Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular simulation”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Neural network potential for Zr–Rh system by machine learning

Zr–Rh metallic glass has enabled its many applications in vehicle parts, sports equipment and so on due to its outstanding performance in mechanical property, but the knowledge of the microstructure determining the superb mechanical property remains yet insufficient. Here, we develop a deep neural network potential of Zr–Rh system by using machine learning, which breaks the dilemma between the accuracy and efficiency in molecular dynamics simulations, and greatly improves the simulation scale in both space and time. Here, the results show that the structural features obtained from the neural network method are in good agreement with the cases in ab initio molecular dynamics simulations. Furthermore, we build a large model of 5400 atoms to explore the influences of simulated size and cooling rate on the melt-quenching process of Zr 77 Rh 23 . Our study lays a foundation for exploring the complex structures in amorphous Zr 77 Rh 23 , which is of great significance for the design and practical application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Investigating Morphology and Diffusion in Simulations of Precise Anion-Conducting Polymers

Using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new class of polymers which are functionalized with quaternary ammonium groups for use as anion exchange membranes. The polymers are precision polyolefins with either a trimethylammonium (p5CNMe3) or a dimethyl-hexyl ammonium (p5CNMe2Hx) pendant group at every fifth carbon along a polyethylene backbone. Simulations are performed at hydration levels of 5, 10, 15, and 20 water molecules per ammonium group. The hydrated polymers form nanoscale, percolated hydrophilic domains (water channels) in the hydrophobic polymer matrix that become wider with increasing water content. Water and hydroxide anion diffusion coefficients also increase with increasing water content. The morphology of the water domains is similar in both polymers, while the diffusion coefficients are somewhat lower in p5CNMe2Hx at fixed water content. Furthermore, the diffusion coefficients in both polymers fall on the same curve as a function of the fractal dimension of the percolated water channels, which appears to be a useful scalar measure of the effects of the nanoscale morphology on water and hydroxide anion transport.

Anion exchange membrane↗
NASA NTRS2003

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗
DOE OSTI2020

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. Finally, this parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Correlation Function Approach for Diffusion in Confined Geometries

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nano channel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

Nanoscale flow, Transport, molecular simulation↗
NASA NTRSDate not supplied

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^8 determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

J Wayne Mullinax↗
NASA NTRSDate not supplied

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^(8) determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

Quantum Computing↗
NASA NTRSDate not supplied

Efficient Implementation for Unitary Coupled Cluster State Preparation for Near-Term Quantum Computers

Unitary coupled cluster theory (UCC) is a common wave function ansatz for quantum simulation of molecular electronic structure using the variational quantum eigenvalue solver (VQE). Even for small molecules using a double-ζ basis, the number of variational parameters required to minimize the electronic energy (i.e., optimize the circuit) is large and beyond the reach of current quantum computers. For example, a circuit simulating C2 using the UCCSD ansatz and the cc-pVDZ basis set with frozen-core will require over 10,000 variational parameters and a Hilbert space of over 10^(8) determinants. To make progress on simulating such molecular systems on near-term quantum computers, we explore how much of the optimization can be approximately prepared with classical simulation while reducing the number of optimization steps performed on a quantum device. Recently, Chen, Cheng, and Freericks [J. Chem. Theory Comput. 2021, 17, 841-847] presented an algorithm for the factorized form of the UCC ansatz that allows for efficient UCC optimizations on classical hardware. We flip the algorithm around and use it to prepare approximate quantum circuits for systems that require a large number of qubits to represent. We will present results from our implementation and discuss strategies for incorporating this implementation for algorithms involving near-term quantum computers.

Quantum Computing↗
NASA NTRS2016

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding supercooling mechanism in sodium sulfate decahydrate phase-change material

Salt hydrate-based phase-change materials are considered promising for future heat storage applications in residential heating/cooling systems. Smooth phase transition from the liquid to solid phase and vice versa is essential for effective heat exchanger; however, supercooling in salt hydrates delays the onset of liquid–solid phase transition. We investigate the molecular level mechanism of supercooling in sodium sulfate decahydrate (SSD). SSD is a complex salt hydrate whose properties are governed by electrostatic forces that include pure Coulombic interactions as well as hydrogen bonds. Experimentally, we examine the importance of a nucleator in reducing supercooling temperatures. We investigated the effect of various mass concentrations of a borax nucleator on a decrease of supercooling temperatures. Molecular dynamics simulation techniques are used to obtain a basic understanding of supercooling in SSD. We observe that by introducing borax as a nucleator, there is a decrease in the supercooling temperature before nucleation. Our molecular dynamics simulations show that long-range electrostatics between sodium and sulfate ion pairs and that with polar water molecules is responsible for delayed nucleation in SSD that results in supercooling, and also, dynamics of charged molecules slows down. The lack of crystallization leads to amorphous structures in supercooled SSD.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗
DOE OSTI2025

Keeping LAMMPS cutting edge

Since its inception 30 years ago, LAMMPS has grown to be a world-class molecular dynamics code and a cornerstone of computational materials science research. This project aimed to keep LAMMPS at the forefront of molecular dynamics simulations by adapting LAMMPS to the latest developments in machine learning technology and hardware. Initially, the project set out to provide a unified implementation of active learning for efficient training data generation in LAMMPS, but the research trajectory pivoted to address more immediate and impactful opportunities. On the hardware side, recent record-breaking molecular dynamics simulations were developed on the Cerebras wafer-scale AI chip, and this project has developed an interface between LAMMPS and the hardware-specific molecular dynamics code to accelerate and simplify development and user adoption. On the software side, PyTorch’s Ahead-of-Time (AOT) compilation features promised increased performance for state-of-the-art equivariant neural network potentials, and this project laid the groundwork for their adoption in LAMMPS, resulting in a nearly 20x acceleration in extreme cases. Combined with a comprehensive benchmark study of LAMMPS across all current exascale systems, this project has reinforced LAMMPS’s role as a versatile, high-performance tool for current and future materials science applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Molecular Lignin Solubility and Structure in Organic Solvents

Lignin, a polymer found in the secondary plant cell wall of terrestrial plants, is the single largest source of renewable aromatics and has attracted considerable attention as a feedstock for potential industrial use. However, the secondary plant cell wall is a crowded environment, and lignin in its native form interacts with other biomass components within a larger network. Application of some organic solvents is known to liberate lignin from this network and creates lignin-rich streams suitable for conversion into target products. Through molecular-scale lignin simulation, we analyze how diverse lignin polymers change their structure in response to varying organic solvent environments. We quantify the relationship between solvent polarity and lignin polymer extension, observing maximum polymer expansion and solvation for solvents with polarity near those of dimethyl sulfoxide. From our observations at the nanoscale, increasing syringyl content within lignin polymers reduces the expansion of the polymer in organic solvent environments and decreases the free energy difference compared to aqueous solvent environments, thereby reducing solubility for high syringyl lignin polymers. The conformational transition rates between lignin polymer shapes increased through a combination of the solvent diffusion constant and polymer extension. The molecular simulations indicate that there is likely no single optimal organic solvent for lignin. Different solvent mixtures have optimal or near-optimal properties in solubilizing lignin polymers, thereby disrupting interactions with other biopolymers.

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

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute–solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (Fe ii Ru iii ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Overall, our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120–200 cm -1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

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

Molecular origin of negative lithium transference in electrolytes with star-shaped multivalent anions

Large-scale molecular dynamics simulations illustrate that highly correlated cation–anion motion leads to negative t 0+ on the order of 1 in lithium electrolytes with star-shaped multivalent anions. Large multivalent anions have gained increasing attention for their potential to improve lithium transference in electrolytes. We employ large-scale molecular dynamics simulations based on the Onsager transport framework to investigate ion transport in a lithium electrolyte with star-shaped multivalent anions. The simulations show that t 0+, the cation transference number with respect to solvent velocity, is negative over a wide range of concentration. This is consistent with experimental data reported previously. The simulation-based Onsager transport coefficients reveal that the magnitudes of the cation–cation, anion–anion, and cation–anion correlations are comparable, a signature of highly correlated motion in the electrolyte. Examination of the cation solvation environment indicates the presence of strong cation–anion association across the entire concentration range, which leads to negative t 0+ on the order of −1. Both simulation and experiment also show that the maximum value of t 0+ reaches 0 when the cation concentration is c + = 0.4 M. This is the concentration at which the anions begin to spatially overlap, and lithium ions serve as dynamic linkers to balance cation–cation and cation–anion correlations. Our results provide molecular-level insights into the origin of transference in multivalent electrolytes.

Fang, Chao↗