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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computation of the Thermal Expansion Coefficient of Graphene with Gaussian Approximation Potentials

Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this work, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

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

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗
DOE OSTI2024

Elucidating the Structure of the Eu‐EDTA Complex in Solution at Various Protonation States

Abstract Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu 3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu‐EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X‐ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu 3+ aqua ion and Eu‐EDTA complexes at pH 3 and 11. This work details how Eu‐EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu 3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation.

Chemistry↗
DOE OSTI2021

Hoobas: A highly object-oriented builder for molecular dynamics

Polydispersity and random sequences are ubiquitous features of polymers, and molecular dynamics simulations can help elucidate the impact of disorder in polymer systems. However, currently available packages for building polymer topologies do not enable the user to include randomness in a straightforward fashion. Here, we introduce Hoobas, a molecular builder package that easily handles polydispersity using a prototype-builder design pattern. This enables fast and easy building of systems comprised of thousands of distinct objects. It is written in the Python programming language, which ensures compatibility with a wide range of molecular dynamics packages and tools, as well as easy integration into most workflows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗
DOE OSTI2022

Large Scale MD to Predict Epitope Regions in HIV Env [Slides]

Highly dense carbohydrates located on the surface of the HIV Env protein play a key role in immune evasion. Such evolutionary adaptation hampers any attempt to obtain a full mechanistic understanding of the role played by the glycans in protecting the virus against an effective immune response. Moreover, and due to their chemical variability, an accurate molecular understanding of the so called “glycan shield” is still limited by the lack of effective resolution of state-of-the-art experimental technics. Here, we have used extensive computational modelling in order to fill this gap, addressing the presence of a large glycan variability as observed experimentally. Based on an automated pipeline, we were able to assemble, set-up and simulate via Molecular dynamics hundreds of different glycosylated Env variants at nearly atomic resolution, recapitulating the glycosylation distributions observed experimentally. Results from these simulations were subjected to machine learning and very accurate prediction of simulation derived glycan shielding areas of each glycan as a function of static sequence features. Such predictive models of per-glycan shielding, incorporating both glycan dynamics and heterogeneity, were used to develop a novel sequence-based glycan shield mapping strategy. Parallel to these studies, we also developed an accurate machine learning approach to predict glycan heterogeneity data using sequence features and found good prediction accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2023

Simulations of Glass Transition and Mechanical Behavior of Off-Stoichiometric Crosslinked Polymers

Here, this work explores the influence of blend composition, network architecture, and hydrogen bonding on the material properties of crosslinked epoxy networks, focusing on the glass transition temperature (T g ) and Young’s modulus (Y). We used coarse-grained molecular dynamics simulations to simulate varying compositions of stiff and flexible components in epoxy monomer blends with varying excess of curative. We find that, without hydrogen bonding, networks of any composition show a monotonically increasing T g with decreasing excess curative, consistent with theory. In contrast, we find that when hydrogen bonding is introduced, the binary blend networks show significant enhancement in T g for lightly crosslinked systems. This result contributes to an explanation of the anomalous T g behavior observed experimentally in these systems. We further find that Y is generally enhanced by hydrogen bonds, especially below T g , demonstrating that hydrogen bonding has a significant influence on mechanical properties and can allow access to other desirable dynamic behavior, especially self-healing.

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

Modeling bicarbonate formation in an alkaline solution with multi-level quantum mechanics/molecular dynamics simulations

Understanding carbonate speciation and how it may be modulated is essential for the advancement of carbon dioxide (CO 2 ) capture and storage technologies, which often rely on the transformation of CO 2 into carbonate, e.g. via the formation of carbonate minerals. To date, few atomic-level, quantum-mechanics-based simulations have been carried out to characterize how carbonic acid (H 2 CO 3 ) and bicarbonate ($HCO^{-}_{3}$) form in aqueous solution, and how pH affects this process. Recently, Martirez and Carter utilized rare-event sampling density functional theory molecular dynamics simulations in combination with multi-level embedded correlated wavefunction theory, thus accounting for both solvent dynamics and electron correlation accurately, to elucidate the mechanism of H 2 CO 3 formation in neutral solution (J. Am. Chem. Soc., 145, 12561, 2023). Here, we perform a complementary simulation using the same method to map out the energetics of $HCO^{-}_{3}$ formation from dissolved CO 2 in basic solution. We find that, as in H 2 CO 3 formation, including water dynamics is important to obtain an accurate prediction of the energetics for the aforementioned reaction. Furthermore, only with MD did we identify the correct pathway for the reaction, in which water – not hydroxide – acts as the initial nucleophile and only at the transition state does it lose a proton.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Effect of an oxide layer on high velocity impact of tantalum particles characterized using molecular dynamics

Micro-cold spray, also commonly referred to as the aerosol deposition method or vacuum kinetic spraying, is a process for depositing nanostructured films of metals and ceramics by impacting solid nanoparticles at high velocity onto a substrate at room temperature and at low pressures. Tantalum films deposited by micro-cold spray are desirable because tantalum is non-reactive to certain molten metals and can therefore serve as a non-reactive barrier coating. However, unlike most other metals deposited by micro-cold spray, tantalum quickly forms an oxide layer when exposed to air. In the deposition of films using traditional cold spray with much larger particles, native oxides must fracture for particle–substrate bonding to occur. Furthermore, the fraction of the particle that is oxidized in micro-cold spray and the strain rates experienced by the particles upon impact are both much larger than in cold spray, suggesting that the role of the oxide on film deposition may be distinct. The effect of this oxide layer on particle deformation and adherence to a substrate under high strain rates is studied using molecular dynamics simulations. The effect of oxide layer thickness, particle diameter, and impact velocity are examined.

42 ENGINEERING↗
DOE OSTI2020

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of adiabatic quantum computing for molecular ground states

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.

Kremenetski, Vladimir↗
DOE OSTI2023

Molecular pathway and structural mechanism of human oncochannel TRPV6 inhibition by the phytocannabinoid tetrahydrocannabivarin

The calcium-selective oncochannel TRPV6 is an important driver of cell proliferation in human cancers. Despite increasing interest of pharmacological research in developing synthetic inhibitors of TRPV6, natural compounds acting at this channel have been largely neglected. On the other hand, pharmacokinetics of natural small-molecule antagonists optimized by nature throughout evolution endows these compounds with a medicinal potential to serve as potent and safe next-generation anti-cancer drugs. Here we report the structure of human TRPV6 in complex with tetrahydrocannabivarin (THCV), a natural cannabinoid inhibitor extracted from Cannabis sativa . We use cryo-electron microscopy combined with electrophysiology, calcium imaging, mutagenesis, and molecular dynamics simulations to identify THCV binding sites in the portals that connect the membrane environment surrounding the protein to the central cavity of the channel pore and to characterize the allosteric mechanism of TRPV6 inhibition. We also propose the molecular pathway taken by THCV to reach its binding site. Our study provides a foundation for the development of new TRPV6-targeting drugs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗
DOE OSTI2020

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗