LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] reg_unknown number of entries in pair/local

Re: [lammps-users] reg_unknown number of entries in pair/local

 From: santhosh mathesan Date: Tue, 22 May 2018 22:45:46 +0300

Thank you Axel.

Is there any other way in LAMMPS to determine the force on atom i due to atom j in the system?

On Tue, May 22, 2018 at 9:28 PM, Axel Kohlmeyer wrote:

On Tue, May 22, 2018 at 12:25 PM, santhosh mathesan wrote:
Dear all,

I am trying to obtain force acting on atom i due to atom j in a gold system during the deformation process. I did used a unit cell (simple system) to check whether pair/local calculates the forces for all the pairs within the system. I am not sure why the number of entries are more than the combination expected from the system.

​because minimum image conditions don't apply for LAMMPS. it creates real copies of the periodic images of your atoms up to the communication cutoff and all pairs (including interactions with periodic images) within the cutoff are considered.

axel.​

units metal
atom_style atomic
boundary p p p

lattice fcc 4.08 orient x 1 0 0  orient y 0 1 0 orient z 0 0 1 origin 0.05 0.05 0.05
region sim_region block 0 1 0 1 0 1 units lattice
create_box 1 sim_region
create_atoms 1 box

pair_style eam/alloy
pair_coeff * * Au_GLJ10_3.eam.alloy Au

dump 1 all atom 100 dump_min.atom
thermo 100
thermo_style custom step pe press vol etotal
min_style cg
minimize 1e-15 1e-7 70000 100000
print "MINIMIZATION DONE"
undump 1

write_data aftermin.txt

dump 2 all atom 1000 dump_for_deform.atom
variable tmp equal "lx"
variable L0 equal \${tmp}
variable strain equal "(lx - v_L0)/v_L0"
variable p1 equal "v_strain"
variable p2 equal "-pxx*0.0001"

variable srate equal 1.0e9
variable srate1 equal "v_srate/1.0e12"

fix 1 all nve
fix 2 all deform 1 x erate \${srate1} units box remap x
thermo_style custom step temp vol density pe ke etotal pxx pyy pzz lx ly lz
thermo          1000
timestep 0.001

compute 1 all property/local patom1 patom2
compute 2 all pair/local dist fx fy fz
dump 1 all local 1 tmp.dump index c_1[1] c_1[2] c_2[1] c_2[2] c_2[3] c_2[4]

run 1
undump 2
unfix 1
unfix 2

The output from dump file is

ITEM: TIMESTEP
1
ITEM: NUMBER OF ENTRIES
108
ITEM: BOX BOUNDS pp pp pp
0 4.08
0 4.08
0 4.08
ITEM: ENTRIES index c_1[1] c_1[2] c_2[1] c_2[2] c_2[3] c_2[4]
1 1 2 2.885 -0.0711631 -0.0711631 0
2 1 3 2.885 -0.0711631 0 -0.0711631
3 1 4 2.885 0 -0.0711631 -0.0711631
4 1 1 4.08 0.0463516 -0 -0
5 1 4 4.99696 0.0165693 0.00828464 0.00828464
6 1 1 5.76999 -0.000857057 0.000857057 -0
7 1 3 4.99696 -0.00828464 0.0165693 0.00828464
8 1 1 4.08 -0 0.0463516 -0
9 1 3 4.99696 0.00828464 0.0165693 0.00828464
10 1 1 5.76999 0.000857057 0.000857057 -0
11 1 2 4.99696 -0.00828464 -0.00828464 0.0165693
12 1 1 5.76999 -0 -0.000857057 0.000857057
13 1 2 4.99696 0.00828464 -0.00828464 0.0165693
14 1 1 5.76999 -0.000857057 -0 0.000857057
15 1 2 4.99696 -0.00828464 0.00828464 0.0165693
16 1 1 4.08 -0 -0 0.0463516
17 1 2 4.99696 0.00828464 0.00828464 0.0165693
18 1 1 5.76999 0.000857057 -0 0.000857057
19 1 1 5.76999 -0 0.000857057 0.000857057

I would like to know why there are more than " 1 2 " atom combination in the output.

Regards,
Santhosh

------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
lammps-users mailing list
lammps-users@...655....net
https://lists.sourceforge.net/lists/listinfo/lammps-users

--
Dr. Axel Kohlmeyer  akohlmey@...24...  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.

--
Sincerely,

Santhosh Mathesan
Postdoctoral Fellow
Nanomechanics Simulation Laboratory
Faculty of Mechanical Engineering
Technion Israel Institute of Technology
Ph.No. +972-559587456 (Israel)
Ph.No. +91-9003555294 (India)