#------------------------------------------------------------------------------- # Stage 2.4: NPT integration for 10 ps with a timestep of 1 fs # Temperature 298.2 K # Pressure 1 atm #------------------------------------------------------------------------------- # Fix shake is needed for water in pcff+ fix shaken all shake .0001 20 50000 a 2 b 1 0 = # of size 2 clusters 0 = # of size 3 clusters 0 = # of size 4 clusters 3248 = # of frozen angles find clusters CPU = 0.002 seconds change_box all triclinic Changing box ... triclinic box = (0 0 0) to (33.78995 48.916 87.2) with tilt (0 0 0) kspace_style pppm 0.00001 reset_timestep 0 thermo_style custom step v_time press vol v_sysdensity temp ebond eangle edihed eimp evdwl ecoul etail elong pe ke thermo ${Nthermo} thermo 0 change_box all triclinic Changing box ... triclinic box = (0 0 0) to (33.78995 48.916 87.2) with tilt (0 0 0) fix 1 movable npt temp 298.2 298.2 100 z 1 1 100 drag 0 mtk yes nreset 200 fix 2 movable ave/time 1 999 1000 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_averages.txt off 1 fix 3 movable ave/time 1 1 1 v_time c_thermo_temp c_thermo_press v_sysvol v_sysdensity v_etotal v_cella v_cellb v_cellc v_cellalpha v_cellbeta v_cellgamma v_pe v_ke v_evdwl v_coulomb v_sxx v_syy v_szz v_syz v_sxz v_sxy file 2.4_instantaneous.txt fix 4 fixed setforce 0.0 0.0 0.0 restart 10000 2.4.restart dump sci all custom 1000 2.4.xyz id mol type q xs ys zs timestep 1 run 10000 PPPM initialization ... WARNING: System is not charge neutral, net charge = 0.0406 (../kspace.cpp:325) using 12-bit tables for long-range coulomb (../kspace.cpp:340) G vector (1/distance) = 0.25951072 grid = 30 40 60 stencil order = 5 estimated absolute RMS force accuracy = 0.0029008773 estimated relative force accuracy = 8.7359058e-06 using double precision KISS FFT 3d grid and FFT values/proc = 8228 2400 generated 0 of 36 mixed pair_coeff terms from sixthpower mixing rule Neighbor list info ... update every 1 steps, delay 0 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 14 ghost atom cutoff = 14 binsize = 7, bins = 5 7 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair lj/class2/coul/long, perpetual attributes: half, newton on pair build: half/bin/newton/tri stencil: half/bin/3d/tri bin: standard SHAKE stats (type/ave/delta/count) on step 0 1 1.00281 0.087494 9744 2 97.5782 10.5185 3248 Per MPI rank memory allocation (min/avg/max) = 13.21 | 13.68 | 14.16 Mbytes Step v_time Press Volume v_sysdensity Temp E_bond E_angle E_dihed E_impro E_vdwl E_coul E_tail E_long PotEng KinEng 0 1e-06 1001615.3 144130.19 1.3687485 420.74834 23.79694 153.05899 0 0 740557.84 -224051.2 -382.57099 -379053.56 137629.93 9747.4161 ERROR: Non-numeric pressure - simulation unstable (../fix_nh.cpp:1061) Last command: run 10000