Bound states in the continuum and Fano resonances in the strong mode coupling regime

The study of resonant dielectric nanostructures with high refractive index is a new research direction in nanoscale optics and metamaterial-inspired nanophotonics. Because of the unique optically-induced electric and magnetic Mie resonances, high-index nanoscale structures are expected to complement or even replace different plasmonic components in a range of potential applications. Here we study strong coupling between modes of a single subwavelength high-index dielectric resonator and analyse the mode transformation and Fano resonances when resonator's aspect ratio varies. We demonstrate that strong mode coupling results in resonances with high quality factors, which are related to the physics of bound states in the continuum when the radiative losses are almost suppressed due to the Friedrich-Wintgen scenario of destructive interference. We explain the physics of these states in terms of multipole decomposition and show that their appearance is accompanied by drastic change of the far-field radiation pattern. We reveal a fundamental link between the formation of the high-quality resonances and peculiarities of the Fano parameter in the scattering cross-section spectra. Our theoretical findings are confirmed by microwave experiments for the scattering of a high-index cylindrical resonators with a tunable aspect ratio. The proposed mechanism of the strong mode coupling in single subwavelength high-index resonators accompanied by resonances with high quality factor helps to extend substantially functionalities of all-dielectric nanophotonics that opens new horizons for active and passive nanoscale metadevices.


I. INTRODUCTION
The physics of resonant structures with strong mode coupling is of the fundamental importance being responsible for a variety of interesting phenomena governing both transport and localization of waves. The modes supported by traditional resonators and microcavities [1] exist due to reflection of waves from the resonator's boundaries under the conditions of constructive interference. To achieve high values of resonator's quality factor (Q factor), one can improve reflectivity by using metals [2,3], photonic bandgap structures [4], or the total internal reflection at glancing angles of incidence in whispering-gallery-mode (WGM) resonators [5]. All those physical mechanisms require large sizes of cavities with a complex design. A more attractive way to confine light is to use destructive interference in the regime of strong mode coupling [6][7][8]. This mechanism is related to the physics of bound states in the continuum (BIC) [9]. It was first proposed in quantum mechanics by Friedrich and Wintgen [10] and then was extended to acoustics [11][12][13] and electrodynamics [14,15]. A true optical BIC is a mathematical abstraction since its realization demands either infinite size of the structure or zero (or infinite) permittivity [16][17][18]. However, the BIC-inspired mechanism of light localization makes possible realization of high-Q states in photonic crystal cavities and slabs [15,17,19], coupled waveguide arrays [20][21][22], dielectric gratings [14], core-shell spherical * ysk@internode.on.net particles [18], dielectric resonators [23][24][25][26], and hybrid plasmonic-photonic systems [27].
By definition, the Q factor of a true BIC is infinite. Hence, in the wave scattering, BIC manifests itself as "a collapse" of the Fano resonance when the width of the resonance vanishes and the Fano feature disappears from the scattering spectrum [28,29]. In practice, both Q factor and the width of the Fano resonance at the frequency of BIC remain finite because of absorption, finite-size samples, roughness, and other imperfections [30]. Remarkably, in terms of the scattering matrix, BIC corresponds to merging of a pole and a zero of the scattering operator on the real axis [31]. Properties of the Fano resonances in the systems with BICs have been considered in several studies [15,17,18].
A conventional device supporting light localization via BIC-inspired mechanism is based on periodic photonic structures [17] or chains of scatterers [32]. For these structures, strong localization can be achieved only for a large number of scatterers because it is governed by their mutual interference. Other implementations of BICs in photonic structures were presented in Refs. [18,25]. In the former work, the design is based on a metallic box with a single dielectric scatterer inside. While the conditions of true BICs can be achieved here by imitating periodic boundary conditions for such a unit cell, the structure is not subwavelength, and it demands using metallic components. The latter work describes a single scatterer but it employs near-zero refraction index constituents which requires large scales of the structure [33], especially, at optical frequencies. Recently BIC-inspired supercavity modes in subwavelength dielectric resonators  Frequency ωr/c without singular permittivity values were proposed theoretically [8]. However, the coexisting of Fano resonance being effect of weak coupling and strong mode coupling underlying BIC is a vital problem of modern photonics.

A A A A A A A A A A A A A A A A A A A A A A A A A A A
In this paper we demonstrate, both theoretically and experimentally, strong light localization and existence of quasi-BICs in the simplest object -a single homogeneous cylindrical subwavelength dielectric resonator in free space. We show that quasi-BICs appear in accord with the Friedrich-Wintgen interference mechanism because of strong coupling between Mie-like and Fabry-Perot-like modes. We develop a novel analytical approach to describe light scattering by finite-size dielectric resonators and reveal close relationships between the peculiarities of quasi-BICs and the critical behaviour of Fano asymmetry parameter in the strong coupling regime. We show that the Fano asymmetry parameter becomes singular at the frequency of the quasi-BIC, and it vanishes when the mode becomes almost dark for far field excitation. We derive an exact form of coupling coefficients between modes and corresponding Rabi frequencies. We analyse effects of material losses and reveal that the mode coupling remains strong even for highly absorptive resonators. Our findings make evident the counterintuitive fact that even a subwavelength dielectric resonator could provide strong light localization.

A. Interplay of Mie and Fabry-Perot modes
We start our study with numerical simulations of the scattering cross-section (SCS) σ of a finite dielectric cylinder with permittivity ε 1 = 80, radius r, and length l placed in vacuum (ε 2 = 1), as shown in Fig. 1a. The n = 0,1 p ≠ 0 p = 0,1,.. spectra are calculated by using the CST Microwave Studio software and T-matrix computations [34,35]. The electric field of the incident wave is assumed to be perpendicular to the axis of the cylinder (see Fig. 1a). To compare σ for cylinders with different aspect ratios, we normalize σ by the projected cross-section of the resonator, S = 2rl. The maps of the normalized SCS σ/S calculated for cylinders with different aspect ratio r/l excited by TM and TE-polarized wave are shown in Figs. 1c and 1d, respectively. According to the standard nomenclature (see, e.g., Ref. [36]), we denote the modes of a cylindrical resonator as TE n,k,p and TM n,k,p , where n, k, p are the indices denoting the azimuthal, radial, and axial wavenumbers, respectively. Generally speaking, distinguishing between TE n,k,p and TM n,k,p modes for a cylinder of a finite length is justified only for n = 0. For other cases, the polarisation is hybrid [37]. In the case of arbitrary n, k, p the mode polarization is mixed. Thus, under the terms "TE" or "TM" we further imply the dominant polarization of the modes.
The low-frequency spectrum of the dielectric cylinder under consideration consists of three types of modes. The modes with the axial index p = 0 and azimuthal index n = 0, 1 demonstrate a small frequency shift with changing r/l. They are formed mainly due to reflection from a side wall of the cylinder, and they could be associated with the Mie resonances of an infinite cylinder (see Figs. 1a and 2a). The modes with the indices p > 0 and n = 0, 1 demonstrate a strong shift to higher frequencies with increasing aspect ratio r/l. They are formed mainly due to reflection from the faces of the cylinder, and they could be associated with the Fabry-Perot modes (see Figs. 1a and 2a). The modes with the azimuthal index n = 2, 3, ... are formed due to the wave incident on the side wall of the cylinder at the angles bigger than the total internal reflection angle, which is about 6.4 degrees for ε 1 = 80. Therefore, they are close in nature to the whispering gallery modes (see Fig. 2a) and their high Q factor is explained by total internal reflection but not by destructive interference as we have for quasi-BIC. Properties of WGMs are well-studied (see, e.g., Refs. [5,38,39]) and further we focus on the Mie-like (TE 1,1,0 ) and Fabry-Perot-like (TM 1,1,1 ) modes. Their electric field distributions are shown in Fig. 1b.
In quantum mechanics, in the simplest case, the system with light-matter interaction is described by a sum of Hamiltonian without interactionĤ 0 and an interaction potentialV (see, e.g., Ref. [40]). The diagonal components ofV are responsible for energy shift and the off-diagonal components are responsible for the coupling. The interaction results in a mixing of the light and matter states and in appearance of an avoided resonance crossing -the characteristic feature of the strong coupling regime [41].
In electromagnetism, due the fact that a resonator is an open system, description of the interaction between the modes becomes more complicated. There are two main approaches describing the interaction between the modes in open system. The first one considers an open system (dielectric cylindrical resonator in our case) as a closed system with non-radiating modes |ϕ a and |ϕ b interacting with a continuum of the radiation modes outside of the resonator in accord with the Friedrich-Wintgen mechanism [10] (Fig. 2b). The difficulty of this method is to correctly define the basis of the non-radiating modes and their coupling constants with the radiation continuum. In the second approach, the resonator is primordially considered as an open non-Hermitian system, characterised by a complex eigenfrequency spectrum. In this approach, a small change of resonator shape could be described as a perturbation δε(r) playing a role of the interaction po-tentialV between modes |ϕ a and |ϕ b (see Refs. [6,42]). In our case, a perturbation δε(r) is responsible for change of the aspect ratio of the cylindrical resonator (Fig. 2c). This method is well-developed for quantum mechanics and electrodynamics [43][44][45][46][47]. It allows to find spectrum, eigenmodes, and interaction constants straightforwardly from the Maxwell's equations (see Appendix B).
For the cylindrical resonator, the strong coupling between the Mie-like and Fabry-Perot-like modes is clearly manifested in the map of the SCS as avoided resonance crossing points (Figs. 1c and 1d). The most pronounced regions of the avoided resonance crossing are marked by red ellipses in Fig. 1d. More detailed analysis shows that in the vicinity of the avoided resonance crossing, the Q factor of one the coupled mode becomes very high that corresponds to the appearance of a quasi-BIC. The dramatical increase of the Q factor is a result of destructive interference between the modes with similar radiation patterns in far field.

B. Analysis of the Fano parameter
The scattering of light by high-index dielectric particles is governed by the Mie resonances of the structure. In the case of highly-symmetric geometries, such as infinite rods [48,49], spheres [50] or core-shell particles [51] it was shown by means of the Mie theory that the scattering cross section represents a series of Fano resonances where each resonance can be described by the Fano formula [52,53]. For other designs of resonators, e.g. finite dielectric cylinders, analytical Mie solution does not exist because the variables of the Maxwell's equations are not separable. However, description of the SCS by the Fano formula is still convenient but, in general, Fano parameters are introduced phenomenologically. The most challenging problem is determination of exact expressions for Fano parameters in concise and clear form. It is worth mentioning that the Fano formula for resonators with complicated geometries can be obtained but for special assumptions [54].
In this paper, we derive an elegant analytical solution for the finite-size cylinder scattering problem by proving rigorously that SCS of a lossless dielectric cylinder irradiated by a plane wave represents the conventional Fano formula. By this approach, we investigate the strong coupling between Mie-like and Fabry-Perot-like modes and reveal the relationship between the peculiarities of the mode spectra, particularly, emergence of quasi-BICs, and the singularities of Fano parameters.
We assume harmonic time dependence of the incident field in the form E inc e −iωt and determine SCS through the extinction cross section which for the lossless case can be calculated using the optical theorem [55]. The main idea of our approach is the expansion of scattered field amplitude into the sum of independent terms where each term corresponds to an eigenmode of the cylinder. This becomes possible by applying the recently developed procedure of the resonant-state expansion that allows for rigorous characterization of eigenmode spectrum of open optical resonators [56]. The cylinder eigenmodes (or resonant states) are treated as self-standing resonator excitations with a complex spectrum describing both the resonant frequencies ω 0 and damping rates γ. Our straightforward but cumbersome calculations (see Supplemental Materials) show that the frequency dependence of the SCS could be rigorously described by the Fano formula and the Fano parameters could be expressed analytically through the material and geometrical parameters of the cavity: Here Ω = (ω −ω 0 )/γ is the relative frequency detuning, q is the Fano asymmetry parameter, A is the smooth amplitude of the peak, ∆ is the resonant phase and I bg is the background contribution describing the non-resonant scattering terms. For ω in the vicinity of one of the eigenfrequencies ω 0 − iγ frequency dispersion of Fano parameters can be neglected. The key parameter of the model which determines both q and A is the coupling coefficient between the electric field of the resonant state E rs and of the incident field E inc The developed approach allows for investigation of evolution of the Fano parameters of coupled modes in parametric space (for different aspect ratio of the cylinder).
As an example, we apply it to study the coupling between TE 1,1,0 and TM 1,1,1 modes which SCS spectrum in the vicinity of the avoided resonance crossing is shown in Fig. 3(a). The calculated resonance positions, Q factor and Fano parameters A and q together with results of extraction of the same data by numerical fitting of the SCS to the Fano formula are shown in Fig 3(b-e). Analysis of Fig. 3 reveals strong correlation between evolution of Fano parameters and quality factor of the high-frequency band. Foremost, q tends to infinity exactly at r/l = 0.543 where the quasi-BIC with high Q factor emerges. Next, at r/l = 0.59 where the SCS shows a narrow deep with a symmetric quasi-Lorentzian antiresonance and q = 0, the peak amplitude A decreases dramatically to a negligible but non-zero value. This means that mode excitation from the far field is strongly suppressed so the mode becomes almost dark. Importantly, we can claim that a quasi-BIC with high Q factor and a dark mode with vanishing intensity manifest themselves at different values of the aspect ratio. This counterintuitive result shows that quasi-BIC differs substantially from a true BIC supported by unbound structures, which is always a dark mode. While true BICs are invisible in the scattering spectrum, quasi-BICs can be tracked easily by controlling peak shapes in the SCS.
To gain deeper insight into a link between formation of quasi-BIC and peculiarities of Fano parameters in SCS spectra, we consider the radiation losses as a perturbation. This approach is natural and justified since we are working in the vicinity of quasi-BIC, where radiation losses are strongly suppressed. Remarkably, nonperturbative extension of this method is developed in Ref. [57].
other on the aspect ratio r/l. (e) Dependence of the first-order corrections to the radiation amplitudes, δDmain and δD other on the aspect ratio r/l. In general, radiation continuum represents a set of independent channels which we label by α. For a single resonator, the independent channels can be attributed to the spherical multipoles. The details of the multipole analysis of quasi-BICs are provided in Sec. II C. The coupling amplitude D α between the resonant state and the radiation continuum mode E α is given by (see Supplementary materials) According to the the reciprocity theorem, the same amplitudes D α determine coupling with the incident field [58]. Following the perturbation theory, we represent each resonant state E rs as sum of a closed resonator mode E (0) and first order correction iδE responsible for the radiation. Therefore, the coupling constant can be also expanded as D α = D (0) α + iδD α . The inverse radiation lifetime γ of the resonant state can be calculated as the sum of radiation losses into all radiation channels (like the Fermi's golden rule in quantum mechanics): Furthermore, the coupling coefficient κ between the resonant state and the incident field (see Eq. 2) can be decomposed into a series of independent contributions of all channels (see Supplemental materials). For most of the eigenmodes their radiation is mainly determined by a single dominant channel, which we denote as D main . The rest small amount of power goes via other channels, which we denote as D other = α D α (see Fig. 4a). The sum is taken over all radiation channels except the main one. In these notations, the critical behaviour of the Fano asymmetry parameter q is determined by the following simple conditions: Inspection of Eq. 4 and Eqs. 5 reveals the crucial role of the main channel in both formation of high-Q modes and critical behavior of q. If radiation to the main channel is completely suppressed (D (0) main + δD main = 0) we get a quasi-BIC, characterized by minimal radiation losses. However, the Fano asymmetry parameter q tends to infinity not exactly at quasi-BIC, but very close to it, when condition 5a is fulfilled. The dark mode (q = 0) is formed under condition 5b, when the radiation into the main and the rest channels compensate each other in the first approximation.
For deeper understanding we focus on the particular example and investigate evolution of coupling coefficients D α for the high-frequency band of the avoided resonance crossing between TE 1,1,0 and TM 1,1,1 modes. In this case the main channel represents electric dipole radiation and other channels are dominated by the magnetic quadrupole radiation. The comparison of the evolution of the inverse radiation lifetime γ, phase ∆ and the amplitudes D (0) and δD for different channels with respect to r/l is shown in Figs. 4b-e, respectively. Figures 4d and 4e show that the radiation to the main channel is suppressed in virtue of coupling between TE 1,1,0 and TM 1,1,1 modes and it is completely vanished at r/l = 0.543 where the quasi-BIC emerges. Importantly, in this regime δD main nulls simultaneously with D (0) main since both of them are proportional to the rate radiation into the main channel. Under further change of the aspect ratio r/l the amplitudes D (0) main and δD main become negative and continue to decrease. At r/l = 0.546, the perturbation δD main compensate δD other and q goes to infinity (Eq. 5a). For further increase the aspect ratio, the unperturbed amplitude D (0) other of the non-dominant channels becomes exactly opposite to D (0) main at r/l = 0.575. Thus, conditions 5b are satisfied and q = 0. Therefore, the quasi-BIC and dark mode are not appeared at the same value of the aspect ratio because of contribution of other non-dominant channels of radiation losses.

C. Multipole analysis
In order to gain deeper insight into the physics of quasi-BIC in a single resonator we illustrate cancelation of its radiation losses through the dominant channel in term of multipoles. The far-field of a single resonator could be expanded into a multipole series of vector spherical harmonics. Each harmonic plays a role of an independent radiation channel. One-to-one correspondence between eigenmodes and spherical multipoles can be established only for spherical resonators. Any mode of other res- onators is always contributed by the infinite number of multipoles. However, in this infinite series it is possible to distinguish the dominant term (dominant channel) making the main contribution to the radiated power. It is possible to show using the group symmetry analysis that for the mode TM 1,1,1 the main contribution to the radiation energy is given by electric dipole moment and the rest part of energy is mainly radiated through the magnetic quadrupole moment. For the aspect ratio r/l of the quasi-BIC, radiation through the dipole channel becomes negligible and the dominant radiation channel is the magnetic quadrupole [see Fig. 6(a)] and the radiation pattern changes dramatically [see Fig. 6(b)]. It is possible to show that all other quasi-BICs demonstrate similar bahaviour in far-field [59]. Therefore, quasi-BICs in single resonators are manifested not only in the scattering spectra as a singularity of the Fano asymmetry parameter but in the far-field since the radiation pattern changes dramatically.

D. Two-band model of strong mode coupling
In this section we discuss the mechanism of strong coupling between modes in a cylindrical resonator and the reason of complete suppression of radiation to the main channel. As it was shown in the previous section, mode coupling is realised not in real, but in parametric space. It means that the cylinder aspect ratio represents the parameter determining the strength of interaction between modes. Therefore, coupling between the modes is governed by the perturbation of cylinder permittivity (see Fig. 2c). Generally, this perturbation mixes all resonant states of the cylinder. New resonant states of the perturbed resonator can be found by means of the resonant-state expansion (see Appendix B) which is a special rigorous technique that allows for careful investigation of spectrum of open systems. However, in the vicinity of an avoided resonance crossing the interaction between modes involves only two of them. Therefore, the general form of the resonant-state expansion can be reduced to a two-band model.
We start from the two of uncoupled cylinder eigenfunctions |ϕ a and |ϕ b with eigenfrequencies ω a − iγ a and ω b − iγ b , respectively. They undergo strong coupling and the perturbed eigenfunction with eigenfrequency ω represents their linear combination |ϕ = C a |ϕ a + C b |ϕ b , where coefficients C a,b are, in general, complex. The twoband model can be written as a generalized eigenvalue problem [56] ω a −iγ a 0 0 Here the perturbation is determined by the symmetric matrix with complex-valued elements which makes the problem non-Hermitian Figure 5 demonstrates perfect coincidence of the exact dispersion of two coupled modes and the two-band model results. Since both modes possess sufficiently low radiative losses, when we change cylinder aspect ratio they couple predominantly in the near-field region inside the resonator. Therefore, when the coupling is maximal, the coefficients become C a = 1, C b = ±1 for lowand high-frequency modes in Fig. 5, respectively. We recall, that |ϕ a and |ϕ b are characterized by similar far-field pattern as they have the same mode symmetry with respect to azimuthal direction and inversion symmetry of structure (see Methods). This means that farfield distribution of radiation of the high-frequency mode in the strong coupling regime is almost suppressed, i.e.
[|ϕ b − |ϕ a ] | r→∞ 0. Therefore, the radiation to the main channel for this value of r/l becomes completely forbidden which explains to formation of high-Q quasi-BIC mode and the dark mode discussed in the previous section.
Conventional method to characterize strength of mode coupling is the Rabi frequency Ω R , which is a half of minimal distance between the dispersion curves of coupled modes (see Fig. 5). For the two resonances to be spectrally separable, the minimum mode-splitting needs to be greater than the sum of the half linewidths of the modes, which is a necessary condition to observe strong coupling [61]. For the avoided resonance crossing under consideration Ω R r/c = 0.018 which is 35 times higher than sum of the half linewidths which clearly manifests the strong coupling of modes.

E. Effect of material losses
For the cylindrical resonator analyzed above we neglect the material losses and take into account only the radiative ones. Here we analyzed the effect of material losses on the quality factor of quasi-BIC and on the conditions of strong coupling. Figure 7a shows the dependence of the total quality factor Q on the aspect ratio for the high-frequency mode in the vicinity of the avoided crossing regime between the modes TE 1,1,0 and TM 1,1,1 at different material loss level. One can see that the Q factor strongly depends on the material losses and could be substantially decreased. In the presence of material losses the total Q tot factor could be estimated as Here Q rad and Q mat are responsible for the radiative and material losses, respectively. Therefore, the results obtained for the lossless cavity are justified and Q tot ≈ Q rad if the radiative losses are dominant (Q rad > Q mat ). Since material losses decrease the Q factor, they affect the strength of mode coupling as well. Figure 7b shows how the Rabi frequency Ω R and the sum of the half linewidths change depending on the level of material losses. Strong coupling regime breaks when the resonance become spectrally inseparable which is realized for Im ε = 6. Therefore, even for relative high absorption the strong coupling can be realized, which is extremely useful for the experimental measurements described in the next Section.

F. Experimental results
Finally, we perform the experimental study to demonstrate the existence of the avoided crossing regime between the TE 1,1,0 and TM 1,1,1 resonances in the microwave frequency range. In the experiment, the plastic cylindrical vessel filled with water is placed in the middle between two antennas. The aspect ratio of the cylindrical resonator is defined by the amount of water. The photo of the experimental setup is shown in Fig. 8a (see Appendix A for details). The resonator is excited by TE polarized electromagnetic wave incident perpendicular to the cylinder axis z (see Fig. 8a). The measured dependence of the SCS of the cylindrical resonator depending on its aspect ratio is shown in Fig. 8b. The results of the numerical simulations taking into account the losses in water are shown in Fig. 8c. One can see that the experimental positions of the resonances are in a good agreement with the real part of eigenfrequencies (marked by white circles) calculated using the resonant state expansion method (see Methods for details). In spite of losses in water, which broaden the resonances, the avoided crossing regime between the TE 1,1,0 and TM 1,1,1 modes and suppression of SCS clearly manifest themselves for the aspect ratio in the range of 0.5 < r/l < 0.6. Discrepancies between the measured and calculated maps of SCS could be explained by not perfect plane wave radiated by a horn antenna and parasitic scattering from the auxiliary equipment (holder of the resonator and plastic cylindrical vessel).

III. DISCUSSION
As we mentioned above, a true BIC is mathematical abstraction and it is not practically implementable. However for periodic photonic structures with quite large number of periods, the radiation of high-Q states could be almost suppressed being much less that other loss mechanisms in the system. Such states are closest to the true BICs. Here we demonstrate that the radiation losses can be substantially suppressed, being much smaller than other losses, via BIC-inspired mechanism even in a single isolated resonator. Therefore, BIC in a finite size periodic structure and quasi-BIC in a single resonator could be indistinguishable in practice if their radiation losses will be strongly suppressed. We believe that the proposed concept of quasi-BIC in a single resonator is more favourable for compact nanophotonic applications and much easier implementable than other designs.
The developed theory predicts that the shape and amplitude of the SCS spectra represents an unambiguous indicator of the quasi-BIC. Namely, both regimes of q = 0 and q → ∞ describe important features inherent to a true BIC. The latter condition q → ∞ practically coincides with the emergence of a quasi-BIC. Therefore, a quasi-BIC could be recognized in the experimental spectra by the naked eye since the asymmetry of the Fano resonance is easily distinguishable without fitting.
The difference in the aspect ratios corresponding to quasi-BIC and singularity of q is determined by the magnitude of radiation losses to non-dominant channels. The difference becomes negligible with increase of the permittivity of the resonator. Even for high-index dielectric nanostructures in the visible and near-IR range, i.e. ε ∼ 10 − 12, the relative difference between aspect ratios corresponding to the singular q and the maximal Q factor is less than 5%. The predicted strong coupling regime and peculiarities of the Fano parameters could be observed in a wide spectral range from the visible to centimetre wavelengths. These results lift the veil on the nature of quasi-BICs in single resonators and emphasize profound relationship between the shape of Fano peak and emergence of bound states in the continuum.
Remarkably, that for BICs in periodic photonic structures, the non-dominant channels are absent because the radiation continuum is discretized [32]. Thus, the leaky modes interact with only one channel which in the simplest case represents the zero-order diffraction [62]. Therefore, leaky states can transform to true BICs, when the radiation to the main channel is completely forbidden. Since for true BIC, the system of Eqs. 5 is satisfied simultaneously, the Fano parameter q becomes illdefined, which corresponds to the collapse of Fano resonance. Furthermore, true BIC is always a dark mode, which can be easily understood from energy consumption arguments -BIC does not radiate at all, thus we are not able to pump it. But for finite-size dielectric cylinders the collapse of Fano resonance is not manifested, instead of this Fano parameter monotonically evolves as cylinder aspect ratio changes and consistently passes the values q = ±∞ and q = 0 where mode exhibits features of a true BIC state.
Recently, the study of resonant dielectric nanostructures has been established as a new research direction in modern nanoscale optics and metamaterial-inspired nanophotonics, due to their optically-induced electric and magnetic Mie resonances [63,64]. However, the Q factor of Mie resonances is about tens that is far from the values achieved in WGM resonators, photonic crystal or Bragg cavities. The proposed mechanism of strong mode coupling in single high-index dielectric resonators accompanied by emergence of quasi-BIC helps to substantially extend functionality of all-dielectric nanophotonics. This opens new horizons for active and passive nanoscale metadevices including low-threshold nanolasers, biosensors, on-chip parametric amplifiers, and nanophotonics quantum circuits.

IV. CONCLUSION
We have demonstrated that a subwavelength homogeneous dielectric resonator can support strongly interacting modes. We have shown that the strong coupling regime is accompanied by the formation of a quasi-BIC when the radiative losses are almost suppressed due to the Friedrich-Wintgen destructive interference. The analysis of the scattering cross-section reveals an abrupt change of the Fano asymmetry parameter from minus to plus infinity in the vicinity of quasi-BICs. Therefore, quasi-BIC manifests itself in scattering spectra by the symmetric Lorentzian shape. Appearance of quasi-BIC is accompanied by drastic change of far-field radiation pattern explained by suppression of the radiation through the dominant multipole. This singularity could be used as an indication of quasi-BIC. In contrast to true BIC, the Fano-resonance feature for quasi-BIC does not vanish completely since it is not completely decoupled from the radiation continuum. We have confirmed our theoretical results in microwave experiment by using a cylindrical resonator filled with water. Our results open new horizons for active and passive optical nanodevices including efficient biosensors, low threshold nanolasers, perfect filters, waveguides, and nanoantennas.  Frequency ωr/c Aspect ratio r / l Quality factor FIG. 9. Complex spectrum of eigenmodes. The spectrum is shown for the modes with the azimuthal index n = 0 (red dotted lines) and n = ±1 (blue dotted lines), which are even with respect to (z → −z) symmetry. Dot sizes are proportional to Q factor. The region of the avoided crossing between modes TE1,1,0 and TM1,1,1 is marked by the green ellipse. Calculations are performed by using the resonant-state expansion method.
V is responsible for the transformation of the sphere into the cylinder. Thus, this operatorV differs from those defined in Figs. 2b and 2c.
The resonant state expansion represents a generalization of the Brillouin-Wigner perturbation theory for non-Hermitian systems [46]. Therefore, numerical accuracy is determined by the size of the basis set N . We choose the basis in such a way that for a given orbital number l, azimuthal number n and parity we select all resonant states with frequencies lying inside the circle |ωR/c| < 10, where R is the radius of the sphere that describes the cylinder. We consider l < 80 which results in N = 1035 that is enough to achieve 99.9% accuracy for calculation of real part of frequencies. Since the perturbation V αβ conserves the axial symmetry and mirror (z → −z) symmetry, we study problem for each azimuthal index n and each parity independently.
The dependence of the complex spectrum of eigenmodes with azimuthal indices n = 0, ±1, which are even with respect to up-down reflection symmetry, vs. the cylinder aspect ratio r/l is shown in Fig. 9 by dotted lines. Dot sizes are proportional to the Q factor. We can clearly observe multiple avoided resonance crossings between modes with the same azimuthal number. In vicinity of an avoided crossing point the Q factor of one of coupled modes dramatically increases, which confirms the results of SCS calculations (Figs. 1c and 3c).
In the high-frequency region, the behavior of interaction between modes becomes more complicated, e.g. Q factors of some of n = ±1 modes remain high in a broad range of parameters x and r/l, as shown in Fig. 9. We explain this phenomenon as strong coupling between three and more eigenmodes with complex spectrum. This broadband high-Q regime will be the subject of our further investigations. In this article we restrict our studies to the mechanism of strong coupling between two eigenmodes because it illustrates the basic peculiarities of spectrum of subwavelength dielectric resonators.