Change of entropy for the one-dimensional ballistic heat equation: Sinusoidal initial perturbation

Acoustic transparency of the chain-chain interface

We study propagation of wave packets through the interface between two dissimilar harmonic chains with on-site potentials (e.g., chains lying on elastic foundations). An expression for the transmission coefficient, relating energies of the incident and transmitted wave packets is derived using two different approaches. Without elastic foundation, the transmission coefficient monotonically decreases with increasing wave frequency. We show that by adding elastic foundations, one can qualitatively change this dependence and make it nonmonotonic or even increasing. Moreover, in some cases, the interface is totally transparent (the transmission coefficient is equal to unity at some frequency) if at least one of the chains has the elastic foundation. Presented results may serve for manipulation of the transmission coefficient and corresponding interfacial thermal resistance in low-dimensional nanosystems.

Transition to thermal equilibrium in a crystal subjected to instantaneous deformation

An adiabatic transition between two equilibrium states corresponding to different stiffnesses in an infinite chain of particles is studied. Initially, the particles have random displacements and random velocities corresponding to uniform initial temperature distributions. An instantaneous change in the parameters of the chain initiates a transitional process. Analytical expressions for the chain temperature as a function of time are obtained from statistical analysis of the dynamic equations. It is shown that the transition process is oscillatory and that the temperature converges non-monotonically to a new equilibrium state, in accordance with what is usually unexpected for thermal processes. The analytical results are supplemented by numerical simulations.

Geometrically nonlinear dynamic model for a hexagonal lattice

It is shown that angular stiffness in the hexagonal lattice model plays a significant role in the geometrical nonlinear terms in the equations of the continuum limit. A geometrically nonlinear discrete model is formulated for the hexagonal lattice by considering the interaction of two sublattices. An asymptotic procedure is developed in order to obtain the nonlinear coupled equations of motion in the continuum limit of the discrete model. An interaction of longitudinal and shear plane strain waves is studied by using the solutions of the obtained equations.

Ballistic resonance and thermalization in the Fermi-Pasta-Ulam-Tsingou chain at finite temperature

We study conversion of thermal energy to mechanical energy and vice versa in an α-Fermi-Pasta-Ulam-Tsingou (FPUT) chain with a spatially sinusoidal profile of initial temperature. We show analytically that coupling between macroscopic dynamics and quasiballistic heat transport gives rise to mechanical vibrations with growing amplitude. This phenomenon is referred to as ballistic resonance. At large times, these mechanical vibrations decay monotonically, and therefore the well-known FPUT recurrence paradox occurring at zero temperature is eliminated at finite temperatures.

Equilibration of energies in a two-dimensional harmonic graphene lattice

We study dynamical phenomena in a harmonic graphene (honeycomb) lattice, consisting of equal particles connected by linear and angular springs. Equations of in-plane motion for the lattice are derived. Initial conditions typical for molecular dynamic modelling are considered. Particles have random initial velocities and zero displacements. In this case, the lattice is far from thermal equilibrium. In particular, initial kinetic and potential energies are not equal. Moreover, initial kinetic energies (and temperatures), corresponding to degrees of freedom of the unit cell, are generally different. The motion of particles leads to equilibration of kinetic and potential energies and redistribution of kinetic energy among degrees of freedom. During equilibration, the kinetic energy performs decaying high-frequency oscillations. We show that these oscillations are accurately described by an integral depending on dispersion relation and polarization matrix of the lattice. At large times, kinetic and potential energies tend to equal values. Kinetic energy is partially redistributed among degrees of freedom of the unit cell. Equilibrium distribution of the kinetic energies is accurately predicted by the non-equipartition theorem. Presented results may serve for better understanding of the approach to thermal equilibrium in graphene. This article is part of the theme issue 'Modelling of dynamic phenomena and localization in structured media (part 2)'.

Thermal equilibration in a one-dimensional damped harmonic crystal

The features for the unsteady process of thermal equilibration ("the fast motions") in a one-dimensional harmonic crystal lying in a viscous environment (e.g., a gas) are under investigation. It is assumed that initially the displacements of all the particles are zero and the particle velocities are random quantities with zero mean and a constant variance, thus, the system is far away from the thermal equilibrium. It is known that in the framework of the corresponding conservative problem the kinetic and potential energies oscillate and approach the equilibrium value that equals a half of the initial value of the kinetic energy. We show that the presence of the external damping qualitatively changes the features of this process. The unsteady process generally has two stages. At the first stage oscillations of kinetic and potential energies with decreasing amplitude, subjected to exponential decay, can be observed (this stage exists only in the underdamped case). At the second stage (which always exists), the oscillations vanish, and the energies are subjected to a power decay. The large- time asymptotics for the energy is proportional to t^{-3/2} in the case of the potential energy and to t^{-5/2} in the case the kinetic energy. Hence, at large values of time the total energy of the crystal is mostly the potential energy. The obtained analytic results are verified by independent numerical calculations.