Charmm36-jul2021.ff virtual sites

GROMACS version:2021.5
GROMACS modification: Yes/No
Here post your question
I was going prepare a simulation and modelling residues with virtual site interactions using the charmm36-jul2021.ff, but with pdb2gmx I have gotten the following error:
“Fatal error:
A dummy mass for type CT3 bonded to type C is required, but no virtual site
database (.vsd) files where found”
Indeed there is no *.vsd file present in the ff directory.
Furthermore, in the charmm36-feb2021.ff release the ffbonded.itp has a section [ constrainttypes ] under which constraints imposed for dummy atoms. This section is also missing from the jul2021 release.
I wonder, would it be safe to copy over the constrainttypes section for the feb2021 to jul2021 and also copy merged.vsd to jul2021.ff/aminoacid.vsd?
Sandor

Dear Sándor,

Virtual site information was omitted from the release you mentioned, but in principle you can copy it from the previous one. The same applies for the constraints. We plan to reintroduce these in an upcoming release of the force field.

The reason: virtual site data and the corresponding are not present in the upstream CHARMM force field distribution. @jalemkul, please correct me if I’m wrong.

Kind regards,

András

Hello Andras,
Yeah, I have thought that I can copy over, the file and constraints, but I wanted to make sure if it would be correct and, sort of, draw attention to it.
thanks for the quick response
Sandor

It should be safe to copy it but just be aware that this file is entirely a GROMACS construct and has nothing to do with the actual CHARMM force field, hence why it got left out of the new port, which was generated in a more automated fashion from the CHARMM force field files. We used to write one as a boilerplate entry from the charmm27.ff file, but we don’t test it in any way.