Abstract:
The present invention concerns a confocal laser scanning microscopy apparatus, comprising: means for emitting a laser beam; means for scanning this laser beam in at least two directions onto an observed sample; means for generating a non linear light signal from the transmitted laser light, these non linear light signal generating means being disposed in the light path between the observed sample and detecting means which are adapted for detecting said non linear light signal.

Description:
[0001]     The present invention claims priority to U.S. 60/539254, which was filed on 26 Jan. 2004. 
     
    
     BACKGROUND OF THE INVENTION  
       [0002]     The present invention concerns a confocal laser scanning microscopy apparatus, notably a transmission confocal laser scanning microscopy apparatus.  
         [0003]     Such microscopy apparatus exist and are for example described in the U.S. Pat. No. 3,013,467, or in “Theory and Practice of Scanning Optical Microscopy”, Academic press, London, 1984, by T. Wilson and C. Sheppard.  
         [0004]     One of the drawbacks of such microscopy apparatus is the need for de-scanning the laser beam since such aparatus comprise a pinhole in the detection path.  
       SUMMARY OF THE INVENTION  
       [0005]     The present invention improves such known microscopy apparatus by providing a confocal laser scanning microscopy apparatus, comprising 
        means for emitting a laser beam;     means for scanning this laser beam in at least two directions onto an observed sample;     means for generating a non linear light signal from the transmitted laser light, these non linear light signal generating means being disposed in the light path between the observed sample and detecting means which are adapted for detecting said non linear light signal.        
 
         [0009]     A confocal laser scanning microscopy apparatus according to the invention may further comprise one or more of the following features: 
        the generating means are second harmonic laser light generating means, and the detecting means are adapted for detecting said second harmonic laser light;     the generating means comprise a second harmonic generating crystal;     the crystal is a Type I, lithium triborate crystal;     the thickness of the crystal along the light path is roughly equal to 200 micrometers;     the lateral dimensions of the crystal are sufficiently large to cover the area scanned by the laser beam;     means for scanning are adapted for scanning the laser in two lateral and one longitudinal directions;     the apparatus further comprises a short-pass filter disposed in the light path immediately behind the non linear light signal generating means such that only non linear light signal is detected;     the non linear light signal generating means comprise a two-photon excited fluorophore;     the non linear light signal generating means comprise a wide-bandgap semiconductor; and     the emitting means are adapted for emitting laser pulses.       
 
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0020]     The invention will be more clearly understood from the following description, given by way of example only, with reference to the accompanying drawings in which: 
         FIG. 1  represents an experimental layout (where MO=microscope objective. Filter transmits only SHG light);      FIG. 2  represents plots of measured &lt;SHG&gt;/SHG 0  and &lt;P&gt; 2 /P 0   2  for z-scan of a 530 nm latex bead (averaged over 5 scans) (Traces are derived from a Gaussian approximation model. Panels illustrate non-averaged SHG (a) and direct P 2  (b) x-y scan of a bead. Scale bars=5_m);      FIG. 3  represents plots of measured &lt;SHG&gt; and &lt;P&gt; 2  for a z-scan of a 170_m thick agarose slab of 1_m latex beads, and theoretical trace derived from Gaussian approximation, normalized to arbitrary units;      FIG. 4  represents x-y images of an onion slice under a 200_m agarose slab of 1_m latex beads using SHG detection (left), and direct P 2  detection (right) (Scale bars=100_m);      FIG. 5  represents a parfocal unit-magnification configuration of a quadratic detection ACM (A scanning laser beam is focused into a sample and re-focused onto a thin SHG crystal. The resulting SHG is isolated (filter) and detected by a photomultiplier tube (PMT));      FIG. 6  represents quadratic detection ACM point spread functions (First order signal responses produced by a point object that is purely absorbing (a: axial; b: radial), or purely phase-shifting (c: axial; radial response is null));      FIG. 7  represents (a) a configuration in which a point object is embedded at a depth L A  within a turbid slab of thickness L and (b) an equivalent configuration in which phase variations provoked by slab are projected into appropriate lens pupils;      FIG. 8  represents afilter function H L  as a function of normalized spatial frequency ξ d  (Ballistic light is attenuated by a factor exp(−Γ L   (φ) (0)) (dotted line). Non-ballistic light is transmitted only below the cut-off frequency ξ 3 dB );      FIG. 9  represents a background ACM signal (closed circles) obtained experimentally when translating a scattering slab in the z direction (slab is centered when z slab =0) (The slab is composed of 1 μm latex beads in agarose (λ=870 nm, I s =39 μm, I* s =490 μm, L=340 μm). Theoretical trace (solid line) from Eq. 35 is shown for comparison);      FIG. 10  represents (a) a Simultaneous TPEF (left) and ACM (right) images of a 500 nm fluorescent latex bead, (b) a corresponding ACM contrast (squares) produced by non-fluorescent beads as a function of penetration depth L A  in a scattering medium (λ=870 nm, I s =126 μm, I* s =1800 μm, L=700 μm) consisting of 1 μm diameter non-fluorescent beads, and sparsely distributed 0.5 μm diameter fluorescent beads, (c) an ACM contrast as a function of depth for a different scattering medium (λ=870 nm, I s =39 μm, I* s =490 μm, L=340 μm) (Solid lines are theoretical traces for ACM contrast decay as a function of depth (no free parameters). For reference we illustrate decays of TPEF signal (circles)).       
 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0031]     The present invention concerns a transmission-mode confocal scanning laser microscope system based on the use of second harmonic generation (SHG) for signal detection. This method exploits the quadratic intensity dependence of SHG to preferentially reveal unscattered signal light and reject out-of-focus scattered background. The SHG crystal plays the role of a virtual pinhole that remains self-aligned without a need for de-scanning.  
         [0032]     Usually, confocal laser scanning microscopy (CLSM) is based on the use of a pinhole in the detection path to provide 3-dimensional image resolution and enhanced background rejection. In the usual CLSM implementation, detected light is de-scanned so that the pinhole effectively tracks the position of the laser focus at the sample. Such de-scanning is readily accomplished in a reflection configuration by retracing the signal path through the laser scanning optics. In a transmission configuration, however, de-scanning is technically much more difficult and typically requires the use of a second synchronized scanning system or of an elaborate beam path to re-direct the transmitted light into the backward direction. We present a simple technique to accomplish self-aligned de-scanning in a transmission CLSM based on signal conversion with a second harmonic generation (SHG) means comprising a crystal.  
         [0033]     In standard transmission CLSM, laser light transmitted through the sample is focused onto a pinhole of area A P  before detection. If the transmitted light has power P and is distributed over a characteristic area A at the pinhole plane, the detected power scales as PA P /A (assuming A P &lt;A). In our method, the pinhole is replaced by a thin nonlinear crystal and only SHG is detected. Because SHG scales quadratically with incident intensity, the resulting signal scales approximately as P 2 /A. In both cases, the detected signal scales inversely with A, implying that out-of-focus light at the aperture (or crystal) plane is rejected. A distinct advantage of using a SHG crystal instead of a pinhole is that it has a large area, allowing it to play the role of an aperture even when the transmitted signal light is not de-scanned. That is, the crystal may be thought of as a self-aligned virtual pinhole.  
         [0034]     We demonstrate the above principle with the experimental setup shown in  FIG. 1 . We use a mode-locked Ti:Sapphire laser (Spectra-Physics) to generate laser pulses at 860 nm wavelength, ˜100 fs duration, and 82 MHz repetition rate, that are focused into a target sample with a 60 — 0.9 N.A. water-immersion objective (Olympus; focal waist w 0 =0.5 μm). The transmitted light is collected with an identical objective and refocused onto a Type I, 200 μm thick lithium triborate crystal (Castech). The total magnification factor, M, from the sample to the crystal is approximately 30, leading to a confocal parameter at the crystal of about 800 μm (ie. the crystal is thin relative to this confocal parameter). The laser power incident on the sample, P 0 , is typically 10 mW. The laser beam is raster scanned in the x-y direction with galvanometer mounted mirrors, and the sample is scanned in the z direction with a motorized translation stage.  
         [0035]     To begin, we consider the signal obtained from a single isolated scatterer, a latex bead, which we scan in the z direction, yielding the axial point-spread-function of our apparatus (z=0 denotes the focal plane). We first observe that, in our imaging configuration, the phase of the scattered light at the crystal plane is approximately in quadrature with that of the unscattered light, independently of the bead position z. This result is a consequence of the cumulative Guoy shifts incurred by both the scattered and unscattered beams before they reach the crystal plane (when z&gt;&gt;w 0   2 /λ the scattered beam incurs no net Guoy shift). To a good approximation, the total intensity incident at the crystal plane is then simply given by the sum of the respective unscattered and scattered intensities: 
 
 I ( r,z )= P   0 ((1−ε z ) W   0 ( r )+ε z   ηW   S ( r,z )),   (1) 
 
 where r is the radius from the optical axis at the crystal plane (we assume cylindrical symmetry), ε z P 0  is the total power scattered by the bead, η is the fraction of this power accepted by the microscope exit pupil (defined here by the collection objective), and W 0,S (r) are flux densities normalized so that 2π∫W 0,S (r)r dr=1. These functions allow us to define the characteristic areas A 1J =(2π∫W 1 (r)W J (r)r dr) −1 . Since the SHG produced by the crystal is proportional to ∫I 2 (r)r dr, we conclude that  
               SHG     SHG   0       =         (     1   -     ɛ   z       )     2     +     2   ⁢       ɛ   z     ⁡     (     1   -     ɛ   z       )       ⁢   η   ⁢       A   00       A     0   ⁢   S           +       ɛ   z   2     ⁢     η   2     ⁢         A   00       A   SS       .                 (   2   )             
 
         [0036]     Several comments are in order. First, A 00 /A 0S  and A 00 /A SS  are smaller than 1, since W 0 (r) corresponds to a diffraction limited intensity profile. Second, it is apparent that as far as scattered light is concerned, A 00 /A 0S  and A 00 /A SS  play an aperturing role similar to that of the microscope exit pupil. The smaller the values of A 00 /A 0S  and A 00 /A SS  both of which depend on z, the less the scattered light contributes to the SHG signal (ie. the more it is rejected). Finally, for purposes of comparison, we note that if the SHG crystal were removed and the power directly detected and squared, the expression for P 2 /P 0   2  would be given by (2) with the replacements A 00 /A 0S →1 and A 00 /A SS →1. In other words, direct detection of power provides no scattered light rejection beyond that of the exit pupil.  
         [0037]      FIG. 2  illustrates SHG/SHG 0  and P 2 /P 0   2  for a z-scan of a 530 nm diameter bead. In both cases, the presence of the bead is recognized as a reduction in unscattered laser power (first term in Eq. 2). This reduction is undermined by the concurrent detection of forward-directed scattered power that is transmitted through the exit pupil, which we refer to as background (second and third terms in Eq. 2). Because background rejection is more efficient with SHG than with direct detection, our method leads to a more highly contrasted bead signal.  
         [0038]     The parameters in Eq. 2 can be roughly estimated in a paraxial approximation by assuming that W 0 (r) and W S (r,z) are Gaussian in profile, leading to A 00 ≈M 2 πw 0   2 , A SS ≈M 2 πw S   2 (1+(λz/π S   2 ) 2 ), and A 0S =(A 00 +A SS )/2, where w S  is the effective bead radius as it appears through the microscope exit pupil. The scattering parameter ε z  is dependent on z since it depends on the laser intensity incident on the bead. Denoting σ as the bead scattering cross-section, then ε z ≈σ/U z , where  
         U   z     =       1   2     ⁢   π   ⁢           ⁢       w   0   2     ⁡     (     1   +       (     λ   ⁢           ⁢     z   /   π     ⁢           ⁢     w   0   2       )     2       )             
 
 is the effective area of the laser beam at the axial position z. The pupil transmission η, on the other hand, is approximately independent of z for small z&#39;s. The following estimates are derived from Mie theory: σ≈π×(0.15 μm) 2 , w S ≈0.83 μm, and η≈0.75. As is evident from a comparison with experimental data, our Gaussian approximation is overly simplistic and cannot account for the observed ringing in the SHG trace, presumably caused by pupil apodization. Nevertheless, it illustrates a salient principle of our microscopy technique, namely that A 00 /A 0S  and A 00 /A SS  are smaller than 1, leading here to an improvement in signal contrast with SHG detection. 
 
         [0039]     To demonstrate that virtual pinhole microscopy with SHG detection also leads to improved out-of-focus background rejection, we acquire a z-stack of x-y scans of a slab of 1 mm latex beads suspended in 0.3% agarose (number concentration N=0.0071 μm −3 ; slab thickness L=170 μm). Since ε z  fluctuates randomly for different x-y-z positions in the slab, we write ε z =&lt;ε z &gt;+δε z  where the brackets refer to the average over an ensemble of x-y scans. If the scattering beads are randomly distributed in the slab and δz is chosen large enough so that ε z  and ε z+δz  are uncorrelated, then &lt;ε z &gt;≈Nσδz and &lt;δε z   2 &gt;≈Nσ 2 δz/U z . Though these last expressions require σ&lt;&lt;U z , meaning their validity breaks down somewhat in the immediate vicinity of the focal plane, we infer that &lt;ε z &gt;,&lt;δε z   2 &gt;&lt;&lt;1 throughout most of the sample. Eq. (2) then leads to the approximation:  
                   〈   SHG   〉       SHG   0       ≈       ∏     z   =       z   slab     -     L   /   2             z   slab     +     L   /   2         ⁢           ⁢     (     1   -     2   ⁢     〈     ɛ   z     〉     ⁢     (     1   -       η   z     ⁢       A   00       A     0   ⁢   S             )       +     
     ⁢       〈     δɛ   z   2     〉     ⁢     (     1   -     2   ⁢     η   z     ⁢       A   00       A     0   ⁢   S           +       η   z   2     ⁢       A   00       A   SS           )         )         ,           (   3   )             
 
 where z slab  is the axial location of the slab center and η z  is no longer assumed to be constant since |z| can be large. Expression 3 is readily evaluated with the substitution  
                   ∏   z             ⁢           ⁢     (     1   -       f   ⁡     (   z   )       ⁢   δ   ⁢           ⁢   z       )       ≈     exp   ⁡     (     -       ∫   z     ⁢       f   ⁡     (   z   )       ⁢   δ   ⁢           ⁢   z         )         ,       yielding   ⁢       〈   SHG   〉       SHG   0         ≈     exp   ⁡     (       -   2     ⁢   N   ⁢           ⁢   σ   ⁢       ∫       z   slab     -     L   /   2           z   slab     +     L   /   2         ⁢       (     1   -       η   z     ⁢           ⁢       A   00       A     0   ⁢   S           -     σ     2   ⁢     U   z           )     ⁢   δ   ⁢           ⁢   z         )         ,           (   4   )             
 
 where only the dominant terms have been kept. Relation 4 can be analytically expressed when using the Gaussian approximation. Again for comparison, we note that in the case of direct detection  
                   〈   P   〉     2       P   0   2       ≈       exp   ⁡     (       -   2     ⁢   N   ⁢           ⁢   σ   ⁢       ∫       z   slab     -     L   /   2           z   slab     +     L   /   2         ⁢       (     1   -     η   z       )     ⁢   δ   ⁢           ⁢   z         )       .             (   5   )             
 
         [0040]     In particular, we observe that SHG detection is sensitive to δε z   2  whereas direct detection is not.  
         [0041]      FIG. 3  illustrates both SHG and direct detection signals, averaged over x-y, for different values of z slab . The qualitative difference in the traces is striking. The large but gradual increase in &lt;P&gt; 2  as the slab approaches the focal plane indicates that a significant fraction of the transmitted power consists of out-of-focus scattered light. This is expected from the fact that the scattering is mostly forward directed (η z ≈0.9 near the focal plane). As is manifest from  FIG. 3 , the slab displacement must be quite large (|z slab |&gt;400 μm) before scattered light is significantly rejected by the exit pupil. In contrast, out-of-focus scattered light is much more efficiently rejected when using SHG detection because A 00 /A 0S  tends towards zero for relatively small displacements from the focal plane (see  FIG. 2 ). When complete rejection is achieved, only unscattered light produces signal, and &lt;SHG&gt; and &lt;P&gt; 2  are both proportional to e −2NσL , which is z-independent. The apparent plateau in the SHG trace stems from the fact that A 00 /A 0S  and δε z   2  are non-negligible only when the slab spans the focal plane. This plateau clearly identifies the slab boundaries, demonstrating the advantage of improved out-of-focus background rejection with SHG detection.  
         [0042]     Finally, for purposes of illustration, we use our virtual pinhole technique to image an onion slice submerged under a 200_m suspension of 1_m latex beads (number concentration N=0.0048 μm −3 ). The &lt;SHG&gt; image and the corresponding &lt;P&gt; 2  image are shown in  FIG. 4 . The former exhibits both a marked improvement in signal contrast and a suppression of speckle noise presumably caused by scattered background.  
         [0043]     In conclusion, we have demonstrated a new implementation of transmitted light CLSM where an SHG crystal serves as a self-aligned virtual pinhole. Because the SHG signal scales inversely with the area of the incident light distribution, it preferentially reveals unscattered (focused) rather than scattered (diffuse) transmitted power. We emphasize that our technique works well provided an adequate supply of unscattered light survives transmission through the sample. The fact that unscattered power decays exponentially with sample thickness imposes limits on the technique&#39;s applicability. In particular, for thick samples, SHG signal from unscattered light can easily be dominated by SHG from scattered background, despite the suppression of the latter by the virtual pinhole effect. We have empirically observed, with samples comprising 1_m beads, that our technique is effective up to sample thicknesses of roughly 3/Nσ (ie. 3 scattering lengths).  
         [0044]     A notable advantage of our technique lies in its ease of implementation, particularly in combination with standard two-photon excited microscopy, which can be operated simultaneously. Finally, we note that our technique is not limited to signal conversion with an SHG crystal. Alternative techniques involving, for example, 2-photon excited fluorophores or wide-bandgap semiconductors could achieve similar virtual pinhole effects.  
         [0045]     We also describe a simple and robust technique for transmission confocal laser scanning microscopy wherein the detection pinhole is replaced by a thin second-harmonic-generation crystal. The advantage of this technique is that self-aligned confocality is achieved without a need for signal de-scanning. We derive the point-spread function of our instrument, and quantify both signal degradation and background rejection when imaging deep within a turbid slab. As an example, we consider a slab whose index of refraction fluctuations exhibit Gaussian statistics. Our model is corroborated by experiment.  
         [0046]     A pulsed infrared laser beam is focused through a sample and then imaged (re-focused) onto the crystal. A short-pass filter is placed immediately behind the crystal such that only second-harmonic generation (SHG) is detected. Because the SHG power is inversely proportional to the effective area of the laser spot incident on the crystal, the crystal acts as a virtual pinhole, producing a large signal only when the laser spot is tightly focused, similarly to a physical pinhole. The notable advantage of this technique is that virtual confocality is ensured regardless of where the laser spot is focused onto the crystal, meaning that fast beam scanning is allowed without any need for elaborate de-scanning. We call such an instrument an auto-confocal microscope (ACM).  
         [0047]     We presented in the first embodiment a cursory description of an ACM based on quadratic detection and valid for thin samples only. Our goal here is to characterize the imaging properties of such an ACM for both thin and thick samples. In this embodiment, we consider a semi-transparent sample and derive the ACM point-spread function (PSF) for both absorbing and phase-shifting point objects. We qualitatively argue that optical sectioning is obtained only to the extent that scattered background is incoherent. Then we will extend our discussion to thick samples, and explicitly quantify the degree to which ACM rejects scattered background—a fundamental property of confocal microscopy. For simplicity, we consider only non-absorbing media, which we characterize by a (real) refractive index auto-correlation function. Finally, we theoretically evaluate the capacity of an ACM to distinguish a localized object of interest embedded within a turbid slab, assuming the refractive index fluctuations in the slab obey Gaussian statistics.  
         [0048]     The basic layout of our ACM is shown in  FIG. 5 . A laser beam, depicted here as a point source, is focused into a sample by a lens of numerical aperture sin α. A second lens re-focuses this focal spot onto a thin nonlinear crystal. We consider the case of a parfocal geometry wherein the lenses are identical and of unit magnification. Generalizations to non-identical lenses or non-unity magnifications are straightforward and will not be considered here. Our goal in this section is to derive the PSF of an ACM, and discuss its capacity for axial sectioning. We consider a semi-transparent sample and begin by deriving the intensity distribution incident on the image plane (ie. on the nonlinear crystal).  
         [0049]     For ease of notation, we drop all scaling constants throughout this paper. Following the usual notational convention, we write the PSF&#39;s of the lenses as  
               h   ⁡     (       v   1     ,     u   1       )       =     ∫       P   ⁡     (     i   ^     )       ⁢     ⅇ       -   ⅈ     ⁢           ⁢       v   1     ·     i   ^           ⁢     exp   ⁡     (     i   ⁢           ⁢       u   1     ⁡     (         1   4     ⁢           ⁢       sin   2     ⁡     (     α   /   2     )         -       ξ   2     /   2       )         )       ⁢     ⅆ     i   ^                   (   1   )             
 
 where we adopt the axial and radial optical units u 1 =4 k z sin 2 (α/2) and v 1 =k ñ sin α, respectively, and k is the wave-vector in the sample medium. We assume the lenses are ideal and possess no aberrations. That is, the coordinates of the lens pupil functions P({circumflex over (l)}) are normalized such that P(ξ≦1)=1 and P(ξ&gt;1)=0. 
 
         [0050]     To determine the SHG power produced by the crystal, we evaluate the electric field at the image plane, given by 
 
 U ( v )=∫ h ( v   1   ,u   1 ) t ( v   1   ,u   1 ) h ( v−v   1   ,−u   1 ) dv   1   du   1    (2) 
 
 where t(v 1 ,u 1 ) is the 3-dimensional object transmission function, and we neglect multiple scattering since we consider here only semi-transparent samples. 
 
         [0051]     We begin by treating the simplest case of a completely transparent sample that produces no scattering. In this case t(v 1 ,u 1 )=δ(u 1 ) and the field distribution at the image plane becomes 
 
 U   0 ( v )=∫ h ( v   1 ,0) h ( v−v   1 ,0) dv   1   =h ( v ,0)   (3) 
 
         [0052]     Accordingly, the intensity distribution at this plane becomes 
 
 I   0 ( v )=| U   0 ( v )| 2   =∫P ({circumflex over (l)} 1 ) P ({circumflex over (l)} 2 ) e   −iv·({circumflex over (l)}     1     −{circumflex over (l)}     2     )   d{circumflex over (l)}   1 d{circumflex over (l)} 2    (4) 
 
         [0053]     Eq. 4 represents a ballistic light distribution, since it is arises from unscattered transmitted laser light only. Making use of the variable changes {circumflex over (l)} c =({circumflex over (l)} 1 +{circumflex over (l)} 2 )/2 and {circumflex over (l)} d ={circumflex over (l)} 1 −{circumflex over (l)} 2 , we note that I 0 (v) is the Fourier transform of the function: 
 
 H   0 ({circumflex over (l)} d )=∫ P ({circumflex over (l)} c +{circumflex over (l)} d  /2) P ({circumflex over (l)} c −{circumflex over (l)} d /2) d{circumflex over (l)}   c    (5) 
 
         [0054]     Eq. 5 is the diffraction limited optical transfer function (OTF) of a simple lens. This is expected since our parfocal two-lens system is equivalent to a single lens when the sample is transparent. The functions I 0 (v) and H 0 ({circumflex over (l)} d ) will play important roles below.  
         [0055]     To derive the PSF in our microscope configuration, we suppose that our sample now contains a single point perturbation located at the position (v ε ,u ε ). That is, we write 7 : 
 
 t ( v   1   ,u   1   ;v   ε   ,u   ε )=δ( u   1 )−εδ( v   1   −v   ε )δ( u   1   −u   ε )   (6) 
 
 where |ε| is the modulus of the transmission perturbation, assumed small. The real part of ε corresponds to an absorption perturbation whereas the imaginary part corresponds to a phase perturbation. For simplicity, we assume that the sample is scanned in 3-dimensions, with the understanding that formally equivalent results are obtained if the beam is scanned instead of the sample. The perturbed intensity distribution at the image plane is 
 
 I ( v;v   ε   ,u   ε )=| U   0 ( v )−ε U   ε ( v;v   ε   ,u   ε )| 2    (7) 
 
 where 
 
 U   ε ( v;v   ε   ,u   ε )= h ( v   ε   ,u   ε ) h ( v−v   ε   ,−u   ε )   (8) 
 
 and, accordingly, resultant SHG power produced by the crystal is 
 
 SHG ( v   ε   ,u   ε )=∫ I   2 ( v;v   ε   ,u   ε ) dv=S   0 +4 Re[εS   1 ( v   ε   ,u   ε )]+  (9) 
 
 expanded only to the first order perturbation in ε. 
 
         [0056]     The zeroth order ballistic component is defined by 
 
 S   0   =∫|U   0 ( v )| 4   dv=∫I   0   2 ( v ) dv=∫H   0   2 (ξ d ) d{circumflex over (l)}   d    (10) 
 
 where {circumflex over (l)} d  is interpreted as a normalized spatial frequency, and the last equality is an expression of Parseval&#39;s theorem. 
 
         [0057]     The first order term, corresponding to the product of a scattered and three ballistic fields, is defined by 
 
 S   1 ( v   ε   ,u   ε )=∫ I   0 ( v ) U   0 ( v ) U*   ε ( v;v   ε   ,u   ε ) dv    (11) 
 
         [0058]     As is apparent from Eq. (11), the function I 0 (v) plays an identical role here as a pinhole transmission function in standard confocal microscopy—hence the appellation “auto-confocal microscopy”for our technique.  
         [0059]      FIG. 6  depicts various PSF&#39;s obtained for purely absorbing or phase shifting perturbations. We recall that the function I 0 (v) represents the distribution of ballistic light at the crystal plane. As defined by Eq. 4, I 0 (v) is the Airy function (J 1 (v)/v) 2  whose effect, as observed from  FIG. 6 , is essentially identical to that of a standard TCLSM. When using an amplitude perturbation, the contrast S 1 /S 0  of our ACM is found to be the same as that of a standard TCLSM whose pinhole radius is 1.65 optical units. This comparison provides a convenient estimate for the effective pinhole size of our ACM.  
         [0060]     The theoretical results shown in  FIG. 6  may be compared with the experimental results of first embodiment, bearing in mind that the point perturbation in this reference (a latex bead) provoked both a phase-shift and an effective absorption, since the light scattered by the bead was partially clipped by the lens pupil.  
         [0061]     We also note that while a pure phase shifting perturbation does not change the total power incident on the image plane, it can, according to  FIG. 6   c  lead to an increase in power transmitted through a finite (but non-zero) size pinhole.  
         [0062]     It is well known that the main advantage of confocal fluorescence microscopy is its capacity for out-of-focus fluorescence background rejection. In particular, a uniformly fluorescent transverse slice produces a signal that scales as u s   −2 , where u s  is its axial distance from the focal plane. Such a scaling law, which is necessary for optical sectioning, applies even in a transmission geometry because of the incoherent (random phase) nature of fluorescence emission.  
         [0063]     However, there is a fundamental difference between TCLSM&#39;s that are based on fluorescence and on transmission. Whereas a fluorescence microscope exhibits a dark background in the absence of a sample, an ACM, in contrast, exhibits a bright background, stemming from the term S 0  in Eq. 9. This background cannot be easily eliminated. Moreover, the capacity of an ACM for optical sectioning is sample dependent. This problem is readily apparent if one considers simple samples such as a uniformly phase-shifting or absorbing transverse slice. The ACM signals produced by either of these samples is independent of u s  and no optical sectioning is possible (this inability to reject a uniform background is sometimes referred to as the “missing-cone” problem).  
         [0064]     However, samples of interest are rarely so simple. If one considers a transverse slice that instead produces locally random phase-shifts or absorptions (about a mean), the signal produced by an ACM then crucially depends on u s . The transmittance of such samples can be written as  
               t   ⁡     (       𝓋   1     ,       u   1     ;     u   s         )       =       δ   ⁡     (     u   1     )       +       ∑   n     ⁢       ɛ   n     ⁢     δ   ⁡     (       𝓋   1     -     𝓋     ɛ   n         )       ⁢     δ   ⁡     (       u   1     -     u   s       )                     (   12   )             
 
 where infinitesimally small area elements are summed, characterized by complex perturbations ε n  that are randomly distributed in phase. Insertion of Eq. 12 into Eq. 9 leads to a cancellation of the S 1  term, leaving the second order term as a sample dependent response. Such a response exhibits optical sectioning since it scales with u s  in the same way as a fluorescence confocal response. In effect, by imposing random phases to ε n  we have mimicked the incoherence of a fluorescence signal. We note that, while our argument assumes that each perturbation ε n  covers an infinitesimally small area, it remains valid even for finite area perturbations, provide these are small relative to the local laser-beam spot size. Hence, even though the optical sectioning may not be as tightly confined as with a standard fluorescence confocal microscope, it remains nonetheless confined since the laser spot-size expands with increasing |u s |. 
 
         [0065]     In practice, samples of interest are often highly scattering, leading to severe limitations on imaging depth. Our goal in this section is to quantify these limitations by extending our above analysis to thick samples. We consider an intermediate regime often encountered in biological imaging wherein light propagating through a sample is neither wholly ballistic nor wholly diffusive. In particular, we consider scattering that is dominantly forward directed. Such scattering arises from samples that provoke local phase variations that do not significantly deflect the light field but nonetheless highly degrade image quality. We adopt the geometry shown in  FIG. 3 . The sample consists of a slab of thickness L, in which a small object of interest is embedded. As in Section II.B, we suppose the object provokes a localized amplitude or phase perturbation whose signal we wish to evaluate. We derive both signal and background as a function of slab position (or, equivalently, object depth). For simplicity, we assume that the object is situated exactly at the focal point, and that the slab medium is non-absorbing, homogeneous, and isotropic. These assumptions allow us to emphasize the main features of our results, though they are not fundamental to our analysis.  
         [0066]     As is apparent from  FIG. 7   a , the sample may be thought of as two adjacent semi-slabs of thicknesses L A  and L B , situated respectively before and after the object plane. By assumption, back-scattered light is neglected and we assume the light traversing these semi-slabs travels from left to right only. The semi-slabs provoke random phase variations in the light field whose effects can be examined separately: semi-slab A de-focuses the light as it propagates to the object plane, while semi-slab B further de-focuses the light as it continues to propagate to the image plane. Bearing this picture in mind, we develop a formalism based on the alternative equivalent geometry shown in  FIG. 7   b , where we project the phase variations provoked respectively by semi-slabs A and B into the pupil functions of the corresponding lenses. In other words, we mimic the de-focusing effects of the semi-slabs by introducing lens aberrations, and write: 
 
P A,B ({circumflex over (l)})→P({circumflex over (l)})e iδφ     A,B     ({circumflex over (l)})    (13) 
 
         [0067]     The statistics of these aberrations must be correctly defined so as to properly match those of the semi-slabs. We will discuss how to define these statistics. For now, we assume the lens aberrations are characterized by their auto-correlation function, which, by assumption of transverse homogeneity and isotropy, is a function only of the distance between the aberration coordinates. We write, for lens A, 
 
Γ φ   (A) (ξ d )=&lt;δφ A ({circumflex over (l)} 1 )δφ A ({circumflex over (l)} 2 )&gt;  (14) 
 
 where {circumflex over (l)} d ={circumflex over (l)} 1 −{circumflex over (l)} 2  and the brackets correspond to an ensemble average, and we assume Γ φ   (A) (ξ d )→0 for ξ d  sufficiently large. A similar equation applies to lens B. Also, since the phase variations provoked by the semi-slabs are assumed to be uncorrelated, then &lt;δφ A ({circumflex over (l)} 1 )δφ B ({circumflex over (l)} 2 )&gt;=0. 
 
         [0068]     Before deriving the signal produced by an isolated perturbation of interest, we derive the associated background in the absence of any specific perturbation. As previously, we must calculate the field U 0  at the image plane. This time, however, we take into account the phase shifts incurred by the light upon propagation through the entire slab thickness. These are δφ L ({circumflex over (l)})=δφ A ({circumflex over (l)})+δφ B ({circumflex over (l)}). By correspondence with Eqs. 1 and 3, we write 
 
U 0 (v)→∫P({circumflex over (l)} 1 )e −iv{circumflex over (l)}     1   e iδφ     L     ({circumflex over (l)}     1     ) d{circumflex over (l)} 1    (15) 
 
 leading to 
 
 S   0   =∫P ({circumflex over (l)} 1 ) P ({circumflex over (l)} 2 ) P ({circumflex over (l)} 3 ) P ({circumflex over (l)} 4   e   −iv·({circumflex over (l)}     2     +{circumflex over (l)}     3     −{circumflex over (l)}     4     )   K   L   1,2,3,4   d{circumflex over (l)}   1   d{circumflex over (l)}   2   d{circumflex over (l)}   3   d   4   dv    (16) 
 
 where we have defined 
 
 K   L   1,2,3,4 =exp [ i (δφ L ({circumflex over (l)} 1 )−δφ L ({circumflex over (l)} 2 )+δφ L ({circumflex over (l)} 3 )−δφ L ({circumflex over (l)} 4 ))]  (17) 
 
         [0069]     Since we are concerned here with a typical background, we perform an ensemble average of K L   1,2,3,4 . By assumption, the slab is thick enough that δφ L  represents a sum of many independent phase variations, and we write  
               〈     K   L     1   ,   2   ,   3   ,   4       〉     =     exp   ⁡     (       -     1   2       ⁢       ∑     i   ,   j       ⁢         (     -   1     )       i   +   j       ⁢     〈         δϕ   L     ⁡     (       i   ^     i     )       ⁢       δϕ   L     ⁡     (       i   ^     j     )         〉           )               (   18   )             
 
 where we have invoked the Central Limit Theorem and made use of the relation  
         〈     exp   ⁡     (     i   ⁢           ⁢   δϕ     )       〉     =     exp   ⁡     (       -     1   2       ⁢     〈     δϕ   2     〉       )           
 
 valid for Gaussian variables. 
 
         [0070]     An integration of Eq. 16 over the variable v imposes the constraint {circumflex over (l)} 1 −{circumflex over (l)} 2 ={circumflex over (l)} 4 −{circumflex over (l)} 3 , leading to the simplification  
               〈     K   L     1   ,   2   ,   3   ,   4       〉     =         H   L   2     ⁡     (     ξ   d     )       ⁢     exp   ⁡     [       -     1   2       ⁢     ∑   ′       ]                 (   19   )             
 
 where 
 
Σ′=&lt;δφ L ({circumflex over (l)} 1 )δφ L ({circumflex over (l)} 3 )&gt;−δφ L ({circumflex over (l)} 1 )δφ L ({circumflex over (l)} 4 &gt;+&lt;δφ L ({circumflex over (l)} 2 )δφ L ({circumflex over (l)} 4 )&lt;−&gt;δφ L ({circumflex over (l)} 2 )δφ L ({circumflex over (l)} 3 )&lt;  (20) 
 
 and we have introduced the transfer function 
 
 H   L (ξ d )=exp(−Γ φ   (L) (0)+Γ φ   (L) (ξ d ))   (21) 
 
         [0071]     The exponent in Eq. 21 is often referred to as (twice) the structure function of the phase variations {circumflex over (l)}. The physical meaning of H L  will be elaborated on below. We note here that if the slab is transparent (or nonexistent), then H L (ξ d )=1 for all ξ d . If, instead, the slab is thick enough to provoke significant phase variations, then H L (ξ d ) rapidly decays from unity at ξ d =0 to a small baseline value exp(&lt;Γ φ   (L) (0) (see  FIG. 4 ). We define below what we mean by “significant” and assume for now that H L  is sufficiently peaked around the origin that H L   2  takes on non-negligible values in Eq. 19 only when ξ d ≈0. As a result, the main contribution in the integration in Eq. 16 comes from the region where &lt;K L   1,2,3,4 &gt;≈H L   2 (ξ d ). We then obtain 
 
&lt; S   0   &gt;≈∫P ({circumflex over (l)} c +{circumflex over (l)} d /2) P ({circumflex over (l)} c +{circumflex over (l)} d /2)  P ({circumflex over (l)}′ c −{circumflex over (l)} 2 /2) P ({circumflex over (l)}′ c −{circumflex over (l)} d /2) H   L   2 (ξ d ) d{circumflex over (l)}   c   d{circumflex over (l)}   c   ′d{circumflex over (l)}   d    (22) 
 
 which, with Eq. 5, simplifies to, 
 
&lt; S   0   &gt;≈∫H   0   2 (ξ d ) H   L   2 (ξ d ) d{circumflex over (l)}   d    (23) 
 
         [0072]     &lt;S 0 &gt; is the average background SHG power obtained when only the slab is taken into account and nothing more (ie. no object of interest lies at the focal center). A comparison of Eq. 23 with Eq. 10 suggests that H L (ξ d ) can be interpreted as a filter function similar to H 0 (ξ d ). By limiting the extent of the spatial frequencies that are transferred to the image plane, H L (ξ d ) provokes a blurring of the focal spot incident on the SHG crystal. Hence, though the presence of the slab does not alter the total power incident on the crystal, it does lead to a reduction in the resultant SHG the crystal produces. The intrinsic sensitivity of nonlinear detection to de-focusing is the basis of ACM background rejection.  
         [0073]     We now derive the signal produced by point object located at the focal center. We use the same formalism developed above for deriving background, but this time we treat the semi-slabs individually. Referring to Eq. 8, and explicitly identifying the respective phase aberrations in lenses A and B, we write, 
 
U ε (v)→∫P({circumflex over (l)} 1 )P({circumflex over (l)} 2 )e −v·{circumflex over (l)}     2   e i(δφ     A     ({circumflex over (l)}     1     )+δφ     B     ({circumflex over (l)}     2     )) d{circumflex over (l)} 1 d{circumflex over (l)} 2    (24) 
 
         [0074]     We will restrict our analysis here to the first order perturbation for both absorption and phase contrasts. This first order signal (Eq. 11) becomes 
 
 S   1   =∫P ({circumflex over (l)} 1 ) P ({circumflex over (l)} 2 ) P ({circumflex over (l)} 3 ) P ({circumflex over (l)} 4 ) P ({circumflex over (l)} 6 ) e   −iv·({circumflex over (l)}     1     −{circumflex over (l)}     2     +{circumflex over (l)}     3     −{circumflex over (l)}     6     )    K   A   1,2,3,4,   K   B   1,2,3,6   d{circumflex over (l)}   1   d{circumflex over (l)}   2   d{circumflex over (l)}   3   d{circumflex over (l)}   4   d{circumflex over (l)}   6   dv    (25) 
 
 where we have used definitions for K A  and K B  similar to Eq. 17 and adjusted our indices in accord with Eq. 18. An integration over the variable v leads to the constraint {circumflex over (l)} d ={circumflex over (l)} 1 −{circumflex over (l)} 2 ={circumflex over (l)} 6 −{circumflex over (l)} 3  and, following the same reasoning as in the previous section, we obtain 
 
&lt; K   B   1,2,3,6   &gt;≈H   B   2 (ξ d )   (26) 
 
&lt; K   A   1,2,3,4   &gt;≈H   A (ξ d ) H   A (ξ′ d )   (27) 
 
 where we have defined {circumflex over (l)}′ d ={circumflex over (l)} 3 −{circumflex over (l)} 4  and have assumed that H B   2 (ξ d ) is non-negligible only for small ξ d , as before, leading to the restriction {circumflex over (l)} 6 ≈{circumflex over (l)} 3 . The signal produced by a localized amplitude perturbation is then given by 
 
&lt; S   1   &gt;≈∫H   0 (ξ′ d ) H   A (ξ′ d ) d{circumflex over (l)}′   d   ∫H   0 (ξ d ) H   B   2 (ξ d ) H   A (ξ d ) d{circumflex over (l)}   d    (28) 
 
         [0075]     We note that Eq. 28 resembles Eq. 23 except that a component of the light transmitted through semi-slab A prior to its interaction with the object has been isolated (first integral). We also remind the reader that S 1  reveals a phase gradient rather than a phase exactly at the focal center (see  FIG. 6   c ).  
         [0076]     Eqs. 23 and 28 are the main results of this section, and represent formal expressions for the background and highest order signal obtained when using a quadratic detection ACM to image inside a thick slab.  
         [0077]     So far, we have made no assumptions on the detailed nature of the phase fluctuations introduced by the slab. We consider here the specific example where these are produced by refractive index fluctuations that obey locally Gaussian statistics. Such statistics are routinely used to describe scattering media, and are particularly convenient because of their tractability. To this end, we define a transverse autocorrelation function for the refractive index fluctuations, 
 
&lt;δ n ( {tilde over (n)}   1 )δ n ( ñ   2 )&gt;=&lt;δ n   2 &gt; exp(−ρ d   2   /l   n   2 )   (29) 
 
 where we have reverted to the lab-frame coordinate system (ñ,z) relative to the focal center, and l n  is a characteristic fluctuation scale, assumed to be the same in all three dimensions. If light propagates an axial distance δz&lt;&lt;l n , it incurs a phase shift kδz. On the other hand, for longer axial distances δz&gt;&gt;l n  then the phase shift is no longer proportional to the propagation distance but instead performs a random walk with step size≈kl n . In this latter case the variance of the phase fluctuations, as opposed to their amplitude, scales linearly with axial propagation distance, and we write 
 
Γ φ   (δz) (ρ d )≈δ zl   n   k   2   &gt;δn ( ñ   1 ) δn ( ñ   2 ) &gt;  (30) 
 
 where ñ d =ñ 1 −ñ 2 , and δz is assumed to be small enough that we may neglect beam convergence or divergence. 
 
         [0078]     As described above, we use the technique of projecting the slab fluctuations into the lens pupils, which requires the coordinate transformation ρ d →ξ d z sin α. Referring to Eq. 21, we obtain then, 
 
 H   δz (ξ d )≈exp(−δzσ φ   2 (1−γ φ   (δz) (ξ d   z  sin α)))   (31) 
 
 where σ φ   2 ≈k 2 l n &lt;δn 2 &gt; is the variance of the phase fluctuations per unit propagation distance, and we define γ 100    (δz) (ρ d )=Γ φ   (δz) (ρ d )/Γ φ   (δz) (0). We note that γ φ   (δz) (ρ d ) is always 
 
 equal to one at the origin, but becomes more and more narrowly peaked as the propagation distance through the slab increases. 
 
         [0079]     To derive the filter function through a thick slab, not just a thin slice, we must take beam convergence or divergence into account. Since the filter functions for sequential slices of thickness δz are assumed to operate independently, we make the approximation 
 
 H   L (ξ d )≈Π L   H   δz (ξ d )   (32) 
 
         [0080]     This last step represents one of the main advantages of our having projected the phase fluctuations from the slab (spatial coordinates) to the lens pupils (frequency coordinates) where the filter functions operate multiplicatively.  
         [0081]     Expression 32 is a product over the entire slab thickness, and can be evaluated by integrating the exponent in Eq. 31. We obtain the approximate expression 
 
 H   L (ξ d )≈exp(− Lσ   φ   2 )+(1−exp(− Lσ   φ   2 ) exp(−ξ d   2 σ φ   2   V/l   n   2 )   (33) 
 
 where we have defined V=∫ L (z sin α) 2 δz, which roughly corresponds to the volume of the laser beam inside the slab (shaded region in  FIG. 3   a ). As described in section III.A and is explicit in Eq. 33, H L (ξ d ) consists of a spatial-frequency-independent baseline (first term) onto which a narrow peak around the origin (second term) is superposed (see  FIG. 8 ). The physical meaning of these terms is as follows: 
 
         [0082]     H L (ξ d ) represents the effect of the slab on the transmitted beam. This effect is two-fold. The baseline term in Eq. 33 is an expression of Lambert&#39;s law and describes the frequency-independent attenuation of the ballistic (non-scattered) light transmitted through the slab. With this interpretation, the scattering mean-free-path (MFP) of the slab is defined as l s =σ φ   −2 . The peak term in Eq. 33 represents the effect of H L (ξ d ) on the rest of the light transmitted through the slab that has been scattered. Whereas very low spatial frequencies are efficiently transmitted, frequencies higher than a cut-off ξ 3 dB ≈l n {square root}{square root over (l s /V)} are severely attenuated. We remind the reader that a diffraction-limited focus requires a transmission of frequencies up to ξ d ≈1. Hence, inasmuch as ξ 3 dB &lt;&lt;1 (we will quantify this below), the second term in H L (ξ d ) leads to a significant blurring of the non-ballistic light at the image plane.  
         [0083]     We now directly evaluate the background produced by the SHG crystal. For convenience, we make two approximations. First, even though H 0 (ξ d ), as defined by Eq. 5, can be expressed analytically, we adopt the much simpler Gaussian beam approximation H 0 (ξ d )≈exp(−ξ d   2 ), which is valid in the paraxial limit. Second, we relate l n  to the more experimentally accessible transport MFP, defined by l* s =l s /(1−&lt; cos θ s &gt;), where θ s  is the deflection angle occasioned by a single scattering event. For Gaussian refractive index fluctuations (Eq. 29), these are approximately related by k 2 l n   2 ≈l* s /l s . As an example, l n  is on the order of a micron for most biological tissues of interest, meaning that the scattering is highly forward directed at optical wavelengths and l* s  is typically 10 to 20 times longer than l s .  
         [0084]     Using Eq. 33 and performing the integral in Eq. 23, we obtain 
 
&lt; S   0   &gt;≈SHG   0 { exp(−2 L/l   s )+(1−exp(−2 L/l   s )) R ( l*   s   ,V )}  (34) 
 
 where SHG 0  corresponds to the SHG power obtained if there were no slab (L=0). As discussed above, the effect of the slab is to convert non-scattered ballistic light into scattered light. The thicker the slab, the more this conversion is complete, and the first and second terms in Eq. 34 correspond to these ballistic and non-ballistic components respectively. However the non-ballistic component is significantly rejected here by the factor  
               R   ⁡     (       l   s   *     ,   V     )       =     (       l   s   *         l   s   *     +       k   2     ⁢   V         )             (   35   )             
 
 This rejection factor is a fundamental consequence of the fact that de-focused non-ballistic light is ineffective in producing SHG. The greater the de-focusing, the greater the rejection, as indicated by the relation R(l* s ,V)≈ξ 3 dB   2 . Moreover, the rejection depends only on the intrinsic slab parameter l* s , and on extrinsic parameters such as slab thickness and position along the optical axis, both of which govern the interaction volume through the geometric relation  
             V   =         V   A     +     V   B       ≈           sin   2     ⁢   α     3     ⁢     (       L   A   3     +     L   B   3       )                 (   36   )             
 
         [0085]     An illustration of R for different V&#39;s is shown in  FIG. 9 . In this experimental example, the slab is thick enough that the ballistic component can be neglected, meaning &lt;S 0 &lt;≈R(l* s ,V). We emphasize that even when the interaction volume is at a minimum here (slab is centered on focal plane), the rejection factor still remains considerably smaller than one, indicating that non-ballistic light is highly de-focused and justifying a posteriori the assumptions that led to Eq. 22. The theoretical fit shown in  FIG. 9  contains no free parameters and is remarkably accurate despite the simplicity of Eq. 35. We note that, for this example and those presented henceforth, V is always large enough that we may use the approximation R(l* s ,V)≈l* s /k 2 V.  
         [0086]     To evaluate the capacity of our ACM to perform deep imaging in a scattering slab, we consider, as previously, the signal produced by a point perturbation of interest located at the focal center. The depth L A  of this perturbation relative to the slab surface is governed by the slab position, which in turn governs L B , V A , V B , and V (only L remains unchanged). Approximating the filter functions in Eq. 28, as was done above to obtain Eq. 33, we arrive at  
               〈     S   1     〉     ∝       {       exp   ⁡     (       -     L   A       /     l   s       )       +       (     1   -     exp   ⁡     (       -     L   A       /     l   s       )         )     ⁢     R   ⁡     (       l   s   *     ,     V   A       )           }     ×     {       exp   ⁡     (       -     (     L   +     L   B       )       /     l   s       )       +       (     1   -     exp   ⁡     (       -     (     L   +     L   B       )       /     l   s       )         )     ⁢     R   ⁡     (       l   s   *     ,     V   +     V   B         )           }               (   37   )             
 
         [0087]     The leftmost bracketed terms in Eq. 37 represents the laser intensity incident exactly at the point object, consisting of ballistic and non-ballistic components. The latter component is diminished by the factor R(l s ,V A ) because of spreading of the non-ballistic light.  
         [0088]     Using the same apparatus as described in the first embodiment we experimentally corroborate the validity of these results with test slabs consisting of 1 μm latex beads embedded in scattering media (themselves composed of latex beads, some of which are fluorescent, in agarose gels). The parameters l s  and l* s  can be prescribed for each slab based on the sizes and concentrations of the beads. Moreover, the parameter l s  can easily be verified by monitoring the average two-photon excited fluorescence (TPEF) signal produced by the fluorescent beads, which is known to decay as exp(−2L A /l s ) to moderate depths.  
         [0089]     Two regimes may be distinguished, based on the relative contributions of ballistic and non-ballistic components in the average SHG signal (Eq. 34). If the ballistic component is dominant (first term in Eq. 34), then &lt;S 1 &gt; is essentially independent of L A  or V A , meaning that the signal produced by a point object of interest, whether absorbing or phase-shifting, remains the same at all depths throughout the slab. This regime is illustrated in  FIG. 10   b.    
         [0090]     If, instead, the non-ballistic component is dominant (second term in Eq. 34), then the amount of ballistic light incident of the SHG crystal is negligible. This should not be confused, however, with the amount of ballistic light incident on the point object itself, which can be much greater and lead to contrast. This second regime is illustrated in  FIG. 10   c . Hence, though only non-ballistic light produces signal in this second regime, high-resolution images can nonetheless be obtained. The signal here decays with object depth. From Eq. 37, we infer the rough scaling law for moderate depths &lt;S 1 &gt;∝ exp(−L A /ζl s ), where the slab MFP has been effectively lengthened by the factor ζ≈1(1−3l s /L), valid for L&gt;l s . This apparent lengthening of the MFP stems directly from the effectiveness of nonlinear detection in rejecting non-ballistic background. In summary, the main advantages of ACM are that it allows fast beam scanning, provides effective background rejection, and can be readily combined with TPEF microscopy. The performance of an ACM is essentially the same as that of a standard TCLSM. The capacity of an ACM for depth penetration depends on the net amount of ballistic light that traverses the slab. If the slab is thin, then both background and signal are independent of depth. If the slab is thick, then background scales roughly inversely with the light-slab interaction volume while the signal decays moderately with depth, in accord with the simple model presented above.