Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Hydroxyl ion absorption in on-chip high-Q resonators

Open Access Open Access

Abstract

Thermal silica is a common dielectric used in all-silicon photonic circuits. Additionally, bound hydroxyl ions (Si-OH) can provide a significant component of optical loss in this material on account of the wet nature of the thermal oxidation process. A convenient way to quantify this loss relative to other mechanisms is through OH absorption at 1380 nm. Here, using ultra-high-quality factor (Q-factor) thermal-silica wedge microresonators, the OH absorption loss peak is measured and distinguished from the scattering loss baseline over a wavelength range from 680 nm to 1550 nm. Record-high on-chip resonator Q-factors are observed for near-visible and visible wavelengths, and the absorption limited Q-factor is as high as 8 billion in the telecom band. Hydroxyl ion content level around 2.4 ppm (weight) is inferred from both Q measurements and by secondary ion mass spectroscopy (SIMS) depth profiling.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

High-quality-factor ($Q$-factor) integrated resonators have become an essential component in nonlinear photonics. Most often, the guided light in these structures has a significant fractional overlap with silica (e.g., all-silica wedge resonators [1,2] and silica-clad ultra-low-loss silicon nitride waveguides [3,4]). It is therefore important to understand the loss limits imposed by the silica used in silicon photonic processing. In addition to interface scattering loss, optical absorption from bound hydroxyl ions (Si-OH) can be a significant component of loss [5], especially since thermal silica is prepared using a process involving water. Bound hydroxyl ions produce a well-known overtone absorption peak near 1380 nm [6] that is used here to measure OH absorption loss in ultrahigh-$Q$ thermal-silica wedge microresonators. Further comparison to scattering loss is made over a wavelength range from 680 nm to 1550 nm. Also, using cavity-enhanced photo-thermal spectroscopy [5] near the 1380-nm band, the OH ion content level is estimated to be 2.4 ppm (weight). This value also agrees with secondary ion mass spectroscopy (SIMS) depth profiling performed on the resonator material. Outside of the 1380-nm band, scattering loss is confirmed by measurement and modeling to be the dominant loss mechanism in the samples tested.

Thermally grown silica wedge whispering-gallery resonator devices were prepared as measurement samples, and featured 8-$\mathrm{\mu}$m thermal oxide thickness with resonator diameter 6.5 mm (10-GHz free spectral range, FSR). The fabrication steps are detailed in Ref. [2] and a photo micrograph of a typical resonator is the inset of Fig. 1. The thermal silica was grown from float-zone (low background doping level) silicon wafers using the wet oxidation method. As a final step, the devices were annealed for 18 hours at 1000$^{\circ }$C in N$_{2}$. The samples were stored and measured in a dry N$_2$ environment to minimize environmental impact on optical loss. Tapered fiber couplers [7,8] were used to couple probe light to the resonator samples.

 figure: Fig. 1.

Fig. 1. Optical loss measured from 680 nm to 1550 nm reveals the existence of the OH absorption peak near 1380 nm. The data are obtained by optical $Q$-factor measurement. The purple solid curve is the predicted loss from surface-roughness scattering (see text). The light blue solid curve is the estimated OH absorption loss in the 1380-nm band (corresponding to 2.4 ppm weight). Dashed contours correspond to different $Q$-factor values. Inset shows typical optical micrograph of the optical resonator used in the measurements.

Download Full Size | PDF

Figure 1 summarizes the optical loss measurements using six of these devices measured at wavelengths from 680 nm to 1550 nm. All resonators feature intrinsic $Q$-factor over 500 million at 1550 nm. The corresponding $Q$-factors are also given in the plot (dashed contours). Significantly, the devices have record-high $Q$ values for on-chip devices in the shorter wavelength bands: 300 million at 685 nm, 500 million at 780 nm, and 600 million at 1064 nm. As an aside, the loss per unit length $\alpha$ (dB/km) is determined from the optical $Q$-factor measurement using the relationship $\alpha =4343\times \frac {2\pi n_g}{\lambda Q_0}\; \textrm {(dB/km)}$, where $Q_0$ is intrinsic $Q$-factor in billions, $n_g$ is the group refractive index and $\lambda$ is the wavelength (nm). The intrinsic $Q$-factor $Q_0$ is obtained from the relationship $1/Q_0= 1/Q_L - 1/Q_c$ by measuring the resonance linewidth with Lorentzian line shape fitting to obtain the loaded $Q$-factor $Q_L$ followed by measurement of transmission depth to infer the coupling $Q$-factor $Q_c$.

A strong increase in loss near 1380 nm is apparent in all six samples, corresponding to the OH absorption band and reaching over 200 dB/km. This absorption quickly decreases for wavelengths above and below 1380 nm. The loss in other spectral regions is believed to be dominated by Rayleigh scattering. The wavelength dependence of this scattering within the resonator mode volume would scale approximately as $\lambda ^{-4}$, and does not fit the data. Modeling of surface scattering is described at the end of the Letter and provides better agreement with the wavelength dependence. This theoretical dependence is given by the purple curve (for the fundamental TE mode) where surface roughness variance (1.9 nm) and correlation length (350 nm) are assumed in the plot [9]. The origin of the increased loss at the shortest wavelengths measured is not known, but possibilities include absorption loss from metallic ion impurities and scattering from material density fluctuations. For example, SIMS data have shown that chromium is a residual contaminant in our processing, and Cr$^{3+}$ ions can contribute 1.6 dB/km/ppbw at 800 nm (peak at 625 nm) to absorption [10].

Cavity-enhanced photothermal spectroscopy (CEPS) [5] was used to further study the OH absorption loss. This method measures the microresonator resonance thermal triangle formation induced by the thermo-optical effect [14] to determine mode volume absorption. Details on this method are provided in Ref. [5]. Figure 2 summarizes the wavelength dependence of the measured absorption loss in both the 1380-nm band and at 1550 nm. The measured absorption near 1380 nm follows reasonably well the OH overtone line shape in silica as determined elsewhere [6] (blue curve superimposed in the plot). To fit the magnitude of the line shape function to the data, an OH content level of 2.4 ppm weight is used. Also, the data near 1550 nm indicate that by further reduction of scattering loss, the existing thermal silica can provide absorption-limited $Q$-factors as high as 8 billion. This value is 2$\times$ larger than measured for wedge resonators in a previous study [5] and is attributed to application of an improved resist cleaning step. The absorption loss measurement results at 1385 nm and 1550 nm are summarized in Table 1 and compared with values from the literature. The measured absorption ratio $\alpha$(1385 nm)/$\alpha$(1550 nm) in the current data is more consistent with the report in Ref. [11]. The high OH concentration in Ref. [6] could possibly account for this difference.

 figure: Fig. 2.

Fig. 2. Absorption loss measurement results. Absorption loss as measured by cavity-enhanced photothermal spectroscopy in the 1380-nm band. Red dots, data from Sample 1 in Fig. 1. Blue solid line, 2.4 ppm (weight) OH content level absorption line shape based on Ref. [6].

Download Full Size | PDF

Tables Icon

Table 1. Summary of Absorption Loss Rates at Both 1385 nm and 1550 nm from Fig. 2 and Taken from the Literaturea

To further confirm the OH ion impurity level and also study its spatial distribution over the resonator oxide thickness, SIMS depth profiling [15] was employed (see Fig. 3). The absolute OH ion concentration in this measurement is obtained from the ratio $^{17}$OH/$^{30}$Si, as calibrated against reference samples [16,17]. The OH concentration remains nearly constant along the depth direction with a value in reasonable agreement with the CEPS results in Fig. 2. This uniformity is believed to result from diffusion-driven equilibrium [11] during the long oxide growth process and the post high-temperature annealing process [2]. The higher OH concentration near the surfaces is consistent with the hydrophilic nature of silica causing water to be adsorbed on the surfaces [1820].

 figure: Fig. 3.

Fig. 3. Hydroxyl ion depth distribution measured using SIMS over the resonator thickness (Sample 1 in Fig. 1). Dashed line indicates 2.4-ppmw level inferred from data in Fig. 2.

Download Full Size | PDF

In summary, using thermal silica wedge microresonators and by measuring optical loss from 680 nm to 1550 nm, the OH ion absorption peak in the 1380-nm band stands out from baseline scattering loss. Absorption-only loss in the 1380-nm band was also measured using cavity enhanced photo-thermal spectroscopy. This gave a maximum absorption loss of 152 dB/km at 1385 nm. Fitting of the absorption spectrum gives an OH content level of 2.4 ppm weight, which is consistent with separate measurements performed on the resonator oxide by SIMS. Outside the OH overtone absorption band, surface scattering was determined to be the dominant loss mechanism. Record high $Q$-factors (limited by scattering) were measured for visible and near-visible wavelengths. Also, with further reduction of scattering loss, $Q$ values approaching 8 billion are feasible at 1550 nm. This value is higher than determined in previous measurements [5]. Overall, this investigation highlights the potential for even lower loss levels and higher $Q$-factors using thermal silica in photonic devices.

Methods: The scattering loss for silica wedge resonators is dominated by surface inhomogeneity scattering rather than bulk inhomogeneity scattering, as evident by the following two experimental observations. First, the optical loss is roughly inversely proportional to the bending radius [2], a characteristic feature for surface roughness scattering [18,21]. Second, with bending radius over 1 mm, the resonances of wedge resonators do not feature a splitting effect [22] even with $Q$-factor over 1 billion [2].

Surface scattering loss can be modeled as dipole radiation of surface inhomogeneities polarized by the local electrical field [9]. The equivalent current density induced by the corresponding dipole moments is

$$\vec J(\vec r)={-}i\omega \epsilon_0(n_0^2-n_1^2) \vec E(\vec r),$$
where $\omega$ is the angular frequency of the light, $\epsilon _0$ is the free space permittivity, $n_0,n_1$ are refractive indices of cavity and cladding materials, and $\vec E(\vec r)$ is the local electrical field. The corresponding radiation field and Poynting vector are given by
$$\vec{A}(\vec{r})=\frac{\mu e^{in_1 k_0 r}}{4\pi r}\oint e^{{-}in_1 k_0\hat{r}\cdot\vec{r'} }J(\vec{r'})h(\vec{r'})dS',$$
$$\mathcal{\vec S}(\vec r)=\hat r\frac{\omega n_1k_0}{2\mu}|\hat r\times \vec A(\vec{r})|^2,$$
where $k_0$ is the wavenumber of the light and $h(\vec {r})dS$ is the volume element of inhomogeneity. When adding up the radiation field from different parts of surface, their coherence is limited by the correlation length. Using the correlation function of surface roughness $\langle h(x)h(y)\rangle =\sigma ^2 R(\frac {x-y}{B})$, and assuming that the correlation length $B$ is small enough to ignore the variation of the field within this length, the Poynting vector field can be reduced to
$$\mathcal{\vec S}(\vec r)=\hat r\frac{\mu \omega n_1k_0}{2(4\pi r)^2}\sigma^2B^2\oint |\vec{J}(\vec{r'})\times \hat r|^2\tilde{R}[\vec\beta(\vec{r'})-n_1k_0\hat r_\|(\vec{r'})]dS',$$
where $\tilde R(\vec k)$ is Fourier transformation of $R(\vec x)$, $\vec \beta (\vec {r'})$ is the propagation constant of cavity mode at $\vec {r'}$, and $\hat r_\|(\vec {r'})$ is the component of $\hat r$ parallel to the surface at $\vec {r'}$.

In the calculation of $\mathcal {\vec S}$ and the radiation power, the parameters of the cavity mode, such as mode profiles and propagation constants, are acquired by numerical simulation (Lumerical MODE). The mode used for calculating the purple curve in Fig. 1 is the fundamental TE mode. The correlation function of surface roughness in the calculation is assumed to be a Gaussian function with the form $R(\frac {x-y}{B})=\frac {1}{B\sqrt {2\pi }}e^{-\frac {(x-y)^2}{2B^2}}$. The parameters of the surface roughness are acquired by atomic force microscopy (AFM) measurements, fitted with a Gaussian-type function, revealing that the lower surface has $\sigma$ of 2 nm, whereas the upper and wedge surfaces have $\sigma$ smaller than 0.5 nm [1,23].

Funding

Defense Advanced Research Projects Agency (HR0011-22-2-0009).

Acknowledgments

The authors acknowledge support provided by The Kavli Nanoscience Institute (KNI) at Caltech and Caltech Microanalysis Center (CMC) at Division of Geological $\&$ Planetary Sciences. The authors thank George R. Rossman and Yunbin Guan for fruitful discussions.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

REFERENCES

1. H. Lee, T. Chen, J. Li, K. Yang, S. Jeon, O. Painter, and K. J. Vahala, Nat. Photonics 6, 369 (2012). [CrossRef]  

2. L. Wu, H. Wang, Q. Yang, Q.-X. Ji, B. Shen, C. Bao, M. Gao, and K. Vahala, Opt. Lett. 45, 5129 (2020). [CrossRef]  

3. W. Jin, Q.-F. Yang, L. Chang, B. Shen, H. Wang, M. A. Leal, L. Wu, M. Gao, A. Feshali, M. Paniccia, K. J. Vahala, and J. E. Bowers, Nat. Photonics 15, 346 (2021). [CrossRef]  

4. M. W. Puckett, K. Liu, N. Chauhan, Q. Zhao, N. Jin, H. Cheng, J. Wu, R. O. Behunin, P. T. Rakich, K. D. Nelson, and D. J. Blumenthal, Nat. Commun. 12, 934 (2021). [CrossRef]  

5. M. Gao, Q.-F. Yang, Q.-X. Ji, H. Wang, L. Wu, B. Shen, J. Liu, G. Huang, L. Chang, W. Xie, S.-P. Yu, S. B. Papp, J. E. Bowers, T. J. Kippenberg, and K. J. Vahala, Nat. Commun. 13, 3323 (2022). [CrossRef]  

6. O. Humbach, H. Fabian, U. Grzesik, U. Haken, and W. Heitmann, J. Non-Cryst. Solids 203, 19 (1996). [CrossRef]  

7. M. Cai, O. Painter, and K. J. Vahala, Phys. Rev. Lett. 85, 74 (2000). [CrossRef]  

8. S. Spillane, T. Kippenberg, O. Painter, and K. Vahala, Phys. Rev. Lett. 91, 043902 (2003). [CrossRef]  

9. T. Barwicz and H. A. Haus, J. Lightwave Technol. 23, 2719 (2005). [CrossRef]  

10. R. Olshansky, Rev. Mod. Phys. 51, 341 (1979). [CrossRef]  

11. G. A. Thomas, B. I. Shraiman, P. F. Glodis, and M. J. Stephen, Nature 404, 262 (2000). [CrossRef]  

12. T. Li, Optical Fiber Communications: Fiber Fabrication (Academic Press, 1985).

13. G. Keiser, Optical Fiber Communications4th edition (McGraw-Hill Education, 2010).

14. T. Carmon, L. Yang, and K. J. Vahala, Opt. Express 12, 4742 (2004). [CrossRef]  

15. R. G. Wilson, F. A. Stevie, and C. W. Magee, Secondary Ion Mass Spectrometry: A Practical Handbook for Depth Profiling and Bulk Impurity Analysis (Wiley, 1989).

16. E. Hauri, J. Wang, J. E. Dixon, P. L. King, C. Mandeville, and S. Newman, Chem. Geol. 183, 99 (2002). [CrossRef]  

17. J. L. Mosenfelder, M. Le Voyer, G. R. Rossman, Y. Guan, D. R. Bell, P. D. Asimow, and J. M. Eiler, Am. Mineral. 96, 1725 (2011). [CrossRef]  

18. M. L. Gorodetsky, A. A. Savchenkov, and V. S. Ilchenko, Opt. Lett. 21, 453 (1996). [CrossRef]  

19. H. Rokhsari, S. Spillane, and K. Vahala, Appl. Phys. Lett. 85, 3029 (2004). [CrossRef]  

20. S. Honari, S. Haque, and T. Lu, Appl. Phys. Lett. 119, 031107 (2021). [CrossRef]  

21. M. L. Gorodetsky, A. D. Pryamikov, and V. S. Ilchenko, J. Opt. Soc. Am. B 17, 1051 (2000). [CrossRef]  

22. T. J. Kippenberg, S. M. Spillane, and K. J. Vahala, Opt. Lett. 27, 1669 (2002). [CrossRef]  

23. H. Lee, T. Chen, J. Li, O. Painter, and K. J. Vahala, Nat. Commun. 3, 867 (2012). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1.
Fig. 1. Optical loss measured from 680 nm to 1550 nm reveals the existence of the OH absorption peak near 1380 nm. The data are obtained by optical $Q$-factor measurement. The purple solid curve is the predicted loss from surface-roughness scattering (see text). The light blue solid curve is the estimated OH absorption loss in the 1380-nm band (corresponding to 2.4 ppm weight). Dashed contours correspond to different $Q$-factor values. Inset shows typical optical micrograph of the optical resonator used in the measurements.
Fig. 2.
Fig. 2. Absorption loss measurement results. Absorption loss as measured by cavity-enhanced photothermal spectroscopy in the 1380-nm band. Red dots, data from Sample 1 in Fig. 1. Blue solid line, 2.4 ppm (weight) OH content level absorption line shape based on Ref. [6].
Fig. 3.
Fig. 3. Hydroxyl ion depth distribution measured using SIMS over the resonator thickness (Sample 1 in Fig. 1). Dashed line indicates 2.4-ppmw level inferred from data in Fig. 2.

Tables (1)

Tables Icon

Table 1. Summary of Absorption Loss Rates at Both 1385 nm and 1550 nm from Fig. 2 and Taken from the Literaturea

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

J ( r ) = i ω ϵ 0 ( n 0 2 n 1 2 ) E ( r ) ,
A ( r ) = μ e i n 1 k 0 r 4 π r e i n 1 k 0 r ^ r J ( r ) h ( r ) d S ,
S ( r ) = r ^ ω n 1 k 0 2 μ | r ^ × A ( r ) | 2 ,
S ( r ) = r ^ μ ω n 1 k 0 2 ( 4 π r ) 2 σ 2 B 2 | J ( r ) × r ^ | 2 R ~ [ β ( r ) n 1 k 0 r ^ ( r ) ] d S ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.