From 6765fd3dbfd18257d25c3f1c27df2958e5991c98 Mon Sep 17 00:00:00 2001 From: Emmanuel Giner Date: Wed, 25 Mar 2020 15:41:34 +0100 Subject: [PATCH] added some on top --- Revised_Manuscript/srDFT_SC.tex | 86 +++++ new/data_on_top/data | 54 +-- new/data_on_top/data_on_top_he.eps | 551 +++++++++++++++++++++++++++++ 3 files changed, 664 insertions(+), 27 deletions(-) create mode 100644 new/data_on_top/data_on_top_he.eps diff --git a/Revised_Manuscript/srDFT_SC.tex b/Revised_Manuscript/srDFT_SC.tex index 9c7ba4c..b0d81aa 100644 --- a/Revised_Manuscript/srDFT_SC.tex +++ b/Revised_Manuscript/srDFT_SC.tex @@ -143,6 +143,7 @@ % effective interaction \newcommand{\twodm}[4]{\mel{\Psi}{\psixc{#4}\psixc{#3} \psix{#2}\psix{#1}}{\Psi}} \newcommand{\murpsi}[0]{\mu_{\wf{}{\Bas}}({\bf r})} +\newcommand{\murcas}[0]{\mu_{\text{CASSCF}}({\bf r})} \newcommand{\murpsibas}[0]{\mu_{\wf{}{\Bas}}({\bf r})} \newcommand{\ntwo}[0]{n_{2}} \newcommand{\ntwohf}[0]{n_2^{\text{HF}}} @@ -267,6 +268,11 @@ \newcommand{\dbr}[1]{d\br{#1}} \newcommand{\PBEspin}{PBEspin} \newcommand{\PBEueg}{PBE-UEG-{$\tilde{\zeta}$}} +\newcommand{\ontopcas}{\langle n_2^{\text{CAS}}(\br{},\br{}) \rangle} +\newcommand{\ontopextrap}{\langle \mathring{n}_{2}^{\text{CAS}}(\br{},\br{}) \rangle} +\newcommand{\ontopcipsi}{\langle n_2^{\text{CIPSI}}(\br{},\br{}) \rangle} +\newcommand{\muaverage}{\langle \murcas \rangle} +\newcommand{\largemu}{E_{c,md}^{\mu \rightarrow \infty}} \newcommand{\LCT}{Laboratoire de Chimie Th\'eorique (UMR 7616), Sorbonne Universit\'e, CNRS, Paris, France} \newcommand{\ISCD}{Institut des Sciences du Calcul et des Donn\'ees, Sorbonne Universit\'e, Paris, France} @@ -681,6 +687,55 @@ The performance of each of these functionals is tested in the following. Note th \end{table*} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +\begin{table*} +\caption{Integral of the on-top pair density in real space at various levels of theory (see text for details) for N$_2$, N, O$_2$ an O in the aug-cc-pVXZ basis sets (X=D,T,Q).} +\begin{ruledtabular} +\begin{tabular}{lrcccccc} +%\begin{tabular}{lrcccccc} + + System & \tabc{Basis set} &\tabc{$\ontopcas$}& \tabc{$\ontopextrap$}& \tabc{$\ontopcipsi$}& \tabc{$\muaverage$} \\ +\hline +\ce{N2} & aug-cc-pVDZ & 0.587712 & 0.329839 & 0.513967 & 0.946 \\ + & aug-cc-pVTZ & 0.591622 & 0.385069 & 0.461386 & 1.328 \\ + & aug-cc-pVQZ & 0.592422 & 0.461393 & 0.419336 & 1.706 \\[0.1cm] + +\ce{N} & aug-cc-pVDZ & 0.172320 & 0.098115 & 0.127420 & 0.910 \\ + & aug-cc-pVTZ & 0.173022 & 0.113155 & 0.111726 & 1.263 \\ + & aug-cc-pVQZ & 0.173071 & 0.123332 & 0.106121 & 1.601 \\[0.1cm] +\hline +\ce{O2} & aug-cc-pVDZ & 1.164428 & 0.707757 & 0.971217 & 1.107 \\ + & aug-cc-pVTZ & 1.166676 & 0.797154 & 0.884406 & 1.545 \\ + & aug-cc-pVQZ & 1.167688 & 0.861134 & 0.841366 & 1.979 \\[0.1cm] + +\ce{O} & aug-cc-pVDZ & 0.513919 & 0.314972 & 0.416040 & 1.080 \\ + & aug-cc-pVTZ & 0.516070 & 0.369056 & 0.376896 & 1.499 \\ + & aug-cc-pVQZ & 0.516288 & 0.383353 & 0.362491 & 1.924 \\[0.1cm] +% & &\multicolumn{4}{c}{Estimated exact:\fnm[1] 665.4} \\[0.2cm] +\hline +\end{tabular} +\end{ruledtabular} +\label{tab:d1} +\end{table*} + +%& \tabc{$\largemu$} +% +%& -0.0520 +%& -0.0210 +%& -0.0104 +% +%& -0.0107 +%& -0.0064 +%& -0.0032 +% +%& -0.0726 +%& -0.0302 +%& -0.0152 +% +%& -0.0318 +%& -0.0138 +%& -0.0069 + \subsection{Computational details} We present potential energy curves of small molecules up to the dissociation limit @@ -775,6 +830,37 @@ We report in Figs.~\ref{fig:N2}, \ref{fig:O2}, and \ref{fig:F2} the potential en Just as in \ce{H10}, the accuracy of the atomization energies is globally improved by adding the basis-set correction and it is remarkable that $\pbeontXi$ and $\pbeontns$ provide again very similar results. The latter observation confirms that the dependence on the on-top pair density allows one to remove the dependence of any kind of spin polarization for a quite wide range of covalent bonds and also for an open-shell system like \ce{O2}. More quantitatively, an error below 1.0 mHa compared to the estimated exact valence-only atomization energy is found for \ce{N2}, \ce{O2}, and \ce{F2} with the aug-cc-pVTZ basis set using the $\pbeontns$ functional, whereas such a feat is far from being reached within the same basis set at the near-FCI level. In the case of \ce{F2} it is clear that the addition of diffuse functions in the double- and triple-$\zeta$ basis sets strongly improves the accuracy of the results, which could have be anticipated due to the strong breathing-orbital effect induced by the ionic valence-bond forms in this molecule. \cite{HibHumByrLen-JCP-94} It should be also noticed that when reaching the aug-cc-pVQZ basis set for \ce{N2}, the accuracy of the atomization energy slightly deteriorates for the $\pbeontXi$ and $\pbeontns$ functionals, but it remains nevertheless more accurate than the estimated FCI atomization energy and very close to chemical accuracy. +\manu{ +The overestimation of the atomization energy appearing for \ce{N2} in large basis sets reveals a kind of unbalanced treatment between the molecule and atoms in favour of the molecular system. +As the integral of the exact on-top pair density is proportional to the correlation energy in the large $\mu$ limit\cite{PazMorGorBac-PRB-06,FerGinTou-JCP-18} (see Eq. \eqref{eq:lim_mularge}), the accuracy of a given approximation to the exact on-top pair density will have a direct influence on the accuracy of the related correlation energy. +To quantify the quality of various flavour of on-top pair densities for a given system and a given basis set $\basis$, we define the following quantities +} +\begin{equation} + \ontopcas = \int \text{d}\br{}\, n_2^{\text{CASSCF}}(\br{},\br{}), +\end{equation} +\begin{equation} + \ontopextrap = \int \text{d}\br{}\, \ntwoextrap(n_2^{\text{CASSCF}}(\br{},\br{}),\murcas), +\end{equation} +\begin{equation} + \ontopcipsi = \int \text{d}\br{}\,n_2^{\text{CIPSI}}(\br{},\br{}), +\end{equation} +%\begin{equation} +% \largemu = \int \text{d}\br{}\, \frac{(-2+\sqrt{2})\sqrt{2\pi}}{3\left(\murcas\right)^3} \ntwoextrap(n_2^{\text{CASSCF}}(\br{},\br{}),\murcas), +%\end{equation} +\begin{equation} + \muaverage = \frac{1}{N_{e}}\int \text{d}\br{}\,n^{\text{CASSCF}}(\br{}) \,\, \murcas +\end{equation} +\manu{ +The quantity $n_2^{\text{CIPSI}}(\br{},\br{})$ is the on-top pair density of the largest CIPSI wave function obtained for each system in a given basis, which contains here at leas $10^7$ Slater determinants. All of these three quantities were computed excluding all contributions from the $1s$ orbitals. +We report in Table \ref{tab:d1} these quantities for N, \ce{N2}, O and \ce{O2} in different basis sets. +From this Table \ref{tab:d1} one can notice that integral of the on-top pair density at the CIPSI level is systematically lower than that at the CASSCF level, which is expected as the short-range correlation, digging the coulomb hone in a given basis set $\basis$ at near FCI level, is missing from the valence CASSCF wave function. +Also, the on-top pair density at the CIPSI level decreases roughly by $20\%$ between the aug-cc-pVDZ and aug-cc-pVQZ, whereas the on-top pair density at the CASSCF level is almost constant with respect to the basis set. +In order to estimate the integral of exact on-top pair density, we take as reference the value of $\ontopcipsi$ in the aug-cc-pVQZ basis set, although it is certainly an upper bound to the exact values. +Regarding $\ontopextrap$, such a quantity is directly linked to the basis set correction in the large $\mu$ limit. +The lowering of $\ontopextrap$ with respect to $\ontopcas$ is noticeable in a given basis set, but such $\ontopextrap$ globally increases when enlarging the basis set. This can be understood easily by remembering Eq. \eqref{eq:def_n2extrap} and realizing that the CASSCF on-top pair density is globally constant with the basis set whereas the value of $\murcas$ globally increases (as evidenced by $\muaverage$). +By comparing $\ontopextrap$ to the $\ontopcipsi$ in the aug-cc-pVQZ basis set, it is quite clear that the error +} + Regarding now the performance of the basis-set correction along the whole potential energy curve, it is interesting to notice that it fails to provide a noticeable improvement far from the equilibrium geometry. Acknowledging that the weak-correlation effects in these regions are dominated by dispersion interactions which are long-range effects, the failure of the present approximations for the complementary functional can be understood easily. Indeed, the whole scheme designed here is based on the physics of correlation near the electron-electron coalescence point: the local range-separation function $\mu(\br{})$ is based on the value of the effective electron-electron interaction at coalescence and the ECMD functionals are suited for short-range correlation effects. Therefore, the failure of the present basis-set correction to describe dispersion interactions is theoretically expected. We hope to report further on this in the near future. diff --git a/new/data_on_top/data b/new/data_on_top/data index f972ff4..b1ec59e 100644 --- a/new/data_on_top/data +++ b/new/data_on_top/data @@ -1,29 +1,29 @@ ========================= N2 - CAS CAS-extrap CIPSI -avdz -N2 0.587712 0.337529 0.513967 -N 0.172320 0.100855 0.127420 -------------------------------------- -avtz -N2 0.591622 0.388739 0.461386 -N 0.173022 0.114365 0.111726 -------------------------------------- -avqz -N2 0.592422 0.423673 0.419336 - 0.173071 0.124562 0.106121 -====================================== -====================================== -O2 - CAS CAS-extrap CIPSI -avdz -O2 1.164428 0.70775731 0.97121 -O 0.513919 0.31497205 0.41604 --------------------------------------- -avtz -O2 1.166676 0.79715452 0.88440 -O 0.516070 0.36905624 0.37689 --------------------------------------- -avqz -O2 1.167688 0.861134291 0.84136 -O 0.516288 0.398760190 0.36249 + CAS CAS-extrap CIPSI CIPSI-extrap +avdz +N2 0.58771 0.32983 0.51397 +N 0.17232 0.09811 0.12742 0.07343 +------------- ---------- ------------------------- +avtz +N2 0.59162 0.38506 0.46139 +N 0.17302 0.11315 0.11173 0.07414 +------------- ---------- ------------------------- +avqz +N2 0.59242 0.42006 0.41934 +N 0.17307 0.12333 0.10612 0.07582 +============= ==================================== +============= ==================================== +O2 + CAS CAS-extrap CIPSI +avdz +O2 1.16443 0.69100 0.97121 +O 0.51391 0.30623 0.41604 0.24934 +-------------------------------------------------- +avtz +O2 1.16668 0.78996 0.88440 +O 0.51607 0.35068 0.37689 0.25768 +-------------------------------------------------- +avqz +O2 1.16769 0.85386 0.84136 +O 0.51628 0.37922 0.36249 0.26283 diff --git a/new/data_on_top/data_on_top_he.eps b/new/data_on_top/data_on_top_he.eps new file mode 100644 index 0000000..5cb0939 --- /dev/null +++ b/new/data_on_top/data_on_top_he.eps @@ -0,0 +1,551 @@ +%!PS-Adobe-3.0 EPSF-3.0 +%%Creator: cairo 1.14.6 (http://cairographics.org) +%%CreationDate: Wed Mar 25 15:41:04 2020 +%%Pages: 1 +%%DocumentData: Clean7Bit +%%LanguageLevel: 2 +%%BoundingBox: 0 0 360 216 +%%EndComments +%%BeginProlog +save +50 dict begin +/q { gsave } bind def +/Q { grestore } bind def +/cm { 6 array astore concat } bind def +/w { setlinewidth } bind def +/J { setlinecap } bind def +/j { setlinejoin } bind def +/M { setmiterlimit } bind def +/d { setdash } bind def +/m { moveto } bind def +/l { lineto } bind def +/c { curveto } bind def +/h { closepath } bind def +/re { exch dup neg 3 1 roll 5 3 roll moveto 0 rlineto + 0 exch rlineto 0 rlineto closepath } bind def +/S { stroke } bind def +/f { fill } bind def +/f* { eofill } bind def +/n { newpath } bind def +/W { clip } bind def +/W* { eoclip } bind def +/BT { } bind def +/ET { } bind def +/pdfmark where { pop globaldict /?pdfmark /exec load put } + { globaldict begin /?pdfmark /pop load def /pdfmark + /cleartomark load def end } ifelse +/BDC { mark 3 1 roll /BDC pdfmark } bind def +/EMC { mark /EMC pdfmark } bind def +/cairo_store_point { /cairo_point_y exch def /cairo_point_x exch def } def +/Tj { show currentpoint cairo_store_point } bind def +/TJ { + { + dup + type /stringtype eq + { show } { -0.001 mul 0 cairo_font_matrix dtransform rmoveto } ifelse + } forall + currentpoint cairo_store_point +} bind def +/cairo_selectfont { cairo_font_matrix aload pop pop pop 0 0 6 array astore + cairo_font exch selectfont cairo_point_x cairo_point_y moveto } bind def +/Tf { pop /cairo_font exch def /cairo_font_matrix where + { pop cairo_selectfont } if } bind def +/Td { matrix translate cairo_font_matrix matrix concatmatrix dup + /cairo_font_matrix exch def dup 4 get exch 5 get cairo_store_point + /cairo_font where { pop cairo_selectfont } if } bind def +/Tm { 2 copy 8 2 roll 6 array astore /cairo_font_matrix exch def + cairo_store_point /cairo_font where { pop cairo_selectfont } if } bind def +/g { setgray } bind def +/rg { setrgbcolor } bind def +/d1 { setcachedevice } bind def +%%EndProlog +%%BeginSetup +%%BeginResource: font DejaVuSans +11 dict begin +/FontType 42 def +/FontName /DejaVuSans def +/PaintType 0 def +/FontMatrix [ 1 0 0 1 0 0 ] def +/FontBBox [ 0 0 0 0 ] def +/Encoding 256 array def +0 1 255 { Encoding exch /.notdef put } for +Encoding 32 /space put +Encoding 45 /hyphen put +Encoding 46 /period put +Encoding 47 /slash put +Encoding 48 /zero put +Encoding 49 /one put +Encoding 50 /two put +Encoding 51 /three put +Encoding 52 /four put +Encoding 53 /five put +Encoding 57 /nine put +Encoding 67 /C put +Encoding 70 /F put +Encoding 73 /I put +Encoding 90 /Z put +Encoding 97 /a put +Encoding 100 /d put +Encoding 101 /e put +Encoding 102 /f put +Encoding 103 /g put +Encoding 104 /h put +Encoding 105 /i put +Encoding 108 /l put +Encoding 110 /n put +Encoding 111 /o put +Encoding 112 /p put +Encoding 114 /r put +Encoding 116 /t put +Encoding 120 /x put +/CharStrings 30 dict dup begin +/.notdef 0 def +/space 1 def +/zero 2 def +/period 3 def +/nine 4 def +/one 5 def +/two 6 def +/three 7 def +/four 8 def +/five 9 def +/i 10 def +/n 11 def +/t 12 def +/e 13 def +/g 14 def +/r 15 def +/a 16 def +/l 17 def +/o 18 def +/f 19 def +/h 20 def +/hyphen 21 def +/p 22 def +/slash 23 def +/Z 24 def +/F 25 def +/C 26 def +/I 27 def +/x 28 def +/d 29 def +end readonly def +/sfnts [ +<0001000000090080000300106376742000691d39000011c8000001fe6670676d7134766a0000 +13c8000000ab676c7966ba1ad50d0000009c0000112c68656164086b92ef0000147400000036 +686865610d9f078b000014ac00000024686d74787fa40f6d000014d0000000786c6f63610000 +fcd4000015480000007c6d617870048b0671000015c400000020707265703b07f100000015e4 +0000056800020066fe96046605a400030007001a400c04fb0006fb0108057f0204002fc4d4ec +310010d4ecd4ec301311211125211121660400fc73031bfce5fe96070ef8f272062900020087 +ffe3048f05f0000b00170023401306a01200a00c91128c18091c0f1e031c151b1810fcecf4ec +310010e4f4ec10ee30012202111012333212111002273200111000232200111000028b9c9d9d +9c9d9d9d9dfb0109fef7fbfbfef701090550fecdfeccfecdfecd0133013301340133a0fe73fe +86fe87fe73018d0179017a018d00000100db000001ae00fe00030011b7008302011900180410 +fcec31002fec3037331523dbd3d3fefe00020081ffe3048705f00018002400584023071f1901 +860019a00aa504a00089161fa01091168c25071c1c21131e0022221c0d1b2510fcece4f4ecec +310010e4f4ec10e6fef5ee10ee111239304016c419c21ac01bc01cc01dc21ec41f07aa12bc12 +e912035d015d37351e01333212130e0123220035340033200011100021222601323635342623 +2206151416e14c9c4bc8d30f3ab26ce0fefb0110e201030111feb1fee54c9c013e889f9f8888 +9f9f1fb82426010d0112565c010febe60116fe73fe86fe9ffe5b1e0297baa2a1bbbba1a2ba00 +000100e10000045a05d5000a004040154203a00402a005810700a009081f061c03001f010b10 +d44bb00f5458b9000100403859ecc4fcec31002fec32f4ecd4ec304b5358592201b40f030f04 +025d3721110535253311211521fe014afe990165ca014afca4aa047348b848fad5aa00000001 +00960000044a05f0001c009e4027191a1b03181c11050400110505044210a111940da0149104 +00a00200100a02010a1c171003061d10fc4bb015544bb016545b4bb014545b58b90003ffc038 +59c4d4ecc0c011123931002fec32f4ecf4ec304b5358071005ed0705ed01b01c101117395922 +0140325504560556077a047a05761b87190704000419041a041b051c74007606751a731b741c +82008619821a821b821ca800a81b115d005d25211521353600373e0135342623220607353e01 +333204151406070600018902c1fc4c73018d33614da7865fd3787ad458e80114455b19fef4aa +aaaa7701913a6d974977964243cc3132e8c25ca5701dfeeb00000001009cffe3047305f00028 +0070402e0015130a86091f862013a0150da00993061ca020932391068c15a329161c13000314 +191c2620101c03141f09062910fc4bb016544bb014545b58b90009ffc03859c4c4d4ecf4ec11 +173939310010ece4f4e4ec10e6ee10ee10ee10ee11123930014009641e611f6120642104005d +011e0115140421222627351e013332363534262b013533323635342623220607353e01333204 +151406033f91a3fed0fee85ec76a54c86dbec7b9a5aeb6959ea39853be7273c959e6010c8e03 +251fc490ddf22525c33132968f8495a67770737b2426b42020d1b27cab0000020064000004a4 +05d50002000d0081401d010d030d0003030d4200030b07a00501038109010c0a001c0608040c +0e10dc4bb00b544bb00d545b58b9000cffc03859d43cc4ec32113931002fe4d43cec32123930 +4b5358071004c9071005c9592201402a0b002a0048005900690077008a000716012b0026012b +0336014e014f0c4f0d5601660175017a0385010d5d005d090121033311331523112311213503 +06fe0201fe35fed5d5c9fd5e0525fce303cdfc33a8fea00160c300000001009effe3046405d5 +001d005e4023041a071186101d1aa00714a010890d02a000810d8c07a41e171c010a031c000a +10061e10fc014bb016544bb014545b58b90010ffc038594bb00f5458b9001000403859c4d4ec +10c4ee310010e4e4f4ec10e6ee10fec410ee1112393013211521113e01333200151400212226 +27351e0133323635342623220607dd0319fda02c582cfa0124fed4feef5ec3685ac06badcaca +ad51a15405d5aafe920f0ffeeeeaf1fef52020cb3130b69c9cb624260000000200c100000179 +061400030007002b400e06be04b100bc020501080400460810fc3cec3231002fe4fcec30400b +1009400950096009700905015d1333112311331523c1b8b8b8b80460fba00614e900000100ba +00000464047b001300364019030900030e0106870e11b80cbc0a010208004e0d09080b461410 +fcec32f4ec31002f3ce4f4c4ec1112173930b46015cf1502015d011123113426232206151123 +1133153e013332160464b87c7c95acb9b942b375c1c602a4fd5c029e9f9ebea4fd870460ae65 +64ef00010037000002f2059e0013003840190e05080f03a9001101bc08870a0b080902040008 +10120e461410fc3cc4fc3cc432393931002fecf43cc4ec3211393930b2af1501015d01112115 +211114163b01152322263511233533110177017bfe854b73bdbdd5a28787059efec28ffda089 +4e9a9fd202608f013e00000000020071ffe3047f047b0014001b007040240015010986088805 +15a90105b90c01bb18b912b80c8c1c1b1502081508004b02120f451c10fcecf4ecc411123931 +0010e4f4ece410ee10ee10f4ee1112393040293f1d701da01dd01df01d053f003f013f023f15 +3f1b052c072f082f092c0a6f006f016f026f156f1b095d71015d0115211e0133323637150e01 +232000111000333200072e0123220607047ffcb20ccdb76ac76263d06bfef4fec70129fce201 +07b802a5889ab90e025e5abec73434ae2a2c0138010a01130143feddc497b4ae9e0000020071 +fe56045a047b000b0028004a4023190c1d0912861316b90f03b92623b827bc09b90fbd1a1d26 +1900080c4706121220452910fcc4ecf4ec323231002fc4e4ece4f4c4ec10fed5ee1112393930 +b6602a802aa02a03015d01342623220615141633323617100221222627351e013332363d010e +0123220211101233321617353303a2a59594a5a59495a5b8fefefa61ac51519e52b5b439b27c +cefcfcce7cb239b8023dc8dcdcc8c7dcdcebfee2fee91d1eb32c2abdbf5b6362013a01030104 +013a6263aa00000100ba0000034a047b001100304014060b0700110b03870eb809bc070a0608 +0008461210fcc4ec3231002fe4f4ecc4d4cc11123930b450139f1302015d012e012322061511 +231133153e0133321617034a1f492c9ca7b9b93aba85132e1c03b41211cbbefdb20460ae6663 +050500000002007bffe3042d047b000a002500bc4027191f0b17090e00a91706b90e1120861f +ba1cb923b8118c170c001703180d09080b1f030814452610fcecccd4ec323211393931002fc4 +e4f4fcf4ec10c6ee10ee11391139123930406e301d301e301f3020302130223f27401d401e40 +1f402040214022501d501e501f50205021502250277027851d871e871f8720872185229027a0 +27f0271e301e301f30203021401e401f40204021501e501f50205021601e601f60206021701e +701f70207021801e801f80208021185d015d0122061514163332363d01371123350e01232226 +353436332135342623220607353e0133321602bedfac816f99b9b8b83fbc88accbfdfb0102a7 +9760b65465be5af3f00233667b6273d9b4294cfd81aa6661c1a2bdc0127f8b2e2eaa2727fc00 +000100c100000179061400030022b7009702010800460410fcec31002fec30400d1005400550 +0560057005f00506015d13331123c1b8b80614f9ec0000020071ffe30475047b000b0017004a +401306b91200b90cb8128c1809120f51031215451810fcecf4ec310010e4f4ec10ee3040233f +197b007b067f077f087f097f0a7f0b7b0c7f0d7f0e7f0f7f107f117b12a019f01911015d0122 +06151416333236353426273200111000232200111000027394acab9593acac93f00112feeef0 +f1feef011103dfe7c9c9e7e8c8c7e99cfec8feecfeedfec7013901130114013800000001002f +000002f8061400130059401c0510010c08a906018700970e06bc0a02130700070905080d0f0b +4c1410fc4bb00a5458b9000b004038594bb00e5458b9000bffc038593cc4fc3cc4c412393931 +002fe432fcec10ee321239393001b640155015a015035d01152322061d012115211123112335 +333534363302f8b0634d012ffed1b9b0b0aebd0614995068638ffc2f03d18f4ebbab000100ba +000004640614001300344019030900030e0106870e11b80c970a010208004e0d09080b461410 +fcec32f4ec31002f3cecf4c4ec1112173930b2601501015d0111231134262322061511231133 +113e013332160464b87c7c95acb9b942b375c1c602a4fd5c029e9f9ebea4fd870614fd9e6564 +ef000001006401df027f028300030011b6009c020401000410dccc310010d4ec301321152164 +021bfde50283a400000200bafe5604a4047b0010001c003e401b1ab9000e14b90508b80e8c01 +bd03bc1d11120b471704000802461d10fcec3232f4ec310010e4e4e4f4c4ec10c4ee30400960 +1e801ea01ee01e04015d2511231133153e013332001110022322260134262322061514163332 +360173b9b93ab17bcc00ffffcc7bb10238a79292a7a79292a7a8fdae060aaa6461febcfef8fe +f8febc6101ebcbe7e7cbcbe7e700000000010000ff4202b205d50003002d4014001a01020102 +1a03000342029f008104020001032fc43939310010f4ec304b5358071005ed071005ed592201 +3301230208aafdf8aa05d5f96d0000000001005c0000051f05d500090090401b031107080708 +11020302420895008103950508030001420400060a10dc4bb009544bb00a545b58b90006ffc0 +3859c4d4e411393931002fecf4ec304b5358071005ed071005ed592201404005020a07180729 +02260738074802470748080905030b08000b16031a08100b2f0b350339083f0b47034a084f0b +55035908660369086f0b770378087f0b9f0b165d005d13211501211521350121730495fc5003 +c7fb3d03b0fc6705d59afb6faa9a04910000000100c90000042305d500090029401206950402 +95008104ad08050107031c00040a10fcec32d4c431002fecf4ec10ee30b20f0b01015d132115 +21112115211123c9035afd700250fdb0ca05d5aafe48aafd370000010073ffe3052705f00019 +0036401a0da10eae0a951101a100ae04951791118c1a07190d003014101a10fcec32ec310010 +e4f4ecf4ec10eef6ee30b40f1b1f1b02015d01152e0123200011100021323637150e01232000 +111000213216052766e782ff00fef00110010082e7666aed84feadfe7a0186015386ed0562d5 +5f5efec7fed8fed9fec75e5fd34848019f01670168019f470000000100c90000019305d50003 +002eb700af02011c00040410fc4bb0105458b9000000403859ec31002fec3001400d30054005 +500560058f059f05065d13331123c9caca05d5fa2b000001003b000004790460000b01434046 +0511060706041103040707060411050401020103110202010b110001000a11090a0101000a11 +0b0a0708070911080807420a070401040800bf05020a0704010408000208060c10d44bb00a54 +4bb00f545b4bb010545b4bb011545b58b90006004038594bb0145458b90006ffc03859c4d4c4 +11173931002f3cec321739304b5358071005ed071008ed071008ed071005ed071005ed071008 +ed071008ed071005ed59220140980a04040a1a04150a260a3d04310a55045707580a660a7601 +7a047607740a8d04820a99049f049707920a900aa601a904af04a507a30aa00a1c0a03040505 +090a0b1a03150515091a0b2903260525092a0b200d3a013903370534073609390b300d490346 +0545094a0b400d590056015902590357055606590756085609590b500d6f0d78017f0d9b0194 +07ab01a407b00dcf0ddf0dff0d2f5d005d09022309012309013309010464fe6b01aad9febafe +bad901b3fe72d9012901290460fddffdc101b8fe48024a0216fe71018f0000020071ffe3045a +06140010001c003840191ab9000e14b905088c0eb801970317040008024711120b451d10fcec +f4ec323231002fece4f4c4ec10c4ee30b6601e801ea01e03015d0111331123350e0123220211 +100033321601141633323635342623220603a2b8b83ab17ccbff00ffcb7cb1fdc7a79292a8a8 +9292a703b6025ef9eca86461014401080108014461fe15cbe7e7cbcbe7e7013500b800cb00cb +00c100aa009c01a600b800660000007100cb00a002b20085007500b800c301cb0189022d00cb +00a600f000d300aa008700cb03aa0400014a003300cb000000d9050200f4015400b4009c0139 +0114013907060400044e04b4045204b804e704cd0037047304cd04600473013303a2055605a6 +0556053903c5021200c9001f00b801df007300ba03e9033303bc0444040e00df03cd03aa00e5 +03aa0404000000cb008f00a4007b00b80014016f007f027b0252008f00c705cd009a009a006f +00cb00cd019e01d300f000ba018300d5009803040248009e01d500c100cb00f600830354027f +00000333026600d300c700a400cd008f009a0073040005d5010a00fe022b00a400b4009c0000 +0062009c0000001d032d05d505d505d505f0007f007b005400a406b80614072301d300b800cb +00a601c301ec069300a000d3035c037103db0185042304a80448008f0139011401390360008f +05d5019a0614072306660179046004600460047b009c00000277046001aa00e904600762007b +00c5007f027b000000b4025205cd006600bc00660077061000cd013b01850389008f007b0000 +001d00cd074a042f009c009c0000077d006f0000006f0335006a006f007b00ae00b2002d0396 +008f027b00f600830354063705f6008f009c04e10266008f018d02f600cd03440029006604ee +00730000140000960000b707060504030201002c2010b002254964b040515820c859212d2cb0 +02254964b040515820c859212d2c20100720b00050b00d7920b8ffff5058041b0559b0051cb0 +032508b0042523e120b00050b00d7920b8ffff5058041b0559b0051cb0032508e12d2c4b5058 +20b0fd454459212d2cb002254560442d2c4b5358b00225b0022545445921212d2c45442d2cb0 +0225b0022549b00525b005254960b0206368208a108a233a8a10653a2d000001000000025999 +fadb104a5f0f3cf5001f080000000000d184f71800000000d184f718f7d6fc4c0e5909dc0000 +0008000000010000000000010000076dfe1d00000efef7d6fa510e5900010000000000000000 +000000000000001e04cd0066028b000005170087028b00db05170081051700e1051700960517 +009c051700640517009e023900c1051200ba0323003704ec007105140071034a00ba04e7007b +023900c104e5007102d1002f051200ba02e30064051400ba02b20000057b005c049a00c90596 +0073025c00c904bc003b05140071000000000000004400000044000000c8000000f0000001c4 +00000234000003340000041c000004d800000598000005e800000660000006dc000007b00000 +0878000008e800000a1400000a5000000af400000b8c00000c0400000c3000000cd000000d1c +00000ddc00000e3000000ec800000f10000010940000112c00010000001e0354002b0068000c +000200100099000800000415021600080004b8028040fffbfe03fa1403f92503f83203f79603 +f60e03f5fe03f4fe03f32503f20e03f19603f02503ef8a4105effe03ee9603ed9603ecfa03eb +fa03eafe03e93a03e84203e7fe03e63203e5e45305e59603e48a4105e45303e3e22f05e3fa03 +e22f03e1fe03e0fe03df3203de1403dd9603dcfe03db1203da7d03d9bb03d8fe03d68a4105d6 +7d03d5d44705d57d03d44703d3d21b05d3fe03d21b03d1fe03d0fe03cffe03cefe03cd9603cc +cb1e05ccfe03cb1e03ca3203c9fe03c6851105c61c03c51603c4fe03c3fe03c2fe03c1fe03c0 +fe03bffe03befe03bdfe03bcfe03bbfe03ba1103b9862505b9fe03b8b7bb05b8fe03b7b65d05 +b7bb03b78004b6b52505b65d40ff03b64004b52503b4fe03b39603b2fe03b1fe03b0fe03affe +03ae6403ad0e03acab2505ac6403abaa1205ab2503aa1203a98a4105a9fa03a8fe03a7fe03a6 +fe03a51203a4fe03a3a20e05a33203a20e03a16403a08a4105a096039ffe039e9d0c059efe03 +9d0c039c9b19059c64039b9a10059b19039a1003990a0398fe0397960d0597fe03960d03958a +410595960394930e05942803930e0392fa039190bb0591fe03908f5d0590bb039080048f8e25 +058f5d038f40048e25038dfe038c8b2e058cfe038b2e038a8625058a410389880b0589140388 +0b03878625058764038685110586250385110384fe038382110583fe0382110381fe0380fe03 +7ffe0340ff7e7d7d057efe037d7d037c64037b5415057b25037afe0379fe03780e03770c0376 +0a0375fe0374fa0373fa0372fa0371fa0370fe036ffe036efe036c21036bfe036a1142056a53 +0369fe03687d036711420566fe0365fe0364fe0363fe0362fe03613a0360fa035e0c035dfe03 +5bfe035afe0359580a0559fa03580a035716190557320356fe03555415055542035415035301 +1005531803521403514a130551fe03500b034ffe034e4d10054efe034d10034cfe034b4a1305 +4bfe034a4910054a1303491d0d05491003480d0347fe0346960345960344fe0343022d0543fa +0342bb03414b0340fe033ffe033e3d12053e14033d3c0f053d12033c3b0d053c40ff0f033b0d +033afe0339fe033837140538fa033736100537140336350b05361003350b03341e03330d0332 +310b0532fe03310b03302f0b05300d032f0b032e2d09052e10032d09032c32032b2a25052b64 +032a2912052a25032912032827250528410327250326250b05260f03250b0324fe0323fe0322 +0f03210110052112032064031ffa031e1d0d051e64031d0d031c1142051cfe031bfa031a4203 +1911420519fe031864031716190517fe031601100516190315fe0314fe0313fe031211420512 +fe0311022d05114203107d030f64030efe030d0c16050dfe030c0110050c16030bfe030a1003 +09fe0308022d0508fe030714030664030401100504fe03401503022d0503fe0302011005022d +0301100300fe0301b80164858d012b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b002b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b +2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b2b1d00> +] def +/f-0-0 currentdict end definefont pop +%%EndResource +%%EndSetup +%%Page: 1 1 +%%BeginPageSetup +%%PageBoundingBox: 0 0 360 216 +%%EndPageSetup +q 4 2 349 207 rectclip q +0 g +0.25 w +0 J +0 j +[] 0.0 d +3.8 M q 1 0 0 -1 0 216 cm +48.051 186.75 m 51.75 186.75 l S Q +q 1 0 0 -1 0 216 cm +345.148 186.75 m 341.449 186.75 l S Q +BT +8 0 0 8 22.570312 26.25 Tm +/f-0-0 1 Tf +( 0.09)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 157.398 m 51.75 157.398 l S Q +q 1 0 0 -1 0 216 cm +345.148 157.398 m 341.449 157.398 l S Q +BT +8 0 0 8 27.660156 55.597656 Tm +/f-0-0 1 Tf +( 0.1)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 128.102 m 51.75 128.102 l S Q +q 1 0 0 -1 0 216 cm +345.148 128.102 m 341.449 128.102 l S Q +BT +8 0 0 8 22.570312 84.898438 Tm +/f-0-0 1 Tf +( 0.11)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 98.75 m 51.75 98.75 l S Q +q 1 0 0 -1 0 216 cm +345.148 98.75 m 341.449 98.75 l S Q +BT +8 0 0 8 22.570312 114.25 Tm +/f-0-0 1 Tf +( 0.12)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 69.398 m 51.75 69.398 l S Q +q 1 0 0 -1 0 216 cm +345.148 69.398 m 341.449 69.398 l S Q +BT +8 0 0 8 22.570312 143.597656 Tm +/f-0-0 1 Tf +( 0.13)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 40.102 m 51.75 40.102 l S Q +q 1 0 0 -1 0 216 cm +345.148 40.102 m 341.449 40.102 l S Q +BT +8 0 0 8 22.570312 172.898438 Tm +/f-0-0 1 Tf +( 0.14)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 10.75 m 51.75 10.75 l S Q +q 1 0 0 -1 0 216 cm +345.148 10.75 m 341.449 10.75 l S Q +BT +8 0 0 8 22.570312 202.25 Tm +/f-0-0 1 Tf +( 0.15)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 186.75 m 48.051 183.051 l S Q +q 1 0 0 -1 0 216 cm +48.051 10.75 m 48.051 14.449 l S Q +BT +8 0 0 8 37.847656 16.949219 Tm +/f-0-0 1 Tf +( 0.15)Tj +ET +q 1 0 0 -1 0 216 cm +90.5 186.75 m 90.5 183.051 l S Q +q 1 0 0 -1 0 216 cm +90.5 10.75 m 90.5 14.449 l S Q +BT +8 0 0 8 82.84375 16.949219 Tm +/f-0-0 1 Tf +( 0.2)Tj +ET +q 1 0 0 -1 0 216 cm +132.949 186.75 m 132.949 183.051 l S Q +q 1 0 0 -1 0 216 cm +132.949 10.75 m 132.949 14.449 l S Q +BT +8 0 0 8 122.746094 16.949219 Tm +/f-0-0 1 Tf +( 0.25)Tj +ET +q 1 0 0 -1 0 216 cm +175.398 186.75 m 175.398 183.051 l S Q +q 1 0 0 -1 0 216 cm +175.398 10.75 m 175.398 14.449 l S Q +BT +8 0 0 8 167.742188 16.949219 Tm +/f-0-0 1 Tf +( 0.3)Tj +ET +q 1 0 0 -1 0 216 cm +217.801 186.75 m 217.801 183.051 l S Q +q 1 0 0 -1 0 216 cm +217.801 10.75 m 217.801 14.449 l S Q +BT +8 0 0 8 207.597656 16.949219 Tm +/f-0-0 1 Tf +( 0.35)Tj +ET +q 1 0 0 -1 0 216 cm +260.25 186.75 m 260.25 183.051 l S Q +q 1 0 0 -1 0 216 cm +260.25 10.75 m 260.25 14.449 l S Q +BT +8 0 0 8 252.59375 16.949219 Tm +/f-0-0 1 Tf +( 0.4)Tj +ET +q 1 0 0 -1 0 216 cm +302.699 186.75 m 302.699 183.051 l S Q +q 1 0 0 -1 0 216 cm +302.699 10.75 m 302.699 14.449 l S Q +BT +8 0 0 8 292.496094 16.949219 Tm +/f-0-0 1 Tf +( 0.45)Tj +ET +q 1 0 0 -1 0 216 cm +345.148 186.75 m 345.148 183.051 l S Q +q 1 0 0 -1 0 216 cm +345.148 10.75 m 345.148 14.449 l S Q +BT +8 0 0 8 337.492188 16.949219 Tm +/f-0-0 1 Tf +( 0.5)Tj +ET +q 1 0 0 -1 0 216 cm +48.051 10.75 297.098 176 re S Q +BT +0 8 -8 0 10.699219 74.636719 Tm +/f-0-0 1 Tf +(integral of the on-top)Tj +8 0 0 8 189.941406 3 Tm +(1/Z)Tj +2.63916 23.862305 Td +[(FCI e)18(xtrapolated on-top)]TJ +ET +0.580392 0 0.827451 rg +q 1 0 0 -1 0 216 cm +310.852 19.102 m 334.949 19.102 l S Q +q 1 0 0 -1 0 216 cm +345.148 179.699 m 203.648 169.301 l 132.949 161.75 l 90.5 157.449 l 62.199 + 154.648 l S Q +q 1 0 0 -1 0 216 cm +342.148 179.699 m 348.148 179.699 l S Q +q 1 0 0 -1 0 216 cm +345.148 176.699 m 345.148 182.699 l S Q +q 1 0 0 -1 0 216 cm +200.648 169.301 m 206.648 169.301 l S Q +q 1 0 0 -1 0 216 cm +203.648 166.301 m 203.648 172.301 l S Q +q 1 0 0 -1 0 216 cm +129.949 161.75 m 135.949 161.75 l S Q +q 1 0 0 -1 0 216 cm +132.949 158.75 m 132.949 164.75 l S Q +q 1 0 0 -1 0 216 cm +87.5 157.449 m 93.5 157.449 l S Q +q 1 0 0 -1 0 216 cm +90.5 154.449 m 90.5 160.449 l S Q +q 1 0 0 -1 0 216 cm +59.199 154.648 m 65.199 154.648 l S Q +q 1 0 0 -1 0 216 cm +62.199 151.648 m 62.199 157.648 l S Q +q 1 0 0 -1 0 216 cm +319.898 19.102 m 325.898 19.102 l S Q +q 1 0 0 -1 0 216 cm +322.898 16.102 m 322.898 22.102 l S Q +0 g +BT +8 0 0 8 264.675781 184.597656 Tm +/f-0-0 1 Tf +(FCI on-top)Tj +ET +0 0.619608 0.45098 rg +q 1 0 0 -1 0 216 cm +310.852 28.398 m 334.949 28.398 l S Q +q 1 0 0 -1 0 216 cm +345.148 22.75 m 203.648 64.25 l 132.949 85.398 l 90.5 97.25 l 62.199 104.699 + l S Q +q 1 0 0 -1 0 216 cm +342.148 19.75 m 348.148 25.75 l S Q +q 1 0 0 -1 0 216 cm +342.148 25.75 m 348.148 19.75 l S Q +q 1 0 0 -1 0 216 cm +200.648 61.25 m 206.648 67.25 l S Q +q 1 0 0 -1 0 216 cm +200.648 67.25 m 206.648 61.25 l S Q +q 1 0 0 -1 0 216 cm +129.949 82.398 m 135.949 88.398 l S Q +q 1 0 0 -1 0 216 cm +129.949 88.398 m 135.949 82.398 l S Q +q 1 0 0 -1 0 216 cm +87.5 94.25 m 93.5 100.25 l S Q +q 1 0 0 -1 0 216 cm +87.5 100.25 m 93.5 94.25 l S Q +q 1 0 0 -1 0 216 cm +59.199 101.699 m 65.199 107.699 l S Q +q 1 0 0 -1 0 216 cm +59.199 107.699 m 65.199 101.699 l S Q +q 1 0 0 -1 0 216 cm +319.898 25.398 m 325.898 31.398 l S Q +q 1 0 0 -1 0 216 cm +319.898 31.398 m 325.898 25.398 l S Q +0 g +q 1 0 0 -1 0 216 cm +48.051 10.75 297.098 176 re S Q +Q Q +showpage +%%Trailer +end restore +%%EOF