Dear all:

I am trying to simulate the water flow upon a solid plate, the model of water is tip4p and solid is cuu3. I arrange the water and cu by read_data command. But after runnig for a while, the system shows as follows:

Lammps (1 Oct 2006)

scanning data file��

2=max bond/atom

1=max angles/atom

Reading data file��

1 by 1 by 1 processor grid

528 atoms

320 bonds

160 angles

Finding 1-2 1-3 1-4 neighbors��

2��max # of 1-2 neighbors

1��max # of 1-3 neighbors

1��max # of 1-4 neighbors

2��max # of special neighbors

Segmentation Fault

Now I am so confused. How can I solve this problem?

The input script is as follows:

units real

dimension 2

boundary p s p

atom_style hybrid atomic full

read_data data.flownew

neighbor 2.0 bin

neigh_modify delay 5

bond_style harmonic

bond_coeff 1 19.514 0.9572

angle_style harmonic

angle_coeff 1 2.385 104.52

group water id<=300

group cu subtract all water

fix 1 water shake 0.0001 20 10 b 1 a 1 t 1 2

����