Thermalization in asymmetric harmonic chains

Abstract

The symmetry of the interparticle interaction potential (IIP) plays a critical role in determining the thermodynamic and transport properties of solids. This study investigates the isolated effect of IIP asymmetry on thermalization. Asymmetry and nonlinearity are typically intertwined. To isolate the effect of asymmetry, we introduce a one-dimensional asymmetric harmonic (AH) model whose IIP possesses asymmetry but no nonlinearity, evidenced by energy-independent vibrational frequencies. Extensive numerical simulations confirm a power-law relationship between thermalization time ($T_{\rm eq}$) and perturbation strength for the AH chain, revealing an exponent larger than the previously observed inverse-square law in the thermodynamic limit. Upon adding symmetric quartic nonlinearity into the AH model, we systematically study thermalization under combined asymmetry and nonlinearity. Matthiessen's rule provides a good estimate of $T_{\rm eq}$ in this case. Our results demonstrate that asymmetry plays a distinct role in enhancing higher-order effects and governing relaxation dynamics.