# Energy minimization does not work

**URL:** https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109
**Category:** User discussions
**Created:** [January 22, 2025, 1:40pm UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109 "2025-01-22T13:40:46Z")
**Posts on this page:** 9
**Page:** 1

<div class="post-metadata">

### Author: ![XiaoDingdang](https://avatars.discourse-cdn.com/v4/letter/x/a9a28c/32.png) [@XiaoDingdang](https://gromacs.bioexcel.eu/u/XiaoDingdang)
#### Post date: [January 22, 2025, 1:40pm UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/1 "2025-01-22T13:40:46Z")

</div>

GROMACS version:  
GROMACS modification: Yes/No  
Here post your question  
[com.top](https://gromacs.bioexcel.eu/uploads/short-url/zexr0bk5i6iZsZ3aBEIlzaG5qVY.top) (22.4 KB)

Non-default thread affinity set probably by the OpenMP library,  
disabling internal thread affinity

Back Off! I just backed up em1.trr to ./#em1.trr.3#

Back Off! I just backed up em1.edr to ./#em1.edr.3#

Steepest Descents:  
Tolerance (Fmax) = 1.00000e+03  
Number of steps = 50000  
Step= 0, Dmax= 1.0e-02 nm, Epot= 1.05915e+34 Fmax= 7.12759e+07, atom= 7991  
Is it because of structural problems that it cannot be operated when energy is minimized?

---

<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: [January 23, 2025, 7:25am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/2 "2025-01-23T07:25:42Z")

</div>

I would think that it’s due to atoms being too close to each other or a bond that is too long. A good start would be to have a look at atom 7991, all atoms in the same residue as atom 7991 and all atoms within 0.4 nm (0.6 if running a coarse-grained simulation) of atom 7991.

---

<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: [January 24, 2025, 6:49am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/4 "2025-01-24T06:49:59Z")

</div>

I’d suggest opening the coordinate file (.gro) in, e.g., VMD. And use the representation options (`Graphics` → `Represenations...`) to show only the atoms you are interested in, e.g.  
`same residue as within 4 of index 7990`.  
This will show all residues with at least one atom within 0.4 nm of atom 7991 (index 7990 in VMD). Then see if all bonds you expect are correctly represented and if there are any atoms in close vicinity. You may want to show the atoms as VdW spheres to make that slightly easier.

---

<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: [January 24, 2025, 6:51am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/5 "2025-01-24T06:51:06Z")

</div>

> [@XiaoDingdang](#):
>
> 7991SOL OW28077 0.539 7.512 7.960  
> 7991SOL HW128078 0.458 7.513 8.018  
> 7991SOL HW228079 0.542 7.595 7.905  
> This is my 7991 atom.

No, that is your residue 7991, containing atoms 28077, 28078 and 28079.

---

<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: [January 24, 2025, 7:06am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/7 "2025-01-24T07:06:07Z")

</div>

That sounds very strange. Do you know if the problem is with VMD or with your gro files?

You’ll probably want some kind of visualization software, so if VMD doesn’t work for you, try another one. Most of them would be able to do similar visualizations, but you’d probably have to check their documentation.

---

<div class="post-metadata">

### Author: ![XiaoDingdang](https://avatars.discourse-cdn.com/v4/letter/x/a9a28c/32.png) [@XiaoDingdang](https://gromacs.bioexcel.eu/u/XiaoDingdang)
#### Post date: [January 24, 2025, 7:10am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/8 "2025-01-24T07:10:50Z")

</div>

Now that I’ve successfully opened it, where exactly should I find that atom?

---

<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: [January 24, 2025, 7:11am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/9 "2025-01-24T07:11:48Z")

</div>

I don’t know where your atoms are located. Did you follow my advice above?

---

<div class="post-metadata">

### Author: ![XiaoDingdang](https://avatars.discourse-cdn.com/v4/letter/x/a9a28c/32.png) [@XiaoDingdang](https://gromacs.bioexcel.eu/u/XiaoDingdang)
#### Post date: [January 24, 2025, 7:36am UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/10 "2025-01-24T07:36:14Z")

</div>

Yes, I followed your advice, but I didn’t look for where the atoms were

---

<div class="post-metadata">

### Author: ![XiaoDingdang](https://avatars.discourse-cdn.com/v4/letter/x/a9a28c/32.png) [@XiaoDingdang](https://gromacs.bioexcel.eu/u/XiaoDingdang)
#### Post date: [February 1, 2025, 1:27pm UTC](https://gromacs.bioexcel.eu/t/energy-minimization-does-not-work/11109/11 "2025-02-01T13:27:53Z")

</div>

Could you please elaborate on that more?
