Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Coarse-Grained Molecular Dynamics Simulation of the Aggregation Properties of Multiheaded Cationic Surfactants in Water

The aggregation property of multiheaded surfactants has been investigated by constant pressure molecular dynamics (MD) simulation in aqueous medium. The model multiheaded surfactants contain more than one headgroup (x = 2, 3, and 4) for a single tail group. This increases the hydrophilic charge progressively over the hydrophobic tail which has dramatic consequences in the aggregation behavior. In particular, we have looked at the change in the aggregation property such as critical micellar concentration (cmc), aggregation number, and size of the micelles for the multiheaded surfactants in water. We find with increasing number of headgroups of the Multiheaded surfactants that the cmc values increase and the aggregation numbers as well as the size of the micelles decrease. These trends are in agreement with the experimental findings as reported earlier with x = 1, 2, and 3. We also predict the aggregation properties of multiheaded surfactant With four headgroups (x = 4) for which no experimental studies exist yet.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

Molecular Dynamics Simulation of the Bio-Adhesion in Molecular Motors

Molecular dynamics (MD) simulation is employed to study the bio-adhesion in F1 ATP molecular motor. Histidine-peptide is widely used as linkage in micro systems because of its strong binding strength to metals. This paper focuses on the adhesion between a synthetic peptide containing 6xHis-tag (Gly-Gly-Lys-Gly-Gly-Lys-Gly-Gly-His-His-His-His-His-His) and metal substrate, which is used to define the position of the F1 ATP molecular motor on the metal substrate. It is shown that the binding strength between histidine and nickel substrate is the strongest, while that of copper is smaller and that of gold substrate is the smallest. From the result of simulation, we find that the stability of adhesion between histidine and the metal substate result of the ringed structure in histidine.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Surface physics at the nano-scale via scanning probe microscopy and molecular dynamics simulations

Probe-based scanning microscopes, such as the STM and the AFM, are used to obtain the topographical and electronic structure maps of material surfaces, and to modify their morphologies on nanoscopic scales. They have generated new areas of research in condensed matter physics and materials science. We will review some examples from the fields of experimental nano-mechanics, nano-electronics and nano-magnetism. These now form the basis of the emerging field of Nano-technology. A parallel development has been brought about in the field of Computational Nano-science, using quantum-mechanical techniques and computer-based numerical modelling, such as the Molecular Dynamics (MD) simulation method. We will report on the simulation of nucleation and growth of nano-phase films on supporting substrates. Furthermore, a theoretical modelling of the formation of STM images of metallic clusters on metallic substrates will also be discussed within the non-equilibrium Keldysh Green function method to study the effects of coherent tunnelling through different atomic orbitals in a tip-sample geometry.

Nanoindentation of polysilicon and single crystal silicon: Molecular dynamics simulation and experimental validation

This paper presents novel advances on the deformation behaviour of polycrystalline and single crystal silicon using molecular dynamics (MD) simulation and validation of the same via nanoindentation experiments. In order to unravel the mechanism of deformation, four simulations were performed: Indentation of polycrystalline silicon substrate with a (i) Berkovich pyramidal and a (ii) spherical (arc) indenter, and indentation of a single crystal silicon substrate with these two indenters. The simulation results reveal that high pressure phase transformation (HPPT) in silicon (Si-I to Si-II phase transformation) occurred in all cases, however, its extent and the manner in which it occurred differed significantly between polycrystalline silicon and single crystal silicon, and was the main driver of differences in nanoindentation deformation behaviour between the two types of silicon. An interesting observation was that in polycrystalline silicon, the HPPT was observed to occur preferentially along the grain boundaries than across the grain boundaries. An automated dislocation extraction algorithm (DXA) revealed no dislocations in the deformation zone, suggesting HPPT to be the primary mechanism in inducing plasticity in silicon.

Diamond machining of silicon: A review of advances in molecular dynamics simulation

Molecular dynamics (MD) simulation has enhanced our understanding about ductile-regime machining of brittle materials such as silicon and germanium. In particular, MD simulation has helped understand the occurrence of brittle–ductile transition due to the high-pressure phase transformation (HPPT), which induces Herzfeld–Mott transition. In this paper, relevant MD simulation studies in conjunction with experimental studies are reviewed with a focus on (i) the importance of machining variables: undeformed chip thickness, feed rate, depth of cut, geometry of the cutting tool in influencing the state of the deviatoric stresses to cause HPPT in silicon, (ii) the influence of material properties: role of fracture toughness and hardness, crystal structure and anisotropy of the material, and (iii) phenomenological understanding of the wear of diamond cutting tools, which are all non-trivial for cost-effective manufacturing of silicon. The ongoing developmental work on potential energy functions is reviewed to identify opportunities for overcoming the current limitations of MD simulations. Potential research areas relating to how MD simulation might help improve existing manufacturing technologies are identified which may be of particular interest to early stage researchers.

Influence of temperature on the anisotropic cutting behaviour of single crystal silicon: A molecular dynamics simulation investigation

Using molecular dynamics (MD) simulation, this paper investigates anisotropic cutting behaviour of single crystal silicon in vacuum under a wide range of substrate temperatures (300 K, 500 K, 750 K, 850 K, 1173 K and 1500 K). Specific cutting energy, force ratio, stress in the cutting zone and cutting temperature were the indicators used to quantify the differences in the cutting behaviour of silicon. A key observation was that the specific cutting energy required to cut the (111) surface of silicon and the von Mises stress to yield the silicon reduces by 25% and 32%, respectively, at 1173 K compared to what is required at 300 K. The room temperature cutting anisotropy in the von Mises stress and the room temperature cutting anisotropy in the specific cutting energy (work done by the tool in removing unit volume of material) were obtained as 12% and 16% respectively. It was observed that this changes to 20% and 40%, respectively, when cutting was performed at 1500 K, signifying a very strong correlation between the anisotropy observed during cutting and the machining temperature. Furthermore, using the atomic strain criterion, the width of primary shear zone was found to vary with the orientation of workpiece surface and temperature i.e. it remains narrower while cutting the (111) surface of silicon or at higher machining temperatures. A major anecdote of the study based on the potential function employed in the study is that, irrespective of the cutting plane or the cutting temperature, the state of the cutting edge of the diamond tool did not show direct diamond to graphitic phase transformation.