Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2000

Modeling Molecules

The molecule modeling method known as Multibody Order (N) Dynamics, or MBO(N)D, was developed by Moldyn, Inc. at Goddard Space Flight Center through funding provided by the SBIR program. The software can model the dynamics of molecules through technology which stimulates low-frequency molecular motions and properties, such as movements among a molecule's constituent parts. With MBO(N)D, a molecule is substructured into a set of interconnected rigid and flexible bodies. These bodies replace the computation burden of mapping individual atoms. Moldyn's technology cuts computation time while increasing accuracy. The MBO(N)D technology is available as Insight II 97.0 from Molecular Simulations, Inc. Currently the technology is used to account for forces on spacecraft parts and to perform molecular analyses for pharmaceutical purposes. It permits the solution of molecular dynamics problems on a moderate workstation, as opposed to on a supercomputer.

Source record↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗
DOE OSTI2024

Early Career Researcher Support for Foundations of Molecular Modeling and Simulation Conference

This project enabled 14 early career researchers, including graduate students and postdocs, to attend the 8th triennial conference on Foundations of Molecular Modeling and Simulation (FOMMS 2022), which took place in Delavan, WI (Lake Geneva area), from July 17 to 21, 2022. Molecular-level modeling and simulation, which serve as critical tools for advancing science and engineering in a myriad of applications, were at the forefront of discussions. These ranged from the design of new materials for various energy applications to insights into the interaction between drug molecules and protein receptors within the body.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Near-physiological-temperature serial crystallography reveals conformations of SARS-CoV-2 main protease active site for improved drug repurposing

The COVID-19 pandemic has resulted in 198 million reported infections and more than 4 million deaths as of July 2021 (covid19.who.int). Research to identify effective therapies for COVID-19 includes: (1) designing a vaccine as future protection; (2) de novo drug discovery; and (3) identifying existing drugs to repurpose them as effective and immediate treatments. To assist in drug repurposing and design, we determine two apo structures of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease at ambient temperature by serial femtosecond X-ray crystallography. We employ detailed molecular simulations of selected known main protease inhibitors with the structures and compare binding modes and energies. The combined structural and molecular modeling studies not only reveal the dynamics of small molecules targeting the main protease but also provide invaluable opportunities for drug repurposing and structure-based drug design strategies against SARS-CoV-2.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Simulation of supersonic molecular beam injection fueling into H-mode plasmas on EAST using BOUT++

Transport dynamic 2D simulations for supersonic molecular beam injection (SMBI) fueling into H-mode deuterium plasmas on the Experimental Advanced Superconducting Tokamak (EAST) is first simulated using a seven-field two-fluid model in the BOUT++ framework. The SMB is assumed to be injected into plasma from the midplane at the low field side with a fixed width and constant molecular flux. The different densities and injection velocities of SMBI are investigated within the upper single-null geometry of EAST. The simulations indicate that the SMBI has a self-shielding effect on molecules' inward transport into the plasma, and the deposition of SMBI leads to a large increase in plasma density and decrease in plasma temperature. There is a velocity threshold for SMB penetrating the pedestal and depositing at the top of density pedestal. However, the deposition point would be back toward the plasma boundary after SMB arriving at the deepest penetration position. Comparing the different molecular injection velocities and densities, the outcomes demonstrate that the depth of deposition is closely related to the injection velocity rather than the molecular density. The simulated results show good agreement with the EAST experiments by comparing the electron density profiles obtained by simulation and experiment, respectively. These results will be helpful for guiding future experiments, as well as the design of the SMBI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Incorporating Electronic Information into Machine Learning Potential Energy Surfaces via Approaching the Ground-State Electronic Energy as a Function of Atom-Based Electronic Populations

Machine learning (ML) approximations to density functional theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present, they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new machine learning approach called "BpopNN" for obtaining efficient approximations to DFT PESs. Conceptually, the methodology is based on approaching the true DFT energy as a function of electron populations on atoms; in practice, this is realized with available density functionals and constrained DFT (CDFT). Additionally, the new approach creates approximations to this function with neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on LinHn clusters.

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

Reticular Materials and AI-Driven Computer Simulations for Seawater Mining of Valuable Metals (Final Technical Report)

This Final Technical Report describes our exploratory efforts that combine reticular materials synthesis (hydrolytically robust metal–organic frameworks, MOFs) with AI‑enabled molecular simulations to develop mechanistic, quantitative design rules for recovering lithium and other alkali-metal ions from highly dilute, competitive aqueous resources (e.g., seawater). The central outcome is a joint experimental–computational study of ion uptake in MOF‑808 (Chemical Science, 2025) that quantifies both thermodynamics and kinetics of Li + , Na + , and K + uptake and identifies how pore size, pore hydration state, dehydration penalties, and pore-window transport barriers govern selectivity. Guided by these insights, we synthesized and tested functionalized MOF‑808 and multivariate MOFs incorporating ion-recognition motifs (including carboxylates and crown-ether linkers) and evaluated uptake in synthetic seawater, highlighting framework topology and pore chemistry as levers for improved Li + /Na + discrimination. We also developed transferable simulation models, enhanced-sampling protocols, and automated workflows that enable systematic screening of porous sorbents.

42 ENGINEERING↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
DOE OSTI2021

Editorial: Advances in computational molecular biophysics

This special issue of BBA is dedicated to Dr. Jeremy C. Smith on the occasion of his 60th birthday and his more than three decades as a practitioner of computational molecular biophysics. The special issue showcases cutting-edge work in computational molecular biophysics. The senior authors of the manuscripts in this issue are Dr. Smith’s previous or present collaborators. Their articles reflect Dr. Smith’s distinguished career integrating experimental and molecular simulation approaches into a unified description of the function, structure, and dynamics of biomolecules. Each of the articles presented here illustrates one aspect of the broad nature of computational molecular biophysics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Simulation studies of polypeptoids using replica exchange with dynamical scaling and dihedral biasing

Polypeptoids differ from polypeptides in that the amide bond can more frequently adopt both cis and trans conformations. The transition between the two conformations requires overcoming a large energy barrier, making it difficult for conventional molecular simulations to adequately visit the cis and trans structures. A replica-exchange method is presented that allows for easy rotations of the amide bond and also an efficient linking to a high temperature replica. The method allows for just three replicas (one at the temperature and Hamiltonian of interest, a second high temperature replica with a biased dihedral potential, and a third connecting them) to overcome the amide bond sampling problem and also enhance sampling for other coordinates. Here our results indicate that for short peptoid oligomers, the conformations can range from all cis to all trans with an average cis/trans ratio that depends on side chain and potential model.

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

Scattering evidence of positional charge correlations in polyelectrolyte complexes

Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral polymers, where the total scattering function exhibits an Ornstein–Zernike form. In spite of consensus among different theoretical predictions, the existence of positional correlations between polyanion and polycation charges has not been confirmed experimentally. Here, we present small-angle neutron scattering profiles where the polycation scattering length density is matched to that of the solvent to extract positional correlations among anionic monomers. The polyanion scattering functions exhibit a peak at the inverse polymer screening radius of Coulomb interactions, q* ≈ 0.2 Å –1 . This peak, attributed to Coulomb repulsions between the fragments of polyanions and their attractions to polycations, is even more pronounced in the calculated charge scattering function that quantifies positional correlations of all polymer charges within the PEC. Screening of electrostatic interactions by adding salt leads to the gradual disappearance of this correlation peak, and the scattering functions regain an Ornstein–Zernike form. Experimental scattering results are consistent with those calculated from the random phase approximation, a scaling analysis, and molecular simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent Cα model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Curated Collection of More than 20,000 Experimentally Reported One-Dimensional Metal–Organic Frameworks

Herein, a collection of more than 20,000 experimentally derived crystal structures for metal–organic frameworks (MOFs) that do not have two- or three-dimensional covalently bonded networks has been developed from the materials available at the Cambridge Crystallographic Data Centre. Of these 20,000 1D MOFs, more than 12,000 structures have been verified to be solvent-free and in exact agreement with the stoichiometry of the synthesized materials. More than 10% of the complete data set comprise materials including two or more distinct metals. The band gaps of more than 12,000 1D MOFs have been computed at the density functional theory-generalized gradient approximation level, finding more than 2000 materials that have a zero band gap. Molecular simulations of CH4 adsorption in a small number of 1D MOFs indicated that adsorbate-induced deformation plays a significant role in determining adsorption isotherms in these materials. As a result, methods that have been used previously for high-throughput predictions of molecular adsorption in 3D MOFs are not suitable for 1D MOFs.

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

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations—molecular dynamics in combination with density functional theory calculations—to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4+) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Nonvariational ADAPT algorithm for quantum simulations

We explore a nonvariational quantum state preparation approach combined with the ADAPT operator selection strategy in the application of preparing the ground state of a desired target Hamiltonian. In this algorithm, energy gradient measurements determine both the operators and the gate parameters in the quantum circuit construction. We compare this nonvariational algorithm with ADAPT-VQE and with feedback-based quantum algorithms in terms of the rate of energy reduction, the circuit depth, and the measurement cost in molecular simulation. We find that, despite using deeper circuits, this new algorithm reaches chemical accuracy at a similar measurement cost to ADAPT-VQE. Since it does not rely on a classical optimization subroutine, it may provide robustness against circuit parameter errors due to imperfect control or gate synthesis.

Tang'S, Ho Lun [Virginia Polytechnic Inst. and Sta↗