# Pressure scaling more than 1%

**URL:** <https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072>\
**Category:** User discussions\
**Created:** [April 30, 2021, 8:12pm UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072 "2021-04-30T20:12:26Z")\
**Posts on this page:** 6\
**Page:** 1

<div class="post-metadata">

**Author:** ![MOH4](https://avatars.discourse-cdn.com/v4/letter/m/a87d85/32.png) [@MOH4](https://gromacs.bioexcel.eu/u/MOH4)\
**Post date:** [April 30, 2021, 8:12pm UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/1 "2021-04-30T20:12:26Z")

</div>

Dear all,

I encountered this error:  
Warning: Pressure scaling more than 1%. This may mean your system  
is not yet equilibrated. Use of Parrinello-Rahman pressure coupling during  
equilibration can lead to simulation instability, and is discouraged.

I saw gromacs explanation to solve this error.

But my question is that, does this error mostly assigned to NPT or prvious steps may bring about this error?

I apperciate your responses in advance,  
Best

---

<div class="post-metadata">

**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\
**Post date:** [April 30, 2021, 8:44pm UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/2 "2021-04-30T20:44:38Z")

</div>

Without knowing what you’ve done to prepare the system, it’s hard to say specifically, but the Parrinello-Rahman barostat is prone to large swings in systems that are not adequately equilibrated already so you may want to try using the Berendsen method for the barostat first, then switching to P-R afterwards.

---

<div class="post-metadata">

**Author:** ![MOH4](https://avatars.discourse-cdn.com/v4/letter/m/a87d85/32.png) [@MOH4](https://gromacs.bioexcel.eu/u/MOH4)\
**Post date:** [April 30, 2021, 9:13pm UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/3 "2021-04-30T21:13:56Z")

</div>

[npt.mdp](https://gromacs.bioexcel.eu/uploads/short-url/ixglKvOYmemNJSFkY8WMTqfMkk2.mdp) (2.7 KB)

Dear Prof. Lemkul,  
Thanks for your responses.

Will the npt.mdp that I enclosed work?  
I am performing the simulation of DNA-ligand in water.  
And, could you please explain how to switch from Berendsen to Par\_R?  
I know only how to use a npt.mdp.

I appreciate your help.  
Best

---

<div class="post-metadata">

**Author:** ![SpongeBob](https://avatars.discourse-cdn.com/v4/letter/s/90db22/32.png) [@SpongeBob](https://gromacs.bioexcel.eu/u/SpongeBob)\
**Post date:** [May 3, 2021, 9:12am UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/4 "2021-05-03T09:12:46Z")

</div>

There is an option “pcoupl = Berendsen” in your npt.mdp, you can use this file to equilebrated the system first, then switching to “pcoupl = Parrinello-Rahman”

---

<div class="post-metadata">

**Author:** ![mirror](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@mirror](https://gromacs.bioexcel.eu/u/mirror)\
**Post date:** [July 5, 2023, 9:03am UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/5 "2023-07-05T09:03:27Z")

</div>

Excuse me, in the 2023 version of gromacs, I first used pcoupl = Berendsen to simulate 100ps in the pre-balance, and then when I used pcoupl = Parrinello-Rahman, an error still occurred: pressure scaling more than 1%. This may mean your system is not yet equilibrated. Use of Parrinello-Rahman pressure coupling during equilibration can lead to simulation instability, and is discouraged. Is this related to the 2023 version of gromacs I use?

---

<div class="post-metadata">

**Author:** ![MagnusL](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/magnusl/32/2380_2.png) [@MagnusL](https://gromacs.bioexcel.eu/u/MagnusL)\
**Post date:** [September 21, 2023, 10:11am UTC](https://gromacs.bioexcel.eu/t/pressure-scaling-more-than-1/2072/7 "2023-09-21T10:11:16Z")

</div>

Equilibrating for 100ps with the Berendsen barostat might not be enough. It depends on the system. However, after the `c-rescale` barostat was added, it is recommended to use that both for equilibration and production simulations.
