Transport in the Frenkel-Kontorova model. III. Thermal conductivity

Approaching the perfect diode limit through a nonlinear interface

We consider a system formed by two different segments of particles, coupled to thermal baths, one at each end, modeled by Langevin thermostats. The particles in each segment interact harmonically and are subject to an on-site potential for which three different types are considered, namely, harmonic, ϕ^{4}, and Frenkel-Kontorova. The two segments are nonlinearly coupled, between interfacial particles, by means of a power-law potential with exponent μ, which we vary, scanning from subharmonic to superharmonic potentials, up to the infinite-square-well limit (μ→∞). Thermal rectification is investigated by integrating the equations of motion and computing the heat fluxes. As a measure of rectification, we use the difference of the currents, resulting from the interchange of the baths, divided by their average (all quantities taken in absolute value). We find that rectification can be optimized by a given value of μ that depends on the bath temperatures and details of the chains. But, regardless of the type of on-site potential considered, the interfacial potential that produces maximal rectification approaches the infinite square well (μ→∞) when reducing the average temperature of the baths. Our analysis of thermal rectification focuses on this regime, for which we complement numerical results with heuristic considerations.

High-heat-flux rectification due to a localized thermal diode

A theoretical implementation of a localized thermal diode with a rectification factor greater than 10^{6} is demonstrated. In reverse thermal bias, extremely low thermal conductivity is achieved through phononic Rayleigh scattering from a finite-depth defect. In forward bias, the diode oscillator escapes the defect and thermal conductivity becomes up to four orders of magnitude higher. The setup provides a minimal model of a localized thermal diode between two identical oscillator chains and opens up a pathway for thermal diode implementations.

Phonon coupling and transport in individual polyethylene chains: a comparison study with the bulk crystal

Using the first-principles-based anharmonic lattice dynamics, we calculate the thermal conductivities (κ) of both bulk and single-chain polyethylene (PE) and characterize the mode-wise phonon transport and scattering channels. A significantly higher room-temperature axial thermal conductivity in single-chain PE (1400 W m^-1 K^-1) is observed compared to bulk PE crystals (237 W m^-1 K^-1). The reduction of scattering phase space caused by the diminished inter-chain van der Waals interactions explains the much larger κ in single-chain PE. Different from many previous studies, the thermal conductivity of single-chain PE is predicted to converge at a chain length of ∼1 mm at 300 K. The convergence is explained by the indirect thermal resistance from momentum-conserving scatterings of long-wavelength phonons. It is also found that longitudinal phonon modes dominate the thermal transport in PE chains, while transverse phonon branches with quadratic dispersions contribute little to κ due to their vanishing group velocities and limited lifetimes in the long wavelength limit. The predicted high κ of bulk crystalline and single-chain PE show great potential for use of polymers in thermal management, and the unveiled phonon transport mechanisms offer guides for their molecule-level design.

Dynamical thermalization of Frenkel-Kontorova model in the thermodynamic limit

We study numerically the process of dynamical thermalization in the Frenkel-Kontorova (FK) model with weak nonlinearity. The total energy has initially equidistributed among some of the lowest frequency linear modes. It is found that the energy transfers continuously to the high-frequency modes and finally evolves towards energy equipartition in the FK model. However, the metastable state, which was found in Fermi-Pasta-Ulam (FPU) model and φ(4) model in a relatively short time scale, is not found in the FK model. We further perform a very accurate systematic study of the equipartition time T(eq) as functions of the particle number N, the nonlinear parameter β, and the energy density ɛ. In the thermodynamic limit, the dependence of T(eq) on β and ɛ is found to display a power law behavior: T(eq)∝β(a)ɛ(b). The exponents a and b are numerically found to be approximately -2.0 and 1.43. This scaling law is also quite different from those of the FPU-β model and φ(4) model.

Boundary-induced instabilities in coupled oscillators

A novel class of nonequilibrium phase transitions at zero temperature is found in chains of nonlinear oscillators. For two paradigmatic systems, the Hamiltonian XY model and the discrete nonlinear Schrödinger equation, we find that the application of boundary forces induces two synchronized phases, separated by a nontrivial interfacial region where the kinetic temperature is finite. Dynamics in such a supercritical state displays anomalous chaotic properties whereby some observables are nonextensive and transport is superdiffusive. At finite temperatures, the transition is smoothed, but the temperature profile is still nonmonotonic.

Nonequilibrium discrete nonlinear Schrödinger equation

We study nonequilibrium steady states of the one-dimensional discrete nonlinear Schrödinger equation. This system can be regarded as a minimal model for the stationary transport of bosonic particles such as photons in layered media or cold atoms in deep optical traps. Due to the presence of two conserved quantities, namely, energy and norm (or number of particles), the model displays coupled transport in the sense of linear irreversible thermodynamics. Monte Carlo thermostats are implemented to impose a given temperature and chemical potential at the chain ends. As a result, we find that the Onsager coefficients are finite in the thermodynamic limit, i.e., transport is normal. Depending on the position in the parameter space, the "Seebeck coefficient" may be either positive or negative. For large differences between the thermostat parameters, density and temperature profiles may display an unusual nonmonotonic shape. This is due to the strong dependence of the Onsager coefficients on the state variables.

Momentum conserving one-dimensional system with a finite thermal conductivity

A one-dimensional system of particles is examined in which even numbered particles are bound to adjacent even particles by harmonic spring forces, while odd particles are free. Even and odd particles collide elastically. This is a momentum conserving modification of the famous "ding-a-ling" model. Molecular-dynamics simulations are carried out and the current power spectra are obtained. The energy current power spectrum has zero slope at low frequencies. This implies that the thermal conductivity κ is finite and independent of system length L , for L sufficiently large. Steady-state simulations provide further evidence that κ is independent of L at large values of L . The relevance of this result to the proof by Prosen and Campbell that momentum conservation with nonvanishing pressure implies an infinite thermal conductivity is discussed.

Heat conduction in the Frenkel-Kontorova model

Heat conduction is an old yet important problem. Since Fourier introduced the law bearing his name almost 200 years ago, a first-principle derivation of this simple law from statistical mechanics is still lacking. Worse still, the validity of this law in low dimensions, and the necessary and sufficient conditions for its validity are far from clear. In this paper we will review recent works on heat conduction in a simple nonintegrable model called the Frenkel-Kontorova model. The thermal conductivity of this model has been found to be finite. We will study the dependence of the thermal conductivity on the temperature and other parameters of the model such as the strength and the periodicity of the external potential. We will also discuss other related problems such as phase transitions and finite-size effects. The study of heat conduction is not only of theoretical interest but also of practical interest. We will show various recent designs of thermal rectifiers and thermal diodes by coupling nonlinear chains together. The study of heat conduction in low dimensions is also important to the understanding of the thermal properties of carbon nanotubes.

Heat conduction in a one-dimensional chain of hard disks with substrate potential

Heat conduction in a one-dimensional chain of equivalent rigid particles in the field of the external on-site potential is considered. The zero diameters of the particles correspond to the integrable case with the divergent heat conduction coefficient. By means of a simple analytical model it is demonstrated that for any nonzero particle size the integrability is violated and the heat conduction coefficient converges. The result of the analytical computation is verified by means of numerical simulation in a plausible diapason of parameters, and good agreement is observed.

Fermi-Pasta-Ulam beta lattice: Peierls equation and anomalous heat conductivity

The Peierls equation is considered for the Fermi-Pasta-Ulam beta lattice. Explicit form of the linearized collision operator is obtained. Using this form the decay rate of the normal-mode energy as a function of wave vector k is estimated to be proportional to k(5/3). This leads to the t(-3/5) long-time behavior of the current correlation function, and, therefore, to the divergent coefficient of heat conductivity. These results are in good agreement with the results of recent computer simulations. Compared to the results obtained through the mode coupling theory our estimations give the same k dependence of the decay rate but a different temperature dependence. Using our estimations we argue that adding a harmonic on-site potential to the Fermi-Pasta-Ulam beta lattice may lead to finite heat conductivity in this model.

Heat conduction in one-dimensional lattices with on-site potential

The process of heat conduction in one-dimensional lattices with on-site potential is studied by means of numerical simulation. Using the discrete Frenkel-Kontorova, phi(4), and sinh-Gordon models we demonstrate that contrary to previously expressed opinions the sole anharmonicity of the on-site potential is insufficient to ensure the normal heat conductivity in these systems. The character of the heat conduction is determined by the spectrum of nonlinear excitations peculiar for every given model and therefore depends on the concrete potential shape and the temperature of the lattice. The reason is that the peculiarities of the nonlinear excitations and their interactions prescribe the energy scattering mechanism in each model. For sine-Gordon and phi(4) models, phonons are scattered at a dynamical lattice of topological solitons; for sinh-Gordon and for phi(4) in a different parameter regime the phonons are scattered at localized high-frequency breathers (in the case of phi(4) the scattering mechanism switches with the growth of the temperature).

Fermi-Pasta-Ulam beta model: boundary jumps, Fourier's law, and scaling

We examine the interplay of surface and volume effects in systems undergoing heat flow. In particular, we compute the thermal conductivity in the Fermi-Pasta-Ulam beta model as a function of temperature and lattice size, and scaling arguments are used to provide analytic guidance. From this we show that boundary temperature jumps can be quantitatively understood, and that they play an important role in determining the dynamics of the system, relating soliton dynamics, kinetic theory, and Fourier transport.

Normal heat conductivity of the one-dimensional lattice with periodic potential of nearest-neighbor interaction

The process of heat conduction in a chain with a periodic potential of nearest-neighbor interaction is investigated by means of molecular dynamics simulation. It is demonstrated that the periodic potential of nearest-neighbor interaction allows one to obtain normal heat conductivity in an isolated one-dimensional chain with conserved momentum. The system exhibits a transition from infinite to normal heat conductivity with the growth of its temperature. The physical reason for normal heat conductivity is the excitation of high-frequency stationary localized rotational modes. These modes absorb the momentum and facilitate locking of the heat flux.