Questions tagged [gromacs]
GROMACS software package to perform molecular dynamics, i.e. simulate the Newtonian equations of motion for systems with hundreds to millions of particles. Primarily designed for biochemical molecules like proteins, lipids and nucleic acids.
89
questions
18
votes
3
answers
657
views
Does anyone know any software for MD simulation which include the effect of magnetic field?
I am interested to generate some configurations from MD simulations which include the presence of strong magnetic field. Does anyone know any software that has some functionality to do this?
- 1,335
16
votes
1
answer
185
views
Why do interaction potential parametrizations of bulk systems not work for corresponding nano-systems?
I am studying several thermodynamics properties of crystalline iron oxides ($\ce{Fe3O4}$, $\gamma-\ce{Fe2O3}$) in their solid phases by MD simulations. At the moment, I am facing many difficulties in ...
- 852
14
votes
2
answers
888
views
Finding parameters for a molecule in various forcefields
Let's say I want to run a simulation for benzene in GROMACS using parameters from a certain force-field. I logged on to the Automated Topology Builder, I looked up benzene, and I found an ITP file for ...
- 4,141
14
votes
1
answer
115
views
How can I decide on the order of different monomers in a co-polymer?
While performing molecular dynamics simulations on a co-polymer, how long should the minimum size of the polymer be, to replicate the experimental properties with a reasonable error?
Within this ...
- 1,407
13
votes
3
answers
4k
views
What does "force field" mean?
In MD simulations we often hear the term force field.
What is it?
We know, for example that, Gromacs provides various types of forcefields.
- 1,674
13
votes
2
answers
511
views
How to find out the total functional form of the forcefield that is being used for a current simulation in GROMACS?
This is a doubt regarding the usage of various forcefields together in MD simulation of hydrates. For Eg., we use TIP4P/ice for water molecules, OPLS-UA for methane and EPM2 for carbon dioxide and so ...
- 1,407
12
votes
1
answer
741
views
Analyzing MD trajectory: molecule orientation
I'm running a MD simulation (in Gromacs) of dye molecules incorporated into polymer to determine how that system responds to external electric field.
There are two order parameters that are used to ...
- 2,363
12
votes
1
answer
2k
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 ...
- 493
11
votes
2
answers
1k
views
How to choose appropriate energy term in GROMACS MD simulation (Lennard-Jones or Coulomb)?
I was performing MD simulation for protein-ligand complex
following GROMACS tutorial.
I compared the energy of two ligands (S1 and S2) to a protein,
using GROMACS LJ-SR (Lennard-Jones Short Range) ...
- 1,674
11
votes
1
answer
1k
views
How do I calculate the radial distribution function from the centre of mass of unit cell in GROMACS?
I am simulating a droplet of water with GROMACS. The droplet consists of about 200 TIP3P water molecules and is placed near the center of a large, mostly empty box (...
- 6,591
11
votes
1
answer
318
views
What does vacuum signify in the simulation box after the simulation?
While carrying out constant volume simulations, I've observed that sometimes when the gaseous molecules come closer from their initial spacing, an empty space is created in the box. As I should ...
- 1,407
11
votes
1
answer
318
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 ...
- 4,141
10
votes
2
answers
894
views
How can I modify the Martini v2 or v3 force field to increase the Lennard Jones potential well depth (epsilon)?
I am trying to run some coarse-grained simulations of an Intrinsically Disordered Protein (IDP) using the Martini force field, and I want to increase the strength of interaction between water and ...
- 423
10
votes
3
answers
2k
views
How to perform "gmx make_ndx" with single command line bypassing the interactive shell?
Upon creating the index file with gmx make_ndx, we are usually
prompted with interactive shell.
For example:
...
- 1,674
10
votes
2
answers
1k
views
Molecular dynamics frame distorted by VMD?
I am running a VMD renderer to view my GROMACS molecular dynamics output.
At one frame, the geometry looks fairly normal like this:
At another frame, it displays a distorted geometry, with the ...
- 1,674
10
votes
1
answer
372
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,066
10
votes
1
answer
134
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
0
answers
382
views
Setup for 4-point flexible water model [closed]
I tried to search online if there are examples or suggestions on how to set up a 4-point flexible water model, such as TIP4P/2005f, qTIP4P/f, TIP4P/$\epsilon$ FLEX, etc., but there's no clear example ...
- 2,262
9
votes
2
answers
539
views
How to simulate peptide solubility using molecular dynamics (GROMACS)?
What's the suggested best practice to simulate peptide solubility?
I'm considering to use Gromacs Tutorial Lysozyme in Water as the model.
Later, I'm considering to use Radius of Gyration as measure ...
- 1,674
9
votes
1
answer
1k
views
How to manage disk space for Gromacs XTC (trajectory) file output
I'm running and MD run with this command:
gmx mdrun -v -deffnm md
It also produces trajectory file with .xtc extension.
...
- 1,674
9
votes
2
answers
2k
views
Large pressure fluctuations in molecular dynamics
I'm performing MD simulations, using NPT ensemble. The problem is that the pressure fluctuates severely and doesn't average at required level.
For example, I set the pressure at 1 bar and it averages ...
- 2,363
9
votes
1
answer
511
views
Is there a way to understand pi-pi interactions through GROMACS?
I have two polymer molecules with benzene rings as a part of them. I would like to understand how the interactions between them influences the solvability of the polymer in a non-aqeous solvent. I ...
- 1,407
9
votes
1
answer
194
views
Computational biology textbook similar to "Computational Materials Science" by June Gunn Lee
I'm looking for a text for computational biology which is similar in structure to a book called Computational Materials Science by June Gunn Lee.
In his book, Lee discusses the theory of MD ...
- 91
9
votes
1
answer
1k
views
Why does my calculation using GROMACS get stuck at step 0?
I have been using the HPCE on our campus for MD simulations using GROMACS.
However, when I use the following script, after the job starts to run on the cluster, gets stuck at Step 0 in the log file.
<...
- 1,407
9
votes
1
answer
653
views
What is the purpose of DSSP when using martinize2?
I am studying an IDP with a lot of helical secondary structure, and trying to find a trajectory that I can use as an ensemble that fits my SAXS data. My starting models are the output from the Robetta ...
- 423
9
votes
1
answer
97
views
Generating topology for charged molecules
I have to generate a topology file for a molecule like Sodium Triphosphate, and I do not find forcefields available in the literature. So as a starting point would like to use engines like LigParGen ...
- 1,407
8
votes
2
answers
799
views
Confused about the application/definition of improper torsions
The only basic definition of improper torsion/improper dihedral I have is - improper dihedral is formed when we have a central atom bonded to the three atoms. Can anyone help me visualize/help to ...
- 1,066
8
votes
1
answer
204
views
How to address the errors on Automated Topology Builder?
I was trying to get the charges and the topology file for this molecule ($\ce{C28H58O7S8}$).
But I get the following error despite using the right number of hydrogens.
...
- 1,407
8
votes
1
answer
482
views
Alternative to CGenFF for generating large ligand topology
I am now preparing a ligand topology following GROMACS protein-ligand tutorial.
However, I have a large ligand (446 atoms).
$ wc -l my_ligand.pdb
446 my_ligand.pdb
...
- 1,674
8
votes
1
answer
143
views
Very slow performance for an OPLS-AA methanol simulation
I am running a methanol simulation with a topology I obtained from Tinker for OPLS-AA, at 298 K and 1 bar. This is my .itp file:
...
- 4,141
8
votes
1
answer
727
views
Viscosity calculation using molecular dynamics simulations
I use the GROMACS software for MD simulations and I am trying to understand the calculation of viscosity using combinations of different ensembles and different methods.
For a non-equilibrium MD ...
- 1,407
8
votes
0
answers
180
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,066
8
votes
0
answers
287
views
How to add double bonds to hexagons when creating a graphene layer in VMD? [closed]
When creating a Graphene layer in VMD using the Nanotube Builder it only generates hexagons with single bonds because of this errors occur when creating the topology for gromacs software using,
...
- 305
8
votes
0
answers
123
views
Calculate FF interaction energy of small molecular dimers [closed]
I have a set of molecular dimers (~25 atoms) for which DFT-D3 interaction energies are calculated by Psi4 and Turbomole program. Now I want to compare these values against Force-Field methods like UFF/...
- 810
8
votes
0
answers
130
views
How to go from tinker parameters to a GROMACS itp file, specifically for dihedrals? [closed]
I am trying to write an itp file for the compound 1,2-diethoxyethance (CCOCCOCC).I used tinker analyze to obtain the parameter file, f2.out, which I have attached to this post. In order to check if my ...
- 4,141
7
votes
2
answers
184
views
Import an initial trajectory into OPS without OpenMM
I'm trying to use openpathsampling to simulate the binding of a ligand to my enzyme. I'm using Gromacs for my normal simulation on the HPC system provided by my university. I created an initial ...
- 143
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
7
votes
1
answer
72
views
Keeping a coarse grained particle off plane (Molecular Dynamics)
I have a coarse-grained particle (bead) that I want to keep off the plane from the plane of other beads which form like a 2D plane. Do people employ improper dihedral for this or any other bonded ...
- 1,066
7
votes
1
answer
233
views
How to go from Tinker parameter file to GROMACS topology file for improper and proper dihedrals?
I have been playing with a polymer called glyme, and I have run Tinker analyze on it to get OPLS-AA parameters for it. I have attached the output parameter file to this post.
My question is, how do I ...
- 4,141
7
votes
1
answer
535
views
7
votes
1
answer
86
views
How do I setup my MSTIS simulation correctly?
I'm trying to setup a Multi State Transition Interface Sampling (MSTIS) simulation to study the ligand binding of my enzyme using OPS and GROMACS. My enzyme acts on small sugar chains and I would like ...
- 143
7
votes
1
answer
679
views
Is it possible to simulate the Raman spectra via Molecular Dynamics (MD) using the GROMACS package?
I'm new to molecular dynamics and DFT. The aim was to study the Raman spectra of a molecule in a medium. As we know the Raman spectra of a molecule change due to the induced polarization in the ...
- 4,448
7
votes
0
answers
148
views
How to determine minimum reasonable MD timing for ligand-protein complex [closed]
I am working on protein-ligand complex MD simulation using GROMACS.
Now I'm running the tutorial.
...
- 1,674
7
votes
0
answers
162
views
How to effectively simulate high density systems in GROMACS? [closed]
I have a system with a dense liquid solvent and a polymer solute. I have inserted all the molecules randomly using gmx insert-molecules and then did the energy ...
- 1,407
7
votes
0
answers
136
views
How to do molecular dynamics calculations and analyse trajectories of small molecules? [closed]
I need to perform some MD calculations and then trajectories of some small molecules analyzed.
What I have is
protein
cofactor (FAD)
small molecule (either single O2 molecule or single Chlorine atom)
...
6
votes
2
answers
1k
views
How to initiate protein-ligand docking with MD simulation tools
I have two PDB files (each corresponding to a ligand and a receptor).
Each downloaded separately from PDB.
I'd like to perform ligand to protein docking process with MD with GROMACS or Amber.
Using ...
- 1,674
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
6
votes
1
answer
1k
views
How to visualize protein-ligand complex MD run using VMD
GROMACS version: gmx, version 2021.4
GROMACS modification: No
I'm following the Protein-Ligand Complex tutorial.
I'd like to visualize the binding using VMD.
I use the output from these commands:
<...
- 1,674
6
votes
1
answer
2k
views
Error in GROMACS for adding ions to a solution
I am following Justin Lemkul's GROMACS simulation given on this website here for GROMACS.
I believe I have followed the process to the T. I would like to mention that in the "Generate Topology&...
- 4,141
6
votes
1
answer
1k
views
How to add molecules to a specific location in periodic box using gromacs?
In gromacs using the commands, command_1
gmx insert-molecules -f run03.pdb -ci CL.pdb -o out.pdb -nmol 119
will distribute CL ions in the system randomly, but how ...
- 305