Circular cladding waveguides in Pr:YAG fabricated by femtosecond laser inscription: Raman, luminescence properties and guiding performance

We report on the fabrication of circular cladding waveguides with cross-section diameters of 60−120 μm in Pr:YAG crystal by applying femtosecond laser inscription. The fabricated waveguides present 2D guidance on the cross-section and fairly low propagation losses. Multiple high-order guiding modes are observed in waveguides with different diameters. Corresponding simulation results reveal the origin of a specific kind of guiding modes. Confocal micro-Raman (μ-Raman) experiments demonstrate the modification effects in femtosecond laser affected areas and ascertain the refractive index induced guiding mechanism. In addition, luminescence emission properties of Pr 3+ ions at waveguide volume are well preserved during the femtosecond laser inscription process, which may result in a potential high-power visible wave-guide laser.


Introduction
In the past decades, femtosecond laser inscription (FLI) as a new method for micromachining 1−2 , has emerged to be an effective technology for optical waveguide fabrication in various transparent dielectric materials. And in recent years, various kinds of waveguides and functional devices have been realized in some kinds of crystals, such as Yb, Na:CaF 2 3 , diamond 4−5 , Nd:GdVO 4 6 , dye polymer 7 , Gallium Lanthanum Sulphide glass 8 and LiNbO 3 9−11 . Nonlinear energy transfer that occurs during the FLI process simultaneously produces localized moderation at the focal spot. The increase, decrease of refractive index or even material plasmatization can occur depending on the utilized femtosecond laser 12 . In addition, waveguide structure with a tubular cladding geometry (unmodified guiding region is surrounded by laser induced low-refractive-index damage tracks) gains growing attention because the guiding area of such a waveguide can still maintain the original material properties after the FLI process 13−19 . Furthermore, tubular cladding waveguide shows higher stability and consistency than the stress induced dual-line and the single-line waveguides, as well as supports orthogonal polarizations in the meantime.
Praseodymium ions (Pr 3+ ) present extremely complex energy level and relatively low photon energy, which make them applicable for laser transitions with different excitation wavelengths ranging from visible to infrared.
During the past few years, laser generations have been reported based on different Pr-doped fluoride materials: waveguide laser in Pr:YLF 20  In this paper, we exhibit the fabrication of circular cladding waveguides with low propagation losses in Pr:YAG crystal by using FLI. By investigating the guiding mode profiles and confocal μ-Raman properties, localized modification effects caused by the femtosecond laser and the guiding mechanism are revealed. In addition, luminescence and guiding properties of the guiding area are also studied.

Waveguide fabrication
The Pr:YAG crystal sample (praseodymium doped yttrium aluminum garnet, Pr:Y 3 Al 5 O 12 , doped by 0.5 at. % Pr 3+ ions, obtained from Atom Optics Co. Ltd, Shanghai, China) used in this work is cut into dimensions of 10 mm × 6 mm × 2 mm with all facets optically polished. During the waveguide fabrication process, an amplified Ti:sapphire laser system (Astrella, Coherent Inc., USA) is used as laser source, which generates linearly polarized pulses (with repetition rate of 1 kHz, central wavelength of 800 nm, pulse width of 35 fs, and single pulse energy of 7 mJ). Watt Pilot motorized attenuator and a neutral density filter are utilized to adjust the laser power precisely, and specifically, single pulse energy of 0.23 μJ is used in this work. The laser beam is focused by a microscope objective (40 ×, NA = 0.6) into the substrate beneath a 10 mm × 6 mm facet (the upper facet). The substrate is placed on a 6-axes precision motorized stage and the scanning speed is set to be 0.8 mm/s. Under the above experiment conditions, one damage track with vertical width of 15 μm and horizontal width of 2.6 μm can be inscribed via a single scan. Arrays of parallel damage tracks are inscribed following the designed geometries (the horizontal separation between ad-jacent tracks is 3 μm) so as to form desired cladding waveguides with different diameters, as shown in Fig. 1. All the waveguides are fabricated along the y-axis of the substrate (i.e., the waveguide lengths are all 10 mm).
Confocal μ-Raman measurement and luminescence characterization Raman properties of the femtosecond laser induced damage tracks and the formed guiding core area are explored by a confocal μ-Raman spectroscopy system (XperRam200, Nanobase, South Korea). A continuous laser beam (532 nm, 30 mW) is focused onto the endface by a microscope objective (40 ×, NA = 0.75). Part of the beam is affected by the Raman scattering effect of the sample and is then collected and transferred to a spectrometer system. With the help of a galvo scanning system, we can accurately gain the Raman properties of every region over the whole waveguide cross-section.
As for the luminescence characterization experiment process, a continuous laser beam at 400 nm is focused onto the end-face of the sample. And the emergent light from another end-face is collected and transferred to a spectrometer system which can provide the luminescence information of irradiated area in real time.

Waveguide characterization and guiding mode simulation
The near field intensity distributions of guiding modes are measured by an end-face coupling system. This system uses a 632.8 nm He-Ne laser as light source. A Glan-Taylor prism and a 1/2 waveplate at 632.8 nm are utilized to generate linear-polarized incident laser beam (i.e. TE and TM, TE: the electric field is parallel to the x-axis; TM: parallel to the z-axis). The linearly polarized beam is then coupled into one waveguide by a microscope objective (40 ×, NA = 0.65) and the transmission light is collected by a long working distance microscope objective (50 ×, NA = 0.55) from another end-face. The profile is subsequently imaged via a CCD camera (WCD-UCD12-1310, DataRay, USA).
The guiding mode simulation process is accomplished by the finite element method (FEM). The geometry structures used for simulation, as shown in Fig. 1(b), are determined by referring to the micrographs and 2D μ-Raman mapping images, as shown in Fig. 1(c) and Fig. 3, respectively. Since the mechanical stress produced by the physical lattice expansion at damage tracks attenuates rapidly over distance and has negligible influence on the waveguide core area, and to simplify the refractive index distribution used in the simulation, we assume a step-index configuration with the refractive index difference of -0.04 at all damage track boundaries.

Results and discussion
Figures 1(a) and 1(b) exhibit the 3D and 2D geometry structures of waveguides, and the 2D one is also used in the guiding mode simulation process. Figure 1(c) shows the cross-section microscope images (optical transmission) of the circular cladding waveguides with different diameters: 120 μm, 100 μm, and 60 μm, respectively. Since our aim of this work is to build waveguides that fit the demand of luminescence and laser generation, we are making waveguides with relatively large diameters. But according to our former experience, cladding waveguides with minimum diameter of 10 μm can be fabricated in such kind of crystal. The damage tracks (of different fabricated waveguides) closest to the upper facet share a depth of 120 μm beneath the upper facet. The refractive index change in the track regions is estimated to be -0.04 by measuring the NA value of the waveguides. No evident damage in the core areas of the fabricated waveguides is oberserved.
Raman properties of one fabricated waveguide with diameter of 120 μm are shown in Fig. 2 and Fig. 3. Raman spectra accurately measured at different positions reveal the variation of lattice after FLI process. As we can see from Fig. 2, the Raman intensity at damage track region suffers a strong quenching without evident variation on the spectral distribution while the Raman spectra of non-processed bulk area and waveguide volume are indistinguishable.
In order to further investigate the Raman spectra variation over the whole waveguide cross-section, the 2D mapping is implemented and the corresponding results are presented in Fig. 3. Intensity, shift and full width at half maximum (FWHM) changes of characteristic peak (259 cm -1 ) are shown as different channels. And significant differences: Raman intensity quenching, blue shift (i.e., wavenumber decrease) and broadening of the emission peak can be observed at the damage track regions. These phenomena indicate the formation of lattice defects/damages and the elongation of bond-lengths (i.e., lattice parameter) due to the expansive stress. The partial dilatation also results in the refractive index decrease at these regions. Therefore, the refractive index induced guiding mechanism of such cladding waveguides can be ascertained. Otherwise, the waveguide volume manifests fairly similar Raman properties with respect to the nonprocessed bulk area, which demonstrates that the lattice structure in the guiding area is well preserved during the FLI process. A further analysis about the Raman spectra of damage tracks and the surrounding areas is introduced here. As mentioned above, significant differences are observed at the damage tracks with respect to the non-processed bulk area. Raman intensity quenching and peak broadening originates from the formation of lattice defects and damages; while for peak blue shift, the elongation of bond-lengths becomes the dominant factor. That is in good agreement with the subjective impressions, which is longer bond-lengths correspond to lower molecular vibration frequencies, closer molecular vibration energy levels, and lower Raman shift for specific peaks. But when we focus on the surrounding areas of damage tracks, a strong red shift (i.e., wavenumber increase) occurs without evident variation on Raman intensity or FWHM with respect to the non-processed bulk area, indicating a homogeneous lattice compression at these areas. In other words, regions modified by the femtosecond laser present lattice dilatation with refractive index decrease while the surrounding areas are squeezed by the modified regions, forming lattice compression with refractive index increase. This phenomenon is in agreement with several published works accomplished on different samples 1,5 . Figure 4 depicts the luminescence emission spectra of the guiding volume of one fabricated 120 μm waveguide and non-processed bulk area. Virtually identical emission spectra (both spectral distribution and intensity) indicate that the luminescence properties of Pr 3+ ions at waveguide volume have been well preserved during the FLI process.
All the fabricated cladding waveguides support both TE and TM polarizations (following the aforesaid definition). Since the diameters of the fabricated waveguide are relatively large (~100 μm) compared to the incident laser wavelength (632.8 nm), large amounts of high-order modes are supported to propagate in the waveguides. Because of that, different modes with different intensity

Fig. 4 | Luminescence emission spectra of waveguide volume (blue dashed line) and non-processed bulk area (red solid line).
Inset shows the output end-face image captured by a color CCD camera while a laser beam at 400 nm has been coupled into the waveguide through another end-face.
The guiding modes with a specific kind of intensity distribution of waveguides with different diameters and the corresponding simulation results are exhibited in Fig. 5. These guiding modes possess horizontal fringes with TE polarization, which is actually the characteristic of LP 1n modes (LP mode: linearly polarized mode). It's worth mentioning that these modes only exist when the incident laser beam is focused on the positions a little bit right above or below the waveguide center. As can be seen in this figure, the experiment and the simulation results match well.
Here, we introduce a further discussion about the origin of guiding modes with this specific kind of intensity distributions. In consideration of the similarity between the observed modes and the typical modes in optical fiber, an ideal step-index circularly-symmetric fiber model that has the same refractive index profile and waveguide core area with the 60 μm cladding waveguide model is established. Simulation result indicates that such fiber supports the propagation of TE 04 mode, as shown in Fig. 6(a). Circularly symmetric fiber supports the propagation of modes with the same electric field distribution but different poloidal rotation angles about the fiber axis, and these modes present nearly indistinguishable effective refractive indexes, which means, they are degenerate. In addition, simulation method that we used (using effective refractive indexes to distinguish modes, mode analysis, FEM) tends to show a combined mode as result instead of several degenerate modes if they have   Since these LP 14 modes are degenerate, they are more inclined to combine and form a circularly symmetric TE 04 mode in simulation, which means, a single LP 14 mode will not emerge under such condition. But if under non-circularly-symmetric condition, such as in elliptical fiber, previous degenerate modes with different polarizations will manifest similar electric-field distributions and slightly different effective refractive indexes, scilicet the removal of degeneracy 31 . To confirm the reliability of this thought, an elliptical fiber model with major axis of 60 μm and minor axis of 56 μm is established. And under such simulation condition, single LP 14 modes can be observed, as shown in Figs. 6(b) and 6(d), which is coincident with the experiment and simulation results of the fabricated 60 μm waveguide. This result demonstrates the non-circular symmetry of the fabricated waveguides, also suggests that with refined design, the fabric-ated waveguides can support the propagation of circularly symmetric modes. These conclusions are also applicable for the waveguides with different diameters.
The losses of the waveguides at 632.8 nm when the incident laser is differently polarized are given in Table 1. The total losses are directly measured by an optical power meter and the end-face coupling system used in the measurement of the intensity distributions of guiding modes while the coupling losses are roughly calculated by the BeamPROP module of Rsoft Photonics CAD Suite. The main simulation parameters are listed here: spot diameter of 10 μm; refractive index difference of -0.04. Therefore, the propagation losses of the waveguides can be directly given and the error induced during the measurements should be less than ±0.2 dB/cm. Since the beam diameter of 10 μm is relatively small compared to the waveguide diameters of 100 and 120 μm, simulation results reveal the same coupling loss value of the waveguides with diameters of 100 and 120 μm. The propagation losses result shows such circular cladding waveguides in Pr:YAG present fairly low propagation losses at 632.8 nm.

Conclusions
In conclusion, we have successfully fabricated circular cladding waveguides which support multi-mode propagation with orthogonal polarizations in Pr:YAG crystal by FLI. The fabricated waveguides show good guiding performance without any deterioration of the luminescence properties of the Pr 3+ ions in the waveguide volume. A specific kind of guiding mode profiles is analyzed in detail, and simulation results suggest the feasibility of mode profile control in such waveguides. By means of the Raman property analysis of damage tracks and the surrounding areas, we have ascertained the modification effects of femtosecond laser and the refractive index induced guiding mechanism. With excellent luminescence performance and fairly low propagation losses, the fabricated waveguides present potential for integrated luminescence and laser generation applied in advances photonics circuits.