compute cna/atom when using read restart

LAMMPS LAMMPS Mailing List Mirror
1

Dear lammps users,

When I tried to read a restart file and apply common neighbor analysis (compute cna/atom), it did not work properly, the return value is not showing any dislocation structure but some grid of the system (which i guess is related by parallel part).
Did any one meet this situation before and how should I fix it?

Thanks

The following is my input file and I am running the system in 128 processors.

units metal
dimension 3
boundary p p p
atom_style atomic

neighbor setting

neighbor 2 bin
neigh_modify delay 1

read_restart restart.700000.equil

EAM potentials

pair_style eam
pair_coeff * * /nfs/12/osu5634/software/lammps/lammps2/potentials/Ni_u3.eam

lattice fcc 3.52 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1

compute CNA all cnasurface/atom 3.00 6 8
compute CENTRO all centro/atom fcc
compute PE all pe/atom
compute KE all ke/atom
compute stress all stress/atom pair

compute new all temp/deform
velocity all create 1E-3 887723 temp new

timestep 0.001
thermo 200
thermo_modify temp new

fix 1 all nve

2

a) try converting the restart file to a data file, read it in,
    and see if you get the same answer - if so, it has nothing
    to do with restart files
b) this is your command:
    compute CNA all cnasurface/atom 3.00 6 8
    what is "cnasurface"? It's not the original LAMMPS
    compute - if it's your own compute, then you're on your own

Steve