Questions tagged [molecular-dynamics]
Molecular dynamics is a method in computational chemistry that is used to study the movement of atoms and molecules.
80
questions
14
votes
4answers
1k views
Why are MD (Molecular Dynamics) simulations performed at room temperature for the purposes of studying proteins?
All of the MD simulations which are cited in the literature seem to be performed at room temperature (~300 Kelvins), while enzymes are usually bioactive at body temperature (~310 Kelvins).
I am ...
10
votes
2answers
286 views
Why aren't the Anton computers used more widely?
The Anton computers are purpose-built by DE Shaw research for molecular dynamics simulations. They achieve performance up to an order of magnitude faster than can be achieved with commodity hardware, ...
10
votes
2answers
1k views
Boltzmann Distribution in Molecular Dynamics Simulation?
The Boltzmann distribution is often used to describe the population of energy states at different temperatures, i.e. something like
$$\dfrac{g_i}{g_j} = \mathrm{exp}\left(-(E_i-E_j)\times \frac{1}{...
10
votes
0answers
121 views
In a molecular dynamics context, is the methyl rotation in propene a symmetric or asymmetric internal rotor?
In a prior question, I asked about the (a)symmetry of the potential energy surface of the methyl rotation of propene. In that context, the kinetic energy of the nuclear motions is of course assumed ...
9
votes
4answers
455 views
Is it true that an evaporating molecule has the same kinetic energy as a molecule in a pot of boiling water?
A molecule on the surface of room-temperature water shoots off the surface of said water, or in other words, it "evaporates." It does so because it gained kinetic energy ${x}$, and ${x}$ was great ...
9
votes
1answer
1k views
How to perform molecular dynamics simulations of charged systems?
The GROMACS official documentation (see here) states that a system with non-zero total charge will yield an error:
System has non-zero total charge
Notifies you that counter-ions may be required for ...
8
votes
1answer
719 views
How to simulate phase separation in a discrete system using LAMMPS?
I would like to simulate phase separation in LAMMPS and need advice what kind of two particle systems and interactions will yield phase separation.
For example, a meta stable phase separation is ...
8
votes
1answer
106 views
How did Noe et al. derive $k_{AB}$ in their PNAS 2009 paper?
This question is regarding Constructing the equilibrium ensemble of folding pathways from short off-equilibrium simulations which uses Markov State Model to find the equilibrium protein folding ...
7
votes
0answers
347 views
Cloud-based molecular modeling software
I’m an IT Acquisitions and Licensing specialist researching well-supported cloud software alternatives to Fujitsu’s molecular modeling SCiGRESS. I am coming in with a tech background, but little to ...
6
votes
1answer
775 views
Why do the peaks in a radial distribution function graph have a probability >1
Why do RDF plots have the probability at a particular distance >1. How can probability be more than one?
I am asking this because I want to understand what do the peaks of different length but at ...
6
votes
1answer
335 views
How to write an entire trajectory into a single PDB file?
I have a trajectory (few frames) that I want to stack together and build a single PDB. I thought either to use VMD or MDAnalysis library.
With VMD if I do ...
6
votes
2answers
601 views
Force Constants for Angles From Gaussian To AMBER: which value should I put in the matrix equation?
I am trying to make a new force field in AMBER for a novel compound including boron (which does not exist in GAFF). For this, I calculated freq and I found bond force constants but I couldn't find ...
6
votes
1answer
77 views
What types of energies can be determined by molecular dynamics in gas phase or explicit solvents?
For molecular dynamics in the gas phase I've seen that electronic energy can be determined and also decomposed in kinetic and potential energy. I've also seen that we can include ZPE. But can we ...
6
votes
0answers
78 views
Solving the Schrödinger equation for a rotating triatomic linear molecule
This source is showing that solving the Schrödinger equation for a triatmoic linear molecule yields the same formula for the rotationaI quantum states $BJ(J+1)$ as for dipoles.
For dipoles, the total ...
6
votes
0answers
96 views
Why the use of MD simulations is less popular for the study of carbohydrates?
In classical all-atom MD simulations, the force fields for proteins, nucleic acids and lipids seem to be more well developed than carbohydrates. It seems to me that the use of MD simulations is not so ...
5
votes
1answer
52 views
Vibrational state dependency of the reaction cross section
I'm reading p.1238 of Mcquarrie's Physical Chemistry: A Molecular Approach textbook.
In the last line of the page, the author explains that the value of the reaction cross section in the constant ...
5
votes
1answer
111 views
Building an average model from MD snapshots
I am performing a molecular dynamics simulation, and I am taking "snapshots" of the system at regular intervals.
The aim is that of finding the preferred geometry of a molecule in solution.
Are you ...
4
votes
1answer
108 views
Tools for Molecular Simulation
I am a high school student and I don't know a whole lot about molecular simulation. I wish to view different conformers of cyclohexane and the changes that occur when substituents are attached (for ...
4
votes
1answer
263 views
Hessian Matrix and Physical directions: Potential Energy Surface
How do the eigenvectors of the Hessian Matrix for potential energy surface correspond directly to the physical directions of translation, rotation and vibration?
I see that they might correspond to ...
4
votes
1answer
55 views
Would there be need for periodic DFT if we are able to get all the forcefields of all the elements in the periodic table in ReaxFF?
Since ReaxFF combines the advantages of DFT and Classical MD, would there be any need for periodic DFT calculations if we can have all the forcefields of all the elements of all the possible systems?
4
votes
1answer
59 views
Keeping the electronic configuration in the ground state is expensive in Ab-initio MD?
In E. Sandré & A. Pasturel's An Introduction to Ab-Initio Molecular Dynamics Schemes
If using small time steps (typically in the order of a few tenth of femto-second), the changes in the ...
4
votes
1answer
205 views
What is the “curse of dimensionality” in molecular dynamics?
The 'curse of dimensionality' or the 'bottleneck problem' in molecular dynamics is explained in page 5 and 6 of Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods by Dominic Marx and Jurg ...
4
votes
0answers
141 views
What computational methods are available for predicting nanoparticle structure?
I'm a quantum chemist, and nearly every system I've ever cared about has been studied by DFT. However, I am currently trying to model the surface interactions on a nanoparticle around 100 nm in ...
4
votes
0answers
185 views
LAMMPS: MD simulation of molecules? [closed]
I need to do molecular dynamics simulation of cyclodextrine molecules for a project and I am searching for an open source solution that can be used. I found LAMMPS to be quite well maintained and ...
3
votes
1answer
63 views
How much time for a water molecule in pure water to reach the top of a glass of water?
I have a curiosity to know how long it will take for a water molecule to move a certain distance in pure water. I can put it simply as how much time for a water molecule to travel the whole length of ...
3
votes
1answer
134 views
Is there a known relation between the diffusion of a molecule and the diffusion of the reactants?
I have a plasma consisting of elements A and B (which have a similar mass). These elements can react to form molecules AB.
$$
\ce{A + B -> AB}
$$
Now if my plasma or gas is expanding based on ...
3
votes
1answer
51 views
Determining the torsion energy profile of butadiene by molecular dynamics
Disclaimer: this question is related to my previous question, although it is much more detailed
Consider the following rotation of butadiene:
I want to determine the relative energy of each dihedral ...
3
votes
1answer
355 views
How can you simulate the effect of a solvent on a polymer in an MD simulation?
I'm doing a molecular dynamics simulation of polymer melts and I want to include the effect of a good/bad solvent on the melt. How can this be done?
My simulation models polymers using a bead-spring ...
3
votes
0answers
19 views
Can someone help me understand the motivation behind and visualize the dihedral potential and the Urey-Bradley potential?
I am studying statistical mechanics and force-fields, and I see a lot of this formula being thrown around with no explanation:
$$U_{\mathrm{dihedral}} = \begin{cases}
k(1+\cos (n\psi + \phi)),&...
3
votes
0answers
28 views
Normal Mode Analysis of proteins
I'm trying to do Normal Mode Analysis for my MD trajectory. I'm following Amber tutorial (http://ambermd.org/tutorials/advanced/tutorial3/py_script/section5.htm) which uses an example of estrogen ...
3
votes
0answers
121 views
Calculate collision density of air
Caluculate the collision density between $\ce{O2}$ and $\ce{N2}$ in air at $T = \pu{10 ^\circ C}$. Assume that the air is an ideal gas and that the molecules are spheres with diameters $\pu{0.36 nm}$ ...
3
votes
1answer
51 views
How to extract force when the center of mass of molecules is put in the box in GROMACS?
To extract the force as in the title I used the gmx trjconv command like below.
...
3
votes
0answers
82 views
Simple sandbox software for an interactive molecular dynamics simulation
I'm not sure if this is the right site to ask a such question. But I'm getting to learn the world of the small.
I know there is software out there like NAMD, and LAMMPS. But these tools are very ...
3
votes
0answers
158 views
Minimum image convention and the atom “seeing” itself
The minimum image convention in molecular dynamics can be stated as:
A particle in the central square must not interact with one of its neighbors in the central square and an image of this neighbor ...
2
votes
1answer
136 views
Calculate force on atoms given the energy of configuration, trajectory of simulation and position of atoms in each frame
Background
I am using a neural network to calculate the potential energy of atoms in a configuration and then adding energy of all atoms to compare it with the true energy of the configuration(label) ...
2
votes
1answer
85 views
Calculation of free energy from potential and kinetic energy
I am having trouble understanding how (and even whether it's possible) to get free energy from kinetic and potential energy.
I have run molecular dynamics simulations (NPT ensemble, CP2K software), ...
2
votes
2answers
976 views
Question during radius of gyration calculation from XYZ trajectoris of atoms
I have molecular dynamic simulation result of 100 polymer molecules (liquid), and I'm trying to calculate the radius of gyration using atomistic mass, molecular mass, and xyz coordinate of atoms. I'm ...
2
votes
1answer
69 views
Understanding water models?
I generally work on fluid flow, so some basic concepts in chemistry are new to me.
Work I'm doing requires use of molecular dynamics (MD) to do classical potential simulations of water–NaCl systems ...
2
votes
1answer
27 views
How can I build the wyoming like montmorillonite unit cell in avogadro for a molecular dynamics simulation?
I am trying to build a unit cell for Wyoming-like montmorillonite, using Avogadro software. I am not sure how can I get access to the atom coordination and how to generate the .xyz file at the end? Is ...
2
votes
1answer
595 views
Radial pair distribution function (VMD). How to define it for a water NaCl system with multiple Na and Cl?
I'm struggling to understand a paper in which someone has presented radial distribution functions (RDFs) for a system that contains multiple Na, Cl, and water atoms/molecules. I don't know exactly ...
2
votes
1answer
230 views
Courses for Molecular Modeling? [closed]
I'm probably doing overkill trying to learn a bit of everything, but here's my story. I have a B.S. in Math with a minor in Microbiology, and I'm currently taking Computer Science courses online to ...
2
votes
0answers
33 views
Molecular dynamics simulation of a protein in acidic medium
I want to perform an MD simulation of a protein under acidic solvent conditions. A quick literature search seems to indicate that people are more interested in the protonation of protein side chains (...
2
votes
0answers
37 views
How are class II force field cross-terms derived?
In some class II force fields, there are cross-terms that characterize bond-bond, bond-angle interactions and others. Example of equations below:
$$\sum_{bond-bond}^{}k_{bb}(b-b_{0})(b^{'}-b^{'}_{0})$$...
2
votes
0answers
59 views
Relationship between reaction quotient and potential energy [closed]
I am trying to calculate the Gibbs free energy in an NVT simulation. Although
$$\Delta G = -RT\ln Q,$$
in my simulation I can't calculate
$$Q = \frac{[\text{Products}]}{[\text{Reactants}]}$$
because I ...
2
votes
0answers
30 views
Difference between docking and scoring?
I'm taking an undergraduate level class on computer based drug design. The class is more for biology students than people with a strong understanding in physics and computer science.
With that said,...
2
votes
0answers
51 views
Do I understand correctly that simulation volume influences relative permittivity calculations? If so, how is the result useful?
I am having a lot of trouble understanding the influence of volume on relative permittivity/static dielectric constant. I'm not a chemist, but am using molecular dynamics to calculate relative ...
2
votes
0answers
62 views
How to Scan Aromatic Dihedrals To Obtain Force Field Parameters
I am trying to generate new force field parameters based on quantum mechanical calculation. I used amber manual "parameter development" title for this and I calculated bond-angle-non bond parameters ...
2
votes
0answers
193 views
Amber Force Field Improper Parameterization Huge Kb Problem
I am trying to parameterize a molecule. I know dihedrals and impropers are the same for Amber Force Field. (Source: http://archive.ambermd.org/201004/0334.html)
So, I did PET scan for improper and I ...
2
votes
0answers
177 views
Amber paramfit “monomer (or ligand)” problem [closed]
I am trying to parameterize a novel molecule which has atoms not included gaff or parm force field. I finished parameterizing vdw-bond and angle force field. Also I did PES Scan for dihedrals and I ...
2
votes
0answers
221 views
How to develop a Reactive Force Field (ReaxFF)? (Molecular simulation)
Currently, I am working on CO2 capture on Hydrotalcites as solid adsorbents. I want to run molecular simulations to analyse the behaviour of the material under different temperatures, pressures and ...
