Questions tagged [molecular-dynamics]
For questions related to the use of or theory behind classical or ab initio molecular dynamics.
444
questions
4
votes
0
answers
32
views
MD simulations of liquid isopropanol predicts wrong thermodynamics proprieties
I am using GROMACS to run molecular dynamics simulations of pure liquids, in the NpT ensemble at room temperature and p = 1 atm. My first test was with isopropanol (also known as 2-propanol). I used ...
1
vote
0
answers
17
views
Thermal Conductance along the lateral and longitudinal direction
Suppose the left and right leads are both copper rod and the central scattering region is one monolayer graphene sheet. There is no chemical bonding formed between Cu and C atoms. The copper and ...
4
votes
1
answer
70
views
Thermal vs dynamical stability
Dynamical stability is checked by a phonon curve. Thermal stability is checked by a molecular dynamics simulation. What is the difference between them?
5
votes
2
answers
335
views
Calculation of pH in Molecular Dynamics simulation?
I have a system that has 5M aluminum hydroxide (10 molecules) ion [Al(OH)4]- and 5M NaOH (10 molecules) in a box of side length 15 angstroms each. I have 113 water molecules in the system, so that the ...
7
votes
0
answers
25
views
Question about Constructing Supercell Layers for MD Simulations for Solid-Solid Phase Transitions
I am a newcomer to MD simulations,
so my questions may be at a basic level,
but I appreciate your understanding.
I would like to simulate solid-solid phase transitions using MD simulations.
To begin ...
5
votes
1
answer
59
views
How can I specify the direction of motion in the colvars file?
I am doing umbrella sampling in NAMD and I want to control the distance in the colvars to be in one direction. Is there any possibility?
2
votes
1
answer
27
views
How can a harmonic constraint be added between two groups in one direction using NAMD?
As a new user of NAMD, I am seeking assistance regarding a specific task. I intend to pull a LIGAND from outside to the center of a lipid using NAMD. To achieve this, I have utilized the colvar module ...
4
votes
1
answer
66
views
How can one create various arrangements of two or more types of atoms (like O, N etc.) with distinct coverages on a metal surface?
How can I modify the procedure described in this answer to achieve adsorption of two or more types of atoms on a surface, each with distinct coverages?
Is there an extension to the described method ...
4
votes
0
answers
35
views
Modelling excited states using classical MD simulations
My question is related to the question here but not exactly the same. The studies which have looked into the excited states of molecules have generally done it by parameterizing the ground-state and ...
2
votes
1
answer
44
views
How to write the REBO force field
I want to use the REBO-CHO force field to study graphene oxide using lammps. But I don't know how to create the input file with the parameters by reading the article or if there is a tutorial or ...
4
votes
1
answer
97
views
Advantages/Disadvantages of Qeq over Drude polarizable model?
For classical Molecular Dynamics (cMD) simulations, the usual methodology is to use force fields (FFs) with fixed point charges (examples include CHARMM, AMBER and Gromos). However, this approach ...
4
votes
2
answers
83
views
Given a .gro file of a monomer, how to polymerize it?
I have a structure (.gro) file of a PolyEthylene Glycol (PEG) molecule as shown below. How do I get the corresponding structure file for PEG-20 i.e., a molecule with a degree of polymerization=20, ...
4
votes
0
answers
71
views
Solvation_analysis subpackage for MDAnalysis multi-atom solute
I have a simulation trajectory from DFTB+ for 1 NaNO3 in 32 H2O molecules. The trajectory is in xyz format and I converted a frame of xyz to pdb using openbabel to serve as topology input.
I am ...
3
votes
1
answer
86
views
Controlling Initial Velocities in NVE Trajectory Simulations with VASP
"I've determined the initial velocities of the surface following its equilibration at 300 K in the NVT ensemble. Now, I aim to conduct AIMD (NVE ensemble)with the metal surface (initially defined ...
2
votes
1
answer
63
views
Velocity autocorrelation function for molecular dynamics trajectory
I have a 1 ns long molecular dynamics trajectory of 884 molecules of water, and I am trying to compute the velocity autocorrelation function using MDAnalysis to analyze the GROMACS output files.
I ...
2
votes
0
answers
42
views
Prerelaxation before AIMD with VASP
I am investigating a bulk structure by ab initio molecular dynamics (AIMD) implemented in VASP.
I was wondering, how to choose the initial geometry (= POSCAR file) to run the AIMD simulation. In the ...
1
vote
1
answer
54
views
How to compute the Elastic Modulus of a molecule using GROMACS?
Say we take the basic example of a PEG molecule (PolyEthylene Glycol). What is the fundamental physics for computing the Elastic Modulus or the Young's Modulus from the Molecular Dynamics (MD) ...
2
votes
2
answers
54
views
My CHARMM IC table generation has caused hydrogen to bond to two carbons. How do I fix this?
As above, I am attempting to generate an IC table for decanoic acid using CHARMM (version 48b1) so I have a minimized geometry. I have attached my build.inp script and the stream file I have written. ...
1
vote
1
answer
73
views
What is the best way to optimize molecules with 8000 atoms. This includes C,H,O,N,Si
My molecule includes C,H,O,N,Si
2
votes
1
answer
19
views
Geometry of H2O in the v=1 vibrational state for AIMD simulation
This is a follow-up to my recent question: Algorithm for finding the v=1 state of an H2O molecule
I want the geometry (atomic positions) of the H2O molecule in the v=1 state. How are initial ...
1
vote
0
answers
28
views
Error with running GROMACS with SLURM
I am running into errors while submitting a job for a basic MD Simulation with HPC on SLURM. My sbatch script reads as follows:
...
1
vote
1
answer
56
views
Equilibrating Surface Temperature Through NVT Ensemble Ab Initio Molecular Dynamics
In my AIMD simulations with the NVT ensemble at T= 300K (Nosé–Hoover thermostat)for a Pd (111) surface (2*2, 6 layers) and vacuum sizes of 14 Å and 18 Å, All other INCAR tags are same(only NGZ for 14Å ...
2
votes
0
answers
23
views
Eliminate the entropic contribution in the free energy of dimerization between two molecules in umbrella sampling
For my set-up to compute the free energy of dimerization between two molecules in a solvent using umbrella sampling, I can see that there is an entropic contribution since I have used distance as a ...
6
votes
1
answer
175
views
Speed of atom from kinetic energy of molecule
What is the velocity of each hydrogen atom in an H2 molecule, given a kinetic energy of hydrogen molecule is 0.12 eV, when expressed in angstroms per femtosecond and assuming equal velocities for both ...
1
vote
1
answer
69
views
How to use martinize.py?
I have been trying to get the CG representation of a basic octane molecule. The pdb file is as follows:
...
3
votes
1
answer
99
views
How to convert a SMILES into a LAMMPS or GROMACS input file?
Are there any dedicated Python tools/packages/scripts that help me convert a SMILES structure into an MD coordinate and topology file (both for GROMACS and just the single one for LAMMPS)?
Demonstrate ...
4
votes
1
answer
38
views
phase average uncertainty in molecular dynamics
I was wondering if I could find some scaling relation of the uncertainity of phase averages of thermodynamic quantities with the number of atoms $N_a$, the number of timesteps $N_t$ and size of the ...
2
votes
0
answers
31
views
Compute the monomeric fraction in MD simulation
I have a system of molecules in two different solvents in which aggregation occurs, but visually in different ways. In one, aggregation occurs rapidly, but in another sparsely. These MD simulations ...
2
votes
1
answer
139
views
Choice of SMASS for AIMD simulations in VASP
How can one verify that the values (0, 1, 2, and 3) in the SMASS parameter correspond appropriately to the NVT ensemble(Nose-Hoover thermostat) in the context of an ab initio molecular dynamics ...
3
votes
1
answer
138
views
having trouble to get the NPT working
I am currently encountering challenges in successfully conducting NPT simulations in LAMMPS for a triclinic cell persistently remain unchanged throughout the NPT simulation.my input file is:
...
6
votes
0
answers
103
views
Why is potential energy used and not total energy when calculating free energy?
I have run a replica-exchange molecular dynamics simulation, using a number, $n$, of replicas initiated at the same number, $n$, of unique temperatures. I have calculated the free energy profile over ...
2
votes
0
answers
31
views
Help with understanding the method of MBD & MBD-NL
I am having a hard time understanding the steps of Tkatchenko's method of many-body dispersion (MBD) and nonlocal many-body dispersion (MBD-NL).
I would greatly appreciate a detailed & clear ...
2
votes
1
answer
45
views
How to make initial trajectory for incident A2 molecule on any surface?
I want to write code to simulate the impingement of any molecule on a surface at different (r, theta, and phi) angles in XYZ format.
4
votes
0
answers
99
views
How to perform a rotation of a structure, specifically throughout molecular dynamics or umbrella sampling simulations?
I have a structure, as depicted in the figure, which is a part of a longer chain. The portion within the red circle is predominantly parallel. My objective is to rotate this part, aiming to orient the ...
6
votes
0
answers
120
views
Free energy calculation as a function of temperature over phase transition
I've run a molecular dynamics simulation (NVT) of a simple Lennard-Jones system of a few particles and I'm trying to calculate the Helmholtz free energy $\Delta F$ as a function of temperature over ...
4
votes
1
answer
35
views
Is water penetration and fatty acid dissolution reasonable when Gromacs modelling of an aqueous lipid bilayer?
as mentioned, I'm building a molecular dynamics simulation of decanoic acid and decanoate (with neutralising sodium cations) in water. I've minimized and equilibrated the simulation up to 295 K and I'...
4
votes
0
answers
176
views
Free energy as a function of temperature over phase transition in molecular dynamics simulation
I've constructed a low density standard MD simulation (in Python), using the Lennard-Jones potential with an Andersen thermostat, of 10 particles in a box V = 10 $\times$ 10 $ \times$ 10 in the NVT ...
4
votes
1
answer
107
views
How do we calculate diffusion coefficient when GROMACS MSD graph with two slopes?
I calculated the MSD using the gmx msd -f run05.part0005.xtc -s run05.tpr -o msd_200-250ns.xvg. As a result, I got a graph with 2 slopes as attached.
In the MD simulation, I followed a simulated ...
6
votes
1
answer
97
views
Lattice constant of an FCC crystal from molecular dynamics simulation
So, I have a system of n number of molecules inside a box which is supposed to be solid. I performed the md simulation of the system using stochastic integrator in a box. I know the lattice of the ...
4
votes
1
answer
65
views
DFPT vs analytical nuclear gradient
I have read some literatures about the gradient wrt. atom positions of the DFT energy, and find different methods one commonly used in MD(analytical nuclear gradient,pulay's paper), and another ...
4
votes
1
answer
201
views
Random velocity generation
How to generate random velocities in the geometry coordinates of a surface (POSCAR file) using the Boltzmann distribution law with Atomic Simulation Environment (ASE) or other codes for Molecular ...
1
vote
0
answers
33
views
LBFGS -- Incomplete Convergence due to BadIntpln [closed]
I am doing energy minimization of a molecular system in TINKER. While doing so, I am getting the error in the title of this question. What might be my problem?
2
votes
1
answer
61
views
Is there a way to force GROMACS to keep a certain density for water?
I am trying to fully populate my GROMACS simulation, but the water density in the bulk is too low.
Here are the steps I used:
Add organic monolayer into 9x9x9 box
Solvate with water
equilibrate t-...
3
votes
1
answer
54
views
Exclusion of Lennard-Jones potential between O-H and H-H in water models
In reading Jorgensen's 1983 paper comparing the properties of 6 water models it was mentioned that there is the inclusion of an L-J potential between pairs of oxygen molecules, but no mention of ...
2
votes
1
answer
150
views
Absorption enthaply at different temprature in VASP software via ab initio molecular dynamics
Can we calculate the absorption enthalpy at different temperature in VASP software via ab initio molecular dynamics? If yes, how?
1
vote
0
answers
96
views
Free energy calculations for lammps
Is there any tutorial or reference available to perform free energy calculations (umbrella sampling; thermodynamic integration; free energy perturbation; wisdom's particle insertion method) etc. with ...
3
votes
0
answers
32
views
How to perform atomistic spin dynamics (ASD) Ab Initio Molecular Dynamics (AIMD) using VASP?
I have tried normal AIMD simulation using VASP for Non magnetic materials. Now, I am interested to perform the molecular dynamics simulations for the finite temperature for Magnetic materials to ...
3
votes
0
answers
24
views
How to get energy gap change at specific high symmetry points during ab initio molecular simulations?
I want to check the energy gap change during finite temperature ab initio molecular simulations. The results may look like this picture:
How can set aimd simulation to get this result in VASP?
1
vote
1
answer
46
views
When optimizing AMBER molecular dynamics parameters for protein-ligand interactions, how do I balance accuracy and efficiency?
Are there established practices for fine-tuning AMBER MD parameters in scenarios like mine, involving specific factors like mutant proteins or non-standard ligands? Has anyone successfully tackled ...
3
votes
1
answer
56
views
Adsorption of Li on GeC (optimization problem)
Actually I am trying to study the adsorption of Li or K on the sheet starting from the hexagon centre the system optimises with Li or (K) adsorbed over one of the Ge atoms and not remain at the centre....