Inverse Scattering Method Solves the Problem of Full Statistics of Nonstationary Heat Transfer in the Kipnis-Marchioro-Presutti Model

Basic Information

• Paper ID: 2112.02474
• Title: Inverse Scattering Method Solves the Problem of Full Statistics of Nonstationary Heat Transfer in the Kipnis-Marchioro-Presutti Model
• Authors: Eldad Bettelheim (Hebrew University), Naftali R. Smith (ENS Paris & Ben-Gurion University), Baruch Meerson (Hebrew University)
• Classification: cond-mat.stat-mech, math.PR, nlin.SI
• Submission Date: December 2021
• Paper Link: https://arxiv.org/abs/2112.02474

Abstract

This paper determines the complete statistical properties of nonstationary heat transfer over long times in the Kipnis-Marchioro-Presutti (KMP) lattice gas model by discovering and exploiting the complete integrability of the fundamental equations of macroscopic fluctuation theory (MFT). These equations are intimately related to the derivative nonlinear Schrödinger equation (DNLS), and the authors employ the Zakharov-Shabat inverse scattering method (ISM) adapted by Kaup and Newell (1978) for DNLS to solve them. For an initially localized thermal pulse, explicit results for the exact large deviation function of heat transport are obtained, and a nontrivial symmetry relation is discovered.

Research Background and Motivation

Core Problem

The core problem addressed in this research is determining the complete statistical properties of matter or energy flow in nonequilibrium macroscopic systems, particularly the nonstationary fluctuation statistics problem. This is a fundamental yet challenging problem in statistical physics.

Problem Significance

1. Theoretical Importance: A fundamental problem in nonequilibrium statistical physics, involving large deviation theory and fluctuation theorems
2. Practical Value: Heat conduction and diffusion processes have widespread applications in materials science, biophysics, and other fields
3. Mathematical Challenge: Exact solution of nonstationary processes is extremely difficult; existing methods are primarily limited to asymptotic analysis

Limitations of Existing Methods

• Steady-State Results: Significant progress has been made for nonequilibrium steady states, but nonstationary fluctuation analysis remains difficult
• Asymptotic Analysis: Previous work primarily obtained asymptotic behavior for small and large deviations, lacking exact solutions
• Numerical Constraints: Monte Carlo simulations have prohibitively high computational costs in the large deviation regime

Research Motivation

The KMP model, as a paradigmatic model for heat conduction, is the first model rigorously proven to satisfy Fourier's heat diffusion law at the coarse-graining level. This work aims to obtain the first exact analytical results for the statistics of nonstationary heat transport in this model.

Core Contributions

1. Discovery of Complete Integrability: First reveals the relationship between the MFT equations of the KMP model and the derivative nonlinear Schrödinger equation (DNLS)
2. Exact Analytical Solution: Obtains the first exact nonstationary large deviation results for the statistics of particle flow in interacting lattice gases under quenched initial conditions using the inverse scattering method
3. Nontrivial Symmetry: Discovers and proves a time-reversal mirror symmetry relation: $v(x,t) = -λu(-x,1-t)$
4. Complete Rate Function: Obtains explicit parametric representation of the large deviation function for heat transport
5. Asymptotic Behavior: Derives exact asymptotic expressions for small deviations ($j \to 0$) and large deviations ($|j| \to 1/2$)

Detailed Methodology

Task Definition

Input: Initial condition as a localized thermal pulse $u(x,t=0) = W\delta(x)$Output: Probability distribution of heat excess $J = \int_0^∞ u(x,T)dx - W/2$ observed on the right half-axis $x > 0$ at time $T$Objective: Solve the large deviation rate function $s(j)$, where $\ln P(J,T,W) \simeq -\sqrt{T}s(J/W)$

Theoretical Framework

Macroscopic Fluctuation Theory (MFT)

Starting from the stochastic Langevin equation: $\partial_t u = \partial_x(\partial_x u + \sqrt{2u}\eta)$

Through saddle-point approximation of the path integral, the MFT equations are obtained: $\partial_t u = \partial_x(\partial_x u + 2u^2 v)$$\partial_t v = \partial_x(-\partial_x v + 2uv^2)$

Boundary Conditions

• Initial Condition: $u(x,0) = \delta(x)$
• Terminal Condition: $v(x,1) = -λ\delta(x)$
• Constraint: $\int_0^∞ u(x,1)dx - 1/2 = j$

Core Technical Innovations

1. Integrability Identification

The key discovery is that the MFT equations are equivalent to the DNLS equation in imaginary time and imaginary space. A Lax pair is constructed: $\partial_x \psi = U(x,t,k)\psi, \quad \partial_t \psi = V(x,t,k)\psi$

where the compatibility condition $\partial_t U - \partial_x V + [U,V] = 0$ for matrices $U$ and $V$ is equivalent to the MFT equations.

2. Inverse Scattering Transform

Define the full-space propagator: $G(t,k) = \lim_{\substack{x \to \infty \\ y \to -\infty}} \begin{pmatrix} e^{ikx/2} & 0 \\ 0 & e^{-ikx/2} \end{pmatrix} T(x,y,t,k) \begin{pmatrix} e^{-iky/2} & 0 \\ 0 & e^{iky/2} \end{pmatrix}$

The time evolution of the scattering amplitude has a simple form: $G(t,k) = \begin{pmatrix} a(0,k) & \tilde{b}(0,k)e^{k^2 t} \\ b(0,k)e^{-k^2 t} & \tilde{a}(0,k) \end{pmatrix}$

3. Boundary Value Problem Solution

By computing $G(t,k)$ at $t=0$ and $t=1$, the key equation is obtained: $ik[Q_+(k) + Q_-(k)] - ikQ_-(k) \times ikQ_+(k) = -λike^{-k^2}$

where $Q_±(k)$ are the Fourier transforms of $v(z,0)$ in the regions $z \gtrless 0$.

4. Exact Solution

The solution is: $ikQ_±(k) = 1 - (1 \pm v_±)e^{\Phi_±(k)}$

where: $\Phi_±(k) = ± \int_{-\infty}^∞ \frac{\ln(1 + iλk'e^{-k'^2})}{k' - k \mp i0^+} \frac{dk'}{2πi}$

Main Results

Parametric Representation of the Rate Function

The relationship between heat excess and Lagrange multiplier: $j(λ) = \frac{1}{4πλ} \int_{-\infty}^∞ \frac{\ln(1 + λ^2k^2e^{-2k^2})}{k^2} dk$

The rate function: $s(λ) = Q_+(0) + \int_{-\infty}^∞ \frac{\text{Li}_2(-λ^2k^2e^{-2k^2})}{8πk^2} dk + \frac{λ}{2}$

where $\text{Li}_2(z)$ is the dilogarithm function.

Asymptotic Behavior

Small Deviation Limit ($j \to 0$)

$s(j \to 0) \simeq \sqrt{8π}j^2$

This yields Gaussian fluctuations with variance scaling as $T^{-1/2}$.

Large Deviation Limit ($|j| \to 1/2$)

$s(|j| \to 1/2) \simeq \frac{4}{3π} \ln^{3/2}\left(\frac{2}{π∆}\sqrt{\ln\frac{2}{π∆}\sqrt{\ln\frac{2}{π∆}\cdots}}\right)$

where $∆ = 1/2 - |j| \ll 1$, exhibiting singular behavior with nested logarithms.

Symmetry Relation

A nontrivial time-reversal symmetry is discovered: $v(x,t) = -λu(-x,1-t)$

This relation mathematically simplifies the solution of the boundary value problem.

Verification and Numerical Results

Monte Carlo Simulation

• Conducted $10^6$ direct Monte Carlo simulations to verify theoretical results
• Parameter settings: $T = 10^2$, number of lattice sites $2L+1 = 51$
• Theory predictions agree highly with numerical results

Numerical Algorithm Verification

The Chernykh-Stepanov iterative algorithm is used to numerically solve the MFT equations, verifying the correctness of the analytical results.

Related Work

Historical Development

1. KMP Model: Proposed by Kipnis, Marchioro, Presutti (1982), the first model rigorously proven to satisfy Fourier's law
2. Macroscopic Fluctuation Theory: Theoretical framework developed by Bertini et al. based on weak noise theory
3. Inverse Scattering Method: Zakharov-Shabat method, Kaup-Newell adaptation for DNLS

Related Research

• Steady-State Statistics: Important work by Derrida et al. on steady-state large deviations
• Step Initial Conditions: Previous research primarily focused on step function initial conditions
• KPZ Equation: Recent breakthroughs in interface growth statistics

Advantages of This Work

1. First Exact Solution: First complete analytical result for quenched initial conditions
2. Elegant Method: More direct compared to Fredholm determinant methods
3. Universality: Method generalizable to other integrable systems

Conclusions and Discussion

Main Conclusions

1. Complete Solution: First complete solution of large deviation statistics for nonstationary heat transport in the KMP model
2. Integrable Structure: Reveals the deep integrable structure underlying macroscopic fluctuation theory
3. Exact Formulas: Obtains explicit parametric representation of the rate function and asymptotic behavior

Limitations

1. Specific Initial Conditions: Results apply only to δ-function initial conditions
2. Long-Time Limit: Theory is valid for $T \gg 1$
3. One-Dimensional Systems: Results limited to one-dimensional lattice systems

Future Directions

1. Extension to SSEP: Similar analysis for the simple symmetric exclusion process
2. General Initial Conditions: Generalization to more general initial conditions
3. Higher Dimensions: Study of integrability in multidimensional systems
4. Experimental Verification: Verification in actual physical systems

In-Depth Evaluation

Strengths

1. Theoretical Breakthrough: Significant theoretical progress in nonequilibrium statistical physics
2. Mathematical Rigor: Application of inverse scattering method is rigorous and innovative
3. Complete Results: Complete description from small to large deviations
4. Comprehensive Verification: Triple verification through theory, numerics, and simulation

Weaknesses

1. Limited Scope: Restricted to specific models and initial conditions
2. Physical Intuition: Physical interpretation of certain mathematical results could be deeper
3. Computational Complexity: Although exact, calculations remain relatively complex

Impact

1. Academic Value: Provides new paradigm for exact solutions in nonequilibrium statistical physics
2. Methodological: Novel application of inverse scattering method in statistical physics
3. Inspirational: Provides insights for large deviation analysis of other integrable systems

Applicable Scenarios

• Exact analysis of one-dimensional heat conduction systems
• Large deviation theory of integrable stochastic processes
• Theoretical research in nonequilibrium statistical physics
• Applications of integrable systems in mathematical physics

References

This paper cites important literature in statistical physics, integrable systems, and large deviation theory, including:

• Work by Derrida et al. on macroscopic fluctuation theory
• Kaup-Newell inverse scattering method for DNLS
• Bertini et al.'s nonequilibrium statistical physics theory
• Recent work on large deviation analysis of the KPZ equation

Summary: This is a paper of significant theoretical importance in nonequilibrium statistical physics. It is the first to use integrable systems theory to exactly solve the nonstationary large deviation problem of a classical statistical model, providing a new research paradigm and exact results for the field.