Optical heating of resonant nanostructures is one of the key issues in modern nanophotonics, being either harmful or desirable effect depending on the applications. Despite a linear regime of light-to-heat conversion being well-studied both for metal and semiconductor resonant systems is generalized as a critical coupling condition, the clear strategy to optimize optical heating upon high-intensity light irradiation is still missing. This work proposes a simple analytical model for such a problem, taking into account material properties changes caused by the heating. It allows us to derive a new general critical coupling condition for the nonlinear case, requiring a counterintuitive initial spectral mismatch between the pumping light frequency and the resonant one. Based on the suggested strategy, we develop an optimized design for efficient nonlinear optical heating, which employs a cylindrical nanoparticle supporting the quasi bound state in the continuum mode (quasi-BIC or so-called ‘super-cavity mode’) excited by the incident azimuthal vector beam. Our approach provides a background for various nonlinear experiments related to optical heating and bistability, where self-action of the intense laser beam can change resonant properties of the irradiated nanostructure.
The strong resonant response of all-dielectric resonant nanosystems in the visible and infrared region along with the diversity of their optical properties opens the way for various applications in nonlinear and laser optics . Remarkably, such low-loss nanostructures can support high-Q modes under proper management of radiative losses  even in single nanocavities on dielectric substrates owing to excitation of super-cavity modes , or in the one- and two-dimensional arrays of resonant nanostructures [4–7], which are often related to the bound state in the continuum  (BIC) or quasi-BIC in the systems of finite size. Even though in many applications, the inherent optical heating of resonators is considered as a parasitic effect, strong thermooptical coefficients of all-dielectric materials can be utilized for developing thermally tunable and reconfigurable nanophotonic devices . Indeed, the achievements of thermo-nanophotonics based on all-dielectric  nanostructures (similarly to thermoplasmonics [11–13]) showed that efficient nano- and micro-scale sources of heat are in strong demand in various nanoscience applications.
From this point of view, the problem of efficient heating of all-dielectric nanostructures requires special optimization, depending on the final application of the nanophotonic design. Based on this approach and employing advanced methods of nanothermometry, the case of linear optical heating of all-dielectric Mie-resonant nanoparticles was successfully described theoretically [14–16] and demonstrated experimentally [17–19]. However, the situation becomes much more complicated once the optical nonlinearity due to thermorefractive effects is taken into account . Recently, there has been a significant progress in the field of nonlinear thermal nanophotonics with resonant all-dielectric  and plasmonic  systems, showing large values of thermal nonlinearity in single resonant nanostructures [23–26]. In this prospective, the problem of efficient heating is tightly connected to maximizing the absorption of resonators [27, 28]. While the critical coupling condition is required for the most optimal linear heating regime , the exact conditions of maximal optical heating in the nonlinear regime, when the elevated temperature drives change of the real and imaginary parts of the refractive index, are yet to be identified.
In this work, we firstly develop a simple analytical model for the nonlinear optical heating of a single-mode resonator supporting efficient light-to-heat conversion and bistable regime of operation, and then employ numerical simulations to propose a realistic design based on the super-cavity. The developed formalism and numerical design dealing with doped silicon cylindrical particles (see Figure 1) allow for the creation of novel thermooptical nanophotonic devices for optical switching and signal processing, as well as useful for the experiments where heat generation in resonant nanostructures should be minimized.
2 Optical heating of a nanoresonator
Optical heating of matter is a rather complex process, which combines free carrier generation, their interaction with light and phonons, and transport of the phonons across the material . However, in the steady-state continuous wave (CW) excitation of subwavelength resonators considered in this paper the situation becomes much simpler. For a large class of nanophotonic structures such as nanoparticles and nanoantennas made of material with high thermal conductivity (metals or semiconductors) as compared to surrounding material (air, glass, and water) the steady-state temperature appears to be homogeneously distributed along the nanostructure  and can be described by the expression 
where P abs is the light power absorbed inside the resonator due nonradiative (ohmic) losses, R is the typical radius of the sphere of the same volume as the nanostructure, η eff is the effective thermal conductivity of surrounding media, and C is a constant which is defined by the shape of the nanoantenna and geometry of the problem.
The problem of efficient heating of nanoresonator is, thus, reduced to maximization of the absorbed power, which has its fundamental limitation . This limit, however, can be reached in the critical coupling regime when scattered and absorbed power are balanced, which paves the way for optimizing the absorption losses . However, despite the critical coupling condition is usually considered for a linear system, it still holds in strongly nonlinear Kerr systems  or two-level systems with saturated absorption [32, 33] allowing for coherent perfect absorption. Here, we discuss the nonlinear critical coupling in the context of optical heating of single nanoresonators. We start by considering a single-mode nonlinear resonator within a coupled-mode theory  as a toy model to identify the main conditions to reach the maximal absorption condition.
2.1 Critical coupling in linear regime
We consider a subwavelength resonator with resonance at the frequency ω 0 having radiative and non-radiative (ohmic) losses rates γ r and γ nr correspondingly. The latter are responsible for light absorption inside the resonator and its consequent heating. Indeed, the absorbed power P abs is proportional to nonradiative losses γ nr and the total electric energy W stored inside the resonator: P abs = 2γ nr W. A single mode resonator can be well-described by the temporal coupled mode theory  and the amplitude of the excited mode can be given by the equation :
where ω 0 is the eigenfrequency of the resonator, γ = γ r + γ nr is the total loss rate, corresponds the amplitude of the incident wave. In the stationary regime under harmonic excitation , the spectral amplitude intensity of the field can be found from Eq. (2):
where Δω = ω 0 − ω is the detuning between the pumping field and the mode frequency. The mode intensity equals the total energy stored inside the resonator W = |a|2. Thus, the absorbed power P abs = 2γ nr|a|2 also has resonant Lorentz spectral profile and at the resonance reaches its maximal value
In these terms, the critical coupling condition manifests itself in equal radiative and non-radiative losses γ r = γ nr. Indeed, once this condition is fulfilled, the absorbed power tends to its maximal value of P abs → |f|2/2.
2.2 Bistability of nonlinear resonator
Due to the thermorefractive effect, both real and imaginary parts of the refractive index of the material may start to depend on temperature resulting in correspondent dispersive and absorptive Kerr-type nonlinearities. One of the main consequences is the appearance of a bistability regime  when the system may have several stationary states. In this section, we will discuss this behavior in more detail as it is directly related to nonlinear critical coupling conditions and optimized optical heating. It is worth mentioning, that we neglect the thermal expansion effect as it gives a much smaller contribution to the resonance shift in comparison with the thermooptical effect (see Table 2 ).
Within a single-mode approximation, the thermooptical Kerr nonlinearity can be accounted for in the first order of perturbation theory through the shift of frequency spectral position ω 0 → ω 0 − α|a|2 and varied nonradiative losses γ nr → γ nr + β|a|2. Here, we assume that both nonlinear coefficients are positive α, β > 0 which results in the redshift of the resonance and increase of losses with temperature increase. The exact expressions are presented in the following Section 2.4. With an account for this the equation for the stationary mode amplitude Eq. (3) will be modified:
Here, we introduce dimensionless parameters for the mode intensity y = α|a|2/γ, the external pump intensity y 0 = αγ r|f|2/γ 3, the relative nonlinear coefficient b = β/α, and normalized frequency detuning . One can notice that in the system there are only three independent parameters: the relative nonlinear coefficient b, normalized detuning frequency , and normalized pump intensity y 0. Interestingly, to analyze the pump intensity y 0 = Δω eff/γ which has the form of the ratio of the effective spectral shift of the resonance under the external pump Δω eff = α|f|2/γ 2 over the spectral width of the resonance γ. Thus, to increase y 0 one can increase the pump field intensity, decrease the total losses in the system, or increase the nonlinear coefficient making the system more sensitive to heating. Alternatively, the equation can be rewritten in a more compact form
Equations (4) and (5) is the central equation describing the state of the nonlinear resonator and its solution y describes the intensity of the field inside the resonator. However, this is a third-order algebraic equation and it may have not a unique solution but also three solutions that correspond to the bistability regime . Indeed, the dependence of the mode intensity y on the pump intensity y 0 is shown in Figure 2(a) for different values of the detuning parameter demonstrating highly nonlinear dependence. Moreover, for the detuning, one can observe the hysteresis behavior. The loop region corresponds to three different solutions and only two of them are stable at the upper and lower parts of the s-type curve  (see Figure 2(a)). For the larger values of detuning the hysteresis will be observed at larger pump intensities y 0 not shown in the plot.
Careful analysis of Eq. (4) shows that there exists a set of parameters that provide the bistability condition. To illustrate that, we plot a map showing the mode intensity y dependence on pump intensity y 0 and detuning (see Figure 2(b)) for fixed value of the relative nonlinear coefficient b = 0.05. The color shows the mode intensity normalized over the maximal value y max in Eq. (5) and discussed below. The hatched area denotes the region of bistability and the upper branch of the solution is depicted on the map. One can see the abrupt drop of the intensity at the right edge of the bistability region corresponding to the drop from the upper to lower branch of the solution (point B in the figure). The points A and B show the intensities corresponding to the switching between the upper and lower branches of the solution (see also Figure 2(a)).
Further analysis shows that the bistability region strongly depends on b, which is the ratio of the imaginary and real parts of thermorefractive constant. Critical point O which denotes the appearance of the bistability region moves upwards on the map with the increase of b as illustrated in Figure 2(c) where the bistability regions are shown for varied parameters b. It turns out that for b = 0, which corresponds to an absence of thermal dependence of ohmic losses, the bistability can appear at the smallest value of the intensity and detuning . With the increase of b the bistability region moves upwards and shrinks in the spectral width.
Such behavior has a clear physical explanation. Large b means that the losses in the system rapidly increase with the amplitude of mode. To maintain bistability one needs to decrease the losses in the “cold” system γ, which means to increase the pump intensity y 0. However, what is less obvious, that the bistability exists only for values of b smaller the . Indeed, once b → b cr the point O goes to infinity along the black dashed line in Figure 2(c). Thus, for a large enough absorptive nonlinear coefficient the bistability can not be observed. One can interpret the existence of the critical value b cr in the following manner: at large β/α, the spectral width of the resonance increases with temperature much faster than the resonant spectral position shifts, and the bistability condition (i.e. the shift of the peak should be larger than its width) simply can not be achieved. This does not allow achieving a bistability regime in a single-mode nonlinear resonator (for more details see Supplementary Information S1).
Finally, it is important to discuss the hysteresis loop parameters such as hysteresis width W h and height H h, which is the x-axis and y-axis distance between the turning points (A and B) in Figure 2(a), respectively. In the case of purely dispersive nonlinearity b = 0, the parameters can be immediately derived in a simple form (see details in Supplementary Information, Section S1):
One can see that the height and width of the hysteresis loop appear to be dependent on the detuning frequency and are mutually dependent. For large frequency detuning values , these quantities behave as and , whereas for close to the critical point values they have the following dependencies on the detuning: H h ∼ δ 1/2 and W h ∼ δ 3/2. Consequently, for small values of detuning near the critical point hysteresis height increases more rapidly than its width. Thus, for potential optical switching applications, it is more prospective to work closer to the critical frequency , where the height of the hysteresis loop is high, while the width is small providing a stronger amplitude difference between the stationary states at lower switching intensities. The graphical plots of these parameters along with additional details on b > 0 case are provided in Supplementary Information materials, Section, Section S1.
2.3 Nonlinear critical coupling
Till now, we have discussed the possible states of the resonator with thermally induced Kerr nonlinearity. The maximal possible mode intensity is provided by the nonlinear critical coupling condition based on the straightforward analysis of Eq. (4). The Lorentz-type of right-hand side of Eq. (5) ensures that the value y can not be greater than y ≤ y max. Moreover, it turns out that there are a certain set of parameters for which the maximal value y max is reached and they are given by the equation
This equation defines the nonlinear critical coupling conditions maximizing the mode intensity y. In the first-order perturbation, one can obtain the proper detuning frequency provides the critical coupling. This condition means that to reach the maximal mode intensity one should have particular detuning between the pump frequency and the “cold” resonant frequency. The stronger the pump intensity the larger should be the detuning. Now, once the b-coefficient becomes stronger the higher-order dependence appears providing that .
The critical coupling regime is seen in Figure 2(b) where the normalized mode intensity is plotted as the function of and y 0. One can see that the maximal value is reached along the black dashed line, which is obtained as the solution of the nonlinear critical coupling equation Eq. (7). One can see that the critical coupling is reached close to the edge of the bistability region at the upper branch of the solution. So, to achieve critical coupling one needs to drive the system in the bistable state at the upper branch.
The nonlinear critical coupling condition formulated in Eq. (7) in dimensionless units can be reformulated in terms of radiative and nonradiative losses, which is convenient for the designing the optical resonators. The total absorption in the regime of maximal heating when y = y max will have a for similar to the classical one with slight modification
Now for any fixed value of γ r the maximal absorption will be observed at
which is an analog of the classical critical coupling condition and has a clear physical meaning, at least in the first order perturbation with respect to b: (i) the system in its “cold” state should be out of the linear critical coupling condition; (ii) once the resonator is pumped with the particular intensity and at spectral detuning in accordance to Eq. (9), in the final “hot” stationary state the radiative and nonradiative losses should be balanced. Indeed, expanding Eq. (7) in the series with respect to b ≪ 1, one obtains that and Thus, the second term in Eq. (9) simply corresponds to the added losses due to the heating of the resonator.
2.4 Nonlinear coefficients α and β
The behavior of the nonlinear resonator strongly depends on the values of dispersive and absorptive nonlinear coefficients α and β which provide the nonlinear coupling. Their amplitude is fully defined by the resonant mode properties and thermorefractive coefficients expressing the thermal origin of nonlinear coupling. One can derive it starting from very general considerations of the absorbed power:
where ɛ 0 is the permittivity of vacuum, ɛ″ is the imaginary part of the complex dielectric permittivity ɛ = ɛ′ + iɛ″, ω is the excitation frequency and E(r) is the electric field amplitude inside of the nanoparticle, and the integration is taken over the volume of the particle. We can introduce electric field through the mode amplitude as , where M(r) is the normalized eigenmode field distribution . For high-Q modes this normalization can be approximated by the expression . The thermorefractive effect in the linear approximation can be introduced as follows:
where n 0 and k 0 are the initial values of complex refractive index at room temperature T = 298 K, n 1 = dn/dT| T=298K and k 1 = dk/dT| T=298K are the linear thermorefracitve coefficients. The imaginary part of the dielectric permittivity also increases linearly with temperature and in the first order approximation we have then and expression (10) could be rewritten as:
Finally, we relate the temperature with the absorbed electromagnetic power ΔT = κP abs through the linear coefficient κ similarly to Eq. (1). This coefficient depends on the particular geometry of the system and can be found with help of exact numerical simulation, for instance. With this, the temperature can be expressed from Eq. (12)
where is proportional to the mode intensity, and is critical temperature increase. One can see that the temperature increase by itself has a nonlinear relation with the mode amplitude. However, as far as the temperature increase is smaller than the critical one the connection can be linearized.
For a dielectric resonator the frequency shift is governed by the change of the refractive index, thus providing Δω 0/ω 0 = −Δn/n. Recalling introduced mode spectral shift from Section 2.2 and temperature dependence of the complex refractive index (11), we obtain
where the last approximation is made for small values of frequency detuning Δω/ω 0 ≪ 1.
Similarly, we can derive the absorptive nonlinear coefficient β. From the amplitude nonradiative losses rate definition:
Using expression (14) for the relation between temperature increase inside of the nanoparticle and mode amplitude, we obtain the nonlinear coefficient β:
which is true if we work in the vicinity of an unperturbed eigenfrequency position. Surprisingly, the relative thermooptical parameter b, which governs the general behavior of the nonlinear resonator, has the following form:
which does not depend on the resonant characteristics of the nanoparticles but only on the material properties and thermorefractive characteristics.
3 Quasi-BIC nanoresonator heating
In this section, we implement the developed approach for optimizing the optical heating of a single nanoresonator. In strong contrast to nanoplasmonic designs, we aim at semiconductor materials, which demonstrate strong thermooptical nonlinear coefficient  and also can demonstrate high-Q resonance which provides strong thermooptical reconfiguration under moderate laser intensities. Moreover, the semiconductor structures allow for precise tuning of losses in a wide spectral range by doping the material with free carriers. Together with the fine engineering of radiative losses with the concept of quasi-BIC modes [3, 38] one can achieve high efficiency of optical heating of a single wavelength scale resonator.
3.1 Radiative and nonradiative losses optimization
Quasi-BIC states are characterized by the destructive interference of radiation in the far-field zone and, hence, efficient localization of electromagnetic energy in the resonator volume . The high-Q states in nanocylinders can be excited with polarized vector beam  enhancing light–matter interaction. Full-wave numerical simulations allowed us to locate the spectral position of high-Q modes silicon nanocylinders on quartz substrates (see Figure 3(a) inset). The modes of cylindrical resonators are classified by their azimuthal number m , which indicates the symmetry of the electromagnetic field E, H ∼ eimφ , where φ is the azimuthal angle to cylinder axes. We investigate the azimuthally symmetrical (m = 0) eigenmodes matching the symmetry of the incident azimuthal vector beam (see Figure 1). The map showing the eigenmodes spectral position as a function of the aspect ratio of the cylinder is shown in Figure 3(a). One can see the formation of a high-Q state due to coupling between two modes in the vicinity of the anticrossing region. Fixing then the incident wavelength λ = 1400 nm we obtain geometrical parameters for the declared quasi-BIC state which are radius ρ = 1238 nm and height h = 885 nm.
Next, we have optimized the nonradiative losses to get close to the critical coupling condition. The doping of semiconductors provides an additional degree of freedom for precise nonradiative losses control. We choose silicon for the resonator material since it possesses zero optical losses in the near-IR region  which gives an opportunity to finely tune the losses in the wide range via modulation of dopant concentration.
For chosen geometry of the resonator, we perform rigorous optimization for nonradiative losses by simulating in numerical commercial software COMSOL Multiphysics temperature increase inside of the nanoparticle excited by the azimuthal vector beam  Figure 3(b). Analysis of the temperature as a function of the imaginary part of the refractive index k responsible for the nonradiative losses rate in the system allows identifying the critical coupling value of optical losses k = 0.0031. The estimated concentration of dopants providing this condition according to the Drude model is expressed as:
where ω is the frequency of the incident radiation; τ e = 1 fs is the electron momenta relaxation time; plasmonic frequency ; m eff = 0.18m e is the effective mass of electrons in the conduction band of c-Si ; ɛ 0 is the permittivity of vacuum, and e is the elementary charge. Assuming that the real part of the refractive index n does not change significantly with the free carrier concentration increase , we immediately obtain donors concentration n d = 4.6 ⋅ 1018 cm−3 corresponding to the optimal value of k = 0.0031 Extracted total loss rate for the resonator with optimized parameters is then γ = γ r + γ nr = 2.4 ⋅ 1012 1/s.
3.2 Numerical modeling of optical heating
We have performed the simulations on the optical heating of the designed nanoresonators demonstrating a huge temperature increase ΔT = 815 K for relatively low incident power flux I = 0.55 mW/μm2 in the linear regime when the thermooptical effects are omitted. However, once the nonlinearity of the system is taken into account, heating is significantly suppressed under the resonant excitation (compare linear and nonlinear regimes in Figure 3(c)). Such a noticeable deviation from the linear trend is associated with the eigenfrequency spectral shift shown in Figure 3(d). Thus, for maximizing heating efficiency in a nonlinear thermorefractive system one needs to obtain the optimal condition concerning both excitation frequency and incident intensity by the basic theory described in Section 2.
In our design of the nanoresonator tuned for IR-region, the relative absorptive nonlinear coefficient b appears to be negligibly low b ≈ 0. Indeed, in the near-IR region crystalline silicon has zero intrinsic optical losses and therefore nonradiative losses in the system are only defined by the concentration of dopants. Since typical values of donors ionization energy are less than thermal energy for the room temperature E ∼ kT, we suppose them to be fully ionized at T = 298 K . Consequently, temperature increase does not influence optical losses and imaginary part of the thermorefractive coefficient can be set to zero k 1 = dk/dT| T=298K = 0. At the same time, the real part of the thermooptical coefficient at the wavelength λ = 1400 nm equals to n 1 = 2 ⋅ 10−4 1/K  which from Eqs. (15)–(17) results into α = 3.25 ⋅ 1028 J−1 s−1 and β = 1.44 ⋅ 1025 J−1 s−1 nonlinear coefficients.
With the given material parameters, we perform full-wave simulations coupled with the heat transfer module in commercial software COMSOL Multiphysics. The thermooptical coupling provides the nonlinear response of the simulated system which may initiate computational difficulties once the bistability regime is reached. In this case, the final state of the iterative numerical solution depends on the initial solutions guess whether the solution is located on the upper, lower, or unstable part of the s-curve shown in Figure 4(a). For that, we were carefully choosing the initial guess and used the obtained solutions as the initial guess for the next set of parameters. The dependence of the heating efficiency denoting the temperature increase per unit incident power ΔT/I 0 is shown in Figure 4(b) and (c) in full analog with the plots shown in Figure 2(b) for a toy model. Figure 4(b) and (c) shows the heating efficiency at the upper and lower branch of the solution correspondingly. One can see that the maximal efficiency of heating is indeed reached at the upper branch by the basic theory discussed in Section 2.2. To reach this solution in the simulations, one can simply start with the high intensity of the laser pump and high temperature initial guess and then gradually decrease the pump power. However, this is not very physical from the experimental point of view, when normally the excitation of the systems starts from the cold regime. In this view, we show several “trajectories” in the parameters space (see Figure 2), which are tested in numerical simulations. The trajectory I stands for the increase of the incident power at the frequency detuning where there is no bistability region. For the trajectory II the detuning frequency is larger and with the increase of the intensity one starts at the lower branch of the solution (see Figure 2) and then by a gradual increase of the pump intensity reaches the upper branch, after that with the decreasing of the temperature one can reach the maximal heating efficiency point A. However, such a route leads to necessary overheating of the structure: reaching the optimal heating point one first need to jump on the upper steady-state branch at higher temperature and then get to the critical point A (see Figure 2(b) and trajectory II there). Thus, alternatively, one can choose the trajectory II′ where one first increases the pump at small detuning and then gradually increases the detuning at constant power reaching again point A. Finally, the shown trajectory III does not reach the supper branch and the heating efficiency is kept at a very low level.
It is also worth noting that at the point of maximal heating efficiency in the nonlinear regime, its value becomes equal to the heating efficiency in the linear regime (see the linear dash-dotted line in Figure 4(a)), which also matches the basic model described in Section 2. However, that occurs only in the case of small dissipative nonlinearity b ≈ 0.
4 Discussion and conclusion
The main idea of this paper is to reveal the key aspects of optical heating of a semiconductor resonator: (i) once the resonator is heated in the CW regime, the efficiency of heating immediately drops due to thermooptical reconfiguration of the resonator; (ii) one can reach high efficiency only the heating is accompanied with proper spectral detuning of the pumping laser from the ‘cold’ resonance of the mode; (iii) The maximal heating efficiency is reached in the bistability regime which appears under certain conditions on the pumping intensity, thermooptical coefficients, and of the resonant mode.
However, despite the discussed design being aimed at the near-IR range, where the intrinsic losses of silicon are vanishing, one can tune the proposed design closer to the visible spectra and balance the radiative losses with the intrinsic losses of silicon, which occur at around λ = 865 nm. Here, the imaginary part of thermorefractive coefficient is non-zero k 1 = 8.5 ⋅ 10−6 1/K , which provides the relative nonlinear coefficient b = 0.015.
It is also worth mentioning that the suggested effect of nonlinear critical coupling can be observed not only in the CW regime, but also under the pulsed excitation once the pulse duration is long enough so that the equilibrium temperature is achieved, i.e. for nanosecond laser pulses. In that case, the quasi-CW regime can be considered and chirped laser pulses could be used to achieve the efficient heating of nanoresonators. Alternatively, the heating effects under the short pico- and femto-second pulses excitation requires more complex models based on the analysis of nonequilibrium dynamics of carriers .
Finally, we would like to provide a brief comparison of the proposed nanophotonic design in terms of heating efficiency to existing analogs of nanoscale optical heaters. Indeed, Figure 4 depicts the map of heating efficiencies as a function of pump intensity and the detuning factor. The maximum values of heating efficiencies reached 1400 K/(mW/μm2), which is the highest value of heating efficiencies for the structures with direct thermal contact with substrates. For a single c-Si nanodisk laser heating on a substrate at magnetic dipole and quadrupole modes the efficiency reaches 150 K/(mW/μm2) , and for a c-Si sphere is about 300 K/(mW/μm2)  with considerably less thermal contact with the substrate than the disk has. The latter results were obtained experimentally, however, the best theoretical values reach 750 K/(mW/μm2) for complex semiconductor structures with quasi-infinite c-Si nanorod covered with a-Si film  and one of the most promising nanostructures for light-to-heat conversion is golden doughnut supporting efficiency values up to 230 K/(mW/μm2) in aqueous media. Nevertheless, it remains questionable of fabrication, experimental feasibility, and temperature detection possibility . On the other hand, thermal nonlinearity driven by pulsed laser heating is also a rather promising approach, but the experimentally demonstrated efficiency for c-Si nanocubes reached only 30 K/(mW/μm2)  mostly due to oil immersion, which drastically increases the thermal conductivity of the nanoresonator’s surrounding medium.
In conclusion, we have developed a new simple theoretical approach to optimization of the resonator optical heating in the nonlinear regime. The proposed design based on the super-cavity mode in doped silicon cylindrical particles has allowed for efficient light-to-heat conversion when the initial spectrum of incident laser is detuned from the initial spectral position of the resonance. Moreover, we have revealed a bistability regime in the optical heating at intensity around 1 mW/μm2. Our results are also helpful for resolving the thermal challenges for all-dielectric resonator-based photonic devices [49, 50], Raman microlasers [51, 52], and nanoscale photo-thermal chemistry and sensing [53–55]. As an outlook, we believe that the developed bistability approach is quite universal and can be further applied not only for various types of nonlinearities based on Kerr effect [56, 57], electron-hole plasma generation [42, 43, 58], and excitonic effects [59–61].
Award Identifier / Grant number: 075-15-2021-589
Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
Acknowledgments: The authors thank Kirill Koshelev, Zarina Sadrieva, Alexey Yulin for fruitful discussions. The research was supported by Priority 2030 Federal Academic Leadership Program. The analytical simulations were supported by Russian Scientific Foundation (proj. 18-72-10140). Numerical simulations were supported by the Ministry of Science and Higher Education of the Russian Federation (proj. 075-15-2021-589).
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
 A. I. Kuznetsov, A. E. Miroshnichenko, M. L. Brongersma, Y. S. Kivshar, and B. Luk’yanchuk, “Optically resonant dielectric nanostructures,” Science, vol. 354, no. 6314, p. aag2472, 2016. https://doi.org/10.1126/science.aag2472.Search in Google Scholar PubMed
 K. Koshelev, G. Favraud, A. Bogdanov, Y. Kivshar, and A. Fratalocchi, “Nonradiating photonics with resonant dielectric nanostructures,” Nanophotonics, vol. 8, no. 5, pp. 725–745, 2019. https://doi.org/10.1515/nanoph-2019-0024.Search in Google Scholar
 K. Koshelev, S. Kruk, E. Melik-Gaykazyan, et al.., “Subwavelength dielectric resonators for nonlinear nanophotonics,” Science, vol. 367, no. 6475, pp. 288–292, 2020. https://doi.org/10.1126/science.aaz3985.Search in Google Scholar PubMed
 E. N. Bulgakov and A. F. Sadreev, “High- Q resonant modes in a finite array of dielectric particles,” Phys. Rev., vol. 99, no. 3, pp. 1–9, 2019. https://doi.org/10.1103/physreva.99.033851.Search in Google Scholar
 E. N. Bulgakov and D. N. Maksimov, “Topological bound states in the continuum in arrays of dielectric spheres,” Phys. Rev. Lett., vol. 118, no. 26, pp. 1–5, 2017. https://doi.org/10.1103/PhysRevLett.118.267401.Search in Google Scholar PubMed
 D. F. Kornovan, R. S. Savelev, Y. Kivshar, and M. I. Petrov, “High-Q localized states in finite arrays of subwavelength resonators,” ACS Photonics, vol. 8, no. 12, pp. 3627–3632, 2021. https://doi.org/10.1021/acsphotonics.1c01262.Search in Google Scholar
 S. I. Azzam and A. V. Kildishev, “Photonic bound states in the continuum: from basics to applications,” Adv. Opt. Mater., vol. 9, no. 1, pp. 16–24, 2021. https://doi.org/10.1002/adom.202001469.Search in Google Scholar
 C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos, and M. Soljačić, “Bound states in the continuum,” Nat. Rev. Mater., vol. 1, no. 9, pp. 1–13, 2016. https://doi.org/10.1038/natrevmats.2016.48.Search in Google Scholar
 S. V. Makarov, A. S. Zalogina, M. Tajik, et al.., “Light-induced tuning and reconfiguration of nanophotonic structures,” Laser Photon. Rev., vol. 11, no. 5, p. 1700108, 2017. https://doi.org/10.1002/lpor.201700108.Search in Google Scholar
 G. P. Zograf, M. I. Petrov, S. V. Makarov, and Y. S. Kivshar, “All-dielectric thermonanophotonics,” Adv. Opt Photon, vol. 13, p. 643, 2021. https://doi.org/10.1364/aop.426047.Search in Google Scholar
 G. Baffou and R. Quidant, “Thermo-plasmonics: using metallic nanostructures as nano-sources of heat,” Laser Photon. Rev., vol. 7, no. 2, pp. 171–187, 2013. https://doi.org/10.1002/lpor.201200003.Search in Google Scholar
 G. Baffou, F. Cichos, and R. Quidant, “Applications and challenges of thermoplasmonics,” Nat. Mater., vol. 19, no. 9, pp. 946–958, 2020. https://doi.org/10.1038/s41563-020-0740-6.Search in Google Scholar PubMed
 S. Danesi, M. Gandolfi, L. Carletti, et al.., “Photo-induced heat generation in non-plasmonic nanoantennas,” Phys. Chem. Chem. Phys., vol. 20, no. 22, pp. 15307–15315, 2018. https://doi.org/10.1039/c8cp01919c.Search in Google Scholar PubMed
 D. Rocco, M. Gandolfi, A. Tognazzi, et al.., “Opto-thermally controlled beam steering in nonlinear all-dielectric metastructures,” Opt Express, vol. 29, no. 23, pp. 37128–37139, 2021. https://doi.org/10.1364/oe.440564.Search in Google Scholar
 O. Pashina, K. Frizyuk, G. Zograf, and M. Petrov, “Thermo-optical reshaping of second-harmonic emission from dimer all-dielectric nanoresonators,” Opt. Lett., vol. 47, pp. 1992–1995, 2022. https://doi.org/10.1364/OL.444348.Search in Google Scholar PubMed
 G. P. Zograf, M. I. Petrov, D. A. Zuev, et al.., “Resonant nonplasmonic nanoparticles for efficient temperature-feedback optical heating,” Nano Lett., vol. 17, no. 5, pp. 2945–2952, 2017. https://doi.org/10.1021/acs.nanolett.7b00183.Search in Google Scholar PubMed
 M. Aouassa, E. Mitsai, S. Syubaev, et al.., “Temperature-feedback direct laser reshaping of silicon nanostructures,” Appl. Phys. Lett., vol. 111, p. 243103, 2017. https://doi.org/10.1063/1.5007277.Search in Google Scholar
 M. Celebrano, D. Rocco, M. Gandolfi, et al.., “Optical tuning of dielectric nanoantennas for thermo-optically reconfigurable nonlinear metasurfaces,” Opt Lett., vol. 46, no. 10, p. 2453, 2021. https://doi.org/10.1364/ol.420790.Search in Google Scholar
 X. Jiang and L. Yang, “Optothermal dynamics in whispering-gallery microresonators,” Light Sci. Appl., vol. 9, pp. 1–15, 2020. https://doi.org/10.1038/s41377-019-0239-6.Search in Google Scholar PubMed PubMed Central
 Yu-L. Tang, Te-H. Yen, K. Nishida, et al.., “Mie-enhanced photothermal/thermo-optical nonlinearity and applications on all-optical switch and super-resolution imaging [Invited],” Opt. Mater. Express, vol. 11, no. 11, p. 3608, 2021. https://doi.org/10.1364/ome.431533.Search in Google Scholar
 C. H. Li, Y. L. Tang, J. Takahara, and S. W. Chu, “Nonlinear heating and scattering in a single crystalline silicon nanostructure,” J. Chem. Phys., vol. 155, no. 20, p. 204202, 2021.10.1063/5.0067251Search in Google Scholar PubMed
 G. J. Huang, H. Y. Cheng, Y. L. Tang, et al.., “Transient super-/sub-linear nonlinearities in silicon nanostructures,” Adv. Opt. Mater., vol. 10, no. 5, pp. 1–10, 2022. https://doi.org/10.1002/adom.202101711.Search in Google Scholar
 Y. S. Duh, Y. Nagasaki, Y. L. Tang, et al.., “Giant photothermal nonlinearity in a single silicon nanostructure,” Nat. Commun., vol. 11, no. 1, pp. 1–9, 2020. https://doi.org/10.1038/s41467-020-17846-6.Search in Google Scholar PubMed PubMed Central
 T. Zhang, Y. Che, K. Chen, et al.., “Anapole mediated giant photothermal nonlinearity in nanostructured silicon,” Nat. Commun., vol. 11, no. 1, pp. 1–9, 2020. https://doi.org/10.1038/s41467-020-16845-x.Search in Google Scholar PubMed PubMed Central
 V. Grigoriev, N. Bonod, J. Wenger, and B. Stout, “Optimizing nanoparticle designs for ideal absorption of light,” ACS Photonics, vol. 2, no. 2, pp. 263–270, 2015. https://doi.org/10.1021/ph500456w.Search in Google Scholar
 A. E. Miroshnichenko and M. I. Tribelsky, “Ultimate absorption in light scattering by a finite obstacle,” Phys. Rev. Lett., vol. 120, no. 3, p. 033902, 2018. https://doi.org/10.1103/PhysRevLett.120.033902.Search in Google Scholar PubMed
 J. Cunha, T.-L. Guo, G. D. Valle, A. N. Koya, R. P. Zaccaria, and A. Alabastri, “Controlling light, heat, and vibrations in plasmonics and phononics,” Adv. Opt. Mater., vol. 225, p. 2020, 2001.10.1002/adom.202001225Search in Google Scholar
 K. N. Reddy, A. V. Gopal, and S. D. Gupta, “Nonlinearity induced critical coupling,” Opt Lett., vol. 38, no. 14, p. 2517, 2013. https://doi.org/10.1364/ol.38.002517.Search in Google Scholar PubMed
 S. Longhi, “Coherent perfect absorption in a homogeneously broadened two-level medium,” Phys. Rev., vol. 83, no. 5, pp. 1–4, 2011. https://doi.org/10.1103/physreva.83.055804.Search in Google Scholar
 J. T. Shen and S. Fan, “Quantum critical coupling conditions for zero single-photon transmission through a coupled atom-resonator-waveguide system,” Phys. Rev., vol. 82, no. 2, pp. 4–7, 2010. https://doi.org/10.1103/physreva.82.021802.Search in Google Scholar
 W. Suh, Z. Wang, and S. Fan, “Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities,” IEEE J. Quant. Electron., vol. 40, no. 10, pp. 1511–1518, 2004. https://doi.org/10.1109/jqe.2004.834773.Search in Google Scholar
 H. Haus, Waves and Fields in Optoelectronics, vol. 402, Englewood Cliffs, NJ 07632, USA, Prentice-Hall, 1984.Search in Google Scholar
 M. B. Doost, W. Langbein, and E. A. Muljarov, “Resonant-state expansion applied to three-dimensional open optical systems,” Phys. Rev., vol. 90, no. 1, p. 013834, 2014. https://doi.org/10.1103/physreva.90.013834.Search in Google Scholar
 A. A. Bogdanov, K. L. Koshelev, P. V. Kapitanova, et al.., “Bound states in the continuum and fano resonances in the strong mode coupling regime,” Adv. Photon., vol. 1, no. 1, p. 016001, 2019. https://doi.org/10.1117/1.ap.1.1.016001.Search in Google Scholar
 S. Gladyshev, K. Frizyuk, and A. Bogdanov, “Symmetry analysis and multipole classification of eigenmodes in electromagnetic resonators for engineering their optical properties,” Phys. Rev. B, vol. 102, no. 7, p. 75103, 2020. https://doi.org/10.1103/physrevb.102.075103.Search in Google Scholar
 M. A. Green, “Self-consistent optical parameters of intrinsic silicon at 300 k including temperature coefficients,” Sol. Energy Mater. Sol. Cell., vol. 92, no. 11, pp. 1305–1310, 2008. https://doi.org/10.1016/j.solmat.2008.06.009.Search in Google Scholar
 K. Sokolowski-Tinten and D. von der Linde, “Generation of dense electron-hole plasmas in silicon,” Phys. Rev. B, vol. 61, no. 4, p. 2643, 2000. https://doi.org/10.1103/physrevb.61.2643.Search in Google Scholar
 S. Makarov, S. Kudryashov, I. Mukhin, et al.., “Tuning of magnetic optical response in a dielectric nanoparticle by ultrafast photoexcitation of dense electron–hole plasma,” Nano Lett., vol. 15, no. 9, pp. 6187–6192, 2015. https://doi.org/10.1021/acs.nanolett.5b02534.Search in Google Scholar PubMed
 B. J. Frey, D. B. Leviton, and T. J. Madison, “Temperature-dependent refractive index of silicon and germanium,” in Optomechanical Technologies for Astronomy, vol. 6273, International Society for Optics and Photonics, 2006, p. 62732J.10.1117/12.672850Search in Google Scholar
 G. E. JellisonJr and F. A. Modine, “Optical functions of silicon at elevated temperatures,” J. Appl. Phys., vol. 76, no. 6, pp. 3758–3761, 1994. https://doi.org/10.1063/1.357378.Search in Google Scholar
 G. P. Zograf, Ye. F. Yu, K. V. Baryshnikova, A. I. Kuznetsov, and S. V. Makarov, “Local crystallization of a resonant amorphous silicon nanoparticle for the implementation of optical nanothermometry,” JETP Lett. (Engl. Transl.), vol. 107, no. 11, pp. 699–704, 2018. https://doi.org/10.1134/s0021364018110140.Search in Google Scholar
 J. González-Colsa, G. Serrera, J. María Saiz, et al.., “Gold nanodoughnut as an outstanding nanoheater for photothermal applications,” Opt Express, vol. 30, no. 1, pp. 125–137, 2022. https://doi.org/10.1364/oe.446637.Search in Google Scholar PubMed
 K. Padmaraju and K. Bergman, “Resolving the thermal challenges for silicon microring resonator devices,” Nanophotonics, vol. 3, nos 4–5, pp. 269–281, 2014. https://doi.org/10.1515/nanoph-2013-0013.Search in Google Scholar
 S. Wang, Y. Li, B. E. Little, et al.., “Athermal third harmonic generation in micro-ring resonators,” Opto-Electron. Adv., vol. 3, no. 12, p. 12200028, 2020. https://doi.org/10.29026/oea.2020.200028.Search in Google Scholar
 D. Agarwal, M.-L. Ren, J. S. Berger, J. Yoo, A. Pan, and R. Agarwal, “Nanocavity-enhanced giant stimulated Raman scattering in si nanowires in the visible light region,” Nano Lett., vol. 19, no. 2, pp. 1204–1209, 2019. https://doi.org/10.1021/acs.nanolett.8b04666.Search in Google Scholar PubMed
 G. P. Zograf, D. Ryabov, V. Rutckaia, et al.., “Stimulated Raman scattering from mie-resonant subwavelength nanoparticles,” Nano Lett., vol. 20, no. 8, pp. 5786–5791, 2020. https://doi.org/10.1021/acs.nanolett.0c01646.Search in Google Scholar PubMed
 M. Caldarola, P. Albella, E. Cortés, et al.., “Non-plasmonic nanoantennas for surface enhanced spectroscopies with ultra-low heat conversion,” Nat. Commun., vol. 6, no. 1, pp. 1–8, 2015. https://doi.org/10.1038/ncomms8915.Search in Google Scholar PubMed PubMed Central
 R. Regmi, J. Berthelot, P. M. Winkler, et al.., “All-dielectric silicon nanogap antennas to enhance the fluorescence of single molecules,” Nano Lett., vol. 16, no. 8, pp. 5143–5151, 2016. https://doi.org/10.1021/acs.nanolett.6b02076.Search in Google Scholar PubMed
 V. A. Milichko, D. A. Zuev, D. G. Baranov, et al.., “Metal-dielectric nanocavity for real-time tracing molecular events with temperature feedback,” Laser Photon. Rev., vol. 12, no. 1, p. 1700227, 2018. https://doi.org/10.1002/lpor.201700227.Search in Google Scholar
 M. R. Shcherbakov, P. P. Vabishchevich, A. S. Shorokhov, et al.., “Ultrafast all-optical switching with magnetic resonances in nonlinear dielectric nanostructures,” Nano Lett., vol. 15, no. 10, pp. 6985–6990, 2015. https://doi.org/10.1021/acs.nanolett.5b02989.Search in Google Scholar PubMed
 G. Grinblat, H. Zhang, M. P. Nielsen, et al.., “Efficient ultrafast all-optical modulation in a nonlinear crystalline gallium phosphide nanodisk at the anapole excitation,” Sci. Adv., vol. 6, no. 34, p. eabb3123, 2020. https://doi.org/10.1126/sciadv.abb3123.Search in Google Scholar PubMed PubMed Central
 I. S. Sinev, K. Koshelev, Z. Liu, et al.., “Observation of ultrafast self-action effects in quasi-bic resonant metasurfaces,” Nano Lett., vol. 21, no. 20, pp. 8848–8855, 2021. https://doi.org/10.1021/acs.nanolett.1c03257.Search in Google Scholar PubMed
 H. M. Gibbs, S. L. McCall, T. N. C. Venkatesan, A. C. Gossard, A. Passner, and W. Wiegmann, “Optical bistability in semiconductors,” Appl. Phys. Lett., vol. 35, no. 6, pp. 451–453, 1979. https://doi.org/10.1063/1.91157.Search in Google Scholar
 H. M. Gibbs, S. S. Tarng, J. L. Jewell, et al.., “Room-temperature excitonic optical bistability in a gaas-gaalas superlattice étalon,” Appl. Phys. Lett., vol. 41, no. 3, pp. 221–222, 1982. https://doi.org/10.1063/1.93490.Search in Google Scholar
 M. A. Masharin, V. A. Shahnazaryan, F. A. Benimetsky, et al.., “Polaron-enhanced polariton nonlinearity in lead halide perovskites,” arXiv preprint arXiv:2201.10265, 2022.10.1021/acs.nanolett.2c03524Search in Google Scholar PubMed
The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph-2022-0074).
© 2022 the author(s), published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.