Dear lammpers,
I used the pair_style eam/fs to define the force between atom Fe (atom1 ) and Atom C(atom 2) in steel, but it always looks like there is no interacting force.
The pair_style commend is list below:
pair_style hybrid eam/fs lj/cut 10.0 # tersoff lcbop
pair_coeff * * eam/fs Fe-C.eam.fs Fe C C C #Fe-C
#pair_coeff * * lcbop C.lcbop NULL C NULL NULL #C-C
#pair_coeff 3 4 lj/cut 0.00759200 3.5 9.7625 #Fe-C
#pair_coeff 3 4 tersoff BNC.tersoff NULL NULL B N #B-N
#pair_coeff 1 1 lj/cut 0.105 3.85 10 #Fe-Fe
pair_coeff 2 4 lj/cut 0.004068 3.367 #N-C
pair_coeff 2 3 lj/cut 0.003311 3.411 #B-C
pair_coeff 1 3 lj/cut 0.00284 3.4 #Fe-B
pair_coeff 1 4 lj/cut 0.00759200 3.5 9.7625 #Fe-N
so what is wrong with my code?the file Fe-C.eam.fs is download from the web site:http://www.ctcms.nist.gov/potentials/Download/FeC-GJA/Fe-C_Hepburn_Ackland.eam.fs
College of Mechanical Engineering, Donghua University
Room 3034, academy building 4, 2999 North Renmin Road, Songjiang District
Shanghai, 201620, China
Tel: 86-159-0160-8286
Email: liuyao13@…8…