Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Peptide Folding and Translocation Across the Water-Membrane Interface

The ability of small peptides to organize at aqueous interfaces was examined by performing a series of large-scale, molecular dynamics computer simulations of several peptides composed of two amino acids, nonpolar leucine (L) and polar glutamine (Q). The peptides differed in size and sequence of the amino acids. Studies on dipeptides LL, LQ, QL and QQ were extended to two heptamers, LQQLLQL and LQLQLQL, designed to maximize interfacial stability of an alpha-helix and a beta-strand, respectively, by exposing polar side chains to water and nonpolar side chains to a nonpolar phase. Finally, a transition of an undecamer, composed entirely of leucine residues, from a disordered structure in water to an alpha-helix in a nonpolar phase representing the interior of the membrane was investigated. Complete folding of a peptide in solution was accomplished for the first time in computer simulations. The simulations revealed several basic principles governing the sequence-dependent organization of peptides at interfaces. Short peptides tend to accumulate at interfaces and acquire ordered structures, providing that they have a proper sequence of polar and nonpolar amino acids. The dominant factor determining the interfacial structure of peptides is the hydrophobic effect, which is manifested at aqueous interfaces as a tendency for polar and nonpolar groups of the solute to segregate into the aqueous and nonpolar phases, respectively. If peptides consist of nonpolar residue's only, they become inserted into the nonpolar phase. As demonstrated by the example of the leucine undecamer, such peptides fold into an alpha-helix as they partition into the nonpolar medium. The folding proceeds through an intermediate, called 3-10-helix, which remains in equilibrium with the alpha-helix. Once in the nonpolar environment, the peptides can readily change their orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with changing external conditions may have provided a simple mechanism of transmitting signals from the environment to the interior of a cell.

Pohorille, Andrew↗
NASA NTRS2002

Analysis of Carbon Nanotube Pull-out from a Polymer Matrix

Molecular dynamics (MD) simulations of carbon nanotube (NT) pull-out from a polymer matrix are carried out. As the NT pull-out develops in the simulation, variations in the displacement and velocities of the NT are monitored. The existence of a carbon-ring-based period in NT sliding during pull-out is identified. Linear trends in the NT velocity-force relation are observed and used to estimate an effective viscosity coefficient for interfacial sliding at the NT/polymer interface. As a result, the entire process of NT pull-out is characterized by an interfacial friction model that is based on a critical pull-out force, and an analog of Newton's friction law used to describe the NT/polymer interfacial sliding.

Frankland, S. J. V.↗
NASA NTRS2003

Structural Basis for "Flip-Flop" Action of Human Pyruvate Dehydrogenase

The derivative of vitamin B1, thiamin pyrophosphate is a cofactor of pyruvate dehydrogenase, a component enzyme of the mitochondrial pyruvate dehydrogenase multienzyme complex that plays a major role in directing energy metabolism in the cell. This cofactor is used to cleave the C(sup alpha)-C(=O) bond of pyruvate followed by reductive acetyl transfer to lipoyl-dihydrolipoamide acetyltransferase. In alpha(sub 2)beta(sub 2)-tetrameric human pyruvate dehydrogenase, there are two cofactor binding sites, each of them being a center of independently conducted, although highly coordinated enzymatic reactions. The dynamic nonequivalence of two, otherwise chemically equivalent, catalytic sites can now be understood based on the recently determined crystal structure of the holo-form of human pyruvate dehydrogenase at 1.95A resolution. The structure of pyruvate dehydrogenase was determined using a combination of MAD phasing and molecular replacement followed by rounds of torsion-angles molecular-dynamics simulated-annealing refinement. The final pyruvate dehydrogenase structure included coordinates for all protein amino acids two cofactor molecules, two magnesium and two potassium ions, and 742 water molecules. The structure was refined to R = 0.202 and R(sub free) = 0.244. Our structural analysis of the enzyme folding and domain assembly identified a simple mechanism of this protein motion required for the conduct of catalytic action.

Ciszak, Ewa↗
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

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 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 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 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 NTRS2007

Interactions of the C-terminal Domain of Human Ku70 with DNA Substrate: A Molecular Dynamics Study

NASA is developing a systems biology approach to improve the assessment of health risks associated with space radiation. The primary toxic and mutagenic lesion following radiation exposure is the DNA double strand break (DSB), thus a model incorporating proteins and pathways important in response and repair of this lesion is critical. One key protein heterodimer for systems models of radiation effects is the Ku(sub 70/80) complex. The Ku70/80 complex is important in the initial binding of DSB ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. The C-terminal domain of Ku70 (Ku70c, residues 559-609), contains an helix-extended strand-helix motif and similar motifs have been found in other nucleic acid-binding proteins critical for DNA repair. However, the exact mechanism of damage recognition and substrate specificity for the Ku heterodimer remains unclear in part due to the absence of a high-resolution structure of the Ku70c/DNA complex. We performed a series of molecular dynamics (MD) simulations on a system with the subunit Ku70c and a 14 base pairs DNA duplex, whose starting structures are designed to be variable so as to mimic their different binding modes. By analyzing conformational changes and energetic properties of the complex during MD simulations, we found that interactions are preferred at DNA ends, and within the major groove, which is consistent with previous experimental investigations. In addition, the results indicate that cooperation of Ku70c with other subunits of Ku(sub 70/80) is necessary to explain the high affinity of binding as observed in experiments.

Hu, Shaowen↗
NASA NTRS2007

Interactions of the SAP Domain of Human Ku70 with DNA Substrate: A Molecular Dynamics Study

NASA is developing a systems biology approach to improve the assessment of health risks associated with space radiation. The primary toxic and mutagenic lesion following radiation exposure is the DNA double strand break (DSB), thus a model incorporating proteins and pathways important in response and repair of this lesion is critical. One key protein heterodimer for systems models of radiation effects is the Ku70/80 complex. The Ku70/80 complex is important in the initial binding of DSB ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. The SAP domain of Ku70 (residues 556-609), contains an a helix-extended strand-helix motif and similar motifs have been found in other nucleic acid-binding proteins critical for DNA repair. However, the exact mechanism of damage recognition and substrate specificity for the Ku heterodimer remains unclear in part due to the absence of a high-resolution structure of the SAP/DNA complex. We performed a series of molecular dynamics (MD) simulations on a system with the SAP domain of Ku70 and a 10 base pairs DNA duplex. Large-scale conformational changes were observed and some putative binding modes were suggested based on energetic analysis. These modes are consistent with previous experimental investigations. In addition, the results indicate that cooperation of SAP with other domains of Ku70/80 is necessary to explain the high affinity of binding as observed in experiments.

Hu, Shaowen↗
NASA NTRS2008

A New Concurrent Multiscale Methodology 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.

Yamakov, Vesselin↗
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 NTRS2009

Filtering Water by Use of Ultrasonically Vibrated Nanotubes

Devices that could be characterized as acoustically driven molecular sieves have been proposed for filtering water to remove all biological contaminants and all molecules larger than water molecules. Originally intended for purifying wastewater for reuse aboard spacecraft, these devices could also be attractive for use on Earth in numerous settings in which there are requirements to obtain potable, medical-grade, or otherwise pure water from contaminated water supplies. These devices could also serve as efficient means of removing some or all water from chemical products . for example, they might be useful as adjuncts or substitutes for stills in the removal of water from alcohols and alcoholic beverages. These devices may be constructed using various materials, such as ceramics, metallics, or polymers, depending on end-use requirements. A representative device of this type (see figure) would include a polymeric disk, about 1 mm in diameter and between 1 and 40 microns thick, within which would be embedded single-wall carbon nanotubes aligned along the thickness axis. The polymeric disk would be part of a unitary polymeric ring assembly. An acoustic transducer in the form of a piezoelectric-film-and-electrode subassembly - typically 9 microns thick and made of poly(vinylidene fluoride) coated with copper 150 nm thick -. would be affixed to the outside of the outer polymeric ring by means of an electrically nonconductive epoxy. The nanotubes would be chosen to have diameters between about 8 and about 13.5 A because water molecules could fit into the nanotubes, but larger molecules could not. Water to be purified would be placed in contact with one face (typically, the upper face) of the filter disk. The surface tension of water is low enough that water molecules should enter and travel along the nanotubes, and computational simulations of molecular dynamics and experimental measurements have shown that the water molecules inside the nanotubes in this size range can be expected to become aligned into helical columns that exhibit properties of both hexagonal ice crystals and liquid water

Gavalas, Lillian Susan↗
NASA NTRS2014

Multiscale Analysis of Structurally-Graded Microstructures Using Molecular Dynamics, Discrete Dislocation Dynamics and Continuum Crystal Plasticity

A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.

Saether, Erik↗
NASA NTRS2019

Emergence of Simple Ion Channels by Self-Assembly of Short Peptides

As protocells transitioned from structures bound by membranes made of simple amphiphiles, such as fatty acids, to structures containing phospholipid membranes, they became nearly impermeable to charged species. Yet, the ability to transport ions across membranes was vital to regulating cellular volume, pH homeostasis, generating energy and sensing the environment. For this reason, life evolved ion channels, protein structures surrounding water-filled pores in the membrane that facilitated ion transport. In contemporary cells, ion channels are complex protein assemblies. Could the same function be carried out, even if less efficiently and selectively, by much simpler peptides at the origins of life? To answer this question, we turned to trichogen (GAIV), an 11-residue peptide that was found to permeabilize membrane, causing leakage of the vesicle content including ions. Since GAIV is synthesized non-gnomically, it represents as excellent model for the earliest, poorly optimized channels. Understanding the mechanism by which short peptides, such as GAIV, mediate ion transport across membranes provides crucial information on the emergence, structure and the mechanism of action of the earliest ion channels. In this study, we used extensive molecular dynamics (MD) simulations to investigate directly spontaneous assembly of GAIV peptide and its structural dynamics, the formation of an induced water pore across the membrane and consequently ion conduction through the membrane.

Wei, Chenyu↗
NASA NTRSDate not supplied

Interfacial thermodynamics of cryogenic fluids: The effect on non-condensable gas on fluid storage

Propellant tanks that contain cryogenic fluids (CFs) in low gravity conditions are usually pressurized with non-condensable (NC) gases for fast extraction. Unfortunately, the presence of NC gases causes CFs to exhibit higher boil-offs compared to pure CF systems. For optimal utilization of the cryogenic fuels, these higher boil-off rates must be minimized. Our goal is to quantify the effects NC gases have on the evaporation and condensation dynamics of CFs to develop strategies that mitigate the higher boil-offs. The hypothesis is that the NC gas accumulates around the CF’s liquid vapor interface in the Knudsen layer, creating a kinetic barrier for both evaporation and condensation. Testing this hypothesis with experimental methods is difficult due to the transient nature of the Knudsen layer, which is only a few nanometers thick. Furthermore, experimental capabilities are limited in low gravity conditions and extremely expensive. Accordingly, we approach this problem from a theoretical perspective. In this study, we employ molecular dynamics (MD) simulations to probe the interfacial mechanisms that affect CF’s evaporation and condensation at varying concentrations of NC gases. Specifically, we use nitrogen (N2) and oxygen (O2) as our CFs and neon (Ne) as our NC gas. Using MD simulations, we show that Ne accumulates at N2 liquid vapor interface across a wide range of Ne concentrations, thereby impeding mass transport of N2. Our simulations allow for direct computation of the molar flux as well as the mass accommodation coefficient (MAC) which can then be used as input parameters to continuum fluid dynamics (CFD) models for optimal storage tank design.

Michael Robert DeLyser↗
NASA NTRSDate not supplied

Predicting Char Yield of High-Temperature Resins

A simulation protocol is developed to predict the char yield of organic resins during high-temperature processing. Such in silico methods can help screen promising new formulations for advanced materials, but previously no chemistry-sensitive technique existed to predict the important experimental value of char yield. The method utilizes a reactive force field (ReaxFF) to model the chemical transformation of precursor monomers into carbonized structures during three processing stages: ramp-up to processing temperatures (~3000 K), pyrolysis, and quenching. Achieving good agreement with experimental char yields requires continuous removal of small byproduct molecules to mimic outgassing, and the application of high pressure to eliminate porosity and encourage graphitization. More than ten different resin chemistries are investigated, including arylacetylenes, cyanate esters, phthalonitriles, and polyimides, representing a diverse group of precursors with respect to initial cyclic content, heteroatoms and reactive groups. The protocol correctly predicts the relative char yield between the investigated chemistries and provides quantitative agreement with experimental values, especially for high char yield resins. The properties of the resins during processing are compared, including outgassing products, morphology of the final chemical configurations, cyclic content and mechanical properties.

polymer simulations, molecular dynamics↗