Solitary waves in nonequilibrium molecular dynamics simulations of heat flow in one-dimensional lattices

Abstract

We study the use of the Evans Nonequilibrium Molecular Dynamics (NEMD) heat flow algorithm for the computation of the heat conductivity in one-dimensional lattices. For the well-known Fermi-Pasta-Ulam (FPU) model, it is shown that when the heat field strength is greater than a certain critical value (which depends on the system size) solitons can be generated in molecular dynamics simulations starting from random initial conditions. Such solitons are stable and travel with supersonic speeds. For smaller heat fields, no solitons are generated in the molecular dynamics simulations; the heat conductivity obtained via the NEMD algorithm increases monotonically with the size of the system.

13 pages + 9 figures

Cited by in corpus (12)

• FPU $\beta$ model: Boundary Jumps, Fourier's Law and Scaling
• Finite thermal conductivity in 1D models having zero Lyapunov exponents
• Energy carriers in the Fermi-Pasta-Ulam $\beta$ lattice: Solitons or Phonons?
• Renormalized waves and discrete breathers in $\beta$-FPU chains
• Energy Relaxation in Nonlinear One-Dimensional Lattices
• Thermal conductivities of one-dimensional anharmonic/nonlinear lattices: renormalized phonons and effective phonon theory
• Breathers and Thermal Relaxation in Fermi-Pasta-Ulam Arrays
• Temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with soft on-site potential
• Solitons as candidates for energy carriers in Fermi-Pasta-Ulam lattices
• Thermal conductivity of one-dimensional lattices with self-consistent heat baths: a heuristic derivation
• Propagation dynamics on the Fermi-Pasta-Ulam lattices
• Chemical potential and variable number of particles control the quantum state: Quantum oscillators as a showcase