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
����