Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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 1,459 records · Page 81

A quantum processor based on coherent transport of entangled atom arrays

The ability to engineer parallel, programmable operations between desired qubits within a quantum processor is key for building scalable quantum information systems. In most state-of-the-art approaches, qubits interact locally, constrained by the connectivity associated with their fixed spatial layout. Here we demonstrate a quantum processor with dynamic, non-local connectivity, in which entangled qubits are coherently transported in a highly parallel manner across two spatial dimensions, between layers of single- and two-qubit operations. Our approach makes use of neutral atom arrays trapped and transported by optical tweezers; hyperfine states are used for robust quantum information storage, and excitation into Rydberg states is used for entanglement generation. We use this architecture to realize programmable generation of entangled graph states, such as cluster states and a seven-qubit Steane code state. Furthermore, we shuttle entangled ancilla arrays to realize a surface code state with thirteen data and six ancillary qubits and a toric code state on a torus with sixteen data and eight ancillary qubits. Finally, we use this architecture to realize a hybrid analogue–digital evolution and use it for measuring entanglement entropy in quantum simulations, experimentally observing non-monotonic entanglement dynamics associated with quantum many-body scars. Realizing a long-standing goal, these results provide a route towards scalable quantum processing and enable applications ranging from simulation to metrology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗

Dynamical Bond Formation in KNi 2 Se 2

Emphanisis, or the appearance out of nothing, has been used to describe the phenomena of spontaneous atom off-centering and dipole formation at elevated temperatures in lead chalcogenides. Here, we provide spectroscopic evidence of spontaneous formation of metal-metal bonds above T~50 K in the layered metal KNi 2 Se 2 . These bonds form zig-zag chains that lower the local symmetry from tetragonal to monoclinic. Energy-resolved pair distribution function measurements exclude a pure phonon origin of our observations, and instead imply the existence of extra, slowly fluctuating, Ni-Ni bonds above T=50 K. Density functional theory calculations support this lower symmetry configuration as an instability of tetragonal KNi 2 Se 2 . We thus demonstrate that the phenomena of emphansis is not limited to local electric dipole formation, but can also be driven by the formation of metal-metal bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗

Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging

Abstract The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D + + OD + intermediates, from other processes leading to the same D + + D + + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1 Σ + and a 1 Δ electronic states of the intermediate OD + ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD + intermediate.

74 ATOMIC AND MOLECULAR PHYSICS↗

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the importance of excited state species in low pressure capacitively coupled plasma argon discharges

In the past three decades, first principles-based fully kinetic particle-in-cell Monte Carlo collision (PIC/MCC) simulations have been proven to be an important tool for the understanding of the physics of low pressure capacitive discharges. However, there is a long-standing issue that the plasma density determined by PIC/MCC simulations shows quantitative deviations from experimental measurements, even in argon discharges, indicating that certain physics may be missing in previous modeling of the low pressure radio frequency (rf) driven capacitive discharges. In this work, we report that the energetic electron-induced secondary electron emission (SEE) and excited state atoms play an important role in low pressure rf capacitive argon plasma discharges. The ion-induced secondary electrons are accelerated by the high sheath field to strike the opposite electrode and produce a considerable number of secondary electrons that lead to additional ionizing impacts and further increase of the plasma density. Importantly, the presence of excited state species even further enhances the plasma density via excited state neutral and resonant state photon-induced SEE on the electrode surface. The PIC/MCC simulation results show good agreement with the recent experimental measurements in the low pressure range (1–10 Pa) that is commonly used for etching in the semiconductor industry. At the highest pressure (20 Pa) and driving voltage amplitudes 250 and 350 V explored here, the plasma densities from PIC/MCC simulations considering excited state neutrals and resonant photon-induced SEE are quantitatively higher than observed in the experiments, requiring further investigation on high pressure discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Curvature-induced hydrophobicity at imogolite–water interfaces

Imogolite, a nanotubular aluminosilicate mineral, is commonly found in volcanic soils, where it exerts a control on carbon dynamics. Synthetic imogolites are used for the removal of contaminants from industrial effluents and are considered for a range of other applications including gas adsorption and functionalised heterogeneous catalysts. In spite of their environmental and industrial relevance, the properties of imogolite–water interfaces remain poorly understood. In this work, an experimental and computational study is presented in which the structure and energetics of water are characterized on the curved external surface of imogolite and the hydrophilicity of this surface is contrasted with that of gibbsite, its planar counterpart. Atomic force spectroscopy experiments show that in spite of their identical surface structure, imogolite has a lower hygroscopicity than gibbsite. Molecular dynamics simulations provide an explanation for this observation: the curvature of imogolite prevents the formation of in-plane H-bonds along the directions of the nanotube circumference, lowering the enthalpy of adsorption of water molecules. The different arrangement of surface H-bonds and the resulting differences in hydration properties also affects the acidity constants of surface hydroxyl groups. This ‘nanotube effect’ may be relevant to other nanotubular systems with high curvatures, potentially impacting their wetting properties, their colloidal stability and their affinity towards hydrophobic organic moieties.

36 MATERIALS SCIENCE↗

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

Titan’s atmosphere is composed mostly of N 2 with a small amount of CH 4 , and so, shock layers around craft entering Titan’s atmosphere will contain a variety of molecules formed from H, C, and N atoms, including the cyanogen radical CN. Sensitivity analysis has shown that the radiative heat flux predicted by computational fluid dynamics (CFD) simulations of Titan entry has up to 14% uncertainty due to the rate coefficients for collisional (de)excitation reactions that control the population of CN in its first and second excited states. The red and violet emission bands from CN’s first and second excited states, respectively, are known to be large sources of radiative heat flux on capsules entering Titan’s atmosphere.[2, 3] So, the simulated population of CN in its first and second excited states is very important, but currently has some inherent uncertainty coming from the data for the rate coefficients for reaction 1. The goal of the present project is to provide improved rate coefficient data for these reactions from first principles quantum chemistry calculations. This work reports on preliminary electronic structure calculations generated at a large number of triatomic geometries of interest, which show multiple avoided crossings at collinear arrangements. This suggests that collisional (de)excitation of CN by N atoms is likely to proceed through these geometries.

Eric Geistfeld↗

The Li + CaF → Ca + LiF chemical reaction under cold conditions

The calcium monofluoride (CaF) molecule has emerged as a promising candidate for precision measurements, quantum simulation, and ultracold chemistry experiments. Inelastic and reactive collisions of laser cooled CaF molecules in optical tweezers have recently been reported and collisions of cold Li atoms with CaF are of current experimental interest. In this paper, we report ab initio electronic structure and full-dimensional quantum dynamical calculations of the Li + CaF → LiF + Ca chemical reaction. The electronic structure calculations are performed using the internally contracted multi-reference configuration-interaction method with Davidson correction (MRCI + Q). An analytic fit of the interaction energies is obtained using a many-body expansion method. Furthermore, a coupled-channel quantum reactive scattering approach implemented in hyperspherical coordinates is adopted for the scattering calculations under cold conditions. Results show that the Li + CaF reaction populates several low-lying vibrational levels and many rotational levels of the product LiF molecule and that the reaction is inefficient in the 1–100 mK regime allowing sympathetic cooling of CaF by collisions with cold Li atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A spherical crystal diffraction imager for Sandia’s Z Pulsed Power Facility

Sandia’s Z Pulsed Power Facility is able to dynamically compress matter to extreme states with exceptional uniformity, duration, and size, which are ideal for investigating fundamental material properties of high energy density conditions. X-ray diffraction (XRD) is a key atomic scale probe since it provides direct observation of the compression and strain of the crystal lattice and is used to detect, identify, and quantify phase transitions. Because of the destructive nature of Z-Dynamic Material Property (DMP) experiments and low signal vs background emission levels of XRD, it is very challenging to detect a diffraction signal close to the Z-DMP load and to recover the data. We have developed a new Spherical Crystal Diffraction Imager (SCDI) diagnostic to relay and image the diffracted x-ray pattern away from the load debris field. The SCDI diagnostic utilizes the Z-Beamlet laser to generate 6.2-keV Mn–He α x rays to probe a shock-compressed material on the Z-DMP load. Finally, a spherically bent crystal composed of highly oriented pyrolytic graphite is used to collect and focus the diffracted x rays into a 1-in. thick tungsten housing, where an image plate is used to record the data.

47 OTHER INSTRUMENTATION↗

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structural distortion and electron redistribution in dual-emitting gold nanoclusters

Deciphering the complicated excited-state process is critical for the development of luminescent materials with controllable emissions in different applications. Here we report the emergence of a photo-induced structural distortion accompanied by an electron redistribution in a series of gold nanoclusters. Such unexpected slow process of excited-state transformation results in near-infrared dual emission with extended photoluminescent lifetime. We demonstrate that this dual emission exhibits highly sensitive and ratiometric response to solvent polarity, viscosity, temperature and pressure. Thus, a versatile luminescent nano-sensor for multiple environmental parameters is developed based on this strategy. Furthermore, we fully unravel the atomic-scale structural origin of this unexpected excited-state transformation, and demonstrate control over the transition dynamics by tailoring the bi-tetrahedral core structures of gold nanoclusters. Overall, this work provides a substantial advance in the excited-state physical chemistry of luminescent nanoclusters and a general strategy for the rational design of next-generation nano-probes, sensors and switches.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Time-resolved Brownian tomography of single nanocrystals in liquid during oxidative etching

Colloidal nanocrystals inherently undergo structural changes during chemical reactions. The robust structure-property relationships, originating from their nanoscale dimensions, underscore the significance of comprehending the dynamic structural behavior of nanocrystals in reactive chemical media. Moreover, the complexity and heterogeneity inherent in their atomic structures require tracking of structural transitions in individual nanocrystals at three-dimensional (3D) atomic resolution. In this study, we introduce the method of time-resolved Brownian tomography to investigate the temporal evolution of the 3D atomic structures of individual nanocrystals in solution. The methodology is applied to examine the atomic-level structural transformations of Pt nanocrystals during oxidative etching. The time-resolved 3D atomic maps reveal the structural evolution of dissolving Pt nanocrystals, transitioning from a crystalline to a disordered structure. Our study demonstrates the emergence of a phase at the nanometer length scale that has received less attention in bulk thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomistic measurement and modeling of intrinsic fracture toughness of two-dimensional materials

Quantifying the intrinsic mechanical properties of two-dimensional (2D) materials is essential to predict the long-term reliability of materials and systems in emerging applications ranging from energy to health to next-generation sensors and electronics. Currently, measurements of fracture toughness and identification of associated atomistic mechanisms remain challenging. Herein, we report an integrated experimental–computational framework in which in-situ high-resolution transmission electron microscopy (HRTEM) measurements of the intrinsic fracture energy of monolayer MoS 2 and MoSe 2 are in good agreement with atomistic model predictions based on an accurately parameterized interatomic potential. Changes in crystalline structures at the crack tip and crack edges, as observed in in-situ HRTEM crack extension tests, are properly predicted. Such a good agreement is the result of including large deformation pathways and phase transitions in the parameterization of the inter-atomic potential. The established framework emerges as a robust approach to determine the predictive capabilities of molecular dynamics models employed in the screening of 2D materials, in the spirit of the materials genome initiative. Moreover, it enables device-level predictions with superior accuracy (e.g., fatigue lifetime predictions of electro- and opto-electronic nanodevices).

Zhang, Xu↗