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.

Filter by
Sorted by
Tagged with
4
votes
0answers
51 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 non equilibrium MD ...
6
votes
0answers
43 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: ...
7
votes
0answers
75 views

Setup for 4-point flexible water model

I tried to search online if there are examples or suggestions on how to setup a 4-point flexible water model, such as TIP4P/2005f, qTIP4P/f, TIP4P/$\epsilon$ FLEX, etc. but there's no clear example on ...
6
votes
0answers
47 views

Back-mapping Martini 3 coarse-grain to all-atom representations

I have a $10\mu s$ trajectory with 10,000 frames that I would like to backmap to atomistic representations so that I can compare my ensemble to experimental measurements. I am using MartiniTools for ...
9
votes
0answers
102 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 ...
10
votes
2answers
417 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 ...
10
votes
1answer
111 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 ...
8
votes
1answer
134 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. ...
7
votes
0answers
76 views

How to go from tinker parameters to a GROMACS itp file, specifically for dihedrals?

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 ...
7
votes
1answer
52 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 ...
3
votes
1answer
44 views

Missing residues in ATB downloaded pdb/rtp files

I have downloaded the diglyme molecule from the ATB server (http://atb.uq.edu.au/molecule.py?molid=34570#panel-md) Two files are available: a PDB file and a topology ITP/RTP file. The name of the ...
8
votes
2answers
239 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 ...
12
votes
1answer
188 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 ...
10
votes
1answer
232 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 (...
11
votes
1answer
117 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 ...
14
votes
2answers
167 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 ...
5
votes
1answer
341 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&...
5
votes
1answer
57 views

How to overcome gromacs make_ndx automatic ranking of atomic indexes?

When working with gmx make_ndx I found very annoying the fact that atomic indexes are automatically ranked in the outputted index ...
12
votes
1answer
430 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 ...
7
votes
0answers
81 views

How to effectively simulate high density systems in GROMACS?

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 ...
9
votes
1answer
117 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 ...
8
votes
1answer
58 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 ...
9
votes
1answer
340 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. <...
6
votes
1answer
139 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 ...
9
votes
1answer
152 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 ...
7
votes
0answers
102 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) ...
15
votes
1answer
161 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 ...
14
votes
1answer
94 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 ...
17
votes
3answers
215 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?
13
votes
2answers
243 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 ...