Search NASA⌕ Search

SEARCH · Search NASA

Results for “force field development”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations

Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.

42 ENGINEERING↗
DOE OSTI2025

Tunable Topological Phonon for Next-generation Quantum Transduction

This project aims to understand the critical factors that determine transduction performance of topological phonons across an oxide perovskite/tungsten diselenide heterojunction. We will investigate mechanisms of their propagation and interfacial coupling using modeling and machine learning approaches. Methods include density functional theory, molecular dynamics, numerical transport simulations, and active learning for building up a training dataset for force field development.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1996

Computational Materials Research

Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.

Hinkley, Jeffrey A.↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational Exploration of High Entropy Alloys as Promising Materials for Future Beam Windows

With the ever-increasing demand for high beam power, the currently used beam-intercepting devices (BIDs) such as targets, and beam windows may not be able to handle the high power required for future accelerator complexes or the lifetime may be reduced drastically. As beam power increases, the damage incurred by BIDs, including thermal shock, fatigue, and irradiation damage, also rises. Therefore, it is imperative to design materials that can withstand high beam power for longer lifetimes. High entropy alloys (HEAs) have emerged as potential alternative materials for designing next-generation BIDs. In this study, we primarily focus on materials for developing beam windows for next-generation accelerator complexes. We propose an integrated approach that combines various computational techniques to study and design new materials. Specifically, we use CALPHAD, density functional theory (DFT), and molecular dynamics (MD) to comprehensively investigate the defect properties of suitable HEAs, offering potential alternatives for future beam windows. We begin by scanning the extensive phase space provided by Cr-Mn-V-Ti-Al-Co HEAs, selecting 8 compositions after evaluating approximately 120,000 unique compositions using CALPHAD. We, then employ DFT-informed machine learning techniques to develop force-field parameters. Finally, MD simulations using these developed force-field parameters will be used to study the effects of radiation damage on the defect and mechanical properties of the selected alloys. This research explains the use of the CALPHAD approach and shows how critical modeling (DFT and MD) is in developing novel material such as HEAs. It also highlights the promising role of machine learning in this field. The results from this study will greatly improve the novel materials development to be used in next-generation accelerator components, leading to higher beam power and longer operational times of BIDs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Development of Field Excavator with Embedded Force Measurement

A semi-intelligent excavation mechanism was developed for use with the NASA-built Centaur 2 rover prototype. The excavator features a continuously rotatable large bucket supported between two parallel arms, both of which share a single pivot axis near the excavator base attached to the rover. The excavator is designed to simulate the collection of regolith, such as on the Moon, and to dump the collected soil into a hopper up to one meter tall for processing to extract oxygen. Because the vehicle can be autonomous and the terrain is generally unknown, there is risk of damaging equipment or using excessive power when attempting to extract soil from dense or rocky terrain. To minimize these risks, it is critical for the rover to sense the digging forces and adjust accordingly. It is also important to understand the digging capabilities and limitations of the excavator. This paper discusses the implementation of multiple strain gages as an embedded force measurement system in the excavator's arms. These strain gages can accurately measure and resolve multi-axial forces on the excavator. In order to validate these sensors and characterize the load capabilities, a series of controlled excavation tests were performed at Glenn Research Center with the excavator at various depths and cut angles while supported by a six axis load cell. The results of these tests are both compared to a force estimation model and used for calibration of the embedded strain gages. In addition, excavation forces generated using two different types of bucket edge (straight vs. with teeth) were compared.

Soils↗
DOE OSTI2020

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

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

Separation of Linear and Cyclic Siloxanes in Pure Silica Zeolites

We present a computational assessment of pure-silica zeolites for separating linear and cyclic siloxanes. We developed a force field (FF) for pure silica zeolites, when combined with our previously developed FF for siloxanes using standard Lorentz–Berthelot combining rules, shows good agreement with dispersion-corrected density functional theory calculations. We used molecular dynamics simulations to investigate diffusion of siloxanes in pure silica zeolites and identified a pure silica zeolite with the structure code FAU that enables kinetic separation of linear and cyclic siloxanes. FAU allows the diffusion of linear siloxanes (L2–L6) while excluding cyclic siloxanes. D4 siloxane does not diffuse in any of the investigated zeolites, eliminating the potential of pure silica zeolites to achieve equilibrium-based separations of linear and cyclic siloxanes.

adsorption↗
DOE OSTI2023

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗
DOE OSTI2024

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Finally, library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Development of a Systematic and Extensible Force Field for Peptoids (STEPs)

Peptoids (N-substituted glycines) are a class of biomimetic polymers that have attracted significant attention due to their accessible synthesis and enzymatic and thermal stability relative to their naturally occurring counterparts (polypeptides). While these polymers provide the promise of more robust functional materials via hierarchical approaches, they present a new challenge for computational structure prediction for material design. The reliability of calculations hinges on the accuracy of interactions represented in the force field used to model peptoids. For proteins, structure prediction based on sequence and de novo design has made dramatic progress in recent years; however, these models are not readily transferable for peptoids. Current efforts to develop and implement peptoid-specific force fields are spread out, leading to replicated efforts and a fragmented collection of parameterized sidechains. Here, we developed a peptoid-specific force field containing 70 different side chains, using GAFF2 as starting point. The new model is validated based on the generation of Ramachandran-like plots from DFT optimization compared against force field reproduced potential energy and free energy surfaces as well as the reproduction of equilibrium cis/trans values for some residues experimentally known to form helical structures. In conclusion, equilibrium cis/trans distributions (Kct) are estimated for all parameterized residues to identify which residues have an intrinsic propensity for cis or trans states in the monomeric state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from a source surface, given the proper boundary conditions. The results of this work; describing the mathematical formalism that was developed, the numerical techniques employed, and the stability criteria developed for these numerical schemes are presented. An analytical solution is used for a test case; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent ( 5%).

Wu, S. T.↗
DOE OSTI2025

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗