Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Recovery of MOF-5 from Extreme High-Pressure Conditions Facilitated by a Modern Pressure Transmitting Medium

Here mechanisms underlying the mechanically induced amorphization of metal-organic frameworks (MOFs) are of current interest, and both high-pressure experimentation and molecular dynamics simulations have been used to reveal the fundamentals of load bearing, deformation, and pressure-induced amorphization (PIA) in these highly porous materials. Unfortunately, MOFs are typically highly susceptibility to amorphization, which limits the conditions under which they can be processed and used. However, their flexible structures can be stabilized to high pressures by incorporating guest species into the framework matrix. In this study, a large molecule pressure transmitting medium (DAPHNE 7575 ) is used as a structure-fortifying guest species to stabilize the prototypical MOF -5 to ultra-high pressures ( >9 GPa) and enable the recovery of crystalline material upon decompression. Structural changes associated with the penetration of the pressure transmitting medium on compression are examined using a combination of high pressure synchrotron powder diffraction and molecular dynamics simulations. This work enhances the understanding of PIA in MOFs while showcasing a potential route for the stabilization of MOFs to surprisingly high-pressures.

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

Temperature and composition dependent screw dislocation mobility in austenitic stainless steels from large-scale molecular dynamics

Abstract Extensive molecular dynamics simulations are performed to determine screw dislocation mobility in austenitic Fe 0.7 Ni x Cr 0.3-x stainless steels as a function of temperature ranging from 100 to 1300 K, resolved shear stress from 30 to 140 MPa, and Ni composition from 0.0 to 30.0 at%. These mobility data are fitted to a linear mobility law with a nonzero stress offset, referred to as the threshold stress. We find that both the linear drag coefficient and the threshold stress increase with Ni composition. The drag coefficient increases with temperature, whereas the threshold stress decreases with temperature. Based on these calculations, we determine fitting functions for the linear solute drag coefficient as a function of temperature and composition. The mobility laws determined in this study may serve to inform dislocation dynamics simulations pertinent to dislocation network evolution at elevated temperatures for a wide composition range of austenitic stainless steels.

Chu, Kevin (ORCID:000000031929710X)↗
DOE OSTI2017

Influence of metal ions intercalation on the vibrational dynamics of water confined between MXene layers

Two-dimensional carbides and nitrides of early transition metals (MXenes) combine high conductivity with hydrophilic surfaces, which make them promising for energy storage, electrocatalysis, and water desalination. Effects of intercalated metal ions on the vibrational states of water confined in Ti 3 C 2 T x MXenes have been explored using inelastic neutron scattering (INS) and molecular dynamics simulations to better understand the mechanisms that control MXenes’ behavior in aqueous electrolytes, water purification and other important applications. Here, we observe INS signal from water in all samples, pristine and with lithium, sodium or potassium ions intercalated between the 2D Ti 3 C 2 T x layers. However, only a small amount of water is found to reside in Ti 3 C 2 T x intercalated with metal ions. Water in pristine Ti 3 C 2 T x is more disordered, with bulk-like characteristics, in contrast to intercalated Ti 3 C 2 T x , where water is more ordered, irrespective of the metal ions used for intercalation. The ordering of the confined water increases with the ion size. Lastly, this finding is further confirmed from molecular dynamics simulation which showed an increase in interference of water molecules with increasing ion size resulting in a concomitant decrease in water mobility, therefore, providing a guidance to tailor MXene properties for energy and environmental applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dueling Backbones: Comparing Peptoid and Peptide Analogues of a Mussel Adhesive Protein

Ensembles of amino acid side chains often dominate the interfacial interactions of intrinsically disordered proteins; however, backbone contributions are far from negligible. Using a combination of nanoscale force measurements and molecular dynamics simulations, we demonstrated with analogous mussel-mimetic adhesive peptides and peptoids both 34 residues long that highly divergent adhesive/cohesive outcomes can be achieved on mica surfaces by altering backbone chemistry only. The Phe, Tyr, and Dopa containing peptoid variants used in this study deposited as dehydrated and incompressible films that facilitated analysis of peptoid side chain contributions to adhesion and cohesion. For example, whereas Phe and Dopa peptoids exhibited similar cohesion, Dopa peptoids were ~3 times more adhesive than Phe peptoids on mica. Compared with the peptides, Phe peptoid achieved only ~20% of Phe containing peptide adhesion, but the Dopa peptoids were >2-fold more adhesive than the Dopa peptides. Cation-π interactions accounted for some but not all of the cohesive interactions. Finally, our results were corroborated by molecular dynamics simulations and highlight the importance of backbone chemistry and the potential of peptoids or peptoid/peptide hybrids as wet adhesives and primers.

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

Imaging covalent bond formation by H atom scattering from graphene

Viewing the atomic-scale motion and energy dissipation pathways involved in forming a covalent bond is a longstanding challenge for chemistry. We performed scattering experiments of H atoms from graphene and observed a bimodal translational energy loss distribution. Using accurate first-principles dynamics simulations, we show that the quasi-elastic channel involves scattering through the physisorption well where collision sites are near the centers of the six-membered C-rings. The second channel results from transient C–H bond formation, where H atoms lose 1 to 2 electron volts of energy within a 10-femtosecond interaction time. This remarkably rapid form of intramolecular vibrational relaxation results from the C atom’s rehybridization during bond formation and is responsible for an unexpectedly high sticking probability of H on graphene.

Science & Technology - Other Topics↗
DOE OSTI2016

Theory of multiple quantum dot formation in strained-layer heteroepitaxy

We develop a theory for the experimentally observed formation of multiple quantum dots (QDs) in strained-layer heteroepitaxy based on surface morphological stability analysis of a coherently strained epitaxial thin film on a crystalline substrate. Using a fully nonlinear model of surface morphological evolution that accounts for a wetting potential contribution to the epitaxial film's free energy as well as surface diffusional anisotropy, we demonstrate the formation of multiple QD patterns in self-consistent dynamical simulations of the evolution of the epitaxial film surface perturbed from its planar state. The simulation predictions are supported by weakly nonlinear analysis of the epitaxial film surface morphological stability. We find that, in addition to the Stranski-Krastanow instability, long-wavelength perturbations from the planar film surface morphology can trigger a nonlinear instability, resulting in the splitting of a single QD into multiple QDs of smaller sizes, and predict the critical wavelength of the film surface perturbation for the onset of the nonlinear tip-splitting instability. The theory provides a fundamental interpretation for the observations of “QD pairs” or “double QDs” and other multiple QDs reported in experimental studies of epitaxial growth of semiconductor strained layers and sets the stage for precise engineering of tunable-size nanoscale surface features in strained-layer heteroepitaxy by exploiting film surface nonlinear, pattern forming phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

In Situ Real-Time Radiographic Study of Thin Film Formation Inside Rotating Hollow Spheres

Hollow spheres with uniform coatings on the inner surface have applications in optical devices, time- or site controlled drug release, heat storage devices, and target fabrication for inertial confinement fusion experiments. The fabrication of uniform coatings, which is often critical for the application performance, requires precise understanding and control over the coating process and its parameters. Here, we report on in-situ real-time radiography experiments that provide critical spatio-temporal information about the distribution of fluids inside hollow spheres during uniaxial rotation. Image analysis and computer fluid dynamics simulations were used to explore the effect of liquid viscosity and rotational velocity on the film uniformity. The data were then used to demonstrate the fabrication of uniform sol-gel chemistry derived porous polymer films inside 2mm inner diameter diamond shells.

aerogel↗
DOE OSTI2017

NetFlow Dynamics

NetFlow Dynamics is a web-accessible analysis environment for simulating dynamic flows of materials on model networks. Performing a simulation requires both the NetFlow Dynamics application and a network model which is a description of the structure of the nodes and edges of a network including the flow capacity of each edge and the storage capacity of each node, and the sources and sinks of the material flowing on the network. NetFlow Dynamics consists of databases for storing network models, algorithms to calculate flows on networks, and a GIS-based graphical interface for performing simulations and viewing simulation results. Simulated flows are dynamic in the sense that flows on each edge of the network and inventories at each node change with time and can be out of equilibrium with boundary conditions. Any number of network models could be simulated using Net Flow Dynamics. To date, the models simulated have been models of petroleum infrastructure. The main model has been the National Transportation Fuels Model (NTFM), a network of U.S. oil fields, transmission pipelines, rail lines, refineries, tank farms, and distribution terminals. NetFlow Dynamics supports two different flow algorithms, the Gradient Flow algorithm and the Inventory Control algorithm, that were developed specifically for the NetFlow Dynamics application. The intent is to add additional algorithms in the future as needed. The ability to select from multiple algorithms is desirable because a single algorithm never covers all analysis needs. The current algorithms use a demand-driven capacity-constrained formulation which means that the algorithms strive to use all available capacity and stored inventory to meet desired flows to sinks, subject to the capacity constraints of each network component. The current flow algorithms are best suited for problems in which a material flows on a capacity-constrained network representing a supply chain in which the material supplied can be stored at each node of the network. In the petroleum models, the flowing materials are crude oil and refined products that can be stored at tank farms, refineries, or terminals (i.e. the nodes of the network). Examples of other network models that could be simulated are currency flowing in a financial network, agricultural products moving to market, or natural gas flowing on a pipeline network.

Corbet Jr., Thomas F↗
DOE OSTI2018

Development of CHARMM-Compatible Force-Field Parameters for Cobalamin and Related Cofactors from Quantum Mechanical Calculations

Corrinoid cofactors such as cobalamin are used by many enzymes and are essential for most living organisms. Therefore, there is broad interest in investigating cobalamin–protein interactions with molecular dynamics simulations. Previously developed parameters for cobalamins are based mainly on crystal structure data. Here, we report CHARMM-compatible force field parameters for several corrinoids developed from quantum mechanical calculations. We provide parameters for corrinoids in three oxidation states, Co 3+ , Co 2+ , and Co 1+ , and with various axial ligands. Lennard-Jones parameters for the cobalt center in the Co(II) and Co(I) states were optimized using a helium atom probe, and partial atomic charges were obtained with a combination of natural population analysis (NPA) and restrained electrostatic potential (RESP) fitting approaches. The Force Field Toolkit was used to optimize all bonded terms. The resulting parameters, determined solely from calculations of cobalamin alone or in water, were then validated by assessing their agreement with density functional theory geometries and by analyzing molecular dynamics simulation trajectories of several corrinoid proteins for which X-ray crystal structures are available. In each case, we obtained excellent agreement with the reference data. In comparison to previous CHARMM-compatible parameters for cobalamin, we observe a better agreement for the fold angle and lower RMSD in the cobalamin binding site. In conclusion, the approach described here is readily adaptable for developing CHARMM-compatible force-field parameters for other corrinoids or large biomolecules.

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

Mechanical behavior of ultralight nickel metamaterial

The mechanical response of ultralight kagomé structures consisting of hollow nickel (Ni) nanotubes and solid Ni nanorods to compression is studied using molecular dynamics simulations. In both kagomé architectures, 16[112] Shockley partial dislocations and twin formation are observed under compression. The structure made from solid nanorods shows deformation near both the nodes and beams of the kagomé lattice. The hollow kagomé architecture has a higher yield point than the solid kagomé structure. The deformation in the hollow nanotube structure is mostly localized in the nodal region for strains less than 11%. At higher strains, the deformation sets in all the struts and nodes of the hollow kagomé lattice. Owing to this two-stage deformation mechanism, the hollow Ni nanotube kagomé structure shows less bending and greater toughness than the solid Ni nanorod kagomé architecture.

Rajak, Pankaj (ORCID:0000000263446056)↗
DOE OSTI2020

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Raman studies of hydrogen trapped in As 4 O 6 ·2H 2 at high pressure and low temperature

In this work, Raman spectroscopic measurements of the arsenolite–hydrogen inclusion compound As 4 O 6 ·2H 2 were performed in diamond anvil cells at high pressure and variable temperature down to 80 K. The experimental results were complemented by ab initio molecular dynamics simulations and phonon calculations. Observation of three hydrogen vibrons in As 4 O 6 ·2H 2 is reported in the entire temperature and pressure range studied (up to 24 GPa). While the experiments performed with protium and deuterium at variable temperatures allowed for the assignment of two vibrons as Q 1 (1) and Q 1 (0) transitions of ortho and para spin isomers of hydrogen trapped in the inclusion compound, the origin of the third vibron could not be unequivocally established. Low-temperature spectra revealed that the lowest-frequency vibron is actually composed of two overlapping bands of A g and T 2g symmetries dominated by H 2 stretching modes as predicted by our previous density functional theory calculations. We observed low-frequency modes of As 4 O 6 ·2H 2 vibrations dominated by H 2 “librations,” which were missed in a previous study. A low-temperature fine structure was observed for the J = 0 → 2 and J = 1 → 3 manifolds of hydrogen trapped in As 4 O 6 ·2H 2 , indicating the lifting of degeneracy due to an anisotropic environment. A non-spherical distribution was captured by molecular dynamics simulations, which revealed that the trajectory of H 2 molecules is skewed along the crystallographic $\langle$111$\rangle$ direction. Last but not least, low-temperature synchrotron powder x-ray diffraction measurements on As 4 O 6 ·2H 2 revealed that the bulk structure of the compound is preserved down to 5 K at 1.6 GPa.

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

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions. Our scheme is based on quadratic response theory and includes orbital relaxation terms neglected in previous approaches. Simultaneous computation of multiple derivative couplings and energy gradients enables efficient multistate nonadiabatic molecular dynamics simulations in conjunction with Tully's fewest switches surface hopping (SH) method. We benchmark the thus obtained multistate TDDFT-SH scheme by simulating ultrafast decay of UV-photoexcited thymine, for which accurate gas-phase data from ultrafast spectroscopy experiments are available. The calculations predict a fast 153 fs decay from the bright S 2 to the dark S 1 excited state, followed by a much slower 14 ps S 1 deactivation to the ground state; statistical uncertainties were estimated using bootstrap sampling. These results agree well with the experimentally observed time constants of 100–200 fs and 5–7 ps, respectively, unlike previous multiconfigurational self-consistent field and second-order algebraic diagrammatic construction calculations. Furthermore, our results support the S 1 -trapping hypothesis. For thymine, the computational cost of a single TDDFT-SH time-step including the lowest 3 states, all couplings and gradients, is ~5 times larger than the cost of a single Born–Oppenheimer dynamics time step for the ground state in our implementation. Thus, ps nonadiabatic dynamics simulations using multistate hybrid TDDFT-SH for systems with up to ~100 atoms are possible without drastic approximations on single workstation nodes. Our implementation will be made available through Turbomole.

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

Challenges in simulating beam dynamics of dielectric laser acceleration

Dielectric Laser Acceleration (DLA) achieves the highest gradients among structure-based electron accelerators. The use of dielectrics increases the breakdown field limit, and thus the achievable gradient, by a factor of at least 10 in comparison to metals. Experimental demonstrations of DLA in 2013 led to the Accelerator on a Chip International Program (ACHIP), funded by the Gordon and Betty Moore Foundation. In ACHIP, our main goal is to build an accelerator on a silicon chip, which can accelerate electrons from below 100 keV to above 1 MeV with a gradient of at least 100 MeV/m. For stable acceleration on the chip, magnet-only focusing techniques are insufficient to compensate the strong acceleration defocusing. Thus, spatial harmonic and Alternating Phase Focusing (APF) laser-based focusing techniques have been developed. Additionally, we have developed the simplified symplectic tracking code DLAtrack6D, which makes use of the periodicity and applies only one kick per DLA cell, which is calculated by the Fourier coefficient of the synchronous spatial harmonic. Due to coupling, the Fourier coefficients of neighboring cells are not entirely independent and a field flatness optimization (similarly as in multi-cell cavities) needs to be performed. The simulation of the entire accelerator on a chip by a Particle In Cell (PIC) code is possible, but impractical for optimization purposes. Finally, we have also outlined the treatment of wake field effects in attosecond bunches in the grating within DLAtrack6D, where the wake function is computed by an external solver.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Frequency–wavenumber spectral analysis of spatio-temporal flows

We propose a fully spatio-temporal approach for identifying spatially varying modes of oscillation in fluid dynamics simulation output by means of multitaper frequency–wavenumber spectral analysis. One-dimensional spectrum estimation has proven to be a valuable tool in the analysis of turbulence data applied spatially to determine the rate of energy transport between spatial scales, or temporally to determine frequencies of oscillatory flows. It also allows for the quantitative comparison of flow characteristics between two scenarios using a standard basis. It has the limitation, however, that it neglects coupling between spatial and temporal structures. Two-dimensional frequency–wavenumber spectral analysis allows one to decompose waveforms into standing or travelling variety. The extended higher-dimensional multitaper method proposed here is shown to have improved statistical properties over conventional non-parametric spectral estimators, and is accompanied by confidence intervals which estimate their uncertainty. Multitaper frequency–wavenumber analysis is applied to a canonical benchmark problem, namely, a direct numerical simulation of von Kármán vortex shedding off a square wall-mounted cylinder with two inflow scenarios with matching momentum-thickness Reynolds numbers $Re_{\unicode[STIX]{x1D703}}\approx 1000$ at the obstacle. Frequency–wavenumber analysis of a two-dimensional section of these data reveals that although both the laminar and turbulent inflow scenarios show a turbulent $-5/3$ cascade in wavenumber ( $\unicode[STIX]{x1D708}$ ) and frequency ( $f$ ), the flow characteristics differ in that there is a significantly more prominent discrete harmonic oscillation near $(f,\unicode[STIX]{x1D708})=(0.2,0.21)$ in wavenumber and frequency in the laminar inflow scenario than the turbulent scenario. This frequency–wavenumber pair corresponds to a travelling wave with velocity near one near the centre path of the vortex street.

Geoga, Christopher J. (ORCID:0000000250604799)↗
DOE OSTI2016

Atomistic Conversion Reaction Mechanism of WO 3 in Secondary Ion Batteries of Li, Na, and Ca

Abstract Intercalation and conversion are two fundamental chemical processes for battery materials in response to ion insertion. The interplay between these two chemical processes has never been directly seen and understood at atomic scale. Here, using in situ HRTEM, we captured the atomistic conversion reaction processes during Li, Na, Ca insertion into a WO 3 single crystal model electrode. An intercalation step prior to conversion is explicitly revealed at atomic scale for the first time for Li, Na, Ca. Nanoscale diffraction and ab initio molecular dynamic simulations revealed that after intercalation, the inserted ion–oxygen bond formation destabilizes the transition‐metal framework which gradually shrinks, distorts and finally collapses to an amorphous W and M x O (M=Li, Na, Ca) composite structure. This study provides a full atomistic picture of the transition from intercalation to conversion, which is of essential importance for both secondary ion batteries and electrochromic devices.

He, Yang↗