Variation of diffusivity with the cation radii in molten salts of superionic conductors containing iodine anion: A molecular dynamics study

A molecular dynamics study of the dependence of diffusivity of the cation on ionic radii in molten AgI is reported. We have employed modified Parinello-Rahman-Vashistha interionic pair potential proposed by Shimojo and Kobayashi.(1) Our results suggest that the diffusivity of the cation exhibits an increase followed by a decrease as the ionic radius is increased. Several structural and dynamical properties are reported.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

Force Biased Molecular Dynamics Simulation Study of Effect of Dendrimer Generation on Interaction with DNA

We have studied the effect of dendrimer generation on the interaction between dsDNA and the PAMAM dendrimer using force biased simulation of dsDNA with three generations of dendrimer: G3, G4, and G5. Our results for the potential of mean force (PMF) and the dendrimer asphericity along the binding pathway, combined with visualization of the simulations, demonstrate that dendrimer generation has a pronounced impact on the interaction. The PMF increases linearly with increasing generation of the dendrimer. While, in agreement with previous results, we see an increase in the extent to which the dendrimer bends the dsDNA with increasing dendrimer generation, we also see that the deformation of the dendrimer is greater with smaller generation of the dendrimer. The larger dendrimer forces the dsDNA to conform to its structure, while the smaller dendrimer is forced to conform to the structure of the dsDNA. Monitoring the number of bound cations at different values of force bias distance shows the expected effect of ions being expelled when the dendrimer binds dsDNA.

Diffusion of pentane isomers in faujasite-type zeolites : NMR and molecular dynamics study

Diffusion of pentane isomers in zeolites NaX has been investigated using pulsed field gradient nuclear magnetic resonance (PFG-NMR) and molecular dynamics (MD) techniques respectively. Temperature and concentration dependence of diffusivities have been studied. The diffusivities obtained from NMR are roughly an order of magnitude smaller than those obtained from MD. The dependence of diffusivity on loading at high temperatures exhibits a type I behavior according to the classification of Karger and Pfeifer 1]. NMR diffusivities of the isomers exhibit the order D(n-pentane) > D(isopentane) > D(neopentane). The results from MD suggest that the diffusivities of the isomers follow the order D(n-pentane) < D(isopentane) < D(neopentane). The activation energies from NMR show E-a(n-pentane) < E-a(isopentane) < E-a(neopentane) whereas those from MD suggest the order E-a(n-pentane) > (isopentane) > E-a(neopentane). The latter follows the predictions of levitation effect whereas those of NMR appears to be due to the presence of defects in the zeolite crystals. The differences between diffusivities estimated by NMR and MD are attributed to the longer time and length scales sampled by the NMR technique, as compared to MD. (C) 2012 Elsevier Inc. All rights reserved.

Structure, molecular dynamics, and stress in a linear polymer

We present a study correlating uniaxial stress in a polymer with its underlying structure when it is strained. The uniaxial stress is significantly influenced by the mean-square bond length and mean bond angle. In contrast, the size and shape of the polymer, typically represented by the end-to-end length, mass ratio, and radius of gyration, contribute negligibly. Among externally set control variables, density and polymer chain length play a critical role in influencing the anisotropic uniaxial stress. Short chain polymers more or less behave like rigid molecules. Temperature and rate of loading, in the range considered, have a very mild effect on the uniaxial stress.

Molecular dynamics simulations of PPI dendrimer-drug complexes

Dendrimeric nanoparticles are potential drug delivery devices which can enhance the solubility of hydrophobic drugs, thus increasing their bioavailability and sustained release action. A quantitative understanding of the dendrimer-drug interactions can give valuable insight into the solubility and release profile of hydrophobic drug molecules in various solvent conditions. Fully atomistic molecular dynamics (MD) simulations have been performed to study the interactions of G5 PPIEDA (G5 ethylenediamine cored poly(propylene imine)) dendrimer and two well known drugs (Famotidine and Indomethacin) at different pH conditions. The study suggested that at low pH the dendrimer-drug complexes are thermodynamically unstable as compared to neutral and high pH conditions. Calculated Potential of Mean Force (PMF) by umbrella sampling showed that the release of drugs from the dendrimer at low pH is spontaneous, median release at neutral pH and slow release at high pH. In addition, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding free energy calculations were also performed at each umbrella sampling window to identify the various energy contributions. To understand the effect of dendrimer chemistry and topology on the solubility and release profile of drugs, this study is extended to explore the solubility and release profile of phenylbutazone drug complexed with G3 poly(amidoamine) and G4 diaminobutane cored PPI dendrimers. The results indicate that the pH-induced conformational changes in dendrimer, ionization states, dendrimer type and pK(a) of the guest molecules influence the free energy barrier and stability of complexation, and thus regulate drug loading, solubility and release.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Performance metrics in a hybrid MPI-OpenMP based molecular dynamics simulation with short-range interactions

We discuss the computational bottlenecks in molecular dynamics (MD) and describe the challenges in parallelizing the computation-intensive tasks. We present a hybrid algorithm using MPI (Message Passing Interface) with OpenMP threads for parallelizing a generalized MD computation scheme for systems with short range interatomic interactions. The algorithm is discussed in the context of nano-indentation of Chromium films with carbon indenters using the Embedded Atom Method potential for Cr-Cr interaction and the Morse potential for Cr-C interactions. We study the performance of our algorithm for a range of MPI-thread combinations and find the performance to depend strongly on the computational task and load sharing in the multi-core processor. The algorithm scaled poorly with MPI and our hybrid schemes were observed to outperform the pure message passing scheme, despite utilizing the same number of processors or cores in the cluster. Speed-up achieved by our algorithm compared favorably with that achieved by standard MD packages. (C) 2013 Elsevier Inc. All rights reserved.

Topological blocking in quantum quench dynamics

We study the nonequilibrium dynamics of quenching through a quantum critical point in topological systems, focusing on one of their defining features: ground-state degeneracies and associated topological sectors. We present the notion of ``topological blocking,'' experienced by the dynamics due to a mismatch in degeneracies between two phases, and we argue that the dynamic evolution of the quench depends strongly on the topological sector being probed. We demonstrate this interplay between quench and topology in models stemming from two extensively studied systems, the transverse Ising chain and the Kitaev honeycomb model. Through nonlocal maps of each of these systems, we effectively study spinless fermionic p-wave paired topological superconductors. Confining the systems to ring and toroidal geometries, respectively, enables us to cleanly address degeneracies, subtle issues of fermion occupation and parity, and mismatches between topological sectors. We show that various features of the quench, which are related to Kibble-Zurek physics, are sensitive to the topological sector being probed, in particular, the overlap between the time-evolved initial ground state and an appropriate low-energy state of the final Hamiltonian. While most of our study is confined to translationally invariant systems, where momentum is a convenient quantum number, we briefly consider the effect of disorder and illustrate how this can influence the quench in a qualitatively different way depending on the topological sector considered.

Conducting polymer-carbon black nanocomposite sensor for volatile organic compounds and correlating sensor response by molecular dynamics

Volatile organic compounds (VOCs) are present in our every day used products such as plastics, cosmetics, air fresheners, paint, etc. The determination of amount of VOC present in atmosphere can be carried out via various sensors. In this work a nanocomposite of a novel thiophene based conducting polymer and carbon black is used as a volatile organic compound sensor. The fabricated 2 lead chemiresistor sensor was tested for vapours of toluene, acetone, cylcohexane, and carbon tetrachloride. The sensor responds to all the vapours, however, exhibit maximum response to toluene vapours. The sensor was evaluated for various concentrations of toluene. The lower limit of detection of the sensor is 15 +/- 10 ppm. The study of the effect of humidity on senor response to toluene showed that the response decreases at higher humidity conditions. The surface morphology of the nanocomposite was characterized by scanning electron microscopy. Diffuse reflectance spectroscopy was used to investigate the absorption of vapours by the nanocomposite film. Contact angle measurements were used to present the effect of water vapour on the toluene response of nanocomposite film. Solubility parameter of the conducting polymer is predicted by molecular dynamics. The sensing behaviour of the conducting polymer is correlated with solubility parameter of the polymer. Dispersion interaction of conducting polymer with toluene is believed to be the reason for the selective response towards toluene. (C) 2014 Elsevier B.V. All rights reserved.

Molecular Dynamics Study of Phonon Screening in Graphene

Phonon interaction with electrons or phonons or with structural defects result in a phonon mode conversion. The mode conversion is governed by the frequency wave-vector dispersion relation. The control over phonon mode or the screening of phonon in graphene is studied using the propagation of amplitude modulated phonon wave-packet. Control over phonon properties like frequency and velocity opens up several wave guiding, energy transport and thermo-electric applications of graphene. One way to achieve this control is with the introduction of nano-structured scattering in the phonon path. Atomistic model of thermal energy transport is developed which is applicable to devices consisting of source, channel and drain parts. Longitudinal acoustic phonon mode is excited from one end of the device. Molecular dynamics based time integration is adopted for the propagation of excited phonon to the other end of the device. The amount of energy transfer is estimated from the relative change of kinetic energy. Increase in the phonon frequency decreases the kinetic energy transmission linearly in the frequency band of interest. Further reduction in transmission is observed with the tuning of channel height of the device by increasing the boundary scattering. Phonon mode selective transmission control have potential application in thermal insulation or thermo-electric application or photo-thermal amplification.

A molecular dynamics study of ambient and high pressure phases of silica: Structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, alpha-PbO2-type, and pyrite-type for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase. (C) 2014 AIP Publishing LLC.

Kinetics of phase separation in polymer mixtures: A molecular dynamics study

We present detailed results from a molecular dynamics (MD) simulation of phase-separation kinetics in polymer mixtures. Our MD simulations naturally incorporate hydrodynamic effects. We find that polymeric phase separation (with dynamically symmetric components) is in the same universality class as segregation of simple fluids: the degree of polymerization only slows down the segregation kinetics. For d = 2 polymeric fluids, the domain growth law is L(t) similar to t(phi) with phi showing a crossover from 1/3 -> 1/2 -> 2/3. For d = 3 polymeric fluids, we see the crossover phi = 1/3 -> 1. Our MD simulations do not yet access the inertial hydrodynamic regime (with L similar to t(2/3)) of phase separation in 3-d fluids. (C) 2014 AIP Publishing LLC.