3.1. Waveguide Meaurement
The photograph of the experimental setup is shown in
Figure 3a. On the input side, light from a tunable laser (1500-1600 nm, T100S-HP, EXFO) and a 2004 nm CW laser (Ultraline 2000, NPI Lasers) can be manually switched to couple into the waveguide. The fiber-chip-fiber alignment is facilitated by a set of fine positioning stages. The 3-axis piezoelectric controller can move the fiber at a precision of ±10 nm along X, Y, and Z directions. A pair of lensed fiber (TSMJ-X-1550-9/125-0.25-7-2.5-14-2-AR, OZ Optics) is adopted for waveguide coupling for better efficiency, as shown in
Figure 3b and 3c. On the output side, an optical component tester (CT440-PDL, EXFO) is used for collection and analysis of the transmitted light from 1500 nm to 1600 nm. For the 2004-nm laser, the output light is coupled to a thermopile optical power meter (3A, Ophir) for characterization.
A standard cut-back method measurement is performed using three sets of identical Ge nanostrip waveguides with different lengths, i.e., 0.5 cm, 1 cm, and 1.5 cm, respectively. First, the Ge waveguide chip is placed between the lensed fiber pair and finely aligned. The power from the tunable laser is set to 10 mW. The fiber polarization controller is finely adjusted to make the transmitted power reach its minimal value at the detector. Since the Ge waveguide can transmit the TM mode with significantly lower loss, the input polarization state giving the minimum transmitted power is considered to be the TE polarization. After scanning the transmitted power under the current polarization state from 1500 nm to 1600 nm, the laser system automatically rotates the polarization direction by 90°, i.e., to the TM polarization, and repeat the wavelength scan. For the second measurement, the input is switched to the 2004-nm laser (also set to 10-mW output power), and the output fiber is connected to the thermopile power meter. The polarization adjustment is repeated, and the transmission values are obtained.
The fiber-to-fiber transmission without the waveguide chip in the middle serves as a reference to rule out the system influence. Linear fitting of the cut-back measurement gives a clear indication of the waveguide propagation loss (slope) and the coupling loss with the fiber (intercept). The results for the 1500-1600 nm range are plotted in
Figure 4a but only for the TM mode. For the TE polarization, we have only detected light power at noise level (-70 dBm) for the 0.5-cm-long waveguide and the loss is expected to be well over 100 dB/cm. As indicated by the grey curve in
Figure 4a, the waveguide propagation loss gradually decreases from 20.56 dB/cm at 1500 nm to 11.75 dB/cm at 1600 nm. At 2004 nm, the transmitted power of TE-polarized light input still stays below the detection limit (-30 dBm) of the chosen power meter, while the linear fitting reveals a propagation loss of 4.89 dB/cm for the TM mode, as shown in
Figure 4b. The TE polarization undergoes a strong absorption and cannot be decisively measured even at 2004 nm, while the TM polarization shows an increasing degree of transparency from 1500 nm to 1600 nm and further to 2004 nm. The Ge nanostrip waveguide with a cross-section of 4 μm × ~40 nm is only demonstrated. The ~40-nm thickness is chosen since the waveguide principle is similar to high-aspect-ratio dielectric waveguide, e.g., silicon nitride waveguide in glass [
42]. Thicker Ge layer will increase the propagation loss. Numerical simulations reveal that smaller Ge nanostrip widths/thicknesses can further weaken the mode confinement and effectively reduce the propagation loss of the TM fundamental mode. We project the propagation loss will be lower than 2 dB/cm including the absorption in the polymer cladding at 2 μm wavelength.
Next, the ellipsometer measurement is performed on the dummy sample to obtain the optical parameters (
n and
κ). Considering the random atom orientations within the amorphous Ge film in the absence of long-range-order, and the fact that a ~40-nm-thick layer already stacks ~150 Ge atoms, the in-plane polarization (ordinary light or “o” light, referring to the TE polarization in the waveguide) and the vertical polarization (extraordinary light or “e” light, referring to the TM polarization in the waveguide) “sees” nearly the same optical constants. Therefore, an isotropic material model is adopted for the ellipsometry. The original data of the test sample is measured by 1650 nm. The refractive index
n and the extinction coefficient
κ is well fitted and plotted as the solid curves in
Figure 5. One Tauc-Lorentz oscillator and two Lorentz oscillators are applied to reshape the imaginary part of the dielectric function, i.e., the extinction coefficient
κ, which ensures a low MSE while avoiding overfitting by excessive oscillators. Both optical constants exhibit similar trends to the previously published reference [
44]. Since the light source of the ellipsometer cuts off at 1650 nm, the optical constants at longer wavelengths need to be extended by mathematical fitting, and the Cauchy model is considered a suitable candidate for such treatment [
45].
The Cauchy model gives
where
λb is usually a value near the shortest wavelength of the fitting region. Data in the 1400-1650 nm wavelength range in the original
n-
κ model is selected as the extension base. The dashed curves in
Figure 5 plots the extension results. The mean square error (MSE) stays on the order of 1×10
-6, indicating a reliable fitting for the model extension.
3.2. n-κ Model Correction
In ellipsometry, the imaginary part of the complex dielectric function is usually fitted by a combination of independent oscillators. For the measurement of thin films with strong optical absorption, the absorbing regions can be described by several symmetric/asymmetric oscillators, while their “tails” stack up in the near-transparent region. The consequence is that κ can only accumulate with the increasing number of oscillators, resulting in an overestimation of the extinction coefficients near the cut-off wavelength. Moreover, this accumulation is hard to tune down by modeling, unless a simpler oscillator combination is adopted, but this often results in a higher MSE. On the other hand, accurate optical parameters are crucial for the precise component design, e.g., a saturable absorber in a laser cavity. The existing Ge n-κ model, especially the extinction coefficient κ, needs confirmation and possibly a correction to match the experimental results.
In the first step, the real part (
n) of the complex index is investigated. Since
n is directly related to the eigenmode field size, deviation of the model can be revealed by comparing the measured mode field area (MFA) with the simulated value using a commercial eigenmode solver. As shown in
Figure 6a, an imaging system is established with a diaphragm, a set of lenses, a polarizer, and a camera. This setup is similar to a commercial beam profiler, which has been extensively used to examine the laser beam waist profile [
46,
47]. The tunable laser from 1500 nm to 1600 nm is coupled to the chip, while the beam spot at the output facet is magnified on the camera through the imaging system. The polarizer direction is first adjusted as parallel to the in-plane polarization direction of the Ge waveguide. Then the polarization controller is finely tuned to make the beam spot completely invisible on camera, reaching the TE mode. Subsequently, the polarizer is rotated by 90°, and the imaged mode field is considered as the pure TM mode. The first row in
Figure 6b plots the captured spots from the waveguide output facet at 1500 nm, 1550 nm and 1600 nm wavelength, respectively. The mode field of a standard single-mode fiber at 1550 nm wavelength is also captured as reference. The second row plots the simulated mode fields at corresponding wavelengths using the
n-
κ model presented in
Figure 5 (for Ge waveguide modes). The MFA of the captured waveguide mode can be calculated as
where the reference MFA
fiber mode is ~85 μm
2 at 1550 nm from the fiber data sheet, and effective pixel stands for the count of pixels in which the captured light power is no smaller than 1/
e2 of the peak value.
Figure 6c gives a comparison of the measured/simulated MFAs at 1500 nm, 1550 nm, and 1600 nm wavelengths, respectively. Each measurement is repeated multiple times and the average value is taken. The deviations between measured/simulated MFAs stay within 5%.
In the second step, the imaginary part
κ is corrected with the previously obtained transmission characteristics of the Ge waveguides. As listed in the first row of
Table 1, six wavelengths with a step size of 20 nm are selected as samples. As presented in the second row, corresponding Ge core absorptions at these wavelengths are extracted from the grey curve (Ge waveguide loss vs.
λ) in
Figure 4a. According to the material datasheet, the polymer cladding exhibits an absorption loss of ~0.35 dB/cm at 1550 nm wavelength. This loss is slightly wavelength dependent from 1500 nm to 1600 nm but remains about two orders of the magnitude smaller than the Ge material absorption. Therefore, we neglect the polymer absorption in the calculation. Since no etching is involved in the definition of the Ge nanostrip, the waveguide scattering loss can also be neglected. Therefore, the waveguide loss arises mainly from the Ge core absorption.
Next, we establish the Ge waveguide model in the commercial eigenmode solver and scan
κ in the Ge material definition at all six sample wavelengths to make the simulated propagation loss as close as possible to the measured value. The given
n of the Ge core for these simulations is listed in the third row of
Table 1, which is directly obtained from the ellipsometer measurement in
Figure 5. Finally, the optimized
κ is listed in the fourth row of
Table 1 and considered to reveal the real
κ at the corresponding wavelength. The relation of
κ vs.
λ from 1500 nm wavelength is fitted by Cauchy mode in Eq. (2). The corrected
κ is plotted as the dashed curve in
Figure 7, which exhibits a significant offset from the original model. The MSE of the Cauchy fitting is on the order of 1×10
-8.
The corrected
n-
κ model is further examined at 2004 nm wavelength. First, a series of polymer channel waveguides with 1.47-1.45 core-cladding index contrast are measured to estimate the polymer absorption at 2004 nm. The cross-sectional design of the waveguide core is 3.5 μm × 3.5 μm. As shown in
Figure 8a, the cut-back fitted propagation loss of the TE and TM mode is 3.27 dB/cm and 3.16 dB/cm, respectively. The sidewall scattering loss of the polymer is estimated to be 0.14 dB/cm at 2004 nm wavelength by Rayleigh scattering principles. Due to similar material compositions, the cladding polymer (ZPU450) and the core polymer (ZPU470) share also similar optical absorption. Therefore, the polymer absorption is determined to be 3.02 dB/cm for the TM polarization at 2004 nm, as summarized in
Figure 8b.
For the Ge waveguide, the dominating field component
Ez of the TM fundamental mode is mainly concentrated in the polymer cladding. Thus, as plotted in
Figure 8b, the absorption for the polymer cladding is already obtained as 3.02 dB/cm, and the absorption of Ge core is calculated to be 1.87 dB/cm at the 2004 nm wavelength. Finally, the corrected material
n-
κ data is imported into the established Ge waveguide model in the commercial eigenmode solver. At 2004 nm wavelength, the solver gives a propagation loss of 1.71 dB/cm (without polymer absorption). The 0.16 dB/cm deviation from the measured value may be attributed to the slight thickness deviation of the Ge layer in the experiment. Nevertheless, the results above demonstrate the accuracy and credibility of the correction method for the optical parameters of Ge. This waveguide-assisted characterization method can also be extended to other thin-film materials in the near-transparent region with evident bandgaps.