**Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging**

Shuanglong Liu, Wei Liu and Hanben Niu

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/59987

### **1. Introduction**

Supercontinuum (SC) generation describes the phenomenon that the output spectrum becomes much wider than the input spectrum when an ultrashort optical pulse passes through a nonlinear medium, such as bulk fused silica, photonic crystal fiber and so on. Since this phenomenon was first discovered by Alfano and Shapiroin in 1970s by passing picosecond pulses through a bulk BK7 glass [1], it has been explored in a wide variety of nonlinear media, including solids, gas, inorganic liquids and various waveguides.

In 1974, an ingenious idea of introducing microstructures into traditional step index optical fibers to modify the guidance properties by Kaiser and Altle [2] established the theoretical foundation for photonic crystal fiber (PCF). And the exploration for fabricating such fibers were never stopped until the advent of an all-silica PCF in Southampton University by Russell and co-workers in 1996 [3]. His birth drew much more interest from various scientific re‐ searchers and initiated a revolution in SC generation.

The PCF used for SC generation often has a solid core surrounded by regularly arranged microscopic air holes running along the fiber length [4], resulting in a surrounding region with lower refractive index compared to that of the central district. Therefore the light will be trapped and guided in the fiber core through total internal reflection. This modified surround‐ ing shows good flexibility to engineer the waveguide properties. It is reported that the zerodispersion wavelength (ZDW) in fused silica can shift to a shorter wavelength than the intrinsic ZDW about 1300 nm by suitable design of the air hole diameter and air filling fraction [5]. Then Ranka *et al.* first reported that the PCF can exhibit anomalous dispersion at visible region [6]. And they obtained an ultrabroadband continuum generation from 400 nm to 1600 nm by injecting a femtosecond pulse into a 75 cm section PCF. In 2000, Wadsworth *et al.* observed soliton effects by pumping a short PCF near the ZDW in the anomalous regime [7]. Accom‐

© 2015 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and eproduction in any medium, provided the original work is properly cited.

panied with other nonlinear effects, the fundamental and second soliton were identified. The phenomenon of soliton self-frequency shift [8] in a tapered air–silica microstructure fiber was reported by Liu *et al.* from 1.3 to 1.65 µm. Numerical modeling of pulse evolution in PCFs was first proposed by Husakou and Herrmann [9] based on reduced Maxwell equations. And the important roles of soliton fission, four-wave mixing in spectral broadening were first discussed in detail. Due to the lack of Raman scattering, much fine structure was lost in the solution. An extended nonlinear Schro¨dinger equation (NLSE) was first applied to the modeling of femtosecond pulse propagation in PCF [10] with higher-order dispersion, self-steepening and Raman effects included. The numerical simulations based on extended NLSE is accurately consistent with experimental results. Based on the previous works, Genty *et al.* quantitatively studied the mechanisms of SC generation when the PCF was pumped in its normal and anomalous regimes [11]. Soon afterwards, continuum generation from 1065 to 1375 nm was demonstrated by pumping a honey fiber with continuous wave [12]. These numerical and experimental studies led to a detailed description of the SC generation in PCFs from the ultraviolet to the infrared in both femtosecond, picosecond regimes and continuous wave.

SC generation in PCFs involves several nonlinear effects such as self-phase modulation, (SPM) cross-phase modulation (XPM), self-steepening and optical shock formation, Raman scatter‐ ing, solition fission, four-wave mixing and so on. Greatest spectral broadening was reported when a PCF was pumped in the anomalous dispersion regime near the ZDW [13]. The broadening mechanism in this case is dominated by soliton dynamics, especially breakup of higher-order solitons into fundamental solitons and non-solitonic radiation, though the SPM effect in initial propagation. The ejected solitons will shifts toward longer wavelength under intrapulse Raman scattering, known as soliton self-frequency shift (SSFS). Spectral broadening in normal dispersion region is mainly due to SPM effect and optical wave breaking. With soliton dynamics suppressed in this region, the generated SC possesses a good temporal structure in spite of a smaller spectral broadening. Under specific pump pulses, continua with broad spectral bandwidth, good spatial coherence, uniform intensity profile and simple temporal structure can be generated in certain PCFs. These continua show significant appli‐ cations in many fields such as optical frequency metrology, optical communication and cellular biology.

Recent years have seen the fast development of cellular biology based on optical microscope and coherent lasers. And many microscopies were developed to observe the fine structures inside biologic cells, such as scanning electron microscope [14] (SEM), scanning tunneling microscope [15] (STM), and atomic-force microscope [16] (AFM). Though they all have very high spatial resolutions, vacuum circumstance needed in SEM makes it not suit for live-cell imaging while STM and AFM are restricted to surface mapping due to the cantilever tip [17]. With the development of several novel microscopic methods, such as photo-activated locali‐ zation microscopy [18] (PALM), stochastic optical reconstruction microscopy [19] (STORM), stimulated emission depletion [20] (STED) microscopy, fluorescence microscopies obtain a spatial resolution about 20 nm and are widely used in cellular biology. But photo-toxicity, photo-bleaching and influence of invasive marker on cells cannot be ignored in vivo imaging, let alone some biomolecule that are difficult to be labeled [21]. Therefore it is urgent to develop a microscopy with nano-scaled resolution to study the process of metabolism in live cells and its response to invasive substances.

panied with other nonlinear effects, the fundamental and second soliton were identified. The phenomenon of soliton self-frequency shift [8] in a tapered air–silica microstructure fiber was reported by Liu *et al.* from 1.3 to 1.65 µm. Numerical modeling of pulse evolution in PCFs was first proposed by Husakou and Herrmann [9] based on reduced Maxwell equations. And the important roles of soliton fission, four-wave mixing in spectral broadening were first discussed in detail. Due to the lack of Raman scattering, much fine structure was lost in the solution. An extended nonlinear Schro¨dinger equation (NLSE) was first applied to the modeling of femtosecond pulse propagation in PCF [10] with higher-order dispersion, self-steepening and Raman effects included. The numerical simulations based on extended NLSE is accurately consistent with experimental results. Based on the previous works, Genty *et al.* quantitatively studied the mechanisms of SC generation when the PCF was pumped in its normal and anomalous regimes [11]. Soon afterwards, continuum generation from 1065 to 1375 nm was demonstrated by pumping a honey fiber with continuous wave [12]. These numerical and experimental studies led to a detailed description of the SC generation in PCFs from the ultraviolet to the infrared in both femtosecond, picosecond regimes and continuous wave.

SC generation in PCFs involves several nonlinear effects such as self-phase modulation, (SPM) cross-phase modulation (XPM), self-steepening and optical shock formation, Raman scatter‐ ing, solition fission, four-wave mixing and so on. Greatest spectral broadening was reported when a PCF was pumped in the anomalous dispersion regime near the ZDW [13]. The broadening mechanism in this case is dominated by soliton dynamics, especially breakup of higher-order solitons into fundamental solitons and non-solitonic radiation, though the SPM effect in initial propagation. The ejected solitons will shifts toward longer wavelength under intrapulse Raman scattering, known as soliton self-frequency shift (SSFS). Spectral broadening in normal dispersion region is mainly due to SPM effect and optical wave breaking. With soliton dynamics suppressed in this region, the generated SC possesses a good temporal structure in spite of a smaller spectral broadening. Under specific pump pulses, continua with broad spectral bandwidth, good spatial coherence, uniform intensity profile and simple temporal structure can be generated in certain PCFs. These continua show significant appli‐ cations in many fields such as optical frequency metrology, optical communication and cellular

Recent years have seen the fast development of cellular biology based on optical microscope and coherent lasers. And many microscopies were developed to observe the fine structures inside biologic cells, such as scanning electron microscope [14] (SEM), scanning tunneling microscope [15] (STM), and atomic-force microscope [16] (AFM). Though they all have very high spatial resolutions, vacuum circumstance needed in SEM makes it not suit for live-cell imaging while STM and AFM are restricted to surface mapping due to the cantilever tip [17]. With the development of several novel microscopic methods, such as photo-activated locali‐ zation microscopy [18] (PALM), stochastic optical reconstruction microscopy [19] (STORM), stimulated emission depletion [20] (STED) microscopy, fluorescence microscopies obtain a spatial resolution about 20 nm and are widely used in cellular biology. But photo-toxicity, photo-bleaching and influence of invasive marker on cells cannot be ignored in vivo imaging, let alone some biomolecule that are difficult to be labeled [21]. Therefore it is urgent to develop

biology.

28 Photonic Crystals

As a new kind of microscopies based on molecular vibrations, coherent anti-Stokes Raman scattering (CARS) microscopy has exhibited their prospective applications in live-cell imaging for its distinct characteristics, including high sensitivity and spatial resolution, label-free chemical specificity, three-dimension sectioning capability [22-24]. In traditional CARS microscopy, the pump and Stokes pulses can resonant only single or few molecular bonds due to the limitation of the spectral bandwidth. Obviously it is not adequate to acquire accurate recognition and mapping of an unknown biomolecule or that with complex components. Although this can be achieved by sequentially tuning the frequency difference of them or synchronizing two ultrashort lasers, it is time-consuming and costs a lot. When a SC pulse generated in a high nonlinear PCF serves as both the pump the Stokes, most distinct vibration modes, perhaps full modes even, can be probed simultaneously as long as different frequency components in the generated SC propagate at about a same speed. And a method named additional probe-beam-induced phonon depletion [25] (APIPD) is proposed to improve the spatial resolution on CARS microscopy. In APIPD method, an additional doughnut beam with different frequency from that of the probe beam is introduced to deplete the phonons on the periphery of point spread function (PSF). The subsequent Gaussian probe beam reacts with the rest phonons, generating an anti-Stokes signal near centre of the PSF. By filtering out the signal induced by doughnut beam, the effective PSF is hence decreased, meaning a higher spatial resolution. And the resolution can be reduced below 100 nm by properly modifying the probe and additional probe beam. These characters make CARS microscopy especially suit for studying the process of metabolism in live cells.

In the latest 10 years, infrared (IR) microscope has been widely used in many industries, such as material analysis and cellular biology, for its ability of nondestructive imaging and molec‐ ular location based on vibrational spectroscopy [26,27]. Traditional IR microscopes are based on either thermal IR sources [28] with low brightness and coherence, such as a globar or Hglamp, or synchrotron radiation [29,30]. The synchrotron radiation is complex, expensive, power consuming and always accompanied with intensity fluctuation in spite of high brightness and coherence. The application of non-silica fibers [31-33] brings about a new IR source with both high brightness and coherence, not to mention the broad bandwidth.

In this chapter, the generalized nonlinear Schro¨dinger equation is briefly deduced, based on the Maxwell's equations. It takes not only the higher-order dispersion into consideration, but also the higher-order nonlinear effects including self-steepening, optical shock formation and intrapulse Raman scattering. Then a predictor-corrector split-step Fourier method is used to simulate the pulse evolutions in both temporal and spectral domain for its high accuracy and fast calculating speed. Based on the fiber structure, the dispersion parameters and nonlinear coefficient can be calculated with finite element method. Then the process of continuum generation in a PCF is described in detail, including pulse evolutions in both temporal and spectral domain. The generated SC has important applications in nonlinear optics, especially CARS microscopy and IR microscopy. Based on the vibrational spectroscopy, CARS micro‐ scopy and IR microscopy has been widely used in cellular biology for their label-free and nondestructive imaging. In terms of the limited spatial resolution in CARS microscopy, an APIPD is proposed to break the diffraction barrier, leading to nano-scaled imaging. Further‐ more, the combination of this CARS nanoscopy and the broadband continuum is of crucial importance to study the fine structures and metabolic dynamics in live cells. Besides, based on the atomic force microscope, FTIR microscopy has realized a spatial resolution below 100nm. Serving a mid-infrared supercontinuum generated with the chalcogenide PCF as pumping source, this nano-FTIR will become a powerful tool for chemical identification of unknown nanostructures.

#### **2. Numerical modeling of supercontinuum generation**

#### **2.1. Nonlinear propagation equation**

The propagation equation describing evolution of laser pulses in optical fibers could be derived from Maxwell's equations [34]

$$
\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial \mathbf{t}} \quad \nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial \mathbf{t}} \quad \nabla \cdot \mathbf{D} = \mathbf{Q}\_{\mathbf{f}\mathbf{f}} \quad \nabla \cdot \mathbf{B} = 0 \tag{1}
$$

where **E** and **B** are electric field vector and magnetic flux density, **H** and **D** are magnetic field vector and electric flux density respectively. **J** and ρ<sup>f</sup> are current density vector and free charge density. They could be related to themselves through the constitutive relations given by

$$\mathbf{D} = \varepsilon\_0 \mathbf{E} + \mathbf{P}, \quad \mathbf{B} = \mu\_0 \mathbf{H} + \mathbf{M} \tag{2}$$

where ε0 and µ0 are vacuum permittivity and permeability, **P** and **M** are the induced electric and magnetic polarizations respectively. In optical fibers, **J**, ρ<sup>f</sup> and **M** all equal to zero, for they are nonmagnetic and absent of free charges.

Based on the Maxwell's equations and constitutive relations, wave equation describing light propagation in optical fibers will be

$$\nabla^2 \mathbf{E} \cdot \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \mu\_0 \frac{\partial^2 \mathbf{P}}{\partial t^2} \tag{3}$$

The induced electric polarization **P** , which in essence expresses the interaction between light and medium, will depend on the electric field strength as

$$\mathbf{P} = \varepsilon\_0 \mathbf{I} \chi^{(1)} \bullet \mathbf{E} + \chi^{(2)} : \mathbf{E} \mathbf{E} + \chi^{(3)} : \mathbf{E} \mathbf{E} \mathbf{E} + \dots \mathbf{J} \tag{4}$$

where χ(k) (k=1,2, 3, ...) is the k-order susceptibility. The linear susceptibility χ(1) often contrib‐ utes to the refractive index and fiber loss, while the second order susceptibility *χ* (2) are related to the second order nonlinearity effects such as second

harmonic generation (SHG), sum frequency generation (SFG). The third order susceptibility *χ* (3) is ofter accompanied by third harmonic generation (THG), four wave mixing (FWM) and so on. The second order susceptibility will be zero, because the fiber is Circular symmetric. So the induced electric polarization is

$$\mathbf{P} = \varepsilon\_0 \mathbf{\bar{L}} \chi^{(1)} \bullet \mathbf{E} + \chi^{(3)} \text{ : } \mathbf{EEE} + \dots \mathbf{J} \tag{5}$$

Also the induced electric polarization can be divided into two parts, the linear component **P***<sup>L</sup>* and nonlinear component **P***NL* given by PL <sup>=</sup>ε0χ(1) <sup>∙</sup>**E**, PNL <sup>=</sup>χ(3) ⋮**EEE**. The relations between linear component, nonlinear component and electric field could expressed as

PL =ε0*∫* -*∞ <sup>t</sup>* χ(1) (*t* - *t '*)∙**E**(*r*, *t '*)*dt '* (6)

$$\mathbf{P}\_{NL} = \varepsilon\_0 \int\_{\omega} dt\_1 \int\_{\omega}^{t} dt\_2 \int\_{\omega} dt\_3 \mathbf{j}\_{\omega}^{\dagger} \chi^{(3)}(t \text{ - } t\_1 \text{ - } t \text{ - } t\_2 \text{ - } t \text{ - } t\_3) \times \mathbf{E}(r, t\_1) \mathbf{E}(r, t\_2) \mathbf{E}(r, t\_3) \tag{7}$$

So Eq. (3) can also be

nondestructive imaging. In terms of the limited spatial resolution in CARS microscopy, an APIPD is proposed to break the diffraction barrier, leading to nano-scaled imaging. Further‐ more, the combination of this CARS nanoscopy and the broadband continuum is of crucial importance to study the fine structures and metabolic dynamics in live cells. Besides, based on the atomic force microscope, FTIR microscopy has realized a spatial resolution below 100nm. Serving a mid-infrared supercontinuum generated with the chalcogenide PCF as pumping source, this nano-FTIR will become a powerful tool for chemical identification of

The propagation equation describing evolution of laser pulses in optical fibers could be derived

where **E** and **B** are electric field vector and magnetic flux density, **H** and **D** are magnetic field

where ε0 and µ0 are vacuum permittivity and permeability, **P** and **M** are the induced electric

Based on the Maxwell's equations and constitutive relations, wave equation describing light

The induced electric polarization **P** , which in essence expresses the interaction between light

∶**EE** + χ(3)

(k=1,2, 3, ...) is the k-order susceptibility. The linear susceptibility χ(1)

utes to the refractive index and fiber loss, while the second order susceptibility *χ* (2)

<sup>∂</sup> <sup>t</sup> , ∇ ∙**D**=ρ<sup>f</sup>

, ∇ ∙**B**=0 (1)

and **M** all equal to zero, for they

are current density vector and free charge

<sup>2</sup> (3)

⋮**EEE** + ... (4)

often contrib‐

are related

**D**=ε0**E** + **P**, **B**=µ0**H** + **M** (2)

∂ **D**

density. They could be related to themselves through the constitutive relations given by

**2. Numerical modeling of supercontinuum generation**

<sup>∂</sup> <sup>t</sup> , ∇ ×**H**=**J** +

unknown nanostructures.

30 Photonic Crystals

**2.1. Nonlinear propagation equation**

<sup>∇</sup> <sup>×</sup>**E**<sup>=</sup> - <sup>∂</sup>**<sup>B</sup>**

are nonmagnetic and absent of free charges.

propagation in optical fibers will be

where χ(k)

vector and electric flux density respectively. **J** and ρ<sup>f</sup>

and magnetic polarizations respectively. In optical fibers, **J**, ρ<sup>f</sup>

and medium, will depend on the electric field strength as

**<sup>P</sup>**=ε<sup>0</sup> <sup>χ</sup>(1)

to the second order nonlinearity effects such as second

<sup>∇</sup>2**<sup>E</sup>** - <sup>1</sup> c2 ∂<sup>2</sup> **E** ∂ t <sup>2</sup> =µ0 ∂<sup>2</sup> **P** ∂ t

∙**E** + χ(2)

from Maxwell's equations [34]

$$\nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \mu\_0 \frac{\partial^2 \mathbf{P}\_{\mathbf{L}}}{\partial t^2} + \mu\_0 \frac{\partial^2 \mathbf{P}\_{NL}}{\partial t^2} \tag{8}$$

Before solving Eq.(8), some assumptions should be made to simplify it. First **P***<sup>N</sup>* <sup>L</sup> will be treated as a smallperturbationto**P**<sup>L</sup> , because thenonlinear changes inthe refractive index is very small, less than 10-6 in fact. Then the input pulse will keep its polarization when it travels in an optical fiber, to ensure the scalar approximation reasonable. At last, the inputfield is treated as a quasimonochromatic light, which is justified when the pulse duration is more than 100 fs.

Supposing the third order susceptibility *χ* (3) has the following form

$$
\chi^{\{3\}}(t \text{ - } t\_1 \text{ - } t \text{ - } t\_2 \text{ - } t\_3) = \chi^{\{3\}} R \{ t \text{ - } t\_1 \} \delta(t \text{ - } t\_2) \delta(t \text{ - } t\_3) \tag{9}
$$

with *R*(*t*) being the nonlinear resoponse function that has been normalized. So the nonlinear induced electric polarization is

$$\mathbf{P}\_{\rm NL} = \varepsilon\_0 \chi^{\rm \tiny{\rm 3}} \mathbf{E}(r, t) \mathbf{j}\_{\rm \tiny{\rm on}}^t \mathbf{R}(t - t\_1) |\mathbf{E}(r, t\_1)|^2 \mathbf{d}t\_1 \tag{10}$$

Considering above equations, the numerical equation describing evolution of an optical pulse when it travels in a single-mode fiber is

$$\frac{\partial \mathbf{A}}{\partial x} + \frac{\alpha \mathbf{A}}{2} \mathbf{\dot{\mathbf{I}}} \text{ - } \sum\_{\mathbf{k} \mathbf{2}}^{\text{i} \times \text{t}} \mathbf{\dot{\beta}}\_{\mathbf{k}} \frac{\partial \mathbf{\dot{A}}}{\partial \mathbf{T}^{\text{k}}} = \mathbf{\dot{r}} \gamma \left\{ \mathbf{1} + \mathbf{i} \mathbf{\tau}\_{\text{shock } \overline{\text{OT}}} \frac{\partial}{\partial \mathbf{T}} \right\} \mathbf{\dot{A}} \begin{bmatrix} \mathbf{\dot{r}} \\ \mathbf{\dot{A}} \mathbf{\dot{T}} \end{bmatrix} \times \left\{ \mathbf{A} \mathbf{\dot{T}} \mathbf{\dot{T}} \right\} \tag{11}$$

This equation is known as the generalized nonlinear Schr*o* ¨dinger equation (NLSE). In Eq. (11), the left side models linear propagation effects, with A being the pulse envelope variation in a retarded time frame, α the fiber loss and βk the k-th order Taylor series expansion of mode propagation constant around the center frequency ω0 as

$$
\beta(\omega) = \beta\_0 + \beta\_1 \{\omega \cdot \omega\_0\} + \frac{1}{2} \beta\_2 \{\omega \cdot \omega\_0\}^2 + \dots \tag{12}
$$

The right side of Eq. (11) describes the nonlinear effects, such as self-steepening, Raman scattering, etc. with γ being the nonlinear coefficient expressed as γ(*ω*0) =*n*2(*ω*0)ω<sup>0</sup> / *cAeff* , where *n*2 is nonlinear refractive index, c is the speed of light in vacuum, *Aeff* is effective mode area defined as *Aeff* =(*∬*|*<sup>F</sup>* (*x*, *<sup>y</sup>*)|<sup>2</sup> *dxdy*)<sup>2</sup> / *∬*|*F* (*x*, *y*)|<sup>4</sup> *dxdy*. Here, *F* (*x*, *y*) is the modal distribution for the fundamental fiber mode. The term of τshock derives from the frequency dependence of the effective mode area *Aeff* and effective index *neff* , associated with effects such as self – steepening and optical shock formation. The relation of the time derivative term τshock and effective mode area *Aeff* , effective index *neff* can be expressed as [35,36]

$$\tau\_{\text{shock}} = \frac{1}{\omega\_0} \text{ -} \left[ \frac{1}{n\_{\text{eff}} \left(\omega\right)} \frac{\partial}{\partial \omega} n\_{\text{eff}} \left(\omega\nu\right) \right]\_{\omega = \omega\_0} \text{ -} \left[ \frac{1}{A\_{\text{eff}} \left(\omega\right)} \frac{\partial}{\partial \omega} A\_{\text{eff}} \left(\omega\nu\right) \right]\_{\omega = \omega\_0} \tag{13}$$

In Eq. (13), *τshock* could approximately equal to 1 / *ω*0 when the spectral broadening is limited to 20 THz. While the spectral broadening increases to 100 THz or more, the second and third term should be taken into account [37].

R(t) is the nonlinear response function with both the electronic and nuclear contributions in‐ cluded. Since the electronic contribution is nearly instantaneous comparing to the nuclear contribution, R(t) could be written as [38-40]

$$\mathbf{R}(t) = \begin{pmatrix} 1 \ -f \ \_R \end{pmatrix} \delta \begin{pmatrix} t \ -t \ \_e \end{pmatrix} + f \ \_Rh \, \_R \begin{pmatrix} t \ 
 \end{pmatrix} \tag{14}$$

with *f <sup>R</sup>* the fractional contribution of the delayed Raman response to nonlinear polarization, *te* the short delay in electronic response, *h <sup>R</sup>*(*t*) the Raman response function. Often *te* is very shot (<1 fs) that it could be neglected. So the nonlinear response function will be R(*t*)=(1 - *f <sup>R</sup>*)*δ*(*t*) + *f <sup>R</sup>h <sup>R</sup>*(*t*). The Raman response function has different forms for fibers made of different materials. For silica fibers, a useful form of Raman response function will be [35]

$$h\_{\,R}(t) = \left(\tau\_1^2 + \tau\_2^2\right) \exp\left(-\frac{t}{\tau\_2}\right) \sin\left(\frac{t}{\tau\_1}\right) \Theta(t) \left/\left(\tau\_1 \tau\_1^2\right)\right.\tag{15}$$

where *τ*1 and *τ*2 are two adjustable parameters, *Ѳ*(*t*) is Heaviside function. In a research in 1989, values of *τ*<sup>1</sup> =12.2 fs and *τ*<sup>1</sup> =32 fs were used to model the profile of Raman response, and the result shows good fit to the actual Raman-gain spectrum of silica fiber [35]. Based the peak Raman gain, the Raman response factor could be calculated, saying about 0.18 [38].

In Eq. (15), single Lorentzian profile was used to approximate the actual Raman gain spectrum which is widely used in investigating ultrafast nonlinear process in optical fibers for its simplicity. However, this model underestimates Raman gain below 10THz, while overesti‐ mating it beyond 15THz, making the Raman induced frequency shift not so accurate. To address this issue, another model of approximating nonlinear response function was proposed as [41]

$$\mathbf{R}(t) = \begin{bmatrix} 1 \ -f\_R \end{bmatrix} \delta(t) + f\_R \left[ (f\_a + f\_c)h\_a(t) + f\_b h\_b(t) \right] \tag{16}$$

with *f <sup>a</sup>* =0.75, *f <sup>b</sup>* =0.21, *f <sup>c</sup>* =0.04, where *ha* (*t*) and *hb* (*t*) have the following forms

$$h\_a(t) = \left(\tau\_1^2 + \tau\_2^2\right) \exp\left(-\frac{t}{\tau\_2}\right) \sin\left(\frac{t}{\tau\_1}\right) \left(\left(\tau\_1 \tau\_1^2\right)\right) \tag{17}$$

$$\mathcal{H}\_b(t) = \frac{2\mathbf{r}\_3 \cdot t}{\mathbf{r}\_3^2} \exp\left(-\frac{t}{\mathbf{r}\_3}\right) \tag{18}$$

when the parameters are *τ*<sup>1</sup> =12.2 fs, *τ*<sup>2</sup> =32 fs, *τ*<sup>3</sup> =96 fs and *f <sup>R</sup>* =0.245, the simulated Raman gain curve matchs the actual Raman gain profile very well in the whole range of 0-15 THz. To make the investigation more approximate to actual Raman response, a more accurate model could be used although it has more complicated forms [42].

#### **2.2. Frequency-resolved optical gating**

This equation is known as the generalized nonlinear Schr*o*

β(*ω*)=*β*<sup>0</sup> + *β*1(*ω* - *ω*0) +

*dxdy*)<sup>2</sup> / *∬*|*F* (*x*, *y*)|<sup>4</sup>

<sup>∂</sup> *<sup>ω</sup> neff* (*ω*) *<sup>ω</sup>*=*ω*<sup>0</sup>

effective mode area *Aeff* , effective index *neff* can be expressed as [35,36]

∂

propagation constant around the center frequency ω0 as

defined as *Aeff* =(*∬*|*<sup>F</sup>* (*x*, *<sup>y</sup>*)|<sup>2</sup>

32 Photonic Crystals

*<sup>τ</sup>shock* <sup>=</sup> <sup>1</sup>

term should be taken into account [37].

contribution, R(t) could be written as [38-40]

*h <sup>R</sup>*(*t*)=(*τ*<sup>1</sup>

<sup>2</sup> <sup>+</sup> *<sup>τ</sup>*<sup>2</sup> 2 )exp(- *<sup>t</sup> τ*2 )*sin*( *<sup>t</sup> τ*1

*<sup>ω</sup>*<sup>0</sup> - <sup>1</sup> *neff* (*ω*)

(11), the left side models linear propagation effects, with A being the pulse envelope variation in a retarded time frame, α the fiber loss and βk the k-th order Taylor series expansion of mode

1

The right side of Eq. (11) describes the nonlinear effects, such as self-steepening, Raman scattering, etc. with γ being the nonlinear coefficient expressed as γ(*ω*0) =*n*2(*ω*0)ω<sup>0</sup> / *cAeff* , where *n*2 is nonlinear refractive index, c is the speed of light in vacuum, *Aeff* is effective mode area

for the fundamental fiber mode. The term of τshock derives from the frequency dependence of the effective mode area *Aeff* and effective index *neff* , associated with effects such as self – steepening and optical shock formation. The relation of the time derivative term τshock and

In Eq. (13), *τshock* could approximately equal to 1 / *ω*0 when the spectral broadening is limited to 20 THz. While the spectral broadening increases to 100 THz or more, the second and third

R(t) is the nonlinear response function with both the electronic and nuclear contributions in‐ cluded. Since the electronic contribution is nearly instantaneous comparing to the nuclear

with *f <sup>R</sup>* the fractional contribution of the delayed Raman response to nonlinear polarization, *te* the short delay in electronic response, *h <sup>R</sup>*(*t*) the Raman response function. Often *te* is very shot (<1 fs) that it could be neglected. So the nonlinear response function will be R(*t*)=(1 - *f <sup>R</sup>*)*δ*(*t*) + *f <sup>R</sup>h <sup>R</sup>*(*t*). The Raman response function has different forms for fibers made of different materials. For silica fibers, a useful form of Raman response function will be [35]

where *τ*1 and *τ*2 are two adjustable parameters, *Ѳ*(*t*) is Heaviside function. In a research in 1989, values of *τ*<sup>1</sup> =12.2 fs and *τ*<sup>1</sup> =32 fs were used to model the profile of Raman response, and the


<sup>∂</sup> *<sup>ω</sup> Aeff* (*ω*) *<sup>ω</sup>*=*ω*<sup>0</sup>

R(*t*)=(1 - *f <sup>R</sup>*)*δ*(*t* - *te*) + *f <sup>R</sup>h <sup>R</sup>*(*t*) (14)

)Ѳ(*t*)/ (*τ*1*τ*<sup>1</sup>

2

) (15)

¨dinger equation (NLSE). In Eq.

<sup>2</sup> *β*2(*ω* - *ω*0)<sup>2</sup> + ... (12)

*dxdy*. Here, *F* (*x*, *y*) is the modal distribution

(13)

In this chapter, we will use the predictor-corrector SSFM proposed by Lee et.al [43] to simulate evolutions of pulses in time and frequency domain. This method shows both fast calculation speed and accurate numerical results. Also the cross-correlation frequency-resolved optical gating method [44-46] (XFROG) is employed to characterize the intensity and phase profiles of the generated continuum. The method describes the pulse structure simultaneously in the time and frequency domain which is intuitive for understanding the dynamics in continuum generation.

In measurement of a continuum based on XFROG, the continuum pulse and a reference pulse should be focused on a nonlinear crystal to generate a sum-frequency signal (SFG signal) after being synchronized in the time domain by a time delay system as shown in figure 1(b). The generated SFG signals under different delay times will be recorded by a spectrometer, with which the XFROG trace could be retrieved. The XFROG trace is

$$I\_{XFORG}(\omega,\ \tau) = \left| \bigcup\_{\omega}^{\star\omega} E\_{\text{sig}}(t,\ \tau) \exp(-i\omega t)dt \right|^2\tag{19}$$

$$I\_{FORG}(\omega,\ \tau) = \left| \bigcup\_{\omega}^{\ast \ast} \mathcal{E}(T)E(T-\tau)\exp(-i\omega t)dt \right|^2\tag{20}$$

$$\text{Im}\,\,^2(\omega) = 1 + \sum\_{j=1}^{m} \frac{\left\|{\,^B\omega^j}\_j\right\|^2}{\left\|{\,^B\omega^j}\_j - {\,^B\omega^j}\_j\right\|}\tag{21}$$

$$B\_1 = 0.6961663, \quad \lambda\_1 = 0.0684043 \text{ }\mu m$$

$$B\_2 = 0.4079426, \quad \lambda\_2 = 0.1162414 \text{ }\mu m$$

$$B\_3 = 0.8974794, \quad \lambda\_3 = 9.8961610 \text{ }\mu m$$

$$\beta(\omega) = n(\omega)\frac{\omega}{c} = \beta\_0 + \beta\_1(\omega \text{ - } \omega\_0) + \frac{1}{2}\beta\_2(\omega \text{ - } \omega\_0)^2 + \dots \tag{22}$$

$$\, \_\beta \beta\_k = \left( \frac{\partial^k \beta(\omega)}{\partial \omega^{\,k}} \right)\_{\omega = \omega\_0} \tag{23}$$

The second dispersion coefficient *β*2 accounts for the group velocity dispersion (GVD), mainly leading to pulse broadening. In fact, another dispersion parameter D is usually used to describe the dispersive effect.

$$\mathbf{D} = -\frac{2\pi c}{\lambda^2} \beta\_2 = -\frac{\lambda}{c} \frac{\partial^2 u}{\partial \lambda^2} \tag{24}$$

Dispersive effect in a PCF is associated with both the material dispersion and waveguide dispersion while it mainly arises from material dispersion in bulk materials. Therefore, in terms of dispersion in a PCF, the fiber structure and index distribution should be both taken into account. Considering a solid-core PCF made of fused silica, the finite element method (FEM) is used to model the fundamental mode distribution for different input wavelengths [52-54]. In the modeling, the air hole diameter *d* and pitch Λ are set to be 1.4 µm and 1.6 µm respectively [6], leading to the cross-section profile of the PCF as shown in figure 3(a). Figure 3(b) shows a two-dimension distribution of the fundamental mode for input wavelength λ=0.4 µm with the material dispersion included. Based on this simulation, the effective index of the PCF *neff* is about 1.462353.

Fig.3 (a) Cross-section model of the PCF; (b) Two-dimension distribution of the fundamental mode **Figure 3.** (a) Cross-section model of the PCF; (b) Two-dimension distribution of the fundamental mode

By calculating every ���� corresponding to different input wavelength, the variation of effective index with wavelength is fitted using MATLAB. Then the GVD parameter �� and dispersion parameter D are obtained with Eq. (24) and the profiles with wavelength are plotted in figure 4 respectively. Obviously the ZDW locates at about 730 nm and there is only one ZDW in the range from visible light to near-infrared. By calculating every *neff* corresponding to different input wavelength, the variation of effective index with wavelength is fitted using MATLAB. Then the GVD parameter *β*2 and dispersion parameter D are obtained with Eq. (24) and the profiles with wavelength are plotted in figure 4 respectively. Obviously the ZDW locates at about 730 nm and there is only one ZDW in the range from visible light to near-infrared.

#### (a) (b) *2.3.2. SC generation with PCF*

2.3.2 SC generation with PCF

����� � �0 fs, neglecting the initial chirp.

In this section, we numerically simulate the SC generation process in the PCF described above based on NLSE in Eq. (11) by PC-SSFM. In the simulation, the more accurate expression for nonlinear response function as in Eq. (18) is used. We first consider the SC generation with

Fig.4 Variation of GVD parameter�� (a) and dispersion parameter D (b) with wavelength

Assume the initial input pulse has a hyperbolic scant profile that could be expressed as

with� being the peak power, �� the pulse duration at 1 �⁄ -intensity point, � the initial chirp parameter. It is worth noting that the full width at half maximum (FWHM) is 2�n�1 � √2� times

The temporal and spectral evolutions of input pulse with propagation distance are depicted in

as large as�� . In this simulation, the peak power and pulse duration are set as � � 10 �W ,

A�0� �� � √�s��� � �

��

� ��� �� ����

���

� � (25)

10��fs�/nm、�� � 2.�4� � 10��fs�/nm、�� � �4.�0� � 10��fs�/nm.

In this section, we numerically simulate the SC generation process in the PCF described above based on NLSE in Eq. (11) by PC-SSFM. In the simulation, the more accurate expression for nonlinear response function as in Eq. (18) is used. We first consider the SC generation with PCF described above with hole diameter 1.4 μm and hole pitch 1.6 μm. Because of the ZDW about 730 nm, pump pulses with central wavelength of 800 nm are used to launch the fiber where the nonlinear parameter is estimated to be0.113 W��m�� . And the Taylor series expansion coefficients are �� � �1.��1 � 10��fs�/nm、�� � 7.�1� � 10��fs�/nm、�� � �1.01� �

figure 4 respectively. Obviously the ZDW locates at about 730 nm and there is only one ZDW in the range from visible light to near-infrared. Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 37

Fig.3 (a) Cross-section model of the PCF; (b) Two-dimension distribution of the fundamental mode

Λ

(a) (b)

�

effective index with wavelength is fitted using MATLAB. Then the GVD parameter �� and dispersion parameter D are obtained with Eq. (24) and the profiles with wavelength are plotted in

By calculating every ���� corresponding to different input wavelength, the variation of

Fig.4 Variation of GVD parameter�� (a) and dispersion parameter D (b) with wavelength **Figure 4.** Variation of GVD parameter *β*2 (a) and dispersion parameter D (b) with wavelength

The second dispersion coefficient *β*2 accounts for the group velocity dispersion (GVD), mainly leading to pulse broadening. In fact, another dispersion parameter D is usually used to describe

> *<sup>λ</sup>* <sup>2</sup> *<sup>β</sup>*<sup>2</sup> <sup>=</sup> - *<sup>λ</sup> c* ∂<sup>2</sup> *n*

Dispersive effect in a PCF is associated with both the material dispersion and waveguide dispersion while it mainly arises from material dispersion in bulk materials. Therefore, in terms of dispersion in a PCF, the fiber structure and index distribution should be both taken into account. Considering a solid-core PCF made of fused silica, the finite element method (FEM) is used to model the fundamental mode distribution for different input wavelengths [52-54]. In the modeling, the air hole diameter *d* and pitch Λ are set to be 1.4 µm and 1.6 µm respectively [6], leading to the cross-section profile of the PCF as shown in figure 3(a). Figure 3(b) shows a two-dimension distribution of the fundamental mode for input wavelength λ=0.4 µm with the material dispersion included. Based on this simulation, the effective index of the PCF *neff* is

Fig.3 (a) Cross-section model of the PCF; (b) Two-dimension distribution of the fundamental mode

**Figure 3.** (a) Cross-section model of the PCF; (b) Two-dimension distribution of the fundamental mode

effective index with wavelength is fitted using MATLAB. Then the GVD parameter �� and dispersion parameter D are obtained with Eq. (24) and the profiles with wavelength are plotted in figure 4 respectively. Obviously the ZDW locates at about 730 nm and there is only one ZDW in

In this section, we numerically simulate the SC generation process in the PCF described above based on NLSE in Eq. (11) by PC-SSFM. In the simulation, the more accurate expression for nonlinear response function as in Eq. (18) is used. We first consider the SC generation with

By calculating every *neff* corresponding to different input wavelength, the variation of effective index with wavelength is fitted using MATLAB. Then the GVD parameter *β*2 and dispersion parameter D are obtained with Eq. (24) and the profiles with wavelength are plotted in figure 4 respectively. Obviously the ZDW locates at about 730 nm and there is only one ZDW in the

Fig.4 Variation of GVD parameter�� (a) and dispersion parameter D (b) with wavelength

Assume the initial input pulse has a hyperbolic scant profile that could be expressed as

with� being the peak power, �� the pulse duration at 1 �⁄ -intensity point, � the initial chirp parameter. It is worth noting that the full width at half maximum (FWHM) is 2�n�1 � √2� times

The temporal and spectral evolutions of input pulse with propagation distance are depicted in

as large as�� . In this simulation, the peak power and pulse duration are set as � � 10 �W ,

A�0� �� � √�s��� � �

��

� ��� �� ����

���

� � (25)

10��fs�/nm、�� � 2.�4� � 10��fs�/nm、�� � �4.�0� � 10��fs�/nm.

In this section, we numerically simulate the SC generation process in the PCF described above based on NLSE in Eq. (11) by PC-SSFM. In the simulation, the more accurate expression for nonlinear response function as in Eq. (18) is used. We first consider the SC generation with PCF described above with hole diameter 1.4 μm and hole pitch 1.6 μm. Because of the ZDW about 730 nm, pump pulses with central wavelength of 800 nm are used to launch the fiber where the nonlinear parameter is estimated to be0.113 W��m�� . And the Taylor series expansion coefficients are �� � �1.��1 � 10��fs�/nm、�� � 7.�1� � 10��fs�/nm、�� � �1.01� �

By calculating every ���� corresponding to different input wavelength, the variation of

<sup>∂</sup> *<sup>λ</sup>* <sup>2</sup> (24)

D= - <sup>2</sup>*π<sup>c</sup>*

Λ

(a) (b)

�

the dispersive effect.

36 Photonic Crystals

about 1.462353.

the range from visible light to near-infrared.

*2.3.2. SC generation with PCF*

range from visible light to near-infrared.

(a) (b)

2.3.2 SC generation with PCF

����� � �0 fs, neglecting the initial chirp.

PCF described above with hole diameter 1.4 µm and hole pitch 1.6 µm. Because of the ZDW about 730 nm, pump pulses with central wavelength of 800 nm are used to launch the fiber where the nonlinear parameter is estimated to be 0.113 W-1 m-1 . And the Taylor series expansion coefficients are *β*<sup>2</sup> <sup>=</sup> - 1.581×10-5 fs<sup>2</sup> / nm, *β*<sup>3</sup> =7.819×10-5 fs<sup>3</sup> / nm, *β*<sup>4</sup> <sup>=</sup> - 1.015×10-4 fs<sup>4</sup> / nm, *<sup>β</sup>*<sup>5</sup> =2.549×10-4 fs<sup>5</sup> / nm, *β*<sup>6</sup> <sup>=</sup> - 4.808×10-6 fs<sup>6</sup> / nm. 2.3.2 SC generation with PCF In this section, we numerically simulate the SC generation process in the PCF described above based on NLSE in Eq. (11) by PC-SSFM. In the simulation, the more accurate expression for nonlinear response function as in Eq. (18) is used. We first consider the SC generation with PCF described above with hole diameter 1.4 μm and hole pitch 1.6 μm. Because of the ZDW

Assume the initial input pulse has a hyperbolic scant profile that could be expressed as about 730 nm, pump pulses with central wavelength of 800 nm are used to launch the fiber where

$$\mathbf{A}(0, t) = \sqrt{P} \operatorname{sech}\left(\frac{t}{t\_0}\right) \exp\left(-\frac{\kappa C t^2}{2t\_0^2}\right) \tag{25}$$

with*P* being the peak power, *t*0 the pulse duration at 1 /*e*-intensity point, *C* the initial chirp parameter. It is worth noting that the full width at half maximum (FWHM) is 2ln(1 + 2) times as large as *t*0. In this simulation, the peak power and pulse duration are set as *P* =10 kW , *tFWHM* =50 fs, neglecting the initial chirp. 10��fs�/nm、�� � 2.�4� � 10��fs�/nm、�� � �4.�0� � 10��fs�/nm. Assume the initial input pulse has a hyperbolic scant profile that could be expressed as A�0� �� � √�s��� � � �� � ��� �� ���� ��� � � (25)

The temporal and spectral evolutions of input pulse with propagation distance are depicted in figure 5. In the initial stage of propagation, the spectral broadening is almost symmetric. And the broadening mainly occurs within this first 6 mm, during which the pulse is com‐ pressed strongly. After that, the spectral broadening becomes strong asymmetric, with distinct peaks on the long and short wavelength sides. The long wavelength components manifest the breakup of input pulse into several sub-pulses, known as soliton fission [52] caused by highorder dispersion, self-steepening effect and Raman scattering. The short wavelength compo‐ nents are related to SPM effect and dispersion wave generation [55,56]. with� being the peak power, �� the pulse duration at 1 �⁄ -intensity point, � the initial chirp parameter. It is worth noting that the full width at half maximum (FWHM) is 2�n�1 � √2� times as large as�� . In this simulation, the peak power and pulse duration are set as � � 10 �W , ����� � �0 fs, neglecting the initial chirp. The temporal and spectral evolutions of input pulse with propagation distance are depicted in

With longer propagation, the SC shows larger spectral broadening mainly on the long wavelength side due to soliton self-frequency shift (SSFS) induced by stimulated Raman scattering, while the spectral broadening on the short wavelength side almost remains unchanged. The spectral broadening process is often accompanied with re-distribution of dispersion wave generation[55,56].

figure 5. In the initial stage of propagation, the spectral broadening is almost symmetric. And the broadening mainly occurs within this first 6 mm, during which the pulse is compressed strongly. After that, the spectral broadening becomes strong asymmetric, with distinct peaks on the long and short wavelength sides. The long wavelength components manifest the breakup of input pulse into several sub-pulses, known as soliton fission[52] caused by high-order dispersion, self-steepening

Fig.5 Temporal (left) and spectral (right) evolution over a propagation distance for input pulse centered at 800 nm with peak power 10kW and pulse duration 50fs **Figure 5.** Temporal (left) and spectral (right) evolution over a propagation distance for input pulse centered at 800 nm with peak power 10kW and pulse duration 50fs

energy, that the long wavelength components possess more while the short wavelength components possess less. This is explained by that: With the blue-shift components as pumping light, Raman gain will amplify the red-shift components effectively, leading to energy's transferring from blue components to red components. This energy transfer describes the red shift of the soliton spectrum with propagation distance. With longer propagation, the SC shows larger spectral broadening mainly on the long wavelength side due to soliton self-frequency shift (SSFS) induced by stimulated Raman scattering, while the spectral broadening on the short wavelength side almost remains unchanged. The spectral broadening process is often accompanied with re-distribution of energy, that the long

To study the detailed characteristics of the generated SC, temporal and spectral slices at 15 cm propagation distance are plotted in figure 6. The time structure and spectrum are both shown on a linear scale and the intensity is normalized. There are three main peaks in the time domain located at 0.46 ps, 1.2 ps, and 3.8 ps, corresponding to the distinct peaks at 870 nm, 953 nm, and 1143 nm respectively in the frequency domain. These three peaks are caused by soliton fission and their central positions in both time and frequency domain are determined by their intrinsic characteristics, SSFS effect and the propagation distance. wavelength components possess more while the short wavelength components possess less. This is explained by that: With the blue-shift components as pumping light, Raman gain will amplify the red-shift components effectively, leading to energy's transferring from blue components to red components. This energy transfer describes the red shift of the soliton spectrum with propagation distance. To study the detailed characteristics of the generated SC, temporal and spectral slices at 15

cm propagation distance are plotted in figure 6. The time structure and spectrum are both shown

separation in spectrum at delay time 1.2 ps. There are still some fine oscillating structures lower than the two peaks in figure (a) ranging from 1.4 ps to 1.8 ps. There exists a little bump ranging from 4.6 ps to 6.5 ps in figure (a) with figure (b) being the magnification of partial structure. This

Fig.7 (a) Oscillating structure ranging from 1.2 ps to 1.7ps (b) Generated dispersion wave ranging from 4.6

370 nm

Coherent anti-stokes Raman scattering (CARS) has proven to be a promising nonlinear

The full quantum picture of CARS process is depicted in Fig.8 (A). |��and|��correspond to the vibration excited and ground state respectively. For simplicity of analysis, a single frequency model is provided. All laser fields-the pump, Stokes, probe, and phonons (or molecules) are quantized. The whole CARS process can be divided into two individual steps as shown in Fig.8. In the first step, a coherent phonon population is directly created by coherent Raman scattering when a pump beam at frequency�� and a Stokes beam at frequency �� with frequency difference �� � �� tuned to a Raman-active molecular vibration �� arrive at the sample simultaneously via a high numerical aperture objective. The process is impulsively stimulated Raman scattering (ISRS), which is considered as the main mechanism of coherent phonon generation[60,61]. The coherent phonons are bosons and do not obey an exclusion principle. They

optical technique that is capable of obtaining high-sensitivity and three dimensional images based on molecular vibrational spectrum without labeling[57,58]. CARS is an example of a four-photon

process. The first two light-matter interactions, one with pump frequency ω�and another with Stokes frequency ω�, set to that of a vibrational resonance, a pump photon is annihilated, correspondingly a Stokes photon and a phonon (i.e., the energy quanta of the material excitation)is generated, the number of phonons is equal to the generated Stokes photons. The second light-matter interactions, the energy quanta of the material excitation, phonons, interact with the probe photons to generate a coherent anti-Stokes signal. In general, both processes occur simultaneously and the whole process can be treated as a four-photon process[59]. The detected

ps to 6.5 ps (c) Spectrogram with 15 cm propagation distance

**3 Supercontinuum used in coherent anti-Stokes Raman scattering nanoscopy** 

**3.1 Theories of CARS process** 

signal strength depends nonlinearly on incident intensity.

Fig.6 Pulse characteristics in (a) time and (b) frequency domain after 15 cm propagation **Figure 6.** Pulse characteristics in (a) time and (b) frequency domain after 15 cm propagation between the two sidelobes. The relation of them is easy to see in figure 7(c) with 370 nm

(a) (b) (c)

bump is related to the dispersion wave generation around 546 nm.

To observe the detailed structure in figure 6(a), magnified profiles of the segments ranging from 1.1 ps to 1.8ps, and 4ps to 7ps are plotted in figure 7(a) and (b) respectively. In figure (a), the ultrafast oscillating structure is associated with two sidelobes on both side of input wavelength. The main peak belongs to a second soliton and the oscillating structure results from beating between the two sidelobes. The relation of them is easy to see in figure 7(c) with 370 nm separation in spectrum at delay time 1.2 ps. There are still some fine oscillating structures lower than the two peaks in figure (a) ranging from 1.4 ps to 1.8 ps. There exists a little bump ranging from 4.6 ps to 6.5 ps in figure (a) with figure (b) being the magnification of partial structure. This bump is related to the dispersion wave generation around 546 nm.

Fig.6 Pulse characteristics in (a) time and (b) frequency domain after 15 cm propagation

Fig.7 (a) Oscillating structure ranging from 1.2 ps to 1.7ps (b) Generated dispersion wave ranging from 4.6 ps to 6.5 ps (c) Spectrogram with 15 cm propagation distance **Figure 7.** (a) Oscillating structure ranging from 1.2 ps to 1.7ps (b) Generated dispersion wave ranging from 4.6 ps to 6.5 ps (c) Spectrogram with 15 cm propagation distance

#### **3 Supercontinuum used in coherent anti-Stokes Raman scattering nanoscopy 3.1 Theories of CARS process 3. Supercontinuum used in coherent anti-Stokes Raman scattering nanoscopy**

#### Coherent anti-stokes Raman scattering (CARS) has proven to be a promising nonlinear **3.1. Theories of CARS process**

energy, that the long wavelength components possess more while the short wavelength components possess less. This is explained by that: With the blue-shift components as pumping light, Raman gain will amplify the red-shift components effectively, leading to energy's transferring from blue components to red components. This energy transfer describes the red

**Figure 5.** Temporal (left) and spectral (right) evolution over a propagation distance for input pulse centered at 800 nm

With longer propagation, the SC shows larger spectral broadening mainly on the long

Fig.5 Temporal (left) and spectral (right) evolution over a propagation distance for input pulse centered at 800 nm with peak power 10kW and pulse duration 50fs

wavelength side due to soliton self-frequency shift (SSFS) induced by stimulated Raman scattering, while the spectral broadening on the short wavelength side almost remains unchanged. The spectral broadening process is often accompanied with re-distribution of energy, that the long wavelength components possess more while the short wavelength components possess less. This is explained by that: With the blue-shift components as pumping light, Raman gain will amplify the red-shift components effectively, leading to energy's transferring from blue components to red components. This energy transfer describes the red shift of the soliton spectrum with propagation

figure 5. In the initial stage of propagation, the spectral broadening is almost symmetric. And the broadening mainly occurs within this first 6 mm, during which the pulse is compressed strongly. After that, the spectral broadening becomes strong asymmetric, with distinct peaks on the long and short wavelength sides. The long wavelength components manifest the breakup of input pulse into several sub-pulses, known as soliton fission[52] caused by high-order dispersion, self-steepening effect and Raman scattering. The short wavelength components are related to SPM effect and

To study the detailed characteristics of the generated SC, temporal and spectral slices at 15 cm propagation distance are plotted in figure 6. The time structure and spectrum are both shown on a linear scale and the intensity is normalized. There are three main peaks in the time domain located at 0.46 ps, 1.2 ps, and 3.8 ps, corresponding to the distinct peaks at 870 nm, 953 nm, and 1143 nm respectively in the frequency domain. These three peaks are caused by soliton fission and their central positions in both time and frequency domain are determined by their

To study the detailed characteristics of the generated SC, temporal and spectral slices at 15

cm propagation distance are plotted in figure 6. The time structure and spectrum are both shown on a linear scale and the intensity is normalized. There are three main peaks in the time domain located at 0.46 ps, 1.2 ps, and 3.8 ps, corresponding to the distinct peaks at 870 nm, 953 nm, and 1143 nm respectively in the frequency domain. These three peaks are caused by soliton fission and their central positions in both time and frequency domain are determined by their intrinsic

Fig.6 Pulse characteristics in (a) time and (b) frequency domain after 15 cm propagation

To observe the detailed structure in figure 6(a), magnified profiles of the segments ranging

from 1.1 ps to 1.8ps, and 4ps to 7ps are plotted in figure 7(a) and (b) respectively. In figure (a), the ultrafast oscillating structure is associated with two sidelobes on both side of input wavelength. The main peak belongs to a second soliton and the oscillating structure results from beating between the two sidelobes. The relation of them is easy to see in figure 7(c) with 370 nm separation in spectrum at delay time 1.2 ps. There are still some fine oscillating structures lower than the two peaks in figure (a) ranging from 1.4 ps to 1.8 ps. There exists a little bump ranging from 4.6 ps to 6.5 ps in figure (a) with figure (b) being the magnification of partial structure. This

Fig.7 (a) Oscillating structure ranging from 1.2 ps to 1.7ps (b) Generated dispersion wave ranging from 4.6

370 nm

Coherent anti-stokes Raman scattering (CARS) has proven to be a promising nonlinear

The full quantum picture of CARS process is depicted in Fig.8 (A). |��and|��correspond to the vibration excited and ground state respectively. For simplicity of analysis, a single frequency model is provided. All laser fields-the pump, Stokes, probe, and phonons (or molecules) are quantized. The whole CARS process can be divided into two individual steps as shown in Fig.8. In the first step, a coherent phonon population is directly created by coherent Raman scattering when a pump beam at frequency�� and a Stokes beam at frequency �� with frequency difference �� � �� tuned to a Raman-active molecular vibration �� arrive at the sample simultaneously via a high numerical aperture objective. The process is impulsively stimulated Raman scattering (ISRS), which is considered as the main mechanism of coherent phonon generation[60,61]. The coherent phonons are bosons and do not obey an exclusion principle. They

optical technique that is capable of obtaining high-sensitivity and three dimensional images based on molecular vibrational spectrum without labeling[57,58]. CARS is an example of a four-photon

process. The first two light-matter interactions, one with pump frequency ω�and another with Stokes frequency ω�, set to that of a vibrational resonance, a pump photon is annihilated, correspondingly a Stokes photon and a phonon (i.e., the energy quanta of the material excitation)is generated, the number of phonons is equal to the generated Stokes photons. The second light-matter interactions, the energy quanta of the material excitation, phonons, interact with the probe photons to generate a coherent anti-Stokes signal. In general, both processes occur simultaneously and the whole process can be treated as a four-photon process[59]. The detected

ps to 6.5 ps (c) Spectrogram with 15 cm propagation distance

**3 Supercontinuum used in coherent anti-Stokes Raman scattering nanoscopy** 

**3.1 Theories of CARS process** 

signal strength depends nonlinearly on incident intensity.

shift of the soliton spectrum with propagation distance.

with peak power 10kW and pulse duration 50fs

(a) (b)

dispersion wave generation[55,56].

38 Photonic Crystals

distance.

intrinsic characteristics, SSFS effect and the propagation distance.

(a) (b)

characteristics, SSFS effect and the propagation distance.

(a) (b) (c)

bump is related to the dispersion wave generation around 546 nm.

**Figure 6.** Pulse characteristics in (a) time and (b) frequency domain after 15 cm propagation

optical technique that is capable of obtaining high-sensitivity and three dimensional images based on molecular vibrational spectrum without labeling[57,58]. CARS is an example of a four-photon process. The first two light-matter interactions, one with pump frequency ω�and another with Stokes frequency ω�, set to that of a vibrational resonance, a pump photon is annihilated, correspondingly a Stokes photon and a phonon (i.e., the energy quanta of the material excitation)is generated, the number of phonons is equal to the generated Stokes photons. The second light-matter interactions, the energy quanta of the material excitation, phonons, interact with the probe photons to generate a coherent anti-Stokes signal. In general, both processes occur simultaneously and the whole process can be treated as a four-photon process[59]. The detected Coherent anti-stokes Raman scattering (CARS) has proven to be a promising nonlinear optical technique that is capable of obtaining high-sensitivity and three dimensional images based on molecular vibrational spectrum without labeling [57,58]. CARS is an example of a four-photon process. The first two light-matter interactions, one with pump frequency ω*L* and another with Stokes frequency ω*S* , set to that of a vibrational resonance, a pump photon is annihilated, correspondingly a Stokes photon and a phonon (i.e., the energy quanta of the material excitation)is generated, the number of phonons is equal to the generated Stokes photons. The second light-matter interactions, the energy quanta of the material excitation, phonons, interact with the probe photons to generate a coherent anti-Stokes signal. In general, both processes occur simultaneously and the whole process can be treated as a four-photon process [59]. The detected signal strength depends nonlinearly on incident intensity.

signal strength depends nonlinearly on incident intensity. The full quantum picture of CARS process is depicted in Fig.8 (A). |��and|��correspond to the vibration excited and ground state respectively. For simplicity of analysis, a single frequency The full quantum picture of CARS process is depicted in Fig.8. |*e* and |*g* correspond to the vibration excited and ground state respectively. For simplicity of analysis, a single frequency model is provided. All laser fields-the pump, Stokes, probe, and phonons (or molecules) are quantized. The whole CARS process can be divided into two individual steps as shown in Fig.

model is provided. All laser fields-the pump, Stokes, probe, and phonons (or molecules) are quantized. The whole CARS process can be divided into two individual steps as shown in Fig.8. In the first step, a coherent phonon population is directly created by coherent Raman scattering when a pump beam at frequency�� and a Stokes beam at frequency �� with frequency difference �� � �� tuned to a Raman-active molecular vibration �� arrive at the sample simultaneously via a high numerical aperture objective. The process is impulsively stimulated Raman scattering (ISRS), which is considered as the main mechanism of coherent phonon generation[60,61]. The coherent phonons are bosons and do not obey an exclusion principle. They 8. In the first step, a coherent phonon population is directly created by coherent Raman scattering when a pump beam at frequency *ωL* and a Stokes beam at frequency *ω<sup>S</sup>* with frequency difference *ω<sup>L</sup>* - *ωS* tuned to a Raman-active molecular vibration *ω<sup>v</sup>* arrive at the sample simultaneously via a high numerical aperture objective. The process is impulsively stimulated Raman scattering (ISRS), which is considered as the main mechanism of coherent phonon generation [60,61]. The coherent phonons are bosons and do not obey an exclusion principle. They are in a non-equilibrium state and in phase. In the second step, the coherent phonons interact with a time-delayed probe beam at frequency *ω<sup>P</sup>* to generate a blue-shift anti-Stokes signal at ω*AS* =ω*<sup>P</sup>* + ω*v* , leading to separation of the CARS signal from the incident laser beams conveniently and efficiently rejecting fluorescence. Delay time between coherent phonons and probe beam is certainly much shorter than the vibrational dephasing time.

**Figure 8.** Energy level diagram of the two individual steps in CARS

The quantum mechanical treatment of the first process, which is a second-order process involving two light-matter interactions, resembles the spontaneous Raman effect. The important difference of the two is that in CARS process, the Stokes light stems from an applied laser field. Spontaneous Raman is a weak effect because the spontaneous interaction through the vacuum field occurs only rarely. This weakness can be overcome when the spontaneous nature of the initial state to final state *i* → *f* transition is eliminated by applying a second field of frequency *ωs*. The treatment is based on the calculation of the transition rate between the states of the molecule, which is described by Fermi's golden rule. The equation (26) predicts the growth of the generated Stokes photons, i.e., phonon number [62,63]

Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 41

$$\frac{d\ln\_{\rm S}^{\rm inc}}{d\mathbf{l}} = \mathbf{N} \bullet \frac{\pi e^{4} \omega\_{L} \cdot \omega\_{\rm S}}{2\varepsilon\_{0}^{2} \hbar^{2} V^{2} \mu\_{L}^{2} \mu\_{\rm S} c} \mid \alpha\_{R} \mid {}^{2}\boldsymbol{n}\_{L}^{\*} \left[\boldsymbol{n}\_{\rm S}^{\*} + \mathbf{1}\right] \delta(\omega\_{\rm S} \cdot \omega\_{\rm L} + \omega\_{\rm v}) \tag{26}$$

Consider a volume *V* , which contains *N* molecules and has a length *l* parallel to the beam propagation direction. The Raman transition polarizability *αR* plays a role similar to the classical polarizability change (∂*α* / ∂*q*). *ω<sup>L</sup>* and *ωS* is the frequencies of pump and Stokes light, respectively. *e* is the electron charge, *ε*0 is the vacuum permittivity, ℏ is Planck's constant, and *δ* is the Dirac delta function.*μL* and *μS* denote the refractive indexes. *nS inc* indicates the increased (decreased) number of Stokes (pump) photons, i.e., the phonon number of energy ℏ*ωv* , *nphonon*. In the CARS process, the Stokes light stems from an applied laser field, and the *nS \** consists of the applied Stokes photon *nS* and increased number *nS inc* that stems from the first two light-matter interaction process. Similarly, *nL \** , amounts to the pump photon number *nL* subtracting the decreased number *nS inc* . With the use of the pump and Stokes laser light sources of sufficiently high power, *nS inc*≪*nL* , *nS* ,the present pump and Stokes number *nL \** <sup>=</sup>*nL* - *nS inc* , *nS \** <sup>=</sup>*nS* <sup>+</sup> *nS inc* , approximately amount to *nL* and *nS* , and the factor 1 in equation (26) can be omitted at the same time. The laser intensities *I <sup>L</sup>* and *I <sup>S</sup>* exist on the focal plane in the form of Gaussian distribution, and we adopt an infinitesimal region in the focal plane where the intensities can be supposed constant *I <sup>L</sup>* =(ℏ *ω<sup>L</sup> c* / *μL V* )*nL* , *IS* =(ℏ *ω<sup>S</sup> c* / *μSV* )*nS* . In an fs-CARS transient, the pump and Stokes beams temporally overlapped, coherently excite the molecular field, the duration of the transitory process denoted by the coherent excitation time is negligible. The homogeneous broadening of the atomic final state converts the delta function to a Lorentzian lineshape. when *ω<sup>S</sup>* <sup>=</sup>*ω<sup>L</sup>* - *<sup>ω</sup>v*, *<sup>δ</sup>*(*ω<sup>S</sup>* - *<sup>ω</sup><sup>L</sup>* <sup>+</sup> *<sup>ω</sup>v*)<sup>→</sup> <sup>Γ</sup> / *<sup>π</sup>* (*ω<sup>S</sup>* - *<sup>ω</sup> <sup>L</sup>* <sup>+</sup> *<sup>ω</sup>v*)2 <sup>+</sup> <sup>Γ</sup><sup>2</sup> <sup>=</sup> <sup>1</sup> *<sup>π</sup>*<sup>Γ</sup> , where Γ denotes the linewidth relating to the dephasing time *T*2 of the physical system, Γ=1 / *T*2.

Accordingly, from the equation (26), we obtain:

8. In the first step, a coherent phonon population is directly created by coherent Raman scattering when a pump beam at frequency *ωL* and a Stokes beam at frequency *ω<sup>S</sup>* with frequency difference *ω<sup>L</sup>* - *ωS* tuned to a Raman-active molecular vibration *ω<sup>v</sup>* arrive at the sample simultaneously via a high numerical aperture objective. The process is impulsively stimulated Raman scattering (ISRS), which is considered as the main mechanism of coherent phonon generation [60,61]. The coherent phonons are bosons and do not obey an exclusion principle. They are in a non-equilibrium state and in phase. In the second step, the coherent phonons interact with a time-delayed probe beam at frequency *ω<sup>P</sup>* to generate a blue-shift anti-Stokes signal at ω*AS* =ω*<sup>P</sup>* + ω*v* , leading to separation of the CARS signal from the incident laser beams conveniently and efficiently rejecting fluorescence. Delay time between coherent phonons and probe beam is certainly much shorter than the vibrational dephasing time.

40 Photonic Crystals

**Figure 8.** Energy level diagram of the two individual steps in CARS

The quantum mechanical treatment of the first process, which is a second-order process involving two light-matter interactions, resembles the spontaneous Raman effect. The important difference of the two is that in CARS process, the Stokes light stems from an applied laser field. Spontaneous Raman is a weak effect because the spontaneous interaction through the vacuum field occurs only rarely. This weakness can be overcome when the spontaneous nature of the initial state to final state *i* → *f* transition is eliminated by applying a second field of frequency *ωs*. The treatment is based on the calculation of the transition rate between the states of the molecule, which is described by Fermi's golden rule. The equation (26) predicts

the growth of the generated Stokes photons, i.e., phonon number [62,63]

$$m\_S^{inc} = \frac{I\_S \mu\_S V}{\hbar \omega\_S \cdot c} \bullet I\_L \text{ g}l \tag{27}$$

where *g* is the gain coefficient, *N*0 is the molecule density.

$$\log g = \mathcal{N}\_0 \bullet \frac{e^4}{2\epsilon\_0^2 \hbar^3 V \,^2 \mu\_L^2 \, ^2 \mu\_S \, c^2} \frac{\omega\_\mathcal{S}}{\mu\_L \, \mu\_\mathcal{S}} \parallel \, \alpha\_\mathcal{R} \parallel \stackrel{2}{\mathcal{T}} \tag{28}$$

Then, considering a molecule with just a single vibrational transition, the coherent phonon number per time unit in the first ISRS step can be written as

$$
\mu\_v = \frac{I\_s \,\mu\_s \, V}{\hbar \omega\_s \, c} \bullet I\_L \, gl \, \tag{29}
$$

where *nv* indicates the coherent phonon number. The interaction volume V, which contains N molecules, is defined by the diffraction-limited focal spot area *Afoc* and the interaction length *l* paralleling to the beam propagation direction. For g∝ *N*0 , it is apparent that the coherent phonon number is linearly dependent on *N*0.

The phonon number subjects to a saturation effect and will not increase infinitely with the growth of incident light powers [64]. The maximum of coherent phonon number per pulse is molecule number in the focal volume. In another word, all of the molecules in the focus are excited to the vibrational state.

$$m\_v^{pulse} = \frac{I\_S \mu\_S \, V}{\hbar \omega\_S \, \text{c}} \bullet I\_L \text{ g} \tau \le \text{N} \tag{30}$$

Here, *nv pulse* is the coherent phonon number per pulse. The magnitude of the spot is roughly estimated as ~10-9 cm<sup>2</sup> , corresponding to a focal volume *V* ~ 8×10-14 cm<sup>3</sup> . The gain coefficient is *<sup>g</sup>* <sup>~</sup> *<sup>N</sup>*<sup>0</sup> ×10-31cm / W. From equation (30), we see that saturation of coherent phonons results in product of the pump and Stokes intensities not exceeding ~1018W2 / cm<sup>4</sup> . Therefore in the first step of a T-CARS process, the pump and Stokes power density at focus is on the same order of ~10<sup>9</sup> W / cm<sup>2</sup> . In principle, if the pump and Stokes intensity exceed such an extreme value, the further increased intensity will not contribute to generation of the coherent phonons any more. In practice, however, energy of the pump and Stokes may be higher than that theoretical estimation due to the spectral broadening of femtosecond laser pulses.

Hereafter, we analyze the second step in CARS process. In the spontaneous Raman anti-Stokes process, the involved lattice vibrations called incoherent phonons are thermally excited with a density function *nv* given by the Bose-Einstein distribution. In this case, oscillatory phases of the phonons are completely random. Thus, individual oscillatory motions are canceled out by each other. Spontaneous Raman scattering is an incoherent linear process and the incoher‐ ent signal is randomly scattered to 4*π* solid angle, then the collected signal is [65]

$$P\_A = N \left(\frac{\partial \sigma}{\partial D}\right)\_A \mathbb{1}\_v I\_L \ \mathfrak{Q}\_{collection} \tag{31}$$

with *PA* representing the anti-Stokes power, *I <sup>L</sup>* the incident laser intensity and Ω*collection* the effective collection angle.( <sup>∂</sup> <sup>σ</sup> <sup>∂</sup> <sup>Ω</sup> )*A* is the spontaneous Raman differential scattering cross sec‐ tion. The scattered signal is proportional to the incident laser intensity and the occupation number of the upper Raman level. Different from spontaneous anti-Stokes Raman scatter‐ ing, the second step in CARS is via introducing probe photons to interact with the coherent phonons, which are in a non-equilibrium state and in phase. Delay time between coherent phonons and probe beam is certainly much shorter than the vibrational dephasing time. Therefore, relaxation of the phonons is omitted, that is, all phonons are considered to have chances to interact with the probe photons and the impact of the delay time on the signal strength can be readily omitted. An important advantage of CARS microscopy over sponta‐ neous Raman scattering microscopy is the fact that CARS signal is emitted into a 10-4 *sr*range, comparing to spontaneous Raman signal randomly scattered into 4*π* solid angle, which provides an important information that scattering cross section ∂*σ<sup>C</sup>* for the interaction of co‐ herent phonons and probe photons is estimated to be enhanced by about five orders of mag‐ nitude [66]. In some sense, coherent process is equivalent to increasing the scattering cross section. The coherent anti-Stokes photons are confined to such a small solid angle, thus a condenser can almost completely collect the forward propagating signal. For low concentra‐ tion samples, the distance between adjacent solute molecules is so great that energy of the coherent phonons cannot effectively transmit from one molecule to another. Accordingly, in the second step, the anti-Stokes signal strength is exclusively proportional to the number of coherent phonons and then linear dependence on the concentration of solute molecules. Comparing with incoherent (spherical) spontaneous Raman process, we calculate the num‐ ber of the detected photons for coherent (directional) CARS in the repetition mode as

where *nv* indicates the coherent phonon number. The interaction volume V, which contains N molecules, is defined by the diffraction-limited focal spot area *Afoc* and the interaction length *l* paralleling to the beam propagation direction. For g∝ *N*0 , it is apparent that the coherent

The phonon number subjects to a saturation effect and will not increase infinitely with the growth of incident light powers [64]. The maximum of coherent phonon number per pulse is molecule number in the focal volume. In another word, all of the molecules in the focus are

*pulse* is the coherent phonon number per pulse. The magnitude of the spot is roughly

. In principle, if the pump and Stokes intensity exceed such an extreme

is *<sup>g</sup>* <sup>~</sup> *<sup>N</sup>*<sup>0</sup> ×10-31cm / W. From equation (30), we see that saturation of coherent phonons results

first step of a T-CARS process, the pump and Stokes power density at focus is on the same

value, the further increased intensity will not contribute to generation of the coherent phonons any more. In practice, however, energy of the pump and Stokes may be higher than that

Hereafter, we analyze the second step in CARS process. In the spontaneous Raman anti-Stokes process, the involved lattice vibrations called incoherent phonons are thermally excited with a density function *nv* given by the Bose-Einstein distribution. In this case, oscillatory phases of the phonons are completely random. Thus, individual oscillatory motions are canceled out by each other. Spontaneous Raman scattering is an incoherent linear process and the incoher‐

with *PA* representing the anti-Stokes power, *I <sup>L</sup>* the incident laser intensity and Ω*collection* the

tion. The scattered signal is proportional to the incident laser intensity and the occupation number of the upper Raman level. Different from spontaneous anti-Stokes Raman scatter‐ ing, the second step in CARS is via introducing probe photons to interact with the coherent phonons, which are in a non-equilibrium state and in phase. Delay time between coherent phonons and probe beam is certainly much shorter than the vibrational dephasing time. Therefore, relaxation of the phonons is omitted, that is, all phonons are considered to have chances to interact with the probe photons and the impact of the delay time on the signal strength can be readily omitted. An important advantage of CARS microscopy over sponta‐

<sup>ℏ</sup>*ω<sup>S</sup> <sup>c</sup>* <sup>∙</sup> *<sup>I</sup> <sup>L</sup> <sup>g</sup><sup>τ</sup>* <sup>≤</sup> *<sup>N</sup>* (30)

<sup>∂</sup> *<sup>Ω</sup>* )*Anv I <sup>L</sup> Ωcollection* (31)

<sup>∂</sup> <sup>Ω</sup> )*A* is the spontaneous Raman differential scattering cross sec‐

. The gain coefficient

. Therefore in the

phonon number is linearly dependent on *N*0.

*nv*

*pulse* <sup>=</sup> *ISμSV*

estimated as ~10-9 cm<sup>2</sup> , corresponding to a focal volume *V* ~ 8×10-14 cm<sup>3</sup>

in product of the pump and Stokes intensities not exceeding ~1018W2 / cm<sup>4</sup>

theoretical estimation due to the spectral broadening of femtosecond laser pulses.

ent signal is randomly scattered to 4*π* solid angle, then the collected signal is [65]

*PA* <sup>=</sup> *<sup>N</sup>* ( <sup>∂</sup> *<sup>σ</sup>*

excited to the vibrational state.

W / cm<sup>2</sup>

effective collection angle.( <sup>∂</sup> <sup>σ</sup>

Here, *nv*

42 Photonic Crystals

order of ~10<sup>9</sup>

$$
\ln \mathbf{n}\_{\text{CARS}}^{low} = \mathbf{n}\_v^{pulse} \bullet \mathbf{n}\_P \bullet \frac{1}{A\_{fc}} \partial \sigma\_\mathbb{C} \bullet \Delta t \, f\_{rep} \tag{32}
$$

In contrast, for high concentration samples, the coherent phonon energy can freely transmit among solute molecules in the focal volume, therefore in the second step, anti-Stokes signal strength isn't simply proportional to the coherent phonon number but also concerns with the molecule number. The overall T-CARS process depends quadratically on concentration, hence the anti-Stokes photon number in the repetition mode is

$$m\_{\rm{CARS}}^{high} = \mathbf{N} \bullet m\_v^{pulse} \bullet m\_P \bullet \frac{1}{A\_{\rm{fc}}} \partial \sigma\_\mathbf{C} \bullet \Delta t \, f\_{rep} \tag{33}$$

where, *nCARS low* and *nCARS high* denote anti-Stokes photon number in the low and high concentration samples, respectively. *nP* is photon number of the probe beam per pulse. The three pulsed nearinfrared excitation beams are tightly focused with a 1.2NA water-dipping lens for deeper penetration depth to the spot area *A foc*. ∆*t* is exposure time. *N* is molecule number in the foci. In short, for *nv pulse* <sup>∝</sup> *<sup>N</sup>*0 , the anti-Stokes signal strength is strictly linear and quadratic on concentration in the low and high concentration samples, respectively, which conflicts with the popular statements of CARS signal consistently quadratic concentration dependence. The linear concentration dependence is especially beneficial to quantitative analysis for low concentration samples [67-70]. In the cross area between the low and high concentration, the observed concentration dependence from the relative experiments is neither quadratic nor linear, i.e., *nCARS* <sup>∝</sup> *<sup>N</sup> <sup>α</sup>*, where 1<*<sup>α</sup>* <2. An experimental proof in Figure 9 shows the depend‐ ence of T-CARS signal strength on the concentration of alcohol in pure water at C-H stretch band [66,71] (2900cm-1). We record the signal intensity at every volume fraction by 20 times to exclude the influence of signal fluctuation. When volume fraction of alcohol gradually increases, concentration dependence transits from strictly linear to sub-quadratic, until quadratic. It is apparent that when volume fraction of alcohol in water is lower than 20%, the theoretical analysis.

concentration dependence is strictly linear. We prepare the aqueous solution of alcohol with smaller volume fraction intervals to obtain a more precise fitting result. This experiment result consists with the theoretical analysis, and more different samples will be prepared to verify the universality of the theoretical analysis. dependence is strictly linear. We prepare the aqueous solution of alcohol with smaller volume fraction intervals to obtain a more precise fitting result. This experiment result consists with the theoretical analysis, and more different samples will be prepared to verify the universality of the

concentration dependence from the relative experiments is neither quadratic nor linear,

i.e.,����� � ��, where � � � � �. An experimental proof in figure 9 shows the dependence of T-CARS signal strength on the concentration of alcohol in pure water at C-H stretch band[66,71] (2900cm-1). We record the signal intensity at every volume fraction by 20 times to exclude the

influence of signal fluctuation. When volume fraction of alcohol gradually increases,

apparent that when volume fraction of alcohol in water is lower than 20%, the concentration

Fig.9 T-CARS intensity for increasing volume fractions of alcohol in water at different volume fraction from 7% to 100%. **Figure 9.** T-CARS intensity for increasing volume fractions of alcohol in water at different volume fraction from 7% to 100%.

#### **3.2. Break the lateral diffraction barrier in CARS microscopy**

**3.2 Break the lateral diffraction barrier in CARS microscopy** 

resolution of around 120nm, etc.

The application of CARS microscopy is limited, nevertheless, by the acquirable lateral spatial resolution, typically around 300nm. With anti-Stokes signal generated from the focal spot, which is then scanned to build up a an image of the sample, a promising way to generate sub-diffraction-limited CARS images is via minimizing the spot to nanometer extent, as in STED[20]. Several schemes have been suggested recently for breaking the diffraction barrier in CARS microscopy, and can be mainly classified as three means[72-74]: by use of an additional saturation laser to populate coupled vibrational levels; a phase-controlled local oscillator generating ring field at the anti-Stokes frequency being out of phase with the induced CARS field in the focal volume; to add structured illumination for achieving widefield CARS microscopy. Whereas, these schemes still rest on the theory, and no experimental implementations have been reported. And the flows of these schemes are obvious: it's difficult to implement the schemes and The application of CARS microscopy is limited, nevertheless, by the acquirable lateral spatial resolution, typically around 300nm. With anti-Stokes signal generated from the focal spot, which is then scanned to build up a an image of the sample, a promising way to generate subdiffraction-limited CARS images is via minimizing the spot to nanometer extent, as in STED [20]. Several schemes have been suggested recently for breaking the diffraction barrier in CARS microscopy, and can be mainly classified as three means [72-74]: by use of an additional saturation laser to populate coupled vibrational levels; a phase-controlled local oscillator generating ring field at the anti-Stokes frequency being out of phase with the induced CARS field in the focal volume; to add structured illumination for achieving widefield CARS microscopy. Whereas, these schemes still rest on the theory, and no experimental implemen‐ tations have been reported. And the flows of these schemes are obvious: it's difficult to implement the schemes and can only obtain resolution enhanced images of the molecule specified by the anti-Stokes Raman spectral signal for its particular chemical-bond, not its broadband even total CARS spectral signals; the second method must know the type of the bond in advance, and the third can merely obtain a resolution of around 120nm, etc.

can only obtain resolution enhanced images of the molecule specified by the anti-Stokes Raman spectral signal for its particular chemical-bond, not its broadband even total CARS spectral signals; We suggest an approach, the so-called additional probe-beam-induced phonon depletion (APIPD) method [25], to visualize the resolution enhancement by introducing an additional

the second method must know the type of the bond in advance, and the third can merely obtain a

(APIPD) method[25], to visualize the resolution enhancement by the introduction of an additional

doughnut probe beam with a wavelength that is different from the Gaussian probe beam, which almost synchronizes with the pump and Stokes beams, to deplete the phonons at the periphery of

We suggest an approach, the so-called additional probe-beam-induced phonon depletion

doughnut probe beam with a wavelength that is different from the Gaussian probe beam, which almost synchronizes with the pump and Stokes beams, to deplete the phonons at the periphery of the focal spot. The Gaussian probe beam that immediately follows, with a delay time less than 1ps, will yield anti-Stokes signal in the center but not at the periphery of the focal spot. The difference of the two probe beams in wavelength leads to the disparity of two anti-Stokes signals in spectra, and the signals at the periphery can be filter out by virtue of a particular filter. Consequently, the effective focal spot can be substantially reduced. Superresolution images can be obtained by scanning the suppressed focal spot.

The rationale of our scheme is to suppress the anti-Stokes signal generation at the periphery of the diffraction-limited spot, i.e., deplete the phonons produced by the pump and Stokes light in this region, then no phonon interact with the probe photons to generate CARS signal here anymore. We adopt an additional probe beam, the first one, with frequency ω*P*1 , which is doughnut-mode, analogy to STED, to achieve this by interacting with phonons to generate useless anti-Stokes signal with frequency ω*A*1 at the periphery which can be filtered out through a specific filter. The probe beam (ω*P*1) is nearly simultaneous with the pump and Stokes beams on the sample. The intensity of the first probe beam (ω*P*1) is strong enough to completely deplete the phonons at the peripheric region, while the immediately followed second probe light (ω*P*2) with Gaussian profile only generates the useful anti-Stokes signal (ω*A*2) at the central part of the spot. The time delay between the two probe beams is certainly much smaller than the vibrational dephasing time. The difference of ω*P*1 and ω*P*2 results in the disparity of ω*A*1 and ω*A*2 , therefore we can select the specific filter to filter out the unwanted anti-Stokes signal of frequency ω*A*1 to obtain the useful signal ω*P*2 at the suppressed extent of the spot. The process of the APIPD is illustrated as in figure 10.

**Figure 10.** The illustration of the process of the APIPD method

concentration dependence is strictly linear. We prepare the aqueous solution of alcohol with smaller volume fraction intervals to obtain a more precise fitting result. This experiment result consists with the theoretical analysis, and more different samples will be prepared to verify

Fig.9 T-CARS intensity for increasing volume fractions of alcohol in water at different volume fraction from 7% to 100%.

**Figure 9.** T-CARS intensity for increasing volume fractions of alcohol in water at different volume fraction from 7% to

**Concentration (Val/V)** 

Sub-quadratic fit Quadratic fit

Experiment Linear fit

resolution, typically around 300nm. With anti-Stokes signal generated from the focal spot, which is then scanned to build up a an image of the sample, a promising way to generate sub-diffraction-limited CARS images is via minimizing the spot to nanometer extent, as in STED[20]. Several schemes have been suggested recently for breaking the diffraction barrier in CARS microscopy, and can be mainly classified as three means[72-74]: by use of an additional

The application of CARS microscopy is limited, nevertheless, by the acquirable lateral spatial resolution, typically around 300nm. With anti-Stokes signal generated from the focal spot, which is then scanned to build up a an image of the sample, a promising way to generate subdiffraction-limited CARS images is via minimizing the spot to nanometer extent, as in STED [20]. Several schemes have been suggested recently for breaking the diffraction barrier in CARS microscopy, and can be mainly classified as three means [72-74]: by use of an additional saturation laser to populate coupled vibrational levels; a phase-controlled local oscillator generating ring field at the anti-Stokes frequency being out of phase with the induced CARS field in the focal volume; to add structured illumination for achieving widefield CARS microscopy. Whereas, these schemes still rest on the theory, and no experimental implemen‐ tations have been reported. And the flows of these schemes are obvious: it's difficult to implement the schemes and can only obtain resolution enhanced images of the molecule specified by the anti-Stokes Raman spectral signal for its particular chemical-bond, not its broadband even total CARS spectral signals; the second method must know the type of the

saturation laser to populate coupled vibrational levels; a phase-controlled local oscillator generating ring field at the anti-Stokes frequency being out of phase with the induced CARS field in the focal volume; to add structured illumination for achieving widefield CARS microscopy. Whereas, these schemes still rest on the theory, and no experimental implementations have been reported. And the flows of these schemes are obvious: it's difficult to implement the schemes and can only obtain resolution enhanced images of the molecule specified by the anti-Stokes Raman spectral signal for its particular chemical-bond, not its broadband even total CARS spectral signals; the second method must know the type of the bond in advance, and the third can merely obtain a

We suggest an approach, the so-called additional probe-beam-induced phonon depletion (APIPD) method [25], to visualize the resolution enhancement by introducing an additional

bond in advance, and the third can merely obtain a resolution of around 120nm, etc.

We suggest an approach, the so-called additional probe-beam-induced phonon depletion

(APIPD) method[25], to visualize the resolution enhancement by the introduction of an additional

doughnut probe beam with a wavelength that is different from the Gaussian probe beam, which almost synchronizes with the pump and Stokes beams, to deplete the phonons at the periphery of

The application of CARS microscopy is limited, nevertheless, by the acquirable lateral spatial

concentration dependence from the relative experiments is neither quadratic nor linear,

i.e.,����� � ��, where � � � � �. An experimental proof in figure 9 shows the dependence of T-CARS signal strength on the concentration of alcohol in pure water at C-H stretch band[66,71] (2900cm-1). We record the signal intensity at every volume fraction by 20 times to exclude the

influence of signal fluctuation. When volume fraction of alcohol gradually increases, concentration dependence transits from strictly linear to sub-quadratic, until quadratic. It is apparent that when volume fraction of alcohol in water is lower than 20%, the concentration dependence is strictly linear. We prepare the aqueous solution of alcohol with smaller volume fraction intervals to obtain a more precise fitting result. This experiment result consists with the theoretical analysis, and more different samples will be prepared to verify the universality of the

the universality of the theoretical analysis.

**3.2 Break the lateral diffraction barrier in CARS microscopy** 

**3.2. Break the lateral diffraction barrier in CARS microscopy**

resolution of around 120nm, etc.

theoretical analysis.

44 Photonic Crystals

100%.

*f <sup>S</sup>* (*r*) and *f <sup>L</sup>* (*r*) are the Gaussian Stokes and pump beam distribution which can be replaced by squared cosine function, corresponding to that in equation (34). *f <sup>P</sup>*<sup>1</sup> (*r*) and *f <sup>P</sup>*<sup>2</sup> (*r*) are the first Gaussian and second doughnut probe beams, respectively, which can be similarly re‐ placed by the squared sine and squared cosine functions. Introducing the peak intensities of *I L max* , *IS max* , *IP*<sup>1</sup> *max* and *IP*<sup>2</sup> *max* representing the pump, Stokes, the first and second probe beams respectively, we have four beams' spatial distribution

$$\begin{aligned} f\_L \ \{ r \} &= I\_L^{\max} \cos^2 \left\{ \pi r \mu\_L \sin \alpha \, / \lambda\_L \, \right\} \quad \{ a \} \\ f\_S \ \{ r \} &= I\_S^{\max} \cos^2 \left\{ \pi r \mu\_S \sin \alpha \, / \lambda\_S \right\} \quad \{ b \} \\ f\_{P1} \ \{ r \} &= I\_{P1}^{\max} \sin^2 \left\{ \pi r \mu\_{P1} \sin \alpha \, / \lambda\_{P1} \right\} \, \{ c \} \\ f\_{P2} \ \{ r \} &= I\_{P2}^{\max} \cos^2 \left\{ \pi r \mu\_{P2} \sin \alpha \, / \lambda\_{P2} \right\} \, \{ d \} \end{aligned} \tag{34}$$

where *α* is the semi-aperture of lens and *μL* , *μS* , *μP*1, *μP*2 are the refractive indexes, *λL* , *λS* , *λP*1 , *λP*2 are the wavelength of the four beams, respectively.

The first doughnut probe light (ω*P*1) with intensity distribution *f <sup>P</sup>*<sup>1</sup> (*r*) on the focal plane, almost synchronizes with the pump and Stokes light on the sample. The second Gaussian probe light (ω*P*2) with intensity distribution *f <sup>P</sup>*<sup>2</sup> (*r*) follows the first doughnut subsequently within 1 ps much smaller than the vibrational dephasing time. Consequently, after the course of the first probe light, the effective point-spread-function (PSF) of the useful anti-Stokes signal (ω*A*2) at the suppressed extent of the spot can be expressed as:

$$\mu^\*(r) = f\_{P2}(r) \left[ n\_v - n\_v \bullet n\_P \bullet \frac{1}{A\_{fc}} \partial \sigma\_C \right]$$

Hence, the detailed expression of *h* (*r*) is

$$h(r) = \begin{vmatrix} h^\*(r), & 1 \geq n\_P \bullet \frac{1}{A\_{fc}} \partial \sigma\_C \\ 0, & \text{others} \end{vmatrix} \tag{35}$$

The first term in the square bracket of the equation (35) denotes the total phonon number, and the second term is the phonon number interacting with the first doughnut probe photons at any position of the focal volume. Here A, the cross section area of the volume V, is *V* / *l*. *N* is the molecule number. When the second term in the square bracket is not larger than the first term, i.e., phonon number interacting with the first doughnut probe beam is not larger than the total phonon number generated by the pump and Stokes beams, the expression *h \**(*r*) is tenable, or else, the phonons have been depleted before the second probe beam arrives and *h* (*r*)=0 , therefore, the PSF in the whole space is a piecewise function expressing as:

$$h(r) = \begin{vmatrix} h \ ^\*(r) \mathbf{1} \ge n\_P \bullet \frac{1}{A\_{fc}} \partial \sigma\_{\mathbb{C}} \\ 0 & \text{else} \end{vmatrix} \tag{36}$$

The full width at half maximum (FWHM) of *h \**(*r*), denoted by ∆*r*, is approximated by a Taylor series expansion of *h \**(*r*) to the second order:

placed by the squared sine and squared cosine functions. Introducing the peak intensities of

cos<sup>2</sup> (*πrμL* sin *<sup>α</sup>* / *<sup>λ</sup><sup>L</sup>* ) (*a*)

cos<sup>2</sup> (*πrμS* sin *<sup>α</sup>* / *<sup>λ</sup><sup>S</sup>* ) (*b*)

*sin* <sup>2</sup> (*πrμP*1*sin <sup>α</sup>* / *<sup>λ</sup>P*1) (*c*)

*cos* <sup>2</sup> (*πrμP*2*sin <sup>α</sup>* / *<sup>λ</sup>P*2) (*d*)

where *α* is the semi-aperture of lens and *μL* , *μS* , *μP*1, *μP*2 are the refractive indexes, *λL* , *λS* ,

synchronizes with the pump and Stokes light on the sample. The second Gaussian probe light

much smaller than the vibrational dephasing time. Consequently, after the course of the first probe light, the effective point-spread-function (PSF) of the useful anti-Stokes signal (ω*A*2) at

(*r*) *nv* - *nv* <sup>∙</sup>*nP* <sup>∙</sup> <sup>1</sup>

*Afoc* ∂*σ<sup>C</sup>*

(*r*), <sup>1</sup>≥*nP* <sup>∙</sup> <sup>1</sup>

0, others

*h* (*r*)=0 , therefore, the PSF in the whole space is a piecewise function expressing as:

*Afoc* ∂*σ<sup>C</sup>* 0 *else*

(*r*)1≥*nP* <sup>∙</sup> <sup>1</sup>

The first term in the square bracket of the equation (35) denotes the total phonon number, and the second term is the phonon number interacting with the first doughnut probe photons at any position of the focal volume. Here A, the cross section area of the volume V, is *V* / *l*. *N* is the molecule number. When the second term in the square bracket is not larger than the first term, i.e., phonon number interacting with the first doughnut probe beam is not larger than the total phonon number generated by the pump and Stokes beams, the expression *h \**(*r*) is tenable, or else, the phonons have been depleted before the second probe beam arrives and

*Afoc* ∂*σ<sup>C</sup>*

representing the pump, Stokes, the first and second probe beams

(*r*) follows the first doughnut subsequently within 1 ps

(34)

(35)

(36)

(*r*) on the focal plane, almost

*I L max* , *IS max* , *IP*<sup>1</sup> *max*

46 Photonic Crystals

 and *IP*<sup>2</sup> *max*

(ω*P*2) with intensity distribution *f <sup>P</sup>*<sup>2</sup>

Hence, the detailed expression of *h* (*r*) is

respectively, we have four beams' spatial distribution

*f <sup>P</sup>*<sup>1</sup>

*f <sup>P</sup>*<sup>2</sup>

the suppressed extent of the spot can be expressed as:

*h* \*

*<sup>h</sup>* (*r*)={*<sup>h</sup>* \*

*<sup>h</sup>* (*r*)={*<sup>h</sup>* \*

(*r*)= *f <sup>P</sup>*<sup>2</sup>

*f <sup>L</sup>* (*r*)= *I <sup>L</sup>*

*f <sup>S</sup>* (*r*)= *IS*

(*r*)= *IP*<sup>1</sup> *max*

(*r*)= *IP*<sup>2</sup> *max*

*λP*1 , *λP*2 are the wavelength of the four beams, respectively.

The first doughnut probe light (ω*P*1) with intensity distribution *f <sup>P</sup>*<sup>1</sup>

*max*

*max*

$$
\Delta r = \sqrt{2} \frac{\lambda}{\pi n \sin \alpha \sqrt{3 + \mathcal{K}}} \approx \frac{0.9}{\sqrt{3 + \mathcal{K}}} \bullet \frac{\lambda}{2 \pi \sin \alpha} \tag{37}
$$

Here, influence of the refractive index difference and wavelength difference of the pump, Stokes and probe beams are negligible, and both of them are uniformly expressed as *λ* and *n*. *K* describes magnitude of the doughnut probe beam. Resembling the STED microscopy, the effective PSF, in which the useful anti-Stokes signal (*ωA*2) is measured, becomes de‐ formed by the unwanted anti-Stokes signal (ω*A*1) generation. Resolution is enhanced via a factor of square root of the nonlinear order *K* , namely depletion level, defined by  *K* = *IP*<sup>1</sup> *max* / *IP*<sup>2</sup> *max*. Apparently, if no doughnut probe beam introduced, namely *K* =0, the equa‐ tion (37) largely reproduces the CARS nonlinear optical microscopy's lateral resolution, which is about 3 times better than that of linear optical microscopy. With the gradually in‐ creasing of *K*, the lateral spot width decreases continually, following a square-root law, i.e., resolution is increased by a factor of square root of K. In principle, almost no limitation to super-resolution can be obtained. The resolution enhancement technique can yield an effec‐ tive focal spot which can be fundamentally reduced, theoretically, to an infinitesimal spot, approximating to the size of a molecule or even further to sub-molecular dimensions. Scan‐ ning the suppressed spots automatically renders images with resolution breaking the Abbe's diffraction barrier. Our approach can only fundamentally reduce the focal spot in the lateral direction.

Taking K≈50 to obtain a lateral resolution offor instance, the lateral resolution attainable using this microscopy is approximately ~40 nm, indicating an approach to 5-fold improvement in x or y direction over the diffraction barrier. The PSFs of the traditional and the resolution enhanced CARS microscopy are displayed in figure 11 exploiting the same parameter.

**Figure 11.** Simulation of the point spread function (PSF) of the proposed resolution enhanced CARS microscopy (green, solid line) contrasting with the traditional system (black, dashed line). The new technique indicates that rough‐ ly approximate to 8-fold improvement in lateral direction over the diffraction barrier.

In the APIPD scheme, a super-continuum (SC) source generated from a photonic crystal fiber (PCF), as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. We adopt 500nm and 800nm as the central wavelength of the first and second probe beam, resulting in absolutely distinguishing the two sets of anti-Stokes signals from each other in spectrum. The unwanted anti-Stokes signal on the periphery can be filter out to achieve resolution enhanced images of the molecule specified by its broadband even total CARS spectral signals not only by the anti-Stokes signal of its particular chemical-bond.

#### **3.3. SC used in CARS nanoscopy**

The traditional CARS device consists of a pump and stokes pulses with a frequency *ωL* and *ω<sup>S</sup>* respectively. And the two beams focus onto a sample after being temporally overlapped. The two waves interact together through the third-order susceptibility to generate a blue-shift anti-stokes signal. CARS being a resonant process, the anti-Stokes signal is strongly enhanced when *ω<sup>L</sup>* - *ω<sup>S</sup>* =Ω<sup>r</sup> .To avoid the use of two synchronized laser oscillators emitting the pump and Stokes frequencies required in the process, several previous studies have demonstrated multiplex CARS system by using photonic crystal fibers (PCF) [4]. Super-continuum (SC) source generated from a PCF, as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. The specific and decisive technical advantage in using PCF for CARS measurements is the possibility to generate both pump and Stokes radiations from a single laser. Figure 12 depicts the simple and multiplex CARS process in an energy diagram, the vibrational (or rotational) levels of energy of molecules being schematized by ωR.

In our scheme, a continuum pulse generated from a PCF, as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. We adopt 500nm and 800nm as the central wavelength of the first and second probe beam, resulting in absolutely distinguishing the two sets of anti-Stokes signals from each other in spectrum. The unwanted anti-Stokes signal on the periphery can be filter out to achieve resolution enhanced images of the molecule specified by its broadband even total CARS spectral signals not only by the anti-Stokes signal of its particular chemical-bond.

Good compatibility with practical application requirements and acquisition time in CARS microscopy is made with a continuum pulse serving as the Stokes which is generated by injecting ultrashort pulses into a PCF with high nonlinearity. Because of its broad spectra range, frequency difference between the pump and Stokes can cover most of the biologically inter‐ esting fingerprint region. Therefore, most distinct vibration modes, perhaps full modes even, can be probed simultaneously by a probe pulse as long as different frequency components in the generated SC propagate at about a same speed, known as broadband CARS spectroscopy.

It is critical to generate a very SC pulse to meet its application in CARS microspectroscopy. To acquire the full CARS spectrum of a certain molecule, the frequency difference should at least cover the fingerprint region from 500 cm-1 to 3000 cm-1 at a same moment. This determines that the generated SC pulses have the following characteristics: (i) the spectral bandwidth is broad

Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 49

**Figure 12.** Single CARS process (left) and multiples CARS process with a broadband Stokes pulse

enough, (ii) various components in SC propagate at about a same speed, (iii) spectral intensity should be large and uniform.

#### *3.3.1. SC generation in traditional PCFs*

In the APIPD scheme, a super-continuum (SC) source generated from a photonic crystal fiber (PCF), as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. We adopt 500nm and 800nm as the central wavelength of the first and second probe beam, resulting in absolutely distinguishing the two sets of anti-Stokes signals from each other in spectrum. The unwanted anti-Stokes signal on the periphery can be filter out to achieve resolution enhanced images of the molecule specified by its broadband even total CARS spectral signals not only by the anti-Stokes signal

The traditional CARS device consists of a pump and stokes pulses with a frequency *ωL* and *ω<sup>S</sup>* respectively. And the two beams focus onto a sample after being temporally overlapped. The two waves interact together through the third-order susceptibility to generate a blue-shift anti-stokes signal. CARS being a resonant process, the anti-Stokes signal is strongly enhanced

and Stokes frequencies required in the process, several previous studies have demonstrated multiplex CARS system by using photonic crystal fibers (PCF) [4]. Super-continuum (SC) source generated from a PCF, as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. The specific and decisive technical advantage in using PCF for CARS measurements is the possibility to generate both pump and Stokes radiations from a single laser. Figure 12 depicts the simple and multiplex CARS process in an energy diagram, the vibrational (or rotational) levels of

In our scheme, a continuum pulse generated from a PCF, as the pump and Stokes pluses, is provided with broad spectral width that allows simultaneous detection over a wide range of Raman shifts. We adopt 500nm and 800nm as the central wavelength of the first and second probe beam, resulting in absolutely distinguishing the two sets of anti-Stokes signals from each other in spectrum. The unwanted anti-Stokes signal on the periphery can be filter out to achieve resolution enhanced images of the molecule specified by its broadband even total CARS

Good compatibility with practical application requirements and acquisition time in CARS microscopy is made with a continuum pulse serving as the Stokes which is generated by injecting ultrashort pulses into a PCF with high nonlinearity. Because of its broad spectra range, frequency difference between the pump and Stokes can cover most of the biologically inter‐ esting fingerprint region. Therefore, most distinct vibration modes, perhaps full modes even, can be probed simultaneously by a probe pulse as long as different frequency components in the generated SC propagate at about a same speed, known as broadband CARS spectroscopy.

It is critical to generate a very SC pulse to meet its application in CARS microspectroscopy. To acquire the full CARS spectrum of a certain molecule, the frequency difference should at least cover the fingerprint region from 500 cm-1 to 3000 cm-1 at a same moment. This determines that the generated SC pulses have the following characteristics: (i) the spectral bandwidth is broad

spectral signals not only by the anti-Stokes signal of its particular chemical-bond.

.To avoid the use of two synchronized laser oscillators emitting the pump

of its particular chemical-bond.

**3.3. SC used in CARS nanoscopy**

energy of molecules being schematized by ωR.

when *ω<sup>L</sup>* - *ω<sup>S</sup>* =Ω<sup>r</sup>

48 Photonic Crystals

Although dispersive property of a PCF leads to different propagating velocities of different frequency components in the generated SC, soliton pulse trains from a PCF have been considered as the Stokes in many CARS spectroscopy schemes [75-77]. In this part, we will detailed discuss the generation and characteristics of the soliton pulse trains used in CARS spectroscopy.

When an ultrashort pulse travels in a nonlinear material such as a PCF, different frequency components in it will propagate at different velocity, because of the dispersive effect, including material dispersion and waveguide dispersion. In the anomalous regime, this dispersive effect manifests blue-shift components moving faster than the red-shift components, leading to the temporal broadening. At the same time, SPM generates a nonlinear chirp which is negative near the leading edge and positive near the trailing edge. Therefore the propagating speed decreases near the leading edge and increases near the trailing edge, making the pulse narrow down. Under certain conditions, the optical envelope will keep its shape and propagating velocity known as optical soliton [78,79].

As discussed in section two, distinct peaks in time domain are obvious. Because the soliton ejected earlier has larger group velocity, higher peak power and shorter durations, it will propagate faster than others. So the separation between these peaks becomes larger with a longer propagation distance. Figure 13 shows the spectrograms of the generated SC when the propagation is 10 cm, 15 cm, 20 cm and 25 cm. The input pulse is central at 800 nm with peak power 10kW and pulse duration 50fs as before. From the figures, the first-order soliton shifts from 2 ps to about 8.5 ps when the propagation changes from 10 cm to 25 cm gradually. Also

power, but quite few.

the time delay between solitions becomes larger with a longer propagation distance by making contrast of the four graphs. For the dispersion wave on the short wavelength side, time delay between different frequency components also becomes larger. The time span is about 2 ps when the propagation distance is 10 cm, while it is about 5 ps with 25 cm travelling distance. The shifting of solitons in time domain is corresponding to the spectral shifting. The central wavelength of the first-order soliton shifts from 1100 nm to 1200 nm with the propagation from 10 cm to 25 cm. when the propagation changes from 10 cm to 25 cm gradually. Also the time delay between solitions becomes larger with a longer propagation distance by making contrast of the four graphs. For the dispersion wave on the short wavelength side, time delay between different frequency components also becomes larger. The time span is about 2 ps when the propagation distance is 10 cm, while it is about 5 ps with 25 cm travelling distance. The shifting of solitons in time domain is corresponding to the spectral shifting. The central wavelength of the first-order soliton shifts from

1100 nm to 1200 nm with the propagation from 10 cm to 25 cm.

duration 50fs as before. From the figures, the first-order soliton shifts from 2 ps to about 8.5 ps

Fig.13 Spectrograms at different propagation distances for input pulse at 800 nm with peak power 10kW and pulse duration 50fs (a) 10 cm; (b) 15 cm; (c) 20 cm; (d) 25 cm; **Figure 13.** Spectrograms at different propagation distances for input pulse at 800 nm with peak power 10kW and pulse duration 50fs (a) 10 cm; (b) 15 cm; (c) 20 cm; (d) 25 cm;

Although soliton trains will shift towards long-wavelength side due to SSFS[80,81], the pulse duration and spectral bandwidth of first-order soliton are nearly invariable. Figure 14 shows the temporal and spectal profiles of the first-order soliton with different propagating distance with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm. The first-order soliton is selected as both the pump and Stokes just for it has shorter duration and wider spectral bandwidth. It is estimated that the soliton has a pulse duration about 30 fs and spectral bandwidth about 47 nm with 10 cm Although soliton trains will shift towards long-wavelength side due to SSFS [80,81], the pulse duration and spectral bandwidth of first-order soliton are nearly invariable. Figure 14 shows the temporal and spectal profiles of the first-order soliton with different propagating distance with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm. The first-order soliton is selected as both the pump and Stokes just for it has shorter duration and wider spectral bandwidth. It is estimated that the soliton has a pulse duration about 30 fs and spectral bandwidth about 47 nm with 10 cm propagation distance. With further propagation, the pulse duration will increase from 30 to 42 fs, accompanied with the spectral bandwidth decreasing from 47 to 41

propagation distance. With further propagation, the pulse duration will increase from 30 to 42 fs, accompanied with the spectral bandwidth decreasing from 47 to 41 nm induced by high-order dispersive and nonlinear effects. Of course the bandwidth is increasing with improving the peak

nm induced by high-order dispersive and nonlinear effects. Of course the bandwidth is increasing with improving the peak power, but quite few.

the time delay between solitions becomes larger with a longer propagation distance by making contrast of the four graphs. For the dispersion wave on the short wavelength side, time delay between different frequency components also becomes larger. The time span is about 2 ps when the propagation distance is 10 cm, while it is about 5 ps with 25 cm travelling distance. The shifting of solitons in time domain is corresponding to the spectral shifting. The central wavelength of the first-order soliton shifts from 1100 nm to 1200 nm with the propagation from

15 cm, 20 cm and 25 cm. The input pulse is central at 800 nm with peak power 10kW and pulse duration 50fs as before. From the figures, the first-order soliton shifts from 2 ps to about 8.5 ps when the propagation changes from 10 cm to 25 cm gradually. Also the time delay between solitions becomes larger with a longer propagation distance by making contrast of the four graphs. For the dispersion wave on the short wavelength side, time delay between different frequency components also becomes larger. The time span is about 2 ps when the propagation distance is 10 cm, while it is about 5 ps with 25 cm travelling distance. The shifting of solitons in time domain is corresponding to the spectral shifting. The central wavelength of the first-order soliton shifts from

Fig.13 Spectrograms at different propagation distances for input pulse at 800 nm with peak power 10kW and pulse duration 50fs (a) 10 cm; (b) 15 cm; (c) 20 cm; (d) 25 cm;

duration and spectral bandwidth of first-order soliton are nearly invariable. Figure 14 shows the temporal and spectal profiles of the first-order soliton with different propagating distance with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm. The first-order soliton is selected as both the pump and Stokes just for it has shorter duration and wider spectral bandwidth. It is estimated that the soliton has a pulse duration about 30 fs and spectral bandwidth about 47 nm with 10 cm propagation distance. With further propagation, the pulse duration will increase from 30 to 42 fs, accompanied with the spectral bandwidth decreasing from 47 to 41 nm induced by high-order dispersive and nonlinear effects. Of course the bandwidth is increasing with improving the peak

Although soliton trains will shift towards long-wavelength side due to SSFS[80,81], the pulse

Although soliton trains will shift towards long-wavelength side due to SSFS [80,81], the pulse duration and spectral bandwidth of first-order soliton are nearly invariable. Figure 14 shows the temporal and spectal profiles of the first-order soliton with different propagating distance with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm. The first-order soliton is selected as both the pump and Stokes just for it has shorter duration and wider spectral bandwidth. It is estimated that the soliton has a pulse duration about 30 fs and spectral bandwidth about 47 nm with 10 cm propagation distance. With further propagation, the pulse duration will increase from 30 to 42 fs, accompanied with the spectral bandwidth decreasing from 47 to 41

**Figure 13.** Spectrograms at different propagation distances for input pulse at 800 nm with peak power 10kW and pulse

10 cm to 25 cm.

50 Photonic Crystals

power, but quite few.

1100 nm to 1200 nm with the propagation from 10 cm to 25 cm.

(a) (b)

(c) (d)

duration 50fs (a) 10 cm; (b) 15 cm; (c) 20 cm; (d) 25 cm;

Fig.14 (a ~ c) temporal and (d ~ f) spectral structure of the first-order soliton with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm **Figure 14.** (a ~ c) temporal and (d ~ f) spectral structure of the first-order soliton with (a) (d) at 10 cm, (b) (e) at 15 cm and (c) (f) at 25 cm

Though this first-order soliton shows good temporal and spectral structures, the limited spectral bandwidth makes it not the best choice for CARS microspectroscopy. That is because the

Recent studies show that SC with broad spectral bandwidth, uniform intensity profile and

good coherence is obtained by launching ultrashort laser pulses into an all-normal-dispersion PCF[82,83]. These PCFs exhibit a convex dispersive profile which is flat near the maximum dispersion wavelength (MDW). And they possess none ZDW point in the visible and near-infrared spectral region, known as all-normal dispersion. Without anomalous dispersion region, all the soliton-related dynamics are suppressed, including breakup of input pulse into several pulse trains. And the generated SC is not so sensitive to fluctuations of input pulse and shot noise of pump

pump and Stokes (soliton pulse) can resonant only few molecular bonds instead of the whole bonds, leading to a CARS spectrum with few Raman distinct peaks. To acquire all the Raman distinct peaks simultaneously, SC with both good temporal structure and broad spectral bandwidth is

needed.

3.2.2 SC with all-normal dispersion PCF

Though this first-order soliton shows good temporal and spectral structures, the limited spectral bandwidth makes it not the best choice for CARS microspectroscopy. That is because the pump and Stokes (soliton pulse) can resonant only few molecular bonds instead of the whole bonds, leading to a CARS spectrum with few Raman distinct peaks. To acquire all the Raman distinct peaks simultaneously, SC with both good temporal structure and broad spectral bandwidth is needed.

#### *3.3.2. SC with all-normal dispersion PCF*

which point the nonlinear coefficient is about 37 W��m��.

single pulse all the time in propagation.

near pulse edges are related to optical wave breaking[85].

Recent studies show that SC with broad spectral bandwidth, uniform intensity profile and good coherence is obtained by launching ultrashort laser pulses into an all-normal-dispersion PCF [82,83]. These PCFs exhibit a convex dispersive profile which is flat near the maximum dispersion wavelength (MDW). And they possess none ZDW point in the visible and nearinfrared spectral region, known as all-normal dispersion. Without anomalous dispersion region, all the soliton-related dynamics are suppressed, including breakup of input pulse into several pulse trains. And the generated SC is not so sensitive to fluctuations of input pulse and shot noise of pump laser [84]. Because SC generation in these PCFs is dominated by the SPM effect and optical wave breaking, the SC pulse can keep characteristic of a single pulse during the propagation. And the SC pulse is easy for dispersion compensation with such smooth phase distribution. These properties make the generated SC in all-normal dispersion PCFs especially suit for time-resolved applications, such as CARS spectroscopy. laser[84]. Because SC generation in these PCFs is dominated by the SPM effect and optical wave breaking, the SC pulse can keep characteristic of a single pulse during the propagation. And the SC pulse is easy for dispersion compensation with such smooth phase distribution. These properties make the generated SC in all-normal dispersion PCFs especially suit for time-resolved applications, such as CARS spectroscopy.

In the following simulations, an all-normal dispersion PCF (from NKT Photonics) with hole diameter 0.5616 µm and hole pitch 1.44 µm is used. Figure 15 shows the cross-section and dispersion profile of the all-normal dispersion PCF respectively. From figure (b), it is clear that this fiber shows all-normal dispersion in the range from visible to near-infrared. And the dispersion parameter D is estimated about -10 ps / nm∙km at wavelength about 1060 nm, at which point the nonlinear coefficient is about 37 W-1 m-1 . In the following simulations, an all-normal dispersion PCF (from NKT Photonics) with hole diameter 0.5616 μm and hole pitch 1.44 μm is used. Figure 15 shows the cross-section and dispersion profile of the all-normal dispersion PCF respectively. From figure (b), it is clear that this fiber shows all-normal dispersion in the range from visible to near-infrared. And the dispersion parameter D is estimated about �10 ����m � �m at wavelength about 1060 nm, at

**Figure 15.** (a) Cross-section model used in simulation and (b) dispersion profile of the all-normal dispersion PCF

Fig.15 (a) Cross-section model used in simulation and (b) dispersion profile of the all-normal dispersion PCF

hyperbolic secant pulses centered at 1060 nm with pulse duration 100 fs and peak power 30 kw are injected into the PCF. From figure 16 (b), the spectral broadening is first concluded at the short wavelength through propagation about 38 mm, then at the long wavelength side about 56mm. The spectrum will gain no more broadening with further propagation, except that the generated SC has a smoother spectral profile. In temporal domain, the pulse duration changes a little in the initial stage, but increases gradually with further propagation. And the input pulse keeps its character of

Temporal and spectral evolutions with propagating distance are shown in figure 16 when

To focus more detailed structures in temporal and spectral domain, evolution slices at

different propagating distances are plotted in figure 17 with linear scales. The most notable feature in time domain lies in the conservation of a single pulse during the whole propagation. Although the pulse is broadening during the propagation, self-steepening effect leads to a steeper trailing edge which is clear in slices within 3 centimeters' propagation. After that, the pulse edges becomes steeper which makes the pulse shape like a rectangular function, because even weak dispersion effect also brings about significant changes of pulse shape under abundant SPM-induced frequency chirp. Oscillating structure near the trailing edge as shown in the left graph with 3cm propagation derives from pulse distortion induced by high-order dispersion. And the oscillations Temporal and spectral evolutions with propagating distance are shown in figure 16 when hyperbolic secant pulses centered at 1060 nm with pulse duration 100 fs and peak power 30 kw are injected into the PCF. From figure 16 (b), the spectral broadening is first concluded at the short wavelength through propagation about 38 mm, then at the long wavelength side about 56mm. The spectrum will gain no more broadening with further propagation, except that the generated SC has a smoother spectral profile. In temporal domain, the pulse duration changes a little in the initial stage, but increases gradually with further propagation. And the input pulse keeps its character of single pulse all the time in propagation.

Though this first-order soliton shows good temporal and spectral structures, the limited spectral bandwidth makes it not the best choice for CARS microspectroscopy. That is because the pump and Stokes (soliton pulse) can resonant only few molecular bonds instead of the whole bonds, leading to a CARS spectrum with few Raman distinct peaks. To acquire all the Raman distinct peaks simultaneously, SC with both good temporal structure and broad

Recent studies show that SC with broad spectral bandwidth, uniform intensity profile and good coherence is obtained by launching ultrashort laser pulses into an all-normal-dispersion PCF [82,83]. These PCFs exhibit a convex dispersive profile which is flat near the maximum dispersion wavelength (MDW). And they possess none ZDW point in the visible and nearinfrared spectral region, known as all-normal dispersion. Without anomalous dispersion region, all the soliton-related dynamics are suppressed, including breakup of input pulse into several pulse trains. And the generated SC is not so sensitive to fluctuations of input pulse and shot noise of pump laser [84]. Because SC generation in these PCFs is dominated by the SPM effect and optical wave breaking, the SC pulse can keep characteristic of a single pulse during the propagation. And the SC pulse is easy for dispersion compensation with such smooth phase distribution. These properties make the generated SC in all-normal dispersion PCFs especially

laser[84]. Because SC generation in these PCFs is dominated by the SPM effect and optical wave breaking, the SC pulse can keep characteristic of a single pulse during the propagation. And the SC pulse is easy for dispersion compensation with such smooth phase distribution. These properties make the generated SC in all-normal dispersion PCFs especially suit for time-resolved

In the following simulations, an all-normal dispersion PCF (from NKT Photonics) with hole diameter 0.5616 µm and hole pitch 1.44 µm is used. Figure 15 shows the cross-section and dispersion profile of the all-normal dispersion PCF respectively. From figure (b), it is clear that this fiber shows all-normal dispersion in the range from visible to near-infrared. And the dispersion parameter D is estimated about -10 ps / nm∙km at wavelength about 1060 nm, at

diameter 0.5616 μm and hole pitch 1.44 μm is used. Figure 15 shows the cross-section and dispersion profile of the all-normal dispersion PCF respectively. From figure (b), it is clear that this fiber shows all-normal dispersion in the range from visible to near-infrared. And the dispersion parameter D is estimated about �10 ����m � �m at wavelength about 1060 nm, at

Fig.15 (a) Cross-section model used in simulation and (b) dispersion profile of the all-normal dispersion PCF

**Figure 15.** (a) Cross-section model used in simulation and (b) dispersion profile of the all-normal dispersion PCF

hyperbolic secant pulses centered at 1060 nm with pulse duration 100 fs and peak power 30 kw are injected into the PCF. From figure 16 (b), the spectral broadening is first concluded at the short wavelength through propagation about 38 mm, then at the long wavelength side about 56mm. The spectrum will gain no more broadening with further propagation, except that the generated SC has a smoother spectral profile. In temporal domain, the pulse duration changes a little in the initial stage, but increases gradually with further propagation. And the input pulse keeps its character of

Temporal and spectral evolutions with propagating distance are shown in figure 16 when

To focus more detailed structures in temporal and spectral domain, evolution slices at

different propagating distances are plotted in figure 17 with linear scales. The most notable feature in time domain lies in the conservation of a single pulse during the whole propagation. Although the pulse is broadening during the propagation, self-steepening effect leads to a steeper trailing edge which is clear in slices within 3 centimeters' propagation. After that, the pulse edges becomes steeper which makes the pulse shape like a rectangular function, because even weak dispersion effect also brings about significant changes of pulse shape under abundant SPM-induced frequency chirp. Oscillating structure near the trailing edge as shown in the left graph with 3cm propagation derives from pulse distortion induced by high-order dispersion. And the oscillations

In the following simulations, an all-normal dispersion PCF (from NKT Photonics) with hole

m-1 .

spectral bandwidth is needed.

52 Photonic Crystals

applications, such as CARS spectroscopy.

single pulse all the time in propagation.

near pulse edges are related to optical wave breaking[85].

*3.3.2. SC with all-normal dispersion PCF*

suit for time-resolved applications, such as CARS spectroscopy.

which point the nonlinear coefficient is about 37 W-1

(a) (b)

which point the nonlinear coefficient is about 37 W��m��.

To focus more detailed structures in temporal and spectral domain, evolution slices at different propagating distances are plotted in figure 17 with linear scales. The most notable feature in time domain lies in the conservation of a single pulse during the whole propagation. Although the pulse is broadening during the propagation, self-steepening effect leads to a steeper trailing edge which is clear in slices within 3 centimeters' propagation. After that, the pulse edges becomes steeper which makes the pulse shape like a rectangular function, because even weak dispersion effect also brings about significant changes of pulse shape under abundant SPMinduced frequency chirp. Oscillating structure near the trailing edge as shown in the left graph with 3cm propagation derives from pulse distortion induced by high-order dispersion. And the oscillations near pulse edges are related to optical wave breaking [85].

Fig.16 Temporal (left) and spectral (right) evolutions over propagation distance for input pulse centered at 1060 nm with peak power 30 kW and pulse duration 100 fs **Figure 16.** Temporal (left) and spectral (right) evolutions over propagation distance for input pulse centered at 1060 nm with peak power 30 kW and pulse duration 100 fs

In spectral domain, the spectral broadening induced by SPM effect is significant with a clear oscillating structure through the full wavelength. The SPM-induced spectral broadening spreads out towards two sides of the pump wavelength with red-shift components near the leading edge and blue-shift components near the trailing edge. And in the later propagation, the oscillating structure is explained by optical wave breaking. The red-shift components near the leading edge transfer faster and catch up with the red components near leading edge, leading to the interference between them. The blue-shift components near the trailing edge are just an opposite to that. This In spectral domain, the spectral broadening induced by SPM effect is significant with a clear oscillating structure through the full wavelength. The SPM-induced spectral broadening spreads out towards two sides of the pump wavelength with red-shift components near the leading edge and blue-shift components near the trailing edge. And in the later propagation, the oscillating structure is explained by optical wave breaking. The red-shift components near the leading edge transfer faster and catch up with the red components near leading edge, leading to the interference between them. The blue-shift components near the trailing edge are just an opposite to that. This oscillating structure in spectral domain and the oscillations near pulse edge in the temporal domain are manifestations of the same phenomenon.

oscillating structure in spectral domain and the oscillations near pulse edges in the temporal

Fig.17 Evolution slices at different propagation distances in temporal (left) and spectral (right) domain

structures, it is useful to study the impact of input pulses on the SC generation, such as pulse

of pulse energy on SC generating. The input pulses are centered at 1060 nm with pulse duration 100 fs and propagation distance 15 cm. Figure 18(a) shows the variations of spectral bandwidth

Although SC generated in all-normal dispersion PCF has good temporal and spectral

First, the central wavelength, pulse duration and fiber length are restricted to show the impact

domain are manifestations of the same phenomenon.

(a) (b)

duration, pulse energy and so on.

domain are manifestations of the same phenomenon.

(a) (b)

Fig.16 Temporal (left) and spectral (right) evolutions over propagation distance for input pulse centered at 1060 nm with peak power 30 kW and pulse duration 100 fs

oscillating structure through the full wavelength. The SPM-induced spectral broadening spreads out towards two sides of the pump wavelength with red-shift components near the leading edge and blue-shift components near the trailing edge. And in the later propagation, the oscillating structure is explained by optical wave breaking. The red-shift components near the leading edge transfer faster and catch up with the red components near leading edge, leading to the interference between them. The blue-shift components near the trailing edge are just an opposite to that. This

In spectral domain, the spectral broadening induced by SPM effect is significant with a clear

Although SC generated in all-normal dispersion PCF has good temporal and spectral **Figure 17.** Evolution slices at different propagation distances in temporal (left) and spectral (right) domain

structures, it is useful to study the impact of input pulses on the SC generation, such as pulse duration, pulse energy and so on. First, the central wavelength, pulse duration and fiber length are restricted to show the impact Although SC generated in all-normal dispersion PCF has good temporal and spectral struc‐ tures, it is useful to study the impact of input pulses on the SC generation, such as pulse duration, pulse energy and so on.

Fig.17 Evolution slices at different propagation distances in temporal (left) and spectral (right) domain

of pulse energy on SC generating. The input pulses are centered at 1060 nm with pulse duration 100 fs and propagation distance 15 cm. Figure 18(a) shows the variations of spectral bandwidth First, the central wavelength, pulse duration and fiber length are restricted to show the impact of pulse energy on SC generating. The input pulses are centered at 1060 nm with pulse duration 100 fs and propagation distance 15 cm. Figure 18(a) shows the variations of spectral bandwidth with different pulse energies ranging from 2 nJ to 10 nJ. It is clear that the generated SC shows larger spectral broadening, much more than 600 nm with pulse energy 4 nJ. For a higher pulse energy, the SC generating is much influenced by the SPM and self-steepening effects, leading to the distinct peaks on each side of the pump wavelength. But the intensity near the pump dips more due to them. The spectral bandwidth is increasing, but the increment becomes small, with pulse energy.

Figure 18(b) shows the impact of pulse duration on generated SC with peak power being a constant 30 kW. The relation between peak power and pulse duration is E=P∙ ∆τ with E being pulse energy, P and ∆τ being peak power and pulse duration at half maximum. It is clear that the spectral broadening changes few with different pulse durations. The spectrum reaches 600 nm when the peak power keeps 30 kW. With larger pulse duration, the oscillating structures

due to them. The spectral bandwidth is increasing, but the increment becomes small, with pulse Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 55

with different pulse energies ranging from 2 nJ to 10 nJ. It is clear that the generated SC shows larger spectral broadening, much more than 600 nm with pulse energy 4 nJ. For a higher pulse energy, the SC generating is much influenced by the SPM and self-steepening effects, leading to the distinct peaks on each side of the pump wavelength. But the intensity near the pump dips more

Fig.18 Spectral broadening with different pulse (a) energies and (b) durations **Figure 18.** Spectral broadening with different pulse (a) energies and (b) durations

additional difficulty to dispersion compensation.

energy.

Fig.16 Temporal (left) and spectral (right) evolutions over propagation distance for input pulse centered at 1060 nm with peak power 30 kW and pulse duration 100 fs

oscillating structure through the full wavelength. The SPM-induced spectral broadening spreads out towards two sides of the pump wavelength with red-shift components near the leading edge and blue-shift components near the trailing edge. And in the later propagation, the oscillating structure is explained by optical wave breaking. The red-shift components near the leading edge transfer faster and catch up with the red components near leading edge, leading to the interference between them. The blue-shift components near the trailing edge are just an opposite to that. This oscillating structure in spectral domain and the oscillations near pulse edges in the temporal

Fig.17 Evolution slices at different propagation distances in temporal (left) and spectral (right) domain

**Figure 17.** Evolution slices at different propagation distances in temporal (left) and spectral (right) domain

structures, it is useful to study the impact of input pulses on the SC generation, such as pulse

Although SC generated in all-normal dispersion PCF has good temporal and spectral struc‐ tures, it is useful to study the impact of input pulses on the SC generation, such as pulse

of pulse energy on SC generating. The input pulses are centered at 1060 nm with pulse duration 100 fs and propagation distance 15 cm. Figure 18(a) shows the variations of spectral bandwidth

First, the central wavelength, pulse duration and fiber length are restricted to show the impact of pulse energy on SC generating. The input pulses are centered at 1060 nm with pulse duration 100 fs and propagation distance 15 cm. Figure 18(a) shows the variations of spectral bandwidth with different pulse energies ranging from 2 nJ to 10 nJ. It is clear that the generated SC shows larger spectral broadening, much more than 600 nm with pulse energy 4 nJ. For a higher pulse energy, the SC generating is much influenced by the SPM and self-steepening effects, leading to the distinct peaks on each side of the pump wavelength. But the intensity near the pump dips more due to them. The spectral bandwidth is increasing, but the increment becomes small,

Figure 18(b) shows the impact of pulse duration on generated SC with peak power being a constant 30 kW. The relation between peak power and pulse duration is E=P∙ ∆τ with E being pulse energy, P and ∆τ being peak power and pulse duration at half maximum. It is clear that the spectral broadening changes few with different pulse durations. The spectrum reaches 600 nm when the peak power keeps 30 kW. With larger pulse duration, the oscillating structures

Although SC generated in all-normal dispersion PCF has good temporal and spectral

First, the central wavelength, pulse duration and fiber length are restricted to show the impact

domain are manifestations of the same phenomenon.

54 Photonic Crystals

(a) (b)

(a) (b)

duration, pulse energy and so on.

with pulse energy.

duration, pulse energy and so on.

In spectral domain, the spectral broadening induced by SPM effect is significant with a clear

in output spectrum are more obvious. Also the dips near pump wavelength are much more clear. Due to optical wave breaking induced four-wave mixing, the sidebands at each short and long wavelength side are significant, and their peaks are moving towards short and long wavelength direction respectively. Figure 18(b) shows the impact of pulse duration on generated SC with peak power being a constant 30 kW. The relation between peak power and pulse duration is � � � � ∆τ with E being pulse energy, P and ∆τ being peak power and pulse duration at half maximum. It is clear that the

Although the generated SC has a larger spectral broadening when the PCF is pumped near the MDW, the impacts of central wavelength of input pulses on SC generation are analyzed. In figure 19, temporal and spectral structures of generated SC with different pump wave‐ lengths are plotted where the pulse duration, pulse energy and propagation distance are fixed to 100 fs, 2 nJ and 15 cm respectively. For 1060 nm pump, the spectrum is approximately symmetric from 800 nm to 1350 nm with clear oscillating structures near the pump wavelength. For pump wavelength far from the MDW, not only the full spectrum is shorter that that with 1060 nm pump, but the spectrum shows a tail towards the MDW with much lower amplitude than the peaks induced by SPM. There are obviously small tails near the leading edge for 1200 nm pump and trailing edge for 900 nm pump in figure (a) which are induced by high-order dispersion. These tails will apply additional difficulty to dispersion compensation. spectral broadening changes few with different pulse durations. The spectrum reaches 600 nm when the peak power keeps 30 kW. With larger pulse duration, the oscillating structures in output spectrum are more obvious. Also the dips near pump wavelength are much more clear. Due to optical wave breaking induced four-wave mixing, the sidebands at each short and long wavelength side are significant, and their peaks are moving towards short and long wavelength direction respectively. Although the generated SC has a larger spectral broadening when the PCF is pumped near the MDW, the impacts of central wavelength of input pulses on SC generation are analyzed. In figure 19, temporal and spectral structures of generated SC with different pump wavelengths are

From above analyses and discussions, the all-normal dispersion PCF should be pump near the maximum-dispersion point with proper peak power and pulse duration to generate an excellent SC suitable for CARS spectroscopy with broad spectral bandwidth, uniform spectral plotted where the pulse duration, pulse energy and propagation distance are fixed to 100 fs, 2 nJ and 15 cm respectively. For 1060 nm pump, the spectrum is approximately symmetric from 800 nm to 1350 nm with clear oscillating structures near the pump wavelength. For pump wavelength

far from the MDW, not only the full spectrum is shorter that that with 1060 nm pump, but the spectrum shows a tail towards the MDW with much lower amplitude than the peaks induced by SPM. There are obviously small tails near the leading edge for 1200 nm pump and trailing edge for 900 nm pump in figure (a) which are induced by high-order dispersion. These tails will apply

bad.

 Fig.19 Pulse shapes with different pump wavelengths in (a) temporal and (b) spectral domain **Figure 19.** Pulse shapes with different pump wavelengths in (a) temporal and (b) spectral domain

intensities and smooth and excellent temporal structure with smooth phase distribution. The spectral intensity profile is much more smooth with shorter pulse duration under the same peak power. Although higher peak power is beneficial to spectral broadening, it often causes strong tails near the pulse edges and bad flatness of spectral intensity profile. And the fiber length is also very important, because a longer fiber length makes a more flat spectral profile while temporal structure becomes bad. From above analyses and discussions, the all-normal dispersion PCF should be pump near the maximum-dispersion point with proper peak power and pulse duration to generate an excellent SC suitable for CARS spectroscopy with broad spectral bandwidth, uniform spectral intensities and smooth and excellent temporal structure with smooth phase distribution. The spectral intensity profile is much more smooth with shorter pulse duration under the same peak power. Although

Based on the basic experimental facilities, pulses centered at 1060 nm with pulse duration 50 fs and pulse energy 2nJ are injected into the all-normal dispersion PCF. The spectrum of output pulse can be expanded to 600~ nm (from ~800 nm to 1400~ nm) after 12 cm propagation with the spectrogram of generated SC being displayed in figure 20(a). This SC shows good spectral continuity and uniformity with tails on both the long and short wavelength sides caused by high-order dispersion. Figure 20(b) demonstrates the temporal distribution of the generated SC pulse with slightly asymmetric structure towards the tailing edge. And the tails on long and short wavelength sides are related to the pedestals near leading and trailing edges. From figure (b), the delay time between different components is less than 2 ps and the temporal phase is predominantly quadratic, which makes the dispersion compensation much easier. higher peak power is beneficial to spectral broadening, it often causes strong tails near the pulse edges and bad flatness of spectral intensity profile. And the fiber length is also very important, because a longer fiber length makes a more flat spectral profile while temporal structure becomes Based on the basic experimental facilities, pulses centered at 1060 nm with pulse duration 50 fs and pulse energy 2nJ are injected into the all-normal dispersion PCF. The spectrum of output pulse can be expanded to 600~ nm (from ~800 nm to 1400~ nm) after 12 cm propagation with the spectrogram of generated SC being displayed in figure 20(a). This SC shows good spectral

To exploit this continuum pulse into CARS spectroscopy, it is necessary to compensate the delay time induced by dispersive effect and other nonlinear effects, leading to a same propa‐ gating speed of different frequency components. In terms of such small dispersion, a pair of prisms made of BK7 glass are used to compensate it. The spectrogram and temporal, spectral profiles after dispersion compensation are plotted in figure 21, with (a) being the spectrogram, (b) the temporal profile and (c) the spectral profile. After dispersion compensation, most components in the generated SC are propagating at a same speed, except the tails on long and short wavelength sides, for the prisms are not effective enough to compensation of high-order dispersion. From figure (b), the pulse duration of generated SC is compressed to about 15 fs from initial ~2 ps in spite of a little oscillation near the pulse edges. Figure (c) shows the spectral profile on a linear scale with bandwidth over 600 nm. And the spectral intensity possesses a good uniformity. continuity and uniformity with tails on both the long and short wavelength sides caused by high-order dispersion. Figure 20(b) demonstrates the temporal distribution of the generated SC pulse with slightly asymmetric structure towards the tailing edge. And the tails on long and short wavelength sides are related to the pedestals near leading and trailing edges. From figure (b), the delay time between different components is less than 2 ps and the temporal phase is predominantly quadratic, which makes the dispersion compensation much easier. To exploit this continuum pulse into CARS spectroscopy, it is necessary to compensate the delay time induced by dispersive effect and other nonlinear effects, leading to a same propagating speed of different frequency components. In terms of such small dispersion, a pair of prisms made

of BK7 glass are used to compensate it. The spectrogram and temporal, spectral profiles after dispersion compensation are plotted in figure 21, with (a) being the spectrogram, (b) the temporal profile and (c) the spectral profile. After dispersion compensation, most components in the generated SC are propagating at a same speed, except the tails on long and short wavelength sides, for the prisms are not effective enough to compensation of high-order dispersion. From figure (b), the pulse duration of generated SC is compressed to about 15 fs from initial ~2 ps in spite of a little oscillation near the pulse edges. Figure (c) shows the spectral profile on a linear scale with

bandwidth over 600 nm. And the spectral intensity possesses a good uniformity.

Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 57

**Figure 20.** (a) Spectrogram and (b) temporal profile of the generated SC pulse

(a) (b) (c)

intensities and smooth and excellent temporal structure with smooth phase distribution. The spectral intensity profile is much more smooth with shorter pulse duration under the same peak power. Although higher peak power is beneficial to spectral broadening, it often causes strong tails near the pulse edges and bad flatness of spectral intensity profile. And the fiber length is also very important, because a longer fiber length makes a more flat spectral profile

maximum-dispersion point with proper peak power and pulse duration to generate an excellent SC suitable for CARS spectroscopy with broad spectral bandwidth, uniform spectral intensities and smooth and excellent temporal structure with smooth phase distribution. The spectral intensity profile is much more smooth with shorter pulse duration under the same peak power. Although higher peak power is beneficial to spectral broadening, it often causes strong tails near the pulse edges and bad flatness of spectral intensity profile. And the fiber length is also very important, because a longer fiber length makes a more flat spectral profile while temporal structure becomes

From above analyses and discussions, the all-normal dispersion PCF should be pump near the

Fig.19 Pulse shapes with different pump wavelengths in (a) temporal and (b) spectral domain

**Figure 19.** Pulse shapes with different pump wavelengths in (a) temporal and (b) spectral domain

Based on the basic experimental facilities, pulses centered at 1060 nm with pulse duration 50 fs and pulse energy 2nJ are injected into the all-normal dispersion PCF. The spectrum of output pulse can be expanded to 600~ nm (from ~800 nm to 1400~ nm) after 12 cm propagation with the spectrogram of generated SC being displayed in figure 20(a). This SC shows good spectral continuity and uniformity with tails on both the long and short wavelength sides caused by high-order dispersion. Figure 20(b) demonstrates the temporal distribution of the generated SC pulse with slightly asymmetric structure towards the tailing edge. And the tails on long and short wavelength sides are related to the pedestals near leading and trailing edges. From figure (b), the delay time between different components is less than 2 ps and the temporal phase is predominantly quadratic, which makes the dispersion compensation much easier.

Based on the basic experimental facilities, pulses centered at 1060 nm with pulse duration

To exploit this continuum pulse into CARS spectroscopy, it is necessary to compensate the delay time induced by dispersive effect and other nonlinear effects, leading to a same propa‐ gating speed of different frequency components. In terms of such small dispersion, a pair of prisms made of BK7 glass are used to compensate it. The spectrogram and temporal, spectral profiles after dispersion compensation are plotted in figure 21, with (a) being the spectrogram, (b) the temporal profile and (c) the spectral profile. After dispersion compensation, most components in the generated SC are propagating at a same speed, except the tails on long and short wavelength sides, for the prisms are not effective enough to compensation of high-order dispersion. From figure (b), the pulse duration of generated SC is compressed to about 15 fs from initial ~2 ps in spite of a little oscillation near the pulse edges. Figure (c) shows the spectral profile on a linear scale with bandwidth over 600 nm. And the spectral intensity possesses a

To exploit this continuum pulse into CARS spectroscopy, it is necessary to compensate the

delay time induced by dispersive effect and other nonlinear effects, leading to a same propagating speed of different frequency components. In terms of such small dispersion, a pair of prisms made of BK7 glass are used to compensate it. The spectrogram and temporal, spectral profiles after dispersion compensation are plotted in figure 21, with (a) being the spectrogram, (b) the temporal profile and (c) the spectral profile. After dispersion compensation, most components in the generated SC are propagating at a same speed, except the tails on long and short wavelength sides, for the prisms are not effective enough to compensation of high-order dispersion. From figure (b), the pulse duration of generated SC is compressed to about 15 fs from initial ~2 ps in spite of a little oscillation near the pulse edges. Figure (c) shows the spectral profile on a linear scale with

predominantly quadratic, which makes the dispersion compensation much easier.

bandwidth over 600 nm. And the spectral intensity possesses a good uniformity.

50 fs and pulse energy 2nJ are injected into the all-normal dispersion PCF. The spectrum of output pulse can be expanded to 600~ nm (from ~800 nm to 1400~ nm) after 12 cm propagation with the spectrogram of generated SC being displayed in figure 20(a). This SC shows good spectral continuity and uniformity with tails on both the long and short wavelength sides caused by high-order dispersion. Figure 20(b) demonstrates the temporal distribution of the generated SC pulse with slightly asymmetric structure towards the tailing edge. And the tails on long and short wavelength sides are related to the pedestals near leading and trailing edges. From figure (b), the delay time between different components is less than 2 ps and the temporal phase is

while temporal structure becomes bad.

(a) (b)

56 Photonic Crystals

good uniformity.

bad.

Fig.20 (a) Spectrogram and (b) temporal profile of the generated SC pulse

Fig.20 (a) Spectrogram and (b) temporal profile of the generated SC pulse

Supposing that the delay time between different frequency components is limited to less than Fig.21 (a) Spectrogram, (b) temporal profile and (c) spectral profile after dispersion compensation **Figure 21.** (a) Spectrogram, (b) temporal profile and (c) spectral profile after dispersion compensation

SC after dispersion compensating shows wonderful temporal structures, broad spectral bandwidth, and uniform spectral intensity, which make it well suited for time-resolved applications, especially CARS microspectroscopy. To make the best of the tails in figure (a), spatial light modulator is suggested for the high-order dispersion compensation[86,87]. Supposing that the delay time between different frequency components is limited to less than 40 fs, it is estimated that the useful components still span over 400 nm from 900 to 1300 nm. This SC after dispersion compensating shows wonderful temporal structures, broad spectral bandwidth, and uniform spectral intensity, which make it well suited for time-resolved applications, especially CARS microspectroscopy. To make the best of the tails in figure (a), spatial light modulator is Supposing that the delay time between different frequency components is limited to less than 40 fs, it is estimated that the useful components still span over 400 nm from 900 to 1300 nm. This SC after dispersion compensating shows wonderful temporal structures, broad spectral bandwidth, and uniform spectral intensity, which make it well suited for time-resolved applications, especially CARS microspectroscopy. To make the best of the tails in figure (a), spatial light modulator is suggested for the high-order dispersion compensation [86,87].

40 fs, it is estimated that the useful components still span over 400 nm from 900 to 1300 nm. This

#### **4 Supercontinuum used in infrared nanoscopy 4. Supercontinuum used in infrared nanoscopy**

#### **4.1 Theory of Infrared absorption spectroscopy 4 Supercontinuum used in infrared nanoscopy 4.1. Theory of Infrared absorption spectroscopy**

suggested for the high-order dispersion compensation[86,87].

and is used extensively to detect functional groups.

and is used extensively to detect functional groups.

The infrared **(**IR**)** absorption spectroscopy is an excellent method for biological analyses. It enables the nondestructive, label-free extraction of biochemical information and images[27,88]. The physical basis of IR absorption is that, when an infrared light interacts with the molecule, the vibrations or rotations in a molecule may cause a net change in the dipole moment. The alternating **4.1 Theory of Infrared absorption spectroscopy**  The infrared **(**IR**)** absorption spectroscopy is an excellent method for biological analyses. It enables the nondestructive, label-free extraction of biochemical information and images[27,88]. The physical basis of IR absorption is that, when an infrared light interacts with the molecule, the The infrared (IR) absorption spectroscopy is an excellent method for biological analyses. It enables the nondestructive, label-free extraction of biochemical information and images [27,88]. The physical basis of IR absorption is that, when an infrared light interacts with the molecule, the vibrations or rotations in a molecule may cause a net change in the dipole moment. The alternating electrical field of the radiation interacts with fluctuations in the dipole

electrical field of the radiation interacts with fluctuations in the dipole moment of the molecule. If the frequency of the radiation matches the vibrational frequency of the molecule, then the radiation will be absorbed[89,90]. Every material has its unique IR absorption spectra that can be used as a "fingerprint" for identification. Absorption in the infrared region results in changes in vibrational and rotational status of the molecules. The absorption frequency depends on the vibrational frequency of the molecule, whereas the absorption intensity depends on how effectively the infrared photon energy can be transferred to the molecule, and on the change in the dipole moment that occurs as a result of molecular vibration[91-93]. In a word, a molecule will absorb infrared light only if the absorption causes a change in the dipole moment. Infrared spectroscopy is widely applied to various samples such as liquid, gas, and solid-state matter to identify and to quantify the unknown materials. It is an effective technique to identify compounds

vibrations or rotations in a molecule may cause a net change in the dipole moment. The alternating electrical field of the radiation interacts with fluctuations in the dipole moment of the molecule. If the frequency of the radiation matches the vibrational frequency of the molecule, then the radiation will be absorbed[89,90]. Every material has its unique IR absorption spectra that can be used as a "fingerprint" for identification. Absorption in the infrared region results in changes in vibrational and rotational status of the molecules. The absorption frequency depends on the vibrational frequency of the molecule, whereas the absorption intensity depends on how effectively the infrared photon energy can be transferred to the molecule, and on the change in the dipole moment that occurs as a result of molecular vibration[91-93]. In a word, a molecule will absorb infrared light only if the absorption causes a change in the dipole moment. Infrared spectroscopy is widely applied to various samples such as liquid, gas, and solid-state matter to identify and to quantify the unknown materials. It is an effective technique to identify compounds moment of the molecule. If the frequency of the radiation matches the vibrational frequency of the molecule, then the radiation will be absorbed [89,90]. Every material has its unique IR absorption spectra that can be used as a "fingerprint" for identification. Absorption in the infrared region results in changes in vibrational and rotational status of the molecules. The absorption frequency depends on the vibrational frequency of the molecule, whereas the absorption intensity depends on how effectively the infrared photon energy can be transferred to the molecule, and on the change in the dipole moment that occurs as a result of molecular vibration [91-93]. In a word, a molecule will absorb infrared light only if the absorption causes a change in the dipole moment. Infrared spectroscopy is widely applied to various samples such as liquid, gas, and solid-state matter to identify and to quantify the unknown materials. It is an effective technique to identify compounds and is used extensively to detect functional groups.

Collecting of IR spectra signal usually occurs in two different ways. One is the conventional dispersive-type spectrometer, which employs a grating or a prism to disperse light into individual frequencies, and a slit placed in front of the detector to determine which frequency to reach the detector [94]. However, dispersive-type spectrometer is outdated today. Modern IR spectrometers are the so-called Fourier transform infrared absorption (FTIR) spectroscopy instruments [95,96]. All the FTIR spectrometer operates on the Fourier transform principle, and refers to the manner in which the data is collected and converted from an interference pattern to a spectrum. FTIR spectrometers have progressively replaced the dispersive instru‐ ments for most applications due to their superior speed and sensitivity. They have greatly extended the capabilities of IR spectroscopy and have been applied to many fields that are very difficult or nearly impossible to study by dispersive instruments [97,98]. FTIR can be utilized to measure some components of an unknown mixture and is currently applied to the analysis of solids, liquids, and gases.

The Basic FTIR System as shown in figure 22, there are mainly four parts: Light source, interferometer, sample compartment, and the detector. The beam splitter divides the incoming infrared beam into two beams. One beam reflects off the fixed mirror, and the other beam reflects off the moving mirror, constituting a Michelson interferometer. After reflected from the two mirrors, the divided beams meet each other again at the beam splitter. Approximately 50% of the light passes through the beam splitter and is reflected back along its path by a fixed mirror, and half of the light is reflected by the same beam splitter. The other 50% fraction of the incident light is reflected onto a moving mirror. Light from the moving mirror returns along its original path and half of the light intensity is transmitted through the beam splitter. The two beams meet again at the same beam splitter and then pass through the sample cell, to the infrared sensor. Thus 25% of the incident light from the source reaches the sensor from the fixed mirror and 25% from the moving mirror. As the path length of the two light beams striking the sensor will differ, there will be destructive and constructive interference.

The beam from the moving mirror has traveled a different distance than the beam from the fixed mirror. When the two beams are combined by the beam splitter, an interference pattern is created, since some of the wavelengths recombine constructively and some destructively. This interference pattern is called an interferogram. This interferogram then goes from the

**Figure 22.** The main components of a FTIR

moment of the molecule. If the frequency of the radiation matches the vibrational frequency of the molecule, then the radiation will be absorbed [89,90]. Every material has its unique IR absorption spectra that can be used as a "fingerprint" for identification. Absorption in the infrared region results in changes in vibrational and rotational status of the molecules. The absorption frequency depends on the vibrational frequency of the molecule, whereas the absorption intensity depends on how effectively the infrared photon energy can be transferred to the molecule, and on the change in the dipole moment that occurs as a result of molecular vibration [91-93]. In a word, a molecule will absorb infrared light only if the absorption causes a change in the dipole moment. Infrared spectroscopy is widely applied to various samples such as liquid, gas, and solid-state matter to identify and to quantify the unknown materials. It is an effective technique to identify compounds and is used extensively to detect functional

Collecting of IR spectra signal usually occurs in two different ways. One is the conventional dispersive-type spectrometer, which employs a grating or a prism to disperse light into individual frequencies, and a slit placed in front of the detector to determine which frequency to reach the detector [94]. However, dispersive-type spectrometer is outdated today. Modern IR spectrometers are the so-called Fourier transform infrared absorption (FTIR) spectroscopy instruments [95,96]. All the FTIR spectrometer operates on the Fourier transform principle, and refers to the manner in which the data is collected and converted from an interference pattern to a spectrum. FTIR spectrometers have progressively replaced the dispersive instru‐ ments for most applications due to their superior speed and sensitivity. They have greatly extended the capabilities of IR spectroscopy and have been applied to many fields that are very difficult or nearly impossible to study by dispersive instruments [97,98]. FTIR can be utilized to measure some components of an unknown mixture and is currently applied to the

The Basic FTIR System as shown in figure 22, there are mainly four parts: Light source, interferometer, sample compartment, and the detector. The beam splitter divides the incoming infrared beam into two beams. One beam reflects off the fixed mirror, and the other beam reflects off the moving mirror, constituting a Michelson interferometer. After reflected from the two mirrors, the divided beams meet each other again at the beam splitter. Approximately 50% of the light passes through the beam splitter and is reflected back along its path by a fixed mirror, and half of the light is reflected by the same beam splitter. The other 50% fraction of the incident light is reflected onto a moving mirror. Light from the moving mirror returns along its original path and half of the light intensity is transmitted through the beam splitter. The two beams meet again at the same beam splitter and then pass through the sample cell, to the infrared sensor. Thus 25% of the incident light from the source reaches the sensor from the fixed mirror and 25% from the moving mirror. As the path length of the two light beams

striking the sensor will differ, there will be destructive and constructive interference.

The beam from the moving mirror has traveled a different distance than the beam from the fixed mirror. When the two beams are combined by the beam splitter, an interference pattern is created, since some of the wavelengths recombine constructively and some destructively. This interference pattern is called an interferogram. This interferogram then goes from the

groups.

58 Photonic Crystals

analysis of solids, liquids, and gases.

beam splitter to the sample, where some energy is absorbed and some is transmitted. The transmitted portion reaches the detector. Every stroke of the moving mirror in the interfer‐ ometer equals one scan of the entire IR source spectrum, and individual scan can be combined to give better representation of the actual absorbance of the sample. Each point in the inter‐ ferogram contains information from each wavelength of light being measured. The detector reads information about every wavelength in range of the IR source simultaneously. To obtain the infrared spectrum, the detector signal is sent to the computer, and an algorithm called a Fourier transform is performed on the interferogram to convert it into a spectrum. This transmittance spectrum can be converted to absorbance. The process of collecting the IR spectrum in an FTIR spectrometer is illustrated in figure 23. FTIR spectrometer simultaneously collects spectral data in a wide spectral range of the IR source [99]. In contract, in a dispersive spectrometer, every wavelength across the spectrum must be measured individually. This is a slow process, and typically only one measurement scan of the sample is made. Accordingly, the superior speed and sensitivity of FTIR is the significant advantage over a conventional dispersive spectrometer which measures intensity over a narrow range of wavelengths at a time [100,101].

The FTIR refers to the manner in which the data is collected and converted from an interference pattern to a spectrum. The FTIR spectrometer operates on a principle called Fourier transform. The mathematical expression of Fourier transform can be expressed as

$$\mathbf{F}(\omega) = \int\_{-\infty}^{+\infty} \mathbf{f}(\mathbf{x}) \mathbf{e}^{\mathrm{i}\omega \infty} d\mathbf{x} \tag{38}$$

And the reverse Fourier transform is

$$\mathbf{f}(\mathbf{x}) = \frac{1}{2\pi} \int\_{-\infty}^{+\infty} \mathbf{f}(\mathbf{x}) \mathbf{e}^{-i\omega \infty} \mathbf{d}\omega \tag{39}$$

**Figure 23.** The process of collecting the IR spectrum in an FTIR spectrometer

where *ω* is angular frequency and *x* is the optical path difference. *F* (*ω*) is the spectrum and *f* (*x*) is called the interferogram. It is clear that if the interferogram *f* (*x*), is determined experimentally by using the Michelson interferometer, the spectrum F(ω) can be obtained by using Fourier transform. The detector receives a signal, i.e., the interferogram, which is a summation of all the interferences resulting from the constructive and destructive interaction between each wavelength component and all the others. The two most popular detectors for a FTIR spectrometer are deuterated triglycine sulfate (pyroelectric detector) and mercury cadmium telluride (photon or quantum detector). The measured signal is sent to the computer where the Fourier transformation takes place.

FTIR spectrometer simultaneously collects spectral data of the IR in a wide spectral range. This confers a significant advantage over a dispersive spectrometer which measures intensity over a narrow range of wavelengths at a time. FTIR spectrometers have progressively replaced dispersive instruments for most applications due to their superior speed and sensitivity, and opens up new applications of infrared spectroscopy

#### **4.2. SC used in IR nanoscopy**

FTIR spectroscopy is a widely used analytical tool for chemical identification of inorganic, organic, and biomedical materials, as well as for exploring conduction phenomena [102-104]. Because of the diffraction limit, however, conventional FTIR cannot achieve nano-scaled resolution [105]. Therefore a FTIR system that allows for infrared-spectroscopic nano-imaging is required. The ability to use a spectrally broad source in a FTIR system makes this technique promising at the nano-scale for the chemical identification of unknown nanostructures. Relying on the use of scattering-type scanning near-field optical microscopy [106,107] (s-SNOM), the diffraction limit of a FTIR system is overcome by the use of a sharp atomic force microscope (AFM) tip acting as an antenna to concentrate the incident radiation to nano-scaled volumes [108-110]. Recording of the tip-scattered light thus yields nano-scaled resolved optical images. With the introduction of s-SNOM, the diffraction-limited resolution in FTIR micro‐ scopy can be overcome by several orders of magnitude. Nano-FTIR has interesting application potential in widely different sciences and technologies.

Traditional FTIR spectrometers are based on either thermal IR sources with low brightness and coherence, such as a globar or Hg-lamp, or synchrotron radiation. The synchrotron radiation is complex, expensive, power consuming and always accompanied with intensity fluctuation in spite of high brightness and coherence [29,111]. The advent of non-silica fibers brings us a new IR source with both high brightness and coherence, not to mention the broad bandwidth [31,33,112,113]. Besides fused silica, the PCFs can also be made of other materials, such as chalcogenide and fluoride which show novel characteristics near mid-infrared, such as high transmissivity. Based on these materials, the PCF can shift its ZDW to mid-infrared region, which makes the SC generation in mid-infrared realized easily [114]. And their intrinsic high IR transmissivity and nonlinear coefficient makes it especially suit for mid-infrared SC generation. With these materials, step-index fibers, tapered fibers and photonic crystal fiber can be used to generate the mid-infrared SC [115,116]. Under properly pumping condition, the generated SC shows an ultra-broad spectrum range from 2 µm to 12 µm. These midinfrared supercontinuum possess a broad spectral bandwidth, higher brightness and coher‐ ence [117]. The nano-FTIR combing with such a mid-infrared supercontinuum source, has interesting application potential in widely different sciences and technologies, ranging from the semiconductor industry to nanogeochemistry.

where *ω* is angular frequency and *x* is the optical path difference. *F* (*ω*) is the spectrum and *f* (*x*) is called the interferogram. It is clear that if the interferogram *f* (*x*), is determined experimentally by using the Michelson interferometer, the spectrum F(ω) can be obtained by using Fourier transform. The detector receives a signal, i.e., the interferogram, which is a summation of all the interferences resulting from the constructive and destructive interaction between each wavelength component and all the others. The two most popular detectors for a FTIR spectrometer are deuterated triglycine sulfate (pyroelectric detector) and mercury cadmium telluride (photon or quantum detector). The measured signal is sent to the computer

FTIR spectrometer simultaneously collects spectral data of the IR in a wide spectral range. This confers a significant advantage over a dispersive spectrometer which measures intensity over a narrow range of wavelengths at a time. FTIR spectrometers have progressively replaced dispersive instruments for most applications due to their superior speed and sensitivity, and

FTIR spectroscopy is a widely used analytical tool for chemical identification of inorganic, organic, and biomedical materials, as well as for exploring conduction phenomena [102-104]. Because of the diffraction limit, however, conventional FTIR cannot achieve nano-scaled resolution [105]. Therefore a FTIR system that allows for infrared-spectroscopic nano-imaging is required. The ability to use a spectrally broad source in a FTIR system makes this technique

where the Fourier transformation takes place.

**Figure 23.** The process of collecting the IR spectrum in an FTIR spectrometer

opens up new applications of infrared spectroscopy

**4.2. SC used in IR nanoscopy**

60 Photonic Crystals

In the following paragraphs, we will discuss the process of SC generation with the widely used PCF made of chalcogenide glass. It is necessary to figure out the refraction index and nonlinear response, because they often vary a lot for different materials. The refraction index of bulk As2Se3 glass is approximate to the Sellmeier equation as in Eq. (21). The approximate profile show good agreement with measured data [118].

$$\tan^2(\lambda) = 1 + \lambda^2 \left[ \frac{A\_0^2}{\lambda^2 \cdot A\_1^2} + \frac{A\_2^2}{\lambda^2 \cdot 19^2} + \frac{A\_4^2}{\lambda^2 \cdot 4A\_1^2} \right] \tag{40}$$

with n(*λ*) being the refraction index, *λ* the input wavelength expressed in microns. The relat‐ ed Sellmeier coefficients are *A*<sup>0</sup> =2.234921, *A*<sup>1</sup> =0.24164, *A*<sup>2</sup> =0.347441, *A*<sup>4</sup> =1.308575. By calcula‐ tion, the variations of refractive index and dispersion parameter with wavelength are shown in figure 24(a) with *n* being the refractive index, *n*g the group index. It is clear the group in‐ dex *n*g decreases with wavelength below about 7.2 µm and increases beyond that wave‐ length point. Figure 24(b) exhibits the variation of dispersion with wavelength, where single ZDW point locates in the region from 2 to 14 µm.

nonlinear response is *f <sup>R</sup>* =0.115 and the delayed Raman response function *h <sup>R</sup>*(*t*) can be ex‐ pressed as [118,119]

$$\ln h\_R(t) = \frac{\pi\_1^2 + \pi\_2^2}{\pi\_1 \pi\_1^2} \exp\left(-\frac{t}{\tau\_2}\right) \sin\left(\frac{t}{\tau\_1}\right) \tag{41}$$

with *τ*<sup>1</sup> =23.1 fs, *τ*<sup>2</sup> =195 fs.

point. Figure 24(b) exhibits the variation of dispersion with wavelength, where single ZDW point

point. Figure 24(b) exhibits the variation of dispersion with wavelength, where single ZDW point

Fig.24 (a) Variation of refractive index � and group index ��with wavelength (b) variation of dispersion parameter with wavelength

Fig.24 (a) Variation of refractive index � and group index ��with wavelength (b) variation of dispersion parameter with wavelength

**Figure 24.** (a) Variation of refractive index *n* and group index *n*gwith wavelength (b) variation of dispersion parameter

about 3.0 μm. The ideal structure used in FEM is plotted as in figure 25(a) with a solid core surrounded by four hexagonal rings of air holes. And figure 25(b) shows the dispersion profile of this PCF with single ZDW about 3.72μm . To obtain a larger spectral broadening, this PCF is pump in the anomalous dispersion regime near ZDW, saying 3.75μm at which the nonlinear coefficient is about 2.7m��W�� . The nonlinear coefficient in As2Se3-based PCF is tens to hundreds times larger than that in fused silica due to the higher nonlinear refractive indexn� about ��� �

tens to hundreds times larger than that in fused silica due to the higher nonlinear refractive

W-1

Considering the As2Se3-based PCF has an air hole diameter about 1.2 µm and a hole pitch about 3.0 µm. The ideal structure used in FEM is plotted as in figure 25(a) with a solid core surrounded by four hexagonal rings of air holes. And figure 25(b) shows the dispersion pro‐ file of this PCF with single ZDW about 3.72 µm. To obtain a larger spectral broadening, this PCF is pump in the anomalous dispersion regime near ZDW, saying 3.75 µm at which the

about 3.0 μm. The ideal structure used in FEM is plotted as in figure 25(a) with a solid core surrounded by four hexagonal rings of air holes. And figure 25(b) shows the dispersion profile of this PCF with single ZDW about 3.72μm . To obtain a larger spectral broadening, this PCF is pump in the anomalous dispersion regime near ZDW, saying 3.75μm at which the nonlinear coefficient is about 2.7m��W�� . The nonlinear coefficient in As2Se3-based PCF is tens to hundreds times larger than that in fused silica due to the higher nonlinear refractive indexn� about ��� �

(a) (b)

(a) (b)

nonlinear coefficient is about 2.7 m-1

index n2 about 2.4×10-17m2 / W.

Fig.25 (a) Cross-section structure in simulation and (b) dispersion profile of the PCF

Before modeling the process of SC generation in this As2Se3-based PCF, it is necessary to fig‐ ure out the nonlinear response function R(*t*)=(1 - *f <sup>R</sup>*)*δ*(*t*) + *f <sup>R</sup>h <sup>R</sup>*(*t*), due to the different mate‐ rials. For chalcogenide glasses, the fractional contribution of the Raman response to the total

Fig.25 (a) Cross-section structure in simulation and (b) dispersion profile of the PCF

**Figure 25.** (a) Cross-section structure in simulation and (b) dispersion profile of the PCF

figure out the nonlinear response functionR��� � �1�������� � ������� , due to the different materials. For chalcogenide glasses, the fractional contribution of the Raman response to the total nonlinear response is �� � 0�11� and the delayed Raman response function����� can be

figure out the nonlinear response functionR��� � �1�������� � ������� , due to the different materials. For chalcogenide glasses, the fractional contribution of the Raman response to the total nonlinear response is �� � 0�11� and the delayed Raman response function����� can be

Before modeling the process of SC generation in this As2Se3-based PCF, it is necessary to

Before modeling the process of SC generation in this As2Se3-based PCF, it is necessary to

In the simulation, the injected pulse has a hyperbolic secant profile with with 200 fs pulse

In the simulation, the injected pulse has a hyperbolic secant profile with with 200 fs pulse

duration and 30 W peak power. The temporal and spectral evolutions of input pulse with

duration and 30 W peak power. The temporal and spectral evolutions of input pulse with

����� � ��

����� � ��

���� � ����

���� � ����

. The nonlinear coefficient in As2Se3-based PCF is

� ��� �� �

� ��� �� �

��

��

� ��� � � ��

� ��� � � ��

� (41)

� (41)

Considering the As2Se3-based PCF has an air hole diameter about 1.2 μm and a hole pitch

Considering the As2Se3-based PCF has an air hole diameter about 1.2 μm and a hole pitch

locates in the region from 2 to 14μm .

locates in the region from 2 to 14μm .

62 Photonic Crystals

(a) (b)

(a) (b)

10���m�/W.

10���m�/W.

with wavelength

0

0

expressed as[118,119]

expressed as[118,119]

with �� � ���1 ��, �� � 1�� ��.

with �� � ���1 ��, �� � 1�� ��.

In the simulation, the injected pulse has a hyperbolic secant profile with with 200 fs pulse duration and 30 W peak power. The temporal and spectral evolutions of input pulse with propagating distance are plotted in figure 26. It is clear that the output spectral bandwidth reaches 1.6µm from 3 µm to about 4.7µm after propagating about 20 cm. Further propagation brings in only little spectral broadening on the long-wavelength side, which is induced by intrapulse Raman scattering. In time domain, strongly temporal compression occurs within about 11 cm. After that, there is an obvious pulse broadening and soliton fission. Because different frequency components often possess different velocities, time delay between different frequencies becomes larger with a longer propagation,. propagating distance are plotted in figure 26. It is clear that the output spectral bandwidth reaches 1.6μm from 3μm to about 4.7μm after propagating about 20 cm. Further propagation brings in only little spectral broadening on the long-wavelength side, which is induced by intrapulse Raman scattering. In time domain, strongly temporal compression occurs within about 11 cm. After that, there is an obvious pulse broadening and soliton fission. Because different frequency components often possess different velocities, time delay between different frequencies becomes larger with a longer propagation,.

Fig.26 Temporal (a) and spectral (b) evolutions over propagation distance for input pulse centered at 3.75 μm with peak power 30 W and pulse duration 200 fs **Figure 26.** Temporal (a) and spectral (b) evolutions over propagation distance for input pulse centered at 3.75 µm with peak power 30 W and pulse duration 200 fs

To elaborate on the dynamics of SC generation, different spectral slices are exhibited in figure 27(a) on a logarithmic scale. All the spectrums in figure 27(a) are cut off at -80 dB. In the initial propagation stage, spectral broadening is approximately symmetrical. And the main mechanism of spectral broadening is SPM effect. The SPM-induced spectral broadening has a notable oscillation structure with the outermost peaks possessing the most intense as shown in figure. After around 12 cm, the spectral broadening becomes asymmetric with a larger broadening on the long wavelength side. The out side peaks on long- and short-wavelength side are related to To elaborate on the dynamics of SC generation, different spectral slices are exhibited in figure 27(a) on a logarithmic scale. All the spectrums in figure 27(a) are cut off at -80 dB. In the initial propagation stage, spectral broadening is approximately symmetrical. And the main mecha‐ nism of spectral broadening is SPM effect. The SPM-induced spectral broadening has a notable oscillation structure with the outermost peaks possessing the most intense as shown in figure. After around 12 cm, the spectral broadening becomes asymmetric with a larger broadening on the long wavelength side. The out side peaks on long- and short-wavelength side are related to soliton fission and dispersive wave generation. With further propagation, spectrum on longwavelength side spans little due to weak Raman scattering.

soliton fission and dispersive wave generation. With further propagation, spectrum on

10 W 30 W

50 W 70 W

Fig.27 Spectral evolutions for (a) selected propagating distance and (b) peak power

Although the generated mid-SC has a bandwidth of 1.6 μm, it is not enough for drive all the vibraions simultaneously. A feasible approach is to increase the peak power of pump pulse or shorten the pulse duration. Figure 27(b) shows the numerical spectrums for different peak powers after 30 cm propagation with the same pulse duration 200 fs. The input peak powers are 10 W, 30 W, 50 W and 70 W for the top-left, top–right, bottom-left and bottom–right figures respectively. It

long-wavelength side spans little due to weak Raman scattering.

(a) (b)

4 cm 8 cm

14 cm 20 cm

30 cm 40 cm

larger with a longer propagation,.

(a) (b)

Fig.26 Temporal (a) and spectral (b) evolutions over propagation distance for input pulse centered at 3.75 μm with peak power 30 W and pulse duration 200 fs

figure 27(a) on a logarithmic scale. All the spectrums in figure 27(a) are cut off at -80 dB. In the initial propagation stage, spectral broadening is approximately symmetrical. And the main mechanism of spectral broadening is SPM effect. The SPM-induced spectral broadening has a notable oscillation structure with the outermost peaks possessing the most intense as shown in figure. After around 12 cm, the spectral broadening becomes asymmetric with a larger broadening

soliton fission and dispersive wave generation. With further propagation, spectrum on

To elaborate on the dynamics of SC generation, different spectral slices are exhibited in

propagating distance are plotted in figure 26. It is clear that the output spectral bandwidth reaches 1.6μm from 3μm to about 4.7μm after propagating about 20 cm. Further propagation brings in only little spectral broadening on the long-wavelength side, which is induced by intrapulse Raman scattering. In time domain, strongly temporal compression occurs within about 11 cm. After that, there is an obvious pulse broadening and soliton fission. Because different frequency components often possess different velocities, time delay between different frequencies becomes

Although the generated mid-SC has a bandwidth of 1.6 μm, it is not enough for drive all the **Figure 27.** Spectral evolutions for (a) selected propagating distance and (b) peak power

long-wavelength side spans little due to weak Raman scattering.

vibraions simultaneously. A feasible approach is to increase the peak power of pump pulse or shorten the pulse duration. Figure 27(b) shows the numerical spectrums for different peak powers after 30 cm propagation with the same pulse duration 200 fs. The input peak powers are 10 W, 30 W, 50 W and 70 W for the top-left, top–right, bottom-left and bottom–right figures respectively. It Although the generated mid-SC has a bandwidth of 1.6 µm, it is not enough for drive all the vibraions simultaneously. A feasible approach is to increase the peak power of pump pulse or shorten the pulse duration. Figure 27(b) shows the numerical spectrums for different peak powers after 30 cm propagation with the same pulse duration 200 fs. The input peak powers are 10 W, 30 W, 50 W and 70 W for the top-left, top–right, bottom-left and bottom–right figures respectively. It is clear that the total spectral bandwidth of generated SC is increasing with peak power. And with higher peak power, shorter propagation distance is need to make the spectral broadening to its maximum.

Fig.27 Spectral evolutions for (a) selected propagating distance and (b) peak power

To study the impact of pulse duration on SC generation, peak power of input pulses is kept an constant 30 W in the following simulations. Figure 28(b) exhibit the spectral evolutions for different pulse durations with figures on top-left, top-right, bottom-left, and bottom-right being 100 fs, 200 fs, 300 fs and 500 fs respectively. For pump pulse with shorter duration, the distinct soliton structure is much more clear as in figure (a) and (b) because of a shorter dispersive length *L <sup>D</sup>*. They also have a smaller spectral broadening compared with that in figure (c) and (d) due to a lower pulse energy. However, larger input pulse is associated with longer dispersive length, larger pulse energy and higher soliton order which result in a less clear soliton structure and more complex spectral profile in spite of a larger spectral broaden‐ ing. In the initial propagation, spectral broadening is more with shorter pulse duration because of the larger SPM-induced chirp and its role in FWM.

The spectral broadening with As2Se3-based PCF is much more than that in silica-glass PCF because of the high nonlinearity. To obtain a continuum with both broad spectral bandwidth and good time structure, a shorter propagation distance is suggested because the weak SSFS effect shifts soliton little towards long-wavelength side with finite propagating distance. And with novel fiber structures, the generated SC has a much larger spectral broadening, such as the mid-infrared SC from 2 to 10 µm by Baili *et al* [118]*.*

propagation, spectral broadening is more with shorter pulse duration because of the larger SPM-induced chirp and its role in FWM. Supercontinuum Generation With Photonic Crystal Fibers and Its Application in Nano-imaging http://dx.doi.org/10.5772/59987 65

is clear that the total spectral bandwidth of generated SC is increasing with peak power. And with higher peak power, shorter propagation distance is need to make the spectral broadening to its

constant 30 W in the following simulations. Figure 28(b) exhibit the spectral evolutions for different pulse durations with figures on top-left, top-right, bottom-left, and bottom-right being 100 fs, 200 fs, 300 fs and 500 fs respectively. For pump pulse with shorter duration, the distinct soliton structure is much more clear as in figure (a) and (b) because of a shorter dispersive length�� . They also have a smaller spectral broadening compared with that in figure (c) and (d) due to a lower pulse energy. However, larger input pulse is associated with longer dispersive length, larger pulse energy and higher soliton order which result in a less clear soliton structure and more complex spectral profile in spite of a larger spectral broadening. In the initial

To study the impact of pulse duration on SC generation, peak power of input pulses is kept an

Fig.28 Evolutions over propagation distance for different pulse durations (a) 100 fs; (b) 200 fs; (c) 300 fs; (d) 500 fs **Figure 28.** Evolutions over propagation distance for different pulse durations (a) 100 fs; (b) 200 fs; (c) 300 fs; (d) 500 fs

The spectral broadening with As2Se3-based PCF is much more than that in silica-glass PCF

#### **5. Conclusions and prospects**

maximum.

propagating distance are plotted in figure 26. It is clear that the output spectral bandwidth reaches 1.6μm from 3μm to about 4.7μm after propagating about 20 cm. Further propagation brings in only little spectral broadening on the long-wavelength side, which is induced by intrapulse Raman scattering. In time domain, strongly temporal compression occurs within about 11 cm. After that, there is an obvious pulse broadening and soliton fission. Because different frequency components often possess different velocities, time delay between different frequencies becomes

> Fig.26 Temporal (a) and spectral (b) evolutions over propagation distance for input pulse centered at 3.75 μm with peak power 30 W and pulse duration 200 fs

figure 27(a) on a logarithmic scale. All the spectrums in figure 27(a) are cut off at -80 dB. In the initial propagation stage, spectral broadening is approximately symmetrical. And the main mechanism of spectral broadening is SPM effect. The SPM-induced spectral broadening has a notable oscillation structure with the outermost peaks possessing the most intense as shown in figure. After around 12 cm, the spectral broadening becomes asymmetric with a larger broadening on the long wavelength side. The out side peaks on long- and short-wavelength side are related to soliton fission and dispersive wave generation. With further propagation, spectrum on

long-wavelength side spans little due to weak Raman scattering.

(a) (b)

4 cm 8 cm

14 cm 20 cm

30 cm 40 cm

spectral broadening to its maximum.

of the larger SPM-induced chirp and its role in FWM.

the mid-infrared SC from 2 to 10 µm by Baili *et al* [118]*.*

**Figure 27.** Spectral evolutions for (a) selected propagating distance and (b) peak power

To elaborate on the dynamics of SC generation, different spectral slices are exhibited in

Fig.27 Spectral evolutions for (a) selected propagating distance and (b) peak power

10 W 30 W

50 W 70 W

Although the generated mid-SC has a bandwidth of 1.6 μm, it is not enough for drive all the vibraions simultaneously. A feasible approach is to increase the peak power of pump pulse or shorten the pulse duration. Figure 27(b) shows the numerical spectrums for different peak powers after 30 cm propagation with the same pulse duration 200 fs. The input peak powers are 10 W, 30 W, 50 W and 70 W for the top-left, top–right, bottom-left and bottom–right figures respectively. It

Although the generated mid-SC has a bandwidth of 1.6 µm, it is not enough for drive all the vibraions simultaneously. A feasible approach is to increase the peak power of pump pulse or shorten the pulse duration. Figure 27(b) shows the numerical spectrums for different peak powers after 30 cm propagation with the same pulse duration 200 fs. The input peak powers are 10 W, 30 W, 50 W and 70 W for the top-left, top–right, bottom-left and bottom–right figures respectively. It is clear that the total spectral bandwidth of generated SC is increasing with peak power. And with higher peak power, shorter propagation distance is need to make the

To study the impact of pulse duration on SC generation, peak power of input pulses is kept an constant 30 W in the following simulations. Figure 28(b) exhibit the spectral evolutions for different pulse durations with figures on top-left, top-right, bottom-left, and bottom-right being 100 fs, 200 fs, 300 fs and 500 fs respectively. For pump pulse with shorter duration, the distinct soliton structure is much more clear as in figure (a) and (b) because of a shorter dispersive length *L <sup>D</sup>*. They also have a smaller spectral broadening compared with that in figure (c) and (d) due to a lower pulse energy. However, larger input pulse is associated with longer dispersive length, larger pulse energy and higher soliton order which result in a less clear soliton structure and more complex spectral profile in spite of a larger spectral broaden‐ ing. In the initial propagation, spectral broadening is more with shorter pulse duration because

The spectral broadening with As2Se3-based PCF is much more than that in silica-glass PCF because of the high nonlinearity. To obtain a continuum with both broad spectral bandwidth and good time structure, a shorter propagation distance is suggested because the weak SSFS effect shifts soliton little towards long-wavelength side with finite propagating distance. And with novel fiber structures, the generated SC has a much larger spectral broadening, such as

larger with a longer propagation,.

64 Photonic Crystals

(a) (b)

In this chapter, we mainly described the generalized nonlinear Schro¨dinger equation which is widely used in simulating the pulse evolution in photonic crystal fibers. In this equation, both the higher-order dispersion and nonlinear effects are included, such as self-steepening, optical shock formation and intrapulse Raman scattering. Then predictor-corrector split-step Fourier method are introduced to solve the nonlinear Schro¨dinger equation. Base the equation and numerical method, process of SC generation in a traditional PCF is analyzed in details. Spectral broadening in anomalous dispersions regime is dominated by soliton dynamics, including soliton fission and soliton self-frequency shift, while it is related to SPM effect and optical wave breaking in normal dispersion region. The generated SC has many applications, such as optical frequency metrology, tunable ultrafast fiber laser and nonlinear microscopy for its broad spectrum, high coherence and uniform intensity. Based on the vibrational spectroscopy, CARS microscopy and IR microscopy has been widely used in cellular biology for their label-free and nondestructive imaging. Although the spatial resolution of CARS microscopy is not such high due to the diffraction limit, the APIPD opens up a quite feasible approach to realize nanoscaled imaging. Furthermore, the combination of this CARS nanoscopy and the supercontin‐ because of the high nonlinearity. To obtain a continuum with both broad spectral bandwidth and good time structure, a shorter propagation distance is suggested because the weak SSFS effect shifts soliton little towards long-wavelength side with finite propagating distance. And with novel fiber structures, the generated SC has a much larger spectral broadening, such as the mid-infrared SC from 2 to 10 μm by Baili *et al*[118]*.* **5 Conclusions and prospects**  In this chapter, we mainly described the generalized nonlinear Schro�dinger equation which is uum generated in an all-normal dispersion PCF is of crucial importance to study the fine structures and metabolic dynamics in live cells. Besides, combining with atomic force micro‐ scope, FTIR microscopy has a spatial resolution below 100nm. This nano-FTIR combing with a mid-infrared supercontinuum source generated from the PCF becoming a powerful tool for chemical identification of unknown nanostructures.

#### **Author details**

Shuanglong Liu, Wei Liu and Hanben Niu\*

\*Address all correspondence to: hbniu@szu.edu.cn

Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, College of Optoelectronic Engineering, Shenzhen University, Shenzhen, P.R.China

#### **References**


[10] Washburn, B. R., Ralph, S. E., Lacourt, P. A., Dudley, J. M. *et al.*, Tunable near-infra‐ red femtosecond soliton generation in photonic crystal fibres, ELECTRONICS LET‐ TERS, 2001, 37(25): 1510-1512

uum generated in an all-normal dispersion PCF is of crucial importance to study the fine structures and metabolic dynamics in live cells. Besides, combining with atomic force micro‐ scope, FTIR microscopy has a spatial resolution below 100nm. This nano-FTIR combing with a mid-infrared supercontinuum source generated from the PCF becoming a powerful tool for

Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, College of Optoelectronic Engineering, Shenzhen University,

Coupling in Glass, Physical Review Letters, 1970, 24(11): 584-587

[4] Russell, P., Photonic crystal fibers, science, 2003, 299(5605): 358-362

photonic crystal fibres at 850 nm. *Electronics Letters* 36, 53-55 (2000).

ca, Bell System Technical Journal, 1974, 53(6): 1021-1039

crystal fibers, Opt. Lett., 1998, 23(21): 1662-1664

[1] Alfano, R. R.Shapiro, S. L., Emission in the Region 4000 to 7000 Å Via Four-Photon

[2] Kaiser, P.Astle, H. W., Low-Loss Single-Material Fibers Made From Pure Fused Sili‐

[3] Knight, J. C., Birks, T. A., Russell, P. S. J. & Atkin, D. M., All-silica single-mode opti‐

[5] Mogilevtsev, D., Birks, T. A. & Russell, P. S. J., Group-velocity dispersion in photonic

[6] Ranka, J. K., Windeler, R. S. & Stentz, A. J., Visible continuum generation in air-silica microstructure optical fibers with anomalous dispersion at 800 nm, Opt. Lett., 2000,

[7] Wadsworth, W. J., Knight, J. C., Ortigosa-Blanch, A., Arriaga, J. *et al.* Soliton effects in

[8] Liu, X., Xu, C., Knox, W. H., Chandalia, J. K. *et al.*, Soliton self-frequency shift in a short tapered air-silica microstructure fiber, Opt. Lett., 2001, 26(6): 358-360

[9] Husakou, A. V.Herrmann, J., Supercontinuum Generation of Higher-Order Solitons by Fission in Photonic Crystal Fibers, Physical Review Letters, 2001, 87(20): 203901

cal fiber with photonic crystal cladding, Opt. Lett., 1996, 21(19): 1547-1549

chemical identification of unknown nanostructures.

Shuanglong Liu, Wei Liu and Hanben Niu\*

\*Address all correspondence to: hbniu@szu.edu.cn

**Author details**

66 Photonic Crystals

Shenzhen, P.R.China

25(1): 25-27

**References**


[38] Mamyshev, P. V.Chernikov, S. V., Ultrashort-pulse propagation in optical fibers, Opt. Lett., 1990, 15(19): 1076-1078

[24] Xia, L., Hui, Z., Xiang-Yun, Z., Shi-An, Z. *et al.*, Experimental Investigation on Selec‐ tive Excitation of Two-Pulse Coherent Anti-Stokes Raman Scattering, Chinese Phys‐

[25] Liu Wei, N. H., Diffraction barrier breakthrough in coherent anti-Stokes Raman scat‐ tering microscopy by additional probe-beam-induced phonon depletion, Physical

[26] Sage, J. T., Zhang, Y., McGeehan, J., Ravelli, R. B. G. *et al.*, Infrared protein crystallog‐ raphy, Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics, 2011, 1814(6):

[27] Gillie, J. K.Hochlowski, J., Infrared Spectroscopy, Analytical Chemistry, 2000, 75(12): [28] Huth, F., Schnell, M., Wittborn, J., Ocelic, N. *et al.*, Infrared-spectroscopic nanoimag‐

[29] Reffner, J. A., Martoglio, P. A. & Williams, G. P., Fourier transform infrared micro‐ scopical analysis with synchrotron radiation: the microscope optics and system per‐

[30] Tobin, M. J., Chesters, M. A., Chalmers, J. M., Rutten, F. J. *et al.*, Infrared microscopy of epithelial cancer cells in whole tissues and in tissue culture, using synchrotron ra‐

[31] Kumar, V., George, A., Knight, J. & Russell, P., Tellurite photonic crystal fiber, Opt.

[32] Sanghera, J. S., Shaw, L. B. & Aggarwal, I. D., Applications of chalcogenide glass op‐

[33] Xia, C., Kumar, M., Kulkarni, O. P., Islam, M. N. *et al.*, Mid-infrared supercontinuum generation to 4.5 µm in ZBLAN fluoride fibers by nanosecond diode pumping, Opt.

[35] Blow, K. J.Wood, D., Theoretical description of transient stimulated Raman scattering in optical fibers, Quantum Electronics, IEEE Journal of, 1989, 25(12): 2665-2673 [36] Karasawa, N., Nakamura, S., Nakagawa, N., Shibata, M. *et al.*, Comparison between theory and experiment of nonlinear propagation for a-few-cycle and ultrabroadband optical pulses in a fused-silica fiber, Quantum Electronics, IEEE Journal of, 2001,

[37] Kibler, B., Dudley, J. M. & Coen, S., Supercontinuum generation and nonlinear pulse propagation in photonic crystal fiber: influence of the frequency-dependent effective

ing with a thermal source, Nature materials, 2011, 10(5): 352-356

formance, Review of Scientific Instruments, 1995, 66(2): 1298-1302

diation, Faraday discussions, 2004, 126(27-39

mode area, Appl. Phys. B, 2005, 81(2-3): 337-342

tical fibers, Comptes Rendus Chimie, 2002, 5(12): 873-883

[34] Diament, P. *Wave transmission and fiber optics*. Vol. 74 (Macmillan, 1990).

Express, 2003, 11(20): 2641-2645

Lett., 2006, 31(17): 2553-2555

37(3): 398-404

ics Letters, 2008, 25(6): 2062

Review A, 2011, 83(2): 023830

760-777

68 Photonic Crystals


[69] Bergner, G., Chatzipapadopoulos, S., Akimov, D., Dietzek, B. *et al.*, Quantitative CARS Microscopic Detection of Analytes and Their Isotopomers in a Two-Channel Microfluidic Chip, Small, 2009, 5(24): 2816-2818

[54] Brechet, F., Marcou, J., Pagnoux, D. & Roy, P., Complete analysis of the characteris‐ tics of propagation into photonic crystal fibers, by the finite element method, Optical

[55] Shibata, N., Braun, R. & Waarts, R., Phase-mismatch dependence of efficiency of wave generation through four-wave mixing in a single-mode optical fiber, Quantum

[56] Cristiani, I., Tediosi, R., Tartara, L. & Degiorgio, V., Dispersive wave generation by solitons in microstructured optical fibers, Opt. Express, 2004, 12(1): 124-135

[57] Cheng, J.-X., Jia, Y. K., Zheng, G. & Xie, X. S., Laser-Scanning Coherent Anti-Stokes Raman Scattering Microscopy and Applications to Cell Biology, Biophysical Journal,

[58] Cheng, J.-X.Xie, X. S., Coherent Anti-Stokes Raman Scattering Microscopy: Instru‐ mentation, Theory, and Applications, The Journal of Physical Chemistry B, 2003,

[59] Horie, K., Ushiki, H. & Winnik, F. M., The Interaction of Light with Materials II, Mo‐

[60] Vallée, F.Bogani, F., Coherent time-resolved investigation of LO-phonon dynamics in

[61] Waltner, P., Materny, A. & Kiefer, W., Phonon relaxation in CdSSe semiconductor quantum dots studied by femtosecond time-resolved coherent anti-Stokes Raman

[62] Penzkofer, A., Laubereau, A. & Kaiser, W., High intensity Raman interactions, Prog‐

[64] El-Diasty, F., Coherent anti-Stokes Raman scattering: Spectroscopy and microscopy,

[65] Portnov, A., Rosenwaks, S. & Bar, I., Detection of particles of explosives via back‐ ward coherent anti-Stokes Raman spectroscopy, Applied Physics Letters, 2008, 93(4):

[66] Begley, R. F., Harvey, A. B. & Byer, R. L., Coherent anti‐Stokes Raman spectroscopy,

[67] Cui, M., Bachler, B. R. & Ogilvie, J. P., Comparing coherent and spontaneous Raman scattering under biological imaging conditions, Opt. Lett., 2009, 34(6): 773-775

[68] Zumbusch, A., Holtom, G. R. & Xie, X. S., Three-Dimensional Vibrational Imaging by Coherent Anti-Stokes Raman Scattering, Physical Review Letters, 1999, 82(20):

lecular Photonics: Fundamentals and Practical Aspects, 177-199

scattering, Journal of Applied Physics, 2000, 88(9): 5268-5271

[63] Loudon, R. *The quantum theory of light*. (Oxford university press, 2000).

GaAs, Physical Review B, 1991, 43(14): 12049

ress in Quantum Electronics, 1979, 6(2): 55-140

Vibrational Spectroscopy, 2011, 55(1): 1-37

Applied Physics Letters, 1974, 25(7): 387-390

Fiber Technology, 2000, 6(2): 181-191

2002, 83(1): 502-509

108(3): 827-840

70 Photonic Crystals


4142-4145

Electronics, IEEE Journal of, 1987, 23(7): 1205-1210


[98] Helm, D., Labischinski, H., Schallehn, G. & Naumann, D., Classification and identifi‐ cation of bacteria by Fourier-transform infrared spectroscopy, Journal of general mi‐ crobiology, 1991, 137(1): 69-79

[83] Hooper, L. E., Mosley, P. J., Muir, A. C., Wadsworth, W. J. *et al.*, Coherent supercon‐ tinuum generation in photonic crystal fiber with all-normal group velocity disper‐

[84] Hartung, A., Heidt, A. M. & Bartelt, H., Design of all-normal dispersion microstruc‐ tured optical fibers for pulse-preserving supercontinuum generation, Opt. Express,

[85] Tomlinson, W., Stolen, R. H. & Johnson, A. M., Optical wave breaking of pulses in

[86] Dorrer, C., Salin, F., Verluise, F. & Huignard, J., Programmable phase control of fem‐ tosecond pulses by use of a nonpixelated spatial light modulator, Opt. Lett., 1998,

[87] Karasawa, N., Li, L., Suguro, A., Shigekawa, H. *et al.*, Optical pulse compression to 5.0 fs by use of only a spatial light modulator for phase compensation, JOSA B, 2001,

[88] Knoll, B.Keilmann, F., Near-field probing of vibrational absorption for chemical mi‐

[89] Chabal, Y. J., Surface infrared spectroscopy, Surface Science Reports, 1988, 8(5):

[91] Chappell, J., Bloch, A., Bryden, W., Maxfield, M. *et al.*, Degree of charge transfer in organic conductors by infrared absorption spectroscopy, Journal of the American

[92] Zhu, X., Suhr, H. & Shen, Y., Surface vibrational spectroscopy by infrared-visible

[93] Kahn, A. H., Theory of the infrared absorption of carriers in germanium and silicon,

[94] Iwata, K.Hamaguchi, H.-O., Construction of a Versatile Microsecond Time-Resolved

[95] VAN DE WEERT, M., Fourier Transform Infrared, Methods for Structural Analysis

[96] Rochow, T. G.Tucker, P. A. in *Introduction to Microscopy by Means of Light, Electrons, X*

[97] Arrondo, J. L. R., Muga, A., Castresana, J. & Goñi, F. M., Quantitative studies of the structure of proteins in solution by Fourier-transform infrared spectroscopy, Prog‐

sion, Opt. Express, 2011, 19(6): 4902-4907

croscopy, Nature, 1999, 399(6732): 134-137

Chemical Society, 1981, 103(9): 2442-2443

Physical Review, 1955, 97(6): 1647

of Protein Pharmaceuticals, 2005, 3(131

*Rays, or Acoustics* 257-264 (Springer, 1994).

[90] Conley, R. T. *Infrared spectroscopy*. (Allyn and Bacon, 1972).

sum frequency generation, Physical Review B, 1987, 35(6): 3047

Infrared Spectrometer, Appl. Spectrosc., 1990, 44(9): 1431-1437

ress in biophysics and molecular biology, 1993, 59(1): 23-56

nonlinear optical fibers, Opt. Lett., 1985, 10(9): 457-459

2011, 19(8): 7742-7749

23(9): 709-711

72 Photonic Crystals

18(11): 1742-1746

211-357


**Chapter 3**
