Dynamic stress response kernels for dislocations and cracks: unified anisotropic Lagrangian formulation

This paper presents a unified anisotropic Lagrangian formulation for dynamic stress response kernels of dislocations and cracks using the Stroh formalism, deriving their space-time and Fourier representations in terms of the prelogarithmic Lagrangian factor and its derivative to enable causality-compliant numerical simulations across all motion regimes.

The Gist

Imagine you are walking through a crowded room. If you walk slowly, people around you just shift slightly to make room, and the disturbance settles down quickly. But if you start running, or even better, if you are a "ghost" moving faster than the speed of sound in that room, the way the crowd reacts changes completely. You create shockwaves, and the energy you expend doesn't just disappear; it radiates away, dragging on you.

This paper is about understanding exactly how that "drag" and those "shockwaves" work for tiny defects inside solid materials, like cracks or dislocations (which are like missing rows of atoms in a crystal).

Here is a breakdown of the paper's big ideas using simple analogies:

1. The Problem: The "Crystal Traffic Jam"

Materials aren't always uniform. Some are like a perfectly smooth pond (isotropic), while others are like a wooden floor with distinct grain (anisotropic).

• The Old Way: Scientists used to study cracks and moving defects assuming the material was like a smooth pond. They had a perfect map (math) for how waves travel in that simple world.
• The New Reality: Real materials (like the metal in a jet engine or the silicon in a chip) are more like the wooden floor. The waves travel differently depending on the direction. The old maps didn't work here, and scientists were flying blind when trying to predict how fast a crack would grow or how a defect would move in these complex materials.

2. The Solution: A Universal "Traffic Controller"

The authors (Pellegrini, Josien, and Chassard) have built a new, universal map. They used a sophisticated mathematical tool called the Stroh formalism (think of it as a high-tech GPS for crystal structures) to figure out how stress waves behave in any direction, at any speed.

They discovered that all the complex math describing these waves boils down to a single, elegant function they call $L(v)$.

• The Analogy: Imagine $L(v)$ is the "speedometer" of the material's resistance. It tells you how much energy is stored in the defect and how much is being lost to the surroundings as it moves.
• The Twist: This function has a secret "imaginary" part. In math, this sounds scary, but in physics, it's the braking system. It represents the energy that leaks away as sound waves (radiation) when the defect moves fast.

3. The Three Speed Zones

The paper covers three distinct regimes of motion, which the authors handle with a single unified theory:

1. Sub-sonic (The Jogger): Moving slower than the speed of sound in the material. The waves ripple out ahead and behind, but the defect is in control.
2. Intersonic (The Speedster): Moving faster than some sound waves but slower than others. This is a chaotic zone where the material behaves strangely, like a boat creating a wake that it can't outrun.
3. Supersonic (The Jet): Moving faster than all sound waves. This creates a "Mach cone" (like a sonic boom). The paper shows how to calculate the stress even in this extreme regime.

4. The "Ghost" Connection: Energy vs. Force

One of the paper's coolest findings is connecting two things that usually seem unrelated: Energy and Force.

• Energy ($L$): How much "fuel" the defect has.
• Force (Stress): How hard the material pushes back.
• The Discovery: The authors proved that if you know the "fuel" function ($L$), you can instantly calculate the "push back" (stress) just by taking its derivative (a math operation that tells you how fast the fuel is changing).
• The Metaphor: It's like knowing that if you know exactly how much gas a car burns at different speeds, you can instantly calculate the air resistance pushing against the car without ever measuring the wind.

5. Why This Matters (The "So What?")

Why do we care about math for moving cracks?

• Safety: If you are designing a bridge or a nuclear reactor, you need to know exactly how fast a crack will grow before it causes a catastrophic failure.
• Technology: In computer chips, tiny defects can ruin performance. Understanding how they move helps engineers make better, faster chips.
• Simulation: The authors designed their math specifically to be used in computer simulations (like a video game physics engine). This means engineers can now run accurate simulations of cracks in complex materials much faster and more reliably than before.

Summary

Think of this paper as writing the instruction manual for the universe's most complex traffic system. Before, we only knew the rules for a flat, empty road. Now, thanks to this paper, we have the rules for driving a race car on a bumpy, twisting mountain track at supersonic speeds. They found a single "magic number" ($L$) that predicts exactly how the road will push back, ensuring we can predict when things will break and how to stop it.

Technical Summary

1. Problem Statement

The paper addresses a significant gap in the elastodynamic theory of moving defects (dislocations and cracks) in anisotropic media.

• Context: While the stress response of moving defects is well-understood in isotropic elasticity (using generalized functions and specific kernels like the Geubelle-Rice model for cracks or the Dynamic Peierls Equation for dislocations), these formulations rely on the prelogarithmic Lagrangian factor $L(v)$ and its derivative.
• The Gap: For anisotropic elasticity, the space-time representations of these stress-response kernels were previously unknown. Existing theoretical studies often rely on isotropic approximations, which fail to capture the complex wave propagation and radiative losses inherent in real crystalline materials.
• Specific Challenge: There is a need to unify the description of radiative stress responses for both cracks and dislocations in anisotropic media, covering subsonic, intersonic, and supersonic regimes, while strictly adhering to causality constraints.

2. Methodology

The authors employ a rigorous mathematical framework combining Fourier Transform (FT) techniques with the Stroh formalism (a powerful method for solving anisotropic elasticity problems).

• Stroh Formalism with Causality: The authors utilize the Stroh formalism to solve the Navier equation for an elastic half-space. Crucially, they implement the Principle of Limiting Absorption (PLA) by replacing the frequency $\omega$ with $\omega + i0^+$. This ensures the solution satisfies causality and correctly handles radiative drag in intersonic/supersonic regimes.
• Fourier Domain Derivation: Instead of solving directly in real space, the problem is formulated in the Fourier domain (wavevector $\mathbf{k}$ and frequency $\omega$). The stress traction $\mathbf{t}$ is related to the displacement discontinuity (slip or crack opening) $\boldsymbol{\eta}$ via a convolution kernel.
• Eigenvalue Analysis: The core of the derivation involves solving a sextic eigenvalue problem to find the Stroh eigenvalues $p_\alpha$ and eigenvectors ($\mathbf{A}_\alpha, \mathbf{L}_\alpha$). The authors categorize these based on their imaginary parts to distinguish between surface-localized modes (subsonic) and bulk radiative modes (intersonic/supersonic).
• Analytic Continuation: The theory relies on the analytic continuation of the prelogarithmic Lagrangian factor $L(v)$ and its derivative (the prelogarithmic impulsion function $p(v) = L'(v)$) into the complex velocity plane.

3. Key Contributions

A. Unified Anisotropic Stress-Response Kernel

The paper derives a unified expression for the stress-response kernel for both dislocations and cracks in anisotropic media.

• Fourier Representation: The traction-discontinuity relationship is expressed as:
  $$\mathbf{t}(\mathbf{k}, \omega) = (2\pi)k \, \mathbf{L}(\hat{\mathbf{k}}, \frac{\omega}{k} + i0^+) \cdot \boldsymbol{\eta}(\mathbf{k}, \omega)$$
  where $\mathbf{L}$ is the Lagrangian energy tensor, derived entirely from Stroh eigenvectors.
• Space-Time Representation: The authors establish the space-time form of the kernel, showing it depends on the derivative of the Lagrangian factor, $p(v) = L'(v)$.

B. The Role of the Lagrangian Factor $L(v)$

A major theoretical insight is that the entire radiative in-plane stress response can be formulated exclusively in terms of:

1. The prelogarithmic Lagrangian factor $L(v)$ (related to energy).
2. Its first derivative $p(v) = L'(v)$ (the impulsion function).
   The paper clarifies that the imaginary part of $L(v + i0^+)$, often associated with radiative dissipation, arises naturally from the causal analytic continuation required to handle the $i0^+$ term.

C. Reformulation of the Weertman Equation

For steady-state motion, the authors derive the stress kernel for the Weertman equation (governing steady dislocation motion) explicitly in terms of $L(v)$. This provides a direct link between the energetic properties of the defect and its mobility law in anisotropic materials.

D. Instantaneous Radiative Term

The paper derives the anisotropic radiative drag coefficient (instantaneous radiative term) for the Dynamic Peierls Equation.

• It is shown that the radiative term corresponds to a local drag force proportional to the time derivative of the slip: $\mathbf{t}_{rad} = -\mathbf{K} \cdot \partial_t \boldsymbol{\eta}$.
• The tensor $\mathbf{K}$ is derived as a limit of the Lagrangian tensor and is shown to be independent of the wavevector direction, depending only on the material density and wave speeds along the normal.

4. Key Results

• Generalized Kernel: The stress response kernel in space-time is given by:
  $$\sigma_r(x, t) = -\frac{\mu}{\pi} \int dx' \int dt' K(x-x', t-t') \frac{\partial \eta}{\partial x'} - \kappa \frac{\mu}{2c} \frac{\partial \eta}{\partial t}$$
  where the non-local kernel $K(x,t)$ is directly related to the impulsion function:
  $$K(x, t) \equiv \theta(t) \frac{1}{2w_0 t^2} \text{Re} \left[ p\left(\frac{x}{t} + i0^+\right) \right]$$
• Regime Coverage: The formulation is valid for subsonic, intersonic, and supersonic velocities. The use of $i0^+$ automatically handles the transition between these regimes, including the formation of Mach cones in supersonic motion.
• Isotropic Limit: The derived anisotropic expressions correctly reduce to known isotropic results (e.g., specific $\kappa$ coefficients for screw, edge, and climb dislocations) when the elastic tensor is isotropic.
• Green's Function Connection: The paper proves that the causal stress-response kernel differs from the purely Lagrangian kernel only by its radiative part, explaining why analytic continuation is necessary to retrieve physical stress fields from energy-based Lagrangians.

5. Significance and Applications

• Numerical Implementation: The theory is specifically tailored for phase-field-type numerical codes that utilize Fast Fourier Transforms (FFT) in space and Laplace transforms in time. By expressing the kernels in terms of $L(v)$ and $p(v)$, the paper provides a computationally efficient route to simulate defect dynamics in anisotropic crystals without solving complex integro-differential equations in real space.
• Physical Insight: It resolves long-standing ambiguities regarding the imaginary part of the Lagrangian factor, demonstrating that it is not merely a mathematical artifact but a direct representation of radiative energy loss and causality.
• Unified Framework: It bridges the gap between dislocation dynamics (Peierls-Nabarro models) and fracture mechanics (cohesive zone models), showing they share the same underlying elastodynamic kernel structure in anisotropic media.
• Mobility Laws: The results allow for the immediate transfer of "Collective-Variable-Approximation" equations of motion (previously limited to isotropic cases) to anisotropic materials, enabling more accurate predictions of defect mobility at ultra-high speeds.

In summary, this paper provides the missing theoretical foundation for modeling fast-moving defects in anisotropic materials, unifying energy-based Lagrangian approaches with causal stress-response kernels, and offering a practical framework for advanced numerical simulations.