One-dimensional superdiffusive heat propagation induced by optical phonon-phonon interactions

Universality classes of thermalization and energy diffusion

In this paper, we show that classical lattices can be classified into two universality classes for thermalization, based solely on the properties of their eigenmodes. This discovery is a consequence of our systematic multiwave quasiresonance analysis, a tool developed to this end. Lattices with extended modes belong to one class that can thermalize within a finite time, even when the nonlinearity strength is very weak, provided the system size is sufficiently large. In contrast, lattices with purely localized modes fall into another class. For these systems, the scaling behavior of thermalization time shifts stepwise from low-order to progressively higher-order quasiresonances as nonlinear strength decreases, implying that thermalization may become unattainable within a reasonable time for sufficiently weak nonlinearity strength. Furthermore, we show that the real-space energy diffusion behavior of the two classes is qualitatively different as well.

Nonequilibrium thermal rectification at the junction of harmonic chains

A thermal diode or rectifier is a system that transmits heat or energy in one direction better than in the opposite direction. We investigate the influence of the distribution of energy among wave numbers on the diode effect for the junction of two dissimilar harmonic chains. An analytical expression for the diode coefficient, characterizing the difference between heat fluxes through the junction in two directions, is derived. It is shown that the diode coefficient depends on the distribution of energy among wave numbers. For an equilibrium energy distribution, the diode effect is absent, while for non-equilibrium energy distributions the diode effect is observed even though the system is harmonic. We show that the diode effect can be maximized by varying the energy distribution and relative position of spectra of the two harmonic chains. Conditions are formulated under which the system acts as an ideal thermal rectifier, i.e., transmits heat only in one direction. The results obtained are important for understanding the heat transfer in heterogeneous low-dimensional nanomaterials.

Equilibration of sinusoidal modulation of temperature in linear and nonlinear chains

The equilibration of sinusoidally modulated distribution of the kinetic temperature is analyzed in the β-Fermi-Pasta-Ulam-Tsingou chain with different degrees of nonlinearity and for different wavelengths of temperature modulation. Two different types of initial conditions are used to show that either one gives the same result as the number of realizations increases and that the initial conditions that are closer to the state of thermal equilibrium give faster convergence. The kinetics of temperature equilibration is monitored and compared to the analytical solution available for the linear chain in the continuum limit. The transition from ballistic to diffusive thermal conductivity with an increase in the degree of anharmonicity is shown. In the ballistic case, the energy equilibration has an oscillatory character with an amplitude decreasing in time, and in the diffusive case, it is monotonous in time. For smaller wavelength of temperature modulation, the oscillatory character of temperature equilibration remains for a larger degree of anharmonicity. For a given wavelength of temperature modulation, there is such a value of the anharmonicity parameter at which the temperature equilibration occurs most rapidly.

Anomalous energy diffusion in two-dimensional nonlinear lattices

Heat transport in one-dimensional (1D) momentum-conserved lattices is generally assumed to be anomalous, thus yielding a power-law divergence of thermal conductivity with system length. However, whether heat transport in a two-dimensional (2D) system is anomalous or not is still up for debate because of the difficulties involved in experimental measurements or due to the insufficiently large simulation cell size. Here we simulate energy and momentum diffusion in the 2D nonlinear lattices using the method of fluctuation correlation functions. Our simulations confirm that energy diffusion in the 2D momentum-conserved lattices is anomalous and can be well described by the Lévy-stable distribution. As is expected, we verify that 2D nonlinear lattices with on-site potentials exhibit normal energy diffusion, independent of the dimension. Contrary to the hypothesis of a 1D system, we further clarify that anomalous heat transport in the 2D momentum-conserved system cannot be corroborated by the momentum superdiffusion any longer. Our findings offer some valuable insights into mechanisms of thermal transport in 2D system.