Simple posterior frequency correction for vibrational spectra from molecular dynamics
NASA Astrophysics Data System (ADS)
Tikhonov, Denis S.
2016-05-01
Vibrational spectra computed from molecular dynamics simulations with large integration time steps suffer from nonphysical frequency shifts of signals [M. Praprotnik and D. Janežič, J. Chem. Phys. 122, 174103 (2005)]. A simple posterior correction technique was developed for compensation of this behavior. It performs through replacement of abscissa in the calculated spectra using following formula: ν corrected = /√{ 2 ṡ (" separators=" 1 - cos ( 2 π ṡ Δ t ṡ ν initial ) ) } 2 π ṡ Δ t , where ν are initial and corrected frequencies and Δt is the MD simulation time step. Applicability of this method was tested on gaseous infrared spectra of hydrogen fluoride and formic acid.
Simple posterior frequency correction for vibrational spectra from molecular dynamics.
Tikhonov, Denis S
2016-05-01
Vibrational spectra computed from molecular dynamics simulations with large integration time steps suffer from nonphysical frequency shifts of signals [M. Praprotnik and D. Janežič, J. Chem. Phys. 122, 174103 (2005)]. A simple posterior correction technique was developed for compensation of this behavior. It performs through replacement of abscissa in the calculated spectra using following formula: νcorrected=2⋅1-cos(2π⋅Δt⋅νinitial)2π⋅Δt, where ν are initial and corrected frequencies and Δt is the MD simulation time step. Applicability of this method was tested on gaseous infrared spectra of hydrogen fluoride and formic acid. PMID:27155626
NASA Astrophysics Data System (ADS)
Costandy, Joseph; Michalis, Vasileios K.; Tsimpanogiannis, Ioannis N.; Stubos, Athanassios K.; Economou, Ioannis G.
2016-03-01
We introduce a simple correction to the calculation of the lattice constants of fully occupied structure sI methane or carbon dioxide pure hydrates that are obtained from classical molecular dynamics simulations using the TIP4PQ/2005 water force field. The obtained corrected lattice constants are subsequently used in order to obtain isobaric thermal expansion coefficients of the pure gas hydrates that exhibit a trend that is significantly closer to the experimental behavior than previously reported classical molecular dynamics studies.
Costandy, Joseph; Michalis, Vasileios K; Tsimpanogiannis, Ioannis N; Stubos, Athanassios K; Economou, Ioannis G
2016-03-28
We introduce a simple correction to the calculation of the lattice constants of fully occupied structure sI methane or carbon dioxide pure hydrates that are obtained from classical molecular dynamics simulations using the TIP4PQ/2005 water force field. The obtained corrected lattice constants are subsequently used in order to obtain isobaric thermal expansion coefficients of the pure gas hydrates that exhibit a trend that is significantly closer to the experimental behavior than previously reported classical molecular dynamics studies. PMID:27036466
NASA Astrophysics Data System (ADS)
Ma, Zhonghua; Zhang, Yanli; Tuckerman, Mark E.
2012-07-01
It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)], 10.1002/jcc.20495. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92 g/cm3 while that obtained with dispersion corrections is 1.07 g/cm3, indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.
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.
Molecular dynamics of large systems with quantum corrections for the nuclei
NASA Astrophysics Data System (ADS)
Gu, Bing; Garashchuk, Sophya
2015-12-01
This paper describes an approximate approach to quantum dynamics based on the quantum trajectory formulation of the Schrödinger equation. The quantum-mechanical effects are incorporated through the quantum potential of the mean-field type, acting on a trajectory ensemble in addition to the classical potential. Efficiency for large systems is achieved by using the quantum corrections for selected degrees of freedom and introduction of empirical friction into the ground-state energy calculations. The classical potential, if needed, can be computed on-the-fly using the Density Functional Tight Binding method of electronic structure merged with the quantum trajectory dynamics code. The approach is practical for a few hundred atoms. Applications include a study of adsorption of quantum hydrogen colliding with the graphene model, C37H15 and a calculation of the ground state of solid 4He simulated by a cell 180-atoms.
Molecular dynamics of large systems with quantum corrections for the nuclei
Gu, Bing; Garashchuk, Sophya
2015-12-31
This paper describes an approximate approach to quantum dynamics based on the quantum trajectory formulation of the Schrödinger equation. The quantum-mechanical effects are incorporated through the quantum potential of the mean-field type, acting on a trajectory ensemble in addition to the classical potential. Efficiency for large systems is achieved by using the quantum corrections for selected degrees of freedom and introduction of empirical friction into the ground-state energy calculations. The classical potential, if needed, can be computed on-the-fly using the Density Functional Tight Binding method of electronic structure merged with the quantum trajectory dynamics code. The approach is practical for a few hundred atoms. Applications include a study of adsorption of quantum hydrogen colliding with the graphene model, C{sub 37}H{sub 15} and a calculation of the ground state of solid {sup 4}He simulated by a cell 180-atoms.
VandeVondele, Joost; Sprik, Michiel
2005-04-01
We have performed density functional theory based molecular dynamics (MD) simulations of the *OH radical in solution using self-interaction corrected (SIC) methods. We use a scheme recently proposed by M. d'Avezac, M. Calandra and F. Mauri [arXiv:cond-mat/0407750] in which a correction is only applied to the spin density within a restricted open shell formulation. In addition to two correction formulas employed within this scheme by M. d'Avezac, M. Calandra and F. Mauri, we propose and test an new empirical form which only introduces a scaled Coulomb term. This new functional leads to good agreement with reference calculations on radical cation dimers and on the hydroxyl water dimer in the gas phase. Applied in ab initio MD simulations, these three SIC methods provide a picture of the *OH solvation that differs qualitatively from the one obtained using the standard generalised gradient approximation (GGA). Hemibonded water, observed in GGA simulations and believed to be an artefact due to self-interaction error, is not present. We find that the *OH acts as a good hydrogen bond donor, but accepts less than two hydrogen bonds on average. These hydrogen bonds are part of a mobile, otherwise quasi-hydrophobic solvation cage. Our results show the potential of this computationally expedient scheme, which might extend the range of problems that can be modelled adequately with density functional theory. PMID:19787955
Correcting for the free energy costs of bond or angle constraints in molecular dynamics simulations
König, Gerhard; Brooks, Bernard R.
2014-01-01
Background Free energy simulations are an important tool in the arsenal of computational biophysics, allowing the calculation of thermodynamic properties of binding or enzymatic reactions. This paper introduces methods to increase the accuracy and precision of free energy calculations by calculating the free energy costs of constraints during post-processing. The primary purpose of employing constraints for these free energy methods is to increase the phase space overlap between ensembles, which is required for accuracy and convergence. Methods The free energy costs of applying or removing constraints are calculated as additional explicit steps in the free energy cycle. The new techniques focus on hard degrees of freedom and use both gradients and Hessian estimation. Enthalpy, vibrational entropy, and Jacobian free energy terms are considered. Results We demonstrate the utility of this method with simple classical systems involving harmonic and anharmonic oscillators, four-atomic benchmark systems, an alchemical mutation of ethane to methanol, and free energy simulations between alanine and serine. The errors for the analytical test cases are all below 0.0007 kcal/mol, and the accuracy of the free energy results of ethane to methanol is improved from 0.15 to 0.04 kcal/mol. For the alanine to serine case, the phase space overlaps of the unconstrained simulations range between 0.15 and 0.9%. The introduction of constraints increases the overlap up to 2.05%. On average, the overlap increases by 94% relative to the unconstrained value and precision is doubled. Conclusions The approach reduces errors arising from constraints by about an order of magnitude. Free energy simulations benefit from the use of constraints through enhanced convergence and higher precision. General Significance The primary utility of this approach is to calculate free energies for systems with disparate energy surfaces and bonded terms, especially in multi-scale molecular mechanics
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.
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.
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.
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.
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
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.
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.
Quantum corrections to inflaton and curvaton dynamics
Markkanen, Tommi; Tranberg, Anders E-mail: anders.tranberg@nbi.dk
2012-11-01
We compute the fully renormalized one-loop effective action for two interacting and self-interacting scalar fields in FRW space-time. We then derive and solve the quantum corrected equations of motion both for fields that dominate the energy density (such as an inflaton) and fields that do not (such as a subdominant curvaton). In particular, we introduce quantum corrected Friedmann equations that determine the evolution of the scale factor. We find that in general, gravitational corrections are negligible for the field dynamics. For the curvaton-type fields this leaves only the effect of the flat-space Coleman-Weinberg-type effective potential, and we find that these can be significant. For the inflaton case, both the corrections to the potential and the Friedmann equations can lead to behaviour very different from the classical evolution. Even to the point that inflation, although present at tree level, can be absent at one-loop order.
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 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
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
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 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)
Self-correcting maps of molecular pathways.
Rzhetsky, Andrey; Zheng, Tian; Weinreb, Chani
2006-01-01
Reliable and comprehensive maps of molecular pathways are indispensable for guiding complex biomedical experiments. Such maps are typically assembled from myriads of disparate research reports and are replete with inconsistencies due to variations in experimental conditions and/or errors. It is often an intractable task to manually verify internal consistency over a large collection of experimental statements. To automate large-scale reconciliation efforts, we propose a random-arcs-and-nodes model where both nodes (tissue-specific states of biological molecules) and arcs (interactions between them) are represented with random variables. We show how to obtain a non-contradictory model of a molecular network by computing the joint distribution for arc and node variables, and then apply our methodology to a realistic network, generating a set of experimentally testable hypotheses. This network, derived from an automated analysis of over 3,000 full-text research articles, includes genes that have been hypothetically linked to four neurological disorders: Alzheimer's disease, autism, bipolar disorder, and schizophrenia. We estimated that approximately 10% of the published molecular interactions are logically incompatible. Our approach can be directly applied to an array of diverse problems including those encountered in molecular biology, ecology, economics, politics, and sociology. PMID:17183692
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%.
Stochastic Dynamics with Correct Sampling for Constrained Systems.
Peters, E A J F; Goga, N; Berendsen, H J C
2014-10-14
In this paper we discuss thermostatting using stochastic methods for molecular simulations where constraints are present. For so-called impulsive thermostats, like the Andersen thermostat, the equilibrium temperature will differ significantly from the imposed temperature when a limited number of particles are picked and constraints are applied. We analyze this problem and give two rigorous solutions for it. A correct general treatment of impulsive stochastic thermostatting, including pairwise dissipative particle dynamics and stochastic forcing in the presence of constraints, is given and it is shown that the constrained canonical distribution is sampled rigorously. We discuss implementation issues such as second order Trotter expansions. The method is shown to rigorously maintain the correct temperature for the case of extended simple point charge (SPC/E) water simulations. PMID:26588119
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
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
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.
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.
Small, Matthew W; Sanchez, Sergio I; Menard, Laurent D; Kang, Joo H; Frenkel, Anatoly I; Nuzzo, Ralph G
2011-03-16
This study describes a prototypical, bimetallic heterogeneous catalyst: compositionally well-defined Ir-Pt nanoclusters with sizes in the range of 1-2 nm supported on γ-Al(2)O(3). Deposition of the molecular bimetallic cluster [Ir(3)Pt(3)(μ-CO)(3)(CO)(3)(η-C(5)Me(5))(3)] on γ-Al(2)O(3), and its subsequent reduction with hydrogen, provides highly dispersed supported bimetallic Ir-Pt nanoparticles. Using spherical aberration-corrected scanning transmission electron microscopy (C(s)-STEM) and theoretical modeling of synchrotron-based X-ray absorption spectroscopy (XAS) measurements, our studies provide unambiguous structural assignments for this model catalytic system. The atomic resolution C(s)-STEM images reveal strong and specific lattice-directed strains in the clusters that follow local bonding configurations of the γ-Al(2)O(3) support. Combined nanobeam diffraction (NBD) and high-resolution transmission electron microscopy (HRTEM) data suggest the polycrystalline γ-Al(2)O(3) support material predominantly exposes (001) and (011) surface planes (ones commensurate with the zone axis orientations frequently exhibited by the bimetallic clusters). The data reveal that the supported bimetallic clusters exhibit complex patterns of structural dynamics, ones evidencing perturbations of an underlying oblate/hemispherical cuboctahedral cluster-core geometry with cores that are enriched in Ir (a result consistent with models based on surface energetics, which favor an ambient cluster termination by Pt) due to the dynamical responses of the M-M bonding to the specifics of the adsorbate and metal-support interactions. Taken together, the data demonstrate that strong temperature-dependent charge-transfer effects occur that are likely mediated variably by the cluster-support, cluster-adsorbate, and intermetallic bonding interactions. PMID:21341654
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 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.
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.
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
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.
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
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
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
Proposal on dynamic correction method for resonance ionization mass spectrometry
NASA Astrophysics Data System (ADS)
Noto, Takuma; Tomita, Hideki; Richter, Sven; Schneider, Fabian; Wendt, Klaus; Iguchi, Tetsuo; Kawarabayashi, Jun
2013-04-01
For high precision and accuracy in isotopic ratio measurement of transuranic elements using laser ablation assisted resonance ionization mass spectrometry, a dynamic correction method based on correlation of ion signals with energy and timing of each laser pulse was proposed. The feasibility of this dynamic correction method was investigated through the use of a programmable electronics device for fast acquisition of the energy and timing of each laser pulse.
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)
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.
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
Software for Correcting the Dynamic Error of Force Transducers
Miyashita, Naoki; Watanabe, Kazuhide; Irisa, Kyouhei; Iwashita, Hiroshi; Araki, Ryosuke; Takita, Akihiro; Yamaguchi, Takao; Fujii, Yusaku
2014-01-01
Software which corrects the dynamic error of force transducers in impact force measurements using their own output signal has been developed. The software corrects the output waveform of the transducers using the output waveform itself, estimates its uncertainty and displays the results. In the experiment, the dynamic error of three transducers of the same model are evaluated using the Levitation Mass Method (LMM), in which the impact forces applied to the transducers are accurately determined as the inertial force of the moving part of the aerostatic linear bearing. The parameters for correcting the dynamic error are determined from the results of one set of impact measurements of one transducer. Then, the validity of the obtained parameters is evaluated using the results of the other sets of measurements of all the three transducers. The uncertainties in the uncorrected force and those in the corrected force are also estimated. If manufacturers determine the correction parameters for each model using the proposed method, and provide the software with the parameters corresponding to each model, then users can obtain the waveform corrected against dynamic error and its uncertainty. The present status and the future prospects of the developed software are discussed in this paper. PMID:25004158
Software for correcting the dynamic error of force transducers.
Miyashita, Naoki; Watanabe, Kazuhide; Irisa, Kyouhei; Iwashita, Hiroshi; Araki, Ryosuke; Takita, Akihiro; Yamaguchi, Takao; Fujii, Yusaku
2014-01-01
Software which corrects the dynamic error of force transducers in impact force measurements using their own output signal has been developed. The software corrects the output waveform of the transducers using the output waveform itself, estimates its uncertainty and displays the results. In the experiment, the dynamic error of three transducers of the same model are evaluated using the Levitation Mass Method (LMM), in which the impact forces applied to the transducers are accurately determined as the inertial force of the moving part of the aerostatic linear bearing. The parameters for correcting the dynamic error are determined from the results of one set of impact measurements of one transducer. Then, the validity of the obtained parameters is evaluated using the results of the other sets of measurements of all the three transducers. The uncertainties in the uncorrected force and those in the corrected force are also estimated. If manufacturers determine the correction parameters for each model using the proposed method, and provide the software with the parameters corresponding to each model, then users can obtain the waveform corrected against dynamic error and its uncertainty. The present status and the future prospects of the developed software are discussed in this paper. PMID:25004158
Coulomb-corrected molecular orbital tomography of nitrogen
NASA Astrophysics Data System (ADS)
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-03-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation.
Coulomb-corrected molecular orbital tomography of nitrogen.
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-01-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation. PMID:27000666
Coulomb-corrected molecular orbital tomography of nitrogen
Zhai, Chunyang; He, Lixin; Lan, Pengfei; Zhu, Xiaosong; Li, Yang; Wang, Feng; Shi, Wenjing; Zhang, Qingbin; Lu, Peixiang
2016-01-01
High-order harmonic generation (HHG) from aligned molecules has provided a promising way to probe the molecular orbital with an Ångström resolution. This method, usually called molecular orbital tomography (MOT) replies on a simple assumption of the plane-wave approximation (PW), which has long been questioned due to that PW approximation is known to be valid in the keV energy region. However, the photon energy is usually no more than 100 eV in HHG. In this work, we experimentally reconstruct the highest occupied molecular orbital (HOMO) of nitrogen (N2) by using a Coulomb-corrected MOT (CCMOT) method. In our scheme, the molecular continuum states are described by a Coulomb wave function instead of the PW approximation. With CCMOT, the reconstructed orbital is demonstrated to agree well with the theoretical prediction and retain the main features of the HOMO of N2. Compared to the PW approximation method, the CCMOT shows a significant improvement in eliminating the artificial structures caused by PW approximation. PMID:27000666
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.
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.
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
Correcting for Purifying Selection: An Improved Human Mitochondrial Molecular Clock
Soares, Pedro; Ermini, Luca; Thomson, Noel; Mormina, Maru; Rito, Teresa; Röhl, Arne; Salas, Antonio; Oppenheimer, Stephen; Macaulay, Vincent; Richards, Martin B.
2009-01-01
There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at ∼15 kya that—unlike the uncorrected clock—matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55–70 kya, 5–20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses. PMID:19500773
Localised distributions and criteria for correctness in complex Langevin dynamics
Aarts, Gert; Giudice, Pietro; Seiler, Erhard
2013-10-15
Complex Langevin dynamics can solve the sign problem appearing in numerical simulations of theories with a complex action. In order to justify the procedure, it is important to understand the properties of the real and positive distribution, which is effectively sampled during the stochastic process. In the context of a simple model, we study this distribution by solving the Fokker–Planck equation as well as by brute force and relate the results to the recently derived criteria for correctness. We demonstrate analytically that it is possible that the distribution has support in a strip in the complexified configuration space only, in which case correct results are expected. -- Highlights: •Characterisation of the equilibrium distribution sampled in complex Langevin dynamics. •Connection between criteria for correctness and breakdown. •Solution of the Fokker–Planck equation in the case of real noise. •Analytical determination of support in complexified space.
Erba, Alessandro; Maul, Jefferson; Civalleri, Bartolomeo
2016-01-31
An ab initio quantum-mechanical theoretical framework is presented to compute the thermal properties of molecular crystals. The present strategy combines dispersion-corrected density-functional-theory (DFT-D), harmonic phonon dispersion, quasi-harmonic approximation to the lattice dynamics for thermal expansion and thermodynamic functions, and quasi-static approximation for anisotropic thermo-elasticity. The proposed scheme is shown to reliably describe thermal properties of the urea molecular crystal by a thorough comparison with experimental data. PMID:26670006
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.
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.
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.
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.
Sergiievskyi, Volodymyr P; Jeanmairet, Guillaume; Levesque, Maximilien; Borgis, Daniel
2014-06-01
Molecular density functional theory (MDFT) offers an efficient implicit-solvent method to estimate molecule solvation free-energies, whereas conserving a fully molecular representation of the solvent. Even within a second-order approximation for the free-energy functional, the so-called homogeneous reference fluid approximation, we show that the hydration free-energies computed for a data set of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by 2-3 orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canonical to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justification to empirical partial molar volume corrections that have been proposed recently. PMID:26273876
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.
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 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
Dynamic shear modulus of glycerol: Corrections due to instrument compliance
NASA Astrophysics Data System (ADS)
Schröter, K.; Hutcheson, S. A.; Shi, X.; Mandanici, A.; McKenna, G. B.
2006-12-01
A recent article by Shi et al. [J. Chem. Phys.123, 174507 (2005)] reports results from mechanical measurements on three simple inorganic glass formers: glycerol, m-toluidine, and sucrose benzoate. The experiments carried out were stress relaxation, aging, and dynamic (all in shear) using a torsional rheometer, an advanced rheometric expansion system (TA Instruments). The original force rebalance transducer (2KFRT) supplied with the system was replaced with a custom-made load cell (Sensotec) that had a capacity of 20000gcm in torque and 5000g in normal force. The replacement of the load cell was done due to the belief that the main source of compliance in this instrument was from the 2KFRT. With this assumption, the authors published their results for the three materials of interest and compared their results with the techniques of Schröter and Donth [J. Chem. Phys.113, 9101 (2000)] for the measurements on glycerol and reported important differences. These differences were disputed by one of the present authors (Schröter), and the present report shows that the results from Schröter and Donth are correct. We show that the reasons have to do with the instrument compliance being greater than originally thought by Shi et al. Here we examine the effects of platen diameter/geometry on the glycerol dynamic moduli, describe a means to correct dynamic data, present a revised comparison of the corrected data with that of Schröter and Donth, and provide a discussion of future work and conclusions.
Dynamically corrected gates for an exchange-only qubit
NASA Astrophysics Data System (ADS)
Hickman, G. T.; Wang, Xin; Kestner, J. P.; Das Sarma, S.
2013-10-01
We provide analytical composite pulse sequences that perform dynamical decoupling concurrently with arbitrary rotations for a qubit coded in the spin state of a triple quantum dot. The sequences are designed to respect realistic experimental constraints such as strictly nonnegative couplings. Logical errors and leakage errors are simultaneously corrected. A short pulse sequence is presented to compensate nuclear noise and a longer sequence is presented to simultaneously compensate nuclear and charge noise. The capability developed in this work provides a clear prescription for combatting the relevant sources of noise that currently hinder exchange-only qubit experiments.
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
Annihilation of craters: Molecular dynamic simulations on a silver surface
Henriksson, K. O. E.; Nordlund, K.; Keinonen, J.
2007-12-15
The ability of silver cluster ions containing 13 atoms to fill in a preexisting crater with a radius of about 28 A ring on a silver (001) target has been investigated using molecular dynamics simulations and the molecular-dynamics-Monte Carlo corrected effective medium potential. The largest lateral distance r between crater and ion was about three times the radius of the preexisting crater, namely, 75 A ring . The results reveal that when r<20 A ring and r>60 A ring the preexisting crater is partially filled in, and for other distances there is a net growth of the crater. The lattice damage created by the cluster ions, the total sputtering yield, the cluster sputtering yield, and simulated transmission electron microscopy images of the irradiated targets are also presented.
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
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.
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.
Leading concentration correction to polymer dynamic self-structure factor
NASA Astrophysics Data System (ADS)
Perico, Angelo; La Ferla, Roberto; Freed, Karl F.
1989-10-01
The discrete chain representation of multiple scattering theory of the concentration dependence of the hydrodynamics of polymer solutions is applied to the calculation of the leading concentration correction to the dynamic structure factor S(k,t) and its first and second cumulants of individual labeled Gaussian chains in theta solutions at nonzero concentrations. Contributions are separated into those from overall translational and internal chain motions as well as couplings between different internal modes and between translation and internal modes, coupling that are introduced by interchain hydrodynamic interactions. The separate contributions are analyzed as a function of k and of t in order to isolate regions where certain contributions are dominant. As expected, short times and larger k tend to favor contributions from internal chain dynamics, while longer times and smaller k make concentration dependent translational effects predominate. Computations for shorter chains are extrapolated to provide the asymptotic long chain behavior.
Susceptibility artefact correction using dynamic graph cuts: application to neurosurgery.
Daga, Pankaj; Pendse, Tejas; Modat, Marc; White, Mark; Mancini, Laura; Winston, Gavin P; McEvoy, Andrew W; Thornton, John; Yousry, Tarek; Drobnjak, Ivana; Duncan, John S; Ourselin, Sebastien
2014-10-01
Echo Planar Imaging (EPI) is routinely used in diffusion and functional MR imaging due to its rapid acquisition time. However, the long readout period makes it prone to susceptibility artefacts which results in geometric and intensity distortions of the acquired image. The use of these distorted images for neuronavigation hampers the effectiveness of image-guided surgery systems as critical white matter tracts and functionally eloquent brain areas cannot be accurately localised. In this paper, we present a novel method for correction of distortions arising from susceptibility artefacts in EPI images. The proposed method combines fieldmap and image registration based correction techniques in a unified framework. A phase unwrapping algorithm is presented that can efficiently compute the B0 magnetic field inhomogeneity map as well as the uncertainty associated with the estimated solution through the use of dynamic graph cuts. This information is fed to a subsequent image registration step to further refine the results in areas with high uncertainty. This work has been integrated into the surgical workflow at the National Hospital for Neurology and Neurosurgery and its effectiveness in correcting for geometric distortions due to susceptibility artefacts is demonstrated on EPI images acquired with an interventional MRI scanner during neurosurgery. PMID:25047865
Dynamic Distortion Correction for Endoscopy Systems with Exchangeable Optics
NASA Astrophysics Data System (ADS)
Stehle, Thomas; Hennes, Michael; Gross, Sebastian; Behrens, Alexander; Wulff, Jonas; Aach, Til
Endoscopic images are strongly affected by lens distortion caused by the use of wide angle lenses. In case of endoscopy systems with exchangeable optics, e.g. in bladder endoscopy or sinus endoscopy, the camera sensor and the optics do not form a rigid system but they can be shifted and rotated with respect to each other during an examination. This flexibility has a major impact on the location of the distortion centre as it is moved along with the optics. In this paper, we describe an algorithm for the dynamic correction of lens distortion in cystoscopy which is based on a one time calibration. For the compensation, we combine a conventional static method for distortion correction with an algorithm to detect the position and the orientation of the elliptic field of view. This enables us to estimate the position of the distortion centre according to the relative movement of camera and optics. Therewith, a distortion correction for arbitrary rotation angles and shifts becomes possible without performing static calibrations for every possible combination of shifts and angles beforehand.
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
Parallel molecular dynamics on a multi signalprocessor system
NASA Astrophysics Data System (ADS)
Scott, W.; Gunzinger, A.; Bäumle, B.; Kohler, P.; Müller, U. A.; Mühll, H.-R. Vonder; Eichenberger, A.; Guggenbühl, W.; Ironmonger, N.; Müller-Plathe, F.; van Gunsteren, W. F.
1993-04-01
This paper gives an overview of a parallel computer architecture called MUSIC (Multi Signalprocessor System with Intelligent Communication), which has been developed at the Swiss Federal Institute of Technology. The current version achieves a peak performance of 3.8 GFlops. We discuss the system software and tools used to program the system and then present our implementation of a molecular dynamics simulation program which uses the architecture of MUSIC in an efficient way. We demonstrate the correctness of our implementation and give measurements of the performance of the system. To the best of our knowledge, MUSIC outperforms the most powerful present-day vector supercomputers.
Danel, J.-F.; Blottiau, P.; Kazandjian, L.; Piron, R.; Torrent, M.
2014-10-15
The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on a semiclassical approximation and possibly on a gradient correction, is a simulation method available at high temperature. For a high-Z element such as lutetium, we examine how orbital-free molecular dynamics applied to the equation of state of a dense plasma can be regarded as the limit of quantum molecular dynamics at high temperature. For the normal mass density and twice the normal mass density, we show that the pressures calculated with the quantum approach converge monotonically towards those calculated with the orbital-free approach; we observe a faster convergence when the orbital-free approach includes the gradient correction. We propose a method to obtain an equation of state reproducing quantum molecular dynamics results up to high temperatures where this approach cannot be directly implemented. With the results already obtained for low-Z plasmas, the present study opens the way for reproducing the quantum molecular dynamics pressure for all elements up to high temperatures.
Shadow correction in high dynamic range images for generating orthophotos
NASA Astrophysics Data System (ADS)
Suzuki, Hideo; Chikatsu, Hirofumi
2011-07-01
High dynamic range imagery is widely used in remote sensing. With the widespread use of aerial digital cameras such as the DMC, ADS40, RMK-D, and UltraCamD, high dynamic range imaging is generally expected for generating minuteness orthophotos in digital aerial photogrammetry. However, high dynamic range images (12-bit, 4,096 gray levels) are generally compressed into an 8-bit depth digital image (256 gray levels) owing to huge amount of data and interface with peripherals such as monitors and printers. This means that a great deal of image data is eliminated from the original image, and this introduces a new shadow problem. In particular, the influence of shadows in urban areas causes serious problems when generating minuteness orthophotos and performing house detection. Therefore, shadow problems can be solved by addressing the image compression problems. There is a large body of literature on image compression techniques such as logarithmic compression and tone mapping algorithms. However, logarithmic compression tends to cause loss of details in dark and/or light areas. Furthermore, the logarithmic method intends to operate on the full scene. This means that high-resolution luminance information can not be obtained. Even though tone mapping algorithms have the ability to operate over both full scene and local scene, background knowledge is required. To resolve the shadow problem in digital aerial photogrammetry, shadow areas should be recognized and corrected automatically without the loss of luminance information. To this end, a practical shadow correction method using 12-bit real data acquired by DMC is investigated in this paper.
Telescopic systems with dynamic nonlinear optical correction for distortions
Vasil'ev, Michail V; Venediktov, Vladimir Yu; Leshchev, Alexey A
2001-01-31
The review of basic achievements in the field of non-linear adaptive optics is presented. In particular, schematics and properties of adaptive optical telescopes considered in which the image distortions introduced by defects of the primary mirror and other optical elements are compensated by nonlinear optical methods. The conventional methods of laser optics, such as phase conjugation and dynamic holography, make it possible both to solve the problems of classical (imaging) optics related to the building of telescopes for imaging remote objects with high resolution, which are based on large, light-weight or sectional mirrors, and create the systems that produce laser beams with the high-quality wave front. The basic designs of such telescopes are considered and the possibilities of corrections for distortions in them are analysed and confirmed by experiments. (review)
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
Frequency-domain correction of sensor dynamic error for step response.
Yang, Shuang-Long; Xu, Ke-Jun
2012-11-01
To obtain accurate results in dynamic measurements it is required that the sensors should have good dynamic performance. In practice, sensors have non-ideal dynamic characteristics due to their small damp ratios and low natural frequencies. In this case some dynamic error correction methods can be adopted for dealing with the sensor responses to eliminate the effect of their dynamic characteristics. The frequency-domain correction of sensor dynamic error is a common method. Using the existing calculation method, however, the correct frequency-domain correction function (FCF) cannot be obtained according to the step response calibration experimental data. This is because of the leakage error and invalid FCF value caused by the cycle extension of the finite length step input-output intercepting data. In order to solve these problems the data splicing preprocessing and FCF interpolation are put forward, and the FCF calculation steps as well as sensor dynamic error correction procedure by the calculated FCF are presented in this paper. The proposed solution is applied to the dynamic error correction of the bar-shaped wind tunnel strain gauge balance so as to verify its effectiveness. The dynamic error correction results show that the adjust time of the balance step response is shortened to 10 ms (shorter than 1/30 before correction) after frequency-domain correction, and the overshoot is fallen within 5% (less than 1/10 before correction) as well. The dynamic measurement accuracy of the balance is improved significantly. PMID:23206091
Frequency-domain correction of sensor dynamic error for step response
NASA Astrophysics Data System (ADS)
Yang, Shuang-Long; Xu, Ke-Jun
2012-11-01
To obtain accurate results in dynamic measurements it is required that the sensors should have good dynamic performance. In practice, sensors have non-ideal dynamic characteristics due to their small damp ratios and low natural frequencies. In this case some dynamic error correction methods can be adopted for dealing with the sensor responses to eliminate the effect of their dynamic characteristics. The frequency-domain correction of sensor dynamic error is a common method. Using the existing calculation method, however, the correct frequency-domain correction function (FCF) cannot be obtained according to the step response calibration experimental data. This is because of the leakage error and invalid FCF value caused by the cycle extension of the finite length step input-output intercepting data. In order to solve these problems the data splicing preprocessing and FCF interpolation are put forward, and the FCF calculation steps as well as sensor dynamic error correction procedure by the calculated FCF are presented in this paper. The proposed solution is applied to the dynamic error correction of the bar-shaped wind tunnel strain gauge balance so as to verify its effectiveness. The dynamic error correction results show that the adjust time of the balance step response is shortened to 10 ms (shorter than 1/30 before correction) after frequency-domain correction, and the overshoot is fallen within 5% (less than 1/10 before correction) as well. The dynamic measurement accuracy of the balance is improved significantly.
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.
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 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
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
Motion correction via nonrigid coregistration of dynamic MR mammography series
NASA Astrophysics Data System (ADS)
Krol, Andrzej; Magri, Alphonso; Unlu, Mehmet; Feiglin, David; Lipson, Edward; Mandel, James; Tillapaugh-Fay, Gwen; Lee, Wei; Coman, Ioana; Szeverenyi, Nikolaus M.
2006-03-01
The objectives of this investigation are to improve quality of subtraction MR breast images and improve accuracy of time-signal intensity curves (TSIC) related to local contrast-agent concentration in dynamic MR mammography. The patients, with up to nine fiducial skin markers (FSMs) taped to each breast, were prone with both breasts suspended into a single well that housed the receiver coil. After a preliminary scan, paramagnetic contrast agent gadopentate digmeglumine (Gd) was delivered intravenously, followed by physiological saline. The field of view was centered over the breasts. We used a gradient recalled echo (GRE) technique for pre-Gd baseline, and five more measurements at 60s intervals. Centroids were determined for corresponding FSMs visible on pre-Gd and any post-Gd images. This was followed by segmentation of breast surfaces in all dynamic-series images, and meshing of all post-Gd breast images. Tetrahedral volume and triangular surface elements were used to construct a finite element method (FEM) model. We used ANSYS TM software and an analogy between orthogonal components of the displacement field and the temperature differences in steady-state heat transfer (SSHT) in solids. The floating images were warped to a fixed image using an appropriate shape function for interpolation from mesh nodes to voxels. To reduce any residual misregistration, we performed surface matching between the previously warped floating image and the target image. Our method of motion correction via nonrigid coregistration yielded excellent differential-image series that clearly revealed lesions not visible in unregistered differential-image series. Further, it produced clinically useful maximum intensity projection (MIP) 3D images.
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.
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.
2015-11-01
In the article by Heuslein et al, which published online ahead of print on September 3, 2015 (DOI: 10.1161/ATVBAHA.115.305775), a correction was needed. Brett R. Blackman was added as the penultimate author of the article. The article has been corrected for publication in the November 2015 issue. PMID:26490278
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 simulation of liquid water: Hybrid density functionals
Todorova, T; Seitsonen, A; Hutter, J; Kuo, W; Mundy, C
2005-09-12
The structure, dynamical and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta functional, four gradient corrected functionals, the local density and Hartree-Fock approximation. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and under-structured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller averaged numbers of hydrogen bonds and similar hydrogen bond populations as pure density functionals. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than from the corresponding pure density functionals.
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.
Validation of Bubble Dynamics Equation for a Nano-scale Bubble via Molecular Dynamics Simulation
NASA Astrophysics Data System (ADS)
Tsuda, S.; Hyodo, H.; Watanabe, S.
2015-12-01
For a validation of the application of conventional bubble dynamics to a nano-scale bubble behaviour, we simulated a nano-scale bubble collapsing or vibration by Molecular Dynamics (MD) method and compared the result with the solution of Rayleigh-Plesset (RP) equation and that of Confined RP (CRP) equation, whose boundary condition was corrected to be consistent with that of MD simulation. As a result, a good coincidence was obtained between MD, RP, and CRP in the case of one-component fluid. In addition, also a good correspondence was obtained particularly in the comparison between MD and CRP in the case of two-component fluid containing non-condensable gas. The present results indicate that conventional bubble dynamics equation can be applied even to a nano-scale tiny bubble.
Accelerated molecular dynamics methods: introduction and recent developments
Uberuaga, Blas Pedro; Voter, Arthur F; Perez, Danny; Shim, Y; Amar, J G
2009-01-01
A long-standing limitation in the use of molecular dynamics (MD) simulation is that it can only be applied directly to processes that take place on very short timescales: nanoseconds if empirical potentials are employed, or picoseconds if we rely on electronic structure methods. Many processes of interest in chemistry, biochemistry, and materials science require study over microseconds and beyond, due either to the natural timescale for the evolution or to the duration of the experiment of interest. Ignoring the case of liquids xxx, the dynamics on these time scales is typically characterized by infrequent-event transitions, from state to state, usually involving an energy barrier. There is a long and venerable tradition in chemistry of using transition state theory (TST) [10, 19, 23] to directly compute rate constants for these kinds of activated processes. If needed dynamical corrections to the TST rate, and even quantum corrections, can be computed to achieve an accuracy suitable for the problem at hand. These rate constants then allow them to understand the system behavior on longer time scales than we can directly reach with MD. For complex systems with many reaction paths, the TST rates can be fed into a stochastic simulation procedure such as kinetic Monte Carlo xxx, and a direct simulation of the advance of the system through its possible states can be obtained in a probabilistically exact way. A problem that has become more evident in recent years, however, is that for many systems of interest there is a complexity that makes it difficult, if not impossible, to determine all the relevant reaction paths to which TST should be applied. This is a serious issue, as omitted transition pathways can have uncontrollable consequences on the simulated long-time kinetics. Over the last decade or so, we have been developing a new class of methods for treating the long-time dynamics in these complex, infrequent-event systems. Rather than trying to guess in advance what
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.
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
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
2015-12-01
In the article by Narayan et al (Narayan O, Davies JE, Hughes AD, Dart AM, Parker KH, Reid C, Cameron JD. Central aortic reservoir-wave analysis improves prediction of cardiovascular events in elderly hypertensives. Hypertension. 2015;65:629–635. doi: 10.1161/HYPERTENSIONAHA.114.04824), which published online ahead of print December 22, 2014, and appeared in the March 2015 issue of the journal, some corrections were needed.On page 632, Figure, panel A, the label PRI has been corrected to read RPI. In panel B, the text by the upward arrow, "10% increase in kd,” has been corrected to read, "10% decrease in kd." The corrected figure is shown below.The authors apologize for these errors. PMID:26558821
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.
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
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-28
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics
NASA Astrophysics Data System (ADS)
Petsev, Nikolai D.; Leal, L. Gary; Shell, M. Scott
2015-01-01
We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.
Counterpoise-corrected interaction energy analysis based on the fragment molecular orbital scheme
NASA Astrophysics Data System (ADS)
Okiyama, Yoshio; Fukuzawa, Kaori; Yamada, Haruka; Mochizuki, Yuji; Nakano, Tatsuya; Tanaka, Shigenori
2011-06-01
Basis set superposition error (BSSE) correction with counterpoise (CP) procedure under the environmental electrostatic potential is newly introduced to interfragment interaction energy (IFIE), which is important for interaction analysis in the fragment molecular orbital method. The CP correction for IFIE is applied to a stacked dimer of base pair and a protein-ligand complex of estrogen receptor and 17β-estradiol with scaled third-order Møller-Plesset perturbation theory. The BSSEs amount to about quarter of IFIE for hydrogen-bonding and electrostatic interactions and half or even more for dispersion interactions. Estimation of IFIE with the CP correction is therefore preferred for the quantitative discussion.
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.
Zhang, Lin; Tang, Ronghong; Bai, Shu; Connors, Natalie K; Lua, Linda H L; Chuan, Yap P; Middelberg, Anton P J; Sun, Yan
2013-05-01
Virus-like particles (VLPs) are highly organized nanoparticles that have great potential in vaccinology, gene therapy, drug delivery, and materials science. However, the application of VLPs is hindered by obstacles in their design and production due to low efficiency of self-assembly. In the present study, all-atom (AA) molecular dynamics (MD) simulations coupled with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method are utilized to examine the molecular interactions in the capsomere of a murine polyomavirus (MPV) VLP. It is found that both low ionic strength and the intracapsomere disulfide bonds are favorable for maintaining a stable capsomere. Simulation results examining the effects of solution conditions on the stabilization of a capsomere were verified by calorimetry experiments. Simulation results of free energy decomposition indicate that hydrophobic interaction is favorable for the formation of a capsomere, whereas electrostatic interaction is unfavorable. With increasing ionic strength, the dominant interaction for the stabilization of a capsomere changes from hydrophobic to electrostatic. By comprehensive analyses, the key amino acid residues (hot spots) in VP1 protein aiding formation of a capsomere in different solution conditions have been identified. These results provide molecular insights into the stabilization of building blocks for VLP and are expected to have implications in their partitioning between the correct and off-pathway reactions in VLP assembly. PMID:23586433
NASA Astrophysics Data System (ADS)
1995-04-01
Seismic images of the Brooks Range, Arctic Alaska, reveal crustal-scale duplexing: Correction Geology, v. 23, p. 65 68 (January 1995) The correct Figure 4A, for the loose insert, is given here. See Figure 4A below. Corrected inserts will be available to those requesting copies of the article from the senior author, Gary S. Fuis, U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025. Figure 4A. P-wave velocity model of Brooks Range region (thin gray contours) with migrated wide-angle reflections (heavy red lines) and migreated vertical-incidence reflections (short black lines) superimposed. Velocity contour interval is 0.25 km/s; 4,5, and 6 km/s contours are labeled. Estimated error in velocities is one contour interval. Symbols on faults shown at top are as in Figure 2 caption.
Molecular Dynamics Simulations of Temperature Equilibration in Dense Hydrogen
Glosli, J; Graziani, F; More, R; Murillo, M; Streitz, F; Surh, M; Benedict, L; Hau-Riege, S; Langdon, A; London, R
2008-02-14
The temperature equilibration rate in dense hydrogen (for both T{sub i} > T{sub e} and T{sub i} < T{sub e}) has been calculated with large-scale molecular dynamics simulations for temperatures between 10 and 300 eV and densities between 10{sup 20}/cc to 10{sup 24}/cc. Careful attention has been devoted to convergence of the simulations, including the role of semiclassical potentials. We find that for Coulomb logarithms L {approx}> 1, Brown-Preston-Singleton [Brown et al., Phys. Rep. 410, 237 (2005)] with the sub-leading corrections and the fit of Gericke-Murillo-Schlanges [Gericke et al., PRE 65, 036418 (2003)] to the T-matrix evaluation of the collision operator, agrees with the MD data to within the error bars of the simulation. For more strongly-coupled plasmas where L {approx}< 1, our numerical results are consistent with the fit of Gericke-Murillo-Schlanges.
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.
Quantum Thermal Bath for Path Integral Molecular Dynamics Simulation.
Brieuc, Fabien; Dammak, Hichem; Hayoun, Marc
2016-03-01
The quantum thermal bath (QTB) method has been recently developed to account for the quantum nature of the nuclei by using standard molecular dynamics (MD) simulation. QTB-MD is an efficient but approximate method when dealing with strongly anharmonic systems, while path integral molecular dynamics (PIMD) gives exact results but in a huge amount of computation time. The QTB and PIMD methods have been combined in order to improve the PIMD convergence or correct the failures of the QTB-MD technique. Therefore, a new power spectral density of the random force within the QTB has been developed. A modified centroid-virial estimator of the kinetic energy, especially adapted to QTB-PIMD, has also been proposed. The method is applied to selected systems: a one-dimensional double-well system, a ferroelectric phase transition, and the position distribution of an hydrogen atom in a fuel cell material. The advantage of the QTB-PIMD method is its ability to give exact results with a more reasonable computation time for strongly anharmonic systems. PMID:26799437
Accelerated Molecular Dynamics studies of He Bubble Growth in Tungsten
NASA Astrophysics Data System (ADS)
Uberuaga, Blas; Sandoval, Luis; Perez, Danny; Voter, Arthur
2015-11-01
Understanding how materials respond to extreme environments is critical for predicting and improving performance. In materials such as tungsten exposed to plasmas for nuclear fusion applications, novel nanoscale fuzzes, comprised of tendrils of tungsten, form as a consequence of the implantation of He into the near surface. However, the detailed mechanisms that link He bubble formation to the ultimate development of fuzz are unclear. Molecular dynamics simulations provide insight into the He implantation process, but are necessarily performed at implantation rates that are orders of magnitudes faster than experiment. Here, using accelerated molecular dynamics methods, we examine the role of He implantation rates on the physical evolution of He bubbles in tungsten. We find that, as the He rate is reduced, new types of events involving the response of the tungsten matrix to the pressure in the bubble become competitive and change the overall evolution of the bubble as well as the subsequent morphology of the tungsten surface. We have also examined how bubble growth differs at various microstructural features. These results highlight the importance of performing simulations at experimentally relevant conditions in order to correctly capture the contributions of the various significant kinetic processes and predict the overall response of the material.
Anharmonic infrared and Raman spectra in Car-Parrinello molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Pagliai, Marco; Cavazzoni, Carlo; Cardini, Gianni; Erbacci, Giovanni; Parrinello, Michele; Schettino, Vincenzo
2008-06-01
The infrared and Raman spectra of naphthalene crystal with inclusion of anharmonic effects have been calculated by adopting the generalized variational density functional perturbation theory in the framework of Car-Parrinello molecular dynamics simulations. The computational approach has been generalized for cells of arbitrary shape. The intermolecular interactions have been analyzed with and without the van der Waals corrections, showing the importance of such interactions in the naphthalene crystal to reproduce the structural, dynamical, and spectroscopic properties.
2016-02-01
Neogi T, Jansen TLTA, Dalbeth N, et al. 2015 Gout classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Ann Rheum Dis 2015;74:1789–98. The name of the 20th author was misspelled. The correct spelling is Janitzia Vazquez-Mellado. We regret the error. PMID:26881284
Local Dynamic Reactive Power for Correction of System Voltage Problems
Kueck, John D; Rizy, D Tom; Li, Fangxing; Xu, Yan; Li, Huijuan; Adhikari, Sarina; Irminger, Philip
2008-12-01
Distribution systems are experiencing outages due to a phenomenon known as local voltage collapse. Local voltage collapse is occurring in part because modern air conditioner compressor motors are much more susceptible to stalling during a voltage dip than older motors. These motors can stall in less than 3 cycles (.05s) when a fault, such as on the sub-transmission system, causes voltage to sag to 70 to 60%. The reasons for this susceptibility are discussed in the report. During the local voltage collapse, voltages are depressed for a period of perhaps one or two minutes. There is a concern that these local events are interacting together over larger areas and may present a challenge to system reliability. An effective method of preventing local voltage collapse is the use of voltage regulation from Distributed Energy Resources (DER) that can supply or absorb reactive power. DER, when properly controlled, can provide a rapid correction to voltage dips and prevent motor stall. This report discusses the phenomenon and causes of local voltage collapse as well as the control methodology we have developed to counter voltage sag. The problem is growing because of the use of low inertia, high efficiency air conditioner (A/C) compressor motors and because the use of electric A/C is growing in use and becoming a larger percentage of system load. A method for local dynamic voltage regulation is discussed which uses reactive power injection or absorption from local DER. This method is independent, rapid, and will not interfere with conventional utility system voltage control. The results of simulations of this method are provided. The method has also been tested at the ORNL s Distributed Energy Communications and Control (DECC) Laboratory using our research inverter and synchronous condenser. These systems at the DECC Lab are interconnected to an actual distribution system, the ORNL distribution system, which is fed from TVA s 161kV sub-transmission backbone. The test results
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..
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.
2016-02-01
In the article by Guessous et al (Guessous I, Pruijm M, Ponte B, Ackermann D, Ehret G, Ansermot N, Vuistiner P, Staessen J, Gu Y, Paccaud F, Mohaupt M, Vogt B, Pechère-Bertschi A, Martin PY, Burnier M, Eap CB, Bochud M. Associations of ambulatory blood pressure with urinary caffeine and caffeine metabolite excretions. Hypertension. 2015;65:691–696. doi: 10.1161/HYPERTENSIONAHA.114.04512), which published online ahead of print December 8, 2014, and appeared in the March 2015 issue of the journal, a correction was needed.One of the author surnames was misspelled. Antoinette Pechère-Berstchi has been corrected to read Antoinette Pechère-Bertschi.The authors apologize for this error. PMID:26763012
Walters, Chris; Errington, Neil; Rowe, Arther J; Harding, Stephen E
2002-01-01
Over recent years the binding ability of the molecular chaperone cpn60 (GroEL14) and its co-chaperone cpn10 (GroES7) has been reported to occur under an assortment of specific conditions from the use of non-hydrolysable ATP analogues (namely adenosine 5'-[gamma-thio]triphosphate) to requiring hydrolysable ATP for any interaction to occur. We have investigated this further using the molecular hydrodynamic methods (hydrodynamic bead modelling, sedimentation-velocity analytical ultracentrifugation and dynamic light-scattering), allowing the process to be followed under physiologically relevant dilute solution conditions, combined with absorption spectrophotometry to determine GroES7-GroEL14 interaction through the rate inhibition of the cpn60's ATPase activity by GroES7. The results found here indicate that the presence of hydrolysable ATP is required to facilitate correct GroES7 interaction with GroEL14 in solution. PMID:12049650
Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides
NASA Technical Reports Server (NTRS)
Sokalski, W. A.; Keller, D. A.; Ornstein, R. L.; Rein, R.
1993-01-01
The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.
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.
Dynamic pressure corrections in a clearance-sealed piston prover for gas flow measurements
NASA Astrophysics Data System (ADS)
Kutin, Jože; Bobovnik, Gregor; Bajsić, Ivan
2013-02-01
The dynamic pressure effects and their corrections in a high-speed, clearance-sealed realization of a piston prover for gas flow measurements are discussed. The experimental results show the deterministic, rather than stochastic, nature of the dynamic pressure conditions and, consequently, the repeatable nature of their influence on the flow measurements. The experimental validation proves the advantage of the polytropic/adiabatic pressure correction model, which was proposed by the authors, as compared with the isothermal pressure correction model. The paper ends with an estimation of the measurement uncertainty related to the pressure corrections using either the adiabatic or isothermal model.
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
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
Molecular-dynamics simulation of hydrogen diffusion in palladium
NASA Astrophysics Data System (ADS)
Li, Yinggang; Wahnström, Göran
1992-12-01
Molecular-dynamics simulations for hydrogen diffusion in Pd are performed for a system consisting of 256 Pd atoms and 8 H atoms at the temperature T=623 K. Under these conditions detailed quasielastic-neutron-scattering (QNS) data are available. For the interatomic interactions we use the embedded-atom method (EAM), which incorporates some essential many-body effects in metals. Based on the EAM approach, the wave-vector dependence of the width of the QNS peak is investigated in detail. It is found that a single electronically adiabatic potential-energy surface cannot reproduce the observed wave-vector dependence. After incorporating the coupling of hydrogen atoms to the low-lying electron-hole pair excitations among the conduction electrons, close agreement with the experimental data is obtained. This is a strong indication that one has to go beyond the Born-Oppenheimer approximation in order to characterize correctly the diffusive motion of hydrogen in metals. To reveal the diffusive behavior in more detail, the residence time distribution and the correlation character in diffusion direction are investigated. We found that including the nonadiabatic corrections reduces the probability for the H atoms to move over several lattice sites without getting trapped in between. As a result, the motion of the H atoms becomes more similar to that assumed in the Chudley-Elliott model, which describes well the QNS data for the wave-vector dependence of the width.
Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations.
Soares, Rosemberg O; Torres, Pedro H M; da Silva, Manuela L; Pascutti, Pedro G
2016-08-01
The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding. PMID:27291071
Automated motion correction based on target tracking for dynamic nuclear medicine studies
NASA Astrophysics Data System (ADS)
Cao, Xinhua; Tetrault, Tracy; Fahey, Fred; Treves, Ted
2008-03-01
Nuclear medicine dynamic studies of kidneys, bladder and stomach are important diagnostic tools. Accurate generation of time-activity curves from regions of interest (ROIs) requires that the patient remains motionless for the duration of the study. This is not always possible since some dynamic studies may last from several minutes to one hour. Several motion correction solutions have been explored. Motion correction using external point sources is inconvenient and not accurate especially when motion results from breathing, organ motion or feeding rather than from body motion alone. Centroid-based motion correction assumes that activity distribution is only inside the single organ (without background) and uniform, but this approach is impractical in most clinical studies. In this paper, we present a novel technique of motion correction that first tracks the organ of interest in a dynamic series then aligns the organ. The implementation algorithm for target tracking-based motion correction consists of image preprocessing, target detection, target positioning, motion estimation and prediction, tracking (new search region generation) and target alignment. The targeted organ is tracked from the first frame to the last one in the dynamic series to generate a moving trajectory of the organ. Motion correction is implemented by aligning the organ ROIs in the image series to the location of the organ in the first image. The proposed method of motion correction has been applied to several dynamic nuclear medicine studies including radionuclide cystography, dynamic renal scintigraphy, diuretic renography and gastric emptying scintigraphy.
2015-05-22
The Circulation Research article by Keith and Bolli (“String Theory” of c-kitpos Cardiac Cells: A New Paradigm Regarding the Nature of These Cells That May Reconcile Apparently Discrepant Results. Circ Res. 2015:116:1216-1230. doi: 10.1161/CIRCRESAHA.116.305557) states that van Berlo et al (2014) observed that large numbers of fibroblasts and adventitial cells, some smooth muscle and endothelial cells, and rare cardiomyocytes originated from c-kit positive progenitors. However, van Berlo et al reported that only occasional fibroblasts and adventitial cells derived from c-kit positive progenitors in their studies. Accordingly, the review has been corrected to indicate that van Berlo et al (2014) observed that large numbers of endothelial cells, with some smooth muscle cells and fibroblasts, and more rarely cardiomyocytes, originated from c-kit positive progenitors in their murine model. The authors apologize for this error, and the error has been noted and corrected in the online version of the article, which is available at http://circres.ahajournals.org/content/116/7/1216.full ( PMID:25999426
NASA Astrophysics Data System (ADS)
1998-12-01
Alleged mosasaur bite marks on Late Cretaceous ammonites are limpet (patellogastropod) home scars Geology, v. 26, p. 947 950 (October 1998) This article had the following printing errors: p. 947, Abstract, line 11, “sepia” should be “septa” p. 947, 1st paragraph under Introduction, line 2, “creep” should be “deep” p. 948, column 1, 2nd paragraph, line 7, “creep” should be “deep” p. 949, column 1, 1st paragraph, line 1, “creep” should be “deep” p. 949, column 1, 1st paragraph, line 5, “19774” should be “1977)” p. 949, column 1, 4th paragraph, line 7, “in particular” should be “In particular” CORRECTION Mammalian community response to the latest Paleocene thermal maximum: An isotaphonomic study in the northern Bighorn Basin, Wyoming Geology, v. 26, p. 1011 1014 (November 1998) An error appeared in the References Cited. The correct reference appears below: Fricke, H. C., Clyde, W. C., O'Neil, J. R., and Gingerich, P. D., 1998, Evidence for rapid climate change in North America during the latest Paleocene thermal maximum: Oxygen isotope compositions of biogenic phosphate from the Bighorn Basin (Wyoming): Earth and Planetary Science Letters, v. 160, p. 193 208.
[Performance analysis and radiometric correction of novel molecular hyperspectral imaging system].
Liu, Hong-Ying; Li, Qing-Li; Gu, Bin; Wang, Yi-Ting; Xue, Yong-Qi
2012-11-01
Integrating molecular imaging technology and hyperspectral technology, a novel molecular hyperspectral imaging (MHSI) system based on AOTF was presented. The system consists of microscope, spectrometer, matrix CCD, image collection card and computer. The system's performance was synthetically evaluated referring every part's performance. The spectral range of the MHSI system is from 550 to 1 000 nm. Two hundred twenty five bands can be continuously captured at a time. The spectral resolution is less than 2 nm. The spatial resolution is about 0.061 5 microm. CCD acquisition speed achieved 2.612 5 s x B(-1) in the integration mode and about 0.11 micros x B(-1) in the non-integration mode. Due to the infection of lamp, a spectral curve extracted directly from the original hyperspectral data can not truly present biochemical character and needs to be corrected. The paper proposes the gray correction coefficient algorithm with spatial dimension and spectral dimension, and gives concrete realization of the algorithm. Taking the sample of leukemia blood, by comparing the single-band images, pseudo-color images and spectra before and after correction, the results indicate the effectiveness of correction algorithm. The corrected data is effective for subsequent analysis. PMID:23387200
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.
Modeling and Bio molecular Self-assembly via Molecular Dynamics and Dissipative Particle Dynamics
NASA Astrophysics Data System (ADS)
Rakesh, L.
2009-09-01
Surfactants like materials can be used to increase the solubility of poorly soluble drugs in water and to increase drug bioavailability. A typical case study will be demonstrated using DPD simulation to model the distribution of anti-inflammatory drug molecules. Computer simulation is a convenient approach to understand drug distribution and solubility concepts without much wastage and costly experiments in the laboratory. Often in molecular dynamics (MD) the atoms are represented explicitly and the equation of motion as described by Newtonian dynamics is integrated explicitly. MD has been used to study spontaneous formation of micelles by hydrophobic molecules with amphiphilic head groups in bulk water, as well as stability of pre-configured micelles and membranes. DPD is a state-of the- art mesoscale simulation, it is a more recent molecular dynamics technique, originally developed for simulating complex fluids but lately also applied to membrane dynamics, hemodynamic in biomedical applications. Such fluids pervade industrial research from paints to pharmaceuticals and from cosmetics to the controlled release of drugs. Dissipative particle dynamics (DPD) can provide structural and dynamic properties of fluids in equilibrium, under shear or confined to narrow cavities, at length- and time-scales beyond the scope of traditional atomistic molecular dynamics simulation methods. Mesoscopic particles are used to represent clusters of molecules. The interaction conserves mass and momentum and as a consequence the dynamics is consistent with Navier-Stokes equations. In addition to the conservative forces, stochastic drive and dissipation is introduced to represent internal degrees of freedom in the mesoscopic particles. In this research, an initial study is being conducted using the aqueous solubilization of the nonsteroidal, anti-inflammatory drug is studied theoretically in micellar solution of nonionic (dodecyl hexa(ethylene oxide), C12E6) surfactants possessing the
Thermal transpiration: A molecular dynamics study
NASA Astrophysics Data System (ADS)
T, Joe Francis; Sathian, Sarith P.
2014-12-01
Thermal transpiration is a phenomenon where fluid molecules move from the cold end towards the hot end of a channel under the influence of longitudinal temperature gradient alone. Although the phenomenon of thermal transpiration is observed at rarefied gas conditions in macro systems, the phenomenon can occur at atmospheric pressure if the characteristic dimensions of the channel is less than 100 nm. The flow through these nanosized channels is characterized by the free molecular flow regimes and continuum theory is inadequate to describe the flow. Thus a non-continuum method like molecular dynamics (MD) is necessary to study such phenomenon. In the present work, MD simulations were carried out to investigate the occurance of thermal transpiration in copper and platinum nanochannels at atmospheric pressure conditions. The mean pressure of argon gas confined inside the nano channels was maintained around 1 bar. The channel height is maintained at 2nm. The argon atoms interact with each other and with the wall atoms through the Lennard-Jones potential. The wall atoms are modelled using an EAM potential. Further, separate simulations were carried out where a Harmonic potential is used for the atom-atom interaction in the platinum channel. A thermally insulating wall was introduced between the low and high temperature regions and those wall atoms interact with fluid atoms through a repulsive potential. A reduced cut off radius were used to achieve this. Thermal creep is induced by applying a temperature gradient along the channel wall. It was found that flow developed in the direction of the increasing temperature gradient of the wall. An increase in the volumetric flux was observed as the length of the cold and the hot regions of the wall were increased. The effect of temperature gradient and the wall-fluid interaction strength on the flow parameters have been studied to understand the phenomenon better.
Thermal transpiration: A molecular dynamics study
T, Joe Francis; Sathian, Sarith P.
2014-12-09
Thermal transpiration is a phenomenon where fluid molecules move from the cold end towards the hot end of a channel under the influence of longitudinal temperature gradient alone. Although the phenomenon of thermal transpiration is observed at rarefied gas conditions in macro systems, the phenomenon can occur at atmospheric pressure if the characteristic dimensions of the channel is less than 100 nm. The flow through these nanosized channels is characterized by the free molecular flow regimes and continuum theory is inadequate to describe the flow. Thus a non-continuum method like molecular dynamics (MD) is necessary to study such phenomenon. In the present work, MD simulations were carried out to investigate the occurance of thermal transpiration in copper and platinum nanochannels at atmospheric pressure conditions. The mean pressure of argon gas confined inside the nano channels was maintained around 1 bar. The channel height is maintained at 2nm. The argon atoms interact with each other and with the wall atoms through the Lennard-Jones potential. The wall atoms are modelled using an EAM potential. Further, separate simulations were carried out where a Harmonic potential is used for the atom-atom interaction in the platinum channel. A thermally insulating wall was introduced between the low and high temperature regions and those wall atoms interact with fluid atoms through a repulsive potential. A reduced cut off radius were used to achieve this. Thermal creep is induced by applying a temperature gradient along the channel wall. It was found that flow developed in the direction of the increasing temperature gradient of the wall. An increase in the volumetric flux was observed as the length of the cold and the hot regions of the wall were increased. The effect of temperature gradient and the wall-fluid interaction strength on the flow parameters have been studied to understand the phenomenon better.
Nanoscale deicing by molecular dynamics simulation
NASA Astrophysics Data System (ADS)
Xiao, Senbo; He, Jianying; Zhang, Zhiliang
2016-07-01
Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice adhesion strength by an aqueous water layer, and provide atomistic details that support previous experimental studies. Our results contribute quantitative comparison of nanoscale adhesion strength of ice on hydrophobic and hydrophilic surfaces, and supply for the first time theoretical references for understanding the mechanics at the atomistic origins of macroscale ice adhesion.Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice
How Dynamic Visualization Technology Can Support Molecular Reasoning
ERIC Educational Resources Information Center
Levy, Dalit
2013-01-01
This paper reports the results of a study aimed at exploring the advantages of dynamic visualization for the development of better understanding of molecular processes. We designed a technology-enhanced curriculum module in which high school chemistry students conduct virtual experiments with dynamic molecular visualizations of solid, liquid, and…
Ruggiero, Michael T; Gooch, Jonathan; Zubieta, Jon; Korter, Timothy M
2016-02-18
The problem of nonlocal interactions in density functional theory calculations has in part been mitigated by the introduction of range-corrected functional methods. While promising solutions, the continued evaluation of range corrections in the structural simulations of complex molecular crystals is required to judge their efficacy in challenging chemical environments. Here, three pyridinium-based crystals, exhibiting a wide range of intramolecular and intermolecular interactions, are used as benchmark systems for gauging the accuracy of several range-corrected density functional techniques. The computational results are compared to low-temperature experimental single-crystal X-ray diffraction and terahertz spectroscopic measurements, enabling the direct assessment of range correction in the accurate simulation of the potential energy surface minima and curvatures. Ultimately, the simultaneous treatment of both short- and long-range effects by the ωB97-X functional was found to be central to its rank as the top performer in reproducing the complex array of forces that occur in the studied pyridinium solids. These results demonstrate that while long-range corrections are the most commonly implemented range-dependent improvements to density functionals, short-range corrections are vital for the accurate reproduction of forces that rapidly diminish with distance, such as quadrupole-quadrupole interactions. PMID:26814572
Emergence of spacetime dynamics in entropy corrected and braneworld models
Sheykhi, A.; Dehghani, M.H.; Hosseini, S.E. E-mail: mhd@shirazu.ac.ir
2013-04-01
A very interesting new proposal on the origin of the cosmic expansion was recently suggested by Padmanabhan [arXiv:1206.4916]. He argued that the difference between the surface degrees of freedom and the bulk degrees of freedom in a region of space drives the accelerated expansion of the universe, as well as the standard Friedmann equation through relation ΔV = Δt(N{sub sur}−N{sub bulk}). In this paper, we first present the general expression for the number of degrees of freedom on the holographic surface, N{sub sur}, using the general entropy corrected formula S = A/(4L{sub p}{sup 2})+s(A). Then, as two example, by applying the Padmanabhan's idea we extract the corresponding Friedmann equations in the presence of power-law and logarithmic correction terms in the entropy. We also extend the study to RS II and DGP braneworld models and derive successfully the correct form of the Friedmann equations in these theories. Our study further supports the viability of Padmanabhan's proposal.
Statistical coarse-graining of molecular dynamics into peridynamics.
Silling, Stewart Andrew; Lehoucq, Richard B.
2007-10-01
This paper describes an elegant statistical coarse-graining of molecular dynamics at finite temperature into peridynamics, a continuum theory. Peridynamics is an efficient alternative to molecular dynamics enabling dynamics at larger length and time scales. In direct analogy with molecular dynamics, peridynamics uses a nonlocal model of force and does not employ stress/strain relationships germane to classical continuum mechanics. In contrast with classical continuum mechanics, the peridynamic representation of a system of linear springs and masses is shown to have the same dispersion relation as the original spring-mass system.
Constant pressure and temperature discrete-time Langevin molecular dynamics.
Grønbech-Jensen, Niels; Farago, Oded
2014-11-21
We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems-a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation. PMID:25416875
Constant pressure and temperature discrete-time Langevin molecular dynamics
Grønbech-Jensen, Niels; Farago, Oded
2014-11-21
We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems—a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.
Molecular Dynamics Simulations of Coulomb Explosion
Bringa, E M
2002-05-17
A swift ion creates a track of electronic excitations in the target material. A net repulsion inside the track can cause a ''Coulomb Explosion'', which can lead to damage and sputtering of the material. Here we report results from molecular-dynamics (MD) simulations of Coulomb explosion for a cylindrical track as a function of charge density and neutralization/quenching time, {tau}. Screening by the free electrons is accounted for using a screened Coulomb potential for the interaction among charges. The yield exhibits a prompt component from the track core and a component, which dominates at higher excitation density, from the heated region produced. For the cases studied, the number of atoms ejected per incident ion, i.e. the sputtering yield Y, is quadratic with charge density along the track as suggested by simple models. Y({tau} = 0.2 Debye periods) is nearly 20% of the yield when there is no neutralization ({tau} {yields} {infinity}). The connections between ''Coulomb explosions'', thermal spikes and measurements of electronic sputtering are discussed.
A molecular dynamics study of dielectric friction
Kurnikova, M.G.; Waldeck, D.H.; Coalson, R.D.
1996-07-01
A molecular dynamics study of the friction experienced by the dye molecule resorufamine rotating in a polar solvent is performed. The validity of simple continuum theories of dielectric friction is tested. It is found that the Alavi{endash}Waldeck theory gives reasonable results for the zero frequency dielectric friction coefficient while the Nee{endash}Zwanzig theory requires an unphysically small cavity radius. A procedure for evaluating the time dependent friction kernel from torques and angular velocities, which enables the contributions to the friction from the van der Waals and Coulomb forces to be evaluated separately, is suggested. This study of a realistic system shows that electrostatic interactions can enhance friction by at least two physical mechanisms. First is a contribution to the friction which arises solely from retardation of the solvent reaction field. Second is a contribution arising from local structural changes of the solvent which are driven by the electrostatic field, i.e., a change in the local viscosity. {copyright} {ital 1996 American Institute of Physics.}
Quantum molecular dynamics simulations of dense matter
Collins, L.; Kress, J.; Troullier, N.; Lenosky, T.; Kwon, I.
1997-12-31
The authors have developed a quantum molecular dynamics (QMD) simulation method for investigating the properties of dense matter in a variety of environments. The technique treats a periodically-replicated reference cell containing N atoms in which the nuclei move according to the classical equations-of-motion. The interatomic forces are generated from the quantum mechanical interactions of the (between?) electrons and nuclei. To generate these forces, the authors employ several methods of varying sophistication from the tight-binding (TB) to elaborate density functional (DF) schemes. In the latter case, lengthy simulations on the order of 200 atoms are routinely performed, while for the TB, which requires no self-consistency, upwards to 1000 atoms are systematically treated. The QMD method has been applied to a variety cases: (1) fluid/plasma Hydrogen from liquid density to 20 times volume-compressed for temperatures of a thousand to a million degrees Kelvin; (2) isotopic hydrogenic mixtures, (3) liquid metals (Li, Na, K); (4) impurities such as Argon in dense hydrogen plasmas; and (5) metal/insulator transitions in rare gas systems (Ar,Kr) under high compressions. The advent of parallel versions of the methods, especially for fast eigensolvers, presage LDA simulations in the range of 500--1000 atoms and TB runs for tens of thousands of particles. This leap should allow treatment of shock chemistry as well as large-scale mixtures of species in highly transient environments.
Integrating influenza antigenic dynamics with molecular evolution
Bedford, Trevor; Suchard, Marc A; Lemey, Philippe; Dudas, Gytis; Gregory, Victoria; Hay, Alan J; McCauley, John W; Russell, Colin A; Smith, Derek J; Rambaut, Andrew
2014-01-01
Influenza viruses undergo continual antigenic evolution allowing mutant viruses to evade host immunity acquired to previous virus strains. Antigenic phenotype is often assessed through pairwise measurement of cross-reactivity between influenza strains using the hemagglutination inhibition (HI) assay. Here, we extend previous approaches to antigenic cartography, and simultaneously characterize antigenic and genetic evolution by modeling the diffusion of antigenic phenotype over a shared virus phylogeny. Using HI data from influenza lineages A/H3N2, A/H1N1, B/Victoria and B/Yamagata, we determine patterns of antigenic drift across viral lineages, showing that A/H3N2 evolves faster and in a more punctuated fashion than other influenza lineages. We also show that year-to-year antigenic drift appears to drive incidence patterns within each influenza lineage. This work makes possible substantial future advances in investigating the dynamics of influenza and other antigenically-variable pathogens by providing a model that intimately combines molecular and antigenic evolution. DOI: http://dx.doi.org/10.7554/eLife.01914.001 PMID:24497547
Molecular chaperone-mediated nuclear protein dynamics.
Echtenkamp, Frank J; Freeman, Brian C
2014-05-01
Homeostasis requires effective action of numerous biological pathways including those working along a genome. The variety of processes functioning in the nucleus is considerable, yet the number of employed factors eclipses this total. Ideally, individual components assemble into distinct complexes and serially operate along a pathway to perform work. Adding to the complexity is a multitude of fluctuating internal and external signals that must be monitored to initiate, continue or halt individual activities. While cooperative interactions between proteins of the same process provide a mechanism for rapid and precise assembly, the inherent stability of such organized structures interferes with the proper timing of biological events. Further prolonging the longevity of biological complexes are crowding effects resulting from the high concentration of intracellular macromolecules. Hence, accessory proteins are required to destabilize the various assemblies to efficiently transition between structures, avoid off-pathway competitive interactions, and to terminate pathway activity. We suggest that molecular chaperones have evolved, in part, to manage these challenges by fostering a general and continuous dynamic protein environment within the nucleus. PMID:24694369
Nanoscale deicing by molecular dynamics simulation.
Xiao, Senbo; He, Jianying; Zhang, Zhiliang
2016-08-14
Deicing is important to human activities in low-temperature circumstances, and is critical for combating the damage caused by excessive accumulation of ice. The aim of creating anti-icing materials, surfaces and applications relies on the understanding of fundamental nanoscale ice adhesion mechanics. Here in this study, we employ all-atom modeling and molecular dynamics simulation to investigate ice adhesion. We apply force to detach and shear nano-sized ice cubes for probing the determinants of atomistic adhesion mechanics, and at the same time investigate the mechanical effect of a sandwiched aqueous water layer between ice and substrates. We observe that high interfacial energy restricts ice mobility and increases both ice detaching and shearing stresses. We quantify up to a 60% decrease in ice adhesion strength by an aqueous water layer, and provide atomistic details that support previous experimental studies. Our results contribute quantitative comparison of nanoscale adhesion strength of ice on hydrophobic and hydrophilic surfaces, and supply for the first time theoretical references for understanding the mechanics at the atomistic origins of macroscale ice adhesion. PMID:27431975
NASA Astrophysics Data System (ADS)
Brennan, John K.; Lísal, Martin; Gubbins, Keith E.; Rice, Betsy M.
2004-12-01
A molecular simulation method to study the dynamics of chemically reacting mixtures is presented. The method uses a combination of stochastic and dynamic simulation steps, allowing for the simulation of both thermodynamic and transport properties. The method couples a molecular dynamics simulation cell (termed dynamic cell) to a reaction mixture simulation cell (termed control cell) that is formulated upon the reaction ensemble Monte Carlo (RxMC) method, hence the term reaction ensemble molecular dynamics. Thermodynamic and transport properties are calculated in the dynamic cell by using a constant-temperature molecular dynamics simulation method. RxMC forward and reverse reaction steps are performed in the control cell only, while molecular dynamics steps are performed in both the dynamic cell and the control cell. The control cell, which acts as a sink and source reservoir, is maintained at reaction equilibrium conditions via the RxMC algorithm. The reaction ensemble molecular dynamics method is analogous to the grand canonical ensemble molecular dynamics technique, while using some elements of the osmotic molecular dynamics method, and so simulates conditions that directly relate to real, open systems. The accuracy and stability of the method is assessed by considering the ammonia synthesis reaction N2+3H2⇔2NH3 . It is shown to be a viable method for predicting the effects of nonideal environments on the dynamic properties (particularly diffusion) as well as reaction equilibria for chemically reacting mixtures.
Dynamic stresses, coulomb failure, and remote triggering: corrected
Hill, David P.
2012-01-01
Dynamic stresses associated with crustal surface waves with 15–30 s periods and peak amplitudes <1 MPa are capable of triggering seismicity at sites remote from the generating mainshock under appropriate conditions. Coulomb failure models based on a frictional strength threshold offer one explanation for instances of rapid‐onset triggered seismicity that develop during the surface‐wave peak dynamic stressing. Evaluation of the triggering potential of surface‐wave dynamic stresses acting on critically stressed faults using a Mohr’s circle representation together with the Coulomb failure criteria indicates that Love waves should have a higher triggering potential than Rayleigh waves for most fault orientations and wave incidence angles. That (1) the onset of triggered seismicity often appears to begin during the Rayleigh wave rather than the earlier arriving Love wave, and (2) Love‐wave amplitudes typically exceed those for Rayleigh waves suggests that the explanation for rapid‐onset dynamic triggering may not reside solely with a simple static‐threshold friction mode. The results also indicate that normal faults should be more susceptible to dynamic triggering by 20‐s Rayleigh‐wave stresses than thrust faults in the shallow seismogenic crust (<10 km) while the advantage tips in favor of reverse faults greater depths. This transition depth scales with wavelength and coincides roughly with the transition from retrograde‐to‐prograde particle motion. Locally elevated pore pressures may have a role in the observed prevalence of dynamic triggering in extensional regimes and geothermal/volcanic systems. The result is consistent with the apparent elevated susceptibility of extensional or transtensional tectonic regimes to remote triggering by Rayleigh‐wave dynamic stresses than compressional or transpressional regimes.
Molecular dynamics in cytochrome c oxidase Moessbauer spectra deconvolution
Bossis, Fabrizio; Palese, Luigi L.
2011-01-07
Research highlights: {yields} Cytochrome c oxidase molecular dynamics serve to predict Moessbauer lineshape widths. {yields} Half height widths are used in modeling of Lorentzian doublets. {yields} Such spectral deconvolutions are useful in detecting the enzyme intermediates. -- Abstract: In this work low temperature molecular dynamics simulations of cytochrome c oxidase are used to predict an experimentally observable, namely Moessbauer spectra width. Predicted lineshapes are used to model Lorentzian doublets, with which published cytochrome c oxidase Moessbauer spectra were simulated. Molecular dynamics imposed constraints to spectral lineshapes permit to obtain useful information, like the presence of multiple chemical species in the binuclear center of cytochrome c oxidase. Moreover, a benchmark of quality for molecular dynamic simulations can be obtained. Despite the overwhelming importance of dynamics in electron-proton transfer systems, limited work has been devoted to unravel how much realistic are molecular dynamics simulations results. In this work, molecular dynamics based predictions are found to be in good agreement with published experimental spectra, showing that we can confidently rely on actual simulations. Molecular dynamics based deconvolution of Moessbauer spectra will lead to a renewed interest for application of this approach in bioenergetics.
Free energy calculations using dual-level Born-Oppenheimer molecular dynamics
NASA Astrophysics Data System (ADS)
Retegan, Marius; Martins-Costa, Marilia; Ruiz-López, Manuel F.
2010-08-01
We describe an efficient and accurate method to compute free energy changes in complex chemical systems that cannot be described through classical molecular dynamics simulations, examples of which are chemical and photochemical reactions in solution, enzymes, interfaces, etc. It is based on the use of dual-level Born-Oppenheimer molecular dynamics simulations. A low-level quantum mechanical method is employed to calculate the potential of mean force through the umbrella sampling technique. Then, a high-level quantum mechanical method is used to estimate a free energy correction on selected points of the reaction coordinate using perturbation theory. The precision of the results is comparable to that of ab initio molecular dynamics methods such as the Car-Parrinello approach but the computational cost is much lower, roughly by two to three orders of magnitude. The method is illustrated by discussing the association free energy of simple organometallic compounds, although the field of application is very broad.
Combined molecular dynamics-spin dynamics simulations of bcc iron
Perera, Meewanage Dilina N; Yin, Junqi; Landau, David P; Nicholson, Don M; Stocks, George Malcolm; Eisenbach, Markus; Brown, Greg
2014-01-01
Using a classical model that treats translational and spin degrees of freedom on an equal footing, we study phonon-magnon interactions in BCC iron with combined molecular and spin dynamics methods. The atomic interactions are modeled via an empirical many-body potential while spin dependent interactions are established through a Hamiltonian of the Heisenberg form with a distance dependent magnetic exchange interaction obtained from first principles electronic structure calculations. The temporal evolution of translational and spin degrees of freedom was determined by numerically solving the coupled equations of motion, using an algorithm based on the second order Suzuki-Trotter decomposition of the exponential operators. By calculating Fourier transforms of space- and time-displaced correlation functions, we demonstrate that the the presence of lattice vibrations leads to noticeable softening and damping of spin wave modes. As a result of the interplay between lattice and spin subsystems, we also observe additional longitudinal spin wave excitations, with frequencies which coincide with that of the longitudinal lattice vibrations.
Singh, Surya Pratap; Gupta, Dwijendra K
2016-04-01
Casein kinase-1 (CK1) isoforms actively participate in the down-regulation of canonical Wnt signaling pathway; however recent studies have shown their active roles in oncogenesis of various tissues through this pathway. Functional loss of two isoforms (CK1-α/ε) has been shown to activate the carcinogenic pathway which involves the stabilization of of cytoplasmic β-catenin. Development of anticancer therapeutics is very laborious task and depends upon the structural and conformational details of the target. This study focuses on, how the structural dynamics and conformational changes of two CK1 isoforms are synchronized in carcinogenic pathway. The conformational dynamics in kinases is the responsible for their action as has been supported by the molecular docking experiments. PMID:26788877
Molecular Dynamics Simulation of Disordered Zircon
Devanathan, Ram; Corrales, Louis R.; Weber, William J.; Chartier, Alain; Meis, Constantin
2004-02-27
The melting of zircon and the amorphous state produced by quenching from the melt were simulated by molecular dynamics using a new partial charge model combined with the Ziegler-Biersack-Littmark potential. The model has been established for the description of the crystalline and aperiodic structures of zircon in order to be used for the simulation of displacement cascades. It provides an excellent fit to the structure, and accounts with convenient precision the mechanical and thermodynamic properties of zircon. The calculated melting temperature is about 2100 K. The activation energy for self-diffusion of ions in the liquid state was determined to be 190-200 kJ/mole. Melt quenching was employed to produce two different disordered states with distinct densities and structures. In the high density disordered state, the zircon structure is intact but the bond angle distributions are broader, 4% of the Si units are polymerized, and the volume swelling is about 8%. In the low density amorphous state, the Zr and Si coordination numbers are lower, and the Zr-O and Si-O bond lengths are shorter than corresponding values for the crystal. In addition, a highly polymerized Si network, with average connectivity of two, is observed in the low density amorphous state. These features have all been experimentally observed in natural metamict zircon. The present findings, when considered in light of experimental radiation effects studies, suggest that the swelling in zircon arises initially from disorder in the zircon crystal, and at high doses the disordered crystal is unable to accommodate the volume expansion and transforms to the amorphous state.
Condensation on nanorods by molecular dynamics
NASA Astrophysics Data System (ADS)
Suh, Donguk; Yasuoka, Kenji
2016-06-01
Many recent experimental studies have been conducted on constructing nanorods and nanowires to use in a wide range of applications. In this study, molecular dynamics is used to directly examine the condensation rate of nanorods and the results are compared with other basic configurations such as cubes or spheres. According to previous studies conducted by Suh and Yasuoka [J. Phys. Chem. B 115, 10631 (2011); 116, 14637 (2012)], a simple change in the configuration of the seed produces a shape effect, where the curvature of the solid seed surface directly affects the growth generating an orderly difference depending on the curvature. Nanoscale cuboids or nanorods were studied to find an aspect ratio effect when condensation occurs on the surface. Various aspect ratios were examined for different nanorod sizes over a wide range of supersaturation ratios. The results show that the growth rate of the nanorod is independent of the supersaturation ratio, which was also observed for the sphere and cube. The growth rate for the rod fell between those of the cube and the sphere, and this is due to an increase in the surface area of the nanorod compared to the cube and curvature effect in comparison with the sphere. A clear size dependence of the seed was observed, which is also similar to the cube and sphere. Furthermore, no aspect ratio influence was seen for the growth rate. This does not mean that the actual amount of condensation is the same for longer seeds, but rather from the definition of the growth rate, the amount of accumulation per unit area is the same for all seed lengths.
NASA Astrophysics Data System (ADS)
Masunov, Artëm E.; Gangopadhyay, Shruba
2015-12-01
New method to eliminate the spin-contamination in broken symmetry density functional theory (BS DFT) calculations is introduced. Unlike conventional spin-purification correction, this method is based on canonical Natural Orbitals (NO) for each high/low spin coupled electron pair. We derive an expression to extract the energy of the pure singlet state given in terms of energy of BS DFT solution, the occupation number of the bonding NO, and the energy of the higher spin state built on these bonding and antibonding NOs (not self-consistent Kohn-Sham orbitals of the high spin state). Compared to the other spin-contamination correction schemes, spin-correction is applied to each correlated electron pair individually. We investigate two binuclear Mn(IV) molecular magnets using this pairwise correction. While one of the molecules is described by magnetic orbitals strongly localized on the metal centers, and spin gap is accurately predicted by Noodleman and Yamaguchi schemes, for the other one the gap is predicted poorly by these schemes due to strong delocalization of the magnetic orbitals onto the ligands. We show our new correction to yield more accurate results in both cases.
Accelerating MP2C dispersion corrections for dimers and molecular crystals
NASA Astrophysics Data System (ADS)
Huang, Yuanhang; Shao, Yihan; Beran, Gregory J. O.
2013-06-01
The MP2C dispersion correction of Pitonak and Hesselmann [J. Chem. Theory Comput. 6, 168 (2010)], 10.1021/ct9005882 substantially improves the performance of second-order Møller-Plesset perturbation theory for non-covalent interactions, albeit with non-trivial computational cost. Here, the MP2C correction is computed in a monomer-centered basis instead of a dimer-centered one. When applied to a single dimer MP2 calculation, this change accelerates the MP2C dispersion correction several-fold while introducing only trivial new errors. More significantly, in the context of fragment-based molecular crystal studies, combination of the new monomer basis algorithm and the periodic symmetry of the crystal reduces the cost of computing the dispersion correction by two orders of magnitude. This speed-up reduces the MP2C dispersion correction calculation from a significant computational expense to a negligible one in crystals like aspirin or oxalyl dihydrazide, without compromising accuracy.
Beam spot shift in a dynamic astigmatism correction type (DQ-DAF) electron gun
NASA Astrophysics Data System (ADS)
Ueda, Y.; Sukeno, M.; Suzuki, H.
This paper presents the results of an investigation to avoid the convergence error in color CRT guns operating with a dynamic focus correction. The dynamic spot shift as well as the spot coma can be avoided at the same time with a proper arrangement of the quadrupole and main lenses with an oblique incidence of the side beam. The design confirmation is also presented.
NASA Astrophysics Data System (ADS)
Grest, Gary S.
2008-03-01
Twenty years ago at the APS March Meeting, Kurt Kremer and I presented the first numerical evidence from computer simulations that the reptation model of Edwards and de Gennes correctly describes the dynamics of entangled linear polymer melts. For chains longer than the entanglement length Ne, the monomers of a chain move predominantly along their own contour. The distinctive signature of reptation dynamics, which we observed, was that on intermediate time scales, the mean squared displacement of a monomer increases with time as t^ 1/4. Though these early simulations were limited to chains of a few Ne, they demonstrated the potential of computer simulations to contribute to our understanding of polymer dynamics. Here I will review the progress over the past twenty years and present an outlook for the future in modeling entangled polymer melts and networks. With present day computers coupled with efficient parallel molecular dynamics codes, it is now possible to follow the equilibrium dynamics of chains of length 10-20Ne from the early Rouse regime to the long time diffusive regime. Result of these simulations support the earlier results obtained on chains of only a few Ne. Further evidence for the tube models of polymer dynamics has been obtained by identifying the primitive path mesh that characterizes the microscopic topological state of the computer- generated conformations of the chains. In particular, the plateau moduli derived on the basis of this analysis quantitatively reproduce experimental data for a wide spectrum of entangled polymer liquids including semi-dilute theta solutions of synthetic polymers, the corresponding dense melts, and solutions of semi-flexible (bio)polymers such as f-actin or suspensions of rodlike viruses. We also find that in agreement with the reptation model, the stress, end-to-end distance and entanglement length of an entangled melt subjected to uniaxial elongation, all relax on the same time scale.
Diamond, Jonathan S; Nashed, Joseph Y; Johansson, Roland S; Wolpert, Daniel M; Flanagan, J Randall
2015-07-22
Numerous studies have shown that people are adept at learning novel object dynamics, linking applied force and motion, when performing reaching movements with hand-held objects. Here we investigated whether the control of rapid corrective arm responses, elicited in response to visual perturbations, has access to such newly acquired knowledge of object dynamics. Participants first learned to make reaching movements while grasping an object subjected to complex load forces that depended on the distance and angle of the hand from the start position. During a subsequent test phase, we examined grip and load force coordination during corrective arm movements elicited (within ∼150 ms) in response to viewed sudden lateral shifts (1.5 cm) in target or object position. We hypothesized that, if knowledge of object dynamics is incorporated in the control of the corrective responses, grip force changes would anticipate the unusual load force changes associated with the corrective arm movements so as to support grasp stability. Indeed, we found that the participants generated grip force adjustments tightly coupled, both spatially and temporally, to the load force changes associated with the arm movement corrections. We submit that recently learned novel object dynamics are effectively integrated into sensorimotor control policies that support rapid visually driven arm corrective actions during transport of hand held objects. Significance statement: Previous studies have demonstrated that the motor system can learn, and make use of, internal models of object dynamics to generate feedforward motor commands. However, it is not known whether such internal models are incorporated into rapid, automatic arm movement corrections that compensate for errors that arise during movement. Here we demonstrate, for the first time, that internal models of novel object dynamics are integrated into rapid corrective arm movements made in response to visuomotor perturbations that, importantly, do
Las Palmeras Molecular Dynamics: A flexible and modular molecular dynamics code
NASA Astrophysics Data System (ADS)
Davis, Sergio; Loyola, Claudia; González, Felipe; Peralta, Joaquín
2010-12-01
Las Palmeras Molecular Dynamics (LPMD) is a highly modular and extensible molecular dynamics (MD) code using interatomic potential functions. LPMD is able to perform equilibrium MD simulations of bulk crystalline solids, amorphous solids and liquids, as well as non-equilibrium MD (NEMD) simulations such as shock wave propagation, projectile impacts, cluster collisions, shearing, deformation under load, heat conduction, heterogeneous melting, among others, which involve unusual MD features like non-moving atoms and walls, unstoppable atoms with constant-velocity, and external forces like electric fields. LPMD is written in C++ as a compromise between efficiency and clarity of design, and its architecture is based on separate components or plug-ins, implemented as modules which are loaded on demand at runtime. The advantage of this architecture is the ability to completely link together the desired components involved in the simulation in different ways at runtime, using a user-friendly control file language which describes the simulation work-flow. As an added bonus, the plug-in API (Application Programming Interface) makes it possible to use the LPMD components to analyze data coming from other simulation packages, convert between input file formats, apply different transformations to saved MD atomic trajectories, and visualize dynamical processes either in real-time or as a post-processing step. Individual components, such as a new potential function, a new integrator, a new file format, new properties to calculate, new real-time visualizers, and even a new algorithm for handling neighbor lists can be easily coded, compiled and tested within LPMD by virtue of its object-oriented API, without the need to modify the rest of the code. LPMD includes already several pair potential functions such as Lennard-Jones, Morse, Buckingham, MCY and the harmonic potential, as well as embedded-atom model (EAM) functions such as the Sutton-Chen and Gupta potentials. Integrators to
Molecular density functional theory for water with liquid-gas coexistence and correct pressure
Jeanmairet, Guillaume Levesque, Maximilien; Sergiievskyi, Volodymyr; Borgis, Daniel
2015-04-21
The solvation of hydrophobic solutes in water is special because liquid and gas are almost at coexistence. In the common hypernetted chain approximation to integral equations, or equivalently in the homogenous reference fluid of molecular density functional theory, coexistence is not taken into account. Hydration structures and energies of nanometer-scale hydrophobic solutes are thus incorrect. In this article, we propose a bridge functional that corrects this thermodynamic inconsistency by introducing a metastable gas phase for the homogeneous solvent. We show how this can be done by a third order expansion of the functional around the bulk liquid density that imposes the right pressure and the correct second order derivatives. Although this theory is not limited to water, we apply it to study hydrophobic solvation in water at room temperature and pressure and compare the results to all-atom simulations. The solvation free energy of small molecular solutes like n-alkanes and hard sphere solutes whose radii range from angstroms to nanometers is now in quantitative agreement with reference all atom simulations. The macroscopic liquid-gas surface tension predicted by the theory is comparable to experiments. This theory gives an alternative to the empirical hard sphere bridge correction used so far by several authors.
MOLECULAR DYNAMICS STUDIES OF MODULAR POLYKETIDE SYNTHASE KETOREDUCTASE STEREOSPECIFICITY
Mugnai, Mauro L.; Shi, Yue; Keatinge-Clay, Adrian T.; Elber, Ron
2015-01-01
Ketoreductases (KRs) from modular polyketide synthases (PKSs) can perform stereospecific catalysis, selecting a polyketide with a D-α or an L-α-methyl substituent for NADPH-mediated reduction. In this report Molecular Dynamics (MD) simulations were performed to investigate the interactions that control stereospecificity. We studied the A1-type KR from the second module of the amphotericin PKS (A1), which is known to be stereospecific for a D-α-methyl-substituted diketide substrate (dkD). MD simulations of two ternary complexes comprised of the enzyme, NADPH, and either the correct substrate, dkD, or its enantiomer (dkL) were performed. The coordinates for the A1/NADPH binary complex were obtained from a crystal structure [Zheng, J. T. et al. (2010) Structure 18, 913–922], and substrates were modeled in the binding pocket in conformations appropriate for reduction. Simulations were intended to reproduce the initial weak binding of the polyketide substrate to the enzyme. Long (tens of nanoseconds) MD simulations show that the correct substrate is retained in a conformation closer to the reactive configuration. Many short (up to a nanosecond) MD runs starting from the initial structures display evidence that Q364, three residues N-terminal to the catalytic tyrosine, forms a hydrogen bond to the incorrect dkL substrate to yield an unreactive conformation that is more favorable than the reactive configuration. This interaction is not as strong for dkD, as the D-α-methyl substituent is positioned between the glutamine and the reactive site. This result correlates with experimental findings [Zheng, J. T. et al. (2010) Structure 18, 913–922] in which a Q364H mutant was observed to lose stereospecificity. PMID:25835227
NASA Astrophysics Data System (ADS)
Kritayakornupong, Chinapong; Hannongbua, Supot
2007-01-01
The structural and dynamical properties of high-spin Ru 2+ in aqueous solution have been theoretically studied using molecular dynamics (MD) simulations. The conventional MD simulation based on pair potentials gives the overestimated average first shell coordination number of 9, whereas the value of 5.9 was observed when the three-body corrected function was included. A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to take into account the many-body effects on the hydration shell structure of Ru 2+. The most important region, the first hydration shell, was treated by ab initio quantum mechanics at UHF level using the SBKJC VDZ ECP basis set for Ru 2+ and the 6-31G ∗ basis sets for water. An exact coordination number of 6 for the first hydration shell was obtained from the QM/MM simulation. The QM/MM simulation predicts the average Ru 2+-O distance of 2.42 Å for the first hydration shell, whereas the values of 2.34 and 2.46 Å are resulted from the pair potentials without and with the three-body corrected simulations, respectively. Several other structural properties representing position and orientation of the solvate molecules were evaluated for describing the hydration shell structure of the Ru 2+ ion in dilute aqueous solution. A mean residence time of 7.1 ps was obtained for water ligands residing in the second hydration shell.
Kronik, Leeor; Tkatchenko, Alexandre
2014-11-18
CONSPECTUS: Molecular crystals are ubiquitous in many areas of science and engineering, including biology and medicine. Until recently, our ability to understand and predict their structure and properties using density functional theory was severely limited by the lack of approximate exchange-correlation functionals able to achieve sufficient accuracy. Here we show that there are many cases where the simple, minimally empirical pairwise correction scheme of Tkatchenko and Scheffler provides a useful prediction of the structure and properties of molecular crystals. After a brief introduction of the approach, we demonstrate its strength through some examples taken from our recent work. First, we show the accuracy of the approach using benchmark data sets of molecular complexes. Then we show its efficacy for structural determination using the hemozoin crystal, a challenging system possessing a wide range of strong and weak binding scenarios. Next, we show that it is equally useful for response properties by considering the elastic constants exhibited by the supramolecular diphenylalanine peptide solid and the infrared signature of water libration movements in brushite. Throughout, we emphasize lessons learned not only for the methodology but also for the chemistry and physics of the crystals in question. We further show that in many other scenarios where the simple pairwise correction scheme is not sufficiently accurate, one can go beyond it by employing a computationally inexpensive many-body dispersive approach that results in useful, quantitative accuracy, even in the presence of significant screening and/or multibody contributions to the dispersive energy. We explain the principles of the many-body approach and demonstrate its accuracy for benchmark data sets of small and large molecular complexes and molecular solids. PMID:24901508
CHARACTERIZING COUPLED CHARGE TRANSPORT WITH MULTISCALE MOLECULAR DYNAMICS
Swanson, Jessica
2011-08-31
This is the final progress report for Award DE-SC0004920, entitled 'Characterizing coupled charge transport with multi scale molecular dynamics'. The technical abstract will be provided in the uploaded report.
Masses, luminosities and dynamics of galactic molecular clouds
NASA Technical Reports Server (NTRS)
Solomon, P. M.; Rivolo, A. R.; Mooney, T. J.; Barrett, J. W.; Sage, L. J.
1987-01-01
Star formation in galaxies takes place in molecular clouds and the Milky Way is the only galaxy in which it is possible to resolve and study the physical properties and star formation activity of individual clouds. The masses, luminosities, dynamics, and distribution of molecular clouds, primarily giant molecular clouds in the Milky Way are described and analyzed. The observational data sets are the Massachusetts-Stony Brook CO Galactic Plane Survey and the IRAS far IR images. The molecular mass and infrared luminosities of glactic clouds are then compared with the molecular mass and infrared luminosities of external galaxies.
Thermostat artifacts in replica exchange molecular dynamics simulations.
Rosta, Edina; Buchete, Nicolae-Viorel; Hummer, Gerhard
2009-01-01
We explore the effects of thermostats in replica exchange molecular dynamics (REMD) simulations. For thermostats that do not produce a canonical ensemble, REMD simulations are found to distort the configuration-space distributions. For bulk water, we find small deviations of the average potential energies, the buildup of tails in the potential energy distributions, and artificial correlations between the energies at different temperatures. If a solute is present, as in protein folding simulations, its conformational equilibrium can be altered. In REMD simulations of a helix-forming peptide with a weak-coupling (Berendsen) thermostat, we find that the folded state is overpopulated by about 10% at low temperatures, and underpopulated at high temperatures. As a consequence, the enthalpy of folding deviates by almost 3 kcal/mol from the correct value. The reason for this population shift is that non-canonical ensembles with narrowed potential energy fluctuations artificially bias toward replica exchanges between low-energy folded structures at the high temperature and high-energy unfolded structures at the low temperature. We conclude that REMD simulations should only be performed in conjunction with thermostats that produce a canonical ensemble. PMID:20046980
Thermostat artifacts in replica exchange molecular dynamics simulations
Rosta, Edina; Buchete, Nicolae-Viorel; Hummer, Gerhard
2009-01-01
We explore the effects of thermostats in replica exchange molecular dynamics (REMD) simulations. For thermostats that do not produce a canonical ensemble, REMD simulations are found to distort the configuration-space distributions. For bulk water, we find small deviations of the average potential energies, the buildup of tails in the potential energy distributions, and artificial correlations between the energies at different temperatures. If a solute is present, as in protein folding simulations, its conformational equilibrium can be altered. In REMD simulations of a helix-forming peptide with a weak-coupling (Berendsen) thermostat, we find that the folded state is overpopulated by about 10% at low temperatures, and underpopulated at high temperatures. As a consequence, the enthalpy of folding deviates by almost 3 kcal/mol from the correct value. The reason for this population shift is that non-canonical ensembles with narrowed potential energy fluctuations artificially bias toward replica exchanges between low-energy folded structures at the high temperature and high-energy unfolded structures at the low temperature. We conclude that REMD simulations should only be performed in conjunction with thermostats that produce a canonical ensemble. PMID:20046980
Hydration structure of salt solutions from ab initio molecular dynamics.
Bankura, Arindam; Carnevale, Vincenzo; Klein, Michael L
2013-01-01
The solvation structures of Na(+), K(+), and Cl(-) ions in aqueous solution have been investigated using density functional theory (DFT) based Car-Parrinello (CP) molecular dynamics (MD) simulations. CPMD trajectories were collected for systems containing three NaCl or KCl ion pairs solvated by 122 water molecules using three different but commonly employed density functionals (BLYP, HCTH, and PBE) with electron correlation treated at the level of the generalized gradient approximation (GGA). The effect of including dispersion forces was analyzed through the use of an empirical correction to the DFT-GGA scheme. Special attention was paid to the hydration characteristics, especially the structural properties of the first solvation shell of the ions, which was investigated through ion-water radial distribution functions, coordination numbers, and angular distribution functions. There are significant differences between the present results obtained from CPMD simulations and those provided by classical MD based on either the CHARMM force field or a polarizable model. Overall, the computed structural properties are in fair agreement with the available experimental results. In particular, the observed coordination numbers 5.0-5.5, 6.0-6.4, and 6.0-6.5 for Na(+), K(+), and Cl(-), respectively, are consistent with X-ray and neutron scattering studies but differ somewhat from some of the many other recent computational studies of these important systems. Possible reasons for the differences are discussed. PMID:23298049
Nuclear structure and reactions in the fermionic molecular dynamics approach
NASA Astrophysics Data System (ADS)
Neff, T.; Feldmeier, H.
2008-05-01
The Fermionic Molecular Dynamics (FMD) model uses Gaussian wave packets as single-particle states. Intrinsic many-body basis states are constructed as Slater determinants which have to be projected on parity, angular momentum and total linear momentum to restore the symmetries of the Hamiltonian. The flexibility of the Gaussian basis allows to economically describe states with shell structures as well as states featuring clustering or halos. We use an effective interaction that is derived from the realistic Argonne V18 interaction by means of the Unitary Correlation Operator Method (UCOM). A phenomenological momentum-dependent two-body correction simulates contributions from missing three-body forces and three-body correlations. We discuss 12C with a special emphasis on the structure of the excited 0+ and 2+ states. We analyze the degree of α-clustering and confirm, taking inelastic electron scattering data into account, the conjecture that the Hoyle state has to be understood as a loosely bound system of alpha particles. We will also present first results on the application of FMD for the calculation of scattering phase shifts in 3He — 4He.
Stresses and elastic constants of crystalline sodium, from molecular dynamics
Schiferl, S.K.
1985-02-01
The stresses and the elastic constants of bcc sodium are calculated by molecular dynamics (MD) for temperatures to T = 340K. The total adiabatic potential of a system of sodium atoms is represented by pseudopotential model. The resulting expression has two terms: a large, strictly volume-dependent potential, plus a sum over ion pairs of a small, volume-dependent two-body potential. The stresses and the elastic constants are given as strain derivatives of the Helmholtz free energy. The resulting expressions involve canonical ensemble averages (and fluctuation averages) of the position and volume derivatives of the potential. An ensemble correction relates the results to MD equilibrium averages. Evaluation of the potential and its derivatives requires the calculation of integrals with infinite upper limits of integration, and integrand singularities. Methods for calculating these integrals and estimating the effects of integration errors are developed. A method is given for choosing initial conditions that relax quickly to a desired equilibrium state. Statistical methods developed earlier for MD data are extended to evaluate uncertainties in fluctuation averages, and to test for symmetry. 45 refs., 10 figs., 4 tabs.
Maintain rigid structures in Verlet based Cartesian molecular dynamics simulations
NASA Astrophysics Data System (ADS)
Tao, Peng; Wu, Xiongwu; Brooks, Bernard R.
2012-10-01
An algorithm is presented to maintain rigid structures in Verlet based Cartesian molecular dynamics (MD) simulations. After each unconstrained MD step, the coordinates of selected particles are corrected to maintain rigid structures through an iterative procedure of rotation matrix computation. This algorithm, named as SHAPE and implemented in CHARMM program suite, avoids the calculations of Lagrange multipliers, so that the complexity of computation does not increase with the number of particles in a rigid structure. The implementation of this algorithm does not require significant modification of propagation integrator, and can be plugged into any Cartesian based MD integration scheme. A unique feature of the SHAPE method is that it is interchangeable with SHAKE for any object that can be constrained as a rigid structure using multiple SHAKE constraints. Unlike SHAKE, the SHAPE method can be applied to large linear (with three or more centers) and planar (with four or more centers) rigid bodies. Numerical tests with four model systems including two proteins demonstrate that the accuracy and reliability of the SHAPE method are comparable to the SHAKE method, but with much more applicability and efficiency.
Hydration structure of salt solutions from ab initio molecular dynamics
Bankura, Arindam; Carnevale, Vincenzo; Klein, Michael L.
2013-01-07
The solvation structures of Na{sup +}, K{sup +}, and Cl{sup -} ions in aqueous solution have been investigated using density functional theory (DFT) based Car-Parrinello (CP) molecular dynamics (MD) simulations. CPMD trajectories were collected for systems containing three NaCl or KCl ion pairs solvated by 122 water molecules using three different but commonly employed density functionals (BLYP, HCTH, and PBE) with electron correlation treated at the level of the generalized gradient approximation (GGA). The effect of including dispersion forces was analyzed through the use of an empirical correction to the DFT-GGA scheme. Special attention was paid to the hydration characteristics, especially the structural properties of the first solvation shell of the ions, which was investigated through ion-water radial distribution functions, coordination numbers, and angular distribution functions. There are significant differences between the present results obtained from CPMD simulations and those provided by classical MD based on either the CHARMM force field or a polarizable model. Overall, the computed structural properties are in fair agreement with the available experimental results. In particular, the observed coordination numbers 5.0-5.5, 6.0-6.4, and 6.0-6.5 for Na{sup +}, K{sup +}, and Cl{sup -}, respectively, are consistent with X-ray and neutron scattering studies but differ somewhat from some of the many other recent computational studies of these important systems. Possible reasons for the differences are discussed.
Hydration structure of salt solutions from ab initio molecular dynamics
NASA Astrophysics Data System (ADS)
Bankura, Arindam; Carnevale, Vincenzo; Klein, Michael L.
2013-01-01
The solvation structures of Na^+, K^+, and Cl^- ions in aqueous solution have been investigated using density functional theory (DFT) based Car-Parrinello (CP) molecular dynamics (MD) simulations. CPMD trajectories were collected for systems containing three NaCl or KCl ion pairs solvated by 122 water molecules using three different but commonly employed density functionals (BLYP, HCTH, and PBE) with electron correlation treated at the level of the generalized gradient approximation (GGA). The effect of including dispersion forces was analyzed through the use of an empirical correction to the DFT-GGA scheme. Special attention was paid to the hydration characteristics, especially the structural properties of the first solvation shell of the ions, which was investigated through ion-water radial distribution functions, coordination numbers, and angular distribution functions. There are significant differences between the present results obtained from CPMD simulations and those provided by classical MD based on either the CHARMM force field or a polarizable model. Overall, the computed structural properties are in fair agreement with the available experimental results. In particular, the observed coordination numbers 5.0-5.5, 6.0-6.4, and 6.0-6.5 for Na^+, K^+, and Cl^-, respectively, are consistent with X-ray and neutron scattering studies but differ somewhat from some of the many other recent computational studies of these important systems. Possible reasons for the differences are discussed.
Johnston, Jennifer M.
2014-01-01
The majority of biological processes mediated by G Protein-Coupled Receptors (GPCRs) take place on timescales that are not conveniently accessible to standard molecular dynamics (MD) approaches, notwithstanding the current availability of specialized parallel computer architectures, and efficient simulation algorithms. Enhanced MD-based methods have started to assume an important role in the study of the rugged energy landscape of GPCRs by providing mechanistic details of complex receptor processes such as ligand recognition, activation, and oligomerization. We provide here an overview of these methods in their most recent application to the field. PMID:24158803
Molecular Dynamics Study of Potassium Azide (KN_3)
NASA Astrophysics Data System (ADS)
Ossowski, M.; Hardy, J. R.
1998-03-01
An ab initio model developed for intermolecular and intramolecular potentials in ionic molecular solids(H. M. Lu and J. R. Hardy, Phys. Rev. B, 42, 8339 (1990)) is employed to study the phase diagram of potassium azide (KN_3). We performed first-principles static structural relaxation, supercell molecular dynamics, lattice- dynamical studies and predict the existence of a high temperature rotationally disordered phase in KN_3. A selected work on other members of the alkali azide family will also be discussed.
HTMD: High-Throughput Molecular Dynamics for Molecular Discovery.
Doerr, S; Harvey, M J; Noé, Frank; De Fabritiis, G
2016-04-12
Recent advances in molecular simulations have allowed scientists to investigate slower biological processes than ever before. Together with these advances came an explosion of data that has transformed a traditionally computing-bound into a data-bound problem. Here, we present HTMD, a programmable, extensible platform written in Python that aims to solve the data generation and analysis problem as well as increase reproducibility by providing a complete workspace for simulation-based discovery. So far, HTMD includes system building for CHARMM and AMBER force fields, projection methods, clustering, molecular simulation production, adaptive sampling, an Amazon cloud interface, Markov state models, and visualization. As a result, a single, short HTMD script can lead from a PDB structure to useful quantities such as relaxation time scales, equilibrium populations, metastable conformations, and kinetic rates. In this paper, we focus on the adaptive sampling and Markov state modeling features. PMID:26949976
Dynamical analysis of highly excited molecular spectra
Kellman, M.E.
1993-12-01
The goal of this program is new methods for analysis of spectra and dynamics of highly excited vibrational states of molecules. In these systems, strong mode coupling and anharmonicity give rise to complicated classical dynamics, and make the simple normal modes analysis unsatisfactory. New methods of spectral analysis, pattern recognition, and assignment are sought using techniques of nonlinear dynamics including bifurcation theory, phase space classification, and quantization of phase space structures. The emphasis is chaotic systems and systems with many degrees of freedom.
Development of the four-body corrected fragment molecular orbital (FMO4) method
NASA Astrophysics Data System (ADS)
Nakano, Tatsuya; Mochizuki, Yuji; Yamashita, Katsumi; Watanabe, Chiduru; Fukuzawa, Kaori; Segawa, Katsunori; Okiyama, Yoshio; Tsukamoto, Takayuki; Tanaka, Shigenori
2012-01-01
The four-body corrected fragment molecular orbital (FMO4) method was implemented at the second-order Møller-Plesset perturbation (MP2) level. A series of accuracy tests relative to the previous two-body and three-body treatments were performed. As expected, FMO4 provided better accuracy in total energies in comparison with the reference values by regular MO calculations. A nonconventional fragmentation by separating main and side chains in amino acid residues was examined for Ala-pentamer and Chignolin, where the four-body corrections were shown to be substantial. A large complex of HIV-1 protease (total 198 residues) with lopinavir was calculated as well. Furthermore, this new FMO scheme was successfully applied to adamantane-shaped clusters with three-dimensional bonding framework.
An Analog Gamma Correction Scheme for High Dynamic Range CMOS Logarithmic Image Sensors
Cao, Yuan; Pan, Xiaofang; Zhao, Xiaojin; Wu, Huisi
2014-01-01
In this paper, a novel analog gamma correction scheme with a logarithmic image sensor dedicated to minimize the quantization noise of the high dynamic applications is presented. The proposed implementation exploits a non-linear voltage-controlled-oscillator (VCO) based analog-to-digital converter (ADC) to perform the gamma correction during the analog-to-digital conversion. As a result, the quantization noise does not increase while the same high dynamic range of logarithmic image sensor is preserved. Moreover, by combining the gamma correction with the analog-to-digital conversion, the silicon area and overall power consumption can be greatly reduced. The proposed gamma correction scheme is validated by the reported simulation results and the experimental results measured for our designed test structure, which is fabricated with 0.35 μm standard complementary-metal-oxide-semiconductor (CMOS) process. PMID:25517692
An analog gamma correction scheme for high dynamic range CMOS logarithmic image sensors.
Cao, Yuan; Pan, Xiaofang; Zhao, Xiaojin; Wu, Huisi
2014-01-01
In this paper, a novel analog gamma correction scheme with a logarithmic image sensor dedicated to minimize the quantization noise of the high dynamic applications is presented. The proposed implementation exploits a non-linear voltage-controlled-oscillator (VCO) based analog-to-digital converter (ADC) to perform the gamma correction during the analog-to-digital conversion. As a result, the quantization noise does not increase while the same high dynamic range of logarithmic image sensor is preserved. Moreover, by combining the gamma correction with the analog-to-digital conversion, the silicon area and overall power consumption can be greatly reduced. The proposed gamma correction scheme is validated by the reported simulation results and the experimental results measured for our designed test structure, which is fabricated with 0.35 μm standard complementary-metal-oxide-semiconductor (CMOS) process. PMID:25517692
Elucidation of molecular dynamics of invasive species of rice
Technology Transfer Automated Retrieval System (TEKTRAN)
Cultivated rice fields are aggressively invaded by weedy rice in the U.S. and worldwide. Weedy rice results in loss of yield and seed contamination. The molecular dynamics of the evolutionary adaptive traits of weedy rice are not fully understood. To understand the molecular basis and identify the i...
Molecular dynamics simulation of interfacial adhesion
Yarovsky, I.; Chaffee, A.L.
1996-12-31
Chromium salts are often used in the pretreatment stages of steel painting processes in order to improve adhesion at the metal oxide/primer interface. Although well established empirically, the chemical basis for the improved adhesion conferred by chromia is not well understood. A molecular level understanding of this behaviour should provide a foundation for the design of materials offering improved adhesion control. Molecular modelling of adhesion involves simulation and analysis of molecular behaviour at the interface between two interacting phases. The present study concerns behaviour at the boundary between the metal coated steel surface (with or without chromium pretreatment) and an organic primer based on a solid epoxide resin produced from bisphenol A and epichlorohydrin. An epoxy resin oligomer of molecular weight 3750 was used as the model for the primer.
Visualizing Functional Motions of Membrane Transporters with Molecular Dynamics Simulations
2013-01-01
Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins. PMID:23298176
Tunable Interfacial Thermal Conductance by Molecular Dynamics
NASA Astrophysics Data System (ADS)
Shen, Meng
We study the mechanism of tunable heat transfer through interfaces between solids using a combination of non-equilibrium molecular dynamics simulation (NEMD), vibrational mode analysis and wave packet simulation. We investigate how heat transfer through interfaces is affected by factors including pressure, interfacial modulus, contact area and interfacial layer thickness, with an overreaching goal of developing fundamental knowledge that will allow one to tailor thermal properties of interfacial materials. The role of pressure and interfacial stiffness is unraveled by our studies on an epitaxial interface between two Lennard-Jones (LJ) crystals. The interfacial stiffness is varied by two different methods: (i) indirectly by applying pressure which due to anharmonic nature of bonding, increases interfacial stiffness, and (ii) directly by changing the interfacial bonding strength by varying the depth of the potential well of the LJ potential. When the interfacial bonding strength is low, quantitatively similar behavior to pressure tuning is observed when the interfacial thermal conductance is increased by directly varying the potential-well depth parameter of the LJ potential. By contrast, when the interfacial bonding strength is high, thermal conductance is almost pressure independent, and even slightly decreases with increasing pressure. This decrease can be explained by the change in overlap between the vibrational densities of states of the two crystalline materials. The role of contact area is studied by modeling structures comprised of Van der Waals junctions between single-walled nanotubes (SWCNT). Interfacial thermal conductance between SWCNTs is obtained from NEMD simulation as a function of crossing angle. In this case the junction conductance per unit area is essentially a constant. By contrast, interfacial thermal conductance between multiwalled carbon nanotubes (MWCNTs) is shown to increase with diameter of the nanotubes by recent experimental studies [1
Tunable Interfacial Thermal Conductance by Molecular Dynamics
NASA Astrophysics Data System (ADS)
Shen, Meng
We study the mechanism of tunable heat transfer through interfaces between solids using a combination of non-equilibrium molecular dynamics simulation (NEMD), vibrational mode analysis and wave packet simulation. We investigate how heat transfer through interfaces is affected by factors including pressure, interfacial modulus, contact area and interfacial layer thickness, with an overreaching goal of developing fundamental knowledge that will allow one to tailor thermal properties of interfacial materials. The role of pressure and interfacial stiffness is unraveled by our studies on an epitaxial interface between two Lennard-Jones (LJ) crystals. The interfacial stiffness is varied by two different methods: (i) indirectly by applying pressure which due to anharmonic nature of bonding, increases interfacial stiffness, and (ii) directly by changing the interfacial bonding strength by varying the depth of the potential well of the LJ potential. When the interfacial bonding strength is low, quantitatively similar behavior to pressure tuning is observed when the interfacial thermal conductance is increased by directly varying the potential-well depth parameter of the LJ potential. By contrast, when the interfacial bonding strength is high, thermal conductance is almost pressure independent, and even slightly decreases with increasing pressure. This decrease can be explained by the change in overlap between the vibrational densities of states of the two crystalline materials. The role of contact area is studied by modeling structures comprised of Van der Waals junctions between single-walled nanotubes (SWCNT). Interfacial thermal conductance between SWCNTs is obtained from NEMD simulation as a function of crossing angle. In this case the junction conductance per unit area is essentially a constant. By contrast, interfacial thermal conductance between multiwalled carbon nanotubes (MWCNTs) is shown to increase with diameter of the nanotubes by recent experimental studies [1
The Computer Simulation of Liquids by Molecular Dynamics.
ERIC Educational Resources Information Center
Smith, W.
1987-01-01
Proposes a mathematical computer model for the behavior of liquids using the classical dynamic principles of Sir Isaac Newton and the molecular dynamics method invented by other scientists. Concludes that other applications will be successful using supercomputers to go beyond simple Newtonian physics. (CW)