QStandardPaths: XDG_RUNTIME_DIR not set, defaulting to '/tmp/runtime-bwaters2' [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [I] RUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:18:08 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 4 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 13 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 1 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=1872 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.041666666666667 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.16666666666667 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 4.1666666667e-02 -0.0000000000e+00 -1.6666666667e-01 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=3.797819e+01 pxpb=0.000000e+00 pypb=5.818505e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 581 need remove atom 584 need remove atom 613 need remove atom 617 need remove atom 620 need remove atom 628 need remove atom 649 need remove atom 653 need remove atom 656 need remove atom 664 need remove atom 685 need remove atom 689 need remove atom 692 need remove atom 700 need remove atom 721 need remove atom 725 need remove atom 728 need remove atom 736 need remove atom 757 need remove atom 761 need remove atom 764 need remove atom 772 need remove atom 793 need remove atom 797 need remove atom 800 need remove atom 808 need remove atom 1077 need remove atom 1106 need remove atom 1113 need remove atom 1142 need remove atom 1149 need remove atom 1178 need remove atom 1185 need remove atom 1214 need remove atom 1221 need remove atom 1250 need remove atom 1257 need remove atom 1286 need remove atom 1293 [I] need removenum=39 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2265165180e-01 9.9871069543e-02 -3.5054641616e-01 dub= -3.3848569130e-01 9.9873514585e-02 3.5734546097e-01 duc= 3.2265255180e-01 -1.0120529100e-01 -3.5054753635e-01 [I] overall tilt ux= -6.6113734310e-01 2.4450420311e-06 7.0789187713e-01 uy= 8.9999954123e-07 -2.0107636054e-01 -1.1201866604e-06 [I] storedr[0]= -3.2903322664e-02 -6.4856338045e-04 -5.4342057636e-02 [I] originally 1872 atoms [I] insert 0 atoms [I] now 1872 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 1872 atoms [I] remove 39 atoms [I] now 1833 atoms [I] NP=1833 n=1833 [I] ASSIGN finalcnfile = dipole_1.0.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.0.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.0.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:18:08 2023 CPU time spent: 0.157843 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:18:44 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 5 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 16 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 1 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=2880 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.033333333333333 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.16666666666667 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 3.3333333333e-02 -0.0000000000e+00 -1.6666666667e-01 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=4.747274e+01 pxpb=0.000000e+00 pypb=7.161237e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 1298 need remove atom 1317 need remove atom 1320 need remove atom 1334 need remove atom 1343 need remove atom 1345 need remove atom 1353 need remove atom 1356 need remove atom 1364 need remove atom 1370 need remove atom 1379 need remove atom 1381 need remove atom 1389 need remove atom 1392 need remove atom 1400 need remove atom 1406 need remove atom 1415 need remove atom 1417 need remove atom 1425 need remove atom 1428 need remove atom 1436 need remove atom 1442 need remove atom 1451 need remove atom 1453 need remove atom 1461 need remove atom 1464 need remove atom 1472 need remove atom 1478 need remove atom 1487 need remove atom 1489 need remove atom 1497 need remove atom 1500 need remove atom 1508 need remove atom 1514 need remove atom 1523 need remove atom 1525 need remove atom 1533 need remove atom 1536 need remove atom 1544 need remove atom 1550 need remove atom 1559 need remove atom 1561 need remove atom 1569 need remove atom 1572 need remove atom 1580 need remove atom 1595 need remove atom 1597 need remove atom 1616 [I] need removenum=48 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2141155831e-01 1.0030255369e-01 -3.4818189825e-01 dub= -3.3711888039e-01 1.0030485216e-01 3.5463137594e-01 duc= 3.2141246976e-01 -1.0160575309e-01 -3.4818301481e-01 [I] overall tilt ux= -6.5853043870e-01 2.2984762590e-06 7.0281327419e-01 uy= 9.1144464465e-07 -2.0190830677e-01 -1.1165637032e-06 [I] storedr[0]= -2.8377296483e-02 -6.5222106525e-04 -4.9845081202e-02 [I] originally 2880 atoms [I] insert 0 atoms [I] now 2880 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 2880 atoms [I] remove 48 atoms [I] now 2832 atoms [I] NP=2832 n=2832 [I] ASSIGN finalcnfile = dipole_1.3.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.3.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.3.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:18:44 2023 CPU time spent: 0.247128 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:19:19 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 5 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 16 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 2 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=5760 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.033333333333333 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.083333333333333 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 3.3333333333e-02 -0.0000000000e+00 -8.3333333333e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=4.747274e+01 pxpb=0.000000e+00 pypb=7.161237e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 2594 need remove atom 2613 need remove atom 2616 need remove atom 2630 need remove atom 2649 need remove atom 2652 need remove atom 2666 need remove atom 2675 need remove atom 2677 need remove atom 2685 need remove atom 2688 need remove atom 2696 need remove atom 2702 need remove atom 2711 need remove atom 2713 need remove atom 2721 need remove atom 2724 need remove atom 2732 need remove atom 2738 need remove atom 2747 need remove atom 2749 need remove atom 2757 need remove atom 2760 need remove atom 2768 need remove atom 2774 need remove atom 2783 need remove atom 2785 need remove atom 2793 need remove atom 2796 need remove atom 2804 need remove atom 2810 need remove atom 2819 need remove atom 2821 need remove atom 2829 need remove atom 2832 need remove atom 2840 need remove atom 2846 need remove atom 2855 need remove atom 2857 need remove atom 2865 need remove atom 2868 need remove atom 2876 need remove atom 2882 need remove atom 2891 need remove atom 2893 need remove atom 2901 need remove atom 2904 need remove atom 2912 need remove atom 2918 need remove atom 2927 need remove atom 2929 need remove atom 2937 need remove atom 2940 need remove atom 2948 need remove atom 2954 need remove atom 2963 need remove atom 2965 need remove atom 2973 need remove atom 2976 need remove atom 2984 need remove atom 2990 need remove atom 2999 need remove atom 3001 need remove atom 3009 need remove atom 3012 need remove atom 3020 need remove atom 3026 need remove atom 3035 need remove atom 3037 need remove atom 3045 need remove atom 3048 need remove atom 3056 need remove atom 3062 need remove atom 3071 need remove atom 3073 need remove atom 3081 need remove atom 3084 need remove atom 3092 need remove atom 3098 need remove atom 3107 need remove atom 3109 need remove atom 3117 need remove atom 3120 need remove atom 3128 need remove atom 3134 need remove atom 3143 need remove atom 3145 need remove atom 3153 need remove atom 3156 need remove atom 3164 need remove atom 3179 need remove atom 3181 need remove atom 3200 need remove atom 3215 need remove atom 3217 need remove atom 3236 [I] need removenum=96 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2141155831e-01 1.0030255369e-01 -3.4818189825e-01 dub= -3.3711888039e-01 1.0030485216e-01 3.5463137594e-01 duc= 3.2141246976e-01 -1.0160575309e-01 -3.4818301481e-01 [I] overall tilt ux= -6.5853043870e-01 2.2984762590e-06 7.0281327419e-01 uy= 9.1144464465e-07 -2.0190830677e-01 -1.1165637032e-06 [I] storedr[0]= -2.8377296483e-02 -6.5222106525e-04 -4.9845081202e-02 [I] originally 5760 atoms [I] insert 0 atoms [I] now 5760 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 5760 atoms [I] remove 96 atoms [I] now 5664 atoms [I] NP=5664 n=5664 [I] ASSIGN finalcnfile = dipole_1.5.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.5.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.5.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:19:19 2023 CPU time spent: 0.483336 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:19:58 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 6 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 19 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 2 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=8208 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.027777777777778 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.083333333333333 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 2.7777777778e-02 -0.0000000000e+00 -8.3333333333e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=5.696729e+01 pxpb=0.000000e+00 pypb=8.503969e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 3097 need remove atom 3101 need remove atom 3104 need remove atom 3112 need remove atom 3133 need remove atom 3137 need remove atom 3140 need remove atom 3148 need remove atom 3169 need remove atom 3173 need remove atom 3176 need remove atom 3184 need remove atom 3205 need remove atom 3209 need remove atom 3212 need remove atom 3220 need remove atom 3241 need remove atom 3245 need remove atom 3248 need remove atom 3256 need remove atom 3277 need remove atom 3281 need remove atom 3284 need remove atom 3292 need remove atom 3313 need remove atom 3317 need remove atom 3320 need remove atom 3328 need remove atom 3349 need remove atom 3353 need remove atom 3356 need remove atom 3364 need remove atom 3385 need remove atom 3389 need remove atom 3392 need remove atom 3400 need remove atom 3421 need remove atom 3425 need remove atom 3428 need remove atom 3436 need remove atom 3457 need remove atom 3461 need remove atom 3464 need remove atom 3472 need remove atom 3493 need remove atom 3497 need remove atom 3500 need remove atom 3508 need remove atom 3529 need remove atom 3533 need remove atom 3536 need remove atom 3544 need remove atom 3565 need remove atom 3569 need remove atom 3572 need remove atom 3580 need remove atom 3601 need remove atom 3605 need remove atom 3608 need remove atom 3616 need remove atom 3637 need remove atom 3641 need remove atom 3644 need remove atom 3652 need remove atom 3673 need remove atom 3677 need remove atom 3680 need remove atom 3688 need remove atom 3709 need remove atom 3713 need remove atom 3716 need remove atom 3724 need remove atom 3745 need remove atom 3760 need remove atom 3781 need remove atom 3796 need remove atom 4490 need remove atom 4497 need remove atom 4526 need remove atom 4533 need remove atom 4562 need remove atom 4569 need remove atom 4598 need remove atom 4605 need remove atom 4634 need remove atom 4641 need remove atom 4670 need remove atom 4677 need remove atom 4706 need remove atom 4713 need remove atom 4742 need remove atom 4749 need remove atom 4778 need remove atom 4785 need remove atom 4814 need remove atom 4821 need remove atom 4850 need remove atom 4857 need remove atom 4886 need remove atom 4893 need remove atom 4922 need remove atom 4929 need remove atom 4958 need remove atom 4965 need remove atom 4994 need remove atom 5001 need remove atom 5030 need remove atom 5037 need remove atom 5066 need remove atom 5073 need remove atom 5102 need remove atom 5109 need remove atom 5138 need remove atom 5174 [I] need removenum=114 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2056269483e-01 1.0059376853e-01 -3.4657730662e-01 dub= -3.3617953780e-01 1.0059596403e-01 3.5278858480e-01 duc= 3.2056361385e-01 -1.0187597663e-01 -3.4657842049e-01 [I] overall tilt ux= -6.5674223264e-01 2.1955042347e-06 6.9936589142e-01 uy= 9.1901394494e-07 -2.0246974515e-01 -1.1138667277e-06 [I] storedr[0]= -2.4894744268e-02 -6.4776761058e-04 -4.4614980217e-02 [I] originally 8208 atoms [I] insert 0 atoms [I] now 8208 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 8208 atoms [I] remove 114 atoms [I] now 8094 atoms [I] NP=8094 n=8094 [I] ASSIGN finalcnfile = dipole_1.6.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.6.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.6.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:19:59 2023 CPU time spent: 0.688097 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:21:13 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 7 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 22 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 2 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=11088 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.023809523809524 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.083333333333333 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 2.3809523810e-02 -0.0000000000e+00 -8.3333333333e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=6.646184e+01 pxpb=0.000000e+00 pypb=9.846701e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 5186 need remove atom 5195 need remove atom 5197 need remove atom 5205 need remove atom 5208 need remove atom 5216 need remove atom 5222 need remove atom 5231 need remove atom 5233 need remove atom 5241 need remove atom 5244 need remove atom 5252 need remove atom 5258 need remove atom 5267 need remove atom 5269 need remove atom 5277 need remove atom 5280 need remove atom 5288 need remove atom 5294 need remove atom 5303 need remove atom 5305 need remove atom 5313 need remove atom 5316 need remove atom 5324 need remove atom 5330 need remove atom 5339 need remove atom 5341 need remove atom 5349 need remove atom 5352 need remove atom 5360 need remove atom 5366 need remove atom 5375 need remove atom 5377 need remove atom 5385 need remove atom 5388 need remove atom 5396 need remove atom 5402 need remove atom 5411 need remove atom 5413 need remove atom 5421 need remove atom 5424 need remove atom 5432 need remove atom 5438 need remove atom 5447 need remove atom 5449 need remove atom 5457 need remove atom 5460 need remove atom 5468 need remove atom 5474 need remove atom 5483 need remove atom 5485 need remove atom 5493 need remove atom 5496 need remove atom 5504 need remove atom 5510 need remove atom 5519 need remove atom 5521 need remove atom 5529 need remove atom 5532 need remove atom 5540 need remove atom 5546 need remove atom 5555 need remove atom 5557 need remove atom 5565 need remove atom 5568 need remove atom 5576 need remove atom 5582 need remove atom 5591 need remove atom 5593 need remove atom 5601 need remove atom 5604 need remove atom 5612 need remove atom 5618 need remove atom 5627 need remove atom 5629 need remove atom 5637 need remove atom 5640 need remove atom 5648 need remove atom 5654 need remove atom 5663 need remove atom 5665 need remove atom 5673 need remove atom 5676 need remove atom 5684 need remove atom 5690 need remove atom 5699 need remove atom 5701 need remove atom 5709 need remove atom 5712 need remove atom 5720 need remove atom 5726 need remove atom 5735 need remove atom 5737 need remove atom 5745 need remove atom 5748 need remove atom 5756 need remove atom 5762 need remove atom 5771 need remove atom 5773 need remove atom 5781 need remove atom 5784 need remove atom 5792 need remove atom 5798 need remove atom 5807 need remove atom 5809 need remove atom 5817 need remove atom 5820 need remove atom 5828 need remove atom 5834 need remove atom 5843 need remove atom 5845 need remove atom 5853 need remove atom 5856 need remove atom 5864 need remove atom 5870 need remove atom 5879 need remove atom 5881 need remove atom 5889 need remove atom 5892 need remove atom 5900 need remove atom 5906 need remove atom 5915 need remove atom 5917 need remove atom 5925 need remove atom 5928 need remove atom 5936 need remove atom 5942 need remove atom 5951 need remove atom 5953 need remove atom 5961 need remove atom 5964 need remove atom 5972 [I] need removenum=132 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.1994518979e-01 1.0080347901e-01 -3.4541704755e-01 dub= -3.3549435150e-01 1.0080559832e-01 3.5145562881e-01 duc= 3.1994611417e-01 -1.0207054037e-01 -3.4541815935e-01 [I] overall tilt ux= -6.5543954129e-01 2.1193100660e-06 6.9687267637e-01 uy= 9.2438448379e-07 -2.0287401938e-01 -1.1117977068e-06 [I] storedr[0]= -2.2301937125e-02 -6.4194480749e-04 -3.9867607879e-02 [I] originally 11088 atoms [I] insert 0 atoms [I] now 11088 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 11088 atoms [I] remove 132 atoms [I] now 10956 atoms [I] NP=10956 n=10956 [I] ASSIGN finalcnfile = dipole_1.9.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.9.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.9.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:21:14 2023 CPU time spent: 0.922977 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:22:36 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 9 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 29 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 2 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=18792 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.018518518518519 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.083333333333333 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 1.8518518519e-02 -0.0000000000e+00 -8.3333333333e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=8.545094e+01 pxpb=0.000000e+00 pypb=1.297974e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 8875 need remove atom 8877 need remove atom 8880 need remove atom 8911 need remove atom 8913 need remove atom 8916 need remove atom 8939 need remove atom 8947 need remove atom 8949 need remove atom 8952 need remove atom 8958 need remove atom 8960 need remove atom 8975 need remove atom 8983 need remove atom 8985 need remove atom 8988 need remove atom 8994 need remove atom 8996 need remove atom 9011 need remove atom 9019 need remove atom 9021 need remove atom 9024 need remove atom 9030 need remove atom 9032 need remove atom 9047 need remove atom 9055 need remove atom 9057 need remove atom 9060 need remove atom 9066 need remove atom 9068 need remove atom 9083 need remove atom 9091 need remove atom 9093 need remove atom 9096 need remove atom 9102 need remove atom 9104 need remove atom 9119 need remove atom 9127 need remove atom 9129 need remove atom 9132 need remove atom 9138 need remove atom 9140 need remove atom 9155 need remove atom 9163 need remove atom 9165 need remove atom 9168 need remove atom 9174 need remove atom 9176 need remove atom 9191 need remove atom 9199 need remove atom 9201 need remove atom 9204 need remove atom 9210 need remove atom 9212 need remove atom 9227 need remove atom 9235 need remove atom 9237 need remove atom 9240 need remove atom 9246 need remove atom 9248 need remove atom 9263 need remove atom 9271 need remove atom 9273 need remove atom 9276 need remove atom 9282 need remove atom 9284 need remove atom 9299 need remove atom 9307 need remove atom 9309 need remove atom 9312 need remove atom 9318 need remove atom 9320 need remove atom 9335 need remove atom 9343 need remove atom 9345 need remove atom 9348 need remove atom 9354 need remove atom 9356 need remove atom 9371 need remove atom 9379 need remove atom 9381 need remove atom 9384 need remove atom 9390 need remove atom 9392 need remove atom 9407 need remove atom 9415 need remove atom 9417 need remove atom 9420 need remove atom 9426 need remove atom 9428 need remove atom 9443 need remove atom 9451 need remove atom 9453 need remove atom 9456 need remove atom 9462 need remove atom 9464 need remove atom 9479 need remove atom 9487 need remove atom 9489 need remove atom 9492 need remove atom 9498 need remove atom 9500 need remove atom 9515 need remove atom 9523 need remove atom 9525 need remove atom 9528 need remove atom 9534 need remove atom 9536 need remove atom 9551 need remove atom 9559 need remove atom 9561 need remove atom 9564 need remove atom 9570 need remove atom 9572 need remove atom 9587 need remove atom 9595 need remove atom 9597 need remove atom 9600 need remove atom 9606 need remove atom 9608 need remove atom 9623 need remove atom 9631 need remove atom 9633 need remove atom 9636 need remove atom 9642 need remove atom 9644 need remove atom 9659 need remove atom 9667 need remove atom 9669 need remove atom 9672 need remove atom 9678 need remove atom 9680 need remove atom 9695 need remove atom 9703 need remove atom 9705 need remove atom 9708 need remove atom 9714 need remove atom 9716 need remove atom 9731 need remove atom 9739 need remove atom 9741 need remove atom 9744 need remove atom 9750 need remove atom 9752 need remove atom 9767 need remove atom 9775 need remove atom 9777 need remove atom 9780 need remove atom 9786 need remove atom 9788 need remove atom 9803 need remove atom 9811 need remove atom 9813 need remove atom 9816 need remove atom 9822 need remove atom 9824 need remove atom 9839 need remove atom 9847 need remove atom 9849 need remove atom 9852 need remove atom 9858 need remove atom 9860 need remove atom 9875 need remove atom 9883 need remove atom 9885 need remove atom 9888 need remove atom 9894 need remove atom 9896 need remove atom 9911 need remove atom 9919 need remove atom 9921 need remove atom 9924 need remove atom 9930 need remove atom 9932 [I] need removenum=174 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2196755306e-01 1.0010998274e-01 -3.4923900536e-01 dub= -3.3773250104e-01 1.0011234751e-01 3.5584498461e-01 duc= 3.2196845943e-01 -1.0142703831e-01 -3.4924012360e-01 [I] overall tilt ux= -6.5970005410e-01 2.3647702652e-06 7.0508398997e-01 uy= 9.0636999395e-07 -2.0153702105e-01 -1.1182354171e-06 [I] storedr[0]= -1.9716472414e-02 -6.6495814173e-04 -2.9405354811e-02 [I] originally 18792 atoms [I] insert 0 atoms [I] now 18792 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 18792 atoms [I] remove 174 atoms [I] now 18618 atoms [I] NP=18618 n=18618 [I] ASSIGN finalcnfile = dipole_2.5.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_2.5.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_2.5.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:22:38 2023 CPU time spent: 1.576538 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:23:56 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 10 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 32 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=34560 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.016666666666667 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.055555555555556 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 1.6666666667e-02 -0.0000000000e+00 -5.5555555556e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=9.494548e+01 pxpb=0.000000e+00 pypb=1.432247e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 14693 need remove atom 14695 need remove atom 14696 need remove atom 14729 need remove atom 14731 need remove atom 14732 need remove atom 14765 need remove atom 14767 need remove atom 14768 need remove atom 14796 need remove atom 14797 need remove atom 14801 need remove atom 14803 need remove atom 14804 need remove atom 14812 need remove atom 14832 need remove atom 14833 need remove atom 14837 need remove atom 14839 need remove atom 14840 need remove atom 14848 need remove atom 14868 need remove atom 14869 need remove atom 14873 need remove atom 14875 need remove atom 14876 need remove atom 14884 need remove atom 14904 need remove atom 14905 need remove atom 14909 need remove atom 14911 need remove atom 14912 need remove atom 14920 need remove atom 14940 need remove atom 14941 need remove atom 14945 need remove atom 14947 need remove atom 14948 need remove atom 14956 need remove atom 14976 need remove atom 14977 need remove atom 14981 need remove atom 14983 need remove atom 14984 need remove atom 14992 need remove atom 15012 need remove atom 15013 need remove atom 15017 need remove atom 15019 need remove atom 15020 need remove atom 15028 need remove atom 15048 need remove atom 15049 need remove atom 15053 need remove atom 15055 need remove atom 15056 need remove atom 15064 need remove atom 15084 need remove atom 15085 need remove atom 15089 need remove atom 15091 need remove atom 15092 need remove atom 15100 need remove atom 15120 need remove atom 15121 need remove atom 15125 need remove atom 15127 need remove atom 15128 need remove atom 15136 need remove atom 15156 need remove atom 15157 need remove atom 15161 need remove atom 15163 need remove atom 15164 need remove atom 15172 need remove atom 15192 need remove atom 15193 need remove atom 15197 need remove atom 15199 need remove atom 15200 need remove atom 15208 need remove atom 15228 need remove atom 15229 need remove atom 15233 need remove atom 15235 need remove atom 15236 need remove atom 15244 need remove atom 15264 need remove atom 15265 need remove atom 15269 need remove atom 15271 need remove atom 15272 need remove atom 15280 need remove atom 15300 need remove atom 15301 need remove atom 15305 need remove atom 15307 need remove atom 15308 need remove atom 15316 need remove atom 15336 need remove atom 15337 need remove atom 15341 need remove atom 15343 need remove atom 15344 need remove atom 15352 need remove atom 15372 need remove atom 15373 need remove atom 15377 need remove atom 15379 need remove atom 15380 need remove atom 15388 need remove atom 15408 need remove atom 15409 need remove atom 15413 need remove atom 15415 need remove atom 15416 need remove atom 15424 need remove atom 15444 need remove atom 15445 need remove atom 15449 need remove atom 15451 need remove atom 15452 need remove atom 15460 need remove atom 15480 need remove atom 15481 need remove atom 15485 need remove atom 15487 need remove atom 15488 need remove atom 15496 need remove atom 15516 need remove atom 15517 need remove atom 15521 need remove atom 15523 need remove atom 15524 need remove atom 15532 need remove atom 15552 need remove atom 15553 need remove atom 15557 need remove atom 15559 need remove atom 15560 need remove atom 15568 need remove atom 15588 need remove atom 15589 need remove atom 15593 need remove atom 15595 need remove atom 15596 need remove atom 15604 need remove atom 15624 need remove atom 15625 need remove atom 15629 need remove atom 15631 need remove atom 15632 need remove atom 15640 need remove atom 15660 need remove atom 15661 need remove atom 15665 need remove atom 15667 need remove atom 15668 need remove atom 15676 need remove atom 15696 need remove atom 15697 need remove atom 15701 need remove atom 15703 need remove atom 15704 need remove atom 15712 need remove atom 15732 need remove atom 15733 need remove atom 15737 need remove atom 15739 need remove atom 15740 need remove atom 15748 need remove atom 15768 need remove atom 15769 need remove atom 15773 need remove atom 15775 need remove atom 15776 need remove atom 15784 need remove atom 15804 need remove atom 15805 need remove atom 15809 need remove atom 15811 need remove atom 15812 need remove atom 15820 need remove atom 15840 need remove atom 15841 need remove atom 15845 need remove atom 15847 need remove atom 15848 need remove atom 15856 need remove atom 15876 need remove atom 15877 need remove atom 15881 need remove atom 15883 need remove atom 15884 need remove atom 15892 need remove atom 15912 need remove atom 15913 need remove atom 15917 need remove atom 15919 need remove atom 15920 need remove atom 15928 need remove atom 15948 need remove atom 15949 need remove atom 15953 need remove atom 15955 need remove atom 15956 need remove atom 15964 need remove atom 15984 need remove atom 15985 need remove atom 15989 need remove atom 15991 need remove atom 15992 need remove atom 16000 need remove atom 16020 need remove atom 16021 need remove atom 16025 need remove atom 16027 need remove atom 16028 need remove atom 16036 need remove atom 16056 need remove atom 16057 need remove atom 16061 need remove atom 16063 need remove atom 16064 need remove atom 16072 need remove atom 16092 need remove atom 16093 need remove atom 16097 need remove atom 16099 need remove atom 16100 need remove atom 16108 need remove atom 16128 need remove atom 16129 need remove atom 16133 need remove atom 16135 need remove atom 16136 need remove atom 16144 need remove atom 16164 need remove atom 16165 need remove atom 16169 need remove atom 16171 need remove atom 16172 need remove atom 16180 need remove atom 16200 need remove atom 16201 need remove atom 16205 need remove atom 16207 need remove atom 16208 need remove atom 16216 need remove atom 16236 need remove atom 16237 need remove atom 16241 need remove atom 16243 need remove atom 16244 need remove atom 16252 need remove atom 16272 need remove atom 16273 need remove atom 16277 need remove atom 16279 need remove atom 16280 need remove atom 16288 need remove atom 16308 need remove atom 16309 need remove atom 16313 need remove atom 16315 need remove atom 16316 need remove atom 16324 need remove atom 16344 need remove atom 16345 need remove atom 16349 need remove atom 16351 need remove atom 16352 need remove atom 16360 need remove atom 16380 need remove atom 16381 need remove atom 16385 need remove atom 16387 need remove atom 16388 need remove atom 16396 need remove atom 16416 need remove atom 16417 need remove atom 16432 need remove atom 16452 need remove atom 16453 need remove atom 16468 need remove atom 16488 need remove atom 16489 need remove atom 16504 [I] need removenum=288 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2141155831e-01 1.0030255369e-01 -3.4818189825e-01 dub= -3.3711888039e-01 1.0030485216e-01 3.5463137594e-01 duc= 3.2141246976e-01 -1.0160575309e-01 -3.4818301481e-01 [I] overall tilt ux= -6.5853043870e-01 2.2984762599e-06 7.0281327419e-01 uy= 9.1144464465e-07 -2.0190830677e-01 -1.1165637032e-06 [I] storedr[0]= -1.8395840399e-02 -6.5804840090e-04 -2.6956368915e-02 [I] originally 34560 atoms [I] insert 0 atoms [I] now 34560 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 34560 atoms [I] remove 288 atoms [I] now 34272 atoms [I] NP=34272 n=34272 [I] ASSIGN finalcnfile = dipole_3.0.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_3.0.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_3.0.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:24:00 2023 CPU time spent: 2.926217 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:28:01 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 12 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 38 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=49248 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.013888888888889 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.055555555555556 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 1.3888888889e-02 -0.0000000000e+00 -5.5555555556e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=1.139346e+02 pxpb=0.000000e+00 pypb=1.700794e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 21600 need remove atom 21601 need remove atom 21605 need remove atom 21607 need remove atom 21608 need remove atom 21616 need remove atom 21636 need remove atom 21637 need remove atom 21641 need remove atom 21643 need remove atom 21644 need remove atom 21652 need remove atom 21672 need remove atom 21673 need remove atom 21677 need remove atom 21679 need remove atom 21680 need remove atom 21688 need remove atom 21708 need remove atom 21709 need remove atom 21713 need remove atom 21715 need remove atom 21716 need remove atom 21724 need remove atom 21744 need remove atom 21745 need remove atom 21749 need remove atom 21751 need remove atom 21752 need remove atom 21760 need remove atom 21780 need remove atom 21781 need remove atom 21785 need remove atom 21787 need remove atom 21788 need remove atom 21796 need remove atom 21816 need remove atom 21817 need remove atom 21821 need remove atom 21823 need remove atom 21824 need remove atom 21832 need remove atom 21852 need remove atom 21853 need remove atom 21857 need remove atom 21859 need remove atom 21860 need remove atom 21868 need remove atom 21888 need remove atom 21889 need remove atom 21893 need remove atom 21895 need remove atom 21896 need remove atom 21904 need remove atom 21924 need remove atom 21925 need remove atom 21929 need remove atom 21931 need remove atom 21932 need remove atom 21940 need remove atom 21960 need remove atom 21961 need remove atom 21965 need remove atom 21967 need remove atom 21968 need remove atom 21976 need remove atom 21996 need remove atom 21997 need remove atom 22001 need remove atom 22003 need remove atom 22004 need remove atom 22012 need remove atom 22032 need remove atom 22033 need remove atom 22037 need remove atom 22039 need remove atom 22040 need remove atom 22048 need remove atom 22068 need remove atom 22069 need remove atom 22073 need remove atom 22075 need remove atom 22076 need remove atom 22084 need remove atom 22104 need remove atom 22105 need remove atom 22109 need remove atom 22111 need remove atom 22112 need remove atom 22120 need remove atom 22140 need remove atom 22141 need remove atom 22145 need remove atom 22147 need remove atom 22148 need remove atom 22156 need remove atom 22176 need remove atom 22177 need remove atom 22181 need remove atom 22183 need remove atom 22184 need remove atom 22192 need remove atom 22212 need remove atom 22213 need remove atom 22217 need remove atom 22219 need remove atom 22220 need remove atom 22228 need remove atom 22248 need remove atom 22249 need remove atom 22253 need remove atom 22255 need remove atom 22256 need remove atom 22264 need remove atom 22284 need remove atom 22285 need remove atom 22289 need remove atom 22291 need remove atom 22292 need remove atom 22300 need remove atom 22320 need remove atom 22321 need remove atom 22325 need remove atom 22327 need remove atom 22328 need remove atom 22336 need remove atom 22356 need remove atom 22357 need remove atom 22361 need remove atom 22363 need remove atom 22364 need remove atom 22372 need remove atom 22392 need remove atom 22393 need remove atom 22397 need remove atom 22399 need remove atom 22400 need remove atom 22408 need remove atom 22428 need remove atom 22429 need remove atom 22433 need remove atom 22435 need remove atom 22436 need remove atom 22444 need remove atom 22464 need remove atom 22465 need remove atom 22469 need remove atom 22471 need remove atom 22472 need remove atom 22480 need remove atom 22500 need remove atom 22501 need remove atom 22505 need remove atom 22507 need remove atom 22508 need remove atom 22516 need remove atom 22536 need remove atom 22537 need remove atom 22541 need remove atom 22543 need remove atom 22544 need remove atom 22552 need remove atom 22572 need remove atom 22573 need remove atom 22577 need remove atom 22579 need remove atom 22580 need remove atom 22588 need remove atom 22608 need remove atom 22609 need remove atom 22613 need remove atom 22615 need remove atom 22616 need remove atom 22624 need remove atom 22644 need remove atom 22645 need remove atom 22649 need remove atom 22651 need remove atom 22652 need remove atom 22660 need remove atom 22680 need remove atom 22681 need remove atom 22685 need remove atom 22687 need remove atom 22688 need remove atom 22696 need remove atom 22716 need remove atom 22717 need remove atom 22721 need remove atom 22723 need remove atom 22724 need remove atom 22732 need remove atom 22752 need remove atom 22753 need remove atom 22757 need remove atom 22759 need remove atom 22760 need remove atom 22768 need remove atom 22788 need remove atom 22789 need remove atom 22793 need remove atom 22795 need remove atom 22796 need remove atom 22804 need remove atom 22824 need remove atom 22825 need remove atom 22829 need remove atom 22831 need remove atom 22832 need remove atom 22840 need remove atom 22860 need remove atom 22861 need remove atom 22865 need remove atom 22867 need remove atom 22868 need remove atom 22876 need remove atom 22896 need remove atom 22897 need remove atom 22901 need remove atom 22903 need remove atom 22904 need remove atom 22912 need remove atom 22932 need remove atom 22933 need remove atom 22937 need remove atom 22939 need remove atom 22940 need remove atom 22948 need remove atom 22968 need remove atom 22969 need remove atom 22973 need remove atom 22975 need remove atom 22976 need remove atom 22984 need remove atom 23004 need remove atom 23005 need remove atom 23009 need remove atom 23011 need remove atom 23012 need remove atom 23020 need remove atom 23040 need remove atom 23041 need remove atom 23045 need remove atom 23047 need remove atom 23048 need remove atom 23056 need remove atom 23076 need remove atom 23077 need remove atom 23081 need remove atom 23083 need remove atom 23084 need remove atom 23092 need remove atom 23112 need remove atom 23113 need remove atom 23117 need remove atom 23119 need remove atom 23120 need remove atom 23128 need remove atom 23148 need remove atom 23149 need remove atom 23153 need remove atom 23155 need remove atom 23156 need remove atom 23164 need remove atom 23184 need remove atom 23185 need remove atom 23189 need remove atom 23191 need remove atom 23192 need remove atom 23200 need remove atom 23220 need remove atom 23221 need remove atom 23225 need remove atom 23227 need remove atom 23228 need remove atom 23236 need remove atom 23256 need remove atom 23257 need remove atom 23261 need remove atom 23263 need remove atom 23264 need remove atom 23272 need remove atom 23292 need remove atom 23293 need remove atom 23297 need remove atom 23299 need remove atom 23300 need remove atom 23308 need remove atom 23328 need remove atom 23329 need remove atom 23333 need remove atom 23335 need remove atom 23336 need remove atom 23344 need remove atom 23364 need remove atom 23365 need remove atom 23369 need remove atom 23371 need remove atom 23372 need remove atom 23380 need remove atom 23400 need remove atom 23401 need remove atom 23405 need remove atom 23407 need remove atom 23408 need remove atom 23416 need remove atom 23436 need remove atom 23437 need remove atom 23441 need remove atom 23443 need remove atom 23444 need remove atom 23452 need remove atom 23472 need remove atom 23473 need remove atom 23477 need remove atom 23479 need remove atom 23480 need remove atom 23488 need remove atom 23508 need remove atom 23509 need remove atom 23513 need remove atom 23515 need remove atom 23516 need remove atom 23524 need remove atom 23544 need remove atom 23545 need remove atom 23549 need remove atom 23551 need remove atom 23552 need remove atom 23560 need remove atom 23580 need remove atom 23581 need remove atom 23585 need remove atom 23587 need remove atom 23588 need remove atom 23596 need remove atom 23616 need remove atom 23617 need remove atom 23621 need remove atom 23623 need remove atom 23624 need remove atom 23632 [I] need removenum=342 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2056269483e-01 1.0059376853e-01 -3.4657730662e-01 dub= -3.3617953780e-01 1.0059596403e-01 3.5278858480e-01 duc= 3.2056361385e-01 -1.0187597663e-01 -3.4657842049e-01 [I] overall tilt ux= -6.5674223264e-01 2.1955042352e-06 6.9936589142e-01 uy= 9.1901394494e-07 -2.0246974515e-01 -1.1138667277e-06 [I] storedr[0]= -1.6431888044e-02 -6.4703129831e-04 -2.2951179875e-02 [I] originally 49248 atoms [I] insert 0 atoms [I] now 49248 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 49248 atoms [I] remove 342 atoms [I] now 48906 atoms [I] NP=48906 n=48906 [I] ASSIGN finalcnfile = dipole_3.5.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_3.5.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_3.5.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:28:06 2023 CPU time spent: 4.206639 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:33:30 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 14 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 45 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=68040 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.011904761904762 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.055555555555556 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 1.1904761905e-02 -0.0000000000e+00 -5.5555555556e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=1.329237e+02 pxpb=0.000000e+00 pypb=2.014098e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 30351 need remove atom 30353 need remove atom 30355 need remove atom 30387 need remove atom 30389 need remove atom 30391 need remove atom 30423 need remove atom 30425 need remove atom 30427 need remove atom 30456 need remove atom 30459 need remove atom 30461 need remove atom 30463 need remove atom 30470 need remove atom 30472 need remove atom 30492 need remove atom 30495 need remove atom 30497 need remove atom 30499 need remove atom 30506 need remove atom 30508 need remove atom 30528 need remove atom 30531 need remove atom 30533 need remove atom 30535 need remove atom 30542 need remove atom 30544 need remove atom 30564 need remove atom 30567 need remove atom 30569 need remove atom 30571 need remove atom 30578 need remove atom 30580 need remove atom 30600 need remove atom 30603 need remove atom 30605 need remove atom 30607 need remove atom 30614 need remove atom 30616 need remove atom 30636 need remove atom 30639 need remove atom 30641 need remove atom 30643 need remove atom 30650 need remove atom 30652 need remove atom 30672 need remove atom 30675 need remove atom 30677 need remove atom 30679 need remove atom 30686 need remove atom 30688 need remove atom 30708 need remove atom 30711 need remove atom 30713 need remove atom 30715 need remove atom 30722 need remove atom 30724 need remove atom 30744 need remove atom 30747 need remove atom 30749 need remove atom 30751 need remove atom 30758 need remove atom 30760 need remove atom 30780 need remove atom 30783 need remove atom 30785 need remove atom 30787 need remove atom 30794 need remove atom 30796 need remove atom 30816 need remove atom 30819 need remove atom 30821 need remove atom 30823 need remove atom 30830 need remove atom 30832 need remove atom 30852 need remove atom 30855 need remove atom 30857 need remove atom 30859 need remove atom 30866 need remove atom 30868 need remove atom 30888 need remove atom 30891 need remove atom 30893 need remove atom 30895 need remove atom 30902 need remove atom 30904 need remove atom 30924 need remove atom 30927 need remove atom 30929 need remove atom 30931 need remove atom 30938 need remove atom 30940 need remove atom 30960 need remove atom 30963 need remove atom 30965 need remove atom 30967 need remove atom 30974 need remove atom 30976 need remove atom 30996 need remove atom 30999 need remove atom 31001 need remove atom 31003 need remove atom 31010 need remove atom 31012 need remove atom 31032 need remove atom 31035 need remove atom 31037 need remove atom 31039 need remove atom 31046 need remove atom 31048 need remove atom 31068 need remove atom 31071 need remove atom 31073 need remove atom 31075 need remove atom 31082 need remove atom 31084 need remove atom 31104 need remove atom 31107 need remove atom 31109 need remove atom 31111 need remove atom 31118 need remove atom 31120 need remove atom 31140 need remove atom 31143 need remove atom 31145 need remove atom 31147 need remove atom 31154 need remove atom 31156 need remove atom 31176 need remove atom 31179 need remove atom 31181 need remove atom 31183 need remove atom 31190 need remove atom 31192 need remove atom 31212 need remove atom 31215 need remove atom 31217 need remove atom 31219 need remove atom 31226 need remove atom 31228 need remove atom 31248 need remove atom 31251 need remove atom 31253 need remove atom 31255 need remove atom 31262 need remove atom 31264 need remove atom 31284 need remove atom 31287 need remove atom 31289 need remove atom 31291 need remove atom 31298 need remove atom 31300 need remove atom 31320 need remove atom 31323 need remove atom 31325 need remove atom 31327 need remove atom 31334 need remove atom 31336 need remove atom 31356 need remove atom 31359 need remove atom 31361 need remove atom 31363 need remove atom 31370 need remove atom 31372 need remove atom 31392 need remove atom 31395 need remove atom 31397 need remove atom 31399 need remove atom 31406 need remove atom 31408 need remove atom 31428 need remove atom 31431 need remove atom 31433 need remove atom 31435 need remove atom 31442 need remove atom 31444 need remove atom 31464 need remove atom 31467 need remove atom 31469 need remove atom 31471 need remove atom 31478 need remove atom 31480 need remove atom 31500 need remove atom 31503 need remove atom 31505 need remove atom 31507 need remove atom 31514 need remove atom 31516 need remove atom 31536 need remove atom 31539 need remove atom 31541 need remove atom 31543 need remove atom 31550 need remove atom 31552 need remove atom 31572 need remove atom 31575 need remove atom 31577 need remove atom 31579 need remove atom 31586 need remove atom 31588 need remove atom 31608 need remove atom 31611 need remove atom 31613 need remove atom 31615 need remove atom 31622 need remove atom 31624 need remove atom 31644 need remove atom 31647 need remove atom 31649 need remove atom 31651 need remove atom 31658 need remove atom 31660 need remove atom 31680 need remove atom 31683 need remove atom 31685 need remove atom 31687 need remove atom 31694 need remove atom 31696 need remove atom 31716 need remove atom 31719 need remove atom 31721 need remove atom 31723 need remove atom 31730 need remove atom 31732 need remove atom 31752 need remove atom 31755 need remove atom 31757 need remove atom 31759 need remove atom 31766 need remove atom 31768 need remove atom 31788 need remove atom 31791 need remove atom 31793 need remove atom 31795 need remove atom 31802 need remove atom 31804 need remove atom 31824 need remove atom 31827 need remove atom 31829 need remove atom 31831 need remove atom 31838 need remove atom 31840 need remove atom 31860 need remove atom 31863 need remove atom 31865 need remove atom 31867 need remove atom 31874 need remove atom 31876 need remove atom 31896 need remove atom 31899 need remove atom 31901 need remove atom 31903 need remove atom 31910 need remove atom 31912 need remove atom 31932 need remove atom 31935 need remove atom 31937 need remove atom 31939 need remove atom 31946 need remove atom 31948 need remove atom 31968 need remove atom 31971 need remove atom 31973 need remove atom 31975 need remove atom 31982 need remove atom 31984 need remove atom 32004 need remove atom 32007 need remove atom 32009 need remove atom 32011 need remove atom 32018 need remove atom 32020 need remove atom 32040 need remove atom 32043 need remove atom 32045 need remove atom 32047 need remove atom 32054 need remove atom 32056 need remove atom 32076 need remove atom 32079 need remove atom 32081 need remove atom 32083 need remove atom 32090 need remove atom 32092 need remove atom 32112 need remove atom 32115 need remove atom 32117 need remove atom 32119 need remove atom 32126 need remove atom 32128 need remove atom 32148 need remove atom 32151 need remove atom 32153 need remove atom 32155 need remove atom 32162 need remove atom 32164 need remove atom 32184 need remove atom 32187 need remove atom 32189 need remove atom 32191 need remove atom 32198 need remove atom 32200 need remove atom 32220 need remove atom 32223 need remove atom 32225 need remove atom 32227 need remove atom 32234 need remove atom 32236 need remove atom 32256 need remove atom 32259 need remove atom 32261 need remove atom 32263 need remove atom 32270 need remove atom 32272 need remove atom 32292 need remove atom 32295 need remove atom 32297 need remove atom 32299 need remove atom 32306 need remove atom 32308 need remove atom 32328 need remove atom 32331 need remove atom 32333 need remove atom 32335 need remove atom 32342 need remove atom 32344 need remove atom 32364 need remove atom 32367 need remove atom 32369 need remove atom 32371 need remove atom 32378 need remove atom 32380 need remove atom 32400 need remove atom 32403 need remove atom 32405 need remove atom 32407 need remove atom 32414 need remove atom 32416 need remove atom 32436 need remove atom 32439 need remove atom 32441 need remove atom 32443 need remove atom 32450 need remove atom 32452 need remove atom 32472 need remove atom 32475 need remove atom 32477 need remove atom 32479 need remove atom 32486 need remove atom 32488 need remove atom 32508 need remove atom 32511 need remove atom 32513 need remove atom 32515 need remove atom 32522 need remove atom 32524 need remove atom 32544 need remove atom 32547 need remove atom 32549 need remove atom 32551 need remove atom 32558 need remove atom 32560 need remove atom 32580 need remove atom 32583 need remove atom 32585 need remove atom 32587 need remove atom 32594 need remove atom 32596 need remove atom 32616 need remove atom 32619 need remove atom 32621 need remove atom 32623 need remove atom 32630 need remove atom 32632 need remove atom 32652 need remove atom 32655 need remove atom 32657 need remove atom 32659 need remove atom 32666 need remove atom 32668 need remove atom 32688 need remove atom 32691 need remove atom 32693 need remove atom 32695 need remove atom 32702 need remove atom 32704 need remove atom 32724 need remove atom 32727 need remove atom 32729 need remove atom 32731 need remove atom 32738 need remove atom 32740 need remove atom 32760 need remove atom 32763 need remove atom 32765 need remove atom 32767 need remove atom 32774 need remove atom 32776 need remove atom 32796 need remove atom 32799 need remove atom 32801 need remove atom 32803 need remove atom 32810 need remove atom 32812 [I] need removenum=405 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2176988181e-01 1.0017861216e-01 -3.4886261558e-01 dub= -3.3751449228e-01 1.0018095346e-01 3.5541291136e-01 duc= 3.2177079000e-01 -1.0149073172e-01 -3.4886373323e-01 [I] overall tilt ux= -6.5928437410e-01 2.3413075771e-06 7.0427552694e-01 uy= 9.0818475146e-07 -2.0166934388e-01 -1.1176498014e-06 [I] storedr[0]= -1.5476138652e-02 -6.6087783565e-04 -1.8535814344e-02 [I] originally 68040 atoms [I] insert 0 atoms [I] now 68040 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 68040 atoms [I] remove 405 atoms [I] now 67635 atoms [I] NP=67635 n=67635 [I] ASSIGN finalcnfile = dipole_4.0.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_4.0.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_4.0.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:33:37 2023 CPU time spent: 5.789666 s [I] MDFrame::initvars() [I] MDFrame::initparser() [I] EXEC setnolog [I] EXEC setoverwrite [I] ASSIGN dirname = output/lammps/in/dipole_ref [W] directory output/lammps/in/dipole_ref already exists [I] RERUN dirname = output/lammps/in/dipole_ref [I] Begin Time = Wed May 31 21:37:07 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1648494526744 [I] ASSIGN makecnspec = -2 [I] ASSIGN makecnspec(1) = -2 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 15 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 48 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 4 [I] ASSIGN makecnspec(11) = 4 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ -2 -1 1 -2 1 1 1 0 4]; basis=[ 0.833333 0 0.166667 0 0 0.944444 0 0.388889 0 1 0.5 0.5 0 0 2 0.777778 0.5 0.0555556 0 3 0.611111 0.5 0.222222 0 4 0.888889 0.5 0.277778 0 5 0.722222 0.5 0.444444 0 6 0.833333 0.5 0.666667 0 7 0.944444 0.5 0.888889 0 8 0.277778 0 0.0555556 0 9 0.555556 0 0.111111 0 10 0.388889 0 0.277778 0 11 0.666667 0 0.333333 0 12 0.5 0 0.5 0 13 0.777778 0 0.555556 0 14 0.611111 0 0.722222 0 15 0.888889 0 0.777778 0 16 0.722222 0 0.944444 0 17 0.0555556 0.5 0.111111 0 18 0.333333 0.5 0.166667 0 19 0.166667 0.5 0.333333 0 20 0.444444 0.5 0.388889 0 21 0.277778 0.5 0.555556 0 22 0.555556 0.5 0.611111 0 23 0.388889 0.5 0.777778 0 24 0.666667 0.5 0.833333 0 25 0 0 0 0 26 0.111111 0 0.222222 0 27 0.222222 0 0.444444 0 28 0.0555556 0 0.611111 0 29 0.333333 0 0.666667 0 30 0.166667 0 0.833333 0 31 0.444444 0 0.888889 0 32 0 0.5 0.5 0 33 0.111111 0.5 0.722222 0 34 0.222222 0.5 0.944444 0 35 ]; [I] makecn: _NP=103680 [I] ASSIGN finalcnfile = perf.cn [I] EXEC writecn [I] FILEOPEN perf.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE perf.cn [I] ASSIGN mkdipole = 3 [I] ASSIGN mkdipole(1) = 2 [I] ASSIGN mkdipole(2) = 0.011111111111111 [I] ASSIGN mkdipole(3) = -0 [I] ASSIGN mkdipole(4) = -0.041666666666667 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28101794975058 [I] ASSIGN mkdipole(9) = -10 [I] ASSIGN mkdipole(10) = 10 [I] ASSIGN mkdipole(11) = -10 [I] ASSIGN mkdipole(12) = 10 [I] ASSIGN mkdipole(13) = 1 [I] EXEC makedipole [I] sb= 1.1111111111e-02 -0.0000000000e+00 -4.1666666667e-02 b= 1.5824247263e+00 0.0000000000e+00 -2.2378865094e+00 [I] px= 1.0000000000e+00 0.0000000000e+00 0.0000000000e+00 [I] py= 0.0000000000e+00 1.0000000000e+00 0.0000000000e+00 [I] pz= 0.0000000000e+00 0.0000000000e+00 1.0000000000e+00 [I] bx=1.5824247263372 by=0 pxpa=1.424182e+02 pxpb=0.000000e+00 pypb=2.148371e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 50131 need remove atom 50134 need remove atom 50136 need remove atom 50167 need remove atom 50170 need remove atom 50172 need remove atom 50203 need remove atom 50206 need remove atom 50208 need remove atom 50239 need remove atom 50242 need remove atom 50244 need remove atom 50265 need remove atom 50267 need remove atom 50275 need remove atom 50278 need remove atom 50280 need remove atom 50286 need remove atom 50301 need remove atom 50303 need remove atom 50311 need remove atom 50314 need remove atom 50316 need remove atom 50322 need remove atom 50337 need remove atom 50339 need remove atom 50347 need remove atom 50350 need remove atom 50352 need remove atom 50358 need remove atom 50373 need remove atom 50375 need remove atom 50383 need remove atom 50386 need remove atom 50388 need remove atom 50394 need remove atom 50409 need remove atom 50411 need remove atom 50419 need remove atom 50422 need remove atom 50424 need remove atom 50430 need remove atom 50445 need remove atom 50447 need remove atom 50455 need remove atom 50458 need remove atom 50460 need remove atom 50466 need remove atom 50481 need remove atom 50483 need remove atom 50491 need remove atom 50494 need remove atom 50496 need remove atom 50502 need remove atom 50517 need remove atom 50519 need remove atom 50527 need remove atom 50530 need remove atom 50532 need remove atom 50538 need remove atom 50553 need remove atom 50555 need remove atom 50563 need remove atom 50566 need remove atom 50568 need remove atom 50574 need remove atom 50589 need remove atom 50591 need remove atom 50599 need remove atom 50602 need remove atom 50604 need remove atom 50610 need remove atom 50625 need remove atom 50627 need remove atom 50635 need remove atom 50638 need remove atom 50640 need remove atom 50646 need remove atom 50661 need remove atom 50663 need remove atom 50671 need remove atom 50674 need remove atom 50676 need remove atom 50682 need remove atom 50697 need remove atom 50699 need remove atom 50707 need remove atom 50710 need remove atom 50712 need remove atom 50718 need remove atom 50733 need remove atom 50735 need remove atom 50743 need remove atom 50746 need remove atom 50748 need remove atom 50754 need remove atom 50769 need remove atom 50771 need remove atom 50779 need remove atom 50782 need remove atom 50784 need remove atom 50790 need remove atom 50805 need remove atom 50807 need remove atom 50815 need remove atom 50818 need remove atom 50820 need remove atom 50826 need remove atom 50841 need remove atom 50843 need remove atom 50851 need remove atom 50854 need remove atom 50856 need remove atom 50862 need remove atom 50877 need remove atom 50879 need remove atom 50887 need remove atom 50890 need remove atom 50892 need remove atom 50898 need remove atom 50913 need remove atom 50915 need remove atom 50923 need remove atom 50926 need remove atom 50928 need remove atom 50934 need remove atom 50949 need remove atom 50951 need remove atom 50959 need remove atom 50962 need remove atom 50964 need remove atom 50970 need remove atom 50985 need remove atom 50987 need remove atom 50995 need remove atom 50998 need remove atom 51000 need remove atom 51006 need remove atom 51021 need remove atom 51023 need remove atom 51031 need remove atom 51034 need remove atom 51036 need remove atom 51042 need remove atom 51057 need remove atom 51059 need remove atom 51067 need remove atom 51070 need remove atom 51072 need remove atom 51078 need remove atom 51093 need remove atom 51095 need remove atom 51103 need remove atom 51106 need remove atom 51108 need remove atom 51114 need remove atom 51129 need remove atom 51131 need remove atom 51139 need remove atom 51142 need remove atom 51144 need remove atom 51150 need remove atom 51165 need remove atom 51167 need remove atom 51175 need remove atom 51178 need remove atom 51180 need remove atom 51186 need remove atom 51201 need remove atom 51203 need remove atom 51211 need remove atom 51214 need remove atom 51216 need remove atom 51222 need remove atom 51237 need remove atom 51239 need remove atom 51247 need remove atom 51250 need remove atom 51252 need remove atom 51258 need remove atom 51273 need remove atom 51275 need remove atom 51283 need remove atom 51286 need remove atom 51288 need remove atom 51294 need remove atom 51309 need remove atom 51311 need remove atom 51319 need remove atom 51322 need remove atom 51324 need remove atom 51330 need remove atom 51345 need remove atom 51347 need remove atom 51355 need remove atom 51358 need remove atom 51360 need remove atom 51366 need remove atom 51381 need remove atom 51383 need remove atom 51391 need remove atom 51394 need remove atom 51396 need remove atom 51402 need remove atom 51417 need remove atom 51419 need remove atom 51427 need remove atom 51430 need remove atom 51432 need remove atom 51438 need remove atom 51453 need remove atom 51455 need remove atom 51463 need remove atom 51466 need remove atom 51468 need remove atom 51474 need remove atom 51489 need remove atom 51491 need remove atom 51499 need remove atom 51502 need remove atom 51504 need remove atom 51510 need remove atom 51525 need remove atom 51527 need remove atom 51535 need remove atom 51538 need remove atom 51540 need remove atom 51546 need remove atom 51561 need remove atom 51563 need remove atom 51571 need remove atom 51574 need remove atom 51576 need remove atom 51582 need remove atom 51597 need remove atom 51599 need remove atom 51607 need remove atom 51610 need remove atom 51612 need remove atom 51618 need remove atom 51633 need remove atom 51635 need remove atom 51643 need remove atom 51646 need remove atom 51648 need remove atom 51654 need remove atom 51669 need remove atom 51671 need remove atom 51679 need remove atom 51682 need remove atom 51684 need remove atom 51690 need remove atom 51705 need remove atom 51707 need remove atom 51715 need remove atom 51718 need remove atom 51720 need remove atom 51726 need remove atom 51741 need remove atom 51743 need remove atom 51751 need remove atom 51754 need remove atom 51756 need remove atom 51762 need remove atom 51777 need remove atom 51779 need remove atom 51787 need remove atom 51790 need remove atom 51792 need remove atom 51798 need remove atom 51813 need remove atom 51815 need remove atom 51823 need remove atom 51826 need remove atom 51828 need remove atom 51834 need remove atom 51849 need remove atom 51851 need remove atom 51859 need remove atom 51862 need remove atom 51864 need remove atom 51870 need remove atom 51885 need remove atom 51887 need remove atom 51895 need remove atom 51898 need remove atom 51900 need remove atom 51906 need remove atom 51921 need remove atom 51923 need remove atom 51931 need remove atom 51934 need remove atom 51936 need remove atom 51942 need remove atom 51957 need remove atom 51959 need remove atom 51967 need remove atom 51970 need remove atom 51972 need remove atom 51978 need remove atom 51993 need remove atom 51995 need remove atom 52003 need remove atom 52006 need remove atom 52008 need remove atom 52014 need remove atom 52029 need remove atom 52031 need remove atom 52039 need remove atom 52042 need remove atom 52044 need remove atom 52050 need remove atom 52065 need remove atom 52067 need remove atom 52075 need remove atom 52078 need remove atom 52080 need remove atom 52086 need remove atom 52101 need remove atom 52103 need remove atom 52111 need remove atom 52114 need remove atom 52116 need remove atom 52122 need remove atom 52137 need remove atom 52139 need remove atom 52147 need remove atom 52150 need remove atom 52152 need remove atom 52158 need remove atom 52173 need remove atom 52175 need remove atom 52183 need remove atom 52186 need remove atom 52188 need remove atom 52194 need remove atom 52209 need remove atom 52211 need remove atom 52219 need remove atom 52222 need remove atom 52224 need remove atom 52230 need remove atom 52245 need remove atom 52247 need remove atom 52255 need remove atom 52258 need remove atom 52260 need remove atom 52266 need remove atom 52281 need remove atom 52283 need remove atom 52291 need remove atom 52294 need remove atom 52296 need remove atom 52302 need remove atom 52317 need remove atom 52319 need remove atom 52327 need remove atom 52330 need remove atom 52332 need remove atom 52338 need remove atom 52353 need remove atom 52355 need remove atom 52363 need remove atom 52366 need remove atom 52368 need remove atom 52374 need remove atom 52389 need remove atom 52391 need remove atom 52399 need remove atom 52402 need remove atom 52404 need remove atom 52410 need remove atom 52425 need remove atom 52427 need remove atom 52435 need remove atom 52438 need remove atom 52440 need remove atom 52446 need remove atom 52461 need remove atom 52463 need remove atom 52471 need remove atom 52474 need remove atom 52476 need remove atom 52482 need remove atom 52497 need remove atom 52499 need remove atom 52507 need remove atom 52510 need remove atom 52512 need remove atom 52518 need remove atom 52533 need remove atom 52535 need remove atom 52543 need remove atom 52546 need remove atom 52548 need remove atom 52554 need remove atom 52569 need remove atom 52571 need remove atom 52579 need remove atom 52582 need remove atom 52584 need remove atom 52590 need remove atom 52605 need remove atom 52607 need remove atom 52615 need remove atom 52618 need remove atom 52620 need remove atom 52626 need remove atom 52641 need remove atom 52643 need remove atom 52651 need remove atom 52654 need remove atom 52656 need remove atom 52662 need remove atom 52677 need remove atom 52679 need remove atom 52687 need remove atom 52690 need remove atom 52692 need remove atom 52698 need remove atom 52713 need remove atom 52715 need remove atom 52723 need remove atom 52726 need remove atom 52728 need remove atom 52734 need remove atom 52749 need remove atom 52751 need remove atom 52759 need remove atom 52762 need remove atom 52764 need remove atom 52770 need remove atom 52785 need remove atom 52787 need remove atom 52795 need remove atom 52798 need remove atom 52800 need remove atom 52806 need remove atom 52821 need remove atom 52823 need remove atom 52831 need remove atom 52834 need remove atom 52836 need remove atom 52842 need remove atom 52857 need remove atom 52859 need remove atom 52867 need remove atom 52870 need remove atom 52872 need remove atom 52878 need remove atom 52893 need remove atom 52895 need remove atom 52903 need remove atom 52906 need remove atom 52908 need remove atom 52914 need remove atom 52929 need remove atom 52931 need remove atom 52939 need remove atom 52942 need remove atom 52944 need remove atom 52950 need remove atom 52965 need remove atom 52967 need remove atom 52975 need remove atom 52978 need remove atom 52980 need remove atom 52986 need remove atom 53001 need remove atom 53003 need remove atom 53011 need remove atom 53014 need remove atom 53016 need remove atom 53022 need remove atom 53037 need remove atom 53039 need remove atom 53047 need remove atom 53050 need remove atom 53052 need remove atom 53058 need remove atom 53073 need remove atom 53075 need remove atom 53083 need remove atom 53086 need remove atom 53088 need remove atom 53094 need remove atom 53109 need remove atom 53111 need remove atom 53119 need remove atom 53122 need remove atom 53124 need remove atom 53130 need remove atom 53145 need remove atom 53147 need remove atom 53155 need remove atom 53158 need remove atom 53160 need remove atom 53166 need remove atom 53181 need remove atom 53183 need remove atom 53191 need remove atom 53194 need remove atom 53196 need remove atom 53202 need remove atom 53217 need remove atom 53219 need remove atom 53227 need remove atom 53230 need remove atom 53232 need remove atom 53238 need remove atom 53253 need remove atom 53255 need remove atom 53263 need remove atom 53266 need remove atom 53268 need remove atom 53274 need remove atom 53289 need remove atom 53291 need remove atom 53299 need remove atom 53302 need remove atom 53304 need remove atom 53310 need remove atom 53325 need remove atom 53327 need remove atom 53335 need remove atom 53338 need remove atom 53340 need remove atom 53346 need remove atom 53361 need remove atom 53363 need remove atom 53371 need remove atom 53374 need remove atom 53376 need remove atom 53382 need remove atom 53397 need remove atom 53399 need remove atom 53407 need remove atom 53410 need remove atom 53412 need remove atom 53418 need remove atom 53433 need remove atom 53435 need remove atom 53443 need remove atom 53446 need remove atom 53448 need remove atom 53454 need remove atom 53469 need remove atom 53471 need remove atom 53479 need remove atom 53482 need remove atom 53484 need remove atom 53490 need remove atom 53505 need remove atom 53507 need remove atom 53515 need remove atom 53518 need remove atom 53520 need remove atom 53526 need remove atom 53541 need remove atom 53543 need remove atom 53551 need remove atom 53554 need remove atom 53556 need remove atom 53562 need remove atom 53577 need remove atom 53579 need remove atom 53598 need remove atom 53613 need remove atom 53615 need remove atom 53634 need remove atom 53649 need remove atom 53651 need remove atom 53670 need remove atom 53685 need remove atom 53687 need remove atom 53706 [I] need removenum=576 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2141155831e-01 1.0030255369e-01 -3.4818189825e-01 dub= -3.3711888039e-01 1.0030485216e-01 3.5463137594e-01 duc= 3.2141246976e-01 -1.0160575309e-01 -3.4818301481e-01 [I] overall tilt ux= -6.5853043870e-01 2.2984762599e-06 7.0281327419e-01 uy= 9.1144464454e-07 -2.0190830677e-01 -1.1165637028e-06 [I] storedr[0]= -1.4894763427e-02 -6.5617513045e-04 -1.7387560220e-02 [I] originally 103680 atoms [I] insert 0 atoms [I] now 103680 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 103680 atoms [I] remove 576 atoms [I] now 103104 atoms [I] NP=103104 n=103104 [I] ASSIGN finalcnfile = dipole_4.5.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_4.5.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_4.5.lammps [W] *********************************************** [W] * H is reoriented to upper triangular matrix! * [W] *********************************************** [I] writevelocity=0 [I] ASSIGN color00 = red [I] ASSIGN color01 = blue [I] ASSIGN color02 = green [I] ASSIGN color03 = magenta [I] ASSIGN color04 = cyan [I] ASSIGN color05 = purple [I] ASSIGN color06 = gray80 [I] ASSIGN color07 = white [I] ASSIGN atomradius = 0.67 [I] ASSIGN bondradius = 0.3 [I] ASSIGN bondlength = 0 [I] ASSIGN atomcolor = cyan [I] ASSIGN highlightcolor = purple [I] ASSIGN bondcolor = red [I] ASSIGN backgroundcolor = gray [I] ASSIGN fixatomcolor = yellow [I] ASSIGN plotfreq = 10 [I] ASSIGN rotateangles = 0 [I] ASSIGN rotateangles(1) = 0 [I] ASSIGN rotateangles(2) = 0 [I] ASSIGN rotateangles(3) = 1.2 [I] ASSIGN win_width = 500 [I] ASSIGN win_height = 500 [I] EXEC openwin [I] EXEC reversergb [I] EXEC alloccolors [W] No window to allocate color for! [I] EXEC rotate [W] No window to rotate for! [I] EXEC saverot [W] No window to rotate for! [I] ASSIGN NIC = 100 [I] ASSIGN NNM = 100 [I] ASSIGN plot_color_axis = 2 [I] EXEC refreshnnlist [I] ASSIGN NCS = 8 [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] ASSIGN finalcnfile = topol.cn [I] EXEC writecn [I] FILEOPEN topol.cn Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE topol.cn [I] ASSIGN plot_color_windows = 2 [I] ASSIGN plot_color_windows(1) = 0.2 [I] ASSIGN plot_color_windows(2) = 4 [I] ASSIGN plot_color_windows(3) = 0 [I] ASSIGN plot_color_windows(4) = 4 [I] ASSIGN plot_color_windows(5) = 20 [I] ASSIGN plot_color_windows(6) = 7 [I] EXEC plot [I] ASSIGN conj_ftol = 1e-07 [I] ASSIGN conj_fevalmax = 3000 [I] ASSIGN conj_fixbox = 1 [I] ASSIGN conj_dfpred = 0.0001 [I] ASSIGN stressmul = 0 [I] EXEC relax [I] Relax [I] Empty potential [I] relax: 1st potential call finished. [I] ############################################################ [I] iteration neval energy (eV) |gradient|^2 (eV^2) [I] ############################################################ 0 1 0.00000000000000e+00 0.00000000000000e+00 conj_f = 0 [I] ASSIGN finalcnfile = relaxed_b01_70_b02_50.lammps [I] EXEC calcentralsymmetry [I] ASSIGN writeall = 1 [I] EXEC writecn [I] FILEOPEN relaxed_b01_70_b02_50.lammps Configuration File for Atom Positions, Format: sx1 sy1 sz1; ... ; sxn syn szn; H(3x3) [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] EXEC writeLAMMPS [I] FILEOPEN relaxed_b01_70_b02_50.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE relaxed_b01_70_b02_50.lammps [I] writevelocity=0 [I] End Time = Wed May 31 21:37:17 2023 CPU time spent: 8.821373 s {"realtime":1513.80,"usertime":2963.66,"systime":28.19,"memmax":142300,"memavg":0}