arxiv:math/ v1 [math.ap] 2 Jul 2006

Transcription

1 arxiv:math/677v1 [math.ap] Jul 6 Relativistic hydrogenic atoms in strong magnetic fields Jean Dolbeault 1, Maria J. Esteban 1, Michael Loss 1 Ceremade (UMR CNRS no. 7534), Université Paris Dauphine, Place de Lattre de Tassigny, Paris Cédex 16, France School of Mathematics, Georgia Institute of Technology, Atlanta, GA 333, USA October 1, 17 Abstract: In the Dirac operator framework we characterize and estimate the ground state energy of relativistic hydrogenic atoms in a constant magnetic field and describe the asymptotic regime corresponding to a large field strength using relativistic Landau levels. We also define and estimate a critical magnetic field beyond which stability is lost. Key words. Dirac-Coulomb Hamiltonian relativistic hydrogen atom constant magnetic field Landau levels min-max levels Hardy inequality selfadjoint operators AMS classification (). 35Q4, 35Q75, 46N5, 81Q1; 34L4, 35P5, 47A5, 47N5, 81V45 1. Introduction In this paper we characterize the ground state energy of hydrogenic atoms in magnetic fields. We deal with fields of large strength in the Dirac operator framework, far away from the perturbative regime. To compute eigenvalues of Dirac operators, the usual min-max principle does not apply. More sophisticated versions of this principle have been established over the last few years, see [11,1,5]. These techniques are powerful enough to provide accurate and efficient algorithms for calculating eigenvalues of Dirac operators [7, 6]. In this paper we demonstrate that they are also flexible enough to coverthe casewith a magneticfield and providereasonableresultsfor ahighly non-perturbative problem, when paired with the right physical insight. The Dirac operator for a hydrogenic atom in the presence of a constant magnetic field B in the x 3 -direction is given by H B ν [ 1 with H B := α i + 1 ] B( x,x 1,) +β, (1) where ν = Zα < 1, Z is the nuclear charge number. The Sommerfeld finestructure constant is α 1/ The energy is measured in units of mc, i.e., the restenergyofthe electron,the length in units of /mc, i.e., the Compton wavelengthdividedbyπ,andthemagneticfieldstrengthb ismeasuredinunits of m c q Tesla. Here m is the mass of the electron, c the speed of light, q the charge of the electron (measured in Coulomb) and is Planck s constant divided by π. It is worth recalling the the earth s magnetic field is of the order of 1 Gauss and 1 Tesla is 1 4 Gauss.

2 Jean Dolbeault, Maria J. Esteban, Michael Loss In (1), α 1, α, α 3 and β are 4 4 complex matrices, whose standard form (in blocks) is ( ) ( ) I σk β =, α I k = (k = 1,,3), σ k ( ) 1 where I = and σ 1 k are the Pauli matrices: σ 1 = ( ) 1, σ 1 = ( ) ( ) i 1, σ i 3 =. 1 The magnetic Dirac operator without the Coulomb potential has essential spectrum (, 1] [1, ) and no eigenvalues in the gap ( 1,1). For ν (,1) the Hamiltonian (1) has the same essential spectrum and eigenvalues in the gap. The ground state energy λ 1 (ν,b) is defined as the smallest eigenvalue in the gap. As the field gets largeenough, one expects that the groundstate energyof the Dirac operator decreases and eventually penetrates the lower continuum. The implication of this for a second quantized model is that electron positron pair creation comes into the picture [18, ]. The intuition for that can be gleaned from the Pauli equation, where the magnetic field tends to lower the energy because of the spin. It is therefore reasonable to define the critical field strength B(ν) as the supremum of the positive B s for which λ 1 (ν,b) is in the gap ( 1,1) for all b (,B). As a function of ν, λ 1 (ν,b) is non-increasing. Hence the function B(ν) is also non-increasing. One of our goals is to give estimates on this critical field as a function of the nuclear charge. Our first result, proved in Section is that this critical field exists and we give some rough estimates in terms of ν: For some C >, 4 B(ν) min 5ν ( 18πν [3ν ] + ), e C/ν. () As a corollary we get the noteworthy result that as ν 1 the critical field B(ν) stays strictly positive. This is somewhat remarkable, since in the case without magnetic field the ground state energy as a function of ν tends to as ν 1 but with an infinite slope. Thus, one might expect very large variations of the eigenvalue at ν = 1 as the magnetic field is turned on, in particular one might naively expect that the ground state energy leaves the gap for small fields B. This is not the case. Moreover, since the hydrogenic Hamiltonian ceases to be selfadjoint at ν = 1 it is hard to visualize how one might arrive at such estimates using standard perturbation theory. Section 3 is devoted to the asymptotics of B(ν) as ν. We define the notion of lowest relativistic Landau level which leads to a one dimensional effective theory. This effective theory can be analyzed in great detail and allows to calculate the ground state energy λ L 1(ν,B) of the magnetic Dirac Coulomb equation (1) in the lowest relativistic Landau level. It is given by the variational problem λ L 1 (ν,b) := inf λ L [f,ν,b], f C (R,C)\{} where λ = λ L [f,ν,b] is defined by ( λ f(z) f (z) ) dz = 1+λ+νa B (z) +(1 νab (z)) f(z) dz, and R R a B (z) := B + se 1 Bs s +z ds.

3 Relativistic hydrogenic atoms in strong magnetic fields 3 The point here is that for B not too small and ν not too large (the precise bounds are given in Theorem 4), λ L 1(ν +ν 3/,B) λ 1 (ν,b) λ L 1(ν ν 3/,B). The one dimensional λ L 1 (ν,b) problem, although not trivial, is simpler to calculate than the λ 1 (ν,b) problem. As a result, in the limit as ν, this new theory yields the first term in the asymptotics of the logarithm of the critical field. In particular we have the following result, lim νlog(b(ν)) = π. ν From a methodological point of view, the ground state energy of the Dirac operator is not given by a minimum problem for the corresponding Rayleigh quotient, but it is a min max in the sense that one decomposes the whole Hilbert space H = H 1 H, maximizes the energy over functions in H and then minimizes over non-zero functions in H 1. While the choice of these Hilbert spaces is not arbitrary, there is some flexibility in choosing them, see [11, 1, 5]. For certain choices, the maximization problem can be worked out almost explicitly leading to a new energy functional for which the ground state energy λ 1 (ν,b) is the minimum. In this sense, the ground state energy of the Dirac operator appears as a minimum of a well defined functional. Both variational characterizations, the min max and the min, are of course equivalent, and our approach depends on the interplay between the two. Our results are different from the work of [] which considered the nonrelativistic hydrogen atom and worked out the asymptotics of the ground state energy as B for every ν >. In our case, however, ν has to stay in the interval [,1) in order that the operator can be defined as a selfadjoint operator. Further, the critical field is always finite and we are interested in estimating it as a function of ν. The similarity with [] comes as we let ν, since then the critical field tends to infinity but the estimates are not the same. While the mathematical methods are the main point of this paper, let us make a few additional remarks about its physical motivation. Spontaneous pair creation in strong external fields, although never experimentally confirmed, has been analyzed by Nenciu [17, 18]. In [17] it was conjectured that by adiabatically switching the potential on and off, there is spontaneous pair creation provided some of the eigenvalues emerging from the negative spectrum crossed eigenvalues emerging from the positive spectrum. This conjecture was partly proved in [18] and[].sincesuchacrossoveroccursinthedirachydrogenicatomwithastrong magnetic field, it is natural to try to estimate the strength of the magnetic fields for which this crossing phenomenon occurs. Note that the unit in which we measure the magnetic field is huge, about Tesla. Sources of gigantic magnetic fields are neutron stars that can carrymagneticfields ofabout 1 9 Tesla.Fields of1 11 Tesla for aneutron star in its gestation and in magnetars are expected, and there is speculation that fields ofupto1 1 Teslamayexistintheinteriorofamagnetar.Thereisaconsiderable literature in this area and an entertaining introduction can be found in [14]. Further, it isexpected that nearthe surfaceofaneutron staratomspersist up to about Z = 4. We show that the critical field at Z = 4 must be larger than Tesla, and preliminary calculations using numerical methods based on Landau levels yield an estimated value of about Tesla. Although improvements on these estimates are currently under investigation we believe it is unlikely that they will yield relevant values for the magnetic field strength. For elements with higher Z, the values for the critical field are much lower. In the case of Uranium (Z = 9), they are sandwiched between Tesla and an estimated value (using Landau levels) of Tesla. Speculations that large magnetic fields facilitate the creation of electron - positron pairs are not new in the physics and astrophysics literature. Clearly,

4 4 Jean Dolbeault, Maria J. Esteban, Michael Loss the Dirac operator coupled to a magnetic field but without electrostatic potential has a gap of mc independent of the magnetic field. It was pointed out in [3, 19] that the anomalous magnetic moment narrows the gap, i.e., it decreases the energy needed for pair production. In lowest nontrivial order the anomalous magnetic energy is proportional to the magnetic field which leads indeed to a narrowing of the gap; in fact the gap closes at a field strength of about Tesla. It was observed in [1], however, that the anomalous magnetic energy depends in a non linear fashion on the external field. Further it is shown that even at field strengths of 1 1 Tesla the gap narrows only a tiny bit, irrelevant for pair production. For a review of these issues the reader may consult [8]. Our contribution is to take into account simultaneously the magnetic field and the Coulomb singularity, in which case no explicit or simple calculations are possible. Of course our analysis only deals with a single electron and a fixed nucleus. A description of the non-relativistic many electron atom under such extreme situations has been given in [9,15,16]. The authors study various limits as the nuclear charge and the magnetic field strength gets large and determine the shape of the atom in these limits. In non-relativistic physics the natural scale for the magnetic field is α m c q.4 15 Tesla, much smaller than the ones under considerations in our paper. As we have mentioned before, for small Z the critical magnetic field is of the order of e π Zα and hence non relativistic physics is sufficient to explain the rough shapes and sizes of atoms even at very high field strengths. It may very well be, however, that for heavy elements and very large fields qualitatively new effects appear that cannot be understood on the basis of non relativistic physics alone. Should such effects occur, then it could make sense to treat the many body relativistic electron problem using the Dirac - Fock approximation.. Ground state and critical magnetic field In this section, we set some notations, establish basic properties and prove Estimate () on the critical magnetic field..1. Min max characterization of the ground state energy. The eigenvalue equation for the Hamiltonian (1) H B ψ ν ψ = λψ (3) is an equation for four complex functions. It is convenient to split ψ as ( ) φ ψ = χ where φ,χ L (R 3 ;C ) are the upper and lower components. Written in terms of φ and χ, (3) is given by Here P B denotes the operator where P B χ+φ ν φ = λφ, (4) P B φ χ ν χ = λχ. (5) P B := iσ ( ia B (x)), A B (x) := B x x 1

5 Relativistic hydrogenic atoms in strong magnetic fields 5 is the magnetic potential associated with the constant magnetic field B(x) :=. B Using (5) we can eliminate the lower component χ in (4). Taking then the inner product with φ we get J[φ,λ,ν,B] =, (6) where J[φ,λ,ν,B] := R 3 ( P B φ 1+λ+ ν ) +(1 λ) φ ν φ d 3 x. Thus, we see that that for any eigenvalue λ ( 1,1) of (1) the corresponding eigenvector leads to a solution of (6). Reciprocally, the functional J can be used to characterize the eigenvalues. For this purpose, a few definitions are useful. The functional J[φ,λ,ν,B] is defined for any B R +, ν (,1), λ 1 and φ C (R3,C ). Further, in order that (6) makes sense we introduce the set A(ν,B) := {φ C (R 3 ) : φ L (R 3 ) = 1, λ J[φ,λ,ν,B] changes sign in ( 1,+ )}. Note that this set might be a priori empty. Finally, since the function J is decreasing in λ, we define λ = λ[φ,ν,b] to be either the unique solution to or λ[φ,ν,b] = 1 if J[φ, 1,ν,B]. J[φ,λ,ν,B] = if φ A(ν,B), Theorem 1. Let B R + and ν (,1). If 1 < this infimum is achieved and inf λ[φ,ν,b] < 1, φ C (R3,C ) λ 1 (ν,b) := inf φ A(ν,B) λ[φ,ν,b] is the lowest eigenvalue of H B ν in the gap of its continuous spectrum, ( 1,1). Proof. This proposition is a consequence of Theorem 3.1 in [5]. The essential assumptions of this theorem are: i) The selfadjointness of H B ν 1 which is proved in the appendix. It is crucial here that < ν < 1. ii) The existence of a direct decomposition of L (R 3 ;C 4 ) as the sum of two subspaces H 1 H such that (x,(h B ν 1 )x) a := sup x H (x,x) < c 1 := inf x H 1 sup y H (x+y,(h B ν 1 )(x+y)) x+y.

6 6 Jean Dolbeault, Maria J. Esteban, Michael Loss Set b := infσ ess (H B ν 1 ) (a,+ ). If c 1 < b then c 1 is the lowest eigenvalue of H B ν 1 in the interval (a,b). In the present case we choose the decomposition ( ) ( ) ( ) φ φ ψ = = + χ χ based on the upper and lower components of the four components spinor ψ. It is easy to see that a = 1. Furthermore the essential spectrum of H B ν 1 is (, 1] [1,+ ) independently of B, see [1]. Hence b = 1. It remains to calculate the supremum in the definition of c 1 as a function of x = ( φ ). Note that the Rayleigh quotient in the definition of c 1 is strictly concave in y = ( χ). Therefore the supremum is uniquely achieved by ( χ[φ] = 1+λ[φ,ν,B]+ ν ) 1 P B φ and its value is λ[φ,ν,b], that is, λ[φ,ν,b] = sup χ C (R3,C ),ψ=( φ χ) ( (HB ) ) ν ψ,ψ. (ψ,ψ) Remarks: 1) The eigenvalue λ 1 (ν,b) can be characterized either as the minimum of the functional λ[φ,ν,b] or as a min-max level of H B ν 1. Both characterizations will be useful in the sequel of this paper. ) Under the assumptions of Theorem 1, we have J[φ,λ,ν,B] φ C (R3,C ) for any λ λ 1 (ν,b). The eigenvalue λ 1 (ν,b) can therefore be interpreted as the best constant in the above inequality. 3) When λ 1 (ν,b) is equal to 1, it belongs to the continuous spectrum and it is not necessarily an eigenvalue of H B ν 1... Basic properties of the ground state energy. Proposition 1. For all B, the function ν λ 1 (ν,b) is monotone nonincreasing on (, 1). Theproofisleft tothe reader.itisaconsequenceofthedefinition ofj[φ,λ,ν,b]. Proposition. For all B, the function ν λ 1 (ν,b) is continuous in the interval (,1) as long as λ 1 (ν,b) ( 1,1). Proof. By Theorem 1, if λ 1 (ν,b) ( 1,1) there exists a function φ ν such that J[φ ν,λ 1 (ν,b),ν,b] =. For any sequence {ν n } n converging to ν, the upper semi-continuity of ν λ 1 (ν,b) holds: lim sup n + λ 1 (ν n,b) limsupλ[φ ν,ν n,b] = λ 1 (ν,b). n + If ν n ν, then λ n := λ 1 (ν n,b) λ 1 (ν,b) and {λ n } n converges to λ. Consider therefore a {ν n } n converging to ν from above. We have to face two cases: First case: λ n > 1 for all n N. Since J[φ ν,λ 1 (ν,b),ν n,b], we know that λ n λ 1 (ν,b). Consider the corresponding eigenfunctions ψ n, such that

7 Relativistic hydrogenic atoms in strong magnetic fields 7 J[φ n,λ n,ν n,b] =, where φ n denotes the upper component of ψ n and assume that φ n L (R 3 ) = 1. By Theorem 1, we have R 3 ( ) P B φ n +(1 λ 1+λ n + νn n ) φ n d 3 ν n x = R φ n d 3 x. 3 Assume that Λ := liminf n + λ 1 (ν n,b) < λ 1 (ν,b). Up to the extraction of a subsequence, assume further that {λ 1 (ν n,b)} n N converges to some value in ( 1,λ 1 (ν,b)) and choose λ (Λ,λ 1 (ν,b)). For n large enough, λ 1 (ν n,b) < λ and ( R 3 P B φ n 1+ λ+ νn ) +(1 λ) φ n d 3 x R 3 ν n φ n. (7) Second case: λ 1 (ν,b) = 1 for all ν > ν. We choose λ ( 1,λ 1 (ν,b)) and find a {φ n } n such that φ n L (R 3 ) = 1 and J[φ n, λ,ν n,b] for n large: (7) also holds. Using the monotonicity of the {ν n } n, which implies the monotonicity of the {λ n } n by Proposition 1, and the fact that in both cases, ν (,1) and λ ( 1,1), we deduce from (7) a uniform bound for the functions φ n : sup n R 3 ( P B φ n + φ n ) d 3 x < +. (8) The proof goes as follows. It is sufficient to prove that R 3 1 φ n d 3 x is uniformly bounded. Let χ be a smooth truncation function such that χ(r) 1 if r [,1), χ(r) if r >, and χ 1. Since R3 φ n R3 d 3 φ n x d 3 x+ 1 φ n R R 3 R3 φ n d 3 x+ 1 R ( ( )) 1 χ d 3 x R with φ n = χ( R 1 )φ n, it is therefore sufficient to prove that R 3 1 φ n d 3 x is uniformly bounded, for some R >, eventually small. Using the estimate we get the following lower bound a (a+b) 1+ε b ε, R3 P B φ n R3 d νn 1+ λ+ 3 P B φ n χ x d νn 1+ λ+ 3 x R3 P B φn (1+ε) ( ) d νn 1+ λ+ 3 x C ε φ n L (R 3 ) where C is a constant which depends on χ L (R + ), B and R. Next, with the same type of arguments, we can write P B φn σ φ n 1+ε 1 ε B φ σ φ n n B R φ n. 1+ε ε

8 8 Jean Dolbeault, Maria J. Esteban, Michael Loss Collecting these estimates, this gives R3 σ φ n (1+ε) ( ) d νn 1+ λ+ 3 ν n x C(ε,R,χ)+ R φ n d 3 x. 3 Because φ n has a compact support in the ball of radius R, if δ = δ(r) > (1+ λ)r/ν n at least for n large enough, then 1+ 1 λ+ νn ν n (1+δ) x B(,R) so that 1 (1+ε) ν n (1+δ) On the other hand, according to [5,4], σ φ n d 3 x C(ε,R,χ)+ R 3 σ φ n d 3 x R 3 R 3 1 φ n d 3 x. R 3 ν n φ n d 3 x. This provides a uniform upper bound on R 3 1 φ n d 3 x if ε and δ are chosen small enough in order that 1 (1+ε) (1+δ) > ν n for n large. This can always be done since ν n converges to ν (,1) and δ(r) can be taken as small as desired for R > sufficiently small. This concludes the proof of (8). Summarizing, we obtain that [ J φ n, 1 ] ( λ+λ 1 (ν,b)),ν,b, for n large enough: hence λ[φ n,ν,b] 1 ( λ+λ 1 (ν,b)) < λ 1 (ν,b), a contradiction. Consider now the effect of a scaling on J. Lemma 1. Let B, λ 1, θ > and φ θ (x) := θ 3/ φ(θx) for any x R 3. Then Aθ B φ θ (x) = θ 5/ [ φ(θx) ia B (θx)φ(θx)], and for any a R, ν (,1), ( J[φ θ,λ,θ a ν,θ B] = θ P B φ R 3 1+λ+ θa+1 ν +(1 λ) φ θa+1 ν φ ) d 3 x. Using this scaling, we prove some properties enjoyed by the function λ 1 (ν,b). Take λ = λ[φ θ,ν,θ B] > 1, θ > 1 and a = in (9). With 1+λ = θ(1+µ), = J[φ θ,λ,ν,θ B] ( ) P B φ = θ R 1+µ+ ν +(1 µ) φ ν 3 φ d 3 x+(1 θ) φ d 3 x. R 3 (9)

9 Relativistic hydrogenic atoms in strong magnetic fields 9 Assuming that φ L (R 3 ) = 1, we get and thus, J[φ,µ,ν,B] = θ 1 θ > λ[φ,ν,b] > µ = λ θ θ 1 θ On the other hand, µ J[φ,µ,ν,B] 1, so that Summarizing, we have the estimate λ[φ θ,ν,θ B] θ λ[φ,ν,b] µ+j[φ,µ,ν,b] = λ θ + θ 1 θ θ 1 θ λ[φ,ν,b] λ[φ θ,ν,θ B] θ.. + θ 1 θ θ > 1. The above estimate, which holds provided λ[φ θ,ν,θ B] > 1 is equivalent to λ[φ,ν,θ B] θ θ 1 θ λ[φ 1/θ,ν,B] λ[φ,ν,θ B] θ + θ 1 θ θ > 1 (1) under the condition λ[φ,ν,θ B] > 1. As a consequence, we have the following result. Proposition 3. For all ν (,1), the function B λ 1 (ν,b) is continuous as long as it takes its values in ( 1,+ ). Moreover λ 1 (ν,θ B) θ θ 1 θ if λ 1 (ν,b) ( 1,+ ) and θ λ 1 (ν,b) λ 1(ν,θ B) + θ 1, (11) θ θ ( ) 1, 1 λ 1(ν,B). As a consequence, B λ 1 (ν,b) is Lipschitz continuous for any ν (,1) and B > such that λ 1 (ν,b) ( 1,+ ): λ 1 1 B λ 1 B λ 1 +1 B. Proof. Choose a ( 1,λ 1 (ν,b)) and take any φ C (R 3,C ). Since J[φ,λ,ν,B] 1, λ an integration on the interval [a,λ[φ,ν,b]] shows that [ ] λ=λ[φ,ν,b] J[φ,a,ν,B] = J[φ,λ,ν,B] λ[φ,ν,b]+a λ=a wherethefirstequalityholdsbydefinitionofλ[φ,ν,b], i.e.j[φ,λ[φ,ν,b],ν,b] =. As a consequence, J[φ,a,ν,B] λ[φ,ν,b] a >. The function θ J[φ,a,ν,θ B] is continuous, so only two cases are possible: First case: J[φ,a,ν,θ B] θ > 1, Second case: there exists a constant θ = θ(a,φ) > 1 such that (i) J[φ,a,ν,θ B] > for any θ (1, θ), (ii) J[φ,a,ν, θ B] = or, equivalently, λ[φ,ν, θ B] = a.

10 1 Jean Dolbeault, Maria J. Esteban, Michael Loss In the second case, by (i), we know that λ[φ,ν,θ B] > a > 1 for any θ (1, θ) and so (1) applies: θλ 1 (ν,b) θλ[φ 1/θ,ν,B] λ[φ,ν,θ B]+θ 1. In the limit case θ = θ, we get θλ 1 (ν,b)+1 θ λ[φ,ν, θ B] = a, which gives the estimate θ 1 a 1 λ 1 (ν,b) =: θ (a). Thus the inequality θλ 1 (ν,b) θλ[φ 1/θ,ν,B] λ[φ,ν,θ B]+θ 1 holds for any θ [,θ (a)] and for any φ C (R 3,C ), which proves the r.h.s. inequality in (11) by letting a 1: lim a 1 θ (a) = 1 λ 1 (ν,b). The l.h.s. inequality is obtained in the same manner. With the appropriate test functions one can prove that λ 1 (ν,b) is always below 1forB large.werecallthatλ 1 (ν,b) = 1meansthatifJ[φ, 1,ν,B] for any φ A(ν,B). Let us give some details. Proposition 4. Let ν (,1). Then for B large enough, λ 1 (ν,b) and there exists B > such that λ 1 (ν,b) = 1 for any B B. Proof. Let us consider B > and the trial function ( ) φ ψ =, where φ = B π e B 4 ( x1 + x ) ( ) f(x3 ) and f C (R,R) is such that f 1 for δ, δ small but fixed, and f L (R) = 1. Note that φ Ker(P B +iσ 3 x3 ) and so, P B φ = iσ 3 x3 φ, Moreover, the state φ is normalized in L (R 3 ). With r =, we can define ( r G B [φ] := ν P Bφ ν r φ ) d 3 x (1) and compute G B [φ] = B = R R + R R + R 3 ( s + x 3 ν ( B 1 s + x 3 ν f (x 3 ) ν f(x 3) s + x 3 )se Bs / dsdx 3 f (x 3 ) ν f(x 3) ) se s / dsdx 3. B 1 s + x 3

11 Relativistic hydrogenic atoms in strong magnetic fields 11 Using B 1 s + x 3 B 1/ s+ x 3 and R R+ f(x 3 ) 1 / B 1 s + x 3 se s ds ds δ e for B 1, we can therefore bound G B [φ] from above by G B [φ] C 1 ν +C ν C 3 νlogb. s B 1/ s+x 3 dx e log(δ B), where C i, i = 1,, 3, are positive constants which depend only on f. For B 1 large enough, G B [φ] + φ and λ 1 (ν,b) = 1, since in this case J[φ, 1,ν,B]..3. The critical magnetic field. Proposition 4 motivates the following definition. Definition 1. Let ν (,1). We define the critical magnetic field as { } B(ν) := inf B > : lim λ 1 (ν,b) = 1. bրb Corollary 1. For all ν (,1), λ 1 (ν,b) < 1 for any B (,B(ν)). Proof. For B = we have λ 1 (ν,) = 1 ν < 1. Given B >, small, by continuity of B λ 1 (ν,b) we know that λ 1 (ν,b) (,1). Let us consider θ (1, B(ν)/B) such that 1 < λ 1 (ν,θ B). This is made possible by Propositions 3 and 4. Then, by Proposition 3, λ 1 (ν,b) θ 1 θ < 1. The computations of Proposition 4 show the existence of a constant C 3 > such that B(ν) e C3/ν, for all ν (,1). This estimate can be made more precise for any ν not too small: Theorem. For all ν (,1), there exists a constant C > such that ( ) 4 18πν B(ν) min 5ν [3ν ], e C/ν. + The proof of this theorem uses Proposition 4. Otherwise it is splitted in two partial results stated in Propositions 5 and 7. Notice that there is a big gap between these lower and upper estimates when ν is small. To try to better understand this problem, in the next section, we will analyze the case when B is large, proving that the 3d definition of B(ν) is actually asymptotically equivalent to a 1d problem related to the lowest relativistic Landau level. More precisely we will prove that when ν is small, B is not too small and λ 1 (ν,b) ( 1,1), the eigenvalue associated with λ 1 (ν,b) ( 1,1) is almost equal to the corresponding eigenvalue in the lowest relativistic Landau level class of functions, see Theorem 4. We will also establish that B(ν) behaves in the limit ν like the upper bound in Theorem and obtain the corresponding value of C, see Theorem 5. Our first partial result is the following

12 1 Jean Dolbeault, Maria J. Esteban, Michael Loss Proposition 5. For any ν ( /3,1), B(ν) 3 πν 3ν. Proof. Consider the trial function ψ = ( φ ) where φ = ( ) 3/4 B e B /4 π ( ) 1, isliketheonechosenintheproofofproposition4,withf(x 3 ) = ( B 1/4e Bx π) 3 3 /4. Here, with the notation r :=, we find G B [φ] = B r x3 ν φ d 3 x 4ν r φ d 3 x [ = (π) 3 1 π π ] B e r / r 5 dr cos θ sinθdθ 4πν e r / rdr ν [ = (π) 3 1 π ] B e r / r 5 dr 4πν e r / rdr 3ν [ ] = (π) 3 1 8π B 3ν 4πν. If ν (/3,1)and B 3 πν 3ν, then G B[φ] = φ and so λ 1 (ν,b) = 1, which proves the Proposition. Proposition 4 shows that for ν > /3 and B large, λ 1 (ν,b) possibly ceases tobeaneigenvalueoftheoperatorh B ν 1.Thiscanbeinterpretedbysaying that for strong magnetic fields, the Coulomb potential does not stabilize the electron. At some level, electron-positron pairs could appear and then Quantum Field Theory (or QED) becomes unavoidable for a correct description of the electron dynamics, see []. Proposition 6. For given ν (,1) and B >, the function θ λ 1 (θ 1 ν,θ B) is monotone nondecreasing as long as it takes its values in ( 1,1) and ν/θ (,1). Proof. Take a = 1 in (9): J [ φ θ,λ,θ 1 ν,θ B ] R3 = (θ P B φ 1) so that for θ < 1, 1+λ+ ν J [ φ θ,λ,θ 1 ν,θ B ] J[φ,λ,ν,B] d 3 x+j[φ,λ,ν,b] at least for 1 θ >, small, so that φ θ A ( θ 1 ν,θ B ) for φ = φ ν such that λ[φ ν,ν,b] = λ 1 (ν,b). This proves that λ 1 ( θ 1 ν,θ B ) λ 1 (ν,b) (13) for 1 θ >, small. By continuation, the property holds as long as the assumptions of Proposition 6 are satisfied. The case θ > 1 follows by multiplying θν and θ B by respectively θ 1 and θ. Corollary. There exists a positive constant Λ such that B(ν) Λ ν as ν ց.

13 Relativistic hydrogenic atoms in strong magnetic fields 13 Proof. Let (ν,b ) be such that B(ν ) > B, i.e., λ 1 (ν,b ) > 1 and take ν (,ν ), θ = ν/ν (,1), B = θ B in (13): 1 < λ 1 (ν,b ) λ 1 (ν,b) = λ 1 ( B /Bν,B). By Proposition 1, this inequality can be extended to This amounts to say that λ 1 ( B /Bν,B) λ 1 (ν,b) ν (, B /Bν ). ν B /Bν = B(ν) B ν ν, which proves the result with Λ = B ν. The constant Λ can be made more precise. The remainder of this section is devoted to the following improvement of Corollary. Proposition 7. For all ν (,1), G B [φ] = R ν P Bφ d 3 ν x 3 R φ d 3 x ν 5B φ d 3 x. 3 R 3 In particular this implies that B(ν) 4 5ν. Proof. Scaling the function φ according to ( ) φ B := B 3/4 φ B 1/ x preserves the L norm, and yields G B [φ B ] = BG 1 [φ], where G B has been defined in (1). Obviously it is sufficient to find a good estimate on the functional G 1. Let us collect some preliminary observations. Recalling that the angular momentum vector L is given by a simple calculation shows that L = i x, ( σ = σ x ) ( r 1 ) r r σ L with r = and r = x r. We also recall that ( ia 1 (x) σ = σ x ) r (σ q(x)) where q(x) = 1 r x 3x 1 x 3 x x 1 +x and ip 1 = σ ia 1 (x) σ

14 14 Jean Dolbeault, Maria J. Esteban, Michael Loss so that we can expand P 1 φ as P 1 φ = ( r 1r ) σ L σ q(x) φ = r φ + 1 r σ Lφ + q φ r φ,σ qφ + φ,σ ( r q)φ + 1 [ ] r φ,σ Lφ + σ Lφ,σ qφ + σ qφ,σ Lφ. r As a last preliminary remark, we notice that r r q = q. Since the vector potential grows linearly, we localize the problem near the origin. To this end consider the function 1 if r R, t(r) = R/r if r R. Since t(r) 1 and ν 1 we get the lower bound t(r)r G 1 [φ] P 1 φ d 3 ν x R ν 3 R r φ d 3 x ( 3 ) ν t(r)r P 1 φ d 3 1 x R 3 R r φ d 3 x 3 ( ) 1 = ν K[φ] R r φ d 3 x 3 where the kinetic part is defined by K[φ] := R 3 t(r)r P 1 φ d 3 x and satisfies K[φ] = t(r)r [ r φ + 1r ] R σ Lφ + q φ d 3 x 3 [ ] + t(r) r r φ,σ qφ +r φ,σ ( r q)φ r φ,σ Lφ d 3 x R 3 [ ] + t(r) σ Lφ,σ qφ + σ qφ,σ Lφ d 3 x. R 3 An integration by parts in the r variable yields K[φ] = t(r)r [ r φ + 1r ] R σ Lφ + q φ d 3 x 3 + t(r) [4 φ,σ qφ + r ] φ,σ Lφ + σ Lφ,σ qφ + σ qφ,σ Lφ d 3 x R 3 [ ] + t (r) r φ,σ qφ + φ,σ Lφ d 3 x, R 3 where we have also used that r r q = q. Consider now the region of integration where r R and denote the corresponding expression by K 1 [φ]. There the derivative of t(r) vanishes and hence collecting terms we find K 1 [φ] = r R [ [ tr r φ + 1 r ] (σ L+1)+σ q φ ]d 1 r φ 3 x + φ,σ qφ d 3 x. r R

15 Relativistic hydrogenic atoms in strong magnetic fields 15 At this point we have decoupled the derivatives with respect to r from the magnetic field or, to be precise, from q. The problem is that the angular momentum is still coupled to the magnetic field. Obviously [ 1 ] r (σ L+1)+σ q φ 1 σ qφ r (σ L+1)φ. Further, since the eigenvalues of σ L+1 are given by ±1, ±..., and we have that [ 1 ] r (σ L+1)+σ q φ (σ L+1)φ L (S ) φ L (S ), L (S ) 1 [ r φ 1 ] L (S ) q φ L (S ) r r φ L (S ), since q(r) r/. For r R, the factor [ 1 r r ] is nonnegative. Since [ 1 r r ] 1 r = r 4 1 and since q(r) r, we obtain the lower bound K 1 [φ] r R [ [ r r r φ ] + 4 φ ]d 3 x. The function r3 4 r is a decreasing function on the interval [, ]. Hence [ R K 1 [φ] r r φ d 3 3 x+ r R 4 ]r R R φ d 3 x, if R <. Next, we look at the contribution to K[φ] of the region where r R, which we denote by K [φ], K [φ] = t + r R r R r R t(r)r t(r) t(r) r [ r φ + 1r σ Lφ + q φ ] d 3 x [4 φ,σ qφ + r φ,σ Lφ + σ Lφ,σ qφ + σ qφ,σ Lφ ] d 3 x [ ] r φ,σ qφ + φ,σ Lφ, using the fact that t = t/r. Collecting the terms, we get K [φ] = t(r)r [ r φ + 1r ] σ Lφ + q φ d 3 x r R + t(r) [3 φ,σ qφ + 1r ] φ,σ Lφ + σ Lφ,σ qφ + σ qφ,σ Lφ d 3 x. r R This can be rewritten as [ K [φ] t(r)r r φ + 1 [ [ σ L+ 1 ] ] +σ q φ r R r + t(r) φ,σ qφ d 3 x. r R ] 1 4r φ d 3 x

16 16 Jean Dolbeault, Maria J. Esteban, Michael Loss Finally we get K [φ] r R and, using q /r 1/, K [φ] t(r)r r φ d 3 x 1 φ d 3 1 x R 4R r R r R r q φ d 3 x, r R [ t(r)r r φ d 3 x R+ 1 ] φ d 3 x. 4R r R Thus we can estimate K[φ] = K 1 [φ]+k [φ] as follows: [ R K[φ] t(r)r r φ d 3 3 x+ R 4 ]r R R 3 Observe that generally [ φ d 3 x R+ 1 ] φ d 3 x. 4R r R t(r)r rφ+ 1 R r φ d 3 x 3 = t(r)r r φ d 3 x t(r) 1 R 3 R r φ d 3 x t (r) φ d 3 x. 3 R 3 Since t (r) and t(r) 1 on [,R), and t (r) = t(r)/r on (R, ), the following estimate holds t(r)r r φ d 3 1 x R 3 R r φ d 3 1 t(r) x 3 R r 3 = and hence r R φ d 3 x+ 1 r φ d 3 x 1 R t (r) φ d 3 x R 3 φ d 3 x 1 [ R K[φ] R r φ d 3 3 [ x 4 ]r R R φ d 3 x R+ 5 ] φ d 3 x. 4R 3 r R Optimizing on R (, ], i.e. using { [R 3 ] [ max R (, min ) 4 R, R 5 ] } [ = R 5 ] 4R 4R R= = 5, 5/ r R we get ( ) 1 G 1 [φ] ν K[φ] R r φ d 3 x ν 5 φ d 3 x. 3 R 3 Hence the condition G B [φ] = BG 1 [φ] φ entails that B(ν) 4 5ν.

17 Relativistic hydrogenic atoms in strong magnetic fields Asymptotics for the critical magnetic field In the large magnetic field limit, the upper component of the eigenfunction corresponding to the lowest energy levels in the gap of Dirac operator with magnetic field H B ν 1 is expected to behave like the eigenfunctions associated to the lowest levels of the Landau operator L B := iσ 1 x1 iσ x σ A B (x), which can alsowewritten asl B = P B +iσ 3 x3 or L B = i( x1 +ibx /)σ 1 i( x ibx 1 /)σ. The goalofthis section is tocomparethe lowestenergylevels of H B ν 1 with its lowest energy levels on a space generated by the lowest energy levels of L B. The asymptotic analysis for the small coupling limit ν + isnotthatsimplebecausethelandaulevelsarenotstableundertheactionofthe kinetic part of the Dirac Hamiltonian. The way out is to choose a representation of H B ν 1 that diagonalizes the kinetic energy in the Dirac Hamiltonian and to project both upper and lower components on the lowest Landau levels Projection on Landau levels. To start with, we observe that P B = L B x 3 and summarize the basic properties of the lowest energy levels of L B. Lemma. [[1], Section 7.1.3] The operator L B in L (R,C ) has discrete spectrum {nb : n N}, each eigenvalue being infinitely degenerate. Moreover the kernel of this operator, that is, the eigenspace corresponding to the eigenvalue, is the set generated by the L -normalized functions φ l := B(l+1)/ πl l! (x +ix 1 ) l e Bs /4 ( 1 ), l N, s = x 1 +x. Next we diagonalize the free magnetic Dirac Hamiltonian. First we write it in the form ( ) ( ) I PB R Q K B = = I+P P B I B, Q R where R and Q are operators acting on spinors, given by R = 1 P B, Q =, I+P B I+P B and satisfy the relation The matrix R +Q = I. ( ) R Q Q R is a reflection matrix and hence has eigenvalues 1 and 1. It can be diagonalized using the matrix ( ) 1 Q R I U =. (I R) I R Q The operator defined by U is unitary and such that ( ) U I+P K B U = B I+PB.

18 18 Jean Dolbeault, Maria J. Esteban, Michael Loss The potential V = 1 r is transformed into the nonnegative operator ( ) P := U p q VU = q. t Here and from now on, we will omit I whenever it is multiplied by a scalarvalued function. If we denote by Z any 4-spinor and decompose it as ( ) X Z =, Y where X and Y are the upper and lower components, in the new representation, the full magnetic Dirac Hamiltonian takes the form ( ) ( ) U H B U = U I+P K B U νp = B p q I+PB ν q. t The Dirac energy for an electronic wave function Z in the electromagnetic potential (V,A) is now K[Z] := (Z,U K B UZ) = E ν [Z] := K[Z] ν(z,p Z), (X, ) I+P B X (Y, ) I+P B Y As we shall see below, in the new representation, restricting the upper and lower components X and Y to the lowest Landau levels makes sense for studying the regime of asymptotically large B. The price we pay for that is that all quantities like R, Q, U, P... depend on B. Denote by Π L the projector on the lowest Landau level, whose image is generated by all functions (x 1,x,x 3 ) φ l (x 1,x )f(x 3 ) l N, f L (R), and define ΠL c := I Π L. Notice that Π L commutes with L B. With the above notations, for all ξ L (R 3,C ), we have that ( ) ( ) ( ) ξ, I+P B ξ = Π L ξ, I+PB Π Lξ + ΠLξ, c I+PB Πc Lξ. (14) Next, we decompose any Z (L (R 3,C)) 4 as Z = ΠZ +Π c Z,. where Π := ( ) ( ) ΠL, Π Π c Π c := L L ΠL c. 3.. Main estimates. From (14), it follows that K[Z] = K[Π Z]+K[Π c Z]. Since the operator P is nonnegative, we also have (Z,P Z) ( 1+ ν ) ( (Π L Z,P Π Z) ) (Π c Z,P Π c Z), ν (Z,P Z) ( 1 ν ) ( (Π L Z,P Π Z)+ 1 1 ) (Π c Z,P Π c Z). ν

19 Relativistic hydrogenic atoms in strong magnetic fields 19 This simply follows from the identities [ (1+ (a+b) a +b ) + ab = inf ν a +( 1+ 1 ) ] b, ν> ν [ (1 (a+b) a +b ) ab = sup ν a +( 1 1 ) ] b. ν> ν The above remarks prove the next proposition. Proposition 8. For all Z C (R3,C 4 ], E ν+ν 3/[ΠZ]+E ν+ ν [Π c Z] E ν [Z] E ν ν 3/[ΠZ]+E ν ν [Π c Z]. The following result will allow us to get rid of the higher Landau levels when looking for the ground state energy, i.e. of the term E ν+ ν [Π c Z] in Proposition 8. Consider ν (,1) such that ( ν + ν) =, i.e. ν.56, and for any ν (, ν), define d(δ) := (1 δ) δ, d ± (ν) := d(δ ± (ν)) with δ ± (ν) := ν ±ν. We have d(δ) > δ < 1 /, d ± (ν) > if ν < ν. Proposition 9. Let B > and δ (1 /). For any Z = ( X ), Z = ( Y), X,Y L (R 3,C ) E δ [Π c Z] d(δ) B Π c L X L (R 3 ). E δ [Π c Z] d(δ) B Π c L Y L (R 3 ). Proof. An elementary computation shows that P B = ( ia B) I+σ B. Using the diamagnetic inequality (see [1]), ( ia B )ψ dx Hardy s inequality, R 3 R 3 R 3 ψ dx, ψ dx 1 R3 ψ 4 dx, and the nonnegativity of σ B+BI, we get P B +(1+B)I = ( ia B ) I+σ B+(1+B) ( 1 4 ) 1 +1 Since the square root is operator monotone, we have B + I+PB PB +(1+B)I

20 Jean Dolbeault, Maria J. Esteban, Michael Loss Now, for any δ >, I+P B δ = (1 δ) I+PB +δ I+PB δ (1 δ) I+PB +δ δ B δ (1 δ) I+PB δ B. On the range of ΠL c the operator I+PB I+L B is bounded from below by 1+B. Hence I+PB δ (1 δ) 1+B δ B, which is equivalent to I+P B δp (1 δ) 1+B δ B. Since as B, and (1 δ) 1+B δ B d(δ) B inf [(1 δ) 1+B δ B] > B> if d(δ) >, the right hand side is positive for any field strength provided that δ < 1 / The restricted problem. Nextweprovethatthegroundstateenergyλ 1 (ν,b) is comparable with the one obtained by restricting it to the lowest Landau level, both in the upper and in the lower components of the wave function, provided that the Coulomb potential is slightly modified. By a result similar to Theorem 1 (also see Theorem 3.1 in [5] in case B = ), for all < B < B(ν), ν (,1), λ 1 (ν,b) is characterized as λ 1 (ν,b) = inf X C (R3,C ) X sup E ν [Z]. (15) Y C (R3,C ) Z L (R 3 =1, ) Z=( X Y) With the notation Z = ( X Y), we define the restricted min-max problem λ L 1(ν,B) := inf X C (R3,C ) Π c L X=, < X L (R 3 ) <1 sup Y C (R3,C ), Z=( X Y) Π c L Y =, Y L (R 3 ) =1 X L (R 3 ) E ν [Z]. We show below that this restricted problem is actually a one-dimensional problem. For this purpose, let us define the function a B : R R+ R given by a B (z) := B + se Bs / s +z ds. and implicitly define µ L as the unique solution of ( f (z) ) R 1+µ L [f,ν,b]+νa µ L [f,ν,b] = B (z) +(1 νab (z)) f(z) dz. f(z) dz R

21 Relativistic hydrogenic atoms in strong magnetic fields 1 Theorem 3. For all B > and ν (,1), λ L 1 (ν,b) = inf µ L [f,ν,b]. f C (R,C)\{} Proof. The definition of λ L 1(ν,B) is equivalent to λ L 1(ν,B) = inf sup φ C (R3,C ) Π L c φ=, φ χ C (R3,C ) Π c L χ=, ψ=(φ χ) ψ L (R 3 ) =1 1 ν r iσ 3 x3 iσ 3 x3 ψ,ψ 1 ν r or λ L 1(ν,B) = inf sup φ C (R3,C ) Π L φ χ C (R3,C ) Π Lχ, ψ=( φ χ) Π ψ L (R 3 =1 ) 1 ν r iσ 3 x3 iσ 3 x3 Π ψ,πψ 1 ν r with the notation r = = x 1 +x +x 3. For any given φ such that Πc L φ =, the supremum in χ is achieved by the function ( χ L [φ] = V L + ( λ L [φ,ν,b]+1 ) ) 1ΠL Π L ( iσ 3 x3 )φ, ν with V L (x) := Π L r Π L and λ L [φ,ν,b] := sup 1 ν r χ C (R3,C ) iσ 3 x3 Π L χ, ψ=( χ) φ iσ 3 x3 Π ψ Π ψ,. 1 ν Πψ L (R 3 ) Πψ L (R 3 ) r Since σ 3 σ3 = σ 3 = I, this yields the expression [ λ L 1 (ν,b) = inf φ C (R3,C ) Π c L φ=, φ L (R 3 ) =1 R 3 x3 φ V L +λ L [φ,ν,b]+1 + (1 ν r ] ) φ dx. Now, with the notations of Lemma, for all l l, for all h : R R, we have R φ l h(r)φ l dx 1dx =, with s = x 1 +x, r = s +x 3, and ( ( λ L 1(ν,B) = inf inf (V L +λ L [φ,ν,b]+1) 1 x3 φ, x3 φ)+ 1 ν ) ) φ,φ. l N φ L (R 3 =1 ) r φ C (R,span(φ l )) A simple calculation shows that for any g C (R,C ), ( (VL )g,g ) L (R,C ) = ( (g,φ l ) L (R,C ) φ l, ν ) r φ l l L (R,C ) a.e. in R x 3, and also that ( φ l, 1 ) r φ l = a B l (x 3) := Bl+1 L (R,C ) l l! + s l+1 e Bs / s +x 3 A simple integration by parts shows that for all l, a B l a B l 1 a.e. When minimizing, only the l = component has therefore to be taken into account. ds.

22 Jean Dolbeault, Maria J. Esteban, Michael Loss Corollary 3. For all ν (,1), the function [,+ ) B λ L 1 (ν,b) is nonincreasing in B. Proof. A simple change of variables shows that a B (z) = + se s / s B +z ds. By Theorem 3 and according to the definition of µ L [f,ν,b], this implies the monotonicity of λ L 1 (ν, ) in [, ). Proposition 1. For all B, the function ν λ L 1 (ν,b) is continuous in the interval (,1) as long as it takes its values in ( 1,1). Moreover, for any ν (,1), as long as λ L 1 (ν,b) takes its values in ( 1,1), there exists a function Z Range(Π) with Z L (R 3,C 4 ) = 1 such that E ν [Z] = λ L 1 (ν,b). The proof is similar to the one of Proposition. The above proposition enables us to define B L (ν) := inf{b > : λ L 1(ν,B) = 1}. Recall that lim ν + d + (ν) = 1/. We are now ready to state and prove the main result of this section Asymptotic results. Theorem 4. Let ν (, ν). For any B ( 1/d + (ν),min { B(ν),B L (ν+ν 3/ ) }), we have ( ) ν +ν 3/,B λ 1 (ν,b) λ L 1(ν ν 3/,B). λ L 1 Notice that the right hand side inequality holds for any B ( 1/d + (ν),b(ν) ). Proof. To prove the upper estimate, we use (15) and notice that λ 1 (ν,b) inf X C (R3,C ) Π c L X=,Π L X sup E ν (Z), Y C (R3,C ) Z L (R 3 =1, ) Z=( X Y) since adding the condition ΠL c X = increases the value of the infimum. Then, by Propositions 8 and 9, λ 1 (ν,b) inf X C (R3,C ) Π c L X=,Π L X inf X C (R3,C ) Π c L X=,Π L X inf X C (R3,C ) Π c L X=,Π L X sup E ν ν 3/[Π Z]+E ν ν [Π c Z] Y C (R3,C ) E ν ν 3/[Π Z] d (ν) B ΠL c sup Y L (R 3 ) Y C (R3,C ) ΠZ L (R 3 ) + Πc L Y L (R 3 ) E sup ν ν 3/[Π Z] Y C (R3,C ) ΠZ = λ L 1 (ν ν3/,b). L (R 3 ) Next, we establish the lower bound. By taking in (15) a smaller maximizing class of functions we decrease the maximum: λ 1 (ν,b) inf X C (R3,C ) X sup E ν [Z]. Y C (R3,C ), Z=( X Y) Π L c Y =, Z L (R 3 ) =1

23 Relativistic hydrogenic atoms in strong magnetic fields 3 Therefore, by Propositions 8 and 9, λ 1 (ν,b) inf X C (R3,C ) X ( sup Y C (R3,C ), Z=( X Y) Π L c Y =, Z L (R 3 ) =1 ) E ν+ν 3/[Π Z]+E ν+ ν [Π c Z], inf X C (R3,C ) X E ν+ν 3/[Π Z]+d + (ν) B ΠL c sup X L (R 3 ) Y C (R3,C ), Z=( ΠZ X Y) L (R 3 ) + Πc L X L (R 3 ) Π L c Y =, Z L (R 3 ) =1 Let us now notice that for every X C (R3,C ), X, E sup ν+ν 3/[Π Z] Y C (R3,C ) ΠZ = λ L [X,ν +ν 3/,B] L Z=( (R 3 ) X Y) is uniquely achieved at some Y L [X] because of the same concavity argument as in the proof of Theorem 1, after noticing that for any B (,B L (ν + ν 3/ )), 1 < λ L 1 (ν +ν3/,b) λ L [X,ν +ν 3/,B] and V L. Recall that λ L 1(ν +ν 3/,B) = inf λ L [X,ν +ν 3/,B]. X C (R3,C ) X Denoting Z L [X] = ( X Y L[X]), for any given X, we find. E ν+ν 3/[Π Z]+d + (ν) B ΠL c sup X L (R 3 ) Y C (R3,C ) ΠZ L (R 3 ) + Πc L X L (R 3 ) E ν+ν 3/(Π Z L[X])+d + (ν) B Π c L X L (R 3 ) ΠZ L [X] L (R 3 ) + Πc L X L (R 3 ) = λ LX,ν +ν 3/,B] ΠZ L [X] L (R 3 ) +d +(ν) B Π c L X L (R 3 ) ΠZ L [X] L (R 3 ) + Πc L X L (R 3 ) λ L 1 (ν +ν3/,b) for B large enough so that d + (ν) B λ L 1 (ν +ν3/,b). As we shall see below, this is always possible. Note indeed that on (, ν), d + (ν). Hence, d + (ν) 1/. Now, by monotonicity (see Corollary 3), λ L 1 (δ,b) λl 1 (δ,1/) for all δ (, ν + ν) = (,1 /) and for all B 1/. Moreover, one can prove very easily that λ L 1(δ,1/) 1 for all δ (,1). Indeed, by Theorem 3, for all B, λ L 1(δ,B) = inf f µ L [f,ν,b]. A simple scaling argument shows that we can make R f dz as small as we wish while keeping R f dz constant. Taking into account the definition of µ L [f,ν,b], this shows that λ L 1 (δ,b) 1 for all B, for all δ (,1). Therefore, for all B d + (ν), d + (ν) B λ L 1 (ν +ν3/,b) holds true. From Theorem 4, we deduce the following Corollary 4. Let ν (, ν). Then B L (ν +ν 3/ ) B(ν) B L (ν ν 3/ ). (16) From Theorem 3 and Corollary 4, better estimates of the critical magnetic strength B(ν) than those of Theorem can be established for ν small.

24 4 Jean Dolbeault, Maria J. Esteban, Michael Loss Theorem 5. The critical strength B(ν) satisfies: lim νlogb(ν) = π. ν Proof. Because of Corollary 4, B(ν) can be estimated using B L (δ) with δ = ν ±ν 3/. This amounts to look for the smallest positive B for which ( f (z) ) λ L (δ,b) := 1+ δa B (z) δab (z) f(z) dz inf f C (R,C)\{} f L (R) =1 R is such that Using the identity λ L (δ,b) = 1. a B (z) = Ba 1 ( Bz ), by the changes of variable and function y(z) := z a 1 (t)dt, f one tranforms the above minimization problem into λ L (δ,b) 1 = B Hence, inf g C (R,C)\{} R g(y) dµ(y)=1 ( z B ) = B 1/4 g(y), R ( 1 δ g (y) δ g(y) ) dy. λ L (δ,b) = 1+ B(λ L (δ,1) 1). (17) For a given δ, let κ = κ(δ) := δ ( 1 λ L (δ,1) ) and µ(y) := 1/a 1 (z(y)). The problem is reduced to look for the first eigenvalue E 1 = E 1 (δ) of the operator y +κ(δ)µ(y), namely to find δ such that The function a 1 satisfies δ = E 1 (δ). a 1 (z) a 1 () = π z R, a 1 (z) 1 as z. z There exists therefore a constant c > such that µ(y) ce y for any y R. To get an upper estimate of E 1, we may now consider the function g 1 (y) := cos(πy/) on ( 1,1) y and the rescaled functions g σ (y) := σ 1/ g 1 ( /σ). 1 1 E 1 (δ) π 4σ +κ g 1 µ(σy)dy π 1 4σ +κc e σ y g 1 dy 1 π 4σ +κc(eσ 1). Optimizing in σ in the above expression, we choose σ = σ(δ) satisfying π = 4κcσ 3 e σ, which implies σ(δ) logκ =: σ δ as δ.

25 Relativistic hydrogenic atoms in strong magnetic fields 5 A Taylor expansion at next order shows that σ(δ) σ δ 3log ( σ δ ), which yields E 1 (δ) π 4σδ (1+o(1)). Next, in order to obtain a lower estimate of E 1 (δ), we consider the function µ δ which is equal to in the interval ( σ δ,σ δ ) and equal to κµ(σ δ ) elsewhere. The function µ δ is positive, but below the function κµ(y). Then observe that the first eigenvalue of the operator y +µ δ(y), that we denote E1 δ, satisfies the equation and as δ goes to + this implies ( ) E1 δσ κµ(σ δ ) E1 δ δ = arctan E1 δ Summarizing, what we have obtained is E 1 (δ) E1 δ = π ( ) 1+o(1). 4σδ π E 1 (δ) = 4 ( logκ(δ) ) ( ) 1+o(1). So, imposing E 1 (δ) = δ, we get κ(δ) = e π δ (1+o(1)). Since by (17),, we get 1 = λ L ( δ,bl (δ) ) = 1+ B L (δ) κ(δ) δ B L (δ) = 4δ e π δ (1+o(1)), which, together with (16), concludes the proof., Appendix 1 : selfadjointness of H B ν 1 With the notation P = i, consider a Dirac operator of the form H = α P +mβ +V (x) defined on (C (R3 )) 4. We assume for instance that (V ) i,j (x) L loc (R3 ). If f C (R3 ), then the following identity holds on (C (R3 )) 4 : H f f H iα f =. (18) Still denote by H the a selfadjoint extension of α P+mβ+V (x) with domain D and let f C (R3 ) be such that f(x) 1, f(x) = 1 for 1, f(x) = for, and f a (x) := f(x/a). If } D := {g f a ψ a 1, ψ D, 1 The following proof was explained to us by George Nenciu to whom we are grateful.

26 6 Jean Dolbeault, Maria J. Esteban, Michael Loss then in all interesting cases, including the case of local Coulomb singularities with ν < 1 (see the characterisation of D in [13]), one has D D, and then by density, Identity (18) implies that for all ψ D, ( H f a f a H iα f a ) ψ =. (19) A standard characterization of selfadjoint operators that we are going to use is the following: Let T be a closed symmetric operator. T is selfadjoint if and only if Ker(T ± i) = {}. The first remark is that H is essentially selfadjoint on D. Suppose indeed that Ψ,(H ± i)(f a ψ) = a 1, ψ D for some Ψ (L (R 3 )) 4. Now from (19), Ψ,(H ± i)f a ψ = Ψ,(H ± i)ψ (1 f a )Ψ,(H ± i)ψ +i Ψ,(α f a )ψ, so that, taking a at fixed ψ, one obtains and then Ψ =. This result also applies to the case Ψ,(H ± i)ψ = ψ D, V = ν, ν < 1. See [13,17] for more details. Define as the matrix norm. Suppose that V 1 is locally L and more precisely satisfies M(R) := sup V 1 (x) < R R +, R lim M(R) =. R Consider on D the operator H = H +V 1. Lemma 3. Under the above assumptions, H is essentially selfadjoint on D. Proof. Let χ R be the characteristic function of the set {x R}. Since χ R V 1 is bounded, H R = H +χ R V 1 is essentially selfadjoint on D. Suppose that there exists Ψ (L (R 3 )) 4 such that Ψ,(H ±i)(f a ψ) = a 1, ψ D and assume that Ψ L (R 3 ) = 1. Using (18), one deduces that < f a Ψ,(H±i)ψ = i Ψ,α f a ψ. Observe now that < f a Ψ,(H±i)ψ = f a Ψ,(H a ±i)ψ, which amounts to f a Ψ,(H a ±i)ψ = i Ψ,α f a ψ. () Since H a is essentially selfadjoint on D, Range(H a ±i) is dense in D and there exists ψ ± a D such that with δ a L (R 3 ) 1/a. Also notice that (H a ±i)ψ ± a = f a Ψ +δ a (1) (H a ±i)ψ ± a L (R 3 ) = H aψ ± a L (R 3 ) + ψ± a L (R 3 ) ψ± a L (R 3 ).

27 Relativistic hydrogenic atoms in strong magnetic fields 7 From () written for ψ ± a and (1), we get f a Ψ L (R 3 ) + f aψ,δ a = i Ψ,α f a ψ ± a so that f a Ψ L (R 3 ) 1 a + f a L (R 3 ) ψ ± a L (R 3 ) 1 a + f a L (R 3 ) (H a ±i)ψ a ± L (R 3 ) 1 ( a + f a L (R 3 ) 1+ 1 ). a For a, f a Ψ L (R 3 ) 1 and f a L (R 3 ), a contradiction: Ψ =. Acknowledgments. We are grateful to George Nenciu for showing us the proof of self-adjointness in the appendix, to P. Pickl and D. Dürr for references. Some of this research has been carried out at the Erwin Schrödinger Institute. M.L. would like to thank the Ceremade for its hospitality. J.D. and M.J.E. acknowledge support from ANR Accquarel project and European Program Analysis and Quantum HPRN-CT # -77. M.L. is partially supported by U.S. National Science Foundation grant DMS c 6 by the authors. This paper may be reproduced, in its entirety, for non-commercial purposes. References 1. J. Avron, I. Herbst, and B. Simon, Schrödinger operators with magnetic fields. I. General interactions, Duke Math. J., 45 (1978), pp J. E. Avron, I. W. Herbst, and B. Simon, Schrödinger operators with magnetic fields. III. Atoms in homogeneous magnetic field, Comm. Math. Phys., 79 (1981), pp V. Canuto and H.-Y. Chiu, Thermodynamic properties of a magnetized fermi gas, Phys. Rev., 173 (1968), pp J. Dolbeault, M. J. Esteban, M. Loss, and L. Vega, An analytical proof of Hardy-like inequalities related to the Dirac operator, J. Funct. Anal., 16 (4), pp J. Dolbeault, M. J. Esteban, and E. Séré, On the eigenvalues of operators with gaps. Application to Dirac operators, J. Funct. Anal., 174 (), pp J. Dolbeault, M. J. Esteban, and E. Séré, A variational method for relativistic computations in atomic and molecular physics, Int. J. Quantum Chemistry, 93 (3), pp J. Dolbeault, M. J. Esteban, E. Séré, and M. Vanbreugel, Minimization methods for the one-particle dirac equation, Phys. Rev. Letters, 85 (), pp R. C. Duncan, Physics in ultra-strong magnetic fields. Review talk given at Fifth Huntsville Gamma-Ray Burst Symposium. 9. I. Fushiki, E. H. Gudmundsson, C. J. Pethick, and J. Yngvason, Matter in a magnetic field in the Thomas-Fermi and related theories, Ann. Physics, 16 (199), pp M. Griesemer, R. T. Lewis, and H. Siedentop, A minimax principle for eigenvalues in spectral gaps: Dirac operators with Coulomb potentials, Doc. Math., 4 (1999), pp (electronic). 11. M. Griesemer and H. Siedentop, A minimax principle for the eigenvalues in spectral gaps, J. London Math. Soc. (), 6 (1999), pp B. Jancovici, Radiative correction to the ground-state energy of an electron in an intense magnetic field, Phys. Rev., 187 (1969), pp M. Klaus and R. Wüst, Characterization and uniqueness of distinguished selfadjoint extensions of Dirac operators, Comm. Math. Phys., 64 (1978/79), pp C. Kouveliotou, R. C. Duncan, and C. Thompson, Magnetars, Scientific American Magazine, (3), pp E. H. Lieb, J. P. Solovej, and J. Yngvason, Asymptotics of heavy atoms in high magnetic fields. I. Lowest Landau band regions, Comm. Pure Appl. Math., 47 (1994), pp , Asymptotics of heavy atoms in high magnetic fields. II. Semiclassical regions, Comm. Math. Phys., 161 (1994), pp

28 8 Jean Dolbeault, Maria J. Esteban, Michael Loss 17. G. Nenciu, Self-adjointness and invariance of the essential spectrum for Dirac operators defined as quadratic forms, Comm. Math. Phys., 48 (1976), pp , Existence of the spontaneous pair creation in the external field approximation of Q.E.D., Commun. Math. Phys., 19 (1987), pp R. F. O Connell, Effect of the anomalous magnetic moment of the electron on spontaneous pair production in a strong magnetic field, Phys. Rev. Lett., 1 (1968), pp P. Pickl, Existence of Spontaneous Pair Creation, PhD thesis, Mathematisches Institut der Ludwig-Maximilians-Universität München, B. Thaller, The Dirac equation, Texts and Monographs in Physics, Springer-Verlag, Berlin, 199.