From c2cb7d7a27716c79de3700d857fe33a841e3aded Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Wed, 27 Mar 2024 17:20:04 +0100 Subject: [PATCH] CuCl OK --- Data/cucl.dat | 15 ++++++ Data/cucl.plt | 30 ++++++++++++ Manuscript/cucl.pdf | Bin 0 -> 9957 bytes Manuscript/stochastic_triples.bib | 4 +- Manuscript/stochastic_triples.tex | 75 ++++++++++++++++-------------- triples.org | 36 ++++++++------ 6 files changed, 110 insertions(+), 50 deletions(-) create mode 100644 Data/cucl.dat create mode 100644 Data/cucl.plt create mode 100644 Manuscript/cucl.pdf diff --git a/Data/cucl.dat b/Data/cucl.dat new file mode 100644 index 0000000..ceda691 --- /dev/null +++ b/Data/cucl.dat @@ -0,0 +1,15 @@ +# R CCSD(T) Stochastic CCSD(T) +#--------------------------------------------------------- + 1.55 -2099.590506349950 -2099.58890791 1.3630E-03 + 1.65 -2099.671184187604 -2099.67175286 1.5710E-03 + 1.75 -2099.720045862965 -2099.71877319 1.3199E-03 + 1.85 -2099.747811193906 -2099.74897746 1.4668E-03 + 1.95 -2099.761752030920 -2099.76232971 1.6105E-03 + 2.05 -2099.766727898670 -2099.76565267 1.5202E-03 + 2.15 -2099.765956694308 -2099.76485609 1.7470E-03 + 2.25 -2099.761562105614 -2099.76237391 1.7474E-03 + 2.35 -2099.754944906474 -2099.75681975 1.9951E-03 + 2.45 -2099.747028328725 -2099.74813718 2.4288E-03 + 2.55 -2099.738443175793 -2099.74031232 2.4057E-03 + 2.65 -2099.729597826175 -2099.72866832 1.6894E-03 + diff --git a/Data/cucl.plt b/Data/cucl.plt new file mode 100644 index 0000000..26055d3 --- /dev/null +++ b/Data/cucl.plt @@ -0,0 +1,30 @@ +#!/usr/bin/env gnuplot + +reset +set grid +set xlabel "Cu-Cl distance (bohr)" +set ylabel "CCSD(T) energy (au)" +set format y "%.2f" + +a0 = 1.8897161646321 +E(r) = De * (1-exp(-a*(r-re)))**2 + E0 + +a = 0.84615 # +/- 0.03216 (3.8%) +re = 3.92539 # +/- 0.01058 (0.2696%) +De = 0.101589 # +/- 0.00932 (9.174%) +E0 = -2099.77 # +/- 0.0008014 (3.817e-05%) + +set xrange [2.7:5.2] +fit E(x) 'cucl.dat' using ($1*a0):3:4 via a, re, De, E0 + +set xrange [3:5.2] +set term pdfcairo enhanced font "Times,14" linewidth 2 rounded size 5.0in, 3.0in +set output 'cucl.pdf' +set pointsize 0.5 +plot \ +'cucl.dat' using ($1*a0):2 pointtype 7 lt 4 title "Full (T)", \ +E(x) title "" lt 3, \ +'cucl.dat' using ($1*a0):3:4 w err pt 0 lt 1 title "1% (T)" + + + diff --git a/Manuscript/cucl.pdf b/Manuscript/cucl.pdf new file mode 100644 index 0000000000000000000000000000000000000000..66edd9e650a3525b5be6dd8751ac126cce54fd53 GIT binary patch literal 9957 zcma)i1zc23_dg)rjfAkEbhFfw(jeU_y>z)OvZQni64FS6lA=;dcXvt%(%lX6Tk!FH zpZ9%!|Ihz-@7_Ci&YYPubIzT)XTN9Y)ud(FxjFbS=__~2t}uB4TmWZtJ4_J~0H+ez z$r@?{;DeJiFaZDnr>w056byl<4rWlWG}yw~5{xM(hUp50fXy5+VdZ33kC3HA)k;|p{HhGN{Fso1!y$9&N^lY0JOaWknYjS z<1k6o2eSuQ6-uh%eN6#;bR&{jR7ReDM#~=}-4L2%B?&T^nxgyQ%ErcSDlocyy*!7* zTU$M(bre8Pum|`+jvq}*%Fnc!&Lf32)Q?Q=MJch$j=WL>kZ%)4w?hre3J%5z5Gqs} zLV~hN)Xul5dTJI#`!>gfx$lo+ql;rq;QZiB2b$rbmUS|r>+uT`Lr6cy|nDyb+Vyc+rY8@;p~%U((>5W z>=+5(2}z&cD@vgAv%S(a!@5m-ihhP8mia8}_V)qT#7mF90Xr~?tk^YY0XE0wbSR2; zljmYFMKu|T%0#+h`ktGn?iv=3sXZi41?NI9pY0f91qtivOO;m)m-Oy?=`)SlJ~v-< z!)R~kn*LrN#p1MPKH5g2lK4(&Dr2?=VJ(Jrwhj^uv$6DHpHt8vo>%5in6XR!zA%~U zI6JU!^-%Yr@tp{sZO8OYaEQsF7fAan(hP-Evn-yk@69n?vRHr=SqjQ-ZDd8mQ-yHO z){CWvp&yyQj(yLLeZ_sh(_HT+shXTk>@{rTnVmHr-0ZeipLs=i+7HqTgUIC2a4qc9 zh6Ep8y_}*kT38!BKQL*`-niYxo7fe*hFl}Bo`u$3b-p?^I~y!ln`|RV?wK#3`1mAb z!GbzdETF$}yL8$|I72p~`UOba!X61Vo8@wy#nn~CCwe~$^fRXzFN@9tuk_+-8SB^CqSB!C=aihZNC;}5K5j7k~OoJ29d>6v5&-Vg2Fqfg+WB|{lvo! zT-V8Zw1tN&@3Z?6hWT--Dat}{QV8V6v9;cLgeB|s;vD*Qi7t;(eDoTt6R{V>^D^YO zRQ{CX-}Pv1(aN0t3T+j^wjTiQT{bzmq;b;P>;y#l)&hBnh%*oUMo4#d8-JvX+BQJ zsth215p2)ETZfJYQi~|=yG4lTo3kWvY4@E^&B|HhvBbf|$k1Om7oAJd?mi#bmlQ>O zqYET8hPu#cikxVxO~@NAPM>0wT8QJ;#nvP88CT_H!&ITZ^p{eA50W#wwhvsFoHC#o z5`zyWWYFESe$l$HeY|Ee6{PKc$a%ufFN{ypPH8{-gi2HmNlP!^&)-K3VS%Fnrm03e ze>SukE=kW%|FF|kQ{&+iR^{99)kjiN(b-Qn@NFjnf#DgJc})~KQ~_CjX204|Yn}w; z!!o<5boq-DT8Y?YuEYC(gk!ot39hpo5IK8qrncBJWiOJ+M}ope;^frLqxyUikq2YE zC%nC9dQfL$1hQ0Z#XMf~zB?W$+c2oDuRj30Jk@Fc2>cEv2_Hyd`LQ3tp(YZf*3>2F z;huO#Wuq2QX_U7Y@ub%cOUc1XW`v`qm{RZi<;ENL@*}68E6v{q79Fita`By!&_J^< zl=H>#d7FEFmTZIMX-M3d1#KuiY&Y#*e446Qk3)TC7^Q(0h)ei@G#<^KP!)9-dwf5n za5H?M3Eww5YD`E|P#+gjq9qVnda=7mR(W|}ast$8H=m z+Y&oVXW{(MkYeD1wUqI&ifat#C#t#v>MpvX*s-PfzJ#PPpdJnQgGVoU$0fi7sqO@F zv7W}T2yMw5Axw_BJtrVq8xXb?s*k0OCgQLXRo`bO(8OmIvFa{CDy3&jR`BukZDA`& z?THoYHx}E!9?k$X0$wRVqqp{U+{N1$GlV~WJ9ysu3sGfZOVf?5LpQXb;3ov(aBu?L zJ6x;eNZIg^FVf@tFS`D@?m!zVP(eSuyoFI5Z$q zZk17hysN$?Nq#v9-2EgjFDmm4;m1T@>d@SH1|kEfZQv>MflS;GiP0Ow?vQWWoUUtt zl_Zs)IkBzl3`1X@l{F3ZKrz#Z=KahY9~d8mO-gOYM$grbtEx;=0(L!#wN9Qbusq#< znbrF)eb%mjkI%j(-v{HZFXADVCr3cQ_ANPdT2p=m2R&*sWo zsqg|x%yd=#26H^Gat*dpB^V9mlHvJ-U{Z#FGFXxN?pc60oqBLN8R)~v*kH}<_cuGY zbz`RH^{vlO*D}6v8FPM5MeUiV&l<!Jqe2ln5-UlJcQ zF|bgKReq2#C4DpO@q2Z0vBSs+uPniLiMom+v>PFwv9*gRPKMV}MvcG4*wEggAX41v;@x9&B5Plk$QCKw(_e(ODtzHIW=C z+TN}@MkKr*l@SNHmDmTQOZEQ+WK{X?V9x%HKRg-jY6~zBRn#o%Us%kZLRM#U?n!Xg z1}~hflG~S-t$qc;R3{iwY}p%HEJ@ywtO+0DT_UwLjq41yJwB;&t`^6=Me~zojKu^y zS^nic!k_N^&O2xGE`__BJpA0;cMNw~;J@6`dq-3T29?u`Wtxszq#@rSGcfI{5Be{xcAE@|*fTMMv;oa7zx zuc2yyU7g(^7GPJvo!T`1HZlQt{=d9?!GCq@fAj9|6szfG4!zgHT|c>{{=?(XNJ4{+ z$$syo@3*Pkjze)sOeSFk;RAxK4E%zEtiE#B_)(g2)hy1>lx=xBmb32SlzZ*rP*P%} zK0uXZoO;UjgjP;+*H01hIJ4%PGmf`y=h58RtVqmiO6+li=h}GR^;)O}zNuieD!`-t4^ka4y)O6#=86t(<>N0o2EN3yxn{I&~jgpSItb2g-O82@xIKd*B%c@1M z-ukDW4>+!@BXOe$Bk5T10@mju(j-sCD2}1XuJ@&>Yd!G!s*daGKt}a^uW;Ms_;mXZ zmx|$>{19B)5a)LfDTFU;g}f*RM);5(&DGb5QikGnUkd^<8| zLfbY1x{NLWPdFQ0S{E%OT!b!@34j@%D@~A%g0z71^QnkD0>vBjxTiewh#vqr#D0=4 z#Kef~OsMo(#0mZoHAEH(pbcVa7S5MI+#!r;e~uxX%Vu&HQfxoh1(6yA)MjNDtXWhO zKVegJ?7(je)HQew$aqcY=?IK13PT7g!94P+jL@2xlc2v9$#p+}xiRo77-m;$f^ znz@i=o8%v&`K1BW2x zkK|R&FQ*WN7rCcUzAe6(df>H)0V6zZ*4*da4z*s;7ewc5QQYUbrZV;ySY(G$`XcT% z5uD?T1*k3nxUImM{WG)QuJL5~>{(A6>=vU;JvQQiV5kn!4_zf}hlg?N=4W($t zAGGw#$PhYa`ocGRaj45$8nPqHyc#MzPZTNQ_F_l7$(OH|l^~>`v}?a~DGor=8^!kY zjWrFhO^|`aRqtwa)IBCjPS;Vy`%DVWI#_(L?!))R9H7{5aB#rFO__4Zi3KS*S+f0Z5&nM=an_ zgY|8!l3%q~5SmX{X?CED`-T3qO_=06o^F+{LiJ+L5@)_@QY{^xjErSW7QV%xD&{N8 zd94b)*f1FiMCVxk|&N^fM*?lf&=SJSV^&HC%M)pZM%UL5K zAwEgfnylkP1T>nMeTTc0m(bVetdoaQ4 z3KI9T#Z$!#ce4H!a^LUq)AD#5qBj^7U8alf*OvUVib4=44$d`{EBaUsZj5>jdciVe3eR5-!p4hFEiI1`f|I;>=XCSO*@iph&RYlIt@t$-rLF=EMaTnsi|C_U7a!76&)~e8m~wKig_` z-AR^?u$X0vX-5iJR6Vgqo$g&rZM={;opPPZ$*ipk?@WV%`jb$xc8hq-<5p7L2ng)Q zBrOM?$4yv)q26xHeU+i)^WzOsotAZ7W*41a4^{0yhvgqPhxG~RhycceIvF~KIfkd_ z*xbW4822z-U(XT>@kKrBkyq@?O9SRo>CJ}~xgHOQ$sdg&y<09kig)IeSx|~Rg+2K8 z`q^u_&sXxkyB;5&kr0q{BL`nfdyMqD&q)yxY}fm&K?mQ>mDGfOc)zZ z`&6y&=V+`>BBoJ&ZgeTBuh~O7n&C~BelGRSEzIt6r}XpcgLH}M+fzzlKFty8+@DG8p<$wH(sWH(o#=WnEUwh|i~lh%q+ zeeID`qWyF|zIMlsn@fvrHN5(zn`V4mw(dBL>*D#n8uu}fIgK=?p}f7_++p>)U6z%M z&E#NMqm{mQ!tmiMQ&8vd=p<ol~sbT0&N0kNT$I@i_DlA~fgoSVg`;>zc9+Bzi z4k*fBeq6aD$Q}Gu`1zC z%`^*x%T5r^+Ds8iV`AQY{WytQmIGW}PfEETVo2#^af+>@>3wwCUc8zY4-pN5qq?2n zeVHdu6upAkgyd^CNicTwYcG=88I|y#yBb=|xm_-jAFULhau_uvuLvAkSt2v$K4Ag3 zzc)rPJ|$ieLt%PeEpAAS?)uWapN+~Rf`)M^+|w8~8&KgawjuB|bDPbNJjLoJDi1lYR&ek3E&OXd$ssLbPR~L--rdYnx2n2e-0w?O3<^ zh;ZOEk{>-wl>8MkQr$DXC2e)@#-ItzkgItOKdcTeq&t&{u1Q49h+n z?8Q9RSkn_eS5;+I)G%mm&5aRWuK1YnljkA}vp%x?L*u>&c7->xnwf!Vr`#5qE2aYT z`^foXsOcB$6wI7qJXpro?8_B{FAtH784L|-^F~7|oIyJR^10zL4RY9VWsRLNpYb$6 zsUP3ka0H-BMY06dJ@~ejU*scVi(oUc@EnV!Xi=rFRP-m-*i2 z<7JFDSF;7vyroEup_U(f%OIWfiR{LZjlW0i(^S<2D^q=6IB1hRmHOEgYl1&V2IA!2 zY^lwjSxj{jo>A;cu$$L7k3W4g=qfVk60e+^(^-uDY9#NK|7)N3X~3n${C%9+SQ=}W zZ>7b%bT=0jH#6xZ91UV6p}XA-;@KQQR=eBJsa{MXes7tJr?ZzZd!`|aIdiBC6(Y|I zi^j?FkVA|)QqN2f zZD6Bw3eu8y5B8R&agKl*~l7eB9e5IaFun4 zq{*i@3Uz(Tbiu9Odo$#pdrao{nIo<&A(vyb(%$k6{mVXOUbUBfPujO&fd)J8kdE#l z%R-vOI#Vi7?`I{ObS>DsDOz8(N(mAj2aJJ={G#6Z)yzxvb91)!DQXb;XxF2H$PQnA z#+}Z+*MKDgE^9WIxz{8H*8}KWgbH!XuyHqHeyJxH`ay$)3L-0>*3?;eyIbz{jFZ(> ziCb$75yz!u-Sj(<>kjnqf=E_f^gntJ#AEz^f`!)9ebP%X%e{jFWD&?Je74C2_9dOq z+vOf#PGT1H5p0M{yNodz`YPnPL-0%M%6s(e02-tTSC{VUYEC) z#~>s$?<5XLN;c~d&`Z=%anhs;oMsEDcoQ`82Lv>r+WcBC)NuP0I3oO=!@I~D`-A4? zW9D#jsG-#{;wcz$MGSH-#sl+5=Y61TbYofdAth27|1)J3wew~%OAIE5)#j`se`(S| zlhe$zVjp^vC~5X`;{HTM{)i_wb?HZ=rD%Ekh33ng0&1x-Idr;P(Sq+A!VEHU`bK97 z=~Q{+FEuo?`Y!YY6BmcW4XA22tpaT@iw4NNNmJ)?&1z}hCtLb{CiswdXxKXf-1w;! zOnEgD@{9{ZyKG|JYHZG>vF`skUg>F|RCR_sX{x75D;XT6%5 zqZ%3-F?j_ARnas*Qj3%om6eCxK529|1}YlIb2hQPzz{>se4jfHyY@0?p?z%AlHx*Z z6y7S0H>{WrP$BebbNz{efPkbSS4oT?@Zu?Qo2dv4O}+`4i$H0_)`2)|V#9&73l#() zTlT*e$h%6lj0vP3!cy^z&%#SbDttvdBjG||G-7+Swt!iKS?0gB&_UVb(_Xni3qx)| zG6{5Uve+lSM)-yO3+WvqvLF}cFI?C@;q`-SqH8R3{7__M{Fa0LM@$~At)6N$Dk0dG z*|gqx1EE;ZOXp)!j$eB2EP4+nTijekc|#}{WNR>-LJWUl8H=5VB-X#mUQv69N}`6* zdiz+Qu!R?fT+oW0jPco&zt6vx3vabG;5NR{&~(c%`I1i#ugP^^q;6mMT?1x$^VWWc zbB!lDZtJ&w>N>o}AUtEuD!nKjdEJ?vvsu}+y>hXw?vO92n5Uon*jE~(G z&J`q%pBdZ6W9OC@9P2bXfC5Rhn^&9zv?YOyp~E5@E-$5yP%b4?CP$2PM6~o#YFg9% zznlPB)KcQ8+fXgv=!+TKiX8`3c3>&PadRpwpKGK*I=ZE$0#wvC*^nNjPmhHIl23HR zLNhGCg>Gm1avKc~cYYZoov7?WJ^NCw42EqGKy*KuSotd%{?JQlZWAD(e=^$dc=5K< z!AIM72}8f~tCoJ}eo+$7haxk<( zMKP1I6(Y)VbMfV356nlj7Dqo1KUvBh zjS9!7c+X2HXt$BbS_RP|mUDLJ$t2RMco|#-s8ZP5-v_Fp`t8)AN@)`VrBa%}gEo3c2bCJ9xkl=4CFJHF4>A`MMZ=3c_O*h;2|c9 znIm+pfLn|7fs+tVoYP0wQ(z%k>a5gQ%+2JUz=|x$szMG)elRLGn;t6jLVY6zr!#5B z01YJkk^j^T*nC0JcPNh>NCb_7DA-(?^O{6Wy zGp;@zMJFzop_sHnrZbLa_Z$y~Qc`Vpn`faKiC~>C6BaG1dB_tFSwGYzC{V-DsiRUE!f9b3u)Z1Ga;AOAhJZxxwjuwCTH>Kgf|&OM_Jve_h%}SB&?a; zA{2w1qrg*Id%BR6)O1+TZeQ5+U;{LlaH80vPg@cNGFaF*V*ZH;x>}*-E~jipIcK?D&5R(LaytcM9af z58Cs^q#2b$<~9`|jL^y_!@W@-6+@^m#(rijv*1_Kzfe?#l3bFg~^$f)3@_4Z5U%^j%+i?C@6ls&M5<*5-_gOmq3n| zJAk6-%-~Z}s^EE(<<%`%0Q9TMC?3u-SxRK9gwUpWZPH6mN3=M79%o9jwJY$P?97jd zK{vV&BWJj=+v_&(3>33Z$Ozru(OAy8p-If&q_JYNtqu5AON{lj=X+42sgtq1prTt_ zNcX^-m(vM4Y5VH$luXkc(??N~qUh3XJf`)*3bEd_tVw&B(6ieUr7T%-j9K&1c0Q~e z)AKEKSu0&Lw_r7uB8AVIUodwOXtz9>{-eo0&K=FCQ)n|;t5dSNZ%X&5fx@`oIIu4S zNLj*div#WdDvO%I16_zFUmNi*BfR~nbc^qSIgtB?5!bu6+Hl!^)wk{}ralDUkde=E zlYBREezd;4@P)aqr!S(8Uz^N7O?_{`zibTz{mazccNp^@y!JOrdlZ>%N+V ztHm8M&MyeO&))U(F3%3+7Pu!!nYqY=ZLMwK;5odRdy*y;?5G2Pe{S!QW_U*aPv&l< za7Y>q1l*z4e@C(JQ2T#s<#!AJfmPp=B+XpG_u&8kS3R~6SE!VY83X|2xv!~Y_E(;p z=Qr~IpUMZ^o9FkmlW~G?hOu?B25{=yI!QRW+WsNjH>7K83AJ&BhX;80?q=%mQxJHE z&cnk6@bvEbO}%^i&BX`T_W$Mj)84;RetrRf0MA_%gOBe{SO5EYR|b!;{69*$0KB}s z_wxUq{<}v!0(Z+3JkW8U{^{o(<$gSZp!Dj>4i#-J?=^T==b|c9On}FEL3A@o@3*ashb+ zxcT`6xIwI3+>BgY%zqd9@8Hy*c#)JDd~P^f-*^0e2!BP3{vBh2LBLj+@ae>h37=Je z9RPklJ|G{!3h);V1iUlf-2rg=L*o&IbNn|=faf24TwL7$3m+G}tN+=aprGJCXu!KU z{@;9DAh;p_t1dV9Kl%9JE5(2F@pHk0U;kMb#LM#!J`f)e{!;zngF?*UvmbK56HU|B y6MUzuyG267*%@whxZuCOcnVHdaNFJ){{B^Rg_=R2_u2*Vfw(d0>19-9G5-e_luy+F literal 0 HcmV?d00001 diff --git a/Manuscript/stochastic_triples.bib b/Manuscript/stochastic_triples.bib index 86d7d53..d2819b4 100644 --- a/Manuscript/stochastic_triples.bib +++ b/Manuscript/stochastic_triples.bib @@ -132,7 +132,7 @@ i@article{watson_2016, month = sep, issn = {1870-249X}, publisher = {Sociedad Qu{\'{\i}}mica de M{\'{e}}xico A.C.}, - url = {https://www.scielo.org.mx/scielo.php?script=sci_arttext&pid=S1870-249X2012000300014} + url = {https://www.scielo.org.mx/scielo.php?script=sci_arttext\&pid=S1870-249X2012000300014} } @article{du_2020, @@ -170,7 +170,7 @@ volume = {118}, number = {24}, pages = {e1797915}, year = {2020}, -publisher = {Taylor & Francis}, +publisher = {Taylor and Francis}, doi = {10.1080/00268976.2020.1797915}, URL = {https://doi.org/10.1080/00268976.2020.1797915}, eprint = {https://doi.org/10.1080/00268976.2020.1797915} diff --git a/Manuscript/stochastic_triples.tex b/Manuscript/stochastic_triples.tex index b2c69bf..94981dc 100644 --- a/Manuscript/stochastic_triples.tex +++ b/Manuscript/stochastic_triples.tex @@ -1,5 +1,6 @@ \documentclass[aip,jcp,reprint,noshowkeys,superscriptaddress]{revtex4-1} -\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,longtable,wrapfig,bbold,siunitx,xspace} +%\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,longtable,wrapfig,bbold,siunitx,xspace} +\usepackage{graphicx,xcolor,physics,siunitx,xspace} \usepackage[version=4]{mhchem} \usepackage[utf8]{inputenc} @@ -13,6 +14,9 @@ urlcolor={red!80!black} } +\DeclareSIUnit{\bohr}{\text{\ensuremath{a_{0}}}} +\DeclareSIUnit{\hartree}{\text{\ensuremath{E_{\textup{h}}}}} + \newcommand{\mc}{\multicolumn} \newcommand{\fnm}{\footnotemark} \newcommand{\fnt}{\footnotetext} @@ -69,7 +73,7 @@ \begin{document} -\title{Stochastically accelerated perturbative triples in coupled cluster calculations} +\title{Stochastically accelerated perturbative triples correction in coupled cluster calculations} % Alphabetic order \author{Yann \surname{Damour}} @@ -104,10 +108,10 @@ This work opens up new avenues for efficient and accurate computations, enabling investigations of complex molecular systems that were previously computationally prohibitive. \bigskip -\begin{center} +%\begin{center} % \boxed{\includegraphics[width=0.5\linewidth]{TOC}} -\end{center} -\bigskip +%\end{center} +%\bigskip \end{abstract} \maketitle @@ -203,41 +207,44 @@ accelerators.\cite{ma_2011,haidar_2015,dinapoli_2014,springer_2018} % - Benzene TZ % - Streptocyanine QZ: Small molecule in a large basis set % - Caffeine def2-svp: Large molecule in a small basis set -% - Vibrational frequency of F2/cc-pvqz +% - Vibrational frequency of CuCl/cc-pvqz %b. Discussion of the obtained results, comparing against other methods % - Measure flops and compare to the peak %c. Analysis of the algorithm's accuracy, efficiency, and scalability %d. Discussion of any observed limitations or challenges -\subsection{Vibrational frequency of \ce{F2}} +\subsection{Vibrational frequency of copper chloride} -In this example, we compute the vibrational frequency of \ce{F2} by -computing the potential energy curve, and fitting it with a Morse -potential +Our methodology proves especially advantageous for scenarios requiring the aggregation of numerous CCSD(T) energies, such as neural network training or the exploration of potential energy surfaces. +In this section, we discuss the application of our novel algorithm within the context of computing vibrational frequencies, specifically through the example of copper chloride (\ce{CuCl}). +A demonstrative application presented here involves the determination of the equilibrium bond length and the computation of the vibrational frequency of \ce{CuCl} using the CCSD(T)/cc-pVQZ level of theory. +The procedure involves determining the CCSD(T) potential energy curve for \ce{CuCl}, followed by its analytical representation through a Morse potential fitting: \begin{equation} - E(r) = D_e \left( 1 - e^{-a (r - r_e)} \right)^2 + E(r_e) +E(r) = D_e \left( 1 - e^{-a (r - r_e)} \right)^2 + E_0 \end{equation} -where $E(r)$ is the energy at distance $r$, $D_e$ is the well depth, -$r_e$ is the equilibrium bond distance, and $a$ is a parameter -controlling the width of the potential well. -The vibrational frequency $\nu$ is calculated as +where $E(r)$ represents the energy at a bond length $r$, $D_e$ the depth of the potential well, $r_e$ the equilibrium bond length, $a$ the parameter defining the potential well's width, and $E_0$ the energy at the equilibrium bond length. The vibrational frequency, $\nu$, is derived as follows: \begin{equation} - \nu = \frac{1}{2 \pi c} \sqrt{\frac{2D_e a^2}{\mu} +\nu = \frac{1}{2 \pi c} \sqrt{\frac{2D_e a^2}{\mu}} \end{equation} -where $\mu$ is the mass of the Fluorine atom, and $c$ is the speed of -light in cm/s. +with $\mu$ denoting the reduced mass of the \ce{CuCl} molecule, and $c$ the speed of light. -% CCSD -%a = 2.2936 +/- 0.006318 (0.2755%) -%De = 0.125888 +/- 0.0005213 (0.4141%) -%re = 1.3893 +/- 0.0003428 (0.02468%) -%E0 = -199.338 +/- 6.422e-05 (3.222e-05%) +\begin{figure} +\includegraphics[width=\columnwidth]{cucl.pdf} +\caption{\label{fig:cucl} CCSD(T) energies of CuCl obtained with the exact CCSD(T) algorithm (dots), the stochastic algorithm using only 1\% of the contributions (error bars), and the Morse potential fitting the points obtained with the stochastic algorithm.} +\end{figure} + +The initial step involved the precise calculation of the CCSD(T) energy across various points along the potential curve. +We froze the six lowest molecular orbitals, specifically the $1s$ orbital of \ce{Cl} and the $1s$, $2s$, and $2p$ orbitals of \ce{Cu}, and correlated 34 electrons within 157 molecular orbitals. +The fitted Morse potential revealed a vibrational frequency of $\nu = \SI{414.7}{\per\centi\meter}$ and an equilibrium bond length of $r_e = \SI{3.92}{\bohr}$, aligning remarkably well with experimental values $\nu = \SI{414}{\per\centi\meter}$ and $r_e = \SI{3.88}{\bohr}$. + +Subsequently, we applied our semi-stochastic algorithm to estimate the perturbative triples correction, utilizing merely 1\% of the total contributions. +This approach yielded a hundredfold acceleration in computational efficiency, achieving statistical uncertainty within the range of \SI{1.3} to \SI{2.5}{\milli\hartree}. +The vibrational frequency and equilibrium distance estimated using this data, $\nu = \SI{415.0}{\per\centi\meter}$ and $r_e = \SI{3.92}{\bohr}$, demonstrated comparable precision to the full computational results. +Figure \ref{fig:cucl} illustrates the potential energy surface of \ce{CuCl}, displaying both the exact CCSD(T) energies and those estimated via the semi-stochastic method. + + +%%% -% CCSD(T) exact -%a = 2.65592 +/- 0.0403 (1.518%) -%De = 0.0718253 +/- 0.001879 (2.617%) -%re = 1.4105 +/- 0.00215 (0.1524%) -%E0 = -199.358 +/- 0.0003179 (0.0001595%) \section{Conclusion} \label{sec:conclusion} @@ -248,20 +255,20 @@ light in cm/s. %=================================================================% -%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%% \acknowledgements{ This work was supported by the European Centre of Excellence in Exascale Computing TREX --- Targeting Real Chemical Accuracy at the Exascale. This project has received funding from the European Union's Horizon 2020 — Research and Innovation program --- under grant agreement No.~952165. This work was performed using HPC resourced from CALMIP (Toulouse) under allocations p18005 and p22001.} -%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%% -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section*{Data availability statement} -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% The data that supports the findings of this study are available within the article and its supplementary material. -%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%% \bibliography{stochastic_triples} -%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%% \end{document} diff --git a/triples.org b/triples.org index 98f4dac..81a539e 100644 --- a/triples.org +++ b/triples.org @@ -2,7 +2,7 @@ * Stochastic formulation -The perturbative correction reads: + The perturbative correction reads:2te \[ E_{(T)} = \sum_{ijkabc} E_{ijk}^{abc} = \sum_{ijkabc} \frac{(4 W_{ijk}^{abc} + @@ -7697,7 +7697,7 @@ re = 3.9 De = 0.1 E0 = -2099.767 set xrange [2.7:5.2] -fit E(x) data using ($1*a0):2 via a, re, De, E0 +fit E(x) data using ($1*a0):2:3 via a, re, De, E0 plot E(x), data using ($1*a0):2:3 w err #+end_src @@ -7717,14 +7717,8 @@ E0 = -2099.77 +/- 0.0008014 (3.817e-05%) #+RESULTS: : 413.5302408975902 - #+CALL:freq(0.895573,0.0854261) - - #+RESULTS: - : 401.3588602143032 - ** CCSD(T) exact - #+name:cucl_ccsdt_ex | 1.50 | -2099.533616067071 | | 1.55 | -2099.590506349950 | | 1.60 | -2099.635662051331 | @@ -7750,6 +7744,20 @@ E0 = -2099.77 +/- 0.0008014 (3.817e-05%) | 2.60 | -2099.734033236772 | | 2.65 | -2099.729597826175 | + #+name:cucl_ccsdt_ex + | 1.55 | -2099.590506349950 | + | 1.65 | -2099.671184187604 | + | 1.75 | -2099.720045862965 | + | 1.85 | -2099.747811193906 | + | 1.95 | -2099.761752030920 | + | 2.05 | -2099.766727898670 | + | 2.15 | -2099.765956694308 | + | 2.25 | -2099.761562105614 | + | 2.35 | -2099.754944906474 | + | 2.45 | -2099.747028328725 | + | 2.55 | -2099.738443175793 | + | 2.65 | -2099.729597826175 | + #+begin_src gnuplot :var data=cucl_ccsdt_ex :results file :file cucl_ccsdt_ex.png reset a0 = 1.8897161646321 @@ -7767,16 +7775,16 @@ plot E(x), data using ($1*a0):2 [[file:cucl_ccsdt_ex.png]] #+begin_example -a = 0.853035 +/- 0.006204 (0.7273%) -re = 3.91994 +/- 0.002032 (0.05183%) -De = 0.100264 +/- 0.001906 (1.901%) -E0 = -2099.77 +/- 0.0001757 (8.367e-06%) +a = 0.840928 +/- 0.009924 (1.18%) +re = 3.92266 +/- 0.003278 (0.08358%) +De = 0.103446 +/- 0.002926 (2.829%) +E0 = -2099.77 +/- 0.0002473 (1.178e-05%) #+end_example - #+CALL:freq(0.853035,0.100264) + #+CALL:freq(0.840928,0.103446) #+RESULTS: - : 414.16742408686565 + : 414.7173810736408 #+begin_src gnuplot :var data=cucl_ccsdt :var data2=cucl_ccsdt_ex :results file :file cucl_ccsdt2.png