From 6afeb95e53eb87f5c9ae8061c07bce3b13a4ecbe Mon Sep 17 00:00:00 2001 From: sjplimp Date: Mon, 4 Nov 2013 17:25:25 +0000 Subject: [PATCH] git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@10944 f3b2605a-c512-4ea7-a41b-209d697bcdaa --- doc/Eqs/pair_comb1.jpg | Bin 15231 -> 19008 bytes doc/Eqs/pair_comb1.tex | 11 +-- doc/Eqs/pair_comb2.jpg | Bin 17322 -> 40980 bytes doc/Eqs/pair_comb2.tex | 22 +++-- doc/Section_commands.html | 18 ++-- doc/Section_commands.txt | 2 +- doc/pair_comb.html | 198 +++++++++++++------------------------- doc/pair_comb.txt | 191 ++++++++++++------------------------ 8 files changed, 159 insertions(+), 283 deletions(-) diff --git a/doc/Eqs/pair_comb1.jpg b/doc/Eqs/pair_comb1.jpg index 0e7fc4e35275a5f15cb2e4631d7625af6a2bc4e3..d92f37be1c2014f37f84fcdf972b7a411aad6911 100644 GIT binary patch literal 19008 zcmb@tbySp5)HgZ^NJ&a49fL67poD;QD9j8!lyoaKAl)Dc5|YCVA<`W~cdB%E4kg`4 zr+9ha`+fI&-?i?#-@R+y=UMCg^Q`CWbN1P1=WqX;`L_(9fGI!~05~`}0Lc9d_=g3^ z1MqS2@bGZ)?+1K*{D%aO2np^VauO1vN0j7LRFveDl#gi{=^j(l(@;`AVS7T)#KgkF z@|cbt$j%I8WM*Oh_X9ve2uQ%?!@;2h{ObTP0&oDh_ayF#{LklpT0GnbIE43K)hPft zxOh0Y4{-5u@rlUsZ~-_EaPa{64=Je#9<%YwQnPCi(g>IZ(sDQn%E5!;=tLB>-9N{p ztALzB@`{=|y5`PbvVL`pKKVH|j!ek>HwPg64?@y_yZ&4d437H$ZO~;Z-M4GZUDrQQ z-+aipFJ>W?h-)`>SsX4V#__WZiy6Y&L;@EaC|O@>|9H*0R*l3w#uNm$JaD_B-^HUW zcPIV-xEkH>;g1v?8`w5x70!A)@uL0#9vd$C=(lx!Y#Y!22l#d#CCrFMDI_wj<5V(y z_HI>u9?2*K|XgzxqD>I;Jzu;|7sA)1?lAFz0XKNKLACdH07~PN9q{ zwq8q4;YPqU*}-2l!+4;h41{M7bN&Gu!A*D3?*sk;3=;nVerxw%m|XkBNR7Xy+h>Yx zE)V8eHBDp!{sYj|K5P>kExSlTMD((sxWwjVbVD?w*>W3QA-GT(9-UrhfK<({O+CL$ zP6M*5p%T6A8ogRZ(b(SlQdI-OWQg-OC}mYuA5( z6GJZtMTf3E6Q4o&cN<4RCX}#dU8+|{?U*c(oJU`w7eEf&(vrBnTX-Mp<-;C%ia)3I zi;0xYA@hRe4w|Q>p0*o|3PQ!;%Q#NMb(Dh)o0{C4&}giEVrFtwMr$aHf)<5pYX3!7 zLkiGBKB&&_*r_=8ec;aK+z+{`^qKE3B5yNon2Xt5x&R9W=$E|*wzRj$@9Izf#h@R88%E1#3j>5|!1rUr7|EZhL zm}f%GgNS@PJj;({B;d;yO`(BXq&eRH%5|#J=96a@H;AbCjt?b9VaaudU^@Yla0qvu%xBilPRZSEt z?PYG*KmCjmjm-f+*{ehYqk;$(<4j7=p3RKgI=x%8ms99b735dE!S#8?5b0Kh(<5ko zX{hHRnz&7*|1I=mtu&483NlODGj0>gWr8*qVCZc9`99vWu{qC70is!~r1bIAe;zTH zB`H9~Ti@Z4;c}Mc`}NwnUK{~)ai(DcWAz<)=5>sa`zFGsQ+WLsJ8gU}ry7tVfWPtT zW^VFMFyt;(`gXcxb|h2A`L}=p<%)2$onL2gCg>XJYn$buM~YhYu=`qRtZYrm7zdVa zjy^%;66=DQc_D0Eu%i~wS6W-ba_uc*otV1Fz}}0iSSe5^%+1cASoe?-!N=PkQffv_ z1k8_)uF*B1t1JftFC*pJaArXmg^Xbc;>DHChx)t_vMmp{XvB7<-+7LJFg| za<&40%Un|$`MqqDMw;r_Rd23{J9#s}(%5ry^VVp`@+o4WVtBh#d3lBz?a5Ts;Y7;G z5pauY9QbqxLQ#I0;#Ca`uEmF*(#HQ&U!gCvsF2w$IBmFY!}h6E>-dG`#!3TXefo-j+P}_tly)u4Y@|~LYSh6FJfZ#Oj;W( z-mm;RP{~w%vV7YK0ugh}9PuiG9btHrt>^9 zY0X)dft%AZ=kkiIJ7RQB5mB9uM0eoSFvZ`hO*JGn$og+So9hT};NgMikHo6kC89!D zWU?fV7EfiT!d^}P53A6e*K5nDXeB*d-A6qAvj(h>`vD5MS8D*?bHP1W3O&wl_M&4fn}`z?0%Zi@!ZL+sxa3$Qn@>A=T~KQ`x5$+?}_ zfM%!DnSHg7|I^5|GTY?O8oUq=)pz3Ll^Y}0LJG+eOz0j*x=Ej#?yCK2%qLe8xMDz9 zB2%>D_xuB;=L|-V4PN~(xt>cz;@jXqG;llQ9(;>|kSk?#wgG z5#bpl5h8ku)Ur3HhQtEJZR*yK4r423HF4LI!t81be%ootN=(!R5%~Q9O*Ck2968CJHt-x$(jCQR?SRSTpSqPF5JBD)khvRN+|Cs%Y{*m>jnKMPQb^2-2 zz%!-6{PKQ;@X{eIn%W-2BTg(uo$eP%I{s`veIhM;En?QHKhjRwT(O8<5o+AVcXCL+ zLJD(>XHAV>Mlw}%KRf&w+gUmXS$Zm9q{MC9a(S{@=Ksp9tK`qo4Ib}vKm|Pm$Ohkb zt!8%;CRGHF?&azQ#_t-#j}4ufwu^yHe~TnX4RX#g?Ky@PvkzU5(H{w)(E8OGjlSX-}KEYYOWS-Pi(0Of18Fj4N1S1@P9uzq^0 zb1mtg>f8QR*F%2#tHdnr*RP-K8c#h){1$vXs31(NC(K8nAH(G=rQS0D2hmw`gI_Rxk_7mjiEjP zxL8^A8OlmUKu1|-+=*yH{gMTFr)U!Tx5(8O{{ejWKNbCaZ5gccb%z4!=bBn59Sc{D z5DLiIq2uiWzh=;&krI4wNHEHSS01(U1%_*Ve#|0<?n#0`#K6X5Khdo@2)ymUEznmlGft588dbv{l zeBH}QD!EYKa@uGf07d4!EFNW+t%!8pj}OG1*FFM9zXq6W(mFQ;lI4XrO~0htsM3)? zuE_*&z}7cM9%L#ucDMC!4=H|VlqNp<*munB{?Pe}`W@qf83v5sT%8k_hZQiG^Vj?l zw_AniV-EkjSiuWte{#g1qKP!^wqbk|y3P_J$Dr0oZ}S55<^dHo(Nu8?PRpyPw?mF1 zD^~9X9jipKni~fDLRn$=mRBID?NMgJB3-B@@3V%r*VuKXH%*_X z`rmP_J5cLVPmQ1QNhIp1)i0C(rg0Lvp=kB%aMT%f|2@3y4`QM0Hdx0lfvrEccb@0A>E)kQlaoIP@OR;Sj^ z;*P#wLete?$ggbSPWot?jiu;g+$gcRpxElT@|?KkAc-K}F*K8Z+f2zL;lJl+h#H}t zdf~NMh!w9+U51XrP)l(Z(=WVElzb|^+-HdAmsD5}!e_aS;7@V*=&6|4Y3a9&ERC+l zp~+k^mAdD_r+IOE^Oy6u!Q>9BN>J$(yz7E*ZSaie9u0z3k0J(km2VMWE|l3=lJ`Ui zoj``RtRs12RTs;%0K!^Y0_o^=5k@&%_H+^t6ytif>~ZuK&X*jUTwNHq-HE^rl~$z> ztf0{uRYZ9Uh>LZJfN#IKb%Mg+32n@>91Xp(B}|t62bHT2=PNBrYRp2i48Iv`Rsu0O zcU-jj@Huj()cj$Lwp0FGc#Lbp*a{5~O+BOu>&-c_9;UCS3>?q!6!+rv1Y`x zI;kwn$lsToC)|FmJK&E-NQT9x2pAM5U4kK>YHqZ%j5BM`=&TpZKPbyzDRnG9k^~Q@ zPO|jJv65Tnf|voN8^?xDvo4u+vHH5t7K0i}@+_>JkYo|xRnI>OxvLg0_Ip(=u@`A4 z>^|IN;^o$?0J5#&TQdW5;#w$IPfB;-Xd*hkj`i34^$n77XFuf<`N6|B3Sw~AtYzXW zfE1C~q{9Hi(JP)~sngSNa4YQImNdItuc_H}xU-alePQ ztB?Nx?It(h8pbkS_Rps*Etyb_gKzO?FTv+O8Dp229#{Ha6<)*9U;bh|m$H{zgl9Q; zdY`!C;a(Zdv0^w$8J{jg;xyUpUwG`dDN=Q3;0@XxkclgYp5aI3m;R1 zo1yr4b0{)(<;pqLu8BWNIC#j`7BOnsIO|rr{d?XcY-((pbI%^3QYikg(4#5Ak*?UY zm4tv_D^P5S@3)(q&eVe>t5QPy44ll~$=#t_;^r07HG3mR4M$za=9TJIxyMCe3ADmd z7_OL*D~=u(jeQ-Mwx(!Yj3AGIqaJWqDFm-J?8MOYxnnPzRE^Q){Nuaw7bd%jJMv8E zu-D!ZU;VAzQn|}6oaP%PR&-o_snhX%uCBGuQAKUqAF*UUL54>qQ6S}~q%!JhY`1zk=--8$9Nln0FYPs)d&x*&D9Yp76W`RwSjYPcc`Js|a4 zVd{Ug5`HQ6>Vn}%)IJ)VZw#3|mY|QzRal$iRZy{vvuu;G^=njQ-X%MXs4oC>b*B@|Rxp zrY(liZs1amMB#g^p34AKdmFO!pJ z+MfrOxG)%)FR%om{k?$EbFPp}`RptE$+^YCufb!RFFQMC4~*y8tK8pHVGFO?^{TU5z}p^)Hz30hv6C+4+(O2L zmO?91ZW`{r?HU%b1%1c5eyi9YbX-%T^(qdJ9G+u_9AEgtJ#jmH^EY&muht?g>zdw+ zm?!!+_h@-iQVw2Fjm?7?aDS;2qt{{pkKv>R;A)KcN(SM~^gqD5{}`%dc(11IF8%01 zWi+<6@%&mBvNx8D5;aUy8h;rI)TVf-Vu1!us3@2lPM#m*L z1UiE;UJ75W{a$;N)Ez$CWTa1sV!D^RoyPdxxQLFQYtx=%@n~vMh*Qz!o{SsBCI1_7IcmAataUb@`B3i_{ zSM$UZi%S@u)5NTcmo}0_7aHf<${F!!oULK@9_tVOPzicI<%_tHa304HtPZT{{^4fS zorOUpx_fZMAt#Ngx#n%}!&Ns`7Oz{p)ZA4TeU9vNf3^kHR+bVW7;ZV0NO(C6!x&C! z_RrNK2b$G52A^89m1MK02(#V=zY*25!rh8x&R2};#fRQ>KpEdx$4Hi`o4jb zZ#&>AhDTu;cl+Iks+B8WT1MDQ$o~Q2iWI!qq>L6GfMQU;)l<(4_j%M>D}8u^CTz=FSonxtF)KCgJ+sTp$AB>l~3wD{Cm*Gfg# zu_oC50$h_LgiP9j#mXUl0_9DEM3yJ-=+rWP?5H&HMm^g&HaFmW!=36D|3b;B3#Rkw zc(Nln^9N+=irmOEH3TW^^8}5Q?)_@%#I|@ywJA)`!!!{b1k<9>O0|zwXt8UZW+&}q zwCezd)bWWi*-Q$l*9d1&r-lGBmp^CWFt0T1CGHw+-^A0-l6blI+J?Tc$Ka52f;1pJ z0A5-7K6lBwa$D7V$=#Q(8Z{_#tLCw4;PX7f3a*A69uA8A zw0ocb4-rWm0e)xpTRoHLwOPJIxhaAWDwNAfTRQ9vFqLoGRq_3m2R7$n9AIIibk!yH z7R?P__bxye1m|x({B%yU6vcD=-CsOeaH?iLNUM;{V63`pnjNk=fVPWaLi%~^8+-O& zc**&uL>5aVIrPCXc6BG>Q5{YP+(_zH!&=RpU6WL$De3d5X2Z3Gmubm!mXb7H=t28o zKQwiTj2Xxf#VOyXl~X_(<6xryd$%)dTHtCd?YqIf{_IFOwmg3-BS-o0vEQ`*$f>Mk zW|>Oh!NVFaq(ijj2s_ql?h!3~rL;u67yb$Kx7RUzOpz6X#x0l7YJrO#C5w;W#aVi} zT5>vPbUt_qcaJ^!s-%LH{GmSvxhBqB*Ew%QoU69>*0c7s_Ejf32+VV;K^uCA)SO_q zZ)GR*V;be)536~VwmTDqDaT1jqm?jfAzF4aN|YX%;lAx4)YVOk;HSBU=1K18Je**nbt7ApJNA(<|_5^TJ#ip8A*IMZj2| z^m_Y!&MDx)c2m21;{(+1N7w6Ua7Gpb{BW2Prg}?-C!bx0NMEr{8GoJJ7+`b9!rX1R zDF3`+=<-;!Y;A@<%w1ZV;hK2#xufx-O}kH>L8d6RxJ~G5XSvV`w9sPAdNy1r*Q##l zo8^Un{>;ghUR23oMfnW{yI)ro6V4h6W9h5LUKx)$#{F@mS%+oODAk!Z>l%K3p2!G= z#FH~?xG`L)Im5==b|%iBn#Dw#PwhVXkQN(e9}K!c*Ri-vI+WaCY7k_?Zdc5fZzM zsKz;f9%(sbEnwSD<+QgG7~e~5B>9UO(rehLWz=SoqtF}?_30emQg8KL#1Yqr*6vFd zEuB?G@~T3+2S2q!si;%&Me-@AGS~J($kk{xw$s&))6z}P=jFOlGv5OpVeLxsW^Ss* z-zI8ksQVyXwo;$8q7q5k?iOF6t)2$95XQ~hE;*PlvUi2mY=K8!W0XvoGH4!Wc;H|c5Hs;8nq$cY$AQhF&t6rg zF4ru#x4bc5d(c)tBEkyzxQ0cBN1;Dt=^wz6zOOR!$6_$#;FV(?zayo^8~sz-Eyk(b zvgi@w@DX^;uCYi-JTPXmB!D4;J*nyrv`yLP1*%uBj$7g(u+dLBo1_1WeuwrG_W}5KzcmS^-^{YFbznC_e0{V z;>3)pmfh?k5eDr<+$Ix8Go`w@8NpO*CW0mtF#{=y5N@E|(#YjZ0dtJS_Sdn3%Y`oa zE>C!(_wv2J(U{Jzv>)X6xGqA%IoiJ1ct^=kIZ_{U#J zNr}RcE9M`5ymib~3z$^1)&XM^mqzTt$9UE`OQtG&)=X0O+uO{c7wG`U(A2m_^o|14 zW1-59FSj9Y5*_Ll8?d>YIPWF$=z438#E6+Q@-561d!!W1V&z&)q{cm8MPGQCuRYh< z_PBb~8D^oVaK;q!5k9YBB2ZOXwM>e0;$VE8V&7IxQ*C=h@xjR$2YW+gOS*6m>*+5xv#V0+azMp{+SyRYiEV~rzYkT!h zTGQTx8F@+viiuf6baZD|0kJFWM3=ROsO^2V(2n$Nuw3_tn@ux?zp*Qe4Zr>JY*9l~8*wF=Ix{RcSY46e_T{5+WMJz}F>;yo$$ z%nh{w6nIP(@XJ@=>ho~zD6AQ`BFk^5--HUcFM7qI9_vU=F*h%(565_NvIZ&qtra<= zNtD@6)1)hkzm9ehSN_`EZ;Qyxw^w@p^2a;llm`%XPU!4hOktsNJJrrZC2Qx2)~q*G zioPFX*YkSx{SfILsQR@~Q7YlrKzoPg#yZcY-%rF9ADcgSdRl9SO@cqTIo4UV?}|hn z`jT_>TnFd>WQH^u(YRgGrqt!efhRdh+$`elNfAm@KXvmyG3VraGq5*}JqlWnH+Saq zpb_5DbW5jY4CV@tf11oEG+u=J1CERl1F)zA#EPRp+WsZ0yaV2Lu@=uU;GT5+u2ryc zm?xPE_dfuBaQawdqfaas_pg1AU$a)*&&eaq27H z`hB9HSNs+vYxzltKpNnXBiPLUlxB&ghfU>e>rD*QZPz!S;{uB_)$dSSdj%}yfu6%5 zAm~yQ;~{CNQ;T2GVO8n0(EJcrOa%#|V%92IV9H3c+F3*GAQop11=X0HrT~6i!&aQ-l1oc^3!f{eB5?7RL?U`fws+vxH6p`MSdE-XKI zP03C1-GxcCx_7b zKQNjYlYm-_l`?|v%s;-93HPpy;k1EU>n=_Y9fbdbyoXU2usV^! zIE=$aoptL|Ko|{!P5Tk=n>E$yl__!G>IYP})}%r!vh-qbIN{9Ee^$BR*_A%q&W<11d%uu*Ogkk}E6PTYR1fOzketrHs6xf4IK(5sK z5!w){#i$XnVBE%v$<1g(UG!e-d0Quf3NR2$(X{a(fvH+G)Uhh6RUZ}vg}IY8&u+uu z+4jQk65c?~A1uDnP1Ed&)FQD5DcX`oAwWEwhVC-JUHMM=c`>h0LxE?teu^?BP^cUMAu z-raoubps_i=QPKj`_k$<=%d{RLN5unj#l28R%3Hx@1-AfZd<~r&mfa| zR4^MUFTg-*7hZ%95r-)$=~n8NZPkC{B=7HoFsWS_!sx@!GQmXF-JII_7g?I#z`x%= z&Z`wT3AVKz_wR48n%H=R5B8`|pt<74p&G2*yTd^q`UxiHQSlbYJ?AxO;Fb!sukTFju2)hJUTT%3?~m) zbRff#?#(rUmY_1_tuG|bvKD+ZesCi&6JFAJ#h)m_~F2qHqfe*nB( zj?*&lq95BEOa#EcIbYwCE)o)H4Trhb=X@?lQQiI?0J! zC%l7tai1KA?8dSPG&^2(?;2NTcB0lk0Lw#lX%DpAn0fNI&-)-{zob2ZcH4baKkwe3aguPahshLZI^Y@rYML()1^@Wq7`P2kfHt*IT z&^4GE+LM7fyA9A~vGPQ_>giX`S)$->e8@W-M^Pwzo@0at4AW^7J1E!tJY~F8@sjxHM?eR7FJ!^qjCz_ zpyIMFNYzuqV^aAIz0NwSx}i?fQs3qXrnnad$p$?78={!%J}Qij1IUH}cLbp*!}Lzf zD#pL}^n|6W7r)=FC4pI8d}B-IfH&Y<1Itv2RxXe5f>pRflWM)Z!9e&986=uooX3@v zPG;&?4}gh@>4CzzGTOy3L$;w^OB}^LZt*;)FqxFIgF7dP2o}gpi05yrEh5r=qEzWX z%K-k>@6eU6|F%WCQQxv%KLWXm`6ODbLV@D4!7)|~%yz26uFN4Q_SFscJHPzmD3967 zBqD;`{op;I7A;RPaC#!o?KIh39#^WQAA!ogLPE=+kLuXPBgKdi9n{?T)Txx&YhEn5 zy8HtKRNanLFK7+K%!W~35Wg3qT1TS{#uHu?M5E$UJ2qdnOC<+yvT0`!YOw%+*b|tM zE8rojT-cYBqxYNzcBU*^SpWK?5+T)ho9CDX-%cVPkm}thjN-i%6#Z{hPUt{|VH-Hk zxqP8iIC#b7S*IHt>N!l=y&ls=O}LV^8a}FXr$r?MJ@WnjE2rIC3@G_jl{G`67sq<$ zO!X%NOVX?L54At54W86uE!c%1TgzCRrNA@b6)+X4dS8jIw0$=Il~7OfiZNw>Y7;!u zMi>lx#98}oa~Mi06y$qTkWOQEFC+g8ahJz;Vb*lUK8tYJ_b^|@4|n8$nom6sErcQa zJ8hU5_n+sG>EpaX9V+UI{jZ z|8p%%QN?uOWm8UMNMWHMJETR~C8gKN%8JJ&94*oB}a#8LEwvjYeG+-%KbSm?|&RoP~(d<8X zI?dLTT^;}NE5fT@iHPxL$OfFjqR~6F6w~9!o6Q`B*gjO*J8$51;zv+$v|}HE*D+9- z#&IT~5C(7g1dKNn!HM6XSo#F!w-?qUT&{D7%+ksq^K)b9w`LBeqt4#T>y$a;qK&uM z>c6bz=^S%G&cAJOMrpIZd3S6C12&^RL5L7aY@|TE%avGs>nB;bo~ObFbtx{7vTNx1 z>A~el_*8chR0ax?QaJUEi2rZ$@c)K8|Ak{ku0v^h8*9!7C;1EOH9YoSL)bKQAD36~ zr%^%1L}X1_TiW=MtdsVk#_ zAStTUV*)#7p5gy`0;pW=s^>)t(y`{}c?iVN&E4uB*6_Dvf_dNF96m0~AyGn@+aG%| zs~R&>ULu@>&JtGhAFa?XV9JRv zO9uKll(_y+4i~@e%aNRX>IugnzO3pSG(0UpX=l`F%`s^D5>R=Hdh}eXZLW;}_x=3W z%lbcfSTbCeD%a+^*!y7Lz7p%OKlr_~CFxMf;DzW-u0GKiHQ|`Y?%B0hC3rFOVegwq zhGoV$NtbwkMdnWb@`U;fITtdq^i}ds(DKjL5^O5EB?(d=#nz^PKGbef$MvqWUV8%{ zNSPB@4PC@qymnN*$5vIHICxY5O&HkNqR073qD6zrNrT0IJdlXZqSk`B(F+tWEM;CC zl&CZ-SA`g?8$Q=~x`KLZg=0cy2*s0F{SMcnh53J^A}EY&&+qH=?x){%F;(#Hj8L%Y z**%&~mAFWglf6)K!rq=1t;VMq@?8$ASE1*v44%)l$Ft4pnaN6CnCG0twy`v?9gLt@ z5~0UL;>;~_4e!rO6D{~i4j&cC?bqJfAVB1XXQco92l&osIE;!>$6+mNR|>Bzjanrqbq#8vNliAD(AY6Za39L65DLKOOv7F? zU#s5gz7cl{*e9prIf(zz(N(d^W~) zEO34#xqMF?VPKEDTEByAv1j6}NjH9FlMB)=DLvz2X!udz|B@LZhcC(f1xbz^Pu0@SKs4uLkm4pH))Y|ehhqXWIBbdy0>RSeaIGS(8 zKbS#b0gLBKh~l2B3pO=C2!bFnb1Ax*T-~ujL`C2rdm#r}Az!ZOm-KQ*!V4BXwdmz@ zYuZvfHhbk7!x~=!p5RSzwB2?WSZ@8(m&1p3+(;OUJ?q^}BBNt~5bftOt!Qj#`H2RB zW9A%%g!Io!@4~xyAPxE7mjGQ5A6sbV^t@ ziqPEug@ORu-8!9Y76R9a^^qZtBZb8DI0ltt%{hHhbTvl;Pd_v++qACy@+!FoI!>hL zdXF`fe(^S~M?`fwn?%o9f}!?T46_ifasS1dgj1J-}| zvvuZ2S)GlLFKKUGu2o>Y7Q1c)6s;Ka!)3JQO3s!vg_A3JqZ>OjNEk*x;l& zy#LV?ds#CoM>;*5&xmNuk!BzMTf6BrtTk+2r};?s-j@?`v^r<494NF+xTe63DI}70 zGTwv&$d`E+QMo*s$?L-g$GpaDrjY}C(HkLB2@4ve@)1D~)>_IWWZ6KFQT;FmTw4|< z$$rj7$9oivhq7ZUe?t^M+YF}gQS10X?#6eRRxULXO`+MV4_dsvY=EN;hYrS_`mK8z zICt# z%0AJ}VWRFBTEPLz_`Od?IU)X4&XtjPaYEA9ab&B0O2en+jQRJ6_^A1kkC^-Z0+X)B zJ~YUgQC%*6PJiqHo`>Sf;^Z1R-qiJy9b#89oANWFx#{74Y=!Xg9ygqsv`eY;g2=!| zjK_&I2WJi6Q6F#eU^3!osIU890j=Vi$(~>AR)hO!p{^_15!5~F;LXMIwygHCs6e{M z`Lx2|m-ecC^{{doC*H5tUWlZgtC(9@Hp&k&USCNX}q63rH_Rb>6RSq7_<%=O4jGq~rSJqM`%{Qsk(rbPN zeQat$W)Uj}R(dXiFZwNS7KR3C&Pov_)z4TftH0b7ZbWf^LL;R7;ckaTaYdn;0wgWb zQ#Z9P_=y|_>+?C#Z@(6;@{(jYCDf~#VLbF#p}{3rQ$@wL1J{}UKb;=6c19TJr2JZW zn$WF%GnHMRrV1seqeA77|B5EjAcjgoA!A%i;8@3m{^$8q;nOnPVKY{?B>O;+k2*g0 z2*(j+n{1VG!i9a;9{eydnK9Bb6>DZ`C+gr`|dsvSb zUwDtlXgB>5dg+&_QNU|Q#j&OrqLJcX)Q6NiuFkQ=A9ZomIy8Np#S-X#%64#js?Tt| zvA;S|YGtynH?aRRVkxhkAO|$ad6{8k){zNWoT8?w2HQvkF0uQ$+95{+_4j8jOVA=D zCCODJyr^VyPL<|eh!f$*XpV3qC;147%m^E6)U^F0xj_<%#`hS3Y>x(c^T(yU`bsJ^ z6pvRetmFfXiLQ0mv&ky4f<7A-OlXFd?b@fvouSLxKF?gxHSUgo4tef7yH{Jt!LE}D ziICOaD?P(7sQ4}Tjl9!H`bcK`UyED2lhW4<7erl$BK`8TtIWOlh0gx35*ewqb{DE& z&dQDdGy!ocCNyx|qgK`z#Zb=nR?4gHXYUR|+FcQfkzZlBrr9Z}gQ7LuTe+lL*3pW; zqN)!2KkV~&-R$0s$X%{XTlxA(Tsm)-GOVvgU-|4B+ih(7*I@aknzUv~9IR$LYlx~V zYc8mRX^O2_gFnTpPRN{}72M^|F|?ho0xKKS=Z(_-C~Qunqq}_!~v-Ce}J?5j4^a~2ztXY+FZt}lv-zDV5CI6 zls!a-HMyGlHdGF3AXmM{`D&qaOLaNWtsjmzU@(fyzc%0odO?wh`VDd0-~A+&vbq-4 zdE={l;N-xUQS?Sj^ow>?PRP;v>UCC8d)_(jLd0$rbG5jO_NaSEVOy{k`t<@2E62C^ zdF*%UtTf{}069sWQi%uy&Ek)S3+r>AoJ6Z8V(}gXY3s^OgcBY|`f-X04G1zT%~y?A zQo=V0?q`T+NoI4mOG#9Sm`Y!Q=rbvD0U!iM>l6e`6sex{g^Z2ebX=XKNo$)qf>fK7 zGcl!NzYGbTIC)ketW1x0NIg2D2OuS#k#qWHD9+lE^f&j-TGF3rQdRXD{PMu%3*zs) zBP3tt9FESy;nOchDW{}M0t~*mp!EI>sIJB@t=&pCE>63x<{xsRfAY(6@nvE-P5FLT)Rb_~Aek&TOQ-v6Lxs7~U3{$+xOG`-tA} z%3c>-??Sr7$-oZ{!k2J5xO2scuTZXlW4rdk^Dl>gMSGu@=Eo-6#rE8y40LJ82bd_Y zBCfVo_U9AxM7w0!e3$C=Eufyecl;{AvVQ>KoqnyMqJW_E_2<%m-CT0v-X}=^o>NEO z9;!{m%+xgeSKZdh=i2Ezx4fznVh8#0)=%a$lWlmdg8WlivVydyG6PpS*r|gOOEC zb=|z1VJb*N)`<}Xi=m04p-R-YEwiM5Ob%NavzeADAwV>@RV0_SN5x zR`r(W=W}8TyIw8f{1Hk8F9inSj7nup`$G0QA_m?K=H@+>Ge8iX!emBpFLCrH)>9Ja z%(0~7wf3Ib`Ll~`3EOS5wdwCYdh5T-W}jtP zn_1C`uGcw&q*q%Ts)fw_qkP^oN1c=E!VHHK+>CBTBz=+R<)}pC(a;3`@a@BWvv))Ou(kEJ)_h-D!oy-eHJ_&CT_=?+~x#TCO_$Z6Rv z?>_+OeZW?6>i?CbDD^>I3it30H7LgrbI8!K#O@d%0I6Q^Fx=h;ZE4pwi7OXs)%oft z26>7$A;ix`6AZ+W>zNrh9GhM$=(t|V4%~`zFX9n_yTTRWi^q>?(;xgrElkXFKl;Hg zqZu4K>t#J758;8Mp!AG2Yg2(1by<_|LlKBgW@R0UX}9*8dbp5Q_CfYefS}QTWT^(l zsM0^9CrAX-d5;924uiO+gJb7Wyq6RgF~1LYwfp%4zxw41?tOh(|F-s(OAR)qyp4aN zHuN2gFcxd86R7S+#ECmRK-N4x40#6Y6N{Qu&A_sQCi66fNtB$C`m291`$IBYxt!Cl zDR*9B>}|`8*n1qVfFyfOpG$f}=+D?`YV{eGEh~Eo^?R6KEt2H|y>mHyq}0^Z1aE{L zpWV>UQf*V-0MF!*^IHpO?jRW_+z-~AJfmXLpcxvy8!c#(x3oYzf{)}2Q~-?Rmsv_C z%jI2K$8PEyRZ{yPYQc@{;!atWU-T4D)ltDY0j4yNO3P$mX*sMbjZ)h&#zydcEoxtb zc6FUnot7Mpsa??jc{=8ndD-S&7H64s`ix7z^(nrGfnO`|gF5vG@D(_7FxLM~ZcRyy zd@o))%0_pR0)^)O1eoQTtBR?z;AH}Sc5)(?Z5m8?dG5 zIAxm*Rt(t-g7P>85)FV%LGTR4h(Yh}px=~YPHu5@WZmz)n9}$}?5%hhyZ~Flha;11 zH)*CUrz}tZ0kkUq0akSs3=j#i<0`qeBe5tp4p~-D0J)joKY)o~iUZTQeaLY0s)&zK z#bMJ5fui5E4l)=7zGHt7-HSUj8IVCYPfod&28w|UICZg@&xQ2b=nBJgBn}qbfX*{) zcIwv?LtF`ZRgs;=@3{sC0-nej`rr3<)H229fNI}4n&M{?Zq!&?Gf%}}Y=r)|(xb6t z-8J8;{G(i9q$cYIDR9~mcw7XOF>z&Cfk#ZG!jW?*6dg4V8|x=6|1@9i$N5ib-+B!= z7h4T#Fn_k@e|JzB^)acoDW*LBc8!hK(*AZnqVAL>%2JQR?_K3;*{wjfAz0JlF?%qX zx;i$JE#~{!XrgLYsHsEprg}geJ*-E6 zYYlZNdIWI^9X9|&fWbgC6(OGW<8t2N-0Aov5fKr2fo$w%fx>;;Dt=87FOT{mH4BJO zG@m?`D!6W?^gede%kI5=QGG)$ew?)*3QExOqplx+4dOUi1+t7|7hmk_y^r_$RWXV+ zNbvO1Zqz;nI;vCE5NKsc!jU-ze;+_kA7+V7#?88vWaoS0e_8IFk{o7xvpxLwox51J zysTH(P52d*tFD}AI%*f`;n8w`q5Lrxm-8Cbo#{tvkZ8?bm0wml(O2r?u!);6ImLXY zJ;y&Y*MZKqgMn*w9uue<;#qdN-qSI9c-xIufC7SiMhqG@o0|!VCmW>6A1TwSr&r^@ z<~7fMksuB}`!drqyMG~^qAMgKAixaM#@OTIW|w&nH74KRnty=zQY>WMdIxL#3@#xR zs1w=E>Qt^CX;^*Lr{&1TQ2k#OwexcKjaGBB_&_f_W+_rV-mA@h1oE9@9Dj6bjdvYJ zi5 z)-bn1zYTxf-qb9$z$ZxOEl>rB3Z;}$dKH5nn8JA!1tCaI;1>#(0#n0YK>ey_*a7-i z7VOmvr8d)}4S1E+QwcT|fP^!HLtrp5LiQa{#UeCB<=SNu^`d-I_2($ROxGb4hYe#M zn4{*s=7OwhrjQ1#yh>PFygrMxyt!60WYH4PpYdt@VqS}~D$CH7+Z~A};(A4c>6IOz z?mHp62d3RSf@4|QaS`!lKML``QSB_ttxn!rWWI}c7++~^j@$UdP z{}2B#Bk=uqAD+wRppgRfb=M=jgv!9+#tFX>e(y2{i>7s%pL-1O(#%-Q@Hu@mpbyj> z?Z-v5VQ2DR{H&(%hYUeTvwl*y@d3ZTX%(I?0h`&Zy6c`lihOSD+FKwse6umw*342r zv+0OK57<9n5{@jNn6nwMgSYO)TyY75lHHA}=z^XNID)XIUCqV%Y-02jyv;+iPg>hO zcr$;ugsm^=_EK8F>E2}`z$^!jXaT!6ynl9_BsXLJ{iP+nHwcIBj_Y) z%)4DRYZy`Lo7CKmE06&IPp(Uz6CpE^_2MWG6%BpK)HO~mldg68=6()n;_r=nTS~|@ z(vn6$cX*}%xB})7ikA64JZ@u(Cayrbu6i*7aGznS(kpbdK_L+9V{qT}Jo7v6$ zfVf&#2G07>90o7N*CG@pQs+vOEYRRJGM~EaEjf~HwQxjwSL-Tr7OH}cUvMrT*GBiP^q O;UCcT|2KAgpZFi~E^e0q literal 15231 zcmdse2T)UA`)24xX`x96NrVJJLFq`97J6vXRHT<6ML;?tARr}#07@?j5Fm6CiXgoT zNQVRvr6^sIq9U^X{@?7*?#ws4vortM-Pvi{Sy zC;)nYAHW3}par0zqNb*%qWN3V(9qD*U8bk|TUah#Vz|u40tB+Lu(7do@pH3t@N%-T zT@ks$Dh1tk?VfQF40C_=|BeN}_Q%#of`M$?Z=^wG0haQ}2} zc`NIn7!*3I7-ZfwJp;aWTUOi2Ik2?8dGg~02|)i(KWc#UY#1Dkh5V;hnb?rCgVIC zHV*(6v5Y2-{^1Pu`(tYpwW96WJ{OtM^eH#EA{5b#`KWJkM5&hL`Hk5yB%S9q+vbPV zGY`LVJ1cmTxU;yBuLfddS&iJ9&5kV3)f9{Ba4q035=- zO2;w;IpjADLN4n%HY@E^DZp6MQ6&-4fSfX&?NdWgQ7fOr!hLlGsLvw{il2TqQ9ZTB zYNG}>&W7;Tp%YU6btQ(Gi*h!Utll#O$iL6gYS5V48b^g6j zdOYCwt5PY?$(x~(o3uggIP^8qfwdKRR)s=0?{IboQaKz=%~wl>#xEi)hgt~tjLN?@ z3BdV4+2Iv{r+@0Id|rfEi@`nle6TI#8tIZ(+tI)8d8Xp4U`bP}5~d(lMY_SAXFmVv zL*_RR*n88aIItEq8&W&cxP@>5SbUBb3K`TYTy3;gWj)iQE?X3vpjuZH6*Nm9RBNW!Oe$&Sz*MeI!~b5)xD65nc++8xMvsRv%NAiQ2A2FlCXy4?T6?vZ0} zNM9CdVLGjaIKzCC5#7vr{r91NWZ;OB0Zw+VB?VVw>}h9jhiEv1F4D}9qK@KUb=?3u zbWvv{M_J+JaJHTEH7<8&t%Mmp*+687Xb5EER;c^M`P96h0vSgL)0Pmy36uP=(}VM z#)=42r>z63j|IwWSr6(2a_L%0KF$@MD|v$J!Ap)48y(n(XUV+jc;$DB?5nT7)ig@_ z09pF(8fbT^KIqV#e*Vj~Dqr!vx4Sz%!A_R*(4m3f ztaaSP+$8TRy|DEYLA1otfb~U%V^m`9Y~Qf(NKvjtj1sq%!JRvtPBg8C;4altQQ>nP z`fA=w#`2h^XiE6qND{9fk2D`U4uOW(LzV~9cGCtt(i{<8a_#Yhvk#s z^lQCrRut%~+xZ+kOV@FCuNQkfi7rm_jrem<^?qVvmdzSg4j@hcn}1{FuA z%i(OlRlN?h@K!HRO|Cf%x0YER9$a%6H@LOrVlu21mkBAi*((^_kMB~Gx|NY(Rgi3u zsa8soda`u`mMvz|6>_@biyySA_yV;>hrkm-s1VR>!@sp#HBKjGop&C6(J8aPGm`I{ zDw%Ft*hd6a+=@Ndv?~3LbAHd8+J8xFp>#VKyWqWld+AFQ{Roh-nA$mci7JYvIAh{$V6+*h&IW50Oiz&6zAmsED_ z!^RK?JhcOt#h+nK>A?DhZsXspK!(BQC!W?XPHF#oGy=Z-lIf8CAWGG8EXvFA0#J~` z^XFaE*{_rF>ezB*;n(5PypiYy^x#PASGB0kHcN=9fYy{8BCxbZ$@ceP=;h5M4t2g< zYjjJ-E+&sPGIOYND0RnxS60`s8Oqf20aKPh*y-WoASakEFN}v{!BFzB`u*u-kac*a zo3+}pcG@|(If_KaNtJnmJfsRdnG8#ih^)`+1eQ%(QyBa2*=N6ZP)FtEp$*ILjl|sI zaLx^hk3hM3?3@`Li7ennL^g-~XTsWZMX_WS&s}+aIc0q3xlap@mO6-=ck|OcS<{q2 z#DgEk+$42_*rJQa15&)?3a1W>$x>GDalXL#b#hnUPCY!osEI(9kc_6F0kN3y`Ssu1 zri}Qr!|POit3;N1?tL}*Q9W6re)d;l|>E8ldoOL%9HPXTavm`{7DJnLo@)| z1_P_&gf~Fw0%@eLUiTfjrL%@&S&<)b;gQ_l)PvrQ!lGrrLM`2yMCmtW|#gPg{JfUk$;P zJTH~%NfBFQUv9^*Dj#m>%gQo=Mk#*Z^*k#dnz6=ux(J+nAulwsKzyJ2oBljdjTPDk zRq!im8VexKH7G?V<4!c%WOo#{|Qwsb2H^uhoWmgn{$b!nWOB zo8Me{wU&}iDBu*C!iquk*4LmvDLrksWbSHkf81**&&CvI#&bMd)O^#1-!{RjepIv` z;$5n|VHfB!3ft0iuLqY#bUKLI;5|c>a}lRAx_(+qA|hUcohwhrtYu?uubTb5{CrOQ{$}OYuleFYcI^($ZDibADyhP@ z&ndV*)U6oHZDB{Tyj8xjl4JjTS?TpHzUBG`1{tt2lTmV*WBFWl;qJ_`uV9K@Eich7 zP_2S%O!K(H{zhd}pY7GAuYL(e@yOK`W*`jhD~)3I!hN~j;W8Z>VDdoFBKPdQ;p<9L zq1~D%rp>R5W9g{qK+@UMH+C+=k?RgW_%goTkkLdE4vR(zim~4NC{}4= z6u`C2o8&+0O4NJi8S^SdAD%6Oz-xTYDvoJx&)D2vv(*$=ehZ=UI08Eb`uwUsvS!EEO92>51a*xXHTbF zc|dya+*0CS&Q&^+g#}{85B0BMy$S+pztMhoYJTez+}X=ZSibTmk3-V^ z)%*KXDzR1mEZ;=5NR7wknhr`SSkRV*r}#;(Vb_QLXOGP%6zjX{ke%C_a2TAPT&VKh z?hx;w;SXHXoh1eFh?8F2o38PBybkBQUVhfiX*#HSescavyn~ zOAAvldQoQ8ZGSZ=h29{|9Xa6pb)BrZh9-2HxvsdHN`!xB}e&Aa8HZNQ5dtHIW(5)6>TsJpc2SS9AF3 z^AJyilh0wh(qO-r!`+&5&*Sz)hhRkE3jmPgP9=W!Pr=|e=#CoDsQ(9j)5<~{WM~~b z2qamdNEYUbIHClwsIr|3W+33=CVr;o7G#lRVL6nWO@IPZRGaow4^ zGV&Tckb~OEX3)4qk%c-dgO?PgDkHF6Bhjk}ctis54ue+1FC#>*w1njHyO+-mZ+vYO zG#FPJ)fdWS@Nm~NiKeBcar;s!Ky-H}8N1$HmWllGL~t z6qb7y?_aLhaUZjMKhS6wvd_x;Ap;pJLO;>DwCcubG5=Vkv+2u;#SsR_5E`0fto%&; zhwGZ0!V*WTK}eMF2&@I>*Ic>Atu$I%+>bo%mmJQvNu>$RV0-m3w*ClVJdAKrWki`E zrEi4A1qoP8pI*s9*@{3x@h*#+6fQFWlRox!h=io*`Qv%zU!M(|aaYNf-ykn#ru{rO zXpZHVv>wZG6vF-GLB(GBJan^24bJh~8ewJ;CC3rNJ8QtDx`j1cQAGg;19!h;?qzd# zusw!W15}llCYUb&UjOoto4IN?iZ4AR1X1@XkuBGY%{0?_U8(-coBGdEF*TMM%2Nr} zUV&L_j5|2+_RZPJqk%#Qg@k#y3x5cGKuhukfQA22$uUOfNMZKZFYvZ(E3PgBx=Nnx znIk&Mz%y@`)KR|lfce7f(a+O_=bAv`p(eZa-mUlf4g1Ku2lvYx8n|A=`L)vU#HiY{ zL6tBZ`A{Ov!Ef+;aMylMZGR=#c)y zcQj6Fv|fpr2VXKzg3rA6px)432(!RHH~g?0e;(Gl?RpSA8e>Y95lj-*jp?Bv#maHi z!eCjAfO1V0d|VM<39Rg<Q@aPvgU((;#Sfs#Qf&bj-Mtf3L>e3Dn0_u8`mFwHzBt@c6mQe_&#@LVAgAV z`z$@ryc${b)!dU2mNf8(CBy2KS?&v_3Or;(Zi*6kKg6f{*||$~rue!3wML7e;(F4- z*PO(_)bbpUll3O}HrdXHfFH88is3-9w8*vsmD#Q~a^2GRh&Rn5Aw(${r~oWQlnqnp z3a^u_x(AQty!z*sUf-%tlROLP@5j_NmwTpCP}m5LP+QVUHjhQkjjFb|Q5kK~AFkPf zSk8-O)R{V)!WT8!nM-OrsAAp;PgiB_4&FUTzUk^VYvBkcNSiU-)`sQQT2)cZ@OOF8-4i*Lga8V1Xe*-Or_a1!pjyl`M z?=ZhSj-V}n@z5XjUIM>%e?b!EyMSh zu|4VAHiG@Ik+hkde5`d~YSuZHG)*MGp|4Cj(w8oNgTAWev$xz10ZA~BIBvqZ-aJGW zE7D$ndh=Br8(NxI*~262{Y|upD1+WF0HBug!(|7~Ck70L{W4eVa3D?s^J*qVPo4fi z0MWlbWl&l6o0C6dFvq%+=uqW2YO=oCn(9eIU0qEeONx}`8X$^A#5Hv&)mH&`;zIxG%^)P+MLXzv{M613VZtX484lD< zjbCiGN*(=)pCsG^~;%<`v%!L}Z0u%a|! zB&LzesXjaXHv$43oK~pKlCZxdkV=D^8|>{O4)|QRBH|nfxq2}6OrS9kFbw*=6ZCJ7 zR{vF%ElsQ*e3Y-K4vji6Kr_hHNP-mE%I02J@C6>k`&NrbH_0@>)P)hpV*UZ%|N6%~ zv=K4S9mXzq4*XN-;?sTsFqV(}h8XZWzry5rzjS?{|8hWq$WdTbRpkspc0yk#V1Odo`J(t8wP<$xT`I0!XTrWi%@ zblUns`|VO#>s$YP-+kr_K#)XrwzE=0NxW9KTDv}oB2R>YN0`zpxZM18n|qk{oa{>A zipq~qPBNQK@hg5r|4*8zSga_8hA_n{i&%!*|7LU_*OqIzB6UZ*ZkxWLxP*ytYwu)f z!jG75r(ps@hCTTH6@?Q#EV*$TV;nRe9eb#Utt+iH*FI4 zApG*IkRBedJ-f46mY_vwcfh|?;eM|cv{q$XklA#$#)k6ysAsAY1Gx9ry?NyO)Am-c z>CZ7PzPrcC>kpkyjq<``$~(}@hF+4g0WDvY_fp3~kjp(EnFd6Mv$m{V4D|Iy)Lm5A zhv)=-!z~7~qK|8OxmFIHnnX(p)6y4!6360R1ZWiqM55Q9<+)KQ0H6FrIu#Dl)0D!(Q<=wb{ zP0NUS#hHm*t$n=Tf=!d2BK-uP< zg7xx`tosf$x>S%u8IzoanYMZTE!^blhhv(~#cO&r@fao6a(_KOjqgJ$ESc7!h-~FBItu7Ds2~HcuS@Kgdvvg{6)b9QY_piut4IBy3LT_0}uS z0V8OnGX(bj&2outBxitD!XvG%v=0oU2uXXYBHM_T;W&POdEH6ZE6NQyB@Eid!Fp8;NXZ3nu;Kg*K z^5vTGgNMtJ{2$ji-`6dEi_v|4RGDorn?G+zDt+r@=WZyG-#pqoqG#N&Og4MrCz(#@!~%Dtl?l&(_|B}D#1_bmlsDN@-OZAzcP!l zn4O6BW1>09wo~6IzLyzj&D}7-E%}wmk$>Jn2F`mWf0h`Ji2Gv^f3-izB0JY5-qb#D zkOPPYUCsoC0~E&|Hrn&GkrrAS-o?L0=I?<|POg`#W^xpmr^irq1CpctX_x;)_$0Ir z++)s9MT-ZpkhBE1P~YtL%$#*+XjFdh&{D?gCUFa9!Hk=;wKjlYwv~@AnFR~HO?o0M z0eTh_>C8(CHnZWK9+_dbFD;o#Z-t|)#Wo6{=`(n5u_Gt=_IU2Ld+l&jPmmM5 zB>sbk1BMeP=5IkyVoBzIHg-=cqwL*e-IaQW^|GBHN7~D?x2KlWO)uS$q@kdJ=nwPh z_yg(TD9=%bY|%cCY^<_mABaZfi=OSoyV&TwKhakB_ygc6Q^3PuZuYh^Qiy)f1BO++ z;0DhSyjlF(M@*^uD1Mh)TA-`d^ zF(E4K!~8kXxGCutBr9t%gwO*rSmj~DHdp#iPV zwV>D-O}NRBKrVA7w~?7@*9K~b=7(QS<}z}KN#O-!G{wX8x09$sVYUK)npAYwrLy)&2u+&Hj%8?{4L{$*pNVXk7k8WF{aTLSoSWe;mQwte^M?0;*L>GQRmik$NOo1 zeS{PZCYTw8F#GVI?Y}7RK;(qb&nWvq-3rFzYa11++X;;y%~p$UXNs+7a*ChpZMM-2 zQ?AW;8TwN$Hfe`~Zv}wWE(H zVYO~&fUOKe#`U^RYui*OqUw*etUIq4b4fHM_-F>cu+@^pvJu*vg$2%gRtGa?ZQ{xh95`W8m>VT_;Is#fQi$18e6)YMBhNW zr&b*E?w^V`4MKiz-;D49&q18RZ+Pi7G=YHSSghJho~nYjnW>{w*ntxEi^NMMrus$t zq>r$naN%8D>|Op^_tiq#oMFxswj%u8d98uHPR_$U;Z~?t5(5-`muXz70hUMw1r^Ft zPgo*QlCV~-`Ul`f!F})V%51&)uOU%&MK`0KB5SSIzY&KSkLES8P(BM!gUCYG98@-D ztVQkX*M|cXCc^X6Fj~y8UuL{ zS4O5&SkJS%k5>}Yn`%=>*}h$YFEyo&amMGNKAC*2Zit^zVf`q0?D5vY+|woI)eh=+ zc9>yPCRjNEmZc9*V4x5;1O2X((W5PYxSdjY&TX(H8DzIMx2#u}fJ!!`oh8F3eBB^= zB)H|80?9%~`vhGf@u&QS*Bj0HhThXl@75*fM0Ck-ZYR^IG?Y;CQDc0>Q0jVSpXxj$ z9jV~=A6vbuzb1Wcgs_+=sLr4@%6O+vLpebY)}YXv*nLQYdZYT?Y{BnzS<6TcMdY`Q zfuo9=UyEJ|5Vd zha_7=BngF+;+k;P=)W2};SlxxMU?~eed2LbReSc&qNgwsg2l=~jz#mDpp0TocCG8k zv(o*?>i36t9@f)e0KQ2f$!h1?B^~I-I7!1r%;ng(Ll-0$5>jn9Pz(NJobiTa~eP_%k5BiEg~KiO`J);X%Khh)D__uf@n z6hXW;DbADx`7@d}SF^48L{#=lf3fxNnG)5ys^Z+tb=?ZU1Huw=jrFKd>BE2!5QQae zXq9)mr{#!1(oe}(rvYxS7gG6mhYil-D4w77#f1{ZQabgl%i1`hi4lHQK8;($hN!GS z4j#EJ{jXd*3V#@FDB3v-@bQ$Ai7S{pyCIrn>$N!x;W7RIJqLaz{i>=U*qQ zCR%5%T$6+x&d?Gu5}V33(wl9x|ImtB;i$oouoyTBUm>_*V8uKG7Ouz>r05a@%r!Y| zppPa@mL<7%;zuzS?>utzX7S#1@CRlg#|Nf*!#tl4qX7dqZDr$<1}7D4*ZVso#t~yn zVK*TDo78vWRx`zMx>8YgfTvxY0nOtn0R-(!%`Mmw^kYV4_{>`Oq9>@JTj)k9CL&xLwKVAcnB2*Fyy8C`T@{W?qg$iNpQHN);_ zA~HoDcDLKxC*-$19(Ue#jkbF1Zm}Z2>-1NVsPy(B-VP6qOcs5z-C*ebY67HwKa3?D{_XpbuFW}x>O$3?#HYrl zJz;)dy)moKOAjspXp-9>p{NT0&jYhd_v2zZ49b!MhPdduMUQ`sG)&bGGgsciIqGZm zQ_Ep5M_2ore6P_udD^m-&pp{Lm)GX-$!`B-buCx{^#SHw_9!2or~mi#Y0)h3iq0Pf z^QW(Qh9M)r&r2HkrLN0Mfd#UhnQEjtn~)^y32K~ii!RNbv2vuBfA5n}ZwQegiExz> zr1%web?iWH{mHYl{whyROXM32Cq-Ic$bf^IlF`7e!vaH=7v|fn`7b0ovtL~yH47Ah z?|z2xMRy$wdHy9DKJqj!wiCD#(gAZ${^R-~QJ}9P;#P-RC6msFwU>i`SDE~Ne25e@ zCUvyMaD@K}Wy_iPhPeCO>VD<6&Xa+ctxJV**X)e6BIZ%O+m0Gdaxg**ExaKJdu1K( z@l~d{s^+m-Q?LoBc2x@a&~A7C%}p9r-LZmJX(a<7Z2o3BnQQ^Y83ppA3nUX427~1&Z z_U}u*W9?Gfi&(jQ9P9p62`2e0-d;~^+u!rK=j2>mU^ZOmmMqu6(Q$E?PH_Pj3}gRn zk)c(7$*Vq@prne(D!sT zX3guK!pbNA=;&a_umqA~-mE-?YsWV8@mcd|?udbJ>G4BAwx+yFTkkMoE`@3b$+q;o&+Ra{mx#l%Q?;uS(lgz9a z0q(!@U;PPG#vLFLZGs;w(Xsl8luH}KbUaJw+xFB^BBl{J!VnTRr*+a*lfiv;~Mw|a&yX}1j}lFsz9|)OEWjGBs;#ywZY`#hi!PVXr|AH zY>ES|9x7aYt7vzCdw&6N-l^usV8;-*4PCvT1c2k% zZoSp#BV4<5Safff5y&;V5pyMV2u4t#)f)NzJUQ=IM)K^ZxluEGx>JJ?)9n$~=(8ld zXClH(BW6%|r!z|b&tUYjGnrUX5M=HJ=e$CbVkNp(Y z)s^cy`Vjo%p1J6kfp{44oU}2v0d6t-Z1}_VgD+CJybf8i+9D*Hpv z)RzjIo`k`TV zUqvy7lg%oU%sKBzdFvIXhA(|BGPa43=aT-c&GI}%O&jn0gnKNz_w;ydjbU|6h%W}E z++HU*-~LR&aLAa(tIHd2ui&J~fp?I;=#Mu5b?EXL=HQ3Cv+*=-o&?r!n=Xs&Y5!5QM@hU(|fR{ zbjaHTVP=#4u)MI@#})H@*D!M(GVggK?|R&%nLn_8r^(gpw|d(u{#0ywG|6pB<)IrO zc{$MG@aC*iLNqf<-bq1DxqmF>8) zzx1eVX*^;X`RSheul#Ybw9w8O>T^-gTx)o*dFZitXNpwe3Z zc3_E4X3koMEKqMop3wE!bi&(FUSFRdym~8ICOiYB56amhd`VG{ozdYIK{VaWR1(~BBSq?sXq4xXw2B{_n z8FE98b%mP9_xdZ0L)0iFol1ek&Z(GXg{YlyuGintSbL1+E$@V_zp=xt7~{Q>)%9|! z^g>R{9h5<+&uD<@sNAxFC+Y$)`SvBES>i}9!0#7mw$Nr-bfGDqbOC^%-^^CBPba|w z&e=U52;GOazn9=oQJY-8Md@>p>d%^e%l9~3YyHiBWc!k7bxzPv=I=8vt6Vt;RQ@_i z6c^5ii(~XanHaziW-ffO!a!rkRiFdY#mk1JJ>BN)Q6`+l>VL(v6;POSm1zO_gCE+d zh_GzLm{7-}Q&!Yx=M!Yw@Y4zSwQvvXOG$%oJy_?u)yjSveLogbz0rh)JmAyhQ`ho< z0n?+=$rKjt9dSVs4#C@uWhSo{h_~a&CE8MYwoWh5yzX)S9OX#_Cn(%}dywaF4liGB`#4+qgQnLaMg;Ub&y9)y=0vL&8WM zUFt)X*;^GimTf1A5fpBTPqM>Xx^{J5ErG9zMgf}%FiTX{RF6yv!$*<&Oo%$)7h(-_ zzINph;S-nWH5F~VGHM8>3HSR+c$OEo#&0!C5al|1;?9PG_qn1*ALb#^(R6A>l{c)K z%F1S^s;`p^nq(2p4y6tbKCrwGVDC(bqySD>)K9Z<;^P4H{u!pA&FiIcId|Zo+z+3A z>8Uw?{=jlBjD+t3Hb|49rvvn?>AShPa69xnt$8X&Czk6IAJt4b8WsIFi$F6=df&Q0 zF8QD~sJD@?ejwbnH6mxW;&y!FtnDdWo@;e^w>xdIS#D>%4^e#^Nf|x*ux0I~^%b7A zq@1mAhSAn-yrqBqZ{5d4^$NV|BBHO?5hF9cE(o`PujLxk@B32>;k4W`{O5K*NEyX{ zBH7wm4Fs>_1_8`^Rtj|P3+d~WAJq~loc_UIGFKah!&Pm)+tQn^V|=JK3Nt+)Ech5@ zIZ5^VKe_~iCdVOeMBgHxbEaSf=|O_f#HptbE=|Q}SBmH*)P35Y7X>#rS%L>I6UX}j zMX^fFAh)5%s`S^iiY76#Q=^7I6X4(#gZbAC*613iQ4~D#-u|;RtbtLZULyYW%pBo) zAy<#^^1T(HJpu}K`B=c1D{)UW#VY#qCDDzF=E#&f8nFu6?v4G!` z-gBTZnZnLAEJ8ywf%E;_3es)CK&EZYjKJ9CxhEAP2F88vM3K*0g|Tsadk_6Buz;H= z?P=Sd)@AQh8+@fj)EkPlYmcme74Iq>J#|}s%SBYXkfmMSnMAFW$Hc*B4Cwm50KH3A z#>ao)rN8Ejv!L0QsOfodb$97W*^?~Z)bC3hh^BH5k_mw#r1wufzVb&C)?M6I6IC2T zKXkS8HW?xJ7u=UD7G;3x0RC8xhfexrD*3#=8+LbVS88)Q!&-pRx42?`1W`mA6zqIU z0PM!l6LI&U^?1!)*Fdtq4RGPd=*FQ?-^aV3ZUb8@yHEH?ngQKlr!5D88Iie}=?d8U zIyk1{ZytPd;~oyA85+>KxRR(9PPWjIFN|Vu&GgeMG8sn zQO)VO%~h9yFa2|?FXnl#YZC1w3@%ePaBd?o#Pq*_v6b=Ka~q`xl>YmnRhIU}Z`sIu z@1YVV1UtW*3m%snTj|i594(vEG0J(-{?LiXF+0{ar##G9pJ8!CU3xTG@_Dn^X0T*G z`|o?7-DOdSJ~1nwm$u?@Z-V-o#pRF3-x~$sMW>Z9QgLHj&7=c?R^ThtDDjFTOe%X0 z+tLiSGq%81v$CX7hu!}~u3>H3UezvG%!3y##mUB#pkjP>!#%9@3_9T;%8+Zsz4|Le za;bVUDogZsmL2taEC!%>PGc6@jFpwE$#isGIxen4RPOO^@AzA6{*-XW>AUagQ}88m z7DO>i+n(#EOIYo`ufua|bbe!!S7h#oi@rkOTAA@uLi(_nFwv!o3V^1CHu3)pv9V@y zbM%6HoZ?@Gn60dkzd2)E%+}@V){N5%X0c*|Yr_tDreKA=(ZdUH*T) zh>me^o3cj>`F1E~#R)1&I#9KobAk%vu&Rey@#6|3e3C_nQ8f3ekq+^R#>0V9f9Y;; zK{zdtnG&9Ujs`RQ!#;B@RZ8=j;m#NG{~P|f`N*)4XW;|8cmwQPT)&R_;X-Tm71(7Rbef2fc8P zYMge)`BmL-5JlsFf!#x$%P9c^w`?0RF7JYUXXPK_Vcjg)kiToB`bW*GUGgtUy-h1kk#1~2l_M3V7_^c zL}8$$r`fL@QZ38v=`l!ou2xFbh!s-fk4^`O8uOQyLkkd`9a|t8n%s9N@cypKEZIop z)rY+_Rs`Fvx0Op1A%^$xV4mfPlv2^z8a`m?v52%;0x75u|5?>}2^7qov9I~4=$vwJ zs#Ym}n)durzDsN4N}OWcU<6OMxT)}Js>bISMY^1BEL*VayBXz=%wsRl;$X06z$v6$ zN7~=?o~4YBcF{vU$#^rqrpSJWgI^7IM^h$W;=5UDZr_XZAy{Qyp>dX&A?rUDjg7@r z=MU00Jn5v~SUR?qm=vsANK zwqqLcOQ&V$I}KNHuBbri@-10$R6k|#j)WVf{u^XJdR9rb7l)BK-jry#WTaxkU3X0RLY=aqZ%5pMu8arFPBu&jSz>oRrd8 z>Ph>4upm<$kgX#;&q;x*2Tl|7-whV8y!_Y5j>$~=c^LQR^sIHp;`N}K-Hd!F4Sx0^ z;knNI>63OUI6J~o&%Qit>b_sU!p;r}$bl6{`J%DHPyb0oW&H2j`^Ri} [E_{ij}^{short} (r_{ij}, q_i, q_j) + E_{ij}^{Coul} (r_{ij}, q_i, q_j)] + \\ +&& E^{polar} (q_i, r_{ij}) + E^{vdW} (r_{ij}) + E^{barr} (q_i) + E^{corr} (r_{ij}, \theta_{jik})] \\ \end{eqnarray*} - \end{document} diff --git a/doc/Eqs/pair_comb2.jpg b/doc/Eqs/pair_comb2.jpg index 206046000af026c893712da6b2073a67c4841aad..c73a174ac27c8281c13d42fe641b273c77c0ac6f 100644 GIT binary patch literal 40980 zcmdSBcU03|^Di2jQlvMLE{W8jB2_vlgrJmwbWpkh5_(5eM0)6<7m;42CRCB$LFteH z(xoE^h#(%G-&*%R@VxI`XZ_ya`QzMVt*o7u@7^3}bNX|bZvt;}0NFS=IsSYCP>>Qh6H*Zn&;tIn1K7^T2N0dpAOQS7uPekP zSBVG-NY8)Or2!BS5fKs+QePn?Jts*-LO@7F3^*rDOUF%rmz05#>87}%ww}I6e8T7$ zv&a+k>KX}W9!VwG>lpMf(89eEbT2UL`YlnHDlz2;IyUbsTK-G}NdF=x4>;?~<%h!J z|IdEtrk{xVU(r32a%kJ5IIaC`kjm}E78EqC+~-XhUVQ;P9M*KX7b5c z{%SERB3*8In z5>(eS?;RCu4(a^cRW;DtQmr?+W%6EmcynN)H1>gE=-QTYwp%5xMLM~G^V=-K7}fO_ zYrtjmpN@}EtvBPY6Q0SPO5MbQeirF#+Gz}y?@hT&3p|&jK$S`#uwz*RDF1IA|Ixba zmyM-j9v{=>B0orf?~AgGv9F?g>FRcKSP$vKV3kkXlHOaja}xh{dH)+i``_8M|M9ST zLliwDv3cunr;r|SkHS`lf~r1}_LB0EuK|6LTsc`kIG}ScT8Ql4qDLT*1rK$cDV426Blga2EF}{)&-w+Tt-;cCdeD&e}4TqplR|^=56K4fhG*VSW*6PISeB= zJq{XC{gQPr%4Tz}u>MK9#=2BOnj9R92Pvc}8Ha4TM8E$v$_nQVp=~uWQe1`^k;KmN zip81NIUEdjf447MXa==|1ksEG_=saN02PqEj(T%^RTJTApRn-(Pm$#QE5cv0L3yia%-OekTQzQ|xo@EweFj#B65Yh`S2z2Da(0 zvetIQe6AgdIIQYbgY&6iu{{t^`ML%QDvQIX=N`DVHWDOW@MLp5h% zraV2~WAeMBvC$>DIJckgMT#}FqgU@*D)f{Y90bi=tK(}ylqI-s7X)Ub+WA1>+@`@c z*WuUA!rGM0(aT0DX})CrJFPwD(b5YIl8eF^tm+>?iPJ>t zj4Zb%jmqia`kbSfzv#XOIw)PKm$9`RP1New)o;IjZn!VW$^S-&Hv%V9JZx>w9{lbr zn+quqH)4?gom(GJlYcB!D*Xr0Z(qKAOj5h5J$;>v8@Mmw@wJo1&*w9` zXkn&a={<k@0R~dVggD`# zG@oS6DcqMFhcf9-4|c*o&dgQv+-~@|r!^2eL(Li7Uiv2LJ|lQnVneSZru=MY*%X=$ z=IlJpo@$U-N1|>dQCfz~LwW^-Z zOIb7qQvU(GAk(Q}#-gG8I=;sw?s(|zRxt|(w^S$)+iyii!HM%k;PrnN6#N085CZo! zWssvU;Vi%NX7t=UMXByZ+I5!CAfKKUHNGxL6Flh&!*O-tdg)+IU723p!}-t9?RRei{a6xVZTz0U2olf6C5hjxtafRsHox~mD5He3!cF=Fus4UUcb#H zp}%fhW(q~dE)Wo#T~4m8aCL7|-#b_3`CQn&rksywTf!w9kx!eC6Dx<EAYU)C4I+2;x*KKfrQ(sEDC@Eugpr}!*6(hYX;wjI}vzZi@d8Dg#BoK|@W$396v z|9a6%oSz2SV$_{ml7}%n(?~H}9Yb*zMB@_EDLQcpcYAzgyx%uagF~dZt}u}Wh7sZ^rWi%LSm5kOpz{v!=hok27vChdOt^TYCIB8Cl;N@#6s`{7Wo`9maqeWmAe( zrGKs^jU@MG54}pn@=0j&_rjCTYpG!;c6MdLl_xFjEOypHrN6ru7x*-Aj82g^l^*Ji zT4(>9i}#hV`TcmfWZk%N*fn8qY?=9>8Tro$% zG(64|{j{rBvhVzs3XND{j>wZ!^-q;M%$snf>vzAY1}C9(Z0&_A!MuqXq9nyuc^_Ls zM){^jrt}^ZH|jAnqWaB3F$^C;qO_k873&KpCok>F3bj)t79w`wx0@T(i~7Q;dIL8- zZgmGNboNBIYMD6h6To4fXK*pb$hB-HnaE`0^h zFFB&24=+$=x}E%iM$UCG_E)D}H$2>uzeVmv=S}%5g{2#LK(u*76&M+57!)RNY3p1} z5F!5q(xn}>OV(gJITy_o>~%S6BF4Iz?SlGQ;UmHN6p>u~5@Zi9EOSG`Jr0)ozcKOR@ctgZ;i~YJXuGm zdKYR39vt5di#5~@6itX-omr3uk#U0_aw~@byvM+Q>6z$$=a%YWOy)FV%oYl(OFNd* zec{+GJym6hG%=XT#z<$odvKmxR(hBD`LB^j9A``?yIbWOa|x^UM)}{6EZq2lW0!#0)5`S0a8m|{k%!rKs$i?_^9Y9B?N%|iW zxWLQ-uU%##kAaI=Hhp_URn)X{DO-X9f{#vEhMZ!Hee^xFMGXP~0U4cXQU=%TGP7!J zJ6!xTHU1k-n3vFRH4W{|>pwGqi$&cxMCex3xycYO`d<-hJ)rN*0K92e(pcBRS-lmg z`cHQxUUM|wVxL$>lg_*B=k143Pf+x9z`&rkkpj(=dZnMQ^iB!@5H*>b|JM>oFY?bp znRkyzXwr!c%EPIuT?`4T-$813`EFNx!a@Z}KwNg!obP&a3!n>6t8sj&a&>-k!O#PF}pmu$ousvXhowD!l{a4jH6BVpsKMDqvzc*;e z=J+0_q6{a%n+1{d@1?-=q4YQL;Bvi9HCpTXR-86WYWpqauVAyYr&1*iRNP@|ZJU&M zR3An-h6N$`c-qLsdAP;dm|qx4P0r+rcYh+6ZnQFW=ck_cUaq$JPusX{W3I{v(5tH#x#wfFA^O@ZhnKPg*tMeJDhOoPJ$IBdj@N ztx_`=+s-so-@_f$r}&WnEG`;9q1aYY;L%^#_=Y{z#iz%s`g`rje#VRVTBf(63v}zD zibiZ^Umau!`h4nhRgC0O$EsBxj-xXIgLRJ951i6zdmOJpl9}ca0ay!2ICz|CIxJ?$ zjn8rP{?Xp?fc2#FHJAqHGvpG|w8AZ<4{n4U9$t)0pa`2%b9M64s`|KCRQO|kLcAxq zwQ2}0U|g=d=>~~kW&I!;te22bRF1qPYsXZ^~SCjD-$oT6-a=`J(ZFGQ3T3E6Q&LHHJ%iV)rR-Y z666)@XSaasQ}$XS;mwX-H2E2YCTBgz!tNzkZ-j`vhT0{24qiDw}-fMbloF@L#Fsup*2oGD&upiMwgIUV{UnTMFdo?N=mre*>FbNnMcKT6jD=Ks8V2fPJ%vXNKEtu!44E^a8NSTJ4d#o_xSMMI+@{4OFkjeQycic{^F>R_oH+IFhh$5^qrh{eNx7r<*f zdRi84eq|S4c`osIRBXIku|A=~a-|1j&E*BoyXSY6X+Rm$^R=f|aynXBREpqd4@~}V&sDDdxQd@5?Pa5Q#c{~F{W0yaUC0{xfM>DkZ4Kg14 zpuC&f(goa9mA|ktt33|s15k;_1``_V^{~djQd>E9=0EainTye)ONBX*zor#YHJTD)y!_+Y3sXt%V;Oy2ScKBAr)VV$uda+t&}wILrCCOy&?cj z+irh>w7(RQuCY>@t7vO*$nr3DBmYLUG3;7+6rqx|uUi34rLL|b zIfOt4SX5P4pTCyCvYV4xpdXPZcI!5^!)nq!HsjqNN$r69}_1l#^ zX1HEKLl(p?bgIa_5z>~Me?UYCZ$HTO8u=QvD{43nVR3+;xIV|$iHN>9UAAokr z$CdTWJOkU-SpO`Iznd^0Z6cU{c>kTfezk;#uT3N%yq$8GAhVRB9kGT+kKh8gG%{KN z^uYPX;NP7L_KZiL4_;YXo#YYZaZ0?KOa}Bm#@ip>^N9jk(^lbWmkH`97wH%j{s1UY zzDE+ezeU>Q7~&4dJ^yaT)hLX{l}O#7O|Vch=-L8JXozuPa#)x-hyy(i*R;5;3JL-k zq6xj!p4hs&YmNxrn^%+9%KslXOD=B4K_p5gc2$*Au7$%4w~LMF*2W!EpsX%%`9E*A z4l`|+%+`=J;3su#BVsyX=;0}Z)`?`;On*IN>q55tvF#n~@?kc?Uel*n zzM4#|Mpg1va`fY>z)c3F+w{W(o_RbwcV^cJn(He)^}Z75dv}U4Exy7FjdvqDy$d1* zrlt!i1%v`~5)!_!I`Iok2Apd=c!QapX-0Cymp#R2RtsZHgG;A?y4!tAWf#l=(qm2v zii=&Y8UnVqSbUPZidc4FCFA)ml%xItmWcxc{J4i)%ZXedgxDu61^yE7 zuIecN5$EFj%d!a`Kj}q@z%jzt;M6ME!bg%DYOoZ2I`5pAdQI-J6-G2lVI z2h|_v5s2f9w79d7pmV*pGAE`ELMojeC6Xzt5C8r8I);~?)13S6yV8~S@A!&eN20Cc zxBUj0t&m5Ea685Fmp-J$awdFMjL&`ffdL^V(fj zE;Ut$%3iJE`!upE4$2)Qwi| zdk$92m1>r|IY)-ymk#ibGrf7qe>ah+62{rT)Yx*(RS4P)v|?Im+ZKtfe%CZ0^8iU- zp)Xu!JEVtA^HtxI5;m`NkUX8Z^~pyGXLAqqxEjsP<{L#uKGV{gwG72bhm~|JKwqPb zz780e6zopHmx*d&!><@#_TQKsy_lk{RkcDCb(g}1z2m@`h5S9n_v(Y?D7=Dgv{A^4 z48xwB`AX4dCetUEpA+^Zc|#Va3zKaki#?HTirQ9urH^Q;z0~yS6=M3_V#X=mQZ6}A zwVpN)D{(xUU{|2=f;8}aHgEc@1q>8@Q@ldv0=K>0UB&7)1G;ukY+~oi3NXf7o(3^& zrYR=kbp3g&P8k!|N&$i>4JTAwe}L)}leM;q;~Oudn&WG}JT<|T7Y9gsL7##n?;Tht6noyjbn)DtXg%vxGa-1vgfXw#>?cX`J zhpF!80oY59$pqV`Ze;is!56KzLOfhmbD-2W>eAr$`gvo`SKjep5^2G{3if9c(At@V zu<2OcE6c#!FZ**iL8j+e)3JkPcKGQ_k3>DUGK%*JCjKXd)7>-U zV-9s&(Yck$cRz3U2yN=j6Lm+c%*B*VMH#bHPJ+ju_rj81GZvJh-QzK@Vw!)F5B@{t zqO4)D8C?;hg^E0N?MMxLZ!C$U6PrDYWRUdT>%ZT9OrlL!Ze{gTG^O+7D#gxaK3&xl zy|fh*jGH||l!BN%MaOJCRU|nD&f>ml_v$f`sR$L%Gwf37(=AMa7U(@|`M9O3`U@U= zWBc7&g}Fd%pMM$1d*j6=7kGkhGZM*$f_(W|Wo$m4-i;@Re`P>N1|vdb{ks1E?wUfo z!lxWO8&Jbv9K%UGuNWYeou|~{hR^F-E9uJxrUtP4K<{GD%VRQbb_@MZFxEPDeC-XY z1Kvi&-#4?&d*L z?FMU%5#velVEflh*Q{cd?@p^fC*J3^dL;*;Z0-|wm!DKs6S-f zvZ*$5DDRgRt7OF7H8iy0A%nB2mh9#bG;4>e?H<6&JA98sy!dRh^v^~SKhuD%_I0R6R=EQmcsFz3(rMvUEc`kkvw34N*<#xfK8Es>ow2*rmABUqOJ! zRj!d4JF4iSl$k@wnZ24iJcmi#gIT#h2xDu zBW5C!`Ygz-f-iH(&u{B%evjr{H|fDUT+OK5lQv01P8O%180{51067@txxI?mFFp&5ic56pX<*WaG+lUCE? zzNb1>P#xNObi!v>hF@~sw7q+UA)a-3{j=d46pXU?H38E$htO6FPegV zdJ_8I@Rx5;S9WAUoUazD-DJLNX!LOiSWL88wAQG-xf*|6ubSav zPl$nHKv-j6e)W5$TWvfRsq>+R zxsvJZBxd5%9xL+m`{~ABLTaAl8*>H-x(~Y%I1a28RfkgvMoA*5Vw^YfZyJXUY& z$Le?jZ&fwoqsBXJw=tP9rQ%Q6EGl9Ic_xH1Mq7{OJ#@qad13(jcRbieGxr%wT)tZ) zS&~+*m=r=<*sPbu;H2~W-#sF6{_YQ-ca?ein%L)lFcg`JwiahMKQyT3r{X%`2P&-y z3jR%^P==YlSvuafTP1I~gg@DCk$2HhCkcab+{ z$*iCDSE&}K=lJdJGh#at5s7FCs@J!uo_SXvh`dUx2}fQH%y_v^EsLCM|L|qf`RKs) zE+TJOuck=s9dQn24{xcqy+IM54LO3Bpw410e7^{)z)L!U#^9^*_veBDome;()%?P-JsIt7mZErFu{c2XmMuXX;vrFB7 zq0-z0IOnZot8q0mP24EdSJg?-HR_X@0_XOvI;%uRRW0>zk>U48dt?>iG*~q`%Nm65 z--p%zn4|LQAcoWz8(QXl4UZcdXCZS8Jt5Y@bxlRHsx;Q5s$2@JKo_DHgG)ciVV7M9 zeb5A(;*4_5l1I0-(){xp%yQPsF@NFmVCFNwu(JTkxo$Q&pZ>1OPBG4h-*F+zc4q3I zl@mD~i8H`?V%Q|^I-<;IeT~bp*SbqV_MJexDUOa%3U12M@T4wHre~&4h}xAjwl@~B_-;nLlgf2yhRjUDC)3<&in9q zSK-9e^>KBk&svT4L>X>92o%@7PuuTI+vISUG!EGm$6V6Jri5!x>q>rm9}0tY^!k_z zsH1jzJ_m(pLejeq6{S-)|KYCR+blTyAobz5*wdq1O}i+CPLnrp7w>Fe>#xOU<(Bu7 z%~aPgl-kK)mSFK)wI_2iMW-s=ozfP_FSOwI_ioGOu_9#n5jFKSSpaWB*t^W#(wa+w z?^?Z5jsw>uqk_+0ZTh0$&v;$)$=JTa%|%w9g&Y z5oMA&S!CT?xlrv%00-Lu;$ZfBJ&U=>xoIoZVN->s?^C8LEIF4 z+J=+7h#p4ke3}m)GxA+q$H0u~CKXKi=#f`#?X3DUq3Pnb&QyG1q9KJhzcfDq`8Oj0O(U-q>l(^Wb1zP7=yQM8EV7QDR(4FU zHi^&fM2?q64)i2ZXmEK0)Q~DY?26C!W*eSdc42#^h0654)zilt=L|-UPVF+`C3~H& zTJiC-27MWC!m@5som+QPDRh*TSCH5(%2vh}+#?(2KLDydJmPmFdYnmsXp6`Pcqj$L)FRmJ zan2tC&%WQv+N_W|BUYHpov0Ls(^b(jaX7gBR!wGDqeJGf1i{!a6rS(= zp+0O;FEqj9iKY8-6o}6#yFM??i*M@>ASH(ju?$9G@p{?iVSb=PnZE}N%F}&0xMj<Ul75crE|vbDYuQEC-OhjiN6q6UO!BXfU#M~W3wvG)b^N7ve|ho$QdYU*`n}lW zVd?QzH{fjYmLKv<3x1=&mi4<@ony!Jc6QV@suqjOz#%>Hzf&wTLv3EN%zU)Tv`JcB zFwE@qbWAS2r4#j4@=Ms*;Bl%4ek@)&m=oOCd<*&tA2j_SAeUfS56AnWhplksGE*2nuMOBl?-r^?6YmwJ`e!V1S zvGkyTXkE3Ynp*PS|M0Tvk?VJ%p+#wI|9Np(K88sof+R zAMNYWScjIKxADq(3P6)YjuTWiK(`0RQ|0M0a!k(=bD2P!&b}cxliH=%&Wf zad5(Whgi4u1G&{J<4sS{w?_W}><1AN4c~;MFsA}DXH9kExEnRk5;p@70WwXhp)P@~ z!29i|kw@1P;@?MCJ9iVnXMo7`&I9!e1TU;k&X$p{xmo#=Q~xRUy*+1jJ^EVZ3LD5p zm@A-yCsb|sN)IqF_B}vA@8E(66!P+{jj*&o61N`GnxDzzhleDyfVU{+U`(1sL8@eU z@k*(ET%tkI+=ac1Bpys%rOU5E&4`g3dmIZthyfXjtzlsMzRo--7(ih6^`tkY{zcAq zO=ZiQ<823o{f+2Zx!^fWdwnSkiWmj64gP5IvHGFF(!mz~xwrg;RhuJhr|S294sPv> zfeP-r9T6jfD_N!^@Sb=?-kc?{sc{ZcsYIh2dexflepH%whu7?MO9=+MoyhkK(=Pc< zSms;f(*#xsE1HA32S8-Dm3NGT0uj0pAn8?MCCn7W8u+DAmf?`ZiCI_?6$8@`~7ZPYohlL~z-hUAblB zi31t6*@vtgemgSF8*L|CjC_=1<*RxZ_9TBQPwHoa?e}!tik&7@^swzVbJ9lD-e{o| z$=bUceZmM{E0bj|I@G4^3mb$EkV`jnIjXpY^a-Z>(eDcENR6GH z6mW1t!k%!EAclFOrT5R=HRUxcp5L*0jCSw;-0l&-QMA*$?FEaelgT-LtFboFtXV_Y zn>sGxJjz>I=J*Br@>Y`wu7Mw-Phj%8qcT!sK~&EZmT}0 zTeY*6NE{6eL^CqvDM|z4hK{{A`Y1YgAGNVt{5Ee|spM+J(x7h(H;Px6Qt%n+q9G(e zexSp{XSTc5t5=`oGM@~31v}EylM`w0rH=G@FS~OQLu7imUZS*>M*5K^;yHOpm(T75 ziqzIeZeZzShPd;vmB8VJ)8c|mX#TwuBY%st=)Zs1z)0BvW!q#&S{8VxK7b5B%!-A|5lU=H;va61RV%i;jjNlwEay{iCPE&mDx7d-u7?1hoGExa`q~ zqN?Z^%RP5XF2~3dv?P{gLB=$fK#% zveXu4_}nRY9Od5_|GH$`Hj)o@U+Sk#oTYp?iF;~*bRzbn=d|hbvai>U9f`G8Fm?ce z+jLS%(u3*3uC!&WQOH4?Dvb@QTkb46eqo;VqQuM#*JDXJqXoJA74#4nQ%-A7ho&`S zYNFOvD(^49KX96F4)#KP*5}G#y87np-h!Nh{Ra!EiXyVhRU$eePpt}ANXp7<+a?^d zx}@Eg89&A@aPek@dzUX)NnDvn-^^_99OM2C6!IfZe9mRlGZze5?*A?^%;R{Q(OCGr z+L^BfOkc8{5q=9E->Ie7tB%Y@wKd1UpG&?xiNl&I9tv<%N)$a74a_1P^{t1+qXHDZ zonP-MuD)c%H_q#Ugd}o3R>^hWrA6Y;b&=OLm`f(F3|tp0eB5Tx>W~QOue0;Fj2sHU zz*Z?r#c2XQH1*gCu)54}_~caJaMe@Krjv~lMkNQjl+fc!baZ5=$9mQm7iiwc%7oVtg!TU7rhBlE;@Z=4Bb5um9`Jhw3_?g_^eM6 zq3~8lTVOC%&E3IiuQ#>29$vbTJwC?5_?Ttrm)k(cmKjD;RkS!Z1Fa6ze+7&yyRu#__K)WF=e1?}RqQ*u52R>1SO&vR^FN??l+`!1go+EBKhB}v zq|6fnwc(+8cmEH7E+DR%E;=^eMf&Fh#|qqwqBxBll{yPOU1`T!RkTAyt1kWJxQ?wDK8lz%_(kNyzX)1hn@7WU>nSOrN`! z6-Lul{1+{{-&Od3eE+awjS0?&SZVWy`RJb1AR z^!6b>7-HEineE1mg-+YToU;rD|kP1K#=3LH4#x!!GUM zxtbJ~YGqDe2h9b{xzx(v;;ER7B4mwaHx5TA@>mhI%CG44u+{A|(ieXIASLQ2A=g^! zdZT$txEq3U4q39jW+;7+n?mW8G?1TyH)*O;QLJ$D?Hl0Tym5O_!f&^2DpQC0M<&Mh zJ0MyEDw0~gc&9?Y7JlM_osVPF^J+*d@m$fwSX==G#2Z~*PJI+z&JcG@L+N{tSuT}# zy5(Lvum1}rjmooDFHWmA!)wRRFN$4Cx11xMe}b+{BJ(egUn&)c{v#UqUjXW_kAJCA z^XCm)zA+c=Ie;>Bd{K&;mVR1)$0+8Dvs&GG8Z2jbr#1aN59VF$dZ?M3QF>?S2jARg z{*+vTPj#F8=j@l2S@%-)5-nNsI^NH1-{3*@BbF=sB2=nBt5i``^(JFBsT9U!ZrZ=HADu&r>7zPk>)Fib4D%83xw?M9!*uDuVWwPozW zy#wQ%mU;5r*h%#*@g%+S8ju<}Q7e~Y#a_zHyq9pC)z3zX{rlSCN{=upwV;R8isZFGX5@T+a)@X zQ2PiLn9Jo3h`!4GB=cgoy~JFjXQTS*mg-JxR9WP(8;@fJ)We8k75OF|$Lf|%_8G$S z4Z#b40C3EqoPn=R%N^x3cC0asDA1XTewo#EZohO}A`ITtCzOnu5xvJ*3AMOsUAYu4 zbM;2-)KdaJ9dnam4?KIw0$bhK(+eAViHW!xctYxJ0~blb5xkl1C-T;x=CsUjtQnVXb2dpnm5s<2MLMWu9!Q=&{`w>H||k_6^zL5{TI&5xtlZi_KT zQsdVI5BLPFf$sY}C%k)e6F1ABO`ZZxgiRyIH9k{i(Q#Vf5S}%kbBTY1L3YDbPBLgb z02x{O@)xD;MckTg@uPdc$t6w%&QVw2_(Yv(vt;N@m8wwyYhAAr4Vmf; zlrH}DD4_z|TfK|D=(>=@aGwY#YHD-!Thyb{2d^1zDDSzW2?+LVkKA=|5k z6Br^coW060-N-j1JTLWKR zSqG*E8{VoX`A2k5Fc1%%acnxju55%rm|kzYJHA+o1-Wr2yE1llG?Uj)Y2XmYJ4+vC zbkIJ%HI|Q@d~47lLg7`3N5A-+@sYFfPk}meyghmL$~K@87mcE;TyNC+e##%f z$dmK6vEQfb-gCQ`e^_I&$8CJ8k~L1u#n2gelWsiFDF;n;Guo>T(Nb}O(|?FSQmWIE+R2z)kqslHdqj}YnwtMtyKh$?O=O10uU#5owdJ#bPD>t%sU?^M!7N6UAU zcVGB7hOu;K;U4BpORehAZ(kU}OIAn*5p*udHcJMp4S0z%F+@s4B>awBoZ^r>D4I_7 zvl1DJ!py(8?Y8D_4w-YTr*G9*4E8s9Vt?5KrJ0b7eP*1$&R&6*a&FxdN2c+e3_iHC zhy|DXy{gWo{{xsZz3)(a)Z+p9Y^J_DM^;zbz7L&JB~O_=CQs9`<8YeZNcW)o1JDyX zzpQKgYht77@*Ti)XD4s)RJHl#h}FxVXp?=WxoexD-FG`2WWLrc0)Fm2UCc53!@$U?W|y zrJKta+4;?>7{$LZ7p=K!X0Y`EoIFJ>EzxAbOzM$sh5TXw+Rtz6(%-MTj0#;r@^7C2brmWOv@aC>a$wV+?z zJW<%pWR}s|sYo1DIRw!P(^u8(n{wK{7^5@yr$yoT$WN|P+^`b_qlx^cMc;e4FSTzJj&)2p+r&WL7N?BrJ6iF%hj zW&hc=2EQ7);qOCjuJUh$=fIlPjjp@+p|5?||3%$j z=t)3I=v8SVC4tbS2m(??T1Wz@AVTN_sY&S4J189j(m|?HL=n(8_w(*@yPoy#b*%O7 z4|{*fF(2lb$y{^I%ypjU|NlGtjJK}}oG8lwklWrH{pUMVv{-K~?)P5fHKGABCM2tD zw54YbFbR_-bDijP=-g%oFls=m}co)%(@&f@Kzk9 zeDwd-5gU`tyBK937vJwX$>)YHjX3xWba-5ENUe|}H@ zC)XKss<>NmP zOSMsp-$65$M4iXtWfHD}R{d>~g8#d)iUohJ6Khv$w8f zGu4pwq-kFNpu4ZXYMjd(OaLV&d7Ct;7gx9s-jlpiKl$#k#KuV@j>nk8UeEgN{9YCW znv4_2;VPftghYVQOG%I)LHa%euo2E~4#Cj-@26I*l5KjvgZ; z5rS>jPo2o*Y*h9~^1F9#vkx9;5Kr`3##azLQ#*7yHh0&+7XoXFi0mCJQxL;X#6P|j z?*I?)KHrVSV#LC2*?f~D8f%ekj5Hy{N%LzhJ~L$L?!9^$1=S$=^x5}nejS5>YAam3 zNXdssEQC$`j!c#;BJd7%Hzmd64K5$?ay{2`wJU=gI=^mZ3FRgk0$?U76&uq0FDW_m zzq@F2YI3hM%uh;mU8x~Iu;C+ID_E|$tUe7315gR)N%QaD3ie~#bvP}Qb=du~yIX(m z>y{>T`!YhGus4jO^Lnd=LNx_b`O_Zh*dh+TT&;3RcdnNyS0MuT$uf75sxs=~02r6@ z7A2ey>XOw={6`9ayZb0;G{nM@u4ilfN9^m7*16l_d8D*AuK9M>64L7axA+_GemQ;i zi77v-%+j3( zdrR&d3BnU_rxULj#2DXiFDH<=7{dXB!3ZOyweMX&zWABTm~n>9?wQ#$?g@AJ9XepX z`gSlFKj+`_r)IE9lEuCJIeh$$3~0gg0mr$dv^N>!`sMjnl;4x_c@LzD0AFKxy6kI= z|Hui}4%WInKD}Ml)o;~?~&4D6@)Bh!_CwX)vZ_6NBW>ckID1TXj?>_K13_gQ{Z z{if3Bmr4Hoo#KaZ#}=--v^fZne#}Tb^{Gke5V^Cs(NdNdm~^Fb35!!0U8l?E6%KDF ztGW2Gw9w)J&cYE&u`}>S@MqJ*?No=y%*#j5(>zlHu~!mF$pSXCS99S|VU}JrrckhH zBXwfJ{M#Q^wAB#IrLnb}zp3&|el1UQN*((1kBr_P$g96N&4?|0S62|b&C9HI7bp2F z)@So;gu>p0?u)kTuYl##)@oKWB4T)`5iDt6TAbo(5BlR)Vp7N#g4@*>d$#>0855+V?|D=O7z|` zr}>{!nTfy43<(!@2b|VtbAV;H5@`m>7c`6ov74u?+bn4xeId+Z#|;lcLqqb>3p{w0 z_xSXvOhi2o6L)>t>36d{ZoJD;k>>U?|9J_Me<=EGVv2nx%QP^SuVfk>6j!e%qZDV~ z_k&O6pVV&U$3OGDD*c&Vw!D(K#R2#9J?wJG`D&&h@(AXih&kV1q)XpzQ7F`Uf1@Mz zH`V4-t5bisKOh;*V^pQC4UVQ4W&vGt${)*cPu*sv9TQaaXP+5o#b^nK&(zY{B7 zJ11BlCbgxw9!>sb>h2M&sFJ{4LRKJ@c77BAv)TCtcb~uhn<}G|>KONRi*dENw`8%r z-vY0dIBn&OK1){DZUN9dpr_Q4{s?K-DOCuH($IT(&UxGoULgvx>=H)CahP<|#D3!< zH2$XQHckC;W9yl*M@A>?Kw(g2G|~$AVYGlFMu=Kgb*UH%?SbB9WT3tkTsU_6uO5?e zC)T;?#Kif#1igSAZt{q^pwvot;@Q{IYs+7sW@mtCcn=~kkox~_`I5YA`c)G!Bvs)T-W28P` z?7=qQBktCv?GU6FD$1FVPLK| zv%V*~o~H3P5fSI7g>}MW?7`MWAuKo0bjj7i*;$wB7oMsV|B-m`0X}W>Za8dX?Rdar zv333Sjc<@wfKiH#594cJeM7s;>4%LSvTdd8nTz-5E#dkGC8PvWkQ@(yu@jSbiwJ|i zsnX6LjXVFh3Y)ofXT`Zi9>VPY(MRM*j_ot5^rBVXfe0O49q)Og2lR|HKR}W_gLQby zBlnNR_A-~-`^akdP_>y8{cp*>q}kl0+QD`mjp|p}8Qbz>>o2wBnE0K{!O}4|DjlOySmagqJQd>$-ZAU>3fJDfGw0`f+e}z>4mvrhPkeOwu ze)%$`0^?ll7jT-#lqgPC(Qb9J0b5>nInP9j<(JC*fPXI_Vfvq)%Or~#d;evx{7SA~ z<;3XW-Qn@*W%8NyTm4~Kbh~;~B~DQR5*p40Z+JPhT>dkvyJQ}p^U(`Chl?Xu2rzxM zAWLf)JY2YV4WLDz=wQGB0jLsYB_oSKQse(OYW`o2%ijq06ifp4etO5M2wtwwS$vr{ z4^wk_l7p$t8*Lb{D!wa@O6`LW2S8IR6aLc+aO=C;LoqQb1ag%5wi>v;3=kC9A-T!=U{!?=QKR>J)Px3&jSEJjx9{&tN%4HVSs{1|-2$VuKlhrIiJyM0EZGdtpL(MjA>0u?zh>*IfE{0<`LMOcrvZ+Ke2Zs zX`{UN$LRi82K8EB)-+G2*zRrAlLa4zRU;G@A;Jk=7R&DOuA-Hy@d~55_eK=EUnkf% zUKVG=l%ueHPpqY)8^mfiPl}f=5bf|8wf95JWlkg3Jg#=>l@GVO7%nQm?fj5AT&==D z|CFSEvN8Rq5oq-kW$ejgE86TtkXb&y*E+`Hd9AR6PR5)g^oNxxN#JDmVI$t+P{XTL z*hYA5;`4;d_V(?VEcRKKJY##ii(C^oI0C~f$c4`a8xNR4Phq}Ued)Rp>9Bl}_S?gX z^XvS`dqv>6H;e#r&nEJAEjg(E57E@b1@EJjh3Sv(G05&RURAUG>wQw1A3knd=F6Wi z&#ulB-$iXG+=kt0D|Y6Mp#JfFy8jH#5hpiUbF$Kb8^ePqmq=koywad&6oV z%_;5ZyE|rxpkntXmnAnBwdpz!Gwwv*^|P(%CTR{%Z8)35(e@okWgST zmg?%_8T^*~@$^#Mfx#y-_sz@kLu^by4hLx~!0H_XWr zfU?cK<zUh7t&R!n^~-p_Q=qTEXW7HM_PyWFsnBJ zlk9$}S=;Kzq%2x+Z(b>1H8<_6Z1sQYZUgK4m?whcF3Jm=v{+;m<^qUwC~x|?Z|IVJ|GeAC%}_AQPZ3`mpD0< zGn34^iGrBAI$-z8*wM?XB-bWudFz|4JdXbLQ@AosKzW|8tO7X>?xq`O#fA~jy9c1< z2rOi^Hl6)f3+6SeidQz-25}J!EH&~s>j@_BEf(F#0&Kk3e?3XgJYher@t{)O1krWj zReaoPmJCe}j-!T+69-=(ziPhz>ji#YHnkRn4qp`jb6=mLL&!gLH0$1bt<@*oY`vg^ z%%&+%Ph4DS+Ju78h0$T+ep_K`(;2G49=D#f;18o6$mA-r?Bn42dD^FKKx#BS?+TNj zmP_nD5Gmvsv^n;_4md~H+VZj2_&4NqNAvMp5?yb5Lg*sWMDWyV0I(_qhN2eV;(B-c z+XzL~^M65tOH|DF4^t-vo7}W3Z+y}Ssz%;#voTwvEF$&JGy92XN(C~MjYX23EX(^O za;#dH)rdMTHo}klG+ZHfR%JZ@bZ-1Z`PDITacA~m@8j9fdrwRNJe?f-EI@3)2P^*E z1OjSNXB9(p^li*9aSUEDwK%o)HRIVCF}1?oDFY^c@CFq{{Q7q|bM3nIKu@D!B*KbVt>FC3qFDjmx{{fKxO_JMVwHWBw0JlCC% z5`@rqf8IBFXMgpZL;KTj&G^vzn}59O;x@S$^<1lz$BHn@zO8$FkFmpvP=6JPhKzfj zr(<$v2EVBuQ4BZ(@hPrnpQ2g=Ie=^IFFP7X_nD!KMO`)!L3>R?W(iwe)R?2mbmNQsaCi?jE&sH#B&I?UQvR3o5xltVoSf z02Ss>Q)t%4|NiZ?JRhG&$@bxQ{~(d*G})^J+_jGuR;X#VNC^QgfTaCUq}kw?crqfx zFJ2N|d)7ID!zn6bg@?3nDW75)uOye+jrLl9ZN2@VpI}(LDgC>l= ztuQ#GAbGwtx-6PWERAFFnbWnl7bqb_%S*mRR|YT`3?w8RQ4d(Z86Z5z?-^#o zx~2Gqa9rQc-?Sb6^V63;`dUslF!4^3P1$#_S{`A_c~}@Ok3}*pGEdAYnYGQ!eP<|&d_`%dUg)N%Z9M{qc1ho_oFz1pCX)sH4(eGRkcm^WtlZ~^ z=&Dm1DWtsh^x>`3U?%C?tS732@6E&;m8;(AWsHzGcg#vSz4YkW+3)!#UBht!m|!2? zH86@7IM#0B35+jlVjYdH{~+qV_+e!j3zL#Wxl8KF^!xxQH=o}joA`FABAfKkM<-vCx8!Nqn?96|=|AC2_#7kBDrl54UdpF0_crR=lZHQ8#=`dRCBq^TClxY8IQl!cVtx^|GVJh6t@xbKhOl7^QWflgtuLv#}A&?nguslAQ$ z|J&?l*UAeo3m(@mrLlhd+8G?@ZL`CAv|N(E7$}j^r)1QgJqtgMx8foUcuE{A& z+ijWDgP)K8jm;7Jikk+4)ptC<+&_mIo+%hy`jQ1z*v{|-&y%aq&#nhN-uwEce9pf1 z-{6`54LAO`&2RZJoF+?ayn3%1pIoV94_YXgAyvq8^PVkU_ zOYzq5*c}b8q36q|w=@5Z&FV=?UbG8-o_svfCgDp2?2{-Xlu-ZBdRR$eo926vo#06W zXiISt5&WM`EZ4W~V@gLIBbO;R+uW0E-?ANgw49h!GmCDzE&E}ugA9U~lb@A->GPQd z8XqX0`%z{vvi|>v&CNM|yy6cIZ}%*ni&H$d8#@_=NoO3E3B#PMSG*{n>>Lw4Kbf7a zGcod}e|BC z24115Gp(R^te9?GDIYC-r}adtdXH1tX{X~*2i3FSwRqP= z_-RIl9vV|vPa(pm%ijVtkQOf<{s`V*wO(=>&)^g;6xBV~tC4skx1$mYgDieE?R(d9 z+il3<*eP^yW#LsdFK__qi$eVj68FkPbKwtZ1zYF~DG-^wpp1kycAPZ#-t;oYMxu&G zID|38qk0jAISAhMfMv-#7EUx}_E!^mFW1W#?pD4E+DYoOo~3k+ntruF+qI+~3uQa| zi=CflXwaX6549dg&rbE0?vTeiHF3J?FBjCGo8)BsGt#n|1fIY3aG_L3NBo4*IweoNW5(1fml5>5_=K|0O5BdpGKn)xPG0RVqdN z7$NDeg$n6}B(7>(KJTo9)O7chmJY?Wk+0jl9=3If6M-zLnXS)?vy)a?zirZpb(QhN z@y7CQ+SGeER5xFbgX6K!Z~uU`EEt?m^Xcomk9~q_3qTiQ$<>x}unz&1nBP>8kE!Ke z;`>sP$>!#378A#-GQ;4C*>{FkzS;h5m+xy6hTqkh!PGo1iHLhH7YTKqsmd9z#;;|5 zUbA@#&GgfkiQ;8sq?mZbx2Ob$P<76K6HA3iBlLq0Dn zP23h+QIGyjwPfo~-I?QXq5i_jnX=Ir2T|9!Ds!S|0ud4t&{cnC$^-%Vz5Lk)Ot0q0 z?`TEbtSH)BiIMG%ZyN^wVkhbnui+DqFDc)K>F+eZseZ#<=`@{E`ohk!Y2Up$+Jqmw z33{tXr4#mcZL@gWvHCk_g{Nr6?!oL(M;dREXXD3;wvOLaKfYK@PcDJm4%TIX0OO3P zC}2cjk*x}A+x{fyNjN^bf|B$M9R~g6yV3Quw8?&J_T6u)0oiI)!Bzd6mEYO1Dbtqi z*Hy-+6T3}kvv_RW)6Bsbxk`iYCQqX4nv`7=94p$Z+i7xNOhU{V=sNt0oX&3pGz)(@ zO(<-o23d*Dr*z-y!eo%G>b)iCvg8cYR5ST}6`M|o3-o18#}Zkuc=}>>tDeoGoe2kb z7v|!imM|Lw*=lTjxJF4#d}#4uSMm_W2RP*u#d6TbN1!dS?`XA4ZR&c5fXp2+O+469 znQ&Hk$~t=9+w@Vm=qfsiT9JmNwoeC0X57mRzh#;+xqrPdUAJPEc_|I^$->iq^sTJ! z>yx&JW4}P=O$m8`AwLrd5wRw-!T0d~mKXOP-FPnUUky_7112gf>i*1>e8V$*34E!~ ze(dDrTS3nfk+WQ)ERqZ3z_!foSt|=9%)p5IyBuzFl|}$FI=D1gLK*lPmoge0nDU$I z^k#gK&#LIJ!>svA0T-#>nQijxgYe%}k*&J$@>#5aKg+_vrz)cpB2|w8FP`j#G@3Yc;9lf>{;s4V@u0unH0V5ZIAvMaB^&E z6~fx>)rgLTF`d?m#~WIuz)L6DPhKTlb?~ympC{_qA{ztVNV^j;-Iqp`uM%uDGNn|1 z#4zZ06^gSKA8p-XFmN1K#kQ)`#)p9o3LO3 z^?4gu-mJv7N@xYi7yC(3Iz(s{m`uMHb8M|n#=W^t%RbJc>-f?$c8M6L7<>a1390;a zzR$}xuY~Y=Z9xU2qEV1>JpXM>^GAwHOUhXfe_RJ@_+>~Mj4OuiZ6zu@kWEeaIi}7v zOVGDQn_wfXpUNYa0OFRlb>`MFhBB$aLnJvAghPf^qEFQuzoehN9*luwd>^if96ta@ zuKBIx=@p6>yx^+4gA-0R|DJf#^RBZzo{lJzh9xi92sYJCEg{WaB@1X`^YZAF3=NqQ zd&HtZduPh8IOw|_%7-*N>{u>ap=#rkq&Svc<2`CVzSROOK@FE)6;h#egE?USwWZdL z1iDYtr-XHuxN@9i`0Ay*Ry=pEJsqt1P@yZJiw=%W64fIW!5l#fj>G6KwkfS9jN81p zw*-hoS;%zotrg2$CF`*t$>QqA`R;1XDAvWuTf8oOHc=`x_#$Zr5%AvE!tIIjvWx`B z!zbR%$y*;%xeM%!zsSVtBG*j6zec)A=vuDSp((GQiKu$6a`KK`L2G2b5Jsd&KN9Sv zIC8~<2%F<#<(JPe;~4W6`*&m z{0*t7w4PYsSmGd9N=`k{B{x$@;LG`f7-SfNybt@7hj3*6n(qSJ?79C z-(uXX5T}i0s@Da2QBip!(sroMb|2@KvTRanm<1L!KL^eZ`MZ6&%tcRA@o~1cNW}I2 zMUy;<4|(XqGpL*Ze|~sPP|N%4v5LxRVC2p@gV)$pS*UE~vq~Dsf@gy5<&DX}ljU!z zme;P$JmgW97etJb$7=_bmI)F-Z;(Aga6B)K37unbtLJ8lTNlu5ilteQW7ZGHb)(4G z4rN9E?pdA{>9If}8KEGk-q0y&uh#?1h+cN?_iY5QUlTQk_(Xs7nrPGvyWY0d!>%f$>XBu$xi5vQw4dN-1w)Tk-xHnZ_8H$@_uf7 z%O6U%UY)t%p64%PaAt#@B%pQcK>%d!H`TTM@$K(ciET+ui`K5y7`z+kzqo1L3i#Uq zim}aF=IO;Iwf~7td)uDqW}tb0BDCwrRc)XfYq+?l9!K61FriYKdDDiD^gf&vO~=+-}-)p zz-4;Qcgo*kcT|2>+5*4VYw9)@8+*^|xn9)?shTH6v9!i4YdxgAKV*cKRJaM8sR((o zw&x${y54vxOe=y|M!_^0G$-!DY})fXEC7qI1ud@&MwpxeFU#k}w8`d$a6W^7z$Zn# zCCvyu*%U?NWz{egi_Wnx^T4z!z+9Ojo{`8n``s)67WvC*4?pHEcP?V+foq!VNXn%G z7^+vvq=p1T<$1joF0vNq7K{;Q1X0@^>gLsbe-!UuJ@C+j5l-*e?QBCK1{OZEMCS$A zT-GMG`u`U-Dz*JDXw?3A*+=Y>qK9ciG2I?BR#A!MQuHDE2{(3d!|Z*)gMv3$;)4i9 zdN0@idP!QnAWD`y_~-g{@1$&(qXKk?I5^J#FdIu_+XVa7(22M*L1WNDx%W--0tPTb zT#C+pepkwSnP#l&P$}d09xNWQq)Aof>?X}vlHMcT)h{_^ao`H`G~T2~$o$U}{B;(Y zD7TMa%$1{x6IjoB)G`4lFsygUD;>|L(YY%u^&hXCU$ga;{!JCH2NG2?_Z3=ra$ejw zK!4{>9Ue7xTm5sH%cSn9Ma894&&~T~N#QF-mB;cEU~GM3Eg`RwtH%4e1riFkkx&C= z4jWL|1h-;EEBt&e`>>1%BuUvbB9GyX6O+*Bb=kyZI9Q^avY395b17yn7FY(-=C}!rSF!|z7Vt$_Gx=9$5RTmlK8b(=b2;+Q{uv< zIkEY`sw}?&t46_jd&|R9Ivgvv@VLARJ zUy{nrUXnh=Ev9JxUIX-q;I`K|;Am#hXwWvi`nsU*&;yMB5w$UIpSyIghUad7Dc ze}{6P&5^7X4vIoOBi;TCx3JpoebqBtlxI55JMhuNQWZK-uDW20u4Xa&N_B!W;}ws zY9RS74LwLTLR!w!3sCjM9`P{caT&|@$J-M%9aS?bd5*M?QAhybFCZ}oR6Ux!@r|hX z3jQ9hXyfqd+qmF#NG)R~-ytfp@*p~FVbKbTt}cmu>V!mq&s5d8?Yfwu>?GGNaPRmf zlyiF#kXqI)Xtt*GuiXm(DC&Gf6R>NRoay1SwKw<75nJD1a;xRoByFk7_+y=Ly8YN& z23ZMRCNaMx68Gq1GWIu>xyMi7%TxDi+pP&Dgk6xPHJ66_c8nqQXE1c=%qZ*n;#YtA zSX@Gh+(eqZq?`{&CLDvP^N;bh&@EQ~aNt!w9Sh^~hne6aQi-4VxqYVZpY7fS`IQ(V z4qvEjJEypkAZa1dglUsFN+ZrdcO9A%_2gG9^#%yB@pUz$-nF3TAb6pp(dyv z_?q=~fe-|HJAOEWmi8?OBNNJ%Sk>=I$V|H?N33wyxFpym8J$#x9(7Mg33V`)`#2aeXhC&tY2wk~dHW)U#T`ZyI2DB0pyc%&dRm)ktb) z*Odw8fhOk7#I4U8u0WLy);?WTce+7Nk|gRet$1OC49Z^digQytw~hA`w@}>#Y3H7! zkD_m>LzoaS0ah;3l8X}>w{@Nke_E`&yHxafmQC;4#)n8= zj=U)Lx~>_f4bXPq(-J}m&p>F+wCth#YsIU#M1jvmY7W~9sSd;?T*BSxeRl&uH0?0$ zu?+XVjMNO}IDJYhHCvEI-4%&3kw-yF7H~}f69^y;Q-hGduspBbE>TJ1(L${!PhUyn zuZf;qDh%>;lI>`5O){x<;Kp+`$omtv` z44UBem#2|Z z8UpDhw_aO60v*PE2YK(QloAvATYY&+)4{{t#TdqmqwJ;~7@b zB9zf_V7$(}9cO^7o)`Q1qL%ZEqORQ?kCeyW!}AY)x{)l^BQqi)RIsGTc_<1@@ljKm zoR>}-zmiB!^;dSOg7 zD4ljeY?uIPrUfxC|H%kHz!~h;)cB6vo7*_9b^6)EYM1YwYgRJBii91Vr#y8wf)153@ONla0 zeH@~sU92S_P>v>wvSxQ*Ic49R+!rXms`&-m;12fuzB$vmK+c|fcNH8M0s+7Np6nyp zO|=T4@+ls3fkb~xp7Y?9QMH~yYsH>KDU{J;8C{7e497GjEnbTBK&c7yU;IBq9GP_O zy7}qnFE$wvr)NBU$MezAwoEc`BR(}`)o2s9_&x;0wLe`^ZT981Y-VTEhD?IhWy>kk zgo#yE^jo*YjhD#zT;eacv-DZZj?9)}ab?;8D6C3&lX)o&nSK) zysP;~THl@BKL0_B2iKzfDNCtM{T|~&`jF?zTuiCs^tSHR^}vp)q!Q*}(He!cN=xr) z8ElxKvk|4=7H93(WHTl8==m&-Ort~4E3NMpoBUOzEm!>I-ckNYukxHJAT7Rtb3CGHIIctI;N^1*dC{ zRZovD;eS&FGRjce1)&27;w*zDdQ~q%CnlR8yeI8KMtzEqBN<`pQwg~E=(Zns|B)< z@;uAcv3HqN<2G)7_v*QRSP#HQ^{vzZF1$|H(uOK<@XE<0_kgx|aJb~u{_WnHpY~}B z&%6C1Y~pA0eYtGw{l>|G=MPd+m{hpRBT;5q!p7ZfrhR4Q6?)h+sorBJ#R~{!!H&ey z&aG&58}B(jm!CXP4rrZ=35*El$9XSvXnAa##s2tIzh*4glPXNEG-*k5&62&UZSS04 zCuw2zp&F-@+6u!6_8Q3MT5W?n9>{VRNVfxxad>e)UY#p^eFtUL5AvML-@&;dMU` z%s^#a5ab(c!o_!lwof|v`g74PDHZh}miN_rUm_!_D%FS8Dea=rCCE|t`|j3K zqj2zOpz(Cf3rI}|K?8gI`rH7r`btq}52e(q%K(jD(20+E3gmUF^-VKhd#Xguz#Z z(v%(E#PN915EMS{p1qlOI z&1N8pFE#K=Y5Z2|$<>*y;j0X{jcmSDK6C=tai-N=*1b|IJAU%0+p1M}u&Bn15Ei}8&qGL>sO|${B)%|jY_P{ zuPrN-jVHyuaopZ{A8dB?|NMeFMPiMR@a8bBY=Ns(4Cl{FZ2=!2v zAK_$4&8!=1MPPz~VP@14;w(q(hDYZf1M3v=yvMnn%zVC6iI$?hxA-d!`B>h%!mVZE z2G3+?5*HXjrk6c0-s+Pco>fb@IrXCvUu1}Y0f#4kQ!(6kmeOdiR%2k+VU9tzGDt78 zAmpEz3-sgDNTl2$8&B_wp2I`goSg%#7(FT>!^rr#{tjCzfgFj;vQF6d?UJm>g_J4; zrvh&V60IFuQyfBsp?X+D^npFb=eT=Tn(Uxp2=aM1GuF8sv>!c(KoD6cJ%~~dt@{{V zuNEP(*xn)GJPBrDLE-|)&n%6x57}8S;wR|+<(fTFbn20fxj0zE6-B<+Yn-S6`I?*@ zn1zo%wZ0gg1P=ynIri<&m|4re11@tS7$Qc)(u6#0&_B=4K2VC>c{Zv)Hg@Uwb}tKSE1`FNBovbjA1s>(zrf12Lw=aOk-{5xO+?C?`g(BNRkkg#@oVmnTDhr z%S3f_G=T?L2~=J-JSow9;l0iyk91CxWp@(x8gP26k-M&BVlAR2J;O3cHVb(^120FD zlYIAt{5-pZV-Siwk3r`}btv6lz)+YX$Sg$07FFe)>WH~GEXs!Q(GA~ zZsLZ|`ESOq&%WWUZ)DOS6n4jwKSD{i6i>JjlK@)(q&Bh4&+>nj5&b_>r6%HVlP>Y) z%H4cnCz<*-dk%3eo|qqG_)|Zf7P6yt$yt{NjdqtRbwV}Wt61n?hpWh>y!qTi#zAes zOWG%-6cT+3XpBY1DA9>&YkFl98B2@bAL-$=uQm#TQ9Z8hO2d8bP*h+97fL1=dLnZd zW@wU}7xKYC{-orv;A(AVaUW)?qPE|LQ+jVBBYMxFu**1{$r9pIC=E|LgQ;V_tEmp{y_+6i;Uz!3 z>5v-SaN8p4f@Fdk9aI9!2!*8#P!5L zreG)`Eg&s0oNGP(*DBB|YMCJJ{JGuK+v&bdq)QI_+>0*e`E?ZnIhK8*Sk~&W%%J+@Rfuq3(1N;6PSF8owb0LH1VyIw&FJucVFIn-#*0 zB!6nM?b3|kGd>B=VsU5H{10bWpfp8dy#~^=oq-%;H$@tM2;Q}NETpc;z5ZV#Z==%V zmE^fqH()MO>(1OIOg+6T_~dkHoF7UtBV8oVc_?}e*_QX0%`rJw4m__=Bs_#V8_m2X zen3m06|#9-vp90{#KphNf0;1%@-!#z*JBb+>@yZ9EQ1leOCPjMwPFh zovB|QSRDVnV|D3e-KoM|ads-GQgy)~Tun^SF^BE7OtP^ig~Nc`rHnaMXCybaCs-wI6;J&PalxP3TzA9`lznayo6kM%lmn+I^dp<%96kc8&7?k2w zIfdZUHhQ6q=oFipCRHidf{_yf$R~j3R?FGzvzPXEngoX|hr0Kf#{mhFfLz{Q6kRcZtW zwX}r3!WnJyj{XuN=!xUi-b&_YM&9Ha2=TJ5YGzjdENoN%?fWSKW>#S@^iGyb zhJd7H!Ow(r=j0qc$FnYPwGOp%KIfNFdUfQHFt-zL6yWXx?k=B+GooS@3M9qyhvJwvm0W$=9jLw)CCI3-k00YIPBBpp4(2wEAmI}_Nn?dURS?n$Xdjj z2BijAU6c@uL{NIc5oZuv-vy!qRZH$oPjHhcg1AE$Z9g~NKG-Nm0ftLFlf@RgM`WT; zGc`27_(54&)XTcd8)zt`p)5EQkpc*qd##96sy!UB~N*Swwx5+ogwV;%Jo7W-F1 zyZcL?Y?ko!4Ctj>rqZqS==q7|Lb)(S(vnYEZrx32DHeEjwQkB!hO_ft7C&JoYU)b< zqnBkFc~yL#9#)?ZC%?j*FQpad%WV$2O*rv@T%UPXHd9TNo$qobKTF_cu<hmW&_Ctm(43+%f{iNBy_dA ziqv#a{*%9YdlB-UmUOxN;X)>4Dt7`O*47DZji)z^o|%`T!iBdh5kL{i3esu4H@$6u z*DBm@#qh^_m(N5xO8Ce^pD&9yluufDi$ro=4U-9#)kM?uk}o#A@2()je~mAAvvkUM=N2Z!!P}!Gk;@<<#~#?!pk}>7 z_JaS>g@L9g?$uomI<>N9!@`XpNDB)d)sLacpJvqBvx3=)Hix!>16Gu7b?|HOL-Tq5 zn=L=KL0-t51?EHyu;Z+0>s_-RFg;$gq|n?v(aywwOoTuV$mV66?j6R!`{~#i`H%eh zkD|O2dwCJvt#H?8Oq^Q+5#wgWyrMh_2X|FQMwUXZ)b7j{@ambm66cpLfp%Ae+PC+1 zt&Nw2Ul-DL)+;T(3=Z$eMWyS1o}wzUb@Z9ftX&xzlJyt3PuweYH9N+=nxIU|i`W=Y z_TtREy{qAIoW-n8{7Xuhl4R7^=uA9k=lrTEJ=yx1|I67%xg6CRhG84}TxcswQ5~(M z-bn^PW`#Ml&s|v(#uTSDa5ny6sXfaqzog4KUJXnWKvI!ZL|U7vgs9tD*IHK^Ll*gJ zHl@6e5YL=rE-u-VK%TFc_Rx++?O`~cSRIatWE_gYfSyfA5K%oM?S8X$ei^P|A|_&R zW7+R&Qm$WtSzjE}B$8EOtYQ+O0Aelf8f(Z}iFT#9uL)L*HWMC+wO6l;TBXxfi(|OZ zgvPOM4-^&X4Fi1W9?DIk?Na5WH{g)AZ0vE`MPb*p9_%#k{AbL27ze}%3RQ|Ge`)t$ zz@r-oYpw&kGeO_nqa&RzYango0s1_4&N=JQ*5u4#LyX5;LOZ--4!IyHY7>HZ_M*2F zemqdsd|6_JgM{;i1G5PxQ;)b3C8_z9)W{*>%h#uecE(Tjl(3c$2vfq%)a!)fV#J1D zZ3sJEXRAr?88dX&%`;_?Ri)@Q^x{KgrPpy?PG`WYu zN0^|cJcn8r#N1{!68-$zDuVSj&bC}!V$vg}EECDkH3klLM5+%Su_Q0dhMI~MSBbAo z=-R;(mnP~>v{RA~+mN;i0DVB3<)b>oazkn+b(4S?BYm>q_tLl{J)#fs>#=~;P=6dk ztJQJ%4UjNxWA$N}3Tk4K$1@|o$4GUFIi!SBPp)s=k!|)9;r%oCQVO3(zWVD>)%;6W z@Jy;AG7u>vCVqIg$Kx58)83BP71*kkh~}&7*xnEQ8PfJO9*Q8ml%+uP%5N&!y0{V5 zV?h75MaS+12eGMlX5YmdWr!mUIvkv>7~!xy2{#$4f=~z($g_N?rHcQT=Akw9F$)1W z+~Un<*ek1ReIC_^h~f=m!|nO}rdrP`e$ijqqulIp+Y+TmyopFYEnDqNdHJL+lE*eM zOBVnE9Rz=4ts^u#_mp3co3z*9e3!QS;e8^J86H{v!FrI?m_$*b)`g1~3H$BH!If>= z*9Ps)CygdwuQb@$QiB_31PPI5f$_L~qN24y>58!Z~f^Kn`1%dWrS_P%0# zaeY$H?4pFgJyr>1-d$ltE@2)%SQF%V`ChYjWxWT5p(o^1RW>Lu=x=F*$k+V>L?!dD zkqqaKo2RGC4gdy7>yu+0-Fll}=*9ctA5Wr4k^3w~HUtArL`;E7DH%BdIPzoI(Y zjp#x$>$)-t$cl4NuaIUgiQLTIisNwlJZ|6`?n#=FoY@}IeBLgFxSdC;W}}{mEL_@` zu|FDrXL7wZC@tH&)c@gfTP@bHp06a*+>zs$StGnA?>*|wbESZ~+8xH}m!$f{Y*(kP z+e1$E+Vyi65i)TZjGZqdP}ON>OU2YRw36`;qGVo#w~;_ex9-|D^`RH9KWOhpD1x+V z4-{v_Sb15TE(Jw991-t@z7ntSfV~&+7Igbe$g-+1uMr$bVx?x-qcr= z_9wA-$}xV}Jz2YFVVXntQY=}JsD=c`i0B4@`0~1_T^V&L0v6hDo)rN4Rvx8q^|VpX z?JAD3#?O1Y?VNOK<(wf29-jrR-7U;x<{KivQuKflV18u4&6?rY(zCIvsmJ{eyPoZH z#Dcu=8%=x&3!fHahNo=?;T>JXh{xKqnyh&?RgjriQ447o%=-kFp|*=FH=`X|yww7U z&UJ)CDLQ~Sdyg~7vv$&_Na0{fy!w%K{~hLnCf(kQo&vS!L0UJ?t!5!{PeZ`2f*Sa{ zkze{ke&KALAY&{y8{5mB1=HBZqkWI{oKI!MjYBpds+FjsnrpJbx(@oBO zzK2k2L{as9Q-t=7Q@jsdj)8iHtgihd$o;c3s!~QW$GN>)JncQiE9K6I)_J|&`6^?H zS6xX5BK)~vd$2GMtVXJ^2*HL!sfmc(jL)4ct$oRpd82n*PE@WVdwS)dw1`f>$Uves zA`mEYivv}?bDr03KAiZG9|8Tv6XjyW(=io0+penfsubJoh6BgG?D1%fQaI#f{T}uvz?TqH>8b?{i_ zL@oejOVw|kK_I8oQ%F$Kx-hxfkjQ!$Fg0EG!J(VBmTc{;`xxPus5V#sR)E0qdXL~#do>`n_KG$ z!-9Z7;85Rm!Y_4--DI}|WfoJd%EJ4icQxM5BrYm_;K~c8WlF|2!#rARAyPaa z|9JDPK^3)}m{1_xpz#e-o<+r{Xp>Zz>m(boR#c6M8}5|`m}Qa}F_-2gwkEP-jXIAV z8Rav?t_C*xa;M#y^)aEor$M@}{hemcsbnRB6hi4q?mbU`Hun;4H~5fty5?BEt+0qo za$~v>D%tnVG;g^~f8cM};O!_rb!#HIRr%p!tyu(_HzFc&g*oCIs7Z1mJ8OrlfdP9B zfiI9x6EiFsyuIzq6?sOZNW+)&ZMaCX%c~YC)n01|X|8J(u~knehmO!R3L?|8?X72l zaN3(oRVmJb{0LT~teS=FK%xI=iy|lL?%Xrg35BS_%T}SGB(67Aat4ad6nYEM7SM|9 zGvn|HPD!WCk8EYL+Farhuxi>ovK&TnO`W)t<|=#y3SQ7x5ZFFLS=Fb=KxCrnO<) zL96+xzaN*WX+FLyd7uYlBmqimfRF*g+suBC!FoyOD@%GWz^(F5^16!80qsHRvK?nT z`UIN@zn0k!i{~oGI6nP$gvv;o@I`D3(og``vjXWmO?>u1eHUtDyq_Lz|7)C5k$Ah= zIHSq&h05pB%43P_suazrfvD3W8uZm66sFa8lglyPA3^q2c_pQbJE{pgOef_`vdq| z&Mi-oQj>lAttL3~UR+6)XpbeG;I0Y4_|N~-Y5o;!J}Cr2B6fr5mo8Zmqx+)Prd-`l zG%PwJQV-PLl7sLNCgg_Mj`xwvSvd<8O{OQzSBsM(Sc<@IDQm&=N#cn~TXS3xt(}em a`{DEOC0N%C5E!iHUZ3Ll-`bv!vcCanGPzO! literal 17322 zcmce-bx>U27AD%bLxA8K8rLAfo!}4>2o~HSxI=If+#7ci+}&LpCpa{X(-5T5;5Pj3 zy)$!Pz4yobG4)PYojO&$`|Q1!d~2<5uYO*7-T=H(l2eoeARr(B6ka}n=M{hq00jvd z85s%X<%WWSf{KQXj`nilVPRol&jO+)?M1HeKD)FT`rAm9U@2LO}+1OOt^KXd%A3k4Mk84VEu{pDNvR{%r=WCSF{ zmno6ZF_4h}2#81kWE6Z+WB>k zFA(JFf?j(oKRuwJPimO0H{?>!3(fdzGYRgm(Ik(>B$D5vD>t^+I?Pa|wkiUW9M)jM z6k2GSVzI-z=2{+g*;Wcn*7R;@4GLXkQdh#_xXr)N#ifwGCROTNN{=xhO;`wzCybts z&9r~&`4Fvl7%9L-13Ud@mSe?DR8eX}i#oenLl+j1c~nJbkYo9fEYjlw0oRQ6{}SOm85S^E zurT|}rST<%0XZqU+1he3WEo9w;a*PjY75{&WDFTY*; zkXvHFvB|Rpnc|ceT*fN$N=}354w-wUC2Uf*kl!U3aAnEF4gI~~Kqig${SCIEZLaZg zd#lX-4A7=Aq3Y%fk#BCT?r=m2+P}vBP8BwYZg*rqjj6LtoZ?bO(o$et-+Ik7a;_-8 z+cJ7anejBipfNSCKW=b%QLiqht()nwyCDIsAkJ6 zkp8Xpf=9my6UrzP^|L5 zVPvkXsSo+053Wh+0^dj<@8 z4Y9GTuAZfQggW}He?wh=x@1tQHW!MeUS#r#iEykV(YaYrUOKX;&sssH8l`EV$%=b} zx)hAtAS_bJaZ0QaxyWf=^( z^sW93a}j0%gd*GFo(QkSsBU=UQiy~+SP{lY?hcUP&?1vAW8hX5pvtMEv-Z4_3nw!8 z)9M?B@jb=BzS^7n!BEc8J@O+qeq4=$AN`YYuFF@r64(`4fPZF1e;`qY^)*PqpuCiv z_WLnYIB2+lXhsOP!9ixv9ZX!6Q6B`%z2)PIKcvZ=)$S$I0rjS4T3F z*;1KuX8%NW4{I_#sLxxRz~}Csx|M5EKqh=Q;<(<+=& z!4&=$zb|XOqBg%sXN%IS@Jo?6U3NHj_$W%EofGvSc)^}fqQyCyZ{mw|X=3DZ%VOMB z2VSDITTC1;(&6M1V?pVfDy3uA;*{&a0r*$2mM$VJ|JZJ>+@cWIi~yNE&f6Bm1q_ps z1QZ<#81j^Y?=*&m5an0>EX-K9hp+Y0Q>Vj&t~qFF<{l5UzdMfy5N_yTQ-JMCzt=>* zVFMuY$V1NU?g=jl0yMs4>4D771+p?ZZi+tYnXNGw*S81ZvtK&YUo(BZ&iR{88cHAh zu8=kovPI_QN>(K~y|&Y-Cz57w1vUT^R3_fBFds4VukN#QuQYXbsiT)Tt-Bs?+ z*TfrC?Powyw!|~wT=f|M)bonwRKhN7m~p}-G>@R=aJKcIV3)Zc7RT|BUug;8*8Za2 z>aXI@egDw}zM_Ulxphqv7}Pa#yr+{IP+g3i`}?V}enDJ~iG0*~aRi(e6{&Jf6ScPe zCVdUqj#_F$@*hr;7sKe_@9L`NM3a_D) z>&fzqG8fzwGocbnKQFnyuNN!aS|qyq?u_u$UG79Gu)kLyK5OTn`Uyoadpwn&4sF)? z7BGHxGqs!=WyN(|^$TOjj7QbFlztkWqub{#&iO<3k6UC% z8L3m{XYbjQ=Z}{4lc#&{E()}zPu1v$yvT@r2$=9_-3>gd16O_{{wgm0AsLQHvRO(uu z0j__`B^buH?Jm-Axg7)iqCmDc`O!q%TOFit}>*)0i`V3K?&xb=p<@`)!F$20$qM9L< zF*T-rY+fnvj>Q#!Q|}%4H%jd;9gf_!4?58GujwyCA{=uziPFjcaP`Yj2qMcx5KYW= zx!1$%6sq-qO=7JzI8$o2c{d4UMcy{=G&o@iN=cZV;pIsjO&w};ngloQK+S+bo!-FE zE7Yj@=>086_k@6f3M9Na%*8%1r(>+V)oATUMNgz(hBABz8zAd+U5hy5;E>?q;?)m9 zvn+_tffoM?XcAxcL_aX>!?J(xwvL-KaVNKBTcNZ&FnVN|p7SzetRT{%Fg3oHNISCi z{H4_lj4;@9tb{eU?+EEt)DN-}xyLGE_K?!hh7;iD@7QX-2^)g$v`$R1cZcfGH8$Eg zu##!30IjrgO!1O0$pr&xR>NtwemDLj&yNq&t>kC6S6hmvx4k5{X`RS@dS-q-_}vtX z@dB0c0v{hpVpkz99~8%2rrC8Vk80@pYg9WM@UP^41j{>{v4_16IrPLx$Tx6&Pi)55^YAlu9*c}5Mcn@{Fp zt-Zy-P0}xkvqkioZnbnbX56g_5<^5O- zij5O4SEVim80?3g@ZAoue|WRQuRcssq$vDiH4AfAb9-z%qsdcP8%BMv?h_lQU2@cf zg^>yp;%TvSzaJ|6RAmoR&$|!Qp+eHhjXrB01E={*MA#)6bY*YWM~`RcD;F<1{XP}F zA==uGa!%*jh1_XaF42#deowDQ{&f04=I%Nsl)TQj>2=FqFZ3qB z(gA{1t>P>CHf{8FT0C7`)Ee*I9DFkrf)0mo#gdwRET#NPq$ytm@IowmAS81hD(c!6 zkFJbX+-V=2B(bg;hCxt)%lRJ0dG|k6?aB`pmF5MHyho#JYZ`^{(fXshW&KP@JUxro ztBzN1;cLyiUFK(L&|i}2qJD+vdXuTIDhJqTV?T@NlC)606(1LkEYQ&z>-QF5fQcru zm=N=F1R%=NmGfQ^g<@Hg`S{Cf^Uxq4jq6+mN?p$+3w@9o_l`HZu&a_2QE|(-wWx0l zm#Zp1h+)&`N!WNv1gYlNhK5G_gn+g6&5eqd-a9&G=+yG2lRF5o^a%QUkm9?~Fi}UY zfH`&=U~_wYY(I&u@asgsGi9LcDo^YTRp)CEcb$lzw7=umAGM=qt&5sByUHlSYnyM= zO3yi5$>Vd0+z13eAK^&kAJ%H@V2~K2Qq;MO+fUI?Zk^#_mxkB0MKSoXai(l3WOkuz z``;yKX-TN5%AZkvDgE!two5ony7EXaF8ebc!6q~yisgZNV%rtjW36JgWAaU$f0k4h z2;V1D(rcOi2&D4JVrc$89)XT6PzKQdjrNFU&+vb~Ex&e&l3(PFbA0X4mz7m$3a?5+ z;0bhaWooIGU49PaH=!_-yhDdUE*VJ`h~48&>tpu}(0B%rU;cYF+UcJdV)7G+rbE+p z?C|d*FO{^@rh$e_)1N#(Bw5>K=RO|J|EDC^G@U9DWxNo$VW?&tnn*m>ym_DM5)78H zJ?2;Iv~w}E9k`jl?2HYFH#Id&E6_!S8W11s9yum+di4kGwzFX&euUCIZoDUT#SR{2D3$uL0nC~7Xl@PG zgoBvKfo4aa%MnFi^+c}!n`2-l9!ox38NV=w*VQP-qiW*gJizf0HMJIPd595$L7(pr zib@SR86(%$lp~<@pWb+y9hISul&o``;4#s!x2>GZO)Y4;AnWN`OXl(dIMvKX9l#`D!u?_SS#v|A3qiv^O|MVt@_M#a_-oL1&!c+Uz7s!~qE1y#@CK&@r4zncBu9j&DVbA?LPJA3&)OR4Z z0KSY<=NbFim5}N~#phdtm{+*r$aY9I{QLjtl}#GVVdI|MIfy{7o#GcmuU-`;w*pTH z==c#zL8EExChIu!GAR0Kie~@Gd7e5?V7!%6>P9m^*Rr`*ilFW>Ig0yH;lcok{yT8a zPPG+bb3k{(k5-~rw7WGZ>>_Rpd|6F4CbA|qyfv`CRrkucQhL=_-!qD|wM7f{-6J$u zvcx~YqzUVYWVhPMK`6E7GA?x|l`x)EbiOh?@?h)IQoI9E)5WLIMnw0En4OgAUEL83x4|h~%rXDrQf$|8?MoXbOz`2K5{ZH}_NDfQTfo@l&E{WD~RU`Q|5pKytL!SXShWk2P$M6bvSVu=XZF7iELBeEk zT&HA!UTj5ELy^eh2syFiV^EdaZ~c=i$3lIX zvm=Yk+8u52iq!KTnm#LF?;dlGwb%7mZ1_7@H}E%c)`iX|kgfx^So|NGGYX%))oAO$ zY7h*ICISpgdpuz*>bPfsyV6U7yt%uz=&Ou}?d`k!Z9xy)Uu)GK7VPt2mQ81TCf5o#k;O^6l&8qtzV` zuS3!=wh_6@pZ&j&NJyfA6p$vgSF$;>mSnPf%xkV%+mjj_ayIiFMzRQByRz1D*Z&DR z7HKF5mK1lUl4|Ze-~gkxZ+k;O$gweY#oA8W(YE{$l%B3r8SF@?T{(XS5H|@B`pkc~ zorR_fqz1YS*QPMfCoS>GD&VnlL3#QxD(rX=fw1G9&YQyMAN- z4#G#<_@^l{YTmv(fza)AmChmwu2=&3#D?~4Etn&@md)#*V&l3?SJu=pR&C?RuCpdL z{bu&Dftxcgq8X6{$kl>V=@Efn#EbiTu^o6M)Yz_5PRhH65ASkb@N0Xp1@JlX@mFK% zl}QvBC+J4;mD9E5Vk*+2&=}+Z6YLse=1)g;cW7;wK}e2vr)fkKHT|@fO!IQuI5l0q z9xg{<~MOdo4OF{(5S0>FSl%Bz%fHoRHC!-wB<4*pO}3F{}`I0YY$mSgghcPv=v z9oOs}zt(Ts)i>r-%Q`^R+#O( z8BLq3qOAde!}?2Rr{@zwwqQlj8gmbDfxb2?81V?6X-tK2T4xQ( z_hN-u&kJVso{Dwf5}~i`q35JsaxOqr#(mx)>kGK*-nJrC7Q8s4ED)FrTX(0)o%4Gz zveBv5g25T|gnowZ{KQvemwT_Z7(7a*jQ9KPI#knVqxBhZuR0X)_u(0Ek(@C)n^Cx< z>N70(u+ZC^BIuMTchvo{!Xb_2Q7IRug{j51;&ptq*N|MKvm2A|iakLX*#Ex#CwB2v zsQL@3NnY{{Nc*IN-Ju%@7z+9#P-tgsZS2C zz3>vk!r~j6u2-3@8RAivXHBtM7ktM3D?!I6?mW`+PugY=LFKXjV+@YH=cLE%(|S9< zhDLDNN2-Q*@rlB7e+$^V7_zy8FPhUhDKa;2tH?%8Kv*-^S7iZ&wpxmepWEf={Nk}Z zscJ-M7j2t58(Z7j8bd$q+f$ik`zXLIKW^B?&?E<)T~~Q8COPQ-&XK`I^|p*Xco!)j z%4qqvu&)yQ^@eFjS3WMF=&8+jhE!;eZBsuWMiplg3neK`AhEMCo)7;)Py)xHf}boy zy5h4pW~MOYTBqgMWr}w$y25e{=;$6j*eBBn4Mh$0E`%R{?<*+YyzixLHmgpk%U>j) zSlJa0iEmD&iErqx8_-ihG8B=fZ2=kfoZh_8bt3R3cwKUQjg;+Ts5i))+b~7y6RR8+ zhc{T$>{;%1U0Z?Y<1~gI)seH{9?fe_6*CuR$V2?7xJlqdjmECqk)ySXv%|hJ*Em?< zQV`=ze|XNuY?O@T3GzkC*Crl=tX@jsf%h6m(vOagk1MNE;M=ikYSW=|8Tf_nj4rR5 zo~)sPqTIu`W46bIrp_e%WayT}sf!-!3@c?k@6Tj_s^~MxX3T%bJbefYg`wJT?jc^$ z_I$DqaYUwEU<;>dv#tmZIG;z!4DDE&Dj-TTGT_{A)H$t}GM#en#M9mtmq%rAItN4j zv+;S~3|tIYWV{V_2&WG;`IdM?qRqv$4OJ7dA~uVWb4+C|NQ#x$mU19{MIX^aXN+OE zHvLQIbu`RRuJ01*%5ax=MbX|SB@DqWb!jl3BxoA5_S3^B!m34Tbl~Au!MFC|v#jsG zNPOkv6sh$R_PF=s8V2P&o$3TNEtB!mQr;@_N?03sivV58#?7Gt+%pXASx{`h*{!`? zbND(HYdOgGVd|BUgf>yYj+e6NW$?0*Dkzj8UV*%y!P$bb2PmgNMH&}>YZ3UugX1?MsXrXY zB2jLAV?@I){HHEb&aRwame)Pi4MmS<@M-p36sg+KT749j=ziXjIivC%rce1cG`Pv_Sc1&~y|1X=ob|tez#S8w6eOvUc zC!LNf{mXb$ytuy8?D`NSp{jZHLf9B$zfB+5WgagRkx8cN@oT=SAoKd-Ytbifk2C%g zYw-1YLV-~_l7DuOw#5D&xbMykrJ9ep$R8@0#lCIUbzxE(!YIwbj z{73n(425U(zQZeXGwvNnQ@iBzn%j`UvBrVCr=KvBM^@=;9A}$}pF2w}KXOSBbzsIP znC3A$4ZJuv9AuuzUJR7?Nt*n^RhMWn0&Km#F?yTUblr`eVyAbtV}DfY!dH=4EZDyr zzrA`2&7bVzN?2NPJDvf*a~eyoi(~V1w;Oxilt6mSn51~C;J>acGh%6snw0Qc2q(H# z<#XCjvRUL#_%Q|rBSC;y{-p)vWPi20C_2RTSa(oJ+_8+(2!v6yO$aV*0kK)xh$()t z9J zQU3*2cYUp^H>?vp&3D0DU}eakkoNml0@Qr!<&!R?`3ztg9Rjw2@XajH0zCUy*of&9 zw!7qJ`8$J$Bz`2C{k!1(Xug^LP|X~X>~n zv*97{m0d?UzE7%(XAdk^=;%a!KQ51f5Jel2A2}|mo<;bg@HcIIzZR#EyX!EVx3wdr z6E-Q%0<2N?(ZF_U2tuGZaab5G|m zq4$ty5!>LAW?(@7XvO-&R@rfXPKeDP3A`~5Rdoib*tfKUEZN~)MVdGJCD9| zm!E(G0j71gPrAFucLXLVm8roYn1v(9!tW9buSlj6Vfy6^ zR$>LM&^%k6;%G44axZ}lo(m1DP6fy#8~?x?aHP4GtTeV^gf~|AC}lf#KRhH%@9SPz z3ff$Ns`N33JHE*PBD`lt`d7VY6zk%u}V+Eo0vP1+V1S+4@t1#b}aBhz-$=Yks{1AZ>p(1 zm|8^XZ^$+~jz`nt%rKvvWaCaNLhwNj7I2+a`_yPwM(-7Oq-7*A#>vGaq`=dPoDHwE z#~LGZ4+iLm-sVS`EuH@S+$si^t8?q>yFK?UPd2aEYe*GI!m)jQmbRlONT*w~dI{8j@MBn^6FOHXApnroY&;wwo0U%5aoP?B4o@$*BGznz22+k0)7 zmBr=+=)aqie}CWma_Zk0l`lA+!&~HqGHLEYhCOg|V>ztq<4!+mWsbe!NtlNLlrik9 zAbTl##LxRijnwtpZ~a8PObvI9McZ@pWT};De&f_L^Su7IYg?PtEu@Z zCLzth1Fgm+L77hLgv&=(mquP8Hg<7N-xN7VIS9(J2b|sSFW^(%xFUcV86*iv)P-Bv zz)jOPK&S~wqy>xL2m~MydOIH~0*NmpO1%lKHzF9vPt;p`)6%=lgl%$2P}WDa{bT>+ z3kn|_5O`&Wok=-LFML7ALZ;?Gu9v9x|D}VLsW;DnTLHM^P;G9`V%XVzRz=v2?BFaq zk=ji#?&eih(A*D_G<2sg*scQcOC44a1+w-klDn0Hd8&mR-mG1H9z(7OQdo0W22CvJ zU~`)hIh|_zV{Ow|S<>uNWKDWjk$Jrqq)WgbJO4H!`Aweab#?C)7ij}- z*Z$zT^*UA=$_t{$FN%G&cXzCKu*xpwIczN?A*+~<`mO{=9<1n4>sJ!)E7G&++UZnV)^Yfe>Y`rZNYUhVvMYN;x|yIX%CHWf$yi-@@W%v zmGEnNZh$AUZ<2W8r6`d7Y~lZ>{FGIndH$$L?ed$~Dcw|#l>5kc1Weyh;j{!If#Pmz zwUq~vF#YdjLr(C-(eV+IlK4T&uZFrX56lF){4o)30Kzli(DD~y<_dVMR$uEGfS*u1 z3Il??=mJ#dO=1+SEah1Fq~8E?Gh`RQ8I7UpgEENyP{-jiP7k?$_uN^FN+Emce_FaRCS4ToeYjHy4dMgBe-R_y#iL!8 zr5-K+@E7>K#%A*k(%#lU(#-K@mT!^W(sr8ou2Ft?Xmk~cb}2(vM;Ay)}a<% z%VBSevfu#|kN%zZD^LLHIXMFeh$NEeSAV!jU3SLS9FjX<*xF^&6FNC2CddOjB7O%x z3?TYjAjpGs{P^X%*fJqI-i#vnxHOzKJ4fbS4V^xwZ>)b8r`-K?ID?QPx_O|*K$54U zP%KP~Rr5P(YJb%=-Rn@jG}oDZHq0DeI<}bb%46EBsQxv)T>1*7^TXJ{SnTR}X*KDXX zPM$_K{|}d}lJ5`YjdpJ@x#d)iZw9aJ7StxvFhI*+yXp?wz51pCL9WQja}te9${ojO zH9O9SQD)4tSgkUW5ZAryob2}WbfvqxOQ6|$A3fZ;4Kw2!_1EaE@VDU+ojIhjNUu#< z8yfdjjq4Z8Qktn_|0TUzB(bWLf$4AInC6jSF2e0SJ{}+krIHt( z_o}>QeD4KH4}v{}xOg42s%m?OFK->#<=sBAA?+|{wrxt1i`5(@)|J1D+ZBvj^Gn?1 zuH3S4UPxFSJCM|3in)Z*EH~ERtprSIQddAQgiB{tLe;%CwwDge@1r7^jBBE%!mwmD7hWcl49iIGESji`ZbW{xnQ&5dJcRUS8WH9AYIX5ntN4n|6x#doa}8p6*CsNW?35cbVh9o=tuy;T-g?4z>FKM_T5yC-d-B*tV`p zW0kmm^p%#Y{-729ER^Ss33)aZF;J%Wx=LjJs7kHGfTzz%MA_DDN3=}uQ*tACjbp~l zZXk;OyrcB|$BYNkgKR_RyvHudLC1FZS%lY<@YSf#i8@No+bM+RU`2mciC-_Ii=@J` zjjzWaPoG^OT{Ql!8`}qJ@w$);c6jA7WhGx@GDX;G-_4VVmLx@ zYy4NyF!HXZn66-jDvkp9lZ8RteQ18S4e1B@knNieZU=`?8Fl@lk=;B!GvDgH@4tN< zF(6@J7;NI#u4yD8!ORGz3Q-65EZ;HT<>sMD-fIdJrL*wp6Z>8P(PEohM=dzb@$6{| zi8gLRBOG-c14I=-w&n?Wc(s+22|FC;4_~&qQ8|ZtR$-2QiX6Zr<%QVI%^Q$dZ~Xhp z9S{PNmE^nz8?Na|qI={;>@~ll6LH!KjmpD18CB!>D8e1hMM=OG-r}Ih|7H7C-4a^C z&0TWWPH@qCIqIMgI=sJdSaV=#h-A3Y)O;Av0R^L<@~_Fo7Jrj618o_gZGRHTr`k!8 ze*atn31oEC|*Ak?Z+;fiNY46$0cr$q|$Q8I7WgQruS>$q@gq(k^CM(FodO$A&0vij7nOdy~lra`7qM|7G z&AS=l)h~K?u7`RjC*JXLQq-eI3Xsqa7YZ_$v{LR7Rxr^BK9qo8h5qSG{$Q+eeYQP7 zl?-zaS|91+yZ%|Qo1V~Qh$&9lnO`nx`9nVrLR$?pyCtB|hj0Yu3s>DFSjz6)H`*7D48)n>Yu9kSax0gjj~ zU{9T#k6aNPyL`W&l&b%Nm>%$TYNi6^nrd@e&thCPrMZ@6LNBd zCes%=^}QcDGFpB6L5rJeE(QpdReEIK91^~KfUpYnEQb~lj4I|k_qmY*UZglHSvG>Lwu zp|0Jb49nea^$0+{QQ%Iq;jlm2ZBaDZs-ou44w9*w)A@u%MSSdn+VG`!!bu$luW;ms zQ+oTV{_X&e0EyH$b*Tkt2?!Kcf8PinA5XKgIzV3fn^??)HasOK2;+4zs6Qvat@&7Ad&Z zEKmLs2uB9x_e9@V2Jm*sV=&?FFK{8mg`Hv4>imIE%lx>Bzz4AFcD84l=g!^VYxxre z3sR~i=&Ll7VU4A8xe%7nDm%Yfd4*S3f*C0u87ZFp1(v0exO?|=@Rl|cms|9kUIbl} zN8(&^QM8=71edkB29x#OwNGRX)2ksa{#ID9ZFh8QBs%OzTbfQEJ$_TWPX|W9+3HGt zrybJAPKT|I@0q8WT`Xa3;4gaShzvW6983X4bwD2GPY+6kisC(wwnilj5A9+@E|hd* zh&9+%!XmmnHPF1~KqXyHNsZ!L;CUN`=>Z?WlH}f+uBj;!g1#VVDw*)>o;4^B?Qg|@ z`}K>elS@;Ew$B$-EfVX4&zIBn!91vjFyqlRzi?x8@r^)PbMW%=_R+7pnug$+1tNGP zDG)teyjMO@BH8M{0OEf`#sAk6Ywi384d*^DYK`RHNV4}s(hf|`H*t#5W*0ddI;z-xbf#_I|lPOY;|ZHPEAiNp%d z8=+KB$5&h8)}()_UJl4M(qRfkv2r2z446Gh_NZ{j+F6$AJN)aH?V2^`V5PlZxRB>F zkEiCxE|!-cev~2t{YMULl)SyXqY0{%F<4hhL_y3D9LTy^kw*Ss;EC@4-U}P$MK3@` zbb5fQkI4rC3SO3b9J$w0D%2~NS{V7I990ylk5)0Q=n5R7*?QNgMb8!VK&!b>tLp5C zUDMBZCQpB){@MyK+^hU5|G*97_Z^~ta>_3uT>n|TcOm(*@8+N6G>EbbsJaRYwOdzq zX<1og>Ind+qM$t3cB|tM`)=COnywpRwhW)0t(czJbA*$P`=8r>y!~pn2lY zJb?|&**}^Sdkp(~_&Y)(-Y;DS3nRAJQ7h|Qzsfm=x-xH*CHi!3ao?FASbp1!bt->I z4j#HC=7ub>fB(9uFD%fWql98Xwi}h!rB0QHq3VYascA@Xzu^XfC}wJ*IxxHCM9XZDF_mnj zlCdoYa7OeOMNHMNQwPfZXkvB@_Wbsbkn~flv!)0w+=slvuHtCS__f9CiM9AU6>`{K zYvmX|#(&b6EAgz1K);(m3EZf-xuzjbd#GIaFyj;PLW$Yn_vq@kDBJ@+k>V5rY&9j6 zDld7(T)~2xPfERJzEPxbJ7TFu+tM8!6J$U6QaU@oW<8y96)W5=@vLjV6t8>Foa1xX z3n4&9{jh%pxOM4qqzq4O7C7(w?c#$X2P~Jo-BM8HH7_3P^c_b8F&_JjZdkbg`U5>WG1PQl_cB6 zR~vhaO#D-3(x6F9-~BADY)F_wU^qOVqLye zBnxfws@v#`bZvIen!KY`CwZQ@)FGRgf@tsk($%Jerjfmo5k(^d_YE!Yxd*YmZID~- z44^mdw}^*xDHxdFGVVM}6yH4w7iX%>`Bt&$CA%$%LT4KA4RVNtN7(bI|D-s^&ME-x z^Gy5YP1{yKS1j!c=leMCK+X=iEpyG>od|R#$_8qzAB!Bj(~aQrl+OU61&5wYEexKB zxjF(J=Ec@!c_c%a4j%11TBdBw$fd-0m;LjS2ZP#ltbnMichnGGolsEGxaGf}dIxuSI?4Gl={>E>CSMCs` z0`DvKDUGU2?+6&VTv8U&Mg4_!jx%hp5=^b(1zwA_J`o9XMA5cA$3Ux|Y!iZ|2NDzCB{#^ziwf&o$c3g0_oHUhR%ed=Q{A zVCK-Oyyt?}U=#*L+Jmo6-`pRld60GeAtL zG8&ecc?LAB6W2%_+&lxQI!9?FE>j5@+8iI}GL04r+SW5BOioE{{e=<}X;R1b?9M#u z*gp;D&nLRnvz6Piv#@zBCR|<#WD1%GFVm$Q*cN}<)s+r)+*=|F*g&x;L(>D_F-0q^ z9{^!3-#H%}91^qxPU-@H`F|@NA>YPW68_-vH1y|>cq}a#hx1R6*Vq7P)ISO?E>!x7 z`r4jSi!Rv!6De*iS{lqnG+sqWAKrbfuaZ5v5J+ohWW#Adap2shO)IwdpgSUBC$Z~^ zqnf@E8!GK2hmo=zlZfHrF(jqQKcUQizqByhXWI6LbWEVFB1FCrL6Pqm32ZqAn~QK8gln zn*P&k*`uBt0y!r{??#0<4cDvmWulf~)pUdFQeIrsCO`70UOSJ5lm0YTBIe3J8maUy zHzYh0IaW>>#UAW8(bigMGE6cE{mfoBvPRO~+^&S+=w`OV6Er1Zd&tzWZ;?2TF|cBF zcCjnoH?D#l9U!Tccm>EcF^C~EK*~qeqCln!W^MNjNSN5>P1xYXsSU#B7nlxlWf=Xw zo;N!BQQ#}#H31$uYI!FgOnbR2%=t;tJ;2#7b!__=>u3V7bEPUW4+8>l0$V6(s_PjO zYGJij3xYpXJJ#RcM;%5xh%1+m2Zl1I-8rXve8yPHhY2I|cNNTAMZ~&FC~#-RF^%qi zkYQxt+Kb$pgno&@_Af3vs4WknfQwb3FI<1oewgA&+A&}pHseVZgmQ5OD>_uZ|Ja=S zU4R?+T3Sk&)y?;2Aua2sMzf3NX!?KIys) zptshn4vQ`vGqO10xfak8m?h(Gs*@!M=lCn~VIGOwtceP2qAPV9#Mu9tI{y33)JG_e zG_=V`c_?T;`&U=k6K%k*hkqkp5(vyRVfWVg%nqI-K$3j}`|ETLMuBg7Shje#iy$kV zi)OCvI6tmW&GOu0seE9(U<-H~Iys5fk|T4J@GHga5x|%yY!=tiy!lbzJm1!?NVj0b zhE@_g^L9-KQgHJiO6;pfbAWrilovLa?fbaqRQ%$g zj}vY>LnPrQwO-|vz9~(JIu#roN0rpFQ7!u(uc!y8VJk{vBb7%*a_f{FeqP79TMp_F z3F&3nY;&zckI!e@v8E(Hh01?I(m&{2VkkAlEz~|{dtcBJ{4`H4$i=~FxaKBDTj%)` z`6p{ZDdY&LGSJox?fcN?j|q#{p*F+3P&M8OzP7dY0vDa2`n(UJe*r|w&7jkq=LH=ph+glKuUIba^Gr(tWD0r{FZH!Bg@KvjFMCH%kDL+tnEw`tt zOR05P&aft33k5mH@yb~_iNyAH(&7)J^rC8TfV-M_<n(BkmkCwR2ZF;NpCCxKhDHLz2oP z>bes5r87?HTU z)82!xOAbv!Z&X}p#@5?~MXJuVAAkB3%oi7%WaSq3^(eJjZ?o7vNwkJ$ZY5m(tkIQGXt_c@LePuxEW%qC`y0ZfCs7Om zVpt42>tm4&2(FQMg3+ViExxivTrRK25)`rohzkf>SdJ7dyKL?_Yc|t~fjYmXjikr5wL2B&Pl*7Jq59zcrntg(kr# zQ@&Yl?hV+S8xe`9b60p5cWnq!p)ewK>3hwn)YE%SPm0+!vDjukHX&Yavv?csZ7$c- z;jY!l>^uBJaPjTC?Pkrd&7l;wUjO>r#!FcXCJ*`ikq>~>1Mp+Z`YRJ^aWjOzWR*gj zwsL=Sdd3!Z)C4-vR@8~3I*QK2bN^J0lBRWH;H=vvJ9}8g{bNYYYgzYAb2~`rSG26{ z;hWrxwxb5aJb<6yrJiUVSZK^i<1EwHeF%xlhpf4~JwqqDZ#4)vo?T@_gmbYq^2hNK zwp^J&8KnRl=mQoWz7BgVngC^}7Er4xc)?4jta%sM383~pT(VohzDDe7J?N$$(V=&C zUhx-mj~u>MZ1NgkIH#wjVd|bW9`|g-=Ep#-5Yrj8c-w);mNbY(aMc7F-P)(>52C;) zFB2x1TrI(KTJ(6hClhHecX&6rW%3nfZ>(Ssr)88+=)p{2e5hrTBzzb_vTXe38r-rC z)uK?o1*S5%@+~f70M~8xr+~^_4$pueyvatVCFMlwZRM*;aTxviS|&NetcD&hz1CU4 zh1`HmrznhPz)*Rm+{Rm-ENtxrQ}%^(Etw zsC0!NLrw5f`5QIUA#dZCWsQpkE%bK(&$OE56IlnVKG?@QlGy9$??Q~DV=})Ddrda= zjOVY>f|G`9@|i8suXZ1dxU)x$DzEf7O$bLQ4YQmFF^OX74HENomq5PtwKFlyb?V7K zGy6;~;^k4T$gJWcnJ$6KX#{jXk{3u?fW-1ML1YKJ0;@a2MoHUQ$m zo&OFY$+5fWoc23St)W6cla&-HBgqWZg3-m%CB6o>ZaJ-5v2!`Nkl>lcw328E!eaSC zG2S!M`fFl;MzzHJ@SR)S?WdjFKanUo_-=CMX}Pkx(zq=soa6O$n7Y)`lKB6ajR0>t z1LFb{R#0WN>OX_}W912VTxGuny-8+1sI%is2&eX8{Z%d7^9xv;zBAOcFdxvpn*cN~ zZ0<8-`O<3EKiAAA-H+Dh&Yr0dK1uF<-8JAb1=1I#uZCMJ7e0COljxH<(=BH9?1?yf p?5IlEqPt>tA3oVY9K diff --git a/doc/Eqs/pair_comb2.tex b/doc/Eqs/pair_comb2.tex index 07822d119c..faa664ef09 100644 --- a/doc/Eqs/pair_comb2.tex +++ b/doc/Eqs/pair_comb2.tex @@ -2,11 +2,21 @@ \begin{document} -\begin{eqnarray*} -U_{ij}^R(r_{ij}) & = & f_{S_{ij}} A_{ij} \exp (-\lambda_{ij} r_{ij}) \\ -U_{ij}^A(r_{ij},q_i,q_j) & = & -f_{S_{ij}} b_{ij} B_{ij} \exp (-\alpha_{ij} r_{ij}) \\ -U_{ij}^I(r_{ij},q_i,q_j) & = & J_{ij} (r_{ij}) q_i q_j \\ -U_{ij}^V(r_{ij}) & = & f_{L_{ij}} (C_{VDW_i}C_{VDW_j})^\frac{1}{2}/r_{ij}^6 \\ -\end{eqnarray*} +\begin{table}[h] +\begin{tabular}{|c|c|c|c|c|c|c|c|} +\hline + & $O$ & $Cu$ & $N$ & $C$ & $H$ & $Ti$ & $Zn$ \\ \hline +$O$ & F & F & F & F & F & F & F \\ \hline +$Cu$ & F & F & P & F & F & P & F \\ \hline +$N$ & F & P & F & M & F & P & P \\ \hline +$C$ & F & F & M & F & F & M & M \\ \hline +$H$ & F & F & F & F & F & M & M \\ \hline +$Ti$ & F & P & P & M & M & F & P \\ \hline +$Zn$ & F & F & P & M & M & P & F \\ \hline +\multicolumn{8}{l}{F: Fully optimized} \\ +\multicolumn{8}{l}{M: Only optimized for dimer molecule} \\ +\multicolumn{8}{l}{P: in Progress but have it from mixing rule} \\ +\end{tabular} +\end{table} \end{document} diff --git a/doc/Section_commands.html b/doc/Section_commands.html index 0e4209b95d..51cbc24029 100644 --- a/doc/Section_commands.html +++ b/doc/Section_commands.html @@ -482,12 +482,12 @@ potentials. Click on the style itself for a full description: lj/long/coul/longlj/long/dipole/longlj/long/tip4p/longlj/smooth lj/smooth/linearlj96/cutlubricatelubricate/poly lubricateUlubricateU/polymeammie/cut -morsenm/cutnm/cut/coul/cutnm/cut/coul/long -peri/lpsperi/pmbperi/vesreax -reboresquaredsoftsw -tabletersofftersoff/modtersoff/zbl -tip4p/cuttip4p/longtri/ljyukawa -yukawa/colloidzbl +morsenb3b/harmonicnm/cutnm/cut/coul/cut +nm/cut/coul/longperi/lpsperi/pmbperi/ves +reaxreboresquaredsoft +swtabletersofftersoff/mod +tersoff/zbltip4p/cuttip4p/longtri/lj +yukawayukawa/colloidzbl

These are pair styles contributed by users, which can be used if @@ -498,9 +498,9 @@ package. awpmd/cutcoul/dieleam/cdedip eff/cutgauss/cutlistlj/cut/dipole/sf lj/sdklj/sdk/coul/longlj/sdk/coul/msmlj/sf -meam/splinemeam/sw/splinenb3b/harmonicreax/c -sph/heatconductionsph/idealgassph/ljsph/rhosum -sph/taitwatersph/taitwater/morristersoff/table +meam/splinemeam/sw/splinereax/csph/heatconduction +sph/idealgassph/ljsph/rhosumsph/taitwater +sph/taitwater/morristersoff/table

These are accelerated pair styles, which can be used if LAMMPS is diff --git a/doc/Section_commands.txt b/doc/Section_commands.txt index 71ee898ae0..fd89ea9609 100644 --- a/doc/Section_commands.txt +++ b/doc/Section_commands.txt @@ -778,6 +778,7 @@ potentials. Click on the style itself for a full description: "meam"_pair_meam.html, "mie/cut"_pair_mie.html, "morse"_pair_morse.html, +"nb3b/harmonic"_pair_nb3b_harmonic.html, "nm/cut"_pair_nm.html, "nm/cut/coul/cut"_pair_nm.html, "nm/cut/coul/long"_pair_nm.html, @@ -818,7 +819,6 @@ package"_Section_start.html#start_3. "lj/sf"_pair_lj_sf.html, "meam/spline"_pair_meam_spline.html, "meam/sw/spline"_pair_meam_sw_spline.html, -"nb3b/harmonic"_pair_nb3b_harmonic.html, "reax/c"_pair_reax_c.html, "sph/heatconduction"_pair_sph_heatconduction.html, "sph/idealgas"_pair_sph_idealgas.html, diff --git a/doc/pair_comb.html b/doc/pair_comb.html index b140b639a7..8d46099968 100644 --- a/doc/pair_comb.html +++ b/doc/pair_comb.html @@ -13,51 +13,50 @@

pair_style comb/omp command

+

pair_style comb3 command +

Syntax:

-
pair_style comb 
+
pair_style comb
+pair_style comb3 keyword 
 

+
keyword = polar
+  polar value = polar_on or polar_off = whether or not to include atomic polarization 
+

+
+
 
Examples:
 

 
pair_style comb
 pair_coeff * * ../potentials/ffield.comb Si
 pair_coeff * * ../potentials/ffield.comb Hf Si O 
 

+
pair_style comb3 polar_off
+pair_coeff * * ../potentials/ffield.comb3 O Cu N C O 
+

 
Description:
 

-
Style comb computes a variable charge COMB (Charge-Optimized
-Many-Body) potential as described in (COMB_1) and
-(COMB_2).  The energy E of a system of atoms 
-is given by
+
Style comb computes the second-generation variable charge COMB 
+(Charge-Optimized Many-Body) potential.  Style comb3 computes 
+the third-generation COMB potential.  These COMB potentials are 
+described in (COMB1) and (COMB2).  Briefly, 
+the total energy ET of a system of atoms is given by
 

 

 

-
where ET is the total potential energy of the system,
-ESi is the self-energy term of atom i,
-Vij is the interatomic potential between the ith and
-jth atoms, rij is the distance of the atoms i and
-j, and qi and qj are charges of the atoms,
-and EBBi is the bond-bending term of atom i.
+
where Eiself is the self-energy of atom i 
+(including atomic ionization energies and electron affinities),
+Eijshort is the bond-order potential between 
+atoms i and j,
+EijCoul is the Coulomb interactions,
+Epolar is the polarization term for organic systems 
+(style comb3 only), 
+EvdW is the van der Waals energy (style comb3 only), 
+Ebarr is a charge barrier function, and 
+Ecorr are angular correction terms.
 

-
The interatomic potential energy Vij consists of four
-components: two-body short-range repulsion,
-URij, many-body short-range attraction,
-UAij, long-range Coulombic electrostatic
-interaction, UIij, and van der Waals energy,
-UVij, which are defined as:
-

-

-

-
The short-range repulsion and attraction are based on the
-Tersoff potential (see the pair_style
-tersoff command); thus for a zero-charge pure
-element system with no van der Waals interaction, the COMB potential
-reduces to Tersoff potential, typically truncated at a short cutoff,
-e.g. 3 to 4 Angstroms.  The long-range Coulombic term uses the Wolf
-summation method described in Wolf, spherically truncated at a
-longer cutoff, e.g. 12 Angstroms.
-

-
The COMB potential is a variable charge potential.  The equilibrium
+
The COMB potentials (styles comb and comb3) are variable 
+charge potentials.  The equilibrium
 charge on each atom is calculated by the electronegativity
 equalization (QEq) method.  See Rick for further details.
 This is implemented by the fix qeq/comb command,
@@ -67,14 +66,11 @@ command has options that determine how often charge equilibration is
 performed, its convergence criterion, and which atoms are included in
 the calculation.
 

-
Only a single pair_coeff command is used with the comb style which
-specifies the COMB potential file with parameters for all needed
-elements.  These are mapped to LAMMPS atom types by specifying N
-additional arguments after the potential file in the pair_coeff
-command, where N is the number of LAMMPS atom types.  The provided
-potential file ffield.comb contains all currently-available COMB
-parameterizations: for Si, Cu, Hf, Ti, O, their oxides and Zr, Zn and
-U metals.
+
Only a single pair_coeff command is used with the comb and comb3 
+styles which specifies the COMB potential file with parameters for 
+all needed elements.  These are mapped to LAMMPS atom types by 
+specifying N additional arguments after the potential file in the 
+pair_coeff command, where N is the number of LAMMPS atom types.  
 

 
For example, if your LAMMPS simulation of a Si/SiO2/
 HfO2 interface has 4 atom types, and you want the 1st and
@@ -93,80 +89,28 @@ This can be used when a comb potential is used as part of the
 hybrid pair style.  The NULL values are placeholders for atom types
 that will be used with other potentials.
 

-
The ffield.comb potential file is in the potentials directory of
-the LAMMPS distribution.  See the pair_coeff doc page
-for alternate ways to specify the path for the potential file.  Lines
-that are not blank or comments (starting with #) define parameters for
-a triplet of elements.  The 49 parameters in a single entry correspond
-to coefficients in the formula above:
+
For style comb, the provided potential file ffield.comb contains 
+all currently-available 2nd generation COMB parameterizations: 
+for Si, Cu, Hf, Ti, O, their oxides and Zr, Zn and U metals.  
+For style comb3, the potential file ffield.comb3 contains all 
+currently-available 3rd generation COMB paramterizations: 
+O, Cu, N, C, H, Ti and Zn.  
+The status of the optimization of the compounds, for example 
+Cu2O, TiN and hydrocarbons, are given in the 
+following table:
 

-
  • element 1 (the center atom in a 3-body interaction)
-
  • element 2 (the atom bonded to the center atom)
-
  • element 3 (the atom influencing the 1-2 bond in a bond-order sense)
-
  • m
-
  • c
-
  • d
-
  • h (cos_theta0 (can be a value -1 or 1))
-
  • n
-
  • beta
-
  • lambda21, lambda2 of element 1 (1/distance units)
-
  • lambda22, lambda2 of element 2 (1/distance units)
-
  • B of element 1 (energy units)
-
  • B of element 2 (energy units)
-
  • R (cutoff, distance units, 0.5*(r_outer + r_inner))
-
  • D (cutoff, distance units, R - r_inner)
-
  • lambda11, lambda1 of element 1 (1/distance units)
-
  • lambda12, lambda1 of element 2 (1/distance units)
-
  • A of element 1 (energy units) 
-
  • A of element 2 (energy units) 
-
  • K_LP_1 (energy units, 1st order Legendre polynomial coefficient)
-
  • K_LP_3 (energy units, 3rd order Legendre polynomial coefficient)
-
  • K_LP_6 (energy units, 6th order Legendre polynomial coefficient)
-
  • A123 (cos_theta, theta = equilibrium MOM or OMO bond angles)
-
  • Aconf (cos_theta, theta = equilibrium MOM or OMO bond-bending coefficient)
-
  • addrep (energy units, additional repulsion)
-
  • R_omiga_a (unit-less scaler for A)
-
  • R_omiga_b (unit-less scaler for B)
-
  • R_omiga_c (unit-less scaler for 0.5*(lambda21+lambda22))
-
  • R_omiga_d (unit-less scaler for 0.5*(lambda11+lambda12))
-
  • QL1 (charge units, lower charge limit for element 1)
-
  • QU1 (charge units, upper charge limit for element 1)
-
  • DL1 (distance units, ion radius of element 1 with charge QL1)
-
  • DU1 (distance units, ion radius of element 1 with charge QU1)
-
  • QL2 (charge units, lower charge limit for element 2)
-
  • QU2 (charge units, upper charge limit for element 2)
-
  • DL2 (distance units, ion radius of element 2 with charge QL2)
-
  • DU2 (distance units, ion radius of element 2 with charge QU2)
-
  • chi (energy units, self energy 1st power term)
-
  • dJ (energy units, self energy 2nd power term)
-
  • dK (energy units, self energy 3rd power term)
-
  • dL (energy units, self energy 4th power term)
-
  • dM (energy units, self energy 6th power term)
-
  • esm (distance units, orbital exponent)
-
  • cmn1 (self energy penalty, rho 1 of element 1)
-
  • cml1 (self energy penalty, rho 1 of element 2)
-
  • cmn2 (self energy penalty, rho 2 of element 1)
-
  • cmn2 (self energy penalty, rho 2 of element 2)
-
  • coulcut (long range Coulombic cutoff, distance units)
-
  • hfocor (coordination term) 
-

-
The parameterization of COMB potentials start with a pure element
-(e.g. Si, Cu) then extend to its oxide and polymorphs
-(e.g. SiO2, Cu2O). For interactions not
-involving oxygen (e.g.  Si-Cu or Hf-Zr), the COMB potential uses a
-mixing rule to generate these parameters.  For furthur details on the
-parameterization and parameters, see the Tersoff
-doc page and the COMB publications (COMB_1) and
-(COMB_2).  For more details on 3-body interaction types
-(e.g. SiSiO vs SiOSi), the mixing rule, and how to generate the
-potential file, please see the Tersoff doc page.
+

+

+
For style comb3, in addition to ffield.comb3, a special 
+parameter file, lib.comb3, 
+that is exclusively used for C/O/H systems, will be automatically 
+loaded if carbon atom is detected in LAMMPS input structure.  
+Keyword polar indicates whether the force field includes the atomic 
+polarization. Since the equilibration of the polarization has not 
+yet been implemented, it can only set polar_off at present. 
 

-
In the potentials directory, the file ffield.comb provides the
-LAMMPS parameters for COMB's Si, Cu, Ti, Hf and their oxides, as well
-as pure U, Zn and Zr metals. This file can be used for pure elements
-(e.g. Si, Zr), binary oxides, binary alloys (e.g.  SiCu, TiZr), and
-complex systems.  Note that alloys and complex systems require all
-3-body entries be pre-defined in the potential file.
+
IMPORTANT NOTE: You can not use potential file ffield.comb with 
+style comb3, nor file ffield.comb3 with style comb.
 

 

 
@@ -198,16 +142,16 @@ more instructions on how to use the accelerated styles effectively.
 two different element types, mixing is performed by LAMMPS as
 described above from values in the potential file.
 

-
This pair style does not support the pair_modify
+
These pair styles does not support the pair_modify
 shift, table, and tail options.
 

-
This pair style does not write its information to binary restart
+
These pair styles do not write its information to binary restart
 files, since it is stored in potential files.  Thus, you
 need to re-specify the pair_style, pair_coeff, and fix
 qeq/comb commands in an input script that reads a
 restart file.
 

-
This pair style can only be used via the pair keyword of the
+
These pair styles can only be used via the pair keyword of the
 run_style respa command.  It does not support the
 inner, middle, outer keywords.
 

@@ -215,15 +159,15 @@ restart file.
 
 
Restrictions:
 

-
This pair style is part of the MANYBODY package.  It is only enabled
+
These pair styles are part of the MANYBODY package.  It is only enabled
 if LAMMPS was built with that package (which it is by default).  See
 the Making LAMMPS section for more info.
 

-
This pair style requires the newton setting to be "on"
+
These pair styles requires the newton setting to be "on"
 for pair interactions.
 

-
The COMB potentials in the ffield.comb file provided with LAMMPS
-(see the potentials directory) are parameterized for metal
+
The COMB potentials in the ffield.comb and ffield.comb3 files provided 
+with LAMMPS (see the potentials directory) are parameterized for metal
 units.  You can use the COMB potential with any LAMMPS
 units, but you would need to create your own COMB potential file with
 coefficients listed in the appropriate units if your simulation
@@ -238,27 +182,19 @@ doesn't use "metal" units.
 

 

 
-
+
 
-
(COMB_1) J. Yu, S. B. Sinnott, S. R. Phillpot, Phys Rev B, 75, 085311 (2007),
+
(COMB1) S. R. Phillpot and S. B. Sinnott, Science 325, 1634-1635 (2009)
 

-
+
 
-
(COMB_2) T.-R. Shan, B. D. Devine, T. W. Kemper, S. B. Sinnott, S. R.
-Phillpot, Phys Rev B, 81, 125328 (2010).
-

-
-
-
(Tersoff) J. Tersoff, Phys Rev B, 37, 6991 (1988).
+
(COMB2) T. Liang, T.-R. Shan, Y.-T. Cheng, B. D. Devine, M. Noordhoek, 
+Y. Li, Z. Lu, S. R. Phillpot, and S. B. Sinnott, Mat. Sci. & Eng: R, 
+in press (DOI: 10.1016/j.mser.2013.07.001)
 

 
 
 
(Rick) S. W. Rick, S. J. Stuart, B. J. Berne, J Chem Phys 101, 6141
 (1994).
 

-
-
-
(Wolf) D. Wolf, P. Keblinski, S. R. Phillpot, J. Eggebrecht, J Chem
-Phys, 110, 8254 (1999).
-

 
diff --git a/doc/pair_comb.txt b/doc/pair_comb.txt
index 95e8ddb682..28eeacc100 100644
--- a/doc/pair_comb.txt
+++ b/doc/pair_comb.txt
@@ -8,52 +8,48 @@
 
 pair_style comb command :h3
 pair_style comb/omp command :h3
+pair_style comb3 command :h3
 
 [Syntax:]
 
-pair_style comb :pre
+pair_style comb
+pair_style comb3 keyword :pre
+
+keyword = {polar}
+  {polar} value = {polar_on} or {polar_off} = whether or not to include atomic polarization :pre
+:ule
 
 [Examples:]
 
 pair_style comb
 pair_coeff * * ../potentials/ffield.comb Si
 pair_coeff * * ../potentials/ffield.comb Hf Si O :pre
+pair_style comb3 polar_off
+pair_coeff * * ../potentials/ffield.comb3 O Cu N C O :pre
 
 [Description:]
 
-Style {comb} computes a variable charge COMB (Charge-Optimized
-Many-Body) potential as described in "(COMB_1)"_#COMB_1 and
-"(COMB_2)"_#COMB_2.  The energy E of a system of atoms 
-is given by
+Style {comb} computes the second-generation variable charge COMB 
+(Charge-Optimized Many-Body) potential.  Style {comb3} computes 
+the third-generation COMB potential.  These COMB potentials are 
+described in "(COMB1)"_#COMB1 and "(COMB2)"_#COMB2.  Briefly, 
+the total energy {ET} of a system of atoms is given by
 
 :c,image(Eqs/pair_comb1.jpg)
 
-where {ET} is the total potential energy of the system,
-{ESi} is the self-energy term of atom {i},
-{Vij} is the interatomic potential between the {i}th and
-{j}th atoms, {rij} is the distance of the atoms {i} and
-{j}, and {qi} and {qj} are charges of the atoms,
-and {EBBi} is the bond-bending term of atom {i}.
+where {Eiself} is the self-energy of atom {i} 
+(including atomic ionization energies and electron affinities),
+{Eijshort} is the bond-order potential between 
+atoms {i} and {j},
+{EijCoul} is the Coulomb interactions,
+{Epolar} is the polarization term for organic systems 
+(style {comb3} only), 
+{EvdW} is the van der Waals energy (style {comb3} only), 
+{Ebarr} is a charge barrier function, and 
+{Ecorr} are angular correction terms.
 
-The interatomic potential energy {Vij} consists of four
-components: two-body short-range repulsion,
-{URij}, many-body short-range attraction,
-{UAij}, long-range Coulombic electrostatic
-interaction, {UIij}, and van der Waals energy,
-{UVij}, which are defined as:
-
-:c,image(Eqs/pair_comb2.jpg)
-
-The short-range repulsion and attraction are based on the
-"Tersoff"_#Tersoff potential (see the "pair_style
-tersoff"_pair_tersoff.html command); thus for a zero-charge pure
-element system with no van der Waals interaction, the COMB potential
-reduces to Tersoff potential, typically truncated at a short cutoff,
-e.g. 3 to 4 Angstroms.  The long-range Coulombic term uses the Wolf
-summation method described in "Wolf"_#Wolf, spherically truncated at a
-longer cutoff, e.g. 12 Angstroms.
-
-The COMB potential is a variable charge potential.  The equilibrium
+The COMB potentials (styles {comb} and {comb3}) are variable 
+charge potentials.  The equilibrium
 charge on each atom is calculated by the electronegativity
 equalization (QEq) method.  See "Rick"_#Rick for further details.
 This is implemented by the "fix qeq/comb"_fix_qeq_comb.html command,
@@ -63,14 +59,11 @@ command has options that determine how often charge equilibration is
 performed, its convergence criterion, and which atoms are included in
 the calculation.
 
-Only a single pair_coeff command is used with the {comb} style which
-specifies the COMB potential file with parameters for all needed
-elements.  These are mapped to LAMMPS atom types by specifying N
-additional arguments after the potential file in the pair_coeff
-command, where N is the number of LAMMPS atom types.  The provided
-potential file {ffield.comb} contains all currently-available COMB
-parameterizations: for Si, Cu, Hf, Ti, O, their oxides and Zr, Zn and
-U metals.
+Only a single pair_coeff command is used with the {comb} and {comb3} 
+styles which specifies the COMB potential file with parameters for 
+all needed elements.  These are mapped to LAMMPS atom types by 
+specifying N additional arguments after the potential file in the 
+pair_coeff command, where N is the number of LAMMPS atom types.  
 
 For example, if your LAMMPS simulation of a Si/SiO2/
 HfO2 interface has 4 atom types, and you want the 1st and
@@ -89,80 +82,28 @@ This can be used when a {comb} potential is used as part of the
 {hybrid} pair style.  The NULL values are placeholders for atom types
 that will be used with other potentials.
 
-The {ffield.comb} potential file is in the {potentials} directory of
-the LAMMPS distribution.  See the "pair_coeff"_pair_coeff.html doc page
-for alternate ways to specify the path for the potential file.  Lines
-that are not blank or comments (starting with #) define parameters for
-a triplet of elements.  The 49 parameters in a single entry correspond
-to coefficients in the formula above:
+For style {comb}, the provided potential file {ffield.comb} contains 
+all currently-available 2nd generation COMB parameterizations: 
+for Si, Cu, Hf, Ti, O, their oxides and Zr, Zn and U metals.  
+For style {comb3}, the potential file {ffield.comb3} contains all 
+currently-available 3rd generation COMB paramterizations: 
+O, Cu, N, C, H, Ti and Zn.  
+The status of the optimization of the compounds, for example 
+Cu2O, TiN and hydrocarbons, are given in the 
+following table:
 
-element 1 (the center atom in a 3-body interaction)
-element 2 (the atom bonded to the center atom)
-element 3 (the atom influencing the 1-2 bond in a bond-order sense)
-m
-c
-d
-h (cos_theta0 (can be a value -1 or 1))
-n
-beta
-lambda21, lambda2 of element 1 (1/distance units)
-lambda22, lambda2 of element 2 (1/distance units)
-B of element 1 (energy units)
-B of element 2 (energy units)
-R (cutoff, distance units, 0.5*(r_outer + r_inner))
-D (cutoff, distance units, R - r_inner)
-lambda11, lambda1 of element 1 (1/distance units)
-lambda12, lambda1 of element 2 (1/distance units)
-A of element 1 (energy units) 
-A of element 2 (energy units) 
-K_LP_1 (energy units, 1st order Legendre polynomial coefficient)
-K_LP_3 (energy units, 3rd order Legendre polynomial coefficient)
-K_LP_6 (energy units, 6th order Legendre polynomial coefficient)
-A123 (cos_theta, theta = equilibrium MOM or OMO bond angles)
-Aconf (cos_theta, theta = equilibrium MOM or OMO bond-bending coefficient)
-addrep (energy units, additional repulsion)
-R_omiga_a (unit-less scaler for A)
-R_omiga_b (unit-less scaler for B)
-R_omiga_c (unit-less scaler for 0.5*(lambda21+lambda22))
-R_omiga_d (unit-less scaler for 0.5*(lambda11+lambda12))
-QL1 (charge units, lower charge limit for element 1)
-QU1 (charge units, upper charge limit for element 1)
-DL1 (distance units, ion radius of element 1 with charge QL1)
-DU1 (distance units, ion radius of element 1 with charge QU1)
-QL2 (charge units, lower charge limit for element 2)
-QU2 (charge units, upper charge limit for element 2)
-DL2 (distance units, ion radius of element 2 with charge QL2)
-DU2 (distance units, ion radius of element 2 with charge QU2)
-chi (energy units, self energy 1st power term)
-dJ (energy units, self energy 2nd power term)
-dK (energy units, self energy 3rd power term)
-dL (energy units, self energy 4th power term)
-dM (energy units, self energy 6th power term)
-esm (distance units, orbital exponent)
-cmn1 (self energy penalty, rho 1 of element 1)
-cml1 (self energy penalty, rho 1 of element 2)
-cmn2 (self energy penalty, rho 2 of element 1)
-cmn2 (self energy penalty, rho 2 of element 2)
-coulcut (long range Coulombic cutoff, distance units)
-hfocor (coordination term) :ul
+:c,image(Eqs/pair_comb2.jpg)
 
-The parameterization of COMB potentials start with a pure element
-(e.g. Si, Cu) then extend to its oxide and polymorphs
-(e.g. SiO2, Cu2O). For interactions not
-involving oxygen (e.g.  Si-Cu or Hf-Zr), the COMB potential uses a
-mixing rule to generate these parameters.  For furthur details on the
-parameterization and parameters, see the "Tersoff"_pair_tersoff.html
-doc page and the COMB publications "(COMB_1)"_#COMB_1 and
-"(COMB_2)"_#COMB_2.  For more details on 3-body interaction types
-(e.g. SiSiO vs SiOSi), the mixing rule, and how to generate the
-potential file, please see the "Tersoff"_pair_tersoff.html doc page.
+For style {comb3}, in addition to ffield.comb3, a special 
+parameter file, {lib.comb3}, 
+that is exclusively used for C/O/H systems, will be automatically 
+loaded if carbon atom is detected in LAMMPS input structure.  
+Keyword {polar} indicates whether the force field includes the atomic 
+polarization. Since the equilibration of the polarization has not 
+yet been implemented, it can only set polar_off at present. 
 
-In the potentials directory, the file {ffield.comb} provides the
-LAMMPS parameters for COMB's Si, Cu, Ti, Hf and their oxides, as well
-as pure U, Zn and Zr metals. This file can be used for pure elements
-(e.g. Si, Zr), binary oxides, binary alloys (e.g.  SiCu, TiZr), and
-complex systems.  Note that alloys and complex systems require all
-3-body entries be pre-defined in the potential file.
+IMPORTANT NOTE: You can not use potential file {ffield.comb} with 
+style {comb3}, nor file {ffield.comb3} with style {comb}.
 
 :line
 
@@ -194,16 +135,16 @@ For atom type pairs I,J and I != J, where types I and J correspond to
 two different element types, mixing is performed by LAMMPS as
 described above from values in the potential file.
 
-This pair style does not support the "pair_modify"_pair_modify.html
+These pair styles does not support the "pair_modify"_pair_modify.html
 shift, table, and tail options.
 
-This pair style does not write its information to "binary restart
+These pair styles do not write its information to "binary restart
 files"_restart.html, since it is stored in potential files.  Thus, you
 need to re-specify the pair_style, pair_coeff, and "fix
 qeq/comb"_fix_qeq_comb.html commands in an input script that reads a
 restart file.
 
-This pair style can only be used via the {pair} keyword of the
+These pair styles can only be used via the {pair} keyword of the
 "run_style respa"_run_style.html command.  It does not support the
 {inner}, {middle}, {outer} keywords.
 
@@ -211,15 +152,15 @@ This pair style can only be used via the {pair} keyword of the
 
 [Restrictions:]
 
-This pair style is part of the MANYBODY package.  It is only enabled
+These pair styles are part of the MANYBODY package.  It is only enabled
 if LAMMPS was built with that package (which it is by default).  See
 the "Making LAMMPS"_Section_start.html#start_3 section for more info.
 
-This pair style requires the "newton"_newton.html setting to be "on"
+These pair styles requires the "newton"_newton.html setting to be "on"
 for pair interactions.
 
-The COMB potentials in the {ffield.comb} file provided with LAMMPS
-(see the potentials directory) are parameterized for metal
+The COMB potentials in the {ffield.comb} and {ffield.comb3} files provided 
+with LAMMPS (see the potentials directory) are parameterized for metal
 "units"_units.html.  You can use the COMB potential with any LAMMPS
 units, but you would need to create your own COMB potential file with
 coefficients listed in the appropriate units if your simulation
@@ -234,20 +175,14 @@ doesn't use "metal" units.
 
 :line
 
-:link(COMB_1)
-[(COMB_1)] J. Yu, S. B. Sinnott, S. R. Phillpot, Phys Rev B, 75, 085311 (2007),
+:link(COMB1)
+[(COMB1)] S. R. Phillpot and S. B. Sinnott, Science 325, 1634-1635 (2009)
 
-:link(COMB_2)
-[(COMB_2)] T.-R. Shan, B. D. Devine, T. W. Kemper, S. B. Sinnott, S. R.
-Phillpot, Phys Rev B, 81, 125328 (2010).
-
-:link(Tersoff)
-[(Tersoff)] J. Tersoff, Phys Rev B, 37, 6991 (1988).
+:link(COMB2)
+[(COMB2)] T. Liang, T.-R. Shan, Y.-T. Cheng, B. D. Devine, M. Noordhoek, 
+Y. Li, Z. Lu, S. R. Phillpot, and S. B. Sinnott, Mat. Sci. & Eng: R, 
+in press (DOI: 10.1016/j.mser.2013.07.001)
 
 :link(Rick)
 [(Rick)] S. W. Rick, S. J. Stuart, B. J. Berne, J Chem Phys 101, 6141
 (1994).
-
-:link(Wolf)
-[(Wolf)] D. Wolf, P. Keblinski, S. R. Phillpot, J. Eggebrecht, J Chem
-Phys, 110, 8254 (1999).