Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Dynamics of charged particles in a Paul radio-frequency quadrupole trap

A molecular-dynamics simulation of hundreds of ions confined in a Paul trap has been performed. The simulation includes the trapped particles' micromotion and interparticle Coulomb interactions. A random walk in velocity was implemented to bring the secular motion to a given temperature which was numerically measured. When the coupling Gamma is large the ions from concentric shells which undergo a quadrupole oscillation at the RF frequency, while the ions within a shell form a 2D hexagonal lattice. Ion clouds at 5 mK show no RF heating for q(z) less than about 0.6, whereas rapid heating is seen for qz = 0.8.

Prestage, J. D.↗
NASA NTRS1993

Calculation of molecular excitation rates

State-to-state collisional excitation rates for interstellar molecules observed by radio astronomers continue to be required to interpret observed line intensities in terms of local temperatures and densities. A problem of particular interest is collisional excitation of water which is important for modeling the observed interstellar masers. In earlier work supported by a different NASA Grant, excitation of water in collisions with He atoms was studied; after many years of successively more refined calculations that problem now seems to be well understood, and discrepancies with earlier experimental data for related (pressure broadening) phenomena are believed to reflect experimental errors. Because of interstellar abundances, excitation by H2, the dominant interstellar species, is much more important than excitation by He, although it has been argued that rates for excitation by these are similar. Under the current grant theoretical study of this problem has begun which is greatly complicated by the additional degrees of freedom which must be included both in determining the interaction potential and also in the molecular scattering calculation. We have now computed the interaction forces for nearly a thousand molecular geometries and are close to having an acceptable global fit to these points which is necessary for the molecular dynamics calculations. Also, extensive modifications have been made to the molecular scattering code, MOLSCAT. These included coding the rotational basis sets and coupling matrix elements required for collisions of an asymmetric top with a linear rotor. A new method for numerical solution of the coupled equations has been incorporated. Because of the long-ranged nature of the water-hydrogen interaction it is necessary to integrate the equations to rather large intermolecular separations, and the integration methods previously available in MOLSCAT are not ideal for such cases. However, the method used by Alexander in his HIBRIDON code is particularly suited for such cases. We have obtained this code and incorporated that part which solves the coupled differential equations as an option in the MOLSCAT program.

Flynn, George↗
NASA NTRS2000

Modeling Molecules

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

Source record↗
NASA NTRS1994

Equilibrium and Dynamics Properties of Poly(oxyethylene) Melts and Related Poly(alkylethers) from Simulations and Ab Initio Calculations

Molecular dynamics simulations of POE melts have been performed utilizing a potential force field parameterized to reproduce conformer energies and rotational energy barriers in dimethoxyethane as determined from ab initio electronic structure calculations. Chain conformations and dimensions of POE from the simulations were found to be in good agreement with predictions of a rotational isomeric state (RIS) model based upon the ab initio conformational. energies. The melt chains were found to be somewhat extended relative to chains at theta conditions. This effect will be discussed in light of neutron scattering experiments which indicate that POE chains are extended in the melt relative to theta solutions. The conformational characteristics of POE chains will also be compared with those of other poly (alkylethers), namely poly(oxymethylene), poly(oxytrimethylene) and poly(oxytetramethylene). Local conformational dynamics were found to be more rapid than in polymethylene. Calculated C-H vector correlation times were found to be in reasonable agreement with experimental values from C-13 NMR spin-lattice relaxation times. The influence of ionic salts on local conformations and dynamics will also be discussed.

Smith, Grant D.↗
NASA NTRS1972

Dynamics of the molecular and atomic mechanisms for the hydrogen-iodine exchange reaction.

The molecular and atomic mechanisms for the hydrogen-iodine exchange reaction are treated theoretically by means of extensive classical trajectories calculated on a reasonable potential energy surface on which the single adjustable parameter is the iodine-core effective charge. The analysis shows the molecular mechanism to be dynamically forbidden, but gives an over-all rate constant for the atomic mechanism that is in agreement with the experimental values. It is indicated that the formation of a weak H2I complex plays an important dynamical role if the atomic mechanism is limited to reactions with collision complexes involving no more than two hydrogen atoms and two iodine atoms. Excellent agreement with experiment is obtained for the rate constant for the recombination I+I+H2 yields I2+H2 and its negative temperature coefficient.

Raff, L. M.↗
NASA NTRS2011

Ndarts

Ndarts software provides algorithms for computing quantities associated with the dynamics of articulated, rigid-link, multibody systems. It is designed as a general-purpose dynamics library that can be used for the modeling of robotic platforms, space vehicles, molecular dynamics, and other such applications. The architecture and algorithms in Ndarts are based on the Spatial Operator Algebra (SOA) theory for computational multibody and robot dynamics developed at JPL. It uses minimal, internal coordinate models. The algorithms are low-order, recursive scatter/ gather algorithms. In comparison with the earlier Darts++ software, this version has a more general and cleaner design needed to support a larger class of computational dynamics needs. It includes a frames infrastructure, allows algorithms to operate on subgraphs of the system, and implements lazy and deferred computation for better efficiency. Dynamics modeling modules such as Ndarts are core building blocks of control and simulation software for space, robotic, mechanism, bio-molecular, and material systems modeling.

Jain, Abhinandan↗
NASA NTRS1998

Computer Simulation of the far Infrared Collision Induced Absorption Spectra of Gaseous CO2

Far infrared collision induced absorption spectra of gaseous CO2 were computed using molecular dynamics simulations. The quadrupole and hexadecapole multipolar induction, through the trace, and the anisotropy of the molecular polarizability were found to be insufficient to represent properly the dipole induction mechanism. For a detailed analysis of the induction process the spectra obtained were decomposed into components resulting from different terms of the induced dipole. Based on this decomposition, all additional overlap contribution for each term was proposed. When spectra were recomputed including such overlap, good agreement between experiment and simulation was achieved over the temperature range at which measurements exist (233-400 K). The use of an anisotropic intermolecular potential was found to be of critical importance for obtaining the right shape of the far wings of the spectra.

Gruszka, Marcin↗
NASA NTRS2016

Computational and Experimental Study of Li-doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Psi). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Psi. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Psi-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Psi and suggest a source of impedance for Li+ to reach the electrode.

Haskins, Justin B.↗
NASA NTRS2010

Interactions of Ku70/80 with Double-Strand DNA: Energetic, Dynamics, and Functional Implications

Space radiation is a proficient inducer of DNA damage leading to mutation, aberrant cell signaling, and cancer formation. Ku is among the first responding proteins in nucleus to recognize and bind the DNA double strand breaks (DSBs) whenever they are introduced. Once loaded Ku works as a scaffold to recruit other repair factors of non-homologous end joining and facilitates the following repair processes. The crystallographic study of the Ku70/80 heterodimer indicate the core structure of this protein shows virtually no conformational change after binding with DNA. To investigate the dynamical features as well as the energetic characteristics of Ku-DNA binding, we conduct multi-nanosecond molecular dynamics simulations of a modeled Ku70/80 structure and several complexes with two 24-bp DNA duplexes. Free energy calculations show significant energy differences between the complexes with Ku bound at DSBs and those with Ku associated at an internal site of a chromosome. The results also reveal detailed interactions between different nucleotides and the amino acids along the DNA-binding cradle of Ku, indicating subtle binding preference of Ku at specific DNA sequences. The covariance matrix analyses along the trajectories demonstrate the protein is stimulated to undergo correlated motions of different domains once bound to DNA ends. Additionally, principle component analyses identify these low frequency collective motions suitable for binding with and translocation along duplex DNA. It is proposed that the modification of dynamical properties of Ku upon binding with DSBs may provide a signal for the further recruitment of other repair factors such as DNA-PKcs, XLF, and XRCC4.

Hu, Shaowen↗
NASA NTRS2003

Atomic-scale Modeling of the Structure and Dynamics of Dislocations in Complex Alloys at High Temperatures

We report on the progress made during the first year of the project. Most of the progress at this point has been on the theoretical and computational side. Here are the highlights: (1) A new code, tailored for high-end desktop computing, now combines modern Accelerated Dynamics (AD) with the well-tested Embedded Atom Method (EAM); (2) The new Accelerated Dynamics allows the study of relatively slow, thermally-activated processes, such as diffusion, which are much too slow for traditional Molecular Dynamics; (3) We have benchmarked the new AD code on a rather simple and well-known process: vacancy diffusion in copper; and (4) We have begun application of the AD code to the diffusion of vacancies in ordered intermetallics.

Daw, Murray S.↗
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 NTRS1990

DNA dynamics in aqueous solution: opening the double helix

The opening of a DNA base pair is a simple reaction that is a prerequisite for replication, transcription, and other vital biological functions. Understanding the molecular mechanisms of biological reactions is crucial for predicting and, ultimately, controlling them. Realistic computer simulations of the reactions can provide the needed understanding. To model even the simplest reaction in aqueous solution requires hundreds of hours of supercomputing time. We have used molecular dynamics techniques to simulate fraying of the ends of a six base pair double strand of DNA, [TCGCGA]2, where the four bases of DNA are denoted by T (thymine), C (cytosine), G (guanine), and A (adenine), and to estimate the free energy barrier to this process. The calculations, in which the DNA was surrounded by 2,594 water molecules, required 50 hours of CRAY-2 CPU time for every simulated 100 picoseconds. A free energy barrier to fraying, which is mainly characterized by the movement of adenine away from thymine into aqueous environment, was estimated to be 4 kcal/mol. Another fraying pathway, which leads to stacking between terminal adenine and thymine, was also observed. These detailed pictures of the motions and energetics of DNA base pair opening in water are a first step toward understanding how DNA will interact with any molecule.

NASA Discipline Number 52-20↗
NASA NTRS1999

Energy Minimization of Molecular Features Observed on the (110) Face of Lysozyme Crystals

Molecular dynamics and energy minimization have been carried out using the program XPLOR to check the plausibility of a model lysozyme crystal surface. The molecular features of the (110) face of lysozyme were observed using atomic force microscopy (AFM). A model of the crystal surface was constructed using the PDB file 193L, and was used to simulate an AFM image. Molecule translations, van der Waals radii, and assumed AFM tip shape were adjusted to maximize the correlation coefficient between the experimental and simulated images. The highest degree of 0 correlation (0.92) was obtained with the molecules displaced over 6 A from their positions within the bulk of the crystal. The quality of this starting model, the extent of energy minimization, and the correlation coefficient between the final model and the experimental data will be discussed.

Perozzo, Mary A.↗
NASA NTRS2003

A Dynamic Pathway for Stone-Wales Bond Rotation on Carbon Nanotubes through Diamond-Like Bonds

A new lower energy barrier with a two-step pathway of Stone-Wales (SW) ,ond rotation on carbon nanotubes (CNTs) is found through molecular dynamics (MD) simulations of CNTs under tension. The first step involves going over to a stable sp3-like metastable configuration with half rotated and partially tilted C-C bond. The second step involves going over to the fully rotated C-C bond with the formation of a SW defect in the nanotube. The energy barrier for this two-step dynamic pathway is significantly lower than the previously known static barrier for in-plane rotation of the C-C bond on a tensile strained (> 4%) CNT.

Wei, Chen-Yu↗
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 NTRS2001

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗
NASA NTRS1993

Point charge representation of multicenter multipole moments in calculation of electrostatic properties

Distributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multi-center multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.

Non-NASA Center↗
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↗