The ground-state electron density of a polaron bound to a Coulomb potential in a homogeneous magnetic field—the transverse coordinates integrated out—converges pointwise and weakly in the strong magnetic field limit to the square of a hyperbolic secant function.

A non-relativistic hydrogen atom in a strong magnetic field interacting with the quantized longitudinal optical modes of an ionic crystal is considered within the framework of Fröhlich’s 1950 polaron model.26 Starting with Platzman’s variational treatment in 1962, the polaron hydrogenic atom has been of interest for describing an electron bound to a donor impurity in a semiconductor.49 Its first rigorous examination, however, came much later in 1988 from Löwen who disproved several longstanding claims about a self-trapping transition.45 

A study of the polaron hydrogenic atom in strong magnetic fields was initiated by Larsen in 1968 for interpreting cyclotron resonance measurements in InSb.35 The model has since been considered in formal analogy to the hydrogen atom in a magnetic field though the latter was understood rigorously again much later in 1981 by Avron et al. who proved several properties including the non-degeneracy of the ground state.3 Whether or not these hydrogenic properties indeed persist when a coupling to a quantized field is turned on remains to be seen.

Polarons are the simplest quantum field theory models, yet their most basic features such as the effective mass, ground-state energy, and wave function cannot be evaluated explicitly. And, quite unfortunately, while several successful theories have been proposed over the years to approximate the energy and effective mass of various polarons, usually they are built entirely on unjustified, even questionable, Ansätze for the wave function (see Refs. 15, 26, and 48); the only exceptions are the rigorous treatments1,2,6–8,13,17,19–23,27,30,42,43,45–47,53,57 including the recent work on the effective mass.14,38,39 The paper provides now for the first time an explicit description of the wave function.

For the polaron hydrogenic atom in a homogeneous magnetic field, its ground-state electron density in the magnetic-field direction is shown to converge in the strong field limit to the square of a hyperbolic secant function: a sharp contrast to the paradigmatic Gaussian variational wave functions (see Ref. 58 and also Ref. 54 and references therein). The explicit limiting function is realized as a density of the minimizer of a one-dimensional problem with a delta-function potential describing the second leading-order term of the ground-state energy (cf. Refs. 5, 17, and 41).

The Fröhlich model is defined by the Hamiltonian

H(B)HB32βx1+N+α2πR3akeikxk+akeikxkdk
(2.1)

acting on the Hilbert space H:=L2(R3)F, where F:=n0snL2(R3) is a symmetric phonon Fock space over L2(R3). The creation and annihilation operators for a phonon mode ak and ak act on F and satisfy [ak,ak]=δ(kk). The energy of the phonon field is described by the operator N=R3akakdk. The kinetic energy of the electron is described by the operator HB32 acting on L2(R3), where HB=j=1,2ij+Ajx2 is the two-dimensional Landau Hamiltonian with the magnetic vector potential Ax1,x2,x3=B/2x2,x1,0 corresponding to a homogeneous magnetic field of strength B ≥ 0 in the x3-direction; the transverse coordinates are denoted by x = (x1, x2). Furthermore, inf spec HB = B. The parameters α ≥ 0 and β > 0 denote the strengths of the Coulombic electron–phonon interaction and the localizing Coulomb potential; the coupling function 1/k is proportional to the square root of the Fourier transform of the Coulomb interaction. The ground-state energy is

E0(B)infΨ,H(B)Ψ:Ψ=1,ΨHA1R3domN,
(2.2)

where HA1R3 is a magnetic Sobolev space of order one. A ground state exists since −i∇ − β|x|−1 has a negative energy bound state in L2R3.29 

Unlike previous treatments, here, the arguments remain valid for all values of the parameters α ≥ 0 and β > 0. First, the large B asymptotics of the ground-state energy is derived; the main result is given as Theorem 2.2. Since the pioneering work of Larsen, the model has been considered only in the perturbative regime αβ, and the ground-state energy E0(B) has been approximated as the hydrogenic energy,

EH(B)inf spec HB32βx1,

with a supposedly small correction from the electron–phonon interaction. The large B asymptotics of the hydrogenic energy was derived rigorously in 1981 by Avron et al.3 using ideas from Refs. 9 and 55,

EH(B)=Bβ24lnB2+β2lnBlnlnBβ2γE/2+ln2lnBβ2lnlnB2+2β2γE/21+ln2lnlnB+O(1)asB,
(2.3)

with γE being the Euler–Mascheroni constant, and the expansion can be carried out to arbitrary order. The first three terms are understood heuristically: For large B, the electron is tightly bound in the transverse plane to the lowest Landau orbit while localized in the magnetic-field direction by a one-dimensional effective Coulomb potential that behaves to leading order like a delta well of strength β ln(B/(ln B)2)12,51,52 [see (4.1) and  Appendix B below]. The electron motion is effectively one-dimensional (cf. Ref. 11). The second and third leading-order terms describe the dominant asymptotic behavior of the ground-state energy of this one-dimensional electron confined along the magnetic field. The pronounced anisotropy in the system is reflected by the characteristic length scales of the electron density in the transverse and the magnetic-field direction 1/B and 1/ln B, respectively.

The above hydrogenic heuristics still apply when a coupling to the phonon field is introduced, i.e., α > 0. For large B, the phonons cannot follow the electron’s rapid motion in the transverse plane and so resign themselves to dressing its entire Landau orbit: not only is the electron again localized in the magnetic-field direction by the one-dimensional effective Coulomb potential but also the electron–phonon coupling function is now proportional to the square root of the Fourier transform of the same effective Coulomb interaction [cf. Refs. 34 and 56 and property (k) in Ref. 51]; the system behaves as a one-dimensional strongly coupled polaron to the leading order with interaction strength α ln(B/(ln B)2) confined along the magnetic field by a delta well of strength β ln(B/(ln B)2), i.e., in the effective one-dimensional model, the electron–phonon coupling is mediated by the magnetic field. The analogous large B asymptotics of the polaron hydrogenic energy is derived to second order.

Theorem 2.1.
LetE0(B) be as defined in(2.2)above. Then,
E0(B)=B+e0lnB2+OlnB3/2 as B, with
(2.4)
e0:=infRφ2dxα2Rφ4dxβφ(0)2:Rφ2dx=1
(2.5)
=148α2+6αβ+12β2.
(2.6)

Here, the second leading-order term describes the dominant asymptotic behavior of the ground-state energy of the effective one-dimensional strongly coupled polaron confined along the magnetic field. It is evaluated explicitly by minimizing a nonlinear functional. Furthermore, the cross term in (2.6) indicates that for large B, the effect of the electron–phonon interaction is not perturbative.

The large B asymptotics for the polaron hydrogenic energy is argued differently from the proof of (2.3) given by Avron et al. and generalizes the result of Frank and Geisinger who proved (2.4)–(2.6) when β = 0 using upper and lower bounds to the ground-state energy.17 Their upper bound is established with a trial wave function. Their lower bound is established by showing that the Hamiltonian when restricted to the lowest Landau level is bounded from below in the sense of quadratic forms by an essentially one-dimensional strong-coupling Hamiltonian; the strategy from Ref. 42 is then used to arrive at the nonlinear minimization problem for the second leading-order term along with lower order error terms (cf. Ref. 28).

For proving Theorem 2.1, Frank and Geisinger's strategy in Ref. 17 applies mutatis mutandis. Here, the argument for the upper bound is simplified using the effective Coulomb potential, given in (4.1), which plays an essential role in the proof of (2.3) by Avron et al. but is conspicuously absent from Ref. 17. Here, the argument for the lower bound uses the bathtub principle37 to extract a delta function from the Coulomb potential in the lowest Landau level; just like in Ref. 17, the error terms are bounded by O((lnB)3/2), but—and this has been demonstrated recently in Ref. 24 for a strongly coupled polaron—the error terms in the lower bound should be much smaller. Moreover, the upper bound suggests that the third leading-order term again analogously to its hydrogenic counterpart in (2.3) is 4e0lnBlnlnB.

The expansion in (2.4) should be carried out to higher order. I conjecture the first six leading-order terms behave analogously to their hydrogenic counterparts in (2.3): these should describe the leading asymptotics for the minimization of a nonlinear functional arising naturally in the proof of the upper bound and given below in (4.3), which is a classical approximation to the Fröhlich model of the strongly coupled one-dimensional polaron confined along the magnetic field: in this classical approximation, the electron–phonon interaction is replaced with the one-dimensional effective Coulomb self-interaction of the electron. The second leading-order term above arises from this classical approximation by estimating the one-dimensional effective Coulomb interaction as a delta interaction of strength ln B. Furthermore, in the seventh leading-order term, there should be an order-one quantum correction to the classical approximation (cf. Refs. 24 and 32). These higher-order asymptotics should be provable with better control of the error terms when arguing the lower bound (cf. Ref. 24) and by making full use of the one-dimensional effective Coulomb potential: the electron–phonon interaction of the one-dimensional Hamiltonian that is derived both here and in Frank and Geisinger's proof of the lower bound is described using an artificial coupling function, given in (5.9) below, when instead it should really be argued that the coupling is proportional to the square-root of the Fourier transform of the effective Coulomb potential in (4.1); then, the strategy from Ref. 42 can be used to arrive at the classical approximation in (4.3) now as a lower bound.

The two-term asymptotics of the ground-state energy achieved in Theorem 2.1 suffices for arguing the main result.

Theorem 2.2.
LetΨ(B)Hbe any approximate ground-state wave function satisfyingΨ(B),H(B)Ψ(B)=E0(B)+o(lnB2). The one-dimensional minimization problem in(2.5)admits up to complex phase a unique minimizer,
ϕ0x3=α+2β8αcoshα+2β4x3+tanh12βα+2β,
and, forW, a sum of a bounded Borel measure on the real line and aLRfunction,
limB1lnBRW(x3)R2Ψ(B)F2x,x3lnBdxdx3=RW(x3)ϕ0(x3)2dx3.
(2.7)

By choosing W as a delta-function potential, pointwise convergence is obtained. When α = 0, the limiting density in (2.7) is β/2exp(β|x3|/2) (cf. Refs. 25 and 50).

For proving Theorem 2.2, the strategy is to add to the Hamiltonian ϵ times the one-dimensional potential W scaled appropriately in the magnetic-field direction (cf. Refs. 4, 18, 33, 40, 41, and 44). For

Hϵ(B)H(B)ϵlnB2WlnBx3
(2.8)

and Eϵ(B) being the corresponding ground-state energy, it is argued vis-à-vis Theorem 2.1 that

Eϵ(B)=B+eϵlnB2+OlnB3/2asB,
(2.9)

with

eϵinfφ2=1Rφ2dxα2Rφ4dxβφ(0)2ϵRW(x)φ2dx.
(2.10)

By the variational principle, Eϵ(B)Ψ(B),Hϵ(B)Ψ(B), and the expectation value on the right-hand side evaluates to

E0(B)+o(lnB2)ϵlnBRWx3R2Ψ(B)F2x,x3lnBdxdx3.

Then, for ϵ > 0,

E0(B)Eϵ(B)ϵlnB21lnBRWx3R2Ψ(B)F2x,x3lnBdxdx3+o1,

and taking the limit B, by Theorem 2.1 and (2.9),

e0eϵϵlim supB1lnBRWx3R2Ψ(B)F2x,x3lnBdxdx3.

For ϵ < 0, the above inequality is reversed with “lim sup” replaced by “lim inf.” Hence, the theorem follows if the quotient on the left-hand side has a limit as ϵ → 0. The map ϵeϵ is differentiable at ϵ = 0 for all W if and only if the one-dimensional problem for the energy e0 in (2.5) admits up to the complex phase a unique minimizer: Uniqueness is established by explicitly solving the corresponding Euler–Lagrange equation (cf. Ref. 41) and by the variational principle and a compactness argument,

limϵ0e0eϵϵ=RW(x3)ϕ0(x3)2dx3.
(2.11)

In Sec. III, the differentiation of the one-dimensional energy (2.11) is proved as Theorem 3.3. In Sec. IV, an upper bound to the ground-state energy is proved as Theorem 4.1. In Sec. V, a lower bound to the ground-state energy is proved as Theorem 5.1, and Theorem 2.1 follows from Theorems 4.1 and 5.1. In Sec. VI, the main result Theorem 2.2 is proved.

The one-dimensional problem for the energy e0 in (2.5) is denoted as

e0infφ2=1E0φ,

with

E0φRφ2dxα2Rφ4dxβφ(0)2.
(3.1)

In the following, Lemmas 3.1, 3.2, and Theorem 3.3 apply for all α ≥ 0 and β > 0.

Lemma 3.1.
The minimization problem in(2.5)for the energye0admits up to the complex phase a unique minimizer,
ϕ0x=α+2β8αcoshα+2β4x+tanh12βα+2β
(3.2)
and
e0=148α2+6αβ+12β2.

When α = 0, the corresponding minimizer in (3.2) is β/2exp(β|x|/2).

Proof.
The existence of a minimizer is shown in  Appendix A. Any minimizer up to multiplication by a complex phase is a nonnegative C2R\0 function in H1R solving the Euler–Lagrange equation −ψ″ − αψ3βδ(x)ψ = −λψ with
λ=e0α2Rψ4dx<0.
(3.3)
Or equivalently, it must solve
ψαψ3=λψforx>0
(3.4)
and satisfy the boundary condition
limϵ0+ψϵψ(ϵ)=βψ(0).
(3.5)
The first integral of (3.4) is
ψ2=α2ψ4+λψ2forx>0.
(3.6)
When α = 0, the lemma is obvious. When α > 0, any nonnegative C2R\{0} solution of (3.6) in H1(R) satisfying the boundary condition in (3.5) must be of the form
ψ=2λα1coshλxτ
(3.7)
for some τ (cf. Ref. 16). The boundary condition and that ‖ψ2 = 1 require
tanhτλ=β2λandα4λ=1+tanhτλ,
respectively. Any minimizer up to complex phase must therefore be of the form in (3.7) with
λ=α+2β42andτ=4α+2βtanh12βα+2β.
The explicit calculation of e0 now follows from (3.3).□

Lemma 3.2.

Ifϕnn=1is a minimizing sequence fore0, thenϕnconverges inH1Rto the minimizerϕ0given in(3.2).

Proof.

By Theorem 7.8 in Ref. 37, ϕnn=1 is also a minimizing sequence for e0. It is argued in  Appendix A that every minimizing sequence for e0 has a subsequence converging in H1(R) to some minimizer. Since ϕ0 is up to complex phase the unique minimizer, every subsequence of ϕnn=1 must converge in H1R to ϕ0.□

Theorem 3.3.
LetWbe a sum of a bounded Borel measure on the real line and aLRfunction. Forϵa real parameter, consider the one-dimensional energy
eϵinfφ2=1Eϵφ,
(3.8)
where
EϵφE0φϵRWxφx2dx
(3.9)
with the functionalE0as given in(3.1). Then, the mapϵeϵis differentiable atϵ = 0 and
ddϵϵ=0eϵ=RW(x)ϕ0(x)2dx,
withϕ0being the minimizer given in(3.2)for the energye0.

Proof.
Writing Wμ + ω, where μ is a signed bounded measure on the real line and ωLR, it follows from Hölder’s inequality, the Sobolev inequality, and completion of the square that for φH1R,
Eϵφφ22φ2α2φ22+ϵμR+βϵωφ2234φ22φ22α2φ22+ϵμR+β2ϵωφ22.
(3.10)
Hence, eϵ>, and for each ϵ, there is a ϕϵH1R satisfying
Eϵϕϵeϵ+ϵ2andϕϵ2=1.
(3.11)
By the variational principle,
e0E0ϕϵ=Eϵϕϵ+ϵRWxϕϵx2dxeϵ+ϵ2+ϵRWxϕϵx2dx
and
eϵEϵϕ0=e0ϵRW(x)ϕ0(x)2dx.
Then, for ϵ > 0,
RWxϕϵx2dxϵeϵe0ϵRW(x)ϕ0(x)2dx.
(3.12)
For ϵ < 0, the inequalities in (3.12) are reversed. It suffices therefore to show for any sequence ϵnn=1, ϵn>0 and ϵn → 0 that
limnRWxϕϵnx2dx=RW(x)ϕ0(x)2dx.
(3.13)
Since eϵ>, the concave map ϵeϵ is continuous. It follows from (3.10), the Sobolev inequality, and (3.11) that ϕϵn<C. Thus,
limnϵnRWxϕϵnx2dx=0,
e0=limneϵnlim supnE0ϕϵnϵnRWxϕϵnx2dxϵn2=lim supnE0ϕϵn,
and ϕϵnn=1 is a minimizing sequence for e0. By Lemma 3.2, ϕϵn converges in H1R to ϕ0, and by Theorem 8.6 in Ref. 37, ϕϵn converges also pointwise uniformly to ϕ0 on bounded sets. The convergence in (3.13) now follows.□

Theorem 4.1.
There is a constantC > 0 such that forB > 1,
E0(B)B+e0lnB24e0lnBlnlnB+ClnB.

Theorem 4.1 will follow from Lemmas 4.2 and 4.3.

Lemma 4.2.
E0(B)E0c(B), whereE0c(B):=infPψ:ψ2=1is the classical Pekar energy withPdenoting the three-dimensional magnetic Pekar functional
ψ,HB32ψL2α2R3R3ψx2ψy2xydxdyβR3ψx2xdx.

Proof.
By the variational principle,
E0(B)infΨ,H(B)Ψ:Ψ=1,Ψ=φ(x)ΦandΦF=E0c(B),
with the equality following from the completion of the square in annihilation and creation operators.□

Lemma 4.3.
There is a constantC > 0 such that forB > 1,
E0cBB+e0lnB24e0lnBlnlnB+ClnB.

Proof.
With the lowest Landau state in the zero angular momentum sector,
γBxB2πexpB4x2,i.e.,HBγB=BγB,
by an elementary calculation,12 
VUB(x3)R2γBx2x2+x32dx=0eux32+2uBdu
(4.1)
and
R2R2γBx2γBy2xy2+x3y32dxdy=12VUBx3y32.
(4.2)
For L > 0 and with ϕ0 being the minimizer for the energy e0, μ(B): = ln B − 2 ln ln B and fB(x3):=μ(B)ϕ0(μ(B)x3) by the variational principle, Lemma B.1, and Corollary B.2,
E0c(B)PγBxfBx3=B+RfB2dx3α8R×RfB(x3)2VUBx3y32fB(y3)2dx3dy3βRVUBx3fB(x3)2dx3B+RfB2dx3α/2μ(B)RfB(x3)4dx3βμ(B)fB(0)2+α/2+βL1+8LfB23/2+GB,L/2fB2,
(4.3)
where GB,L:=2lnL+2lnlnB+20euln1u+2BL2+1uduln2. Now choosing L = 1/ln B, it can be verified that GB,L/2<C for B > 1. Since ‖fB′‖2μ(B)‖ϕ0′‖2 and
RfB2dx3α/2μ(B)fB4dx3βμ(B)fB(0)2=μ(B)2e0,
(4.4)
the lemma follows.□

Corollary 4.4.
LetWbe a sum of a bounded Borel measure on the real line and aL(R)function. Forϵa real parameter, Eϵ(B) the ground-state energy of the HamiltonianHϵ(B)in(2.8) and eϵ the one-dimensional energy in(2.10) there is a constantC > 0 such that forB > 1,
Eϵ(B)B+eϵlnB2+ClnBlnlnB+ClnB.

Proof.
By the estimate in (3.10), eϵ>. Then, with the functional Eϵ as given in (3.9) for B > 1, there is a ϕBH1R, ‖ϕB2 = 1, satisfying
EϵϕB<eϵ+1/lnBandϕB2<C.
(4.5)
For L > 0 and with gBx3:=lnBϕBlnBx3 by trivial modifications to Lemma B.1 and Corollary B.2, the arguments in the Proofs of Lemmas 4.2 and 4.3 apply mutatis mutandis, and EϵBEϵc(B) with
Eϵc(B)infψ2=1PψϵlnB2R3WlnBx3ψx,x32dxdx3B+RgB2dx3α/2lnBRgB4dx3βlnBgB(0)2ϵlnB2RWlnBx3gBx32dx3+α/2+βL1+8LgB23/2+G̃B,L/2gB2,
where G̃B,L:=2lnL+20euln1u+2BL2+1uduln2. Now choosing L = 1/ln B, it can be verified that |G̃B,L/2|<2lnlnB+C for B > 1. Since gB2=lnBϕB2, the corollary follows from (4.5) and scaling as in (4.4).□

Theorem 5.1.
There is a constantC > 0 such that forBC,
E0(B)B+e0lnB2ClnB3/2.

The Proof of Theorem 5.1 is given at the end of Subsection V B.

Below in Lemma 5.2, the Coulomb potential restricted to the lowest Landau level is approximated as a delta well of strength ln(B/(ln B)2). First, some remarks about the notation are in order. For ΨHA1R3F, its electron density in the magnetic-field direction is denoted as Ψ̃2x3, i.e.,

Ψ̃x3R2ΨF2x,x3dx1/2.

By Hölder’s inequality, 3Ψ̃23Ψ and Ψ̃H1R. Furthermore, the integral operator P0B acting on L2R2 with kernel

P0B(x,y)B2πeB4xy2eiB2x1y2x2y1
(5.1)

is the projection onto the lowest Landau level, i.e., the ground state of the Landau Hamiltonian HB, and P>B:=1P0B. The operators P0B1 and P0B11 acting on L2R2L2R and L2R2L2RF, respectively, are also denoted P0B.

Lemma 5.2.
LetL > 0. ForB > 1 andΨHA1(R3)F,
R3P0BΨF2x,x3x2+x32dxdx3lnB2lnlnBΨ̃(0)2L1Ψ̃22+8L3Ψ̃23/2Ψ̃21/2+DB,L3Ψ̃2Ψ̃2;
DB,L2lnL+2lnlnB+22BL2+1+1+2ln1+2BL2+1ln2.
(5.2)

Proof.
By Hölder’s inequality,
P0BΨF2x,x3B2πΨ̃x32,
(5.3)
and since P0B is a projection,
R3P0BΨF2x,x3dxdx3RΨ̃x32dx3.
(5.4)
It follows from the bathtub principle37,10 that the maximum of the expression
R2Gx,x3x2+x3dx
over all functions G satisfying the conditions (5.3) and (5.4) above is attained by the function
Gmaxx,x3=B2πΨ̃x32when xR,0when x>R,
where R=2/B. Therefore,
RR2P0BΨF2x,x3x2+x32dxdx3Rx2/BGmaxx,x3x2+x32dxdx3=RVLB(x3)Ψ̃x32dx3,
with
VLB(x3)22B+x32+x3.
(5.5)
The lemma follows from Corollary B.3 in  Appendix B.□

Lemma 5.2 is used to prove the following corollary.

Corollary 5.3.
Letτ > 0. There is a constantC > 0 such that forBCandΨHA1R3F,
3Ψ2τΨ,P0Bx1P0BΨτ24lnB2+τ2lnBlnlnBClnBΨ2.

Proof.
Recall that for ϵ > 0, 3Ψ̃2Ψ̃2ϵ3Ψ̃22+ϵ1Ψ̃22. With DB,L as given in (5.2), under the assumption that 2τϵDB,L<1/2 and denoting μ(B)=lnB2lnlnB, it follows from Lemma 5.2 that
3Ψ2τΨ,P0Bx1P0BΨ12τϵDB,L3Ψ̃22τμ(B)Ψ̃(0)2τL1+ϵ1DB,LΨ̃22+τϵDB,L3Ψ̃228τLΨ̃21/23Ψ̃23/2τ2μ(B)2412τϵDB,LΨ̃22τ432L2DB,L3ϵ3+L1+ϵ1DB,LΨ̃22
τ24τ3ϵDB,Lμ(B)2432τL2DB,L3ϵ3τLτDB,LϵΨ̃22.
At the second inequality, 3Ψ̃22ζ(Ψ̃(0))2(1/4)ζ2Ψ̃22 for ζ ≥ 0 and ax2bx3/2 ≥−(27/256)b4/a3 for a > 0, b ≥ 0, and x ≥ 0 are used. Now choosing ϵ = 1/ln B and L = 1/ln B, the above assumption that 2τϵDB,L<1/2 can be verified to be true for BC. The lemma follows.□

Corollary 5.3 is used to prove Lemma 5.6 in Subsection V B. The arguments for Lemma 5.2 and Corollary 5.3 are used to prove Lemma 5.12 in Subsection V C.

Below in Lemma 5.8, a Hamiltonian restricted to the lowest Landau level is derived as a lower bound. First, an ultraviolet cutoff is needed. With K>0 and ΓK:=kR3:maxk,k3K, the cutoff Hamiltonian is denoted as

hKco18απKHB32+12ΓKakakdk+12N+α2πΓKakeikx|k|+akeikx|k|dk.

Lemma 5.4.

IfK>8α/π, thenH(B)hKcoβx11/4.

Proof.

The lemma follows from Lemma 5.1 in Ref. 17.□

Lemma 5.4 is used to prove the following lemma.

Lemma 5.5.
There exists a constantC > 0 such that forBC,
H(B)18αlnB2πBHB+128αlnB2πB32βx1+12NClnB2.

Proof.
Let K=B/lnB2 and K3=16αlnB/π. By completion of the square in annihilation and creation operators,
|k3|K3|k|K12akak+α2πak|k|eikx+α2πak|k|eikxdkdk3α2π2|k3|K3|k|K1|k|2dkdk3=απ0K3lnK2+k32k32dk3=απK3lnK2K32+1+2KarctanK3KαπK3lnK2K32+1+2+23K32K2ClnB2,
(5.6)
valid for some constant C and BC. The argument corrects a mistake in the Proof of Lemma 5.2 in Ref. 17.
Denoting ΛB=k,k3ΓK:K3k3KandkK, it can be argued as in the Proof of Lemma 5.2 in Ref. 17 that for large B,
α2πΛBakkeikx+akkeikxdk123212ΛBakakdk+12.
(5.7)
The lemma follows from Lemma 5.4 and the estimates (5.6) and (5.7).□

Lemma 5.5 and Corollary 5.3 from Subsection V A are used to prove the following lemma.

Lemma 5.6.
There exists a constantC > 0 such that forBCand allΨHA1R3domNwith ‖Ψ‖ = 1,
Ψ,H(B)ΨB+B2P>BΨ2+12Ψ,NΨClnB2.

Proof.
Recall that by the diamagnetic inequality for τ > 0,
Ψ,HB32ΨτΨ,β|x|1Ψ41β2τ2Ψ2.
Then, with 0 < η < 1, A>1, and θ:=128αlnB2πB for B large,
θΨ,HB32ΨΨ,β|x|1Ψ=θP0BΨ,HB32P0BΨP0BΨ,β|x|1P0BΨ+θ1ηP>BΨ,HB32P>BΨ+θηP>BΨ,HB32P>BΨP>BΨ,β|x|1P>BΨP0BΨ,β|x|1P>BΨP>BΨ,β|x|1P0BΨθP0BΨ,HB32P0BΨP0BΨ,β|x|1P0BΨ+θ1ηP>BΨ,HB32P>BΨA2β24θη1P>BΨ2+A1P>BΨ,β|x|1P>BΨ2P0BΨ,β|x|1P0BΨ12P>BΨ,β|x|1P>BΨ12θP0BΨ,HBP0BΨ+θP0BΨ,32P0BΨA/(A1)P0BΨ,β|x|1P0BΨ+θ1ηP>BΨ,HB32P>BΨA2β24θη1P>BΨ2θB+θP0BΨ,32P0BΨA/(A1)P0BΨ,β|x|1P0BΨ+θBP>BΨ2+θB3BηA2β24θη1P>BΨ2.
Choosing η = 1/4 and A=Bθ/2β for B large A/(A1)<2 and by Corollary 5.3, there exists some C > 0 such that for BC,
θP0BΨ,32P0BΨA/(A1)P0BΨ,β|x|1P0BΨClnB2.
The lemma now follows from Lemma 5.5.□

The following observation is immediate from Theorem 4.1: For every M>e0 and large B, there exist wave functions ΨHA1(R3)domN satisfying

Ψ,H(B)ΨB+MlnB2andΨ=1.
(5.8)

The estimate in (5.8) and Lemma 5.6 are used to prove the following corollary.

Corollary 5.7.
For everyMR, there exists a constantCM > 0 such that forBCMand allΨHA1(R3)domNsatisfying (5.8),
P>BΨ2CMlnB2B1andΨ,NΨCMlnB2.

Proof.

The corollary follows from Lemma 5.6.□

Corollary 5.7 and Lemma 5.4 are used to prove Lemma 5.8.

Lemma 5.8.
LetK>8α/πand1<A<B/lnB. Denotingκ=18α/πKfor everyMR, there exists a constantCM > 0 such that forBCMand allΨHA1(R3)domNsatisfying(5.8),
Ψ,H(B)ΨP0BΨ,hKcoβ1/1A1x1P0BΨ+κBP>BΨ2CMlnB2KB1+KB1CMκ1lnB1/4.

Proof.
It can be argued that as in the Proof of Lemma 5.6 with 0 < η < 1,
κΨ,HB32ΨΨ,β|x|1ΨκP0BΨ,HB32P0BΨA/A1P0BΨ,β|x|1P0BΨ+κBP>BΨ2+κ2B3BηA2β24ηκ1P>BΨ2.
It now follows from Lemma 5.4 that
Ψ,H(B)ΨP0BΨ,hKcoβ1/1A1x1P0BΨ+κBP>BΨ2+κ2B3BηA2β24ηκ1P>BΨ214+P>BΨ,ΓKakak+α2πakeikx|k|+α2πakeikx|k|dkP>BΨ+P0BΨ,α2πΓKakeikx|k|+akeikx|k|dkP>BΨ+P>BΨ,α2πΓKakeikx|k|+akeikx|k|dkP0BΨ.
By completion of the square in annihilation and creation operators,
ΓKakak+α2πak|k|eikx+α2πak|k|eikxdkα4π2ΓK1|k|2dkdk3=α2ln(2)+π4πK,
and by Corollary 5.7 for BCM,
P>BΨ,ΓKakak+α2πakeikx|k|+α2πakeikx|k|dkP>BΨα2ln(2)+π4πKP>BΨ2CMKB1lnB2.
Furthermore, it can be argued as in the Proof of Lemma 5.4 in Ref. 17 and using Corollary 5.7 that for BCM,
P0BΨ,α2πΓKak|k|eikxdkP>BΨ
CKP>BΨN+1P0BΨCMKB1lnB2.
The remaining interaction terms are estimated similarly. Choosing η = 2/3, the lemma follows from Corollary 5.7.□

The Hamiltonian P0BhKcoβ1/1A1x1P0B realized as a lower bound in Lemma 5.8 is estimated from below.

Proposition 5.9.
There exists a constantC > 0 such that forB, κ, andKsatisfyingBC,ClnB1/2κC1lnB,KB, and some1<A<B/lnB,
P0BhKcoβ1/1A1x1P0BκB+κ1lnB2e0Cκ1/2lnB3/2C1+κ2lnBP0B.

Proposition 5.9 and Lemma 5.8 are used to prove Theorem 5.1.

Proof of Theorem 5.1.
Fix M>e0. For large B by Theorem 4.1, there exist wave functions satisfying (5.8). It suffices to argue the desired lower bound on Ψ,H(B)Ψ with those wave functions. By Lemma 5.8 and Proposition 5.9,
Ψ,H(B)ΨκB+κ1lnB2e0Cκ1/2lnB3/2C1+κ2lnBP0BΨ2ClnB2KB1+KB1Cκ1lnBC,
with κ=18α/πK. Choosing K=BlnB4/3 and since P0BΨ1,
Ψ,H(B)ΨB+e0lnB2ClnB3/2,
yielding the claimed lower bound.□

Proposition 5.9 is proved in Subsection V C.

1. Reduction to one dimension

Below in Lemma 5.10, a one-dimensional Hamiltonian is derived as a lower bound. In Ref. 17, the authors consider a one-dimensional Hamiltonian with 0<K3K and 1KK,

h1dκ132+|k3|K3âk3âk3dk3+α2π|k3|K3νk3âk3eik3x3+âk3eik3x3dk3,

acting on L2(R3)F with κ1:=κ8α/πK301+t1exptK32/2Bdt and κ as in the statement of Proposition 5.9, and

âk31νk31kKakkeikxdk,
ν(k3)1kKk2dk12=πlnK2+k32ln1+k3212,
(5.9)

satisfying [âk3,âk3]=δk3k3 and âk3,âk3=[âk3,âk3]=0 for k3,k3R.

Lemma 5.10.
Denotingκ2=κ2απ1K3K2,
P0BhKcoβ1/1A1x1P0Bκ2BP0B+P0Bh1dβ1/1A1P0Bx1P0BP0B1+α2P0B.

Proof.

The lemma follows from Lemmas 6.1–6.4 in Ref. 17. The Proof of Lemma 6.3 in Ref. 17 uses an incorrect vector operator for cutting off high modes in the k-direction (cf. Ref. 43); the mistake can be fixed, and the lemma stands true.□

2. Localization and decomposition

Below in Lemma 5.11, the one-dimensional Hamiltonian is bounded from below by a Hamiltonian describing an electron interacting with only finitely many modes. In Ref. 17, the authors decompose the mode space into M intervals, indexed with b, each of length 2K3/M, and consider for uR and 0 < γ < 1 the block Hamiltonian

hγ(u)κ132+b1γAb(u)*Ab(u)+α2πV(b)Ab(u)eikbx3+Ab(u)*eikbx3,

with kb being a mode in block b and the block creation and annihilation operators Ab(u)* and Ab(u) acting on FL2R3, where

Ab(u)1V(b)bνk3eik3kbuâk3dk3,

Vb:=bνk32dk312, satisfying

[Ab(u),Ab(u)*]=δbb,[Ab(u),Ab(u)]=[Ab(u)*,Ab(u)*]=0for all blocksb,b.

Furthermore, localizing the electron in the x3-direction on an interval of length J, the following lemma is proved.

Lemma 5.11.
ForχC0R,χ2=1being a nonnegative function supported on the interval [−1/2, 1/2] and forJ > 0, denotingχuJ(x3)=J1/2χJ1x3u,
h1dβ1/1A1P0Bx1P0BRχuJhγ(u)β1/1A1P0Bx1P0BχuJduαK32J24π2γM2Rχ22J2,
with
R|k3|K3νk32dk3=π|k3|K3lnK+k32ln1+k32dk3.

Proof.

The lemma follows from Lemma 6.5 in Ref. 17.□

3. Error estimates

The block Hamiltonian is bounded from below. Similarly, as in Refs. 17, 28, and 42 representing the block creation and annihilation operators by coherent state integrals and completing the square, it follows for a suitably chosen kb that

hγ(u)β1/1A1P0Bx1P0BIM,
(5.10)

where

Iinfϕ2=1κ13ϕ22α4π21γRνk32R3eik3x3ϕx,x32dx2dk3β1A1R3x1P0Bϕx,x32dx

will be estimated from below in Lemma 5.12. Combining (5.10) with Lemma 5.11,

h1dβ1/1A1P0Bx1P0BIMαK32J24π2γM2Rχ22J2.

Now, it follows from Lemma 5.10 that for some constant C > 0,

P0BhKcoβ1/1A1x1P0BκBP0B+IP0BCK3BK2+M+K32J2γM2R+1J2P0B.
(5.11)

The argument for Lemma 5.2 from Subsection V A is used to prove Lemma 5.12.

Lemma 5.12.
For anyL > 0 andϵ > 0 and withD(B,L)as given in(5.2)assuming
2lnK1γ=μ(B)1A1andκ̃1κ14βϵD(B,L)lnKμ(B)1γ>0
(5.12)
withμ(B): = ln B − 2 ln ln B,
I4lnK2e0κ̃11γ22βlnKμ(B)1γ432L2D(B,L)3ϵ3+1L+D(B,L)ϵ.

Proof.
For ϕL2R3, ‖ϕ2 = 1, and ϕ̃(x3):=R2ϕx,x32dx1/2,
κ13ϕ22α4π21γRνk32R3eik3x3ϕx,x32dx2dk3β1A1R3P0Bϕx,x32xdxκ13ϕ22αlnK1γRR2ϕx,x32dx2dx3β1A1R3P0Bϕx,x32xdxκ12βϵD(B,L)1A13ϕ̃22αlnK1γϕ̃44βμ(B)1A1ϕ̃(0)2β1A1432L2D(B,L)3ϵ3+1L+D(B,L)ϵ4lnK2e0κ̃11γ22βlnKμ(B)1γ432L2D(B,L)3ϵ3+1L+D(B,L)ϵ.
At the first inequality,
νk32=πlnK2+k32ln1+k322πlnK
(5.13)
is used along with Plancherel’s identity. At the second inequality, the bathtub principle in Lemma 5.2 and the argument for Corollary 5.3 apply mutatis mutandis. At the third inequality, the assumptions in (5.12) are used.□

From Lemma 5.12 and further assuming that

κκ̃1κ2andγ12,
(5.14)

it can be seen that there is a constant C > 0 such that

I4κ1lnK2e0ClnK2κκκ̃1κ+γ+lnKμ(B)L2DB,L3ϵ3+1L+DB,Lϵ.

With the above bound and (5.11), the argument in Ref. 17 applies mutatis mutandis, and choosing J2=κ1/5K33/5lnK3/5, M=[J2], and γ=κ4/5K33/5lnK7/5 yields

P0BhKcoβ1/1A1x1P0BκB+4lnK2κe0P0BCκ1κκ1lnK2+κ1/5K33/5lnK3/5+BK3K2P0BClnKμ(B)L2ϵ3DB,L3+1L+DB,Lϵ+ϵDB,LlnK2κ2P0B.

It is shown in Ref. 17 choosing K=B1/2 and K3=κ1/2lnB3/2 that

κ1κκ1lnK2+κ1/5K33/5lnK3/5+BK3K2Cκ1/2lnB3/2.

Now choosing L = 1/ln B and ϵ = 1/ln B, since 1DB,LC for BC,

lnKμ(B)L2ϵ3DB,L3+1L+DB,Lϵ+ϵDB,LlnK2κ2C1+κ2lnB.

With the above choice of parameters and the assumptions ClnB1/2κC1lnB and KB, the conditions in (5.12) and (5.14) and that 0<K3K, 1KK and 1<A<B/lnB are verified, and

P0BhKcoβ1/1A1x1P0BκB+κ1lnB2e0Cκ1/2lnB3/2C1+κ2lnBP0B,

thereby concluding the Proof of Proposition 5.9.

Corollary 5.13.
LetWbe a sum of a bounded Borel measure on the real line and aL(R)function. Forϵa real parameter andEϵ(B) andeϵas in Corollary 4.4, there is a constantC > 0 such that forBC,
Eϵ(B)B+eϵlnB2ClnB3/2.

Proof.

The above arguments apply mutatis mutandis (see also Ref. 31).□

Proof of Theorem 2.2.
Let Eϵ(B) be the ground-state energy of the Hamiltonian Hϵ(B) in (2.8) and let the one-dimensional energy eϵ be as defined in (2.10). The large B asymptotics of Eϵ(B) in (2.9) follows from Corollaries 4.4 and 5.13. As explained in the Introduction, it follows from the variational principle, Theorem 2.1, and (2.9) that for ϵ > 0,
e0eϵϵlim supB1lnBRWx3R2Ψ(B)F2x,x3lnBdxdx3,
and for ϵ < 0,
e0eϵϵlim infB1lnBRWx3R2Ψ(B)F2x,x3lnBdxdx3.
Theorem 2.2 now follows from Theorem 3.3.□

This paper is based on the work supported by the NSF (Grant No. DMS-1600560). I thank Michael Loss for several valuable suggestions, Rupert Frank for the general strategy used in the Proof of Theorem 2.2, and the referee for helpful comments about the presentation.

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

In the following, Theorem A.1 applies for all α ≥ 0 and β > 0.

Theorem A.1.

If a sequenceψnn=1,ψn2=1satisfieslimnE0ψn=e0with the functionalE0as given in(3.1), then there exists a subsequence{ψnk}k=1and someψH1Rsuch thatψ2=1,E0ψ=e0, andψnkψH10ask.

Proof.
For φH1R, it can be argued as in the Proof of Theorem 3.3 that
E0φ34φ22φ22α2φ22+β2.
Furthermore, since ψn is a minimizing sequence, E0ψn<e0+1 for n large. Then, ψnH1<C and there exists a subsequence {ψnk}k=1 and some ψH1R such that ψnk converges to ψ weakly in H1R.
Step 1 (compactness): it is argued that the subsequence {ψnk}k=1 satisfies
δ>0,R<s.t.ψnkL2x<R2>1δ.
(A1)
Essential to the argument is the binding inequality e0<eT, where eT:=infφ2=1ETφ and ETφ:=Rφ2dx(α/2)Rφ4dx is the translation-invariant problem admitting when α > 0 a symmetric decreasing minimizer ϕT.17 When α = 0, the binding inequality is obvious, and when α > 0,
e0E0ϕT=ETϕTβϕT2(0)=eTβϕT2(0)<eT.
Moreover, it should be noted that
E0φe0φ22andETφeTφ22whenφ21.
(A2)
Also a quadratic partition of unity is chosen, χ2+χ̃21, where 0 ≤ χ ≤ 1 is a smooth function with χ(x) = 1 when |x| < 1/2 and χ(x) = 0 when |x| > 1. Denoting χR=χR1, it follows from (A2) that
E0ψnk=E0χRψnk+ETχ̃RψnkαRχR2χ̃R2ψnk4dxRψnk2χR2+χ̃R2dxe0eTχRψnk22+eTαRχR2χ̃R2ψnk4dxRψnk2χR2+χ̃R2dx.
(A3)
Since χ,χ̃ have bounded derivatives,
Rψnk2χR2+χ̃R2dx<CR2
for some C > 0. Furthermore, with DR:=R/2|x|R,
RχR2χ̃R2ψnk4dxψnkL4DR4andψnkL4DR4ψL4DR4
by Rellich-Kondrashov, so the first term in (A3) can also be made arbitrarily small with R chosen to be large enough uniformly in k. Hence, for any δ > 0, there is some R such that for all k,
E0ψnke0eTχRψnk22+eTδ(eTe0)/2.
(A4)
Since {ψnk}k=1 is a minimizing sequence for e0, for large k,
E0ψnke0+δ(eTe0)/2.
(A5)
Compactness now follows from (A4) and (A5).
Step 2 (weak limit is a minimizer): by Rellich–Kondrashov and the compactness property in (A1),
ψnkψ20andψ2=1.
(A6)
Since ψnH1<C, by Sobolev and Hölder’s inequalities,
R|ψnk|4ψ4dxCψnkψ20.
Furthermore, by Theorem 8.6 in Ref. 37, ψnk(0)ψ0, so
α2R|ψnk|4dx+β|ψnk0|2α2Rψ4dx+βψ02.
(A7)
Then, since lim infkψnk2ψ2, e0=limkE0(ψnk)E0ψe0 and E0ψ=e0.
Step 3 (convergence in H1(R)): from (A7),
limkψnk22=limkE0(ψnk)+α2ψnk44+βψnk(0)2=e0+α2ψ44+βψ02=ψ22,
and since ψnkψ in H1, ψnkψ20. Strong convergence in H1 now follows from (A6).□

Recall the effective Coulomb potential VUB from (4.1).

Lemma B.1.
For anyL > 0 andϕH1(R), one has forB > 1,
RVUBxϕx2dxlnB2lnlnBϕ(0)2L1ϕ22+8Lϕ23/2ϕ21/2+GB,Lϕ2ϕ2;
GB,L2lnL+2lnlnB+20euln1u+2BL2+1uduln2.

Proof.
Writing RVUBxϕx2 as
|ϕ(0)|2|x|LVUB(x)+|x|LVUB(x)|ϕ(x)|2|ϕ(0)|2+|x|LVUB(x)|ϕ(x)|2,
it is possible to bound
|x|LVUBxϕ(x)2dxL1|x|L|ϕ(x)|2dx,
(B1)
|x|LVUBxϕ(x)2ϕ(0)2dx8Lϕ23/2ϕ21/2
(B2)
and to evaluate the integral
|x|LVUBxdx=lnB2lnlnB+GB,L.
To arrive at the bound in (B2), the following inequalities are used:
ϕ(x)ϕ(0)xϕ2andϕ2ϕ2ϕ2.
(B3)
The lemma now follows from (B1), (B2), and the rightmost inequality of (B3).□

Corollary B.2.
For anyL > 0 andϕH1(R), one has forB > 1,
R×Rϕ(x)212VUBxy2ϕy2dxdylnB2lnlnBϕ44L1ϕ24+8Lϕ23/2ϕ25/2+GB,L/2ϕ2ϕ23
withGB,Las above.

Proof.

The corollary follows from Lemma B.1.□

Now, recall the potential VLB from (5.5).

Corollary B.3.
For anyL > 0 andϕH1(R), one has for B > 1,
RVLB(x)ϕx2dxlnB2lnlnBϕ02L1ϕ22+8Lϕ23/2ϕ21/2+D(B,L)ϕ2ϕ2;
DB,L2lnL+2lnlnB+22BL2+1+1+2ln1+2BL2+1ln2.

Proof.
Evaluating the integral
x<LVLB(x)dx=lnB2lnlnB+DB,L,
the argument follows the Proof of Lemma B.1.□

1.
I.
Anapolitanos
and
M.
Hott
, “
Asymptotic behavior of the ground state energy of a Fermionic Fröhlich multipolaron in the strong coupling limit
,” arXiv:1601.05272 (
2016
).
2.
I.
Anapolitanos
and
B.
Landon
, “
The ground state energy of the multipolaron in the strong coupling limit
,”
Lett. Math. Phys.
103
(
12
),
1347
1366
(
2013
).
3.
J. E.
Avron
,
I. W.
Herbst
, and
B.
Simon
, “
Schrödinger operators with magnetic fields. III. Atoms in homogeneous magnetic field
,”
Commun. Math. Phys.
79
,
529
572
(
1981
).
4.
B.
Baumgartner
, “
The Thomas-Fermi-Theory as result of a strong-coupling-limit
,”
Commun. Math. Phys.
47
(
3
),
215
219
(
1976
).
5.
B.
Baumgartner
,
J. Ph.
Solovej
, and
J.
Yngvason
, “
Atoms in strong magnetic fields: The high field limit at fixed nuclear charge
,”
Commun. Math. Phys.
212
(
3
),
703
724
(
2000
).
6.
R. D.
Benguria
and
G. A.
Bley
, “
Improved results on the no-binding of bipolarons
,”
J. Phys. A: Math. Theor.
45
(
4
),
045205
(
2012
).
7.
G. A.
Bley
, “
Absence of binding in the Nelson and piezoelectric polaron models
,”
Rev. Math. Phys.
31
(
2
),
1950006
(
2019
).
8.
G. A.
Bley
and
L. E.
Thomas
, “
Estimates on functional integrals of quantum mechanics and non-relativistic quantum field theory
,”
Commun. Math. Phys.
350
(
1
),
79
103
(
2017
).
9.
R.
Blankenbecler
,
M. L.
Goldberger
, and
B.
Simon
, “
The bound states of weakly coupled long-range one-dimensional quantum Hamiltonians
,”
Ann. Phys.
108
(
1
),
69
78
(
1977
).
10.
F.
Bonetto
and
M.
Loss
, “
Entropic chaoticity for the steady state of a current carrying system
,”
J. Math. Phys.
54
,
103303
(
2013
).
11.
R.
Brummelhuis
and
P.
Duclos
, “
Effective Hamiltonians for atoms in very strong magnetic fields
,”
J. Math. Phys.
47
,
032103
(
2006
).
12.
R.
Brummelhuis
and
M. B.
Ruskai
, “
One-dimensional models for atoms in strong magnetic fields, II: Anti-symmetry in the Landau levels
,”
J. Stat. Phys.
116
(
1
),
547
570
(
2004
).
13.
M. D.
Donsker
and
S. R. S.
Varadhan
, “
Asymptotics for the polaron
,”
Commun. Pure Appl. Math.
36
(
4
),
505
528
(
1983
).
14.
W.
Dybalski
and
H.
Spohn
, “
Effective mass of the polaron–revisited
,”
Ann. Henri Poincaré
21
,
1573
1594
(
2020
).
15.
R.
Evrard
,
E.
Kartheuser
, and
J.
Devreese
, “
Theory of cyclotron resonance of Fröhlich polarons in high magnetic fields
,”
Phys. Status Solidi B
41
(
1
),
431
438
(
1970
).
16.
R. L.
Frank
,
Ground States of Semi-Linear PDEs
, Lecture Notes from the “Summer School on Current Topics in Mathematical Physics,” preprint; available at http://www.math.caltech.edu/~rlfrank/luminy140202.pdf.
17.
R. L.
Frank
and
L.
Geisinger
, “
The ground state energy of a polaron in a strong magnetic field
,”
Commun. Math. Phys.
338
(
1
),
1
29
(
2015
).
18.
R. L.
Frank
,
K.
Merz
,
H.
Siedentop
, and
B.
Simon
, “
Proof of the strong Scott conjecture for Chandrasekhar atoms
,” arXiv:1907.04894v2 (
2019
).
19.
R. L.
Frank
,
E. H.
Lieb
, and
R.
Seiringer
, “
Binding of polarons and atoms at threshold
,”
Commun. Math. Phys.
313
,
405
424
(
2012
).
20.
R. L.
Frank
,
E. H.
Lieb
,
R.
Seiringer
, and
L. E.
Thomas
, “
Bipolaron and N-polaron binding energies
,”
Phys. Rev. Lett.
104
,
210402
(
2010
).
21.
R. L.
Frank
,
E. H.
Lieb
,
R.
Seiringer
, and
L. E.
Thomas
, “
Stability and absence of binding of multi-polaron systems
,”
Publ. Math. L’IHÉS
113
(
1
),
39
67
(
2011
).
22.
R. L.
Frank
,
E. H.
Lieb
,
R.
Seiringer
, and
L. E.
Thomas
, “
Binding, stability, and non-binding of multi-polaron systems
,” in
Proceedings of the QMath11 Conference
, edited by
P.
Exner
(
World Scientific
,
Singapore
,
2011
), pp.
21
32
.
23.
R. L.
Frank
,
E. H.
Lieb
,
R.
Seiringer
, and
L. E.
Thomas
, “
Ground state properties of multi-polaron systems
,” in
XVIIth International Congress on Mathematical Physics, Proceedings of the ICMP Held in Aalborg, August 6–11, 2012
, edited by
A.
Jensen
(
World Scientific
,
Singapore
,
2013
), pp.
477
485
.
24.
R. L.
Frank
and
R.
Seiringer
, “
Quantum corrections to the Pekar asymptotics of a strongly coupled polaron
,” arXiv:1902.02489v2 (
2019
).
25.
R.
Froese
and
R.
Waxler
, “
The spectrum of a hydrogen atom in an intense magnetic field
,”
Rev. Math. Phys.
6
,
699
832
(
1995
).
26.
H.
Fröhlich
,
H.
Pelzer
, and
S.
Zienau
, “
Properties of slow electrons in polar materials
,”
London, Edinburgh, Dublin Philos. Mag. J. Sci.
41
,
221
242
(
1950
).
27.
B.
Gerlach
and
H.
Löwen
, “
Analytical properties of polaron systems or: Do polaronic phase transitions exist or not?
,”
Rev. Mod. Phys.
63
(
1
),
63
90
(
1991
).
28.
R.
Ghanta
, “
Exact ground state energy of the 1D strong-coupling polaron
,” Junior thesis,
Princeton
,
2012
.
29.
M.
Griesemer
,
E. H.
Lieb
, and
M.
Loss
, “
Ground states in non-relativistic quantum electrodynamics
,”
Invent. Math.
145
(
3
),
557
595
(
2001
).
30.
M.
Griesemer
and
J. S.
Møller
, “
Bounds on the minimal energy of translation invariant N-polaron systems
,”
Commun. Math. Phys.
297
(
1
),
283
297
(
2010
).
31.
M.
Griesemer
and
D.
Wellig
, “
The strong-coupling polaron in static electric and magnetic fields
,”
J. Phys. A: Math. Theor.
46
,
425202
(
2013
).
32.
E. P.
Gross
, “
Strong coupling polaron theory and translational invariance
,”
Ann. Phys.
99
(
1
),
1
29
(
1976
).
33.
A.
Iantchenko
,
E. H.
Lieb
, and
H.
Siedentop
, “
Proof of a conjecture about atomic and molecular cores related to Scott’s conjecture
,”
J. Reine Angew. Math.
472
,
177
195
(
1996
).
34.
E. A.
Kochetov
,
H.
Leschke
, and
M. A.
Smondyrev
, “
Diagrammatic weak-coupling expansion for the magneto-polaron energy
,”
Z. Phys. B
89
,
177
186
(
1992
).
35.
D. M.
Larsen
, “
Shallow donor levels of InSb in a magnetic field
,”
J. Phys. Chem. Solids
29
,
271
280
(
1968
).
36.
T. D.
Lee
,
F. E.
Low
, and
D.
Pines
, “
The motion of slow electrons in a polaron crystal
,”
Phys. Rev.
90
(
2
),
297
302
(
1953
).
37.
E. H.
Lieb
and
M.
Loss
,
Analysis
, Graduate Studies in Mathematics Vol. 14, 2nd ed. (
American Mathematical Society
,
Providence, RI
,
2001
).
38.
E. H.
Lieb
and
R.
Seiringer
, “
Equivalence of two definitions of the effective mass of a polaron
,”
J. Stat. Phys.
154
,
51
57
(
2014
).
39.
E. H.
Lieb
and
R.
Seiringer
, “
Divergence of the effective mass of a polaron in the strong coupling limit
,”
J. Stat. Phys.
180
,
22
33
(
2020
).
40.
E. H.
Lieb
and
B.
Simon
, “
The Thomas-Fermi theory of atoms, molecules and solids
,”
Adv. Math.
23
,
22
116
(
1977
).
41.
E. H.
Lieb
,
J. P.
Solovej
, and
J.
Yngvason
, “
Asymptotics of heavy atoms in high magnetic fields: I. Lowest Landau band regions
,”
Commun. Pure Appl. Math.
47
,
513
591
(
1994
).
42.
E. H.
Lieb
and
L. E.
Thomas
, “
Exact ground state energy of the strong-coupling polaron
,”
Commun. Math. Phys.
183
(
3
),
511
519
(
1997
).
43.
E. H.
Lieb
and
K.
Yamazaki
, “
Ground-state energy and the effective mass of the polaron
,”
Phys. Rev.
111
(
3
),
728
733
(
1958
).
44.
E. H.
Lieb
and
H.-T.
Yau
, “
The Chandrasekhar theory of stellar collapse as the limit of quantum mechanics
,”
Commun. Math. Phys.
112
(
1
),
147
174
(
1987
).
45.
H.
Löwen
, “
Analytical behavior of the ground-state energy and pinning transitions for a bound polaron
,”
J. Math. Phys.
29
(
6
),
1505
1513
(
1988
).
46.
H.
Löwen
, “
Spectral properties of an optical polaron in a magnetic field
,”
J. Math. Phys.
29
(
6
),
1498
1504
(
1988
).
47.
T.
Miyao
and
H.
Spohn
, “
The bipolaron in the strong coupling limit
,”
Ann. Henri Poincaré
8
,
1333
1370
(
2007
).
48.
S. I.
Pekar
, “
Research in electron theory of crystals
,” Technical Report No. AEC-tr-5575,
United States Atomic Energy Commission
,
Washington, DC
,
1963
.
49.
P. M.
Platzman
, “
Ground-state energy of bound polarons
,”
Phys. Rev.
125
(
6
),
1961
1965
(
1962
).
50.
A. R. P.
Rau
and
L.
Spruch
, “
Energy levels of hydrogen in magnetic fields of arbitrary strength
,”
Astrophys. J.
207
,
671
679
(
1976
).
51.
M. B.
Ruskai
and
E.
Werner
, “
Study of a class of regularizations of 1/x using Gaussian integrals
,”
SIAM J. Math. Anal.
32
(
2
),
435
463
(
2000
).
52.
L. I.
Schiff
and
H.
Snyder
, “
Theory of the quadratic Zeeman effect
,”
Phys. Rev.
55
,
59
63
(
1939
).
53.
R.
Seiringer
, “
The polaron at strong coupling
,”
Rev. Math. Phys.
33
,
2060012
(
2020
).
54.
J. M.
Shi
,
F. M.
Peeters
, and
J. T.
Devreese
, “
Magnetopolaron effect on shallow donor states in GaAs
,”
Phys. Rev. B
48
(
8
),
5202
(
1993
).
55.
B.
Simon
, “
The bound state of weakly coupled Schrödinger operators in one and two dimensions
,”
Ann. Phys.
97
(
2
),
279
288
(
1976
).
56.
M. A.
Smondyrev
,
E. A.
Kochetov
,
G.
Verbist
,
F. M.
Peeters
, and
J. T.
Devreese
, “
Equivalence of 3D bipolarons in a strong magnetic field to 1D bipolarons
,”
Europhys. Lett.
19
(
6
),
519
524
(
1992
).
57.
H.
Spohn
, “
Effective mass of the polaron: A functional integral approach
,”
Ann. Phys.
175
(
2
),
278
318
(
1987
).
58.
I.
Zorkani
,
R.
Belhissi
, and
E.
Kartheuser
, “
The ground state energy of a bound polaron in the presence of a magnetic field
,”
Phys. Status Solidi B
197
,
411
419
(
1996
).