Surface Topological Plexcitons: Strong Coupling in a Bi 2 Se 3 Topological Insulator Nanoparticle-Quantum Dot Molecule

: Strong coupling of quantum states with electromagnetic modes of topological matter offer an interesting platform for the exploration of new physics and applications. In this work, we report a novel hybrid mode, a surface topological plexciton, arising from strong coupling between the surface topological plasmon mode of a Bi 2 Se 3 topological insulator nanoparticle and the exciton of a two-level quantum emitter. We study the power absorption spectrum of the system by working within the dipole and rotating-wave approximations, using a density matrix approach for the emitter, and a classical dielectric-function approach for the topological-insulator nanoparticle. We show that a Rabi-type splitting can appear in the spectrum suggesting the presence of strong coupling. Furthermore, we study the dependence of the splitting on the separation of the two nanoparticles as well as the dipole moment of the quantum emitter. These results can be useful for exploring exotic phases of matter, furthering research in topological insulator plasmonics, as well as for applications in the far-infrared and quantum computing.


Introduction
Strong coupling in light-matter interactions has long been a topic of significant interest, both because of its unique features, as well as its potential importance to a number of current and future applications [1,2].In the strong coupling regime, the rate of coherent exchange of energy overcomes the losses in the coupled light-matter system, and a new hybrid energy structure emerges.In this modified energy spectrum, we observe the emergence of characteristic splittings between the original states, with energy separations related to the coupling strength of the interaction.The strong coupling regime is characterized by an increased coherence and is often associated with the formation of new hybrid, polaritonic-like states, which have drastically different features when compared with the non-interacting states.Among others, strong coupling is expected to be of importance in nonlinear optics and sensing [1], quantum information processing [2], thresholdless lasing [3] and even in chemistry [4].
Plasmonic states, thanks to their intrinsic strong enhancement and subwavelength confinement of the electric field facilitate strong coupling in light-matter interactions [5,6].Because of their unique features, the use of plasmonic states has been a staple strategy for a variety of applications where an enhanced light-matter interaction is important.Moreover, plasmonic states have also been employed as a means for the creation of systems with enhanced performance and unique behaviour [6,7].The potential of plasmonics has been explored for quite some time in the development of all-optical circuits, improved sensors and imaging, metamaterials and theranostics [8,9].In addition, the evolving field of quantum plasmonics, which considers the quantization of electrons and/or the optical fields involved, promises advances and applications in areas such as integrated nanophotonics, energy harvesting, sensing and quantum information technology [7,10].
One particular research direction in quantum plasmonics has been the occurence of strong coupling between excitons and plasmonic excitations [11].Strong coupling has already been reported for the coupling of surface plasmon polaritons with different types of emitters, such as J-aggregates [12][13][14][15][16][17][18][19][20], dye molecules [21][22][23] and quantum dots [24,25].The new hybrid mode that arises from this strong coupling, the plexciton, is expected to be interesting for a variety of applications, such as light harvesting, light emitting devices and optical communications [26].Plexcitons may also offer a novel platform for exploring exotic phases of matter and many-body quantum phenomena, as well as for controlling nanoscale energy flows, the latter of which could be of potential importance to light-harvesting and all-optical circuit architectures [27].Another conceivable area of application would be quantum information technology, as strong coupling is an essential component of future systems utilizing light-matter interactions.
While plasmons are most often associated with metals, topological insulators could also offer a new interesting playground for plasmonics and plexcitonics [28].Topological insulators are a novel class of topological matter, which feature unprecedented phenomena and behaviour, such as spin-momentum locking and gapless, chiral conducting surface states protected by time reversal symmetry.These surface states can be robust against non-magnetic perturbations and can also be protected from backscattering because of T-symmetry, rendering them effectively lossless at low temperatures where lattice interactions are essentially absent.Topological insulators are being studied, among others, for their potential applications in spintronics and quantum computation, as well as for the exploration of exotic phases of matter.[29][30][31][32] Previously, there has been work exploring alternatives to metals for plasmonics [33], which are associated with considerable ohmic losses-something that can be quite limiting for applications [10].Classes of materials that have been considered include transparent conducting oxides (TCOs) [34], transition-metal nitrides [35] and 2D materials such as graphene [36].In all these cases, the materials support plasmonic excitations due to the presence of some form of free carriers, which result in a negative real permittivity response.
Recently, there has been considerable work in exploring the possibility of plasmonic excitations in bulk topological insulators (TIs) [37][38][39] as well as in topological insulator nanoparticles (TINPs) [40].These new investigations open up a new world of possibilities in applications utilizing plasmonics, strong coupling and exotic phases of matter.Recently, theoretical results from Siroki et al. [41] about a nanoparticle made of Bi 2 Se 3 , a single Dirac-cone topological insulator [42], may present yet another link between plasmonics and topological insulators.Their results suggest the existence of a novel collective excitation in a Bi 2 Se 3 nanoparticle when it is subjected to light, which the authors called a surface topological plasmon mode (SToP mode).This excitation seems to arise due to the coupling of light with electrons in the surface states and phonons in the nanoparticle.Furthermore, additional theoretical work [43] suggests that the electron surface states involved are also symmetry protected, similar to the case of a bulk crystal.This mode seems to be characterized by a comparatively high absorption rate, as well as a quite narrow bandwidth.This interesting combination of features could potentially make this excitation a good candidate in achieving strong coupling with other excitations, such as excitons in quantum emitters.
In this work, we computationally investigate the possibility of achieving strong coupling between this SToP mode and the exciton of a simple two-level quantum emitter.We considered a system comprising of a spherical TINP made of Bi 2 Se 3 , separated by some distance from a two-level quantum emitter, coupled via an external sinusoidal electric field, with polarization parallel to the main axis of the system.
To model the behaviour of the interacting nanoparticles, we followed a quantum master equation-density matrix approach for the quantum emitter, while we used a classical dielectric function approach for the TINP.By working within the dipole approximation and later applying the rotating wave approximation, we derive a set of equations describing the evolution of the quantum emitter's states.The solutions from these equations are then used in estimating the absorption of the interacting TINP-QE system.
By studying the absorption spectrum of the system, we have found that splitting can occur, which suggests that our system is operating in the strong coupling regime.This, in turn, suggests that a new coherent state is created, a localized surface topological plexciton, which results from the strong coupling between the SToP mode in the TINP and the exciton of the quantum emitter.
The paper is organized as follows.We first present a concise overview of the effective permittivity model for a Bi 2 Se 3 TINP subjected to illumination, as calculated in [41], reproducing and presenting the most relevant results for our investigation.This section also features a brief discussion on the various excitations of the irradiated TINP.In Section 2 we present in detail the model we used to describe the absoprtion of a single semiconductor quantum dot(SQD)-TINP dimer.The central equations and expressions, as well as the general strategy for studying the system are also discussed there.In Section 3 we present numerical results demonstrating the emergence of Rabi-like splitting in the absorption spectra for the interacting SToP and exciton mode.We also show how the splitting (and thus the coupling) is affected by the separation of the two nanoparticles as well as the dipole moment of the quantum emitter.We then conclude with a brief discussion on the results, the method employed, as well as on possible future research directions and applications.

Topological Insulator Nanoparticle Subjected to Light
Before we move on to the interacting TINP-QE system, it is very helpful to first study the behaviour of a single Bi 2 Se 3 nanoparticle interacting with light.In the first part of this section, we present the permittivity model we used for the TINP in our numerical calculations.In the second part, there is a brief presentation and discussion on the excitations of the nanoparticle, focusing on the surface topological plasmon (SToP) mode of the TINP.

Permittivity Model
Based on the model for a spherical topological insulator developed by Imura et al. in [44], Siroki et al. have derived the energy structure for a spherical Bi 2 Se 3 nanoparticle by using time-dependent perturbation theory [41].This model should be valid for energies close to the Dirac point [44,45] and for relatively weak field intensities [41].
According to this model, the surface states of the Bi 2 Se 3 nanoparticle become discrete, at energies ±A/α, ±2 A/α, . . .relative to the Dirac point.Here, A = 3.0 eV•Å is a constant obtained from density functional theory calculations [45], while α is the radius of the nanoparticle.In Figure 1 we can see a simple diagram depicting these discrete surface states, when only states below the Dirac point are occupied.This is the case that will be presented in Section 4. Similarly to [41], the bulk dielectric function model we employ is of a simple Lorentzian form This model includes the contributions from three different excitations in the bulk material: α and β transverse phonons as well as free carriers f due to defects in the crystal structure.
As shown in the same paper using perturbation theory, the absorption cross section of a Bi 2 Se 3 nanoparticle in the quasi-static limit is given by where the term δ R appears because of the interaction between the surface states and the external field.In the special case where we consider the Fermi level below the Dirac point (E = 0), this term is given by It is important to note here that the expression of δ R is only valid for sufficiently weak electric fields, a fact which we have taken into account for our numerical calculations.It should also be noted that the applied electric field is assumed to be perpendicular to the c axis of the crystal.
It is also interesting to note that Equation ( 2) is reminiscent of the classical equation for the absorption cross section of a dielectric sphere.We can thus approximate the dielectric function of the spherical TINP using the simple expression This approximation is used extensively in the model we used for the calculation of the absorption of the system, as described in Section 3.

The SToP Mode
In [41], using Equation (2), the absorption cross section of a Bi 2 Se 3 TINP was calculated, showing a novel excitation which the authors called a SToP mode (surface topological plasmon mode).The peak corresponding to this mode seems to be characterized by a comparatively high absorption rate, as well as a quite narrow bandwidth.This interesting combination of features could potentially make this excitation a good candidate for achieving strong coupling with quantum emitters.
For the sake of completeness, we reproduce the power absorption of the nanoparticle by using Equation (29) (see Section 3), which, for a TINP in a vacuum, yields an equivalent expression to Equation (2) (with the difference of a constant factor): In Figure 2 we can see the excitations of the TINP for varying values of its radius.Comparing with the case where we ignore the effect of the surface states (a), we can clearly see that including them (b) gives rise to an additional excitation, the SToP mode.The proposed mechanism for the appearance of the SToP mode in [41] is that the surface states act as mediators of the interaction between light and the α phonon, which normally does not absorb.The second graph suggests that for relatively small TINP radii, the localized surface phonon polariton (LSPhP) couples strongly to the SToP mode, resulting in splitting behaviour.In Figure 3, the SToP mode, along with the other excitations of the TINP, are shown for a nanoparticle with a radius of 10 nm.When we ignore the effects of the surface states (red graph), we see the same three excitations that are also visible in Figure 2a, namely, from left to right: the LSPP, the β phonon and the LSPhP.However, when the surface states are taken into account (green graph), a new peak emerges, the SToP.It is exactly this case that we present in more detail in Section 4.
A short note should be made here on how we have determined the nature of these excitations.By comparing the positions of the peaks with the values of the dielectric function, we can conclude that while the β phonon peak is a simple resonance, the SToP mode, the LSPP and the LSPhP all occur approximately when ℜ[ϵ TINP ] = −2 and ℑ[ϵ TINP ] is relatively small, suggesting that these excitations are Fröhlich-type resonances.This hints at the SToP mode being a plasmonic mode, a collective electromagnetic excitation occuring near the surface of the nanoparticle.Another indication for this is the association with free carriers, namely the electrons in the surface states, analogously to the case of the LSPP with the free carriers in the bulk (see Equation ( 1)).At this point, we should also note why we have chosen to present this particular case of nanoparticle radius in this work.Firstly, it seems that around this nanoparticle radius there is a significant interaction between the surface states and the LSPhP (see Figure 3).Furthermore, comparisons with tight-binding calculations presented in [41] suggest that this nanoparticle radius could be a suitable lowest bound for the validity of this model.This bound is also interesting as the energy level spacing increases with decreasing radius, and a greater spacing could potentially mean that the topological states are more robust to thermal (scattering) effects at finite temperatures (see also discussion in [41]).Considering all of the above, this lowest bound appears to be a very interesting case to present in this article.

Interacting TINP-QE
To investigate whether there exists a strong coupling regime for a TINP-SQD dimer, we derive a formula for the absorption of the dimer.Then, by calculating the corresponding absorption spectra, the presence of splittings or avoided-crossing areas imply the occurence of strong coupling between excitations of the nanoparticles.This section includes a short presentation of the model and the expressions describing the absorption of the interacting QE-TINP system.The interaction between a two level quantum emitter and a classical nanoparticle is theoretically studied in the manner of [46].

The Model
We consider a spherical TINP of radius α, interacting with a semiconductor quantum dot, with the two nanoparticles separated by a distance R. We also assume the two nanoparticles are placed in a matrix with dielectric constant ϵ B and the system is subjected to an external polarized field of the form E = E 0 cos(ωt).
We further assume that the quantum emitter can be described as a simple two-level system with transition energy hω 0 , transition dipole moment µ and dielectric constant ϵ S .We will treat the quantum dot quantum mechanically using a density matrix approach, while we will consider the TINP as a classical spherical particle, described by a dielectric function ϵ TINP .The calculation setup is shown in Figure 4. Our system consists of an interacting topological insulator nanoparticle (TINP) made of Bi 2 Se 3 and a two-level quantum emitter (QE, here denoted SQD, for semiconductor quantum dot), and the system is irradiated by a harmonic exciting field of intensity I and frequency ω.The TINP has a radius α and is described by a permittivity ε TINP , while the QE is described by a transition dipole moment µ and a transition frequency ω 0 .

Dipole Approximation
We use the following Lindbladian master equation to model the time evolution of the states of the quantum emitter [47,48]: In this expression, Γ ij models the relaxation rates of the states, while ρ is the density matrix of the system at a state of equilibrium.It should be noted that Γ ij acts on the (i, j) element of (ρ − ρ), and we assume that Γ ij = 1/τ 0 for i = j, while Γ ij = 1/T 20 for i ̸ = j, where τ 0 , T 20 are the characteristic relaxation times of the density matrix.More specifically, in this model, τ 0 is a depopulation (decay) time affecting the diagonal elements and T 20 a dephasing time affecting the off-diagonal elements of the density matrix.
The Hamiltonian of the two level system in the dipole approximation is written as The electric field on the quantum emitter, E SQD , can be written as the sum of the externally applied field and the field induced on the emitter by the TINP where ϵ effS = 2ϵ B +ϵ S 3ϵ B is the effective permittivity of the SQD (see [49]) and s α = −2(1) when the applied field is parallel (perpendicular) to the major axis of the system.The polarization of the TINP, P TINP , will then be given by the Hermitian expression where γ = ϵ TINP −ϵ B ϵ TINP +2ϵ B (see [49]) and Ẽ(±) TINP are the positive and negative frequency parts of the electric field on the TINP.In turn, the electric field on the TINP, E TINP , will be given as the sum of the external plus the induced field from the SQD The polarization of the SQD can then be found by using the non-diagonal elements of its density matrix [48] P SQD = µ(ρ 12 + ρ 21 ) Factoring out the high-frequency time dependence of the off-diagonal terms of the density matrix we get (10) Substituting P SQD into E TINP and then subsequently into P TINP , we can use the resulting expression to find E SQD in terms of ρ12 and ρ21 where we have defined the quantities (as in [46])

Interaction Picture and Rotating Wave Approximation
Before we apply the rotating wave approximation, we will first transform the Hamiltonian to the interaction picture.We thus write the Hamiltonian in the form where H 0,S = hω 0 α † α is the unperturbed Hamiltonian of the two-level system while H 1,S = −µE SQD (α + α † ) describes the perturbation due to the external field as well as the induced field from the TINP.To change the Hamiltonian to the interaction picture we apply the unitary transformation U ≡ e iH 0,S t/h Now, we need to apply an appropriate transformation to the master equation in order to make use of the transformed Hamiltonian.The most obvious method would be by applying the same transformation It can easily be shown that the transformed master equation can be rewritten in the simpler form where If we choose the representation g = (1 0) T and e = (0 1) T , then it can easily be shown that the first term is while the last term in the master equation becomes The Hamiltonian in the interaction picture is It is easy to see that H 0,I = H 0,S , whereas the second term is non-trivial Therefore we have hω 0 (21) and substituting for the expression of E SQD and by applying the RWA, we finally find that

Time Evolution of the QE States
Now we can substitute the expressions we have calculated above into our transformed master equation eventually giving us a set of four coupled differential equations: If we define ρ12 = A + iB, ρ21 = A − iB and ∆ = ρ 11 − ρ 22 then we can find three coupled differential equations for A, B and ∆: The solutions of these coupled equations will be used to calculate the absorption rate of the system, as described in the following section.

Absorption of the System
For the absorption of the system, we assume that we can write the total absorption as the sum of two components: one contributed by the quantum emitter and one by contributed by the TINP, i.e., For the quantum emitter, the absorption can be calculated using the simple expression which of course can be written in terms of ∆ For the TINP, we can use where j is the current density, whereas E inside TINP is the electric field inside the nanoparticle.By substituting Equation ( 9) in Equation ( 8) and by using Equation ( 10) with ρ12,21 = A ± iB, we can write the field in the TINP in the form E TINP = E C cos(ωt) − E S sin(ωt), where The electric field inside the TINP can then be written in the convenient form with ϵ effT = (2ϵ B + ϵ TINP )/3ϵ B being the effective dielectric function of the TINP.
For the current density, we use j = dP TINP /dt V , where V is the volume of the spherical nanoparticle.Using E TINP = E C cos(ωt) − E S sin(ωt) in Equation (7), we can write By calculating j and substituting in the expression for Q TINP , Equation ( 29), we can eventually get the expression To summarize, we have found the steady-state solution of the density matrix, expressed through the quantities A(t), B(t) and ∆(t).Using these steady-state values, we can then calculate the absorption rates Q SQD , Q TINP of the SQD and TINP, respectively.

Numerical Results and Discussion
In the previous sections, we have presented an overview of the excitations of a Bi 2 Se 3 TINP of radius α = 10 nm (Section 2.2), as well as of the model we have employed for the interacting TINP-quantum emitter system (Section 3).In this section, after first presenting our choice of parameters and giving a brief description of our methodology in studying the system, we show how the interaction between the SToP mode of the TINP and the exciton of the QE gives rise to a topological plexciton.

Problem Parameters
We assume a spherical Bi 2 Se 3 TINP of radius α = 10 nm, modeled by the dielectric function described in Section 2.1.We include the effect of the surface states in the response of the nanoparticle, expressed through the perturbation term δ R , which, assuming the Fermi level is below the Dirac point, is given by Equation (3).
We assume that the quantum emitter has a constant dielectric function ϵ S = 6ϵ 0 and that its dipole moment is µ = 0.65 e•nm, though we allow its value to vary in some cases in order to show its effect on the splitting.For the relaxation times τ 0 , T 20 we have assumed a linear dependence on the transition frequency: τ 0 = 0.8 ns • ω 0 /2.5 eV, T 20 = 0.3 ns • ω 0 /2.5 eV, where ω 0 is the exciton transition frequency and 0.8 ns, 0.3 ns the relaxation times at ω 0 = 2.5 eV.All these correspond to typical values for CdSe-based quantum dots [46, [50][51][52][53][54].For the SToP mode at a TINP radius a = 10 nm (≈36 meV), this provides relaxation times on the order of ps (τ 0 ≈ 11.5 ps, T 20 ≈ 4.3 ps).Comparing with [55], concerning lifetimes of SQD intersublevel transitions in the terahertz range, this seems to be a decent (and maybe a quite conservative) approximation.Lastly, for the initial conditions of the density matrix we have taken A(0) = B(0) = ∆(0) = 0, which is a simple mixed state of the ground and excited states of the quantum emitter.
As for the other parameters, the separation distance R could, in principle, take any value greater than the sum of the radii of the two nanoparticles.In our study, however, we used the radius of the TINP as the lowest limit to the separation, to identify the limits of the splitting when the emitter is of a very small (negligible) size.We have assumed a field intensity I = 1 W/cm 2 , with the electric field applied along the main axis of the system (so s α = 2).For the background material we have assumed that ϵ B = ϵ 0 .

Method of Study
In our previous analysis of Section 3, we have applied the dipole and rotating-wave approximations to derive equations to describe our system.Therefore, for any choice of ω 0 , we only study the absorption spectrum for exciting field energies around this value.For this system, we generally choose ω 0 close to the SToP peak for a TINP of radius α = 10 nm, which is around ≈36 meV.

Strong Coupling between the SToP Mode and the Exciton
Next, we study the conditions under which strong coupling occurs between the SToP mode of the TINP and the exciton of the SQD.More specifically, we observe splitting behaviour, which suggests operation in the strong coupling regime and the creation of a novel hybrid excitation, a (localized) surface topological plexciton.
In Figure 5 we see the effect of decreasing the separation distance between the nanoparticles.Here, we assume µ = 0.65 e•nm and exciton transition frequency close to the SToP mode.On each graph, the components contributed by each nanoparticle are also visible.It is clear that the splitting increases with decreasing distance, reaching a maximum of ≈2.75 meV at the limiting case R = 10 nm.The strong coupling is also evident from the fact that the splitting is visible in both components of the absorption.
In Figures 6 and 7 we see the splitting when we vary the transition frequency of the quantum emitter for specific values of R and µ.In the same figures, other than the total absorption of the system (top graph), we can also see the components contributed by the TINP (bottom left) and by the quantum emitter (bottom right).Two cases are presented: in Figure 6, for R = 20 nm, and in Figure 7 for R = 10 nm, for the same value of the dipole moment µ = 0.65 e•nm.By comparing the two figures, we can see that once again the same trend is visible.Namely, when we decrease the separation distance, we observe a marked increase in the splitting, and thus in the coupling between the two excitations.As in the previous case, we also identify the same splitting behaviour appearing in both components of the absorption, implying the increased coherence of the state.The figure shows the power absorption spectrum Q (10 −10 W) for decreasing values of R (nm), for a constant emitter dipole moment µ = 0.65 e•nm.The horizontal axis is the difference between the exciting frequency ω and the QE transition frequency ω 0 , with ω 0 ≈ 36 meV in order to match the spectral position of the SToP mode.On each graph the absorption components contributed by the quantum emitter (smooth line) and the TINP (textured line) are also visible below the total absorption line, allowing us to observe that the splitting also appears in the absorption components of the system.These results suggest that for a sufficiently small separation of the two nanoparticles we can achieve strong coupling in the system.
The results when varying µ with constant R follow the inverse trend of varying R with constant µ, i.e., the splitting increases by increasing the dipole moment µ.Both of these trends can also be summarized by the graphs in Figures 8 and 9.In Figure 8, we concisely see the effect of varying the separation distance R with constant µ = 0.65 e•nm.We can see that the splitting in the total absorption emerges at values around 20 nm, while it becomes visible in both components at values below 16 nm.These graphs suggest a maximum splitting of about 2 meV at the limiting case R = 10 nm.
In Figure 9, we see the effect of varying µ with R = 10 nm.In this case, we see that the splitting emerges for quite small values of the dipole moment, apparently around µ ≈ 0.1 e•nm.These graphs suggest a splitting of almost 3 meV at µ = 0.7 e•nm for the limiting case of touching nanoparticles.If we use a value of 2 meV as indicative of the possible size of the coupling and compare it with the exciton transition frequency, we obtain the ratio: Ω/ω 0 = 2 meV/36 meV ≈ 0.056, which should place the system well within the strong coupling regime according to common classification schemes [56,57].

Outlook
At this point, it might be conducive to highlight a few limitations of the model and indicate a few possibilities for future research.A first thing to notice is the perturbative approach in modeling the dielectric function, preventing us from studying the system beyond the weak-field limit.Therefore, developing a non-perturbative approach could allow for a richer exploration of the interaction.Another thing to consider would be working outside the dipole and rotating-wave approximations, or even attempting a fully quantum treatment of the problem, in order to obtain a fuller and more accurate description.
Moreover, an important assumption that was made in this model was that the system operates at a very low temperature, allowing us to largely ignore phonon interactions in the nanoparticles and from any surrounding medium.At higher temperatures, however, we would likely need to somehow take these effects into account.However this is not crucial for the study presented here since recent theoretical work [43] seems to support the symmetry protection of the surface states in TINPs, which is expected to limit phonon interactions by preventing backscattering.Moreover, recent experiments at room temperature [58] have pointed out that the low-temperature effects discussed here can realistically persist at higher temperatures.
In terms of applications, these results may be useful for the exploration of novel, exotic phases of matter or to further research in the area of topological insulator plasmonics, but also for applications in the far-infrared and quantum computing.For example, we could think of sensing applications in the far infrared by studying the effect of the surrounding or adjacent material or the influence of the separation distance on the properties of the system.We could also consider using ensembles of QEs and/or ensembles of TINPs instead of a simple nanoparticle dimer.One possible example would be an implementation of a "spaser" making use of the SToP mode in the TINP and using a set of surrounding quantum emitters as the active medium [59].
It might also be interesting to consider placing similar systems in (narrow linewidth) optical resonators, for example to possibly suppress interactions with phonons in the nanoparticles and in the surrounding medium [60,61].As was noted in [41], this system could be of interest in quantum information applications because of the ability to adjust the Fermi level of the TINP and the different response of the surface states to circularly and linearly polarized light.Since the topological nature of the surface states [43] has been experimentally demonstrated for TINPs [58], the discretized states of the TINP could then be used as basis states for a qubit that is more resistant to decoherence.Finally, the ability to control and manipulate the strong coupling in the reported system may facilitate the creation and preservation of quantum coherence and entanglement, both crucial aspects in quantum computing.

Conclusions
In summary, we report a novel hybrid mode, a surface topological plexciton, arising from strong coupling between the exciton of a two-level quantum emitter and the surface topological plasmon (SToP) mode of a Bi 2 Se 3 topological insulator nanoparticle (TINP) [41].Our results suggest that the plexcitonic modes may exhibit a vacuum Rabi splitting of at least 2 meV for a TINP with a radius of 10 nm (SToP mode ∼36 meV) and a quantum emitter with dipole moment µ = 0.65 e•nm, which corresponds to a ratio Ω/ω 0 = 2 meV/36 meV ≈ 0.056, suggesting that our system is operating in the strong coupling regime.The presence of strong coupling in this novel system may open up new possibilities for exploring exotic phases of matter and plasmonics in topological insulators, as well as for applications in the far-infrared technology and in quantum computing.

Figure 1 .
Figure 1.The discrete energy levels of a Bi 2 Se 3 nanoparticle subjected to light [41].Here we assume that only energy levels below the Dirac point are occupied.

Figure 2 .
The excitations of the Bi 2 Se 3 nanoparticle: Power absorption heatmap (10 −10 W) as a function of the exciting energy ω (meV) (vertical axis) and the nanoparticle radius α (nm) (horizontal axis).(a) Excitations of the TINP when we ignore the effect of the surface states.From top to bottom, the visible excitations are: the localized surface phonon polariton (LSPhP), the β phonon and the LSPP.(b) Excitations of the TINP when we include the surface states.Here, we focus on radii below 40 nm, where the interaction between the LSPhP (blue line) and the surface states (red line) is clearly visible.We can now see an additional excitation, the SToP mode, located between the LSPhP and the β phonon.

Figure 3 .
Figure 3.The excitations of a 10 nm Bi 2 Se 3 nanoparticle.Red dotted line: Excitations of the TINP when we ignore the surface states.Only three excitations are visible in this case (left to right): LSPP, β phonon and a localized surface phonon polariton.Green continuous line: Excitations of the TINP when we include the effect of the surface states.In this case we see an additional excitation, the SToP mode.In particular we can see (from left to right): LSPP, β phonon, SToP mode and the localized surface phonon polariton.

Figure 4 .
Figure 4.Our system consists of an interacting topological insulator nanoparticle (TINP) made of Bi 2 Se 3 and a two-level quantum emitter (QE, here denoted SQD, for semiconductor quantum dot), and the system is irradiated by a harmonic exciting field of intensity I and frequency ω.The TINP has a radius α and is described by a permittivity ε TINP , while the QE is described by a transition dipole moment µ and a transition frequency ω 0 .

Figure 5 .
Figure 5. From weak to strong coupling: emergence of splitting as we vary the separation distance R.The figure shows the power absorption spectrum Q (10 −10 W) for decreasing values of R (nm), for a constant emitter dipole moment µ = 0.65 e•nm.The horizontal axis is the difference between the exciting frequency ω and the QE transition frequency ω 0 , with ω 0 ≈ 36 meV in order to match the spectral position of the SToP mode.On each graph the absorption components contributed by the quantum emitter (smooth line) and the TINP (textured line) are also visible below the total absorption line, allowing us to observe that the splitting also appears in the absorption components of the system.These results suggest that for a sufficiently small separation of the two nanoparticles we can achieve strong coupling in the system.

Figure 6 .
Beginning of splitting for small separation distances.The distance between the particles is set to R = 20 nm, and the QE dipole moment is set to µ = 0.65 e•nm, for which parameters a splitting begins to emerge.The heatmaps show the power absorption (10 −10 W) as a function of the exciting light frequency ω (meV) (vertical axis) and the quantum emitter transition frequency ω 0 (meV) (horizontal axis).From the top down, the figures show (a) the total absorption of the two particles, (b) the TINP component and (c) the QE component of absorption.

Figure 7 .
Emergence of splitting for small separation distances.The distance between the particles is set to R = 10 nm, and the QE dipole moment is set to µ = 0.65 e•nm, for which parameters a splitting has formed.The heatmaps show the power absorption (10 −10 W) as a function of the exciting light frequency ω (meV) (vertical axis) and the quantum emitter transition frequency ω 0 (meV) (horizontal axis).From the top down, the figures show (a) the total absorption of the two particles, (b) the TINP component and (c) the QE component of absorption.In the total absorption shown in figure (a), the horizontal red dotted line ω ≈ 36 meV is the position of the SToP mode, while the diagonal green dotted line corresponds to when the exciting frequency matches the QE resonance, ω = ω 0 .

Figure 8 .
Separation distance and splitting.The heatmaps show the power absorption (10 −10 W) in the system as a function of the exciting light frequency ω (meV) (vertical axis) and the separation between the two particles R (nm) (horizontal axis).A splitting arises as the separation distance R between the particles becomes small enough.The QE transition frequency was set to ω 0 ≈ 36 meV in order to coincide with the SToP mode position, and the dipole moment was set to µ = 0.65 e•nm.From the top down, the figures show (a) the total absorption, (b) the TINP component and (c) the QE component of absorption.

Figure 9 .
QE dipole moment and splitting.The heatmaps show the power absorption (10 −10 W) in the system as a function of the exciting light frequency ω (meV) (vertical axis) and the dipole moment of the QE µ (e • nm) (horizontal axis), for a separation distance R = 10 nm.A splitting arises and grows as the dipole moment is incresed.We note that The QE transition frequency was set to ω 0 ≈ 36 meV in order to coincide with the SToP mode position.From the top down, the figures show (a) the total absorption, (b) the TINP component and (c) the QE component of absorption.