Thermo-Optical Sensitivity of Whispering Gallery Modes in As2S3 Chalcogenide Glass Microresonators

Glass microresonators with whispering gallery modes (WGMs) have a lot of diversified applications, including applications for sensing based on thermo-optical effects. Chalcogenide glass microresonators have a noticeably higher temperature sensitivity compared to silica ones, but only a few works have been devoted to the study of their thermo-optical properties. We present experimental and theoretical studies of thermo-optical effects in microspheres made of an As2S3 chalcogenide glass fiber. We investigated the steady-state and transient temperature distributions caused by heating due to the partial thermalization of the pump power and found the corresponding wavelength shifts of the WGMs. The experimental measurements of the thermal response time, thermo-optical shifts of the WGMs, and heat power sensitivity in microspheres with diameters of 80–380 µm are in a good agreement with the theoretically predicted dependences. The calculated temperature sensitivity of 42 pm/K does not depend on diameter for microspheres made of commercially available chalcogenide fiber, which may play an important role in the development of temperature sensors.


Introduction
Dielectric microresonators with whispering gallery modes (WGMs) have plenty of diversified applications [1], among which sensing applications play a significant role [2][3][4][5][6][7][8][9][10][11]. Microresonators with certain advantages and disadvantages are produced from different materials using special technologies [12]. Among dielectric microresonators, glass microresonators are widely used, and silica ones are the most common since they can be manufactured in a reproducible manner from standard telecommunication fibers [13][14][15][16][17][18]. However, the use of other non-silica glasses with very different properties is also actively attracting the attention of different research groups since this allows expanding the boundaries of the possibilities and applications for microresonators and improving the characteristics of the devices based on them.
Chalcogenide glasses have been actively investigated for implementing nonlinear optical effects, since such glasses have huge cubic nonlinearities, ultrawide transparency ranges, and physical and chemical properties suitable for producing optical microresonators. For example, Raman lasing [19,20], optical parametric oscillations [21], the generation of Raman-Kerr optical frequency combs [22], and Brillouin generation [23] have been attained in chalcogenide microresonators.
Effects based on thermo-optical nonlinearity [24] have also been observed in chalcogenide microresonators. For instance, a laser wavelength in an Nd-doped chalcogenide microsphere was tuned when the pump power was changed; this fact was explained by the influence of nonlinear and/or thermal effects [25]. In [22], the tunability of a Raman wavelength was obtained with a change in temperature due to shifts of the WGM resonances.
The start of Raman lasing was attained in [26] with the use of an auxiliary laser diode for microsphere heating to control the WGM shifts. The measurement of temperature sensitivity based on the same WGM wavelength shift effect with increasing ambient temperature in a chalcogenide microsphere was demonstrated in [11]. However, there have been almost no systematic studies of thermo-optical effects, including the investigation of temperature distributions in a microresonator during pump thermalization, thermal response times and the influence of these factors on the resonant wavelengths of WGMs, or on the influence of geometrical parameters on chalcogenide microresonators. However, such studies are very important for sensing applications, primarily for measuring ambient temperature as chalcogenide microresonators have higher sensitivity compared to, for example, silica ones [11]. Such research is useful in the search for an optimal design, choosing optimal parameters for sensor elements, and studying their limiting capabilities. Moreover, thermooptical control using external heating or heating with an auxiliary laser diode can be used to quickly switch on/off nonlinear optical or laser generation and tuning the generated or transformed wavelengths.
Here we present detailed experimental and theoretical studies of thermo-optical effects in As 2 S 3 microspheres. We investigated the temperature distributions caused by heating due to partial thermalization of the pump power and the corresponding wavelength shifts of WGMs. The temperature sensitivity and heat power sensitivity were studied.

Experimental Materials and Methods
The simplified experimental scheme is shown in Figure 1a. We performed studies with microspheres of different diameters, d, in the 80-380 µm range, made of a commercially available fiber based on As 2 S 3 chalcogenide glass. Images of the microspheres used in all our experiments are presented in Figure 1b. The manufacturing method was based on melting the end of the fiber using a resistive microheater with the formation of a sphere under the action of surface tension forces. The details are given in our previous article [20]. CW pump laser radiation was coupled into the microsphere through a fiber taper made of a commercially available silica telecommunication fiber by heating and stretching [16,27]. We used a commercially available tunable telecommunication C-band laser (Pure Photonics, PLCC550-180-60) with a linewidth of 10 kHz and adjustable power up to 18 dBm (63 mW) supporting the frequency sweeping mode. The loaded quality factors (Q-factors) of microspheres were measured, and thermo-optical dynamic processes occurring during thermalization of the radiation power were investigated with the use of this laser. The central wavelength of the pump laser in all our experiments was set at about 1.55 µm.
The Q-factors of the microspheres were measured by the dynamic resonance scanning method, when the pump laser frequency was swept linearly at a rate of 10 GHz/s, and the resonant dip was recorded using an oscilloscope (Tektronix MSO-64 6-BW-2500, Beaverton, OR USA). The typical Q-factors for all samples, regardless of their diameter, were Q~10 6 . An example of the measured data and the corresponding Lorentz fit are shown in Figure 1c. The pump power was attenuated to hundreds of nW to avoid nonlinear effects during measurements of the resonant dip.
In the study of thermo-optical effects, the pump laser power was set to be sufficiently high (~20-40 mW). Additionally, low-power broadband radiation from an erbium-doped fiber source of amplified spontaneous emission (ASE) was coupled with the microsphere to monitor WGM wavelengths. The radiation transformed in the microsphere was extracted through the same taper and then divided into two diagnostic channels. In one channel, the spectra were measured using an optical spectrum analyzer (OSA, Yokogawa AQ6370D, Tokyo, Japan), which allowed us to determine the WGM wavelength shifts. In the other channel, the signal was received by a photodetector (PD, Thorlabs PDA015C) and an oscilloscope, which made it possible to measure the transmitted power. When detuning of the pump laser frequency was changed, the absorbed power changed too, which could be estimated from the variations in the transmitted power. In these experiments, we did not observe nonlinear optical conversion of the pump frequency due to the Kerr and Raman nonlinearities, so we assumed that the absorbed power corresponded to the thermalized power. The pump laser was controlled by a computer. We processed and analyzed experimental data using specially developed algorithms and computer codes written on their basis. channels. In one channel, the spectra were measured using an optical spectrum a (OSA, Yokogawa AQ6370D,Tokyo, Japan), which allowed us to determine the wavelength shifts. In the other channel, the signal was received by a photodetect Thorlabs PDA015C) and an oscilloscope, which made it possible to measure th mitted power. When detuning of the pump laser frequency was changed, the ab power changed too, which could be estimated from the variations in the tran power. In these experiments, we did not observe nonlinear optical conversion pump frequency due to the Kerr and Raman nonlinearities, so we assumed that sorbed power corresponded to the thermalized power. The pump laser was contro a computer. We processed and analyzed experimental data using specially dev algorithms and computer codes written on their basis.

Theoretical Model and Methods
To theoretically study thermo-optical processes and explain experimental and dependences, we developed and applied a theoretical model based on find temperature distributions in As2S3 glass microspheres during pump power ther tion, followed by calculating the WGM eigenfrequencies and finding the wav shifts of WGMs Δλ when changing the system parameters. In our earlier work, we similar model to study thermo-optical effects in tellurite glass microspheres [28].
The microresonator geometry used in the numerical simulations is shown in 2a. In our model, we assumed that pump thermalization occurs uniformly in a near the sphere equator and the effective area Seff and volume Veff of this region ar to the effective area and volume of the fundamental TE mode of the resonator, tively, at a central wavelength of about 1.55 μm (Figure 2b). To calculate Seff and found the fields of the fundamental modes using expressions from [29].

Theoretical Model and Methods
To theoretically study thermo-optical processes and explain experimental results and dependences, we developed and applied a theoretical model based on finding the temperature distributions in As 2 S 3 glass microspheres during pump power thermalization, followed by calculating the WGM eigenfrequencies and finding the wavelength shifts of WGMs ∆λ when changing the system parameters. In our earlier work, we used a similar model to study thermo-optical effects in tellurite glass microspheres [28].
The microresonator geometry used in the numerical simulations is shown in Figure 2a. In our model, we assumed that pump thermalization occurs uniformly in a region near the sphere equator and the effective area S eff and volume V eff of this region are equal to the effective area and volume of the fundamental TE mode of the resonator, respectively, at a central wavelength of about 1.55 µm (Figure 2b). To calculate S eff and V eff , we found the fields of the fundamental modes using expressions from [29]. To determine the steady-state and dynamic temperature distributions caused by partial thermalization of the pump power, we solved (by the finite element method) the heat equation numerically [30,31]: where ΔT is the temperature increase relative to the ambient temperature, ρ is the As2S3 glass density, cp is the As2S3 glass heat capacity at constant pressure, k is the As2S3 glass thermal conductivity, q is heat flow, and Q is the heat source uniformly distributed over the volume Veff ( Figure 2). The natural convection in air was set as a boundary condition on all surfaces, except for the end of the fiber opposite to the microsphere, where the temperature was fixed (although the boundary condition on this surface had practically no effect on the temperature distributions inside the microsphere). The boundary condition of natural convection was [30]: where nnorm is a normalized outward-pointing vector, h ≈ 2 kair/d for a sphere, and kair is the air thermal conductivity (a more precise expression for h used in our numerical simulations is given in Appendix A) [30]. Constants describing the material properties of As2S3 glass microspheres and other parameters are listed in Table 1.  To determine the steady-state and dynamic temperature distributions caused by partial thermalization of the pump power, we solved (by the finite element method) the heat equation numerically [30,31]: where ∆T is the temperature increase relative to the ambient temperature, ρ is the As 2 S 3 glass density, c p is the As 2 S 3 glass heat capacity at constant pressure, k is the As 2 S 3 glass thermal conductivity, q is heat flow, and Q is the heat source uniformly distributed over the volume V eff ( Figure 2). The natural convection in air was set as a boundary condition on all surfaces, except for the end of the fiber opposite to the microsphere, where the temperature was fixed (although the boundary condition on this surface had practically no effect on the temperature distributions inside the microsphere). The boundary condition of natural convection was [30]: where n norm is a normalized outward-pointing vector, h ≈ 2 k air /d for a sphere, and k air is the air thermal conductivity (a more precise expression for h used in our numerical simulations is given in Appendix A) [30]. Constants describing the material properties of As 2 S 3 glass microspheres and other parameters are listed in Table 1.
To find the eigenfrequencies of the fundamental TE modes, we used the characteristic equation obtained from Maxwell's equations based on the well-known approach [29]: where k 0 = 2π/λ, c is the speed of light in vacuum, λ is the wavelength, d is the microsphere diameter, n is the wavelength-dependent As 2 S 3 glass refractive index (defined in Table 1), l is the polar WGM index, J x is the Bessel function of order x, H x (1) is the Hankel function of the 1st kind of order x, and the prime means the derivative with respect to the argument (nk 0 d/2 or k 0 d/2). Equation (4) has multiple roots corresponding to different radial WGM indices (q = 1 for the fundamental WGM). We numerically solved Equation (4) with allowance for the As 2 S 3 glass dispersion and the dependence of n and d on the temperature increase where n 0 is a linear refractive index; ∆n = (dn/dT)∆T mode and ∆d are changes in the refractive indexes of the As 2 S 3 glass and microsphere diameters, respectively, under heating; dn/dT and ε are the thermo-optical and thermal expansion coefficients, respectively, of the As 2 S 3 glass; and ∆T mode and ∆T av are the temperature increases averaged over the regions corresponding to the effective WGM size and the overall the sphere, respectively. Therefore, in our model we assumed that the temperature increase averaged over the sphere affects the thermal expansion, but ∆n is important only in the region where the fundamental WGMs are located.

Theoretical Study of Temperature Distributions
First of all, we studied the temperature distributions arising in microspheres of different diameters upon partial thermalization of the pump power in the region corresponding to the effective size of the fundamental WGM at a wavelength of about 1.55 µm. As an example, the steady-state temperature distribution in a microsphere with a diameter of 140 µm is shown in Figure 3a. The maximum temperature is reached in the region where the heating source is located and decreases when moving away from the source. When the pump is switched on at the time t = 0, the microsphere begins to heat up; the curves for the average temperature increases ∆T mode and ∆T av are shown in Figure 3b. After a few hundred ms, the average temperatures reach steady-state values. After the pump is switched off at the moment of time t off , marked by the vertical dotted gray line, the microsphere begins to cool down, and the average temperatures relax to the ambient temperature (Figure 3b). The time dependence of the temperature decrease averaged over the sphere is well approximated by an exponential function with a characteristic time t 0 ('fit' in Figure 3b): We numerically simulated the dynamic cooling of microspheres of different diameters. The relaxation time t 0 as a function of d is plotted in Figure 3c. Note that t 0 is well approximated by the function t 0 = C 1 ·d 2 , where C 1 = 5.2 × 10 −6 s·µm −2 ('fit' in Figure 3c). Moreover, we found an analytical solution to the problem of cooling an ideal sphere with reasonable approximations, which gives an exponential decay for the average temperature and the dependence t 0~d 2 . The found analytical solution allowed us to independently estimate the proportionality coefficient as 5.0 × 10 −6 s·µm −2 , which agrees very well with the numerically found coefficient C 1 . The analytical solution is given in Appendix A. estimate the proportionality coefficient as 5.0 × 10 −6 s⸱ μm −2 , which agrees very well with the numerically found coefficient C1. The analytical solution is given in Appendix A. As for the temperature averaged over the WGM region, its dynamics are more complex, and within the framework of our model, are not described by a single decaying exponential dependence on time. Namely, when the pump is switched off, ΔTmode first rapidly decreases to a value of about ΔTav (t = toff) in time << t0, and then both average values ΔTav and ΔTmode decrease exponentially with a characteristic time t0 (Figure 3b).
Next, we calculated steady-state temperature distributions. The average temperatures increases ΔTav and ΔTmode as functions of a thermalized power, P, and a microsphere As for the temperature averaged over the WGM region, its dynamics are more complex, and within the framework of our model, are not described by a single decaying exponential dependence on time. Namely, when the pump is switched off, ∆T mode first rapidly decreases to a value of about ∆T av (t = t off ) in time << t 0 , and then both average values ∆T av and ∆T mode decrease exponentially with a characteristic time t 0 (Figure 3b). Next, we calculated steady-state temperature distributions. The average temperatures increases ∆T av and ∆T mode as functions of a thermalized power, P, and a microsphere diameter, d, are presented in Figure 3d,e, respectively. Note that, for fixed d, these dependencies are almost linear (∆T~P). The corresponding examples for d = 140 µm are plotted in Figure 3f. For fixed P, an average temperature increase is inversely proportional to the microsphere diameter. The corresponding examples for P = 1 mW are plotted in Figure 3g.

Theoretical Study of Steady-State Wavelength Shifts of the WGMs
Subsequently, we numerically simulated the wavelength shifts of the WGMs for a uniform increase in the temperature of the microsphere, which occurs with an increase in the ambient temperature. In this case, we assumed that the WGMs can be detected with the help of radiation of a very low power which does not affect the temperature distribution. The sensitivity ∆λ/∆T, as a function of microsphere diameter, is shown in Figure 4. It is seen that the sensitivity is practically independent of the diameter, which is important for sensing. This peculiarity can be explained by using the resonance condition l·λ = πdn eff , where n eff is the WGM effective refractive index, which is primarily a function of λ but also weakly depends on the structure. Using Equations (5) and (6), an approximate relation between ∆λ and ∆T can be obtained:

Theoretical Study of Steady-State Wavelength Shifts of the WGMs
Subsequently, we numerically simulated the wavelength shifts of th uniform increase in the temperature of the microsphere, which occurs with the ambient temperature. In this case, we assumed that the WGMs can be the help of radiation of a very low power which does not affect the temp bution. The sensitivity Δλ/ΔT, as a function of microsphere diameter, is sh 4. It is seen that the sensitivity is practically independent of the diameter portant for sensing. This peculiarity can be explained by using the reson l•λ = πdneff, where neff is the WGM effective refractive index, which is prima of λ but also weakly depends on the structure. Using Equations (5) and ( mate relation between Δλ and ΔT can be obtained: Therefore, Δλ/ΔT is almost a constant value for modes near a give Namely, when using As2S3 chalcogenide microspheres having diameter tens to several hundreds of microns, the predicted sensitivity is ~42 pm small perturbations of temperature sensitivity are observed for microres ferent sizes (Figure 4). The explanation is that the eigenfrequencies and cor found by solving Equation (4)   Next, we studied the wavelength shifts of the WGMs occurring as a perature increase due to the thermalization of the pump power. The num lated Δλ as a function of P and d is shown in Figure 5a. For a microsphere Δλ is an almost linear function of P. The example for d = 140 μm is plotte Moreover, at a fixed power, Δλ is inversely proportional to d (Figure 5c). Therefore, ∆λ/∆T is almost a constant value for modes near a given wavelength. Namely, when using As 2 S 3 chalcogenide microspheres having diameters from several tens to several hundreds of microns, the predicted sensitivity is~42 pm/K. However, small perturbations of temperature sensitivity are observed for microresonators of different sizes (Figure 4). The explanation is that the eigenfrequencies and corresponding n eff found by solving Equation (4) are slightly varied for the fundamental WGMs (nearest to 1.55 µm) in microspheres with different diameters.
Next, we studied the wavelength shifts of the WGMs occurring as a result of temperature increase due to the thermalization of the pump power. The numerically simulated ∆λ as a function of P and d is shown in Figure 5a. For a microsphere of a fixed size, ∆λ is an almost linear function of P. The example for d = 140 µm is plotted in Figure 5b. Moreover, at a fixed power, ∆λ is inversely proportional to d (Figure 5c).

Experimental and Theoretical Study of Temporal Dynamics of the Wavelength Shifts of the WGMs
We investigated the temporal dynamics of the WGM wavelength shifts after the pump was switched off. In this experiment, a steady-state temperature distribution was established under the action of a CW pump and the WGM wavelengths had certain values. Then, the pump was switched off, after which the microsphere cooled down and the WGM wavelengths shifted to the short-wavelength range.
We measured the temporal dynamics of Δλ for all samples of the produced microresonators. An example of the measurements for a microresonator with a diameter of 140 μm is shown in Figure 6a. The pump was switched off at t = 0. The decreasing function Δλ(t) can be approximated by the exponential dependence on time We performed direct simulations in the framework of the model described in Section 2.2 and obtained the numerical dependence of Δλ for t after the pump was switched off. An example of Δλ(t) for d = 140 μm is plotted in Figure 6b. Experimental and numerical results agree very well (compare Figure 6a,b). However, in our theoretical model, Δλ(t) is not an ideal exponential function because Δλ depends on both ΔTav and ΔTmode. The explanation based on the results obtained in Section 3.1 is that ΔTav is well-fitted by an exponential dependence, but ΔTmode is a more complex function. Consequently, the simulated values of tΔλ are slightly lower than t0 (the difference is about 10%). The ex-

Experimental and Theoretical Study of Temporal Dynamics of the Wavelength Shifts of the WGMs
We investigated the temporal dynamics of the WGM wavelength shifts after the pump was switched off. In this experiment, a steady-state temperature distribution was established under the action of a CW pump and the WGM wavelengths had certain values. Then, the pump was switched off, after which the microsphere cooled down and the WGM wavelengths shifted to the short-wavelength range.
We measured the temporal dynamics of ∆λ for all samples of the produced microresonators. An example of the measurements for a microresonator with a diameter of 140 µm is shown in Figure 6a. The pump was switched off at t = 0. The decreasing function ∆λ(t) can be approximated by the exponential dependence on time ∆λ = ∆λ| t=0 · exp(−t/t ∆λ ).

Experimental and Theoretical Study of Temporal Dynamics of the Wavelength Shifts of the WGMs
We investigated the temporal dynamics of the WGM wavelength shifts after the pump was switched off. In this experiment, a steady-state temperature distribution was established under the action of a CW pump and the WGM wavelengths had certain values. Then, the pump was switched off, after which the microsphere cooled down and the WGM wavelengths shifted to the short-wavelength range.
We measured the temporal dynamics of Δλ for all samples of the produced microresonators. An example of the measurements for a microresonator with a diameter of 140 μm is shown in Figure 6a. The pump was switched off at t = 0. The decreasing function Δλ(t) can be approximated by the exponential dependence on time We performed direct simulations in the framework of the model described in Section 2.2 and obtained the numerical dependence of Δλ for t after the pump was switched off. An example of Δλ(t) for d = 140 μm is plotted in Figure 6b. Experimental and numerical results agree very well (compare Figure 6a,b). However, in our theoretical model, Δλ(t) is not an ideal exponential function because Δλ depends on both ΔTav and ΔTmode. The explanation based on the results obtained in Section 3.1 is that ΔTav is well-fitted by an exponential dependence, but ΔTmode is a more complex function. Consequently, the simulated values of tΔλ are slightly lower than t0 (the difference is about 10%). The ex- Figure 6. Experimental (a) and numerically simulated (b) dependences of thermo-optical shift in WGMs ∆λ at time t. Thermalized power is 0.9 mW for t < 0 and the pump is switched off at t = 0. 'Fit' is approximation of ∆λ by exponential decay with characteristic time t ∆λ for t > 0; (c) t ∆λ dependence on microsphere diameter obtained by experimental measurements, by direct simulations and fitted by t ∆λ = C 2 ·d 2 , where C 2 = 4.6 × 10 −6 s·µm −2 .
We performed direct simulations in the framework of the model described in Section 2.2 and obtained the numerical dependence of ∆λ for t after the pump was switched off. An example of ∆λ(t) for d = 140 µm is plotted in Figure 6b. Experimental and numerical results agree very well (compare Figure 6a,b). However, in our theoretical model, ∆λ(t) is not an ideal exponential function because ∆λ depends on both ∆T av and ∆T mode . The explanation based on the results obtained in Section 3.1 is that ∆T av is well-fitted by an exponential dependence, but ∆T mode is a more complex function. Conse-quently, the simulated values of t ∆λ are slightly lower than t 0 (the difference is about 10%). The experimentally measured and numerically simulated t ∆λ values for different microsphere diameters are shown in Figure 6c. The corresponding approximation t ∆λ = C 2 ·d 2 , where C 2 = 4.6 × 10 −6 s·µm −2 is also plotted.

Experimental and Theoretical Study of ∆λ/P
Finally, we measured the heat power sensitivity ∆λ/P for all produced microspheres (Figure 7). The simulations demonstrate that ∆λ/P is proportional to 1/d. Note that the presented experimental data and numerical dependences agree very well. The main sources of uncertainties in our measurements are the OSA resolution of 0.02 nm, its limited scanning speed, and some drifts in the photodetector signal used for power measurements. For the microspheres with a high density of WGMs, wavelength tracking errors may introduce additional uncertainty into the WGM thermal shift measurements larger than the OSA resolution.

Experimental and Theoretical Study of Δλ/P
Finally, we measured the heat power sensitivity Δλ/P for all produ (Figure 7). The simulations demonstrate that Δλ/P is proportional to 1 presented experimental data and numerical dependences agree very sources of uncertainties in our measurements are the OSA resolution of ited scanning speed, and some drifts in the photodetector signal used urements. For the microspheres with a high density of WGMs, wavele rors may introduce additional uncertainty into the WGM thermal shi larger than the OSA resolution.

Discussion
We produced chalcogenide glass microspheres with diameters of typical Q-factors of 10 6 from a commercially available As2S3 fiber and i thermo-optical properties experimentally and theoretically. The theoret based on numerical simulations of the heat equation using the finite-ele experimental-like geometry. Average temperature increases were calc into account aimed at finding the eigenfrequencies of the WGMs an shifts caused by thermal expansion of the microsphere and the tempera of the refractive index. We investigated the steady-state and transient tributions caused by heating due to the partial thermalization of the p found the corresponding wavelength shifts of the WGMs. We showed average steady-state temperature increases are proportional to the th power and inversely proportional to the microsphere diameter. We de mentally, numerically, and even analytically that the thermal resp switching off the pump is proportional to the microsphere diameter coefficient is about 5 × 10 −6 ⸱ s μm −2 . We demonstrated through numerical simulations that for As2S3 crospheres the temperature sensitivity Δλ/ΔT = 42 pm/K does not depe

Discussion
We produced chalcogenide glass microspheres with diameters of 80-380 µm and typical Q-factors of 10 6 from a commercially available As 2 S 3 fiber and investigated their thermo-optical properties experimentally and theoretically. The theoretical studies were based on numerical simulations of the heat equation using the finite-element method for experimental-like geometry. Average temperature increases were calculated and taken into account aimed at finding the eigenfrequencies of the WGMs and thermo-optical shifts caused by thermal expansion of the microsphere and the temperature dependence of the refractive index. We investigated the steady-state and transient temperature distributions caused by heating due to the partial thermalization of the pump power and found the corresponding wavelength shifts of the WGMs. We showed theoretically that average steady-state temperature increases are proportional to the thermalized pump power and inversely proportional to the microsphere diameter. We determined experimentally, numerically, and even analytically that the thermal response time after switching off the pump is proportional to the microsphere diameter squared, and the coefficient is about 5 × 10 −6 s·µm −2 .
We demonstrated through numerical simulations that for As 2 S 3 chalcogenide microspheres the temperature sensitivity ∆λ/∆T = 42 pm/K does not depend on diameters ranging from several tens to several hundreds of microns, and a simple explanation in terms of the effective refractive index of the mode was provided. Note that this sensitivity is higher than that for microspheres made from many other glasses. For example, for a gallium-germanium-antimony-sulfide chalcogenide glass microsphere with a diameter of about 100 µm, the measured sensitivity was 28 pm/K [11]. For silica microspheres, the experimental sensitivity at room temperature was about 10 pm/K and was almost independent of the diameter [10]. In addition, although we did not directly measure ∆λ/ ∆T for an external increase in ambient temperature, we measured heat power sensitivity ∆λ/P and obtained good agreement with the theoretically predicted dependence ∆λ/P~1/d. Therefore, this fact confirms that our numerical model is well calibrated and gives reliable results.
The microresonators studied here can serve as microheaters and sensors simultaneously, while the measurement of the WGM shifts provides quite accurate temperature readings. Chalcogenide glasses ensure extended transparency in the mid-IR range, where the spectral fingerprints of many important chemical substances are located. The chalcogenide glass-based microresonators may be used as multipurpose devices combining heating/temperature-measuring functionalities and spectral sensing of their environment. Currently, fibers based on chalcogenide As 2 S 3 glass are commercially available, so microsphere manufacturing based on them can be a routine, fast, and inexpensive process. The independence of temperature sensitivity and microsphere diameter means that it is not necessary to control the diameter of the samples for such sensors. Thus, the demonstrated results are quite universal and may be practically significant.