Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lower temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.
Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lowermore » temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.« less
Surface electron density models for accurate ab initio molecular dynamics with electronic friction
NASA Astrophysics Data System (ADS)
Novko, D.; Blanco-Rey, M.; Alducin, M.; Juaristi, J. I.
2016-06-01
Ab initio molecular dynamics with electronic friction (AIMDEF) is a valuable methodology to study the interaction of atomic particles with metal surfaces. This method, in which the effect of low-energy electron-hole (e-h) pair excitations is treated within the local density friction approximation (LDFA) [Juaristi et al., Phys. Rev. Lett. 100, 116102 (2008), 10.1103/PhysRevLett.100.116102], can provide an accurate description of both e-h pair and phonon excitations. In practice, its applicability becomes a complicated task in those situations of substantial surface atoms displacements because the LDFA requires the knowledge at each integration step of the bare surface electron density. In this work, we propose three different methods of calculating on-the-fly the electron density of the distorted surface and we discuss their suitability under typical surface distortions. The investigated methods are used in AIMDEF simulations for three illustrative adsorption cases, namely, dissociated H2 on Pd(100), N on Ag(111), and N2 on Fe(110). Our AIMDEF calculations performed with the three approaches highlight the importance of going beyond the frozen surface density to accurately describe the energy released into e-h pair excitations in case of large surface atom displacements.
Wijma, Hein J; Marrink, Siewert J; Janssen, Dick B
2014-07-28
Computational approaches could decrease the need for the laborious high-throughput experimental screening that is often required to improve enzymes by mutagenesis. Here, we report that using multiple short molecular dynamics (MD) simulations makes it possible to accurately model enantioselectivity for large numbers of enzyme-substrate combinations at low computational costs. We chose four different haloalkane dehalogenases as model systems because of the availability of a large set of experimental data on the enantioselective conversion of 45 different substrates. To model the enantioselectivity, we quantified the frequency of occurrence of catalytically productive conformations (near attack conformations) for pairs of enantiomers during MD simulations. We found that the angle of nucleophilic attack that leads to carbon-halogen bond cleavage was a critical variable that limited the occurrence of productive conformations; enantiomers for which this angle reached values close to 180° were preferentially converted. A cluster of 20-40 very short (10 ps) MD simulations allowed adequate conformational sampling and resulted in much better agreement to experimental enantioselectivities than single long MD simulations (22 ns), while the computational costs were 50-100 fold lower. With single long MD simulations, the dynamics of enzyme-substrate complexes remained confined to a conformational subspace that rarely changed significantly, whereas with multiple short MD simulations a larger diversity of conformations of enzyme-substrate complexes was observed. PMID:24916632
Lee, M.W.; Meuwly, M.
2013-01-01
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories.
Lee, Myung Won; Meuwly, Markus
2013-12-14
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories. PMID:24170171
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-14
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online. PMID:25770527
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-14
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10{sup 3}-10{sup 5} molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.
NASA Astrophysics Data System (ADS)
Schwörer, Magnus; Lorenzen, Konstantin; Mathias, Gerald; Tavan, Paul
2015-03-01
Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 103-105 molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat.
Liu, Jian; Li, Dezhang; Liu, Xinzijian
2016-07-14
We introduce a novel simple algorithm for thermostatting path integral molecular dynamics (PIMD) with the Langevin equation. The staging transformation of path integral beads is employed for demonstration. The optimum friction coefficients for the staging modes in the free particle limit are used for all systems. In comparison to the path integral Langevin equation thermostat, the new algorithm exploits a different order of splitting for the phase space propagator associated to the Langevin equation. While the error analysis is made for both algorithms, they are also employed in the PIMD simulations of three realistic systems (the H2O molecule, liquid para-hydrogen, and liquid water) for comparison. It is shown that the new thermostat increases the time interval of PIMD by a factor of 4-6 or more for achieving the same accuracy. In addition, the supplementary material shows the error analysis made for the algorithms when the normal-mode transformation of path integral beads is used. PMID:27421393
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
NASA Astrophysics Data System (ADS)
Liu, Jian; Li, Dezhang; Liu, Xinzijian
2016-07-01
We introduce a novel simple algorithm for thermostatting path integral molecular dynamics (PIMD) with the Langevin equation. The staging transformation of path integral beads is employed for demonstration. The optimum friction coefficients for the staging modes in the free particle limit are used for all systems. In comparison to the path integral Langevin equation thermostat, the new algorithm exploits a different order of splitting for the phase space propagator associated to the Langevin equation. While the error analysis is made for both algorithms, they are also employed in the PIMD simulations of three realistic systems (the H2O molecule, liquid para-hydrogen, and liquid water) for comparison. It is shown that the new thermostat increases the time interval of PIMD by a factor of 4-6 or more for achieving the same accuracy. In addition, the supplementary material shows the error analysis made for the algorithms when the normal-mode transformation of path integral beads is used.
Lippert, Ross A; Predescu, Cristian; Ierardi, Douglas J; Mackenzie, Kenneth M; Eastwood, Michael P; Dror, Ron O; Shaw, David E
2013-10-28
In molecular dynamics simulations, control over temperature and pressure is typically achieved by augmenting the original system with additional dynamical variables to create a thermostat and a barostat, respectively. These variables generally evolve on timescales much longer than those of particle motion, but typical integrator implementations update the additional variables along with the particle positions and momenta at each time step. We present a framework that replaces the traditional integration procedure with separate barostat, thermostat, and Newtonian particle motion updates, allowing thermostat and barostat updates to be applied infrequently. Such infrequent updates provide a particularly substantial performance advantage for simulations parallelized across many computer processors, because thermostat and barostat updates typically require communication among all processors. Infrequent updates can also improve accuracy by alleviating certain sources of error associated with limited-precision arithmetic. In addition, separating the barostat, thermostat, and particle motion update steps reduces certain truncation errors, bringing the time-average pressure closer to its target value. Finally, this framework, which we have implemented on both general-purpose and special-purpose hardware, reduces software complexity and improves software modularity. PMID:24182003
Chocholousová, Jana; Feig, Michael
2006-04-30
Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs. PMID:16518883
Accurate path integral molecular dynamics simulation of ab-initio water at near-zero added cost
NASA Astrophysics Data System (ADS)
Elton, Daniel; Fritz, Michelle; Soler, José; Fernandez-Serra, Marivi
It is now established that nuclear quantum motion plays an important role in determining water's structure and dynamics. These effects are important to consider when evaluating DFT functionals and attempting to develop better ones for water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads, which is typically 32. Here we introduce a method whereby PIMD can be incorporated into a DFT molecular dynamics simulation at virtually zero cost. The method is based on the cluster (many body) expansion of the energy. We first subtract the DFT monomer energies, using a custom DFT-based monomer potential energy surface. The evolution of the PIMD beads is then performed using only the more-accurate Partridge-Schwenke monomer energy surface. The DFT calculations are done using the centroid positions. Various bead thermostats can be employed to speed up the sampling of the quantum ensemble. The method bears some resemblance to multiple timestep algorithms and other schemes used to speed up PIMD with classical force fields. We show that our method correctly captures some of key effects of nuclear quantum motion on both the structure and dynamics of water. We acknowledge support from DOE Award No. DE-FG02-09ER16052 (D.E.) and DOE Early Career Award No. DE-SC0003871 (M.V.F.S.).
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Filizola, Marta
2009-01-01
For years conventional drug design at G-protein coupled receptors (GPCRs) has mainly focused on the inhibition of a single receptor at a usually well-defined ligand-binding site. The recent discovery of more and more physiologically relevant GPCR dimers/oligomers suggests that selectively targeting these complexes or designing small molecules that inhibit receptor-receptor interactions might provide new opportunities for novel drug discovery. To uncover the fundamental mechanisms and dynamics governing GPCR dimerization/oligomerization, it is crucial to understand the dynamic process of receptor-receptor association, and to identify regions that are suitable for selective drug binding. This minireview highlights current progress in the development of increasingly accurate dynamic molecular models of GPCR oligomers based on structural, biochemical, and biophysical information that has recently appeared in the literature. In view of this new information, there has never been a more exciting time for computational research into GPCRs than at present. Information-driven modern molecular models of GPCR complexes are expected to efficiently guide the rational design of GPCR oligomer-specific drugs, possibly allowing researchers to reach for the high-hanging fruits in GPCR drug discovery, i.e. more potent and selective drugs for efficient therapeutic interventions. PMID:19465029
NASA Astrophysics Data System (ADS)
Yogurtcu, Osman N.; Johnson, Margaret E.
2015-08-01
The dynamics of association between diffusing and reacting molecular species are routinely quantified using simple rate-equation kinetics that assume both well-mixed concentrations of species and a single rate constant for parameterizing the binding rate. In two-dimensions (2D), however, even when systems are well-mixed, the assumption of a single characteristic rate constant for describing association is not generally accurate, due to the properties of diffusional searching in dimensions d ≤ 2. Establishing rigorous bounds for discriminating between 2D reactive systems that will be accurately described by rate equations with a single rate constant, and those that will not, is critical for both modeling and experimentally parameterizing binding reactions restricted to surfaces such as cellular membranes. We show here that in regimes of intrinsic reaction rate (ka) and diffusion (D) parameters ka/D > 0.05, a single rate constant cannot be fit to the dynamics of concentrations of associating species independently of the initial conditions. Instead, a more sophisticated multi-parametric description than rate-equations is necessary to robustly characterize bimolecular reactions from experiment. Our quantitative bounds derive from our new analysis of 2D rate-behavior predicted from Smoluchowski theory. Using a recently developed single particle reaction-diffusion algorithm we extend here to 2D, we are able to test and validate the predictions of Smoluchowski theory and several other theories of reversible reaction dynamics in 2D for the first time. Finally, our results also mean that simulations of reactive systems in 2D using rate equations must be undertaken with caution when reactions have ka/D > 0.05, regardless of the simulation volume. We introduce here a simple formula for an adaptive concentration dependent rate constant for these chemical kinetics simulations which improves on existing formulas to better capture non-equilibrium reaction dynamics from dilute
Bai, Fang; Liao, Sha; Gu, Junfeng; Jiang, Hualiang; Wang, Xicheng; Li, Honglin
2015-04-27
Metalloproteins, particularly zinc metalloproteins, are promising therapeutic targets, and recent efforts have focused on the identification of potent and selective inhibitors of these proteins. However, the ability of current drug discovery and design technologies, such as molecular docking and molecular dynamics simulations, to probe metal-ligand interactions remains limited because of their complicated coordination geometries and rough treatment in current force fields. Herein we introduce a robust, multiobjective optimization algorithm-driven metalloprotein-specific docking program named MpSDock, which runs on a scheme similar to consensus scoring consisting of a force-field-based scoring function and a knowledge-based scoring function. For this purpose, in this study, an effective knowledge-based zinc metalloprotein-specific scoring function based on the inverse Boltzmann law was designed and optimized using a dynamic sampling and iteration optimization strategy. This optimization strategy can dynamically sample and regenerate decoy poses used in each iteration step of refining the scoring function, thus dramatically improving both the effectiveness of the exploration of the binding conformational space and the sensitivity of the ranking of the native binding poses. To validate the zinc metalloprotein-specific scoring function and its special built-in docking program, denoted MpSDockZn, an extensive comparison was performed against six universal, popular docking programs: Glide XP mode, Glide SP mode, Gold, AutoDock, AutoDock4Zn, and EADock DSS. The zinc metalloprotein-specific knowledge-based scoring function exhibited prominent performance in accurately describing the geometries and interactions of the coordination bonds between the zinc ions and chelating agents of the ligands. In addition, MpSDockZn had a competitive ability to sample and identify native binding poses with a higher success rate than the other six docking programs. PMID:25746437
NASA Astrophysics Data System (ADS)
Osei-Kuffuor, Daniel; Fattebert, Jean-Luc
2014-01-01
We present the first truly scalable first-principles molecular dynamics algorithm with O(N) complexity and controllable accuracy, capable of simulating systems with finite band gaps of sizes that were previously impossible with this degree of accuracy. By avoiding global communications, we provide a practical computational scheme capable of extreme scalability. Accuracy is controlled by the mesh spacing of the finite difference discretization, the size of the localization regions in which the electronic wave functions are confined, and a cutoff beyond which the components of the overlap matrix can be omitted when computing selected elements of its inverse. We demonstrate the algorithm's excellent parallel scaling for up to 101 952 atoms on 23 328 processors, with a wall-clock time of the order of 1 min per molecular dynamics time step and numerical error on the forces of less than 7×10-4 Ha/Bohr.
Osei-Kuffuor, Daniel; Fattebert, Jean-Luc
2014-01-01
We present the first truly scalable first-principles molecular dynamics algorithm with O(N) complexity and controllable accuracy, capable of simulating systems with finite band gaps of sizes that were previously impossible with this degree of accuracy. By avoiding global communications, we provide a practical computational scheme capable of extreme scalability. Accuracy is controlled by the mesh spacing of the finite difference discretization, the size of the localization regions in which the electronic wave functions are confined, and a cutoff beyond which the components of the overlap matrix can be omitted when computing selected elements of its inverse. We demonstrate the algorithm's excellent parallel scaling for up to 101 952 atoms on 23 328 processors, with a wall-clock time of the order of 1 min per molecular dynamics time step and numerical error on the forces of less than 7x10^{-4} Ha/Bohr.
Ladd, A.J.C.
1988-08-01
The basic methodology of equilibrium molecular dynamics is described. Examples from the literature are used to illustrate how molecular dynamics has been used to resolve theoretical controversies, provide data to test theories, and occasionally to discover new phenomena. The emphasis is on the application of molecular dynamics to an understanding of the microscopic physics underlying the transport properties of simple fluids. 98 refs., 4 figs.
NASA Astrophysics Data System (ADS)
Osei-Kuffuor, Daniel; Fattebert, Jean-Luc
2014-03-01
We present a truly scalable First-Principles Molecular Dynamics algorithm with O(N) complexity and fully controllable accuracy, capable of simulating systems of sizes that were previously impossible with this degree of accuracy. By avoiding global communication, we have extended W. Kohn's condensed matter ``nearsightedness'' principle to a practical computational scheme capable of extreme scalability. Accuracy is controlled by the mesh spacing of the finite difference discretization, the size of the localization regions in which the electronic wavefunctions are confined, and a cutoff beyond which the components of the overlap matrix can be omitted when computing selected elements of its inverse. We demonstrate the algorithm's excellent parallel scaling for up to 100,000 atoms on 100,000 processors, with a wall-clock time of the order of one minute per molecular dynamics time step. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Sangiovanni, D. G.; Hellman, O.; Alling, B.; Abrikosov, I. A.
2016-03-01
We revisit the color-diffusion algorithm [Aeberhard et al., Phys. Rev. Lett. 108, 095901 (2012), 10.1103/PhysRevLett.108.095901] in non equilibrium ab initio molecular dynamics (NE-AIMD) and propose a simple efficient approach for the estimation of monovacancy jump rates in crystalline solids at temperatures well below melting. Color-diffusion applied to monovacancy migration entails that one lattice atom (colored atom) is accelerated toward the neighboring defect site by an external constant force F. Considering bcc molybdenum between 1000 and 2800 K as a model system, NE-AIMD results show that the colored-atom jump rate kNE increases exponentially with the force intensity F , up to F values far beyond the linear-fitting regime employed previously. Using a simple model, we derive an analytical expression which reproduces the observed kNE(F ) dependence on F . Equilibrium rates extrapolated by NE-AIMD results are in excellent agreement with those of unconstrained dynamics. The gain in computational efficiency achieved with our approach increases rapidly with decreasing temperatures and reaches a factor of 4 orders of magnitude at the lowest temperature considered in the present study.
NASA Astrophysics Data System (ADS)
Genova, Alessandro; Ceresoli, Davide; Pavanello, Michele
2016-06-01
In this work we achieve three milestones: (1) we present a subsystem DFT method capable of running ab-initio molecular dynamics simulations accurately and efficiently. (2) In order to rid the simulations of inter-molecular self-interaction error, we exploit the ability of semilocal frozen density embedding formulation of subsystem DFT to represent the total electron density as a sum of localized subsystem electron densities that are constrained to integrate to a preset, constant number of electrons; the success of the method relies on the fact that employed semilocal nonadditive kinetic energy functionals effectively cancel out errors in semilocal exchange-correlation potentials that are linked to static correlation effects and self-interaction. (3) We demonstrate this concept by simulating liquid water and solvated OH• radical. While the bulk of our simulations have been performed on a periodic box containing 64 independent water molecules for 52 ps, we also simulated a box containing 256 water molecules for 22 ps. The results show that, provided one employs an accurate nonadditive kinetic energy functional, the dynamics of liquid water and OH• radical are in semiquantitative agreement with experimental results or higher-level electronic structure calculations. Our assessments are based upon comparisons of radial and angular distribution functions as well as the diffusion coefficient of the liquid.
NASA Astrophysics Data System (ADS)
Meng, Qingyong; Chen, Jun; Zhang, Dong H.
2015-09-01
The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMD is capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures.
Willow, Soohaeng Yoo; Salim, Michael A; Kim, Kwang S; Hirata, So
2015-01-01
A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690
Willow, Soohaeng Yoo; Salim, Michael A.; Kim, Kwang S.; Hirata, So
2015-01-01
A direct, simultaneous calculation of properties of a liquid using an ab initio electron-correlated theory has long been unthinkable. Here we present structural, dynamical, and response properties of liquid water calculated by ab initio molecular dynamics using the embedded-fragment spin-component-scaled second-order many-body perturbation method with the aug-cc-pVDZ basis set. This level of theory is chosen as it accurately and inexpensively reproduces the water dimer potential energy surface from the coupled-cluster singles, doubles, and noniterative triples with the aug-cc-pVQZ basis set, which is nearly exact. The calculated radial distribution function, self-diffusion coefficient, coordinate number, and dipole moment, as well as the infrared and Raman spectra are in excellent agreement with experimental results. The shapes and widths of the OH stretching bands in the infrared and Raman spectra and their isotropic-anisotropic Raman noncoincidence, which reflect the diverse local hydrogen-bond environment, are also reproduced computationally. The simulation also reveals intriguing dynamic features of the environment, which are difficult to probe experimentally, such as a surprisingly large fluctuation in the coordination number and the detailed mechanism by which the hydrogen donating water molecules move across the first and second shells, thereby causing this fluctuation. PMID:26400690
Accurate Molecular Polarizabilities Based on Continuum Electrostatics
Truchon, Jean-François; Nicholls, Anthony; Iftimie, Radu I.; Roux, Benoît; Bayly, Christopher I.
2013-01-01
A novel approach for representing the intramolecular polarizability as a continuum dielectric is introduced to account for molecular electronic polarization. It is shown, using a finite-difference solution to the Poisson equation, that the Electronic Polarization from Internal Continuum (EPIC) model yields accurate gas-phase molecular polarizability tensors for a test set of 98 challenging molecules composed of heteroaromatics, alkanes and diatomics. The electronic polarization originates from a high intramolecular dielectric that produces polarizabilities consistent with B3LYP/aug-cc-pVTZ and experimental values when surrounded by vacuum dielectric. In contrast to other approaches to model electronic polarization, this simple model avoids the polarizability catastrophe and accurately calculates molecular anisotropy with the use of very few fitted parameters and without resorting to auxiliary sites or anisotropic atomic centers. On average, the unsigned error in the average polarizability and anisotropy compared to B3LYP are 2% and 5%, respectively. The correlation between the polarizability components from B3LYP and this approach lead to a R2 of 0.990 and a slope of 0.999. Even the F2 anisotropy, shown to be a difficult case for existing polarizability models, can be reproduced within 2% error. In addition to providing new parameters for a rapid method directly applicable to the calculation of polarizabilities, this work extends the widely used Poisson equation to areas where accurate molecular polarizabilities matter. PMID:23646034
NASA Astrophysics Data System (ADS)
Meng, Qingyong; Chen, Jun; Zhang, Dong H.
2016-04-01
To fast and accurately compute rate coefficients of the H/D + CH4 → H2/HD + CH3 reactions, we propose a segmented strategy for fitting suitable potential energy surface (PES), on which ring-polymer molecular dynamics (RPMD) simulations are performed. On the basis of recently developed permutation invariant polynomial neural-network approach [J. Li et al., J. Chem. Phys. 142, 204302 (2015)], PESs in local configuration spaces are constructed. In this strategy, global PES is divided into three parts, including asymptotic, intermediate, and interaction parts, along the reaction coordinate. Since less fitting parameters are involved in the local PESs, the computational efficiency for operating the PES routine is largely enhanced by a factor of ˜20, comparing with that for global PES. On interaction part, the RPMD computational time for the transmission coefficient can be further efficiently reduced by cutting off the redundant part of the child trajectories. For H + CH4, good agreements among the present RPMD rates and those from previous simulations as well as experimental results are found. For D + CH4, on the other hand, qualitative agreement between present RPMD and experimental results is predicted.
Multiscale reactive molecular dynamics
NASA Astrophysics Data System (ADS)
Knight, Chris; Lindberg, Gerrick E.; Voth, Gregory A.
2012-12-01
Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system.
Multiscale reactive molecular dynamics
Knight, Chris; Lindberg, Gerrick E.; Voth, Gregory A.
2012-01-01
Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system. PMID:23249062
Developing accurate molecular mechanics force fields for conjugated molecular systems.
Do, Hainam; Troisi, Alessandro
2015-10-14
A rapid method to parameterize the intramolecular component of classical force fields for complex conjugated molecules is proposed. The method is based on a procedure of force matching with a reference electronic structure calculation. It is particularly suitable for those applications where molecular dynamics simulations are used to generate structures that are therefore analysed by electronic structure methods, because it is possible to build force fields that are consistent with electronic structure calculations that follow classical simulations. Such applications are commonly encountered in organic electronics, spectroscopy of complex systems and photobiology (e.g. photosynthetic systems). We illustrate the method by parameterizing the force fields of a molecule used in molecular semiconductors (2,2-dicyanovinyl-capped S,N-heteropentacene or DCV-SN5), a polymeric semiconductor (thieno[3,2-b]thiophene-diketopyrrolopyrrole TT-DPP) and a chromophore embedded in a protein environment (15,16-dihydrobiliverdin or DBV) where several hundreds of parameters need to be optimized in parallel. PMID:26349916
Accelerated molecular dynamics methods
Perez, Danny
2011-01-04
The molecular dynamics method, although extremely powerful for materials simulations, is limited to times scales of roughly one microsecond or less. On longer time scales, dynamical evolution typically consists of infrequent events, which are usually activated processes. This course is focused on understanding infrequent-event dynamics, on methods for characterizing infrequent-event mechanisms and rate constants, and on methods for simulating long time scales in infrequent-event systems, emphasizing the recently developed accelerated molecular dynamics methods (hyperdynamics, parallel replica dynamics, and temperature accelerated dynamics). Some familiarity with basic statistical mechanics and molecular dynamics methods will be assumed.
Le, Hung M; Dinh, Thach S; Le, Hieu V
2011-10-13
The singlet-triplet transformation and molecular dissociation of ozone (O(3)) gas is investigated by performing quasi-classical molecular dynamics (MD) simulations on an ab initio potential energy surface (PES) with visible and near-infrared excitations. MP4(SDQ) level of theory with the 6-311g(2d,2p) basis set is executed for three different electronic spin states (singlet, triplet, and quintet). In order to simplify the potential energy function, an approximation is adopted by ignoring the spin-orbit coupling and allowing the molecule to switch favorably and instantaneously to the spin state that is more energetically stable (lowest in energy among the three spin states). This assumption has previously been utilized to study the SiO(2) system as reported by Agrawal et al. (J. Chem. Phys. 2006, 124 (13), 134306). The use of such assumption in this study probably makes the upper limits of computed rate coefficients the true rate coefficients. The global PES for ozone is constructed by fitting 5906 ab initio data points using a 60-neuron two-layer feed-forward neural network. The mean-absolute error and root-mean-squared error of this fit are 0.0446 eV (1.03 kcal/mol) and 0.0756 eV (1.74 kcal/mol), respectively, which reveal very good fitting accuracy. The parameter coefficients of the global PES are reported in this paper. In order to identify the spin state with high confidence, we propose the use of a pattern-recognition neural network, which is trained to predict the spin state of a given configuration (with a prediction accuracy being 95.6% on a set of testing data points). To enhance the prediction effectiveness, a buffer series of five points are validated to confirm the spin state during the MD process to gain better confidence. Quasi-classical MD simulations from 1.2 to 2.4 eV of total internal energy (including zero-point energy) result in rate coefficients of singlet-triplet transformation in the range of 0.027 ps(-1) to 1.21 ps(-1). Also, we find very
Towards Accurate Molecular Modeling of Plastic Bonded Explosives
NASA Astrophysics Data System (ADS)
Chantawansri, T. L.; Andzelm, J.; Taylor, D.; Byrd, E.; Rice, B.
2010-03-01
There is substantial interest in identifying the controlling factors that influence the susceptibility of polymer bonded explosives (PBXs) to accidental initiation. Numerous Molecular Dynamics (MD) simulations of PBXs using the COMPASS force field have been reported in recent years, where the validity of the force field in modeling the solid EM fill has been judged solely on its ability to reproduce lattice parameters, which is an insufficient metric. Performance of the COMPASS force field in modeling EMs and the polymeric binder has been assessed by calculating structural, thermal, and mechanical properties, where only fair agreement with experimental data is obtained. We performed MD simulations using the COMPASS force field for the polymer binder hydroxyl-terminated polybutadiene and five EMs: cyclotrimethylenetrinitramine, 1,3,5,7-tetranitro-1,3,5,7-tetra-azacyclo-octane, 2,4,6,8,10,12-hexantirohexaazazisowurzitane, 2,4,6-trinitro-1,3,5-benzenetriamine, and pentaerythritol tetranitate. Predicted EM crystallographic and molecular structural parameters, as well as calculated properties for the binder will be compared with experimental results for different simulation conditions. We also present novel simulation protocols, which improve agreement between experimental and computation results thus leading to the accurate modeling of PBXs.
NASA Astrophysics Data System (ADS)
Kapil, V.; VandeVondele, J.; Ceriotti, M.
2016-02-01
The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modeling the quantum nature of the nuclei, the cost is still prohibitive when combined with advanced electronic structure methods. Here we present how multiple time step integrators can be combined with ring-polymer contraction techniques (effectively, multiple time stepping in imaginary time) to reduce virtually to zero the overhead of modelling nuclear quantum effects, while describing inter-atomic forces at high levels of electronic structure theory. This is demonstrated for a combination of MP2 and semi-local DFT applied to the Zundel cation. The approach can be seamlessly combined with other methods to reduce the computational cost of path integral calculations, such as high-order factorizations of the Boltzmann operator or generalized Langevin equation thermostats.
Integration methods for molecular dynamics
Leimkuhler, B.J.; Reich, S.; Skeel, R.D.
1996-12-31
Classical molecular dynamics simulation of a macromolecule requires the use of an efficient time-stepping scheme that can faithfully approximate the dynamics over many thousands of timesteps. Because these problems are highly nonlinear, accurate approximation of a particular solution trajectory on meaningful time intervals is neither obtainable nor desired, but some restrictions, such as symplecticness, can be imposed on the discretization which tend to imply good long term behavior. The presence of a variety of types and strengths of interatom potentials in standard molecular models places severe restrictions on the timestep for numerical integration used in explicit integration schemes, so much recent research has concentrated on the search for alternatives that possess (1) proper dynamical properties, and (2) a relative insensitivity to the fastest components of the dynamics. We survey several recent approaches. 48 refs., 2 figs.
Meng, Qingyong Chen, Jun Zhang, Dong H.
2015-09-14
The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMD is capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures.
Kapil, V; VandeVondele, J; Ceriotti, M
2016-02-01
The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modeling the quantum nature of the nuclei, the cost is still prohibitive when combined with advanced electronic structure methods. Here we present how multiple time step integrators can be combined with ring-polymer contraction techniques (effectively, multiple time stepping in imaginary time) to reduce virtually to zero the overhead of modelling nuclear quantum effects, while describing inter-atomic forces at high levels of electronic structure theory. This is demonstrated for a combination of MP2 and semi-local DFT applied to the Zundel cation. The approach can be seamlessly combined with other methods to reduce the computational cost of path integral calculations, such as high-order factorizations of the Boltzmann operator or generalized Langevin equation thermostats. PMID:26851912
Nonequilibrium molecular dynamics
Hoover, W.G. . Dept. of Applied Science Lawrence Livermore National Lab., CA )
1990-11-01
The development of nonequilibrium molecular dynamics is described, with emphasis on massively-parallel simulations involving the motion of millions, soon to be billions, of atoms. Corresponding continuum simulations are also discussed. 14 refs., 8 figs.
A sampling of molecular dynamics
NASA Astrophysics Data System (ADS)
Sindhikara, Daniel Jon
The sheer vastness of the number of computations required to simulate a biological molecule puts incredible pressure on algorithms to be efficient while maintaining sufficient accuracy. This dissertation summarizes various projects whose purposes address the large span of types of problems in molecular dynamics simulations of biological systems including: increasing efficiency, measuring convergence, avoiding pitfalls, and an application and analysis of a biological system. Chapters 3 and 4 deal with an enhanced sampling algorithm called "replica exchange molecular dynamics" which is designed to speed-up molecular dynamics simulations. The optimization of a key parameter of these simulations is analyzed. In these successive projects, it was found conclusively that maximizing "exchange attempt frequency" is the most efficient way to run a replica exchange molecular dynamics simulation. Chapter 5 describes an enhanced metric for convergence in parallel simulations called the normalized ergodic measure. The metric is applied to several properties for several replica exchange simulations. Advantages of this metric over other methods are described. Chapter 6 describes the implementation and optimization of an enhanced sampling algorithm similar to replica exchange molecular dynamics called multicanonical algorithm replica exchange molecular dynamics. The algorithm was implemented into a biomolecular simulation suite called AMBER. Additionally several parameters were analyzed and optimized. In Chapter 7, a pitfall in molecular dynamics is observed in biological systems that is caused by negligent use of a simulation's "thermostat". It was found that if the same pseudorandom number seed were used for multiple systems, they eventually synchronize. In this project, synchronization was observed in biological molecules. Various negative effects including corruption of data are pointed out. Chapter 8 describes molecular dynamics simulation of NikR, a homotetrameric nickel
Substructured multibody molecular dynamics.
Grest, Gary Stephen; Stevens, Mark Jackson; Plimpton, Steven James; Woolf, Thomas B. (Johns Hopkins University, Baltimore, MD); Lehoucq, Richard B.; Crozier, Paul Stewart; Ismail, Ahmed E.; Mukherjee, Rudranarayan M. (Rensselaer Polytechnic Institute, Troy, NY); Draganescu, Andrei I.
2006-11-01
We have enhanced our parallel molecular dynamics (MD) simulation software LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, lammps.sandia.gov) to include many new features for accelerated simulation including articulated rigid body dynamics via coupling to the Rensselaer Polytechnic Institute code POEMS (Parallelizable Open-source Efficient Multibody Software). We use new features of the LAMMPS software package to investigate rhodopsin photoisomerization, and water model surface tension and capillary waves at the vapor-liquid interface. Finally, we motivate the recipes of MD for practitioners and researchers in numerical analysis and computational mechanics.
Superposition State Molecular Dynamics.
Venkatnathan, Arun; Voth, Gregory A
2005-01-01
The ergodic sampling of rough energy landscapes is crucial for understanding phenomena like protein folding, peptide aggregation, polymer dynamics, and the glass transition. These rough energy landscapes are characterized by the presence of many local minima separated by high energy barriers, where Molecular Dynamics (MD) fails to satisfy ergodicity. To enhance ergodic behavior, we have developed the Superposition State Molecular Dynamics (SSMD) method, which uses a superposition of energy states to obtain an effective potential for the MD simulation. In turn, the dynamics on this effective potential can be used to sample the configurational free energy of the real potential. The effectiveness of the SSMD method for a one-dimensional rough potential energy landscape is presented as a test case. PMID:26641113
Open boundary molecular dynamics
NASA Astrophysics Data System (ADS)
Delgado-Buscalioni, R.; Sablić, J.; Praprotnik, M.
2015-09-01
This contribution analyzes several strategies and combination of methodologies to perform molecular dynamic simulations in open systems. Here, the term open indicates that the total system has boundaries where transfer of mass, momentum and energy can take place. This formalism, which we call Open Boundary Molecular Dynamics (OBMD), can act as interface of different schemes, such as Adaptive Resolution Scheme (AdResS) and Hybrid continuum-particle dynamics to link atomistic, coarse-grained (CG) and continuum (Eulerian) fluid dynamics in the general framework of fluctuating Navier-Stokes equations. The core domain of the simulation box is solved using all-atom descriptions. The CG layer introduced using AdResS is located at the outer part of the open box to make feasible the insertion of large molecules into the system. Communications between the molecular system and the outer world are carried out in the outer layers, called buffers. These coupling preserve momentum and mass conservation laws and can thus be linked with Eulerian hydro- dynamic solvers. In its simpler form, OBMD allows, however, to impose a local pressure tensor and a heat flux across the system's boundaries. For a one component molecular system, the external normal pressure and temperature determine the external chemical potential and thus the independent parameters of a grand-canonical ensemble simulation. Extended ensembles under non-equilibrium stationary states can also be simulated as well as time dependent forcings (e.g. oscillatory rheology). To illustrate the robustness of the combined OBMD-AdResS method, we present simulations of star-polymer melts at equilibrium and in sheared flow.
A robust and accurate formulation of molecular and colloidal electrostatics.
Sun, Qiang; Klaseboer, Evert; Chan, Derek Y C
2016-08-01
This paper presents a re-formulation of the boundary integral method for the Debye-Hückel model of molecular and colloidal electrostatics that removes the mathematical singularities that have to date been accepted as an intrinsic part of the conventional boundary integral equation method. The essence of the present boundary regularized integral equation formulation consists of subtracting a known solution from the conventional boundary integral method in such a way as to cancel out the singularities associated with the Green's function. This approach better reflects the non-singular physical behavior of the systems on boundaries with the benefits of the following: (i) the surface integrals can be evaluated accurately using quadrature without any need to devise special numerical integration procedures, (ii) being able to use quadratic or spline function surface elements to represent the surface more accurately and the variation of the functions within each element is represented to a consistent level of precision by appropriate interpolation functions, (iii) being able to calculate electric fields, even at boundaries, accurately and directly from the potential without having to solve hypersingular integral equations and this imparts high precision in calculating the Maxwell stress tensor and consequently, intermolecular or colloidal forces, (iv) a reliable way to handle geometric configurations in which different parts of the boundary can be very close together without being affected by numerical instabilities, therefore potentials, fields, and forces between surfaces can be found accurately at surface separations down to near contact, and (v) having the simplicity of a formulation that does not require complex algorithms to handle singularities will result in significant savings in coding effort and in the reduction of opportunities for coding errors. These advantages are illustrated using examples drawn from molecular and colloidal electrostatics. PMID:27497538
A robust and accurate formulation of molecular and colloidal electrostatics
NASA Astrophysics Data System (ADS)
Sun, Qiang; Klaseboer, Evert; Chan, Derek Y. C.
2016-08-01
This paper presents a re-formulation of the boundary integral method for the Debye-Hückel model of molecular and colloidal electrostatics that removes the mathematical singularities that have to date been accepted as an intrinsic part of the conventional boundary integral equation method. The essence of the present boundary regularized integral equation formulation consists of subtracting a known solution from the conventional boundary integral method in such a way as to cancel out the singularities associated with the Green's function. This approach better reflects the non-singular physical behavior of the systems on boundaries with the benefits of the following: (i) the surface integrals can be evaluated accurately using quadrature without any need to devise special numerical integration procedures, (ii) being able to use quadratic or spline function surface elements to represent the surface more accurately and the variation of the functions within each element is represented to a consistent level of precision by appropriate interpolation functions, (iii) being able to calculate electric fields, even at boundaries, accurately and directly from the potential without having to solve hypersingular integral equations and this imparts high precision in calculating the Maxwell stress tensor and consequently, intermolecular or colloidal forces, (iv) a reliable way to handle geometric configurations in which different parts of the boundary can be very close together without being affected by numerical instabilities, therefore potentials, fields, and forces between surfaces can be found accurately at surface separations down to near contact, and (v) having the simplicity of a formulation that does not require complex algorithms to handle singularities will result in significant savings in coding effort and in the reduction of opportunities for coding errors. These advantages are illustrated using examples drawn from molecular and colloidal electrostatics.
Molecular Dynamics of Acetylcholinesterase
Shen, T Y.; Tai, Kaihsu; Henchman, Richard H.; Mccammon, Andy
2002-06-01
Molecular dynamics simulations are leading to a deeper understanding of the activity of the enzyme acetylcholinesterase. Simulations have shown how breathing motions in the enzyme facilitate the displacement of substrate from the surface of the enzyme to the buried active site. The most recent work points to the complex and spatially extensive nature of such motions and suggests possible modes of regulation of the activity of the enzyme.
A general, accurate procedure for calculating molecular interaction force.
Yang, Pinghai; Qian, Xiaoping
2009-09-15
The determination of molecular interaction forces, e.g., van der Waals force, between macroscopic bodies is of fundamental importance for understanding sintering, adhesion and fracture processes. In this paper, we develop an accurate, general procedure for van der Waals force calculation. This approach extends a surface formulation that converts a six-dimensional (6D) volume integral into a 4D surface integral for the force calculation. It uses non-uniform rational B-spline (NURBS) surfaces to represent object surfaces. Surface integrals are then done on the parametric domain of the NURBS surfaces. It has combined advantages of NURBS surface representation and surface formulation, including (1) molecular interactions between arbitrary-shaped objects can be represented and evaluated by the NURBS model further common geometries such as spheres, cones, planes can be represented exactly and interaction forces are thus calculated accurately; (2) calculation efficiency is improved by converting the volume integral to the surface integral. This approach is implemented and validated via its comparison with analytical solutions for simple geometries. Calculation of van der Waals force between complex geometries with surface roughness is also demonstrated. A tutorial on the NURBS approach is given in Appendix A. PMID:19596335
Accurate Evaluation Method of Molecular Binding Affinity from Fluctuation Frequency
NASA Astrophysics Data System (ADS)
Hoshino, Tyuji; Iwamoto, Koji; Ode, Hirotaka; Ohdomari, Iwao
2008-05-01
Exact estimation of the molecular binding affinity is significantly important for drug discovery. The energy calculation is a direct method to compute the strength of the interaction between two molecules. This energetic approach is, however, not accurate enough to evaluate a slight difference in binding affinity when distinguishing a prospective substance from dozens of candidates for medicine. Hence more accurate estimation of drug efficacy in a computer is currently demanded. Previously we proposed a concept of estimating molecular binding affinity, focusing on the fluctuation at an interface between two molecules. The aim of this paper is to demonstrate the compatibility between the proposed computational technique and experimental measurements, through several examples for computer simulations of an association of human immunodeficiency virus type-1 (HIV-1) protease and its inhibitor (an example for a drug-enzyme binding), a complexation of an antigen and its antibody (an example for a protein-protein binding), and a combination of estrogen receptor and its ligand chemicals (an example for a ligand-receptor binding). The proposed affinity estimation has proven to be a promising technique in the advanced stage of the discovery and the design of drugs.
Multisurface Adiabatic Reactive Molecular Dynamics.
Nagy, Tibor; Yosa Reyes, Juvenal; Meuwly, Markus
2014-04-01
Adiabatic reactive molecular dynamics (ARMD) simulation method is a surface-crossing algorithm for modeling chemical reactions in classical molecular dynamics simulations using empirical force fields. As the ARMD Hamiltonian is time dependent during crossing, it allows only approximate energy conservation. In the current work, the range of applicability of conventional ARMD is explored, and a new multisurface ARMD (MS-ARMD) method is presented, implemented in CHARMM and applied to the vibrationally induced photodissociation of sulfuric acid (H2SO4) in the gas phase. For this, an accurate global potential energy surface (PES) involving 12 H2SO4 and 4 H2O + SO3 force fields fitted to MP2/6-311G++(2d,2p) reference energies is employed. The MS-ARMD simulations conserve total energy and feature both intramolecular H-transfer reactions and water elimination. An analytical treatment of the dynamics in the crossing region finds that conventional ARMD can approximately conserve total energy for limiting cases. In one of them, the reduced mass of the system is large, which often occurs for simulations of solvated biomolecular systems. On the other hand, MS-ARMD is a general approach for modeling chemical reactions including gas-phase, homogeneous, heterogeneous, and enzymatic catalytic reactions while conserving total energy in atomistic simulations. PMID:26580356
Interactive molecular dynamics
NASA Astrophysics Data System (ADS)
Schroeder, Daniel V.
2015-03-01
Physics students now have access to interactive molecular dynamics simulations that can model and animate the motions of hundreds of particles, such as noble gas atoms, that attract each other weakly at short distances but repel strongly when pressed together. Using these simulations, students can develop an understanding of forces and motions at the molecular scale, nonideal fluids, phases of matter, thermal equilibrium, nonequilibrium states, the Boltzmann distribution, the arrow of time, and much more. This article summarizes the basic features and capabilities of such a simulation, presents a variety of student exercises using it at the introductory and intermediate levels, and describes some enhancements that can further extend its uses. A working simulation code, in html5 and javascript for running within any modern Web browser, is provided as an online supplement.
Microsatellites Are Molecular Clocks That Support Accurate Inferences about History
Mullikin, James C.; Patterson, Nick; Reich, David E.
2009-01-01
Microsatellite length mutations are often modeled using the generalized stepwise mutation process, which is a type of random walk. If this model is sufficiently accurate, one can estimate the coalescence time between alleles of a locus after a mathematical transformation of the allele lengths. When large-scale microsatellite genotyping first became possible, there was substantial interest in using this approach to make inferences about time and demography, but that interest has waned because it has not been possible to empirically validate the clock by comparing it with data in which the mutation process is well understood. We analyzed data from 783 microsatellite loci in human populations and 292 loci in chimpanzee populations, and compared them with up to one gigabase of aligned sequence data, where the molecular clock based upon nucleotide substitutions is believed to be reliable. We empirically demonstrate a remarkable linearity (r2 > 0.95) between the microsatellite average square distance statistic and sequence divergence. We demonstrate that microsatellites are accurate molecular clocks for coalescent times of at least 2 million years (My). We apply this insight to confirm that the African populations San, Biaka Pygmy, and Mbuti Pygmy have the deepest coalescent times among populations in the Human Genome Diversity Project. Furthermore, we show that microsatellites support unbiased estimates of population differentiation (FST) that are less subject to ascertainment bias than single nucleotide polymorphism (SNP) FST. These results raise the prospect of using microsatellite data sets to determine parameters of population history. When genotyped along with SNPs, microsatellite data can also be used to correct for SNP ascertainment bias. PMID:19221007
Accurate dynamics in an azimuthally-symmetric accelerating cavity
NASA Astrophysics Data System (ADS)
Appleby, R. B.; Abell, D. T.
2015-02-01
We consider beam dynamics in azimuthally-symmetric accelerating cavities, using the EMMA FFAG cavity as an example. By fitting a vector potential to the field map, we represent the linear and non-linear dynamics using truncated power series and mixed-variable generating functions. The analysis provides an accurate model for particle trajectories in the cavity, reveals potentially significant and measurable effects on the dynamics, and shows differences between cavity focusing models. The approach provides a unified treatment of transverse and longitudinal motion, and facilitates detailed map-based studies of motion in complex machines like FFAGs.
Introduction to Accelerated Molecular Dynamics
Perez, Danny
2012-07-10
Molecular Dynamics is the numerical solution of the equations of motion of a set of atoms, given an interatomic potential V and some boundary and initial conditions. Molecular Dynamics is the largest scale model that gives unbiased dynamics [x(t),p(t)] in full atomistic detail. Molecular Dynamics: is simple; is 'exact' for classical dynamics (with respect to a given V); can be used to compute any (atomistic) thermodynamical or dynamical properties; naturally handles complexity -- the system does the right thing at the right time. The physics derives only from the interatomic potential.
Molecular Dynamics Calculations
NASA Technical Reports Server (NTRS)
1996-01-01
The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two
Molecular adsorption at Pt(111). How accurate are DFT functionals?
Gautier, Sarah; Steinmann, Stephan N; Michel, Carine; Fleurat-Lessard, Paul; Sautet, Philippe
2015-11-21
Molecular chemisorption at a metal surface is a key step for many processes, such as catalysis, electrochemistry, surface treatment, tribology and friction. Modeling with density functional theory is largely used on these systems. From a detailed comparison with accurate micro-calorimetric data on ten systems (involving ethylene, cyclohexene, benzene, naphthalene, CO, O2, H2, methane, ethane), we study the accuracy, for chemisorption on Pt(111), of five exchange-correlation functionals including one generalized gradient approximation functional (PBE) and four functionals that take into account van der Waals interactions (optPBE-vdW, optB86b-vdW, BEEF-vdW, PBE-dDsC). If the functionals used provide very similar geometries and electronic structures, as shown by projected density of states, they give strikingly different results for the adsorption energy of molecules on Pt(111). Among the set of chemisorption data, the lowest mean absolute deviations (MAD) are obtained with the optPBE-vdW and PBE-dDsC functionals (∼0.2 eV) while PBE and optB86b-vdW give twice larger MAD (∼0.45 eV). BEEF-vdW is intermediate with a MAD of 0.33 eV. For laterally π-bound unsaturated hydrocarbons (cyclohexene, benzene, naphthalene) the PBE and the BEEF-vdW functionals are severally under-bound, while optPBE-vdW and PBE-dDsC provide a good match with experiments. Hence both the incorporation of van der Waals dispersive forces and the choice of the exchange functional have a key influence on the chemisorption energy. Vertically bound ethylidyne and CO are in contrast over-bound with all functionals, the best agreement being obtained with BEEF-vdW. None of the selected functionals hence provides a universally accurate treatment of chemisorption energies. PMID:26455444
Molecular dynamics simulations
Alder, B.J.
1985-07-01
The molecular dynamics computer simulation discovery of the slow decay of the velocity autocorrelation function in fluids is briefly reviewed in order to contrast that long time tail with those observed for the stress autocorrelation function in fluids and the velocity autocorrelation function in the Lorentz gas. For a non-localized particle in the Lorentz gas it is made plausible that even if it behaved quantum mechanically its long time tail would be the same as the classical one. The generalization of Fick's law for diffusion for the Lorentz gas, necessary to avoid divergences due to the slow decay of correlations, is presented. For fluids, that generalization has not yet been established, but the region of validity of generalized hydrodynamics is discussed. 20 refs., 5 figs.
VMD: visual molecular dynamics.
Humphrey, W; Dalke, A; Schulten, K
1996-02-01
VMD is a molecular graphics program designed for the display and analysis of molecular assemblies, in particular biopolymers such as proteins and nucleic acids. VMD can simultaneously display any number of structures using a wide variety of rendering styles and coloring methods. Molecules are displayed as one or more "representations," in which each representation embodies a particular rendering method and coloring scheme for a selected subset of atoms. The atoms displayed in each representation are chosen using an extensive atom selection syntax, which includes Boolean operators and regular expressions. VMD provides a complete graphical user interface for program control, as well as a text interface using the Tcl embeddable parser to allow for complex scripts with variable substitution, control loops, and function calls. Full session logging is supported, which produces a VMD command script for later playback. High-resolution raster images of displayed molecules may be produced by generating input scripts for use by a number of photorealistic image-rendering applications. VMD has also been expressly designed with the ability to animate molecular dynamics (MD) simulation trajectories, imported either from files or from a direct connection to a running MD simulation. VMD is the visualization component of MDScope, a set of tools for interactive problem solving in structural biology, which also includes the parallel MD program NAMD, and the MDCOMM software used to connect the visualization and simulation programs. VMD is written in C++, using an object-oriented design; the program, including source code and extensive documentation, is freely available via anonymous ftp and through the World Wide Web. PMID:8744570
NASA Astrophysics Data System (ADS)
Heaps, Charles W.; Mazziotti, David A.
2016-08-01
Quantum molecular dynamics requires an accurate representation of the molecular potential energy surface from a minimal number of electronic structure calculations, particularly for nonadiabatic dynamics where excited states are required. In this paper, we employ pseudospectral sampling of time-dependent Gaussian basis functions for the simulation of non-adiabatic dynamics. Unlike other methods, the pseudospectral Gaussian molecular dynamics tests the Schrödinger equation with N Dirac delta functions located at the centers of the Gaussian functions reducing the scaling of potential energy evaluations from O ( N 2 ) to O ( N ) . By projecting the Gaussian basis onto discrete points in space, the method is capable of efficiently and quantitatively describing the nonadiabatic population transfer and intra-surface quantum coherence. We investigate three model systems: the photodissociation of three coupled Morse oscillators, the bound state dynamics of two coupled Morse oscillators, and a two-dimensional model for collinear triatomic vibrational dynamics. In all cases, the pseudospectral Gaussian method is in quantitative agreement with numerically exact calculations. The results are promising for nonadiabatic molecular dynamics in molecular systems where strongly correlated ground or excited states require expensive electronic structure calculations.
Floating orbital molecular dynamics simulations.
Perlt, Eva; Brüssel, Marc; Kirchner, Barbara
2014-04-21
We introduce an alternative ab initio molecular dynamics simulation as a unification of Hartree-Fock molecular dynamics and the floating orbital approach. The general scheme of the floating orbital molecular dynamics method is presented. Moreover, a simple but sophisticated guess for the orbital centers is provided to reduce the number of electronic structure optimization steps at each molecular dynamics step. The conservation of total energy and angular momentum is investigated in order to validate the floating orbital molecular dynamics approach with and without application of the initial guess. Finally, a water monomer and a water dimer are simulated, and the influence of the orbital floating on certain properties like the dipole moment is investigated. PMID:24600690
Accurate measurements of dynamics and reproducibility in small genetic networks
Dubuis, Julien O; Samanta, Reba; Gregor, Thomas
2013-01-01
Quantification of gene expression has become a central tool for understanding genetic networks. In many systems, the only viable way to measure protein levels is by immunofluorescence, which is notorious for its limited accuracy. Using the early Drosophila embryo as an example, we show that careful identification and control of experimental error allows for highly accurate gene expression measurements. We generated antibodies in different host species, allowing for simultaneous staining of four Drosophila gap genes in individual embryos. Careful error analysis of hundreds of expression profiles reveals that less than ∼20% of the observed embryo-to-embryo fluctuations stem from experimental error. These measurements make it possible to extract not only very accurate mean gene expression profiles but also their naturally occurring fluctuations of biological origin and corresponding cross-correlations. We use this analysis to extract gap gene profile dynamics with ∼1 min accuracy. The combination of these new measurements and analysis techniques reveals a twofold increase in profile reproducibility owing to a collective network dynamics that relays positional accuracy from the maternal gradients to the pair-rule genes. PMID:23340845
Fang, Tao; Li, Wei; Gu, Fangwei; Li, Shuhua
2015-01-13
We extend the generalized energy-based fragmentation (GEBF) approach to molecular crystals under periodic boundary conditions (PBC), and we demonstrate the performance of the method for a variety of molecular crystals. With this approach, the lattice energy of a molecular crystal can be obtained from the energies of a series of embedded subsystems, which can be computed with existing advanced molecular quantum chemistry methods. The use of the field compensation method allows the method to take long-range electrostatic interaction of the infinite crystal environment into account and make the method almost translationally invariant. The computational cost of the present method scales linearly with the number of molecules in the unit cell. Illustrative applications demonstrate that the PBC-GEBF method with explicitly correlated quantum chemistry methods is capable of providing accurate descriptions on the lattice energies and structures for various types of molecular crystals. In addition, this approach can be employed to quantify the contributions of various intermolecular interactions to the theoretical lattice energy. Such qualitative understanding is very useful for rational design of molecular crystals. PMID:26574207
Grebner, Christoph; Becker, Johannes; Weber, Daniel; Bellinger, Daniel; Tafipolski, Maxim; Brückner, Charlotte; Engels, Bernd
2014-09-15
The presented program package, Conformational Analysis and Search Tool (CAST) allows the accurate treatment of large and flexible (macro) molecular systems. For the determination of thermally accessible minima CAST offers the newly developed TabuSearch algorithm, but algorithms such as Monte Carlo (MC), MC with minimization, and molecular dynamics are implemented as well. For the determination of reaction paths, CAST provides the PathOpt, the Nudge Elastic band, and the umbrella sampling approach. Access to free energies is possible through the free energy perturbation approach. Along with a number of standard force fields, a newly developed symmetry-adapted perturbation theory-based force field is included. Semiempirical computations are possible through DFTB+ and MOPAC interfaces. For calculations based on density functional theory, a Message Passing Interface (MPI) interface to the Graphics Processing Unit (GPU)-accelerated TeraChem program is available. The program is available on request. PMID:25056524
Molecular photoionization dynamics
Dehmer, Joseph L.
1982-05-01
This program seeks to develop both physical insight and quantitative characterization of molecular photoionization processes. Progress is briefly described, and some publications resulting from the research are listed. (WHK)
Langevin stabilization of molecular dynamics
NASA Astrophysics Data System (ADS)
Izaguirre, Jesús A.; Catarello, Daniel P.; Wozniak, Justin M.; Skeel, Robert D.
2001-02-01
In this paper we show the possibility of using very mild stochastic damping to stabilize long time step integrators for Newtonian molecular dynamics. More specifically, stable and accurate integrations are obtained for damping coefficients that are only a few percent of the natural decay rate of processes of interest, such as the velocity autocorrelation function. Two new multiple time stepping integrators, Langevin Molly (LM) and Brünger-Brooks-Karplus-Molly (BBK-M), are introduced in this paper. Both use the mollified impulse method for the Newtonian term. LM uses a discretization of the Langevin equation that is exact for the constant force, and BBK-M uses the popular Brünger-Brooks-Karplus integrator (BBK). These integrators, along with an extrapolative method called LN, are evaluated across a wide range of damping coefficient values. When large damping coefficients are used, as one would for the implicit modeling of solvent molecules, the method LN is superior, with LM closely following. However, with mild damping of 0.2 ps-1, LM produces the best results, allowing long time steps of 14 fs in simulations containing explicitly modeled flexible water. With BBK-M and the same damping coefficient, time steps of 12 fs are possible for the same system. Similar results are obtained for a solvated protein-DNA simulation of estrogen receptor ER with estrogen response element ERE. A parallel version of BBK-M runs nearly three times faster than the Verlet-I/r-RESPA (reversible reference system propagator algorithm) when using the largest stable time step on each one, and it also parallelizes well. The computation of diffusion coefficients for flexible water and ER/ERE shows that when mild damping of up to 0.2 ps-1 is used the dynamics are not significantly distorted.
Spectroscopy and molecular dynamics in nonpolar fluids
NASA Astrophysics Data System (ADS)
Everitt, Karl Frederick
This thesis considers the mechanisms by which molecular dynamics in nonpolar liquids influences solvation dynamics and vibrational energy relaxation. We use semiclassical molecular dynamics simulations to calculate photon echo signals for two simple fluids. We demonstrate that two new observables are directly related to the relevant molecular quantity, the frequency- frequency time correlation function (TCF), in contrast to the commonly measured 3PEPS, which cannot be simply related to this TCF at short times. We also present a semianalytic photon echo theory, based on an ansatz which determines the full time dependence from the short time expansion coefficients of the TCF. We demonstrate that this theory accurately predicts most photon echo observables, even when the theory's gaussian approximation is not accurate. We also consider vibrational energy relaxation (VER) in liquid oxygen. Using semiclassical molecular dynamics simulations and an intermolecular potential from the literature, we evaluate the required quantity (the spectral density of a certain force-force TCF) using the same ansatz described above. We demonstrate numerically that this procedure is accurate. Approximately relating this semiclassical rate to the fully quantum mechanical VER rate, using one of the more accurate ``quantum corrections'' available in the literature, yields a result which is in order-of-magnitude agreement with the experimental VER rate. We also calculate the VER rate for liquid oxygen/argon mixtures. The rotations of the solvent near a vibrationally excited molecule, and of that molecule itself, have important consequences for the short-time dynamics of the force-force TCF. We propose a simple statistical model which quantitatively explains the mole- fraction dependence of the observed VER rate. Next, we demonstrate that a newly-developed model for oxygen very accurately describes the liquid, by comparing to experimental measures of microscopic structure and dynamics. We also
Dynamic-domain-decomposition parallel molecular dynamics
NASA Astrophysics Data System (ADS)
Srinivasan, S. G.; Ashok, I.; Jônsson, Hannes; Kalonji, Gretchen; Zahorjan, John
1997-05-01
Parallel molecular dynamics with short-range forces can suffer from load-imbalance problems and attendant performance degradation due to density variations in the simulated system. In this paper, we describe an approach to dynamical load balancing, enabled by the Ādhāra runtime system. The domain assigned to each processor is automatically and dynamically resized so as to evenly distribute the molecular dynamics computations across all the processors. The algorithm was tested on an Intel Paragon parallel computer for two and three-dimensional Lennard-Jones systems containing 99 458 and 256000 atoms, respectively, and using up to 256 processors. In these benchmarks, the overhead for carrying out the load-balancing operations was found to be small and the total computation time was reduced by as much as 50%.
Romanov, V N; Cygan, R T; Myshakin, E M
2012-06-21
Naturally occurring clay minerals provide a distinctive material for carbon capture and carbon dioxide sequestration. Swelling clay minerals, such as the smectite variety, possess an aluminosilicate structure that is controlled by low-charge layers that readily expand to accommodate water molecules and, potentially, CO2. Recent experimental studies have demonstrated the efficacy of intercalating CO2 in the interlayer of layered clays, but little is known about the molecular mechanisms of the process and the extent of carbon capture as a function of clay charge and structure. A series of molecular dynamics simulations and vibrational analyses have been completed to assess the molecular interactions associated with incorporation of CO2 and H2O in the interlayer of montmorillonite clay and to help validate the models with experimental observation. An accurate and fully flexible set of interatomic potentials for CO2 is developed and combined with Clayff potentials to help evaluate the intercalation mechanism and examine the effect of molecular flexibility onthe diffusion rate of CO2 in water.
Accurate and predictive antibody repertoire profiling by molecular amplification fingerprinting
Khan, Tarik A.; Friedensohn, Simon; de Vries, Arthur R. Gorter; Straszewski, Jakub; Ruscheweyh, Hans-Joachim; Reddy, Sai T.
2016-01-01
High-throughput antibody repertoire sequencing (Ig-seq) provides quantitative molecular information on humoral immunity. However, Ig-seq is compromised by biases and errors introduced during library preparation and sequencing. By using synthetic antibody spike-in genes, we determined that primer bias from multiplex polymerase chain reaction (PCR) library preparation resulted in antibody frequencies with only 42 to 62% accuracy. Additionally, Ig-seq errors resulted in antibody diversity measurements being overestimated by up to 5000-fold. To rectify this, we developed molecular amplification fingerprinting (MAF), which uses unique molecular identifier (UID) tagging before and during multiplex PCR amplification, which enabled tagging of transcripts while accounting for PCR efficiency. Combined with a bioinformatic pipeline, MAF bias correction led to measurements of antibody frequencies with up to 99% accuracy. We also used MAF to correct PCR and sequencing errors, resulting in enhanced accuracy of full-length antibody diversity measurements, achieving 98 to 100% error correction. Using murine MAF-corrected data, we established a quantitative metric of recent clonal expansion—the intraclonal diversity index—which measures the number of unique transcripts associated with an antibody clone. We used this intraclonal diversity index along with antibody frequencies and somatic hypermutation to build a logistic regression model for prediction of the immunological status of clones. The model was able to predict clonal status with high confidence but only when using MAF error and bias corrected Ig-seq data. Improved accuracy by MAF provides the potential to greatly advance Ig-seq and its utility in immunology and biotechnology. PMID:26998518
Accurate and predictive antibody repertoire profiling by molecular amplification fingerprinting.
Khan, Tarik A; Friedensohn, Simon; Gorter de Vries, Arthur R; Straszewski, Jakub; Ruscheweyh, Hans-Joachim; Reddy, Sai T
2016-03-01
High-throughput antibody repertoire sequencing (Ig-seq) provides quantitative molecular information on humoral immunity. However, Ig-seq is compromised by biases and errors introduced during library preparation and sequencing. By using synthetic antibody spike-in genes, we determined that primer bias from multiplex polymerase chain reaction (PCR) library preparation resulted in antibody frequencies with only 42 to 62% accuracy. Additionally, Ig-seq errors resulted in antibody diversity measurements being overestimated by up to 5000-fold. To rectify this, we developed molecular amplification fingerprinting (MAF), which uses unique molecular identifier (UID) tagging before and during multiplex PCR amplification, which enabled tagging of transcripts while accounting for PCR efficiency. Combined with a bioinformatic pipeline, MAF bias correction led to measurements of antibody frequencies with up to 99% accuracy. We also used MAF to correct PCR and sequencing errors, resulting in enhanced accuracy of full-length antibody diversity measurements, achieving 98 to 100% error correction. Using murine MAF-corrected data, we established a quantitative metric of recent clonal expansion-the intraclonal diversity index-which measures the number of unique transcripts associated with an antibody clone. We used this intraclonal diversity index along with antibody frequencies and somatic hypermutation to build a logistic regression model for prediction of the immunological status of clones. The model was able to predict clonal status with high confidence but only when using MAF error and bias corrected Ig-seq data. Improved accuracy by MAF provides the potential to greatly advance Ig-seq and its utility in immunology and biotechnology. PMID:26998518
Accurate molecular classification of cancer using simple rules
Wang, Xiaosheng; Gotoh, Osamu
2009-01-01
Background One intractable problem with using microarray data analysis for cancer classification is how to reduce the extremely high-dimensionality gene feature data to remove the effects of noise. Feature selection is often used to address this problem by selecting informative genes from among thousands or tens of thousands of genes. However, most of the existing methods of microarray-based cancer classification utilize too many genes to achieve accurate classification, which often hampers the interpretability of the models. For a better understanding of the classification results, it is desirable to develop simpler rule-based models with as few marker genes as possible. Methods We screened a small number of informative single genes and gene pairs on the basis of their depended degrees proposed in rough sets. Applying the decision rules induced by the selected genes or gene pairs, we constructed cancer classifiers. We tested the efficacy of the classifiers by leave-one-out cross-validation (LOOCV) of training sets and classification of independent test sets. Results We applied our methods to five cancerous gene expression datasets: leukemia (acute lymphoblastic leukemia [ALL] vs. acute myeloid leukemia [AML]), lung cancer, prostate cancer, breast cancer, and leukemia (ALL vs. mixed-lineage leukemia [MLL] vs. AML). Accurate classification outcomes were obtained by utilizing just one or two genes. Some genes that correlated closely with the pathogenesis of relevant cancers were identified. In terms of both classification performance and algorithm simplicity, our approach outperformed or at least matched existing methods. Conclusion In cancerous gene expression datasets, a small number of genes, even one or two if selected correctly, is capable of achieving an ideal cancer classification effect. This finding also means that very simple rules may perform well for cancerous class prediction. PMID:19874631
Molecular dynamics and protein function
Karplus, M.; Kuriyan, J.
2005-01-01
A fundamental appreciation for how biological macromolecules work requires knowledge of structure and dynamics. Molecular dynamics simulations provide powerful tools for the exploration of the conformational energy landscape accessible to these molecules, and the rapid increase in computational power coupled with improvements in methodology makes this an exciting time for the application of simulation to structural biology. In this Perspective we survey two areas, protein folding and enzymatic catalysis, in which simulations have contributed to a general understanding of mechanism. We also describe results for the F1 ATPase molecular motor and the Src family of signaling proteins as examples of applications of simulations to specific biological systems. PMID:15870208
Accurate Langevin approaches to simulate Markovian channel dynamics
NASA Astrophysics Data System (ADS)
Huang, Yandong; Rüdiger, Sten; Shuai, Jianwei
2015-12-01
The stochasticity of ion-channels dynamic is significant for physiological processes on neuronal cell membranes. Microscopic simulations of the ion-channel gating with Markov chains can be considered to be an accurate standard. However, such Markovian simulations are computationally demanding for membrane areas of physiologically relevant sizes, which makes the noise-approximating or Langevin equation methods advantageous in many cases. In this review, we discuss the Langevin-like approaches, including the channel-based and simplified subunit-based stochastic differential equations proposed by Fox and Lu, and the effective Langevin approaches in which colored noise is added to deterministic differential equations. In the framework of Fox and Lu’s classical models, several variants of numerical algorithms, which have been recently developed to improve accuracy as well as efficiency, are also discussed. Through the comparison of different simulation algorithms of ion-channel noise with the standard Markovian simulation, we aim to reveal the extent to which the existing Langevin-like methods approximate results using Markovian methods. Open questions for future studies are also discussed.
Hall, G.E.; Prrese, J.M.; Sears, T.J.; Weston, R.E.
1999-05-21
The goal of this research is the understanding of elementary chemical and physical processes important in the combustion of fossil fuels. Interest centers on reactions involving short-lived chemical intermediates and their properties. High-resolution high-sensitivity laser absorption methods are augmented by high temperature flow-tube reaction kinetics studies with mass spectrometric sampling. These experiments provide information on the energy levels, structures and reactivity of molecular flee radical species and, in turn, provide new tools for the study of energy flow and chemical bond cleavage in the radicals in chemical systems. The experimental work is supported by theoretical and computational work using time-dependent quantum wavepacket calculations that provide insights into energy flow between the vibrational modes of the molecule.
ADAPTIVE MULTILEVEL SPLITTING IN MOLECULAR DYNAMICS SIMULATIONS*
Aristoff, David; Lelièvre, Tony; Mayne, Christopher G.; Teo, Ivan
2014-01-01
Adaptive Multilevel Splitting (AMS) is a replica-based rare event sampling method that has been used successfully in high-dimensional stochastic simulations to identify trajectories across a high potential barrier separating one metastable state from another, and to estimate the probability of observing such a trajectory. An attractive feature of AMS is that, in the limit of a large number of replicas, it remains valid regardless of the choice of reaction coordinate used to characterize the trajectories. Previous studies have shown AMS to be accurate in Monte Carlo simulations. In this study, we extend the application of AMS to molecular dynamics simulations and demonstrate its effectiveness using a simple test system. Our conclusion paves the way for useful applications, such as molecular dynamics calculations of the characteristic time of drug dissociation from a protein target. PMID:26005670
Molecular dynamics investigation of nanoscale cavitation dynamics
NASA Astrophysics Data System (ADS)
Sasikumar, Kiran; Keblinski, Pawel
2014-12-01
We use molecular dynamics simulations to investigate the cavitation dynamics around intensely heated solid nanoparticles immersed in a model Lennard-Jones fluid. Specifically, we study the temporal evolution of vapor nanobubbles that form around the solid nanoparticles heated over ps time scale and provide a detail description of the following vapor formation and collapse. For 8 nm diameter nanoparticles we observe the formation of vapor bubbles when the liquid temperature 0.5-1 nm away from the nanoparticle surface reaches ˜90% of the critical temperature, which is consistent with the onset of spinodal decomposition. The peak heat flux from the hot solid to the surrounding liquid at the bubble formation threshold is ˜20 times higher than the corresponding steady state critical heat flux. Detailed analysis of the bubble dynamics indicates adiabatic formation followed by an isothermal final stage of growth and isothermal collapse.
Dynamic fracture toughness determined using molecular dynamics
Swadener, J. G.; Baskes, M. I.; Nastasi, Michael Anthony,
2004-01-01
Molecular dynamics (MD) simulations of fracture in crystalline silicon are conducted in order to determine the dynamic fracture toughness. The MD simulations show how the potential energy released during fracture is partitioned into surface energy, energy stored in defects and kinetic energy. First, the MD fracture simulations are shown to produce brittle fracture and be in reasonable agreement with experimental results. Then dynamic hcture toughness is calculated as the sum of the surface energy and the energy stored as defects directly from the MD models. Models oriented to produce fracture on either (111) or (101) planes are used. For the (101) fracture orientation, equilibrium crack speeds of greater than 80% of the Rayleigh wave speed are obtained. Crack speeds initially show a steep increase with increasing energy release rate followed by a much more gradual increase. No plateau in crack speed is observed for static energy release rates up to 20 J/m{sup 2}. At the point where the change in crack speed behavior occur, the dynamic fracture toughness (J{sub d}) is still within 10% of two times the surface energy (2{gamma}{sub 0}) and changing very slowly. From these MD simulations, it appears that the change in crack speed behavior is due to a change in the kinetic energy generation during dynamic fracture. In addition, MD simulations of facture in silicon with defects were conducted. The addition of defects increases the inelastic dissipation and the energy stored in defects.
Fiftieth anniversary of molecular dynamics
NASA Astrophysics Data System (ADS)
Melker, Alexander I.
2007-04-01
The history of computer application in physics for solving nonlinear problems is considered. Examples from different branches of condensed matter physics (nonlinear vibrations of anharmonic chains of atoms, dynamics of radiation damage of crystals, deformation and fracture of crystals) are given. A new line of investigation and the results obtained in the field of computer simulation of physical processes realized in the department of metal physics and computer technologies in materials science are considered. This line incorporates both a study of self-organization and properties of new materials (fullerenes, carbon nanotubes) and biological objects by molecular dynamics technique as well as the development of new computer simulation methods.
Available Instruments for Analyzing Molecular Dynamics Trajectories
Likhachev, I. V.; Balabaev, N. K.; Galzitskaya, O. V.
2016-01-01
Molecular dynamics trajectories are the result of molecular dynamics simulations. Trajectories are sequential snapshots of simulated molecular system which represents atomic coordinates at specific time periods. Based on the definition, in a text format trajectory files are characterized by their simplicity and uselessness. To obtain information from such files, special programs and information processing techniques are applied: from molecular dynamics animation to finding characteristics along the trajectory (versus time). In this review, we describe different programs for processing molecular dynamics trajectories. The performance of these programs, usefulness for analyses of molecular dynamics trajectories, strong and weak aspects are discussed. PMID:27053964
Molecular dynamics on vector computers
NASA Astrophysics Data System (ADS)
Sullivan, F.; Mountain, R. D.; Oconnell, J.
1985-10-01
An algorithm called the method of lights (MOL) has been developed for the computerized simulation of molecular dynamics. The MOL, implemented on the CYBER 205 computer, is based on sorting and reformulating the manner in which neighbor lists are compiled, and it uses data structures compatible with specialized vector statements that perform parallel computations. The MOL is found to reduce running time over standard methods in scalar form, and vectorization is shown to produce an order-of-magnitude reduction in execution time.
Molecular dynamics for dense matter
NASA Astrophysics Data System (ADS)
Maruyama, Toshiki; Watanabe, Gentaro; Chiba, Satoshi
2012-08-01
We review a molecular dynamics method for nucleon many-body systems called quantum molecular dynamics (QMD), and our studies using this method. These studies address the structure and the dynamics of nuclear matter relevant to neutron star crusts, supernova cores, and heavy-ion collisions. A key advantage of QMD is that we can study dynamical processes of nucleon many-body systems without any assumptions about the nuclear structure. First, we focus on the inhomogeneous structures of low-density nuclear matter consisting not only of spherical nuclei but also of nuclear "pasta", i.e., rod-like and slab-like nuclei. We show that pasta phases can appear in the ground and equilibrium states of nuclear matter without assuming nuclear shape. Next, we show our simulation of compression of nuclear matter which corresponds to the collapsing stage of supernovae. With the increase in density, a crystalline solid of spherical nuclei changes to a triangular lattice of rods by connecting neighboring nuclei. Finally, we discuss fragment formation in expanding nuclear matter. Our results suggest that a generally accepted scenario based on the liquid-gas phase transition is not plausible at lower temperatures.
Numerical methods for molecular dynamics
Skeel, R.D.
1991-01-01
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
Scalable Molecular Dynamics with NAMD
Phillips, James C.; Braun, Rosemary; Wang, Wei; Gumbart, James; Tajkhorshid, Emad; Villa, Elizabeth; Chipot, Christophe; Skeel, Robert D.; Kalé, Laxmikant; Schulten, Klaus
2008-01-01
NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This paper, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Next, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, e.g., the Tcl scripting language. Finally, the paper provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu. PMID:16222654
Local Refinements in Classical Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Fackeldey, Konstantin; Weber, Marcus
2014-03-01
Quantum mechanics provide a detailed description of the physical and chemical behavior of molecules. However, with increasing size of the system the complexity rises exponentially, which is prohibitive for efficient dynamical simulation. In contrast, classical molecular dynamics procure a coarser description by using less degrees of freedom. Thus, it seems natural to seek for an adequate trade-off between accurateness and computational feasibility in the simulation of molecules. Here, we propose a novel method, which combines classical molecular simulations with quantum mechanics for molecular systems. For this we decompose the state space of the respective molecule into subsets, by employing a meshfree partition of unity. We show, that this partition allows us to localize an empirical force field and to run locally constrained classical trajectories. Within each subset, we compute the energy on the quantum level for a fixed number of spatial states (ab initio points). With these energy values from the ab initio points we have a local scattered data problem, which can be solved by the moving least squares method.
Better, Cheaper, Faster Molecular Dynamics
NASA Technical Reports Server (NTRS)
Pohorille, Andrew; DeVincenzi, Donald L. (Technical Monitor)
2001-01-01
Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.
Molecular Dynamics Simulations of Polymers
NASA Astrophysics Data System (ADS)
Han, Jie
1995-01-01
Molecular dynamics (MD) simulations have been undertaken in this work to explore structures and properties of polyethylene (PE), polyisobutylene (PIB), atactic polypropylene (aPP) and atactic polystyrene (aPS). This work has not only demonstrated the reliability of MD simulations by comparing results with available experiments, but more importantly has revealed structure-property relationships on a molecular level for these selected polymers. Structures of these amorphous polymers were characterized by radial distribution functions (RDFs) or scattering profiles, and properties of the polymers studied were pressure-volume -temperature (PVT) equation of state, enthalpy, cohesive energy, the diffusion coefficient of methane in the polymer, and glass transition temperature. Good agreement was found for these structures and properties between simulation and experiment. More importantly, the scientific understanding of structure-property relationships was established on a molecular level. In the order of aPP (PE), PIB and aPS, with the chain surface separation or free volume decreasing, the density increases and the diffusion coefficient decreases. Therefore, the effects of changes or modifications in the chemical structure of monomer molecules (substituting pendent hydrogen with methyl or phenyl) on polymeric materials performance were attributed to the effects of molecular chain structure on packing structure, which, in turn, affects the properties of these polymers. Local chain dynamics and relaxation have been studied for bulk PE and aPS. Cooperative transitions occur at second-neighbor bonds for PE, and first-neighbor bonds for aPS due to the role of side groups. The activation energy is a single torsional barrier for overall conformational transitions, and is single torsional barrier plus locally "trapped" barrier for relaxation. Temperature dependence is Arrhenius for transition time, and is WLF for relaxation time. The mean correlation times derived from
Accurate Determination of Membrane Dynamics with Line-Scan FCS
Ries, Jonas; Chiantia, Salvatore; Schwille, Petra
2009-01-01
Here we present an efficient implementation of line-scan fluorescence correlation spectroscopy (i.e., one-dimensional spatio-temporal image correlation spectroscopy) using a commercial laser scanning microscope, which allows the accurate measurement of diffusion coefficients and concentrations in biological lipid membranes within seconds. Line-scan fluorescence correlation spectroscopy is a calibration-free technique. Therefore, it is insensitive to optical artifacts, saturation, or incorrect positioning of the laser focus. In addition, it is virtually unaffected by photobleaching. Correction schemes for residual inhomogeneities and depletion of fluorophores due to photobleaching extend the applicability of line-scan fluorescence correlation spectroscopy to more demanding systems. This technique enabled us to measure accurate diffusion coefficients and partition coefficients of fluorescent lipids in phase-separating supported bilayers of three commonly used raft-mimicking compositions. Furthermore, we probed the temperature dependence of the diffusion coefficient in several model membranes, and in human embryonic kidney cell membranes not affected by temperature-induced optical aberrations. PMID:19254560
An Accurate and Dynamic Computer Graphics Muscle Model
NASA Technical Reports Server (NTRS)
Levine, David Asher
1997-01-01
A computer based musculo-skeletal model was developed at the University in the departments of Mechanical and Biomedical Engineering. This model accurately represents human shoulder kinematics. The result of this model is the graphical display of bones moving through an appropriate range of motion based on inputs of EMGs and external forces. The need existed to incorporate a geometric muscle model in the larger musculo-skeletal model. Previous muscle models did not accurately represent muscle geometries, nor did they account for the kinematics of tendons. This thesis covers the creation of a new muscle model for use in the above musculo-skeletal model. This muscle model was based on anatomical data from the Visible Human Project (VHP) cadaver study. Two-dimensional digital images from the VHP were analyzed and reconstructed to recreate the three-dimensional muscle geometries. The recreated geometries were smoothed, reduced, and sliced to form data files defining the surfaces of each muscle. The muscle modeling function opened these files during run-time and recreated the muscle surface. The modeling function applied constant volume limitations to the muscle and constant geometry limitations to the tendons.
Accurate direct Eulerian simulation of dynamic elastic-plastic flow
Kamm, James R; Walter, John W
2009-01-01
The simulation of dynamic, large strain deformation is an important, difficult, and unsolved computational challenge. Existing Eulerian schemes for dynamic material response are plagued by unresolved issues. We present a new scheme for the first-order system of elasto-plasticity equations in the Eulerian frame. This system has an intrinsic constraint on the inverse deformation gradient. Standard Godunov schemes do not satisfy this constraint. The method of Flux Distributions (FD) was devised to discretely enforce such constraints for numerical schemes with cell-centered variables. We describe a Flux Distribution approach that enforces the inverse deformation gradient constraint. As this approach is new and novel, we do not yet have numerical results to validate our claims. This paper is the first installment of our program to develop this new method.
Radiation in molecular dynamic simulations
Glosli, J; Graziani, F; More, R; Murillo, M; Streitz, F; Surh, M
2008-10-13
Hot dense radiative (HDR) plasmas common to Inertial Confinement Fusion (ICF) and stellar interiors have high temperature (a few hundred eV to tens of keV), high density (tens to hundreds of g/cc) and high pressure (hundreds of Megabars to thousands of Gigabars). Typically, such plasmas undergo collisional, radiative, atomic and possibly thermonuclear processes. In order to describe HDR plasmas, computational physicists in ICF and astrophysics use atomic-scale microphysical models implemented in various simulation codes. Experimental validation of the models used to describe HDR plasmas are difficult to perform. Direct Numerical Simulation (DNS) of the many-body interactions of plasmas is a promising approach to model validation but, previous work either relies on the collisionless approximation or ignores radiation. We present a new numerical simulation technique to address a currently unsolved problem: the extension of molecular dynamics to collisional plasmas including emission and absorption of radiation. The new technique passes a key test: it relaxes to a blackbody spectrum for a plasma in local thermodynamic equilibrium. This new tool also provides a method for assessing the accuracy of energy and momentum exchange models in hot dense plasmas. As an example, we simulate the evolution of non-equilibrium electron, ion, and radiation temperatures for a hydrogen plasma using the new molecular dynamics simulation capability.
The "Collisions Cube" Molecular Dynamics Simulator.
ERIC Educational Resources Information Center
Nash, John J.; Smith, Paul E.
1995-01-01
Describes a molecular dynamics simulator that employs ping-pong balls as the atoms or molecules and is suitable for either large lecture halls or small classrooms. Discusses its use in illustrating many of the fundamental concepts related to molecular motion and dynamics and providing a three-dimensional perspective of molecular motion. (JRH)
Accurate boundary conditions for exterior problems in gas dynamics
NASA Technical Reports Server (NTRS)
Hagstrom, Thomas; Hariharan, S. I.
1988-01-01
The numerical solution of exterior problems is typically accomplished by introducing an artificial, far field boundary and solving the equations on a truncated domain. For hyperbolic systems, boundary conditions at this boundary are often derived by imposing a principle of no reflection. However, waves with spherical symmetry in gas dynamics satisfy equations where incoming and outgoing Riemann variables are coupled. This suggests that natural reflections may be important. A reflecting boundary condition is proposed based on an asymptotic solution of the far field equations. Nonlinear energy estimates are obtained for the truncated problem and numerical experiments presented to validate the theory.
Accurate boundary conditions for exterior problems in gas dynamics
NASA Technical Reports Server (NTRS)
Hagstrom, Thomas; Hariharan, S. I.
1988-01-01
The numerical solution of exterior problems is typically accomplished by introducing an artificial, far-field boundary and solving the equations on a truncated domain. For hyperbolic systems, boundary conditions at this boundary are often derived by imposing a principle of no reflection. However, waves with spherical symmetry in gas dynamics satisfy equations where incoming and outgoing Riemann variables are coupled. This suggests that natural reflections may be important. A reflecting boundary condition is proposed based on an asymptotic solution of the far-field equations. Nonlinear energy estimates are obtained for the truncated problem and numerical experiments presented to validate the theory.
Dynamic pseudos: How accurate outside their parent case?
Ekrann, S.; Mykkeltveit, J.
1995-12-31
If properly constructed, dynamic pseudos allow the parent solution from which they were derived to be exactly reproduced, in a certain well-defined sense, in a subsequent coarse grid simulation. The paper reports extensive numerical experimentation, in 1D homogeneous and heterogeneous media, to determine the performance of pseudos when used outside their parent case. The authors perturb fluid viscosities and injection rate, as well as realization. Parent solutions are produced analytically, via a generalization of the Buckley-Leverett technique, as are true solutions in off-parent cases. Capillarity is neglected in these experiments, while gravity is sometimes retained in order to force rate sensitivity.
A concurrent multiscale micromorphic molecular dynamics
Li, Shaofan Tong, Qi
2015-04-21
In this work, we have derived a multiscale micromorphic molecular dynamics (MMMD) from first principle to extend the (Andersen)-Parrinello-Rahman molecular dynamics to mesoscale and continuum scale. The multiscale micromorphic molecular dynamics is a con-current three-scale dynamics that couples a fine scale molecular dynamics, a mesoscale micromorphic dynamics, and a macroscale nonlocal particle dynamics together. By choosing proper statistical closure conditions, we have shown that the original Andersen-Parrinello-Rahman molecular dynamics is the homogeneous and equilibrium case of the proposed multiscale micromorphic molecular dynamics. In specific, we have shown that the Andersen-Parrinello-Rahman molecular dynamics can be rigorously formulated and justified from first principle, and its general inhomogeneous case, i.e., the three scale con-current multiscale micromorphic molecular dynamics can take into account of macroscale continuum mechanics boundary condition without the limitation of atomistic boundary condition or periodic boundary conditions. The discovered multiscale scale structure and the corresponding multiscale dynamics reveal a seamless transition from atomistic scale to continuum scale and the intrinsic coupling mechanism among them based on first principle formulation.
Nonadiabatic Molecular Dynamics with Trajectories
NASA Astrophysics Data System (ADS)
Tavernelli, Ivano
2012-02-01
In the mixed quantum-classical description of molecular systems, only the quantum character of the electronic degrees of freedom is considered while the nuclear motion is treated at a classical level. In the adiabatic case, this picture corresponds to the Born-Oppenheimer limit where the nuclei move as point charges on the potential energy surface (PES) associated with a given electronic state. Despite the success of this approximation, many physical and chemical processes do not fall in the regime where nuclei and electrons can be considered decoupled. In particular, most photoreactions pass through regions of the PES in which electron-nuclear quantum interference effects are sizeable and often crucial for a correct description of the phenomena. Recently, we have developed a trajectory-based nonadiabatic molecular dynamics scheme that describes the nuclear wavepacket as an ensemble of particles following classical trajectories on PESs derived from time-dependent density functional theory (TDDFT) [1]. The method is based on Tully's fewest switches trajectories surface hopping (TSH) where the nonadiabatic coupling elements between the different potential energy surfaces are computed on-the-fly as functionals of the ground state electron density or, equivalently, of the corresponding Kohn-Sham orbitals [2]. Here, we present the theoretical fundamentals of our approach together with an extension that allows for the direct coupling of the dynamics to an external electromagnetic field [3] as well as to the external potential generated by the environment (solvent effects) [4]. The method is applied to the study of the photodissociation dynamics of simple molecules in gas phase and to the description of the fast excited state dynamics of molecules in solution (in particular Ruthenium (II) tris(bipyridine) in water). [4pt] [1] E. Tapavicza, I. Tavernelli, U. Rothlisberger, Phys. Rev. Lett., 98, (2007) 023001. [0pt] [2] Tavernelli I.; Tapavicza E.; Rothlisberger U., J. Chem
Molecular dynamics simulation of benzene
NASA Astrophysics Data System (ADS)
Trumpakaj, Zygmunt; Linde, Bogumił B. J.
2016-03-01
Intermolecular potentials and a few models of intermolecular interaction in liquid benzene are tested by Molecular Dynamics (MD) simulations. The repulsive part of the Lennard-Jones 12-6 (LJ 12-6) potential is too hard, which yields incorrect results. The exp-6 potential with a too hard repulsive term is also often used. Therefore, we took an expa-6 potential with a small Gaussian correction plus electrostatic interactions. This allows to modify the curvature of the potential. The MD simulations are carried out in the temperature range 280-352 K under normal pressure and at experimental density. The Rayleigh scattering of depolarized light is used for comparison. The results of MD simulations are comparable with the experimental values.
Buckybomb: Reactive Molecular Dynamics Simulation.
Chaban, Vitaly V; Fileti, Eudes Eterno; Prezhdo, Oleg V
2015-03-01
Energetic materials, such as explosives, propellants, and pyrotechnics, are widely used in civilian and military applications. Nanoscale explosives represent a special group because of the high density of energetic covalent bonds. The reactive molecular dynamics (ReaxFF) study of nitrofullerene decomposition reported here provides a detailed chemical mechanism of explosion of a nanoscale carbon material. Upon initial heating, C60(NO2)12 disintegrates, increasing temperature and pressure by thousands of Kelvins and bars within tens of picoseconds. The explosion starts with NO2 group isomerization into C-O-N-O, followed by emission of NO molecules and formation of CO groups on the buckyball surface. NO oxidizes into NO2, and C60 falls apart, liberating CO2. At the highest temperatures, CO2 gives rise to diatomic carbon. The study shows that the initiation temperature and released energy depend strongly on the chemical composition and density of the material. PMID:26262672
Molecular dynamics of membrane proteins.
Woolf, Thomas B.; Crozier, Paul Stewart; Stevens, Mark Jackson
2004-10-01
Understanding the dynamics of the membrane protein rhodopsin will have broad implications for other membrane proteins and cellular signaling processes. Rhodopsin (Rho) is a light activated G-protein coupled receptor (GPCR). When activated by ligands, GPCRs bind and activate G-proteins residing within the cell and begin a signaling cascade that results in the cell's response to external stimuli. More than 50% of all current drugs are targeted toward G-proteins. Rho is the prototypical member of the class A GPCR superfamily. Understanding the activation of Rho and its interaction with its Gprotein can therefore lead to a wider understanding of the mechanisms of GPCR activation and G-protein activation. Understanding the dark to light transition of Rho is fully analogous to the general ligand binding and activation problem for GPCRs. This transition is dependent on the lipid environment. The effect of lipids on membrane protein activity in general has had little attention, but evidence is beginning to show a significant role for lipids in membrane protein activity. Using the LAMMPS program and simulation methods benchmarked under the IBIG program, we perform a variety of allatom molecular dynamics simulations of membrane proteins.
Sezen, Melda; Register, Jeffrey T; Yao, Yao; Glisic, Branko; Loo, Yueh-Lin
2016-06-01
The polarity and the magnitude of polyaniline's gauge factor are tuned through structural modification. Combining conducting polymers with gauge factors of opposite polarities yields an accurate temperature sensor, even when deployed under dynamic strains. PMID:27061270
Zhou, Nengji; Chen, Lipeng; Huang, Zhongkai; Sun, Kewei; Tanimura, Yoshitaka; Zhao, Yang
2016-03-10
By employing the Dirac-Frenkel time-dependent variational principle, we study the dynamical properties of the Holstein molecular crystal model with diagonal and off-diagonal exciton-phonon coupling. A linear combination of the Davydov D1 (D2) ansatz, referred to as the "multi-D1 ansatz" ("multi-D2 ansatz"), is used as the trial state with enhanced accuracy but without sacrificing efficiency. The time evolution of the exciton probability is found to be in perfect agreement with that of the hierarchy equations of motion, demonstrating the promise the multiple Davydov trial states hold as an efficient, robust description of dynamics of complex quantum systems. In addition to the linear absorption spectra computed for both diagonal and off-diagonal cases, for the first time, 2D spectra have been calculated for systems with off-diagonal exciton-phonon coupling by employing the multiple D2 ansatz to compute the nonlinear response function, testifying to the great potential of the multiple D2 ansatz for fast, accurate implementation of multidimensional spectroscopy. It is found that the signal exhibits a single peak for weak off-diagonal coupling, while a vibronic multipeak structure appears for strong off-diagonal coupling. PMID:26871592
Shapiro like steps reveals molecular nanomagnets' spin dynamics
NASA Astrophysics Data System (ADS)
Abdollahipour, Babak; Abouie, Jahanfar; Ebrahimi, Navid
2015-09-01
We present an accurate way to detect spin dynamics of a nutating molecular nanomagnet by inserting it in a tunnel Josephson junction and studying the current voltage (I-V) characteristic. The spin nutation of the molecular nanomagnet is generated by applying two circularly polarized magnetic fields. We demonstrate that modulation of the Josephson current by the nutation of the molecular nanomagnet's spin appears as a stepwise structure like Shapiro steps in the I-V characteristic of the junction. Width and heights of these Shapiro-like steps are determined by two parameters of the spin nutation, frequency and amplitude of the nutation, which are simply tuned by the applied magnetic fields.
Dynamical Localization in Molecular Systems.
NASA Astrophysics Data System (ADS)
Wang, Xidi
In the first four chapters of this thesis we concentrate on the Davydov model which describes the vibrational energy quanta of Amide I bonds (C=O bonds on the alpha -helix) coupled to the acoustic phonon modes of the alpha-helix backbone in the form of a Frohlich Hamiltonian. Following a brief introduction in chapter one, in chapter two we formulate the dynamics of vibrational quanta at finite temperature by using coherent state products. The fluctuation-dissipation relation is derived. At zero temperature, in the continuum limit, we recover the original results of Davydov. We also achieve good agreement with numerical simulations. In chapter three, the net contraction of the lattice is calculated exactly at any temperature, and its relation to the so -call "topological stability" of the Davydov soliton is discussed. In the second section of the chapter three we calculate the overtone spectra of crystalline acetanilide (according to some opinions ACN provides experimental evidence for the existence of Davydov solitons). Good agreement with experimental data has been obtained. In chapter four we study the self-trapped vibrational excitations by the Quantum Monte Carlo technique. For a single excitation, the temperature dependence of different physical observables is calculated. The quasi-particle which resembles the Davydov soliton has been found to be fairly narrow using the most commonly used data for the alpha -helix; at temperatures above a few Kelvin, the quasi-particle reaches its smallest limit (extends over three sites), which implies diffusive motion of the small polaron-like quasi-particle at high temperatures. For the multi-excitation case, bound pairs and clusters of excitations are found at low temperatures; they gradually dissociate when the temperature of the system is increased as calculated from the density-density correlation function. In the last chapter of this thesis, we study a more general model of dynamical local modes in molecular systems
NASA Astrophysics Data System (ADS)
Barone, Luciano Maria; Simonazzi, Riccardo; Tenenbaum, Alexander
1995-09-01
We have studied portability, efficiency and accuracy of a standard Molecular Dynamics simulation on the SIMD parallel computer APE100. Computing speed performance and physical system size range have been analyzed and compared with those of a conventional computer. Short range and long range potentials have been considered, and the comparative advantage of different simulation approaches has been assessed. For long range potentials, APE turns out to be faster than a conventional computer; large systems can be conveniently simulated using either the cloning approach (up to ˜ 10 5 particles) or a domain decomposition with the systolic method. In the case of short range potentials and systems with diffusion (like a liquid), APE is convenient only when using a large number of processors. In a special case (a crystal without diffusion), a specific domain decomposition technique with frames makes APE advantageous for intermediate and large systems. Using the latter technique we have studied in detail the effect of different numerical error sources, and compared the accuracy of APE with that of a conventional computer.
Dynamics of riboswitches: Molecular simulations.
Sanbonmatsu, Karissa Y
2014-10-01
Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches. PMID:24953187
Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations
Gottwald, Fabian; Karsten, Sven; Ivanov, Sergei D. Kühn, Oliver
2015-06-28
Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into a few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation, which can be rigorously derived by means of a linear projection technique. Within this framework, a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here, we discuss that this task is more naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.
Time-Dependent Molecular Reaction Dynamics
Oehrn, Yngve
2007-11-29
This paper is a brief review of a time-dependent, direct, nonadiabatic theory of molecular processes called Electron Nuclear Dynamics (END). This approach to the study of molecular reaction dynamics is a hierarchical theory that can be applied at various levels of approximation. The simplest level of END uses classical nuclei and represents all electrons by a single, complex, determinantal wave function. The wave function parameters such as average nuclear positions and momenta, and molecular orbital coefcients carry the time dependence and serve as dynamical variables. Examples of application are given of the simplest level of END to ion-atom and ion-molecule reactions.
Hansen, Katja; Biegler, Franziska; Ramakrishnan, Raghunathan; Pronobis, Wiktor; von Lilienfeld, O. Anatole; Müller, Klaus -Robert; Tkatchenko, Alexandre
2015-06-04
Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstratemore » prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. The same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies.« less
Hansen, Katja; Biegler, Franziska; Ramakrishnan, Raghunathan; Pronobis, Wiktor; von Lilienfeld, O. Anatole; Müller, Klaus -Robert; Tkatchenko, Alexandre
2015-06-04
Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstrate prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. The same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies.
Hele, Timothy J. H.; Willatt, Michael J.; Muolo, Andrea; Althorpe, Stuart C.
2015-05-21
We recently obtained a quantum-Boltzmann-conserving classical dynamics by making a single change to the derivation of the “Classical Wigner” approximation. Here, we show that the further approximation of this “Matsubara dynamics” gives rise to two popular heuristic methods for treating quantum Boltzmann time-correlation functions: centroid molecular dynamics (CMD) and ring-polymer molecular dynamics (RPMD). We show that CMD is a mean-field approximation to Matsubara dynamics, obtained by discarding (classical) fluctuations around the centroid, and that RPMD is the result of discarding a term in the Matsubara Liouvillian which shifts the frequencies of these fluctuations. These findings are consistent with previous numerical results and give explicit formulae for the terms that CMD and RPMD leave out.
Application of the G-JF discrete-time thermostat for fast and accurate molecular simulations
NASA Astrophysics Data System (ADS)
Grønbech-Jensen, Niels; Hayre, Natha Robert; Farago, Oded
2014-02-01
A new Langevin-Verlet thermostat that preserves the fluctuation-dissipation relationship for discrete time steps is applied to molecular modeling and tested against several popular suites (AMBER, GROMACS, LAMMPS) using a small molecule as an example that can be easily simulated by all three packages. Contrary to existing methods, the new thermostat exhibits no detectable changes in the sampling statistics as the time step is varied in the entire numerical stability range. The simple form of the method, which we express in the three common forms (Velocity-Explicit, Störmer-Verlet, and Leap-Frog), allows for easy implementation within existing molecular simulation packages to achieve faster and more accurate results with no cost in either computing time or programming complexity.
Modeling the Hydrogen Bond within Molecular Dynamics
ERIC Educational Resources Information Center
Lykos, Peter
2004-01-01
The structure of a hydrogen bond is elucidated within the framework of molecular dynamics based on the model of Rahman and Stillinger (R-S) liquid water treatment. Thus, undergraduates are exposed to the powerful but simple use of classical mechanics to solid objects from a molecular viewpoint.
Molecular Dynamics Simulations of Simple Liquids
ERIC Educational Resources Information Center
Speer, Owner F.; Wengerter, Brian C.; Taylor, Ramona S.
2004-01-01
An experiment, in which students were given the opportunity to perform molecular dynamics simulations on a series of molecular liquids using the Amber suite of programs, is presented. They were introduced to both physical theories underlying classical mechanics simulations and to the atom-atom pair distribution function.
Molecular ions, Rydberg spectroscopy and dynamics
Jungen, Ch.
2015-01-22
Ion spectroscopy, Rydberg spectroscopy and molecular dynamics are closely related subjects. Multichannel quantum defect theory is a theoretical approach which draws on this close relationship and thereby becomes a powerful tool for the study of systems consisting of a positively charged molecular ion core interacting with an electron which may be loosely bound or freely scattering.
Shaughnessy, M C; Jones, R E
2016-02-01
We develop and demonstrate a method to efficiently use density functional calculations to drive classical dynamics of complex atomic and molecular systems. The method has the potential to scale to systems and time scales unreachable with current ab initio molecular dynamics schemes. It relies on an adapting dataset of independently computed Hellmann-Feynman forces for atomic configurations endowed with a distance metric. The metric on configurations enables fast database lookup and robust interpolation of the stored forces. We discuss mechanisms for the database to adapt to the needs of the evolving dynamics, while maintaining accuracy, and other extensions of the basic algorithm. PMID:26669825
Parallel Molecular Dynamics Program for Molecules
Energy Science and Technology Software Center (ESTSC)
1995-03-07
ParBond is a parallel classical molecular dynamics code that models bonded molecular systems, typically of an organic nature. It uses classical force fields for both non-bonded Coulombic and Van der Waals interactions and for 2-, 3-, and 4-body bonded (bond, angle, dihedral, and improper) interactions. It integrates Newton''s equation of motion for the molecular system and evaluates various thermodynamical properties of the system as it progresses.
Fermionic Molecular Dynamics for Nuclear Dynamics and Thermodynamics
NASA Astrophysics Data System (ADS)
Hasnaoui, K. H. O.; Chomaz, Ph; Gulminelli, F.
A new Fermionic Molecular Dynamics (FMD) model based on a Skyrme functional is proposed in this paper. After introducing the basic formalism, some first applications to nuclear structure and nuclear thermodynamics are presented.
Molecular dynamics: A stitch in time
NASA Astrophysics Data System (ADS)
Deupi, Xavier
2014-01-01
Lengthy molecular dynamics simulations of complex systems at the atomic scale usually require supercomputers. Now, by stitching together many shorter independent simulations run 'in the cloud', this requirement has been circumvented, allowing two milliseconds of the dynamics of a G-protein-coupled receptor to be simulated.
NASA Astrophysics Data System (ADS)
Sellers, Michael; Lisal, Martin; Brennan, John
2015-06-01
Investigating the ability of a molecular model to accurately represent a real material is crucial to model development and use. When the model simulates materials in extreme conditions, one such property worth evaluating is the phase transition point. However, phase transitions are often overlooked or approximated because of difficulty or inaccuracy when simulating them. Techniques such as super-heating or super-squeezing a material to induce a phase change suffer from inherent timescale limitations leading to ``over-driving,'' and dual-phase simulations require many long-time runs to seek out what frequently results in an inexact location of phase-coexistence. We present a compilation of methods for the determination of solid-solid and solid-liquid phase transition points through the accurate calculation of the chemical potential. The methods are applied to the Smith-Bharadwaj atomistic potential's representation of cyclotrimethylene trinitramine (RDX) to accurately determine its melting point (Tm) and the alpha to gamma solid phase transition pressure. We also determine Tm for a coarse-grain model of RDX, and compare its value to experiment and atomistic counterpart. All methods are employed via the LAMMPS simulator, resulting in 60-70 simulations that total 30-50 ns. Approved for public release. Distribution is unlimited.
Ab initio molecular dynamics: concepts, recent developments, and future trends.
Iftimie, Radu; Minary, Peter; Tuckerman, Mark E
2005-05-10
The methodology of ab initio molecular dynamics, wherein finite-temperature dynamical trajectories are generated by using forces computed "on the fly" from electronic structure calculations, has had a profound influence in modern theoretical research. Ab initio molecular dynamics allows chemical processes in condensed phases to be studied in an accurate and unbiased manner, leading to new paradigms in the elucidation of microscopic mechanisms, rationalization of experimental data, and testable predictions of new phenomena. The purpose of this work is to give a brief introduction to the technique and to review several important recent developments in the field. Several illustrative examples showing the power of the technique have been chosen. Perspectives on future directions in the field also will be given. PMID:15870204
Reinbolt, Jeffrey A.; Haftka, Raphael T.; Chmielewski, Terese L.; Fregly, Benjamin J.
2013-01-01
Variations in joint parameter values (axis positions and orientations in body segments) and inertial parameter values (segment masses, mass centers, and moments of inertia) as well as kinematic noise alter the results of inverse dynamics analyses of gait. Three-dimensional linkage models with joint constraints have been proposed as one way to minimize the effects of noisy kinematic data. Such models can also be used to perform gait optimizations to predict post-treatment function given pre-treatment gait data. This study evaluates whether accurate patient-specific joint and inertial parameter values are needed in three-dimensional linkage models to produce accurate inverse dynamics results for gait. The study was performed in two stages. First, we used optimization analyses to evaluate whether patient-specific joint and inertial parameter values can be calibrated accurately from noisy kinematic data, and second, we used Monte Carlo analyses to evaluate how errors in joint and inertial parameter values affect inverse dynamics calculations. Both stages were performed using a dynamic, 27 degree-of-freedom, full-body linkage model and synthetic (i.e., computer generated) gait data corresponding to a nominal experimental gait motion. In general, joint but not inertial parameter values could be found accurately from noisy kinematic data. Root-mean-square (RMS) errors were 3° and 4 mm for joint parameter values and 1 kg, 22 mm, and 74,500 kg*mm2 for inertial parameter values. Furthermore, errors in joint but not inertial parameter values had a significant effect on calculated lower-extremity inverse dynamics joint torques. The worst RMS torque error averaged 4% bodyweight*height (BW*H) due to joint parameter variations but less than 0.25% BW*H due to inertial parameter variations. These results suggest that inverse dynamics analyses of gait utilizing linkage models with joint constraints should calibrate the model’s joint parameter values to obtain accurate joint
Dynamic molecular crystals with switchable physical properties.
Sato, Osamu
2016-06-21
The development of molecular materials whose physical properties can be controlled by external stimuli - such as light, electric field, temperature, and pressure - has recently attracted much attention owing to their potential applications in molecular devices. There are a number of ways to alter the physical properties of crystalline materials. These include the modulation of the spin and redox states of the crystal's components, or the incorporation within the crystalline lattice of tunable molecules that exhibit stimuli-induced changes in their molecular structure. A switching behaviour can also be induced by changing the molecular orientation of the crystal's components, even in cases where the overall molecular structure is not affected. Controlling intermolecular interactions within a molecular material is also an effective tool to modulate its physical properties. This Review discusses recent advances in the development of such stimuli-responsive, switchable crystalline compounds - referred to here as dynamic molecular crystals - and suggests how different approaches can serve to prepare functional materials. PMID:27325090
Molecular potentials and relaxation dynamics
NASA Astrophysics Data System (ADS)
Karo, A. M.
1981-03-01
The use of empirical pseudopotentials, in evaluating interatomic potentials, provides and inexpensive and convenient method for obtaining highly accurate potential curves and permits the modeling of core-valence correlation, and the inclusion of relativistic effects when these are significant. As an example, recent calculations of the chi and the a 3 Sigma + states of LiH, NaH, KH, RbH, and CsH and the chi 2 Sigma + states of their anions are discussed. Pseudopotentials, including core polarization terms, have been used to replace the core electrons, and this has been coupled with the development of compact, highly-optimized basis sets for the corresponding one- and two-electron atoms. Comparisons of the neutral potential curves with experiment and other ab initio calculations show good agreement (within 1000/cm over most of the potential curves) with the difference curves being considerably more accurate.
Molecular dynamics simulations: advances and applications
Hospital, Adam; Goñi, Josep Ramon; Orozco, Modesto; Gelpí, Josep L
2015-01-01
Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.
Molecular dynamics simulations of large macromolecular complexes
Perilla, Juan R.; Goh, Boon Chong; Cassidy, C. Keith; Liu, Bo; Bernardi, Rafael C.; Rudack, Till; Yu, Hang; Wu, Zhe; Schulten, Klaus
2015-01-01
Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone. PMID:25845770
NASA Astrophysics Data System (ADS)
Mills, Andrew A.; Ford, Kyle B.; Kreckel, Holger; Perera, Manori; Crabtree, Kyle N.; McCall, Benjamin J.
2009-06-01
With the advent of Herschel and SOFIA, laboratory methods capable of providing molecular rest frequencies in the terahertz and sub-millimeter regime are increasingly important. As of yet, it has been difficult to perform spectroscopy in this wavelength region due to the limited availability of radiation sources, optics, and detectors. Our goal is to provide accurate THz rest frequencies for molecular ions by combining previously recorded microwave transitions with combination differences obtained from high precision mid-IR spectroscopy. We are constructing a Sensitive Resolved Ion Beam Spectroscopy setup which will harness the benefits of kinematic compression in a molecular ion beam to enable very high resolution spectroscopy. This ion beam is interrogated by continuous-wave cavity ringdown spectroscopy using a home-made widely tunable difference frequency laser that utilizes two near-IR lasers and a periodically-poled lithium niobate crystal. Here, we report our efforts to optimize our ion beam spectrometer and to perform high-precision and high-accuracy frequency measurements using an optical frequency comb. footnote
A large catalog of accurate distances to molecular clouds from PS1 photometry
Schlafly, E. F.; Rix, H.-W.; Martin, N. F.; Green, G.; Finkbeiner, D. P.; Bell, E. F.; Burgett, W. S.; Chambers, K. C.; Hodapp, K. W.; Kaiser, N.; Magnier, E. A.; Tonry, J. L.; Draper, P. W.; Metcalfe, N.; Price, P. A.
2014-05-01
Distance measurements to molecular clouds are important but are often made separately for each cloud of interest, employing very different data and techniques. We present a large, homogeneous catalog of distances to molecular clouds, most of which are of unprecedented accuracy. We determine distances using optical photometry of stars along lines of sight toward these clouds, obtained from PanSTARRS-1. We simultaneously infer the reddenings and distances to these stars, tracking the full probability distribution function using a technique presented in Green et al. We fit these star-by-star measurements using a simple dust screen model to find the distance to each cloud. We thus estimate the distances to almost all of the clouds in the Magnani et al. catalog, as well as many other well-studied clouds, including Orion, Perseus, Taurus, Cepheus, Polaris, California, and Monoceros R2, avoiding only the inner Galaxy. Typical statistical uncertainties in the distances are 5%, though the systematic uncertainty stemming from the quality of our stellar models is about 10%. The resulting catalog is the largest catalog of accurate, directly measured distances to molecular clouds. Our distance estimates are generally consistent with available distance estimates from the literature, though in some cases the literature estimates are off by a factor of more than two.
NASA Astrophysics Data System (ADS)
Garrison, Stephen L.
2005-07-01
The combination of molecular simulations and potentials obtained from quantum chemistry is shown to be able to provide reasonably accurate thermodynamic property predictions. Gibbs ensemble Monte Carlo simulations are used to understand the effects of small perturbations to various regions of the model Lennard-Jones 12-6 potential. However, when the phase behavior and second virial coefficient are scaled by the critical properties calculated for each potential, the results obey a corresponding states relation suggesting a non-uniqueness problem for interaction potentials fit to experimental phase behavior. Several variations of a procedure collectively referred to as quantum mechanical Hybrid Methods for Interaction Energies (HM-IE) are developed and used to accurately estimate interaction energies from CCSD(T) calculations with a large basis set in a computationally efficient manner for the neon-neon, acetylene-acetylene, and nitrogen-benzene systems. Using these results and methods, an ab initio, pairwise-additive, site-site potential for acetylene is determined and then improved using results from molecular simulations using this initial potential. The initial simulation results also indicate that a limited range of energies important for accurate phase behavior predictions. Second virial coefficients calculated from the improved potential indicate that one set of experimental data in the literature is likely erroneous. This prescription is then applied to methanethiol. Difficulties in modeling the effects of the lone pair electrons suggest that charges on the lone pair sites negatively impact the ability of the intermolecular potential to describe certain orientations, but that the lone pair sites may be necessary to reasonably duplicate the interaction energies for several orientations. Two possible methods for incorporating the effects of three-body interactions into simulations within the pairwise-additivity formulation are also developed. A low density
Molecular potentials and relaxation dynamics
Karo, A.M.
1981-03-27
The use of empirical pseudopotentials, in evaluating interatomic potentials, provides an inexpensive and convenient method for obtaining highly accurate potential curves and permits the modeling of core-valence correlation, and the inclusion of relativistic effects when these are significant. As an example, recent calculations of the chi/sup 1/..sigma../sup +/ and a/sup 3/..sigma../sup +/ states of LiH, NaH, KH, RbH, and CsH and the chi/sup 2/..sigma../sup +/ states of their anions are discussed. Pseudopotentials, including core polarization terms, have been used to replace the core electrons, and this has been coupled with the development of compact, highly-optimized basis sets for the corresponding one- and two-electron atoms. Comparisons of the neutral potential curves with experiment and other ab initio calculations show good agreement (within 1000 cm/sup -1/ over most of the potential curves) with the difference curves being considerably more accurate.
Accurate force fields and methods for modelling organic molecular crystals at finite temperatures.
Nyman, Jonas; Pundyke, Orla Sheehan; Day, Graeme M
2016-06-21
We present an assessment of the performance of several force fields for modelling intermolecular interactions in organic molecular crystals using the X23 benchmark set. The performance of the force fields is compared to several popular dispersion corrected density functional methods. In addition, we present our implementation of lattice vibrational free energy calculations in the quasi-harmonic approximation, using several methods to account for phonon dispersion. This allows us to also benchmark the force fields' reproduction of finite temperature crystal structures. The results demonstrate that anisotropic atom-atom multipole-based force fields can be as accurate as several popular DFT-D methods, but have errors 2-3 times larger than the current best DFT-D methods. The largest error in the examined force fields is a systematic underestimation of the (absolute) lattice energy. PMID:27230942
Spinelli, Orietta; Rambaldi, Alessandro; Rigo, Francesca; Zanghì, Pamela; D'Agostini, Elena; Amicarelli, Giulia; Colotta, Francesco; Divona, Mariadomenica; Ciardi, Claudia; Coco, Francesco Lo; Minnucci, Giulia
2015-01-01
The diagnostic work-up of acute promyelocytic leukemia (APL) includes the cytogenetic demonstration of the t(15;17) translocation and/or the PML-RARA chimeric transcript by RQ-PCR or RT-PCR. This latter assays provide suitable results in 3-6 hours. We describe here two new, rapid and specific assays that detect PML-RARA transcripts, based on the RT-QLAMP (Reverse Transcription-Quenching Loop-mediated Isothermal Amplification) technology in which RNA retrotranscription and cDNA amplification are carried out in a single tube with one enzyme at one temperature, in fluorescence and real time format. A single tube triplex assay detects bcr1 and bcr3 PML-RARA transcripts along with GUS housekeeping gene. A single tube duplex assay detects bcr2 and GUSB. In 73 APL cases, these assays detected in 16 minutes bcr1, bcr2 and bcr3 transcripts. All 81 non-APL samples were negative by RT-QLAMP for chimeric transcripts whereas GUSB was detectable. In 11 APL patients in which RT-PCR yielded equivocal breakpoint type results, RT-QLAMP assays unequivocally and accurately defined the breakpoint type (as confirmed by sequencing). Furthermore, RT-QLAMP could amplify two bcr2 transcripts with particularly extended PML exon 6 deletions not amplified by RQ-PCR. RT-QLAMP reproducible sensitivity is 10−3 for bcr1 and bcr3 and 10−2 for bcr2 thus making this assay particularly attractive at diagnosis and leaving RQ-PCR for the molecular monitoring of minimal residual disease during the follow up. In conclusion, PML-RARA RT-QLAMP compared to RT-PCR or RQ-PCR is a valid improvement to perform rapid, simple and accurate molecular diagnosis of APL. PMID:25815362
Clustering Molecular Dynamics Trajectories for Optimizing Docking Experiments
De Paris, Renata; Quevedo, Christian V.; Ruiz, Duncan D.; Norberto de Souza, Osmar; Barros, Rodrigo C.
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
Clustering molecular dynamics trajectories for optimizing docking experiments.
De Paris, Renata; Quevedo, Christian V; Ruiz, Duncan D; Norberto de Souza, Osmar; Barros, Rodrigo C
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
PyVCI: A flexible open-source code for calculating accurate molecular infrared spectra
NASA Astrophysics Data System (ADS)
Sibaev, Marat; Crittenden, Deborah L.
2016-06-01
The PyVCI program package is a general purpose open-source code for simulating accurate molecular spectra, based upon force field expansions of the potential energy surface in normal mode coordinates. It includes harmonic normal coordinate analysis and vibrational configuration interaction (VCI) algorithms, implemented primarily in Python for accessibility but with time-consuming routines written in C. Coriolis coupling terms may be optionally included in the vibrational Hamiltonian. Non-negligible VCI matrix elements are stored in sparse matrix format to alleviate the diagonalization problem. CPU and memory requirements may be further controlled by algorithmic choices and/or numerical screening procedures, and recommended values are established by benchmarking using a test set of 44 molecules for which accurate analytical potential energy surfaces are available. Force fields in normal mode coordinates are obtained from the PyPES library of high quality analytical potential energy surfaces (to 6th order) or by numerical differentiation of analytic second derivatives generated using the GAMESS quantum chemical program package (to 4th order).
Molecular dynamics simulations of a lithium/sodium carbonate mixture.
Ottochian, Alistar; Ricca, Chiara; Labat, Frederic; Adamo, Carlo
2016-03-01
The diffusion and ionic conductivity of Li x Na1-x CO3 salt mixtures were studied by means of Molecular Dynamics (MD) simulations, using the Janssen and Tissen model (Janssen and Tissen, Mol Simul 5:83-98; 1990). These salts have received particular attention due to their central role in fuel cells technology, and reliable numerical methods that could perform as important interpretative tool of experimental data are thus required but still lacking. The chosen computational model nicely reproduces the main structural behaviour of the pure Li2CO3, Na2CO3 and K2CO3 carbonates, but also of their Li/K and Li/Na mixtures. However, it fails to accurately describe dynamic properties such as activation energies of diffusion and conduction processes, outlining the need to develop more accurate models for the simulation of molten salt carbonates. PMID:26897519
Accurate and molecular-size-tolerant NMR quantitation of diverse components in solution
Okamura, Hideyasu; Nishimura, Hiroshi; Nagata, Takashi; Kigawa, Takanori; Watanabe, Takashi; Katahira, Masato
2016-01-01
Determining the amount of each component of interest in a mixture is a fundamental first step in characterizing the nature of the solution and to develop possible means of utilization of its components. Similarly, determining the composition of units in complex polymers, or polymer mixtures, is crucial. Although NMR is recognized as one of the most powerful methods to achieve this and is widely used in many fields, variation in the molecular sizes or the relative mobilities of components skews quantitation due to the size-dependent decay of magnetization. Here, a method to accurately determine the amount of each component by NMR was developed. This method was validated using a solution that contains biomass-related components in which the molecular sizes greatly differ. The method is also tolerant of other factors that skew quantitation such as variation in the one-bond C–H coupling constant. The developed method is the first and only way to reliably overcome the skewed quantitation caused by several different factors to provide basic information on the correct amount of each component in a solution. PMID:26883279
Molecular dynamic simulations of ocular tablet dissolution.
Ru, Qian; Fadda, Hala M; Li, Chung; Paul, Daniel; Khaw, Peng T; Brocchini, Steve; Zloh, Mire
2013-11-25
Small tablets for implantation into the subconjunctival space in the eye are being developed to inhibit scarring after glaucoma filtration surgery (GFS). There is a need to evaluate drug dissolution at the molecular level to determine how the chemical structure of the active may correlate with dissolution in the nonsink conditions of the conjunctival space. We conducted molecular dynamics simulations to study the dissolution process of tablets derived from two drugs that can inhibit fibrosis after GFS, 5-fluorouracil (5-FU) and the matrix metalloprotease inhibitor (MMPi), ilomastat. The dissolution was simulated in the presence of simple point charge (SPC) water molecules, and the liquid turnover of the aqueous humor in the subconjunctival space was simulated by removal of the dissolved drug molecules at regular intervals and replacement by new water molecules. At the end of the simulation, the total molecular solvent accessible surface area of 5-FU tablets increased by 60 times more than that of ilomastat as a result of tablet swelling and release of molecules into solution. The tablet dissolution pattern shown in our molecular dynamic simulations tends to correlate with experimental release profiles. This work indicates that a series of molecular dynamic simulations can be used to predict the influence of the molecular properties of a drug on its dissolution profile and could be useful during preformulation where sufficient amounts of the drug are not always available to perform dissolution studies. PMID:24073784
Molecular Scale Dynamics of Large Ring Polymers
NASA Astrophysics Data System (ADS)
Gooßen, S.; Brás, A. R.; Krutyeva, M.; Sharp, M.; Falus, P.; Feoktystov, A.; Gasser, U.; Pyckhout-Hintzen, W.; Wischnewski, A.; Richter, D.
2014-10-01
We present neutron scattering data on the structure and dynamics of melts from polyethylene oxide rings with molecular weights up to ten times the entanglement mass of the linear counterpart. The data reveal a very compact conformation displaying a structure approaching a mass fractal, as hypothesized by recent simulation work. The dynamics is characterized by a fast Rouse relaxation of subunits (loops) and a slower dynamics displaying a lattice animal-like loop displacement. The loop size is an intrinsic property of the ring architecture and is independent of molecular weight. This is the first experimental observation of the space-time evolution of segmental motion in ring polymers illustrating the dynamic consequences of their topology that is unique among all polymeric systems of any other known architecture.
Dynamic signature of molecular association in methanol
NASA Astrophysics Data System (ADS)
Bertrand, C. E.; Self, J. L.; Copley, J. R. D.; Faraone, A.
2016-07-01
Quasielastic neutron scattering measurements and molecular dynamics simulations were combined to investigate the collective dynamics of deuterated methanol, CD3OD. In the experimentally determined dynamic structure factor, a slow, non-Fickian mode was observed in addition to the standard density-fluctuation heat mode. The simulation results indicate that the slow dynamical process originates from the hydrogen bonding of methanol molecules. The qualitative behavior of this mode is similar to the previously observed α-relaxation in supercooled water [M. C. Bellissent-Funel et al., Phys. Rev. Lett. 85, 3644 (2000)] which also originates from the formation and dissolution of hydrogen-bonded associates (supramolecular clusters). In methanol, however, this mode is distinguishable well above the freezing transition. This finding indicates that an emergent slow mode is not unique to supercooled water, but may instead be a general feature of hydrogen-bonding liquids and associating molecular liquids.
Dynamic signature of molecular association in methanol.
Bertrand, C E; Self, J L; Copley, J R D; Faraone, A
2016-07-01
Quasielastic neutron scattering measurements and molecular dynamics simulations were combined to investigate the collective dynamics of deuterated methanol, CD3OD. In the experimentally determined dynamic structure factor, a slow, non-Fickian mode was observed in addition to the standard density-fluctuation heat mode. The simulation results indicate that the slow dynamical process originates from the hydrogen bonding of methanol molecules. The qualitative behavior of this mode is similar to the previously observed α-relaxation in supercooled water [M. C. Bellissent-Funel et al., Phys. Rev. Lett. 85, 3644 (2000)] which also originates from the formation and dissolution of hydrogen-bonded associates (supramolecular clusters). In methanol, however, this mode is distinguishable well above the freezing transition. This finding indicates that an emergent slow mode is not unique to supercooled water, but may instead be a general feature of hydrogen-bonding liquids and associating molecular liquids. PMID:27394112
Molecular scale dynamics of large ring polymers.
Gooßen, S; Brás, A R; Krutyeva, M; Sharp, M; Falus, P; Feoktystov, A; Gasser, U; Pyckhout-Hintzen, W; Wischnewski, A; Richter, D
2014-10-17
We present neutron scattering data on the structure and dynamics of melts from polyethylene oxide rings with molecular weights up to ten times the entanglement mass of the linear counterpart. The data reveal a very compact conformation displaying a structure approaching a mass fractal, as hypothesized by recent simulation work. The dynamics is characterized by a fast Rouse relaxation of subunits (loops) and a slower dynamics displaying a lattice animal-like loop displacement. The loop size is an intrinsic property of the ring architecture and is independent of molecular weight. This is the first experimental observation of the space-time evolution of segmental motion in ring polymers illustrating the dynamic consequences of their topology that is unique among all polymeric systems of any other known architecture. PMID:25361284
Semiclassical guided optimal control of molecular dynamics
Kondorskiy, A.; Mil'nikov, G.; Nakamura, H.
2005-10-15
An efficient semiclassical optimal control theory applicable to multidimensional systems is formulated for controlling wave packet dynamics on a single adiabatic potential energy surface. The approach combines advantages of different formulations of optimal control theory: quantum and classical on one hand and global and local on the other. Numerical applications to the control of HCN-CNH isomerization demonstrate that this theory can provide an efficient tool to manipulate molecular dynamics of many degrees of freedom by laser pulses.
An accurate dynamical electron diffraction algorithm for reflection high-energy electron diffraction
NASA Astrophysics Data System (ADS)
Huang, J.; Cai, C. Y.; Lv, C. L.; Zhou, G. W.; Wang, Y. G.
2015-12-01
The conventional multislice method (CMS) method, one of the most popular dynamical electron diffraction calculation procedures in transmission electron microscopy, was introduced to calculate reflection high-energy electron diffraction (RHEED) as it is well adapted to deal with the deviations from the periodicity in the direction parallel to the surface. However, in the present work, we show that the CMS method is no longer sufficiently accurate for simulating RHEED with the accelerating voltage 3-100 kV because of the high-energy approximation. An accurate multislice (AMS) method can be an alternative for more accurate RHEED calculations with reasonable computing time. A detailed comparison of the numerical calculation of the AMS method and the CMS method is carried out with respect to different accelerating voltages, surface structure models, Debye-Waller factors and glancing angles.
Shapiro like steps reveals molecular nanomagnets’ spin dynamics
Abdollahipour, Babak; Abouie, Jahanfar Ebrahimi, Navid
2015-09-15
We present an accurate way to detect spin dynamics of a nutating molecular nanomagnet by inserting it in a tunnel Josephson junction and studying the current voltage (I-V) characteristic. The spin nutation of the molecular nanomagnet is generated by applying two circularly polarized magnetic fields. We demonstrate that modulation of the Josephson current by the nutation of the molecular nanomagnet’s spin appears as a stepwise structure like Shapiro steps in the I-V characteristic of the junction. Width and heights of these Shapiro-like steps are determined by two parameters of the spin nutation, frequency and amplitude of the nutation, which are simply tuned by the applied magnetic fields.
Danshita, Ippei; Polkovnikov, Anatoli
2010-09-01
We study the quantum dynamics of supercurrents of one-dimensional Bose gases in a ring optical lattice to verify instanton methods applied to coherent macroscopic quantum tunneling (MQT). We directly simulate the real-time quantum dynamics of supercurrents, where a coherent oscillation between two macroscopically distinct current states occurs due to MQT. The tunneling rate extracted from the coherent oscillation is compared with that given by the instanton method. We find that the instanton method is quantitatively accurate when the effective Planck's constant is sufficiently small. We also find phase slips associated with the oscillations.
Optimizing replica exchange moves for molecular dynamics.
Nadler, Walter; Hansmann, Ulrich H E
2007-11-01
We sketch the statistical physics framework of the replica exchange technique when applied to molecular dynamics simulations. In particular, we draw attention to generalized move sets that allow a variety of optimizations as well as new applications of the method. PMID:18233794
Molecular dynamics calculations of nuclear stimulated desorption
Glikman, E.; Kelson, I. ); Doan, N.V. )
1991-09-01
Molecular dynamics calculations of nuclear stimulated desorption are carried out for a palladium crystal containing radioactive palladium atoms. The total desorption probability from various sites are computed, as well as the angular distribution of the desorbing atoms. The implications of the results to different experimental scenarios are discussed.
Reaction dynamics in polyatomic molecular systems
Miller, W.H.
1993-12-01
The goal of this program is the development of theoretical methods and models for describing the dynamics of chemical reactions, with specific interest for application to polyatomic molecular systems of special interest and relevance. There is interest in developing the most rigorous possible theoretical approaches and also in more approximate treatments that are more readily applicable to complex systems.
Molecular Dynamics Simulations of Graphene Oxide Frameworks
Zhu, Pan; Sumpter, Bobby G; Meunier, V.; Nicolai, Adrien
2013-01-01
We use quantum mechanical calculations to develop a full set of force field parameters in order to perform molecular dynamics simulations to understand and optimize the molecular storage properties inside Graphene Oxide Frameworks (GOFs). A set of boron-related parameters for commonly used empirical force fields is determined to describe the non-bonded and bonded interactions between linear boronic acid linkers and graphene sheets of GOF materials. The transferability of the parameters is discussed and their validity is quantified by comparing quantum mechanical and molecular mechanical structural and vibrational properties. The application of the model to the dynamics of water inside the GOFs reveals significant variations in structural flexibility of GOF depending on the linker density, which is shown to be usable as a tuning parameter for desired diffusion properties.
Molecular dynamics simulations of weak detonations.
Am-Shallem, Morag; Zeiri, Yehuda; Zybin, Sergey V; Kosloff, Ronnie
2011-12-01
Detonation of a three-dimensional reactive nonisotropic molecular crystal is modeled using molecular dynamics simulations. The detonation process is initiated by an impulse, followed by the creation of a stable fast reactive shock wave. The terminal shock velocity is independent of the initiation conditions. Further analysis shows supersonic propagation decoupled from the dynamics of the decomposed material left behind the shock front. The dependence of the shock velocity on crystal nonlinear compressibility resembles solitary behavior. These properties categorize the phenomena as a weak detonation. The dependence of the detonation wave on microscopic potential parameters was investigated. An increase in detonation velocity with the reaction exothermicity reaching a saturation value is observed. In all other respects the model crystal exhibits typical properties of a molecular crystal. PMID:22304055
Molecular dynamics of PLK1 during mitosis
Schmucker, Stephane; Sumara, Izabela
2014-01-01
Polo-like kinase 1 (PLK1) is a key regulator of eukaryotic cell division. During mitosis, dynamic regulation of PLK1 is crucial for its roles in centrosome maturation, spindle assembly, microtubule–kinetochore attachment, and cytokinesis. Similar to other members of the PLK family, the molecular architecture of PLK1 protein is characterized by 2 domains—the kinase domain and the regulatory substrate-binding domain (polo-box domain)—that cooperate and control PLK1 function during mitosis. Mitotic cells employ many layers of regulation to activate and target PLK1 to different cellular structures in a timely manner. During the last decade, numerous studies have shed light on the precise molecular mechanisms orchestrating the mitotic activity of PLK1 in time and space. This review aims to discuss available data and concepts related to regulation of the molecular dynamics of human PLK1 during mitotic progression. PMID:27308323
Approximate but accurate quantum dynamics from the Mori formalism: I. Nonequilibrium dynamics.
Montoya-Castillo, Andrés; Reichman, David R
2016-05-14
We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics. PMID:27179468
Approximate but accurate quantum dynamics from the Mori formalism: I. Nonequilibrium dynamics
NASA Astrophysics Data System (ADS)
Montoya-Castillo, Andrés; Reichman, David R.
2016-05-01
We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics.
Thompson, Andrew R; Binder, Benjamin P; McCaffrey, Jesse E; Svensson, Bengt; Thomas, David D
2015-01-01
While EPR allows for the characterization of protein structure and function due to its exquisite sensitivity to spin label dynamics, orientation, and distance, these measurements are often limited in sensitivity due to the use of labels that are attached via flexible monofunctional bonds, incurring additional disorder and nanosecond dynamics. In this chapter, we present methods for using a bifunctional spin label (BSL) to measure muscle protein structure and dynamics. We demonstrate that bifunctional attachment eliminates nanosecond internal rotation of the spin label, thereby allowing the accurate measurement of protein backbone rotational dynamics, including microsecond-to-millisecond motions by saturation transfer EPR. BSL also allows for accurate determination of helix orientation and disorder in mechanically and magnetically aligned systems, due to the label's stereospecific attachment. Similarly, labeling with a pair of BSL greatly enhances the resolution and accuracy of distance measurements measured by double electron-electron resonance (DEER). Finally, when BSL is applied to a protein with high helical content in an assembly with high orientational order (e.g., muscle fiber or membrane), two-probe DEER experiments can be combined with single-probe EPR experiments on an oriented sample in a process we call BEER, which has the potential for ab initio high-resolution structure determination. PMID:26477249
Structural and dynamical properties of hot dense matter by a Thomas-Fermi-Dirac molecular dynamics
NASA Astrophysics Data System (ADS)
Lambert, F.; Clérouin, J.; Mazevet, S.
2006-09-01
We use a model combining, in a consistent way, orbital-free density functional theory (OF-DFT) and molecular dynamics (MD), to compute the thermodynamical, structural and dynamical properties of Fe and Au plasmas at conditions relevant to astrophysics and inertial confinement fusion (ICF). The newly developed parallel numerical scheme presented here allows to propagate hundreds of particles and to obtain accurate transport properties. This allows us to investigate the validity of the commonly used one-component plasma (OCP) model in predicting the pair correlation, the diffusion and viscosity coefficients for these two high-temperature high-density plasmas.
Accurate calculation of binding energies for molecular clusters - Assessment of different models
NASA Astrophysics Data System (ADS)
Friedrich, Joachim; Fiedler, Benjamin
2016-06-01
In this work we test different strategies to compute high-level benchmark energies for medium-sized molecular clusters. We use the incremental scheme to obtain CCSD(T)/CBS energies for our test set and carefully validate the accuracy for binding energies by statistical measures. The local errors of the incremental scheme are <1 kJ/mol. Since they are smaller than the basis set errors, we obtain higher total accuracy due to the applicability of larger basis sets. The final CCSD(T)/CBS benchmark values are ΔE = - 278.01 kJ/mol for (H2O)10, ΔE = - 221.64 kJ/mol for (HF)10, ΔE = - 45.63 kJ/mol for (CH4)10, ΔE = - 19.52 kJ/mol for (H2)20 and ΔE = - 7.38 kJ/mol for (H2)10 . Furthermore we test state-of-the-art wave-function-based and DFT methods. Our benchmark data will be very useful for critical validations of new methods. We find focal-point-methods for estimating CCSD(T)/CBS energies to be highly accurate and efficient. For foQ-i3CCSD(T)-MP2/TZ we get a mean error of 0.34 kJ/mol and a standard deviation of 0.39 kJ/mol.
2012-01-01
A natural bond orbital (NBO) analysis of unpaired electron spin density in metalloproteins is presented, which allows a fast and robust calculation of paramagnetic NMR parameters. Approximately 90% of the unpaired electron spin density occupies metal–ligand NBOs, allowing the majority of the density to be modeled by only a few NBOs that reflect the chemical bonding environment. We show that the paramagnetic relaxation rate of protons can be calculated accurately using only the metal–ligand NBOs and that these rates are in good agreement with corresponding rates measured experimentally. This holds, in particular, for protons of ligand residues where the point-dipole approximation breaks down. To describe the paramagnetic relaxation of heavy nuclei, also the electron spin density in the local orbitals must be taken into account. Geometric distance restraints for 15N can be derived from the paramagnetic relaxation enhancement and the Fermi contact shift when local NBOs are included in the analysis. Thus, the NBO approach allows us to include experimental paramagnetic NMR parameters of 15N nuclei as restraints in a structure optimization protocol. We performed a molecular dynamics simulation and structure determination of oxidized rubredoxin using the experimentally obtained paramagnetic NMR parameters of 15N. The corresponding structures obtained are in good agreement with the crystal structure of rubredoxin. Thus, the NBO approach allows an accurate description of the geometric structure and the dynamics of metalloproteins, when NMR parameters are available of nuclei in the immediate vicinity of the metal-site. PMID:22329704
Nonadiabatic molecular dynamics simulations: synergies between theory and experiments.
Tavernelli, Ivano
2015-03-17
Recent developments in nonadiabatic dynamics enabled ab inito simulations of complex ultrafast processes in the condensed phase. These advances have opened new avenues in the study of many photophysical and photochemical reactions triggered by the absorption of electromagnetic radiation. In particular, theoretical investigations can be combined with the most sophisticated femtosecond experimental techniques to guide the interpretation of measured time-resolved observables. At the same time, the availability of experimental data at high (spatial and time) resolution offers a unique opportunity for the benchmarking and the improvement of those theoretical models used to describe complex molecular systems in their natural environment. The established synergy between theory and experiments can produce a better understanding of new ultrafast physical and chemical processes at atomistic scale resolution. Furthermore, reliable ab inito molecular dynamics simulations can already be successfully employed as predictive tools to guide new experiments as well as the design of novel and better performing materials. In this paper, I will give a concise account on the state of the art of molecular dynamics simulations of complex molecular systems in their excited states. The principal aim of this approach is the description of a given system of interest under the most realistic ambient conditions including all environmental effects that influence experiments, for instance, the interaction with the solvent and with external time-dependent electric fields, temperature, and pressure. To this end, time-dependent density functional theory (TDDFT) is among the most efficient and accurate methods for the representation of the electronic dynamics, while trajectory surface hopping gives a valuable representation of the nuclear quantum dynamics in the excited states (including nonadiabatic effects). Concerning the environment and its effects on the dynamics, the quantum mechanics/molecular
2012-01-01
In the present work, we employ excited state accelerated ab initio molecular dynamics (A-AIMD) to efficiently study the excited state energy landscape and photophysical topology of a variety of molecular systems. In particular, we focus on two important challenges for the modeling of excited electronic states: (i) the identification and characterization of conical intersections and crossing seams, in order to predict different and often competing radiationless decay mechanisms, and (ii) the description of the solvent effect on the absorption and emission spectra of chemical species in solution. In particular, using as examples the Schiff bases formaldimine and salicylidenaniline, we show that A-AIMD can be readily employed to explore the conformational space around crossing seams in molecular systems with very different photochemistry. Using acetone in water as an example, we demonstrate that the enhanced configurational space sampling may be used to accurately and efficiently describe both the prominent features and line-shapes of absorption and emission spectra. PMID:22904696
Exciton dynamics in perturbed vibronic molecular aggregates
Brüning, C.; Wehner, J.; Hausner, J.; Wenzel, M.; Engel, V.
2015-01-01
A site specific perturbation of a photo-excited molecular aggregate can lead to a localization of excitonic energy. We investigate this localization dynamics for laser-prepared excited states. Changing the parameters of the electric field significantly influences the exciton localization which offers the possibility for a selective control of this process. This is demonstrated for aggregates possessing a single vibrational degree of freedom per monomer unit. It is shown that the effects identified for the molecular dimer can be generalized to larger aggregates with a high density of vibronic states. PMID:26798840
Exciton dynamics in perturbed vibronic molecular aggregates.
Brüning, C; Wehner, J; Hausner, J; Wenzel, M; Engel, V
2016-07-01
A site specific perturbation of a photo-excited molecular aggregate can lead to a localization of excitonic energy. We investigate this localization dynamics for laser-prepared excited states. Changing the parameters of the electric field significantly influences the exciton localization which offers the possibility for a selective control of this process. This is demonstrated for aggregates possessing a single vibrational degree of freedom per monomer unit. It is shown that the effects identified for the molecular dimer can be generalized to larger aggregates with a high density of vibronic states. PMID:26798840
Exploiting molecular dynamics in Nested Sampling simulations of small peptides
NASA Astrophysics Data System (ADS)
Burkoff, Nikolas S.; Baldock, Robert J. N.; Várnai, Csilla; Wild, David L.; Csányi, Gábor
2016-04-01
Nested Sampling (NS) is a parameter space sampling algorithm which can be used for sampling the equilibrium thermodynamics of atomistic systems. NS has previously been used to explore the potential energy surface of a coarse-grained protein model and has significantly outperformed parallel tempering when calculating heat capacity curves of Lennard-Jones clusters. The original NS algorithm uses Monte Carlo (MC) moves; however, a variant, Galilean NS, has recently been introduced which allows NS to be incorporated into a molecular dynamics framework, so NS can be used for systems which lack efficient prescribed MC moves. In this work we demonstrate the applicability of Galilean NS to atomistic systems. We present an implementation of Galilean NS using the Amber molecular dynamics package and demonstrate its viability by sampling alanine dipeptide, both in vacuo and implicit solvent. Unlike previous studies of this system, we present the heat capacity curves of alanine dipeptide, whose calculation provides a stringent test for sampling algorithms. We also compare our results with those calculated using replica exchange molecular dynamics (REMD) and find good agreement. We show the computational effort required for accurate heat capacity estimation for small peptides. We also calculate the alanine dipeptide Ramachandran free energy surface for a range of temperatures and use it to compare the results using the latest Amber force field with previous theoretical and experimental results.
Concise NMR approach for molecular dynamics characterizations in organic solids.
Aliev, Abil E; Courtier-Murias, Denis
2013-08-22
Molecular dynamics characterisations in solids can be carried out selectively using dipolar-dephasing experiments. Here we show that the introduction of a sum of Lorentzian and Gaussian functions greatly improve fittings of the "intensity versus time" data for protonated carbons in dipolar-dephasing experiments. The Lorentzian term accounts for remote intra- and intermolecular (1)H-(13)C dipole-dipole interactions, which vary from one molecule to another or for different carbons within the same molecule. Thus, by separating contributions from weak remote interactions, more accurate Gaussian decay constants, T(dd), can be extracted for directly bonded (1)H-(13)C dipole-dipole interactions. Reorientations of the (1)H-(13)C bonds lead to the increase of T(dd), and by measuring dipolar-dephasing constants, insight can be gained into dynamics in solids. We have demonstrated advantages of the method using comparative dynamics studies in the α and γ polymorphs of glycine, cyclic amino acids L-proline, DL-proline and trans-4-hydroxy-L-proline, the Ala residue in different dipeptides, as well as adamantane and hexamethylenetetramine. It was possible to distinguish subtle differences in dynamics of different carbon sites within a molecule in polymorphs and in L- and DL-forms. The presence of overall molecular motions is shown to lead to particularly large differences in dipolar-dephasing experiments. The differences in dynamics can be attributed to differences in noncovalent interactions. In the case of hexamethylenetetramine, for example, the presence of C-H···N interactions leads to nearly rigid molecules. Overall, the method allows one to gain insight into the role of noncovalent interactions in solids and their influence on the molecular dynamics. PMID:23879450
Molecular dynamics at constant Cauchy stress
NASA Astrophysics Data System (ADS)
Miller, Ronald E.; Tadmor, Ellad B.; Gibson, Joshua S.; Bernstein, Noam; Pavia, Fabio
2016-05-01
The Parrinello-Rahman algorithm for imposing a general state of stress in periodic molecular dynamics simulations is widely used in the literature and has been implemented in many readily available molecular dynamics codes. However, what is often overlooked is that this algorithm controls the second Piola-Kirchhoff stress as opposed to the true (Cauchy) stress. This can lead to misinterpretation of simulation results because (1) the true stress that is imposed during the simulation depends on the deformation of the periodic cell, (2) the true stress is potentially very different from the imposed second Piola-Kirchhoff stress, and (3) the true stress can vary significantly during the simulation even if the imposed second Piola-Kirchhoff is constant. We propose a simple modification to the algorithm that allows the true Cauchy stress to be controlled directly. We then demonstrate the efficacy of the new algorithm with the example of martensitic phase transformations under applied stress.
Molecular dynamics at constant Cauchy stress.
Miller, Ronald E; Tadmor, Ellad B; Gibson, Joshua S; Bernstein, Noam; Pavia, Fabio
2016-05-14
The Parrinello-Rahman algorithm for imposing a general state of stress in periodic molecular dynamics simulations is widely used in the literature and has been implemented in many readily available molecular dynamics codes. However, what is often overlooked is that this algorithm controls the second Piola-Kirchhoff stress as opposed to the true (Cauchy) stress. This can lead to misinterpretation of simulation results because (1) the true stress that is imposed during the simulation depends on the deformation of the periodic cell, (2) the true stress is potentially very different from the imposed second Piola-Kirchhoff stress, and (3) the true stress can vary significantly during the simulation even if the imposed second Piola-Kirchhoff is constant. We propose a simple modification to the algorithm that allows the true Cauchy stress to be controlled directly. We then demonstrate the efficacy of the new algorithm with the example of martensitic phase transformations under applied stress. PMID:27179471
Molecular dynamics study of cyclohexane interconversion
NASA Astrophysics Data System (ADS)
Wilson, Michael A.; Chandler, David
1990-12-01
Classical molecular dynamics calculations are reported for one C 6H 12 molecule in a bath of 250 CS 2 molecules at roomtemperature and liquid densities of 1.0, 1.3, 1.4 and 1.5 g/cm 3. The solvent contribution to the free energy of activation for the chair-boat isomerization has been determined to high accuracy. The transmission coefficient and reactive flux correlation functions have also been computed. The results obtained agree with earlier conclusions drawn from RISM integral equation calculations and stochastic molecular dynamics calculations. Namely, the solvent effect on the rate manifests a qualitative breakdown of transition state theory and the RRKM picture of unimolecular kinetics. Analysis of the activated trajectories indicate a significant degree of quasiperiodicity.
Molecular dynamics studies of polyurethane nanocomposite hydrogels
NASA Astrophysics Data System (ADS)
Strankowska, J.; Piszczyk, Ł.; Strankowski, M.; Danowska, M.; Szutkowski, K.; Jurga, S.; Kwela, J.
2013-10-01
Polyurethane PEO-based hydrogels have a broad range of biomedical applicability. They are attractive for drug-controlled delivery systems, surgical implants and wound healing dressings. In this study, a PEO based polyurethane hydrogels containing Cloisite® 30B, an organically modified clay mineral, was synthesized. Structure of nanocomposite hydrogels was determined using XRD technique. Its molecular dynamics was studied by means of NMR spectroscopy, DMA and DSC analysis. The mechanical properties and thermal stability of the systems were improved by incorporation of clay and controlled by varying the clay content in polymeric matrix. Molecular dynamics of polymer chains depends on interaction of Cloisite® 30B nanoparticles with soft segments of polyurethanes. The characteristic nanosize effect is observed.
New faster CHARMM molecular dynamics engine
Hynninen, Antti-Pekka; Crowley, Michael F
2014-01-01
We introduce a new faster molecular dynamics (MD) engine into the CHARMM software package. The new MD engine is faster both in serial (i.e., single CPU core) and parallel execution. Serial performance is approximately two times higher than in the previous version of CHARMM. The newly programmed parallelization method allows the MD engine to parallelize up to hundreds of CPU cores. PMID:24302199
Nanoindentation of Zr by molecular dynamics simulation
NASA Astrophysics Data System (ADS)
Lu (芦子哲), Zizhe; Chernatynskiy, Aleksandr; Noordhoek, Mark J.; Sinnott, Susan B.; Phillpot, Simon R.
2015-12-01
Molecular dynamics simulations of nanoindentation are used to study the deformation behaviors of single crystal Zr for four different surface orientations. The comparison of results for two different potentials, an embedded atom method potential and a charged optimized many body potential, reveals the influence of stable and unstable stacking fault energy on dislocation behaviors under nanoindentation. The load-displacement curve, hardness and deformation behaviors of the various surface orientations Zr are compared and the elastic and plastic deformation behaviors are analyzed.
Molecular dynamics modelling of solidification in metals
Boercker, D.B.; Belak, J.; Glosli, J.
1997-12-31
Molecular dynamics modeling is used to study the solidification of metals at high pressure and temperature. Constant pressure MD is applied to a simulation cell initially filled with both solid and molten metal. The solid/liquid interface is tracked as a function of time, and the data are used to estimate growth rates of crystallites at high pressure and temperature in Ta and Mg.
Molecular crowding and protein enzymatic dynamics.
Echeverria, Carlos; Kapral, Raymond
2012-05-21
The effects of molecular crowding on the enzymatic conformational dynamics and transport properties of adenylate kinase are investigated. This tridomain protein undergoes large scale hinge motions in the course of its enzymatic cycle and serves as prototype for the study of crowding effects on the cyclic conformational dynamics of proteins. The study is carried out at a mesoscopic level where both the protein and the solvent in which it is dissolved are treated in a coarse grained fashion. The amino acid residues in the protein are represented by a network of beads and the solvent dynamics is described by multiparticle collision dynamics that includes effects due to hydrodynamic interactions. The system is crowded by a stationary random array of hard spherical objects. Protein enzymatic dynamics is investigated as a function of the obstacle volume fraction and size. In addition, for comparison, results are presented for a modification of the dynamics that suppresses hydrodynamic interactions. Consistent with expectations, simulations of the dynamics show that the protein prefers a closed conformation for high volume fractions. This effect becomes more pronounced as the obstacle radius decreases for a given volume fraction since the average void size in the obstacle array is smaller for smaller radii. At high volume fractions for small obstacle radii, the average enzymatic cycle time and characteristic times of internal conformational motions of the protein deviate substantially from their values in solution or in systems with small density of obstacles. The transport properties of the protein are strongly affected by molecular crowding. Diffusive motion adopts a subdiffusive character and the effective diffusion coefficients can change by more than an order of magnitude. The orientational relaxation time of the protein is also significantly altered by crowding. PMID:22476233
Monoamine transporters: insights from molecular dynamics simulations
Grouleff, Julie; Ladefoged, Lucy Kate; Koldsø, Heidi; Schiøtt, Birgit
2015-01-01
The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson’s disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na+ symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors. PMID:26528185
Optimally designed fields for controlling molecular dynamics
NASA Astrophysics Data System (ADS)
Rabitz, Herschel
1991-10-01
This research concerns the development of molecular control theory techniques for designing optical fields capable of manipulating molecular dynamic phenomena. Although is has been long recognized that lasers should be capable of manipulating dynamic events, many frustrating years of intuitively driven laboratory studies only serve to illustrate the point that the task is complex and defies intuition. The principal new component in the present research is the recognition that this problem falls into the category of control theory and its inherent complexities require the use of modern control theory tools largely developed in the engineering disciplines. Thus, the research has initiated a transfer of the control theory concepts to the molecular scale. Although much contained effort will be needed to fully develop these concepts, the research in this grant set forth the basic components of the theory and carried out illustrative studies involving the design of optical fields capable of controlling rotational, vibrational and electronic degrees of freedom. Optimal control within the quantum mechanical molecular realm represents a frontier area with many possible ultimate applications. At this stage, the theoretical tools need to be joined with merging laboratory optical pulse shaping capabilities to illustrate the power of the concepts.
Structure and dynamics of complex liquid water: Molecular dynamics simulation
NASA Astrophysics Data System (ADS)
S, Indrajith V.; Natesan, Baskaran
2015-06-01
We have carried out detailed structure and dynamical studies of complex liquid water using molecular dynamics simulations. Three different model potentials, namely, TIP3P, TIP4P and SPC-E have been used in the simulations, in order to arrive at the best possible potential function that could reproduce the structure of experimental bulk water. All the simulations were performed in the NVE micro canonical ensemble using LAMMPS. The radial distribution functions, gOO, gOH and gHH and the self diffusion coefficient, Ds, were calculated for all three models. We conclude from our results that the structure and dynamical parameters obtained for SPC-E model matched well with the experimental values, suggesting that among the models studied here, the SPC-E model gives the best structure and dynamics of bulk water.
Dynamic Maintenance and Visualization of Molecular Surfaces
Bajaj, C L; Pascucci, V; Shamir, A; Holt, R J; Netravali, A N
2004-12-16
Molecular surface computations are often necessary in order to perform synthetic drug design. A critical step in this process is the computation and update of an exact boundary representation for the molecular surface (e.g. the Lee-Richards surface). In this paper they introduce efficient techniques for computing a molecular surface boundary representation as a set of NURBS (non-uniform rational B-splines) patches. This representation introduces for molecules the same geometric data structure used in the solid modeling community and enables immediate access to a wide range of modeling operations and techniques. Furthermore, this allows the use of any general solid modeling or visualization system as a molecular modeling interface. However, using such a representation in a molecular modeling environment raises several efficiency and update constraints, especially in a dynamic setting. For example, changes in the probe radius result in both geometric and topological changes to the set of patches. The techniques provide the option of trading accuracy of the representation for the efficiency of the computation, while still tracking the changes in the set of patches. In particular, they discuss two main classes of dynamic updates: one that keeps the topology of the molecular configuration fixed, and a more complicated case where the topology may be updated continuously. In general the generated output surface is represented in a format that can be loaded into standard solid modeling systems. It can also be directly triangulated or rendered, possibly at different levels of resolution, by a standard graphics library such as OpenGL without any additional effort.
The 2011 Dynamics of Molecular Collisions Conference
Nesbitt, David J.
2011-07-11
The Dynamics of Molecular Collisions Conference focuses on all aspects of molecular collisions--experimental & theoretical studies of elastic, inelastic, & reactive encounters involving atoms, molecules, ions, clusters, & surfaces--as well as half collisions--photodissociation, photo-induced reaction, & photodesorption. The scientific program for the meeting in 2011 included exciting advances in both the core & multidisciplinary forefronts of the study of molecular collision processes. Following the format of the 2009 meeting, we also invited sessions in special topics that involve interfacial dynamics, novel emerging spectroscopies, chemical dynamics in atmospheric, combustion & interstellar environments, as well as a session devoted to theoretical & experimental advances in ultracold molecular samples. Researchers working inside & outside the traditional core topics of the meeting are encouraged to join the conference. We invite contributions of work that seeks understanding of how inter & intra-molecular forces determine the dynamics of the phenomena under study. In addition to invited oral sessions & contributed poster sessions, the scientific program included a formal session consisting of five contributed talks selected from the submitted poster abstracts. The DMC has distinguished itself by having the Herschbach Medal Symposium as part of the meeting format. This tradition of the Herschbach Medal was first started in the 2007 meeting chaired by David Chandler, based on a generous donation of funds & artwork design by Professor Dudley Herschbach himself. There are two such awards made, one for experimental & one for theoretical contributions to the field of Molecular Collision Dynamics, broadly defined. The symposium is always held on the last night of the meeting & has the awardees are asked to deliver an invited lecture on their work. The 2011 Herschbach Medal was dedicated to the contributions of two long standing leaders in Chemical Physics, Professor
High temperature phonon dispersion in graphene using classical molecular dynamics
Anees, P. Panigrahi, B. K.; Valsakumar, M. C.
2014-04-24
Phonon dispersion and phonon density of states of graphene are calculated using classical molecular dynamics simulations. In this method, the dynamical matrix is constructed based on linear response theory by computing the displacement of atoms during the simulations. The computed phonon dispersions show excellent agreement with experiments. The simulations are done in both NVT and NPT ensembles at 300 K and found that the LO/TO modes are getting hardened at the Γ point. The NPT ensemble simulations capture the anharmonicity of the crystal accurately and the hardening of LO/TO modes is more pronounced. We also found that at 300 K the C-C bond length reduces below the equilibrium value and the ZA bending mode frequency becomes imaginary close to Γ along K-Γ direction, which indicates instability of the flat 2D graphene sheets.
Improved dynamic compensation for accurate cutting force measurements in milling applications
NASA Astrophysics Data System (ADS)
Scippa, A.; Sallese, L.; Grossi, N.; Campatelli, G.
2015-03-01
Accurate cutting-force measurements appear to be the key information in most of the machining related studies as they are fundamental in understanding the cutting processes, optimizing the cutting operations and evaluating the presence of instabilities that could affect the effectiveness of cutting processes. A variety of specifically designed transducers are commercially available nowadays and many different approaches in measuring cutting forces are presented in literature. The available transducers, though, express some limitations since they are conditioned by the vibration of the surrounding system and by the transducer's natural frequency. These parameters can drastically affect the measurement accuracy in some cases; hence an effective and accurate tool is required to compensate those dynamically induced errors in cutting force measurements. This work is aimed at developing and testing a compensation technique based on Kalman filter estimator. Two different approaches named "band-fitting" and "parallel elaboration" methods, have been developed to extend applications of this compensation technique, especially for milling purpose. The compensation filter has been designed upon the experimentally identified system's dynamic and its accuracy and effectiveness has been evaluated by numerical and experimental tests. Finally its specific application in cutting force measurements compensation is described.
Coarse-grained red blood cell model with accurate mechanical properties, rheology and dynamics.
Fedosov, Dmitry A; Caswell, Bruce; Karniadakis, George E
2009-01-01
We present a coarse-grained red blood cell (RBC) model with accurate and realistic mechanical properties, rheology and dynamics. The modeled membrane is represented by a triangular mesh which incorporates shear inplane energy, bending energy, and area and volume conservation constraints. The macroscopic membrane elastic properties are imposed through semi-analytic theory, and are matched with those obtained in optical tweezers stretching experiments. Rheological measurements characterized by time-dependent complex modulus are extracted from the membrane thermal fluctuations, and compared with those obtained from the optical magnetic twisting cytometry results. The results allow us to define a meaningful characteristic time of the membrane. The dynamics of RBCs observed in shear flow suggests that a purely elastic model for the RBC membrane is not appropriate, and therefore a viscoelastic model is required. The set of proposed analyses and numerical tests can be used as a complete model testbed in order to calibrate the modeled viscoelastic membranes to accurately represent RBCs in health and disease. PMID:19965026
Polymer Fluid Dynamics: Continuum and Molecular Approaches.
Bird, R B; Giacomin, A J
2016-06-01
To solve problems in polymer fluid dynamics, one needs the equations of continuity, motion, and energy. The last two equations contain the stress tensor and the heat-flux vector for the material. There are two ways to formulate the stress tensor: (a) One can write a continuum expression for the stress tensor in terms of kinematic tensors, or (b) one can select a molecular model that represents the polymer molecule and then develop an expression for the stress tensor from kinetic theory. The advantage of the kinetic theory approach is that one gets information about the relation between the molecular structure of the polymers and the rheological properties. We restrict the discussion primarily to the simplest stress tensor expressions or constitutive equations containing from two to four adjustable parameters, although we do indicate how these formulations may be extended to give more complicated expressions. We also explore how these simplest expressions are recovered as special cases of a more general framework, the Oldroyd 8-constant model. Studying the simplest models allows us to discover which types of empiricisms or molecular models seem to be worth investigating further. We also explore equivalences between continuum and molecular approaches. We restrict the discussion to several types of simple flows, such as shearing flows and extensional flows, which are of greatest importance in industrial operations. Furthermore, if these simple flows cannot be well described by continuum or molecular models, then it is not necessary to lavish time and energy to apply them to more complex flow problems. PMID:27276553
Molecular dynamics simulation of triclinic lysozyme in a crystal lattice.
Janowski, Pawel A; Liu, Chunmei; Deckman, Jason; Case, David A
2016-01-01
Molecular dynamics simulations of crystals can enlighten interpretation of experimental X-ray crystallography data and elucidate structural dynamics and heterogeneity in biomolecular crystals. Furthermore, because of the direct comparison against experimental data, they can inform assessment of molecular dynamics methods and force fields. We present microsecond scale results for triclinic hen egg-white lysozyme in a supercell consisting of 12 independent unit cells using four contemporary force fields (Amber ff99SB, ff14ipq, ff14SB, and CHARMM 36) in crystalline and solvated states (for ff14SB only). We find the crystal simulations consistent across multiple runs of the same force field and robust to various solvent equilibration schemes. However, convergence is slow compared with solvent simulations. All the tested force fields reproduce experimental structural and dynamic properties well, but Amber ff14SB maintains structure and reproduces fluctuations closest to the experimental model: its average backbone structure differs from the deposited structure by 0.37Å; by contrast, the average backbone structure in solution differs from the deposited by 0.65Å. All the simulations are affected by a small progressive deterioration of the crystal lattice, presumably due to imperfect modeling of hydrogen bonding and other crystal contact interactions; this artifact is smallest in ff14SB, with average lattice positions deviating by 0.20Å from ideal. Side-chain disorder is surprisingly low with fewer than 30% of the nonglycine or alanine residues exhibiting significantly populated alternate rotamers. Our results provide helpful insight into the methodology of biomolecular crystal simulations and indicate directions for future work to obtain more accurate energy models for molecular dynamics. PMID:26013419
Application of optimal prediction to molecular dynamics
Barber IV, John Letherman
2004-12-01
Optimal prediction is a general system reduction technique for large sets of differential equations. In this method, which was devised by Chorin, Hald, Kast, Kupferman, and Levy, a projection operator formalism is used to construct a smaller system of equations governing the dynamics of a subset of the original degrees of freedom. This reduced system consists of an effective Hamiltonian dynamics, augmented by an integral memory term and a random noise term. Molecular dynamics is a method for simulating large systems of interacting fluid particles. In this thesis, I construct a formalism for applying optimal prediction to molecular dynamics, producing reduced systems from which the properties of the original system can be recovered. These reduced systems require significantly less computational time than the original system. I initially consider first-order optimal prediction, in which the memory and noise terms are neglected. I construct a pair approximation to the renormalized potential, and ignore three-particle and higher interactions. This produces a reduced system that correctly reproduces static properties of the original system, such as energy and pressure, at low-to-moderate densities. However, it fails to capture dynamical quantities, such as autocorrelation functions. I next derive a short-memory approximation, in which the memory term is represented as a linear frictional force with configuration-dependent coefficients. This allows the use of a Fokker-Planck equation to show that, in this regime, the noise is {delta}-correlated in time. This linear friction model reproduces not only the static properties of the original system, but also the autocorrelation functions of dynamical variables.
Numerical methods for molecular dynamics. Progress report
Skeel, R.D.
1991-12-31
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
NASA Astrophysics Data System (ADS)
Tao, Jianmin; Rappe, Andrew M.
2016-01-01
Due to the absence of the long-range van der Waals (vdW) interaction, conventional density functional theory (DFT) often fails in the description of molecular complexes and solids. In recent years, considerable progress has been made in the development of the vdW correction. However, the vdW correction based on the leading-order coefficient C6 alone can only achieve limited accuracy, while accurate modeling of higher-order coefficients remains a formidable task, due to the strong non-additivity effect. Here, we apply a model dynamic multipole polarizability within a modified single-frequency approximation to calculate C8 and C10 between small molecules. We find that the higher-order vdW coefficients from this model can achieve remarkable accuracy, with mean absolute relative deviations of 5% for C8 and 7% for C10. Inclusion of accurate higher-order contributions in the vdW correction will effectively enhance the predictive power of DFT in condensed matter physics and quantum chemistry.
Symmetry Reduced Dynamics of Charged Molecular Strands
NASA Astrophysics Data System (ADS)
Ellis, David C. P.; Gay-Balmaz, François; Holm, Darryl D.; Putkaradze, Vakhtang; Ratiu, Tudor S.
2010-09-01
The equations of motion are derived for the dynamical folding of charged molecular strands (such as DNA) modeled as flexible continuous filamentary distributions of interacting rigid charge conformations. The new feature is that these equations are nonlocal when the screened Coulomb interactions, or Lennard-Jones potentials between pairs of charges, are included. The nonlocal dynamics is derived in the convective representation of continuum motion by using modified Euler-Poincaré and Hamilton-Pontryagin variational formulations that illuminate the various approaches within the framework of symmetry reduction of Hamilton’s principle for exact geometric rods. In the absence of nonlocal interactions, the equations recover the classical Kirchhoff theory of elastic rods. The motion equations in the convective representation are shown to arise by a classical Lagrangian reduction associated to the symmetry group of the system. This approach uses the process of affine Euler-Poincaré reduction initially developed for complex fluids. On the Hamiltonian side, the Poisson bracket of the molecular strand is obtained by reduction of the canonical symplectic structure on phase space. A change of variables allows a direct passage from this classical point of view to the covariant formulation in terms of Lagrange-Poincaré equations of field theory. In another revealing perspective, the convective representation of the nonlocal equations of molecular strand motion is transformed into quaternionic form.
Stochastic Event-Driven Molecular Dynamics
Donev, Aleksandar Garcia, Alejandro L.; Alder, Berni J.
2008-02-01
A novel Stochastic Event-Driven Molecular Dynamics (SEDMD) algorithm is developed for the simulation of polymer chains suspended in a solvent. SEDMD combines event-driven molecular dynamics (EDMD) with the Direct Simulation Monte Carlo (DSMC) method. The polymers are represented as chains of hard-spheres tethered by square wells and interact with the solvent particles with hard-core potentials. The algorithm uses EDMD for the simulation of the polymer chain and the interactions between the chain beads and the surrounding solvent particles. The interactions between the solvent particles themselves are not treated deterministically as in EDMD, rather, the momentum and energy exchange in the solvent is determined stochastically using DSMC. The coupling between the solvent and the solute is consistently represented at the particle level retaining hydrodynamic interactions and thermodynamic fluctuations. However, unlike full MD simulations of both the solvent and the solute, in SEDMD the spatial structure of the solvent is ignored. The SEDMD algorithm is described in detail and applied to the study of the dynamics of a polymer chain tethered to a hard-wall subjected to uniform shear. SEDMD closely reproduces results obtained using traditional EDMD simulations with two orders of magnitude greater efficiency. Results question the existence of periodic (cycling) motion of the polymer chain.
Electronic continuum model for molecular dynamics simulations.
Leontyev, I V; Stuchebrukhov, A A
2009-02-28
A simple model for accounting for electronic polarization in molecular dynamics (MD) simulations is discussed. In this model, called molecular dynamics electronic continuum (MDEC), the electronic polarization is treated explicitly in terms of the electronic continuum (EC) approximation, while the nuclear dynamics is described with a fixed-charge force field. In such a force-field all atomic charges are scaled to reflect the screening effect by the electronic continuum. The MDEC model is rather similar but not equivalent to the standard nonpolarizable force-fields; the differences are discussed. Of our particular interest is the calculation of the electrostatic part of solvation energy using standard nonpolarizable MD simulations. In a low-dielectric environment, such as protein, the standard MD approach produces qualitatively wrong results. The difficulty is in mistreatment of the electronic polarizability. We show how the results can be much improved using the MDEC approach. We also show how the dielectric constant of the medium obtained in a MD simulation with nonpolarizable force-field is related to the static (total) dielectric constant, which includes both the nuclear and electronic relaxation effects. Using the MDEC model, we discuss recent calculations of dielectric constants of alcohols and alkanes, and show that the MDEC results are comparable with those obtained with the polarizable Drude oscillator model. The applicability of the method to calculations of dielectric properties of proteins is discussed. PMID:19256627
Molecular dynamics studies of aromatic hydrocarbon liquids
McLaughlin, E.; Gupta, S.
1990-01-01
This project mainly involves a molecular dynamics and Monte Carlo study of the effect of molecular shape on thermophysical properties of bulk fluids with an emphasis on the aromatic hydrocarbon liquids. In this regard we have studied the modeling, simulation methodologies, and predictive and correlating methods for thermodynamic properties of fluids of nonspherical molecules. In connection with modeling we have studied the use of anisotropic site-site potentials, through a modification of the Gay-Berne Gaussian overlap potential, to successfully model the aromatic rings after adding the necessary electrostatic moments. We have also shown these interaction sites should be located at the geometric centers of the chemical groups. In connection with predictive methods, we have shown two perturbation type theories to work well for fluids modeled using one-center anisotropic potentials and the possibility exists for extending these to anisotropic site-site models. In connection with correlation methods, we have studied, through simulations, the effect of molecular shape on the attraction term in the generalized van der Waals equation of state for fluids of nonspherical molecules and proposed a possible form which is to be studied further. We have successfully studied the vector and parallel processing aspects of molecular simulations for fluids of nonspherical molecules.
Molecular dynamics in high electric fields
NASA Astrophysics Data System (ADS)
Apostol, M.; Cune, L. C.
2016-06-01
Molecular rotation spectra, generated by the coupling of the molecular electric-dipole moments to an external time-dependent electric field, are discussed in a few particular conditions which can be of some experimental interest. First, the spherical-pendulum molecular model is reviewed, with the aim of introducing an approximate method which consists in the separation of the azimuthal and zenithal motions. Second, rotation spectra are considered in the presence of a static electric field. Two particular cases are analyzed, corresponding to strong and weak fields. In both cases the classical motion of the dipoles consists of rotations and vibrations about equilibrium positions; this motion may exhibit parametric resonances. For strong fields a large macroscopic electric polarization may appear. This situation may be relevant for polar matter (like pyroelectrics, ferroelectrics), or for heavy impurities embedded in a polar solid. The dipolar interaction is analyzed in polar condensed matter, where it is shown that new polarization modes appear for a spontaneous macroscopic electric polarization (these modes are tentatively called "dipolons"); one of the polarization modes is related to parametric resonances. The extension of these considerations to magnetic dipoles is briefly discussed. The treatment is extended to strong electric fields which oscillate with a high frequency, as those provided by high-power lasers. It is shown that the effect of such fields on molecular dynamics is governed by a much weaker, effective, renormalized, static electric field.
NASA Astrophysics Data System (ADS)
Niklasson, Anders; Cawkwell, Marc
2012-02-01
Born-Oppenheimer molecular dynamics (BOMD) based on density functional theory offers a very accurate quantum mechanical approach to atomistic simulations that is more reliable and general compared to classical MD. Unfortunately, BOMD simulations are often limited by a high computational cost or by problems such as unbalanced phase space trajectories, numerical instabilities and a systematic long-term energy drift. These problems become particularly severe in combination with reduced complexity or linear scaling algorithms that are necessary for the study of large systems. We have recently taken some steps toward a new generation of first principles MD, which combines some of the best features of regular BOMD and Car-Parrinello MD, while avoiding their most serious shortcomings. The new dynamics is given in terms of an extended Lagrangian (XL), where auxiliary extended electronic degrees of freedom are added to the nuclear part. Our framework enables accurate geometric integration of both the nuclear and electronic degrees of freedom that provide a time-reversible and energy conserving dynamics on the ground state BO potential energy surface that is stable also under approximate SCF convergence. XL-BOMD provides a surprisingly simple and general framework for atomistic simulations
Thermostability of Enzymes from Molecular Dynamics Simulations.
Zeiske, Tim; Stafford, Kate A; Palmer, Arthur G
2016-06-14
Thermodynamic stability is a central requirement for protein function, and one goal of protein engineering is improvement of stability, particularly for applications in biotechnology. Herein, molecular dynamics simulations are used to predict in vitro thermostability of members of the bacterial ribonuclease HI (RNase H) family of endonucleases. The temperature dependence of the generalized order parameter, S, for four RNase H homologues, from psychrotrophic, mesophilic, and thermophilic organisms, is highly correlated with experimentally determined melting temperatures and with calculated free energies of folding at the midpoint temperature of the simulations. This study provides an approach for in silico mutational screens to improve thermostability of biologically and industrially relevant enzymes. PMID:27123810
8B structure in Fermionic Molecular Dynamics
NASA Astrophysics Data System (ADS)
Henninger, K. R.; Neff, T.; Feldmeier, H.
2015-04-01
The structure of the light exotic nucleus 8B is investigated in the Fermionic Molecular Dynamics (FMD) model. The decay of 8B is responsible for almost the entire high- energy solar-neutrino flux, making structure calculations of 8B important for determining the solar core temperature. 8B is a proton halo candidate thought to exhibit clustering. FMD uses a wave-packet basis and is well-suited for modelling clustering and halos. For a multiconfiguration treatment we construct the many-body Hilbert space from antisymmetrised angular-momentum projected 8-particle states. First results show formation of a proton halo.
Molecular dynamics simulations of dense plasmas
Collins, L.A.; Kress, J.D.; Kwon, I.; Lynch, D.L.; Troullier, N.
1993-12-31
We have performed quantum molecular dynamics simulations of hot, dense plasmas of hydrogen over a range of temperatures(0.1-5eV) and densities(0.0625-5g/cc). We determine the forces quantum mechanically from density functional, extended Huckel, and tight binding techniques and move the nuclei according to the classical equations of motion. We determine pair-correlation functions, diffusion coefficients, and electrical conductivities. We find that many-body effects predominate in this regime. We begin to obtain agreement with the OCP and Thomas-Fermi models only at the higher temperatures and densities.
Molecular beam studies of reaction dynamics
Lee, Y.T.
1987-03-01
Purpose of this research project is two-fold: (1) to elucidate detailed dynamics of simple elementary reactions which are theoretically important and to unravel the mechanism of complex chemical reactions or photo chemical processes which play an important role in many macroscopic processes and (2) to determine the energetics of polyatomic free radicals using microscopic experimental methods. Most of the information is derived from measurement of the product fragment translational energy and angular distributions using unique molecular beam apparati designed for these purposes.
Molecular dynamics simulation of ice XII
NASA Astrophysics Data System (ADS)
Borzsák, István; Cummings, Peter T.
1999-02-01
Molecular dynamics simulations have been performed on the newly discovered metastable ice XII. This new crystalline ice phase [C. Lobban, J.L. Finney, W.F. Kuhs, Nature (London) 391 (1998) 268] is proton-disordered. Thus 90 possible configurations of the unit cell can be constructed which differ only in the orientations of the water molecules. The simulation used the TIP4P potential model for water at constant temperature and density. About one-quarter of the initial configurations did not melt in the course of the simulation. This result is supportive of the experimental structure and also demonstrates the ability of this water model to study ice phases.
Crystallization of nickel nanoclusters by molecular dynamics
NASA Astrophysics Data System (ADS)
Chamati, H.; Gaminchev, K.
2012-12-01
We investigated the melting properties of bulk nickel and the crystallization of nickel nanocrystals via molecular dynamics using a potential in the framework of the second moment approximation of tight-binding theory. The melting behavior was simulated with the hysteresis approach by subsequently heating and cooling gradually the system over a wide range of temperatures. The crystallization of nickel nanoclusters consisting of 55, 147 and 309 atoms was achieved after repeatedly annealing and quenching the corresponding quasicrystals several times to avoid being trapped in a local energy minimum. The time over which the global minimum was reached was found to increase with the cluster size.
Molecular dynamics at constant temperature and pressure
NASA Astrophysics Data System (ADS)
Toxvaerd, S.
1993-01-01
Algorithms for molecular dynamics (MD) at constant temperature and pressure are investigated. The ability to remain in a regular orbit in an intermittent chaotic regime is used as a criterion for long-time stability. A simple time-centered algorithm (leap frog) is found to be the most stable of the commonly used algorithms in MD. A model of N one-dimensional dimers with a double-well intermolecular potential, for which the distribution functions at constant temperature T and pressure P can be calculated, is used to investigate MD-NPT dynamics. A time-centered NPT algorithm is found to sample correctly and to be very robust with respect to volume scaling.
Charge transport network dynamics in molecular aggregates.
Jackson, Nicholas E; Chen, Lin X; Ratner, Mark A
2016-08-01
Due to the nonperiodic nature of charge transport in disordered systems, generating insight into static charge transport networks, as well as analyzing the network dynamics, can be challenging. Here, we apply time-dependent network analysis to scrutinize the charge transport networks of two representative molecular semiconductors: a rigid n-type molecule, perylenediimide, and a flexible p-type molecule, [Formula: see text] Simulations reveal the relevant timescale for local transfer integral decorrelation to be [Formula: see text]100 fs, which is shown to be faster than that of a crystalline morphology of the same molecule. Using a simple graph metric, global network changes are observed over timescales competitive with charge carrier lifetimes. These insights demonstrate that static charge transport networks are qualitatively inadequate, whereas average networks often overestimate network connectivity. Finally, a simple methodology for tracking dynamic charge transport properties is proposed. PMID:27439871
Accurate rotor loads prediction using the FLAP (Force and Loads Analysis Program) dynamics code
Wright, A.D.; Thresher, R.W.
1987-10-01
Accurately predicting wind turbine blade loads and response is very important in predicting the fatigue life of wind turbines. There is a clear need in the wind turbine community for validated and user-friendly structural dynamics codes for predicting blade loads and response. At the Solar Energy Research Institute (SERI), a Force and Loads Analysis Program (FLAP) has been refined and validated and is ready for general use. Currently, FLAP is operational on an IBM-PC compatible computer and can be used to analyze both rigid- and teetering-hub configurations. The results of this paper show that FLAP can be used to accurately predict the deterministic loads for rigid-hub rotors. This paper compares analytical predictions to field test measurements for a three-bladed, upwind turbine with a rigid-hub configuration. The deterministic loads predicted by FLAP are compared with 10-min azimuth averages of blade root flapwise bending moments for different wind speeds. 6 refs., 12 figs., 3 tabs.
Combining molecular dynamics with mesoscopic Green's function reaction dynamics simulations
NASA Astrophysics Data System (ADS)
Vijaykumar, Adithya; Bolhuis, Peter G.; ten Wolde, Pieter Rein
2015-12-01
In many reaction-diffusion processes, ranging from biochemical networks, catalysis, to complex self-assembly, the spatial distribution of the reactants and the stochastic character of their interactions are crucial for the macroscopic behavior. The recently developed mesoscopic Green's Function Reaction Dynamics (GFRD) method enables efficient simulation at the particle level provided the microscopic dynamics can be integrated out. Yet, many processes exhibit non-trivial microscopic dynamics that can qualitatively change the macroscopic behavior, calling for an atomistic, microscopic description. We propose a novel approach that combines GFRD for simulating the system at the mesoscopic scale where particles are far apart, with a microscopic technique such as Langevin dynamics or Molecular Dynamics (MD), for simulating the system at the microscopic scale where reactants are in close proximity. This scheme defines the regions where the particles are close together and simulated with high microscopic resolution and those where they are far apart and simulated with lower mesoscopic resolution, adaptively on the fly. The new multi-scale scheme, called MD-GFRD, is generic and can be used to efficiently simulate reaction-diffusion systems at the particle level.
Molecular dynamics simulation of liquid sulfur dioxide.
Ribeiro, Mauro C C
2006-05-01
A previously proposed model for molecular dynamics (MD) simulation of liquid sulfur dioxide, SO(2), has been reviewed. Thermodynamic, structural, and dynamical properties were calculated for a large range of thermodynamic states. Predicted (P,V,T) of simulated system agrees with an elaborated equation of state recently proposed for liquid SO(2). Calculated heat capacity, expansion coefficient, and isothermal compressibility are also in good agreement with experimental data. Calculated equilibrium structure agrees with X-ray and neutron scattering measurements on liquid SO(2). The model also predicts the same (SO(2))(2) dimer structure as previously determined by ab initio calculations. Detailed analysis of equilibrium structure of liquid SO(2) is provided, indicating that, despite the rather large dipole moment of the SO(2) molecule, the structure is mainly determined by the Lennard-Jones interactions. Both single-particle and collective dynamics are investigated. Temperature dependency of dynamical properties is given. The MD results are compared with previous findings obtained from the analysis of inelastic neutron scattering spectra of liquid SO(2), including wave-vector dependent structural relaxation, tau(k), and viscosity, eta(k). PMID:16640437
Frank, Martin; Gutbrod, Peter; Hassayoun, Chokri; von Der Lieth, Claus-W
2003-10-01
Molecular dynamics is a rapidly developing field of science and has become an established tool for studying the dynamic behavior of biomolecules. Although several high quality programs for performing molecular dynamic simulations are freely available, only well-trained scientists are currently able to make use of the broad scientific potential that molecular dynamic simulations offer to gain insight into structural questions at an atomic level. The "Dynamic Molecules" approach is the first internet portal that provides an interactive access to set up, perform and analyze molecular dynamic simulations. It is completely based on standard web technologies and uses only publicly available software. The aim is to open molecular dynamics techniques to a broader range of users including undergraduate students, teachers and scientists outside the bioinformatics field. The time-limiting factors are the availability of free capacity on the computing server to run the simulations and the time required to transport the history file through the internet for the animation mode. The interactive access mode of the portal is acceptable for animations of molecules having up to about 500 atoms. PMID:12908101
Steered Molecular Dynamics Methods Applied to Enzyme Mechanism and Energetics.
Ramírez, C L; Martí, M A; Roitberg, A E
2016-01-01
One of the main goals of chemistry is to understand the underlying principles of chemical reactions, in terms of both its reaction mechanism and the thermodynamics that govern it. Using hybrid quantum mechanics/molecular mechanics (QM/MM)-based methods in combination with a biased sampling scheme, it is possible to simulate chemical reactions occurring inside complex environments such as an enzyme, or aqueous solution, and determining the corresponding free energy profile, which provides direct comparison with experimental determined kinetic and equilibrium parameters. Among the most promising biasing schemes is the multiple steered molecular dynamics method, which in combination with Jarzynski's Relationship (JR) allows obtaining the equilibrium free energy profile, from a finite set of nonequilibrium reactive trajectories by exponentially averaging the individual work profiles. However, obtaining statistically converged and accurate profiles is far from easy and may result in increased computational cost if the selected steering speed and number of trajectories are inappropriately chosen. In this small review, using the extensively studied chorismate to prephenate conversion reaction, we first present a systematic study of how key parameters such as pulling speed, number of trajectories, and reaction progress are related to the resulting work distributions and in turn the accuracy of the free energy obtained with JR. Second, and in the context of QM/MM strategies, we introduce the Hybrid Differential Relaxation Algorithm, and show how it allows obtaining more accurate free energy profiles using faster pulling speeds and smaller number of trajectories and thus smaller computational cost. PMID:27497165
Dynamic Shear Modulus of Polymers from Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Byutner, Oleksiy; Smith, Grant
2001-03-01
In this work we describe the methodology for using equilibrium molecular dynamics simulations (MD) simulations to obtain the viscoelastic properties of polymers in the glassy regime. Specifically we show how the time dependent shear stress modulus and frequency dependent complex shear modulus in the high-frequency regime can be determined from the off-diagonal terms of the stress-tensor autocorrelation function obtained from MD trajectories using the Green-Kubo method and appropriate Fourier transforms. In order to test the methodology we have performed MD simulations of a low-molecular-weight polybutadiene system using quantum chemistry based potential functions. Values of the glassy modulus and the maximum loss frequency were found to be in good agreement with experimental data for polybutadiene at 298 K.
An Efficient Time-Stepping Scheme for Ab Initio Molecular Dynamics Simulations
NASA Astrophysics Data System (ADS)
Tsuchida, Eiji
2016-08-01
In ab initio molecular dynamics simulations of real-world problems, the simple Verlet method is still widely used for integrating the equations of motion, while more efficient algorithms are routinely used in classical molecular dynamics. We show that if the Verlet method is used in conjunction with pre- and postprocessing, the accuracy of the time integration is significantly improved with only a small computational overhead. We also propose several extensions of the algorithm required for use in ab initio molecular dynamics. The validity of the processed Verlet method is demonstrated in several examples including ab initio molecular dynamics simulations of liquid water. The structural properties obtained from the processed Verlet method are found to be sufficiently accurate even for large time steps close to the stability limit. This approach results in a 2× performance gain over the standard Verlet method for a given accuracy. We also show how to generate a canonical ensemble within this approach.
Exact dynamic properties of molecular motors
NASA Astrophysics Data System (ADS)
Boon, N. J.; Hoyle, R. B.
2012-08-01
Molecular motors play important roles within a biological cell, performing functions such as intracellular transport and gene transcription. Recent experimental work suggests that there are many plausible biochemical mechanisms that molecules such as myosin-V could use to achieve motion. To account for the abundance of possible discrete-stochastic frameworks that can arise when modeling molecular motor walks, a generalized and straightforward graphical method for calculating their dynamic properties is presented. It allows the calculation of the velocity, dispersion, and randomness ratio for any proposed system through analysis of its structure. This article extends work of King and Altman ["A schematic method of deriving the rate laws of enzyme-catalyzed reactions," J. Phys. Chem. 60, 1375-1378 (1956)], 10.1021/j150544a010 on networks of enzymatic reactions by calculating additional dynamic properties for spatially hopping systems. Results for n-state systems are presented: single chain, parallel pathway, divided pathway, and divided pathway with a chain. A novel technique for combining multiple system architectures coupled at a reference state is also demonstrated. Four-state examples illustrate the effectiveness and simplicity of these methods.
Structure and Dynamics of Cellulose Molecular Solutions
NASA Astrophysics Data System (ADS)
Wang, Howard; Zhang, Xin; Tyagi, Madhusudan; Mao, Yimin; Briber, Robert
Molecular dissolution of microcrystalline cellulose has been achieved through mixing with ionic liquid 1-Ethyl-3-methylimidazolium acetate (EMIMAc), and organic solvent dimethylformamide (DMF). The mechanism of cellulose dissolution in tertiary mixtures has been investigated by combining quasielastic and small angle neutron scattering (QENS and SANS). As SANS data show that cellulose chains take Gaussian-like conformations in homogenous solutions, which exhibit characteristics of having an upper critical solution temperature, the dynamic signals predominantly from EMIMAc molecules indicate strong association with cellulose in the dissolution state. The mean square displacement quantities support the observation of the stoichiometric 3:1 EMIMAc to cellulose unit molar ratio, which is a necessary criterion for the molecular dissolution of cellulose. Analyses of dynamics structure factors reveal the temperature dependence of a slow and a fast process for EMIMAc's bound to cellulose and in DMF, respectively, as well as a very fast process due possibly to the rotational motion of methyl groups, which persisted to near the absolute zero.
Exact dynamic properties of molecular motors.
Boon, N J; Hoyle, R B
2012-08-28
Molecular motors play important roles within a biological cell, performing functions such as intracellular transport and gene transcription. Recent experimental work suggests that there are many plausible biochemical mechanisms that molecules such as myosin-V could use to achieve motion. To account for the abundance of possible discrete-stochastic frameworks that can arise when modeling molecular motor walks, a generalized and straightforward graphical method for calculating their dynamic properties is presented. It allows the calculation of the velocity, dispersion, and randomness ratio for any proposed system through analysis of its structure. This article extends work of King and Altman ["A schematic method of deriving the rate laws of enzyme-catalyzed reactions," J. Phys. Chem. 60, 1375-1378 (1956)] on networks of enzymatic reactions by calculating additional dynamic properties for spatially hopping systems. Results for n-state systems are presented: single chain, parallel pathway, divided pathway, and divided pathway with a chain. A novel technique for combining multiple system architectures coupled at a reference state is also demonstrated. Four-state examples illustrate the effectiveness and simplicity of these methods. PMID:22938213
Gaussian Molecular Dynamics in Imaginary and Real Time
NASA Astrophysics Data System (ADS)
Georgescu, Ionut; Mandelshtam, Vladimir
2010-03-01
The variational Gaussian wavepacket (VGW) method can be used to estimate the equilibrium density matrix by propagating Gaussian wavepackets in imaginary time [1,2]. It has proven to be practically accurate and computationally less expensive than the path integral methods. We compare the VGW method to the Feynman-Kleinert approximation (FKA), which has comparable computational cost. Although both methods are variational, they utilize different variational principles: In FKA the partition function is optimized, while in VGW it is the imaginary-time-dependent wave packet. We show that the VGW method is more accurate for a wide variety of systems. The differences are particularly important when thermodynamic properties, such as heat capacity, are of main interest. Moreover, unlike the case of FKA, in the VGW method the imaginary frequencies do not arise. In the spirit of the Centroid Molecular Dynamics the VGW method has also been extended to simulate the real-time dynamics, e.g., it can be used to estimate the Kubo-transformed quantum time correlation functions. The latter are exact in the high-temperature and harmonic limits. [1] P. Frantsuzov and V.A. Mandelshtam, J. Chem. Phys 121, 9247 (2004)[2] C. Predescu, P. Frantsuzov and V.A. Mandelshtam J. Chem. Phys 122, 154305 (2005)
Earthquake Rupture Dynamics using Adaptive Mesh Refinement and High-Order Accurate Numerical Methods
NASA Astrophysics Data System (ADS)
Kozdon, J. E.; Wilcox, L.
2013-12-01
Our goal is to develop scalable and adaptive (spatial and temporal) numerical methods for coupled, multiphysics problems using high-order accurate numerical methods. To do so, we are developing an opensource, parallel library known as bfam (available at http://bfam.in). The first application to be developed on top of bfam is an earthquake rupture dynamics solver using high-order discontinuous Galerkin methods and summation-by-parts finite difference methods. In earthquake rupture dynamics, wave propagation in the Earth's crust is coupled to frictional sliding on fault interfaces. This coupling is two-way, required the simultaneous simulation of both processes. The use of laboratory-measured friction parameters requires near-fault resolution that is 4-5 orders of magnitude higher than that needed to resolve the frequencies of interest in the volume. This, along with earlier simulations using a low-order, finite volume based adaptive mesh refinement framework, suggest that adaptive mesh refinement is ideally suited for this problem. The use of high-order methods is motivated by the high level of resolution required off the fault in earlier the low-order finite volume simulations; we believe this need for resolution is a result of the excessive numerical dissipation of low-order methods. In bfam spatial adaptivity is handled using the p4est library and temporal adaptivity will be accomplished through local time stepping. In this presentation we will present the guiding principles behind the library as well as verification of code against the Southern California Earthquake Center dynamic rupture code validation test problems.
Dynamic transitions in molecular dynamics simulations of supercooled silicon
NASA Astrophysics Data System (ADS)
Mei, Xiaojun; Eapen, Jacob
2013-04-01
Two dynamic transitions or crossovers, one at a low temperature (T* ≈ 1006 K) and the other at a high temperature (T0 ≈ 1384 K), are shown to emerge in supercooled liquid silicon using molecular dynamics simulations. The high-temperature transition (T0) marks the decoupling of stress, density, and energy relaxation mechanisms. At the low-temperature transition (T*), depending on the cooling rate, supercooled silicon can either undergo a high-density-liquid to low-density-liquid (HDL-LDL) phase transition or experience an HDL-HDL crossover. Dynamically heterogeneous domains that emerge with supercooling become prominent across the HDL-HDL transition at 1006 K, with well-separated mobile and immobile regions. Interestingly, across the HDL-LDL transition, the most mobile atoms form large prominent aggregates while the least mobile atoms get spatially dispersed akin to that in a crystalline state. The attendant partial return to spatial uniformity with the HDL-LDL phase transition indicates a dynamic mechanism for relieving the frustration in supercooled states.
Atomistic molecular dynamic simulations of multiferroics.
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L
2012-08-10
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings. PMID:23006300
Atomistic Molecular Dynamic Simulations of Multiferroics
NASA Astrophysics Data System (ADS)
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L.
2012-08-01
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings.
The Relation between Approximation in Distribution and Shadowing in Molecular Dynamics
NASA Astrophysics Data System (ADS)
Tupper, Paul
2009-01-01
Molecular dynamics refers to the computer simulation of a material at the atomic level. An open problem in numerical analysis is to explain the apparent reliability of molecular dynamics simulations. The difficulty is that individual trajectories computed in molecular dynamics are accurate for only short time intervals, whereas apparently reliable information can be extracted from very long-time simulations. It has been conjectured that long molecular dynamics trajectories have low-dimensional statistical features that accurately approximate those of the original system. Another conjecture is that numerical trajectories satisfy the shadowing property: they are close over long time intervals to exact trajectories but with different initial conditions. We prove that these two views are actually equivalent to each other, after we suitably modify the concept of shadowing. A key ingredient of our result is a general theorem that allows us to take random elements of a metric space that are close in distribution and embed them in the same probability space so that they are close in a strong sense. This result is similar to the Strassen-Dudley theorem except that a mapping is provided between the two random elements. Our results on shadowing are motivated by molecular dynamics but apply to the approximation of any dynamical system when initial conditions are selected according to a probability measure.
Role of Molecular Dynamics and Related Methods in Drug Discovery.
De Vivo, Marco; Masetti, Matteo; Bottegoni, Giovanni; Cavalli, Andrea
2016-05-12
Molecular dynamics (MD) and related methods are close to becoming routine computational tools for drug discovery. Their main advantage is in explicitly treating structural flexibility and entropic effects. This allows a more accurate estimate of the thermodynamics and kinetics associated with drug-target recognition and binding, as better algorithms and hardware architectures increase their use. Here, we review the theoretical background of MD and enhanced sampling methods, focusing on free-energy perturbation, metadynamics, steered MD, and other methods most consistently used to study drug-target binding. We discuss unbiased MD simulations that nowadays allow the observation of unsupervised ligand-target binding, assessing how these approaches help optimizing target affinity and drug residence time toward improved drug efficacy. Further issues discussed include allosteric modulation and the role of water molecules in ligand binding and optimization. We conclude by calling for more prospective studies to attest to these methods' utility in discovering novel drug candidates. PMID:26807648
Molecular dynamics simulations on the melting of gold nanoparticles
NASA Astrophysics Data System (ADS)
Qiao, Zhiwei; Feng, Haijun; Zhou, Jian
2014-01-01
Molecular dynamics is employed to study the melting of bulk gold and gold nanoparticles. PCFF, Sutton-Chen and COMPASS force fields are adopted to study the melting point of bulk gold and we find out that the Sutton-Chen force field is the most accurate model in predicting the melting point of bulk gold. Consequently, the Sutton-Chen force field is applied to study the melting points of spherical gold nanoparticles with different diameters. Variations of diffusion coefficient, potential energy and translational order parameter with temperature are analyzed. The simulated melting points of gold nanoparticles are between 615∼1115 K, which are much lower than that of bulk gold (1336 K). As the diameter of gold nanoparticle drops, the melting point also descends. The melting mechanism is also analyzed for gold nanoparticles.
Molecular dynamics of shock loading of metals with defects
Belak, J.F.
1997-12-31
The finite rise time of shock waves in metals is commonly attributed to dissipative or viscous behavior of the metal. This viscous or plastic behavior is commonly attributed to the motion of defects such as dislocations. Despite this intuitive understanding, the experimental observation of defect motion or nucleation during shock loading has not been possible due to the short time scales involved. Molecular dynamics modeling with realistic interatomic potentials can provide some insight into defect motion during shock loading. However, until quite recently, the length scale required to accurately represent a metal with defects has been beyond the scope of even the most powerful supercomputers. Here, the author presents simulations of the shock response of single defects and indicate how simulation might provide some insight into the shock loading of metals.
Fast and accurate modeling of molecular atomization energies with machine learning.
Rupp, Matthias; Tkatchenko, Alexandre; Müller, Klaus-Robert; von Lilienfeld, O Anatole
2012-02-01
We introduce a machine learning model to predict atomization energies of a diverse set of organic molecules, based on nuclear charges and atomic positions only. The problem of solving the molecular Schrödinger equation is mapped onto a nonlinear statistical regression problem of reduced complexity. Regression models are trained on and compared to atomization energies computed with hybrid density-functional theory. Cross validation over more than seven thousand organic molecules yields a mean absolute error of ∼10 kcal/mol. Applicability is demonstrated for the prediction of molecular atomization potential energy curves. PMID:22400967
NASA Astrophysics Data System (ADS)
Xin, Cui; Di-Yu, Zhang; Gao, Chen; Ji-Gen, Chen; Si-Liang, Zeng; Fu-Ming, Guo; Yu-Jun, Yang
2016-03-01
We demonstrate that the interference minima in the linear molecular harmonic spectra can be accurately predicted by a modified two-center model. Based on systematically investigating the interference minima in the linear molecular harmonic spectra by the strong-field approximation (SFA), it is found that the locations of the harmonic minima are related not only to the nuclear distance between the two main atoms contributing to the harmonic generation, but also to the symmetry of the molecular orbital. Therefore, we modify the initial phase difference between the double wave sources in the two-center model, and predict the harmonic minimum positions consistent with those simulated by SFA. Project supported by the National Basic Research Program of China (Grant No. 2013CB922200) and the National Natural Science Foundation of China (Grant Nos. 11274001, 11274141, 11304116, 11247024, and 11034003), and the Jilin Provincial Research Foundation for Basic Research, China (Grant Nos. 20130101012JC and 20140101168JC).
Osmosis : a molecular dynamics computer simulation study
NASA Astrophysics Data System (ADS)
Lion, Thomas
Osmosis is a phenomenon of critical importance in a variety of processes ranging from the transport of ions across cell membranes and the regulation of blood salt levels by the kidneys to the desalination of water and the production of clean energy using potential osmotic power plants. However, despite its importance and over one hundred years of study, there is an ongoing confusion concerning the nature of the microscopic dynamics of the solvent particles in their transfer across the membrane. In this thesis the microscopic dynamical processes underlying osmotic pressure and concentration gradients are investigated using molecular dynamics (MD) simulations. I first present a new derivation for the local pressure that can be used for determining osmotic pressure gradients. Using this result, the steady-state osmotic pressure is studied in a minimal model for an osmotic system and the steady-state density gradients are explained using a simple mechanistic hopping model for the solvent particles. The simulation setup is then modified, allowing us to explore the timescales involved in the relaxation dynamics of the system in the period preceding the steady state. Further consideration is also given to the relative roles of diffusive and non-diffusive solvent transport in this period. Finally, in a novel modification to the classic osmosis experiment, the solute particles are driven out-of-equilibrium by the input of energy. The effect of this modification on the osmotic pressure and the osmotic ow is studied and we find that active solute particles can cause reverse osmosis to occur. The possibility of defining a new "osmotic effective temperature" is also considered and compared to the results of diffusive and kinetic temperatures..
Molecular dynamics simulations of solvated yeast tRNA(Asp).
Auffinger, P; Louise-May, S; Westhof, E
1999-01-01
Transfer RNA molecules are involved in a variety of biological processes, implying complex recognition events with proteins and other RNAs. From a structural point of view, tRNAs constitute a reference system for studying RNA folding and architecture. A deeper understanding of their structural and functional properties will derive from our ability to model accurately their dynamical behavior. We present the first dynamical model of a fully neutralized and solvated tRNA molecule over a 500-ps time scale. Starting from the crystallographic structure of yeast tRNA(Asp), the 75-nucleotide molecule was modeled with 8055 water molecules and 74 NH4+ counterions, using the AMBER4.1 program and the particle mesh Ewald (PME) method for the treatment of long-range electrostatic interactions. The calculations led to a dynamically stable model of the tRNA molecule. During the simulation, all secondary and tertiary base pairs are maintained while a certain lability of base triples in the tRNA core is observed. This lability was interpreted as resulting from intrinsic factors associated with the "weaker" hydrogen bonding patterns seen in these base triples and from an altered ionic environment of the tRNA molecule. Calculated thermal factors are used to compare the dynamics of the tRNA in solution and in the crystal. The present molecular dynamics simulation of a complex and highly charged nucleic acid molecule attests to the fact that simulation methods are now able to investigate not only the dynamics of proteins, but also that of large RNA molecules. Thus they also provide a basis for further investigations on the structural and functional effects of chemical and posttranscriptionally modified nucleotides as well as on ionic environmental effects. PMID:9876122
Fiber lubrication: A molecular dynamics simulation study
NASA Astrophysics Data System (ADS)
Liu, Hongyi
Molecular and mesoscopic level description of friction and lubrication remains a challenge because of difficulties in the phenomenological understanding of to the behaviors of solid-liquid interfaces during sliding. Fortunately, there is the computational simulation approach opens an opportunity to predict and analyze interfacial phenomena, which were studied with molecular dynamics (MD) and mesoscopic dynamics (MesoDyn) simulations. Polypropylene (PP) and cellulose are two of most common polymers in textile fibers. Confined amorphous surface layers of PP and cellulose were built successfully with xenon crystals which were used to compact the polymers. The physical and surface properties of the PP and cellulose surface layers were investigated by MD simulations, including the density, cohesive energy, volumetric thermal expansion, and contact angle with water. The topology method was employed to predict the properties of poly(alkylene glycol) (PAG) diblock copolymers and Pluronic triblock copolymers used as lubricants on surfaces. Density, zero shear viscosity, shear module, cohesive energy and solubility parameter were predicted with each block copolymer. Molecular dynamics simulations were used to study the interaction energy per unit contact area of block copolymer melts with PP and cellulose surfaces. The interaction energy is defined as the ratio of interfacial interaction energy to the contact area. Both poly(proplene oxide) (PPO) and poly(ethylene oxide) (PEO) segments provided a lipophilic character to both PP and cellulose surfaces. The PPO/PEO ratio and the molecular weight were found to impact the interaction energy on both PP and cellulose surfaces. In aqueous solutions, the interaction energy is complicated due to the presence of water and the cross interactions between the multiple molecular components. The polymer-water-surface (PWS) calculation method was proposed to calculate such complex systems. In a contrast with a vacuum condition, the presence
Development of semiclassical molecular dynamics simulation method.
Nakamura, Hiroki; Nanbu, Shinkoh; Teranishi, Yoshiaki; Ohta, Ayumi
2016-04-28
Various quantum mechanical effects such as nonadiabatic transitions, quantum mechanical tunneling and coherence play crucial roles in a variety of chemical and biological systems. In this paper, we propose a method to incorporate tunneling effects into the molecular dynamics (MD) method, which is purely based on classical mechanics. Caustics, which define the boundary between classically allowed and forbidden regions, are detected along classical trajectories and the optimal tunneling path with minimum action is determined by starting from each appropriate caustic. The real phase associated with tunneling can also be estimated. Numerical demonstration with use of a simple collinear chemical reaction O + HCl → OH + Cl is presented in order to help the reader to well comprehend the method proposed here. Generalization to the on-the-fly ab initio version is rather straightforward. By treating the nonadiabatic transitions at conical intersections by the Zhu-Nakamura theory, new semiclassical MD methods can be developed. PMID:27067383
Molecular Dynamics Simulation of Shock Induced Detonation
NASA Astrophysics Data System (ADS)
Tomar, Vikas; Zhou, Min
2004-07-01
This research focuses on molecular dynamics (MD) simulation of shock induced detonation in Fe2O3+Al thermite mixtures. A MD model is developed to simulate non-equilibrium stress-induced reactions. The focus is on establishing a criterion for reaction initiation, energy content and rate of energy release as functions of mixture and reinforcement characteristics. A cluster functional potential is proposed for this purpose. The potential uses the electronegativity equalization to account for changes in the charge of different species according to local environment. Parameters in the potential are derived to fit to the properties of Fe, Al, Fe2O3, and Al2O3. NPT MD simulations are carried out to qualitatively check the energetics of the forward (Fe2O3+Al) as well as backward (Al2O3+Fe) thermite reactions. The results show that the potential can account for the energetics of thermite reactions.
Assessing Molecular Dynamics Simulations with Solvatochromism Modeling.
Schwabe, Tobias
2015-08-20
For the modeling of solvatochromism with an explicit representation of the solvent molecules, the quality of preceding molecular dynamics simulations is crucial. Therefore, the possibility to apply force fields which are derived with as little empiricism as possible seems desirable. Such an approach is tested here by exploiting the sensitive solvatochromism of p-nitroaniline, and the use of reliable excitation energies based on approximate second-order coupled cluster results within a polarizable embedding scheme. The quality of the various MD settings for four different solvents, water, methanol, ethanol, and dichloromethane, is assessed. In general, good agreement with the experiment is observed when polarizable force fields and special treatment of hydrogen bonding are applied. PMID:26220273
Nonequilibrium molecular dynamics: The first 25 years
Hoover, W.G. |
1992-08-01
Equilibrium Molecular Dynamics has been generalized to simulate Nonequilibrium systems by adding sources of thermodynamic heat and work. This generalization incorporates microscopic mechanical definitions of macroscopic thermodynamic and hydrodynamic variables, such as temperature and stress, and augments atomistic forces with special boundary, constraint, and driving forces capable of doing work on, and exchanging heat with, an otherwise Newtonian system. The underlying Lyapunov instability of these nonequilibrium equations of motion links microscopic time-reversible deterministic trajectories to macroscopic time-irreversible hydrodynamic behavior as described by the Second Law of Thermodynamics. Green-Kubo linear-response theory has been checked. Nonlinear plastic deformation, intense heat conduction, shockwave propagation, and nonequilibrium phase transformation have all been simulated. The nonequilibrium techniques, coupled with qualitative improvements in parallel computer hardware, are enabling simulations to approximate real-world microscale and nanoscale experiments.
Molecular Dynamics Simulations of Water Evaporation
NASA Astrophysics Data System (ADS)
Wen, Chengyuan; Grest, Gary; Cheng, Shengfeng
2015-03-01
The evaporation of water from the liquid/vapor interface is studied via large-scale molecular dynamics simulations for systems of more than a million atoms at 550K and 600K. The TIP4P-2005 water model whose liquid/vapor surface tension is in excellent agreement with experiments is used. Evaporative cooling at the interface is observed from temperature profiles determined from both translational and rotational kinetic energy. During evaporation, the density of water is slightly enhanced near the liquid-vapor interface. The velocity distribution of water molecules in the vapor phase during evaporation at various distances relative to the interface fit a Maxwell-Boltzmann distribution. While our results indicate an imbalance between evaporating and condensing water molecules, local thermal equilibrium is found to hold in addition to mechanical equilibrium. Department of Physics, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA.
Cluster production within antisymmetrized molecular dynamics
NASA Astrophysics Data System (ADS)
Ono, Akira
2016-06-01
Clusters are quite important at various situations in heavy-ion collisions. Antisymmetrized molecular dynamics was improved to take into account the correlations to form light clusters, such as deuterons and α particles, and light nuclei composed of several clusters. The momentum fluctuations of emitted particles are also taken into account by a simple method. Formation of fragments and light clusters in a wide range of heavy-ion collisions was well described with a single set of model parameters. Fragmentation in a proton induced reaction was also well reproduced by introducing cluster correlations. Calculated results demonstrate strong impacts of clusters in various observables including those usually regarded as probes of the density dependence of symmetry energy.
Cell list algorithms for nonequilibrium molecular dynamics
NASA Astrophysics Data System (ADS)
Dobson, Matthew; Fox, Ian; Saracino, Alexandra
2016-06-01
We present two modifications of the standard cell list algorithm that handle molecular dynamics simulations with deforming periodic geometry. Such geometry naturally arises in the simulation of homogeneous, linear nonequilibrium flow modeled with periodic boundary conditions, and recent progress has been made developing boundary conditions suitable for general 3D flows of this type. Previous works focused on the planar flows handled by Lees-Edwards or Kraynik-Reinelt boundary conditions, while the new versions of the cell list algorithm presented here are formulated to handle the general 3D deforming simulation geometry. As in the case of equilibrium, for short-ranged pairwise interactions, the cell list algorithm reduces the computational complexity of the force computation from O(N2) to O(N), where N is the total number of particles in the simulation box. We include a comparison of the complexity and efficiency of the two proposed modifications of the standard algorithm.
Molecular-dynamics simulations of lead clusters
NASA Astrophysics Data System (ADS)
Hendy, S. C.; Hall, B. D.
2001-08-01
Molecular-dynamics simulations of nanometer-sized lead clusters have been performed using the Lim-Ong-Ercolessi glue potential [Surf. Sci. 269/270, 1109 (1992)]. The binding energies of clusters forming crystalline (fcc), decahedron and icosahedron structures are compared, showing that fcc cuboctahedra are the most energetically favored of these polyhedral model structures. However, simulations of the freezing of liquid droplets produced a characteristic form of surface-reconstructed ``shaved'' icosahedron, in which atoms are absent at the edges and apexes of the polyhedron. This arrangement is energetically favored for 600-4000 atom clusters. Larger clusters favor crystalline structures. Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect fcc Wulff particle, containing a number of parallel stacking faults. The effects of temperature on the preferred structure of crystalline clusters below the melting point have been considered. The implications of these results for the interpretation of experimental data is discussed.
Dielectrophoresis of nanocolloids: a molecular dynamics study.
Salonen, E; Terama, E; Vattulainen, I; Karttunen, M
2005-10-01
Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer-sized colloidal particles. We introduce a simplified model for a polarizable nanoparticle, consisting of a large charged macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the particle at different combinations of temperature and electric field strength. In accord with linear response theory, the particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the mobility further results in an apparent weak variation of the DEP displacements with temperature. PMID:16195818
Dynamics of dewetting at the nanoscale using molecular dynamics.
Bertrand, E; Blake, T D; Ledauphin, V; Ogonowski, G; Coninck, J De; Fornasiero, D; Ralston, J
2007-03-27
Large-scale molecular dynamics simulations are used to model the dewetting of solid surfaces by partially wetting thin liquid films. Two levels of solid-liquid interaction are considered that give rise to large equilibrium contact angles. The initial length and thickness of the films are varied over a wide range at the nanoscale. Spontaneous dewetting is initiated by removing a band of molecules either from each end of the film or from its center. As observed experimentally and in previous simulations, the films recede at an initially constant speed, creating a growing rim of liquid with a constant receding dynamic contact angle. Consistent with the current understanding of wetting dynamics, film recession is faster on the more poorly wetted surface to an extent that cannot be explained solely by the increase in the surface tension driving force. In addition, the rates of recession of the thinnest films are found to increase with decreasing film thickness. These new results imply not only that the mobility of the liquid molecules adjacent to the solid increases with decreasing solid-liquid interactions, but also that the mobility adjacent to the free surface of the film is higher than in the bulk, so that the effective viscosity of the film decreases with thickness. PMID:17328565
Molecular dynamics simulations of microscale fluid transport
Wong, C.C.; Lopez, A.R.; Stevens, M.J.; Plimpton, S.J.
1998-02-01
Recent advances in micro-science and technology, like Micro-Electro-Mechanical Systems (MEMS), have generated a group of unique liquid flow problems that involve characteristic length scales of a Micron. Also, in manufacturing processes such as coatings, current continuum models are unable to predict microscale physical phenomena that appear in these non-equilibrium systems. It is suspected that in these systems, molecular-level processes can control the interfacial energy and viscoelastic properties at the liquid/solid boundary. A massively parallel molecular dynamics (MD) code has been developed to better understand microscale transport mechanisms, fluid-structure interactions, and scale effects in micro-domains. Specifically, this MD code has been used to analyze liquid channel flow problems for a variety of channel widths, e.g. 0.005-0.05 microns. This report presents results from MD simulations of Poiseuille flow and Couette flow problems and addresses both scaling and modeling issues. For Poiseuille flow, the numerical predictions are compared with existing data to investigate the variation of the friction factor with channel width. For Couette flow, the numerical predictions are used to determine the degree of slip at the liquid/solid boundary. Finally, the results also indicate that shear direction with respect to the wall lattice orientation can be very important. Simulation results of microscale Couette flow and microscale Poiseuille flow for two different surface structures and two different shear directions will be presented.
Molecular dynamics simulations of supramolecular polymer rheology
NASA Astrophysics Data System (ADS)
Li, Zhenlong; Djohari, Hadrian; Dormidontova, Elena E.
2010-11-01
Using equilibrium and nonequilibrium molecular dynamics simulations, we studied the equilibrium and rheological properties of dilute and semidilute solutions of head-to-tail associating polymers. In our simulation model, a spontaneous complementary reversible association between the donor and the acceptor groups at the ends of oligomers was achieved by introducing a combination of truncated pseudo-Coulombic attractive potential and Lennard Jones repulsive potential between donor, acceptor, and neighboring groups. We have calculated the equilibrium properties of supramolecular polymers, such as the ring/chain equilibrium, average molecular weight, and molecular weight distribution of self-assembled chains and rings, which all agree well with previous analytical and computer modeling results. We have investigated shear thinning of solutions of 8- and 20-bead associating oligomers with different association energies at different temperatures and oligomer volume fractions. All reduced viscosity data for a given oligomer length can be collapsed into one master curve, exhibiting two power-law regions of shear-thinning behavior with an exponent of -0.55 at intermediate ranges of the reduced shear rate β and -0.8 (or -0.9) at larger shear rates. The equilibrium viscosity of supramolecular solutions with different oligomer lengths and associating energies is found to obey a power-law scaling dependence on oligomer volume fraction with an exponent of 1.5, in agreement with the experimental observations for several dilute or semidilute solutions of supramolecular polymers. This implies that dilute and semidilute supramolecular polymer solutions exhibit high polydispersity but may not be sufficiently entangled to follow the reptation mechanism of relaxation.
Nguyen, Thuy-Diem; Schmidt, Bertil; Zheng, Zejun; Kwoh, Chee-Keong
2015-01-01
De novo clustering is a popular technique to perform taxonomic profiling of a microbial community by grouping 16S rRNA amplicon reads into operational taxonomic units (OTUs). In this work, we introduce a new dendrogram-based OTU clustering pipeline called CRiSPy. The key idea used in CRiSPy to improve clustering accuracy is the application of an anomaly detection technique to obtain a dynamic distance cutoff instead of using the de facto value of 97 percent sequence similarity as in most existing OTU clustering pipelines. This technique works by detecting an abrupt change in the merging heights of a dendrogram. To produce the output dendrograms, CRiSPy employs the OTU hierarchical clustering approach that is computed on a genetic distance matrix derived from an all-against-all read comparison by pairwise sequence alignment. However, most existing dendrogram-based tools have difficulty processing datasets larger than 10,000 unique reads due to high computational complexity. We address this difficulty by developing two efficient algorithms for CRiSPy: a compute-efficient GPU-accelerated parallel algorithm for pairwise distance matrix computation and a memory-efficient hierarchical clustering algorithm. Our experiments on various datasets with distinct attributes show that CRiSPy is able to produce more accurate OTU groupings than most OTU clustering applications. PMID:26451819
Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method
NASA Technical Reports Server (NTRS)
Blech, Richard A.
2002-01-01
One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.
Molecular beam studies of reaction dynamics
Lee, Y.T.
1993-12-01
The major thrust of this research project is to elucidate detailed dynamics of simple elementary reactions that are theoretically important and to unravel the mechanism of complex chemical reactions or photochemical processes that play important roles in many macroscopic processes. Molecular beams of reactants are used to study individual reactive encounters between molecules or to monitor photodissociation events in a collision-free environment. Most of the information is derived from measurement of the product fragment energy, angular, and state distributions. Recent activities are centered on the mechanisms of elementary chemical reactions involving oxygen atoms with unsaturated hydrocarbons, the dynamics of endothermic substitution reactions, the dependence of the chemical reactivity of electronically excited atoms on the alignment of excited orbitals, the primary photochemical processes of polyatomic molecules, intramolecular energy transfer of chemically activated and locally excited molecules, the energetics of free radicals that are important to combustion processes, the infrared-absorption spectra of carbonium ions and hydrated hydronium ions, and bond-selective photodissociation through electric excitation.
Molecular beam studies of reaction dynamics
Lee, Yuan T.
1991-03-01
The major thrust of this research project is to elucidate detailed dynamics of simple elementary reactions that are theoretically important and to unravel the mechanism of complex chemical reactions or photochemical processes that play important roles in many macroscopic processes. Molecular beams of reactants are used to study individual reactive encounters between molecules or to monitor photodissociation events in a collision-free environment. Most of the information is derived from measurement of the product fragment energy, angular, and state distributions. Recent activities are centered on the mechanisms of elementary chemical reactions involving oxygen atoms with unsaturated hydrocarbons, the dynamics of endothermic substitution reactions, the dependence of the chemical reactivity of electronically excited atoms on the alignment of excited orbitals, the primary photochemical processes of polyatomic molecules, intramolecular energy transfer of chemically activated and locally excited molecules, the energetics of free radicals that are important to combustion processes, the infrared-absorption spectra of carbonium ions and hydrated hydronium ions, and bond-selective photodissociation through electric excitation.
Accelerated electronic structure-based molecular dynamics simulations of shock-induced chemistry
NASA Astrophysics Data System (ADS)
Cawkwell, Marc
2015-06-01
The initiation and progression of shock-induced chemistry in organic materials at moderate temperatures and pressures are slow on the time scales available to regular molecular dynamics simulations. Accessing the requisite time scales is particularly challenging if the interatomic bonding is modeled using accurate yet expensive methods based explicitly on electronic structure. We have combined fast, energy conserving extended Lagrangian Born-Oppenheimer molecular dynamics with the parallel replica accelerated molecular dynamics formalism to study the relatively sluggish shock-induced chemistry of benzene around 13-20 GPa. We model interatomic bonding in hydrocarbons using self-consistent tight binding theory with an accurate and transferable parameterization. Shock compression and its associated transient, non-equilibrium effects are captured explicitly by combining the universal liquid Hugoniot with a simple shrinking-cell boundary condition. A number of novel methods for improving the performance of reactive electronic structure-based molecular dynamics by adapting the self-consistent field procedure on-the-fly will also be discussed. The use of accelerated molecular dynamics has enabled us to follow the initial stages of the nucleation and growth of carbon clusters in benzene under thermodynamic conditions pertinent to experiments.
Molecular dynamics studies on nanoscale gas transport
NASA Astrophysics Data System (ADS)
Barisik, Murat
Three-dimensional molecular dynamics (MD) simulations of nanoscale gas flows are studied to reveal surface effects. A smart wall model that drastically reduces the memory requirements of MD simulations for gas flows is introduced. The smart wall molecular dynamics (SWMD) represents three-dimensional FCC walls using only 74 wall Molecules. This structure is kept in the memory and utilized for each gas molecule surface collision. Using SWMD, fluid behavior within nano-scale confinements is studied for argon in dilute gas, dense gas, and liquid states. Equilibrium MD method is employed to resolve the density and stress variations within the static fluid. Normal stress calculations are based on the Irving-Kirkwood method, which divides the stress tensor into its kinetic and virial parts. The kinetic component recovers pressure based on the ideal gas law. The particle-particle virial increases with increased density, while the surface-particle virial develops due to the surface force field effects. Normal stresses within nano-scale confinements show anisotropy induced primarily by the surface force-field and local variations in the fluid density near the surfaces. For dilute and dense gas cases, surface-force field that extends typically 1nm from each wall induces anisotropic normal stress. For liquid case, this effect is further amplified by the density fluctuations that extend beyond the three field penetration region. Outside the wall force-field penetration and density fluctuation regions the normal stress becomes isotropic and recovers the thermodynamic pressure, provided that sufficiently large force cut-off distances are utilized in the computations. Next, non-equilibrium SWMD is utilized to investigate the surface-gas interaction effects on nanoscale shear-driven gas flows in the transition and free molecular flow regimes. For the specified surface properties and gas-surface pair interactions, density and stress profiles exhibit a universal behavior inside the
Długosz, Maciej; Antosiewicz, Jan M
2015-07-01
Proper treatment of hydrodynamic interactions is of importance in evaluation of rigid-body mobility tensors of biomolecules in Stokes flow and in simulations of their folding and solution conformation, as well as in simulations of the translational and rotational dynamics of either flexible or rigid molecules in biological systems at low Reynolds numbers. With macromolecules conveniently modeled in calculations or in dynamic simulations as ensembles of spherical frictional elements, various approximations to hydrodynamic interactions, such as the two-body, far-field Rotne-Prager approach, are commonly used, either without concern or as a compromise between the accuracy and the numerical complexity. Strikingly, even though the analytical Rotne-Prager approach fails to describe (both in the qualitative and quantitative sense) mobilities in the simplest system consisting of two spheres, when the distance between their surfaces is of the order of their size, it is commonly applied to model hydrodynamic effects in macromolecular systems. Here, we closely investigate hydrodynamic effects in two and three-body systems, consisting of bead-shell molecular models, using either the analytical Rotne-Prager approach, or an accurate numerical scheme that correctly accounts for the many-body character of hydrodynamic interactions and their short-range behavior. We analyze mobilities, and translational and rotational velocities of bodies resulting from direct forces acting on them. We show, that with the sufficient number of frictional elements in hydrodynamic models of interacting bodies, the far-field approximation is able to provide a description of hydrodynamic effects that is in a reasonable qualitative as well as quantitative agreement with the description resulting from the application of the virtually exact numerical scheme, even for small separations between bodies. PMID:26068580
Molecular dynamics simulation of amorphous indomethacin.
Xiang, Tian-Xiang; Anderson, Bradley D
2013-01-01
Molecular dynamics (MD) simulations have been conducted using an assembly consisting of 105 indomethacin (IMC) molecules and 12 water molecules to investigate the underlying dynamic (e.g., rotational and translational diffusivities and conformation relaxation rates) and structural properties (e.g., conformation, hydrogen-bonding distributions, and interactions of water with IMC) of amorphous IMC. These properties may be important in predicting physical stability of this metastable material. The IMC model was constructed using X-ray diffraction data with the force-field parameters mostly assigned by analogy with similar groups in Amber-ff03 and atomic charges calculated with the B3LYP/ccpVTZ30, IEFPCM, and RESP models. The assemblies were initially equilibrated in their molten state and cooled through the glass transition temperature to form amorphous solids. Constant temperature dynamic runs were then carried out above and below the T(g) (i.e., at 600 K (10 ns), 400 K (350 ns), and 298 K (240 ns)). The density (1.312 ± 0.003 g/cm(3)) of the simulated amorphous solid at 298 K was close to the experimental value (1.32 g/cm(3)) while the estimated T(g) (384 K) was ~64 degrees higher than the experimental value (320 K) due to the faster cooling rate. Due to the hindered rotation of its amide bond, IMC can exist in different diastereomeric states. Different IMC conformations were sufficiently sampled in the IMC melt or vapor, but transitions occurred rarely in the glass. The hydrogen-bonding patterns in amorphous IMC are more complex in the amorphous state than in the crystalline polymorphs. Carboxylic dimers that are dominant in α- and γ-crystals were found to occur at a much lower probability in the simulated IMC glasses while hydrogen-bonded IMC chains were more easily identified patterns in the simulated amorphous solids. To determine molecular diffusivity, a novel analytical method is proposed to deal with the non-Einsteinian behavior, in which the temporal
GAS-PHASE MOLECULAR DYNAMICS: VIBRATIONAL DYNAMICS OF POLYATOMIC MOLECULES
MUCKERMAN,J.T.
1999-06-09
The goal of this research is the understanding of elementary chemical and physical processes important in the combustion of fossil fuels. Interest centers on reactions and properties of short-lived chemical intermediates. High-resolution, high-sensitivity, laser absorption methods are augmented by high-temperature, flow-tube reaction kinetics studies with mass-spectrometric sampling. These experiments provide information on the energy levels, structures and reactivity of molecular free radical species and, in turn, provide new tools for the study of energy flow and chemical bond cleavage in radicals involved in chemical systems. The experimental work is supported by theoretical studies using time-dependent quantum wavepacket calculations, which provide insight into energy flow among the vibrational modes of polyatomic molecules and interference effects in multiple-surface dynamics.
Gas-Phase Molecular Dynamics: Vibrational Dynamics of Polyatomic Molecules
Muckerman, J.T.
1999-05-21
The goal of this research is the understanding of elementary chemical and physical processes important in the combustion of fossil fuels. Interest centers on reactions and properties of short-lived chemical intermediates. High-resolution, high-sensitivity, laser absorption methods are augmented by high- temperature, flow-tube reaction kinetics studies with mass-spectrometic sampling. These experiments provide information on the energy levels, structures and reactivity of molecular free radical species and in turn, provide new tools for the study of energy flow and chemical bond cleavage in the radicals involved in chemical systems. The experimental work is supported by theoretical studies using time-dependent quantum wavepacket calculations, which provide insight into energy flow among the vibrational modes of polyatomic molecules and interference effects in multiple-surface dynamics.
Internal Coordinate Molecular Dynamics: A Foundation for Multiscale Dynamics
2015-01-01
Internal coordinates such as bond lengths, bond angles, and torsion angles (BAT) are natural coordinates for describing a bonded molecular system. However, the molecular dynamics (MD) simulation methods that are widely used for proteins, DNA, and polymers are based on Cartesian coordinates owing to the mathematical simplicity of the equations of motion. However, constraints are often needed with Cartesian MD simulations to enhance the conformational sampling. This makes the equations of motion in the Cartesian coordinates differential-algebraic, which adversely impacts the complexity and the robustness of the simulations. On the other hand, constraints can be easily placed in BAT coordinates by removing the degrees of freedom that need to be constrained. Thus, the internal coordinate MD (ICMD) offers an attractive alternative to Cartesian coordinate MD for developing multiscale MD method. The torsional MD method is a special adaptation of the ICMD method, where all the bond lengths and bond angles are kept rigid. The advantages of ICMD simulation methods are the longer time step size afforded by freezing high frequency degrees of freedom and performing a conformational search in the more important low frequency torsional degrees of freedom. However, the advancements in the ICMD simulations have been slow and stifled by long-standing mathematical bottlenecks. In this review, we summarize the recent mathematical advancements we have made based on spatial operator algebra, in developing a robust long time scale ICMD simulation toolkit useful for various applications. We also present the applications of ICMD simulations to study conformational changes in proteins and protein structure refinement. We review the advantages of the ICMD simulations over the Cartesian simulations when used with enhanced sampling methods and project the future use of ICMD simulations in protein dynamics. PMID:25517406
Tosso, Rodrigo D; Andujar, Sebastian A; Gutierrez, Lucas; Angelina, Emilio; Rodríguez, Ricaurte; Nogueras, Manuel; Baldoni, Héctor; Suvire, Fernando D; Cobo, Justo; Enriz, Ricardo D
2013-08-26
A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper. PMID:23834278
Molecular Detection of Foodborne Pathogens: A Rapid and Accurate Answer to Food Safety.
Mangal, Manisha; Bansal, Sangita; Sharma, Satish K; Gupta, Ram K
2016-07-01
Food safety is a global health concern. For the prevention and recognition of problems related to health and safety, detection of foodborne pathogen is of utmost importance at all levels of food production chain. For several decades, a lot of research has been targeted at the development of rapid methodology as reducing the time needed to complete pathogen detection tests has been the primary goal of food microbiologists. With the result, food microbiology laboratories now have a wide array of detection methods and automated technologies such as enzyme immunoassay, polymerase chain reaction, and microarrays, which can cut test times considerably. Nucleic acid amplification strategies and advances in amplicon detection methodologies have been the key factors in the progress of molecular microbiology. A comprehensive literature survey has been carried out to give an overview in the field of foodborne pathogen detection. In this paper, we describe the conventional methods, as well as recent developments in food pathogen detection, identification, and quantification, with a major emphasis on molecular detection methods. PMID:25830555
In situ structure and dynamics of DNA origami determined through molecular dynamics simulations
Yoo, Jejoong; Aksimentiev, Aleksei
2013-01-01
The DNA origami method permits folding of long single-stranded DNA into complex 3D structures with subnanometer precision. Transmission electron microscopy, atomic force microscopy, and recently cryo-EM tomography have been used to characterize the properties of such DNA origami objects, however their microscopic structures and dynamics have remained unknown. Here, we report the results of all-atom molecular dynamics simulations that characterized the structural and mechanical properties of DNA origami objects in unprecedented microscopic detail. When simulated in an aqueous environment, the structures of DNA origami objects depart from their idealized targets as a result of steric, electrostatic, and solvent-mediated forces. Whereas the global structural features of such relaxed conformations conform to the target designs, local deformations are abundant and vary in magnitude along the structures. In contrast to their free-solution conformation, the Holliday junctions in the DNA origami structures adopt a left-handed antiparallel conformation. We find the DNA origami structures undergo considerable temporal fluctuations on both local and global scales. Analysis of such structural fluctuations reveals the local mechanical properties of the DNA origami objects. The lattice type of the structures considerably affects global mechanical properties such as bending rigidity. Our study demonstrates the potential of all-atom molecular dynamics simulations to play a considerable role in future development of the DNA origami field by providing accurate, quantitative assessment of local and global structural and mechanical properties of DNA origami objects. PMID:24277840
An ab initio molecular dynamics study of the roaming mechanism of the H2+HOC+ reaction
NASA Astrophysics Data System (ADS)
Yu, Hua-Gen
2011-08-01
We report here a direct ab initio molecular dynamics study of the p-/o-H2+HOC+ reaction on the basis of the accurate SAC-MP2 potential energy surface. The quasi-classical trajectory method was employed. This work largely focuses on the study of reaction mechanisms. A roaming mechanism was identified for this molecular ion-molecule reaction. The driving forces behind the roaming mechanism were thoroughly investigated by using a trajectory dynamics approach. In addition, the thermal rate coefficients of the H2+HOC+ reaction were calculated in the temperature range [25, 300] K and are in good agreement with experiments.
Restoring electronic coherence/decoherence for a trajectory-based nonadiabatic molecular dynamics
NASA Astrophysics Data System (ADS)
Zhu, Chaoyuan
2016-04-01
By utilizing the time-independent semiclassical phase integral, we obtained modified coupled time-dependent Schrödinger equations that restore coherences and induce decoherences within original simple trajectory-based nonadiabatic molecular dynamic algorithms. Nonadiabatic transition probabilities simulated from both Tully’s fewest switches and semiclassical Ehrenfest algorithms follow exact quantum electronic oscillations and amplitudes for three out of the four well-known model systems. Within the present theory, nonadiabatic transitions estimated from statistical ensemble of trajectories accurately follow those of the modified electronic wave functions. The present theory can be immediately applied to the molecular dynamic simulations of photochemical and photophysical processes involving electronic excited states.