From 216db0b5f7665f24be364128c8fcff66a7675a7d Mon Sep 17 00:00:00 2001 From: SotaYoshida Date: Mon, 29 Jan 2024 11:23:01 +0900 Subject: [PATCH 1/3] Correction to pi-gamma and number of mesh points for SFR integrals --- src/chiEFTint/dict_LECs.jl | 2 +- src/chiEFTint/main_chiEFTint.jl | 13 +++-- src/chiEFTint/pionexchange.jl | 83 ++++++++++++++++++++++++-------- src/chiEFTint/struct_const_io.jl | 13 +++-- test/runtests.jl | 11 +++-- 5 files changed, 89 insertions(+), 33 deletions(-) diff --git a/src/chiEFTint/dict_LECs.jl b/src/chiEFTint/dict_LECs.jl index 31932a2b..7fee48ec 100644 --- a/src/chiEFTint/dict_LECs.jl +++ b/src/chiEFTint/dict_LECs.jl @@ -191,7 +191,7 @@ function dict_emn500n4lo() dLECs["D_3P2"] = 5.342585336 dLECs["hD_3S1"] = -2.949089421 dLECs["D_3S1"] = -20.793199632 - dLECs["hD_3SD1"] = 1.3545478412 + dLECs["hD_3SD1"] = 1.345478412 dLECs["D_3SD1"] = 2.176852098 dLECs["D_3D1"] = -6.01826561 dLECs["D_1D2"] = -1.545851484 diff --git a/src/chiEFTint/main_chiEFTint.jl b/src/chiEFTint/main_chiEFTint.jl index 820071b4..0fa98d03 100644 --- a/src/chiEFTint/main_chiEFTint.jl +++ b/src/chiEFTint/main_chiEFTint.jl @@ -25,8 +25,12 @@ function make_chiEFTint(;is_show=false,itnum=1,writesnt=true,nucs=[],optimizer=" BE_d_bare = Calc_Deuteron(chiEFTobj,to;io=io) @timeit to "renorm." SRG(chiEFTobj,to) BE_d_srg = Calc_Deuteron(chiEFTobj,to;io=io) - println("E(2H): bare = ",@sprintf("%8.5f", BE_d_bare), - " srg = ", @sprintf("%8.5f", BE_d_srg), " Diff.", @sprintf("%8.3e", BE_d_bare - BE_d_srg)) + if chiEFTobj.params.srg + println("E(2H): bare = ",@sprintf("%9.6f", BE_d_bare), + " srg = ", @sprintf("%9.6f", BE_d_srg), " Diff.", @sprintf("%8.3e", BE_d_bare - BE_d_srg)) + else + println("E(2H): bare = ",@sprintf("%9.6f", BE_d_bare)) + end HFdata = prepHFdata(nucs,ref,["E"],corenuc) if do_svd; target_LSJ = [[0,0,0,0],[0,2,1,1],[1,1,1,0],[2,2,0,2]]; svd_vmom(chiEFTobj,target_LSJ); end @@ -34,8 +38,8 @@ function make_chiEFTint(;is_show=false,itnum=1,writesnt=true,nucs=[],optimizer=" if write_vmom target_LSJ = [[0,0,1,1],[1,1,1,0],[1,1,0,1],[1,1,1,1],[0,0,0,0],[0,2,1,1],[3,3,1,3]] write_onshell_vmom(chiEFTobj,2,target_LSJ;label="pn"); write_onshell_vmom(chiEFTobj,3,target_LSJ;label="nn") - momplot(chiEFTobj,2,target_LSJ; fnlabel=ifelse(chiEFTobj.params.srg,"srg","bare")) - momplot(chiEFTobj,3,target_LSJ; fnlabel=ifelse(chiEFTobj.params.srg,"srg","bare")) + #momplot(chiEFTobj,2,target_LSJ; fnlabel=ifelse(chiEFTobj.params.srg,"srg","bare")) + #momplot(chiEFTobj,3,target_LSJ; fnlabel=ifelse(chiEFTobj.params.srg,"srg","bare")) end if do2n3ncalib #calibrate 2n3n LECs by HFMBPT @@ -83,6 +87,7 @@ function construct_chiEFTobj(do2n3ncalib,itnum,optimizer,MPIcomm,io,to;fn_params f_ss!(opfs[2]);f_c!(opfs[3]) ## prep. Gauss point for integrals ts, ws = Gauss_Legendre(-1,1,40) + ts, ws = Gauss_Legendre(-1,1,96) ## prep. for TBMEs infos,izs_ab,nTBME = make_sp_state(params;io=io) println(io,"# of channels 2bstate ",length(infos)," #TBME = $nTBME") diff --git a/src/chiEFTint/pionexchange.jl b/src/chiEFTint/pionexchange.jl index 6d2dd677..15c65478 100644 --- a/src/chiEFTint/pionexchange.jl +++ b/src/chiEFTint/pionexchange.jl @@ -5,7 +5,7 @@ calc. One-pion exchange potential in the momentum-space Reference: R. Machleidt, Phys. Rev. C 63 024001 (2001). """ -function OPEP(chiEFTobj,to;pigamma=false) +function OPEP(chiEFTobj,to;pigamma=true) ts = chiEFTobj.ts; ws = chiEFTobj.ws; xr = chiEFTobj.xr; V12mom = chiEFTobj.V12mom dict_pwch = chiEFTobj.dict_pwch; lsjs = chiEFTobj.lsjs; tllsj = chiEFTobj.tllsj @@ -34,13 +34,13 @@ function OPEP(chiEFTobj,to;pigamma=false) nfac = 1.0 / (xdwn * ydwn) ree = 1.0 /sqrt(ex*ey) *freg(x,y,4) f_sq!(opfs,xdwn,ydwn) - t_fc = fff * ree *hc3 * nfac + t_fc = fff *hc3 * nfac * ree tVs .= 0.0 cib_lsj_opep(t_fc,opfs,x,y,mpi02,1,J,pnrank,ts,ws,tVs,QLdict) - if pnrank == 2 ## pp/nn + if pnrank == 2 # charged pion exchange cib_lsj_opep(t_fc,opfs,x,y,mpipm2,2,J,pnrank,ts,ws,tVs,QLdict;additive=true) if pigamma - t_fc_pig = t_fc + t_fc_pig = t_fc cib_lsj_opep(t_fc,opfs,x,y,mpipm2,2,J,pnrank,ts,ws,tVs,QLdict,pigamma;additive=true,factor_pig=t_fc_pig) end end @@ -49,7 +49,7 @@ function OPEP(chiEFTobj,to;pigamma=false) tl,tlp,S,tJ = lsj[idx] if pnrank%2 == 1 && (tl+S+1)%2 != 1;continue;end V12idx = get(tdict,tllsj,-1) - if V12idx == -1;continue;end + if V12idx == -1;continue;end V12mom[V12idx][i,j] += tVs[idx] end end @@ -59,25 +59,35 @@ function OPEP(chiEFTobj,to;pigamma=false) return nothing end +""" +pi-gamma correction Eq.(4) of U. van Kolck et al., Phys. Rev. Lett. 80, 4386 (1998). + +Note that 1/(1+beta^2) is ommited, since it is also included in the OPEP contribution, +i.e. one needs to consider only the difference other than this factor. +""" function fac_pig(beta,c5=0.0) return - (1.0-beta)^2 / (2*beta^2) * log(1+beta) +(1.0+beta)/(2*beta) -2.0*c5 end """ Ref: R. Machleidt, Phys. Rev. C 63, 024001 (2001). + +For pi-gamma correction term, which is introduced in e.g. EM500 interaction, +one needs to evaluate integrals in (B11)~ terms in an explicit manner. """ -function cib_lsj_opep(fac_in,opfs,x,y,mpi2,nterm,J,pnrank,ts,ws,tVs,QLdict,pigamma=false;additive=false,factor_pig=1.0) +function cib_lsj_opep(fac_in,opfs,x,y,mpi2,nterm,J,pnrank,ts,ws,tVs,QLdict,pigamma=false;additive=false,factor_pig=1.0) x2 = x^2; y2 = y^2; z = (mpi2+x2+y2) / (2*x*y) QJ = QJm1 = 0.0 - nfac = fac_in + nfac = fac_in q2s = Float64[ ] + qfac = 1.0 if pigamma q2s = zeros(Float64,length(ts)) - nfac = ifelse(J==0,0.0,factor_pig) for (i,t) in enumerate(ts) q2 = x2 + y2 -2.0*x*y*t beta = q2 / mpi2 - q2s[i] = (fsalpha/pi) * fac_pig(beta) + nfac = fac_in * fsalpha/pi + q2s[i] = fac_pig(beta) end QJ = QL_numeric_fac(z,J,ts,ws,q2s) if J>0;QJm1=QL_numeric_fac(z,J-1,ts,ws,q2s);end @@ -85,10 +95,42 @@ function cib_lsj_opep(fac_in,opfs,x,y,mpi2,nterm,J,pnrank,ts,ws,tVs,QLdict,pigam QJ = QL(z,J,ts,ws,QLdict) if J>0;QJm1=QL(z,J-1,ts,ws,QLdict);end end + + # Using Eq.(B17)-(B20), which is valid for OPEP (not pi-gamma correction term) IJ0 = nfac * QJ - IJ1 = nfac * (z * QJ -delta(J,0)) #Eq. (B18) - IJ2 = nfac * (J*z* QJ + QJm1) /(J+1) #Eq. (B19) - IJ3 = nfac * sqrt(J/(J+1)) * (z* QJ - QJm1) #Eq. (B20) + IJ1 = nfac * (z * QJ - delta(J,0) ) + IJ2 = nfac * (J*z* QJ + QJm1) /(J+1) + IJ3 = nfac * sqrt(J/(J+1)) * (z* QJ - QJm1) + + if pigamma + QJ = QJm1 = 0.0 + for i=1:length(ws) + t = ts[i] + q2 = x2 + y2 -2.0*x*y*t + denom = (mpi2 +q2) + pj = Legendre(J,t) + QJ += pj * q2s[i] * ws[i] / denom + if J > 0 + QJm1 += Legendre(J-1,t) * q2s[i] * ws[i] / denom + end + end + QJ *= x * y + QJm1 *= x * y + IJ0 = nfac * QJ + QJ_ = 0.0 + for i=1:length(ws) + t = ts[i] + q2 = x2 + y2 -2.0*x*y*t + denom = mpi2 +q2 + pj = Legendre(J,t) + QJ_ += ts[i]* pj * q2s[i] * ws[i] / denom + end + QJ_ *= x * y + IJ1 = nfac * QJ_ + IJ2 = nfac * ( J * QJ_ + QJm1)/ (J+1) + IJ3 = nfac * sqrt(J/(J+1)) * (QJ_ - QJm1) + end + #Eq. (B28) v1 = opfs[1] * IJ0 + opfs[2] *IJ1 v2 = opfs[3] * IJ0 + opfs[4] *IJ2 @@ -103,6 +145,7 @@ function cib_lsj_opep(fac_in,opfs,x,y,mpi2,nterm,J,pnrank,ts,ws,tVs,QLdict,pigam v34 = -sqrt(J*(J+1)) * (v3-v4) v56 = sqrt(J*(J+1)) * (v5+v6) d2j1 = 1.0/(2*J+1) + if nterm == 1 phase = ifelse(pnrank==2,-1.0,1.0) tVs[1] = additive_sum(additive,tVs[1],v1 *phase) @@ -119,11 +162,11 @@ function cib_lsj_opep(fac_in,opfs,x,y,mpi2,nterm,J,pnrank,ts,ws,tVs,QLdict,pigam tVs[1] = additive_sum(additive,tVs[1],ttis * v1) tVs[2] = additive_sum(additive,tVs[2],ttit * v2) tVs[3] = additive_sum(additive,tVs[3],d2j1 * ((J+1)* (ttis*v3) + J*(ttis*v4)-(ttis*v56))) - tVs[4] = additive_sum(additive,tVs[4],d2j1 * ( J*(v3*ttis) + (J+1)*(ttis*v4) +(ttis*v56))) + tVs[4] = additive_sum(additive,tVs[4],d2j1 * ttis * ( J*v3 + (J+1)*v4 +v56) ) tVs[5] = additive_sum(additive,tVs[5],-d2j1 * ((ttis*v34)-(J+1)*(ttis*v5)+J*(ttis*v6))) tVs[6] = additive_sum(additive,tVs[6],-d2j1 * ((ttis*v34)+J*(ttis*v5)-(J+1)*(ttis*v6))) end - return nothing + return nothing end function additive_sum(TF::Bool,retv,inv) @@ -206,7 +249,7 @@ end """ tpe(chiEFTobj,to::TimerOutput) -calc. two-pion exchange terms up to N3LO(EM) or N4LO(EMN) +Function to calculate two-pion exchange terms up to N3LO(EM) or N4LO(EMN) The power conting schemes for EM/EMN are different; The ``1/M_N`` correction terms appear at NNLO in EM and at N4LO in EMN. @@ -485,7 +528,7 @@ function Vs_term(chi_order,LoopObjects,w,tw2,q2,k2,Lq,Aq,nd_mpi,r_d145,Fpi2;EMN= else f_N3LO_Vs = gA4 / (32.0 * pi^2 * Fpi4) f_N3LO_2l_Vs = -gA2 * r_d145 /(32.0*pi^2 *Fpi4) - # There is type in EM review (D.11) 5/8 -> 3/8. + # There is typo in EM review (D.11) 5/8 -> 3/8. # See also N. Kaiser, Phys. Rev. C 65 (2002) 017001. tmp_s += Lq * (k2 + 3.0*q2/8.0 + nd_mpi4 /w2) *f_N3LO_Vs tmp_s += w2 * Lq * f_N3LO_2l_Vs @@ -533,7 +576,7 @@ function Ws_term(chi_order,LoopObjects,w,q2,Lq,Aq,nd_mpi,c4,Fpi2;EMN=false,useMa f_N3LO_a_Ws = c4^2 / (96.0*pi^2 *Fpi4) tmp_s += Lq * w2 * f_N3LO_a_Ws f_N3LO_b_Ws = - c4 / (192.0*pi^2 *Fpi4) - tmp_s += Lq * ( gA2 *(16.0*nd_mpi2 + 7.0 * q2) -w2) * f_N3LO_b_Ws + tmp_s += Lq * ( gA2 *(16.0*nd_mpi2 + 7.0 * q2) -w2) * f_N3LO_b_Ws #ci/MN f_N3LO_c_Ws = 1.0/ (1536.0* pi^2 * Fpi4) f_N3LO_2l_Ws= gA4 / (2048.0 * pi^2 * Fpi6) if EMN @@ -603,7 +646,7 @@ function Vc_term(chi_order,w,tw2,q2,Lq,Aq,nd_mpi,c1,c2,c3,Fpi2,LoopObjects;EMN=f tmp_s += Lq *( brak6 * f_N3LO_f6_Vc) ### 2-loop corrections if EMN - obj = LoopObjects.n3lo; ImV = obj.ImVc; tmp_s += Integral_ImV(obj.mudomain,q2,obj,ImV) + obj = LoopObjects.n3lo; ImV = obj.ImVc; tmp_s += Integral_ImV(obj.mudomain,q2,obj,ImV) else ## EM eq.(D.9) Mm2cor = Lq *(2.0*nd_mpi2^4 / w^4 + 8.0*nd_mpi2^3 / w2 -q2^2 -2.0*nd_mpi2^2 ) + nd_mpi2^3 /(2.0 *w2) @@ -712,7 +755,7 @@ function Vls_term(chi_order,w,tw2,q2,Lq,Aq,nd_mpi,c2,Fpi2;EMN=false) end if chi_order >= 3 f_N3LO_a_Vls = c2 *gA2 /(8.0* pi^2 * Fpi4) - tmp_s += f_N3LO_a_Vls * w2 * Lq + tmp_s += f_N3LO_a_Vls * w2 * Lq f_N3LO_b_Vls = gA4 /(4.0* pi^2 * Fpi4) if !EMN tmp_s += f_N3LO_b_Vls * Lq * (11.0/32.0 * q2 + nd_mpi4 / w2) @@ -734,7 +777,7 @@ function Wls_term(chi_order,w,q2,Lq,Aq,nd_mpi,c4,Fpi2;EMN=false) if chi_order >= 3 f_N3LO_a_Wls = -c4 /(48.0*pi^2 * Fpi4) f_N3LO_b_Wls = 1.0 /(256.0*pi^2 * Fpi4) - tmp_s += f_N3LO_a_Wls*Lq * (gA2 *(8.0*nd_mpi2 + 5.0*q2) + w2) + tmp_s += f_N3LO_a_Wls*Lq * (gA2 *(8.0*nd_mpi2 + 5.0*q2) + w2) #ci/MN if !EMN tmp_s += f_N3LO_b_Wls*Lq * (16.0*gA2 *(nd_mpi2 + 3.0/8.0 *q2) + 4.0/3.0 * gA4 *(4.0*nd_mpi4 /w2 -11.0/4.0 * q2 -9.0*nd_mpi2) -w2) diff --git a/src/chiEFTint/struct_const_io.jl b/src/chiEFTint/struct_const_io.jl index dac95d53..89bd4f5a 100644 --- a/src/chiEFTint/struct_const_io.jl +++ b/src/chiEFTint/struct_const_io.jl @@ -713,6 +713,14 @@ end plot nn potential in partial wave over relative momentum space """ function momplot(chiEFTobj,pnrank,llpSJ_s;fnlabel="",ctext="",fpath="",write_hdf5=true) + cpnrank = "pp" + if pnrank == 1 + cpnrank = "pp" + elseif pnrank == 2 + cpnrank = "pn" + elseif pnrank == 3 + cpnrank = "nn" + end xr = chiEFTobj.xr_fm V12mom = chiEFTobj.V12mom tdict = chiEFTobj.dict_pwch[pnrank] @@ -735,7 +743,6 @@ function momplot(chiEFTobj,pnrank,llpSJ_s;fnlabel="",ctext="",fpath="",write_hdf end for i =1:n_mesh; tv[i] = V[i,i]; end if fpath != ""; tfdat = [xf,yfs[vidx]];end - #pw_plt(tx,xr,V,tv,pnrank,fnlabel;fdat=tfdat) if write_hdf5 fn = "vmom_"*fnlabel*"_chiEFT_"*tx*"_"*string(pnrank)*".h5" h5open(fn,"w") do file @@ -753,8 +760,8 @@ end function momplot_from_file(pnrank,tx) fn_bare = "vmom_bare_chiEFT_"*tx*"_"*string(pnrank)*".h5" fn_srg = "vmom_srg_chiEFT_"*tx*"_"*string(pnrank)*".h5" - if !isfile(fn_bare);println("file $fn_bare is not found");return nothing;end - if !isfile(fn_srg);println("file $fn_srg is not found");return nothing;end + if !isfile(fn_bare);println("file $fn_bare is not found, you need to run srg=false calculations to create 'vmom' files");return nothing;end + if !isfile(fn_srg);println("file $fn_srg is not found , you need to run srg=true calculations to create 'vmom' files");return nothing;end file = h5open(fn_bare,"r"); x_b = read(file,"x"); V_b = read(file,"V"); V_b_diag = diag(V_b); close(file) file = h5open(fn_srg,"r"); x_s = read(file,"x"); V_s = read(file,"V"); V_s_diag = diag(V_s); close(file) tls = ["pp", "pn", "nn"] diff --git a/test/runtests.jl b/test/runtests.jl index 45ee4738..220725d2 100644 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -1,6 +1,7 @@ -using NuclearToolkit +# using NuclearToolkit using Test - +include("src/NuclearToolkit.jl") +using .NuclearToolkit @testset "NuclearToolkit.jl" begin @testset "generate NN potential" begin @@ -14,7 +15,7 @@ using Test sntf = "tbme_em500n3lo_barehw20emax2.snt.bin" ## HF-MBPT from snt/snt.bin @testset "HFMBPT results under bare EM500,hw20,e2,nmesh50" begin - Eref = [1.47561,-5.80114,0.39324] + Eref = [1.493, -5.805, 0.395] HFobj1 = hf_main(nucs,sntf,hw,emax;return_obj=true,verbose=true) Es1 = [HFobj1.E0, HFobj1.EMP2, HFobj1.EMP3] println("Eref $Eref") @@ -35,7 +36,7 @@ using Test @test ((HFobj1.E0-HFobj2.E0)^2 + (HFobj1.EMP2-HFobj2.EMP2)^2 + (HFobj1.EMP3-HFobj2.EMP3)^2) < 1.e-6 end end - Eref = -4.06623902 + Eref = -4.05224909 @testset "IMSRG results under bare EM500,hw20,e2,nmesh50" begin IMSRGobj = hf_main(nucs,sntf,hw,emax;doIMSRG=true,return_obj=true) Es = IMSRGobj.H.zerobody[1] @@ -47,7 +48,7 @@ using Test @test abs(Eref - Es[1]) < 1.e-6 end @testset "shell model calculation" begin - Eref = [ -10.737, -8.426] + Eref = [ -10.720, -8.410] vs_sntf = "vsimsrg_p-shell_coreHe4refHe4_He4_hw20e2_Delta0.0.snt"; n_eigen=2;targetJ=[] Ens = main_sm(vs_sntf,"Be8",n_eigen,targetJ) @test ((Eref[1]-Ens[1])^2 + (Eref[2] - Ens[2])^2) < 1.e-6 From 5c4140b49462d64b7763edc39ef95be84f1aa367 Mon Sep 17 00:00:00 2001 From: SotaYoshida Date: Mon, 29 Jan 2024 11:31:16 +0900 Subject: [PATCH 2/3] test code updated --- test/runtests.jl | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/test/runtests.jl b/test/runtests.jl index 220725d2..c25ed43f 100644 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -1,7 +1,5 @@ -# using NuclearToolkit +using NuclearToolkit using Test -include("src/NuclearToolkit.jl") -using .NuclearToolkit @testset "NuclearToolkit.jl" begin @testset "generate NN potential" begin From cad273ec727e77021ecc64995c09dabfb1a98ac2 Mon Sep 17 00:00:00 2001 From: SotaYoshida Date: Mon, 29 Jan 2024 02:34:50 +0000 Subject: [PATCH 3/3] Update log_sample_script.txt --- example/log_sample_script.txt | 555 +++++++++++++++++----------------- 1 file changed, 277 insertions(+), 278 deletions(-) diff --git a/example/log_sample_script.txt b/example/log_sample_script.txt index e99ac81c..9150331d 100644 --- a/example/log_sample_script.txt +++ b/example/log_sample_script.txt @@ -1,294 +1,293 @@ size of dWS (jmax 9 lmax 40 e2max 8 Nnmax 20): dtri 4.46 MB dcgm0 1.11 MB d6j_int 1.11 MB d6j_lj 0.28 MB d9j_lsj 1.76 MB dictHOB 0.44 MB -E(2H): bare = -2.23001 srg = -2.23001 Diff.1.487e-10 +E(2H): bare = -2.224578 srg = -2.224578 Diff.1.484e-10 size of dWS (jmax 9 lmax 40 e2max 8 Nnmax 20): dtri 4.46 MB dcgm0 1.11 MB d6j_int 1.11 MB d6j_lj 0.28 MB d9j_lsj 1.76 MB dictHOB 0.44 MB -target: O16 Ref. => Z=8 N=8 E: -132.9506 = E1b 363.70 + E2b -496.65 ( -55.16 -367.32 -74.18) -EMP2 -23.70628 1b -0.000 pp -1.437 pn -20.815 nn -1.454 -EMP3 0.37645 pp -0.111 hh -1.895 ph 2.383 -E_HF -132.95058 E_MBPT(3) = -156.2804 Eexp: -127.619 +target: O16 Ref. => Z=8 N=8 E: -132.7334 = E1b 363.65 + E2b -496.38 ( -55.15 -367.06 -74.17) +EMP2 -23.70205 1b -0.000 pp -1.437 pn -20.810 nn -1.455 +EMP3 0.37877 pp -0.110 hh -1.895 ph 2.384 +E_HF -132.73341 E_MBPT(3) = -156.0567 Eexp: -127.619 parameters in optional_parameters.jl will be used. def-by-run d6j_lj done! 4191 step: s E0 ||Omega_1|| ||Omega_2|| ||Eta_1|| ||Eta_2|| Ncomm. nwritten - 0 0.000 -132.95057781 0.000000e+00 0.000000e+00 1.884675e-16 1.410292e+00 0 0 - 1 0.500 -150.64837609 9.423374e-17 7.051459e-01 4.457285e-02 6.465559e-01 10 0 - 2 1.000 -154.78704838 2.228643e-02 3.232779e-01 5.387093e-02 3.441643e-01 18 1 - 3 1.500 -156.02709507 2.693546e-02 1.720822e-01 5.200791e-02 1.968159e-01 25 2 - 4 2.000 -156.46234611 5.291574e-02 2.697514e-01 4.609815e-02 1.192463e-01 34 2 - 5 2.500 -156.63786066 2.304908e-02 5.962317e-02 3.915110e-02 7.603387e-02 40 3 - 6 3.000 -156.71754528 4.261375e-02 9.740456e-02 3.246552e-02 5.082054e-02 47 3 - 7 3.500 -156.75734646 5.883469e-02 1.223106e-01 2.655233e-02 3.554179e-02 55 3 - 8 4.000 -156.77877263 7.209319e-02 1.393224e-01 2.156017e-02 2.604640e-02 63 3 - 9 4.500 -156.79104484 8.284970e-02 1.513341e-01 1.745834e-02 2.007233e-02 71 3 - 10 5.000 -156.79848127 9.154892e-02 1.600942e-01 1.414576e-02 1.629448e-02 79 3 - 11 5.500 -156.80324373 9.858457e-02 1.666945e-01 1.150104e-02 1.388250e-02 87 3 - 12 6.000 -156.80646858 1.042894e-01 1.718341e-01 9.405852e-03 1.230019e-02 95 3 - 13 6.500 -156.80877527 1.089366e-01 1.759706e-01 7.754598e-03 1.120726e-02 104 3 - 14 7.000 -156.81051145 1.127468e-01 1.794075e-01 6.457349e-03 1.039706e-02 113 3 - 15 7.500 -156.81187735 1.158958e-01 1.823489e-01 5.439592e-03 9.749738e-03 122 3 - 16 8.000 -156.81299135 1.185223e-01 1.849332e-01 4.640660e-03 9.198211e-03 131 3 - 17 8.500 -156.81392512 1.207286e-01 1.872546e-01 4.015938e-03 8.705218e-03 139 3 - 18 9.000 -156.81472443 1.226094e-01 1.893775e-01 3.520493e-03 8.251996e-03 147 3 - 19 9.500 -156.81541866 1.242303e-01 1.913467e-01 3.124629e-03 7.827798e-03 155 3 - 20 10.000 -156.81602782 1.256424e-01 1.931932e-01 2.805005e-03 7.426899e-03 163 3 - 21 10.500 -156.81656612 1.268858e-01 1.949382e-01 2.543557e-03 7.046196e-03 171 3 - 22 11.000 -156.81704414 1.279916e-01 1.965967e-01 2.326485e-03 6.683924e-03 179 3 - 23 11.500 -156.81747005 1.289842e-01 1.981790e-01 2.143354e-03 6.338987e-03 187 3 - 24 12.000 -156.81785041 1.298828e-01 1.996926e-01 1.986340e-03 6.010600e-03 195 3 - 25 12.500 -156.81819065 1.307022e-01 2.011427e-01 1.849605e-03 5.698119e-03 203 3 - 26 13.000 -156.81849533 1.314544e-01 2.025332e-01 1.728808e-03 5.400952e-03 211 3 - 27 13.500 -156.81876840 1.321487e-01 2.038672e-01 1.620718e-03 5.118518e-03 219 3 - 28 14.000 -156.81901328 1.327926e-01 2.051470e-01 1.522924e-03 4.850238e-03 227 3 - 29 14.500 -156.81923296 1.333922e-01 2.063745e-01 1.433620e-03 4.595530e-03 235 3 - 30 15.000 -156.81943009 1.339523e-01 2.075515e-01 1.351441e-03 4.353810e-03 243 3 - 31 15.500 -156.81960704 1.344770e-01 2.086796e-01 1.275347e-03 4.124497e-03 251 3 - 32 16.000 -156.81976590 1.349696e-01 2.097602e-01 1.204533e-03 3.907021e-03 259 3 - 33 16.500 -156.81990854 1.354330e-01 2.107947e-01 1.138373e-03 3.700822e-03 267 3 - 34 17.000 -156.82003317 5.691865e-04 1.850411e-03 1.074658e-03 3.506997e-03 272 4 - 35 17.500 -156.82014499 1.106509e-03 3.603905e-03 1.015349e-03 3.323020e-03 277 4 - 36 18.000 -156.82024531 1.614168e-03 5.265407e-03 9.598947e-04 3.148469e-03 282 4 - 37 18.500 -156.82033533 2.094090e-03 6.839628e-03 9.078782e-04 2.982918e-03 287 4 - 38 19.000 -156.82041609 2.547997e-03 8.331069e-03 8.589724e-04 2.825941e-03 292 4 - 39 19.500 -156.82048854 2.977444e-03 9.744018e-03 8.129099e-04 2.677124e-03 297 4 - 40 20.000 -156.82055356 3.383857e-03 1.108256e-02 7.694672e-04 2.536064e-03 302 4 - 41 20.500 -156.82061188 3.768545e-03 1.235056e-02 7.284533e-04 2.402376e-03 307 4 - 42 21.000 -156.82066421 4.132725e-03 1.355172e-02 6.897020e-04 2.275686e-03 312 4 - 43 21.500 -156.82071117 4.477530e-03 1.468954e-02 6.530664e-04 2.155638e-03 317 4 - 44 22.000 -156.82075329 4.804017e-03 1.576733e-02 6.184151e-04 2.041891e-03 323 4 - 45 22.500 -156.82079108 5.113180e-03 1.678825e-02 5.856293e-04 1.934120e-03 329 4 - 46 23.000 -156.82082498 5.405952e-03 1.775529e-02 5.546002e-04 1.832017e-03 335 4 - 47 23.500 -156.82085540 5.683211e-03 1.867127e-02 5.252279e-04 1.735286e-03 341 4 - 48 24.000 -156.82088269 5.945786e-03 1.953889e-02 4.974198e-04 1.643648e-03 347 4 - 49 24.500 -156.82090717 6.194459e-03 2.036069e-02 4.710897e-04 1.556836e-03 353 4 - 50 25.000 -156.82092913 6.429969e-03 2.113909e-02 4.461569e-04 1.474599e-03 359 4 - 51 25.500 -156.82094883 6.653014e-03 2.187636e-02 4.225460e-04 1.396697e-03 365 4 - 52 26.000 -156.82096651 6.864257e-03 2.257469e-02 4.001858e-04 1.322902e-03 371 4 - 53 26.500 -156.82098237 7.064321e-03 2.323613e-02 3.790095e-04 1.253000e-03 377 4 - 54 27.000 -156.82099659 7.253799e-03 2.386261e-02 3.589540e-04 1.186785e-03 383 4 - 55 27.500 -156.82100935 7.433251e-03 2.445599e-02 3.399597e-04 1.124063e-03 389 4 - 56 28.000 -156.82102080 7.603207e-03 2.501800e-02 3.219703e-04 1.064652e-03 395 4 - 57 28.500 -156.82103107 7.764171e-03 2.555032e-02 3.049325e-04 1.008377e-03 401 4 - 58 29.000 -156.82104028 7.916617e-03 2.605449e-02 2.887961e-04 9.550725e-04 407 4 - 59 29.500 -156.82104854 8.060996e-03 2.653201e-02 2.735132e-04 9.045825e-04 413 4 - 60 30.000 -156.82105595 8.197734e-03 2.698429e-02 2.590387e-04 8.567587e-04 419 4 - 61 30.500 -156.82106260 8.327237e-03 2.741266e-02 2.453300e-04 8.114605e-04 425 4 - 62 31.000 -156.82106857 8.449887e-03 2.781838e-02 2.323465e-04 7.685550e-04 431 4 - 63 31.500 -156.82107392 8.566046e-03 2.820265e-02 2.200499e-04 7.279160e-04 437 4 - 64 32.000 -156.82107872 8.676057e-03 2.856660e-02 2.084039e-04 6.894239e-04 443 4 - 65 32.500 -156.82108302 8.780246e-03 2.891131e-02 1.973741e-04 6.529656e-04 449 4 - 66 33.000 -156.82108688 8.878922e-03 2.923778e-02 1.869279e-04 6.184337e-04 455 4 - 67 33.500 -156.82109035 8.972374e-03 2.954699e-02 1.770345e-04 5.857267e-04 461 4 - 68 34.000 -156.82109345 9.060881e-03 2.983985e-02 1.676646e-04 5.547481e-04 467 4 - 69 34.500 -156.82109624 9.144704e-03 3.011722e-02 1.587905e-04 5.254069e-04 473 4 - 70 35.000 -156.82109874 9.224090e-03 3.037991e-02 1.503861e-04 4.976166e-04 479 4 - 71 35.500 -156.82110098 9.299275e-03 3.062872e-02 1.424264e-04 4.712953e-04 485 4 - 72 36.000 -156.82110300 9.370480e-03 3.086436e-02 1.348880e-04 4.463654e-04 491 4 - 73 36.500 -156.82110480 9.437916e-03 3.108754e-02 1.277486e-04 4.227536e-04 497 4 - 74 37.000 -156.82110642 9.501784e-03 3.129891e-02 1.209870e-04 4.003901e-04 503 4 - 75 37.500 -156.82110787 9.562271e-03 3.149910e-02 1.145833e-04 3.792090e-04 509 4 - 76 38.000 -156.82110917 9.619556e-03 3.168870e-02 1.085185e-04 3.591479e-04 515 4 - 77 38.500 -156.82111034 9.673809e-03 3.186827e-02 1.027747e-04 3.401476e-04 521 4 - 78 39.000 -156.82111139 9.725191e-03 3.203834e-02 9.733496e-05 3.221521e-04 527 4 - 79 39.500 -156.82111233 9.773854e-03 3.219941e-02 9.218310e-05 3.051083e-04 533 4 - 80 40.000 -156.82111317 9.819940e-03 3.235197e-02 8.730392e-05 2.889658e-04 539 4 - 81 40.500 -156.82111393 9.863588e-03 3.249645e-02 8.268300e-05 2.736771e-04 545 4 - 82 41.000 -156.82111461 9.904925e-03 3.263328e-02 7.830665e-05 2.591970e-04 551 4 - 83 41.500 -156.82111522 9.944074e-03 3.276288e-02 7.416194e-05 2.454828e-04 557 4 - 84 42.000 -156.82111576 9.981151e-03 3.288562e-02 7.023662e-05 2.324940e-04 563 4 - 85 42.500 -156.82111625 1.001627e-02 3.300186e-02 6.651905e-05 2.201922e-04 569 4 - 86 43.000 -156.82111669 1.004952e-02 3.311196e-02 6.299826e-05 2.085412e-04 575 4 - 87 43.500 -156.82111708 1.008102e-02 3.321622e-02 5.966383e-05 1.975066e-04 581 4 - 88 44.000 -156.82111744 1.011085e-02 3.331498e-02 5.650588e-05 1.870556e-04 587 4 - 89 44.500 -156.82111775 1.013910e-02 3.340850e-02 5.351509e-05 1.771575e-04 593 4 - 90 45.000 -156.82111804 1.016585e-02 3.349708e-02 5.068259e-05 1.677831e-04 599 4 - 91 45.500 -156.82111829 1.019119e-02 3.358097e-02 4.800002e-05 1.589046e-04 605 4 - 92 46.000 -156.82111852 1.021519e-02 3.366042e-02 4.545944e-05 1.504959e-04 611 4 - 93 46.500 -156.82111873 1.023792e-02 3.373567e-02 4.305333e-05 1.425320e-04 617 4 - 94 47.000 -156.82111891 1.025944e-02 3.380693e-02 4.077457e-05 1.349895e-04 623 4 - 95 47.500 -156.82111907 1.027983e-02 3.387442e-02 3.861643e-05 1.278461e-04 629 4 - 96 48.000 -156.82111922 1.029913e-02 3.393835e-02 3.657252e-05 1.210806e-04 635 4 - 97 48.500 -156.82111936 1.031742e-02 3.399889e-02 3.463679e-05 1.146731e-04 641 4 - 98 49.000 -156.82111947 1.033473e-02 3.405622e-02 3.280352e-05 1.086046e-04 647 4 - 99 49.500 -156.82111958 1.035113e-02 3.411052e-02 3.106728e-05 1.028572e-04 653 4 - 100 50.000 -156.82111968 1.036667e-02 3.416195e-02 2.942294e-05 9.741393e-05 659 4 - 101 50.500 -156.82111976 1.038138e-02 3.421066e-02 2.786563e-05 9.225868e-05 665 4 - 102 51.000 -156.82111984 1.039531e-02 3.425678e-02 2.639075e-05 8.737623e-05 671 4 - 103 51.500 -156.82111991 1.040850e-02 3.430047e-02 2.499393e-05 8.275213e-05 677 4 - 104 52.000 -156.82111997 1.042100e-02 3.434185e-02 2.367105e-05 7.837272e-05 683 4 - 105 52.500 -156.82112003 1.043283e-02 3.438103e-02 2.241818e-05 7.422506e-05 689 4 - 106 53.000 -156.82112008 1.044404e-02 3.441814e-02 2.123163e-05 7.029688e-05 695 4 - 107 53.500 -156.82112012 1.045465e-02 3.445329e-02 2.010788e-05 6.657657e-05 701 4 - 108 54.000 -156.82112016 1.046471e-02 3.448658e-02 1.904361e-05 6.305314e-05 707 4 - 109 54.500 -156.82112020 1.047423e-02 3.451811e-02 1.803567e-05 5.971616e-05 713 4 - 110 55.000 -156.82112023 1.048324e-02 3.454796e-02 1.708108e-05 5.655577e-05 719 4 - 111 55.500 -156.82112026 1.049178e-02 3.457624e-02 1.617701e-05 5.356263e-05 725 4 - 112 56.000 -156.82112028 1.049987e-02 3.460302e-02 1.532079e-05 5.072788e-05 731 4 - 113 56.500 -156.82112031 1.050753e-02 3.462839e-02 1.450990e-05 4.804315e-05 737 4 - 114 57.000 -156.82112033 1.051479e-02 3.465241e-02 1.374192e-05 4.550050e-05 743 4 - 115 57.500 -156.82112035 1.052166e-02 3.467516e-02 1.301459e-05 4.309241e-05 749 4 - 116 58.000 -156.82112036 1.052816e-02 3.469670e-02 1.232575e-05 4.081176e-05 755 4 - 117 58.500 -156.82112038 1.053433e-02 3.471711e-02 1.167338e-05 3.865181e-05 761 4 - 118 59.000 -156.82112039 1.054016e-02 3.473643e-02 1.105553e-05 3.660617e-05 767 4 - 119 59.500 -156.82112041 1.054569e-02 3.475474e-02 1.047038e-05 3.466878e-05 773 4 - 120 60.000 -156.82112042 1.055092e-02 3.477207e-02 9.916210e-06 3.283393e-05 779 4 - 121 60.500 -156.82112043 1.055588e-02 3.478849e-02 9.391367e-06 3.109619e-05 785 4 - 122 61.000 -156.82112043 1.056058e-02 3.480403e-02 8.894303e-06 2.945041e-05 791 4 - 123 61.500 -156.82112044 1.056502e-02 3.481876e-02 8.423548e-06 2.789173e-05 797 4 - 124 62.000 -156.82112045 1.056923e-02 3.483271e-02 7.977709e-06 2.641554e-05 803 4 - 125 62.500 -156.82112046 1.057322e-02 3.484591e-02 7.555467e-06 2.501748e-05 809 4 - 126 63.000 -156.82112046 1.057700e-02 3.485842e-02 7.155573e-06 2.369341e-05 815 4 - 127 63.500 -156.82112047 1.058058e-02 3.487027e-02 6.776845e-06 2.243942e-05 821 4 - 128 64.000 -156.82112047 1.058397e-02 3.488149e-02 6.418163e-06 2.125179e-05 827 4 - 129 64.500 -156.82112048 1.058717e-02 3.489211e-02 6.078464e-06 2.012701e-05 833 4 - 130 65.000 -156.82112048 1.059021e-02 3.490218e-02 5.756745e-06 1.906177e-05 839 4 - 131 65.500 -156.82112048 1.059309e-02 3.491171e-02 5.452054e-06 1.805290e-05 845 4 - 132 66.000 -156.82112049 1.059582e-02 3.492073e-02 5.163490e-06 1.709743e-05 851 4 - 133 66.500 -156.82112049 1.059840e-02 3.492928e-02 4.890199e-06 1.619253e-05 857 4 - 134 67.000 -156.82112049 1.060084e-02 3.493738e-02 4.631372e-06 1.533551e-05 863 4 - 135 67.500 -156.82112049 1.060316e-02 3.494505e-02 4.386245e-06 1.452386e-05 869 4 - 136 68.000 -156.82112049 1.060535e-02 3.495231e-02 4.154091e-06 1.375516e-05 875 4 - 137 68.500 -156.82112050 1.060743e-02 3.495919e-02 3.934225e-06 1.302715e-05 881 4 - 138 69.000 -156.82112050 1.060940e-02 3.496570e-02 3.725996e-06 1.233767e-05 887 4 - 139 69.500 -156.82112050 1.061126e-02 3.497187e-02 3.528788e-06 1.168467e-05 893 4 - 140 70.000 -156.82112050 1.061302e-02 3.497771e-02 3.342018e-06 1.106624e-05 899 4 - 141 70.500 -156.82112050 1.061469e-02 3.498324e-02 3.165133e-06 1.048054e-05 905 4 - 142 71.000 -156.82112050 1.061628e-02 3.498848e-02 2.997610e-06 9.925842e-06 911 4 - 143 71.500 -156.82112050 1.061777e-02 3.499345e-02 2.838954e-06 9.400498e-06 917 4 - 144 72.000 -156.82112050 1.061919e-02 3.499815e-02 2.688695e-06 8.902960e-06 923 4 - 145 72.500 -156.82112050 1.062054e-02 3.500260e-02 2.546389e-06 8.431754e-06 929 4 - 146 73.000 -156.82112051 1.062181e-02 3.500681e-02 2.411615e-06 7.985487e-06 935 4 - 147 73.500 -156.82112051 1.062302e-02 3.501081e-02 2.283974e-06 7.562840e-06 941 4 - 148 74.000 -156.82112051 1.062416e-02 3.501459e-02 2.163089e-06 7.162562e-06 947 4 - 149 74.500 -156.82112051 1.062524e-02 3.501817e-02 2.048602e-06 6.783469e-06 953 4 - 150 75.000 -156.82112051 1.062626e-02 3.502156e-02 1.940174e-06 6.424440e-06 959 4 - 151 75.500 -156.82112051 1.062723e-02 3.502477e-02 1.837486e-06 6.084413e-06 965 4 - 152 76.000 -156.82112051 1.062815e-02 3.502781e-02 1.740232e-06 5.762383e-06 971 4 - 153 76.500 -156.82112051 1.062902e-02 3.503070e-02 1.648126e-06 5.457397e-06 977 4 - 154 77.000 -156.82112051 1.062985e-02 3.503342e-02 1.560895e-06 5.168553e-06 983 4 - 155 77.500 -156.82112051 1.063063e-02 3.503601e-02 1.478281e-06 4.894996e-06 989 4 - 156 78.000 -156.82112051 1.063137e-02 3.503846e-02 1.400039e-06 4.635918e-06 995 4 - 157 78.500 -156.82112051 1.063207e-02 3.504077e-02 1.325938e-06 4.390552e-06 1001 4 - 158 79.000 -156.82112051 1.063273e-02 3.504297e-02 1.255760e-06 4.158172e-06 1007 4 - 159 79.500 -156.82112051 1.063336e-02 3.504505e-02 1.189295e-06 3.938092e-06 1013 4 - 160 80.000 -156.82112051 1.063395e-02 3.504702e-02 1.126349e-06 3.729660e-06 1019 4 - 161 80.500 -156.82112051 1.063452e-02 3.504888e-02 1.066734e-06 3.532259e-06 1025 4 - 162 81.000 -156.82112051 1.063505e-02 3.505065e-02 1.010275e-06 3.345306e-06 1031 4 - 163 81.500 -156.82112051 1.063555e-02 3.505232e-02 9.568032e-07 3.168249e-06 1037 4 - 164 82.000 -156.82112051 1.063603e-02 3.505390e-02 9.061620e-07 3.000562e-06 1043 4 - 165 82.500 -156.82112051 1.063648e-02 3.505540e-02 8.582011e-07 2.841750e-06 1049 4 - 166 83.000 -156.82112051 1.063691e-02 3.505683e-02 8.127787e-07 2.691344e-06 1055 4 - 167 83.500 -156.82112051 1.063732e-02 3.505817e-02 7.697603e-07 2.548899e-06 1061 4 - 168 84.000 -156.82112051 1.063771e-02 3.505945e-02 7.290189e-07 2.413992e-06 1067 4 - 169 84.500 -156.82112051 1.063807e-02 3.506065e-02 6.904337e-07 2.286226e-06 1073 4 - 170 85.000 -156.82112051 1.063841e-02 3.506180e-02 6.538908e-07 2.165222e-06 1079 4 - 171 85.500 -156.82112051 1.063874e-02 3.506288e-02 6.192820e-07 2.050623e-06 1085 4 - 172 86.000 -156.82112051 1.063905e-02 3.506390e-02 5.865050e-07 1.942089e-06 1091 4 - 173 86.500 -156.82112051 1.063934e-02 3.506487e-02 5.554627e-07 1.839299e-06 1097 4 - 174 87.000 -156.82112051 1.063962e-02 3.506579e-02 5.260635e-07 1.741950e-06 1103 4 - 175 87.500 -156.82112051 1.063989e-02 3.506667e-02 4.982203e-07 1.649753e-06 1109 4 - 176 88.000 -156.82112051 1.064013e-02 3.506749e-02 4.718507e-07 1.562436e-06 1115 4 - 177 88.500 -156.82112051 1.064037e-02 3.506827e-02 4.468768e-07 1.479741e-06 1121 4 - 178 89.000 -156.82112051 1.064059e-02 3.506901e-02 4.232248e-07 1.401422e-06 1127 4 - 179 89.500 -156.82112051 1.064081e-02 3.506971e-02 4.008245e-07 1.327248e-06 1133 4 - 180 90.000 -156.82112051 1.064101e-02 3.507038e-02 3.796099e-07 1.257000e-06 1139 4 - 181 90.500 -156.82112051 1.064120e-02 3.507100e-02 3.595181e-07 1.190471e-06 1145 4 - 182 91.000 -156.82112051 1.064138e-02 3.507160e-02 3.404897e-07 1.127462e-06 1151 4 - 183 91.500 -156.82112051 1.064155e-02 3.507216e-02 3.224685e-07 1.067789e-06 1157 4 - 184 92.000 -156.82112051 1.064171e-02 3.507270e-02 3.054010e-07 1.011273e-06 1163 4 - 185 92.500 -156.82112051 1.064186e-02 3.507320e-02 2.892369e-07 9.577492e-07 1169 4 - 186 93.000 -156.82112051 1.064200e-02 3.507368e-02 2.739283e-07 9.070580e-07 1175 4 + 0 0.000 -132.73340718 0.000000e+00 0.000000e+00 1.050065e-16 1.410541e+00 0 0 + 1 0.500 -150.42715474 5.250326e-17 7.052707e-01 4.464929e-02 6.465989e-01 10 0 + 2 1.000 -154.56414557 2.232464e-02 3.232995e-01 5.395684e-02 3.441616e-01 18 1 + 3 1.500 -155.80358359 2.697842e-02 1.720808e-01 5.208556e-02 1.968026e-01 25 2 + 4 2.000 -156.23862273 5.299755e-02 2.697434e-01 4.616261e-02 1.192332e-01 34 2 + 5 2.500 -156.41406698 2.308130e-02 5.961658e-02 3.920241e-02 7.602368e-02 40 3 + 6 3.000 -156.49373083 4.267163e-02 9.739283e-02 3.250542e-02 5.081370e-02 47 3 + 7 3.500 -156.53352793 5.891257e-02 1.222953e-01 2.658294e-02 3.553780e-02 55 3 + 8 4.000 -156.55495503 7.218644e-02 1.393048e-01 2.158347e-02 2.604452e-02 63 3 + 9 4.500 -156.56722918 8.295468e-02 1.513152e-01 1.747599e-02 2.007182e-02 71 3 + 10 5.000 -156.57466729 9.166283e-02 1.600747e-01 1.415906e-02 1.629469e-02 79 3 + 11 5.500 -156.57943088 9.870525e-02 1.666746e-01 1.151103e-02 1.388297e-02 87 3 + 12 6.000 -156.58265638 1.044152e-01 1.718142e-01 9.413302e-03 1.230062e-02 95 3 + 13 6.500 -156.58496333 1.090663e-01 1.759507e-01 7.760114e-03 1.120750e-02 104 3 + 14 7.000 -156.58669949 1.128795e-01 1.793877e-01 6.461390e-03 1.039707e-02 113 3 + 15 7.500 -156.58806521 1.160306e-01 1.823291e-01 5.442510e-03 9.749531e-03 122 3 + 16 8.000 -156.58917892 1.186587e-01 1.849134e-01 4.642728e-03 9.197811e-03 131 3 + 17 8.500 -156.59011234 1.208663e-01 1.872347e-01 4.017372e-03 8.704660e-03 139 3 + 18 9.000 -156.59091127 1.227480e-01 1.893574e-01 3.521453e-03 8.251313e-03 147 3 + 19 9.500 -156.59160512 1.243695e-01 1.913264e-01 3.125237e-03 7.827018e-03 155 3 + 20 10.000 -156.59221389 1.257821e-01 1.931726e-01 2.805356e-03 7.426046e-03 163 3 + 21 10.500 -156.59275183 1.270257e-01 1.949173e-01 2.543721e-03 7.045289e-03 171 3 + 22 11.000 -156.59322950 1.281317e-01 1.965754e-01 2.326516e-03 6.682978e-03 179 3 + 23 11.500 -156.59365508 1.291245e-01 1.981573e-01 2.143291e-03 6.338013e-03 187 3 + 24 12.000 -156.59403514 1.300231e-01 1.996705e-01 1.986210e-03 6.009609e-03 195 3 + 25 12.500 -156.59437510 1.308425e-01 2.011201e-01 1.849430e-03 5.697118e-03 203 3 + 26 13.000 -156.59467952 1.315946e-01 2.025102e-01 1.728602e-03 5.399946e-03 211 3 + 27 13.500 -156.59495235 1.322888e-01 2.038437e-01 1.620491e-03 5.117513e-03 219 3 + 28 14.000 -156.59519700 1.329326e-01 2.051230e-01 1.522684e-03 4.849237e-03 227 3 + 29 14.500 -156.59541648 1.335321e-01 2.063501e-01 1.433373e-03 4.594536e-03 235 3 + 30 15.000 -156.59561343 1.340920e-01 2.075266e-01 1.351190e-03 4.352826e-03 243 3 + 31 15.500 -156.59579021 1.346166e-01 2.086542e-01 1.275094e-03 4.123526e-03 251 3 + 32 16.000 -156.59594891 1.351091e-01 2.097343e-01 1.204281e-03 3.906064e-03 259 3 + 33 16.500 -156.59609140 1.355723e-01 2.107684e-01 1.138123e-03 3.699880e-03 267 3 + 34 17.000 -156.59621591 5.690614e-04 1.849940e-03 1.074412e-03 3.506075e-03 272 4 + 35 17.500 -156.59632761 1.106260e-03 3.602974e-03 1.015107e-03 3.322119e-03 277 4 + 36 18.000 -156.59642782 1.613798e-03 5.264024e-03 9.596567e-04 3.147590e-03 282 4 + 37 18.500 -156.59651774 2.093601e-03 6.837806e-03 9.076451e-04 2.982060e-03 287 4 + 38 19.000 -156.59659841 2.547391e-03 8.328818e-03 8.587443e-04 2.825105e-03 292 4 + 39 19.500 -156.59667079 2.976725e-03 9.741349e-03 8.126873e-04 2.676309e-03 297 4 + 40 20.000 -156.59673572 3.383025e-03 1.107948e-02 7.692501e-04 2.535272e-03 302 4 + 41 20.500 -156.59679398 3.767605e-03 1.234709e-02 7.282419e-04 2.401605e-03 307 4 + 42 21.000 -156.59684625 4.131680e-03 1.354787e-02 6.894962e-04 2.274937e-03 312 4 + 43 21.500 -156.59689315 4.476381e-03 1.468531e-02 6.528663e-04 2.154911e-03 317 4 + 44 22.000 -156.59693522 4.802768e-03 1.576274e-02 6.182207e-04 2.041185e-03 323 4 + 45 22.500 -156.59697297 5.111834e-03 1.678330e-02 5.854404e-04 1.933436e-03 329 4 + 46 23.000 -156.59700683 5.404511e-03 1.774999e-02 5.544169e-04 1.831353e-03 335 4 + 47 23.500 -156.59703721 5.681678e-03 1.866565e-02 5.250501e-04 1.734643e-03 341 4 + 48 24.000 -156.59706446 5.944164e-03 1.953294e-02 4.972474e-04 1.643025e-03 347 4 + 49 24.500 -156.59708891 6.192750e-03 2.035443e-02 4.709226e-04 1.556233e-03 353 4 + 50 25.000 -156.59711084 6.428177e-03 2.113253e-02 4.459951e-04 1.474016e-03 359 4 + 51 25.500 -156.59713052 6.651142e-03 2.186951e-02 4.223893e-04 1.396133e-03 365 4 + 52 26.000 -156.59714817 6.862305e-03 2.256756e-02 4.000342e-04 1.322357e-03 371 4 + 53 26.500 -156.59716401 7.062294e-03 2.322872e-02 3.788628e-04 1.252473e-03 377 4 + 54 27.000 -156.59717821 7.251698e-03 2.385494e-02 3.588122e-04 1.186276e-03 383 4 + 55 27.500 -156.59719095 7.431079e-03 2.444806e-02 3.398226e-04 1.123572e-03 389 4 + 56 28.000 -156.59720238 7.600967e-03 2.500983e-02 3.218378e-04 1.064178e-03 395 4 + 57 28.500 -156.59721264 7.761864e-03 2.554191e-02 3.048045e-04 1.007919e-03 401 4 + 58 29.000 -156.59722184 7.914246e-03 2.604586e-02 2.886725e-04 9.546313e-04 407 4 + 59 29.500 -156.59723009 8.058563e-03 2.652316e-02 2.733939e-04 9.041571e-04 413 4 + 60 30.000 -156.59723749 8.195242e-03 2.697523e-02 2.589236e-04 8.563487e-04 419 4 + 61 30.500 -156.59724413 8.324687e-03 2.740339e-02 2.452190e-04 8.110655e-04 425 4 + 62 31.000 -156.59725008 8.447281e-03 2.780891e-02 2.322394e-04 7.681745e-04 431 4 + 63 31.500 -156.59725543 8.563386e-03 2.819299e-02 2.199467e-04 7.275495e-04 437 4 + 64 32.000 -156.59726022 8.673346e-03 2.855676e-02 2.083044e-04 6.890711e-04 443 4 + 65 32.500 -156.59726452 8.777486e-03 2.890129e-02 1.972783e-04 6.526260e-04 449 4 + 66 33.000 -156.59726837 8.876113e-03 2.922759e-02 1.868356e-04 6.181070e-04 455 4 + 67 33.500 -156.59727183 8.969520e-03 2.953664e-02 1.769456e-04 5.854123e-04 461 4 + 68 34.000 -156.59727493 9.057982e-03 2.982934e-02 1.675790e-04 5.544458e-04 467 4 + 69 34.500 -156.59727772 9.141762e-03 3.010655e-02 1.587082e-04 5.251163e-04 473 4 + 70 35.000 -156.59728021 9.221107e-03 3.036911e-02 1.503069e-04 4.973372e-04 479 4 + 71 35.500 -156.59728245 9.296252e-03 3.061777e-02 1.423502e-04 4.710268e-04 485 4 + 72 36.000 -156.59728446 9.367419e-03 3.085328e-02 1.348148e-04 4.461074e-04 491 4 + 73 36.500 -156.59728626 9.434819e-03 3.107633e-02 1.276781e-04 4.225057e-04 497 4 + 74 37.000 -156.59728788 9.498651e-03 3.128757e-02 1.209193e-04 4.001520e-04 503 4 + 75 37.500 -156.59728933 9.559104e-03 3.148765e-02 1.145182e-04 3.789803e-04 509 4 + 76 38.000 -156.59729063 9.616356e-03 3.167713e-02 1.084560e-04 3.589284e-04 515 4 + 77 38.500 -156.59729179 9.670579e-03 3.185659e-02 1.027147e-04 3.399369e-04 521 4 + 78 39.000 -156.59729284 9.721930e-03 3.202656e-02 9.727725e-05 3.219498e-04 527 4 + 79 39.500 -156.59729378 9.770564e-03 3.218753e-02 9.212768e-05 3.049142e-04 533 4 + 80 40.000 -156.59729462 9.816623e-03 3.233998e-02 8.725071e-05 2.887796e-04 539 4 + 81 40.500 -156.59729537 9.860244e-03 3.248437e-02 8.263192e-05 2.734985e-04 545 4 + 82 41.000 -156.59729605 9.901555e-03 3.262112e-02 7.825763e-05 2.590257e-04 551 4 + 83 41.500 -156.59729666 9.940680e-03 3.275063e-02 7.411490e-05 2.453185e-04 557 4 + 84 42.000 -156.59729720 9.977734e-03 3.287329e-02 7.019148e-05 2.323365e-04 563 4 + 85 42.500 -156.59729769 1.001283e-02 3.298945e-02 6.647576e-05 2.200412e-04 569 4 + 86 43.000 -156.59729813 1.004606e-02 3.309947e-02 6.295674e-05 2.083965e-04 575 4 + 87 43.500 -156.59729852 1.007754e-02 3.320367e-02 5.962401e-05 1.973678e-04 581 4 + 88 44.000 -156.59729887 1.010734e-02 3.330235e-02 5.646770e-05 1.869227e-04 587 4 + 89 44.500 -156.59729919 1.013558e-02 3.339581e-02 5.347848e-05 1.770302e-04 593 4 + 90 45.000 -156.59729947 1.016231e-02 3.348432e-02 5.064751e-05 1.676611e-04 599 4 + 91 45.500 -156.59729973 1.018763e-02 3.356815e-02 4.796640e-05 1.587878e-04 605 4 + 92 46.000 -156.59729996 1.021161e-02 3.364754e-02 4.542722e-05 1.503840e-04 611 4 + 93 46.500 -156.59730016 1.023433e-02 3.372273e-02 4.302246e-05 1.424248e-04 617 4 + 94 47.000 -156.59730034 1.025583e-02 3.379394e-02 4.074500e-05 1.348869e-04 623 4 + 95 47.500 -156.59730051 1.027620e-02 3.386139e-02 3.858810e-05 1.277478e-04 629 4 + 96 48.000 -156.59730066 1.029550e-02 3.392526e-02 3.654539e-05 1.209865e-04 635 4 + 97 48.500 -156.59730079 1.031377e-02 3.398575e-02 3.461081e-05 1.145830e-04 641 4 + 98 49.000 -156.59730091 1.033107e-02 3.404304e-02 3.277864e-05 1.085184e-04 647 4 + 99 49.500 -156.59730101 1.034746e-02 3.409730e-02 3.104346e-05 1.027747e-04 653 4 + 100 50.000 -156.59730111 1.036298e-02 3.414869e-02 2.940014e-05 9.733501e-05 659 4 + 101 50.500 -156.59730120 1.037768e-02 3.419735e-02 2.784381e-05 9.218318e-05 665 4 + 102 51.000 -156.59730127 1.039160e-02 3.424344e-02 2.636986e-05 8.730400e-05 671 4 + 103 51.500 -156.59730134 1.040478e-02 3.428709e-02 2.497394e-05 8.268304e-05 677 4 + 104 52.000 -156.59730140 1.041727e-02 3.432844e-02 2.365192e-05 7.830664e-05 683 4 + 105 52.500 -156.59730146 1.042909e-02 3.436759e-02 2.239988e-05 7.416186e-05 689 4 + 106 53.000 -156.59730151 1.044029e-02 3.440467e-02 2.121412e-05 7.023644e-05 695 4 + 107 53.500 -156.59730155 1.045090e-02 3.443979e-02 2.009113e-05 6.651878e-05 701 4 + 108 54.000 -156.59730159 1.046094e-02 3.447304e-02 1.902759e-05 6.299788e-05 707 4 + 109 54.500 -156.59730163 1.047046e-02 3.450454e-02 1.802035e-05 5.966333e-05 713 4 + 110 55.000 -156.59730166 1.047946e-02 3.453437e-02 1.706643e-05 5.650527e-05 719 4 + 111 55.500 -156.59730169 1.048800e-02 3.456263e-02 1.616300e-05 5.351436e-05 725 4 + 112 56.000 -156.59730172 1.049608e-02 3.458938e-02 1.530740e-05 5.068175e-05 731 4 + 113 56.500 -156.59730174 1.050373e-02 3.461472e-02 1.449709e-05 4.799906e-05 737 4 + 114 57.000 -156.59730176 1.051098e-02 3.463872e-02 1.372967e-05 4.545837e-05 743 4 + 115 57.500 -156.59730178 1.051784e-02 3.466145e-02 1.300288e-05 4.305215e-05 749 4 + 116 58.000 -156.59730180 1.052434e-02 3.468298e-02 1.231457e-05 4.077330e-05 755 4 + 117 58.500 -156.59730181 1.053050e-02 3.470336e-02 1.166269e-05 3.861506e-05 761 4 + 118 59.000 -156.59730182 1.053633e-02 3.472267e-02 1.104531e-05 3.657106e-05 767 4 + 119 59.500 -156.59730184 1.054185e-02 3.474096e-02 1.046062e-05 3.463525e-05 773 4 + 120 60.000 -156.59730185 1.054708e-02 3.475827e-02 9.906884e-06 3.280190e-05 779 4 + 121 60.500 -156.59730186 1.055204e-02 3.477467e-02 9.382457e-06 3.106559e-05 785 4 + 122 61.000 -156.59730187 1.055673e-02 3.479021e-02 8.885791e-06 2.942119e-05 791 4 + 123 61.500 -156.59730187 1.056117e-02 3.480492e-02 8.415416e-06 2.786383e-05 797 4 + 124 62.000 -156.59730188 1.056538e-02 3.481885e-02 7.969942e-06 2.638890e-05 803 4 + 125 62.500 -156.59730189 1.056936e-02 3.483204e-02 7.548048e-06 2.499204e-05 809 4 + 126 63.000 -156.59730189 1.057313e-02 3.484454e-02 7.148488e-06 2.366912e-05 815 4 + 127 63.500 -156.59730190 1.057671e-02 3.485637e-02 6.770079e-06 2.241622e-05 821 4 + 128 64.000 -156.59730190 1.058009e-02 3.486758e-02 6.411701e-06 2.122965e-05 827 4 + 129 64.500 -156.59730191 1.058330e-02 3.487820e-02 6.072295e-06 2.010588e-05 833 4 + 130 65.000 -156.59730191 1.058633e-02 3.488825e-02 5.750855e-06 1.904159e-05 839 4 + 131 65.500 -156.59730191 1.058921e-02 3.489777e-02 5.446430e-06 1.803365e-05 845 4 + 132 66.000 -156.59730192 1.059193e-02 3.490679e-02 5.158121e-06 1.707905e-05 851 4 + 133 66.500 -156.59730192 1.059451e-02 3.491532e-02 4.885073e-06 1.617499e-05 857 4 + 134 67.000 -156.59730192 1.059695e-02 3.492341e-02 4.626480e-06 1.531877e-05 863 4 + 135 67.500 -156.59730192 1.059927e-02 3.493107e-02 4.381575e-06 1.450789e-05 869 4 + 136 68.000 -156.59730192 1.060146e-02 3.493833e-02 4.149634e-06 1.373992e-05 875 4 + 137 68.500 -156.59730193 1.060353e-02 3.494520e-02 3.929972e-06 1.301261e-05 881 4 + 138 69.000 -156.59730193 1.060550e-02 3.495170e-02 3.721937e-06 1.232379e-05 887 4 + 139 69.500 -156.59730193 1.060736e-02 3.495786e-02 3.524914e-06 1.167144e-05 893 4 + 140 70.000 -156.59730193 1.060912e-02 3.496370e-02 3.338322e-06 1.105362e-05 899 4 + 141 70.500 -156.59730193 1.061079e-02 3.496923e-02 3.161606e-06 1.046850e-05 905 4 + 142 71.000 -156.59730193 1.061237e-02 3.497446e-02 2.994245e-06 9.914351e-06 911 4 + 143 71.500 -156.59730193 1.061387e-02 3.497942e-02 2.835744e-06 9.389538e-06 917 4 + 144 72.000 -156.59730193 1.061528e-02 3.498411e-02 2.685632e-06 8.892506e-06 923 4 + 145 72.500 -156.59730194 1.061663e-02 3.498856e-02 2.543467e-06 8.421783e-06 929 4 + 146 73.000 -156.59730194 1.061790e-02 3.499277e-02 2.408828e-06 7.975978e-06 935 4 + 147 73.500 -156.59730194 1.061910e-02 3.499676e-02 2.281316e-06 7.553772e-06 941 4 + 148 74.000 -156.59730194 1.062024e-02 3.500053e-02 2.160553e-06 7.153914e-06 947 4 + 149 74.500 -156.59730194 1.062132e-02 3.500411e-02 2.046184e-06 6.775223e-06 953 4 + 150 75.000 -156.59730194 1.062235e-02 3.500750e-02 1.937868e-06 6.416577e-06 959 4 + 151 75.500 -156.59730194 1.062332e-02 3.501071e-02 1.835286e-06 6.076916e-06 965 4 + 152 76.000 -156.59730194 1.062423e-02 3.501374e-02 1.738135e-06 5.755235e-06 971 4 + 153 76.500 -156.59730194 1.062510e-02 3.501662e-02 1.646126e-06 5.450582e-06 977 4 + 154 77.000 -156.59730194 1.062592e-02 3.501935e-02 1.558988e-06 5.162056e-06 983 4 + 155 77.500 -156.59730194 1.062670e-02 3.502193e-02 1.476462e-06 4.888803e-06 989 4 + 156 78.000 -156.59730194 1.062744e-02 3.502437e-02 1.398305e-06 4.630014e-06 995 4 + 157 78.500 -156.59730194 1.062814e-02 3.502669e-02 1.324285e-06 4.384924e-06 1001 4 + 158 79.000 -156.59730194 1.062880e-02 3.502888e-02 1.254184e-06 4.152808e-06 1007 4 + 159 79.500 -156.59730194 1.062943e-02 3.503096e-02 1.187793e-06 3.932979e-06 1013 4 + 160 80.000 -156.59730194 1.063002e-02 3.503292e-02 1.124917e-06 3.724787e-06 1019 4 + 161 80.500 -156.59730194 1.063059e-02 3.503478e-02 1.065369e-06 3.527615e-06 1025 4 + 162 81.000 -156.59730194 1.063112e-02 3.503655e-02 1.008973e-06 3.340880e-06 1031 4 + 163 81.500 -156.59730194 1.063162e-02 3.503822e-02 9.555630e-07 3.164030e-06 1037 4 + 164 82.000 -156.59730194 1.063210e-02 3.503980e-02 9.049799e-07 2.996542e-06 1043 4 + 165 82.500 -156.59730194 1.063255e-02 3.504130e-02 8.570745e-07 2.837920e-06 1049 4 + 166 83.000 -156.59730194 1.063298e-02 3.504272e-02 8.117050e-07 2.687694e-06 1055 4 + 167 83.500 -156.59730194 1.063339e-02 3.504406e-02 7.687371e-07 2.545421e-06 1061 4 + 168 84.000 -156.59730194 1.063377e-02 3.504533e-02 7.280437e-07 2.410679e-06 1067 4 + 169 84.500 -156.59730194 1.063414e-02 3.504654e-02 6.895045e-07 2.283069e-06 1073 4 + 170 85.000 -156.59730194 1.063448e-02 3.504768e-02 6.530053e-07 2.162214e-06 1079 4 + 171 85.500 -156.59730194 1.063481e-02 3.504876e-02 6.184383e-07 2.047757e-06 1085 4 + 172 86.000 -156.59730194 1.063512e-02 3.504979e-02 5.857011e-07 1.939359e-06 1091 4 + 173 86.500 -156.59730194 1.063541e-02 3.505076e-02 5.546968e-07 1.836698e-06 1097 4 + 174 87.000 -156.59730194 1.063569e-02 3.505167e-02 5.253337e-07 1.739472e-06 1103 4 + 175 87.500 -156.59730194 1.063595e-02 3.505254e-02 4.975250e-07 1.647393e-06 1109 4 + 176 88.000 -156.59730194 1.063620e-02 3.505337e-02 4.711884e-07 1.560188e-06 1115 4 + 177 88.500 -156.59730194 1.063643e-02 3.505415e-02 4.462458e-07 1.477599e-06 1121 4 + 178 89.000 -156.59730194 1.063666e-02 3.505489e-02 4.226237e-07 1.399382e-06 1127 4 + 179 89.500 -156.59730194 1.063687e-02 3.505559e-02 4.002519e-07 1.325305e-06 1133 4 + 180 90.000 -156.59730194 1.063707e-02 3.505625e-02 3.790645e-07 1.255150e-06 1139 4 + 181 90.500 -156.59730194 1.063726e-02 3.505688e-02 3.589986e-07 1.188708e-06 1145 4 + 182 91.000 -156.59730194 1.063744e-02 3.505747e-02 3.399949e-07 1.125784e-06 1151 4 + 183 91.500 -156.59730194 1.063761e-02 3.505803e-02 3.219971e-07 1.066190e-06 1157 4 + 184 92.000 -156.59730194 1.063777e-02 3.505857e-02 3.049521e-07 1.009751e-06 1163 4 + 185 92.500 -156.59730194 1.063792e-02 3.505907e-02 2.888094e-07 9.562998e-07 1169 4 Starting VS-IMSRG flow step: s E0 ||Omega_1|| ||Omega_2|| ||Eta_1|| ||Eta_2|| Ncomm. nwritten - 0 93.000 -156.82112051 0.000000e+00 0.000000e+00 3.229158e-01 6.378179e+00 1175 5 - 1 93.500 -156.82112051 1.614579e-01 3.189090e+00 9.227151e-02 2.761175e+00 1185 5 - 2 94.000 -156.82112051 4.613575e-02 1.380588e+00 1.389716e-01 1.441046e+00 1193 6 - 3 94.500 -156.82112051 6.948582e-02 7.205232e-01 1.184896e-01 8.269515e-01 1201 7 - 4 95.000 -156.82112051 5.924480e-02 4.134758e-01 8.893113e-02 5.161785e-01 1208 8 - 5 95.500 -156.82112051 4.446556e-02 2.580892e-01 6.482493e-02 3.480325e-01 1215 9 - 6 96.000 -156.82112051 3.241246e-02 1.740163e-01 4.762813e-02 2.502454e-01 1221 10 - 7 96.500 -156.82112051 5.621124e-02 2.977482e-01 3.579587e-02 1.888736e-01 1229 10 - 8 97.000 -156.82112051 1.789794e-02 9.443679e-02 2.761892e-02 1.475082e-01 1235 11 - 9 97.500 -156.82112051 3.169891e-02 1.677967e-01 2.182843e-02 1.179198e-01 1241 11 - 10 98.000 -156.82112051 4.259421e-02 2.259722e-01 1.759313e-02 9.577239e-02 1247 11 - 11 98.500 -156.82112051 5.136303e-02 2.728047e-01 1.439169e-02 7.863840e-02 1254 11 - 12 99.000 -156.82112051 7.195846e-03 3.931920e-02 1.189991e-02 6.506699e-02 1259 12 - 13 99.500 -156.82112051 1.314460e-02 7.181908e-02 9.914422e-03 5.413659e-02 1265 12 - 14 100.000 -156.82112051 1.809915e-02 9.881703e-02 8.303855e-03 4.522856e-02 1271 12 - 15 100.500 -156.82112051 2.224708e-02 1.213313e-01 6.980349e-03 3.790600e-02 1277 12 - 16 101.000 -156.82112051 2.573216e-02 1.401637e-01 5.882715e-03 3.184809e-02 1283 12 - 17 101.500 -156.82112051 2.866758e-02 1.559548e-01 4.966578e-03 2.681195e-02 1289 12 - 18 102.000 -156.82112051 3.114431e-02 1.692216e-01 4.198543e-03 2.260932e-02 1295 12 - 19 102.500 -156.82112051 3.323662e-02 1.803857e-01 3.552696e-03 1.909163e-02 1301 12 - 20 103.000 -156.82112051 3.500574e-02 1.897928e-01 3.008432e-03 1.614001e-02 1307 12 - 21 103.500 -156.82112051 3.650263e-02 1.977283e-01 2.549069e-03 1.365835e-02 1313 12 - 22 104.000 -156.82112051 3.776984e-02 2.044287e-01 2.160925e-03 1.156828e-02 1319 12 - 23 104.500 -156.82112051 3.884307e-02 2.100906e-01 1.832674e-03 9.805459e-03 1325 12 - 24 105.000 -156.82112051 3.975235e-02 2.148784e-01 1.554884e-03 8.316818e-03 1331 12 - 25 105.500 -156.82112051 4.052298e-02 2.189293e-01 1.319662e-03 7.058360e-03 1337 12 - 26 106.000 -156.82112051 4.117629e-02 2.223586e-01 1.120388e-03 5.993494e-03 1343 12 - 27 106.500 -156.82112051 4.173029e-02 2.252629e-01 9.514930e-04 5.091694e-03 1349 12 - 28 107.000 -156.82112051 4.220019e-02 2.277237e-01 8.082886e-04 4.327428e-03 1356 12 - 29 107.500 -156.82112051 4.259886e-02 2.298095e-01 6.868218e-04 3.679298e-03 1363 12 - 30 108.000 -156.82112051 4.293716e-02 2.315779e-01 5.837569e-04 3.129333e-03 1370 12 - 31 108.500 -156.82112051 4.322431e-02 2.330777e-01 4.962771e-04 2.662418e-03 1377 12 - 32 109.000 -156.82112051 4.346809e-02 2.343501e-01 4.220026e-04 2.265824e-03 1384 12 - 33 109.500 -156.82112051 4.367510e-02 2.354298e-01 3.589212e-04 1.928813e-03 1391 12 - 34 110.000 -156.82112051 4.385091e-02 2.363463e-01 3.053310e-04 1.642322e-03 1398 12 - 35 110.500 -156.82112051 4.400025e-02 2.371243e-01 2.597919e-04 1.398690e-03 1405 12 - 36 111.000 -156.82112051 4.412714e-02 2.377850e-01 2.210844e-04 1.191437e-03 1412 12 - 37 111.500 -156.82112051 4.423496e-02 2.383462e-01 1.881759e-04 1.015080e-03 1419 12 - 38 112.000 -156.82112051 4.432660e-02 2.388229e-01 1.601912e-04 8.649718e-04 1426 12 - 39 112.500 -156.82112051 4.440450e-02 2.392279e-01 1.363885e-04 7.371735e-04 1433 12 - 40 113.000 -156.82112051 4.447072e-02 2.395721e-01 1.161388e-04 6.283445e-04 1440 12 - 41 113.500 -156.82112051 4.452704e-02 2.398646e-01 9.890843e-05 5.356499e-04 1447 12 - 42 114.000 -156.82112051 4.457493e-02 2.401133e-01 8.424462e-05 4.566827e-04 1454 12 - 43 114.500 -156.82112051 4.461566e-02 2.403247e-01 7.176295e-05 3.893980e-04 1461 12 - 44 115.000 -156.82112051 4.465031e-02 2.405045e-01 6.113702e-05 3.320584e-04 1468 12 - 45 115.500 -156.82112051 4.467979e-02 2.406573e-01 5.208959e-05 2.831867e-04 1475 12 - 46 116.000 -156.82112051 4.470487e-02 2.407874e-01 4.438509e-05 2.415268e-04 1482 12 - 47 116.500 -156.82112051 4.472622e-02 2.408980e-01 3.782334e-05 2.060101e-04 1489 12 - 48 117.000 -156.82112051 4.474438e-02 2.409921e-01 3.223416e-05 1.757273e-04 1496 12 - 49 117.500 -156.82112051 4.475984e-02 2.410721e-01 2.747285e-05 1.499045e-04 1503 12 - 50 118.000 -156.82112051 4.477300e-02 2.411402e-01 2.341636e-05 1.278827e-04 1510 12 - 51 118.500 -156.82112051 4.478421e-02 2.411982e-01 1.996000e-05 1.091008e-04 1517 12 - 52 119.000 -156.82112051 4.479374e-02 2.412476e-01 1.701471e-05 9.308090e-05 1524 12 - 53 119.500 -156.82112051 4.480186e-02 2.412895e-01 1.450471e-05 7.941587e-05 1531 12 - 54 120.000 -156.82112051 4.480878e-02 2.413253e-01 1.236549e-05 6.775881e-05 1538 12 - 55 120.500 -156.82112051 4.481466e-02 2.413557e-01 1.054214e-05 5.781407e-05 1545 12 - 56 121.000 -156.82112051 4.481968e-02 2.413816e-01 8.987909e-06 4.932969e-05 1552 12 - 57 121.500 -156.82112051 4.482395e-02 2.414037e-01 7.662995e-06 4.209089e-05 1559 12 - 58 122.000 -156.82112051 4.482758e-02 2.414224e-01 6.533495e-06 3.591454e-05 1566 12 - 59 122.500 -156.82112051 4.483068e-02 2.414384e-01 5.570534e-06 3.064452e-05 1573 12 - 60 123.000 -156.82112051 4.483331e-02 2.414520e-01 4.749516e-06 2.614769e-05 1580 12 - 61 123.500 -156.82112051 4.483556e-02 2.414636e-01 4.049488e-06 2.231051e-05 1587 12 - 62 124.000 -156.82112051 4.483747e-02 2.414734e-01 3.452595e-06 1.903615e-05 1594 12 - 63 124.500 -156.82112051 4.483910e-02 2.414818e-01 2.943626e-06 1.624200e-05 1601 12 - 64 125.000 -156.82112051 4.484048e-02 2.414890e-01 2.509617e-06 1.385762e-05 1608 12 - 65 125.500 -156.82112051 4.484166e-02 2.414951e-01 2.139518e-06 1.182289e-05 1615 12 - 66 126.000 -156.82112051 4.484267e-02 2.415002e-01 1.823913e-06 1.008654e-05 1622 12 - 67 126.500 -156.82112051 4.484352e-02 2.415046e-01 1.554774e-06 8.604821e-06 1629 12 - 68 127.000 -156.82112051 4.484425e-02 2.415084e-01 1.325260e-06 7.340408e-06 1636 12 - 69 127.500 -156.82112051 4.484487e-02 2.415116e-01 1.129537e-06 6.261450e-06 1643 12 - 70 128.000 -156.82112051 4.484540e-02 2.415143e-01 9.626320e-07 5.340772e-06 1650 12 - 71 128.500 -156.82112051 4.484585e-02 2.415166e-01 8.203041e-07 4.555181e-06 1657 12 - 72 129.000 -156.82112051 4.484623e-02 2.415185e-01 6.989383e-07 3.884890e-06 1664 12 - 73 129.500 -156.82112051 4.484655e-02 2.415202e-01 5.954516e-07 3.313011e-06 1671 12 - 74 130.000 -156.82112051 4.484683e-02 2.415216e-01 5.072152e-07 2.825132e-06 1678 12 - 75 130.500 -156.82112051 4.484706e-02 2.415228e-01 4.319873e-07 2.408955e-06 1685 12 - 76 131.000 -156.82112051 4.484726e-02 2.415238e-01 3.678562e-07 2.053984e-06 1692 12 - 77 131.500 -156.82112051 4.484743e-02 2.415247e-01 3.131915e-07 1.751261e-06 1699 12 - 78 132.000 -156.82112051 4.484757e-02 2.415254e-01 2.666027e-07 1.493143e-06 1706 12 - 79 132.500 -156.82112051 4.484770e-02 2.415260e-01 2.269041e-07 1.273107e-06 1713 12 - 80 133.000 -156.82112051 4.484780e-02 2.415266e-01 1.930844e-07 1.085586e-06 1720 12 - 81 133.500 -156.82112051 4.484789e-02 2.415270e-01 1.642813e-07 9.258295e-07 1727 12 + 0 92.500 -156.59730194 0.000000e+00 0.000000e+00 3.231086e-01 6.378496e+00 1169 5 + 1 93.000 -156.59730194 1.615543e-01 3.189248e+00 9.220145e-02 2.760810e+00 1179 5 + 2 93.500 -156.59730194 4.610072e-02 1.380405e+00 1.389228e-01 1.440777e+00 1187 6 + 3 94.000 -156.59730194 6.946140e-02 7.203887e-01 1.184664e-01 8.267959e-01 1195 7 + 4 94.500 -156.59730194 5.923319e-02 4.133980e-01 8.892307e-02 5.160971e-01 1202 8 + 5 95.000 -156.59730194 4.446154e-02 2.580486e-01 6.482492e-02 3.479942e-01 1209 9 + 6 95.500 -156.59730194 3.241246e-02 1.739971e-01 4.763196e-02 2.502307e-01 1215 10 + 7 96.000 -156.59730194 5.621315e-02 2.977220e-01 3.580123e-02 1.888713e-01 1223 10 + 8 96.500 -156.59730194 1.790061e-02 9.443566e-02 2.762460e-02 1.475122e-01 1229 11 + 9 97.000 -156.59730194 3.170442e-02 1.677977e-01 2.183386e-02 1.179267e-01 1235 11 + 10 97.500 -156.59730194 4.260242e-02 2.259770e-01 1.759806e-02 9.578036e-02 1241 11 + 11 98.000 -156.59730194 5.137366e-02 2.728140e-01 1.439603e-02 7.864655e-02 1248 11 + 12 98.500 -156.59730194 7.198015e-03 3.932328e-02 1.190366e-02 6.507487e-02 1253 12 + 13 99.000 -156.59730194 1.314864e-02 7.182711e-02 9.917645e-03 5.414399e-02 1259 12 + 14 99.500 -156.59730194 1.810480e-02 9.882878e-02 8.306606e-03 4.523539e-02 1265 12 + 15 100.000 -156.59730194 2.225410e-02 1.213465e-01 6.982693e-03 3.791223e-02 1271 12 + 16 100.500 -156.59730194 2.574034e-02 1.401821e-01 5.884712e-03 3.185375e-02 1277 12 + 17 101.000 -156.59730194 2.867675e-02 1.559759e-01 4.968281e-03 2.681706e-02 1283 12 + 18 101.500 -156.59730194 3.115432e-02 1.692453e-01 4.199999e-03 2.261392e-02 1289 12 + 19 102.000 -156.59730194 3.324734e-02 1.804117e-01 3.553943e-03 1.909575e-02 1295 12 + 20 102.500 -156.59730194 3.501708e-02 1.898209e-01 3.009502e-03 1.614370e-02 1301 12 + 21 103.000 -156.59730194 3.651449e-02 1.977582e-01 2.549990e-03 1.366165e-02 1307 12 + 22 103.500 -156.59730194 3.778214e-02 2.044601e-01 2.161718e-03 1.157121e-02 1313 12 + 23 104.000 -156.59730194 3.885576e-02 2.101235e-01 1.833359e-03 9.808073e-03 1319 12 + 24 104.500 -156.59730194 3.976537e-02 2.149125e-01 1.555476e-03 8.319141e-03 1325 12 + 25 105.000 -156.59730194 4.053629e-02 2.189646e-01 1.320174e-03 7.060420e-03 1331 12 + 26 105.500 -156.59730194 4.118984e-02 2.223948e-01 1.120832e-03 5.995318e-03 1337 12 + 27 106.000 -156.59730194 4.174405e-02 2.253001e-01 9.518776e-04 5.093307e-03 1343 12 + 28 106.500 -156.59730194 4.221413e-02 2.277616e-01 8.086222e-04 4.328853e-03 1350 12 + 29 107.000 -156.59730194 4.261296e-02 2.298480e-01 6.871112e-04 3.680555e-03 1357 12 + 30 107.500 -156.59730194 4.295140e-02 2.316170e-01 5.840080e-04 3.130440e-03 1364 12 + 31 108.000 -156.59730194 4.323867e-02 2.331173e-01 4.964951e-04 2.663392e-03 1371 12 + 32 108.500 -156.59730194 4.348255e-02 2.343902e-01 4.221918e-04 2.266679e-03 1378 12 + 33 109.000 -156.59730194 4.368964e-02 2.354703e-01 3.590853e-04 1.929564e-03 1385 12 + 34 109.500 -156.59730194 4.386553e-02 2.363871e-01 3.054734e-04 1.642980e-03 1392 12 + 35 110.000 -156.59730194 4.401494e-02 2.371655e-01 2.599154e-04 1.399266e-03 1399 12 + 36 110.500 -156.59730194 4.414188e-02 2.378264e-01 2.211915e-04 1.191942e-03 1406 12 + 37 111.000 -156.59730194 4.424975e-02 2.383878e-01 1.882687e-04 1.015521e-03 1413 12 + 38 111.500 -156.59730194 4.434143e-02 2.388647e-01 1.602716e-04 8.653565e-04 1420 12 + 39 112.000 -156.59730194 4.441937e-02 2.392699e-01 1.364581e-04 7.375090e-04 1427 12 + 40 112.500 -156.59730194 4.448563e-02 2.396142e-01 1.161990e-04 6.286368e-04 1434 12 + 41 113.000 -156.59730194 4.454197e-02 2.399068e-01 9.896052e-05 5.359043e-04 1441 12 + 42 113.500 -156.59730194 4.458989e-02 2.401556e-01 8.428964e-05 4.569038e-04 1448 12 + 43 114.000 -156.59730194 4.463064e-02 2.403671e-01 7.180184e-05 3.895901e-04 1455 12 + 44 114.500 -156.59730194 4.466531e-02 2.405470e-01 6.117058e-05 3.322251e-04 1462 12 + 45 115.000 -156.59730194 4.469480e-02 2.406999e-01 5.211853e-05 2.833312e-04 1469 12 + 46 115.500 -156.59730194 4.471990e-02 2.408300e-01 4.441003e-05 2.416519e-04 1476 12 + 47 116.000 -156.59730194 4.474125e-02 2.409407e-01 3.784480e-05 2.061183e-04 1483 12 + 48 116.500 -156.59730194 4.475943e-02 2.410348e-01 3.225261e-05 1.758208e-04 1490 12 + 49 117.000 -156.59730194 4.477489e-02 2.411149e-01 2.748870e-05 1.499852e-04 1497 12 + 50 117.500 -156.59730194 4.478806e-02 2.411831e-01 2.342995e-05 1.279522e-04 1504 12 + 51 118.000 -156.59730194 4.479927e-02 2.412411e-01 1.997164e-05 1.091606e-04 1511 12 + 52 118.500 -156.59730194 4.480881e-02 2.412904e-01 1.702467e-05 9.313235e-05 1518 12 + 53 119.000 -156.59730194 4.481694e-02 2.413325e-01 1.451321e-05 7.946004e-05 1525 12 + 54 119.500 -156.59730194 4.482385e-02 2.413682e-01 1.237274e-05 6.779665e-05 1532 12 + 55 120.000 -156.59730194 4.482974e-02 2.413987e-01 1.054830e-05 5.784644e-05 1539 12 + 56 120.500 -156.59730194 4.483476e-02 2.414246e-01 8.993140e-06 4.935734e-05 1546 12 + 57 121.000 -156.59730194 4.483903e-02 2.414466e-01 7.667421e-06 4.211444e-05 1553 12 + 58 121.500 -156.59730194 4.484267e-02 2.414654e-01 6.537229e-06 3.593457e-05 1560 12 + 59 122.000 -156.59730194 4.484576e-02 2.414814e-01 5.573673e-06 3.066150e-05 1567 12 + 60 122.500 -156.59730194 4.484840e-02 2.414950e-01 4.752145e-06 2.616206e-05 1574 12 + 61 123.000 -156.59730194 4.485065e-02 2.415066e-01 4.051679e-06 2.232263e-05 1581 12 + 62 123.500 -156.59730194 4.485256e-02 2.415164e-01 3.454412e-06 1.904633e-05 1588 12 + 63 124.000 -156.59730194 4.485419e-02 2.415248e-01 2.945122e-06 1.625053e-05 1595 12 + 64 124.500 -156.59730194 4.485557e-02 2.415320e-01 2.510840e-06 1.386473e-05 1602 12 + 65 125.000 -156.59730194 4.485675e-02 2.415380e-01 2.140508e-06 1.182880e-05 1609 12 + 66 125.500 -156.59730194 4.485776e-02 2.415432e-01 1.824706e-06 1.009142e-05 1616 12 + 67 126.000 -156.59730194 4.485862e-02 2.415476e-01 1.555401e-06 8.608835e-06 1623 12 + 68 126.500 -156.59730194 4.485934e-02 2.415514e-01 1.325746e-06 7.343688e-06 1630 12 + 69 127.000 -156.59730194 4.485996e-02 2.415546e-01 1.129905e-06 6.264115e-06 1637 12 + 70 127.500 -156.59730194 4.486049e-02 2.415573e-01 9.629018e-07 5.342922e-06 1644 12 + 71 128.000 -156.59730194 4.486094e-02 2.415596e-01 8.204924e-07 4.556907e-06 1651 12 + 72 128.500 -156.59730194 4.486132e-02 2.415615e-01 6.990596e-07 3.886268e-06 1658 12 + 73 129.000 -156.59730194 4.486164e-02 2.415632e-01 5.955187e-07 3.314108e-06 1665 12 + 74 129.500 -156.59730194 4.486192e-02 2.415646e-01 5.072389e-07 2.826006e-06 1672 12 + 75 130.000 -156.59730194 4.486215e-02 2.415658e-01 4.319772e-07 2.409658e-06 1679 12 + 76 130.500 -156.59730194 4.486235e-02 2.415668e-01 3.678205e-07 2.054560e-06 1686 12 + 77 131.000 -156.59730194 4.486252e-02 2.415677e-01 3.131373e-07 1.751751e-06 1693 12 + 78 131.500 -156.59730194 4.486266e-02 2.415684e-01 2.665364e-07 1.493582e-06 1700 12 + 79 132.000 -156.59730194 4.486278e-02 2.415690e-01 2.268312e-07 1.273526e-06 1707 12 + 80 132.500 -156.59730194 4.486289e-02 2.415695e-01 1.930099e-07 1.086012e-06 1714 12 + 81 133.000 -156.59730194 4.486297e-02 2.415700e-01 1.642094e-07 9.262894e-07 1721 12 Operator:Rp2 - HF point proton radius 2.049501 charge radius 2.211009 => HF+PT 2.229791 -IMSRG point proton radius 2.099186 charge radius 2.257141 + HF point proton radius 2.049643 charge radius 2.211140 => HF+PT 2.229936 +IMSRG point proton radius 2.099363 charge radius 2.257306 Mg24 Z,N=(12,12) c(8,8) v(4,4) mdim: 28503 ( 4.45 ) J [0.0, 2.0, 2.0, 3.0, 4.0, 4.0, 5.0, 1.0, 6.0, 2.0] -En. -109.0654 -107.7170 -106.0557 -104.8003 -104.7485 -103.2378 -101.1769 -100.7040 -100.1376 -100.1193 -Ex. 0.0000 1.3484 3.0097 4.2651 4.3170 5.8276 7.8886 8.3615 8.9278 8.9461 +En. -108.9515 -107.6036 -105.9415 -104.6864 -104.6365 -103.1245 -101.0649 -100.5969 -100.0279 -100.0124 +Ex. 0.0000 1.3479 3.0101 4.2651 4.3150 5.8270 7.8866 8.3546 8.9236 8.9391