Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic systems”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Digital quantum simulation of cavity quantum electrodynamics: insights from superconducting and trapped ion quantum testbeds

We explore the potential for hybrid development of quantum hardware where currently available quantum computers simulate open cavity quantum electrodynamical (CQED) systems for applications in optical quantum communication, simulation and computing. Our simulations make use of a recent quantum algorithm that maps the dynamics of a singly excited open Tavis–Cummings model containing N atoms coupled to a lossy cavity. We report the results of executing this algorithm on two noisy intermediate-scale quantum computers: a superconducting processor and a trapped ion processor, to simulate the population dynamics of an open CQED system featuring N = 3 atoms. By applying technology-specific transpilation and error mitigation techniques, we minimize the impact of gate errors, noise, and decoherence in each hardware platform, obtaining results which agree closely with the exact solution of the system. These results can be used as a recipe for efficient and platform-specific quantum simulation of cavity–emitter systems on contemporary and future quantum computers.

cavity QED↗
NASA NTRS1972

An application of a space technique for pollution detection.

Description of the hazardous gas detection system (HGDS) used during Apollo-Saturn launch tests for identifying and analyzing extremely minute amounts of gases escaping during cryogenic loading. These gases included hydrogen, helium, nitrogen, oxygen, nitrogen tetroxide, and hydrazine N2H3 radical. The heart of the HGDS is a mass spectrometer operating under high vacuum and capable to measure both qualitatively and quantitatively the presence and percentage of any gas whose molecular weight is between 2 and 150 atomic mass units. It is suggested that such a system could be utilized to detect environmental pollutants.

Hall, J. J.↗
DOE OSTI2021

Vibrational Properties of Pristine and Lithium-Intercalated Black Phosphorous under High-Pressure

Structural evolution of Li-intercalated and pristine black phosphorous (BP) under high-pressure (up to ≈8 GPa) is studied using in situ Raman spectroscopy. Even though both materials show a monotonic blueshift of the out-of-plane vibrational mode (A 1 g ) with pressure, Li-intercalated BP do not show a blueshift until a threshold pressure (2.4 GPa) is reached to compensate the structural expansion caused by intercalation. However, the in-plane modes (B 2g and A 2 g) in each sample respond differently. In the mid-pressure region, they both show redshifts which in Li-intercalated BP is also followed by abrupt blueshifts. Such behavior indicates pressure-induced structural reorganizations inside the material. Computational modeling reveals the existence of a process of P-P bond breaking and reforming in the system due to the redistribution of intercalated Li atoms under pressure. This work shows the significance of combined effect of pressure and intercalation on structural changes in the search for new phases of BP and other two-dimensional (2D) materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the integration of molecular dynamics, data science, and experiments for studying solvent effects on catalysis

Computational workflows that combine molecular dynamics (MD) simulations and emerging data-centric (DC) methods can accelerate the screening and analysis of solvent systems experimentally and computationally. Here, MD simulations provide atomic positions and velocities of reactant, solvent, and catalyst materials that can be manipulated into data representations that in turn can be used by DC techniques to conduct predictive modeling, feature extraction, and experimental design. For liquid-phase catalytic applications, emerging DC techniques such as Convolutional and Graph Neural Networks (CNN/GNN), Topological Data Analysis (TDA), and Active Learning (AL) can leverage MD and experimental data to quickly predict solvent effects on reaction outcomes. For instance, in recent studies, 3D solvent environments obtained with MD have been exploited by CNNs to predict experimental reaction rates for homogeneous acid-catalyzed lignocellulosic processes. In this perspective, we discuss basic principles of DC methods and how these can be combined with MD to enable high-throughput screening of solvent selection for diverse catalysis applications.

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

Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains

Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monte Carlo simulations for understanding the transport properties of metallic nanowires

We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

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

Intelligent Agents for the Optimization of Atomic Layer Deposition

Atomic layer deposition (ALD) is a highly controllable thin film synthesis approach with applications in computing, energy, and separations. The flexibility of ALD means that it can access a massive chemical catalogue; however, this chemical and process diversity results in significant challenges in determining processing parameters that result in stable and uniform film growth with minimal precursor consumption. In situ measurements of the ALD growth per cycle (GPC) can accelerate process development but it still requires expert intuition and time-consuming trial and error to identify acceptable processing parameters. This procedure is made more difficult by the presence of experimental noise in the GPC values and the complexity of ALD surface chemistries. Here, a need exists for efficient optimization approaches capable of autonomously determining processing conditions resulting in optimal ALD film growth. In this work, we present the development of three optimization strategies and compare their performance in optimizing four simulated ALD processes. Furthermore, the effect of noise in the GPC measurements on optimization convergence is studied.

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

How Precisely Can Individual Molecules Be Analyzed? A Case Study on Locally Quantifying Forces and Energies Using Scanning Probe Microscopy

Recent advances in scanning probe microscopy methodology have enabled the measurement of tip-sample interactions with picometer accuracy in all three spatial dimensions, thereby providing a detailed site-specific and distance-dependent picture of the related properties. This paper explores the degree of detail and accuracy that can be achieved in locally quantifying probe-molecule interaction forces and energies for adsorbed molecules. Toward this end, cobalt phthalocyanine (CoPc), a promising CO 2 reduction catalyst, was studied on Ag(111) as a model system using low-temperature, ultrahigh vacuum noncontact atomic force microscopy. Data were recorded as a function of distance from the surface, from which detailed three-dimensional maps of the molecule's interaction with the tip for normal and lateral forces as well as the tip-molecule interaction potential were constructed. The data were collected with a CO molecule at the tip apex, which enabled a detailed visualization of the atomic structure. Determination of the tip-substrate interaction as a function of distance allowed isolation of the molecule-tip interactions; when analyzing these in terms of a Lennard-Jones-type potential, the atomically resolved equilibrium interaction energies between the CO tethered to the tip and the CoPc molecule could be recovered. Interaction energies peaked at less than 160 meV, indicating a physisorption interaction. As expected, the interaction was weakest at the aromatic hydrogens around the periphery of the molecule and strongest surrounding the metal center. In conclusion, the interaction, however, did not peak directly above the Co atom but rather in pockets surrounding it.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice thermal transport in two-dimensional alloys and fractal heterostructures

Abstract Engineering thermal transport in two dimensional materials, alloys and heterostructures is critical for the design of next-generation flexible optoelectronic and energy harvesting devices. Direct experimental characterization of lattice thermal conductivity in these ultra-thin systems is challenging and the impact of dopant atoms and hetero-phase interfaces, introduced unintentionally during synthesis or as part of deliberate material design, on thermal transport properties is not understood. Here, we use non-equilibrium molecular dynamics simulations to calculate lattice thermal conductivity of $${\mathrm {(Mo|W)Se_2}}$$ ( Mo | W ) Se 2 monolayer crystals including $${\mathrm {Mo}}_{1-x}{\mathrm {W}}_x{\mathrm {Se_2}}$$ Mo 1 - x W x Se 2 alloys with substitutional point defects, periodic $${\mathrm {MoSe_2}|\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures with characteristic length scales and scale-free fractal $${\mathrm {MoSe_2}}|{\mathrm {WSe_2}}$$ MoSe 2 | WSe 2 heterostructures. Each of these features has a distinct effect on phonon propagation in the crystal, which can be used to design fractal and periodic alloy structures with highly tunable thermal conductivities. This control over lattice thermal conductivity will enable applications ranging from thermal barriers to thermoelectrics.

42 ENGINEERING↗
DOE OSTI2023

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. Here, in order to gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li + , Na + and K + . The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

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

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2021

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

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

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Driving and characterizing nucleation of urea and glycine polymorphs in water

Crystal nucleation is relevant across the domains of fundamental and applied sciences. However, in many cases, its mechanism remains unclear due to a lack of temporal or spatial resolution. To gain insights into the molecular details of nucleation, some form of molecular dynamics simulations is typically performed; these simulations, in turn, are limited by their ability to run long enough to sample the nucleation event thoroughly. To overcome the timescale limits in typical molecular dynamics simulations in a manner free of prior human bias, here, we employ the machine learning-augmented molecular dynamics framework “reweighted autoencoded variational Bayes for enhanced sampling (RAVE).” We study two molecular systems—urea and glycine—in explicit all-atom water, due to their enrichment in polymorphic structures and common utility in commercial applications. From our simulations, we observe multiple back-and-forth nucleation events of different polymorphs from homogeneous solution; from these trajectories, we calculate the relative ranking of finite-sized polymorph crystals embedded in solution, in terms of the free-energy difference between the finite-sized crystal polymorph and the original solution state. We further observe that the obtained reaction coordinates and transitions are highly nonclassical.

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

Quantum optimal control of ten-level nuclear spin qudits in Sr 87

We study the ability to implement unitary maps on states of the I = 9 / 2 nuclear spin in Sr 87 , a d = 10 dimensional (qudecimal) Hilbert space, using quantum optimal control. Through a combination of nuclear spin resonance and a tensor ac Stark shift, by solely modulating the phase of a radio-frequency magnetic field, the system is quantum controllable. Alkaline-earth-metal atoms, such as Sr 87 , have a very favorable figure of merit for such control due to narrow intercombination lines and the large hyperfine splitting in the excited states. We numerically study the quantum speed limit, optimal parameters, and the fidelity of arbitrary state preparation and full SU(10) maps, including the presence of decoherence due to optical pumping induced by the light-shifting laser. We also study the use of robust control to mitigate some dephasing due to inhomogeneities in the light shift. We find that with an rf Rabi frequency of Ω rf and 0.5% inhomogeneity in the the light shift we can prepare an arbitrary Haar-random state in a time T = 4.5 π / Ω rf with average fidelity ( F ψ ) = 0.9992 , and an arbitrary Haar-random SU(10) map in a time T = 24 π / Ω rf with average fidelity ( F U ) = 0.9923 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗