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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS2002

On the Boundary Condition for Water at a Hydrophobic, Dense Surface

We study the no-slip boundary conditions for water at a hydrophobic (graphite) surface using non-equilibrium molecular-dynamics simulations. For the planar Couette flow, we find a slip length of 64 nm at 1 bar and 300 K, decreasing with increasing system pressure to a value of 31 nm at 1000 bar. Changing the properties of the interface to from hydrophobic to strongly hydrophilic reduces the slip to 14 nm. Finally, we study the flow of water past an array of carbon nanotubes mounted in an inline configuration with a spacing of 16.4 x 16.4 nm. For tube diameters of 1.25 and 2.50 nm we find drag coefficients in good agreement with the macroscopic, Navier-Stokes values. For carbon nanotubes, the no-slip condition is valid to within the definition of the position of the interface.

Walther, J. H.↗

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

Theory of Transport of Long Polymer Molecules through Carbon Nanotube Channels

A theory of transport of long chain polymer molecules through carbon nanotube (CNT) channels is developed using Fokker-Planck equation and direct molecular dynamics (MD) simulations. The mean transport or translocation time tau is found to depend on the chemical potential energy, entropy and diffusion coefficient. A power law dependence tau approx. N(sup 2)is found where N is number of monomers in a molecule. For 10(exp 5)-unit long polyethylene molecules, tau is estimated to be approx. 1micro-s. The diffusion coefficient of long polymer molecules inside CNTs, like that of short ones, are found to be few orders of magnitude larger than in ordinary silicate based zeolite systems.

Wei, Chenyu↗

NASA NTRS2003

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

NASA NTRS2003

Discrimination of Single Base Pair Differences Among Individual DNA Molecules Using a Nanopore

The protein toxin alpha-hemolysin form nanometer scale channels across lipid membranes. Our lab uses a single channel in an artificial lipid bilayer in a patch clamp device to capture and examine individual DNA molecules. This nanopore detector used with a support vector machine (SVM) can analyze DNA hairpin molecules on the millisecond time scale. We distinguish duplex stem length, base pair mismatches, loop length, and single base pair differences. The residual current fluxes also reveal structural molecular dynamics elements. DNA end-fraying (terminal base pair dissociation) can be observed as near full blockades, or spikes, in current. This technique can be used to investigate other biological processes dependent on DNA end-fraying, such as the processing of HIV DNA by HIV integrase.

Vercoutere, Wenonah↗

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 NTRS2003

Detection of Strand Cleavage And Oxidation Damage Using Model DNA Molecules Captured in a Nanoscale Pore

We use a biological nano-scale pore to distinguish among individual DNA hairpins that differ by a single site of oxidation or a nick in the sugar-phosphate backbone. In earlier work we showed that the protein ion channel alpha-hemolysin can be used as a detector to distinguish single-stranded from double-stranded DNA, single base pair and single nucleotide differences. This resolution is in part a result of sensitivity to structural changes that influence the molecular dynamics of nucleotides within DNA. The strand cleavage products we examined here included a 5-base-pair (5-bp) hairpin with a 5-prime five-nucleotide overhang, and a complementary five-nucleotide oligomer. These produced predictable shoulder-spike and rapid near-full blockade signatures, respectively. When combined, strand annealing was monitored in real time. The residual current level dropped to a lower discrete level in the shoulder-spike blockade signatures, and the duration lengthened. However, these blockade signatures had a shorter duration than the unmodified l0bp hairpin. To test the pore sensitivity to nucleotide oxidation, we examined a 9-bp hairpin with a terminal 8-oxo-deoxyguanosine (8-oxo-dG), or a penultimate 8-oxo-dG. Each produced blockade signatures that differed from the otherwise identical control 9bp hairpins. This study showed that DNA structure is modified sufficiently by strand cleavage or oxidation damage at a single site to alter in a predictable manner the ionic current blockade signatures produced. This technique improves the ability to assess damage to DNA, and can provide a simple means to help characterize the risks of radiation exposure. It may also provide a method to test radiation protection.

Vercoutere, W.↗

NASA NTRS2004

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗

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 NTRS2003

Conductance Oscillations in Squashed Carbon Nanotubes

A combination of molecular dynamics and electrical conductance calculations are used to probe the electromechanical properties of squashed metallic carbon nanotubes. We find that the conductance and bandgap of armchair nanotubes show oscillations upon squashing. The physical origin of these oscillations is attributed to interaction of carbon atoms with a fourth neighbor. Squashing of armchair and zigzag nanotubes ultimately leads to metallic behavior.

Mehrez, H.↗

NASA NTRS2003

Structure, Function, Self-Assembly and Origin of Simple Membrane Proteins

Integral membrane proteins perform such essential cellular functions as transport of ions, nutrients and waste products across cell walls, transduction of environmental signals, regulation of cell fusion, recognition of other cells, energy capture and its conversion into high-energy compounds. In fact, 30-40% of genes in modem organisms codes for membrane proteins. Although contemporary membrane proteins or their functional assemblies can be quite complex, their transmembrane fragments are usually remarkably simple. The most common structural motif for these fragments is a bundle of alpha-helices, but occasionally it could be a beta-barrel. In a series of molecular dynamics computer simulations we investigated self-organizing properties of simple membrane proteins based on these structural motifs. Specifically, we studied folding and insertion into membranes of short, nonpolar or amphiphatic peptides. We also investigated glycophorin A, a peptide that forms sequence-specific dimers, and a transmembrane aggregate of four identical alpha-helices that forms an efficient and selective voltage-gated proton channel was investigated. Many peptides are attracted to water-membrane interfaces. Once at the interface, nonpolar peptides spontaneously fold to a-helices. Whenever the sequence permits, peptides that contain both polar and nonpolar amino also adopt helical structures, in which polar and nonpolar amino acid side chains are immersed in water and membrane, respectively. Specific identity of side chains is less important. Helical peptides at the interface could insert into the membrane and adopt a transmembrane conformation. However, insertion of a single helix is unfavorable because polar groups in the peptide become completely dehydrated upon insertion. The unfavorable free energy of insertion can be regained by spontaneous association of peptides in the membrane. The first step in this process is the formation of dimers, although the most common are aggregates of 4-7 helices. The helices could arrange themselves such that they formed pores capable of transporting ions and small molecules across membranes. Stability of transmembrane aggregates of simple proteins is often only marginal and, therefore, it can be regulated by environmental signals or small sequence modifications in the region of interhelical interactions. A key step in the earliest evolution of membrane proteins was the emergence of selectivity for specific substrates. Many channels could become selective if one or only a few properly chosen amino acids are properly placed along the channel, acting as filters or gates. This is a convenient evolutionary solution because it does not require imposing conditions on the whole sequence.

Pohorille, Andrew↗

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 NTRS1995

Scattering of water from the glycerol liquid-vacuum interface

Molecular dynamics calculations of the scattering of D2O from the glycerol surface at different collision energies are reported. The results for the trapping probabilities and energy transfer are in good agreement with experiments. The calculations demonstrate that the strong attractive forces between these two strongly hydrogen bonding molecules have only a minor effect on the initial collision dynamics. The trapping probability is influenced to a significant extent by the repulsive hard sphere-like initial encounter with the corrugated surface and, only at a later stage, by the efficiency of energy flow in the multiple interactions between the water and the surface molecules.

NASA Center ARC↗

NASA NTRS2004

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗

NASA NTRS2005

Applications of Quantum Chemistry to the Study of Carbon Nanotubes

For several years, scientists at NASA Ames have been studying the properties of carbon nanotubes using various experimental and computational methods. In this talk, I will compare different strategies for using quantum chemistry calculations to describe the electronic structure, deformation and chemical functionalization of single wall carbon nanotubes (SWNT) and the physisorption of small molecules on nanotube surfaces. The SWNT can be treated as an infinite (periodic) or finite length carbon cylinder or as a polycyclic aromatic hydrocarbon (PAH) molecule with an imposed curvature maintained by external constraints (as if it were cut out of the SWNT surface). Calculations are carried out using DFT and MP2 methods and a variety of atomic orbital basis sets from minimal (STO-3G) to valence triple zeta. The optimal approach is based on the particular SWNT property of interest. Examples to be discussed include: nanotube fluorination and other functionalization reactions; coating of nanotubes by water vapor and low-molecular weight organic molecules; and the nature of the interface between SWNT and liquids such as water and amines. In many cases, the quantum chemistry calculations are used to parameterize or validate force fields for molecular dynamics simulations. The results of these calculations have helped explain experimental data and contributed to the design of novel materials and sensors based on carbon nanotubes. Some of this research is described in the following papers:

Jaffe, Richard L.↗

NASA NTRS2004

Modeling the interaction Between Ethylene Diamine and Water Films on the Surface of a Carbon Nanotube

It has been observed that a carbon nanotube (CNT) AFM tip coated with ethylene diamine (EDA) penetrates the liquid water-air interface more easily than an uncoated nanotube tip. The EDA coating remains intact through repeated cycles of dipping and removal. In order to understand the physical basis for this observation, we use ab initio quantum chemistry calculations to study the EDA-CNT-water interaction and to parameterize a force field describing this system. Molecular dynamics (MD) simulations are carried out for EDA-water mixtures and an EDA-coated carbon nanotube immmed in water. These simulations are similar to our earlier MD study that characterized the CNT-water interface. The attractive CNT-EDA and CNT-water interactions arise primarily from van der Waals forces, and the EDA-EDA, EDA-water and water-water interactions are mainly due to hydrogen bond formation. The binding energ of single EDA molecule to the nanotube is nearly three times larger than the corresponding value found for water (4.3 versus 1.5 kcal mol, respectively). The EDA molecules readily stick to and diffuse along the CNT surface. As a resulf mixing of the EDA and water films does not occur on the timescale of the MD simulations. The EDA film reduces the hydrophobicity of the nanotube surface and acts like a prototypical surfactant in stabilizing the suspension of carbon nanotubes in water. For this presentation, we use the MD simulations to determine how the presence of the carbon nanotube surface perturbs the properties of EDA-water mixtures.

Jaffe, Richard L.↗

NASA NTRS2005

Prediction of Mechanical Properties of Polymers With Various Force Fields

The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide.

Odegard, Gregory M.↗

NASA NTRS2005

Multi-Scale Modeling of Cross-Linked Nanotube Materials

The effect of cross-linking single-walled carbon nanotubes on the Young's modulus of a nanotube-reinforced composite is modeled with a multi-scale method. The Young's modulus is predicted as a function of nanotube volume fraction and cross-link density. In this method, the constitutive properties of molecular representative volume elements are determined using molecular dynamics simulation and equivalent-continuum modeling. The Young's modulus is subsequently calculated for cross-linked nanotubes in a matrix which consists of the unreacted cross-linking agent. Two different cross-linking agents are used in this study, one that is short and rigid (Molecule A), and one that is long and flexible (Molecule B). Direct comparisons between the predicted elastic constants are made for the models in which the nanotubes are either covalently bonded or not chemically bonded to the cross-linking agent. At a nanotube volume fraction of 10%, the Young's modulus of Material A is not affected by nanotube crosslinking, while the Young's modulus of Material B is reduced by 64% when the nanotubes are cross-linked relative to the non-cross-linked material with the same matrix.

Frankland, S. J. V.↗