Questions tagged [molecular-dynamics]

For questions related to the use of or theory behind classical or ab initio molecular dynamics.

Filter by
Sorted by
Tagged with
18
votes
3answers
554 views

What are the types of ab initio Molecular Dynamics?

I am new to the world of Molecular Dynamics, and am curious to know exactly what is considered to be an ab initio Molecular Dynamics (AIMD) method, and how do they work? The seminal work by Roberto ...
19
votes
5answers
2k views

What software is available to do molecular dynamics on Windows?

Is there a Windows-based molecular dynamics simulation software which is easy to use?
9
votes
1answer
113 views

What are the types of Quantum Molecular Dynamics (QMD)?

In similar spirit to recent questions on Quantum Monte Carlo, ab-initio Molecular Dynamics, types of SCF, and others, I would like to ask: What are the types of Quantum Molecular Dynamics (QMD)? As I ...
14
votes
2answers
239 views

What are the parameters that need to be calculated to design the structure of a molecule?

Before performing Molecular Dynamics on a co-polymer that doesn't have a ready-to-use structure anywhere, the monomer fragment was built in Avogadro software with 0 net charge, B3LYP theory and 6-31G(...
35
votes
2answers
380 views

How should one choose the time step in a molecular dynamics integration?

Choosing too small of a timestep leads to an unrealistic simulation time, whereas too big of a timestep leads to the system not being represented correctly (or, in the case of an algorithm like SHAKE, ...
11
votes
1answer
207 views

What are good resources to learn Materials Modeling? [closed]

I wonder if there are any books or resources to start learning Materials modelling using methods density functional theory quantum monte carlo molecular dynamics
11
votes
1answer
98 views

Selection of appropriate Langevin damping parameter for MD of solid metal

I am performing an MD simulation involving solid metal. The system is in the NVT ensemble, and I am using the Langevin thermostat (fix Langevin in LAMMPS) to enforce this. I have previously used the ...
10
votes
2answers
244 views

What are the factors to consider before docking 2 proteins?

With no background in biological science, I need help with docking 2 proteins together. This issue is with regards to SARS-CoV-2. The 2 proteins are Spike glycoprotein and membrane protein. What are ...
19
votes
2answers
725 views

Does increasing Lennard-Jones cutoff means higher accuracy?

I am currently trying to simulate a large box of QM/MM water molecules with SPC/Fw model and GGA DFT functional (a cube of 40 angstrom sides). To validate the QM and MM simulations separately, I am ...
14
votes
1answer
938 views

How to rectify the error “Lost atoms: original 28437 current 3659” in LAMMPS?

I have modeled a complex structure of a protein and a Carbon nano-cone (used as an indenter to the protein) using PACKMOL. The LAMMPS data file of the complex structure consists of 4 atom types (C O N ...
13
votes
0answers
123 views

Empirically determining thermostat damping factor

Is there a way to determine an appropriate thermostat damping factor given a timestep size and a numerical integration method? I am computing a set of simulations in LAMMPS, using ReaxFF, tracking ...
12
votes
1answer
80 views

Rule-of-thumb for Morse potential cutoff in molecular dynamics?

Is there a rule-of-thumb for selecting a cutoff for the Morse potential in molecular dynamics simulation? For instance, is it typically some multiple of the equilibrium distance?
12
votes
1answer
314 views

Compute RDF between COMs of two selections (or groups) using GROMACS

I have been using gromacs since years but I am always confused when I want to perform the following task: compute the RDF (radial distribution function) between the ...
11
votes
1answer
78 views

Unable to recreate data for benzene from paper in GROMACS 2020.5

I have been trying to simulate benzene in OPLSS-AA forcefield as per this paper: https://pubs.acs.org/doi/10.1021/ct2002122. I have been at a loss, since I believe I have everything running as per the ...
10
votes
1answer
150 views

Can I compute the free energy difference in (temporally) nearby micro-states using the Zwanzig equation for Free Energy Perturbation?

I wanted to understand the Free Energy Perturbation in an NVT ensemble. Wikipedia explains it as: $$ \tag{1} \Delta F = F_1 - F_0 = -\kappa T \ln \left<\exp\left( \frac{-(E_1-E_0)}{\kappa T}\right)...
8
votes
1answer
123 views

What is the “overcoherence problem” in Ehrenfest dynamics?

I am a novice in Ehrenfest dynamics and would like to know why Ehrenfest dynamics suffers from the "overcoherence problem"?