Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamic 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Solid-state Precipitation in Molecular Dynamics: kMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗
NASA NTRS2020

Ultrafast Melting at Metal-Semiconductor Interfaces: A Molecular Dynamics Study

Quantifying the thermodynamics of interfacial melting is important for several applications such as laser machining and designing radiation hardened electronic devices. Due to the short time and length scales that are involved in the process, molecular simulations provide an ideal path for understating the interfacial melting phenomena. In this work we investigate tungsten (W) – silicon carbide (SiC) interfaces which are relevant in modern high-power electronic devices. Classical molecular dynamics simulations along with the two-temperature model (TTM) were used to study the radiation induced melting of the interface while a neural network-based machine learning framework was used to track the transient re-crystallization process. The damage endured by the semiconductor at the interface following the radiation event is discussed in detail while the re-crystallization of the metal is described in terms of the classical nucleation theory.

Ashwin Ravichandran↗
NASA NTRSDate not supplied

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics↗
NASA NTRS2014

Modeling Materials: Design for Planetary Entry, Electric Aircraft, and Beyond

NASA missions push the limits of what is possible. The development of high-performance materials must keep pace with the agency's demanding, cutting-edge applications. Researchers at NASA's Ames Research Center are performing multiscale computational modeling to accelerate development times and further the design of next-generation aerospace materials. Multiscale modeling combines several computationally intensive techniques ranging from the atomic level to the macroscale, passing output from one level as input to the next level. These methods are applicable to a wide variety of materials systems. For example: (a) Ultra-high-temperature ceramics for hypersonic aircraft-we utilized the full range of multiscale modeling to characterize thermal protection materials for faster, safer air- and spacecraft, (b) Planetary entry heat shields for space vehicles-we computed thermal and mechanical properties of ablative composites by combining several methods, from atomistic simulations to macroscale computations, (c) Advanced batteries for electric aircraft-we performed large-scale molecular dynamics simulations of advanced electrolytes for ultra-high-energy capacity batteries to enable long-distance electric aircraft service; and (d) Shape-memory alloys for high-efficiency aircraft-we used high-fidelity electronic structure calculations to determine phase diagrams in shape-memory transformations. Advances in high-performance computing have been critical to the development of multiscale materials modeling. We used nearly one million processor hours on NASA's Pleiades supercomputer to characterize electrolytes with a fidelity that would be otherwise impossible. For this and other projects, Pleiades enables us to push the physics and accuracy of our calculations to new levels.

Supercomputing↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
NASA NTRS2000

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).

NASA Center ARC↗
NASA NTRS2012

Development and Application of Interatomic Potentials for Ultra High Temperature Ceramics (UHTC): ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are characterized by high melting point, good strength, and reasonable oxidation resistance. These materials are of interest for use as sharp leading edges for hypersonic vehicles among other applications. Progress in computational modeling of UHTCs has been limited in part due to the absence of suitable interatomic potentials. We present a Tersoff style parameterization of such potentials for ZrB2 and HfB2 appropriate for atomistic simulations. Parameters are fit to data generated from ab initio calculations. The accuracy of the potentials is assessed against further ab initio data. As a first non ]trivial application, molecular dynamics simulations are performed to evaluate the thermal conductivity of single crystals and the thermal resistance of high symmetry grain boundaries.

Lawson, John W.↗
NASA NTRS2009

Prediction of Material Properties of Nanostructured Polymer Composites Using Atomistic Simulations

Atomistic models of epoxy polymers were built in order to assess the effect of structure at the nanometer scale on the resulting bulk properties such as elastic modulus and thermal conductivity. Atomistic models of both bulk polymer and carbon nanotube polymer composites were built. For the bulk models, the effect of moisture content and temperature on the resulting elastic constants was calculated. A relatively consistent decrease in modulus was seen with increasing temperature. The dependence of modulus on moisture content was less consistent. This behavior was seen for two different epoxy systems, one containing a difunctional epoxy molecule and the other a tetrafunctional epoxy molecule. Both epoxy structures were crosslinked with diamine curing agents. Multifunctional properties were calculated with the nanocomposite models. Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between the carbon nanotube and the surrounding epoxy matrix. These estimated values were used in a multiscale model in order to predict the thermal conductivity of a nanocomposite as a function of the nanometer scaled molecular structure.

Hinkley, J.A.↗
NASA NTRS1999

Computer Simulation Studies of Ion Channel Gating: Characteristics of the M2 Channel of Influenza-A Virus in a Phospholipid Bilayer

The 25 amino acids long, transmembrane fragment of the Influenza virus M2 protein forms a homotetrameric channel that transports protons across lipid bilayers. It has been postulated that high efficiency and selectivity of this process is due to gating by four histidine residues that occlude the channel lumen in the closed state. Two mechanisms of gating have been postulated. In one mechanism, the proton is "shuttled" through the gate by attaching to the delta nitrogen atom on the extracellular side of the imidazole ring, followed by the release of the proton attached to the epsilon nitrogen atom on the opposite side. In the second mechanism, the four histidines move away from each other due to electrostatic repulsion upon protonation, thus opening the gate sufficiently that a wire of water molecules can penetrate the gate. Then, protons are transported by "hopping" along the wire. In this paper, both mechanisms are evaluated in a series of molecular dynamics simulations by investigating stability of different protonation states of the channel that are involved in these mechanisms. For the shuttle mechanism, these are states with all epsilon protonated histidines, one biprotonated residue or one histidine protonated in the delta position. For the gate opening mechanism, this is the state in which all four histidines are biprotonated. In addition, a state with two biprotonated histidines is considered. For each system, composed of the protein channel embedded in phospholipid bilayer located between two water lamellae, a molecular dynamics trajectory of approximately 1.3 ns (after equilibration) was obtained. It is found that the states involved in the shuttle mechanism are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfers involved in the shuttle. In contract, the fully biprotonated state, implicated in the gate opening mechanism, is not stable and the channel looses its structural integrity. If only two histidines are biprotonated the channel deforms but remains intact with the gate mostly closed. In summary, the results of this study lend support to the shuttle mechanism but not to the gate opening mechanism of proton gating in M2.

Schweighofer, Karl J.↗
NASA NTRS1993

Molecular dynamics of interface rupture

Several situations have been studied in which a fluid-vapor or fluid-fluid interface ruptures, using molecular dynamics simulations of 3000 to 20,000 Lennard-Jones molecules in three dimensions. The cases studied are the Rayleigh instability of a liquid thread, the burst of a liquid drop immersed in a second liquid undergoing shear, and the rupture of a liquid sheet in an extensional flow. The late stages of the rupture process involve the gradual withdrawal of molecules from a thinning neck, or the appearance and growth of holes in a sheet. In all cases, it is found that despite the small size of the systems studied, tens of angstroms, the dynamics is in at least qualitative accord with the behavior expected from continuum calculations, and in some cases the agreement is to within tens of percent. Remarkably, this agreement occurs even though the Eulerian velocity and stress fields are essentially unmeasurable - dominated by thermal noise. The limitations and prospects for such molecular simulation techniques are assessed.

Koplik, Joel↗
NASA NTRS1992

Terraced spreading of simple liquids on solid surfaces

We have studied the spreading of liquid drops on a solid surface by molecular-dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor, and solid. We consider both spherically symmetric atoms and diatomic molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observe a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with substantial molecular diffusion both within and between layers. The quantitative behavior resembles recent experimental findings, but the detailed dynamics differ. In particular, the layers exhibit an unusual spreading law, where their radii vary in time as R-squared approximately equal to log10t, which disagrees with experiments on polymeric liquids as well as recent calculations.

Yang, Ju-Xing↗
NASA NTRS2016

Yield and Failure Behavior Investigated for Cross-Linked Phenolic Resins Using Molecular Dynamics

Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.

Molecular dynamics↗
NASA NTRS2016

Ab Initio Simulations of Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. We show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing phase transformation temperatures is discussed.

Simulations↗
NASA NTRSDate not supplied

Validation of Accurate Cryogenic Fluid Vapor-Liquid Boundary Conditions Via Molecular Simulations

Vapor-liquid interfaces drive many important phenomena in cryogenic fluid management, including heat transfer, evaporation, and capillary flow. Design of cryogenic fluid systems, such as propellant storage, requires accurate predictions of fluid behavior, including evaporation rates. Many models have been proposed for heat and mass transfer at vapor liquid interfaces, but the accuracy of these models in the context of cryogenic fluids has not been performed. We use molecular dynamics simulations, which allow for nanometer scale resolution of fluid phenomena, to evaluate the accuracy of a variety of vapor-liquid boundary conditions at evaporating and condensing interfaces. We find that an anisotropic temperature distribution is a critical ingredient for accurate prediction of intensive evaporation and condensation.

Daniel Vigil↗
NASA NTRS2001

Simulation of Peptides at Aqueous Interfaces

Behavior of peptides at water-membrane interfaces is of great interest in studies on cellular transport and signaling, membrane fusion, and the action of toxins and antibiotics. Many peptides, which exist in water only as random coils, can form sequence-dependent, ordered structures at aqueous interfaces, incorporate into membranes and self-assembly into functional units, such as simple ion channels. Multi -nanosecond molecular dynamics simulations have been carried out to study the mechanism and energetics of interfacial folding of both non-polar and amphiphilic peptides, their insertion into membranes and association into higher-order structures. The simulations indicate that peptides fold non-sequentially, often through a series of amphiphilic intermediates. They further incorporate into the membrane in a preferred direction as folded monomers, and only then aggregate into dimers and, possibly, further into "dimers of dimers".

Pohorille, Andrew↗
NASA NTRS2016

Yield and Failure Behavior Investigated for Cross-linked Phenolic Resins Using Molecular Dynamics

Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on the glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.

Stress-strain curves↗
NASA NTRS2016

Ab Initio Simulations of Temperature Dependent Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. In particular, we show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing these phase transformations is discussed.

nickel titanium↗