From 44292e02a3fd40a70832a75a2fe21916850ee39c Mon Sep 17 00:00:00 2001 From: Daniel Pozsar Date: Fri, 8 Nov 2024 14:40:21 +0100 Subject: [PATCH] deleted old tests and run pre commit hook --- Fe3GeTe2_benchmark_on_15k_300eset.log | 283 ----------------- Fe3GeTe2_benchmark_on_15k_300eset.pickle | Bin 7650 -> 0 bytes ...eTe2_benchmark_on_15k_300eset_orb_test.log | 281 ----------------- ...2_benchmark_on_15k_300eset_orb_test.pickle | Bin 7519 -> 0 bytes ...Te2_benchmark_on_15k_300eset_orb_test2.log | 283 ----------------- ..._benchmark_on_15k_300eset_orb_test2.pickle | Bin 7774 -> 0 bytes ...Te2_benchmark_on_15k_300eset_orb_test3.log | 288 ------------------ ..._benchmark_on_15k_300eset_orb_test3.pickle | Bin 7520 -> 0 bytes Fe3GeTe2_final_test.pickle | Bin 7839 -> 0 bytes Fe3GeTe2_notebook.pickle | Bin 7821 -> 0 bytes all_atoms.txt | 88 ------ documentation.md | 60 ---- out with HK = HK + Ef * SK.txt | 239 --------------- src/grogu_magn/core.py | 20 ++ src/grogu_magn/grogu.py | 21 -- src/grogu_magn/io.py | 20 ++ src/grogu_magn/magnetism.py | 20 ++ test.py | 1 - 18 files changed, 60 insertions(+), 1544 deletions(-) delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset.log delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset.pickle delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test.log delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test.pickle delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.log delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.pickle delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test3.log delete mode 100644 Fe3GeTe2_benchmark_on_15k_300eset_orb_test3.pickle delete mode 100644 Fe3GeTe2_final_test.pickle delete mode 100644 Fe3GeTe2_notebook.pickle delete mode 100644 all_atoms.txt delete mode 100644 documentation.md delete mode 100644 out with HK = HK + Ef * SK.txt diff --git a/Fe3GeTe2_benchmark_on_15k_300eset.log b/Fe3GeTe2_benchmark_on_15k_300eset.log deleted file mode 100644 index 675fdb4..0000000 --- a/Fe3GeTe2_benchmark_on_15k_300eset.log +++ /dev/null @@ -1,283 +0,0 @@ -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Setup done. Elapsed time: 1.0880625 s -================================================================================================================================================================ -Hamiltonian and exchange field rotated. Elapsed time: 1.55488075 s -================================================================================================================================================================ -Site and pair dictionaries created. Elapsed time: 1.603311625 s -================================================================================================================================================================ -k set created. Elapsed time: 1.618691125 s -================================================================================================================================================================ -Rotations done perpendicular to quantization axis. Elapsed time: 1.848489583 s -================================================================================================================================================================ -Starting matrix inversions -Total number of k points: 225 -Number of energy samples per k point: 300 -Total number of directions: 3 -Total number of matrix inversions: 202500 -The shape of the Hamiltonian and the Greens function is 84x84=7056 -Memory taken by a single Hamiltonian is: 0.015625 KB -Expected memory usage per matrix inversion: 0.5 KB -Expected memory usage per k point for parallel inversion: 450.0 KB -Expected memory usage on root node: 98.876953125 MB -================================================================================================================================================================ -Calculated Greens functions. Elapsed time: 295.988709083 s -================================================================================================================================================================ -Magnetic entities integrated. -Pairs integrated. -Magnetic parameters calculated. -##################################################################### GROGU OUTPUT ############################################################################# -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Atomic information: ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[atom index]Element(orbitals) x [Ang] y [Ang] z [Ang] Sx Sy Sz Q Lx Ly Lz Jx Jy Jz ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(2) -7.339158738013707e-06 4.149278510690423e-06 11.657585837928032 - -[4]Fe(2) -7.326987662162937e-06 4.158274523275774e-06 8.912422537596708 - -[5]Fe(2) 1.8954667088117545 1.0943913231921656 10.285002698393109 - -================================================================================================================================================================ -Exchange [meV] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Magnetic entity1 Magnetic entity2 [i j k] d [Ang] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(2) [4]Fe(2) [0 0 0] d [Ang] 2.745163300331324 -Isotropic: -39.91775489990159 -DMI: [ 5.22846405e-03 -1.25316789e-05 -2.62478203e-06] -Symmetric-anisotropy: [-2.32809260e+00 1.41154602e-06 -9.40957764e-09 1.41154602e-06 - -1.35443302e+00 5.63570534e-07 -9.40957764e-09 5.63570534e-07 - 3.68252563e+00] -J: [-4.22458475e+01 1.41154602e-06 -9.40957764e-09 1.41154602e-06 - -4.12721879e+01 5.63570534e-07 -9.40957764e-09 5.63570534e-07 - -3.62352293e+01] -Energies for debugging: -array([[-3.52462334e-02, 5.22790048e-06, -5.22902762e-06, - -3.50238913e-02], - [-3.72242252e-02, 1.25410885e-08, -1.25222693e-08, - -3.69712117e-02], - [-4.75204846e-02, -4.03632805e-09, 1.21323601e-09, - -4.75204833e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.03697121, -0.04752048, -0.03524623]) -Test J_xx = E(y,z) = E(z,y) --0.0369712117186039 -0.04752048328815303 - -[3]Fe(2) [5]Fe(2) [0 0 0] d [Ang] 2.5835033632437767 -Isotropic: -64.2035893374817 -DMI: [ 3.34203129e+00 -5.80278227e+00 5.97617073e-04] -Symmetric-anisotropy: [ 0.90424825 0.05359555 -0.07541288 0.05359555 0.5167559 -0.0366015 - -0.07541288 -0.0366015 -1.42100415] -J: [-6.32993411e+01 5.35955501e-02 -7.54128755e-02 5.35955501e-02 - -6.36868334e+01 -3.66014961e-02 -7.54128755e-02 -3.66014961e-02 - -6.56245935e+01] -Energies for debugging: -array([[-6.59077017e-02, 3.37863278e-03, -3.30542979e-03, - -6.63454261e-02], - [-6.53414853e-02, 5.87819515e-03, -5.72736940e-03, - -6.56325743e-02], - [-6.10282408e-02, -5.29979330e-05, -5.41931672e-05, - -6.09661079e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06563257, -0.06102824, -0.0659077 ]) -Test J_xx = E(y,z) = E(z,y) --0.06563257426370037 -0.06096610790875501 - -[4]Fe(2) [5]Fe(2) [0 0 0] d [Ang] 2.583501767937866 -Isotropic: -64.1990652847126 -DMI: [-3.36554719e+00 5.81679250e+00 5.98938514e-04] -Symmetric-anisotropy: [ 0.9024007 0.05359252 0.07068773 0.05359252 0.51528204 0.03236462 - 0.07068773 0.03236462 -1.41768273] -J: [-6.32966646e+01 5.35925196e-02 7.06877343e-02 5.35925196e-02 - -6.36837832e+01 3.23646213e-02 7.06877343e-02 3.23646213e-02 - -6.56167480e+01] -Energies for debugging: -array([[-6.58989484e-02, -3.39791181e-03, 3.33318257e-03, - -6.63392511e-02], - [-6.53345476e-02, -5.88748024e-03, 5.74610477e-03, - -6.56271418e-02], - [-6.10283154e-02, -5.29935811e-05, -5.41914581e-05, - -6.09661874e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06562714, -0.06102832, -0.06589895]) -Test J_xx = E(y,z) = E(z,y) --0.06562714176880605 -0.06096618740895869 - -[3]Fe(2) [5]Fe(2) [-1 -1 0] d [Ang] 2.5834973202859075 -Isotropic: -64.40746869580421 -DMI: [-6.66837514e+00 3.56807199e-04 -5.24962006e-04] -Symmetric-anisotropy: [ 2.81689694e-01 -9.99033947e-05 1.00133116e-04 -9.99033947e-05 - 1.23775856e+00 9.79547273e-02 1.00133116e-04 9.79547273e-02 - -1.51944825e+00] -J: [-6.41257790e+01 -9.99033947e-05 1.00133116e-04 -9.99033947e-05 - -6.31697101e+01 9.79547273e-02 1.00133116e-04 9.79547273e-02 - -6.59269169e+01] -Energies for debugging: -array([[-6.51429577e-02, -6.76632987e-03, 6.57042041e-03, - -6.54030227e-02], - [-6.67108762e-02, -4.56940315e-07, 2.56674082e-07, - -6.71912559e-02], - [-6.09363975e-02, -4.25058611e-07, 6.24865400e-07, - -6.10603021e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06719126, -0.0609364 , -0.06514296]) -Test J_xx = E(y,z) = E(z,y) --0.06719125593394114 -0.061060302070081784 - -[4]Fe(2) [5]Fe(2) [-1 -1 0] d [Ang] 2.583495745338251 -Isotropic: -64.40900602860754 -DMI: [ 6.66852978e+00 -7.33794176e-04 -5.22013428e-04] -Symmetric-anisotropy: [ 2.81831507e-01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - 1.23919957e+00 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -1.52103108e+00] -J: [-6.41271745e+01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - -6.31698065e+01 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -6.59300371e+01] -Energies for debugging: -array([[-6.51430169e-02, 6.76649535e-03, -6.57056420e-03, - -6.54030961e-02], - [-6.67170573e-02, 6.60963771e-07, -8.06624581e-07, - -6.71939352e-02], - [-6.09365168e-02, -4.22426970e-07, 6.21599885e-07, - -6.10604139e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06719394, -0.06093652, -0.06514302]) -Test J_xx = E(y,z) = E(z,y) --0.06719393518677816 -0.06106041385741498 - -[3]Fe(2) [5]Fe(2) [-1 0 0] d [Ang] 2.583541444641373 -Isotropic: -64.18690345424889 -DMI: [ 3.36476076e+00 5.80269172e+00 -7.78587899e-05] -Symmetric-anisotropy: [ 0.9015663 -0.05345052 0.07536033 -0.05345052 0.51805291 -0.03239989 - 0.07536033 -0.03239989 -1.41961921] -J: [-6.32853372e+01 -5.34505248e-02 7.53603270e-02 -5.34505248e-02 - -6.36688505e+01 -3.23998894e-02 7.53603270e-02 -3.23998894e-02 - -6.56065227e+01] -Energies for debugging: -array([[-6.58846583e-02, 3.39716065e-03, -3.33236087e-03, - -6.63246913e-02], - [-6.53283870e-02, -5.87805205e-03, 5.72733140e-03, - -6.56193992e-02], - [-6.10130098e-02, 5.33726660e-05, 5.35283836e-05, - -6.09512751e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.0656194 , -0.06101301, -0.06588466]) -Test J_xx = E(y,z) = E(z,y) --0.06561939924718083 -0.06095127506840987 - -[4]Fe(2) [5]Fe(2) [-1 0 0] d [Ang] 2.5835398672184064 -Isotropic: -64.18749916411404 -DMI: [-3.34143011e+00 -5.81613454e+00 -8.04864431e-05] -Symmetric-anisotropy: [ 0.90467906 -0.05344778 -0.07078134 -0.05344778 0.5154147 0.03662271 - -0.07078134 0.03662271 -1.42009377] -J: [-6.32828201e+01 -5.34477814e-02 -7.07813383e-02 -5.34477814e-02 - -6.36720845e+01 3.66227144e-02 -7.07813383e-02 3.66227144e-02 - -6.56075929e+01] -Energies for debugging: -array([[-6.58936437e-02, -3.37805282e-03, 3.30480739e-03, - -6.63310120e-02], - [-6.53215421e-02, 5.88691588e-03, -5.74535321e-03, - -6.56142142e-02], - [-6.10131569e-02, 5.33672949e-05, 5.35282678e-05, - -6.09514260e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06561421, -0.06101316, -0.06589364]) -Test J_xx = E(y,z) = E(z,y) --0.06561421416763057 -0.06095142603633385 - -[4]Fe(2) [5]Fe(2) [-2 0 0] d [Ang] 5.951322298958084 -Isotropic: 4.59955059881751 -DMI: [ 0.3863217 0.8953163 -0.65738037] -Symmetric-anisotropy: [-0.02046716 -0.0444248 -0.03977157 -0.0444248 0.43165605 -0.07442493 - -0.03977157 -0.07442493 -0.41118888] -J: [ 4.57908344 -0.0444248 -0.03977157 -0.0444248 5.03120664 -0.07442493 - -0.03977157 -0.07442493 4.18836172] -Energies for debugging: -array([[ 0.00472963, 0.00046075, -0.0003119 , 0.00466184], - [ 0.00364709, -0.00085554, 0.00093509, 0.0036577 ], - [ 0.00540057, -0.00061296, 0.00070181, 0.00550046]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.0036577 , 0.00540057, 0.00472963]) -Test J_xx = E(y,z) = E(z,y) -0.003657703314451758 0.005500463555586455 - -[4]Fe(2) [5]Fe(2) [-3 0 0] d [Ang] 9.638732176310562 -Isotropic: -0.23103792475126567 -DMI: [ 0.09123182 0.33044193 -0.36639877] -Symmetric-anisotropy: [-0.08681631 0.00193669 -0.00538302 0.00193669 -0.07518869 -0.03449947 - -0.00538302 -0.03449947 0.162005 ] -J: [-0.31785424 0.00193669 -0.00538302 0.00193669 -0.30622661 -0.03449947 - -0.00538302 -0.03449947 -0.06903292] -Energies for debugging: -array([[ 9.75938878e-06, 1.25731290e-04, -5.67323532e-05, - -1.49392424e-04], - [-1.47825232e-04, -3.25058913e-04, 3.35824948e-04, - -2.23552488e-04], - [-4.63060805e-04, -3.68335462e-04, 3.64462081e-04, - -4.12155989e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-2.23552488e-04, -4.63060805e-04, 9.75938878e-06]) -Test J_xx = E(y,z) = E(z,y) --0.00022355248815234167 -0.00041215598873809247 - -================================================================================================================================================================ -Runtime information: -Total runtime: 295.33235375 s ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Initial setup: 0.09973091699999992 s -Hamiltonian conversion and XC field extraction: 0.467 s -Pair and site datastructure creatrions: 0.048 s -k set cration and distribution: 0.015 s -Rotating XC potential: 0.230 s -Greens function inversion: 294.140 s -Calculate energies and magnetic components: 0.332 s diff --git a/Fe3GeTe2_benchmark_on_15k_300eset.pickle b/Fe3GeTe2_benchmark_on_15k_300eset.pickle deleted file mode 100644 index 9fb6305cb6fdd55e99d763464888c6747ddbcc0c..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 7650 zcmcgx30xD$8s~_D9EyNKQB;Z+LA`%52Xr=XtAz4TWB$$KKTGVPmqo^#^ zDk>ty0}oK2N3q#>JoTvv$f-V&3QAGrcnIG3W|N&1IV8Sof4?vP{dRU{{bgVE*%n~{I3gZRQ zLKa8tBoIdNqMZCtmpSMHbcR!mP#}zEC88oHv5@T)&62nBLZjz*)I zVqOx8F|=hAHX0q>pRgFZ>oKmmOedm;s=1c^Y1;qC(RXgW2pvv?Bq9Ewyxxtktn!_B zgKnTMw~KQA>X(8lFJnC2SD07t6>@gmx~E%RX~f{9vXL>YNC7J0v6*1r@+3S|Jp(8# zi7*CZfXNk&F?3`jgrabsgcZ%?2{&iVmOIx5;-KIBq|3)k+33($1Kxd?6TAkwVyE=W65+A z5|hBSzE73Nd{IQwd3JjHl54ecZ>~W%EwYTtiSf5M9&D$KBTJ@_3Q$ob5x#y}VCVpK zfi|4aTrjr0|J0ba`-1L-YjpC9qk*ynuNFqXI*K(SXecMzqOY)wdlO!dP$* zajZ~*mWV_#i4qc}B*u{#PhtXzi6kbGm`vhI62BvH6^Y-Im_p)e64#KpmP950b;{_u z!su~RKNh*`J{)~UQLX%>^%pmrV#}y49^HnM%WYJV<8JBuD0I~U^GaX1wFv~OeCCvk z1a1m_vy1`)=DOZ?ZaWCrHi2E~9S|t)wV~B_KMdHA!2&i149FkYTwnpf7{KNMn-44y zSP-yaV1os0Gq{I1(p8|fg~Y8SZX+>+M2y5t5`Q9bJBd3;{F%g^B<>^1Gky59~xc)E3>CIOp>$(KkXv}p|Pp})KV z*e&04C^1r-R@I2jT>>z*$-Hx0m^N*O<0d5%88cpJM<5dVy4*`3fyS15_gbOV^7{XCgrs*LHc(_W(6x>rZx z9ogoy@z@!vW$o#Zm#)5Y=&!5~@4ArRiEF(adRo6JM_ms6ec{^gOQfnSDAV@srze=3 z4Zl2Lvm#20B-3Ln8sO8jz)1$_?Y{ z9A5=w1xEq>))4#dG2?s`N+l7w`Zhz0TMlw{ML&r}1I=*nsUG#p#=;Wae9C|tfpcjP zjbuFx#X}zy;PX5l9a?>qk4%swW8RtZ@$BW?jk<3OaB9gicaL8S@D_`*?jo}UB&A?H zB~u!Y$Cc>TF^UUt%1_ImTOh~riTe$$>3=uGb8cKVJkB%1gKOyHJKT>eHNt(OqC7fd z+b%uhV3j(!Eq%uil~Lyvs<*SX8L_zX%)2f&A91X+v+@Z%i9`SK{OVxR2-+rJjLXGH;lg=B^YPSkE8@XF3!lzp`OwFVtM69@P`>Cf^ zg~uJelHj}`1!M#nKtIbO-u_8-kODG-<lX94OOvuLuiAdwg6B2dAn{~B_EG7rCxRpsOK_!X%YP{o3C@3EuT4Ds|gmqH{!6@Ix6*E@TuLbvfR zD}OW6Koz>d>Cy3Z%}8jv!{#KQ3fFHv+1o8osD0bnWeVr->JZc4jwV}f%T|{||F(a@ z=!#>5tK#i~4DH-z1BG2=6AksXpq(1!cGxFd6IEQ<(r)$YmmtMJ47Jm;-PQwDKt^x` z(EsA#k(&u0s$$SwL3~vAJm5m3)G(PV+QfTvwR)=pzT0Z9is%D365Ng!;LA4|)aaZ^ z#vk7;znR&Tgh+-f_S57A`1P&1SG&v#@qad7xiIzc9(k5Y+NNyJt<)0lQg@L{AzpgS z&~3QOaXigbxIw%=OCH(z&yAa2; z(%A|RyRBbUhFGb6mxDG>?-h3-oVhW3(#_5xupIguC}}Sb2E80KS3KE~k{WfyN^xuW zugI+MnaEvyU&u+9*{XUU%@uqD<|QNjV8uWT6MafYoXZSSKt{b8w%&-_6_WgExuT~c zi1Yqbap#a%+55x2R*+B66?>N2$GZbn%x`L^JDo{J*8K4YdcG+M_p&LioeflR5v3hCrHXM*d+nxaqKenM{-OKhu-nM>^R(d4$5`R3YpTk(SeU7O$NYwti{%~2wXmOe zer0}6T@L-0vU|rTem%G<#5JuWGmco%Qm3bMuQBw;nGkJP+HN1_iVby*o?V6-}`(Y;mi>AB)3|N4eV@Lo}V8co;esv6sDMqZWtvzWt zgNhtG6%cbV))kL_&J-~$#YoM8qb(_QF^DZk6jY%<53h{-;nAqvKqS>OZ|%jd7$oGv z7UQOIaY*CtJ0ogF`Qnp2Y;2~y_e1QsVJi;&7>6h+liwhRAMDuPhMJO;Kk^M?erfO6 zuucr8pzphCwuA?e9xpJ9(QP(>rM5_uDBN*jowJT)Qb`FL|0WCTY5{l!(4 z4ZkZ)lz*FO`VI0<16BA2FK^VFWBZ?~B70bN8&d z(K)FSF@7hlG!s{Apo$ZpU6_AvQH7cjWB?h*j1jM$`E&Vz44a?ncuhP~sO$;+^Q9+} zf{lk5RJpNfAeOL15~k9hgS2MzHiy)FBQpRyD$-bZwViv)Gt4AhHr$`d5^$JdbU#?~ zvhMo+FMb~-8!zTbP-P+UuNsrXV-x$4EU+~xhPS>t;9%h(Z<87DnV@m?)yMXS>2AYh zmSATP6@g7plqnLbb~DvA-(~Ph$DevGvyK#@sDK$E6)0N=dlPYBFBbGeHbxVDU#8f- Vv4Qw!o}<3D&Sri4j(n-q@n87URvQ2S diff --git a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test.log b/Fe3GeTe2_benchmark_on_15k_300eset_orb_test.log deleted file mode 100644 index 0759261..0000000 --- a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test.log +++ /dev/null @@ -1,281 +0,0 @@ -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset_orb_test.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Setup done. Elapsed time: 1.162749041 s -================================================================================================================================================================ -Hamiltonian and exchange field rotated. Elapsed time: 1.565606958 s -================================================================================================================================================================ -Site and pair dictionaries created. Elapsed time: 1.599692791 s -================================================================================================================================================================ -k set created. Elapsed time: 1.621649791 s -================================================================================================================================================================ -Rotations done perpendicular to quantization axis. Elapsed time: 1.848495125 s -================================================================================================================================================================ -Starting matrix inversions -Total number of k points: 225 -Number of energy samples per k point: 300 -Total number of directions: 3 -Total number of matrix inversions: 202500 -The shape of the Hamiltonian and the Greens function is 84x84=7056 -Memory taken by a single Hamiltonian is: 0.015625 KB -Expected memory usage per matrix inversion: 0.5 KB -Expected memory usage per k point for parallel inversion: 450.0 KB -Expected memory usage on root node: 98.876953125 MB -================================================================================================================================================================ -Calculated Greens functions. Elapsed time: 341.658981125 s -================================================================================================================================================================ -Magnetic entities integrated. -Pairs integrated. -Magnetic parameters calculated. -##################################################################### GROGU OUTPUT ############################################################################# -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset_orb_test.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Atomic information: ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[atom index]Element(orbitals) x [Ang] y [Ang] z [Ang] Sx Sy Sz Q Lx Ly Lz Jx Jy Jz ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(1) -7.339158738013707e-06 4.149278510690423e-06 11.657585837928032 - -[4]Fe(2) -7.326987662162937e-06 4.158274523275774e-06 8.912422537596708 - -[5]Fe(2) 1.8954667088117545 1.0943913231921656 10.285002698393109 - -================================================================================================================================================================ -Exchange [meV] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Magnetic entity1 Magnetic entity2 [i j k] d [Ang] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(1) [4]Fe(2) [0 0 0] d [Ang] 2.745163300331324 -Isotropic: 0.34783476777916644 -DMI: [5.53740563e-04 1.42358020e-07 2.27549752e-08] -Symmetric-anisotropy: [ 3.04085587e-04 6.89183595e-08 -2.49215900e-07 6.89183595e-08 - -4.93936735e-04 5.24663366e-05 -2.49215900e-07 5.24663366e-05 - 1.89851148e-04] -J: [ 3.48138853e-01 6.89183595e-08 -2.49215900e-07 6.89183595e-08 - 3.47340831e-01 5.24663366e-05 -2.49215900e-07 5.24663366e-05 - 3.48024619e-01] -Energies for debugging: -array([[ 3.46973270e-04, 5.01274226e-07, -6.06206899e-07, - 3.50103963e-04], - [ 3.49075968e-04, 1.06857880e-10, 3.91573921e-10, - 3.51700101e-04], - [ 3.44577700e-04, -4.61633843e-11, -9.16733348e-11, - 3.44577606e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.0003517 , 0.00034458, 0.00034697]) -Test J_xx = E(y,z) = E(z,y) -0.000351700100715513 0.00034457760601775604 - -[3]Fe(1) [5]Fe(2) [0 0 0] d [Ang] 2.5835033632437767 -Isotropic: 0.2133083205964182 -DMI: [ 9.77395336e-03 -1.71505848e-02 -2.50531184e-06] -Symmetric-anisotropy: [ 0.00063826 -0.00205346 0.00155526 -0.00205346 0.00126536 0.00091559 - 0.00155526 0.00091559 -0.00190362] -J: [ 0.21394658 -0.00205346 0.00155526 -0.00205346 0.21457368 0.00091559 - 0.00155526 0.00091559 0.2114047 ] -Energies for debugging: -array([[ 2.09552786e-04, 8.85836233e-06, -1.06895444e-05, - 2.12830226e-04], - [ 2.13256611e-04, 1.55953293e-05, -1.87058404e-05, - 2.13948891e-04], - [ 2.16317132e-04, 2.05095493e-06, 2.05596555e-06, - 2.13944277e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00021395, 0.00021632, 0.00020955]) -Test J_xx = E(y,z) = E(z,y) -0.0002139488907894629 0.0002139442771288117 - -[4]Fe(2) [5]Fe(2) [0 0 0] d [Ang] 2.583501767937866 -Isotropic: -64.1990652847126 -DMI: [-3.36554719e+00 5.81679250e+00 5.98938514e-04] -Symmetric-anisotropy: [ 0.9024007 0.05359252 0.07068773 0.05359252 0.51528204 0.03236462 - 0.07068773 0.03236462 -1.41768273] -J: [-6.32966646e+01 5.35925196e-02 7.06877343e-02 5.35925196e-02 - -6.36837832e+01 3.23646213e-02 7.06877343e-02 3.23646213e-02 - -6.56167480e+01] -Energies for debugging: -array([[-6.58989484e-02, -3.39791181e-03, 3.33318257e-03, - -6.63392511e-02], - [-6.53345476e-02, -5.88748024e-03, 5.74610477e-03, - -6.56271418e-02], - [-6.10283154e-02, -5.29935811e-05, -5.41914581e-05, - -6.09661874e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06562714, -0.06102832, -0.06589895]) -Test J_xx = E(y,z) = E(z,y) --0.06562714176880605 -0.06096618740895869 - -[3]Fe(1) [5]Fe(2) [-1 -1 0] d [Ang] 2.5834973202859075 -Isotropic: 0.21288578271842049 -DMI: [-2.00783982e-02 1.01211533e-06 1.05202503e-06] -Symmetric-anisotropy: [ 2.14788484e-03 6.88789785e-07 -6.77046017e-07 6.88789785e-07 - -3.60175342e-04 -1.77759079e-03 -6.77046017e-07 -1.77759079e-03 - -1.78770950e-03] -J: [ 2.15033668e-01 6.88789785e-07 -6.77046017e-07 6.88789785e-07 - 2.12525607e-01 -1.77759079e-03 -6.77046017e-07 -1.77759079e-03 - 2.11098073e-01] -Energies for debugging: -array([[ 2.12888094e-04, -1.83008074e-05, 2.18559890e-05, - 2.12284870e-04], - [ 2.09308053e-04, -3.35069315e-10, 1.68916135e-09, - 2.12557778e-04], - [ 2.12766345e-04, 3.63235247e-10, -1.74081482e-09, - 2.17509557e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00021256, 0.00021277, 0.00021289]) -Test J_xx = E(y,z) = E(z,y) -0.0002125577781844703 0.00021750955693244014 - -[4]Fe(2) [5]Fe(2) [-1 -1 0] d [Ang] 2.583495745338251 -Isotropic: -64.40900602860754 -DMI: [ 6.66852978e+00 -7.33794176e-04 -5.22013428e-04] -Symmetric-anisotropy: [ 2.81831507e-01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - 1.23919957e+00 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -1.52103108e+00] -J: [-6.41271745e+01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - -6.31698065e+01 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -6.59300371e+01] -Energies for debugging: -array([[-6.51430169e-02, 6.76649535e-03, -6.57056420e-03, - -6.54030961e-02], - [-6.67170573e-02, 6.60963771e-07, -8.06624581e-07, - -6.71939352e-02], - [-6.09365168e-02, -4.22426970e-07, 6.21599885e-07, - -6.10604139e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06719394, -0.06093652, -0.06514302]) -Test J_xx = E(y,z) = E(z,y) --0.06719393518677816 -0.06106041385741498 - -[3]Fe(1) [5]Fe(2) [-1 0 0] d [Ang] 2.583541444641373 -Isotropic: 0.21320567745060778 -DMI: [9.89526228e-03 1.71523371e-02 5.08375662e-07] -Symmetric-anisotropy: [ 0.00072503 0.0020527 -0.00155485 0.0020527 0.0009888 0.00102463 - -0.00155485 0.00102463 -0.00171383] -J: [ 0.21393071 0.0020527 -0.00155485 0.0020527 0.21419448 0.00102463 - -0.00155485 0.00102463 0.21149184] -Energies for debugging: -array([[ 2.09743450e-04, 8.87062962e-06, -1.09198949e-05, - 2.12089719e-04], - [ 2.13240238e-04, -1.55974916e-05, 1.87071826e-05, - 2.13932607e-04], - [ 2.16299244e-04, -2.05219386e-06, -2.05321061e-06, - 2.13928806e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00021393, 0.0002163 , 0.00020974]) -Test J_xx = E(y,z) = E(z,y) -0.00021393260657957431 0.0002139288059246158 - -[4]Fe(2) [5]Fe(2) [-1 0 0] d [Ang] 2.5835398672184064 -Isotropic: -64.18749916411404 -DMI: [-3.34143011e+00 -5.81613454e+00 -8.04864431e-05] -Symmetric-anisotropy: [ 0.90467906 -0.05344778 -0.07078134 -0.05344778 0.5154147 0.03662271 - -0.07078134 0.03662271 -1.42009377] -J: [-6.32828201e+01 -5.34477814e-02 -7.07813383e-02 -5.34477814e-02 - -6.36720845e+01 3.66227144e-02 -7.07813383e-02 3.66227144e-02 - -6.56075929e+01] -Energies for debugging: -array([[-6.58936437e-02, -3.37805282e-03, 3.30480739e-03, - -6.63310120e-02], - [-6.53215421e-02, 5.88691588e-03, -5.74535321e-03, - -6.56142142e-02], - [-6.10131569e-02, 5.33672949e-05, 5.35282678e-05, - -6.09514260e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06561421, -0.06101316, -0.06589364]) -Test J_xx = E(y,z) = E(z,y) --0.06561421416763057 -0.06095142603633385 - -[4]Fe(2) [5]Fe(2) [-2 0 0] d [Ang] 5.951322298958084 -Isotropic: 4.59955059881751 -DMI: [ 0.3863217 0.8953163 -0.65738037] -Symmetric-anisotropy: [-0.02046716 -0.0444248 -0.03977157 -0.0444248 0.43165605 -0.07442493 - -0.03977157 -0.07442493 -0.41118888] -J: [ 4.57908344 -0.0444248 -0.03977157 -0.0444248 5.03120664 -0.07442493 - -0.03977157 -0.07442493 4.18836172] -Energies for debugging: -array([[ 0.00472963, 0.00046075, -0.0003119 , 0.00466184], - [ 0.00364709, -0.00085554, 0.00093509, 0.0036577 ], - [ 0.00540057, -0.00061296, 0.00070181, 0.00550046]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.0036577 , 0.00540057, 0.00472963]) -Test J_xx = E(y,z) = E(z,y) -0.003657703314451758 0.005500463555586455 - -[4]Fe(2) [5]Fe(2) [-3 0 0] d [Ang] 9.638732176310562 -Isotropic: -0.23103792475126567 -DMI: [ 0.09123182 0.33044193 -0.36639877] -Symmetric-anisotropy: [-0.08681631 0.00193669 -0.00538302 0.00193669 -0.07518869 -0.03449947 - -0.00538302 -0.03449947 0.162005 ] -J: [-0.31785424 0.00193669 -0.00538302 0.00193669 -0.30622661 -0.03449947 - -0.00538302 -0.03449947 -0.06903292] -Energies for debugging: -array([[ 9.75938878e-06, 1.25731290e-04, -5.67323532e-05, - -1.49392424e-04], - [-1.47825232e-04, -3.25058913e-04, 3.35824948e-04, - -2.23552488e-04], - [-4.63060805e-04, -3.68335462e-04, 3.64462081e-04, - -4.12155989e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-2.23552488e-04, -4.63060805e-04, 9.75938878e-06]) -Test J_xx = E(y,z) = E(z,y) --0.00022355248815234167 -0.00041215598873809247 - -================================================================================================================================================================ -Runtime information: -Total runtime: 340.88071287500003 s ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Initial setup: 0.11835845800000011 s -Hamiltonian conversion and XC field extraction: 0.403 s -Pair and site datastructure creatrions: 0.034 s -k set cration and distribution: 0.022 s -Rotating XC potential: 0.227 s -Greens function inversion: 339.810 s -Calculate energies and magnetic components: 0.266 s diff --git a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test.pickle b/Fe3GeTe2_benchmark_on_15k_300eset_orb_test.pickle deleted file mode 100644 index caec7936f2781c39778cc84a8d97d2741b32013f..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 7519 zcmcgx30MDzRLxb$*N#P|sE?dM$1To@7DO+OaC5{sb#B8qA zP9%=vN7;EJjy}j-#KBG|7KsJy1Vmyd6?5zaY?&j=)oq%(+%|$6fuWy>>TBIojJJQbfiQOTt+)L@Kk5;d7z|2!;-qaw1sq(8LlxB9gIX ze6a{aM~LLYm;_snSc2FJEs?MjFw}r0K{#TOR3?#gWLPGK>L_o7G(=q0Etwco zBNg+cg=61f!I(9Z#?(s1c!N_A(_jwFVL>MeF$>Wx@Wn}o4Zt@wmQ zc&+!5L=hYP$nEu~9ulgh(QZDeI*LrY2lb*oM8} z^mB6?yGY==Bf@IK)JU*>#@T`cPLbgHH2H#rX_3HkUv2$~q6kpr?Du%okqDq%7h*pz zXKM2#OiiBj!`AiTR3xtVw4M%!=mu>Dv@f8!Lz@YW4h=wKK$``PXp^_s#5E);`KKyZ z&xu$)2e^8?w2!Ua1rM$sv-qCki21~d#+XWa%k-|nldg?ZtsF1o!>Q0xFPMA!yiJWw zhKy^jHm3;uM@Yvq?N&pE)4K7q+nyRSHjAvT)YKX>c;|F)-h$X|!F2jtFdsM|Uud(T z`9Whsn*+@sS^%^_XhG2WTdfkp<2n@HSDVkU_giCakgj>N4bW|6p!#B373 zCviK8J4oC~Vh)M9B<>=y+k!a{Y{9$@gC1TF=>k*j9wg+Zw1L!;{r8(TSq@~uEHT9r zSp<#2Oe5nlM>7lw(4qL#Mz;uNd<*8R)X+72S zxefzc2U8jz6ti`=SOiZgVI#e(Q(F~I!n#J$YUm{{AmEnOkBybN^py#_hWRYWQv~{4 z81`G63wRlvdN;M!1r#`SiDj(2Gvn*dPs)4rdgdj)86woc6^yj=X*b;B0`lSyI5vfy zFR4wi?R@G`Q?lXC3Y54`Ob)8!$Kp$226UGoK6fSURI^7FP!dPA{8%Bq%*-G-Gcelu zX6^LY`_mIg@cc5!)=QSqHI@vw`1*K;g7Q1{pA1=U$G~%JCT?q7Js42slpoe*+ZOu3 zT1A8A)k{xACP>8QCjTEUj+x*AytCps;iK2k-R6wcyD18#D+fUH>#e5pDNC^I^pCcC zfohq1LM99jT?} z+P^83^*Vd0r&fiZ);8Vj@NzXHiUGxVnAWVMiE-t9GW_n)`U!zU;8eU*1?|6~3g-6A zSDHOlL4^w#IR>bTjT=JhOYarn>$h#yYP}JS#~RUpwGW5}JkI&H?>h?c=t!T+)5i;N z@#`_m6>qI4yGHY&9)X`{QR-mfEkZ z1x6)@lT0_~s>><=s&B!_>pwp4R~4OeGiiCv`lVKpjSR%ggmH9)^8t@sPvu(iKB|f{ z8Es~-eh4h>iy>`&uJgK^fu$59Y6O%&ao?ee_((Vv?^MxmtRP-wyC2=3FI22(mFDDW z^i&0v+YVe6-t(}E`3*LpqVcasTkP8v!YOs1H9vk%*RT6(;;3sjz;0Gel8e{|c@@Jsv~#sTiUNGYEmw(5~X8!nYA!eqoaUsWQLUS_Y^9hAs<)8rluN8w{_)xGAwb; z(bRV;2Na{84C7Lkc6#Xiy;#v*5yW|Ks;JrTUitPwk0;3Aj}<$YTF1E@F2Ma8+dS=l zNdhU49w8?i6LI%(SMJSH6yT>2&!lC>h4|dd^rO3{?xDX-XSm!NdQkCmo%!S($3mQ5 zyvy0hF(1#|=NUV5;dX_0xiCBVb{@TAgv07U&-Y`+Si2mnsb&L+6|cX)qxH?8tKh;( zM$qS@&G6OLRo601hO2$sgY_+^uC;@6q1&>*G(N5_r~H=6+xeFBU-qjCY4y`#nTO06 zY15LsQU?3tjh~&J^TKSFs@_LqMN07Zx#RVMO8a71H~p;Z$t{7U6eDT`lpnEk*@@Jo ze=k;0YM@a0yd(1p>(ljTO#bs&(bDv3>eZDtz*O#N7@qxrp7f-&W%I*_ihvyxD$xoX zAT3;7WCinzWp|6@O^aF;8$;ffAKcMS7k~N2lmqhy_uV^|R4kkVnu^Y^Hgqjh94)&N z@9tKn@G0?-di2OEiX%E}g8c>%E0VX~TGdhJ1OAL!YH|Ha7+A2_dH1S>g*fG>{P5y% zsACZD+&Q#u=lLLYIpyPRT;Y7sCv&Gu(M8byP#!7jfHgLy?~6Axqelz2LtYt``*ik8i4%2S`&4h-;Wj139CBa4jw?Ph!t6_ z9@pLs+_PHRRxdq>pB*&g;_1+x*JgClrO7{oO$$E|9ORd zE&C@>tY7N%#pE-`JWWb_7)U z9`zHRigCgPf3YM?Qz@w4d$=XJRtUy(C4txV=HNGC*ZwggF94*u?OAi`xex@GWf(S& zjRg%mYKByg@WLldA2-hCtv9gZg)ZOwO)OASCbP{0_p`RPB1R0Ii|FSS

iS%Q2yeo8}+}J4pshb zLcW>ytBN-RP=!}eSi>hijHZQ+$EyP-M{MzT#yj*JZJae-fG6h5nCd?jUs9Z`bw8Aj|5myDbl&*s`0Y{b zWBVSs0xy^0tFCK)p_*V)6%o1dTlQ?ei1+&r@_NE>gDkh7agU3WB8L^Tz?ow~7n2QQ z)&9dNw;yDHBlwLSVM8D9--T1rYo7a(*>I_2E|xc?@;+&+MHgdYAmrPkmWz5tkL}=}ncqGC2Nw zfrdS5A?3>uWg+pe8jH*45Z?gV@N=M4eQVZ&xXbrPp!(4)*tq)YpFW0TGeS^P_%#lZ zz%PynOCnZ%K~&d#lPS(AH_$}QBP9qTVnxVB%GTkxf>`(=6!t?l#*i=i*6zZ&2I8N2 QuKL=ti*CMILb=@bU%E*zfB*mh diff --git a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.log b/Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.log deleted file mode 100644 index fb7e858..0000000 --- a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.log +++ /dev/null @@ -1,283 +0,0 @@ -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Setup done. Elapsed time: 1.35647775 s -================================================================================================================================================================ -Hamiltonian and exchange field rotated. Elapsed time: 1.7424465 s -================================================================================================================================================================ -Site and pair dictionaries created. Elapsed time: 1.774876708 s -================================================================================================================================================================ -k set created. Elapsed time: 1.79503 s -================================================================================================================================================================ -Rotations done perpendicular to quantization axis. Elapsed time: 2.036713416 s -================================================================================================================================================================ -Starting matrix inversions -Total number of k points: 225 -Number of energy samples per k point: 300 -Total number of directions: 3 -Total number of matrix inversions: 202500 -The shape of the Hamiltonian and the Greens function is 84x84=7056 -Memory taken by a single Hamiltonian is: 0.015625 KB -Expected memory usage per matrix inversion: 0.5 KB -Expected memory usage per k point for parallel inversion: 450.0 KB -Expected memory usage on root node: 98.876953125 MB -================================================================================================================================================================ -Calculated Greens functions. Elapsed time: 377.287089208 s -================================================================================================================================================================ -Magnetic entities integrated. -Pairs integrated. -Magnetic parameters calculated. -##################################################################### GROGU OUTPUT ############################################################################# -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Atomic information: ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[atom index]Element(orbitals) x [Ang] y [Ang] z [Ang] Sx Sy Sz Q Lx Ly Lz Jx Jy Jz ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(all) -7.339158738013707e-06 4.149278510690423e-06 11.657585837928032 - -[4]Fe(2) -7.326987662162937e-06 4.158274523275774e-06 8.912422537596708 - -[5]Fe(2) 1.8954667088117545 1.0943913231921656 10.285002698393109 - -================================================================================================================================================================ -Exchange [meV] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Magnetic entity1 Magnetic entity2 [i j k] d [Ang] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(all) [4]Fe(2) [0 0 0] d [Ang] 2.745163300331324 -Isotropic: -40.87189112563518 -DMI: [ 6.16487380e-03 3.16430556e-06 -1.59991820e-06] -Symmetric-anisotropy: [-2.32043212e+00 2.53801348e-06 -1.08593736e-06 2.53801348e-06 - -1.35046068e+00 7.74258719e-05 -1.08593736e-06 7.74258719e-05 - 3.67089280e+00] -J: [-4.31923232e+01 2.53801348e-06 -1.08593736e-06 2.53801348e-06 - -4.22223518e+01 7.74258719e-05 -1.08593736e-06 7.74258719e-05 - -3.72009983e+01] -Energies for debugging: -array([[-3.62152318e-02, 6.08744793e-06, -6.24229967e-06, - -3.59883911e-02], - [-3.81867648e-02, -2.07836820e-09, 4.25024293e-09, - -3.79283356e-02], - [-4.84563126e-02, -4.13793168e-09, -9.38095273e-10, - -4.84563108e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.03792834, -0.04845631, -0.03621523]) -Test J_xx = E(y,z) = E(z,y) --0.03792833564278811 -0.0484563108453157 - -[3]Fe(all) [5]Fe(2) [0 0 0] d [Ang] 2.5835033632437767 -Isotropic: -62.028936345638364 -DMI: [ 3.23051343e+00 -5.60964351e+00 5.89430234e-04] -Symmetric-anisotropy: [ 0.90811936 0.04789656 -0.07353011 0.04789656 0.52317596 -0.03533139 - -0.07353011 -0.03533139 -1.43129532] -J: [-6.11208170e+01 4.78965641e-02 -7.35301053e-02 4.78965641e-02 - -6.15057604e+01 -3.53313894e-02 -7.35301053e-02 -3.53313894e-02 - -6.34602317e+01] -Energies for debugging: -array([[-6.37479933e-02, 3.26584482e-03, -3.19518204e-03, - -6.41681572e-02], - [-6.31724700e-02, 5.68317362e-03, -5.53611341e-03, - -6.34538327e-02], - [-5.88433636e-02, -4.73071338e-05, -4.84859943e-05, - -5.87878012e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06345383, -0.05884336, -0.06374799]) -Test J_xx = E(y,z) = E(z,y) --0.06345383273317398 -0.058787801246131074 - -[4]Fe(2) [5]Fe(2) [0 0 0] d [Ang] 2.583501767937866 -Isotropic: -64.1990652847126 -DMI: [-3.36554719e+00 5.81679250e+00 5.98938514e-04] -Symmetric-anisotropy: [ 0.9024007 0.05359252 0.07068773 0.05359252 0.51528204 0.03236462 - 0.07068773 0.03236462 -1.41768273] -J: [-6.32966646e+01 5.35925196e-02 7.06877343e-02 5.35925196e-02 - -6.36837832e+01 3.23646213e-02 7.06877343e-02 3.23646213e-02 - -6.56167480e+01] -Energies for debugging: -array([[-6.58989484e-02, -3.39791181e-03, 3.33318257e-03, - -6.63392511e-02], - [-6.53345476e-02, -5.88748024e-03, 5.74610477e-03, - -6.56271418e-02], - [-6.10283154e-02, -5.29935811e-05, -5.41914581e-05, - -6.09661874e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06562714, -0.06102832, -0.06589895]) -Test J_xx = E(y,z) = E(z,y) --0.06562714176880605 -0.06096618740895869 - -[3]Fe(all) [5]Fe(2) [-1 -1 0] d [Ang] 2.5834973202859075 -Isotropic: -62.23514735929866 -DMI: [-6.44473800e+00 3.69898564e-04 -5.24592826e-04] -Symmetric-anisotropy: [ 2.91447198e-01 -1.04541004e-04 1.02420326e-04 -1.04541004e-04 - 1.23896201e+00 9.61259121e-02 1.02420326e-04 9.61259121e-02 - -1.53040921e+00] -J: [-6.19437002e+01 -1.04541004e-04 1.02420326e-04 -1.04541004e-04 - -6.09961853e+01 9.61259121e-02 1.02420326e-04 9.61259121e-02 - -6.37655566e+01] -Energies for debugging: -array([[-6.29770048e-02, -6.54086391e-03, 6.34861208e-03, - -6.32308800e-02], - [-6.45541084e-02, -4.72318891e-07, 2.67478238e-07, - -6.50151564e-02], - [-5.87614907e-02, -4.20051821e-07, 6.29133830e-07, - -5.88722439e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06501516, -0.05876149, -0.062977 ]) -Test J_xx = E(y,z) = E(z,y) --0.06501515642010514 -0.0588722439032051 - -[4]Fe(2) [5]Fe(2) [-1 -1 0] d [Ang] 2.583495745338251 -Isotropic: -64.40900602860754 -DMI: [ 6.66852978e+00 -7.33794176e-04 -5.22013428e-04] -Symmetric-anisotropy: [ 2.81831507e-01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - 1.23919957e+00 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -1.52103108e+00] -J: [-6.41271745e+01 -9.95864576e-05 7.28304048e-05 -9.95864576e-05 - -6.31698065e+01 -9.79655745e-02 7.28304048e-05 -9.79655745e-02 - -6.59300371e+01] -Energies for debugging: -array([[-6.51430169e-02, 6.76649535e-03, -6.57056420e-03, - -6.54030961e-02], - [-6.67170573e-02, 6.60963771e-07, -8.06624581e-07, - -6.71939352e-02], - [-6.09365168e-02, -4.22426970e-07, 6.21599885e-07, - -6.10604139e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06719394, -0.06093652, -0.06514302]) -Test J_xx = E(y,z) = E(z,y) --0.06719393518677816 -0.06106041385741498 - -[3]Fe(all) [5]Fe(2) [-1 0 0] d [Ang] 2.583541444641373 -Isotropic: -62.01268695598744 -DMI: [ 3.25190510e+00 5.60946137e+00 -6.89689002e-05] -Symmetric-anisotropy: [ 0.90589405 -0.04774783 0.07347635 -0.04774783 0.523936 -0.0317999 - 0.07347635 -0.0317999 -1.42983004] -J: [-6.11067929e+01 -4.77478280e-02 7.34763489e-02 -4.77478280e-02 - -6.14887510e+01 -3.17998960e-02 7.34763489e-02 -3.17998960e-02 - -6.34425170e+01] -Energies for debugging: -array([[-6.37256875e-02, 3.28370499e-03, -3.22010520e-03, - -6.41494027e-02], - [-6.31593465e-02, -5.68293772e-03, 5.53598502e-03, - -6.34406363e-02], - [-5.88280992e-02, 4.76788591e-05, 4.78167969e-05, - -5.87729495e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06344064, -0.0588281 , -0.06372569]) -Test J_xx = E(y,z) = E(z,y) --0.06344063629910189 -0.058772949517250955 - -[4]Fe(2) [5]Fe(2) [-1 0 0] d [Ang] 2.5835398672184064 -Isotropic: -64.18749916411404 -DMI: [-3.34143011e+00 -5.81613454e+00 -8.04864431e-05] -Symmetric-anisotropy: [ 0.90467906 -0.05344778 -0.07078134 -0.05344778 0.5154147 0.03662271 - -0.07078134 0.03662271 -1.42009377] -J: [-6.32828201e+01 -5.34477814e-02 -7.07813383e-02 -5.34477814e-02 - -6.36720845e+01 3.66227144e-02 -7.07813383e-02 3.66227144e-02 - -6.56075929e+01] -Energies for debugging: -array([[-6.58936437e-02, -3.37805282e-03, 3.30480739e-03, - -6.63310120e-02], - [-6.53215421e-02, 5.88691588e-03, -5.74535321e-03, - -6.56142142e-02], - [-6.10131569e-02, 5.33672949e-05, 5.35282678e-05, - -6.09514260e-02]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-0.06561421, -0.06101316, -0.06589364]) -Test J_xx = E(y,z) = E(z,y) --0.06561421416763057 -0.06095142603633385 - -[4]Fe(2) [5]Fe(2) [-2 0 0] d [Ang] 5.951322298958084 -Isotropic: 4.59955059881751 -DMI: [ 0.3863217 0.8953163 -0.65738037] -Symmetric-anisotropy: [-0.02046716 -0.0444248 -0.03977157 -0.0444248 0.43165605 -0.07442493 - -0.03977157 -0.07442493 -0.41118888] -J: [ 4.57908344 -0.0444248 -0.03977157 -0.0444248 5.03120664 -0.07442493 - -0.03977157 -0.07442493 4.18836172] -Energies for debugging: -array([[ 0.00472963, 0.00046075, -0.0003119 , 0.00466184], - [ 0.00364709, -0.00085554, 0.00093509, 0.0036577 ], - [ 0.00540057, -0.00061296, 0.00070181, 0.00550046]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.0036577 , 0.00540057, 0.00472963]) -Test J_xx = E(y,z) = E(z,y) -0.003657703314451758 0.005500463555586455 - -[4]Fe(2) [5]Fe(2) [-3 0 0] d [Ang] 9.638732176310562 -Isotropic: -0.23103792475126567 -DMI: [ 0.09123182 0.33044193 -0.36639877] -Symmetric-anisotropy: [-0.08681631 0.00193669 -0.00538302 0.00193669 -0.07518869 -0.03449947 - -0.00538302 -0.03449947 0.162005 ] -J: [-0.31785424 0.00193669 -0.00538302 0.00193669 -0.30622661 -0.03449947 - -0.00538302 -0.03449947 -0.06903292] -Energies for debugging: -array([[ 9.75938878e-06, 1.25731290e-04, -5.67323532e-05, - -1.49392424e-04], - [-1.47825232e-04, -3.25058913e-04, 3.35824948e-04, - -2.23552488e-04], - [-4.63060805e-04, -3.68335462e-04, 3.64462081e-04, - -4.12155989e-04]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-2.23552488e-04, -4.63060805e-04, 9.75938878e-06]) -Test J_xx = E(y,z) = E(z,y) --0.00022355248815234167 -0.00041215598873809247 - -================================================================================================================================================================ -Runtime information: -Total runtime: 376.47007308300005 s ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Initial setup: 0.10374624999999993 s -Hamiltonian conversion and XC field extraction: 0.386 s -Pair and site datastructure creatrions: 0.032 s -k set cration and distribution: 0.020 s -Rotating XC potential: 0.242 s -Greens function inversion: 375.250 s -Calculate energies and magnetic components: 0.436 s diff --git a/Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.pickle b/Fe3GeTe2_benchmark_on_15k_300eset_orb_test2.pickle deleted file mode 100644 index 776b91616813a88a886b363f6f618a6987d1c218..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 7774 zcmcgx30M=!x+Me@Sw=(<+)yqmhzQ8y3h5{eB2hrlaRWn48WPB65>`=B)B&dzk#^jL zaZqu?h>W6WDx%{ydL0xHm#4U)C?XEZxO1zM?g)y&_04U0qeDs!sh?ok}%W z7iO!cxfbG%n6Z>2aYq)7MPH8dO)x{AB1wu8S@Z>?@r-nBDjlPF<3T||fz9yJ{IK*C z_`D2X8*OKP`lXS#I84SYW#|{A<#+_{NOU$$g7b~}W&{m>_m(8AVX;&(p9N2PyLhbey@(Du=Yo_VX-M{r?Yasj4UnnkxS<*-{lz z6>Iwz#SCJh>~U6i%rKTGfU^2fIazQeRvuKzB7ZnQI5*CZT0_R&|`cW-Z&9`Yw(*ZyBy8iFK!@gruBq>~`dHxZk zRQe72!`E1s?hxbl{SIZNw;9OvlkB|j9a?O={ArdalW1HV)Z_f^| z4?lNX<-<4r$c!#y@I(RD2SOY)is4{Js2G(+6GCaD1eQM3F)YJvar#R)MF zn4iSl>}bGN!*_@)Qww@)NX#N}Es5DA;v}vk@n;g(lemGzjU;X&aWjcqNZd-|HWG74 z%q8&`5?c(ISLX)I-z=o2JnSPf-1S~k?$Q^?@)Nu7Jk1)^nE@MwtHKpgG$xC`d+Ixg zA646SCq>hkm5=aA;SkenEnlsjNn^UZtkR%nz);Lk$|068xaPyzhUs$zaH|C&hFG5H zp>f={X0>*{*pRLK0*I0^1D-&x!1+d47de+BW^mFkt=tZr*3tf&nSPz&zUidjHkTVJzxO_d7NYGSRWj9c z)T&=Wz%Bg%&J8H&V&5#!w%iaWDQXT<&ru612(;&(Y8&>EOyihduV3ZI<7Ip5!0 zW!<;ry)S&OX@p@QTOfx!zuN?Lmi?W9nx8ik3JV8_BzXA5T~knZX?#b;sy^4N0O{v3IC~D4UHM{e9qi@SZL5Kw;29BkRm6V=E|v3t?W zUfThcPWdas4?SI1)1C1X3WmrTjgw|Evn3i=3Tz_dV^#lCH`Nbi=2r*jx4{hYr z>7s6bM16c_keXsdF`)eVhd-Y@t`1UDjHqFNp4BQMI{(I?@U}DD%j&IDw>I~A%Bm;(Pirwgw} z9%!Ws`lG2SCrm>+Rs}sIQWr@Zbk zT4dSmLI)BMw&f9+9tR0+I;Osi@#1Eb^7CVkos#3->|uQX^Vr@6r;}> zZUZk)ez^PJMT?f2Anx1J#r56ZWzBn9ts>tNEw)8CCXC##1XF5X__`ihfGmCd7(HH_ z47~0Byggo}1f{6&g1MFl!L-YaqV2p`%cHt<%C9fGd0@~F9jiiK^}I*+zCO(K5vd=SngoE_ zfoCWG-KR6ESQ@%H*WM^Z{WXU7N1q*ad|j}bVnnro@}suRJ+^$ow?zv|4HUMF7LKPL zoU#5VqQ(7(6**xcaMo`qJSeJ(LHw6h{b|aH0lgOI%(D?H!JrQ-o4tl9LDP)Pw;%iy zfI0KKJvIB{kwt{KH%G(7fVrix!j~C~T<<{W|ByUafWUEz#oVfcp=k;ojovzM{#2`(nXQAGW+I zu!u!$Lq{Iz0rwW~Yy0KY7AwIntVd&YzeL1s4sEG4Jpo)ks8?kT_ZHkqfmi$c8~_jR z@B%B}njpMqcl_s@n<58ZHY9iPY8@?NlNP;iHt9sPh!;C7v1_bDyf@UpSQ4=pob5LD z(ox$&K>4pc^56I`szz20&f|3|3Q&sC3FdXrs7vFRE3v=y1pRkz7P5 zHBi`66~=XE#SZ_rZuuqhEvo33dA2WyEl2vCOFLxZ&Oi!^0!5`#Pau5Bma%N)NY$?W z^{Le&#GWS$E;pJCD&m(tw8{%Y(#PyfD{T@Xp(odv)!M??vGscQDl0!Qc(k3}ux5Y6 zfj@Kpt{>wO4Q28>WY3)q>l;uDl5#uWA(p3d`ptZW0|T_VsY^Pc{nN}2>C!2`33qdT8gfHBy6t?2)L${!x#(yW@-+R{9P!RNWLff))f2is z?1U;FAt}ha5DVtl7^>0}x4#@<$)p%jEuj24l@$-}s?9ZjpLF&+_>x}y*b7C_eLRpBYQ3?*B__tU{V!P zxrys`uDt}bf4sD#inoJ%nr;|#7)-b{G`dF#8#z-Paw*j`Ugy_1-@La5IRq-UMw&m_ z{R^PdDPKA?Apt*rS}VzRRPmtv(wo0$A3^SCzkU2_a0O!aA+DmQ{8}eeQP}(Bl%q4s zb&Mzm6k}0edD^IrSHH@z?wP*#pnV55HN~M=S6Li9>JXutLsetBf+JI~HU328rHPxX ztSAgKhG#f&QtIieqJPmNyOd$y1+Y0{9($&*Kb7+CZg>1#(h_WdT%bTTiNwEaY@UEi zoPcuRVW^yX^*#E~Rx%*->#~P3iFbBOaau!+OZVmq9EDQ*jh1&jIY4_ZG)xYku47q2~c5-I9`fKY(2y) zMM5#cliG^Jk%CBDPt;*98iY=_jTVcqu#ctOaCR(6 zVu=71$q<=9EK*?OMDpmEI2*25g4#sOg)#vmksxsjY!q99a>XL4Od{vX6x$S-w(3Q| zL&Q_RvQ5F#NLIMhIEsym00o0ZV`(NU_#-iGo-8f~C9=>3PKvZtwO13PxhVYn{Jgut zrTSrM$P}D3cVl&D{{2fMZxz6e8A;H7q@{|b3I@^HG_iuO$JZxl@C|yC7!Qi2iut=> z^_{9%`ihlAA%o|pa`i2R$UStfzPX89s2ot{n+=_Lz$&MFkZs>ARQ>-CY^ka!3aTpq zFWFKRQ5CEE7R9upfb114157)LCjhe%qsIanYoQ#mOyT_!hz{HnZ}?87VcNu9jKcS^ zCZ%h@WiX?>_ zRL{GV3BP)LQQ@xFP9=HUYoRs`D-H}N&PL7R zc`J|!(4HBfg~-Iw3Kp1LLIsPC8Gu_ z+6p9Vkfwq(4I~?oY(cUE$sVNXAQ2ofLurga#O8=&KTx$1s2a@7%m9#yxB}In5kz7z zi6JB|ATgB0g(NN_aWRQYNL)(dG7{M&hLMPn$RUwSB9BCrL{)>dWJmvb_BY9L;;*Z8CU7OHVXFkMuHN+O8i z8=wRh9jFVGVOLzLZ+z#eK66j0m-StM-|gD`l}?*JbJk^hMq{f!6aQTKsOp|Rb3?-I zJCy)c8X;oKa+Vff!lLt~{YxiCJwDBy;RKpE3nXWdW`pDck}F6MNH9oDkmi6yl*wP+ zzZ_}^W57GawOkDv5)!2(%1D%xxPru$B*u~$M`ApQt4K^BF_Fa8B(5PbiNv)et|M_h zi7Nids?l?#MvtHRp^>}n#?xh;uE&p?PPy_&ObyiH)HTww_A~X!@wfEv3eBDYbJ-l6 z{s$PS+OMq75%^u;+htl{z#LP|zj+P@Y^%uP>P;|E{PK_LZoNNX+F-zF0}U7Um|CO*TG<#C;^@lK2ye`$_E8Ni&E}vIi3Xa3^`{`#-!M*ac6s zZHUWU_Z&_>d*JTV^e=|fNjp25jvCmyDuur1;nDG~O6KGIamwAU%GBd)bK~7yl_g0A z*Oig5viyZ#asxnuoE*ipV-O;kgVhfx9HS!wu;~C0Nd#)m06BM{nWXNHe)|ft)v9j9 zw0HukOu^U1K9X_~AtC|W6@H|M0!N zb;Ei=#mH>sWKP_I$Y(d8d%wJ~%51rT&+O8__)|-kGIzoD)A|NEO2;HO@3avTqD}tHXSL3jO&X6IvqW5 z%1(H|c~yF-L#fgEQIYN+**$*B&f?IztL5s*gV8(&sELY)Bb^IBTOKKgSGmVrZ=CxC zZ{0r0CSR)zS2cjrGNV=WqfMI z3iS#Jgm)D?%v`tLUAdZNm0rCM>9Jqw&bFHLDFWBooVBVC3WO;?Wp`$Tg}7aeR0)*d z(Z3!_7wPRBhwjI#l(zP7YQ8;~2Dg>lxNf&BgjA&84Un%#6)dNSB@>0sJ&Y&@6l0|? z8iGli2W9w4VPi+$T+kInd{L`{_P&I3h>ML_QBsDZG}X88W!!hXu>%UzJ^C$^no2V zXlX_QRdhNy$>MaLlCtrHaapfT!r(PWT*gmyUfAQOY<|vd9ggnMsyN~@sv~z|0(`gZ&CGF? zX-eD9o73)nQZ@`#OcCseQbslp$}r(%O=ojV^PmiuTx4Wi@_bJfL&gf?Alv)uy}nSf zqD`8yPou9YsQu*dRdIj4xbz*(7yjL^_2TQQHt754IdSa^AL1#`uGJ<*zI12t?6wXXqwyW$0LIs#v~yzmF)No$9zmEik{*0iyBtF6Z8 zTV?OZif6mzFIJ($P(|4O=s1%LOJS16zRr-#)qzefxiE7s=9LiK_+!OuDbkD5kbhHe_4t(n>O(=b#~I=^tuGFkJW z498e6$k|wu>J^Vf<|wFc|f0OQ+pVuYBMt=Y(KvkvUXjyj~s=^BFEbUQjajHy=2a z_i>&!P=&FPQB^}m$61%n@p+f4JJ0GJ)iv1CM-{sY63-nCABHN%uJ&&`b2sfydi$Qv|xlk&!) zNqQ4rm*OL`nTy;Xh7Lm&#&pq|sFlrwGSrv6tCa#(P>lLAJeWDU-PZ3tRSY{;P->u1 zwce5Y3dZ@yvhn|Yta!1ee(H{wKJeS5ji$`3`}l^>4X6D(^j9co_mVB{U?`tCr(QZ2;|tT?=5xn%<=_Q>Y8H6D$icU1 zVLIKB@5hSnuI{|6!&yP5E?j*ub9XS@*=Ryvoa7FN9{PQ{@jgUJ*@ro2FWqlj21i{9 z-jV6Fu*Xl?59>{x=Xegif?Rd1`}or21m%gzKbqwpPE#^})NOCn*^dvsuds%KP!(O|5<;ng z!rrRTYrZJ5s2czN8RVatLQWO4D~{wBk27cWmqua0In zK&N(?H%;HKKITzpko{LXrp?HLltZ3D7KKYHO~0M3yj-!tyYfd5(q zpU+o@LdyBe9#$&;y- zjB}zhan8`HVAC2GZ#A2%M4vKj8`oCh9apvm#CFhTK=Az$|%uJyYp-haYYNo53fWu9t+IWDY|W7o5#gFZcxw%5cY&+yC&3wA&;S+ zcF+$e!EI*PXfHM*;;}>BJt@!I7G}YzGa}4PDv+TnC-JWunr+2xt^2hetg^5*CsWl9&SlFT4*{L0FYq ziU)@xRY66t+TB={sudJb@vM}BTB{V5T2K6ElO04rNuU1zzt1yy_uFr0=G~chzWFwr zV78gCW9g5@gadA+6sQEzmJ4b{^5V(axt1&d* z-?_=fkC4a(GQI{=Yn+u5QIrfL@E&SW1V0uOg-U|SH3E%9At&%&a&5FS&RL{TVb0N7 znMNW|sRVHZ-kGn$L<+fDqtc2rL@I%s>u&`8$-@n|q!K*VD#D8`B-Ri?gaePwGhIc9 zJK^Txnm8p!S22-riFIoYUQ?XyPWb!#Pk9Rt{SS|&A%yF^`WOc7-7lLtC4d|6slr;3 zRugjx2YOK13PNlpwx(%_yEZ4$8#D`>dDpPvOy8_Y%mTWR57&mZTDC&hZogK`+H@^U z3n+^}ZNGAxUD|(?ZObaK|NjG9*cGIp@A7}i7Ip+XHjFKXn?y6Z+tKNfu_Ofd750!EzW#?aYx2BT3 z6*<4T2Nm0``eV~*5BlM>c=u>Qq#V;oM0~JrB^n9Vv;xco8bvh01B**W@J#To3YAcz z0SjLu5648{w$^I10&4|4gJkr(mD%<)I1B@%hlAn>$_P-LKye1e1r!b_BS9Gj3e6Gk zrdCSie4!%tBVG3dU4xamBM4NYk6;6^LKqBXa5{rC7@W!AEC#<~a5jT;7@W)CJO=p; zhA}8$P{^Q&!EgpK2K5hQrV&Ka!;D+bcAe{kIgA=c2yoNbxOi}_^{CR3Pf7=!M@P2K zIr9rSL2T9N5m}CeNBi|Y5I9snj{@8hlVhq#divUEfoB5L1=`>YaX5Z@H}wR)c?pfm z1Ni5<$fzMFo}gvxxlZGC0L?!^{yg#om8PbYKLV(y%jUJ!Cd10(r!4p0Y3C_c)xPYW zi)ASDoQ9dR=JDw`FHvhSxhRehTjCwmB7sbx0ufv+^|EAz)kV~9bV7;MAj<1w*nFOu zSj96Dt6M(>t*-cWUJhK)#TTHsfifNxcTgsPf`CGI?*Yn0P-yPrC9T_mQc;3?=wrSC zSSkk93~CtEG8n_)0tRCljAJmK!G#PiVsJ5oOBh_rU;=~77+lWa3I_H3SLszU8mdIR z_G1b;=NWFvJN65CXn_6s-<9P^9k;R5u!_$N3K1u_9tzVYVaynjqbLG`*w$dXLk{p( zUmMFzfMxrZ+4%!uBu=nD_7Vu9TTsSluNh&8x361~Xj1^%HK?f~uu!GN=i={1wV zEC#a~%wcdhgUv>K3_ao_!MgnTi2GOv-n|^$h&s7c$L(1D7+v*4?wvn24K`-P*M?SL zwO(8*WODK3G#;0_nO^F+J%~&FJmiol$(Kucq60nBv%a7@KVBY~15i)yBO~sSm|t;e z+a$CqWaij&Hb~nep3`4={QN=0u--7-Tq&TpGqC@I`~f!+NI=>JKuT{_ra@ZK2i85E z0x?6@3TT?RSvXkogkhGil{+}~d1mgR?^6u)?9eRifm5h$tZjjs86sC$%3d?40qP>^ z89LaSUByDBkEU2_ZnLOq-?`4VUsyyPs>v5VUr)j?L!5AB8=$2Di~PuH{5hwH~Y zmiP$>*H7Dy1nmh${L(Pgi|!R%FwHajphnSkd?A^E1fjy-s~~ zq3dq4V2<<8A?vwl%#jYgR^VLp%%bwCv+w1SXM*B)+|MW^xAm5FRK|1BtECY~XZdka zyG*1{j$bmlYUH%Y*k#G&!2Y#;Y^QPcyP?zqB#8#rY1Yxsh+oo9^z~M63{W3k zroTA*b&nwGvX5^~z6l2v-boKj{WTDUc70Cjwdvn_ScLWE8^wZzzP9p8CS1r9xon08zw9g*42v(LG7$qMW?i zz~?nPkz22(YUhvNi~N`nT=M5ia6My$#o9zJ>YzKedH=jT#4kQ@oo`AJ;&nddgkS*| z<-Yu@Fa{{Z%f@cuV6S9Sc-d*G$LeHq&$Ei#qgpFu(4XF`_Wx>}LJl{aP&_U8lFH~h z&VM>k$fvuuZLk^{)MOub>_7akk}1@tVFtsXXjrJSn)1bqL5chp2MlhPfxPmQ%Jip zgnqAYe!yLamNbkZ_3G_etd9z3+vAr0STfBse4&;ApRj}d+66SYF)KOK4=&qdB zB`CP)-mqOZxtUxvKcVX8gbCfq2=eR_X)PN$P**CX)^X7gkC5^E1Gwme<4c7zqI!@W ztCsfNf20@UylBgfu|BP0$d=3BPIap?P9bN?(|kXhH;EG36lga;oKIe?^UNHw5J#Ym zm-SvXyDEq}xALLI>zGgq`k}2HSK}Z{nrnGuE6={2r%nOP8(t91& zZ!+!15GTu*;}fq`qL3Hl0qrx#`#hdiT8Tnl&=2jh;*#RCZ@`#<*)=wXG*<_G-j+hH z=en1_-P>Z9`R6gjDnfBCA1I_-Vd;;GRYyo&TC7{b{6ezv{3i}^Kp|%4uh@sPxu`g4 zS#e=ZC*-d4_m6p_ZO8`xR(p987tMOIadZw)$fe~oGGi-kN%i?NNn_9TL27@mp`2SQ zWS5g`=LJT_kk>Onrq%?NYjg?9{*>8OHvFO-k2(;&P<24)`974Sf&T~#(5Q;)S zw8Jvbq&xZ>r4Zqbkq4JlvZzZvht^th4v>dSo+Un+SY(_+y5z4zmFl>2VI86fL-RAMePBiekXaq zd)f23oW10r2L5g>*kw{W^okz@@{rgbVzU2Vxk%s^pF?Kb3Xzn4!w+``dB`UzIUU~k za?v-JR=fW_CmE@DZtuP_G8wTHo0u0*Z52Zj7SFnM*GPptoN%J9^jaX*@#yu`s<$qb zeU%e0+2WH>=!f<>PS=3EKa9#lb`Omlw7LI4bd%fY{_08& z`n>k6x9A(fIE8o~mv-p_@{kYlOmunCSW!`lLS8L-7T|l`TDEO9lFFK@~%LzuRey0{$yFS*9znzxyezXGeVM)=eNI8Mk|w%aXD_p>DCIF zuz1DPCDq35G8*4$@&0u|lv&-;o;_I{>W^1B=Y2RqD71@j3><#rhA*{YX7>#NKYBL# zp?&M%+@xbSjZ(Bg`a#<@W&`reEVtVu%8L!F@ZH<)B=udl^`FnA^nZu6^R{cv%H5OP&zYJ& zdGDO+iaW^K#0}gqTLOW0ZSw1h2httUiMnNnDzuGFerUVJiN8>sA6H29mIUqD zdnhU_=6ub4&OsM%uX22F(#bf5_@=%Jczt*F2YC+rH7BUyKnMkS!RLT>SZ%j1`zvM` z6oTOvDlPbi5{>D<%37&40+oiZw|&uuD6=_@p+j*i@Nr11gjb&h53vnC>4tal;|t{B z{Fz=p(DUX_z@%AE596PzB^pfcr2nephf75Cr&R&?%BqI9Hrgex47=AAw~FF}#+$B& zFWes5y%g^QK5}9z@b^2$S1Al1W?@Ya_SmtJm2S`Q0g);Ulk+3Aa((OIGfNEkkPP}^ d8uM6~qpA@8Zc2Y-2yeQ!H#W%oexX+D{2!BWcQ*h4 diff --git a/Fe3GeTe2_notebook.pickle b/Fe3GeTe2_notebook.pickle deleted file mode 100644 index 6fe024042121840754dc59adca6a7f0899e34dc3..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 7821 zcmcgx30M=?+77D-DwIVi6`|r4Z~+t-w8%t>phOl0lqwL0WFUc%ge(e(3MvA11cXtv z;)aHT;B_q*6*<8TcZ8}{Y+XUEdR5T2Zoq$%%%BD;>3{$K&hvbE=giDG-<$Kj=R0TS zr0cI=F_fk{mJ-gmu?Q8TiI@}_h{e&{9 zx^Y|A?prK@P>RI}g$b@A9y>wl*+RmOWhFxsig}noib{Dx0fAcyWQn2_SGG`$xhBf^ zQXVQ6qbUSFfF;J*LV-joma(NoA%PpH9)#Wm9Q7lG1e2CY_~^tCS;SnznW@9n%Otoa zxB*9+BEpnaba|dcPPY2491O^8BcfmvT!=y=xm>eG)b*Fv!)gey_cpJAAV?9Vq zh)BX&*{M20f@{njsHDL){g{L`Y?f-~L&NGbRkON_7AqV1>)KQ;%~q7P`>v%~TUiT} zLCReFzAN|Ya>QrZYF45A{}CvReu4y(5==|nAc4b!qYl?i zhBDTNV?au9jc`K=8|9;7C`+y-t{eN2Fu688W$9zoFdY-kz3ieORzq>)gHV}Nmj9mJ~Oih z-?;Q=$RIk^GZd~;UGD^G*BHak`xtlMcZzqPA>Vh+YaqyD2utsK&8;4q`DERAy0S-Q z_!o(2oB&h(V}X<>XOkNMw`9ArT`{^@sGO zXq-wYMhiV6eKF^8E(D6}C8wl9@$~FfWx9oNx6%1`RZ+y5TeVZU#+~%IMsH5yOzXYD zPO3PfxDh76#Bs{-_0j@USFIiQLW{QJPJ7#t&HGZZU6FC?pY4b^Hp%{{%68kW%sQ*I)({vqR@*Zla27P#eP^zZ1PF-Jd*zCmBn_k+DN6WT0j{?M4v zWf}S}?S}qMr+&Q670}pyiXej>G~I3rQqMTulrEcP>Kf%8H-7 zVsr1&PakG{=_dnn>|?lr2vttvaMc4>O}H+~gNscF@dQ%`*PAP2|3%_aF8s(Y>j6rd zxIPEw95UZrZ*(4G^M?Ulm2bFkY&pZhMdjU(#u?W8V9{5n&S=jXNiuhP6|l_88+7lu z8B!enD$DM;87K;lNO3Q(0uJw+PA#ad0zp@1md_R(^4>XTps1l_e=Yyb!}1XgRiK%* zyy~p)LD1D+-`0|pU3*|&sKuk%*|kx>VO|&Ebro*|e-=*ym%FM&xucDn#WoEi*Cu-Y ztdL&1!$3XuoWJ+l(wP`Wqlf2()txWZrN3)f_SMLXPRO074nfDSS^?S!ev+W$=USB; z9Y!xLQS;1X&aw+HBKIL#W~}N38|H#db0k}nO-sm{{=L=gUYi{2C$+A|Llv6dsj}0I zD3>v*p$dwPCOf@sZrBW`c?u;<)zt8oN#nkQQ^Wh-vrb<^h}YAh(bT<}_HciLUWNW& z%77YheaIl9@RECxtsEGbS)xzYh9WvemxzufA%Hu0)~%J};Nx(#py;UmA0H&;0V@xv`e#m;h-4C2kY zS?>+s6*{(rq1HgWrSf`@S0PFQ;SXHbL|m>w2I=`_L^7covC%^jB2)`ok) ztrf?L9ILmhU3}GRQXH@cc4pe9GMY@fhapS{_`Su=3Z%^{KHFNnYJ>} zm#ITP6F+ml+ar7JL~z~hZ2PlfhC#%AzgBWp`iF^^4otW?-w276A|>836$R) zK0c_qLG@ohWP%rxZM$fiBKT-R&{717;ppX29lOAhMBnSh9&B(xG7@jHib0+@@ynqT_qlfL~pygZBD{JfXYEAK1Wr1r`YcHNj-BUAN4*bv|b>FOQ7?P2H#bQ*_e)Vwvd?xU>_v6zh=+ezZ`^U37{sCU? zh&@z1O-m7!EYdkscZQ+R^iP#b9eCcupa>{Fn*5H~G^Fk;_(TL+X96nt6l!Dw=Mycb zt^e~(aJa}VK=7p;*gB?V@x-CXtfrV(Z_kGS13{Q;3S1LxsVSIatS<*X2TX!8a!_z~ zX{)ZdI2yTLdm>{UoF;~5yu+SNm4mdvEu9w*I3jOmnN8Vo$PwA)mPm|o>LG%Qf}$eB z9okI?k9aPIPOX(l?w?jy-rcMRNhc~B)2E&QEf330hQ6poR@(h6?FuHe5kb?hOUxFR zzC;$C)NRPBq|qHRHnlhq?^Wx<_`WiMq~&VB#Z3%_rhlq@YE?qvy=piekbE@xx!zyZ z%Kl$w0xI|vYKTDh>FAe(HAV3EGeK~>;K#901j+RlIbW^~1x^Q&BT*4^Y%R z=(MiduFZ7NYWu*=7Jmu;nG`>9y}ynkI%0iShxn~_BAD`Lwe3>4Fx2!K4vzu{%!l0Nx_JV_TtBO%HRF*1yBK?|H%5Zf zUe@8?wnqcm*(cX3?o0ynJN-U@KL=AG2GCktIT5V@0fcSOM?LvU_S(Sh=f{Sexm z{Qaqm-gwuZb`L6hY=q?JZ$4L2rUyK8Pd%P={C;i0`~?}dp@vB1Jx6PsA54G^!*<_} z8IFiz`h@VO6ZH|5bU%pTHS1`g!Rw9cZz_F=K+nA~@O#@;fbpi2Iqq&DGN-=Oi04=a zcHMliwlzlz=KV0#{NDW~+KAv+Bd_-g*K%aX>5UanPRPN$-?sTN3l3{1f@gX1K+)B; zJ^WLAwvLXxxS3yyP<%A`@j7OO(r4zXx199FoGFLmdg?vTh*+bg)ega@~GlsA>1s912z zL%ETP!aY?9_2~P=;3vCpVYqPu3pP$&?J#b9-}uc0+y-uGVq*B+8)J!u>TOu6CY@Fj zP_SVXirdABF-*XUl?ha>!~G4oTO@&g$j0ay9yKmt@F?ZieNOi^&ldXDusWH{^?v|N CbuW= 3.9 -- Operating System: Platform independent - -## Dependencies -The package requires the following Python packages: -- numpy -- scipy -- sisl -- netcdf4 -- openmpi -- mpi4py - -## Installation -The package uses the `hatchling` build system. Installation details will be provided in the project's repository. - -## Project Links -- Homepage: https://gitea.vo.elte.hu/et209d/grogu -- Documentation: https://gitea.vo.elte.hu/et209d/grogu -- Repository: https://gitea.vo.elte.hu/et209d/grogu -- Issue Tracker: https://gitea.vo.elte.hu/et209d/grogu - -## Development Status -The project is currently in Alpha status (Development Status :: 3 - Alpha), indicating it's in early development phase and may undergo significant changes. - -## Testing -The project uses pytest for testing with the following configurations: -- Python path is set to "src/grogu_magn/" -- Includes doctest modules -- Uses importlib import mode - -## Keywords -- DFT -- physics -- grogu -- magnetic interactions - -## License -This project is licensed under the MIT License. - -## Command Line Interface -The package provides a command-line interface through the `grogu` command: -```bash -grogu = "grogu:main" -``` - -## Project Structure -The main package code is located in the `src/grogu_magn/` directory, following modern Python packaging standards. diff --git a/out with HK = HK + Ef * SK.txt b/out with HK = HK + Ef * SK.txt deleted file mode 100644 index 321b5ad..0000000 --- a/out with HK = HK + Ef * SK.txt +++ /dev/null @@ -1,239 +0,0 @@ -Magnetic entities integrated. -Pairs integrated. -Magnetic parameters calculated. -##################################################################### GROGU OUTPUT ############################################################################# -================================================================================================================================================================ -Input file: -/Users/danielpozsar/Downloads/nojij/Fe3GeTe2/monolayer/soc/lat3_791/Fe3GeTe2.fdf -Output file: -./Fe3GeTe2_notebook.pickle -Number of nodes in the parallel cluster: 1 -================================================================================================================================================================ -Cell [Ang]: -[[ 3.79100000e+00 0.00000000e+00 0.00000000e+00] - [-1.89550000e+00 3.28310231e+00 0.00000000e+00] - [ 1.25954923e-15 2.18160327e-15 2.05700000e+01]] -================================================================================================================================================================ -DFT axis: -[0 0 1] -Quantization axis and perpendicular rotation directions: -[1 0 0] --» [array([0, 1, 0]), array([0, 0, 1])] -[0 1 0] --» [array([1, 0, 0]), array([0, 0, 1])] -[0 0 1] --» [array([1, 0, 0]), array([0, 1, 0])] -================================================================================================================================================================ -Parameters for the contour integral: -Number of k points: 15 -k point directions: xy -Ebot: -13 -Eset: 300 -Esetp: 1000 -================================================================================================================================================================ -Atomic information: ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[atom index]Element(orbitals) x [Ang] y [Ang] z [Ang] Sx Sy Sz Q Lx Ly Lz Jx Jy Jz ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(2) -7.339158738013707e-06 4.149278510690423e-06 11.657585837928032 - -[4]Fe(2) -7.326987662162937e-06 4.158274523275774e-06 8.912422537596708 - -[5]Fe(1) 1.8954667088117545 1.0943913231921656 10.285002698393109 - -================================================================================================================================================================ -Exchange [meV] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Magnetic entity1 Magnetic entity2 [i j k] d [Ang] ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -[3]Fe(2) [4]Fe(2) [0 0 0] d [Ang] 2.745163300331324 -Isotropic: 3.7364077022979965 -DMI: [1.50680376e-04 6.63325925e-06 2.14458459e-07] -Symmetric-anisotropy: [ 7.51357435e-02 -7.11950712e-07 3.59136771e-08 -7.11950712e-07 - 7.51324470e-02 3.13082991e-08 3.59136771e-08 3.13082991e-08 - -1.50268191e-01] -J: [ 3.81154345e+00 -7.11950712e-07 3.59136771e-08 -7.11950712e-07 - 3.81154015e+00 3.13082991e-08 3.59136771e-08 3.13082991e-08 - 3.58613951e+00] -Energies for debugging: -array([[ 3.58613658e-03, 1.50649068e-07, -1.50711685e-07, - 3.58823154e-03], - [ 3.58614245e-03, -6.66917293e-09, 6.59734558e-09, - 3.58823758e-03], - [ 4.03484876e-03, 9.26409172e-10, 4.97492253e-10, - 4.03484931e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00358824, 0.00403485, 0.00358614]) -Test J_xx = E(y,z) = E(z,y) -0.003588237582808751 0.004034849308800343 - -[3]Fe(2) [5]Fe(1) [0 0 0] d [Ang] 2.5835033632437767 -Isotropic: 1.3665069721027028 -DMI: [ 6.01665372e-02 -1.04123577e-01 -2.32184533e-06] -Symmetric-anisotropy: [ 0.00046905 0.00297291 -0.00117676 0.00297291 -0.00302165 -0.0006331 - -0.00117676 -0.0006331 0.0025526 ] -J: [ 1.36697602e+00 2.97290780e-03 -1.17675699e-03 2.97290780e-03 - 1.36348532e+00 -6.33095832e-04 -1.17675699e-03 -6.33095832e-04 - 1.36905957e+00] -Energies for debugging: -array([[ 1.36869655e-03, 6.07996330e-05, -5.95334413e-05, - 1.36354141e-03], - [ 1.36942260e-03, 1.05300334e-04, -1.02946820e-04, - 1.36707747e-03], - [ 1.36342924e-03, -2.97522965e-06, -2.97058595e-06, - 1.36687457e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136708, 0.00136343, 0.0013687 ]) -Test J_xx = E(y,z) = E(z,y) -0.0013670774733037568 0.0013668745662554472 - -[4]Fe(2) [5]Fe(1) [0 0 0] d [Ang] 2.583501767937866 -Isotropic: 1.3665590525954592 -DMI: [-6.01423084e-02 1.04104179e-01 -2.13066078e-06] -Symmetric-anisotropy: [ 0.0004341 0.00297286 0.00136618 0.00297286 -0.00305201 0.00083159 - 0.00136618 0.00083159 0.00261791] -J: [1.36699315e+00 2.97285931e-03 1.36618160e-03 2.97285931e-03 - 1.36350704e+00 8.31592395e-04 1.36618160e-03 8.31592395e-04 - 1.36917697e+00] -Energies for debugging: -array([[ 1.36884061e-03, -6.09739008e-05, 5.93107160e-05, - 1.36356111e-03], - [ 1.36951332e-03, -1.05470361e-04, 1.02737998e-04, - 1.36708800e-03], - [ 1.36345297e-03, -2.97498997e-06, -2.97072865e-06, - 1.36689830e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136709, 0.00136345, 0.00136884]) -Test J_xx = E(y,z) = E(z,y) -0.0013670880049909515 0.001366898297231859 - -[3]Fe(2) [5]Fe(1) [-1 -1 0] d [Ang] 2.5834973202859075 -Isotropic: 1.3666275581659808 -DMI: [-1.20163250e-01 -3.54573594e-06 6.70818449e-07] -Symmetric-anisotropy: [-4.78793224e-03 -5.40296355e-06 4.35410484e-06 -5.40296355e-06 - 2.17809565e-03 1.46954942e-03 4.35410484e-06 1.46954942e-03 - 2.60983659e-03] -J: [ 1.36183963e+00 -5.40296355e-06 4.35410484e-06 -5.40296355e-06 - 1.36880565e+00 1.46954942e-03 4.35410484e-06 1.46954942e-03 - 1.36923739e+00] -Energies for debugging: -array([[ 1.36990885e-03, -1.21632799e-04, 1.18693700e-04, - 1.36893010e-03], - [ 1.36856594e-03, -8.08368896e-10, -7.89984078e-09, - 1.36187154e-03], - [ 1.36868121e-03, 6.07378200e-09, 4.73214510e-09, - 1.36180771e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136187, 0.00136868, 0.00136991]) -Test J_xx = E(y,z) = E(z,y) -0.00136187154036587 0.0013618077114880157 - -[4]Fe(2) [5]Fe(1) [-1 -1 0] d [Ang] 2.583495745338251 -Isotropic: 1.3666148028343932 -DMI: [ 1.20163013e-01 -6.22809007e-06 4.59622024e-07] -Symmetric-anisotropy: [-4.76680949e-03 -5.40568363e-06 4.15669294e-06 -5.40568363e-06 - 2.20342253e-03 -1.46945257e-03 4.15669294e-06 -1.46945257e-03 - 2.56338696e-03] -J: [ 1.36184799e+00 -5.40568363e-06 4.15669294e-06 -5.40568363e-06 - 1.36881823e+00 -1.46945257e-03 4.15669294e-06 -1.46945257e-03 - 1.36917819e+00] -Energies for debugging: -array([[ 1.36992100e-03, 1.21632465e-04, -1.18693560e-04, - 1.36894224e-03], - [ 1.36843538e-03, 2.07139714e-09, -1.03847830e-08, - 1.36187518e-03], - [ 1.36869421e-03, 5.86530565e-09, 4.94606160e-09, - 1.36182081e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136188, 0.00136869, 0.00136992]) -Test J_xx = E(y,z) = E(z,y) -0.0013618751762170837 0.0013618208104796187 - -[3]Fe(2) [5]Fe(1) [-1 0 0] d [Ang] 2.583541444641373 -Isotropic: 1.3663664851465724 -DMI: [6.01418914e-02 1.04129666e-01 1.11624086e-06] -Symmetric-anisotropy: [ 0.00044657 -0.00296955 0.00117251 -0.00296955 -0.00303728 -0.00083176 - 0.00117251 -0.00083176 0.00259071] -J: [ 1.36681306e+00 -2.96954847e-03 1.17250768e-03 -2.96954847e-03 - 1.36332921e+00 -8.31759183e-04 1.17250768e-03 -8.31759183e-04 - 1.36895719e+00] -Energies for debugging: -array([[ 1.36864715e-03, 6.09736506e-05, -5.93101323e-05, - 1.36338366e-03], - [ 1.36926723e-03, -1.05302174e-04, 1.02957159e-04, - 1.36692240e-03], - [ 1.36327476e-03, 2.97066471e-06, 2.96843223e-06, - 1.36670371e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136692, 0.00136327, 0.00136865]) -Test J_xx = E(y,z) = E(z,y) -0.0013669224001181225 0.0013667037105967244 - -[4]Fe(2) [5]Fe(1) [-1 0 0] d [Ang] 2.5835398672184064 -Isotropic: 1.3663569094972243 -DMI: [-6.01699640e-02 -1.04105528e-01 -7.75522025e-07] -Symmetric-anisotropy: [ 0.00045093 -0.00296959 -0.00137041 -0.00296959 -0.0030254 0.00063314 - -0.00137041 0.00063314 0.00257447] -J: [ 1.36680784e+00 -2.96958730e-03 -1.37041361e-03 -2.96958730e-03 - 1.36333151e+00 6.33143106e-04 -1.37041361e-03 6.33143106e-04 - 1.36893138e+00] -Energies for debugging: -array([[ 1.36852649e-03, -6.08031071e-05, 5.95368209e-05, - 1.36338712e-03], - [ 1.36933627e-03, 1.05475941e-04, -1.02735114e-04, - 1.36691071e-03], - [ 1.36327590e-03, 2.96881177e-06, 2.97036282e-06, - 1.36670497e-03]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([0.00136691, 0.00136328, 0.00136853]) -Test J_xx = E(y,z) = E(z,y) -0.0013669107098080865 0.001366704969722615 - -[4]Fe(2) [5]Fe(1) [-2 0 0] d [Ang] 5.951322298958084 -Isotropic: -0.0034746104317761807 -DMI: [-0.00564163 0.00190664 0.00218617] -Symmetric-anisotropy: [ 1.71572275e-04 -1.64462154e-04 -6.17520368e-05 -1.64462154e-04 - 5.31812409e-05 2.26715191e-04 -6.17520368e-05 2.26715191e-04 - -2.24753515e-04] -J: [-3.30303816e-03 -1.64462154e-04 -6.17520368e-05 -1.64462154e-04 - -3.42142919e-03 2.26715191e-04 -6.17520368e-05 2.26715191e-04 - -3.69936395e-03] -Energies for debugging: -array([[-3.70537194e-06, -5.86835019e-06, 5.41491981e-06, - -3.68176164e-06], - [-3.69335596e-06, -1.84488317e-06, 1.96838725e-06, - -3.62869748e-06], - [-3.16109674e-06, 2.35062932e-06, -2.02170502e-06, - -2.97737884e-06]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([-3.62869748e-06, -3.16109674e-06, -3.70537194e-06]) -Test J_xx = E(y,z) = E(z,y) --3.628697478791367e-06 -2.977378835723539e-06 - -[4]Fe(2) [5]Fe(1) [-3 0 0] d [Ang] 9.638732176310562 -Isotropic: 0.00557283408339648 -DMI: [-0.00306844 0.00323139 0.00243485] -Symmetric-anisotropy: [ 0.00095486 -0.00049048 -0.00010491 -0.00049048 -0.00065003 -0.00015335 - -0.00010491 -0.00015335 -0.00030483] -J: [ 0.00652769 -0.00049048 -0.00010491 -0.00049048 0.0049228 -0.00015335 - -0.00010491 -0.00015335 0.00526801] -Energies for debugging: -array([[ 5.39229066e-06, -2.91508841e-06, 3.22178614e-06, - 4.70355358e-06], - [ 5.14372081e-06, -3.12647464e-06, 3.33630142e-06, - 6.55373590e-06], - [ 5.14204978e-06, 2.92533702e-06, -1.94436895e-06, - 6.50165378e-06]]) -J_ii for debugging: (check if this is the same as in calculate_exchange_tensor) -array([6.55373590e-06, 5.14204978e-06, 5.39229066e-06]) -Test J_xx = E(y,z) = E(z,y) -6.553735898142947e-06 6.501653779660145e-06 - -================================================================================================================================================================ -Runtime information: -Total runtime: 350.25778825 s ----------------------------------------------------------------------------------------------------------------------------------------------------------------- -Initial setup: 0.13192529200000003 s -Hamiltonian conversion and XC field extraction: 0.516 s -Pair and site datastructure creatrions: 0.060 s -k set cration and distribution: 0.045 s -Rotating XC potential: 0.247 s -Greens function inversion: 349.010 s -Calculate energies and magnetic components: 0.248 s diff --git a/src/grogu_magn/core.py b/src/grogu_magn/core.py index 23e4dcd..3d8591e 100644 --- a/src/grogu_magn/core.py +++ b/src/grogu_magn/core.py @@ -1,3 +1,23 @@ +# Copyright (c) [2024] [Daniel Pozsar] +# +# Permission is hereby granted, free of charge, to any person obtaining a copy +# of this software and associated documentation files (the "Software"), to deal +# in the Software without restriction, including without limitation the rights +# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell +# copies of the Software, and to permit persons to whom the Software is +# furnished to do so, subject to the following conditions: +# +# The above copyright notice and this permission notice shall be included in all +# copies or substantial portions of the Software. +# +# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR +# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, +# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE +# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER +# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, +# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE +# SOFTWARE. + import numpy as np from numpy.linalg import inv diff --git a/src/grogu_magn/grogu.py b/src/grogu_magn/grogu.py index b01b08c..21ea471 100644 --- a/src/grogu_magn/grogu.py +++ b/src/grogu_magn/grogu.py @@ -18,27 +18,6 @@ # OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE # SOFTWARE. -# Copyright (c) [2024] [Daniel Pozsar] -# -# Permission is hereby granted, free of charge, to any person obtaining a copy -# of this software and associated documentation files (the "Software"), to deal -# in the Software without restriction, including without limitation the rights -# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell -# copies of the Software, and to permit persons to whom the Software is -# furnished to do so, subject to the following conditions: -# -# The above copyright notice and this permission notice shall be included in all -# copies or substantial portions of the Software. -# -# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR -# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, -# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE -# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER -# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, -# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE -# SOFTWARE. - - import os os.environ["OMP_NUM_THREADS"] = "1" # export OMP_NUM_THREADS=1 diff --git a/src/grogu_magn/io.py b/src/grogu_magn/io.py index 0ef5529..38d2410 100644 --- a/src/grogu_magn/io.py +++ b/src/grogu_magn/io.py @@ -1,3 +1,23 @@ +# Copyright (c) [2024] [Daniel Pozsar] +# +# Permission is hereby granted, free of charge, to any person obtaining a copy +# of this software and associated documentation files (the "Software"), to deal +# in the Software without restriction, including without limitation the rights +# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell +# copies of the Software, and to permit persons to whom the Software is +# furnished to do so, subject to the following conditions: +# +# The above copyright notice and this permission notice shall be included in all +# copies or substantial portions of the Software. +# +# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR +# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, +# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE +# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER +# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, +# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE +# SOFTWARE. + from argparse import ArgumentParser from pickle import dump, load diff --git a/src/grogu_magn/magnetism.py b/src/grogu_magn/magnetism.py index 289cc55..9206577 100644 --- a/src/grogu_magn/magnetism.py +++ b/src/grogu_magn/magnetism.py @@ -1,3 +1,23 @@ +# Copyright (c) [2024] [Daniel Pozsar] +# +# Permission is hereby granted, free of charge, to any person obtaining a copy +# of this software and associated documentation files (the "Software"), to deal +# in the Software without restriction, including without limitation the rights +# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell +# copies of the Software, and to permit persons to whom the Software is +# furnished to do so, subject to the following conditions: +# +# The above copyright notice and this permission notice shall be included in all +# copies or substantial portions of the Software. +# +# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR +# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, +# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE +# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER +# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, +# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE +# SOFTWARE. + import numpy as np diff --git a/test.py b/test.py index adb7e53..21ea471 100644 --- a/test.py +++ b/test.py @@ -18,7 +18,6 @@ # OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE # SOFTWARE. - import os os.environ["OMP_NUM_THREADS"] = "1" # export OMP_NUM_THREADS=1