|
1.IntroductionThe calculation of electromagnetic fields in focal regions has been the subject of an immense amount of work over, at minimum, the last century. A formalism allowing the solution of this problem has been developed by Ignatowsky,1 Luneburg,2 and Richards and Wolf.3 This formalism has been applied to a variety of focusing problems, from studying the structure of tightly focused fields obtained by focusing linearly polarized light3 to more complex applications such as focusing Gauss–Laguerre4 and azimuthally and radially polarized5 beams, for example. The formalism has also been applied to focusing through stratified media (see, for example6–8). This is an important extension to focusing theory due to the large number of applications where this is relevant, such as image formation in high numerical aperture confocal microscopes,9 optical coherence tomography (OCT),10 and photoacoustic imaging,11 to name just a few. Many contributions have been made to the problem of focusing vector beams through stratified media (see, for example7,8,12–15). A barrier to researchers in the field putting these solutions to use is their implementation in computer software. Although the integral equations are more detailed than complex, their implementation may appear daunting to researchers from outside the field. Furthermore, even slight modifications to the beam type require a good understanding of the underlying formalism to correctly modify the equations governing the form of the focused field. One of the reasons for this is that the computational complexity of calculating the form of focused fields has hitherto dictated that as much as possible of the computation be performed analytically. The result of this is that it can be difficult for nonspecialists to modify the theory for their application and the final equations can lose connection with the starting point of the focusing formalism, which can be understood conceptually. Calculating the field at each point in the focal region requires independent evaluation of a two-dimensional (2-D) diffraction integral. Early calculations of fields focused through a stratified medium employed an Intel Pentium processor operating at a frequency of 90 MHz.12 Numerical solution of the integrals presented in that work took advantage of the rotational symmetry of the optical system. In particular, integration about the axis of revolution of the optical system (i.e., the optical axis) was performed analytically, leaving only one-dimensional (1-D) integrals to be numerically evaluated. Aside from the elegance of the analytic solution, this was also necessary to keep the computational cost of evaluating the focused fields feasible. The downside of using this semianalytic approach is that the process of making the required coordinate system transformations and manipulations results in a series of equations that are not readily physically interpreted and are not easily modifiable to treat different beam types or optical systems. The results contained within this paper were calculated on an Intel Xeon processor with 12 cores operating at a frequency of 2.3 GHz, thus reflecting a significant increase in processing power over the past two decades. This increase in computational power makes it feasible to solve the problem of tight focusing in the presence of a stratified medium largely numerically. This carries several advantages. First, it enables a more direct mapping between the equations governing the solution and the numerical implementation, making the solution easier to extend to alternative beam types and optical systems. Second, the focusing problem can be considered as a vectorial angular spectrum of plane waves, thus allowing the focusing and stratified media problems to be naturally decoupled. This provides for a very convenient way of calculating the electromagnetic field throughout a stratified medium, not just the transmitted field as has previously been the case.6 The aim of this work is to provide a formalism, along with freely available computer software,16 which is able to be used, with little prior knowledge, to calculate a range of focusing problems, possibly involving a stratified medium. Furthermore, the formalism is structured in a manner that prioritizes conceptual understanding of the code rather than computational efficiency. This allows the formalism to be extended by users with an understanding of the underlying theory and not necessarily a detailed knowledge of scientific computing, which is why we have presented this work in the form of a tutorial. Furthermore, it is our intention that this tutorial, along with the associated computer code, will enable researchers from a broad range of fields to easily solve problems related to high numerical aperture focusing. In the rest of this paper, we introduce the mathematical basis of our method for calculating the focused field in the vicinity of a stratified medium that is adapted from that previously proposed by van de Nes et al.7 We begin with the Debye–Wolf integral before considering the matrix formalism for treating the propagation of plane waves within a stratified medium. We then show how these two formalisms can be combined before considering the numerical implementation of the resulting formalism. We illustrate the use of the model through three examples: image formation in OCT, the focusing of cylindrical vector beams through a dielectric interface, and a simulation of aberrations typically experienced in confocal microscopy. We note that the stratified media focusing code, along with code to generate the results presented in this paper, is freely available for download.16 2.Mathematical FoundationThroughout this paper, we represent the electric field vector at position , as a time-harmonic quantity, , such that the time-dependent field is given as where is the angular frequency of time-harmonic oscillation and refers to the real part. We shall also use the convention whereby a three-dimensional (3-D) vector, for example, can be written as , where . This allows, for example, for the transverse and axial components in Cartesian space to be conveniently expressed independently. We shall denote unit vectors with a hat, i.e., , and where unit vectors are plotted (such as in Fig. 1), it is understood that no attention is paid to scale.2.1.Debye–Wolf IntegralThe field at position , in the vicinity of the focus of a lens of arbitrary numerical aperture (NA) and focal length , which satisfies the sine condition, can be found by evaluating the Debye–Wolf integral1–3 as The axial position [i.e., in Eq. (2) and Fig. 1] of the point of observation is measured from the location of the nominal focus. The lateral position of the point of observation [i.e., in Eq. (2) and Fig. 1] is measured from the optical axis. The point of observation is assumed to be located within a homogeneous region of refractive index . As shown in Fig. 1, upon multiplication by , the 2-D vector uniquely defines the position within the aperture from where a typical ray, assumed to propagate parallel to the optical axis, originates before being refracted by the lens. After refraction, the ray propagates parallel to . The field on a Gaussian reference sphere located in the exit pupil of the configuration shown in Fig. 1 is denoted as . The Gaussian reference sphere is shown in Fig. 1 and we note explicitly that corresponds to the field at the location on the reference sphere intersected by a line parallel to and intersecting the nominal focus. It is important to note that as the aperture is located in the back focal plane of the lens, the exit pupil appears at infinity as viewed from the space containing the focal region. It can be shown17 that this means that can, in fact, be found using geometrical optics.The first task in evaluating Eq. (2) is thus to find, using geometrical optics, for an arbitrary field in the back focal plane, which we denote by , where the superscript stands for the focal plane. It is assumed that has no axial component. We consider a ray propagating in the sample space, parallel to unit vector . We define as the angle that this ray makes with the optical axis, as shown in Fig. 1(a). This ray propagates in a particular meridional plane that makes an angle to the axis in the sample space as shown in Fig. 1(c). Suppose , where and are the unit vectors in the - and -directions, respectively. The generalized Jones matrix formalism18 provides a systematic way to find from . We start by defining a coordinate system transformation matrix , which transforms from a global Cartesian coordinate system to one where the and axes have been rotated by an angle around the axis. We also define a lens rotation matrix that rotates a vector by an angle of about the axis of the coordinate system. These matrix transformations are defined as18 It is then straightforward to show that and furthermore that18 We emphasize that when we write an equation involving one or more of , , and , it is implied that the values and are determined by such that Eq. (4) is satisfied.2.2.Stratified Medium Matrix Formalism for Plane WavesAs will be explained in further detail in this section, the Debye–Wolf integral expresses the electromagnetic focusing problem as a superposition of plane waves. In this section, we thus consider the interaction of plane waves with stratified media. 2.2.1.Preliminary considerationsWe follow a stratified medium matrix formalism as proposed by Azzam and Bashara.19 We consider a stratified medium as shown in Fig. 2 with interfaces located at , therefore possessing up to different media having refractive indices , the first and last layers being infinite half-spaces. We allow each refractive index potentially to be complex valued as and assume that both and are positive. We assume that the first medium is lossless and so is purely real valued. We decompose the electric field vector of waves propagating in the positive -directions onto transverse electric (TE) and transverse magnetic (TM) components with unit vectors and , respectively. We follow a similar convention for waves propagating in the negative -direction, except we interchange the superscript + with −. 2.2.2.Incident, reflected, and transmitted fieldsWe describe a plane wave, incident from the objective lens, propagating in the first medium as where we decompose into TE and TM components denoted by and , respectively. Following the convention of Azzam and Bashara,19 we define these vectors as where is the unit vector parallel with the optical axis. Without loss of generality, we henceforth assume that each plane wave component is either TE or TM polarized, thus allowing us to drop the TE and TM subscripts from Eq. (7). Consider such a plane wave incident upon the stratified medium shown in Fig. 2. The wave incident upon the first interface, at the first interface, takes the form , where we have added the superscript − to to indicate that this is the limiting value of the field as approaches from the negative -direction (see Fig. 2). Note that the -dependence of is expressed entirely by the scalar . Following the matrix formalism presented by Azzam and Bashara,19 we start by defining , , and as shown in Fig. 2, meaning that we have defined plane waves propagating in the positive and negative -directions, evaluated on both sides of the interface at . We can then define a matrix relationship between , , , and as19 where is implicitly defined and and are the Fresnel transmission and reflection coefficients, respectively, for the interface at . Next, we consider relating the fields within the second medium using a similar matrix formalism. To do this, consider plane wave components at the beginning and end of the second medium, i.e., , , , and . We can then define a relationship between their plane wave amplitudes, , , , and as19 where and is implicitly defined. It is then possible to construct a matrix that relates the fields at and as where is obtained by successive application of Eqs. (9) and (10) as We will assume in this work that , meaning that all light reaches the stratified medium after transmission through the objective lens, thus giving The only thing remaining is to find the vectors and , the polarization vectors in media 0 and , respectively. This is trivial in the TE case as they remain fixed in all media. In the TM case, the transverse component of the wave propagation vector in each medium (i.e., and ) is first calculated using the Snell’s law (see Sec. 2.2.5). The polarization vectors may then be found according to Eq. (8). The reflected field at the first interface and transmitted field at the last interface can thus be found as and , respectively. Note that the dependence upon is entirely embedded within the matrix formalism. In general, this step must be performed for the TE and TM cases independently and the results superimposed to obtain the total field. We have thus shown how to calculate the reflected and transmitted fields at the first and last interfaces, respectively, for an arbitrary stratified medium and an incident plane wave.2.2.3.Internal fields at internal interfacesSuppose it is desired to know the field at an interface in the interior of the stratified medium, say, for example, the field at , for some positive integer . As in Sec. 2.2.2, we can construct matrix relationships relating the field components at to the incident, reflected, and transmitted fields as where and This means that Eq. (14) can be used in conjunction with Eq. (13) to yield or Eq. (15) can be used in conjunction with Eq. (13) to yield2.2.4.Calculating fields between interfacesThere are many instances where it is necessary to evaluate the field away from an interface, either internal or external. One approach to this would be to first calculate the field due to the plane wave in question at an interface adjacent to the plane of constant in which the field must be evaluated. The field at the desired location can then be propagated from the adjacent interface using the plane wave’s propagation vector obtained using Snell’s law. There is a certain numerical efficiency to this approach, however, it requires keeping track of propagation vectors and plane wave phase terms. An alternative approach that is consistent with the matrix formalism discussed above is to introduce a “phantom” interface, which is an interface between two layers with identical refractive indices. In particular, suppose we wish to evaluate the field at some location that lies in the layer having refractive index . The field at can be found by following a similar procedure to that outlined in Sec. 2.2.3, except that we introduce an additional interface at . Then, we evaluate and instead of and , respectively, and apply either Eq. (18) or Eq. (19) to evaluate and . 2.2.5.Conducting media and total internal reflectionIn the case where media have a complex-valued refractive index, or when total internal reflection occurs, the wave vector of a plane wave becomes complex valued. This does not pose any difficulties as Snell’s law and the Fresnel reflection and transmission equations remain formally valid (see, for example19–21). Snell’s law can be generalized in the notation of this work by defining as the transverse component of a plane wave’s propagation vector in the ’th medium. The propagation vector in this medium will thus be given by . We can then write Snell’s law as where it is assumed that, as in Sec. 2.2.2, is the transverse component of the incident plane wave in a homogeneous medium of refractive index . The axial component of the plane wave in the ’th medium may then be found as .2.3.Stratified Medium Matrix Formalism for Focused IlluminationWe have shown in Sec. 2.1 how focusing of electromagnetic waves can be described as a superposition of plane waves using the Debye–Wolf integral. We have also shown in Sec. 2.2 how the field at any location within or outside a stratified medium can be evaluated for an arbitrary incident plane wave. Thus, by virtue of the linearity of the system, the forward or back propagating field anywhere within, or outside, the stratified medium may be found according to where and are evaluated using Eqs. (13), (18), or (19) as appropriate, depending on the value of . We find and using Eq. (8) combined with knowledge of in the appropriate medium as obtained from Eq. (20).2.4.Numerical ImplementationEquation (21) must be evaluated numerically. This is done by changing the variables of integration from to as defined by Eq. (4). Evaluation of the Jacobian of this transformation results in . The resulting integral is then evaluated using Gauss–Legendre numerical integration. The Gauss–Legendre technique allows for a definite integral to be evaluated as22 where are the zeros of the ’th-order Legendre polynomial, are the associated weights, and is the number of points used to evaluate the integration. A detailed explanation regarding the calculation of both the and is given by Davis et al.22 The 2-D integration of Eq. (21) is evaluated by integrating into and separately.3.Examples3.1.A-Scan Formation in Optical Coherence TomographyThe inverse scattering problem for OCT has been addressed previously in a 1-D setting using the theory of stratified media for plane waves.23 Although we do not consider the inverse scattering problem here, rather the forward problem, we demonstrate how selected problems with practical significance can be addressed using the model developed in this paper. We consider two stratified media as shown in Fig. 3. The stratified medium in Fig. 3(a) is used to demonstrate how the so-called confocal function of an OCT system can be calculated and that of Fig. 3(b) shows how dispersion of water changes the axial resolution of an OCT system. The simulation presented here is idealized and valid for both spectral domain and swept source OCT systems in the sense that a range of wavenumbers spanning to is uniformly sampled. We consider a fiber-based system with an objective having . The resulting spot size of the beam is determined by the combination of collimator, objective lens and mode field diameter of the single-mode fiber. This is integrated into our model by specifying the functions and as where is derived to be consistent with a previously characterized OCT system.24 The A-scan formation model is based upon a recently published full-wave model of OCT image formation.24 In brief, for both stratified media in Fig. 3, the incident beam that would be observed in the absence of a stratified medium is calculated just in front of the first refractive index discontinuity yielding . The field scattered back by the stratified medium is also calculated yielding . The modal coefficient describing the scattered field coupled back in to the fiber is given as25 a similar quantity, representing the field coupled back into the fiber from the reference arm is evaluated by repeating the above calculation for a single refractive index interface between air and a nonphysical material having refractive index . The OCT A-scan is then evaluated as where represents the effective shape of the OCT system’s spectrum. As already mentioned, full details of this model are available in a previous publication.24The first stratified medium considered shows how the confocal function of an OCT system can be calculated using a single A-scan. A refractive index distribution was modeled as shown in Fig. 3(a), where slabs, 100-nm wide, of material with refractive index of are spaced every in air. The narrow slabs act as weak reflectors that do not significantly perturb the illumination beam. The magnitude of the A-scan is plotted in Fig. 4(a), which shows a peak (plotted in blue) for each refractive index discontinuity in the stratified medium. These peaks lie on an envelope, plotted in red, which is evaluated as where the subscript denotes that is evaluated in free space. This example does not produce significant multiple reflections due to the low refractive index contrast between the reflecting slabs and the background refractive index .The second example of A-scan simulation involves a layer of water sandwiched between layers of glass as shown in Fig. 3(b). For this example, we assume that the glass layers are infinitely thick. The refractive index of the glass was assumed to be 1.525 at all wavelengths within the spectrum of the OCT system. The refractive index of water was found by interpolating the data published by Segelstein.26 An A-scan is plotted in Fig. 4(b) on a logarithmic vertical axis. The horizontal axis is scaled by the group refractive index of water obtained by numerically evaluating at the central wavelength of the OCT system’s spectrum, where is the phase refractive index. The A-scan of Fig. 4(b) shows three peaks separated by a physical distance of . The peak at corresponds to the reflection from the interface between glass and water at in Fig. 3(b). Similarly, the peak at is from the interface between water and glass at in Fig. 3(b). The third peak is due to light that has been multiply reflected and has undergone two complete round-trips of the water layer. The impact of dispersion due to water on the axial response of an OCT system can be studied by aligning each of these three peaks upon the same axial location as is shown in Fig. 4(c). These plots demonstrate that the peak corresponding to the first reflection does not exhibit dispersion as the A-scan is symmetric about . The two subsequent peaks that correspond to light propagating greater distances in water show how dispersion broadens axial response of an OCT system. 3.2.Cylindrical Vector Beams Focusing Through a Dielectric InterfaceThere are many examples in the literature of where the properties of unconventional beams focused through dielectric interfaces are studied (e.g., Refs. 4,7, and 27–29). In these, and many other related papers, much effort is devoted to evaluating Eq. (2) partially analytically, for a nonstandard form of such as may result when focusing, for example, a Bessel–Gauss beam. In such papers, is evaluated for a specific beam in the aperture of the focusing lens, , using the generalized Jones matrix formalism18 outlined in Sec. 2.1. In the absence of a stratified medium, the resulting value of is substituted into Eq. (2) before changing the variables of integration from to . The integration over is performed analytically using the identities for some integer , where is the ’th order Bessel function of the first kind. This results in an expression for the complex amplitude of the electric field in the vicinity of the focus in terms of the so-called “”-functions, which are functions of that are evaluated by performing a numerical integration for each unique combination of and . The significance of this procedure is that when a new beam type is to be modeled, it is generally necessary to derive the equations that determine the field from scratch. This process is made more difficult when a stratified medium is considered. In this example, we demonstrate first how a previously published result can be obtained, using the formalism presented in this paper, without the need to perform any analytic derivations. In particular, we consider the focusing of a radially polarized Bessel–Gauss beam originally considered by Biss and Brown27 who calculated how a beam of the form incident upon the aperture of a focusing lens, interacts with an air interface located at the focus. The lens was assumed to have a numerical aperture of 1.4 when focused into a medium of refractive index .We consider three different scenarios: focusing into a homogeneous medium with refractive index , through an interface at the focus between refractive indices and and, finally, through a wide layer of refractive index embedded in a refractive index . Scenario 1, depicted in Figs. 5(a) and 5(d), shows the profile of the beam in a homogeneous medium. The second scenario [Figs. 5(b) and 5(e)] replicates what was published by van de Nes et al.30 The final result [Figs. 5(c) and 5(f)] that we show as an example is where, instead of having an interface between and , there is a wide layer of refractive index . The subwavelength width of the layer allows for the evanescent waves to couple into the region causing the field after the final interface to more closely resemble that of the homogeneous case. This example is intended to demonstrate how complex focusing problems can be studied using this formalism in a much simpler way than was previously possible. 3.3.Simulation of Aberrations Arising in Confocal MicroscopyFor the final demonstration of how this formalism can be used, we consider an example from the literature of where a similar model has been applied.12 In particular, we study the illumination point spread function (PSF) that arises when performing confocal microscopy with a dry objective, with a sample embedded in water, covered by a glass cover slip as shown in Fig. 6. This particular example, considered previously by Török and Varga,12 employs a dry objective of numerical aperture 0.9, a cover slip of refractive index and a probe fluorescent molecule situated in water below the glass/water interface. The fluorescent molecule is merely a way of depicting the location at which we evaluate the focused field, we do not consider the process of fluorescence. When focusing from air through a glass cover slip, spherical aberration will occur if aberration compensation is not implemented. Although an analytic approach was derived by Török and Varga,12 we will use the stratified medium matrix formalism employed in this paper for this purpose. Suppose that we have an arbitrarily complex stratified medium that we wish to compensate for. We must calculate the phase aberration suffered by each plane wave components, and , appearing in Eq. (21), due to the stratified medium. This phase aberration can be found, as a function of , using the matrix that describes the entire stratified medium in Eq. (12). In particular, we know the ratio between the amplitude of the plane wave incident upon the stratified medium and that of the transmitted wave is given by Eq. (13) as We then use the matrix formalism again to calculate the ratio of amplitudes for the same incident plane wave when all refractive indices are set to that of air, i.e., the aberration free case. We denote this ratio . We thus define the aberration caused by the stratified medium by as where arg represents the argument of a complex number. Although is often real, it may be a complex value if evanescent waves are produced, which cannot be compensated for.We used this approach to calculate the rotationally symmetric aberration produced by the glass cover slip as plotted in Fig. 7(b). Another aspect of the model presented in this paper is the ability to model aberration functions with ease. In particular, it is straightforward to incorporate Zernike fringe polynomials into the model. For this work, we consider only the rotationally symmetric Zernike fringe polynomials defined as where . As these polynomials are orthogonal on the unit circle, we can decompose the aberration onto these functions using numerical integration. The result of this is shown in Fig. 7(a), which shows the contribution due to each individual polynomial. However, since we consider only the first six rotationally symmetric polynomials, aberration compensation described by a sum of these polynomials cannot match, exactly, the actual aberration. This error is plotted in Fig. 7(c) for the case where between three and six polynomials are considered.Now that the aberration introduced by the glass cover slip has been calculated, it can be compensated for by applying a phase modulation to the beam incident upon the aperture of the focusing lens as where, in this example, will be evaluated in two different ways. One way is by directly setting . The alternative way is by setting , for values of between 3 and 6, where is the weight associated with polynomial evaluated using numerical integration as .The unaberrated illumination PSF, evaluated as if the probe fluorescent molecule was in air, is plotted in Fig. 8(f). All other PSFs in Fig. 8 are normalized by the peak of this PSF. The PSFs in Figs. 8(a)–8(e) are calculated with the combination of glass cover slip and fluorescent molecule embedded in water. Aberration compensation is applied using an increasing number of Zernike fringe polynomials beginning with three in Fig. 8(a) (i.e., ) and going up to six in Fig. 8(d) (i.e., ). Significant spherical aberration is present in each of cases (a), (c), and (d) due to the fluorescent molecule being located in water rather than air. It is interesting, however, that case (b), where four Zernike fringe polynomials are considered, is substantially less aberrated. Analysis shows that the combination of aberration due to compensation, combined with propagation through the cover slip and water, results in less overall aberration in the water. Although not considered in this example, we emphasize that noncircularly symmetric aberration polynomials can be just as easily treated within this model. 4.Discussion and ConclusionsThe model presented in this paper will be useful across many fields of optics, particularly in optical imaging. The strength of the model is that its components that deal with the beam type, stratified medium and focusing are decoupled from one another. This has some important ramifications. One is that this means that incident beams and stratified media of arbitrary complexity can be modeled just as simply as a free-space focusing problem. Another consequence is that the implementation in software closely mirrors the theory itself, making the model amenable to use and modification by others. This approach is quite deliberate as it is intended that this model be available for use by researchers in the field, freeing them from the need to implement their own bespoke focusing code when a new, slightly different, application arises. Examples including OCT, cylindrical vector beams, and high numerical aperture focusing have been considered. Each of these examples is based upon a result that has already been published, plus some new results. These examples have been chosen to demonstrate the breadth of applications, which this model is able to address. There are several extensions that we intend to make to this model. Perhaps the most straightforward is the calculation of the magnetic fields, which is trivially evaluated as the magnetic field associated with a plane wave can be directly evaluated from the electric field amplitude and polarization, propagation vector and the refractive index. Consideration of a tilted stratified medium is also possible and of practical significance. Modeling of anisotropic layers within the stratified medium would also be of high practical importance. There are many improvements that can be made to the computational efficiency of the code. Little has been done to date in this regard as the emphasis of the code was on readability. However, there is scope for making the existing code more efficient. Furthermore, as the code is written in MATLAB it should be straightforward to execute it on a graphics processing unit using MATLAB’s parallel computing toolbox. AcknowledgmentsP.M. was supported by a Royal Society University Research Fellowship. P.M. acknowledges fruitful discussions with Professor Peter Török. ReferencesV. S. Ignatowsky,
“Diffraction by a lens of arbitrary aperture,”
Trans. Optical Inst. Petrograd, 1
(4), 1
–36
(1919). Google Scholar
R. Luneburg, Mathematical Theory of Optics, University of California Press, Berkeley and Los Angeles
(1966). Google Scholar
B. Richards and E. Wolf,
“Electromagnetic diffraction in optical systems ii. Structure of the image field in an aplanatic system,”
Proc. R. Soc. (London) A, 253 358
–379
(1959). https://doi.org/10.1098/rspa.1959.0200 Google Scholar
P. Török and P. R. T. Munro,
“The use of Gauss–Laguerre vector beams in STED microscopy,”
Opt. Express, 12 3605
–3617
(2004). https://doi.org/10.1364/OPEX.12.003605 OPEXFF 1094-4087 Google Scholar
K. S. Youngworth and T. G. Brown,
“Focusing of high numerical aperture cylindrical-vector beams,”
Opt. Express, 7 77
–87
(2000). https://doi.org/10.1364/OE.7.000077 OPEXFF 1094-4087 Google Scholar
P. Török, P. Varga and G. R. Booker,
“Electromagnetic diffraction of light focused through a planar interface between materials of mismatched refractive indices—structure of the electromagnetic field I,”
J. Opt. Soc. Am. A, 12
(10), 2136
–2144
(1995). https://doi.org/10.1364/JOSAA.12.002136 JOAOD6 0740-3232 Google Scholar
A. S. van de Nes et al.,
“Calculation of the vectorial field distribution in a stratified focal region of a high numerical aperture imaging system,”
Opt. Express, 12
(7), 1281
–1293
(2004). https://doi.org/10.1364/OPEX.12.001281 OPEXFF 1094-4087 Google Scholar
R. Kant,
“Vector diffraction problem of focusing a category 1 aberrated wavefront through a multilayered lossless medium,”
J. Mod. Opt., 51
(3), 343
–366
(2004). https://doi.org/10.1080/09500340408235528 JMOPEW 0950-0340 Google Scholar
O. Haeberlé,
“Focusing of light through a stratified medium: a practical approach for computing microscope point spread functions: Part II: confocal and multiphoton microscopy,”
Opt. Commun., 235
(1–3), 1
–10
(2004). https://doi.org/10.1016/j.optcom.2004.02.068 OPCOB8 0030-4018 Google Scholar
P. R. T. Munro, A. Curatolo and D. D. Sampson,
“Full wave model of image formation in optical coherence tomography applicable to general samples,”
Opt. Express, 23 2541
–2556
(2015). https://doi.org/10.1364/OE.23.002541 OPEXFF 1094-4087 Google Scholar
E. Z. Y. Zhang and P. C. Beard,
“Ultrahigh-sensitivity wideband Fabry–Perot ultrasound sensors as an alternative to piezoelectric PVDF transducers for biomedical photoacoustic detection,”
Proc SPIE, 5320 222
–229
(2004). https://doi.org/10.1117/12.531160 Google Scholar
P. Török and P. Varga,
“Electromagnetic diffraction of light focused through a stratified medium,”
Appl. Opt., 36
(11), 2305
–2312
(1997). https://doi.org/10.1364/AO.36.002305 APOPAI 0003-6935 Google Scholar
I. Gregor and J. Enderlein,
“Focusing astigmatic Gaussian beams through optical systems with a high numerical aperture,”
Opt. Lett., 30 2527
–2529
(2005). https://doi.org/10.1364/OL.30.002527 OPLEDP 0146-9592 Google Scholar
D. G. Flagello, T. Milster and A. E. Rosenbluth,
“Theory of high-NA imaging in homogeneous thin films,”
J. Opt. Soc. Am. A, 13 53
–64
(1996). https://doi.org/10.1364/JOSAA.13.000053 JOAOD6 0740-3232 Google Scholar
D. G. Flagello and T. Milster,
“3D modelling of high numerical aperture imaging in thin films,”
Proc. SPIE, 1625 246
–261
(1992). https://doi.org/10.1117/12.58952 PSISDG 0277-786X Google Scholar
P. R. T. Munro,
“Multi-layer vectorial focusing,”
(2018) http://prtmunro.net September ). 2018). Google Scholar
P. Török,
“An imaging theory for advanced, high numerical aperture optical microscopes,”
(2003). Google Scholar
P. Higdon, P. Török and T. Wilson,
“Imaging properties of high aperture multiphoton fluorescence scanning microscopes,”
J. Microsc., 193
(2), 127
–141
(1999). https://doi.org/10.1046/j.1365-2818.1999.00448.x JMICAR 0022-2720 Google Scholar
R. Azzam and N. Bashara, Ellipsometry and Polarized Light, 3rd ed.Elsevier, North Holland, Amsterdam
(1992). Google Scholar
C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles, Wiley Interscience, New York
(1983). Google Scholar
J. Stratton, Electromagnetic Theory, McGraw-Hill, New York
(1941). Google Scholar
P. R. Philip, J. Davis and W. Rheinbolt, Methods of Numerical Integration, Academic Press Limited(1984). Google Scholar
O. Bruno and J. Chaubell,
“One-dimensional inverse scattering problem for optical coherence tomography,”
Inverse Prob., 21 499
–524
(2005). https://doi.org/10.1088/0266-5611/21/2/006 INPEEY 0266-5611 Google Scholar
P. R. T. Munro,
“Three-dimensional full wave model of image formation in optical coherence tomography,”
Opt. Express, 24 27016
–27031
(2016). https://doi.org/10.1364/OE.24.027016 OPEXFF 1094-4087 Google Scholar
P. R. T. Munro,
“Exploiting data redundancy in computational optical imaging,”
Opt. Express, 23 30603
–30617
(2015). https://doi.org/10.1364/OE.23.030603 OPEXFF 1094-4087 Google Scholar
D. J. Segelstein,
“The complex refractive index of water,”
University of Missouri–Kansas City,
(1981). Google Scholar
D. P. Biss and T. G. Brown,
“Cylindrical vector beam focusing through a dielectric interface,”
Opt. Express, 9
(10), 490
–497
(2001). https://doi.org/10.1364/OE.9.000490 OPEXFF 1094-4087 Google Scholar
K. Hu, Z. Chen and J. Pu,
“Generation of super-length optical needle by focusing hybridly polarized vector beams through a dielectric interface,”
Opt. Lett., 37 3303
–3305
(2012). https://doi.org/10.1364/OL.37.003303 OPLEDP 0146-9592 Google Scholar
D. Khoptyar et al.,
“Tight focusing of laser beams in a -microcavity,”
Opt. Express, 16 9907
–9917
(2008). https://doi.org/10.1364/OE.16.009907 OPEXFF 1094-4087 Google Scholar
A. van de Nes et al.,
“Cylindrical vector beam focusing through a dielectric interface: comment,”
Opt. Express, 12 967
–969
(2004). https://doi.org/10.1364/OPEX.12.000967 OPEXFF 1094-4087 Google Scholar
BiographyPeter R. T. Munro received his BSc (computer science) and BEng (electrical engineering) degrees from the University of Western Australia in 2000. He completed his PhD in the Department of Physics at Imperial College in 2006. Currently, he holds a Royal Society University research fellowship focused on performing quantitative imaging using a combination of x-ray, optical, and photoacoustic techniques. |