Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical molecular dynamics”

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 NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
NASA NTRS2002

Computational Nanomechanics of Carbon Nanotubes and Composites

Nanomechanics of individual carbon and boron-nitride nanotubes and their application as reinforcing fibers in polymer composites has been reviewed with interplay of theoretical modeling, computer simulations and experimental observations. The emphasis in this work is on elucidating the multi-length scales of the problems involved, and of different simulation techniques that are needed to address specific characteristics of individual nanotubes and nanotube polymer-matrix interfaces. Classical molecular dynamics simulations are shown to be sufficient to describe the generic behavior such as strength and stiffness modulus but are inadequate to describe elastic limit and nature of plastic buckling at large strength. Quantum molecular dynamics simulations are shown to bring out explicit atomic nature dependent behavior of these nanoscale materials objects that are not accessible either via continuum mechanics based descriptions or through classical molecular dynamics based simulations. As examples, we discus local plastic collapse of carbon nanotubes under axial compression and anisotropic plastic buckling of boron-nitride nanotubes. Dependence of the yield strain on the strain rate is addressed through temperature dependent simulations, a transition-state-theory based model of the strain as a function of strain rate and simulation temperature is presented, and in all cases extensive comparisons are made with experimental observations. Mechanical properties of nanotube-polymer composite materials are simulated with diverse nanotube-polymer interface structures (with van der Waals interaction). The atomistic mechanisms of the interface toughening for optimal load transfer through recycling, high-thermal expansion and diffusion coefficient composite formation above glass transition temperature, and enhancement of Young's modulus on addition of nanotubes to polymer are discussed and compared with experimental observations.

Srivastava, Deepak↗
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 NTRS2001

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗
NASA NTRS2020

Charged Particle Dynamics in the Lunar Environment

As the Artemis Program streams forward, organizations of scientists and engineers across the country have been coming together to solve the complex network of problems to once again achieve the milestone of successfully touching down on the moon. This project is no exception, and it has been an honor to work with the Electrostatics and Surface Physics Laboratory (ESPL), a lab within the Exploration Research and Technology Programs’ Spaceport Technologies Office (UB-G) located at the National Aeronautics and Space Administration at Kennedy Space Center (NASA KSC). The authors and their mentor James R. Phillips III3, in tandem with researchers from the Astrodynamics and Space Robotics Laboratory (ASRL) at the University of Central Florida (UCF), have been working on creating a state-of-the-art (SOA) granular gas dynamics model for particulate contamination prevention and mitigation purposes. In essence, the underlying objective for this project is to more accurately model the resulting electrodynamic interactions between lunar regolith grains with applications to dust mitigation and rocket engine plume surface interactions. To achieve this goal, the team has been expanding upon existing open source classical molecular dynamics code developed by Sandia National Laboratories (SNL). The following reports on the details of the problem at hand as well as the contributions that the authors have made towards resolution, including but not limited to the encoding of physical attributes and interactions for non-spherical polydisperse particle distributions within various bed geometries and preemptive data analysis implementations. Significant data analysis processes were utilized, and several original algorithms were created to perform critical evaluations, resulting in only a 0.006% error in discrepancy.

Annelisa B Esparza↗
NASA NTRS1999

Virtual Prototyping for Carbon Nanotubes and Their Compositions

Carbon nanotubes (CNT) are of interest as constituents-reinforcements for composites, probes for scanning microscopy, quantum wires for electronics, and building blocks for nanomachinery. Their small size makes direct experimental measurements difficult. Therefore, the critical answers about their properties should be obtained based on adequate modeling, virtual prototyping, using the methods of classical molecular dynamic (MD), in combination with the analytical methods and continuum mechanics.

Yakobson, Boris I.↗
NASA NTRS2001

Thermodynamic Behavior of Nano-sized Gold Clusters on the (001) Surface

We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Paik, Sun M.↗
NASA NTRS2000

Computational Nanotechnology of Materials, Devices, and Machines: Carbon Nanotubes

The mechanics and chemistry of carbon nanotubes have relevance for their numerous electronic applications. Mechanical deformations such as bending and twisting affect the nanotube's conductive properties, and at the same time they possess high strength and elasticity. Two principal techniques were utilized including the analysis of large scale classical molecular dynamics on a shared memory architecture machine and a quantum molecular dynamics methodology. In carbon based electronics, nanotubes are used as molecular wires with topological defects which are mediated through various means. Nanotubes can be connected to form junctions.

Srivastava, Deepak↗
NASA NTRS2001

Computational Nanotechnology of Materials, Electronics and Machines: Carbon Nanotubes

This report presents the goals and research of the Integrated Product Team (IPT) on Devices and Nanotechnology. NASA's needs for this technology are discussed and then related to the research focus of the team. The two areas of focus for technique development are: 1) large scale classical molecular dynamics on a shared memory architecture machine; and 2) quantum molecular dynamics methodology. The areas of focus for research are: 1) nanomechanics/materials; 2) carbon based electronics; 3) BxCyNz composite nanotubes and junctions; 4) nano mechano-electronics; and 5) nano mechano-chemistry.

Srivastava, Deepak↗
NASA NTRS2001

Thermal Expansion and Diffusion Coefficients of Carbon Nanotube-Polymer Composites

Classical molecular dynamics (MD) simulations employing Brenner potential for intra-nanotube interactions and van der Waals forces for polymer-nanotube interface have been used to investigate thermal expansion and diffusion characteristics of carbon nanotube-polyethylene composites. Addition of carbon nanotubes to polymer matrix is found to significantly increase the glass transition temperature Tg, and thermal expansion and diffusion coefficients in the composite above Tg. The increase has been attributed to the temperature dependent increase of the excluded volume for the polymer chains, and the findings could have implications in the composite processing, coating and painting applications.

Wei, Chengyu↗
NASA NTRS1998

Nanoscale Etching and Indentation of Silicon(001) Surface with Carbon Nanotube Tips

The possibility of nanoscale etching and indentation of Si(001)(2x1) surface by (8,0) and (10,10) carbon nanotube tips is demonstrated, for the first time, by classical molecular dynamics simulations employing Tersoff's many-body potential for a mixed C/Si/Ge system. In the nanotube tip barely touching the surface scenario atomistic etching is observed, where as in the nanoindentation scenario nanotube tip penetrates the surface without much hindrance. The results are explained in terms of the relative strength of C-C, C-Si, and Si-Si bonds.

Dzegilenko, Fendor N.↗
NASA NTRS1998

Nanoscale Etching and Indentation of Silicon Surfaces with Carbon Nanotubes

The possibility of nanolithography of silicon and germanium surfaces with bare carbon nanotube tips of scanning probe microscopy devices is considered with large scale classical molecular dynamics (MD) simulations employing Tersoff's reactive many-body potential for heteroatomic C/Si/Ge system. Lithography plays a key role in semiconductor manufacturing, and it is expected that future molecular and quantum electronic devices will be fabricated with nanolithographic and nanodeposition techniques. Carbon nanotubes, rolled up sheets of graphene made of carbon, are excellent candidates for use in nanolithography because they are extremely strong along axial direction and yet extremely elastic along radial direction. In the simulations, the interaction of a carbon nanotube tip with silicon surfaces is explored in two regimes. In the first scenario, the nanotubes barely touch the surface, while in the second they are pushed into the surface to make "nano holes". The first - gentle scenario mimics the nanotube-surface chemical reaction induced by the vertical mechanical manipulation of the nanotube. The second -digging - scenario intends to study the indentation profiles. The following results are reported in the two cases. In the first regime, depending on the surface impact site, two major outcomes outcomes are the selective removal of either a single surface atom or a surface dimer off the silicon surface. In the second regime, the indentation of a silicon substrate by the nanotube is observed. Upon the nanotube withdrawal, several surface silicon atoms are adsorbed at the tip of the nanotube causing significant rearrangements of atoms comprising the surface layer of the silicon substrate. The results are explained in terms of relative strength of C-C, C-Si, and Si-Si bonds. The proposed method is very robust and does not require applied voltage between the nanotube tips and the surface. The implications of the reported controllable etching and hole-creating for nanolithography on silicon are discussed in detail.

Dzegilenko, Fedor N.↗
NASA NTRS2015

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗
NASA NTRSDate not supplied

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗
NASA NTRS1998

Carbon Nanotube Based Molecular Electronics and Motors: A View from Classical and Quantum Dynamics Simulations

The tubular forms of fullerenes popularly known as carbon nanotubes are experimentally produced as single-, multiwall, and rope configurations. The nanotubes and nanoropes have shown to exhibit unusual mechanical and electronic properties. The single wall nanotubes exhibit both semiconducting and metallic behavior. In short undefected lengths they are the known strongest fibers which are unbreakable even when bent in half. Grown in ropes their tensile strength is approximately 100 times greater than steel at only one sixth the weight. Employing large scale classical and quantum molecular dynamics simulations we will explore the use of carbon nanotubes and carbon nanotube junctions in 2-, 3-, and 4-point molecular electronic device components, dynamic strength characterization for compressive, bending and torsional strains, and chemical functionalization for possible use in a nanoscale molecular motor. The above is an unclassified material produced for non-competitive basic research in the nanotechnology area.

Srivastava, Deepak↗
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 NTRS2020

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗