From e48f77e0bf6089168e5a96af2ea07e49445aa0cd Mon Sep 17 00:00:00 2001 From: Emmanuel Giner Date: Sun, 24 Mar 2019 11:44:16 +0100 Subject: [PATCH] added some data --- Data_DFT/VDZ/G2_cc-pvdz.dat | 67 ++++++++++++++++++++++++++++++++ Data_DFT/VDZ/data.py | 52 +++++++++++++++++++++++++ Data_DFT/VDZ/data_G2_VDZ.tar | Bin 0 -> 20480 bytes Data_DFT/VDZ/data_HF_HF_vdz_lda | 67 ++++++++++++++++++++++++++++++++ Data_DFT/VDZ/data_HF_HF_vdz_pbe | 67 ++++++++++++++++++++++++++++++++ Data_DFT/VDZ/data_HF_vdz_lda | 67 ++++++++++++++++++++++++++++++++ Data_DFT/VDZ/data_HF_vdz_pbe | 67 ++++++++++++++++++++++++++++++++ 7 files changed, 387 insertions(+) create mode 100644 Data_DFT/VDZ/G2_cc-pvdz.dat create mode 100755 Data_DFT/VDZ/data.py create mode 100644 Data_DFT/VDZ/data_G2_VDZ.tar create mode 100644 Data_DFT/VDZ/data_HF_HF_vdz_lda create mode 100644 Data_DFT/VDZ/data_HF_HF_vdz_pbe create mode 100644 Data_DFT/VDZ/data_HF_vdz_lda create mode 100644 Data_DFT/VDZ/data_HF_vdz_pbe diff --git a/Data_DFT/VDZ/G2_cc-pvdz.dat b/Data_DFT/VDZ/G2_cc-pvdz.dat new file mode 100644 index 0000000..8d8f26b --- /dev/null +++ b/Data_DFT/VDZ/G2_cc-pvdz.dat @@ -0,0 +1,67 @@ +Be -14.61741101 +BeH -15.18923334 +C -37.76050114 +C2H2 -77.1109468 +C2H4 -78.35649991 +C2H6 -79.5834817 +CH -38.38025056 +CH2_1A1 -39.02292896 +CH2_3B1 -39.04192746 +CH3 -39.71633284 +CH3Cl -499.44091892 +CH4 -40.38755496 +CN -92.49316793 +CO -113.05534345 +CO2 -188.149784 +CS -435.60599183 +Cl -459.59921683 +Cl2 -919.26412521 +ClF -559.19444863 +ClO -534.57426372 +F -99.52770333 +F2 -199.09818574 +H -0.4992784 +H2CO -114.22021595 +H2O -76.24165057 +H2O2 -151.19382739 +H2S -398.8641923 +H3COH -115.42044081 +H3CSH -438.05482386 +HCN -93.19016189 +HCO -113.57564188 +HCl -460.25501712 +HF -100.22866223 +HOCl -535.22969696 +Li -7.43263686 +Li2 -14.90139263 +LiF -107.15596609 +LiH -8.01472826 +N -54.47835898 +N2 -109.27677437 +N2H4 -111.56718685 +NH -55.09157748 +NH2 -55.73304031 +NH3 -56.4026668 +NO -129.59990735 +Na -161.85420851 +Na2 -323.73377653 +NaCl -621.59018413 +O -74.90996889 +O2 -149.98775604 +OH -75.55967664 +P -340.79165231 +P2 -681.72627866 +PH2 -342.00929124 +PH3 -342.63840859 +S -397.60087998 +S2 -795.33180604 +SO -472.6631306 +SO2 -547.72059833 +Si -288.91550937 +Si2 -577.92881827 +Si2H6 -581.6110694 +SiH2_1A1 -290.13799899 +SiH2_3B1 -290.10837062 +SiH3 -290.74828195 +SiH4 -291.39324391 +SiO -364.08184458 diff --git a/Data_DFT/VDZ/data.py b/Data_DFT/VDZ/data.py new file mode 100755 index 0000000..695e5a1 --- /dev/null +++ b/Data_DFT/VDZ/data.py @@ -0,0 +1,52 @@ +import sys +import os + + +system = [] +ev = [] +hf_lda = [] +val_lda = [] +basis=sys.argv[1] +func=sys.argv[2] +print basis +print func + + +filepath = 'G2_cc-p'+basis+'.dat' + +with open(filepath, "r") as fp: + cipsi = [] + for line in fp: + a=line.split() + cipsi.append(line.split()) + system.append(a[0]) + ev.append(a[1]) + + +#print cipsi + +filepath = 'data_HF_'+basis+'_'+func +with open(filepath, "r") as fp2: + for line in fp2: + a=line.split() + # print a[0], a[1] + hf_lda.append(a[1]) + val_lda.append(a[2]) + +file_pbe = open("G2_CIPSI_"+basis+"_"+func+"_full.dat","w+") +count=0 +for e in ev: + file_pbe.write(system[count] +' '+str(float(e)+float(hf_lda[count]))+'\n') + count += 1 + + +file_pbe_val = open("G2_CIPSI_"+basis+"_"+func+"_valence.dat","w+") +count=0 +for e in ev: + if (count==1): + print float(e), float(val_lda[count]) , float(hf_lda[count]) + + file_pbe_val.write(system[count] +' '+''+str(float(e)+float(val_lda[count]))+'\n') + count += 1 + + diff --git a/Data_DFT/VDZ/data_G2_VDZ.tar b/Data_DFT/VDZ/data_G2_VDZ.tar new file mode 100644 index 0000000000000000000000000000000000000000..bdd8115854c7ac467302878ed7363b3b3a09ff35 GIT binary patch literal 20480 zcmeI(-Oe0GZU$f-&sFFKR@+tmQ{4wSSPL+8u#gSt0t_Wq0wZ=1c%8)}SD%Mf)jwl> za_sdNCo$&3V&$IZbhB6_i|zSy#M;c_kVdcmb|1te*ON3_g}sF{Hwoy z{hwcb`}O-5zxr?Rw@lNE_}AZ${~c1EUZioF=4JM`ycp70QhM?7Z@>D@KjQ)a@s}Uo ze@~^KA^Lk1$>rsDKYa7uhkr=(JWhGcX?j(LWgFAYt@lsdTE_J@4ln=diCb$)d6{m% zM?3vq*Yxbc{)C|-=v(}kX&ECAe_6IF|MHI={(HBu&%i^RlJ-h-FTrtCB)w>LPQR$GmZnDpJzs zJL6hKmR)4avsAS;Ro5YT#x`xY*Sn_6Qd05U_8QNc)-9EhK>jNIDRq6z=$x5zZc^iz zmMrfLRcicfm-5lddFTGDEe+dNhIrYyKWEvt!NRwyu`_g9QrdDzy#BVyOhZ|wWnQQH zlr$th?@NVDU-z(D^f)c++N>)ksJ+nUZC~}YW&-mtHQO!4i8*b|`$ZR6<~$Y=)-+vN zu`o8vebLRI*@e}{p&HI+4BL?Cb$k7;Jz^}=$|7vuj6c||d6)_Gb&Pq6wQj>u1!k`s z$9ZJgZ~HrzxukJPY!nj7fl}5j=b6Ok-4@fhrO{IAOL$$Zdf3?F+rE;?TaxprOC$@M zsg`g32`OdPF5}W{Rfe>!_F%q!zArRy3$5n4DZz{^V;^t3gS1%n|8*J1)$L6_tm%Do zD$_hy5$EAx6~jbiw>Lj7PP`ru%d!lk9jFq!9>X$b+HdtCjKL4eJkTTW7Ee=}w^_Pt z2DKelR`QMSe9^T{>o^Z%>h4SrGQ+r3nRi`=N%_LIO-)_yzRY|n-QFDjGf%Uuu8_=t=AgTGq0N$kcRMi&~&eHs;N-{UbjYNPGNFbCJ(yn(R9xt<$5dO%3TlAm)i0 zU+mAA#xyT1*PanJNdpUvw?BW={KL|wd9~k5RrI(b(9JrOi7}KaP2#vBvFkB!i$}0r zJR_yTHf2hM#O?yF=)u_aEyvZZhno^;V|AWf@8?3+W)8QvT?soyHR05)B{Nrfj;d64=rVIxUtVK&f%0A8O`Y z#<|6;pjjD2%)Z^XwXB)V*Pg+)$ng>=f%nbec`AAHlA(G=9z?v;V57zRzShb9lyzED z)xw6!h|)_Uv2W`fYZbV;+KWD%iOX)cu7y-G7B z^Ut><43rH%9hwq@1##g~ynWf1AwMoF&uNCpgH*)Dx7(L}71D1ON!!|X%x+F%d|Gc` z?$*v530~6@&%m|EK@ylr^kmOaQcYfaoC!qI64-5z=V7!=-)YNqblCHi+q~;aRDUQ% zTC9Gxid6rb*W0^psA;erl+0}z+ZJb%r4ZQHQJq50tNXK^ObCbxS1q z$3m76*dM`X=6Nf*C4+Htm@JF`Vc(rX+og#b|)=3Xg@^sDl>ZkGlzyIThr@!}0K8^GBXXF26Zuq}E@c%Ed@!u=^ z5&q8!77$WT;Qx9{2u>g4|NgBct>;Im{k@y#~qf#BgJ8?n7zmU(ELM`5B`QYdDizLij_Q8a=I>BRs1C#zvw(M(X? zs0t!*hySAhQL*l|1em3!gjBC}%K+EGZ0hn>-W2FrAvBTdWs4Q4FGjV0l^Agk6>NmW zUL7mx(rOiL;qZT$YR0i%@xOAkEWM0}Euvo9=)B|q87;0@EFHh8a)OEoc;D*Z09vyR z*jjmu@PL7qD$uLGyyjW`sv-?d6PCgU%Ys7a)hjTso;oOmqTp7emka)1gVMpaG#>!4 zH=t7eb>M%9*w>OD4e=n~WSM*Vp3HG|mUWiDx88M{XOY>4N?xvJxP-p^@Y2ME~kl5!(1V?2-TG zKl6$H@jeJc8M7^wap%Pb+@XJA0_U~7fTM&(Bad9rKeWn`ySV_*vf!<@EBfD15X(-d zfJWiR0``Lb6N2Z6BXH<>NldxUS8gJF8ZlI4E}XA3V12$#K;+foKOL zCn3^J-qCsztIIEVOoENMY9^1qY z^v~?>3~OiP;dji5{S2{wj7nc9aA~fookji{NIpb5`gcxDarW#YlGuVL34);e8Q?N% zJFsb8S=@A@iQ8wdn|?13{0aR3Y$E`1-2QC*UlRN;%*O}*|1Es%WBfnA8nE5?VYlx4 z|GC~m@;-t8+qaZM^5t)#_v0PUV}Tw2KU=Urq1GAypTPgu{(p=AIpqibf5Mc{{?UQ| zCtX#@yW^ST>Wcq015|;%fv6+iu1@@~NK`xavsaO1u*-pckN=~p10n123@+&(t3B@U zzq)QbS}<@k{G{T#!~dK1zt%*nz13_GX26US|7#B__WqDqhhhWEILFMspWaW|40#_ckM-=x++MLQ0K(|Y8mwn^b-qLN`bdh zKqtAJ1ZuT!Pxe_}5c zFs(JG4*ahLWac6Y)M^Xu99402?@#>iEp}8jUM27*8NSfW-f^b8>Q$OmEeeo{I{{UX z{=XtHnvJ>D$o4e`6^4fY1Fj!`^#RA#l>yL1EwKCCOo91Qd}5IT-*>)69S`dEbqs-v z!*N0ClL6p8QP&O?lo)e#K4J~-`~NYt0nb-WW{$51X6)mx|L=Q9f&&xi%PF#T?oS4Q zF{-+K2qlOHXvs-^?f>I5pvUH5xH;g5&KCwCcoN0)+RUo+3@+0B!T^|p04{pxQBCKa z0Z>l{z{r4s9i8B|Gn(7rXOIZj{y&Ks0}&!T%vuscYT^RsiG;nX@n391Nm{%-Z`~T5))swH9h8NeXg8?vA!s(>MdS8|ZE<-$?{eQcU zyyy%7N)i$WlkOM*jZqk+G%bS1kw0wirT-62qBENkqs@ozSp3remrQVued-5U<~17h zR|dd*2~GA15>cZMFJWH${~`caQ&+MXM__$y+Q9&r@(5nBu~1?cVaD5yDO~#hLWsLp zFLgH9FaY|>07!Z$C4tzDXB_m|G;lTiSr*RH>x9x^`2HbIs z)E*4LzzG80_rzlW!1;qJJ{f?B(N&UYSw$2^j^%kf1K{nrJlsWlgcFqIQMm*<1AsqU z;**vj4LV>0qE#Lem;S%Q#3V<1itu4u|4-(6jLOR~(KC#Vfs%+(7WmvIcATR!c%&!? zLyM=NB27kx(mMl?F>GeuumnBiLT%(nh+O*rSTvMjGcV>)>X(%d20${!)!CD#PS7LH zr{U3*{6yvW-{b#Z-xWOT|9a9MBeMIsdpG}I%m>?4|2*dZKe1=e`SHj2e|;4*IGPY2 zn;2!;`f37pL2u$hoAK${Qv8N{Rw@~`2Pg{zxMxI{G)9j`2Pu0KKn-p z{s&$F=3c7?lxbI2{GSFTuj93zy&3}T>cs!EN>nGb9p8|wqr6UjU+4cF$TL>HJ|Ze3 zL!Z$;zV!b+*XY31PIzhn!g@J(_8LyTrdN30(aq0iJS?OLpc>?&NM|R+Ub*F-=PY_7<(ZunR6aPp5 ze*ue{Rlz&6$xvc21RcFA`rt##S7e27iXUp|uwU^1wwR{IBsV?dzZHdATs!{viMkt# zsVikw4K2SsJMlk93^JuPufQ2&*ipOwlKy|dIaO`EA@~|4EdX>u|0YcA2Y^s@U`F5s zOuM50XllaSebPNTgJX~(MEbmimha}dXa66GfnzXnf0_Odun|q&`XbN^C=WjCivHs# zGJrSwtr-Laq`cAR3;NeYX2gaWY?_b&M(J{nj{t_~qyD1`1$5A6^%n#?`4f#_@9*eH zM>u}w5^ve}#lW_H=RSArrvk`5Q1gQRO_s+LH*TgH7U#vRxV#JcuYFev)}=)MKYmXT zKM&~GMMRE$mJMNMSzbB`mgzo?mK5Eg~Zi?{$qlG`wY#)#6SQp z;}Qq@FEO_eyVoqf8CeL1=nnnsPsZp$-4uAtF&M4-EBZGFVK`<_9JoDG@^|RpHUx3l z)I?}_qQ>|Y{Y$98e_h2RJFybzO3k@XrZj^C;nx4@5`=-&i|+Z}rQ-Irg=I4TbG9}R!x z2?!eR08R0XF-~wn|2z|F?)(8)m)>Qs9sTS7+YY}OYNnDL0%Jq8m;S%!3*~k^As?Jy z63vPJ8O|@2gihn)P@ggo>!0X9`u~!$@wB?@W)y;FjOqSPkIKiUs?3YD`+#V$upo2y zU5BIJS0RGh1A@c$bzyj*f7o9-6hzeyF;;@2sPw=N-8pt}FO7?dAusC|xgc;?!eoj` z6tgB(3BPaw^O0c$`Xj=4jxbwVCE)JG_VNA;`q%%L$&IhH(*$9F+BCmU=6Z