Crossover from ballistic to normal heat transport in the ϕ^{4} lattice: If nonconservation of momentum is the reason, what is the mechanism?

Antipersistent energy-current correlations in strong long-ranged Fermi-Pasta-Ulam-Tsingou-type models

We study heat transfer in one-dimensional Fermi-Pasta-Ulam-Tsingou-type systems with long-range (LR) interactions. The strength of the LR interaction between two lattice sites decays as a power σ of the inverse of their distance. We focus on the strong LR regime (0≤σ≤1) and show that the thermal transport behaviors are remarkably nuanced. Specifically, we observe that the antipersistent (negative) energy current correlation in this regime is intricately dependent on σ, displaying a nonmonotonic variation. Notably, a significant qualitative change occurs at σ_{c}=0.5, where with respect to other σ values the correlation shows a minimum negative value. Furthermore, our findings also demonstrate that within the long-time range considered, these antipersistent correlations will eventually vanish for certain σ>0.5. The underlying mechanisms behind these intriguing phenomena are related to the crossover of two diverse space-time scaling properties of equilibrium heat correlations and the various scattering processes of phonons and discrete breathers.

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.

Subdiffusive energy transport and antipersistent correlations due to the scattering of phonons and discrete breathers

While there are many physical processes showing subdiffusion and some useful particle models for understanding the underlying mechanisms have been established, a systematic study of subdiffusive energy transport is still lacking. Here we present convincing evidence that, in the range of system size investigated, the energy subdiffusion can take place in a Hamiltonian lattice system with both harmonic nearest-neighbor and anharmonic long-range interactions. In particular, we show that the interaction range dependence of antipersistent energy-current correlations are relevant to this special type of energy subdiffusion. The underlying mechanisms are related to the various scattering processes of phonons and discrete breathers. Our results shed light on understanding the extremely slow energy transport.

Modulational Instability of Delocalized Modes in fcc Copper

Delocalized nonlinear vibrational modes (DNVMs) are exact solutions of the equations of motion, and therefore, DNVMs exist at any vibration amplitude and do not depend on interaction potentials. For the first time, modulation instability of four one-component three-dimensional DNVMs is studied in a single crystal of fcc copper with the use of methods of molecular dynamics. DNVMs frequencies, evolution of stresses, kinetic and potential energies, and heat capacity depending on the oscillation amplitudes are analyzed. It is found that all four DNVMs are characterized by a hard-type anharmonicity. Modulation instability of DNVMs results in a formation of chaotic discrete breathers (DBs) with frequency above the upper edge of the phonon spectrum of the crystal. The lifetime of chaotic DBs is found to be in the range of 30-100 ps. At low-oscillation frequencies, longer-lived DBs are formed. The growth of modulation instability leads to an increase in mechanical stresses and a decrease in the heat capacity of the crystal. The results obtained in this work enrich our understanding of the influence of the modulation instability of DNVMs on the properties of metals.

Anisotropic thermal conductivity of the nanoparticles embedded GaSb thin film semiconductor

The prior theoretical model shows that GaSb is one of the few non-alloy semiconductors showing phonons ballistic effect in the thermal conductivity. However, no previous literature had been reported on the experimental measurements on the quasi-ballistic thermal transport of the GaSb thin film. In this paper, we employed the time-domain thermoreflectance (TDTR) to study the thermal transport of nanoparticles embedded GaSb thin film. Our measurements results provide first experimental evidence to verify the quasi-ballistic effect in the thermal transport of the GaSb thin film. The apparent cross-plane thermal conductivity of pure GaSb sample drops ∼15% when the pump laser modulation frequency is increased from 0.8 MHz to 10 MHz at room temperature. To further understand the thermal transport mechanism, Tempered Lévy analysis is employed to study the quasi-ballistic effect of the GaSb thin film. The model shows that GaSb thin film thermal transport has a superdiffusion exponent, [Formula: see text] = 1.51 ± 0.23 and Lévy-Fourier transition length, r _LF = 0.19 ± 0.13 µm. Both obtained values via Tempered Lévy indicates the quasi-ballistic transport phenomena in GaSb thin film. However, this frequency dependence of the cross-plane thermal conductivity will disappear in the presence of the 3%-20% ErSb nanoparticles. Another thermal transport mechanism, i.e. anisotropic thermal transport, can be observed in GaSb thin film. The ratio of in- to cross-plane thermal conductivity varies from ∼0.2 to ∼0.7 in the 0%-20% ErSb nanoparticles volume concentrations. Detailed temperature dependence of the in-plane thermal conductivity of ErSb:GaSb samples with 0%-20% are also included in the paper for the understanding of the scattering mechanism in the thin film thermal transport. With enhanced understanding of the quasi-ballistic and anisotropic thin film thermal transport, our results might improve the thermal management efficiency of the GaSb devices.

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.

Heat conduction in two-dimensional momentum-conserving and -nonconserving gases

Compared to that for two-dimensional (2D) lattices, our understanding of heat conduction in 2D gases is still limited. Here we study heat conduction behavior of 2D gas systems with momentum-conserving and -nonconserving interparticle interactions by using the nonequilibrium and equilibrium molecular dynamics methods. For the momentum-conserving system, we find that when the dimensionality of the system is changed from 2D to quasi-one-dimensional (quasi-1D), the heat conductivity κ diverges with the system size L as κ∼lnL (the theoretical prediction for 2D systems) for a short L and shows, in the thermodynamic limit, a tendency to κ∼L^{1/3} like that predicted in 1D fluids. This suggests that the dimensional-crossover effect of heat conduction exists in 2D systems with conserved momentum. In contrast, for the momentum-nonconserving system, as L increases, finite heat conductivity independent of L is observed. These findings are in agreement with the predictions given by hydrodynamic theory and thus further confirm the validity of the theory in 2D gases.

Universal scaling for recovery of Fourier's law in low-dimensional solids under momentum conservation

Dynamic renormalization group (RG) of fluctuating viscoelastic equations is investigated to clarify the cause for numerically reported disappearance of anomalous heat conduction (recovery of Fourier's law) in low-dimensional momentum-conserving systems. RG flow is obtained explicitly for simplified two model cases: a one-dimensional continuous medium under low pressure and incompressible viscoelastic medium of arbitrary dimensions. Analyses of these clarify that the inviscid fixed point of contributing the anomalous heat conduction becomes unstable under the RG flow of nonzero elastic-wave speeds. The dynamic RG analysis further predicts a universal scaling of describing the crossover between the growth and saturation of observed heat conductivity, which is confirmed through the numerical experiments of Fermi-Pasta-Ulam β (FPU-β) lattices.

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.

Geometry-induced local thermal current from cold to hot in a classical harmonic system

The second law of thermodynamics requires the overall thermal current to flow from hot to cold. However, it does not forbid a local thermal current from flowing from cold to hot. By coupling a harmonic system of three masses connected by a few springs to two Langevin reservoirs at different temperatures, a local atypical thermal current is found to flow from cold to hot in the steady state while the overall thermal current is still from hot to cold. The direction of the local thermal current can be tuned by the mass, spring constant, and system-reservoir coupling. The local thermal current can vanish if the parameters are tuned to proper values. We also consider nonlinear effect from the system-substrate coupling and find that the local atypical thermal current survives in the presence of the nonlinear potential. Moreover, the local atypical thermal current is robust against asymmetry of the system-reservoir coupling, inhomogeneity of the nonlinear potential, and additions of more masses and springs. In molecular or nanomechanical systems where the setup may find its realization, the direction of the local thermal current may be controlled by mechanical or electromagnetic means, which may lead to applications in information storage.

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

One-dimensional anomalous heat propagation is usually characterized by a Lévy walk superdiffusive spreading function with two side peaks located on the fronts due to the finite velocity of acoustic phonons. In the case when the acoustic phonons vanish, e.g., due to the phonon-lattice interactions such that the system's momentum is not conserved, the side peaks will disappear and a normal Gaussian diffusive heat-propagating behavior will be observed. Here we show that there exists another new type of superdiffusive, non-Gaussian heat propagation but without side peaks in a typical nonacoustic, momentum-nonconserving system. It implies that thermal transport in this system disobeys the Fourier law, in clear contrast with the existing theoretical predictions. The underlying mechanism is related to an effect of optical phonon-phonon interactions. These findings may open a new avenue for further exploring thermal transport in low dimensions.

Discrete breathers assist energy transfer to ac-driven nonlinear chains

A one-dimensional chain of pointwise particles harmonically coupled with nearest neighbors and placed in sixth-order polynomial on-site potentials is considered. The power of the energy source in the form of single ac driven particle is calculated numerically for different amplitudes A and frequencies ω within the linear phonon band. The results for the on-site potentials with hard and soft anharmonicity types are compared. For the hard-type anharmonicity, it is shown that when the driving frequency is close to (far from) the upper edge of the phonon band, the power of the energy source normalized to A^{2} increases (decreases) with increasing A. In contrast, for the soft-type anharmonicity, the normalized power of the energy source increases (decreases) with increasing A when the driving frequency is close to (far from) the lower edge of the phonon band. Our further demonstrations indicate that in the case of hard (soft) anharmonicity, the chain can support movable discrete breathers (DBs) with frequencies above (below) the phonon band. It is the energy source quasiperiodically emitting moving DBs in the regime with driving frequency close to the DB frequency that induces the increase of the power. Therefore, our results here support the mechanism that the moving DBs can assist energy transfer from the ac driven particle to the chain.