GROMACS version: 2021
GROMACS modification: No
I’m running a simulation with a unit cell with parameters: -box 3.9970 2.6950 67.7900 -angles 89.24 94.59 105.58. This creates an issue with ‘Triclinic box too skewed’ as |zx| > 0.5 |xx| (|-5.423| > 0.5 * 3.997). As I understand it, this shouldn’t be an issue as we could add the a vector to the c vector to then reduce the value of |zx|, and Gromacs seems to do this after using genion.
However after viewing trajectories in vmd, the system does not seem to follow the correct unit cell and gmx solvate does not seem to place the solvent inside the correct box (instead it seems to follow a unit cell in which zx is equal to 0). Is this just an issue with the way gromacs is outputting the trajectories and during the simulation it will follow the correct unit cell? Or is the triclinic box being too skewed mean the simulation is using the wrong unit cell?
Thanks for any help.