Energy conservation and pressure relaxation in an extended two-temperature model for copper with an electron temperature-dependent interaction potential

This paper presents an implementation of an electron temperature-dependent interaction potential for copper within a two-temperature model molecular dynamics framework, introducing a specific algorithm to ensure energy conservation and addressing pressure relaxation effects caused by electron temperature gradients following laser irradiation.

The Gist

Imagine you are a conductor leading a massive orchestra of copper atoms. Usually, this orchestra plays a steady, predictable tune. But what happens when you blast them with a super-fast, intense laser pulse?

This paper is about how we simulate that chaotic moment on a computer, and how the authors fixed a few major "bugs" in the way we've been doing it until now.

Here is the story of the paper, broken down into simple concepts:

1. The Setup: The Two-Temperature Orchestra

When a laser hits copper, it doesn't heat everything up evenly.

• The Electrons (The Violins): They get hit first. They get super excited, hot, and energetic almost instantly.
• The Atoms (The Cellos): They are heavy and slow. They don't know the laser is there yet. They stay cool for a split second.

This creates a "Two-Temperature" situation: The violins are screaming (hot electrons), but the cellos are still playing a slow, calm melody (cool atoms). Eventually, the violins pass their energy to the cellos, and the whole orchestra heats up.

2. The Problem: The Rules of the Game Keep Changing

In the past, computer simulations treated the copper atoms like they were made of a rigid, unchanging material. The "interaction potential" (the invisible glue holding the atoms together) was fixed.

The Reality: When the electrons get super hot, they actually change the glue!

• Bond Hardening: Imagine the copper atoms are holding hands. When the electrons get excited, they actually squeeze their hands tighter. The material gets harder and more resistant to breaking.
• The Reference Shift: The "zero point" of energy changes. It's like if you suddenly decided that sea level is now 100 feet higher. If your computer simulation doesn't update its map, it thinks the water is rising when it's actually just the map that changed.

The Bug: Previous simulations didn't account for these changes. They kept using the old "rigid glue" rules while the electrons were screaming. This led to errors in how much energy was conserved and how the material reacted.

3. The Solution: A Smarter Energy Accountant

The authors (Simon and Johannes) built a new, smarter system for the computer simulation.

The Energy Conservation Fix:
Think of the simulation as a bank account.

• Old Way: The laser deposits money (energy). The electrons spend some, the atoms spend some. But because the "glue" changed strength, the bank thought money appeared out of nowhere or disappeared. The books didn't balance.
• New Way: The authors created a strict accounting algorithm. Every time the electrons get hotter, the computer recalculates the "cost" of holding the atoms together. If the glue gets stronger (bond hardening), the system knows that some of the laser energy is being "spent" just to tighten that glue, rather than heating the atoms up. This ensures the total energy in the system never magically changes.

The Pressure Fix (The "Blast Force"):
When the electrons get hot, they push outward, creating pressure.

• Old Way: Scientists used to manually add a "kick" or a "blast force" to the atoms to simulate this push. It was like a director shouting, "Hey, push harder!" It worked, but it was a clumsy, artificial fix that sometimes broke the energy rules.
• New Way: The authors showed that if you use the new "temperature-dependent glue," the pressure happens naturally. The atoms push each other apart because the rules of their interaction changed. You don't need to shout "push!"; the physics does it for you. It's like realizing the orchestra naturally gets louder when the violins get excited, rather than having to manually turn up the volume knob.

4. The Results: What Happens When You Get It Right?

When they ran the simulation with these new, smarter rules on a copper target:

• Slower Melting: Because the "glue" gets harder when hot (bond hardening), the copper resists melting longer. It's like the material puts up a better fight.
• Less Ejection: In the old simulations, the artificial "blast force" would often rip individual atoms off the surface immediately. In the new simulation, the material holds together better. The atoms don't fly off as easily; they melt more slowly and deeply.
• Cooler Electrons: Because some energy is used to tighten the bonds, the electrons don't get quite as hot as the old models predicted.

The Big Picture Analogy

Imagine you are trying to predict how a balloon pops when you poke it with a hot needle.

• The Old Model: You assume the rubber stays the same thickness and strength. You calculate that the needle will pop it instantly.
• The New Model: You realize that as the needle heats the rubber, the rubber actually gets stiffer and tougher right where the heat is. So, the needle has to push harder and longer before it pops. The balloon doesn't explode as fast, and the hole is different.

Why Does This Matter?

This isn't just about copper; it's about precision. If we want to use lasers to manufacture microchips, clean medical tools, or cut metal with extreme precision, we need to know exactly how the material will react.

By fixing the "accounting" (energy conservation) and letting the physics happen naturally (pressure relaxation) instead of forcing it, the authors have given scientists a much more accurate map for navigating the chaotic world of ultra-fast laser interactions. They showed that when you treat the material with respect for its changing nature, the simulation tells a truer story.

Technical Summary

Here is a detailed technical summary of the paper "Energy conservation and pressure relaxation in an extended two-temperature model for copper with an electron temperature-dependent interaction potential" by Simon Kümmel and Johannes Roth.

1. Problem Statement

When materials like copper are subjected to ultra-short laser pulses, they undergo strong electronic excitation, creating a non-equilibrium state where electrons and the lattice (atoms) have different temperatures. The Two-Temperature Model (TTM) coupled with Molecular Dynamics (MD) is the standard framework for simulating this process.

However, existing TTM-MD simulations often rely on static interatomic potentials (IPs) that do not account for changes in electronic excitation. Recent Density Functional Theory (DFT) calculations have shown that interatomic forces, bond strengths, and reference energies change significantly with electron temperature ($T_e$).

• The Core Issue: When using an electron temperature-dependent potential, standard energy conservation schemes fail because the potential energy landscape shifts as $T_e$ changes.
• The Pressure Issue: Electron temperature gradients generate pressure gradients. Previous studies often artificially added a "blast force" term to account for this, but this approach often violates energy conservation and introduces arbitrary assumptions about density dependence.

The paper addresses the need for a rigorous implementation of electron temperature-dependent potentials that strictly enforces energy conservation and naturally reproduces pressure relaxation without artificial force terms.

2. Methodology

The authors developed a modified TTM-MD framework for copper incorporating the following key methodological advancements:

A. Electron Temperature-Dependent Interaction Potential

The study utilizes a potential developed in previous work [20] where the interaction strength and reference energy depend on $T_e$.

• Reference Energy Correction: The potential development adopted single-atom energies at specific $T_e$ as reference energies. Since these energies change with excitation, the authors derived a correction scheme to ensure the total energy (lattice + electron) remains consistent.
• Polynomial Fitting:
  • Electron energy ($E_e$) and single-atom energy ($E_{single}$) were fitted to second-order polynomials of $T_e$.
  • The total free energy of the lattice was fitted to a second-order polynomial: $E_{tot} = a_{str} + b(k_B T_e) + c(k_B T_e)^2$.
  • Crucially, the coefficients ($a_{str}, b, c$) were found to be density-dependent and were fitted to sixth-order polynomials of density ($\rho$).

B. Energy Conservation Algorithm

To enforce energy conservation during TTM-MD time steps (where energy is exchanged between electrons and the lattice):

1. Energy Balance: The algorithm calculates the total energy change ($\Delta E$) in a Finite Difference (FD) cell due to laser absorption or heat diffusion.
2. Iterative Correction: It solves a quadratic equation to find the final equilibrium temperature ($T_f$) where the sum of the new electron energy and the new lattice free energy equals the initial total energy plus $\Delta E$.
   • The equation solved is: $A(k_B T_f)^2 + B(k_B T_f) + C = 0$.
   • Special handling (edge cases) is applied for temperatures above 1.2 eV, where the potential behavior changes.

C. Pressure Relaxation and Blast Force Comparison

• Natural Pressure: The authors demonstrate that the electron temperature-dependent IP naturally generates pressure gradients due to $T_e$ gradients, eliminating the need for artificial force terms.
• Blast Force Comparison: They compared their natural approach against the traditional "blast force" method ($\vec{F}_{blast} = -\nabla p_e / n_i$), which is often added as an external force term. The study highlights that the blast force method violates energy conservation and assumes uniform force distribution within an FD cell, whereas the IP method correctly localizes forces at interfaces.

D. Simulation Setup

• System: A copper target (2000 unit cells) irradiated by a 1056 nm laser pulse (680 fs duration, 3.8 J/m² fluence).
• Parameters: Used DFT-derived electron heat capacity, electron-phonon coupling, and optical properties.
• Phase Detection: A Local Order Parameter (LOP) was used to distinguish between solid and liquid states dynamically.
• Comparisons: Four scenarios were simulated:
  1. Static IP (no blast force).
  2. Static IP + Blast Force.
  3. Dynamic IP (no energy conservation).
  4. Dynamic IP + Energy Conservation.

3. Key Contributions

1. Rigorous Energy Conservation Scheme: Presented a novel algorithm that solves for the equilibrium temperature in TTM-MD simulations using electron temperature-dependent potentials, ensuring strict energy conservation even as the potential landscape shifts.
2. Natural Pressure Relaxation: Demonstrated that electron temperature-dependent potentials inherently reproduce the pressure gradients caused by $T_e$ gradients, rendering the artificial "blast force" term unnecessary and physically inconsistent.
3. Density-Dependent Parameterization: Provided a comprehensive set of fit parameters for the potential energy coefficients as a function of density, allowing for accurate simulations across a wide range of material states (solid to ablated plasma).
4. Systematic Comparison: Conducted large-scale simulations comparing the new approach against existing methods, quantifying the differences in material dynamics.

4. Key Results

• Electron Temperature Suppression: Enforcing energy conservation with the dynamic potential results in significantly lower electron temperatures compared to simulations without conservation or with static potentials. This is because part of the absorbed energy is consumed to adjust the changing interaction strength (bond hardening).
• Delayed Melting and Ablation:
  • The lower electron temperature leads to a slower increase in lattice temperature.
  • The dynamic potential predicts bond hardening at high excitation, effectively raising the melting point.
  • Result: The ablation process is significantly delayed, and the total amount of ablated material is reduced compared to static potential models.
• Pressure Wave Dynamics:
  • Blast Force Model: Produces sharp pressure gradients at the surface, leading to the immediate ejection of single atoms or small clusters.
  • Dynamic IP (Conserved): Produces smoother pressure gradients. The pressure wave propagating into the sample is weaker because less energy is available for thermal expansion due to the energy conservation constraint.
• Validity of Blast Force: The study confirms that while the blast force approximates the magnitude of forces, it fails to capture the correct spatial distribution (localizing forces at interfaces) and violates energy conservation.

5. Significance

This work provides a critical correction to the simulation of laser-matter interactions. By enforcing energy conservation in the context of excitation-dependent potentials, the authors show that previous models may have overestimated electron temperatures and ablation rates.

• Physical Accuracy: The study establishes that "bond hardening" and energy redistribution are dominant factors in the early stages of laser ablation, which static potentials miss.
• Methodological Standard: The proposed energy conservation algorithm and the rejection of artificial blast forces in favor of naturally emerging pressure gradients set a new standard for accurate TTM-MD simulations.
• Implications: These findings suggest that for accurate prediction of laser processing thresholds (melting, ablation depth) in metals like copper, simulations must account for the thermodynamic cost of changing interatomic potentials, not just the heat transfer.

In conclusion, the paper demonstrates that a self-consistent, energy-conserving approach with electron temperature-dependent potentials yields a more physically realistic description of laser-induced material dynamics, characterized by delayed phase transitions and reduced material removal compared to conventional methods.