Lammps exit without error message printed out

Dear Lammps Users,

I encountered an error that I could not identify any causes of it. Basically, I would like to simulate a system with spherical particles surrounded by water. The following is my input file. It exited automatically after creating the atoms without any error message printing out. I am pretty new to Lammps, and hope to get some ideas from you.

Here are the input script

Triclinic box water/sphere

units lj
boundary p p p
atom_style sphere

create geometry

lattice custom 0.85 a1 10.0 0.0 0.0 a2 0.0 10.0 0.0 a3 0.0 0.0 10.0 basis 0.0

0.0 0.0 basis 0.5 0.5 0.5
region box prism 0.0 10 0.0 10 0.0 10 0.0 0.0 0.0
create_box 2 box
create_atoms 1 box

set group all type/fraction 2 0.1 95392 # fraction 0.1
set type 1 mass 1.0
set type 2 mass 4.91 # calculated from dia=1.7

velocity all create 0.5 87287

Potential

pair_style hybrid/overlay lj/gromacs 1.91 2.5 lj/cut 1.12246 dpd/tstat 0.5 0.5 2.5

482748
pair_coeff 1 1 lj/gromacs 1.0 1.0 1.91 2.5
pair_coeff 2 2 lj/cut 14.19 1.7 1.12246
pair_modify shift yes
pair_coeff * * dpd/tstat 1.0 2.5

Integration

fix NVE_sys all nve/sphere
compute tempmainrun all temp/profile 0 1 0 z 5

neighbor and communicate

neighbor 1 multi
neigh_modify delay 0

Run and dump

communicate single vel yes
timestep 0.003
thermo 1
thermo_style custom step temp atoms pe ke etotal epair
dump 1 all atom 50 dump.flow

Your help will be greatly appreciated!
Best Regards,
Ling Qi

It exited because it came to the end of your input script.
There is no "run" command to run a simulation in
your script.

Why are you using atom_style sphere for finite size
particles? You are not using a pair style that requires
such particles, like a granular or colloidal pair style. Nor
are you creating any finite size particles by setting
their radii to non-zero. I think you should just use
atom_style atomic. The sigma setting in a potential
like lj sets the effective size of the particle.

Steve