Questions tagged [molecular-dynamics]

Molecular dynamics is a method in computational chemistry that is used to study the movement of atoms and molecules.

Filter by
Sorted by
Tagged with
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 ...