A molecular dynamics simulation of thermalization of crystalline lattice with harmonic interaction

This paper utilizes molecular dynamics simulations of a harmonic lattice to reveal the complex thermalization process, highlighting distinct relaxation rates for velocity components, power-law frequency proliferation, concurrent defect generation, and two-stage out-of-plane fluctuation behaviors driven by broken symmetry.

The Gist

Imagine a giant, perfectly round trampoline made of a grid of tiny, bouncy springs connecting thousands of marbles. This is the "drum" in our story. In a perfect, frozen world, these marbles would sit still in their neat triangular patterns. But in our real world, things are never perfectly still.

This paper is like a high-speed movie camera watching what happens when we give this trampoline a sudden, random shake. The scientists wanted to understand how a chaotic, jiggly mess eventually settles down into a calm, "thermal" state (what we call heat).

Here is the story of what they found, explained simply:

1. The Two-Step Dance of Relaxation

When you shake the drum, the marbles start moving wildly. The scientists noticed something funny about how they calm down:

• The Flat Dance (Longitudinal): The marbles moving along the surface of the drum (side-to-side) calm down very quickly. It's like a dancer who stops spinning almost immediately.
• The Bouncing Dance (Transverse): The marbles moving up and down (bouncing off the trampoline) take much longer to settle. They keep bouncing around like a rubber ball that just won't stop.

The Analogy: Imagine a crowd of people in a room. If you tell them to stop walking side-to-side, they stop almost instantly. But if you tell them to stop jumping up and down, it takes much longer because the floor is bouncy. The "heat" (temperature) of the system is determined by how hard you shook it at the start, and surprisingly, the temperature goes up with the square of your shake strength.

2. The Frequency Explosion (The Music of Chaos)

As the marbles jiggle, they create vibrations, like notes on a guitar string.

• The Beginning: At first, there is only one main "note" or frequency of vibration.
• The Explosion: As time goes on, the geometry of the triangular grid causes these notes to mix. A low note and a high note crash together to create a new note. Then those new notes mix again.
• The Pattern: The number of different notes (frequencies) doesn't just grow; it explodes in a specific mathematical pattern (a "power law"). It's like a snowball rolling down a hill, picking up more snow (new frequencies) at an accelerating rate.

The Analogy: Think of a single drop of ink falling into a glass of water. At first, it's just one spot. Then it swirls and creates streaks, then clouds, then a complex, swirling galaxy of colors. The paper found that the "swirling" of these vibrations follows a strict, predictable rule, even though it looks chaotic.

3. The "Glitch" in the Grid (Topological Defects)

The triangular grid is very orderly. But as the shaking gets stronger, the grid starts to get "glitchy."

• The Defects: In a perfect triangle, every marble has exactly 6 neighbors. When things get too chaotic, some marbles end up with 5 neighbors (like a pentagon) and some with 7 (like a heptagon). These are called "defects."
• The Connection: The scientists found a magical link: The moment the "frequency explosion" (the music getting crazy) happens, the "defects" (the glitches in the grid) also suddenly appear in huge numbers. It's as if the music and the structural damage are two sides of the same coin.

4. The Two-Stage Bounce (Symmetry Breaking)

Finally, they looked at how the drum bounces up and down.

• Stage 1 (Gentle): When the shake is weak, the drum bounces up and down fairly evenly. It's symmetric.
• Stage 2 (Strong): When the shake gets stronger, something interesting happens. The drum starts to favor bouncing more in one direction (say, up) than the other (down).
• The Break: This is called "breaking the up-down symmetry." The paper shows that this happens at a very specific, low threshold of shaking. Once this symmetry breaks, the way the drum fluctuates changes its mathematical rule.

The Analogy: Imagine a seesaw. At first, it balances perfectly. But if you push it just a little bit harder than a certain point, it stops balancing and starts tilting permanently to one side, even though you aren't pushing it that way anymore.

Why Does This Matter?

This might sound like just a math game with marbles, but it's actually a window into how the universe works.

• From Order to Chaos: It helps us understand how things go from being still to being hot and chaotic.
• Materials Science: It explains why materials might crack or melt when they get hot.
• The Hidden Rules: Even in a system that looks completely random and chaotic, there are hidden, beautiful mathematical rules (like the power laws) that govern the behavior.

In short, the paper tells us that even when you shake a crystal until it's a mess, it doesn't just turn into random noise. It follows a complex, rhythmic dance of frequencies and structural changes that we can now predict and understand.

Technical Summary

1. Problem Statement

The paper addresses the fundamental scientific question of how many-body systems achieve thermal equilibrium through the thermalization process. While macroscopic thermodynamics describes the final state, the microscopic dynamical mechanisms—bridging classical dynamics and thermodynamics—remain complex. Specifically, the study investigates:

• How a crystalline lattice loses memory of its initial state and explores permissible states.
• The role of geometric nonlinearity in a system with purely harmonic (linear) interaction potentials.
• The emergence of statistical regularities, topological defects, and symmetry breaking during the transition from a disturbed state to thermal equilibrium.

2. Methodology

The authors employ Molecular Dynamics (MD) simulations based on a specific model and numerical approach:

• Model System: A circular, drum-like configuration of a 2D triangular lattice.
  • Particles: $N$ point particles connected by identical linear springs (harmonic potential $V_0(r) = \frac{1}{2}k_0(r-\ell_0)^2$).
  • Boundary Conditions: Anchored (clamped) boundary; particles in an outer annulus are fixed in space.
  • Dimensionality: Particles move in 3D space (allowing out-of-plane deformations), though 2D-constrained cases are also analyzed for comparison.
• Initial Conditions: The lattice starts in mechanical equilibrium. Thermalization is initiated by imposing a random velocity vector $\vec{v}_{ini} = v_0 \vec{\xi}$ on each particle, where $v_0$ is the disturbance strength and $\vec{\xi}$ is a random vector.
• Numerical Integration: The equations of motion are solved using the Verlet integration algorithm with a time step $dt = 10^{-3}\tau_0$. The system is treated as a dissipationless Hamiltonian system, ensuring energy conservation.
• Analysis Techniques:
  • Velocity Space: Analysis of longitudinal ($v_{\parallel}$, in-plane) and transverse ($v_{\perp}$, out-of-plane) velocity distributions against the Maxwell-Boltzmann distribution.
  • Spectral Analysis: Fourier transformation of the kinetic energy curve to track the proliferation of dominant frequencies over time.
  • Coordinate Space: Delaunay triangulation to identify topological defects (disclinations) and analysis of out-of-plane deformation fluctuations.
  • Theoretical Modeling: A discrete theoretical model was constructed to simulate frequency mixing and proliferation.

3. Key Contributions and Results

A. Anisotropic Velocity Relaxation

• Distinct Relaxation Rates: The study reveals that the longitudinal component of velocity ($v_{\parallel}$) relaxes significantly faster than the transverse component ($v_{\perp}$).
• Mechanism: In-plane motion leads to higher efficiency in velocity relaxation compared to out-of-plane motion.
• Temperature Scaling: Both components eventually converge to a Maxwell distribution characterized by a common temperature $T$. The temperature scales with the square of the initial disturbance strength: $T \propto v_0^2$.

B. Nonlinear Proliferation of Frequencies

• Power Law Growth: The number of dominant frequencies ($\Omega(t)$) in the system grows according to a power law with respect to observation time $t$ (after an initial transient period):
  $$\Omega(t) - \Omega(t^*) = C(t - t^*)^\alpha$$
• Tunable Exponent: The exponent $\alpha$ is not constant; it depends on the disturbance strength $v_0$.
  • $\alpha$ ranges from $\approx 1.1$ to $2.9$ as $v_0$ varies.
  • A critical disturbance strength ($v_0^* \approx 0.6$) triggers a rapid rise in $\alpha$, indicating a transition to more efficient frequency proliferation.
• Geometric Origin: The frequency mixing (sum/difference of frequencies) and drift are driven by the intrinsic geometric nonlinearity of the triangular lattice configuration, even though the interaction potential is harmonic.
• Theoretical Model: A discrete model incorporating frequency mixing, doubling, and halving rules successfully reproduces the observed power-law behavior and non-unity exponents.

C. Concurrent Proliferation of Topological Defects

• Correlation: There is a direct correlation between the rapid proliferation of frequencies and the emergence of topological defects (disclinations).
• Critical Threshold: Both phenomena occur concurrently above a critical disturbance strength ($v_0 \approx 0.5 - 0.6$).
• Defect Types: The study identifies the formation of dislocations (pairs of 5- and 7-fold disclinations) and quadrupoles. The unbinding of these defects leads to the destruction of crystalline order.
• Finite Size Stability: Unlike infinite 2D systems where the Hohenberg-Mermin-Wagner theorem predicts the absence of long-range order, the finite-size drum system maintains crystalline order at low $v_0$ even over long simulation times ($10^4 \tau_0$).

D. Two-Stage Out-of-Plane Fluctuations and Symmetry Breaking

• Two-Stage Behavior: The mean squared transverse displacement $\langle h^2 \rangle$ exhibits two distinct regimes governed by different power laws as $v_0$ increases.
• Scaling Relation: $\langle h^2 \rangle \propto T^\beta$, where the exponent $\beta$ changes across a transition point $v_{0,b} \approx 0.023$.
  • Regime 1 ($v_0 < v_{0,b}$): $\beta \approx 0.65$ ($k_1/2$).
  • Regime 2 ($v_0 > v_{0,b}$): $\beta \approx 0.49$ ($k_2/2$).
• Symmetry Breaking: This transition coincides with the breaking of up-down symmetry in the out-of-plane deformations. While the system is symmetric at very low $v_0$, a measurable asymmetry ($\langle h_+ \rangle - \langle h_- \rangle \neq 0$) emerges at $v_{0,b}$, triggering the change in fluctuation susceptibility.

4. Significance and Implications

• Harmonic vs. Nonlinear Dynamics: The work demonstrates that geometric nonlinearity alone (arising from the lattice configuration) is sufficient to induce complex nonlinear dynamics, frequency mixing, and thermalization, even in the absence of anharmonic interaction potentials.
• Thermalization Mechanism: It provides a microscopic view of how deterministic chaos and frequency proliferation lead to statistical equilibrium, bridging the gap between classical dynamics and thermodynamics.
• Material Stability: The findings offer insights into the mechanical stability of crystalline materials in thermal environments and the onset of 2D melting (topological phase transitions).
• Symmetry Breaking: The identification of a two-stage fluctuation behavior linked to up-down symmetry breaking in a clamped membrane system advances the understanding of fluctuation dynamics in constrained elastic systems.

In conclusion, this paper elucidates the rich dynamical adaptations of many-body systems to external disturbances, showing that thermalization is a multi-stage process involving distinct relaxation rates, power-law frequency proliferation, and symmetry-breaking fluctuations.