The action of a nanoscopic spherically symmetric refractive index profile on a focused Gaussian beam may easily be envisaged as the action of a phase-modifying element, i.e. a lens: Rays traversing the inhomogeneous refractive index field n(r) collect an additional phase along their trajectory which advances or retards their phase with respect to the unperturbed ray. This lens-like action has long been understood as being the mechanism behind the signal of thin sample photothermal absorption measurements [Appl. Opt. 34, 41–50 (1995)], [Jpn. J. Appl. Phys. 45, 7141–7151 (2006)], where a cylindrical symmetry and a different lengthscale is present. In photothermal single (nano-)particle microscopy, however, a complicated, though prediction-wise limited, electrodynamic scattering treatment was established [Phys. Rev. B 73, 045424 (2006)] during the emergence of this new technique. Our recent study [ACS Nano, DOI: 10.1021/nn300181h] extended this approach into a full ab-initio model and showed for the first time that the mechanism behind the signal, despite its nanoscopic origin, is also the lens-like action of the induced refractive index profile only hidden in the complicated guise of the theoretical generalized Mie-like framework. The diffraction model proposed here yields succinct analytical expressions for the axial photothermal signal shape and magnitude and its angular distribution, all showing the clear lens-signature. It is further demonstrated, that the Gouy-phase of a Gaussian beam does not contribute to the relative photothermal signal in forward direction, a fact which is not easily evident from the more rigorous EM treatment. The presented model may thus be used to estimate the signal shape and magnitude in photothermal single particle microscopy.
© 2012 Optical Society of America
CorrectionsMarkus Selmke, Marco Braun, and Frank Cichos, "Nano-lens diffraction around a single heated nano particle: errata," Opt. Express 21, 25344-25345 (2013)
Photothermal lens spectroscopy (PLS) has become a valuable tool in the study of solids and liquids [1–3]. Recent publications include the study of non-linear effects  and nanoparticles in solution [3, 5]. Many authors have focused on the theoretical description of the thin-sample slab geometry which is common in such macroscopic lensing experiments, often providing numerical equations which may be used to obtain absorption coefficients [6, 7]. In all these models the thermal lens induced originates from the absorbed power of a heating laser which constitutes a spatially extended cylindrically symmetric heat source in the heat equation. The solution obtained is then used to study the effect on the propagation of a probing laser beam. In thermal lens spectroscopy (TL) the probing beam is coaxial (possibly offset) with the heating beam [7, 8] while in beam deflection spectroscopy [9, 10] (BDS) the probing and the heating beams are aligned perpendicularly to each other. Both methods have in principle the same sensitivity .
In contrast to these macroscopic techniques, a new microscopic approach has been developed in the 1990’s by Harada and Kitamori [12, 13] termed photothermal lens microscopy. The developement of single particle photothermal microscopy [14–16] followed and detects a very different kind of thermal lens with a modified standard confocal fluorescence microscope: Instead of an axially symmetric refractive index profile a spherically symmetric profile n(r), created by the point-like heat-source which is the absorbing nano-particle, is probed. Also, instead of a profile that decays on the length-scale of the heating beam focus, in PT single particle microscopy a lens is probed which decays to half its value on the length-scale of the nanoparticle. On the other hand, the profile extends infinitely as 1/r and thus misses a characteristic length-scale. This is the reason, why the description presented in this paper and within the recent electrodynamic (EM) study  provide no evidence for a role-play of the Gouy-phase which is otherwise important for the probing of small scatterers as shown by Hwang and Moerner . While models for spherical absorbers have been put forward [19, 20], these were numerical in nature and targeted for the μm-size regime. The theoretical description of the nanoscopic photothermal lens has been first given by Berciaud et al.  in a scattering treatment considering an extinction mechanism. Our recent ab-initio theoretical description of the electrodynamic problem has shown, however, that instead a simple lensing mechanism is responsible for the photothermal signal of nanoscopic absorbers showing a clear angular thermal diffraction signature and a double-lobe lens signature in axial scans. We have also demonstrated the importance of both the interference- (extinction) and scatter-term in forward detection appearing in the general situation which complicated an intuitive understanding within such a scattering framework.
It is the aim of this paper to show that both worlds, the vast literature on macroscopic thermal lens spectroscopy and the recent emerging tool of single particle thermal lens microscopy are very similar. The diffraction picture which has been a successful tool in the first domain will be shown to yield analytical and easily tractable expressions in the nano-scopic domain. They will allow for a quantitative assessment of absorption cross-sections of single nano-object based on standard photothermal measurements and will be able to explain the main phenomena of photothermal microscopy qualitatively as well as quantitatively, while providing an intuitive picture of the working mechanism. The quality of the simple model is checked against the more elaborate electromagnetic model within the extended scattering description. Axial scans and angular patterns of single heated nano-particles will be described and compared. The angular diffraction pattern will be shown to explain the signal inversion observed for the first time for a single nano-particle upon the introduction of an inverse aperture in the detection path.
2. The thermal lens
The lens to be considered in single particle generated nano-lens experiments such as photothermal microscopy originates from the absorption of electromagnetic power provided by a focused laser beam by a single nano-particle. The absorbing particle can be treated as a point-like heat-source, yielding the steady-state temperature profile T(r) = T0 + ΔTR/r which decays with the inverse distance r from the particle of radius R to the unperturbed ambient temperature T0. In the case of modulated heating, as utilized in the lock-in approach common to photothermal single particle microscopy, the assumption remains valid as long as the modulation frequencies used remain below Ωm ≈ 1MHz, depending on the thermal conductivity κ and the heat capacity C per unit volume of the medium [17, 21]. Only in case of higher modulation frequencies where the thermal length scale Rth = [2κ/(CΩm)]1/2 becomes comparable to the diffraction limit of the focusing optics the detailed form of the time-dependent temperature profile needs to be considered. As the signal of photothermal single particle microscopy has been shown to decrease considerably above this limiting frequency , we postpone a detailed discussion of this issue to a different publication.
By the temperature dependence of the surrounding mediums’ refractive index n(T), quanti-fied by its thermorefractive coefficient (dn/dT), a corresponding refractive index profile n(r) is established as indicated in Fig. 1(a):Fig. 1(b)). Together with the thermal conductivity κ and the radius R of the particle this controls the induced temperature and thus the contrast Δn of the lens Eq. (1) decays to half its maximum perturbation Δn at a distance of r = 2R, making it a nanoscopic object. Nonetheless, it has infinite extent. It will only be limited by the probing beam which will for most cases be a focused beam of diffraction-limited extent. As Δn will quantify the photothermal single particle signal, the knowledge of the heating beam intensity and the thermal properties of the embedding bulk material will then allow the determination of the particles’ absorption cross-section σabs. It is thus necessary to obtain expressions for the photothermal single particle signal from the induced refractive contrast Δn. This is the purpose of the next section.
3. The diffraction integral
The Frauenhofer diffraction integral for circular apertures in the Fresnel-grade approximation connects the complex field-amplitude of the probe beam Ua (ρ) in the aperture plane at z = 0 to the complex field amplitude U (r,z) in the image plane at a distance z (see Fig. 1(b)) [22,23]. A further factor is included which represents the collected phase Δχ for each component-wave as a result of the thermal lens (see Fig. 1(a)) .
To simplify matters, we will consider an axially symmetric scenario only, i.e. the refractive index profile will be symmetric with respect to the optical axis and the steady-state or low-frequency limit (as defined above) of the temperature profile 1/r will be assumed. This corresponds to the case of a heated nano-particle positioned along the optical axis, see Fig. 1(a).
The phase advance Δχ (ρ) for a single ray passing the lens immersed in a sample slab of thickness 2L at a distance ρ from the optical axis can be approximated via a straight ray calculation already utilized for collimated beams in reference  and similarly in reference :
The field amplitude in the aperture plane will be taken to be the complex field of the probing Gaussian beam  having a focus displaced by zp from the center of the heated particle:
To put photothermal single particle microscopy on a quantitative footing, we have introduced the relative PT signal Φ . It is the change in intensity I ∝ |U|2 within the image-plane relative to the constant much larger background of the unperturbed field:
Through direct computation it can be shown that for small particles R < 100 nm the relative signal is not affected much by the geometrical diffraction of the inverse aperture disc of radius R. The diffraction-integral may be extended to R = 0 without an appreciable change in the result. We find the following analytic expression on the optical axis (r = 0), corresponding to the intensity-change detected in a closed aperture scenario:Fig. 2(a)) obtained with this expression may have been anticipated from the close similarity of the situation presented here and the thermal lens model used in thin sample slab geometries [6, 7, 27]. In these cases, the parabolic refractive index profile approximation for example yields an optimal probing beam offset relative to the sample of one confocal distance, zp = zR. Here, the maximum signal is obtained when the heated particle is offset by about zp = 0.7zR. For the off-axis image-plane signal, where the particle is still assumed to be positioned on the optical axis, the radial coordinate r in the image-plane was transformed to an angular coordinate via tan(θ) = r/z. The expression for the signal (Eq. (6)) detected under an angle θ for r ≠ 0 then reads: Eq. (7) since tan(θ) → 0 and 1F1 (c, 1, 0) = 1 for any number c. For a probing beam positioned behind the lens (zp < 0), the angular pattern described by this equation shows a peak towards the center and an annular dip at larger angles (see Fig. 2(b)). For the case of a probing beam being positioned in front of the lens (zp > 0), the angular pattern simply changes its sign, relative to the previous scenario. Overall, the energy of the beam is only redistributed by the action of the lens, i.e. an integration of Φ (θ,zp) sin (θ) from 0 to π/2 gives zero. This corresponds to a geometrical optics picture where bundles of rays would be refracted to form a changed intensity pattern behind the lens [24, 28]. From this result we may, similar to the treatment given in , compute the total detected photothermal signal if a finite angular detection domain, determined by the numerical aperture (NA) of the collecting microscope objective, is used. In this case it is not the modulation of the intensity detected on-axis or under a certain angle θ, but rather one needs to integrate the angular spectrum given in Eq. (8) over the angular detection domain: Eq. (8) instead of the true Gaussian intensity. Eq. (7) and Eq. (9) shows a change in sign and a dispersion-like signal if Δn is changed according to the axial shifting of the absorbing particle with respect to the local heating laser. This is just what one would expect for the collected signal behind a (modifiable) diverging lens if probed by a focused beam and the lens were moved along the beam axis across the focus while its focal length were to change in such a way that it has its smallest |f| close to the probe beam focus and gets weaker, i.e. |f| larger, for large offsets. A further consequence of the given angular pattern of the photothermal signal is the possibility of an inversion of the photothermal signal. By using an annular aperture, only the annular region described by Eq. (8) may be collected, and the total collected signal as described by Eq. (9) will be of opposite sign as compared to the use of a usual circular aperture (Objective NA). Indeed, the reversal of sign in the photothermal signal could both be observed experimentally for R = 60 nm AuNPs and confirmed in the calculations using the generalized Lorenz-Mie theory (see section 4).
Equations (7), (8) and (9) present the main results of this paper and give the background normalized intensity change Φ = ΔI/I detected either on the optical axis or detected under a finite angle θ with respect to the optical axis upon the introduction of the lens n(r), Eq. (1), on the optical axis. Hereby, the particle/lens was assumed to be displaced by a distance zp with respect to the probing beam-waist position. In this approach the Gouy-phase terms, while present in the fields in the aperture plane, cancel each other in the relative photothermal signal since they do not depend on the integration variable ρ. The found expressions which are exemplarily visualized in the plots of Fig. 2 allow the study of (finite NA) zp-scans in photothermal microscopy setups similar to the thin sample slab studies in [4, 29–33]. To this end, the collection angle depends on the collecting microscope objectives’ numerical aperture through θmax = arcsin (NA/n0), while the illuminating objective determines the beam waist(s).
4. Comparison to rigorous vectorial electromagnetic (EM) model
To compare the found results Eq. (7) and (8) to a more rigorous solution of the problem, a full vectorial electromagnetic treatment will be used in the following. Therefore, the problem at hand may be expressed as the scattering of a shaped incident probe field interacting with a multilayered scatterer. The scatterer described by Eq. (1) may be viewed as an unbounded gradient refractive index lens (GRIN). The scattering of such an object has been studied in the literature and a publicly available C-code is attainable through reference  providing the scattering coefficients of the GRIN aN+1 and bN+1 when the scatterer is discretized into N concentric spherical shells as shown in Fig. 3(b). The appropriate scattering framework is the generalized Lorenz-Mie theory (GLMT), which is applicable to any spherically symmetric scatterer  positioned in an arbitrarily shaped beam [36–38]. We have therefore adopted and modified it to yield an analytical expression for the power 𝒫d of the total electromagnetic field (Et, Ht) detected in the far-field under a scattering-angle θ. This treatment has been verified by exhaustive comparison to single gold nanoparticle scattering and photothermal microscopy and is described in more detail in our reference . The quantity of interest is the projected total fields’ Poynting-vector St under a certain angle θ with respect to the optical axis (Fig. 3(a)):
Within the GLMT formalism the total field is mathematically decomposed into series expressions of the incidence field Ei and Hi with complex expansion coefficients gn, and an outgoing scattered field Es and Hs described by complex scatter coefficients an and bn. Thereby, the total field consists of the incident field and the scattered field, Et = Ei + Es and Ht = Hi + Hs. The scattered far-field is commonly expressed via the scattering amplitudes S1,2 through , , , . These amplitudes are:39]: Eq. (11), with the integration over the azimuthal angle ϕ carried out in the far-field can now be expressed as differential cross-sections in any forward/backward polar angle θ : d𝒫ext (θ) = −dσext (θ)I0, d𝒫sca (θ) = dσsca (θ)I0 and d𝒫inc (θ) = dσinc (θ)I0, where the Gaussian beam focus intensity was used: 40]. The above expressions may now be used to obtain the angular spectrum of the difference signal, i.e. the PT signal analogously defined to Eq. (6).
To obtain the relative PT signal collected with a finite detection aperture, the integrated differential cross-sections dσsca,ext have to be used instead of the differential ones. This corresponds then to Eq. (8) in the diffraction model and to Eq. (9) when integrated accordingly. The predictions of both models are displayed in Fig. 4. The best agreement is found in the low-focusing regime (ω0 ≫ λ). In the strong focusing case, the GLMT may be expected to deviate since the beam shape coefficients used here rely on the low focusing expansion of the paraxial field, while within the diffraction treatment the straight ray phase advance approximation, Eq. (4), becomes less accurate. Both frameworks agree for typical experimental parameters within a factor of the order of unity. A direct scaling-free comparison of the on-axis diffraction formula Eq. (7) and the Gaussian GLMT prediction with Eq. (21) is shown in Fig. 4(b), where the solid thin black line corresponds to Eq. (7) and the red dashed line corresponds to Eq. (21). The qualitative agreement is illustrated by further zp-scan examples with finite numerical collection apertures (see Fig. 4(c), NA = 0.3 and NA = 0.75), i.e. Eq. (9), scaled by 1.5 to match the GLMT predictions. The angular pattern in Fig. 4(a) of the photothermal signal has two consequences: The relative change in intensity as compared to the unperturbed beam, i.e. the relative photothermal signal, is maximal, if the detection is on-axis or in a small angular detection domain around the forward direction (see Fig. 4(b) and 4(c)). Further, the relative as well as the absolute PT signal will decrease if the angular collection domain extends across the zero-crossing of the angular pattern. Also, both frameworks predict a signal inversion for the collection of an annular region of detection angles (Fig. 4(d)), for example for [θmin, θmax] = [21°, 31°]. To test this prediction, a central beam-stop experiment has been carried out which is detailed in the following section.
Notably, the diffraction framework Eq. (9) also predicts the dependence of the spatial photothermal signal shape on the axial displacement of the heating and the detection beam focus Δzf as observed experimentally in our recent study . Once again, the calculations have been carried out in both frameworks. In Fig. 5, the GLMT calculations are shown by the solid lines, while the diffration framework corresponds to the dashed-solid lines. The numerical aperture of the detection lens for these calculations is NA = 0.75, resulting in a maximum collection angle of θmax = 30°. The focal displacement Δzf that is included in Eq. (2) has been varied from −3zR,d to 3zR,d. The effect in the photothermal signal can be seen in Fig. 5(a). Depending on Δzf, either the positive signal is enhanced while the negative one is decreased (for a positive Δzf) or vice versa. However, the position of the zero-intersection, when the detection beam is focused to the center of the thermal lens, is not influenced, which is covered by both frameworks. By scaling the curves obtained by diffraction theory again by a constant factor of 1.6, the overall shape matches very well with the curves predicted by GLMT. Also, the dependence of the maximum and the minimum signal on the focal displacement Δzf shows a reasonable agreement (see Fig. 5(b)), when the scaling factor of order unity is included. The influence of Δzf on the maxima and minima position is depicted in Fig. 5(c). The axial distance between maximum and minimum is the shortest for Δzf = 0. In this region the predictions of both theories are nearly identical. Only for a large focal displacement Δzf > 1.5zR,d, the deviation becomes noticable.
In summary, the comparison of the simple and intuitive diffraction-picture, as put forward in the preceding section, and the rigorous vectorial treatment of a Gaussian beam scattered by the GRIN lens shows a near-perfect qualitative agreement, while a quantitative agreement within a factor of order unity is found.
5. Signal inversion of an axial single particle scan
Having demonstrated the equivalence of the results obtainable within the generalized Lorenz-Mie framework and in the simple diffraction model, the following part of the paper will show a direct consequence of the above findings. Whereas the collection of the probe-beam in an angular domain around the forward-direction will show a dispersion-like lens-signature as displayed in Fig. 2, the collection of an annular region can invert the detected signal. Indeed, the measurement of the photothermal signal of a single gold nanoparticle of R = 30 nm radius shows this effect. To this end, a central beam-stop was introduced into the detection beam path of a photothermal microscope as described in detail in . The collimated beam of diameter Do = 10 mm was thereby reduced to an annular ring of Di = 9 mm inner diameter corresponding to an angle of θmin = arcsin (NA Di/ (Don0)) = 27° as given by the Abbe sine condition and the used numerical aperture of the detection objective. The illumination objective was a high-NA oil immersion objective and the resulting probe and heating beam waists were ω0,d = 281 nm and ω0,h = 233 nm, respectively. The beams were offset in axial direction by Δzf = 350 nm to achieve a symmetric signal configuration with the aberrated beams. The observed z- and xz-scans (Fig. 6(a)) of the signal measured correspond very well to the calculated scans (Fig. 6(b)) and clearly show the inversion of the signal. In addition to a pure change in sign of the signal, the zero-crossing shifts. Although these details are only reconcilable with the predictions of the exact GLMT description (see supplement of ), which includes the details of the aberrated beams, the basic principle and the prediction of the effect is understood based on the angular diffraction pattern described in section 3.
While a close relation of Gaussian beam diffraction and scattering has been shown by J. Lock et al.  for the case of spherical dielectric particles, the more complicated scattering approach was thus far the only theoretical approach for the emerging technique of single particle photothermal microscopy. Although this ansatz may be used, the signature of a simple lensing mechanism (e.g. in zp-scans) suggests a more intuitive model. Employing the scalar diffraction formalism common in thin sample slab absorption spectroscopy we have demonstrated that the signal obtained in photothermal single (nano-)particle microscopy can indeed be understood as the signal of a phase-modifying element, i.e. a lens, despite the microscopic origin of the latter. In contrast to the former, the spherical symmetry of the heated single-particle lens allows a direct analytical evaluation of the relative photothermal signal. The analytical model presented here could be verified in shape and absolute value when compared to experiments and a full electromagnetic vectorial treatment as published by the authors earlier . Further, it could be shown, that the Gouy-phase of a Gaussian detection beam does not contribute to the relative photothermal signal in forward direction. The understanding of the mechanism and angular distribution of the photothermal signal was then shown to explain the observed signal inversion upon the collection of an angular domain corresponding to the outer angles only. This introduces a simple and intuitive model for single particle absorption measurements and allows the quantitative assessment of nano-particle absorption cross sections.
Financial support by the DFG research unit 877 and the graduate school BuildMoNa as well as funding by the European Union and the Free State of Saxony is acknowledged.
References and links
1. X. S. Xie, S. J. Lu, W. Min, S. S. Chong, and G. R. Holtom, “Label-free imaging of heme proteins with two-photon excited photothermal lens microscopy,” Appl. Phys. Lett. 96, 113701 (2010). [CrossRef]
2. S. Sinha, A. Ray, and K. Dasgupta, “Solvent dependent nonlinear refraction in organic dye solution,” J. Appl. Phys. 87, 3222–3226 (2000). [CrossRef]
3. A. V. Brusnichkin, D. A. Nedosekin, M. A. Proskurnin, and V. P. Zharov, “Photothermal lens detection of gold nanoparticles: Theory and experiments,” Appl. Spectrosc. 61, 1191–1201 (2007). [CrossRef] [PubMed]
4. G. Battaglin, P. Calvelli, E. Cattaruzza, F. Gonella, R. Polloni, G. Mattei, and P. Mazzoldi, “Z-scan study on the nonlinear refractive index of copper nanocluster composite silica glass,” Appl. Phys. Lett. 78, 3953–3955 (2001). [CrossRef]
7. J. Moreau and V. Loriette, “Confocal dual-beam thermal-lens microscope: Model and experimental results,” Jpn. J. Appl. Phys. 45, 7141–7151 (2006). [CrossRef]
9. M. Harada, T. Kitamori, and T. Sawada, “Phase signal of optical beam deflection from single microparticles -theory and experiment,” J. Appl. Phys. 73, 2264–2271 (1993). [CrossRef]
10. J. Q. Wu, T. Kitamori, and T. Sawada, “Theory of optical beam deflection for single microparticles,” J. of Appl. Phys. 69, 7015–7020 (1991). [CrossRef]
12. M. Harada, K. Iwamotok, T. Kitamori, and T. Sawada, “Photothermal microscopy with excitation and probe beams coaxial under the microscope and its application to microparticle ananlysis,” Anal. Chem. 65, 2938–2940 (1993). [CrossRef]
13. K. Uchiyama, A. Hibara, H. Kimura, T. Sawada, and T. Kitamori, “Thermal lens microscope,” Jpn. J. Appl. Phys. 39, 5316–5322 (2000). [CrossRef]
15. S. Berciaud, L. Cognet, G. Blab, and B. Lounis, “Photothermal heterodyne imaging of individual nonfluorescent nanoclusters and nanocrystals,” Phys. Rev. Lett. 93, 257402 (2004). [CrossRef]
18. J. Hwang and W. E. Moerner, “Interferometry of a single nanoparticle using the gouy phase of a focused laser beam,” Opt. Commun. 280, 487–491 (2007). [CrossRef]
19. G. C. K. Chen, M. Andika, and S. Vasudevan, “Excitation temporal pulse shape and probe beam size effect on pulsed photothermal lens of single particle,” J. Opt. Soc. Am. B 27, 796–805 (2010). [CrossRef]
20. M. Harada, M. Shibata, T. Kitamori, and T. Sawada, “Application of coaxial beam photothermal microscopy to the analysis of a single biological cell in water,” Anal. Chim. Acta. 299, 343–347 (1995). [CrossRef]
21. S. Berciaud, D. Lasne, G. Blab, L. Cognet, and B. Lounis, “Photothermal heterodyne imaging of individual metallic nanoparticles: Theory versus experiment,” Phys. Rev. B 73, 045424 (2006). [CrossRef]
22. Z. L. Horvath and Z. Bor, “Focusing of truncated gaussian beams,” Opt. Commun. 222, 51–68 (2003). [CrossRef]
23. S. Teng, T. Zhou, and C. Cheng, “Fresnel diffraction of truncated gaussian beam,” Optik 118, 435–439 (2007). [CrossRef]
24. A. A. Vigasin, “Diffraction of light by absorbing inclusions in solids,” Kvant. Elektron. (Moscow) [Sov. J. Quantum Electron.] 4, 662–666 (1977).
25. G. Baffou, P. Bon, J. Savatier, J. Polleux, M. Zhu, M. Merlin, H. Rigneault, and S. Monneret, “Thermal Imaging of Nanostructures by Quantitative Optical Phase Analysis,” ACS Nano , DOI: [CrossRef] [PubMed]
26. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics (John Wiley and Sons, Inc., 1991). [CrossRef]
28. M. Selmke, M. Braun, and F. Cichos, “Photonic Rutherford Scattering,” in preparation (2012).
29. R. Escalona, “Comparative study between interferometric and z-scan techniques for thermal lensing characterization,” Opt. Commun. 281, 1323–1330 (2008). [CrossRef]
30. S. M. Mian, S. B. McGee, and N. Melikechi, “Experimental and theoretical investigation of thermal lensing effects in mode-locked femtosecond z-scan experiments,” Opt. Commun. 207, 339–345 (2002). [CrossRef]
32. A. Gnoli, A. M. Paoletti, G. Pennesi, G. Rossi, and M. Righini, “High-accuracy z-scan measurements of the optical nonlinearity of bis-phthalocyanines,” J. Porphyrins Phthalocyanines 11, 481–486 (2007). [CrossRef]
33. R. Polloni, B. F. Scremin, P. Calvelli, E. Cattaruzza, G. Battaglin, and G. Mattei, “Metal nanoparticles-silica composites: Z-scan determination of non-linear refractive index,” J. Non-Cryst. Solids 322, 300–305 (2003). [CrossRef]
34. O. Pena and U. Pal, “Scattering of electromagnetic radiation by a multilayered sphere,” Comput. Phys. Commun. 180, 2348–2354 (2009). [CrossRef]
36. G. Gouesbet, G. Grehan, and B. Maheu, “Scattering of a gaussian-beam by a mie scatter center using a bromwich formalism,” J. Optics-Nouvelle Revue D Optique 16, 83–93 (1985).
37. G. Gouesbet, B. Maheu, and G. Grehan, “Light-scattering from a sphere arbitrarily located in a gaussian-beam, using a bromwich formulation,” J. Opt. Soc. Am. A 5, 1427–1443 (1988). [CrossRef]
38. G. Gouesbet, J. Lock, and G. Grehan, “Generalized lorenz-mie theories and description of electromagnetic arbitrary shaped beams: Localized approximations and localized beam models, a review,” J. Quant. Spectrosc. Radiat. Transfer 112, 1–27 (2011). [CrossRef]
40. G. H. Meeten, “Computation of s1–s2 in mie scattering-theory,” J. Phys. D: Appl. Phys. 17, L89–L91 (1984). [CrossRef]
41. J. A. Lock and E. A. Hovenac, “Diffraction of a gaussian-beam by a spherical obstacle,” Am. J. Phys. 61, 698–707 (1993). [CrossRef]