arxiv: v1 [quant-ph] 15 Nov 2018

Transcription

1 Bayesian estimation of switching ates fo blinking emittes axiv:8.6627v [quant-ph] 5 Nov 28 Jemy Geody,, 2 Lachlan J Roges,, 2, Cameon M Roges, 3 Thomas Volz,, 2 and Alexei Gilchist, 2 Depatment of Physics and Astonomy, Macquaie Univesity, Sydney, NSW 222, Austalia. 2 ARC Cente of Excellence fo Engineeed Quantum Systems. 3 Launceston Chuch Gamma School, Mowbay Heights, Tasmania, 7248, Austalia. (Dated: Novembe 9, 28 Single quantum emittes of light ae valuable esouces fo engineeed quantum systems. They can function as obust single-photon geneatos, allow optical contol of single spins, povide eadout capabilities fo atomic-scale sensos, and povide intefaces between stationay and flying qubits. Envionmental factos can lead to single emittes exhibiting blinking, wheeby the fluoescence level switches between on and off states. Detailed chaacteisation of this blinking behaviou including detemining the switching ates, is often a poweful way to gain undestanding about the undelying physical mechanisms. Bight low-noise emittes lead to high-contast blinking, and simple thesholds can be used to extact the on and off intevals fom a time-seies of photon counts (and hence detemine the switching ates. Howeve, such appoaches become difficult fo emittes fluoescing at low levels, high noise, o fast switching ates. We develop a Bayesian appoach capable of infeing switching ates diectly fom the time-seies. This is able to deal with high levels of noise and fast switching in fluoescence taces. Moeove, the Bayesian infeence also yields a obust pictue of the paamete uncetainties, poviding a benefit also fo bight emittes. Finally, ou method is applicable to a boad ange of systems that do show behaviou analogous to a single blinking emitte. I. INTRODUCTION A wide vaiety of natual and man-made systems exhibit stochastic switching behaviou. They switch discetely between distinct states, and although the states ae not diectly obsevable thei occuence is coelated with a vaiable that may be measued. The incomplete knowledge about the system and its inteaction with the envionment often limits ou undestanding of the undelying causes of the blinking. A canonical example of a typical system showing blinking is povided by quantum light emittes, which ae valuable esouces fo engineeed quantum systems. Most quantum light emittes exhibit intemittency in thei fluoescence unde continuous excitation [], including diamond colou centes [2 6], quantum dots [7 ], nanowies [8, 2], nanoods [2] and oganic semiconductos [3]. In these quantum emittes, the fluoescence blinking is typically caused by physical mechanisms that fluctuate andomly (as opposed to peiodically and without egad fo thei histoy. The applications of single quantum emittes ae limited by the pesence of this blinking in thei fluoescence. It is impotant to undestand obseved blinking statistics in ode to analyse the cause of the blinking. The key insight into the physical mechanism often can be distilled fom the switching ates, and yet these ates ae often the hadest to extact fom the aw time seies of counts in paticula fo small light signals o noisy data. Blinking typically leads to step-like switches in the fluoescence time tace, as illustated in Fig.. Most common method that used to analyse the on and off lachlan.j.oges@quantum.diamonds states fom a binking time seies is theshold analysis [2, 4, 5]. Given a theshold, the emitte is taken to be in the on state if the intensity I n in an n th time bin, satisfies I n > I theshold and it is off if I n I theshold. Once the on and off states ae diffeentiated then the duations of states can be found fom the binned data seies [5, 6] and the appopiate paamete distibutions can be fitted to this data. Howeve, the choice of theshold intensity and binning time can significantly influence the statistics of the on and off states [6, 7]. These ae often essentially abitay choices. What is moe, this method is useful only if the on and off states ae clealy diffeentiated fom the blinking time seies as fo the simulated data in Fig.. The use of the theshold technique to analyse the on and off states becomes difficult fo emittes fluoescing at low intensity levels o when thee is a high backgound noise as depicted in Fig. 2 (a. In both these cases bief on and off intevals ae difficult to distinguish fom the noise spikes. A simila poblem aises if the blinking switches at ates high elative to the detecto inteval as shown in an exteme case in Fig. 2 (b. Hee we pesents a method utilising Bayesian statistical infeence to find the blinking ates diectly fom the photon-count data seies. The method woks fo blinking quantum light emittes independent of thei bightness o speed of switching. This method does not equie a theshold intensity level to diffeentiate between the on and off states. We model the blinking emittes on hidden, discete and continuous-time Makov pocesses. The accumulated photon counts ove a detecto inteval obscues o hides the undelying pocess which is assumed to be memoyless. The Bayesian infeence on this model then allows us to detemine the undelying switching ates of these hidden states. The infeence is able to tackle

2 <latexit sha_base64="x9yxdfeuw7gurlvjappoiaslw=">aaac3hicbvhljtmwfhxdawivgviymurikokqoivgoknygazlrqqnxjukvv2q7ojlxjhiw4infav/wd/gzioaedklo7opt5h995kc2zdmv4drtdu3p95+bufo/+g4epdo8ennvgeinrhflphvyypmke8ccpntkiik/nq/a7n3+ixjilt9glauifi+acns2a6xxm6/lhko7tvva+ylygqds6mr+n/hsio2gheos6yvluybemy4dthcdfyqogsyulbdtzkseuhac3vbfetxqzwkuikgt234+j8epsylzcktdvfswametrtrbwxnaxd7hxldb9a4xzuyzvi3jkpygbroohykd8vcntouol4jaihhytrmvmcauldsujdmfcgkzbojlgznzfhzkktmku/bdwzq6xncrzixpnx7podfkv2s/hpxupaj+vxd4tn/j/luujo+lbofgngegu/jcpfs97774cwfz+k4+/aqox67vfbeoqeo5coq6/rmxqpttaeectrd/ql/y4+rl+jb9h3s2k2v55ggyv/fwhmvjqdq==</latexit> <latexit sha_base64="x9yxdfeuw7gurlvjappoiaslw=">aaac3hicbvhljtmwfhxdawivgviymurikokqoivgoknygazlrqqnxjukvv2q7ojlxjhiw4infav/wd/gzioaedklo7opt5h995kc2zdmv4drtdu3p95+bufo/+g4epdo8ennvgeinrhflphvyypmke8ccpntkiik/nq/a7n3+ixjilt9glauifi+acns2a6xxm6/lhko7tvva+ylygqds6mr+n/hsio2gheos6yvluybemy4dthcdfyqogsyulbdtzkseuhac3vbfetxqzwkuikgt234+j8epsylzcktdvfswametrtrbwxnaxd7hxldb9a4xzuyzvi3jkpygbroohykd8vcntouol4jaihhytrmvmcauldsujdmfcgkzbojlgznzfhzkktmku/bdwzq6xncrzixpnx7podfkv2s/hpxupaj+vxd4tn/j/luujo+lbofgngegu/jcpfs97774cwfz+k4+/aqox67vfbeoqeo5coq6/rmxqpttaeectrd/ql/y4+rl+jb9h3s2k2v55ggyv/fwhmvjqdq==</latexit> <latexit sha_base64="x9yxdfeuw7gurlvjappoiaslw=">aaac3hicbvhljtmwfhxdawivgviymurikokqoivgoknygazlrqqnxjukvv2q7ojlxjhiw4infav/wd/gzioaedklo7opt5h995kc2zdmv4drtdu3p95+bufo/+g4epdo8ennvgeinrhflphvyypmke8ccpntkiik/nq/a7n3+ixjilt9glauifi+acns2a6xxm6/lhko7tvva+ylygqds6mr+n/hsio2gheos6yvluybemy4dthcdfyqogsyulbdtzkseuhac3vbfetxqzwkuikgt234+j8epsylzcktdvfswametrtrbwxnaxd7hxldb9a4xzuyzvi3jkpygbroohykd8vcntouol4jaihhytrmvmcauldsujdmfcgkzbojlgznzfhzkktmku/bdwzq6xncrzixpnx7podfkv2s/hpxupaj+vxd4tn/j/luujo+lbofgngegu/jcpfs97774cwfz+k4+/aqox67vfbeoqeo5coq6/rmxqpttaeectrd/ql/y4+rl+jb9h3s2k2v55ggyv/fwhmvjqdq==</latexit> <latexit sha_base64="x9yxdfeuw7gurlvjappoiaslw=">aaac3hicbvhljtmwfhxdawivgviymurikokqoivgoknygazlrqqnxjukvv2q7ojlxjhiw4infav/wd/gzioaedklo7opt5h995kc2zdmv4drtdu3p95+bufo/+g4epdo8ennvgeinrhflphvyypmke8ccpntkiik/nq/a7n3+ixjilt9glauifi+acns2a6xxm6/lhko7tvva+ylygqds6mr+n/hsio2gheos6yvluybemy4dthcdfyqogsyulbdtzkseuhac3vbfetxqzwkuikgt234+j8epsylzcktdvfswametrtrbwxnaxd7hxldb9a4xzuyzvi3jkpygbroohykd8vcntouol4jaihhytrmvmcauldsujdmfcgkzbojlgznzfhzkktmku/bdwzq6xncrzixpnx7podfkv2s/hpxupaj+vxd4tn/j/luujo+lbofgngegu/jcpfs97774cwfz+k4+/aqox67vfbeoqeo5coq6/rmxqpttaeectrd/ql/y4+rl+jb9h3s2k2v55ggyv/fwhmvjqdq==</latexit> <latexit sha_base64="m8g3ktsdoyg8zpycqydgz5mkc=">aaac3nicbvhljtmwfhxdawivgviymuribkjicalwacq2ladpojoiicqn47zw/yhsbqy3jdbfiwgu/hp/gbneyrjucydlrucfn6n5bojbl6z/r9g6zdu3tq7hd+5e+/+g/2dhydwn4aycdvcm7mslbncsynjtcz2jcqpwcn5eptz/9xizlwh27dckcqvf55ycc9qh8iihus9hvs32k3sc9ov3admabg3qahyw+pnxmjaskucfwdslae2kfozjvdafx3ljwqbqvwibcotj+wepwyxq+ai8h6mmaciszrdtp6pgzwfr4k4yugehtibthyty6dspo2yov6kbn39dtfzvjwokxgtng4gdxt26cmuno6sawbqebgnyuyoc4snc4v+ylkcqni8n+ckp2lyawgiwkhm/fkyy6wxnccln/e7dmu2yszfelg7jpzc4b5k/+f8y253xbd4twmj22hgfh4ut7fvugpnstnixef8yoxyzufeeoyeoueioffoeldr2ickflob/qffkcfo6/rt+j7htqabf48qookfv4dzcbqfw==</latexit> <latexit sha_base64="m8g3ktsdoyg8zpycqydgz5mkc=">aaac3nicbvhljtmwfhxdawivgviymuribkjicalwacq2ladpojoiicqn47zw/yhsbqy3jdbfiwgu/hp/gbneyrjucydlrucfn6n5bojbl6z/r9g6zdu3tq7hd+5e+/+g/2dhydwn4aycdvcm7mslbncsynjtcz2jcqpwcn5eptz/9xizlwh27dckcqvf55ycc9qh8iihus9hvs32k3sc9ov3admabg3qahyw+pnxmjaskucfwdslae2kfozjvdafx3ljwqbqvwibcotj+wepwyxq+ai8h6mmaciszrdtp6pgzwfr4k4yugehtibthyty6dspo2yov6kbn39dtfzvjwokxgtng4gdxt26cmuno6sawbqebgnyuyoc4snc4v+ylkcqni8n+ckp2lyawgiwkhm/fkyy6wxnccln/e7dmu2yszfelg7jpzc4b5k/+f8y253xbd4twmj22hgfh4ut7fvugpnstnixef8yoxyzufeeoyeoueioffoeldr2ickflob/qffkcfo6/rt+j7htqabf48qookfv4dzcbqfw==</latexit> <latexit sha_base64="m8g3ktsdoyg8zpycqydgz5mkc=">aaac3nicbvhljtmwfhxdawivgviymuribkjicalwacq2ladpojoiicqn47zw/yhsbqy3jdbfiwgu/hp/gbneyrjucydlrucfn6n5bojbl6z/r9g6zdu3tq7hd+5e+/+g/2dhydwn4aycdvcm7mslbncsynjtcz2jcqpwcn5eptz/9xizlwh27dckcqvf55ycc9qh8iihus9hvs32k3sc9ov3admabg3qahyw+pnxmjaskucfwdslae2kfozjvdafx3ljwqbqvwibcotj+wepwyxq+ai8h6mmaciszrdtp6pgzwfr4k4yugehtibthyty6dspo2yov6kbn39dtfzvjwokxgtng4gdxt26cmuno6sawbqebgnyuyoc4snc4v+ylkcqni8n+ckp2lyawgiwkhm/fkyy6wxnccln/e7dmu2yszfelg7jpzc4b5k/+f8y253xbd4twmj22hgfh4ut7fvugpnstnixef8yoxyzufeeoyeoueioffoeldr2ickflob/qffkcfo6/rt+j7htqabf48qookfv4dzcbqfw==</latexit> <latexit sha_base64="m8g3ktsdoyg8zpycqydgz5mkc=">aaac3nicbvhljtmwfhxdawivgviymuribkjicalwacq2ladpojoiicqn47zw/yhsbqy3jdbfiwgu/hp/gbneyrjucydlrucfn6n5bojbl6z/r9g6zdu3tq7hd+5e+/+g/2dhydwn4aycdvcm7mslbncsynjtcz2jcqpwcn5eptz/9xizlwh27dckcqvf55ycc9qh8iihus9hvs32k3sc9ov3admabg3qahyw+pnxmjaskucfwdslae2kfozjvdafx3ljwqbqvwibcotj+wepwyxq+ai8h6mmaciszrdtp6pgzwfr4k4yugehtibthyty6dspo2yov6kbn39dtfzvjwokxgtng4gdxt26cmuno6sawbqebgnyuyoc4snc4v+ylkcqni8n+ckp2lyawgiwkhm/fkyy6wxnccln/e7dmu2yszfelg7jpzc4b5k/+f8y253xbd4twmj22hgfh4ut7fvugpnstnixef8yoxyzufeeoyeoueioffoeldr2ickflob/qffkcfo6/rt+j7htqabf48qookfv4dzcbqfw==</latexit> <latexit sha_base64="3hzqsigydsirhg9lg+m+5i5kjii=">aaac3xicbvhlbtqwfpwert4tbbky4iq2dckoqwqpxyscxspy2arkmbxzoxo/idqiget2ialnva/ad/g5moujptkk6ovf4hn7ypw69l7yi6dfvo3z3de/h9bw8fpd7bf3jidwmomattdkwtlbfzs47gq7qwdwqp2wi7fd/3tz8xytwxw9wsklbqfm4pueb9iq/zkjmyvbo9jb2nfeftqnygqes6mu2p/usvpokylebk5jwuibem4fczhcd5yvgndwoknygqqlqfsblceqftmcd9trabzlzouh9phlyco8yy85xdpdpxawuszvb2nnaz5e39aanm78twq7qxjff4lmjcbo425buokgusdwaqapiyg6tkyoc7snm4vu6baslbvm5cg/jqubs5alqrlykufcpbm9yllhosm/gth6txtdjspvfqdknztcn8+f/p2ab6c74nu8wygr7blypw8xj5n23wuk2jumyfdxidt+tb7+lnqex6bui6a6rb/qezogiit6gx6h39es+hp9i75fsapr+s9tnkjo5z85uol</latexit> <latexit sha_base64="3hzqsigydsirhg9lg+m+5i5kjii=">aaac3xicbvhlbtqwfpwert4tbbky4iq2dckoqwqpxyscxspy2arkmbxzoxo/idqiget2ialnva/ad/g5moujptkk6ovf4hn7ypw69l7yi6dfvo3z3de/h9bw8fpd7bf3jidwmomattdkwtlbfzs47gq7qwdwqp2wi7fd/3tz8xytwxw9wsklbqfm4pueb9iq/zkjmyvbo9jb2nfeftqnygqes6mu2p/usvpokylebk5jwuibem4fczhcd5yvgndwoknygqqlqfsblceqftmcd9trabzlzouh9phlyco8yy85xdpdpxawuszvb2nnaz5e39aanm78twq7qxjff4lmjcbo425buokgusdwaqapiyg6tkyoc7snm4vu6baslbvm5cg/jqubs5alqrlykufcpbm9yllhosm/gth6txtdjspvfqdknztcn8+f/p2ab6c74nu8wygr7blypw8xj5n23wuk2jumyfdxidt+tb7+lnqex6bui6a6rb/qezogiit6gx6h39es+hp9i75fsapr+s9tnkjo5z85uol</latexit> <latexit sha_base64="3hzqsigydsirhg9lg+m+5i5kjii=">aaac3xicbvhlbtqwfpwert4tbbky4iq2dckoqwqpxyscxspy2arkmbxzoxo/idqiget2ialnva/ad/g5moujptkk6ovf4hn7ypw69l7yi6dfvo3z3de/h9bw8fpd7bf3jidwmomattdkwtlbfzs47gq7qwdwqp2wi7fd/3tz8xytwxw9wsklbqfm4pueb9iq/zkjmyvbo9jb2nfeftqnygqes6mu2p/usvpokylebk5jwuibem4fczhcd5yvgndwoknygqqlqfsblceqftmcd9trabzlzouh9phlyco8yy85xdpdpxawuszvb2nnaz5e39aanm78twq7qxjff4lmjcbo425buokgusdwaqapiyg6tkyoc7snm4vu6baslbvm5cg/jqubs5alqrlykufcpbm9yllhosm/gth6txtdjspvfqdknztcn8+f/p2ab6c74nu8wygr7blypw8xj5n23wuk2jumyfdxidt+tb7+lnqex6bui6a6rb/qezogiit6gx6h39es+hp9i75fsapr+s9tnkjo5z85uol</latexit> <latexit sha_base64="3hzqsigydsirhg9lg+m+5i5kjii=">aaac3xicbvhlbtqwfpwert4tbbky4iq2dckoqwqpxyscxspy2arkmbxzoxo/idqiget2ialnva/ad/g5moujptkk6ovf4hn7ypw69l7yi6dfvo3z3de/h9bw8fpd7bf3jidwmomattdkwtlbfzs47gq7qwdwqp2wi7fd/3tz8xytwxw9wsklbqfm4pueb9iq/zkjmyvbo9jb2nfeftqnygqes6mu2p/usvpokylebk5jwuibem4fczhcd5yvgndwoknygqqlqfsblceqftmcd9trabzlzouh9phlyco8yy85xdpdpxawuszvb2nnaz5e39aanm78twq7qxjff4lmjcbo425buokgusdwaqapiyg6tkyoc7snm4vu6baslbvm5cg/jqubs5alqrlykufcpbm9yllhosm/gth6txtdjspvfqdknztcn8+f/p2ab6c74nu8wygr7blypw8xj5n23wuk2jumyfdxidt+tb7+lnqex6bui6a6rb/qezogiit6gx6h39es+hp9i75fsapr+s9tnkjo5z85uol</latexit> <latexit sha_base64="viyeeg5o5at2rgmwz3gbjwwjqfi=">aaac23icbvfna9taef2x6n6lbthxtskqk9geoh2fak99jhchacsyuasbn4p8tuqsfv9trb6awhxt/v/rf9ov4kjkz2dh8ebte8xmwxnmbjl8hqv37t67/2dvyfjo8zonz/ypnp8zwike6q4uclgorm4sqyy/g8gii5dgtvx+6/vqzasouplxggsbs8kwjildqvck8mu/hytjpk9of6qbezfmn84prn7xstbeolevgzcxnaluoc2jhfy5o3bgugkltjwwh2cvg7acaxnxxrlhqozx5kegikth/prw88uulpf2tnuzgrmiy9ai9mo2z3ovk23qyxi/dfy2tdwjtomjr8miqqntvvdgnpkba86nf9aiuos3guysl6ksamtv5qugnulnucglxzj9monbsuveaqtnpda9yoq6/jbpe2ei+wts5ug+av/h7ptdktdm3clki9c6f/ulp9n3wvk2tpnx+ukwpj7a3h6pvcsvyvusknfkmhwkj2rckfmrh+qx+ruwdfgw/d9whqmnn9ekeefp/8bbcbpmq==</latexit> <latexit sha_base64="viyeeg5o5at2rgmwz3gbjwwjqfi=">aaac23icbvfna9taef2x6n6lbthxtskqk9geoh2fak99jhchacsyuasbn4p8tuqsfv9trb6awhxt/v/rf9ov4kjkz2dh8ebte8xmwxnmbjl8hqv37t67/2dvyfjo8zonz/ypnp8zwike6q4uclgorm4sqyy/g8gii5dgtvx+6/vqzasouplxggsbs8kwjildqvck8mu/hytjpk9of6qbezfmn84prn7xstbeolevgzcxnaluoc2jhfy5o3bgugkltjwwh2cvg7acaxnxxrlhqozx5kegikth/prw88uulpf2tnuzgrmiy9ai9mo2z3ovk23qyxi/dfy2tdwjtomjr8miqqntvvdgnpkba86nf9aiuos3guysl6ksamtv5qugnulnucglxzj9monbsuveaqtnpda9yoq6/jbpe2ei+wts5ug+av/h7ptdktdm3clki9c6f/ulp9n3wvk2tpnx+ukwpj7a3h6pvcsvyvusknfkmhwkj2rckfmrh+qx+ruwdfgw/d9whqmnn9ekeefp/8bbcbpmq==</latexit> <latexit sha_base64="viyeeg5o5at2rgmwz3gbjwwjqfi=">aaac23icbvfna9taef2x6n6lbthxtskqk9geoh2fak99jhchacsyuasbn4p8tuqsfv9trb6awhxt/v/rf9ov4kjkz2dh8ebte8xmwxnmbjl8hqv37t67/2dvyfjo8zonz/ypnp8zwike6q4uclgorm4sqyy/g8gii5dgtvx+6/vqzasouplxggsbs8kwjildqvck8mu/hytjpk9of6qbezfmn84prn7xstbeolevgzcxnaluoc2jhfy5o3bgugkltjwwh2cvg7acaxnxxrlhqozx5kegikth/prw88uulpf2tnuzgrmiy9ai9mo2z3ovk23qyxi/dfy2tdwjtomjr8miqqntvvdgnpkba86nf9aiuos3guysl6ksamtv5qugnulnucglxzj9monbsuveaqtnpda9yoq6/jbpe2ei+wts5ug+av/h7ptdktdm3clki9c6f/ulp9n3wvk2tpnx+ukwpj7a3h6pvcsvyvusknfkmhwkj2rckfmrh+qx+ruwdfgw/d9whqmnn9ekeefp/8bbcbpmq==</latexit> <latexit sha_base64="viyeeg5o5at2rgmwz3gbjwwjqfi=">aaac23icbvfna9taef2x6n6lbthxtskqk9geoh2fak99jhchacsyuasbn4p8tuqsfv9trb6awhxt/v/rf9ov4kjkz2dh8ebte8xmwxnmbjl8hqv37t67/2dvyfjo8zonz/ypnp8zwike6q4uclgorm4sqyy/g8gii5dgtvx+6/vqzasouplxggsbs8kwjildqvck8mu/hytjpk9of6qbezfmn84prn7xstbeolevgzcxnaluoc2jhfy5o3bgugkltjwwh2cvg7acaxnxxrlhqozx5kegikth/prw88uulpf2tnuzgrmiy9ai9mo2z3ovk23qyxi/dfy2tdwjtomjr8miqqntvvdgnpkba86nf9aiuos3guysl6ksamtv5qugnulnucglxzj9monbsuveaqtnpda9yoq6/jbpe2ei+wts5ug+av/h7ptdktdm3clki9c6f/ulp9n3wvk2tpnx+ukwpj7a3h6pvcsvyvusknfkmhwkj2rckfmrh+qx+ruwdfgw/d9whqmnn9ekeefp/8bbcbpmq==</latexit> <latexit sha_base64="iwms4tcu3pbwp6oyxl6ifbqxp/g=">aaac3hicbvhlbtnafj2yvzgvfpzsbiwkvpfticekvwldskhngyk2fvkkowyl82mqcgdhbuqmy66gw/hp/gbxm6q6qrxguno3dpn6n5bac6ss9o/g+jo3xv3h+w9jb89fvl2f7b8xoakpoiciuzlgcszmtdosy43ssdqvrcxpat6/dnvfim5lfbaokweg2zwrcomzmwgdxs5juj+kw7qvgmwdespo+nb4e8xu6qwvdcwdpjlmpxnmaci5z6oc5qszwqfsxoo4ghpkl8jxxglhp+ki976kmaegszdpn5/gbwmzwxjnwpmmd23mcyeavehzettxkvepvubv6hbjjutaosxafna46dwu2y8iwzshxfbwdesdaajkswqfxyavxi+poiudomqiy4cdz2rqc5iltjdvzfcgkuk5wlmsf6tr+x6ht5ddnkvx2cvb2iflnw+lv/4/5dozvinabrjrhbvv43dxbpu+u+akh2bpmpvyljn8uln9hnqjxqo3kepvsh6ji7qcbhesx6hx5hx6ofxlcs2nbps/lcvoqt/u7vqgw==</latexit> <latexit sha_base64="iwms4tcu3pbwp6oyxl6ifbqxp/g=">aaac3hicbvhlbtnafj2yvzgvfpzsbiwkvpfticekvwldskhngyk2fvkkowyl82mqcgdhbuqmy66gw/hp/gbxm6q6qrxguno3dpn6n5bac6ss9o/g+jo3xv3h+w9jb89fvl2f7b8xoakpoiciuzlgcszmtdosy43ssdqvrcxpat6/dnvfim5lfbaokweg2zwrcomzmwgdxs5juj+kw7qvgmwdespo+nb4e8xu6qwvdcwdpjlmpxnmaci5z6oc5qszwqfsxoo4ghpkl8jxxglhp+ki976kmaegszdpn5/gbwmzwxjnwpmmd23mcyeavehzettxkvepvubv6hbjjutaosxafna46dwu2y8iwzshxfbwdesdaajkswqfxyavxi+poiudomqiy4cdz2rqc5iltjdvzfcgkuk5wlmsf6tr+x6ht5ddnkvx2cvb2iflnw+lv/4/5dozvinabrjrhbvv43dxbpu+u+akh2bpmpvyljn8uln9hnqjxqo3kepvsh6ji7qcbhesx6hx5hx6ofxlcs2nbps/lcvoqt/u7vqgw==</latexit> <latexit sha_base64="iwms4tcu3pbwp6oyxl6ifbqxp/g=">aaac3hicbvhlbtnafj2yvzgvfpzsbiwkvpfticekvwldskhngyk2fvkkowyl82mqcgdhbuqmy66gw/hp/gbxm6q6qrxguno3dpn6n5bac6ss9o/g+jo3xv3h+w9jb89fvl2f7b8xoakpoiciuzlgcszmtdosy43ssdqvrcxpat6/dnvfim5lfbaokweg2zwrcomzmwgdxs5juj+kw7qvgmwdespo+nb4e8xu6qwvdcwdpjlmpxnmaci5z6oc5qszwqfsxoo4ghpkl8jxxglhp+ki976kmaegszdpn5/gbwmzwxjnwpmmd23mcyeavehzettxkvepvubv6hbjjutaosxafna46dwu2y8iwzshxfbwdesdaajkswqfxyavxi+poiudomqiy4cdz2rqc5iltjdvzfcgkuk5wlmsf6tr+x6ht5ddnkvx2cvb2iflnw+lv/4/5dozvinabrjrhbvv43dxbpu+u+akh2bpmpvyljn8uln9hnqjxqo3kepvsh6ji7qcbhesx6hx5hx6ofxlcs2nbps/lcvoqt/u7vqgw==</latexit> <latexit sha_base64="iwms4tcu3pbwp6oyxl6ifbqxp/g=">aaac3hicbvhlbtnafj2yvzgvfpzsbiwkvpfticekvwldskhngyk2fvkkowyl82mqcgdhbuqmy66gw/hp/gbxm6q6qrxguno3dpn6n5bac6ss9o/g+jo3xv3h+w9jb89fvl2f7b8xoakpoiciuzlgcszmtdosy43ssdqvrcxpat6/dnvfim5lfbaokweg2zwrcomzmwgdxs5juj+kw7qvgmwdespo+nb4e8xu6qwvdcwdpjlmpxnmaci5z6oc5qszwqfsxoo4ghpkl8jxxglhp+ki976kmaegszdpn5/gbwmzwxjnwpmmd23mcyeavehzettxkvepvubv6hbjjutaosxafna46dwu2y8iwzshxfbwdesdaajkswqfxyavxi+poiudomqiy4cdz2rqc5iltjdvzfcgkuk5wlmsf6tr+x6ht5ddnkvx2cvb2iflnw+lv/4/5dozvinabrjrhbvv43dxbpu+u+akh2bpmpvyljn8uln9hnqjxqo3kepvsh6ji7qcbhesx6hx5hx6ofxlcs2nbps/lcvoqt/u7vqgw==</latexit> <latexit sha_base64="igpjzynwalso/xktct5wfpv5ks=">aaac23icbvfna9taef/ujuyq95qjwfhoykgipxdopccu4jhgctnaj53f+yf2vw2usqfesi899nl+l/6p/puufbciowmljzdv32nmyoozy5pk7ydyefdwepdvfdj2fpx+wfhf4yvwuki6q4pclgorm4sgyy/gygii5dgulx/a/vgzasouplecgsbc8nmjilfhp8xittpfjzjhfw2ddaisq6z6chgtz5ttbyolevgzcrnklsoc2jhfy5xbcugsfthuwh2cvk7hcbxfxxrluqpjx5kegikphvprw88m4vmsvsnbdsxpscqxqxe6zwtp9nster9vult5++lhsmqtijpbdc85pfvubuamyustxhgdvzi8wsvqqk3fajhlvkzkcjczji8uelandkhueaypzeulii7qq3cztuo2hdpndtckzu4xe2wxmnzf/9/5htplvtmxdgbbnuxyu9bdpn++7ds6yyzom/h8enj+va75ii8jm9jst6ru/krnjerowrjfpbf5hdqbf+db8h3w2kwwp95sxov/pwhkczppw==</latexit> <latexit sha_base64="igpjzynwalso/xktct5wfpv5ks=">aaac23icbvfna9taef/ujuyq95qjwfhoykgipxdopccu4jhgctnaj53f+yf2vw2usqfesi899nl+l/6p/puufbciowmljzdv32nmyoozy5pk7ydyefdwepdvfdj2fpx+wfhf4yvwuki6q4pclgorm4sgyy/gygii5dgulx/a/vgzasouplecgsbc8nmjilfhp8xittpfjzjhfw2ddaisq6z6chgtz5ttbyolevgzcrnklsoc2jhfy5xbcugsfthuwh2cvk7hcbxfxxrluqpjx5kegikphvprw88m4vmsvsnbdsxpscqxqxe6zwtp9nster9vult5++lhsmqtijpbdc85pfvubuamyustxhgdvzi8wsvqqk3fajhlvkzkcjczji8uelandkhueaypzeulii7qq3cztuo2hdpndtckzu4xe2wxmnzf/9/5htplvtmxdgbbnuxyu9bdpn++7ds6yyzom/h8enj+va75ii8jm9jst6ru/krnjerowrjfpbf5hdqbf+db8h3w2kwwp95sxov/pwhkczppw==</latexit> <latexit sha_base64="igpjzynwalso/xktct5wfpv5ks=">aaac23icbvfna9taef/ujuyq95qjwfhoykgipxdopccu4jhgctnaj53f+yf2vw2usqfesi899nl+l/6p/puufbciowmljzdv32nmyoozy5pk7ydyefdwepdvfdj2fpx+wfhf4yvwuki6q4pclgorm4sgyy/gygii5dgulx/a/vgzasouplecgsbc8nmjilfhp8xittpfjzjhfw2ddaisq6z6chgtz5ttbyolevgzcrnklsoc2jhfy5xbcugsfthuwh2cvk7hcbxfxxrluqpjx5kegikphvprw88m4vmsvsnbdsxpscqxqxe6zwtp9nster9vult5++lhsmqtijpbdc85pfvubuamyustxhgdvzi8wsvqqk3fajhlvkzkcjczji8uelandkhueaypzeulii7qq3cztuo2hdpndtckzu4xe2wxmnzf/9/5htplvtmxdgbbnuxyu9bdpn++7ds6yyzom/h8enj+va75ii8jm9jst6ru/krnjerowrjfpbf5hdqbf+db8h3w2kwwp95sxov/pwhkczppw==</latexit> <latexit sha_base64="igpjzynwalso/xktct5wfpv5ks=">aaac23icbvfna9taef/ujuyq95qjwfhoykgipxdopccu4jhgctnaj53f+yf2vw2usqfesi899nl+l/6p/puufbciowmljzdv32nmyoozy5pk7ydyefdwepdvfdj2fpx+wfhf4yvwuki6q4pclgorm4sgyy/gygii5dgulx/a/vgzasouplecgsbc8nmjilfhp8xittpfjzjhfw2ddaisq6z6chgtz5ttbyolevgzcrnklsoc2jhfy5xbcugsfthuwh2cvk7hcbxfxxrluqpjx5kegikphvprw88m4vmsvsnbdsxpscqxqxe6zwtp9nster9vult5++lhsmqtijpbdc85pfvubuamyustxhgdvzi8wsvqqk3fajhlvkzkcjczji8uelandkhueaypzeulii7qq3cztuo2hdpndtckzu4xe2wxmnzf/9/5htplvtmxdgbbnuxyu9bdpn++7ds6yyzom/h8enj+va75ii8jm9jst6ru/krnjerowrjfpbf5hdqbf+db8h3w2kwwp95sxov/pwhkczppw==</latexit> 2 2 Count Time step on off theshold 2 Feqency FIG.. Simulated blinking timetace with clea on and off states. The dotted line epesents a theshold, which is used to distinguish the on and off states. (a Count (b Count Time step Fequency FIG. 2. (a Simulated blinking time tace with highe backgound noise but same switching ates as in Fig.. (b A high ate of switching events elative to the detecto inteval makes it hade even to see the blinking ates anging anywhee fom low ates to vey high ates, whee multiple switch events occu within a given detecto inteval. We have modeled thee behavious fo the state elative to the detecto inteval. In section II we analyse simplest blinking pocess using discete time Makov pocess, whee the emitte is supposed to eithe stay on o off duing each detecto inteval and switching of state can occu only at the inteval boundaies. The exteme case of continuous switching within the detecto intevals is analysed in section III using a continuous time Makov pocess. In eal expeiments only a finite numbe of switching events ae likely to occu within any given detecto inteval athe than no switching o infinite switching events. Such a pocess is discussed in section IV using a discete time Makov pocess and consideing subintevals inside the detecto inteval. II. DISCRETE-TIME MARKOV CHAIN SINGLE-STEP MODEL A. Model desciption The citical step in making an infeence on the blinking ates fom the obseved fluoescence data is constucting a model of how the counts aise given the stuctue of the poblem. We conside that the emitte has only two elevant states that affect the fluoescence levels: an on state and an off state, and thee is some mechanism that causes the emitte to switch between these states. We model the switching as a Makov pocess, whee the cuent state is sufficient to detemine the futue dynamics. Fo simplicity, we begin by consideing the state to be fixed fo the entie detecto inteval and it can only switch at the boundaies with switch-on and switch-off pobabilities α and β espectively, whee the subscipt epesents the allowed numbe of state changes pe detecto inteval. This notation anticipates a discete-time Makov chain which allows fo multiple switching events pe detecto inteval, which is discussed in late sections. The Makov chain fo this moe geneal model is depicted on the left in Fig. 3, whee d = coesponds to the single switching oppotunity case. a d b d On Off a d b d FIG. 3. Makov pocess modeling the evolution of the state. Left: a discete time Makov pocess with d time steps pe detection inteval. In this case α d and β d ae the switching pobabilities. Right: a continuous time Makov pocess whee α and β ae switching ates. The blinking model also contains the fluoescence and backgound ates λ and µ. If the emitte is in an offstate then the detecto can still egiste counts (known as dak-counts due to noise pocesses in the detecto o electonics. These counts ae well modeled by a Poissonian pocess with a ate µ. While fluoescing, it is assumed that the emitte is being diven by a stong classical pump so the statistics of the photon counts ae also well modelled by a Poissonian pocess (i.e. shot noise with an oveall ate λ that includes detecto inefficiencies. The obseved counts ae then Poissonian with a ate of eithe µ o µ + λ since the combination of two independent Poissonian pocesses is again Poissonian. Given the state s t at the beginning of time-step t, λ and µ the pobability of seeing c t counts is { P(ct µ s P (c t s t λ µ I = t = ( P(c t µ + λ s t = whee s t = epesents the on state and s t = epesents the off state fo inteval t, I tags the single inteval model as the backgound infomation, and P(n λ = n λ exp( λ/n! is the pobability of n given a Poisson distibution with ate λ. Clealy, if we knew the state of the emitte at each data point it would be a simple matte to infe the switching ates. Unfotunately these states ae not diectly obsevable, and wost still, because the states ae not obseved a On Off b

3 <latexit sha_base64="mys89qmsy7tmfxhlnwnn4jsy8=">aaac4xicbvhlbtnafj2yr4t5tbbkm2ahssgylss6gddskhnwym2ouvjjblhu7mnsi48whsebswbobh+a/+hebpcbpluy6ovfmob3vuudtp7yc6dfvo3z3de/h9bw8fpd7bf3litgmzhzejjt2whepnb+hqmnpastbvzkfvov3xf/e7dogh2my5qxcuzazaqddftpjm2bqi8pvs79zc9jh2lfdbtkk5cqvrn9gd/iqlhjeiamqtnxllay9mcrcek93fcni7xwbyw52nmorq49dybxhewxfj6zp6hkagxv3z9kn6+jiwuzoznjyntgfxnea5trvuhaebpxktfxg3odspw6lpbtlvkyrfa3tnkanwnkgchkamcvcajsdgwwgya9xoflnzpqcpw2lyoifz2vbsnbzyzpsq8lfhx7wwwsfbbx+c2hxpnfnnym8vb2sfmw2l5/8/5pvpahbdauwqj223sfh4tnmfbfbst7mmh28uy+g5+3yjlwghg3pjd8oeckrfh5il8il/i74hfx6nvfcatry/xlki6+q822uv7</latexit> <latexit sha_base64="mys89qmsy7tmfxhlnwnn4jsy8=">aaac4xicbvhlbtnafj2yr4t5tbbkm2ahssgylss6gddskhnwym2ouvjjblhu7mnsi48whsebswbobh+a/+hebpcbpluy6ovfmob3vuudtp7yc6dfvo3z3de/h9bw8fpd7bf3litgmzhzejjt2whepnb+hqmnpastbvzkfvov3xf/e7dogh2my5qxcuzazaqddftpjm2bqi8pvs79zc9jh2lfdbtkk5cqvrn9gd/iqlhjeiamqtnxllay9mcrcek93fcni7xwbyw52nmorq49dybxhewxfj6zp6hkagxv3z9kn6+jiwuzoznjyntgfxnea5trvuhaebpxktfxg3odspw6lpbtlvkyrfa3tnkanwnkgchkamcvcajsdgwwgya9xoflnzpqcpw2lyoifz2vbsnbzyzpsq8lfhx7wwwsfbbx+c2hxpnfnnym8vb2sfmw2l5/8/5pvpahbdauwqj223sfh4tnmfbfbst7mmh28uy+g5+3yjlwghg3pjd8oeckrfh5il8il/i74hfx6nvfcatry/xlki6+q822uv7</latexit> <latexit sha_base64="mys89qmsy7tmfxhlnwnn4jsy8=">aaac4xicbvhlbtnafj2yr4t5tbbkm2ahssgylss6gddskhnwym2ouvjjblhu7mnsi48whsebswbobh+a/+hebpcbpluy6ovfmob3vuudtp7yc6dfvo3z3de/h9bw8fpd7bf3litgmzhzejjt2whepnb+hqmnpastbvzkfvov3xf/e7dogh2my5qxcuzazaqddftpjm2bqi8pvs79zc9jh2lfdbtkk5cqvrn9gd/iqlhjeiamqtnxllay9mcrcek93fcni7xwbyw52nmorq49dybxhewxfj6zp6hkagxv3z9kn6+jiwuzoznjyntgfxnea5trvuhaebpxktfxg3odspw6lpbtlvkyrfa3tnkanwnkgchkamcvcajsdgwwgya9xoflnzpqcpw2lyoifz2vbsnbzyzpsq8lfhx7wwwsfbbx+c2hxpnfnnym8vb2sfmw2l5/8/5pvpahbdauwqj223sfh4tnmfbfbst7mmh28uy+g5+3yjlwghg3pjd8oeckrfh5il8il/i74hfx6nvfcatry/xlki6+q822uv7</latexit> <latexit sha_base64="mys89qmsy7tmfxhlnwnn4jsy8=">aaac4xicbvhlbtnafj2yr4t5tbbkm2ahssgylss6gddskhnwym2ouvjjblhu7mnsi48whsebswbobh+a/+hebpcbpluy6ovfmob3vuudtp7yc6dfvo3z3de/h9bw8fpd7bf3litgmzhzejjt2whepnb+hqmnpastbvzkfvov3xf/e7dogh2my5qxcuzazaqddftpjm2bqi8pvs79zc9jh2lfdbtkk5cqvrn9gd/iqlhjeiamqtnxllay9mcrcek93fcni7xwbyw52nmorq49dybxhewxfj6zp6hkagxv3z9kn6+jiwuzoznjyntgfxnea5trvuhaebpxktfxg3odspw6lpbtlvkyrfa3tnkanwnkgchkamcvcajsdgwwgya9xoflnzpqcpw2lyoifz2vbsnbzyzpsq8lfhx7wwwsfbbx+c2hxpnfnnym8vb2sfmw2l5/8/5pvpahbdauwqj223sfh4tnmfbfbst7mmh28uy+g5+3yjlwghg3pjd8oeckrfh5il8il/i74hfx6nvfcatry/xlki6+q822uv7</latexit> <latexit sha_base64="jpqnipstojkke9laqm9h7jz6r8=">aaac23icbvhlbhmxfhwmpmwasuszcaiiq3reirulhvswbapasplrtedxmt+dgypuxb9yodysocdfxl/6n/g2cape7skk6ovf4hn7q5ozy7psehdt3lv/4ohuo/jxk6fpnu/th5wawhcrrxpc8qmjqzsuewwu7pakbvjyoq+wntj++onowju/sqqalgivkcbabmpmps6+w366l2bdktkg+asng6jqf7g6tipkgjqlsegzetnnw2dkati5z6oc4aq2sgshqvwmpcvl5hivantfwtfmkp54ekefqu7pupj+8ccwsmssdnrjx/acqbizelvqtp5ms9esd/umjz/kb2tdwopjddb84ynvixtpiz5ryvgoaigzhtiscgwziwbjqtkvrakbcuaksoof4nivhoscxrf+5gswnuexxotun33lonpltkzs/wxyuxwj+27b49f9jvplutxdfuwat3in2pg4xx5v33qan+rbnq/zlfxcx37xfqsvuzveuah6ah9rsdohahaol/od/obldh36ef8yaddz/xqberb//azqj6w4=</latexit> <latexit sha_base64="jpqnipstojkke9laqm9h7jz6r8=">aaac23icbvhlbhmxfhwmpmwasuszcaiiq3reirulhvswbapasplrtedxmt+dgypuxb9yodysocdfxl/6n/g2cape7skk6ovf4hn7q5ozy7psehdt3lv/4ohuo/jxk6fpnu/th5wawhcrrxpc8qmjqzsuewwu7pakbvjyoq+wntj++onowju/sqqalgivkcbabmpmps6+w366l2bdktkg+asng6jqf7g6tipkgjqlsegzetnnw2dkati5z6oc4aq2sgshqvwmpcvl5hivantfwtfmkp54ekefqu7pupj+8ccwsmssdnrjx/acqbizelvqtp5ms9esd/umjz/kb2tdwopjddb84ynvixtpiz5ryvgoaigzhtiscgwziwbjqtkvrakbcuaksoof4nivhoscxrf+5gswnuexxotun33lonpltkzs/wxyuxwj+27b49f9jvplutxdfuwat3in2pg4xx5v33qan+rbnq/zlfxcx37xfqsvuzveuah6ah9rsdohahaol/od/obldh36ef8yaddz/xqberb//azqj6w4=</latexit> <latexit sha_base64="jpqnipstojkke9laqm9h7jz6r8=">aaac23icbvhlbhmxfhwmpmwasuszcaiiq3reirulhvswbapasplrtedxmt+dgypuxb9yodysocdfxl/6n/g2cape7skk6ovf4hn7q5ozy7psehdt3lv/4ohuo/jxk6fpnu/th5wawhcrrxpc8qmjqzsuewwu7pakbvjyoq+wntj++onowju/sqqalgivkcbabmpmps6+w366l2bdktkg+asng6jqf7g6tipkgjqlsegzetnnw2dkati5z6oc4aq2sgshqvwmpcvl5hivantfwtfmkp54ekefqu7pupj+8ccwsmssdnrjx/acqbizelvqtp5ms9esd/umjz/kb2tdwopjddb84ynvixtpiz5ryvgoaigzhtiscgwziwbjqtkvrakbcuaksoof4nivhoscxrf+5gswnuexxotun33lonpltkzs/wxyuxwj+27b49f9jvplutxdfuwat3in2pg4xx5v33qan+rbnq/zlfxcx37xfqsvuzveuah6ah9rsdohahaol/od/obldh36ef8yaddz/xqberb//azqj6w4=</latexit> <latexit sha_base64="jpqnipstojkke9laqm9h7jz6r8=">aaac23icbvhlbhmxfhwmpmwasuszcaiiq3reirulhvswbapasplrtedxmt+dgypuxb9yodysocdfxl/6n/g2cape7skk6ovf4hn7q5ozy7psehdt3lv/4ohuo/jxk6fpnu/th5wawhcrrxpc8qmjqzsuewwu7pakbvjyoq+wntj++onowju/sqqalgivkcbabmpmps6+w366l2bdktkg+asng6jqf7g6tipkgjqlsegzetnnw2dkati5z6oc4aq2sgshqvwmpcvl5hivantfwtfmkp54ekefqu7pupj+8ccwsmssdnrjx/acqbizelvqtp5ms9esd/umjz/kb2tdwopjddb84ynvixtpiz5ryvgoaigzhtiscgwziwbjqtkvrakbcuaksoof4nivhoscxrf+5gswnuexxotun33lonpltkzs/wxyuxwj+27b49f9jvplutxdfuwat3in2pg4xx5v33qan+rbnq/zlfxcx37xfqsvuzveuah6ah9rsdohahaol/od/obldh36ef8yaddz/xqberb//azqj6w4=</latexit> <latexit sha_base64="c949cdtaqfzwpdcjzzpbnfr7ju=">aaac4xicbvhlbtqwfpwert4tbdsxhahsweurhwcggzze6bavjnlxekbw+jhadmhi/qhdvwxysief4t/4mzpiduzvzklo3opz9g9t6wety5n/w2ie/cfpnzafhq/fvl2fod3rfhvtegshhvqpvteiwtxlgr46w8owkkvgj+xiu9s/+cqm5voduwxfcglzxwecggtuysdn7hayveo+mlokg7tvamicuqofudtnyhf/oppvkylebo5jwmiaem4fczhcv5bvgfdwjwfygqqtvsrlcwqxftcc99thabzlzoumg9phnykz4pk4yugo7tmbthypy6bspe6yxv6on9svgq6q2jffb4jmtcbo43zjemono4sawbqebgnzmwqf3ya5w9oqkufnm7w4mekahibai5yqi58x7bghncjcq3nczvohsa7lzjkttdsoumbttml/6/zfbt3fgn3m7acobi+n9hc5ou+7cy6ziumh5mv75odj6vbbaa+9rm8rqr/qafqmdteiuxsofqlf6e9eo8vokvpxi4gqz8vua+ix9caov6</latexit> <latexit sha_base64="c949cdtaqfzwpdcjzzpbnfr7ju=">aaac4xicbvhlbtqwfpwert4tbdsxhahsweurhwcggzze6bavjnlxekbw+jhadmhi/qhdvwxysief4t/4mzpiduzvzklo3opz9g9t6wety5n/w2ie/cfpnzafhq/fvl2fod3rfhvtegshhvqpvteiwtxlgr46w8owkkvgj+xiu9s/+cqm5voduwxfcglzxwecggtuysdn7hayveo+mlokg7tvamicuqofudtnyhf/oppvkylebo5jwmiaem4fczhcv5bvgfdwjwfygqqtvsrlcwqxftcc99thabzlzoumg9phnykz4pk4yugo7tmbthypy6bspe6yxv6on9svgq6q2jffb4jmtcbo43zjemono4sawbqebgnzmwqf3ya5w9oqkufnm7w4mekahibai5yqi58x7bghncjcq3nczvohsa7lzjkttdsoumbttml/6/zfbt3fgn3m7acobi+n9hc5ou+7cy6ziumh5mv75odj6vbbaa+9rm8rqr/qafqmdteiuxsofqlf6e9eo8vokvpxi4gqz8vua+ix9caov6</latexit> <latexit sha_base64="c949cdtaqfzwpdcjzzpbnfr7ju=">aaac4xicbvhlbtqwfpwert4tbdsxhahsweurhwcggzze6bavjnlxekbw+jhadmhi/qhdvwxysief4t/4mzpiduzvzklo3opz9g9t6wety5n/w2ie/cfpnzafhq/fvl2fod3rfhvtegshhvqpvteiwtxlgr46w8owkkvgj+xiu9s/+cqm5voduwxfcglzxwecggtuysdn7hayveo+mlokg7tvamicuqofudtnyhf/oppvkylebo5jwmiaem4fczhcv5bvgfdwjwfygqqtvsrlcwqxftcc99thabzlzoumg9phnykz4pk4yugo7tmbthypy6bspe6yxv6on9svgq6q2jffb4jmtcbo43zjemono4sawbqebgnzmwqf3ya5w9oqkufnm7w4mekahibai5yqi58x7bghncjcq3nczvohsa7lzjkttdsoumbttml/6/zfbt3fgn3m7acobi+n9hc5ou+7cy6ziumh5mv75odj6vbbaa+9rm8rqr/qafqmdteiuxsofqlf6e9eo8vokvpxi4gqz8vua+ix9caov6</latexit> <latexit sha_base64="c949cdtaqfzwpdcjzzpbnfr7ju=">aaac4xicbvhlbtqwfpwert4tbdsxhahsweurhwcggzze6bavjnlxekbw+jhadmhi/qhdvwxysief4t/4mzpiduzvzklo3opz9g9t6wety5n/w2ie/cfpnzafhq/fvl2fod3rfhvtegshhvqpvteiwtxlgr46w8owkkvgj+xiu9s/+cqm5voduwxfcglzxwecggtuysdn7hayveo+mlokg7tvamicuqofudtnyhf/oppvkylebo5jwmiaem4fczhcv5bvgfdwjwfygqqtvsrlcwqxftcc99thabzlzoumg9phnykz4pk4yugo7tmbthypy6bspe6yxv6on9svgq6q2jffb4jmtcbo43zjemono4sawbqebgnzmwqf3ya5w9oqkufnm7w4mekahibai5yqi58x7bghncjcq3nczvohsa7lzjkttdsoumbttml/6/zfbt3fgn3m7acobi+n9hc5ou+7cy6ziumh5mv75odj6vbbaa+9rm8rqr/qafqmdteiuxsofqlf6e9eo8vokvpxi4gqz8vua+ix9caov6</latexit> <latexit sha_base64="mselso9exuww44o+tlbpzgu4=">aaac3hicbvhlbhmxfhwgvxlelsy7myyqweuziyryvbhwasmjrqpr3hsaz4mbi9oddjhiu38c38b39tzzrinarxsnr7ve5uvewledwpem/qxtn737d/yexo8ep3n6bp/g+anvtafsrlxqzlyczyinnlcctaudanzcnzwj62/bovzfiuylbv6yqsfb8zim4qi3nlicoljddt9jh2hxebdkgjghtx9odwv8y7swtdkqwnpjllauama4tgxzcuxqyyqgkiwpgir8pt2pvacwbwxxvcu8cvcczlzpupo9fb2ag59qepxzu2j4tsgvxsgzktnu9ytt6kdvp3rcnvvtum6hxqvbbyadwuc8+4ydsjdqbadq+jyboea9sflczesw9uswlqpdsgj9krumeqivgsxbu+4ivc53gpmmi6tr+x6ht5ddnkvx2cvb2iflnq/ly/8d8o9z35b2auy2j8b7ofw8277vljjnhk6zd6/ty4+bg6/hw7rk/qgzegdokkfdeaiyoeukc/z/os/qj+hn9upzgg82ff6hxeuvvo/qhw==</latexit> <latexit sha_base64="mselso9exuww44o+tlbpzgu4=">aaac3hicbvhlbhmxfhwgvxlelsy7myyqweuziyryvbhwasmjrqpr3hsaz4mbi9oddjhiu38c38b39tzzrinarxsnr7ve5uvewledwpem/qxtn737d/yexo8ep3n6bp/g+anvtafsrlxqzlyczyinnlcctaudanzcnzwj62/bovzfiuylbv6yqsfb8zim4qi3nlicoljddt9jh2hxebdkgjghtx9odwv8y7swtdkqwnpjllauama4tgxzcuxqyyqgkiwpgir8pt2pvacwbwxxvcu8cvcczlzpupo9fb2ag59qepxzu2j4tsgvxsgzktnu9ytt6kdvp3rcnvvtum6hxqvbbyadwuc8+4ydsjdqbadq+jyboea9sflczesw9uswlqpdsgj9krumeqivgsxbu+4ivc53gpmmi6tr+x6ht5ddnkvx2cvb2iflnq/ly/8d8o9z35b2auy2j8b7ofw8277vljjnhk6zd6/ty4+bg6/hw7rk/qgzegdokkfdeaiyoeukc/z/os/qj+hn9upzgg82ff6hxeuvvo/qhw==</latexit> <latexit sha_base64="mselso9exuww44o+tlbpzgu4=">aaac3hicbvhlbhmxfhwgvxlelsy7myyqweuziyryvbhwasmjrqpr3hsaz4mbi9oddjhiu38c38b39tzzrinarxsnr7ve5uvewledwpem/qxtn737d/yexo8ep3n6bp/g+anvtafsrlxqzlyczyinnlcctaudanzcnzwj62/bovzfiuylbv6yqsfb8zim4qi3nlicoljddt9jh2hxebdkgjghtx9odwv8y7swtdkqwnpjllauama4tgxzcuxqyyqgkiwpgir8pt2pvacwbwxxvcu8cvcczlzpupo9fb2ag59qepxzu2j4tsgvxsgzktnu9ytt6kdvp3rcnvvtum6hxqvbbyadwuc8+4ydsjdqbadq+jyboea9sflczesw9uswlqpdsgj9krumeqivgsxbu+4ivc53gpmmi6tr+x6ht5ddnkvx2cvb2iflnq/ly/8d8o9z35b2auy2j8b7ofw8277vljjnhk6zd6/ty4+bg6/hw7rk/qgzegdokkfdeaiyoeukc/z/os/qj+hn9upzgg82ff6hxeuvvo/qhw==</latexit> <latexit sha_base64="mselso9exuww44o+tlbpzgu4=">aaac3hicbvhlbhmxfhwgvxlelsy7myyqweuziyryvbhwasmjrqpr3hsaz4mbi9oddjhiu38c38b39tzzrinarxsnr7ve5uvewledwpem/qxtn737d/yexo8ep3n6bp/g+anvtafsrlxqzlyczyinnlcctaudanzcnzwj62/bovzfiuylbv6yqsfb8zim4qi3nlicoljddt9jh2hxebdkgjghtx9odwv8y7swtdkqwnpjllauama4tgxzcuxqyyqgkiwpgir8pt2pvacwbwxxvcu8cvcczlzpupo9fb2ag59qepxzu2j4tsgvxsgzktnu9ytt6kdvp3rcnvvtum6hxqvbbyadwuc8+4ydsjdqbadq+jyboea9sflczesw9uswlqpdsgj9krumeqivgsxbu+4ivc53gpmmi6tr+x6ht5ddnkvx2cvb2iflnq/ly/8d8o9z35b2auy2j8b7ofw8277vljjnhk6zd6/ty4+bg6/hw7rk/qgzegdokkfdeaiyoeukc/z/os/qj+hn9upzgg82ff6hxeuvvo/qhw==</latexit> <latexit sha_base64="j/erh+nujcgcuc4ckwgxr+epxc=">aaac2xicbvhlatwwfnw49r9je2yg7wmnvgmczdo2myxqysca2w7vgmxgjh5hklkmirxclmyy6ax+m/9g/qexmiz7jbchh3knzupfwdwfgpunfuxtv/oohj3yex+epnv+ynfv5ylrszqhrx+qwgqzmtdgkz5fssrrezelpvfzu9u8vqdzmywo7amglychznbgwgzqqqbn+upuk47qvva2ynujquo6me6m/5uyrvlbpcqdjiixtbovawy49xfctoy2qjawok4bhskk8gnwgdv4jvyekauapqgqklcp5zh7wizw3olw5mw9+zacyqxkehzedpnnsdevevao38ohjmnq2lktwezvuocldpvcmauosxwuarlmwgibnoihysm+4lpquukanlmyuclhilb7ngnmku/vcwplyxxczzqxun33lonfltkys/wxyufwj+27b88/9jvplutxdltwat3lh2pg4xzzbvuwo8ngwjmv+8nhx/xtd9b9ba9rxn6ga7rz3sejogghq7rl/q7kqlvy/o6kyajdz/xqfbrt//av7b6pw=</latexit> <latexit sha_base64="j/erh+nujcgcuc4ckwgxr+epxc=">aaac2xicbvhlatwwfnw49r9je2yg7wmnvgmczdo2myxqysca2w7vgmxgjh5hklkmirxclmyy6ax+m/9g/qexmiz7jbchh3knzupfwdwfgpunfuxtv/oohj3yex+epnv+ynfv5ylrszqhrx+qwgqzmtdgkz5fssrrezelpvfzu9u8vqdzmywo7amglychznbgwgzqqqbn+upuk47qvva2ynujquo6me6m/5uyrvlbpcqdjiixtbovawy49xfctoy2qjawok4bhskk8gnwgdv4jvyekauapqgqklcp5zh7wizw3olw5mw9+zacyqxkehzedpnnsdevevao38ohjmnq2lktwezvuocldpvcmauosxwuarlmwgibnoihysm+4lpquukanlmyuclhilb7ngnmku/vcwplyxxczzqxun33lonfltkys/wxyufwj+27b88/9jvplutxdltwat3lh2pg4xzzbvuwo8ngwjmv+8nhx/xtd9b9ba9rxn6ga7rz3sejogghq7rl/q7kqlvy/o6kyajdz/xqfbrt//av7b6pw=</latexit> <latexit sha_base64="j/erh+nujcgcuc4ckwgxr+epxc=">aaac2xicbvhlatwwfnw49r9je2yg7wmnvgmczdo2myxqysca2w7vgmxgjh5hklkmirxclmyy6ax+m/9g/qexmiz7jbchh3knzupfwdwfgpunfuxtv/oohj3yex+epnv+ynfv5ylrszqhrx+qwgqzmtdgkz5fssrrezelpvfzu9u8vqdzmywo7amglychznbgwgzqqqbn+upuk47qvva2ynujquo6me6m/5uyrvlbpcqdjiixtbovawy49xfctoy2qjawok4bhskk8gnwgdv4jvyekauapqgqklcp5zh7wizw3olw5mw9+zacyqxkehzedpnnsdevevao38ohjmnq2lktwezvuocldpvcmauosxwuarlmwgibnoihysm+4lpquukanlmyuclhilb7ngnmku/vcwplyxxczzqxun33lonfltkys/wxyufwj+27b88/9jvplutxdltwat3lh2pg4xzzbvuwo8ngwjmv+8nhx/xtd9b9ba9rxn6ga7rz3sejogghq7rl/q7kqlvy/o6kyajdz/xqfbrt//av7b6pw=</latexit> <latexit sha_base64="j/erh+nujcgcuc4ckwgxr+epxc=">aaac2xicbvhlatwwfnw49r9je2yg7wmnvgmczdo2myxqysca2w7vgmxgjh5hklkmirxclmyy6ax+m/9g/qexmiz7jbchh3knzupfwdwfgpunfuxtv/oohj3yex+epnv+ynfv5ylrszqhrx+qwgqzmtdgkz5fssrrezelpvfzu9u8vqdzmywo7amglychznbgwgzqqqbn+upuk47qvva2ynujquo6me6m/5uyrvlbpcqdjiixtbovawy49xfctoy2qjawok4bhskk8gnwgdv4jvyekauapqgqklcp5zh7wizw3olw5mw9+zacyqxkehzedpnnsdevevao38ohjmnq2lktwezvuocldpvcmauosxwuarlmwgibnoihysm+4lpquukanlmyuclhilb7ngnmku/vcwplyxxczzqxun33lonfltkys/wxyufwj+27b88/9jvplutxdltwat3lh2pg4xzzbvuwo8ngwjmv+8nhx/xtd9b9ba9rxn6ga7rz3sejogghq7rl/q7kqlvy/o6kyajdz/xqfbrt//av7b6pw=</latexit> <latexit sha_base64="fhlzcdekggu24jrugggsnptok=">aaac23icbvfnb9naen24qiv5asuri8fc4htzviu4vndhwksmqrrbxgzsvfzd2t3trxcpxfdxdjuv4l/4n/w9onuppjwe3x9tzntvpwzmyr/b8hog4epdvceh+epnv+yv/g8myowlmcucwvpi/bigcsr5zzjuevrhalx3g5/nt2x9rg6bkqvvwahysdznfkynxnqal2hhuh8nw6sabukaxctdzmdwz/8pmitubpkqdjjmls2aibbrnl6miwwwqjewwkyc7uokcjwgl2o4jtyuu898vccqfm3xgueuuzwtrx2j9po47toyewzivk2w9zwavje/two7/atfarulvg+yj6yk2lfm6arw7yakhmfsixoagavgwzijvvcgjweanbtdkiytnibcc4vxj9wrjtj7ikx3gflgmyuyzxd/yk7ve7k5h/v/x2wz3wx5octghotxohv3i6ed9tcjyn2syfjmkjz+ub79hxpe35bjyxtytd6teziilczjl3jdfgdf8d34efy8lqad9z+xpffb9t8pyom</latexit> <latexit sha_base64="fhlzcdekggu24jrugggsnptok=">aaac23icbvfnb9naen24qiv5asuri8fc4htzviu4vndhwksmqrrbxgzsvfzd2t3trxcpxfdxdjuv4l/4n/w9onuppjwe3x9tzntvpwzmyr/b8hog4epdvceh+epnv+yv/g8myowlmcucwvpi/bigcsr5zzjuevrhalx3g5/nt2x9rg6bkqvvwahysdznfkynxnqal2hhuh8nw6sabukaxctdzmdwz/8pmitubpkqdjjmls2aibbrnl6miwwwqjewwkyc7uokcjwgl2o4jtyuu898vccqfm3xgueuuzwtrx2j9po47toyewzivk2w9zwavje/two7/atfarulvg+yj6yk2lfm6arw7yakhmfsixoagavgwzijvvcgjweanbtdkiytnibcc4vxj9wrjtj7ikx3gflgmyuyzxd/yk7ve7k5h/v/x2wz3wx5octghotxohv3i6ed9tcjyn2syfjmkjz+ub79hxpe35bjyxtytd6teziilczjl3jdfgdf8d34efy8lqad9z+xpffb9t8pyom</latexit> <latexit sha_base64="fhlzcdekggu24jrugggsnptok=">aaac23icbvfnb9naen24qiv5asuri8fc4htzviu4vndhwksmqrrbxgzsvfzd2t3trxcpxfdxdjuv4l/4n/w9onuppjwe3x9tzntvpwzmyr/b8hog4epdvceh+epnv+yv/g8myowlmcucwvpi/bigcsr5zzjuevrhalx3g5/nt2x9rg6bkqvvwahysdznfkynxnqal2hhuh8nw6sabukaxctdzmdwz/8pmitubpkqdjjmls2aibbrnl6miwwwqjewwkyc7uokcjwgl2o4jtyuu898vccqfm3xgueuuzwtrx2j9po47toyewzivk2w9zwavje/two7/atfarulvg+yj6yk2lfm6arw7yakhmfsixoagavgwzijvvcgjweanbtdkiytnibcc4vxj9wrjtj7ikx3gflgmyuyzxd/yk7ve7k5h/v/x2wz3wx5octghotxohv3i6ed9tcjyn2syfjmkjz+ub79hxpe35bjyxtytd6teziilczjl3jdfgdf8d34efy8lqad9z+xpffb9t8pyom</latexit> <latexit sha_base64="fhlzcdekggu24jrugggsnptok=">aaac23icbvfnb9naen24qiv5asuri8fc4htzviu4vndhwksmqrrbxgzsvfzd2t3trxcpxfdxdjuv4l/4n/w9onuppjwe3x9tzntvpwzmyr/b8hog4epdvceh+epnv+yv/g8myowlmcucwvpi/bigcsr5zzjuevrhalx3g5/nt2x9rg6bkqvvwahysdznfkynxnqal2hhuh8nw6sabukaxctdzmdwz/8pmitubpkqdjjmls2aibbrnl6miwwwqjewwkyc7uokcjwgl2o4jtyuu898vccqfm3xgueuuzwtrx2j9po47toyewzivk2w9zwavje/two7/atfarulvg+yj6yk2lfm6arw7yakhmfsixoagavgwzijvvcgjweanbtdkiytnibcc4vxj9wrjtj7ikx3gflgmyuyzxd/yk7ve7k5h/v/x2wz3wx5octghotxohv3i6ed9tcjyn2syfjmkjz+ub79hxpe35bjyxtytd6teziilczjl3jdfgdf8d34efy8lqad9z+xpffb9t8pyom</latexit> <latexit sha_base64="bokjaclam5dgjgxnc5cbux+3rl4=">aaac33icbvfnb9naen2y2k+wjhywbcqoev2vamofvw4fqlpk9lwnn5slx2w9odg4k7z26icwcu8e/4h/wbm6q6qqjft5u7mpmqlsjhmv4drbdu37l7b+9+/odho8dp9g+enjtwmanzehjzytwxanpyhq8vpaclcv5gfv6n3xp/verngn+c65qwcprylwqadlbtzw6dqa4p+tp+k47qvuguydujipo5nb6m/xdywrngntijzezbwwlzgutdjfrwxjemsbusevuddjha+agac9q+gi8hajzaduo7sq2h9htv4gz4wx4wmkpttwauxcwlvb2xm67v5h3ttlgy8lvuh6wa5zldbizsnltbf5lyxnkdqdaaijuxybfhigcaf5p+zuqvcqcz7pyaqojesj9mlhwpwhhvbzzivttf4hydem7lukkxufgdlnzi5bli8+p9xsp2odftwc2hpfrwunm3fdxectszzos4+hizh7za33yppyuvymmtkdtkih8gxmrjgdplbfphfeurfo2/r9ytpnn8euygff38b/pz64g=</latexit> <latexit sha_base64="bokjaclam5dgjgxnc5cbux+3rl4=">aaac33icbvfnb9naen2y2k+wjhywbcqoev2vamofvw4fqlpk9lwnn5slx2w9odg4k7z26icwcu8e/4h/wbm6q6qqjft5u7mpmqlsjhmv4drbdu37l7b+9+/odho8dp9g+enjtwmanzehjzytwxanpyhq8vpaclcv5gfv6n3xp/verngn+c65qwcprylwqadlbtzw6dqa4p+tp+k47qvuguydujipo5nb6m/xdywrngntijzezbwwlzgutdjfrwxjemsbusevuddjha+agac9q+gi8hajzaduo7sq2h9htv4gz4wx4wmkpttwauxcwlvb2xm67v5h3ttlgy8lvuh6wa5zldbizsnltbf5lyxnkdqdaaijuxybfhigcaf5p+zuqvcqcz7pyaqojesj9mlhwpwhhvbzzivttf4hydem7lukkxufgdlnzi5bli8+p9xsp2odftwc2hpfrwunm3fdxectszzos4+hizh7za33yppyuvymmtkdtkih8gxmrjgdplbfphfeurfo2/r9ytpnn8euygff38b/pz64g=</latexit> <latexit sha_base64="bokjaclam5dgjgxnc5cbux+3rl4=">aaac33icbvfnb9naen2y2k+wjhywbcqoev2vamofvw4fqlpk9lwnn5slx2w9odg4k7z26icwcu8e/4h/wbm6q6qqjft5u7mpmqlsjhmv4drbdu37l7b+9+/odho8dp9g+enjtwmanzehjzytwxanpyhq8vpaclcv5gfv6n3xp/verngn+c65qwcprylwqadlbtzw6dqa4p+tp+k47qvuguydujipo5nb6m/xdywrngntijzezbwwlzgutdjfrwxjemsbusevuddjha+agac9q+gi8hajzaduo7sq2h9htv4gz4wx4wmkpttwauxcwlvb2xm67v5h3ttlgy8lvuh6wa5zldbizsnltbf5lyxnkdqdaaijuxybfhigcaf5p+zuqvcqcz7pyaqojesj9mlhwpwhhvbzzivttf4hydem7lukkxufgdlnzi5bli8+p9xsp2odftwc2hpfrwunm3fdxectszzos4+hizh7za33yppyuvymmtkdtkih8gxmrjgdplbfphfeurfo2/r9ytpnn8euygff38b/pz64g=</latexit> <latexit sha_base64="bokjaclam5dgjgxnc5cbux+3rl4=">aaac33icbvfnb9naen2y2k+wjhywbcqoev2vamofvw4fqlpk9lwnn5slx2w9odg4k7z26icwcu8e/4h/wbm6q6qqjft5u7mpmqlsjhmv4drbdu37l7b+9+/odho8dp9g+enjtwmanzehjzytwxanpyhq8vpaclcv5gfv6n3xp/verngn+c65qwcprylwqadlbtzw6dqa4p+tp+k47qvuguydujipo5nb6m/xdywrngntijzezbwwlzgutdjfrwxjemsbusevuddjha+agac9q+gi8hajzaduo7sq2h9htv4gz4wx4wmkpttwauxcwlvb2xm67v5h3ttlgy8lvuh6wa5zldbizsnltbf5lyxnkdqdaaijuxybfhigcaf5p+zuqvcqcz7pyaqojesj9mlhwpwhhvbzzivttf4hydem7lukkxufgdlnzi5bli8+p9xsp2odftwc2hpfrwunm3fdxectszzos4+hizh7za33yppyuvymmtkdtkih8gxmrjgdplbfphfeurfo2/r9ytpnn8euygff38b/pz64g=</latexit> <latexit sha_base64="iwcvmymum7wvdqvstgdze8pl9b=">aaac23icbvhlahsxfjunfattv5iuu5l2kbqkzjqumvon2memcbz2duylijmcgavjcravusjdddnp+s/4jfpnxiwmnqucw7lh53dvwojm2y6gc+/+g4e7j+lht54+e763f3bqvkmjhrhflt6wfdojbzzjk9qzufuxe6paf2v74gmdldyxq5qwahaszrkbg6gxmt7dvvpxpons66sbydxiexop7ud66kmsknonisdszmcfbbog2jhdq47hodk2blgfbxq8xenle6wae7dn+wtnnosoarbtem68xzyjjczzk5enimhdtzamhmslrb2xqazv5l3twbnhb+oxrmo2lktwezruewjwupmtfni+soaijqfrjydhqidruncm/eiueyjkkg+gktxcjaltln86fucjbwd4dlfhe4fsuh+s3tdl8bnfqdon5bcpif+p+wa6d4v7qkccfa+zhchg/edxuc5kocdfgx9+nrx/xtd9fl9bq9rrgdoipgr2jesjoix6hp+hvvebfox/rzxtpnfj/eyf6ff3+b48j6ww=</latexit> <latexit sha_base64="iwcvmymum7wvdqvstgdze8pl9b=">aaac23icbvhlahsxfjunfattv5iuu5l2kbqkzjqumvon2memcbz2duylijmcgavjcravusjdddnp+s/4jfpnxiwmnqucw7lh53dvwojm2y6gc+/+g4e7j+lht54+e763f3bqvkmjhrhflt6wfdojbzzjk9qzufuxe6paf2v74gmdldyxq5qwahaszrkbg6gxmt7dvvpxpons66sbydxiexop7ud66kmsknonisdszmcfbbog2jhdq47hodk2blgfbxq8xenle6wae7dn+wtnnosoarbtem68xzyjjczzk5enimhdtzamhmslrb2xqazv5l3twbnhb+oxrmo2lktwezruewjwupmtfni+soaijqfrjydhqidruncm/eiueyjkkg+gktxcjaltln86fucjbwd4dlfhe4fsuh+s3tdl8bnfqdon5bcpif+p+wa6d4v7qkccfa+zhchg/edxuc5kocdfgx9+nrx/xtd9fl9bq9rrgdoipgr2jesjoix6hp+hvvebfox/rzxtpnfj/eyf6ff3+b48j6ww=</latexit> <latexit sha_base64="iwcvmymum7wvdqvstgdze8pl9b=">aaac23icbvhlahsxfjunfattv5iuu5l2kbqkzjqumvon2memcbz2duylijmcgavjcravusjdddnp+s/4jfpnxiwmnqucw7lh53dvwojm2y6gc+/+g4e7j+lht54+e763f3bqvkmjhrhflt6wfdojbzzjk9qzufuxe6paf2v74gmdldyxq5qwahaszrkbg6gxmt7dvvpxpons66sbydxiexop7ud66kmsknonisdszmcfbbog2jhdq47hodk2blgfbxq8xenle6wae7dn+wtnnosoarbtem68xzyjjczzk5enimhdtzamhmslrb2xqazv5l3twbnhb+oxrmo2lktwezruewjwupmtfni+soaijqfrjydhqidruncm/eiueyjkkg+gktxcjaltln86fucjbwd4dlfhe4fsuh+s3tdl8bnfqdon5bcpif+p+wa6d4v7qkccfa+zhchg/edxuc5kocdfgx9+nrx/xtd9fl9bq9rrgdoipgr2jesjoix6hp+hvvebfox/rzxtpnfj/eyf6ff3+b48j6ww=</latexit> <latexit sha_base64="iwcvmymum7wvdqvstgdze8pl9b=">aaac23icbvhlahsxfjunfattv5iuu5l2kbqkzjqumvon2memcbz2duylijmcgavjcravusjdddnp+s/4jfpnxiwmnqucw7lh53dvwojm2y6gc+/+g4e7j+lht54+e763f3bqvkmjhrhflt6wfdojbzzjk9qzufuxe6paf2v74gmdldyxq5qwahaszrkbg6gxmt7dvvpxpons66sbydxiexop7ud66kmsknonisdszmcfbbog2jhdq47hodk2blgfbxq8xenle6wae7dn+wtnnosoarbtem68xzyjjczzk5enimhdtzamhmslrb2xqazv5l3twbnhb+oxrmo2lktwezruewjwupmtfni+soaijqfrjydhqidruncm/eiueyjkkg+gktxcjaltln86fucjbwd4dlfhe4fsuh+s3tdl8bnfqdon5bcpif+p+wa6d4v7qkccfa+zhchg/edxuc5kocdfgx9+nrx/xtd9fl9bq9rrgdoipgr2jesjoix6hp+hvvebfox/rzxtpnfj/eyf6ff3+b48j6ww=</latexit> 3 s t 2 c t a s t c t b s t c t+ The ates of fluoescence and backgound counts, λ and µ can be added then emoved by maginalisation, P (α β c I = dλ dµ P (α β I P ( c λ µ α β I P ( c I (6 = dλ dµp (α β λ µ I P ( c Ω. P ( c I (7 l µ FIG. 4. A patial Bayesian netwok epesenting the joint pobability distibution of poblem paametes. The patten of nodes and connections inside the cental gayed egion epesenting detection inteval t ae epeated fo each data value. We ae estimating paametes α, β, λ, and µ given the obseved counts c t. Fo the discete-time models α and β epesent the pobabilities α d and β d espectively. Enough of the full netwok is dawn to be able to easily detemine the vaiable independencies. In the case of a single step ove the detecto inteval, s t is the state ove the entie inteval and the dashed aows ae absent. In all othe cases s t and s t ae the states at the boundaies of the detecto inteval. the switching pobabilities become dependent on the entie histoy of the count data. This makes the infeence consideably moe involved. We can summaize the dependencies amongst the vaiables by the Bayesian netwok (BN [8] shown in Fig. 4. The BN can be used to detemine conditional independencies between the poblem vaiables, using the popety of d-sepaability. In the BN given in Fig. 4, s t and s t epesents the state of the emitte at the boundaies of the t th detecto inteval, but fo the single-step model I the initial state extends fo the whole duation so the dashed aows in the figue should be omitted. That is Eq. ( depends only on s t (this will change fo multistep models. The key vaiable independencies ae the following α β λ µ I (2 s t s u s t Ω (u < t (3 c t c u, s v s t Ω (u t, v t (4 whee Ω = α β λ µ I fo compactness. We will make use of these independencies in the following section. B. Bayesian Infeence Ultimately, the quantity we want to infe is P (α β c I. Using Bayes ule, the posteio pobability distibution of α and β can be witten as, P (α β c I = P (α β IP ( c α β I P ( c I (5 Note that if µ and λ ae known then the integals in (7 ae unnecessay and this will be tue in all the models we develop. Given (2 we can facto P (α β λ µ I = P (α I P (β I P (λ I P (µ I. We will go futhe and take them all as constant ove some initial ange so that P (α β c I = dλ dµp ( c Ω (8 N and with the nomalisation facto N to be detemined at the end. As stated ealie, if all the states s s s... s n whee known togethe with Ω the pobability of all the counts would be easily detemined. We can use the same tick, adding these paametes then maginalising ove them so that, P ( c Ω = s P ( c s Ω = s P ( c s Ω P ( s Ω (9 whee s = s... s n. Using the independencies (4, P ( c s Ω can be simplified to P ( c s Ω = P (c t s t Ω, ( t= and each tem is detemined by (. The emaining P ( s Ω tem cannot be simply factoised ove the states despite being a Makov chain, as obseving Ω intoduces possible dependencies. We can howeve expand using the poduct ule and simplify by making use of the independencies in (3: P ( s Ω = P (s t s t Ω P (s Ω, ( t= whee we ve chosen to expand in tempoal ode. Finally, the infeence becomes P (α β c I = dλ dµ P (c t s t s t Ω N s t= P (s Ω (2 whee we ve used (4 to wite the joint distibution between c t and s t. The poblem with (2 is that the sum ove s contains 2 N tems each of which has N poducts. This will apidly

4 99% 4 become intactable as the size of the data gows. Fotunately it s possible to ewite (2 as a single tem with N 2 2 matix poducts. Conside the following matix: [ ] P (ct s R t = t = s t = Ω P (c t s t = s t = Ω, P (c t s t = s t = Ω P (c t s t = s t = Ω (3 and vecto then s t= D = [ ] P (s = Ω P (s = Ω (4 P (c t s t s t Ω P (s Ω = [ ] R td (5 t= whee the poduct on the RHS is in deceasing t to the ight. The equivalence in (5 is eadily veified by expanding a few tems out. The matix multiplication will sum ove the columns of R t which is a summation ove the states s t. Each successive multiplication sums ove anothe state and the final multiplication by the ow vecto [ ] will sum ove s N. With this equivalence, the infeence eads P (α β c I = N dλ dµ[ ] R td, (6 t= which is easily computable. Note that cae must be taken to avoid undeflow o oveflow in calculating the matix poducts. As a demonstation of the algoithm, Fig. 5 (a shows the cedible egions of α and β fo the data given in Fig. 2 (a. The maginal distibutions of all unknown paametes (α, β, λ and µ ae given in Fig. 5 (b. It should be noted that the infeence contains no appoximations and makes use of all the data. The key point of contention with obseved data is if the model is ealistic to the system. If it is, then the infeence faithfully conveges on the undelying values egadless of what they ae. Fo example in Fig. 6 we geneated count data assuming high pobabilities of switching (α =.8 and β =.9 so the data looks vey unlike the usual blinking tace. The infeence quickly conveges on the coect pobabilities despite this (the 99% cedible egion is given. In the following sections we develop models that assume the emitte is fee to switch states within a detection window, aguably a moe geneal situation. III. CTMC MODEL In the pevious section the emitte could only switch states at the boundaies of the detection inteval. Though this is a hidden Makov model, the infeence was simple as the state is only obfuscated by the Poissonian distibution of counts. In this section we go to the (b.2 P( P(..2. (a % 9% P( P( FIG. 5. (a Infeence of the switching pobabilities on simulated data given in Fig. 2 (a. The contous contain the 5%, 9%, and 99% cedible egions. The tue value is shown as ed dot. (b Pobability distibutions of α and β and maginal distibutions of λ and µ. The tue values of the unknown paametes ae given in ed ba. othe exteme and assume the emitte can change state at any point in time, possibly multiple times in a single detection inteval. The pevious single-inteval model is a good appoximation to the continuous time model in the limit that the switching ates ae vey small. In that limit the occasional switch duing a detecto inteval intoduces negligible eo and most intevals ae eithe entiely on o entiely off. A. Model The emitte can switch at any time, so that the state is a function of time, i.e. s(t. This situation is modeled by a continuous time Makov chain (CTMC as depicted on the ight of Fig. 3. At any given time, the state can still only have one of two values but can switch abitaily often in any time inteval. The detection events ae as befoe with the detectos epoting the accumulated count fom a time window of constant duation T. Fo claity we shall take time to be expessed in units of the detecto inteval so that T =. In a continuous-time Makov chain the time spent in any state is exponentially distibuted, and since thee ae only two states, the next state is deteministic. The

5 5 (a 2 Counts (b Time step Fequency B. Infeence The basic infeence is as befoe, we have a list of n detecto counts c and we want to infe α and β while optionally maginalising ove the ates µ and λ leading an equivalent expession to Eq. (23 but fo α and β : P ( α β c I c = dλ dµp ( c Ω c (23 N % 99% FIG. 6. (a Time tace with high switching pobability of α =.8 and β =.9, simulated using the DTMC single-step model. (b The posteio distibution fo switching pobabilities infeed using single inteval model. chaacteistic lifetime of the off-intevals will be invesely popotional to the ate of switching out of the off state which will be the switch-on ate α. Similaly the onintevals ae invesely popotional to the switch-off ate β. In this model, these two paametes α and β, ae the paametes we want to infe, eplacing α and β. The CTMC can be solved by the fowad Kolmogoov equations. Repesenting the off state by the vecto [ ] T and the on state by [ ] T the tansitions dynamics ae encoded by the equation whee dp (t dt Q = = QP (t (7 [ ] α α. (8 β β This has the fomal solution P (t = exp(qt. Fom this solution we can detemine the tansition pobabilities P ab (t fom state a to state b in time t: P (t = P (t = P (t = P (t = ( β + α e (α+ βt, (9 α + β ( α e (α+ βt, (2 α + β ( β e (α+ βt, (2 α + β ( α + β e (α+ βt. (22 α + β whee Ω c = α β λ µ I c simila to peviously, and the subscipt denotes the continuous-time model, and N is suitably edefined. To get taction on the poblem, we add in the states at the boundaies of the detection inteval s j = s(t j, whee the n th detection inteval goes fom t n to t n and t n = n since T =. P ( c Ω c = s P ( c s Ω c P ( s Ω c (24 Note that in this situation the states s j at the boundaies ae not enough to specify the state duing the detection inteval, unlike in the pevious case whee the state was constant duing the whole inteval. The indepencencies ae now and we can simplify P ( c Ω c = s = s c n c m s n s n Ω c m n (25 s n s m s n Ω c m < n (26 P (c n s n s n Ω c P (s n s n Ω c P (s Ω c n= (27 P (c n s n s n Ω c P (s Ω c. (28 n= This expession can again be conveted into an efficient matix equation identical in fom to ealie: P (c t s n s n Ω c P (s Ω c = [ ] R nd s n= n= (29 with [ ] P (cn s R n = n = s n = Ω c P (c n s n = s n = Ω c, P (c n s n = s n = Ω c P (c n s n = s n = Ω c (3 and D = [ ] P (s = Ω c P (s = Ω c (3 All that emains is to calculate P (c n s n s n Ω c. If we knew f, the faction of the inteval fo which the emitte

6 was in the on state, then the counts would be given by a weighted Poissonian. That is, P (c n s n s n Ω c = and dfp (c n s n s n f Ω c P (s n f s n Ω c, (32 P (c n s n s n f Ω c = P(c n µ + fλ. (33 Now, we need to find an expession fo P (s n f s n Ω c. Fo bevity, We ll intoduce the notation R ab (f P (s n = b f s n = a Ω c. This pobability accounts fo all the possible histoies whee the detection inteval spent a faction f in state on and ended in the state s n given it stated in state s n. The diffeent histoies can be decomposed by diffeing numbes of switch events. Conside the case whee the state stats and ends in the off state, i.e. R (f. Clealy an odd numbe of switch events will not be consistent with the bounday states and hence will have pobability zeo. With zeo switch events the entie inteval is off and the pobability is just the exponential distibution δ(fe α, whee δ is the Diac delta function. Two switch events patitions the inteval into thee egions l, l 2, and l 3 with states off on off and pobability Summing ove all switch events gives R (f =δ(fe α + e β( f β f f k ( f k (k!k! α k β k (37 k= =δ(fe α e α( f (f βf α β f J (2 (f f α β (38 whee J (x is the Bessel function of the fist kind. A simila agument yields the othe possible faction pobabilities: R (f = α J (2 (f f α β e α( f βf (39 R (f = β J (2 (f f α β e α( f βf (4 R (f =δ( fe β + e α( f f βf α β f J (2 (f f α β. (4 So in summay we have 6 f dl e αl α e βl 2 β e αl3 = α β e α( f βf (34 whee we have used l 2 = f and l +l 3 = f. In geneal the inteval will be patitioned into a set of on states of total duation f and a set of off states of duation f so that the exponentials combine and will not depend on the duations l j. Fo example with six switches thee ae seven egions with thee switch-on and thee switch-off events, and the pobability is 3 α 3 βe α( f βf f f l2 dl 2 dl 4 }{{} f A f l f l 3 dl dl 3 dl 5 (35 } {{ } B The integals in A calculate the volume of a twodimensional simplex of side f, and those in B the volume of a thee dimensional simplex of side f. Since the volume of an n-dimensional simplex of side f is f n /n!, Eq. (35 evaluates to f 2 ( f 3 α 3 β 3e α( f βf. (36 2!3! P (c n s n s n Ω c = dfp(c n µ + fλr sn s n (f (42 which is an integal that needs to be done numeically. The data in Fig. 2 (b was used as an example of the infeence and is pesented in Fig. 7. The initial ange of paametes was assumed to be equally likely anywhee in α, β 8, and λ, µ 8, and a 4-dimensional gid whee each ange was divided into 7 values was used in evaluating the posteio pobability in (29 fo each data point to allow fo tempoay nomalization and contol of undeflow. C. Relation between single inteval model and continuous-time model The single inteval model and the continuous time model make diffeent assumptions about the undelying switching mechanism, nevetheless if the switching ates α and β ae low, we would expect the single inteval model to yield vey simila esults. This is because the detecto intevals that contain intenal switching events will be elatively ae and will not contibute significantly to the infeence. Additionally, the discete time model assumes that the emitte emains in a constant state ove each detecto inteval, which is a easonable appoximation fo small switching ates. It might be supposed that, if the ates ae small, the single inteval model could be used as a moe computationally efficient appoximation

7 7 (b. P( P(..2. (a P( P( FIG. 7. (a Infeence of the switching ates on simulated data given in Fig. 2 (b. The contous contain the 5%, 9%, and 99% cedible egions. The tue value is shown as ed dot. (b Pobability distibutions of α and β and maginal distibutions of λ and µ. The tue values of the unknown paametes ae given in ed ba. to the full infeence. We now exploe the validity of this appoximation. The two models can be connected by the time spent in a given state. Conside a time inteval t = T/d which is a faction of the measuement time. In the continuoustime model the pobability of a state having a lifetime t is exponentially distibuted, e.g P ( t P (s(t = : t t = exp( α t. So fo a model whee the state is fixed fo duation t, we can elate the pobability of switching state to the pobability of the state not having lifetime t, e.g. α d = P ( t. Consequently, α d = e αt/d (43 α d = e βt/d (44 Assuming small switching ates in the continuous model we can expand to fist ode in the switching ates. Fo Bessel function of the fist kind, ( z m J m (z z (45 Γ(m + 2 whee, Γ(x is the Gamma function. Hence at a low ate of switching, we can expand fist ode in α, β to obtain R (f δ(f( α, R (f δ( f( β, R (f α and R (f β. Futhemoe, fom (43 and (44 with d =, α α and β β. With this appoximations (42 yields, whee, P (c n Ω c ( α P(c n µ (46 P (c n Ω c ( β P(c n µ + λ (47 P (c n Ω c α c n (48 P (c n Ω c β c n (49 c n = dfp(c n µ + fλ, (5 is the pobability of obtaining c n counts aveaged ove all duations spent on in the inteval. The fist two pobabilities ae identical to those obtained in the single inteval model. The last two diffe, as in the single inteval model we made the abitay decision to take the switch event as happening at the ends of the intevals, wheeas in the coespondence between models we allow the switch to occu anywhee. Howeve, detecto intevals that contain a switch event become ae when the ates ae small. So the continuous time model becomes indistinguishable fom the single step model in this limit as expected. The single inteval model will loses accuacy as the ate of switching inceases and this can be exploed numeically. In Fig. 8 we summaise such an investigation by plotting the diffeence between the tue switching ate (fo data simulated fom a continuous-switching model and the esult of a single-inteval model infeence. The colou map epesents the Eucledian distance fom the tue value to the maximum value fom the posteio pobability distibution. The black egion coesponds to a small distance, and is theefoe the most applicable egion fo the single-inteval infeence. As the (continuous ates incease the single-inteval model stats failing. An appoximate theshold fo applicability of the single-step model is when α β <., which coesponds oughly to a elative eo 25%. The continuous time model pefoms well thoughout this entie egion (and beyond as shown in Fig. 7, even though the taditional theshold technique cannot even be applied to the data, as illustated in Fig. 2(b. A gid of 5 5 simulations wee un to poduce the coloumap in Fig. 8, and the posteio distibutions fo switching ates ae shown fo fou illustative cases. Note that the scales ae diffeent fo each contou plot. The lowe left plot shows a good convegence of the infeence to the tue value, whee the ates ae small. As both ates incease up the diagonal of the plane, the plots on the uppe left and uppe ight show the inceasing deviation of the infeence fom the tue value. In the top ight the switching ates ae so high that thee is an aveage of two switching events pe detecto inteval. The posteio distibution is tightly conveged since thee ae a lage numbe of switching events, but it is significantly diveged fom the tue value. The lowe ight plot consides a situation whee the on ate is vey high and off

8 Infeence eo (Euclidean dist FIG. 8. Accuacy of the single-step infeence to data simulated using a continuous model with vaious switching ates. The infeence eo taken as the Euclidean distance between the tue α, β and the mode of the posteio distibution. Fo small switching pobabilities the infeence gives good ageement with the known paametes, but fo highe switching pobabilites the infeence is less accuate. The cyan line indicates an appoximate theshold fo applicability of the single-step model of α β <., which coesponds appoximately to a elative eo 25%. ate is low. The convegence is not good fo this case, since the asymmety of switching ates esults in vey few switching off events the data emains in the on state fo most of the time. The single-inteval infeence povides a good estimate fo the slow ate, but damatically undeestimates the fast ate. IV. DTMC MULTI-STEP MODEL A. Recusion elation In the CTMC model of the pevious section the emitte may switch state at any time, and we deived a final integal ove the distibution of counts that has to be evaluated numeically. In geneal this may be difficult o inefficient, and in such cases the continuous model may be appoximated by a discete switching model that povides multiple possible switching points within each detecto inteval. Conside the n th detecto inteval, having bounday states s n and s n ; and obseved counts c n. We have aleady noted that this distibution is independent of the photon counts in adjacent detecto intevals, and so the infeence developed within this inteval will be equally valid fo all intevals in the measuement. Futhemoe, this holds within a detecto inteval also; the count distibution in a given subinteval is independent of adjacent subintevals given knowledge of the states at the boundaies of the sub-intevals. We intoduce the notation ζ x to denote the state of the emitte at faction x within the detecto inteval ( x, whee ζ = s n and ζ = s n. Also the counts in a factional inteval fom time u to time v will be denoted c u,v. In paticula, we intoduce the mid-point state ζ /2 and then maginalize ove it in ode to identify a useful ecuence elationship fo the photon count distibution. This intemediate state beaks the oiginal inteval into two subintevals, and clealy thee ae only c n + possible assignments fo the subintevals (c,, c 2 2, to give the total counts: {(c n,, (c n,,..., (, c n }. In the following discussion fequent distinction will be made between the fou sub-cases of the photon count distibution coesponding with the possible stat and finish state possibilities ζ u and ζ u+t of an inteval of length t. To make this easie we intoduce the following notation: f i,j (t, κ P (c u,u+t =κ, ζ u+t =j ζ u =i, Ω (5 Using the intoduced intemediate state ζ /2 gives ise to the expession m f i,j (t, κ = f i, (t/2, c d f,j (t/2, κ c d + c d = m f i, (t/2, c d f,j (t/2, κ c d (52 c d = which is the sum of the two discete convolutions of count pobability distibutions coesponding to the two possible mid-inteval states. This offes a fundamental ecusion elation between the pobability distibutions of counts ove successive halvings of the inteval. Fo a given inteval t, f i,j (t, is function of counts (a slice though a joint pobability distibution. Since in any given expeiment thee will be a maximum count in an inteval, the function f i,j (t, can be ecoded as a single dimensional aay. In light of this Eq. (52 shows that values of the functions f i,j (t, can be calculated numeically fom knowledge of the function values fo half the inteval, f i,j (t/2,.