Towards -Symmetric Optical Dimer Fabrication without a Light-Absorbing Material

We consider an approach to engineer an optical dimer of particles operating in the spectral region near the dipolar resonance that exhibits parity–time symmetry-like features. Both particles are assumed to be made of a gain medium with the same refractive index and extinction coefficient. We suggest introducing a gain– loss contrast by altering the radiative loss of the particles through changing their shape. To demonstrate our approach, we consider a dimer of infinite filled and hollow cylinders. We demonstrate that a larger hollow diameter leads to a stronger radiative decay. Then we find the parameters of a dimer that has an exceptional point at a real frequency and exhibits two real eigenfrequencies when the gain–loss contrast is decreased

During the last decade, a lot of research has been conducted in the field of non-Hermitian photonics [1][2][3][4]. The unintuitive physics that arises from the non-Hermiticity, such as the existence of real Hamiltonian eigenvalues in the systems that obey paritytime ( ) symmetry, and the exceptional points where the Hamiltonian becomes degenerate (nondiagonalizable), leads to a plethora of optical effects, such as single mode lasing in a two-resonator setup [5], unidirectional transmission of light without reflection [6,7], violation of reciprocity [8], and splitting of eigenfrequencies according to the root law when the degeneracy is lifted [9].
These effects have been successfully demonstrated in systems that are large compared to the wavelength, such as coupled waveguides, coupled ring and disk resonators, and flat layered structures. The non-Hermiticity is introduced by loss and gain from the medium absorption and stimulated emission under pumping, respectively [10], or by applying a spatially structured optical pumping to the gain medium [11]. To fabricate non-Hermitian metasurfaces with unique properties based on the symmetry and exceptional points, one needs to use resonant particles with a size comparable to, or smaller than the wavelength. In this case, application of spatially-structured optical pumping is hindered by the diffraction limit, while fabrication of metamaterials with a unit cell consisting of gain and lossy media requires a complex technological process.
Here, we consider the idea of shaping the resonators made of a gain medium to achieve a gain-loss PT PT contrast by increasing the radiative loss, allowing a -symmetric metasurface design solely from a gain medium operating under uniform pumping. To demonstrate the approach, we chose the simplest system to analyze, that is a dimer of an infinite cylinder made of a dielectric gain medium, and an infinite hollow cylinder made of the same material, which is shown in Fig. 1. We consider the TE polarization, where the magnetic field vector is oriented along the cylinder axes. We demonstrate that the radiative loss becomes stronger when a hole is introduced, therefore it becomes possible to achieve a gain-loss contrast in a dimer made of the same material. Then we find the geometry of a dimer consisting of a cylinder and a hollow cylinder, that corresponds to an exceptional point, and analyze its eigenfrequencies and scattering spectra in the vicinity of the exceptional point.  -symmetric and exceptional point photonics is conventionally studied in the framework of non-Hermitian two-level models well known from cavity quantum electrodynamics (QED), by assigning "levels" to the localized excitations in the resonators or waveguides [12]. Systems that are large compared to the wavelength, have been studied comprehensively and are well described by those models. However, despite the simplicity of the two-level model, rigorous treatment of the underlying electromagnetic wave propagation problem often appears to be quite complicated. The reason behind this complexity is that, in contrast to cavity quantum electrodynamics, where different kind of particles are coupled (e.g., excited electron and a photon being the levels and the coupling is carried by inter-level quantum transitions), in coupled resonators, levels and the coupling are essentially the same kind of excitation, namely, the electromagnetic field. This, in turn, leads to two consequences that complicate the problem.

OPTICS AND LASER PHYSICS
The first one is that the energy exchange between the resonators is carried by electromagnetic waves involving retarded potentials. The coupling, in turn, becomes subject to these retardation effects. In photonics, the effective two-level Hamiltonian is defined as the time evolution operator of the system, so it becomes dependent on its own eigenvalue, i.e., the frequency. However, if the resonators are large compared to the wavelength, it is possible to treat them as waveguides, which have translation symmetry allowing to introduce a wave vector, which becomes the eigenvalue of the system master equation, instead of frequency.
The second complication is the existence of the radiative loss, making the losses and coupling inseparable. In waveguide-like structures, such as resonators operating on whispering gallery modes, the radiative losses are usually negligible, with the main loss channel usually being either scattering or absorption, allowing the coupled system to be treated by a twolevel model. In contrast, when the size of the resonators is comparable to the wavelength, the radiative loss remains intertwined with the coupling [13].
Despite these difficuities, it is possible to formulate a two-level model using an approach based on the multiple scattering theory for coupled photonic resonators operating on the dipolar Mie resonance, however the effective Hamiltonian turns out to be dependent on its own eigenvalue [13]: (1) where is the dimer complex eigenfrequency, are the scattering poles of each scatterer (a cylinder and a hollow cylinder) in the dimer and are the respec- tive oscillator strengths, is the Hankel function of the first kind which describes the dipole field in the 2D space, c is the lightspeed and d is the distance between the centers of the cylinders.
To describe the electromagnetic waves around the cylinders, we use the basis of multipoles that are described by the Bessel and Hankel functions. In this case, scattering on a cylinder is described by the Lorenz-Mie coefficients. For infinite cylinders of radius , the scattering Lorenz-Mie coefficient can be expressed as (2) where is the cylinder refractive index, l is the azimutal number, which is equal to zero for dipoles, and . The scattering poles are the points where the denominator of the 2D Lorenz-Mie dipole scattering coefficient becomes zero, i.e., . By solving this equation for the , we obtain the refractive index of the cylinder, which has a scattering pole at the frequency . After that, the Lorenz-Mie coefficient can be approximated as (3) in the vicinity of the scattering pole, where is the oscillator strength.
Next, we consider a hollow cylinder, which can be treated as a core-shell particle. By matching the boundary conditions on the interfaces, one may arrive at the following formula for the scattering Lorenz-Mie coefficient of a hollow cylinder with a shell refractive index , external radius and internal radius : 11 22 21 ( ) ( ) ( ) ( ) 11 11 12 21 = , To find the scattering poles, we once again equate the denominator of the scattering coefficient to zero: (5) By solving this equation for , we can obtain the shell refractive index and the hollow radius , which correspond to a hollow cylinder having a scattering pole at the frequency . Then the scattering coefficient may be approximated in the vicinity of the pole in the same manner as in Eq. (3).
Let us study the dependence of the radiative loss with the hollow radius. Figure 2a,b show the dependences of the eigenfrequency and the eigenstate damping rate of hollow cylinders with various refractive indices on the diameter of the hollow, while Fig. 2c shows the corresponding oscillator strengths . The losses can be seen to increase with increasing the hollow diameter. To distinguish between the increase of the radiative loss and the decrease of the gain medium volume, gain media ( , shown with solid and dashed lines), zero-loss zero-gain media ( , shown with dotted lines) and absorptive media ( , shown with dash-dotted lines) have been studied. In all three cases, the losses increase while the medium volume decreases, therefore this increase can be attributed to the radiative loss.
The increase in the radiative loss is accompanied by a frequency shift, as can be seen from Fig. 2a. As we are interested in introducing a gain-loss contrast without a frequency shift, a scale transformation should be applied to the hollow cylinder to alleviate the frequency shift.
Let us now analyze the eigenfrequencies of a dimer composed of a cylinder and a hollow cylinder. The Hamiltonian is dependent on its eigenvalues, making the equation for the eigenfrequency non-quadratic. Therefore, to find the dimer eigenfrequencies, one has to search for the zeros of the function , which is defined as follows (6) We search for the exceptional point, which is the point where two or more eigenfrequencies coalesce. Using a Taylor expansion, one may easily demonstrate that such a point satisfies the following conditions f f Solving this equation system for the scattering poles , one arrives at the following expression (8) where and . Using Eq. (8), one may find the scattering poles of the particles, which, placed at the distance d, will form a dimer that is at the exceptional point, i.e., only has one eigenmode at the frequency instead of a pair of eigenmodes at different frequencies.
Using Eq. (8), we find the scattering poles of a dimer that is expected to have an exceptional point at the frequency (which corresponds to the vacuum wavelength 530 nm and the cylinder radius nm). We also fix the distance between ω , the centers at , which leads to a purely imaginary difference between the scattering poles, necessary for the symmetry regime. The following scattering poles are then found from Eq. (8): , which corresponds to , and . In order to match the frequencies and , the hollow cylinder has been upscaled, so that its external radius would become , where is the radius of the filled cylinder. The hollow radius corresponding to the scattering pole at is . The oscillator strengths of the cylinders then become and . These values have been substituted into the effective Hamiltonian at Eq. (1). Then, to check if the system is indeed at the exceptional point, a symmetric frequency detuning was introduced into the scattering poles given by Eq. (8): , and . The eigenfrequencies iγ ± of the Hamiltonian were then sought for, as the dependence of the complex detuning . The resulting eigenfrequency sheets, shown in Fig. 3, display the topology resemblant of the Riemann surface of a complex square root function, which is a characteristic fea- ω ω − = δω ture of an exceptional point [1,3,4]. We also note that, when only the gain-loss contrast is varied ( ), a behavior similar to the symmetry is observed. When the gain-loss contrast is decreased compared to its value at the exceptional point, the dimer has two purely real eigenfrequencies. Increasing the gain-loss contrast, in turn, leads to the symmetry breakinglike behavior, where a damped ( ) and an amplifying ( ) modes exist at the same frequency ( ). We, however, note, that the dimer itself does not possess the symmetry in the rigorous sense: because of a different geometry of the scatterers (a filled and a hollow cylinder), the composition of a parity and a time inversion transformations would yield a different system compared to the dimer considered.
Finally, we consider the scattering spectra on the dimer made of a gain medium with a refractive index and the extinction coefficient , consisting of a cylinder with a radius and a hollow cylinder with the external radius . To obtain the scattering spectra, we use the multiple scattering theory described elsewhere [13][14][15][16][17][18][19][20][21][22][23][24], considering only the dipolar contribution. The excitation is a plane wave with an incidence direction normal to the plane passing through the axes of the cylinders. As the dimer is made of a gain material and exists at a radiant regime, we consider the outflowing energy, which we define as the extinction cross-section taken with the negative sign. The spectra are shown in Fig. 4 as a dependence of the hole radius. Decreasing the hole size leads to a decrease in the gain-loss contrast, which is accompanied by a small eigenfrequency shift. This leads to a stronger coupling compared to the exceptional point, which leads to a frequency splitting. Increasing the hollow radius corresponds to an increase in the gain-loss contrast, which leads to the symmetry breaking. In this regime, the eigenmodes exist separately in the cylinder, which encompasses a radiant mode, and in the hollow cylinder, which has a decaying mode at the same frequency.

CONCLUSIONS
We have considered infinite dielectric cylinders with hollows and demonstrated that the radiative losses increase with increasing hollow diameter. This dependence allowed us to achieve a gain-loss contrast between a hollow cylinder and a cylinder without a hole made of the same material. Using the dipole approximation, we derived the equations for the geometric parameters of a dimer consisting of a cylinder and a hollow cylinder, that has an exceptional point at a given frequency. By selecting a purely real exceptional point frequency and a distance that corresponds to a purely imaginary difference between the scattering poles of the solid and hollow cylinder, we were able to construct a dimer with a symmetry-like behavior. With decreasing gain-loss contrast, this dimer displays two purely real eigenfrequencies, while increasing the gain-loss contrast leads to the appearance of a damping and an amplifying mode at the same frequency. We have analyzed the scattering spectra of the dimers of the solid cylinder and a cylinder with varied hollow diameter and demonstrated that decreasing the gain-loss contrast by decreasing the hollow diameter leads to emergence of a splitting of the resonance in the spectra.
Our findings might be useful for designing the metasurfaces operating almost at the symmetry and exceptional point regimes. We admit that the negative imaginary part of the refractive index as large as , that is required for our dimer to exhibit the symmetry-like behavior, can not be realized in optical gain media, which are currently limited by . Such a strong gain is required to compensate the strong radiative decay of the dipole resonances. At the same time, the radiative losses, if not compensated, will overwhelm the loss contrast, making the symmetry effects impossible to observe. Therefore, considering the geometric resonances of a higher order, such as quadrupole or octopole states, which are characterized by a higher quality factor, seems to be a way to alleviate this problem. In particular, halide perovskites might be a suitable gain medium, as the compensation of the radiative losses in quadrupole and octopole states by gain, leading to single-particle lasing, has been experimentally demonstrated on this platform [25].
We considered infinite cylinders due to the simplicity of analysis. To fabricate real structures one PT would need particles with finite sizes, such as finite cylinders. Hollow core waveguides can be used for implementation in the terahertz range, as their technology of fabrication with the diameters as small as 10 μm is well established [26]. In optical range one may use hollow cylinders fabricated by electron beam lithography, which may have a radius as small as 100 nm. Heights of such cylinders can exceed the radius by several times [27]. However, in this case the resonance frequency is shifted from the value predicted by Mie theory (about 10% for the cylinder height being equal to three times the radius [28]) due to the coupling between the Mie modes and the Fabry-Perot modes in the cylinder, so one would have to take extra steps to treat this shift.

TOWARDS
-SYMMETRIC OPTICAL DIMER FABRICATION PT 809