LAMMPS (7 Jan 2022)
KOKKOS mode is enabled (../kokkos.cpp:105)
will use up to 0 GPU(s) per node
log 2.1_Initialize.out
#-------------------------------------------------------------------------------
# Stage 2.1: Initialize LAMMPS run for 3-d periodic
#-------------------------------------------------------------------------------
units real
boundary p p p
atom_style full
box tilt large
read_data structure.dat
Reading data file ...
orthogonal box = (0 0 0) to (12.5214 12.5214 10.6652)
3 by 4 by 3 MPI processor grid
reading atoms ...
216 atoms
Finding 1-2 1-3 1-4 neighbors ...
special bond factors lj: 0 0 0
special bond factors coul: 0 0 0
0 = max # of 1-2 neighbors
0 = max # of 1-3 neighbors
0 = max # of 1-4 neighbors
1 = max # of special neighbors
special bonds CPU = 0.040 seconds
read_data CPU = 0.111 seconds
pair_style reax/c NULL checkqeq no safezone 10.0 mincap 1000
pair_coeff * * parameters.dat C O Si
variable R equal 0.00198722
variable sysvol equal vol
variable sysmass equal mass(all)/6.0221367e+023
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
216 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
#-------------------------------------------------------------------------------
# Stage 2.2: Minimization
#-------------------------------------------------------------------------------
min_style cg
min_modify dmax 0.05 line quadratic
reset_timestep 0
thermo_style custom step fmax fnorm press vol v_sysdensity v_sxx v_syy v_szz v_syz v_sxz v_sxy pe v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_sum_fx v_sum_fy v_sum_fz
WARNING: New thermo_style command, previous thermo_modify settings will be lost (../output.cpp:888)
dump sci all custom 1000 2.2.xyz id mol type q xs ys zs
thermo 100
fix 2_2_qeq all qeq/reax 1 0.0 10.0 1.0e-6 reax/c
minimize 0.0 1.0 1000 10000
CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE
Your simulation uses code contributions which should be cited:
- pair reaxff command:
- fix qeq/reaxff command:
The log file lists these citations in BibTeX format.
CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE
WARNING: Using 'neigh_modify every 1 delay 0 check yes' setting during minimization (../min.cpp:187)
Neighbor list info ...
update every 1 steps, delay 0 steps, check yes
max neighbors/atom: 2000, page size: 100000
master list distance cutoff = 12
ghost atom cutoff = 12
binsize = 6, bins = 3 3 2
2 neighbor lists, perpetual/occasional/extra = 2 0 0
(1) pair reax/c/kk, perpetual
attributes: half, newton off, ghost, kokkos_device
pair build: half/bin/ghost/kk/device
stencil: full/ghost/bin/3d
bin: kk/device
(2) fix qeq/reax/kk, perpetual, copy from (1)
attributes: half, newton off, ghost, kokkos_device
pair build: copy/kk/device
stencil: none
bin: none
Setting up cg/kk style minimization ...
Unit style : real
Current step : 0
Per MPI rank memory allocation (min/avg/max) = 23.35 | 23.44 | 23.59 Mbytes
Step Fmax Fnorm Press Volume v_sysdensity v_sxx v_syy v_szz v_syz v_sxz v_sxy PotEng v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_sum_fx v_sum_fy v_sum_fz
0 14459.142 24633.239 6078855.7 1672.1483 4.0817993 -4883986.3 -4834873.8 -8517707.1 18397.659 99448.595 269323.83 -11617.161 12.5214 12.5214 10.6652 1.5707963 1.5707963 1.5707963 6.9931243e-06 5.680964e-06 1.43384e-06