nb3b/harmonic parameters are in real units

This commit is contained in:
Axel Kohlmeyer
2020-06-29 13:24:13 -04:00
parent 3eda735d29
commit 5c08ce0818
2 changed files with 111 additions and 108 deletions

View File

@ -1,14 +1,17 @@
---
lammps_version: 5 May 2020
date_generated: Wed May 27 06:42:57 202
epsilon: 1.0e-12
lammps_version: 15 Jun 2020
date_generated: Mon Jun 29 13:21:51 202
epsilon: 1e-12
prerequisites: ! |
pair nb3b/harmonic
pre_commands: ! |
variable newton_pair delete
variable newton_pair index on
variable units delete
variable units index real
post_commands: ! |
change_box all x final 0 4.2 y final 0 4.2 z final 0 4.2
velocity all scale 100.0
input_file: in.manybody
pair_style: nb3b/harmonic
pair_coeff: ! |
@ -18,139 +21,139 @@ natoms: 64
init_vdwl: 2654.82155652424
init_coul: 0
init_stress: ! |-
-5.1017158429026651e+01 -1.3611903256981696e+01 6.4629061686009067e+01 -8.7135857294938191e+01 -1.1070515375917164e+02 7.5452837495497633e+01
-5.1017158429026708e+01 -1.3611903256981552e+01 6.4629061686009067e+01 -8.7135857294937580e+01 -1.1070515375917185e+02 7.5452837495497562e+01
init_forces: ! |2
1 -1.0877972789090633e+01 2.7818609977822151e+00 5.3286346737471801e-01
2 2.3844946347198697e+00 3.8120293270515626e+00 -1.7852915817427366e+01
3 -5.5060972496922012e+01 -1.4649767245154461e+01 -8.0050176965474122e+01
4 1.7259636030947714e+01 9.3366797180178907e+00 1.5880959275683665e+01
1 -1.0877972789090633e+01 2.7818609977822115e+00 5.3286346737471402e-01
2 2.3844946347198643e+00 3.8120293270515857e+00 -1.7852915817427341e+01
3 -5.5060972496921963e+01 -1.4649767245154461e+01 -8.0050176965474009e+01
4 1.7259636030947703e+01 9.3366797180178764e+00 1.5880959275683665e+01
5 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
6 1.7794992633465359e+01 -5.3419222655618320e+01 1.6321598833864002e+01
7 -6.7476300345434893e+01 6.2036681257356090e+00 -7.5139089756926083e+01
8 -2.4916802467563383e+02 3.9421473842157511e+00 -3.9134655098942254e+01
9 -4.2488750688310093e+00 4.0608847565742447e+00 -1.2332028245418186e+01
10 -5.6132393889097036e+00 1.7823276348533055e+00 1.2023878783019342e+00
11 2.7495846047096700e+02 3.6847059220182459e+01 -3.5747872835095954e+01
12 -8.5777111027672873e+01 -6.6703893761473649e+00 3.6593219992613264e+01
13 2.8150260504826656e+00 4.6881843084850182e+01 -2.8703190101594132e+00
14 4.3310482261409405e+01 5.7711396678423561e+01 5.5582267493210210e+01
6 1.7794992633465419e+01 -5.3419222655618618e+01 1.6321598833863909e+01
7 -6.7476300345423340e+01 6.2036681257305277e+00 -7.5139089756916974e+01
8 -2.4916802467563477e+02 3.9421473842150050e+00 -3.9134655098942197e+01
9 -4.2488750688310146e+00 4.0608847565742368e+00 -1.2332028245418179e+01
10 -5.6132393889096601e+00 1.7823276348532602e+00 1.2023878783018827e+00
11 2.7495846047096711e+02 3.6847059220182388e+01 -3.5747872835096082e+01
12 -8.5777111027672959e+01 -6.6703893761473418e+00 3.6593219992613299e+01
13 2.8150260504829605e+00 4.6881843084849855e+01 -2.8703190101594780e+00
14 4.3310482261409398e+01 5.7711396678423561e+01 5.5582267493210239e+01
15 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
16 1.2606972666826963e+00 -2.9949919839800152e+01 -3.8329654261193085e+01
17 -5.6603372959879501e+01 4.4199755708912090e+00 1.6646523859302970e+01
18 2.9145047807552785e+01 -3.8254292709284357e+01 -7.8685029062057863e+00
19 -3.9052132885118897e+02 -3.4874531498136278e+01 2.4335359036689007e+02
20 3.0523826465323211e+01 -2.7070245767455226e+01 5.8408159705563616e+00
16 1.2606972666821876e+00 -2.9949919839796344e+01 -3.8329654261188118e+01
17 -5.6603372959879628e+01 4.4199755708913484e+00 1.6646523859303024e+01
18 2.9145047807552785e+01 -3.8254292709284357e+01 -7.8685029062057987e+00
19 -3.9052132885118920e+02 -3.4874531498136022e+01 2.4335359036689030e+02
20 3.0523826465323282e+01 -2.7070245767455212e+01 5.8408159705563483e+00
21 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
22 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
23 7.4409325483379803e+01 -6.1338064797828309e+01 -4.0331673797102184e+01
23 7.4409325483379916e+01 -6.1338064797828324e+01 -4.0331673797102226e+01
24 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
25 -1.5434301681368330e+01 -1.8496773706788222e+01 -4.1692996127143651e+00
26 -7.2440757908713049e-01 -1.5411185566412183e+01 3.6729264463045794e+00
27 1.4857759801010695e+02 -8.7246598153787403e+01 3.5575180405196285e+00
28 1.1863298792067371e+02 8.2702860857456869e+01 1.5436298352562847e+00
29 4.7621783523283590e+01 2.8039369004096750e+00 2.7605524532590053e+01
30 2.8830156773978769e+02 -3.2784354352079447e+01 1.0018854518372478e+01
25 -1.5434301681368334e+01 -1.8496773706788261e+01 -4.1692996127144077e+00
26 -7.2440757908713937e-01 -1.5411185566412181e+01 3.6729264463045856e+00
27 1.4857759801010724e+02 -8.7246598153787659e+01 3.5575180405192732e+00
28 1.1863298792067371e+02 8.2702860857456841e+01 1.5436298352562812e+00
29 4.7621783523283653e+01 2.8039369004096772e+00 2.7605524532590049e+01
30 2.8830156773978769e+02 -3.2784354352079447e+01 1.0018854518372489e+01
31 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
32 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
33 -5.8839244877108614e+00 4.3985789859461818e+01 2.3192635019620489e+01
34 -8.7966629741582576e+00 -3.7772557274042633e+00 -5.0822108474305443e+00
35 -1.1045699552194318e+02 -2.5865892104728207e+01 -3.3353488725325462e+01
36 -3.9928053482362905e+01 -3.7650338080860749e+01 -8.6463197758339003e+00
33 -5.8839244877086143e+00 4.3985789859458791e+01 2.3192635019619907e+01
34 -8.7966629741582505e+00 -3.7772557274042629e+00 -5.0822108474305479e+00
35 -1.1045699552194318e+02 -2.5865892104728044e+01 -3.3353488725325391e+01
36 -3.9928053482362998e+01 -3.7650338080860742e+01 -8.6463197758338932e+00
37 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
38 3.7625199613228766e-01 -5.9409197479738708e+01 1.1096821440484145e+01
38 3.7625199613237004e-01 -5.9409197479738559e+01 1.1096821440484220e+01
39 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
40 -4.2079999725616096e+02 1.5030064155408525e+01 2.9502576524714826e+01
41 -7.1831348401232091e+00 2.3501863961013282e+01 3.3211487338838609e+01
42 -5.4333378031078610e+00 -3.8222885761439054e+00 5.2740780894998984e-01
40 -4.2079999725616096e+02 1.5030064155408574e+01 2.9502576524714936e+01
41 -7.1831348401232287e+00 2.3501863961013250e+01 3.3211487338838566e+01
42 -5.4333378031078574e+00 -3.8222885761439134e+00 5.2740780894998363e-01
43 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
44 6.3343303260433999e+01 4.9488759873343774e+01 5.0212324296502423e+01
45 4.2425954317559700e+02 5.6875319339115556e+01 -2.5998498159230604e+02
46 1.5374524942979843e+02 -2.4634003877741762e+01 -5.9674041391425958e-01
44 6.3343303260434006e+01 4.9488759873343746e+01 5.0212324296502395e+01
45 4.2425954317559712e+02 5.6875319339115613e+01 -2.5998498159230616e+02
46 1.5374524942979843e+02 -2.4634003877741691e+01 -5.9674041391425869e-01
47 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
48 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
49 2.4664457425766684e-01 8.5221588811519311e-03 7.1570435255873655e+00
50 -2.9771445332463109e+01 -1.3053973386895201e+00 9.5074220439575612e+00
51 -1.8696262048502865e+02 -3.0046250420275893e+01 -1.9739208373694105e+01
52 -2.3539906256917567e+02 8.1697006804398669e+01 -1.1372898509058719e+02
49 2.4664457425767483e-01 8.5221588811448257e-03 7.1570435255873699e+00
50 -2.9771445332463266e+01 -1.3053973386896178e+00 9.5074220439575896e+00
51 -1.8696262048502857e+02 -3.0046250420275893e+01 -1.9739208373694158e+01
52 -2.3539906256917558e+02 8.1697006804398626e+01 -1.1372898509058722e+02
53 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
54 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
55 -9.9888474306294608e+00 -6.3056149949578383e+01 -3.8789913572864435e+01
56 -1.5108125115180235e+02 8.7451440758083976e+01 6.2279587307098822e+01
57 -2.4707820063785416e+01 8.1430106606145536e+00 3.7927211164970971e+00
58 -7.9324636824924024e+00 1.1055979589866887e+00 -1.1979321661997941e+01
59 8.8647131566710783e+01 3.8727979707847247e+01 4.8159748997763351e+01
60 2.3037697796149806e+02 -4.6916609535445453e+00 9.6808251082403842e+01
61 8.6364330939177634e+01 1.4895642916188269e+01 9.8020771081533500e+00
62 4.1476164742574909e+01 -9.7738882325904193e+00 2.2124574239200477e+01
55 -9.9888474306294714e+00 -6.3056149949578362e+01 -3.8789913572864407e+01
56 -1.5108125115180263e+02 8.7451440758084161e+01 6.2279587307099106e+01
57 -2.4707820063785420e+01 8.1430106606145536e+00 3.7927211164970944e+00
58 -7.9324636824924042e+00 1.1055979589866993e+00 -1.1979321661997956e+01
59 8.8647131566696885e+01 3.8727979707853429e+01 4.8159748997752480e+01
60 2.3037697796149948e+02 -4.6916609535455187e+00 9.6808251082401242e+01
61 8.6364330939177592e+01 1.4895642916188271e+01 9.8020771081532168e+00
62 4.1476164742574909e+01 -9.7738882325904157e+00 2.2124574239200474e+01
63 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
64 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
run_vdwl: 2586.23310804073
run_vdwl: 2496.16387090594
run_coul: 0
run_stress: ! |-
-4.6723183908900717e+02 1.7852419684221977e+02 2.8870764224678692e+02 -7.3768673952427577e+01 8.1285715133964914e+01 -2.2175376931373922e+01
-3.1050965336040292e+01 -6.3663930994022884e+01 9.4714896330062999e+01 -1.1897654655567628e+02 -6.4811369677402881e+00 1.2347862730997569e+02
run_forces: ! |2
1 -1.1750670180272866e+00 1.3131514736699824e+01 -3.8826576290188388e+00
2 2.3357180987384143e+00 2.6172640820169175e+01 -1.3859108643114411e+00
3 6.0695590546373353e+01 2.0642136613920513e+01 5.4996351260125280e+00
4 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
1 -6.6214436114187141e+00 6.8003032240472616e-01 8.1610480878961345e-01
2 6.8952756210409225e+00 6.9800495995491607e+00 -2.1085653305920978e+01
3 -2.5923266064876653e+01 -7.2254055664041772e-01 -6.8249224259047452e+01
4 1.7514735996364870e+01 9.1382329672789950e+00 1.4609547038415045e+01
5 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
6 2.4546655526930735e+01 -3.1101905626646428e+01 3.8841426269193406e+01
7 -4.9616838878159811e+01 -1.0992414702422360e+01 4.8515728874677269e+00
8 -5.6453384170560540e+01 4.7672434067033826e+00 -1.9612102502967304e+00
9 -1.8392638626001174e+00 1.1828167793968992e+01 -2.0356506218224546e+01
10 -8.6013994200124220e+00 -1.0417698459793890e+01 5.4845710971785913e+00
11 6.0381617656756774e+01 -1.6620349247888100e+01 2.2225135196909434e+01
12 5.0979659960037345e+01 -2.2641974764565148e+01 8.6466127914636193e+00
6 1.8311337599848436e+01 -5.2737007463348519e+01 1.6786691829845285e+01
7 -7.5126152872375002e+01 3.9208090425206592e+00 -7.8501700045948425e+01
8 -1.7883231274220702e+02 1.2382156266165076e+01 -5.9370199178248981e+01
9 5.5713719759234035e+00 5.2160432479756871e+00 -7.2308807479301063e+00
10 -8.8969852437999064e-01 -1.4444413139727268e+00 1.1573923441213608e+01
11 1.4439507123610309e+02 -4.2519782552988204e+01 -3.4698021577236396e+01
12 -7.9891827746142440e+01 -1.3089531661279894e+01 3.3544238612989560e+01
13 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
14 9.6203506994347148e+00 2.0490558201590559e+01 3.8564561780893207e+00
14 1.5112474793527387e+01 3.7694869755527407e+01 3.7356461110106324e+01
15 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
16 4.2696341100597985e+00 -6.4394247758381127e+01 -2.3287529910455991e+00
17 -1.7809619943717419e+01 -1.9261779466882487e+00 2.6614972406661607e+01
18 -3.0062559069631483e+01 -4.6813849188661194e+01 2.0071390221802680e+01
19 2.2231227996316667e+00 5.5240594291246383e+00 -4.0911170180556759e+01
20 1.2156858565323886e+02 3.9830629488572413e+01 3.7706116674646424e+01
16 8.2461187714168827e-01 -3.3645319303810005e+01 -3.7521462215646530e+01
17 -6.1699650496788180e+01 1.2415735742737171e+01 1.3686739582320218e+01
18 2.8937858224261120e+01 -3.2763761650136210e+01 -2.1828920543524459e+00
19 7.5423787139990210e+01 1.5671749716356253e+01 -7.2118194981816799e+01
20 3.1376794755628680e+01 -3.1236998918287441e+01 7.5194815218724624e+00
21 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
22 -1.6364437371629357e+01 3.2268779082011257e+01 -2.5552249817970370e+01
23 -8.3597998911963472e+01 -1.7654897749235602e+01 7.1337191308940682e+00
22 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
23 6.8299963765568918e+01 -5.5707235162344126e+01 -3.8157782755725137e+01
24 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
25 1.8771431821434961e+00 3.7184795487210792e+01 -6.1577014025376728e+00
26 -2.3472845707924179e+00 -4.2517776291726975e+01 1.4877047560362083e+01
27 9.9311610633439727e+01 5.8342547785154601e+01 6.2490107493557900e+01
28 8.1997428592624331e+01 4.9305082539955812e+01 -3.1924492771504966e+01
29 1.5755014250094833e+01 1.6333818401948033e+01 2.0837131374062341e+01
30 4.6444961759604695e+01 2.9380895149140301e+01 -1.2549027721695549e+02
25 -1.4057244393584686e+01 -2.0289162180649946e+01 -6.8170959892909968e+00
26 -2.0187519714808841e+00 -1.7816926619477194e+01 3.6965530794144814e+00
27 1.3550009286410088e+02 -6.7838442433020745e+01 4.3133759785537102e+01
28 9.6888228678225971e+01 8.8823964822113297e+01 -5.0921310684595458e+00
29 5.0758136978827828e+01 -8.8983484466933565e+00 2.4899348744763654e+00
30 2.4524050422835563e+02 -1.3660139575889508e+01 -5.9330622488468725e+00
31 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
32 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
33 -5.5005500969351843e+01 -4.4888400262284573e+00 -7.9881718186484276e+00
34 1.2391358513429662e+01 4.2205224695342682e+00 -8.7307512727212266e+00
35 -6.8048459171953510e+01 -5.0746781744565169e+01 1.6313940341424583e+01
36 1.7865132435199360e+00 -4.6549184761244362e+01 -1.8272256554692827e+00
33 -5.9389712314058594e+00 5.1469485825476355e+01 1.5282542616314375e+01
34 -4.5946066700985710e+00 -4.4047708825320973e+00 -2.5896753139020205e+00
35 -7.6885241750207925e+01 -3.7858402496175358e+00 -1.0416208866746173e+01
36 -4.3251004252238630e+01 -3.0812732079587814e+01 1.5493207265310204e+01
37 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
38 2.9421387114598204e+01 -4.2086146965794974e+01 1.6800774864890816e+01
38 1.5022834606957978e+00 -6.0344797002288928e+01 9.7293521036630022e+00
39 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
40 -3.8166935284463143e+01 4.1676208829419494e+01 -7.1239426844135672e+00
41 8.9777269850365204e-01 4.4289814981844543e+01 2.5485173899838237e+01
42 1.1863601054760339e+01 -5.5963639046028719e+00 -2.6542953642884349e+01
43 -1.4318530885017411e+01 2.2057393766498969e+01 4.4372607097293105e+01
44 1.4461973377279693e+02 -4.1284603574668154e+01 3.8323143295598186e+01
45 1.2862655952927490e+01 5.8051830292127509e+00 5.6062109819492129e+01
46 3.9950844327128916e+00 1.0816632720372052e+01 2.7359382513314188e+01
40 -2.7902725378002401e+02 1.1520394360420082e+02 7.6029058865840060e+00
41 1.9875813960392859e+00 2.5898055298981063e+01 2.5255670110917048e+01
42 -1.2584608623742763e+00 -5.6180254661161735e+00 4.3355672941135257e+00
43 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
44 6.1017414877125326e+01 2.3847994346393680e+01 4.3262891796747915e+01
45 -5.7359830701947253e+01 3.4172133537980798e+01 8.2107294758634893e+01
46 1.1199817492468900e+02 -2.3757850899635525e+01 1.5384138774546539e+01
47 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
48 -8.2471941292632351e+00 7.2384082709877987e+01 7.6275073773084081e+00
49 7.0874669396804322e+00 1.0710819490993529e+01 6.8487083926871106e+00
50 -1.3702452391153543e+01 -3.1315676795828317e+01 -2.9390721004259010e+01
51 -5.1841862825969486e+01 -6.5229295766089322e+01 6.4345194792685618e+01
52 -1.0076730856559496e+02 -3.1323724404286470e+01 -2.7740629900930994e+01
48 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
49 -6.7826664559677354e+00 -1.8594048126803564e-01 5.8956046460290992e+00
50 -1.9821181472416256e+01 -1.8120466552580972e+00 4.4346376430197623e+00
51 -1.4173233779399769e+02 -4.2613815950203332e+01 -1.2166481368092634e+00
52 -1.8872790570220141e+02 8.0898507623818261e+01 -1.1989280050987685e+02
53 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
54 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
55 -1.1783270126548655e+02 -1.9472584883618183e-01 -3.6574242640823016e+01
56 -5.3306408717222062e+01 -4.8054342798707381e+01 -3.9204734285983008e+01
57 -1.1225984573690962e+01 2.0902821304606945e+01 -1.0842828823499091e+00
58 -1.5639600973783883e+01 3.7055729017065011e+00 -1.9828404272923528e+01
59 5.4139156830996429e+00 5.9778757589813964e+01 -2.3883800709664158e+01
60 4.0289832362530653e+01 -4.1876270741423269e+00 -1.8019718733793713e+01
61 -2.9803789778276474e+01 -1.2846624364406747e+00 -5.0425625243893258e+01
62 8.0660338957777533e-01 5.4686577469802673e+00 -1.3547627168270241e+01
63 -8.5708358255591914e+00 -3.0150178540572348e+01 -1.4967860711132145e+01
64 9.0239994663437773e-01 5.5410990098449542e-01 4.1571851677482030e+00
55 -6.3070794314782628e+00 -3.9493489491808347e+01 -3.5654394269934343e+01
56 -1.5536342282055168e+02 6.8759186611823679e+01 1.9484632984497214e+01
57 -1.6311185567740008e+01 5.5975547175367124e+00 -1.4425077069506711e+00
58 -6.4235602672067103e+00 2.6970192907549695e+00 -1.0923186469627634e+01
59 9.2546558740452369e+01 4.2666948660770792e+01 4.7533444080031025e+01
60 1.4559994425694987e+02 -3.5129424621390534e+01 1.1583916891031559e+02
61 5.9875382725562169e+01 6.8646156842044093e+00 2.3254650701763939e+00
62 3.9267471066687051e+01 -1.0670715066325148e+01 1.9913762076436925e+01
63 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
64 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00
...

View File

@ -189,7 +189,7 @@ TEST_F(PotentialFileReaderTest, GWZBL)
TEST_F(PotentialFileReaderTest, Nb3bHarmonic)
{
if (!verbose) ::testing::internal::CaptureStdout();
lmp->input->one("units metal");
lmp->input->one("units real");
PotentialFileReader reader(lmp, "MOH.nb3b.harmonic", "NB3B Harmonic");
if (!verbose) ::testing::internal::GetCapturedStdout();