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 289 records · Page 16

The Structural Properties of Vapor Deposited Water Ice and Astrophysical Implications

Films of vapor deposited water ice at low temperature (T<30 K) show a number of interesting structural changes during a gradual warmup. We would like to talk about the structure of the low temperature high density amorphous form of water ice, the process of crystallization, and some recent work on the morphological changes of water ice films at high temperature. The studies of the high density amorphous form are from in-situ electron microscopy as well as numerical simulations of molecular dynamics and have lead to new insights into the physical distinction between this high density amorphous form and the low density amorphous form. For the process of crystallization, we propose a model that describes the crystallization of water ice from the amorphous phase to cubic ice in terms of the nucleation of small domains in the ice. This model agrees well with the behavior of water ice in our electron microscopy studies and finds that pure water above the glass transition is a strong liquid. In more recent work, we have concentrated on temperatures above the crystallization temperature and we find interesting morphological changes related to the decrease in viscosity of the amorphous component in the cubic crystalline regime. Given enough time, we would like to put these results in an astrophysical context and discuss some observed features of the frost on interstellar grains and the bulk ice in comets.

Jenniskens, P.↗

Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes

Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.

Jaffe, Richard↗

Comparison of Oxygen Gauche Effects in Poly(Oxyethylene) and Poly(ethylene terephtylene) Based on Quantum Chemistry Calculations

The so-called oxygen gauche effect in poly(oxyethylene) (POE) and its model molecules such as 1,2-dimethoxyethane (DME) and diglyme (CH3OC2H4OC2H4OCH3) is manifested in the preference for gauche C-C bond conformations over trans. This has also been observed for poly(ethylene terephthalate) (PET). Our previous quantum chemistry calculations demonstrated that the large C-C gauche population in DME is due, in part, to a low-lying tg +/- g+ conformer that exhibits a substantial 1,5 CH ... O attraction. New calculations will be described that demonstrate the accuracy of the original quantum chemistry calculations. In addition, an extension of this work to model molecules for PET will be presented. It is seen that the C-C gauche preference is much stronger in 1,2 diacetoxyethane than in DME. In addition, there exist low-lying tg +/- g+/- and g+/-g+/-g+/- conformers that exhibit 1,5 CH ... O attractions involving the carbonyl oxygens. It is expected that the -O-C-C-O- torsional properties will be quite different in these two polymers. The quantum chemistry results are used to parameterize rotational isomeric states models (RIS) and force fields for molecular dynamics simulations of these polymers.

Jaffe, Richard↗

Anisotropic Plasticity of BN Nanotubes

Plastic collapse of compressed BN nanotubes are investigated and compared with carbon nanotubes of similar nature. Using a generalized tight-binding molecular dynamics (GTBMD) method for system containing B, N and C atoms we compute stiffness and plastic collapse of BN and C nanotubes under axial compression. For small compressional strain, BN nanotubes are found to be about 92% as stiff as similar C nanotubes. Due to BN bond buckling effect, however, the elastic limit of BN nanotubes is found to be more than C nanotubes. A route to plasticity is explored in which we find that at elastic limit the accumulated strain is released by a local plastic deformation of the nanotube. The mechanism of strain release and the resulting plastic deformation, however. are anisotropic in nature. The strain is released preferentially towards N as leading edge of a buckled BN bond and the tube, compressed at both ends, plastically collapses preferentially towards one end. Details of the anisotropic plasticity and prospective applications will be discussed in this presentation.

Madhu, Menon↗

Proton Transport

The transport of protons across membranes is an essential process for both bioenergetics of modern cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATP), synthesized from adenosine diphosphate. ATP, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this Study was: how the same process can be accomplished with the aid of similar but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC) which is a good model of the biological membranes forming cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M$_2$ protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M$_2$ protein is 97 amino acids in length, but a fragment 25 amino acids long. which contains a transmembrane domain of 19 amino acids flanked by three amino acids on each side. is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules. and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small ions, such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel with the histidine residues in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel loses its structural integrity and one helix moves away from the remaining three. Additional information is contained within the original extended abstract.

Pohorille, Andrew↗

Carbon Nanotube Based Molecular Electronics

Carbon nanotubes and the nanotube heterojunctions have recently emerged as excellent candidates for nanoscale molecular electronic device components. Experimental measurements on the conductivity, rectifying behavior and conductivity-chirality correlation have also been made. While quasi-one dimensional simple heterojunctions between nanotubes with different electronic behavior can be generated by introduction of a pair of heptagon-pentagon defects in an otherwise all hexagon graphene sheet. Other complex 3- and 4-point junctions may require other mechanisms. Structural stability as well as local electronic density of states of various nanotube junctions are investigated using a generalized tight-binding molecular dynamics (GDBMD) scheme that incorporates non-orthogonality of the orbitals. The junctions investigated include straight and small angle heterojunctions of various chiralities and diameters; as well as more complex 'T' and 'Y' junctions which do not always obey the usual pentagon-heptagon pair rule. The study of local density of states (LDOS) reveal many interesting features, most prominent among them being the defect-induced states in the gap. The proposed three and four pointjunctions are one of the smallest possible tunnel junctions made entirely of carbon atoms. Furthermore the electronic behavior of the nanotube based device components can be taylored by doping with group III-V elements such as B and N, and BN nanotubes as a wide band gap semiconductor has also been realized in experiments. Structural properties of heteroatomic nanotubes comprising C, B and N will be discussed.

Srivastava, Deepak↗

Thermal Expansion and Diffusion Coefficients of Carbon Nanotube-Polymer Composites

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

Wei, Chengyu↗

Fullerene Derived Molecular Electronic Devices

The carbon Nanotube junctions have recently emerged as excellent candidates for use as the building blocks in the formation of nanoscale electronic devices. While the simple joint of two dissimilar tubes can be generated by the introduction of a pair of heptagon-pentagon defects in an otherwise perfect hexagonal grapheme sheet, more complex joints require other mechanisms. In this work we explore structural and electronic properties of complex 3-point junctions of carbon nanotubes using a generalized tight-binding molecular-dynamics scheme.

Menon, Madhu↗

Carbon Nanotubes: Molecular Electronic Components

The carbon Nanotube junctions have recently emerged as excellent candidates for use as the building blocks in the formation of nanoscale molecular electronic networks. While the simple joint of two dissimilar tubes can be generated by the introduction of a pair of heptagon-pentagon defects in an otherwise perfect hexagonal graphene sheet, more complex joints require other mechanisms. In this work we explore structural characteristics of complex 3-point junctions of carbon nanotubes using a generalized tight-binding molecular-dynamics scheme. The study of pi-electron local densities of states (LDOS) of these junctions reveal many interesting features, most prominent among them being the defect-induced states in the gap.

Srivastava, Deepak↗

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↗

Observation and Modeling of Single Wall Carbon Nanotube Bend Junctions

Single wall carbon nanotube (SWNT) bends, with diameters from approx. 1.0 to 2.5 nm and bend angles from 18 deg. to 34 deg., are observed in catalytic decomposition of hydrocarbons at 600 - 1200 C. An algorithm using molecular dynamics simulation (MD) techniques is developed to model these structures that are considered to be SWNT junctions formed by topological defects (i.e. pentagon-heptagon pairs). The algorithm is used to predict the tube helicities and defect configurations for bend junctions using the observed tube diameters and bend angles. The number and arrangement of the defects at the junction interfaces are found to depend on the tube helicities and bend angle. The structural and energetic calculations using the Brenner potential show a number of stable junction configurations for each bend angle with the 34 deg. bends being more stable than the others. Tight binding calculations for local density of state (LDOS) and transmission coefficients are carried out to investigate electrical properties of the bend junctions.

Han, Jie↗

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

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

Dzegilenko, Fendor N.↗

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.↗

Research in Computational Astrobiology

We present results from several projects in the new field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. We have developed a procedure for calculating long-range effects in molecular dynamics using a plane wave expansion of the electrostatic potential. This method is expected to be highly efficient for simulating biological systems on massively parallel supercomputers. We have perform genomics analysis on a family of actin binding proteins. We have performed quantum mechanical calculations on carbon nanotubes and nucleic acids, which simulations will allow us to investigate possible sources of organic material on the early earth. Finally, we have developed a model of protobiological chemistry using neural networks.

Chaban, Galina↗

Collaborative Simulation Grid: Multiscale Quantum-Mechanical/Classical Atomistic Simulations on Distributed PC Clusters in the US and Japan

A multidisciplinary, collaborative simulation has been performed on a Grid of geographically distributed PC clusters. The multiscale simulation approach seamlessly combines i) atomistic simulation backed on the molecular dynamics (MD) method and ii) quantum mechanical (QM) calculation based on the density functional theory (DFT), so that accurate but less scalable computations are performed only where they are needed. The multiscale MD/QM simulation code has been Grid-enabled using i) a modular, additive hybridization scheme, ii) multiple QM clustering, and iii) computation/communication overlapping. The Gridified MD/QM simulation code has been used to study environmental effects of water molecules on fracture in silicon. A preliminary run of the code has achieved a parallel efficiency of 94% on 25 PCs distributed over 3 PC clusters in the US and Japan, and a larger test involving 154 processors on 5 distributed PC clusters is in progress.

Kikuchi, Hideaki↗

The M2 Proton Channel of Influenza Virus: How Does It Work?

The transport of protons across membranes is an essential process for both bioenergetics of modem cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATR), synthesized from adenosine diphosphate. ATR, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this study was: how the same process can be accomplished with the aid of similar, but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC), which is a good model of the biological membranes focusing cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M2 protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M2 protein is 97 amino acids in length, but a fragment 25 amino acids long, which contains a transmembrane domain of 19 amino acids flanked by 3 amino acids on each side, is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules, and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional, important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel, in which the histidine residues were in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel looses its structural integrity and one helix moves away from the remaining three. This result indicates that such a mechanism of proton transport is unlikely. The simulations revealed that translocation along a network of water molecules in the channel and tautomerization of the histidine residues in the M2 proteins in the most likely mechanism of proton transport. The results not only explain how a remarkably simple protein system can efficiently aid in the formation of proton gradients across cell walls, but also suggest how this system can be genetically re-engineered to become a directional, reversible proton pump. Such a pump can provide energy to laboratory-built models of simple cellular systems. If they were successfully constructed it would greatly advance our understanding of the beginnings of life and find important applications in medicine and pharmacology.

Pohorille, Andrew↗

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗

Tensile Yielding of Multi-Wall Carbon Nanotube

The tensile yielding of multiwall carbon nanotubes (MWCNTs) has been studied using Molecular Dynamics simulations and a Transition State Theory based model. We find a strong dependence of the yielding on the strain rate. A critical strain rate has been predicted above/below which yielding strain of a MWCNT is larger/smaller than that of the corresponding single-wall carbon nanotubes. At experimentally feasible strain rate of 1% /hour and T = 300K, the yield strain of a MWCNT is estimated to be about 3-4 % higher than that of an equivalent SWCNT (Single Wall Carbon Nanotube), in good agreement with recent experimental observations.

Wei, Chenyu↗