Nonlinear Dirac Equations
We construct nonlinear extensions of Dirac's relativistic electron equation that preserve its other desirable properties such as locality, separability, conservation of probability and Poincaré invariance. We determine the constraints that the nonlinear term must obey and classify the resultant...
Saved in:
| Published in: | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Date: | 2009 |
| Main Authors: | , |
| Format: | Article |
| Language: | English |
| Published: |
Інститут математики НАН України
2009
|
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/149176 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Nonlinear Dirac Equations / Wei Khim Ng, Rajesh R. Parwani // Symmetry, Integrability and Geometry: Methods and Applications. — 2009. — Т. 5. — Бібліогр.: 46 назв. — англ. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1860030789034442752 |
|---|---|
| author | Wei Khim Ng Rajesh R. Parwani |
| author_facet | Wei Khim Ng Rajesh R. Parwani |
| citation_txt | Nonlinear Dirac Equations / Wei Khim Ng, Rajesh R. Parwani // Symmetry, Integrability and Geometry: Methods and Applications. — 2009. — Т. 5. — Бібліогр.: 46 назв. — англ. |
| collection | DSpace DC |
| container_title | Symmetry, Integrability and Geometry: Methods and Applications |
| description | We construct nonlinear extensions of Dirac's relativistic electron equation that preserve its other desirable properties such as locality, separability, conservation of probability and Poincaré invariance. We determine the constraints that the nonlinear term must obey and classify the resultant non-polynomial nonlinearities in a double expansion in the degree of nonlinearity and number of derivatives. We give explicit examples of such nonlinear equations, studying their discrete symmetries and other properties. Motivated by some previously suggested applications we then consider nonlinear terms that simultaneously violate Lorentz covariance and again study various explicit examples. We contrast our equations and construction procedure with others in the literature and also show that our equations are not gauge equivalent to the linear Dirac equation. Finally we outline various physical applications for these equations.
|
| first_indexed | 2025-12-07T16:52:12Z |
| format | Article |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 5 (2009), 023, 20 pages
Nonlinear Dirac Equations
Wei Khim NG and Rajesh R. PARWANI
Department of Physics, National University of Singapore, Kent Ridge, Singapore
E-mail: g0500437@nus.edu.sg, parwani@nus.edu.sg
URL: http://staff.science.nus.edu.sg/∼parwani/
Received November 12, 2008, in final form February 23, 2009; Published online February 27, 2009
doi:10.3842/SIGMA.2009.023
Abstract. We construct nonlinear extensions of Dirac’s relativistic electron equation that
preserve its other desirable properties such as locality, separability, conservation of proba-
bility and Poincaré invariance. We determine the constraints that the nonlinear term must
obey and classify the resultant non-polynomial nonlinearities in a double expansion in the
degree of nonlinearity and number of derivatives. We give explicit examples of such nonlinear
equations, studying their discrete symmetries and other properties. Motivated by some pre-
viously suggested applications we then consider nonlinear terms that simultaneously violate
Lorentz covariance and again study various explicit examples. We contrast our equations
and construction procedure with others in the literature and also show that our equations
are not gauge equivalent to the linear Dirac equation. Finally we outline various physical
applications for these equations.
Key words: nonlinear Dirac equation; Lorentz violation
2000 Mathematics Subject Classification: 81P05; 81Q99; 83A05
1 Introduction
When Schrödinger obtained his wave equation to realise de Broglie’s speculation about the wave
nature of particles, he used a number of heuristic arguments and assumed the simplest possibility,
that of linearity of the equation [1]. Fortunately that assumption led to very good agreement with
experiment and till today no deviations from quantum linearity have been detected eventhough
a few low energy experiments have attempted to observe them [2, 3, 4, 5]. Currently the main
interest in nonlinear Schrödinger equations is that they appear, in form, as approximations in
optics and condensed matter [6, 7].
When Dirac generalised Schrödinger’s equation to the relativistic domain, he too kept lineari-
ty. Nonlinear versions of Dirac’s equation have been studied for various purposes since then.
Heisenberg’s proposal [8] was in the context of field theory and was motivated by the question of
mass. In the quantum mechanical context, nonlinear Dirac equations have been used as effective
theories in atomic, nuclear and gravitational physics [9, 10, 11]. Some of the simpler versions
have been analysed rigorously [12].
Although there is as yet no evidence for fundamental quantum nonlinearities, their absence
is seen as a puzzle by several authors and requires an understanding [13, 14, 15, 16, 17]. Based
on an extrapolation of some information theoretic arguments at the non-relativistic level, it was
proposed in [18] that perhaps quantum linearity might be intimately tied to Lorentz invariance
and that the possible violation of the latter at a fundamental level might lead to quantum
nonlinearities. If true, then perhaps the appropriate regime to seek such inter-related violations
would be at high energies or at very short distances.
Since quantum nonlinearities, if they exist, must be very small, the best place to search for
them is where they might show up at leading order, not masked by other corrections. Thus one
hopes to detect the nonlinearities at the quantum mechanical level, rather than as supplements to
mailto:g0500437@nus.edu.sg
mailto:parwani@nus.edu.sg
http://staff.science.nus.edu.sg/~parwani/
http://dx.doi.org/10.3842/SIGMA.2009.023
2 W.K. Ng and R.R. Parwani
loop effects in field theory. Neutrinos are therefore an ideal probe of such potential nonlinearities
as they are weakly interacting and so not affected much by field theory corrections.
Indeed, neutrino oscillations were suggested in [18] as one place where quantum nonlinearities
might be relevant and a heuristic study was conducted using a provisional nonlinear Dirac
equation. That equation was very complicated and it did not conserve probability.
In this paper we discuss Dirac equations, at the quantum mechanical level, which preserve
all the others desirable features such as conservation of probability. We intend to use these
equations to study not just neutrino oscillations but also various other high-energy phenomenon
which are briefly discussed in the last section.
However it is possible that our equations might also be relevant as approximate equations,
for use either in particle physics or condensed matter physics, and we discuss this also in the
concluding section.
As there are various obstacles to generalising the non-relativistic information theory approach
of [18] to the relativistic domain, we proceed in a different manner here. We write the nonlinear
equation as
(iγµ∂µ −m+ F )ψ = 0, (1.1)
where F is a function of the wavefunction ψ, its adjoint and their derivatives1. We begin
by requiring, just as for F = 0, that equation (1.1) be local, Poincaré covariant, conserves
probability and is separable for multi-particle states. The constraints on F are then solved in an
expansion procedure to be detailed in Section 2.5.1. That is, we implement a systematic scheme
to construct a large class of nonlinear extensions of the Dirac equation.
The constraints we adopt are similar to those used in understanding non-relativistic quantum
theory in [19, 20]. There it was deduced that the Schrödinger equation is the unique single
universal parameter (~) extension of classical ensemble dynamics. Although the speed of light, c,
is a universal parameter for relativistic dynamics, it already appears at the classical level and
plays the role of converting the dimensions of space to those of time. One expects that further
extensions of quantum theory either at the non-relativistic or relativistic level would involve
other universal parameters, for example a universal length.
Our approach and most of our results differ from previous constructs of nonlinear Dirac
equations in the literature. Most studies [9, 10, 11] do not impose separability, which is a strong
constraint that leads to non-polynomiality of F . In [21] separability was imposed in a somewhat
different manner from what we do here, but more importantly the authors of [21] only considered
nonlinear Dirac equations that are obtained from the linear Dirac equation through a process
of gauge-completion: thus their class of equations is more restrictive than ours. Some further
contrasts of our procedure compared to others [22, 23] is that we allow derivatives of the wave-
function in F , and also study nonlinearities which violate one or more of the discrete P, C, T
symmetries as such cases are expected to be phenomenologically relevant.
Furthermore, proceeding with the suggestion of [18], we also construct versions of (1.1) that
are simultaneously Lorentz violating and nonlinear: such equations have also not been studied
before in general; however we note that one example of such an equation, without derivatives in
the nonlinearity, has been studied in [24, 25], motivated by anisotropic space-times [26].
The rest of the paper is structured as follows: In Section 2 we discuss and make explicit
the various constraints on the nonlinear term F ; we note that the class of nonlinearities we
consider can also be motivated without imposing separability and so are potentially useful also
as effective equations at low energies. The simplest examples of such equations are discussed in
Section 3 followed by their plane-wave solutions and the corresponding dispersion relations in
1But we do not consider F ’s that have free derivatives acting to the right on the final ψ of the equation (1.1).
So our nonlinearity is a matrix in spinor space with spacetime dependent coefficients.
Nonlinear Dirac Equations 3
Section 4. In Section 5 we study examples of F that simultaneously violate Lorentz covariance.
In Section 6 we illustrate more complicated examples of the nonlinear equations and also discuss
the alternative approach whereby the nonlinear equations are obtained from a Lagrangian. In
Section 7 we explain how to distinguish our nonlinearities from those that may be obtained from
the linear equation through a nonlinear gauge transformation. A summary and outlook is in
Section 8.
Although the evolution equation (1.1) has been modified, we keep the usual kinematical
structure of quantum mechanics; some arguments, that fundamental nonlinear quantum theories
are intrinsically pathological, are discussed in the final section. The conventions we use are
similar to those in the textbook [27]; unless stated, our discussion is representation independent.
Although we work in 3 + 1 dimensional flat spacetime with metric gµν = (1,−1,−1,−1), some
effects of gravity could possibly be encoded in an effective nonlinearity; we do not study in
this paper explicit couplings to gravity though this might yield some interesting consequences
as seen for the linear Dirac equation [28].
2 Constraints
The usual, linear, quantum-mechanical (“first quantised”) Dirac equation has many appealing
properties which we will mostly preserve so as to achieve a minimal deformation. Later in
Sections 5 and 6, we discuss the possibility of further extensions motivated by physical considera-
tions.
We now list and explain the various constraints that we are going to impose on the nonlinear
Dirac equation and hence on F in (1.1).
2.1 Locality
We continue to assume that physics, as described by the wavefunction ψ, is accurately captured
by a local evolution equation: that is we require F to depend only on ψ, its conjugate and their
derivatives all evaluated at a single point x. Note that F below is in general a matrix in spinor
space though later we will specialise to various cases, such as F proportional to the identity
matrix.
Notice that we demand locality of the equations of motion rather than of a Lagrangian.
This means that some of our equations might not be obtainable from a local Lagrangian. One
could of course implement a construction procedure similar to that described below at the
local Lagrangian level: we illustrate this in Section 6 and discuss the relative advantages and
disadvantages.
2.2 Poincaré invariance
Under the Poincaré transformation x′ = Λx + a the linear Dirac equation is covariant if the
wavefunction transforms as [27]
ψ′(x′) = S(Λ)ψ(x) = ψ
(
Λ−1(x′ − a)
)
,
where S−1(Λ)γµΛν
µS(Λ) = γν . Explicitly we have S(Λ) = exp(− i
4σαβω
αβ), with ωαβ the
transformation parameters. If we demand that the nonlinear equation (1.1) be covariant under
the same transformations then we obtain the following constraint,
S−1(Λ)F ′S(Λ) = F,
where F ′ is the Poincaré transformed F ; recall that F is a function depending on ψ̄, ψ and their
derivatives.
4 W.K. Ng and R.R. Parwani
2.3 Hermiticity
In quantum mechanics we usually require the Hamiltonian to be Hermitian so as to guarantee
reality of eigenvalues. Rewriting the nonlinear Dirac equation in Hamiltonian form we have,
i
∂
∂t
ψ = (HD − βF )ψ,
where β = γ0 and HD is the linear Dirac Hamiltonian. Since H†
D = HD, thus we also impose2
γ0F †γ0 = F. (2.1)
2.3.1 Current conservation
In terms of the familiar adjoint ψ̄ = ψ†γ0, the linear Dirac equation has the conserved current
jµ = ψ̄γµψ, (2.2)
which allows ψ†ψ to be interpreted as a probability density. The divergence of the same expres-
sion (2.2) in the nonlinear theory is
∂µj
µ = ψ̄
(
iF − iγ0F †γ0
)
ψ, (2.3)
which vanishes due to the Hermiticity condition (2.1).
Thus requiring Hermiticity of the Hamiltonian also ensures conservation of (2.2). On the
other hand, in some future applications, we may want to consider non-Hermitian Hamiltonians
that model open systems. Then the right-hand-side of (2.3) can be used to measure leakage
from the system.
2.3.2 Chiral current
For completeness we also discuss the chiral current, for which the expression in the linear theory
is jµ
5 = ψ̄γµγ5ψ. Using the nonlinear equations of motion, we obtain
∂µj
µ
5 = −iψ̄
(
γ0F †γ0γ5 + γ5F
)
ψ + 2imψ̄γ5ψ.
For the usual chiral current to be conserved in the massless, m → 0, limit of the nonlinear
equation, we require
γ5F + γ0F †γ0γ5 = 0,
which, on using the Hermiticity condition (2.1), simplifies to
{F, γ5} = 0.
2.4 Universality
The usual Dirac equation has the property, as all linear equations do, that it is invariant under
a rescaling of the wavefunction, ψ → λψ. In quantum mechanics such a condition allows
solutions of the equation to be freely normalised, which is not only convenient but also sometimes
demanded for an interpretation of measurements [13, 14, 15, 16].
We would like our nonlinear generalisation to preserve the same scale-invariance property,
which one may motivate with alternative reasoning as follows. We desire equations that are
2Recall, we are adopting the standard kinematical structure of quantum mechanics, in particular the standard
inner product. See also the first footnote.
Nonlinear Dirac Equations 5
as universal as possible. So, for example, the equation should have the same form whether it
describes a single particle or a system of particles. More specifically, the parameters describing
the strength of the nonlinearity F should not be dependent on the number of particles in the
system, just as Planck’s constant ~ is universal in the multiparticle Schrödinger equation.
If ψ represents the wavefunction for a N -particle state, then the normalisation of probability
implies that the dimension of ψ depends on N , just as in the non-relativistic case [19, 20], and
so the dimension of F would then be N dependent in general. We can avoid this conclusion by
requiring that F have the above-mentioned scaling property
F (λψ) = F (ψ), (2.4)
where we mean that the wavefunction and its conjugate are all scaled by the same factor λ on
the left-hand-side. Equation (2.4) implies that F must be non-polynomial,
F ∼ F (A/B), (2.5)
where A, B have equal factors of the wavefunction.
2.5 Separability
The usual Dirac equation may be used to describe a collection of particles and is separable for
independent subsystems. It seems useful to have this separability property also for our nonlinear
generalisation. However as we will explain in a later section, one may omit the separability
constraint in favour of other arguments which result in similar forms for the eventual F ’s, and
those forms anyway become separable with a suitable interpretation of the multiparticle states.
Thus with the same structure for F we can use the equation for fundamental, phenomenological
or effective dynamics.
Let us review separability first for the linear Dirac equation so as to motivate suitable defi-
nitions of ψ̄ and jµ for many-body systems. In the multi-time formalism [29, 30, 31], which
preserves manifest Poincaré invariance, the many-body linear Dirac equation for non-interacting
particles may be written as∑
s
(iγµ
s ∂µ,s −ms)ψ = 0,
where
ψ = ψ1 ⊗ ψ2 ⊗ · · · ⊗ ψs︸︷︷︸
sth site
⊗ · · · ,
γµ
s = 1⊗ 1⊗ · · · ⊗ γµ︸︷︷︸
sth site
⊗1⊗ · · · ,
ms = 1⊗ 1⊗ · · · ⊗ m(s)︸︷︷︸
sth site
⊗1⊗ · · ·
and s labels the particle. Consider explicitly the two-particles case,[(
iγµ∂µ,1 −m(1)
)
ψ1
]
⊗ ψ2 + ψ1 ⊗
[(
iγµ∂µ,2 −m(2)
)
ψ2
]
= 0. (2.6)
Let φ1 and φ2 be arbitrary single particle wavefunctions for the two independent variables 1, 2.
Then multiplying by φ̄1 ⊗ φ̄2/(φ̄1ψ1)(φ̄2ψ2) on the left of (2.6), we have
φ̄1
(
iγµ∂µ,1 −m(1)
)
ψ1
φ̄1ψ1
⊗ 1 + 1⊗
φ̄2
(
iγµ∂µ,2 −m(2)
)
ψ2
φ̄2ψ2
= 0. (2.7)
6 W.K. Ng and R.R. Parwani
The result is clearly separable in that solutions of the individual single particle Dirac equations
satisfy the two-particle equation and vice-versa.
Furthermore it is easy to show that if ψ̄ for a many-body system is defined to be ψ̄ =
ψ̄1 ⊗ ψ̄2 ⊗ · · · ⊗ ψ̄s ⊗ · · · then the two-particle adjoint equation that follows from (2.6) will have
the same form as the one-particle equation, and since form-invariance is in the spirit of the
universality criteria of Section 2.4, this justifies our definition.
Now consider, as an example, the expression for the multi-particle current jµ. Multiply (2.6)
from the left by ψ̄1 ⊗ ψ̄2, multiply the adjoint of (2.6) from the right by ψ1 ⊗ ψ2, and take the
difference to get[
∂µ,1j
µ
(1)
]
⊗ ψ̄2ψ2 + ψ̄1ψ1 ⊗
[
∂µ,2j
µ
(2)
]
= 0 (2.8)
⇒
2∑
s=1
∂µ,sj
µ
s = 0,
where the current is defined to be
jµ
s = ψ̄1ψ1 ⊗ ψ̄2ψ2 ⊗ · · · ⊗ jµ
(s)︸︷︷︸
sth site
⊗ψ̄s+1ψs+1 · · · .
Multiplying (2.8) by
(
ψ̄1ψ1 ⊗ ψ̄2ψ2
)−1 gives
∂µ,1j
µ
(1)
ψ̄1ψ1
⊗ 1 + 1⊗
∂µ,2j
µ
(2)
ψ̄2ψ2
= 0.
Thus conservation of individual currents implies the conservation of the two-particle current and
vice-versa.
Similarly the definition
γ5
s = 1⊗ 1⊗ · · · ⊗ γ5︸︷︷︸
sth site
⊗1⊗ · · · (2.9)
allows the multiparticle chiral current to be defined and in the massless limit conservation of
individual chiral currents implies the conservation of the two-particle chiral current and vice-
versa.
2.5.1 Structure of F
We would like our nonlinear equation to be separable in this minimal sense: for a wavefunction
which is the product of two independent states, the composite equation should decompose into
two independent equations3. Looking at the expressions (2.6), (2.7) we see that for the nonlinear
equation (1.1) to be separable as such, we require F to decompose as
F (ψ1 ⊗ ψ2) = F (ψ1)⊗ 1 + 1⊗ F (ψ2)
for a state made up of two independent particles (constraints of this type have been studied
before for non-relativistic systems in [33]). Equation (2.5) and the examples above suggest that
this can be achieved if we have the structure
F
(
N
D
)
∼ N
D
→ N1
D1
⊗ 1 + 1⊗ N2
D2
.
Thus for a product state we require N → N1 ⊗D2 +D1 ⊗N2 while D → D1 ⊗D2.
3We note that other implementations of separability might lead to more constraints, see for example [32].
Nonlinear Dirac Equations 7
Requiring N and D to be separately Poincaré invariant, we see that the only functional of ψ
that would decompose as required for D is ψ̄ψ and powers thereof. Thus our nonlinear term F
can be a sum of terms of the form
N
(
ψ̄, ψ
)(
ψ̄ψ
)n , (2.10)
subject to the other constraints that have yet to be imposed.
Our deduction of (2.10) has been somewhat heuristic and so the reader may prefer to think
of it as an ansatz within which we discuss our equations.
As mentioned earlier, the separability condition is appropriate for fundamental equations that
describe an arbitrary collection of particles. However if the nonlinearities are an approximate
description of an underlying dynamics, as effective equations attempt to do, then the universality
and separability arguments do not seem appropriate. However even then one may motivate the
structure (2.10) as follows. Generally, for slowly varying fields, one may perform a gradient
expansion for F when seeking local equations,
F ∼ N0
D0
+
N1
D1
+
N2
D2
+ · · ·+ Ni
Di
+ · · · , (2.11)
where the Ni’s depend on the wavefunction and contain exactly i derivatives. The Di’s also
depend on the wavefunction but do not contain any derivatives.
Now in most nonlinear Scrödinger or Dirac equations the nonlinear terms break the scale
invariance, ψ → λψ, present in the linear theory. That is, typically the nonlinearities make the
equations sensitive to the amplitude of the fields thus giving rise to very interesting phenomena.
However it is possible to have nonlinearities that preserve the scale invariance of the linear theory
and though the effects are then likely to be milder, they can still lead be novel and interesting
effects [34, 35]. So if we focus on such “soft” nonlinearities, and also impose Lorentz invariance,
then (2.11) is included in the form (2.10). Indeed, as we shall verify later, even without imposing
separability at the outset, separability of the resultant structures appears to be possible with
consistent definitions of the multi-particle states.
In summary, we will discuss in this paper the class of nonlinearities of the form (2.10) by
looking at several cases corresponding to a specific degree of nonlinearity, n = 1, 2, . . . , and
a derivative expansion of the numerator.
We remark that the scale-invariant nonlinearities (2.10) we introduce here might also be
interesting for future quantum field theory investigations: these nonlinearities correspond to
Lagrangians that are still naively power-counting renormalisable.
2.6 Discrete symmetries
The Standard Model of particle physics encodes both parity and CP violation as these are em-
pirically observed facts. Thus in our nonlinear equation we find it interesting to allow violation
of individual symmetries. However in line with general theorems [27, 36] on local, Hermitian,
Lorentz covariant theories, we do find by explicit verification that our specific examples preserve
the combined PCT invariance although we do not impose it.
The discrete symmetry operators are the same as in the linear theory [27], and they place
constraints on the nonlinear term F so that the nonlinear equation (1.1) is form invariant (similar
to the discussion in Section 2.2).
Ignoring unobservable phases, the representation independent parity operator is P̂ = γ0 and
parity invariance requires
P̂−1FP P̂ ≡ F,
8 W.K. Ng and R.R. Parwani
where FP is the parity transformed F . Charge conjugation invariance is achieved if
Ĉ−1FC Ĉ ≡ F ∗,
where FC is the charge conjugated F and Ĉ = iγ2 in the Dirac–Pauli representation. The
time-reversal invariance constraint on F is
T̂ −1FT T̂ ≡ F ∗,
FT being the time reversed F and T̂ = iγ1γ3 in the Dirac–Pauli representation4.
Under the combined PCT transformation, Θ̂, the nonlinear Dirac equation in invariant if
Θ̂−1FΘΘ̂ ≡ F,
where FΘ is the PCT transformed F . The representation independent form for Θ̂ is proportional
to γ5.
3 Explicit examples of nonlinear equations with F ∝ I, n = 1
We found earlier in Section 2.5.1 that F has the form
N
(
ψ̄, ψ
)(
ψ̄ψ
)n , (3.1)
where the number of factors of the wavefunction in the numerator is 2n.
In the absence of other dynamical fields, Poincaré invariance requires spacetime indices of
matrices like γµ to be contracted among themselves or with derivatives ∂µ. We will assume that
the spinor indices of ψ and ψ̄ are contracted in the natural way with ψ̄ acting like a row vector
and ψ a column vector, for example N ∼ Aψ̄BψC where A, B, C are matrices in spinor space.
In this Section we restrict the explicit discussion to the important case where F is proportional
to the identity matrix I in spinor space,
F = fI (3.2)
and so the nonlinearity f may be thought of as a spacetime dependent mass. This choice is
motivated by our interest in neutrino oscillations. We also consider here only the lowest order
of nonlinearity, n = 1. In Section 6, we discuss some other types of F .
Current conservation for the case (3.2) simply amounts to the statement that f is a real
function of the wavefunction,
f = f∗.
3.1 No derivatives
In the absence of derivatives, the most general structure of the nonlinear term with F ∝ I and
n = 1 is given by
F =
ψ̄Aψ
ψ̄ψ
, (3.3)
4We remind the reader that we are treating the wavefunction as a classical object rather than as a Grassmann
variable. Thus the denominator of (2.10) obtains a negative sign when performing a transpose operation, such as
occurs in charge conjugation. Such negative signs mutually cancel for our scale-invariant nonlinearities (2.10).
Nonlinear Dirac Equations 9
where A is a matrix. In the absence of other fields which carry spacetime indices we must
therefore have
A = aI + ibγ5,
where a, b are constants. The a term is clearly equivalent to a mass term in the linear equation
and so may be ignored in the following discussion. Notice that the form A = ibγ5 in (3.3),
which is a consequence of Lorentz invariance, also automatically satisfies the PCT invariance
condition.
As for individual discrete symmetries, using the equations of Section 2.6, we see that the
term with b 6= 0 preserves C invariance but breaks parity. Time-reversal invariance requires b
to be purely imaginary, which conflicts with the requirement from current conservation which
requires b to be real.
We thus conclude that our simplest nonlinear equation, with F ∝ I and n = 1,
F = iε
ψ̄γ5ψ
ψ̄ψ
.
unavoidably breaks P and CP, something that is surely intriguing from the perspective of particle
physics phenomenology. We have indicated the small nonlinearity parameter by ε.
Note that the multiparticle version of the above equation is separable, so that does not impose
additional constraints. Nonlinear Dirac equations without derivatives in the nonlinear part have
been studied in [22, 23] and (3.3) is a special case of the equations studied there.
3.1.1 Lorentz vs PCT invariance
Let us discuss the situation whereby the PCT invariance is imposed on (3.3) first. Then we
find, using Θ̂ ∝ γ5, that we require
[A, γ5] = 0,
which is satisfied if A has the form
A = aI + bγ5 + cµνσµν .
If there are no other dynamical fields other than the wavefunction, then cµν can only be a con-
stant background field, thus explicitly breaking Lorentz invariance. Indeed, explicitly imple-
menting Lorentz invariance of (3.3) gives
S(Λ)−1AS(Λ) ≡ A,
which for the infinitesimal case gives [A, σαβ ] = 0. This only allows
A = aI + bγ5
as we argued earlier.
In other words, we can have PCT invariance even if we give up Lorentz invariance, which
again is consistent with general results in the literature [27, 36].
3.2 One derivative
The most general form of F is now given by the linear combination of the following two terms,(
∂µψ̄
)
AγµBψ
ψ̄ψ
,
ψ̄CγµD∂µψ
ψ̄ψ
.
10 W.K. Ng and R.R. Parwani
As in the no derivative case, Lorentz covariance requires that both A,B be proportional to
a linear combination of I, γ5 and so we may write
F =
(
∂µψ̄
)
(aI + ibγ5)γµψ
ψ̄ψ
+
ψ̄(cI − idγ5)γµ∂µψ
ψ̄ψ
,
a result which also satisfies PCT invariance. Hermiticity of this F , and hence current conserva-
tion, is satisfied if we have c = a∗ and d = b∗. Clearly parity invariance is violated if b 6= 0; in
that case C invariance requires b to be purely imaginary while T invariance requires b to be real.
The constant a is not constrained by parity but both C and T invariance separately require a
to be purely imaginary.
Let us consider the special case where each of the discrete symmetries is individually pre-
served: b = 0 and a = iε with ε a real parameter that controls the strength of the nonlinearity.
Then we may write, using explicitly the on-shell current conservation condition,
F = 2iε
(
∂µψ̄
)
γµψ
ψ̄ψ
= −2iε
ψ̄γµ∂µψ
ψ̄ψ
. (3.4)
For ε small, one may simplify F in (3.4) by solving the nonlinear Dirac equation (1.1) itera-
tively. To leading order (iγµ∂µ −m)ψ = 0 which when used in F gives F = −2εm. Thus to
leading order in ε the nonlinearity (3.4) is just a mass shift.
We remark that just as in the no derivative case, we could have imposed PCT invariance
first and obtained cases which violate Lorentz covariance. However we defer further discussion
of Lorentz violating cases to a later section.
3.3 Two derivatives
There are well-known problems in constructing Lorentz covariant higher-derivative first-quanti-
sed theories. Consider a normalised state,
1 =
∫
d3xψ†ψ.
Applying ∂
∂t to both sides gives
0 =
∫
d3x(ψ̇†ψ + ψ†ψ̇).
Now, if the evolution is second-order in time, then one can specify ψ(0, x) and ψ̇(0, x) indepen-
dently and that would mean that the right-hand-side of the above equation need not be zero in
general, leading to a contradiction.
However, in our nonlinear equations, Hermiticity and hence current conservation are ensured
by construction and so the above-mentioned problem does not occur. This of course does not
guarantee that all other physical quantities will be well-behaved, but it is plausible that that is
the case if the higher-order terms are treated perturbatively.
The general structure of the two-derivative nonlinear term, F ∝ I, without embedded γ
matrices is
F =
a
(
∂µ∂
µψ̄
)
ψ + bψ̄∂µ∂
µψ + c
(
∂µψ̄
)
(∂µψ)
ψ̄ψ
.
Each numerator/denominator term is separately Poincaré and PCT invariant. However while
each term is also separately parity invariant, C or T invariance requires all the coefficients a, b, c
to be real.
Current conservation, F = F † implies that b = a∗ and c = c∗. Thus we conclude that for a
not real, both C and T (or CP) are violated.
Nonlinear Dirac Equations 11
4 Plane-wave solutions and dispersion relations
We wish to construct plane-wave solutions to the nonlinear equations of the previous section.
As in the case for the linear theory, we require the solutions to be simultaneous eigenstates of
momentum and energy. Let us clarify what this means in the nonlinear theory.
Although we allow the equations to be nonlinear, we keep the fundamental commutation
relation between the position and momentum operators, [x̂, p̂] = i~. Thus in the Schrödinger
representation we have p̂ = −i~∂ and the momentum eigenvalue is given by p̂ψp = pψp. Likewise
the energy-eigenvalue equation is given by i~∂tψE = EψE .
With Lorentz covariance preserved, the method to find plane-wave solutions is similar to the
linear case. We seek solutions of the form5
ψ(x, t) = e−ik.xu(k) (4.1)
with kµ a four vector.
The dispersion relations will be covariantly modified from that of the linear theory. Consider
the nonlinear Dirac equation,
i∂tψ = [iα · ∂ + βm− βF (ψ)]ψ (4.2)
for the case F = fI where αi = γ0γi. Substituting the plane wave ansatz into the above
equation, squaring this and re-arranging gives
ψ†k2ψ = ψ† [m− f(kµ)]2 ψ. (4.3)
Thus we have,
k2 =
[
m− f
(
k2
)]2
. (4.4)
(Since equation (4.3) is covariant, then f must be also covariant.)
The solution of (4.4) requires the explicit form for f , the nonlinear term. It may also require
the explicit form for the plane wave solutions which we discuss next. Note that from the above
expression, one may view the effect of the nonlinearity for plane wave states as giving rise to an
effective mass.
Assume m 6= 0. Then in the rest frame we have from (4.1), (4.2),
Mu = [βm− βF (u)]u, (4.5)
where the rest energy has been labelled by M > 0.
For the case F ∝ I, the rest frame Hamiltonian is therefore proportional to γ0 = β and the
eigenstates are as in the linear theory [27]. These can then be boosted as usual to obtain the
general solutions. The net result is similar to the usual spinor solutions of the linear theory but
with the effective mass M in place of the bare mass m,
E2 = k2 +M2.
The expression for M in terms of m and the nonlinear parameters can be determined by sub-
stituting the rest frame spinors into (4.5).
5We have set ~ = 1.
12 W.K. Ng and R.R. Parwani
4.1 Perturbative method
The procedure of boosting rest frame solutions is valid if Lorentz invariance is a symmetry of
the theory. If we relax the constraint of Lorentz invariance6, we will not be able to use this
method to find the energy dispersion relations. Thus we will now introduce a method to obtain
the energy dispersion relation, to leading order in the nonlinearity, even if we do not know the
exact plane wave solutions to the theory.
From (4.3), we have
E2 = k2 +m2 − 2mf + f2.
Since the nonlinear term will contain a small nonlinearity parameter ε, we can explicitly factor
it out. That is, f = εf̃ . Then to leading order, we have
E2 = k2 +m2 − 2εmf̃.
Now we assume the following,
kµ = k(0)µ +O(ε), u(k) = u(0)
(
k(0)
)
+O(ε),
where k(0) and u(0)(k(0)) are the usual 4-momentum and u’s for the linear theory. Thus to
leading order in ε we have
E2 =
(
k2 +m2
)
− 2εmf̃
(
u(0)
(
k(0)
))
+O
(
ε2
)
=
(
k2 +m2
)
− 2εmf̃
(
u(0)(k)
)
+O
(
ε2
)
.(4.6)
Note that in the last step, we have replaced k(0) by k. This is alright because we are dropping
terms that are order ε2 or higher.
The perturbative method allows us to find corrections to the linear theorie’s energy dispersion
relation. We only need to substitute linear plane wave solutions into the nonlinear term. Note
that the above method works only for the massive theory. If we consider the massless limit then
we might need to keep terms that are of order ε2.
4.2 Example
We look at an explicit example corresponding to F ∝ I and n = 1 with two derivatives, and
obtain the corresponding expression for the effective mass M for plane wave states. Although
one can work covariantly with the expressions (4.4), it is faster to work in the rest frame, that
is by using (4.5).
Consider the nonlinear term when each of the discrete symmetries is preserved,
F =
a∂µ∂µ
(
ψ̄ψ
)
+ (c− 2a)
(
∂µψ̄
)
(∂µψ)
ψ̄ψ
.
Substituting the plane wave solution, the first term drops out leaving
F = (c− 2a)M2 ≡ εM2.
Thus
M2 =
(
m− εM2
)2
. (4.7)
6Here we refer to violation of particle Lorentz invariance while keeping observer Lorentz invariance. This can
be done by introducing background fields, see Section 5.
Nonlinear Dirac Equations 13
Taking the square root and solving we get
M =
∓1±
√
1 + 4εm
2ε
. (4.8)
For the rest energy to be real, we need ε ≥ − 1
4m . Let us consider the case where ε > 0. Then
since we have taken M > 0 by convention, only the following two of the four solutions in (4.8)
are physical:
M =
∓1 +
√
1 + 4εm
2ε
.
In the limit ε� 1, we have
M =
{
m− εm2,
1
ε +m− εm2.
(4.9)
There are therefore two legitimate positive energy solutions for 0 < ε � 1. This is because
the equation (4.7) is a quartic equation instead of the usual quadratic which arises when only
first-order derivatives appear in the Dirac equation. The first possibility in (4.9) represents
a perturbation to the usual rest mass and is seen also in the direct perturbative approach
of (4.6). It results in the dispersion relation
E2 ∼= k2 +m2 − 2εm3.
The other solution in (4.9) represents a non-perturbative mass generation that exists even
when m→ 0.
5 Lorentz violating nonlinear equations
There are various ways of motivating the study of Lorentz violating theories. For example, at
short distances space might not be smooth and so dynamical equations might require higher-
spatial derivatives to adequately describe the situation. However if one still restricts the time
derivatives to first or second order, to avoid potential causality problems, then clearly one has
to give up on Lorentz covariance.
We will consider nonlinear terms F which simultaneously violate Lorentz invariance [18].
The Lorentz violation will be implemented via constant background fields: in the terminology
of [37, 38] our equations will preserve the observer Lorentz covariance but break the particle
Lorentz symmetry which involves boosting the particles and local fields but not background
fields [37, 38].
In this part of the paper we illustrate some of the possibilities rather than work out all cases
as this becomes tedious and is better left for specific applications.
As Lorentz violation is constrained by phenomenology to be small [39, 40], we may use
perturbative methods to determine the corrected dispersion relations.
5.1 An example: no derivatives
If the Lorentz violation is described by background vector fields, then for F ∝ I and n = 1 we
may write
F1 = Aµ
ψ̄γµψ
ψ̄ψ
+Bµ
ψ̄γ5γ
µψ
ψ̄ψ
, (5.1)
where Aµ and Bµ are constants; current conservation requires them to be real.
14 W.K. Ng and R.R. Parwani
Under a PCT transformation of the spinors alone in (5.1) we have F → −F . Thus we have
here our first example of PCT violation associated with Lorentz violation. However it is possible
to maintain PCT while still violating Lorentz covariance. Consider
F2 = Aαβ
ψ̄σαβψ
ψ̄ψ
+ iBαβ
ψ̄γ5σ
αβψ
ψ̄ψ
,
where Aαβ and Bαβ are real background tensor fields. Both current conservation and PCT
invariance are satisfied in this case.
The dispersion relation for perturbed plane waves can be obtained using the perturbative
method, equation (4.6). For example, for the case
F = Aµ
ψ̄γµψ
ψ̄ψ
we get F = A·k
m . Thus E2 = k2 +m2 − 2A · k. Notice the correction is O(k).
6 Other cases
In this section we look at some other examples of nonlinear equations within the class (3.1) such
as those with higher nonlinearities, n ≥ 2, or with F ∝ γµ. We also discuss the Lagrangian
approach and some examples of nonlinearities outside the class (3.1).
6.1 Lorentz invariant equation with F ∝ I, n = 2
For simplicity we consider here only cases where there are no derivatives in F . An example is
given by
F = ε
(
ψ̄γ5ψ
)2(
ψ̄ψ
)2 .
It is Poincaré invariant and invariant under each of the discrete symmetries while Hermiticity
requires ε to be real. It is easy to verify, using the definition from Section 2.5 that F is separable.
6.2 Lorentz violating equation with F ∝ γµ, n = 1
Here we consider an F that is proportional to γµ. Such terms will allow the chiral current to be
conserved, as discussed in Section 2.3.2. If we exclude derivatives then the simplest possibility
is to let the Lorentz index of the gamma matrix contract with that of the background field Aµ,
F = iAµγ
µ ψ̄γ5ψ
ψ̄ψ
.
Hermiticity requires the background field to be real. This F individually breaks all the discrete
symmetries and is PCT odd! It is separable.
6.3 Equations from a Lagrangian
There are both advantages and disadvantages in using a Lagrangian approach. Firstly, a local
equation does not necessarily have a local Lagrangian. Also, even though a Lagrangian might
be simple, the resultant equations of motion might look complicated. On the other hand,
it is probably easier to discuss conservation laws corresponding to symmetries starting from
Nonlinear Dirac Equations 15
a Lagrangian. Another possible advantage of a Lagrangian approach will appear after we look
at an example.
Consider the Lagrangian density
L =
i
2
[
ψ̄γµ (∂µψ)−
(
∂µψ̄
)
γµψ
]
−mψ̄ψ + LNL.
Suppose, for simplicity, LNL contains no derivatives. Then the equation of motion will reduce
to
iγµ∂µψ −mψ +
∂LNL
∂ψ̄
= 0,
which is a similar to (1.1) and so we label the last nonlinear term here as FE.O.M.ψ. As an
example, using,
LNL =
(
ψ̄Aψ
) (
ψ̄Bψ
)
ψ̄ψ
gives
FE.O.M.ψ =
∂LNL
∂ψ̄
=
(
ψ̄Aψ
ψ̄ψ
)
Bψ +
(
ψ̄Bψ
ψ̄ψ
)
Aψ −
(
ψ̄Aψ
) (
ψ̄Bψ
)(
ψ̄ψ
)2 ψ.
Thus we see that a n = 1 nonlinearity in the Lagrangian will introduce a mixture of n = 1, 2
terms into the equations of motion. This then might be one advantage of the Lagrangian
approach: it generates constrained complexity from simplicity.
7 Gauge inequivalence
It is possible to generate a nonlinear equation from the linear Dirac equation through a nonlinear
gauge transformation [21]. The transformed equation is equivalent to the original equation in
the sense that the probability density is an invariant. Here we show that the nonlinear terms
we have investigated in this paper cannot be obtained by performing a gauge transformation on
the linear Dirac equation, and so represent genuine and distinct nonlinear structures.
We define the following gauge transformation.
ψ → ψ′(x) = eiθ(x)ψ(x),
where θ(x) is a function of ψ̄’s and ψ’s. In general, we will treat θ(x) as a 4 × 4 matrix7. We
require that the probability to be invariant under the gauge transformation,
ψ†ψ → ψ′†ψ′ =
(
eiθψ
)†
eiθψ ≡ ψ†ψ
and so θ† = θ.
Under an infinitesimal gauge transformation of the linear Dirac equation we get
(1− iθ) (iγµ∂µ −m) (1 + iθ)ψ ' 0,
(iγµ∂µ −m)ψ + [θ, γµ] ∂µψ − γµ (∂µθ)ψ ' 0. (7.1)
We wish to identify the θ dependent terms with the nonlinearity F in our nonlinear Dirac
equation (1.1) so we set
Fψ = [θ, γµ] ∂µψ − γµ (∂µθ)ψ.
7For ease of notation, we will often suppress the x-dependence in θ and ψ.
16 W.K. Ng and R.R. Parwani
Thus
ψ̄Fψ = ψ̄ [θ, γµ] ∂µψ − ψ̄γµ (∂µθ)ψ. (7.2)
We note that equation (7.2) is not symmetric in ∂µψ and so this representation of F is not
Hermitian. In order to obtain a symmetric equation, we will repeat the above steps on the
adjoint Dirac equation (this also removes any ambiguity when taking the adjoint of ∂µ).
For the adjoint equation, we have for an infinitesimal gauge transformation,
0 = ψ̄
(
iγµ←−∂µ +m
)
+
(
∂µψ̄
) [
γ0θγ0γµ − γµθ
]
+ ψ̄γ0 (∂µθ) γ0γµ + imψ̄
(
θ − γ0θγ0
)
. (7.3)
Now the adjoint of (1.1) is, upon using the Hermiticity constraint (2.1),
ψ̄ (iγµ∂µ +m)− ψ̄F = 0.
Thus comparing with (7.3) we label
ψ̄F = −
(
∂µψ̄
) [
γ0θγ0γµ − γµθ
]
− ψ̄γ0 (∂µθ) γ0γµ − imψ̄
(
θ − γ0θγ0
)
. (7.4)
Multiplying (7.4) by ψ from the right and adding to (7.2) gives
2ψ̄Fψ = ψ̄ [θ, γµ] ∂µψ −
(
∂µψ̄
) [
γ0θγ0γµ − γµθ
]
ψ
− ψ̄
[
∂µ
(
γµθ + γ0θγ0γµ
)]
ψ − imψ̄
(
θ − γ0θγ0
)
ψ. (7.5)
The left hand-side is Hermitian if the constraint (2.1) on F is applied. But the adjoint of the
right-hand-side is(
∂µψ̄
) (
γµγ0θγ0 − γ0θγ0γµ
)
ψ − ψ̄
(
γµθ − γ0θγ0γµ
)
∂µψ
− ψ̄
[
∂µ
(
γ0θγ0γµ + γµθ
)]
ψ + imψ̄
(
γ0θγ0 − θ
)
ψ. (7.6)
Comparing (7.5) and (7.6), we require
γ0θγ0 = θ ≡
[
θ, γ0
]
= 0.
Then (7.5) becomes
2ψ̄Fψ = ψ̄ [θ, γµ] ∂µψ −
(
∂µψ̄
)
[θ, γµ]ψ − ψ̄ [∂µ {θ, γµ}]ψ. (7.7)
So far we have deduced two constraints on θ,
θ = θ†,
[
θ, γ0
]
= 0,
coming respectively from the invariance of the probability density and Hermiticity. These are
necessary constraints for a nonlinear equation generated by gauge transformation to be equiva-
lent to a theory of our general class, but one must still check if any candidate solution, θ, is
actually a solution, that is, sufficiency is not guaranteed by (7.7).
7.1 Lorentz invariant case
We will now look at the constraint from Poincaré invariance. Recall that we need S−1F ′S ≡ F
under ψ → ψ′ = Sψ. The l.h.s. of (7.7) is clearly invariant while the r.h.s. transforms into
ψ̄S−1
[
θ′, γµ
]
Λν
µS∂νψ −
(
∂νψ̄
)
S−1Λν
µ
[
θ′, γµ
]
Sψ − ψ̄
[
∂νS
−1Λν
µ
{
θ′, γµ
}
S
]
ψ. (7.8)
Nonlinear Dirac Equations 17
Comparing (7.8) with (7.7), we get[
S−1θ′S, γν
]
≡ [θ, γν ] ,
{
S−1θ′S, γν
}
≡ {θ, γν} .
Thus we have the constraint
S−1θ′S = θ,
which for an infinitesimal Lorentz transformation gives
θ′ − i
4
ωab
[
θ′, σab
]
= θ.
Therefore in total we have 3 constraints,
constraint 1: θ† = θ,
constraint 2:
[
θ, γ0
]
= 0,
constraint 3: θ′ − i
4
ωab
[
θ′, σab
]
= θ.
From constraint 2, θ must be proportional to I or γ0. If θ ∝ I, then all constraints are satisfied
but for θ ∝ γ0, we cannot satisfy constraint 3: Let θ = gγ0, where g is a scalar function of the
wavefunctions. Then the Poincaré transformed θ is given by θ′ = g′γ0. Substituting this into
the left-hand-side of constraint 3, we get
g′γ0 − i
4
g′ωab
[
γ0, σab
]
= g′γ0 − i
4
g′ωab (γa − γb) .
Since ωab (γa − γb) is non-zero, the result is not proportional to γ0 and so θ ∝ γ0 does not satisfy
constraint 3. Thus we conclude that θ can only be proportional to I.
Hence with θ ∝ I, equation (7.7) becomes
ψ̄Fψ = −1
2
ψ̄ [∂µ {θ, γµ}]ψ = −ψ̄ (∂µθ) γµψ = −jµ∂µθ. (7.9)
Consider the specific case where F is proportional to I. Writing F = fI, we deduce from (7.9)
that
f = −(∂µθ) jµ
ψ̄ψ
. (7.10)
Remember that θ is a function of ψ̄’s and ψ’s, and recall our condition (2.4): we see therefore
that θ must be invariant under a scaling of the wavefunction. As long as the nonlinearities
cannot be expressed in the form shown in (7.10), we can be sure that they cannot be obtained
by performing a gauge transformation on the linear Dirac equation. In particular we conclude
that the Lorentz covariant nonlinear Dirac equations we have explicitly studied in this paper
are not gauge equivalent to the linear Dirac equation.
Now consider the class of nonlinearities where F is proportional to γµ. We let F = fµγ
µ,
where fµ are functions of ψ̄’s and ψ’s. Then (7.9) becomes
fµψ̄γ
µψ = fµj
µ = −jµ∂µθ. (7.11)
Therefore if fµ cannot be expressed as a total derivative of a scale-invariant θ function like (7.11)
then those nonlinear structures proportional to γµ cannot be obtained from the linear Dirac
equation by a gauge transformation. In particular the cases we considered in Section 6.2 are
safe.
18 W.K. Ng and R.R. Parwani
7.2 Lorentz violating cases
Finally let us consider the case where F is Lorentz violating. We have constructed our Lorentz
violating terms by introducing a constant background field Aµ (independent of the wavefunc-
tion). We may write F as AµG
µ where Gµ is the nonlinear factor which may be proportional
to I, γµ etc.
Could the Lorentz violating examples we have considered be obtained by a nonlinear gauge
transformation of the linear Dirac equation with or without Lorentz violation? The linear Dirac
equation to start with would now be of the form
(iγµ∂µ −m)ψ + LV ψ = 0,
where LV is a state-independent Lorentz violating term, if it is not zero (we assume that LV
does not have free derivatives that act to the right on ψ). Gauge transforming this equation with
a state-dependent but Hermitian θ ∝ I can generate at most Lorentz covariant nonlinearities. So
consider the other possibility, θ ∝ γ0. Then one would generate Lorentz violating nonlinearties
and on the right-hand-side of (7.1) there would be an additional term ∼ [LV, γ0]. Now if we
write θ = θ̄γ0, (7.7) becomes
ψ̄Fψ = θ̄ψ†γi∂iψ − θ̄
(
∂iψ
†)γiψ −
(
∂
∂t
θ̄
)
ψ̄ψ. (7.12)
The first observation is that in order to write the right-hand-side in covariant form we need to
introduce background tensor (for the first two terms) and vector (for the last term) fields. Also
from the structural form of our F (3.1), we see by comparing both sides of (7.12) that θ must
be invariant under scaling of the wavefunction. The examples we have explicitly discussed in
this paper therefore do not fall under the category of nonlinearities described by (7.12). For
example, with F = fI, (7.12) becomes
f = θ̄
[
ψ†γi∂iψ −
(
∂iψ
†) γiψ
ψ̄ψ
]
− ˙̄θ, (7.13)
which means having at least n = 2 and a simultaneous use of tensor and vector fields: these
are necessary conditions for the nonlinearity to be obtained through a Lorentz violating gauge
transformation of the usual linear Dirac equation.
8 Discussion
In [18] it was suggested that fundamental quantum nonlinearities might be related to potential
Lorentz violation [37, 38, 39, 40]. This current paper is a step towards a quantitative study of the
suggestions in [18]. We have discussed a framework for systematically constructing nonlinear
Dirac equations, at the quantum mechanical level, that satisfy other conventional properties
such as Hermiticity, Poincaré invariance and ψ → λψ invariance although, as shown, even those
can be relaxed.
We gave several examples of such equations, different in structure from those studied pre-
viously in the literature, and discussed their properties. We also demonstrated that our equations
were not gauge equivalent to the linear Dirac equation. More explicit examples of our class of
nonlinear Dirac equations may be found in [41] and their non-relativistic limit is studied in [42].
As mentioned in Section 1, one application of such equations is to study neutrino oscilla-
tions [43] which would be an ideal probe of quantum nonlinearities, with of without a simultane-
ous Lorentz violation [18]. Other examples we hope to study with the nonlinear equations are CP
violation and dark matter/energy. In this regard, it would be useful to obtain non-plane-wave
Nonlinear Dirac Equations 19
solutions to our nonlinear equations, similar to what has been done for simpler polynomial-type
nonlinear Dirac equations in [22, 23].
A number of authors had argued that nonlinear quantum evolution of states within the stan-
dard kinematical framework of quantum theory would lead to pathologies. However, on closer
examination, such attempts at “no go” theorems were seen to require one or more assump-
tions that are not very obvious on physical grounds; for detailed critiques and citations to the
literature the interested reader is referred to [44, 45].
We have kept open the possibility that the nonlinearities we proposed might be fundamental,
effective or only phenomenological. Of course there is less contention if the nonlinearities are
only an approximate representation of more complex underlying dynamics; in any case, from
a Wilsonian perspective, one deals in physics with a sequence of approximate theories.
Effective or phenomenological nonlinear equations are quite common in the non-relativistic
domain [6, 7] and there are also a few examples of phenomenological relativistic nonlinear equa-
tions [9, 10, 11]. As another possibility of the latter case, we note that some condensed matter
systems have (linear) relativistic-looking equations for their quasi-particles [46]: these are surely
approximations to nonlinear equations.
References
[1] Schrödinger E., Collected papers on wave mechanics, Blackie & Son, London, 1928.
[2] Bolinger J.J., Heinzen D.J., Itano W.M., Gilbert S.L., Wineland D.J., Test of the linearity of quantum
mechanics by rf spectroscopy of the 9Be+ ground state, Phys. Rev. Lett. 63 (1989), 1031–1034.
[3] Chupp T., Hoare R., Coherence in freely precessing 21Ne and a test of linearity of quantum mechanics, Phys.
Rev. Lett. 64 (1990), 2261–2264.
[4] Walsworth R.L., Silvera I.F., Mattison E.M., Vessot R.F.C., Test of the linearity of quantum mechanics in
an atomic system with a hydrogen maser, Phys. Rev. Lett. 64 (1990), 2599–2602.
[5] Majumder P.K., Venema B.J., Lamoreaux S.K., Heckel B.R., Test of the linearity of quantum mechanics in
optically pumped 201Hg, Phys. Rev. Lett. 65 (1990), 2931–2934.
[6] Akhmediev N.N., Ankiewicz A., Solitons: nonlinear pulses and beams, Chapman and Hall, 1997.
[7] Pethick C., Smith H., Bose–Einstein condensation in dilute gases, Cambridge University Press, 2001.
[8] Heisenberg W., Quantum theory of fields and elementary particles, Rev. Mod. Phys. 29 (1957), 269–278.
[9] Ionescu D.C., Reinhardt R., Muller B., Greiner W., Nonlinear extensions of the Dirac equation and their
implications in QED, Phys. Rev. A 38 (1988), 616–620.
[10] Zecca A., Dirac equation in space-time with torsion, Internat. J. Theoret. Phys. 41 (2002), 421–428.
[11] Esteban M.J., Sere E., An overview on linear and nonlinear Dirac equations, Discrete Contin. Dyn. Syst. 8
(2002), 381–397.
[12] Esteban M.J., Sere E., Stationary states of the nonlinear Dirac equation: a variational approach, Comm.
Math. Phys. 171 (1995), 323–350.
[13] Bialynicki-Birula I., Mycielski J., Nonlinear wave mechanics, Ann. Physics 100 (1976), 62–93.
[14] Kibble T.W.B., Relativistic models of nonlinear quantum mechanics, Comm. Math. Phys. 64 (1978), 73–82.
[15] Weinberg S., Testing quantum mechanics, Ann. Physics 194 (1989), 336–386.
[16] Doebner H.D., Goldin G.A., Nattermann P., Gauge Transformations in quantum mechanics and the
unification of nonlinear Schrödinger equations, J. Math. Phys. 40 (1999), 49–63, quant-ph/9709036.
[17] Parwani R., Why is Schrödinger’s equation linear?, Braz. J. Phys. 35 (2005), 494–496, quant-ph/0412192.
[18] Parwani R., An information theoretic link between space-time symmetries and quantum linearity, Ann.
Physics 315 (2005), 419–452, hep-th/0401190.
[19] Parwani R., Information measures for inferring quantum mechanics and its deformations, J. Phys. A: Math.
Gen. 38 (2005), 6231–6237, quant-ph/0408185.
[20] Parwani R., A physical axiomatic approach to Schrödinger’s equation, Internat. J. Theoret. Phys. 45 (2006),
1901–1933, quant-ph/0508125.
http://arxiv.org/abs/quant-ph/9709036
http://arxiv.org/abs/quant-ph/0412192
http://arxiv.org/abs/hep-th/0401190
http://arxiv.org/abs/quant-ph/0408185
http://arxiv.org/abs/quant-ph/0508125
20 W.K. Ng and R.R. Parwani
[21] Doebner H.D., Zhdanov R., Nonlinear Dirac equations and nonlinear gauge transformations,
quant-ph/0304167.
[22] Fushchich W.I., Zhdanov R.Z., Symmetry and exact solutions of nonlinear spinor equations, Phys. Rep. 172
(1989), 123–174.
[23] Fushchych W.I., Zhdanov R.Z., Symmetries and exact solutions of nonlinear Dirac equations, Mathematical
Ukraina Publisher, Kyiv, 1997, math-ph/0609052.
[24] Bogoslovsky G.Yu., Lorentz symmetry violation without violation of relativistic symmetry, Phys. Lett. A
350 (2006), 5–10, hep-th/0511151.
[25] Bogoslovsky G.Yu., Goenner H.F., Concerning the generalized Lorentz symmetry and the generalization of
the Dirac equation, Phys. Lett. A 323 (2004), 40–47, hep-th/0402172.
[26] Bogoslovsky G.Yu., A special-relativistic theory of the locally anisotropic space-time, Nuovo Cimento B 40
(1977), 99–134.
[27] Itzykson C., Zuber J.B., Quantum field theory, McGraw-Hill International Book Co., New York, 1980.
[28] Kostelecky V.A., Gravity, Lorentz violation, and the standard model, Phys. Rev. D 69 (2004), 105009,
20 pages, hep-th/0312310.
[29] Bohm D., Hiley B.J., The undivided universe. An ontological interpretation of quantum theory, Routledge,
London, 2003.
[30] Moshinsky M., Nikitin A., The Many body problem in relativistic quantum mechanics, Rev. Mexicana F́ıs.
50 (2005), 66–73, hep-ph/0502028.
[31] Goldstein S., Tumulka R., Opposite arrows of time can reconcile relativity and nonlocality, Classical Quan-
tum Gravity 20 (2003), 557–564, quant-ph/0105040.
[32] Svetlichny G., On linearity of separating multiparticle differential Schrödinger operators for identical partic-
les, J. Math. Phys. 45 (2004), 959–964, quant-ph/0308001.
[33] Goldin G.A., Svetlichny G., Nonlinear Schrödinger equations and the separation property, J. Math. Phys.
35 (1994), 3322–3332.
[34] Parwani R., Tan H.S., Exact solutions of a non-polynomially nonlinear Schrödinger equation, Phys. Lett. A
363 (2007), 197–201, quant-ph/0605123.
[35] Parwani R., Tabia G., Universality in an information-theoretic motivated nonlinear Schrödinger equation,
J. Phys. A: Math. Gen. 40 (2007), 5621–5635, quant-ph/0607222.
[36] Greenberg O.W., CPT violation implies violation of Lorentz invariance, Phys. Rev. Lett. 89 (2002), 231602,
4 pages, hep-ph/0201258.
[37] Colladay D., Kostelecky V.A., CPT violation and the standard model, Phys. Rev. D 55 (1997), 6760–6774,
hep-ph/9703464.
[38] Lehnert R., Dirac theory within the standard-model extension, J. Math. Phys. 45 (2004), 3399–3412,
hep-ph/0401084.
[39] Jacobson T., Liberate S., Mattingly D., Lorentz violation at high energy: concepts, phenomena, and astro-
physical constraints, Ann. Physics 321 (2006), 150–196, astro-ph/0505267.
[40] Vucetich H., Testing Lorentz invariance violation in quantum gravity theories, gr-qc/0502093.
[41] Ng W.K., Parwani R., Nonlinear Dirac equations, arXiv:0707.1553v1.
[42] Ng W.K., Parwani R., Nonlinear Schrödinger–Pauli equations, arXiv:0807.1877.
[43] Ng W.K., Parwani R., Probing quantum nonlinearities through neutrino oscillations, arXiv:0805.3015.
[44] Kent A., Nonlinearity without superluminality, Phys. Rev. A 72 (2005), 012108, 4 pages, quant-ph/0204106.
[45] Holman M., On arguments for linear quantum dynamics, quant-ph/0612209.
[46] Geim A.K., Novoselov K.S., The rise of graphene, Nature Materials 6 (2007), 183–191, cond-mat/0702595.
http://arxiv.org/abs/quant-ph/0304167
http://arxiv.org/abs/math-ph/0609052
http://arxiv.org/abs/hep-th/0511151
http://arxiv.org/abs/hep-th/0402172
http://arxiv.org/abs/hep-th/0312310
http://arxiv.org/abs/hep-ph/0502028
http://arxiv.org/abs/quant-ph/0105040
http://arxiv.org/abs/quant-ph/0308001
http://arxiv.org/abs/quant-ph/0605123
http://arxiv.org/abs/quant-ph/0607222
http://arxiv.org/abs/hep-ph/0201258
http://arxiv.org/abs/hep-ph/9703464
http://arxiv.org/abs/hep-ph/0401084
http://arxiv.org/abs/astro-ph/0505267
http://arxiv.org/abs/gr-qc/0502093
http://arxiv.org/abs/0707.1553v1
http://arxiv.org/abs/0807.1877
http://arxiv.org/abs/0805.3015
http://arxiv.org/abs/quant-ph/0204106
http://arxiv.org/abs/quant-ph/0612209
http://arxiv.org/abs/cond-mat/0702595
1 Introduction
2 Constraints
2.1 Locality
2.2 Poincaré invariance
2.3 Hermiticity
2.3.1 Current conservation
2.3.2 Chiral current
2.4 Universality
2.5 Separability
2.5.1 Structure of F
2.6 Discrete symmetries
3 Explicit examples of nonlinear equations with F \propto I, n = 1
3.1 No derivatives
3.1.1 Lorentz vs PCT invariance
3.2 One derivative
3.3 Two derivatives
4 Plane-wave solutions and dispersion relations
4.1 Perturbative method
4.2 Example
5 Lorentz violating nonlinear equations
5.1 An example: no derivatives
6 Other cases
6.1 Lorentz invariant equation with F \propto I, n = 2
6.2 Lorentz violating equation with F \propto \gamma_\mu, n = 1
6.3 Equations from a Lagrangian
7 Gauge inequivalence
7.1 Lorentz invariant case
7.2 Lorentz violating cases
8 Discussion
References
|
| id | nasplib_isofts_kiev_ua-123456789-149176 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-0659 |
| language | English |
| last_indexed | 2025-12-07T16:52:12Z |
| publishDate | 2009 |
| publisher | Інститут математики НАН України |
| record_format | dspace |
| spelling | Wei Khim Ng Rajesh R. Parwani 2019-02-19T18:13:18Z 2019-02-19T18:13:18Z 2009 Nonlinear Dirac Equations / Wei Khim Ng, Rajesh R. Parwani // Symmetry, Integrability and Geometry: Methods and Applications. — 2009. — Т. 5. — Бібліогр.: 46 назв. — англ. 1815-0659 2000 Mathematics Subject Classification: 81P05; 81Q99; 83A05 https://nasplib.isofts.kiev.ua/handle/123456789/149176 We construct nonlinear extensions of Dirac's relativistic electron equation that preserve its other desirable properties such as locality, separability, conservation of probability and Poincaré invariance. We determine the constraints that the nonlinear term must obey and classify the resultant non-polynomial nonlinearities in a double expansion in the degree of nonlinearity and number of derivatives. We give explicit examples of such nonlinear equations, studying their discrete symmetries and other properties. Motivated by some previously suggested applications we then consider nonlinear terms that simultaneously violate Lorentz covariance and again study various explicit examples. We contrast our equations and construction procedure with others in the literature and also show that our equations are not gauge equivalent to the linear Dirac equation. Finally we outline various physical applications for these equations. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Nonlinear Dirac Equations Article published earlier |
| spellingShingle | Nonlinear Dirac Equations Wei Khim Ng Rajesh R. Parwani |
| title | Nonlinear Dirac Equations |
| title_full | Nonlinear Dirac Equations |
| title_fullStr | Nonlinear Dirac Equations |
| title_full_unstemmed | Nonlinear Dirac Equations |
| title_short | Nonlinear Dirac Equations |
| title_sort | nonlinear dirac equations |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/149176 |
| work_keys_str_mv | AT weikhimng nonlineardiracequations AT rajeshrparwani nonlineardiracequations |