Large temperature-up-jump simulations of a binary Lennard-Jones system

This paper investigates the validity of the Tool-Narayanaswamy material-time concept for describing the physical aging of a binary Lennard-Jones liquid following large temperature up-jumps, finding that the theory works well for smaller jumps but fails for larger ones, thereby confirming its limitation to systems near equilibrium while suggesting future research into quantity-specific or locally defined material times.

The Gist

The Big Picture: The "Glassy" Problem

Imagine you have a jar of honey. If you leave it on the shelf, it slowly hardens and settles into a rigid state. This is similar to what happens in glasses (like window glass or plastic). They aren't quite solid crystals, but they aren't liquids either. They are "stuck" in a messy, frozen state.

When you change the temperature of this glassy material, it doesn't snap back to a new state instantly. It takes a long, slow time to "relax" or settle down. This slow settling process is called physical aging.

Scientists have a theory called the Tool-Narayanaswamy (TN) theory to predict how this aging happens. The theory suggests that instead of using a normal clock (seconds, minutes, hours), we should use a special "Material Time" clock.

The Analogy: Think of a normal clock as a metronome ticking at a steady beat. But when you age a glass, the "beat" of the material speeds up or slows down depending on how tired or energetic it is. The "Material Time" is like a custom stopwatch that speeds up when the material is moving fast and slows down when it's sluggish, so that the aging process looks the same on this special clock, no matter what the temperature is.

What Did This Paper Do?

The researchers wanted to test if this "Material Time" idea works when the temperature change is huge.

Usually, scientists test this with small temperature changes (like moving a glass from a cool room to a warm room). But in this paper, they simulated a massive temperature jump.

• The Setup: They took a simulated liquid made of two types of particles (like big red balls and small blue balls) that was very cold and sluggish.
• The Jump: They suddenly heated it up significantly (like taking a frozen block of ice and throwing it into a hot oven).
• The Goal: They wanted to see if the "Material Time" clock could still make sense of the chaos, or if the system was too far from equilibrium for the theory to work.

The "Triangle" Test

Before they could use the Material Time clock, they had to prove it was valid. They used a mathematical rule called the "Triangular Relation."

The Analogy: Imagine you are walking through a foggy forest.

1. You look at a tree at time A.
2. You look at the same tree at time B.
3. You look at it at time C.

The "Triangular Relation" says: If you know how much the tree moved between A and B, and how much it moved between B and C, you can perfectly predict how much it moved between A and C. It's like a triangle where if you know two sides, the third is fixed.

The Result: They checked this for the potential energy of the particles. It worked! The triangle held up. This meant they could confidently define a "Material Time" based on the energy of the system.

The Main Experiment: Does the Clock Work?

Once they had their "Material Time" clock, they watched five different things happen as the system aged:

1. Energy: How much energy the particles had.
2. Movement: How far particles moved (like a drunk person stumbling).
3. Structure: How the particles were arranged.
4. Clumping: How much the particles moved together in groups.
5. Weirdness: How much the movement deviated from a normal pattern.

They ran two scenarios:

• Scenario A: A medium-sized temperature jump.
• Scenario B: A huge temperature jump (from very cold to warm).

The Results:

• For the medium jump: The "Material Time" worked beautifully. When they plotted the data using the special clock, all the messy curves collapsed into a single, neat line. It was like taking a tangled ball of yarn and smoothing it out perfectly.
• For the huge jump: The "Material Time" failed. The curves didn't collapse into one line. They were still messy and spread out.

Why Did It Fail?

The paper concludes that the "Material Time" theory works best when the system is never too far from equilibrium (not too far from being "settled").

The Metaphor:
Imagine a classroom.

• Small Jump: The teacher asks the class to switch from "quiet reading" to "group discussion." Everyone adjusts their behavior smoothly. You can predict the noise level easily.
• Huge Jump: The teacher suddenly screams "FREEZE!" and then immediately yells "RUN AROUND THE ROOM!" The students are confused, some run, some hide, some laugh. The behavior is chaotic and different for every student. You can't use a single "classroom clock" to predict what everyone is doing because the system is too chaotic.

In the simulation, the huge temperature jump created Dynamic Heterogeneity. This means some parts of the material were moving fast while others were still stuck in the "slow" mode. Because different parts of the material were aging at different rates, a single "Global Clock" (Material Time) couldn't describe the whole system anymore.

The Takeaway

1. The Theory is Robust (mostly): The "Material Time" concept is a powerful tool for understanding how glasses age, but it has limits.
2. Extreme Conditions Break the Rules: If you shock the system too hard (a massive temperature jump), the simple "one clock fits all" idea breaks down.
3. Future Questions: The authors ask: Maybe we need local clocks? Instead of one clock for the whole glass, maybe we need a different clock for the fast-moving parts and a different clock for the slow-moving parts?

In short: The "Material Time" is a great GPS for navigating the aging of glass, but if you drive off the road into a massive storm (a huge temperature jump), the GPS loses its signal, and you need a more complex map to find your way.

Technical Summary

1. Problem Statement

Physical aging in glasses is a complex, nonlinear phenomenon where a system slowly relaxes toward equilibrium after a temperature change. The Tool-Narayanaswamy (TN) formalism is a widely used theoretical framework that attempts to linearize this process by introducing a material time ($\xi$). According to TN theory, if laboratory time ($t$) is replaced by material time, all aging trajectories for a given final temperature should collapse onto a single master curve, regardless of the initial state.

While the TN formalism works well for small temperature perturbations (near equilibrium), its validity is challenged by large temperature up-jumps. In such cases, the system starts from a highly non-equilibrium state with a vastly different relaxation time than the final state. The central question of this paper is: Does the single-parameter TN material-time concept hold for extreme, large temperature up-jumps in a binary Lennard-Jones system? Specifically, can a single global material time defined by potential energy accurately describe the aging dynamics of diverse structural and dynamic quantities?

2. Methodology

The authors employed large-scale molecular dynamics (MD) simulations using the GPU-optimized RUMD code.

• System: A modified binary Kob-Andersen Lennard-Jones (mBLJ) mixture consisting of 80% large particles (A) and 20% small particles (B). The model parameters were slightly adjusted (shifted-force cutoffs) to suppress crystallization, allowing for long simulations of the glassy state.
• Protocol:
  • Systems were equilibrated at low temperatures ($T = 0.43$ and $T = 0.37$).
  • Large Temperature Up-Jumps: The temperature was instantaneously raised to a final temperature of $T = 0.48$.
  • Monitoring: Five distinct quantities were tracked during the aging process:
    1. Potential energy time-autocorrelation ($C_{uu}$).
    2. Self-intermediate scattering function ($F_s$).
    3. Mean-square displacement ($\langle \Delta r^2 \rangle$).
    4. Dynamic susceptibility ($\chi_4$), measuring dynamic heterogeneity.
    5. Non-Gaussian parameter ($\alpha_2$), measuring deviations from Gaussian diffusion.
• Material Time Definition: The material time $\xi(t)$ was derived strictly from the potential-energy time-autocorrelation function. The definition relies on the triangular relation (proposed by Cugliandolo and Kurchan), which states that for any three times $t_1 < t_2 < t_3$, the correlation $C_{13}$ is a unique function of $C_{12}$ and $C_{23}$. If this relation holds, a unique material time exists.

3. Key Contributions and Results

A. Validation of the Triangular Relation

The authors first tested whether the triangular relation holds for the potential energy during large up-jumps.

• Finding: The triangular relation was found to be well-obeyed (with only minute deviations) for the potential energy, even for the largest jump ($0.37 \to 0.48$).
• Implication: This validates the mathematical consistency required to define a global material time $\xi$ based on potential energy. The aging process, despite being highly nonlinear, effectively takes place "close to equilibrium" in terms of the potential energy landscape.

B. Testing the Collapse of Aging Quantities

Using the material time $\xi$ derived from potential energy, the authors re-plotted the time-autocorrelation functions of all five monitored quantities against the material-time difference ($\xi_2 - \xi_1$).

• Small Jump ($0.43 \to 0.48$):
  • The collapse was successful for most quantities (potential energy, $F_s$, and mean-square displacement).
  • The dynamic susceptibility ($\chi_4$) and non-Gaussian parameter ($\alpha_2$) showed poorer collapse but were significantly improved compared to laboratory time plots.
• Large Jump ($0.37 \to 0.48$):
  • The collapse was significantly worse. While the introduction of material time eliminated much of the "spread" seen in laboratory time, a unique master curve was not achieved for any of the five quantities.
  • Specific Anomalies:
    • The potential energy autocorrelation showed residual slope differences, suggesting the global clock is breaking down.
    • Divergent Behavior: $\chi_4$ and $\alpha_2$ behaved oppositely. $\chi_4$ curves shifted toward later times (slower relaxation) compared to equilibrium, while $\alpha_2$ curves at short waiting times relaxed faster than those at longer waiting times. This suggests that the onset of flow events ($\alpha_2$) and the spatial correlation of those events ($\chi_4$) evolve differently under extreme non-equilibrium conditions.

4. Significance and Discussion

The study provides critical insights into the limits of the TN formalism:

1. Limits of the Global Clock: The results confirm that the TN formalism works best for systems that are never very far from equilibrium. For large up-jumps, a single global material time defined by potential energy is insufficient to describe the full complexity of the aging dynamics.
2. Role of Dynamic Heterogeneity: The failure of the collapse for the large jump suggests that dynamic heterogeneity plays a dominant role. In extreme up-jumps, different regions of the glass may age at different rates (some regions may remain "frozen" while others become active), invalidating the assumption of a single global clock.
3. Future Directions: The paper raises two pivotal questions for future research:
   • Does the collapse improve if each physical quantity is allowed its own material time (multi-parameter aging)?
   • Can the formalism be saved by generalizing material time to be locally defined (a "local clock") to account for spatial dynamic heterogeneity, rather than assuming a global time scale?

Conclusion

This paper demonstrates that while the triangular relation holds for potential energy in large temperature up-jumps, the resulting single-parameter material time fails to universally collapse the aging dynamics of other structural and dynamic observables. This highlights the breakdown of the standard TN formalism under extreme non-equilibrium conditions and points toward the necessity of incorporating local dynamic heterogeneity or multi-parameter descriptions to accurately model glass aging in such regimes.