Foldy–Wouthuysen transformation

"FW transformation" redirects to here.

The Foldy–Wouthuysen transform is widely used in high energy physics. It was historically formulated by Leslie Lawrance Foldy and Siegfried Adolf Wouthuysen in 1949 to understand the nonrelativistic limit of the Dirac equation, the equation for the spin-1/2 particles.^[1][2][3][4] A detailed general discussion of the Foldy–Wouthuysen-type transformations in particle interpretation of relativistic wave equations is in Acharya and Sudarshan (1960).^[5]

A canonical transform

Foldy and Wouthuysen made use of a canonical transform that has now come to be known as the Foldy–Wouthuysen transformation. A brief account of the history of the transformation is to be found in the obituaries of Foldy and Wouthuysen^{[6] [7]} and the biographical memoir of Foldy. ^[8] Before their work, there was some difficulty in understanding and gathering all the interaction terms of a given order, such as those for a Dirac particle immersed in an external field. With their procedure the physical interpretation of the terms was clear, and it became possible to apply their work in a systematic way to a number of problems that had previously defied solution.^[9][10] The Foldy–Wouthuysen transform was extended to the physically important cases of the spin-0 and the spin-1 particles,^[11] and even generalized to the case of arbitrary spins.^[12]

The Foldy–Wouthuysen (FW) transformation (after Lesley L. Foldy and Siegfried A. Wouthuysen) is a unitary transformation on a fermion wave function of the form:

$\psi \to \psi '=U\psi$

(1)

where the unitary operator is the 4 × 4 matrix:

$U=e^{{\beta {\mathbf {\alpha }}\cdot {\hat {p}}\theta }}=\cos \theta +\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin \theta =e^{{{\mathbf {\gamma }}\cdot {\hat {p}}\theta }}=\cos \theta +{\mathbf {\gamma }}\cdot {\hat {p}}\sin \theta$

(2)

Above, ${\hat {p^{i}}}\equiv p^{i}/|p|$ is the unit vector oriented in the direction of the fermion momentum. The above are related to the Dirac matrices by β = γ⁰ and αⁱ = γ⁰γⁱ, with i = 1, 2, 3. A straightforward series expansion applying the commutativity properties of the Dirac matrices demonstrates that (2) above is true. The inverse

$U^{{-1}}=e^{{-\beta {\mathbf {\alpha }}\cdot {\hat {p}}\theta }}=\cos \theta -\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin \theta \,,$

so it is clear that U⁻¹U = I, where I is a 4 × 4 identity matrix.

Foldy–Wouthuysen transformation of the Dirac Hamiltonian for a free fermion

This transformation is of particular interest when applied to the free-fermion Dirac Hamiltonian operator

${\hat {H}}_{0}\equiv \alpha \cdot p+\beta m$

in bi-unitary fashion, in the form:

${\hat {H}}_{0}\to {\hat {H}}'_{0}\equiv U{\hat {H}}_{0}U^{{-1}}=U(\alpha \cdot p+\beta m)U^{{-1}}=(\cos \theta +\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin \theta )(\alpha \cdot p+\beta m)(\cos \theta -\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin \theta )$

(3)

Using the commutativity properties of the Dirac matrices, this can be massaged over into the double-angle expression:

${\hat {H}}'_{0}=(\alpha \cdot p+\beta m)(\cos \theta -\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin \theta )^{{2}}=(\alpha \cdot p+\beta m)e^{{-2\beta {\mathbf {\alpha }}\cdot {\hat {p}}\theta }}=(\alpha \cdot p+\beta m)(\cos 2\theta -\beta {\mathbf {\alpha }}\cdot {\hat {p}}\sin 2\theta )$

(4)

This factors out into:

${\hat {H}}'_{0}=\alpha \cdot p\left(\cos 2\theta -{\frac {m}{|p|}}\sin 2\theta \right)+\beta (m\cos 2\theta +|p|\sin 2\theta )$

(5)

Choosing a particular representation: Newton–Wigner

Clearly, the FW transformation is a continuous transformation, that is, one may employ any value for θ which one chooses. Now comes the distinct question of choosing a particular value for θ, which amounts to choosing a particular transformed representation.

One particularly important representation, is that in which the transformed Hamiltonian operator ${\hat {H}}'_{0}$ is diagonalized. Clearly, a completely diagonalized representation can be obtained by choosing θ such that the α · p term in (5) is made to vanish. Such a representation is specified by defining:

$\tan 2\theta \equiv |p|/m$

(6)

so that (5) is reduced to the diagonalized (this presupposes that β is taken in the Dirac-Pauli representation (after Paul Dirac and Wolfgang Pauli) in which it is a diagonal matrix):

${\hat {H}}'_{0}=\beta (m\cos 2\theta +|p|\sin 2\theta )$

(7)

By elementary trigonometry, (6) also implies that:

$\sin 2\theta =|p|/{\sqrt {m^{2}+|p|^{2}}}$ and $\cos 2\theta =m/{\sqrt {m^{2}+|p|^{2}}}$

(8)

so that using (8) in (7) now leads following reduction to:

${\hat {H}}'_{0}=\beta {\sqrt {m^{2}+|p|^{2}}}$

(9)

This calculation can be examined in further detail in the following link.

Prior to Foldy and Wouthuysen publishing their transformation, it was already known that (9) is the Hamiltonian in the Newton–Wigner (NW) representation (named after Theodore Duddell Newton and Eugene Wigner) of the Dirac equation. What (9) therefore tells us, is that by applying a FW transformation to the Dirac-Pauli representation of Dirac's equation, and then selecting the continuous transformation parameter θ so as to diagonalize the Hamiltonian, one arrives at the NW representation of Dirac's equation, because NW itself already contains the Hamiltonian specified in (9). See this link

If one considers an "on shell" mass—fermion or otherwise—given by m² = p^σp_σ, and employs a Minkowski metric tensor for which diag(η) = (+1, −1, −1, −1), it should be apparent that the expression

$p^{0}={\sqrt {m^{2}+{|p|}^{2}}}$

is equivalent to the E ≡ p⁰ component of the energy-momentum vector p^μ, so that (9) is alternatively specified rather simply by ${\hat {H}}'_{0}=\beta E$.

Correspondence between the Dirac-Pauli and Newton–Wigner representations, for a fermion at rest

Now let us consider a fermion "at rest", which we may define in this context as a fermion for which |p| = 0. From (6) or (8), this means that cos(2θ) = 1, so that θ = 0, ±π, ±2π and, from (2), that the unitary operator U = ±I. Therefore, any operator O in the Dirac-Pauli representation upon which we perform a bi-unitary transformation, will be given, for an "at rest" fermion, by:

$O\to O'\equiv UOU^{{-1}}=(\pm I)(O)(\pm I)=O$

(10)

Contrasting the original Dirac-Pauli Hamiltonian Operator

${\hat {H}}_{0}\equiv \alpha \cdot p+\beta m$

with the NW Hamiltonian (9), we do indeed find the |p| = 0 "at rest" correspondence:

${\hat {H}}_{0}={\hat {H}}'_{0}=\beta m$

(11)

The velocity operator in the Dirac-Pauli representation

Now, let us consider the velocity operator. To obtain this operator, we must commute the Hamiltonian operator ${\hat {H}}_{0}$ with the canonical position operators $x_{i}$, i.e., we must calculate ${\hat {v_{i}}}\equiv i[{\hat {H}}_{0},x_{i}]$. One good way to approach this calculation, is to start by writing the scalar rest mass m as $m=\gamma ^{0}{\hat {H}}_{0}+\gamma ^{j}p_{j}$, and then to mandate that the scalar rest mass commute with the $x_{i}$. Thus, we may write:

$0=[m,x_{i}]=[(\gamma ^{0}{\hat {H}}_{0}+\gamma ^{j}p_{j}),x_{i}]=[\gamma ^{0}{\hat {H}}_{0},x_{i}]+i\gamma _{i}$

(12)

where we have made use of the Heisenberg canonical commutation relationship $[x_{i},p_{j}]=-i\eta _{{ij}}$ to reduce terms. Then, multiplying from the left by $\gamma ^{0}$ and rearranging terms, we arrive at:

${\frac {{\hat {dx_{i}}}}{dt}}={\hat {v_{i}}}\equiv i[{\hat {H}}_{0},x_{i}]=\alpha _{i}$

(13)

Because the canonical relationship $i[{\hat {H}}_{0},{\hat {v}}_{i}]\neq 0$, the above provides the basis for computing an inherent, non-zero acceleration operator, which specifies the oscillatory motion known as Zitterbewegung.

The velocity operator in the Newton–Wigner representation

In the Newton–Wigner representation, we now wish to calculate ${\hat {v_{i}}}'\equiv i[{\hat {H}}'_{0},x_{i}]$. If we use the result at the very end of section 2 above, ${\hat {H}}'_{0}=\beta p_{0}$, then this can be written instead as:

${\hat {v_{i}}}'\equiv i[{\hat {H}}'_{0},x_{i}]=i\beta [p_{0},x_{i}]$

(14)

Using the above, we need simply to calculate $[p_{0},x_{i}]$, then multiply by iβ.

The canonical calculation proceeds similarly to the calculation in section 4 above, but because of the square root expression in $p^{0}={\sqrt {m^{2}+|p|^{2}}}$, one additional step is required.

First, to accommodate the square root, we will wish to require that the scalar square mass m² commute with the canonical coordinates x_i, which we write as:

$0\equiv [m^{2},x_{i}]=[(p^{0}p_{0}+p^{j}p_{j}),x_{i}]=[p^{0}p_{0},x_{i}]+2ip_{i}$

(15)

where we again use the Heisenberg canonical relationship $[x_{i},p_{j}]=-i\eta _{{ij}}$. Then, we need an expression for $[p_{0},x_{i}]$ which will satisfy (15). It is straightforward to verify that:

$i[p_{0},x_{i}]={\frac {p_{i}}{p^{0}}}=v_{i}$

(16)

will satisfy (15) when again employing $[x_{i},p_{j}]=-i\eta _{{ij}}$. Now, we simply return the iβ factor via (14), to arrive at:

${\frac {{\hat {dx_{i}}}'}{dt}}={\hat {v_{i}}}'\equiv i[{\hat {H}}'_{0},x_{i}]=\beta {\frac {p_{i}}{p^{0}}}=\beta v_{i}$

(17)

This is understood to be the velocity operator in the Newton–Wigner representation. Because:

$i[{\hat {H}}'_{0},{\hat {v}}_{i}']=i[\beta p_{0},\beta v_{i}]=0$

(18)

it is commonly thought that the Zitterbewegung motion arising out of (12), vanishes when a fermion is transformed into the Newton–Wigner representation.

The velocity operators for a fermion at rest

Now, let us compare equations (13) and (17) for a fermion "at rest", defined earlier in section 3 as a fermion for which |p| = 0. Here, (13) remains:

${\hat {v_{i}}}\equiv i[{\hat {H}}_{0},x_{i}]=\alpha _{i}$

(19)

while (17) becomes:

${\hat {v_{i}}}'\equiv i[{\hat {H}}'_{0},x_{i}]=\beta {\frac {p_{i}}{p^{0}}}=0$

(20)

In equation (10) we found that for a fermion at rest, O′ = O for any operator. One would expect this to include:

${\hat {v_{i}}}'={\hat {v_{i}}}$

(21)

however, equations (19) and (20) for a |p| = 0 fermion appear to contradict (21).

Similar alternatives - perturbative schemes

Starting with the one-particle Dirac equation written earlier with ${\hat {H}}_{0}$ and rewritten here as:

$(c{\vec {\alpha }}\cdot {\vec {p}}+\beta m_{o}c^{2}+IV)\Phi =E\Phi$

where I = I₄ is the 4 × 4 unit matrix. This Hamiltonian is rewritten, namely divided into two parts:

${\hat {H}}\equiv {\hat {H}}_{0}+\lambda {\hat {H}}_{1}$

where

${\hat {H}}_{0}=c{\vec {\alpha }}\cdot {\vec {p}}+\beta m_{o}c^{2}+{\frac {1}{2}}(I+\beta )V$

and

${\hat {H}}_{1}={\frac {1}{\alpha ^{2}}}{\hat {H'}}_{1}={\frac {1}{2\alpha ^{2}}}(I-\beta )V$

where α ≈ 1/137 is the Fine-structure constant (not to be confused with the Dirac alpha matrices). Letting

$\phi ^{0}=A\left({\begin{array}{c}\Psi _{1}^{0}\\\Psi _{2}^{0}\end{array}}\right)$

into the zero order equation for ${\hat {H}}_{0}$ and using a particular but known representation of the Dirac operators, yields:

$\Psi _{2}^{0}=c~{\frac {{\vec {\sigma }}\cdot {\vec {p}}}{(E_{0}+m_{0}c^{2})}}\Psi _{1}^{0}$

where σ_i are the 2 × 2 Pauli matrices. Note that the potential V does not appear in the equation above. The equation for the other spinor is:

$c\sigma \cdot {\vec {p}}~\Psi _{2}^{0}+V\Psi _{1}^{0}=\epsilon _{0}\Psi _{1}^{0}$

where $\epsilon _{0}=E_{0}-m_{0}c^{2}$. Eliminating $\Psi _{2}^{0}$ gives:

$\left\{{\frac {p^{2}}{2m_{0}}}+\left(1+{\frac {\epsilon _{0}}{2m_{0}c^{2}}}\right)(V-\epsilon _{0})\right\}\Psi _{1}^{0}\equiv 0.$

This is simply the non-relativistic equation for a system with a re-normalized potential and energy eigenvalue:

$E'=\left(1+{\frac {\epsilon _{0}}{2m_{0}c^{2}}}\right)\epsilon _{0}\quad {\mbox{and}}\quad V'=\left(1+{\frac {\epsilon _{0}}{2m_{0}c^{2}}}\right)V$

The higher-order corrections can be obtained by conventional perturbation theory. This is known as Moore's decoupling technique. Though it resembles the FW transformation, it is computationally and conceptually much simpler. Though misunderstood at first, in part because the fine structure constant appears in both the equations and the order parameter λ requiring care in the "bookkeeping" of the perturbative scheme, Moore's decoupling technique was vindicated for the (relativistic) hydrogen atom using conventional Rayleigh Schrödinger perturbation theory and computer algebra and proven to converge to the correct solution.^[13]

It has been applied successfully to relativistic calculations on Alkali metals and represents one of many relativistic perturbative schemes investigated by Werner Kutzelnigg.^[14][15]

Other applications

The powerful machinery of the Foldy–Wouthuysen transform originally developed for the Dirac equation has found applications in many situations such as acoustics, and optics.

It has found applications in very diverse areas such as atomic systems^[16][17] synchrotron radiation^[18] and derivation of the Bloch equation for polarized beams.^[19]

The applications of the Foldy–Wouthuysen transformation in acoustics is very natural; comprehensive and mathematically rigorous accounts.^[20][21][22]

In the traditional scheme the purpose of expanding the light optics Hamiltonian

$\widehat {H}=-\left(n^{2}(r)-\widehat {p}_{\perp }^{2}\right)^{{1/2}}$

in a series using

$\left({\widehat {p}_{\perp }^{2}}/{n_{0}^{2}}\right)$

as the expansion parameter is to understand the propagation of the quasiparaxial beam in terms of a series of approximations (paraxial + nonparaxial). Similar is the situation in the case of the charged-particle optics. Let us recall that in relativistic quantum mechanics too one has a similar problem of understanding the relativistic wave equations as the nonrelativistic approximation plus the relativistic correction terms in the quasirelativistic regime. For the Dirac equation (which is first order in time) this is done most conveniently using the Foldy–Wouthuysen transformation leading to an iterative diagonalization technique. The main framework of the newly developed formalisms of optics (both light optics and the charged-particle optics) is based on the transformation technique of the Foldy–Wouthuysen theory which casts the Dirac equation in a form displaying the different interaction terms between the Dirac particle and an applied electromagnetic field in a nonrelativistic and easily interpretable form.

In the Foldy–Wouthuysen theory the Dirac equation is decoupled through a canonical transformation into two two-component equations: one reduces to the Pauli equation^[23] in the nonrelativistic limit and the other describes the negative-energy states. It is possible to write a Dirac-like matrix representations of Maxwell's equations. In such a matrix form the Foldy–Wouthuysen can be applied.^{[24] [25] [26] [27] [28] [29]}

There is a close algebraic analogy between the Helmholtz equation (governing scalar optics) and the Klein–Gordon equation; and the matrix-form of the Maxwell's equations (governing vector optics) and the Dirac equation. So, it is natural to use the powerful machinery of standard quantum mechanics (particularly, the Foldy–Wouthuysen transform) in analyzing these systems.

The suggestion to employ the Foldy–Wouthuysen Transformation technique in the case of the Helmholtz equation was mentioned in the literature as a remark.^[30]

It was only in the recent works, that this idea was exploited to analyze the quasiparaxial approximations for specific beam optical system.^[31] The Foldy–Wouthuysen technique is ideally suited for the Lie algebraic approach to optics. With all these plus points, the powerful and ambiguity-free expansion, the Foldy–Wouthuysen Transformation is still little used in optics. The technique of the Foldy–Wouthuysen Transformation results in what we call as the non-traditional prescriptions of Helmholtz optics^[32] and Maxwell optics^[33] respectively. The non-traditional approaches give rise to very interesting wavelength-dependent modifications of the paraxial and aberrating behaviour. The non-traditional formalism of Maxwell optics provides a unified framework of light beam optics and polarization. The non-traditional prescriptions of light optics are in close analogy with the quantum theory of charged-particle beam optics.^{[34][35][36][37]} In optics, it has enabled to see the deeper connections in the wavelength-dependent regime between light optics and charged-particle optics (see Electron optics).^{[38] [39]}

See also

• Relativistic quantum mechanics

Notes

1. Foldy, L. L. and Wouthuysen, S. A. (1950). On the Dirac Theory of Spin 1/2 Particles and its Non-Relativistic Limit. Physical Review, 78, 29-36.
2. Foldy, L. L. (1952). The Electromagnetic Properties of the Dirac Particles. Physical Review, 87, (5), 682-693.
3. Pryce, M. H. L. (1948). The mass-centre in the restricted theory of relativity and its connexion with the quantum theory of elementary particles. Proceedings of the Royal Society of London, Series A: Mathematical and Physical Sciences, A195, 62-81.
4. Tani, S. (1951). Connection between particle models and field theories. I. The case spin 1/2. Progress of Theoretical Physics, 6, 267-285.
5. Acharya, R., and Sudarshan, E. C. G. (1960). Front Description in Relativistic Quantum Mechanics. Journal of Mathematical Physics, 1, 532-536.
6. Brown, R. W., Krauss, L. M., and Taylor, P. L. (2001). Obituary of Leslie Lawrence Foldy. Physics Today, 54 (12), 75.
7. Leopold, H. (1997). Obituary of Siegfried A Wouthuysen. Physics Today, 50, (11), 89.
8. Foldy, L. L. (2006). Origins of the FW Transformation: A Memoir, Appendix G, pp. 347-351 in Physics at a Research University, Case Western Reserve University 1830-1990, Ed. William Fickinger.
9. Bjorken, J. D., and Drell, S. D. (1964). Relativistic Quantum Mechanics. McGraw-Hill, New York, San Francisco.
10. Costella, J. P., and McKellar, B. H. J. (1995). The Foldy-Woutuysen transformation. .
11. Case, K. M. (1954). Some generalizations of the Foldy–Wouthuysen transformation. Physical Review, 95, 1323-1328.
12. Jayaraman, J. (1975). A note on the recent Foldy–Wouthuysen transformations for particles of arbitrary spin. J. Phys. A: Math. Gen., 8, L1-L4.
13. T.C. Scott, R.A. Moore, G.J. Fee, M.B. Monagan and E.R. Vrscay, "Perturbative Solutions of Quantum Mechanical Problems by Symbolic Computation", J. Comp. Phys., 87, 366-395 (1990).
14. W. Kutzelnigg,"Perturbation theory of relativistic corrections. II. Analysis and classification of known and other possible methods.", Z. Phys. D 15, 27 (1990).
15. W. Kutzelnigg,"Perturbation theory of relativistic effects", in Relativistic Electronic Structure Theory, Part I, P.Schwerdtfeger ed. Elsevier (2002).
16. Asaga, T., Fujita, T., and Hiramoto M. (2000). Foldy–Wouthuysen transformed EDM operator in atomic system. .
17. Pachucki, K. (2004). Higher-order effective Hamiltonian for light atomic systems, .
18. Lippert, M., Bruckel, Th., Kohler, Th., and Schneider, J. R. (1994). High-Resolution Bulk Magnetic Scattering of High-Energy Synchrotron Radiation. Europhys. Lett. 27(7), 537-541.
19. Heinemann, K., and Barber, D. P. (1999). The semiclassical Foldy–Wouthuysen transformation and the derivation of the Bloch equation for spin-1/2 polarized beams using Wigner functions, , Chen P (ed.): Proceedings of the 15th Advanced ICFA Beam Dynamics Workshop on Quantum Aspects of Beam Physics, 4–9 January 1998, Monterey, California, USA, World Scientific, Singapore.
20. Fishman, L. (1992). Exact and operator rational approximate solutions of the Helmholtz, Weyl composition equation in underwater acoustics-the quadratic profile, J. Math. Phys., 33, (5), 1887-1914.
21. Fishman, L. (2004). One-way wave equation modeling in two-way wave propagation problems, in: B. Nilsson, L. Fishman (Eds.), Mathematical Modelling of Wave Phenomena 2002, Mathematical Modelling in Physics, Engineering and Cognitive Sciences, vol. 7, Vaxjo University Press, Vaxjo, Sweden, pp. 91-111.
22. Wurmser, D. (2004). A parabolic equation for penetrable rough surfaces: using the Foldy–Wouthuysen transformation to buffer density jumps. Annals of Physics, 311, 53-80.
23. Osche, G. R. (1977). Dirac and Dirac-Pauli equation in the Foldy–Wouthuysen representation, Physical Review D, 15 (8) 2181-2185.
24. Bialynicki-Birula, I. (1996b). Photon wave function. in Progress in Optics, Vol. XXXVI, Emil Wolf. (ed.), Elsevier, Amsterdam, 245-294.
25. Sameen Ahmed Khan. (2002). Maxwell Optics: I. An exact matrix representation of the Maxwell equations in a medium. E-Print: http://arXiv.org/abs/physics/0205083/.
26. Sameen Ahmed Khan. (2005). An Exact Matrix Representation of Maxwell's Equations. Physica Scripta, 71(5), 440-442.
27. Laporte, O., and Uhlenbeck, G. E. (1931). Applications of spinor analysis to the Maxwell and Dirac Equations. Physical Review, 37, 1380-1397.
28. Majorana, E. (1974). (unpublished notes), quoted after Mignani, R., Recami, E., and Baldo, M. About a Diraclike Equation for the Photon, According to Ettore Majorana. Lett. Nuovo Cimento, 11, 568-572.
29. Moses, E. (1959).Solutions of Maxwell's equations in terms of a spinor notation: the direct and inverse problems. Physical Review, 113(6), 1670-1679.
30. Fishman, L., and McCoy, J. J. (1984).Derivation and Application of Extended Parabolic Wave Theories. Part I. The Factored Helmholtz Equation. Journal of Mathematical Physics, 25, 285-296.
31. Sameen Ahmed Khan, Ramaswamy Jagannathan and Rajiah Simon, (2002). Foldy–Wouthuysen transformation and a quasiparaxial approximation scheme for the scalar wave theory of light beams. .
32. Sameen Ahmed Khan, (2005). Wavelength-dependent modifications in Helmholtz Optics. International Journal of Theoretical Physics, 44 (1), 95-125. (Kluwer Academic Publishers, ).
33. Sameen Ahmed Khan, (2006). Wavelength-Dependent Effects in Light Optics. in New Topics in Quantum Physics Research, Editors: Volodymyr Krasnoholovets and Frank Columbus, Nova Science Publishers, New York, . pp. 163-204.
34. Jagannathan, R. Simon, E. C. G. Sudarshan and N. Mukunda, Quantum theory of magnetic electron lenses based on the Dirac equation, Physics Letters A, 134, 457-464 (1989).
35. Jagannathan, Quantum theory of electron lenses based on the Dirac equation, Physical Review A, 42, 6674-6689 (1990).
36. R. Jagannathan and S. A. Khan, Quantum theory of the optics of charged particles, Advances in Imaging and Electron Physics, Editors: Peter W. Hawkes, B. Kazan and T. Mulvey, (Academic Press, Logo, San Diego, 1996), Vol. 97, 257-358 (1996).
37. Conte, M., Jagannathan, R., Khan, S. A., and Pusterla, M. (1996). Beam optics of the Dirac particle with anomalous magnetic moment. Particle Accelerators, 56, 99-126.
38. Sameen Ahmed Khan, (2006). The Foldy–Wouthuysen Transformation Technique in Optics. Optik-International Journal for Light and Electron Optics.117, (10), pp. 481-488.
39. Sameen Ahmed Khan, The Foldy–Wouthuysen Transformation Technique in Optics, in Advances in Imaging and Electron Physics, Volume 152, pp. 49–78 (Elsevier, 2008).