Normal heat conduction in one-dimensional momentum conserving lattices with asymmetric interactions

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.

Too Close to Integrable: Crossover from Normal to Anomalous Heat Diffusion

Energy transport in one-dimensional chains of particles with three conservation laws is generically anomalous and belongs to the Kardar-Parisi-Zhang dynamical universality class. Surprisingly, some examples where an apparent normal heat diffusion is found over a large range of length scales were reported. We propose a novel physical explanation of these intriguing observations. We develop a scaling analysis that explains how this may happen in the vicinity of an integrable limit, such as, but not only, the famous Toda model. In this limit, heat transport is mostly supplied by quasiparticles with a very large mean free path ℓ. Upon increasing the system size L, three different regimes can be observed: a ballistic one, an intermediate diffusive range, and, eventually, the crossover to the anomalous (hydrodynamic) regime. Our theoretical considerations are supported by numerical simulations of a gas of diatomic hard-point particles for almost equal masses and of a weakly perturbed Toda chain. Finally, we discuss the case of the perturbed harmonic chain, which exhibits a yet different scenario.

Energy and spin diffusion in the one-dimensional classical Heisenberg spin chain at finite and infinite temperatures

The energy and spin diffusion behaviors in the one-dimensional classical Heisenberg spin chain have been systematically investigated using the equilibrium diffusion method. The spatiotemporal autocorrelation functions for energy and spin are calculated at finite and infinite temperatures. As conserved quantities, the spreading of excess energy and spin can be used to determine their actual diffusion behaviors. At low temperatures, the energy diffusion shows almost ballistic behavior, and spin shows superdiffusion behavior for finite chain size. For energy diffusion, normal diffusion behavior can be obtained when the temperature is higher than 0.75. For spin diffusion, normal diffusion behavior is observed at infinite temperature.

Energy and momentum diffusion in one-dimensional periodic and asymmetric nonlinear lattices with momentum conservation

The energy diffusion in one-dimensional (1D) momentum conserving nonlinear lattices usually exhibits anomalous superdiffusion, except the coupled rotator lattice with symmetric and periodic interacting potential which has normal energy diffusion corresponding to normal heat conduction. For nonperiodic 1D lattices with momentum conservation, it has been argued that the asymmetric potential can induce normal heat conduction. Later results indicate the observed normal behavior might be the finite size effect and the anomalous behavior will appear in the thermodynamical limit. Here we propose asymmetric and periodic 1D nonlinear lattices with momentum conservation. The energy and momentum diffusion behaviors will be investigated in detail and the same normal diffusion behaviors for both energy and momentum can be observed. These results confirm that the periodicity is the key for normal transport behavior in 1D momentum conserving lattices, whether the potential is symmetric or asymmetric.

Temperature-dependent divergence of thermal conductivity in momentum-conserving one-dimensional lattices with asymmetric potential

In this study we used a nonequilibrium simulation method to investigate the temperature dependent divergence of thermal conductivity in a one-dimensional momentum conserving system with an asymmetric double well nearest-neighbor interaction potential. We show that across all temperatures thermal conductivity exhibits power-law divergence with the chain length and the value of the divergence exponent (α) depends on the temperature of the system. At low and high temperatures α reaches close to ∼0.5 and ∼0.33, respectively. Whereas in the intermediate temperature the divergence of thermal conductivity with the chain length saturates with α∼0.07. Subsequent analysis showed that the estimated value of α in the intermediate temperature may not have reached its thermodynamic limit. Further calculations of local α revealed that its approach towards the thermodynamic limit is crucially dependent on the temperature of the system. At low and high temperatures local α reaches its thermodynamic limits in shorter chain lengths. On the contrary, in the case of intermediate temperature its progress towards the asymptotic limit is nonmonotonic.

Solitons as candidates for energy carriers in Fermi-Pasta-Ulam lattices

Currently, effective phonons (renormalized or interacting phonons) rather than solitary waves (for short, solitons) are regarded as the energy carriers in nonlinear lattices. In this work, by using the approximate soliton solutions of the corresponding equations of motion and adopting the Boltzmann distribution for these solitons, the average velocities of solitons are obtained and are compared with the sound velocities of energy transfer. Excellent agreements with the numerical results and the predictions of other existing theories are shown in both the symmetric Fermi-Pasta-Ulam-β lattices and the asymmetric Fermi-Pasta-Ulam-αβ lattices. These clearly indicate that solitons are suitable candidates for energy carriers in Fermi-Pasta-Ulam lattices. In addition, the root-mean-square velocity of solitons can be obtained from the effective phonons theory.

Observing golden-mean universality class in the scaling of thermal transport

We address the issue of whether the golden-mean [ψ=(sqrt[5]+1)/2≃1.618] universality class, as predicted by several theoretical models, can be observed in the dynamical scaling of thermal transport. Remarkably, we show strong evidence that ψ appears to be the scaling exponent of heat mode correlation in a purely quartic anharmonic chain. This observation seems to somewhat deviate from the previous expectation and we explain it by the unusual slow decay of the cross correlation between heat and sound modes. Whenever the cubic anharmonicity is included, this cross correlation gradually dies out and another universality class with scaling exponent γ=5/3, as commonly predicted by theories, seems recovered. However, this recovery is accompanied by two interesting phase transition processes characterized by a change of symmetry of the potential and a clear variation of the dynamic structure factor, respectively. Due to these transitions, an additional exponent close to γ≃1.580 emerges. All this evidence suggests that, to gain a full prediction of the scaling of thermal transport, more ingredients should be taken into account.

Response and correlation functions of nonlinear systems in equilibrium states

In this paper, we study systematically a serial of correlation functions in some one-dimensional nonlinear lattices. Due to the energy conservation law, they are implicitly interdependent. Various transport coefficients are thus also connected. In the studies of the autocorrelations of local energy density and of local heat current, a general relation between diverging heat conduction and super heat diffusion has been proposed recently. We clarify that such a relation is valid only in systems without temperature pressure. In those with temperature pressure, a constant but nontrivial term appears. This term explains a previously observed fact that heat diffusion in such systems is always ballistic but heat conduction can diverge very slowly. Such a result not only disproves the existence of any general relation between diverging heat conduction and super heat diffusion, but it also breaks the long-term presumption that ballistic heat conduction and diffusion always coexist.

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

Anomalous (non-Fourier) heat transport is no longer just a theoretical issue since it has been observed experimentally in a number of low-dimensional nanomaterials, such as SiGe nanowires, carbon nanotubes, and others. To understand these anomalous behaviors, exploring the microscopic origin of normal (Fourier) heat transport is a fascinating theoretical topic. However, this issue has not yet been fully understood even for one-dimensional (1D) model chains, in spite of a great amount of thorough studies done to date. From those studies, it has been widely accepted that the conservation of momentum is a key ingredient to induce anomalous heat transport, while momentum-nonconserving systems usually support normal heat transport where Fourier's law is valid. But if the nonconservation of momentum is the reason, what is the underlying microscopic mechanism for the observed normal heat transport? Here we carefully revisit a typical 1D momentum-nonconserving ϕ^{4} model, and we present evidence that the mobile discrete breathers, or, in other words, the moving intrinsic localized modes with frequency components above the linear phonon band, can be responsible for that.

Scattering of lattice solitons and decay of heat-current correlation in the Fermi-Pasta-Ulam-α-β model

As is well known, solitons can be excited in nonlinear lattice systems; however, whether, and if so, how, this kind of nonlinear excitation can affect the energy transport behavior is not yet fully understood. Here we study both the scattering dynamics of solitons and heat transport properties in the Fermi-Pasta-Ulam-α-β model with an asymmetric interparticle interaction. By varying the asymmetry degree of the interaction (characterized by α), we find that (i) for each α there exists a momentum threshold for exciting solitons from which one may infer an α-dependent feature of probability of presentation of solitons at a finite-temperature equilibrium state and (ii) the scattering rate of solitons is sensitively dependent on α. Based on these findings, we conjecture that the scattering between solitons will cause the nonmonotonic α-dependent feature of heat conduction. Fortunately, such a conjecture is indeed verified by our detailed examination of the time decay behavior of the heat current correlation function, but it is only valid for an early time stage. Thus, this result may suggest that solitons can have only a relatively short survival time when exposed in a thermal environment, eventually affecting the heat transport in a short time.

Long-time instability in the Runge-Kutta algorithm for a Nosé-Hoover heat bath of a harmonic chain and its stabilization

In this paper, we investigate the Runge-Kutta algorithm for the Nosé-Hoover heat bath of a harmonic chain. The Runge-Kutta algorithm is found to be unstable in long-time calculations, with the system temperature growing exponentially. The growth rate increases if time step size is chosen larger. By analyzing the Fourier spectra in both space (wave number) and time (frequency), we discover that the growth is caused by spurious energy accumulation, particularly at the largest wave number. Such accumulation may be explained by von Neumann analysis for an infinite chain, with the nonlinear heat bath being ignored. Furthermore, we propose to add a filter to remove excessive energy, which effectively stabilizes the algorithm.

Anomalous temperature-dependent heat transport in one-dimensional momentum-conserving systems with soft-type interparticle interaction

We numerically investigate the heat transport problem in a one-dimensional momentum-conserving lattice with a soft-type (ST) anharmonic interparticle interaction. It is found that with the increase of the system's temperature, while the introduction of ST anharmonicity softens phonons and decreases their velocities, this type of nonlinearity like its hard type (HT) counterpart, can still not be able to fully damp the longest wavelength phonons. Therefore, a usual anomalous temperature dependence of heat transport with certain scaling properties similarly to those shown in the Fermi-Pasta-Ulam-β-like systems with HT interactions can be seen. Our detailed examination from simulations verifies this temperature-dependent behavior well.

Heat conduction in a chain of colliding particles with a stiff repulsive potential

One-dimensional billiards, i.e., a chain of colliding particles with equal masses, is a well-known example of a completely integrable system. Billiards with different particle masses is generically not integrable, but it still exhibits divergence of a heat conduction coefficient (HCC) in the thermodynamic limit. Traditional billiards models imply instantaneous (zero-time) collisions between the particles. We relax this condition of instantaneous impact and consider heat transport in a chain of stiff colliding particles with the power-law potential of the nearest-neighbor interaction. The instantaneous collisions correspond to the limit of infinite power in the interaction potential; for finite powers, the interactions take nonzero time. This modification of the model leads to a profound physical consequence-the probability of multiple (in particular triple) -particle collisions becomes nonzero. Contrary to the integrable billiards of equal particles, the modified model exhibits saturation of the heat conduction coefficient for a large system size. Moreover, the identification of scattering events with triple-particle collisions leads to a simple definition of the characteristic mean free path and a kinetic description of heat transport. This approach allows us to predict both the temperature and density dependencies for the HCC limit values. The latter dependence is quite counterintuitive-the HCC is inversely proportional to the particle density in the chain. Both predictions are confirmed by direct numerical simulations.

Pressure-induced recovery of Fourier's law in one-dimensional momentum-conserving systems

We report the two typical models of normal heat conduction in one-dimensional momentum-conserving systems. They show the Arrhenius and the non-Arrhenius temperature dependence. We construct the two corresponding phenomenologies, transition-state theory of thermally activated dissociation and the pressure-induced crossover between two fixed points in fluctuating hydrodynamics. Compressibility yields the ballistic fixed point, whose scaling is observed in Fermi-Pasta-Ulam (FPU) β lattices.

Heat conduction and energy diffusion in momentum-conserving one-dimensional full-lattice ding-a-ling model

The ding-a-ling model is a kind of half lattice and half hard-point-gas (HPG) model. The original ding-a-ling model proposed by Casati et al. does not conserve total momentum and has been found to exhibit normal heat conduction behavior. Recently, a modified ding-a-ling model which conserves total momentum has been studied and normal heat conduction has also been claimed. In this work, we propose a full-lattice ding-a-ling model without hard point collisions where total momentum is also conserved. We investigate the heat conduction and energy diffusion of this full-lattice ding-a-ling model with three different nonlinear inter-particle potential forms. For symmetrical potential lattices, the thermal conductivities diverges with lattice length and their energy diffusions are superdiffusive signaturing anomalous heat conduction. For asymmetrical potential lattices, although the thermal conductivity seems to converge as the length increases, the energy diffusion is definitely deviating from normal diffusion behavior indicating anomalous heat conduction as well. No normal heat conduction behavior can be found for the full-lattice ding-a-ling model.

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.

Anomalous dynamical scaling in anharmonic chains and plasma models with multiparticle collisions

We study the anomalous dynamical scaling of equilibrium correlations in one-dimensional systems. Two different models are compared: the Fermi-Pasta-Ulam chain with cubic and quartic nonlinearity and a gas of point particles interacting stochastically through multiparticle collision dynamics. For both models-that admit three conservation laws-by means of detailed numerical simulations we verify the predictions of nonlinear fluctuating hydrodynamics for the structure factors of density and energy fluctuations at equilibrium. Despite this, violations of the expected scaling in the currents correlation are found in some regimes, hindering the observation of the asymptotic scaling predicted by the theory. In the case of the gas model this crossover is clearly demonstrated upon changing the coupling constant.

Deviation from the Maxwell-Cattaneo law: Role of asymmetric interparticle interactions

Nonstationary heat conduction in a few one-dimensional nonlinear lattices is studied numerically based on the Maxwell-Cattaneo (MC) law. We simulate the relaxation process and calculate the magnitudes of the temperature oscillation A(T)(t) and the local heat current oscillation A(j)(t). A phase difference between A(T)(t) and A(j)(t) is observed, which not only verifies the existence of the time lag τ in the MC law but also provides a better way of determining the critical wavelength L(*) that separates between oscillatory and diffusive relaxation modes. However, clear deviations from the MC law are observed. Not only do the decay exponents differ from the theoretical expectations, but, more importantly, suboscillation in the diffusive regime, which is not expected by the MC law, is found in the lattices with asymmetric interactions as well. These findings imply that higher-order effects must be considered in order to well describe the nonstationary heat conduction process in these systems.

Effects of interaction symmetry on delocalization and energy transport in one-dimensional disordered lattices

We study effects of interaction symmetry in one-dimensional, momentum-conserving disordered lattices. It is found that asymmetric and symmetric interparticle interactions may result in significant difference: localized modes can be delocalized by very weak asymmetric interactions but survive much stronger symmetric interactions. Moreover, in the delocalization regime, asymmetric and symmetric interactions also have qualitatively different effects on transport: the former (the latter) may lead to a fast decaying (slow power-law decaying) heat current correlation function and in turn a convergent (divergent) heat conductivity. A method for detecting delocalization in systems at a nonzero temperature is proposed as well.

Thermoelectricity of interacting particles: a numerical approach

A method for computing the thermopower in interacting systems is proposed. This approach, which relies on Monte Carlo simulations, is illustrated first for a diatomic chain of hard-point elastically colliding particles and then in the case of a one-dimensional gas with (screened) Coulomb interparticle interaction. Numerical simulations up to N>10^{4} particles confirm the general theoretical arguments for momentum-conserving systems and show that the thermoelectric figure of merit increases linearly with the system size.

Super heat diffusion in one-dimensional momentum-conserving nonlinear lattices

Heat diffusion processes in various one-dimensional total-momentum-conserving nonlinear lattices with symmetric interaction and asymmetric interaction are systematically studied. It is revealed that the asymmetry of interaction largely enhances the heat diffusion; while according to our existing studies for heat conduction in the same lattices, it slows the divergence of heat conductivity in a wide regime of system size. These findings violate the proposed relations that connect anomalous heat conduction and super heat diffusion. The generality of those expectations is thus questioned.

Renormalized phonons in nonlinear lattices: A variational approach

We propose a variational approach to study renormalized phonons in momentum-conserving nonlinear lattices with either symmetric or asymmetric potentials. To investigate the influence of pressure for phonon properties, we derive an inequality which provides both the lower and upper bound of the Gibbs free energy as the associated variational principle. This inequality is a direct extension to the Gibbs-Bogoliubov inequality. Taking the symmetry effect into account, the reference system for the variational approach is chosen to be harmonic with an asymmetric quadratic potential which contains variational parameters. We demonstrate the power of this approach by applying it to one-dimensional nonlinear lattices with a symmetric or asymmetric Fermi-Pasta-Ulam-type potential. For a system with a symmetric potential and zero pressure, we recover existing results. For other systems which are beyond the scope of existing theories, including those having symmetric potential and pressure and those having the asymmetric potential with or without pressure, we also obtain accurate sound velocity.

Temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with a soft on-site potential

The temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with soft on-site potential (soft-KG) are investigated systematically. Similarly to the previously studied hard-KG lattices, the existence of renormalized phonons is also confirmed in soft-KG lattices. In particular, the temperature dependence of the renormalized phonon frequency predicted by a classical field theory is verified by detailed numerical simulations. However, the thermal conductivities of soft-KG lattices exhibit the opposite trend in temperature dependence in comparison with those of hard-KG lattices. The interesting thing is that the temperature-dependent thermal conductivities of both soft- and hard-KG lattices can be interpreted in the same framework of effective phonon theory. According to the effective phonon theory, the exponents of the power-law dependence of the thermal conductivities as a function of temperature are only determined by the exponents of the soft or hard on-site potentials. These theoretical predictions are consistently verified very well by extensive numerical simulations.

Localization of low-frequency oscillations in single-walled carbon nanotubes

In the framework of the continuum shell theory, we analytically predict a new phenomenon: the weak localization of optical low-frequency oscillations in carbon nanotubes. We clarify the origin of the localization by means of the concept of the limiting phase trajectory and confirm the obtained analytical results by molecular dynamics simulations of simply supported carbon nanotubes. The performed analysis contributes to the new universal approach to the treatment of nonstationary resonant processes.

Nonintegrability and the Fourier heat conduction law

We study in momentum-conserving systems, how nonintegrable dynamics may affect thermal transport properties. As illustrating examples, two one-dimensional (1D) diatomic chains, representing 1D fluids and lattices, respectively, are numerically investigated. In both models, the two species of atoms are assigned two different masses and are arranged alternatively. The systems are nonintegrable unless the mass ratio is one. We find that when the mass ratio is slightly different from one, the heat conductivity may keep significantly unchanged over a certain range of the system size and as the mass ratio tends to one, this range may expand rapidly. These results establish a new connection between the macroscopic thermal transport properties and the underlying dynamics.

Temperature dependence of heat conduction in the Fermi-Pasta-Ulam-β lattice with next-nearest-neighbor coupling

We show numerically that introducing the next-nearest-neighbor interactions (of appropriate strength) into the one-dimensional (1D) Fermi-Pasta-Ulam-β (FPU-β) lattice can result in an unusual, nonmonotonic temperature dependent divergence behavior in a wide temperature range, which is in clear contrast to the universal divergence manner independent of temperature as suggested previously in the conventional 1D FPU-β models with nearest-neighbor (NN) coupling only. We also discuss the underlying mechanism of this finding by analyzing the temperature variations of the properties of discrete breathers, especially that with frequencies having the intraband components. The results may provide useful information for establishing the connection between the macroscopic heat transport properties and the underlying dynamics in general 1D systems with interactions beyond NN couplings.

Numerical test of hydrodynamic fluctuation theory in the Fermi-Pasta-Ulam chain

Recent work has developed a nonlinear hydrodynamic fluctuation theory for a chain of coupled anharmonic oscillators governing the conserved fields, namely, stretch, momentum, and energy. The linear theory yields two propagating sound modes and one diffusing heat mode, all three with diffusive broadening. In contrast, the nonlinear theory predicts that, at long times, the sound mode correlations satisfy Kardar-Parisi-Zhang scaling, while the heat mode correlations have Lévy-walk scaling. In the present contribution we report on molecular dynamics simulations of Fermi-Pasta-Ulam chains to compute various spatiotemporal correlation functions and compare them with the predictions of the theory. We obtain very good agreement in many cases, but also some deviations.

Cross-correlations between phonon modes in anharmonic oscillator chains: role in heat transport

We have computed current-current correlation functions in chains of anharmonic oscillators described by various models (FPU-β, FPU-αβ, ϕ4), considering both the total current and the currents associated with individual phonon modes, which are important in view of the Green-Kubo relation for heat conductivity. Our simulations show that, contrary to the common hypothesis, there are, under some circumstances, significant correlations between neighboring modes. These cross-mode correlations are the dominant contribution to the conductivity in the low anharmonicity regime. The inverse of the timescale over which they are significant, 1/τc, is related to the anharmonicity level in a way similar to the largest Lyapunov exponent, suggesting that the two quantities are related. Cross-mode correlations exist in both anomalous and regular heat-conducting systems although we are unable to observe a transition to the independent-mode regime in the latter case.

Thermal conductivity of molecular chains with asymmetric potentials of pair interactions

We provide molecular-dynamics simulation of heat transport in one-dimensional molecular chains with different interparticle pair potentials. We show that the thermal conductivity is finite in the thermodynamic limit in chains with the potentials that allow for bond dissociation. The Lennard-Jones, Morse, and Coulomb potentials are such potentials. The convergence of the thermal conductivity is provided by phonon scattering on the locally strongly stretched loose interatomic bonds at low temperature and by the many-particle scattering at high temperature. On the other hand, chains with a confining pair potential, which does not allow for bond dissociation, possess anomalous thermal conductivity, diverging with the chain length. We emphasize that chains with a symmetric or asymmetric Fermi-Pasta-Ulam potential or with combined potentials, containing a parabolic and/or a quartic confining potential, all exhibit anomalous heat transport.

Finite-size effects on current correlation functions

We study why the calculation of current correlation functions (CCFs) still suffers from finite-size effects even when the periodic boundary condition is taken. Two important one-dimensional, momentum-conserving systems are investigated as examples. Intriguingly, it is found that the state of a system recurs in the sense of microcanonical ensemble average, and such recurrence may result in oscillations in CCFs. Meanwhile, we find that the sound mode collisions induce an extra time decay in a current so that its correlation function decays faster (slower) in a smaller (larger) system. Based on these two unveiled mechanisms, a procedure for correctly evaluating the decay rate of a CCF is proposed, with which our analysis suggests that the global energy CCF decays as ∼ t(-2/3) in the diatomic hard-core gas model and in a manner close to ∼ t(-1/2) in the Fermi-Pasta-Ulam-β model.

Mechanical control of heat conductivity in molecular chains

We discuss a possibility to control heat conductivity in molecular chains by means of external mechanical loads. To illustrate such possibilities we consider first well-studied one-dimensional chain with degenerate double-well potential of the nearest-neighbor interaction. We consider varying lengths of the chain with fixed number of particles. Number of possible energetically degenerate ground states strongly depends on the overall length of the chain, or, in other terms, on average length of the link between neighboring particles. These degenerate states correspond to mechanical equilibria; therefore, one can say that formation of such structures mimics a process of plastic deformation. We demonstrate that such modification of the chain length can lead to quite profound (almost fivefold) reduction of the heat conduction coefficient. Even more profound effect is revealed for a model with a single-well nonconvex potential. It is demonstrated that in a certain range of constant external forcing, this model becomes effectively double-well and has a multitude of possible states of equilibrium for fixed value of the external load. Due to this degeneracy, the heat-conduction coefficient can be reduced by two orders of magnitude. We suggest a mechanical model of a chain with periodic double-well potential, which allows control of the heat transport. The models considered may be useful for description of heat transfer in biological macromolecules and for control of the heat transport in microsystems. The possibility of the heat transport control in more realistic three-dimensional systems is illustrated by simulation of a three-dimensional model of polymer α-helix. In this model, the mechanical stretching also brings about the structural inhomogeneity and, in turn, to essential reduction of the heat conductivity.

Validity of Fourier's law in one-dimensional momentum-conserving lattices with asymmetric interparticle interactions

We have numerically studied heat conduction in a few one-dimensional momentum-conserving lattices with asymmetric interparticle interactions by the nonequilibrium heat bath method, the equilibrium Green-Kubo method, and the heat current power spectra analysis. Very strong finite-size effects are clearly observed. Such effects make the heat conduction obey a Fourier-like law in a wide range of lattice lengths. However, in yet longer lattice lengths, the heat conductivity regains its power-law divergence. Therefore, the power-law divergence of the heat conductivity in the thermodynamic limit is verified, as is expected by many existing theories.

Heat transport enhanced by optical phonons in one-dimensional anharmonic lattices with alternating bonds

In lattice systems, the effects of optical phonons on heat transport are usually neglected due to their relatively small group velocities compared with acoustic phonons, or even assumed to be negative because introducing optical phonons may simultaneously reduce the group velocities of acoustic phonons. In order to well understand the role played by optical phonons, we propose a one-dimensional anharmonic lattice model with alternating interactions, where the optical phonons can be conveniently tuned. We find that in contrast to previous studies, the optical phonons (in coordination with the nonlinearities) can enhance heat transport in the thermodynamical limit, suggesting that optical phonons can also play an active role. The underlying mechanism is related to the effects of two kinds of nonlinear excitations, i.e., the optical and the gap discrete breathers (DBs). These DBs release energy and in turn facilitate heat transport.