Intermediate Thermal Equilibrium Stages in Molecular Dynamics Simulations of two Bodies in Contact

This study utilizes classical molecular dynamics simulations of argon-based two- and three-region models to analyze the intermediate fluctuations, correlations, and temperature distributions that characterize the process of heat conduction leading to thermal equilibrium, thereby providing a detailed microscopic perspective on the Zeroth Law of Thermodynamics.

The Gist

Imagine you have two cups of coffee: one is scalding hot, and the other is ice cold. If you put them next to each other, eventually, they will both settle at a lukewarm temperature. This is the basic idea of the Zeroth Law of Thermodynamics: things in contact will eventually reach the same temperature.

But here is the question this paper asks: What happens while they are getting there? Does the heat flow smoothly like water in a pipe, or is it a chaotic, bumpy ride?

The authors of this paper used a super-powerful computer simulation (like a microscopic movie camera) to watch this process happen in real-time, atom by atom. They didn't just look at the start and the finish; they zoomed in on the messy middle part.

Here is the story of their experiment, explained simply:

The Setup: The "Thermal Dance Floor"

Imagine a long hallway divided into rooms.

• Scenario A (The Simple Dance): They put a hot crowd of people (Argon atoms) on the left and a cold crowd on the right. Between them, they built three thin, transparent walls made of "graphene" (a super-thin material). These walls let heat pass through but keep the people in their own rooms.
• Scenario B (The Complicated Dance): They added a tiny, cramped middle room between the hot and cold sides. This middle room had fewer people and started at a medium temperature.

They turned off the air conditioning (the external temperature control) and let the crowds interact naturally, watching how the heat traveled from the hot side to the cold side.

The Findings: It's Not a Straight Line

1. The Simple Case (Two Rooms)

In the first scenario, the heat transfer was like a smooth slide. The hot side cooled down, the cold side warmed up, and they met in the middle. It was fast, predictable, and the temperature settled down quickly. It was a "textbook" example of reaching equilibrium.

2. The Complicated Case (Three Rooms)

When they added that tiny middle room, things got weird and interesting.

• The "Bottleneck" Effect: The middle room acted like a busy traffic intersection. The heat had to flow from the left, stop at the middle, and then go to the right. This slowed everything down.
• The "Bumpy Ride": Instead of a smooth slide, the temperature in the middle room went up and down like a rollercoaster. It was chaotic. The authors found that the middle room kept "jittering" with energy long after the sides had settled down.
• The "Double Peak" Mystery: If you looked at the temperature history of the middle room, it didn't just go straight to the answer. It seemed to get stuck in a "halfway" state for a while, creating a temporary, local balance before finally giving in to the global balance.

The Analogy: The Relay Race vs. The Solo Run

Think of the two scenarios like a relay race:

• Scenario A is a solo runner sprinting from start to finish. They run fast and stop when they cross the line.
• Scenario B is a relay race with a tiny, nervous middle runner. The first runner passes the baton to the middle runner, who fumbles it a bit, runs in circles, and then passes it to the final runner. The whole race takes longer, and the middle runner is sweating and shaking (fluctuating) much more than the others.

Why Does This Matter?

Usually, we think of the Zeroth Law as a simple rule: "If A touches B, and B touches C, then A and C are the same temperature." We assume this happens instantly and smoothly.

This paper shows that reality is messier.

• Heat isn't just a flow; it's a conversation. The atoms are constantly bumping into each other, exchanging energy in bursts.
• Size matters. The tiny middle room was more "jittery" because it had fewer atoms to share the load. It's like trying to balance a seesaw with one person versus a whole team; the single person wobbles much more.
• The path to peace is bumpy. Before a system becomes perfectly stable, it often gets stuck in "local" states where things look balanced for a moment, but aren't quite there yet.

The Bottom Line

The authors used advanced math (like a tool called GARCH, which is usually used to predict stock market volatility) to measure how "nervous" the temperature was. They found that adding a middle layer makes the system take 2.5 times longer to calm down.

In everyday terms: If you want your house to reach a comfortable temperature, don't just think about the thermostat. Think about the walls, the furniture, and the air pockets in between. The journey to a comfortable temperature is full of little bumps, delays, and temporary pauses, not just a straight line. The Zeroth Law tells us where we end up, but this paper shows us the chaotic, fascinating dance we have to do to get there.

Technical Summary

1. Problem Statement

The study addresses a gap in the classical formulation of the Zeroth Law of Thermodynamics. While the law states that if two systems are in thermal equilibrium with a third, they are in equilibrium with each other, it typically describes only the final state. It rarely addresses the dynamic, non-equilibrium process leading to that state, particularly regarding:

• The temporal evolution of temperature fluctuations.
• The role of intermediate regions (metastable states) in heat transport.
• The statistical nature of thermal relaxation before global equilibrium is achieved.

The authors aim to investigate the microscopic mechanisms of heat conduction and the "intermediate stages" of equilibration in ideal gas systems separated by heat-conducting walls.

2. Methodology

The researchers employed Classical Molecular Dynamics (MD) simulations to model the behavior of Argon atoms.

• Simulation Setup:

  • Software: LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator).
  • System: A 3D box ($40 \times 40 \times 40$ Å) containing Argon atoms modeled via the Lennard-Jones (LJ) potential.
  • Walls: Separating walls were simulated as fixed graphene sheets (Carbon atoms) interacting via LJ potentials.
    • Internal walls: Designed as diathermal (allowing heat transfer) by reducing interaction energy parameters.
    • External walls: Designed as adiabatic (insulating) to isolate the system from the environment.
  • Ensembles: Initial equilibration in the NVT ensemble (constant Number, Volume, Temperature), followed by the main simulation in the NVE ensemble (constant Number, Volume, Energy) to observe spontaneous thermal transfer.

• Experimental Configurations:

  • Case 1 (Two-Region): Two regions (Left/Right) with 400 atoms each, separated by three thin walls. Initial temperatures: 100 K (Left) and 500 K (Right).
  • Case 2 (Three-Region): Three regions (Left/Middle/Right). Lateral regions have 400 atoms; the central region has 100 atoms (1/8th the volume). Initial temperatures: 100 K (Left), 300 K (Middle), and 500 K (Right).

• Analytical Tools:

  • GARCH Model: Used to analyze conditional temperature volatility (time-dependent variance) to distinguish between different statistical regimes of thermal fluctuation.
  • Correlation Functions: Calculated to determine thermal relaxation times ($\tau$) based on the decay of temperature fluctuations ($C(t) \sim e^{-t/\tau}$).
  • Frequency Distributions: Analyzed to identify unimodal vs. bimodal temperature states.

3. Key Results

A. Temperature Evolution and Equilibration Speed

• Case 1 (Two Regions): Temperatures converged almost exponentially to an average value (~300 K) with minimal fluctuations. Equilibrium was reached relatively quickly.
• Case 2 (Three Regions): The presence of the intermediate region introduced complex behavior. The central region exhibited significant temperature variations and acted as a heat transport mediator, delaying the overall equilibration process.
• Relaxation Time:
  • Case 1 average relaxation time: 662 time steps.
  • Case 2 lateral regions: 1669 time steps (approx. 2.5 times slower than Case 1).
  • Case 2 central region: 216 time steps (fastest to stabilize locally, but creates a bottleneck for the whole system).

B. Volatility and Fluctuations (GARCH Analysis)

• Volatility Clustering: In Case 2, the central region showed "volatility clusters" (intermittent high-variance events) at specific time intervals (e.g., 4000–4500 and 6000–8000 steps), indicating non-uniform energy exchange.
• Susceptibility: The smaller central region (fewer atoms) was more susceptible to energy fluctuations, maintaining high volatility even after the lateral regions stabilized.
• Convergence: Lateral regions in Case 2 reached low volatility intensity later than in Case 1.

C. Distribution and Transitivity

• Bimodality: In Case 2, the temperature frequency distribution for lateral regions showed bimodal peaks (a main peak and a secondary peak) before reaching the final global equilibrium.
• Metastability: This bimodality suggests the existence of transient "quasi-equilibrium" states (local maxima in the free energy landscape) where thermal synchronization occurs at different time scales. This challenges the instantaneous transitivity implied by a static reading of the Zeroth Law.

D. Correlations

• Anti-correlation: In Case 2, the lateral regions exhibited strong anti-correlation (heating of one coincides with cooling of the other).
• Mediator Role: The central region showed slight correlation with the cold zone and anti-correlation with the hot zone, confirming its role as a thermal intermediary that buffers the direct exchange between the ends.

4. Key Contributions

1. Dynamic View of the Zeroth Law: The paper moves beyond the static definition of thermal equilibrium to characterize the kinetic pathway and intermediate stages required to satisfy the law.
2. Quantification of Relaxation: It provides quantitative metrics (relaxation times, volatility clusters) for how geometric configurations (specifically the addition of an intermediate region) impact heat transport efficiency.
3. Statistical Mechanics Application: The novel application of the GARCH model to MD data allows for the rigorous characterization of conditional temperature volatility, distinguishing between stable and intermittent thermal regimes.
4. Identification of Metastable States: The study identifies and characterizes "local maxima" in the equilibration process, demonstrating that systems can temporarily stabilize in inhomogeneous thermal states before achieving global homogeneity.

5. Significance

This work is significant for non-equilibrium thermodynamics and microscopic heat transport theory. It demonstrates that the path to thermal equilibrium is not a simple, monotonic decay but a complex process involving:

• Bottlenecks: Intermediate regions can significantly delay equilibration despite having an initial temperature close to the final equilibrium.
• Fluctuation-Driven Dynamics: The time to reach equilibrium is governed by the persistence of thermal fluctuations and the statistical properties of the subsystems (e.g., size-dependent volatility).
• Practical Implications: These findings are relevant for designing nanoscale thermal management systems, understanding energy dissipation in cellular biophysics, and modeling glass transitions where systems get trapped in local energy minima.

In conclusion, the authors argue that the Zeroth Law describes the destination of a system, but the journey is governed by dynamic fluctuations, geometric constraints, and transient metastable states that require detailed statistical analysis to fully understand.