From f58715d53c55902ee65edc27117d22bf1bc8ad13 Mon Sep 17 00:00:00 2001 From: Pierre-Francois Loos Date: Thu, 30 May 2019 11:49:37 +0200 Subject: [PATCH] CH2 --- Manuscript/CH2.pdf | Bin 46225 -> 46222 bytes Manuscript/Ex-srDFT.bib | 6 +- Manuscript/Ex-srDFT.tex | 242 ++++++++++++---------------------------- 3 files changed, 73 insertions(+), 175 deletions(-) diff --git a/Manuscript/CH2.pdf b/Manuscript/CH2.pdf index b15f2a761d539abb71fb45009e9b303810289545..414ed788443c4b2de2bb7f499ca54e83cf822ef7 100644 GIT binary patch delta 11675 zcmZvCbyQr>mUSSwySqD0(~WzAyF+j%cyRCF?(P=cH4xky4X(l6LvRR^;l26Znwk8* z{-f5d-TRzrR7lf~nX5$=rY;Ck(m4jz8gb{CV@@3zofv9hzp=~-3Ay`oq@OY{3O zAKGHbAMH`MRDDR<`#dQlmBTEc-z?jrz(MiXldeg5ekR_dMBF_QgO_I`4oLTz9dnDc z?&-y#%kcg*LuY~^&rtgFXM?y-$k#_Z|1!K~g^=KKqf2A;qxK)4*v+>lFYo7rhRxZT zcdE@qr8w))5-`v@mUiz9Miq#yl|qM+jb_ApTF+XtWz3CW=Hd!8}FeMsi~jnm6{|VZm=m`)-ke zIthbkyO4$m&YZ8}bz9VDKM3~$;|N`QGzF|V?c}tx=_JR;9x>tXoee8Ot@~f)&#q0X zGd9mIH;En7Qd#^|jNsHjSt;nU-=6jYyYqEw5uZ1{+(nIPBoE76lR(Tp7$4_~^beFA zr1#=gE4nvC3Efn8SA5{e&+lfc>8FrCKr1~ua9rN`SGOQk1@CIF%_9(x7!BuJhtC=-S!&5Jw6sm zk)sg)?6*I&#>8uKYmyY@y649nH(EqrUczz&Y>zHi&Q~_h9&=x=SGuQtYWdR)BTKe_ z5#edhDp7}%r2regBz7t96?JUVjoT@K=`zc#QEFH|(QTva9Ir8Jqnp-Yql-dbw5Y_sdyx9PtrzPm{G06%C%c}_ zPpaC-D(R_pazCd*h4ik#u#rs|!ft^MPhGHJ?8vFZqrMVAiAmVe+ZDctS~_&DQGBlDJY+DvMLdQwoSdSMnF%*_^F$Rq=>G?Cr5gtIi{vE z5YC8p5=eOV!?_d^U_|oR)-zY;!o|ukHIa@_DIPOKbO4Ln>{loEQyIQQS-&Q}mO{we zhoyD=E^zE_N7P&V7}M;lmA5cxuDC#>RH8)Le$KYkIdwHz^fWQ+&#$mQGPOomcbJg7 zQqk+TZ9<0|MmODmu1g_LXOrhvNi&^KFAQk|7wd zSZ{-52=b*cI%^)^5nTl3Z=!ZHPkTnac$e;KWsy-PuZ|1d?C5HCn7fU$t{;*_YNL;8 z8GzwxNp?G<6yZ6#A8kT1&2+0hN66As1uQonwfFB!F~_H`zV* zUML+-i-W zF_Hx{tNQ~+z23b?vyd8qZ;PZ){~_=Ze|7j`1}{Bw&|2OCd8T4p{D(#=lRXPTu>pNw zQ#r&s-5ot{_O5a`VVPbfD^gNL-W<80Vw-wSBbCU8h2Wh5ePCI&hAkF2CcMx?emP+g zsi?&G;(Jo~n^G2vO$CbxrxH`_2~C=T~fvx##Ih+j{PC}-<~Yk{z4dP zG5MusZ#OV}1qYek=YnpKRBI!tQlz(bj~pRYadWmCSplSOy=Cq(-{BnRraRPDHD#zg zl)!}7mHZ4Rf)3FmesRz-^{2!1A$QLBp#+%x(Kwe4Wo&|&3c8Y3augma!`V4k;LWen zCW+N3-rDJFq}n z6-d&Z2W53!1c`WqiCLf{Nbig!T;5E#DjUx6geGnD*Bh(#NT#}9WAnE--_`?orZc(H zSNnH~eI#Ws-CSU_MoaL`_7!L}NixZGJm?iG0y|l~ciirRz(*@zkv+sJ!bdR%(osBK z7-ZkG^YeogbiIR2ZDV5Amv|7tlemz0sZG|^2W{nuu20WdXn74!G3TV(Y_coL7aj*G zii&3hD(jL({8gK}N<1X9syj&fsPskObmmBRkA=HToVE_{&=o#5XsbSOEg7fk3))`B zMXcb%B17MEv&*K>kJJx_Egx$zwliTg;1!Ho!~@9b#a7QP-UszF%HF5pe1YNCH5q!r_b09JP3z_-*-XoNmdbdit(=?qxp+ig zp~K*%MHp^m<4=;8Sn1Z=z`B0O`1gQ5XaU5@1wmyM48O`jI;q4NFIgG##vfq+AN=3qE=iFxKTxfQy#vhZdhhWWwu$?Lo8+A)}BD zD4%+UVs6&nmjr!r*ZrHDL^30Vp2T)y5pO$R(!CNOITfhjB{o-&C<(6=NkFPjktf`cc}heX)o z^cf%D_%7psknUQy5sbaGNUEz7J#k!Gmh7rFLZM&18+~s_&aC)0Cz>DFH$_VLH~sgB z#$LYK$~u3ti|l_S%Dg(WUzl-bK3~N~W|YyNEn4ov?Rj)RdD~L}?WLy2Dx9(JZJ?t% zyM<+@G;?Czb<{e2ktSan*r<7tCW@skdQ1>TJHKR7l(-Kxb z%HbGx>-Np?sbG{<4^Aa?U?>^`dul>gkrnU+PzFt-nLO)Zu_6lj+xbgGBzDrpQ=Im!V@i6yU(zTJ7>HcZU*|jqzp3T~PB1VQN56gA1cpQv^_qjPw;Vt#P8G?uuf2ewZA<=q<*SR4ypEru1`?Y?G$x>og{+rxB zQS;gG%X&e>6)0~0$#4n2dVjVqH7x7H)rBoD@1RQ?f;G2|NT;u|)>ZxHa`nDP(}e2P z*}E_MoMpePZu)ISh`Vf}f-dFyS-XB~%=vl%o!f)?3*?LHjSdtcDMe;UhDs2eb9)0s zVy3sS&BrPsnu09fLle+OV>jQQr?ck4`i@iKvL*`$T$0^z9r!yHk{haX*{t8kApUMb zvQN%Pn34PjpQ(O*y~%~170hs#`Pj?({#U%J1u5V%R6F`DlHR=vES@Ea!sk!#cQTM% zO~GTZ=6GR~_-t=`&?jS{^%w>a?2Z6LvlPeOg)+Utmm3iQkgsnpCKrVZ@$YQj5Cgz|Q3eD-4X`Jg&2 zi>3vprHN%hWoA^9XZDksq6~tIH3CKl9C{Vl-SWf7)gLwQ272ni?tApPThf#?7d?aD zWbCa@R9kqWyy4_WNp>s7AtXI6bAgDwAsfb1(Ab_%`&;r3$K5wf8(HGSl~?My=Og%P zC^=wmO;AXXnPVdcpJl9xNlQp{4(X(_j(Qj&p6&4}4hEQt{B$@8ajaco)1>43MQrEZ zrY>$dDqNl0@1Ok}^>Fjl#N+dQT~CULJ}4Wm8sX?6ot>y(+LfwHA!Z8e!hxe*F|CkU zW}Xm7;GmXZs4E)FJau3QY}|M~cF%j=Eo2~UyKC#&n#v&8*3o zV!2tpcJy!9`qax|c{0kG$KuKuHA$F1WsfFhaItR1vay3uAjK|r+XR`_C~1381cN8F zwwSlbAvSLH37(0|C4HU54uy}#;ZCd_BoK!_T&*Kf|$A{G-7#SOaDH{pn!FB#?T zrwUW3QXSql?OcfrU}uQ)k$FlX8!U|AyQ4+R+Au!MLbRJ2ICCHApL#p=RB)9rJ~Wr! z2T=;$quvMUoL0Y-jQ5oy1;5T&iGuGe7Fw;Y(Sm*j9r(isIeFn{8lgPraZ7f4{yexs zq_j^MB}5x~HPra0y*_$)hBLmQ!uP!{-GRcMjxnQezgsQ6=DoXG4!YwkJZynX^gM01 zQJO+9Est48M!RPa73+Q>f_vEE`r`S5RNjsEQM-;bvi0ekQ}b8{Gayr#AY247sz%Ek zCZ28I^w$V@VsssGznEc(I0y_>Axh*V3h~NY zBna;7%exdR!jg=rcf>+a#L1TM8Z7L+7nXvUK+C=^hOn0fx=k(?DfI>Gwas#=x{4Cy zAWnZBN)v5B*$JUOT21={qh<1Q>8+Y_1By%ks@yYIWNITy_t;omJ&KE?k;T5S)|RPO zT8H&Qh^8mB2DS26vPqEUeVav|C4QHU?g|c3s@q6qD!SIzL8@P-rsrrwFN@GlBWfk& zm*m)QCAqJ5G8K+X-V?9zxb@P$bgEyZW~=!yJtx6UAB&KB1cBC8Mj?jQR(ay<5Ji8M zWw2I9?KNHZDim!(sI^7h&r#LKJ^j(--)gjE|53|0#H;F44}H}oPBYL6;UrmT(P}p^ z_+@JSK__G~q(m8_8zS$oe`(6k`5F(lHo?zjJ~VzE4nj2VZTMIA?*LWuhG4ro6lG!_ zosd#A(pOr+zhta+Y0+Cv41N_Fyvk?vyvlG}zA9LextrjNm@HVmriV|P)CplK*C${| zLH-GiNJ;wO`~oVe1}xo=u!o^T>Me>x5>m&=2)gK)Z+jF-EiDUMw_N}mRU%< z;SeuJ2|4F^Q*orAC{U1G?>C^IO`ft=*K&v>=Iu$X=MzPrPlT8TJ zO0zd7wu+WbV4kr~-|8TvLQ)#8GPd@HXwbvp#A}czVg$wKCHcz(F(=nnS#>i{hVgff z+bD_?#_~-i>DzraftFJ1Y1xh=_i60jHKDmpc4iY2{0@gqm`E7`aa9F_gEl2?__M{?hM3YGW`jiL>KH8$mGAAb}$zIj!#Vh%UQ^^v6(;?fl zu*_~w*@eO(m23WNm_ilNRL97ImGx?Ugo{%Xb##!D#Vjz)RDcuj+;>cA0d-Dw=;zR| zhTZ!{%J|b(u(`tMy&FV4>yi;{ai9*(vQW3)Rp80uMwlOJme0*{X4z>B#{NjjPKMPY zLyjuonN?we70A1#**75(m1zD2FZu&@yPHXWN6W)eV$}aG8<6n%K85bn(#Zu| zk83og{Kl;}owLUWxMYFnO=8No;Elw1$-6R3?{LohZCx2O~@V^*sVnjwl{;S zjXyOomA`kyd0fnCP8D98{MD#x;2#Px#th z&E*8u>p_ecb?m!X?(NlaJ|oXjGJkTNf$dfJS^d%3tY$s41m|%{7&#?zO>fpbHq?hT zh9X*J+#z_QR~DDCuYU%KWMH?koptZ+O39qTZ14AI0etH6mHnFbYvMi=j{}UT{+}I2 zYver4G3sPIE@6FSseYf@uXJ3L8Qo3CgBl+9+fpFHYZ7O>Ch9e#XCU)krNJ!K4|U|I zW@cfUvb-F_KiJ<*Vw6$>+GnTi8yEE8Q;X<2I@_NRRB>v#94d_MIJ6@fDG-?D`Eo3FiAr!+n7TKRyM_zyN1j)o-nUbJT)w0aqgXrm-FY7F)Ij;IlEtY) zo^Kg4I1kEwBgv~YMWTFt7tY91=h>iaW;;5!U%-Fu*+A*-pP+d!2(yh$4?_|#7-lfR z@rf@O&sC4Uvpb0I0?Jkt&cNb?a+rEoGjyv5D!{ux#EW1E_H6myWt1a+enW2xMM%J4H=313mV`w}nAnaR9g2^a{!sc}nk4-LyE*mhM#Vp#&q zt-tx%cOgL&8O>#*)AISaQH~IXrr!=qlA^6mALHCo*PA4`pZU=&26hyiu@7xq9|KUm z%SS_YB?(WZnux$FCu*vNMzv~m4-~U@Vta%H@|4D9z7HGx{C%m&Gv*g%iYN*Ura_4K z=*C+m{h8(58nlOQVFJl`8QcjiS$3T-Se_wr0^e&_9hs?!^I2CvgMBk4EGY_{WoN11 zF;OSMnw)Ohsg0|dG|{@wC|Xv_X-+T|QADeKL3i|ZVdK)4DIR7<&d~UF5jKpC*0*8` z5lSWM4hD0!6tNML(9k48wutjlApM>vvfN?nu^K$&#Zsh7}h5V-~MKFI*t&_3oUFVfkZ zZ@*`0XozPf2B{5*Jv~PBdx}y2kR1|mg-D5D|9*t76bT6iU*+DWft{4n zM`k1Ns+$}0VIAwMhr$ls9EemK#;kyZMB#CL@OwP9xa^c5k-QdIs-Zuqvy*6!eJCgAJ;_a94-R7F_D zNIpr-X*i4;S&Uz6bB#>B#0S6s-pHjBa79f3Pl}M&gQ$%hfseVfi*=TQrmN3IzzQJ2 zIh(|in410>n9p8Tb=&Ti8p;Xdq%LGWGrU~&WHTUl6v-uysf{mlF4jmMyucJ)OH>Wf zv9}dm)XVOwKlp=~ywE=Ma^(}ZcDNxr+YKpzbkh=Y0!-x78q>3JOasr`>WVQeDMBhx zx(@o(g2rgWR5W`zcltiG#LRg+nR^DJJ#Gbc2cqi@N8NY8_~~x@KXJP*VdgzO_%6XZ z-(PDVsSR)k{|ms6UK$FNU-j}1h2TDJC3Oelw|7NjZoGr@2_<-(4ie2H)T_fTDrsH9 zMxWyzSc;=~!XIBc9}j>T=AQfaY;_57(q>iLCzR~pbn?f;r;e}#V((LF%oi-t78OD{ z290D&+rZ7C4+E(Fr!-vtRg&48L!rC6{&wOX`8A7nnRNH{7I@?kL?Y1z`rf$aW2~p8 zS8Ed|*`O1S>|Z4tNU0~4>hnH^f)uvhzVUkGZHDn{II6qdN5Q*7|8z<{6{z#ijs_z| zoF^Thk;0o^4xy6f z5E|MQe?p8&nD}K%E6U*6oKN`SG=rpR%5)OXb;?vZ#RW2DYQ$!>qI~xH5USGMK&Rg7 z^jn?jlNRrM5Nzz%{SYP?>6`I zbxlsDd8qPv&Ze}2uQJJma;F)`jViA)F0VuKr9^~Zb8?CbwmZ0Y%pG0ol(`*&_BPod z2pi(Vz{VcEM2KEVhxfH|OZ)pbEYnH2Vyc1M`eTW3;hh>(>RaSWk$2&K91Xrw;e$kv zI<3}`i`0}gxFe`R({Q3wq! z1sX~nrt1~mArJntJ>xBG$Q=Wl@3!$el~o|7Ua%IO0dsJ+{wuwWJpVEiAgXWWUNFgVdIQ*r3i?yz7(%}yXa z`AJG!BC5p5HzJV>YfshbHW}N~Cy3#o^_l+Q2f-cWn(sNRx_Y?Cl4d)1LK9(Q_=^kB zl&WUrVrvKek<2WW$lHCBL9dbfna$QjFkJ4;3{HCDSGbxc5PXffqBxVJ-lp6FL9iMK z0rO--9-6k$lfowFBPXAJYs<~9Fw~iWU1~70?aCyX|LR)i3h)wsn}GdCAon;G1|z0r zrA=_E+6`2%wnb8Lq_FJdJ&_weYEc31r#Z$$_E_(hG_hkjf=!`f`aZ9wX51&tSjUZ1 z-#_Je1cHx{;mkxbRngFF`$@fEJaz@0$F=mj&=lmch!F+?hqu=T^?PJEMN!7;_v8z3 zo*OZC<8b;&Y@c;Z+wCjiAa>MY297_H1_biVe51bd7r6xJZ5DdCvj-OP#)*je>qee@ zsylGvj=(X25f4{cXQguY*wBLE`6Ymz3J*5ovFHa<3=Y?!c|FoJa=dv^Q={%T4N&;% z_EUCY%#*vi(S-VAum9kpD>Xc7l^n9Gcv8D+Q+}3ycbT2eQVUU0GsI&^Dl|mIom9z< zv~6`kCB?pumoRzHpOw3@&cqA70ne%*_PpcBV5wHJV#J*h>q9U-W1e4}`N218vwoKK zp9Kv@>{}~=m}?HAabT;4f?TeR#FFC%N6n?W3tb7VrZ@sc`(KSs%fi(Zg8C`}~UudMm9me=jICF*j zAXBN?Mon5;K4R!v+wG5u1X?T^2i_sMf+=bG*LxX^m3u26MhGp9AUBc$LfdDFPomeR zDLo@Z^=J(SU6t0<2+QY6XutrRAR=hUm(M2w@T3PNUphjMG~aX|Q|(AHosvNbgUiSd z@Wm2?h?*oO;!fAxE#r{`t#q8o1|-m6_p7%$!pOrv#@#-{#t@k=MTj$C9b;lycB8p< z8VLZhWx~aaQ6Ve3;o>2+7V=VPux!#(XaO#mnbb~-P~V-fuvs+zJsi-x4#KlWKv3(~ zvacPVLZ=VFrt$t0Egi!mOEa0?!|1p{o>_3XCVJS8J(7Kyw<RT%fDi~irwDNs7B3^s2jo2oD;0HX~)*F7lx81A#@C3nM=Krx>>j8VaX!?xzNv!Q0 zO09I1ftrAxxPIvqpR`d_jTg2Sgiwi*wJl|QvBXANI&a23Wbsw(c0of(s}9#p(O4wo zdR_mW&&Q~tpK_N3afA4B_S}SPZP>uL;(~r6TA>zlr4vhO8|RJ6Hf$sMc6O(>Exakm zMakUDpGpu(8J1RlyVUCSt-58el|-#uQjcYmC>)bP6I7$gGTBDP^@mEq(neFR-@OQA zW`n6bt!ubiY+Z@x5>{GO)MkS$?pkcQdynbpJbTD(`AEW1kJgHr0wz$F4^2?9(0%!1-9oWZq)p z?)N1Il@}hcUHUc$Iu0}y7JVz*=A&w-(_+_3$A;QTr2G+{Zoa)Xa*7&71)X!=@ zQ&+oNcreAs9E$VSnw z58~b_GD1^n2;Bv^TZSc_{cngiG$0j0{U2{J;!MLW$23`!qp7o(;B_JW2 zuVvJ0Jb?sb=Yy8A7_!;SWK49zuHPH6%>uB4qKU>LifsmTI?9$qoWRT*vK&#gI&Hk8 zi^G~nd3rzL*F3TGmX%X98G_8k_S=!TQjPj!D+Y9smJ&jUqV(NFD&(J_Tc zmTdbAyLftGN15jOEsk?4GZ_w-*fm237pD%Jx)BQlgpRPH_k4ea6$<# zNu|3ZM?_41w$vycqIu}E0y(MdW#N(|mEW`&zexC~iM+qB-=~H`;=pWAONCZ!$f5)? zrEgugVGjE7%?2)gE@1~mLZhGp{;^(}&{zxgl-*MuU)GY+re3PxN00NfU%3|9a+P^R8k_6K9m7mN#uwHjA1zk1Vl9=uZM4e_at5IRZTgx7u$ zBYOx1D%}BWMmG&)v4Peo2pRd_lWlJZ){5URF6rj!-|RZ?H;0F?iJG5~D z?x>gPs~t6MGIWq!otb>*>R5xVK~g;Jso^4>$0c;IMz-C+z0tOi#7`lLA@a#VdvL}u zhzFPi-j9GB&=$|W@n?f0!x&&n#vDJjW)ych9Fu@tXZGcr-nVAKPpWia$dd>GspcGjcW? z*?;T}JNlJ+`tpaCmw*#dmd3&W^AF|9ANV0GzuzcHf-F=p+>QDy(IH6O04^?WZXht} zS1rN+Ofc#o%-mFfe|@MNs8IfW!NtYJ$NM)1-~;g|wbU^}z<+M@^zcI-FDiZ|wJD7_X1pY@mydd6xw8P5>_xK!w>*~ix15GkB)&LuK!B}0DwRcH~(Lq0D-*!Ka2a*WaCh0ARkq&j*ML z0Q&oUfVjc`IMlBgKZwfg|Nd59cNp;Z-o0YLzdHee!2c{71m^okEEvECND^pLMgpVI K(@QE#q5K!X>u*#5 delta 11711 zcmZvCbyOVR(rq9RWN>$Pml+0k3GQxz83-2KZE$yj1b24}uE7Gq-Gh_hBuKEwefPWT zeK)_?f7I&UXIFLYI_Io)s#nJy!6zQUC#b*_Yw5ZY*F|A6?2iEoqhMSjnb>c( zJU5>#y4S;=uX>E1ZWkXfYo9ZVicOy;TTGvno}WJ4hr^0CZ^pzN_A*xYdQ83B7gk^# zji!b|9Ge^;ntqOlZF0btUmQbRwd!wd1qoyZ3tBm&Qv7T38L@2c-F~^w$={Vha*2sW ze{8%x+I&uN&Zb|aWSK6VCdB(wqO=JdYo72rIx)C*Wk%X#Fm^78&rE|or;Tko2pLA5 zs1x}3UXV3-HOV^^rG4YJ6nnP_IfI>uW>T`q@M%(;{YVq8@e?~X{xE|?p1XfhiIg~M zpfO}P>>p$Du2|Up59Hg9oS-7XaBGVoN1Ms);dMz)twri&N+#7Cx93%Q%9D&4wFh^v zGEcg)Uq^G^?C|an$|AW4S6gYMvILi8?2OkxlHQEZ&aSa`j0YR~cg(5e&S#)(aA2BT zxLcRc4F8g$+H277O9(Mo&QU;q3gYRa_?0VvM(9z_6rVT(C2TlRX(@-vwW&0D#BRj$ zJU%|vGs*sGHz+6CyJ*t1OuUob3r8)6@4Puh%+H7a_Nl1aFR;mj1^;9CHNMiTJ*32+ zYO84PjkL73jbmhXOfF9@29wIg9IP1un-)xTzXi(`WUfaGV9%?v^qO&nB^8gH;V!*i zE&?Odvu+wkpS|x+R@YRYZ{nAHV>Phtg=@r(S4i6 z)j`{l5Q}@U;B$q=sDy@crG@cjg$45)4GI#OVNSOD-#2PnAN*+6a_j9kX5DJ%dVMX{ zC&vV5b8lyKTmrd&IMVUlrccAM*`IB-B?7wPfdSsvfcn}xV&0iWM zN?w{f9mGX!4DrJV^UCUbG?JHov>~SmP%v^7;L?kbf}*omj&&8=)q+=64B{m(3(Uq; z@Ezv_CuEp@fQ`n|+p;R9Zey0OhZycIkYz++)zv3^dxBbJ>kg2bRJYN})MBgE@}OT= z*cHFL)^F+$k{62Y&5CceCS=Nj;?~)So;$79+b4e#MLWTa&&_=$httd_0&|8mFRVKW zb4PyG%9&QGzcpbyR-ll3U7NIDZlhKHEPD9ei$io?u>85FXED2-e3j)?5(^TM7=Z~HMv_bo!r&)=7D2ylk5IDft@<@h zY|Cy!vu;ZC8Gx0NuUL{#$+OgJQx+7G+1cBuZinHJ$la0$v$VB`aIptD@5hWD*lmZ{ zV?O1`p2YK0#q+Z&?9cG#(}O=BGumRL-3)S-3LzV;aN#waPe|0Zc5s+5upb${ukwwh zhz@#RblBuUWmoFm!@&v?p6_=sJ&E76>Dk9Mo8+C~q|8sHfx?whX;m})2oA33gOzVA zDr-oBbkQ$2220@6HvQ%$ztvY6bV4zd)X8z(o&YWnvq<}g>9y90nuCfuq`$RR3i3VY znS}9Sq$ELv=;QjuYFMpLHMYDqu@)cp@gZsrDeYCeow-D`e%quVVO4+vSPT$IHw~WKBxGq~J6Oa`VwYJ7_{;fr@+Shxc$c>9o3DHZ8s1-L;v;rdH z`bg&78D}`o-TEk%hL9)}NtyQm$w7n~b14&;t`nT4y|$cOem(?Vb!bSlW;;l2I#8l4 zy;r-*=O;-RE$hN}G<;lXniV%xeB9P3cX;z|H76$F38R+8+Tnq@Vpe_Dsi`-YAx6O0 z6}#U8j@V8D%Jqb{Sf;Vkqd*8XIP86o3ii96)~HH@=1s*72-eycR3BIrKl0&6&(z?+ z%vQ93N*!UtNB8r)BgZ&p%7Zjz^bH|u+80Eis!oWeNSTSUQC1i^h!=`bcO3*p z{c}T>{F@v+wA(GWlfvTfIu7X|9N8IprM5TRc2*j}wp2-1(n^lXM`f};=X#kE=2#(w zS!4D{`wG_Mm1;HOdc5p!g4N}$@zF-c&S@%7mz$bR3_BxTt-z@-&LY=|TuU8h&(f+y zlI=X$r}bWx2pmGpGBk?3V=C$L%B(NebdE19d2^utokp)rYTX1-pvC#RUY2h8Ku0 z4D!f8zOYM) z`UO*_yW?2G_*^9@eN*fGy4(pSVXNPc_qX_4SL=dM*Y@n zntx+;+p&%}gLK@8BtF0>f18FwhOVX^axDuE&qEpR^O)@K)SrRK8&yxtvb)kn)_T^9 zz4Dy$ezIG1!n`TVnZxu1(pyEP7AUgYA4BHFN8*tt^rx+tPxBd#roy8Ld3~>Z@%>46 zlq_T?H6R1>781hv9J#E>ZdtFh#FtJE!~6QU6Ov>uc$NZikh~wcVFRY>#w48Cx^K=o zy^l=vFqxblr9@Q9hUaYJ8@M;h#V=g?Hu3sOmNq3Ua-ESG6D3Qs^4#5{7GtmI5f!+L zZ*Hml`rrXD)1OJ?Upet5_7^_~G01B1>$`AzZUy&*A|la+!XX+XW<|2rx?E5@apdd z3o>X9t{yrxA%NxT3hp)z8^%J&A+L9j035EuZ_%rQ^&9Rtk{JMmFtgHpL5lrEIl^J7 zmmQKW)TqnuUuq(2U&<`+`DlrvxFar**FqpAcPR)A9DUV6;>^gUe#nFu&0VmC4-lt| z+CSd+JXp^>JTPhp0_P6FNDEe|Yiyh|(rbITrjna&s>sm$r2gytP9Lf*ei5AbDXVWJAY=CT=nCI<~K%4cBxm@{!KrJ6(~2XS$J(OBi^*5jX7rwB%Y~$b+qTx zkfECAUU^4e?K0zp}{sujftZDz-;ahRzb&R88UVxwv7 z9WE86C4cdO87@k^`5GWzEvRKSAE(f0Mi<;b%Reo&|F!fLCGW@=R12RxVL89?Hp(ok z@I|ed46y%BpeVZ}sfA^cYodMNKAwn(r_w-b{%F^}|2Y6#7^=1WCZ_W{{?4oZFwu{% zLO;F1eMdI!$pD|N_pRX-@!)qC#IaY711^Q5HKsqJwqWSOM3rl0%@Y{<0tvWVNEQrI z4UsgST>Mc~PV;5kRh#@38pXA|EK*5#+f1P>wDhdFaS*{dO}I&)>5plZMLE;5kfheGrdwYybQ z5H{YMyicV&ZTFM+++Z2I*LDutRry6Jm3GxTLm3HZ`%`FSo+u=Gw#Va5Gp)_W;Sqsm zw#UAkAF0!GiKZRuCN!$r9qount3nC^fJpHPLK_g-(2om6a6fX4ex@!IUbWRW zC3}NUx!%k2;eArRLdyCjcsO#^*`Mf%jNY2&b^TOa7G6@DSmxQJjb~TrT{HxDn@cWCkW# zs3thfV4uHS5$nYmGBwDVvdWJtR5m{-yDe{-iF#w%8hIK_-Cz9$6A}CYfbUGIkGhCz zEgJo8NMi+Q+^B7Q5Mk-r2kk==s1*F6PVweaBDByKBCPIZrA(C{PLutrXokK$FF2q! zeZ?iONAYMynzdCwO{)Hv9G{;GsU)mO!0k4$25sYQZuxpv!9B1|ej&z&v@0VE6^=D} zcW)_Bh>~WN1ZBp@EvrNVJducr5_Ss(#v#3Oz<6P@)4Ebdlx&ZDuT!W&BIokX&+7JXW2cDMP9ak z%$)l>w*GrUx274@X=E|aQ&_lXNZVxw1PevEjtd3U$*$F(z@Mh}dR&shb+pm0--zo; z|1g9FS(F-qNE{*GF@>oRt01xUxge`QL3oU1-0O1Ja0*o!r?Wo`(>NDwp0C&vou(}> z9E&&CpWxfaK7QZXl1D)B^zooAl<94ZPq3g>(mTShtOV#oXHEcKSy*QmN2dgS)qNL* z8xF(j*X+2+f}iz=!V%pf-4BVTM4V{y2mTM3=Kz=Cr+Duq({$2(?H`yQPSE|odwAV( z`sy4DzzH1rZ{Zc}w)l#Ub=TP0tV2l}||MySG&F28q=Bg8HjIs7z?XfpK=XN0m zx1ps~q!0bP{v9uaO*JRA?mO*UgSSoRfZS-++)Dr~?V)vdu;1i1^nJy2Kc}$fRNsr4 zcnLc1+(LeCrkB=fv$XYZL8`=jFSPC1eoo)%h9NS?f%Yw#4vvTZ#%b%9H0w`Z+;^H^ zJW*Cf2RMB>`hABurJeQ~94nL@!8TQ$hXg=s{D)Ua+_W*VfhYN5nZfC-hp^b6Xm=R1x+uF5S65+YwNZaiE%Es$ZnxL%c!L-1(S2nLsiSZkfP zTfZ5hIYQ9qcGLeB|4{L+m~{wG|7{%=HNh`PzT_BvbSk+Mo&g66J4~!%xLK0tfr&>< zjIqfG$NJ8Us|G(wW(^9s4CevA8XN^01F|DD&U zTa$220lOsKj&hJDY`aPxS;eRcr3U&zl1Zj;D|&5#sdZ+CMYt-2v&|#}QMbX7oLQ!` zrQAAxaYDivzJ*o47o#|0sRoKI^pOTCRyD@ixkPW-k%XwJY8_cpvO7JOlIr+Ffj~uy zdApodUl_d@xlAb_QJr=f54f70r+tdbw{{Yuq?r`Gbs z5hUea190F;gRQ?!?i58*shD9P)}^=ldP$zDm;y#ABW1p!`{~a9drL)Nl&lr zqqqgejt0+!)b*aZB3E5qJ^I;DOvCPNBW1#I7YLg56?W?;k;%ghvUsmElxg8(wX4*U z&HLtPxLGkLvw?M|F@)5blHCrgV2d_dpkRox;1J=_~I zIj2+lAk1(*KP~h8F3ZCOk<~+@z$Ip`L&ktF@5qR@L=aKJGDM4U%pqW`S*Ll=7{1cJ zyEBy#?m!+h{^nyXqF&lEuFTz%8+yS_lr4>4u(zz(agEXr;pBbZ09iD;50U$T-Cfqw zA1A>mF?fwmj3}8E*#yCKp1QPWMDKw+_Z?H+(C}_aiy^zPQ25}G=`X}Qe6V-G@cx8R z?@r!KKnHJ}x7|<@9Kl4v6g0pPYw7mDOqkfeRZFguknj4BRKQ6~f68}^?R14yI7!Au z*b>i^+ckywUWl^>aA7H%J2i1s*cU@KS?S>A&H2P9B^qK!6#|8SQX^6soca;32)g*8> z>)+x#dLehy=6@zclXV<^T92Cb1BRpG?gWZu6rYa8Z5qOT6xohMWy?99^~~gH_Oikd z$vFHK#buNX*6ItmRBy9}^6T~uvaMdl*`C_fSCzeho{}@a# z9a1+=Wko(*kutk>`@NnW`Mp>2XJQ(dV*`IWR=?G^8B5hp0K;jP3W>5@&m~@X12flm zB!xT;_E)8sOTEIPyK7jCagu8FgBsRky~ax>Z>1BR=m8D%$98A+wpj~kb8j38)Wo2n1$)o zzO+X2VI5H&aXelvHbW7QJoQ$7rK9|vPj8%};OcA&>+|Q~0%Sq!erDEk#9Tc1S~dC7 zyA(tYM7r0mp`>+Ts+uF?Oj31oCJm|>@+{pbeU=#a5WV)%D(&*G6c%&AhEYx5 zV;$n?_ffm>1n8+}6~Rwf0x=RqFr%381_a`*ZHx;N8Ch=7ILLg*RlT}6hzfH~+UMl8 z{#GkxnXs zcfN$ydh6$J8soBzsSul(&+v%UWg<#%EfECY=UytA0TBXxW_kGnrZhk zknE9pR*6&y_HRe&4p5Ns@zv688{|nTeco;cUG?zN-fdug_E6e!odc5U!tE52@L4AJ zlH$X9)=Y%63F4dagM4{oH8k}T9&W)!^iUty)B5l`rs!IlS>exTP6%tim6yxWnsGzn z#NJWRAC_Aah~3}Im5Z5F4R*%1QE zhGLSb3B{DtQZ>&;El1&H#`}jO#nClXF$%#xq1I*`j&JAy3Ai+`XGT`M9L2B#3FT1X zS$DEgz&a#ClSQ^ZbXG146_0F9+EhxHTn;I-%o2tcV~bAQ8!Tiob`18fR3UJ@AM5)E zB8cCozvhCj&EIM_#MF$`FC&4(%(kbmvU@Xz<~mY7_Mx^Qh3li+OZ}nvYKEu;2nWFe z83S+g>2#{-^Oolw3r4J}!iuv?PT|Hi8>ExEm%_mTO%7QgPK4{> zjT8C8Q<~Ac>Ild1N;m29K%E9}t9q^8Hi=0gOUQKqf?BAWSI!pnNxBhwl`j(HX%Ct& zE3sK$PDNTGUt6y+fU6Q8V*78SR=>h_v$B7(PkdjxdVcNva~<-0*EMYa>n!(ZvuCPj zGclTwzrJ{^TuB9syYHrINdw8{r<6A0W}|s|-(I3rS;baby*c~p(ATv#u}cK@>w%ZH zg|nqMI%?JsyP~^w{(OuPCd%9LkwM98Z_63;hS>604)Y1s@>r@qQA90$G64b;^S#tP zhb*RDQ2Sch=1`CL|Js^-KcZ0>kyF82^Ycl-?K1knA5rSMk0==YN?0oK@vzUdFVXLu z^^`_fCh6rVcQVN@#;vpF22J_85D6B+Xjpb1A{u@LRmKw5k9O}%Hl`a&Zg#{KUVdp{ zlG*3*nSM5BvLNTM;Znqn2OTiSd@BMkhp7HDJ*4+R>%~7rf=O z%S#B}a(jF&&lQCXUzJUeZgo2sbTZw|Y$T__mJ2EdKAuHe)raUzh#P;wszLH6UNmzR-F{m;wdL`bs< zoR67z-UwXBEOP8F5yn=M-mJHTnZ_yuz}Ia{znH;%4Uiz+Xz%=?oN_4u=vs*z>BTT} z&MoOG>I>kT>Ed?UExsH2ml}xbGAO^jTIquvNd<|^L zo3R{a;h`nC)KO9n(L*p+^~l6k+OZUJ(ApN4y>p#*wkSEyH)b$m+hV=ssiBepQDpBt zpu4OL0`+)zo;#xyzqCRbIen*sG_A9v!m_5W(se&>x8s*BjUT*5`NrkT#Z=sV^cS7_lH%OH8I2<5Lm$y{!+@P-hvP*ccmG-QCYXDTd3E*$ z;^O{9N>z1aNJ0HYCgs6^TePk1j>5u9hz1Z`1O2AVB;#{cci9Q`wRSWbC4pD9wnkAJ zHH6nS$bi>vD53gT!0t$I3y@?c2@ZEVr{7Vu-BKVZ7AdAzx|hiyDA{577BO1`brSqY zqxRCc{7{Q}`&h@)1NB9jht)v2=QIt1%2@TW(2FAqN}%dudIA}LK^<;;J*E=#S*>#h70N^k5Nr2rP6UwEg3~1TlHukqc`w1 zV($x-7w{SGPtB4yM&z_?oxlTzEqv#*ET0uThctmW+voKbZTN+R0a%U4nY=k3=~v_+ zR-0NLA;MNf4)5qi8YXYTqu0A>#jRXOKA7NSCTaq#FUWg6@aAbZAMi)K8WkceB)+97 z>(|G~aci(iZ48YZ1H`*awBFWRKvy@~8c!PSD%^C~=-Bn@ZZkM`-vw%y2Xa_B8%7bh zCIg65@Bym?ke2JyCD_~gso|P|a(mjvb&`q(^6UMehBZ$MO9CA>EWRva*3>jaud9x_ zchOqOCn{9GkQG%#?BKG=QB(b+1% z(HWb??0@an|I9mPCmU_yq4*+EeCZcE=^#2W=fZg=gp7Z@Lobx165gl&SY=*Fu5;#S zo<<7ObKzcr=PDE(fBl*~eG*So>tseRL?}ChZ=E)sE&&AG+;-sLdLf%nc!icEssUOJX%_W>Z1o>5;i(SQX$Fg8_UK z%JN?5vz2+n3Z0?-6K1{=5i3siaC8F5lB~NABNX=Ukv~!0h9M41vLR5v+n|~7Sq@Bp zvfAXFLY55sSE&3~aAv6EYKx?$uLmMSsI=PaazbHLTDJu)Pz35;XM}Yu{+m2f_(doh zzrQE{dV(g01o(^3v?xF-!YA7MjPyN`p9wAU(ShAq^l>nIb$KRhIu!y_$db+%Qpe`H z5|NMd*a-RO-y9e1KW8DTSg@ViEFMSp5OAnm>B@RM85Fb*n#ycF8B`HcJsI3p62N5K z<62<1KeGcKoAtzG$kgn6og!wv_a&`@i#|V4?C#!Ju@0Ku zmsQRXxtK6r7sw`$p2lq1Arfk7>O0A)gIuKe3JEqs(20d~>Kh>!Y@1Vdb8U-|(CR@-;}RYyKBCmZzdDy?gMW7}c}+zfm1u^csB9lqFqoFK=TE$$GPdH= zo3wdlz3>XFdqcYDObGPh6lw|N-9C#p57_~=nJWURh5@}Yv2OfTNMSv;;(gUoKK__A zcY;BA<8U&FisF5HfiCDwt@N;3IslV6JA9NRCR5Z)L~$wEY;WGYy!@S*UPn`OCbyB_ z@i-DTy6qnaQquy0Dlg;c^X4dkB3}G5;^p2yEyd#*)BSBIRBune(B~Wu+g*@68mTJq zjT#8|ZOGb32;{u0?lPR&j_|nn;x~)ZP&|U4Gm*gV+?-`9JLIuQZUWQVb9}vR1k*TppMsY#2&?G0P(T8+@>urPq({x@!KT7?8eCf{!$%>C%B z>RmMbA>R|0vs?ClG2Eys>cSBcohY`_<8V_)2Jfh^x*-pYyD&Om@g!c%cK=#;pHDQ= zoc&Y^ZQ#^vKzj_KZDr9dg+YaZ6AP4Il>N zh3jQInI)097bz7vu!e)zhNS_qk#~OW>?pD9AnTUFmc;3+^BPIJTX8xt!O#J)&)LMd zCl$hIOJOby(;w(hX7^r5I)?_j3W_sQII0+NTe9^qa*7{ATCJ*f4qkpFVbkZ+mY+x( z{A}*6K#C^Wns~E;ai%|{-+;(Scd?uzh^KRW6TZ(}>fje(6@$aJoP>G3eTKz3ccPHF zNXXg}E$!rOl-rjyIGcBMl=WWJb4J6eLEKq1;f~cA#vZ(Rqojw)|Mu3V1oh)M0~#ieDARP>?DG$aRV z^#wKvO!R0_f%fV$h)*Hv!A7?))A$sXzhmPRbqGmZ`V*Tst9yzt4Iwd2a%#}K zZ|Uxf-Fkw%-kv8^P$7j&hfVy!G&22$JOA^1L0<%g=jTM)=(J#-#D<#i=lGLncynd= zz=~ur3pE^H<4cbo%D@BW;Q{i1_zaj)si>$pl^vkcX08@gOgvoN0$iZ~j3xu>V9b0} zT>rYL?5R-y?+rKi3*~Q&iyy?3+*rp5mNo0rPYJPr&~d56sW|k9hoC0{<-- z#>K_W590lYVjdv(e-pX5UM2v@^Oq7H9`OGU;{*N8mNk0 zz&{h@=i&VakDrI{9|^vQ`zIy<{EuMRe^v=EHxJK$!2g8-|2fAX;NO+uM>$8Qi3b%#TV>FQY1p`hNhH`jdPB diff --git a/Manuscript/Ex-srDFT.bib b/Manuscript/Ex-srDFT.bib index 8f1e26f..1de5e6a 100644 --- a/Manuscript/Ex-srDFT.bib +++ b/Manuscript/Ex-srDFT.bib @@ -1,7 +1,7 @@ %% This BibTeX bibliography file was created using BibDesk. %% http://bibdesk.sourceforge.net/ -%% Created for Pierre-Francois Loos at 2019-05-29 23:30:47 +0200 +%% Created for Pierre-Francois Loos at 2019-05-30 11:26:43 +0200 %% Saved with string encoding Unicode (UTF-8) @@ -12712,9 +12712,9 @@ Year = {2018}, Bdsk-Url-1 = {https://doi.org/10.1021/acs.jctc.8b00591}} -@article{LooBogSceCafJAc-JCTC-19, +@article{LooBogSceCafJac-JCTC-19, Author = {Loos, Pierre-Fran{\c c}ois and Boggio-Pasqua, Martial and Scemama, Anthony and Caffarel, Michel and Jacquemin, Denis}, - Date-Modified = {2019-04-07 14:02:34 +0200}, + Date-Modified = {2019-05-30 11:26:43 +0200}, Doi = {10.1021/acs.jctc.8b01205}, Journal = {J. Chem. Theory Comput.}, Number = {3}, diff --git a/Manuscript/Ex-srDFT.tex b/Manuscript/Ex-srDFT.tex index 52b32ec..68dec94 100644 --- a/Manuscript/Ex-srDFT.tex +++ b/Manuscript/Ex-srDFT.tex @@ -333,11 +333,13 @@ This computationally-lighter functional will be refered to as PBE. In the present study, we compute the ground- and excited-state energies, one-electron and on-top densities with a selected CI method known as CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively). \cite{HurMalRan-JCP-73, GinSceCaf-CJC-13, GinSceCaf-JCP-15} The total energy of each state is obtained via an efficient extrapolation procedure of the sCI energies designed to reach near-FCI accuracy. \cite{QP2} These energies will be labeled exFCI in the following. +Using near-FCI excitation energies (within a given basis set) has the indisputable advantage to remove the error inherent to the WFT method. +Indeed, in the present case, the only source of error on the excitation energies is due to basis set incompleteness. We refer the interested reader to Refs.~\onlinecite{HolUmrSha-JCP-17, SceGarCafLoo-JCTC-18, LooSceBloGarCafJac-JCTC-18, SceBenJacCafLoo-JCP-18, LooBogSceCafJac-JCTC-19, QP2} for more details. The one-electron and on-top densities are computed from a very large CIPSI expansion containing several million determinants. All the RS-DFT and exFCI calculations have been performed with {\QP}. \cite{QP2} For the numerical quadratures, we employ the SG-2 grid. \cite{DasHer-JCC-17} -Except for methylene for which FCI/TZVP geometries have been taken from Ref.~\onlinecite{SheLeiVanSch-JCP-98}, the other geometries have been extracted from Refs.~\onlinecite{LooSceBloGarCafJac-JCTC-18, LooBogSceCafJAc-JCTC-19} and have been obtained at the CC3/aug-cc-pVTZ level of theory. +Except for methylene for which FCI/TZVP geometries have been taken from Ref.~\onlinecite{SheLeiVanSch-JCP-98}, the other geometries have been extracted from Refs.~\onlinecite{LooSceBloGarCafJac-JCTC-18, LooBogSceCafJac-JCTC-19} and have been obtained at the CC3/aug-cc-pVTZ level of theory. For the sake of completeness, they are also reported in the {\SI}. Frozen-core calculations are systematically performed and defined as such: a \ce{He} core is frozen from \ce{Li} to \ce{Ne}, while a \ce{Ne} core is frozen from \ce{Na} to \ce{Ar}. The frozen-core density-based correction is used consistently with the frozen-core approximation in WFT methods. @@ -365,12 +367,18 @@ We have also computed these adiabatic energies at the exFCI/AV5Z level and used \end{equation} These results are illustrated in Fig.~\ref{fig:CH2} and reported in Table \ref{tab:CH2} alongside reference values from the literature obtained with various approaches. \cite{ChiHolAdaOttUmrShaZim-JPCA-18, SheLeiVanSch-JCP-98, JenBun-JCP-88, SheLeiVanSch-JCP-98, ZimTouZhaMusUmr-JCP-09} +Figure \ref{fig:CH2} clearly shows that, for the double-$\zeta$ basis, the exFCI adiabatic energies are far from being chemically accurate with errors as high as 0.015 eV. +From triplet-$\zeta$ onward, the exFCI excitation energies are chemically-accurate though. + %%% TABLE 1 %%% + \begin{turnpage} \begin{squeezetable} \begin{table*} \caption{ - Total energies $E$ (in hartree) and adiabatic transition energies $\Ead$ (in eV) of excited states of methylene for various methods and basis sets.} + Total energies $E$ (in hartree) and adiabatic transition energies $\Ead$ (in eV) of excited states of methylene for various methods and basis sets. + The value in parenthesis is an estimate on the last digit of the extrapolation error. + The relative difference with respect to the exFCI/CBS result is reported in square brackets.} \label{tab:CH2} \begin{ruledtabular} \begin{tabular}{llddddddd} @@ -386,64 +394,64 @@ These results are illustrated in Fig.~\ref{fig:CH2} and reported in Table \ref{t & \tabc{$E$ (a.u.)} & \tabc{$\Ead$ (eV)} \\ \hline exFCI & AVDZ & -39.04846(1) - & -39.03225(1) & 0.441 - & -38.99203(1) & 1.536 - & -38.95076(1) & 2.659 \\ + & -39.03225(1) & 0.441 [+0.053] + & -38.99203(1) & 1.536 [+0.146] + & -38.95076(1) & 2.659 [+0.154] \\ & AVTZ & -39.08064(3) - & -39.06565(2) & 0.408 - & -39.02833(1) & 1.423 - & -38.98709(1) & 2.546 \\ + & -39.06565(2) & 0.408 [+0.020] + & -39.02833(1) & 1.423 [+0.034] + & -38.98709(1) & 2.546 [+0.042] \\ & AVQZ & -39.08854(1) - & -39.07402(2) & 0.395 - & -39.03711(1) & 1.399 - & -38.99607(1) & 2.516 \\ + & -39.07402(2) & 0.395 [+0.007] + & -39.03711(1) & 1.399 [+0.010] + & -38.99607(1) & 2.516 [+0.012] \\ & AV5Z & -39.09079(1) - & -39.07647(1) & 0.390 - & -39.03964(3) & 1.392 - & -38.99867(1) & 2.507 \\ - & CBS & -39.09111 - & -39.07682 & 0.388 - & -39.04000 & 1.390 - & -38.99904 & 2.504 \\ - \\ - exFCI+LDA & AVDZ & -39.07450(1) - & -39.06213(1) & 0.337 - & -39.02233(1) & 1.420 - & -38.97946(1) & 2.586 \\ - & AVTZ & -39.09099(3) - & -39.07779(2) & 0.359 - & -39.04051(1) & 1.374 - & -38.99859(1) & 2.514 \\ - & AVQZ & -39.09319(1) - & -39.07959(2) & 0.370 - & -39.04267(1) & 1.375 - & -39.00135(1) & 2.499 \\ - \\ - exFCI+PBE & AVDZ & -39.07282(1) - & -39.06150(1) & 0.308 - & -39.02181(1) & 1.388 - & -38.97873(1) & 2.560 \\ - & AVTZ & -39.08948(3) - & -39.07639(2) & 0.356 - & -39.03911(1) & 1.371 - & -38.99724(1) & 2.510 \\ - & AVQZ & -39.09247(1) - & -39.07885(2) & 0.371 - & -39.04193(1) & 1.375 - & -39.00066(1) & 2.498 \\ + & -39.07647(1) & 0.390 [+0.001] + & -39.03964(3) & 1.392 [+0.002] + & -38.99867(1) & 2.507 [+0.003] \\ + & CBS & -39.09141 + & -39.07715 & 0.388 + & -39.04034 & 1.390 + & -38.99939 & 2.504 \\ \\ exFCI+PBEot & AVDZ & -39.06924(1) - & -39.05651(1) & 0.347 - & -39.01777(1) & 1.401 - & -38.97698(1) & 2.511 \\ + & -39.05651(1) & 0.347 [-0.042] + & -39.01777(1) & 1.401 [+0.011] + & -38.97698(1) & 2.511 [+0.007] \\ & AVTZ & -39.08805(3) - & -39.07430(2) & 0.374 - & -39.03742(1) & 1.378 - & -38.99652(1) & 2.491 \\ + & -39.07430(2) & 0.374 [-0.014] + & -39.03742(1) & 1.378 [-0.012] + & -38.99652(1) & 2.491 [-0.013] \\ & AVQZ & -39.09189(1) - & -39.07795(2) & 0.379 - & -39.04124(1) & 1.378 - & -39.00044(1) & 2.489 \\ + & -39.07795(2) & 0.379 [-0.009] + & -39.04124(1) & 1.378 [-0.011] + & -39.00044(1) & 2.489 [-0.016] \\ + \\ + exFCI+PBE & AVDZ & -39.07282(1) + & -39.06150(1) & 0.308 [-0.080] + & -39.02181(1) & 1.388 [-0.002] + & -38.97873(1) & 2.560 [+0.056] \\ + & AVTZ & -39.08948(3) + & -39.07639(2) & 0.356 [-0.032] + & -39.03911(1) & 1.371 [-0.019] + & -38.99724(1) & 2.510 [+0.006] \\ + & AVQZ & -39.09247(1) + & -39.07885(2) & 0.371 [-0.017] + & -39.04193(1) & 1.375 [-0.015] + & -39.00066(1) & 2.498 [-0.006] \\ + \\ + exFCI+LDA & AVDZ & -39.07450(1) + & -39.06213(1) & 0.337 [-0.051] + & -39.02233(1) & 1.420 [+0.030] + & -38.97946(1) & 2.586 [+0.082] \\ + & AVTZ & -39.09099(3) + & -39.07779(2) & 0.359 [-0.029] + & -39.04051(1) & 1.374 [-0.016] + & -38.99859(1) & 2.514 [+0.010] \\ + & AVQZ & -39.09319(1) + & -39.07959(2) & 0.370 [-0.018] + & -39.04267(1) & 1.375 [-0.015] + & -39.00135(1) & 2.499 [-0.005] \\ \\ SHCI\fnm[1] & AVQZ & -39.08849(1) & -39.07404(1) & 0.393 @@ -473,13 +481,14 @@ These results are illustrated in Fig.~\ref{fig:CH2} and reported in Table \ref{t \fnt[5]{References \onlinecite{SheLeiVanSch-JCP-98, JenBun-JCP-88}.} \end{table*} \end{squeezetable} + \end{turnpage} %%% %%% %%% %%% FIG 1 %%% \begin{figure} \includegraphics[width=\linewidth]{CH2} \caption{Error in adiabatic excitation energies $\Ead$ (in eV) of methylene for various basis sets and methods. - The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal}). + The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal} or 0.043 eV). See Table \ref{tab:CH2} for raw data.} \label{fig:CH2} \end{figure} @@ -497,7 +506,7 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac \begin{table*} \caption{ Vertical absorption energies $\Eabs$ (in eV) of excited states of ammonia, carbon dimer, water and ethylene for various methods and basis sets. - The TBEs have been extracted from Refs.~\onlinecite{LooSceBloGarCafJac-JCTC-18, LooBogSceCafJAc-JCTC-19} on the same geometries. + The TBEs have been extracted from Refs.~\onlinecite{LooSceBloGarCafJac-JCTC-18, LooBogSceCafJac-JCTC-19} on the same geometries. See the {\SI} for raw data.} \begin{ruledtabular}{} \begin{tabular}{lllddddddddddddd} @@ -553,33 +562,6 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac & 0.11 & 0.02 & 0.00 \\ \\ -% Hydrogen chloride& ${}^1\Sigma \ra {}^1\Pi$ & CT\fnm[2] & 7.86 & -0.04 & -0.02 & 0.02 -% & 0.13 & 0.06 & 0.06 -% & 0.11 & 0.04 & 0.05 -% & 0.10 & 0.05 & 0.06 -% \\ -% \\ -% Hydrogen sulfide & $1\,^{1}A_1 \ra 1\,^{1}A_2$ & Ryd. & 6.10 & 0.00 & 0.08 & 0.05 -% & 0.15 & 0.12 & 0.07 -% & 0.14 & 0.11 & 0.07 -% & 0.14 & 0.11 & 0.07 -% \\ -% & $1\,^{1}A_1 \ra 1\,^{1}B_1$ & Ryd. & 6.29 & 0.00 & -0.05 & 0.00 -% & -0.12 & 0.01 & 0.03 -% & -0.14 & 0.00 & 0.03 -% & -0.14 & 0.01 & 0.03 -% \\ -% & $1\,^{1}A_1 \ra 1\,^{3}A_2$ & Ryd. & 5.74 & 0.01 & 0.07 & 0.05 -% & 0.18 & 0.12 & 0.08 -% & 0.20 & 0.13 & 0.08 -% & 0.19 & 0.13 & 0.08 -% \\ -% & $1\,^{1}A_1 \ra 1\,^{3}B_1$ & Ryd. & 5.94 & -0.04 & -0.05 & -0.01 -% & 0.07 & 0.02 & 0.03 -% & 0.09 & 0.03 & 0.03 -% & 0.07 & 0.04 & 0.04 -% \\ -% \\ Water & $1\,^{1}A_1 \ra 1\,^{1}B_1$ & Ryd. & 7.70 & -0.17 & -0.07 & -0.02 & 0.01 & 0.00 & 0.02 & -0.02 & -0.01 & 0.00 @@ -611,32 +593,6 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac & 0.06 & 0.03 & 0.04 \\ \\ -% Acetylene & $1\,^{1}\Sigma_{g}^{+} \ra 1\,^{1}\Sigma_{u}^{-}$ & Val. & 7.10 & 0.10 & 0.00 -% & 0.07 & 0.00 -% & 0.11 & 0.00 -% & 0.11 & 0.00 -% \\ -% & $1\,^{1}\Sigma_{g}^{+} \ra 1\,^{1}\Delta_{u}$ & Val. & 7.44 & 0.07 & 0.00 -% & 0.04 & -0.01 -% & 0.12 & 0.02 -% & 0.11 & 0.02 -% \\ -% & $1\,^{1}\Sigma_{g}^{+} \ra 1\,^{3}\Sigma_{u}^{+}$ & Val. & 5.56 & -0.06 & -0.03 -% & 0.07 & 0.02 -% & 0.04 & 0.00 -% & 0.02 & 0.00 -% \\ -% & $1\,^{1}\Sigma_{g}^{+} \ra 1\,^{3}\Delta_{u}$ & Val. & 6.40 & 0.06 & 0.00 -% & 0.10 & 0.02 -% & 0.14 & 0.03 -% & 0.12 & 0.03 -% \\ -% & $1\,^{1}\Sigma_{g}^{+} \ra 1\,^{3}\Sigma_{u}^{-}$ & Val. & 7.09 & 0.05 & -0.01 -% & 0.08 & 0.00 -% & 0.16 & 0.04 -% & 0.14 & 0.04 -% \\ -% \\ Ethylene & $1\,^{1}A_{1g} \ra 1\,^{1}B_{3u}$ & Ryd. & 7.43 & -0.12 & -0.04 & & -0.05 & -0.01 & & -0.04 & -0.01 & @@ -668,55 +624,9 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac & 0.05 & 0.04 & \\ \\ -% Formaldehyde& $1\,^{1}A_{1} \ra 1\,^{1}A_{2}$ & Val. & 3.97 & 0.02 & 0.01 & -% & 0.05 & 0.02 & -% & 0.03 & 0.02 & -% & 0.02 & 0.01 & -% \\ -% & $1\,^{1}A_{1} \ra 1\,^{1}B_{2}$ & Ryd. & 7.30 & -0.19 & -0.07 & -% & 0.00 & 0.00 & -% & -0.02 & 0.00 & -% & -0.04 & 0.00 & -% \\ -% & $1\,^{1}A_{1} \ra 2\,^{1}B_{2}$ & Ryd. & 8.14 & -0.10 & -0.01 & -% & 0.09 & 0.07 & -% & 0.08 & 0.06 & -% & 0.05 & 0.06 & -% \\ -% & $1\,^{1}A_{1} \ra 2\,^{1}A_{1}$ & Ryd. & 8.27 & -0.15 & -0.04 & -% & 0.03 & 0.04 & -% & 0.02 & 0.03 & -% & 0.00 & 0.03 & -% \\ -% & $1\,^{1}A_{1} \ra 1\,^{3}A_{2}$ & Val. & 3.58 & 0.00 & 0.00 & -% & 0.09 & 0.05 & -% & 0.11 & 0.06 & -% & 0.07 & 0.04 & -% \\ -% & $1\,^{1}A_{1} \ra 1\,^{3}A_{1}$ & Val. & 6.07 & 0.03 & 0.01 & -% & 0.13 & 0.04 & -% & 0.15 & 0.05 & -% & 0.11 & 0.04 & -% \\ -% & $1\,^{1}A_{1} \ra 1\,^{3}B_{2}$ & Ryd. & 7.14 & -0.19 & -0.08 & -% & 0.01 & 0.01 & -% & 0.02 & 0.01 & -% & -0.01 & 0.00 & -% \\ -% & $1\,^{1}A_{1} \ra 2\,^{3}B_{2}$ & Ryd. & 7.96 & -0.09 & -0.02 & -% & 0.13 & 0.08 & -% & 0.14 & 0.08 & -% & 0.10 & 0.07 & -% \\ -% & $1\,^{1}A_{1} \ra 1\,^{3}A_{1}$ & Ryd. & 8.15 & -0.14 & -0.05 & -% & 0.07 & 0.05 & -% & 0.07 & 0.04 & -% & 0.04 & 0.04 & -% \\ \end{tabular} \end{ruledtabular} \fnt[1]{Doubly-excited states of $(\pi,\pi) \ra (\si,\si)$ character.} -% \fnt[2]{CT stands for charge transfer.} \end{table*} \end{squeezetable} %%% %%% %%% @@ -725,7 +635,7 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac \begin{figure} \includegraphics[width=\linewidth]{H2O} \caption{Error in vertical excitation energies (in eV) of water for various basis sets and methods. - The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal}). + The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal} or 0.043 eV). See the {\SI} for raw data.} \label{fig:H2O} \end{figure} @@ -735,7 +645,7 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac \begin{figure} \includegraphics[width=\linewidth]{NH3} \caption{Error in vertical excitation energies (in eV) of ammonia for various basis sets and methods. - The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal}). + The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal} or 0.043 eV). See the {\SI} for raw data.} \label{fig:NH3} \end{figure} @@ -746,13 +656,13 @@ Water \cite{CaiTozRei-JCP-00, RubSerMer-JCP-08, LiPal-JCP-11, LooSceBloGarCafJac \label{sec:C2} %======================= It is also interesting to study doubly-excited states. \cite{AbrShe-JCP-04, AbrShe-CPL-05, Var-JCP-08, PurZhaKra-JCP-09, AngCimPas-MP-12, BooCleThoAla-JCP-11, Sha-JCP-15, SokCha-JCP-16, VarRoc-PTRSMPES-18} -In the carbon dimer, these valence states are of $(\pi,\pi) \ra (\si,\si)$ character and they have been recently studied with state-of-the-art methods. \cite{LooBogSceCafJAc-JCTC-19} +In the carbon dimer, these valence states are of $(\pi,\pi) \ra (\si,\si)$ character and they have been recently studied with state-of-the-art methods. \cite{LooBogSceCafJac-JCTC-19} %%% FIG 4 %%% \begin{figure} \includegraphics[width=\linewidth]{C2} \caption{Error in vertical excitation energies $\Eabs$ (in eV) for two doubly-excited states of the carbon dimer for various basis sets and methods. - The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal}). + The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal} or 0.043 eV). See the {\SI} for raw data.} \label{fig:C2} \end{figure} @@ -766,26 +676,14 @@ In the carbon dimer, these valence states are of $(\pi,\pi) \ra (\si,\si)$ chara Ethylene is an interesting molecules as it contains both valence and Rydberg excited states. \cite{SerMarNebLinRoo-JCP-93, WatGwaBar-JCP-96, WibOliTru-JPCA-02, BarPaiLis-JCP-04, Ang-JCC-08, SchSilSauThi-JCP-08, SilSchSauThi-JCP-10, SilSauSchThi-MP-10, Ang-IJQC-10, DadSmaBooAlaFil-JCTC-12, FelPetDav-JCP-14, ChiHolAdaOttUmrShaZim-JPCA-18} -%\begin{figure} -% \includegraphics[width=\linewidth]{C2H2} -% \caption{Error in vertical excitation energies (in eV) of acetylene for various basis sets and methods.} -% \label{fig:C2H2} -%\end{figure} - \begin{figure} \includegraphics[width=\linewidth]{C2H4} \caption{Error in vertical excitation energies $\Eabs$ (in eV) of ethylene for various basis sets and methods. - The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal}). + The green region corresponds to chemical accuracy (i.e., error below 1 {\kcal} or 0.043 eV). See the {\SI} for raw data.} \label{fig:C2H4} \end{figure} -%\begin{figure} -% \includegraphics[width=\linewidth]{CH2O} -% \caption{Error in vertical excitation energies $\Eabs$ (in eV) of formaldehyde for various basis sets and methods.} -% \label{fig:CH2O} -%\end{figure} - %%%%%%%%%%%%%%%%%%%%%%%% \section{Conclusion} \label{sec:ccl} @@ -796,7 +694,7 @@ We are currently investigating the performance of the present basis set correcti %%%%%%%%%%%%%%%%%%%%%%%% \section*{Supporting Information Available} %%%%%%%%%%%%%%%%%%%%%%%% -See {\SI} for geometries and additional information (including total energies). +See {\SI} for geometries and additional information (including energetic correction of the various functionals). %%%%%%%%%%%%%%%%%%%%%%%% \begin{acknowledgements}