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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Electromechanics of graphene spirals

Among the most fascinating nanostructure morphologies are spirals, hybrids of somewhat obscure topology and dimensionality with technologically attractive properties. Here, we investigate mechanical and electromechanical properties of graphene spirals upon elongation by using density-functional tight-binding, continuum elasticity theory, and classical force field molecular dynamics. It turns out that electronic properties are governed by interlayer interactions as opposed to strain effects. The structural behavior is governed by van der Waals interaction: in its absence spirals unfold with equidistant layer spacings, ripple formation at spiral perimeter, and steadily increasing axial force; in its presence, on the contrary, spirals unfold via smooth local peeling, complex geometries, and nearly constant axial force. These electromechanical trends ought to provide useful guidelines not only for additional theoretical investigations but also for forthcoming experiments on graphene spirals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗
DOE OSTI2013

Atomistic simulation of electron and phonon transport in nano-devices

In this paper, the electrical, thermal, and coupled electro-thermal properties of silicon nanowire field-effect transistors (NWFETs) are investigated through advanced numerical simulations. An atomistic quantum transport solver based on the Non-equilibrium Green's Function (NEGF) formalism is utilized as modeling tool. The electron characteristics are treated within an empirical tight-binding basis, while a modified valence-force-field method is used to describe phonons. Device simulations going beyond the ballistic limit of transport are performed: electron-phonon interactions with an equilibrium and out-of-equilibrium phonon population are considered as well as an-harmonic phonon-phonon scattering. It is found that ultrascaled Si NWFETs with a gate length of 15 nm operate at about 50% of their ballistic limit. This poor performance can be explained by large self-heating effects that drastically increase the lattice temperature and thus reduce the drain current.

Luisier, Mathieu↗
DOE OSTI2025

3D Continuous Forcing Dataset from 3D Constrained Variational Analysis at SGP

The continuous 3D large-scale forcing (VARANAL3D) data set derived from 3D constrained variational analysis (3DCVA) extends the conventional constrained variational analysis method by incorporating multiple sub-columns within the analysis domain. This advancement introduces spatial variability into the large-scale forcing fields, thereby enriching the data set’s applicability. The VARANAL3D data set spans from 2004 to 2018 and covers a region of 5˚×4.5˚ domain around the ARM SGP site. The analysis domain is divided into 10×9 sub-columns with 0.5˚ resolution. The 3D large-scale forcing data provides necessary variables to drive and evaluate single-column models (SCM), cloud-resolving models (CRM) ,and large-eddy simulations (LES), as well as information for testing model sensitivity to spatial variability of the large-scale forcing data, facilitating more rigorous testing and refinement of physical processes in SCM/CRM/LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Direct measurement of optical force induced by near-field plasmonic cavity using dynamic mode AFM

Plasmonic nanostructures have attracted much attention in recent years because of their potential applications in optical manipulation through near-field enhancement. Continuing experimental efforts have been made to develop accurate techniques to directly measure the near-field optical force induced by the plasmonic nanostructures in the visible frequency range. In this work, we report a new application of dynamic mode atomic force microscopy (DM-AFM) in the measurement of the enhanced optical force acting on a nano-structured plasmonic resonant cavity. The plasmonic cavity is made of an upper gold-coated glass sphere and a lower quartz substrate patterned with an array of subwavelength gold disks. In the near-field when the sphere is positioned close to the disk array, plasmonic resonance is excited in the cavity and the induced force by a 1550 nm infrared laser is found to be increased by an order of magnitude compared with the photon pressure generated by the same laser light. Lastly, the experiment demonstrates that DM-AFM is a powerful tool for the study of light induced forces and their enhancement in plasmonic nanostructures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Conformational Order in Aggregates of Conjugated Polymers

With the abundant variety and increasing chemical complexity of conjugated poly-friers proliferating the field of organic semiconductors, it has become increasingly important to correlate the polymer molecular structure with its mesoscale conformational and morphological attributes. For instance, it is unknown which combinations of chemical moieties and periodicities predictably produce mesoscale ordering. Interestingly) not all ordered morphologies result in efficient devices. In this work we have parametrized accurate classical force-fields and used these to compute the conformational and aggregation characteristics of single strands of common conjugated polymers. Molecular dynamics trajectories are shown to reproduce experimentally observed polymeric ordering, concluding that efficient organic photovoltaic devices span a range of polymer conformational classes, and suggesting that the solution-phase morphologies have far-reaching effects. Encouragingly, these simulations indicate that despite the wide-range of conformational classes present in successful devices, local molecular ordering, and not long-range crystallinity, appears to be the necessary requirement for efficient devices. Finally, we examine what makes a "good" solvent for conjugated polymers, concluding that dispersive pi-electron solvent-polymer interactions, and not the electrostatic potential of the backbone interacting with the solvent, are what primarily determine a polymer's solubility in a particular solvent, and consequently its morphological characteristics.

Jackson, Nicholas E.↗
DOE OSTI2019

Molecular polarizability in open ensemble simulations of aqueous nanoconfinements under electric field

Molecular polarization at aqueous interfaces involves fast degrees of freedom that are often averaged-out in atomistic-modeling approaches. The resulting effective interactions depend on a specific environment, making explicit account of molecular polarizability particularly important in solutions with pronounced anisotropic perturbations, including solid/liquid interfaces and external fields. Our work concerns polarizability effects in nanoscale confinements under electric field, open to an unperturbed bulk environment. We model aqueous molecules and ions in hydrophobic pores using the Gaussian-charge-on-spring BK3-AH representation. This involves nontrivial methodology developments in expanded ensemble Monte Carlo simulations for open systems with long-ranged multibody interactions and necessitates further improvements for efficient modeling of polarizable ions. Structural differences between fixed-charge and polarizable models were captured in molecular dynamics simulations for a set of closed systems. Our open ensemble results with the BK3 model in neat-aqueous systems capture the ∼10% reduction of molecular dipoles within the surface layer near the hydrophobic pore walls in analogy to reported quantum mechanical calculations at water/vapor interfaces. The polarizability affects the interfacial dielectric behavior and weakens the electric-field dependence of water absorption at pragmatically relevant porosities. We observe moderate changes in thermodynamic properties and atom and charged-site spatial distributions; the Gaussian distribution of mobile charges on water and ions in the polarizable model shifts the density amplitudes and blurs the charge-layering effects associated with increased ion absorption. The use of polarizable force field indicates an enhanced response of interfacial ion distributions to applied electric field, a feature potentially important for in silico modeling of electric double layer capacitors.

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

Toward transferable interatomic van der Waals interactions without electrons: The role of multipole electrostatics and many-body dispersion

We estimate polarizabilities of atoms in molecules without electron density, using a Voronoi tesselation approach instead of conventional density partitioning schemes. The resulting atomic dispersion coefficients are calculated, as well as many-body dispersion effects on intermolecular potential energies. We also estimate contributions from multipole electrostatics and compare them to dispersion. We assess the performance of the resulting intermolecular interaction model from dispersion and electrostatics for more than 1300 neutral and charged, small organic molecular dimers. Applications to water clusters, the benzene crystal, the anti-cancer drug ellipticine—intercalated between two Watson-Crick DNA base pairs, as well as six macro-molecular host-guest complexes highlight the potential of this method and help to identify points of future improvement. The mean absolute error made by the combination of static electrostatics with many-body dispersion reduces at larger distances, while it plateaus for two-body dispersion, in conflict with the common assumption that the simple 1/R{sup 6} correction will yield proper dissociative tails. Overall, the method achieves an accuracy well within conventional molecular force fields while exhibiting a simple parametrization protocol.

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

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Sensing polarization effects through the analysis of the effective C{sub 6} dispersion coefficients in NaCl solutions

Density functional theory based ab initio molecular dynamics is used to obtain microscopic details of the interactions in sodium chloride solutions. By following the changes in the atomic C{sub 6} coefficients under the Tkatchenko-Scheffler’s scheme, we were able to identify two different coordination situations for the Cl{sup −} ion with significant different capabilities to perform dispersion interactions. This capability is enhanced when the ion-ion distance corresponds to the contact ion-pair situation. Also, the oxygen and hydrogen atoms of the water molecules change their aptitudes to interact through van der Waals like terms when they are close to the cation region of the ion-pair. These results have interesting implications on the design of force fields to model electrolyte solutions.

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

Kohn–Sham accuracy from orbital-free density functional theory via Δ-machine learning

Here, we present a Δ-machine learning model for obtaining Kohn–Sham accuracy from orbital-free density functional theory (DFT) calculations. In particular, we employ a machine-learned force field (MLFF) scheme based on the kernel method to capture the difference between Kohn–Sham and orbital-free DFT energies/forces. We implement this model in the context of on-the-fly molecular dynamics simulations and study its accuracy, performance, and sensitivity to parameters for representative systems. We find that the formalism not only improves the accuracy of Thomas–Fermi–von Weizsäcker orbital-free energies and forces by more than two orders of magnitude but is also more accurate than MLFFs based solely on Kohn–Sham DFT while being more efficient and less sensitive to model parameters. We apply the framework to study the structure of molten Al 0.88 Si 0.12 , the results suggesting no aggregation of Si atoms, in agreement with a previous Kohn–Sham study performed at an order of magnitude smaller length and time scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Flow of long chain hydrocarbons through carbon nanotubes (CNTs)

The pressure-driven flow of long-chain hydrocarbons in nanosized pores is important in energy, environmental, biological, and pharmaceutical applications. This paper examines the flow of hexane, heptane, and decane in carbon nanotubes (CNTs) of pore diameters 1–8 nm using molecular dynamic simulations. Enhancement of water flow in CNTs in comparison to rates predicted by continuum models has been well established in the literature. Our work was intended to observe if molecular dynamic simulations of hydrocarbon flow in CNTs produced similar enhancements. We used the OPLS-AA force field to simulate the hydrocarbons and the CNTs. Our simulations predicted the bulk densities of the hydrocarbons to be within 3% of the literature values. Molecular sizes and shapes of the hydrocarbon molecules compared to the pore size create interesting density patterns for smaller sized CNTs. We observed moderate flow enhancements for all the hydrocarbons (1–100) flowing through small-sized CNTs. For very small CNTs the larger hydrocarbons were forced to flow in a cork-screw fashion. As a result of this flow orientation, the larger molecules flowed as effectively (similar enhancements) as the smaller hydrocarbons.

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

Molecular polarizability of water from local dielectric response theory

Here, we propose a fully ab initio theory to compute the electron density response under the perturbation in the local field. This method is based on our recently developed local dielectric response theory [Phys. Rev. B 92, 241107(R), 2015], which provides a rigorous theoretical framework to treat local electronic excitations in both nite and extended systems beyond the commonly employed dipole approximation. We have applied this method to study the electronic part of the molecular polarizability of water in ice Ih and liquid water. Our results reveal that the crystal field of the hydrogen-bond network has strong anisotropic effects, which significantly enhance the out-of-plane component and suppress the in-plane component perpendicular to the bisector direction. The contribution from the charge transfer is equally important, which increases the isotropic molecular polarizability by 5-6%. Our study provides new insights into the dielectric properties of water, which form the basis to understand electronic excitations in water and to develop accurate polarizable force fields of water.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Phonon-limited low-field mobility in silicon: Quantum transport vs. linearized Boltzmann Transport Equation

We propose to check and validate the approximations made in dissipative quantum transport (QT) simulations solved in the Non-equilibrium Green's Function formalism by comparing them with the exact solution of the linearized Boltzmann Transport Equation (LB) in the stationary regime. For that purpose, we calculate the phonon-limited electron and hole mobility in bulk Si and ultra-scaled Si nanowires for different crystal orientations 〈100〉, 〈110〉, and 〈111〉. In both QT and LB simulations, we use the same sp{sup 3}d{sup 5}s{sup *} tight-binding model to describe the electron/hole properties and the same valence-force-field approach to account for the phonon properties. It is found that the QT simplifications work well for electrons, but are less accurate for holes, where a renormalization of the phonon scattering strength is proved useful to improve the results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2016

Complexation Enhancement Drives Water-to-Oil Ion Transport: A Simulation Study

We address the structures and energetics of ion solvation in aqueous and organic solutions to understand liquid-liquid ion transport. Atomistic molecular dynamics (MD) simulations with polarizable force field are performed to study the coordination transformations driving lanthanide (Ln(III)) and nitrate ion transport between aqueous and an alkylamide-oil solution. An enhancement of the coordination behavior in the organic phase is achieved in contrast with the aqueous solution. In particular, the coordination number of Ce3+ increases from 8.9 in the aqueous to 9.9 in the organic solutions (from 8 in the aqueous to 8.8 in the organic systems for Yb3+). Moreover, the local coordination environ ment changes dramatically. Potential of mean force calculations show that the Ln(III)-ligand coordination interaction strengths follow the order of Ln(III-)nitrate> Ln(III)-water>Ln(III)-DMDBTDMA. They increase 2-fold in the lipophilic environment in comparison to the aqueous phase, and we attribute this to the shedding of the outer solvation shell. Our findings highlight the importance of outer sphere interactions on the competitive solvation energetics that cause ions to migrate between immiscible phases; an essential ingredient for advancing important applications such as rare earth metal separations. Some open questions in simulating the coordination behavior of heavy metals are also addressed.

Qiao, Baofu↗