This report is the final documentation for the one-year LDRD project 226360: Simulated X-ray Diffraction and Machine Learning for Optimizing Dynamic Experiment Analysis. As Sandia has successfully developed in-house X-ray diffraction tools for study of atomic structure in experiments, it has become increasingly important to develop computational analysis methods to support these experiments. When dynamically compressed lattices and orientations are not known a priori, the identification requires a cumbersome and sometimes intractable search of possible final states. These final states can include phase transition, deformation and mixed/evolving states. Our work consists of three parts: (1) development of an XRD simulation tool and use of traditional data science methods to match XRD patterns to experiments; (2) development of ML-based models capable of decomposing and identifying the lattice and orientation components of multicomponent experimental diffraction patterns; and (3) conducting experiments which showcase these new analysis tools in the study of phase transition mechanisms. Our target material has been cadmium sulfide, which exhibits complex orientation-dependent phase transformation mechanisms. In our current one-year LDRD, we have begun the analysis of high-quality c-axis CdS diffraction data from DCS and Thor experiments, which had until recently eluded orientation identification.
Abstract: Here, we report the high pressure phase and morphology behavior of ordered anatase titanium dioxide (TiO2) nanocrystal arrays. One-dimensional TiO2 nanorods and nanorices were synthesized and self-assembled into ordered mesostructures. Their phase and morphological transitions at both atomic scale and mesoscale under pressure were studied using in situ synchrotron wide- and small-angle x-ray scattering (WAXS and SAXS) techniques. At the atomic scale, synchrotron WAXS reveals a pressure-induced irreversible amorphization up to 35 GPa in both samples but with different onset pressures. On the mesoscale, no clear phase transformations were observed up to 20 GPa by synchrotron SAXS. Intriguingly, sintering of TiO2 nanorods at mesoscale into nano-squares or nano-rectangles, as well as nanorices into nanowires, were observed for the first time by transmission electron microscopy. Such pressure-induced nanoparticle phase-amorphization and morphological changes provide valuable insights for design and engineering structurally stable nanomaterials. Impact statement: The high pressure behavior of nanocrystals (NCs) continues to be of interest, as previous studies have demonstrated that an externally applied pressure can serve as an efficient tool to induce structural phase transitions of NC assemblies at both the atomic scale and mesoscale without altering any chemistry by manipulating NC interatomic and interparticle distances. In addition, the high pressure generated deviatoric stress has been proven to be able to force adjacent NCs to connect and fuse into new crystalline nanostructures. Although the atomic structural evolution of TiO2 NCs under pressure has been widely investigated in the past decades, open questions remain regarding the mesoscale phase transition and morphology of TiO2 NC assemblies as a function of pressure. Therefore, in this work, systemic high pressure experiments on ordered arrays of TiO2 nanorods and nanorices were conducted by employing wide/small angle x-ray scattering techniques. The sintering of TiO2 assemblies at mesoscale into various nanostructures under pressure were revealed by transmission electron microscopy. Overall, this high pressure work fills the current gap in research on the mesoscale phase behavior of TiO2 assemblies. The observed morphology tunability attained by applying pressure opens new pathways for engineering nanomaterials and optimizing their collective properties through mechanical compression stresses. Graphical abstract: [Figure not available: see fulltext.].
Pulsed-power generators using the magnetic loading technique are able to produce well-controlled continuous ramp compression of condensed matter for high-pressure equation-of-state studies. X-ray diffraction (XRD) data from dynamically compressed samples provide direct measurements of the elastic compression of the crystal lattice, onset of plastic flow, strength-strain rate dependence, structural phase transitions, and density of crystal defects such as dislocations. Here, we present a cost effective, compact X-ray source for XRD measurements on pulsed-power-driven ramp-loaded samples. This combination of magnetically-driven ramp compression of materials with single, short-pulse XRD diagnostic will be a powerful capability for the dynamic materials community. The success in fielding this new XRD diagnostic dramatically improves our predictive capability and understanding of rate-dependent behavior at or near phase transition. As Sandia plans the next-generation pulse-power driver platform, a key element needed to deliver new state-of-the-art experiments will be having the necessary diagnostic tools to probe new regimes and phenomena. These diagnostics need to be as versatile, compact, and portable as they are powerful. The development of a platform-independent XRD diagnostic gives Sandia researchers a new window to study the microstructure and phase dynamics of materials under load. This project has paved the way for phase transition research in a variety of materials with mission interest.
We describe recent efforts to improve our predictive modeling of rate-dependent behavior at, or near, a phase transition using molecular dynamics simulations. Cadmium sulfide (CdS) is a well-studied material that undergoes a solid-solid phase transition from wurtzite to rock salt structures between 3 and 9 GPa. Atomistic simulations are used to investigate the dominant transition mechanisms as a function of orientation, size and rate. We found that the final rock salt orientations were determined relative to the initial wurtzite orientation, and that these orientations were different for the two orientations and two pressure regimes studied. The CdS solid-solid phase transition is studied, for both a bulk single crystal and for polymer-encapsulated spherical nanoparticles of various sizes.
Recent dynamic compression experiments [M. D. Knudson, Science 348, 1455 (2015)10.1126/science.aaa7471; P. M. Celliers, Science 361, 677 (2018)10.1126/science.aat0970] have observed the insulator-metal transition in dense liquid deuterium, but with an approximately 95-GPa difference in the quoted pressures for the transition at comparable estimated temperatures. It was claimed in the latter of these two papers that a very large latent heat effect on the temperature was overlooked in the first, requiring correction of those temperatures downward by a factor of 2, thereby putting both experiments on the same theoretical phase boundary and reconciling the pressure discrepancy. We have performed extensive path-integral molecular dynamics calculations with density functional theory to directly calculate the isentropic temperature drop due to latent heat in the insulator-metal transition for dense liquid deuterium and show that this large temperature drop is not consistent with the underlying thermodynamics.
Magnetically launched flyer plates were used to investigate the shock response of beryllium between 90 and 300 GPa. Solid aluminum flyer plates drove steady shocks into polycrystalline beryllium to constrain the Hugoniot from 90 to 190 GPa. Multilayered copper/aluminum flyer plates generated a shock followed by an overtaking rarefaction which was used to determine the sound velocity in both solid and liquid beryllium between 130 and 300 GPa. Disappearance of the longitudinal wave was used to identify the onset of melt along the Hugoniot and measurements were compared to density functional theory calculations to explore the proposed hcp-bcc transition at high pressure. The onset of melt along the Hugoniot was identified at ∼205GPa, which is in good agreement with theoretical predictions. These results show no clear indication of an hcp-bcc transition prior to melt along the beryllium Hugoniot. Rather, the shear stress, determined from the release wave profiles, was found to gradually decrease with stress and eventually vanish at the onset of melt.
Measurement of the Hugoniot and sound velocity provides information on the bulk modulus and Grüneisen parameter of a material at extreme conditions. The capability to launch multilayered (copper/aluminum) flyer plates at velocities in excess of 20 km/s with the Sandia Z accelerator has enabled high-precision sound-velocity measurements at previously inaccessible pressures. For these experiments, the sound velocity of the copper flyer must be accurately known in the multi-Mbar regime. Here we describe the development of copper as an absolutely calibrated sound-velocity standard for high-precision measurements at pressures in excess of 400 GPa. Using multilayered flyer plates, we performed absolute measurements of the Hugoniot and sound velocity of copper for pressures from 500 to 1200 GPa. These measurements enabled the determination of the Grüneisen parameter for dense liquid copper, clearly showing a density dependence above the melt transition. Combined with earlier data at lower pressures, these results constrain the sound velocity as a function of pressure, enabling the use of copper as a Hugoniot and sound-velocity standard for pressures up to 1200 GPa.
Measurement of the window interface velocity is a common technique for investigating the dynamic response materials at high strain rates. However, these measurements are limited in pressure to the range where the window remains transparent. The most common window material for this application is lithium fluoride, which under single shock compression becomes opaque at ∼200 GPa. To date, no other window material has been identified for use at higher pressures. Here, we present a Lagrangian technique to calculate the interface velocity from a continuously measured shock velocity, with application to quartz. The quartz shock front becomes reflective upon melt, at ∼100 GPa, enabling the use of velocity interferometry to continuously measure the shock velocity. This technique overlaps with the range of pressures accessible with LiF windows and extends the region where wave profile measurements are possible to pressures in excess of 2000 GPa. We show through simulated data that the technique accurately reproduces the interface velocity within 20% of the initial state, and that the Lagrangian technique represents a significant improvement over a simple linear approximation.
In recent years, α-quartz has been used prolifically as an impedance matching standard in shock wave experiments in the multi-Mbar regime (1 Mbar = 100 GPa = 0.1 TPa). This is due to the fact that above ∼90-100 GPa along the principal Hugoniot α-quartz becomes reflective, and thus, shock velocities can be measured to high precision using velocity interferometry. The Hugoniot and release of α-quartz have been studied extensively, enabling the development of an analytical release model for use in impedance matching. However, this analytical release model has only been validated over a range of 300-1200 GPa (0.3-1.2 TPa). Here, we extend this analytical model to 200-3000 GPa (0.2-3 TPa) through additional α-quartz Hugoniot and release measurements, as well as first-principles molecular dynamics calculations.
Materials in nuclear and conventional weapons can reach multi-megabar pressures and 1000s of degree temperatures on timescales ranging from microseconds to nanoseconds. Understanding the response of complex materials under these conditions is important for designing and assessing changes to nuclear weapons. In the next few decades, a major concern will be evaluating the behavior of aging materials and remanufactured components. The science to enable the program to underwrite decisions quickly and confidently on use, remanufacturing, and replacement of these materials will be critical to NNSA’s new Stockpile Responsiveness Program. Material response is also important for assessing the risks posed by adversaries or proliferants. Dynamic materials research, which refers to the use of high-speed experiments to produce extreme conditions in matter, is an important part of NNSA’s Stockpile Stewardship Program.