Interferometric superlocalization of two incoherent optical point sources
Abstract
A novel interferometric method – SLIVER (Super Localization by Image inVERsion interferometry) – is proposed for estimating the separation of two incoherent point sources with a mean squared error that does not deteriorate as the sources are brought closer. The essential component of the interferometer is an image inversion device that inverts the field in the transverse plane about the optical axis, assumed to pass through the centroid of the sources. The performance of the device is analyzed using the CramérRao bound applied to the statistics of spatiallyunresolved photon counting using photon numberresolving and onoff detectors. The analysis is supported by MonteCarlo simulations of the maximum likelihood estimator for the source separation, demonstrating the superlocalization effect for separations well below that set by the Rayleigh criterion. Simulations indicating the robustness of SLIVER to mismatch between the optical axis and the centroid are also presented. The results are valid for any imaging system with a circularly symmetric pointspread function.
Department of Electrical and Computer Engineering,
National University of Singapore, 4 Engineering Drive 3,
Singapore 117583
Department of Physics,
National University of Singapore, 2 Science Drive 3, Singapore 117551
elernair@nus.edu.sg
(110.2970) Image detection systems; (350.5730) Resolution; (110.0180) Microscopy; (030.5290) Photon statistics; (180.2520) Fluorescence Microscopy; (180.3170) Interference Microscopy; (180.1790) Confocal Microscopy; (180.5810) Scanning Microscopy.
References
 [1] Lord Rayleigh, “XXXI. Investigations in optics, with special reference to the spectroscope,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 8, 261–274 (1879).
 [2] M. Born and E. Wolf, Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light (Cambridge University, 1999).
 [3] E. Bettens, D. Van Dyck, A. Den Dekker, J. Sijbers, and A. Van den Bos, “Modelbased twoobject resolution from observations having counting statistics,” Ultramicroscopy 77, 37–48 (1999).
 [4] S. Ram, E. S. Ward, and R. J. Ober, “Beyond Rayleigh’s criterion: a resolution measure with application to singlemolecule microscopy,” Proceedings of the National Academy of Sciences of the United States of America 103, 4457–4462 (2006).
 [5] S. V. Aert, A. den Dekker, D. V. Dyck, and A. van den Bos, “Highresolution electron microscopy and electron tomography: resolution versus precision,” Journal of Structural Biology 138, 21 – 33 (2002).
 [6] H. L. Van Trees, Detection, Estimation, and Modulation Theory : Part I (WileyInterscience 1st Ed, 2001), 1st ed.
 [7] M. Tsang, R. Nair, and X.M. Lu, “Quantum theory of superresolution for two incoherent optical point sources,” (2015). http://arxiv.org/abs/1511.00552.
 [8] S. Weisenburger and V. Sandoghdar, “Light microscopy: an ongoing contemporary revolution,” Contemporary Physics 56, 123–143 (2015).
 [9] J. Chao, S. Ram, E. S. Ward, and R. J. Ober, “Ultrahigh accuracy imaging modality for superlocalization microscopy,” Nat Meth 10, 335–338 (2013).
 [10] M. Tsang, “Quantum limits to optical pointsource localization,” Optica 2, 646–653 (2015).
 [11] C. W. Helstrom, Quantum Detection and Estimation Theory (Academic Press, 1976).
 [12] A. S. Holevo, Probabilistic and Statistical Aspects of Quantum Theory (Edizioni Della Normale, Pisa, Italy, 2011).
 [13] M. G. Paris, “Quantum estimation for quantum technology,” International Journal of Quantum Information 7, 125–137 (2009).
 [14] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University, 1995).
 [15] J. H. Shapiro, “The quantum theory of optical communications,” IEEE Journal of Selected Topics in Quantum Electronics 15, 1547 –1569 (2009).
 [16] J. W. Goodman, Statistical Optics (John Wiley & Sons, 1985).
 [17] N. Sandeau and H. Giovannini, “Increasing the lateral resolution of 4pi fluorescence microscopes,” J. Opt. Soc. Am. A 23, 1089–1095 (2006).
 [18] K. Wicker and R. Heintzmann, “Interferometric resolution improvement for confocal microscopes,” Optics Express 15, 12206–12216 (2007).
 [19] K. Wicker, S. Sindbert, and R. Heintzmann, “Characterisation of a resolution enhancing image inversion interferometer,” Optics Express 17, 15491–15501 (2009).
 [20] D. Weigel, R. Foerster, H. Babovsky, A. Kiessling, and R. Kowarschik, “Enhanced resolution of microscopic objects by image inversion interferometry,” Optics Express 19, 26451–26462 (2011).
 [21] D. Weigel, H. Babovsky, A. Kiessling, and R. Kowarschik, “Investigation of the resolution ability of an image inversion interferometer,” Optics Communications 284, 2273–2277 (2011).
 [22] D. Weigel, T. Elsmann, H. Babovsky, A. Kiessling, and R. Kowarschik, “Combination of the resolution enhancing image inversion microscopy with digital holography,” Optics Communications 291, 110 – 115 (2013).
 [23] D. Weigel, H. Babovsky, A. Kiessling, and R. Kowarschik, “Imaging properties of different types of microscopes in combination with an image inversion interferometer,” Optics Communications 332, 301–310 (2014).
 [24] D. Weigel, H. Babovsky, A. Kiessling, and R. Kowarschik, “Widefield microscopy with infinite depth of field and enhanced lateral resolution based on an image inverting interferometer,” Optics Communications 342, 102 – 108 (2015).
 [25] J. W. Goodman, Introduction to Fourier Optics (Roberts and Company Publishers, 2005).
 [26] S. Kay and Y. C. Eldar, “Rethinking biased estimation,” IEEE Signal Processing Magazine 25, 133–136 (2008).
 [27] R. S. Kennedy, “A nearoptimum receiver for the binary coherent state channel,” Tech. Rep. 108, MIT RLE Quarterly Progress Report (1973).
 [28] S. J. Dolinar, “An optimum receiver for the binary coherent state channel,” Tech. Rep. 111, MIT RLE Quarterly Progress Report (1973).
 [29] R. Nair, S. Guha, and S.H. Tan, “Realizable receivers for discriminating coherent and multicopy quantum states near the quantum limit,” Physical Review A 89, 032318 (2014).
 [30] U. Kubitscheck (Ed.), Fluorescence Microscopy: from Principles to Biological Applications (John Wiley & Sons, 2013).
Appendix A Introduction
Rayleigh’s criterion for resolution of two incoherent point sources, which asserts that a minimum separation between the sources equal to the diffractionlimited spot size is necessary for them to be resolvable, has perhaps been the most influential resolution criterion in the history of optics despite its heuristic character [1, 2]. Based as it is on the eye as a detection instrument, Rayleigh’s criterion neglects the possibility of better resolution using light of greater intensity or using a longer observation period. Pioneering studies on resolution limits as a function both of the source separation and the mean number of observed photons were made by Bettens et al. in [3] and more recently by Ram et al. in [4] (See also [5] in the context of electron microscopy). In [3], it was shown that for a given mean photon number, ideal continuum imageplane photon counting can locate the centroid, i.e., the midpoint, of the two radiating sources with a finite precision (depending on the size of the diffractionlimited spot in the image plane and on the number of photons detected) no matter how small their separation is. Based on the CramérRao (CR) lower bound of estimation theory [6], it was also shown in [3, 4] that any (unbiased) estimate of the separation between the sources based on imageplane photon counting must suffer a divergent mean squared error for a given mean photon number as the separation tends to zero. This phenomenon was dubbed Rayleigh’s curse in [7] as it suggests a fundamental limitation in resolving incoherent point sources even when the role of the average detected photon number is taken into account. In the past few decades, several farfield superresolution techniques that circumvent the Rayleigh limit have revolutionized microscopy (see [8] for a review), but these rely on placing restrictions on the emissions from the sources and do not fundamentally challenge the criterion (or Rayleigh’s curse) for independently and incoherently radiating sources. Other developments, e.g.,[9], use clever techniques to approach the performance of ideal continuum imageplane photon counting using imperfect detectors.
In [7], following preliminary work in [10], the problem of resolving two incoherent point sources was approached anew from the perspective of quantum estimation theory using the quantum CramérRao (QCR) bound [11, 12, 13]. This bound provides a fundamental limit to the accuracy of estimating the source separation optimized over all possible measurement techniques allowed by quantum mechanics. Under a weaksource assumption, it was found in [7] that the QCR lower bound on the minimum mean squared error (MSE) of estimating the separation of two point sources is independent of that separation. Further, a linear opticsbased measurement – spatialmode demultiplexing (SPADE) – was proposed and shown in principle to approach the quantum bound for all values of the separation. These results are in stark contrast to the performance of imageplane photon counting mentioned above, as the divergent behavior of the minimum MSE with decreasing separation – Rayleigh’s curse – is conspicuously absent.
In this paper, we revisit the resolution problem from the point of view of the semiclassical theory of photodetection [14, 15, 16]. For sources emitting thermal radiation, for linear propagation of light through the imaging system, and for linearoptics processing followed by photon counting, the semiclassical theory of photodetection may be used to obtain the same results as a fully quantum analysis would give [14, 15]. This is because the thermal state has a positive representation, and its counting statistics may be obtained by statistical averaging over those of coherent states via Mandel’s rule [14]. The SLIVER scheme of this paper as well as the schemes of [7] satisfy the above necessary conditions for the applicability of the semiclassical theory. We emphasize that, for the state of radiation and measurements considered here and in [7], “semiclassical” does not imply “approximate” – the fully quantum and semiclassical treatments are in quantitative agreement.
Besides its greater familiarity to the optics community, a semiclassical analysis of the resolution problem offers several advantages. Firstly, the analysis can be carried out for arbitrary source strengths. Indeed, the work in [3, 4] assumes the sources are weak enough so that the counting statistics are nearly Poissonian. Similarly, the quantum analysis in [7] assumed, for mathematical tractability, that the sources are weak enough so that the state of light could be assumed to be confined to the zero and onephoton subspaces. Thus, a semiclassical treatment will incorporate the effects of multiphoton events neglected in [3, 4] and in the quantum analysis of [7]. Such an analysis would also reveal whether the superresolution effect predicted in [7] for weak sources persists as the source strength is increased. A second advantage of the semiclassical approach is that intuitions from semiclassical optics can be harnessed to both understand the reasons for superresolution as well as to design new measurement schemes which can also be analyzed semiclassically provided that they only involve linear optics and photodetection, possibly along with the use of other classical sources of light.
Accordingly, our contributions in this paper are threefold:

We set up the problem of resolving two incoherently radiating point sources in the framework of semiclassical photodetection theory. The sources can be of arbitrary strength and imaging systems with twodimensional circularlysymmetric pointspread functions are studied. The thermal source model is applicable to a wide range of physical scenarios ranging from optical astronomy to fluorescence microscopy.

We propose a new interferometric scheme for estimating source separation that we call SLIVER (Super Localization by Image inVERsion interferometry) and that yields finite resolution for arbitrarily small values of the source separation and for arbitrary source strengths. Devices employing image inversion have been proposed, studied and demonstrated previously in the microscopy literature for various related applications [17, 18, 19, 20, 21, 22, 23, 24], but the fundamental limits on the capability of such devices for improving lateral resolution seem not to have been explored before. Through an analysis of the photodetection statistics of SLIVER in the semiclassical framework, we show both using the CramérRao bound and through explicit simulation that our method manifestly alleviates Rayleigh’s curse.

We offer an explanation for the superlocalization effect of SLIVER at small values of separation in the language of estimation theory applied to the photodetection statistics for thermal light.
Appendix B Source and system model
Consider the depiction in Fig. 1 of the focused image of two incoherent point sources in the image plane of an imaging system coordinatized by . We assume that the imaging system is spatially invariant and coordinates have been rescaled so that the image is of unit magnification [25]. We also assume that the (possibly complexvalued) normalized amplitude pointspread function (PSF) satisfying
(1) 
is circularly symmetric, so that depends only on the magnitude for all . It is straightforward to make spatially invariant and circularly symmetric in an imaging system using two lenses [25]. The prototypical example of such a PSF is of course the Airy disk pattern resulting from a circular pupil. The images of the sources 1 and 2 in the image plane are assumed to be centered at respectively, giving rise in the image plane to a combined field with complex amplitude
(2) 
in normalized units of . We will assume throughout this paper, except in Section 5.3, that the centroid of the sources, i.e., their midpoint, has already been located, perhaps by imageplane photon counting [4]. It is taken to lie at the origin of the image plane. In Sec. 5.3, we consider the effect of a small error in centroid localization. Our initial results suggest that our measurement scheme is robust to such perturbations.
The pair of dimensionless complex numbers are the source amplitudes. The thermal and mutually incoherent nature of the sources dictates that the source amplitudes are circularcomplex Gaussian random variables having the first and second moments [16]:
(3)  
(4)  
(5)  
(6) 
for ranging over the two sources. According to the above relations, the real and imaginary parts of the comprise four statistically independent zeromean Gaussians of variance (for the components of ) or (for the components of ). Eqs. (3)(6) are as appropriate for completely incoherent sources with respective strengths (i.e., mean photon number) and , which can assume any values.
Given that the centroid has been located, the parameter of interest is the separation between the two sources. The angle between the line joining the sources and the axis is another unknown parameter, but the measurement method proposed here works regardless of the value of that angle, which need not be known beforehand in practice. We will thus focus just on the separation and apply the tools of singleparameter estimation theory in our analysis.
Appendix C The SLIVER measurement and its statistics
The SLIVER measurement, illustrated schematically in Fig. 2 and described in detail below, essentially involves the separation of the imageplane field of Eq. (2) into its symmetric and antisymmetric parts with respect to inversion in the image plane about the optical axis. By definition, the symmetric and antisymmetric parts of are given respectively by
(7) 
In writing the righthand sides of Eqs. (7), we have used the fact that is circularly symmetric and have defined the sum and difference
(8)  
(9) 
of the amplitudes of the two sources. The relations
(10) 
hold for these new random variables, which are thus also identically distributed circularcomplex Gaussians and statistically independent of each other when .
The field components of Eq. (7) may be obtained in spatially separated planes using an image inversion interferometer (III) that works by splitting the input field using a 5050 beamsplitter, spatially inverting the field from one output about the optical axis and recombining the two beams at a second 5050 beamsplitter. The two beams output from the second beamsplitter then have the field patterns of Eq. (7). The entire setup is thus essentially a balanced MachZehnder interferometer with an inversion device in one arm (Fig. 2). The critical element in such setups is the inversion device – a device that implements the transformation . Inversion devices have been proposed and demonstrated for various applications in the fluorescence microscopy community and a considerable literature on image inversion microscopy exists [17, 18, 19, 20, 21, 22, 23, 24]. However, to the best of our knowledge, a fundamental statistical analysis of the kind we make here of the capabilities of an III for lateral resolution has not been done (See the discussion in Section 6). The spatial inversion itself can be accomplished in a variety of ways, e.g., using a arrangement with lenses [18], with a 3D setup using plane mirrors [19], or using two concave mirrors [20]. We will be concerned here only with the outputs (7) of such a device rather than its detailed implementation, which may vary according to the application.
Suppose now that the output beams from the III are directed to different photodetectors. Conditioned on the values of and , the intensity patterns and of the fields on their respective detector planes are
(11)  
(12) 
where
(13) 
is an interference term. Conditioned on the values of and , semiclassical detection theory dictates that the photocounts on the two detector planes are independent inhomogeneous Poisson processes with the rate functions and respectively [14, 16, 15]. It follows that the spatiallyunresolved integrated photocounts and at each detector are Poisson random variables with the means
(14)  
(15) 
where
(16)  
(17)  
(18) 
The circular symmetry of has been used to get Eq. (17), and we have retained only the separation as argument. According to Eqs. (10), both and are exponentially distributed with mean , so the photocounts and integrated over the photodetector surfaces are distributed according to BoseEinstein distributions [14, 16, 15] with means
(19)  
(20) 
for . If , the photocounts are also statistically independent because and are independent in that case.
The SLIVER method consists of optimal statistical processing of the sequence of observed photocount measurements of and (or) over independent shots using the same sample to obtain a good estimate of the separation – see Fig. 2 and Secs. 4 and 5. As we show below, the photocount from the antisymmetric component of the input field is much more informative in this regard.
Appendix D CramérRao bounds for the separation
d.1 Numberresolved photon counting
Given the BoseEinstein distribution
(21) 
of numberresolved photon counting of the inversionsymmetric field, we can calculate the fundamental CramérRao (CR) lower bound of estimation theory [6] on the mean squared error (MSE) of any unbiased estimator of that is a function of the photocount . An unbiased estimator is one that satisfies for all in some interval of interest. The statistical expectations here are taken over the distribution of . The CR bound is itself the reciprocal of the Fisher information on the parameter defined as
(22) 
The superscript indicates that the Fisher information pertains to photon counting of the inversionsymmetric field. This quantity turns out to equal
(23)  
(24)  
(25)  
(26) 
Here we have defined
(27) 
In a similar fashion, we find the Fisher information of photon counting of the inversionantisymmetric field to be
(28) 
The CR lower bounds for these individual measurements then respectively read
(29)  
(30) 
For a fixed , both bounds contain a term that displays the familiar shotnoise scaling with respect to the combined source strength . Of greater interest here, however, is the behavior of the bounds as the separation reduces to zero for fixed source strengths. Since and as , we have
(31) 
so that any unbiased estimator based on bulk direct detection of the inversionsymmetric field suffers a divergent mean squared error for separations . Thus, Rayleigh’s curse that plagues spatiallyresolved imageplane photon counting remains in effect here [3, 4]. However, note that
(32) 
Here we have used
(33)  
(34)  
(35) 
where is the squared spectral width of the PSF. Remarkably therefore, direct detection of the inversionantisymmetric field does not suffer Rayleigh’s curse and a limiting MSE
(36) 
is potentially achievable for small , where is the number of independent measurements, and
(37) 
is the the mean number of photons emitted over the shots. The limiting Fisher information value is similar to the quantum limit obtained in [7] and achieved by the SPADE and binary SPADE schemes therein. An advantage of SLIVER over those schemes is that there is no need to couple the imageplane field into a waveguide(s) with mode profile(s) tailored to the PSF.
If , and are statistically independent so that the two Fisher information terms of eqs. (26) and (28) may be added to give the total Fisher information. The reciprocal of this quantity is then the CR bound for any unbiased estimate of based on both photocounts. In the remainder of the paper, however, we will only consider measurements in the antisymmetric output port.
d.2 Onoff detection of the inversionantisymmetric field
Consider bucket direct detection at the antisymmetric port with an onoff detector, i.e., a detector that only distinguishes between no photon and one or more photons. The noclick probability and click probability are respectively
(38)  
(39) 
The Fisher information for this binary measurement can be shown to be
(40) 
which is expected as onoff detection is a coarser measurement than photon counting. However, it is interesting that onoff detection also evades Rayleigh’s curse for small and the Fisher information attains the optimal value at . Note also that the Fisher information of Eq. (40) is negligibly different from that of Eq. (28) for . These behaviors are expected as multiphoton arrivals in the imageplane are rare in the regime while multiphoton detections at the antisymmetric output port are rare in the regime, i.e., the regime of . Equations (28) and (40) indicate that the two antisymmetric detection methods should exhibit superresolution for arbitrary values of source strength beyond the regime.
We illustrate the above results for the circular Gaussian PSF
(41) 
The PSFdependent quantities appearing in the CR bound then become:
(42)  
(43)  
(44) 
The Fisher information quantities of Eqs. (26), (28), and (40) pertaining to photon counting at the symmetric and antisymmetric ports, and onoff detection at the antisymmetric port of the interferometer are plotted as a function of the separation in Figures 34. The source strength is photons and photons in Figures 3 and 4 respectively. The plots are normalized to the maximum value of attained at . We see that the information obtainable from the antisymmetric detection methods is maximum at and decreases thereafter. The information from detection at the symmetric port is zero for and reaches a peak around and decreases thereafter.
Appendix E MaximumLikelihood (ML) estimates of : MonteCarlo analysis
To further support the predictions from the CR bounds, we present MonteCarlo simulations of the MSE for onoff detection and photon counting at the antisymmetric output port of the interferometer. The circular Gaussian PSF of Eq. (41) is assumed and the simulation results are independent of the halfwidth .
e.1 Onoff detection in the antisymmetric output port
Consider direct detection of the output in detection windows using an onoff detector. The ensuing measurement record is a bit string , with if the detector did not fire in the th detection window and if it did. The maximum likelihood (ML) estimator [6] for is then:
(45) 
Here, is the total number of clicks observed in the detector during the measurement, which is a sufficient statistic for generating the estimate. The second case above is necessary because the equation for the maximum likelihood estimate of has no solution if , in which case the estimate is set to an arbitrary value. As increases, the probability of such “large deviations” goes to zero.
Figure 5 shows the simulated MSE of the ML estimate for a source strength of photons. The MSE is shown scaled relative to the minimum value of the CR bound for values of , 200 and . As the separation increases beyond the halfwidth, the curves for different approach the CramérRao bound for onoff detection in the antisymmetric output port. We see that the ML estimate actually beats the CR bound for small . This is because the ML estimate is biased for finite , as may be expected from the highly nonlinear nature of the estimator of Eq. (45). We recall here that bias and MSE are independent characteristics of an estimator. Many minimum MSE estimators are biased [6], and techniques for reducing MSE by introducing bias are wellknown in the signal processing literature [26]. Since MSE is the relevant metric for most applications, the behavior of the MSE observed here is a welcome – though intriguing – feature. For all values of the separation, the performance – when not superior to that predicted by the CR bound – is within a factor of 2 of it. The behavior of the simulated MSE of SLIVER with respect to the CR bound is very similar to that reported for the SPADE and binary SPADE schemes in [7].
For PSFs other than the circular Gaussian, the ML estimator takes a different form than that in Eq. (45), but depends on the measurement result only through the quantity .
e.2 Photon counting in the antisymmetric output port
For numberresolved photon counting in the antisymmetric output port, the measurement record is now a vector of nonnegative integers. The total number of photons counted – – is a sufficient statistic and the ML estimator for is:
(46) 
where the last value is again an arbitrary assignment needed in the event that a solution to the equation for the ML estimate for does not exist.
Figure 6 shows the results for the MSE of the ML estimate for a source strength photons, well outside the regime. For such a source, an onoff detector would need a very large to collect useful statistics because of detector saturation. The simulated MSE is shown scaled relative to the minimum value of the CR bound for values of 50, 100 and 150. The ML estimate again beats the CR bound for small owing to the biasedness of the estimator of Eq. (46). However, the curves for each closely approach the CR bound for photon counting in the antisymmetric output port for separations . For all values of the separation, the observed performance – when not superior to that predicted by the CR bound – is within a factor of 2 of it.
For PSFs other than the circular Gaussian, the ML estimator takes a different form than that in Eq. (46), but depends on the measurement result only through the quantity .
e.3 Effect of centroid mismatch: Onoff detection in the antisymmetric output port
We have so far assumed that the centroid has been located at the optical axis, so that the vector from the origin of the image plane to the midpoint of is zero (Fig. 1). A detailed study of the operation and performance of SLIVER when this assumption is removed will be given elsewhere. Here, we briefly evaluate the performance of onoff detection in the antisymmetric port for nonzero . It can be shown that the noclick and click probability are given by
(47)  
where
(48)  
Consider the case of parallel to the separation vector (Fig. 7). For the Gaussian PSF of Eq. (41), let so that is a normalized misalignment factor. For such a configuration, and for fixed with , simulation results for the MSE of the estimator of Eq. (45) are shown in Fig. 8. Note that the estimator used cannot depend on the centroid position, as this is unknown in practice. We have assumed photons.
Comparing to Fig. 5, we see that the misalignment has little effect on the MSE for separations . For smaller separations, the MSEs in Fig. 8 are greater than those in Fig. 5 and the MSE does not go to zero for zero separation. However, the MSEs still beat the CR bound for the aligned SLIVER system, suggesting they also beat the CR bound – yet to be determined – for the misaligned system.
Figure 9 displays the simulated MSEs for the same source and system parameters, with and for the misalignment factor values and 0.4. We see that the performance degrades gracefully with increasing misalignment, remains within an order of magnitude of the CR bound of the aligned system, and does not display the divergent behavior of imageplane photon counting [3, 4, 7].
A detailed study of SLIVER under misalignment – including derivations of the formulas in this section, CramérRao bounds, and performance analyses – will be given elsewhere. Here, we simply note the important point that the MSE of a misaligned SLIVER system is still finite for arbitrarily small separations, and Rayleigh’s curse has not returned.
Appendix F Discussion and conclusions
Motivated by the greater familiarity and facility of analysis afforded by the semiclassical theory of photodetection over the full quantum theory, we have set up the problem of resolving two incoherently radiating point sources in this framework. This theory is quantitatively exact for such sources and for measurements involving only linearoptics processing prior to photodetection. The sources are assumed to be of arbitrary strength and circularlysymmetric PSF’s directly related to real imaging systems are considered. Since the analysis in [7] was largely confined to a onedimensional model of the pointspread function with a weaksource assumption, the model used here is of interest both theoretically and in practice (See also the remarks in Sec. V therein). We have not made full semiclassical analyses of the binary SPADE and SPADE protocols of [7] here, as these are rather involved and will be given elsewhere. Preliminary calculations indicate that the results from [7] for agree with those from the semiclassical analysis to .
We have also proposed and analyzed a novel interferometric scheme – the SLIVER method of Sec. 3 – for estimating the separation between the sources once their centroid has been located. The scheme alleviates Rayleigh’s curse and provides superresolution at all values of source strengths. The essential ingredient of the scheme is an image inversion interferometer (III) – a device that separates the symmetric and antisymmetric parts of the imageplane field relative to inversion about the centroid of the sources. The SLIVER method, unlike the binary SPADE and SPADE schemes, does not require coupling the imageplane field into one or more waveguides with mode profiles tailored to the PSF of the imaging system. As such, it appears to be more readily implementable in practice.
An explanation of the superlocalization effect of SLIVER for estimating may be given in terms of photon statistics as follows. If both sources were superimposed at the centroid (so that ), the field in the antisymmetric output port of the III is identically zero despite their incoherence. If is small but nonzero, the mean photon number is finite but still small as per Eq. (20). The Fisher information on the mean of a BoseEinstein distribution equals and hence is very large in this region of . This extreme sensitivity to changes in is partially offset by the weak sensitivity of on via which is close to zero for (cf. the term in square brackets in Eq. (25)). These two effects compensate each other to give a finite value of the Fisher information for separations close to zero. In the symmetric port, on the other hand, the mean output is large at and the Fisher information is correspondingly small, so that the second effect dominates and superlocalization is not obtained. Similar remarks apply to onoff detection in the two output ports. It is rather remarkable that appropriate linearoptics processing followed by spatiallyunresolved onoff detection can alleviate the divergent MSE behavior plaguing ideal spatiallyresolved imageplane photon counting for closely separated sources. Somewhat intriguingly, receivers using the same toolkit of operations (linearoptics processing with bulk onoff detection) also approach the quantum limits of binary [27, 28] and ary [29] communication with laser (coherentstate) light.
As mentioned in Sec. 3, various imaging modalities employing image inversion have been proposed and demonstrated previously in the microscopy literature [17, 18, 19, 20, 21, 22, 23, 24]. The importance of the antisymmetric output of the III for improving lateral resolution has also been noticed, e.g., in [18]. However, the claimed resolution enhancements in these works are in terms of the reduction of effective widths of the PSF by fixed factors, which in turn are related to the intensity patterns at the two outputs of the III (or their difference). Thus, from a statistical viewpoint, it appears that these works focus on the spatial structure of the first moment of the photocount random process (i.e, the intensity) at the two output channels of the interferometer and on the information about the separation extractable from it. In contrast, we have focused on spatially unresolved detection at the antisymmetric output port of the III and on the CramérRao bound as the fundamental limit to the mean squared error of any (unbiased) estimate of the separation. By considering statistically optimal estimators, we can exploit the information available in the full probability distribution of the measurement rather than just its first moment. Such subtler estimators are necessitated by the limited number of photons that can be collected in most imaging applications. The maximum likelihood estimator, which approaches the CramérRao bound in the limit of large [6], was shown to yield performance in close agreement to the bound and sometimes exceeding it.
In order to focus on the spatial aspects of the resolution problem, we have, as in [7], suppressed the explicit temporal dependence of the field in our analysis. In effect, we have assumed that the light emitted by the sources is in a single quasimonochromatic temporal mode in each detection window. Such a situation can be easily realized using a pseudothermal light source [16], but is in general unrealistic for many astronomical and biological imaging situations in which the coherence time of the emitted radiation is typically much shorter than the duration of the detection window. An approximate treatment of such cases can be given using the concept of the number of coherence cells [16]. This approach effectively redefines the number of shots according to the coherence properties of the light, with exact results obtainable in certain cases. These issues will be explored in detail elsewhere.
The effects of loss and nonunity quantum efficiency of the detectors can be incorporated in our analysis by appropriate scaling of the source strengths. Since spatially homogeneous losses simply scale the semiclassical field amplitudes, thermal states remain thermal states under such attenuation.Thus, a system resolving sources of strength using a detector with quantum efficiency can be modeled as a system with source strength and an ideal detector.
We have assumed through most of this paper that the centroid of the sources has been located at the optical axis. In optical astronomy applications, the centroid position may be available from previous observational records, or else a long observation time should be available to determine the centroid by, say, imageplane photon counting. In microscopy applications, the time available for measurements may be more limited, e.g., because of photobleaching in fluorescent samples[30], and lack of knowledge of the centroid position can be a significant impediment to implementing SLIVER. To address this problem, one can first perform imageplane photon counting using a portion of the available light to determine the centroid position to within a small fraction of the PSF width before implementing SLIVER. The simulations in Sec. 5.3 indicate that the MSE of separation estimation should still be well below that achievable using imageplane photon counting. A more general study of the effects of centroid mismatch is best carried out in the framework of multiparameter estimation theory and will be given elsewhere.
Apart from exploring these issues, various extensions and refinements of the schemes of this paper suggest themselves. For example, spatiallyresolved detection in one or both output ports of the interferometer can only improve the performance of SLIVER. In any case, the results here should hopefully spur investigations into these issues in order to reap the superlocalization gain promised by the SLIVER method.
Acknowledgments
This work is supported by the Singapore National Research Foundation under NRF Grant No. NRFNRFF201107.