Application aborted error

What is causing this error?

C:\Users\Sadat Abdullah\LAMMPS 64-bit 29Aug2024-MSMPI\bin> mpiexec -np 8 lmp -in in.relax
LAMMPS (29 Aug 2024 - Update 1)

job aborted:
[ranks] message

[0-2] terminated

[3] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 3

[4] terminated

[5] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 5

[6] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 6

[7] terminated

---- error analysis -----

[3,5-6] on LAPTOP-S968HDJK
lmp aborted the job. abort code 1

---- error analysis -----

This is the input script:

#  This is the control script for LAMMPS

echo 		both
#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run
#-------------------------------------------------------------------------------

units           metal
atom_style      charge
boundary        p p p
processors      * * * grid numa

read_data       Bi2Te3.data
replicate       8 8 4

#	group           watom type 1
#	group           seatom type 2
neighbor        2.0 bin
neigh_modify    delay 0 every 1 check yes

pair_style hybrid/overlay morse 30.0 coul/long 11.0
pair_coeff * * coul/long
pair_coeff 1 2 morse 0.975 1.285 3.089
pair_coeff 1 3 morse 0.582 1.257 3.251      
pair_coeff 2 2 morse 0.076 1.675 3.642      
pair_coeff 1 1 morse 0.085 2.212 4.203
pair_coeff 2 3 morse 0.807 0.731 4.497      
pair_coeff 3 3 morse 0.066 2.876 4.312      
kspace_style    pppm 1.0e-5


thermo_style    custom step temp pe ke etotal press vol
thermo          50
dump            1 all atom 100 dump.min
min_style       cg
minimize        1.0e-8 1.0e-8 5000 10000

undump          1

velocity        all create 10.00 156467 mom yes rot yes dist gaussian
fix             1 all npt temp 10.00 10.00 100 iso 0.0 0.0 1000
dump            2 all atom 50000 dumpheat.npt
timestep        0.001
run             100000
unfix           1
# heat from 10 to 100K
velocity        all create 100.00 156467 mom yes rot yes dist gaussian
fix             2a all  nve
run             20000
velocity        all create 100.00 156467 mom yes rot yes dist gaussian
run             20000
velocity        all create 100.00 156467 mom yes rot yes dist gaussian
run             20000
velocity        all create 100.00 156467 mom yes rot yes dist gaussian
run             20000
velocity        all create 100.00 156467 mom yes rot yes dist gaussian
run             20000
unfix           2a
fix             2b all npt temp 100.00 100.00 100 iso 0.0 0.0 1000
run             200000
unfix           2b

undump          2
# final relaxation run 
fix             5 all nve
dump            4 all atom 50000 dump300.nve
run             300000
unfix           5
write_restart   Bi2Te3.restart

Impossible to say from this information. Please try running without ‘mpiexec -np 8’ and report back what output you get.

If that is successful, then try to add the flags -echo screen -nb and report back what the output is.

Without the data file to try reproduce what you see on a different platform, there is not much to do.

C:\Users\Sadat Abdullah\LAMMPS 64-bit 29Aug2024-MSMPI\bin>lmp -in in.relax
LAMMPS (29 Aug 2024 - Update 1)
OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
  using 1 OpenMP thread(s) per MPI task
Loaded 0 plugins from D:\LAMMPS 64-bit 29Aug2024-MSMPI\plugins
#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run
#-------------------------------------------------------------------------------

units           metal
atom_style      charge
boundary        p p p
processors      * * * grid numa

read_data       Bi2Te3.data
Reading data file ...
  triclinic box = (0 0 0) to (4.4181304 3.8262132 32.257299) with tilt (-2.2090652 0 0)
WARNING: Triclinic box skew is large. LAMMPS will run inefficiently. (src/domain.cpp:221)
  1 by 1 by 1 MPI processor grid
  1 by 1 by 1 core grid within node
  reading atoms ...
  15 atoms
  read_data CPU = 0.002 seconds
replicate       8 8 4
Replication is creating a 8x8x4 = 256 times larger system...
  triclinic box = (0 0 0) to (35.345044 30.609706 129.0292) with tilt (-17.672522 0 0)
WARNING: Triclinic box skew is large. LAMMPS will run inefficiently. (src/domain.cpp:221)
  1 by 1 by 1 MPI processor grid
  1 by 1 by 1 core grid within node
  3840 atoms
  replicate CPU = 0.001 seconds

#       group           watom type 1
#       group           seatom type 2
neighbor        2.0 bin
neigh_modify    delay 0 every 1 check yes

pair_style hybrid/overlay morse 30.0 coul/long 11.0
pair_coeff * * coul/long
pair_coeff 1 2 morse 0.975 1.285 3.089
pair_coeff 1 3 morse 0.582 1.257 3.251
pair_coeff 2 2 morse 0.076 1.675 3.642
pair_coeff 1 1 morse 0.085 2.212 4.203
pair_coeff 2 3 morse 0.807 0.731 4.497
pair_coeff 3 3 morse 0.066 2.876 4.312
kspace_style    pppm 1.0e-5


thermo_style    custom step temp pe ke etotal press vol
thermo          50
dump            1 all atom 100 dump.min
min_style       cg
minimize        1.0e-8 1.0e-8 5000 10000

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

Your simulation uses code contributions which should be cited:
- Type Label Framework: https://doi.org/10.1021/acs.jpcb.3c08419
The log file lists these citations in BibTeX format.

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

PPPM initialization ...
  using 12-bit tables for long-range coulomb (src/kspace.cpp:342)
  G vector (1/distance) = 0.24879243
  grid = 20 27 54
  stencil order = 5
  estimated absolute RMS force accuracy = 0.00010939985
  estimated relative force accuracy = 7.5973993e-06
  using double precision FFTW3
  3d grid and FFT values/proc = 54162 29160
Generated 0 of 3 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update: every = 1 steps, delay = 0 steps, check = yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 32
  ghost atom cutoff = 32
  binsize = 16, bins = 4 2 9
  2 neighbor lists, perpetual/occasional/extra = 2 0 0
  (1) pair morse, perpetual
      attributes: half, newton on
      pair build: half/bin/atomonly/newton/tri
      stencil: half/bin/3d/tri
      bin: standard
  (2) pair coul/long, perpetual, trim from (1)
      attributes: half, newton on, cut 13
      pair build: trim
      stencil: none
      bin: none
Setting up cg style minimization ...
  Unit style    : metal
  Current step  : 0
ERROR on proc 0: Neighbor list overflow, boost neigh_modify one (src/npair_bin.cpp:248)
Last command: minimize        1.0e-8 1.0e-8 5000 10000

job aborted:
[ranks] message

[0] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 0

This is the data file

# LAMMPS data file written by OVITO Basic 3.10.6

15 atoms
3 atom types

0.0 4.41813045 xlo xhi
0.0 3.8262132069 ylo yhi
0.0 32.25729904 zlo zhi
-2.209065225 0.0 0.0 xy xz yz

Masses

1 208.98038  # Bi3+
2 127.6  # Te2-
3 127.6

Atoms  # charge

1 1 +0.38 0.0 0.0 12.8006565589
2 1 +0.38 0.0 0.0 19.4566424811
3 1 +0.38 2.209065225 1.2754044023 23.5530895722
4 1 +0.38 2.209065225 1.2754044023 30.2090754945
5 1 +0.38 4.4408920985e-16 2.5508088046 2.0482235455
6 1 +0.38 4.4408920985e-16 2.5508088046 8.7042094678
7 2 -0.26 0.0 0.0 6.9936321034
8 3 -0.24 0.0 0.0 25.2636669366
9 2 -0.26 2.209065225 1.2754044023 17.7460651167
10 2 -0.26 2.209065225 1.2754044023 3.75880091
11 3 -0.24 4.4408920985e-16 2.5508088046 28.49849813
12 2 -0.26 4.4408920985e-16 2.5508088046 14.5112339233
13 2 -0.26 0.0 0.0 0.0
14 3 -0.24 2.209065225 1.2754044023 10.7524330133
15 2 -0.26 4.4408920985e-16 2.5508088046 21.5048660267

C:\Users\Sadat Abdullah\LAMMPS 64-bit 29Aug2024-MSMPI\bin>lmp -in in.relax -echo screen -nb
LAMMPS (29 Aug 2024 - Update 1)
OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
  using 1 OpenMP thread(s) per MPI task
Loaded 0 plugins from D:\LAMMPS 64-bit 29Aug2024-MSMPI\plugins
#  This is the control script for LAMMPS

echo            both
#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run
#-------------------------------------------------------------------------------

units           metal
atom_style      charge
boundary        p p p
processors      * * * grid numa

read_data       Bi2Te3.data
Reading data file ...
  triclinic box = (0 0 0) to (4.4181304 3.8262132 32.257299) with tilt (-2.2090652 0 0)
WARNING: Triclinic box skew is large. LAMMPS will run inefficiently. (src/domain.cpp:221)
  1 by 1 by 1 MPI processor grid
  1 by 1 by 1 core grid within node
  reading atoms ...
  15 atoms
  read_data CPU = 0.001 seconds
replicate       8 8 4
Replication is creating a 8x8x4 = 256 times larger system...
  triclinic box = (0 0 0) to (35.345044 30.609706 129.0292) with tilt (-17.672522 0 0)
WARNING: Triclinic box skew is large. LAMMPS will run inefficiently. (src/domain.cpp:221)
  1 by 1 by 1 MPI processor grid
  1 by 1 by 1 core grid within node
  3840 atoms
  replicate CPU = 0.001 seconds

#       group           watom type 1
#       group           seatom type 2
neighbor        2.0 bin
neigh_modify    delay 0 every 1 check yes

pair_style hybrid/overlay morse 30.0 coul/long 11.0
pair_coeff * * coul/long
pair_coeff 1 2 morse 0.975 1.285 3.089
pair_coeff 1 3 morse 0.582 1.257 3.251
pair_coeff 2 2 morse 0.076 1.675 3.642
pair_coeff 1 1 morse 0.085 2.212 4.203
pair_coeff 2 3 morse 0.807 0.731 4.497
pair_coeff 3 3 morse 0.066 2.876 4.312
kspace_style    pppm 1.0e-5


thermo_style    custom step temp pe ke etotal press vol
thermo          50
dump            1 all atom 100 dump.min
min_style       cg
minimize        1.0e-8 1.0e-8 5000 10000

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

Your simulation uses code contributions which should be cited:
- Type Label Framework: https://doi.org/10.1021/acs.jpcb.3c08419
The log file lists these citations in BibTeX format.

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

PPPM initialization ...
  using 12-bit tables for long-range coulomb (src/kspace.cpp:342)
  G vector (1/distance) = 0.24879243
  grid = 20 27 54
  stencil order = 5
  estimated absolute RMS force accuracy = 0.00010939985
  estimated relative force accuracy = 7.5973993e-06
  using double precision FFTW3
  3d grid and FFT values/proc = 54162 29160
Generated 0 of 3 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update: every = 1 steps, delay = 0 steps, check = yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 32
  ghost atom cutoff = 32
  binsize = 16, bins = 4 2 9
  2 neighbor lists, perpetual/occasional/extra = 2 0 0
  (1) pair morse, perpetual
      attributes: half, newton on
      pair build: half/bin/atomonly/newton/tri
      stencil: half/bin/3d/tri
      bin: standard
  (2) pair coul/long, perpetual, trim from (1)
      attributes: half, newton on, cut 13
      pair build: trim
      stencil: none
      bin: none
Setting up cg style minimization ...
  Unit style    : metal
  Current step  : 0
ERROR on proc 0: Neighbor list overflow, boost neigh_modify one (src/npair_bin.cpp:248)
Last command: minimize        1.0e-8 1.0e-8 5000 10000

job aborted:
[ranks] message

[0] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 0

Here is the error that is terminating your run. The default neighbor list settings in LAMMPS are sufficient for the simulation of typical molecular systems are normal densities and with a typical interaction cutoff of 10-15 Angstrom. You are, however, using a cutoff of 30 Angstrom for the Morse potential and thus the number of neighbors is much larger (remember that the number of neighbors grows O(N^3) with the cutoff). So you either need to reduce the cutoff or increase the size of the neighbor list storage using neigh_modify.

1 Like