Home Which magnetic fields support a zero mode?
Article
Licensed
Unlicensed Requires Authentication

Which magnetic fields support a zero mode?

  • Rupert L. Frank ORCID logo and Michael Loss ORCID logo EMAIL logo
Published/Copyright: April 28, 2022

Abstract

This paper presents some results concerning the size of magnetic fields that support zero modes for the three-dimensional Dirac equation and related problems for spinor equations. It is a well-known fact that for the Schrödinger equation in three dimensions to have a negative energy bound state, the 3/2 norm of the potential has to be greater than the Sobolev constant. We prove an analogous result for the existence of zero modes, namely that the 3/2 norm of the magnetic field has to greater than twice the Sobolev constant. The novel point here is that the spinorial nature of the wave function is crucial. It leads to an improved diamagnetic inequality from which the bound is derived. While the results are probably not sharp, other equations are analyzed where the results are indeed optimal.

Funding statement: Partial support through U.S. National Science Foundation grants DMS-1363432 and DMS-1954995 (Rupert L. Frank) and DMS-1856645 (Michael Loss) and through the Deutsche Forschungsgemeinschaft (German Research Foundation) through Germany’s Excellence Strategy EXC-2111-390814868 (Rupert L. Frank) is acknowledged.

A Some computations involving the Dirac matrices

The construction of zero modes in higher dimensions is more complicated and, as mentioned before, was accomplished by Dunne and Min [11] using information about the Dirac equation on the sphere. The advantage of their construction is that it delivers automatically the dimension of the zero mode space. If one is satisfied with less information, then there is, we believe, a simpler way to construct the Dunne–Min zero modes. Moreover, it gives the opportunity to get acquainted with some of the properties of the Dirac matrices. The basic idea is due to Adolf Hurwitz in his posthumously published paper “Über die Komposition der quadratischen Formen” [25]. In this paper he gave a complete classification of matrices γj, j=1,,d, satisfying the relations

γjγk+γkγj=2δij.

For our purpose we shall assume the matrices γj to be self-adjoint in the space N with the usual inner product.

Theorem A.1.

Let d=2ν+1 or d=2ν, and consider the N×N hermitian matrices γj, j=1,,d, satisfying

γiγj+γjγi=2δij.

Then N=2ν and, if γj is another set of 2ν×2ν hermitian matrices satisfying the same relations, then there exists a 2ν×2ν unitary matrix A such that γj=A*γjA for j=1,,d.

The proof proceeds by reducing the γ matrices to a unitarily equivalent, but canonical set of matrices using an inductive procedure. We omit the proof and refer to [25].

Corollary A.2.

Let R be a d×d orthogonal matrix and define

γj=k=1dRjkγk.

Then there exists a unitary matrix A such that for all j=1,,d one has γj=A*γjA.

The computation with γ matrices can be sometimes tedious and the following framework called “second quantization” is quite helpful.

In the remainder of this section, we assume that d=2ν+1 is odd.

We single out the matrix γ1 and define the “annihilation” and “creation” operators

cj:=12(γ2j+iγ2j+1),cj*=12(γ2j-iγ2j+1),j=1,2,,ν,

so that

γ2j=cj+cj*,γ2j+1=1i(cj-cj*).

One easily checks that

cjcj*+cj*cj=I,cj2=cj*2=0

and, for k,

ckc+cck=0,ck*c*+c*ck*=0,ckc*+c*ck=0.

Note that the matrix γ1 is not involved in these definitions.

Lemma A.3.

There exists a vector ϕC2ν, a vacuum, such that ϕ=1 and

cjϕ=0,j=1,,ν.

Proof.

Since c12=0, it is clear that there exists ϕ0 such that c1ϕ=0. Let k be the first index such that ckϕ0. Setting ψ=ckϕ, we see because of the commutation relations that cjψ=0, j=1,,k-1, and ckψ=ck2ϕ=0. Thus, replacing ϕ by ψ, we have that cjψ=0, j=1,,k. Continuing in this fashion, we have a vector ϕ such that ckϕ=0 for all k=1,,ν. ∎

Lemma A.4.

Let β=(β1,,βν) be a sequence with βj{0,1}, j=1,,ν. Then the vector

|β=ck1*β1ckν*βνϕ

is non-zero if and only if the indices k1,,kν are all distinct. In this case the vector is normalized. Moreover, the vectors |β form an orthonormal basis in C2ν.

In view of this lemma, we will sometimes denote ϕ=|0.

Proof.

If one or more of the indices are not distinct, then by commuting the various operators results in a square of one of the ci*, which is zero. Hence we may assume that the indices k1,,kν are all distinct. We also may assume that β1=1 because otherwise ck1β1=I and we may move on to the next index. We have

ck1*β1ckν*βνϕ2=(|0,ckνβνck1β1ck1*β1ckν*βνϕ)

and using ck1β1ck1*β1=I-ck1*β1ck1β1, we find

(|0,ckνβνck1β1ck1*β1ckν*βνϕ)=(|0,ckνβνck2β2ck2*β2ckν*βνϕ)
-(|0,ckνβνck1*β1ck1β1ckν*βνϕ).

The second term on the right side vanishes because the indices are distinct and thus ck1β1 either commutes or anti-commutes with all the matrices on the right and once it hits ϕ it yields zero. In this fashion we may move the annihilation matrices to the right and obtain that this state is normalized. Incidentally this also makes it clear that the state vanishes if two indices are the same on account of the fact that cj2=cj*2=0. From this argument it also follows that for ββ,

(|β,|β)=0

and hence we have 2ν orthonormal vectors which constitute an orthonormal basis. ∎

Lemma A.5.

The vacuum is unique (up to a constant phase).

Proof.

Suppose that v is another vacuum, i.e., v=1 and for all α=1,,ν,

cαv=0.

We may assume that 0|v=0. Then

(v,ck1*β1ckν*βνϕ)

is always zero and therefore, by Lemma A.4, v=0, which is a contradiction. ∎

We note that γ1ϕ satisfies the same properties as ϕ, namely, γ1ϕ=1 and

cαγ1ϕ=-γ1cαϕ=0

for all α. By the uniqueness result of Lemma A.5 there is a θ such that γ1ϕ=eiθϕ. Since γ1 is self-adjoint, we have eiθ=±1. In case it is -1, we can change the sign of γ1 without changing the commutation relations and arrive at the same relation with +1. Hence we may adopt the convention that γ1ϕ=ϕ.

The point about introducing this formalism is the following result.

Lemma A.6.

Introduce a (2ν+1)×(2ν+1) matrix ω with entries

ωα,β={0𝑖𝑓α=1𝑜𝑟β=1𝑜𝑟α=β,0|iγαγβ|0𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒.

Then

ω=diag(0,-iσ2,,-iσ2),

where the zero is a number and there are ν 2×2-blocks -iσ2.

Proof.

Since ω is skew and vanishes on the diagonal, it suffices to compute ωα,β when α<β. Moreover, since γ1|0=|0 and γ1 is selfadjoint, we have ω1,β=0 for all β>1. For the remaining entries, we need to distinguish whether α and β are even or odd. When both are even, we have for 1j<kν,

ω2j,2k=0|i(cj+cj*)(ck+ck*)|0=i0|cjck*|0=0.

Similarly, when both are odd, we have for 1j<kν,

ω2j+1,2k+1=-0|i(cj-cj*)(ck-ck*)|0=-i0|cjck*|0=0.

Next, we consider α is even and β is odd. If α=β-1, we get

ω2j,2j+1=0|(cj+cj*)(cj-cj*)|0=-1.

Otherwise, for 1j<kν,

ω2j,2k+1=0|(cj+cj*)(ck-ck*)|0=0.

Finally, we have the case where α is odd and β is even. For 1j<kν, we get

ω2j+1,2k=0|(cj-cj*)(ck+ck*)|0=0.

This proves the claimed formula for the entries of the matrix ω. ∎

After these preliminaries we discuss now an alternative approach to the Dunne–Min generalization [11] of [30]. The following example is relevant. It is the higher-dimensional analog of choice for the vector potential in [30]. Consider

(A.1)((1+iγx)η,γ(1+iγx)η),

where η2ν is normalized. Recall that for d=3 the γ matrices are the Pauli matrices and there is the well-know identity

|(η,ση)|2=|η|4.

This leads to the identity

σ(η,ση)η=η,

which is very useful for constructing zero modes. It turns out that this identity also holds for d=5, but not in higher dimensions. In particular, it does not hold for (A.1) for general η.

Things simplify considerably if we choose the constant spinor η to be the vacuum ϕ. We compute

(1-iγx)γj(1+iγx)=γj-iγxγj+iγjγx+γxγjγx
=γj-2ixγγj+ix(γjγ+γγj)-(xγ)2γj
+xγx(γjγ+γγj)
=γj-2ixγγj+2ixj-|x|2γj+2xγxj
=(1-|x|2)γj+2xγxj-2ixγγj+2ixj.

Taking expectation we get

((1+iγx)ϕ,γj(1+iγx)ϕ)=(1-|x|2)(ϕ,γjϕ)+2x(ϕ,γϕ)xj-2kj,1xk(ϕ,iγkγjϕ).

Since γ1ϕ=ϕ, we find that (ϕ,γjϕ)=0, j1. Hence we have that for this particular state

((1+iγx)ϕ,γ1(1+iγx)ϕ)=(1-|x|2+2x12)|ϕ|2.

For the component j1 we find

((1+iγx)ϕ,γj(1+iγx)ϕ)=(2x1xj+2[ωx]j)|ϕ|2.

Here ω is the (2ν+1)×(2ν+1) skew matrix introduced above. We introduce the field

Uj(x):=11+|x|2((1+iγx)ϕ,γj(1+iγx)ϕ)={1-|x|2+2x121+|x|2if j=1,2x1xj+2[ωx]j1+|x|2if j1.

This can be written more concisely as

U(x)=11+|x|2((1-|x|2)e1+2(e1x)x+2ωx),

where [ωx]k=j=1dωkjxj. A straightforward computation shows that

|U(x)|2=1(1+|x|2)2((1-|x|2)e1+2(e1x)x+2ωx)2
=1(1+|x|2)2((1-|x|2)2+4(e1x)2|x|2+4|ωx|2+4(1-|x|2)(e1x)2).

Since |ωx|2=(ωx,ωx)=(x,ωTωx)=j=2dxj2, we get

|U(x)|2=1.

In other words, the vector

U(x)=(1+ixγ(1+|x|2)1/2ϕ,γ1+ixγ(1+|x|2)1/2ϕ)

is a unit vector. Now consider the self adjoint matrix

M:=U(x)γ,

whose square is |U(x)|2=1. Hence the eigenvalues of M are ±1. Moreover,

(1+ixγ(1+|x|2)1/2ϕ,M1+ixγ(1+|x|2)1/2ϕ)=|U(x)|2=1

and hence we have that

M1+ixγ(1+|x|2)1/2ϕ=1+ixγ(1+|x|2)1/2ϕ.

If we set

ψ:=1+ixγ(1+|x|2)d2ϕ,

then a simple computation yields

-iγψ=d1+|x|2ψ

and if we define

A(x):=d1+|x|2U(x),

then

-iγψ=γAψ

and we have constructed our zero modes.

B Generalization of the spin-field interaction term to arbitrary dimensions

Squaring the Dirac equation yields

[γ(-i-A)]2=jkγjγk(-ij-Aj)(-ik-Ak)
=jγjγj(-ij-Aj)(-ij-Aj)
+jkγjγk(-ij-Aj)(-ik-Ak).

We have

jγjγj(-ij-Aj)(-ij-Aj)=(-i-A)2

and

jkγjγk(-ij-Aj)(-ik-Ak)
=12jkγjγk(-ij-Aj)(-ik-Ak)+12jkγkγj(-ik-Ak)(-ij-Aj)
=12jkγjγk[(-ij-Aj)(-ik-Ak)-(-ik-Ak)(-ij-Aj)]
=i2jkγjγk[jAk-kAj].

For each fixed xd, the matrix Bjk:=jAk-kAj is an antisymmetric matrix and there is an orthogonal matrix R (depending on x) such that

RTBR=D,

where

D={diag(D1iσ2,,Dνiσ2,0)ifd=2ν+1is odd,diag(D1iσ2,,Dνiσ2)ifd=2νis even.

Here there are ν 2×2 blocks iσ2 and, if d is odd, an additional 1×1 “block” consisting of the number 0. For instance, in five dimensions

[0D1000-D10000000D2000-D20000000].

Since the trace of B is zero, we have

jγjγjBjj=0.

Hence

i2jkγjγkBjk=i2jkγjγkBjk
=i2αβjkγjγkRjαDαβRkβ
=i2αβ(jγjRjα)(kγkRkβ)Dαβ.

If we set

Γβ:=jγjRjβ,

then, according to Corollary A.2, there exists a unitary matrix U such that

Γα=U*γαU,α=1,,d,

and we can write

i2jkγjγkBjk=U*i2αβγαγβDαβU=U*i[γ1γ2D1+γ3γ4D2++γ2ν-1γ2νDν]U.

The matrices γ1γ2 and γ3γ4, etc., are skew symmetric, commute with each other and we can simultaneously diagonalize them by a unitary matrix V, that is,

γ2k-1γ2k=-iV*Σ2k-1,2kV,k=1,,ν,

with diagonal matrices Σ2k-12k. As (γiγj)2=-1, the eigenvalues of Σ2k-1,2k must be ±1. Thus, all things considered, we get

i2jkγjγkBjk=(VU)*[Σ12D1+Σ34D3+Σ2ν-1,2νDν](VU),

where the matrices Σi,i+1 are diagonal and have ±1 in the diagonal. Thus, if ψ is a spinor, we have that

|ψ,i2jkγjγkBjkψ||ψ|2k=1ν|Dk|.

This fits with the three-dimensional case where ν=1 and |D1|=|B|.

Moreover, we have

k=1ν|Dk|ν12(k=1ν|Dk|2)12=ν1/2(j<k|Bjk|2)12.

The last identity comes from the fact that conjugation by an orthogonal matrix R does not change the Hilbert–Schmidt norm of the matrix D=RTBR.

Acknowledgements

The authors would like to thank H. Kovarik and M. Lewin for helpful remarks.

References

[1] C. Adam, B. Muratori and C. Nash, Zero modes of the Dirac operator in three dimensions, Phys. Rev. D (3) 60 (1999), no. 12, Article ID 125001. 10.1103/PhysRevD.60.125001Search in Google Scholar

[2] Y. Aharonov and A. Casher, Ground state of a spin-12charged particle in a two-dimensional magnetic field, Phys. Rev. A (3) 19 (1979), no. 6, 2461–2462. 10.1103/PhysRevA.19.2461Search in Google Scholar

[3] S. Avramska-Lukarska, D. Hundertmark and H. Kovařík, Absence of positive eigenvalues of magnetic Schrödinger operators, preprint (2020), https://arxiv.org/abs/2003.07294. Search in Google Scholar

[4] A. A. Balinsky and W. D. Evans, On the zero modes of Pauli operators, J. Funct. Anal. 179 (2001), no. 1, 120–135. 10.1006/jfan.2000.3670Search in Google Scholar

[5] A. A. Balinsky, W. D. Evans and R. T. Lewis, Sobolev, Hardy and CLR inequalities associated with Pauli operators in 3, J. Phys. A 34 (2001), no. 5, L19–L23. 10.1088/0305-4470/34/5/101Search in Google Scholar

[6] R. D. Benguria and H. Van Den Bosch, A criterion for the existence of zero modes for the Pauli operator with fastly decaying fields, J. Math. Phys. 56 (2015), no. 5, Article ID 052104. 10.1063/1.4920924Search in Google Scholar

[7] W. Borrelli and R. L. Frank, Sharp decay estimates for critical Dirac equations, Trans. Amer. Math. Soc. 373 (2020), no. 3, 2045–2070. 10.1090/tran/7958Search in Google Scholar

[8] W. Borrelli, A. Malchiodi and R. Wu, Ground state Dirac bubbles and Killing spinors, Comm. Math. Phys. 383 (2021), no. 2, 1151–1180. 10.1007/s00220-021-04013-1Search in Google Scholar

[9] D. M. J. Calderbank, P. Gauduchon and M. Herzlich, Refined Kato inequalities and conformal weights in Riemannian geometry, J. Funct. Anal. 173 (2000), no. 1, 214–255. 10.1006/jfan.2000.3563Search in Google Scholar

[10] L. Cossetti, L. Fanelli and D. Krejčiřík, Absence of eigenvalues of Dirac and Pauli Hamiltonians via the method of multipliers, Comm. Math. Phys. 379 (2020), no. 2, 633–691. 10.1007/s00220-020-03853-7Search in Google Scholar

[11] G. V. Dunne and H. Min, Abelian zero modes in odd dimensions, Phys. Rev. D 78 (2008), no. 6, Article ID 067701. 10.1103/PhysRevD.78.067701Search in Google Scholar

[12] D. M. Elton, The local structure of zero mode producing magnetic potentials, Comm. Math. Phys. 229 (2002), no. 1, 121–139.10.1007/s00220-002-0679-2Search in Google Scholar

[13] L. Erdős and J. P. Solovej, Semiclassical eigenvalue estimates for the Pauli operator with strong non-homogeneous magnetic fields. II. Leading order asymptotic estimates, Comm. Math. Phys. 188 (1997), no. 3, 599–656. 10.1007/s002200050181Search in Google Scholar

[14] L. Erdős and J. P. Solovej, Semiclassical eigenvalue estimates for the Pauli operator with strong nonhomogeneous magnetic fields. I. Nonasymptotic Lieb–Thirring-type estimate, Duke Math. J. 96 (1999), no. 1, 127–173. 10.1215/S0012-7094-99-09604-7Search in Google Scholar

[15] L. Erdős and J. P. Solovej, The kernel of Dirac operators on 𝕊3 and 3, Rev. Math. Phys. 13 (2001), no. 10, 1247–1280. 10.1142/S0129055X01000983Search in Google Scholar

[16] M. J. Esteban, M. Lewin and E. Séré, Dirac-Coulomb operators with general charge distribution II. The lowest eigenvalue, Proc. Lond. Math. Soc. (3) 123 (2021), no. 4, 345–383. 10.1112/plms.12396Search in Google Scholar

[17] P. M. N. Feehan, A Kato–Yau inequality and decay estimate for eigenspinors, J. Geom. Anal. 11 (2001), no. 3, 469–489. 10.1007/BF02922015Search in Google Scholar

[18] R. L. Frank and M. Loss, A sharp criterion for zero modes of the Dirac equation, preprint (2022), https://arxiv.org/abs/2201.03610. Search in Google Scholar

[19] R. L. Frank and M. Loss, Existence of optimizers in a Sobolev inequality for vector fields, Ars Inveniendi Anal. (2022), Paper No. 1. Search in Google Scholar

[20] T. Friedrich, Dirac operators in Riemannian geometry, Grad. Stud. Math. 25, American Mathematical Society, Providence, RI 2000. 10.1090/gsm/025Search in Google Scholar

[21] J. Fröhlich, E. H. Lieb and M. Loss, Stability of Coulomb systems with magnetic fields. I. The one-electron atom, Comm. Math. Phys. 104 (1986), no. 2, 251–270. 10.1007/BF01211593Search in Google Scholar

[22] M. P. Fry, Paramagnetism, zero modes, and mass singularities in qed in 1+1, 2+1, and 3+1 dimensions, Phys. Rev. D 55 (1997), 968–972. 10.1103/PhysRevD.55.968Search in Google Scholar

[23] O. Hijazi, A conformal lower bound for the smallest eigenvalue of the Dirac operator and Killing spinors, Comm. Math. Phys. 104 (1986), no. 1, 151–162. 10.1007/BF01210797Search in Google Scholar

[24] O. Hijazi, Première valeur propre de l’opérateur de Dirac et nombre de Yamabe, C. R. Acad. Sci. Paris Sér. I Math. 313 (1991), no. 12, 865–868. Search in Google Scholar

[25] A. Hurwitz, Über die Komposition der quadratischen Formen, Math. Ann. 88 (1922), no. 1–2, 1–25. 10.1007/978-3-0348-4160-3_46Search in Google Scholar

[26] E. Jannelli and S. Solimini, Concentration estimates for critical problems, Ric. Mat. 48 (1999), 233–257. Search in Google Scholar

[27] H. Kalf, T. Okaji and O. Yamada, The Dirac operator with mass m00: Non-existence of zero modes and of threshold eigenvalues, Doc. Math. 20 (2015), 37–64. 10.4171/dm/485Search in Google Scholar

[28] H. Kalf, T. Okaji and O. Yamada, Erratum for “The Dirac operator with mass m00: Non-existence of zero modes and of threshold eigenvalues”, Doc. Math. 24 (2019), 1361–1363. 10.4171/dm/706Search in Google Scholar

[29] E. H. Lieb and M. Loss, Analysis, 2nd ed., Grad. Stud. Math. 14, American Mathematical Society, Providence, RI 2001. Search in Google Scholar

[30] M. Loss and H.-T. Yau, Stability of Coulomb systems with magnetic fields. III. Zero energy bound states of the Pauli operator, Comm. Math. Phys. 104 (1986), no. 2, 283–290. 10.1007/BF01211595Search in Google Scholar

[31] C. Ross and B. J. Schroers, Magnetic zero-modes, vortices and Cartan geometry, Lett. Math. Phys. 108 (2018), no. 4, 949–983. 10.1007/s11005-017-1023-2Search in Google Scholar PubMed PubMed Central

[32] Y. Saitō and T. Umeda, The zero modes and zero resonances of massless Dirac operators, Hokkaido Math. J. 37 (2008), no. 2, 363–388. 10.14492/hokmj/1253539560Search in Google Scholar

[33] Y. Saitō and T. Umeda, A sequence of zero modes of Weyl–Dirac operators and an associated sequence of solvable polynomials, Spectral theory, function spaces and inequalities, Oper. Theory Adv. Appl. 219, Birkhäuser/Springer, Basel (2012), 197–209. 10.1007/978-3-0348-0263-5_11Search in Google Scholar

[34] B. v. Sz. Nagy, Über Integralungleichungen zwischen einer Funktion und ihrer Ableitung, Acta Univ. Szeged. Sect. Sci. Math. 10 (1941), 64–74. Search in Google Scholar

[35] H. Weyl, The method of orthogonal projection in potential theory, Duke Math. J. 7 (1940), 411–444. 10.1215/S0012-7094-40-00725-6Search in Google Scholar

Received: 2021-01-29
Revised: 2021-11-29
Published Online: 2022-04-28
Published in Print: 2022-07-01

© 2022 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 22.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/crelle-2022-0015/html?lang=en
Scroll to top button