Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation

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.

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.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Ab initio Molecular Dynamics Simulations of the Structural Response of Solids to Ultrashort Laser and XUV Pulses

The theoretical model and underlying physics described in this thesis are about the interaction of femtosecond-laser and XUV pulses with solids. The key to understand the basics of such interaction is to study the structural response of the materials after laser interaction. Depending on the laser characteristics, laser-solid interaction can result in a wide range of structural responses such as solid-solid phase transitions, vacuum phonon squeezing, ultrafast melting, generation of coherent phonons, etc. During my research work, I have modeled the systems irradiated by low-, medium- and high-laser intensities, and studied different types of structural dynamics of solids at various laser fluences.

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Grid warehousing of molecular dynamics protein unfolding data

With the increasing awareness of protein folding disorders, the explosion of genomic information, and the need for efficient ways to predict protein structure, protein folding and unfolding has become a central issue in molecular sciences research. Molecular dynamics computer simulations are increasingly employed to understand the folding and unfolding of proteins. Running protein unfolding simulations is computationally expensive and finding ways to enhance performance is a grid issue on its own. However, more and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. This paper describes efforts to provide a grid-enabled data warehouse for protein unfolding data. We outline the challenge and present first results in the design and implementation of the data warehouse.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Assembling a xylanase-lichenase chimera through all-atom molecular dynamics simulations

Multifunctional enzyme engineering can improve enzyme cocktails for emerging biofuel technology. Molecular dynamics through structure-based models (SB) is an effective tool for assessing the tridimensional arrangement of chimeric enzymes as well as for inferring the functional practicability before experimental validation. This study describes the computational design of a bifunctional xylanase-lichenase chimera (XylLich) using the xynA and bglS genes from Bacillus subtilis. In silico analysis of the average solvent accessible surface area (SAS) and the root mean square fluctuation (RMSF) predicted a fully functional chimera, with minor fluctuations and variations along the polypeptide chains. Afterwards, the chimeric enzyme was built by fusing the xynA and bglS genes. XylLich was evaluated through small-angle X-ray scattering (SAXS) experiments, resulting in scattering curves with a very accurate fit to the theoretical protein model. The chimera preserved the biochemical characteristics of the parental enzymes, with the exception of a slight variation in the temperature of operation and the catalytic efficiency (k cat/Km). The absence of substantial shifts in the catalytic mode of operation was also verified. Furthermore, the production of chimeric enzymes could be more profitable than producing a single enzyme separately, based on comparing the recombinant protein production yield and the hydrolytic activity achieved for XylLich with that of the parental enzymes. © 2013 Elsevier B.V. All rights reserved.

Influence of Molecular Dynamics on the Dielectric Properties of Poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) -Based Devices

This paper uses Nuclear Magnetic Resonance (NMR) and Differential Scanning Calorimetry (DSC) techniques to study the molecular relaxations and phase transitions in poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) (F8BT), which has been extensively studied as the active thin film in organic devices. Besides the identification of the glass transition, beta relaxation and crystal-to-crystal phase transition, we correlate such phenomena with dielectric and transport mechanisms in diodes with F8BT as the active layer. The beta relaxation has been assigned to a transition at about 210 K measured by H-1 and C-13 solid state NMR, and can be attributed to local motions in the side chains. The glass transition has been detected by DSC and H-1 NMR. Dielectric spectroscopy (DS) carried out at low frequencies on diodes made from F8BT show two peaks which are coincident with the above transitions. This allowed us to correlate the electrical changes in the film with the onset of specific molecular motions. In addition, DS indicates a third peak related with a crystal-to-crystal phase transition. Finally, these transitions were correlated with changes in the carrier mobility recorded in thin films and published recently.

Molecular dynamics of poly(Ethylene Glycol) intercalated in clay, studied using 13C solid-state NMR

In this study, Cross-Polarization Magic-angle Spinning CP/MAS, 2D Exchange, Centerband-Only Detection of Exchange (CODEX), and Separated-Local-Field (SLF) NMR experiments were used to study the molecular dynamics of poly(ethylene glycol) (PEG) inside Hectorite/PEG intercalation compounds in both single- and double-layer configurations. The results revealed that the overall amplitude of the motions of the PEG chain in the single-layer configuration is considerably smaller than that observed for the double-layer intercalation compound. This result indicates that the effect of having the polymer chain interacting with both clay platelets is to produce a substantial decrease in the motional amplitudes of those chains. The presence of these dynamically restricted segments might be explained by the presence of anchoring points between the clay platelets and the PEG oxygen atoms, which was induced by the Na+ cations. By comparing the PEG motional amplitudes of the double-layered nanocomposites composed of polymers with different molecular weights, a decrease in the motional amplitude for the smaller PEG chain was observed, which might also be understood using the presence of anchoring points.

Molecular dynamics simulations of sheared brush-like systems

Mit Hilfe von Molekulardynamik-Simulationen untersuchen wir bürstenartige Systeme unter guten Lösungsmittelbedingungen. Diese Systeme sind, dank ihren vielfältigen Beschaffenheiten, die von Molekularparametern und äußeren Bedingungen abhängig sind, wichtig für viele industrielle Anwendungen. Man vermutet, dass die Polymerbürsten eine entscheidende Rolle in der Natur wegen ihrer einzigartigen Gleiteigenschaften spielen. Ein vergröbertes Modell wird verwendet, um die strukturellen und dynamischen Eigenschaften zweier hochkomprimierter Polymerbürsten, die eine niedrige Reibung aufweisen, zu untersuchen. Allerdings sind die Lubrikationseigenschaften dieser Systeme, die in vielen biologischen Systemen vorhanden sind, beeinflußt. Wir untersuchen so-genannte "weiche Kolloide", die zwischen den beiden Polymerbürsten eingebettet sind, und wie diese Makroobjekte auf die Polymerbürsten wirken.rnrnNicht-Gleichgewichts-Molekulardynamik-Simulationen werden durchgeführt, in denen die hydrodynamischen Wechselwirkungen durch die Anwendung des DPD-Thermostaten mit expliziten Lösungsmittelmolekülen berücksichtigt werden. Wir zeigen, dass die Kenntnis der Gleichgewichtseigenschaften des Systems erlaubt, dynamische Nichtgleichgewichtsigenschaften der Doppelschicht vorherzusagen.rnrnWir untersuchen, wie die effektive Wechselwirkung zwischen kolloidalen Einschlüßen durch die Anwesenheit der Bürsten (in Abhängigkeit der Weichheit der Kolloide und der Pfropfdichte der Bürsten) beeinflußt wird. Als nächsten Schritt untersuchen wir die rheologische Antwort von solchen komplexen Doppelschichten auf Scherung. Wir entwickeln eine Skalen-Theorie, die die Abhängigkeit der makroskopischen Transporteigenschaften und der lateralen Ausdehnung der verankerten Ketten von der Weissenberg Zahl oberhalb des Bereichs, in dem die lineare Antwort-Theorie gilt, voraussagt. Die Vorhersagen der Theorie stimmen gut mit unseren und früheren numerischen Ergebnissen und neuen Experimenten überein. Unsere Theorie bietet die Möglichkeit, die Relaxationszeit der Doppelschicht zu berechnen. Wenn diese Zeit mit einer charakteristischen Längenskala kombiniert wird, kann auch das ''transiente'' (nicht-stationäre) Verhalten beschrieben werden.rnrnrnWir untersuchen die Antwort des Drucktensors und die Deformation der Bürsten während der Scherinvertierung für grosse Weissenberg Zahlen. Wir entwickeln eine Vorhersage für die charakteristische Zeit, nach der das System wieder den stationären Zustand erreicht.rnrnrnElektrostatik spielt eine bedeutende Rolle in vielen biologischen Prozessen. Die Lubrikationseigenschaften der Polymerbürsten werden durch die Anwesenheit langreichweitiger Wechselwirkungen stark beeinflusst. Für unterschiedliche Stärken der elektrostatischen Wechselwirkungen untersuchen wir rheologische Eigenschaften der Doppelschicht und vergleichen mit neutralen Systemen. Wir studieren den kontinuierlichen Übergang der Systemeigenschaften von neutralen zu stark geladenen Bürsten durch Variation der Bjerrumlänge und der Ladungsdichte.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Molecular Dynamics Simulations of the d(TAT) Triple Helix

Molecular dynamics (MD) simulations have been used to study the dynamical and time-averaged characteristics of the DNA triple helix d(T)10âd(A)10âd(T)10. The structures sampled during the trajectory resemble closely the B-type model for the DNA triplex proposed on the basis of NMR data, although there are some subtle differences. Alternative P- and A-type conformations for the triplex, suggested from X-ray experiments, are not predicted to contribute significantly to the structure of the DNA triplex in solution. Comparison with the best available experimental data supports the correctnes of the MD-generated structures. The analysis of the collected data gives a detailed picture of the characteristics of triple-helix DNA. A new and interesting pattern of hydration, specific for triplex DNA, is an important observation. The results suggest that molecular dynamics can be useful for the study of novel nucleic acid structures.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.