Surface energy and boundary layers for a chain of atoms at low temperature

TL;DR

This paper investigates the low-temperature behavior of a chain of atoms with Lennard-Jones interactions, focusing on surface energy, boundary layers, and large deviations of Gibbs measures, revealing zero-temperature limits and correlation decay bounds.

Contribution

It establishes large deviation principles for Gibbs measures at low temperature, characterizes surface energy functionals, and links zero-temperature limits to boundary layers and free energy corrections.

Findings

Gibbs measures satisfy large deviations principles with specific energy functionals.

Surface correction converges to zero-temperature surface energy.

Bulk measures can be approximated by Gaussian measures.

Abstract

We analyze the surface energy and boundary layers for a chain of atoms at low temperature for an interaction potential of Lennard-Jones type. The pressure (stress) is assumed small but positive and bounded away from zero, while the temperature $\beta^{-1}$ goes to zero. Our main results are: (1) As $\beta \to \infty$ at fixed positive pressure $p>0$, the Gibbs measures $\mu_\beta$ and $\nu_\beta$ for infinite chains and semi-infinite chains satisfy path large deviations principles. The rate functions are bulk and surface energy functionals $\overline{\mathcal{E}}_{\mathrm{bulk}}$ and $\overline{\mathcal{E}}_\mathrm{surf}$. The minimizer of the surface functional corresponds to zero temperature boundary layers. (2) The surface correction to the Gibbs free energy converges to the zero temperature surface energy, characterized with the help of the minimum of $\overline{\mathcal{E}}_\mathrm{surf}$. (3) The bulk Gibbs measure and Gibbs free energy can be approximated by their Gaussian counterparts. (4) Bounds on the decay of correlations are provided, some of them uniform in $\beta$.

Equations

$${\mathcal{E}}_N(z_1,\dots ,z_{N-1})=\sum _{\genfrac{}{}{0ex}{}{1\le i<j\le N}{|i-j|\le m}}v(z_i+\cdots +z_{j-1})+p\sum _{j=1}^{N-1}{z}_j.$$

$$E_N=\underset{z_1,\dots ,z_{N-1}>0}{\inf }{\mathcal{E}}_N(z_1,\dots ,z_{N-1})$$

$$e_{0}=pa+\sum_{k=1}^{m}v(ka)=\min_{r>0}\Bigl{(}pr+\sum_{k=1}^{m}v(kr)\Bigr{)}.$$

$$h(z_1,\dots ,z_m)$$

$$\displaystyle\mathcal{E}_{\rm surf}\bigl{(}(z_{j})_{j\in\mathbb{N}}\bigr{)}$$

$$e_{\mathrm{surf}}=2\underset{\mathcal{D}}{\min }{\mathcal{E}}_{\mathrm{surf}}-pa-\sum _{k=1}^{m}kv(ka).$$

$$\displaystyle\mathcal{E}_{\mathrm{bulk}}\bigl{(}(z_{j})_{j\in\mathbb{Z}}\bigr{)}$$

$$\displaystyle=\sum_{j=-\infty}^{\infty}\sum_{k=1}^{m}\bigl{(}v(z_{j}+\cdots+z_{j+k-1})-v(ka)+\delta_{1k}p(z_{j}-a)\bigr{)}.$$

$${\mathcal{E}}_{\mathrm{bulk}}((z_j{)}_{j\in \mathbb{Z}})$$

$$+\mathcal{W}(\cdots z_{-1}{z}_0\mid z_1{z}_2\dots ),$$

$${\mathbb{Q}}_N^{(\beta )}(A)=\frac{1}{Q_N(\beta )}{\int }_A{e}^{-\beta {\mathcal{E}}_N(z_1,\dots ,z_{N-1})}\mathrm{d}{z}_1\cdots \mathrm{d}{z}_{N-1}$$

$$Q_N(\beta )={\int }_{{\mathbb{R}}_{+}^{N-1}}{e}^{-\beta {\mathcal{E}}_N(z_1,\dots ,z_{N-1})}\mathrm{d}{z}_1\cdots \mathrm{d}{z}_{N-1}.$$

$$\underset{N\to \infty }{\lim }{\int }_{{\mathbb{R}}_{+}^{N-1}}f(z_1,\dots ,z_k)\mathrm{d}{\mathbb{Q}}_N^{(\beta )}(z_1,\dots ,z_{N-1})={\int }_{{\mathbb{R}}_{+}^{\mathbb{N}}}f(z_1,\dots ,z_k)\mathrm{d}{\nu }_{\beta }((z_j{)}_{j\ge 1}).$$

$$\underset{N\to \infty }{\lim }{\int }_{{\mathbb{R}}_{+}^{N-1}}f(z_{i_N+1},\dots ,z_{i_N+k})\mathrm{d}{\mathbb{Q}}_N^{(\beta )}(z_1,\dots ,z_{N-1})={\int }_{{\mathbb{R}}_{+}^{\mathbb{Z}}}f(z_1,\dots ,z_k)\mathrm{d}{\mu }_{\beta }((z_j{)}_{j\ge 1}).$$

$$\overline{\mathcal{E}}_{\mathrm{surf}}\bigl{(}(z_{j})_{j\in\mathbb{N}}\bigr{)}=\begin{cases}\mathcal{E}_{\mathrm{surf}}\bigl{(}(z_{j})_{j\in\mathbb{N}}\bigr{)},&\quad(z_{j})_{j\in\mathbb{N}}\in\mathcal{D},\\ \infty,&\quad\text{else}.\end{cases}$$

$$\underset{\beta \to \infty }{\mathrm{limsup}}\frac{1}{\beta }\log {\nu }_{\beta }(A)$$

$$\underset{\beta \to \infty }{\mathrm{liminf}}\frac{1}{\beta }\log {\nu }_{\beta }(O)$$

$$-\frac{1}{\beta}\log\mu_{\beta}\bigl{(}{\operatorname{e}}^{-\beta\mathcal{V}}\bigr{)}\to\inf_{\mathcal{D}}(\mathcal{E}_{\mathrm{bulk}}+\mathcal{V})\quad(\beta\to\infty).$$

$$g(\beta)=-\lim_{N\to\infty}\frac{1}{\beta N}\log Q_{N}(\beta),\quad g_{\mathrm{surf}}(\beta)=\lim_{N\to\infty}\Bigl{(}-\frac{1}{\beta}\log Q_{N}(\beta)-Ng(\beta)\Bigr{)},$$

$$\underset{\beta \to \infty }{\lim }g(\beta )=e_0,\underset{\beta \to \infty }{\lim }{g}_{\mathrm{surf}}(\beta )=e_{\mathrm{surf}}.$$

$$\ell (\beta )={\int }_{{\mathbb{R}}_{+}^{\mathbb{Z}}}{z}_0\mathrm{d}{\mu }_{\beta }((z_j{)}_{j\in \mathbb{Z}}).$$

$$\mu_{\beta}\bigl{(}\{(z_{j})_{j\in\mathbb{Z}}\mid z_{0}\geq r\}\bigr{)}\leq C\exp(-\beta pr)$$

$$\lim_{\beta\to\infty}\ell(\beta)=a=\mathrm{argmin}\bigl{(}pr+\sum_{k=1}^{m}v(kr)\bigr{)}.$$

$$g(\beta)=-\frac{1}{\beta}\log\Bigl{(}\int_{0}^{\infty}{\operatorname{e}}^{-\beta[v(r)+p_{\beta}r]}\mathrm{d}r\Bigr{)},\quad\ell(\beta)=\frac{\int_{0}^{\infty}r\exp(-\beta[v(r)+p_{\beta}r])\mathrm{d}r}{\int_{0}^{\infty}\exp(-\beta[v(r)+p_{\beta}r])\mathrm{d}r}.\\$$

$$\underset{\beta \to \infty }{\lim }g(\beta )=\min (e_0^0,-\nu ).$$

$$\int_{\mathbb{R}^{\mathbb{Z}}}s_{i}s_{j}\mathrm{d}\mu^{\mathrm{Gauss}}\bigl{(}(s_{k})_{k\in\mathbb{Z}}\bigr{)}=(\mathcal{H}^{-1})_{ij}$$

$$\lim_{\beta\to\infty}\int_{\mathbb{R}^{n}}\Bigl{|}\beta^{-n/2}\rho_{n}^{(\beta)}\bigl{(}a+\beta^{-1/2}s_{1},\ldots,a+\beta^{-1/2}s_{n}\bigr{)}-\rho_{n}^{\mathrm{Gauss}}(s_{1},\ldots,s_{n})\Bigr{|}\mathrm{d}s_{1}\ldots\mathrm{d}s_{n}=0.$$

$$\lim_{\beta\to\infty}\int_{\mathbb{R}^{\mathbb{Z}}}f\bigl{(}\sqrt{\beta}(z_{j}-a)_{j\in\mathbb{Z}}\bigr{)}\mathrm{d}\mu_{\beta}\bigl{(}(z_{j})_{j\in\mathbb{Z}}\bigr{)}=\int_{\mathbb{R}^{\mathbb{Z}}}f\mathrm{d}\mu^{\mathrm{Gauss}}.$$

$$g(\beta )=e_0-\frac{1}{\beta }\log \sqrt{\frac{2\pi }{\beta (\det C{)}^{1/d}}}+o({\beta }^{-1})$$

$$g(\beta )=e_0-\frac{1}{\beta }\log \sqrt{\frac{2\pi }{\beta c}}+\sum _{j=1}^n{a}_j{\beta }^{-j/2}+O({\beta }^{-(n+1)/2})(\beta \to \infty )$$

$${\mathcal{I}}_{\mathrm{eff}}^{(\beta )}(k)=-\frac{1}{\beta }\log {\mu }_{\beta }(e^{-\beta ({\mathcal{V}}_0+{\mathcal{V}}_k)})+\frac{1}{\beta }\log {\mu }_{\beta }(e^{-\beta {\mathcal{V}}_0})+\frac{1}{\beta }\log {\mu }_{\beta }(e^{-\beta {\mathcal{V}}_k})$$

$${\mathcal{I}}_{\mathrm{eff}}^{(\infty )}(k)=\inf ({\mathcal{E}}_{\mathrm{bulk}}+{\mathcal{V}}_0+{\mathcal{V}}_k)-\inf ({\mathcal{E}}_{\mathrm{bulk}}+{\mathcal{V}}_0)-\inf ({\mathcal{E}}_{\mathrm{bulk}}+{\mathcal{V}}_k),$$

Full text

Keywords: atomistic models of elasticity - surface energy and boundary layers - semi-classical limit of transfer operators - uniform decay of correlations - path large deviations for stationary processes.

MSC2010 classification: 82B21, 74B20, 74G65, 60F10.

1. 1 Introduction

2. 2 Main results

3. 2.1 Zero temperature

4. 2.2 Small positive temperature

5. 2.3 Gaussian approximation

6. 2.4 Decay of correlations

7. 3 Energy estimates

8. 3.1 Bulk periodicity

9. 3.2 Surface energy

10. 3.3 A fixed point equation

11. 4 Gibbs measures for the infinite and semi-infinite chains

12. 4.1 Transfer operator

13. 4.2 Ergodicity

14. 4.3 Thermodynamic limit

15. 5 Large deviations as $\beta\to\infty$

16. 5.1 A tightness estimate

17. 5.2 Gibbs free energy in the bulk

18. 5.3 Large deviations principles for $\nu_{\beta}$ and $\mu_{\beta}$

19. 5.4 Surface corrections to the Gibbs free energy

20. 6 Gaussian approximation

21. 6.1 Decomposition of the energy. Choice of transfer operator

22. 6.2 Some properties of the transfer operator

23. 6.3 Gaussian transfer operator

24. 6.4 Perturbation theory

25. 6.5 Proof of Theorems 2.7, 2.8 and  2.11

26. 7 A Brascamp-Lieb type covariance estimate for $m=\infty$

1. Introduction

The purpose of the present article is to analyze the low-temperature behavior for a one-dimensional chain of atoms that interact via a Lennard-Jones type potential. The model is atomistic and in terms of the Gibbs measures of classical statistical mechanics. Two limiting procedures are at play: the zero-temperature limit, for which the inverse temperature $\beta$ goes to infinity, and the thermodynamic limit, where the number of particles $N$ and the system size go to infinity. The order of the limits matters. When the zero-temperature limit is taken before the $N\to\infty$ limit, the analysis of Gibbs measures is replaced by energy minimization, leading to variational models of non-linear elasticity. We perform instead the zero-temperature limit after the thermodynamic limit. The zero-temperature limit for infinite systems is far from trivial, see [vER07, CGU11, CH10] and the discussion in [BRS10].

For the one-dimensional Lennard-Jones interaction, it is known that energy minimizers (ground states) converge to a periodic lattice [GR79] (“crystallization”). For one-dimensional systems with pair potentials that decay faster than $1/r^{2}$ it is well-known that, in contrast, at positive temperature, no matter how small, there is no crystallization [BL15]. Nevertheless, some quantities can be approximated well by their zero-temperature counterpart. For the bulk free energy this is to be expected, for other quantities such as surface corrections this is already more subtle. For the decay of correlations, it is a priori not even clear what the zero-temperature counterpart should be; we propose a natural candidate, see Eqs (2.10) and (2.11).

At zero temperature, surface corrections and boundary layers have been studied, for example, in order to better understand variational models of fracture, see e.g. [BC07, SSZ11] and the references therein. Fracture might be expected for elongated chains, forced to stretch beyond their preferred length. At small positive temperature, large interparticle distances correspond to low pressure (stress) $p=p_{\beta}\to 0$. We address this regime in a subsequent work and focus here on the elastic regime of positive pressure $p>0$, though the case of small pressure $p_{\beta}\to 0$ is discussed in some comments.

Our main results come in four parts. They are listed in Sections 2.1–2.4 and proven in Sections 3–7. At zero temperature, we extend the result on bulk periodicity from [GR79] to a more general class of potentials and positive pressure, see Theorem 2.1. We prove the existence of bounded surface corrections, and characterize them with the help of an energy functional $\mathcal{E}_{\mathrm{surf}}$ for semi-infinite chains (Theorem 2.2).

At positive temperature, we prove large deviations principles for the Gibbs measures $\mu_{\beta}$ and $\nu_{\beta}$ on $\mathbb{R}_{+}^{\mathbb{Z}}$ and $\mathbb{R}_{+}^{\mathbb{N}}$ (product topology) as $\beta\to\infty$ at fixed $p>0$ (Theorem 2.4). The speed is $\beta$ and the respective rate functions are energy functionals $\overline{\mathcal{E}}_{\mathrm{bulk}}$ and $\overline{\mathcal{E}}_{\mathrm{surf}}-\min\overline{\mathcal{E}}_{\mathrm{surf}}$ whose minimizers are, respectively, the periodic bulk ground state and the zero-temperature boundary layer. The convergence of positive-temperature surface corrections to their zero-temperature counterpart is addressed in Theorem 2.5. These results are intimately related to path large deviations for Markov processes and Hamilton-Jacobi-Bellman equations [FK06], semi-classical analysis [Hel02], and a more direct approach to low-temperature expansions [SL17]. We remark that our results are valid for long range interactions which in particular are not assumed to have superlinear growth at infinity. The large deviations principle is complemented by a result on Gaussian approximations for the bulk Gibbs measure and the Gibbs free energy, valid for finite interaction range $m$ (Theorems 2.7 and 2.8).

Finally we study the temperature-dependence of correlations and informally discuss how correlations connect with effective interactions of defects and the decay of boundary layers. Theorem 2.9 provides a priori estimates that hold for all $\beta,p>0$. In Theorem 2.11 we show that for finite $m$ and small positive pressure $p$, the decay of correlations is exponential with a rate of decay that stays bounded as $\beta\to\infty$—the associated Markov chain has a spectral gap bounded away from zero. This uniform estimate is proven with perturbation theory for the transfer operator. For infinite $m$, we provide instead a uniform estimate for restricted Gibbs measures (Proposition 2.10), which follows from the convexity of the energy (in a neighborhood of the periodic gound state) and techniques from the realm of Brascamp-Lieb inequalities [Hel02]. At vanishing pressure $p_{\beta}\to 0$ or fixed high pressure $p>0$, the spectral gap might become exponentially small because of fracture or metastable wells [BdH15] in non-convex energy landscapes.

Bringing statistical mechanics into atomistic models of crystals and elasticity has a rich tradition [BH98, Wei02, BCF86, Pen02]. Modern developments include: the study of gradient Gibbs measures [FS97] with sophisticated tools such as renormalization groups and cluster expansions [AKM16], random walk representations [BFS82], and Witten Laplacians [Hel02]; scaling limits and gradient Young-Gibbs measures [Pre09, KL14, Run15]; the extension of approximation schemes, e.g., the quasi-continuum method, to positive temperature [BLBLP10, TM11]. In addition, there have been some inroads into the open problem of proving crystallization in the form of orientational order for two-dimensional models [Aum15, HMR14].

To the best of our knowledge, all of the aforementioned mathematical literature, notably on Gibbs gradient measures, is limited to potentials with a superlinear growth at infinity. This is in stark contrast with the decay to zero typically imposed in statistical mechanics of point particles [Rue69]. We work with potentials $v(r)\to 0$, an additional linear term $pr$ enters because we work in the constant pressure ensemble, which is the most convenient ensemble for one-dimensional systems [Rue69, Section 5.6.6]. As a consequence, the by now classical combination of Bakry-Émery estimates and Holley-Stroock perturbation principle, see [Men14] and the references therein, becomes potentially more delicate. We use instead estimates on energy penalties, some aspects of which might generalize to higher-dimensional models.

Another aspect that might generalize to higher dimension concerns the large deviations principle. The existence of a large deviations principle for the Gibbs measure as $\beta\to\infty$, proven using a exponential tightness and fixed point equation for the measure, amounts to the construction of an infinite volume energy functional that vanishes on ground states only. In higher dimension, the role of the fixed point equation is taken by DLR-conditions named after Dobrushin, Lanford, Ruelle [Geo11] and the proof of a large deviations principle reduces to the investigation of a higher-dimensional analogue of a Bellman equation. The theory of the latter, for non-unique ground states, might mirror possible intricacies of the zero-temperature limit of Gibbs measure described in [vER07].

Finally we remark that the results of this work allow for a detailed analysis of typical atomic configurations at low temperature and low density. In [JKST19] we will in particular prove that, when the density is strictly smaller than the density of the ground state lattice, a system with $N$ particles fills space by alternating approximately crystalline domains (“clusters”) with empty domains (“cracks”). The number of domains is of the order of $N\exp(-\beta e_{\rm surf}/2)$ with $e_{\rm surf}$ the surface energy from Theorem 2.2 below.

2. Main results

2.1. Zero temperature

Let $v:(0,\infty)\to\mathbb{R}$ be a pair potential, $m\in\mathbb{N}\cup\{\infty\}$ a truncation parameter and $p\geq 0$ the pressure. At zero temperature we allow for $p=0$, at positive temperature we impose $p>0$. The Gibbs energy at zero temperature and pressure $p$ for a system of $N$ particles with positions $x_{1}<\ldots<x_{N}$ and interparticle spacings $z_{j}=x_{j+1}-x_{j}$, $j=1,\ldots,N-1$, is

[TABLE]

The parameter $m$ restricts the range of the interaction: $m=2$ corresponds to a next-nearest neighbor interaction. This section deals with the minimization problem

[TABLE]

in the limit $N\to\infty$. Throughout we assume that the following assumption holds.

Assumption 1**.**

The pair potential $v:(0,\infty)\to\mathbb{R}\cup\{+\infty\}$ is equal to $+\infty$ on $(0,r_{\rm hc}]$ for some $r_{\rm hc}\geq 0$ and a $C^{2}$ function on $(r_{\rm hc},\infty)$. There exist $r_{\rm hc}<z_{\min}<z_{\max}<2z_{\min}$ and $\alpha_{1},\alpha_{2}>0$, $s>2$ such that the following holds.

• (i)

Shape of $v$: $z_{\max}$ is the unique minimizer of $v$ and satisfies $v(z_{\max})<0$. $v$ is decreasing on $(0,z_{\max})$ and increasing and non-positive on $(z_{\max},\infty)$.

• (ii)

Growth of $v$: $v(z)\geq-\alpha_{1}z^{-s}$ for all $z>0$ and $v(z)+v(z_{\max})-2\alpha_{1}\sum_{n=2}^{\infty}(nz)^{-s}>0$ for all $z<z_{\min}$.

• (iii)

Shape of $v^{\prime\prime}$: $v^{\prime\prime}$ is decreasing on $[z_{\min},z_{\max}]$ and increasing and non-positive on $[2z_{\min},\infty)$.

• (iv)

Growth of $v^{\prime\prime}$: $v^{\prime\prime}(z)\geq-\alpha_{2}z^{-s-2}$ for all $z>r_{\rm hc}$ and $v^{\prime\prime}(z_{\max})+\sum_{n=2}^{\infty}n^{2}v^{\prime\prime}(nz_{\min})>0$.

The assumption is satisfied, for example, by the Lennard-Jones potential $v(r)=r^{-12}-r^{-6}$. As we will see, parts (i) and (ii) of the assumption guarantee that energy minimizers at $p=0$ have interparticle spacings $z_{j}$ in $(z_{\min},z_{\max})$, parts (iii) and (iv) ensure that $\mathcal{E}_{N}$ is uniformly strictly convex in $(z_{\min},z_{\max})^{N-1}$; moreover the Hessian $\mathrm{D}^{2}\mathcal{E}_{N}$ is diagonally dominant with positive diagonal entries and negative off-diagonal entries.

Assumption 2**.**

The pressure $p$ satisfies $0\leq p<p^{*}$ with $p^{*}:=\frac{|v(z_{\max})|}{z_{\max}}$.

At positive temperature we shall assume in addition that $p>0$, $r_{hc}>0$, and for some results we need $\lim_{r\searrow r_{\mathrm{hc}}}v(r)=\infty$. The next theorem is the adaptation of a similar result by Gardner and Radin [GR79]. It is proven in Section 3.1.

Theorem 2.1** (Bulk properties).**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p\in[0,p^{*})$ as in Assumption 2.

1. (a)

For every $N\geq 2$, the map $\mathcal{E}_{N}:\mathbb{R}_{+}^{N-1}\to\mathbb{R}$ has a unique minimizer $(z_{1}^{\scriptscriptstyle{({N}})},\ldots,z_{N-1}^{\scriptscriptstyle{({N}})})$. The mimizer has all its spacings $z_{j}$ in $[z_{\min},z_{\max}]$. 2. (b)

As $j,N\to\infty$ along $N-j\to\infty$, we have $z_{j}^{\scriptscriptstyle{({N}})}\to a$ where $a\in(z_{\min},z_{\max}]$ is the unique minimizer of $\mathbb{R}_{+}\ni r\mapsto pr+\sum_{k=1}^{m}v(kr)$. 3. (c)

The limit $e_{0}=\lim_{N\to\infty}(E_{N}/N)<0$ exists and is given by

[TABLE]

Let $\mathcal{D}_{0}\subset(r_{\rm hc},\infty)^{\mathbb{N}}$ be the space of sequences $(z_{j})_{j\in\mathbb{N}}$ with none or at most finitely many elements different from $a$. Define

[TABLE]

When $m=\infty$, $h((z_{j})_{j\in\mathbb{N}})$ is a function of the whole sequence. $\mathcal{E}_{\mathrm{surf}}$ is the Gibbs energy of a semi-infinite chain, with additive constant chosen in such a way that at spacings $z_{j}\equiv a$ the Gibbs energy is zero; $h(z_{1},z_{2},\ldots)$ represents the interaction of the left-most particle with everybody else. Let $\mathcal{D}=\{(z_{j})_{j\in\mathbb{N}}\in(r_{\rm hc},\infty)^{\mathbb{N}}\mid\sum_{j=1}^{\infty}(z_{j}-a)^{2}<\infty\}$ be the space of square summable strains.

Theorem 2.2** (Surface energy).**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p\in[0,p^{*})$ as in Assumption 2. Equip $\mathcal{D}$ with the $\ell^{2}$-metric. Then

1. (a)

$\mathcal{E}_{\mathrm{surf}}$* extends to a continuous functional on $\mathcal{D}$.* 2. (b)

On $\mathcal{D}\cap[z_{\min},z_{\max}]^{\mathbb{N}}$ it is strictly convex. 3. (c)

$\mathcal{E}_{\mathrm{surf}}$* has a unique minimizer. The minimizer lies in $\mathcal{D}\cap[z_{\min},z_{\max}]^{\mathbb{N}}$.* 4. (d)

The limit $e_{\mathrm{surf}}=\lim_{N\to\infty}(E_{N}-Ne_{0})$ exists and is given by

[TABLE]

The theorem is proven in Section 3.2. Note that $-pa-\sum_{k=1}^{\infty}kv(ka)$ is the surface energy for a clamped chain with all spacings equal to $a$ and $\mathcal{E}_{\mathrm{surf}}$ encodes the effect of boundary layers. $\mathcal{E}_{\mathrm{surf}}$ is multiplied by $2$ because finite chains have two ends. We note that $\min\mathcal{E}_{\mathrm{surf}}$ is exactly the boundary layer energy introduced by Braides and Cicalese [BC07]; Braides and Cicalese dealt with the special case $m=2$ of next-nearest neighbor interactions but more general potentials. For finite $m\geq 2$, see  [SS18, Theorem 4.2].

For later purpose we also define a bulk functional

[TABLE]

It is defined, a priori, on the space $\mathcal{D}_{0}^{+}$ of positive bi-infinite sequences $(z_{j})_{j\in\mathbb{Z}}\in(r_{\rm hc},\infty)^{\mathbb{Z}}$ that have at most finitely many elements $z_{j}\neq a$. Denoting the space of square summable strains $\mathcal{D}^{+}=\{(z_{j})_{j\in\mathbb{Z}}\in(r_{\rm hc},\infty)^{\mathbb{Z}}\mid\sum_{j\in\mathbb{Z}}(z_{j}-a)^{2}<\infty\}$, an analysis similar to the one for the surface functional yields the following result.

Proposition 2.3** (Limiting bulk properties).**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p\in[0,p^{*})$ as in Assumption 2. Equip $\mathcal{D^{+}}$ with the $\ell^{2}$-metric. Then

1. (a)

$\mathcal{E}_{\mathrm{bulk}}$* extends to a continuous functional on $\mathcal{D}^{+}$.* 2. (b)

On $\mathcal{D}^{+}\cap[z_{\min},z_{\max}]^{\mathbb{N}}$ it is strictly convex. 3. (c)

The unique minimizer of $\mathcal{E}_{\mathrm{bulk}}$ is the constant sequence $(\ldots,a,a,\ldots)$. The minimum value is $\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)=0$. 4. (d)

For every $(z_{j})_{j\in\mathbb{Z}}\in\mathcal{D}^{+}$ one has

[TABLE]

where $\mathcal{W}(\cdots z_{-1}z_{0}\mid z_{1}z_{2}\ldots):=\sum_{\genfrac{}{}{0.0pt}{}{j\leq 0,k\geq 1}{|k-j|\leq m-1}}v(z_{j}+\cdots+z_{k})$ is the total interaction between the left and right half-infinite chain.

2.2. Small positive temperature

Next we analyze infinite volume Gibbs measures on $\mathbb{R}_{+}^{\mathbb{N}}$ and $\mathbb{R}_{+}^{\mathbb{Z}}$ in the limit $\beta\to\infty$. We focus on fixed positive $p\in(0,|v(z_{\max})|/z_{\mathrm{\max}})$ but comment on vanishing $p=p_{\beta}\to 0$ at the end of the section. Let $\mathbb{Q}_{N}^{\scriptscriptstyle{({\beta}})}$ be the probability measure on $\mathbb{R}_{+}^{N-1}$ defined by

[TABLE]

where

[TABLE]

Standard arguments (see Section 4) show there is a uniquely defined probability measure $\nu_{\beta}$ on the product space $\mathbb{R}_{+}^{\mathbb{N}}$ such that for every $k\in\mathbb{N}$, every bounded continuous test function $f\in C_{b}(\mathbb{R}_{+}^{k})$,

[TABLE]

Similarly, there is a uniquely defined probabilty measure $\mu_{\beta}$ on $\mathbb{R}_{+}^{\mathbb{Z}}$ such that for all local test functions $f$ as above, and all sequences $i_{N}$ with $i_{N}\to\infty$ and $N-i_{N}\to\infty$,

[TABLE]

Moreover the measure $\mu_{\beta}$ is shift-invariant and mixing. The measure $\mu_{\beta}$ describes the bulk behavior of a semi-infinite chain, the measure $\nu_{\beta}$ is the equilibrium measure for a semi-infinite chain and encodes the probability distribution of boundary layers.

Our first result is a large deviations principle for the equilibrium measure $\nu_{\beta}$ as $\beta\to\infty$. The rate function is a suitable extension of $\mathcal{E}_{\mathrm{surf}}$: define $\overline{\mathcal{E}}_{\mathrm{surf}}:\ \mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}\cup\{\infty\}$ by

[TABLE]

In the same way $\mathcal{E}_{\mathrm{bulk}}$ extends to a map $\overline{\mathcal{E}}_{\mathrm{bulk}}$ from $\mathbb{R}_{+}^{\mathbb{Z}}$ to $\mathbb{R}\cup\{\infty\}$. Both $\mathbb{R}_{+}^{\mathbb{N}}$ and $\mathbb{R}_{+}^{\mathbb{Z}}$ are equipped with the product topology.

Theorem 2.4**.**

Fix $p\in(0,p^{*})$ and $m\in\mathbb{N}\cup\{\infty\}$. Assume that $r_{\mathrm{hc}}>0$ and $\lim_{r\searrow r_{\mathrm{hc}}}v(r)=\infty$. Then as $\beta\to\infty$, the equilibrium measures $(\nu_{\beta})_{\beta>0}$ and $(\mu_{\beta})_{\beta>0}$ satisfy large deviations principles with speed $\beta$ and respective rate functions $\overline{\mathcal{E}}_{\mathrm{surf}}-\min\mathcal{E}_{\mathrm{surf}}$ and $\overline{\mathcal{E}}_{\mathrm{bulk}}$. The rate functions are good, i.e., lower semi-continuous with compact level sets.

The theorem is proven in Section 5.3. The large deviations principle for $\nu_{\beta}$ says that for every closed set $A\subset\mathbb{R}_{+}^{\mathbb{N}}$ and every open set $O\subset\mathbb{R}_{+}^{\mathbb{N}}$ (product topology)

[TABLE]

It is essential that we work in the product topology. Indeed we shall later see that $\nu_{\beta}$ is mixing, therefore for every $\varepsilon>0$, the measure $\nu_{\beta}$ gives full mass $1$ to sequences $(z_{j})_{j\in\mathbb{N}}$ that have infinitely many spacings $|z_{j}-a|>\varepsilon$. Thus for every ball $O=\{(z_{j})\in\mathbb{R}_{+}^{\mathbb{N}}\mid\sum_{j=1}^{\infty}(z_{j}-a)^{2}<\delta\}$, we have $\nu_{\beta}(O)=0$ hence $\beta^{-1}\log\nu_{\beta}(O)=-\infty$, to be contrasted with the lower bound in Eq. (2.5).

Another consequence concerns the evaluation of the Gibbs energies of localized defects: suppose that because of some impurity, the energy is not $\mathcal{E}_{N}$ but $\mathcal{E}_{N}+\mathcal{V}$, where $\mathcal{V}$ is, say, continuous in the product topology, localized in the bulk, and bounded from below. Then by Varadhan’s lemma [DZ98], as $\beta\to\infty$, the effective Gibbs energy converges to the zero temperature energy of the defect,

[TABLE]

Surface energies occur as a specific type of defect, when $\mathcal{V}$ cancels all interactions between two half-infinite chains (see Proposition 4.9(a)), which leads to the following theorem. Define

[TABLE]

the Gibbs free energy $g(\beta)$ per particle in the bulk and the surface correction $g_{\mathrm{surf}}(\beta)$.

Theorem 2.5**.**

Fix $p\in(0,p^{*})$ and $m\in\mathbb{N}\cup\{\infty\}$. The limits (2.6) exist. If in addition $r_{\mathrm{hc}}>0$ and $\lim_{r\searrow r_{\mathrm{hc}}}v(r)=\infty$, then the bulk and surface Gibbs energy approach their zero-temperature counterparts when $\beta\to\infty$:

[TABLE]

This proves that the thermodynamic limit and the zero temperature limit can be exchanged, which is non-trivial (and in fact, fails when the pressure goes to zero too fast, see below).

One last consequence of Theorem 2.4 concerns the distribution of spacings and the pressure-density (or stress-strain) relation. The Gibbs free energy and our partition functions correspond to an ensemble where the overall length of the system is not fixed, but instead may fluctuate with a law that depends on the pressure—high pressures $p$ favor compressed states. In the thermodynamic limit $N\to\infty$, though, the average spacing between particles becomes a well-defined quantity, given by

[TABLE]

By the contraction principle [DZ98, Theorem 4.2.1], the distribution of $z_{0}$ under $\mu_{\beta}$ satisfies a large deviations principle with good rate function $w(z)=\inf\{\overline{\mathcal{E}}_{\mathrm{bulk}}((z_{j})_{j\in\mathbb{Z}})\mid(z_{j})_{j\in\mathbb{Z}}\in\mathbb{R}_{+}^{\mathbb{Z}},\,z_{0}=z\}$. The unique minimizer of $w(z)$ is the ground state spacing $a$. Lemma 5.1 implies that the distribution of spacings has exponential tails

[TABLE]

for some $\beta$-independent constant $C$.

Corollary 2.6**.**

Under the assumptions of Theorem 2.5, we have

[TABLE]

In particular, for large $\beta$, we have $\ell(\beta)<a_{0}$ where $a_{0}$ is the minimizer of the zero-stress Cauchy-Born energy density $\sum_{k}v(kr)$. Conversely, spacings $\ell(\beta)>a_{0}$ (elongated chains) imply vanishing pressure $p=p_{\beta}\to 0$. This is clearly apparent for nearest neighbor interactions ($m=1$, Takahashi nearest neighbor gas [Tak42, LM66]), for which

[TABLE]

Comments on vanishing pressure. We add a superscript to indicate that zero-temperature quantities are evaluated at $p=0$. When $p=p_{\beta}\to 0$ slower than any exponential, it is still true that $g(\beta)\to e_{0}^{0}$. When $\beta p_{\beta}=\exp(-\beta\nu)$ with $\nu>0$, one can show with [JKM15, Jan12] that

[TABLE]

At pressures vanishing faster than $\exp(-\beta|e_{0}^{0}|)$, the most likely configurations have very large spacings (dilute gas phase, $\ell(\beta)\to\infty$) and the previous results no longer apply. For $\liminf\frac{1}{\beta}\log(\beta p_{\beta})>e_{0}^{0}$, we expect that large deviations principles with rate functions $\overline{\mathcal{E}}_{\mathrm{bulk}}^{0}$ and $\overline{\mathcal{E}}_{\mathrm{surf}}^{0}-\min\overline{\mathcal{E}}_{\mathrm{surf}}^{0}$ still hold (in fact our proofs still show weak large deviations principles). However rate functions have non-compact level sets and exponential tightness is lost. Moreover large spacings may contribute to the average (2.7) and Corollary 2.6 need no longer be true, thus allowing for spacings $\ell(\beta)\to\ell>a_{0}$.

2.3. Gaussian approximation

Here we complement the large deviations result by a Gaussian approximation. This section deals with finite $m$ and the bulk measure $\mu_{\beta}$ only. Remember $d=m-1$. We will see that the Hessian of $\mathcal{E}_{\mathrm{bulk}}$ at $(\ldots,a,a,\ldots)$ is associated with a positive-definite, bounded operator $\mathcal{H}$ in $\ell^{2}(\mathbb{Z})$. It is represented by a doubly-infinite matrix $(\mathcal{H}_{ij})_{i,j\in\mathbb{Z}}$ that is diagonally dominant. Write $(\mathcal{H}^{-1})_{ij}$ for the matrix elements of the inverse operator and let $\mu^{\mathrm{Gauss}}$ be the uniquely defined measure on $\mathbb{R}^{\mathbb{Z}}$, equipped with the product topology and its associated Borel $\sigma$-algebra, such that

[TABLE]

for all $i,j\in\mathbb{Z}$, and every finite-dimensional marginal of $\mu^{\mathrm{Gauss}}$ is a multi-dimensional Gaussian distribution. Equivalently, $\mu^{\mathrm{Gauss}}$ is the distribution of a Gaussian process $(N_{j})_{j\in\mathbb{Z}}$ with mean zero and covariance $\mathbb{E}[N_{i}N_{j}]=(\mathcal{H}^{-1})_{ij}$. More concrete expressions for the probability density functions of $nd$-dimensional marginals of $\mu^{\mathrm{Gauss}}$ are provided in Proposition 6.17 below.

In the following we identify the measure $\mu_{\beta}$ on $\mathbb{R}_{+}^{\mathbb{Z}}$ with the measure ${\mathchoice{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.5mu\mathrm{l}}{1\mskip-5.0mu\mathrm{l}}}_{\mathbb{R}_{+}^{\mathbb{Z}}}\mu_{\beta}$ on $\mathbb{R}^{\mathbb{Z}}$. We exclude the trivial case $m=1$.

Theorem 2.7**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and $r_{\mathrm{hc}}>0$. Then for every $n\in\mathbb{N}$, the $n$-dimensional marginals of $\mu_{\beta}$ and $\mu^{\mathrm{Gauss}}$ have probability density functions $\rho_{n}^{(\beta)}$ and $\rho_{n}^{\mathrm{Gauss}}$, and

[TABLE]

It follows that the distribution of the spacings, suitably rescaled, converges locally to the Gaussian measure $\mu^{\mathrm{Gauss}}$: for every bounded function $f:\mathbb{R}^{Z}\to\mathbb{R}$ that depends on finitely many spacings $z_{j}$ only (bounded cylinder functions), we have

[TABLE]

For example, in the limit $\beta\to\infty$, the distribution of a single spacing $z_{j}$ is approximately normal, with mean $a$ and variance $\beta^{-1}(\mathcal{H}^{-1})_{ii}$. We expect that Theorem 2.7 stays true for $m=\infty$ but a proof or disproof is beyond the scope of this article.

The next theorem says that the Gibbs free energy is close to the Gibbs free energy of the approximate Gaussian model.

Theorem 2.8**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and $r_{\mathrm{hc}}>0$. The Gibbs free energy satisfies, as $\beta\to\infty$,

[TABLE]

where $d=m-1$ and $C$ is a $d\times d$ positive-definite matrix.

The matrix $C$ is introduced in Eq. (6.18), see also Lemma 6.7, it is a function of the Hessian of the energy.

Remark* (Gaussian approximation and semi-classical expansions).*

If $v$ is smooth and $p>0$ is fixed, the Gibbs energy should admit an asymptotic expansion of form

[TABLE]

to arbitrarily high order $n$, for some $c>0$ and coefficients $a_{j}\in\mathbb{R}$. The first correction comes from a Gaussian approximation of the partition function (harmonic crystal), see Section 6, with the constant $c$ capturing the asymptotic behavior of the determinant of the Hessian around the energy minimum. Higher order corrections correspond to anharmonic effects. A similar expansion holds for $g_{\mathrm{surf}}(\beta)$. Rigorous results for finite $m$ are derived with semi-classical analysis [Hel02, Møl01, BM03] which build on the analogy with the $\hbar\to 0$ limit from quantum mechanics. For $m=2$ and potentials with superlinear growth at infinity, independent results are given in [SL17].

2.4. Decay of correlations

Suppose that two defects change the energy functional from $\mathcal{E}_{\mathrm{bulk}}$ to $\mathcal{E}_{\mathrm{bulk}}+\mathcal{V}_{0}+\mathcal{V}_{k}$, where we assume for simplicity that $\mathcal{V}_{0}$ and $\mathcal{V}_{k}$ depend on $z_{0}$ and $z_{k}$ alone. For large $k$, we may expect that the Gibbs energies are approximately additive, i.e.,

[TABLE]

should be small when the defects are far apart. $\mathcal{I}_{\mathrm{eff}}^{\scriptscriptstyle{({\beta}})}(k)$ represents an effective interaction between the defects. In the study of systems with many defects it is important to understand how fast the effective interaction decreases at large distances. Some intuition is gained from the zero-temperature counterpart

[TABLE]

however in general the limits $\beta,k\to\infty$ cannot be interchanged and a full study of (2.10) for large $k$ requires techniques beyond variational calculus.

A closely related problem is about the localization of changes induced by a defect: at zero temperature, if $(z_{j})_{j\in\mathbb{Z}}$ is a minimizer of $\mathcal{E}_{\mathrm{bulk}}+\mathcal{V}_{0}$, how fast does $z_{k}$ converge to the ground state spacing $a$ as $k\to\pm\infty$? On a similar note, how fast does $z_{k}\to a$ for a minimizer of the surface energy $\mathcal{E}_{\mathrm{surf}}$ (decay of boundary layers)? At positive temperature, the question is about the speed of convergence, for test functions $f:\mathbb{R}_{+}^{k}\to\mathbb{R}$, in

[TABLE]

as $i\to\infty$. Here $f_{i}((z_{j})_{j\in\mathbb{Z}}):=f(z_{i},\ldots,z_{i+k-1})$, so that $f_{n+i}=f_{i}\circ\tau^{n}$ when $\tau$ denotes the left shift on $\mathbb{R}_{+}^{\mathbb{Z}}$. These questions naturally lead to the investigation of the decay of correlations. We start with a general result which holds for all $\beta,p>0$.

Theorem 2.9**.**

Assume $m\in\mathbb{N}\cup\{\infty\}$ and $p>0$. There exist $c,C>0$ such that for all $\beta,p>0$, $k\in\mathbb{N}$, and bounded $f,g:\mathbb{R}_{+}^{k}\to\mathbb{R}$,

[TABLE]

When $m$ is finite and $k=m-1$, we have the stronger bound

[TABLE]

The theorem is proven in Section 4.2. When $m$ is finite, it implies exponential decay of correlations as $n\to\infty$, however the rate $-\log(1-{\operatorname{e}}^{-c\beta})$ can be exponentially small for large $\beta$. When $m$ is infinite, Theorem 2.9 implies algebraic decay of correlations: for $q=\lfloor n^{\varepsilon}\rfloor$ and sufficiently large $n$, $(1-{\operatorname{e}}^{-c\beta})^{q}$ is negligible compared to $\beta(q/n)^{s-2}$ and we find that as $n\to\infty$

[TABLE]

Better bounds are available for restricted Gibbs measures. Let $\tilde{\mu}_{\beta}^{\scriptscriptstyle{({N}})}$ be the measure $\mathbb{Q}_{N}^{\scriptscriptstyle{({\beta}})}$ conditioned on $[z_{\mathrm{min}},z_{\mathrm{max}}]^{N-1}$ and $\tilde{\mu}_{\beta}$ the probability measure on $[z_{\mathrm{min}},z_{\mathrm{max}}]^{\mathbb{Z}}$ obtained from the thermodynamic limit of $\tilde{\mu}_{\beta}^{\scriptscriptstyle{({N}})}$.

Proposition 2.10**.**

Let $m\in\mathbb{N}\cup\{\infty\}$. There exists $c>0$ such that for all $\beta,p>0$, smooth $f,g:\mathbb{R}_{+}\to\mathbb{R}$, and $i\neq j$,

[TABLE]

Remark*.*

When $m$ is finite, the uniform algebraic decay for the restricted Gibbs measure is replaced with uniform exponential decay $\exp(-\gamma|j-i|)$ with $\beta$-independent $\gamma>0$.

The proposition is proven in Section 7. It follows from the uniform convexity of the energy (Lemma 3.3) and known results from the realm of Brascamp-Lieb, Poincaré and Log-Sobolev inequalities. Proposition 2.10 differs from the estimate (2.12) in two ways: there is no exponentially large prefactor $\exp(c\beta)$, and the rate of algebraic decay is $1/n^{s}$ instead of $1/n^{s-2}$. Exponentially large prefactors are absent because the energy landscape has no local minimum. The improved algebraic decay $1/n^{s}$ arises, roughly, because the Gibbs measure is comparable to a Gaussian measure whose covariance is the inverse of the energy’s Hessian near the minimum, and instead of the tails of $v(r)$, it is the tails of $v^{\prime\prime}(r)$ that count.

We suspect that for large $\beta$ and small pressure, these improvements should carry over to the full Gibbs measure $\mu_{\beta}$, but we have proofs for interactions involving finitely many neighbors only.

Theorem 2.11**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and $r_{\mathrm{hc}}>0$. There exists $\gamma>0$ such that for all sufficiently large $\beta$, suitable $C(\beta)$, all $n\in\mathbb{N}$, and all $f,g:\mathbb{R}_{+}^{d}\to\mathbb{R}$, we have

[TABLE]

If $m=2$, we can pick $C(\beta)=1$.

The theorem is proven in Section 6 with perturbation theory for compact integral operators in $L^{2}(\mathbb{R}^{d})$. When $m=2$, the relevant operators are self-adjoint and spectral norms and operator norms coincide, leading to improved statements. We conclude with a few comments.

Lagrangian vs. Eulerian point of view. The theorems above formulate decay of correlations in terms of labelled spacings, which in the language of continuum mechanics is a Lagrangian viewpoint. On the other hand, in statistical mechanics of point particles it is more common to deal with unlabelled particles (Eulerian viewpoint) and correlations are between portions of space rather than labelled interparticle distances. The difference between the two approaches becomes quite clear for nearest neighbor interactions ($m=1$, see Eq. (2.8)), for which the spacings are i.i.d. with probability density $q_{\beta}(r)$ proportional to $\exp(-\beta[v(r)+p_{\beta}r])$. Because of the independence of spacings, correlations in terms of spacings vanish, $\mu_{\beta}(f_{0}g_{n})-\mu_{\beta}(f_{0})\mu_{\beta}(g_{n})=0$. On the other hand, the two-point function $\rho_{2}(0,x)$111Intuitively, $\rho_{2}(0,x)$ represents the probability for having one particle at [math] and one particle at $x$. Rigorously, $\rho_{2}(x_{1},x_{2})=\rho_{2}(0,x_{2}-x_{1})$ and for every $A$, $\int_{A}\rho_{2}(x_{1},x_{2})\mathrm{d}x_{1}\mathrm{d}x_{2}$ is the average number $\langle N_{A}(N_{A}-1)\rangle$ of ordered particle pairs in $A$. studied in statistical mechanics of particles is a sum over the number of particles contained in $(0,x]$,

[TABLE]

with $q_{\beta}^{*k}$ the $n$-fold convolution of $q_{\beta}$ with itself. It is a well-known fact from renewal theory [Fel71, Chapter XI] that

[TABLE]

but in general the difference is non-zero finite for $x$—in fact changing $q_{\beta}$ the convergence as $x\to\infty$ can be arbitrarily slow, even though correlations of labelled interparticle spacings vanish identically. One should keep this difference in mind when browsing the literature.

Path-large deviations, non-linear semi-groups, Bellman equation. For $m=2$, we may view $\mu_{\beta}$ as the law of a stationary Markov chain with state space $\mathbb{R}_{+}$ and transition kernel $P_{\beta}$ defined in Eq. (6.6). Theorem 2.4 is a path-large deviations result for the Markov chain. Path large deviations are often investigated with the help of non-linear semi-groups and Hamilton-Jacobi-Bellman equations [FK06]. In our context, a natural non-linear semi-group is

[TABLE]

and for sufficiently smooth $f$ we have a convergence of the form

[TABLE]

where $u$ solves

[TABLE]

Similar equations, motivated by quantum mechanics and geometric optics, appear in semi-classical analysis [Hel02, Eq. (5.4.4)]. Proposition 3.9 below provides an infinite-$m$ ersatz and is instrumental in the proof of Theorem 2.4.

Vanishing pressure. When $\beta p=\beta p_{\beta}\to 0$ faster than $\exp(-\beta|e_{0}^{0}|)$ (see (2.9)), the Gibbs measure should no longer be comparable to a Gaussian. Instead, it should be close to the ideal gas measure, for which spacings are i.i.d. exponentially distributed with parameter $\beta p_{\beta}$, and we may again expect uniform exponential decay of correlations (for finite $m$). When $\beta p_{\beta}\to 0$ at a speed comparable to $\exp(-\beta|e_{0}^{0}|)$, we should instead expect an exponentially small spectral gap: the Markov chain has two metastable wells, one corresponding to the optimal spacing $a$ and another well at infinity. The exponentially small spectral gap is associated with the fracture of the chain of atoms, in the spirit of “fracture as a phase transition” [Tru96].

3. Energy estimates

In this section we analyze the variational problems arising at zero temperature. Throughout the section we assume that $p\in[0,p^{*})$ as in Assumption 2.

3.1. Bulk periodicity

Lemma 3.1**.**

Every minimizer of $\mathcal{E}_{N}:\mathbb{R}_{+}^{N-1}\to\mathbb{R}$ lies in $[z_{\min},z_{\max}]^{N-1}$.

Proof.

Let $z_{1},\ldots,z_{N-1}>0$. If $z_{j}>z_{\max}$ for some $j$, define a new configuration by shrinking $z_{j}$ to $z_{\max}$, leaving all other spacings unchanged: $z^{\prime}_{i}=z_{i}$ for $i\neq j$ and $z^{\prime}_{j}=z_{\max}$. Since $z_{\rm max}$ is a strict minimizer of $v$ and $r\mapsto v(r)$ increases on $[z_{\max},\infty)$, shrinking the bonds decreases $\mathcal{E}_{N}$ strictly and the original configuration could not have been a minimizer.

If some interparticle spacing is smaller than $z_{\min}$, we remove a particle and reattach it to one end of the chain as follows. Assume $b:=\min(z_{1},\ldots,z_{N-1})<z_{\min}$ and let $j\in\{1,\ldots,N-1\}$ with $z_{j}=b$. Let $x_{1}=0$ and $x_{i}=z_{1}+\cdots+z_{i-1}$, $i=2,\ldots,N$ be associated particle positions. Thus $x_{j+1}-x_{j}=z_{j}=b$ and $x_{i+1}-x_{i}\geq b$ for all $i$. The interaction of $x_{j}$ with all other particles is

[TABLE]

For finite $m$ we note that, if $v(z_{j-i}+\ldots+z_{j-i+m})>0$ for an $i\in\{1,\ldots,\min\{m,j-1\}\}$, then $v(z_{j-i}+\ldots+z_{j-i+m})<v(z_{j-i}+\ldots+z_{j-1})$ by Assumption 1(i). Removing the particle $x_{j}$ thus leads to a configuration of $N$ atoms whose energy has decreased by at least

[TABLE]

The last inequality holds because of Assumption 1(ii) and $b<z_{\rm min}$. We define a new configuration by attaching the removed particle to either end of the chain at a distance $r=z_{\rm max}$. Since $v(z_{\max})+pz_{\max}<0$ by Assumption 2, this decreases $\mathcal{E}_{N}$ further, so overall the new configuration has strictly smaller energy, and the original sequence of spacings cannot be a minimizer of $\mathcal{E}_{N}$. ∎

At zero pressure, it is a well-known fact that the $N$-particle energy is subadditive, $E_{N+M}\leq E_{N}+E_{M}$. Indeed placing two $N$,$M$-particle minimizers side by side with large mutual distance, because of $v(r)\to 0$ at $r\to\infty$, yields an $N+M$-particle configuration with energy $\leq E_{N}+E_{M}$. Positive pressure penalizes large mutual distances between two consecutive blocks, so the construction has to be modified.

Lemma 3.2**.**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p\in[0,p^{*})$. Then $E_{N+M-1}\leq E_{N}+E_{M}$ for all $N,M\in\mathbb{N}$, and the limit $e_{0}=\lim E_{N}/N$ exists and satisfies $E_{N}\geq(N-1)e_{0}$ for all $N\in\mathbb{N}$.

Proof.

Let $z\in(r_{\mathrm{hc}},\infty)^{N-1}$ and $w\in(r_{\mathrm{hc}},\infty)^{M-1}$ be minimizers of $\mathcal{E}_{N}$ and $\mathcal{E}_{M}$ respectively. Define $y\in(r_{\mathrm{hc}},\infty)^{M+N-2}$ by concatenating $z$ and $w$. By Lemma 3.1, all spacings are in $[z_{\min},z_{\max}]$. Therefore interactions that involve bonds from both blocks are for spacings $\geq 2z_{\min}>z_{\max}$, hence negative, and

[TABLE]

As a consquence, $a_{n}:=E_{n+1}$ is subadditive. By Fekete’s subadditive lemma, the limit $e_{0}=\lim a_{n}/n=\lim E_{n}/n$ exists and is equal to the infimum of $a_{n}/n$, hence $E_{N}\geq(N-1)e_{0}$. Notice that $e_{0}>-\infty$ since

[TABLE]

(In the terminology of statistical mechanics, the energy is stable [Rue69, Chapter 3.2].) ∎

The next lemma in particular shows that $\mathcal{E}_{N}$ is uniformly convex on $[z_{\min},z_{\max}]^{N-1}$. For later purposes, we state and prove this on a slightly larger set.

Lemma 3.3**.**

There are constants $\varepsilon,\eta,C>0$ such that for all $m,N,N_{1},N_{2}\in\mathbb{N}$ with $N_{1}<N_{2}\leq N$, and $z=(z_{1},\ldots,z_{N-1})\in[z_{\min},\infty]^{N-1}$ with $z_{j}\leq z_{\max}+\varepsilon$ for $N_{1}\leq j\leq N_{2}-1$, the Hessian of $\mathcal{E}_{N}$ at $z$ satisfies

[TABLE]

for all $\zeta\in\mathbb{R}^{N-1}$. Moreover, the submatrix $(\partial_{i}\partial_{j}\mathcal{E}_{N}(z))_{N_{1}\leq i.j\leq N_{2}-1}$ of the Hessian has strictly positive diagonal entries $\partial_{i}^{2}\mathcal{E}_{N}(z)>0$ and non-positive off-diagonal entries $\partial_{i}\partial_{j}\mathcal{E}_{N}(z)\leq 0$. In particular, this matrix is monotone.

Note that the Hessian is independent of the pressure $p$.

Proof.

Let $\mathcal{L}$ be the collection of discrete intervals $\{i,\ldots,j-1\}\subset\{1,\ldots,N-1\}$ of length $j-i\leq m$. Then for all $i,j$

[TABLE]

For $i\neq j$ and $i,j\in L$ we have $\sum_{j\in L}z_{j}\geq 2z_{\min}$ hence $v^{\prime\prime}(\sum_{L}z_{j})\leq 0$; it follows that the off-diagonal entries of the Hessian are non-positive. Next we show that the row-sums are bounded from below by some constant $\eta>0$ if $N_{1}\leq i\leq N_{2}-1$.

[TABLE]

Assumption 1 guarantees that $\eta>0$ for $\varepsilon>0$ sufficiently small. Thus row sums are positive, off-diagonal matrix elements non-positive, and consequently diagonal elements positive. Moreover, with $C=2\max\{v^{\prime\prime}(r)\mid r\in[z_{\min},z_{\max}+\varepsilon]\}$ the diagonal elements are bounded from above by $\frac{C}{2}$. The proof of the lemma is then completed with the help of standard arguments, for example every eigenvalue of $(\partial_{i}\partial_{j}\mathcal{E}_{N}(z))_{N_{1}\leq i.j\leq N_{2}-1}$ lies in a Gershgorin circle with center $\partial_{i}^{2}\mathcal{E}_{N}$ and radius $\sum_{j\neq i}|\partial_{i}\partial_{j}\mathcal{E}_{N}|$. In particular, $(\partial_{i}\partial_{j}\mathcal{E}_{N}(z))_{N_{1}\leq i,j\leq N_{2}-1}$ is an M-matrix and thus monotone. ∎

Proof of Theorem 2.1.

(a) By Lemma 3.1 minimizers lie in the compact set $[z_{\min},z_{\max}]^{N-1}$. On that set the Hessian of $\mathcal{E}_{N}$ is positive definite because of Lemma 3.3, so $\mathcal{E}_{N}$ is strictly convex and the minimzer is unique.

(b) The convergence $z_{j}^{\scriptscriptstyle{({N}})}\to a$ as $j,N\to\infty$ along $N-j\to\infty$, where $a\in[z_{\min},z_{\max}]$ is the unique minimizer of $\mathbb{R}_{+}\ni r\mapsto pr+\sum_{k=1}^{m}v(kr)$, with the help of Lemma 3.3 is a straightforward adaptation of the corresponding proof in [GR79] and will be omitted. By Assumption 1(ii) we even have $a>z_{\min}$. We remark that the proof in [GR79] also shows that $\max\{z_{j}^{\scriptscriptstyle{({N+1}})},\ z_{j+1}^{\scriptscriptstyle{({N+1}})}\}\leq z_{j}^{\scriptscriptstyle{({N}})}$ for $j=1,\ldots,N-1$. This in turn implies that the convergence is in fact uniform away from a boundary layer of vanishing volume fraction.

(c) This observation in combination with Lemma 3.2 yields (c). Note that $e_{0}<0$ since $e_{0}\leq pz_{\max}+\sum_{k=1}^{\infty}v(kz_{\max})\leq pz_{\max}+v(z_{\max})<0$ by Assumptions 1 and 2.

∎

Notice that also $a<z_{\max}$ except for the exceptional cases in which only nearest neighbors interact, i.e. $m=1$ or $v(z)=0$ for $z\geq 2z_{\max}$, and the pressure vanishes.

3.2. Surface energy

Proposition 3.4**.**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p\geq 0$. Then

[TABLE]

Proof.

For simplicity we write down the proof for $m=\infty$; the proof when $m\in\mathbb{N}$ is completely analogous. Fix $k\geq 2$ and $\varepsilon>0$. Let $n_{1},n_{2}\in\mathbb{N}$ with $n_{2}\geq k$ and $N=n_{1}+n_{2}+1$. Let $z=(z_{-n_{1}},\ldots,z_{n_{2}-1})\in[z_{\min},z_{\max}]^{n_{1}+n_{2}}$ be the spacings of the $N$-particle ground state, labelled by $j=-n_{1},\ldots,n_{2}-1$ rather than $1,\ldots,N-1$. Choosing $n_{1}$ and $n_{2}$ large enough we may assume $\sum_{j=0}^{k-1}|z_{j}-a|^{2}\leq\varepsilon$. Since the Hessian has matrix norm uniformly bounded from above (Lemma 3.3), changing the spacings $z_{0},\ldots,z_{k-1}$ to $a$ increases the energy by $C\varepsilon$ at most thus

[TABLE]

We decompose the energy of the modified configuration as $A_{N}+B_{N}+C_{N}+D_{N}$ where

[TABLE]

where $\mathcal{W}$ gathers interactions that involve bonds from two consecutive blocks. The term $D_{N}$ represents the interactions between the left and right blocks. It satisfies

[TABLE]

which goes to zero as $k\to\infty$. Next we subtract $Ne_{0}$ from $\mathcal{E}_{N}$ and distribute it as $Ne_{0}=n_{1}e_{0}+(k+1)e_{0}+(n_{2}-k)e_{0}$ over the first three sums. The middle block contributes

[TABLE]

as $k\to\infty$. For the first block, we notice that

[TABLE]

Indeed the only missing piece are negative interactions between the left block and the right tail of a semi-infinite chain. The contribution of the right block $C_{N}$ is estimated in a similar way. We combine the estimates and let first $n_{1},n_{2}\to\infty$, then $k\to\infty$, and finally $\varepsilon\to 0$ and find

[TABLE]

For the upper bound, we take approximate minimizers of $\mathcal{E}_{\mathrm{surf}}$ and glue them together to an $N$-particle configuration by assigning them to the left and right boundaries, with spacings $a$ in between. This yields an $N$-particle configuration with energy $\mathcal{E}_{N}(z)-Ne_{0}\leq 2\inf_{\mathcal{D}_{0}}\mathcal{E}_{\mathrm{surf}}-\sum_{n=1}^{\infty}nv(na)+O(\varepsilon)$, and the required upper bound follows. ∎

Next we extend $\mathcal{E}_{\mathrm{surf}}$ to the space $\mathcal{D}\subset(r_{\rm hc},\infty)^{\mathbb{N}}$ of sequences with $\sum_{j=1}^{\infty}(z_{j}-a)^{2}<\infty$.

Lemma 3.5**.**

Let $m\in\mathbb{N}\cup\{\infty\}$. Let $\beta_{j}=\sum_{k=j+1}^{m}(k-j)v^{\prime}(ka)$, $j=1,\ldots,m-1$. Then for all $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$, we have

[TABLE]

The right-hand side is absolutely convergent for all $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}$.

Proof.

Let $\gamma_{j}=z_{j}-a$. Using $e_{0}=\sum_{k=1}^{m}v(ka)$, we have

[TABLE]

The equilibrium condition $p+\sum_{k=1}^{m}kv^{\prime}(ka)=0$ yields

[TABLE]

and the alternate expression for $\mathcal{E}_{\mathrm{surf}}$ follows. Next consider $(\gamma_{j})\in\ell^{2}(\mathbb{N})$ with $\gamma_{j}>r_{\rm hc}-a$ for all $j\in\mathbb{N}$. Under Assumption 1 the derivatives behave as $v^{\prime\prime}(r)=O(r^{-s-2})$ and $v^{\prime}(r)=O(r^{-s-1})$ as $r\to\infty$ with $s>2$. It follows that $\varepsilon_{j}:=\sum_{k=1}^{\infty}kv^{\prime}(ka)$ decays like $\int_{ja}^{\infty}r\times r^{-s-1}\mathrm{d}r=O(j^{-s+1})$ so that $\sum_{j=1}^{\infty}\varepsilon_{j}^{2}<\infty$. The Cauchy-Schwarz inequality then shows

[TABLE]

for some suitable $m$-independent constant $c$. In particular, when $m=\infty$ the sum $\sum_{j}\beta_{j}\gamma_{j}$ is absolutely convergent. In order to show that the double sum over $k$ and $j$ in Eq. (3.2) is absolutely convergent, we proceed with estimates analogous to Lemma 3.3. Assume first that all spacings $z_{j}=\gamma_{j}+a$ are larger than $z_{\rm min}$. Set $\sup_{r\geq z_{\min}}|v^{\prime\prime}(r)|=c_{1}$ and note that, by Assumption 1(iii) for all $k\geq 2$, $\sup_{r\geq kz_{\min}}|v^{\prime\prime}(r)|\leq|v^{\prime\prime}(kz_{\min})|$. Hence

[TABLE]

More generally, if $(\gamma_{j})\in\ell^{2}(\mathbb{N})\cap(r_{\rm hc}-a,\infty)^{\mathbb{N}}$, then $\gamma_{j}\to 0$ and because of $a>z_{\rm\min}$, there is an $i\in\mathbb{N}$ such that $z_{j}\geq z_{\min}$ for all $j\geq i$. Let $\varepsilon=\min\{|z_{j}|\mid j=1,\ldots,i\}$. Summands with $j\geq i$ can be estimated as before. For $j\leq i$ and $k\geq i+2$, we proceed as before as well, except that we replace $v^{\prime\prime}(kz_{\min})$ by $v^{\prime\prime}((k-i)z_{\min}+i\varepsilon)$. This leaves a finite sum over $j\leq i,k\leq i+2$ and overall, the sum is absolutely convergent. ∎

Lemma 3.6**.**

The map $\mathcal{D}\to\mathbb{R}$, (z_{j})\mapsto\mathcal{E}_{\mathrm{surf}}\bigl{(}(z_{j})_{j\in\mathbb{N}}\bigr{)} defined by (3.2) is continuous.

Proof.

Let $z,z^{\scriptscriptstyle{({1}})},z^{\scriptscriptstyle{({2}})},\ldots$ be sequences in $\mathcal{D}$ such that $z^{\scriptscriptstyle{({n}})}-z\to 0$ in $\ell^{2}(\mathbb{N})$. As $\lim_{i\to\infty}\sum_{j\geq i}(\gamma^{\scriptscriptstyle{({n}})}_{j})^{2}=0$ uniformly in $n$, the estimates above show that for every $\varepsilon>0$, we can find $i\in\mathbb{N}$ such that the sum over $\{(j,k)\mid j\geq i\text{ or }k\geq i\}$ contributes to $\mathcal{E}_{\mathrm{surf}}(\gamma^{\scriptscriptstyle{({n}})})$ and $\mathcal{E}_{\mathrm{surf}}(\gamma)$ an amount bounded by $\varepsilon$. In the remaining finite sum the continuity of $v(r)$ allows us to pass to the limit. The proof is easily concluded with an $\varepsilon/3$ argument. ∎

Lemma 3.7**.**

The restriction of $\mathcal{E}_{\mathrm{surf}}$ to $\mathcal{D}\cap[z_{\min},z_{\max}+\varepsilon]^{\mathbb{N}}$ is strictly convex and satisfies

[TABLE]

for suitable $m$-independent constants $\varepsilon,c_{1},c_{2}>0$.

Proof.

The proof of the convexity is similar to Lemma 3.3 and therefore omitted. For the coercivity, consider first $m=\infty$. Let $\gamma_{j}=z_{j}-a$, $\gamma_{j}^{\scriptscriptstyle{({n}})}=\gamma_{j}{\mathchoice{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.5mu\mathrm{l}}{1\mskip-5.0mu\mathrm{l}}}_{\{j\leq n\}}$ the truncated strain, and $z_{j}^{\scriptscriptstyle{({n}})}=a+\gamma_{j}^{\scriptscriptstyle{({n}})}$. Then

[TABLE]

thus

[TABLE]

where $C=-\sum_{k,\ell=1}^{\infty}v(\ell z_{\min}+ka)<\infty$. Next we cut and paste $(z_{1},\ldots,z_{n})$ into the middle of a large ground state chain: let $k_{1},k_{2}\in\mathbb{N}$ with $k_{2}\geq n+1$, $N=k_{2}+k_{1}+1$ and $(z_{-k_{1}+1}^{\scriptscriptstyle{({N}})},\ldots,z_{k_{2}}^{\scriptscriptstyle{({N}})})$ the spacings of the $N$-particle ground state. Let $z^{\prime}=(z_{-k_{1}+1}^{\scriptscriptstyle{({N}})},\ldots,z_{0}^{\scriptscriptstyle{({N}})},z_{1},\ldots,z_{n},z_{n+1}^{\scriptscriptstyle{({N}})},\ldots,z_{k_{2}}^{\scriptscriptstyle{({N}})})$. A Taylor expansion of $\mathcal{E}_{N}$ around the minmizer $z^{\scriptscriptstyle{({N}})}$ together with Lemma 3.3 and Theorem 2.1 yields

[TABLE]

On the other hand, let $C_{1}=\sum_{\ell=2}^{\infty}\ell|v(\ell z_{\min})|$ be a bound for interactions between blocks and remember $E_{k}\geq ke_{0}$ by Lemma 3.2 and $e_{0}\leq 0$. Then

[TABLE]

We combine with Eq. (3.3) and let first $k_{1},k_{2}\to\infty$, then $n\to\infty$, and conclude that $\frac{\eta}{2}\sum_{j=1}^{\infty}\gamma_{j}^{2}\leq\mathcal{E}_{\mathrm{surf}}(z)+C_{2}$ with the help of Lemma 3.6. This proves the coercivity in the case $m=\infty$. The proof for finite $m$ is similar. ∎

Lemma 3.8**.**

The surface energy $\mathcal{E}_{\mathrm{surf}}$ has a unique minimizer in $\mathcal{D}$. The minimizer is in $\mathcal{D}\cap[z_{\min},z_{\max}]^{\mathbb{N}}$.

Proof.

We proceed as in Section 3.1. Let $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}$. If one of the $z_{j}$’s is larger than $z_{\max}$, we can define a new configuration by shrinking this spacing to $z_{\max}$, leaving all other configurations unchanged. This decreases $\mathcal{E}_{\mathrm{surf}}$. If one of the $z_{j}$’s is smaller than $z_{\min}$, let $b$ be the smallest among them, and $j\in\mathbb{N}$ with $b=z_{j}$. Then we can define a new configuration by removing a participating particle and possibly shrinking a bond, i.e., $(z_{1},z_{2},\ldots)\mapsto(z_{1},z_{2},\ldots,z_{j-1},\min(z_{j}+z_{j+1},z_{\max}),z_{j+2},\ldots)$. Since $e_{0}\leq 0$, just as in Lemma 3.1, we see that this decreases the energy. Repeating these steps if necessary, the initial configuration is mapped to a new one that has strictly lower energy and all spacings in $[z_{\min},z_{\max}]$.

The existence of a minimizer now follows from the coercivity proven in Lemma 3.7, the compactness of $[z_{\min},z_{\max}]^{\mathbb{N}}\cap\mathcal{D}$ with respect to the weak $\ell^{2}$-convergence (shifted by $(a,a,\ldots)$) and the weak lower semicontinuity of $\mathcal{E}_{\mathrm{surf}}$ on that set due to Lemmas 3.6 and 3.7. The minimizer is unique because of the strict convexity from Lemma 3.7. ∎

Proof of Theorem 2.2..

Clear from Lemmas 3.6, 3.7, 3.8 and Proposition 3.4. ∎

Proof of Proposition 2.3..

In complete analogy to Lemma 3.5 we obtain

[TABLE]

for all $(z_{j})_{j\in\mathbb{Z}}\in\mathcal{D}^{+}_{0}$, and as in Lemma 3.6, we see that (3.4) defines a continuous map $\mathcal{D}^{+}\to\mathbb{R}$. The proof of strict convexity, even on $[z_{\min},z_{\max}+\varepsilon]^{\mathbb{Z}}\cap\mathcal{D}^{+}$ for some $\varepsilon>0$, is again similar to Lemma 3.3. As in Lemma 3.8 we have that $\mathcal{E}_{\mathrm{bulk}}$ has a unique minimizer in $\mathcal{D}$, which lies in $\mathcal{D}\cap[z_{\min},z_{\max}]^{\mathbb{N}}$. Since $a\in(z_{\min},z_{\max}]$ and $\partial_{i}\mathcal{E}_{\mathrm{bulk}}((z_{j})_{j\in\mathbb{Z}})=0$ for every $i\in\mathbb{Z}$ by (3.4), the minimizer of $\mathcal{E}_{\mathrm{bulk}}$ is $(\ldots,a,a,\ldots)$. Clearly, $\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)=0$. Finally, the formula connecting $\mathcal{E}_{\mathrm{bulk}}$ and $\mathcal{E}_{\mathrm{surf}}$ is clear on $\mathcal{D}^{+}_{0}$ and follows on $\mathcal{D}^{+}$ by approximation. ∎

3.3. A fixed point equation

In the following we assume that $v$ has a hard core:

Assumption 3**.**

$r_{\rm hc}>0$ and $v(r)\to\infty$ as $r\searrow r_{\rm hc}$.

We extend $h$, defined by (2.1) on $(r_{\rm hc},\infty)^{\mathbb{N}}$, to $\mathbb{R}_{+}^{\mathbb{N}}$ by setting

[TABLE]

Our main aim in this subsection is to obtain the following characterisation of $\overline{\mathcal{E}}_{\mathrm{surf}}$, cf. (2.4).

Proposition 3.9**.**

Let $I=\overline{\mathcal{E}}_{\mathrm{surf}}-\min\mathcal{E}_{\mathrm{surf}}$. Then $I$ is the unique lower semi-continuous solution (product topology) of the equation

[TABLE]

such that $\min I=0$ and $I=\infty$ if $z_{j}\leq r_{\mathrm{hc}}$ for one of the $z_{j}$’s.

Note that, by induction, (3.6) is equivalent to

[TABLE]

for all $k\in\mathbb{N}$ and $z=(z_{j})_{j\in\mathbb{N}}\in\mathbb{R}_{+}^{\mathbb{N}}$. (Observe that $h(z)>-\infty$ for all $z\in\mathbb{R}_{+}^{\mathbb{N}}$ by the decay assumption on $v$ and $r_{\rm hc}>0$.)

We begin with a technical auxiliary result.

Lemma 3.10**.**

If $z_{1},z_{2},\ldots>0$ and $\bar{c}<\infty$ are such that

[TABLE]

then $z=(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}$. Moreover, any $z\in\mathcal{D}$ satisfies

[TABLE]

Proof.

Let $\varepsilon_{0}<\min(a-z_{\min},z_{\max}-a)$. The partial sum $\sum_{j=1}^{k}h(z_{j},\ldots,z_{j+m-1})$ is equal to the energy $\mathcal{E}_{k+1}(z_{1},\ldots,z_{k})$ plus an interaction

[TABLE]

(the inner sum being [math] if $m+j-1<k+1$) which is bounded from below by

[TABLE]

By adding $n_{1}$ and $n_{2}$ spacings $a$ to the left and right respectively, we may view $z$ as a block of spacings in an $N$-particle configuration where $N=n_{1}+n_{2}+k+1$. Let $\hat{z}=(a,\ldots,a,z_{1},\ldots,z_{k},a,\ldots,a)$. The new configuration satisfies

[TABLE]

for some suitable constant $C$ that depends on $r_{\mathrm{hc}}$, $\bar{c}$ and $v$ only. Let $z^{\scriptscriptstyle{({N}})}$ be the $N$-particle ground state with spacings labelled by $j=-n_{1}+1,\ldots,k+n_{2}$ rather than $1,\ldots,N-1$. Since $\mathcal{E}_{N}(z^{\scriptscriptstyle{({N}})})=E_{N}\geq Ne_{0}$ by Lemma 3.2 and $e_{0}\leq 0$, we get

[TABLE]

Suppose that all spacings $z_{j}$ are in $[z_{\min},z_{\max}]$. We use a Taylor approximation around the minimizer $z^{\scriptscriptstyle{({N}})}$, apply Lemma 3.3 and Theorem 2.1, and obtain.

[TABLE]

Letting $k\to\infty$ we obtain an upper bound for the $\ell^{2}$-norm of $(z_{j}-a)_{j\in\mathbb{N}}$. If there are $z_{j}$ with $z_{j}<z_{\min}$ or $z_{j}>z_{\max}$, we modify the configuration $z_{1},\ldots,z_{k}$ without increasing its energy as in the proof of Lemma 3.1 to obtain $z^{\prime}_{1},\ldots,z^{\prime}_{k}$. When we shrink bonds $z_{j}>z_{\max}$ to $z^{\prime}_{j}=z_{\max}$, leaving all other spacings unchanged, both $z^{\prime}_{j}$ and $z_{j}$ are strictly larger than $\varepsilon_{0}$ so the truncated $\ell^{2}$-norm \sum_{j=1}^{k}\min\bigl{(}(z_{j}-a)^{2},\varepsilon_{0}^{2}) is unaffected.

On the other hand suppose $z_{i}=\min(z_{j})<z_{\min}$. Then we remove the particle $x_{i}$, reattach it a distance $z_{\max}$ to the left of the $k$-particle block. This effects the change

[TABLE]

on the $\ell^{2}$-norm. Both $|z_{i}-a|$ and $|z_{\max}-a|$ are larger than $\varepsilon_{0}$, moreover

[TABLE]

So the truncated $\ell^{2}$-norm increases by at most $\varepsilon_{0}^{2}$. Let $n$ be the number of times this step has to be performed. Iterating we arrive at a configuration $z^{\prime\prime}_{1},\ldots,z^{\prime\prime}_{k}\in[z_{\min},z_{\max}]$ with

[TABLE]

and $\mathcal{E}_{k+1}(z^{\prime\prime})\leq\mathcal{E}_{k+1}(z)-n\delta$ for some $\delta>0$, cf. (3.1). Making $\varepsilon_{0}$ smaller if necessary we may assume $\varepsilon_{0}^{2}<\delta$. We combine with Eq. (3.8) for $\hat{z}^{\prime\prime}$ and $C^{\prime\prime}=C-n\delta$ and obtain

[TABLE]

We let $k\to\infty$ and find that the truncated $\ell^{2}$-norm of $(z_{j})_{j\in\mathbb{N}}$ is finite. It follows in particular that there are only finitely many spacings $|z_{j}-a|\geq\varepsilon_{0}$, and $(z_{j}-a)_{j\in\mathbb{N}}$ is square summable. This establishes the first assertion.

In order to show the convergence of the partial sums to $\mathcal{E}_{\mathrm{surf}}$, first observe that $\mathcal{E}_{\mathrm{surf}}$ satisfies (3.7) for $I=\mathcal{E}_{\mathrm{surf}}$. This is clear for $z\in\mathcal{D}_{0}$ and follows for general $z\in\mathcal{D}$ by continuity. If $z\in\mathcal{D}$, the sequence of shifts $((z_{j})_{j\geq k})_{k\in\mathbb{N}}$ converges to $(\ldots,a,a,\ldots)$ strongly and thus

[TABLE]

as $k\to\infty$. ∎

We have actually proven the following: for sufficiently small $\varepsilon_{0}>0$, suitable $c_{1},c_{2}>0$, and all $(z_{j})_{j\in\mathbb{N}}\in\mathbb{R}_{+}^{\mathbb{N}}$,

[TABLE]

Proof of Proposition 3.9..

Let $I=\overline{\mathcal{E}}_{\mathrm{surf}}-\min\mathcal{E}_{\mathrm{surf}}$. Observe that $I$ satisfies (3.6). This is clear for $z\in\mathcal{D}_{0}$ and for $z\notin\mathcal{D}$. For the remaining $z$ it follows from Lemma 3.6. We now show that $I$ is lower semi-continuous with respect to pointwise convergence. Without loss we suppose that $z^{\scriptscriptstyle{({n}})}\in\mathcal{D}$ converges to $z\in[r_{\rm hc},\infty)^{\mathbb{N}}$ pointwise with $I(z^{\scriptscriptstyle{({n}})})\leq\bar{c}<\infty$ for some constant $\bar{c}>0$. Passing to a subsequence (not relabelled) we may furthermore assume that $\liminf_{n\to\infty}I(z^{\scriptscriptstyle{({n}})})=\lim_{n\to\infty}I(z^{\scriptscriptstyle{({n}})})$. Fix an $\varepsilon>0$ such that the estimate in Lemma 3.7 is satisfied. By (3.9)

[TABLE]

for some uniform constant $C>0$ since $z_{\min}<a\leq z_{\max}$. For given $N\in\mathbb{N}$ we denote by $j_{n}$ the first index $j\geq N$, if existent, with $z^{\scriptscriptstyle{({n}})}_{j}\notin[z_{\min},z_{\max}+\varepsilon]$. Passing to a further subsequence (not relabelled) and choosing $N$ sufficiently large we may achieve that either such indices do not exist or that $j_{n}\to\infty$ as $n\to\infty$. In both cases we get that $z_{j}\in[z_{\min},z_{\max}+\varepsilon]$ for $j\geq N$. In particular, $z_{j}>r_{\rm hc}$ for $j\geq N$.

In the second case we define new configurations $\tilde{z}^{\scriptscriptstyle{({n}})}$ by applying the procedure detailed in the proof of Lemma 3.8 to the tails $(z^{\scriptscriptstyle{({n}})}_{j})_{j\geq N}$ shrinking the bonds $z^{\scriptscriptstyle{({n}})}_{j}>z_{\max}+\varepsilon$, $j\geq N$, and deleting particles $x^{\scriptscriptstyle{({n}})}_{j+1}$ if $z^{\scriptscriptstyle{({n}})}_{j}<z_{\min}$, $j\geq N$, so that

[TABLE]

In the first case we simply set $\tilde{z}^{\scriptscriptstyle{({n}})}=z^{\scriptscriptstyle{({n}})}$. Since $j_{n}\to\infty$ in the second case, we still have $\tilde{z}^{\scriptscriptstyle{({n}})}\to z$ pointwise.

By (3.7) with $k=N-1$ we have

[TABLE]

From the decay properties of $v$ and $z^{\scriptscriptstyle{({n}})}_{j}\geq r_{\rm hc}>0$ it is easy to see that, for any $j\in\mathbb{N}$, $h(z^{\scriptscriptstyle{({n}})}_{j},z^{\scriptscriptstyle{({n}})}_{j+1},\ldots)$ converges to $h(z_{j},z_{j+1},\ldots)$. Since $I(z^{\scriptscriptstyle{({n}})})\leq\bar{c}$ and $I\geq 0$, from Assumption 3 we also get $z_{j}>r_{\rm hc}$ for $j=1,\ldots,N-1$. So

[TABLE]

In particular, $I((\tilde{z}^{\scriptscriptstyle{({n}})}_{j})_{j\geq N})\leq C$ and so Lemma 3.7 implies that $z\in\mathcal{D}$ and $\tilde{z}^{\scriptscriptstyle{({n}})}-z\rightharpoonup 0$ in $\ell^{2}$ by coercivity and hence that

[TABLE]

by convexity. Summarizing we obtain

[TABLE]

Suppose, conversely, that a lower semi-continuous $I:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}\cup\{+\infty\}$ satisfies (3.6) with $\min I=0$ and $I(z)=\infty$ if $z_{j}\leq r_{\mathrm{hc}}$ for some $j$. We first note that, since $I\geq 0$, for any $z$ with $I(z)<\infty$ one has

[TABLE]

by (3.7) and so $z\in\mathcal{D}$ by Lemma 3.10. It thus suffices to show that

[TABLE]

for all $z\in\mathcal{D}$.

If $z\in\mathcal{D}$, then $\mathcal{E}_{\mathrm{surf}}(z)$ is indeed finite by Lemma 3.5. We have \lim_{k\to\infty}\sum_{j=1}^{k}\big{(}h(z_{j},\ldots,z_{j+m-1})-e_{0}\big{)}=\mathcal{E}_{\mathrm{surf}}(z) by Lemma 3.10. Since the sequence of shifts $((z_{j})_{j\geq k})_{k\in\mathbb{N}}$ converges to $(a,a,\ldots)$ pointwise as $k\to\infty$, taking the $\liminf$ in (3.7) yields

[TABLE]

Note that, as $I\not\equiv\infty$, this inequality also shows that $I(a,a,\ldots)<\infty$.

For the reverse inequality, by choosing $k$ large enough in (3.7) we first see that (3.10) holds true for all $z\in\mathcal{D}_{0}$. We denote by $z^{(N)}$ the truncation with $z^{(N)}_{j}=z_{j}$ for $j\leq N$ and $z^{(N)}_{j}=a$ for $j\geq N+1$. Since $z^{(N)}\to z$ pointwise and $z^{(N)}-z\to 0$ in $\ell^{2}$ as $N\to\infty$, lower semi-continuity of $I$ and strong continuity of $\mathcal{E}_{\mathrm{surf}}$ (see Lemma 3.6) give

[TABLE]

where we have used that $z^{(N)}\in\mathcal{D}_{0}$ for all $N$. ∎

We now restrict to the case $m<\infty$. Let $d=m-1$. By (3.7) with $k=d$ we have

[TABLE]

for any $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}$, where

[TABLE]

Taking the infimum over $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$, with fixed $z_{1},\ldots,z_{d}$ and setting

[TABLE]

(recall Lemma 3.6) leads to

[TABLE]

In Chapter 6 we will need the following estimate.

Lemma 3.11**.**

Set $A_{\varepsilon}=[z_{\min},z_{\max}+\varepsilon]^{d}$ and $B_{\varepsilon}=\mathbb{R}_{+}^{d}\setminus A_{\varepsilon}$. Then, for any $\varepsilon>0$ there exists a $\delta>0$ such that

[TABLE]

for all $x\in A_{\varepsilon}$.

Proof.

Suppose $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$ is such that $(z_{1},\ldots,z_{d})\in A_{\varepsilon}$, in particular, $z_{j}\geq z_{\min}$ for $j=1,\ldots,d$. If $(z_{d+1},z_{d+2},\ldots)\notin[z_{\min},z_{\max}+\varepsilon]^{\mathbb{N}}$ we construct a new configuration $(z^{\prime}_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$ without changing the first $d$ spacings similarly as in the proofs of Lemma 3.1 and 3.8.

If $z_{i}>z_{\max}+\varepsilon$, we define $(z^{\prime}_{j})_{j\in\mathbb{N}}$ by setting $z^{\prime}_{j}=z_{j}$ for $j\neq i$ and $z^{\prime}_{i}=z_{\max}$. Then

[TABLE]

Now assume $b=\min\{z_{d+1},z_{d+2},\ldots\}<z_{\min}$. We choose an $i\geq d+1$ with $z_{i}=b$ and define $(z^{\prime}_{j})_{j\in\mathbb{N}}$ by setting $z^{\prime}_{j}=z_{j}$ for $j<i$, $z^{\prime}_{i}=\min\{z_{i}+z_{i+1},z_{\max}\}$ and $z^{\prime}_{j}=z_{j+1}$ for $j>i$. As in Lemmas 3.1 and 3.8 (in particular using that $e_{0}\leq 0$), we see that

[TABLE]

The estimates (3.12) and (3.13) show that, for any $(z_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$ with $(z_{1},\ldots,z_{d})\in A_{\varepsilon}$ and $(z_{d+1},\ldots,z_{2d})\in B_{\varepsilon}$ there is a $(z^{\prime}_{j})_{j\in\mathbb{N}}\in\mathcal{D}_{0}$ with $(z^{\prime}_{1},\ldots,z^{\prime}_{d})=(z_{1},\ldots,z_{d})$ such that

[TABLE]

where \delta=\min\big{\{}v(z_{\max}+\varepsilon)-v(z_{\max}),\ 2\alpha_{1}\sum_{n=m+1}^{\infty}(nz_{\min})^{-s}\big{\}}>0. Using (3.11) we arrive at

[TABLE]

The claim now follows by taking the infimum over $(z_{j})_{j\in\mathbb{N}}$ with fixed $(z_{1},\ldots,z_{d})$ conditioned on $(z_{d+1},\ldots,z_{2d})\in B_{\varepsilon}$. ∎

A simpler proof gives the following estimate that will also be needed in Chapter 6.

Lemma 3.12**.**

For any $\varepsilon>0$ there exists a $\delta>0$ such that $\mathcal{E}_{\mathrm{bulk}}(z)\geq\delta$ for all $z\in\mathcal{D}^{+}\setminus[z_{\min},z_{\max}+\varepsilon]^{\mathbb{Z}}$.

Proof.

By continuity we may assume that $z=(z_{j})_{j\in\mathbb{Z}}\in\mathcal{D}^{+}_{0}\setminus[z_{\min},z_{\max}+\varepsilon]^{\mathbb{Z}}$. If $z_{i}>z_{\max}+\varepsilon$, we define $z^{\prime}=(z^{\prime}_{j})_{j\in\mathbb{Z}}$ by setting $z^{\prime}_{j}=z_{j}$ for $j\neq i$ and $z^{\prime}_{i}=z_{\max}$. Then

[TABLE]

If $b=\min\{z_{j}:j\in\mathbb{Z}\}<z_{\min}$. We choose the smallest $i$ with $z_{i}=b$ and define $z=(z^{\prime}_{j})_{j\in\mathbb{N}}$ by setting $z^{\prime}_{j}=z_{j}$ for $j<i$, $z^{\prime}_{i}=\min\{z_{i}+z_{i+1},z_{\max}\}$ and $z^{\prime}_{j}=z_{j+1}$ for $j>i$. As in (3.13) we get

[TABLE]

This concludes the proof. ∎

4. Gibbs measures for the infinite and semi-infinite chains

Here we prove the existence of $\nu_{\beta}$, $\mu_{\beta}$, $g(\beta)$, $g_{\mathrm{surf}}(\beta)$ and check that $\mu_{\beta}$ is shift-invariant and mixing, hence ergodic; the results and methods are fairly standard. In addition, we provide an a priori estimate on the decay of correlations with explicit analysis of the $\beta$-dependence (Theorem 4.4) which to the best of our knowledge is new. The results from this section need only very little on the pair potential: we only use that $v$ has a hard core and that $v(r)=O(1/r^{s})$, for large $r$, with $s>2$. The technical assumption of a hard core frees us from superstability estimates [LP76, Rue76]. The decay of the potential ensures that the infinite volume Gibbs measure is unique, see e.g. [Geo11, Chapter 8.3] and [Pap84a, Pap84b, Kle85].

We follow the classical treatment of one-dimensional systems with transfer operators. For compactly supported pair potentials with a hard core (or, in our case, when $m$ is chosen finite), the transfer operators are integral operators in $L^{2}(\mathbb{R}_{+}^{m-1},\mathrm{d}x)$ [Rue69, Chapter 5.6], see Section 6. For long-range interactions, the transfer operator (also known as Ruelle operator or Ruelle-Perron-Frobenius operator) acts instead from the left on functions of infinitely many variables, and from the right on measures [Rue68, GMS70, Rue78]. The formalism of transfer operators keeps being developed in the context of dynamical systems and ergodic theory [Bal00b, Bal00a].

For the decay of correlations, we adapt [Pol00] to the present context of continuous unbounded spins and carefully track the $\beta$-dependence in the bounds. In Section 5.3, transfer operators will also help us investigate the large deviations behavior of the Gibbs measures; notably the eigenvalue equation from Lemma 4.1 translates into a fixed point equation for the rate function (see Lemma 5.4).

The results of this section hold for all $m\in\mathbb{N}\cup\{\infty\}$ and $\beta,p>0$; the additional condition $p<p^{*}$ is not needed.

4.1. Transfer operator

For $j\in\mathbb{Z}$ and $z_{j},z_{j+1},\ldots>0$ we abbreviate $h_{j}=h(z_{j},z_{j+1},\ldots)$, cf. (2.1) and (3.5). The transfer operator acts on functions as

[TABLE]

The dual action on measures is defined by $(\mathcal{L}_{\beta}^{*}\nu)(f)=\nu(\mathcal{L}_{\beta}f)$ and is given by

[TABLE]

Lemma 4.1**.**

There exist $\lambda_{0}(\beta)>0$ and a probability measure $\nu_{\beta}$ on $\mathbb{R}_{+}^{\mathbb{N}}$ such that

[TABLE]

Moreover $\nu_{\beta}((r_{\mathrm{hc}},\infty)^{\mathbb{N}})=1$ and the pair $(\nu_{\beta},\lambda_{0}(\beta))$ is unique.

We will show in Proposition 4.9 that $\nu_{\beta}$ is the measure satisfying (2.2). The non-compactness of $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ forms an obstacle to the application of a Schauder-Tychonoff fixed point theorem for the map $\nu\mapsto\mathcal{L}_{\beta}^{*}\nu/\nu(\mathcal{L}_{\beta}\mathbf{1})$, see e.g. [Rue68, Proposition 2]. It might be possible to remove the obstacle using tightness estimates, but we prefer to follow a different route and exploit the known uniqueness of infinite volume Gibbs measures [Geo11, Chapter 8.3] instead.

Proof.

Let $\nu$ be a probability measure on $\mathbb{R}_{+}^{\mathbb{N}}$, $\lambda:=\nu(\mathcal{L}_{\beta}\mathbf{1})$, and $\tilde{\nu}:=\frac{1}{\lambda}\mathcal{L}_{\beta}^{*}\nu$. We show that if $\nu$ is a Gibbs measure, then $\tilde{\nu}$ is a Gibbs measure as well. Let us first introduce the kernels needed to formulate that $\nu$ is a Gibbs measure. By [Geo11, Theorem 1.33] it is enough to look at one-point kernels. Pick $k\in\mathbb{N}$. For $z^{\prime}_{k}>0$ and $z=(z_{j})_{j\in\mathbb{N}}\in\mathbb{R}_{+}^{\mathbb{N}}$, let

[TABLE]

where sum runs over discrete intervals $J=\{i,\ldots,\ell-1\}\subset\mathbb{N}$. Further define the kernel

[TABLE]

where $A\subset\mathbb{R}_{+}^{\mathbb{N}}$ and $N_{k}(z)=\int_{0}^{\infty}{\operatorname{e}}^{-\beta H_{k}(z^{\prime}_{k}\mid z)}\mathrm{d}z^{\prime}_{k}$. The kernel acts on functions and measures in the usual way, in particular $(\gamma_{k}{\mathchoice{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.5mu\mathrm{l}}{1\mskip-5.0mu\mathrm{l}}}_{A})(z)=\gamma_{k}(z,A)$. Notice that $\gamma_{k}^{2}f=\gamma_{k}f$ for all $f$. Indeed $\gamma_{k}f$ yields a function where $z_{k}$-dependence has been integrated out, and integrating it against the probability measure $\gamma_{k}(z,\cdot)$ does not change its value. Replacing $\mathbb{N}$ with $\mathbb{N}_{0}$, we define in a completely analogous fashion conditional energies $H_{k}^{0}$ and kernels \gamma_{k}^{0}\bigl{(}(z_{j})_{j\in\mathbb{N}_{0}},B\bigr{)}.

Suppose that $\nu$ is a Gibbs measure, i.e., $\nu\gamma_{k}=\nu$ for all $k\in\mathbb{N}$. Let $f:\mathbb{R}_{+}^{\mathbb{N}_{0}}\to\mathbb{R}_{+}$ be a measurable test function. Treat $\tilde{\nu}=\lambda^{-1}\mathcal{L}_{\beta}^{*}\nu$ as a measure on $\mathbb{R}_{+}^{\mathbb{N}_{0}}$. We check that $\tilde{\nu}(\gamma_{k}^{0}f)=\tilde{\nu}(f)$ for all $k\in\mathbb{N}_{0}$. For $k\in\mathbb{N}$, this property is inherited from the Gibbsianness of $\nu$: we have

[TABLE]

Set $\tilde{f}:=\gamma_{k}^{0}f$. Note $\tilde{f}=(\gamma_{k}^{0})\tilde{f}$. Therefore

[TABLE]

hence $\tilde{\nu}(f)=\tilde{\nu}(\tilde{f})=\tilde{\nu}(\gamma_{k}^{0}f)$. For $k=0$, the required property follows from the definition of $\tilde{\nu}$. Notice $H_{0}^{0}=h_{0}$ and

[TABLE]

Let $\tilde{f}=\gamma_{0}^{0}f$. Then

[TABLE]

The previous identities hold for all non-negative test functions $f$, consequently $\tilde{\nu}\gamma_{k}^{0}=\tilde{\nu}$ for all $k\in\mathbb{N}_{0}$ and $\tilde{\nu}$ is a Gibbs measure as well.

By [Geo11, Theorem 8.39], the Gibbs measure $\nu$ exists and is unique. Treating $\nu$ and $\tilde{\nu}$ both as measures on $\mathbb{R}_{+}^{\mathbb{N}}$, we must therefore have $\nu=\tilde{\nu}$, i.e., the unique Gibbs measure is an eigenmeasure of $\mathcal{L}_{\beta}^{*}$ and in particular, there exists an eigenmeasure. Conversely, let $\nu=\frac{1}{\lambda}\mathcal{L}_{\beta}^{*}\nu$ be an eigenmeasure. Arguments similar to the investigation of $\tilde{\nu}$ given above, based on the iterated fixed point equation $\nu=\frac{1}{\lambda^{k}}{\mathcal{L}_{\beta}^{*}}^{k}\nu$, show that $\nu\gamma_{j}=\nu$ for all $j=1,\ldots,k$ and all $k$, hence for all $j$. Every eigenmeasure is a Gibbs measure. Since the latter is unique, the eigenmeasure is unique as well. Finally, since $v(z_{j})=\infty$ for $z_{j}\leq r_{\mathrm{hc}}$, the eigenmeasure $\nu=\frac{1}{\lambda^{k}}{\mathcal{L}_{\beta}^{*}}^{k}\nu$ must satisfy $\nu(\exists j\in\{1,\ldots,k\}:\,z_{j}\leq r_{\mathrm{hc}})=0$. This holds for all $k\in\mathbb{N}$, hence $\nu((r_{\mathrm{hc}},\infty)^{\mathbb{N}})=1$. ∎

Let $\nu_{\beta}^{-}$ be the probability measure on $\mathbb{R}_{+}^{\{\ldots,-1,0\}}$ obtained by flipping $\nu_{\beta}^{+}=\nu_{\beta}$, i.e., $\nu_{\beta}^{-}$ is the image of $\nu_{\beta}^{+}=\nu_{\beta}$ under the map $(z_{k})_{k\in\mathbb{N}}\mapsto(z_{1-\ell})_{\ell\leq 0}$. The measures $\nu_{\beta}^{\pm}$ represent equilibrium measures for the left and right half-infinite chains. Let

[TABLE]

be the total interaction between left and right half-infinite chains, cf. Proposition 2.3(d). We abbreviate the shifted versions as $\mathcal{W}_{\ell}=\mathcal{W}(\cdots z_{\ell}\mid z_{\ell+1}\cdots)$. Define $\varphi_{\beta}(z_{1},z_{2},\ldots)$ by

[TABLE]

Thus $\varphi_{\beta}(z_{1},z_{2},\ldots)$ represents an averaged contribution to the Boltzmann weight from the left half-infinite chain.

Lemma 4.2**.**

We have $\mathcal{L}_{\beta}\varphi_{\beta}=\lambda_{0}(\beta)\varphi_{\beta}$ and $\nu_{\beta}(\varphi_{\beta})=1$.

Proof.

The normalization is obvious, for the eigenvalue equation let $c_{\beta}=\nu_{\beta}^{-}\otimes\nu_{\beta}^{+}(\exp(-\beta\mathcal{W}_{0}))$ and use the eigenvalue equation for $\nu_{\beta}^{\pm}$

[TABLE]

See also [Rue78, Section 5.12]. ∎

Define the operator

[TABLE]

so that $\mathcal{S}_{\beta}\mathbf{1}=\mathbf{1}$ and $\mathcal{S}_{\beta}^{*}(\varphi_{\beta}\nu_{\beta}^{+})=\varphi_{\beta}\nu_{\beta}^{+}$. Let $\mu_{\beta}$ be the probability measure on $\mathbb{R}_{+}^{\mathbb{Z}}$ given by

[TABLE]

We will show in Proposition 4.9 that $\mu_{\beta}$ is the measure satisfying (2.3). Notice that for every bounded measurable function $f$ that depends on right-chain variables $z_{1},z_{2},\ldots$ only,

[TABLE]

Let $\tau:\mathbb{R}_{+}^{\mathbb{Z}}\to\mathbb{R}_{+}^{\mathbb{Z}}$ be the shift $(\tau z)_{j}=z_{j+1}$.

Lemma 4.3**.**

• (a)

$\mu_{\beta}$* is shift-invariant.*

• (b)

For all $f,g:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}_{+}$ and all $n\in\mathbb{N}$, we have $\mu_{\beta}(f(g\circ\tau^{n}))=\mu_{\beta}((\mathcal{S}_{\beta}^{n}f)g)$.

The proof is standard [Rue78] and therefore omitted. The lemma can be rephrased as follows: let $(Z_{n})_{n\in\mathbb{Z}}$ be a stochastic process with law $\mu_{\beta}$, defined on some probability space $(\Omega,\mathcal{F},\mathbb{P})$. Then $(Z_{n})_{n\in\mathbb{Z}}$ is stationary, and

[TABLE]

Our next task is to show that the process is not only stationary but in fact ergodic and to estimate the decay of correlations.

4.2. Ergodicity

Bounds on correlations are most conveniently expressed with the help of variations, semi-norms that quantify how much a function depends on faraway variables. Notice that $\nu_{\beta}((r_{\rm hc},\infty)^{\mathbb{N}})=\mu_{\beta}((r_{\rm hc},\infty)^{\mathbb{Z}})=1$. Let $f:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}$ be a function and $n\in\mathbb{N}$. The $n$th variation of $f$ on $(r_{\rm hc},\infty)^{\mathbb{N}}$ is

[TABLE]

When $n=0$ the constraint on initial values is empty, $\operatorname{var}_{0}(f)$ is sometimes called the oscillation of $f$ [Geo11, Eq. (8.2)]. The oscillation vanishes if and only $f$ is constant. Notice that $\operatorname{var}_{k}(h)$ decays algebraically: for $k\in\mathbb{N}$, as $v(r)=O(r^{-s})$,

[TABLE]

It follows that the variation is summable, $\sum_{k=1}^{\infty}\operatorname{var}_{k}(h)<\infty$. Set

[TABLE]

Notice that for all $q\in\mathbb{N}_{0}$, $C_{q}$ is independent of $\beta$ and $p$. In fact the pressure only enters the oscillation $\operatorname{var}_{0}(h)$. By a slight abuse of notation we identify a function $f:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}$ with the function $f_{1}:\mathbb{R}_{+}^{\mathbb{Z}}\to\mathbb{R}_{+}$, $(z_{j})_{j\in\mathbb{Z}}\mapsto f((z_{j})_{j\in\mathbb{N}})$ and write $\mu_{\beta}(f)$ instead of $\mu_{\beta}(f_{1})$. The results of this subsection hold for all $p>0$.

Theorem 4.4**.**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p>0$. The measure $\mu_{\beta}$ is mixing with respect to shifts, i.e., $\mu_{\beta}(f(g\circ\tau^{n}))\to\mu_{\beta}(f)\mu_{\beta}(g)$ as $n\to\infty$, for all $f,g\in L^{1}(\mathbb{R}_{+}^{\mathbb{Z}},\mu_{\beta})$. Moreover for $\gamma(\beta)=\exp(-3\beta C_{0})$ and all bounded $f,g:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}$, $q,n\in\mathbb{N}$, $N\geq qn$,

[TABLE]

We prove Theorem 4.4 with Pollicott’s method of conditional expectations [Pol00]. For alternative approaches, see [Sar02] and the references therein. The principal idea is the following: for $n\in\mathbb{N}$, $f\in L^{1}(\mathbb{R}_{+}^{\mathbb{N}},\varphi_{\beta}\nu_{\beta})$ let $\Pi_{n}f$ be the projection

[TABLE]

onto the subspace of functions that depend on the first $n$ coordinates only, i.e., $\operatorname{var}_{n}(f)=0$. In terms of the stationary process $(Z_{n})_{n\in\mathbb{Z}}$ with law $\mu_{\beta}$,

[TABLE]

Notice that

[TABLE]

where $||\cdot||_{1}$ is the $L^{1}(\mathbb{R}_{+}^{\mathbb{N}},\varphi_{\beta}\nu_{\beta})$ norm. Let $q,n\in\mathbb{N}$. Then

[TABLE]

The difference enclosed in parentheses represents a truncation error; it is made small by choosing $n$ large. On the subspace of mean-zero functions, the truncated operator $\mathcal{S}^{n}\Pi_{n}$ satisfies a contraction property uniformly in $n$ (Lemma 4.7), and $(\mathcal{S}_{\beta}^{n}\Pi_{n})^{q}$ goes to zero exponentially fast as $q\to\infty$.

Lemma 4.5**.**

We have $\operatorname{var}_{q}(\log\varphi_{\beta})\leq\beta C_{q}$ for all $q\in\mathbb{N}_{0}$ and $\beta,p>0$.

Proof.

Let $q\in\mathbb{N}_{0}$, $(z_{j})_{j\in\mathbb{Z}},(z^{\prime}_{j})_{j\in\mathbb{Z}}\in(r_{\rm hc},\infty)^{\mathbb{Z}}$ such that $z_{j}=z^{\prime}_{j}$ for all $j\leq q$. Then

[TABLE]

and $\nu_{\beta}^{-}(\exp(-\beta\mathcal{W}_{0}))\leq\exp(\beta C_{q})\nu_{\beta}^{-}(\exp(-\beta\mathcal{W}_{0}^{\prime}))$. The claim then follows from the definition (4.1) of the invariant function. ∎

Lemma 4.6**.**

Let $f:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}$ be a bounded function. Then $n,k\in\mathbb{N}_{0}$,

[TABLE]

Proof.

Let $g=\sum_{j=1}^{n}h_{j}-\beta^{-1}\log[\lambda_{0}^{n}(\beta)\varphi_{\beta}]+\beta^{-1}\log\varphi_{\beta}\circ\tau^{n}$ on $(r_{\rm hc},\infty)^{\mathbb{N}}$ and $g\equiv\infty$ on $\mathbb{R}_{+}^{\mathbb{N}}\setminus(r_{\rm hc},\infty)^{\mathbb{N}}$ so that

[TABLE]

Pick $z,z^{\prime}\in(r_{\rm hc},\infty)^{\mathbb{N}}$ so that $z_{j}=z^{\prime}_{j}$ for $j=1,\ldots,n+k$. Then

[TABLE]

We integrate out $z_{1},\ldots,z_{n}$, observe $\int\exp(-\beta g)\mathrm{d}z_{1}\cdots\mathrm{d}z_{n}=\mathcal{S}_{\beta}^{n}\mathbf{1}=\mathbf{1}$, and deduce

[TABLE]

To conclude, we note

[TABLE]

∎

Lemma 4.7**.**

Let $f\in L^{1}(\mathbb{R}_{+}^{\mathbb{N}},\varphi_{\beta}\nu_{\beta})$ such that $\nu_{\beta}(f\varphi_{\beta})=0$. Then for all $n\geq 1$ and $\gamma(\beta)=\exp(-3\beta C_{0})$

[TABLE]

Proof.

We adapt [Rue68, Proposition 3]. Consider first a non-negative function $f$ that depends on $z_{1},\ldots,z_{n}$ only, i.e., $\operatorname{var}_{n}(f)=0$. Let $k\geq 0$ $z,z^{\prime}$ such that $z_{j}=z^{\prime}_{j}$ for $j=1,\ldots,n$ and $g(z_{1},z_{2},\ldots)$ as in the proof of Lemma 4.6. Then

[TABLE]

By Inequality (4.5) with $k=0$ we have $\operatorname{var}_{n}(g)\leq 3C_{0}$, uniformly in $n$. Thus $\mathcal{S}_{\beta}^{n}f(z)\leq\exp(-3\beta C_{0})(\mathcal{S}_{\beta}^{n}f)(z^{\prime})$ for all $z,z^{\prime}\in(r_{\rm hc},\infty)^{\mathbb{N}}$. For non-negative $f$ with $f=\Pi_{n}f$ we have by Lemma 4.3

[TABLE]

Next let $f$ with $\operatorname{var}_{n}(f)=0$ and $\mu_{\beta}(f)=0$. Then $\mu_{\beta}(f_{+})=\mu_{\beta}(f_{-})$ and

[TABLE]

We integrate against $\mu_{\beta}$, use $\mu_{\beta}(\mathcal{S}_{\beta}^{n}|f|)=\mu_{\beta}(|f|)=||f||_{1}$, and find $||\mathcal{S}_{\beta}^{n}f||_{1}\leq(1-\gamma(\beta))||f||_{1}$. This holds for every local function $\operatorname{var}_{n}(f)=0$ with $\mu_{\beta}(f)=0$. For general $f$, we may apply the bound to $\Pi_{n}f$ and use $\mu_{\beta}(\Pi_{n}f)=\mu_{\beta}(f)=0$ and $\mu_{\beta}(|\Pi_{n}f|)\leq\mu_{\beta}(|f|)$, and we are done. ∎

Lemma 4.8**.**

Let $f\in L^{1}(\mathbb{R}_{+}^{\mathbb{N}},\varphi_{\beta}\nu_{\beta})$ be a bounded map with $\nu_{\beta}(f\varphi_{\beta})=0$. Then for all $q,n\in\mathbb{N}$,

[TABLE]

Proof.

A telescope summation, the triangle inequality, and Lemma 4.7 yield

[TABLE]

where in the second step we use that \nu_{\beta}((\mathcal{S}_{\beta}^{n}\Pi_{n})^{i}\bigl{(}\mathcal{S}_{\beta}^{n}\Pi_{n}-\mathcal{S}_{\beta}^{n}\bigr{)}(\mathcal{S}_{\beta}^{n})^{q-k-1}f\varphi_{\beta})=\nu_{\beta}(f\varphi_{\beta})=0 for $i=1,\ldots,k$ by Lemma 4.3 and the third step follows from |\mathcal{S}_{\beta}^{n}\bigl{(}\Pi_{n}-\mathrm{id}\bigr{)}(\mathcal{S}_{\beta}^{n})^{q-k-1}f|\leq\mathcal{S}_{\beta}^{n}|\bigl{(}\Pi_{n}-\mathrm{id}\bigr{)}(\mathcal{S}_{\beta}^{n})^{q-k-1}f| and Lemma 4.3. By Eq. (4.4) and Lemma 4.6, this can be further estimated as

[TABLE]

Proof of Theorem 4.4.

Let $f,g:\mathbb{R}_{+}^{\mathbb{N}}\to\mathbb{R}$ be bounded functions and $q,n\in\mathbb{N}$, $N\geq qn$. Using Eq. (4.2) and Lemmas 4.7 and 4.8, we get

[TABLE]

since $||\mathcal{S}_{\beta}||_{1}\leq 1$. The explicit estimate on the decay of correlations follows. That $\mu_{\beta}$ is mixing then follows from standard approximation arguments. ∎

Proof of Theorem 2.9.

The estimate for infinite $m$ is an immediate consequence of Theorem 2.9. For finite $m$ and $n=m-1$, the truncation error in Lemma 4.8 for a function $f:\mathbb{R}_{+}^{n}\to\mathbb{R}$ actually vanishes since $\operatorname{var}_{n}(f)=0$ and $C_{n}=0$. The bound simplifies accordingly. ∎

4.3. Thermodynamic limit

Proposition 4.9**.**

Let $m\in\mathbb{N}\cup\{\infty\}$ and $p>0$.

1. (a)

The Gibbs free energy and its surface correction defined by the limits (2.6) exist and are given by

[TABLE] 2. (b)

Eqs. (2.2) and (2.3) hold true.

Proof.

We compute

[TABLE]

Let $\mathcal{W}_{0n}=\sum_{j\leq 0}\sum_{k\geq n+1}v(z_{j}+\cdots+z_{k})$. We note

[TABLE]

and with (4.3) deduce

[TABLE]

Now $\mathcal{W}_{0n}=O(n^{-(s-2)})\to 0$ uniformly on $(r_{\rm hc},\infty)^{\mathbb{Z}}$. By Theorem 4.4, $\mu_{\beta}(\exp(\beta[\mathcal{W}_{0}+\mathcal{W}_{n}]))=\mu_{\beta}(f(f\circ\tau^{n}))\to\mu_{\beta}(f)^{2}$ where $f=\exp(\beta\mathcal{W}_{0})$. Consequently as $n\to\infty$

[TABLE]

from which part (a) of the lemma follows. A computation analogous to Eq. (4.6) shows that for every local test function $f\in C_{b}(\mathbb{R}_{+}^{k})$,

[TABLE]

Part (b) of the lemma then follows from Theorem 4.4. ∎

5. Large deviations as $\beta\to\infty$

Here we analyze the behavior of the bulk and surface Gibbs measures $\mu_{\beta}$ and $\nu_{\beta}$ and of the energies $g(\beta)$ and $g_{\mathrm{surf}}(\beta)$. The large deviations result for the surface measure $\nu_{\beta}$ is a consequence of the eigenvalue equation from Lemma 4.1, exponential tightness, and the uniqueness of the solution to the fixed point equation in Proposition 3.9. Since the bulk measure is absolutely continuous with respect to the product measure of two independent half-infinite chains (Eq. (4.2) and Proposition 4.9(b)), we may go from the surface to the bulk measure with the help of Varadhan’s integral lemma [DZ98, Chapter 4.3]. The asymptotic behavior of $e_{\mathrm{surf}}(\beta)$ is based on the representation from Proposition 4.9(a).

5.1. A tightness estimate

The following estimate will help us prove that the infinite-volume measure $\nu_{\beta}$ is exponentially tight (see the proof of Lemma 5.3) which enters the proof of Theorem 2.4.

Lemma 5.1**.**

For all $\beta,p>0$, $N\in\mathbb{N}$, $k\in\{1,\ldots,N-1\}$, and $r\geq 0$, we have

[TABLE]

Proof.

Fix $k\in\mathbb{N}$ and $r\geq 0$. For $z=(z_{1},\ldots,z_{N-1})\in\mathbb{R}_{+}^{N-1}$ with $z_{k}\geq z_{\max}+r$ we define a new configuration $z^{\prime}$ by setting $z^{\prime}_{k}=z_{k}-r$ and leaving all other spacings unchanged. This decreases the Gibbs energy by an amount at least

[TABLE]

A change of variables thus yields

[TABLE]

and the proof of the lemma is easily concluded. ∎

5.2. Gibbs free energy in the bulk

Lemma 5.2**.**

Let $\beta\to\infty$ at fixed $p$. Then

[TABLE]

Proof of Lemma 5.2.

The relation between $g(\beta)$ and $\lambda_{0}(\beta)$ has been proven in Proposition 4.9. We proceed with an upper bound for $Q_{N}(\beta)$ and $\lambda_{0}(\beta)$. For $z=(z_{1},\ldots,z_{N-1})$, define $z^{\prime}$ by $z^{\prime}_{j}=\min(z_{\max},z_{j})$. Revisiting the proof of Lemma 3.1, we see that

[TABLE]

It follows that

[TABLE]

and

[TABLE]

whence $\beta^{-1}\log\lambda_{0}(\beta)\leq-e_{0}+O(\beta^{-1})$. For a lower bound, we let $\bar{z}\in[z_{\min},z_{\max}]^{N-1}$ be the minimizer of $\mathcal{E}_{N}$ and choose $0<\varepsilon<a-z_{\min}$ so small that by Lemma 3.3

[TABLE]

for every $z\in\times_{j=1}^{N-1}[\bar{z}_{j}-\varepsilon,\bar{z}_{j}+\varepsilon]$. We get

[TABLE]

This yields

[TABLE]

and $\beta^{-1}\log\lambda_{0}(\beta)\geq-e_{0}+O(\beta^{-1}\log\beta)$. ∎

5.3. Large deviations principles for $\nu_{\beta}$ and $\mu_{\beta}$

Here we prove Theorem 2.4.

Lemma 5.3**.**

Every sequence $\beta_{j}\to\infty$ has a subsequence along which $(\nu_{\beta_{j}})_{j\in\mathbb{N}}$ satisfies a large deviations principle with speed $\beta_{j}$ and some good rate function.

Remark*.*

If $p=p_{\beta}\to 0$, we lose exponential tightness and only know that every sequence $(\nu_{\beta_{j}})$ has a subsequence along which it satisfies a weak large deviations principle [DZ98, Lemma 4.1.23], which means that the upper bound in (2.5) is required to hold for compact sets rather than closed sets.

Proof.

The lemma is a consequence of exponential tightness. Let $n\in\mathbb{N}_{0}$. Define $K_{n}=\times_{j=1}^{\infty}[0,z_{\max}+n+j]$. $K_{n}$ is compact in the product topology. Passing to the limit $N\to\infty$ in Lemma 5.1, we find

[TABLE]

for all $k\in\mathbb{N}$ and $r\geq 0$. Therefore

[TABLE]

It follows that the family of measures $(\nu_{\beta})_{\beta\geq 1}$ is exponentially tight, i.e., for every $M>0$, we can find a compact subset $K\subset\mathbb{R}_{+}^{\mathbb{N}}$ such that $\limsup_{\beta\to\infty}\frac{1}{\beta}\log\nu_{\beta}(K^{\mathrm{c}})\leq-M$. $\mathbb{R}_{+}^{\mathbb{N}}$ endowed with the product topology is separable and metrizable and therefore has a countable base. Lemma 4.1.23 in [DZ98] applies and yields the claim. ∎

Lemma 5.4**.**

Suppose that Assumption 3 holds true and assume that along some subsequence $(\beta_{j})$ the measure $\nu_{\beta_{j}}$ satisfies a large deviations principle with good rate function $I(z_{1},z_{2},\ldots)$. Then $I$ satisfies

[TABLE]

on $\mathbb{R}_{+}^{\mathbb{N}}$. In particular, $I((z_{j})_{j\in\mathbb{N}})=\infty$ if $z_{j}\leq r_{\mathrm{hc}}$ for some $j\in\mathbb{N}$.

Proof.

Write $\beta$ instead of $\beta_{j}$. We will see that the fixed point equation for $I$ follows from the eigenvalue equation in Lemma 4.1 and the asymptotics of the principal eigenvalue provided in Lemma 5.2. According to these,

[TABLE]

for any $n\in\mathbb{N}$ where the $o(1)$-term comes from $\log\lambda_{0}^{n}(\beta)=-\beta[ne_{0}+o(1)]$ and is independent of $(z_{j})_{j\in\mathbb{N}}$.

We first show that $I$ can only be finite on $(r_{\rm hc},\infty)^{\mathbb{N}}$. Fix $n\in\mathbb{N}$ and for $\varepsilon>0$ consider the open set $O_{\varepsilon}=\{z\in\mathbb{R}^{\mathbb{N}}\mid 0<z_{n}<r_{\rm hc}+\varepsilon\}$. A repeated application of Lemma 4.1 and Lemma 5.2 give

[TABLE]

Let $-C$ be a lower bound for $-e_{0}+v(z_{\max})+\sum_{k=2}^{\infty}v(z_{1}+\cdots+z_{k})$ on $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$. Then

[TABLE]

and

[TABLE]

Hence

[TABLE]

It follows that

[TABLE]

Since $n$ was arbitrary we have shown that $I\equiv\infty$ on $\mathbb{R}_{+}^{\mathbb{N}}\setminus(r_{\rm hc},\infty)^{\mathbb{N}}$. In particular, as $\nu_{\beta}$ satisfies a large deviations principle on $\mathbb{R}_{+}^{\mathbb{N}}$ with rate function $I$, the same large deviations principle holds on $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$.

We now establish another (weak) large deviations principle on $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$. Let $K\subset(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ be a (relatively) closed set and $[\alpha,b]\subset(r_{\mathrm{hc}},\infty)$ a compact interval. Then (5.1) with $n=1$ yields

[TABLE]

Write $f_{\beta}(z_{1};K)$ for the inner integral. As $h$ is bounded from below and for every fixed $z_{1}>r_{\mathrm{hc}}$, $(z_{2},z_{3},\ldots)\mapsto h(z_{1},z_{2},\ldots)$ is continuous n $(r_{\rm hc},\infty)^{\mathbb{N}}$ with respect to the product topology, we deduce from Varadhan’s lemma [DZ98, Chapter 4.3] that

[TABLE]

for all $z_{1}\in[\alpha,b]$. Next we note that for all $(z_{j})_{j\in\mathbb{N}}\in(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$, $z^{\prime}_{1}>r_{\mathrm{hc}}$, and suitable $C>0$,

[TABLE]

For $z_{1},z^{\prime}_{1}$ bounded away from $r_{\mathrm{hc}}$ we may exploit that the derivative of $v$ is bounded and drop the first term, making $C$ larger if need be. Plugging these estimates into the definition of $f_{\beta}(z_{1},K)$, we find that for some $C_{\alpha}>0$ and all $\beta>0$,

[TABLE]

It follows that the upper bound (5.2) is uniform on compact subsets of $(r_{\mathrm{hc}},\infty)$ and

[TABLE]

A similar argument shows that for all $b>\alpha>r_{\mathrm{hc}}$ and all (relatively) open subsets $O\subset(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$,

[TABLE]

Taking monotone limits, the latter inequality is seen to extend to $\alpha=r_{\mathrm{hc}}$ and $b=\infty$. It follows that $(\nu_{\beta})$, as a family of probability measures on $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$, satisfies a weak large deviations principle with rate function $J=h_{1}-e_{0}+I(z_{2},\ldots)$. (It is indeed sufficient to consider product sets. This is easy to see for the lower bound: If $U\subset(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ is open, then for any $\varepsilon>0$ one finds $\bar{z}\in(\alpha,b)\times O\subset U$ with h(\bar{z}_{1},\bar{z}_{2},\ldots)-e_{0}+I(\bar{z}_{2},\bar{z}_{3},\ldots)-\varepsilon\leq\inf_{z\in U}\bigl{(}h(z_{1},z_{2},\ldots)-e_{0}+I(z_{2},z_{3},\ldots)\bigr{)}, from which it follows that (5.4) holds for $U$ instead of $(\alpha,b)\times O$. The upper bound for a general compact $V\subset(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ is obtained by covering, for given $\varepsilon>0$, $V\subset\bigcup_{i=1}^{N_{\varepsilon}}(\alpha_{x_{i}},b_{x_{i}})\times B_{\delta(x_{i})}(x_{i})$, where for each $x\in V$, $b_{x}>\alpha_{x}>r_{\mathrm{hc}}$ and $\delta(x)>0$ are chosen such that h(x_{1},x_{2},\ldots)-e_{0}+I(x_{2},x_{3},\ldots)-\varepsilon\leq\inf_{z\in(\alpha_{x},b_{x})\times B_{\delta(x)}(x)}\bigl{(}h(z_{1},z_{2},\ldots)-e_{0}+I(z_{2},z_{3},\ldots)\bigr{)}. This is possible since $I$ is lower semicontinuous. With the help of (5.3) we can now deduce that (5.3) holds for $V$ instead of $[\alpha,b]\times K$.)

Since $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ is a Polish space, the rate function in a weak large deviations principle is uniquely defined [DZ98, Chapter 4.1], hence $J=I$ on $(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$. To finish the proof it remains to observe that also $J=I$ on $\mathbb{R}_{+}^{\mathbb{N}}\setminus(r_{\mathrm{hc}},\infty)^{\mathbb{N}}$ because both $I$ and $h$ are equal to $\infty$ on that set. ∎

Proof of Theorem 2.4.

The large deviations principle for $\nu_{\beta}$ with good rate function $\overline{\mathcal{E}}_{\mathrm{surf}}-\min\mathcal{E}_{\mathrm{surf}}$ is an immediate consequence of Lemmas 5.3 and 5.4 and Proposition 3.9. As a consequence, $\nu^{-}_{\beta}\otimes\nu^{+}_{\beta}$ satisfies a deviations principle with good rate function $(z_{j})_{j\in\mathbb{Z}}\mapsto\overline{\mathcal{E}}_{\mathrm{surf}}(z_{1},z_{2},\ldots)+\overline{\mathcal{E}}_{\mathrm{surf}}(z_{0},z_{-1},\ldots)-2\min\mathcal{E}_{\mathrm{surf}}$ on $\mathbb{R}_{+}^{\mathbb{Z}}$ and on $[r_{\rm hc},\infty)^{\mathbb{Z}}$, The large deviations principle for $\mu_{\beta}$ thus follows from Eq. (4.2), Lemmas 4.3.4 and 4.3.6 in [DZ98], $\min\overline{\mathcal{E}}_{\mathrm{bulk}}=0$ and

[TABLE]

by Proposition 2.3, and the observation that $\mathcal{W}_{0}$ is continuous on $[r_{\rm hc},\infty)^{\mathbb{Z}}$. ∎

5.4. Surface corrections to the Gibbs free energy

Proof of Theorem 2.5.

The statements about $g(\beta)$ have already been proven in Lemma 5.2. For $g_{\mathrm{surf}}(\beta)$, we start from the formula in Proposition 4.9(a), to which we apply Lemma 5.2, Theorem 2.4 and Varadhan’s lemma. This yields

[TABLE]

But now for $(z_{j})$ with $\sum_{j\in\mathbb{Z}}(z_{j}-a)^{2}<\infty$

[TABLE]

with $e_{\mathrm{clamp}}=-pa-\sum\limits_{k=1}^{\infty}k\,v(ka)$. So

[TABLE]

6. Gaussian approximation

Here we prove Theorems 2.7 and 2.8 on the Gaussian approximation to the bulk measure $\mu_{\beta}$ when $m$ is finite. We start from a standard idea, namely perturbation theory for transfer operators [Hel02], however we need to put some work into a good choice of transfer operator as the standard symmetrized choice (6.2) does not work well. This aspect is explained in more detail in Section 6.1. Throughout this section $m$ satisfies $2\leq m<\infty$. Remember $d=m-1$.

6.1. Decomposition of the energy. Choice of transfer operator

For finite $m$, the treatment with transfer operators from Section 4.1 can be considerably simplified: instead of an operator that acts on functions of infinitely many variables, the transfer operator becomes an integral operator in $L^{2}(\mathbb{R}^{d})$ ($L^{2}$ space with respect to Lebesgue measure). There are several possible choices, corresponding each to an additive decomposition of the energy. Let $V(z_{1},\ldots,z_{d}):=\mathcal{E}_{m}(z_{1},\ldots,z_{d})$ and

[TABLE]

Let us block variables as $x_{j}=(z_{dj+1},\ldots,z_{dj+d})$. Then for $(z_{j})_{j\in\mathbb{Z}}\in\mathcal{D}_{0}^{+}$ we have

[TABLE]

with only finitely many non-zero summands. By Proposition 2.3 the sum extends to $\mathcal{D}^{+}$ by continuity. The transfer operator associated with the representation (6.1) is the integral operator with kernel $\exp(-\beta[V(x)+W(x;y)])$; it is clearly related to the $d$-th power of the transfer operator $\mathcal{L}_{\beta}$ from Section 4.1. The analysis is simpler for a symmetrized operator with kernel

[TABLE]

which has the advantage of being Hilbert-Schmidt: The pressure term present in $V(x)$ and $V(y)$ ensures that $T_{\beta}(x,y)$ decays exponentially fast when $|x|+|y|\to\infty$ so that $\int_{\mathbb{R}^{2d}}T_{\beta}(x,y)^{2}\mathrm{d}x\mathrm{d}y<\infty$. The transfer operator $T_{\beta}$ corresponds to a rewriting of  (6.1),

[TABLE]

For the analysis of the limit $\beta\to\infty$, we would like to have a transfer operator that concentrates in some sense around the optimal spacings so that we may approximate it with a Gaussian operator. When $m\geq 3$, unfortunately, the function $(x,y)\mapsto\tfrac{1}{2}V(x)+W(x;y)+\tfrac{1}{2}V(y)$ need not have its minimum at $(x,y)=(\boldsymbol{a},\boldsymbol{a})$, with $\boldsymbol{a}=(a,\ldots,a)\in\mathbb{R}^{d}$. Therefore we introduce yet another variant of the transfer operator: we look for a function $\widehat{H}(x,y)$ such that

[TABLE]

and $\widehat{H}(x,y)\geq\widehat{H}(\boldsymbol{a},\boldsymbol{a})=0$, and work with the kernel

[TABLE]

By a slight abuse of notation we use the same letter for the integral operator

[TABLE]

in $L^{2}(\mathbb{R}^{d})$. The function $\widehat{H}$ is defined as follows. Set

[TABLE]

and

[TABLE]

Remember

[TABLE]

Lemma 6.1**.**

Assume $2\leq m<\infty$, $p\in[0,p^{*})$, and $r_{\mathrm{hc}}>0$. Then:

1. (a)

For all $x,y\in(r_{\mathrm{hc}},\infty)^{d}$, we have $\widehat{H}(x,y)\geq\widehat{H}(\boldsymbol{a},\boldsymbol{a})=0$. 2. (b)

The function $g(x):=\frac{1}{2}[u(x)-u(\sigma x)]$ is bounded, and we have

[TABLE] 3. (c)

$\widehat{H}(x,y)=\widehat{H}(\sigma y,\sigma x)$* for all $x,y\in(r_{\mathrm{hc}},\infty)^{d}$.*

Proof.

One easily checks

[TABLE]

which yields

[TABLE]

For $x=y=\boldsymbol{a}$, we have $H(\boldsymbol{a},\boldsymbol{a})=w(\boldsymbol{a})$ hence $\widehat{H}(\boldsymbol{a},\boldsymbol{a})=0$. This proves part (a) of the lemma. The symmetry in part (c) is immediate from the reversal symmetry of $\mathcal{E}_{\mathrm{bulk}}$. For (b), we note that

[TABLE]

the formula for $\widehat{H}$ follows. Because of

[TABLE]

and $V(\sigma x)=V(x)$, $C:=\sup_{(x,y)\in(r_{\mathrm{hc}},\infty)^{2d}}|W(x,y)-W(\sigma x,y)|<\infty$, we have

[TABLE]

The roles of $x$ and $\sigma x$ can be exchanged, hence $u(x)-u(\sigma x)$ is bounded. ∎

6.2. Some properties of the transfer operator

Lemma 6.2**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and and $r_{\mathrm{hc}}>0$. Then:

1. (a)

The kernels $K_{\beta}$ and $T_{\beta}$ are related as follows:

[TABLE] 2. (b)

The operator $K_{\beta}$ is a Hilbert-Schmidt operator in $L^{2}(\mathbb{R}^{d})$, and the kernel has the symmetry $K_{\beta}(x,y)=K_{\beta}(\sigma y,\sigma x)$.

The lemma follows from Lemma 6.1, the elementary proofs are omitted.

By the Krein-Rutman theorem [KR48], [Dei85, Chapter 6], the operator norm $||K_{\beta}||=:\Lambda_{0}(\beta)$ is a simple eigenvalue of $K_{\beta}$, the associated eigenfunction $\phi_{\beta}$ can be chosen strictly positive on $(r_{\mathrm{hc}},\infty)^{d}$, and the other eigenvalues of $K_{\beta}$ have absolute value strictly smaller than $\Lambda_{0}(\beta)$, i.e.,

[TABLE]

By Lemma 6.2(b), the function $\phi_{\beta}\circ\sigma$ is a left eigenfunction of $K_{\beta}$:

[TABLE]

Let $\Pi_{\beta}$ be the rank-one projection in $L^{2}(\mathbb{R}^{d})$ given by

[TABLE]

Then $K_{\beta}\Pi_{\beta}=\Lambda_{0}(\beta)\Pi_{\beta}=\Pi_{\beta}K_{\beta}$ and an induction over $n\in\mathbb{N}$ shows

[TABLE]

Since $\Lambda_{1}(\beta)$ is nothing else but the spectral radius of $K_{\beta}-\Lambda_{0}(\beta)\Pi_{\beta}$, it follows that

[TABLE]

The spectral properties of $K_{\beta}$ are related to the Gibbs free energy and the Gibbs measure as follows.

Lemma 6.3**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and $r_{\mathrm{hc}}>0$. Then:

1. (a)

The Gibbs free energy is given by $g(\beta)=e_{0}-\frac{1}{\beta d}\log\Lambda_{0}(\beta)$. 2. (b)

The $nd$-dimensional marginals of the bulk Gibbs measure $\mu_{\beta}$ have probability density function

[TABLE]

with $c=\langle\phi_{\beta},\phi_{\beta}\circ\sigma\rangle$. 3. (c)

For all $\varepsilon>0$ and all bounded $f,g:\mathbb{R}^{d}\to\mathbb{R}$, writing f_{0}\bigl{(}(z_{j})_{j\in\mathbb{Z}}\bigr{)}:=f(z_{0},\ldots,z_{d-1}) and g_{n}\bigl{(}(z_{j})_{j\in\mathbb{Z}}\bigr{)}:=g(z_{nj},\ldots,z_{nj+d-1}), we have

[TABLE]

with some constant $C_{\varepsilon}(\beta)$ that does not depend on $f$, $g$, or $n$. If $m=2$, we can pick $\varepsilon=0$ and $C_{0}=1$.

Proof of Lemma 6.3.

For $N=nd+1$, the partition function $Q_{N}(\beta)$ is given by

[TABLE]

For the second identity we have used Lemma 6.2(a). The function $u-u\circ\sigma$ is bounded by Lemma 6.1(b) and $\exp(-\beta V)$ is integrable because $V(z_{1},\ldots,z_{d})=\mathcal{E}_{m}(z_{1},\ldots,z_{d})$ grows linearly when $|z_{j}|\to\infty$. Therefore $F_{\beta}:=\exp(-\beta V/2-\beta[u-u\circ\sigma]/2)$ and $F_{\beta}\circ\sigma$ are in $L^{2}(\mathbb{R}^{d})$, and as $n\to\infty$,

[TABLE]

It follows that

[TABLE]

which proves part (a) of the lemma. The standard proof of part (b) is omitted (compare [Hel02, Chapter 4]). For (c), we use the formula for the $(n+1)d$- dimensional marginal provided by (b). Let us choose multiplicative constants in such a way that $c=\langle\phi_{\beta},\phi_{\beta}\circ\sigma\rangle=1$. Then

[TABLE]

Eq. (6.4) yields

[TABLE]

where $||\cdot||$ refers to the $L^{2}$-norm for functions and the operator norm for the operator. We further bound $||g\phi_{\beta}||\leq||g||_{\infty}||\phi_{\beta}||$ and $||f(\phi_{\beta}\circ\sigma)||\leq||f||_{\infty}||\phi_{\beta}||$ and conclude with (6.5). If $m=2$, the operators are symmetric, hence the operator norm is the same as the spectral radius and the estimates simplify accordingly. ∎

Remark* (Associated Markov chain).*

Define the kernel

[TABLE]

on $(r_{\mathrm{hc}},\infty)^{d}$. Then $P_{\beta}$ is a Markov kernel with invariant measure $\rho_{\beta}(x)\mathrm{d}x$ where

[TABLE]

If in the bulk Gibbs measure $\mu_{\beta}$ we group spacing in blocks as $x_{n}=(z_{dn},\ldots,z_{dn+d-1})$, we obtain a probability measure on $(r_{\mathrm{hc}},\infty)^{d}$. This measure is exactly the distribution of the two-sided stationary Markov chain $(X_{j})_{j\in\mathbb{Z}}$ with state space $\mathbb{R}^{d}$, transition kernel $P_{\beta}$, and initial law $\mathcal{L}(X_{0})=\rho_{\beta}(x)\mathrm{d}x$.

6.3. Gaussian transfer operator

Here we introduce the Gaussian counterpart to the transfer operator $K_{\beta}$ and study its spectral properties. We start from the quadratic approximation to the bulk energy $\mathcal{E}_{\mathrm{bulk}}$. The differentiability of $\mathcal{E}_{\mathrm{bulk}}$ in a neighborhood of the constant sequence $z_{j}\equiv a$ is checked in Lemma 6.11 below, for the definition of the Gaussian transfer operator we only need the infinite matrix of partial derivatives at $(\ldots,a,a,\ldots)$.

In the following we block variables as $x_{j}=(z_{dj},\ldots,z_{dj+d-1})$ for $z=(z_{j})_{j\in\mathbb{Z}}$ and $\xi_{j}=(\zeta_{dj},\ldots,\zeta_{dj+d-1})$ for $\zeta=(\zeta_{j})_{j\in\mathbb{Z}}$. Remember the decomposition (6.1). Set $\boldsymbol{a}=(a,\ldots,a)\in\mathbb{R}^{d}$ and define the $d\times d$ matrices

[TABLE]

We note the following relations:

[TABLE]

The Hessian $\mathrm{D}^{2}\mathcal{E}_{\mathrm{bulk}}$ at $(\ldots,a,a,\ldots)$ is a doubly infinite, band-diagonal matrix with block form

[TABLE]

Note that Lemma 3.3 implies that $\mathrm{D}^{2}\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)$ is positive definite. We look for a quadratic form $\mathcal{Q}(x,y)$ on $\mathbb{R}^{2d}$ that is positive-definite and satisfies

[TABLE]

One candidate choice could be

[TABLE]

but it is not easily related to $\widehat{H}(x,y)$. We make a different choice which mimicks the definition of $\widehat{H}(x,y)$ and show later that this amounts to picking the Hessian of $\widehat{H}(x,y)$ (see Lemma 6.12 below).

We introduce the quadratic counterparts to the functions $H(x,y)$, $w(x)$, and $\widehat{H}(x,y)$ from Section 6.2. Remember the bulk Hessian from (6.9). Since it is positive-definite, there exist uniquely defined positive-definite matrices $M\in\mathbb{R}^{2d\times 2d}$ and $N\in\mathbb{R}^{d\times d}$ such that

[TABLE]

for all $x,y\in\mathbb{R}^{d}$. The quadratic forms associated with $M$ and $N$ are the Gaussian counterparts to the functions $H(x,y)$ and $w(x)$, respectively. Finally set

[TABLE]

and

[TABLE]

We will see in the proof of Lemma 6.12 that $M$, $N$ and $\widehat{M}$ are the Hessians of $H$ at $(\boldsymbol{a},\boldsymbol{a})$, $w$ at $\boldsymbol{a}$ and $\widehat{H}$ at $(\boldsymbol{a},\boldsymbol{a})$, respectively. The relation between $\mathcal{Q}$ and $\widehat{\mathcal{Q}}(x,y)$ is clarified in Lemma 6.7 below. We are going to work with the kernel

[TABLE]

and the associated integral operator $(G_{\beta}f)(x)=\int_{\mathbb{R}^{d}}G_{\beta}(x,y)f(y)\mathrm{d}y$. In Section 6.4 we show that $G_{\beta}$ is a good approximation for $K_{\beta}$, here we study the operator $G_{\beta}$ on its own. Clearly it is enough to understand the integral operator $G$ with kernel

[TABLE]

since $G$ and $G_{\beta}$ are related by the change of variables $x\mapsto\sqrt{\beta}(x-\boldsymbol{a})$, see Eq. (6.21) below.

Lemma 6.4**.**

Assume $2\leq m<\infty$, $p\in[0,p^{*})$. Then the quadratic form $\widehat{\mathcal{Q}}$ is positive-definite: $\widehat{\mathcal{Q}}(x,y)\geq\varepsilon(|x|^{2}+|y|^{2})$ for some $\varepsilon>0$ and all $(x,y)\in\mathbb{R}^{2d}$.

Proof.

First we show that $\widehat{M}$ is positive semi-definite, by an argument similar to Lemma 6.2(a). Define

[TABLE]

Clearly

[TABLE]

hence

[TABLE]

for all $(x,y)\in\mathbb{R}^{d}\times\mathbb{R}^{d}$ and $\widehat{M}$ is positive semi-definite. Next let $(x_{0},y_{0})\in\mathbb{R}^{d}\times\mathbb{R}^{d}$ be a zero of the quadratic form associated with $\widehat{M}$. Then by (6.13), the function $y\mapsto F(x_{0},y)$ must be minimal at $y=y_{0}$, hence $\nabla_{y}F(x_{0},y)=0$. Similarly, the function $y\mapsto F(x,y_{0})$ must be minimal at $x=x_{0}$, hence $\nabla_{x}F(x_{0},y_{0})=0$. Thus $(x_{0},y_{0})$ is a critical point of $F$. But $F$ is strictly convex because $M$ is positive-definite, therefore the critical point $(x_{0},y_{0})$ is a global minimizer of $F$ which yields $(x_{0},y_{0})=0$. It follows that $\widehat{M}$ is positive-definite. ∎

It follows from Lemma 6.4 that $\int_{\mathbb{R}^{2d}}G(x,y)^{2}\mathrm{d}x\mathrm{d}y<\infty$, hence $G$ is Hilbert-Schmidt with strictly positive integral kernel and Krein-Rutman theorem is applicable. So we may ask for its principal eigenvalue and eigenvector and its spectral gap. It is natural to look for a Gaussian eigenfunction.

Lemma 6.5**.**

Let $F$ be a positive-definite, symmetric $d\times d$ matrix. Then the following two statements are equivalent:

1. (i)

$\phi(x):=\exp(-\tfrac{1}{2}\langle x,Fx\rangle)$* is an eigenfunction of $G$.* 2. (ii)

The function $x\mapsto\langle x,Fx\rangle$ satisfies the quadratic Bellman equation

[TABLE]

Proof.

The proof is by a straightforward completion of squares: write

[TABLE]

with $d\times d$ -matrices $\widehat{M}_{j}$. The diagonal blocks $\widehat{M}_{1}$ and $\widehat{M}_{3}$ are positive-definite because $\widehat{M}$ is positive-definite, therefore $\widehat{M}_{3}+F$ is positive-definite as well. Then

[TABLE]

It follows that

[TABLE]

and

[TABLE]

Therefore (i) and (ii) hold true if and only if $F$ solves

[TABLE]

In particular, (i) and (ii) are equivalent. ∎

In Lemma 6.7 below we check that $M$ is of the form

[TABLE]

for some positive-definite $d\times d$ matrix $C$.

Lemma 6.6**.**

The principal eigenvalue of $G$ is $\sqrt{(2\pi)^{d}/\det C}$ and the principal eigenfunction is $\exp(-\tfrac{1}{2}\langle x,\tfrac{1}{2}Nx\rangle)$ (up to scalar multiples).

Proof.

A close look at our definitions shows that $F:=\frac{1}{2}N$ solves (6.14) (it is positive-definite because $N$ is). Indeed, by the definition of $\widehat{\mathcal{Q}}$, $\widehat{M}$, we have

[TABLE]

Therefore, by Lemma 6.5, the function $\phi(x)=\exp(-\frac{1}{4}\langle x,Nx\rangle)$ is an eigenfunction of $G$. The matrix $\widehat{M}_{3}+F$ in (6.15) is equal to $(C-\tfrac{1}{2}N)+F=C$, and we find that the principal eigenvalue of $G$ is $\sqrt{(2\pi)^{d}/\det C}$. ∎

In order to identify the block $C$ in (6.16), we introduce the quadratic analogue to the function $u(x)$. Let $A$ and $B$ be the $d\times d$ matrices from (6.7) and $A_{1}:=V_{xx}(\boldsymbol{a})+W_{xx}(\boldsymbol{a},\boldsymbol{a})$. The infinite matrix $(\partial_{i}\partial_{j}\mathcal{E}_{\mathrm{surf}}(a,a,\ldots))_{i,j\in\mathbb{N}}$ is band-diagonal with block structure

[TABLE]

The matrix differs from the bulk Hessian (6.9) by the upper left corner $A_{1}$: we have

[TABLE]

By a reasoning similar to Lemma 3.3, the Hessian of $\mathcal{E}_{\mathrm{surf}}$ is positive-definite. Therefore there is a uniquely defined positive-definite $d\times d$-matrix $D$ such that

[TABLE]

for all $x\in\mathbb{R}^{d}$. (Analogous arguments as in the proof of Lemma 6.12 show that $D$ is the Hessian of $u$ at $\boldsymbol{a}$.) Set

[TABLE]

and

[TABLE]

(remember the symmetries (6.8)).

Lemma 6.7**.**

The matrix $C$ solves

[TABLE]

and Eq. (6.16) holds true. Moreover

[TABLE]

Proof.

Clearly

[TABLE]

hence

[TABLE]

by a completion of squares similar to the proof of Lemma 6.5. We add $W_{yy}(\boldsymbol{a},\boldsymbol{a})$ to both sides, remember (6.17), and obtain the equation for $C$. It is easy to see that

[TABLE]

which proves (6.16). Furthermore,

[TABLE]

hence,

[TABLE]

Let us check that the two expressions for $N$ are indeed identical, and that $\sigma N\sigma=N$. Combining with (6.17) and (6.19), the two expressions for $N$ become

[TABLE]

and

[TABLE]

The two expressions are indeed equal, and from the end formula and (6.8) we read off that $\sigma N\sigma=N$. Actually

[TABLE]

which is the analogue of $w(x)=u(x)+u(\sigma x)-V(x)$.

Now we compute $\widehat{M}$. The off-diagonal blocks of $\widehat{M}$ are the same as those of $M$. The upper left diagonal block is

[TABLE]

A similar computation yields the lower right block. Altogether we find

[TABLE]

and the lemma follows. ∎

Finally we come back to the $\beta$-dependent operator $G_{\beta}$.

Proposition 6.8**.**

Assume $2\leq m<\infty$ and $p\in[0,p^{*})$. The principal eigenvalue of $G_{\beta}$ is

[TABLE]

and the normalized, positive principal eigenfunction is

[TABLE]

Proof.

Let $U_{\beta}:L^{2}(\mathbb{R}^{d})\to L^{2}(\mathbb{R}^{d})$ be the unitary operator given by

[TABLE]

We have

[TABLE]

hence

[TABLE]

and the principal eigenvalue and eigenfunction of $G_{\beta}$ are obtained from those of $G$ in Lemma 6.6 by straightforward transformations. ∎

Remark*.*

When $m=2$, all eigenvalues and eigenfunctions of $G$ (hence $G_{\beta}$) can be computed explicitly, and the eigenfunctions are expressed with Hermite polynomials. See [Hel02, Section 5.2] on the harmonic Kac operator.

6.4. Perturbation theory

Remember the unitary operator $U_{\beta}$ from (6.20) and the relation $G_{\beta}=\beta^{-d/2}U_{\beta}^{*}GU_{\beta}$. The main technical result of this section is the following.

Proposition 6.9**.**

Assume $2\leq m<\infty$, $p\in(0,p^{*})$, and $r_{\mathrm{hc}}>0$. We have $||\beta^{d/2}(K_{\beta}-G_{\beta})||=||G-\beta^{d/2}U_{\beta}K_{\beta}U_{\beta}^{*}||\to 0$ as $\beta\to\infty$.

Before we come to the proof of the proposition, we state a corollary on the principal eigenvalue and eigenfunction. Remember the quantities $\Lambda_{0}(\beta)$, $\Lambda_{1}(\beta)$, $\phi_{\beta}$ defined before Lemma 6.3. We choose multiplicative constants so that $||\phi_{\beta}||=1$. Let $\lambda_{j}^{\mathrm{Gauss}}$, $j\in\mathbb{N}_{0}$, be an enumeration of the eigenvalues of $G$ with $\lambda_{0}^{\mathrm{Gauss}}=||G||$ and

[TABLE]

Corollary 6.10**.**

Under the assumptions of Proposition 6.9: Let $\Lambda_{0}^{\mathrm{Gauss}}(\beta)$ and $\phi^{\mathrm{Gauss}}_{\beta}(x)$ be as in Proposition 6.8. Then as $\beta\to\infty$,

[TABLE]

and

[TABLE]

The corollary follows from Proposition 6.9 and standard perturbation theory for compact operators [RS78]. The proof of Proposition 6.9 builds on several lemmas. First we show that $\mathcal{E}_{\mathrm{bulk}}$ is $C^{2}$ in a neighborhood of its global minimizer.

Lemma 6.11**.**

The mapping $\mathcal{E}_{\mathrm{bulk}}$ is $C^{2}$ in some open neighborhood in $\mathcal{D}^{+}$ of the constant sequence $(\ldots,a,a,\ldots)$.

Proof.

Note that

[TABLE]

defines a $C^{2}$ function in a neighborhood of $(a,\ldots,a)\in\mathbb{R}^{d}\times\mathbb{R}^{d}$ which vanishes for $(z_{1},\ldots,z_{2d})=(a,\ldots,a)$. Moreover, using that $(\ldots,a,a,\ldots)$ minimizes $\mathcal{E}_{\mathrm{bulk}}$ on $\mathcal{D}_{0}^{+}$ and so $\partial_{x_{j}}\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)=0$, we see that also

[TABLE]

For all $z\in\mathcal{D}_{0}^{+}$ the derivative of $\mathcal{E}_{\mathrm{bulk}}$ at $z$ is given by

[TABLE]

for all $\zeta\in\ell^{2}(\mathbb{Z})$ with $\zeta_{j}=0$ for all but finitely many $j$. So

[TABLE]

Since

[TABLE]

for $z,z^{\prime}\in\mathcal{D}^{+}$ in a neighborhood of $(\ldots,a,a,\ldots)$ with a uniform constant $C$, the right hand side of (6.22) extends to a uniformly continuous function there. Writing

[TABLE]

for $z,z^{\prime}\in\mathcal{D}_{0}^{+}$, a standard approximation argument shows that indeed $\mathcal{E}_{\mathrm{bulk}}$ is $C^{1}$ in a neighborhood of $(\ldots,a,a,\ldots)$ also in $\mathcal{D}^{+}$ with $\mathrm{D}\mathcal{E}_{\mathrm{bulk}}$ given by (6.22). In fact, $\mathcal{E}_{\mathrm{bulk}}$ is even $C^{2}$ on a neighborhood of $(\ldots,a,a,\ldots)$ in $\mathcal{D}^{+}$ and

[TABLE]

This follows similarly as above by extending the derivative of $\mathrm{D}\mathcal{E}_{\mathrm{bulk}}$, where we now use that the mappings $\mathbb{R}^{d}\times\mathbb{R}^{d}\times\mathbb{R}^{d}\to\mathbb{R}$, $(x,x^{\prime},x^{\prime\prime})\mapsto V_{xx}(x^{\prime})+W_{xx}(x^{\prime},x^{\prime\prime})+W_{yy}(x,x^{\prime})$ and $\mathbb{R}^{d}\times\mathbb{R}^{d}\to\mathbb{R}$, $(x,x^{\prime})\mapsto W_{xy}(x,x^{\prime})$ are uniformly continuous in a neighborhood of $x=x^{\prime}=x^{\prime\prime}=(a,\ldots,a)$ and so $\mathrm{D}^{2}\mathcal{E}_{\mathrm{bulk}}$ extends to a continuous mapping from a neighborhood of $(\ldots,a,a,\ldots)$ to $L(\ell^{2}(\mathbb{Z}))$ (the space of bounded linear operators on $\ell^{2}(\mathbb{Z})$) given by (6.23). ∎

Next we show that $\widehat{M}$ is in fact the Hessian of $\widehat{H}$.

Lemma 6.12**.**

Assume $2\leq m<\infty$, $p\in[0,p^{*})$, and $r_{\mathrm{hc}}>0$. We have $\widehat{H}(x,y)\geq\widehat{H}(\boldsymbol{a},\boldsymbol{a})=0$ for all $x,y\in\mathbb{R}_{+}^{d}$, moreover as $x,y\to\boldsymbol{a}$,

[TABLE]

The lemma leaves open whether $(\boldsymbol{a},\boldsymbol{a})$ is the unique global minimizer of $\widehat{H}$.

Proof.

The first part of the lemma has already been proven in Lemma 6.2(a). With $M\in\mathbb{R}^{2d\times 2d}$, $N\in\mathbb{R}^{d\times d}$ as in (6.10) and (6.11) we let $\widehat{M}$ as in (6.12). It remains to show that $\mathrm{D}^{2}\widehat{H}(\boldsymbol{a},\boldsymbol{a})=\widehat{M}$. Since, for a suitable $\varepsilon>0$, $\mathcal{E}_{\mathrm{bulk}}$ is convex on $\mathcal{D}^{+}\cap[z_{\min},z_{\max}+\varepsilon]^{\mathbb{Z}}$, see (the proof of) Proposition 2.3, Lemma 3.12 shows that there is a unique function on a neighborhood of $(\boldsymbol{a},\boldsymbol{a})$ in $\mathbb{R}^{d}\times\mathbb{R}^{d}$ with values in $\mathbb{R}^{-\mathbb{N}}\times\mathbb{R}^{\mathbb{N}}$, $(x,y)\mapsto\tilde{z}=(z_{-},z_{+})=(z_{-}(x,y),z_{+}(x,y))$ such that

[TABLE]

As $\mathrm{D}^{2}\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)$ is positive definite, the implicit function theorem shows that this mapping is $C^{1}$ and satisfies

[TABLE]

as well as

[TABLE]

The latter identity implies

[TABLE]

so that $H$ is indeed $C^{2}$ near $(\ldots,a,a,\ldots)$ and

[TABLE]

In particular, since $\tilde{z}(\boldsymbol{a},\boldsymbol{a})=(\ldots,a,a,\ldots)$,

[TABLE]

The same analysis applied to the quadratic approximation $\ell^{2}(\mathbb{Z})\to\mathbb{R}$, $z\mapsto\frac{1}{2}\langle z,\mathrm{D}^{2}\mathcal{E}_{\mathrm{bulk}}(\ldots,a,a,\ldots)z\rangle$ leads to

[TABLE]

too. So we have $\mathrm{D}^{2}H(\boldsymbol{a},\boldsymbol{a})=M$. A completely analogous reasoning gives $\mathrm{D}^{2}w(a,\ldots,a)=N$ and it follows that $\mathrm{D}^{2}\widehat{H}(\boldsymbol{a},\boldsymbol{a})=\widehat{M}$. ∎

Lemma 6.13**.**

Assume $2\leq m<\infty$. For some $c_{2}>0$ and all $(z_{1},\ldots,z_{2d})\in(r_{\mathrm{hc}},\infty)^{2d}$,

[TABLE]

Proof.

Since the pair potential $v$ is bounded from below, we have for some constant $c>0$

[TABLE]

In combination with Lemma 6.1 this yields the claim. ∎

In order to estimate $||K_{\beta}-G_{\beta}||$, we split the configuration space into a neighborhood $\mathcal{A}\supset B_{\delta}(\boldsymbol{a})$ of $\boldsymbol{a}$ and its complement $\mathcal{B}=\mathbb{R}^{d}\setminus\mathcal{A}$ and treat blocks separately. For $U\subset\mathbb{R}^{d}$, we write $\mathbf{1}_{U}$ for the multiplication operator with the indicator function ${\mathchoice{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.0mu\mathrm{l}}{1\mskip-4.5mu\mathrm{l}}{1\mskip-5.0mu\mathrm{l}}}_{U}$.

Lemma 6.14**.**

Suppose that $\mathcal{A}\subset\mathbb{R}^{d}$ is compact, contains an open neighborhood of $\boldsymbol{a}$, and is such that $\widehat{H}(x,y)>0$ for all $(x,y)\in\mathcal{A}\times\mathcal{A}\setminus\{(\boldsymbol{a},\boldsymbol{a})\}$. Then

[TABLE]

Proof.

By Lemma 6.12, for every $\varepsilon>0$, there is a $\delta>0$ such that for all $s,t\in\mathbb{R}^{d}$ with $|s|\leq\delta$ and $|t|\leq\delta$, we have

[TABLE]

Choosing $\delta>0$ small enough we may assume without loss of generality that $B_{\delta}(\boldsymbol{a})\subset\mathcal{A}$. We estimate

[TABLE]

for some $k>0$. On $\mathcal{A}^{2}\setminus B_{\delta}(\boldsymbol{a})^{2}$, the function $\widehat{H}$ stays bounded away from [math], therefore

[TABLE]

A similar estimate clearly holds true for $G_{\beta}$ as well. Hence

[TABLE]

This holds true for every $\varepsilon>0$, so the left-hand side converges to zero. Since operator norms are bounded by Hilbert-Schmidt norms, the lemma follows. ∎

Lemma 6.15**.**

Assume that $\mathcal{B}\subset\mathbb{R}^{d}$ is such that $\mathrm{dist}(\boldsymbol{a},\mathcal{B})>0$ and $\mathcal{B}$ is invariant under reversals, $\sigma(\mathcal{B})=\mathcal{B}$. Then $||\mathbf{1}_{\mathcal{B}}K_{\beta}\mathbf{1}_{\mathcal{B}}||=O({\operatorname{e}}^{-\beta\delta})\to 0$.

Proof.

We may view $K_{\beta}^{\mathcal{B}}=\mathbf{1}_{\mathcal{B}}K_{\beta}\mathbf{1}_{\mathcal{B}}$ as an operator in $L^{2}(\mathcal{B},\mathrm{d}x)$. The Krein-Rutman theorem is applicable and shows that $\lambda=||K_{\beta}^{\mathcal{B}}||$ is a simple eigenvalue and there exists an eigenfunction $\psi$ that is strictly positive on $\mathcal{B}\cap(r_{\mathrm{hc}},\infty)^{d}$. Because of the symmetry $\widehat{H}(\sigma y,\sigma x)=\widehat{H}(x,y)$, the function $\psi\circ\sigma$ is a left eigenfunction. Moreover for all $f,g\in L^{2}(\mathcal{B},\mathrm{d}x)$, we have

[TABLE]

so for all strictly positive functions $f,g\in L^{2}(\mathcal{B},\mathrm{d}x)$,

[TABLE]

We choose $f(y)=\exp(-\beta\widehat{H}(\boldsymbol{a},y))$ and $g(x)=\exp(-\beta\widehat{H}(x,\boldsymbol{a}))$. The scalar product becomes

[TABLE]

with $x_{0}=x_{n+2}=\boldsymbol{a}$. By Lemma 6.1(b) , remembering $u(\boldsymbol{a})=0$, we have

[TABLE]

Define $(z_{1},\ldots,z_{(n+1)d})=(x_{1},\ldots,x_{n+1})$ and for $j\in\mathbb{Z}\setminus\{1,\ldots,(n+1)d\}$, $z_{j}=a$. Then we recognize

[TABLE]

where the constant depends on $e_{0}$, $d$, and $V(a)$ alone. As $z_{1},\ldots,z_{(n+1)d}$ stay bounded away from $a$, we obtain

[TABLE]

for some $\delta,c>0$ and all $n\in\mathbb{N}$ and $x_{1},\ldots,x_{n+1}\in\mathcal{B}$. It follows that $||K_{\beta}^{\mathcal{B}}||=\lambda\leq{\operatorname{e}}^{-\beta\delta}$. ∎

Lemma 6.16**.**

Suppose that $\mathcal{A}\subset\mathbb{R}^{d}$ and $\mathcal{B}=\mathbb{R}^{d}\setminus\mathcal{A}$ are such that

[TABLE]

for some $\delta>0$ and all $x\in\mathcal{A}$, $y\in\mathcal{B}$. Assume also that $\mathcal{A}$ is invariant under reversals, $\sigma(\mathcal{A})=\mathcal{A}$. Then

[TABLE]

Proof.

Revisiting the proof of Lemma 6.1, we see that

[TABLE]

Eqs. (6.25), (6.3) and (6.24) show that $\widehat{H}(x,y)\geq\delta/2$ for all $x\in\mathcal{A}$ and $y\in\mathcal{B}$. This estimate together with the growth estimate from Lemma 6.13 shows

[TABLE]

hence $||\mathbf{1}_{\mathcal{A}}K_{\beta}\mathbf{1}_{\mathcal{B}}||\to 0$. The estimate on $||\mathbf{1}_{\mathcal{B}}K_{\beta}\mathbf{1}_{\mathcal{A}}||$ follows from the symmetry $K_{\beta}(\sigma y,\sigma x)=K_{\beta}(x,y)$. ∎

Proof of Proposition 6.9.

Let $\varepsilon>0$, $\mathcal{A}_{\varepsilon}:=[z_{\min},z_{\max}+\varepsilon]^{d}$, and $\mathcal{B}=\mathbb{R}^{d}\setminus\mathcal{A}$. The sets $\mathcal{A}$ and $\mathcal{B}$ are clearly invariant under reversals, moreover $z_{\min}<a\leq z_{\max}$ by Theorem 2.1(b), so $a$ is in the interior of $\mathcal{A}$ and bounded away from $\mathcal{B}$. Thus $\mathcal{A}$ and $\mathcal{B}$ satisfy the assumptions of Lemmas 6.14 and 6.15. By Lemma 3.11, they also satisfy the condition (6.24) from Lemma 6.16. By the triangle inequality,

[TABLE]

The first term on the right-hand side, multiplied by $\beta^{d/2}$, goes to zero by Lemma 6.14. For the second term, we estimate

[TABLE]

and conclude from Lemmas 6.15 and 6.16 that d $\beta^{d/2}||K_{\beta}-1_{\mathcal{A}}K_{\beta}\mathbf{1}_{\mathcal{A}}||\to 0$. Bounding Hilbert-Schmidt norms, it is straightforward to check that $||\beta^{d/2}(G_{\beta}-\mathbf{1}_{\mathcal{A}}G_{\beta}\mathbf{1}_{\mathcal{A}})||\to 0$ as well, and the proof is complete. ∎

6.5. Proof of Theorems 2.7, 2.8 and  2.11

Proof of Theorem 2.8.

Combining Lemma 6.3(a) and Corollary 6.10, we obtain

[TABLE]

Proof of Theorem 2.11.

The theorem is an immediate consequence of Lemma 6.3(c) and Corollary 6.10. ∎

For the proof of Theorem 2.7, we first express the marginals of $\mu^{\mathrm{Gauss}}$ in terms of the matrices $A$ and $B$ from Eq. (6.7) and the matrix $C$ from (6.18). We group variables in blocks $x_{j}\in\mathbb{R}^{d}$ as usual and view $\mu^{\mathrm{Gauss}}$ as a measure on $(\mathbb{R}^{d})^{\mathbb{Z}}$.

Proposition 6.17**.**

Under the assumptions of Theorem 2.7, the distributions of $x_{0}=(z_{0},\ldots,z_{d-1})$, $(x_{0},x_{1})$, and $(x_{0},\ldots,x_{n})$ ($n\geq 2$) under $\mu^{\mathrm{Gauss}}$ have probability density functions proportional to

1. (a)

$\exp(-\frac{1}{2}\beta\langle x_{0},(\sigma C\sigma-BC^{-1}B^{T})x_{0}\rangle)$, 2. (b)

$\exp(-\tfrac{1}{2}\beta[\langle\sigma x_{0},C\sigma x_{0}\rangle-2\langle x_{0},Bx_{1}\rangle+\langle x_{1},Cx_{1}\rangle])$, 3. (c)

$\exp(-\tfrac{1}{2}(\langle\sigma x_{0},(C-\tfrac{1}{2}A)\sigma x_{0}\rangle+\sum_{i=0}^{n-1}\mathcal{Q}(x_{i},x_{i+1})+\langle x_{n-1},(C-\tfrac{1}{2}A)x_{n-1}\rangle))$**

respectively.

Proof.

We recall a standard fact on marginals of multivariate Gaussians and Schur complements. Suppose we are given a positive-definite $(n+k)\times(n+k)$-matrix in block form

[TABLE]

where $\mathcal{H}_{1},\mathcal{H}_{2},\mathcal{H}_{3}$ are $n\times n$, $n\times k$ and $k\times k$ matrices, respectively. Think of $\mathcal{H}$ as the Hessian of the energy. Consider the Gaussian measure on $\mathbb{R}^{n+k}$ with covariance matrix $\mathcal{H}^{-1}$ and probability density function

[TABLE]

Then for all $x\in\mathbb{R}^{n}$,

[TABLE]

with $\mathcal{M}=\mathcal{H}_{1}-\mathcal{H}_{2}\mathcal{H}_{3}^{-1}\mathcal{H}_{2}^{T}$ the Schur complement of $\mathcal{H}_{3}$ in $\mathcal{H}$. The inverse $\mathcal{M}^{-1}$ is equal to the upper left block of $\mathcal{H}^{-1}$. Another characterization is provided by a completion of squares, similar to the proof of Lemma 6.5: we have

[TABLE]

Now let $\mathcal{H}=(\mathcal{H}_{ij})_{i,j\in\mathbb{Z}}$ be the Hessian of $\mathcal{E}_{\mathrm{bulk}}$ at $(\ldots,a,a,\ldots)$. By definition of $\mu^{\mathrm{Gauss}}$, the distribution of $(z_{1},\ldots,z_{n})$ is Gaussian with mean zero and covariance matrix $(\mathcal{H}^{-1})_{i,j=1,\ldots,n}$. Let $\mathcal{M}=(\mathcal{M}_{ij})_{0\leq i,j\leq n-1}$ be the $n\times n$-matrix defined by $\mathcal{M}^{-1}=(\mathcal{H}^{-1})_{0\leq i,j\leq n-1}$. It is not difficult to check that the considerations above generalize to the infinite matrices at hand, hence for all $z_{0},\ldots,z_{n-1}\in\mathbb{R}$,

[TABLE]

Eq. (6.27) provides a variational description of the covariance matrix $\mathcal{M}^{-1}$ of the $n$-dimensional marginal of $\mu^{\mathrm{Gauss}}$. For $n=2d=2(m-1)$, with $x_{0}=(z_{0},\ldots,z_{d-1})$ and $x_{1}=(z_{d},\ldots,z_{2d-1})$, Eq. (6.27) shows $\mathcal{M}=M$, by the definition (6.10) of $M$. Combining with (6.16) we get

[TABLE]

This proves part (b) of the lemma. The proof of (c) is similar. Part (a) follows from (b) and a relation similar to (6.26). ∎

Proof of Theorem 2.7.

It is enough to treat the $nd$-dimensional marginals with $n\geq 2$. Let $\phi_{\beta}$ be the principal eigenfunction of $K_{\beta}$, with multiplicative constant chosen so that $\langle\phi_{\beta}\circ\sigma,\phi_{\beta}\rangle=1$. Set $\tilde{\phi}_{\beta}(x):=U_{\beta}\phi_{\beta}(x)=\beta^{-1/4}\phi_{\beta}(\boldsymbol{a}+\beta^{-1/2}x)$ and

[TABLE]

By Lemma 6.3, the probability density $\rho_{nd}^{\scriptscriptstyle{({\beta}})}$ for $(x_{1},\ldots,x_{n})\in\mathbb{R}^{nd}$ satisfies

[TABLE]

By Proposition 6.17, the analogous representation for the Gaussian density $\rho_{nd}^{\mathrm{Gauss}}$ is

[TABLE]

with $\tilde{G}(x,y)=(\lambda_{0}^{\mathrm{Gauss}})G(x,y)$ and $\phi^{\mathrm{Gauss}}(x)\propto\exp(-\frac{1}{2}\langle x,\frac{1}{2}Nx\rangle)$ the principal eigenfunction of $G$, normalized so that $\langle\phi^{\mathrm{Gauss}}\circ\sigma,\phi^{\mathrm{Gauss}}\rangle=1$. It follows that

[TABLE]

Using $\tilde{K}_{\beta}\tilde{\phi}_{\beta}=\tilde{\phi}_{\beta}$ and $\tilde{G}^{*}(\phi^{\mathrm{Gauss}}\circ\sigma)=\phi^{\mathrm{Gauss}}\circ\sigma$, we get

[TABLE]

which goes to zero by Proposition 6.9 (see also Corollary 6.10). ∎

7. A Brascamp-Lieb type covariance estimate for $m=\infty$

Here we prove Proposition 2.10. Key to the proof is a matrix lower bound $A$ for the Hessian of $\mathcal{E}_{N}$. For Gaussian measures with probability density proportional to $\exp(-\frac{\beta}{2}\langle z,Az\rangle)$ and test functions $f_{i}=z_{i}$, $g_{j}=z_{j}$, we end up estimating the covariance $C_{ij}=([\beta A]^{-1})_{ij}$. We follow [Men14], see also [OR07].

Proof of Proposition 2.10.

Revisiting the proof of Lemma 3.3, we obtain bounds on matrix elements of the Hessian. Let $N\in\mathbb{N}$, $z\in[z_{\mathrm{min}},z_{\mathrm{max}}]^{N-1}$. For $1\leq i<j\leq N-1$ we have

[TABLE]

with

[TABLE]

For $1\leq i\leq N-1$ we also have

[TABLE]

by Assumption 1(iv). Moreover

[TABLE]

again by Assumption 1(iv). Let $A_{N}$ be the $(N-1)\times(N-1)$-matrix with diagonal $\rho$ and off-diagonal entries $-\kappa_{|j-i|}$; notice that $\eta,\kappa_{j-i},\rho$ do not depend on $N$. $A_{N}$ is symmetric and positive-definite.

The previous estimates together with [Men14, Remark 2.6] show that the energy $\mathcal{E}_{N}$ satisfies the assumptions of [Men14, Theorem 2.3 and Proposition 3.5]. It follows that for all smooth $f,g:\mathbb{R}_{+}\to\mathbb{R}$,

[TABLE]

Let $X_{1},X_{2},\ldots$ be i.i.d. random variables with law

[TABLE]

and $S_{n}=X_{1}+\cdots+X_{n}$. We may decompose $A_{N}$ as $\rho\mathrm{Id}$ plus an off-diagonal matrix, write a Neumann series for the inverse, and find that for $i<j$

[TABLE]

Clearly

[TABLE]

By (7.1), we have $\mathbb{P}(X_{r}=\ell)\leq C/|\ell|^{s}$ for some constant $C>0$. Following [Men14, Proposition 3.5] we may estimate, for each $m\in\mathbb{N}$,

[TABLE]

Similar estimates apply to other $r$. Combining with (7.3) we find

[TABLE]

It follows that

[TABLE]

Notice that the series is convergent. The bound is plugged into the estimate (7.2) and the proposition follows by passing to the limit $N\to\infty$. ∎

Acknowledgments

We thank Nils Berglund, Andrew Duncan, André Schlichting, and Martin Slowik for helpful discussions.