#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run for 3-d periodic
#-------------------------------------------------------------------------------
units real
boundary p p p
atom_style full
pair_style lj/class2/coul/long 12
pair_modify mix sixthpower
pair_modify tail yes
bond_style class2
angle_style class2
dihedral_style class2
improper_style class2
special_bonds lj 0.0 0.0 1.0 coul 0.0 0.0 1.0
box tilt large
read_data structure.dat
Reading data file ...
orthogonal box = (0 0 0) to (38.8222 38.2464 46.9556)
4 by 2 by 4 MPI processor grid
reading atoms ...
3120 atoms
scanning bonds ...
1 = max bonds/atom
scanning angles ...
6 = max angles/atom
scanning impropers ...
4 = max impropers/atom
reading bonds ...
1824 bonds
reading angles ...
1680 angles
reading impropers ...
768 impropers
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 1
special bond factors coul: 0 0 1
4 = max # of 1-2 neighbors
3 = max # of 1-3 neighbors
4 = max # of special neighbors
special bonds CPU = 0.001 seconds
read_data CPU = 0.022 seconds
include parameters.dat
bond_coeff 1 1.5325 277.56 -433.25 2922.8 # o1S- s4o-
bond_coeff 2 0.97 563.28 -1428.2 1902.1 # h* o*
angle_coeff 1 109.47 87.637 -98.528 174.89 # o1S- s4o- o1S-
angle_coeff 2 103.7 49.84 -11.6 -8 # h* o* h*
dihedral_coeff 1 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 # xxx xxx xxx xxx
improper_coeff 1 0 0 # * s_ * *
angle_coeff 1 bb 130.45 1.5325 1.5325 # o1S- s4o- o1S-
angle_coeff 2 bb -9.5 0.97 0.97 # h* o* h*
angle_coeff 1 ba 0.0 0.0 1.5325 1.5325 # o1S- s4o- o1S-
angle_coeff 2 ba 22.35 22.35 0.97 0.97 # h* o* h*
dihedral_coeff 1 mbt 0.0000 0.0000 0.0000 0.0 # xxx xxx xxx xxx
dihedral_coeff 1 ebt 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0 0.0 # xxx xxx xxx xxx
dihedral_coeff 1 at 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0 0.0 # xxx xxx xxx xxx
dihedral_coeff 1 aat 0.0000 0.0 0.0 # xxx xxx xxx xxx
dihedral_coeff 1 bb13 0.0 0.0 0.0 # xxx xxx xxx xxx
improper_coeff 1 aa 0.0 0.0 0.0 120.0 120.0 120.0 # * s_ * *
pair_coeff 1 1 0.05899 5.1728 # Cl -
pair_coeff 2 2 0 1.098 # hw -
pair_coeff 3 3 0.4691 3.292 # K+ -
pair_coeff 4 4 0.2951 2.3959 # mg+ -
pair_coeff 5 5 0.148 3.7824 # o* -
pair_coeff 6 6 0.1436 3.6173 # o1S- -
pair_coeff 7 7 0.24 4.08 # s4o- -
neighbor 2.0 bin
neigh_modify delay 0 every 1 check yes
kspace_style pppm 0.00001
variable R equal 0.00198722
variable sysvol equal vol
variable sysmass equal mass(all)/6.0221367e+23
variable sysdensity equal v_sysmass/v_sysvol/1.0e-24
variable coulomb equal ecoul+elong
variable etotal equal etotal
variable pe equal pe
variable ke equal ke
variable evdwl equal evdwl
variable epair equal epair
variable ebond equal ebond
variable eangle equal eangle
variable edihed equal edihed
variable eimp equal eimp
variable lx equal lx
variable ly equal ly
variable lz equal lz
variable Nthermo equal 0
variable cella equal lx
variable cellb equal sqrt(ly*ly+xy*xy)
variable cellc equal sqrt(lz*lz+xz*xz+yz*yz)
variable cellalpha equal acos((xy*xz+ly*yz)/(v_cellb*v_cellc))
variable cellbeta equal acos(xz/v_cellc)
variable cellgamma equal acos(xy/v_cellb)
variable p equal press
variable pxx equal pxx
variable pyy equal pyy
variable pzz equal pzz
variable pyz equal pyz
variable pxz equal pxz
variable pxy equal pxy
variable sxx equal -pxx
variable syy equal -pyy
variable szz equal -pzz
variable syz equal -pyz
variable sxz equal -pxz
variable sxy equal -pxy
variable fmax equal fmax
variable fnorm equal fnorm
variable time equal step*dt+0.000001
variable surfacetension equal 0.5*v_lz*(0.5*(v_sxx+v_syy)-v_szz)
thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke
thermo_modify flush yes
#
# Set up the fixed and movable groups
#
group movable union all
3120 atoms in group movable
group fixed subtract all movable
0 atoms in group fixed
compute sum_f1 movable reduce sum fx fy fz
variable sum_fx equal c_sum_f1[1]
variable sum_fy equal c_sum_f1[2]
variable sum_fz equal c_sum_f1[3]
log 2.2_Minimization.out