Dear all,
I compute stress/atom with two free surfaces normal to z axis, x and y directions are periodic.
I think, if z is free, atoms are free to move in z direction to bring sigma_zz (stress/atom_zz) to zero. Since in continuum, for free surfaces, Cauchy stress is indeed zero everywhere, from surface to bulk.
I accept that sigma_xx and sigma_yy are not zero because surface stress, but sigma_zz should be zero everywhere.
I also realize if I sum over all blue numbers (see below) in dump file, I do get a value that is close to zero, but theory says it should be zero everywhere.
I am not quite sure what virial stress that LAMMPS calculates is, but I do get the knowledge from somewhere says virial stress is equivalent to Cauchy stress in continuum.
Could someone explain to me where I am wrong?
Best,
Weina
Input file
…………………………….
boundary p p s
lattice fcc 4.08 origin 0 0 0 orient x 1 -1 0 orient y 1 1 0 orient z 0 0 1
region mybox block 0 1 0 1 -1 1 units lattice
create_box 1 mybox
create_atoms 1 region mybox
compute 1 all stress/atom
dump 1 all custom 5000 001stress.dump.* id type x y z c_1[1] c_1[2] c_1[3]
minimize 1e-100 1e-120 10000000 100000000
…………………………….
Dump file
…………………………….
ITEM: ATOMS id type x y z c_1[1] c_1[2] c_1[3]
1 1 2.885 -7.22443e-17 -3.96282 1.5637e+06 1.5637e+06 195121
2 1 4.32749 1.4425 -2.05164 -251372 -251372 -262882
3 1 -5.85237e-17 -3.05374e-18 -3.96282 1.5637e+06 1.5637e+06 195121
4 1 2.9288e-17 2.885 -3.96282 1.5637e+06 1.5637e+06 195121
5 1 1.4425 1.4425 -2.05164 -251372 -251372 -262882
6 1 2.885 2.885 -3.96282 1.5637e+06 1.5637e+06 195121
7 1 4.32749 4.32749 -2.05164 -251372 -251372 -262882
8 1 1.4425 4.32749 -2.05164 -251372 -251372 -262882
9 1 2.885 -2.14746e-16 -2.397e-16 -9253.95 -9253.95 135522
10 1 4.32749 1.4425 2.05164 -251372 -251372 -262882
11 1 -1.97155e-16 6.30718e-17 -2.25604e-16 -9253.95 -9253.95 135522
12 1 -1.15288e-16 2.885 1.09957e-16 -9253.95 -9253.95 135522
13 1 1.4425 1.4425 2.05164 -251372 -251372 -262882
14 1 2.885 2.885 -1.79755e-16 -9253.95 -9253.95 135522
15 1 4.32749 4.32749 2.05164 -251372 -251372 -262882
16 1 1.4425 4.32749 2.05164 -251372 -251372 -262882
17 1 2.885 1.09667e-16 3.96282 1.5637e+06 1.5637e+06 195121
18 1 -1.56657e-16 3.50824e-17 3.96282 1.5637e+06 1.5637e+06 195121
19 1 -6.18136e-17 2.885 3.96282 1.5637e+06 1.5637e+06 195121
20 1 2.885 2.885 3.96282 1.5637e+06 1.5637e+06 195121
…………………………….