Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Nanomechanics of Carbon and CxByNz Nanotubes: Via a Quantum Molecular Dynamics Method

Nanomechanics of single-wall C, BN and BC$_3$ and B doped C nanotubes under axial compression and tension are investigated through a generalized tight-binding molecular dynamics (GTBMD) and {\it ab-initio} electronic structure methods. The dynamic strength of BN, BC$_3$ and B doped C nanotubes for small axial strain are comparable to each other. The main difference is in the critical strain at which structural collapse occurs. For example, even a shallow doping with B lowers the value of critical strain for C nanotubes. The critical strain for BN nanotube is found to be more than that for the similar C nanotube. Once the structural collapse starts to occur we find that carbon nanotubes irreversibly go into plastic deformation regime via the formation of tetrahedral (four-fold coordinated) bonds at the location of sharp pinches or kinks. This finding is considerably different from the classical MD (molecular dynamics) simulation results known so far. The energetics and electronic densities of states of the collapsed structures, investigated with {\it ab-initio) methods, will also be discussed.

Srivastava, Deepak↗
NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗
NASA NTRS2007

Local Elastic Constants for Epoxy-Nanotube Composites from Molecular Dynamics Simulation

A method from molecular dynamics simulation is developed for determining local elastic constants of an epoxy/nanotube composite. The local values of C11, C33, K12, and K13 elastic constants are calculated for an epoxy/nanotube composite as a function of radial distance from the nanotube. While the results possess a significant amount of statistical uncertainty resulting from both the numerical analysis and the molecular fluctuations during the simulation, the following observations can be made. If the size of the region around the nanotube is increased from shells of 1 to 6 in thickness, then the scatter in the data reduces enough to observe trends. All the elastic constants determined are at a minimum 20 from the center of the nanotube. The C11, C33, and K12 follow similar trends as a function of radial distance from the nanotube. The K13 decreases greater distances from the nanotube and becomes negative which may be a symptom of the statistical averaging.

Frankland, S. J. V.↗
NASA NTRS1991

Disordering and H-embrittlement of Pb(110) surface using embedded atom method and molecular dynamics

The dynamics and structure of the crystal melt interface are of great importance in crystal growth, melting, and other properties of the solid-liquid interface. The Embedded Atom Method (EAM), a semiempirical method developed by Daw and Baskes for determining the energetics of atoms in a bulk environment, is a useful method in the computer simulation. The EAM functions of PB in conjunction with the MD code are employed to study melting and disordering of the Pb(100) surface. The EAM functions are fitted to the bulk experimental data at zero temperature, and it is not clear if it will behave correctly at higher temperatures. In fact, it is shown that use of a similar type of potential at higher temperatures might result in a negative coefficient of thermal expansion. The primary purpose here is to perform MD simulation combined with the EAM functions of PB to study melting of the Pb(100) surface. This provides an accurate test of the EAM functions at higher temperatures. In particular, we calculate the two dimensional structure factor, the two dimensional distribution function, and density profile as a function of temperature up to the melting point of lead.

Karimi, Majid↗
NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
NASA NTRS1984

Models of Rapid Solidification

Laser annealing studies provide much information on various consequences of rapid solidification, including the trapping of impurities in the crystal, the generation of vacancies and twins, and on the fundamental limits to the speed of the crystal-melt interface. Some results obtained by molecular dynamics methods of the solidification of a Lennard-Jones liquid are reviewed. An indication of the relationship of interface speed to undercooling for certain materials can be derived from this model. Ising model simulations of impurity trapping in silicon are compared with some of the laser annealing results. The consequences of interface segregation and atomic strain are discussed.

Gilmer, G. H.↗
NASA NTRS2000

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing energy and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following questions about these proteins: (1) How do small proteins self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations identify the mechanisms by which protons move through the gate.

Pohorille, Andrew↗
NASA NTRS2001

Molecular Dynamics Studies of Structure and Functions of Water-Membrane Interfaces

A large number of essential cellular processes occur at the interfaces between water and membranes. The selectivity and dynamics of these processes are largely determined by the structural and electrical properties of the water-membrane interface. We investigate these properties by the molecular dynamics method. Over the time scales of the simulations, the membrane undergoes fluctuations described by the capillary wave model. These fluctuations produce occasional thinning defects in the membrane which provide effective pathways for passive transport of ions and small molecules across the membrane. Ions moving through the membrane markedly disrupt its structure and allow for significant water penetration into the membrane interior. Selectivity of transport, with respect to ionic charge, is determined by the interfacial electrostatic potential. Many small molecules. of potential significance in catalysis, bioenergetics and pharmacology, are shown to bind to the interface. The energetics and dynamics of this process will be discussed.

Pohorille, Andrew↗
NASA NTRS2003

Scalability of a Low-Cost Multi-Teraflop Linux Cluster for High-End Classical Atomistic and Quantum Mechanical Simulations

Scalability of a low-cost, Intel Xeon-based, multi-Teraflop Linux cluster is tested for two high-end scientific applications: Classical atomistic simulation based on the molecular dynamics method and quantum mechanical calculation based on the density functional theory. These scalable parallel applications use space-time multiresolution algorithms and feature computational-space decomposition, wavelet-based adaptive load balancing, and spacefilling-curve-based data compression for scalable I/O. Comparative performance tests are performed on a 1,024-processor Linux cluster and a conventional higher-end parallel supercomputer, 1,184-processor IBM SP4. The results show that the performance of the Linux cluster is comparable to that of the SP4. We also study various effects, such as the sharing of memory and L2 cache among processors, on the performance.

Kikuchi, Hideaki↗
NASA NTRS2003

Vacancy Mediated Mechanism of Nitrogen Substitution in Carbon Nanotubes

Nitrogen substitution reaction in a graphene sheet and carbon nanotubes of different diameter are investigated using the generalized tight-binding molecular dynamics method. The formation of a vacancy in curved graphene sheet or a carbon nanotube is found to cause a curvature dependent local reconstruction of the surface. Our simulations and analysis show that vacancy mediated N substitution (rather than N chemisorption) is favored on the surface of nanotubes with diameter larger than 8 nm. This predicted value of the critical minimum diameter for N incorporation is confirmed by experimental results presented.

Srivastava, Deepak↗
NASA NTRS1998

Structure of Boron Nitride Nanotubes: Tube Closing Vs. Chirality

The structure of boron nitride nanotubes is investigated using a generalized tight-binding molecular dynamics method. It is shown that dynamic relaxation results in a wavelike or "rippled" surface in which the B atoms rotate inward and the N atoms move outward, reminiscent of the surface relaxation of the III-V semiconductors. More importantly, the three different morphologies of the tube closing with flat, conical and amorphous ends, as observed in experiments, are shown to be directly related to the tube chiralities. The abundance of flat end tubes observed in experiments is, thus, shown to be an indication of the greater stability of "zig-zag" BN tubes over the "arm-chair" tubes under experimental conditions.

Srivastava, Deepak↗
NASA NTRS2000

Computer Modeling of the Earliest Cellular Structures and Functions

In the absence of extinct or extant record of protocells (the earliest ancestors of contemporary cells). the most direct way to test our understanding of the origin of cellular life is to construct laboratory models of protocells. Such efforts are currently underway in the NASA Astrobiology Program. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs for molecules that perform proto-cellular functions. Many of these functions, such as import of nutrients, capture and storage of energy. and response to changes in the environment are carried out by proteins bound to membrane< We will discuss a series of large-scale, molecular-level computer simulations which demonstrate (a) how small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes, (b) how these peptides aggregate to form membrane-spanning structures (eg. channels), and (c) by what mechanisms such aggregates perform essential proto-cellular functions, such as proton transport of protons across cell walls, a key step in cellular bioenergetics. The simulations were performed using the molecular dynamics method, in which Newton's equations of motion for each item in the system are solved iteratively. The problems of interest required simulations on multi-nanosecond time scales, which corresponded to 10(exp 6)-10(exp 8) time steps.

Pohorille, Andrew↗
NASA NTRS2006

Computer Modeling of Protocellular Functions: Peptide Insertion in Membranes

Lipid vesicles became the precursors to protocells by acquiring the capabilities needed to survive and reproduce. These include transport of ions, nutrients and waste products across cell walls and capture of energy and its conversion into a chemically usable form. In modem organisms these functions are carried out by membrane-bound proteins (about 30% of the genome codes for this kind of proteins). A number of properties of alpha-helical peptides suggest that their associations are excellent candidates for protobiological precursors of proteins. In particular, some simple a-helical peptides can aggregate spontaneously and form functional channels. This process can be described conceptually by a three-step thermodynamic cycle: 1 - folding of helices at the water-membrane interface, 2 - helix insertion into the lipid bilayer and 3 - specific interactions of these helices that result in functional tertiary structures. Although a crucial step, helix insertion has not been adequately studied because of the insolubility and aggregation of hydrophobic peptides. In this work, we use computer simulation methods (Molecular Dynamics) to characterize the energetics of helix insertion and we discuss its importance in an evolutionary context. Specifically, helices could self-assemble only if their interactions were sufficiently strong to compensate the unfavorable Free Energy of insertion of individual helices into membranes, providing a selection mechanism for protobiological evolution.

Rodriquez-Gomez, D.↗
NASA NTRS2008

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗
NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗
NASA NTRS2013

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗
NASA NTRS2011

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗