# Problem in the umbrella histrograms **URL:** \ **Category:** User discussions\ **Tags:** umbrella-sampling\ **Created:** [September 15, 2025, 8:37pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647 "2025-09-15T20:37:56Z")\ **Posts on this page:** 6\ **Page:** 1
**Author:** ![Eduardo1](https://avatars.discourse-cdn.com/v4/letter/e/838e76/32.png) [@Eduardo1](https://gromacs.bioexcel.eu/u/Eduardo1)\ **Post date:** [September 15, 2025, 8:37pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/1 "2025-09-15T20:37:56Z")
GROMACS version: 2019.4 GROMACS modification: No Dear all, After pulling a polymer through a membrane using steered molecular dynamics simulation, I selected the frames of this trajectory using the pull-x.xvg output as the initial configurations for the umbrella sampling simulations. A part of the list of these 111 frames x COM distance (z axis) is shown below: ```auto 1 10.768 2 10.464 3 10.202 4 9.919 5 9.587 6 9.313 7 8.984 8 8.689 9 8.411 10 8.113 ``` However, after running umbrella sampling simulation and the wham, I noticed that the first histogram (black one below) is not located (centralized) at 10.768 nm as shown in the list above. The umbrella histograms are shown below: ![histogram-problem](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/f/febc7d5614e90f339bbdd9edf4d7130b0dc76c63.png) Should I correct the position of this histogram by increasing the force constant in the umbrella sampling simulation (I used k = 1000 kJ/mol/nm)? I also have many gaps in these histogram despite the fact that the COM distances of each configuration were equally spaced by 0.3 nm as shown in the list above ( a part of this list of 111 configurations). Can the number of bins used in the wham calculation correct this problem or it is a question of increase the force constant in the US simulations? Sincerely, Eduardo. ---
**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\ **Post date:** [September 16, 2025, 7:45am UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/2 "2025-09-16T07:45:37Z")
Your force constant is far too low. Note that pulling a large molecule through a membrane is a very challenging calculation. You might not get convergence with the simulation times available. ---
**Author:** ![Eduardo1](https://avatars.discourse-cdn.com/v4/letter/e/838e76/32.png) [@Eduardo1](https://gromacs.bioexcel.eu/u/Eduardo1)\ **Post date:** [September 16, 2025, 1:13pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/3 "2025-09-16T13:13:07Z")
Thank you for your response. If I increase the force constant in the US simulations, can the histogram gap problem get worse? Do you have a suggestion on what this increase should be? ---
**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\ **Post date:** [September 16, 2025, 1:34pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/4 "2025-09-16T13:34:16Z")
The distributions move towards the center of the potential when increasing the force constant. But you will need many more umbrella. The width of the distributions is proportional to 1/sqrt(k). ---
**Author:** ![Eduardo1](https://avatars.discourse-cdn.com/v4/letter/e/838e76/32.png) [@Eduardo1](https://gromacs.bioexcel.eu/u/Eduardo1)\ **Post date:** [September 16, 2025, 2:05pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/5 "2025-09-16T14:05:42Z")
Thank you! I will increase the force constant and the number of umbrella in the regions with gaps. Should the increase in the force constant be done in the pulling simulation (steered molecular dynamics) or in the umbrella sampling simulations? ---
**Author:** ![hess](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/hess/32/416_2.png) [@hess](https://gromacs.bioexcel.eu/u/hess)\ **Post date:** [September 17, 2025, 3:35pm UTC](https://gromacs.bioexcel.eu/t/problem-in-the-umbrella-histrograms/12647/6 "2025-09-17T15:35:08Z")
Strictly speaking only in the US simulations. But also, or maybe even more, in a steered simulation there is the risk that the system will lag and be far from the center of the umbrella.