# Umbrella sampling for two sheets

GROMACS version: 2019.4
GROMACS modification: No

I am calculating the PMF for two sheets in solution with different interlayer distances (see figures below). The problem is when sheets are in a certain range (between 0.5 to 0.65 nm) they will move toward each other so the distance will be decreased to 0.4 nm so there will not be any window for this range. The results are shown in the figure below.

The question is how can I sample from this range which is not stable and sheets will absorb each other?

More info: I restrained one sheet, the other sheet is free but there’s an umbrella potential between two sheets. Th epull_coord1_k is 1000. Changing k doesn’t help too. Using “pull coord1_type = constraint” will deform the free sheet heavily.

The procedure is something like this figure.

Hello, have you solved this question? I also met this question and want some answers.
I pull one sheet along z axis, when I select some conf.gro to execute US, the upper sheet would decrese heavily so that go to the original position. I can’t calculate the PMF when this situation appeared.
Thank you

To my knowledge, you can’t calculate the PMF along with this reaction coordinate with US. You should consider other reaction coordinates such as :

Regards,

Ehsan

I have read some articles about peeling one sheet to achieve exfoliation and calculating the PMF along the z axis. They usually use ABF from lammps. Why US from GROMACS can not achieve this calculation?
https://pubs.rsc.org/en/content/articlelanding/2019/CP/C9CP03654G#!divAbstract
Best regards

I’m not familiar with ABF, and I’m not sure about the reason behind this but I think it has something to do with the force required to exfoliate sheets in this direction. Whatever the reason is, I couldn’t calculate PMF either with US or AWH method.

Thanks for your patiently answer. I would try the calculation using US again by adjusting the force.

Best regards

You’re welcome, I don’t think changing the force can solve the problem. I guess this is not the appropriate reaction coordinate. BTW, if you could solve the issue, please let me know.

Thanks.

Hello, I change the force and the upper sheet can not move to the original position. You can see my US histograms, but the PMF still can not calculate in this way. Have you met this situation?
Thanks a lot.

I’m not sure if you “position restrained” the bottom sheet or top sheet but whatever you did you shouldn’t restrain both of them because then you have no movement in your system and no PMF. What you have to do is restraining all the atoms in one sheet and then defining an umbrella potential between the COM of the two sheets. I recommend defining another reaction coordinate because I couldn’t get any reasonable answer with this reaction coordinate. Not all the reaction coordinates are feasible.

Regards,

Ehsan