Position restraints and barostat in alchemical free energy double-system/single-box setup

User discussions
, ,
1

GROMACS version: 2023.2
GROMACS modification: plumed_2.10.0_dev

I’m performing calculation protein stability change upon charge modifying mutation. For this reason, I prepared a double-system/single-box setup to maintain constant charge. As suggested in the related chapter, I’m applying position restraints to the COM atom for each of the two protein copies. At the equilibration stage, I get the following warnings, depending on a barostat I use:

For Parrinello-Rahman

You are combining position restraints with Parrinello-Rahman pressure
coupling, which can lead to instabilities. If you really want to combine
position restraints with pressure coupling, we suggest to use Berendsen
pressure coupling instead.

For Berendsen

The Berendsen barostat does not generate any strictly correct ensemble,
and should not be used for new production simulations (in our opinion).
We recommend using the C-rescale barostat instead.

And no barostat-related warning for C-rescale.
So, which is the better choice of barostat for this type of calculations?

Thank you in advance,
Sofya

2

We would always advice to use the C-rescale barostat.

We forgot to update this error message when deprecating the Berendsen barostat.

2 Likes
3

Thank you!