Mass transport properties of Pu/DT mixtures from orbital free molecular dynamics simulations
Kress, Joel David; Ticknor, Christopher; Collins, Lee A.
2015-09-16
Mass transport properties (shear viscosity and diffusion coefficients) for Pu/DT mixtures were calculated with Orbital Free Molecular Dynamics (OFMD). The results were fitted to simple functions of mass density (for ρ=10.4 to 62.4 g/cm^{3}) and temperature (for T=100 up to 3,000 eV) for Pu/DT mixtures consisting of 100/0, 25/75, 50/50, and 75/25 by number.
Transport properties of lithium hydride at extreme conditions from orbital-free molecular dynamics
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Ticknor, C.; Kress, J. D.; Collins, L. A.; Lambert, F.
2013-02-01
We have performed a systematic study of lithium hydride (LiH), using orbital-free molecular dynamics, with a focus on mass transport properties such as diffusion and viscosity by extending our previous studies at the lower end of the warm, dense matter regime to cover a span of densities from ambient to 10-fold compressed and temperatures from 10 eV to 10 keV. We determine analytic formulas for self- and mutual-diffusion coefficients, and viscosity, which are in excellent agreement with our molecular dynamics results, and interpolate smoothly between liquid and dense plasma regimes. In addition, we find the orbital-free calculations begin to agree with the Brinzinskii-Landau formula above about 250 eV at which point the medium becomes fully ionized. A binary-ion model based on a bare Coulomb interaction within a neutralizing background with the effective charges determined from a regularization prescription shows good agreement above about 100 eV with the orbital-free results. Finally, we demonstrate the validity of a pressure-based mixing rule in determining the transport properties from the pure-species quantities.
Liquid Be, Ca and Ba. An orbital-free ab-initio molecular dynamics study
Rio, B. G. del; González, L. E.
2015-08-17
Several static and dynamic properties of liquid beryllium (l-Be), liquid calcium (l-Ca) and liquid barium (l-Ba) near their triple point have been evaluated by the orbital-free ab initio molecular dynamics method (OF-AIMD), where the interaction between valence electrons and ions is described by means of local pseudopotentials. These local pseudopotentials used were constructed through a force-matching process with those obtained from a Kohn-Sham ab initio molecular dynamics study (KS-AIMD) of a reduced system with non-local pseudopotentials. The calculated static structures show good agreement with the available experimental data, including an asymmetric second peak in the structure factor which has been linked to the existence of a marked icosahedral short-range order in the liquid. As for the dynamic properties, we obtain collective density excitations whose associated dispersion relations exhibit a positive dispersion.
NASA Astrophysics Data System (ADS)
Chakraborty, Debajit; Karasiev, Valentin; Trickey, Samuel
Aluminum is frequently used in warm-dense matter (WDM) experiments. However, experimental diagnostic limitations make computational exploration of the Al liquid-vapor transition important. The elevated temperaure and low-density make ab initio molecular dynamics (AIMD) with Kohn-Sham (KS) density functional theory (DFT) searches for the divergent compressibility extremely time consuming. Orbital free DFT (OFDFT) in principle is a cost-effective alternative. Here we report on calculations utilizing the PROFESS@QuantumEspresso interface to explore suitable pseudo-potentials, the limitations of our wholly constraint-based VT84F non-ineracting free-energy functional as exposed in the low-density regime, and possible extensions or extrapolations via tunable non-interacting free energy functionals. Work supported by U.S. Dept. of Energy, Grant DE-SC0002139.
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.
NASA Astrophysics Data System (ADS)
Chen, Mohan; Xia, Junchao; Huang, Chen; Dieterich, Johannes M.; Hung, Linda; Shin, Ilgyou; Carter, Emily A.
2015-05-01
Orbital-free density functional theory (OFDFT) is a linear-scaling first-principles quantum mechanics method used to calculate the ground-state energy of a given system. Here we present a new version of PRinceton Orbital-Free Electronic Structure Software (PROFESS) with new features. First, PROFESS 3.0 provides a set of new kinetic energy density functionals (KEDFs) which are designed to model semiconductors or transition metals. Specifically, PROFESS 3.0 includes the Huang-Carter (HC) KEDF [1], a density decomposition method with fixed localized electronic density [2], the Wang-Govind-Carter (WGC) decomposition KEDF [3], and the Enhanced von Weizsäcker (EvW)-WGC KEDF [4]. Other major new functions are included, such as molecular dynamics with different statistical mechanical ensembles and spin-polarized density optimizers.
Danel, J-F.; Kazandjian, L.; Zerah, G.
2009-06-15
A form of the linear mixing rule involving the equality of excess pressures is tested with various mole fractions and various types of orbital-free molecular dynamics simulations. For all the cases considered, this mixing rule yields, within statistical error, the pressure of a mixture of helium and iron obtained by a direct simulation. In an attempt to interpret the robustness of the mixing rule, we show that it can be derived from thermodynamic stability if the system is regarded as a mixture of independent effective average atoms. The success of the mixing rule applied with equations of state including various degrees of approximation leads us to suggest its use in the thermodynamic domain where quantum molecular dynamics can be implemented.
NASA Astrophysics Data System (ADS)
Danel, J.-F.; Kazandjian, L.
2015-01-01
We test two isothermal-isobaric mixing rules, respectively based on excess-pressure and total-pressure equilibration, applied to the equation of state of a dense plasma. While the equation of state is generally known for pure species, that of arbitrary mixtures is not available so that the validation of accurate mixing rules, that implies resorting to first-principles simulations, is very useful. Here we consider the case of a plastic with composition C2H3 and we implement two complementary ab initio approaches adapted to the dense plasma domain: quantum molecular dynamics, limited to low temperature by its computational cost, and orbital-free molecular dynamics, that can be implemented at high temperature. The temperature and density range considered is 1-10 eV and 0.6-10 g/cm 3 for quantum molecular dynamics, and 5-1000 eV and 1-10 g/cm 3 for orbital-free molecular dynamics. Simulations for the full C2H3 mixture are the benchmark against which to assess the mixing rules, and both pressure and internal energy are compared. We find that the mixing rule based on excess-pressure equilibration is overall more accurate than that based on total-pressure equilibration; except for quantum molecular dynamics and a thermodynamic domain characterized by very low or negative excess pressures, it gives pressures which are generally within statistical error or within 1% of the exact ones. Besides, its superiority is amplified in the calculation of a principal Hugoniot.
Burakovsky, Leonid; Kress, Joel D.; Collins, Lee A.
2012-05-31
Mass transport properties for LiD-U mixtures were calculated using a pressure matching mixture rule for the mixing of LiD and of U properties simulated with Orbital Free Molecular Dynamics (OFMD). The mixing rule was checked against benchmark OFMD simulations for the fully interacting three-component (Li, D, U) system. To obtain transport coefficients for LiD-U mixtures of different (LiD){sub x}U{sub (1-x)} compositions as functions of temperature and mixture density is a tedious task. Quantum molecular dynamics (MD) simulations can be employed, as in the case LiD or U. However, due to the presence of the heavy constituent U, such simulations proceed so slowly that only a limited number of numerical data points in the (x, {rho}, T) phase space can be obtained. To finesse this difficulty, transport coefficients for a mixture can be obtained using a pressure-matching mixing rule discussed. For both LiD and U, the corresponding transport coefficients were obtained earlier from quantum molecular dynamics simulations. In these simulations, the quantum behavior of the electrons was represented using an orbital free (OF) version of density functional theory, and ions were advanced in time using classical molecular dynamics. The total pressure of the system, P = nk{sub B}T/V + P{sub e}, is the sum of the ideal gas pressure of the ions plus the electron pressure. The mass self-diffusion coefficient for species {alpha}, D{sub {alpha}}, the mutual diffusion coefficient for species {alpha} and {beta}, D{alpha}{beta}, and the shear viscosity, {eta}, are computed from the appropriate autocorrelation function. The details of similar QMD calculations on LiH are described in Ref. [1] for 0.5 eV < T < 3 eV, and in Ref. [2] for 2 eV < T < 6 eV.
NASA Astrophysics Data System (ADS)
Romere, P. O.
1982-03-01
A proposed configuration for a Space Operations Center is presented in its eight stages of buildup. The on orbit aerodynamic force and moment characteristics were calculated for each stage based upon free molecular flow theory. Calculation of the aerodynamic characteristics was accomplished through the use of an orbital aerodynamic computer program, and the computation method is described with respect to the free molecular theory used. The aerodynamic characteristics are presented in tabulated form for each buildup stage at angles of attack from 0 to 360 degrees and roll angles from -60 to +60 degrees. The reference altitude is 490 kilometers, however, the data should be applicable for altitudes below 490 kilometers down to approximately 185 kilometers.
NASA Technical Reports Server (NTRS)
Romere, P. O.
1982-01-01
A proposed configuration for a Space Operations Center is presented in its eight stages of buildup. The on orbit aerodynamic force and moment characteristics were calculated for each stage based upon free molecular flow theory. Calculation of the aerodynamic characteristics was accomplished through the use of an orbital aerodynamic computer program, and the computation method is described with respect to the free molecular theory used. The aerodynamic characteristics are presented in tabulated form for each buildup stage at angles of attack from 0 to 360 degrees and roll angles from -60 to +60 degrees. The reference altitude is 490 kilometers, however, the data should be applicable for altitudes below 490 kilometers down to approximately 185 kilometers.
Rio, Beatriz G del; González, Luis E
2014-11-19
We have performed a comprehensive study of the properties of liquid Be, Ca and Ba, through the use of orbital free ab initio simulations. To this end we have developed a force-matching method to construct the necessary local pseudopotentials from standard ab initio calculations. The structural magnitudes are analyzed, including the average and local structures and the dynamic properties are studied. We find several common features, like an asymmetric second peak in the structure factor, a large amount of local structures with five-fold symmetry, a quasi-universal behaviour of the single-particle dynamic properties and a large degree of positive dispersion in the propagation of collective density fluctuations, whose damping is dictated by slow thermal relaxations and fast viscoelastic ones. Some peculiarities in the dynamic properties are however observed, like a very high sound velocity and a large violation of the Stokes-Einstein relation for Be, or an extremely high positive dispersion and a large slope in the dispersion relation of shear waves at the onset of the wavevector region where they are supported for Ba. PMID:25347355
NASA Astrophysics Data System (ADS)
Xu, Fuming; Wang, Bin; Wei, Yadong; Wang, Jian
2013-10-01
Orbital-free density functional theory (OFDFT) replaces the wavefunction in the kinetic energy by an explicit energy functional and thereby speeds up significantly the calculation of ground state properties of the solid state systems. So far, the application of OFDFT has been centered on closed systems and less attention is paid on the transport properties in open systems. In this paper, we use OFDFT and combine it with non-equilibrium Green's function to simulate equilibrium electronic transport properties in silicon nanostructures from first principles. In particular, we study ac transport properties of a silicon atomic junction consisting of a silicon atomic chain and two monoatomic leads. We have calculated the dynamic conductance of this atomic junction as a function of ac frequency with one to four silicon atoms in the central scattering region. Although the system is transmissive with dc conductance around 4 to 5 e2/h, capacitive-like behavior was found in the finite frequency regime. Our analysis shows that, up to 0.1 THz, this behavior can be characterized by a classic RC circuit consisting of two resistors and a capacitor. One resistor gives rise to dc resistance and the other one accounts for the charge relaxation resistance with magnitude around 0.2 h/e2 when the silicon chain contains two atoms. It was found that the capacitance is around 5 aF for the same system.
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.
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.
Structural and dynamical properties of hot dense matter by a Thomas-Fermi-Dirac molecular dynamics
NASA Astrophysics Data System (ADS)
Lambert, F.; Clérouin, J.; Mazevet, S.
2006-09-01
We use a model combining, in a consistent way, orbital-free density functional theory (OF-DFT) and molecular dynamics (MD), to compute the thermodynamical, structural and dynamical properties of Fe and Au plasmas at conditions relevant to astrophysics and inertial confinement fusion (ICF). The newly developed parallel numerical scheme presented here allows to propagate hundreds of particles and to obtain accurate transport properties. This allows us to investigate the validity of the commonly used one-component plasma (OCP) model in predicting the pair correlation, the diffusion and viscosity coefficients for these two high-temperature high-density plasmas.
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
NASA Astrophysics Data System (ADS)
Danel, J.-F.; Kazandjian, L.; Piron, R.
2016-04-01
Of the two approaches of density-functional theory molecular dynamics, quantum molecular dynamics is limited at high temperature by computational cost whereas orbital-free molecular dynamics, based on an approximation of the kinetic electronic free energy, can be implemented in this domain. In the case of deuterium, it is shown how orbital-free molecular dynamics can be regarded as the limit of quantum molecular dynamics at high temperature for the calculation of the equation of state. To this end, accurate quantum molecular dynamics calculations are performed up to 20 eV at mass densities as low as 0.5 g /cm3 and up to 10 eV at mass densities as low as 0.2 g /cm3 . As a result, the limitation in temperature so far attributed to quantum molecular dynamics is overcome and an approach combining quantum and orbital-free molecular dynamics is used to construct an equation of state of deuterium. The thermodynamic domain addressed is that of the fluid phase above 1 eV and 0.2 g /cm3 . Both pressure and internal energy are calculated as functions of temperature and mass density, and various exchange-correlation contributions are compared. The generalized gradient approximation of the exchange-correlation functional, corrected to approximately include the influence of temperature, is retained and the results obtained are compared to other approaches and to experimental shock data; in parts of the thermodynamic domain addressed, these results significantly differ from those obtained in other first-principles investigations which themselves disagree. The equations of state of hydrogen and tritium above 1 eV and above, respectively, 0.1 g /cm3 and 0.3 g /cm3 , can be simply obtained by mass density scaling from the results found for deuterium. This ab initio approach allows one to consistently cover a very large domain of temperature on the domain of mass density outlined above.
Sjostrom, Travis; Crockett, Scott
2015-09-02
The liquid regime equation of state of silicon dioxide SiO_{2} is calculated via quantum molecular dynamics in the density range of 5 to 15 g/cc and with temperatures from 0.5 to 100 eV, including the α-quartz and stishovite phase Hugoniot curves. Below 8 eV calculations are based on Kohn-Sham density functional theory (DFT), and above 8 eV a new orbital-free DFT formulation, presented here, based on matching Kohn-Sham DFT calculations is employed. Recent experimental shock data are found to be in very good agreement with the current results. Finally both experimental and simulation data are used in constructing a new liquid regime equation of state table for SiO_{2}.
Sjostrom, Travis; Crockett, Scott
2015-09-02
The liquid regime equation of state of silicon dioxide SiO2 is calculated via quantum molecular dynamics in the density range of 5 to 15 g/cc and with temperatures from 0.5 to 100 eV, including the α-quartz and stishovite phase Hugoniot curves. Below 8 eV calculations are based on Kohn-Sham density functional theory (DFT), and above 8 eV a new orbital-free DFT formulation, presented here, based on matching Kohn-Sham DFT calculations is employed. Recent experimental shock data are found to be in very good agreement with the current results. Finally both experimental and simulation data are used in constructing a newmore » liquid regime equation of state table for SiO2.« less
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.
Petascale Orbital-Free Density Functional Theory Enabled by Small-Box Algorithms.
Chen, Mohan; Jiang, Xiang-Wei; Zhuang, Houlong; Wang, Lin-Wang; Carter, Emily A
2016-06-14
Orbital-free density functional theory (OFDFT) is a quantum-mechanics-based method that utilizes electron density as its sole variable. The main computational cost in OFDFT is the ubiquitous use of the fast Fourier transform (FFT), which is mainly adopted to evaluate the kinetic energy density functional (KEDF) and electron-electron Coulomb interaction terms. We design and implement a small-box FFT (SBFFT) algorithm to overcome the parallelization limitations of conventional FFT algorithms. We also propose real-space truncation of the nonlocal Wang-Teter KEDF kernel. The scalability of the SBFFT is demonstrated by efficiently simulating one full optimization step (electron density, energies, forces, and stresses) of 1,024,000 lithium (Li) atoms on up to 65,536 cores. We perform other tests using Li as a test material, including calculations of physical properties of different phases of bulk Li, geometry optimizations of nanocrystalline Li, and molecular dynamics simulations of liquid Li. All of the tests yield excellent agreement with the original OFDFT results, suggesting that the OFDFT-SBFFT algorithm opens the door to efficient first-principles simulations of materials containing millions of atoms. PMID:27145175
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.
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
Sheppard, Daniel; Kress, Joel D; Crockett, Scott; Collins, Lee A; Desjarlais, Michael P
2014-12-01
The shock Hugoniot for lithium 6 deuteride ((6)LiD) was calculated via first principles using Kohn-Sham density-functional theory molecular dynamics (KSMD) for temperatures of 0.5-25 eV. The upper limit of 25 eV represents a practical limit where KSMD is no longer computationally feasible due to the number of electronic bands which are required to be populated. To push the Hugoniot calculations to higher temperatures we make use of orbital-free density-functional theory molecular dynamics (OFMD). Thomas-Fermi-Dirac-based OFMD gives a poor description of the electronic structure at low temperatures so the initial state is not well defined. We propose a method of bootstrapping the Hugoniot from OFMD to the Hugoniot from KSMD between 10 and 20 eV, where the two methods are in agreement. The combination of KSMD and OFMD allows construction of a first-principles Hugoniot from the initial state to 1000 eV. Theoretical shock-compression results are in good agreement with available experimental data and exhibit the appropriate high-temperature limits. We show that a unified KSMD-OFMD Hugoniot can be used to assess the quality of the existing equation-of-state (EOS) models and inform better EOS models based on justifiable physics. PMID:25615229
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
Chakraborty, Debdutta; Kar, Susmita; Chattaraj, Pratim Kumar
2015-12-21
The orbital free density functional theory and the single density equation approach are formally equivalent. An orbital free density based quantum dynamical strategy is used to study the quantum-classical correspondence in both weakly and strongly coupled van der Pol and Duffing oscillators in the presence of an external electric field in one dimension. The resulting quantum hydrodynamic equations of motion are solved through an implicit Euler type real space method involving a moving weighted least square technique. The Lagrangian framework used here allows the numerical grid points to follow the wave packet trajectory. The associated classical equations of motion are solved using a sixth order Runge-Kutta method and the Ehrenfest dynamics is followed through the solution of the time dependent Schrodinger equation using a time dependent Fourier Grid Hamiltonian technique. Various diagnostics reveal a close parallelism between classical regular as well as chaotic dynamics and that obtained from the Bohmian mechanics. PMID:26033095
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)
Danel, J-F; Kazandjian, L; Piron, R
2016-04-01
Of the two approaches of density-functional theory molecular dynamics, quantum molecular dynamics is limited at high temperature by computational cost whereas orbital-free molecular dynamics, based on an approximation of the kinetic electronic free energy, can be implemented in this domain. In the case of deuterium, it is shown how orbital-free molecular dynamics can be regarded as the limit of quantum molecular dynamics at high temperature for the calculation of the equation of state. To this end, accurate quantum molecular dynamics calculations are performed up to 20 eV at mass densities as low as 0.5g/cm^{3} and up to 10 eV at mass densities as low as 0.2g/cm^{3}. As a result, the limitation in temperature so far attributed to quantum molecular dynamics is overcome and an approach combining quantum and orbital-free molecular dynamics is used to construct an equation of state of deuterium. The thermodynamic domain addressed is that of the fluid phase above 1 eV and 0.2g/cm^{3}. Both pressure and internal energy are calculated as functions of temperature and mass density, and various exchange-correlation contributions are compared. The generalized gradient approximation of the exchange-correlation functional, corrected to approximately include the influence of temperature, is retained and the results obtained are compared to other approaches and to experimental shock data; in parts of the thermodynamic domain addressed, these results significantly differ from those obtained in other first-principles investigations which themselves disagree. The equations of state of hydrogen and tritium above 1 eV and above, respectively, 0.1g/cm^{3} and 0.3g/cm^{3}, can be simply obtained by mass density scaling from the results found for deuterium. This ab initio approach allows one to consistently cover a very large domain of temperature on the domain of mass density outlined above. PMID:27176421
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%.
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.
Equation of state of dense plasmas with pseudoatom molecular dynamics
NASA Astrophysics Data System (ADS)
Starrett, C. E.; Saumon, D.
2016-06-01
We present an approximation for calculating the equation of state (EOS) of warm and hot dense matter that is built on the previously published pseudoatom molecular dynamics (PAMD) model of dense plasmas [Starrett et al., Phys. Rev. E 91, 013104 (2015), 10.1103/PhysRevE.91.013104]. While the EOS calculation with PAMD was previously limited to orbital-free density functional theory (DFT), the new approximation presented here allows a Kohn-Sham DFT treatment of the electrons. The resulting EOS thus includes a quantum mechanical treatment of the electrons with a self-consistent model of the ionic structure, while remaining tractable at high temperatures. The method is validated by comparisons with pressures from ab initio simulations of Be, Al, Si, and Fe. The EOS in the Thomas-Fermi approximation shows remarkable thermodynamic consistency over a wide range of temperatures for aluminum. We calculate the principal Hugoniots of aluminum and silicon up to 500 eV. We find that the ionic structure of the plasma has a modest effect that peaks at temperatures of a few eV and that the features arising from the electronic structure agree well with ab initio simulations.
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.
Augmented Lagrangian formulation of orbital-free density functional theory
Suryanarayana, Phanish Phanish, Deepa
2014-10-15
We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic Orbital-Free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a sequence of minimization problems without any constraint, thereby making it amenable to powerful unconstrained optimization algorithms. Further, we develop a parallel implementation of this approach for the Thomas–Fermi–von Weizsacker (TFW) kinetic energy functional in the framework of higher-order finite-differences and the conjugate gradient method. With this implementation, we establish that the Augmented Lagrangian approach is highly competitive compared to the penalty and Lagrange multiplier methods. Additionally, we show that higher-order finite-differences represent a computationally efficient discretization for performing OF-DFT simulations. Overall, we demonstrate that the proposed formulation and implementation are both efficient and robust by studying selected examples, including systems consisting of thousands of atoms. We validate the accuracy of the computed energies and forces by comparing them with those obtained by existing plane-wave methods.
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.
Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lower temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.
Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
Sjostrom, Travis; Daligault, Jerome
2014-10-10
Here, we develop and implement a new quantum molecular dynamics approximation that allows fast and accurate simulations of dense plasmas from cold to hot conditions. The method is based on a carefully designed orbital-free implementation of density functional theory. The results for hydrogen and aluminum are in very good agreement with Kohn-Sham (orbital-based) density functional theory and path integral Monte Carlo calculations for microscopic features such as the electron density as well as the equation of state. The present approach does not scale with temperature and hence extends to higher temperatures than is accessible in the Kohn-Sham method and lowermore » temperatures than is accessible by path integral Monte Carlo calculations, while being significantly less computationally expensive than either of those two methods.« less
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
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
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
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.
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.
Nonadiabatic Molecular Dynamics with Trajectories
NASA Astrophysics Data System (ADS)
Tavernelli, Ivano
2012-02-01
In the mixed quantum-classical description of molecular systems, only the quantum character of the electronic degrees of freedom is considered while the nuclear motion is treated at a classical level. In the adiabatic case, this picture corresponds to the Born-Oppenheimer limit where the nuclei move as point charges on the potential energy surface (PES) associated with a given electronic state. Despite the success of this approximation, many physical and chemical processes do not fall in the regime where nuclei and electrons can be considered decoupled. In particular, most photoreactions pass through regions of the PES in which electron-nuclear quantum interference effects are sizeable and often crucial for a correct description of the phenomena. Recently, we have developed a trajectory-based nonadiabatic molecular dynamics scheme that describes the nuclear wavepacket as an ensemble of particles following classical trajectories on PESs derived from time-dependent density functional theory (TDDFT) [1]. The method is based on Tully's fewest switches trajectories surface hopping (TSH) where the nonadiabatic coupling elements between the different potential energy surfaces are computed on-the-fly as functionals of the ground state electron density or, equivalently, of the corresponding Kohn-Sham orbitals [2]. Here, we present the theoretical fundamentals of our approach together with an extension that allows for the direct coupling of the dynamics to an external electromagnetic field [3] as well as to the external potential generated by the environment (solvent effects) [4]. The method is applied to the study of the photodissociation dynamics of simple molecules in gas phase and to the description of the fast excited state dynamics of molecules in solution (in particular Ruthenium (II) tris(bipyridine) in water). [4pt] [1] E. Tapavicza, I. Tavernelli, U. Rothlisberger, Phys. Rev. Lett., 98, (2007) 023001. [0pt] [2] Tavernelli I.; Tapavicza E.; Rothlisberger U., J. Chem
Molecular dynamics simulation of benzene
NASA Astrophysics Data System (ADS)
Trumpakaj, Zygmunt; Linde, Bogumił B. J.
2016-03-01
Intermolecular potentials and a few models of intermolecular interaction in liquid benzene are tested by Molecular Dynamics (MD) simulations. The repulsive part of the Lennard-Jones 12-6 (LJ 12-6) potential is too hard, which yields incorrect results. The exp-6 potential with a too hard repulsive term is also often used. Therefore, we took an expa-6 potential with a small Gaussian correction plus electrostatic interactions. This allows to modify the curvature of the potential. The MD simulations are carried out in the temperature range 280-352 K under normal pressure and at experimental density. The Rayleigh scattering of depolarized light is used for comparison. The results of MD simulations are comparable with the experimental values.
Buckybomb: Reactive Molecular Dynamics Simulation.
Chaban, Vitaly V; Fileti, Eudes Eterno; Prezhdo, Oleg V
2015-03-01
Energetic materials, such as explosives, propellants, and pyrotechnics, are widely used in civilian and military applications. Nanoscale explosives represent a special group because of the high density of energetic covalent bonds. The reactive molecular dynamics (ReaxFF) study of nitrofullerene decomposition reported here provides a detailed chemical mechanism of explosion of a nanoscale carbon material. Upon initial heating, C60(NO2)12 disintegrates, increasing temperature and pressure by thousands of Kelvins and bars within tens of picoseconds. The explosion starts with NO2 group isomerization into C-O-N-O, followed by emission of NO molecules and formation of CO groups on the buckyball surface. NO oxidizes into NO2, and C60 falls apart, liberating CO2. At the highest temperatures, CO2 gives rise to diatomic carbon. The study shows that the initiation temperature and released energy depend strongly on the chemical composition and density of the material. PMID:26262672
Molecular dynamics of membrane proteins.
Woolf, Thomas B.; Crozier, Paul Stewart; Stevens, Mark Jackson
2004-10-01
Understanding the dynamics of the membrane protein rhodopsin will have broad implications for other membrane proteins and cellular signaling processes. Rhodopsin (Rho) is a light activated G-protein coupled receptor (GPCR). When activated by ligands, GPCRs bind and activate G-proteins residing within the cell and begin a signaling cascade that results in the cell's response to external stimuli. More than 50% of all current drugs are targeted toward G-proteins. Rho is the prototypical member of the class A GPCR superfamily. Understanding the activation of Rho and its interaction with its Gprotein can therefore lead to a wider understanding of the mechanisms of GPCR activation and G-protein activation. Understanding the dark to light transition of Rho is fully analogous to the general ligand binding and activation problem for GPCRs. This transition is dependent on the lipid environment. The effect of lipids on membrane protein activity in general has had little attention, but evidence is beginning to show a significant role for lipids in membrane protein activity. Using the LAMMPS program and simulation methods benchmarked under the IBIG program, we perform a variety of allatom molecular dynamics simulations of membrane proteins.
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
NASA Astrophysics Data System (ADS)
Xia, Junchao; Carter, Emily A.
2016-03-01
We propose a simple density decomposition formalism within orbital-free (OF) density functional theory (DFT) based on the Wang-Govind-Carter-decomposition (WGCD) kinetic energy density functional (KEDF). The resulting simple-WGCD (sWGCD) KEDF provides efficient density optimization, full cell relaxation, reasonable bulk properties for various materials compared to both the original OFDFT-WGCD and the Kohn-Sham (KS) DFT values, and has various numerical benefits including more stable convergence and lower computational cost (twice as fast as the WGCD KEDF). We also study amorphous (a-) Li-Si alloys with KSDFT and OFDFT using the Huang-Carter (HC), WGCD, and sWGCD KEDFs. The a-Li-Si alloy samples are prepared with the anneal-and-quench method using NVT molecular dynamics simulations. We report structural properties, equilibrium volumes, bulk moduli, and alloy formation energies for each a-alloy. The HC, WGCD, and sWGCD KEDFs within OFDFT all predict accurate equilibrium volumes compared against KSDFT benchmarks. The HC KEDF bulk moduli agree with KSDFT benchmarks whereas the WGCD/sWGCD KEDFs generally overestimate the bulk moduli, especially for alloys with low Li concentrations. All three KEDFs show limited ability to predict alloy formation energies, which indicates the lack of transferability of these KEDFs among such systems and motivates future developments in OFDFT and KEDF formalisms.
Quantum molecular dynamics simulations of transport properties in liquid and dense-plasma plutonium
Kress, J. D.; Cohen, James S.; Kilcrease, D. P.; Horner, D. A.; Collins, L. A.
2011-02-15
We have calculated the viscosity and self-diffusion coefficients of plutonium in the liquid phase using quantum molecular dynamics (QMD) and in the dense-plasma phase using orbital-free molecular dynamics (OFMD), as well as in the intermediate warm dense matter regime with both methods. Our liquid metal results for viscosity are about 40% lower than measured experimentally, whereas a previous calculation using an empirical interatomic potential (modified embedded-atom method) obtained results 3-4 times larger than the experiment. The QMD and OFMD results agree well at the intermediate temperatures. The calculations in the dense-plasma regime for temperatures from 50 to 5000 eV and densities about 1-5 times ambient are compared with the one-component plasma (OCP) model, using effective charges given by the average-atom code inferno. The inferno-OCP model results agree with the OFMD to within about a factor of 2, except for the viscosity at temperatures less than about 100 eV, where the disagreement is greater. A Stokes-Einstein relationship of the viscosities and diffusion coefficients is found to hold fairly well separately in both the liquid and dense-plasma regimes.
Danel, J-F; Kazandjian, L; Zérah, G
2012-06-01
Computations of the self-diffusion coefficient and viscosity in warm dense matter are presented with an emphasis on obtaining numerical convergence and a careful evaluation of the standard deviation. The transport coefficients are computed with the Green-Kubo relation and orbital-free molecular dynamics at the Thomas-Fermi-Dirac level. The numerical parameters are varied until the Green-Kubo integral is equal to a constant in the t→+∞ limit; the transport coefficients are deduced from this constant and not by extrapolation of the Green-Kubo integral. The latter method, which gives rise to an unknown error, is tested for the computation of viscosity; it appears that it should be used with caution. In the large domain of coupling constant considered, both the self-diffusion coefficient and viscosity turn out to be well approximated by simple analytical laws using a single effective atomic number calculated in the average-atom model. PMID:23005237
Non-Periodic Finite-Element Formulation of Orbital-Free Density Functional Theory
Gavini, V; Knap, J; Bhattacharya, K; Ortiz, M
2006-10-06
We propose an approach to perform orbital-free density functional theory calculations in a non-periodic setting using the finite-element method. We consider this a step towards constructing a seamless multi-scale approach for studying defects like vacancies, dislocations and cracks that require quantum mechanical resolution at the core and are sensitive to long range continuum stresses. In this paper, we describe a local real space variational formulation for orbital-free density functional theory, including the electrostatic terms and prove existence results. We prove the convergence of the finite-element approximation including numerical quadratures for our variational formulation. Finally, we demonstrate our method using examples.
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.
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.
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.
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.
NASA Astrophysics Data System (ADS)
Marqués, M.; González, D. J.; González, L. E.
2016-07-01
The melting curve of sodium for a pressure range up to 100 GPa has been evaluated by the orbital free ab initio molecular dynamics method. This method uses the electronic density as the basic variable combined with an approximate electronic kinetic energy functional and a local ionic pseudopotential and makes it possible to perform simulations with a large number of particles and for long simulation times. The calculated melting curve shows a maximum melting temperature at a pressure around 30 GPa followed by a steep decrease up to 100 GPa. For various pressures and temperatures we have evaluated several static properties, including average and local structure, electronic properties, like the electron localization function (ELF), and dynamic properties, both single-particle and collective ones, from which some transport coefficients are deduced. Despite the accurate reproduction of the available experimental data, we do not observe any indication of an early transition from a bcc-like to an fcc-like liquid, as suggested previously by other authors, but rather pressure-induced change in the variation of icosahedral-like order and bcc-like order, with no sign of fcc-like structures in the whole liquid range studied. We also consider the evolution of the ELF within this type of local arrangement upon pressurization. In the dynamic realm, we find an enlarged wave-vector region where atomic collisions play an important role in the dynamic properties of the system as pressure is increased and temperature decreased along the melting line, leading to a peculiar behavior of the dynamic properties.
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
Molecular dynamic simulations of ocular tablet dissolution.
Ru, Qian; Fadda, Hala M; Li, Chung; Paul, Daniel; Khaw, Peng T; Brocchini, Steve; Zloh, Mire
2013-11-25
Small tablets for implantation into the subconjunctival space in the eye are being developed to inhibit scarring after glaucoma filtration surgery (GFS). There is a need to evaluate drug dissolution at the molecular level to determine how the chemical structure of the active may correlate with dissolution in the nonsink conditions of the conjunctival space. We conducted molecular dynamics simulations to study the dissolution process of tablets derived from two drugs that can inhibit fibrosis after GFS, 5-fluorouracil (5-FU) and the matrix metalloprotease inhibitor (MMPi), ilomastat. The dissolution was simulated in the presence of simple point charge (SPC) water molecules, and the liquid turnover of the aqueous humor in the subconjunctival space was simulated by removal of the dissolved drug molecules at regular intervals and replacement by new water molecules. At the end of the simulation, the total molecular solvent accessible surface area of 5-FU tablets increased by 60 times more than that of ilomastat as a result of tablet swelling and release of molecules into solution. The tablet dissolution pattern shown in our molecular dynamic simulations tends to correlate with experimental release profiles. This work indicates that a series of molecular dynamic simulations can be used to predict the influence of the molecular properties of a drug on its dissolution profile and could be useful during preformulation where sufficient amounts of the drug are not always available to perform dissolution studies. PMID:24073784
Molecular Scale Dynamics of Large Ring Polymers
NASA Astrophysics Data System (ADS)
Gooßen, S.; Brás, A. R.; Krutyeva, M.; Sharp, M.; Falus, P.; Feoktystov, A.; Gasser, U.; Pyckhout-Hintzen, W.; Wischnewski, A.; Richter, D.
2014-10-01
We present neutron scattering data on the structure and dynamics of melts from polyethylene oxide rings with molecular weights up to ten times the entanglement mass of the linear counterpart. The data reveal a very compact conformation displaying a structure approaching a mass fractal, as hypothesized by recent simulation work. The dynamics is characterized by a fast Rouse relaxation of subunits (loops) and a slower dynamics displaying a lattice animal-like loop displacement. The loop size is an intrinsic property of the ring architecture and is independent of molecular weight. This is the first experimental observation of the space-time evolution of segmental motion in ring polymers illustrating the dynamic consequences of their topology that is unique among all polymeric systems of any other known architecture.
Dynamic signature of molecular association in methanol
NASA Astrophysics Data System (ADS)
Bertrand, C. E.; Self, J. L.; Copley, J. R. D.; Faraone, A.
2016-07-01
Quasielastic neutron scattering measurements and molecular dynamics simulations were combined to investigate the collective dynamics of deuterated methanol, CD3OD. In the experimentally determined dynamic structure factor, a slow, non-Fickian mode was observed in addition to the standard density-fluctuation heat mode. The simulation results indicate that the slow dynamical process originates from the hydrogen bonding of methanol molecules. The qualitative behavior of this mode is similar to the previously observed α-relaxation in supercooled water [M. C. Bellissent-Funel et al., Phys. Rev. Lett. 85, 3644 (2000)] which also originates from the formation and dissolution of hydrogen-bonded associates (supramolecular clusters). In methanol, however, this mode is distinguishable well above the freezing transition. This finding indicates that an emergent slow mode is not unique to supercooled water, but may instead be a general feature of hydrogen-bonding liquids and associating molecular liquids.
Dynamic signature of molecular association in methanol.
Bertrand, C E; Self, J L; Copley, J R D; Faraone, A
2016-07-01
Quasielastic neutron scattering measurements and molecular dynamics simulations were combined to investigate the collective dynamics of deuterated methanol, CD3OD. In the experimentally determined dynamic structure factor, a slow, non-Fickian mode was observed in addition to the standard density-fluctuation heat mode. The simulation results indicate that the slow dynamical process originates from the hydrogen bonding of methanol molecules. The qualitative behavior of this mode is similar to the previously observed α-relaxation in supercooled water [M. C. Bellissent-Funel et al., Phys. Rev. Lett. 85, 3644 (2000)] which also originates from the formation and dissolution of hydrogen-bonded associates (supramolecular clusters). In methanol, however, this mode is distinguishable well above the freezing transition. This finding indicates that an emergent slow mode is not unique to supercooled water, but may instead be a general feature of hydrogen-bonding liquids and associating molecular liquids. PMID:27394112
Molecular scale dynamics of large ring polymers.
Gooßen, S; Brás, A R; Krutyeva, M; Sharp, M; Falus, P; Feoktystov, A; Gasser, U; Pyckhout-Hintzen, W; Wischnewski, A; Richter, D
2014-10-17
We present neutron scattering data on the structure and dynamics of melts from polyethylene oxide rings with molecular weights up to ten times the entanglement mass of the linear counterpart. The data reveal a very compact conformation displaying a structure approaching a mass fractal, as hypothesized by recent simulation work. The dynamics is characterized by a fast Rouse relaxation of subunits (loops) and a slower dynamics displaying a lattice animal-like loop displacement. The loop size is an intrinsic property of the ring architecture and is independent of molecular weight. This is the first experimental observation of the space-time evolution of segmental motion in ring polymers illustrating the dynamic consequences of their topology that is unique among all polymeric systems of any other known architecture. PMID:25361284
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.
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.
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 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 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
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
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
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.
Molecular dynamics study of cyclohexane interconversion
NASA Astrophysics Data System (ADS)
Wilson, Michael A.; Chandler, David
1990-12-01
Classical molecular dynamics calculations are reported for one C 6H 12 molecule in a bath of 250 CS 2 molecules at roomtemperature and liquid densities of 1.0, 1.3, 1.4 and 1.5 g/cm 3. The solvent contribution to the free energy of activation for the chair-boat isomerization has been determined to high accuracy. The transmission coefficient and reactive flux correlation functions have also been computed. The results obtained agree with earlier conclusions drawn from RISM integral equation calculations and stochastic molecular dynamics calculations. Namely, the solvent effect on the rate manifests a qualitative breakdown of transition state theory and the RRKM picture of unimolecular kinetics. Analysis of the activated trajectories indicate a significant degree of quasiperiodicity.
Molecular dynamics 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
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
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.
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
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
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.
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
Orbital-free density functional theory implementation with the projector augmented-wave method
Lehtomäki, Jouko; Makkonen, Ilja; Harju, Ari; Lopez-Acevedo, Olga; Caro, Miguel A.
2014-12-21
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
Orbital-free density functional theory implementation with the projector augmented-wave method
NASA Astrophysics Data System (ADS)
Lehtomäki, Jouko; Makkonen, Ilja; Caro, Miguel A.; Harju, Ari; Lopez-Acevedo, Olga
2014-12-01
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
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.
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.
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
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.
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.
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.
Crystallization of nickel nanoclusters by molecular dynamics
NASA Astrophysics Data System (ADS)
Chamati, H.; Gaminchev, K.
2012-12-01
We investigated the melting properties of bulk nickel and the crystallization of nickel nanocrystals via molecular dynamics using a potential in the framework of the second moment approximation of tight-binding theory. The melting behavior was simulated with the hysteresis approach by subsequently heating and cooling gradually the system over a wide range of temperatures. The crystallization of nickel nanoclusters consisting of 55, 147 and 309 atoms was achieved after repeatedly annealing and quenching the corresponding quasicrystals several times to avoid being trapped in a local energy minimum. The time over which the global minimum was reached was found to increase with the cluster size.
Molecular dynamics 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.
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 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 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 dynamics at constant temperature and pressure
NASA Astrophysics Data System (ADS)
Toxvaerd, S.
1993-01-01
Algorithms for molecular dynamics (MD) at constant temperature and pressure are investigated. The ability to remain in a regular orbit in an intermittent chaotic regime is used as a criterion for long-time stability. A simple time-centered algorithm (leap frog) is found to be the most stable of the commonly used algorithms in MD. A model of N one-dimensional dimers with a double-well intermolecular potential, for which the distribution functions at constant temperature T and pressure P can be calculated, is used to investigate MD-NPT dynamics. A time-centered NPT algorithm is found to sample correctly and to be very robust with respect to volume scaling.
Charge transport network dynamics in molecular aggregates.
Jackson, Nicholas E; Chen, Lin X; Ratner, Mark A
2016-08-01
Due to the nonperiodic nature of charge transport in disordered systems, generating insight into static charge transport networks, as well as analyzing the network dynamics, can be challenging. Here, we apply time-dependent network analysis to scrutinize the charge transport networks of two representative molecular semiconductors: a rigid n-type molecule, perylenediimide, and a flexible p-type molecule, [Formula: see text] Simulations reveal the relevant timescale for local transfer integral decorrelation to be [Formula: see text]100 fs, which is shown to be faster than that of a crystalline morphology of the same molecule. Using a simple graph metric, global network changes are observed over timescales competitive with charge carrier lifetimes. These insights demonstrate that static charge transport networks are qualitatively inadequate, whereas average networks often overestimate network connectivity. Finally, a simple methodology for tracking dynamic charge transport properties is proposed. PMID:27439871
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.
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
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
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
Sjostrom, Travis; Daligault, Jérôme
2015-12-01
We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014);] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method. Our results for warm dense aluminum at solid density are inconsistent with the recent experimental results reported by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)]. PMID:26764850
Ionic and electronic transport properties in dense plasmas by orbital-free density functional theory
NASA Astrophysics Data System (ADS)
Sjostrom, Travis; Daligault, Jérôme
2015-12-01
We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014), 10.1103/PhysRevLett.113.155006;] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method. Our results for warm dense aluminum at solid density are inconsistent with the recent experimental results reported by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015), 10.1103/PhysRevLett.115.115001].
Dynamic transitions in molecular dynamics simulations of supercooled silicon
NASA Astrophysics Data System (ADS)
Mei, Xiaojun; Eapen, Jacob
2013-04-01
Two dynamic transitions or crossovers, one at a low temperature (T* ≈ 1006 K) and the other at a high temperature (T0 ≈ 1384 K), are shown to emerge in supercooled liquid silicon using molecular dynamics simulations. The high-temperature transition (T0) marks the decoupling of stress, density, and energy relaxation mechanisms. At the low-temperature transition (T*), depending on the cooling rate, supercooled silicon can either undergo a high-density-liquid to low-density-liquid (HDL-LDL) phase transition or experience an HDL-HDL crossover. Dynamically heterogeneous domains that emerge with supercooling become prominent across the HDL-HDL transition at 1006 K, with well-separated mobile and immobile regions. Interestingly, across the HDL-LDL transition, the most mobile atoms form large prominent aggregates while the least mobile atoms get spatially dispersed akin to that in a crystalline state. The attendant partial return to spatial uniformity with the HDL-LDL phase transition indicates a dynamic mechanism for relieving the frustration in supercooled states.
Atomistic molecular dynamic simulations of multiferroics.
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L
2012-08-10
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings. PMID:23006300
Atomistic Molecular Dynamic Simulations of Multiferroics
NASA Astrophysics Data System (ADS)
Wang, Dawei; Weerasinghe, Jeevaka; Bellaiche, L.
2012-08-01
A first-principles-based approach is developed to simulate dynamical properties, including complex permittivity and permeability in the GHz-THz range, of multiferroics at finite temperatures. It includes both structural degrees of freedom and magnetic moments as dynamic variables in Newtonian and Landau-Lifshitz-Gilbert (LLG) equations within molecular dynamics, respectively, with the couplings between these variables being incorporated. The use of a damping coefficient and of the fluctuation field in the LLG equations is required to obtain equilibrated magnetic properties at any temperature. No electromagnon is found in the spin-canted structure of BiFeO3. On the other hand, two magnons with very different frequencies are predicted via the use of this method. The smallest-in-frequency magnon corresponds to oscillations of the weak ferromagnetic vector in the basal plane being perpendicular to the polarization while the second magnon corresponds to magnetic dipoles going in and out of this basal plane. The large value of the frequency of this second magnon is caused by static couplings between magnetic dipoles with electric dipoles and oxygen octahedra tiltings.
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
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.
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
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.
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.
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.
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
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.
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.
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.
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.
Introducing PROFESS: A new program for orbital-free density functional theory calculations
NASA Astrophysics Data System (ADS)
Ho, Gregory S.; Lignères, Vincent L.; Carter, Emily A.
2008-12-01
We present PROFESS (PRinceton Orbital-Free Electronic Structure Software), a new software package that performs orbital-free density functional theory (OF-DFT) calculations. OF-DFT is a first principles quantum mechanics method primarily for condensed matter that can be made to scale linearly with system size. We describe the implementation of energy, force, and stress functionals and the methods used to optimize the electron density under periodic boundary conditions. All electronic energy and potential terms scale linearly while terms involving the ions exhibit quadratic scaling in our code. Despite the latter scaling, the program can treat tens of thousands of atoms with quantum mechanics on a single processor, as we demonstrate here. Limitations of the method are also outlined, the most serious of which is the accuracy of state-of-the-art kinetic energy functionals, which limits the applicability of the method to main group elements at present. Program summaryProgram title: PROFESS Catalogue identifier: AEBN_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBN_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 35 933 No. of bytes in distributed program, including test data, etc.: 329 924 Distribution format: tar.gz Programming language: Fortran 90 Computer: Intel with ifort; AMD Opteron with pathf90 Operating system: Linux RAM: Problem dependent, but 2 GB is sufficient for up to 10,000 ions Classification: 7.3 External routines: FFTW ( http://www.fftw.org), MINPACK-2 Nature of problem: Given a set of coordinates describing the initial ion positions under periodic boundary conditions, recovers the ground state energy, electron density, ion positions, and cell lattice vectors predicted by orbital-free density functional theory. Except for computation of the
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.
Orbital-free density functional theory at the nanoscale and beyond: Algorithms and applications
NASA Astrophysics Data System (ADS)
Hung, Linda
Orbital-free density functional theory (OFDFT) is a quantum mechanical theory where the total energy of a system depends only on the electron density and an external potential. While not as well-known as its relative, the wavefunction-based Kohn-Sham density functional theory (KSDFT), OFDFT has been established as a fast yet accurate theory for modeling materials -- like aluminum -- where the valence electron density is nearly homogeneous. In this thesis, we describe improvements to the implementation of OFDFT that allow the study of materials at the nanoscale and beyond. We approximate the structure factor that appears in ion-electron and ion-ion terms with B-splines, so that the scaling of their computations is reduced from quadratic to quasilinear ( O(N log(N)), where N is some measure of system size). We also review the treatment of optimization variables, implement a norm-conserving line search, and use preconditioners to improve optimization algorithms. These preconditioners have limited success for simulation cells with regions of low electron density, but can speed up optimization 3-4 times in other cells. Many of these improvements are included in a new version of the PRinceton Orbital-Free Electronic Structure Software (PROFESS). With all energy, potential, force, and stress computations now parallelized and scaling quasilinearly, more than a million Al atoms can be explicitly simulated on 192 cores. Using this PROFESS 2.0 code, we study Al crack tips using a quasi-two-dimensional setup with ~20 nm x 20 nm samples. Although both the embedded-atom method (EAM) and OFDFT can reproduce KSDFT material properties, the mechanism of plasticity at crack tips in simulation can occur through different slip planes, depending on the choice of OFDFT or the EAM. We also model [111]-oriented Al nanowires with 1-8 nm diameter and length up to ~20 nm. During tensile loading from equilibrium, the iv axial displacements of surface atoms (not noticed in past nanowire
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
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
Orbital-free density functional theory study of crystalline Li-Si alloys
NASA Astrophysics Data System (ADS)
Xia, Junchao; Carter, Emily A.
2014-05-01
Li-Si interactions are of great interest currently due to the potential use of silicon anodes in Li-ion batteries. As a first step toward eventual nanoscale characterization of lithiation of silicon, here we study the crystalline Li-Si alloys LiSi, Li12Si7, Li7Si3, Li13Si4, Li15Si4, and Li22Si5 using orbital-free density functional theory (OFDFT). The recently proposed Wang-Govind-Carter decomposition (WGCD) and Huang-Carter (HC) kinetic energy density functionals (KEDFs) are used to evaluate the electron kinetic energy. Both KEDFs predict accurate cell lattice vectors, equilibrium volumes, bulk moduli, and ground-state densities when compared to Kohn-Sham density functional theory (KSDFT) benchmarks. Elastic constants and alloy formation energies calculated with the WGCD KEDF also agree reasonably well with KSDFT. Finally, Li atom adsorption energies on the Si(100) - 2 × 1 surface are calculated as a simple initial test of the Li-Si mixing process during lithiation of silicon. The OFDFT adsorption energies again are fairly close to KSDFT values. The results in this work demonstrate the accuracy of the WGCD and HC KEDFs for materials with mixed covalent-metallic character and their considerable transferability under different chemical environments. Because of its quasilinear scaling, coupled with the level of accuracy shown here, OFDFT appears quite promising for large-scale simulation of such materials phenomena.
Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory
NASA Astrophysics Data System (ADS)
Ghosh, Swarnava; Suryanarayana, Phanish
2016-02-01
We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state and forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms.
ATLAS: A real-space finite-difference implementation of orbital-free density functional theory
NASA Astrophysics Data System (ADS)
Mi, Wenhui; Shao, Xuecheng; Su, Chuanxun; Zhou, Yuanyuan; Zhang, Shoutao; Li, Quan; Wang, Hui; Zhang, Lijun; Miao, Maosheng; Wang, Yanchao; Ma, Yanming
2016-03-01
Orbital-free density functional theory (OF-DFT) is a promising method for large-scale quantum mechanics simulation as it provides a good balance of accuracy and computational cost. Its applicability to large-scale simulations has been aided by progress in constructing kinetic energy functionals and local pseudopotentials. However, the widespread adoption of OF-DFT requires further improvement in its efficiency and robustly implemented software. Here we develop a real-space finite-difference (FD) method for the numerical solution of OF-DFT in periodic systems. Instead of the traditional self-consistent method, a powerful scheme for energy minimization is introduced to solve the Euler-Lagrange equation. Our approach engages both the real-space finite-difference method and a direct energy-minimization scheme for the OF-DFT calculations. The method is coded into the ATLAS software package and benchmarked using periodic systems of solid Mg, Al, and Al3Mg. The test results show that our implementation can achieve high accuracy, efficiency, and numerical stability for large-scale simulations.
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…
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.
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
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.
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
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.}
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.
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.
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.
Improving the orbital-free density functional theory description of covalent materials
NASA Astrophysics Data System (ADS)
Zhou, Baojing; Ligneres, Vincent L.; Carter, Emily A.
2005-01-01
The essential challenge in orbital-free density functional theory (OF-DFT) is to construct accurate kinetic energy density functionals (KEDFs) with general applicability (i.e., transferability). During the last decade, several linear-response (LR)-based KEDFs have been proposed. Among them, the Wang-Govind-Carter (WGC) KEDF, containing a density-dependent response kernel, is one of the most accurate that still affords a linear scaling algorithm. For nearly-free-electron-like metals such as Al and its alloys, OF-DFT employing the WGC KEDF produces bulk properties in good agreement with orbital-based Kohn-Sham (KS) DFT predictions. However, when OF-DFT, using the WGC KEDF combined with a recently proposed bulk-derived local pseudopotential (BLPS), was applied to semiconducting and metallic phases of Si, problems arose with convergence of the self-consistent density and energy, leading to poor results. Here we provide evidence that the convergence problem is very likely caused by the use of a truncated Taylor series expansion of the WGC response kernel. Moreover, we show that a defect in the ansatz for the first-order reduced density matrix underlying the LR KEDFs limits the accuracy of these KEDFs. By optimizing the two free parameters involved in the WGC KEDF, the two-body Fermi wave vector mixing parameter γ and the reference density ρ* used in the Taylor expansion, OF-DFT calculations with the BLPS can achieve semiquantitative results for nine phases of bulk silicon. These new parameters are recommended whenever the WGC KEDF is used to study nonmetallic systems.
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
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.
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.
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
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.
NASA Astrophysics Data System (ADS)
Xia, Junchao
Orbital-free (OF) density functional theory (DFT) is a powerful and numerically efficient first principles quantum mechanics method. Its application has contributed to understanding a diverse set of materials properties in recent decades. However, most previous studies were confined to simple metals. In this thesis, we focus on extending OFDFT to describe covalently-bonded materials and aiming for a balance between accuracy and efficiency. We first apply OFDFT to study diatomic molecules, with the Huang-Carter (HC) kinetic energy density functional (KEDF). OFDFT predicts reasonable equilibrium bond lengths, bond dissociation energies, and vibrational frequencies compared to Kohn-Sham (KS) DFT benchmarks. This work indicates significant progress of OFDFT in describing molecules. However, we find that the HC KEDF is computationally expensive and thus inapplicable for large-scale simulations. Consequently, we propose an electron density decomposition formalism for covalent materials. Based on local density information, the total density is decomposed into localized and delocalized electron densities, which are then described by different KEDF models separately. The resulting Wang--Govind--Carter-decomposition (WGCD) KEDF gives accurate properties for bulk semiconductors and isolated molecules. Furthermore, it offers far superior numerical efficiency compared to the previous HC KEDF. We then test the HC and WGCD KEDFs on Li-Si alloys and obtain accurate structures and bulk properties. The OFDFT Li adsorption energies on the Si(100) surface are also close to KSDFT values. OFDFT is thus promising to study mechanical properties of Li-Si alloys and the mixing mechanism during lithiation and delithiation processes. We next focus on single-point KEDFs for localized densities and pointwise quantities including the local kinetic energy density (KED) and the electron localization function (ELF). Based on a transferable correlation between the reduced density and the KED
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.
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)
Identifying the mechanisms of polymer friction through molecular dynamics simulation.
Dai, Ling; Minn, M; Satyanarayana, N; Sinha, Sujeet K; Tan, V B C
2011-12-20
Mechanisms governing the tribological behavior of polymer-on-polymer sliding were investigated by molecular dynamics simulations. Three main mechanisms governing frictional behavior were identified. Interfacial "brushing" of molecular chain ends over one another was observed as the key contribution to frictional forces. With an increase of the sliding speed, fluctuations in frictional forces reduced in both magnitude and periodicity, leading to dynamic frictional behavior. While "brushing" remained prevalent, two additional irreversible mechanisms, "combing" and "chain scission", of molecular chains were observed when the interfaces were significantly diffused. PMID:22044344
Temperature dependence of protein hydration hydrodynamics by molecular dynamics simulations.
Lau, E Y; Krishnan, V V
2007-07-18
The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.
Reactive Molecular Dynamics Simulations at the Petascale (Invited)
NASA Astrophysics Data System (ADS)
Nakano, A.
2013-12-01
We are developing a divide-conquer-recombine algorithmic framework into a metascalable (or 'design once, scale on new architectures') parallelization scheme to perform large spatiotemporal-scale reactive molecular dynamics simulations. The scheme has achieved parallel efficiency well over 0.9 on 786,432 IBM BlueGene/Q processors for 8.5 trillion-atom molecular dynamics and 1.9 trillion electronic degrees-of-freedom quantum molecular dynamics in the framework of density functional theory. Simulation results reveal intricate interplay between photoexcitation, mechanics, flow, and chemical reactions at the nanoscale. Specifically, we will discuss atomistic mechanisms of: (1) rapid hydrogen production from water using metallic alloy nanoparticles; (2) molecular control of charge transfer, charge recombination, and singlet fission for efficient solar cells; and (3) mechanically enhanced reaction kinetics in nanobubbles and nanojets.
Interfacial Molecular Searching Using Forager Dynamics
NASA Astrophysics Data System (ADS)
Monserud, Jon H.; Schwartz, Daniel K.
2016-03-01
Many biological and technological systems employ efficient non-Brownian intermittent search strategies where localized searches alternate with long flights. Coincidentally, molecular species exhibit intermittent behavior at the solid-liquid interface, where periods of slow motion are punctuated by fast flights through the liquid phase. Single-molecule tracking was used here to observe the interfacial search process of DNA for complementary DNA. Measured search times were qualitatively consistent with an intermittent-flight model, and ˜10 times faster than equivalent Brownian searches, suggesting that molecular searches for reactive sites benefit from similar efficiencies as biological organisms.
Dynamics of molecular superrotors in an external magnetic field
NASA Astrophysics Data System (ADS)
Korobenko, Aleksey; Milner, Valery
2015-08-01
We excite diatomic oxygen and nitrogen to high rotational states with an optical centrifuge and study their dynamics in an external magnetic field. Ion imaging is employed to directly visualize, and follow in time, the rotation plane of the molecular superrotors. The two different mechanisms of interaction between the magnetic field and the molecular angular momentum in paramagnetic oxygen and non-magnetic nitrogen lead to qualitatively different behaviour. In nitrogen, we observe the precession of the molecular angular momentum around the field vector. In oxygen, strong spin-rotation coupling results in faster and richer dynamics, encompassing the splitting of the rotation plane into three separate components. As the centrifuged molecules evolve with no significant dispersion of the molecular wave function, the observed magnetic interaction presents an efficient mechanism for controlling the plane of molecular rotation.
First principles molecular dynamics without self-consistent field optimization
Souvatzis, Petros; Niklasson, Anders M. N.
2014-01-28
We present a first principles molecular dynamics approach that is based on time-reversible extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) construction are required in each integration time step. The proposed dynamics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents a natural starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.
Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber
Luchko, Tyler; Gusarov, Sergey; Roe, Daniel R.; Simmerling, Carlos; Case, David A.; Tuszynski, Jack; Kovalenko, Andriy
2010-01-01
We present the three-dimensional molecular theory of solvation (also known as 3D-RISM) coupled with molecular dynamics (MD) simulation by contracting solvent degrees of freedom, accelerated by extrapolating solvent-induced forces and applying them in large multi-time steps (up to 20 fs) to enable simulation of large biomolecules. The method has been implemented in the Amber molecular modeling package, and is illustrated here on alanine dipeptide and protein G. PMID:20440377
VUV studies of molecular photofragmentation dynamics
White, M.G.
1993-12-01
State-resolved, photoion and photoelectron methods are used to study the neutral fragmentation and ionization dynamics of small molecules relevant to atmospheric and combustion chemistry. Photodissociation and ionization are initiated by coherent VUV radiation and the fragmentation dynamics are extracted from measurements of product rovibronic state distributions, kinetic energies and angular distributions. The general aim of these studies is to investigate the multichannel interactions between the electronic and nuclear motions which determine the evolution of the photoexcited {open_quotes}complex{close_quotes} into the observed asymptotic channels.
Visualizing global properties of a molecular dynamics trajectory.
Zhou, Hao; Li, Shangyang; Makowski, Lee
2016-01-01
Molecular dynamics (MD) trajectories are very large data sets that contain substantial information about the dynamic behavior of a protein. Condensing these data into a form that can provide intuitively useful understanding of the molecular behavior during the trajectory is a substantial challenge that has received relatively little attention. Here, we introduce the sigma-r plot, a plot of the standard deviation of intermolecular distances as a function of that distance. This representation of global dynamics contains within a single, one-dimensional plot, the average range of motion between pairs of atoms within a macromolecule. Comparison of sigma-r plots calculated from 10 ns trajectories of proteins representing the four major SCOP fold classes indicates diversity of dynamic behaviors which are recognizably different among the four classes. Differences in domain structure and molecular weight also produce recognizable features in sigma-r plots, reflective of differences in global dynamics. Plots generated from trajectories with progressively increasing simulation time reflect the increased sampling of the structural ensemble as a function of time. Single amino acid replacements can give rise to changes in global dynamics detectable through comparison of sigma-r plots. Dynamic behavior of substructures can be monitored by careful choice of interatomic vectors included in the calculation. These examples provide demonstrations of the utility of the sigma-r plot to provide a simple measure of the global dynamics of a macromolecule. PMID:26522428
Molecular Mechanotransduction: how forces trigger cytoskeletal dynamics
NASA Astrophysics Data System (ADS)
Ehrlicher, Allen
2012-02-01
Mechanical stresses elicit cellular reactions mediated by chemical signals. Defective responses to forces underlie human medical disorders, such as cardiac failure and pulmonary injury. Despite detailed knowledge of the cytoskeleton's structure, the specific molecular switches that convert mechanical stimuli into chemical signals have remained elusive. Here we identify the actin-binding protein, filamin A (FLNa) as a central mechanotransduction element of the cytoskeleton by using Fluorescence Loss After photoConversion (FLAC), a novel high-speed alternative to FRAP. We reconstituted a minimal system consisting of actin filaments, FLNa and two FLNa-binding partners: the cytoplasmic tail of ß-integrin, and FilGAP. Integrins form an essential mechanical linkage between extracellular and intracellular environments, with ß integrin tails connecting to the actin cytoskeleton by binding directly to filamin. FilGAP is a FLNa-binding GTPase-activating protein specific for Rac, which in vivo regulates cell spreading and bleb formation. We demonstrate that both externally-imposed bulk shear and myosin II driven forces differentially regulate the binding of integrin and FilGAP to FLNa. Consistent with structural predictions, strain increases ß-integrin binding to FLNa, whereas it causes FilGAP to dissociate from FLNa, providing a direct and specific molecular basis for cellular mechanotransduction. These results identify the first molecular mechanotransduction element within the actin cytoskeleton, revealing that mechanical strain of key proteins regulates the binding of signaling molecules. Moreover, GAP activity has been shown to switch cell movement from mesenchymal to amoeboid motility, suggesting that mechanical forces directly impact the invasiveness of cancer.
Molecular dynamics insights into human aquaporin 2 water channel.
Binesh, A R; Kamali, R
2015-12-01
In this study, the first molecular dynamics simulation of the human aquaporin 2 is performed and for a better understanding of the aquaporin 2 permeability performance, the characteristics of water transport in this protein channel and key biophysical parameters of AQP2 tetramer including osmotic and diffusive permeability constants and the pore radius are investigated. For this purpose, recently recovered high resolution X-ray crystal structure of` the human aquaporin 2 is used to perform twenty nanosecond molecular dynamics simulation of fully hydrated tetramer of this protein embedded in a lipid bilayer. The resulting water permeability characteristics of this protein channel showed that the water permeability of the human AQP2 is in a mean range in comparison with other human aquaporins family. Finally, the results reported in this research demonstrate that molecular dynamics simulation of human AQP2 provided useful insights into the mechanisms of water permeation and urine concentration in the human kidney. PMID:26489820
Hydrolysis of Al3+ from constrained molecular dynamics.
Ikeda, Takashi; Hirata, Masaru; Kimura, Takaumi
2006-02-21
We investigated the hydrolysis reactions of Al(3+) in AlCl(3) aqueous solution using the constrained molecular dynamics based on the Car-Parrinello molecular-dynamics method. By employing the proton-aluminum coordination number as a reaction coordinate in the constrained molecular dynamics the deprotonation as well as dehydration processes are successfully realized. From our free-energy difference of DeltaG(0) approximately 8.0 kcal mol(-1) the hydrolysis constant pK(a1) is roughly estimated as 5.8, comparable to the literature value of 5.07. We show that the free-energy difference for the hydrolysis of Al(3+) in acidic conditions is at least 4 kcal mol(-1) higher than that in neutral condition, indicating that the hydrolysis reaction is inhibited by the presence of excess protons located around the hydrated ion, in agreement with the change of the predominant species by pH. PMID:16497053
Theoretical analysis of dynamic processes for interacting molecular motors
NASA Astrophysics Data System (ADS)
Teimouri, Hamid; Kolomeisky, Anatoly B.; Mehrabiani, Kareem
2015-02-01
Biological transport is supported by the collective dynamics of enzymatic molecules that are called motor proteins or molecular motors. Experiments suggest that motor proteins interact locally via short-range potentials. We investigate the fundamental role of these interactions by carrying out an analysis of a new class of totally asymmetric exclusion processes, in which interactions are accounted for in a thermodynamically consistent fashion. This allows us to explicitly connect microscopic features of motor proteins with their collective dynamic properties. A theoretical analysis that combines various mean-field calculations and computer simulations suggests that the dynamic properties of molecular motors strongly depend on the interactions, and that the correlations are stronger for interacting motor proteins. Surprisingly, it is found that there is an optimal strength of interactions (weak repulsion) that leads to a maximal particle flux. It is also argued that molecular motor transport is more sensitive to attractive interactions. Applications of these results for kinesin motor proteins are discussed.
Studying Interactions by Molecular Dynamics Simulations at High Concentration
Fogolari, Federico; Corazza, Alessandra; Toppo, Stefano; Tosatto, Silvio C. E.; Viglino, Paolo; Ursini, Fulvio; Esposito, Gennaro
2012-01-01
Molecular dynamics simulations have been used to study molecular encounters and recognition. In recent works, simulations using high concentration of interacting molecules have been performed. In this paper, we consider the practical problems for setting up the simulation and to analyse the results of the simulation. The simulation of beta 2-microglobulin association and the simulation of the binding of hydrogen peroxide by glutathione peroxidase are provided as examples. PMID:22500085
Replica exchange molecular dynamics simulations of amyloid peptide aggregation
NASA Astrophysics Data System (ADS)
Cecchini, M.; Rao, F.; Seeber, M.; Caflisch, A.
2004-12-01
The replica exchange molecular dynamics (REMD) approach is applied to four oligomeric peptide systems. At physiologically relevant temperature values REMD samples conformation space and aggregation transitions more efficiently than constant temperature molecular dynamics (CTMD). During the aggregation process the energetic and structural properties are essentially the same in REMD and CTMD. A condensation stage toward disordered aggregates precedes the β-sheet formation. Two order parameters, borrowed from anisotropic fluid analysis, are used to monitor the aggregation process. The order parameters do not depend on the peptide sequence and length and therefore allow to compare the amyloidogenic propensity of different peptides.
Trillion-atom molecular dynamics becomes a reality
Kadau, Kai; Germann, Timothy C
2008-01-01
By utilizing the molecular dynamics code SPaSM on Livermore's BlueGene/L architecture, consisting of 212 992 IBM PowerPC440 700 MHz processors, a molecular dynamics simulation was run with one trillion atoms. To demonstrate the practicality and future potential of such ultra large-scale simulations, the onset of the mechanical shear instability occurring in a system of Lennard-Jones particles arranged in a simple cubic lattice was simulated. The evolution of the instability was analyzed on-the-fly using the in-house developed massively parallel graphical object-rendering code MD{_}render.
State-to-state dynamics of molecular energy transfer
Gentry, W.R.; Giese, C.F.
1993-12-01
The goal of this research program is to elucidate the elementary dynamical mechanisms of vibrational and rotational energy transfer between molecules, at a quantum-state resolved level of detail. Molecular beam techniques are used to isolate individual molecular collisions, and to control the kinetic energy of collision. Lasers are used both to prepare specific quantum states prior to collision by stimulated-emission pumping (SEP), and to measure the distribution of quantum states in the collision products by laser-induced fluorescence (LIF). The results are interpreted in terms of dynamical models, which may be cast in a classical, semiclassical or quantum mechanical framework, as appropriate.
Energy conserving, linear scaling Born-Oppenheimer molecular dynamics.
Cawkwell, M J; Niklasson, Anders M N
2012-10-01
Born-Oppenheimer molecular dynamics simulations with long-term conservation of the total energy and a computational cost that scales linearly with system size have been obtained simultaneously. Linear scaling with a low pre-factor is achieved using density matrix purification with sparse matrix algebra and a numerical threshold on matrix elements. The extended Lagrangian Born-Oppenheimer molecular dynamics formalism [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] yields microcanonical trajectories with the approximate forces obtained from the linear scaling method that exhibit no systematic drift over hundreds of picoseconds and which are indistinguishable from trajectories computed using exact forces. PMID:23039583
Mitochondrial dynamics: molecular mechanisms and the role in the heart.
Jazbutyte, V
2010-04-01
Mitochondria are dynamic organelles which actively move along the cytoskeleton within the cell, change their shape and undergo fusion and fission. The heart is a metabolically active organ with high energy demands and rich in mitochondria. Mitochondria not only supply the heart with the high energy compound, adenosine triphosphate (ATP), but also actively participate in cell signaling and apoptotic events and communicate with the cytosol. Recent advantages in molecular biology and imaging techniques helped to study mitochondrial dynamics directly in the cell and under real time conditions. In this review, I will briefly summarize current knowledge about molecular machinery mediating mitochondrial fusion/ fission, its link to apoptosis and cardiovascular disease. PMID:20440252
How Dynamic Visualization Technology can Support Molecular Reasoning
NASA Astrophysics Data System (ADS)
Levy, Dalit
2012-11-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 gas. They interact with the visualizations and carry out inquiry activities to make and refine connections between observable phenomena and atomic level processes related to phase change. The explanations proposed by 300 pairs of students in response to pre/post-assessment items have been analyzed using a scale for measuring the level of molecular reasoning. Results indicate that from pretest to posttest, students make progress in their level of molecular reasoning and are better able to connect intermolecular forces and phase change in their explanations. The paper presents the results through the lens of improvement patterns and the metaphor of the "ladder of molecular reasoning," and discusses how this adds to our understanding of the benefits of interacting with dynamic molecular visualizations.
Scaling properties of the kinetic energy density of atoms - towards an orbital-free meta-GGA
NASA Astrophysics Data System (ADS)
Cancio, Antonio; Redd, Jeremy
The scaling properties of atoms, combining periodicity with gradual increase in density, make a fruitful probe of relationships in density functional theory, and have driven advances in understanding the exchange and correlation energy. Although focus is normally upon the properties of integrated energies, insights can be generated from studying energy density functions as well. We visualize the behavior of the positive-definite kinetic energy density (KED) in closed-shell atoms, in comparison to invariant quantities based upon the gradient and Laplacian of the density. The latter are potential variables for constructing orbital-free functionals for the KE and can be used for analyzing the electronic structure of atoms and molecules. We notice a striking fit of the KED within the core of any atom to a gradient expansion model using both the gradient and the Laplacian, but one different from that derived from first principles for a slowly-varying electron gas. Correlated with this feature, we notice unexpected structure to the KED near the nucleus that cannot be explained simply by the von Weizsacker model, as is often presumed. These unexpected features provide potential insights for developing better orbital-free meta-GGA models for the kinetic energy.
Multiple time step integrators in ab initio molecular dynamics
Luehr, Nathan; Martínez, Todd J.; Markland, Thomas E.
2014-02-28
Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.
Enhanced Sampling Techniques in Molecular Dynamics Simulations of Biological Systems
Bernardi, Rafael C.; Melo, Marcelo C. R.; Schulten, Klaus
2014-01-01
Background Molecular Dynamics has emerged as an important research methodology covering systems to the level of millions of atoms. However, insufficient sampling often limits its application. The limitation is due to rough energy landscapes, with many local minima separated by high-energy barriers, which govern the biomolecular motion. Scope of review In the past few decades methods have been developed that address the sampling problem, such as replica-exchange molecular dynamics, metadynamics and simulated annealing. Here we present an overview over theses sampling methods in an attempt to shed light on which should be selected depending on the type of system property studied. Major Conclusions Enhanced sampling methods have been employed for a broad range of biological systems and the choice of a suitable method is connected to biological and physical characteristics of the system, in particular system size. While metadynamics and replica-exchange molecular dynamics are the most adopted sampling methods to study biomolecular dynamics, simulated annealing is well suited to characterize very flexible systems. The use of annealing methods for a long time was restricted to simulation of small proteins; however, a variant of the method, generalized simulated annealing, can be employed at a relatively low computational cost to large macromolecular complexes. General Significance Molecular dynamics trajectories frequently do not reach all relevant conformational substates, for example those connected with biological function, a problem that can be addressed by employing enhanced sampling algorithms. PMID:25450171
Chain networking revealed by molecular dynamics simulation
NASA Astrophysics Data System (ADS)
Zheng, Yexin; Tsige, Mesfin; Wang, Shi-Qing
Based on Kremer-Grest model for entangled polymer melts, we demonstrate how the response of a polymer glass depends critically on the chain length. After quenching two melts of very different chain lengths (350 beads per chain and 30 beads per chain) into deeply glassy states, we subject them to uniaxial extension. Our MD simulations show that the glass of long chains undergoes stable necking after yielding whereas the system of short chains is unable to neck and breaks up after strain localization. During ductile extension of the polymer glass made of long chain significant chain tension builds up in the load-bearing strands (LBSs). Further analysis is expected to reveal evidence of activation of the primary structure during post-yield extension. These results lend support to the recent molecular model 1 and are the simulations to demonstrate the role of chain networking. This work is supported, in part, by a NSF Grant (DMR-EAGER-1444859)
Electron-phonon interaction within classical molecular dynamics
NASA Astrophysics Data System (ADS)
Tamm, A.; Samolyuk, G.; Correa, A. A.; Klintenberg, M.; Aabloo, A.; Caro, A.
2016-07-01
We present a model for nonadiabatic classical molecular dynamics simulations that captures with high accuracy the wave-vector q dependence of the phonon lifetimes, in agreement with quantum mechanics calculations. It is based on a local view of the e -ph interaction where individual atom dynamics couples to electrons via a damping term that is obtained as the low-velocity limit of the stopping power of a moving ion in a host. The model is parameter free, as its components are derived from ab initio-type calculations, is readily extended to the case of alloys, and is adequate for large-scale molecular dynamics computer simulations. We also show how this model removes some oversimplifications of the traditional ionic damped dynamics commonly used to describe situations beyond the Born-Oppenheimer approximation.
Electron-phonon interaction within classical molecular dynamics
Tamm, A.; Samolyuk, G.; Correa, A. A.; Klintenberg, M.; Aabloo, A.; Caro, A.
2016-07-14
Here, we present a model for nonadiabatic classical molecular dynamics simulations that captures with high accuracy the wave-vector q dependence of the phonon lifetimes, in agreement with quantum mechanics calculations. It is based on a local view of the e-ph interaction where individual atom dynamics couples to electrons via a damping term that is obtained as the low-velocity limit of the stopping power of a moving ion in a host. The model is parameter free, as its components are derived from ab initio-type calculations, is readily extended to the case of alloys, and is adequate for large-scale molecular dynamics computermore » simulations. We also show how this model removes some oversimplifications of the traditional ionic damped dynamics commonly used to describe situations beyond the Born-Oppenheimer approximation.« less
Numerical methods for molecular dynamics. Progress report
Skeel, R.D.
1991-12-31
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
Shapiro like steps reveals molecular nanomagnets' spin dynamics
NASA Astrophysics Data System (ADS)
Abdollahipour, Babak; Abouie, Jahanfar; Ebrahimi, Navid
2015-09-01
We present an accurate way to detect spin dynamics of a nutating molecular nanomagnet by inserting it in a tunnel Josephson junction and studying the current voltage (I-V) characteristic. The spin nutation of the molecular nanomagnet is generated by applying two circularly polarized magnetic fields. We demonstrate that modulation of the Josephson current by the nutation of the molecular nanomagnet's spin appears as a stepwise structure like Shapiro steps in the I-V characteristic of the junction. Width and heights of these Shapiro-like steps are determined by two parameters of the spin nutation, frequency and amplitude of the nutation, which are simply tuned by the applied magnetic fields.
Imaging the molecular dynamics of dissociative electron attachment to water
Adaniya, Hidihito; Rudek, B.; Osipov, Timur; Haxton, Dan; Weber, Thorsten; Rescigno, Thomas N.; McCurdy, C.W.; Belkacem, Ali
2009-10-19
Momentum imaging experiments on dissociative electron attachment to the water molecule are combined with ab initio theoretical calculations of the angular dependence of the quantum mechanical amplitude for electron attachment to provide a detailed picture of the molecular dynamics of dissociation attachment via the two lowest energy Feshbach resonances. The combination of momentum imaging experiments and theory can reveal dissociation dynamics for which the axial recoil approximation breaks down and thus provides a powerful reaction microscope for DEA to polyatomics.
Coarse-Grained Molecular Dynamics: Dissipation Due to Internal Modes
Rudd, R E
2001-12-21
We describe progress on the issue of pathological elastic wave reflection in atomistic and multiscale simulation. First we briefly review Coarse-Grained Molecular Dynamics (CGMD). Originally CGMD was formulated as a Hamiltonian system in which energy is conserved. This formulation is useful for many applications, but recently CGMD has been extended to include generalized Langevin forces. Here we describe how Langevin dynamics arise naturally in CGMD, and we examine the implication for elastic wave scattering.
Input File Creation for the Molecular Dynamics Program LAMMPS.
Energy Science and Technology Software Center (ESTSC)
2001-05-30
The program creates an input data file for the molecular dynamics program LAMMPS. The input file created is a liquid mixture between two walls explicitly composed of particles. The liquid molecules are modeled as a bead-spring molecule. The input data file specifies the position and topology of the starting state. The data structure of input allows for dynamic bond creation (cross-linking) within the LAMMPS code.
Electron trapping in amorphous silicon: A quantum molecular dynamics study
Yang, Lin H.; Kalia, R.K.; Vashishta, P.
1990-12-01
Quantum molecular dynamics (QMD) simulations provide the real-time dynamics of electrons and ions through numerical solutions of the time-dependent Schrodinger and Newton equations, respectively. Using the QMD approach we have investigated the localization behavior of an excess electron in amorphous silicon at finite temperatures. For time scales on the order of a few picoseconds, we find the excess electron is localized inside a void of radius {approximately}3 {Angstrom} at finite temperatures. 12 refs.
Plastic dislocation motion via nonequilibrium molecular and continuum dynamics
Hoover, W.G.; Ladd, A.J.C.; Hoover, N.E.
1980-09-29
The classical two-dimensional close-packed triangular lattice, with nearest-neighbor spring forces, is a convenient standard material for the investigation of dislocation motion and plastic flow. Two kinds of calculations, based on this standard material, are described here: (1) Molecular Dynamics simulations, incorporating adiabatic strains described with the help of Doll's Tensor, and (2) Continuum Dynamics simulations, incorporating periodic boundaries and dislocation interaction through stress-field superposition.
A fast recursive algorithm for molecular dynamics simulation
NASA Technical Reports Server (NTRS)
Jain, A.; Vaidehi, N.; Rodriguez, G.
1993-01-01
The present recursive algorithm for solving molecular systems' dynamical equations of motion employs internal variable models that reduce such simulations' computation time by an order of magnitude, relative to Cartesian models. Extensive use is made of spatial operator methods recently developed for analysis and simulation of the dynamics of multibody systems. A factor-of-450 speedup over the conventional O(N-cubed) algorithm is demonstrated for the case of a polypeptide molecule with 400 residues.
NASA Astrophysics Data System (ADS)
Shimizu, Futoshi; Kimizuka, Hajime; Kaburaki, Hideo
2002-08-01
A new parallel computing environment, called as ``Parallel Molecular Dynamics Stencil'', has been developed to carry out a large-scale short-range molecular dynamics simulation of solids. The stencil is written in C language using MPI for parallelization and designed successfully to separate and conceal parts of the programs describing cutoff schemes and parallel algorithms for data communication. This has been made possible by introducing the concept of image atoms. Therefore, only a sequential programming of the force calculation routine is required for executing the stencil in parallel environment. Typical molecular dynamics routines, such as various ensembles, time integration methods, and empirical potentials, have been implemented in the stencil. In the presentation, the performance of the stencil on parallel computers of Hitachi, IBM, SGI, and PC-cluster using the models of Lennard-Jones and the EAM type potentials for fracture problem will be reported.
Superspreading: molecular dynamics simulations and experimental results
NASA Astrophysics Data System (ADS)
Theodorakis, Panagiotis; Kovalchuk, Nina; Starov, Victor; Muller, Erich; Craster, Richard; Matar, Omar
2015-11-01
The intriguing ability of certain surfactant molecules to drive the superspreading of liquids to complete wetting on hydrophobic substrates is central to numerous applications that range from coating flow technology to enhanced oil recovery. Recently, we have observed that for superspreading to occur, two key conditions must be simultaneously satisfied: the adsorption of surfactants from the liquid-vapor surface onto the three-phase contact line augmented by local bilayer formation. Crucially, this must be coordinated with the rapid replenishment of liquid-vapor and solid-liquid interfaces with surfactants from the interior of the droplet. Here, we present the structural characteristics and kinetics of the droplet spreading during the different stages of this process, and we compare our results with experimental data for trisiloxane and poly oxy ethylene surfactants. In this way, we highlight and explore the differences between surfactants, paving the way for the design of molecular architectures tailored specifically for applications that rely on the control of wetting. EPSRC Platform Grant MACIPh (EP/L020564/).
Molecular circuits for dynamic noise filtering.
Zechner, Christoph; Seelig, Georg; Rullan, Marc; Khammash, Mustafa
2016-04-26
The invention of the Kalman filter is a crowning achievement of filtering theory-one that has revolutionized technology in countless ways. By dealing effectively with noise, the Kalman filter has enabled various applications in positioning, navigation, control, and telecommunications. In the emerging field of synthetic biology, noise and context dependency are among the key challenges facing the successful implementation of reliable, complex, and scalable synthetic circuits. Although substantial further advancement in the field may very well rely on effectively addressing these issues, a principled protocol to deal with noise-as provided by the Kalman filter-remains completely missing. Here we develop an optimal filtering theory that is suitable for noisy biochemical networks. We show how the resulting filters can be implemented at the molecular level and provide various simulations related to estimation, system identification, and noise cancellation problems. We demonstrate our approach in vitro using DNA strand displacement cascades as well as in vivo using flow cytometry measurements of a light-inducible circuit in Escherichia coli. PMID:27078094
Probing Molecular Dynamics by Laser-Induced Backscattering Holography.
Haertelt, Marko; Bian, Xue-Bin; Spanner, Michael; Staudte, André; Corkum, Paul B
2016-04-01
We use differential holography to overcome the forward scattering problem in strong-field photoelectron holography. Our differential holograms of H_{2} and D_{2} molecules exhibit a fishbonelike structure, which arises from the backscattered part of the recolliding photoelectron wave packet. We demonstrate that the backscattering hologram can resolve the different nuclear dynamics between H_{2} and D_{2} with subangstrom spatial and subcycle temporal resolution. In addition, we show that attosecond electron dynamics can be resolved. These results open a new avenue for ultrafast studies of molecular dynamics in small molecules. PMID:27081975
Probing Molecular Dynamics by Laser-Induced Backscattering Holography
NASA Astrophysics Data System (ADS)
Haertelt, Marko; Bian, Xue-Bin; Spanner, Michael; Staudte, André; Corkum, Paul B.
2016-04-01
We use differential holography to overcome the forward scattering problem in strong-field photoelectron holography. Our differential holograms of H2 and D2 molecules exhibit a fishbonelike structure, which arises from the backscattered part of the recolliding photoelectron wave packet. We demonstrate that the backscattering hologram can resolve the different nuclear dynamics between H2 and D2 with subangstrom spatial and subcycle temporal resolution. In addition, we show that attosecond electron dynamics can be resolved. These results open a new avenue for ultrafast studies of molecular dynamics in small molecules.
Ultrafast dynamics in isolated molecules and molecular clusters
NASA Astrophysics Data System (ADS)
Hertel, I. V.; Radloff, W.
2006-06-01
During the past decade the understanding of photo-induced ultrafast dynamics in molecular systems has improved at an unforeseen speed and a wealth of detailed insight into the fundamental processes has been obtained. This review summarizes our present knowledge on ultrafast dynamics in isolated molecules and molecular clusters evolving after excitation with femtosecond pulses as studied by pump-probe analysis in real time. Experimental tools and methods as well as theoretical models are described which have been developed to glean information on primary, ultrafast processes in photophysics, photochemistry and photobiology. The relevant processes are explained by way of example—from wave packet dynamics in systems with a few atoms all the way to internal conversion via conical intersections in bio-chromophores. A systematic overview on characteristic systems follows, starting with diatomic and including larger organic molecules as well as various types of molecular clusters, such as micro-solvated chromophore molecules. For conciseness the focus is on molecular systems which remain unperturbed by the laser pulses—apart from the excitation and detection processes as such. Thus, only some aspects of controlling and manipulating molecular reactions by shaped and/or very intense laser pulses are discussed briefly for particularly instructive examples, illustrating the perspectives of this prospering field. The material presented in this review comprises some prototypical examples from earlier pioneering work but emphasizes studies from recent years and covers the most important and latest developments until January 2006.
Optimizing legacy molecular dynamics software with directive-based offload
NASA Astrophysics Data System (ADS)
Michael Brown, W.; Carrillo, Jan-Michael Y.; Gavhane, Nitin; Thakkar, Foram M.; Plimpton, Steven J.
2015-10-01
Directive-based programming models are one solution for exploiting many-core coprocessors to increase simulation rates in molecular dynamics. They offer the potential to reduce code complexity with offload models that can selectively target computations to run on the CPU, the coprocessor, or both. In this paper, we describe modifications to the LAMMPS molecular dynamics code to enable concurrent calculations on a CPU and coprocessor. We demonstrate that standard molecular dynamics algorithms can run efficiently on both the CPU and an x86-based coprocessor using the same subroutines. As a consequence, we demonstrate that code optimizations for the coprocessor also result in speedups on the CPU; in extreme cases up to 4.7X. We provide results for LAMMPS benchmarks and for production molecular dynamics simulations using the Stampede hybrid supercomputer with both Intel® Xeon Phi™ coprocessors and NVIDIA GPUs. The optimizations presented have increased simulation rates by over 2X for organic molecules and over 7X for liquid crystals on Stampede. The optimizations are available as part of the "Intel package" supplied with LAMMPS.
Reasoning with Atomic-Scale Molecular Dynamic Models
ERIC Educational Resources Information Center
Pallant, Amy; Tinker, Robert F.
2004-01-01
The studies reported in this paper are an initial effort to explore the applicability of computational models in introductory science learning. Two instructional interventions are described that use a molecular dynamics model embedded in a set of online learning activities with middle and high school students in 10 classrooms. The studies indicate…
Nanoscale probing of dynamics in local molecular environments.
Atkin, Joanna M; Sass, Paul M; Teichen, Paul E; Eaves, Joel D; Raschke, Markus B
2015-11-19
Vibrational spectroscopy can provide information about structure, coupling, and dynamics underlying the properties of complex molecular systems. While measurements of spectral line broadening can probe local chemical environments, the spatial averaging in conventional spectroscopies limits insight into underlying heterogeneity, in particular in disordered molecular solids. Here, using femtosecond infrared scattering scanning near-field optical microscopy (IR s-SNOM), we resolve in vibrational free-induction decay (FID) measurements a high degree of spatial heterogeneity in polytetrafluoroethylene (PTFE) as a dense molecular model system. In nanoscopic probe volumes as small as 10(3) vibrational oscillators, we approach the homogeneous response limit, with extended vibrational dephasing times of several picoseconds, that is, up to 10 times the inhomogeneous lifetime, and spatial average converging to the bulk ensemble response. We simulate the dynamics of relaxation with a finite set of local vibrational transitions subject to random modulations in frequency. The combined results suggest that the observed heterogeneity arises due to static and dynamic variations in the local molecular environment. This approach thus provides real-space and real-time visualization of the subensemble dynamics that define the properties of many functional materials. PMID:26528865
Relating Soil Organic Matter Dynamics to its Molecular Structure
Technology Transfer Automated Retrieval System (TEKTRAN)
Our understanding of the dynamics of soil organic matter (SOM) must be integrated with a sound knowledge of it biochemical complexity. The molecular structure of SOM was determined in 98% sand soils to eliminate the known protective effects of clay on the amount and turnover rate of the SOM constitu...
Clustering Molecular Dynamics Trajectories for Optimizing Docking Experiments
De Paris, Renata; Quevedo, Christian V.; Ruiz, Duncan D.; Norberto de Souza, Osmar; Barros, Rodrigo C.
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
Clustering molecular dynamics trajectories for optimizing docking experiments.
De Paris, Renata; Quevedo, Christian V; Ruiz, Duncan D; Norberto de Souza, Osmar; Barros, Rodrigo C
2015-01-01
Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for the k-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand. PMID:25873944
Molecular vibrational dynamics in PMMA studied by femtosecond CARS
NASA Astrophysics Data System (ADS)
Zhao, Yang; Zhang, Sheng; Zhou, Boyang; Fan, Rongwei; Chen, Deying; Zhang, Zhonghua; Xia, Yuanqin
2014-11-01
The ultrafast molecular vibrational dynamics in PMMA sheets is studied by femtosecond time-resolved coherent anti-Stokes Raman spectroscopy at room temperature. The C-H stretch modes at 2870 cm-1 and 3008 cm-1 in PMMA sheets are excited and detected. The coherence relaxation times and beat wavenumbers of the Raman modes are obtained.