Search NASA⌕ Search

SEARCH · Search NASA

Results for “force fields”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

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

Surface structure of linear nanopores in amorphous silica: Comparison of properties for different pore generation algorithms

In this work, we compare the surface structure of linear nanopores in amorphous silica (a-SiO 2 ) for different versions of “pore drilling” algorithms (where the pores are generated by the removal of atoms from the preformed bulk a-SiO 2 ) and for “cylindrical resist” algorithms (where a-SiO 2 is formed around a cylindrical exclusion region). After adding H to non-bridging O, the former often results in a moderate to high density of surface silanol groups, whereas the latter produces a low density. The silanol surface density for pore drilling can be lowered by a final dehydroxylation step, and that for the cylindrical resist approach can be increased by a final hydroxylation step. In this respect, the two classes of algorithms are complementary. We focus on the characterization of the chemical structure of the pore surface, decomposing the total silanol density into components corresponding to isolated and vicinal mono silanols and geminal silanols. The final dehyroxylation and hydroxylation steps can also be tuned to better align some of these populations with the target experimental values.

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

Effect of water models on structure and dynamics of lignin in solution

Lignin, a major biomass component, can be an excellent source for different monomers in the polymer industry. However, the complex and heterogeneous structure of lignin poses a significant challenge for designing energy-efficient processes for depolymerization. As many proposed depolymerization processes are solvothermal, it is essential to understand the structure and dynamics of lignin in solution, in particular aqueous solution. Here, we utilize molecular dynamics simulations to understand the effect of water models on the structure and dynamics of different model lignin oligomers (softwood and hardwood) as a function of temperature. We have examined three different water models: TIP3P, TIP4P/Ew, and flexible SPC/Fw. We find that the diffusion constant of lignin oligomers in an aqueous solution differs significantly depending on the water model used. The diffusion constant of lignin in the TIP3P water model is almost twice as large as that in SPC/Fw and TIP4P/Ew. The softwood and hardwood oligomers adopt an extended structure in TIP3P water compared to SPC/Fw and TIP4P/Ew. Given the different levels of sensitivity of transport and structural properties of aqueous lignin on water models, it is important to take these into account when discussing results from a specific water model.

09 BIOMASS FUELS↗
DOE OSTI2021

First principles reactive simulation for equation of state prediction

The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.

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

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

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

Can a deep-learning model make fast predictions of vacancy formation in diverse materials?

The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.

2D materials↗
DOE OSTI2023

Detailed studies of the processes in low energy H irradiation of Li and Li-compound surfaces

We have used a combination of pico-to-nano temporal/spatial scale computational physics and chemistry modeling of plasma–material interfaces in the tokamak fusion plasma edges to unravel the evolving characteristics, not readily accessible by empirical means, of lithium-, oxygen-, and hydrogen-containing materials of plasma-facing components under irradiation by hydrogen and its isotopes. In the present calculation, amorphous lithium compound surfaces containing oxygen, Li 2 O, and LiOH were irradiated by 1-100 eV particles at incident angles on the surface ranging from perpendicular to almost grazing angles. Consequential surface processes, reflection, retention, and sputtering were studied at “the same footing” and compared to earlier results from amorphous Li and LiH surfaces. In conclusion, the critical role of charging dynamics of lithium, oxygen, and hydrogen atoms in the surface chemistry during hydrogen-fuel irradiation was found to drive the kinetics and dynamics of these surfaces in unexpected ways that ultimately could have profound effects on fusion plasma confinement behavior and surface erosion.

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

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

The influence of counterion structure identity on conductivity, dynamical correlations, and ion transport mechanisms in polymerized ionic liquids

We used equilibrium and non-equilibrium atomistic simulations to probe the influence of anion chemistry on the true conductivity, dynamical correlations, and ion transport mechanisms in polymeric ionic liquids. An inverse correlation was found between anion self-diffusivities, ionic mobilities, and the anion size for spherical anions. While some larger asymmetric anions had higher diffusivities than smaller spherical anions, their diffusivities and mobilities did not exhibit a direct correlation to the anion volumes. Here, the conductivity and anion dynamical correlations also followed the same trends as displayed by the diffusivity and mobility of anions. All the systems we examined displayed positively correlated motion among anions, suggesting a contribution that enhances the conductivity beyond the ideal Nernst–Einstein value. Analysis of ion transport mechanisms demonstrated very similar hopping characteristics among the spherical anions despite differences in their sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2019

Water follows polar and nonpolar protein surface domains

The conformation of water around proteins is of paramount importance, as it determines protein interactions. Although the average water properties around the surface of proteins have been provided experimentally and computationally, protein surfaces are highly heterogeneous. Therefore, it is crucial to determine the correlations of water to the local distributions of polar and nonpolar protein surface domains to understand functions such as aggregation, mutations, and delivery. By using atomistic simulations, we investigate the orientation and dynamics of water molecules next to 4 types of protein surface domains: negatively charged, positively charged, and charge-neutral polar and nonpolar amino acids. The negatively charged amino acids orient around 98% of the neighboring water dipoles toward the protein surface, and such correlation persists up to around 16 Å from the protein surface. The positively charged amino acids orient around 94% of the nearest water dipoles against the protein surface, and the correlation persists up to around 12 Å. The charge-neutral polar and nonpolar amino acids are also orienting the water neighbors in a quantitatively weaker manner. A similar trend was observed in the residence time of the nearest water neighbors. These findings hold true for 3 technically important enzymes (PETase, cytochrome P450, and organophosphorus hydrolase). Our results demonstrate that the water−amino acid degree of correlation follows the same trend as the amino acid contribution in proteins solubility, namely, the negatively charged amino acids are the most beneficial for protein solubility, then the positively charged amino acids, and finally the charge-neutral amino acids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Time-resolved particle image velocimetry

Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. Lastly, a survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparison of three DWM-based wake models at above-rated wind speeds

In this study we investigate three mid-fidelity wind turbine wake models based on the dynamic wake meandering (DWM) model principle, and compare their performance with a reference dataset, produced with large-eddy simulations using the actuator line model. The models are compared with respect to flow field, power, and loads on a row of four 5MW reference turbines experiencing above-rated wind conditions. In general, the DWM models show fairly good agreement with large-eddy simulation for the time-averaged flow fields, blade forces and power, with increasing differences along the turbine row. Also when comparing fatigue loads of blade root moments, the differences between the models increase further into the row, with deviations up to 25 % of the reference case. However, while the development in blade root moment fatigue along the turbine row is predominantly driven by the energy content at the frequency corresponding to the turbine's rotational period (1 P ) for the DWM models, the large-eddy simulation results suggest that the key drivers for the blade root and tower loads are the increase in meandering and energy at higher frequencies (> 1 P ) deeper into the turbine row. For the tower loads, the DWM models highly underestimate the fatigue for the waked turbines. From these results, we suggest priorities for future model developments so that robust model implementations can be used in wind farm design and operation.

17 WIND ENERGY↗