From 80d4bae49fbaf5616f688a773933b0d3553796a3 Mon Sep 17 00:00:00 2001 From: Wojciech Kozlowski Date: Wed, 3 Aug 2016 20:35:53 +0100 Subject: [PATCH] Finished chapter 5 --- Chapter4/chapter4.tex | 1 + Chapter5/Figs/comp.pdf | Bin 0 -> 39788 bytes Chapter5/Figs/figure3.pdf | Bin 0 -> 52511 bytes Chapter5/chapter5.tex | 720 +++++++++++++++++++++++++++++++++++++- Preamble/preamble.tex | 5 +- References/references.bib | 62 ++++ thesis.tex | 2 +- 7 files changed, 781 insertions(+), 9 deletions(-) create mode 100644 Chapter5/Figs/comp.pdf create mode 100644 Chapter5/Figs/figure3.pdf diff --git a/Chapter4/chapter4.tex b/Chapter4/chapter4.tex index 6e5cdf7..f0fbc60 100644 --- a/Chapter4/chapter4.tex +++ b/Chapter4/chapter4.tex @@ -280,6 +280,7 @@ correlations itself. \mynote{introduce citations from PRX/PRA above} \section{The Master Equation} +\label{sec:master} A quantum trajectory is stochastic in nature, it depends on the exact timings of the quantum jumps which are determined randomly. This makes diff --git a/Chapter5/Figs/comp.pdf b/Chapter5/Figs/comp.pdf new file mode 100644 index 0000000000000000000000000000000000000000..897475fde33ac4b68c7bd6dd7c18d5d50f5930c7 GIT binary patch literal 39788 zcmZs?1#lcek2Y*Oju~T)nVFfHnHkpC-qZ|EKmg9!#mUUr7S1F4!biTAzTNSl$9Gt}hoG%CvIGx-2Q!(FcXTw(0Su1~KQI9l7EmM}HVHt{| zhQB`Af|oD7XCB^v_Z55n-v1abk6wKOv#>X2^b3Xq|$;|v!oJNOA`!d(tDAgCu_T{$OxQh)L4DgZ} zZeVv@!Yz-|SvGLg5&eT3G$HKDW1%3lOiq0@A%Hf>38d>O&X$&SETe~9@#`vAK3l$;zXRjfm%GA_BQ?z24@&4yuPoSJLrZLXw9ohiN-_Ov=&SiZ zvg!3SwPOsdemqb&<54#^ftis*Tj5s=l_Zijr!MCw^T(@9Y8Keg#_PzQDODUwRhOic zFF6%2-8NHMe^s&m!7g+9rJ}isAf@G>G}(Wo8DvU8&X%lHsW_0TI6PC4A^0ao4TLLR zURYPN;#RSG`IFoX0TL5P(TQEs@x=hfj67KZE>29)G?cnZ1HvE))7Q{Rws0M z6F7pnJ3bs8{kvM`2YYq9mf~nwjXHxu zbXVg`2c<6`Z+FXFHPhwT#1j31UZ?-7;X))f3nJEy5vn=PmE2`ZKAbV`_(a5w-G`to zz7qa42Jyke4)2ZPAW=R}!SqebY*F-3KWO!s=F6My-H5^3o5Zx@f!5lJ1crhs+aC*E zZHukj92f%G=p82WC*-|NShrr?)4v5+fQ+5 zi(Vo*cSc)T!08ZoyiOE9rJs{!)7qk@y*^tQE4FhpJkGvno?CW{w*rJv9FKQFuiPAs z>>3-sdZmxB1)93SqZKkwYY5t&zu9htSl{$0>yY ze*r^6kNOluRPj$L*B00VPFZ$1GDdPm6!L3T)0`3R>~{4|yQhKfG8>L}q8C02s<;+| zIgs>E6=u4wE2p_xv2Pu!*ty<+Dm zlv{B<4Uea}vCy;f5z(CO7ducaXWpkZ?(zSwHvfc;)I_^jrG!3o0C{VA}ZW z;ihF44S;%=j&ChMhRzOl|11Kym+@@ssvIp^o9jttc0OR8c-fl1!~Lc?I_58YQGRnQ zEtk43N97#O?rK;e%Ll;BJohQ4uAHI>KmeEHONaa*;JX!mXZ1V>Vp;w9C8QArKhTyH zUx?R+B3m`mtMlB`ABQ=?Gi|-mo1~22Sz+VLFM#zqcNsbDJw8V2-FLTMusssaBfgGm zd%pg;7d??uituRo2^^0-qjaWldd10fjqgfGxiNLG>2BBn?7lbiso^OODIpC#_L}J> zNYh42qWHDjTPS0G(t*h-DU5WtAvN#?`6KSZYcc}_tL2#(O2n@Nki>f=JPTaR4^xzZ zU%e~l>_S-&c!JN5iGva4bMTMhS*ls($3jO?ntl-JMeoAx-n>SNS%`vt2!p|Qli#JV z*gWj(nF-{_#@+ok{JEw+$PEYm+lLfP)J(3D4tR|Par(8uqdki)hX~^LR|uloV1jsX zigqs{wDxDCRA*Kf-Z1+wM(+N^7@ihKBk5XLc|1-43J0&RT+CGN?uz*Uuu+Y746v0; z1EkkH+Q;>%#z2-d_ZJ$pzVHQg@H$>5pOj$~0kOkha{i4YABnJIn}wQzw-Q5)Mn_wB`n17Rl{w#O;;J357Z_9_l>!r^Go`yJeL{Ic0>UBTdkC&KmuN z_4nLOyc~x#{pv}H*5>{s-Vrjg{5uB&!53nIh(GD;n489|Da7seVY=K-Y2snp)S-6? zEnV$F1OJH1Mq*?;_ zG02UWd|s9nV2_=wqhuXj`&Gc2YwcyvJPtV1fWELb+T0X>vtfEva%nVQ!i(5xi`pcr zXWvTx%t`3f=H>*B7-MnupgIL-F?-whqBde@rwLa_PRK^@#<#gmlqn3DZgfFNjfPHc zJXBmJnQJ!lA2rcX2n_p@tvP}T_^q5oNaeG@CCs4+JOJmx5v$EQtfGwvXBPKOUD6BI zmVQvAT!I)>UJeM|G$jYN*i`Z5w#b92{S1LwAu|8)Ul`;JvGcGq3+O z%NCK4U4)18!&mkpZLtBj)H^}a+`6^9U`<5lC=Nb?-DSR*;(=ap((z(2fXLJ2x+7gQ|sF?N(@Z%QJX? zp&jD9B$UhLxnZRE^)sduCviOHn9)<)vAy)rj{Lb%ss5BTn16LNxGvfNAlR-EQL3m6 z`Jl<-J^Q;0KIPIB?fUUgi0xsNJ90O=d1yx96l3J(g4hvD0}M+VUOau7P#1BmZ(bB1 zMJO>PE>ArewgYgpfm8D#+rh2%{4PKU(e*3P-{#j~7!(TDNyo+CuLqmkm8Vti@D|u6 z9KbgMoi1*T?T6|1u9;yI;$SG4`h{VBxJz&v#Co=@2DBf^;l8oNs)Uv>{gpws*4E+d zrDz46y?`VSHc@N=2$3S&qh7|~@-x_s$QIR+?E}qoT)=1OYE$GUZSFkg$ZGvcHwQ5X z6PE_9S8Gq%wu@DM(@xC?gMW}Zg95ma3%Q9ZA8LjJ$>>j5?vTDtrXeT+!k^dEm?TJR zrD);&?{0JfgJ50@yXPV!3sdRv8T{VHu`7*vq2x^@uFau))33TUGSe|)u(`Qt1txk_sniK=0NHbObBj-`eTN8 zRhoSy?>t@2Y%QDxcMg{8Y-LdN?iB?uHs^VaP)N&c+by^f9GX2u2Lsz`6gfygYySt-A{;9GSHdUTC|9A-lfs;U<=#of?n5iasf# z$>7I1;#qE}TnOce@a%Mb$tu90!P!?XXX74hLjpc%y6~{pkFF0~=xb?SvM46l6{VK6;ISkBi8p?@)~&6&xKGvTux%M4-tO;%EVP zV1N3kMU4pW31xq4!;188ql!4JO|jwxDm07KRx#n0f*lD8oW;T0!DFAlYTtCw7k%Tm z<`6JR=hYT`uuijVXC|aXM|rTaX%zf&9k5(w-4JF{usF+MOFDaqx^nQ8MRlP6j6?_k*hD4OI+&Mc0(UC}T4Bl&{B6%$2J;RJ;4u3n?#C~+ zlN|#c*60UOIFWN}N-O48fO~Nsfm-c1CU_D|TAL>{18T|H!Z%+j;S;R|a_UE$__v9{ z+#$ep@xX{CoLoHv8)X+fP@#6gKV91E{W?1AVzN}+z{Qm{JJ!@7Kd~z6qTeWk$s;mp{;4!LfG6J-i12+%KM7?!k#nPt zle1ub=TUu>v!3~!U^ANkHydwUBVZ7sk#Jeb2HMC%S7;#S3U`)kqy93>{$}smk2J=T zdh|0A=HRB+{^s^xaAS$`Q@{Kw1w|)UabGuvb*O`CZExpPZN~Q zUwof}Q~T*sc%ONC^&tAv-k;IWqcx>$5ZyK-=FX3}s%Y>pZeLLY} z9uMcl;9s%#y3d1;#Ilb-f`_`#rJuQBA7v(lo~!mSxO+c(h3~QV>pl@g0tp{A$#@#| z9nJjgE}41QZxwqfBM!=KCd8yjI8>y+MA)Ide6$R=ooPS0@W5z=@dRc*eb=(s%yrc? z4s3isAY=VjxUy;!q7Dy#{PwXRjio$0Y-#e=o=(5&hds84nz&OZ9!VkOaOV|?!1hcr zn+$>8)Sd-Pt+>xY$BhEv77d0n*kri_%`&d!*$m%#*l)!V(XhkcEHcs09=1|!;^Vs`A z#rOj|%&>P%U(kw>S2Tur*mo_HKt&>B--D+CTXc4>d>)cyq1t6(n%cJ%dvnO0_(H3Z zksATREdYu}ER%}9t)>Bf@hDLxxP$yKOL)pDFB)=EKbqnR(1?3blWotr?Fa?iQYQHy*;+A(?19lbSz^WxHgNiG#k`dUQJ zd2_tUHPMn;Gzs*L3#5oEFW4+pW^C5dmkOKNWaIU_2BsR*Y_lT45hj?`F?J>51u+iP zj!1>_De_?$GubIl=k!g149j)IDrtbC9pBT(tJDDrJKpi>&YoTC*Wewj@A5dwgQxR+ z)2`WB4Z|K;9wkORX%47LLl@g+S3%HaX^z7Iq={8^LHw7RK~O>5+B_%4gJCcTV`1US zM?|?_nG26pGP$I3%orQx{N7JRY;X;A5`vAK{coy>>5%xGla)hs7j1-pP7sJBK|d{3 zSS?6z$ADM*aG+GRAOhc>Mh-vDNE33LEl9^}IjgU5_Y0{2Q6Fs7p*N1GfYFaG=J0h? zW5P~E$TN^B(yeno%a%tnze3i~)6G`2z4{x$6{*q&GhP@hlDQMhNaLfzD%ZsJbI^0D zeu9iu4uZMeAkb}t8p5CF-LFQq7?1_NlCif@?hQB zJ#aWSY@kTs%noT0qtv_B1K}l_y}@UTD2S^j(5)&^lhy9t!fm75gHd|hb&?{ij_6CFiwosHw>m)E_Eu7v`-C|wUJ~et>>3t=%Pb#aV;1 zmKTI#w132epFq$6`{UQWus#E^_=%24kM`iyG-HNfSXD@cY&7wYZOocAwUN{Q_H+^IxyHtc4ysBjW>F%!P&s|KO zgSIg+815pQFP=Pgt4U+aZ1&*FG4!26Tvo>zi6$iZH6h6W%bUjllsZ{9HR1(vj@|VM zh4NYQAs8Ro3C*DN-7QRXI+A1J3iT8`y;Ap?Fey# zLKj`amiDh=Flc5~bh@M#O5%FLuhtqThkHk(Yis>97gE8$Ma+o{u(*RQAcCfTxBVSb zSn#SF+c)p1|I1u`JdN70yUZo(MoX`u8FaTjUOqiplDp7`e$t)gi-mP{V+WD`r%~Ie z@`LE@z^Luvj5I!pxyg20T+Mn!PRNse-Keqd3LUCT=dMY?3kb-eZHllf)X6X}C z?~Kfr&=vEYcR7xOd{bauE>aw6GmwJeqm(5pg9;U&Y!Age>BI7kF09=aaxv5kWKn=! zPao7}x-B)wMkel~52yHW#_3{wsAMKS)l??#sy|myNfU6KOL;7o>L>5r@xJK#-Y_rq zM9CL9i|jBN2*}&kTKPeaW@Et5HW-?{#kfX3dyx*7vHEtOVe2>7o zq1fZ;+hQBwHAzCWB2%f4vT&yKu2KV<9gw1{FK5!xx0c-wnv&uWkv;jQh$ioYTVJS# zj#f<9c2NSZEpR!|%(EKL1e)HuQAuL3`nuoA8n5k2B;uy80$SgWoCb3APD>o`zB`3G z;yV4zAob*Fh!^y`EQy%ryE!igDx7N{Fcq3auA!LfO-Lcs(Rx>j6-V3jH>c>Q3#5^R z&~B7?vqxlHCS{afxIyCA(=T1bxu4+gUc@!nlDWPG30nrwN!9)Hpfg zu%-z7!bQ98{M|k+yK%KGWBX-dEcsBY-9oezhG|;A=p~_a@$^T~Nq7<_#0s_GwYPGG0ez2%Xv=HFGw9d5q*g`u?0925u3>s~({! zDt054dQL=;=id<3!(9^O+gcjXemP5Xi)<#|O%s{7D@p_&p1;T0OPZwQq{0W=^PK;z zsE8{0vyfOc&#!Yi1W3%>9S{V2imZ!%x^ zhdo2O?|`x)N%q`z0mEQ*F|w37=V>T|NVt-bIA%x(r!Xmy-%UA;>GIEu@2b-+CIX() zk0Ft{`)?%uCnX0OV)R8ZNc2tNYpTzKUiDO>mD+~9StHg!%h!k(uVMY$ggn*lDc3@b zDHG?A>9{T5l!WJR6I1xGdcoG5H9ChvpPjkgVpn&)d`=fF??j=Flx zppZdRQ9m6@ZZF0wumx>T$gCTbxcv?*I||06za73a#G8wLRsYQeJ-49OCv4t~4`(3H zYxA9$W$1Xr{Bk~kYDocT((9)eJ>>QB15X>cAc;s}qzoT*V|5Wye{QysiIn)Is#QQ` zhRT^tQL?xsf&QKrIEuT8X^wz;t@4b{AvSa)py7&S0=tv3A5Bvlz7VXY&Xd5xYaVE= zv)7&R-^+$*{&wu0gJft$k_gvb!UCAQUZBl}DW(U^vhnT{q+EE;I|ngG%1t6b#64snU>saeea>DYk#Uc%HuRvjc1 zKkunG_d|xs{a|)33L4{Q&@2imu@PWJD?|!wE!ditw?H5*nN*I5ds_|y=)-|#)5?05 z*uB5qZ;1%mdg-RS2C?H2V5Cs$+}X@e%l7PH~@y&|cV770Spi#?#H zj(+3V(*b^&qhW6*sQ`=>vRG4ti*Y}3CdX&V(pws^h>JkPg})-pcz45;tT#N=M?5*N zIz{D-rzpwRf4oSY$y7L(kmFBkrXn3VlZGcE!!$HPdb?TiJ{9mSo@C--NsUb~kaDS? z5kXwGOH(wxf6MlZ*Cy{amc~L%ZJd~$*GO*b*5+0^=2;gf3u-*!Ae*!M_7MBg%CP4k z#izxeTZ(%1VNdWaA!~l$-m#LKN&2-Z-{^%exv-B$NiZIFt#Rng6L>?53FLgCmJrUq zodomHmk~iK(g zrTi(L3c-of_D9J?5K!lETm|s+ocfZXne2>!<7^I-ZMu#pgnzzH)x#w)$e=A~oDSW$qjT%+ zG8|V*gZ?i}K>}m;MW0S)_zIBpj&c#ZQ6Y1T6NI%5Htc3s0*Dd8H%Ejxs1Zh1XcI?ycs^hPTlvr=aY_g%wU_gjS;f~ zlgh0rp-|g?A^PB!(T2zb9Ws6-e8K#=Yrf^4k(>6Fm|}rNWWWl{c?o*N8mELqqJWLR ze&Ii_6h_QxxRCfT5KGd$FI9f8r3g3`$ojELvDb?}rdE zgW(={Ee z##Dx#3h&ZNJxICw5zbyYh5eBHI0iTzyd)QU-6Xsll$DypI~1d{Be0Lw<2VK~XC@e4 zS9`>rHU`D#vydI!(6V}4*vu%DFAdq8{ZcKh?QUrZ_Nun-Ff)Z2lB7{^_xgKQ)D-Kf zTT9JN9a2r4&6)>7`4*2{74taSl7h3Bw{IH&yed<5lAK-sJTfWC!gOT9%4S8*fM2 zNo}1p6koV~DPLAc+&&Ve2yQTL9K>cR3PUUk+%~H&=C4SNE=hRALrrz|>o;GQX$UtLlXiu{IKm%%z z7;Q79G@rD-)uF8ro&3tC+b<`MWW<$9bvzx?9$*Ny9A+rN;*K6**&0tBJ8Rz`Bk=Cr zVIjhX`PCE^-jwHqp2-{_S%3o%n_16^`YTc((EQvaxP=^)q#AgP7#B^j3bX}d7BUkQ zNTpwi?FUr>^Jhq~QYzRyM)GfVt8po}MbO%yb8;)4c5x&TnZ4IUc5y5yBmcCcE6}=8 z7j>Q4(0Y#SYy^l9%2VO`-p@~wGm7m4*c&=z+s~$Gp+t6+Jw`+VEvhg@H0Q4P zleU~gHg%aJD!10PlwEZUhNej8*eo6-=@jiq(ysD-nh8Z(pJj?r>4S_rit(R^G98WFiMxdN8zPQ<(+2=N=EX{zPH* zT0MR$jyfp8w;NHX@bmSPLDmc}a|P|-|cdy4gw#zCbE(RQvL5U4k>B6%euY}Pkn z{BIAhM=AtboI8OEGRIlM|2$kg!5SLdf>l`G14PAJ>t_WT&6D_%wy*ZH5y-J&iRpD^ z^FAFZuxr#i#1uhP^mO@3vp78Nx)`U4xR_w6h{p_^j=;0}?Ht{o-LOY^K1m=eg=XlG z?9flW*h_hQkF&WUuG4B8fef@dPwPx9^!(un3odf=?Rpx8QEi_|;qZHukw)8yR0vpN zyK_u|{s#6;ZnKS%8zJOMI$L}GHd4eDq<#wgZ!W9wJ=#Q_cz;>KLf9~L6U?w=gCj1Q zMUU${`1zr6)vXgqd!(d-D-)^+yKqmW&>#RABteK39a;Iv{J)}@V~yk3E}T2Jq5SOt zI})pU&5cp)$N)Qm+#9yuen02mAfc%9cwXw(8FqpeU0AqFY}4Y%qDqs3?_OUmUT$B> zJ^I!{SVQz-l%w<|7@bf;e(+)EPn^AU=b!f*q9WtOD>mX&lWxN*mO6Hk5)7h$liRoK z$H4G`TwEQ7(FYZqth<(a`zO^PJsRPmg_a9`^e zc9vwf@y7>}hk_RKDrhm=p34XK6ebzrcH$NME9E3(wK!t`YKXBV{*z+M=f)zIL(ZGN zOXkrO4wJxtZUbPp3yR=<-qk3>t<-uRQ5H5ew5NXVTML63B!+|?CN7G^3r36m(P|6{ z=*|5!hpsd*mBZ2SN??7uBSi4B3=<FGMw8(o+r0$|@rFL$N^_0qCP?ICK)dUa5Ji<#IuYYZBBJ`NiGYrb{g|nE z!c3{z+FKZ@QUlC!O${M-wJ>5>)6zA>!!(z5vh6Hq%1A^0k_EW%upN$U|4K92c5=)r z+Y|?#!UkiYB`w3Dk+~wmK`>>8WeT9vd&;?A>sNeHYq;5Q9|{bJS`2c76xSRU;fRGJ zd1(j`;h0g#ram{}YF{0WZpdS%So2^qNsxD8)pi>Qmn6pBr*_d3pGhgGz>Y{E46&gj z8yndYv7~*4X{e7vg-c7fr|=}k4a=YbiFw)+<02D$g<8`hT!DeQ*s(C80eU#pVw36l zxKdhn&|vjAE}<;CTgrYn-C2;GRRvEHMc9vDoa{9DG0r8C_-k1f4T0aBQoge^If&3& zj$3o&k26)}fD*SBNk3OCg`VMB1Z%8*tWx`dH|qraJ{N*I&M40yZ8Tt#GZ#<$*;g{r zvfy_Nna(jTGDRIcL5;#ipgkV6m+(={capb;kO9sagKT;_=4=Bx;-WESwHD?hR!iVK z*y5T6)^wSqRR2{9*XeAem8js>E#3ISBr;R7&x*!HI93;0x90J8?U@S^7?4#q=Yg$Cz?(1vKHfOi(m0ow&(@t2)CJ z=n6Y>HK&;sh_3hZ26C1mEPVu22X(jYO7ui1d!?gy@w(x7vz+*mV!||R7Wso(9sWj( zTLUymp~MwM@b(lGDg=6fu9735Dw4E;$Eb1X0Ix>P$pJoqm+$F^a9YSnzG?lB0u`?a zKF`shAMQc1ywCHwJ7>T2JVup8qYYb>JVw`|l0x+nj&c6w;aoIV*x%b_;14p*38w|5 zRKiDt3=L(OFW}u}amIg<$#H9t^z+A36jz@Em_m#rnfV?q5YDlz)ovZJ%M;^L2fT<0fjKT1x1c?jm=GTJ_CI^BLZCMYYYS*c;5Qf>H7+6} zq%H7^kE7R1*x2tkg;XfR7sht&peIl)V|E!HRewD(I{&5D1_GLG$gx|sb6A3f8U{l8 zU3P-X31v3Wrc1HmD@$+S&uEKTiI#>a;Y-cP_I2TOs0K(Skknx;Vh;Q6dZd6zj{c9d z5*u?l4|Z{4!Z1gg7p;g>_!kc*cI)^!=qGOBvSC|??!>q-p9xBs_SZ~#B{y}315jv} zYuS~?MoW;H!}Q@OaXXX$`S%6e`I+GW@YeTRgXhOYPpSzO`_J(9TRlO^1{ztd*c7Sw zxFXmU4Y$36{diJAS+RdRfrtPSE5>cT09O3LzA=tLpw~GKS^`3$JQbtQV(7OwlJE8Xb$4so!2 za8Fv{&T-veSJYil4*+QcuJpu)w#wt9#YD83Ng&;|_R6wFb<%tH@1Z#fJ{%lymXfLb_F2C(OYbd;4dPY8v@jeBIxR=WW5}l)oAPDMaa`a>B zcBj^lrgQ|_GlKo3M{ZJW02B%oiyrW)j+zWO_5R)biNRP`3V0@I!%=CKec+!e{L~{n z;Bat_Gk%GtRgCr#rd_O2T}jdI!>oaT^WG>ud`(r15=W~EFL!>K7EOLZ!bw7bv!1$W z`qDOvTt?J=vkO~M&=sE7C5?MWzbky8C$l$E&5c3O>suA$X;K6o8Va7 z6~cNwnV@z$(8xTS;+*YJYK?gG{16PR4bOfQ&o9Q=uQK8n<;Zz$JEzkqtk*+Hf5tiAYEyfOmfc4<5I;UZ>$Eh`%lSqFW+=V=ZlW4a?o8f@ZNV#-XYB+0>D-~E zR8ulovP1(_?J*8lj4}FqNv4Iqy283b&8fbrc*F16CPaj)=b8Nc0BwIC{|v`_@2sva zmerKex2U5o;^ro`=N~sIJc;TeYIgjaj~d>k+($( ze=}%+kRI@DM^CEm+*|6bR3E5sE4~I{j}$>3k(8b-stdZB4kT4T>XBx2>^bvEzYf8T zo=w^t!7#*REXejku~L^aad=k{vD34S-YjXEjn8^w#H^D&7MJ;y^UuT3*}MH+{xJj7 zV<$e{{dv~F-zyuIpLml^%^x*MdfKNPU}My)w?2NP@AW6QLtn3p^H7t*VYLV);A)bq zZv0D+{esIdKtCoCPi!THcueMt0Zu4F&AHwOGlIq9FRk!C<$|=R99TUblr+t6@qcJf z)j;{%7pSxI)pr?Iu4Ok1@Y%heys*zloDTX_-m0evaCI`mUP%xkN(zYtD4NZ_0{W<7)AjZVnZ`p`tAQQmr zFcgR!o*aLVfg}6-$Lf?C7g@aFsn9Fp$vg^c*?zUG!B^gkvC7a};*6@w(5HnCrrK6L znnAnQL3TMcSxau0U%QtEbJ<>v=rDB!B=ERfrJp9kW{3s4Um$M$Q4{@o_|SP_QWW-k zft~Ek{SfzvxPY2gyO$1gnK8iVD{Dx&o%Hq#)zugJkS1oD_OGrmEFa^?k*yZAEx z?l%!*_BNBs;XB9My#;7APlK7U2e?J#H|!z!lTGiu&tPD8E*K5AOumDaP5#o7>NV(x zZMDfizyc+8&1WXzM`tQO69%cOwnt_aj~_NU^43IIEqd#GfmWDu#G0W(CfL>vX+GJN zW78%d67V&@PDnwhB`kGKde-(qs5yRBZJTc=F%)leEAcT0-A5spUXh#G-92R_-2-22 z$9MW?SiVt1scIaGyn>P}ASlGIpb*!CLOc!%@$Jh^sMo-{n~;7SMRlO z{LcEY&8C8>36cC2*QgXMvc5R7@Utti$mX}(-vo@_hMdYJZze~*bw0S;g)l9Q z0^Mz?Rqclp{F=izqYscu`i<^cc2}w@Q`TZ>PwwG+x{yR!og{o{7iHt(~zi zKbLP?xo0jq6s7H)%ETHy%FNI?lck(S)drYgdZ9g^Eav!}>~tvk)FZTtJPcFjNf^KC zyeNnFKn5P;uE=VZKp(%TgwC~ydAE+j&vV+SMf7lSUZxxaJz+5UwtKp7y$hT=;S-Oe zf)fDAa-CtD(g~4IbJeNiKf41d@Ay3`&3d*^hNL?fLhpX+Pf>BsEEd1a?3z$W!pG#P z#R{xHDBRhao(HEZc`2sWVi1NL$P#h&__kZq@WRB_PNGbtM%ONT6g+nX&(kKNbK^*0dmiWKDoxj#UAY#otx9IZx?eZ%o^NmEDD^t7P5T&vSG$Od z(;m8{o$M`kICZt!Lp=z!j4}$>hArQ*y`&;+FD*khWOtOuxHc%) z3l!-M15l(N-Fg7p<=r_I8?2wz;uWmImlbi4Wnxdb_xES&<|M&s zM%dE=FJ#8A6>AR~L+b(*ug0`7XgY3^5T@a1lF^$)WsgR1Gx+HljK( z5#UM7KxwW?yiT%WGfacBDzid^+T956;h&Rwy=$`OUrOICQMPD%eVd0_|Lz)K-%e6O8bk(3@ZWSqh#8bU{8*e z9}Z#s!6dSmhl{)(-F9J8jb6Z1?z`))brFkFkj@P6>X3o6X#*0Rx}os*{%@=kVTxB% zlwC;xT={s2#GGeMYn52@wHH-|Mp5ENw@)j1B#qL@;(@K;r4%2odd{Kv9`eIVeuD+L z&Ir&z*pgTABjg1m!syZD@rt~f$qt!R0J;%N+gI24QO9>OG4|(ke z>x9?|LN-{6ln%AX+oR`Ej<(C$N(?O=za>$l&1>KhrDzR`xugkMYN=^JzO0&Ci6lb< zG@fdAF~$fn z(}KGM%a*Po;*@u{!QkS{mNz#o>tA!mTr@2cvG2waIESM4w`olS?c?J5dg1M<`@&UIoewRM>euHu9lT3RMjpo`ZqHhcYb!Ws+%E`l|D2pN4JsA#*PM$dtDw85@Bj$Z2#{Ms+fs!+91ZI zQ%~jrdwpG21qG%6{myl>d291hO@h|5%B7y!WcS0(sT5sw7BpFk6kQjMXALE$c{>HB z0qHO!+UpjAg2l}C30VV#Y-FO;YH>76*iE8Vx?$4{MLIPPXXjCn#ZzJMpK0 zPN(dP%22tf4+iFC&7n&17{`Mp>7B_K$3tEmmuyU0*Wz<* zWCqj6LwLC~5_H;}%_f%yp!lq;b_tpcLV#b6>R=Gn;f3d!pJ_?Ji8G9 z(XYd5Cu@B{&rW^3TE-FjwtBy8Z+~1%lEjQ{*PN~YjWxQn67<#nEd7{vul70J{x||X zOvZMzpLW%=-D2AkY*_WV`O2BA+&OE1zb*KX>J-;royF5};dPjIz_8+~_q6Soj$D`H zd0kh*_i`DJa_~LJZcGvFwB_b+G+P3Rug+eB|M}_j)-&7d5-SzvT?JDbO?Xa8fPGU8 z1{*}(UxQ*~R|R^N`*i+t0gGAG#rkjvdeA1*i*8D z+FS*I-5*80uv)kpqu0MHrHg4C7OGLd%t#&_eE=l8jqdpKjf(uAdty4rmlm(GQe`jN zd-e*&G2@&u)*e8yZ-I@)3qTfA+PP7~Y;o|vw`4r&Dz$PX#Gh;I6+q?v_J_CLE(%Bf*S}mnQ>RV7&c?5g z-jN}vNhp4Y1NVr#em{&GVD6Q%O7|PMRWVCv!$ZVkJD9$zVU=3G*BYifyXtuR#g#4VPej&%(o5UdUTuwqcu9s$i@q z=QRBc`XT(l$#!s@{R+YMug$^>ZvQJ=I%j*P88=P`u%l3~Uf@ky*;Aa~n>%nCwv1bk z@}=GPg;RAoJJ@6JZRGF!+q<*R4rES_!51VXq;rch!}VD>3%=qXtOA+7@3)-W`lE*1 zQ9eR0ygH2rud`2qL725qd6$q^jQ915$b{WCf7a~ScVIpQCMu0JB1oa(*e)zAEQNrHN<$rOT z(=*XC5wZVI-}1i+{*jo8IQ}n*8N|i;e@QGLF0TJeVg+$={|^bI5cK|U!5V5oHx*G{ zw{ISd1^3<27uewyuo!Srft>l4#Twy}^62@r8W=pYr7AGuVsyXggy=d%^+_5|kcJ>2 zIkr7n-Tr2is7txz@*nda`#x~htc`;_fM&!Cg`7F{M~p4+6B`!)^ln3XDO5Lw1o0` z@A&tXbACdY@5Mt_KT)?@l#e!DR1(x_2v$>oDYF+7)f6Y&m6*0lKW_+EG*qLCPr75s z4@>L<^9)+0QX(g|T!lPWYVGNor;X)8Fe8Mp$g8o4pOgvd`0x4G)ymY%Ys#;F5o%xQ z68=HuFOB=-?)1jVIn4P}ce(=XqPBJYZQj?k^JAVcgj+yuWy$qHmCwF4_vvOsqgAlw zxa_Uy?>3V6gkW#~9cFikUhz_nX|5%p;`5K;Y_Y&<*O!mgqnwz#LROq?}Sq5AJ;b#jx@k+hMt!bqSp?TLPUx&K3(c6%TI&5 z>hF}2Gu1J_kXm(WXL6pL+!%n@Y%5jXGYMh{sALRT?X^*35xL={wmsA?7DY_3WfHB#K z#f=HDyv{RT%@eZWmJgi%2)MG0EsS*NAGZ3UZ@XGxBZ5ahAw#BZv3ro3mOTdCHFb`h zk#spWzArrAEZGEBOb)^8oSjRlb?Yt9yoBgQbjq7*tLGw9%5DIJmNGjxu!?#8FP=B( z*YTR{EXgm2nS30j*hrK{yEgWrI#4Yyh38~8EX-#-l%0GR2#+qmzR5(mXJlyGAGs{1 zPxo&T0ts1BA37|$%UiLsuy{6flOe^|d6d(7e6bi3@eGN37S;bp-CF=vvSsVK zxVsl#K;abb?(T)VyHmKkJB7O!?(PnSySqCSZVzhrmOi~t_j&KWxN&c+SP|q*Su@AX z4E}$7BWHfJjSP+LJNJt>O9yBo4i5Q*p))XtGwCP5UA+dyRtfX$emADngK_;re zpiXKF*x!pUUvN#>Wf8blH(eBIfUhN{zPJj^*puI3YE7xBe!p6{y!JD0B2bsRze9E= z5WN3l!A7ofYMgpYmBV=%H!@w$S9x$snHhU!L0aE}MRU_g!r3ve2mA#K^@Ae#Gk%>< zz-3cO!HzzGUjb&gi7Te?0)MMeYqt$&=JYziF1H^No!@c4mc)ibeZx5ElByCW#?kqT zp@wBdr)U`Dw_P5%tDcg2%LBjBGBv{B);2=W&%S&n{pPMv!o%+8sUULt$qLK1zRoGg zmMoEm!?4vI-YALSkABvW=)s)+t<|*MD6^o5V?GxSi9MFzs)(ehbu;(~@^CRx35P;1 zE0OQY7{%vT8Y_0j#9$aGi=q`_j-mMIaZ^}*`E#!0IH_U>EpdYRiJ^@J<3^nfG&=K0 zaO=ZV-4>vn+2Bv6TjxyeavX-1XGiB_+Jx0`8`RX02U;K3k?QL%S@n*p1T3Hw>NgW| zvN1b4i_KHaZj&Qd##W>4md~>`BYe|o+#2RPdSWQgZ*<`5@_q=)kkrweQ#8e>Gs(4v z@lRr(I-@|N>`O<5OuZH#Bj4tgeKC)VaJ^yQt9P>}!caSMGM?&l6ws|fq^ioRni0G=;3u22C{iPPWJltQW?ESE3?#I3UgkrSQ< zI1|KewWMhN64Yz{GTR#4E5VO2lNcMYx(ktpBy|g*z@KnE%o-#WHNuM;$_Qv(ybEz1 zQrL?c8*tBW3n&%+tr=RlW4H_P9CafVz)6q=2xtQW;%z6N4Ws&Qb#(s>;N(bIP~(rI zB&7YBl$f|{hFU0}45@oX=C#73W!eqz2AHmovU{e*X^{?X!~7m#yM^nE%0rOkm2`=w zt@=a?6GmKICl_7lZ9`I3)_bf^LY7HBnZDfQ6@^kewoR(tvvZ~gmj`mU{5)C`A;>w#i3*wVSN(2E!e1h5^BW19&Hr!}{+ia8tS zVJ$S2DKO{FuT|#CcrfwzgbVxDqJ0Ed$uh=yY~!nr zK!eXy4!S>YmMX@j>saE+g|AOn#AePoH|OIUiA?B?5f8V#&uQ~fWwbPRt)@meQ zJ!ZdJ`E=|IemRDJ8zU22J9~a(9oyf-AffY*W!kqf`K5#ZDf(}g<`3yEX!$lZl8L1e zKDCmGC9kEO$?pxn_ER#^w>P%?hYfhQ0{^#t=ondl8}v8(VWnVc@;(m~KJ9PL<P_jEw*6B1)B3tgw|(+Z_!ZTQxoNUqp_`b}akNCRdu(>ZPzr zz9Kmq&ZSKA%tZjxi)L1@UiW$NOBBjD@mxYPs1|X^T*wImDJb+Hv+2o`5%CFdpn_~> zQ9prhXTghn;qV#?SWD7im)i9DvJuJIZg*-m%DJ~}cb0}K2Qfj#w&p6I)^u+p z0l^A)B%!x8xwB^Zbj|WJl<#MWmg=j=){l5;51&e`yY{v^)l&GU)ve}mm-@#ocPoS?_HsSjaw14QI{}teZNy5nr?KQ3 ztz1SrvsE&4Z_b+)IAh;TyR-PXd8<26U(00wd^UU$dAjE7vH8+E3py)mEV!AfbQ!PO z@aX{^oz&*0ML3X(Z5teRhGQ=}y6A>=_a@{MEqNs^rF zH_1RIO31|&S^&^`^~i*#!t#MB&fV7tlCIoS1bWr3pE;6LJEbY*1$k`hH)WR{NOjmD z)<`~VljU3Vzc9blAL_R?gTRT#%85Nl#gmYYdo&FQ4 z5FqijqEe)_xD->X+P>wp^@h9!HC-Trx;|E++*H{mfq()Mno>-hKfYr&5xt&L=n0ui z`Z-JZ7MTvGh{YLyE1Q^cF^R#qmKuYCXLB}bcQ%GFq?zJgj%zc7_-xHvMi1-`2e(5# z3rJPG8V=Zw5$!KX{!r62)Yr}4reGltvTu!t{sz=xGoddm*QJYdm-al=lrGRlOF|T| zk23n6Vu89HoO}Mjz{)M(`GHw2Wok;1Z!U%LxI`j!KmCyCg>G~jm^IP0ab7fE?9*eP zH;ud^+;^?J8^{g28?NAL`68{{ArtfwYsOYIbvZ^($c!;pkTh(vN*K5Rc488j>xGlm4qB{0w3Iu}*J3e99>Zznv=xzr1iBt$Z4D8~19VYgtt6caRPY%OI z_NF}zGPN$_)--@w#n{Ac_36aQD(cDM;adx5qDeW#gH#|Rt zJ8G8dIA10&7XhR9Lde1=H1T{FN;-YAtC#|+R>Mzmhu1{d+!a@C*xmk3%?{K^+jPSoG9_8oM(A{xs&)Y4(L{2|3yL&| z96Oh!g@4$MMi*d!embEKUf~+joo_94#|}sF{6WECE0jbE;ikAKJn!3Bu`S#P2KDVu zpaATKQj8usi7Z$$XSzjWchIReLh#cT2pcNuuqN-w=HyloyHIcCaU4gqw}atX?V@Jk zGh2N%9aH_;g%rr?=p=~S8CUY=s9amNg~Iyxk*(>pYiJ!QIK8cr(_{+Q5fWU(vtibv z0})HDf&vD~z2wBi6wnE-(atIqQt-oN%fOQ_2`R%u#po5Hi>I2sI^Mt!^fkKouQTr0 zTp1MAOc4>BM`<_gv~22OX-^uX+!l^dOcwMvX0n*6Qxy zEvyA0$hQ|A%zwUCD-oE!hM>K+^-FxPknDFXQc?~s#6r|PbR_7RE_P{VH#~0nz{)CC z7Q?F9*mo|r=uK1ip!eY{vPyf$t_Nu4Q{w+6SqJzAxe4T9_60SPalbU9KD~? zuLpe4U*HqY`G`Ej#IxPeEc^&g&2i!ghPBz_;|7^9r36SfmYGhwY?!SR?9_hLR+5Em zG_7tkGRmLYM(gR$`8yhUt?e!<_{J7o?QK}s^E|8>ffq@H%alw43ipBou)l+Uf_(Y7uyj&6w8$y zB+HzdwMH%Z1~4o_B}+Cbk$<%fx!ODn7^?lsv-@m7{C2TdhQ~vuiWBYpBU=|Ml|CXE zHa|)@5W0C|vCGn+{>KlebB*iMf*Ar!K2O`b>c;-vq2j=}pa26DthornEN)s4&qZ3c zFP)Hj-qqun*eSlgXCR78z;8FUA)guhdE|v^@nGqZ3Fe-OA}%c=Pr2AKDIKDU4Ci z@9-tR_0g4<$hxsfAIvrk(9p-%%7;|cV2)`_b2lBPh^3Xfra^32Y#~?bZ^Fy;Wu;n zTLV2lBLm|?5oP2+Xj!eZ~U@4KJet!|ey= z-3g<|`vQo~aejh0?%;SazR;p_)(r03y#df((T>=Gq|g6i&vOM1G>-)|iSYwJYLpfj zq9R4ecfgCAD(piCmeB3M0sD40t5_RAJcgVKl}wG zU2pw{dz$I7d|VMDzsY-U|3ZW`0v4c5>m*V~Hk2Qr7E}pH0K~qKGNKGj5a&857dSWI ztN>^hUVT}TfHW8Y&f5TSnNamK#CM=Q5X(etpJHCa36ZjW&S5)^V)=Nyx;K5IpFDxM zRUu%&GI{jj*1;e6?Ry*@=3(U$%!QExK`#LCYPdKQ0jTcWTrd4=d3452bk@2ZqkMP( z1o3$(`(lKm_s=r@L4Co`qcM=>0Zg<_L9}TCSOAQCfC)VDfjk^O%19CIn>oK`*5t{5 z0I)2Ac5>5V2bB3J{s2}3`Y;RVu`PU}ZH-R{A-JM^o@P#RCvbdw1wRuZRGj%fgFo+o3e`8SvB|Y9O1}emuA{0K}8X8VT)cEmI9Z z04~?poQd+78bhSu$%!UJgoK_HVFiz-0~44CEfghyr`v_0Ll=+>%{E`5qy*(Pyi=bOv9H}9JBOibqTc5U#-5MZ-7rcGJnv7Ix-x2#WRxhNs z_Zp#uwp_Cn8W=Fi_^S?P_l4FEJhqQ^DwN7-T5*u|_w|zg>)t;2@I-`$dTWq4o}9TX zDv9xvdQ4W-xz~YW^r0_0G5AMYL(x!7OAj6j;RYHuD%Z~iBo4a>#QicyamB~CkE9U= zFpoPAgnlCiyTil=y1132_w^>~<}eFeTJ26-VlPY8f$X3CQk_Z^NR}`f6>Z6{b?1tJ z?A@N3wt;Q-%1jG((O#3s6TQ*SsC@0La+eInab;T|_~)J@+8Ve%Xl>AWOsQHK%&u$L z>Sd_`nxb5jlUpz@06f&A!^rea-vY%1Om!jW>LK2D&*og*jY8$97 zVa!Hu^{Or%6~#g-)|_~H{o&~C{wnCLo{>~OOX6bf>S{vb-{|UE_cv;#kZU4> zr43mn7#ZJL=uT(HW2S$MS8}X5bDLermHF&eNg!!AsdlWf9bb@V6p^;d^oUeM%09uI zhK=k-rI9r}Ry|Qv+c97x-_z%sKhbpuYnQAMYe002)`}(KYu^Vpl%LoU!Jd(qv5yJ! zQ@V}sWcBvGC4)X>q%7VuP%5ZYjY*T8(a0jI+;XYB?XDxq``*2a)~QrMrhBVA#Rl$}#76Hs^jL>MA3<%@(0Uux2^= z!O_Z#Y3`gd8)K#19y4e^+yeflUWL+nh2AhX0dVuV!NL-;QwPN`*URyoO5gAswRAt1 zPnt-&>$Ru+!4*yZfg^ndrAU~<&hbmK{CWjtT!$m+(o>uR{FUJR$@G(AV;iZ8{ip%g z8-1m);APIKuGhKQKuJWL?I@D&y_~p@ON6tZn+a4NBxuaa(yUKxTx5>S zi$GgRbZw~feZ`QuZczNPnKdq>Q*E4uG<`=jW)qTG^L((E588|#PJpox_5-<0PRP&8 z;)>dUQXz3zsESsf??oPMJ-d|^u*)9*)m(-K?>k3pdleVWybRY!)@|6D_GMS#hm7MB z#4Qaowwlkqtp_mbmOJA{SLBTzt+TipO_3ISgrD3demFU(DxA5N~|vSjVAu8CrC$-ZTLbvJQ(#BO(zlfD5Zg_kh5 zz6=+Cd)c+L-QaUD#588mS+W8ud>BKJQZAu}lRc;XX{9Gd*rM3MBl`wAp;amsO#zz8 zo=CMmmqDE{Z!Lz3Ta@`gGF2Y0m90+HAvf&^joJW^CmE(ML6BJ_L$8X@aa}$B>)ou=vx`UVfokhZ6d20s#{SroCt&Ih3HktT|2u0isH9LU}N5Ih@s zl#<$3M#TzApy?4hUjCA5fRTXA)+$R&B4rjHV!GU_Il<+$rla9rATww-HM{VY>~X@M z!ZaNeMYSzRE{f>_%TuYtT#f+9hOAYu(ZWy#WI7}eSrawm?(I=%Z`JMPanay7WG;@g zUe;(!6G!SI*S(r592sGf5p0e|Ra}PJmGoP|63aS?B+ZZc^GV#WBvLz3BFSdQ*8_yl z-i4jM;6l`vMVSkyRZk`9v#K-l6>aGeOUB2pY78vv>0eZ*`>!}PWOGER17S@-Tf(32 zpXr+1r@%tmA6B63ZBF-{xm-_od5piRp8Z@ZE>oJ1xF@)GXwkd+t^vKcA-Z5VeNCx8 z+D33jPz>E~rYmS!w=7of8nIBZ?0rkwgG_Wo8?FTg-OS-QXD?mksx?#)u`J6n#gTXS zXi~hf9dGZSg@yW5*aMeiPAko5YGZl8m3E7ldZn-SuxLVOh!rfnB?LqMsi^qqGT+FQ zXKSVOa$&{8FzJJSB@XS4QZ6JovkG&W+-Z@!5#E$pxp;}PXRe^Zt(vq-zN})j-=KNV zz6#N5aFu6EZb`dN#*K=3iwe%>PU#&eHuV&1P6A7pR_)qt*H`B>uawmSY>xzU;Nk~z z_Su%>gpf!^zOS+5{ptoID~1v~MEW%ceJiQ1!z4uwkQ7wSOcq-fS-ZVC482Jqmj!0j zzG{{vN1&@by1CczVSzIO3F;Aw{kD9#h@V#RwsD$r1!CemY3kb(^jnfw5tYj}I=Bc= zn4YjUoJN%K96vrb&NjigD)k!Xr4|3^*Ttl7o~esNe?^0}T;a%7D2_$PuI(C5FJ#k` zA`24f3Vy(zMa9BsM9EyPRQTNETaq-;hHI37+Iz1b^h|1!#OjTBZ+%#PvH*Vhyb7l; zyeD${tmo4Bs7ddpkg}_Mw6v=wRl}Uev&B0Q9>NZhGf0r#JTEXZ-b;UkqHFinz?c2P z+ZFB9Nw5XK{PRg%vqPuafB`Z}tQxyof|I8hB8mJO)DWv+&T9ht<}2ijG(waKesHf2 zNVKvy1Q)b^clnc#k|D?@_lzbK?!P7=ewnm?Bp~SNS^j{MbicF8zq~Jm?oSZ*cZTkt z(m682mOc0|Emx>qmhcpp%C(%CA~QEJ{CK8|@YI~P#ZqgUrSDS)Z)uz8)pcBxAfC*7&qwHwvrG!xCVhO)n{`WZR&o}8 zB7!E--rV@Y0*HL8%mR4u`29@$)vVwtwxwz*$TT!MpTd<&Z^yf&d+E>@enDZB+Pu{Z z{%>OaJJ5Vf2+%S8ne6#{hyUO5UpY}>C1n8$DH99buXb|)Bcue_=R}? zHJ?Mv%<_h4Ss30pD?KYdEAu~4AN7CBci%SoOH$yS@4j6g@41XW`R?zV{Aaw2PxoJo z5r5Ak{#j!Ao^$y11M+{+-M733?XRXk^BVu@^#6OjOUv}@Q13tUXK>#};w{(qTPo(y zyvYA2yi4~M;op!iBO??3U%&rvK=6Aa=>MP3dk_9R(tGcI<#GNBcYodHe=QNE`|Gg$ z1MbrOQ8)XyY25cq;vt4{swp39B=xMJlSv9EQoU*=`~|dS6A3P z%Tic^>8s7y8OLk&x+Sm|RnA;uzxgF2_@8@igY!&~7l9yBqasP)WyuxOGsW{E5SaFk z0S5)Pd9PPG?_)tgBCxm_Ih?Qr1IziR-MQ^8|IT(1-`Fm>qJY59FJ&4E1&Z9kP#^cS zaZ0JQ8lFOTA|rX*$pjwrpcC3unw>2V$1om0zN|Gid?|By-3Nc>vwy;Q7!bs`IX=Ne z$+p9ISnnNHDz9%N>^Rf+3>LRi9NKJg0()@YzZl**8s6mXUm7nSpSg_B+mR+N;PQI0 z*k3~h+eC_vcY>X$4eEKkC`1wZ?ApLQlOlDUuSl@DMiDTNaIwg}su?`vhf>~kQZ}!V zCE()bsV)RdrvK%Rik1-RAQ4HdS=eCN1kXPks|JOKx--Afp+D|-y6=!%R7G&OeorBvNKmge>Rk~kGd)Oaz=kUmcN3DO+dF?n9!FJ+T}iOI>? zI=Fa-yR9+B2m-Aog{9jM6{NUkdJJ90J9I!M8r5^KBD=w+^c+o2 z&Mm173ycwzWA`lg1UP|u`ZKVV73G#H$rY7%GK^z?o+0)0k<}|cC|BIFqnE$V0Vi~L z%1{mscd5rvFBOk>t-6<=I5x-=X;wJlkPI3FV!|j(6)rf|Om4`!^q~gAPd3%AXsJpe zdT}DkB5}SxKJsAg>>xT@eBNON-l4bJ{*&3MaMDAxnii9=D2~WC;6nQ06}F%XHz9D@ zZ)1n}eYTNso0O6`H5JRnz97S6`tXW)HkwJ$!b#phgfsk79Mux}ikv=-Hfen=n5~z6fW9q@#OL-_@<_PW%$W(AQ5|dj3(+?+%>GV>&tuOiKd06R z`J(G1C1DG!LfG%IOw{g7VB<{QujsL|J+&9^S3vCxSD@(PmtLC2I1-CFWyFB#AG*p<#owfj2a4 z>>&JxhQsxbenZ29omJap(1&-5ffkByXqc~fuTm8ARBEfD1Mrzl)#maaXjl&|BEsS* z?J}8mQ9W$r31^hsfDwB9haSZaken0~(_4k;ad#UBPLC!fS1;GwL)1676Cc+_C57?x z8%zbK6(rZw7fYr-jGTzbdr==Tev|^102Tqum@=cQal^&0ar0(n>Y(_cBcluF`(c+R z(O^{6tNMoyO+Pl+)eoI>!QC)mw_sKn^K^VO&|ViOM_93_=o92iuj1&t_QzJEPiSb1-#Z$wrp5ZyWXwN2XAr`>hY@9rXF z*YjERc7Cm=I{4h>ZiRyTk!pCY-O+4J=%D;awE9eK13hNKEur7;VzEz$HlBf5?DdJ% zBp|)+730oJHMq5Vjo@yUbJ<12LBMZnRjLiTz z;tJk}*RiK>ZT96j{xA=|IVZHf@Y@)B4HbeI)-bt;^bOiqtMN-x#AH;~?V}$GBHX@6 z0zb(2uwFq0qN9UzkOPb>IShCUy_4a`LM^$6PDnETytn7oU7Sxx-ek+ieV`0B5|@D^ z(YG?uxnzXVyz;YXb9)vNy=d6*!sNI4#{0x;q2B659Zk<~3>zBEY?P#T8yjx#`z$it z{Hma5S6sBW(-{5_T{_o-EpHlV@b^BYL zs2Cr=yvqMlo#?-si2no1M1LvNze-m6zp}^wj*$KvDEx2Y^uHGh{~F+bfx^r*Z@T@j zK=}9D>^~W&{|bcv&3XQ5eEuU4{uS+i0m5(K`kx;)qyL-IaQc5=8va)wLHkyJ{OfZ5 z9SAcqGqe1q;ctMjv(XonrCBO`dxEPS&7`stehVgNG{wbnCF|bb}V1l-A zVnQ)q`nKHK=LA%1~-%dTUbYprKvuCJ?33VL5=((O)E%t|>OpUD)3EI4sh*!^Gy zAOVFhAkHtv0g`3VvEUs!Y1Gq%(KputmwmGF6&2k{p}_9p(RxfFc$@eUt%E)L2Hz$k zo~FtK5LyX1Vs`E1{ap>#E#*$6s z<;e+z4b~)-y1Wek!UxjR+Y5-5;YsHO)Pr*f06zahY#r*trwyTm4uSm%a$v2C24fn{ z+Uvf7tqEA0ojqv93YMP+W)p^#S4T!k$&aMt{PmIevcJ&d;CfS%J! z00jBWtwL`%9y3-P0tj2y+g<^e_%|Hd^8nf5mAHUerbmFYv(Jra)Rzgxl7v=5B5LpLn630y?ePJ=Z)!X{o^t=1v0zq-F14Q|c? zJ3z4f*ekb|tbhhtj|JzC4&vq@>*;<)CnE_w2Xk`-OE0ct>f{GT_(s`ZR6r zZO=YmjZB)1SGp~3f^dA#(~;M-7$_(=HDAhQz?(54U^f}ZPle1qQdl8XU!MY1UWs7V z9T13MpaP8vm$Ta2u>c~yf;wJ_u5q&JMy-RIA4@Y{*YT~J+uOfMBU*-5_;%p{@FL0r zgbDb#^E}lvp`aiFx`BEvod93F=z#;cwOE53IzX)-dS8ND-jeU%RzZgXxs@Rhq3GOV zy&@w4-=*_K1+@d%c)tUKD^bJY?DQ3r)P}~iU!m= zni)22os@_M4qDKE!Kw-YHK9j(TP_P#JnL0kfD>b4-?nX#7?ODJ7dpuVD8Rjq;yPonjYUj zi&{Mf-p7Ye!sO2q$ksc2lDJPBQIt^Wzmhvq)SGT-z3ya%K5?Z1AE?~yBv1xbkhm59C!i)X0j84ttE;NBxDG_4&Tt-lnYIwV8JLRD)qNiH3P3BJbTu>$u9i7 z$tfF`r*RnE$m&M67ogxg5LSz#J?- zGNTy8bo=v6C;dd%n0c~pC4tAM-0Pxm->bE2cEL_{?|J0#r@2b+*FGz;BgpCxM>>PN zl{uL1s&g`EmdPXRq_{9&v6Zu{b%Q?{oVvEiW#l{6Do+&#tW%ooh*mp}iox!_%B9DyQ z5i@3)Z8be~&UYcl7t>G0yp|?q@21B<(nrJK^fmd2>tciT4AEl1YC-o z-gwm4dGe?g#`9hO6KVx zOUrYJk1|0}smYxtcX5%UKaWc!l+*c2*YL1FB=&P&fysFt%H_eJCqNPg!;$C7>ykV-nP$EPsDja-3E z|BWyM%Geg91+j~h>}B_H%Bf`8xuv~xJ>3^@F=NA{iEqRzNPu74rUEPXk=0X}WmQ3- z8^oNY>w|b$I@x-~2MKHCBV?_L*Jxi{X0nf-lMm#q<)S~tn(BK;43Ketrgc#FSo$28 zyVw)=g)d%9QxfYdzr(|jjV7|ZlJH!izFd$a!(h~tH`Z|L6;PbUym1&ZnR9jHP6BC~ zUj)TkOXFCf{^ZZR+V}zuvY%@NA^U*cDcabW+cE`x8>(2Ki zBnPJV&%Jz~B?R{xQ`_w^!tJCoQxKB{_%8IWtbOU5JkDx2nxE!~7Y|EqVw!P|qg246 zOAX=o2v=(|S8SeKrvW%cKZa9*uBdoIwHi0upLYuiB{4=G`du3H^1`QA(p(#HTk(wz zs)8eMvf)24gzN&BHI!+L3Qm+i99L52p4tb-@e%Inm!p`jAJXBV4@%Fb<(2uEl@T82 z?f?O4(2aW%r<6cZt`E;)n?;1+%DC~kdYlmsgwZ)`qd%Hf?~(&YGJqVN1n$Yx-akk! zSCRXfoYK}<}SKmhS$@hmN*J4jzq6Ww9+M@3pd{7m!Olg7iE{- z-ilS$Dl`yW{S5^d39&UM%CG{`(qop~A0a263s_rta+NLco{_R~yYkuMnUYAAWWJk6 zUxfy!JqKd@)jGRJN+#($o=mt)+~Z~4J|U&x#5@|Fi*4mK8Wtq+**ouw4tsx~)40*0 zWf}UiQS3RiWj^fM+ls-6y#<0Cc~DRi>mncmW`2BMswbn5gzjIO6`*cF?Gt}BcS<>D zAkf6NHI@D|c|@rik$xfu3Nw_Z|MdI$_e}>#9OF2_@eKC5$@BG;u9ZA7YT7`WAE zDKy)mt!`?1ovJFamXq->R!$(8m<&WU^T_?i#XF|O)Uuj;dgvu<4pt=5HjPpryo4Sm z6r@2;kCq(a${LWhE$Md4yK;6xG=qae#Ot<1zYM1^nUbWqSoe!-uF#;CQr|fYmU%=W zh!|51-R!mu?5Te}^w1kpUEdPr$}%n@jvRTZ_9dSygI(+3ss?PdH41dtm~>$AZk$fc zw;^9(*U_OiMkC|D=d`F~iqePo6&Ww$quE^eL2r{jYy%YbvocmVBV(h;$ILg`!1#;u zPbn9B^bw7b(Ui)#k?W~ZT4;v`7Or_Fa1iI7l@BoUW9dnm6b9#Y_)BnRqH{5N64}F% z8b+~(n&G-1AM50q1wMJcakoCu+7Ep-j?IwLgoGCH>a)v!^P-!*LT7Ppg-t026sfDt zJBh2nbr5!NHGNy*=rJJDqIq2TS*0<1;mzl1H)E5kcN0J8E1R0x)L3R>ehM+Mc=65s zgj8x4F>)i^@e?!px#7whI4yo7W+F$%_SJ2zK78cDV$V@{Ii+y8-|iI z=LREqva?wo9`T^s$Q_A160CHQan09=@TJ0s!pF%e{>oLc8V|kbACR8-elt8svWW() zr31s)nR3J#AS8$Z<}lW@3;eYLB=93YD5PLr(F7(Rbo%XtFt!()AJOYEM|vu51`cVi z^9`}S>jGvsegeB&ZtYstSd;YI-T1!aaK`+6%+n&PUrv-fDtUTJA_${Bc-8o^;zkiF zlduA!>GCEtmhXqv1&OgW&gXM}#f%jy3bQe9{2ek=7loHNLkH(UsS5ky>O0h%Z$IqF z)C)gGrk8T!9iI356|2o9^1$`vN+z#mcm*dTSU^A2Pnj}eUyb4c$Bmyc#_OV?Wx_625PiPf%#hHV`hRjJ#`k@a{n?XItvk@ zEUQ}{llO&Wf1}tGxT64GrmvmZb4Fm7IrvHqwG%hf0jENJ#Ay&DPN*>#@g(i)85Z~l z(QmdwKB({~T#IO?rfVC>8pE6E;tadu-JMr!B3EwO?21^$1Gi-o# zT*Jtzc5Zz%HIh=)QCpj7y!_0D1fG+7-r1VyQm-MI*Vsf>B{G9^_9%awkvp??Rjc-+ zUwjay)yfB5YN;k5q43hbd>TS`F&q1`aDG-Fso#q16S!0_5;q+JkOVBYTXlKp_YF7t zOGt~y*5N9Jn9@rpF(cRgVePg`n<6rjg-=xE@Ta?^3Wabl#p`3RtWQ9DC%yVwc<71S zbCdEKeU%(pIFo)bak@L0wl=mQ!gnk%B`ru8G^;2hMbMrr0?iq@F{Q7mhPq zFlC_D;QKRxt`)-?680#4NiV1hzres1ew|*C54PM~PS?#<)=sQL@(Q%0S`G~eE!3x? zZU};OA5L09?QXi{z^dV(1Bv5vw2WazrIh0}~ix^woq)=sNseT-Or# zs5QaP3S235^JwsXqzC|JV9)v-LBvcG7hG@-bYBn;YB7lkrztfYT(QUzQz=vn@YPyF zCXc{W4$yYbkDf6QGw%$6rX+o~=S;=n%a9>MaI(|$gPyt)U%QNTrrqKk?dvK(Z0PTL zscki1jH615GFD+le@&ST85RbwkrzreON3Ai8Dy*pv@@!Te7QG5Wdtd&=WhBC~}5%V38D6b*9E2y_OZMh0Z`LHbPb%LNn*VL!d+rHy zUkz4Dn!eyRWfX^N(nXi;_wAd)qwPu*=)mHgkgekUN{rh%KL*A=)r$|NhO+aYf=7_5 zm_Ehtch7p~B-_{xMRVm_m2WG8F&8(1^jahk7-Hq_2X%8-O*)TEF6u#yr=yFl1-8{G znO1*rC`m`9uYm4dF!{+2WEK%M5#4W|+FePx6{=&mk_;KME zALRBt3KEWx8p`jFdmi__>h~h>f@c`i`OPC)fux*U{HRVtR+eV`BCgY{qZDRN^1j%2 zref*+K<+}Pb^wPd5%XAzd{h_}ZJ>&eKKSoMqs##(tOq$ zvPG=S)m(;I##C7V)aUZN?+dfNI}Od!k|oh~i`^tnRWfvXX{SZ17O)1qk?mh;Mpga|` z&gEI~(57#q@Dxpnt?_ji_R+b54R>UHTg5c-Q|NZRK&pRr?eMbdfdBUE8WeBw5b$>q zWqvgjl)ayujn4bbV1|bY%N;jdcRu~05xG?usO&2{6v$e=!z__1--*91w1X6Qw`qj+ zu+0>K%NdMXoE5Jl`rpM0(;uRg8Z~fn_6-k*2Y0?rRqi0)fI>;@pI!Fau|OiZi?lq{ z8?(><0)t&F7lT@Ws>b&=_d6r0-Lm^IuO@ra=J)WERlRSR1Wo!}kA&wa#mu`^FWF4V zC7m>A8+lWF!eMA3;xg+~z1``WbzkgG81{!3DxOpGTd^3Th=+(e%+XlR1ZYM^jB(Uf zsYBa!KN+XK#PChcFqwgMlcdagfI+C=;6iDOU>IzEAHe=*@R8Yg8dB`|v;V|VyIo8D z7J1Ktc=ttL9v$V+ShkZc3fiZfq`#M4@iosPLWDRVaE-fq`7OOW;4Soqki)V$4P z^iu;_mn%kGVS?rpmD{cBt*Tu20GQG!!vO17Rkk3GMel9gCkt=EY+CB!P~oYau~hDw zErXU9TaI82ZJ}_b8a1M?zE%30f}~`+NrrXUKYG-t%lu4p4OejD5EFP6^X)NJ@b@(b zV*Lbetho9=^mqFtv#qO8%3a`H){7rO_#yMvI8 z%d2(Ugncp22&+_)k|dN%Q-^)!C7nvA%jK@qPxt4)<+m*j?31LT*vJ9w?TWL|!!o=i zdZ?9Cw?Vw`!y$I5De&-@-j;SZl!l%sm;6{V6MPD@RO+!)dv$RaR%aGD{9S(qk)oJo z&D1sjxwq{q4f$r}z>93_1DNi;+@Nyd&LN*WvTuQp;k1wBb8Lz_n14d`_T!k3(74)i zd=7z-Xr({s*qkc#oQkI3H}W}`(LdsG_(qA^jRh4r__pGb;6hVG=fb_ zB6t!r-gG}9yCfys~c$mTE) zdNu6#5|)xLnoO}H4U6@u@@i1&0GK0qblhz-&|O({S8Xm{Ez_ISjbAJkujRHTnQQLM zAg^BGlweGP$x9{{t}|tZX@35pqi@(P_#atpLV9Z{r&M@ge~1}{d>^rgqp zDNg7AY3ZsnC*q-fD*4Gp=hagMp2ahX^#EYpZ`88#`VC7BxI5hcm3T-POX zNjFxq7A5)-Qk$fQDSYjg7_OXy<4E;b>8SH=XluK2 zXv0zY38>lq1{3iXRwr$v6T8dgPpfxp8b8YmSDNNrBk9S9c|kK)ZO*1XHzq6U2BJ?7 zl=1FhYE^Dpban?oHpQ67+4+8X2gzuopn_uCjZK$QrrmIM!?$1LX3vc3uQSvkpV!d! zh98@fhg8@qr{;eh-*a8%T+G@M__6oeX;scLdxu@^ zEyf;Cp0PZ7GhWfO;zF)^Qqb(!oqn`y6t~0NXo_p%`z#ZmaxZF`W@X3W;Jzror>_Sw zP>xOh;i>D)kb8{Ch7SFiteftbR=3H>ylS?>d&}kWNyYV>t)s8Tv7FNAr6kRFdxyPK zA54nol5u*z z!Lh~MZdDaSrd^}OwN$4+w~(9c_~{RhK8o^;Y%tX8^BOy0^72S%LuKXUru5p9t}1?j zoRgJLsz-xy>Pe1eEF*Wsd6eEm2Ues^fU7d9=4|xWZgS#|-lg|F>~#&Hm(l7 zl@NFPAAL`!jjK&Fx+}~c7BYjJ7@^>#lo#bJ9aFUC`O ziLwFpmKsRqY@NSVZr*#k{Y}%ks^Op8QLiYL`i8bz&8ub3s&4$nU89I|SgO1o2nrz& z)mY26P#soB1Z85HTb?EQmR~;N8~d|z(gkxzZt0x=`o_im3>w3EKK*efZo1!ftFfaN(sW84zx2}*=636|w_xsMmZ%dDl z{kmd)e&j3uy!1RxuJ0OZ8(LvV+TJ|&K!!@Vn)KUu?d#>HTBLk7`qMwdU@Y+O@PDO= zL>%@bgt5SIFNRq@5yXEFvoye5yh_6UqV1j3YH8Kobv8Nf1#Nc9oQ@^iG}aVrsIf6fC;NW=(lX!&;@TA{JXFfKy!mYROpjU?jK!~f>8Y#9 zr7;GMtCR8MRb(PrePVTZ9?89nRu;W_L6kf=|{3QmDLr8W!3It#zorK5!o#qTK z5e$fMHbTuAB(6DIDETiicmWr&m|h`$cE*A;ei4KJ($xPFfyWE4T#SGhbp8KAz#~+b zkb)Pd--RqX@>7CA#|pQVe3A)Ba!!bRS!lR(3iL;`1Oh>`>xYldmKlZ!nHM<|V1wPl z_eju5TULM{EYd?6crb$*EWh2TMRGiNEDVSNFaU}HS2zcTgzZD^g9s=C0?uT6Alncz z>9KbF1nUgGL$Efyq0pq~!eaV|?FfL!6#gzaMIbX|1LC(ffC~u#5yXNx_#M$00M`M4 zu#U>~|5C}PF2I{K7Urj-%99$+B|x2TTI@`J6DEP>(OEif_wSu+TM=E`KUQG&H6ydxqy7?j|CBDiaZ Wvi(?W!79T-Bmh*<)3c-5EBptfHALkA literal 0 HcmV?d00001 diff --git a/Chapter5/Figs/figure3.pdf b/Chapter5/Figs/figure3.pdf new file mode 100644 index 0000000000000000000000000000000000000000..9c712a5142b8f8268eb7e21d9b30090c13794583 GIT binary patch literal 52511 zcmaI7V|XRa^Y|Uxc210q?LDz=+vdg_+qO^aWaDIG+qP|Nvb#_A{{DY&pBHn^)zjTI z(^b>e)%CfWN=ZzDk(G%Pf$GOq?K=WHfCb=aVuK(c0AQ9gcd&G`0&slZQb7O!0L&8B zc5db_pZ|8oZsuao#`XxFIa}Sbb)38<na&SBax_w+%1l2EhH zy#eoJy~u7=0VFA$H9(DdvX5zf#P@1kE{6YUqAo0zP%%EgD(`_t*Sx3*trjGKRdGfv*M zjTdkJK3Ls2;Jb;xxNQ0ZfBE#0U-Ad`53I+!UQSa9Bl}_3OT5R^$jgY=Q+&q@uxE-c z6Er;Hx2%2sTQH5@RK-cL+CRlWi9P|Ow)It_%C-WSrc$FoKd$rd&QGCdk7^Kvwv7IW zlzUW6XfMniI-9cXxL5pi&$BFtuY5>^^@8t6mp+t~c ze$|0@;glpCx3Ct2!$l8ej=!9Hav{gYLg@h1F+YrCg^{vlG~Jk##^ zUHM$SRGxlXyw|f(z1LHn&ZuaS+^DNEfmzuqk=a0fJiWj#ky%4MDxJ1DBE9@PJcxpm za^aiMuSHK0&&7LXMDCYj|4Y@}YFCRm6t;loiiTxvT~d+t6zWj(@f^1uE|Kp29}+c) zp{3jIGb_N!96Aw++$P8A+_f)Y55H~iT~P+x0+otHWFg&1aS`=NPdMLrOC|hpL@mzv zLMOWs+rCmWVg^=&u~1zbOxu&Y7Pc7AeiR!taiJ_K$MN#QUW6jy&c9!ZGV3`3$0d2! zO{==UKR$3+P1c}Gb&tie!Dl!Acy3piFo$#d;Ok&7ma=i!o$*r-zUB+JKL^=VHHdUU zDL^C^tv|n;619?Ld{;gC`QDA*Kzy3>9GbzGD5Pt(>w)efJh~`EOAdu}e{sCI4YoyU zS#i1S9!(-|?PklQ0XX@ie3HRI=uh@Sw82+KK}|rPHD7Fuc@0~G1{l(q>W#vkEtA-O z=GEgKG0~Bg%w|?MJ(qyX7PUtv;2oFCEm8QR=Nh(XGrB8QPPzih;sv+aIEf&Fnatl^mw6o*c?xz#>Geyj!1Dt$!m=37xfR#Xcwfu7i9yJ7fO%qDQ?;y%BFIO?T___F z@>4+CG*kTYfSV-Hszs~ePUI?8I|Z;|bXCteAFAG^98h3nCg&cKI~s&w!?oidYiKu6 zo0GO_+dUbOtO^c-9{?{=NlO|D2VX_KnWMbt!-?XBCDOwh2H`t@Q1@RE1WG-$!BwNQLP0!kZeweHXCr&&P$lZF>>JF zHj|c-b+73{P`iF(fgIj6AD!jN@u6rs&OahqUKl`|{4#k|=*9OI{{|0HV~E*46TR;; zj-t+?`3w!CzrDa~`tkr6!?Pu`*E8#0mcKBAC>u|)&Mpu_2wx@4p1(`bJNtnTk*T+_ zIn2p=Vn{=cx}G)h9Wx{}{I0>B?d;?adqociF2UgbAzCxBv|ROlw@P`QD}(YJkwQI; zub^Y2GdZ#>^fa=%p!iV^`6V1qSprJwR3$1VJaKJ$^o5!U$!La;@=(%T8r%0^w>sA= zb$Gk#^Fnnnj54W={qD4BXF-lo8iBu!}xlZXNGI$nnB#2X68f{eNwL#KsX3V||!jr9Og_vB~qP(6Ph(MkiS2s@H1{ z!;p)$1YNsOGgEZg<53fscT%F!NN=%I;#R$q)(b8N#sjv_Y?SBI1# z*fd`l`*f&=9bYvZTh#4kGAngYmVL^9SLn?$<9T4lrJ|VF|z_ z&~j8gTo36uIft6+w5NY~JFj2tj4mGw@>Z)g-H7Z>@wdw{5yf!P@p3YmhUmfVlZ7R= zaVeHH&XNjWmi8(cKAu*uvtL`Fo@e#5Sdo8~ulZs03+u{51jb@U?ZqK>*-zEJfQw>S z5o!SsJD;-SvU4$v?!-|A4$4&Z=o_??lw$@v_JRwfS_u#E0|}dcp0EH@>2az?##{7E zf*tJ)nO_}3gEAh8Z)eH-$fUZI3nCK%iU(%uu`+D4^BL=s(iLh#=zDb%Aq@>g0-^Vb zSsfP%PlE1`-w!zb+K31zy~gkmXSYQhW118*f#oHO*Mn6m{ADt(j+e1KWuKq;vda(n3WW@IQFAr8ew#&!)0N23juIBVph@3R68 zpz`zf$GKdDZ58n^c;9XE9#F^lfW#BKcoy2&5vt{Rsy>aw)NBPl(|3)mg7_+{1-~=k zGDZWgwHMt$Z0d;5uMVT`wm41?w65K1v`RD;8%ROKk+BGWDjgs7jFpXlOIN41t=U{tnP)$(}9t;*FIlL}qOA%0Aiv&^uTC@z-ww5faqJO@J)9jb?P zn~e~Q)SD^ZJm=Ci**-wKz@dLjDVi+!ytK0%!K3^4ZvOzYNUhEKVxcxTB0ttlTxW>MV~yc~Qid5O`GnsC$KYs|=GvNR*zm>O3oo}B zGvqtRJZvt_vzS%&1MD1omYxP4Lxb5!q*(p-_9=1qe+7+g_$Or(BCIhXgcB{w)3WG* zyLM}yX_k`6fg!p&{Iq~cweI0yL+XIPnxCiRY+tZJA9!#Z#lqJqTnGr{7AE3H#t#o7 zJEh0_dWE&xW7OR3t~?%WKJ|+Vek{{zD@;$z&Q)kYBRn7jM()Vwi^>{!FBfd+@F|_| zhV{qXd?q!V7(SX!oNIChngzU^W{pnU)@%yCjVB2fJIyz4Ds_ogCFS2N_b~c7gr3G(V<9CBL>5Eo6%(r30X7wMTqRP#{KVjv zn2nUs3Q1rUhM0LE0YaknPbqDlm|)C;Lx)XhFsstyi>8h#z!^g48kxBcGY4aIW?IDi zP4iJlJ5MG&(EQ;z!$ZeSl4PP=k;f3XZmXRTq9nir`Q=|u%tkl74e7=Fy zq4+5;#JXoTweAOn;!<__*pzn1#TU#oJ(pF@QLRIABLrjh`3}a#4n?}Z2O8p#4_HuZ zBDr8s&?@q~N#iUhb&>hHVV>QD-v|NRZeg(oa-%ZBi})5(dfV}?SuBhAyH#buE)J~o z_$WIT0dh6Q8qM#ezSp5$MOIr?WCOnH^X1dO4|fBt+TjEB)aJ;BL)CcsP6U#nAVef- zw)Z)*xf*%-wsiQv@$w$$Y(3H{(r}E_8X$Y|y+kJ)!s+Lt=`1I>Ny+6)zgnQ_=PIqm z_UP6s@MzaMeAYeuoHv;slS5((YqsM58v8S2HYfAs($L*6dH!{vQNg$Ei-wmI|Y7~^<+qi({a{@yy#@B(Or{Rb3i;Y5Ui48?!EhU z6xe=d-n45*vIBw_FQ4(G!6P|XQQ1OFyCc%j!h#plhDZSLAY)Mr{*4=5r-bk8n~*&| zA0OFU6#_CJADzNjwDrX7T>G~Ltbj+RCB5Z+4+ZnIxar(%=JC{+hTaE~=oG)r#uzg?7VZbKWvmQANvU6qY%jT$=%kg(FzjWWa_MqveltZoJ9*pmC5a~2W{|@)G^qj)PIRnxa{YIl? zbnMKA2J9El%)-?Ort_w1k+W2Y+@Zv;l`G6nOY71c$RUge9V1P#48UKk9Uf9=q zWLPh-S-zL8llx{=QxM$Zx(e#r6Ss$pl}+_F=Rw z&gaE;TL)vXBFV&aS3leaO8=7Dq6CDOLN%X*=1*2rvFNXov=Lumq(g1|&y6<|3h zkS`K&h2{j~aq%t=0;hLO<>Elbr$%$qbX9eR!sz{hjkYBFgE|bzu&O;c&gek+Gxa(&5M!4?5!O)pU?*LquQKK@9(;Q642FS69=^Z9Iu?SPH_SA=Ju z(imU;9RXzv3Ku+f2IVk1Tbm6G^IvDty_n`5j?rnfO|1B`^VAsT3C2!6l;cLWsz?h9 zCJBC&25N%xu`@|#W{f5SYQ$+c3!p{6iLQ(uSv?({7(JP^+(TdqnnL$~Bj{ z5x9>3MzqPvHJpCP9L0eqkyA<59IJD;$O{SttMBbO3bY#s??5hc50R=?t`R`9%1h$N zgPk|w^NNDS(x8*ioNw@1PC9V|&IqRBo6r&C=|wUqo>_4qwT>h$%#pwJx$>a_c-Bp& zvt+ZueRg!*iiphjat_XZLg_$`*KVjC05Mz+Fn*3#LNt z1b9F5mEUdS*|hDD=tJw)G3`rnrHz2$$}YSr?q(r zkyb1K(#Iy@r2>sxdOVnC(t#EmdvFU>3H6{A0ey(l)a7(GQZB#~r!aWh@KcgafF3U> zOs7)#N5n^xb?+`m?$11OB}LhoXsp1kGIoV7y?b#Q+*(*JIiCg(WFecEsaEp-OwaoHCWTjK6@imar zN$nRtb@-r(G`Opq8Oz2y4Ga_j(#SVOJWws!d;PYlB^i6WrOz~~`zbb0u81E_8$um^ znnS(0%>;=ja0|cNBCXZjB9WmYj(I}i_h5w{F$Q(bZ2)fCU=tgjnqFk#);=alL?oCD z?;==5v0bk;LQNT;v&qcPtdY{SJy-iX$~sBOs}cE9TZn`=n{@fqpq#G&Bdth>5Nm=m zoflMQ^nk2gXWiS93&tV)stWNPa)*g#;ZP<Qx|slz5+u_^Kd z%LWP@bSJUe(4kDheYj9#8x^&>Hy1ij0cYeKG{0$T z`}xaoGHn7Tc6{6^n1%<|h?{ku7@V)1w9oYY*Nd}UIeg+x9+@US@mXgbwf;Q2;b+n+ z?(MuGN;*dX^Bq%zD%;(i#ES+`Aeej2amc+3#qu%*bmwWaN^LoMRjk4}gm?-YXQ z$v|@LlEs=#;sERfxEfI$ljFwn*9T@)>Ghlyaoo8nIEG-=y3Vk(Epl)UJY2%l?O8I3 zxAUA?lq!&y)wCzDe%vp&BeoNzZJ}T9xiFxZmYR~}4ekFiE@4-q4&`nmYRL@eZVRDV z7w8_k8fVn&Hs?-{@N!B0;WCiZ z`+tuIF_b5dw0`EX-7A`9{{e&G@llo0K+q)Uby$Tqfng~T*c*JT;6pB3AsRPX{h7+g zqCE5d4g{0MuT1)){>coon*08N0=jYvd%^PjQ5A&a%WjfIMKih~yC>WZh%uFeS?y7V z)X$3S{`UVWps)Vwpq_tnf*|%YJ!VfSGrA>gb3K_)I1m(G5cA_7b?9^oHjsf%0=W?L zJbg95@A2QRGngA+Z9nj1deT`*Krt${P5Z|Obnb6d#uL&z7$9Bg%tvMQOLLe&NM03R z;;$Q_y0sncl<{dq&xm23-BEbSiLHvhk7^& zMs#K-Oo*RA!eyKP$gn2loA#@Fh4aMyYlDaLRazE?%ZzM|%glvLJ1z*P9oIc{Gg>8j zGrB|MTD(I1T3ny(@?3ZjfM+TKzssX=r zH{};DjfC`gU5Gf>>yIzbZ92ba&?E>&Ft!jKyK2>?7(px6u5*29y(S%oV|S^7LSV&THe_i|%kxAgHJ`Ou;L}=K7*H?~n-$YR6`k{Oo zS8Lpk9-HYYcdT|T6&wtn#Xc`a*ByvNG7~hRhB>_W!>{!VF<5WC<7^LyF09~_{kJ+6 zFgXvRRr+xxN6~GGO$#zZltTd#*sdLsWF&*V&U+P1_o_W8_^N9XdYBM^Qt*v>BW_$& zazKp;PaLdGkKGmd7(zBoS5pTZINCEga9x>D4`1@A8+Okcb@AXPW1`Ca92QqM{mQ3LroOQ=IR3{_xZ(uzc_-tRU2i+akw|QC> zh>&!hJs3{pL8_ZEm{IC@+j*z`2>$dg8=&2ufyX<#OE!lW6P1z!axL?(Jc8+XqX8!W zbY!$aRf5Ls!4utZG@ZMW#jpi(xI3=A1U@8@PGtV5*N6pj5N5*^LYa8|&48?}x~%RR z*jU}71M#T!EoJ-w8XPAeohdF#v{qCLaR< z4m5%rO|aT?MmG?mO6LVb;zg5>@A(qiy@y5rySd&tK@A*g4jw~x-XfGnvXdUtYMQB9?oa&EV3!4i9L}Rb_xv}0b0%x5bN+1t~>`8{f&xWp*23bBl)>`dH_rQJmn*axbLx6L?x-;R}Afg z4QwNYP7BQKUDT8s+Y$mW7z{Qhy0b?*bmhGifA1z+ST;^YO(J6m#g;0p0cFTxm+igo zuyl6Z6Dj53erl5;JUej8bIn%;49TnsKv?K1&t%Pmf-KaJ!0E#AvB<+NDPSPB(89Xva~{pxvE&ECvOT;YmHbjGL`*nKJh zj?HHS-8R@JzAQK5-FK}_(f%Xvv?5g&|_7foMpNl_nWL=V7rvk3YMTs^TDNp&Uy zM)iUTA431iYc7vl<=|>N#tnN^T{DOsz5zvX)D&Xj%rX4^7M)V*f8K?K#{2kQ0Y=^| z29mJL{l=&!qCTa$bqOQ3W=Fh6iuQmsRBtmu4qhZz{%GB7<0YO8+s90wJnndv#s>6L z{W2Tpr@AK~mIy)u>xd=KoeA!z`(4d$ULZwBT&d&b4y3eY`G5;ik>5fM21FF3&{OVI zx%39oB7j4+;3QPyuXyDCf+|Ts1sJ5Pe3V9PyJaJ%>fUwGmBB_-m!eF&skrGXIE_@; z%T*m`oAP*NX?uvxzSUCFf{BK6f?Hn_vWwl?B{xa4fhV%9W~)R8t`1&|GlWI?BFW#^ zs_idh_u>W}p>pcTlSP7?V4oqMLnh$?UVZNPkG(a=u>oUmaY&&>(OT}yaBzckXjHjR zwFh8oJ?^P606T{Ovz8a8_DMvX3dqob2^UP7AR=IdyUpqPvlBf*Rk0>%Y~d30AiYRj zf5nPhFcKQ2pGb5#Mt>BR9z=;75`cv@KuY{lbVMkSRuTe^0g3+adrefwwRPxIi9WJ+ zfO=@jxkm~3Hod7Xs>>e$S!Y~vBRTY{F}YMWgh-kCOEJRleQiOoELGel3V5w;LBJ|% zBL|PU(fMVi1im!TNT_#tRekUyf(ZK=%Q{etVN3CchZ)R0%J#%OOn{RTvnAfO4t=vi z$_{LR1z3Eq?Z!Z+>1gkD;cQ$pP|6yQbDc2}Dngm06J+4No#}?oi|sZnc$x7mr%Tr9 zL^VIFyPvid0o7Yt=O9*OFJ|wFVAS~+F2O+_T;$0&WYYunoad6MYx+H2E%C_xTZUYj zLDSZQ#y%XHD%S2vAR6B1j$!qxb|g{>!ZD;gHs0haUZ3ugV|w+^YHI(sCw(aage>AH zi2Xa{T8z5;L?xJrJ&|-nt3g#shy}@r-Rz&IjM5f-ios zs~!es=7nU&$wBK3h>1ILcPCfBOX~01&l}$D-Md`;vj2pb&}* zbr+2=Mo$!q%L10LV+_bf$8d7XOV)w|nyQf%1BlA* zM+2n7VO(7A&JU4>Zv;wKQGTvj12B=IuG%K~z##n6hP%MbtGJ1aFnUFgfNh82`sQH3 zC~CCx_JO9VkqRgSfJ`NES&}O^Q05hXv{AqjN9EZigmc(lTLuVl-*8$m0^~m?ISNtJSd%S|9r)(f4!-kwfKz%A72cPbuNLae)n6Ssnq$6xGvA%dGA zD{88FaZ>h%Gr_?C4jcjO8))roF!8!@vjGRq&>{3FZ0b+wxG^Gg-6*nA_@i0jaNQ_P z<$+2y-6)tV5ny}Gun8WsN%%@Z%-Bio?3b6|XvCrTZKW>dlj=wXd5FLfa#$2mV32`- zuMzOu=R6$5q+UUyElfzqoGTOsEV|b9^A1LZVQp270-?dCh=GD3r$F*?kl96xz(qnv2&&)l5O(bTW~h+Z25$x;x1-D;{~||7rMC^ zWFd@C=;JUejZOV%$DCV#VvIia+WCmk&SRE~4(#)jzPfA~G}|tCzxbHR$Sfm+*dnU^ zh?AA{4YBA-{p^`IHxycr=s*^`G0Rt!2&BF6}^&d#1t&o&Mw0C0&K$q`xC6^YOijxf3EJF8TRv@f`rCk7!%ziOx)GTM1M@^!& zJa(dX2$~XhayO(cMClz5j~ni_>J6#WYzTACuE$e4m!8eb|11I>y(U@~cWYj@9OR08 zeFwCQ4y)?GRWB*tK2U*zkFfVlj|;>+^S-k11)CPfH#!J=&%}1muYoHN{#3Lj1De6& z_qV>m_L&FfUyM}qdvnW+F|qf{k+cB54(Jd~@8an7E(d)pd9zbS&p1VOj%vEff;b3JRj<eXq=+DEmSHYj>9lKknGY6UyN%osz@LWpJOf}s= z&aeh?Vq4tUvl<}v+dYo{U*YgP{XW$+Pe#&wH!2d%KYDIf6b2Tzrasn2tsho2#;qKl zG;{fD;U^GkLnciyemuPd$_|a_M`lDwJ+O)4B~5${(6^HY-te9E6!cC)=9trCs=?{i z)S2}#+0#I__cR`3~9ZU&cI!ZXO=a&g=xfbHe zI#zY~zn318g}Qzl$ZF+ibK2eWE#J>a>64zB(IS#ZJW4+8pG8cgs%nTs;}k?tKQ4ux zh=d9l3!5Kw0(xDe&cNDV$scIPD4{V@RMd!4mm(xb^M6 zNRo?ETecR5b{Wafbj5=y8%)`YdvR|OitMb``-q{ejgrW>oRpc&j_B+!-uW)ok3p{( z641#!ps=Cn%kHYJ7kAuo7@5jt$4P(R;v%jyggIfS3wVZCU!v>EP;6u*bEe7=Xvz)~ zYv!+XQZT1z{Vvr%3TF{L-a5Bi??*s}M7UTt+%2F`O*z=Z>%P0sTZy)(QvBUkqCJ1@m)`r>7cn+@-_)vdwwNU9X zTaK{v+e!`Ybt9Z;a*@Y5NK4jlp867Ep=_50W)+8{vHd0&53IV07jtVh-nj}k5@)3! zthH9!by|+A{LD%!qX(PMz34k%WU;yl3SoaEeQS!lf6Ioud9Y}xCPs6K>+&msS8B%d zd#ff9zzv}>Y6V4EkOnS*w_1#h=1EJ=X)BdAo<<8UM>7VL{jyR--)n6R7FYhxmWXAC znC>p~;zz_{ls^IkIsC<)en}qBaECOKwsHMR=9^esAuWZHLUOodgLG%mnaVtiL0c61 z2Qb2!fExWr>j{nN@5f!q$imzN>;is${}emZ))_WlY@g)tjfxn{r-QRmi}v&vl~9b1 zqC{2WTg#+C%obO9=u95E=%>UtSg6%+7gc1LT1M(Tt&Qk+uCjbVdlDxb2~L+n$e6>m zr`PsY9pl;ZVg`l+>_DxjK4Pd#{EBk@i4l`u<2fPUqcXt`1&W0;1$swlR8d~>>E3gB z;lFkDklBOtyK*?4UO)@?&25zc0{@hg2w>LVR2;rXy@~zi?`EC}4az3;bI5qS|2#csmaUhs=K-eRMIYGJdK`IXL-AyOevvg); zTA4gOz&NS&bUK+1t}QXst5BQ{LB7y$l4r)a{-lhJwBFb|(nXgEV&`NJ=Gtp5NIwpr zgI1N>`qA=YXaPM_Y!$f}p~Ok=2t9T;DGW}eD++&FdRSPWe{9noQsP-GG3aDolWaH!gkY)3WrCk0Ea~+55+=I`n z&~0+2s=S{NRZGBCnoaviT%)!%g6#%C%}pi9{9{M##aDz6{ntuuuB26*LaLzf)q^do zQ9O8Z*;87|Eo8>G&zruNr=iQSYaFI`!oIM#-4UxN7Dn=F0nT_e!Pe}OVu$H5ARWMmO&`slry)lFeyHVSIS zA36l>lQsx^F$W6A6A?~_+pzGE+>e1jji{b|>)iOISef!D zU{d3SH{cxn254hE&mU&})6T7+klM!PPmTP}c>*N{7j94Wn7AP~usB4l zVR$Kt>uxdA4;8WxxEbsJ4;WCh-Mn@S48W_NNr6(&7ya9A)$Q%qtclcH=SqgU|FGwS#dI^_AA5u*B1%Wtg~k`+8Dt z;$Vi#hRA4N(hy=aWieZ|#QCcf1;w5*>bcdBWUc7!Pp^%jJWRXyrlyXY4{5~ox$^DLD!|USWDr}r z9si35#EQTIXUdBnvQ^ilG$lSg?Wt(?kA=0js~E;k@_;%mOJr6DG3>rNt3jR)XyQJX z;s-w^7VmsmgyL^l6GglgDBr(Rs?arN?v!j&tzLOz?Vt;QGs4~sY^FT9avUimw^N_! zO>BM<5=zZezxn8vR&QP528-x1L*#+2Z>2O|`N1V@EKr52C`pMhds4k$05-wpLER#l zjHGwwy(4Yj>Iyc83hq0^XA|K}UInKiPH}jYW(HLkJNbVrpDu6e)&c_`8d#m5*q`Nk zv1)4(@}IIqKkeQic~5CA5tLOG@bHUtn-Zlw4DVjmYPCsa!M85VuqUI86Ys@ z#Dy~hneTzkE ze}x`YP8l_iS1cWz`+o6L7bmy7hwN#Z9L=Wpnq^YEkUi}G>4Gbc9#>ecP!H45rX3$ggrrS~Qt z$aFkT(HUFG*bkqu{KS!mo4)K14CB>oGAW(aHdw5tZM?pGt5ECqquV?4t28oI*fobY zSsvwOU^6kx>8P8jc$B;tUO*2^xxXq@WmWz=+e)4%8&)_U95h%sK)oEv7+~qfE_2Ga zQ{IQ4YR~gUd?>9NCejhI5B;5k2d}|mj>;F?=9VA=csM-3zO*h&I!IUp)0y2uWO1(! zZBUQP{$59qiCC(sM_@do2*OUIkX56y)Ilb6#l%k>)2JL0%|&^QoKE{+%qr$W*eR?Q z7@{Gi?r3%@^3!rY(3gFZs?7DaU?T0)f7GVn`2qe5TK^Xhgg0%TeZ~WEpYcHdsXHR8 z{x31&+HI%95&~Y&4B@%T3JKKCr|zQ6HOVMYBMvQd$Bmduc=Vh z_$Ks)59%5N!*VjHlNIPZeb;Gby+gCpC2V<1@nsZA+qXPs*uO+VTUi~Vm0_?&xhtVl z`5GI{B`=Xi*gI_^?o2BB+kSA&n{)DF_c;{b+XH!A@?qd9xUnoc+-!$lB*CIMSh^c% zLCEOeD1m%4*v^{kB}Gy7gkdK+XF89ggNqd7xD(G$5P8#MFWiSdT(T(5I2|te!D0To z2z7IKGv|W~=$PB@-D-soIBO9?3BUD5SIz(OWKIqVwVp=gy@grg!J$*6A>^p84YCB& zm9rQSCf1quvHFmbAg!5`mCB~O`+ukG7Y~!HY?e9;?)Q$}ZQJjjyt*toEY%KtR2=Sq zO}xLxFZ%Rm26UY#RBmvqHmezR>uqEfUpnUpa$~m0mnA@7D+xcR?$5+HOzfOZ+D2%7 zrZ>HaSpr#dXBo_W{$IiYb)Jf$cwf@!bF|k+&0ADv_bxry3}HhcU#t^YTw0;CcHzr& zUl4f7R|xVdRt-K_81AKW&sduB#)U(%gr!LSNSM{$XVK(lMVM*&3=^j96H_r*YGB)S zZA9R^%Tk1cXz;MZkL>vq|5Q5&5G#|gu$eT4$-IBrNWi^sFS$NXQAwykdbgL}P#^!z z-eD8aU{cYSnQFgj!2*px?LhP3{-0Ss`JQ=xsA&sha7Eg2 zDZt~Ie;(8Un4&9e(H``w!1PUka0HotQXEJ6vCHt(r{IV81C@sZK{4!K#yHu6RSVOjQoj0MoI0ho|7LJb<3;hhZAzeE21T(|q^pEgbL>^IbqulnV`-Wz z+MWMk8VKg6UVR$dA!xBCnkenFjoV%@B9Bg|)8nn9s+m6Y$^QAeq@;n&URdy-e=JT# zYNTsfi0@t`J@$e*=8-VX^>{+1cq7)`aLGKH0g<;ptwFs7%f*cD=z`(j$)Wqyw14N9 ztbdr}`xwZ)+Dx#09x(FlHvdO0`?_hLu)ppsAFYIbtap-*iTngI;}oWAe_mtar=iQ@ zcpwvWi3)RcKKXXI;>AX|;##c`WgG1fWqKUr6rIm2>a$6j_^I0wIZ_-`^z6m>f93z{ zCQANxTTM(yWOj;T84y_cv5uuwB-tbe3E7xZ-Wq*q#WSs@$v^I|y$rcR&HVTV8Kd&bo7-?_*BXu-!H z5r*LYosPiVKCpNLY)(M!HJs);F{6h_GD^3ld!-U>oJhZv{X7UqD>G+XF~lh*I+{g4yyPv}ZiUammc##ht3dVlfbe zhCr$MB=Z2%6=FtV3jxh&Ny5gvWcLM#QGyK5G(Lj0?5@8?oH7*P)3J(9;`Hth{sD!H|%{M>fIz(R(}0_K8=8i01AyHLGQ_hbXa zzV7Chd)8Y%`0P1%-Pc6yxW}wkM=B2HX|sO(tX8WrFmK5lV z+QF~b@=%h&2nSwDXTdnck;>7v>*bNSve?r#R$?EM#)$sWqpje3y?z8Gye_XA65Zzz z)Z}qF0uISRXJ&HMOae=iylG3rHy)dGe?&e;M(^0iiW#Wf>f~dHdU|5Cmgi+wrYMCa zSGx0rEd_kiks}i)2u2(2`ZKfCBCe!Cv3eeZl1!cSgPsj)dz<4Dm%iT;_C9bD9U+4f z-vy*x51ks($0dUzBPe)~PiU0Uw{9~6Sr@y+KAMA%o#u4@K@eOOqoPIH2 zb$Ml@)f<-m($=s)i9c8;)1*?&p4UAQHaKdIES}H#M4eM1J!HHy;6`tM*hnGdvbMC< zvnXR=hTZpXxG{;c{|&cq)Q^Ocq+&_Z4?|9jXb;<;Q8KRe&y#zwsmkk*7$QYZqV~W{ zKJ0_zbhVS4C5xZ2gerfXD-PEsjGF-!v|Zke<)7-UCu4Qc4nFGl!JM8i*{`nxXg+(=t4A#xNh9HEy;4;!mQV*WF$UC9^bFQZQX4Rp>c;zhXekp`*+^K+nmuEe}9 zQu}k&E_vj0CmnS}yk;=cuuE+r&KMNoW<@>epuzezM~V@*l8$R5O&u|oG+u!F`Z)HpUJi9z zF7?k#(SWb<=C2!z_w`SDJ+cm;#SjbxE+5=OJV2}Ud>pa&|UmlCWQyama$oz*nbQSuweF(h0TF9Ke+Hq4y1Hs-$ z*PRNaBSIwgv`sipP?T!^3^IM9WCWzoo(Dr3HwH21V60OWfCpQ?1f7Zee~v+eLIB2+ zv;U}6&+hi!M8Riu_m#X~f1&aacX(zAMlzixdy~aHKCZd%H?RFI1KGujjh5Xl2)n5u z7^SH_ciPiJQUmitBg$?$AlXY6v%|@*p^qlWr@9|z#9VHqVRaWmgK#sUzONYJM1uo^ zt1S`zgNxn}_HYR1g;huX7H?9FBXR($RXR@s#2h>XwW1^2sU}OusA*_tz3C0l?Qh=z z@*}Eg1HW=>8@u`pd1n?m`Z!aYV5T3e$UhIz8rxEpKMtjl_!`>Jsf+eQ#nZ>=%2GH4 z@#9ZBqQ7wgn5kbo_WL^!+cP8lLPq4BuIP+ahf9-j+U(_w(`tJI{&MMn8KzWnQA>kkY==Ap051SkMS4!yN=spLjBLr44XO^ z#&hce`X(3O&hWtygyBmI%dH6qX8Yx1k@hav(n9UyP`%c!j*E5Hdb>KizF!~6{^I9~ zmZ6@LHHjsm#aCEt3upe+yGwRXMyIKxYT*|nFU$sV zB>&~dxO;PvPFj`U_yEi>uZ?}7*0{l$5y#G}J6tmhMUScVt2;XO)+sieptbW0Mn+v` zRueDhZjke*ot`e;+PeSQNpY$9Cue3;Sf#1#vKT{HL4TU8(%A)L@I=pkQB_#T@2h_g zO054G`~RjM?z@k^v+;W9t(;h=R1O16IYa;osO9W9^M45TIw z&2Cb21#|hQcUP?p04%OHOk8~qEKK>W|9GxWse)h1CAPHEJWO%4o+k#dJRzpc1J^;`*tZTGyI5lo$FZb__FdR^&&cz_@EwdEiEchqJ`KMk)m_c zLZtJojc{OKR4BD7yXB;B!PGd1#A8h{t4?IT%N#m7k=Q;Y(oQWz`E!~-IFX9st%<`9 zH$@HSoLND=glxUoSabqh_z+#)7bce;S(SzXQ&PDIKpBob)>d6YDHj@LFl=T=36s1%B#!C^I`@wz1XE3} z{ZfxvbRyhLN72UhbzZru24cvp2x?cn4X(8X-^Zhb;g*{AX5f#Kd5I_Sbc^bAC-xN_ zeVzWfsO)itaKaS@OM!PgdDD$!SB=6mlGwp0I6k{6sN9}YHrEMFeYfQVme%5b5GEpf zL-=z)`)Sk^++ZJlFNeucIAkDU*#F1eTL#3@Y-^)H&_Hmv;O;KL-95MkcXxM}0Kwe} z?k+)tySoRs;40TN+#U^`9}>KU=AP)BI>-T+?VBTtB@R(?JMcwT$-7lbem1y~d#RV$X`*e#^2l z-1B*1q|HUAi|ecs$EL=8JA?f!z5e7l>blgd-P%;YsjrIk6M5rmO@mv;P44(pU73|v z!C=YkHd#q!MV)z7TrD}* zGgMlQ4_PX&&vv&(BOB{HR_ryLO@(6(2!y0e15Kt!W5k!moTfH^cN6vkeX9nh4P$ngDQ`%S)fJDGs4kM3{_8~ zb(>Td>h8axO4nT!Us#a+O+zGLP``H5w;jX8MO_i?bY7(Mf2Jxq*%PPhr5_5ncVc7j;)Bg?|TyB7DHqFE^6!> zsi0yc3nh3Q$<=n|*PYdyHOntrS0VG`6RoelKFk*0?2Q^bZ#^m070K{8VyG?52ql4u zJwLU%r_OALHi3t6-#M)k=)@qkuhsZEt?6Re8+k6pT~-wv;C_tIFabx8 zc-Wz3GgCP7Zk?56(d?l;%=JK-9y$I{i0f?Hs5ER=@>W8s;lSgnw7QgS3Kt#COdw*& zlg?QI@JZpQaJ_hI!UZ)K$TJur@SKz%|#7Y^w zQDkS?v`uY>m%miR%7h0!LgRgiOhn63DXu<|vsb)9ees3STwc+aFu5qd2aZI=6Uhn& zmAt@|ny>oR#YrDW@Zq^e{GGYB-5OhQrgQs~c+7{XQmyf#MbA|>Jb;leVxPsCG;8tQ z<-p6(>k0K0#L)UHOo`#li@LE+l{?-IMu_0`{QRWw=u4J6n2ugg_GnIM&Iv4<1znUh zZ3&Bz*JMG#A?Zq7y+|W=>zrhAkEod1gy`!?*QPNa-7oGg^k6}zCi+puisJZ?8&T2w zHDU~S^K<(3Cqcy)Pnu(R-G!uh35)B zR_zU#vEt*=7;fdNn@JRo2|MvAaG;9-2fA6=eHCf0TN--3_F{`#Sn9o*!t#TpkW%Tu zI|W^Bvyz81>%pZ%jVNuKo>BdxZd%hkrc5tMI`n#DD`ov!F6^+d_l^rw&|F?l>SCXM z4fKdV2D$=pp#OUKr9*pYWEE6`zQu%eg08QyY>iCz>6jguPX+qN!V`5j=_!}`-+k3l z8dRTF-oyhOeGsps&ju2x<&C7W;Ef!07mZ@LkcevdP=U@9rcL``F0MAC5L8XQY#YE; z%n;1|#YX7zlR{;y2DUlx>nJZ0DZlT~2-QFWCh=8@7&yvD<$n5?rd9^VjG?=-*R$;x5>dggSjTeDedV~P;6M@Jy}T26FE7rB)@z=; zXQQY@zplAl2_B|~#GUBa{Frk>WzE`Ofg*6UIluV*N1ZlxG#>DZ_0i&qE{$!CE#AFE zm(8h;K5en&N^BhytV)|!XTE=nt^SOAvaO+jKsN0o;`>&OpJTA!^P}gZICWaD#<1DP z9x877A7*oPTXZtwE#}k zF`emrfF5+QNjsFAZNAG(h?4*)#ko$tyK8-(4PlO@B_5>Jug|xg8Ntxi_gVHEiM3;z z*;I6)YezE`K963H5uML&05>1pnDv*Fn`ywkoX^v!4M3n$rDHC3mv`Ux@o5RE@)?)D zgw|7NbBBG`ySW7LjNtR|a_jT*xL*4WaK5}V*JboBCTo0SZkZHju9Z@~=zvYNFDBV} zL3J9|rQ&Qzmg|h)oAk)`{?7aK0d=E}qG;k}Cad$|3OJwk>I`r{?gP-`nmOP7X5an} zFpu|uR_i<8Dp)}CP;uJrMTT`u?*7XG1PJ7mJM*zoQu8!b$%$?$N-^IXEgGbV!X=+P+XJ$`v( zMS0hL>s>}q!#1<-)6S146kV8BO4Rc}IsZ*f_8aufaxPD9IE;F$HMO21wJ5}A8KSJ4 z2}hNTc9kOg0h>MTCx@x)mBfD8Rr4Q#KjKSd^D{FSK>>>o14)g?p4VGrxNLZ3DEgJ@ z*Jp9%ZP6YbdmNsrqUBLIcHLZPL(SVVFxRDXeT)0!-=CZ2gI=8RLS(ePonMk4Pw6Wo z9BwkA1-(5okj|uZqilCn%Oi%g^RvD)_RgbM+Ktg}mkV{K@jUdL0rond_x5YGuFlp0 z_XwT$N9!-UOSSgL$LgbKohxyd0oFY5$L72Js-Q)w@{gu161(MD5bTv#-fCw>#GNf! zVyiaXQxpB^mRS%>x-a)m_p06k=9LShYwJV102|#TzCv~G8v{O{nVIS{oA|u+lZ+tU zx1;pUKiWK$?-C@NJB^8QrJmJNOuu)2s=;-Dbw~ANYfe(&8-J-xcFpqWp^Uq2^#0yk zdva$vijZoa<;u~@biJ7M7_nBpqGGsePxox=r?Q44otB4odAtSj9Q|aTImB&+EPGgFc^wsJGy{G<^En?WlXG zpU!;9h^a^}SJY)IKD)o4et7EjxAE~L-*n;gYvp@BWxoQvOv~+9D0A^XC|U(IFGDf) zdm-U|-S%0$eZF()6qfCSQk}bvOE-~s;2@7-yjfIPOJ?(ks1Uub{E&rnf4sB3U;jf6 zaKQcvmhSXDJCO+AepalC1L$WePlY>7DCroNuS4|q=ojq`(klf`45fnZxHufcW}UA9^DS!Jl;0zuZ*At z$T?(ApVOcC+#&bxPm{F(N(Xu>&!fg*PriXS=_P9R<{e5heMUaf&7H^i zadBa03gd~AsC#Wv%7_@)d@ zw*=oKzyfZ1@64y#ChD^W{aqi+!)|}oaA|Z&;n)dc;#pDzQxU9Z?Jc8A;U!yjr{3YP z;hhhlPl9)L&`hxVo6?2b-t&zA{U_qPSh}BdvYLD*X3}SkW*jDfq`ss1zU9u`n{|EO z@jj>e;QQ!n=NbKtrz=IT=xAe~8iY#fYNA7)5TyR20qyEtPwXh4ClYFYjm09&$RUOx zw^(?AXg}Xx3oCkuyC)ZV0}zyBI`$3;w(`4HA+tN6k*ipqwsisVv7>8D1mGXyqamjh z@$5EMh|xJ7VO=-!F!oqST?x`{bep&1@13+UH&j_swaF|{g$2AsqhsV~ybB^71ypw+cy(3xTVKrH#>w^{ z75$z5FX^8xr3B2?Uf?(+ELNS z!N$qnz{r8%wNDiPS*^c2=`VkIjDHvGKRL{M&FHshJz!(8Gc*3G)NcV30V@*|`#%KC z1gtEq|Keir#LYxm$zq;XM3_GujA)AOt+<>ZLVFlNyA;wpL_rdg*uXU;frEBp@h>{U z7X+bip(#Z2RVk-LV%M4l#!E9FGj(rY9{U@aMQ<`KDp=>0_G$;dL$eoU=ETO)_G4EB z3v34wbnBV%Lmij`F?e83eIVQhG{(CRa1e;lFp}2R#3c9-fO$AcYZMOjklyS9Kzylx zP|tawFBO~|1!Ij1Xo40Z2mu&>Z@vg=!M816XiS>>EMd5h*i}8B$|##v7-aF)NV@Uq z-#xj5yGFf}z6S>f_YotIr$DW`xkGqU+-)Zgj=gao0z!EBR!lz~chj^DE?W>fPhK1c z#P8C!W&o5yQ?uzLEEon3!Ubd|mLfkMr5|lWW%IGVYzYB443sw`EY}@EvIB z@@jh_c?qK+MC?L03xF!ygeFP7Y~&(@YkUl=*LS6v!m9#-MM*>h$@28(StWA)h8`g5 z_9Q#+$Jq@U`xY*&@9CUE4f-if`D1_68y!JJ2}*u(VM#b&e&HMtW#qSaPeRaljSP}(G# z1YTyF02{!28(_4WT!V^Sva|d_msr_BAji9m^h}Nq6kKQ%gm$A>kpdH)qZ(~$uX9*0 zT0CK(6ZJ*4kn(9?;aT&%^DgUM!UQxbE~0F z76m)t+`7}vP65O)#PKYB+1b5E+NbEr)!#o zk`G+NEjqE7JOs6s5{}pO(YIaEdEfKcjgPz&YHNAv47AO1E+*z^}YI_}y=>E(J&u1^+)Oqw=$oa&wUidN_Cn=Km$ zy(uoJ$46m>5Yca#GL;9E$zHbPbuijfpVubSP&q?(vmN#eNY9T**hh>$XRHN%!3b(S znx{(j{1lUEFU!FU+w372MO~u)b15-XoF0~S`h!xQyKNWbo$(hdXg=xo{82ihE`0es z>**j9yq2K(r^AI{I5UHgGcFGrT91ohl@4aza<&)@sv(C&Q|B@yuclXj$1i}8ae zog_w)_pWrP%Rzx}k7z&DP2{B}49`3>)YJL;6K%F^t!0F2^9YDb9%W2&d(+%Xp^?~v z5#iwFY@cM`<4kW>zLT*mjNp_%hXF<(1Z%hj_!bmTsTUsX13vuJN>rz6TCMkL`D2em zc~Iw5g&9X%Aoo>c0<@I5tK@VugO8d$&+SOJSUX7rL&o;twcO6tbfL(^GjzCMmwNM~ z(L$UWfs#fbMRaFqhHN zyvJe@!m>S*v(GB6-`JhR*12V>r>{Myo>!Lp6(zK^@P=!xy^Z_!4_wuaXNd1t=>^Lo zEqXMO{Kn|F)E>%g8V;TIBFuCzIu_{A_-RJ#2BsTz;8Bv$JV{BAL-cFfk$+4~y9#sn z(=+kZ^CHx=IsEL^cpuo)-cbhz-##z=@O<6B=JXktE{d{Z+u93`#6^Z%z}H>&Tt_{E zPE%>{CD=9Vjj}x5jkGew?Cjc;g-Mt*Em@94qi%g$g^@;?1g}`7>e9)zwJ)5>%%>tt zq3wvKFSUbX#%y&UOn`>?PkIy@g#l&ldO9N<*9ZkMiZdi`sJ}5U#S{>L; z%~nf22*cQEL8deur5K`L*UKSipb*sxO+^Nm&mJ>sqC9g@%XS>+uYOvwHjiBtoR-#= zj*NO~M{MO<$8eod!VaIm7tXU2e~C*SYfLxkrTJdOiur?cO!VhbEo9WYl!lv)uDPOa zfh*L43L;~;w1Qr(p~6Fq0q&9};?*o_U*jaE%qY0RV@`=JdkS-O;GCTpJB*4GHjD(=Xq?XnW#* zsdD)cBgCg~?mljhaF%{UqypnkaYsq^wl{8_x8&moe#RNcX@?k}`U4txExJS=g6dKb zPFK}7T93It?t!c9isY?sY4f*eACVtTeWu19k(?Ct3+$&omLS9t+LscUzB^^t`g|3D zBOMaR@suc!8$Alj{-LJO;C29l797%>(<(XkWgB&yDI-rU zF<-2Nn9{#+H9e{XDVXMXaE1n@fBt(nRS)2;tAm~O-UllTFv=$B`&!7V_rwxou85CpS5&62n09vqgL_MaEW3@P;1+JGd!t zE!j60YEDAG^@#Ewc_-e^jnpP~XEgM#cyo%5`P~}(PU^7&=;Q0s6>|2Q z+Q6xw0tWk-%Ea?aW_jDiN+5M8xoM;(rKjRd7Er(N>{Fmjp0j56{piqNK`C4->h<(< z&`d`jhVXzB?d1lG$_*=b>py;LjCEEK+Jb@FFdTF2Ji@S_I-oYMg0SBWvFLx)dN;i^ zC&1~s5#p2-p%A=~%n@s|236Elx_I@5JxqF+R@Jc2^JkFcqSEu}d11-n1 zBJF$A?M4WGmUtO@LR?BuABQ&cCV@OjGw9B~yywv(nj^PT0p_!2YG@IY-t9$;&~Js> z6a4Q#nGi;6_b(@;b_#XXxtcym_4uhP^p7QkB;gh&)nC@l1n%MV(|WZZ*>t7huvCs{ zX`J|1BXegz^4dP&Wye8{IPRE5A7pMtKtdjLfPaG^4sTXetc+`XkP@$t!A^gQ z-z#O2IkZ8e$lhN4(-zCqRCc_hMo}sEu-S24{A~5zw)d>k=`LaiX|VzbN((s!q&=GEe* zbyuWI70mqbiP`!9%qby=$4f6vv&}~T<5_?+#;%VLsr#Q#pJmPA-Kc*j7Ok*nY05IKwI6 zI$@E}d+t?XbB+bri-R84@4MX_&Fl(}b^2vn#_wbz&c$*D;L!0*_E0*b*R7@D-Kue!?381SfrA4<~-Wzp-g=m<>phY|L2d#1f_dUU)QLGeNu+j@6tI z@bpEe`cBuO5R;*(Luua!QQ#}`GZo=#z?~H88e&0r=X+c5_Qxu9*Z9nXlADB^Ap1|S zmv;K{L+3pSyk_AMT&=LzvzE-kNF%%0-k_kCm1}3h)rRlISm4JrzfXp<55|01O#2e& z1%na@js4Y!eNOgFG4fq7GHb2bm-phwXnx3VuaSSYLhw9KtYYxBDHI)21oO$k#G=R_ zVm}&wCobR*6h{qQxm>x-f2K=l(6NZyZ6` z&Rc! z$C|qT!{cjIK8q#V75sE$@AC2d44FIm>1EEbSgxI}cCLKGyFHvI zL$hDFp^ywWjSwC_|D;n9cCU`^XYN?*vFg6~5H?Ps3CYt=av0Zjn3H?xIA6%Eyd%1{ z3@OfejCbr2+M@i{-58ouUzPk_H(APoWBut}Uj%}5b(kU{=ljLOX?I9FcW(^KU(g_7nR0)JbI}>~gS+FuE62DdPDmR=9>)q)I=_<9ZXH7@v&f zobG;S&zf(g?$&+puCE5xi+1b4Fu5;Mb{$})P26y#CT%;O62LpGD_-2LBI&4>uBz4{ zboY5)*yiO~`&bYOi-fBV%6e3{UgGDTajok!oWx+FhF$rg^B2u;Dr6;-dVOGIXz)Il z4&l|6J%Za3F#(!SDU$4-BTX9B&qk2=X*tWEI316QgbZrB1TP%-QyW|E0#z|cZ!gRw zG)Q*1A|?9FcA${a*U)dRlGdZ=SMs5EEXK2T4d!>abSc5zOU6FLUo`A!?-!}pq{KFc z46axP;DHxe_y}T~t>Tf!kxqrpk-(0!REDOVz|)2wcSa+~Jx?7hI^KI8v{IUC(yaz| znK`;oU`@*6MNNoF+|b_1Rf)~0et8=Ez|$c=%yOf();>;lpPdJ`jdLS2?0N#By#IM| z-a31yI8nuPLFVJf^7TfaaN0vgu|?$`5ZU!D@pdu)8-1-1w1EM8Qa+Ir%akq-?eG_> z3^s0o1JbQMk+8k`Pp&9H_l6mFaOX_kF^i>3m^LyX#=5>cy+A2D0{xOk);|T&WCc3a&SaC zR*DrXS#MNsBCXU`i1t&T&pfnC&TILQ7J~|q&6+B(cLGJ{LrWYf8@|dnO_3^;!QtM$ zGg3hqi#Vrq@l(23qSv}`tOxTF!V3S6WQLuoGDouun8BmGuoEy*jO!`qE^~@9|4ty& z-(H#j+6l=NwM$<*X#P%PNNp6CTIQ+dl&7R?lR=8%$|} zZMPiL6-qR}`F5cGjjnmC5wCr^hk7msf^!tvxp`!{RhpsyLa2}4EKYklg4yC(RfgOD zZ5qm_^|<;^n98?loW)%EfdLx0X&6agYes1iX70LBk90}pM{1JR24jZ=TiptOc88N( zes`|HXPx_Kb3LhAvRu!e6@Vo^zM>l(iu<9%n6A#V%857p9@T*2nPo>iZ)N6hn*eHp$PrqG6y`cSWmo)y&fja5jSrl zab&&8F50V8*K<2v0{0Hd`-^O-uE}CtGexae8TQn#3yjxUv_8xU?QfOYPgEc@V-Qg< z7OB`{bayYwW=TrO!Af6F$YiYDDp_vH?j)5LcowI)y=>X0WgOke&q-I^4`|hwI_zwZ z$L*;tiC5l|_c}&=?@O^ek!xcXZx~bKByz%FG`y_9xI@>{wF0eJ!UaJnUJg1w^;XB~ z$aLJ!K3Fg#so}{)-t#lrSJk%W!n^*;E)JE^EX?i5{c)(PhDlb9vSc>3JH>FU(k*Y6 z>Z%+`#AxU#rSZs9g6P%)+e|YC9yzNc=M)t+<|cD6pgq7@VFWw9Zj=4p z($mbrf|KFTBE*kq*P=6oFSeF;HH~4iK$q5DQ;l+-+p^or3Ulv*d&-U+9nz<=z8FcD zz^e)ypecILw6=>1ci1u=|@y=Nf2w+-v?Sb6Nb-I$2NdJ&l=9QaVK&Ip_jJJj|v23-$x=@ze%yVYey zm^!JQ<-Ok@UfoOUgcic%OxOWnShzqKHX2)Z6if#uZyEd{qf#cNHG2vd= zu;8~@3HiCGFPTpndtL!PjoyWdV1BHoAKK?*bHdTsXS$&|^;v4Ileo^DX?Vc^J~{~i>T^E9E=Y(F z{R>Pn%sIa)SlEG0#)WvkE(%Hj0;&4w$N$P>{o+0TAhOt)IGF#Ue*Pr5Ua6ly_7r0L zi-r1MHx<&38@3K)fD^gB52q=LKsuZFjHt-ETJ=My;Mq4enNkr(K@a`O>%^7iOACmu ztS24{Ij75bMjBg6TS(u}szEz@wyqYf08q{@BYUFZf{YR8$8Z~F3wmUh+sBq~E6sLD zGQ{x+T3U@&!#(lZkgkY|Lp!u;(DAL`d*+)>W^~&DepoAHN|Oy!-J!cWbw0Dq-$_zq zGWxE~mSapGhwoN+zJKsGS&37pDbvA$Uee}gFe#DbL+&M?Fzt zt|Z>)3O}6D>MVRH*jr;c^gqH-2KqrS%5UEHcc)v_j8ny`Fbk#L3%~2~B%lWnLvoEQ zmrPbqJ)K=K%_G_>U&QrdO*)c6-Foa`$W`1uf$d3_{HMNq?T=|!x8yKR|Sn-I{enpq22JDB|`_?3{VnW3Yp!yk0%D`ooMy;xb<2>yJ( z>bxHRce;Pd*?w`i>`Y7qOu&9(W%yJ6D*5x6jqQ~%{_kGDb9s&V?|C!*O5^p+>k$Vt zknsGYgs*8azlLAW0rUMc?bq-xA{a>M+NfBYy>3EA!1#NQvVYLcuf_TEhx_L>=Pv^L zKVEbG$;+}ZvT^(??W(G6iK~j%;c_!@G+ix+o+N#+op+#`yXg}9taOA@`ilmqF z7hqeU!d4y3t!o@%)8c5z^vCXA9S&zM<2?H-4rduiiqMlpAB6>rhnj8g<<6eZMF@9d z5&JE7hv`f6_`LZjg|{^IqG;%{IB$n{J8=VFEC@BmGdUl>^xeLg`apT9UOYZb339C8 z=bc@@oeaD9>y*?fx#MzxGcSF?O*P1HJ5)^5^(sn-oqeZ z8%cF15^e$@+1?5D42AuC3lJZtNpL_s+u3B{z%Tu{TZenbCXgx zfPk)YoggO4_M|RqC|vMqJS9=2n|yO=yR;aQ-S#aBYl&VpfAGhvaO(t0PL9?Vv4*As zNwd8=lf^{!Yt%2pnIZ-qC8mLsp=YJVqHp)j5K<3Qa|uu1f7*+xDcTTbRi0720jr}p zwa{l7G38)*vwkZR=eFJ`kj9*|GGS%?r)+unr?()VBZD^1YV<`3J{jTzRiKumTp==6e;( zNDf-F;syia_6C$V#zyfPu^AzOB0eIt|=1QfuQnF*_ zBrWReDR+sfXelAEbh+=8OUy&spzRU^2UAFlx(8+GDL>cuTr8oM-(tzW!6Y47&auLq zDctt`sHureI?iM@JgyF5KL1%i!IT%5UOog;fzsR|YjYE%VOxpYf_Y*o2dx;zxMeGW z!Z6$4t)%CjGOPTIUp`u~lw8Ko+?}T;KK>Kfn?5k7qY^WlL*eqtx`{QEM!tOF$%D+x z+%Z`iogeZhnXsrvN5WQ5m3RR08v|10#TpE|OT==ezOBnT&f)$l8Oyl`H|WDmO;2Ga z7y$XC`@Z6;2`Dk4yQSeUshy#UBvz$L4EkpVNphaUIy{x!kEn*QWbn9XdJpWaEhML) z10|4?LI)qTzo)f+F4XCri>H^*iY6KM-bvYo$0JqaEz2Vr7IRGa)T|_CmRJ~_s>b^@ z4{kW-pg3ihj<%QwfP{O`fx8It=52IU-3IjVs%>DwJ25r*E~G?tj&6Iz=&C1hT>WOK z=qgX?(#fVwd+cdvh#@146^{J{;ut%HUDn;pgdXR~;u#62EXU=|(9!*w5GHwg2r8C2 z{+sx*&&N370o2x`LmApJ7WI5L4$!#jE%O`=;|=lltv>)}p6D^}tHp8$jNTKuf!Xy; zo;(ZB1`|WLD`L1rqMM!t6c<;${qb0Zn#)%jZXV}Crs&$Hkx@rp{e>rep5v@a`&M#s zKH~m%8~Q-(dYisbzBWy_7~QHsw6u!VvsDw#f{g+uO^@=eD|3vDhefa?g~KSO^zlzT z-VG)Wij!n3IX#PJm$P_Zgp?&p>!xpQGBE7vVvU=@WVT&UBdZ32O@R@;nO#eZ$@*Y+ z0=>0ez}dD7^?5X5s}U{K^K%)19Pb=#VivrB?js?FD5s+&C# z=6Gz`lgKtaZch=}lP+T$R6f_T`6;=~(2Ke)MnthMnbO%<@p_VbXpSMtJYFqh^L#bq z(CzEhqB+0Ofln+DcatB$Mr%mXrlZ3t(CH;BN%6oAfq?%JBOe|QC_gl~z1PC3w{x^q zw;8T{{9(()w(Y{{JzPPajW0&LV;8uhx}4tBQjOpd{uBJ>(*f@Z3WtS}4FV6$hZ3{* zy^LcTFm5U+O+{>1j#>_!Q;c0EKffZ%y^#hLhs>HOWwO3U9Svz6}D$5-5%_OB~!(l}K+-(mR z!Zni<4s4GBI;BTXc3LTlMO7_Fi2sTa8$u)W2xds1A|h23`|XwRqIGfz z1rJ?xAGOD$nX8OE0Ah^c*3qWR`14aEjSnDR)ukc{&@pFn(*{=Gd=VUEchGm;zNMDc ztn)I@?X|~qW(?D2T1iqFe)Saa!mt*Vk^CjrteSamBMvcOON!t7Xt4dGA?scv5%D_yS-h}G2b*?&MguhzsL7BdsW z-?R1@U+2XAo8`guYP|fOs`j5|gaOT#zqU}-j_tJWC4lR9$?V>C;`mlT9j>gM1oxa* z=mp|06FviC|2?*mK6YB=5TxPaww9NsvgZoX6Y_p+GCJYgy|1*z$^x=9%Oe3DUUuXC z$NFVDG8pv@J6&d4AzKC0cb9oA^CR5pnGCW&2|)Ne7?B?tdJc}b@a47>rx0jSnJ6Z2 zeRqbMFyeJa1qpCG@vDA%t&KxGDtVrUWm{bDRv8jcrACsn>67JrITYZYhkbijIgi)7 zbKIz(!%ZVuAc^>J@m;ETMMi;qrgE)(F{-)A{L$3;E7JQP9W|&*1xxo3EwFY+|FtcD zwa=ez$;A05`e6F6a|Zv_K>lU>zKVXM4`Bf@xsMXmGWw3DMppWcW(N8|2;%UUlBm9w znWY;6#lM}W{H^c{n*3pxzb5%7HemuzXryIiWo9Q}Vq@kYU}9lr``^$BJLkWn6DHPQ znf#4TUX%GH`rqLR6L6N`ubRsRV_{$HaLR*u(6lmFA}RpY!>-Z@WC zHrMYBPLG7QcA0AnQW9-E%L9rJw13>B`{h1X@9E?kFJxyj@9VToyl`HieYdQCJYA1ZM5rq@bAX{Be zG?gM-LcZc!O_?v2{Cvr$o5eo7!H%3_bOm`Am4Lg5i;{{6(M`17GQr}GZk%S?%@8S? zS6x2K&^omo_@fdq|9A`5 zT^J$UuPGXb?Hc4R^dmpaKgS|yGWhrtLDjb3&UE=Rdh1I z;Ci)4x9KLtA`_i8ZuWv(wGXOiHo&*zj_g<$xcL)%gFM1A`vz;uxNx9Z0|6vqD!ArZ zPUgD|9JV@<(I_({8J43AJnf_MvW0dFH5^~e`ry8;P^Id62Ef4m4(Dpu)U0_-j@kH5K z%r@xSFyvIEC5I6e2vIxD?xG%;Cl}it`)+omQ<)}|R4uicNeV5r`UB#L8?gsK52PBL zXw~+%>q_0@6g>?aN$D^nN}Pkt^?fSs0S9eUN-(B)>Yvfn;2%Z|jfSoxE~fOy>pzS{ z;o31G?fWKv9Hp54s6gKC%Y-nE>}R@5lRrQ7CBgY?SuaZ@HS-q*1Uf1J`j(vFDmKLfBL;vM_F&XznU0W(;l(BN{yWE7Y zlbqjU-l%cPu~s``Dx%~^n){{9F_0o-CHJ>{UA!>rOZti>Fxm6v(bI6D_1o8 z(rs6*Py%ru>Zn;D5b-$UW0bAKL{(b@l`FMKNy_!uy-`zIYY&)l);xuOK7J#%^PVkw zyLKXCZi!?PoQE1#VdCNJ2WeO{%KKexi@%U$8jZzqf}8GQF>P8-c69p+mh=#gZv_l}S?AV;Im)-)ft_H7&H*&K@g5QlF;IaQ-s~L?it#!=kNoQ-M2ywV%)P1^h9Z9|%YGutN^3rbKayRyI!;j!EnAtqs+ zv3=tY%_;{Pm|Bz)|AwO%7gA)8AsGlhdPcAyIl6Hn#eir7(*g5CPdbt>5P5VN!M};& zgo|b#0Ok0G3pvgh^6d(UMZg^3LDCK3-(+#ZNsj2-4t;z>ZI)>)aPlSU`nj(Brb1rVkT8END6}LC z28c+OBb(5e>MXE+Lt{+8)Bz^>7-DRh-cZx)brqHGC-{AzY zO$a;?J=fiQz#Afvq*cBNzXM$oFJ zWn@_YjFa7ME|bnTh+z5=00<8l*!_R7rC+U`zies7zgH9e2Ns6u{{agl-)i2=fOve3 zrl}JySAZt0UhD3RxH1OTP{4n25i9s}(n|pXZ@^QwDdfn8 zXH}q-)*wM7$34LzWEuN?@Ev=ttRVVD6)o02Zxcle{Ei7>k~B_wrWNjm0En+AdBjc+ zc+uG~(4iG;!TWSZ!{R_iy-BoI)r_^Dkp$|Vr&gO@YhxleE*JHrfx&iO!c3D}adiSy z`ht@xMyQSzk;s_Gd>uT>dCaOZp4v%AU%lLwBe6yEc;$z`D<3;fA7^@i-gX(2QmHqS z!Z#0Gfu4yJCOe~b>Q?qR%l?NZ{gV#>J|boU+Un+y3HqWll-1dOHR z<^UvdB&>~XUK4y}Xn+)x84y1Mk(~f=-JH=cJooxAU(Vjp$o`cwp^z{Hl0;^XZhxIp zaHx>1NRo}mH)t^PQ3EYq6eJ9M zRI8M&KR02M(3{l>HKA#>bMO$oZ(!gqFsIlA{Qg0D-&O-GpaaEjI!w-L&zjH1JpgO2 zvj&c{sSaMXa~iCj4=N)+ms(b88W=39kVhcCjgI$_!_WA0X?ipr`D?K>IdwVXcR|U6 zca4-miE#J<0w5T0d^bR0gE7H9AX;@-`Z#>(gp#Ad|DAK3)zN=w6;o?OvG7K1p4 zEH#6uzYw1Ly8CD7RgCFj1cJ-IfI8dzFM*t-dupAm`VnYqR-M?qdlQSYfK&|bbr?R`&@l~Ta}<^%YFX7{hS905xi!R7_nljtt}Sn z?Pc1QUC4QpFH##q>PR?U!zIK1`fSr1kWLJ0{!P8KkUcL~z+)Yrww?}M2LMzko35d= zz^?(nqV0iBo*kj;Yj@TlziKiyBAC5B)Elq&YF+Dfc$lE5dng+y+bAIZtG;?RI38>3 zF7sD*5CDNoK(GK^!%K&c*_{yKMhHsRJl+lrB+&-Gp3k}u`3n~4-MwFfJP!B-wk4Dk z=u`O%NC2Tve<}eu?%_L@-EB8a^4PQ6vr6hKD1?P>t#r>9*?BO|0Kc|N0t(!l5ORp7 z`xmm^9lbSvUF#QeL4Ipc1$_d(EME|LX*I+j00hDffFp=&y}$&3>TI1vMVo*jdm&5v zL*!W}7wDz?%Yg;KvW;%-8$ii;DTENeRTo0y{cy6eXSXLprYFY*ASU=33Xo(47*BuD zlzM5LT1S75dHdx2j717;o|*UNPvC6*)QL6W^77C@Zidf)?`>+zHg8XC>ICpKZIqd zHU>(G z^JmU9XH5DiFPM?}W<-qoiGEbS?;h_sRH^`?s~p5mWNTuwO4_w}50;f1%3P z^}NPkJ1ZQB57LX?jND81OkGZ9tm@5k*@(R?uz)EYL!R5s?7-ezqz!kZx0?2WmHxTQ zzI4w8uCUPRE^^}9rkRvUUD_nmnO&E`kyMSBh>WW0Mn!#T0M8EXB>uz00j&sd8lIhL z+7D=cN7Qai^x}b(4Kwv;Wa7c*p>?`{Jflc_$hwOto6y` z=6*j)sW&l&baIw81|!en(xor2dgmZ3@lcm8I;WxwiEs6Lw#^*AR(FZwNQYsV}lmBgMz;e1Kh65eGem%4qaVRMM&^y`AX{uL{e#tI~5 zJrha}4Xo2qW`=3JC_T>Iq7i`W)Wl^bP9k(3<*1S6mU*=*WU22>7Qp`n5*AE?UiSx_C@ls~pYw|>6!sIk|M%Yli{jAe5F5rD}wDMOv z35cK3Uc2vV_PMo(c`M7lMtB%z$LCr)#Pj?+#vfacU8C5%vOWvi3zG@z zoV?zctVX#aOSC;mtQv93*vC$gfjjk*VI$gfFl_t6AT$%;9d1uVVa;B7je5aTzCBc{ zXAcoAl-68ReUil(??RmQ9@ET@FiaGZ@Nkl6agPMNCUthV8B0A+CmE@bma)PJ1<1t1 zt8dt(==~@)q{|i7$HX`?h%Cudex9#U_lPlQdR)vT-K^U+hns$a8pvC;!pXJTSK=YD0rEoaEiez!@tDNeNaShKND(xpo2%b91HlWU99a0m;F5X)N z9;lIKa!dZzpkr<$U6kZ$!*D|qmpo2#EW8#9_fPW2S;XUxv006WRPI>XPqb~SnhZoE zoZvU&Atw}1%q3nA^;2rmXccVPYHozjayO`j1F>DdT)5(JbP(@wT9=v_z+T{S7^}cI z;bds8VU=@%k_Jd-dk@^lkz#>e1$7fplO1aj7%RR$dpfBX>gX_T$`(3_UL+C|DNrGt zTQ`*ie$GV%Vwu3+h33!*p=5$t4Zoxu~4a>~<>cj@=Cw6x)%+xtCKXpfI zVy-1cN}X=8wmLw*sAm1DRuy9i-*f-S?SAyxLurxVRoBz!4g48)Em*vK_sy4Mp*CM! zy=b9<10A2vUgJOVFa?L=uGCZ`JcKfmzGL`wF3Z8dPR4I}Ez2?`T{x{w+`xR6wcdfA!El$2j&7U< zIL^EUPJ@*?A?IXRqx$YqABL%%z5$FzvA!Rsc_B#RqqKly(W-F5rJiF4=iJ0ZIyQ-t zcxy1?lpnE)0?u0gCK_P4@R24GSM`qJEg;`E`$~^XKQL9W{iQ*eY#kebC#+ry)uH#W zosz}4vHZqtxum*x1s=dnaJLCSiYj^V4m{JEX7dfj7H4< zppQ~d^X^EdahvwWx+GVQ*++)RueYB9NtSEi)Cgl{;g`NZKe2m9)fj8E#e!+)aZ6&|a+u>2SJc$(Ow&No9L z$%VVl+zlbgOY@dhS>T8@Y^j%Th$7u4R^O~I%u=T~aM?P>Jd{P8n>w8L0{!jh>BhoR z6{N}Fpu#_)pC;nznR&@43QD`Sq>iw=WrxR{C8QG3D!koMn>Ym}@1126+I{DFxeHEo z2Iv-`#mj&Ygc%7lnp2|LvvAM2ez|P>I6IusQtK#ZO<*&GG3-QU+WP`3YB!0>(FyK~ zqAE%*|57Y-s&6)Du|kz8rFAfo5!0&x7j=w}x?+>N>sn^EU*t8B;gr^g1q}H2)0t)% zRk~wI4*`-AEI$Mw zqW3b6cW>gMpB`zkVrn}m6jl{4UQ|MU6Nzp_v29}O;*pKO7E0W9e)o~k8eOnEoy3gh zi!DRQC?+aHly-*Wn+(7m-jRl;J0fKMfc#j&Nlr~=->dVt+1xr5R%I?mI_Rg>G=S56pFu}dgTo$nPg=+9Zm6ZDIXPb6`}W3nb?&=cXBHpdOqafmJxHAW zgm%%+&>BV%W!MGJuXZgxkfu?yZi;SrSk~87FDG*-JDZ8H_(2i^f$6J)9Q+_VD)!Fo z*wWjt&&d(aDiKHC&@$V$IaRD6KHb23_;J=YOvxnxcIbO1UfI34quy>sz#Mf`Ca|q{ z))U%yZyu{fQ$}GhlJ|M{R47lGs?UD-G6Si&fe~1gL;Ku zhpG)W7UT4Krsg&;g{SOD(8X04Ff=An)C$w*bavsyNK}I#ce6ysFSD)Q%$1<8Wt)WX zuyw)Y*(|n6Dn7Jn*+8@H%ELs{R5N9Cooc-rbo{6%g-r+beld@aTNRMk^g=9jR> zmLZw!p$x~P9Fbw$`gMo_bJdW&>LkFo$4rf^lu27^Z^bmdwunYjXMF#*e&NGg4*{BE z*E_;{nS}G1+G*)_=yKlAM(~&)&I-Tep$<;4zPDFlD|5^~whupY-JYUu@MES#; zsA?{>EO-=P)NXMjbM?%ohQL}epr%;oX_BKFOHx_y@h|-vN-W9 z9i2XG6StZ?hgFwhX2PicxOai-G|mFWcdCiFCbmTPHY0OjU*lL3lY)=+j+A2R<|xFL zf{q&E{XWf29h}t5A;NPMA-5xKD2p%Lt`?v%Fce3M@~sS9iQ3_u3oCrqIg> zT`H3v?E$0BSEP3c$;{YS-376f{dWYU6}&k}n`v*blACY5R?pCP`w1j@^K-I7S>n4) z~9&=*=hy7f(Ui3xj~u9~aDFedhjkcRKuzFt+rN zsbdEI<+&5v1_`ViTKt8Fdefb0g9jy~_r4UvyYm)l-4FbN7O=Z;V?#w1C{7huEa}K+ zcAa|@`@sX9wHdf^d}IkIZg|GJU1=AJ?hO_}$&!1gj_`(7qvMHmx2hZEN05&)s%slm z^H3zVQmJf6UJ8&s1Uc7sGUX0)Yht@|*<${>A}ypwNpmHr_{#v>!=NJBrD5i*rP)Jr zYGPwBECC~mP+DrmhEF~@(%C8Vi^VPcEtcs)_P*i*Dz@iz#KMa}*iUr%ueGNdieTLb zXucJy-e*Z=M@kN*FLZsgan5{`G^p93T}d|h=^N4jP{(ls8Qo`Z-(a_U#1%P*VV~)) zfm6YIUcTrZW(!Ny6bB4~`M^*$T!U!i$YGveYj~)LNy_l^4@PW;;QTW;gS2BWnmiA0 zcaBqRX@hrC`L<=7gJPldWKCWo=b-nytJIv!i(78UI!Vq%8)%Pqh)P3vea1xHF=s?O z2OhC1mn>0!>W-VLn8Q%O{a`k2@G$euGDHo9LJpU8F{Fw6IP8_n`sfJ>w7#_yr0-+h zbiXsgdH%y*+ugx?(`S>BJ*TOmxjQI^@8NE4a08NPyYzHvnQv*n0xDCqij15_vtD|K zRXi{@&gk;=grq3(n0g7;-y1}V-PvIhMnKmQ@~lF|Qhs)Z-{}k7e0=v=5%g zgO9rnnQzkKVSIFPTx5%5cRN61>mOM$*7?ya=>e>mRtnKem1ob^>&BVE!l@j4A9KO1 z1+&4Kc-mq5FaYxW;BH8w>OAgIV`d#T3I?u97z30|0Q}79>O59&?JQUT7b=Wk%Wn^# z5F%`*59ve36mUH;5nio&LKDa(b3Z__AI@6mPkT|E>%`;4(E~7Izv9N{b|?tfJ?uqs zYYfQ7f3n_%pcrBktwCS*+KGa2<2M;C&>0=V4?eiRy)g=hF;C>Rq6;wg$!kMmJ<{h5dA$1;CD0xXBkg5wrCq|?_gjj^ zjm+tTqPXJ>-o(m4CuF?4>|3*?$LX`9qKvLC-_vj?^SZq)5G$^+ zvR0Conu64@W(&RGQU?#TaJ!|UA9yQiKn!DfW3kWH;^%axmqqgcUVLkBuT(DKaLW?k zD-*&fj|tNne1FgRw8JbuzQL>}N(#51G3N3tO2I=C(;+1)Gu3)fN(;S-MqR|1bbE00 zTS%F0G1t8@vR_}vr!}OlxFKerh}^u&i9A%>o+agk*g=B8gl0mIJomXP6t%?HT8i|0 zsRI>nNi$;wm$oMm4_OS^fMfe87lsZMoM0Pi;X5zGv8q0Hb8)^M+3R~_L9#W12-zvR z2q9}Tm;P}!mz5p414?bht209#bMryace87vQ62J*Q|V%YYd>+N_Kkn)2}htx0Z}$+UOeoJlcW61s2K zl=~~zGdtr)FGGi+jAa{+XL~L4q86RYS9pyJD6=J3{nPwRlA1=RE>x?cjV`)JU~KRc zJPk64IUuh2mw0fnjU-vdec-~=;kbt^c^kyR%!tFlrMKTqHvs}VXHB;;BWY2Ol1{Xm zy(M$@_ET7#+sJQnA5eRIx!w9ANo)$Wa@rx%$sV6&*Wl3aK#nkq4t-RcSEuIX~ z)VP|8Fh9ji48@W7**{f1dk-Sd^r;60@%9ZOQ^vhdCQR~9xig?9wI)saeTJpM0F+{{ z1x>)lrb)Qkg#V1)T{_B)ABhe6-89YO|T!`zZ|UU?UK?aoYv$DSPZ-vGs; zsmWz~^@BHf+ur&?Td> ztC!|Akt?Sz!9U3oji+-Ju$P-8go~cNHn4otC&SBq@{VuzG)(3dmS!LDlriUe40Q}& z`GsaSmQQEjX(FFvq9~mTh2B;JnQOy?Y)fB24urNxYH=RMWW`ZgxJ;3MfS>JSR$=2b zM5erFnz9`{*r+$P(ZJ-H^M_{1uPlO>m1~BYe1Ngv2!~r_|=K~ zij;j0zB-6QsST4pzx#G&^c-?XI7h|koDaLzlpwB+X&JWJo0&F%$26zR{NAHwpbI7e zo>!f+C(MR}cq-WV!vOu5OKsEGZDapPMAM|FEKnScEi~FL=6o~EuRi7#^j#=pRC;rwsyWIjEw)OMV-$O6KAeVjPfY3H z`{99RLkf~6vf37A%7g98!`@~dl&eW*x85>4jtI=zO z%QSu)1aBmm_VBqeS5<0R+64-lbGfF+jl>07(ZLiLbl9@Ah&*@Bzr82XFIZV; z2Q2u$foEJ;og{HUVEBTWq)KO-<)+umI)jD3A44_aelvH4N8jX@V^=n`cU>)ZT&X@^ zWk7aJ#lKxHvf7pON~*zmE+^{qXYVWm&ci;~y~t6cXj4O7Ma$A~Mnz$-3HkR|JLoLp6IxiE%3_SjK<R^2;0f(#Y4x%I>EAH?8Zl4HM;OY669_gHN|gZiHc9=7NoR*m}E7)V{M~ZFTGld zE6w=TK83B^yai4unmip(`#0V^xEvgkgQO-<=cJb|VR!vu4|B0SMDXZZ=Ni)w;Ro$2 zmx`-L@i)fQ;IHRalBylwdX(PmBhx?evNta5(~~}BVTI8f_t!Jx1gN{v4nitYiHv5h zMSf&)hsE17ADO-!P&=(s_Uv1gZ6QOkpSiNG@~(X8uyZY2IASnwKnPa0^i||uy_6d0 z(j0Fdvt4Rp9dJs#SH!L892^RkUxiR#Qf}M?h3pI|N}D+K9=Eig`jqMO$uK9p@(^1n z;i~7&^0EA>QLUuDv=BmEHI`XI=*w^8q4*6OD_q&bp^7$k2WM$+GB4+>!YP3wc88-eU8Ho4NDq?0eVLs9l2hkbYMSy4vZ`VjC*IUzC+|4!+7Eo z+?R&xht$Xn!Wa5tsn8W7t+R%@po5U9_w0 zb%PYQmD**;ly%_Edg5r()LbjU|GW#Td^>R7oR_?eVej-*)Fbq*FNat7fj1l(;qgm_ z9MahZGRqnD^gM%I+^4jXCC2t-QBHiWR>P8U2b_9pE_7t9MowY0pn}PsPkx;j z_FK_ip+-5pW6|hwU$NZJ;4 zAeI)c`D1p*-~-+=cul6M%a5tECEr8d+00&&Df5mA*u&;iyYLUw)k|gXQ{6cZ1l5=` zEMXvx6f_ELOkGl0q2TRA+h7pCe)j`wWf`t z+(=SqpI0Mdg{_T$Y1nqBRvGa^+x^USj5xFAt^Yv{*7E*+jrJ#whvu-4a$-qb?u3cy zE_vS+6<1M*dE*sw>DS$+I8-vY+Z0jVzDW^7?6aliq)JvDrSVKkMcbMjp1ozJm@!4Y zh1;;_fwIb7dG%0FXbl|e8j?;LY!OX$8x5l;3Yh<>oL;gNK`1Vr5`#R_pjFsgVFDtl zDSS|{jVD9)t)vNIC*e%KZnBCIwED>Q(zZ+ckU{O6@UF<=9i}1DWZfz3C19_64-Cs_ zVmz6J8t6sFpGg>eELi`Wj|o;wm)>QBW$7MCHA*4hx#SL6@b&z>tG|I#W%}ipoQEzzW$B6fe^9`p8A5-2`MQF16XOUUyt}^Ko?SkgR zx&{**nXB$Oo2f609^h5Bjz%K7tjD}pK{(;ZeJlj)ipzf9;-NV)Z>!Y$g6L;FgOhlt zx|4iDutqLwwLVEv^n@oN%++#UC7yJtBtM>miD(5?VZN?hXhPF7V1L9GC~;q&@~@}V zsc-mXMkyd1rEgi&v993$zJF1?zc&-2S7s;nZsUpW&Ax+|sCE%or0Q&WF*>?cyh;Fe z5s$rV>ZdPcShAJRUJDs?N*6)aBtXDk00O##7_LlX=>RZ=rLDmyndVMxTa$I}pT~nP z=COvh7RxAy!Og3w_1R%2{0fXh%2UT*dF>x$O-uzd#k^;9bR+A$w^C)wt!Xn|AXiK$ zN|p=m3kwWZgP~Q(^bdTopn#xonR(_gsplm?NMBn|C(rEEC7@uVaS0qqKsCa92&?o` zsXp`HuEszX&GYGLVo#oV1~@MwQs;0f}do?HPa`jVy*w(74zbD@4!8&@3cbWujm8n#@ z;m1#^J2cpBFs~(2@c zVs!p^^c{ryS%szcE0WCz$Ug4GfrD^|{x0tI(A>NC4-n~%x4q7}RIu;!{4k>f^0vIt zG9sbcVIZQ~2l~i9N~W&?HycH%4tx6j^gT3H1|qP{wB!^EXYk@*lELA?&{W5|mQcs| z%w?6i%4227H|`=UxHWxb6~wo7*FCOHvi4>7MDPu3;rp;Ro$H?uc;VQ|laPiDGh2C+ zp%GN=#30-^k(PSXR%RJqiTMOxQzN_4%QJPGY*qF+L_!>);x@i}WI^{jIh#@C#SzEy z)8_N*C$J|u=ot<_K|6>690*s4>$a~)(Af7^;yTY-pTMB?MezTH``jOnQ6M)VPG$h} zAO3S7qW5o+QY?Q$!2iogspu~@?I13Mt2bn~6ZE}Q8hn5s$!nf-zlMir;#v65tI2v! zS6WJ`UQmw+ZlSoBiYm^eL0UmTM@w= ze5ip)<7Uwu1x9&yt)bq(_+n~YdQJ9>$ny%ei4E5CcQ2iYhz~V0Nv;FZQP9hT4~;MhDcW*sBAc4Y9$rDIGM54&QzfsGW=Jj0 zG{P^X#Z_|pgFF8gsoL?n+1RJ-)u7Q0`8%?wjmKSl{a=j*0D(_ev@ zS$Zd!x!zHyMdg*8kiV0Y)QY>gFCOqE==R$x)Iv?wiJ)Q^Q8luKR#U|eEM-i zh}6&-d$^&(OuZ@BgwzSm!|`$^e>}WJx$8r06KibuIQ3GeFYE)<$bLkQtem~=q-ZPq z$I$#(G$q#nmwUup)}2yjNFV0mDt>U=e)SdWmP8hET=P=3eXR)(=ehPr*5WdrGgH2L zj>CSyknA_tMy+_$uUEs=VdA9PhssS_{QIZ(87;!GZ_YIN3X4t$+YN%1N?%mpx)`fG zlKDPt4SLmBowrf42i+O?D3^OrYalO3B^Gv~=Tv@@n_bXZhhZyGLBEqeJ&3-HAT_

uU?Ft;l=4ul-WLJIvQ`IMf7p$%b(_dDEj~7_OLm;m*Zv z(06-T-3MzNYiVI;an;q6vEG%u(H}S;kHHw5AM1N_$8V~)R~nnFMYxUCL!jj*w<}18 z&Pe&xSKk$%(RlOvG1vBS{KK&j&c;51p6Vh*de(z^E5GA)YM&I}P+vKBHFw>r5VmqrZN+JNx{=#*_G_xz(I>_7 zmgf(=IgjG}c#e(3jHiUBgzFb#Cl?Gl)7RC;*92Si;YI^j^)-Y~3g0fDBCVdDK5nzS zrJHUmBo4UF&$p@b%>&L>eU>l$=SOPk7lsH)UNgUzaA}oR7lPN-XC-YPYdQWv%`N2c z&e&As=nZvz!bB+;UF(>+HC=W=j7AL3Tq_+-N7=eUuQT<5v@%TtHrA{jg`i` zfsTcwlFEji-siKXwy;OFn-q;9PVYE?!F%;H-8dwtx_Ne~+iVgMT*Y!tkO5L1;*0ay z3Mmt^HT~+i9fya<4mI?mZ_wHRK72pRJU0utL{f+k!`qH5YbpVA^S(EbTQx}e(0pvkeQ>Z zvM^K7fMapCk%C~3>T0a$Dia;Ttau?A+Q;SX*0Sr?4OL{K!TT{iqR+$ItpjgK%?v)@ z%+_vRZ3nTC6){QHnjVu8Be5EVMn2AZAwMh9!b@W5ef+m7GMT`bK(lpDItH1Ss;ux@ z9Uq1G5d;^NLj5HIR-TYGk97;g={~fpzYDO!iMbeb7+S2h7)frA>+n&UvapYiil^#< zmA2~k58qubDO|6ET02QQEczzC(JNtwKQ<9GgG-`-Gl4H#0pm_lujtA{mN|mY*sId8 zPFk<&dQT?Ypd{uHegq-mBk0;ESu&Yi6u}cfs{5YQQeU7IC58r%xE?VgYGMVt*Mvi; z0x}}X%Ab^Mx*a9vB}f$Bu;aW>U1T`NrG>rmc6AHS(})y9MI>=G2;T6^&gTi@lOQL{ ze4Fs0{9#adR~ECTCioKBMcKwzHuk!iv7&|LYDuj0l1%HBZyoJj<=zSDadR5uT1Of@ zBYDo;scy*8Ji&FVJEuNb7D|f>H*o}7voGThpMpbRU%gxM%h8jeC!q*1GbQoocC071 zhSr+xa1TF2oUcwl(<9wF?m!uPzt^4lDtc4Qhy-i={m$T5)#)Z16xx|6of(K%5~*yC zPeGh{(2>xC)WsQ-s%eOI&k&^sjV<5JXm-83>&O+(wF{vQ;4P7ijK&qOHOLlfWy1#` zMM!=n*Cv}Inv!jYe-{MTmJWHDkSs!rXRLFev#}H?bjoXlfk~Z(ISu1Q0Vh7wRA07{ zN(}`@-#lf2Sf!JlJ8zxv+zOZYMKy6sm}+6S3=X9H0M)#0mk1*R?`ho(RZn7P+@02k zLv4>@OiYXVKKha{vCS_zJAUPLK(@s%NC9GjMR@sAK(pd2ZK5{ zj9VM=ExhgV=2G)u#lwF!rBxlm*F}vJdh!yzr2qp>J`Ci&^Gpa4lZ~>(=AbzG6>Rtn zj|9=L6g99_e+`*KDDCjsFmm7<89XrRHDv7Z%`L^0@$8;rfpvyG$P_yeZGw#o=e0}X zjJTxztOgd!9q4^w*>FIu&>|c1D%}Z1A+kQ`^YE_2t6+~b^hB;rlB^6n6a~3G$4wpM zlT-yB5^%J!@iZQpE6F&FNB-7Fr>mnIQwvXilR~hgi*MacEsx;6R}7h2>2f|mEO^#N zi0XxhcoC$o8JsBk?Yt(7^Uu=wF&%LEc*1DuMeb#+90-Q zzuo<8jLWd2Ru{>wqj1QQKm6e-FceC)VRK>6nCYdbE^T&3n@yLIJ zTuLc~{Z{Y|qxoCNc}BZGTkA~Ab60u3um56Q|3^p{XsN`^%=u5?9?Ku`N`Hy*08tD6 z+XxEzew%h8GHwS;f5a^k@nqGF6 zJ`}sMJ@<{&K}c#<&m(2<35V71=;d;;64Z)6(z%(W179a2s3Bin;9gHgZsn({&#wD= zl=Tqg`Qq653IL`q$152L6bXCp+Bk_)ac-957FTfbd^9e$e-9juC0Y2E`k0b)uAtmWbI8+(;p zr=6Uun6#Peo*LMZ~ ziTD7qApTs;{x=LKB_&Z=5t08%;Mf13T#6ry`oFI7{~BBLr(S+qYySn8;y>x=-&VT~nty?m^BYn0y z7I%jWD@hV8{h}oOnXILde;as=?uV%al8+&YRH&$LgDB6#7qfsUB;O#QJA_wtiIDuE zY~lrs+mBOUbIY#1QC{u$X*uR@;oSa`?r7qVF60!Vq95`&r>3jaUDA??ETpQJuB};b zw`)r4%LDH%6I8MZ*RQM9aI@X-+cv{`-{Ngwu6ca7T390~~IndPPsi45TL zI2p54^PoF{+TqcCIxVa&P)FPT7&tPXSfHkweb5Z+jnkl6-&lA0cyOXC@Pv783k&W8 zN@Sv|KWo@}Iz$|2AjK8!*^OaBtc$dmqcn2*Bby%m1fhiC1@b`T}-Pp&gSQ_)L}0Th4S?YS+4ESnvx+d;tnoC z9t;Jy^yRf~a=k)4Cv`mr4wuUc(W08(JA9Z}Kp5rb7w=Uv^`5gG&mhuP3}@lPnce$B z9gSL@QmJRkOa_LUp`;n;tFwbCVWGIH=m%fC7h6&j2PNS{1<|H zu>(+k!dYvKF9iLhive;ZTK>WbYA3x8+EQE$3rG?HkWyw)KKUL-^mQTPnxRT_?RLon z+Amb)Se3z9#P5o%4405-T`r7PxSx0B%H3d>aelZgfUE#IyWKcf$)zQbDR(2w#z zl$VKA1oZZTsC2?0DqVl!qC^m6C(DQvHh`S26EXMPLp9*{hLY5A!Sib3Kz2b-09yUvv+>^|6I5SWNcZOW3-b*IP^@5`{ zHWPz{yw^mWVNyx6A%x$lbn^oC?z7rmUky`dY{Jh-cRXt#f{ZX9hub`S^9dAMC4!zq z&7Ml&y-rl@s>U55K^42HlJkxz?3?fpkTzT$Q6Wz7&1=yO(OGQ>2YApW=k>eanP;0L zmm_Sz0kcVa9kAU^>b1lv`n#Ru^lsiq&?S&^NmW!iPkcR<#LrRpkrq#<%N%>_aZ30# zIBiA4hM7za@V{z&Z zJ5`KHx?}0=Zz7+V_6H1T?>I{Bm-N3xyy64~T_CMMRO@W=j7rhS*|UW z(bQ~xfY^d~rO=b=?yx(dJuGV8foubHPG55F|}W?5UX5LbG+%z zcQZlUtaJ8B?W+9Ft3;ZUgR5qe8?DOQi4!OF#1qfA%bfB9ez$55`Foy;xgStEsE~8I zOM@GfOPsW72fH^J;2D!uglin8#bJ#FI15Ut?pz=1@pm~^LOoswVLjV=M)po_&lyLH z4-#36Zx9b_g?tO}84&oD2Lz%s8d7+p1Y$a*aKaCPM5akB2;Vfe0N;|^iS!WPRSujy z8v;Qf+_wPNc#ha!8JykEbW9N6MRCFx4Ct4A*j-Ugu`C)n6S%;Zc~0N;0I>-jM=leq zbpQu~(uFb^Tw?uX&{c^(pwcoqbBakfp_Aq|XYGE&yiid!*pmuNy4?hE>nfIZq_If??h&p9D|_+YU*>4F=-+%HY~dU$y00htkP;}JBBdkb`_ibV~>SnQil3eK@FSv2^Lb_=#JF)$GSw|HT!!hrnLuRBH z)}gJ9X&~(dNwuRHLGO0*v?;{A$eP*UL`@HzL3QxF~`_*`-%YaAx)g={E~2)brP6rD|4Y4XG#3N zK*EHK`({RdI;upLUFL)iasfOwoyxSBelm(YQfspVq`* zAP_({P7uWV=lgFEh#%>X{{YbayMBL{|4M=Q4;}vVoan!#K>SeqZxo0h)Bdlc^L@(y z{U*OBV|M&usl75c`p%R10sU6|2kKDR))|D#VEPN}0{F+?!_Q>YjKU7a`XC-W5tW*) zm5GDCHPK%vi=fj#1||jo1Axc`R6^R^#@UU?jf0JdiVa9_VD3cy9m(imXa-`d&=D~M zm{_TQA{pf!Y>k`^jX`z({QE}}=C(E>-*aG7feNw#*qK?cVoQ$o} ze=71PED-qDeUhJJC9DrZ&fA*)P{a?F{39#%pBK3~7@NQXh=6RcfPcP-*jQOXrkV-S zPZ=u{69?!E`1^~<<{vU9P7XHEZo=PWpg+FnLHL^tRPOr|_cs|cJMcTA{%d%~<6Zt< zj)nR6`dFC1Bl3SK$NHOISXfv&e;p4N77q5`WI$#P&_=@F>I1PrnSSjXNcNkqfIv`( zf4vq6_^ls6ARF-4dj_%tzHcV{y+jDC3&(Hw0($;_eGY&efZy6+1u%WzH2GT_ zAQ|hg{RPcf_Fu;hB;)+`{#XI*?7zvFm|1_90a^cYe@+hipsWrKKQhfLnR^(2pR3;w ulND`korr$S`JZ3s-x)%nljR?qCqK@X^&Omk%x>22-zl)<T8vEy|k! literal 0 HcmV?d00001 diff --git a/Chapter5/chapter5.tex b/Chapter5/chapter5.tex index ec8210f..b7eadc8 100644 --- a/Chapter5/chapter5.tex +++ b/Chapter5/chapter5.tex @@ -83,7 +83,9 @@ non-Hermitian Hamiltonian thus extending the notion of quantum Zeno dynamics into the realm of non-Hermitian quantum mechanics joining the two paradigms. -\section{Large-Scale Dynamics due to Weak Measurement} +\section{Quantum Measurement Induced Dynamics} + +\subsection{Large-Scale Dynamics due to Weak Measurement} We start by considering the weak measurement limit when photon scattering does not occur frequently compared to the tunnelling rate @@ -542,7 +544,7 @@ described by this model, is that the width does in fact stay roughly constant. It is only in the later stages when the oscillations reach maximal amplitude that the width becomes visibly reduced. -\section{Three-Way Competition} +\subsection{Three-Way Competition} Now it is time to turn on the inter-atomic interactions, $U/J^\mathrm{cl} \ne 0$. As a result the atomic dynamics will change @@ -787,8 +789,6 @@ into the strong measurement regime in Fig. \ref{fig:squeezing} as the $U$-dependence flattens out with increasing measurement strength as the $\gamma/U \gg 1$ regime is reached. -\section{Quantum Zeno Dynamics} - \subsection{Emergent Long-Range Correlated Tunnelling} When $\gamma \rightarrow \infty$ the measurement becomes @@ -1085,10 +1085,716 @@ which as we prove can successfully compete with other short-range processes in many-body systems. This opens promising opportunities for future research. -\subsection{Non-Hermitian Dynamics in the Quantum Zeno Limit} +\section{Non-Hermitian Dynamics in the Quantum Zeno Limit} +In the previous section we provided a rather high-level analysis of +the strong measurement limit in our quantum gas model. We showed that +global measurement in the strong, but not projective, limit leads to +correlated tunnelling events which can be highly delocalised. Multiple +examples for different optical geometries and measurement operators +demonstrated the incredible felixbility and potential in engineering +dynamics for ultracold gases in an optical lattice. We also claimed +that the behaviour of the system is described by the Hamiltonian given +in Eq. \eqref{eq:hz}. Having developed a physical and intuitive +understanding of the dynamics in the quantum Zeno limit we will now +provide a more rigorous, low-level and fundamental understanding of +the process. +\subsection{Suppression of Coherences in the Density Matrix} -\subsection{Steady-State of the Non-Hermitian Hamiltonian} +At this point we deviate from the quantum trajectory approach and we +resort to a master equation as introduced in section +\ref{sec:master}. We do this, because we have seen that the emergent +long-range correlated tunnelling is a feature of all trajectories and +mostly depends on the geometry of the measurement. Therefore, a +general approach starting from an unconditioned state should be able +to reveal these features. However, we will later make use of the fact +that we are in possession of a measurement record and obtain a +conditioned state. Furthermore, we first consider the most general +case of an open system subject to a quantum measurement and only limit +ourselves to the quantum gas model later on. This demonstrates that +the dynamics we observed in the previous section are a feature of +measurement rather than our specific model. -\section{Conclusions} \ No newline at end of file +As introduced in section \ref{sec:master} we consider a state +described by the density matrix $\hat{\rho}$ whose isolated behaviour +is described by the Hamiltonian $\H_0$ and when measured the jump +operator $\c$ is applied to the state at each detection +\cite{MeasurementControl}. The master equation describing its time +evolution when we ignore the measurement outcomes is given by +\begin{equation} + \dot{\hat{\rho}} = -i [ \H_0 , \hat{\rho} ] + \c \hat{\rho} \cd - \frac{1}{2}( + \cd\c \hat{\rho} + \hat{\rho} \cd\c ). +\end{equation} +We also define $\c = \lambda \op$ and $\H_0 = \nu \h$. The exact +definition of $\lambda$ and $\nu$ is not so important as long as these +coefficients can be considered to be some measure of the relative size +of these operators. They would have to be determined on a case-by-case +basis, because the operators $\c$ and $\H_0$ may be unbounded. If +these operators are bounded, one can simply define them such that +$||\op|| \sim O(1)$ and $||\h|| \sim O(1)$. If they are unbounded, one +possible approach would be to identify the relevant subspace of which +dynamics we are interested in and scale the operators such that the +eigenvalues of $\op$ and $\h$ in this subspace are $\sim O(1)$. + +We will once again use projectors $P_m$ which have no effect on states +within a degenerate subspace of $\c$ ($\op$) with eigenvalue $c_m$ +($o_m$), but annihilate everything else. For convenience we will also +use the following definition $\hat{\rho}_{mn} = P_m \hat{\rho} P_n$. +Note that these are submatrices of the density matrix, which in +general are not single matrix elements. Therefore, we can write the +master equation that describes this open system as a set of equations +\begin{equation} +\label{eq:master} + \dot{\hat{\rho}}_{mn} = -i K P_m \left[ \h \sum_r \hat{\rho}_{rn} + - \sum_r \hat{\rho}_{mr} \h \right] P_n + \lambda^2 \left[ o_m + o_n^* - \frac{1}{2} \left( |o_m|^2 + |o_n|^2 \right) \right] \hat{\rho}_{mn}, +\end{equation} +where the first term describes coherent evolution whereas the second +term causes dissipation. + +First, note that for the density submatrices for which $m = n$, +$\hat{\rho}_{mm}$, the dissipative term vanishes. This means that +these submatrices are subject to coherent evolution only and do not +experience losses and they are thus decoherence free subspaces. It is +crucial to note that these submatrices are simply the density matrices +of the individual degenerate Zeno subspaces. Interestingly, any state +that consists only of these decoherence free subspaces, i.e.~ +$\hat{\rho} = \sum_m \hat{\rho}_{mm}$, and that commutes with the +Hamiltonian, $[\hat{\rho}, \hat{H}_0] = 0$, will be a steady state. +This can be seen by substituting this ansatz into +Eq. \eqref{eq:master} which yields $\dot{\hat{\rho}}_{mn} = 0$ for all +$m$ and $n$. These states can be prepared dissipatively using known +techniques \cite{diehl2008}, but it is not required that the state be +a dark state of the dissipative operator as is usually the case. + +Second, we consider a large detection rate, $\lambda^2 \gg \nu$, for +which the coherences, i.e.~ the density submatrices $\hat{\rho}_{mn}$ +for which $m \ne n$, will be heavily suppressed by dissipation. We can +adiabatically eliminate these cross-terms by setting +$\dot{\hat{\rho}}_{mn} = 0$, to get +\begin{equation} +\label{eq:intermediate} +\hat{\rho}_{mn} = \frac{\nu}{\lambda^2} \frac{i P_m \left[ \h \sum_r \hat{\rho}_{rn} - \sum_r \hat{\rho}_{mr} \h \right] P_n } {o_m o_n^* - \frac{1}{2} \left( |o_m|^2 + |o_n|^2 \right)} +\end{equation} +which tells us that they are of order $\nu/\lambda^2 \ll +1$. Therefore, the resulting density matrix will be given by +$\hat{\rho} \approx \sum_m \hat{\rho}_{mm}$ which consists solely of +the individual Zeno subspace density matrices. One can easily recover +the projective Zeno limit by considering $\lambda \rightarrow \infty$ +when all the subspaces completely decouple. This is exactly the +$\gamma \rightarrow \infty$ limit discussed in the previous +section. However, we have seen that it is crucial we only consider, +$\lambda^2 \gg \nu$, but not infinite. If the subspaces do not +decouple completely, then transitions within a single subspace can +occur via other subspaces in a manner similar to Raman transitions. In +Raman transitions population is transferred between two states via a +third, virtual, state that remains empty throughout the process. By +avoiding the infinitely projective Zeno limit we open the option for +such processes to happen in our system where transitions within a +single Zeno subspace occur via a second, different, Zeno subspace even +though the occupation of the intermediate states will remian +negligible at all times. + +A single quantum trajectory results in a pure state as opposed to the +density matrix and in general, there are many density matrices that +have non-zero and non-negligible $m = n$ submatrices, +$\hat{\rho}_{mm}$, even when the coherences are small. They correspond +to a mixed states containing many Zeno subspaces and it is not clear +what the pure states that make up these density matrices are. However, +we note that for a single pure state the density matrix can consist of +only a single diagonal submatrix $\hat{\rho}_{mm}$. To understand +this, consider the state $| \Phi \rangle$ and take it to span exactly +two distinct subspaces $P_a$ and $P_b$ ($a \ne b$). This wavefunction +can thus be written as +$| \Phi \rangle = P_a | \Phi \rangle + P_b | \Phi \rangle$. The +corresponding density matrix is given by +\begin{equation} + \hat{\rho}_\Psi = P_a | \Phi \rangle \langle \Phi | P_a + P_a | \Phi + \rangle \langle \Phi | P_b + P_b | \Phi \rangle \langle \Phi | P_a + + P_b | \Phi \rangle \langle \Phi | P_b. +\end{equation} +If the wavefunction has significant components in both subspaces then +in general the density matrix will not have negligible coherences, +$\hat{\rho}_{ab} = P_a | \Phi \rangle \langle \Phi | P_b$. A density +matrix with just diagonal components must be in either subspace $a$, +$| \Phi \rangle = P_a | \Phi \rangle$, or in subspace $b$, +$| \Phi \rangle = P_b | \Phi \rangle$. Therefore, a density matrix of +the form $\hat{\rho} = \sum_m \hat{\rho}_{mm}$ without any cross-terms +between different Zeno subspaces can only be composed of pure states +that each lie predominantly within a single subspace. However, because +we will not be dealing with the projective limit, the wavefunction +will in general not be entirely confined to a single Zeno subspace. We +have seen that the coherences are of order $\nu/\lambda^2$. This would +require the wavefunction components to satisfy +$P_a | \Phi \rangle \approx O(1)$ and +$P_b | \Phi \rangle \approx O(\nu/\lambda^2)$ (or vice-versa). This in +turn implies that the population of the states outside of the dominant +subspace (and thus the submatrix $\hat{\rho}_{bb}$) will be of order +$\langle \Phi | P_b^2 | \Phi \rangle \approx +O(\nu^2/\lambda^4)$. Therefore, these pure states, even though they +span multiple Zeno subspaces, cannot exist in a meaningful coherent +superposition in this limit. This means that a density matrix that +spans multiple Zeno subspaces has only classical uncertainty about +which subspace is currently occupied as opposed to the uncertainty due +to a quantum superposition. This is anlogous to the simple qubit +example we considered in section \ref{sec:master}. + +\subsection{Quantum Measurement vs. Dissipation} + +This is where quantum measurement deviates from dissipation. If we +have access to a measurement record we can infer which Zeno subspace +is occupied, because we know that only one of them can be occupied at +any time. We have seen that since the density matrix cross-terms are +small we know \emph{a priori} that the individual wavefunctions +comprising the density matrix mixture will not be coherent +superpositions of different Zeno subspaces and thus we only have +classical uncertainty which means we can resort to clasical +probability methods. Each individual experiment will at any time be +predominantly in a single Zeno subspace with small cross-terms and +negligible occupations in the other subspaces. With no measurement +record our density matrix would be a mixture of all these +possibilities. We can try and determine the Zeno subspace around which +the state evolves in a single experiment from the number of +detections, $m$, in time $t$. + +The detection distribution on time-scales shorter than dissipation (so +we can approximate as if we were in a fully Zeno regime) can be +obtained by integrating over the detection times \cite{mekhov2009pra} +to get +\begin{equation} + P(m,t) = \sum_n \frac{[|c_n|^2 t]^m} {m!} e^{-|c_n|^2 t} \mathrm{Tr} (\rho_{nn}). +\end{equation} +For a state that is predominantly in one Zeno subspace, the +distribution will be approximately Poissonian (up to +$O(\nu^2 / \lambda^4)$, the population of the other +subspaces). Therefore, in a single experiment we will measure +$m = |c_0|^2t \pm \sqrt{|c_0|^2t}$ detections (note, we have assumed +$|c_0|^2 t$ is large enough to approximate the distribution as +normal. This is not necessary, we simply use it here to not have to +worry about the asymmetry in the deviation around the mean value). The +uncertainty does not come from the fact that $\lambda$ is not +infinite. The jumps are random events with a Poisson +distribution. Therefore, even in the full projective limit we will not +observe the same detection trajectory in each experiment even though +the system evolves in exactly the same way and remains in a perfectly +pure state. + +If the basis of $\c$ is continuous (e.g. free particle position or +momentum) then the deviation around the mean will be our upper bound +on the deviation of the system from a pure state evolving around a +single Zeno subspace. However, continuous systems are beyond the scope +of this work and we will confine ourselves to discrete systems. Though +it is important to remember that continuous systems can be treated +this way, but the error estimate (and thus the mixedness of the state) +will be different. + +For a discrete system it is easier to exclude all possibilities except +for one. The error in our estimate of $|c_0|^2$ in a single experiment +decreases as $1/\sqrt{t}$ and thus it can take a long time to +confidently determine $|c_0|^2$ to a sufficient precision this +way. However, since we know that it can only take one of the possible +values from the set $\{|c_n|^2 \}$ it is much easier to instead +exclude all the other values. + +In an experiment we can use Bayes' theorem to infer the state of our +system as follows +\begin{equation} + p(c_n = c_0 | m) = \frac{ p(m | c_n = c_0) p(c_n = c_0) }{ p(m) }, +\end{equation} +where $p(x)$ denotes the probability of the discrete event $x$ and +$p(x|y)$ the conditional probability of $x$ given $y$. We know that +$p(m | c_n = c_0)$ is simply given by a Poisson distribution with mean +$|c_0|^2 t$. $p(m)$ is just a normalising factor and $p(c_n = c_0)$ is +our \emph{a priori} knowledge of the state. Therefore, one can get the +probability of being in the right Zeno subspace from +\begin{align} + p(c_n & = c_0 | m) = \frac{ p_0(c_n = c_0) \frac{ \left( |c_0|^2 t + \right)^{2m} } {m!} e^{-|c_0|^2 t}} {\sum_n p_0(c_n) \frac{ + \left( |c_n|^2 t \right)^{2m} } {m!} e^{-|c_n|^2 t}} \nonumber \\ + & = p_0(c_n = c_0) \left[ \sum_n p_0(c_n) \left( \frac{ + |c_n|^2 } { |c_0|^2 } \right)^{2m} e^{\left( |c_0|^2 - + |c_n|^2 \right) t} \right]^{-1}, +\end{align} +where $p_0$ denotes probabilities at $t = 0$. In a real experiment one +could prepare the initial state to be close to the Zeno subspace of +interest and thus it would be easier to deduce the state. Furthermore, +in the middle of an experiment if we have already established the Zeno +subspace this will be reflected in these \emph{a priori} probabilities +again making it easier to infer the correct subspace. However, we will +consider the worst case scenario which might be useful if we don't +know the initial state or if the Zeno subspace changes during the +experiment, a uniform $p_0(c_n)$. + +This probability is a rather complicated function as $m$ is a +stochastic quantity that also increases with $t$. We want it to be as +close to $1$ as possible. In order to devise an appropriate condition +for this we note that in the first line all terms in the denominator +are Poisson distributions of $m$. Therefore, if the mean values +$|c_n|^2 t$ are sufficiently spaced out, only one of the terms in the +sum will be significant for a given $m$ and if this happens to be the +one that corresponds to $c_0$ we get a probability close to +unity. Therefore, we set the condition such that it is highly unlikely +that our measured $m$ could be produced by two different distributions +\begin{align} + \sqrt{|c_0|^2 t} \ll ||c_0|^2 - |c_n|^2| t, \forall n \ne 0 \\ + \sqrt{|c_n|^2 t} \ll ||c_0|^2 - |c_n|^2| t, \forall n \ne 0 +\end{align} +The left-hand side is the standard deviation of $m$ if the system was +in subspace $P_0$ or $P_n$. The right-hand side is the difference in +the mean detections between the subspace $n$ and the one we are +interested in. The condition becomes more strict if the subspaces +become less distinguishable as it becomes harder to confidently +determine the correct state. Once again, using $\c = \lambda \hat{o}$ +where $\hat{o} \sim O(1)$ we get +\begin{equation} + t \gg \frac{1}{\lambda^2} \frac{|o_{0,n}|^2} {(|o_0|^2 - |o_n|^2|)^2}. +\end{equation} +Since detections happen on average at an average rate of order +$\lambda^2$ we only need to wait for a few detections to satisfy this +condition. Therefore, we see that even in the worst case scenario of +complete ignorance of the state of the system we can very easily +determine the correct subspace. Once it is established for the first +time, the \emph{a priori} information can be updated and it will +become even easier to monitor the system. + +However, it is important to note that physically once the quantum +jumps deviate too much from the mean value the system is more likely +to change the Zeno subspace (due to measurement backaction) and the +detection rate will visibly change. Therefore, if we observe a +consistent detection rate it is extremely unlikely that it can be +produced by two different Zeno subspaces so in fact it is even easier +to determine the correct state, but the above estimate serves as a +good lower bound on the necessary detection time. + +Having derived the necessary conditions to confidently determine which +Zeno subspace is being observed in the experiment we can make another +approximation thanks to measurement which would be impossible in a +purely dissipative open system. If we observe a number of detections +consistent with the subspace $P_m = P_0$ we can set +$\hat{\rho}_{mn} \approx 0$ for all cases when both $m \ne 0$ and +$n \ne 0$ leaving our density matrix in the form +\begin{equation} + \label{eq:approxrho} + \hat{\rho} = \hat{\rho}_{00} + \sum_{r\ne0} (\hat{\rho}_{0r} + + \hat{\rho}_{r0}). +\end{equation} +We can do this, because the other states are inconsistent with the +measurement record. We know from the previous section that the system +must lie predominantly in only one of the Zeno subspaces and when that +is the case, $\hat{\rho}_{0r} \approx O(\nu/\lambda^2)$ and for +$m \ne 0$ and $n \ne 0$ we have +$\hat{\rho}_{mn} \approx O(\nu^2/\lambda^4)$. Therefore, this amounts +to keeping first order terms in $\nu/\lambda^2$ in our approximation. + +This is a crucial step as all $\hat{\rho}_{mm}$ matrices are +decoherence free subspaces and thus they can all coexist in a mixed +state decreasing the purity of the system without +measurement. Physically, this means we exclude trajectories in which +the Zeno subspace has changed (measurement isn't fully projective). By +substituting Eq. \eqref{eq:intermediate} into Eq. \eqref{eq:master} we +see that this happens at a rate of $\nu^2 / \lambda^2$. However, since +the two measurement outcomes cannot coexist any transition between +them happens in discrete transitions (which we know about from the +change in the detection rate as each Zeno subspace will correspond to +a different rate) and not as continuous coherent evolution. Therefore, +we can postselect in a manner similar to Refs. \cite{otterbach2014, + lee2014prx, lee2014prl}, but our requirements are significantly more +relaxed - we do not require a specific single trajectory, only that it +remains within a Zeno subspace. Furthermore, upon reaching a steady +state, these transitions become impossible as the coherences +vanish. This approximation is analogous to optical Raman transitions +where the population of the excited state is neglected. Here, we can +make a similar approximation and neglect all but one Zeno subspace +thanks to the additional knowledge we gain from knowing the +measurement outcomes. + +\subsection{The Non-Hermitian Hamiltonian} + +Rewriting the master equation using $\c = c_0 + \delta \c$, where +$c_0$ is the eigenvalue corresponding to the eigenspace defined by the +projector $P_0$ which we used to obtain the density matrix in +Eq. \eqref{eq:approxrho}, we get +\begin{equation} + \label{eq:finalrho} + \dot{\hat{\rho}} = -i \left( \H_\mathrm{eff} \hat{\rho} - \hat{\rho} + \H_\mathrm{eff}^\dagger \right) + \delta \c \hat{\rho} \delta \cd, +\end{equation} +\begin{equation} + \label{eq:Ham} + \H_\mathrm{eff} = \H_0 + i \left( c_0^*\c - \frac{|c_0|^2}{2} - \frac{\cd\c}{2} \right). +\end{equation} +The first term in Eq. \eqref{eq:finalrho} describes coherent evolution +due to the non-Hermitian Hamiltonian $\H_\mathrm{eff}$ and the second +term is decoherence due to our ignorance of measurement outcomes. When +we substitute our approximation of the density matrix +$\hat{\rho} = \hat{\rho}_{00} + \sum_{r\ne0} (\hat{\rho}_{0r} + +\hat{\rho}_{r0})$ into Eq. \eqref{eq:finalrho}, the last term +vanishes, $\delta \c \hat{\rho} \delta \cd = 0$. This happens, because +$\delta \c P_0 \hat{\rho} = \hat{\rho} P_0 \delta \c^\dagger = 0$. The +projector annihilates all states except for those with eigenvalue +$c_0$ and so the operator $\delta \c = \c - c_0$ will always evaluate +to $c_0 - c_0 = 0$. Recall that we defined +$\hat{\rho}_{mn} = P_m \hat{\rho} P_n$ which means that every term in +our approximate density matrix contains the projector $P_0$. However, +it is important to note that this argument does not apply to other +second order terms in the master equation, because some terms only +have the projector $P_0$ applied from one side, +e.g.~$\hat{\rho}_{0m}$. The term $\delta \c \hat{\rho} \delta \cd$ +applies the fluctuation operator from both sides so it does not matter +in this case, but it becomes relevant for terms such as +$\delta \cd \delta \c \hat{\rho}$. It is important to note that this +term does not automatically vanish, but when the explicit form of our +approximate density matrix is inserted, it is in fact zero. Therefore, +we can omit this term using the information we gained from +measurement, but keep other second order terms, such as +$\delta \cd \delta \c \rho$ in the Hamiltonian which are the origin of +other second-order dynamics. This could not be the case in a +dissipative system. + +Ultimately we find that a system under continuous measurement for +which $\lambda^2 \gg \nu$ in the Zeno subspace $P_0$ is described by +the deterministic non-Hermitian Hamiltonian $\H_\mathrm{eff}$ in +Eq. \eqref{eq:Ham} and thus obeys the following Schr\"{o}dinger +equation +\begin{equation} + i \frac{\mathrm{d} | \Psi \rangle}{\mathrm{d}t} = \left[\H_0 + i \left( + c_0^*\c - \frac{|c_0|^2}{2} - \frac{\cd\c}{2} \right) \right] | + \Psi \rangle. +\end{equation} +Of the three terms in the parentheses the first two represent the +effects of quantum jumps due to detections (which one can think of as +`reference frame' shifts between different degenerate eigenspaces) and +the last term is the non-Hermitian decay due to information gain from +no detections. It is important to emphasize that even though we +obtained a deterministic equation, we have not neglected the +stochastic nature of the detection events. The detection trajectory +seen in an experiment will have fluctuations around the mean +determined by the Zeno subspace, but there simply are many possible +measurement records with the same outcome. This is just like the fully +projective Zeno limit where the system remains perfectly pure in one +of the possible projections, but the detections remain randomly +distributed in time. + +One might then be concerned that purity is preserved even though we +might be averaging over many trajectories within this Zeno +subspace. We have neglected the small terms $\hat{\rho}_{m,n}$ +($m,n \ne 0$) which are $O(\nu^2/\lambda^4)$ and thus they are not +correctly accounted for by our approximation. This means that we have +an $O(\nu^2/\lambda^4)$ error in our density matrix. The purity +given by +\begin{equation} + \mathrm{Tr}(\hat{\rho}^2) = \mathrm{Tr}(\hat{\rho}^2_{00} + \sum_{m \ne + 0} \hat{\rho}_{0m}\hat{\rho}_{m0}) + \mathrm{Tr}(\sum_{m,n\ne0} + \hat{\rho}_{mn} \hat{\rho}_{nm}) +\end{equation} +where the second term contains the terms not accounted for by our +approximation thus introduces an $O(\nu^4/\lambda^8)$ +error. Therefore, this discrepancy is negligible in our +approximation. The pure state predicted by $\H_\mathrm{eff}$ is only +an approximation, albeit a good one, and the real state will be mixed +to a small extent. Whilst perfect purity within the Zeno subspace +$\hat{\rho}_{00}$ is expected due to the measurement's strong +decoupling effect, the nearly perfect purity when transitions outside +the Zeno subspace are included is a nontrivial result. Similarly, in +Raman transitions the population of the neglected excited state is +also non-zero, but negligible. Furthermore, this equation does not +actually require the adiabatic elimination used in +Eq. \eqref{eq:intermediate} (we only used it to convince ourselves +that the coherences are small) and such situations may be considered +provided all approximations remain valid. In a similar way the limit +of linear optics is derived from the physics of a two-level nonlinear +medium, when the population of the upper state is neglected and the +adiabatic elimination of coherences is not required. + +\subsection{Non-Hermitian Dynamics in Ultracold Gases} + +We finally return to our quantum gas model inside of a cavity. We +start by considering the simplest case of a global multi-site +measurement of the form $\hat{D} = \hat{N}_K = \sum_i^K \n_i$, where +the sum is over $K$ illuminated sites. The effective Hamiltonian +becomes +\begin{equation} + \label{eq:nHH2} + \hat{H}_\mathrm{eff} = \hat{H}_0 - i \gamma \left( \delta \hat{N}_K \right)^2, +\end{equation} +where $ \delta \hat{N}_K = \hat{N}_K - N^0_K$ and $N^0_K$ is the Zeno +subspace eigenvalue. It is now obvious that continuous measurement +squeezes the fluctuations in the measured quantity, as expected, and +that the only competing process is the system's own dynamics. + +In this case, if we adiabatically eliminate the density matrix +cross-terms and substitute Eq. \eqref{eq:intermediate} into +Eq. \eqref{eq:master} for this system we obtain an effective +Hamiltonian within the Zeno subspace defined by $N_K$ +\begin{equation} + \H_\varphi = P_0 \left[ \H_0 - i \frac{J^2}{\gamma} + \sum_\varphi \sum_{\substack{\langle i \in \varphi, j \in \varphi^\prime + \rangle \\ \langle k \in \varphi^\prime, l \in \varphi + \rangle}} b^\dagger_i b_j b^\dagger_k b_l \right] P_0, +\end{equation} +where $\varphi$ denotes a set of sites belonging to a single mode and +$\varphi^\prime$ is the set's complement (e.g. odd and even or +illuminated and non-illuminated sites) and $P_0$ is the projector onto +the eigenspace with $N_K^0$ atoms in the illuminated area. We focus on +the case when the second term is not only significant, but also leads +to dynamics within a Zeno subspace that are not allowed by +conventional quantum Zeno dynamics accounted for by the first +term. The second term represents second-order transitions via other +subspaces which act as intermediate states much like virtual states in +optical Raman transitions. This is in contrast to the conventional +understanding of the Zeno dynamics for infinitely frequent projective +measurements (corresponding to $\gamma \rightarrow \infty$) where such +processes are forbidden \cite{facchi2008}. Thus, it is the weak +quantum measurement that effectively couples the states. Note that +this is a special case of the equation in Eq. \eqref{eq:hz} which can +be obtained by considering a more general two mode setup. + +\subsection{Small System Example} + +To get clear physical insight, we initially consider three atoms in +three sites and choose our measurement operator such that +$\hat{D} = \n_2$, i.e.~only the middle site is subject to measurement, +and the Zeno subspace defined by $n_2 = 1$. Such an illumination +pattern can be achieved with global addressing by crossing two beams +and placing the nodes at the odd sites and the antinodes at even +sites. This means that $P_0 \H_0 P_0 = 0$. However, the +first and third sites are connected via the second term. Diagonalising +the Hamiltonian reveals that out of its ten eigenvalues all but three +have a significant negative imaginary component of the order $\gamma$ +which means that the corresponding eigenstates decay on a time scale +of a single quantum jump and thus quickly become negligible. The three +remaining eigenvectors are dominated by the linear superpositions of +the three Fock states $|2,1,0 \rangle$, $|1, 1, 1 \rangle$, and +$|0,1,2 \rangle$. Whilst it is not surprising that these components +are the only ones that remain as they are the only ones that actually +lie in the Zeno subspace $n_2 = 1$, it is impossible to solve the full +dynamics by just considering these Fock states alone as they are not +coupled to each other in $\hat{H}_0$. The components lying outside of +the Zeno subspace have to be included to allow intermediate steps to +occur via states that do not belong in this subspace, much like +virtual states in optical Raman transitions. + +An approximate solution for $U=0$ can be written for the +$\{|2,1,0 \rangle, |1,1,1 \rangle, |0,1,2 \rangle\}$ subspace by +multiplying each eigenvector with its corresponding time evolution +\begin{equation} + | \Psi(t) \rangle \propto \left( \begin{array}{c} + z_1 + \sqrt{2} z_2 e^{-6 J^2 t / \gamma} + z_3 e^{-12 J^2 t / \gamma} \\ + -\sqrt{2} \left(z_1 - z_3 e^{-12 J^2 t / \gamma} \right) \\ + z_1 - \sqrt{2} z_2 e^{-6 J^2 t / \gamma} + z_3 e^{-12 J^2 t / + \gamma} \\ + \end{array} + \right), \nonumber +\end{equation} +where $z_i$ denote the overlap between the eigenvectors and the +initial state, $z_i = \langle v_i | \Psi (0) \rangle$, with +$| v_1 \rangle = (1, -\sqrt{2}, 1)/2$, +$| v_2 \rangle = (1, 0, -1)/\sqrt{2}$, and +$| v_3 \rangle = (1, \sqrt{2}, 1)/2$. The steady state as +$t \rightarrow \infty$ is given by +$| v_1 \rangle = (1, -\sqrt{2}, 1)/2$. This solution is illustrated in +Fig. \ref{fig:comp} which clearly demonstrates dynamics beyond the +canonical understanding of quantum Zeno dynamics as tunnelling occurs +between states coupled via a different Zeno subspace. + +\begin{figure}[hbtp!] + \includegraphics[width=\linewidth]{comp} + \caption[Fock State Populations in a Zeno + Subspace]{Populations of the Fock states in the Zeno subspace + for $\gamma/J = 100$ and initial state $| 2,1,0 \rangle$. It + is clear that quantum Zeno dynamics occurs via Raman-like + processes even though none of these states are connected in + $\hat{H}_0$. The dynamics occurs via virtual intermediate + states outside the Zeno subspace. The system also tends to a + steady state which minimises tunnelling effectively + suppressing fluctuations. The lines are solutions to the + non-Hermitian Hamiltonian, and the dots are points from a + stochastic trajectory calculation.\label{fig:comp}} +\end{figure} + +\subsection{Steady State of non-Hermitian Dynamics} + +A distinctive difference between Bose-Hubbard model ground states and +the final steady state, +$| \Psi \rangle = [|2,1,0 \rangle - \sqrt{2} |1,1,1\rangle + +|0,1,2\rangle]/2$, is that its components are not in phase. Squeezing +due to measurement naturally competes with inter-site tunnelling which +tends to spread the atoms. However, from Eq. \eqref{eq:nHH2} we see +the final state will always be the eigenvector with the smallest +fluctuations as it will have an eigenvalue with the largest imaginary +component. This naturally corresponds to the state where tunnelling +between Zeno subspaces (here between every site) is minimised by +destructive matter-wave interference, i.e.~the tunnelling dark state +defined by $\hat{T} |\Psi \rangle = 0$, where +$\hat{T} = \sum_{\langle i, j \rangle} \bd_i b_j$. This is simply the +physical interpretation of the steady states we predicted for +Eq. \eqref{eq:master}. Crucially, this state can only be reached if +the dynamics aren't fully suppressed by measurement and thus, +counter-intuitively, the atomic dynamics cooperate with measurement to +suppress itself by destructive interference. Therefore, this effect is +beyond the scope of traditional quantum Zeno dynamics and presents a +new perspective on the competition between a system's short-range +dynamics and global measurement backaction. + +We now consider a one-dimensional lattice with $M$ sites so we extend +the measurement to $\hat{D} = \N_\text{even}$ where every even site is +illuminated. The wavefunction in a Zeno subspace must be an +eigenstate of $\c$ and we combine this with the requirement for it to +be in the dark state of the tunnelling operator (eigenstate of $\H_0$ +for $U = 0$) to derive the steady state. These two conditions in +momentum space are +\begin{equation} + \hat{T} | \Psi \rangle = \sum_{\text{RBZ}} \left[ \bd_k b_k - + \bd_{q} b_{q} \right] \cos(ka) |\Psi \rangle = 0, \nonumber +\end{equation} +\begin{equation} + \Delta \N |\Psi \rangle = \sum_{\text{RBZ}} \left[ \bd_k b_{-q} + + \bd_{-q} b_k \right] | \Psi \rangle= \Delta N |\Psi \rangle, \nonumber +\end{equation} +where $b_k = \frac{1}{\sqrt{M}} \sum_j e^{i k j a} b_j$, +$\Delta \hat{N} = \hat{D} - N/2$, $q = \pi/a - k$, $a$ is the lattice +spacing, $N$ the total atom number, and we perform summations over the +reduced Brillouin zone (RBZ), $-\pi/2a < k \le \pi/2a$, as the +symmetries of the system are clearer this way. Now we define +\begin{equation} +\hat{\alpha}_k^\dagger = \bd_k \bd_q - \bd_{-k} \bd_{-q}, +\end{equation} +\begin{equation} +\hat{\beta}_\varphi^\dagger = \bd_{\pi/2a} + \varphi \bd_{-\pi/2a}, +\end{equation} +where $\varphi = \Delta N / | \Delta N |$, which create the smallest +possible states that satisfy the two equations for $\Delta N = 0$ and +$\Delta N \ne 0$ respectively. Therefore, by noting that +\begin{align} + \left[ \hat{T}, \hat{\alpha}_k^\dagger \right] & = 0, \\ + \left[ \hat{T}, \hat{\beta}_\varphi^\dagger \right] & = 0, \\ + \left[ \Delta \N, \hat{\alpha}_k^\dagger \right] & = 0, \\ + \left[ \Delta \N, \hat{\beta}_\varphi^\dagger \right] & = \varphi + \hat{\beta}_\varphi^\dagger, +\end{align} +we can now write the equation for the $N$-particle steady state +\begin{equation} + \label{eq:ss} + | \Psi \rangle \propto \left[ \prod_{i=1}^{(N - |\Delta N|)/2} + \left( \sum_{k = 0}^{\pi/2a} \phi_{i,k} \hat{\alpha}_k^\dagger + \right) \right] \left( \hat{\beta}_\varphi^\dagger \right)^{| + \Delta N |} | 0 \rangle, \nonumber +\end{equation} +where $\phi_{i,k}$ are coefficients that depend on the trajectory +taken to reach this state and $|0 \rangle$ is the vacuum state defined +by $b_k |0 \rangle = 0$. Since this a dark state (an eigenstate of +$\H_0$) of the atomic dynamics, this state will remain stationary even +with measurement switched-off. Interestingly, this state is very +different from the ground states of the Bose-Hubbard Hamiltonian, it +is even orthogonal to the superfluid state, and thus it cannot be +obtained by cooling or projecting from an initial ground state. The +combination of tunnelling with measurement is necessary. + +\begin{figure}[hbtp!] + \includegraphics[width=\linewidth]{figure3} + \caption[Non-Hermitian Steady State]{A trajectory simulation + for eight atoms in eight sites, initially in + $|1,1,1,1,1,1,1,1 \rangle$, with periodic boundary + conditions and $\gamma/J = 100$. (a), The fluctuations in + $\c$ where the stochastic nature of the process is clearly + visible on a single trajectory level. However, the general + trend is captured by the non-Hermitian Hamiltonian. (b), The + local density variance. Whilst the fluctuations in the + global measurement operator decrease, the fluctuations in + local density increase due to tunnelling via states outside + the Zeno subspace. (c), The momentum distribution. The + initial Fock state has a flat distribution which with time + approaches the steady state distribution of two identical + and symmetric distributions centred at $k = \pi/2a$ and + $k = -\pi/2a$.\label{fig:steady}} +\end{figure} + +In order to prepare the steady state one has to run the experiment and +wait until the photocount rate remains constant for a sufficiently +long time. Such a trajectory is illustrated in Fig. \ref{fig:steady} +and compared to a deterministic trajectory calculated using the +non-Hermitian Hamiltonian. It is easy to see from +Fig. \ref{fig:steady}(a) how the stochastic fluctuations around the +mean value of the observable have no effect on the general behaviour +of the system in the strong measurement regime. By discarding these +fluctuations we no longer describe a pure state, but we showed how +this only leads to a negligible error. Fig. \ref{fig:steady}(b) shows +the local density variance in the lattice. Not only does it grow +showing evidence of tunnelling between illuminated and non-illuminated +sites, but it grows to significant values. This is in contrast to +conventional quantum Zeno dynamics where no tunnelling would be +allowed at all. Finally, Fig. \ref{fig:steady}(c) shows the momentum +distribution of the trajectory. We can clearly see that it deviates +significantly from the initial flat distribution of the Fock +state. Furthermore, the steady state does not have any atoms in the +$k=0$ state and thus is orthogonal to the superfluid state as +discussed. + +To obtain a state with a specific value of $\Delta N$ postselection +may be necessary, but otherwise it is not needed. The process can be +optimised by feedback control since the state is monitored at all +times \cite{ivanov2014}. Furthermore, the form of the measurement +operator is very flexible and it can easily be engineered by the +geometry of the optical setup \cite{elliott2015, mazzucchi2016} which +can be used to design a state with desired properties. + +\section{Conclusions} + +In this chapter we have demonstrated that global quantum measurement +backaction can efficiently compete with standard local processes in +many-body systems. This introduces a completely new energy and time +scale into quantum many-body research. This is made possbile by the +ability to structure the spatial profile of the measurement on a +microscopic scale comparable to the lattice period without the need +for single site addressing. The extreme flexibility of the setup +considered allowed us to effectively tailor long-range entanglement +and correlations present in the system. We showed that the competition +between the global backaction and usual atomic dynamics leads to the +production of spatially multimode macroscopic superpositions which +exhibit large-scale oscillatory dynamics which could be used for +quantum information and metrology. We subsequently demonstrated that +when on-site atomic interactions are introduced the dynamics become +much more complicated with different regimes of behaviour where +measurement and interactions can either compete or cooperate. In the +strong measurement regime we showed that conventional quantum Zeno +dynamics can be realised, but more interestingly, by considering a +strong, but not projective, limit of measurement we observe a new type +of nonlocal dynamics. It turns out that a global measurement scheme +leads to correlations between spatially separated tunnelling events +which conserve the Zeno subspace via Raman-like processes which would +be forbidden in the canonical fully projective limit. We subsequently +presented a rigorous analysis of the underlying process of this new +type of quantum Zeno dynamics in which we showed that in this limit +quantum trajectories can be described by a deterministic non-Hermitian +Hamiltonian. In contrast to previous works, it is independent of the +underlying system and there is no need to postselect a particular +exotic trajectory \cite{lee2014prx, lee2014prl}. Finally, we have +shown that the system will always tend towards the eigenstate of the +Hamiltonian with the best squeezing of the observable and the atomic +dynamics, which normally tend to spread the distribution, cooperates +with measurement to produce a state in which tunnelling is suppressed +by destructive matter-wave interference. A dark state of the +tunnelling operator will have zero fluctuations and we provided an +expression for the steady state which is significantly different from +the ground state of the Hamiltonian. This is in contrast to previous +works on dissipative state preparation where the steady state had to +be a dark state of the measurement operator \cite{diehl2008}. + +Such globally paired tunnelling due to a fundamentally new phenomenon, +global quantum measurement backaction, can enrich the physics of +long-range correlated systems beyond relatively short-range +interactions expected from standard dipole-dipole interactions +\cite{sowinski2012, omjyoti2015}. These nonlocal high-order processes +entangle regions of the optical lattice that are disconnected by the +measurement. Using different detection schemes, we showed how to +tailor density-density correlations between distant lattice +sites. Quantum optical engineering of nonlocal coupling to +environment, combined with quantum measurement, can allow the design +of nontrivial system-bath interactions, enabling new links to quantum +simulations~\cite{stannigel2013} and thermodynamics~\cite{erez2008} +and extend these directions to the field of non-Hermitian quantum +mechanics, where quantum optical setups are particularly +promising~\cite{lee2014prl}. Importantly, both systems and baths, +designed by our method, can be strongly correlated systems with +internal long-range entanglement. diff --git a/Preamble/preamble.tex b/Preamble/preamble.tex index 5106e1d..06df970 100644 --- a/Preamble/preamble.tex +++ b/Preamble/preamble.tex @@ -211,4 +211,7 @@ \renewcommand{\a}{a} % in case we decide to put hats on \renewcommand{\c}{\hat{c}} \newcommand{\cd}{\hat{c}^\dagger} -\renewcommand{\b}[1]{\mathbf{#1}} \ No newline at end of file +\renewcommand{\b}[1]{\mathbf{#1}} +\newcommand{\op}{\hat{o}} +\newcommand{\h}{\hat{h}} +\newcommand{\N}{\hat{N}} \ No newline at end of file diff --git a/References/references.bib b/References/references.bib index 792c0b3..061bb2e 100644 --- a/References/references.bib +++ b/References/references.bib @@ -978,3 +978,65 @@ doi = {10.1103/PhysRevA.87.043613}, year = {2005}, publisher = {American Physical Society}, } +@article{otterbach2014, + title = {Dissipative Preparation of Spatial Order in Rydberg-Dressed + Bose-Einstein Condensates}, + author = {Otterbach, Johannes and Lemeshko, Mikhail}, + journal = {Phys. Rev. Lett.}, + volume = {113}, + issue = {7}, + pages = {070401}, + numpages = {6}, + year = {2014}, + month = {Aug}, + publisher = {American Physical Society}, +} +@article{lee2014prx, + title = {Heralded Magnetism in Non-Hermitian Atomic Systems}, + author = {Lee, Tony E. and Chan, Ching-Kit}, + journal = {Phys. Rev. X}, + volume = {4}, + issue = {4}, + pages = {041001}, + numpages = {13}, + year = {2014}, + month = {Oct}, + publisher = {American Physical Society}, +} +@article{lee2014prl, + title = {Entanglement and Spin Squeezing in Non-Hermitian Phase Transitions}, + author = {Lee, Tony E. and Reiter, Florentin and Moiseyev, Nimrod}, + journal = {Phys. Rev. Lett.}, + volume = {113}, + issue = {25}, + pages = {250401}, + numpages = {5}, + year = {2014}, + month = {Dec}, + publisher = {American Physical Society}, +} +@article{stannigel2013, + title = {{Constrained Dynamics via the Zeno Effect in Quantum + Simulation: Implementing Non-Abelian Lattice Gauge + Theories with Cold Atoms}}, + author = {Stannigel, K. and Hauke, P. and Marcos, D. and Hafezi, + M. and Diehl, S. and Dalmonte, M. and Zoller, P.}, + journal = {Phys. Rev. Lett.}, + volume = {112}, + issue = {12}, + pages = {120406}, + numpages = {6}, + year = {2014}, + month = {Mar}, + publisher = {American Physical Society}, + doi = {10.1103/PhysRevLett.112.120406} +} +@article{erez2008, + title={{Repulsively bound atom pairs in an optical lattice}}, + author={Erez, N. and Gordon, G. and Nest, M. and Kurizki, G.}, + journal={Nature}, + volume={452}, + pages = {724}, + year={2008}, + publisher={Nature Publishing Group} +} diff --git a/thesis.tex b/thesis.tex index 23ab763..1770e82 100644 --- a/thesis.tex +++ b/thesis.tex @@ -1,7 +1,7 @@ % ******************************* PhD Thesis Template ************************** % Please have a look at the README.md file for info on how to use the template -\documentclass[a4paper,12pt,times,numbered,print,chapter]{Classes/PhDThesisPSnPDF} +\documentclass[a4paper,12pt,times,numbered,print]{Classes/PhDThesisPSnPDF} % ****************************************************************************** % ******************************* Class Options ********************************