---
author:
- Cédric Bernardin
title: Superdiffusion of energy in Hamiltonian systems perturbed by a conservative noise
---

Introduction
============

Transport properties of one-dimensional Hamiltonian systems consisting of coupled oscillators on a lattice have been the subject of many theoretical and numerical studies, see the review papers [@BLR; @D; @LLP].  Despite many efforts, our knowledge of the fundamental mechanisms necessary and/or sufficient to have a normal diffusion remains very limited. 

Consider a $1$-dimensional chain of oscillators indexed by $x \in {{\mathbb Z}}$, whose formal Hamiltonian is given by $${{{\mathcal H}}} = \sum_{x \in {{\mathbb Z}}} \left[ \frac{p_x^2}{2} +V(r_x) \right],$$ where $r_x=q_{x+1}-q_x$ is the “deformation” of the lattice, $q_x$ being the displacement of the atom $x$ from its equilibrium position and $p_x$ its momentum.  The interaction potential $V$ is a smooth positive function growing at infinity fast enough.  The energy $e_x$ of atom $x \in {{\mathbb Z}}$ is defined by $$e_x = \frac{p_x^2}{2} + V(r_x).$$

Our goal is to understand the macroscopic energy diffusion properties for the corresponding Hamiltonian dynamics $$\cfrac{dr_x}{dt} = p_{x+1} -p_x, \quad \cfrac{dp_x}{dt} = V'(r_{x}) -V'(r_{x-1}), \quad x\in {{\mathbb Z}}.$$ Under suitable conditions on $V$, the infinite dynamics is well defined for a large class of initial conditions. 

Apart from the total energy $\sum_x e_x$, observe that the total momentum $ \sum_x p_x$ and the total deformation $\sum_x r_x$ of the lattice are formally conserved.  This is a consequence of the following microscopic continuity equations: $$\begin{aligned}
&&\frac{de_x}{dt} + \nabla [j_{x-1,x}^e]=0, \quad j_{x,x+1}^e= -p_{x+1} V' (r_x),\\
&&\frac{dp_x}{dt} + \nabla[ -V' (r_{x-1})]=0,\\
&&\frac{dr_x}{dt} + \nabla[ -p_x ]=0.\end{aligned}$$ The function $j_{x,x+1}^e$ is the current of energy going from $x$ to $x+1$. The main open problem ([@LP], [@Sz]) concerning the foundation of statistical mechanics based on classical mechanics is precisely to show that the three quantities above are the only quantities which are conserved by the dynamics.  In some sense, it means that the dynamics, evolving on the manifold defined by fixing the total energy, the total momentum and the total deformation, is ergodic.  Of course, the last sentence does not make sense since we are in infinite volume and $\sum_x e_x, \sum_x p_x$ and $\sum_x r_x$ are typically infinite.  Nevertheless, an alternative meaningful definition will be proposed and discussed in Section \[sec:ergodicity\]. 

Numerical simulations provide a strong evidence of the fact that one dimensional chains of anharmonic oscillators conserving momentum are [[^1]]{} superdiffusive.  It shall be noticed that there is no explanation of this, apart from heuristic considerations, and that some models which do not conserve momentum can also display anomalous diffusion of energy (see [@GDL]). 

An interesting area of current research consists in studying this problem for hybrid models where a stochastic perturbation is superposed to the deterministic evolution.  Even if the problem is considerably simplified, several open challenging questions can be addressed for these systems.  The first benefit of the introduction of stochasticity in the models is to guarantee the ergodicity that we are not able to show for purely deterministic systems.  The added noise must be carefully chosen in order not to destroy the conservation laws we are interested in.  In particular, the noise shall conserve energy.  But we will consider a noise conserving also some of the other quantities conserved by the underlying Hamiltonian dynamics, e.g. the momentum, the deformation or any linear combination of them. 

The paper is organized as follows.  In Section \[sec:ergodicity\] we discuss the problem of the ergodicity of the infinite dynamics mentioned above and the possible stochastic perturbations we can add to the deterministic dynamics to obtain ergodic dynamics.  In Section \[sec:harmonic\] we review some results obtained in the context of harmonic chains perturbed by a conservative noise and we discuss the case of anharmonic chains in the last section. 

Ergodicity {#sec:ergodicity}
==========

Let us first generalize the models introduced above ([@BS]).  Let $U$ and $V$ be smooth positive potentials growing at infinity fast enough and let ${{{\mathcal H}}}:={{{\mathcal H}}}_{U,V}$ be the Hamiltonian $${{{\mathcal H}}}_{U,V} =\sum_{x \in {{\mathbb Z}}} \left[ U(p_x) +V(r_x) \right].$$ The corresponding Hamiltonian dynamics satisfy $$\label{eq:Ham-gen}
\frac{dr_x}{dt}= U'(p_{x+1}) -U' (p_{x}), \quad \frac{dp_x}{dt} =V'(r_x) -V'(r_{x-1}), \quad x \in {{\mathbb Z}}.$$ The energy of particle $x$ is defined by $e_x = U(p_x) + V(r_x)$. The three formal quantities $\sum_x e_x$, $\sum_x r_x$ and $\sum_x p_x$ are conserved by the dynamics.  The fundamental question we address in this section is: are they the only ones?  In finite volume, i.e. replacing the lattice ${{\mathbb Z}}$ by a finite box $\Lambda$, this would correspond to the usual notion of ergodicity for Hamiltonian flows with a finite number of degrees of freedom.  But since we consider the dynamics in infinite volume the notion of conserved quantity has to be properly defined.  The way we follow to attack the problem is to detect the existence of a non-trivial conserved quantity through the existence of a non-trivial invariant state for the infinite dynamics. 

Let $\Omega=({{\mathbb R}}\times {{\mathbb R}})^{{{\mathbb Z}}}$ be the phase space of the dynamics and let us denote a typical configuration by $\omega=(r,p) \in \Omega$. For simplicity we assume that for any $(\beta, \lambda, \lambda' ) \in (0, +\infty) \times {{\mathbb R}}\times {{\mathbb R}}$, the partition function $$Z(\beta, \lambda, \lambda')= \int_{{{\mathbb R}}\times {{\mathbb R}}} e^{-\beta[U(a) +V(b)] -\lambda b - \lambda' a } da \,  db$$ is finite.  Let $\mu_{\beta, \lambda,\lambda'}$ be the product Gibbs measures on $\Omega$ defined by $$d\mu_{\beta, \lambda,\lambda'} (\omega) = \prod_{x \in {{\mathbb Z}}} \frac{1}{Z(\beta, \lambda,\lambda')} \exp \left[ - \beta [U(p_x) + V(r_x)] - \lambda r_x - \lambda' p_x \right] dr_x dp_x.$$ We assume that (\[eq:Ham-gen\]) is well defined for a subset $\Omega_{\beta, \lambda,\lambda'}$ of full measure with respect to $\mu_{\beta, \lambda,\lambda'}$, that the latter is invariant for (\[eq:Ham-gen\]), and that it is possible to define a strongly continuous semigroup in ${{{\mathbb L}}}^2 (\mu_{\beta,\lambda,\lambda'})$ with formal generator $${{{\mathcal A}}}_{U,V} = \sum_{x \in {{\mathbb Z}}} \left[(U'(p_{x+1}) -U'(p_x)) \partial_{r_x} + (V'(r_x) -V'(r_{x-1})) \partial_{p_x} \right].$$ All that can be proved under suitable assumptions on $U$ and $V$ ([@FFL], [@BObook]). 

In order to explain what is meant by ergodicity of the infinite volume dynamics we need to introduce some notation.  For any topological space $X$ equipped with its Borel $\sigma$-algebra we denote by ${{{\mathcal P}}} (X)$ the convex set of probability measures on $X$. The relative entropy $H(\nu|\mu)$ of $\nu \in {{{\mathcal P}}} (X)$ with respect to $\mu \in {{{\mathcal P}}} (X)$ is defined as $$H(\nu | \mu) = \sup_{\phi} \left\{ \int \phi \, d\nu - 
\log \left( \int e^{\phi} \, d\mu \right) \right\},$$ where the supremum is carried over all bounded measurable functions $\phi$ on $X$.

Let $\theta_x, x \in {{\mathbb Z}}$, be the shift by $x$: $(\theta_x \omega)_z=\omega_{x+z}$. For any function $g$ on $\Omega$, $\theta_x g$ is the function such that $(\theta_xg)(\omega) 
= g(\theta_x \omega)$. For any probability measure $\mu \in {{{\mathcal P}}} (\Omega)$, $\theta_x \mu \in {{{\mathcal P}}} (\Omega)$ is the probability measure such that, for any bounded function $g: \Omega \to {{\mathbb R}}$, it holds $\int_\Omega g \, d (\theta_x \mu)= \int_\Omega \theta_x g \, d\mu$. If $\theta_x \mu = \mu$ for any $x$ then $\mu$ is said to be translation invariant. 

If $\Lambda $ is a finite subset of ${{\mathbb Z}}$ the marginal of $\mu \in {{{\mathcal P}}} (\Omega)$ on ${{\mathbb R}}^{\Lambda}$ is denoted by $\mu |_{\Lambda}$. The relative entropy of $\nu \in {{{\mathcal P}}} (\Omega)$ with respect to $\mu \in {{{\mathcal P}}} (\Omega)$ in the box $\Lambda$ is defined by $H(\nu |_{\Lambda} \, | \, \mu |_{\Lambda} )$ and is denoted by $H_{\Lambda} (\nu| \mu)$. We say that a translation invariant probability measure $\nu \in {{{\mathcal P}}} (\Omega)$ has finite entropy density (with respect to $\mu$) if there exists a finite positive constant $C$ such that for any finite $\Lambda \subset {{\mathbb Z}}$, $H_{\Lambda} (\nu | \mu) \le C | \Lambda|$. In fact, if this condition is satisfied, then the limit $$\overline{H} (\nu|\mu)=\lim_{|\Lambda| \to \infty} \frac{H_{\Lambda} (\nu | \mu) }{|\Lambda|}$$ exists and is finite (see [@FFL]).  It is called the entropy density of $\nu$ with respect to $\mu$.

We are now in position to define ergodicity. 

\[def:ergo\] We say that the infinite volume dynamics with infinitesimal generator ${{{\mathcal A}}}_{U,V}$ is *ergodic* if the following claim is true: If $\nu \in {{{\mathcal P}}} (\Omega)$ is a probability measure invariant by translation, invariant by the dynamics generated by ${{{\mathcal A}}}_{U,V}$ and with finite entropy density with respect to $\mu_{1,0,0}$, then $\nu$ is a mixture of the $\mu_{\beta,\lambda,\lambda'}, \beta>0, \lambda, \lambda' \in {{\mathbb R}}$.

In the harmonic case ($U(z)=V(z)=z^2/2$) and for the Toda lattice ($U(z)=z^2/2$, $V(z) =e^{-z} +z -1$), the infinite system is completely integrable and an infinite number of conserved quantities can be explicitly written.  It follows that they are not ergodic in the sense above.  Nevertheless we expect that for a very large class of potentials, the Hamiltonian dynamics are ergodic and that these two cases are exceptional. 

In order that the infinite dynamics enjoy good ergodic properties, we superpose to the deterministic evolution a stochastic noise. 

Given a sequence $u= (u_y)_{y \in {{\mathbb Z}}} \in {{\mathbb R}}^{{{\mathbb Z}}}$ and a site $x\in {{\mathbb Z}}$, we denote by $u^x$ (resp.  $u^{x,x+1}$) the sequence defined by $(u^x)_y =u_y$ if $y \ne x$ and $(u^x)_x = -u_x$ (resp.  $(u^{x,x+1})_y = u_y$ if $y \ne x,x+1$, $(u^{x,x+1})_{x}= u_{x+1}$ and $(u^{x,x+1})_{x+1} = u_x$).  We consider the following noises (jump processes) whose generators are defined by their actions on functions $f: \Omega \to {{\mathbb R}}$ according to:

1.   $({{{\mathcal S}}}_{flip}^p f)(r,p)= \sum_{x} \left[ f(r,p^x) -f(r,p)\right]$.

2.   $({{{\mathcal S}}}_{flip}^r f)(r,p)= \sum_{x} \left[ f(r^x,p) -f(r,p)\right]$.

3.   $({{{\mathcal S}}}_{ex}^p f)(r,p)= \sum_{x} \left[ f(r,p^{x,x+1}) -f(r,p)\right]$.

4.   $({{{\mathcal S}}}_{ex}^r f)(r,p)= \sum_{x} \left[ f(r^{x,x+1},p) -f(r,p)\right]$.

If $U$ is even then the noise ${{{\mathcal S}}}_{flip}^p$ conserves the energy, the deformation but not the momentum; if $U$ is odd the noise has little interest for us since the energy conservation is destroyed.  Similarly, if $V$ is even the the noise ${{{\mathcal S}}}_{flip}^r$ conserves the energy and the momentum but not the deformation.  The noises ${{{\mathcal S}}}_{ex}^p$ and ${{{\mathcal S}}}_{ex}^r$ conserve the energy, the deformation and the momentum. 

Let now $\gamma>0$ and denote by ${{{\mathcal L}}}$ the generator of the infinite Hamiltonian dynamics generated by ${{{\mathcal A}}}_{U,V}$ perturbed by one of the previous noise ${{{\mathcal S}}}$ with intensity $\gamma$, i.e. ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}$.

The dynamics generated by ${{{\mathcal L}}}$ is ergodic in the sense that if $\nu \in {{{\mathcal P}}} (\Omega)$ is a probability measure invariant by translation, invariant by the dynamics generated by ${{{\mathcal L}}}$ and with finite entropy density with respect to $\mu_{1,0,0}$, then it holds:

1.   If $U$ even and ${{{\mathcal S}}}={{{\mathcal S}}}_{flip}^{p}$ then $\nu$ is a mixture of the $\mu_{\beta,\lambda,0}$;

2.   If $V$ is even and ${{{\mathcal S}}}={{{\mathcal S}}}_{flip}^{r}$ then $\nu$ is a mixture of the $\mu_{\beta,0,\lambda'}$.

3.   If ${{{\mathcal S}}}={{{\mathcal S}}}_{ex}^{p}$ or ${{{\mathcal S}}}={{{\mathcal S}}}_{ex}^{r}$ then $\nu$ is a mixture of the $\mu_{\beta,\lambda,\lambda'}$.

The main motivation to establish such a theorem is that by using Yau’s relative entropy method ([@Y]) in the spirit of Olla-Varadhan-Yau ([@OVY]), it is possible to show that if the infinite volume dynamics is ergodic then the propagation of local equilibrium holds in the hyperbolic time scale, before the appearance of the shocks.  As a consequence, the dynamics has a set of compressible Euler equations as hydrodynamic limits ([@BObook], [@BS]).  Observe that this is true also for the deterministic dynamics so that the rigorous derivation of the Euler equations from the first principles of the mechanics in the smooth regime is “reduced” to prove that the dynamics generated by ${{{\mathcal A}}}_{U,V}$ is ergodic. 

Harmonic chains  {#sec:harmonic}
================

Role of the conservation of momentum and deformation
----------------------------------------------------

We consider here the specific (harmonic) case $V(z)=U(z)=z^2/2$. The dynamics is then linear and can be solved analytically using Fourier transform.  Let us introduce a new macroscopic variable $\eta \in {{\mathbb R}}^{{{\mathbb Z}}}$ defined from $(p,r) \in \Omega$ by setting $$\label{eq:eta}
\eta_{2x}=r_x, \quad \eta_{2x+1}=p_{x+1}, \quad x\in {{\mathbb Z}}.$$ Then, the Hamiltonian dynamics can be rewritten in the form $$\label{eq:etadyn}
\frac{d\eta_x}{dt} = V'(\eta_{x+1}) -V' (\eta_{x-1}), \quad x \in {{\mathbb Z}}.$$

We introduce the $k$th mode ${\widehat \eta} (k, \cdot)$ for $k \in {{{\mathbb T}}} = 
\mathbb{R}/{{\mathbb Z}}$, the one-dimensional torus of length $1$: $${\widehat \eta} (t,k) =\sum_{x \in {{\mathbb Z}}} \eta_x (t) \, e^{2 {\mathrm{i}}\pi k x}.$$ Then, the equations of motion are equivalent in the sense of distributions to the following decoupled system of first order differential equations: $$\frac{d{\widehat \eta}}{dt} (t,k)  = {\mathrm{i}}\omega (k) \, {\widehat \eta} (t,k),$$ where the dispersion relation $\omega (k)$ reads $$\omega (k) =- 2 \sin (2\pi k),$$ and the group velocity $v_{\rm g}$ is $$v_{\rm g} (k)  = \omega' (k) = -4\pi \cos (2\pi k).$$ By inverting the Fourier transform, the solution can be written as $$\eta_{x} (t) = \int_{{{\mathbb T}}} {\widehat \eta} (t,k) \, e^{- 2 {\mathrm{i}}\pi k x} \, dk.$$ If the initial configuration $\eta (0)$ is in $\ell_2$ the well defined energy of the $k$th mode $$E_k (t)= \frac{1}{4\pi} |{\widehat \eta} (t,k)|^2 = E_k (0)$$ is conserved by the time evolution, and the total energy current ${\tilde J}^e= \sum_{x \in {{\mathbb Z}}} j^{e}_{x,x+1}$ takes the simple form $${\tilde J}^e = \int_{{{\mathbb T}}} v_{\rm g}(k) E_k \, dk.$$

We interpret the waves ${\widehat \eta} (k,t)$ as fictitious particles (phonons in solid state physics).  In the absence of nonlinearities, they travel the chain without scattering.  The diffusion of energy is then said to be ballistic.  If the potential is non-quadratic, it may be expected that the nonlinearities produce a scattering responsible for the diffusion of the energy.  Nevertheless, the conservation of the deformation and of the momentum implies that $\sum_x (r_x +p_x)$ is conserved $$\label{eq:consj}
{\widehat \eta} (t, 0) = {\widehat \eta} (0,0).$$ The identity (\[eq:consj\]) is valid even if $U \ne V$ and $U,V$ are not quadratic.  It means that the $0$th mode is not scattered at all and crosses the chain ballistically.  In fact, the modes with small wave number $k$ do not experience a strong scattering and they therefore contribute to the observed anomalous diffusion of energy. 

It is usually explained that momentum conservation plays a major role in the anomalous diffusion of energy but it is clear that the deformation conservation plays exactly the same role as momentum and that it is the conservation of their sum which is the real ingredient producing anomalous diffusion of energy (see Theorem \[th:harmflip\] and Theorem \[th:anharmflip\]). 

Green-Kubo formula
------------------

The signature of an anomalous diffusion of energy can be seen at the level of the Green-Kubo formula.  When transport of energy is normal, meaning that the macroscopic equations such as the Fourier’s law or heat equation hold, the transport coefficient appearing in these equations can be expressed by the famous Green-Kubo formula.  In order to define the latter we need to introduce some notations.  Since the discussion about the Green-Kubo formula is not restricted to the harmonic case we go back to a generic anharmonic model in the rest of the Subsection. 

Recall that the probability measures $\mu_{\beta, \lambda, \lambda'}$ form a family of invariant probability measures for the infinite dynamics generated by ${{{\mathcal A}}}_{U,V}$. The following thermodynamic relations (which are valid since we assumed that the partition function $Z$ is well defined on $(0,+\infty) \times {\mathbb R} \times {\mathbb R}$) relate the chemical potentials $\beta, \lambda, \lambda'$ to the mean energy $e$, the mean deformation $u$, the mean momentum $\pi$ under $\mu_{\beta,\lambda,\lambda'}$: $$\begin{aligned}
\label{eq:tr}
&&e(\beta,\lambda, \lambda')= \mu_{\beta,\lambda,\lambda'} (U(p_x) +V(r_x))= -\partial_{\beta}\Big(\log Z(\beta,\lambda, \lambda')\Big),\\
&&u(\beta,\lambda,\lambda') =\mu_{\beta,\lambda,\lambda'} (r_x)= -\partial_{\lambda} \Big(\log Z(\beta,\lambda, \lambda')\Big), \\
&&\pi(\beta,\lambda,\lambda') =\mu_{\beta,\lambda,\lambda'} (p_x)= -\partial_{\lambda'} \Big(\log Z(\beta,\lambda,\lambda')\Big).\end{aligned}$$ These relations can be inverted by a Legendre transform to express $\beta$, $\lambda$ and $\lambda'$ as a function of $e$, $u$ and $\pi$. Define the thermodynamic entropy $S:(0,+\infty) \times {{\mathbb R}}\times {{\mathbb R}}\to [0,+\infty]$ as $$S(e,u,\pi)= \inf_{ \lambda, \lambda' \in {{\mathbb R}}^2, \beta >0} \Big\{ \beta e + \lambda u + \lambda' \pi 
+ \log Z (\beta,\lambda, \lambda') \Big\}.$$ Let ${{{\mathcal U}}}$ be the convex domain of $(0,+\infty) \times {{\mathbb R}}\times {{\mathbb R}}$ where $S(e,u,\pi) <+ \infty$ and $\mathring{{{\mathcal U}}}$ its interior.  Then, for any $(e,u,\pi):=(e(\beta, \lambda,\lambda'),u(\beta,\lambda, \lambda'), \pi(\beta,\lambda,\lambda')) 
\in {\mathring{{{\mathcal U}}}}$, the parameters $\beta,\lambda,\lambda'$ can be obtained as $$\label{eq:14}
\beta = (\partial_e S ) (e,u,\pi), 
\qquad 
\lambda= (\partial_u S) (e,u,\pi),
\qquad
\lambda'= (\partial_\pi S)(e,u,\pi)$$

These thermodynamic relations allow us to parameterize the Gibbs states by the average values of the conserved quantities $(e,u,\pi)$ rather than by the chemical potentials $(\beta,\lambda,\lambda')$. Thus, we denote by $\nu_{e,u,\pi}$ the Gibbs measure $\mu_{\beta,\lambda,\lambda'}$ where $(e,u,\pi)$ are related to $(\beta,\lambda,\lambda')$ by (\[eq:14\]).  Let $J^e:=J^e (e,u,\pi)= \nu_{e,u,\pi} (j^e_{x,x+1})$ be the average of the energy current $j^{e}_{x,x+1}= -U' (p_x) V'(r_x)$ and define the normalized energy current ${\hat j}^e_{x,x+1}$ by $${\hat j}^e_{x,x+1} = j^e_{x,x+1} - J^e  - (\partial_e J^e) (e_x -e) -(\partial_u J^e) (r_x -u) - (\partial_\pi J^e) (p_x -\pi).$$ The normalized energy current is the part of the centered energy current which is orthogonal in ${{{\mathbb L}}}^2 (\nu_{e,u,\pi})$ to the space spanned by the conserved quantities. 

Up to multiplicative thermodynamic parameters (see [@Sp] for details) that we neglect to simplify the notations, the Green-Kubo formula [[^2]]{} is nothing but $$\kappa (e,u,\pi) := \int_0^{\infty} \; \sum_{x \in {{\mathbb Z}}} {\mathbb E}_{\nu_{e,u,\pi}} \left[ {\hat j}^e_{x,x+1} (\omega(t)) \; {\hat j}_{0,1}^e (\omega(0)) \right] \; dt$$ where ${\mathbb E}_{\nu_{e,u,\pi}}$ denotes the expectation corresponding to the law of the infinite volume dynamics $(\omega(t))_{t\ge 0}$ generated by ${{{\mathcal A}}}_{U,V}$ with initial condition $\omega (0)$ distributed according to the equilibrium Gibbs measure $\nu_{e,u,\pi}$. The definition of $\kappa (e,u,\pi)$ is formal but the way we adopt to give it a mathematically well posed definition is to introduce a small parameter $z>0$ and define $\kappa (e,u,\pi)$ as $$\label{eq:GK}
\kappa (e,u,\pi) = \limsup_{z \to 0} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal A}}}_{U,V})^{-1} {\hat j}^e_{0,1} \gg_{e,u,\pi}$$ where the inner-product $\ll \cdot, \cdot \gg_{e,u,\pi}$ is defined for local square integrable functions $f,g : \Omega \to {{\mathbb R}}$ by $$\ll f \, , \, g \gg_{e,u,\pi} \; = \; \sum_{x \in {{\mathbb Z}}} \left[ \left( \int f \theta_x g \, d\nu_{e,u, \pi} \right) - \left(\int f d\nu_{e,u,\pi} \right) \left(\int g d\nu_{e,u,\pi} \right) \right].$$ Since $(z-{{{\mathcal A}}}_{U,V})^{-1} {\hat j}^e_{0,1}$ is not a local function, the term on the RHS of (\[eq:GK\]) has to be interpreted in the Hilbert space obtained by the completion of the space of local bounded functions with respect to the inner product $\ll \cdot, \cdot \gg_{e,u,\pi}$.

The superdiffusion (resp.  normal diffusion) of energy corresponds to an infinite (resp.  finite) value for $\kappa (e,u,\pi)$. In order to study the superdiffusion, it is of interest to estimate the time decay of the autocorrelation of the normalized current $$C (t):= C_{e,u,\pi} (t) = \sum_{x \in {{\mathbb Z}}}  {\mathbb E}_{\nu_{e,u,\pi}} \left[ {\hat j}^e_{x,x+1} (\omega(t)) \; {\hat j}_{0,1}^e (\omega(0)) \right].$$ It is in general easier to estimate the behavior of the Laplace transform $L(z)=\int_{0}^{\infty} e^{-z t} C(t) dt$ as $z\to 0$. Roughly, if $L(z) \sim z^{- \delta}$ for some $\delta \ge 0$ then $C(t) \sim t^{\delta -1}$ as $t \to + \infty$. Observe also that $$L(z)= \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal A}}}_{U,V})^{-1} {\hat j}_{0,1} \gg_{e,u,\pi}.$$

Harmonic chain perturbed by a conservative stochastic noise
-----------------------------------------------------------

We consider now the particular case $U(z)=V(z)=z^2 /2$ and study the Green-Kubo formula for the perturbed dynamics generated by ${{{\mathcal L}}} ={{{\mathcal A}}}_{U,V} +\gamma {{{\mathcal S}}}$ where ${{{\mathcal S}}}$ is one of the noises introduced in Section \[sec:ergodicity\].  Since, depending of the form of the noise, the momentum conservation law (resp.  deformation conservation law) can be suppressed, the corresponding Green-Kubo formula shall be modified by setting $\pi=0$ and $\partial_\pi J^e=0$ (resp.  $u=0$ and $\partial_u J^e =0$). 

We have the following theorem which shows that if momentum conservation law or deformation conservation law is destroyed by the noise then a normal behavior occurs. 

\[th:harmflip\] Let $U$ and $V$ be quadratic potentials. 

1.   Consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{flip}^p$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,u, 0}$$ exists, is finite and strictly positive and can be explicitly computed. 

2.   Consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{flip}^r$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,0,\pi}$$ exists, is finite and strictly positive and can be explicitly computed. 

It shall be noticed that the second statement is a direct consequence of the first one since the process of the second item is equal to the first one by the transformation $$r_x \to p_x,\quad p_{x} \to r_{x-1}.$$ However, the interest of the second statement is to show that [*even if momentum is conserved*]{}, a normal diffusion of energy occurs.  This is because the deformation is no longer conserved. 

The following theorem shows that if the noise added conserves momentum [*and*]{} deformation then the situation is very different since an anomalous diffusion of energy is observed. 

Let $U$ and $V$ be quadratic potentials. 

1.   Consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{ex}^p$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} z^{1/2} \; \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,u, \pi}$$ exists, is finite and strictly positive and can be explicitly computed. 

2.   Consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{ex}^r$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} z^{1/2} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,u,\pi}$$ exists, is finite and strictly positive and can be explicitly computed. 

In particular, in each of the previous case the Green-Kubo formula yields an infinite conductivity. 

 Anharmonic chains
==================

We consider now the anharmonic case.  For deterministic chains generated by ${{{\mathcal A}}}_{U,V}$ we expect usually a superdiffusive behavior of the energy.  If a noise ${{{\mathcal S}}}$ is superposed to the dynamics, we expect that transport is normal for ${{{\mathcal S}}}={{{\mathcal S}}}_{flip}^p$ and ${{{\mathcal S}}}={{{\mathcal S}}}_{flip}^r$ and superdiffusive if ${{{\mathcal S}}}={{{\mathcal S}}}_{ex}^p$ or ${{{\mathcal S}}}={{{\mathcal S}}}_{ex}^r$.

The following theorem generalizes Theorem \[th:harmflip\] to the anharmonic case showing that a noise destroying momentum conservation law or deformation conservation law produces normal transport.  This shows that, also in the anharmonic case, momentum conservation alone is not responsible of anomalous diffusion of energy but that deformation conservation law plays a similar role. 

\[th:anharmflip\] Let $U$ and $V$ be smooth potentials such that there exists a constant $c>0$ such that $$c \le U '' \le c^{-1}, \quad c \le V'' \le c^{-1}.$$

1.   Assume $U$ even and consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{flip}^p$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,u,0}$$ exists and is finite. 

2.   Assume $V$ even and consider the system generated by ${{{\mathcal L}}}={{{\mathcal A}}}_{U,V} + \gamma {{{\mathcal S}}}_{flip}^r$, $\gamma>0$. Then the following limit $$\lim_{z \to 0} \ll {\hat j}^e_{0,1}\, , \, (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,0,\pi}$$ exists and is finite. 

The second statement is a direct consequence of the first one by the symmetry argument evoked for Theorem \[th:harmflip\].  The upper bounds on $U''$ and $V''$ are here to assure the existence of the infinite volume dynamics. 

For simplicity assume that $u=0$ and $\beta:=\beta(e,u,0)=1$. The first statement has been proved in [@BO] in the particular case $U(z)=z^2/2$. The generalization to a non quadratic smooth even potential $U$ is straightforward.  In [@BO], since $U(z)=z^2/2$, we used Hermite polynomials which are orthogonal w.r.t.  the Gaussian measure $d\mu (z) = (2\pi)^{-1/2} \exp\{ -z^2 /2\} dz$. In the present case, the only difference is that we have to replace the Hermite basis by any orthogonal polynomial basis $\{P_n \}_{n \ge 0}$ with respect to the probability measure ${{{\mathcal N}}}^{-1} \exp (- U(z) ) dz$ (with ${{{\mathcal N}}}$ a normalization constant) which satisfies $P_n$ odd if $n$ odd and even otherwise.  Then the proof is exactly the same. 

It would be now of interest to show that if we perturb the dynamics generated by ${{{\mathcal A}}}_{U,V}$ by ${{{\mathcal S}}}_{ex}^p$ or by ${{{\mathcal S}}}_{ex}^r$ then anomalous diffusion of energy occurs [[^3]]{}.  This is an open question and as far as we know the only result going in this direction has been obtained in [@BG]. 

The model considered in [@BG] is the dynamics generated by ${{{\mathcal A}}}_{U,V}$ with $U=V$ taking the particular form $V(z)=e^{-z} +z -1$, perturbed by a noise ${{{\mathcal S}}}$ which conserves energy and $\sum_{x \in {{\mathbb Z}}} (r_x + p_x)$. More exactly, let us rewrite the Hamiltonian dynamics (\[eq:Ham-gen\]) by using the variable $\eta:=(\eta_x)_{x \in {{\mathbb Z}}} \in {{\mathbb R}}^{{{\mathbb Z}}}$ defined by (\[eq:eta\]).  Then we get the equations of motion given by (\[eq:etadyn\]).  With these new variables, the total energy is $2 \sum_{x} V(\eta_x)$, the total deformation is $\sum_{x} \eta_{2x}$ and the total momentum is $\sum_{x} \eta_{2x+1}$. The noise ${{{\mathcal S}}}$ superposed to the dynamics acts on local functions $f: {{\mathbb R}}^{{{\mathbb Z}}} \to {{\mathbb R}}$ according to $$({{{\mathcal S}}} f) (\eta) = \sum_{x \in {{\mathbb Z}}} \left[ f(\eta^{x,x+1}) -f(\eta) \right].$$

Observe that the noise conserves the energy, destroys the momentum and the deformation conservation laws but conserves $\sum_{x} \eta_x = \sum_{x} (p_x +r_x)$, which as explained above is the quantity (that we call the “volume” to follow the terminology used in [@BG]) responsible of the anomalous diffusion of energy.  Since we have now only two conserved quantities (the energy and the volume), the Gibbs states of the perturbed dynamics are given by $\{ \mu_{\beta, \lambda, \lambda} \}_{\beta>0, \lambda}$ or equivalently by $\{ \nu_{e, \pi, \pi} \, ; \, e>0, \pi \}$. The normalized energy current is given by $${\hat j}^e_{x,x+1} (\eta)= -2 V' (\eta_x) V' (\eta_{x+1}) +2 \tau^2 +2 \partial_e (\tau^2) \, (2V(\eta_x) -e) + 2 \partial_\pi (\tau^2) \, (\eta_x -\pi)$$ with $\tau:=\tau(e,\pi) = \int V'(\eta_x) d\nu_{e,\pi,\pi}$.

Let $(e,\pi) \in (0,+\infty) \times {{\mathbb R}}$ such that $\nu_{e,\pi,\pi}$ is well defined.  Consider the dynamics with generator ${{{\mathcal L}}}={{{\mathcal A}}}_{exp} +\gamma {{{\mathcal S}}}$, $\gamma>0$, where $${{{\mathcal A}}}_{exp} = \sum_x (V'(\eta_{x+1}) -V' (\eta_{x-1})) \partial_{\eta_x},$$ and $V(z)=e^{-z}+z-1$. Then there exists a constant $c>0$ such that for any $z>0$ $$c z^{-1/4} \; \le \; \ll {\hat j}^e_{0,1} , (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,\pi,\pi} \le c^{-1} z^{-1/2}.$$ It follows that the Green-Kubo formula of the energy transport coefficient yields an infinite value. 

We expect that the system above belongs to the KPZ class so that $ \ll {\hat j}^e_{0,1} , (z-{{{\mathcal L}}})^{-1} {\hat j}^e_{0,1} \gg_{e,\pi,\pi}$ should diverge like $z^{-1/3}$. In the present state of the art no robust technique is available to show such result apart from the non-rigorous (but powerful) mode-coupling theory ([@MS], [@Sp2], [@VB]).  A second open problem is to generalize the previous theorem to other interaction potentials $V$. Numerical simulations have been reported in [@BS]. 

[^1]: See however the coupled-rotor model which displays normal behavior (see [@LLP], Section 6.4).  This is probably due to the fact that the position space is compact. 

[^2]: The transport coefficient is in fact a matrix whose size is the number of conserved quantities.  Since we are interested in the energy diffusion, we only consider the entry corresponding to the energy-energy flux. 

[^3]: However, if $U$ or $V$ is bounded, like for the rotors model, we expect that diffusion is normal.