From b4acacf5cbc90d593a9455da1cd85e627415bde4 Mon Sep 17 00:00:00 2001 From: Axel Kohlmeyer Date: Tue, 14 Sep 2021 16:40:42 -0400 Subject: [PATCH] add minimal example and a unit test input --- examples/reaxff/in.reaxff.rdx-acks2 | 53 ++++++ .../reaxff/log.31Aug21.reaxff.rdx-acks2.g++.1 | 153 +++++++++++++++ .../reaxff/log.31Aug21.reaxff.rdx-acks2.g++.4 | 153 +++++++++++++++ .../tests/atomic-pair-reaxff-acks2.yaml | 176 ++++++++++++++++++ 4 files changed, 535 insertions(+) create mode 100644 examples/reaxff/in.reaxff.rdx-acks2 create mode 100644 examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.1 create mode 100644 examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.4 create mode 100644 unittest/force-styles/tests/atomic-pair-reaxff-acks2.yaml diff --git a/examples/reaxff/in.reaxff.rdx-acks2 b/examples/reaxff/in.reaxff.rdx-acks2 new file mode 100644 index 0000000000..541095081e --- /dev/null +++ b/examples/reaxff/in.reaxff.rdx-acks2 @@ -0,0 +1,53 @@ +# ReaxFF potential for RDX system + +units real + +atom_style charge +read_data data.rdx + +pair_style reaxff control.reax_c.rdx +pair_coeff * * ffield.reax C H O N + +compute reax all pair reaxff + +variable eb equal c_reax[1] +variable ea equal c_reax[2] +variable elp equal c_reax[3] +variable emol equal c_reax[4] +variable ev equal c_reax[5] +variable epen equal c_reax[6] +variable ecoa equal c_reax[7] +variable ehb equal c_reax[8] +variable et equal c_reax[9] +variable eco equal c_reax[10] +variable ew equal c_reax[11] +variable ep equal c_reax[12] +variable efi equal c_reax[13] +variable eqeq equal c_reax[14] + +neighbor 2.5 bin +neigh_modify every 10 delay 0 check no + +fix 1 all nve +fix 2 all acks2/reaxff 1 0.0 10.0 1.0e-6 reaxff + +variable nqeq equal f_2 + +thermo 10 +thermo_style custom step temp epair etotal press & + v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa & + v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq v_nqeq + +timestep 1.0 + +#dump 1 all atom 10 dump.reaxff.rdx + +#dump 2 all image 25 image.*.jpg type type & +# axes yes 0.8 0.02 view 60 -30 +#dump_modify 2 pad 3 + +#dump 3 all movie 25 movie.mpg type type & +# axes yes 0.8 0.02 view 60 -30 +#dump_modify 3 pad 3 + +run 100 diff --git a/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.1 b/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.1 new file mode 100644 index 0000000000..42d84167a3 --- /dev/null +++ b/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.1 @@ -0,0 +1,153 @@ +LAMMPS (31 Aug 2021) + using 1 OpenMP thread(s) per MPI task +# ReaxFF potential for RDX system + +units real + +atom_style charge +read_data data.rdx +Reading data file ... + orthogonal box = (35.000000 35.000000 35.000000) to (48.000000 48.000000 48.000000) + 1 by 1 by 1 MPI processor grid + reading atoms ... + 21 atoms + read_data CPU = 0.000 seconds + +pair_style reaxff control.reax_c.rdx +WARNING: Ignoring inactive control parameter: simulation_name (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: energy_update_freq (src/REAXFF/reaxff_control.cpp:97) +WARNING: Support for writing native trajectories has been removed after LAMMPS version 8 April 2021 (src/REAXFF/reaxff_control.cpp:113) +WARNING: Ignoring inactive control parameter: traj_title (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_info (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_forces (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_velocities (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: bond_info (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: angle_info (src/REAXFF/reaxff_control.cpp:97) +pair_coeff * * ffield.reax C H O N +Reading potential file ffield.reax with DATE: 2010-02-19 + +compute reax all pair reaxff + +variable eb equal c_reax[1] +variable ea equal c_reax[2] +variable elp equal c_reax[3] +variable emol equal c_reax[4] +variable ev equal c_reax[5] +variable epen equal c_reax[6] +variable ecoa equal c_reax[7] +variable ehb equal c_reax[8] +variable et equal c_reax[9] +variable eco equal c_reax[10] +variable ew equal c_reax[11] +variable ep equal c_reax[12] +variable efi equal c_reax[13] +variable eqeq equal c_reax[14] + +neighbor 2.5 bin +neigh_modify every 10 delay 0 check no + +fix 1 all nve +fix 2 all acks2/reaxff 1 0.0 10.0 1.0e-6 reaxff + +variable nqeq equal f_2 + +thermo 10 +thermo_style custom step temp epair etotal press v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq v_nqeq + +timestep 1.0 + +#dump 1 all atom 10 dump.reaxff.rdx + +#dump 2 all image 25 image.*.jpg type type # axes yes 0.8 0.02 view 60 -30 +#dump_modify 2 pad 3 + +#dump 3 all movie 25 movie.mpg type type # axes yes 0.8 0.02 view 60 -30 +#dump_modify 3 pad 3 + +run 100 + +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: + +@Article{Aktulga12, + author = {H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y. Grama}, + title = {Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques}, + journal = {Parallel Computing}, + year = 2012, + volume = 38, + pages = {245--259} +} + +- fix acks2/reaxff command: + +@Article{O'Hearn2020, + author = {K. A. O'Hearn, A. Alperen, and H. M. Aktulga}, + title = {Fast Solvers for Charge Distribution Models on Shared Memory Platforms}, + journal = {SIAM J. Sci. Comput.}, + year = 2020, + volume = 42, + pages = {1--22} +} + +CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE + +Neighbor list info ... + update every 10 steps, delay 0 steps, check no + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 12.5 + ghost atom cutoff = 12.5 + binsize = 6.25, bins = 3 3 3 + 2 neighbor lists, perpetual/occasional/extra = 2 0 0 + (1) pair reaxff, perpetual + attributes: half, newton off, ghost + pair build: half/bin/newtoff/ghost + stencil: full/ghost/bin/3d + bin: standard + (2) fix acks2/reaxff, perpetual, copy from (1) + attributes: half, newton off, ghost + pair build: copy + stencil: none + bin: none +Per MPI rank memory allocation (min/avg/max) = 15.42 | 15.42 | 15.42 Mbytes +Step Temp E_pair TotEng Press v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq v_nqeq + 0 0 -1808.3917 -1808.3917 28736.41 -2958.4712 79.527715 0.31082031 0 98.589783 25.846176 -0.18034154 0 16.709078 -9.1620736 938.43732 -2.8201336e-06 0 0.00098778735 13 + 10 1304.7509 -1914.0924 -1836.3081 -21153.752 -2726.219 -17.446777 0.19121437 0 58.532415 2.5619735 -77.337316 0 14.042901 -5.8126196 837.38787 -6.2534642e-06 0 0.0069050076 8 + 20 427.93813 -1859.4968 -1833.9848 -23712.924 -2777.9198 1.1772924 0.029117162 0 88.192077 0.1543307 -53.152357 0 28.462512 -9.6520395 863.2094 -1.1899783e-06 0 0.0027021013 9.5 + 30 600.85197 -1867.8807 -1832.0602 -15720.955 -2765.2737 1.4937773 0.06185725 0 80.013841 0.090367845 -47.085245 0 16.406666 -11.121903 857.53144 2.0736158e-06 0 0.0022173663 11 + 40 1108.9095 -1898.9259 -1832.817 -38935.482 -2778.2032 2.0997015 0.14353694 0 78.798725 0.26421692 -58.91856 0 14.943876 -10.029332 851.97505 4.5556609e-07 0 1.3535608e-05 13.5 + 50 1155.316 -1890.9414 -1822.0658 7674.2712 -2804.003 8.5334231 0.15421192 0 78.75008 2.2698852 -77.319055 0 31.922369 -8.9637444 877.71732 6.9388694e-05 0 -0.0028728853 10 + 60 1237.1418 -1895.0813 -1821.3276 854.31642 -2722.5859 -14.2286 0.22149986 0 44.37587 9.5571735 -73.477237 0 18.566827 -4.8577223 847.34568 1.6885606e-05 0 0.0011149828 12 + 70 1025.7868 -1882.2967 -1821.1431 25212.869 -2890.8508 54.451033 0.29038763 0 112.30176 2.4666006 -64.1052 0 19.268525 -11.520433 895.4047 5.9278993e-06 0 -0.0032206344 9.5 + 80 619.97124 -1857.6805 -1820.7202 -12086.828 -2740.4855 -5.2006679 0.12319993 0 99.775674 2.7800075 -77.812022 0 25.298145 -7.562536 845.40484 3.766017e-05 0 -0.0017469014 9.5 + 90 1169.6095 -1893.5597 -1823.832 -6196.2424 -2816.571 34.143182 0.24873152 0 91.015691 0.81251232 -69.865835 0 17.93908 -10.502528 859.2142 1.3276274e-05 0 0.006324974 9.5 + 100 1032.79 -1882.6603 -1821.0893 6839.1293 -2811.285 8.1859322 0.22109056 0 69.13601 8.7302413 -74.363952 0 22.220038 -7.7705968 902.27101 6.5173803e-05 0 -0.005154992 7.5 +Loop time of 0.118719 on 1 procs for 100 steps with 21 atoms + +Performance: 72.777 ns/day, 0.330 hours/ns, 842.328 timesteps/s +99.8% CPU use with 1 MPI tasks x 1 OpenMP threads + +MPI task timing breakdown: +Section | min time | avg time | max time |%varavg| %total +--------------------------------------------------------------- +Pair | 0.09767 | 0.09767 | 0.09767 | 0.0 | 82.27 +Neigh | 0.0055032 | 0.0055032 | 0.0055032 | 0.0 | 4.64 +Comm | 0.00021768 | 0.00021768 | 0.00021768 | 0.0 | 0.18 +Output | 0.00023556 | 0.00023556 | 0.00023556 | 0.0 | 0.20 +Modify | 0.014979 | 0.014979 | 0.014979 | 0.0 | 12.62 +Other | | 0.0001125 | | | 0.09 + +Nlocal: 21.0000 ave 21 max 21 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Nghost: 546.000 ave 546 max 546 min +Histogram: 1 0 0 0 0 0 0 0 0 0 +Neighs: 1102.00 ave 1102 max 1102 min +Histogram: 1 0 0 0 0 0 0 0 0 0 + +Total # of neighbors = 1102 +Ave neighs/atom = 52.476190 +Neighbor list builds = 10 +Dangerous builds not checked +Total wall time: 0:00:00 diff --git a/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.4 b/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.4 new file mode 100644 index 0000000000..74f4038548 --- /dev/null +++ b/examples/reaxff/log.31Aug21.reaxff.rdx-acks2.g++.4 @@ -0,0 +1,153 @@ +LAMMPS (31 Aug 2021) + using 1 OpenMP thread(s) per MPI task +# ReaxFF potential for RDX system + +units real + +atom_style charge +read_data data.rdx +Reading data file ... + orthogonal box = (35.000000 35.000000 35.000000) to (48.000000 48.000000 48.000000) + 1 by 2 by 2 MPI processor grid + reading atoms ... + 21 atoms + read_data CPU = 0.000 seconds + +pair_style reaxff control.reax_c.rdx +WARNING: Ignoring inactive control parameter: simulation_name (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: energy_update_freq (src/REAXFF/reaxff_control.cpp:97) +WARNING: Support for writing native trajectories has been removed after LAMMPS version 8 April 2021 (src/REAXFF/reaxff_control.cpp:113) +WARNING: Ignoring inactive control parameter: traj_title (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_info (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_forces (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: atom_velocities (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: bond_info (src/REAXFF/reaxff_control.cpp:97) +WARNING: Ignoring inactive control parameter: angle_info (src/REAXFF/reaxff_control.cpp:97) +pair_coeff * * ffield.reax C H O N +Reading potential file ffield.reax with DATE: 2010-02-19 + +compute reax all pair reaxff + +variable eb equal c_reax[1] +variable ea equal c_reax[2] +variable elp equal c_reax[3] +variable emol equal c_reax[4] +variable ev equal c_reax[5] +variable epen equal c_reax[6] +variable ecoa equal c_reax[7] +variable ehb equal c_reax[8] +variable et equal c_reax[9] +variable eco equal c_reax[10] +variable ew equal c_reax[11] +variable ep equal c_reax[12] +variable efi equal c_reax[13] +variable eqeq equal c_reax[14] + +neighbor 2.5 bin +neigh_modify every 10 delay 0 check no + +fix 1 all nve +fix 2 all acks2/reaxff 1 0.0 10.0 1.0e-6 reaxff + +variable nqeq equal f_2 + +thermo 10 +thermo_style custom step temp epair etotal press v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq v_nqeq + +timestep 1.0 + +#dump 1 all atom 10 dump.reaxff.rdx + +#dump 2 all image 25 image.*.jpg type type # axes yes 0.8 0.02 view 60 -30 +#dump_modify 2 pad 3 + +#dump 3 all movie 25 movie.mpg type type # axes yes 0.8 0.02 view 60 -30 +#dump_modify 3 pad 3 + +run 100 + +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: + +@Article{Aktulga12, + author = {H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y. Grama}, + title = {Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques}, + journal = {Parallel Computing}, + year = 2012, + volume = 38, + pages = {245--259} +} + +- fix acks2/reaxff command: + +@Article{O'Hearn2020, + author = {K. A. O'Hearn, A. Alperen, and H. M. Aktulga}, + title = {Fast Solvers for Charge Distribution Models on Shared Memory Platforms}, + journal = {SIAM J. Sci. Comput.}, + year = 2020, + volume = 42, + pages = {1--22} +} + +CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE + +Neighbor list info ... + update every 10 steps, delay 0 steps, check no + max neighbors/atom: 2000, page size: 100000 + master list distance cutoff = 12.5 + ghost atom cutoff = 12.5 + binsize = 6.25, bins = 3 3 3 + 2 neighbor lists, perpetual/occasional/extra = 2 0 0 + (1) pair reaxff, perpetual + attributes: half, newton off, ghost + pair build: half/bin/newtoff/ghost + stencil: full/ghost/bin/3d + bin: standard + (2) fix acks2/reaxff, perpetual, copy from (1) + attributes: half, newton off, ghost + pair build: copy + stencil: none + bin: none +Per MPI rank memory allocation (min/avg/max) = 11.41 | 12.58 | 13.79 Mbytes +Step Temp E_pair TotEng Press v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq v_nqeq + 0 0 -1808.3956 -1808.3956 28736.41 -2958.4712 79.527715 0.31082031 0 98.589783 25.846176 -0.18034154 0 16.709078 -9.1620736 938.43732 -5.6770858e-06 0 -0.0028392419 14 + 10 1304.7509 -1914.1049 -1836.3206 -21153.752 -2726.219 -17.446777 0.19121437 0 58.532415 2.5619735 -77.337316 0 14.042901 -5.8126196 837.38787 5.0978181e-06 0 -0.005556783 10 + 20 427.93812 -1859.4876 -1833.9756 -23712.923 -2777.9198 1.1772926 0.029117161 0 88.192077 0.1543307 -53.152357 0 28.462513 -9.6520395 863.2094 2.1862491e-05 0 0.011844766 6 + 30 600.85195 -1867.8815 -1832.061 -15720.949 -2765.2737 1.4937783 0.061857248 0 80.013842 0.090367838 -47.085244 0 16.406666 -11.121903 857.53144 -1.2004549e-06 0 0.0014350102 11 + 40 1108.9094 -1898.8905 -1832.7815 -38935.484 -2778.2032 2.0996992 0.14353693 0 78.798728 0.26421691 -58.91856 0 14.943875 -10.029332 851.97505 4.1455385e-06 0 0.03548355 8 + 50 1155.3161 -1890.9404 -1822.0649 7674.2794 -2804.0031 8.5334252 0.15421192 0 78.75008 2.2698849 -77.319055 0 31.92237 -8.9637447 877.71732 3.7570219e-05 0 -0.0018876928 7 + 60 1237.1418 -1895.0785 -1821.3248 854.31669 -2722.5859 -14.228606 0.22149986 0 44.375869 9.5571747 -73.477236 0 18.566825 -4.8577216 847.34567 -6.0652199e-06 0 0.0039746806 9 + 70 1025.7867 -1882.2715 -1821.118 25212.883 -2890.8508 54.451049 0.29038764 0 112.30178 2.4666008 -64.105192 0 19.268524 -11.520433 895.40472 7.5942987e-05 0 0.021807003 5.5 + 80 619.97151 -1857.6791 -1820.7188 -12086.851 -2740.4854 -5.2006737 0.12319996 0 99.775675 2.7800093 -77.812024 0 25.29814 -7.5625349 845.40485 -3.2701271e-06 0 -0.00028718671 11 + 90 1169.6099 -1893.5491 -1823.8214 -6196.2462 -2816.571 34.143181 0.24873153 0 91.015695 0.81251269 -69.865839 0 17.939079 -10.502527 859.2142 2.297465e-05 0 0.016894027 6.5 + 100 1032.7901 -1882.6521 -1821.081 6839.0982 -2811.2849 8.1859195 0.22109056 0 69.136035 8.730243 -74.363951 0 22.220028 -7.7705949 902.27102 6.4245187e-06 0 0.0030953783 10 +Loop time of 0.140636 on 4 procs for 100 steps with 21 atoms + +Performance: 61.435 ns/day, 0.391 hours/ns, 711.054 timesteps/s +95.2% CPU use with 4 MPI tasks x 1 OpenMP threads + +MPI task timing breakdown: +Section | min time | avg time | max time |%varavg| %total +--------------------------------------------------------------- +Pair | 0.042596 | 0.057261 | 0.073915 | 4.7 | 40.72 +Neigh | 0.0024769 | 0.003423 | 0.0044286 | 1.2 | 2.43 +Comm | 0.0078962 | 0.025491 | 0.040466 | 7.3 | 18.13 +Output | 0.00030279 | 0.00066626 | 0.0017195 | 0.0 | 0.47 +Modify | 0.05257 | 0.0536 | 0.054564 | 0.3 | 38.11 +Other | | 0.0001949 | | | 0.14 + +Nlocal: 5.25000 ave 14 max 0 min +Histogram: 1 0 2 0 0 0 0 0 0 1 +Nghost: 355.500 ave 433 max 282 min +Histogram: 1 0 0 0 2 0 0 0 0 1 +Neighs: 303.000 ave 784 max 0 min +Histogram: 1 0 1 1 0 0 0 0 0 1 + +Total # of neighbors = 1212 +Ave neighs/atom = 57.714286 +Neighbor list builds = 10 +Dangerous builds not checked +Total wall time: 0:00:00 diff --git a/unittest/force-styles/tests/atomic-pair-reaxff-acks2.yaml b/unittest/force-styles/tests/atomic-pair-reaxff-acks2.yaml new file mode 100644 index 0000000000..9fc6e54858 --- /dev/null +++ b/unittest/force-styles/tests/atomic-pair-reaxff-acks2.yaml @@ -0,0 +1,176 @@ +--- +lammps_version: 31 Aug 2021 +date_generated: Tue Sep 14 16:33:13 2021 +epsilon: 2e-10 +skip_tests: omp +prerequisites: ! | + pair reaxff + fix acks2/reaxff +pre_commands: ! | + echo screen + variable newton_pair delete + variable newton_pair index on + atom_modify map array + units real + atom_style charge + lattice diamond 3.77 + region box block 0 2 0 2 0 2 + create_box 3 box + create_atoms 1 box + displace_atoms all random 0.1 0.1 0.1 623426 + mass 1 1.0 + mass 2 12.0 + mass 3 16.0 + set type 1 type/fraction 2 0.5 998877 + set type 2 type/fraction 3 0.5 887766 + set type 1 charge 0.00 + set type 2 charge 0.01 + set type 3 charge -0.01 + velocity all create 100 4534624 loop geom +post_commands: ! | + fix qeq all acks2/reaxff 1 0.0 8.0 1.0e-20 reaxff +input_file: in.empty +pair_style: reaxff NULL checkqeq yes +pair_coeff: ! | + * * ffield.reax.mattsson H C O +extract: ! "" +natoms: 64 +init_vdwl: -3296.3503506624793 +init_coul: -46582.941410876396 +init_stress: ! |- + -4.0889791539749363e+03 -5.3386746000293797e+03 -4.0808531437530951e+03 8.4337212661084749e+01 1.2848022290076187e+02 2.1468330620396631e+02 +init_forces: ! |2 + 1 -1.1929105250740845e+02 -1.4227215649327911e+02 2.1690003483475269e+02 + 2 -1.3331393160584241e+02 -2.4072728673874801e+02 -2.5617184179287540e+02 + 3 -1.2277758113579102e+02 3.0537119398406725e+02 5.2041769044942647e+01 + 4 4.3936404332100686e+01 -5.5511409696656926e+01 1.1097037460922759e+02 + 5 1.7864066365499838e+02 6.4034214571442291e+01 -1.0759738771465670e+02 + 6 3.5224339213063365e+02 -3.4719844127558417e+02 2.6425030914356228e+02 + 7 -2.0955545981628606e+02 -1.5858086260185297e+02 -5.7110672086393009e+01 + 8 -1.7044115435176795e+02 -5.3258138982765971e+01 1.2346092050629277e+02 + 9 -4.3715452075923416e+01 -3.8277423788683727e+01 6.3239179557158693e+01 + 10 9.3379874863494152e+01 5.9621377499868792e+01 9.2692428302505945e+01 + 11 -2.9198145425949332e+02 3.4467679177423884e+02 4.2632171085782431e+02 + 12 -2.6868067466378778e+02 -6.8644694781243447e+02 3.6910245078857537e+02 + 13 7.2123761164950440e+02 5.3620750239907272e+02 -3.8334401482414739e+02 + 14 8.4304627105934017e+01 -7.2466826898455409e+01 -8.7680132417144804e+01 + 15 -8.6225588839900468e+01 1.2225376793640189e+02 3.1577248706704783e+01 + 16 9.5393618380056623e+01 9.2886438384296099e+01 -7.3037388977053901e+00 + 17 1.3836060168005312e+02 -2.7473812676174262e+01 -9.4295524163369763e+01 + 18 -1.0165081413937811e+02 -1.1015152146781125e+02 5.1410985169799659e+01 + 19 1.5206789621670023e+02 2.5283948687709136e+02 -1.3894819646210772e+02 + 20 -1.1131767660064460e+02 -3.4867570304165632e+01 -2.2712139843531563e+01 + 21 -9.3180727628491212e+01 -7.7585741763401899e+01 -1.0055132802024193e+02 + 22 1.5388786773142479e+02 -1.2742120442767296e+02 5.9418228396501462e+01 + 23 1.4353594998513574e+02 1.9241855678142144e+01 -8.1940929642015490e+01 + 24 -1.2522859809648111e+02 3.4415263825860251e+02 -1.1666304527016791e+02 + 25 -6.1721877485216332e+01 2.1180528773074695e+02 -7.5185905744284099e+01 + 26 2.9047739438739399e+02 -2.0857900661856581e+02 2.4684212082797956e+02 + 27 1.9541979561640971e+01 -1.1873091999325840e+01 2.3080404132665154e+02 + 28 1.4779678412475869e+02 1.5072902909656133e+02 3.0786573270627889e+02 + 29 -5.7906765009047632e+01 5.2017884794099068e+01 6.2397336175147302e+00 + 30 -5.0895158598354335e+01 -2.5651607581450094e+02 -2.3443114587348992e+02 + 31 -1.5901188482871649e+02 1.7515333469043051e+02 -1.3486355088745134e+02 + 32 -1.3360496387491185e+02 -1.6026008719190622e+02 -6.3974117051860412e+02 + 33 -2.5865357552595057e+01 -1.3593003904804139e+01 3.6914416827045400e+01 + 34 -5.5345811893104553e+01 7.9787476336364392e+01 -1.1478329621445167e+02 + 35 6.2660444245581886e+01 -3.0334076804276236e+02 2.3657913727895405e+01 + 36 5.6628731130607559e+02 6.3924255616951359e+02 -6.4337876635896555e+02 + 37 -3.2625798296558662e+02 -3.4996432828486286e+02 4.6462264485478067e+02 + 38 -1.7064714626235613e+02 5.4874767598319613e+02 -5.7882903400837307e+01 + 39 2.2473038304436571e+01 -6.2855043226973365e+01 6.4935240565245138e+00 + 40 -2.0126181967178340e+02 1.9610382824263672e+02 4.3026826548397452e+01 + 41 -1.3922156023283114e+02 6.4923943342784469e+01 -4.9530732454693187e+01 + 42 2.5995311603257232e+02 6.2417650092087072e+01 -9.7757113311038410e+01 + 43 2.9060357406136120e+02 1.4859563620459687e+02 1.2902634448684583e+02 + 44 7.2822452489664343e+01 -9.5648852111255877e+01 -8.7605260909388264e+01 + 45 -7.6599938131930799e+01 2.7335310025430903e+02 1.4131077433184396e+02 + 46 -2.7864489804675839e+02 -3.3973466156303232e+02 -7.8954068435531696e+01 + 47 -5.5399367185908162e+01 -1.7022266452251142e+02 -3.6356135964982343e+01 + 48 4.5552031101078306e+01 -1.9615434812741569e+02 1.3353869041480218e+01 + 49 -2.6504355894047563e+02 3.8696727838168464e+01 1.6916938099889404e+02 + 50 4.3639038114698417e+02 -3.5529826932462339e+02 2.2030242322174141e+02 + 51 2.5201585479728806e+01 -5.5088430398826756e+02 2.9465735389668356e+02 + 52 -9.0610418024348732e+01 7.6451159625662612e+01 -1.0407170461264690e+02 + 53 -7.3230804863728633e+02 4.9918863718008072e+02 -3.7838006824966345e+02 + 54 -4.7813763806167337e+01 1.4821304034781713e+02 1.0972324440296121e+02 + 55 9.5154177757014878e+01 -1.6691953935970690e+02 6.0033792912059681e+01 + 56 1.3771822987048927e+02 -2.7596944143424753e+02 -4.7071002706580586e+01 + 57 7.7412666503504113e+01 -4.7147589918747506e+01 -7.4398816462734942e+01 + 58 -6.6332114011965132e+01 1.0893333226816128e+02 1.4604164580109380e+00 + 59 -1.4727436711601135e+02 1.3495995435479597e+02 1.3790547616349306e+02 + 60 4.4129860556328424e+01 8.8352155376666673e+01 -1.4754960139395871e+02 + 61 1.4664747834652729e+00 -1.6970260796882511e+02 -9.4741886350021787e+01 + 62 4.4268359315177207e+01 4.4244449542813292e+01 1.2469468703392482e+02 + 63 1.7715967912563977e+02 -8.9272317987559674e+01 -1.4306028518778581e+02 + 64 4.5068920113622809e+01 1.1297321948286621e+02 6.4572008812611145e+01 +run_vdwl: -3296.3469885397826 +run_coul: -823492.7373897503 +run_stress: ! |- + -1.8585481565254432e+05 -1.6813801222973305e+05 -1.4502245083109321e+05 -5.8156385069933705e+04 -2.9152676301583182e+04 4.2004485688893306e+04 +run_forces: ! |2 + 1 -3.1538860418527776e+02 -2.1667070105808438e+02 2.4327387701805449e+01 + 2 1.6469401642403250e+03 -2.2638495114953002e+03 8.1885571019561985e+03 + 3 -9.5233405445336075e+02 -1.3724167062887382e+03 5.1002346300387226e+03 + 4 -1.9086024577541921e+04 2.6665514538832704e+04 -1.5472292269485462e+04 + 5 7.8578011077546921e+02 1.1543491971265603e+03 -2.9601286588036862e+02 + 6 -4.1699480975660663e+03 6.3051224921793291e+03 -3.3124016041183961e+04 + 7 -1.7016494722270763e+03 -1.6238373161421100e+03 -3.1498701157264741e+03 + 8 -8.7621187948259321e+03 -5.8937249360767862e+02 1.0215462943849338e+03 + 9 2.6988978990627977e+04 -7.5719531974766342e+03 2.4402547340884375e+04 + 10 -5.5005889662574373e+02 2.8664451881339041e+03 3.6505720022004916e+03 + 11 -6.7136518028847252e+02 1.6951905338011699e+01 1.7613233778368631e+01 + 12 -4.1634897368756356e+02 -8.5427900773396550e+02 1.8957750739417386e+03 + 13 -2.4216201154307335e+04 -1.5001197373926550e+04 -1.3022247607934969e+04 + 14 -1.0405638578539543e+04 2.0054631376981383e+03 -1.7796805319319501e+04 + 15 1.1107583611850200e+03 1.3946144565444863e+04 2.5796292910159968e+03 + 16 7.0815040550639987e+03 -1.4994388002155864e+02 5.3941167060157622e+03 + 17 -1.0906710352768299e+02 -3.7797021439898720e+02 -7.1043712362237557e+02 + 18 3.1335355909249570e+03 -1.2979911718656356e+04 4.0028269316329097e+03 + 19 3.0219327888397208e+02 1.3132901589285169e+02 -2.2586289218389567e+02 + 20 -9.9495344251336473e+02 -1.0100323982979886e+04 7.2449256267058790e+03 + 21 -4.1626605209794234e+03 1.8207407400795291e+04 -1.9931886311902479e+04 + 22 -1.3809401079344286e+02 -1.3699694825434080e+03 4.9000482399450652e+02 + 23 -8.4097849417170528e+03 -6.7421987083314089e+03 2.1165500301472144e+04 + 24 6.6770947596674314e+03 7.0334161079873766e+03 -4.2395467771461354e+03 + 25 -3.0766262766235326e+02 4.5208369485829849e+01 -2.4215550554801229e+02 + 26 4.3282796665842488e+03 3.9725663253263883e+03 -1.5104557062597873e+03 + 27 3.4305125101384992e+03 -2.0757974258609988e+02 2.2419823347029745e+04 + 28 -1.2042796161721953e+03 2.8601773568967531e+02 -3.7347872364058077e+02 + 29 -1.4221968667948236e+03 -1.6884513031798073e+02 1.4718645011377384e+03 + 30 5.6087854923317345e+03 -1.3929352331590022e+04 8.9229119079296124e+03 + 31 1.8012567809349953e+03 1.7963340895216925e+03 -6.0693227081363002e+03 + 32 -2.4212994769744628e+02 -1.4663787679151267e+02 -2.6941110807653297e+02 + 33 -8.0205637538401052e+03 -6.7997183454421474e+03 7.0360939390233079e+03 + 34 -5.6825090133159119e+03 -4.9410619181051680e+02 -2.5378668650973555e+03 + 35 2.8749496957898991e+04 -4.6373022885745434e+02 -6.1483933048064773e+03 + 36 3.4096922298128607e+02 1.6757417075507853e+03 -8.8894223458424244e+02 + 37 -3.7343729694303511e+02 4.0656193398515529e+02 -1.8738927422246070e+02 + 38 1.4904299646376625e+04 -1.2064340896314498e+04 2.6539478453790794e+04 + 39 -7.0532532206553833e+03 -2.3957404157713672e+03 -1.1447343311718747e+04 + 40 6.1842420953100736e+03 -3.6659635754801893e+03 2.1054957749221794e+03 + 41 -8.1984840065986391e+03 -1.0365942326176664e+04 1.3952656987833076e+04 + 42 5.9577914748549119e+02 -1.7775021048376038e+02 -1.2021290039743464e+04 + 43 -1.9747460399160684e+03 1.6750061311272725e+03 2.8167237669329461e+03 + 44 -1.4678732439406349e+02 -2.4848519095612715e+04 4.7216919473597045e+02 + 45 4.5832362832616300e+03 1.2867485447504534e+04 1.1425823083324687e+04 + 46 2.1681087046653184e+02 -4.4934132470777809e+02 2.3087937120670163e+03 + 47 -5.2078060850766360e+02 -1.0208437402520989e+03 -2.2033873154011167e+03 + 48 2.6765788351065014e+04 -4.6344248657866183e+03 -2.9668864467185216e+01 + 49 -3.3609820362098151e+04 7.0960668414334923e+03 2.0836704420875976e+03 + 50 5.2223928486543973e+03 -2.9199713953576970e+03 4.6245497550815480e+03 + 51 5.5712197650885059e+01 -1.4394835567396435e+03 -1.5189574303433574e+01 + 52 6.0412308047232736e+03 -4.5508173248590192e+03 -2.8441095710784916e+03 + 53 2.9534447696061790e+02 -1.3661813387871148e+03 -4.0034798093494584e+02 + 54 -6.1569823475168255e+02 -4.7361494916845299e+02 -2.9928489978624157e+02 + 55 7.4978541085184744e+01 -1.5875365627631319e+02 1.2957750158518593e+02 + 56 -1.0429337846817589e+04 1.1766963241565751e+04 -4.9827771650037830e+03 + 57 1.1099912931100564e+04 -5.5053980025296833e+03 -1.1685512724787503e+04 + 58 -1.9431048462192582e+03 5.6326023332911245e+03 -2.6172065262884835e+03 + 59 -5.3914957888616245e+03 6.1370263012888418e+03 9.9937202300412118e+03 + 60 -1.5282247767913028e+04 5.2869249261205459e+03 -8.0462797163142568e+03 + 61 8.2234507406447065e+03 -6.6972733225279626e+00 -4.5676623487624911e+03 + 62 -1.2440498519086455e+03 1.5259685674863772e+03 -6.4071229464774854e+02 + 63 5.2477591133845465e+02 -1.1707580235553110e+03 -5.2912221162139249e+03 + 64 1.1950180636484571e+04 2.2135788612431359e+04 -8.1931421390373671e+03 +...