Gordon decomposition

Last updated May 11, 2022

In mathematical physics, the Gordon decomposition^[1] (named after Walter Gordon) of the Dirac current is a splitting of the charge or particle-number current into a part that arises from the motion of the center of mass of the particles and a part that arises from gradients of the spin density. It makes explicit use of the Dirac equation and so it applies only to "on-shell" solutions of the Dirac equation.

Original statement

For any solution $\psi$ of the massive Dirac equation,

(i\gamma ^{\mu }\nabla _{\mu }-m)\psi =0,

the Lorentz covariant number-current $j^{\mu }={\bar {\psi }}\gamma ^{\mu }\psi$ may be expressed as

{\bar {\psi }}\gamma ^{\mu }\psi ={\frac {i}{2m}}({\bar {\psi }}\nabla ^{\mu }\psi -(\nabla ^{\mu }{\bar {\psi }})\psi )+{\frac {1}{m}}\partial _{\nu }({\bar {\psi }}\Sigma ^{\mu \nu }\psi ),

where

\Sigma ^{\mu \nu }={\frac {i}{4}}[\gamma ^{\mu },\gamma ^{\nu }]

is the spinor generator of Lorentz transformations, and

{\bar {\psi }}=\psi ^{\dagger }\gamma ^{0}

is the Dirac adjoint.

The corresponding momentum-space version for plane wave solutions $u(p)$ and ${\bar {u}}(p')$ obeying

(\gamma ^{\mu }p_{\mu }-m)u(p)=0

{\bar {u}}(p')(\gamma ^{\mu }p'_{\mu }-m)=0,

is

{\bar {u}}(p')\gamma ^{\mu }u(p)={\bar {u}}(p')\left[{\frac {(p+p')^{\mu }}{2m}}+i\sigma ^{\mu \nu }{\frac {(p'-p)_{\nu }}{2m}}\right]u(p)~,

where

\sigma ^{\mu \nu }=2\Sigma ^{\mu \nu }.

Proof

One sees that from Dirac's equation that

{\bar {\psi }}\gamma ^{\mu }(m\psi )={\bar {\psi }}\gamma ^{\mu }(i\gamma ^{\nu }\nabla _{\nu }\psi )

and, from the adjoint of Dirac's equation,

({\bar {\psi }}m)\gamma ^{\mu }\psi =((\nabla _{\nu }{\bar {\psi }})(-i\gamma ^{\nu }))\gamma ^{\mu }\psi .

Adding these two equations yields

{\bar {\psi }}\gamma ^{\mu }\psi ={\frac {i}{2m}}({\bar {\psi }}\gamma ^{\mu }\gamma ^{\nu }\nabla _{\nu }\psi -(\nabla _{\nu }{\bar {\psi }})\gamma ^{\nu }\gamma ^{\mu }\psi ).

From Dirac algebra, one may show that Dirac matrices satisfy

\gamma ^{\mu }\gamma ^{\nu }=\eta ^{\mu \nu }-i\sigma ^{\mu \nu }=\eta ^{\nu \mu }+i\sigma ^{\nu \mu }.

Using this relation,

{\bar {\psi }}\gamma ^{\mu }\psi ={\frac {i}{2m}}({\bar {\psi }}(\eta ^{\mu \nu }-i\sigma ^{\mu \nu })\nabla _{\nu }\psi -(\nabla _{\nu }{\bar {\psi }})(\eta ^{\mu \nu }+i\sigma ^{\mu \nu })\psi ),

which amounts to just the Gordon decomposition, after some algebra.

Utility

The second, spin-dependent, part of the current coupled to the photon field, $-A_{\mu }j^{\mu }$ yields, up to an ignorable total divergence,

-{\frac {e\hbar }{2mc}}\partial _{\nu }A_{\mu }{\bar {\psi }}\sigma ^{\nu \mu }\psi =-{\frac {e\hbar }{2mc}}{\tfrac {1}{2}}F_{\mu \nu }{\bar {\psi }}\sigma ^{\mu \nu }\psi ,

that is, an effective Pauli moment term, $-(e\hbar /2mc){\vec {B}}\cdot \psi ^{\dagger }{\vec {\sigma }}\psi$ .

Massless generalization

This decomposition of the current into a particle number-flux (first term) and bound spin contribution (second term) requires $m\neq 0$ .

If one assumed that the given solution has energy ${\textstyle E={\sqrt {|{\mathbf {k} }|^{2}+m^{2}}}}$ so that ${\textstyle \psi (\mathbf {r} ,t)=\psi ({\mathbf {r} })e^{-iEt}}$ , one might obtain a decomposition that is valid for both massive and massless cases.^[2]

Using the Dirac equation again, one finds that

{\mathbf {j} }\equiv e{\bar {\psi }}{\boldsymbol {\gamma }}\psi ={\frac {e}{2iE}}\left(\psi ^{\dagger }\nabla \psi -(\nabla \psi ^{\dagger })\psi \right)+{\frac {e}{E}}(\nabla \times {\mathbf {S} }).

Here ${\boldsymbol {\gamma }}=(\gamma ^{1},\gamma ^{2},\gamma ^{3})$ , and ${\mathbf {S} }=\psi ^{\dagger }{\hat {\mathbf {S} }}\psi$ with $({\hat {S}}_{x},{\hat {S}}_{y},{\hat {S}}_{z})=(\Sigma ^{23},\Sigma ^{31},\Sigma ^{12}),$ so that

{\hat {\mathbf {S} }}={\frac {1}{2}}\left[{\begin{matrix}{\boldsymbol {\sigma }}&0\\0&{\boldsymbol {\sigma }}\end{matrix}}\right],

where ${\boldsymbol {\sigma }}=(\sigma _{x},\sigma _{y},\sigma _{z})$ is the vector of Pauli matrices.

With the particle-number density identified with $\rho =\psi ^{\dagger }\psi$ , and for a near plane-wave solution of finite extent, one may interpret the first term in the decomposition as the current ${\mathbf {j} }_{\rm {free}}=e\rho {\mathbf {k} }/E=e\rho {\mathbf {v} }$ , due to particles moving at speed ${\mathbf {v} }={\mathbf {k} }/E$ .

The second term, ${\mathbf {j} }_{\rm {bound}}=(e/E)\nabla \times {\mathbf {S} }$ is the current due to the gradients in the intrinsic magnetic moment density. The magnetic moment itself is found by integrating by parts to show that

{\boldsymbol {\mu }}{\stackrel {\rm {}}{=}}{\frac {1}{2}}\int {\mathbf {r} }\times {\mathbf {j} }_{\rm {bound}}\,d^{3}x={\frac {1}{2}}\int {\mathbf {r} }\times \left({\frac {e}{E}}\nabla \times {\mathbf {S} }\right)\,d^{3}x={\frac {e}{E}}\int {\mathbf {S} }\,d^{3}x~.

For a single massive particle in its rest frame, where $E=m$ , the magnetic moment reduces to

{\boldsymbol {\mu }}_{\rm {Dirac}}=\left({\frac {e}{m}}\right){\mathbf {S} }=\left({\frac {eg}{2m}}\right){\mathbf {S} }.

where $|{\mathbf {S} }|=\hbar /2$ and $g=2$ is the Dirac value of the gyromagnetic ratio.

For a single massless particle obeying the right-handed Weyl equation, the spin-1/2 is locked to the direction ${\hat {\mathbf {k} }}$ of its kinetic momentum and the magnetic moment becomes^[3]

{\boldsymbol {\mu }}_{\rm {Weyl}}=\left({\frac {e}{E}}\right){\frac {\hbar {\hat {\mathbf {k} }}}{2}}.

Angular momentum density

For both the massive and massless cases, one also has an expression for the momentum density as part of the symmetric Belinfante–Rosenfeld stress–energy tensor

T_{\rm {BR}}^{\mu \nu }={\frac {i}{4}}({\bar {\psi }}\gamma ^{\mu }\nabla ^{\nu }\psi -(\nabla ^{\nu }{\bar {\psi }})\gamma ^{\mu }\psi +{\bar {\psi }}\gamma ^{\nu }\nabla ^{\mu }\psi -(\nabla ^{\mu }{\bar {\psi }})\gamma ^{\nu }\psi ).

Using the Dirac equation one may evaluate $T_{\rm {BR}}^{0\mu }=({\mathcal {E}},{\mathbf {P} })$ to find the energy density to be ${\mathcal {E}}=E\psi ^{\dagger }\psi$ , and the momentum density,

{\mathbf {P} }={\frac {1}{2i}}\left(\psi ^{\dagger }(\nabla \psi )-(\nabla \psi ^{\dagger })\psi \right)+{\frac {1}{2}}\nabla \times {\mathbf {S} }.

If one used the non-symmetric canonical energy-momentum tensor

T_{\rm {canonical}}^{\mu \nu }={\frac {i}{2}}({\bar {\psi }}\gamma ^{\mu }\nabla ^{\nu }\psi -(\nabla ^{\nu }{\bar {\psi }})\gamma ^{\mu }\psi ),

one would not find the bound spin-momentum contribution.

By an integration by parts one finds that the spin contribution to the total angular momentum is

\int {\mathbf {r} }\times \left({\frac {1}{2}}\nabla \times {\mathbf {S} }\right)\,d^{3}x=\int {\mathbf {S} }\,d^{3}x.

This is what is expected, so the division by 2 in the spin contribution to the momentum density is necessary. The absence of a division by 2 in the formula for the current reflects the $g=2$ gyromagnetic ratio of the electron. In other words, a spin-density gradient is twice as effective at making an electric current as it is at contributing to the linear momentum.

Spin in Maxwell's equations

Motivated by the Riemann–Silberstein vector form of Maxwell's equations, Michael Berry ^[4] uses the Gordon strategy to obtain gauge-invariant expressions for the intrinsic spin angular-momentum density for solutions to Maxwell's equations.

He assumes that the solutions are monochromatic and uses the phasor expressions $\mathbf {E} =\mathbf {E} (\mathbf {r} )e^{-i\omega t}$ , $\mathbf {H} =\mathbf {H} (\mathbf {r} )e^{-i\omega t}$ . The time average of the Poynting vector momentum density is then given by

{\begin{aligned}\langle \mathbf {P} \rangle &={\frac {1}{4c^{2}}}[{\mathbf {E} }^{*}\times {\mathbf {H} }+{\mathbf {E} }\times {\mathbf {H} }^{*}]\\&={\frac {\epsilon _{0}}{4i\omega }}[{\mathbf {E} }^{*}\cdot (\nabla {\mathbf {E} })-(\nabla {\mathbf {E} }^{*})\cdot {\mathbf {E} }+\nabla \times ({\mathbf {E} }^{*}\times {\mathbf {E} })]\\&={\frac {\mu _{0}}{4i\omega }}[{\mathbf {H} }^{*}\cdot (\nabla {\mathbf {H} })-(\nabla {\mathbf {H} }^{*})\cdot {\mathbf {H} }+\nabla \times ({\mathbf {H} }^{*}\times {\mathbf {H} })].\end{aligned}}

We have used Maxwell's equations in passing from the first to the second and third lines, and in expression such as ${\mathbf {H} }^{*}\cdot (\nabla {\mathbf {H} })$ the scalar product is between the fields so that the vector character is determined by the $\nabla$ .

As

\mathbf {P} _{\text{tot}}=\mathbf {P} _{\text{free}}+{\mathbf {P} }_{\text{bound}},

and for a fluid with intrinsic angular momentum density $\mathbf {S}$ we have

\mathbf {P} _{\text{bound}}={\frac {1}{2}}\nabla \times \mathbf {S} ,

these identities suggest that the spin density can be identified as either

\mathbf {S} ={\frac {\mu _{0}}{2i\omega }}{\mathbf {H} }^{*}\times {\mathbf {H} }

or

\mathbf {S} ={\frac {\epsilon _{0}}{2i\omega }}{\mathbf {E} }^{*}\times {\mathbf {E} }.

The two decompositions coincide when the field is paraxial. They also coincide when the field is a pure helicity state – i.e. when ${\mathbf {E} }=i\sigma c{\mathbf {B} }$ where the helicity $\sigma$ takes the values $\pm 1$ for light that is right or left circularly polarized respectively. In other cases they may differ.

Related Research Articles

In particle physics, the Dirac equation is a relativistic wave equation derived by British physicist Paul Dirac in 1928. In its free form, or including electromagnetic interactions, it describes all spin-1⁄2 massive particles such as electrons and quarks for which parity is a symmetry. It is consistent with both the principles of quantum mechanics and the theory of special relativity, and was the first theory to account fully for special relativity in the context of quantum mechanics. It was validated by accounting for the fine details of the hydrogen spectrum in a completely rigorous way.

The Klein–Gordon equation is a relativistic wave equation, related to the Schrödinger equation. It is second-order in space and time and manifestly Lorentz-covariant. It is a quantized version of the relativistic energy–momentum relation. Its solutions include a quantum scalar or pseudoscalar field, a field whose quanta are spinless particles. Its theoretical relevance is similar to that of the Dirac equation. Electromagnetic interactions can be incorporated, forming the topic of scalar electrodynamics, but because common spinless particles like the pions are unstable and also experience the strong interaction the practical utility is limited.

In quantum field theory, the Dirac spinor is the spinor that describes all known fundamental particles that are fermions, with the possible exception of neutrinos. It appears in the plane-wave solution to the Dirac equation, and is a certain combination of two Weyl spinors, specifically, a bispinor that transforms "spinorially" under the action of the Lorentz group.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics. The Hamilton–Jacobi equation is particularly useful in identifying conserved quantities for mechanical systems, which may be possible even when the mechanical problem itself cannot be solved completely.

In theoretical physics, the Rarita–Schwinger equation is the relativistic field equation of spin-3/2 fermions. It is similar to the Dirac equation for spin-1/2 fermions. This equation was first introduced by William Rarita and Julian Schwinger in 1941.

In atomic physics, the electron magnetic moment, or more specifically the electron magnetic dipole moment, is the magnetic moment of an electron caused by its intrinsic properties of spin and electric charge. The value of the electron magnetic moment is approximately −9.284764×10⁻²⁴ J/T. The electron magnetic moment has been measured to an accuracy of 7.6 parts in 10¹³.

In differential geometry, the four-gradient $is the four-vector analogue of the gradient from vector calculus.$

In physics, the Majorana equation is a relativistic wave equation. It is named after the Italian physicist Ettore Majorana, who proposed it in 1937 as a means of describing fermions that are their own antiparticle. Particles corresponding to this equation are termed Majorana particles, although that term now has a more expansive meaning, referring to any fermionic particle that is its own anti-particle.

In differential geometry and mathematical physics, a spin connection is a connection on a spinor bundle. It is induced, in a canonical manner, from the affine connection. It can also be regarded as the gauge field generated by local Lorentz transformations. In some canonical formulations of general relativity, a spin connection is defined on spatial slices and can also be regarded as the gauge field generated by local rotations.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— $in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.$

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

In quantum mechanics, the Pauli equation or Schrödinger–Pauli equation is the formulation of the Schrödinger equation for spin-½ particles, which takes into account the interaction of the particle's spin with an external electromagnetic field. It is the non-relativistic limit of the Dirac equation and can be used where particles are moving at speeds much less than the speed of light, so that relativistic effects can be neglected. It was formulated by Wolfgang Pauli in 1927.

In mathematical physics, spacetime algebra (STA) is a name for the Clifford algebra Cl_1,3(R), or equivalently the geometric algebra G(M⁴). According to David Hestenes, spacetime algebra can be particularly closely associated with the geometry of special relativity and relativistic spacetime.

The intent of this article is to highlight the important points of the derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids.

In physics, and specifically in quantum field theory, a bispinor, is a mathematical construction that is used to describe some of the fundamental particles of nature, including quarks and electrons. It is a specific embodiment of a spinor, specifically constructed so that it is consistent with the requirements of special relativity. Bispinors transform in a certain "spinorial" fashion under the action of the Lorentz group, which describes the symmetries of Minkowski spacetime. They occur in the relativistic spin-1/2 wave function solutions to the Dirac equation.

In physics, relativistic quantum mechanics (RQM) is any Poincaré covariant formulation of quantum mechanics (QM). This theory is applicable to massive particles propagating at all velocities up to those comparable to the speed of light c, and can accommodate massless particles. The theory has application in high energy physics, particle physics and accelerator physics, as well as atomic physics, chemistry and condensed matter physics. Non-relativistic quantum mechanics refers to the mathematical formulation of quantum mechanics applied in the context of Galilean relativity, more specifically quantizing the equations of classical mechanics by replacing dynamical variables by operators. Relativistic quantum mechanics (RQM) is quantum mechanics applied with special relativity. Although the earlier formulations, like the Schrödinger picture and Heisenberg picture were originally formulated in a non-relativistic background, a few of them also work with special relativity.

In mathematical physics, the Belinfante–Rosenfeld tensor is a modification of the energy–momentum tensor that is constructed from the canonical energy–momentum tensor and the spin current so as to be symmetric yet still conserved.

In physics, particularly in quantum field theory, the Weyl equation is a relativistic wave equation for describing massless spin-1/2 particles called Weyl fermions. The equation is named after Hermann Weyl. The Weyl fermions are one of the three possible types of elementary fermions, the other two being the Dirac and the Majorana fermions.

In probability theory and statistics, the normal-inverse-Wishart distribution is a multivariate four-parameter family of continuous probability distributions. It is the conjugate prior of a multivariate normal distribution with unknown mean and covariance matrix.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

References

↑ W. Gordon (1928). "Der Strom der Diracschen Elektronentheorie". Z. Phys. 50: 630–632. Bibcode:1928ZPhy...50..630G. doi:10.1007/BF01327881.
↑ M.Stone (2015). "Berry phase and anomalous velocity of Weyl fermions and Maxwell photons". International Journal of Modern Physics B. 30: 1550249. arXiv: 1507.01807 . doi:10.1142/S0217979215502495.
↑ D.T.Son, N.Yamamoto (2013). "Kinetic theory with Berry curvature from quantum field theories". Physical Review D. 87: 085016. arXiv: 1210.8158 . Bibcode:2013PhRvD..87h5016S. doi:10.1103/PhysRevD.87.085016.
↑ M.V.Berry (2009). "Optical currents". J. Opt. A. 11: 094001 (12 pages). Bibcode:2009JOptA..11i4001B. doi:10.1088/1464-4258/11/9/094001.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] W. Gordon (1928). "Der Strom der Diracschen Elektronentheorie". Z. Phys. 50: 630–632. Bibcode:1928ZPhy...50..630G. doi:10.1007/BF01327881.

[stone-2] M.Stone (2015). "Berry phase and anomalous velocity of Weyl fermions and Maxwell photons". International Journal of Modern Physics B. 30: 1550249. arXiv: 1507.01807 . doi:10.1142/S0217979215502495.

[son-3] D.T.Son, N.Yamamoto (2013). "Kinetic theory with Berry curvature from quantum field theories". Physical Review D. 87: 085016. arXiv: 1210.8158 . Bibcode:2013PhRvD..87h5016S. doi:10.1103/PhysRevD.87.085016.

[berry-4] M.V.Berry (2009). "Optical currents". J. Opt. A. 11: 094001 (12 pages). Bibcode:2009JOptA..11i4001B. doi:10.1088/1464-4258/11/9/094001.

[1]

[2]

[3]

[4]