Molecular dynamics simulation (MDS) to study nanoscale machining processes

Molecular Dynamics Simulations (MDS) are constantly being used to make important contributions to our fundamental understanding of material behaviour, at the atomic scale, for a variety of thermodynamic processes. This chapter shows that molecular dynamics simulation is a robust numerical analysis tool in addressing a range of complex nanofinishing (machining) problems that are otherwise difficult or impossible to understand using other methods. For example the mechanism of nanometric cutting of silicon carbide is influenced by a number of variables such as machine tool performance, machining conditions, material properties, and cutting tool performance (material microstructure and physical geometry of the contact) and all these variables cannot be monitored online through experimental examination. However, these could suitably be studied using an advanced simulation based approach such as MDS. This chapter details how MD simulation can be used as a research and commercial tool to understand key issues of ultra precision manufacturing research problems and a specific case was addressed by studying diamond machining of silicon carbide. While this is appreciable, there are a lot of challenges and opportunities in this fertile area. For example, the world of MD simulations is dependent on present day computers and the accuracy and reliability of potential energy functions [109]. This presents a limitation: Real-world scale simulation models are yet to be developed. The simulated length and timescales are far shorter than the experimental ones which couples further with the fact that contact loading simulations are typically done in the speed range of a few hundreds of m/sec against the experimental speed of typically about 1 m/sec [17]. Consequently, MD simulations suffer from the spurious effects of high cutting speeds and the accuracy of the simulation results has yet to be fully explored. The development of user-friendly software could help facilitate molecular dynamics as an integral part of computer-aided design and manufacturing to tackle a range of machining problems from all perspectives, including materials science (phase of the material formed due to the sub-surface deformation layer), electronics and optics (properties of the finished machined surface due to the metallurgical transformation in comparison to the bulk material), and mechanical engineering (extent of residual stresses in the machined component) [110]. Overall, this chapter provided key information concerning diamond machining of SiC which is classed as hard, brittle material. From the analysis presented in the earlier sections, MD simulation has helped in understanding the effects of crystal anisotropy in nanometric cutting of 3C-SiC by revealing the atomic-level deformation mechanisms for different crystal orientations and cutting directions. In addition to this, the MD simulation revealed that the material removal mechanism on the (111) surface of 3C-SiC (akin to diamond) is dominated by cleavage. These understandings led to the development of a new approach named the “surface defect machining” method which has the potential to be more effective to implement than ductile mode micro laser assisted machining or conventional nanometric cutting.

Molecular dynamics simulation investigation on the plastic flow behaviour of silicon during nanometric cutting

Molecular dynamics (MD) simulation was carried out to acquire an in-depth understanding of the flow behaviour of single crystal silicon during nanometric cutting on three principal crystallographic planes and at different cutting temperatures. The key findings were that (i) the substrate material underneath the cutting tool was observed for the first time to experience a rotational flow akin to fluids at all the tested temperatures up to 1200 K. (ii) The degree of flow in terms of vorticity was found higher on the (1 1 1) crystal plane signifying better machinability on this orientation in accord with the current pool of knowledge (iii) an increase in the machining temperature reduces the springback effect and thereby the elastic recovery and (iv) the cutting orientation and the cutting temperature showed significant dependence on the location of the stagnation region in the cutting zone of the substrate.

Interactions between pyrazole derived enantiomers and Chiralcel OJ: Prediction of enantiomer absolute configurations and elution order by molecular dynamics simulations

The separation of enantiomers and confirmation of their absolute configurations is significant in the development of chiral drugs. The interactions between the enantiomers of chiral pyrazole derivative and polysaccharide-based chiral stationary phase cellulose tris(4-methylbenzoate) (Chiralcel OJ) in seven solvents and under different temperature were studied using molecular dynamics simulations. The results show that solvent effect has remarkable influence on the interactions. Structure analysis discloses that the different interactions between two isomers and chiral stationary phase are dependent on the nature of solvents, which may invert the elution order. The computational method in the present study can be used to predict the elution order and the absolute configurations of enantiomers in HPLC separations and therefore would be valuable in development of chiral drugs.

Self-healing of cracks during ductile regime machining of silicon: Insights from molecular dynamics simulation

During nanoindentation and ductile-regime machining of silicon, a phenomenon known as “self-healing” takes place in that the microcracks, microfractures, and small spallings generated during the machining are filled by the plastically flowing ductile phase of silicon. However, this phenomenon has not been observed in simulation studies. In this work, using a long-range potential function, molecular dynamics simulation was used to provide an improved explanation of this mechanism. A unique phenomenon of brittle cracking was discovered, typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves being filled by plastically flowing silicon during cutting. This observation is supported by the direct imaging. The simulated X-ray diffraction analysis proves that in contrast to experiments, Si-I to Si-II (beta tin) transformation during ductile-regime cutting is highly unlikely and solid-state amorphisation of silicon caused solely by the machining stress rather than the cutting temperature is the key to its brittle-ductile transition observed during the MD simulations

A disappearing act performed by magnesium: the nucleotide exchange mechanism of Ran GTPase by Quantum mechanics/Molecular mechanics studies

The Ran GTPase protein is a guanine nucleotide-binding protein (GNBP) with an acknowledged profile in cancer onset, progression and metastases. The complex mechanism adopted by GNBPs in exchanging GDP for GTP is an intriguing process and crucial for Ran viability. The successful completion of the process is a fundamental aspect of propagating downstream signalling events. QM/MM molecular dynamics simulations were employed in this study to provide a deeper mechanistic understanding of the initiation of nucleotide exchange in Ran. Results indicate significant disruption of the metal-binding site upon interaction with RCC1 (the Ran guanine nucleotide exchange factor), overall culminating in the prominent shift of the divalent magnesium ion. The observed ion drifting is reasoned to occur as a consequence of the complex formation between Ran and RCC1 and is postulated to be a critical factor in the exchange process adopted by Ran. This is the first report to observe and detail such intricate dynamics for a protein in Ras superfamily.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Structural basis of substrate selectivity in the glycerol-3-phosphate: phosphate antiporter GlpT.

Major facilitators represent the largest superfamily of secondary active transporter proteins and catalyze the transport of an enormous variety of small solute molecules across biological membranes. However, individual superfamily members, although they may be architecturally similar, exhibit strict specificity toward the substrates they transport. The structural basis of this specificity is poorly understood. A member of the major facilitator superfamily is the glycerol-3-phosphate (G3P) transporter (GlpT) from the Escherichia coli inner membrane. GlpT is an antiporter that transports G3P into the cell in exchange for inorganic phosphate (Pi). By combining large-scale molecular-dynamics simulations, mutagenesis, substrate-binding affinity, and transport activity assays on GlpT, we were able to identify key amino acid residues that confer substrate specificity upon this protein. Our studies suggest that only a few amino acid residues that line the transporter lumen act as specificity determinants. Whereas R45, K80, H165, and, to a lesser extent Y38, Y42, and Y76 contribute to recognition of both free Pi and the phosphate moiety of G3P, the residues N162, Y266, and Y393 function in recognition of only the glycerol moiety of G3P. It is the latter interactions that give the transporter a higher affinity to G3P over Pi.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (Pi) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.

Correlated electron-ion dynamics in metallic systems

In this paper we briefly discuss the problem of simulating non-adiabatic processes in systems that are usefully modelled using molecular dynamics. In particular we address the problems associated with metals, and describe two methods that can be applied: the Ehrenfest approximation and correlated electron-ion dynamics (CEID). The Ehrenfest approximation is used to successfully describe the friction force experienced by an energetic particle passing through a crystal, but is unable to describe the heating of a wire by an electric current. CEID restores the proper heating.

Liquid Structure and Dynamics of Aqueous Isopropanol over gamma-Alumina

The liquid structures of thin films of aqueous solutions of 0, 7, 19, 50, and 100 mol % isopropanol above O/Al-terminated gamma-alumina surfaces have been investigated by means of classical molecular dynamics simulations. The structuring effect of the oxide oil the liquid mixtures is strong and heavily dependent on the local structure of the oxide. Two distinct re-ions are found oil the oxide Surface characterized by the degree of coordination of Al atoms. Above octahedral Al atoms, water and isopropanol molecules adsorb via the oxygen atoms to maximize the electrostatic interaction, whereas above tetrahedral Al sites the solvent molecules adsorb via hydrogen atoms with the oxygen atoms away front the surface. More mobility is found in the second layer compared with the first; however, its structure is still influenced significantly by the orientation of molecules in the first adsorbed layer. Qualitatively, the displacement of water from the surface by the adsorption of isopropanol occurs with 2.6 Water molecules lost for every alcohol molecule present based on the effective surface areas of the two species calculated from the pure simulations.

A molecular model for drug binding to tandem repeats of telomeric G-quadruplexes.

The extreme 3'-ends of human telomeres consist of 150–250 nucleotides of single-stranded DNA sequence together with associated proteins. Small-molecule ligands can compete with these proteins and induce a conformational change in the DNA to a four-stranded quadruplex arrangement, which is also no longer a substrate for the telomerase enzyme. The modified telomere ends provide signals to the DNA-damage-response system and trigger senescence and apoptosis. Experimental structural data are available on such quadruplex complexes comprising up to four telomeric DNA repeats, but not on longer systems that are more directly relevant to the single-stranded overhang in human cells. The present paper reports on a molecular modelling study that uses Molecular Dynamics simulation methods to build dimer and tetramer quadruplex repeats. These incorporate ligand-binding sites and are models for overhang–ligand complexes.

Comparison of Dynamics of Extracellular Accesses to the β(1) and β(2) Adrenoceptors Binding Sites Uncovers the Potential of Kinetic Basis of Antagonist Selectivity.

From the molecular mechanism of antagonist unbinding in the ß(1) and ß(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of ß(1) -selective Esmolol and ß(2) -selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the ß(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the ß(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the ß(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the ß(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface. © 2012 John Wiley & Sons A/S.