# Segmentation falult (core dumped

**URL:** https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120
**Category:** User discussions
**Tags:** grompp
**Created:** [May 10, 2024, 4:28am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120 "2024-05-10T04:28:07Z")
**Posts on this page:** 8
**Page:** 1

<div class="post-metadata">

### Author: ![bishnu](https://avatars.discourse-cdn.com/v4/letter/b/43a26b/32.png) [@bishnu](https://gromacs.bioexcel.eu/u/bishnu)
#### Post date: [May 10, 2024, 4:28am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/1 "2024-05-10T04:28:07Z")

</div>

GROMACS version:  
GROMACS modification: Yes/No  
Here post your question  
hello i am using gromacs 24 and calculating diffusion of a drug molecule byut i encountered this segmentation fault error during energy minimization can anyone help me?  
gmx mdrun -nt 8 -v -nb gpu -gpu\_id 0 -s em.tpr -deffnm em  
:-) GROMACS - gmx mdrun, 2024.1 (-:

Executable: /usr/local/gromacs/bin/gmx  
Data prefix: /usr/local/gromacs  
Working dir: /home/administrator/Desktop/  
Command line:  
gmx mdrun -nt 8 -v -nb gpu -s em.tpr -deffnm em

Back Off! I just backed up em.log to ./#em.log.4#  
Reading file em.tpr, VERSION 2024.1 (single precision)  
1 GPU selected for this run.  
Mapping of GPU IDs to the 1 GPU task in the 1 rank on this node:  
PP:0  
PP tasks will do (non-perturbed) short-ranged interactions on the GPU  
PP task will update and constrain coordinates on the CPU  
Using 1 MPI thread  
Using 8 OpenMP threads

NOTE: The number of threads is not equal to the number of (logical) cpus  
and the -pin option is set to auto: will not pin threads to cpus.  
This can lead to significant performance degradation.  
Consider using -pin on (and -pinoffset in case you run multiple jobs).

Back Off! I just backed up em.trr to ./#em.trr.4#

Back Off! I just backed up em.edr to ./#em.edr.4#

Steepest Descents:  
Tolerance (Fmax) = 1.50000e+02  
Number of steps = 5000000  
Step= 0, Dmax= 1.0e-03 nm, Epot= 1.98031e+33 Fmax= 8.15762e+06, atom= 14  
Segmentation fault (core dumped)

---

<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: [May 10, 2024, 8:34am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/2 "2024-05-10T08:34:04Z")

</div>

There are probably atoms overlapping causing very high forces and the extremely high potential energy in the system.

---

<div class="post-metadata">

### Author: ![bishnu](https://avatars.discourse-cdn.com/v4/letter/b/43a26b/32.png) [@bishnu](https://gromacs.bioexcel.eu/u/bishnu)
#### Post date: [May 11, 2024, 4:16pm UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/3 "2024-05-11T16:16:03Z")

</div>

thank you for your suggestion i added h-bond constraints in mdp file and it worked perfectly.

---

<div class="post-metadata">

### Author: ![DikshaVats10](https://avatars.discourse-cdn.com/v4/letter/d/a8b319/32.png) [@DikshaVats10](https://gromacs.bioexcel.eu/u/DikshaVats10)
#### Post date: [June 4, 2024, 6:20am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/4 "2024-06-04T06:20:42Z")

</div>

hi MagnusL, i am trying to minimize a small molecule with 73 atom, during enrgy minimization after certain point, 1 hydrogen atom is getting attached to additional oxgyen atom present nearby. i have check the positions of these two atoms in em.gro file, its very close. is it overlapping of atoms?

---

<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: [June 4, 2024, 1:03pm UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/5 "2024-06-04T13:03:42Z")

</div>

Hi,  
When you say that “1 hydrogen atom is getting attached to additional oxygen atom”, I presume you mean that if you visualize the molecule in a software with automatic bond detection there is a bond shown between them, right? If so, that does not mean that it is getting attached, or that they are overlapping. But it is still worrying that energy minimization would put the atoms that close to each other. That indicates a very strained, or unfavorable initial geometry.

It would be improbable that energy minimization would make two atoms actually overlap. The energy barrier for that would be very high and in that case there are definitely other problems in your structure.

---

<div class="post-metadata">

### Author: ![DikshaVats10](https://avatars.discourse-cdn.com/v4/letter/d/a8b319/32.png) [@DikshaVats10](https://gromacs.bioexcel.eu/u/DikshaVats10)
#### Post date: [June 4, 2024, 2:32pm UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/6 "2024-06-04T14:32:37Z")

</div>

@MagnusL Thanks for the clarity.  
yes, may be in visualization it is appearing as a bond due to closeness of two atoms.  
if i ignore this visualization and proceed further then after minimizing energy to possible minimum value, nvt run is showing that “atoms are rotating more than 30 degree(segmentation fault, core dumped)”, it means my system is not stable? how to resolve this?

---

<div class="post-metadata">

### Author: ![DikshaVats10](https://avatars.discourse-cdn.com/v4/letter/d/a8b319/32.png) [@DikshaVats10](https://gromacs.bioexcel.eu/u/DikshaVats10)
#### Post date: [June 5, 2024, 8:51am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/7 "2024-06-05T08:51:49Z")

</div>

@jalemkul could you please help ij this?

---

<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: [June 5, 2024, 9:20am UTC](https://gromacs.bioexcel.eu/t/segmentation-falult-core-dumped/9120/8 "2024-06-05T09:20:24Z")

</div>

If you get a crash in the NVT run, after a “successful” EM run, you have probably ended up in a local unfavourable configuration. Try to see what atoms are involved in the messages you get during NVT, before crashing. Have a look at those atoms, and atoms nearby, in the initial structure and try to resolve all overlaps or unfavourable bond configurations.
