Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular 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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Phenolic Polymer Interactions with Water and Ethylene Glycol Solvents

Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final post-cured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through ab initio techniques and molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H2O. Three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H2O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H2O, and other similar solvents.

Haskins, Justin B.↗
DOE OSTI2021

A molecular dynamics study of domain switching in BiFeO 3 nanofilm under DC electric field

The performance of ferroelectric materials is deeply rooted in the domain switching process triggered by an applied electric field. To understand the effect of controlling factors on the domain switching characteristics, nanoscale studies are essential. In this work, an atomistic-scale computational framework is developed to study the domain switching process of bismuth ferrite (BiFeO 3 ) nanofilms. Through a series of molecular dynamics simulations, the effects of duration and magnitude of an electric field, domain size and nanofilm thickness on the domain switching behavior are investigated. The evolution of domains under different magnitudes and durations of the applied electric field is studied, and the results of unswitched and switched cases are summarized in a switching map. Analysis of various nanofilms reveals that by increasing the film thickness, the switching time is increased. Additionally, the nanofilms with different domain sizes and configurations are examined, and the minimum domain size for complete reversal switching is obtained. Our results are helpful to predict and tune the switching process in ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gradient microstructure induced shear band constraint, delocalization, and delayed failure in CuZr nanoglasses

The mechanical properties of metallic glasses (MGs) can be significantly tuned by introducing a nanoscale glass-glass interface microstructure. Here, we design three types of Cu 64 Zr 36 MG seamless gradient microstructures by combining glassy grain sizes from 3 to 15 nm. We perform tensile loading molecular dynamics simulations on the named gradient nanoglasses to characterize their mechanical behavior. In all cases, results indicate the generation of gradient plasticity following the gradient in grain size, i.e. local plasticity increases from large to small grain size regions. Loading parallel to the gradient direction results in initial broad necking followed by catastrophic failure. In contrast, loading perpendicular to the gradient direction results in diffused shear band propagation from large to small grain size regions, delaying the generation of a dominant shear band and further delocalizing plastic deformation. In conclusion, the results demonstrate the synergistic effects generated by heterogeneous gradient designs in the mechanical behavior of metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

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

Decomposition of liquid-liquid extraction organic phase structure into critical and pre-peak contributions

Solution nanostructure induced by amphiphile self-association plays an important role in various chemical and physical processes, including liquid-liquid extraction (LLE). Small angle scattering techniques are an essential means of probing this structure, but, due to non-uniqueness of scattering patterns, their interpretation is not trivial. Here, in this study, we decompose small angle x-ray scattering (SAXS) patterns for a range of binary LLE organic phases into two components: Ornstein-Zernike and pre-peak contributions resulting from, respectively, critical concentration fluctuations and packing of the amphiphilic extractant molecules in the nonpolar aliphatic diluent. While we previously applied Ornstein-Zernike scattering models to similar organic phases, including the pre-peak explicitly in the fit allows us to measure weak fluctuations at high extractant concentration and simultaneously obtain information on the position and intensity of the correlation peak related to amphiphile packing. Scattering patterns and partial structure factors calculated from molecular dynamics simulations support this interpretation. We demonstrate how this minimal scattering model describes nanostructuring for a large number of extractant types over their entire extractant/diluent composition ranges, suggesting simple and universal behavior. The straightforward assignment of these structural contributions facilitates the comparison of these features between different classes of extractant molecules and will enable further studies of organic phase aggregation. Applicability to more complex, process-relevant organic phases is illustrated with post-contact organic phases containing significant quantities of extracted lanthanide nitrate salt, which we find are readily described by this model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Healing of nanocracks by collision cascades in nickel

We use molecular dynamics simulations to investigate the interaction between collision cascades and pre-existing nano-cracks in nickel (Ni). We find that collision cascades cause the cracks to heal whenever a cascade initiates a thermal spike whose core overlaps with the crack. Moreover, crack healing affects the distribution of radiation-induced defects generated during a collision cascade: defect clusters, such as dislocations and stacking fault tetrahedra, dominate in the presence of crack healing, while isolated point defects dominate in its absence. Furthermore, this work indicates that irradiation may reduce the number of pre-existing crack-like flaws in metallic components, potentially improving their mechanical resilience.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of Interfacial Structure on Heterogeneous Nucleation of Amorphous Carbonates

For this work, classical molecular dynamics simulations were performed to provide physical insight into the impact of interfacial structure on the heterogeneous nucleation of amorphous calcium carbonate (ACC, CaCO 3 ·H 2 O) and amorphous magnesium carbonate (AMC, MgCO 3 ·H 2 O) by using α-quartz as a model substrate. Interfacial structure and energies were computed for ACC and AMC in contact with the (100), (001), and (101) α-quartz surfaces. The simulations showed α-quartz surfaces drew water molecules out of the carbonate nuclei to form a partial hydration layer. The formation of a partial hydration layer and its disruption to the ACC/AMC structure meant the α-quartz–ACC/AMC interfaces were not energetically favored relative to separate α-quartz–water and ACC/AMC–water interfaces and, thus, homogeneous ACC/AMC nucleation was favored over heterogeneous nucleation. The CMD simulations hence provided an atomic-level explanation for a reported nonclassical growth mechanism whereby carbonate minerals grow via homogeneous nucleation and subsequent surface attachment of amorphous intermediates.

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

Molecular Tuning of Ether Cosolvent Chemistry for High-Voltage Sodium-Ion Batteries

Ethers as electrolyte cosolvents in sodium-ion batteries (SIBs) provide favorable Na + solvation and interfacial properties, but their low oxidative stability limits their use in high-voltage SIBs. Herein, we address this limitation via molecular tuning of ether cosolvents for high-voltage (4.2 V) hard carbon || NaNi 0.33 Fe 0.33 Mn 0.33 O 2 full cells. Tetrahydropyran (THP) is functionalized with a nitrile group to form tetrahydropyran-4-carbonitrile (THPCN). To delineate the effect of nitrile functionalization and benchmark ether against a conventional carbonate, THP, THPCN, and diethyl carbonate (DEC) are evaluated as cosolvents with ethylene carbonate. Nitrile functionalization lowers the HOMO energy of the ether, extends the electrolyte stability window, and alters Na⁺ solvation. Spectroscopic techniques and molecular dynamics simulations reveal that THPCN exhibits predominantly aggregate-dominated solvation (95.1 %) with weakened Na + -solvent interactions, producing the most anion-rich environment relative to DEC and THP cosolvents. THPCN-modified solvation promotes the formation of highly conductive, fluorine-enriched interphases that suppress parasitic reactions. Pouch full cells with THPCN sustained ~ 600 cycles at 4.2 V, outperforming THP and DEC. Operando gas analysis reveals that THPCN reduces CO 2 generation by 45% and H 2 generation by 30% relative to THP. Furthermore, the findings demonstrate nitrile functionalization as a molecular design strategy to stabilize ethers and enable high-voltage SIBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Uncertainty-driven dynamics for active learning of interatomic potentials

Machine learning (ML) models, if trained to data sets of high-fidelity quantum simulations, produce accurate and efficient interatomic potentials. Active learning (AL) is a powerful tool to iteratively generate diverse data sets. In this approach, the ML model provides an uncertainty estimate along with its prediction for each new atomic configuration. If the uncertainty estimate passes a certain threshold, then the configuration is included in the data set. Here we develop a strategy to more rapidly discover configurations that meaningfully augment the training data set. The approach, uncertainty-driven dynamics for active learning (UDD-AL), modifies the potential energy surface used in molecular dynamics simulations to favor regions of configuration space for which there is large model uncertainty. The performance of UDD-AL is demonstrated for two AL tasks: sampling the conformational space of glycine and sampling the promotion of proton transfer in acetylacetone. The method is shown to efficiently explore the chemically relevant configuration space, which may be inaccessible using regular dynamical sampling at target temperature conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Water transport through monolayer fullerene membrane

Water transport through nanoporous materials is important in water treatment, desalination, and nanofiltration. Two-dimensional (2D) membranes such as porous graphene have been explored for high-permeance water transport. However, water transport through a new class of 2D membranes based on two-dimensional covalently linked fullerene monolayers has not been fully explored. Here we use classical molecular dynamics simulations to investigate both vapor and liquid water transport through a monolayer fullerene membrane. We find that a quasi-tetragonal phase fullerene membrane possesses the right pore size and geometry that allows fast water vapor transport (∼ 50 g m −2 day −1 Pa −1 ) and water liquid transport (∼ 2.0 g m −2 day −1 Pa −1 ). Furthermore, simulation of sea water transport through the fullerene membrane shows 100 % salt rejection. The much faster vapor transport rate is attributed to the funnel-shaped pore and the optimal size that allows free rotation of water molecules permeating through, while the slower liquid transport is due to the need to desolvate a water molecule to break its hydrogen-bond network across the hydrophobic pore. This work shows the great potential of using monolayer fullerene membranes as 2D membranes for fast and selective water transport.

Desalination↗
DOE OSTI2021

Ultrahigh thermal conductivity in three-dimensional covalent organic frameworks

Covalent organic frameworks (COFs) are a new, growing class of crystalline polymers with exceptional characteristics such as low density, mechanical strength, and flexibility. Taking three-dimensional (3D) COF-300 derivatives as an example, we apply molecular dynamics simulations to show that ultrahigh thermal conductivity (>15 W/(mK)) is achieved at a pore size of d = 0.63 nm. This is the first observation of such an outstanding intrinsic thermal conductivity of a 3D polymer in all directions. Further temperature-dependent thermal conductivity calculations and X-ray diffraction (XRD) patterns highlight the importance of maintaining microscopic order for a large thermal conductivity. We also demonstrate that the tunability and anisotropy in COFs provide sufficient room to engineer and obtain desired thermal conductivity. These findings provide a fundamental understanding of structure-thermal conductivity relation in COFs and shine a light on high-thermal-conductivity COFs for thermal management applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

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

Bicontinuous nanoporous design induced homogenization of strain localization in metallic glasses

Bicontinuous nanoporous metallic glasses (MG) synergize the outstanding properties of MGs and open-cell nanoporous materials. The low-density and high-specific-surface-area of bicontinuous nanoporous structures have the potential to enhance the applicability of MGs in catalysis, sensors, and lightweight structural designs. In this work, we report molecular dynamics simulations of tensile loading deformation and failure of bicontinuous nanoporous Cu 64 Zr 36 MG with 55% porosity and 4.4 nm ligament size. Results indicate an anomalous mechanical behavior featuring delocalized plastic deformation preceding ductile failure. The deformation follows two mechanisms: i) Necking of ligaments aligned with the loading direction and ii) progressive alignment of randomly oriented ligaments. Failure occurs at 0.16 strain, following massive rupture of ligaments. This work indicates that a bicontinuous nanoporous design is able to effectively delocalize strain localization in a MG due to a combination of size effect on the ductility of MGs resulting in nano ligaments necking and progressive asynchronous alignment of ligaments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗