Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of NCERC Critical Experiments with ENDF/B Nuclear Data Libraries

Nuclear data (ND) libraries are the backbone of the nuclear industry, as they are the collections of tabulated probabilities that define sub-atomic particle interactions with matter. In the areas of criticality safety and experiments, no evaluated nuclear data files (ENDF) are more important than those containing neutron cross section data. It is upon these files, and accompanying radiation transport codes, that practitioners are enabled to safely design both subcritical and critical systems. Likewise, in a symbiotic fashion, it is the same critical assemblies which are used primarily to validate that the cross sections are correct. The ENDF/B library, the United States’ national library maintained by the National Nuclear Data Center (NNDC) at Brookhaven National Laboratory (BNL), is soon releasing a new version, ENDF/B-VIII.1. Prior to the official release, several beta versions of the library were prepared and tested in simulation suites. The work presented here are results comparing the newest ENDF/B beta release (ENDF/B-VIII.1b3) and ENDF/B-VIII.0 with recent experiments done at the National Criticality Experiments Research Center (NCERC) as well as correlated experiments from the Los Alamos Critical Experiments Facility (LACEF). These NCERC and LACEF experiments were performed in part to provide validation for various cross sections that were identified as insufficient in the ENDF/B-VIII.0. In particular, lead, copper, fluorine, and chlorine, as well as the major actinides, were targeted from the last decade of critical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Atomic arrangements in an amorphous CoFeB ribbon extracted via an analysis of radial distribution functions

Here we discuss the atomic structure of amorphous ferromagnetic FeCoB alloys, which are used widely in spintronics applications. Specifically, we obtain the pair-distribution functions for various atomic pairs based on high-energy x-ray diffraction data taken from an amorphous Co 20 Fe 61 B 19 specimen. We start our reverse Monte Carlo cycles to determine the disordered structure with a two-phase model in which a small amount of cobalt is mixed with Fe 23 B 6 as a second phase. The structure of the alloy is found to be heterogeneous, where the boron atoms drive disorder through the random occupation of the atomic network. Our analysis also indicates the presence of small cobalt clusters that are embedded in the iron matrix and percolating the latter throughout the structure. This morphology can explain the enhanced spin polarization observed in amorphous magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Oriented Electrostatic Effects on O 2 and CO 2 Reduction by a Polycationic Iron Porphyrin

Next-generation energy technologies require improved methods for rapid and efficient chemical-to-electrical energy transformations. One new approach has been to include atomically-positioned, electrostatic motifs in molecular catalysts to stabilize high-energy, charged intermediates. For example, an iron porphyrin bearing four cationic, ortho-N,N,N-trimethylanilinium groups (o-[N(CH 3 ) 3 ] + ) has recently been used to catalyze the complex, multi-step O 2 and CO 2 reduction reactions (ORR and CO 2 RR) with fast rates and at low overpotentials. The success of this catalyst is attributed, at least in part, to specific charge-charge interactions between the atomically-positioned o-[N(CH 3 ) 3 ] + groups and the bound substrate. However, by nature of the mono-ortho substitution pattern, there are four possible atropisomers of this metalloporphyrin and thus four unique electrostatic environments. This work reports that each of the four individual atropisomers catalyzes both the ORR and CO 2 RR with fast rates and low overpotentials. The maximum turnover frequencies (TOF max ) vary among the atropisomers, by a factor of 60 for the ORR and a factor of 5 for CO 2 RR. For the ORR, the αβαβ isomer is the fastest and has the highest overpotential, while for the CO 2 RR, the order is reversed and the αααα isomer is the fastest and has the highest overpotential. The role of charge-positioning is complex and can affect more than a single step such as CO 2 binding. These data offer a first-of-a-kind perspective on atomically positioned charge and highlight the significance of high charge density, rather than orientation, on the thermodynamics and kinetics of multistep molecular electrochemical transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2–3 × 10 −4 (Δ E / E ). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (Δ E / E = 3 × 10 −3 ) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Erasure conversion in a high-fidelity Rydberg quantum simulator

Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices and for the quest towards fault-tolerant quantum computation. Rydberg arrays have emerged as a prominent platform in this context with impressive system sizes and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors with single-atom resolution, a form of erasure error conversion. However, two-qubit entanglement fidelities in Rydberg atom arrays have lagged behind competitors and this type of erasure conversion is yet to be realized for matter-based qubits in general. Here we demonstrate both erasure conversion and high-fidelity Bell state generation using a Rydberg quantum simulator. When excising data with erasure errors observed via fast imaging of alkaline-earth atoms, we achieve a Bell state fidelity of $\ge 0.997{1}_{-13}^{+10}$ , which improves to $\ge 0.998{5}_{-12}^{+7}$ when correcting for remaining state-preparation errors. We further apply erasure conversion in a quantum simulation experiment for quasi-adiabatic preparation of long-range order across a quantum phase transition, and reveal the otherwise hidden impact of these errors on the simulation outcome. Our work demonstrates the capability for Rydberg-based entanglement to reach fidelities in the 0.999 regime, with higher fidelities a question of technical improvements, and shows how erasure conversion can be utilized in NISQ devices. These techniques could be translated directly to quantum-error-correction codes with the addition of long-lived qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Oxygen Defect Configurations in Single-Phase UO 2.15

Hyperstoichiometric UO 2.15 was characterized by neutron total scattering at high temperature in the single-phase UO 2+x region of the U/O phase diagram. The diffraction data confirmed a single-phase fluorite structure at high temperature. Analysis of the short-range data showed that the same structural model does not fit the pair distribution functions well. Instead, structural models containing specific configurations of oxygen defect clusters best represent the local atomic arrangement. In conclusion, prevalent defect clusters previously proposed were fit to the experimental data, and moderately distorted oxygen cuboctahedra hypothesized by recent molecular dynamics simulations fit the data most accurately.

neutron pair-distribution function analysis↗
DOE OSTI2023

Structure of amorphous materials in the NASICON system Na 1+x Ti 2 Si x P 3−x O 12

Abstract The structure of glasses in the sodium (Na) super-ionic conductor (NASICON) system Na 1 + x Ti 2 Si x P 3 − x O 12 with x = 0.8 and x = 1.0 was explored by combining neutron and high-energy x-ray diffraction with 29 Si, 31 P and 23 Na solid-state nuclear magnetic resonance (NMR) spectroscopy. The 29 Si magic angle spinning (MAS) NMR spectra reveal that the silica component remains fully polymerized in the form of Si 4 units, i.e. the silicon atoms are bound to four bridging oxygen atoms. The 31 P{ 23 Na} rotational echo adiabatic passage double resonance (REAPDOR) NMR data suggest that the 31 P MAS NMR line shape originates from four-coordinated P n units, where n = 1, 2 or 3 is the number of bridging oxygen atoms per phosphorus atom. These sites differ in their 31 P- 23 Na dipolar coupling strengths. The results support an intermediate range order scenario of a phosphosilicate mixed network-former glass in which the phosphate groups selectively attract the Na + modifier ions. Titanium takes a sub-octahedral coordination environment with a mean Ti–O coordination number of 5.17(4) for x = 0.8 and 4.86(4) for x = 1.0. A mismatch between the P–O and Si–O bond lengths of 8% is likely to inhibit the incorporation of silicon into the phosphorus sites of the NASICON crystal structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modification of a force field for molecular dynamics simulations of silicon etching by chlorine atoms

Here, a modified classical molecular dynamics (MD) force field is presented for the interaction between silicon (Si) and chlorine (Cl). The original version of the force field is shown to significantly overestimate the probability of Si etching by thermal Cl atoms. However, the modified force field corrects this problem and results in generally good agreement with experimental data. Further, it is shown that while the modification of the force field improves the prediction of Si spontaneous etching with Cl atoms, it does not degrade predictions of atomic-layer etching of Si with Cl[Formula: see text] molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Resilience of the Aurivillius structure upon La and Cr doping in a Bi 5 Ti 3 FeO 15 multiferroic

Here, combining the experimental techniques of high-resolution X-ray diffraction, magnetometry, specific heat measurement, and X-ray photoelectron, Raman and dielectric spectroscopy techniques, we have studied the influence of La and Cr doping on the crystal structure and magnetism of the room temperature Aurivillius multiferroic Bi 5 Ti 3 FeO 15 by investigating the physical properties of (Bi 4 La)Ti 3 FeO 15 and Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 . The parent (Bi 5 Ti 3 FeO 15 ) and the doped ((Bi 4 La)Ti 3 FeO 15 and Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 ) compounds crystallize in the A2 1 am space group, which is confirmed through our analysis of high-resolution synchrotron X-ray diffraction data obtained on phase-pure polycrystalline powders. We determined the oxidation states of the metal atoms in the studied compounds as Fe 3+ , Ti 4+ , Cr 3+ , and La 3+ through the analysis of X-ray photoelectron spectroscopy data. The magnetic susceptibilities of the three compounds are marked by the absence of a long-range ordered ground state, but dominated by superparamagnetic clusters with dominant antiferromagnetic interactions. This signature of short-range magnetism is also seen in specific heat as a low temperature enhancement which is suppressed upon the application of external magnetic fields up to 8 T. Our dielectric spectroscopy experiments showed that the three studied compounds display similar features in the dielectric constant measured as a function of frequency. However, upon doping La at the Bi site, the width of the ferroelectric hysteresis loop increases for (Bi 4 La)Ti 3 FeO 15 compared to that of the parent compound Bi 5 Ti 3 FeO 15 , and with Cr doping, Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 becomes a leaky dielectric. The resilience of the Aurivillius crystal structure towards doping of La at the Bi site and Cr at the Fe site is clearly seen in the bulk properties of magnetic susceptibility, specific heat and the average crystal structure. The relevance of changes in the local structure is evident from our Raman spectroscopy and X-ray pair distribution function studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accelerating error correction in tomographic reconstruction

Abstract Spurred by recent advances in detector technology and X-ray optics, upgrades to scanning-probe-based tomographic imaging have led to an exponential growth in the amount and complexity of experimental data and have created a clear opportunity for tomographic imaging to approach single-atom sensitivity. The improved spatial resolution, however, is highly susceptible to systematic and random experimental errors, such as center of rotation drifts, which may lead to imaging artifacts and prevent reliable data extraction. Here, we present a model-based approach that simultaneously optimizes the reconstructed specimen and sinogram alignment as a single optimization problem for tomographic reconstruction with center of rotation error correction. Our algorithm utilizes an adaptive regularizer that is dynamically adjusted at each alternating iteration step. Furthermore, we describe its implementation in a software package targeting high-throughput workflows for execution on distributed-memory clusters. We demonstrate the performance of our solver on large-scale synthetic problems and show that it is robust to a wide range of noise and experimental drifts with near-ideal throughput.

Ali, Sajid (ORCID:0000000321864636)↗
DOE OSTI2024

Characterization of Fuel Cladding Chemical Interaction on a High Burnup U-10Zr Metallic Fuel via Electron Energy Loss Spectroscopy Enhanced by Machine Learning

Fuel cladding chemical interaction (FCCI) is one of the main performance limiting factors for metallic nuclear fuels. The interaction destabilizes the martensitic microstructure and deteriorates mechanical properties of HT-9 cladding. The detection of low atomic number elements (Z<10) and overlapping of elemental peaks can be problematic in interpreting energy dispersive X-ray spectroscopy (EDS) data. Electron energy loss spectroscopy (EELS) provides precise elemental edge energy values and can detect elements with a low atomic number. This work utilizes EELS to study the distribution of lanthanides and light elements at the interaction region. The sample was prepared from the FCCI region of a U-10Zr (wt.%) solid fuel with HT-9 cladding, irradiated to a burnup of 13.2 at.%. Processing the EELS data included three major steps: 1) enhance the signal to noise ratio by denoising the spectrum with principal component analysis (PCA) method, removing background and performing deconvolution; 2) identify chemical elements with core energy loss edges; 3) confirm different phases using a popular machine learning method, K-means. This work presents qualitative assessment of lanthanides and light elements like carbon (C) and oxygen (O) enhanced by the application of machine learning algorithms. By comparing with EDS elemental maps, EELS provides higher resolution chemical maps, reveals the distribution of carbon at the interaction region supporting the formation of zirconium carbide, a rind-like microstructure feature that was proposed to mitigate the chemical interaction. Furthermore, the plasmon peak map was also found to indicate an energy shift associated with the formation of phases/compounds. K-means clustering method was used on the processed electron energy loss (EEL) spectrum to automatically reveal different phases. The resulting clustered maps from K-means clustering align well with elemental maps confirming certain phases, especially Fe-Ce and Zr-C, in the FCCI region.

EELS↗
DOE OSTI2024

Enhancing Docking Accuracy with PECAN2, a 3D Atomic Neural Network Trained without Co-Complex Crystal Structures

Decades of drug development research have explored a vast chemical space for highly active compounds. The exponential growth of virtual libraries enables easy access to billions of synthesizable molecules. Computational modeling, particularly molecular docking, utilizes physics-based calculations to prioritize molecules for synthesis and testing. Nevertheless, the molecular docking process often yields docking poses with favorable scores that prove to be inaccurate with experimental testing. To address these issues, several approaches using machine learning (ML) have been proposed to filter incorrect poses based on the crystal structures. However, most of the methods are limited by the availability of structure data. Here, we propose a new pose classification approach, PECAN2 (Pose Classification with 3D Atomic Network 2), without the need for crystal structures, based on a 3D atomic neural network with Point Cloud Network (PCN). The new approach uses the correlation between docking scores and experimental data to assign labels, instead of relying on the crystal structures. We validate the proposed classifier on multiple datasets including human mu, delta, and kappa opioid receptors and SARS-CoV-2 Mpro. Our results demonstrate that leveraging the correlation between docking scores and experimental data alone enhances molecular docking performance by filtering out false positives and false negatives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Atomic Energy Accuracy of Neural Network Potentials: Harnessing Pretraining and Transfer Learning

Machine learning-based interatomic potentials (MLIPs) have transformed the prediction of potential energy surfaces (PESs), achieving accuracy comparable to ab initio calculations. However, atomic energy predictions, often assumed to lack physical meaning, remain underexplored. In this study, we demonstrate that inaccuracies in atomic energy predictions reduce the robustness and transferability of Neural Network Potentials (NNPs) and atomic energy error can be masked in total energy predictions due to error cancellation. Here, we validate this finding using challenging configurations involving deformation and failure under tensile loading. By pretraining atomic energy predictions using empirical potentials and applying transfer learning with density functional theory (DFT) data, we achieve notable improvements in the accuracy of total energy, forces, and stress predictions. Furthermore, this approach enhances the robustness and transferability of NNPs, emphasizing the importance of atomic energy predictions in developing high-quality and reliable MLIPs.

Active Learning↗
DOE OSTI2022

Time-resolved CO 2 , CO, and N 2 vibrational population measurements in Ns pulse discharge plasmas

Abstract Time-resolved CO 2 and N 2 vibrational populations and translational-rotational temperature are measured in a CO 2 –N 2 plasma sustained by a ns pulse discharge burst in plane-to-plane geometry. Time-resolved, absolute number density of CO generated in the plasma is also inferred from the experimental data. CO 2 and CO vibrational populations are measured by mid-IR, tunable quantum cascade laser absorption spectroscopy, and N 2 vibrational populations are measured by the ns broadband vibrational CARS. Transient excitation of N 2 and CO 2 asymmetric stretch vibrational energy modes is detected during the discharge burst. The time-resolved rate of CO generation does not correlate with N 2 or CO 2 ( ν 3 ) vibrational temperatures, indicating that CO 2 dissociation via the vibrational excitation is insignificant at the present conditions. The rate of CO generation decreases gradually during the discharge burst. The estimated specific energy cost of the CO product is close to that of N atoms in pure nitrogen, measured previously at similar operating conditions. Comparison of the experimental data with the kinetic modeling analysis indicates that CO 2 dissociation in collisions with electronically excited N 2 molecules is the dominant channel of CO generation at the present conditions, although the inferred CO yield in these processes is significantly lower than 1. The effect of vibrational energy transfer between N 2 and CO 2 on the plasma chemical processes is insignificant. The kinetic model underpredicts a rapid reduction of the N 2 and CO 2 ( ν 3 ) vibrational temperatures during the later half of the discharge burst and in the afterglow. V–T relaxation of N 2 by N and O atoms generated in the ns pulse discharge plasma does not affect the vibrational relaxation rate in a significant way. However, rapid V–T relaxation of CO 2 by O atoms has a significant effect on the relaxation rate. The difference between the experimental data and the modeling predictions may be due to the unknown scaling of the CO 2 –O V–T rates with the vibrational quantum number.

Physics↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗