Investigation of the long wavelength limit of soliton self-frequency shift in a silica fiber

We explore the long wavelength limit of soliton self-frequency shift in silica-based fibers experimentally and using numerical simulation. We found that the longest wavelength soliton generated by soliton self-frequency shift is approximately 2500 nm because the soliton loses its energy rapidly at wavelength beyond 2400 nm due to material absorption by silica and water. We demonstrate 1580-2520 nm wavelength-tunable, high-pulse energy soliton generation using soliton self-frequency shift in a large-mode-area silica fiber pumped by a compact fiber source. Soliton pulses with pulse width of ~100 fs and pulse energy up to 73 nJ were obtained. Second harmonic generation of the solitons enables a wavelength-tunable femtosecond source from 950 nm to 1260 nm, with pulse energy up to 21 nJ. Using such energetic pulses, we demonstrate in vivo two-photon excited fluorescence imaging of vasculature and neurons in a mouse brain at wavelength beyond the tuning range of a modelocked Ti:Sapphire lasers. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement OCIS codes: (320.7110) Ultrafast nonlinear optics; (060.2380) Fiber optics sources and detectors; (110.0180) Microscopy. References and links 1. X. Liu, C. Xu, W. H. Knox, J. K. Chandalia, B. J. Eggleton, S. G. Kosinski, and R. S. Windeler, “Soliton selffrequency shift in a short tapered air-silica microstructure fiber,” Opt. Lett. 26(6), 358–360 (2001). 2. A. C. Judge, O. Bang, B. J. Eggleton, B. T. Kuhlmey, E. C. Magi, R. Pant, and C. M. de Sterke, “Optimization of the soliton self-frequency shift in a tapered photonic crystal fiber,” J. Opt. Soc. Am. B 26(11), 2064–2071 (2009). 3. D. V. Skryabin, F. Luan, J. C. Knight, and P. S. Russell, “Soliton self-frequency shift cancellation in photonic crystal fibers,” Science 301(5640), 1705–1708 (2003). 4. K. Wang and C. Xu, “Tunable high-energy soliton pulse generation from a large-mode-area fiber and its application to third harmonic generation microscopy,” Appl. Phys. Lett. 99(7), 071112 (2011). 5. Y. Tang, L. G. Wright, K. Charan, T. Wang, C. Xu, and F. W. Wise, “Generation of intense 100 fs solitons tunable from 2 to 4.3μm in fluoride fiber,” Optica 3(9), 948–951 (2016). 6. J. P. Gordon, “Theory of the soliton self-frequency shift,” Opt. Lett. 11(10), 662–664 (1986). 7. F. M. Mitschke and L. F. Mollenauer, “Discovery of the soliton self-frequency shift,” Opt. Lett. 11(10), 659–661 (1986). 8. B. Zysset, P. Beaud, and W. Hodel, “Generation of optical solitons in the wavelength region 1.37–1.49 μm,” Appl. Phys. Lett. 50(16), 1027–1029 (1987). 9. N. Nishizawa and T. Goto, “Compact system of wavelength-tunable femtosecond soliton pulse generation using optical fibers,” IEEE Photonics Technol. Lett. 11(3), 325–327 (1999). 10. A. B. Fedotov, A. A. Voronin, I. V. Fedotov, A. A. Ivanov, and A. M. Zheltikov, “Powerful wavelength-tunable ultrashort solitons in a solid-core photonic-crystal fiber,” Opt. Lett. 34(6), 851–853 (2009). 11. M. C. Chan, S. H. Chia, T. M. Liu, T. H. Tsai, M. C. Ho, A. A. Ivanov, A. M. Zheltikov, J. Y. Liu, H. L. Liu, and C. K. Sun, “1.2to 2.2-μm tunable Raman soliton source based on a Cr:Forsterite laser and a photoniccrystal fiber,” IEEE Photonics Technol. Lett. 20(11), 900–902 (2008). 12. M. Y. Koptev, E. A. Anashkina, A. V. Andrianov, V. V. Dorofeev, A. F. Kosolapov, S. V. Muravyev, and A. V. Kim, “Widely tunable mid-infrared fiber laser source based on soliton self-frequency shift in microstructured tellurite fiber,” Opt. Lett. 40(17), 4094–4097 (2015). 13. E. A. Anashkina, A. V. Andrianov, M. Yu Koptev, S. V. Muravyev, and A. V. Kim, “Generating femtosecond optical pulses tunable from 2 to 3 μm with a silica-based all-fiber laser system,” Opt. Lett. 39(10), 2963–2966 (2014). Vol. 26, No. 15 | 23 Jul 2018 | OPTICS EXPRESS 19637


Introduction
Soliton sources based on soliton self-frequency shift (SSFS) in optical fibers with anomalous dispersion provide superb pulse quality, femtosecond pulse duration and wide spectral tunability [1][2][3][4][5][6][7]. When propagating in an optical fiber, soliton pulses are subject to stimulated Raman scattering (SRS). The longer-wavelength portion can experience Raman amplification at the expense of the optical power in the shorter-wavelength part of the spectrum. At wavelength regimes > 1300 nm, SSFS in standard single-mode fibers and index-guided photonic crystal fibers has been demonstrated [8][9][10][11]. Recently, soliton wavelengths up to 2.65 μm, 3 μm, 3.42 μm and 4.3 μm were realized in microstructured tellurite fibers [12], germanium-doped silica fibers [13], chalcogenide fibers [14] and fluoride fibers [5], respectively. However, these solitons suffer from low pulse energy. One way to increase the soliton pulse energy is to increase the effective mode area Aeff in the fiber, because the pulse energy of the Raman-shifted soliton in a fiber is proportional to the effective mode area of the fiber [15]. In previous reports, wavelength-tunable solitons from 1580 nm to 2130 nm with pulse energy up to 44.8 nJ in large-mode-area (LMA) fibers with A eff = 530 μm 2 were demonstrated [4,15]. However, solitons at longer wavelengths in LMA fibers with high pulse energy have not been demonstrated. In addition, current LMA fibers or photonic crystal rods are all based on silica glass. In this paper, we explore the long wavelength limit of silicabased fibers using SSFS. One of our main motivations and practical applications of pushing the long wavelength limit of SSFS in LMA silica fibers is to create a convenient femtosecond source for long-wavelength multiphoton microscopy (MPM) at wavelength range beyond that of a mode-locked Ti:Sapphire laser.
MPM enables non-invasive, three-dimensional imaging deep within scattering biological tissue. It has been a valuable tool for biological and medical imaging applications [15][16][17][18][19][20][21][22][23][24]. Soliton sources have been successfully applied in MPM. A photonic-crystal-fiber-based soliton source that accommodates the requirements of multi-modality biophotonic imaging was demonstrated [25]. Using a fiber-based energetic soliton source, in vivo hippocampal three-photon excited fluorescence imaging to a depth of 1.35 mm in an adult mouse brain has been demonstrated at the 1700 nm spectral window [22]. Using second harmonic generation, excitation sources have been demonstrated in a variety of fibers, e.g., an LMA fiber for twophoton imaging at 775-950 nm [26], a highly nonlinear germanosilicate fiber for coherent anti-Stokes Raman scattering imaging at 850-1100 nm [27], and an LMA fiber for harmonic generation microscopy at 1150 nm [28]. Recently, the long-wavelength spectral window beyond 1060 nm has attracted a lot of attention in MPM because of the reduced tissue scattering at long wavelengths [29], and significant effort has been made to generate energetic pulses at 1060 nm to 1700 nm in fiber-based systems [30][31][32], including higher-order-mode fibers [33][34][35][36][37]. However, energetic (e.g., > 10 nJ pulse energy), wavelength-tunable soliton sources at the spectral window between 1060 nm and 1300 nm, where high-energy pulses are difficult to generate by a mode-locked laser (e.g., the Ti:Sapphire laser), have not yet been demonstrated. The generation of energetic femtosecond pulses between 1060 nm and 1300 nm today typically requires an optical parametric oscillator (OPO) or an optical parametric amplifier (OPA).
In this paper, we explore the long wavelength limit of SSFS in silica-based fibers, and demonstrate SSFS in a silica LMA fiber with an Aeff ~760 μm 2 . The generated solitons are wavelength tunable from 1580 nm to 2520 nm, with a pulse energy up to 73 nJ and a pulse width ~100 fs. We perform second harmonic generation (SHG) of the solitons (data shown from 950 to 1260 nm), with pulse energy up to 21 nJ. The frequency-doubled soliton source has large wavelength tuning range, short pulse duration, and high pulse energy. Finally, in vivo two-photon fluorescence imaging of vasculature and neurons in a mouse brain is shown as one of the potential applications by using 1150 nm as the excitation wavelength.

Experimental results
The experimental setup is shown in Fig. 1. The first pump source in our experiment is a fiberbased femtosecond laser (FLCPA-02C, Calmar), delivering 412-fs and linearly polarized pulses with 0.66-MHz repetition rate at 1550 nm. SSFS is performed in a 2-m LMA fiber (LMA-PM-40-FUD, NKT Photonics, which is similar to DC-200/40-PZ-Si, NKT Photonics) with an A eff of ~760 ± 100 μm 2 at 1550 nm. The coupling efficiency is measured to be ~64% (measured at the low input power of ~36 mW, using a C-coated lens with a focal length of 50 mm). With the increase of the input pump power, solitons form and continuously shift from 1580 nm to 2520 nm. To separate the most red-shifted soliton from the residual input and other solitons, long-pass or band-pass filters are used. Pulses at the wavelength between 950 nm and 1260 nm are then obtained by SHG of the soliton source, using a frequency-doubling crystal BiBO (θ = 8°, φ = 0°, L = 1.5 mm, Castech). A short-pass filter is then used to separate the SHG from the solitons. To demonstrate the potential of this source for longwavelength two-photon microscopy, pulses at 1150 nm are sent into a microscopy and used for in vivo imaging of vasculature and neurons in a mouse brain.

Soliton generation
Numerical results of launching 1 μJ and 412 fs Gaussian-shaped pulses at 1550 nm into the fundamental mode of a 2-m-long LMA fiber are shown in Fig. 2. The numerical simulation is performed in the same manner as described in [15]. The dispersion of the LMA fiber is assumed to be identical to silica glass, since the waveguide dispersion for the LMA fiber is negligible. The linear propagation loss α is set to be 0.1 dB/m. Since the output pulse of the pump source has some side lobes, which don't contribute to the SSFS process, the pulse energy used in the numerical simulation is scaled down to 65% of the pulse energy used in the experiments to better match the experimental data. The spectral and temporal evolution of the pulses along the fiber exhibit the main feature of SSFS, as shown in Fig. 2. The input pulse excites a higher-order soliton and experiences soliton-effect compression. On reaching its minimum pulse duration, the pulse breaks up into a series of fundamental solitons, under the perturbation of higher-order dispersion and Raman scattering. These solitons experience dispersion and frequency shift toward lower frequencies (intrapulse SRS), and eventually separate spectrally and temporally. A comparison of numerical simulations and experimental results is shown in Fig. 3. We first characterize the solitons. Experimental results and numerical simulations for a number of input pulse energies are shown in Figs. 3(a) and 3(b), respectively. The soliton continuously shifts to longer wavelengths as input pulse energy increases and thus, the wavelength of the most red-shifted soliton can be tuned by adjusting the input pulse energy, as shown in Fig. 3(a). Note that the input pulse energies given in Fig. 3 are already scaled by the coupling efficiency (64%). The center wavelengths of the most red-shifted solitons as a function of the input pulse energy into the LMA fiber are shown in Fig. 4(a). The longest center wavelength is 2520 nm. The experimental results (blue solid line with squares) and numerical simulations (red dashed line with stars) are both shown in Fig. 4(a). The pulse energy of the most red-shifted soliton is measured by using long-pass (LP, 1630 nm, 1720 nm, 1950 nm, Semrock) or band-pass (BP, 2000-2250 nm, 2250-2500 nm, 2500-2750 nm, Thorlabs) filters and a power meter, taking into account of the transmission of the filters and the collimation lens. The measured soliton energies as a function of the soliton wavelength are shown in Fig. 4(b). The data points represented by the solid squares are directly measured using filters and power meters. Because of the lack of LP or BP filters, for the data points represented by the solid circles, we first measured the total power of the entire filtered spectra using a power meter, and then calculated the power of the most red-shifted soliton using the measured spectra (i.e., by calculating the power underneath the soliton spectrum as a fraction of the power of the entire filtered spectrum). The general trend that the soliton pulse energy increases as its center wavelength increases is observed for wavelength below 2400 nm. The pulse energy has a small dip at ~2200 nm and a rapid drop for > 2400 nm, which can be explained by absorption by silica glass and the OH-vibrational overtones of water content in the fiber (see next section). For the most red-shifted solitons, the bandwidths as a function of the soliton wavelength are shown in Fig. 4(c), and the measured second-order interferometric autocorrelation traces are shown in Fig. 3(c). The pulse widths as a function of the soliton wavelength are shown in Fig. 4(d). In particular, there is a good agreement between the numerical simulations (red stars) and the transform-limited (TL) pulse widths (blue solid circles) calculated from the measured bandwidths [ Fig. 4(c)], assuming a sech 2 intensity profile. The measured pulse widths (black solid lines with squares) are somewhat higher than the TL pulse widths. This pulse broadening may be due to dispersion of the filters and the collimation lens and the decreasing bandwidth at long wavelengths (see next section).

The long wavelength limit of SSFS in silica-based fibers
To explore the long wavelength limit of SSFS in silica-based fibers, we conduct further analysis and experiments. To obtain cleaner solitons, we employed an optical parametric amplifier (Monaco-Opera F, Coherent) as the pump source, with shorter pulse width of 70 fs at the wavelength of 1550 nm. The coupling efficiency is measured as 61%. The rest of the experimental configuration is the same as before. The LMA fiber is kept straight to minimize bending loss during these experiments. The measured spectra of SSFS after the LMA fiber are shown in Fig. 5. The center wavelengths of the most red-shifted solitons as a function of the input pulse energy into the LMA fiber are shown in Fig. 6(a). The longest center wavelength is again 2520 nm. Compared to Fig. 4(a), less input pulse energy is required for shifting to the same wavelength due to the shorter and cleaner input pulse. Since there is no overlap in the spectrum between the most red-shifted soliton and the secondary solitons, we can clearly observe that the bandwidth of the most red-shifted soliton starts to narrow after 2350 nm, as shown in Fig. 6(b). This is accompanied by the rapid drop in pulse energy and the increase in pulse width from 2400 nm to 2520 nm, as shown in Figs. 4(b) and 4(d).  We performed numerical simulations with increasing loss at the long wavelength, to investigate the impact of light absorption by silica and residue water content on SSFS. The loss due to silica absorption starts around 2300 nm and increases for longer wavelengths, whereas the loss due to water absorption starts around 1800 nm and increases for longer wavelengths as well [38]. The first loss curve we used in the numerical simulation is similar to the loss of bulk dry silica [38], as shown in Fig. 6(c). The bandwidth and pulse width of the simulation results are shown in Figs. 6(b) and 6(d), respectively, which significantly underestimates the decrease in pulse bandwidth when compared to the experimental results. This result shows that absorption by silica glass alone cannot account for the long wavelength limit of SSFS, and other loss mechanisms, such as OH-absorption, may also contribute. The evidence of OH-absorption is manifested by the experimentally measured soliton pulse energy and bandwidth. Both show a dip between 2200 nm and 2400 nm [ Figs. 4(b) and 6(b)], which coincides with the OH-absorption peak around 2200 nm [38]. The second loss curve we used in our numerical simulation also starts around 2300 nm but increases faster than the first loss curve, as shown in Fig. 6(c). The bandwidth and pulse width of the simulation results are shown in Figs. 6(b) and 6(d), respectively. The agreement between experimental results and numerical simulation results with the second loss curve is much improved. Therefore, SSFS in the LMA fiber is most likely limited by absorption of both silica and the residue water content. However, even without water absorption, the longest wavelength of SSFS in silica-based fibers will still be limited by the absorption of silica glass to approximately 2400 nm in order not to suffer significant loss of soliton energy and bandwidth [see simulation results in Figs. 6(b) and 6(d)].
These conclusions on the long wavelength limit are further supported by SSFS experiments performed using a low-OH single mode optical fiber (experimental fiber from Corning, the parameters are similar to those of SMF-28 except for the low-OH feature. The attenuation due to OH at 1380 nm is about 0.026 dB/km). The measured spectra of SSFS after the low-OH SMF are shown in Fig. 7(a). The center wavelengths of the most red-shifted solitons as a function of the input pulse energy are shown in Fig. 7(b). The longest center wavelength is again ~2500 nm. The bandwidth of the most red-shifted solitons is shown in Fig. 7(c), which again starts to decrease around 2300 nm. Interestingly, the soliton bandwidth doesn't have the dip between 2200 nm and 2400 nm as shown in Figs. 4(b) and 6(b), confirming that absorption due to the water content in the low-OH fiber is negligible.

Second harmonic generation
We then characterize the SHG pulses generated from the fiber-based soliton source (Fig. 4). The measured spectra of the SHG pulses are shown in Fig. 8(a). The bandwidths of the spectra centered from 950 nm to 1260 nm vary from 15 nm to 19 nm. The measured secondorder interferometric autocorrelation traces are shown in Fig. 8(b). For SHG pulses from 950 nm to 1260 nm, the deconvolved pulse widths [ Fig. 8(b)] increase from 77 fs to 180 fs (using the deconvolution factor of 1.54 for sech 2 intensity profile), following a similar trend as the soliton pulse width. These results are longer than the calculated transform-limited pulse durations assuming sech 2 profiles (from 52 fs to 106 fs), which is likely due to group velocity mismatch [26]. The stability of the SHG pulses at 1150 nm is characterized by measuring the two-photon current generated in a GaAsP photodiode when excited by the source. The measured root-mean-square (RMS) fluctuation of the two-photon signal is 3.274% (from 0.003 kHz to 50 kHz). The pulse energy of the SHG pulses at 1150 nm is measured to be 21 nJ, corresponding to a conversion efficiency of 44.6% for the fundamental wavelength at 2300 nm. Higher conversion efficiency has been demonstrated using the same crystal for SHG at other wavelengths [26,39]. Assuming similar SHG efficiency, the pulse energy for 950 nm to 1260 nm will be in the range of 13 nJ to 32 nJ, with the peak pulse energy at ~1200 nm. To achieve the highest conversion efficiency, however, the orientation of the SHG crystal should be adjusted to the phase-matching angle for each input wavelength. This may produce slightly different beam pointing. The beam pointing problem can be solved by using a periodically poled lithium niobate (PPLN) for SHG.

In vivo two-photon excited fluorescence imaging
One of our main motivations and practical applications of pushing the long wavelength limit of SSFS in LMA silica fibers is to create a convenient femtosecond source for longwavelength MPM at wavelength range beyond that of a mode-locked Ti:Sapphire laser. To demonstrate the applications of our source (using the fiber-based femtosecond laser and SSFS), we imaged vasculature and neurons in a mouse brain in vivo. The mouse is 3-month old. The overlying skull was removed by performing a craniotomy before imaging. The vasculature was labelled by a retro-orbital injection of dextran-coupled Texas-Red dye (Invitrogen). The excitation wavelength is 1150 nm. Mouse brain vasculature down to 420 μm below the surface of the brain [Figs. 9(a)-9(d)] can be obtained with a maximum optical power on the brain surface of 2 mW (i.e., 3 nJ). We also imaged red fluorescent protein (RFP) labelled neurons in a B6.Cg-Tg(Thy-Brainbow1.0)HLich/J mouse, with a depth of 264 μm [Figs. 9(e)-9(h)].  Figure 9 showed that the fiber-based femtosecond source and SSFS can provide energetic femtosecond pulses for two-photon imaging at wavelengths currently only reachable by OPO or OPA. While the pulse energy is adequate, the average power, and therefore the repetition rate, demonstrated here is too low for optimum two-photon imaging [39]. The average power is limited by the fiber-based femtosecond laser used in this work. However, it is relatively straightforward to scale up the average power (i.e., scaling up the repetition rate) of fiberbased femtosecond laser at 1550 nm while maintaining the same pulse energy [39]. For twophoton deep tissue imaging in the mouse brain, the optimum repetition rate for many common fluorophores is between 1 and 10 MHz [39]. This is a result of the trade-off between the average power at the sample surface and the peak power at the focus. The low repetition rate for deep imaging potentially limits the ability for fast imaging of dynamic phenomena since a minimum of one excitation pulse per pixel is required. Our demonstration shows that SSFS in a silica-based LMA fiber, together with SHG, is a possible approach for obtaining energetic (> 20 nJ), wavelength-tunable femtosecond pulses at > 1060 nm for long-wavelength twophoton microscopy. On the other hand, material absorption in a silica-based fiber appears to limit the longest wavelength of such an approach to ~1260 nm, which only partially overlaps with the long wavelength window of 1300 nm for deep tissue imaging (~1200 nm to 1350 nm) [22]. Further optimizations of the fiber design, e.g., perhaps with heavy Germanium doping, will be needed to push the long wavelength limit of this approach and expand the spectral coverage.

Conclusions
In summary, we explore the long wavelength limit of SSFS in silica-based fiber experimentally and using numerical simulation. We found that the longest wavelength soliton generated by SSFS is approximately 2500 nm. The soliton loses its energy rapidly at wavelength beyond 2400 nm due to absorption by silica and residue water content. We demonstrate a 1580-2520 nm wavelength-tunable, high-pulse energy soliton source. Efficient SHG from 950 to 1260 nm is obtained, and its potential application for in vivo two-photon fluorescence microscopy is demonstrated. By increasing the repetition rate and further improving the pulse quality (e.g., reducing the pulse width and the pedestal) of the energetic fiber-based femtosecond laser, the demonstrated method has the potential to provide a low-cost, wavelength-tunable, energetic femtosecond source for long-wavelength two-photon microscopy at the spectral range currently only reachable with OPOs or OPAs.