Supergrowth and sub-wavelength object imaging

We further develop the concept of supergrowth [Jordan, Quantum Stud.: Math. Found. $\textbf{7}$, 285-292 (2020)], a phenomenon complementary to superoscillation, defined as the local amplitude growth rate of a function being higher than its largest wavenumber. We identify the superoscillating and supergrowing regions of a canonical oscillatory function and find the maximum values of local growth rate and wavenumber. Next, we provide a quantitative comparison of lengths and relevant intensities between the superoscillating and the supergrowing regions of a canonical oscillatory function. Our analysis shows that the supergrowing regions contain intensities that are exponentially larger in terms of the highest local wavenumber compared to the superoscillating regions. Finally, we prescribe methods to reconstruct a sub-wavelength object from the imaging data using both superoscillatory and supergrowing point spread functions. Our investigation provides an experimentally preferable alternative to the superoscillation based superresolution schemes and is relevant to cutting-edge research in far-field sub-wavelength imaging.


I. INTRODUCTION
As we strive to engineer optical systems with better imaging capabilities, much attention has been devoted to surpassing the Rayleigh resolution limit in diffraction limited optical devices [1,2].To achieve superresolution, methods such as evanescent field based techniques [3,4] and negative refractive index materials [5] have been proposed.To that end, the efficacy of superoscillatory spot generation for label free far-field based superresolution imaging is well established.Superoscillation (SO) refers to the phenomena of local oscillation frequency of a function being faster than its fastest Fourier component [6][7][8][9][10][11]. Thus, in optical systems, superoscillation can be utilized to generate sub-wavelength hot-spots and thus beat the Rayleigh resolution limit [12][13][14].Sub-wavelength spots have been realized at the focal plane of a microscope objective using optical eigenmode approach implemented with a spatial light modulator [15].Also, superoscillatory spot properties of radially polarized Laguerre-Gaussian beams in a confocal laser scanning microscopy setup have been studied [16,17].Optimization of superoscillatory lenses for sub-diffraction limit optical needle generation has been investigated as well [18].
In principle, it is possible to design arbitrarily small superoscillatory optical spots.However, superoscillation is necessarily accompanied by enhanced side lobes.This leads to poor quality in imaging and unrealistic constraints on the dynamic ranges of the detectors.Several recent ventures try to address this problem through numerical or design based approaches.Simulations for simultaneous optimization of superoscillatory spot size and their relative intensities compared to the side lobes have been performed [19].Elimination of sidelobes along a particular dimension by introducing moonlike apertures has been demonstrated [20].
In this work, we show how the problem of enhanced sidelobes can be circumnavigated by utilizing supergrowing functions, a concept proposed by Jordan [21].The phenomenon of supergrowth (SG) is analogous to superoscillation.While superoscillation pertains to the local oscillation rate of a function, supergrowth occurs when the local growth rate of the amplitude is higher than the highest wave number in the Fourier space of a bandlimited function.A large local growth rate enables enhanced spatial resolution.The idea is analogous to evanescent wave imaging microscopes, but is applicable in the farfield.In our analysis, we consider a canonical singleparameter bandlimited oscillatory function and locate its SO and SG regions.In addition to the entire SO and SG areas, we identify the near maxima regions for both phenomena.We provide analytical estimates for the intensities in these regions and find that the amount of light in the SG areas is exponentially higher compared to the regions associated with SO.Finally, we present two parallel schemes to reconstruct an incoherently illuminated subwavelength object with SO and SG point spread functions (PSF).We numerically compare our approaches to object reconstruction using a bandlimited sinc PSF.As expected, superresolved object reconstruction is achieved in the former two cases.
The past two decades have witnessed considerable progress in the superoscillation related research.The far-field nature of superoscillatory fields have been shown by proving that subwavelength structures generated by a diffraction grating can survive farther than the evanes-cent waves [22].This investigation has inspired further studies into Schrödinger equation based evolution of superoscillatory waves [23][24][25].The correspondence between superoscillations and weak values [26,27] is well known.On a related note, super-phenomena in arbitrary quantum observables [28] have been proposed.Towards the more implementational side, numerical optimization of the energy ratio between superoscillatory region and total signal has been studied [29].Also new methods for generating superoscillatory functions have been proposed [30][31][32].In comparison supergrowth is a very recent concept.Jordan [21] showed that it is possible to access superresolving features using supergrowth.Spherical Bessel function based method for systemic generation of SO/SG functions and general approximation scheme using bandlimited functions have been prescribed [33].Our analysis draws inspiration from these works and solidifies the benefits of implementing supergrowth based superresolution imaging in practice.
This article is organized as follows.In Sec.II, we briefly describe the phenomena of SO and SG using a canonical oscillatory function.Sec.II A compares the length of SO and SG regions for the chosen function and Sec.II B compares their intensities.In Sec.III, we present the schemes to reconstruct a sub-wavelength object using both superoscillatory and supergrowing spots.In Sec.IV we discuss the implications of our findings.We conclude in Sec.V.

II. PROPERTIES OF SO/SG FUNCTIONS
For our analysis, we consider the function [22] f (x) = (cos parameterized by a positive real number a, which sets its SO and SG properties, and a natural number N , which gives an upper bound to the Fourier wavenumber k max = N for the function.This corresponds to a shortest wavelength of oscillation of λ min = 2π/N .In the following sections we characterize the SO and SG regions of f (x) and look at the intensities within.

A. SO/SG regions
For any complex valued function f (x), the local wavenumber is defined as k(x) = Im ∂ x log f (x), and the local growth rate is κ(x) = Re ∂ x log f (x).If f (x) is bandlimited with highest wavenumber k max , the function is said to be superoscillating (supergrowing) at For f (x) in Eq. ( 1), local wavenumber and growth rate take the form [21]  and For √ 2 − 1 < a < √ 2 + 1, this function shows SO, but there is no SG.For a > √ 2 + 1, the function displays both SO and SG behavior.We will show in Sec.II B that larger values of a show more SG but less SO, so we will consider only such cases.The function shows SG behavior (κ(x) ≥ N ) in x ∈ [arctan y l , arctan y h ], with The largest growth rate occurs when x = ±x sg = ± arctan(a −1 ), giving a maximum growth rate of For large values of a, this approaches ±N a/2 rapidly.

B. Intensity comparison
In imaging applications, SO/SG spots with higher intensities are ideal since that makes superresolution imaging less susceptible to the influence from sidelobes.In this section, we quantify the amount of light in f (x) in regions showing SO/SG behaviors.For simplicity, we consider the function f (x) to be the intensity point spread function of an incoherent imaging device [34].The total intensity within a single period of f (x) can be defined as For large a and N , the total intensity can be approximated as (see Appendix.A) We can further simplify the expression above using Stirling's approximation for factorials We see indeed that I 0 increases exponentially in N .From Fig. 1, it is clear that the most of this light is available away from the origin.This is a manifestation of pronounced side-lobes around a superoscillatory spot.However, as the subsequent analysis will reveal, SG based imaging can be highly advantageous in this respect.We see at the superoscillatory region x = 0, while the region of superoscillation (where , the useful range is more restrictive.This is where the oscillations are at their fastest, and the amplitude of oscillations is approximately constant.Noting that ln f ≈ iaN x + 1 2 N (a 2 − 1)x 2 , near x = 0, we get a more restrictive range x ∈ [−1/ N (a 2 − 1), 1/ N (a 2 − 1)].We can find the amount of power in the range by integrating |f | 2 over this restrictive interval to find (Appendix.A), We note that the intensity in this region is exponentially suppressed compared to the total intensity in Eqs.(7), (8).
Taking a more relaxed view and using the whole region of superoscillations (even if it becomes impractical to use), we find that integrating over the range [− arctan 1 √ a , arctan 1 √ a ] gives an intensity of exponentially smaller than Eq. ( 8) by a N .Let us now consider the SG region characterized by Eq. ( 4).We find, for a 1, κ has the expansion around x sg of κ ≈ (aN/2)(1 − (a 2 /2)(x − x sg ) 2 ) to quadratic order.Consequently, we get considerable supergrowth of approximately exponential form in the region x ∈ [x sg − 1/a, x sg + 1/a].Despite the fact that the phase of f oscillates several times in this region, it is of no concern to us, since we are using only the magnitude (intensity) for the imaging.
Even better exponential fits can be obtained by reducing the range by (say) a factor of 2, but at the cost of reducing the amount of light that is used.The intensity in this reduced region is which to order 1 a is, It is interesting to note that the intensities in restricted regions Eqs. ( 9) and ( 12) both have 1/a dependence, instead of the exponential dependence of total intensities in Eqs. ( 7) and (10).This is due to the shrinking of the restricted regions as a is increased.The coefficient of 1/a in the Eq. ( 12) is much larger than the corresponding coefficient in Eq. ( 9).Therefore, the near maxima SG region contains more intensity compared to the near maxima SO region.The intensity in the total SG region (κ > 1 and specified by Eq. ( 4)) can be approximated as We see that the intensity of the SG region is exponentially higher than that of the SO region by a N , but exponentially smaller than Eq. ( 8) by (1/2) N .Additionally, the great advantage is that we can profitably use the entire SG range, but only a small fraction of the SO range.This is because the amplitude varies too much once x leaves the more restricted range, as well as the fact that lower values of k become mixed in with the high values (even if still above N ).The SG region has none of those drawbacks.In Fig. 2, we show quantitative comparisons of SO and SG properties of our chosen function for different values of the parameters a and N .We see that SG intensity dominates SO intensity for most of the parameter range, except for a very small region in the bottom left of the plot.This signifies that for supergrowth imaging can provide a sufficient intensity to obtain bright images, overcoming signal-to-noise ratio issues of superoscillation.We also can compare the lengths of the SO and SG regions with varying a in the bottom panel of Fig. 2. Both of these lengths are independent of N .For large values of a, the length of the SO region ∝ 2a −1/2 and approaches 0 asymptotically.However, the length of the SG region ∝ π/4 − 2a −2 and rapidly approaches π/4 = 1.25λ min asymptotically for large a, which can also be seen in the bottom panel of Fig. 1.

III. SUPERRESOLUTION IMAGING AND OBJECT RECONSTRUCTION
In this section, we present schemes for the detection of sub-wavelength features using SO and SG spots in a diffraction limited optical system.For a simple onedimensional model of incoherent imaging, we define the function O(x) for the object intensity, and S(u) for the image intensity.They are related via the intensity PSF of the imaging system, χ(x), which indicates the image created from a point source.Imaging theory dictates so the image created is a convolution of the object with the PSF [1,35].
We consider the problem of reconstructing O(x) from experimental imaging data S(u) [36].Now, the Fourier transforms 1 of the above functions are related by If the imaging system is illuminated by light of wavelength λ 0 = 2π/k 0 and has an exit pupil numerical aperture of NA, the PSF and therefore S(k) in Eq. ( 15) is bandlimited by k 0 NA.Inverting the above equation Õ(k) = S(k)/ χ(k) leads to a loss of information of the higher spatial frequency components and places a constraint on our ability to resolve features smaller than λ 0 .However, in the following investigation we show how isolating the SO/SG region of the function in Eq. ( 1) and using that as the PSF can help us resolve sub-wavelength features.This could be done in a confocal imaging setup, where the object is illuminated by a SO/SG optical field spot.As we will see, the nature of SO/SG field spot dictates the near origin features of the image, thereby providing an access to the subwavelength features in the object.In a confocal setup with one objective, the objective NA and wavelength will determine the extent of the spatial filter function that results from the aperture in front of the detector.A pinhole with a deeply subwavelength diameter will still map to a diameter of ∼ λ N A on the object.NA mismatch between illumination and collection is one way to allow for narrow filter functions with respect to SO/SG regions.
In the following analysis, we assume the object has an extent x ∈ [−l/2, l/2] with l < λ 0 .We can also scan an extended object gradually by filtering length l at a time.For further simplicity, we assume NA = 1 and the object is symmetric (i.e., O(x) = O(−x)).Also, in the sample cases we consider, the amplitude PSFs are bandlimited by N .The subsequent analysis separately shows how 1 Here adopt the convention that the Fourier transform of a func- e −ikx g(x), such that the inverse e ikx g(k).
Reconstruction of two sub-wavelength peaks using the methods in Sec.III.The first panel shows the object, which we assume to be two Gaussian peaks separated by a distance of ∼ 0.08λmin.Throughout the simulations, we assume N = 10.The second panel shows the PSF (Ref (x)) 2 for a = 6 (solid red) and its approximation (green dashed) cos 2 60x near the origin.The third panel shows the PSF |f (x)| 2 (solid red) and its approximation (green dashed) e 119.2x centered at x = xsg for a = 12.The fourth panel shows object reconstruction with a bandlimited sinc PSF using Eq. ( 15).The fourth and fifth panels show the results from reconstruction using the methods laid out in Secs.III A and III B respectively.We see the superior sub-wavelength feature detection ability of our reconstruction scheme.

SO or SG based imaging can help us reconstruct O(x).
We assume, for SO based imaging, the intensity PSF is (Ref (x)) 2 while for SG based imaging the relevant PSF is |f (x)| 2 .

A. Object reconstruction with SO spots
We consider the real part of Eq. ( 1) as the PSF superoscillating at origin with local rate k = aN N .Thus near the origin, for a length larger than the object length l, χ(u − x) ∼ cos 2 k(u − x) is a good approximation for the PSF and Eq. ( 14) can be written as where for symmetrical objects Õ(2k) = Õ(−2k) and C is a constant.The subscript k on S expresses the fact that the observed image intensity depends the chosen PSF.We can invert Eq. ( 16) or measure the image intensity peak reduction to calculate Õ(2k).Thus, Eq. ( 16) provides a prescription for the direct measurement of the high spatial frequency Fourier coefficients of the object.
We can perform a series of measurements with different values of k, i.e. a series of different PSFs, to map out the Fourier transform of the object.If k SO denotes the highest value of the local wavenumber in experiment, the reconstructed object is B. Object reconstruction with SG spots In this case, we consider the function in Eq. ( 1) shifted by x sg such that near origin the PSF can be approximated as χ(u−x) ∼ e κ(u−x) , where κ = N (a−1/a) N .Then, we have, Assuming a Fourier expansion of Eq. ( 18) transforms at u = 0 to (20) Note as n increases, the contributions from higher order terms decrease.In practice the sum on the right hand side will be a good estimate for the left hand side if we consider terms till n = κl.Again we consider a series of different SG PSFs to measure the the left hand side of Eq. ( 20) as a function of κ.We can then solve for the coefficients numerically to reconstruct the object O SG (x).
An alternative approach would be to approximate the image intensity in Eq. ( 18) as a Laplace transform of the object.By performing intensity measurements for different supergrowing PSFs, we gain access to the Laplace transform at different values of κ.The object then can be reconstructed by performing a numerical inverse Laplace transform [37].
One drawback of these schemes is the need to perform imaging using multiple PSFs.In practice, we can limit our experiment to a few different PSFs and use interpolation to approximate Õ for the intermediate values.
Interestingly, the SO based reconstruction relies on intensity in a small region near origin.Therefore, the method's effectiveness heavily depends on the detector's ability to register small changes in the intensity.On the other hand, the SG based approach does not have this limitation because we only look at the image intensity at a particular point.In experiments, of course, both of these methods will be limited by noise.Having more intensity, SG gives much better SNR.
We show the effectiveness of the schemes in Fig. 3.The top panel shows a sub-wavelength object which is reconstructed using a bandlimited sinc function and inverting Eq. ( 15), and also using SO/SG based schemes described in this section.The second and third panel compare the PSFs for SO and SG imaging for highest values of a (6 and 12 respectively) with their approximation.As we have already seen in Fig. 2, the length of the SO region shrinks with increasing a.This restricts the values of a we can use in practice for SO based reconstruction, while the SG based approach does not have this limitation.In the next panel we immediately see that the cardinal sine function is unable to resolve the two peaks in the reconstructed object as expected.However, as the last two panels show, both the SO and SG PSFs are able to reconstruct the peaks.In the SO based approach, the results of reconstruction are limited by our ability to generate long superoscillatory features for higher values of a. Whereas, in SG based approach, the numerical accuracy in determining the coefficients c n in Eq. ( 19) limits the reconstruction.These shortcomings in both cases lead to slightly shifted estimation of the peaks.

IV. DISCUSSION
It was previously shown [21] that the function in Eq. ( 1) is capable of superresolution imaging using the superoscillation and supergrowth using the function over a full period.For the case of two nearby point sources, putting the point sources in the "sweet spot" of the SO or SG region (maximum k(x) or κ(x) respectively) can give superresolution of the two sources.However, due to the nature of the function, any apparatus aiming to use this method would require extremely high dynamic range, both on the illumination and detection side, which limits its utility.We overcome this limitation on the detection side by applying spatial filtering of sweet spot and scanning the image.While both behaviors can resolve features better than the Fourier bound, the SG phenomena benefits from exponentially higher intensity of the illumination point spread function, whereas using SO behavior would require incredibly sensitive detectors due to the exponentially suppressed amplitude.There is also the concern of hitting the diffraction limit for the SO PSF, due to the shrinking length of the region as the parameter a is increased, which is not an issue for the SG case.This demonstrates the superiority of SG spots as a means to achieve superresolution.

V. CONCLUSIONS AND OUTLOOK
We describe supergrowth as a concept analogous to superoscillations, except we evaluate the local growth rate instead of the local wave number of a function.We characterize the supergrowing and superoscillating regions of a canonical oscillatory function as well as provide an-alytical approximations for the energy inside total supergrowing, total superoscillating, near maxima supergrowing and near maxima superoscillating regions.Our analysis reveals that the supergrowing regions can contain intensity that is exponentially larger in terms of the highest local wavenumber compared to that of the superoscillating regions.These results indicate that superresolution imaging using supergrowing spots could be more advantageous compared to superoscillation based superresolution imaging.Finally, we numerically show that the SO and SG based superresolution imaging is able to reconstruct objects beyond the diffraction limit, thereby demonstrating the efficacy of our schemes.These findings highlight supergrowth as a potentially superior far-field superresolution scheme.Physically, superoscillation corresponds to rapidly oscillating regions with small amplitude.However, supergrowth, due to high growth rate, connects the smaller amplitude regions with higher amplitude ones.Therefore, it is not surprising that supergrowing regions can contain significantly more light while also giving enhanced resolution.The analytical approximations for intensities characterize the parametric dependence of the light within different regions.These expressions can be useful for designing optimized superoscillatory or supergrowing lenses.Furthermore, the object reconstruction schemes provide the framework for experimental implementation of supergrowth and superoscillation based superresolution imaging.
Our results offer a variety of avenues for potential future investigations.The analysis presented here is restricted to a simple 1-d oscillatory function.Performing similar calculations for more general and higher dimensional functions could have greater significance for experimental implementation of supergrowth imaging.On that note, the natural next step is to investigate ways to generate optimized supergrowing spots.Another relevant problem is to explore whether simultaneous occurrence of superoscillations and supergrowth could be leveraged for better resolution with suppressed side lobe intensity.Apart from the prospective theoretical ventures, a lab realization of supergrowth based superresolution imaging is an imminent experimental challenge.Therefore, our analysis serves as the groundwork for a novel superresolution and object reconstruction scheme with a multitude of scopes for further prospects. where ˆ1 0 dy y 2N (1 + y 2 ) N +1 , (A25) can be evaluated using the recurrence relation and F 0 = π/4.We can find the asymptotic value of F N by transforming the integral in Eq. (A25) with y = tan θ.Then, F N = ´π/4 0 dθ sin 2N θ.The integral can be expressed in terms of incomplete Beta function as ).As N increases, the integrand approaches zero.We approximate sin θ ≈ 1

Figure 1 .
Figure 1.The solid orange (blue dot-dashed) curve in the upper panel shows the magnitude (real part) of the function defined in Eq. (1) in log scale for a = 6 and N = 10.The vertical grid in the top panel has a spacing equal to the smallest wavelength allowed by the bandlimit.Clearly near x = 0 the function oscillates on a scale smaller than the smallest wavelength possible.This is also apparent from the bottom panel where the the local rates of oscillation from Eq. (2) (blue dotdashed line) and growth from Eq. (3) (solid orange line) and are shown as functions of position.The black dashed lines indicate the bandlimit.In both panels, the shaded regions indicate where the function is superoscillating (blue dots) and supergrowing (orange hatching).

Fig. 1
shows the behavior of the function in Eq. (1) (top panel) and the local rates of oscillation and growth in Eqs.(2) and (3) (bottom panel) for a = 6 and N = 10.The function is superoscillating (k(x) ≥ N ) in x ∈ [arctan(1/ √ a), arctan(−1/ √ a)], shaded orange in Fig. 1 (bottom panel).The largest local wavenumber is k(x = 0) = aN , and it is clear (from top panel) that near origin the function oscillates on a scale shorter than the shortest wavelength.

Figure 2 .
Figure 2. (Top) Ratio of intensities of total SG and SO regions for different values of a and N .With the increase of either a or N , the intensity within the SG region becomes increasingly higher.The dashed black line on the bottom left corner shows the parameter values with equal SO and SG intensities.The green asterisk identifies the region of parameters chosen for object reconstruction scheme presented in Sec.III.(Bottom) Comparison of the lengths of SO and SG regions as a function of a.While the SO region shrinks with the amount of superoscillation (i.e.a), the SG length approaches 1.25λmin.Note that the full range is N independent.