Effective mass of the polaron -- revisited

TL;DR

This paper investigates the energy-momentum relation of the Fröhlich polaron, establishing conditions under which the inverse effective mass is positive and coincides with the diffusion constant, with implications for polaron models across coupling regimes.

Contribution

It combines spectral theory and the central limit theorem to analyze the effective mass of the polaron, extending results to models with ultraviolet cut-off and all coupling constants.

Findings

Inverse effective mass is positive outside an intermediate coupling range.

Inverse effective mass equals the diffusion constant.

Results apply to polaron models with ultraviolet cut-off.

Abstract

Properties of the energy-momentum relation for the Fr\"ohlich polaron are of continuing interest, especially for large values of the coupling constant. By combining spectral theory with the available results on the central limit theorem for the polaron path measure we prove that, except for an intermediate range of couplings, the inverse effective mass is strictly positive and coincides with the diffusion constant. Such a result is established also for polaron-type models with a suitable ultraviolet cut-off and for arbitrary values of the coupling constant. We point out a slightly stronger variant of the central limit theorem which would imply that the energy-momentum relation has a unique global minimum attained at zero momentum.

Equations

$$H=\tfrac{1}{2}p^{2}+\int_{\mathbb{R}^{d}}\mathrm{d}k\,\omega(k)a^{*}(k)a(k)+\sqrt{\alpha}\int_{\mathbb{R}^{d}}\mathrm{d}k\frac{\hat{v}(k)}{\sqrt{2\omega(k)}}\big{(}\mathrm{e}^{\mathrm{i}kx}a(k)+\mathrm{e}^{-\mathrm{i}kx}a^{*}(k)\big{)}.$$

$$P=p+P_{\mathrm{f}},P_{\mathrm{f}}={\int }_{{\mathbb{R}}^d}\mathrm{d}kk{a}^{\ast }(k)a(k),$$

$$H=\Pi^{*}\big{(}\int^{\oplus}_{\mathbb{R}^{d}}\mathrm{d}P\,H(P)\big{)}\Pi,$$

$$H(P)=\tfrac{1}{2}(P-P_{\mathrm{f}})^{2}+\int_{\mathbb{R}^{d}}\mathrm{d}k\,\omega(k)a^{*}(k)a(k)+\sqrt{\alpha}\int_{\mathbb{R}^{d}}\mathrm{d}k\,g(k)\big{(}a(k)+a^{*}(k)\big{)}$$

$$E(P)=\inf \mathrm{spec}(H(P)).$$

$$E(P)\ge E(0).$$

$$(m_{\mathrm{eff}}{)}^{-1}=E_{\mathrm{r}}^{\prime \prime }(0).$$

$$(m_{\mathrm{eff}}{)}^{-1}={\sigma }^2.$$

$$E_{\mathrm{r}}^{\prime \prime }(0)>0.$$

$$E(P)>E(0),P\ne 0.$$

$${\tilde{G}}_{T_{-},T_{+}}(k,t)=⟨{\varphi }_{-},{\mathrm{e}}^{-T_{-}H}{\mathrm{e}}^{\mathrm{i}kx}{\mathrm{e}}^{-tH}{\mathrm{e}}^{-\mathrm{i}kx}{\mathrm{e}}^{-T_{+}H}{\varphi }_{+}⟩,t,T_{-},T_{+}\ge 0.$$

$${\tilde{G}}_{T_{-},T_{+}}(k,t)={\int }_{{\mathbb{R}}^d}\mathrm{d}P\overline{{\hat{\phi }}_{-}(P)}{\hat{\phi }}_{+}(P)⟨\Omega ,{\mathrm{e}}^{-T_{-}H(P)}{\mathrm{e}}^{-tH(P+k)}{\mathrm{e}}^{-T_{+}H(P)}\Omega ⟩.$$

$$\mathbb{E}^{\mathrm{G}}\big{(}u(x,t)u(x^{\prime},t^{\prime})\big{)}=\int_{\mathbb{R}^{d}}\mathrm{d}k\frac{1}{(2\pi)^{d}}\frac{1}{2\omega(k)}\mathrm{e}^{\mathrm{i}k(x-x^{\prime})}\mathrm{e}^{-\omega(k)|t-t^{\prime}|}.$$

$$\displaystyle\tilde{G}_{T_{-},T_{+}}(k,t)=\mathbb{E}^{\mathrm{W}}\!\!\times\!\mathbb{E}^{\mathrm{G}}\Big{(}\varphi_{-}(q(-T_{-}))\varphi_{+}(q(T_{+}+t))\mathrm{e}^{-\mathrm{i}k(q(t)-q(0))}$$

$$\displaystyle\hskip 50.0pt\times\exp\Big{[}\sqrt{\alpha}\int_{-T_{-}}^{T_{+}+t}\mathrm{d}s\int_{\mathbb{R}^{d}}\mathrm{d}x\,v(x)u(q(s)-x,s)\Big{]}\Big{)},$$

$$\displaystyle\tilde{G}_{T_{-},T_{+}}(k,t)=\mathbb{E}^{\mathrm{W}}\Big{(}\varphi_{-}(q(-T_{-}))\varphi_{+}(q(T_{+}+t))\mathrm{e}^{-\mathrm{i}k(q(t)-q(0))}$$

$$\displaystyle\hskip 50.0pt\times\exp\Big{[}\tfrac{1}{2}\alpha\int_{-T_{-}}^{T_{+}+t}\mathrm{d}s\int_{-T_{-}}^{T_{+}+t}\mathrm{d}s^{\prime}\,W(q(s)-q(s^{\prime}),s-s^{\prime})\Big{]}\Big{)}$$

$$W(x,t)={\int }_{{\mathbb{R}}^d}\mathrm{d}k|g(k){|}^2{\mathrm{e}}^{\mathrm{i}kx}{\mathrm{e}}^{-\omega (k)|t|}.$$

$$W^{\mathrm{Fr}}(x,t)=|x{|}^{-1}{\mathrm{e}}^{-|t|},$$

$$W_{\kappa }(x,t)=|x{|}^{-1}\frac{2}{\pi }\arctan (\kappa |x|){\mathrm{e}}^{-|t|},$$

$$G_{T_{-},T_{+}}(k,t)={\tilde{G}}_{T_{-},T_{+}}(k,t)/{\tilde{G}}_{T_{-},T_{+}}(0,t).$$

$$\int_{\mathbb{R}^{d}}\mathrm{d}k|g(k)|^{2}\big{(}\sum_{j=1,2,3}\omega^{-j}\big{)}<\infty,\quad\int_{\mathbb{R}^{d}}\mathrm{d}k|g(k)|^{2}|k|^{2}\big{(}\sum_{j=2,4}\omega^{-j}\big{)}<\infty.$$

$$\displaystyle\tilde{G}_{0,T}(k,t)=\mathbb{E}^{\mathrm{W}}\Big{(}\delta(q(0))\mathrm{e}^{-\mathrm{i}kq(t)}\exp\Big{[}\tfrac{1}{2}\alpha\int_{0}^{T+t}\mathrm{d}s\int_{0}^{T+t}\mathrm{d}s^{\prime}W(q(s)-q(s^{\prime}),s-s^{\prime})\Big{]}\Big{)}$$

$$\underset{T\to \infty }{\lim }{G}_{0,T}(k,t)=G_{0,\infty }(k,t).$$

$$\underset{ϵ\to 0}{\lim }{G}_{0,\infty }(ϵk,{ϵ}^{-2}t)={\mathrm{e}}^{-\frac{1}{2}{\sigma }^2{k}^{2}t}$$

$${\tilde{G}}_{0,T}(k,t)=⟨{\varphi }_{-},{\mathrm{e}}^{\mathrm{i}kx}{\mathrm{e}}^{-tH}{\mathrm{e}}^{-\mathrm{i}kx}{\mathrm{e}}^{-TH}{\varphi }_{+}⟩$$

$$={\int }_{{\mathbb{R}}^d}\mathrm{d}P\delta (P)⟨\Omega ,{\mathrm{e}}^{-tH(P+k)}{\mathrm{e}}^{-TH(P)}\Omega ⟩=⟨\Omega ,{\mathrm{e}}^{-tH(k)}{\mathrm{e}}^{-TH(0)}\Omega ⟩,$$

$$G_{0,\infty }(k,t)=⟨\Omega ,{\mathrm{e}}^{-t(H(k)-E(0))}{\psi }_0⟩/⟨\Omega ,{\psi }_0⟩.$$

$$\tilde{G}_{0,0}(k,t)=\mathbb{E}^{\mathrm{W}}\Big{(}\delta(q(0))\mathrm{e}^{-\mathrm{i}kq(t)}\exp\Big{[}\tfrac{1}{2}\alpha\int_{0}^{t}\mathrm{d}s\int_{0}^{t}\mathrm{d}s^{\prime}W(q(s)-q(s^{\prime}),s-s^{\prime})\Big{]}\Big{)},$$

$$G_{0,0}(k,t)=⟨\Omega ,{\mathrm{e}}^{-tH(k)}\Omega ⟩/⟨\Omega ,{\mathrm{e}}^{-tH(0)}\Omega ⟩.$$

$$\underset{ϵ\to 0}{\lim }{G}_{0,0}(ϵk,{ϵ}^{-2}t)={\mathrm{e}}^{-\frac{1}{2}{\sigma }^2{k}^{2}t},$$

$$\displaystyle H(P)=\frac{1}{2}(P-P_{\mathrm{f}})^{2}+H_{\mathrm{f}}+\sqrt{\alpha}\int_{\mathbb{R}^{d}}\mathrm{d}k\,g(k)\big{(}a(k)+a^{*}(k)\big{)},$$

Full text

Effective mass of the polaron - revisited

Wojciech Dybalski, Herbert Spohn

Zentrum Mathematik and Physik Department, TUM,

Boltzmannstraße 3, 85747 Garching, Germany.

dybalski@ma.tum.de, spohn@tum.de

Abstract. Properties of the energy-momentum relation for the Fröhlich polaron are of continuing interest, especially for large values of the coupling constant. By combining spectral theory with the available results on the central limit theorem for the polaron path measure we prove that, except for an intermediate range of couplings, the inverse effective mass is strictly positive and coincides with the diffusion constant. Such a result is established also for polaron-type models with a suitable ultraviolet cut-off and for arbitrary values of the coupling constant. We point out a slightly stronger variant of the central limit theorem which would imply that the energy-momentum relation has a unique global minimum attained at zero momentum.

08.8.2019

1 Introduction

Polaron refers to an electron interacting with the lattice vibrations of a polar crystal, see [1, 7, 22] as a guide to the physics literature. In the conventional approximations the quantum Hamiltonian reads

[TABLE]

We use units in which the bare electron mass equals one. $x,p$ are position and momentum of the electron in $\mathbb{R}^{d}$, $a^{*}(k),a(k)$ are the creation and annihilation operators of a free scalar Bose field over $\mathbb{R}^{d}$ with commutation relations $[a(k),a^{*}(k^{\prime})]=\delta(k-k^{\prime})$, $\omega$ is the dispersion relation of the Bose field, $\omega\geq 0$, continuous, and $\omega(Rk)=\omega(k)$ for all rotations $R$. The form factor $\hat{v}$ is assumed to be real rotation invariant and has the inverse Fourier transform111We use the convention $\hat{v}(k):=\frac{1}{(2\pi)^{d/2}}\int_{\mathbb{R}^{d}}dx\,\mathrm{e}^{-\mathrm{i}kx}v(x)$ and occasionally write $(Fv)(k):=\hat{v}(k)$. $\hat{v}$. $v(x)$ physically describes the smearing of the interaction between the electron and the Bose field. It is a standard convention to call $g=\hat{v}/\sqrt{2\omega}$ the coupling function. Finally $\alpha$ is the coupling constant, $\alpha\geq 0$. Formally $H$ acts on the Hilbert space $\mathcal{H}=L^{2}(\mathbb{R}^{d})\otimes\mathcal{F}$ with $\mathcal{F}$ the Fock space of the Bose field. The coupling between field and particle is translation invariant and hence the total momentum

[TABLE]

is conserved.

The Fröhlich polaron corresponds to the specific choice $d=3$, $\omega(k)=1$, and $g(k)=(\sqrt{2}\pi|k|)^{-1}$. In particular $\|g\|_{2}=\infty$ because of ultraviolet divergence. Since the strong coupling physics is dominated by the large $|k|$ behavior of the coupling function, no ultraviolet cut-off can be afforded and a separate discussion is required, see Section 4. The acoustic polaron corresponds to $\omega(k)=|k|$ and other variations can be found in the physics literature. As common practice [22, 14], we thus keep $d,\omega,g$ general for a while and add further assumptions on the way. For the purpose of this introductory discussion $\|g\|_{2}<\infty$ is assumed. Under precise conditions to be stated in Section 3, $H$ is a self-adjoint operator and has the fiber decomposition

[TABLE]

since $P$ is conserved. (Here $\Pi:=F\mathrm{e}^{\mathrm{i}P_{\mathrm{f}}x}$, where $F$ is the Fourier transform from $x$ to $P$ variable). The fiber Hamiltonian reads

[TABLE]

and acts on $\mathcal{F}$. The energy-momentum relation, $E(P)$, is the bottom of the spectrum of $H(P)$,

[TABLE]

By construction $E(RP)=E(P)$ for all rotations $R$ and hence $E(P)=E_{\mathrm{r}}(|P|)$. Also, in generality,

[TABLE]

As a widely accepted definition, the effective mass is the inverse of the curvature of $E(P)$ at $P=0$, which by rotation invariance means

[TABLE]

A long-standing open problem is to analyse the effective mass of the Fröhlich polaron in the strong coupling regime [19, 22]. As the coupling is increased, more and more bosons are bound to the electron and one would expect the effective mass to increase with $\alpha$, presumably to diverge in the limit. Recent progress has been achieved by Lieb and Seiringer [12], who by functional analytic methods prove that indeed $m_{\mathrm{eff}}(\alpha)$ diverges as $\alpha\to\infty$. A mathematically orthogonal approach is to study the polaron path measure, as originally introduced by Feynman [4], see also [5, 21]. Mathematically this corresponds to a standard Brownian motion with a Gibbsian-like weight which depends only on the increments. Thus one expects to still observe diffusive behavior on large scales, however with an effective diffusion constant, $\sigma^{2}$. In other words, one conjectures the validity of a central limit theorem (CLT) for such weighted Brownian motion. In fact, the CLT has become now available for a large class of polaron models, including the assertion that $\sigma^{2}>0$. We refer to Section 2 for more details.

As argued in [22], and presumably before, effective mass and diffusion constant should be related as

[TABLE]

However, at the time the reasoning was based on considering second moments for the position of the weighted Brownian motion, a piece of information which is not so easily available from current CLT proofs. For us this by itself is a convincingly enough reason to reconsider the case. As an extra bonus, apparently not noted before, the conventional CLT also yields spectral information about properties apparently unaccessible by current functional analytic techniques. To explain this point in more detail we assume the lower bound $\omega(k)\geq c_{0}>0$ and further conditions as stated in Sections 3 and 4. Then $H(P)$ has a unique ground state with energy $E(P)$ for $|P|<\gamma$ for some $\gamma>0$. Furthermore the continuum edge of $H(P)$ is strictly larger than $E(P)$ in this ball. Possible eigenvalues have a finite multiplicity and can accumulate only at the continuum edge [6, 14]. In particular, $E_{\mathrm{r}}$ is real analytic in $P$ and $\sqrt{\alpha}$. From perturbation theory, requiring $\alpha$ to be sufficiently small, one then infers that

[TABLE]

But for larger $\alpha$ it is difficult to exclude $E(P)$ to have vanishing curvature at $P=0$. Under to be stated conditions we will establish the identity (1.9). Hence from $\sigma^{2}>0$ one concludes that $m_{\mathrm{eff}}<\infty$ for all $\alpha$.

A related issue is the long-standing (physically obvious) conjecture

[TABLE]

The weaker property (1.6) follows from a Kato inequality for $H$ [6]. So the real issue is to exclude points $P^{*}\neq 0$ at which $E(P^{*})=E(0)$. Property 1.10 is claimed in [7], Statement 2. In the proof on p. 78, the authors argue with the degeneracy of the ground state of $H$. But, under the common assumptions, $H$ has no ground state at all. To have a ground state would require the set $\{P\in\mathbb{R}^{d}\,|\,E(P)=E(0)\}$ to have non-zero Lebesgue measure. In Section 5 we will explain how (1.10) follows from a CLT with yet to be studied boundary conditions. Alternatively one might invoke a suitable large deviation result in the context of the available boundary conditions.

Our paper is organized as follows. In Section 2 we explain the connection to the probabilistic CLT and discuss recent results of interest in our context. In Sections 3 and 4 we show the relations (1.8), (1.9) for polaron-type models with the UV cut-off $\|g\|_{2}<\infty$ and for the Fröhlich polaron, respectively. In Section 5 we study the functional analytic side of a CLT with general two-sided pinning.

2 Probabilistic approach and central limit theorem

We choose the boundary states $\phi_{\pm}=\varphi_{\pm}\otimes\Omega\in\mathcal{H}$ with $\Omega$ the Fock vacuum and define

[TABLE]

Using the direct inegral decompostion (1.3), (1.4), see also formula (B.1), one obtains the identity

[TABLE]

In spirit of the Feynman-Kac formula for a Schrödinger operator, the semigroup $\mathrm{e}^{-tH}$, $t\geq 0$, can be written as a weighted average with respect to a Gaussian measure. For the particle trajectories we introduce the Wiener measure $\mathbb{P}^{\mathrm{W}}$ with expectation $\mathbb{E}^{\mathrm{W}}$. The continuous paths of the Wiener measure are denoted by $q(t)$. The Bose field maps to the Gaussian process $u(x,t)$ whose path measure is denoted by $\mathbb{P}^{\mathrm{G}}$, with expectation $\mathbb{E}^{\mathrm{G}}$. The Gaussian process has mean zero, is stationary in space-time, and is uniquely defined through its covariance

[TABLE]

Then

[TABLE]

which makes more explicit how $v(x)$ smears the field $u$ relative to the position of the particle. If $\|g\|_{2}<\infty$, the term in the square brackets is a well-defined Gaussian random variable with respect to $\mathbb{P}^{\mathrm{G}}$. The Gaussian average $\mathbb{E}^{\mathrm{G}}$ can be carried out explicitly leading to

[TABLE]

with

[TABLE]

Note that $W$ is real, continuous, rotation invariant in $x$, and $|W(x,t)|\leq\|g\|^{2}_{2}$. In particular the integrand under the double time integral appearing in (2) is pathwise bounded and continuous.

The Fröhlich polaron is the special case $d=3$, $\omega=1$, and $g(k)=(\sqrt{2\pi}|k|)^{-1}$, thereby defining the Hamiltonian $H^{\mathrm{Fr}}$, for which self-adjointness is established in [6, 8]. The kernel $W$ of the Fröhlich polaron is given by

[TABLE]

which is no longer bounded. Still the factor $\exp[\cdot]$ in (2) is integrable [3]. To establish the validity of the basic identity (2) for $H^{\mathrm{Fr}}$ one introduces the cut-off coupling $g_{\kappa}(k)=(\sqrt{2}\pi|k|)^{-1}\mathrm{e}^{-\frac{1}{\kappa}|k|}$, thereby defining the Hamiltonian $H_{\kappa}$. The strong limit $\lim_{\kappa\to\infty}\mathrm{e}^{-tH_{\kappa}}=\mathrm{e}^{-tH^{\mathrm{Fr}}}$ is established in [6, 13], which controls the left side of (2). On the right side $W$ is replaced by

[TABLE]

which increases monotonously to $W^{\mathrm{Fr}}(x,t)$. Thus by monotonicity the right hand side of (2) converges to the corresponding expression with kernel $W^{\mathrm{Fr}}(x,t)$ and hence (2) remains valid for the Fröhlich polaron.

In (2) the reference process is a standard Brownian motion over the time interval $[-T_{-},t+T_{+}]$. The Brownian motion is pinned by the function $\varphi_{-}$ at the left border and by $\varphi_{+}$ at the right one. The Brownian path is weighted by the exponential of the double time integral involving $W$. Note that the weight depends only on the increments. To have a probability measure we have to normalize by the partition function $\tilde{G}_{T_{-},T_{+}}(0,t)$. The difference $q(t)-q(0)$ is the Brownian motion increment over the time interval $[0,t]$. Of interest is its characteristic function, i.e. the Fourier transform of the corresponding probability density function. Altogether this leads to the normalized characteristic function

[TABLE]

Depending on the precise set-up, one then has to establish the limits $T_{-},T_{+}\to\infty$ followed by the CLT which requires $t\to\infty$.

In the probabilistic literature, two distinct boundary conditions have been studied and we discuss them one by one. In both cases $T_{-}=0$, $\varphi_{-}(x)=\delta(x)$, and $\varphi_{+}(x)=1$, of which the latter two have to be approximated by a suitable sequence of $L^{2}$ functions. We set $T_{+}=T$ in the sequel.

In [2] and the follow-up by Gubinelli [10] the authors require the conditions

[TABLE]

They consider $\tilde{G}_{0,T}$ of the form

[TABLE]

and establish the limit

[TABLE]

The CLT is proved, thus ensuring the limit

[TABLE]

for some $\sigma>0$. In fact, the stronger functional CLT is established, see [2, Theorem 1.1].

It is instructive to rewrite the expectation values from above in the language of operators as in (2.1), (2.2), with the result

[TABLE]

where we used $\hat{\varphi}_{-}(P)=(2\pi)^{-d/2}$, and $\hat{\varphi}_{+}(P)=(2\pi)^{d/2}\delta(P)$. For polaron-type models treated in Section 3 below $H(0)$ has a spectral gap and a unique ground state $\psi_{0}$, thus by the spectral theorem

[TABLE]

More recently, Mukerjee and Varadhan studied the CLT under weaker conditions than imposed in [2, 10]. Their starting formula is

[TABLE]

hence $T=0$, which one recognizes as a particular case of (2) and

[TABLE]

In [16, Theorem 4.2] the CLT of the following form is established for the Fröhlich polaron,

[TABLE]

for some $\sigma>0$, with the restriction $\alpha\in[0,\alpha_{0})\cup(\alpha_{1},\infty)$ for some $0<\alpha_{0}<\alpha_{1}<\infty$.

The functional CLT is not touched upon. In the related study [17] the strong coupling limit and its relation to the Pekar process are investigated. Mukherjee [18] also starts from (2.16) and considers a general weight function $W$, for which he requires $|W(x,t)|\leq C(1+|t|)^{-(2+\delta)}$ for some $C,\delta>0$. In particular, this condition covers the polaron whenever $g\in L^{2}$. [ In the currently posted version in addition $W\geq 0$ is required. As communicated to us by the author this condition can be dropped.] In [18, Theorem 2.1] the conventional CLT of the form (2.18) is proved for arbitrary $\alpha\geq 0$.

Physically one is also interested in the behavior of $E(P)$ away from the origin. Starting from (2.17), instead of $k=\mathcal{O}(\epsilon)$ and would have to consider $k=P=\mathcal{O}(1)$, which probabilistically is a problem of large deviations. In Section 5 we explore a different approach by starting from (2.1) with general square-integrable boundary functions $\varphi_{\pm}$ in the limit $T_{\pm}\to\infty$, but still invoking a CLT.

3 Polaron-type models with a UV cut-off

In this section we show that $\sigma>0$ appearing in the CLT (2.18) coincides with the effective mass for a large class of polaron-type Hamiltonians with a UV cut-off. It is convenient to start from a family of the fiber Hamiltonians of the form

[TABLE]

where $H_{\mathrm{f}}=\int_{\mathbb{R}^{d}}\mathrm{d}k\,\omega(k)a^{*}(k)a(k)$, $P_{\mathrm{f}}=\int_{\mathbb{R}^{d}}\mathrm{d}k\,k\,a^{*}(k)a(k)$. Further asumptions are listed in

Condition C. (i) $g\in L^{2}(\mathbb{R}^{d})$ is real and rotation invariant. The coupling constant $\alpha\geq 0$ is arbitrary.

(ii) $\omega(k)\geq c_{0}>0$, $\omega$ is continuous, rotation invariant, and sub-additive in the sense that

[TABLE]

Then, by the Kato-Rellich theorem, $H(P)$ are self-adjoint, semi-bounded operators on the domain $D(P_{\mathrm{f}}^{2}+H_{\mathrm{f}})$ which is independent of $P$. By the direct integral formula (1.3) one obtains a Hamiltonian of the form (1.1). Under the above assumptions, the HVZ theorem for these models was shown in [6, 14]. All the properties below can be found in [14] except for part 0 for which we refer to [9] or [24, Section 15.2], and part 6 which can be found in [15]. We refer to [14] for a discussion of the literature.

Lemma 3.1

[14, 6]* Assume Condition $C$ and define $E(P)=\inf\,\mathrm{spec}(H(P))$, $E_{\mathrm{ess}}(P)=\inf\mathrm{spec}_{\mathrm{ess}}(H(P))$. Then the following statements hold true:*

$E(0)\leq E(P)$* for all $P\in\mathbb{R}^{d}$.* 2. 1.

$E_{\mathrm{ess}}(P)=\mathrm{inf}_{k\in\mathbb{R}^{d}}(E(P-k)+\omega(k))$. 3. 2.

The interval $\mathcal{I}_{0}:=\{\,P\in\mathbb{R}^{d}\,|\,E(P)<E_{\mathrm{ess}}(P)\,\}$ is non-empty and contains a neighbourhood of any global minimum of $|P|\mapsto E_{\mathrm{r}}(|P|)$. 4. 3.

$E(P)$* is an isolated, simple eigenvalue for $P\in\mathcal{I}_{0}$.* 5. 4.

If $\omega$ is bounded, then, for any $\overline{E}\in\mathbb{R}$, $\lim_{|P|\to\infty,E(P)\leq\overline{E}}(E_{\mathrm{ess}}(P)-E(P))=0$.

In general, $E(P)\geq c_{1}\omega(P)+c_{2}$ for some $c_{1}>0$. 6. 5.

For $P\in\mathcal{I}_{0}$ we have $|\langle\Omega,\psi_{P}\rangle|>0$, where $\psi_{P}$ is the ground state of $H(P)$. 7. 6.

$\mathcal{I}_{0}\ni P\mapsto E(P)$* and $P\mapsto|\psi_{P}\rangle\langle\psi_{P}|$ are real analytic functions.*

Let us comment briefly on the proof of properties 1–5 and the role of various assumptions on $\omega$. We consider the thresholds

[TABLE]

Assuming only that $\omega$ is continuous, bounded and massive (i.e. $\omega(k)\geq c_{0}>0$), Theorem 2.1 of [14] gives $E_{\mathrm{ess}}(P)=\inf_{n\geq 1}E^{(n)}(P)$. As sub-additivity of $\omega$ clearly gives monotonicity of thresholds, with this additional assumption one obtains property 1 of Lemma 3.1 above. As pointed out in [14], it is clear from this relation, and from the fact that $\omega$ is massive, that if $E(P)<\inf_{P^{\prime}}E(P^{\prime})+\inf_{k^{\prime}}\omega(k^{\prime})$ then $E(P)<E_{\mathrm{ess}}(P)$, which gives property 2 of Lemma 3.1. Clearly, the spectrum below $E_{\mathrm{ess}}(P)$ consists at most of eigenvalues of finite multiplicity with $E_{\mathrm{ess}}(P)$ as the only possible accumulation point. Thus $E(P)$ is an isolated eigenvalue, which is simple by Theorem 2.3 of [14]. Thus we obtain property 3 of Lemma 3.1. For the first part of property 4 and property 5 we refer to Theorems 2.3 and 2.4 of [14]. The second part of property 4 can be found in [6] (see also [24, Section 15.2, property (v)]).

As for part 6, we note that for $\xi$ in the resolvent set of $H(P_{0})$ the function $P\mapsto(H(P)-\xi)^{-1}$ can be expanded around any $P_{0}\in\mathbb{R}^{d}$ as in formula (A.2). The real analyticity of the eigenprojections $\mathcal{I}_{0}\ni P\mapsto|\psi_{P}\rangle\langle\psi_{P}|$ follows immediately via the Cauchy formula. (We note that by a suitable choice of the phase, we can ensure that $\mathcal{I}_{0}\ni P\mapsto\psi_{P}$ is norm-continuous, which is the property we will need below). Since $|P|\mapsto H(\hat{P}|P|)$ is a real analytic family of self-adjoint operators in the sense of [11, Chapter VII, §1, §3] and $P\mapsto E(P)$ is a rotation invariant function, we obtain by [11, Chapter VII, §3] that $\mathcal{I}_{0}\ni P\mapsto E(P)$ is real analytic.

Now we are ready to state and prove our main result concerning polaron-type models with a UV cut-off.

Theorem 3.2

Consider polaron-type models satisfying Condition C. Then, for all $\alpha\geq 0$,

[TABLE]

Proof: The proof relies on the CLT as stated in (2.18). We consider the expression

[TABLE]

The spectral calculus immediately gives $\lim_{\epsilon\to 0}\langle\Omega,\mathrm{e}^{-\frac{t}{\epsilon^{2}}(H(0)-E(0))}\Omega\rangle=\langle\Omega,\psi_{0}\rangle\langle\psi_{0},\Omega\rangle$. Concerning the numerator in (3.5), we obtain

[TABLE]

where in the leading term (3.6) we used the analyticity of $P\mapsto E(P)$ near zero (see Lemma 3.1) and we noted that the expression in (3.7) tends to zero as $\epsilon\to 0$ by the spectral calculus. $\Box$

We remark that a priori the diffusion constant obtained from the CLT of the characteristic function (2.17) could differ from the one of (2.13). Our result implies that they agree under Condition C and assumption (2.10).

4 The Fröhlich polaron

Let $H_{\kappa}(P)$ be the polaron Hamiltonian (3.1) with $d=3$, $\omega\equiv 1$ and $g(k)=\frac{\chi_{[0,\kappa]}(|k|)}{\sqrt{2}\pi|k|}$, where $\kappa$ is the UV cut-off222We use here a different UV cut-off than in the discussion in Section 2. However, the limiting Fröhlich Hamiltonians $H^{\mathrm{Fr}}(P)$ are the same, as one can infer from [13, Proposition A.4] and the strong convergence of the Gross transform.. Explicitly, it has the form

[TABLE]

where $N_{\mathrm{f}}$ is the number operator. It is well known, that this sequence of Hamiltonians converges in the norm resolvent sense as $\kappa\to\infty$ to the limiting Fröhlich Hamiltonian $H^{\mathrm{Fr}}(P)$. Also, the sequence of the full Hamiltonians $H_{\kappa}=\Pi^{*}\int^{\oplus}\mathrm{d}P\,H_{\kappa}(P)\Pi$ converges in the norm resolvent sense to $H^{\mathrm{Fr}}=\Pi^{*}\int^{\oplus}\mathrm{d}PH^{\mathrm{Fr}}(P)\,\Pi$, cf. [8] and references therein. Making use of these approximation properties, Lemma 3.1 and further results from [13], it is easy to establish the following:

Lemma 4.1

Let $E(P)=\inf\,\mathrm{spec}(H^{\mathrm{Fr}}(P))$ and $E_{\mathrm{ess}}(P)=\inf\mathrm{spec}_{\mathrm{ess}}(H^{\mathrm{Fr}}(P))$. Then the following statements hold true:

$E(0)\leq E(P)$* for all $P\in\mathbb{R}^{3}$.* 2. 1.

$E_{\mathrm{ess}}(P)=E(0)+1$. 3. 2.

The interval $\mathcal{I}_{0}:=\{\,P\in\mathbb{R}^{d}\,|\,E(P)<E_{\mathrm{ess}}(P)\,\}$ contains a neighbourhood of any global minimum of $|P|\mapsto E_{\mathrm{r}}(|P|)$. 4. 3.

$E(P)$* is an isolated, simple eigenvalue for $P\in\mathcal{I}_{0}$.* 5. 4.

All global minima of $|P|\mapsto E_{\mathrm{r}}(|P|)$ are contained in a compact set. 6. 5.

For $P\in\mathcal{I}_{0}$ we have $|\langle\Omega,\psi_{P}\rangle|>0$, where $\psi_{P}$ is the ground state of $H(P)$. 7. 6.

$\mathcal{I}_{0}\ni P\mapsto E(P)$* and $P\mapsto|\psi_{P}\rangle\langle\psi_{P}|$ are real analytic functions.*

Let us comment on the proofs of the above properties. It is a general consequence of the strong resolvent convergence that for any eigenvalue $\lambda$ of $H(P)$ there exists an approximating sequence $\lambda_{\kappa}\to\lambda$ of eigenvalues of $H_{\kappa}(P)$ [20, Theorem VIII.24]. Therefore, part 0 of Lemma 4.1 follows from part 0 of Lemma 3.1. Next, by [13, Proposition A.4], $\lim_{\kappa\to\infty}E_{\mathrm{ess},\kappa}(P)=E_{\mathrm{ess}}(P)$, where $E_{\mathrm{ess},\kappa}(P)$ is the bottom of the essential spectrum of $H_{\kappa}(P)$. Now part 1 of Lemma 4.1 follows from part 1 of Lemma 3.1 applied to the case of $\omega\equiv 1$. (Alternatively, one can refer to [23, Section IV]). Parts 2 and 3 of Lemma 4.1 follow from parts 1 and 2 of the same lemma, considering that the proof of [13, Theorem 6.4] gives the uniqueness of the ground state whenever it exists, also outside of the ball $|P|<\sqrt{2}$ from the statement of the theorem. Concerning part 4, suppose by contradiction that there is a sequence $P_{\ell}$, $\ell\in\mathbb{N}$, s.t. $E(P_{\ell})=E(0)$ and $|P_{\ell}|\to\infty$. We pick a function $f\in C_{0}^{\infty}(\mathbb{R})$ supported in a ball around $E(0)$ of radius strictly smaller than $1$ and s.t. $0\leq f\leq 1$ and $f(E(0))=1$. Then, by the norm resolvent convergence of $H_{\kappa}$ and [20, Theorem VIII.20] we have

[TABLE]

which is a contradiction. Here in the third step we choose $\ell_{\kappa}$ so large that the spectrum of $H_{\kappa}(P_{\ell})$ is outside of the support of $f$, which is possible by Lemma 3.1, part 4. Part 5 of Lemma 4.1 is a consequence of the strict positivity statement in [13, Theorem 6.4], where again we can disregard the restriction $|P|<\sqrt{2}$, considering the structure of the proof. Part 6 is proven analogously as the corresponding part of Lemma 3.1, given the input from Appendix A.

Now we come to our main result concerning the Fröhlich polaron.

Theorem 4.2

Consider the Fröhlich polaron. Then, for all $\alpha\in[0,\alpha_{0})\cup(\alpha_{1},\infty)$ for some $0<\alpha_{0}<\alpha_{1}<\infty$,

[TABLE]

Proof: The claim follows from the CLT stated in (2.18) by the same steps as in the proof of Theorem 3.2. Instead of Lemma 3.1, Lemma 4.1 is used. $\Box$

5 A CLT for two-sided pinning

We return to the set-up of equation (2.1) with square-integrable boundary functions $\varphi_{\pm}=\varphi$, $T_{\pm}=T$, and $T\to\infty$. A probabilistic study of this variant does not seem to be available in the literature and we focus on the functional analytic side. It will be more transparent to work in a general framework, which includes the polaron models discussed so far, but many more, e.g. systems with a non-quadratic energy momentum relation for the electron.

Let $\mathcal{H}$, $\mathfrak{H}$ be Hilbert spaces and $\Pi:\mathcal{H}\to L^{2}(\mathbb{R}^{d},\mathfrak{H})$ a unitary. For any $\phi\in\mathcal{H}$ we have the corresponding representation $\phi=\Pi^{*}\int^{\oplus}_{\mathbb{R}^{d}}\mathrm{d}P\,\phi_{P}$. For any $k\in\mathbb{R}^{d}$ we define the unitary $U(k)$ by its action on such vectors $\phi$

[TABLE]

Furthermore, we are interested in self-adjoint operators $H$ on a domain $D(H)\subset\mathcal{H}$ which have the representation

[TABLE]

Here $\mathbb{R}^{d}\ni P\mapsto H(P)$ is a real analytic family of positive operators with domains $D(H(P))\subset\mathfrak{H}$, as stated more precisely in the standing assumption 0 below. Furthermore, we note that for any bounded Borel function $f$

[TABLE]

In this section we impose the following standing assumptions:

The family $P\mapsto H(P)$ is real analytic in the sense that for any $P_{0}\in\mathbb{R}^{d}$ and any $\xi\notin\mathrm{spec}(H(P_{0}))$ there exists a real neighbourhood $N_{P_{0}}$ of $P_{0}$ s.t. $\xi\notin\mathrm{spec}(H(P))$ for any $P\in N_{P_{0}}$ and $N_{P_{0}}\ni P\mapsto(H(P)-\xi)^{-1}$ is real analytic. As a consequence, for any $\hat{P}$ on the unit sphere $|P|\mapsto H(\hat{P}|P|)$ is a real analytic family of unbounded operators in the sense of [11, Chapter VII, §1]. Another consequence of this property and of the Helfer-Sjöstrand method of almost analytic extensions is the strong continuity of $\mathbb{R}^{d}\ni P\mapsto\mathrm{e}^{-tH(P)}$, which will be used in the proofs below. 2. 1.

The function $P\mapsto E(P):=\mathrm{infspec}(H(P))$ is rotation invariant and we write as before $E(P)=E_{\mathrm{r}}(|P|)$. $E$ attains its global minima in the sets $M_{\ell}=\{\,P\in\mathbb{R}^{d}\,|\,|P|=Q_{\ell}\,\}$, $\ell=0,1,2\ldots,L$, where $0\leq Q_{0}<Q_{1}<\ldots<Q_{L}$ and $L$ finite. Also, we assume $E_{\mathrm{r}}(Q_{\ell})=0$. 3. 2.

$E$ is analytic in sets $\tilde{M}_{\ell}=\{\,P\in\mathbb{R}^{d}\,|\,|P|\in\Delta_{\ell}\,\}$, where $\Delta_{\ell}$ is a neighbourhood of $Q_{\ell}$. Then we have E_{\mathrm{r}}({\color[rgb]{0,0,0}Q_{\ell}+R})\sim R^{n_{\ell}} for small $R$ and some $n_{\ell}\in\mathbb{N}$, $n_{\ell}\geq 2$. 4. 3.

For $P\in\tilde{M}_{\ell}$, $E(P)$ are simple eigenvalues and the corresponding family of projections $\tilde{M}_{\ell}\ni P\mapsto|\psi_{P}\rangle\langle\psi_{P}|$ is strongly continuous. (It easily follows that $P\mapsto\psi_{P}$ can be chosen strongly continuous in $\tilde{M}_{\ell}$ by a suitable choice of the phases, which is what will be used below). 5. 4.

There exist vectors $\phi\in\mathcal{H}$ such that $\tilde{M}_{\ell}\ni P\mapsto|\langle\phi_{P},\psi_{P}\rangle|$ are continuous and non-zero on $M_{\ell}$.

The above assumptions hold, in particular, for models of Sections 3 and 4 as shown in the following proposition. The proof is postponed to Appendix B.

Proposition 5.1

For the polaron-type models (3.1) satisfying Condition C and the Fröhlich polaron (4.1) the following properties hold true:

1. (a)

The models satisfy the standing assumptions 0,1,2,3 above. 2. (b)

Let $\phi=\varphi\otimes\Omega\in\mathcal{H}$ be s.t. $\varphi\in L^{2}(\mathbb{R}^{d})$, $\hat{\varphi}\in C(\mathbb{R}^{d})$ and $\hat{\varphi}(p)>0$ for all $p\in\mathcal{I}_{0}$. For such $\phi$ assumption 4 above holds.

Coming back to the general framework, we note that assumptions 2, 3 follow from 0,1 and analytic perturbation theory [11, Chapter VII, §3] if $E(P)$ are simple, isolated eigenvalues for $P\in\tilde{M}_{\ell}$. However, our discussion in this section does not require spectral gaps above $E(P)$. Also the standard relation $E(0)\leq E(P)$, $P\in\mathbb{R}^{d}$, for polaron type models, cf. Lemmas 3.1, 4.1, does not follow from the standing assumptions above. However, with additional input which we now explain, we will obtain not only this relation, but even $E(0)<E(P)$, $P\neq 0$, for $d\geq 2$.

For the two-sided boundary condition the properly normalized characteristic function reads

[TABLE]

for $\phi\in\mathcal{H}$ and $t,T\geq 0$. By the spectral theorem, the denominator above is different from zero for any finite $T$. Furthermore, for $\phi$ as in assumption 4, the limits

[TABLE]

exist under our standing assumptions. The explicit expressions are provided in Proposition 5.4 and Lemma 5.5 below. We expect that the latter limit has the form suggested by the CLT.

Conjecture 5.2

There exists $\phi\in\mathcal{H}$ as in assumption 4 above, such that the CLT of the form

[TABLE]

holds true for some $\sigma>0$.

The consequences of this conjecture for models satisfying the above standing assumptions are collected in the following theorem.

Theorem 5.3

Suppose that Conjecture 5.2 holds true for some $\phi\in\mathcal{H}$ as in assumption 4 and $\sigma>0$. Then, for $d\geq 2$, there is a global minimum at zero (i.e. $Q_{0}=0$). Furthermore,

1. (a)

$\sigma^{2}=(\partial^{2}_{|P|}E_{\mathrm{r}})(0)$, 2. (b)

$E(P)>E(0)$* for $P\neq 0$.*

For $d=1$ we obtain that $\sigma^{2}=(\partial^{2}_{|P|}E_{\mathrm{r}})(Q_{\ell})$ for $\ell=0,1,2,\ldots L$ and $0\leq Q_{0}<Q_{1}<\ldots<~{}Q_{L}$.

We stated a minimal conjecture as required for Theorem 5.3 to hold. In fact, the CLT should be in force for a large set of boundary functions, e.g. those satisfying 4. of the standing assumptions. Our theorem then asserts that the diffusion constant is always given by $\sigma^{2}=(\partial^{2}_{|P|}E_{\mathrm{r}})(0)$.

The observation behind Theorem 5.3 is fairly elementary and can be grasped most easily for the polaron models underlying (2.2). We note that the $P$-integral has the weight $|\hat{\varphi}(P)|^{2}>0$ for $\varphi$ as in Proposition 5.1(b). Thus in the limit $T\to\infty$ the $P$-integral concentrates on the set of global minima $\{P\in\mathbb{R}^{d}\,|\,E(P)=E(0)\}$. If the CLT would hold, the limit expression must have come only from $P=0$ and hence $E(P)>E(0)$ for $P\neq 0$.

The actual proof is more involved and the remaining part of this section is devoted to proving Theorem 5.3. We start with two auxiliary results, which do not rely on Conjecture 5.2.

Proposition 5.4

The following statements hold:

If $Q_{0}=0$ and it is the only global minimum of $E$, then

[TABLE] 2. 2.

If $Q_{0}=0$ and there are other global minima at $Q_{\ell}>0$, $\ell=1,2,\ldots,L$, we set $n:=\mathrm{max}_{\ell\neq 0}(n_{\ell})$ (see assumption 2) and distinguish the following cases:

1. (a)

For $n>\frac{n_{0}}{d}$

[TABLE]

where $C_{\ell}\neq 0$ and the sums extend only over $\ell>0$ s.t. $n_{\ell}=n$. 2. (b)

For $n=\frac{n_{0}}{d}$

[TABLE]

where $c_{0},c_{\ell}\neq 0$ and the sums extend only over $\ell>0$ s.t. $n_{\ell}=n$. 3. (c)

For $n<\frac{n_{0}}{d}$

[TABLE] 3. 3.

If $0<Q_{0}<Q_{1}<\ldots<Q_{L}$, for $L\geq 0$, we obtain

[TABLE]

where $C_{\ell}\neq 0$ and the sum extends over $\ell$ s.t. $n_{\ell}=\bar{n}:=\max_{\ell^{\prime}}{n_{\ell^{\prime}}}$.

For $d=1$ the angular integrations above amount to summation over $\hat{P}=\pm 1$.

Proof: In the fiber representation expression (5.4) has the following form

[TABLE]

We denote the spectral measure of $H(P)$ by $M_{P}(\,\cdot\,)$ and choose $\delta>0$ s.t. $E(P)\leq\delta$ implies that $P\in\tilde{M}:=\bigcup_{\ell=0}^{L}\tilde{M}_{\ell}$. By spectral calculus, we have

[TABLE]

Therefore, it suffices to study

[TABLE]

Hence, setting ${G}^{(1)}_{T}(k,t):=\tilde{G}^{(1)}_{T}(k,t)/\tilde{G}^{(1)}_{T}(0,t)$,

[TABLE]

We write $\hat{P}:=P/|P|$ and move on to polar coordinates in $P$ and $P^{\prime}$ integrations:

[TABLE]

where we also used that $P\mapsto E(P)$ is rotation invariant.

Let us first consider a possible global minimum at zero. Since $E$ is analytic near zero, we have that $E_{\mathrm{r}}(|P|)\sim|P|^{n_{0}}$ in this region for some $n_{0}\in\mathbb{N}_{0}$, $n_{0}\geq 2$. Thus an elementary analysis gives for the numerator in (5.16)

[TABLE]

for some $C_{0}\neq 0$. An analogous formula holds for the denominator in (5.16)

[TABLE]

Let us now analyse a global minimum at $Q_{\ell}\neq 0$. By analyticity, we have that E_{\mathrm{r}}({\color[rgb]{0,0,0}Q_{\ell}+R})\sim R^{n_{\ell}} near $R=0$ for some $n_{\ell}\in\mathbb{N}_{0}$, $n_{\ell}\geq 2$. In this case we obtain for the numerator in (5.16)

[TABLE]

for some $C_{\ell}\neq 0$. For the denominator in (5.16) we get in this case

[TABLE]

By substituting (5.17)–(5.20) back to formula (5.16) and considering the different cases from the statement of the proposition, we complete the proof. $\Box$

Lemma 5.5

Suppose that $Q\geq 0$ is a global minimum of $|P|\mapsto E_{\mathrm{r}}(|P|)$. Then the following relations hold:

[TABLE]

Proof: Suppose that $Q=0$. By shifting the vector $\psi_{0}=\psi_{\epsilon k}+(\psi_{0}-\psi_{\epsilon k})$, we obtain

[TABLE]

where $R(\epsilon)$ is an error term which tends to zero as $\epsilon\to 0$ by assumption 3.

Concerning the case $Q>0$, we shift the vector as follows $\psi_{Q\hat{P}}=\psi_{Q\hat{P}+\epsilon k}+(\psi_{Q\hat{P}}-\psi_{Q\hat{P}+\epsilon k})$. This gives

[TABLE]

which completes the proof. $\Box$

Proof of Theorem 5.3: We start from the case $d\geq 2$. Suppose, by contradiction, that there is no global minimum at zero. Then, from the last part of Proposition 5.4 and Lemma 5.5 we obtain

[TABLE]

We denote $x^{2}:=tk^{2}\sigma^{2}/2$, $f_{\ell}(\hat{P}):=C_{\ell}|\langle\phi_{Q_{\ell}\hat{P}},\psi_{Q_{\ell}\hat{P}}\rangle|^{2}$, $m^{-1}_{\ell}:=\partial^{2}_{|P|}E(Q_{\ell})$. This gives

[TABLE]

By averaging both sides w.r.t. the group of rotations we can assume that the functions $f_{\ell}$ are constant and non-zero. Suppose first that all $m_{\ell}^{-1}$ are zero. Then we immediately obtain a contradiction by taking $x^{2}\to\infty$. Now suppose that some333We note as an aside, that if some $m_{\ell}^{-1}>0$ then all $m_{\ell}^{-1}>0$ by definition of $n$. $m_{\ell_{1}}^{-1}>0$. Then we obtain from (5.26)

[TABLE]

As before, we obtain a contradiction by taking $x^{2}\to\infty$, due to the fact that the functions $f_{\ell}$ are constant and non-zero.

Next, we prove part (b). Suppose, by contradiction, that there are several global minima in addition to the global minimum at zero. Let us assume first that case (a) from Proposition 5.4 occurs, that is $n>\frac{n_{0}}{d}$. With the help of Lemma 5.5, we obtain from (5.8)

[TABLE]

We obtain a contradiction by repeating the steps (5.25)–(5.27) above.

Now let us assume that case (b) of Proposition 5.4 occurs, that is $n=\frac{n_{0}}{d}$. Since $d\geq 2$ and $n\geq 2$ we obtain that $n_{0}>2$ which implies $(\partial^{2}_{|P|}E_{\mathrm{r}})(0)=0$. With the help of Lemma 5.5, we obtain from (LABEL:b) using the notation introduced above

[TABLE]

Due to the presence of the non-zero terms involving $c_{0}|\langle\phi_{0},\psi_{0}\rangle|^{2}$, we immediately obtain a contradiction by taking $x^{2}\to\infty$.

Finally, suppose that we are in case (c) of Proposition 5.4, that is $n<\frac{n_{0}}{d}$. Also in this case we have $n_{0}>2$ which implies $(\partial^{2}_{|P|}E_{\mathrm{r}})(0)=0$. By equation (5.10) and Lemma 5.5, we obtain $\mathrm{e}^{-x^{2}}=1$ which is immediately a contradiction. This concludes the proof of part (b) of the theorem.

Given that we have only one global minimum, we can apply formula (5.7). Together with Lemma 5.5, we obtain

[TABLE]

which gives part (a) of the theorem.

For $d=1$, the reasoning above requires several modifications. From the assumption that there is no global minimum at zero we obtain via (5.26) that $m^{-1}_{\ell}=\sigma^{2}$ for $\ell=0,1,\ldots,L$, which is what we wanted to prove.

Now suppose that there is a global minimum at zero and possibly some non-zero global minima. In the case $n>\frac{n_{0}}{d}$ from the fact that $n_{0}\geq 2$ we conclude that $n>2$, hence $m^{-1}_{\ell}=0$ for $\ell\neq 0$. In this situation formula (5.26) gives directly a contradiction.

In the case $n=\frac{n_{0}}{d}$ we distinguish two sub-cases. First, for $n=n_{0}>2$ we have $m_{\ell}^{-1}=0$ for all $\ell$, including $\ell=0$, and thus formula (5.29) gives a contradiction. Second, for $n=n_{0}=2$ we have $m_{\ell}^{-1}\neq 0$ and thus formula (5.29) has to be rewritten as follows

[TABLE]

where angular integration denotes now summation over $\hat{P}=\pm 1$. Clearly, we avoid a contradiction iff $m_{\ell}^{-1}=\sigma^{2}$ for $\ell=0,1,\ldots,L$. (It is important here that if $n=\mathrm{max}_{\ell\neq 0}(n_{\ell})=2$ then $n_{\ell}=2$ for all $\ell\neq 0$).

In the case $n<\frac{n_{0}}{d}$ we obtain a contradiction as before. Thus in the case $d=1$ the assumption that there are several global minima of $E_{\mathrm{r}}$ led us to the conclusion that the inverses of their effective masses $\partial^{2}_{|P|}E_{\mathrm{r}}(Q_{\ell})$, $\ell=0,1,\ldots L$, must all be equal to $\sigma^{2}$. $\Box$

Acknowledgements. We thank Fumio Hiroshima, Tadahiro Miyao, Chiranjib Mukherjee, and Jacob Schach Møller for helpful discussions. This work was partially supported by the Deutsche Forschungsgemeinschaft (DFG) within the grant DY107/2-1.

Appendix A Analyticity of the Fröhlich polaron in total momentum

In this appendix we verify that $P\mapsto H^{\mathrm{Fr}}(P)$ is a real analytic family as specified in the standing assumption 0 of Section 5.

Lemma A.1

Suppose that $\xi\notin\mathrm{spec}(H(P_{0}))$. Then $\xi\notin\mathrm{spec}(H^{\mathrm{Fr}}(P))$ for $P$ in a neighbourhood $N_{P_{0}}$ of $P_{0}$. The function $N_{P_{0}}\ni P\mapsto(H^{\mathrm{Fr}}(P)-\xi)^{-1}$ is real analytic.

Proof: First, suppose that $\xi\notin\mathrm{spec}(H_{\kappa}(P))$ for $\kappa$ sufficiently large and note that on $D(P_{\mathrm{f}}^{2}+N_{\mathrm{f}})$

[TABLE]

Consequently, for $P$ in a small neighbourhood $N_{P_{0}}$ of $P_{0}$, which a priori may depend on $\kappa$, the series on the r.h.s. below converges and defines the inverse of $(H_{\kappa}(P)-\xi)$:

[TABLE]

To eliminate the dependence of $N_{P_{0}}$ on $\kappa$, we show in Lemma A.2 below, that

[TABLE]

uniformly in $\kappa$.

Now we intend to take the limit $\kappa\to\infty$ on both sides of (A.2). By [20, Theorem VIII.23], if $\xi\notin\mathrm{spec}(H^{\mathrm{Fr}}(P))$ then $\xi\notin\mathrm{spec}(H_{\kappa}{(P)})$ for $\kappa$ sufficiently large and $(H_{\kappa}{(P)}-\xi)^{-1}\to(H^{\mathrm{Fr}}{(P)}-\xi)^{-1}$ in norm. As the same is true for $P$ replaced with $P_{0}$ and we can use (A.3) to exchange the limit $\kappa\to\infty$ with summation in (A.2). Thus the proof is complete. $\Box$

Lemma A.2

The following bounds hold uniformly in $\kappa$:

[TABLE]

Proof: First, we recall some material from [13, Appendix A], referring there for more details. Let

[TABLE]

where $K$ is chosen sufficiently large (depending on $\alpha$ but not on $\kappa$ or $P$) as specified above Lemma A.3 of [13]. Let $H^{\mathrm{free}}(P)$ denote the Hamiltonians (4.1) with $\alpha=0$. Then the Gross-transformed Hamiltonians \tilde{H}_{\kappa}{\color[rgb]{0,0,0}(P)}:=e^{T_{K,\kappa}}H_{\kappa}(P)e^{-T_{K,\kappa}} are self-adjoint operators which converge in the norm-resolvent sense to a limiting Hamiltonian \tilde{H}({\color[rgb]{0,0,0}P}) [13, Proposition A.4]. As stated in the proof of this latter proposition, H^{\mathrm{free}}(P)\leq C^{\prime}(\tilde{H}_{\kappa}{\color[rgb]{0,0,0}(P)}+C), with $C,C^{\prime}$ independent of $\kappa$. Hence,

[TABLE]

Next, we can write on $D(H^{\mathrm{free}}(P))$

[TABLE]

where the last term is actually zero by symmetry. Noting that $\int\mathrm{d}k\,\,|k\beta_{K,\kappa}(k)|^{2}\leq c$, uniformly in $\kappa$, we have $\|a^{(*)}(k\beta_{K,\kappa})(1+N_{\mathrm{f}})^{-1/2}\|\leq c$ uniformly in $\kappa$. Therefore, by estimates (A.6),

[TABLE]

is uniformly bounded in $\kappa$. This concludes the proof. $\Box$

Appendix B Proof of Proposition 5.1

Let us consider first polaron-type models satisfying Condition C. For the standing assumption 0 we refer to the discussion of Lemma 3.1, part 6. By part 2 of Lemma 3.1, $\mathcal{I}_{0}$ contains neighbourhoods of all the global minima of $|P|\mapsto E_{\mathrm{r}}(|P|)$. Thus, considering other items of this lemma, it suffices to show that there is a finite number of such minima to complete the proof of Proposition 5.1 (a). To this end, we first note that by part 1 of Lemma 3.1, $P\mapsto E_{\mathrm{ess}}(P)$ is bounded if $\omega$ is bounded. Now by part 4 of the same lemma, all global minima must be localized in a compact set. (In particular, $P\mapsto E(P)$ cannot be a constant function in an open set, hence everywhere). Now suppose there is a finite accumulation point $(Q_{*},E_{\mathrm{r}}(Q_{*}))$ of the set of global minima. Since the spectrum is closed, the corresponding sphere belongs to the spectrum. By Lemma 3.1 parts 2 and 3, for $|P|=Q_{*}$, $E(P)$ is an isolated eigenvalue of $H(P)$. Thus, by analyticity, this eigenvalue can be continued to some neighbourhood of $Q_{*}$ in the radial direction. As any neighbourhood contains infinitely many global minima, we conclude that $P\mapsto E(P)$ is constant near $|P|=Q_{*}$, hence everywhere, which is a contradiction.

Let us now move on to part (b). For the models in question we have $\Pi=F\mathrm{e}^{\mathrm{i}P_{\mathrm{f}}x}$, where $F$ is the unitary Fourier transform from $x$ to $P$ variable. We decompose $\varphi\otimes\Omega=\Pi^{*}\int^{\oplus}_{\mathbb{R}^{d}}\mathrm{d}P\phi_{P}$. Thus to find $\phi_{P}$ we compute

[TABLE]

which gives $\phi_{P}=\hat{\varphi}(P)\Omega$. The function $P\mapsto|\langle\phi_{P},\psi_{P}\rangle|=|\hat{\varphi}(P)|\,|\langle\Omega,\psi_{P}\rangle|$ is continuous by continuity of $\hat{\varphi}$ and analyticity of $P\mapsto|\psi_{P}\rangle\langle\psi_{P}|$. It is non-zero by our assumption on $\hat{\varphi}$ and by part 5 of Lemma 3.1.

Concerning the Fröhlich polaron, standing assumption 0 is proven in Appendix A. The other claims are verified as above, using Lemma 4.1.