Energy transfer mechanisms in a dipole chain: From energy equipartition to the formation of breathers

Equilibria and dynamics of two coupled chains of interacting dipoles

We explore the energy transfer dynamics in an array of two chains of identical rigid interacting dipoles. Varying the distance b between the two chains of the array, a crossover between two different ground-state (GS) equilibrium configurations is observed. Linearizing around the GS configurations, we verify that interactions up to third nearest neighbors should be accounted to accurately describe the resulting dynamics. Starting with one of the GS, we excite the system by supplying it with an excess energy ΔK located initially on one of the dipoles. We study the time evolution of the array for different values of the system parameters b and ΔK. Our focus is hereby on two features of the energy propagation: the redistribution of the excess energy ΔK among the two chains and the energy localization along each chain. For typical parameter values, the array of dipoles reaches both the equipartition between the chains and the thermal equilibrium from the early stages of the time evolution. Nevertheless, there is a region in parameter space (b,ΔK) where even up to the long computation time of this study, the array does neither reach energy equipartition nor thermalization between chains. This fact is due to the existence of persistent chaotic breathers.

Chaos and thermalization in a classical chain of dipoles

We explore the connection between chaos, thermalization, and ergodicity in a linear chain of N interacting dipoles. Starting from the ground state, and considering chains of different numbers of dipoles, we introduce single site excitations with excess energy ΔK. The time evolution of the chaoticity of the system and the energy localization along the chain is analyzed by computing, up to a very long time, the statistical average of the finite-time Lyapunov exponent λ(t) and the participation ratio Π(t). For small ΔK, the evolution of λ(t) and Π(t) indicates that the system becomes chaotic at approximately the same time as Π(t) reaches a steady state. For the largest considered values of ΔK the system becomes chaotic at an extremely early stage in comparison with the energy relaxation times. We find that this fact is due to the presence of chaotic breathers that keep the system far from equipartition and ergodicity. Finally, we show numerically and analytically that the asymptotic values attained by the participation ratio Π(t) fairly correspond to thermal equilibrium.