Far-field intensity signature of sub-wavelength microscopic objects

Information about microscopic objects with features smaller than the diffraction limit is almost entirely lost in a far-field diffraction image but could be partly recovered with data completition techniques. Any such approach critically depends on the level of noise. This new path to superresolution has been recently investigated with the use of compressed sensing and machine learning. We demonstrate a two-stage technique based on deconvolution and genetic optimization which enables the recovery of objects with features of 1/10 of the wavelength. We indicate that l1-norm based optimization in the Fourier domain unrelated to sparsity is more robust to noise than its l2-based counterpart. We also introduce an extremely fast general-purpose restricted domain calculation method for Fourier transform based iterative algorithms operating on sparse data.


Introduction
Superresolving optical microscopy becomes increasingly important in medical and biological applications, in nanotechnology, material science etc.Recent advances in computational imaging [1,2] and deep learning [3] reopen the question of how much resolution can be enhanced by data completition methods [4,5,6,7,8].While scanning near-field optical * rafalk@fuw.edu.plmicroscopy [9] as well as techniques based on fluorescence microscopy [10] allow to reach a deeply sub-wavelength resolution down to the order of several nanometers, the resolution of classical optical imaging is restricted by the Abbe diffraction limit.Fluorescence microscopy gave rise to methods that overcome this limit such as stochastic optical reconstruction microscopy [11], photo-activated localization microscopy [12] or stimulated emission depletion [13,14].These techniques bypass the diffraction limit which refers to two-point resolution or the width of the point spread function but not to the localization precision of point sources [15], whereas optical detection of isolated sub-wavelength objects remains a part of these methods.By purely optical means the resolution may only insignificantly exceed the Abbe limit [16,17].Respective methods range from using a high refractive index immersion liquid and a high-numerical aperture (NA>1) objective on which the diffraction limit actually depends, through structured illumination [18] and deconvolution techniques [19,20,21], up to development of novel optical set-ups such as for 4PI confocal optics [22,23].
Resolution depends on the use of the available degrees of freedom of the imaging system [24,25] given in terms of its space-bandwidth product and independent polarization channels [26,27,16].It may be increased by using a priori knowledge about the object combined with signal modulation and reconstruction techniques [28,5].A well known encoding technique that enables to extend the measured spatial spectrum of the specimen by copying high frequency information to lower frequencies is based on the use of gratings [24,25,29,29].More generally, modulation may involve structured illumination [18], speckle pattern projection [30,31], or structured illumination varied on sub-wavelength scale [29,32,33,34].A potential novel promising approach to superresolution is based on superoscillations [35,36,37].Imaging with nanospheres has been also shown to overcome the diffraction limit [38].Finally, digital processing with deep neural networks may enhance the spatial resolution of regular microscopic images slightly beyond Abbe's limit [4].
In this paper we make use of far-field intensity pattern obtained under coherent illumination to deduce information on the shape and location of subwavelength sized objects.The specimen is placed in slits of a sub-wavelength binary metallic grating.Interference of the Fourier spectra of the object(s) and the grating enhances far-field intensity modulation introduced by the object(s).This approach resembles object recovery in deconvolution microscopy [19].It is also similar to the interscale mixing microscopy (IMM) [39,8], although we do not share the opinion [8] that a binary grating nonoverlapping with the object could introduce mixing of the evanescent spectrum of the object into the far-field pattern.In fact this concept is not upheld in [7].As compared to IMM, we use a different singal reconstruction method and consider a 2D situation.The far-field diffraction intensity pattern is captured and processed to obtain information on the object.We apply a two-step image reconstruction method for sub-wavelength objects and demonstrate that signal recovery from the farfield diffraction pattern remains feasible in 2D despite the significant drop in the fraction of diffracted light within the diffractive pattern observed when a 2D situation is compared against 1D.Numerical recovery of sub-wavelength objects from their far-field interference signatures is computationally challenging and depends on using further assumptions, for instance on the sparsity of the objects.The overall resolution and amount of details that may be recovered from the non-evanescent field is strongly limited by the signal-to-noise ratio. 2 Recovery of sub-wavelength objects from a far-field interference pattern Reconstructing an object from the intensity measurement of its far-field interference pattern is an ambiguous inverse problem.The proposed method consists of two parts.The first stage is derived from the framework of deconvolution microscopy.Then we apply an original optimization procedure with a genetic algorithm using a criterion evaluated with a restricted-domain Fourier transform.
Figure 1 shows a sample object placed together with a binary grating and introduces respective denotations for the geometric features.The field from the grating interferes with that from the object(s) placed within the slits of the grating.Intensity distribution is measured in the Fourier plane of the objective.The overall far field intensity is is the optical transfer function (OTF) and its squared modulus is the modulation transfer function (MTF) of the optical system, O(x, y) is the object field, and R(x, y) is the reference field created by the grating.The caret denotes the 2D Fourier transform.OTF is a low-pass function with a cut-off at |k| < N A • k 0 , where NA is the numerical aperture and k 0 is the wavenumber.When the object O fits in the grating slits, it is not modulated by the grating R and the far field contains a superposition of the two spectra filtered independently by the same OTF.In effect, the far-field holds no information about the spatial spectrum of the object above the cut-off and the role of the grating in the measurement is other than to shift the evanescent spectrum below the cut-off.For small objects, the interference term For instance for 2D square objects and a 2D grating with square masks the reinforcement of the intensity signal due to interference is on the order of the squared ratio of their surfaces and may be substantial.Figure 2 illustrates this enhancement for a 2D 25 × 25 rectangular grating with Λ = 275 nm, w/Λ = 50% at the wavelength of λ = 532 nm and for N A = 1.49when the size of the object is 15 nm.The interference contribution from this deeply subwavelength object is clearly a measurable correction to the far-field intensity of the grating.The interference pattern also encodes the phase of Ô as interference fringes.This justifies the use of the grating in our set-up.At the same time, the interference mechanism enhances noise in a similar way as object information.
Recording the far field interference pattern I(k x , k y ) resembles recording a Fourier hologram with a reference beam R(k x , k y ).The denotation R underlines the role of a reference beam in the holographic recording played by the far-field diffraction pattern of the grating.Similar as in deconvolution microscopy we want to recover O(x, y) from the interference pattern.For a 2D diffraction grating R can be written as where P (k x , k y ) is the Fourier transform of the elementary cell of the grating, P = const − w 2 sinc (k x w/2) sinc (k y w/2), and Ŝ/ Q depends on the distribution of cells in the finite sized grating, Ŝ(k x , k y ) = sin (k x N Λ/2)sin (k y N Λ/2), and Q(k x , k y ) = sin (k y Λ/2) sin (k x Λ/2).Here N , w, and Λ denote the number of grating periods in each dimension, width of a square grating mask, and the grating pitch, respectively, as indicated in Fig. 1.The object may be approximately recovered with the deconvolution formula where the inverse Fourier transform is calculated over part of the domain for which The result is diffraction limited.For a sub-wavelength sized object all fine details are lost, and closely positioned objects can not be resolved, but the locations of the objects may be approximately determined restricting the vast domain that will be examined in further computationally intensive optimization.The deconvolution formula (2) is used by us only to estimate the region of pixels where the object(s) may be located.
The second stage of the algorithm consists of minimizing an error function dependent on the hypotheti-cal object location and shape using the actual far-field intensity measurement.We consider both the 1 and 2 norms to construct the criterion, where O hyp is the tested hypothesis of the objects shape, and is the corresponding far-field intensity which may be compared to the actual measurement I.
In Eq. ( 3), p = 1, 2 decides upon the use of 1 or 2 norm, and Ω denotes a subset of spatial frequencies selected adaptively based on signal to noise ratio.Optimization is computationally intensive, and may only be successful with additional constraints imposed on the object O(x, y).We assume that O is sparse, and is located in the region first estimated in the deconvolution stage.The size of this region considered here is on the order of 1000 pixels, while the number of spatial frequencies in Ω is on the order of 10 2 −10 4 (we used the value of 300 most of the time).Optimization is further simplified for binary objects of a known and fixed or parametrized shape.
In practice Eq. ( 3) has a large number of local minima and can not be minimized with gradient descend methods.In Fig. 4 we show a typical shape of F p near the global minimum, as a function of two parameters corresponding to the possible size of the object when its location is already known.The minimum is broad, and some local minima are also present.The situation becomes a lot more complicated when more degrees of freedom need to be included in optimization.In practice the number of local minima makes optimization difficult.
For this reason, we have used a genetic algorithm for finding the minimum, which had the additional advantage of the ease to encode additional structural information of the objects.The genetic algorithm was also more robust than Gerchberg-Saxton type iterative optimization.
Finally, it is crucial to minimize the evaluation time of the criterion.Equation ( 3) is formulated in the Fourier domain, but the constraints on O(x, y) about object location, sparsity, shape or parametrization can be only easily specified in the image domain.Thus a 2D Fourier transform has to be calculated every time, when we want to calculate the criterion.It would extremely inefficient to work with dense zeropadded matrices with high resolution sampling and use the FFT algorithm for this purpose.Instead we are calculating the discrete Fourier transform directly over small subsets of signal and spectral domains.This approach was faster by two to three orders of magnitude than using FFT from a highly optimized FFTW package included in Matlab.The details of the genetic algorithm are described in the next section and the details of the restricted domain Fourier transform in the Appendix.

Computational imaging beyond the diffraction limit
In order to minimize the cost function defined in Eq. ( 3) with respect to positions and sizes of objects we utilize a dedicated genetic algorithm.We consider a population of solutions, each of which contains information about positions and sizes for a set of objects.
The objects have a discrete pixelated structure and are located within the area of 1000 pixels selected with the deconvolution formula (2).The initial population is generated randomly with a uniform distribution of objects over the allowed area.Then, in each algorithm iteration, a new population is generated with the better half of the solutions preserved, and the other half regenerated with three genetic operators, which are mutation applied to object position, mutation applied to object size, and crossover.The genetic algorithm can be easily adapted to work with parametrized objects.In this work we also analyze arbitrarily oriented line-shaped objects fully characterized by the coordinates of their ends and their width.
Figure 5(a) presents an example of three subwavelength objects reconstructed with the genetic algorithm.The objects and their preliminary localization with deconvolution are shown in Fig. 3.The algorithm does not always converge to the exact solution  2).Isolated objects can be identified and localized but the image is diffraction limited and contains a spurious mirror image.The deconvolved image is the result of the first stage object recovery, which is later refined by numerical optimization beyond the diffraction limit.but clearly allows to achieve a localization accuracy on the order of λ/10 and to resolve objects located at distances smaller than the diffraction limit.In Fig. 5(b) we show a similar example but there is a single linear object parametrized with 5 coefficients representing the coordinates of the two ends and the thickness.The genetic algorithm is aware of the possibility to parametrize the object and the coefficients are recovered without error.Optimization with the genetic algorithm is also shown in the supplementary materials (See Visualisation 1).

Influence of noise on image reconstruction
Presence of noise in the far-field intensity measurement has a profound influence on the possibility of object reconstruction and on its quality.As a measure of noise intensity we will use the ratio of its standard deviation σ n to the intensity of the far-field averaged over the measured region < I > (the same normalization of far-field contributions was used in Fig. 2).
Noise is complex and affects the complex field addi-   tively.Due to the presence of noise, optimization of criteria F 1 and F 2 is not equivalent, and there is no guarantee that either of them has still a global minimum for the correct hypothesis about the microscopic object shapes and locations.Criterion F 1 based on the 1 norm is more robust to noise than F 2 .We note that this observation can not be simply attributed to linking object sparsity to the value of the 1 norm because the criteria are formulated in the Fourier rather than object domain, and only sparse objects are considered in the genetic algorithm.In Fig. 6 we compare the two criteria calculated in the presence of noise for the hypothesis about object shape or location near the correct values.The global minimum to the optimization problem becomes shallower with the increasing level of noise.By comparing ∂F i /∂σ n with ∂F i /∂ , where represents location error or size error, we notice that F 2 criterion loses sensitivity to these errors much faster than F 1 when noise is present.
A fundamental question is up to what level of noise σ n , the criteria F 1 and F 2 retain the global minima at the correct location or in its close vicinity.When this level is exceeded, it will no longer be possible to recover the objects, independently of the computational resources available.In practice though, for more complicated problems it may be extremely difficult to find the global minimium, and especially a flat and shallow minimum may be easily omitted by numerical optimization.
We will now focus on localizing a single binaryvalued linearly-shaped sub-wavelength sized object positioned on the discrete rectangular pixel grid (such as shown in Fig. 5(b)).A regular object of this kind could be interesting for data-storage or security applications, where a simple far-field interference measurement provides more information than is possible to see with an optical microscope at the same wavelength.An object that is fully parametrized by 5 integer numbers is a convenient test example since for these number of degrees of freedom we are able to compare the operation of an iterative algorithm with a brute-force optimization of criteria F 1 and F 2 over the entire parameter space.Such a comparison, calculated at various noise levels, is presented in Fig. 7.The fidelity of the object reconstruction with values between 0 and 100% is defined as where S obj represents the actual surface of the objects and S err represents the incorrectly identified object surface (including both omitted and erroneously attributed areas).Fidelity of 100% corresponds to a perfectly identified object, while the fidelity of 0% usually signifies that less than half of the object pixels have been correctly identified.Two important conclusions may be drawn from the results in Fig. 7.The first is that the 1 -norm based optimization is more robust to noise than the 2 based optimization.The second is that sensitivity to noise puts a severe limitation to the possibility of sub-wavelength sized object recovery from far-field intensity information.Since we have tested the full possible parameter space for a rather idealistic object parametrization, we do not expect any optimization method, including methods of compressed sensing or deep learning, to overcome these limitations.At the same time, the recovery is feasible, if only the level of signal-to-noise ratio is sufficiently large, which may be potentially achieved by temporal or spatial signal averaging, limiting the aperture, improving experimental stability, or other noise-reduction techniques.

Conclusion
We have examined the possibility to recover geometrical information on sub-wavelength sized objects, not limited to their location, from a far-field interference pattern obtained under coherent illumination.The far-field signature of microscopic objects considered by us does not include spatial frequencies beyond the cut-off, i.e. corresponding to evanescent waves or spatial frequencies lost in classical imaging optics with an objective having a given numerical aperture (with N A = 1.49assumed in the presented results).
We have proposed a two-step object recovery algorithm, with the first diffraction-limited step based on deconvolution, and the second based on numerical optimization.This second step involves the use of a genetic algorithm and a restricted-domain Fourier transform (described in the Appendix) the purpose of which is to speed-up calculations of the far-field for sparse objects when only a limited part of the Fourier coefficients are required.Overall, the method makes it possible to recover sub-wavelength sized objects with sizes on the order of λ/10 from far-field information, although for more complicated scenes, the method becomes computationally intensive and does not always converge to the correct result.Sub-wavelength object recovery is only possible at a very low level of noise.Even with a known object parametrization, a small level of noise introduces false minima to the cost functions making the true hypothesis impossible to distinguish from alternatives.This limitation concerns object recovery with both 1 -norm and 2 -based criterion functions and is unlikely to be mitigated by the use of sophisticated optimization frameworks such as compressive sensing or deep learning.At the same time, the robustness of 1 -norm based criterion to noise is considerably better than that of the 2 -based counterpart.Compressive sensing heavily relies on 1 norm as a measure of signal sparsity which leads to convex computationally tractable optimization criteria.Perhaps noise robustness properties of the 1 norm are in effect underestimated.Our optimization of 1 -norm applied to dense Fourier domain information is probably original.
It is unlikely that information recovery from a far field signature could become a significant alternative to existing superresolving wide-field microscopic techniques.Prospect applications may be connected with security purposes, where a specific microscopic object impossible to be analyzed by classical optical microscopy could be detected, localized and verified by the examination of the far-field diffraction pattern.Another area of interest is high density data storage where the geometry of sub-wavelength sized features could be used to enhance the amount of information readable without shortening the laser wavelength.

Figure 1 :Figure 2 :
Figure 1: (a) Subwavelength-sized object (in red) placed in a slit of a 2D binary grating.(b) Far field intensity expressed in dB (and obtained for an object shown in Fig.3(a).The dynamic range of the measurement is limited to 30 dB which results in masking overexposed low-frequency information.

Figure 3 :
Figure 3: Recovery of sub-wavelength sized objects from a far-field interference pattern by deconvolution: (a) Three 15 nm objects placed on the grating (length expressed in the units of diffraction limit).The corresponding far-field interference pattern is shown in shown in Fig. 1(b) (b) Objects recovered from the far-field interference pattern using Eq.(2).Isolated objects can be identified and localized but the image is diffraction limited and contains a spurious mirror image.The deconvolved image is the result of the first stage object recovery, which is later refined by numerical optimization beyond the diffraction limit.

Figure 4 :
Figure 4: Criterion F 1 calculated in the case when two square objects with sizes 53 × 53 nm and 23 × 23 nm are located in close vicinity.The tested hypothesis assumes that we know the exact object locations and the size of the larger object.F 1 is analyzed as a function of the size of the smaller object.This example shows that although the global minimum appears for the correct hypothesis about the size of the object, other local minima also exist.

Figure 5 :
Figure 5: Recovery of sub-wavelength sized objects from a far-field interference pattern with a genetic algorithm (red).The potential areas where the objects may be located (marked in yellow) were initially determined by deconvolution The actual object locations are shown in green.(a) Three 15 nm objects found approximately (See.Fig. 3(b.);(b) single linear object parametrized with 5 coefficients found exactly.Optimization with the genetic algorithm is also shown in Visualisation 1.

Figure 6 :
Figure 6: Sensitivity of 1 and 2 norm based criteria to noise and distortions.Criteria F 1 and F 2 are calculated in the presence of noise when the actual object consists of a single 30 nm square spot (a),(b) or of two such objects separated by 30 nm (c).The change of the shape of the minimum with noise indicates that 1 -based criterion is more robust to noise.(a) the tested hypothesis assumes a correct size but includes a position error; (b) the tested hypothesis assumes a correct position and proportions but includes a scale error.c) the tested hypothesis assumes correct positions and sizes but unknown separation.

Figure 7 :
Figure 7: Reconstruction fidelity of a sub-wavelength linear object found from far-field interference pattern in the presence of noise by optimization of the criterion defined in Eq. (3) with either p = 2 (a),(c) or with p = 1 (b),(d).An example of the object is shown in Fig. 5(b).(a),(b) results corresponding to a global minimum obtained by brute-force optimization, (c),(d) results obtained with a genetic algorithm.The object consists of a single 33 nm-thick line with length varied between 30 and 120 nm.Results are extremely sensitive to noise and indicate a better noise robustness obtained with 1 norm than with a 2 norm.