[lammps-users] (no subject)

Hi everyone,i try to create a binary alloy lattice and get its rdf, so i input as follows:

units metal
boundary p p p
atom_style atomic
lattice bcc 4.117
region box block 0 10 0 10 0 10
create_box 2 box
create_atoms 2 region box
pair_style meam
pair_coeff * * library.meam Cu Al cuzr.meam Cu Al
neighbor 2.0 bin
neigh_modify delay 5
group Cu type 1
compute 1 all temp
compute 2 all pressure 1
compute 3 all etotal/atom 1
fix 1 all npt 1400.0 1400.0 10.0 xyz 0.0 0.0 1000.0
fix 2 all rdf 1000 rdf.out 100 1 1 1 2 2 1 2 2
thermo 20
thermo_modify temp 1
timestep 0.001
dump 1 all atom 10000 dump.meam
reset_timestep 0
run 3000

when i run it ,the output of rdf is

TIMESTEP 0
r, g(1,1,r), ncoord(1,1,r), g(1,2,r), ncoord(1,2,r), g(2,1,r), ncoord(2,1,r), g(2,2,r), ncoord(2,2,r)
0.020000 nan nan nan nan nan nan 0.000000 0.000000
0.060000 nan nan nan nan nan nan 0.000000 0.000000
0.100000 nan nan nan nan nan nan 0.000000 0.000000
0.140000 nan nan nan nan nan nan 0.000000 0.000000
0.180000 nan nan nan nan nan nan 0.000000 0.000000
0.220000 nan nan nan nan nan nan 0.000000 0.000000
0.260000 nan nan nan nan nan nan 0.000000 0.000000
0.300000 nan nan nan nan nan nan 0.000000 0.000000
0.340000 nan nan nan nan nan nan 0.000000 0.000000
0.380000 nan nan nan nan nan nan 0.000000 0.000000
0.420000 nan nan nan nan nan nan 0.000000 0.000000
0.460000 nan nan nan nan nan nan 0.000000 0.000000
0.500000 nan nan nan nan nan nan 0.000000 0.000000
0.540000 nan nan nan nan nan nan 0.000000 0.000000
0.580000 nan nan nan nan nan nan 0.000000 0.000000
0.620000 nan nan nan nan nan nan 0.000000 0.000000
0.660000 nan nan nan nan nan nan 0.000000 0.000000
0.700000 nan nan nan nan nan nan 0.000000 0.000000
0.740000 nan nan nan nan nan nan 0.000000 0.000000
0.780000 nan nan nan nan nan nan 0.000000 0.000000
0.820000 nan nan nan nan nan nan 0.000000 0.000000

and in file log.lammps,it shows

0 atoms in group Cu

the first type atoms Cu weren’t put in a group. which parameter shall i change and how shall i do?
thanks for your time.

sincerely yours

syl

my email: [email protected]…414…

The "create_box 2 box" command will create a bcc lattice of atoms,
all of type 2. If you want an alloy lattice, you need to use other
options on the "lattice" command and define a basis set of 2 atoms.
Also the create_box command has options to specify which atom types
to assign to which basis atoms.

Steve