# Pressure coupling

**URL:** https://gromacs.bioexcel.eu/t/pressure-coupling/2088
**Category:** User discussions
**Tags:** mdp-parameters
**Created:** [May 3, 2021, 10:34pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088 "2021-05-03T22:34:37Z")
**Posts on this page:** 7
**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: [May 3, 2021, 10:34pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/1 "2021-05-03T22:34:37Z")

</div>

Dear all,

I have used deform option as in the following:  
**deform = 0 0 0 7.0 0 0**

But, for pressure coupling I am not sure if I have written it correctly.  
**pcoupl = Parrinello-Rahman**  
**pcoupltype = anisotropic**  
**tau\_p = 2.0**  
**ref\_p = 0.0 0.0 0.0 1.0 0.0 0.0**  
**compressibility = 0.0 0.0 0.0 46e-5 0.0 0.0**

Could you please tell me if it is correct?  
And, if not, could you please direct me to write the correct one?

I appreciate your assistance in advance.  
Best

---

<div class="post-metadata">

### Author: ![MichelePellegrino](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/michelepellegrino/32/628_2.png) [@MichelePellegrino](https://gromacs.bioexcel.eu/u/MichelePellegrino)
#### Post date: [May 4, 2021, 8:36am UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/2 "2021-05-04T08:36:20Z")

</div>

Hi MOH4,

it’s correct syntax-wise: have you tried to run it? Does it work?

Best,  
Michele

---

<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: [May 4, 2021, 5:41pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/3 "2021-05-04T17:41:31Z")

</div>

Hi Michele,

Not yet. I am using another option along with this.  
But, I was not sure if what I had written (particularly zero values) were correct.

Could you please tell me what the zero values of compressibility mean?

Best

---

<div class="post-metadata">

### Author: ![MichelePellegrino](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/michelepellegrino/32/628_2.png) [@MichelePellegrino](https://gromacs.bioexcel.eu/u/MichelePellegrino)
#### Post date: [May 4, 2021, 6:02pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/4 "2021-05-04T18:02:47Z")

</div>

The definition of isotropic compressibility is the ratio between the derivative of volume w.r.t. pressure and total volume.  
In a non-isotropic case you may have a different compressibility value in each direction (or rather: for each component of the pressure tensor), however the physical interpreation kinda remains the same (‘how much’ the specific deformation is sensitive to pressure differences).  
Saying that compressibility is zero means that the box does not get deformed in that direction. For example, if you have a cubic box and compressibilty is non-zero only along ZZ, then the box would only be strained/compressed along z.

---

<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: [May 4, 2021, 7:11pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/5 "2021-05-04T19:11:27Z")

</div>

Many thanks for your informative response.

In this regard, if I write **deform = 0 0 0 7.0 0 0** , it means the box will be deformed along angle between xy directions as it is LyXY. Is it correct?

---

<div class="post-metadata">

### Author: ![MichelePellegrino](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/michelepellegrino/32/628_2.png) [@MichelePellegrino](https://gromacs.bioexcel.eu/u/MichelePellegrino)
#### Post date: [May 5, 2021, 6:49am UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/6 "2021-05-05T06:49:11Z")

</div>

Hi,

the off-diagon terms refer to the vectors b(x) c(x) c(y) in case of a triclinic box; please have a look at [Periodic boundary conditions — GROMACS 2021 documentation](https://manual.gromacs.org/2021/reference-manual/algorithms/periodic-boundary-conditions.html) and [Molecular dynamics parameters (.mdp options) — GROMACS 2021 documentation](https://manual.gromacs.org/documentation/2021/user-guide/mdp-options.html).  
I personally never utilized off-diagonal deformations, so I am sorry but that is the only indication I am able to provide.

Cheers,  
Michele

---

<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: [May 5, 2021, 1:40pm UTC](https://gromacs.bioexcel.eu/t/pressure-coupling/2088/7 "2021-05-05T13:40:32Z")

</div>

Hi Michele,  
Many thanks for your all responses.  
They are really helpful.

For the last question, could you please tell me if you know the range of value that I can give to deform?  
I mean what is the maximum value for deform input components?  
In the other words, is this possible to write e.g. deform = 0 0 0 77.0 0 0 or values more than this?

I apologize for asking many questions.  
Many thanks in advance,  
Best
