Questions tagged [molecular-dynamics]
For questions related to the use of or theory behind classical or ab initio molecular dynamics.
401
questions
6
votes
2
answers
2k
views
How to calculate diffusion coefficient from MSD graph using GROMACS?
In GROMACS I have used the command below to calculate mean square displacement (MSD) over time,
...
- 305
1
vote
1
answer
445
views
I have a force in eV/Angstroms. How do I convert it so that the speed I calculate is in Angstroms/second?
How can I express speed in Angstroms/second after application of a force that's in eV/Angstrom (given the mass, original velocity and time elapsed)?
- 564
6
votes
1
answer
498
views
Is it possible to build a force field that suits all elements based on VASP's machine learning result?
I have tried to use VASP's machine learning force field calculation during running molecular dynamics simulation with a supercell including some elements of Ti, O, Cu, etc. It does increase the speed ...
- 1,977
3
votes
0
answers
83
views
What are the similarities between MD and MC in the case of protein simulation? [closed]
As far as I understand, both MD and MC simulation methods use Force Fields.
Are there any other similarities between them?
- 2,352
4
votes
0
answers
98
views
Radius of gyration and aggregates [closed]
I have performed a MD simulation using GROMACS for a system of 100 caffeine molecules in water. It is expected caffeine form aggregates.
I have calculated the radius of gyration for the simulation ...
- 41
3
votes
0
answers
396
views
Calculating the Diffusion coffecients from Trajectory file in ASE [closed]
I plan to study the Diffusion-Coefficient of my simulation outputs, for this, I want to use ASE.
Sadly I am not entirely sure of how to proceed, has anyone done this and if yes, please share your ...
- 949
8
votes
1
answer
375
views
How to run 2nd generation CPMD in CP2K?
I believe CP2K by default implements Born–Oppenheimer MD. I have seen many recent papers that use 2nd generation Car–Parrinello MD using CP2K. But I couldn't find any direct way of switching to 2nd-...
- 1,581
3
votes
1
answer
392
views
Domain decomposition error while applying bonded restraints - GROMACS
GROMACS version:version 2021.3-MODIFIED
GROMACS modification: Yes
Hi, I have a small system with 87 residues of a protein. And I want to maintain the tertiary structure of the protein, so I came up ...
- 1,291
5
votes
0
answers
179
views
Regarding use of GCMC to add water molecules
I am working on a goethite system, and wish to fill the unit cell with water molecules for AIMD simulation using quantum espresso. For this, some papers[1],[2] indicated the use of Grand Canonical ...
- 911
3
votes
1
answer
186
views
Is the energy due to external electric field included in gromacs energy?
I have a Gromacs NVT simulation with a constant external electric field using the option "electric-field-z = 1.5 0 0 0".
When calling gmx energy -f (...).edr to calculate the potential ...
- 131
7
votes
1
answer
265
views
Simulating autoionization of small molecules using DFT or AIMD
From what I'm reading, simulation of auto-ionization in liquid substances is hard to do using gas phase DFT or AIMD [1][2][3]. I am using ORCA to do my work and I must confess that I'm shocked to see ...
- 327
3
votes
0
answers
41
views
How may I define the pressures of two sides of a membrane in Materials Studio? [closed]
To study SiO2, polyamide and two-layer SiO2-polyamide as the membranes in the reverse osmosis desalination of saline water; I used the Materials Studio software to model the systems. My question is:
...
4
votes
1
answer
1k
views
Why does the temperature fluctuate severely during the AIMD using an NVT ensemble?
I am trying to raise the temperature of my model and keep the temperature at 300K, with an NVT ensemble and Nose-Hoover thermostat.
I am doing this job with VASP, here is the ...
- 1,977
3
votes
0
answers
416
views
Atomic Simulation Environment (ASE) is unable to read a big trajectory file [closed]
The problem that I face is the ASE software can't read the trajectory file and throws this error ValueError: buffer is smaller than requested size. My assumption is ...
- 949
4
votes
1
answer
114
views
How can we test the force field (parameter) for a complex?
I built a force field for a complex. How can I check if the force field is working well or no?
I did a short molecular dynamics using AMBER only for the complex using the parameters I built. It worked ...
- 1,285
6
votes
2
answers
718
views
Parameterising a Lennard-Jones interaction: Which atoms should I scan to build potential energy surface?
I have a question related to parameterize Lennard-Jones interaction to get sigma $\sigma$ and epsilon $\epsilon$ to include in amber force field to do the Molecular Dynamics.
I searched a lot, they ...
- 1,285
4
votes
0
answers
84
views
Periodic boundary conditions or spherical potential? [closed]
I was wondering whether you know if there is any sort of literature discussing whether periodic boundary conditions or spherical potentials could be better or worse for QM/MM or MM simulations, their ...
- 967
5
votes
1
answer
107
views
Where to find/ to get Tersoff and Stillinger Weber Potentials for MD simulation using LAMMPS
I'm trying to calculate the thermal properties of Silicene with hydrogen and halogen atoms. My problem is that, there's seems no available potential for Tersoff and Stillinger Weber Potential. and I'...
- 129
3
votes
1
answer
294
views
Calculate the atomic charge for transition metal using RESP?
I calculate the atomic charge using restricted electrostatic potential (RESP) and compute with the Hartree Fock method. This method is OK for organic compounds, but for complexes that contain ...
- 1,285
8
votes
0
answers
179
views
Young's modulus from MD simulation [closed]
I have a coarse-grained model of a surface which is almost planar/flat in 2D. Now, I want to compute some of its mechanical properties to test it. Young's modulus is one of the properties that come to ...
- 1,031
10
votes
1
answer
364
views
Contact Angle (water droplet on a surface)
I have built a CG model of a surface on which I want to place a water droplet (some water molecules), and calculate the contact angle. Has anyone done similar work with gromacs (MD simulation)? First, ...
- 1,031
4
votes
0
answers
58
views
Gromacs generate forces using .pdb topol.top and .mdp file [closed]
Given the start configuration (.pdb), topology file (topol.top) and run file(water.mdp).
I want to generate the forces on atoms. I tried, using:
...
- 1,083
7
votes
0
answers
248
views
What does the shape of the radial distribution function(RDF) indicate? [closed]
From my understanding, generally in the RDF plot for liquids, the value of $g(r)$ goes to a maximum near the first solvation shell and then goes to a minimum and then slowly tapers off to the value of ...
- 1,407
3
votes
2
answers
243
views
How to quantify the lattice distortion from XRD?
In real crystals, the atoms may not always sit at their ideal lattice sites. How to quantify this lattice distortion?
Based on my understanding from X-ray diffraction books, the Radial Distribution ...
- 633
7
votes
1
answer
202
views
Is it incorrect to use a NVT ensemble for Non-equilibrium MD simulation?
I am getting some conflicting results on this topic while searching for literature online, so I am asking this question here with respect to a simplified experimental setup. If we are trying to ...
- 2,262
7
votes
0
answers
171
views
Diffusion coefficient as a function of distance from the position autocorrelation function [closed]
I am trying to understand if my idea for the computation of diffusion coefficient as a function of distance is right or not. So, the diffusion coefficient as a function of distance is given by,
$$D(z) ...
- 1,031
6
votes
1
answer
165
views
Understanding derivation of discretized FIRE algorithm
Introduction/Preamble
@SusiLehtola's answer to Basics of numerical energy minimization techniques used in molecular dynamics? mentions conjugate gradients, BFGS for energy minimization, Metropolis ...
- 1,874
6
votes
0
answers
57
views
Basics of configuration initialization strategies used in molecular dynamics? [duplicate]
This answer to Basics of numerical energy minimization techniques used in molecular dynamics? mentions
conjugate gradients, BFGS for energy minimization
and
Metropolis Monte Carlo (and the) the ...
- 1,874
10
votes
1
answer
295
views
Basics of numerical energy minimization techniques used in molecular dynamics?
The question below describes my plan to make a basic molecular dynamics calculation using a Python script rather than a canned, self-contained program.
There seems to be three parts:
a model of the ...
- 1,874
8
votes
0
answers
99
views
DIY molecular dynamics for Xenes on crystal surfaces; where can I get applicable open-source force field parameters that I can use in my scripts?
The linked question(s) below describes a DIY project to simulate how Xenes (honeycomb nets of atoms like graphehe) behave on single crystal surfaces. Behaviors include rotation, strain, and heigh ...
- 1,874
17
votes
4
answers
4k
views
I am considering to purchase a GPU. What calculations would need FP32, and what would need FP64?
I would like to know how to choose the best GPU on the basis of its FP32 or FP64 efficiency, depending on different types of simulations.
In brief I would like to know what types of simulations work ...
- 1,479
6
votes
1
answer
107
views
How do biaxial potentials like the Gay-Berne potential work?
I am trying to understand and implement the Gay-Berne potential, but I do not quite understand how one goes about constructing them.
I am following this paper to parse through the idea.
I understand ...
- 4,141
7
votes
1
answer
529
views
6
votes
0
answers
116
views
How to calculate and plot the electrostatic forces in a .pdb file [closed]
Given a pdb file that looks like the following:
...
- 1,083
16
votes
1
answer
320
views
What are the advantages of (semi)-empirical force fields over Machine Learning Potentials?
I am currently working with ReaxFF, an empirical reactive force field that can describe chemical bond forming and breaking. The main advantage over ab initio methods are of course the greatly ...
- 263
15
votes
3
answers
2k
views
Can DFT simulate chemical reactions?
I was wondering, can you simulate the chemical reaction between two or more molecules using DFT? If not, what numerical method is usually used for this?
8
votes
1
answer
218
views
Parameterize a complex that contains Transition metal (Build Force Field)?
It is a general question. I know, but I need some help.
I am interested to parametrize a complex that contains transition metal like Ru and Pt etc.
Furthermore, I will use this Force Field to do ...
- 1,285
8
votes
1
answer
118
views
What software should I use for an MD study of a temperature-dependent crystallographic phase transition?
I want to study the crystallographic phase transitions of crystals with increasing temperature. Which freely available MD package should I choose, which has easily available tutorials/example-input-...
- 181
7
votes
1
answer
2k
views
How to generate a .tpr file from a .pdb or .xyz file in gromacs
I have a .pdb file for a single frame of water in a box. I need to calculate and analyse the electrostatic field for choosing a viable atom in this particular ...
- 1,083
6
votes
1
answer
56
views
Is it possible to simulate the intercalation of hydrogen nanobubbles?
Is it possible to simulate the intercalation of hydrogen nanobubbles into 2D layered materials using first principles calculations such as Density Functional Theory? Would this require a multi-scale ...
- 119
7
votes
0
answers
70
views
Can Markov State Models based on short-term non-equilibrium datasets also capture the dynamics towards equilibrium? [closed]
I am non-specialist who is new to the field of computational physics/biology. I have a few questions about the relevance and effectiveness of Markov State Models (MSM) in capturing the dynamics ...
- 171
13
votes
0
answers
445
views
Applicability of GFN-FF force field and other xTB methods for dynamics of ions in water [closed]
I am working on a project involving water clusters with ions in them. One of the problems that I ran into quite early on was that good forcefield parameters are not available for highly charged ions ...
- 6,516
10
votes
1
answer
357
views
Very high negative pressure in MD simulation of liquid quenching
I am trying to prepare turbostratic carbon structure from an initial random carbon sample with a density of 1.8gcm$^{-3}$. The initial sample image is attached here:
The problem is the high negative ...
- 103
10
votes
1
answer
133
views
What does the slope in the viscosity vs shear acceleration graph from non-equilibrium MD signify?
In finding out the viscosity of a system with non-equilibrium MD using GROMACS, an acceleration is given to the system due to which, the velocity profile changes. This change in energy is used to ...
- 1,407
10
votes
1
answer
155
views
Generalized Born implicit solvent models and pairwise descreening
It seems almost all the common implementations of Generalized Born (GB) implicit solvent models (in Amber, OpenMM etc) use the "pairwise descreening" approximation, which is described in ...
- 801
5
votes
0
answers
86
views
Define stress-strain for a multi-wall carbon nanotubes bundle [closed]
I want to simulate a MWCNT bundle in which the central MWCNT is pulled out of the bundle.
How to define the stress/strain curve of such system?
I am building my system with scikit-nano and using ...
- 667
5
votes
0
answers
123
views
Is there a way to set different special bonds for different atom groups in the same structure in LAMMPS? [closed]
I am running some simulations on a metal organic framework and I need to set different 1-4 coul and lj parameters for different ...
- 51
6
votes
1
answer
108
views
How to check for phase separation in polymer simulations?
Suppose I am running a molecular dynamics simulation of a polymer solution that exhibits upper critical solution temperature. From standard solution thermodynamics, we expect to see a polymer-rich ...
- 4,141
8
votes
0
answers
254
views
Contribution in action from collisions for a gas?
So if I write the action for a gas:
$$ S = \int (T -U) dt \tag{1}$$
where $T$ is the kinetic energy and $U$ is potential energy.
I suspect there is a constant term (which does not affect the equations ...
- 415
6
votes
0
answers
477
views
LAMMPS fix print function [closed]
First of all, I'm very new to LAMMPS so I'm sorry if my question is trivial.
I'm trying to model a gas of atoms that float in a square box, they obey an NVE integration and only feel a shifted Lennard-...
- 61