forked from lijiext/lammps
82 lines
2.2 KiB
Plaintext
82 lines
2.2 KiB
Plaintext
# Pure Si crystal, qeq off, calculates stress-strain, no out.*.cfg
|
|
|
|
units metal
|
|
atom_style charge
|
|
dimension 3
|
|
boundary p p p
|
|
|
|
lattice diamond 5.431
|
|
region box block 0 4 0 4 0 4
|
|
create_box 1 box
|
|
create_atoms 1 box
|
|
mass 1 29.0
|
|
change_box triclinic
|
|
|
|
velocity all create 1.0 277387
|
|
|
|
pair_style comb
|
|
pair_coeff * * ffield.comb Si
|
|
|
|
neighbor 0.5 bin
|
|
neigh_modify every 20 delay 0 check no
|
|
|
|
fix 1 all box/relax aniso 0.0 vmax 0.001
|
|
timestep 0.00020
|
|
|
|
thermo_style custom step temp etotal pe evdwl ecoul &
|
|
#lx ly lz vol pxx pyy pzz
|
|
thermo_modify norm yes
|
|
thermo 1
|
|
minimize 1.0e-14 1.0e-20 1000 10000
|
|
min_modify dmax 0.2 line quadratic
|
|
|
|
unfix 1
|
|
run 1
|
|
|
|
### copy lines after this to any input script for elastic calculations ###
|
|
## Elastic constants calculations: strain box, measure box stress
|
|
## strain x, measure s_x, s_y, s_z, s_yz:
|
|
## calculates C11, C12, C13 and C14
|
|
|
|
fix 2 all deform 1 x scale 1.0001 remap x
|
|
compute perfx all stress/atom pair
|
|
compute fx all reduce sum &
|
|
c_perfx[1] c_perfx[2] c_perfx[3] &
|
|
c_perfx[4] c_perfx[5] c_perfx[6]
|
|
thermo_style custom step lx ly lz xy xz yz &
|
|
c_fx[1] c_fx[2] c_fx[3] c_fx[4] c_fx[5] c_fx[6]
|
|
run 10
|
|
|
|
## strain z, measure s_z: calculates C33
|
|
|
|
fix 2 all deform 1 z scale 1.0001 remap x
|
|
compute perfz all stress/atom pair
|
|
compute fz all reduce sum &
|
|
c_perfz[1] c_perfz[2] c_perfz[3] &
|
|
c_perfz[4] c_perfz[5] c_perfz[6]
|
|
thermo_style custom step lx ly lz xy xz yz &
|
|
c_fz[1] c_fz[2] c_fz[3] c_fz[4] c_fz[5] c_fz[6]
|
|
run 10
|
|
|
|
## strain yz, measure s_yz: calculates C44
|
|
|
|
fix 2 all deform 1 yz erate 0.0001 remap x
|
|
compute perfyz all stress/atom pair
|
|
compute fyz all reduce sum &
|
|
c_perfyz[1] c_perfyz[2] c_perfyz[3] &
|
|
c_perfyz[4] c_perfyz[5] c_perfyz[6]
|
|
thermo_style custom step lx ly lz xy xz yz &
|
|
c_fyz[1] c_fyz[2] c_fyz[3] c_fyz[4] c_fyz[5] c_fyz[6]
|
|
run 10
|
|
|
|
## strain xy, measure s_xy: calculates C66
|
|
|
|
fix 2 all deform 1 xy erate 0.0001 remap x
|
|
compute perfxy all stress/atom pair
|
|
compute fxy all reduce sum &
|
|
c_perfxy[1] c_perfxy[2] c_perfxy[3] &
|
|
c_perfxy[4] c_perfxy[5] c_perfxy[6]
|
|
thermo_style custom step lx ly lz xy xz yz &
|
|
c_fxy[1] c_fxy[2] c_fxy[3] c_fxy[4] c_fxy[5] c_fxy[6]
|
|
run 10
|