The Quantum Mechanical Rotation Operators of Spins 5/2 to 7/2

Abstract

With the creation of logic gates and algorithms for quantum computers and entering our lives, it is predicted that great developments will take place in this area and important efforts are made. Spin rotation processors are quantum mechanical rotation processors and have no classic counterparts. The rotation operators of spin 1/2 are well known and can be found in related textbooks. But rotation operators of other spins greater than 1/2 can be found numerically by evaluating the series expansions of exponential operator obtained from Schr?dinger equation, by evaluation of Wigner-d formula or by recently established expressions in polynomial forms discussed in the text. In a previous paper, rotation operators for spins 1/2 to 2 were published. In this work, explicit symbolic expressions of x, y and z components of rotation operators for spin 5/2, 3 and 7/2 are worked via exponential operator for each element of related spin operators and utilizing simple linear curve fitting process. The procedures gave out exact expressions of each element of the rotation operators.

Share and Cite:

Kocakoç, M. and Tapramaz, R. (2021) The Quantum Mechanical Rotation Operators of Spins 5/2 to 7/2. Journal of Applied Mathematics and Physics, 9, 856-863. doi: 10.4236/jamp.2021.95058.

1. Introduction

Quantum mechanical rotation operators in explicit forms are the essential tools of microscopic systems. Wigner [1] introduced expressions of rotation matrices or Wigner-d matrices for orbital angular momenta with the basis of standard Euler angles. Real rotation operators for total angular momenta of 1/2, 1, 3/2 and 2 were generated from Wigner-d formula in some classical textbooks of quantum mechanics [2] - [8] . A recently published paper on rotation operators by Curtright et al. (2014) and Curtright and Van Kortryk [9] [10] give rotation operator expressions in polynomial form for all angular momenta in Cartesian components. In order to find out the rotation operators in a matrix form, one has to sum up the polynomial terms given, which includes powers of related angular momentum operators.

Pulsed nuclear magnetic resonance, (pulsed-NMRnuclear magnetic resonance), pulsed electron paramagnetic resonance (pulsed-EPRElectron paramagnetic resonance) and pulsed electron nuclear double resonance (pulsed- ENDORElectron nuclear double resonance) spectroscopies, however, utilizes rotation operators in rotating coordinate system or laboratory coordinate system where the spins are polarized along a definite orientation by a known external magnetic field. This direction is defined as z axis and a series of magnetic pulses are applied consequently along laboratory x and/or y axes to rotate the polarized spins around related axes [11] [12] [13] [14] [15] . Besides the spectroscopic applications, spin based quantum-computing systems in which the pulsed magnetic resonance techniques constructed the pioneering pulse arrays, utilizes rotation operators intensively in rotating Cartesian coordinate system [13] [14] .

In this work explicit rotation operator expressions of angular momenta 5/2 to 7/2 are constructed from exponential operators given in equality below for x, y and z components of angular momenta for a series of angles θ between interval θ 1 to θ N , and fitting the obtained values to suitable functions by linear curve fitting procedure.

2. Rotation Operator Expression

We start rewriting the exponential rotation operator expression,

R ^ α = exp ( i θ p J ^ α ) , α = x , y , z (1)

derived from time dependent Schrödinger equation for a rotating magnetic field pulse B1 given as i ψ t = H ^ ψ where H ^ = g β B 1 J ^ α is pulse Hamiltonian applied to a spin system polarized along laboratory z axis by an external magnetic field, where J ^ α ( α = x , y , z ) is Pauli spin matrices. J ^ α can be nuclear spin I or electron spin S, or coupled S and I systems, and g and β have conventional meanings. The rotation angle θ p is defined for rotating pulse duration t p as θ p = g β B 1 t p / = ω p t p .

The exponential function for nuclear spin in equation above can be rewritten by Euler expansion,

R α = exp ( i θ p J ^ α ) = cos ( θ p J ^ α ) + i sin ( θ p J ^ α ) , α = x , y , z (2)

where;

cos ( θ p I ^ α ) = I 1 2 ! θ p 2 J α 2 + 1 4 ! θ p 4 J α 4 1 6 ! θ p 6 J α 6 + = n = 0 1 ( 2 n ) ! θ p 2 n J α 2 n sin ( θ p I ^ α ) = 1 1 ! θ p 1 J α 1 1 3 ! J α 3 θ p 3 + 1 5 ! J α 5 θ p 5 = n = 0 1 ( 2 n + 1 ) ! θ p 2 n + 1 J α 2 n + 1 (3)

where I is the unit matrix with the same dimension as J.

Evaluation of power series given in above requires intensive calculation due to the powers of spin operator matrices. The rotation operators for spin-1/2 systems are rather easy because the elements of the powers of Pauli spin matrices are either zero or unity multiplied by a coefficient. As an example, Pauli matrices for spin-1/2 and corresponding explicit rotation operators are borrowed from textbooks [11] [12] [13] and are given in Equation,

I ^ x = 1 2 [ 0 1 1 0 ] , I ^ y = i 2 [ 0 1 1 0 ] , I ^ z = 1 2 [ 1 0 0 1 ] R ^ x ( θ ) = [ c i s i s c ] , R ^ y ( θ ) = [ c s s c ] , R ^ z ( θ ) = [ z 0 0 z * ] (4)

where c = cos ( θ / 2 ) , s = sin ( θ / 2 ) , z = c + i s .

One of the ways of obtaining explicit expressions of the rotation operators for the spins greater than 1/2 is two-step numerical calculation. In the first step the sine and cosine series given in powers of spin operator matrices are summed up numerically for each element of a spin matrix for angles between e.g. 0˚ and 360˚ with intervals of 5˚, and with the highest possible precision. The results for each angle are found to be real and imaginary elements of rotation matrices.

In the second step, variations of each element of the rotation matrices against rotation angles are fitted to a linear function. The exact fitting functions r i j ( θ ) , found after some trials, are determined and given in equation,

r i j ( θ ) = k = 1 K ξ k cos K k ( θ 2 ) sin k 1 ( θ 2 ) , K = 2 J + 1 ; i , j = 1 , 2 , 3 , , K (5)

where K is the number of terms of the function, θ is rotation angle around x, y or z axis and J is the value of spin (nuclear, electronic or coupled spins). ξ k is the coefficient of k’th term of linear fitting function which is determined by linear curve fitting process. Goodness of all fitting processes was controlled by the value r which is known as regression coefficient or spectral radius, or visually on simultaneous plots of original and fitted curves. All operator matrices obtained were tested by comparing to corresponding operators obtained from Wigner-d formula [1] [2] [3] [4] [5] and operators in polynomial forms [9] .

3. Results and Discussion

The calculations performed for each element of rotation matrices for all components of spins-5/2 to 7/2 produce the operators given below. The symbols c and s are given in Equation 4. It should be noted, just for economy of calculations, that the rotation matrices can be symmetric or anti-symmetric and therefore calculations of only upper triangular matrices will suffice.

Explicit rotation operators for spin-5/2, 3 and 7/2 are given below together with the recommended fit functions. Although all terms of each fitting function are given, most of the coefficients are found to be zero and therefore a few terms will be left as given below. For precise results, only several terms including nonzero coefficients should be taken in fitting processes.

3.1. Rotation Operators of Spin-5/2

Suggested fitting function for rotation operators:

r i j = ξ 1 c 5 + ξ 2 c 4 s + ξ 3 c 3 s 2 + ξ 4 c 2 s 3 + ξ 5 c s 4 + ξ 6 s 5 (6)

Rotation operators were given as,

R ^ x ( θ ) = [ r 11 i r 12 r 13 i r 14 r 15 i r 16 i r 12 r 22 i r 23 r 24 i r 25 r 15 r 13 i r 23 r 33 i r 34 r 24 i r 14 i r 14 r 24 i r 34 r 33 i r 23 r 13 r 15 i r 25 r 24 i r 23 r 22 i r 12 i r 16 r 15 i r 14 r 13 i r 12 r 11 ]

R ^ y ( θ ) = [ r 11 r 12 r 13 r 14 r 15 r 16 r 12 r 22 r 23 r 24 r 25 r 15 r 13 r 23 r 33 r 34 r 24 r 14 r 14 r 24 r 34 r 33 r 23 r 13 r 15 r 25 r 24 r 23 r 22 r 12 r 16 r 15 r 14 r 13 r 12 r 11 ] (7)

R ^ z ( φ ) = [ z 11 0 0 0 0 0 0 z 22 0 0 0 0 0 0 z 33 0 0 0 0 0 0 z 33 * 0 0 0 0 0 0 z 22 * 0 0 0 0 0 0 z 11 * ]

The elements of operators were found as given below;

r 11 = c 5 r 22 = c 5 4 c 3 s 2 r 12 = 5 c 4 s r 23 = 18 c 2 s 3 8 c 4 s r 13 = 10 c 3 s 2 r 24 = 18 c 3 s 2 8 c s 4 r 14 = 10 c 2 s 3 r 25 = 4 c 2 s 3 s 5 r 15 = 10 c s 4 r 33 = c 5 4 c 3 s 2 + 3 c s 4 r 16 = s 5 r 34 = 3 c 4 s + 6 c 2 s 3 s 5

z 11 = ( c 5 10 c 3 s 2 + 5 c s 4 ) + i ( 5 c 4 s 10 c 2 s 3 + s 5 ) z 22 = ( c 5 2 c 3 s 2 3 c s 4 ) + i ( 3 c 4 s + 2 c 2 s 3 s 5 ) z 33 = ( c 5 + 2 c 3 s 2 + c s 4 ) + i ( c 4 s + 2 c 2 s 3 + s 5 ) (8)

3.2. Rotation Operators of Spin-3

Suggested fitting function for rotation operators:

r i j = ξ 1 c 6 + ξ 2 c 5 s + ξ 3 c 4 s 2 + ξ 4 c 3 s 3 + ξ 5 c 2 s 4 + ξ 6 c s 5 + ξ 7 s 6 (9)

Rotation operators were found to be,

R ^ x ( θ ) = [ r 11 i r 12 r 13 i r 14 r 15 i r 16 r 17 i r 12 r 22 i r 23 r 24 i r 25 r 26 i r 16 r 13 i r 23 r 33 i r 34 r 35 i r 25 r 15 i r 14 r 24 i r 34 r 44 i r 34 r 24 i r 14 r 15 i r 25 r 35 i r 34 r 33 i r 23 r 13 i r 16 r 26 i r 25 r 24 i r 23 r 22 i r 12 r 17 i r 16 r 15 i r 14 r 13 i r 12 r 11 ]

R ^ y ( θ ) = [ r 11 r 12 r 13 r 14 r 15 r 16 r 17 r 12 r 22 r 23 r 24 r 25 r 26 r 16 r 13 r 23 r 33 r 34 r 35 r 25 r 15 r 14 r 24 r 34 r 44 r 34 r 24 r 14 r 15 r 25 r 35 r 34 r 33 r 23 r 13 r 16 r 26 r 25 r 24 r 23 r 22 r 12 r 17 r 16 r 15 r 14 r 13 r 12 r 11 ] (10)

R ^ z ( φ ) = [ z 11 0 0 0 0 0 0 0 z 22 0 0 0 0 0 0 0 z 33 0 0 0 0 0 0 0 z 44 0 0 0 0 0 0 0 z 33 * 0 0 0 0 0 0 0 z 22 * 0 0 0 0 0 0 0 z 11 * ]

The elements of operators were found as given below;

r 11 = c 6 r 24 = 30 ( c 4 s 2 + c 2 s 4 ) r 12 = 6 c 5 s r 25 = 40 c 3 s 3 10 c s 5 r 13 = 15 c 4 s 2 r 26 = 5 c 2 s 4 + s 6 r 14 = 20 c 3 s 3 r 33 = c 6 8 c 4 s 2 + 6 c 2 s 4 r 15 = 15 c 2 s 4 r 34 = 108 c 3 s 3 12 ( c 5 s + c s 5 ) r 16 = 6 c s 5 r 35 = 6 c 4 s 2 + 8 c 2 s 4 s 6 r 17 = s 6 r 44 = c 6 9 ( c 4 s 2 + c 2 s 4 ) s 6 r 22 = c 6 5 c 4 s 2 r 23 = 10 c 5 s + 40 c 3 s 3

z 11 = ( c 6 15 c 4 s 2 + 15 c 2 s 4 s 6 ) + i ( 6 c 5 s 20 c 3 s 3 + 6 c s 5 ) z 22 = ( c 6 5 c 4 s 2 5 c 2 s 4 + s 6 ) + 4 i ( c 5 s c s 5 ) z 33 = ( c 6 + c 4 s 2 c 2 s 4 s 6 ) + 2 i ( c 5 s + 2 c 3 s 3 + c s 5 ) z 44 = 1 (11)

3.3. Rotation Operators of Spin-7/2 (Octet)

Suggested fitting function for rotation operators:

r i j = ξ 1 c 7 + ξ 2 c 6 s + ξ 3 c 5 s 2 + ξ 4 c 4 s 3 + ξ 5 c 3 s 4 + ξ 6 c 2 s 5 + ξ 7 c s 6 + ξ 8 s 7 (12)

Rotation operators are given as,

R ^ x ( θ ) = [ r 11 i r 12 r 13 i r 14 r 15 i r 16 r 17 r 18 i r 12 r 22 i r 23 r 24 i r 25 r 26 i r 27 r 17 r 13 i r 23 r 33 i r 34 r 35 i r 36 r 26 i r 16 i r 14 r 24 i r 34 r 44 i r 45 r 35 i r 25 r 15 r 15 i r 25 r 35 i r 45 r 44 i r 34 r 24 i r 14 i r 16 r 26 i r 36 r 35 i r 34 r 33 i r 23 r 13 r 17 i r 27 r 26 i r 25 r 24 i r 23 r 22 i r 12 i r 18 r 17 i r 16 r 15 i r 14 r 13 i r 12 r 11 ]

R ^ y ( θ ) = [ r 11 r 12 r 13 r 14 r 15 r 16 r 17 r 18 r 12 r 22 r 23 r 24 r 25 r 26 r 27 r 17 r 13 r 23 r 33 r 34 r 35 r 36 r 26 r 16 r 14 r 24 r 34 r 44 r 45 r 35 r 25 r 15 r 15 r 25 r 35 r 45 r 44 r 34 r 24 r 14 r 16 r 26 r 36 r 35 r 34 r 33 r 23 r 13 r 17 r 27 r 26 r 25 r 24 r 23 r 22 r 12 r 18 r 17 r 16 r 15 r 14 r 13 r 12 r 11 ] (13)

R ^ z ( φ ) = [ z 11 0 0 0 0 0 0 0 0 z 22 0 0 0 0 0 0 0 0 z 33 0 0 0 0 0 0 0 0 z 44 0 0 0 0 0 0 0 0 z 44 * 0 0 0 0 0 0 0 0 z 33 * 0 0 0 0 0 0 0 0 z 22 * 0 0 0 0 0 0 0 0 z 11 * ]

The elements of operators were found as given below;

r 11 = c 7 r 12 = 7 c 6 s , r 13 = 21 c 5 s 2 r 14 = 35 c 4 s 3 r 15 = 35 c 3 s 4 r 16 = 21 c 2 s 5 r 17 = 7 c s 6 r 18 = s 7 r 22 = c 7 6 c 5 s 2 r 23 = 12 c 6 s + 75 c 4 s 3 r 24 = 45 c 5 s 2 + 80 c 3 s 4 r 25 = 80 c 4 s 3 45 c 2 s 5 r 26 = 75 c 3 s 4 12 c s 6 r 27 = 6 c 2 s 5 + s 7 r 33 = c 7 10 c 5 s 2 c 3 s 4 r 34 = 15 c 6 s 240 c 4 s 3 + 60 c 2 s 5 r 35 = 60 c 5 s 2 240 c 3 s 4 + 15 c s 6 r 36 = 10 c 3 s 4 10 s 5 c 2 + s 7 r 44 = c 7 12 s 2 + 18 c 3 s 4 4 c s 6 r 45 = 4 c 6 s 18 c 4 s 3 + 12 c 2 s 5 s 7

z 11 = ( c 7 21 c 5 s 2 + 35 c 3 s 4 7 c s 6 ) + i ( 7 c 6 s 35 c 4 s 3 21 c 2 s 5 s 7 ) z 22 = ( c 7 9 c 5 s 2 5 c 3 s 4 + 5 c s 6 ) + i ( 5 c 6 s 5 c 4 s 3 9 c 2 s 5 + s 7 ) z 33 = ( c 7 c 5 s 2 5 c 3 s 4 3 c s 6 ) + i ( 3 c 6 s + 5 c 4 s 3 + c 2 s 5 s 7 ) z 44 = ( c 7 + 3 c 5 s 2 + 3 c 3 s 4 + c s 6 ) + i ( c 6 s + 3 c 4 s 3 + 3 c 2 s 5 + s 7 ) (14)

4. Conclusions

The electron spin can form a dual thread, the total spin is 1, and it can form a triplet state or qutrite system. In addition, carbon nanotubes or fullerens will enable three or more paired electron systems. Therefore, it seems necessary to establish the theoretical foundations of large spin systems. Since EPR spectroscopy can work in different spin systems, quantum mechanical rotation for some systems other than spin 1/2, some basic passages and spin processors need to be created. Quantum rotation processors for spin 5/2, 3 and 7/2 were found in this study.

Quantum mechanical rotation operators in explicit forms were constructed by numerical evaluations of exponential rotation operator for spins-5/2 to 7/2. The exponential rotation operator was expanded in power series and numerical values for each operator element at angles consequently between definite intervals, and then the values were fitted to linear trigonometric polynomials of sines and cosines. The rotation operators obtained can be used to determine rotations of corresponding spin around three Cartesian coordinates.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Wigner, E.P. (1959) Group Theory and Its Applications to the Quantum Mechanics of Atomic Spectra. Academic Press, Los Alamos.
[2] Messiah, A. (1967) Quantum Mechanics. Vol. 2, North Holland Publishing Co., Orsay.
[3] Sakurai, J.J. and Napolitano, J.J. (2011) Modern Quantum Mechanics. Cambridge University Press, Cambridge.
[4] Morrison, M.A. and Parker G.A. (1987) A Guide to Rotations in Quantum Mechanics. Australian Journal of Physics, 40, 465-498.
https://doi.org/10.1071/PH870465
[5] Blanca, M.A., Flórez, M. and Bermejo, M. (1997) Evaluation of the Rotation Matrices in the Basis of Real Spherical Harmonics. Journal of Molecular Structure: THEOCHEM, 419, 19-27.
https://doi.org/10.1016/S0166-1280(97)00185-1
[6] Dachsel, H. (2006) Fast and Accurate Determination of the Wigner Rotation Matrices in the Fast Multipole Method. Journal of Chemical Physics, 124, Article ID: 144115.
https://doi.org/10.1063/1.2194548
[7] Gimbutas, Z. and Greengard, L. (2009) A Fast and Stable Method for Rotating Spherical Harmonic Expansions. Journal of Computational Physics, 228, 5621-562.
https://doi.org/10.1016/j.jcp.2009.05.014
[8] Aubert, G. (2013) An Alternative to Wigner D-Matrices for Rotating Real Spherical Harmonics. AIP Advances, 3, Article ID: 062121.
https://doi.org/10.1063/1.4811853
[9] Curtright, T.L. Fairlie, D.B. and Zachos, C.K. (2014) mCompact Formula for Rotations as Spin Matrix Polynomials. Symmetry, Integrability and Geometry: Methods and Applications, 10, 1-15.
https://doi.org/10.3842/SIGMA.2014.084
[10] Curtright, T.L. and Van Kortryk T.S. (2014) On Rotations as Spin Matrix Polynomials. Journal of Physics A: Mathematical and Theoretical, 48, Article No. 025202.
https://doi.org/10.1088/1751-8113/48/2/025202
[11] Fukushima, E. and Roeder, S.B.W. (1993) Experimental Pulse NMR: A Nuts and Bolts Approach. CRC Press, Boca Raton.
[12] Rule, G.S. and Kelvin Hitchens, T. (2006) Fundamentals of Protein NMR Spectroscopy. Springer, Netherlands.
[13] Oliviera, I., Sarthour Jr., R., Bonagamba, T., Azevedo, E. and Freitas, J.C.C. (2007) NMR Quantum Information Processing. Elsevier, Amsterdam.
[14] Jones, J.A. (2001) NMR Quantum Computation. Progress in Nuclear Magnetic Resonance Spectroscopy, 38, 325-360.
https://doi.org/10.1016/S0079-6565(00)00033-9
[15] Schweiger, A. and Jeschke, G. (2001) Principles of Pulse Electron Paramagnetic Resonance. Oxford University Press, Oxford.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.