Femtosecond X-Ray Diffraction Studies of the Reversal of the Microstructural Effects of Plastic Deformation during Shock Release of Tantalum

The thermodynamic-pathway-determined microstructure evolution of copper under shock compression

Shock-induced structural transformations in copper exhibit notable directional dependence and anisotropy, but the mechanisms that govern the responses of materials with different orientations are not yet well understood. In this study, we employ large-scale non-equilibrium molecular dynamics simulations to investigate the propagation of a shock wave through monocrystal copper and analyse the structural transformation dynamics in detail. Our results indicate that anisotropic structural evolution is determined by the thermodynamic pathway. A shock along the [Formula: see text] orientation causes a rapid and instantaneous temperature spike, resulting in a solid-solid phase transition. Conversely, a liquid metastable state is observed along the [Formula: see text] orientation due to thermodynamic supercooling. Notably, melting still occurs during the [Formula: see text]-oriented shock, even if it falls below the supercooling line in the thermodynamic pathway. These results highlight the importance of considering anisotropy, the thermodynamic pathway and solid-state disordering when interpreting phase transitions induced by shock. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. Our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

Spallation Characteristics of Single Crystal Aluminum with Copper Nanoparticles Based on Atomistic Simulations

In this study, the effects of Cu nanoparticle inclusion on the dynamic responses of single crystal Al during shockwave loading and subsequent spallation processes have been explored by molecular dynamics simulations. At specific impact velocities, the ideal single crystal Al will not produce dislocation and stacking fault structure during shock compression, while Cu inclusion in an Al–Cu nanocomposite will lead to the formation of a regular stacking fault structure. The significant difference of a shock-induced microstructure makes the spall strength of the Al–Cu nanocomposite lower than that of ideal single crystal Al at these specific impact velocities. The analysis of the damage evolution process shows that when piston velocity u_p ≤ 2.0 km/s, due to the dense defects and high potential energy at the interface between inclusions and matrix, voids will nucleate preferentially at the inclusion interface, and then grow along the interface at a rate of five times faster than other voids in the Al matrix. When u_p ≥ 2.5 km/s, the Al matrix will shock melt or unloading melt, and micro-spallation occurs; Cu inclusions have no effect on spallation strength, but when Cu inclusions and the Al matrix are not fully diffused, the voids tend to grow and coalescence along the inclusion interface to form a large void.

Semi-Automated Creation of Density Functional Tight Binding Models through Leveraging Chebyshev Polynomial-Based Force Fields

Density functional tight binding (DFTB) is an attractive method for accelerated quantum simulations of condensed matter due to its enhanced computational efficiency over standard density functional theory (DFT) approaches. However, DFTB models can be challenging to determine for individual systems of interest, especially for metallic and interfacial systems where different bonding arrangements can lead to significant changes in electronic states. In this regard, we have created a rapid-screening approach for determining systematically improvable DFTB interaction potentials that can yield transferable models for a variety of conditions. Our method leverages a recent reactive molecular dynamics force field where many-body interactions are represented by linear combinations of Chebyshev polynomials. This allows for the efficient creation of multi-center representations with relative ease, requiring only a small investment in initial DFT calculations. We have focused our workflow on TiH₂ as a model system and show that a relatively small training set based on unit-cell-sized calculations yields a model accurate for both bulk and surface properties. Our approach is easy to implement and can yield reliable DFTB models over a broad range of thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

Fingerprinting shock-induced deformations via diffraction

During the various stages of shock loading, many transient modes of deformation can activate and deactivate to affect the final state of a material. In order to fundamentally understand and optimize a shock response, researchers seek the ability to probe these modes in real-time and measure the microstructural evolutions with nanoscale resolution. Neither post-mortem analysis on recovered samples nor continuum-based methods during shock testing meet both requirements. High-speed diffraction offers a solution, but the interpretation of diffractograms suffers numerous debates and uncertainties. By atomistically simulating the shock, X-ray diffraction, and electron diffraction of three representative BCC and FCC metallic systems, we systematically isolated the characteristic fingerprints of salient deformation modes, such as dislocation slip (stacking faults), deformation twinning, and phase transformation as observed in experimental diffractograms. This study demonstrates how to use simulated diffractograms to connect the contributions from concurrent deformation modes to the evolutions of both 1D line profiles and 2D patterns for diffractograms from single crystals. Harnessing these fingerprints alongside information on local pressures and plasticity contributions facilitate the interpretation of shock experiments with cutting-edge resolution in both space and time.

Revisiting the Power Law Characteristics of the Plastic Shock Front under Shock Loading

Under uniaxial shock compression, the steepness of the plastic shock front usually exhibits power law characteristics with the Hugoniot pressure, also known as the "Swegle-Grady law." In this Letter, we show that the Swegle-Grady law can be described better by a third power law rather than the classical fourth power law at the strain rate between 10^{5}-10^{7}  s^{-1}. A simple dislocation-based continuum model is developed, which reproduced the third power law and revealed very good agreement with recent experiments of multiple types of metals quantitatively. New insights into this unusual macroscopic phenomenon are presented through quantifying the connection between the macroscopic mechanical response and the collective dynamics of dislocation assembles. It is found that the Swegle-Grady law results from the particular stress dependence of the plasticity behaviors, and that the difference between the third power scaling and the classical fourth power scaling results from different shock dissipative actions.

Development of shock-dynamics study with synchrotron-based time-resolved X-ray diffraction using an Nd:glass laser system

The combination of high-power laser and synchrotron X-ray pulses allows us to observe material responses under shock compression and release states at the crystal structure on a nanosecond time scale. A higher-power Nd:glass laser system for laser shock experiments was installed as a shock driving source at the NW14A beamline of PF-AR, KEK, Japan. It had a maximum pulse energy of 16 J, a pulse duration of 12 ns and a flat-top intensity profile on the target position. The shock-induced deformation dynamics of polycrystalline aluminium was investigated using synchrotron-based time-resolved X-ray diffraction (XRD) under laser-induced shock. The shock pressure reached up to about 17 GPa with a strain rate of at least 4.6 × 10⁷ s^-1 and remained there for nanoseconds. The plastic deformation caused by the shock-wave loading led to crystallite fragmentation. The preferred orientation of the polycrystalline aluminium remained essentially unchanged during the shock compression and release processes in this strain rate. The newly established time-resolved XRD experimental system can provide useful information for understanding the complex dynamic compression and release behaviors.

Unraveling the Role of Interfaces on the Spall Failure of Cu/Ta Multilayered Systems

Molecular dynamics (MD) simulations are carried out to investigate the effects of the type and spacing of FCC/BCC interfaces on the deformation and spall behavior. The simulations are carried out using model Cu/Ta multilayers with six different types of interfaces. The results suggest that interface type can significantly affect the structure and intensity of the incoming shock wave, change the activated slip systems, alter dislocation slip and twinning behavior, affect where and how voids are nucleated during spallation and the resulting spall strength. Moreover, the above aspects are significantly affected by the interface spacing. A transition from homogeneous to heterogeneous dislocation nucleation occurs as the interface spacing is decreased to 6 nm. Depending on interface type and spacing, damage (voids) nucleation and spall failure is observed to occur not only at the Cu/Ta interfaces, but also in the weaker Cu layer interior, or even in the stronger Ta layer interior, although different mechanisms underlie each of these three distinct failure modes. These findings point to the fact that, depending on the combination of interface type and spacing, interfaces can lead to both strengthening and weakening of the Cu/Ta multilayered microstructures.

Nonisentropic Release of a Shocked Solid

We present molecular dynamics simulations of shock and release in micron-scale tantalum crystals that exhibit postbreakout temperatures far exceeding those expected under the standard assumption of isentropic release. We show via an energy-budget analysis that this is due to plastic-work heating from material strength that largely counters thermoelastic cooling. The simulations are corroborated by experiments where the release temperatures of laser-shocked tantalum foils are deduced from their thermal strains via in situ x-ray diffraction and are found to be close to those behind the shock.

Extreme Hardening of Pb at High Pressure and Strain Rate

We study the high-pressure strength of Pb and Pb-4wt%Sb at the National Ignition Facility. We measure Rayleigh-Taylor growth of preformed ripples ramp compressed to ∼400  GPa peak pressure, among the highest-pressure strength measurements ever reported on any platform. We find agreement with 2D simulations using the Improved Steinberg-Guinan strength model for body-centered-cubic Pb; the Pb-4wt%Sb alloy behaves similarly within the error bars. The combination of high-rate, pressure-induced hardening and polymorphism yield an average inferred flow stress of ∼3.8  GPa at high pressure, a ∼250-fold increase, changing Pb from soft to extremely strong.

Deducing density and strength of nanocrystalline Ta and diamond under extreme conditions from X-ray diffraction

In situ X-ray diffraction with advanced X-ray sources offers unique opportunities for investigating materials properties under extreme conditions such as shock-wave loading. Here, Singh's theory for deducing high-pressure density and strength from two-dimensional (2D) diffraction patterns is rigorously examined with large-scale molecular dynamics simulations of isothermal compression and shock-wave compression. Two representative solids are explored: nanocrystalline Ta and diamond. Analysis of simulated 2D X-ray diffraction patterns is compared against direct molecular dynamics simulation results. Singh's method is highly accurate for density measurement (within 1%) and reasonable for strength measurement (within 10%), and can be used for such measurements on nanocrystalline and polycrystalline solids under extreme conditions (e.g. in the megabar regime).