A framework for far-field infrared absorption microscopy beyond the diffraction limit

A framework is proposed for infrared (IR) absorption microscopy in the far-field with a spatial resolution below the diffraction limit. The sub-diffraction resolution is achieved by pumping a transient contrast in the population of a selected vibrational mode with IR pulses that exhibit alternating central minima and maxima, and by probing the corresponding absorbance at the same wavelength with adequately delayed Gaussian pulses. Simulations have been carried out on the basis of empirical parameters emulating patterned thin films of octadecyltrichlorosilane and a resolution of 250 nm was found when probing the CH2 stretches at 3.5 μm with pump energies less than ten times the vibrational saturation threshold. ©2012 Optical Society of America OCIS codes: (170.0180) Microscopy; (300.1030) Absorption; (300.6340) Spectroscopy, infrared; (300.6500) Spectroscopy, time-resolved; (300.6360) Spectroscopy, laser; (300.6390) Spectroscopy, molecular. References and links 1. W. Min, C. W. Freudiger, S. Lu, and X. S. Xie, “Coherent nonlinear optical imaging: beyond fluorescence microscopy,” Annu. Rev. Phys. Chem. 62(1), 507–530 (2011). 2. A. Volkmer, L. D. Book, and X. S. Xie, “Time-resolved coherent anti-Stokes Raman scattering microscopy: imaging based on Raman free induction decay,” Appl. Phys. Lett. 80(9), 1505–1507 (2002). 3. M. Jurna, J. P. Korterik, C. Otto, J. L. Herek, and H. L. Offerhaus, “Background free CARS imaging by phase sensitive heterodyne CARS,” Opt. Express 16(20), 15863–15869 (2008). 4. C. L. Evans, E. O. Potma, M. Puoris’haag, D. Côté, C. P. Lin, and X. S. Xie, “Chemical imaging of tissue in vivo with video-rate coherent anti-Stokes Raman scattering microscopy,” Proc. Natl. Acad. Sci. U.S.A. 102(46), 16807–16812 (2005). 5. M. Balu, G. Liu, Z. Chen, B. J. Tromberg, and E. O. Potma, “Fiber delivered probe for efficient CARS imaging of tissues,” Opt. Express 18(3), 2380–2388 (2010). 6. I. Toytman, K. Cohn, T. Smith, D. Simanovskii, and D. Palanker, “Non-scanning CARS microscopy using widefield geometry,” Proc. SPIE 6442, 64420D, 64420D-7 (2007). 7. C. W. Freudiger, W. Min, B. G. Saar, S. Lu, G. R. Holtom, C. He, J. C. Tsai, J. X. Kang, and X. S. Xie, “Labelfree biomedical imaging with high sensitivity by stimulated Raman scattering microscopy,” Science 322(5909), 1857–1861 (2008). 8. G. Romero, E. Rojas, I. Estrela-Lopis, E. Donath, and S. E. Moya, “Spontaneous confocal Raman microscopy: a tool to study the uptake of nanoparticles and carbon nanotubes into cells,” Nanoscale Res. Lett. 6(1), 429 (2011). 9. H. Kim, C. A. Michaels, G. W. Bryant, and S. J. Stranick, “Comparison of the sensitivity and image contrast in spontaneous Raman and coherent Stokes Raman scattering microscopy of geometry-controlled samples,” J. Biomed. Opt. 16(2), 021107 (2011). 10. G. L. Carr, “Resolution limits for infrared microspectroscopy explored with synchrotron radiation,” Rev. Sci. Instrum. 72(3), 1613–1619 (2001). 11. D. McNaughton, “Synchrotron infrared spectroscopy in biology and biochemistry,” Aust. Biochem. 36, 55–58 (2005). 12. H.-Y. N. Holman, R. Miles, Z. Hao, E. Wozei, L. M. Anderson, and H. Yang, “Real-time chemical imaging of bacterial activity in biofilms using open-channel microfluidics and synchrotron FTIR spectromicroscopy,” Anal. Chem. 81(20), 8564–8570 (2009). #173487 $15.00 USD Received 30 Jul 2012; revised 18 Sep 2012; accepted 21 Sep 2012; published 20 Dec 2012 (C) 2012 OSA 31 December 2012 / Vol. 20, No. 28 / OPTICS EXPRESS 29694 13. E. Stavitski, M. H. F. Kox, I. Swart, F. M. F. de Groot, and B. M. Weckhuysen, “In situ synchrotron-based IR microspectroscopy to study catalytic reactions in zeolite crystals,” Angew. Chem. Int. Ed. Engl. 47(19), 3543– 3547 (2008). 14. P. Dumas, G. D. Sockalingum, and J. Sulé-Suso, “Adding synchrotron radiation to infrared microspectroscopy: what’s new in biomedical applications?” Trends Biotechnol. 25(1), 40–44 (2007). 15. E. Levenson, P. Lerch, and M. C. Martin, “Spatial resolution limits for synchrotron-based infrared spectromicroscopy,” Infra. Phys. Tech. 51(5), 413–416 (2008). 16. P. Dumas and L. Miller, “The use of synchrotron infrared microspectroscopy in biological and biomedical investigations,” Vib. Spectrosc. 32(1), 3–21 (2003). 17. H.-Y. N. Holman, H. A. Bechtel, Z. Hao, and M. C. Martin, “Synchrotron IR spectromicroscopy: chemistry of living cells,” Anal. Chem. 82(21), 8757–8765 (2010). 18. G. Ellis, G. Santoro, M. A. Gómez, and C. Marco, “Synchrotron IR microspectroscopy: opportunities in polymer science,” IOP Conf. Ser.: Mater. Sci. Eng. 14, 012019 (2010). 19. M. J. Nasse, M. J. Walsh, E. C. Mattson, R. Reininger, A. Kajdacsy-Balla, V. Macias, R. Bhargava, and C. J. Hirschmugl, “High-resolution Fourier-transform infrared chemical imaging with multiple synchrotron beams,” Nat. Methods 8(5), 413–416 (2011). 20. F. Huth, M. Schnell, J. Wittborn, N. Ocelic, and R. Hillenbrand, “Infrared-spectroscopic nanoimaging with a thermal source,” Nat. Mater. 10(5), 352–356 (2011). 21. F. Lu and M. A. Belkin, “Infrared absorption nano-spectroscopy using sample photoexpansion induced by tunable quantum cascade lasers,” Opt. Express 19(21), 19942–19947 (2011). 22. B. Knoll and F. Keilmann, “Near-field probing of vibrational absorption for chemical microscopy,” Nature 399(6732), 134–137 (1999). 23. S. W. Hell, M. Dyba, and S. Jakobs, “Concepts for nanoscale resolution in fluorescence microscopy,” Curr. Opin. Neurobiol. 14(5), 599–609 (2004). 24. S. W. Hell and J. Wichmann, “Breaking the diffraction resolution limit by stimulated emission: stimulatedemission-depletion fluorescence microscopy,” Opt. Lett. 19(11), 780–782 (1994). 25. T. A. Klar, S. Jakobs, M. Dyba, A. Egner, and S. W. Hell, “Fluorescence microscopy with diffraction resolution barrier broken by stimulated emission,” Proc. Natl. Acad. Sci. U.S.A. 97(15), 8206–8210 (2000). 26. E. Rittweger, K. Y. Han, S. E. Irvine, C. Eggeling, and S. W. Hell, “STED microscopy reveals crystal colour centres with nanometric resolution,” Nat. Photonics 3(3), 144–147 (2009). 27. D. Wildanger, R. Medda, L. Kastrup, and S. W. Hell, “A compact STED microscope providing 3D nanoscale resolution,” J. Microsc. 236(1), 35–43 (2009). 28. E. Rittweger, D. Wildanger, and S. W. Hell, “Far-field fluorescence nanoscopy of diamond color centers by ground state depletion,” Europhys. Lett. 86(1), 14001 (2009). 29. S. W. Hell and M. Kroug, “Ground-state-depletion fluorescence microscopy: a concept for breaking the diffraction resolution limit,” Appl. Phys. B 60(5), 495–497 (1995). 30. S. Bretschneider, C. Eggeling, and S. W. Hell, “Breaking the diffraction barrier in fluorescence microscopy by optical shelving,” Phys. Rev. Lett. 98(21), 218103 (2007). 31. J. Kwon, Y. Lim, J. Jung, and S. K. Kim, “New sub-diffraction-limit microscopy technique: dual-point illumination AND-gate microscopy on nanodiamonds (DIAMOND),” Opt. Express 20(12), 13347–13356 (2012). 32. M. G. L. Gustafsson, “Nonlinear structured-illumination microscopy: wide-field fluorescence imaging with theoretically unlimited resolution,” Proc. Natl. Acad. Sci. U.S.A. 102(37), 13081–13086 (2005). 33. E. H. Rego, L. Shao, J. J. Macklin, L. Winoto, G. A. Johansson, N. Kamps-Hughes, M. W. Davidson, and M. G. L. Gustafsson, “Nonlinear structured-illumination microscopy with a photoswitchable protein reveals cellular structures at 50-nm resolution,” Proc. Natl. Acad. Sci. U.S.A. 109(3), E135–E143 (2012). 34. M. J. Rust, M. Bates, and X. Zhuang, “Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM),” Nat. Methods 3(10), 793–796 (2006). 35. B. Huang, W. Wang, M. Bates, and X. Zhuang, “Three-dimensional super-resolution imaging by stochastic optical reconstruction microscopy,” Science 319(5864), 810–813 (2008). 36. W. P. Beeker, P. Gross, C. J. Lee, C. Cleff, H. L. Offerhaus, C. Fallnich, J. L. Herek, and K.-J. Boller, “A route to sub-diffraction-limited CARS Microscopy,” Opt. Express 17(25), 22632–22638 (2009). 37. W. P. Beeker, C. J. Lee, K.-J. Boller, P. Groß, C. Cleff, C. Fallnich, H. L. Offerhaus, and J. L. Herek, “Spatially dependent Rabi oscillations: an approach to sub-diffraction-limited coherent anti-Stokes Raman-scattering microscopy,” Phys. Rev. A 81(1), 012507 (2010). 38. M. Fushitani, “Applications of pump-probe spectroscopy,” Annu. Rep. Prog. Chem. C 104, 272–297 (2008). 39. J. Cabanillas-Gonzalez, G. Grancini, and G. Lanzani, “Pump-probe spectroscopy in organic semiconductors: monitoring fundamental processes of relevance in optoelectronics,” Adv. Mater. (Deerfield Beach Fla.) 23(46), 5468–5485 (2011). 40. L. Carroll, P. Friedli, P. Lerch, J. Schneider, D. Treyer, S. Hunziker, S. Stutz, and H. Sigg, “Ultra-broadband infrared pump-probe spectroscopy using synchrotron radiation and a tuneable pump,” Rev. Sci. Instrum. 82(6), 063101 (2011). 41. T. Watanabe, M. Fujii, Y. Watanabe, N. Toyama, and Y. Iketaki, “Generation of a doughnut-shaped beam using a spiral phase plate,” Rev. Sci. Instrum. 75(12), 5131–5135 (2004). 42. D. Wildanger, J. Bückers, V. Westphal, S. W. Hell, and L. Kastrup, “A STED microscope aligned by design,” Opt. Express 17(18), 16100–16110 (2009). 43. J. Keller, A. Schönle, and S. W. Hell, “Efficient fluorescence inhibition patterns for RESOLFT microscopy,” Opt. Express 15(6), 3361–3371 (2007). #173487 $15.00 USD Received 30 Jul 2012; revised 18 Sep 2012; accepted 21 Sep 2012; published 20 Dec 2012 (C) 2012 OSA 31 December 2012 / Vol. 20, No. 28 / OPTICS EXPRESS 29695 44. X. Hao, C. Kuang, T. Wang, and X. Liu, “Effects of polarization on the de-excitation dark focal spot in STED microscopy,” J. Opt. 12(11), 115707 (2010). 45. D. A. Guzonas, M. L. Hair, and C. P. Tripp, “Infrared spectra of monolayers adsorbed on mica,” Appl. Spectros. 44(2), 290–293 (1990). 46. A. L. Harris, L. Rothberg, L. Dhar, N. J. Levinos, and L. H. Dubois, “Vibrational energy relaxation of a polyatomic adsorbate on a metal surface: methyl thiolate (CH3S) on Ag(111),” J


Introduction
Developments in biomedical and nano-materials' research rely increasingly on state-of-the-art instruments capable of structural imaging and spectral analysis at sub-cellular and nanometer resolution.Because vibrational modes of molecules are dependent on chemical bonds and local molecular arrangements, vibrational spectroscopies reveal molecular fingerprints, the mapping of which produces label-free imaging modalities, such as coherent anti-Stokes Raman spectroscopy (CARS) [1][2][3][4][5][6], stimulated Raman scattering [7], spontaneous Raman [8,9], and Fourier-transform infrared absorption spectroscopy (IRAS) [10][11][12][13][14][15][16][17][18][19].However, the spatial resolution of these methods is restricted by diffraction to values between λ and λ/2.This imposes a severe limitation especially for IRAS given that in the mid-infrared region the wavelength λ spans from 3 to 20 μm.The resolution for IRAS imaging in the far-field thus remains poor when compared to the resolution that can be achieved using near-field scanning probe approaches (50-100 nm) [20][21][22].Conversely, the latter techniques suffer from practical limitations arising from the scanning of a material probe a few nanometers above the surface of the sample.Achieving high resolution in far-field imaging thus remains of prime interest.
Nanoscale and sub-diffraction resolution have been extensively demonstrated in far-field fluorescence microscopy [23], with reversible optical fluorescence transition (RESOLFT) [24][25][26][27][28][29][30][31], non-linear structured illumination [32,33], and single molecule localization microscopies [34,35].However, these methods generally require the introduction of fluorescent markers in the samples and do not provide label-free local chemical information.Novel chemical imaging tools with high resolution are thus sought for and, in view of that, schemes for sub-diffraction vibrational microscopies based on far-field CARS have been proposed [36,37].For these schemes, it is envisaged to control the spatial extension of the measured sample volume by suppressing the vibrational Raman coherence through an incoherent process involving the excitation of a neighboring vibration [36] or by periodically depleting the ground-state through Rabi oscillations [37].These concepts are expected to improve the spatial resolution down to ca. 65 nm [37], which would represent a breakthrough in CARS microscopy.
In this paper, we discuss and model a framework for IRAS imaging with spatial resolution beyond the diffraction limit.In the proposed methodology, the effective point-spread function (PSF) is confined by generating a spatial contrast in the depletion of the ground-state of a vibrational mode of interest using a scheme of IR pumping and IR probing, which results in a PSF proportional to the local IR absorbance.Here below, the method is termed as vibrationaldepletion IR (VD-IR) microscopy, and it is illustrated by simulating PSF and images of patterned thin films of octadecyltrichlorosilane adsorbed on an IR transparent dielectric substrate.

Vibrational depletion IR microscopy
Figure 1(a) illustrates the concepts of IR absorption and pump-probe spectroscopy [38][39][40], where both pump and probe pulses are tuned to the wavelength of a vibrational transition.A pulse of low intensity undergoes partial absorption through a sample of two level oscillators and the absorption of each IR photon leads to the excitation of one oscillator from the ground state 0 to the excited state 1 (left).The competition between stimulated absorption and emission processes prevents the population of 1 to exceed that of 0 , so that a pulse of very high intensity excites at most half the oscillators, leading to the vibrational saturation of the sample (middle).In the pump-probe experiment, the same pulse of very high intensity is immediately followed by a delayed probe pulse whose absorption is reduced, for delays shorter than the lifetime of 1 (right).In other words, this second pulse probes the groundstate population in the sample and regions of the sample that have been irradiated by the pump are virtually transparent to the probe for short delays, where the ground-state remains depleted.
A suitable PSF for VD-IR imaging is realized by alternating the spatial definition of the pump, by using a first one characterized by a central anti-node and a second by a central node, as shown in Fig. 1(b).The absorption of the probe is minimized or nullified over its entire spatial extension in the first case, thus creating a reference signal, and absorption of the probe is confined to the nodal region of the pump in the second case.The difference in integrated pulse intensity gives the probe absorption confined to the center of the nodal pump, where the sample population 1 remains below saturation, and is defined here as the VD-IR signal.Increase in pump intensities reduces the spatial extension of the unsaturated center and thus increases further the localization of the VD-IR PSF, since sample regions that are saturated for both pump patterns do not contribute to the differential signal.
Although various geometry are possible, the spatial distribution of the nodal pump intensity is set to that of a first-order Gaussian which can be readily generated passing a zeroth-order Gaussian beam through a vortex phase plate [41][42][43][44].The latter modify the phase in the plane normal to the direction of propagation by a factor i e θ , where θ is the polar angle in the plane, effectively shaping the wavefront into a spiral, so that destructive interferences occur at the beam center.The nodal pump is below referred to as vortex.For simplicity, the anti-nodal pump is defined as a Gaussian, concentric with the vortex pump as well as with the Gaussian probe.

Formalism
To study the feasibility of the proposed concept and predict realistic PSF, VD-IR images of a CH 2 stretch vibration at ca. 3.5 μm in octadecyltrichlorosilane were simulated, taking into account of empirical values of the absorbance and relaxation time of the CH stretch modes [45][46][47][48][49][50][51][52][53][54][55].The change in intensity due to stimulated absorption/emission [56], along the coordinate z defined as the axis of propagation of the light pulses, is readily inferred from where k defines the stimulated emission/absorption coefficient, ρ is the density of oscillators, N Δ is the difference in relative population density between the levels 1 and 0 , ( , ) I t r is the local intensity at a given time t and sample position ( , , ) x y z = r , with ( , ) x y the coordinates in the sample plane set normal to the direction of propagation of the light.h and c are the Planck constant and the speed of light in vacuum.N Δ is determined by where Γ is the deexcitation rate of 1 , and N the relative population density of the oscillators in the excited state 1 at a given time and position in the sample.
The sample was modeled in space as an ensemble of voxels each containing a collection of independent oscillators, and a code was written in Matlab to solve the system of equations.The latter was developed using the backward Euler approach and independently for each set of ( , ) x y coordinates.The solution is readily found by iteration with boundary conditions 0 ( , ) 0 describing the temporal evolution of the intensity impinging the sample at 0 z for the copropagative pump and probe pulses.
The departing probe pulse energy Σ was computed by integrating the intensity in the sample plane ( , ) x y and over time, at the coordinate z exceeding the sample thickness.IR absorption images and sub-diffraction VD-IR images were computed by repeating the calculation whilst varying the relative position of the sample with respect to the pulses.The IR absorption and VD-IR signals are defined by 0 0 ( ) (%) 100, ( ) and by (%) 100, where 0 ( ) is the probe pulse energy incident on the sample, and Gauss Σ and vortex Σ are the probe pulse energy transmitted through the sample following a Gaussian or a vortex pump pulse, respectively.The calculation neglects scattering and diffraction in the sample, which is justified since its thickness is of a few microns maximum in transmission IRAS microscopy and thus typically shorter than the wavelength, and since the whole transmitted intensity is integrated over the sample plane.
In the focal plane, the spatio-temporal intensity profile of the Gaussian pulses is of circular symmetry and defined by

Empirical considerations
Aiming to a realistic prediction of the VD-IR microscopy performances, the fwhm of the Gaussian pulses were adjusted to those expected with objectives of NA 0.7 and 0.85.The first value is chosen because it is the largest NA reported for reflective objective, such as those used in synchrotron IR absorption microscopy [10,58], and the second value because aspherical lenses at 3.5 μm with a NA of 0.85 are commercially available.At the diffraction limit, experiments have verified that an objective with NA 0.7 focuses a Gaussian beam with λ at 3.5 μm into a spot of fwhm ca.2.3 μm [10].A fwhm of ca.1.9 μm is thus readily expected for a NA of 0.85.The beam waists 0 w are adjusted to reproduce these fwhms for both the Gaussian pulses and the same waists are used for the vortex pulses.
Pump-probe spectroscopy of vibrational modes is best achieved with picosecond-long pulses [55], affording suitable time resolution without compromising too much on the spectral resolution, and a duration of 1.0 ps was then chosen here for all pulses.The spatial distribution of the intensity of Gaussian and vortex pulses are shown in Fig. 2(a), along with their radial profile computed for an energy of 1.0 nJ.The sample emulates a thin film of octadecyltrichlorosilane self-assembled on a fully transparent dielectric substrate (see Fig. 2(b)), for which experimental measurements reveals ca.0.5% of absorption by the CH 2 stretch modes in homogeneous films [45].This value is reproduced with Eq. ( 1) and Eq. ( 2) for k equals 1.4 nm 2 .fJ−1 and ρ 31.5 nm −3 , accounting for a uniform thickness of 2.0 nm confined into a single pixel, an area of 27 Å 2 per molecule, and 17 CH 2 moieties per molecule.
A sketch of the optical path for the Gaussian probe and the alternating Gaussian and vortex pumps is presented in Fig. 2(c).The pump pulses are focused on the octadecyltrichlorosilane sample and followed by the delayed Gaussian probe pulse which is separated after transmission through the sample by a time-gate and eventually integrated by a detector.In line with Eq. ( 4), the VD-IR signal is built as the difference between the transmitted probe pulse energies at the detector with Gaussian and vortex pumps.Images are generated by scanning the sample in the beam paths.
The intensity at which the saturation of N occurs has been computed by observing the evolution of the population of 1 in a single pixel as a function of the intensity of a pulse of duration 1.0 ps (see Fig. 3).The data reveals a saturation threshold of ca.1.1 × 10 2 kW.μm −2 for the emulated sample which is consistent with earlier observations [55].From Fig. 2(a) and Eq. ( 6), it is seen that energies above 10 nJ are required to achieve saturation with the vortex pump.These energies and intensities are readily available in the mid-IR with tunable picosecond sources based on non-linear optical conversion and pumped by pulsed lasers.

Point-spread functions
The PSF for the IR absorption and VD-IR microscopies have been simulated by computing the (%) IR and (%) VD IR − images for a sample where the molecules are found in a single pixel of dimension 25 × 25 nm 2 , and where all others pixels are left transparent (i.e., k , Γ , and ρ are all null).Since pumps and probe are here always of circular symmetry, only plots of the radial evolution of the PSF are shown.The Gaussian and vortex pumps were set to the same energy and waist, and varied from 1.0 nJ to 1.0 μJ, and the Gaussian probe at 0.1 nJ.The PSF are shown in Fig. 4(a) for a NA of 0.85 and a pump-probe delay of 2.0 ps.
The PSF for the diffraction limited IR absorption image is shown as a reference and reproduces the expected fwhm of 1.9 μm, commensurate with the dimension of the Gaussian probe in the sample plane.From the concept illustrated in Fig. 1(b), the VD-IR PSF is expected to be narrower and this is confirmed by the simulated data presented in Fig. 4(a) and Fig. 4(b).With the pumps at 1.0 nJ, thus about an order of magnitude below the saturation threshold, the VD-IR PSF is already reduced down to 1.05 μm.This value is about half the radial distance between two intensity maxima in the vortex pattern (2.2 μm) and is below-thediffraction-limit (1.9 μm).The improvement in resolution is reminiscent of structured illumination microscopy that exploits sinusoidal illumination patterns to double the spatial resolution in wide-field microscopy [59].Increasing the pump energy to values that are sufficient for the intensity to exceed locally the saturation threshold leads to a further improvement in resolution with a fwhm of 250 nm at 0.1 μJ (i.e., for a maximum vortex intensity ca. 8 times above the saturation threshold), and below 100 nm for higher energies.Figure 4(b) highlights the influence of the dimensions of the vortex profile.Replacing the objective by one of NA 0.7 increases the diameter of the ring of maximum intensity from 2.2 μm to 2.85 μm, and the PSF fwhm is scaled up accordingly, with for example a value of 390 nm at 0.1 μJ.The value also reflects the decrease in local intensity at same pulse energy due to the increase in irradiated area.
Although the PSF for VD-IR microscopy exhibits mainly a point-like geometry, negative contributions at large radial displacements are seen when the excited population N in the outer rim of the probed area differs between vortex and Gaussian pumps (see Fig. 4(a) at 1.0 nJ).However, these negative contributions vanish for large pump energy when the entire probed area is effectively saturated for both pumps.
The magnitude of the VD-IR signal when the pixel of molecules coincides with the center of the beams is presented in Fig. 4(c pulse energies, with NA at 0.85.For short delays such as 2.0 ps and energies close or above the saturation threshold, the VD-IR signal corresponds to the same percentage of intensity measured at the detector as for the IR absorption microscopy.Figure 4(c) illustrates the exponential decrease of the VD-IR signal when the pump-probe delay t Δ is increased.The decay occurs with a rate equivalent to the decay rate Γ of the excited oscillators (i.e., 0.1 ps −1 ), which is understandable given that at longer delays the probe interrogates a population that is exponentially decreased (see inset in Fig. 3).IR and VD-IR images at 3.5 μm have been simulated with the pump set at 0.10 μJ for an array of nine domains of octadecyltrichlorosilane (200 × 200 nm 2 ).As shown in Fig. 5, the pattern is clearly resolved and the VD-IR signal is then ca.0.4 × 10 -2% of the total energy fluency at the detector.The signal-to-noise ratio required for measuring a molecular domain of 200 × 200 nm 2 remains thus of the order of 10 −5 -10 −6 (i.e., for a VD-IR magnitude a decade above noise level), a value that remains acceptable in view of the noise-equivalent power in modern IR detectors (e.g., 0.2 pW .Hz −1/2 for mercury cadmium telluride) [20].Indeed, using such a detector in conjunction with lock-in amplification, Fourier-transform thermal-source IR SNOM imaging has been demonstrated with a spatial resolution below 100 nm and an integration time of 40 ms per pixel, in a configuration where the scattered power was only 10 −8 that of the power in the reference arm [20].A comparison between Fig. 5(d) and Fig. 4(a) shows that the VD-IR signal scales linearly with the molecular domain dimensions, with maximum contrasts of ca.0.4 × 10 -2% and 0.7 × 10 -4% for domains of 200 × 200 nm 2 and 25 × 25 nm 2 , respectively.Figure 4(a) also shows that for a given molecular domain the signal contrast for VD-IR and IR microscopy is the same at short pump-probe delays and that both methods have the same detection capability.Thus, although the instrumental noise (i.e., laser source, time-gate, detector, etc) will ultimately limit the sensibility of VD-IR, the signal contrasts are found sufficient in the sub-micron and nanoscale regimes.

Damage threshold
Akin to the RESOLFT techniques affording nanoscale resolution in fluorescence microscopy [24][25][26][27][28][29][30] the resolution for VD-IR is theoretically bound only by the pump energy, and the latter can be, practically, increased up to the damage threshold of the sample.A brief survey of past IR pump-probe experiments of CH stretch modes in organic materials and of the damage threshold in inorganic materials is thus useful to assess the potential performance of VD-IR.Because pulsed lasers are necessary to achieve peak intensities above the saturation threshold, the averaged power that will be received by the sample is relatively moderate.Assuming a standard repetition rate of 1.0 kHz, the power here is at most 1.0 mW, since 1.0 μJ is the highest pulse energy considered in this study.This is a relatively modest power and similar to the values used in near-field probe imaging of polymer beads with a laser power of ca. 1 mW focused within a 10 μm beam spot [22].The power values are also similar to those used in synchrotron IRAS microscopy of samples including biological cells where the broadband beam has been focused without damage on the same areas on samples [60,61].
Bearing in mind the typical temporal structure of a synchrotron source [11], the much shorter pulses required for VD-IR imply higher peak intensities and their values with respect to the damage thresholds in organic and inorganic materials must be considered.The average energy density of the vortex pump at the sample is 5 nJ.μm −2 for a pulse of 0.10 μJ that is focused on an area of radius ca.2.5 μm (see Fig. 2(a), NA = 0.85) and for which a resolution of about λ/14 is calculated (see Fig. 4).The energy density expected for 0.10 μJ is thus about an order of magnitude below the damage threshold of dielectric substrates such as CaF 2 and silicate exposed to ps IR pulses [62,63].Moreover, sum-frequency generation (SFG) images of the CH stretch modes in organic monolayers [64] and IR pump-SFG probe spectra of CH stretch mode on diamond [48] have been successfully acquired with an energy density of ca. 10 nJ.μm −2 , suggesting that the damage threshold is close to or above that value.It seems thus reasonable to assume that VD-IR affords spatial resolutions better than 250 nm (λ/14) for these materials and that a resolution below 100 nm is also possible (see Fig. 4(b)).
For molecules in a solution or at the interface with another a brief survey of IR pump-probe spectroscopy investigations reveals that saturation of CH in organic molecules is typically made with pulse energy densities of the order of several nJ.μm −2 [52,53,65].For VD-IR microscopy, the pump energy will thus be preferably limited to values ca.0.10 μJ resulting in a spatial resolution of the order of 250 nm or higher.For IR nanosecond pulses, the ablation threshold of soft tissue has been estimated at ca. 20 nJ.μm −2 [66] and ablation experiment supports the negligible heating of tissue at low repetition rates [67], with complete cooling between pulses.

Conclusions
A conceptual framework has been proposed for imaging the IR absorbance of samples at subdiffraction spatial resolutions and simulation of the associated PSF has been discussed.The pulse energies required to pattern the vibrational depletion to achieve imaging beyond the diffraction limit have been determined on the basis of empirical parameters.It has been found that, for energy density of several nJ.μm −2 i.e., roughly 10 times above the saturation threshold, the fwhm of the VD-IR PSF can be as low as λ/14, which highlights the possibility of IRAS imaging beyond the diffraction limit.For monochromatic IR pulses, these energy densities are estimated to be below the damage threshold for thin organic films on dielectric or metallic substrate, and close to for soft tissue and organic molecules at soft interfaces.Chemical imaging based on the measure of the IR absorption by vibrational modes in the farfield is thus possible at sub-micron resolution and down to the nanoscale.

Fig. 1 .
Fig. 1.(a) Illustration of the quantized IR absorption in a two-level system (left), of the saturation where the excited population is at maximum (middle), and of the pump-probe sequence where the population is saturated by a first pulse and probed by a second delayed pulse (right).(b) Description of the differential pumping scheme where the VD-IR PSF is created by alternating nodal and anti-nodal pump profiles.
where 0 t is a reference time before irradiation of the sample, and 0 ( , , , ) I x y z t #173487 -$15.00USD Received 30 Jul 2012; revised 18 Sep 2012; accepted 21 Sep 2012; published 20 Dec 2012 (C) 2012 OSA r is the radial coordinate in the plane normal to the direction of propagation, and where the full width at half maximum (fwhm) of the Gaussian is defined by 0 2 ln(2) w , with 0 w the Gaussian waist, and the temporal pulse duration by 0 2 ln(2)τ . 0Gauss h is a constant adjusted to reproduce pulse energies ranging from 1.0 nJ to 1.0 μJ for the Gaussian pump pulses and 0.1 nJ for the probe.tΔ is zero for the pump pulses and set to a finite value for the probe, marking the pump probe delay.The vortex pulse intensity in the focal plane is written[57] adjusted to the desired vortex pump energy.

Fig. 2 .
Fig. 2. (a) Spatial distribution (scale bar: 2.0 μm) of the intensity of Gaussian and vortex pulses computed for an energy of 1.0 nJ, a duration of 1.0 ps, and a waist of 1.6 μm (Gaussian fwhm: 1.9 μm), along with their corresponding radial profile.(b) Schematic of an octadecyltrichlorosilane layer on a dielectric IR transparent substrate.(c) Sketch of the optical path for the Gaussian probe and the alternating Gaussian and vortex pumps.BS, SPL, L, TG, D stand for beam splitter, sample, lens, time-gate, and detector.

Fig. 3 .
Fig. 3. Local evolution of the population N as a function of the intensity.A saturation threshold of 1.1 × 10 2 kW.μm −2 is determined. (inset) Evolution of N with time at the intensity peak of a Gaussian pulse centered at 0 ps with energy of 1 nJ (continuous line), 10 nJ (shortest dash), 0.10 μJ (long dash), and 1.0 μJ (longest dash).

Fig. 5 .
Fig. 5. (a) Representation of an array of nine domains of octadecyltrichlorosilane (200 × 200 nm 2 ).(b) CH 2 stretch IR image (3.5 μm) computed with a pixel of 50 × 50 nm 2 and a NA of 0.85.(c) CH 2 stretch VD-IR image computed for a pump energy of 0.1 μJ.Scale bar is 0.5 μm in all images.(d) IR and VD-IR magnitudes across the octadecyltrichlorosilane domains.