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 24 14:27:58 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [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) = 11 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=1320 [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.05 [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.28133100175789 [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= 5.0000000000e-02 0.0000000000e+00 1.6666666667e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=3.165200e+01 pxpb=0.000000e+00 pypb=4.923896e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 697 need remove atom 699 need remove atom 701 need remove atom 703 need remove atom 705 need remove atom 707 need remove atom 709 need remove atom 711 need remove atom 713 need remove atom 715 need remove atom 717 need remove atom 719 need remove atom 721 need remove atom 723 need remove atom 725 need remove atom 727 need remove atom 729 need remove atom 731 need remove atom 733 need remove atom 735 need remove atom 737 need remove atom 739 need remove atom 741 need remove atom 743 need remove atom 745 need remove atom 747 need remove atom 749 need remove atom 751 need remove atom 753 need remove atom 755 need remove atom 757 need remove atom 761 need remove atom 765 [I] need removenum=33 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2392161994e-01 9.9402924774e-02 3.5290045102e-01 dub= -3.3987859857e-01 9.9405512757e-02 -3.6004071582e-01 duc= 3.2392250855e-01 -1.0076817378e-01 3.5290157447e-01 [I] overall tilt ux= -6.6380021851e-01 2.5879836478e-06 -7.1294116683e-01 uy= 8.8861152953e-07 -2.0017109856e-01 1.1234478368e-06 [I] storedr[0]= -1.3168128498e-02 6.4204306515e-02 1.5473952645e-02 [I] originally 1320 atoms [I] insert 0 atoms [I] now 1320 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 1320 atoms [I] remove 33 atoms [I] now 1287 atoms [I] NP=1287 n=1287 [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 24 14:27:58 2023 CPU time spent: 0.115414 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 24 14:28:17 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 11 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 12 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=1584 [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.045454545454545 [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.28133100175789 [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.5454545455e-02 0.0000000000e+00 1.6666666667e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=3.481720e+01 pxpb=0.000000e+00 pypb=5.371522e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 756 need remove atom 760 need remove atom 764 need remove atom 768 need remove atom 770 need remove atom 772 need remove atom 774 need remove atom 776 need remove atom 778 need remove atom 780 need remove atom 782 need remove atom 784 need remove atom 786 need remove atom 788 need remove atom 790 need remove atom 792 need remove atom 794 need remove atom 796 need remove atom 798 need remove atom 800 need remove atom 802 need remove atom 804 need remove atom 806 need remove atom 808 need remove atom 810 need remove atom 812 need remove atom 814 need remove atom 816 need remove atom 818 need remove atom 820 need remove atom 822 need remove atom 824 need remove atom 826 need remove atom 830 need remove atom 834 need remove atom 838 [I] need removenum=36 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2327083290e-01 9.9634158119e-02 3.5164366830e-01 dub= -3.3916557339e-01 9.9636672261e-02 -3.5859934701e-01 duc= 3.2327172781e-01 -1.0098285229e-01 3.5164479010e-01 [I] overall tilt ux= -6.6243640629e-01 2.5141422774e-06 -7.1024301531e-01 uy= 8.9491163269e-07 -2.0061701040e-01 1.1217955762e-06 [I] storedr[0]= -1.2773008182e-02 5.8299974165e-02 1.3958256689e-02 [I] originally 1584 atoms [I] insert 0 atoms [I] now 1584 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 1584 atoms [I] remove 36 atoms [I] now 1548 atoms [I] NP=1548 n=1548 [I] ASSIGN finalcnfile = dipole_1.1.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.1.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.1.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 24 14:28:17 2023 CPU time spent: 0.136875 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 24 14:28:41 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [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) = 13 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [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.28133100175789 [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.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=3.798240e+01 pxpb=0.000000e+00 pypb=5.819149e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 975 need remove atom 979 need remove atom 983 need remove atom 985 need remove atom 987 need remove atom 989 need remove atom 991 need remove atom 993 need remove atom 995 need remove atom 997 need remove atom 999 need remove atom 1001 need remove atom 1003 need remove atom 1005 need remove atom 1007 need remove atom 1009 need remove atom 1011 need remove atom 1013 need remove atom 1015 need remove atom 1017 need remove atom 1019 need remove atom 1021 need remove atom 1023 need remove atom 1025 need remove atom 1027 need remove atom 1029 need remove atom 1031 need remove atom 1033 need remove atom 1035 need remove atom 1037 need remove atom 1039 need remove atom 1041 need remove atom 1043 need remove atom 1045 need remove atom 1047 need remove atom 1049 need remove atom 1051 need remove atom 1053 need remove atom 1055 [I] need removenum=39 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2271996287e-01 9.9828368577e-02 3.5058523978e-01 dub= -3.3856055964e-01 9.9830819138e-02 -3.5738503760e-01 duc= 3.2272086302e-01 -1.0116314237e-01 3.5058636010e-01 [I] overall tilt ux= -6.6128052252e-01 2.4505612572e-06 -7.0797027738e-01 uy= 9.0014622822e-07 -2.0099151095e-01 1.1203107246e-06 [I] storedr[0]= -1.2416001239e-02 5.3370686223e-02 1.2689397321e-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.2.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.2.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.2.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 24 14:28:41 2023 CPU time spent: 0.160282 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 24 14:29:19 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 13 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 14 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=2184 [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.038461538461538 [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.28133100175789 [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.8461538462e-02 0.0000000000e+00 1.6666666667e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=4.114760e+01 pxpb=0.000000e+00 pypb=6.266776e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 1056 need remove atom 1058 need remove atom 1060 need remove atom 1062 need remove atom 1064 need remove atom 1066 need remove atom 1068 need remove atom 1070 need remove atom 1072 need remove atom 1074 need remove atom 1076 need remove atom 1078 need remove atom 1080 need remove atom 1082 need remove atom 1084 need remove atom 1086 need remove atom 1088 need remove atom 1090 need remove atom 1092 need remove atom 1094 need remove atom 1096 need remove atom 1098 need remove atom 1100 need remove atom 1102 need remove atom 1104 need remove atom 1106 need remove atom 1108 need remove atom 1110 need remove atom 1112 need remove atom 1114 need remove atom 1116 need remove atom 1118 need remove atom 1120 need remove atom 1122 need remove atom 1124 need remove atom 1126 need remove atom 1128 need remove atom 1130 need remove atom 1132 need remove atom 1134 need remove atom 1136 need remove atom 1138 [I] need removenum=42 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2224765447e-01 9.9993764848e-02 3.4968165306e-01 dub= -3.3804077796e-01 9.9996160131e-02 -3.5634807738e-01 duc= 3.2224855904e-01 -1.0131667078e-01 3.4968277203e-01 [I] overall tilt ux= -6.6028843243e-01 2.3952830663e-06 -7.0602973044e-01 uy= 9.0456236940e-07 -2.0131043563e-01 1.1189759368e-06 [I] storedr[0]= -1.2095649310e-02 4.9195265179e-02 1.1609552437e-02 [I] originally 2184 atoms [I] insert 0 atoms [I] now 2184 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 2184 atoms [I] remove 42 atoms [I] now 2142 atoms [I] NP=2142 n=2142 [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 24 14:29:19 2023 CPU time spent: 0.176260 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 24 14:30:18 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [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) = 15 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 3 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=2520 [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.035714285714286 [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.28133100175789 [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.5714285714e-02 0.0000000000e+00 1.6666666667e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=4.431280e+01 pxpb=0.000000e+00 pypb=6.714403e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 1309 need remove atom 1311 need remove atom 1313 need remove atom 1315 need remove atom 1317 need remove atom 1319 need remove atom 1321 need remove atom 1323 need remove atom 1325 need remove atom 1327 need remove atom 1329 need remove atom 1331 need remove atom 1333 need remove atom 1335 need remove atom 1337 need remove atom 1339 need remove atom 1341 need remove atom 1343 need remove atom 1345 need remove atom 1347 need remove atom 1349 need remove atom 1351 need remove atom 1353 need remove atom 1355 need remove atom 1357 need remove atom 1359 need remove atom 1361 need remove atom 1363 need remove atom 1365 need remove atom 1367 need remove atom 1369 need remove atom 1371 need remove atom 1373 need remove atom 1375 need remove atom 1377 need remove atom 1379 need remove atom 1381 need remove atom 1383 need remove atom 1385 need remove atom 1387 need remove atom 1389 need remove atom 1391 need remove atom 1393 need remove atom 1397 need remove atom 1401 [I] need removenum=45 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2183822978e-01 1.0013630142e-01 3.4890125272e-01 dub= -3.3758942528e-01 1.0013864822e-01 -3.5545227396e-01 duc= 3.2183913812e-01 -1.0144896958e-01 3.4890237050e-01 [I] overall tilt ux= -6.5942765506e-01 2.3468070244e-06 -7.0435352668e-01 uy= 9.0833669264e-07 -2.0158527099e-01 1.1177735824e-06 [I] storedr[0]= -1.1808723122e-02 4.5614307595e-02 1.0678103738e-02 [I] originally 2520 atoms [I] insert 0 atoms [I] now 2520 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 2520 atoms [I] remove 45 atoms [I] now 2475 atoms [I] NP=2475 n=2475 [I] ASSIGN finalcnfile = dipole_1.4.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.4.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.4.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 24 14:30:19 2023 CPU time spent: 0.209975 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 24 14:31:50 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [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) = 16 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 4 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=3840 [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.125 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.2500000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=4.747800e+01 pxpb=0.000000e+00 pypb=7.162030e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 1856 need remove atom 1860 need remove atom 1864 need remove atom 1868 need remove atom 1872 need remove atom 1874 need remove atom 1876 need remove atom 1878 need remove atom 1880 need remove atom 1882 need remove atom 1884 need remove atom 1886 need remove atom 1888 need remove atom 1890 need remove atom 1892 need remove atom 1894 need remove atom 1896 need remove atom 1898 need remove atom 1900 need remove atom 1902 need remove atom 1904 need remove atom 1906 need remove atom 1908 need remove atom 1910 need remove atom 1912 need remove atom 1914 need remove atom 1916 need remove atom 1918 need remove atom 1920 need remove atom 1922 need remove atom 1924 need remove atom 1926 need remove atom 1928 need remove atom 1930 need remove atom 1932 need remove atom 1934 need remove atom 1936 need remove atom 1938 need remove atom 1940 need remove atom 1942 need remove atom 1944 need remove atom 1946 need remove atom 1948 need remove atom 1950 need remove atom 1952 need remove atom 1954 need remove atom 1956 need remove atom 1958 need remove atom 1960 need remove atom 1962 need remove atom 1964 need remove atom 1966 need remove atom 1968 need remove atom 1970 need remove atom 1972 need remove atom 1974 need remove atom 1976 need remove atom 1978 need remove atom 1980 need remove atom 1982 need remove atom 1986 need remove atom 1990 need remove atom 1994 need remove atom 1998 [I] need removenum=64 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2147992018e-01 1.0026040114e-01 3.4822046000e-01 dub= -3.3719383798e-01 1.0026270511e-01 -3.5467065198e-01 duc= 3.2148083178e-01 -1.0156414759e-01 3.4822157668e-01 [I] overall tilt ux= -6.5867375816e-01 2.3039670258e-06 -7.0289111198e-01 uy= 9.1159870003e-07 -2.0182454873e-01 1.1166873630e-06 [I] storedr[0]= -1.1551536733e-02 4.2510132224e-02 9.8655733282e-03 [I] originally 3840 atoms [I] insert 0 atoms [I] now 3840 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 3840 atoms [I] remove 64 atoms [I] now 3776 atoms [I] NP=3776 n=3776 [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 24 14:31:51 2023 CPU time spent: 0.317493 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 24 14:32:55 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 16 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 17 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 4 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=4352 [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.03125 [I] ASSIGN mkdipole(3) = 0 [I] ASSIGN mkdipole(4) = 0.125 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.1250000000e-02 0.0000000000e+00 1.2500000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=5.064320e+01 pxpb=0.000000e+00 pypb=7.609657e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 2243 need remove atom 2247 need remove atom 2251 need remove atom 2255 need remove atom 2257 need remove atom 2259 need remove atom 2261 need remove atom 2263 need remove atom 2265 need remove atom 2267 need remove atom 2269 need remove atom 2271 need remove atom 2273 need remove atom 2275 need remove atom 2277 need remove atom 2279 need remove atom 2281 need remove atom 2283 need remove atom 2285 need remove atom 2287 need remove atom 2289 need remove atom 2291 need remove atom 2293 need remove atom 2295 need remove atom 2297 need remove atom 2299 need remove atom 2301 need remove atom 2303 need remove atom 2305 need remove atom 2307 need remove atom 2309 need remove atom 2311 need remove atom 2313 need remove atom 2315 need remove atom 2317 need remove atom 2319 need remove atom 2321 need remove atom 2323 need remove atom 2325 need remove atom 2327 need remove atom 2329 need remove atom 2331 need remove atom 2333 need remove atom 2335 need remove atom 2337 need remove atom 2339 need remove atom 2341 need remove atom 2343 need remove atom 2345 need remove atom 2347 need remove atom 2349 need remove atom 2351 need remove atom 2353 need remove atom 2355 need remove atom 2357 need remove atom 2359 need remove atom 2361 need remove atom 2363 need remove atom 2365 need remove atom 2367 need remove atom 2369 need remove atom 2371 need remove atom 2373 need remove atom 2375 need remove atom 2377 need remove atom 2379 need remove atom 2381 need remove atom 2383 [I] need removenum=68 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2116371996e-01 1.0036941710e-01 3.4762135172e-01 dub= -3.3684429245e-01 1.0037168295e-01 -3.5398269566e-01 duc= 3.2116463441e-01 -1.0166531962e-01 3.4762246743e-01 [I] overall tilt ux= -6.5800801241e-01 2.2658451815e-06 -7.0160404739e-01 uy= 9.1444543637e-07 -2.0203473672e-01 1.1157028118e-06 [I] storedr[0]= -1.1320507376e-02 3.9794031618e-02 9.1499832221e-03 [I] originally 4352 atoms [I] insert 0 atoms [I] now 4352 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 4352 atoms [I] remove 68 atoms [I] now 4284 atoms [I] NP=4284 n=4284 [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 24 14:32:56 2023 CPU time spent: 0.372135 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 24 14:33:53 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 17 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 18 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 4 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=4896 [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.029411764705882 [I] ASSIGN mkdipole(3) = 0 [I] ASSIGN mkdipole(4) = 0.125 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.9411764706e-02 0.0000000000e+00 1.2500000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=5.380840e+01 pxpb=0.000000e+00 pypb=8.057284e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 2384 need remove atom 2386 need remove atom 2388 need remove atom 2390 need remove atom 2392 need remove atom 2394 need remove atom 2396 need remove atom 2398 need remove atom 2400 need remove atom 2402 need remove atom 2404 need remove atom 2406 need remove atom 2408 need remove atom 2410 need remove atom 2412 need remove atom 2414 need remove atom 2416 need remove atom 2418 need remove atom 2420 need remove atom 2422 need remove atom 2424 need remove atom 2426 need remove atom 2428 need remove atom 2430 need remove atom 2432 need remove atom 2434 need remove atom 2436 need remove atom 2438 need remove atom 2440 need remove atom 2442 need remove atom 2444 need remove atom 2446 need remove atom 2448 need remove atom 2450 need remove atom 2452 need remove atom 2454 need remove atom 2456 need remove atom 2458 need remove atom 2460 need remove atom 2462 need remove atom 2464 need remove atom 2466 need remove atom 2468 need remove atom 2470 need remove atom 2472 need remove atom 2474 need remove atom 2476 need remove atom 2478 need remove atom 2480 need remove atom 2482 need remove atom 2484 need remove atom 2486 need remove atom 2488 need remove atom 2490 need remove atom 2492 need remove atom 2494 need remove atom 2496 need remove atom 2498 need remove atom 2500 need remove atom 2502 need remove atom 2504 need remove atom 2506 need remove atom 2508 need remove atom 2510 need remove atom 2512 need remove atom 2514 need remove atom 2516 need remove atom 2518 need remove atom 2520 need remove atom 2522 need remove atom 2524 need remove atom 2526 [I] need removenum=72 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2088262077e-01 1.0046593709e-01 3.4709006086e-01 dub= -3.3653319954e-01 1.0046816881e-01 -3.5337252500e-01 duc= 3.2088353772e-01 -1.0175488941e-01 3.4709117567e-01 [I] overall tilt ux= -6.5741582031e-01 2.2317102131e-06 -7.0046258586e-01 uy= 9.1695101273e-07 -2.0222082650e-01 1.1148073390e-06 [I] storedr[0]= -1.1112373315e-02 3.7397911168e-02 8.5145716471e-03 [I] originally 4896 atoms [I] insert 0 atoms [I] now 4896 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 4896 atoms [I] remove 72 atoms [I] now 4824 atoms [I] NP=4824 n=4824 [I] ASSIGN finalcnfile = dipole_1.7.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.7.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.7.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 24 14:33:53 2023 CPU time spent: 0.413397 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 24 14:36:00 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 18 [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) = 0 [I] ASSIGN makecnspec(11) = 4 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=5472 [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.125 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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 1.2500000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=5.697360e+01 pxpb=0.000000e+00 pypb=8.504911e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 2817 need remove atom 2819 need remove atom 2821 need remove atom 2823 need remove atom 2825 need remove atom 2827 need remove atom 2829 need remove atom 2831 need remove atom 2833 need remove atom 2835 need remove atom 2837 need remove atom 2839 need remove atom 2841 need remove atom 2843 need remove atom 2845 need remove atom 2847 need remove atom 2849 need remove atom 2851 need remove atom 2853 need remove atom 2855 need remove atom 2857 need remove atom 2859 need remove atom 2861 need remove atom 2863 need remove atom 2865 need remove atom 2867 need remove atom 2869 need remove atom 2871 need remove atom 2873 need remove atom 2875 need remove atom 2877 need remove atom 2879 need remove atom 2881 need remove atom 2883 need remove atom 2885 need remove atom 2887 need remove atom 2889 need remove atom 2891 need remove atom 2893 need remove atom 2895 need remove atom 2897 need remove atom 2899 need remove atom 2901 need remove atom 2903 need remove atom 2905 need remove atom 2907 need remove atom 2909 need remove atom 2911 need remove atom 2913 need remove atom 2915 need remove atom 2917 need remove atom 2919 need remove atom 2921 need remove atom 2923 need remove atom 2925 need remove atom 2927 need remove atom 2929 need remove atom 2931 need remove atom 2933 need remove atom 2935 need remove atom 2937 need remove atom 2939 need remove atom 2941 need remove atom 2943 need remove atom 2945 need remove atom 2947 need remove atom 2949 need remove atom 2951 need remove atom 2953 need remove atom 2955 need remove atom 2957 need remove atom 2959 need remove atom 2961 need remove atom 2965 need remove atom 2969 need remove atom 2973 [I] need removenum=76 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2063108718e-01 1.0055198985e-01 3.4661569066e-01 dub= -3.3625454980e-01 1.0055419082e-01 -3.5282765675e-01 duc= 3.2063200636e-01 -1.0183474126e-01 3.4661680465e-01 [I] overall tilt ux= -6.5688563699e-01 2.2009729784e-06 -6.9944334741e-01 uy= 9.1917296818e-07 -2.0238673112e-01 1.1139900923e-06 [I] storedr[0]= -1.0924256623e-02 3.5268654975e-02 7.9463062421e-03 [I] originally 5472 atoms [I] insert 0 atoms [I] now 5472 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 5472 atoms [I] remove 76 atoms [I] now 5396 atoms [I] NP=5396 n=5396 [I] ASSIGN finalcnfile = dipole_1.8.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_1.8.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_1.8.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 24 14:36:01 2023 CPU time spent: 0.463770 s gzip: dipole_1.8.lammps: file size changed while zipping -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- [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 24 14:36:02 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 19 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 20 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 5 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=7600 [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.026315789473684 [I] ASSIGN mkdipole(3) = 0 [I] ASSIGN mkdipole(4) = 0.1 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.6315789474e-02 0.0000000000e+00 1.0000000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=6.013880e+01 pxpb=0.000000e+00 pypb=8.952537e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 3700 need remove atom 3704 need remove atom 3708 need remove atom 3712 need remove atom 3716 need remove atom 3720 need remove atom 3722 need remove atom 3724 need remove atom 3726 need remove atom 3728 need remove atom 3730 need remove atom 3732 need remove atom 3734 need remove atom 3736 need remove atom 3738 need remove atom 3740 need remove atom 3742 need remove atom 3744 need remove atom 3746 need remove atom 3748 need remove atom 3750 need remove atom 3752 need remove atom 3754 need remove atom 3756 need remove atom 3758 need remove atom 3760 need remove atom 3762 need remove atom 3764 need remove atom 3766 need remove atom 3768 need remove atom 3770 need remove atom 3772 need remove atom 3774 need remove atom 3776 need remove atom 3778 need remove atom 3780 need remove atom 3782 need remove atom 3784 need remove atom 3786 need remove atom 3788 need remove atom 3790 need remove atom 3792 need remove atom 3794 need remove atom 3796 need remove atom 3798 need remove atom 3800 need remove atom 3802 need remove atom 3804 need remove atom 3806 need remove atom 3808 need remove atom 3810 need remove atom 3812 need remove atom 3814 need remove atom 3816 need remove atom 3818 need remove atom 3820 need remove atom 3822 need remove atom 3824 need remove atom 3826 need remove atom 3828 need remove atom 3830 need remove atom 3832 need remove atom 3834 need remove atom 3836 need remove atom 3838 need remove atom 3840 need remove atom 3842 need remove atom 3844 need remove atom 3846 need remove atom 3848 need remove atom 3850 need remove atom 3852 need remove atom 3854 need remove atom 3856 need remove atom 3858 need remove atom 3860 need remove atom 3862 need remove atom 3864 need remove atom 3866 need remove atom 3868 need remove atom 3870 need remove atom 3872 need remove atom 3874 need remove atom 3876 need remove atom 3878 need remove atom 3880 need remove atom 3882 need remove atom 3884 need remove atom 3886 need remove atom 3888 need remove atom 3890 need remove atom 3892 need remove atom 3894 need remove atom 3896 need remove atom 3898 need remove atom 3902 need remove atom 3906 need remove atom 3910 need remove atom 3914 need remove atom 3918 [I] need removenum=100 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2040468919e-01 1.0062918810e-01 3.4618955990e-01 dub= -3.3600352313e-01 1.0063136125e-01 -3.5233814255e-01 duc= 3.2040561035e-01 -1.0190637299e-01 3.4619067314e-01 [I] overall tilt ux= -6.5640821232e-01 2.1731539245e-06 -6.9852770245e-01 uy= 9.2115665307e-07 -2.0253556109e-01 1.1132417413e-06 [I] storedr[0]= -1.0753657636e-02 3.3364237619e-02 7.4348848030e-03 [I] originally 7600 atoms [I] insert 0 atoms [I] now 7600 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 7600 atoms [I] remove 100 atoms [I] now 7500 atoms [I] NP=7500 n=7500 [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 24 14:36:02 2023 CPU time spent: 0.638079 s gzip: dipole_1.9.lammps: file size changed while zipping -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- [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 24 14:36:03 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 20 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 21 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 5 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=8400 [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.025 [I] ASSIGN mkdipole(3) = 0 [I] ASSIGN mkdipole(4) = 0.1 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.5000000000e-02 0.0000000000e+00 1.0000000000e-01 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=6.330400e+01 pxpb=0.000000e+00 pypb=9.400164e+01 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 4303 need remove atom 4307 need remove atom 4311 need remove atom 4315 need remove atom 4319 need remove atom 4321 need remove atom 4323 need remove atom 4325 need remove atom 4327 need remove atom 4329 need remove atom 4331 need remove atom 4333 need remove atom 4335 need remove atom 4337 need remove atom 4339 need remove atom 4341 need remove atom 4343 need remove atom 4345 need remove atom 4347 need remove atom 4349 need remove atom 4351 need remove atom 4353 need remove atom 4355 need remove atom 4357 need remove atom 4359 need remove atom 4361 need remove atom 4363 need remove atom 4365 need remove atom 4367 need remove atom 4369 need remove atom 4371 need remove atom 4373 need remove atom 4375 need remove atom 4377 need remove atom 4379 need remove atom 4381 need remove atom 4383 need remove atom 4385 need remove atom 4387 need remove atom 4389 need remove atom 4391 need remove atom 4393 need remove atom 4395 need remove atom 4397 need remove atom 4399 need remove atom 4401 need remove atom 4403 need remove atom 4405 need remove atom 4407 need remove atom 4409 need remove atom 4411 need remove atom 4413 need remove atom 4415 need remove atom 4417 need remove atom 4419 need remove atom 4421 need remove atom 4423 need remove atom 4425 need remove atom 4427 need remove atom 4429 need remove atom 4431 need remove atom 4433 need remove atom 4435 need remove atom 4437 need remove atom 4439 need remove atom 4441 need remove atom 4443 need remove atom 4445 need remove atom 4447 need remove atom 4449 need remove atom 4451 need remove atom 4453 need remove atom 4455 need remove atom 4457 need remove atom 4459 need remove atom 4461 need remove atom 4463 need remove atom 4465 need remove atom 4467 need remove atom 4469 need remove atom 4471 need remove atom 4473 need remove atom 4475 need remove atom 4477 need remove atom 4479 need remove atom 4481 need remove atom 4483 need remove atom 4485 need remove atom 4487 need remove atom 4489 need remove atom 4491 need remove atom 4493 need remove atom 4495 need remove atom 4497 need remove atom 4499 need remove atom 4501 need remove atom 4503 need remove atom 4505 need remove atom 4507 need remove atom 4509 need remove atom 4511 need remove atom 4513 need remove atom 4515 need remove atom 4517 need remove atom 4519 [I] need removenum=105 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2019983947e-01 1.0069882984e-01 3.4580466711e-01 dub= -3.3577620855e-01 1.0070097770e-01 -3.5189595668e-01 duc= 3.2020076241e-01 -1.0197099004e-01 3.4580577967e-01 [I] overall tilt ux= -6.5597604802e-01 2.1478589916e-06 -6.9770062379e-01 uy= 9.2293827214e-07 -2.0266981988e-01 1.1125542890e-06 [I] storedr[0]= -1.0598421747e-02 3.1650983439e-02 6.9720450219e-03 [I] originally 8400 atoms [I] insert 0 atoms [I] now 8400 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 8400 atoms [I] remove 105 atoms [I] now 8295 atoms [I] NP=8295 n=8295 [I] ASSIGN finalcnfile = dipole_2.0.lammps [I] EXEC writeLAMMPS [I] FILEOPEN dipole_2.0.lammps data file MD++ -> LAMMPS dump file LAMMPS --> MD++ [I] WRITEFILE dipole_2.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 24 14:36:04 2023 CPU time spent: 0.702126 s gzip: dipole_2.0.lammps: file size changed while zipping -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- [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 24 14:36:04 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 25 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 27 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 6 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=16200 [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.02 [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.28133100175789 [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.0000000000e-02 0.0000000000e+00 8.3333333333e-02 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=7.913000e+01 pxpb=0.000000e+00 pypb=1.208593e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms [I] need removenum=0 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2247507687e-01 9.9914253576e-02 3.5011629486e-01 dub= -3.3829117963e-01 9.9916675563e-02 -3.5684690776e-01 duc= 3.2247597931e-01 -1.0124286631e-01 3.5011741448e-01 [I] overall tilt ux= -6.6076625649e-01 2.4219871667e-06 -7.0696320262e-01 uy= 9.0244424250e-07 -2.0115711988e-01 1.1196256924e-06 [I] storedr[0]= -1.0372794534e-02 2.5387761721e-02 4.1431642316e-03 [I] originally 16200 atoms [I] insert 0 atoms [I] now 16200 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 16200 atoms [I] remove 0 atoms [I] now 16200 atoms [I] NP=16200 n=16200 [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 24 14:36:09 2023 CPU time spent: 1.368306 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 24 14:36:27 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 30 [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) = 0 [I] ASSIGN makecnspec(11) = 7 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=26880 [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.071428571428571 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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 7.1428571429e-02 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=9.495600e+01 pxpb=0.000000e+00 pypb=1.432406e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 12768 need remove atom 12772 need remove atom 12776 need remove atom 12780 need remove atom 12784 need remove atom 12788 need remove atom 12792 need remove atom 12796 need remove atom 12798 need remove atom 12800 need remove atom 12802 need remove atom 12804 need remove atom 12806 need remove atom 12808 need remove atom 12810 need remove atom 12812 need remove atom 12814 need remove atom 12816 need remove atom 12818 need remove atom 12820 need remove atom 12822 need remove atom 12824 need remove atom 12826 need remove atom 12828 need remove atom 12830 need remove atom 12832 need remove atom 12834 need remove atom 12836 need remove atom 12838 need remove atom 12840 need remove atom 12842 need remove atom 12844 need remove atom 12846 need remove atom 12848 need remove atom 12850 need remove atom 12852 need remove atom 12854 need remove atom 12856 need remove atom 12858 need remove atom 12860 need remove atom 12862 need remove atom 12864 need remove atom 12866 need remove atom 12868 need remove atom 12870 need remove atom 12872 need remove atom 12874 need remove atom 12876 need remove atom 12878 need remove atom 12880 need remove atom 12882 need remove atom 12884 need remove atom 12886 need remove atom 12888 need remove atom 12890 need remove atom 12892 need remove atom 12894 need remove atom 12896 need remove atom 12898 need remove atom 12900 need remove atom 12902 need remove atom 12904 need remove atom 12906 need remove atom 12908 need remove atom 12910 need remove atom 12912 need remove atom 12914 need remove atom 12916 need remove atom 12918 need remove atom 12920 need remove atom 12922 need remove atom 12924 need remove atom 12926 need remove atom 12928 need remove atom 12930 need remove atom 12932 need remove atom 12934 need remove atom 12936 need remove atom 12938 need remove atom 12940 need remove atom 12942 need remove atom 12944 need remove atom 12946 need remove atom 12948 need remove atom 12950 need remove atom 12952 need remove atom 12954 need remove atom 12956 need remove atom 12958 need remove atom 12960 need remove atom 12962 need remove atom 12964 need remove atom 12966 need remove atom 12968 need remove atom 12970 need remove atom 12972 need remove atom 12974 need remove atom 12976 need remove atom 12978 need remove atom 12980 need remove atom 12982 need remove atom 12984 need remove atom 12986 need remove atom 12988 need remove atom 12990 need remove atom 12992 need remove atom 12994 need remove atom 12996 need remove atom 12998 need remove atom 13000 need remove atom 13002 need remove atom 13004 need remove atom 13006 need remove atom 13008 need remove atom 13010 need remove atom 13012 need remove atom 13014 need remove atom 13016 need remove atom 13018 need remove atom 13020 need remove atom 13022 need remove atom 13024 need remove atom 13026 need remove atom 13028 need remove atom 13030 need remove atom 13032 need remove atom 13034 need remove atom 13036 need remove atom 13038 need remove atom 13040 need remove atom 13042 need remove atom 13044 need remove atom 13046 need remove atom 13048 need remove atom 13050 need remove atom 13052 need remove atom 13054 need remove atom 13056 need remove atom 13058 need remove atom 13060 need remove atom 13062 need remove atom 13064 need remove atom 13066 need remove atom 13068 need remove atom 13070 need remove atom 13072 need remove atom 13074 need remove atom 13076 need remove atom 13078 need remove atom 13080 need remove atom 13082 need remove atom 13084 need remove atom 13086 need remove atom 13088 need remove atom 13090 need remove atom 13092 need remove atom 13094 need remove atom 13096 need remove atom 13098 need remove atom 13100 need remove atom 13102 need remove atom 13104 need remove atom 13106 need remove atom 13108 need remove atom 13110 need remove atom 13112 need remove atom 13114 need remove atom 13116 need remove atom 13118 need remove atom 13120 need remove atom 13122 need remove atom 13124 need remove atom 13126 need remove atom 13128 need remove atom 13130 need remove atom 13132 need remove atom 13134 need remove atom 13136 need remove atom 13138 need remove atom 13140 need remove atom 13142 need remove atom 13144 need remove atom 13146 need remove atom 13148 need remove atom 13150 need remove atom 13152 need remove atom 13154 need remove atom 13156 need remove atom 13158 need remove atom 13160 need remove atom 13162 need remove atom 13164 need remove atom 13166 need remove atom 13168 need remove atom 13170 need remove atom 13172 need remove atom 13174 need remove atom 13176 need remove atom 13178 need remove atom 13180 need remove atom 13182 need remove atom 13184 need remove atom 13186 need remove atom 13188 need remove atom 13190 need remove atom 13192 need remove atom 13194 need remove atom 13196 need remove atom 13198 need remove atom 13200 need remove atom 13202 need remove atom 13204 need remove atom 13206 need remove atom 13208 need remove atom 13210 need remove atom 13212 need remove atom 13214 need remove atom 13218 need remove atom 13222 need remove atom 13226 need remove atom 13230 need remove atom 13234 need remove atom 13238 need remove atom 13242 [I] need removenum=224 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2147992018e-01 1.0026040114e-01 3.4822046000e-01 dub= -3.3719383798e-01 1.0026270511e-01 -3.5467065198e-01 duc= 3.2148083178e-01 -1.0156414759e-01 3.4822157668e-01 [I] overall tilt ux= -6.5867375816e-01 2.3039670248e-06 -7.0289111198e-01 uy= 9.1159870003e-07 -2.0182454873e-01 1.1166873639e-06 [I] storedr[0]= -9.8846743802e-03 2.0998489968e-02 3.1778138526e-03 [I] originally 26880 atoms [I] insert 0 atoms [I] now 26880 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 26880 atoms [I] remove 224 atoms [I] now 26656 atoms [I] NP=26656 n=26656 [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 24 14:36:30 2023 CPU time spent: 2.237314 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 24 14:38:38 2023 [I] ASSIGN latticestructure = body-centered-cubic [I] ASSIGN latticeconst = 3.1651999652386 [I] ASSIGN makecnspec = 0 [I] ASSIGN makecnspec(1) = 0 [I] ASSIGN makecnspec(2) = 1 [I] ASSIGN makecnspec(3) = 35 [I] ASSIGN makecnspec(4) = -1 [I] ASSIGN makecnspec(5) = 1 [I] ASSIGN makecnspec(6) = 0 [I] ASSIGN makecnspec(7) = 37 [I] ASSIGN makecnspec(8) = 1 [I] ASSIGN makecnspec(9) = 1 [I] ASSIGN makecnspec(10) = 0 [I] ASSIGN makecnspec(11) = 8 [I] EXEC makecn [I] makecrystal [I] my_unitcell h=[ 0 -1 1 0 1 1 1 0 0]; basis=[ 0.5 0.5 0 0 0 0 0 0 0 1 0.5 0 0.5 0 2 0 0.5 0.5 0 3 ]; [I] makecn: _NP=41440 [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.014285714285714 [I] ASSIGN mkdipole(3) = 0 [I] ASSIGN mkdipole(4) = 0.0625 [I] ASSIGN mkdipole(5) = -0.01 [I] ASSIGN mkdipole(6) = -0.2499 [I] ASSIGN mkdipole(7) = 0.251 [I] ASSIGN mkdipole(8) = 0.28133100175789 [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.4285714286e-02 0.0000000000e+00 6.2500000000e-02 b= 1.5825999826e+00 0.0000000000e+00 2.2381343592e+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.5825999826193 by=0 pxpa=1.107820e+02 pxpb=0.000000e+00 pypb=1.656219e+02 pypa=0.000000e+00 [I] allocate storedr [I] Need to remove atoms need remove atom 20419 need remove atom 20423 need remove atom 20427 need remove atom 20431 need remove atom 20435 need remove atom 20439 need remove atom 20443 need remove atom 20447 need remove atom 20449 need remove atom 20451 need remove atom 20453 need remove atom 20455 need remove atom 20457 need remove atom 20459 need remove atom 20461 need remove atom 20463 need remove atom 20465 need remove atom 20467 need remove atom 20469 need remove atom 20471 need remove atom 20473 need remove atom 20475 need remove atom 20477 need remove atom 20479 need remove atom 20481 need remove atom 20483 need remove atom 20485 need remove atom 20487 need remove atom 20489 need remove atom 20491 need remove atom 20493 need remove atom 20495 need remove atom 20497 need remove atom 20499 need remove atom 20501 need remove atom 20503 need remove atom 20505 need remove atom 20507 need remove atom 20509 need remove atom 20511 need remove atom 20513 need remove atom 20515 need remove atom 20517 need remove atom 20519 need remove atom 20521 need remove atom 20523 need remove atom 20525 need remove atom 20527 need remove atom 20529 need remove atom 20531 need remove atom 20533 need remove atom 20535 need remove atom 20537 need remove atom 20539 need remove atom 20541 need remove atom 20543 need remove atom 20545 need remove atom 20547 need remove atom 20549 need remove atom 20551 need remove atom 20553 need remove atom 20555 need remove atom 20557 need remove atom 20559 need remove atom 20561 need remove atom 20563 need remove atom 20565 need remove atom 20567 need remove atom 20569 need remove atom 20571 need remove atom 20573 need remove atom 20575 need remove atom 20577 need remove atom 20579 need remove atom 20581 need remove atom 20583 need remove atom 20585 need remove atom 20587 need remove atom 20589 need remove atom 20591 need remove atom 20593 need remove atom 20595 need remove atom 20597 need remove atom 20599 need remove atom 20601 need remove atom 20603 need remove atom 20605 need remove atom 20607 need remove atom 20609 need remove atom 20611 need remove atom 20613 need remove atom 20615 need remove atom 20617 need remove atom 20619 need remove atom 20621 need remove atom 20623 need remove atom 20625 need remove atom 20627 need remove atom 20629 need remove atom 20631 need remove atom 20633 need remove atom 20635 need remove atom 20637 need remove atom 20639 need remove atom 20641 need remove atom 20643 need remove atom 20645 need remove atom 20647 need remove atom 20649 need remove atom 20651 need remove atom 20653 need remove atom 20655 need remove atom 20657 need remove atom 20659 need remove atom 20661 need remove atom 20663 need remove atom 20665 need remove atom 20667 need remove atom 20669 need remove atom 20671 need remove atom 20673 need remove atom 20675 need remove atom 20677 need remove atom 20679 need remove atom 20681 need remove atom 20683 need remove atom 20685 need remove atom 20687 need remove atom 20689 need remove atom 20691 need remove atom 20693 need remove atom 20695 need remove atom 20697 need remove atom 20699 need remove atom 20701 need remove atom 20703 need remove atom 20705 need remove atom 20707 need remove atom 20709 need remove atom 20711 need remove atom 20713 need remove atom 20715 need remove atom 20717 need remove atom 20719 need remove atom 20721 need remove atom 20723 need remove atom 20725 need remove atom 20727 need remove atom 20729 need remove atom 20731 need remove atom 20733 need remove atom 20735 need remove atom 20737 need remove atom 20739 need remove atom 20741 need remove atom 20743 need remove atom 20745 need remove atom 20747 need remove atom 20749 need remove atom 20751 need remove atom 20753 need remove atom 20755 need remove atom 20757 need remove atom 20759 need remove atom 20761 need remove atom 20763 need remove atom 20765 need remove atom 20767 need remove atom 20769 need remove atom 20771 need remove atom 20773 need remove atom 20775 need remove atom 20777 need remove atom 20779 need remove atom 20781 need remove atom 20783 need remove atom 20785 need remove atom 20787 need remove atom 20789 need remove atom 20791 need remove atom 20793 need remove atom 20795 need remove atom 20797 need remove atom 20799 need remove atom 20801 need remove atom 20803 need remove atom 20805 need remove atom 20807 need remove atom 20809 need remove atom 20811 need remove atom 20813 need remove atom 20815 need remove atom 20817 need remove atom 20819 need remove atom 20821 need remove atom 20823 need remove atom 20825 need remove atom 20827 need remove atom 20829 need remove atom 20831 need remove atom 20833 need remove atom 20835 need remove atom 20837 need remove atom 20839 need remove atom 20841 need remove atom 20843 need remove atom 20845 need remove atom 20847 need remove atom 20849 need remove atom 20851 need remove atom 20853 need remove atom 20855 need remove atom 20857 need remove atom 20859 need remove atom 20861 need remove atom 20863 need remove atom 20865 need remove atom 20867 need remove atom 20869 need remove atom 20871 need remove atom 20873 need remove atom 20875 need remove atom 20877 need remove atom 20879 need remove atom 20881 need remove atom 20883 need remove atom 20885 need remove atom 20887 need remove atom 20889 need remove atom 20891 need remove atom 20893 need remove atom 20895 need remove atom 20897 need remove atom 20899 need remove atom 20901 need remove atom 20903 need remove atom 20905 need remove atom 20907 need remove atom 20909 need remove atom 20911 need remove atom 20913 need remove atom 20915 need remove atom 20917 need remove atom 20919 need remove atom 20921 need remove atom 20923 need remove atom 20925 need remove atom 20927 need remove atom 20929 need remove atom 20931 need remove atom 20933 need remove atom 20935 need remove atom 20937 need remove atom 20939 need remove atom 20941 need remove atom 20943 need remove atom 20945 need remove atom 20947 need remove atom 20949 need remove atom 20951 need remove atom 20953 need remove atom 20955 need remove atom 20957 need remove atom 20959 need remove atom 20961 need remove atom 20963 need remove atom 20965 need remove atom 20967 need remove atom 20969 need remove atom 20971 need remove atom 20973 need remove atom 20975 need remove atom 20977 need remove atom 20979 need remove atom 20981 need remove atom 20983 need remove atom 20985 need remove atom 20987 need remove atom 20989 need remove atom 20991 need remove atom 20993 need remove atom 20995 need remove atom 20997 need remove atom 20999 need remove atom 21001 need remove atom 21003 need remove atom 21005 need remove atom 21007 need remove atom 21009 need remove atom 21011 need remove atom 21013 need remove atom 21015 need remove atom 21017 need remove atom 21019 need remove atom 21021 need remove atom 21023 [I] need removenum=296 minx = -10 maxx = 10 miny = -10 maxy = 10 [I] dua= 3.2075345757e-01 1.0051016268e-01 3.4684634804e-01 dub= -3.3639014469e-01 1.0051237863e-01 -3.5309260110e-01 duc= 3.2075437566e-01 -1.0179592866e-01 3.4684746243e-01 [I] overall tilt ux= -6.5714360225e-01 2.2159491027e-06 -6.9993894913e-01 uy= 9.1809436370e-07 -2.0230609135e-01 1.1143895524e-06 [I] storedr[0]= -9.5366048495e-03 1.7873370226e-02 2.4578185147e-03 [I] originally 41440 atoms [I] insert 0 atoms [I] now 41440 atoms [I] shiftbox to accommodate plastic strain makedipole: store=0 [I] free storedr [I] originally 41440 atoms [I] remove 296 atoms [I] now 41144 atoms [I] NP=41144 n=41144 [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 24 14:38:42 2023 CPU time spent: 3.512334 s {"realtime":727.51,"usertime":1267.83,"systime":10.08,"memmax":138244,"memavg":0}