Q6 icon indicating copy to clipboard operation
Q6 copied to clipboard

Question about `number of paths` parameter in Q's forcefield.prm definition

Open David-Araripe opened this issue 1 year ago • 0 comments
trafficstars

I'm working on translating some forcefield amino acid parameters to Q's format, but I have a question about one of the parameters used in the [torsions] section of the .prm file. According to the Qmanual, the 8th column is the "number of paths", which seems to be defined as:

is the number of ways that a two-atom torsion can be defined, i:e: the product of the number of atoms bonded to the two middle atoms. It is used to distribute the force over all the atoms involved.

However, this doesn't seem to match the values. I would understand it if all the carbon-carbon bonds would be defaulted to 9 but that doesn't seem to be the case as you can see here: https://github.com/qusers/Q6/blob/master/ff/amber14sb/qamber14.prm#L759-L762

Further, I can see that whenever wildcard atom flags are present, these values seem to be able to go higher than 1.0. However, for impropers where you have the 4 atom types defined, this variable always seem to default to 1.0, as seen here:

https://github.com/qusers/Q6/blob/master/ff/amber14sb/qamber14.prm#L856-L861

I was looking at the other forcefield definitions to understand this, but for OPLS2015 it seems to be always left as 1.0, as seen here:

https://github.com/qusers/Q6/blob/master/ff/oplsaam2015/qoplsaa.prm#L749-L767

I then tried looking at the Fortran code to see how was this #paths variable used but I couldn't find it. Is this some sort of deprecated parameter? And is it fine if I would define a Amber-based forcefield using this #paths parameter as 1.0 for all the torsions defined?

David-Araripe avatar Jun 24 '24 08:06 David-Araripe