METHOD FOR DETECTING A REVERSIBLY PHOTOSWITCHABLE CHEMICAL SPECIES IN A SAMPLE

Abstract

A method for detecting a reversibly photoswitchable chemical species in a sample, includes the steps of: a) illuminating the sample with light suitable to be absorbed by the chemical species triggering a reaction affecting an optical property of the chemical species, the first light being periodically-modulated at a fundamental modulation frequency; b) measuring the evolution of the optical property; c) extracting at least one of an in-phase component at a frequency which is an even multiple of the fundamental modulation frequency; and a quadrature component at a frequency which is an odd multiple of the fundamental modulation frequency of a signal representing the evolution; and d) using the extracted component or components for detecting the chemical species. An apparatus for carrying out the method is also provided.

Claims

1. A method for detecting a reversibly photoswitchable chemical species in a sample, comprising the steps of: a) illuminating the sample (S) with a first light at a first wavelength suitable to be absorbed by the chemical species triggering a reaction affecting at least one optical property of the chemical species, said first light being periodically-modulated at a fundamental modulation frequency; and b) measuring the evolution of the optical property of the chemical species; further comprising the steps of: c) extracting at least one of: an in-phase component at a frequency which is an even multiple, different from zero, of the fundamental modulation frequency; and a quadrature component at a frequency which is an odd multiple of the fundamental modulation frequency, greater than the fundamental modulation frequency itself, of a signal representing said evolution; and d) using the extracted component or components for detecting the chemical species.

2. The method of claim 1, wherein the average intensity of the light and the fundamental modulation frequency are chosen so as: either to maximize the amplitude of said or at least one said extracted component; or to minimize an interfering signal, at the frequency of said or at least one said extracted component, from a different chemical species in the sample.

3. The method of claim 1, wherein step a) further comprises illuminating the sample with a second light at a second wavelength, different from the first wavelength and suitable to be absorbed by the chemical species triggering said or another reaction affecting said optical property of the chemical species, said second light having a constant intensity.

4. The method of claim 1, wherein step a) further comprises illuminating the sample with a second light at a second wavelength, different from the first wavelength and suitable to be absorbed by the chemical species triggering said or a different reaction affecting said optical property of the chemical species, said second light being periodically modulated at the fundamental modulation frequency, in phase opposition with the first light.

5. The method of claim 3, wherein the ratio of the average intensities of the first and second light and the fundamental modulation frequency are chosen so as: either to maximize the amplitude of said or at least one said extracted component; or to minimize an interfering signal, at the frequency of said or at least one said extracted component, from a different chemical species in the sample.

6. The method of claim 1, wherein steps a), b) and c) are repeated a plurality of times for different illumination conditions, corresponding to different values of the light intensity or intensities and of the fundamental modulation frequency and wherein step d) comprises using a plurality of signal components corresponding to said different illumination conditions for discriminating between a plurality of chemical species.

7. The method of claim 6, wherein step d) comprises applying a machine learning or unmixing method to said plurality of signal components for estimating absolute or relative concentration of said chemical species.

8. The method of claim 1, any one of the preceding claims wherein light is modulated sinusoidally at the fundamental modulation frequency.

9. The method of claim 1, light is modulated by a sinusoid at the fundamental modulation frequency multiplied by a pulse train at a repetition frequency which is a multiple of, and at least ten times larger than, the fundamental modulation frequency.

10. The method of claim 9, wherein steps a) and b) are carried out by light scanning microscopy.

11. The method of claim 1, wherein said chemical species is fluorescent and said optical property whose evolution is measured is the intensity of fluorescence emission.

12. An apparatus for carrying out a method according to claim 1, comprising: at least a first controlled light source (LS1) configured for illuminating a sample with a first light at a first wavelength, said first light being periodically-modulated at a fundamental modulation frequency; a light detector (CAM) configured for measuring the evolution of an optical property of the sample; and a data processing device (DPD); wherein the data processing device is configured for extracting at least one of an in-phase component at a frequency which is an even multiple, different from zero, of the fundamental modulation frequency and a quadrature component at a frequency which is an odd multiple, greater than the fundamental modulation frequency itself, of the fundamental modulation frequency of a signal representing said evolution; and for using the extracted component or components for detecting a chemical species in the sample.

13. The apparatus of claim 12, further comprising a second controlled light source (LS2) configured for illuminating the sample with a second light at a second wavelength, different from the first wavelength, said second light having a constant intensity.

14. The apparatus of claim 12, further comprising a second controlled light source (LS2) configured for illuminating the sample with a second light at a second wavelength, different from the first wavelength, said second light being periodically modulated at the fundamental modulation frequency, in phase opposition with the first light.

15. The apparatus of claim 12, further comprising a scanning device for illuminating the sample.

Description

[0036] Additional features and advantages of the present invention will become apparent from the subsequent description, taken in conjunction with the accompanying drawings, wherein:

[0037] FIG. 1 illustrates the principle at the basis of the invention;

[0038] FIGS. 2A-2L show the theoretical dependence of concentrations (2A-2H) and fluorescence intensities (2I-2L) on dimensionless light intensity and modulation frequency;

[0039] FIGS. 3A-3F show the theoretical dependence of concentrations (3A-3C) and fluorescence intensity (3D-3F) of Dronpa-2 on control parameters (ratio of light intensities and modulation frequency);

[0040] FIGS. 4A-4F, 5A-5F, 6A-6F, and 7A-7I illustrate the discriminatory power of the inventive method;

[0041] FIG. 8 demonstrates the inherent optical sectioning power of the inventive method;

[0042] FIG. 9 is a schematic representation of an apparatus according to an embodiment of the invention;

[0043] FIG. 10 illustrates the light modulation scheme used in some embodiments of the invention; and

[0044] FIG. 11 is a schematic representation of an apparatus according to another embodiment of the invention.

[0045] On FIGS. 2A-2L, respectively associated to positive and negative values, the isodensity curves 95% (thick solid lines), 50% (medium solid lines), 10% (thin solid lines), 0% (dotted lines), −10% (thin dashed lines), −50% (medium dashed lines), −95% (thick dashed lines) refer to absolute extrema.

[0046] On FIGS. 3A-6F, respectively associated to positive and negative values, the isodensity curves 95% (solid lines), 50% (solid lines), 10% (solid lines), 0% (dotted lines), −10% (dashed lines), −50% (dashed lines), −95% (dashed lines) refer to absolute extrema.

[0047] The inventive method will now be described with reference to the case where the detected chemical species are reversibly photoswitchable fluorophores, however—as explained above—it also applies to non-fluorescent reversibly photoswitchable species, provided that they have different measurable properties (e.g. fluorescence emission, but also optical absorption, Raman signal, etc.) in their different states.

[0048] FIG. 1 illustrates the general principle at the basis of the invention. A sample S contains a reversibly photoswitchable species (“label”) P which, under the effect of light, can switch between a first state 1 and a second state 2 having different brightness (central part of the figure). In the example of FIG. 1 (see the left part of the figure), the sample is illuminated at a first wavelength λ.sub.1 with a sinusoidally-modulated light intensity I.sub.1 and, simultaneously, at a second wavelength λ.sub.2 with a constant light intensity I.sub.2 (alternatively, I.sub.2 could be zero, or be sinusoidally modulated at the same frequency as I.sub.1 but in phase opposition; moreover, light may be modulated periodically but not sinusoidally, e.g. by a square wave). In general, the intensity of the fluorescence emission varies periodically, but not sinusoidally; it can therefore be decomposed into harmonic components. More precisely, for each harmonic order there is an “in-phase” component In exhibiting the phase of the exciting light and an “out-of-phase” component Out exhibiting a 90° phase delay with respect to the exciting light. The amplitudes of the harmonic components may be computed e.g. by Fourier-transforming the detected intensity I.sub.F(t) of the fluorescence emission. The harmonics decomposition of the fluorescence emission is illustrated in the right part of FIG. 1.

[0049] A detailed theoretical description of the dynamic behavior of the reversibly photoswitchable fluorescent label P illuminated with light of intensity I(t) involving two components I.sub.1(t) and I.sub.2(t) at wavelengths λ.sub.1 and λ.sub.2 will be provided, under the assumption that this dynamic behavior is reliably described by the two-state exchange

[00001]$\begin{array}{cc}1\begin{array}{c}{k}_{12}\left(t\right)\\ \rightleftharpoons \\ k_{21}\left(t\right)\end{array}2& \left(1\right)\end{array}$

[0050] where the thermodynamically most stable state 1 is photochemically converted into the thermodynamically less stable state 2 at rate constant

k.sub.12(t)=σ.sub.12,1I.sub.1(t)+σ.sub.12,2I.sub.2(t)   (2)

[0051] from which it can relax back to the initial state 1 either by a photochemically- or a thermally-driven process at rate constant

k.sub.21(t)=σ.sub.21,1I.sub.1(t)+σ.sub.21,2I.sub.2(t0+k.sub.21 .sup.Δ  (3)

[0052] where σ.sub.12,1I.sub.1(t), σ.sub.12,2I.sub.2(t), σ.sub.21,1I.sub.1(t), σ.sub.21,2I.sub.2(t), and k.sub.21.sup.Δ respectively provide photochemical and thermal contributions to the rate constants. In that case, the molecular action cross-sections for photoisomerization σ.sub.12,1 and σ.sub.21,1 ( at λ.sub.10, σ.sub.12,2 and σ.sub.21,2 (at λ.sub.2), and the thermal rate constant k.sub.21 .sup.Δ fully define the dynamic behavior of the reversibly photoswitchable fluorescent label.

[0053] We consider that the system is submitted to a modulation of illumination obeying

I(t)=I.sub.1.sup.0[1+α.sub.1h.sub.1(ωt)]+I.sub.2.sup.0[1+α.sub.2h.sub.2(ωt)]  (4)

[0054] where α.sub.1 and h.sub.1(ωt) (respectively α.sub.2 and h.sub.2(ùt)) designate the amplitude of the light modulation at the wavelength λ.sub.1 (respectively λ.sub.2) and a periodic function with fundamental radial frequency ω respectively. We assume that the system is either uniformly illuminated or that it can be considered homogeneous at any time of its evolution.

[0055] Then we rely on the two-state exchange (1) to write the equation governing the concentration evolution

[00002]$\begin{array}{cc}\frac{d2}{\mathrm{dt}}=k_{12}\left(t\right)1-k_{21}\left(t\right)2& \left(5\right)\end{array}$$\mathrm{with}$$\begin{array}{cc}1+2=P_{tot}& \left(6\right)\end{array}$

where P.sub.tot is the total concentration. The steady value of the concentration 2 is

[00003]$\begin{array}{cc}{2}^0=\left[K_{12}^0/\left(1+K_{12}^0\right)\right]P_{\mathrm{tot}}& \left(7\right)\end{array}$$\begin{array}{cc}{K}_{12}^0=\frac{k_{12}^0}{k_{21}^0},& \left(8\right)\end{array}$

is the equilibrium constant, and k.sub.12.sup.0 and k.sub.21.sup.0

k.sub.12.sup.0=σ.sub.12,1I.sub.1.sup.0+σ.sub.12,2I.sub.2.sup.0   (9)

k.sub.21.sup.0=σ.sub.21,1I.sub.1.sup.0+σ.sub.21,2I.sub.2.sup.0+k.sub.21.sup.Δ,   (10)

are the rate constants associated with reaction (1) upon illuminating at light intensities I.sub.1.sup.0 and I.sub.2.sup.0:

For

[0056]
f=(2−2.sup.0)N/(P.sub.tot(α.sub.1+α.sub.2))   (11)

where N is the number of modulated light components of the illumination (N=1 or 2), the differential equation governing the temporal evolution of the concentrations of 1 and 2 (Eq. (5)) becomes

[00004]$\begin{array}{cc}\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{N}{1+\frac{{\alpha }_2}{{\alpha }_1}}\frac{K_{12}^0}{{\left(1+K_{12}^0\right)}^2}{\Delta }_{12,1}^0-{\alpha }_1\left({\sigma }_{12,1}+{\sigma }_{21,1}\right)I_1^0{\tau }_{12}^{0}f\left(\theta x\right)\right]h_1\left(\theta x\right)+\left[\frac{N\frac{{\alpha }_2}{{\alpha }_1}}{1+\frac{{\alpha }_2}{{\alpha }_1}}\frac{K_{12}^0}{{\left(1+K_{12}^0\right)}^2}{\Delta }_{12,2}^0-{\alpha }_2\left({\sigma }_{12,2}+{\sigma }_{21,2}\right)I_2^0{\tau }_{12}^{0}f\left(\theta x\right)\right]h_2\left(\theta x\right)& \left(12\right)\end{array}$$\begin{array}{cc}\mathrm{where}:x=\frac{t}{{\tau }_{12}^0}& \left(13\right)\end{array}$$\begin{array}{cc}\theta =\omega {\tau }_{12}^0,& \left(14\right)\end{array}$

[0057] with τ.sub.12.sup.0 the relaxation time of the reaction (1):

[00005]$\begin{array}{cc}{\tau }_{12}^0=\frac{1}{k_{12}^0+k_{21}^0}& \left(15\right)\end{array}$

[0058] and Δ.sub.12,1.sup.0 (respectively Δ.sub.12,2.sup.0) the difference of the relative contributions of the average of the modulated light at wavelength λ.sub.1 (I.sub.1.sup.0) (respectively λ.sub.2 (I.sub.2.sup.0)) to drive the transition from 1 to 2 (respectively from 2 to 1):

[00006]$\begin{array}{cc}{\Delta }_{12,1}^0=\frac{{\sigma }_{12,1}{I}_1^0}{{\sigma }_{12,1}{I}_1^0+{\sigma }_{12,2}{I}_2^0}-\frac{{\sigma }_{21,1}{I}_1^0}{{\sigma }_{21,1}{I}_1^0+{\sigma }_{21,2}{I}_2^0+k_{21}^{\Delta }}& \left(16\right)\end{array}$$\begin{array}{cc}{\Delta }_{12,2}^0=\frac{{\sigma }_{12,2}{I}_2^0}{{\sigma }_{12,1}{I}_1^0+{\sigma }_{12,2}{I}_2^0}-\frac{{\sigma }_{21,2}{I}_2^0}{{\sigma }_{21,1}{I}_1^0+{\sigma }_{21,2}{I}_2^0+k_{21}^{\Delta }}.& \left(17\right)\end{array}$

[0059] This equation is subsequently solved in three distinct cases: [0060] 1. Case 1. We first consider the case of reversibly photoswitchable species whose forward and backward conversions are both photochemically driven using two different wavelengths with only one of them (λ.sub.1) being modulated (N=1, α.sub.2=0). In the following, this case is addressed by making the following assumptions validated in most reported systems: [0061] a. The transition from the state 1 to the state 2 is assumed to be controlled by the illumination at the wavelength λ.sub.1 (σ.sub.12,2I.sub.2.sup.0<<σ.sub.12,1I.sub.1.sup.0) whereas the transition from the state 2 to the state 1 is driven by the illumination at the wavelength λ.sub.2 (σ.sub.21,1I.sub.1.sup.0+k.sub.21.sup.Δ<<σ.sub.21,2I.sub.2.sup.0), which means that we consider to be in a regime of high enough light intensity I.sub.2.sup.0. [0062] b. The illumination at the wavelength λ.sub.1 contributes more to the rate constants from the state 1 to the state 2 than from the state 2 to the state 1 (σ.sub.21,1<<σ.sub.12,1). Thus Eqs. (8), (15), and (16) become:

[00007]$\begin{array}{cc}{K}_{12}^{0,a}=\frac{{\sigma }_{12,1}{I}_1^0}{{\sigma }_{21,2}{I}_2^0}& \left(18\right)\end{array}$$\begin{array}{cc}{\tau }_{12}^{0,a}=\frac{1}{{\sigma }_{12,1}{I}_1^0+{\sigma }_{21,2}{I}_2^0}& \left(19\right)\end{array}$$\begin{array}{cc}{\Delta }_{12,1}^{0,a}=1& \left(20\right)\end{array}$

[0063] Under the latter assumptions, Eq. (12) yields:

[00008]$\begin{array}{cc}\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{K_{12}^{0,a}}{{\left(1+K_{12}^{0,a}\right)}^2}-{\alpha }_1\frac{K_{12}^{0,a}}{1+K_{12}^{0,a}}f\left(\theta x\right)\right]h_1\left(\theta x\right)& \left(21\right)\end{array}$ [0064] 2. Case 2. We subsequently consider the case of reversibly photoswitchable species whose forward conversion is photochemically driven, whereas the backward reaction involves thermal control. We correspondingly adopt I.sub.2.sup.0=0 (N=1, α.sub.2=0) and rely on reported results to further assume that: [0065] a. The illumination at the wavelength λ.sub.1 contributes more to the rate constants from the state 1 to the state 2 than from the state 2 to the state 1 (σ.sub.21,1<<σ.sub.12,1). [0066] b. The transition from the state 2 to the state 1 is driven by the thermal contribution (k.sub.21.sup.Δ>>σ.sub.21,1I.sub.1.sup.0). Thus Eqs. (8), (15), and (16) become:

[00009]$\begin{array}{cc}{K}_{12}^{0,b}=\frac{{\sigma }_{12,1}{I}_1^0}{k_{21}^{\Delta }}& \left(22\right)\end{array}$$\begin{array}{cc}{\tau }_{12}^{0,b}=\frac{1}{{\sigma }_{12,1}{I}_1^0+k_{21}^{\Delta }}& \left(23\right)\end{array}$$\begin{array}{cc}{\Delta }_{12,1}^{0,b}=1& \left(24\right)\end{array}$

[0067] Under the latter assumptions, Eq. (12) yields:

[00010]$\begin{array}{cc}\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{K_{12}^{0,b}}{{\left(1+K_{12}^{0,b}\right)}^2}-{\alpha }_1\frac{K_{12}^{0,b}}{1+K_{12}^{0,b}}f\left(\theta x\right)\right]h_1\left(\theta x\right)& \left(25\right)\end{array}$ [0068] 3. Case 3. We eventually consider the case of reversibly photoswitchable species whose forward and backward conversions are both photochemically driven using two different wavelengths which are both modulated (N=2). In addition to the assumptions (a) and (b) of the case 1, we further adopt the following assumptions to address this case: [0069] a. The illumination at the wavelength λ.sub.2 contributes more to the rate constants from the state 2 to the state 1 than from the state 1 to the state 2 (σ.sub.12,1<<σ.sub.21,2). Thus Eqs (8), (15), (16) and (17) become:

[00011]$\begin{array}{cc}{K}_{12}^{0,c}=\frac{{\sigma }_{12,1}{I}_1^0}{{\sigma }_{21,2}{I}_2^0}& \left(26\right)\end{array}$$\begin{array}{cc}{\tau }_{12}^{0,c}=\frac{1}{{\sigma }_{12,1}{I}_1^0+{\sigma }_{21,2}{I}_2^0}& \left(27\right)\end{array}$$\begin{array}{cc}{\Delta }_{12,1}^{0,c}=1& \left(28\right)\end{array}$$\begin{array}{cc}{\Delta }_{12,2}^{0,c}=-1& \left(29\right)\end{array}$ [0070] b. The periodic functions h.sub.1(θx) and h.sub.2(θx) are opposite (h.sub.2(θx)=−h.sub.1(θx)). [0071] Under the latter assumptions, Eq. (12) yields:

[00012]$\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{2K_{12}^{0,c}}{{\left(1+K_{12}^{0,c}\right)}^2}-\frac{{\alpha }_1{K}_{12}^{0,c}-{\alpha }_2}{1+K_{12}^{0,c}}f\left(\theta x\right)\right]h_1\left(\theta x\right)$

[0072] The final equations (21) and (25) are similar. In the following, we correspondingly adopt the following general equations for the first two cases:

[00013]$\begin{array}{cc}\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{K_{12}^0}{{\left(1+K_{12}^0\right)}^2}-\alpha \frac{K_{12}^0}{1+K_{12}^0}f\left(\theta x\right)\right]h\left(\theta x\right)& \left(30\right)\end{array}$

and for the third case;

[00014]$\begin{array}{cc}\frac{\mathrm{df}\left(\theta x\right)}{\mathrm{dx}}=-f\left(\theta x\right)+\left[\frac{2K_{12}^0}{{\left(1+K_{12}^0\right)}^2}-\frac{{\alpha }_1{K}_{12}^0-{\alpha }_2}{1+K_{12}^0}f\left(\theta x\right)\right]h\left(\theta x\right)& \left(31\right)\end{array}$

[0073] Beyond the relaxation time τ.sub.12.sup.0, one enters into a permanent regime in which f(θx) is a continuous periodic function. Then f(θx) can be expressed as a Fourier series

[00015]$\begin{array}{cc}f\left(\theta x\right)=a_0+\underset{n=1}{\overset{+\infty }{.Math.}}\left[a_n\cos \left(n\theta x\right)+b_n\sin \left(n\theta x\right)\right]& \left(32\right)\end{array}$

[0074] where a.sub.n and b.sub.n designate the amplitudes of the n-th components of the Fourier series. The a.sub.n and b.sub.n terms can be extracted from Eqs. (31), (32) upon identifying the amplitudes of the components of the same order (harmonic balance).

[0075] Then we obtain the following expressions of the concentrations in 1 and 2:

[00016]$\begin{array}{cc}2=2^0+\frac{{\alpha }_1+{\alpha }_2}{N}\left\{a_0+\underset{n=1}{\overset{+\infty }{.Math.}}\left[a_n\cos \left(n\theta x\right)+b_n\sin \left(n\theta x\right)\right]\right\}{P}_{\mathrm{tot}}& \left(33\right)\end{array}$$\begin{array}{cc}1=1^0-\frac{{\alpha }_1+{\alpha }_2}{N}\left\{a_0+\underset{n=1}{\overset{+\infty }{.Math.}}\left[a_n\cos \left(n\theta x\right)+b_n\sin \left(n\theta x\right)\right]\right\}{P}_{\mathrm{tot}}& \left(34\right)\end{array}$

The following explicit notations are introduced:

=b.sub.n   (35)

=a.sub.n   (36)

where in means in phase and out in quadrature with light modulation and norm stands for scaling of the amplitudes. and respectively represent the in-phase and the out-of-phase parts of the n-order component of the normalized (with respect to P.sub.tot) concentration of the 2 state. Similarly, it is possible to define =− and =−

[0076] Many observables (e.g. absorbance, electrophoretic mobility, . . . ) from the photoswitchable probe P vary linearly with the concentrations in 1 and 2. Therefore the harmonic components of the time dependency of these observables are expressed (apart from a proportionality factor) by the equations above for the harmonic components of 1 and 2.

[0077] Fluorescence emission, which is a hugely important observable for practical application of the present invention e.g. to biology, has a more complex time dependence:

I.sub.F(t)=[Q.sub.1,11(t)].sub.1(t)+[Q.sub.1,21(t)+Q.sub.2,22(t)]I.sub.2(t)   (37)

[0078] where Q.sub.i,j is the brightness of state I (i=1, 2) when illuminated at wavelength λ.sub.j (j=1, 2).

[0079] Like concentrations, intensity of fluorescence emission may be decomposed into in-phase and out-of-phase harmonic components:

I.sub.F(t)=+Σ.sub.n=1.sup.∞[.sup. sin(nθx)+ cos (nθx)]  (38)

[0080] As the inventive method uses harmonic components of the fluorescence intensity (or of other observables) to detect and/or discriminate reversibly photoswitchable species, in preferred embodiments of the invention the modulation of the light intensity is sinusoidal, so that no harmonic component is introduced by the linear part of the dynamical response of the species. Hereafter, then, the particular cases of sinusoidal modulation at one or two wavelengths will be considered in detail. The invention, however, is not limited to these preferred cases, and non-sinusoidal light modulation may also be used. A particularly interesting case is that of a square-wave (i.e. on-off) light intensity modulation, which is particularly easy to implement. As the Fourier development of a square wave only comprises odd harmonics, when an on-off light intensity modulation is applied, it is advantageous to perform the detection/identification of the reversibly photoswitchable species by only taking into account even harmonics of the relevant observable(s).

[0081] Light intensity sinusoidally modulated at one wavelength can be expressed as:

I(t)=I.sub.1.sup.0[1+α sin(ωt)]+I.sub.2.sup.0   (39)

[0082] Injecting (38) into (30) and then extracting the harmonic components according to (32) allows the numerical calculation of the a.sub.n and b.sub.n terms (amplitudes of the out-of-phase and in-phase harmonic components, n being the harmonic order).

[0083] Truncation at the fifth order (n=5) has been found to be sufficient to yield convergence: the dependencies of a.sub.0, a.sub.1, b.sub.1, a.sub.2 and b.sub.2 do not significantly evolve beyond the third order (n=3), whereas the convergences of a.sub.3, b.sub.3 and a.sub.4, b.sub.4 are observed beyond the fourth (n=4) and fifth order (n=5) respectively.

[0084] The amplitudes a.sub.i and b.sub.i do not all share the same sign : a.sub.0, a.sub.1, b.sub.3, a.sub.4 and b.sub.4 are negative whereas b.sub.1, a.sub.2, b.sub.2 and a.sub.3 are positive. The odd order out-of-phase (a.sub.2k+1, with k∈) and even order in-phase (b.sub.2k, with k∈*) amplitudes exhibit an optimum at coordinates (k.sub.12.sup.0,R,θ.sup.R) in the (K.sub.12.sup.0,θ) space, whose position shifts respectively toward higher and smaller K.sub.12.sup.0 and θ values when k increases; otherwise said, they show a resonant behavior (note that the resonant behavior of the a.sub.1 term had already been identified, and is exploited by the OPIOM technique mentioned above).

[0085] Moreover, the bandwidth of the resonance decreases, at a larger rate along θ. As it will be discussed further, this allows improving the resolving power of the inventive method compared with the OPIOM technique.

[0086] In contrast, the other amplitudes (a.sub.2k, with k∈ and b.sub.2k+1, with k∈) exhibit a ridge shape; they can be optimized with respect to K.sub.12.sup.0 but only display a threshold behavior with respect to θ. These components are discarded in preferred embodiments of the inventive method.

[0087] Another interesting result is that the intensity of the harmonic components decreases, but rather slowly—roughly by a factor of 2 whenever the order n increases by 1. This means that, in most cases of practical interest, at least the first few harmonics (up to the third or the fourth) will be usable.

[0088] The inventors have found the following approximate analytical expressions for the amplitudes of the in-phase even harmonics and of the out-of-phase odd harmonics:

[0089] For the out-of-phase components at odd orders: n=2k+1; k∈:

[00017]$\begin{array}{cc}{a}_{2k+1}={\left(-1\right)}^{k+1}{A}_{2k+1}{\alpha }^{2k}\frac{{\left(K_{12}^0\right)}^{2k+1}}{{\left(1+K_{12}^0\right)}^{2k+2}}\frac{\theta P^k\left[{\theta }^2\right]}{{.Math.}_{i=1}^{2k+1}\left[1+i^2{\theta }^2\right]}& \left(40\right)\end{array}$

[0090] For the in-phase components at evenorders: n=2k; k∈*:

[00018]$\begin{array}{cc}{b}_{2k}={\left(-1\right)}^{k+1}{A}_{2k}{\alpha }^{2k-1}\frac{{\left(K_{12}^0\right)}^{2k}}{{\left(1+K_{12}^0\right)}^{2k+1}}\frac{\theta P^{k-1}\left[{\theta }^2\right]}{{.Math.}_{i=1}^{2k}\left[1+i^2{\theta }^2\right]}& \left(41\right)\end{array}$

[0091] where A.sub.n is a positive amplitude factor and P.sup.n a polynomial of order n, with P.sup.n[θ.sup.2]˜1 for θ≲1.

[0092] If properties proportional to the concentrations 1 and 2 are used for detecting and/or discriminating photoswitchable species, it is useful to work with normalized observables:

[00019]$\begin{array}{cc}{O}_{j,\mathrm{norm}}\left(t\right)=\frac{Q_{1,j}1\left(t\right)+Q_{2,j}2\left(t\right)}{\left(Q_{2,j}-Q_{1,j}\right)\alpha P_{\mathrm{tot}}}& \left(42\right)\end{array}$

[0093] where Q.sub.1,j and Q.sub.2,j designate the specific normalized observables associated to the observation at the wavelength λ.sub.j for the states 1 and 2 respectively. Harmonics decomposition gives:

O.sub.j,norm(t)=+Σ.sub.n=1.sup.∞[ sin(nθx)+ cos(nθx)]  (43)

=b.sub.n   (44)

=a.sub.n   (45)

[0094] Similarly, when the observable used for detecting/identifying reversibly photoswitchable species is fluorescence emission, it is useful to define a normalized fluorescence intensity:

[00020]$\begin{array}{cc}{I}_{F,\mathrm{norm}}\left(t\right)=\frac{\left[Q_{1,1}1\left(t\right)+Q_{2,1}2\left(t\right)\right]I_1\left(t\right)+\left[Q_{1,2}1\left(t\right)+Q_{2,2}2\left(t\right)\right]I_2\left(t\right)}{\left[\left(Q_{2,1}-Q_{1,1}\right)I_1^0+\left(Q_{2,2}-Q_{1,2}\right)I_2^0\right]\alpha P_{\mathrm{tot}}}.& \left(46\right)\end{array}$

[0095] Harmonics decomposition gives:

I.sub.F,norm(t)=+Σ.sub.n=1.sup.∞[ sin(nθx)+ cos(nθx)].   (47)

[0096] with:

[00021]$=\begin{array}{cc}\frac{K_{12}^0\left[1+Q_{{\lambda }_1}{K}_{12}^0\right]}{\left[K_{12}^0\left(Q_{{\lambda }_1}-1\right)+\mathrm{\Pi \varsigma }\left(Q_{{\lambda }_2}-1\right)\right]\left(1+K_{12}^0\right)}+b_1+\frac{\alpha }{2}\frac{1}{1+\frac{\xi }{K_{12}^0}}\left(2a_0-a_2\right)& \left(48\right)\end{array}$$\begin{array}{cc}=a_1+\frac{\alpha }{2}\frac{1}{1+\frac{\xi }{K_{12}^0}}{b}_2& \left(49\right)\end{array}$$\begin{array}{cc}=b_n+\frac{\alpha }{2}\frac{1}{1+\frac{\xi }{K_{12}^0}}\left(a_{n-1}-a_{n+1}\right)& \left(50\right)\end{array}$$\begin{array}{cc}=a_n-\frac{\alpha }{2}\frac{1}{1+\frac{\xi }{K_{12}^0}}\left(b_{n-1}-b_{n+1}\right)& \left(51\right)\end{array}$$\mathrm{for}n>1,\mathrm{with}\varsigma =\frac{{\sigma }_{121}}{{\sigma }_{212}},\Pi =\frac{Q_{1,2}}{Q_{1,1}},\xi =\mathrm{\Pi \varsigma }\frac{Q_{{\lambda }_2}-1}{Q{\lambda }_1-1}\mathrm{and}{Q}_{{\lambda }_i}=\frac{Q_{2,i}}{Q_{1,i}}.$

[0097] Light intensity sinusoidally modulated at two wavelengths with phase opposition can be expressed as:

I(t)=.sub.1.sup.0[1+α.sub.1 sin(ωt)]+I.sub.2.sup.0[1−α.sub.2 sin(ωt)]  (52)

[0098] Injecting (52) into (31) and then extracting the harmonic components according to (32) allows the numerical calculation of the a.sub.n and b.sub.n terms (amplitudes of the out-of-phase and in-phase harmonic components, n being the harmonic order).

[0099] Truncation at the fifth order (n=5) has been found to be sufficient to yield convergence: the dependencies of a.sub.0, a.sub.1, b.sub.1, a.sub.2 and b.sub.2 do not significantly evolve beyond the third order (n=3), whereas the convergences of a.sub.3, b.sub.3 and a.sub.4, b.sub.4 are observed beyond the fourth (n=4) and fifth order (n=5) respectively.

[0100] The amplitudes a.sub.i and b.sub.i do not all share the same sign: a.sub.1 and b.sub.3 are negative, b.sub.1 and a.sub.3 are positive whereas a.sub.0, a.sub.2, b.sub.2, a.sub.4 and b.sub.4 exhibit positive and negative components. The odd order amplitudes (a.sub.2k+1 and b.sub.2k+k1, with k∈) exhibit a symmetry axis in K.sub.12.sup.0=1 and the even order amplitudes (a.sub.2k b.sub.2k, with k∈) exhibit an antisymmetry axis in K.sub.12.sup.0=1. The odd order out-of-phase (a.sub.2k+1, with k∈) and even order in-phase (b.sub.2k, with k∈*) amplitudes exhibit at least an optimum in the (K.sub.12.sup.0,θ) space. More specifically, a.sub.1 exhibits an optimum whereas the superior orders (a.sub.2k+1 and b.sub.2k, k>0) exhibit two optima, whose positions shift respectively toward higher and smaller K.sub.12.sup.0 (symmetry in K.sub.12.sup.0=1) and smaller 0 values when k increases; otherwise said, they show a resonant behavior (note that the resonant behavior of the a.sub.1 term had already been identified, and is exploited by the Speed OPIOM technique mentioned above).

[0101] Moreover, the bandwidth of the resonance decreases, at a larger rate along θ. As it will be discussed further, this allows improving the resolving power of the inventive method compared with the OPIOM and Speed OPIOM techniques.

[0102] In contrast, the other amplitudes (a.sub.2k, with k∈ and b.sub.2k+1, with k∈) exhibit at least a ridge shape; b.sub.1 exhibits a ridge shape whereas the superior orders (b.sub.2k+1 and a.sub.2k, k>0) exhibit two ridge shapes which can be optimized with respect to K.sub.12.sup.0 but only exhibit a threshold behavior with respect to θ. These components are discarded in preferred embodiments of the inventive method.

[0103] With respect to the observations performed in the case of the single light modulation and by restricting analysis to the terms exhibiting at least a resonance: [0104] The evolution of the position and bandwidth of the resonance along the θ axis is similar; [0105] The position and bandwidth of the resonance along the K.sub.12.sup.0 axis shift faster when increasing the order n of the concentration terms; [0106] The decrease of the harmonics intensity with increasing harmonic order is similar; [0107] The optimal amplitudes in the (K.sub.12.sup.0,θ) space are the same except for a.sub.1, which is twice higher. [0108] The inventors have found the following approximate analytical expressions for the amplitudes of the in-phase even harmonics and of the out-of-phase odd harmonics:

[0109] For the out-of-phase components at odd orders: n=2k+1; k∈:

[00022]$\begin{array}{cc}{a}_{2k+1}={\left(-1\right)}^{k+1}{A}_{2k+1}\frac{{K_{12}^0\left({\alpha }_1{K}_{12}^0-{\alpha }_2\right)}^{2k}}{{\left(1+K_{12}^0\right)}^{2k+2}}\frac{\theta P^k\left[{\theta }^2\right]}{{\Pi }_{i=1}^{2k+1}\left[1+i^2{\theta }^2\right]}& \left[53\right]\end{array}$

[0110] For the in-phase componeents at even orders: n=2k; k∈*:

[00023]$\begin{array}{cc}{b}_{2k}={\left(-1\right)}^{k+1}{A}_{2k}\frac{{K_{12}^0\left({\alpha }_1{K}_{12}^0-{\alpha }_2\right)}^{2k-1}}{{\left(1+K_{12}^0\right)}^{2k+1}}\frac{\theta P^{k-1}\left[{\theta }^2\right]}{{.Math.}_{i=1}^{2k}\left[1+i^2{\theta }^2\right]}& \left(54\right)\end{array}$

[0111] where A.sub.n is a positive amplitude factor and P.sup.n a polynomial of order n, with P.sup.n[θ.sup.2]−1 for θ≲1.

[0112] If properties proportional to the concentrations 1 and 2 are used for detecting and/or discriminating photoswitchable species, it is useful to work with normalized observables:

[00024]$\begin{array}{cc}{O}_{j,\mathrm{norm}}\left(t\right)=\frac{Q_{1,j}1\left(t\right)+Q_{2,j}2\left(t\right)}{\left(Q_{2,j}-Q_{1,j}\right)\frac{{\alpha }_1+{\alpha }_2}{2}{P}_{\mathrm{tot}}}& \left(55\right)\end{array}$

[0113] where Q.sub.i,j and Q.sub.2,j designate the specific normalized observables associated to the observation at the wavelength λ.sub.j for the states 1 and 2 respectively. Harmonic decomposition gives:

O.sub.j,norm(t)=+Σ.sub.n=1.sup.∞[ sin(nθx)+ cos(nθx)]  (56)

=a.sub.n   (57)

=b.sub.n   (58)

[0114] Similarly, when the observable used for detecting/identifying reversibly photoswitchable species is fluorescence emission, it is useful to define a normalized fluorescence intensity:

[00025]$\begin{array}{cc}{I}_{F,norm}\left(t\right)=\frac{\left[Q_{1,1}1\left(t\right)+Q_{2,1}2\left(t\right)\right]I_1\left(t\right)+\left[Q_{1,2}1\left(t\right)+Q_{2,2}2\left(t\right)\right]I_2\left(t\right)}{\left[\left(Q_{2,1}-Q_{1,1}\right)I_1^0+\left(Q_{2,2}-Q_{1,2}\right)I_2^0\right]\frac{{\alpha }_1+{\alpha }_2}{2}{P}_{tot}}& \left(59\right)\end{array}$

[0115] Harmonic dcomposition gives:

[00026]$\begin{array}{cc}{I}_{F,\mathrm{norm}}\left(t\right)=+{.Math.}_{n=1}^{\infty }\left[\sin \left(n\theta x\right)+\cos \left(n\theta x\right)\right]& \left(60\right)\end{array}$

[0116] with:

[00027]$\begin{array}{cc}=\frac{2\left[{\alpha }_1{K}_{12}^0\left(1+Q_{{\lambda }_1}{K}_{12}^0\right)-{\alpha }_2\mathrm{\Pi \varsigma }\left(1+Q_{{\lambda }_2}{K}_{12}^0\right)\right]}{\left[K_{12}^0\left(Q_{{\lambda }_1}-1\right)+\mathrm{\Pi \varsigma }\left(Q_{{\lambda }_2}-1\right)\right]\left({\alpha }_1+{\alpha }_2\right)\left(1+K_{12}^0\right)}+b_1+\frac{1}{2}\frac{{\alpha }_1-{\alpha }_2\frac{\xi }{K_{12}^0}}{1+\frac{\xi }{K_{12}^0}}\left(2a_0-a_2\right)& \left(61\right)\end{array}$$\begin{array}{cc}=a_1+\frac{1}{2}\frac{{\alpha }_1-{\alpha }_2\frac{\xi }{K_{12}^0}}{1+\frac{\xi }{K_{12}^0}}{b}_2& \left(62\right)\end{array}$$\begin{array}{cc}=b_n+\frac{1}{2}\frac{{\alpha }_1-{\alpha }_2\frac{\xi }{K_{12}^0}}{1+\frac{\xi }{K_{12}^0}}\left(a_{n-1}-a_{n+1}\right)& \left(63\right)\end{array}$$\begin{array}{cc}=a_n-\frac{1}{2}\frac{{\alpha }_1-{\alpha }_2\frac{\xi }{K_{12}^0}}{1+\frac{\xi }{K_{12}^0}}\left(b_{n-1}-b_{n+1}\right)& \left(64\right)\end{array}$$\mathrm{for}n>1,\mathrm{with}\varsigma =\frac{{\sigma }_{121}}{{\sigma }_{212}},\Pi =\frac{Q_{1,2}}{Q_{1,1}},\xi =\mathrm{\Pi \varsigma }\frac{Q_{{\lambda }_2}-1}{Q{\lambda }_1-1}\mathrm{and}{Q}_{{\lambda }_i}=\frac{Q_{2,i}}{Q_{1,i}}.$

[0117] The equations above allow computing maps of the different harmonic components of the fluorescence intensity and of the concentration of state 2 as a function of the control parameters: I.sub.1, I.sub.2, ω, or suitable functions thereof. FIGS. 2A-2L show such maps in the (K.sub.12.sup.0,θ) space. It is recalled that K.sub.12.sup.0 is a dimensionless parameter proportional to I.sub.1.sup.0 if I.sub.2.sup.0=0 or I.sub.1.sup.0/I.sub.2.sup.0 otherwise and θ is a dimensionless parameter proportional to ω. Interestingly, when the harmonic components are expressed as a function of K.sub.12.sup.0 and θ, they do not explicitly depend on the photochemical properties of the reversibly photoswitchable species (this is possible because K.sub.12.sup.0 and θ do depend on them). As a consequence, the maps of FIGS. 2A-2L are “universal”: they depend on the illumination conditions (single wavelength, double wavelength with single modulation, double wavelength with double modulation in phase opposition), but not on a specific reversibly photoswitchable species. In all cases, dual-wavelength illumination with either single modulation or dual modulation in phase opposition, the modulation is sinusoidal and with large amplitude (α=1). Single-wavelength illumination is treated as a specific case of dual-wavelength illumination with single modulation, wherein the intensity of the constant component is zero. More specifically: [0118] FIG. 2A is a map of for single-wavelength modulation; [0119] FIG. 2B is a map of .sub.norm.sup.2,in for single-wavelength modulation; [0120] FIG. 2C is a map of for single-wavelength modulation; [0121] FIG. 2D is a map of for single-wavelength modulation; [0122] FIG. 2E is a map of for dual-wavelength modulation; [0123] FIG. 2F is a map of for dual-wavelength modulation; [0124] FIG. 2G is a map of for dual-wavelength modulation; [0125] FIG. 2H is a map of for dual-wavelength modulation; [0126] FIG. 2I is a map of for single-wavelength modulation; [0127] FIG. 2J is a map of for single-wavelength modulation; [0128] FIG. 2K is a map of for single-wavelength modulation; [0129] FIG. 2L is a map of fo single wavelength modulation.

[0130] In the case of single wavelength modulation involving two wavelengths, the fluorescence maps displayed in FIGS. 2I-2L correspond to the case where the modulated wavelength switches the species from a bright to a dark state associated with a significant brightness at the modulated wavelength. When the modulated wavelength switches the species from a dark to a bright state, the fluorescence maps differ only by the sign from the maps displayed in FIGS. 2I-2L. For a bright state with two different brightnesses at the two wavelengths, qualitatively different fluorescence maps are obtained for light modulation at the wavelength associated with the smallest brightness. In this case, the fluorescence maps are close to the concentration maps.

[0131] It can be seen that the first order amplitude is well-approximated by the out-of-phase response obtained with both light modulations of small amplitude.([Quérard 2015]; [Quérard 2017]) In both single and dual wavelength light modulation, it exhibits a similar symmetrical peak with respect to K.sub.12and θ with a broader bandwidth along K.sub.12.sup.0 than along θ. Its resonance at K.sub.12.sup.0=1 and θ=1 optimizes both the composition shift occurring with light modulation (maximized when the forward and backward reactions driven by illumination occur at the same rate so that K.sub.12.sup.0=1) and the phase lag of the response to light modulation (π/2 when the radial frequency of the light modulation ω is matched with the exchange relaxation time τ.sub.12.sup.0 so that θ=1). Interestingly at resonance is twice higher with dual wavelength than with single wavelength light modulation. In contrast to , single and dual wavelength modulation generate different maps of the higher terms of the concentration modulation. With single wavelength modulation, the position of the resonance shifts towards higher and lower values for K.sub.11.sup.0 and θ respectively from to . Along the K.sub.12.sup.0-axis, approximate analytical expressions yield resonance for K.sub.12.sup.0,R=n where n designates the harmonics order. In contrast, a simple analytical expression of the resonance along the θ-axis could not be derived. Numerical computation showed that the resonant θ.sup.R value dropped from 1 as n increases. With dual wavelength modulation, not anymore one but two resonances are observed at higher orders n in the (K.sub.12.sup.0,θ) space with either a symmetry or anti-symmetry relationship with respect to the K.sub.12.sup.0=1 axis. The position of these resonances shifts symmetrically towards higher and lower values for K.sub.12.sup.0 (with |K.sub.12.sup.0,R|=n±√{square root over (n.sup.2−1)}), and towards lower values for θ from to . For both modulated illuminations, the resonance peak becomes more and more asymmetric and its bandwidth is reduced along the K.sub.11.sup.0 and θ axes when going from to .

[0132] As per fluorescence, essentially exhibits the resonance of but with a lower amplitude and with a slight change of shape. In contrast to and exhibit either one (if λ.sub.2=405 nm is modulated) or two (if λ.sub.1=480 nm is modulated) resonances in the (K.sub.12.sup.0,θ) space. Interestingly if λ.sub.1=480 nm is modulated, and possess one line of vanishing amplitude, and possesses two lines of vanishing amplitude. The highest peak experiences a more pronounced drop of the bandwidth upon increasing n than the resonance observed for the corresponding concentration terms. In particular, the drop is more pronounced along the θ.sub.axis than the K.sub.12.sup.0 axis.

[0133] The dimensionless (K.sub.12.sup.0,θ) space is appropriate for theoretical computations and analysis since all reversibly photoswitchable species exhibit the same behavior in this space. However the illumination control parameters (ω,I.sub.1.sup.9,I.sub.2.sup.0) are more relevant for applications in analysis or imaging. Indeed, in such a space, the positions of the resonance peaks are species-dependent, and therefore allow discriminating different reversibly photoswitchable species. In the case of single-wavelength illumination, a suitable parameter space will be (ω,I.sub.1.sup.0), while in the case of double-wavelength illumination (ω/I.sub.1.sup.0,I.sub.2.sup.0/I.sub.1.sup.0) will be preferably used.

[0134] This will be illustrated with the help of FIGS. 3A-3F, 4A-4F, 5A-5F and 6A-6F.

[0135] FIGS. 3A-3F show the maps of the oscillating concentration and fluorescence intensity components at different orders n=1, 2, 3, and 4 for the reversibly photoswitchable fluorescent protein Dronpa-2 submitted to dual-wavelength illumination in the (ω/I.sub.1.sup.0,I.sub.2.sup.0/I.sub.1.sup.0) space (a single-wavelength illumination which can be considered as a particular case of dual-wavelength illumination with I.sub.2.sup.0=0 in the (K.sub.12.sup.0,θ) space, is not accounted for by these maps). Dronpa-2 is switched to a dark state by blue light (480 nm) and returns to the initial bright state by violet light (405 nm), a behavior known as negative photochromism.

[0136] More precisely: [0137] FIG. 3A is a map of the oscillating concentration components (“out” for n=1, 3: “in” for n=2, 4) for single-wavelength modulation with λ.sub.1=480 nm (modulated) and λ.sub.2=405 nm (constant intensity); [0138] FIG. 3B is a map of the oscillating concentration components (“out” for n=1, 3: “in” for n=2, 4) for single-wavelength modulation with λ.sub.2=405 nm (modulated) and λ.sub.1=480 nm (constant intensity); [0139] FIG. 3C is a map of the oscillating concentration components (“out” for n=1, 3: “in” for n=2, 4) for dual-wavelength modulation with λ.sub.1=480 nm and λ.sub.2=405 nm; [0140] FIG. 3D is a map of the oscillating fluorescence intensity components (“out” for n=1, 3: “in” for n=2) for single-wavelength modulation with λ.sub.1=480 nm (modulated) and λ.sub.2=405 nm (constant intensity); [0141] FIG. 3E is a map of the oscillating fluorescence intensity components (“out” for n=1, 3: “in” for n=2) for single-wavelength modulation with λ.sub.2=405 nm (modulated) and λ.sub.1=480 nm (constant intensity); [0142] FIG. 3F is a map of the oscillating fluorescence intensity components (“out” for n=1, 3: “in” for n=2) for dual-wavelength modulation with λ.sub.1=480 nm and λ.sub.2=405 nm.

[0143] Intensity is set at I.sub.1.sup.0=1 Ein.Math.s.sup.−1.Math.m.sup.−2 (250 W/cm.sup.2) and the photochemical properties of Dronpa-2 are expressed by the following parameters: σ.sub.12,1=196 m.sup.2.Math.mol.sup.−1, σ.sub.21,2=413 m.sup.2.Math.mol.sup.−1 (for λ.sub.1=480 nm and λ.sub.2=405 nm).

[0144] We denote by K.sub.12,n.sup.0,R and θ.sub.n.sup.R the coordinates of the resonance evidenced for the concentration amplitude of the n-th order. Beyond evidencing the behavior already discussed above, the figures clearly show a diagonal shift of the resonance peaks. For a single modulation, in both evaluated cases of light modulation (λ.sub.1=480 nm and λ.sub.2=405 nm), resonance respectively shifts toward lower values of I.sub.2.sup.0/I.sub.1.sup.0 and ω/I.sub.1.sup.0 when the order n of the amplitude is increased due to the corresponding evolution of K.sub.12,n.sup.0,R and θ.sub.n.sup.R. For dual modulation, resonances display two different shifts, respectively toward higher and lower values of I.sub.2.sup.0/I.sub.1.sup.0 and lower values of ω/I.sub.1 .sup.0 when the order n of the amplitude is increased due to the corresponding evolution of K.sub.12,n.sup.0,R and θ.sub.n.sup.R.

[0145] The general theory developed above allows determining the following expression for the resonance conditions in two different limiting cases

[0146] 1. The reversible photoconversion of the fluorophore is driven by the wavelength λ.sub.1 (backward reaction is thermally driven). Using Eqs. (14), (22), (23), the resonance condition becomes:

[00028]$\begin{array}{cc}{I}_1^0=K_{12,n}^{0,R}\frac{k_{21}^{\Delta }}{{\sigma }_{12,1}}& \left(65\right)\end{array}$$\begin{array}{cc}\omega ={\theta }_n^R\left({\sigma }_{12,1}{I}_1^0+k_{21}^{\Delta }\right)& \left(66\right)\end{array}$

[0147] Both control parameters Ihd 1.sup.0 and ω are now fixed at resonance of the amplitude of the n-th order:

[00029]$\begin{array}{cc}{I}_1^0=K_{12,n}^{0,R}\frac{k_{21}^{\Delta }}{{\sigma }_{12,1}}& \left(67\right)\end{array}$$\begin{array}{cc}\omega ={\theta }_n^R\left(1+K_{12,n}^{0,R}\right)k_{21}^{\Delta }& \left(68\right)\end{array}$

[0148] In particular, Eqs. (67), (68) enable to derive the shift of the resonance peak between the orders n and m, in logarithmic scale:

[00030]$\begin{array}{cc}\log \frac{I_{1,n}^0}{I_{1,m}^0}=\log \frac{K_{12,n}^{0,R}}{K_{12,m}^{0,R}}& \left(69\right)\end{array}$$\begin{array}{cc}\log \frac{{\omega }_n}{{\omega }_m}=\log \frac{{\theta }_n^R}{{\theta }_m^R}+\log \frac{1+K_{12,n}^{0,R}}{1+K_{12,m}^{0,R}}& \left(70\right)\end{array}$

[0149] 2. The reversible photoconversion of the fluorophore is driven by two wavelengths λ.sub.1 and λ.sub.2. Using Eqs. (14), (18), (19) the resonance conditions become:

[00031]$\begin{array}{cc}\frac{I_2^0}{I_1^0}=\frac{1}{K_{12,n}^{0,R}}\frac{{\sigma }_{12,1}}{{\sigma }_{21,2}}& \left(71\right)\end{array}$$\begin{array}{cc}\omega ={\theta }_n^R\left({\sigma }_{12,1}{I}_1^0+{\sigma }_{21,2}{I}_2^0\right)& \left(72\right)\end{array}$

[0150] Once the intensity I.sub.1.sup.0 satisfies the condition of neglecting thermal contribution to the backward reaction, both control parameters I.sub.2.sup.0 and ω are fixed to optimize the resonant amplitude of the n-th order:

[00032]$\begin{array}{cc}{I}_2^0=\frac{1}{K_{12,n}^{0,R}}\frac{{\sigma }_{12,1}}{{\sigma }_{21,2}}{I}_1^0& \left(73\right)\end{array}$$\begin{array}{cc}\omega ={\theta }_n^R\left(1+\frac{1}{K_{12,n}^{0,R}}\right){\sigma }_{12,1}{I}_1^0.& \left(74\right)\end{array}$

[0151] Eqs. (73), (74) enable to derive the shift of the resonance peak between the orders n and m, in logarithmic scale:

[00033]$\begin{array}{cc}\log \frac{\frac{I_{2,n}^0}{I_{1,n}^0}}{\frac{I_{2,m}^0}{I_{1,m}^0}}=\log \frac{K_{12,m}^{0,R}}{K_{12,n}^{0,R}}& \left(75\right)\end{array}$$\begin{array}{cc}\log \frac{\frac{{\omega }_n}{I_{1,n}^0}}{\frac{{\omega }_m}{I_{1,m}^0}}=\log \frac{{\theta }_n^R}{{\theta }_m^R}+\log \frac{1+\frac{1}{K_{12,n}^{0,R}}}{1+\frac{1}{K_{12,m}^{0,R}}}& \left(76\right)\end{array}$

[0152] FIGS. 4A-4F, 5A-5F, 6A-6F display the dependence of the amplitudes of the oscillating concentration ( and ) and of the fluorescence emission ( and ) for eight investigated reversibly photoswitchable fluorescent proteins (RSFPs) Dronpa, Dronpa-2, Dronpa-3, rsFastLime, rsEGFP2, and Padron in the space of the control parameters. Dronpa, Dronpa-3, rsFastLime, rsEGFP2 exhibit negative photochromism, like Dronpa-2, while Padron exhibits positive photochromism (blue light switches it to a bright state and violet light switches it back to a dark state). All these species absorb light at 405 nm and 480 nm and emit spectrally similar fluorescence, making them hard to discriminate by prior art methods. In order to avoid overly cluttered figures, maps for different subsets of fluorophores have been displayed on different groups of figures: Dronpa, Dronpa-2, Dronpa-3 and Padron on FIGS. 4A-4F; Dronpa, Dronpa-2, rsFastLime and Padron on FIGS. 5A-5F; Dronpa, Dronpa-2, rsEGFP2 and Padron on FIGS. 6A-6F.

[0153] More precisely: [0154] FIGS. 4A, 5A and 6A are maps of the fundamental oscillating fluorescence intensity components for single-wavelength modulation with λ.sub.1=480 nm (modulated) and λ.sub.2=405 nm (constant intensity); [0155] FIGS. 4B, 5B and 6B are maps of the fundamental oscillating fluorescence intensity components for dual-wavelength modulation with λ.sub.1=480 nm and λ.sub.2=405 nm; [0156] FIGS. 4C, 5C and 6C are maps of the second harmonic fluorescence intensity components for single-wavelength modulation with λ.sub.1=480 nm (modulated) and λ.sub.2=405 nm (constant intensity); [0157] FIGS. 4D, 5D and 6D are maps of the second harmonic fluorescence intensity components for dual-wavelength modulation with λ.sub.1=480 nm and λ.sub.2=405 nm; [0158] FIGS. 4E, 5E and 6E are maps of the third harmonic fluorescence intensity components for single-wavelength modulation with λ.sub.1=480 nm (modulated) and λ.sub.2=405 nm (constant intensity); and [0159] FIGS. 4F, 5F and 6F are maps of the third harmonic fluorescence intensity components for dual-wavelength modulation with λ.sub.1=480 nm and λ.sub.2=405 nm. [0160] As expected from the resonance conditions (73), (74) the positions of the RSFPs resonances are scattered in the space of the control parameters since RSFPs exhibit distinct photochemical parameters (e.g. cross-sections for photoswitching). It is worthy of note that the resonance coordinates associated with a modulation at λ.sub.1 =480 nm with or without a modulation at λ.sub.2=405 nm are different from the coordinates associated with a single modulation at λ.sub.2=405 nm. Scattering is more pronounced for a single modulation at λ.sub.1=480 nm; [0161] In contrast, as anticipated from K.sub.12,n.sup.0,R and θ.sub.n.sup.0,R values (see Eqs. (74), (76)) the shift of the resonance position at different orders n in the space of the control parameters does not depend on the RSFPs nor on the modulated wavelength. Hence the relative position of the RSFP resonances does not evolve with the order n. Note that the n-dependence of the resonance amplitude is the same in the (I.sub.2.sup.0/I.sub.1.sup.0,ω/I.sub.1.sup.0) and (K.sub.12.sup.0,θ) spaces; [0162] For all modulation cases, the resonance bandwidth of the concentration amplitude slowly decreases from the first to the fourth order. The effect is more visible for a single modulation at λ.sub.1=480 nm or dual modulation and from the third order on; [0163] The behavior of the fluorescence intensity amplitude for a single modulation at λ.sub.2=405 nm is similar to the concentration amplitude. In contrast, the decrease of the bandwidth is more visible when the order increases with a single modulation at λ.sub.1=480 nm. The second order brings already a better RSFP discrimination than the first order, and the map is clarified from the third order on (for instance Dronpa stands out). For the dual modulation, the second order display broader bandwidth than the first order, but the third order brings a better RSFP discrimination than the first order. However at all orders the single modulation at λ.sub.1=480 nm brings more RSFP discrimination than the dual modulation; [0164] The figures clearly show that for the single modulation at λ.sub.1=480 nm the decreased resonance overlapping of the third order amplitude of the fluorescence can be used to discriminate four RSFPs: Dronpa, Dronpa-2, Dronpa-3 or rsFastLime, and Padron which exhibits a fluorescence amplitude of opposite sign with respect to the three others. [0165] The utilization of high order responses of the fluorescence emission to the light modulation is favorable to discriminate distinct RSFPs. However, since the amplitude of the n-order response drops with n, a compromise has to be found between the bandwidth (higher orders enable us a better discrimination) and the amplitude (lower orders exhibit a higher amplitude). [0166] Unlike the fundamental components used in OPIOM and Speed OPIOM, harmonic components exhibit, in addition to resonant peaks, zero-amplitude lines which allow selectively “turning off” the contribution of a RSFP at a particular order. Therefore, in order to discriminate between two species, a user has the choice of maximizing the contribution(s) from one species, extinguishing the contribution(s) from another one or finding a suitable compromise to maximize the dynamic contrast. Several oscillating (fundamental and harmonic) contributions may be simultaneously used for detection. Several reversibly photoswitchable species may be discriminated by performing several measurements with different illumination conditions e.g. resonant for the different species to be detected. [0167] Whatever the observable used, signal is proportional to the concentration of the detected species. This can be a drawback, as the out-of-resonance signal of an abundant species may interfere with the resonant signal of the less abundant species to be detected. However, as explained above, the existence of zero-amplitude lines allows “turning off” the interfering signal by a suitable choice of the illumination parameters. Moreover, this makes it possible to estimate absolute or relative concentration of reversibly photoswitchable species from several measurements performed in different illumination conditions. For instance, this can be performed by machine learning (e.g. neural networks, regression analysis . . . ), using samples of known composition for training, or unmixing algorithms.

[0168] FIGS. 7A-7I demonstrate the discriminative power of the inventive method applied to epifluorescence microscopy.

[0169] Four chambers of a microfluidic device were respectively filled with solutions of Dronpa, rsFastLime, Dronpa-2, and EGFP; on the figures, the chambers are labeled D, FL, D2 and E, respectively. Dronpa, rsFastLime and Dronpa-2 are spectrally similar RSFPs, which the inventors failed to discriminate using Speed OPIOM, while EGFP is a spectrally similar non-photochemically active fluorophores. The microfluidic device was illuminated using sinusoidally modulating light at λ.sub.1=480 nm and constant light at λ.sub.2=405 nm at the resonance conditions for the three RSFPs. More precisely: [0170] FIGS. 7A, 7B and 7C show the images obtained by extracting the quadrature fundamental component of the fluorescence intensity, at the resonance conditions for Dronpa (FIG. 7A), Dronpa-2 (FIG. 7B) and rsFastLime (FIG. 7C). Whereas the image recorded at Dronpa-2 resonance could almost eliminate the interfering contributions of the two other RSFPs, the images recorded at resonances of Dronpa and rsFastLime encountered strong interferences from rsFastLime and Dronpa-2 respectively. In fact, such a result was anticipated from FIG. 5A, which evidences that significantly overlaps for the three RSFPs. Hence, use of the first-order component alone (i.e. Speed OPIOM) is not appropriate to discriminate the three RSFPs. [0171] FIGS. 7D, 7E and 7F show the images obtained by extracting the second-order in-phase component of the fluorescence intensity, at the resonance conditions for Dronpa (FIG. 7D), Dronpa-2 (FIG. 7E) and rsFastLime (FIG. 7F). It can be seen that the EGFP signal was again efficiently removed from all the images and, most importantly, that the three RSFPs can be selectively targeted at their respective resonance. Beyond the drop of the resonance bandwidth, it can be estimated that this achievement probably originates from the 0 lines observed in the discrimination map of the second harmonic component (see FIG. 7C). Indeed the 0 lines of rsFastLime and Dronpa-2 respectively go through resonance of Dronpa and rsFastLime. Hence the second harmonic component images at resonances of Dronpa and rsFastLime eliminate the interferences with rsFastLime and Dronpa-2 respectively observed using the first-order contribution. [0172] FIGS. 7G, 7H and 7I show the images obtained by extracting the third-order quadrature component of the fluorescence intensity, at the resonance conditions for Dronpa (FIG. 7G), Dronpa-2 (FIG. 7H) and rsFastLime (FIG. 7I). The theoretical maps suggest a behavior much similar to the one obtained with the second-order component, but the signal-to-noise ratio proved too low to reach any definitive conclusion.

[0173] Since a reversibly photoswitchable species possesses singular kinetic properties of photoswitching (and correspondingly its own resonant conditions), when the spatial profile of modulated illumination (I.sub.1,I.sub.2) is not homogeneous, the resonance conditions for a specific harmonic (or fundamental) contribution can only be met at restricted locations, which filters out the fluorescence contributions from out-of-resonance regions. Hence the inventive method is endowed with intrinsic optical sectioning properties. The plot of FIG. 8 shows the evolution of the spatially-integrated signal intensity with the thickness of the sample (triangles correspond to averaged intensity, squares to the first-order quadrature component and disks to the second order in-phase contribution). It can be seen that, beyond 15 μm, the first- and second-order components reach a plateau, demonstrating that their responses originate from a spatially restricted section of the sample. The average fluorescence intensity, instead, continues increasing (albeit sub-linearly) with sample thickness.

[0174] It will be easily understood that the sectioning power of the inventive method may be very useful for imaging applications, e.g. in biology.

[0175] FIG. 9 shows the block diagram and optical layout of an apparatus for carrying out the invention according to an embodiment of the invention. This apparatus is based on an epi-fluorescence microscope.

[0176] The apparatus integrates two colored LEDs: LS1 centered at 480 nm and LS2, centered at 405 nm, as excitation sources, in order to generate a large focal field around the focal plane. Optical filters OF1, OF2 are used to reduce the spectral width of the generated light. Each LED is triggered by a LED driver (DC 4104, Thorlabs, NJ, US), not represented. In turn the led drivers are controlled by a waveform generator WFG, introducing the desired intensity modulation(s).

[0177] Light from each source is first collimated by lens L1 (AC254-200-A, f=200 mm, Thorlabs, NJ, US) and then reduced by a quasi-afocal system consisting of two lenses (L2: LA1134-A, f=60 mm, L3: LC2679-A, f=−30 mm, Thorlabs, Newton, NJ). The reduced beams are combined by a dichroic mirror (T425LPXR, Chroma, Bellows Falls, VT) and then focused by a convergent lens L4 (LA1131-A, f=50 mm, Thorlabs, NJ, US). The focal plane of L4 is conjugated to the focal plane of the objective (UPLSAP060XW, 60×, NA 1.2, Olympus Corporation, Tokyo, Japan), with the help of the tube lens system L5 and L6 (AC254-200-A, f=200 mm, Thorlabs, NJ, US). The objective along with the tube lenses L5 and L6 gives a total magnitude of 100:3, forming an illuminated area of 15 .sub.lim at the sample. The position of the imaging plan after L4 can be effectively changed by slightly moving the element L3 back and forth along the optical axis, which makes it possible to defocus each wavelength independently from the focal plane of the objective. According to software simulation, a displacement of the element L3 by 1 mm gives rise to a defocus of 4 μm at the focal plane of the objective. The fluorescence signal from the sample is separated from the excitation lights with another dichroic mirror (T505LPXR, Chroma, Bellows Falls, VT) and then recorded by a CCD camera CAM (Luca-R Andor Technology, Belfast, UK) serving as light detector.

[0178] A data processing device such as a computer DPD controls the operation of the waveform generation, in order to determine the desired illumination conditions. Moreover, the computer receives and processes the signal generated by the camera, representative of the fluorescent light intensity, and processes them by extracting the oscillating components used for detecting reversibly photoswitchable species in the sample S.

[0179] In many applications, the use of a scanning microscope may be preferred to that of a wide-field microscope as that of FIG. 9. However, Illumination is by essence discontinuous (pulsed) in scanning microscopy, while the invention requires a continuous, preferably sinusoidal modulation (or at least a square-wave modulation with a duty-cycle near to 0.5). However, it is possible to define experimental conditions for which illumination can be nevertheless considered sinusoidal. We define the time t.sub.c as the duration of a light pulse on a pixel and t.sub.d the time between two pulses. A rectangular wave of period t.sub.c+t.sub.d varying between 0 and 1 is decomposed into Fourier series according to:

[00034]$\begin{array}{cc}R\left(t\right)=\frac{t_c}{t_c+t_d}+{.Math.}_{n=1}^{\infty }\frac{2}{n\pi }\sin \left(\frac{n\Omega t_c}{2}\right)\cos \left(n\Omega t\right)& \left(77\right)\end{array}$

[0180] where Ω=2π/(t.sub.c+t.sub.d) is the fundamental angular frequency. We choose an illumination at the wavelength λ.sub.i with i=1,2, modeled by the product of the rectangular wave R(t) and a sinusoidal envelope of angular frequency ω oscillating around the mean value I.sub.i.sup.0′,

I.sub.i(t)=I.sub.i.sup.0′R(t)(1+α sin(ωt))

[0181] where α is the amplitude of oscillation. Imposing Ω=kω where k is an integer, we find:

[00035]$\begin{array}{cc}{I}_i\left(t\right)=I_i^{0^{\prime }}\frac{t_c}{t_c+t_d}\left(1+\alpha \sin \left(\omega t\right)\right)+I^{0^{\prime }}{.Math.}_{n=1}^{\infty }\frac{1}{n\pi }\sin \left(\frac{\mathrm{nk}\omega t_c}{2}\right)[2\cos \left(\mathrm{nk}\omega t\right)+\alpha \sin \left(\mathrm{nk}+1\right)\omega t)+\alpha \sin \left(\left(\mathrm{nk}-1\right)\omega t\right)]& \left(78\right)\end{array}$

[0182] The desired term I.sub.i.sup.0)(1+α sin(ωt)) with I.sub.i.sup.0=I.sub.i.sup.0′t.sub.c/(t.sub.c+t.sub.d) at the fundamental frequency is not polluted by the first harmonics obtained for n=1 provided that k is sufficiently large. Adapting the inventive method to scanning microscopy requires Ω>>ω (more precisely, ≥10ω), where Ω is the frequency of successive scans on a pixel and ω is the frequency of the sinusoidal envelope of the illumination.

[0183] FIG. 10 shows a pulsed intensity profile suitable for carrying out the invention. The result is valid for light modulation at a single wavelength or two wavelengths. The mean light intensity perceived by a pixel is I.sub.i.sup.0=I.sub.i.sup.0′t.sub.c/(t.sub.c+t.sub.d). The theoretical computations of the concentration and fluorescence amplitudes obtained for wide-field illumination are valid using this definition of I.sub.i.sup.0 for the mean light intensity. In addition, light intensities delivered by a scanning microscope are stronger than in wide-field microscopy. The condition t.sub.c≤τ.sub.12.sup.0 where τ.sub.12.sup.0 depends on light intensity ensures that kinetics is controlled by photochemistry and correctly described by a two-state model. Moreover thermal steps ignored by the two-state model occur on a fast time scale τ.sup.Δ. The time τ.sub.d must be larger than x for the thermal steps to have time to relax between two pulses.

[0184] Before beginning modulation the system has to be prepared in the desired steady state associated with the mean light intensities I.sub.1.sup.0 and I.sub.2.sup.0 used during modulation. Preparation includes a series of pulses at the two wavelengths of duration t.sub.c≤τ.sub.12.sup.0 which imposes a condition on the mean light intensities such that only fractions of 1 and 2 react during t.sub.c. The pulses must by separated by a time t.sub.d larger than σ.sup.Δ for the thermal steps to relax between two pulses. The number p of pulses is such that pt.sub.c≥5τ.sub.12.sup.0 for the steady-state to be reached. Superimposing modulation during the preparation has the same effect and has the advantage of keeping the same protocol during the entire experiment. In this case, at least the first p pulses are necessary to prepare the system in the desired steady-state.

[0185] FIG. 11 schematically represents a scanning microscope SM coupled to a data processing device such as a computer DPD which controls the light intensity profiles and processes signals representing collected light intensity by extracting the oscillating components used for detecting reversibly photoswitchable species in the sample S.

[0186] The invention is not limited to the case of a microscope—scanning or wide-field—of observing the sample. For instance, sample microscope SM may be replaced by a more general scanning optical device. For instance, WO2015/075209 describes a macroscopic scanning imaging device suitable for carrying out the OPIOM method, comprising an electro-optical modulator for modulating the intensity of a laser beam and a pair of steerable mirror for scanning a distant target, e.g. a ground surface area, with the modulated beam. The sample is supposed to contain at least one reversibly photoswitchable fluorescent species, and its fluorescence emission is remotely collected by an objective of a camera. This apparatus can easily be applied to the implementation of the inventive method.

[0187] Up to now, only the case of resonant detection—i.e. where the illumination intensities and angular frequency are selected to match the resonant condition of at least one oscillating component of an observable—have been considered. However, the amplitudes of the Fourier decomposition of the fluorescence response to light modulation constitute the kinetic footprint of a given resonant photoswitchable species, and the amplitudes which vanish far from resonance conditions also offer good discrimination properties.

[0188] For each resonant photoswitchable species two sets of experimental conditions are used. In the case of modulation at λ.sub.1=480 nm, the first-order amplitude displays one resonance R.sub.1 and the second-order amplitude has a resonance close to R.sub.1 and another one, R.sub.2, of weaker amplitude. The experiments give access to the out-of-phase first-order amplitudes (R.sub.1) and (R.sub.2), and the in-phase second-order amplitudes (R.sub.1) and (R.sub.2) evaluated at the resonance conditions R.sub.1 and R.sub.2, respectively.

[0189] It was chosen to scale the amplitudes by the non-resonant term (R.sub.2) in order to eliminate their dependence on concentration, brightness, and the collection factor of the objective. In this work the kinetic footprint is reduced to the set (l.sub.1,l.sub.2,l.sub.3) with l.sub.1= (R.sub.1)/ (R.sub.2), l.sub.2= (R.sub.1)/ (R.sub.2), and l.sub.3= (R.sub.2)/ (R.sub.2) in a 3-dimension space but it can be generalized to an m-dimension space if amplitudes of higher order are detectable. Theoretical computation or experiments performed for a given RSF i lead to the reference kinetic footprint {l.sub.ki} with k=1, . . . ,m. The kinetic footprint M{l.sub.kj} of an unknown resonant photoswitchable species j is compared to the reference i using the distance:

d.sub.ij=√{square root over (Σ.sub.k=1.sup.m(l.sub.ki−l.sub.kj).sup.2)}   (79)

Discrimination between two resonant photoswitchable species is possible if their distance is larger than the cutoff distance d.sub.c imposed by the experimental accuracy Ad.sub.ii on the distance d.sub.ij induced by the uncertainty Δl.sub.ki on the kinetic footprint {l.sub.ki}. Differentiating Eq. (79), we find

[00036]$\begin{array}{cc}\Delta d_{\mathrm{ij}}=\frac{1}{d_{\mathrm{ij}}}{.Math.}_{k=1}^m.Math.l_{\mathrm{ki}}-l_{\mathrm{kj}}.Math.\Delta .Math.l_{\mathrm{ki}}-l_{\mathrm{kj}}.Math..& \left(80\right)\end{array}$

Introducing the maximum uncertainty M=max.sub.k,i(Δl.sub.ki) on the determination of the kinetic footprint {l.sub.ki}, we have Δ|l.sub.ki−l.sub.ki|≤2M. Hence Eq. (80) reads

Δd.sub.ij≤2MX   (81)

with X.sup.2=1+2Σ.sub.k=1.sup.mΣ.sub.k′=k+1.sup.m|l.sub.ki−l.sub.kj∥l.sub.k+i−l.sub.k′j|/d.sub.ij.sup.2. Regardless of the value of |l.sub.ki−l.sub.kj|, the inequality:

2Σ.sub.k=1.sup.mΣ.sub.k′=k+1.sup.m|l.sub.ki−l.sub.kj∥l.sub.k′i−l.sub.k′j|≤(m−1)Σ.sub.k=1.sup.m(l.sub.ki−l.sub.kj).sup.2   (82)

is observed, leading to X≤√{square root over (m)}. Using Eq. (81) we find Δd.sub.ij≤d.sub.c where the cutoff distance obeys

d.sub.c=2M√{square root over (m)}   (83)

[0190] The cutoff distance between two resonant photoswitchable species depends on two parameters, the number m of discriminating dimensions and the uncertainty M on the determination of the kinetic footprint. An imaged RSF i is recognized as the resonant photoswitchable species j if the distance d.sub.ij is smaller than d.sub.c.

[0191] Equation (79) only constitutes an example, and the discriminating distance may be defined otherwise.

[0192] Like resonant detection, a plurality of nonresonant measurements may be used, preferably in conjunction with a machine-learning algorithm, for estimating absolute or relative concentrations of reversibly photoswitchable species from several measurements performed in different illumination conditions.

[0193] The invention has been described with reference to a number of examples which, however, are not limiting. Different embodiments are possible; for instance, illumination at more than two wavelengths may be used; if light at more than one wavelength is modulated, different wavelengths may correspond to different modulation frequencies and/or have different phase relationships.

[0194] The invention specifically relates to a plurality of methods using harmonic components of the measured signal (i.e. components having an angular frequency n.Math.ω with n is an integer greater than 1 and ω is the “fundamental” modulation frequency of the illuminating light), either alone or together with the fundamental component of the signal. However, it is also possible to use, in all the embodiments described above, the out-of-phase fundamental component of the signal alone, without harmonic components. In this case, use of an on-off modulation of the illuminating light may be preferred.

[0195] This is particularly interesting in the case of pulsed illumination, and more particularly of a scanning microscopy, or more generally a scanning optical device, as described above in reference to FIGS. 10 and 11.

[0196] Therefore, it is hereby disclosed a method for detecting a reversibly photoswitchable chemical species in a sample, comprising the steps of:

[0197] a) illuminating the sample with a first light at a first wavelength suitable to be absorbed by the chemical species triggering a reaction affecting at least one optical property of the chemical species, said first light being periodically-modulated at a fundamental modulation frequency;

[0198] b) measuring the evolution of the optical property of the chemical species;

[0199] c) extracting a quadrature component at the frequency of the fundamental modulation frequency of a signal representing said evolution; and

[0200] d) using the extracted component for detecting the chemical species;

[0201] wherein light is modulated by a sinusoid at the fundamental modulation frequency multiplied by a pulse train at a repetition frequency which is a multiple of, and at least ten times larger than, the fundamental modulation frequency.

[0202] According to different, advantageous embodiments of such a method: [0203] The average intensity of the light and the fundamental modulation frequency may be chosen so to maximize the amplitude of said extracted component. [0204] Step a) may further comprise illuminating the sample with a second light at a second wavelength, different from the first wavelength and suitable to be absorbed by the chemical species triggering said or another reaction affecting said optical property of the chemical species, said second light either having a constant intensity or being periodically modulated at the fundamental modulation frequency in phase opposition with the first light (or by a sinusoid at the fundamental modulation frequency and in phase opposition with the first light, multiplied by a pulse train at a repetition frequency which is a multiple of, and at least ten times larger than, the fundamental modulation frequency). In this case, the ratio of the average intensities of the first and second light and the fundamental modulation frequency are preferably chosen so as to maximize the amplitude of said or at least one said extracted component. [0205] Steps a), b) and c) may be repeated a plurality of times for different illumination conditions, corresponding to different values of the light intensity or intensities and of the fundamental modulation frequency and step d) may then comprise using a plurality of signal components corresponding to said different illumination conditions for discriminating between a plurality of chemical species. [0206] A machine learning or unmixing method may be applied to said plurality of signal components for estimating absolute or relative concentration of said chemical species.

REFERENCES

[0207] [Lakowicz 1992] J. R. Lakowicz, H. Szmacinski, K. Nowaczyk, K. W. Berndt, and M. L. Johnson. Fluorescence lifetime imaging. Anal. Biochem., 202:316-330, 1992.

[0208] [Marriott 2008] G. Marriott, S. Mao, T. Sakata, J. Ran, D. K. Jackson, C. Petchprayoon, T. J. Gomez, E. Warp, 0. Tulyathan, H. L. Aaron, E. Y. Isacoff, and Y. Yan. Optical lock-in detection imaging microscopy for contrast enhanced imaging in living cells. Proc. Natl. Acad. Sci. U. S. A., 105:17789-17794, 2008.

[0209] [Richards 2010] C. I. Richards, J.-C. Hsiang, and R. M. Dickson. Synchronously amplified fluorescence image recovery (SAFIRe). J. Phys. Chem. B, 114:660-665, 2010.

[0210] [Widengren 2010] J.Widengren. Fluorescence-based transient state monitoring for biomolecular spectroscopy and imaging. J. R. Soc. Interface, 7:1135-1144, 2010.

[0211] [Quérard 2015] J. Qu6rard, T.-Z. Markus, M.-A. Plamont, C. Gauron, P. Wang, A. Espagne, M. Volovitch, S. Vriz,V. Croquette, A. Gautier, T. Le Saux, and L. Jullien. Photoswitching kinetics and phase-sensitive detection add discriminative dimensions for selective fluorescence imaging. Angew. Chem. Int. Ed., 127:2671-2675, 2015.

[0212] [Quérard 2017] J. Querard, R. Zhang, Z. Kelemen, M.-A. Plamont, X. Xie, R. Chouket, I. Roemgens, Y. Korepina, S. Albright, E. Ipendey, M. Volovitch, H. L. Sladitschek, P. Neveu, L. Gissot, A. Gautier, J.-D. Faure, V. Croquette, T. Le Saux, and L. Jullien. Resonant out-of-phase fluorescence microscopy and remote imaging overcome spectral limitations. Nat. Comm., 8:969, 2017.

[0213] [Rowe 2002] H. M. Rowe, S. P. Chan, J. N. Demas, and B. A. DeGraff. Elimination of Fluorescence and Scattering Backgrounds in Luminescence Lifetime Measurements Using Gated-Phase Fluorometry, Anal. Chem. 2002, 74 , 4821-4827.

[0214] [Jenkins 2015] P. Jenkins, M. A. Naivar, and J. P. Houston. Toward the measurement of multiple fluorescence lifetimes in flow cytometry: maximizing multi-harmonic content from cells and microspheres. J. Biophotonics 2015 November; 8 (11-12): 908-917.