Abstract
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.
Original language | English (US) |
---|---|
Article number | 22A525 |
Journal | Journal of Chemical Physics |
Volume | 137 |
Issue number | 22 |
DOIs | |
State | Published - Dec 14 2012 |
ASJC Scopus subject areas
- Physics and Astronomy(all)
- Physical and Theoretical Chemistry
Access to Document
Other files and links
Link to the citations in Scopus
(Video) DOE CSGF 2020: Multiscale Reactive Molecular Dynamics Applied to Influenza A M2 Reveals Critical...
Fingerprint
Dive into the research topics of 'Multiscale reactive molecular dynamics'. Together they form a unique fingerprint.
View full fingerprint
Cite this
- APA
- Standard
- Harvard
- Vancouver
- Author
- BIBTEX
- RIS
Knight, C., Lindberg, G. E., & Voth, G. A. (2012). Multiscale reactive molecular dynamics. Journal of Chemical Physics, 137(22), [22A525]. https://doi.org/10.1063/1.4743958
Multiscale reactive molecular dynamics. / Knight, Chris; Lindberg, Gerrick E.; Voth, Gregory A.
In: Journal of Chemical Physics, Vol. 137, No. 22, 22A525, 14.12.2012.
Research output: Contribution to journal › Article › peer-review
Knight, C, Lindberg, GE & Voth, GA 2012, 'Multiscale reactive molecular dynamics', Journal of Chemical Physics, vol. 137, no. 22, 22A525. https://doi.org/10.1063/1.4743958
Knight C, Lindberg GE, Voth GA. Multiscale reactive molecular dynamics. Journal of Chemical Physics. 2012 Dec 14;137(22):22A525. doi: https://doi.org/10.1063/1.4743958
Knight, Chris ; Lindberg, Gerrick E. ; Voth, Gregory A. / Multiscale reactive molecular dynamics. In: Journal of Chemical Physics. 2012 ; Vol. 137, No. 22.
@article{44d0a46ef9c34c3c9de22949100f7b1a,
title = "Multiscale reactive molecular dynamics",
abstract = "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.",
author = "Chris Knight and Lindberg, {Gerrick E.} and Voth, {Gregory A.}",
note = "Funding Information: This research was supported in part by the National Science Foundation (NSF, Grant No. CHE-1214087) and the National Institutes of Health (NIH, Grant No. R01-GM053148). We also acknowledge support from the Department of Defense Multidisciplinary University Research Initiative through the (U.S.) Army Research Office (USARO Grant No. W911NF-10-1-0520), the (U.S.) Department of Energy (DOE) under Contract No. DE-AC02-06CH11357, and an Argonne National Laboratory (ANL) Computational Science Postdoctoral Fellowship to C.K. The computations in this work were supported in part by a grant of computer time from the U.S. Department of Defense (DOD) High Performance Computing Modernization Program at the Navy, Engineer Research and Development Center (ERDC), and Air Force Research Laboratory (AFRL) DOD Supercomputing Resource Centers. These computations were also supported in part by the National Science Foundation Teragrid and Extreme Science and Engineering Discovery Environment (XSEDE) computing resources provided by the Texas Advanced Computing Center under Grant No. MCA94P017.",
year = "2012",
month = dec,
day = "14",
doi = "https://doi.org/10.1063/1.4743958",
language = "English (US)",
volume = "137",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "22",
}
TY - JOUR
T1 - Multiscale reactive molecular dynamics
AU - Knight, Chris
AU - Lindberg, Gerrick E.
AU - Voth, Gregory A.
N1 - Funding Information: This research was supported in part by the National Science Foundation (NSF, Grant No. CHE-1214087) and the National Institutes of Health (NIH, Grant No. R01-GM053148). We also acknowledge support from the Department of Defense Multidisciplinary University Research Initiative through the (U.S.) Army Research Office (USARO Grant No. W911NF-10-1-0520), the (U.S.) Department of Energy (DOE) under Contract No. DE-AC02-06CH11357, and an Argonne National Laboratory (ANL) Computational Science Postdoctoral Fellowship to C.K. The computations in this work were supported in part by a grant of computer time from the U.S. Department of Defense (DOD) High Performance Computing Modernization Program at the Navy, Engineer Research and Development Center (ERDC), and Air Force Research Laboratory (AFRL) DOD Supercomputing Resource Centers. These computations were also supported in part by the National Science Foundation Teragrid and Extreme Science and Engineering Discovery Environment (XSEDE) computing resources provided by the Texas Advanced Computing Center under Grant No. MCA94P017.
PY - 2012/12/14
Y1 - 2012/12/14
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=84870938810&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84870938810&partnerID=8YFLogxK
U2 - https://doi.org/10.1063/1.4743958
DO - https://doi.org/10.1063/1.4743958
M3 - Article
C2 - 23249062
SN - 0021-9606
VL - 137
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
IS - 22
M1 - 22A525
ER -
FAQs
What is the largest MD simulation? ›
The largest system that contains 1.6 billion atoms was simulated using MD with a performance of 8.30 ns/day on Fugaku supercomputer. It extends the available size and time of MD simulations to answer unresolved questions of biomacromolecules in a living cell.
What are the limitations of MD simulations? ›A limitation of classical force field-based MD is the restriction to covalent complexes, with exclusion of chemical reactions. The very important appli- cations to reactions in the condensed phase, including enzyme reactions and catalysis in general, need extension with dynamic behaviour of non-covalent intermediates.
Is molecular dynamics accurate? ›Density functional theory (DFT)-based ab-initio molecular dynamics (AIMD) is accurate but computational cost limits its applications to small systems.
What is the time scale for MD simulation? ›To ensure numerical stability, the time steps in an MD simulation must be short, typically only a few femtoseconds (10–15 s) each. Most of the events of biochemical interest—for example, functionally important structural changes in proteins—take place on timescales of nanoseconds, microseconds, or longer.
Why is MD more competitive than do? ›While this may come as a shock to you, the primary reason behind this fact is that there are far fewer accredited DO programs (37) than accredited MD programs (155) in the U.S. In other words, because more MD programs exist, you are statistically more likely to get into an MD program vs. a DO program.
How many steps are in MD simulation? ›“Number of steps for the MD simulation”: 50000. “EM tolerance”: 1000. “Maximum step size”: 0.01. “Generate detailed log”: Yes.
What are at least 3 disadvantages of modelling and simulation? ›- Mistakes may be made in the programming or rules of the simulation or model.
- The cost of a simulation model can be high.
- The cost of running several different simulations may be high.
- Time may be needed to make sense of the results.
Simulation is not precise. It is not an optimization process and does not yield an answer but merely provides a set of the system's responses to different operating conditions. In many cases this lack of precision is difficult to measure. A good simulation model may be very expensive.
What is the weakness of simulation training? ›There is always room for error. No matter how accurate the learning simulation is, there is always some scope for error and doubt when it comes to the re-creation of real life scenarios. The biggest drawback of using simulation is maintenance and updates can be costly.
What is a good RMSD value molecular dynamics? ›An RMSD value ≤ 2 Å is fairly good.
What is a good RMSD in molecular dynamics? ›
rmsd between 2 and 3 Å is normal/average* rmsd between 3 and 4 Å means differences in protein folding. rmsd higher than 4 means problems!
Which molecular drawing is the most accurate? ›A 3-D space-filling model shows connectivity of atoms and the shape of the molecule. These work well for smaller molecules, but can be a challenge to interpret since atoms are so close together. This is the most accurate representation of what a molecule looks like.
What is the typical time step for molecular dynamics simulation? ›time step size is about 0.0333 to 0.01 of the smallest vibrational period in the simulation.
What is MD vs MC simulation? ›Both Classical Monte Carlo (MC) and Classical Molecular Dynamics (MD) simulations are used to perform simulations of ensembles of molecules. These MC calculations are calculating thermodynamic properties via an ensemble average, while the MD simulations are doing so via a time average.
How do you continue MD simulation? ›- Call the program with gmx.
- Select the mdrun command.
- The -v makes the command verbose. It is useful to make clearer what the program is actually doing while running.
- The -deffnm option is followed by the prefix of the tpr file we are using. ...
- We will continue the simulation from the md.
Do MDs Look Down on DOs? In practice, DOs and MDs work side by side and are respected equally by the majority of those in medicine. The consensus in most hospitals and residency programs is that they don't care if you're a DO or MD.
What is the easiest MD to become? ›- Family Medicine. Average Step 1 Score: 216.1. ...
- Psychiatry. Average Step 1 Score: 223.1. ...
- Physical Medicine and Rehabilitation. Average Step 1 Score: 224.2. ...
- Pediatrics. Average Step 1 Score: 226.4. ...
- Internal Medicine (Categorical) Average Step 1 Score: 231.4. ...
- Anesthesiology.
In the United States, an MD degree is typically more well-respected than a DO. That does not mean a physician with either degree is actually better or worse than the other. There are more MDs than DOs, and because of this standardized acceptance of MDs, they often are considered slightly more reputable.
What are the 5 steps of a simulation? ›Step 1→ | Define the problem or system you intended to simulate. |
---|---|
Step 2→ | Formulate the model you intend to use. |
Step 3→ | Test the model; compare its behaviour with the behaviour of the actual problem. |
Step 4→ | Identify and collect the data needed to test the model. |
Step 5→ | Run the simulation |
Molecular dynamics (MD) and related methods are close to becoming routine computational tools for drug discovery. Their main advantage is in explicitly treating structural flexibility and entropic effects.
What are MD trajectories? ›
Abstract. 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.
What are the 4 types of models in simulation? ›- Monte Carlo method.
- Agent-based modeling.
- Discrete event simulation.
- System dynamic modeling.
Simulation requires manpower and it is a time-consuming process. Simulation results are difficult to translate. It requires experts to understand. Simulation process is expensive.
What are the three methods in simulation modeling? ›To date, there exist three methods: System Dynamics • Discrete Event Modeling • Agent-Based Modeling The choice of method should be based on the system being modeled and the purpose of the modeling – though often it is most heavily influenced by the background or available tool set of the modeler.
When should you avoid simulation? ›Waiting to Simulate
Perhaps the biggest mistake is to wait until the design is complete before simulating. If simulation is only performed late in the design cycle then it is almost not worth doing. Instead of accelerating the process, it actually delays the design process.
SUMMARY OF FAILURES
Failure to define an achievable goal 2. Incomplete mix of essential skills 3. Inadequate level of user participation 4. Inappropriate levels of detail -^ , 5.
Typical reasons why simulation projects fail include the following: Failure to state clear objectives at the outset. Failure to involve individuals affected by outcome. Overrunning budget and time constraints.
What are the five common pitfalls to successful simulation? ›Users of simulation methods might encounter the following five pitfalls: distraction, complexity, implementation, interpretation, and acceptance.
What are two limitations of the simulation? ›- A good simulation model may be very expensive. Simulation often requires a significant amount of computer time and is therefore expensive.
- Simulation generates a way of evaluating solutions but does not generate solutions themselves.
...
Main disadvantages of simulation include:
- Expensive to build a simulation model.
- Expensive to conduct simulation.
- Sometimes it is difficult to interpret the simulation results.
What is a bad RMSD value? ›
You especially see RMSD as a metric of accuracy when PDB structures are redocked, wherein the same ligand is docked and reference is the file from the PDB (or its PDB Redo derivative), so RMSD greater than 2 Å is bad as it would not match the density, while under 1 is really good —note that in tests about 50% and 75% ...
Is higher or lower RMSD better? ›On average, smaller r.m.s.d. values are associated with protein structure pairs at better resolution and the r.m.s.d. values tend to increase if the two proteins that are compared have been refined at different resolutions.
What does RMSD tell you? ›The rmsd value gives the average deviation between the corresponding atoms of two proteins: the smaller the rmsd, the more similar the two structures.
What is a good TM score? ›TM-score has the value in (0,1], where 1 indicates a perfect match between two structures. Following strict statistics of structures in the PDB, scores below 0.17 correspond to randomly chosen unrelated proteins whereas structures with a score higher than 0.5 assume generally the same fold in SCOP/CATH.
What does a high RMSF value mean? ›Higher RMSF values indicate greater flexibility during the MD simulation.
What is a good molecular docking score? ›Examples of RMSD for docked ligands (gray) with respect to reference ligand at the crystal structures (green) for illustrating good (RMSD ≤ 2.0 Å), acceptable (RMSD > 2.0 Å and <3.0 Å) and bad (RMSD ≥ 3.0 Å) solutions in each target protein.
What is the most accurate chemical analysis? ›pH measurement is the most significant test performed in laboratories as many of the physical, chemical, and biological processes are dependent on pH. The pH measurement using pH meter gives the most accurate results.
What is the most accurate method to determine molecular weight? ›Gel permeation chromatography (GPC) Gel permeation chromatography is also called size exclusion chromatography. It is widely used method to determine high molecular weight distribution.
Which one is best method for determination of molecular weight? ›Osmotic pressure method is best for the determination of molecular weight of polymers and protiens since osmotic pressure though very small, is measurable.
What is the difference between molecular simulation and molecular dynamics? ›In a word, Molecular mechanics (MM) usually refers to running molecular dynamics simulations with a specific force field developed for the related molecules or system. Molecular dynamics (MD) refers to solving of newton's equations, and it doesn't depend on the fore field parameters.
Which software is used for molecular dynamics? ›
Abalone, is a general purpose molecular modeling program focused on the dynamics of biopolymers.
How do I choose a time step size? ›If there is an event that is trying to be captured, the time step size needs to be small enough to capture the event. Example: A chip or heater is on for 30 seconds and off for 10 seconds. The time step size needs to capture the smaller event, so 2-3 seconds would be sufficient.
How long does it take to do a MD simulation? ›To ensure numerical stability, the time steps in an MD simulation must be short, typically only a few femtoseconds (10–15 s) each. Most of the events of biochemical interest—for example, functionally important structural changes in proteins—take place on timescales of nanoseconds, microseconds, or longer.
What is the biggest MD simulation? ›The largest system that contains 1.6 billion atoms was simulated using MD with a performance of 8.30 ns/day on Fugaku supercomputer. It extends the available size and time of MD simulations to answer unresolved questions of biomacromolecules in a living cell.
Are molecular dynamics simulations accurate? ›Force field-based classical molecular dynamics (CMD) is efficient but its potential energy surface (PES) prediction error can be very large. Density functional theory (DFT)-based ab-initio molecular dynamics (AIMD) is accurate but computational cost limits its applications to small systems.
How do I extend my MD simulation in Desmond? ›- Choose Applications → Desmond → Molecular Dynamics.
- Choose Import from file in the Model system section.
- Click Browse and navigate to the checkpoint file ( . cpt ).
- Increase the total simulation time. This is the only setting you should adjust.
- Start the job.
NVT: Simulation of system based on constant number (N), constant-volume (V), and constant-temperature (T);
How does simulation improve patient care? ›Simulation also can be used to streamline protocols or processes without involving patients. Multidisciplinary simulation exercises in the practice setting can be used to identify latent threats and improve readiness.
Which is the best molecular dynamics software? ›VASP, Quantum Espresso, CASTEP, CPMD and ABINIT are the most popular ab initio molecular dynamics software used for calculating and simulating properties of a wide range of materials.
What is the difference between molecular docking and MD simulation? ›Molecular docking is used to study the interaction of two biological molecules (Protein and drug or protein-protein intreactions). MD simulation will give you more clear idea in relation to time that with respect to time what will happen with these interacting molecules.
What is the MD simulation method? ›
Molecular Dynamics (MD) simulations is a computational method that employs Newton's laws to evaluate the motions of water, ions, small molecules, and macromolecules or more complex systems, for example, whole viruses, to reproduce the behavior of the biological environment, including water molecules and lipid membranes ...
What is Amber in MD simulation? ›Assisted Model Building with Energy Refinement (AMBER) is a family of force fields for molecular dynamics of biomolecules originally developed by Peter Kollman's group at the University of California, San Francisco. AMBER is also the name for the molecular dynamics software package that simulates these force fields.
What is the difference between NAMD and GROMACS? ›GROMACS and NAMD each adopts different temperature control schemes: GROMACS mainly uses Nose-Hoover and Berendsen thermostat while NAMD mainly uses Langevin dynamics or temperature coupling scheme. We performd NVT simulations using each of these temperature control schemes to compare their effects.
Which software is used for molecular dynamic simulation? ›LAMMPS is a classical molecular dynamics code with a focus on materials modeling.
What does RMSD mean in molecular dynamics? ›RMSD. RMSD, or root-mean-square deviation, is a standard measure of structural distance between coordinates. It measures the average distance between a group of atoms (e.g. backbone atoms of a protein).
What does RMSD stand for in molecular dynamics simulation? ›However, a very common technique is the root mean square deviation (RMSD). The RMSD is defined as the spatial difference between two static structures: (1) Here, N denotes the number of atoms, i the current atom, rX the target structure, and rY the reference structure.
What is RMSD in md? ›Root mean square deviation (RMSD) is used for measuring the difference between the backbones of a protein from its initial structural conformation to its final position.
What is a good RMSD value for MD simulation? ›An RMSD value ≤ 2 Å is fairly good.
What are the 3 modern simulation methods? ›Simulation systems include discrete event simulation, process simulation and dynamic simulation.
Why do we do equilibration in MD simulation? ›In molecular dynamics (MD) simulations, an equilibration phase is used to bring the system to the desired conditions (e.g. temperature and pressure) before collecting data in more extensive production phases.
How do you analyze MD simulation results? ›
- Do the RMSD first.
- Then do the RMSF calculations.
- Do the PCA. If you want, study the theory of Normal Mode Analysis and do the calculations.
- There comes one of my favourite tools, g_mmpbsa. ...
- H-bond analyses.
- For a protein-ligand system, you can do Umbrella sampling to find out the delG value. ...
- Plot a contact matrix.
Both are very good for the Molecular Dynamics (MD).
What does having amber energy mean? ›Amber is thought to help absorb negative energy and to release bright, soothing energy, helping to calm nerves and enliven disposition like a mental sunny day. The different colors of amber are often used on the chakras with corresponding colors to facilitate opening and cleansing.