diff --git a/Be2_cc-pvtz/fit_eq.sh b/Be2_cc-pvtz/fit_eq.sh index 3493cf1..490d30c 100755 --- a/Be2_cc-pvtz/fit_eq.sh +++ b/Be2_cc-pvtz/fit_eq.sh @@ -156,14 +156,14 @@ paste det.dat a.dat D.dat > det_aD.dat grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat grep 'FCI' det_aD.dat | sed 's/FCI/1.0E9 /g' >> det_aD_FCI.dat -grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat -grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat -grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' >> det_aD_CI.dat +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat +#grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat +grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' > det_aD_CI.dat grep 'hfCISDT ' det_aD.dat | sed 's/hfCISDT /3 /g' >> det_aD_CI.dat grep 'hfCISDTQ' det_aD.dat | sed 's/hfCISDTQ/4 /g' >> det_aD_CI.dat -grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat -grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' >> det_aD_ooCI.dat +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat +grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' > det_aD_ooCI.dat grep 'ooCISD ' det_aD.dat | sed 's/ooCISD /2 /g' >> det_aD_ooCI.dat grep 'ooCISDT ' det_aD.dat | sed 's/ooCISDT /3 /g' >> det_aD_ooCI.dat @@ -174,8 +174,8 @@ grep 'hfCIs4 ' det_aD.dat | sed 's/hfCIs4/4 /g' >> det_aD_CIs.dat grep 'ooCIs0 ' det_aD.dat | sed 's/ooCIs0/0 /g' > det_aD_ooCIs.dat grep 'ooCIs2 ' det_aD.dat | sed 's/ooCIs2/2 /g' >> det_aD_ooCIs.dat -grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat -grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' >> det_aD_CIo.dat +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat +grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' > det_aD_CIo.dat grep 'hfCIo1.5' det_aD.dat | sed 's/hfCIo1.5/1.5 /g' >> det_aD_CIo.dat grep 'hfCIo2 ' det_aD.dat | sed 's/hfCIo2/2 /g' >> det_aD_CIo.dat grep 'hfCIo2.5' det_aD.dat | sed 's/hfCIo2.5/2.5 /g' >> det_aD_CIo.dat @@ -183,8 +183,8 @@ grep 'hfCIo3 ' det_aD.dat | sed 's/hfCIo3/3 /g' >> det_aD_CIo.dat grep 'hfCIo3.5' det_aD.dat | sed 's/hfCIo3.5/3.5 /g' >> det_aD_CIo.dat grep 'hfCIo4 ' det_aD.dat | sed 's/hfCIo4/4 /g' >> det_aD_CIo.dat -grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat -grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' >> det_aD_ooCIo.dat +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat +grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' > det_aD_ooCIo.dat grep 'ooCIo1.5' det_aD.dat | sed 's/ooCIo1.5/1.5 /g' >> det_aD_ooCIo.dat grep 'ooCIo2 ' det_aD.dat | sed 's/ooCIo2/2 /g' >> det_aD_ooCIo.dat grep 'ooCIo2.5' det_aD.dat | sed 's/ooCIo2.5/2.5 /g' >> det_aD_ooCIo.dat @@ -255,14 +255,14 @@ paste det.dat xe.dat D.dat > det_xe.dat grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat grep 'FCI' det_xe.dat | sed 's/FCI/1E9 /g' >> det_xe_FCI.dat -grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat -grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat -grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' >> det_xe_CI.dat +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat +#grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat +grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' > det_xe_CI.dat grep 'hfCISDT ' det_xe.dat | sed 's/hfCISDT /3 /g' >> det_xe_CI.dat grep 'hfCISDTQ' det_xe.dat | sed 's/hfCISDTQ/4 /g' >> det_xe_CI.dat -grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat -grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' >> det_xe_ooCI.dat +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat +grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' > det_xe_ooCI.dat grep 'ooCISD ' det_xe.dat | sed 's/ooCISD /2 /g' >> det_xe_ooCI.dat grep 'ooCISDT ' det_xe.dat | sed 's/ooCISDT /3 /g' >> det_xe_ooCI.dat @@ -273,8 +273,8 @@ grep 'hfCIs4 ' det_xe.dat | sed 's/hfCIs4/4 /g' >> det_xe_CIs.dat grep 'ooCIs0 ' det_xe.dat | sed 's/ooCIs0/0 /g' > det_xe_ooCIs.dat grep 'ooCIs2 ' det_xe.dat | sed 's/ooCIs2/2 /g' >> det_xe_ooCIs.dat -grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat -grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' >> det_xe_CIo.dat +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat +grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' > det_xe_CIo.dat grep 'hfCIo1.5' det_xe.dat | sed 's/hfCIo1.5/1.5 /g' >> det_xe_CIo.dat grep 'hfCIo2 ' det_xe.dat | sed 's/hfCIo2/2 /g' >> det_xe_CIo.dat grep 'hfCIo2.5' det_xe.dat | sed 's/hfCIo2.5/2.5 /g' >> det_xe_CIo.dat @@ -282,8 +282,8 @@ grep 'hfCIo3 ' det_xe.dat | sed 's/hfCIo3/3 /g' >> det_xe_CIo.dat grep 'hfCIo3.5' det_xe.dat | sed 's/hfCIo3.5/3.5 /g' >> det_xe_CIo.dat grep 'hfCIo4 ' det_xe.dat | sed 's/hfCIo4/4 /g' >> det_xe_CIo.dat -grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat -grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' >> det_xe_ooCIo.dat +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat +grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' > det_xe_ooCIo.dat grep 'ooCIo1.5' det_xe.dat | sed 's/ooCIo1.5/1.5 /g' >> det_xe_ooCIo.dat grep 'ooCIo2 ' det_xe.dat | sed 's/ooCIo2/2 /g' >> det_xe_ooCIo.dat grep 'ooCIo2.5' det_xe.dat | sed 's/ooCIo2.5/2.5 /g' >> det_xe_ooCIo.dat diff --git a/Be2_cc-pvtz/fit_eq_0.sh b/Be2_cc-pvtz/fit_eq_0.sh new file mode 100755 index 0000000..3f3e7b1 --- /dev/null +++ b/Be2_cc-pvtz/fit_eq_0.sh @@ -0,0 +1,292 @@ +#!/bin/bash + +####################################################################### + +if [ -f "fit.log" ]; then + rm "fit.log" +fi + +gnuplot fit_eq.gnu + +####################################################################### +grep ^a fit.log | grep = | sed 's/=//g' | cut -c -30 | + sed 's/a0 /FCI /g' | + sed 's/a1 /RHF /g' | + sed 's/a2 /hfCISD /g' | + sed 's/a3 /hfCISDT /g' | + sed 's/a4 /hfCISDTQ /g' | + sed 's/a5 /hfCIs0 /g' | + sed 's/a6 /hfCIs2 /g' | + sed 's/a7 /hfCIs4 /g' | + sed 's/a8 /hfCIo1 /g' | + sed 's/a9 /hfCIo1.5 /g' | + sed 's/a10/hfCIo2 /g' | + sed 's/a11/hfCIo2.5 /g' | + sed 's/a12/hfCIo3 /g' | + sed 's/a13/hfCIo3.5 /g' | + sed 's/a14/hfCIo4 /g' | + sed 's/a15/ooCISD /g' | + sed 's/a16/ooCISDT /g' | + sed 's/a17/ooCIs0 /g' | + sed 's/a18/ooCIs2 /g' | + sed 's/a19/ooCIo1 /g' | + sed 's/a20/ooCIo1.5 /g' | + sed 's/a21/ooCIo2 /g' | + sed 's/a22/ooCIo2.5 /g' | + sed 's/a23/ooCIo3 /g' > a.dat + +grep ^D fit.log | grep = | sed 's/=//g' | cut -c -30 | + sed 's/D0 /FCI /g' | + sed 's/D1 /RHF /g' | + sed 's/D2 /hfCISD /g' | + sed 's/D3 /hfCISDT /g' | + sed 's/D4 /hfCISDTQ /g' | + sed 's/D5 /hfCIs0 /g' | + sed 's/D6 /hfCIs2 /g' | + sed 's/D7 /hfCIs4 /g' | + sed 's/D8 /hfCIo1 /g' | + sed 's/D9 /hfCIo1.5 /g' | + sed 's/D10/hfCIo2 /g' | + sed 's/D11/hfCIo2.5 /g' | + sed 's/D12/hfCIo3 /g' | + sed 's/D13/hfCIo3.5 /g' | + sed 's/D14/hfCIo4 /g' | + sed 's/D15/ooCISD /g' | + sed 's/D16/ooCISDT /g' | + sed 's/D17/ooCIs0 /g' | + sed 's/D18/ooCIs2 /g' | + sed 's/D19/ooCIo1 /g' | + sed 's/D20/ooCIo1.5 /g' | + sed 's/D21/ooCIo2 /g' | + sed 's/D22/ooCIo2.5 /g' | + sed 's/D23/ooCIo3 /g' > D.dat + +paste a.dat D.dat > aD.dat + +grep 'FCI' aD.dat | sed 's/FCI/-0.1 /g' > aD_FCI.dat +grep 'FCI' aD.dat | sed 's/FCI/4.1 /g' >> aD_FCI.dat + +grep ^a fit.log | grep = | sed 's/=//g' | cut -c -30 > a0.dat +grep 'a0 ' a0.dat | sed 's/a0 /FCI /g' > a.dat +grep 'a1 ' a0.dat | sed 's/a1 /RHF /g' >> a.dat +grep 'a1 ' a0.dat | sed 's/a1 /hfCIS /g' >> a.dat +grep 'a2 ' a0.dat | sed 's/a2 /hfCISD /g' >> a.dat +grep 'a3 ' a0.dat | sed 's/a3 /hfCISDT /g' >> a.dat +grep 'a4 ' a0.dat | sed 's/a4 /hfCISDTQ /g' >> a.dat +grep 'a5 ' a0.dat | sed 's/a5 /hfCIs0 /g' >> a.dat +grep 'a6 ' a0.dat | sed 's/a6 /hfCIs2 /g' >> a.dat +grep 'a7 ' a0.dat | sed 's/a7 /hfCIs4 /g' >> a.dat +grep 'a8 ' a0.dat | sed 's/a8 /hfCIo1 /g' >> a.dat +grep 'a9 ' a0.dat | sed 's/a9 /hfCIo1.5 /g' >> a.dat +grep 'a10 ' a0.dat | sed 's/a10/hfCIo2 /g' >> a.dat +grep 'a11 ' a0.dat | sed 's/a11/hfCIo2.5 /g' >> a.dat +grep 'a12 ' a0.dat | sed 's/a12/hfCIo3 /g' >> a.dat +grep 'a13 ' a0.dat | sed 's/a13/hfCIo3.5 /g' >> a.dat +grep 'a14 ' a0.dat | sed 's/a14/hfCIo4 /g' >> a.dat +grep 'a1 ' a0.dat | sed 's/a1 /ooCIS /g' >> a.dat +grep 'a15 ' a0.dat | sed 's/a15/ooCISD /g' >> a.dat +grep 'a16 ' a0.dat | sed 's/a16/ooCISDT /g' >> a.dat +grep 'a17 ' a0.dat | sed 's/a17/ooCIs0 /g' >> a.dat +grep 'a18 ' a0.dat | sed 's/a18/ooCIs2 /g' >> a.dat +grep 'a19 ' a0.dat | sed 's/a19/ooCIo1 /g' >> a.dat +grep 'a20 ' a0.dat | sed 's/a20/ooCIo1.5 /g' >> a.dat +grep 'a21 ' a0.dat | sed 's/a21/ooCIo2 /g' >> a.dat +grep 'a22 ' a0.dat | sed 's/a22/ooCIo2.5 /g' >> a.dat +grep 'a23 ' a0.dat | sed 's/a23/ooCIo3 /g' >> a.dat + +grep ^D fit.log | grep = | sed 's/=//g' | cut -c -30 > D0.dat +grep 'D0 ' D0.dat | sed 's/D0 /FCI /g' > D.dat +grep 'D1 ' D0.dat | sed 's/D1 /RHF /g' >> D.dat +grep 'D1 ' D0.dat | sed 's/D1 /hfCIS /g' >> D.dat +grep 'D2 ' D0.dat | sed 's/D2 /hfCISD /g' >> D.dat +grep 'D3 ' D0.dat | sed 's/D3 /hfCISDT /g' >> D.dat +grep 'D4 ' D0.dat | sed 's/D4 /hfCISDTQ /g' >> D.dat +grep 'D5 ' D0.dat | sed 's/D5 /hfCIs0 /g' >> D.dat +grep 'D6 ' D0.dat | sed 's/D6 /hfCIs2 /g' >> D.dat +grep 'D7 ' D0.dat | sed 's/D7 /hfCIs4 /g' >> D.dat +grep 'D8 ' D0.dat | sed 's/D8 /hfCIo1 /g' >> D.dat +grep 'D9 ' D0.dat | sed 's/D9 /hfCIo1.5 /g' >> D.dat +grep 'D10 ' D0.dat | sed 's/D10/hfCIo2 /g' >> D.dat +grep 'D11 ' D0.dat | sed 's/D11/hfCIo2.5 /g' >> D.dat +grep 'D12 ' D0.dat | sed 's/D12/hfCIo3 /g' >> D.dat +grep 'D13 ' D0.dat | sed 's/D13/hfCIo3.5 /g' >> D.dat +grep 'D14 ' D0.dat | sed 's/D14/hfCIo4 /g' >> D.dat +grep 'D1 ' D0.dat | sed 's/D1 /ooCIS /g' >> D.dat +grep 'D15 ' D0.dat | sed 's/D15/ooCISD /g' >> D.dat +grep 'D16 ' D0.dat | sed 's/D16/ooCISDT /g' >> D.dat +grep 'D17 ' D0.dat | sed 's/D17/ooCIs0 /g' >> D.dat +grep 'D18 ' D0.dat | sed 's/D18/ooCIs2 /g' >> D.dat +grep 'D19 ' D0.dat | sed 's/D19/ooCIo1 /g' >> D.dat +grep 'D20 ' D0.dat | sed 's/D20/ooCIo1.5 /g' >> D.dat +grep 'D21 ' D0.dat | sed 's/D21/ooCIo2 /g' >> D.dat +grep 'D22 ' D0.dat | sed 's/D22/ooCIo2.5 /g' >> D.dat +grep 'D23 ' D0.dat | sed 's/D23/ooCIo3 /g' >> D.dat + +if [[ $(grep 'a0 ' a0.dat) ]]; then grep '0 ' det_FCI.dat | cut -c 4- > det.dat; fi +if [[ $(grep 'a1 ' a0.dat) ]]; then grep '0 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a1 ' a0.dat) ]]; then grep '1 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a2 ' a0.dat) ]]; then grep '2 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a3 ' a0.dat) ]]; then grep '3 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a4 ' a0.dat) ]]; then grep '4 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a5 ' a0.dat) ]]; then grep '0 ' det_CIs.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a6 ' a0.dat) ]]; then grep '2 ' det_CIs.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a7 ' a0.dat) ]]; then grep '4 ' det_CIs.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a8 ' a0.dat) ]]; then grep '1 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a9 ' a0.dat) ]]; then grep '1.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a10 ' a0.dat) ]]; then grep '2 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a11 ' a0.dat) ]]; then grep '2.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a12 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a13 ' a0.dat) ]]; then grep '3.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a14 ' a0.dat) ]]; then grep '4 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a1 ' a0.dat) ]]; then grep '1 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a15 ' a0.dat) ]]; then grep '2 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a16 ' a0.dat) ]]; then grep '3 ' det_CI.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a17 ' a0.dat) ]]; then grep '0 ' det_CIs.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a18 ' a0.dat) ]]; then grep '2 ' det_CIs.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a19 ' a0.dat) ]]; then grep '1 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a20 ' a0.dat) ]]; then grep '1.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a21 ' a0.dat) ]]; then grep '2 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a22 ' a0.dat) ]]; then grep '2.5 ' det_CIo.dat | cut -c 4- >> det.dat; fi +if [[ $(grep 'a23 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det.dat; fi + +paste det.dat a.dat D.dat > det_aD.dat + +grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat +grep 'FCI' det_aD.dat | sed 's/FCI/1.0E8 /g' >> det_aD_FCI.dat + +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat +#grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat +grep 'hfCISD ' det_aD.dat | sed 's/hfCISD /2 /g' > det_aD_CI.dat +grep 'hfCISDT ' det_aD.dat | sed 's/hfCISDT /3 /g' >> det_aD_CI.dat +grep 'hfCISDTQ' det_aD.dat | sed 's/hfCISDTQ/4 /g' >> det_aD_CI.dat + +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCI.dat +#grep 'ooCIS ' det_aD.dat | sed 's/ooCIS /1 /g' >> det_aD_ooCI.dat +grep 'ooCISD ' det_aD.dat | sed 's/ooCISD /2 /g' > det_aD_ooCI.dat +grep 'ooCISDT' det_aD.dat | sed 's/ooCISDT/3 /g' >> det_aD_ooCI.dat + +grep 'hfCIs0 ' det_aD.dat | sed 's/hfCIs0/0 /g' > det_aD_CIs.dat +grep 'hfCIs2 ' det_aD.dat | sed 's/hfCIs2/2 /g' >> det_aD_CIs.dat +grep 'hfCIs4 ' det_aD.dat | sed 's/hfCIs4/4 /g' >> det_aD_CIs.dat + +grep 'ooCIs0 ' det_aD.dat | sed 's/ooCIs0/0 /g' > det_aD_ooCIs.dat +grep 'ooCIs2 ' det_aD.dat | sed 's/ooCIs2/2 /g' >> det_aD_ooCIs.dat + +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CIo.dat +grep 'hfCIo1 ' det_aD.dat | sed 's/hfCIo1/1 /g' > det_aD_CIo.dat +grep 'hfCIo1.5' det_aD.dat | sed 's/hfCIo1.5/1.5 /g' >> det_aD_CIo.dat +grep 'hfCIo2 ' det_aD.dat | sed 's/hfCIo2/2 /g' >> det_aD_CIo.dat +grep 'hfCIo2.5' det_aD.dat | sed 's/hfCIo2.5/2.5 /g' >> det_aD_CIo.dat +grep 'hfCIo3 ' det_aD.dat | sed 's/hfCIo3/3 /g' >> det_aD_CIo.dat +grep 'hfCIo3.5' det_aD.dat | sed 's/hfCIo3.5/3.5 /g' >> det_aD_CIo.dat +grep 'hfCIo4 ' det_aD.dat | sed 's/hfCIo4/4 /g' >> det_aD_CIo.dat + +#grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_ooCIo.dat +grep 'ooCIo1 ' det_aD.dat | sed 's/ooCIo1/1 /g' > det_aD_ooCIo.dat +grep 'ooCIo1.5' det_aD.dat | sed 's/ooCIo1.5/1.5 /g' >> det_aD_ooCIo.dat +grep 'ooCIo2 ' det_aD.dat | sed 's/ooCIo2/2 /g' >> det_aD_ooCIo.dat +grep 'ooCIo2.5' det_aD.dat | sed 's/ooCIo2.5/2.5 /g' >> det_aD_ooCIo.dat +grep 'ooCIo3 ' det_aD.dat | sed 's/ooCIo3/3 /g' >> det_aD_ooCIo.dat + +####################################################################### + +####################################################################### +grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 > xe0.dat + +grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 | + sed 's/xe0 /FCI /g' | + sed 's/xe1 /RHF /g' | + sed 's/xe2 /hfCISD /g' | + sed 's/xe3 /hfCISDT /g' | + sed 's/xe4 /hfCISDTQ /g' | + sed 's/xe5 /hfCIs0 /g' | + sed 's/xe6 /hfCIs2 /g' | + sed 's/xe7 /hfCIs4 /g' | + sed 's/xe8 /hfCIo1 /g' | + sed 's/xe9 /hfCIo1.5 /g' | + sed 's/xe10/hfCIo2 /g' | + sed 's/xe11/hfCIo2.5 /g' | + sed 's/xe12/hfCIo3 /g' | + sed 's/xe13/hfCIo3.5 /g' | + sed 's/xe14/hfCIo4 /g' | + sed 's/xe15/ooCISD /g' | + sed 's/xe16/ooCISDT /g' | + sed 's/xe17/ooCIs0 /g' | + sed 's/xe18/ooCIs2 /g' | + sed 's/xe19/ooCIo1 /g' | + sed 's/xe20/ooCIo1.5 /g' | + sed 's/xe21/ooCIo2 /g' | + sed 's/xe22/ooCIo2.5 /g' | + sed 's/xe23/ooCIo3 /g' > xe.dat + +grep ^xe fit.log | grep = | sed 's/=//g' | cut -c -30 > xe0.dat +grep 'xe0 ' xe0.dat | sed 's/xe0 /FCI /g' > xe.dat +grep 'xe1 ' xe0.dat | sed 's/xe1 /RHF /g' >> xe.dat +grep 'xe1 ' xe0.dat | sed 's/xe1 /hfCIS /g' >> xe.dat +grep 'xe2 ' xe0.dat | sed 's/xe2 /hfCISD /g' >> xe.dat +grep 'xe3 ' xe0.dat | sed 's/xe3 /hfCISDT /g' >> xe.dat +grep 'xe4 ' xe0.dat | sed 's/xe4 /hfCISDTQ /g' >> xe.dat +grep 'xe5 ' xe0.dat | sed 's/xe5 /hfCIs0 /g' >> xe.dat +grep 'xe6 ' xe0.dat | sed 's/xe6 /hfCIs2 /g' >> xe.dat +grep 'xe7 ' xe0.dat | sed 's/xe7 /hfCIs4 /g' >> xe.dat +grep 'xe8 ' xe0.dat | sed 's/xe8 /hfCIo1 /g' >> xe.dat +grep 'xe9 ' xe0.dat | sed 's/xe9 /hfCIo1.5 /g' >> xe.dat +grep 'xe10 ' xe0.dat | sed 's/xe10/hfCIo2 /g' >> xe.dat +grep 'xe11 ' xe0.dat | sed 's/xe11/hfCIo2.5 /g' >> xe.dat +grep 'xe12 ' xe0.dat | sed 's/xe12/hfCIo3 /g' >> xe.dat +grep 'xe13 ' xe0.dat | sed 's/xe13/hfCIo3.5 /g' >> xe.dat +grep 'xe14 ' xe0.dat | sed 's/xe14/hfCIo4 /g' >> xe.dat +grep 'xe1 ' xe0.dat | sed 's/xe1 /ooCIS /g' >> xe.dat +grep 'xe15 ' xe0.dat | sed 's/xe15/ooCISD /g' >> xe.dat +grep 'xe16 ' xe0.dat | sed 's/xe16/ooCISDT /g' >> xe.dat +grep 'xe17 ' xe0.dat | sed 's/xe17/ooCIs0 /g' >> xe.dat +grep 'xe18 ' xe0.dat | sed 's/xe18/ooCIs2 /g' >> xe.dat +grep 'xe19 ' xe0.dat | sed 's/xe19/ooCIo1 /g' >> xe.dat +grep 'xe20 ' xe0.dat | sed 's/xe20/ooCIo1.5 /g' >> xe.dat +grep 'xe21 ' xe0.dat | sed 's/xe21/ooCIo2 /g' >> xe.dat +grep 'xe22 ' xe0.dat | sed 's/xe22/ooCIo2.5 /g' >> xe.dat +grep 'xe23 ' xe0.dat | sed 's/xe23/ooCIo3 /g' >> xe.dat + +paste det.dat xe.dat D.dat > det_xe.dat + +grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat +grep 'FCI' det_xe.dat | sed 's/FCI/1E8 /g' >> det_xe_FCI.dat + +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat +#grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat +grep 'hfCISD ' det_xe.dat | sed 's/hfCISD /2 /g' > det_xe_CI.dat +grep 'hfCISDT ' det_xe.dat | sed 's/hfCISDT /3 /g' >> det_xe_CI.dat +grep 'hfCISDTQ' det_xe.dat | sed 's/hfCISDTQ/4 /g' >> det_xe_CI.dat + +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCI.dat +#grep 'ooCIS ' det_xe.dat | sed 's/ooCIS /1 /g' >> det_xe_ooCI.dat +grep 'ooCISD ' det_xe.dat | sed 's/ooCISD /2 /g' > det_xe_ooCI.dat +grep 'ooCISDT' det_xe.dat | sed 's/ooCISDT/3 /g' >> det_xe_ooCI.dat + +grep 'hfCIs0 ' det_xe.dat | sed 's/hfCIs0/0 /g' > det_xe_CIs.dat +grep 'hfCIs2 ' det_xe.dat | sed 's/hfCIs2/2 /g' >> det_xe_CIs.dat +grep 'hfCIs4 ' det_xe.dat | sed 's/hfCIs4/4 /g' >> det_xe_CIs.dat + +grep 'ooCIs0 ' det_xe.dat | sed 's/ooCIs0/0 /g' > det_xe_ooCIs.dat +grep 'ooCIs2 ' det_xe.dat | sed 's/ooCIs2/2 /g' >> det_xe_ooCIs.dat + +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CIo.dat +grep 'hfCIo1 ' det_xe.dat | sed 's/hfCIo1/1 /g' > det_xe_CIo.dat +grep 'hfCIo1.5' det_xe.dat | sed 's/hfCIo1.5/1.5 /g' >> det_xe_CIo.dat +grep 'hfCIo2 ' det_xe.dat | sed 's/hfCIo2/2 /g' >> det_xe_CIo.dat +grep 'hfCIo2.5' det_xe.dat | sed 's/hfCIo2.5/2.5 /g' >> det_xe_CIo.dat +grep 'hfCIo3 ' det_xe.dat | sed 's/hfCIo3/3 /g' >> det_xe_CIo.dat +grep 'hfCIo3.5' det_xe.dat | sed 's/hfCIo3.5/3.5 /g' >> det_xe_CIo.dat +grep 'hfCIo4 ' det_xe.dat | sed 's/hfCIo4/4 /g' >> det_xe_CIo.dat + +#grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_ooCIo.dat +grep 'ooCIo1 ' det_xe.dat | sed 's/ooCIo1/1 /g' > det_xe_ooCIo.dat +grep 'ooCIo1.5' det_xe.dat | sed 's/ooCIo1.5/1.5 /g' >> det_xe_ooCIo.dat +grep 'ooCIo2 ' det_xe.dat | sed 's/ooCIo2/2 /g' >> det_xe_ooCIo.dat +grep 'ooCIo2.5' det_xe.dat | sed 's/ooCIo2.5/2.5 /g' >> det_xe_ooCIo.dat +grep 'ooCIo3 ' det_xe.dat | sed 's/ooCIo3/3 /g' >> det_xe_ooCIo.dat + +####################################################################### + +#epspdf fit_eq.eps +#okular fit_eq.pdf diff --git a/H4_cc-pvdz/pes_ooCIs2.dat b/H4_cc-pvdz/pes_ooCIs2.dat index ae43735..8e6871a 100644 --- a/H4_cc-pvdz/pes_ooCIs2.dat +++ b/H4_cc-pvdz/pes_ooCIs2.dat @@ -31,3 +31,43 @@ 3.4 -2.06260088 3.5 -2.05451700 3.6 -2.04724897 +3.7 -2.04075589 +3.8 -2.03498844 +3.9 -2.02989402 +4.0 -2.02541500 +4.1 -2.02149231 +4.2 -2.01807018 +4.3 -2.01509383 +4.4 -2.01251351 +4.5 -2.01028245 +4.6 -2.00835935 +4.7 -2.00670418 +4.8 -2.00528242 +4.9 -2.00406073 +5.0 -2.00301776 +5.1 -2.00213173 +5.2 -2.00137685 +5.3 -2.00073304 +5.4 -2.00018601 +5.5 -1.99972098 +5.6 -1.99932659 +5.7 -1.99899193 +5.8 -1.99870814 +5.9 -1.99846786 +6.0 -1.99826413 +6.1 -1.99809160 +6.2 -1.99794551 +6.3 -1.99782184 +6.4 -1.99771717 +6.5 -1.99762841 +6.6 -1.99755320 +6.7 -1.99748942 +6.8 -1.99743528 +6.9 -1.99738933 +7.0 -1.99735029 +7.5 -1.99722636 +8.0 -1.99716905 +8.5 -1.99714119 +9.0 -1.99712732 +9.5 -1.99712026 +10.0 -1.99711687 diff --git a/H4_cc-pvdz/pes_ooCIs4.dat b/H4_cc-pvdz/pes_ooCIs4.dat new file mode 100644 index 0000000..15de768 --- /dev/null +++ b/H4_cc-pvdz/pes_ooCIs4.dat @@ -0,0 +1,73 @@ +1.0 -1.89416406 +1.1 -2.03006332 +1.2 -2.12225090 +1.3 -2.18372568 +1.35 -2.20580385 +1.4 -2.22325188 +1.45 -2.23675483 +1.5 -2.24689119 +1.55 -2.25415055 +1.6 -2.25895296 +1.65 -2.26165601 +1.7 -2.26256902 +1.75 -2.26195628 +1.8 -2.26004722 +1.85 -2.25704006 +1.9 -2.25310597 +2.0 -2.24303439 +2.1 -2.23080806 +2.2 -2.21716895 +2.3 -2.20268609 +2.4 -2.18779608 +2.5 -2.17283680 +2.6 -2.15806971 +2.7 -2.14369739 +2.8 -2.12987642 +2.9 -2.11672522 +3.0 -2.10433183 +3.1 -2.09275639 +3.2 -2.08203592 +3.3 -2.07218609 +3.4 -2.06320395 +3.5 -2.05507098 +3.6 -2.04775524 +3.7 -2.04121517 +3.8 -2.03540184 +3.9 -2.03026136 +4.0 -2.02573774 +4.1 -2.02177422 +4.2 -2.01831523 +4.3 -2.01530714 +4.4 -2.01269952 +4.5 -2.01044552 +4.6 -2.00850205 +4.7 -2.00683018 +4.8 -2.00539489 +4.9 -2.00416498 +5.0 -2.00311285 +5.1 -2.00221416 +5.2 -2.00144765 +5.3 -2.00079468 +5.4 -2.00023911 +5.5 -1.99976688 +5.6 -1.99936588 +5.7 -1.99902562 +5.8 -1.99873710 +5.9 -1.99849257 +6.0 -1.99828540 +6.1 -1.99810994 +6.2 -1.99796134 +6.3 -1.99783549 +6.4 -1.99772890 +6.5 -1.99763859 +6.6 -1.99756205 +6.7 -1.99749714 +6.8 -1.99744206 +6.9 -1.99739529 +7.0 -1.99735552 +7.5 -1.99722213 +8.0 -1.99716834 +8.5 -1.99714166 +9.0 -1.99712801 +9.5 -1.99711922 +10.0 -1.99711666 diff --git a/H8_cc-pvdz/pes_CISDTQ.dat b/H8_cc-pvdz/pes_CISDTQ.dat index 734b0b8..9d6a519 100644 --- a/H8_cc-pvdz/pes_CISDTQ.dat +++ b/H8_cc-pvdz/pes_CISDTQ.dat @@ -49,10 +49,10 @@ 5.2 -3.93653352 5.3 -3.93724673 5.4 -3.93790329 -#5.5 -3.87997693 +5.5 -3.93851531 #5.6 -3.97676976 -#5.7 -3.89997467 -#5.8 -3.94015306 +5.7 -3.93961604 +5.8 -3.94012202 #5.9 -3.94038278 6.0 -3.94106391 #6.1 -3.93960957 @@ -61,9 +61,9 @@ #6.4 -4.03895963 #6.5 -3.93781661 6.6 -3.94353415 -#6.7 -3.88251916 +6.7 -3.94390984 #6.8 -3.94486071 -#6.9 -3.94131095 +6.9 -3.94462960 #7.0 -3.94076026 7.5 -3.94652098 8.0 -3.94774171 diff --git a/H8_cc-pvdz/pes_ooCISD.dat b/H8_cc-pvdz/pes_ooCISD.dat index bbb19d3..7870ef7 100644 --- a/H8_cc-pvdz/pes_ooCISD.dat +++ b/H8_cc-pvdz/pes_ooCISD.dat @@ -44,3 +44,5 @@ 4.7 -3.79187796 4.8 -3.77798402 4.9 -3.76474094 +5.0 -4.00645858 +5.1 -4.00462218 diff --git a/H8_cc-pvdz/pes_ooCIo2.5_1.dat b/H8_cc-pvdz/pes_ooCIo2.5_1.dat new file mode 100644 index 0000000..6f29c2e --- /dev/null +++ b/H8_cc-pvdz/pes_ooCIo2.5_1.dat @@ -0,0 +1 @@ +1.0 -3.44375043 diff --git a/H8_cc-pvdz/pes_ooCIo2_1.dat b/H8_cc-pvdz/pes_ooCIo2_1.dat new file mode 100644 index 0000000..ff03d27 --- /dev/null +++ b/H8_cc-pvdz/pes_ooCIo2_1.dat @@ -0,0 +1,4 @@ +1.0 -3.44264191 +1.1 -3.80394822 +1.2 -4.05321075 +1.3 -4.22379281 diff --git a/HF_cc-pvdz/pes_ooCIo3.dat b/HF_cc-pvdz/pes_ooCIo3.dat index eacd838..8404756 100644 --- a/HF_cc-pvdz/pes_ooCIo3.dat +++ b/HF_cc-pvdz/pes_ooCIo3.dat @@ -1,12 +1,6 @@ -0.5 -99.51299514 -0.55 -99.76770787 -0.6 -99.93901362 0.65 -100.05395002 -0.7 -100.12857177 -0.75 -100.17623091 -0.8 -100.20523177 -0.85 -100.22126897 -0.9 -100.22645291 +0.85 -100.22127021 +0.9 -100.22731775 0.95 -100.22684695 1.0 -100.22189555 1.05 -100.21391078 diff --git a/ethylene_cc-pvdz/fit_eq.gnu b/ethylene_cc-pvdz/fit_eq.gnu index 4885b04..3e416a9 100644 --- a/ethylene_cc-pvdz/fit_eq.gnu +++ b/ethylene_cc-pvdz/fit_eq.gnu @@ -40,31 +40,31 @@ set ylabel 'Energy (Hartree)' # check if file exists file_exists(file) = int(system("[ -f '".file."' ] && echo '1' || echo '0'")) -D0=0.0;a0=1.0;xe0=1.0;C0=0.0 -D1=0.0;a1=1.0;xe1=1.0;C1=0.0 -D2=0.0;a2=1.0;xe2=1.0;C2=0.0 -D3=0.0;a3=1.0;xe3=1.0;C3=0.0 -D4=0.0;a4=1.0;xe4=1.0;C4=0.0 -D5=0.0;a5=1.0;xe5=1.0;C5=0.0 -D6=0.0;a6=1.0;xe6=1.0;C6=0.0 -D7=0.0;a7=1.0;xe7=1.0;C7=0.0 -D8=0.0;a8=1.0;xe8=1.0;C8=0.0 -D9=0.0;a9=1.0;xe9=1.0;C9=0.0 -D10=0.0;a10=1.0;xe10=1.0;C10=0.0 -D11=0.0;a11=1.0;xe11=1.0;C11=0.0 -D12=0.0;a12=1.0;xe12=1.0;C12=0.0 -D13=0.0;a13=1.0;xe13=1.0;C13=0.0 -D14=0.0;a14=1.0;xe14=1.0;C14=0.0 -D15=0.0;a15=1.0;xe15=1.0;C15=0.0 -D16=0.0;a16=1.0;xe16=1.0;C16=0.0 -D17=0.0;a17=1.0;xe17=1.0;C17=0.0 -D18=0.0;a18=1.0;xe18=1.0;C18=0.0 -D19=0.0;a19=1.0;xe19=1.0;C19=0.0 -D20=0.0;a20=1.0;xe20=1.0;C20=0.0 -D21=0.0;a21=1.0;xe21=1.0;C21=0.0 -D22=0.0;a22=1.0;xe22=1.0;C22=0.0 -D23=0.0;a23=1.0;xe23=1.0;C23=0.0 -D24=0.0;a24=1.0;xe24=1.0;C24=0.0 +D0=0.3;a0=1.0;xe0=1.0;C0=0.0 +D1=0.3;a1=1.0;xe1=1.0;C1=0.0 +D2=0.3;a2=1.0;xe2=1.0;C2=0.0 +D3=0.3;a3=1.0;xe3=1.0;C3=0.0 +D4=0.3;a4=1.0;xe4=1.0;C4=0.0 +D5=0.3;a5=1.0;xe5=1.0;C5=0.0 +D6=0.3;a6=1.0;xe6=1.0;C6=0.0 +D7=0.3;a7=1.0;xe7=1.0;C7=0.0 +D8=0.3;a8=1.0;xe8=1.0;C8=0.0 +D9=0.3;a9=1.0;xe9=1.0;C9=0.0 +D10=0.3;a10=1.0;xe10=1.0;C10=0.0 +D11=0.3;a11=1.0;xe11=1.0;C11=0.0 +D12=0.3;a12=1.0;xe12=1.0;C12=0.0 +D13=0.3;a13=1.0;xe13=1.0;C13=0.0 +D14=0.3;a14=1.0;xe14=1.0;C14=0.0 +D15=0.3;a15=1.0;xe15=1.0;C15=0.0 +D16=0.3;a16=1.0;xe16=1.0;C16=0.0 +D17=0.3;a17=1.0;xe17=1.0;C17=0.0 +D18=0.3;a18=1.0;xe18=1.0;C18=0.0 +D19=0.3;a19=1.0;xe19=1.0;C19=0.0 +D20=0.3;a20=1.0;xe20=1.0;C20=0.0 +D21=0.3;a21=1.0;xe21=1.0;C21=0.0 +D22=0.3;a22=1.0;xe22=1.0;C22=0.0 +D23=0.3;a23=1.0;xe23=1.0;C23=0.0 +D24=0.3;a24=1.0;xe24=1.0;C24=0.0 # function for fitting diff --git a/ethylene_cc-pvdz/fit_eq.sh b/ethylene_cc-pvdz/fit_eq.sh index 3493cf1..e6f572f 100755 --- a/ethylene_cc-pvdz/fit_eq.sh +++ b/ethylene_cc-pvdz/fit_eq.sh @@ -154,7 +154,7 @@ if [[ $(grep 'a24 ' a0.dat) ]]; then grep '3 ' det_CIo.dat | cut -c 4- >> det. paste det.dat a.dat D.dat > det_aD.dat grep 'FCI' det_aD.dat | sed 's/FCI/1 /g' > det_aD_FCI.dat -grep 'FCI' det_aD.dat | sed 's/FCI/1.0E9 /g' >> det_aD_FCI.dat +grep 'FCI' det_aD.dat | sed 's/FCI/1.0E11/g' >> det_aD_FCI.dat grep 'RHF' det_aD.dat | sed 's/RHF/0 /g' > det_aD_CI.dat grep 'hfCIS ' det_aD.dat | sed 's/hfCIS /1 /g' >> det_aD_CI.dat @@ -253,7 +253,7 @@ grep 'xe24 ' xe0.dat | sed 's/xe24/ooCIo3 /g' >> xe.dat paste det.dat xe.dat D.dat > det_xe.dat grep 'FCI' det_xe.dat | sed 's/FCI/1 /g' > det_xe_FCI.dat -grep 'FCI' det_xe.dat | sed 's/FCI/1E9 /g' >> det_xe_FCI.dat +grep 'FCI' det_xe.dat | sed 's/FCI/1E11/g' >> det_xe_FCI.dat grep 'RHF' det_xe.dat | sed 's/RHF/0 /g' > det_xe_CI.dat grep 'hfCIS ' det_xe.dat | sed 's/hfCIS /1 /g' >> det_xe_CI.dat diff --git a/ethylene_cc-pvdz/pes_ooCIS.dat b/ethylene_cc-pvdz/pes_ooCIS.dat index a088f9e..f3e7756 100644 --- a/ethylene_cc-pvdz/pes_ooCIS.dat +++ b/ethylene_cc-pvdz/pes_ooCIS.dat @@ -2,17 +2,17 @@ 1.6 -77.29332495 1.7 -77.53168130 1.8 -77.70802692 -1.9 -77.83677134 -2.0 -77.92890141 -2.1 -77.99223605 -2.2 -78.03465007 +1.9 -77.84159820 +2.0 -77.93610860 +2.1 -78.00252424 +2.2 -78.04732377 2.3 -78.07553049 -2.4 -78.07339214 -2.5 -78.07667227 +2.4 -78.09104792 +2.5 -78.09691316 2.52053 -78.09714826 -2.6 -78.07260476 +2.6 -78.09548975 2.7 -78.08861735 -2.8 -78.04937855 +2.8 -78.07772924 2.9 -78.06394457 3.0 -78.04814002 3.1 -78.03100749 @@ -25,4 +25,22 @@ 3.8 -77.90893039 3.9 -77.89353420 4.0 -77.87879484 -4.2 -77.64245355 +4.2 -77.85124327 +4.4 -77.82613461 +4.6 -77.80327734 +4.8 -77.78246546 +5.0 -77.76349711 +5.2 -77.74618578 +5.4 -77.73036272 +5.6 -77.71587711 +5.8 -77.70259501 +6.0 -77.69039650 +6.2 -77.76688989 +6.4 -77.76929692 +6.6 -77.77128725 +6.8 -77.77294261 +7.0 -77.77432459 +7.2 -77.77548143 +7.4 -77.77645064 +7.6 -77.77726332 +7.8 -77.77794502 diff --git a/ethylene_cc-pvdz/pes_ooCIs0.dat b/ethylene_cc-pvdz/pes_ooCIs0.dat index 44bda83..5671b84 100644 --- a/ethylene_cc-pvdz/pes_ooCIs0.dat +++ b/ethylene_cc-pvdz/pes_ooCIs0.dat @@ -1,3 +1,27 @@ +1.5 -77.04903276 +1.6 -77.37200507 +1.7 -77.61389011 +1.8 -77.79366129 +1.9 -77.92573191 +2.0 -78.02113169 +2.1 -78.08832259 +2.2 -78.13390033 +2.3 -78.16289654 +2.4 -78.17926547 +2.5 -78.18608883 +2.52053 -78.18653736 +2.6 -78.18576775 +2.7 -78.18017151 +2.8 -78.17075380 +2.9 -78.15864566 +3.0 -78.14472842 +3.1 -78.12968920 +3.2 -78.11406141 +3.3 -78.09825578 +3.4 -78.08258197 +3.5 -78.06727022 +3.6 -78.05248331 +3.7 -78.03834858 3.8 -78.02492837 3.9 -78.01226506 4.0 -78.00037586