MOLECULAR BEAMS AND SURFACE INTERACTIONS
Gas dynamics - molecular beams and surface interactions
SURFACE REACTION↗
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.
Free-Energy Calculations. A Mathematical Perspective
Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.
free energy↗
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↗
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↗
A molecular model for ice nucleation and growth, attachment 1
The quantum mechanical technique is used to study ionic, configurational, and impurity defects in the ice surface. In addition to static calculations of the energetics of the water monomer-ice surface interactions, molecular dynamics studies were initiated. The calculations of the monomer-ice surface interaction, molecular dynamics studies were initiated. The calculations of monomer-ice surface interactions indicate that many adsorption sites exist on the ice surfaces and that the barriers between bonding sites are relatively low. Bonding on the prism face of ice is preferentially above lattice sites.
Plummer, P. L. M.↗
Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential
The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.
molecular dynamics↗
Building and Breaking Carbon Composites with REACTER
Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.
Molecular Dynamics↗
Computer Simulation of the Forces Acting on a Submerged Polystyrene Probe as it Approaches the Succinonitrile Melt-Solid Interface
A Modeling approach to simulate both mesoscale and microscopic forces acting in a typical AFM experiment is presented. A mesoscale level interaction between the cantilever tip and the sample surface is primarily described by the balance of attractive Van der Waals and repulsive forces. Ultimately, the goal is to measure the forces between a particle and the crystal-melt interface. Two modes of AFM operation are considered in this paper - a stationary and a "tapping" one. The continuous mechanics approach to model tip-surface interaction is presented. At microscopic levels, tip contamination and details of tip-surface interaction are modeled using a molecular dynamics approach for the case of polystyrene - succinonitrile contact. Integration of the mesoscale model with a molecular dynamic model is discussed.
Bune, Andris V.↗
Interaction of a sodium ion with the water liquid-vapor interface
Molecular dynamics results are presented for the density profile of a sodium ion near the water liquid-vapor interface at 320 K. These results are compared with the predictions of a simple dielectric model for the interaction of a monovalent ion with this interface. The interfacial region described by the model profile is too narrow and the profile decreases too abruptly near the solution interface. Thus, the simple model does not provide a satisfactory description of the molecular dynamics results for ion positions within two molecular diameters from the solution interface where appreciable ion concentrations are observed. These results suggest that surfaces associated with dielectric models of ionic processes at aqueous solution interfaces should be located at least two molecular diameters inside the liquid phase. A free energy expense of about 2 kcal/mol is required to move the ion within two molecular layers of the free water liquid-vapor interface.
NASA Center ARC↗
Mechanism of unassisted ion transport across membrane bilayers
To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.
NASA Discipline Exobiology↗
Dynamics of Molecular Hydrogen in Hypersaline Microbial Mars
Early Earth microbial communities that centered around the anaerobic decomposition of organic molecular hydrogen as a carrier of electrons, regulator of energy metabolism, and facilitator of syntroph'c microbial interactions. The advent of oxygenic photosynthetic organisms added a highly dynamic and potentially dominant term to the hydrogen economy of these communities. We have examined the daily variations of hydrogen concentrations in cyanobacteria-dominated microbial mats from hypersaline ponds in Baja California Sur, Mexico. These mats bring together phototrophic and anaerobic bacteria (along with virtually all other trophic groups) in a spatially ordered and chemically dynamic matrix that provides a good analog for early Earth microbial ecosystems. Hydrogen concentrations in the photic zone of the mat can be three orders of magnitude or more higher than in the photic zone, which are, in turn, an order of magnitude higher than in the unconsolidated sediments underlying the mat community. Within the photic zone, hydrogen concentrations can fluctuate dramatically during the diel (24 hour day-night) cycle, ranging from less than 0.001% during the day to nearly 10% at night. The resultant nighttime flux of hydrogen from the mat to the environment was up to 17% of the daytime oxygen flux. The daily pattern observed is highly dependent on cyanobacterial species composition within the mat, with Lyngbya-dominated systems having a much greater dynamic range than those dominated by Microcoleus; this may relate largely to differing degrees of nitrogen-fixing and fermentative activity in the two mats. The greatest H2 concentrations and fluxes were observed in the absence of oxygen, suggesting an important potential feedback control in the context of the evolution of atmospheric composition. The impact of adding this highly dynamic photosynthetic term to the hydrogen economy of early microbial ecosystems must have been substantial. From an evolutionary standpoint, the H2 generated in mats could have represented a very important new source of electrons and energy - but one that could not be harnessed without substantial adaptation to the highly variable chemistry of the mat surface. In addition, the emergent chemistry of anaerobic communities is often highly dependent on ambient hydrogen concentrations, so that incorporation of these communities into photosynthetic mats could have significantly affected the composition and flux of reduced "biosignature' gases to the environment.
Hoehler, Tori M.↗
Calculation of Non-Bonded Forces Due to Sliding of Bundled Carbon Nanotubes
An important consideration for load transfer in bundles of single-walled carbon nanotubes is the nonbonded (van der Waals) forces between the nanotubes and their effect on axial sliding of the nanotubes relative to each other. In this research, the non-bonded forces in a bundle of seven hexagonally packed (10,10) single-walled carbon nanotubes are represented as an axial force applied to the central nanotube. A simple model, based on momentum balance, is developed to describe the velocity response of the central nanotube to the applied force. The model is verified by comparing its velocity predictions with molecular dynamics simulations that were performed on the bundle with different force histories applied to the central nanotube. The model was found to quantitatively predict the nanotube velocities obtained from the molecular dynamics simulations. Both the model and the simulations predict a threshold force at which the nanotube releases from the bundle. This force converts to a shear yield strength of 10.5-11.0 MPa for (10,10) nanotubes in a bundle.
Frankland, S. J. V.↗
Computer simulation of surface and film processes
All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.
Tiller, W. A.↗
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.↗
Computational Study on Full-length Human Ku70 with Double Stranded DNA: Dynamics, Interactions and Functional Implications
The Ku70/80 heterodimer is the first repair protein in the initial binding of double-strand break (DSB) ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. In this study we constructed a full-length human Ku70 structure based on its crystal structure, and performed 20 ns conventional molecular dynamic (CMD) simulations on this protein and several other complexes with short DNA duplexes of different sequences. The trajectories of these simulations indicated that, without the topological support of Ku80, the residues in the bridge and C-terminal arm of Ku70 are more flexible than other experimentally identified domains. We studied the two missing loops in the crystal structure and predicted that they are also very flexible. Simulations revealed that they make an important contribution to the Ku70 interaction with DNA. Dislocation of the previously studied SAP domain was observed in several systems, implying its role in DNA binding. Targeted molecular dynamic (TMD) simulation was also performed for one system with a far-away 14bp DNA duplex. The TMD trajectory and energetic analysis disclosed detailed interactions of the DNA-binding residues during the DNA dislocation, and revealed a possible conformational transition for a DSB end when encountering Ku70 in solution. Compared to experimentally based analysis, this study identified more detailed interactions between DNA and Ku70. Free energy analysis indicated Ku70 alone is able to bind DNA with relatively high affinity, with consistent contributions from various domains of Ku70 in different systems. The functional implications of these domains in the processes of Ku heterodimerization and DNA damage recognition and repair can be characterized in detail based upon this analysis.
Hu, Shaowen↗
Tensile Strength of Carbon Nanotubes Under Realistic Temperature and Strain Rate
Strain rate and temperature dependence of the tensile strength of single-wall carbon nanotubes has been investigated with molecular dynamics simulations. The tensile failure or yield strain is found to be strongly dependent on the temperature and strain rate. A transition state theory based predictive model is developed for the tensile failure of nanotubes. Based on the parameters fitted from high-strain rate and temperature dependent molecular dynamics simulations, the model predicts that a defect free micrometer long single-wall nanotube at 300 K, stretched with a strain rate of 1%/hour, fails at about 9 plus or minus 1% tensile strain. This is in good agreement with recent experimental findings.
Wei, Chen-Yu↗
Plasticity and Kinky Chemistry of Carbon Nanotubes
Since their discovery in 1991, carbon nanotubes have been the subject of intense research interest based on early predictions of their unique mechanical, electronic, and chemical properties. Materials with the predicted unique properties of carbon nanotubes are of great interest for use in future generations of aerospace vehicles. For their structural properties, carbon nanotubes could be used as reinforcing fibers in ultralight multifunctional composites. For their electronic properties, carbon nanotubes offer the potential of very high-speed, low-power computing elements, high-density data storage, and unique sensors. In a continuing effort to model and predict the properties of carbon nanotubes, Ames accomplished three significant results during FY99. First, accurate values of the nanomechanics and plasticity of carbon nanotubes based on quantum molecular dynamics simulations were computed. Second, the concept of mechanical deformation catalyzed-kinky-chemistry as a means to control local chemistry of nanotubes was discovered. Third, the ease of nano-indentation of silicon surfaces with carbon nanotubes was established. The elastic response and plastic failure mechanisms of single-wall nanotubes were investigated by means of quantum molecular dynamics simulations.
Srivastava, Deepak↗
Multiscale Modeling for the Analysis for Grain-Scale Fracture Within Aluminum Microstructures
Multiscale modeling methods for the analysis of metallic microstructures are discussed. Both molecular dynamics and the finite element method are used to analyze crack propagation and stress distribution in a nanoscale aluminum bicrystal model subjected to hydrostatic loading. Quantitative similarity is observed between the results from the two very different analysis methods. A bilinear traction-displacement relationship that may be embedded into cohesive zone finite elements is extracted from the nanoscale molecular dynamics results.
Glaessgen, Edward H.↗