Fixing bugs in FF14SB port for Gromacs
Bugs? The newly released amber14sb made modifications for protein, DNA and RNA. This force field could be ported into gromacs as I already introduced previous. However, there are some small bugs in the improper dihedral angles of HID, HIE and HIP. So when generating tpr file, errors would be prompted. When run grompp -f npt.mdp -p topol.top -c protein.gro -o product.tpr -n index, errors occur: ERROR 1 [file topolff14sb_Protein_chain_B.itp, line 208041]: No default Improper Dih. types The reason is that the parameters for some dihedral angles are missing. So we need to fix this. How to fix? Add the following two lines in the ffbonded.itp after the 730th lines: NA CW CC CT 4 180.00 4.60240 2 ;;; HID force field added NA CV CC CT 4 180.00 4.60240 2 ;;; HID force fie...