Quantitative chemical imaging with background-free multiplex coherent anti-Stokes Raman scattering by dual-soliton Stokes pulses.

Coherent anti-Stokes Raman microscopy (CARS) is a quantitative, chemically specific, and label-free optical imaging technique for studying inhomogeneous systems. However, the complicating influence of the nonresonant response on the CARS signal severely limits its sensitivity and specificity and especially limits the extent to which CARS microscopy has been used as a fully quantitative imaging technique. On the basis of spectral focusing mechanism, we establish a dual-soliton Stokes based CARS microspectroscopy and microscopy scheme capable of quantifying the spatial information of densities and chemical composition within inhomogeneous samples, using a single fiber laser. Dual-soliton Stokes scheme not only removes the nonresonant background but also allows robust acquisition of multiple characteristic vibrational frequencies. This all-fiber based laser source can cover the entire fingerprint (800-2200 cm-1) region with a spectral resolution of 15 cm-1. We demonstrate that quantitative degree determination of lipid-chain unsaturation in the fatty acids mixture can be achieved by the characterization of C = C stretching and CH2 deformation vibrations. For microscopy purposes, we show that the spatially inhomogeneous distribution of lipid droplets can be further quantitatively visualized using this quantified degree of lipid unsaturation in the acyl chain for contrast in the hyperspectral CARS images. The combination of compact excitation source and background-free capability to facilitate extraction of quantitative composition information with multiplex spectral peaks will enable wider applications of quantitative chemical imaging in studying biological and material systems.


Introduction
Raman microscopy has been extremely successful in visualizing the chemical changes and distribution of species within soft and complex materials, including polymeric and biological samples, as a powerful label-free technique [1,2]. It especially provides a high level of sensitivity and specificity to identify disease related chemical changes [3]. The conventional spontaneous Raman imaging methods bear two major weaknesses, including weak signal intensity and an often strong fluorescence background in the same spectral range. The discrimination of subtly different states of cells and tissues relying on multiple Raman peaks within the weakly scattering fingerprint region (<1,800 cm −1 ) generally requires long integration time in the order of seconds per spectrum [4], seriously limiting the use of Raman spectroscopy in high-speed imaging. To bolster the inherently weak Raman scattering process, coherent Raman techniques have been developed, which coherently populate selected vibrational states of molecules through their nonlinear response to multiple pulsed laser fields. Successful techniques for microscopy include coherent anti-Stoke Raman scattering (CARS) microscopy [5] and stimulated Raman scattering (SRS) microscopy [6]. In both cases, video-rate imaging of single Raman bands can be achieved with narrow spectral bandwidth [7,8], but suffers from limitations in breadth and speed of laser tuning rates. With respect to the multiplex or broadband CARS or SRS spectroscopy and imaging, there is considerably more information contained throughout an entire spectrum than at just a single frequency within that spectrum [9,10]. This advantage is more useful for qualitative discrimination of many overlapping peaks, in general, and more amenable to quantitative analysis of variations in the spatial distribution of different molecules in the sample.
In order to further reconcile the bandwidth and resolution conflict of broadband lasers in the broadband (multiplex) CARS, Raman spectrum can be acquired by scanning the time delay between the two pulses when they are appropriately chirped. This is known as the spectral focusing approach [11], which concentrates most of the optical power into a single Raman vibrational mode, thereby providing a flexible way to optimize the spectral resolution and signal-to-noise ratio. However, a well-known background from nonresonant (NR) CARS contributions severely limits the sensitivity and specificity of multiplex CARS, which does not carry chemically specific information. The NR background can further distort the resonant signal of interest, which has limited the extent to which CARS microscopy has been used as a fully quantitative imaging technique. It especially renders quantitative CARS hyperspectral imaging on target samples exceedingly difficult in the fingerprint spectral region where molecules present their unique vibrational signature with lower cross section than the frequently used C-H stretch [9]. Some approaches have been devised to minimize the NR background contribution to the measured multiplex CARS spectrum [12][13][14][15][16][17]. Frequency modulation (FM) CARS, has proven to be an effective way to eliminate the nonresonant background, relying on the different spectral shapes of the resonant CARS and NR signals [15][16][17]. If one can perform a differential measurement with slightly shifted excitation frequencies, the nonresonant background can be eliminated in the FM-CARS scheme. In the previously demonstrated narrow FM-CARS microscopy methods, one needs either two different pump lasers that are switched by a Pockels cell [15] or a sophisticated optical parametric oscillator (OPO) that can modulate the frequency of the output pulses [16]. With respect to the broadband CARS, Chen et al. firstly proposed a new FM-CARS version by combination of fast frequency modulation of time-delay between pump and Stokes pulses and lock-in signal detection for quantitative analysis of unsaturation composition in mixtures of fatty acids [17]. Hence, sophisticated laser systems or pulse-scan system were required to perform the necessary fast frequency modulation in beam-scanning microscopy for the abovementioned methods, yielding a complex and costly setup. The challenge is then to find a way of suppressing the NR signals of multiplex CARS, while beating the inaccessibility and complexity of conventional modulation system. Very recently, we established a new all-fibergenerated dual-soliton Stokes based scheme for background-free CARS under the spectral focusing mechanism [18], without any additional modulation. This scheme, termed as DS-CARS, takes full advantage of the all-fiber-based source technology which could be used in less favorable environments while being compact and inexpensive.
In this article, we establish a dual-soliton Stokes based multiplex and background-free CARS method capable of quantifying the spatial information of densities and chemical composition within inhomogeneous samples. We show that our method not only removes the NR background but also generates Raman-like CARS signals at multiple characteristic vibrational frequencies in the fingerprint region. We demonstrate that the degree of unsaturation of mixture lipid droplets can be determined with a well-defined correlation between the vibration ratio of C = C/CH 2 and the chemical concentrations. We further demonstrated that, the spatially inhomogeneous distribution of lipid droplets can be quantitatively visualized using the degree of acyl chain unsaturation for contrast in the hyperspectral CARS images. We believe that the powerful combination compact excitation source and background-free capability can greatly simplify multiplex CARS spectroscopy measurements and extend the practicality of quantitative CARS microscopy for studying material and biological systems.

Dual-soliton pulse generation
The current gold-standard laser system for CARS microscopy is synchronized picosecond (ps) mode-locked solid-state oscillators [19], or synchronously pumped ps OPOs [20]. Free space lasers generally require a stable environment or active feedback, whereas an all-fiber source could be used in less favorable environments while being more compact and inexpensive. Here we show, by exploring the potential of high birefringence of polarizationmaintaining photonic crystal fiber (PM-PCF), that a new version of fiber-generated Stokes pulse can be harnessed for background-free CARS spectroscopy and microscopy [18]. This birefringence has the potential to simultaneously generate two continua (two soliton pulses) with orthogonal polarizations, allowing for an extra degree of freedom in tuning the properties of the supercontinuum. Governed by the coupled vector nonlinear Schrödinger equations [21], numerically simulated results in Fig. 1 qualitatively display the generated spectra and soliton pulse profiles. When the input polarization isn't parallel to the fast or slow axes, contour plots of the spectra provides an unambiguous proof that the generated spectrum is a linear combination of two continua generated separately along the two principal axes. Indeed, spectral components along one axis are not coupled to those in the other axis. This is attributable to the different dispersion characteristics of the two eigenpolarizations, thus leading to the temporal walk-off of this dual-soliton τ 0 . It is worth noting that both frequency separation and temporal interval between dual-soliton pulses can be easily adjusted by rotating the input polarization. Soliton 2 shifts just about 1.5nm (~10 cm −1 ) relative to soliton 1, towards to shorter wavelength, and the temporal interval is about 1.3ps. Due to the balance between anomalous dispersion and the Kerr nonlinearity [22], both soliton 1 and 2 are nearly transform-limited.

Spectral focusing CARS
In the CARS process, molecules are excited by two input light fields, a pump pulse and a Stokes pulse, whose frequency difference is chosen to match the frequency of Raman active transitions. If transform limited femtosecond pulses are used, the signal size increases at the expense of reduced spectral resolution. A straightforward way, known as spectral focusing method is to equally apply a linear temporal chirp on both pump and Stokes pulses to resolve the conflict. For linearly chirped pulses under the stationary phase approximation, the pump and Stokes electric field can be described by: where E P and E S are the complex field envelopes and E P0 and E S0 represent the absolute field amplitude. ω P0 and ω S0 are the center frequency of pump and Stokes pulses and Δt is the delay between them. α is the common chirp parameter. The frequency components of the two beams are thereby stretched in time, thus leading to the instantaneous frequencies: ω P (t) = ω P0 + 2αt and ω S (t) = ω S0 + 2α(t-Δt), respectively. Let us consider the nonlinear vibrational excitation field, A(t), by chirped pump and Stokes pulses [11,17,23]: Obviously, equal linear chirp on both pump and Stokes pulses creates a constant, narrowed, instantaneous frequency difference (IFD) as Ω IFD = ω p0 -ω s0 + 2αΔt, which drives a single Raman coherence [23][24][25][26]. It also offers a route at tuning the IFD within the spectral width of the chirped pulses by simply adjusting their relative time delay, allowing to perform broadband CARS spectroscopy without tuning the laser sources. Most importantly, one can achieve the best compromise between spectral resolution and signal strength for CARS microscopy when the time durations of the pump and Stokes pulses are similar to the vibrational dephasing time of molecules (a few picoseconds) [23].

Dual-soliton CARS
The temporal walk-off dual-soliton pulses can be intrinsically applied to the spectral focusing scheme, which performs CARS spectroscopy by sequentially sweeping the interpulse delay between pump and Stokes pulses. In the spectral focusing mechanism, a single Raman level can be excited by multiple pairs of pump and Stokes frequency components, accompanying with NR nonspecific four-wave-mixing signal at the same frequency as schematically illustrated in Fig. 2(a) and 2(b). The total CARS response can be described as the sum of resonant third-order nonlinear polarization (3) R P and the NR one (3) NR P . The detected CARS signal I CARS (ω as ) at anti-Stokes frequency (ω as ) is therefore proportional to [9]: (3) NR P is real, chemically-nonspecific typically assumed to be slowly varying in frequency (although not necessarily constant). As described by the last term of Eq. (4), the non-resonant background electric field distorts the CARS spectrum due to its constructive and destructive interference with the resonant vibrational contribution, , on the red and blue sides of the Raman peak, respectively [26]. The resonant nonlinear polarization (3) R ( ) as P ω can be further written in a more tractable form [9]: R as s as p as pr as where  and ⊗ are the cross-correlation and convolution operations, respectively; E p (ω as ), E s (ω as ) and E pr (ω as ) are the pump, Stokes and probe electric field, respectively; (3) R χ is the nonlinear susceptibility for the resonant components. ω p , ω s and ω pr are the pump, Stokes and probe frequencies, respectively;A j , Ω j , and Γ j are the amplitude, frequency, and line width of the j th vibrational mode and N is the number of scatterers per unit volume. In the Dual-soliton CARS scheme, the Raman level Ω can be excited twice while scanning the interpulse delay between the pump and the Stokes beam shown in Fig. 2(a). So, we can acquire CARS signals at different time-delay: I CARS1 (ω as , ω as1 ) and I CARS2 (ω as , ω as2 ). By combining the time and frequency picture in Fig. 2(b), the effect of the frequency separation between the two soliton pulses on the vibrational coherence and especially on the signal generation can be promptly understood. As revealed by Fig. 2(b), the delay interval τ 1 is not equal to the temporal interval of two TL soliton pulses τ 0 , since soliton 2 shifts to higher frequency by a slight amount of Δω. Then the excellent linear correlation between the time delay and frequency components of the chirped pump (or probe) pulse ensures the same amount of shift between ω pr1 and ω pr2 , further leading to the frequency difference of the two anti-Stokes signals ω as2 -ω as1 = Δω as illustrated in Fig. 2(a). This frequency modulation with two excitation wavelengths shifted by an amount comparable to the bandwidth of the measured Raman bands has been demonstrated as an effective way to eliminate the large backgrounds [17,24]. Dual-soliton CARS scheme is therefore naturally endowed with this capability to suppress the NR signals. In the scenario of weak or dilute analytes, We can then exhibit a difference CARS spectrum between I CARS1 (ω as , ω as1 ) and I CARS2 (ω as , ω as2 ).
CARS as NR R as as R as as The DS-CARS response is then linear in concentration amplified by the NR background without any offset from the strong NR contribution. We will further elaborate this linear correlation with simulations in Fig. 3. Two simulated CARS signals generated by soliton 1 and 2 with an appropriate frequency separation are both distorted by the nonresonant background in Fig. 3(a). However, the output lineshape of the difference spectrum between the two distorted CARS signals looks close to that of spontaneous Raman scattering. More importantly, the DS-CARS become linearly proportional to sample concentration. This linear dependence can be clearly demonstrated in the case of different Raman peak width (5, 15 and 10 cm −1 ) and frequency separations (4, 8, 10, 12 and 16 cm −1 ) in (b) and (c), respectively. This is due to the fact that our DS-CARS measures the cross term , which allows a broader application of quantitative DS-CARS spectroscopy and microscopy. Besides, a smaller frequency separation between soliton 1 and 2 generally leads to the increase of the spectral resolution of the DS-CARS, yet the decrease of the intensity as shown in Fig. 3(c). In order to maximize the spectral resolution of DS-CARS signals with minimum loss of its peak intensity, the frequency separation between the two Stokes pulses should be matched with the Raman mode to be probed.

Experimental setup
The implementation of DS-CARS proceeds as presented in Fig. 4(a) [18]. A home-built Ybdoped fiber laser with an output of 2W average power at 100 MHz repetition rate was used as the primary source. This laser output was split into the pump and the Stokes paths. For efficient Stokes generation, about 650mW is coupled into a 192 cm highly nonlinear PM-PCF (coupling efficiency~30%).The continuum redshifts to about 1280nm so as to cover most of the fingerprint region from 1150 to 1800 cm −1 , by adjusting the input power to the PCF using a half wave-plate and polarizing beam splitter combination. By rotating the input polarization, the generated spectrum at the fiber output actually consists of a superposition of the spectra generated independently by the two eigenpolarization modes, which is termed as Soliton 1 (red line) and Soliton 2 (green line) in Fig. 4(b). These results are in good agreement with the simulated contour plots in Fig. 1. The temporal walk-off (τ 0 ) and the frequency separation (ω) between two soliton pulses are 1.4ps and 10cm −1 , respectively. SF57 glass rods with different length are added to ensure that the same amount of linear positive chirp are applied to both pump and Stokes beam. After traveling through these high-index glass rod, the pulse durations of pump and Stokes pulses are measured to be 2.2ps and 1.4ps, respectively with an autocorrelator (APE-Berlin). It is worth noting that both the frequency separation (ω) and the degree of linear chirp are chosen considering the Raman line width of the mode to be probed. A motorized delay line is inserted into the pump beam to scan the delay between the two pulse trains. These two beams are overlapped in space by a dichroic mirror and sent into a custom-made scanning microscope with a near-infrared optimized excitation objective lens (NA = 0.65, LCPLN50XIR, Olympus) and the collection objective is the same as the excitation one. The signal at anti-Stokes frequency is detected in the forward direction by means of an imaging spectrograph (IsoPlane160, Princeton Instruments) attached with a backilluminated, deep depletion CCD (PIXIS1024BR, Princeton Instruments). For imaging, the sample is raster scanned. Most of CARS measurements in this work are performed with the pump and Stokes powers of 40.0mW and 12.0mW (I soliton 1 = 6.5mW, I soliton 2 = 5.5mW), respectively. Prominent spontaneous Raman peaks of ethanol are used to calibrate the linear coefficient between the optical delay and Raman wavenumber. Linear fitting of these Raman peaks to corresponding optical delays in DS-CARS measurements results in a calibration curve which can be used to determine Raman shift at any interpulse delay (Fig. 4(c)).

Sample preparation
For multiplex DS-CARS acquisition, a kind of polymer film, polyethylene glycol terephthalate (PET) with the thickness of 200μm was directly used. The fatty acids used in this work are fish oil (Nature's Bounty, USA), olive oil (Betis, Spainsh) and oleic acid (Sigma-Aldrich) and used without further purification. For quantitative identification of lipid unsaturation, different proportions of fish oil were directly incorporated into olive oil. For quantitative DS-CARS imaging, oil droplets are prepared in the following way [17]: 0.1 g of each oil is mixed with 0.1 g of Triton X-100 (Sigma-Aldrich) and stirred for 10 min, and 10 mL of water is added to the mixture, which is stirred further for 2 h. Separately prepared olive and fish oil droplets are mixed afterward. For all the liquid samples, a drop of solution was pipetted inside a 120μm thick imaging spacer (GraceTMBio-Lab SecureSealTM) glued on a glass slide (1 mm thick) in order to create a chamber, which was sealed by a second coverslip (0.13mm thick) on top.

Multiplex DS-CARS spectra of polymer film
We first address the background-free multiplex CARS spectra in the dual-soliton scheme. Figure 5(a) shows the measured CARS spectrograms of PET film with the thickness of 200μm. We are bringing clear evidences that the two soliton Stokes pulses together with one pump pulse can excite a Raman level subsequently while scanning the interpulse delay between the pump and the Stokes beam, also as illustrated in Fig. 2(a) The true CARS spectra are extracted from the three-dimensional spectrograms as the amplitude along the time-delay by integrating the total power of each anti-Stokes spectrum from the spectrometer. As shown in Fig. 5(b), CARS signals produced from both soliton 1 and 2 present distorted and typically dissipative spectral lineshape, which is ascribed to the strong NR backgrounds and wellunderstood in the CARS community. In the DS-CARS scheme the pump pulse scans the two Stokes soliton pulses at different delay-time and the delay interval is about τ 1 in (b). It is worth noting that τ 1 is not equal to the temporal interval of two soliton pulses τ 0 , because the frequency separation of soliton 1 and 2 translates to difference of excited Raman frequency. After shifting the CARS spectra of soliton 2 in the time domain by the amount of τ 0 , the underlying broad NR background can be subtracted from the difference between these two CARS signals, and its output is almost background-free and more resolvable (Fig. 5(c)). DS-CARS spectra of PET is almost identical to its spontaneous counterparts except for relative intensity differences at low and high wavenumbers due to decreased pulse overlap [10]. This intensity variation is always constant and can be removed by calibration against the known spontaneous Raman signals if needed. Figure 5(d) shows the intensity calibrated DS-CARS spectra, in which the spectral positions and FWHMs of different Raman resonance match very well with the spontaneous spectrum, allowing quantitative analysis in CARS spectroscopy and microscopy.

Quantitative chemical identification of fatty acids
Understanding the diverse biological functions of lipid droplets and their role in the development of metabolic diseases requires the number of double C = C bonds within lipid acyl chains, i.e., the degree of unsaturation. Clearly, quantitative spatial and temporal information on the degree of lipid unsaturation would greatly contribute to our understanding of metabolic diseases and the influence of nutritional factors [28]. Oleic acid is a monounsaturated omega-9 fatty acid with 18-carbon chain and one C = C double bond (18:1). According to the manufacture label of the olive oil used in this experiment, it contains 14, 7, and 79% of saturated, polyunsaturated, and monounsaturated fatty acid, respectively. The major component in olive oil is oleic acid. For the fish oil, 30% of its content is a mixture of eicosapentanoic acid (20:5), docosahexanoic acid (22:6) and other polyunsaturated fatty acids according to the manufacture's label. Figure 6 shows the DS-CARS and spontaneous Raman spectra of fish oil, olive oil and oleic acid. The originally measured CARS spectra by soliton1 and 2 are severely distorted by the nonresonant background in Fig. 6(a), while DS-CARS produces Raman-like CARS signals in (b). Since these spectra contain rich information from multiple characteristic vibration peaks such as olefinic C = C and methylene CH 2 deformations, background-free DS-CARS is more likely to identify molecular structures by quantitative spectral analysis. Besides, difference based DS-CARS is even more resolved than the spontaneous Raman spectra in the crowded Raman region at 1260 and 1300 cm −1 . The spectral resolution of DS-CARS can be estimated at 15 cm −1 from the two Raman peaks. In Fig. 6(b), the 1260 and 1300 cm −1 peaks are assigned to in-plane olefinic hydrogen bending (1260 cm −1 ) and methylene twisting deformation (1300 cm −1 ) [17].The former increase with the numbers of C = C double bonds; the latter increase with the number of C-C single bonds. The resonances at 1650 and 1450 cm −1 originated from the v(C = C) stretching and δ(CH 2 ) deformation vibrations, respectively. The signal at 1650 cm −1 was directly proportional to the concentration of C = C bonds, and the1450 cm −1 resonance can be used as a measure of total fatty-acid concentration [26]. Obviously, the peak ratios of 1260 and 1300 cm −1 (I 1260 /I 1300 ) and 1450 and 1650 cm −1 ((I 1450 /I 1650 )) provide a quantitative measure for the degree of lipid-chain unsaturation: C = C/CH 2 and CH 2 /C = C, respectively. Low unsaturation level of olive oil shows I 1260 /I 1300 <1 and I 1450 /I 1650 >1. On the other hand, high unsaturation level of fish oil exhibits I 1260 /I 1300 >1 and I 1450 /I 1650 <1 as shown in Fig. 7(a).
In order to verify the linear dependence on the sample concentration experimentally, as we demonstrated in Fig. 4, we further test the variation of peak ratios (I 1260 /I 1300 and I 1450 /I 1650 ) with the fish fraction in the fish/olive mixture oil droplets. On the basis of the theoretical analysis in Eq. (7), we can express the peak ratio in the mixture oil as  (11). It is apparent from Fig. 7(b) that a good correlation of two peak ratios (I 1260 /I 1300 and I 1450 /I 1650 ) and the fish fraction in fish/olive mixture can be established. More importantly, when the two experimental data are compared with theoretical fitting curves of Eq. (8) and (9), it is noticeable the proximity between experimental results and theoretical plot. Indeed, the calculated average relative error between the theoretical and experimental results are 0.007% and 0.016% and the maximum difference are of 3.3% and 2.9% in the entire range of concentration variation for I 1260 /I 1300 and I 1450 /I 1650 , respectively. These features provide rich chemical information about acyl chains of fatty acids and lipids to determine the unsaturation level. They serve as quantitative means of chemical selectivity in oil droplet imaging in the next sections. Background-free multiplex CARS imaging can be used to quantitatively map the spatial variations in the level of lipid unsaturation and acyl-chain order of different oil droplets. Figure 8(a) and 8(b) show the DS-CARS images of a mixture of olive and fish oil droplets at 1260 and 1300 cm −1 , respectively. The main difference between their chemical structures is the degree of unsaturation in their acyl chains. Olive oil consists of oleic acid (18:1) primarily (79%), while fish oil contains a large amount (30%) of polyunsaturated fatty acids such as EPA (20:5) and DHA (22:6). Note that both oils possess vibrational peaks at 1260 and 1300 cm −1 . Because of the high unsaturation level, one can see larger signals from fish oil droplets than the olive ones at 1260 cm −1 at droplet A in Fig. 8(a). On the other hand, the olive oil droplet B is a little brighter than the fish ones in Fig. 8(b). However, brightness contrast in microscopy may lead to false chemical identification. For example, the droplet marked with arrows in Fig. 8(a) is a fish oil droplet positioned vertically off from the laser focus. Its apparent signal strength at 1260 cm −1 is similar to that of the olive oil droplet. The olive oil droplet marked with arrows in Fig. 8(b) also suffers from this problem. So, spatial and temporal information of lipid droplets can be easily obfuscated by the qualitative visual inspection. Background-free multiplex CARS has the potential to be used to quantify the degree of lipid-chain unsaturation and then the two oil droplets in microscopy, using the peak ratio (I 1260 /I 1300 ). We compare the CARS images at 1260 cm −1 (Fig. 8(a)) and 1300 cm −1 (Fig.  8(b)) and take the ratio image as shown in Fig. 8(c). It can quantitatively and clearly distinguish these two different oil droplets. We find that the intensity ratios of 1260 and 1300 cm −1 peaks (I 1260 /I 1300 ) are ∼1.6 and ∼0.48 for fish and olive oils, respectively. Therefore, one can also clearly identify each droplet by the peak intensity ratio like I 1260 /I 1300 , avoiding false interpretation of local chemical components. It is worth noting that the spot by the blue arrow, which is actually a fish droplets shows green color with low signal strength in (g), because it has small size and is positioned slightly off from the laser focus. The validity of the background-free multiplex CARS microscopy can be seen from the quantitative images of both Fig. 8(c) and 8(g) and the results are in good agreement with the previously reported work [17,28].
We can also implement microspectroscopy while taking CARS images. Figure 8(d) and 8(h) show DS-CARS spectra taken at the droplets A, B and C, D. We position the laser focus at the point of interest and scan the time delay to obtain these spectra. Thus, we are able to not only acquire chemical images at several important frequencies in real-time but also perform microspectroscopy at a few positions of interest. Combination of multiple frequency imaging and microspectroscopy will be a powerful tool in the study of complex samples such as cells and tissues, offering more chemical sensitivity and selectivity.

Conclusion
The combination of the high specificity of background-free CARS with rich information from multiple characteristic vibration peak in multiplex CARS spectra allows for excellent capability in quantitative label-free imaging of different molecules. The realization of such a powerful combination is obtained through dual-soliton pulses based on a single fiber laser source under the spectral focusing mechanism. All-fiber laser source for CARS application is advantageous over the free space lasers in terms of its simplicity, flexibility and lower cost. Dual-soliton generation just takes advantage of a piece of PM-PCF fiber without any change of laser cavity or extra modulation. The experimental implementation of DS-CARS is straightforward and could be used in less favorable environments while being more compact. All these features are unique to the DS-CARS scheme and differ from previously reported efforts aimed at suppressing the NR background signals in CARS. In addition, advanced spectral focusing was proposed and established recently by further exploring the potential of this technique [25,26], and we believe that this improvement will benefit the DS-CARS method.
Since the DS-CARS method removes the nonresonant background and measures Ramanlike vibrational signals in both microspectroscopy and microscopy, we can perform label-free imaging at multiple characteristic vibrational frequencies to obtain more information about the distributions of multiple chemical species. As we demonstrate in this work, the combination of ratio imaging and in situ background-free microspectroscopy is a powerful tool for study of inhomogeneous samples. The strength of C = C double bond stretching over CH 2 deformation vibrations (I 1260 /I 1300 and I 1450 /I 1650 ) can quantify the degree of lipid-chain unsaturation in the fatty acids mixture, benefiting the reveal of the function of lipid droplets in health and disease. Besides, multiplex DS-CARS imaging is further identified to quantify different oil droplets by the peak intensity ratio. This distinct vibrational markers (I 1260 /I 1300 and I 1450 /I 1650 ) can be used to quantitatively map the spatial variations in the level of lipid unsaturation and acyl-chain order of different lipid droplets within a single cell, or even within individual lipid droplets. Although the DS-CARS capability has been demonstrated by use of a spectrometer with a CCD, it is a parallel process and significantly limit the imaging speed. However, enlarging the temporal interval between the dual-soliton pulses by a longer PCF can allow the utilization of a single detector, such as APD or PMT to accelerate the imaging speed. Besides, our method allows intrinsically feasible tunability of excitation for CARS based on the dual-soliton supercontinuum generation, and we can obtain fingerprint vibrational spectroscopy and images over a large spectral range (800-2200 cm −1 ) with a spectral resolution of 15 cm −1 . This vibrational region has profound chemical information and can help to explore complicated biochemical processes in cells and tissues.