Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: topological insulators ( tis ) are amongst the most actively investigated systems in condensed matter physics @xcite . in reality , there is evidence for their existence in two @xcite and three @xcite spatial dimensions . due to bulk - boundary correspondence , non - trivial topological states of matter have edge states at their boundaries with peculiar transport and optical properties . for instance , the two - dimensional ( 2d ) , time - reversal symmetric quantum spin hall state  that is realized in hg(cd)te quantum wells ( qws )  is known to come along with helical edge states that are protected against elastic backscattering of non - magnetic impurities @xcite . however , not only the edge state physics of these systems is interesting but also the 2d bulk physics bears exciting novelties . the low - energy excitations of hg(cd)te qws are described by a model  the bernevig - hughes - zhang ( bhz ) model @xcite  that interpolates between the limiting cases of schrdinger and dirac fermions . this interplay between schrdinger and dirac physics constitutes an opportunity for new phenomena to emerge . we have , for instance , recently discovered collective charge excitations at zero doping , i.e. intrinsic plasmons , in this system which are absent in both separate limits @xcite .    in this article , we complement our study of the screening properties and the collective charge excitations of hg(cd)te qws on the basis of random phase approximation ( rpa ) , and hence present a comprehensive analysis of its polarization function in the static and full dynamic limit , at zero and finite doping . continuously tuning the parameters of the bhz model , we reproduce the limits of pure dirac and pure schrdinger fermions and explore intermediate regimes , in order to understand how analogies and differences emerge . we support our numerical calculations of the polarization functions with analytical expressions derived by f - sum rules . in the static limit , we calculate the screening properties due to the intrinsic system and at finite doping , analyzing the induced charge density ( with friedel oscillations ) in response to a charged impurity . different to the dirac fermion system graphene , where static screening in the intrinsic limit is momentum independent and can therefore be absorbed into an effective dielectric constant @xcite , the bhz model shows a significant momentum dependence that translates into a finite extent of the induced charge density . in the dynamic limit , we are particularly interested in a better understanding of the plasmon excitations of this system away from zero doping where we previously found a new plasmon due to the interplay between schrdinger and dirac fermion physics  @xcite . at finite doping , under certain conditions specified below that are e.g. applicable to hg(cd)te qws , we find a coexistence between this novel ( interband ) plasmon and an ordinary ( intraband ) plasmon . both plasmons can be rather weakly damped by single - particle excitations and should therefore be observable . interestingly , the two plasmons respond to the topology of the bandstructure with a distinctive behavior . they seem to merge one into the other in a normal insulating phase , while they remain clearly resolved when the system realizes a topological insulator . generally , rpa is known to provide reliable predictions at large densities and in systems with a large number of fermionic degrees of freedom . while its validity was indeed questioned for the intrinsic dirac limit , where the system is unable to screen the coulomb interaction and strong renormalization effects are expected  @xcite , rpa has been shown to yield a quantitative description of many - body effects in graphene  @xcite . it has been widely used for the study of plasmons in the dirac model , including various forms of ( multilayer ) graphene and ti surface states , see ref .   for a comprehensive review . closely related to our work , the intraband plasmons of black phosphorous have been studied on the basis of rpa and an extended version of the bhz model including anisotropy  @xcite . a similar study has been done for mos@xmath0  @xcite . our article is organized as follows . in sec . [ sec_model ] , we introduce the bhz model and present the general formalism we employ to calculate the static and dynamical dielectric function and the induced charge density . the nature of the nontrivial pseudospin , the origin of possible interband plasmons , experimentally relevant parameters and the different contributions to the f - sum rule are also discussed here . subsequently , in sec . [ sec_undoped ] , we present the static screening properties , the dynamical excitation spectrum ( new interband plasmon ) and the f - sum rule in the undoped regime . here we revisit and go beyond the results from ref .  . [ sec_doped ] , this analysis is extended to the case of finite doping where inter- and intraband excitations equally matter . we begin by discussing the ability of the bhz model to interpolate between dirac and schrdinger physics . afterwards , we have a closer look at parameters which are experimentally relevant for hg(cd)te qws , see sec . [ sec : hg(cd)te quantum wells ] . in this limit , we find a coexistence of inter- and intraband plasmons occuring for energies and momenta which are suitable for raman or electron loss spectroscopy . we close this chapter by investigating the influence of a non - trivial topology on the plasmonic excitation spectrum . finally , in sec .  [ sec_con ] , a conclusion and a brief outlook are given . the bhz hamiltonian  @xcite for a two - dimensional electron gas ( 2deg ) near the @xmath1-point has the form @xmath2 here @xmath3 are the pauli matrices associated with the band - pseudospin degree of freedom ( band @xmath4 and @xmath5 in hg(cd)te quantum wells ( qws ) ) , @xmath6 , @xmath7 with @xmath8 . the system possesses time - reversal symmetry and @xmath9 is block diagonal in the kramer s partner or spin degree of freedom . restricting ourselves to the block @xmath10 , the results can be extended to the other one by applying the time reversal operator . @xmath10 describes fermions with intermediate properties between a dirac and a conventional 2deg system . the off - diagonal term ( @xmath11 parameter ) is typical for a dirac system ( @xmath12 in graphene ) , with @xmath13 the dirac mass ( corresponding to a gap of @xmath14 ) . we consider positive and negative masses , where the latter one corresponds to an inversion of the bandstructure and the system is topologically non - trivial  @xcite . for simplicity , we restrict ourselves to a bandstructure with a minimum at the @xmath1-point , which limits the mass to @xmath15 . in analogy to a 2deg , the diagonal elements bear kinetic energy elements which preserve ( @xmath16 parameter ) and break ( @xmath17 parameter ) particle - hole ( p - h ) symmetry ( @xmath18 for schrdinger fermions with @xmath19 the quasi - particle mass ) .    the eigenstates of eq . ( [ eq : hamiltonian ] ) are described by the following dispersion and pseudospin @xmath20 with @xmath21 for valence and conduction band . note that we consider electrons to be perfectly localized on the 2d x - y plane and therefore we neglect the real shapes of the envelope functions due to the quantum confinement along @xmath22 direction  @xcite . the bhz model is characterized by intrinsic scales for momentum , @xmath23 , and energy , @xmath24 , which reflect the interpolating character of the model between dirac ( @xmath11 parameter ) and schrdinger ( @xmath16 parameter ) system . fermi momentum @xmath25 and chemical potential @xmath26 provide externally tunable momentum and energy scales , which we call fermi scales in the following . we expect the ratio between fermi and intrinsic scales to govern the physics of this system . we therefore define the dimensionless quantities @xmath27 where we set @xmath28 in the following . @xmath29 is defined to be the energy to the wave vector @xmath30 , such that @xmath31 if @xmath32 . for @xmath33 , we therefore expect intermediate physics , while in the limit @xmath34 ( @xmath35 ) the dirac ( 2deg ) physics should be recovered . the linear response of an homogeneous system to an external applied potential is described by the density - density generalized susceptibility or retarded polarization function @xmath36 . this response comprises two main phenomena : screening , described by the real part @xmath37 $ ] , and dissipation by single - particle excitations ( spes ) , given by the imaginary part @xmath38 $ ] . the polarization function in rpa yields the expression @xmath39 with @xmath40 , @xmath41 a positive infinitesimal , @xmath42 for spin degeneracy , @xmath43 and @xmath44 the fermi - dirac function with @xmath45 and @xmath46 the boltzmann constant . in the following we will assume zero temperature , @xmath47 . the overlap factor is given by @xmath48.\label{eq : f factor}\ ] ] eq . ( @xmath49 ) implies that @xmath50 is only a function of the reduced dimensionless variables @xmath51 and @xmath52 and parametrically depends on @xmath53 , @xmath54 and @xmath30 .      in the massless dirac limit ( @xmath55 ) , eigenspinors are characterized by their helicity and consequently the overlap factor @xmath56 only depends on the angle @xmath57 between @xmath58 and @xmath59 . it is strictly one ( zero ) for states with the same ( opposite ) helicity .    in the bhz model , the quadratic terms have the effect of turning the pseudospin of the eigenstates out of plane in opposite directions for conduction and valence bands at large @xmath51 , see fig . [ fig : overlap ] .    , and a ti phase ( b ) , @xmath60 . the bands are separated by an additional @xmath61 for better illustration of the pseudospin . [ fig : overlap],title="fig:",width=158 ] , and a ti phase ( b ) , @xmath60 . the bands are separated by an additional @xmath61 for better illustration of the pseudospin . [ fig : overlap],title="fig:",width=158 ]    this results in a decay of the overlap factor down to @xmath62 in the limit of a conventional 2deg system ( @xmath63 or @xmath64 ) . a finite mass @xmath65 has a similar effect , but in the limit of @xmath66 . the pseudospin turns in the same ( opposite ) direction as for the quadratic term for positive ( negative ) mass , see fig . [ fig : overlap ] . this has the direct consequence that for a normal insulator ( ni ) phase the interband overlap factor is reduced , while it is increased for a ti phase . on the contrary , a positive ( negative ) mass enhances ( diminishes ) the intraband overlap factor . this picture is also confirmed in section  [ f - sum rule introduction ] by calculating the f - sum rule . the bare coulomb interaction @xmath67 in an electron gas is modified by screening into the effective interaction @xmath68 . there , screening is described by the dynamical dielectric function . employing dimensionless units , it acquires the form @xmath69 where we have introduced the interaction strength parameter @xmath70 ( effective dirac fine structure constant  @xcite ) and the dimensionless function @xmath71 @xmath72 in graphene one finds  @xcite @xmath73 , while in hg(cd)te qws it is of the order @xmath74  @xcite . here , @xmath75 is the background dielectric constant , accounting for screening of internal electronic shells , while @xmath76 gives the dynamic screening due to electrons in the low energy bands . zeros of @xmath77 describe a density - density ( longitudinal ) perturbation of the system that it is able to sustain itself , which forms a collective mode called plasmon . it is defined by @xmath78 with the plasma frequency @xmath79 , and the finite imaginary part @xmath80 accounts for the possible damping due to single - particle excitations  @xcite . the dissipation of the interacting system , including both single - particle excitation and the plasmon mode , is then described by the imaginary part of the interacting polarization function @xmath81 . in order to compare to the non - interacting one , we will plot the normalized functions @xmath82 , \ \pi^{im}\equiv \im\left[\pi^{r}\right ] , \ \pi^{re}\equiv \re\left[\pi^{r}\right ] \nonumber\end{aligned}\ ] ] in the following , with @xmath83 .      in rpa , eq . ( [ eq : epsi ] ) characterizes the screening of the interaction between two electrons exchanging momentum @xmath51 and energy @xmath52 , by the creation of electron - hole pairs in the electron gas with the same momentum @xmath51 . if these pairs are resonant in energy @xmath84 , they correspond to a physical process leading to dissipation and a lowering of the coulomb interaction - described by the imaginary part of the polarization function , eq . ( [ eq : pi_par ] ) . when @xmath85 , we have only virtual electron - hole pairs , which either still screen the interaction , if @xmath86 , or even enhance it ( antiscreening effect ) , if @xmath87 . these effects depend on the energy of the created pair , for @xmath88 one finds antiscreening , while @xmath89 leads to a screening of the bare coulomb interaction . this can be directly seen from the definition of the polarization function , eq . ( [ eq : pi_par ] ) . for every allowed excitation , the real part of the integrand in eq . ( [ eq : pi_par ] ) becomes @xmath90}{\omega^2-\omega_{eh } \left[\boldsymbol{\tilde{x}},\boldsymbol{\tilde{x}'}\right]^{2 } } , \label{eq : integrand}\ ] ] with @xmath91 ( @xmath92 ) for intraband ( interband ) excitations . therefore every process with energy less than @xmath52 increases @xmath93 , lowering @xmath94 and thus increasing the interaction . in the intrinsic dirac system within rpa one finds @xmath95 for all energies @xmath52 where electron - hole excitations are allowed  @xcite . thus the screening effect of virtual excitations with @xmath89 cancels exactly with the one from excitations with @xmath88 , such that the only screening comes from the resonant process @xmath84 . in the bhz model , the high energy excitations become less likely as the electron and the hole band get decoupled for large @xmath52 . additionally their excitation energy is higher as in the dirac case for the same momentum @xmath51 , leading to an additional reduction of their influence on @xmath93 due to the lorentzian in eq . ( [ eq : integrand ] ) . further , low energy excitations become more important , as processes are allowed that where forbidden in the dirac system by helicity ( see sec . [ sec : bhz model no mass no d ] for details ) . combining these effects , one finds the virtual excitations which increase the coulomb interaction , @xmath88 , dominating for larger frequency @xmath52 , leading to an increased effective interaction and the possibility of intrinsic plasmons in the bhz model  @xcite . more mathematically speaking , the described effects alter the high energy behaviour of @xmath96 from a decay like @xmath97 in the dirac case to a @xmath98 decay in the bhz model , as is shown in sec . [ sec : expansion of pi ] . taking the kramers - kronig relation @xmath99 one finds directly that the real part of the polarization changes sign for @xmath100 , but not for @xmath101 . in more general terms , one can expect intrinsic interband plasmons to appear in all models for which @xmath96 decays faster as @xmath97 for high energies . the static limit of the polarization function is obtained by sending @xmath102 at finite momentum @xmath51 in eq . ( [ eq : pi_par ] ) . in this limit we can easily analyze the response of the system to the application of a static ( or sufficiently slowly varying ) external potential . an important physical problem of this kind is the screening of a charged impurity by the electronic system . the static polarization is a strictly real function , that we define as @xmath103 in a multiband system , like the bhz model , it is useful to separate the contributions to the static polarization coming from the intrinsic neutral system , @xmath104 ( obtained for @xmath105 ) , and the contribution due to a finite charge density , @xmath106 ( finite @xmath26 ) . consistently with the notation of eq . ( [ eq : pi_static ] ) , the dielectric function , eq . ( [ eq : epsi ] ) , can therefore be rearranged into @xmath107.\end{aligned}\ ] ]    from the static dielectric constant we can find the induced charge density in response to a test charge @xmath108 placed at the origin . the variation of the electronic charge density in momentum space corresponds to @xmath109 , where @xmath110 is given by  @xcite @xmath111}-1= \\   & = & n_r(x ) + n_0(x ) + n_\mu(x ) . \nonumber\end{aligned}\ ] ] here the induced charge density can be seen as a sum of three contributions of different physical nature . the first is due to the background polarization @xmath112 ( high energy polarization of the system ) , the second to the intrinsic polarization @xmath113 ( polarization of the natural system ) and the third to the polarization of the finite charge density in the system @xmath114 , with @xmath115    in real space , the density fluctuation ( using physical dimensional units ) is given by @xmath116 with @xmath117 the zero - th order bessel function .      including coulomb interaction , we now have a 4-dimensional parameter space consisting of @xmath53 , @xmath54 , @xmath30 and @xmath70 . this parameter space will be explored systematically in the following . while the exploration of the different physical behaviors featured by the bhz model in different regions of this parameter space has a clear theoretical significance , we want to stress that our discussion is also relevant for experiments . in particular , realistic parameters for hg(cd)te qw structures  @xcite are roughly @xmath118 , @xmath119 , @xmath120 and masses @xmath13 with absolute values up to several mev . the interaction strength is around @xmath121 with an average @xmath122 from the cdte substrate ( @xmath123 ) and hgte ( @xmath124 ) . considering the experimental acceptable damping rate for plasmons , we refer to experiments on the surface states of a 3d ti  @xcite . there , plasmons with a ratio of @xmath125 are perfectly resolvable . the f - sum rule for the polarization function provides the total spectral weight of all excitations in the system . it is identical for the interacting and noninteracting system , as the interaction conserves the number of particles . thus the sum rule is a powerful tool to check our numerics . additionally , it offers a deeper insight concerning the shift of spectral weight between the inter and intra spes as well as the different plasmons in the system . the f - sum rule is defined by  @xcite @xmath126=g_s\left\langle 0\left| \left[\left[n_{\mathbf{q}},h^0\right],n_{\mathbf{q}}^{\dagger}\right]\right|0\right\rangle\ \label{eq : f - sum def}\ ] ] with the density operator @xmath127 and the hamiltonian @xmath128 with @xmath10 as defined in eq . ( [ eq : hamiltonian ] ) . @xmath129 is a spinor associated with the band - pseudospin degree of freedom ( band @xmath4 and @xmath5 in hg(cd)te qws ) . the spin degree of freedom enters via the degeneracy factor @xmath130 . for the calculation we follow the steps outlined in the appendix of ref .  , where the f - sum rule for the dirac model is obtained . for the bhz model the computational steps are the same , therefore we only present important intermediate results and differences to the dirac limit . the commutator in eq . ( [ eq : f - sum def ] ) is given by @xmath131,n_{\mathbf{q}}^{\dagger}\right]&= & \sum_\mathbf{k}\left(\psi^\dag_{\mathbf{k}}h^0_{\mathbf{k},\mathbf{q}}\psi_{\mathbf{k}}-\psi^\dag_{\mathbf{k}+\mathbf{q}}h^0_{\mathbf{k}+\mathbf{q},\mathbf{q}}\psi_{\mathbf{k}+\mathbf{q}}\right)\nonumber\\ & -&2q^2\sum_\mathbf{k}\psi^\dag_{\mathbf{k}+\mathbf{q}}\left(d\sigma_0+b\sigma_z\right)\psi_{\mathbf{k}+\mathbf{q}}\label{eq : f - sum commutator}\end{aligned}\ ] ] with @xmath132 . a simple shift of the momentum sums in eq .  ( [ eq : f - sum commutator ] ) would put the first line to zero , but this is not allowed . in the same way as in the dirac system , the operators are unbounded and one has to work with a large momentum cutoff @xmath133 . while in the dirac limit one finds simply @xmath134 and the second line of eq . ( [ eq : f - sum commutator ] ) would be zero , now the latter gives rise to a contribution depending on the chemical potential , as one would expect for a 2deg . the sums in eq . ( [ eq : f - sum commutator ] ) are then converted into integrals and solved in the limit of large @xmath133 . care has to been taken when converting the momentum cutoff @xmath133 into the frequency cutoff @xmath135 , such that both integrals cover the same phase space .      for a pure dirac system one find the f - sum rule  @xcite @xmath136=-\frac{g_s q^{2}\lambda}{16 } \label{eq : dirac f - sum rule}\ ] ] where the cutoff @xmath135 is needed as the dirac spectrum is unbounded . in a 2deg system one finds @xmath137=\frac{g_s}{4 } \left(b\pm d\right ) k_{f}^{2}q^{2}=-\frac{\pi nq^{2}}{2 m } \label{eq:2deg f - sum rule}\ ] ] with @xmath138 the electron density and @xmath139\neq0 $ ] only over a finite range of @xmath140 . similar to a dirac system , the bhz spectrum is unbounded which complicates the evaluation of the sum rule and makes it necessary to introduce a high - energy cutoff @xmath141 . we find approximately for @xmath142 @xmath143 \label{eq : bhz f - sum rule}\\ =   & \frac{g_{s}}{8}x^{2}\biggl[\ln\left(\frac{2\lambda e^{-1 - 2\xi_{m}+2\left|\omega_{f}\right|}}{1 + 2x_{f}^{2}\left(1+\gamma\xi_{d}\right)+2\xi_{m}+2\left|\omega_{f}\right|}\right ) \nonumber \\    + & \frac{1-x^{2}+4\xi_{m}}{\lambda}-\frac{2x^{4}+\left(1 + 4\xi_{m}\right)^{2}-4x^{2}\left(2 + 7\xi_{m}\right)}{4\lambda^{2}}\biggl ] \nonumber \\    + & \mathcal{o}\left(\frac{\xi_{d}}{\lambda^{2}}\right)+\mathcal{o}\left(\frac{1}{\lambda^{3}}\right ) \nonumber \end{aligned}\ ] ] with @xmath144 $ ] and euler s number @xmath145 , so the leading order term diverges logarithmically with @xmath146 . this is due to the fact that @xmath147 $ ] decays like @xmath98 for @xmath148 , and not as @xmath97 as for a dirac system . the sum rule is exact up to order @xmath149 ( @xmath150 ) for finite ( zero ) @xmath151 . the f - sum rules for bhz , dirac and 2deg models are always proportional to @xmath152 in the leading order , but otherwise distinct from one another . taking the limit @xmath63 in the bhz result , eq . ( [ eq : bhz f - sum rule ] ) , gives the 2deg case , eq . ( [ eq:2deg f - sum rule ] )  ) , the same is not possible for the limit @xmath153 , as there the defined cutoff @xmath154 would go to zero . ] . we begin our discussion of eq . ( [ eq : bhz f - sum rule ] ) by comparing the contributions from the different orders @xmath155 , @xmath156 and @xmath157 . in the limit of @xmath158 we find @xmath159 and @xmath160 , thus the ratio @xmath161 determines the importance of higher order corrections for @xmath162 . we take @xmath163 for the cutoff in the following . already for @xmath164 and a maximal momentum @xmath165 , the corrections of order @xmath156 are 2% of order @xmath155 , while contributions of order @xmath157 are smaller than 0.1% . a modest cutoff @xmath166 works best for comparing eq . ( [ eq : bhz f - sum rule ] ) to numerical data , as the latter one is only given over a finite range of @xmath52 . a larger @xmath146 makes it necessary to extrapolate the data , providing a source for errors . next , we investigate changes to the f - sum rule and therefore to the total spectral weight by varying the mass . the influence of a finite mass is studied in fig . [ fig : f - sum rule ]     of the lowest order f - sum rule including mass over the one without mass . @xmath169 and @xmath32 . ( b ) ratio @xmath170 of the lowest order f - sum rule including finite doping over the one without doping , @xmath171 ( @xmath172 ) as a black , solid ( red , dashed ) line . @xmath173 and @xmath174 . @xmath164 and @xmath175 in both plots . [ fig : f - sum rule],title="fig:",width=158 ]   of the lowest order f - sum rule including mass over the one without mass . @xmath169 and @xmath32 . ( b ) ratio @xmath170 of the lowest order f - sum rule including finite doping over the one without doping , @xmath171 ( @xmath172 ) as a black , solid ( red , dashed ) line . @xmath173 and @xmath174 . @xmath164 and @xmath175 in both plots . [ fig : f - sum rule],title="fig:",width=158 ]    \(a ) for @xmath169 and @xmath32 . a positive mass lowers the f - sum rule , while a negative mass increases it linearly . this is a direct consequence from the change of the overlap factor : a negative mass enhances the coupling between the two bands , while a positive mass diminishes it , as in the latter case the pseudospins do not match . it is also consistent with the increase in the optical conductivity observed in the undoped limit with negative mass  @xcite . last , we consider the effects of finite doping . it blocks interband transitions close to the dirac point , but due to the small density of states , these transitions carry only a small spectral weight . on the other hand , doping enables intraband transitions , which carry a large spectral weight due the combined effects of larger overlap factor , density of states and smaller excitation energies compared to interband transitions . therefore , a finite doping usually increases the f - sum rule , as seen in fig . [ fig : f - sum rule ] ( b ) , where we plot @xmath170 for positive ( black , solid line ) and negative ( red , dashed line ) doping with @xmath174 and @xmath176 . a finite @xmath54 adds a term @xmath177 to the leading order of the f - sum rule , @xmath178 ( @xmath179 ) for positive ( negative ) doping . it can be seen as an increased ( decreased ) contribution from the 2deg part of the spectrum , eq . ( [ eq:2deg f - sum rule ] ) , and leads to the slight decrease of the f - sum rule for negative doping in panel ( b ) .      in order to compare the importance of different excitations in the system , one should compare their spectral weight and thus their contribution to the f - sum rule . the latter has the benefit of being independent of the coulomb interaction strength and the position of the excitation peaks , in contrast to the polarization function @xmath180 . as an example , we assume that the excitation spectrum , @xmath180 , is governed by a single plasmonic peak following a lorentzian shape with width @xmath1 and peak height @xmath181 . then the f - sum rule is proportional to @xmath182 . the value of this integral should be independent of @xmath70 and thus of @xmath183 . therefore we find @xmath184 , such that the peak height of a resonance in @xmath180 naturally has to scale with @xmath185 to fulfill the f - sum rule . we conclude that the importance of a resonance in @xmath180 should be judged by its spectral weight , which can be estimated by multiplying the peak height with its position @xmath183 . the relevant width of the peak is given by @xmath186 , with @xmath1 being the width of the resonance in @xmath180 . in this section , we focus on an intrinsic ( undoped limit @xmath105 ) bhz model system . first , we analyze the static polarization function and the static screening properties . then we consider the long wavelength limit of the dynamical polarization function , providing an analytical expansion . finally , we add some complementary arguments elucidating the origin of the new interband plasmon ( absent both in the dirac and 2deg cases ) , whose appearence for the intrinsic bhz model has been proposed in ref .  . in order to set a reference with a closely related and analytically solvable model , we discuss the static intrinsic polarization for a massive dirac limit , given by  @xcite @xmath187        \overset{m\rightarrow0}{\longrightarrow } \frac{-g q}{16 a } ,    \label{eq : pi_dirac_static}\end{aligned}\ ] ] where the index @xmath62 stands for intrinsic limit @xmath105 , @xmath188 account for possible spin and band degeneracy , and @xmath189 . when the dirac system is massless ( @xmath190 ) , @xmath191 is a linear function of the momentum @xmath192 . a finite dirac mass suppresses the polarization for @xmath193 , where @xmath191 shows a super - linear behavior . for @xmath194 , the mass is negligible instead and the result of the massless limit is reproduced . the static polarization function of the bhz model is simply obtained by direct numerical evaluation of eq . ( [ eq : pi_par ] ) at zero frequency . in fig . [ fig : pi_0_static ] , we show @xmath104 calculated for a particle - hole symmetric bhz system ( @xmath195 ) . note that we obtain the massless dirac case in the limit @xmath196 ( and therefore @xmath197 ) , where @xmath198 . a finite @xmath16 parameter determines a fundamental qualitative change with respect to a dirac system . indeed , @xmath104 reaches a maximum at @xmath199 and then decays as @xmath200 for @xmath201 as shown in the inset of fig . [ fig : pi_0_static ] .     and @xmath195.,width=302 ] a finite and positive dirac mass @xmath13 leads to a general suppression of the polarization function with respect to the massless case . in the region @xmath202 ( where quadratic terms are less important ) , @xmath104 resembles the massive dirac case , with a super - linear increase in the region @xmath203 , due to the suppression of the interband overlap factor determining a reduction of the polarization at small momentum . for intermediate values @xmath204 , analogously to the massive dirac limit , @xmath104 is approximatively linear in @xmath51 . considering larger momenta @xmath205 , the behavior is dominated by the quadratic terms and the polarization eventually vanish for @xmath206 . in general , the interplay of quadratic terms and a finite dirac mass shifts the maximum of @xmath104 . when the dirac mass @xmath13 is negative ( topological insulator phase ) , we observe a less pronounced suppression of the polarization for @xmath207 , with respect to a massive dirac system ( normal insulator ) with equal modulus of @xmath13 . moreover , on the contrary to the @xmath208 case , @xmath104 is enhanced at large @xmath51 with respect to the massless , particle - hole symmetric limit . this behavior is due to the enhanced overlap factor between electron and hole bands in the topological insulator phase .    in fig . [ fig : d](a - c ) ,     value.,width=302 ]    we analyze the effects of a finite value of the parameter @xmath151 in the bhz model , for @xmath209 , @xmath62 and @xmath210 . a finite @xmath151 breaks particle - hole symmetry by changing the effective masses of conduction and valence bands . we only found quantitative changes to @xmath104 , which is progressively reduced for increasing @xmath151 .      in a massless dirac system , where the static polarization is linear in @xmath192 [ eq .  ( [ eq : pi_dirac_static ] ) ] , the dielectric function is a constant @xmath211 therefore the intrinsic polarization contribution can be absorbed into an effective background dielectric constant @xmath212 . as a consequence , a test charge @xmath213 , placed at the origin , induces a screening electronic density @xmath214= z e \left ( \frac{1-\varepsilon}{\varepsilon}\right),\ ] ] which in real space corresponds to a screening image charge [ a fraction @xmath215 of the external one ] placed exactly at the same position @xmath216 note that the screening charge only due to the electronic system ( without background contribution ) is a fraction @xmath217 of the external one .    in a massive dirac system , the large @xmath192 behavior of @xmath191 reproduces the massless limit and therefore a screening charge given by eq . ( [ eq : charge_dirac ] ) is also developed at vanishing distances @xmath218 in response to an external test charge . however , in the long wave length limit ( @xmath219 ) @xmath191 has a superlinear behavior and thus @xmath220 . thus an induced charge density of the same sign as the external charge is developed at finite distances  @xcite [ summing up to @xmath221 , so that the test charge feels only the background screening over long distances , as expected in an insulator . for the bhz model , we find similar to eq . ( [ constant screening dirac ] ) @xmath222 in the long wavelength limit , but @xmath223 . in order to understand this , we discuss next the induced charge density in real space for the bhz model . it is given by @xmath224 with @xmath225 a natural charge density constant of the model . we note that @xmath226 is proportional to @xmath227 and @xmath70 , but @xmath226 has an additional dependence on @xmath70 ( and thus on @xmath11 ) through its integrand . it also parametrically depends on @xmath167 and @xmath151 through @xmath104 and @xmath228 .    in fig . [ fig : n_int ] ,     due to a test charge in the intrinsic limit of the bhz model for @xmath229 . the plot is invariant under a change of @xmath16 parameter and only depend on the effective fine structure constant @xmath70 . in the inset , @xmath230 calculated for @xmath231 and finite dirac mass.,width=302 ]    we plot the induced charge density @xmath226 in real space for @xmath229 with different values of @xmath70 . opposite to a dirac system , the induced charge density has a finite extent over a distance of the order of @xmath232 , which is clearly related to the decay of @xmath233 at large wavevector due to the presence of quadratic @xmath16 terms . @xmath226 decays at large distances as @xmath234 . an electron far away from this induced charge , @xmath235 , does not see the finite extent of it and is therefore screened in the same way as in the dirac system , leading to the similarity of eqs . ( [ constant screening dirac ] ) and ( [ constant screening bhz ] ) . in the opposite limit where the electron sits on top of the induced charge , @xmath236 , it does not feel it at all , resulting in no screening besides @xmath237 .    in the inset of fig . [ fig : n_int ] , we study the effect of a finite dirac mass term . with a finite @xmath167 , the induced density ( as in the case of pure dirac systems ) shows a qualitatively different behavior . @xmath226 changes sign for sufficiently large @xmath218 , ensuring a vanishing total induced charge . from a quantitative point of view , a finite negative ( positive ) @xmath13 enhances ( suppresses ) the features of @xmath226 , due to its effect on the interband overlap factor . an analytic discussion of the polarization function is only possible in the limit @xmath197 . here , we focus on the limit of vanishing mass @xmath238 to extract an analytic formula of the plasmon dispersion . an expansion of @xmath239 in @xmath51 gives , for @xmath240 @xmath241 where one finds an @xmath98 behaviour with an additional logarithmic correction for the real part in the high frequency limit . calculating the plasmon dispersion by performing an expansion of eq . ( [ eq : plasmon equation ] ) up to second order in @xmath242 , one finds the linear dispersion @xmath243 which is only valid for sufficiently large @xmath70 , such that the conditions @xmath244>0 $ ] and @xmath240 are fulfilled . the linearity of the dispersion follows from eq . ( [ eq : pi undoped expanded ] ) only by inclusion of the damping via @xmath1 . without the substitution @xmath245 , @xmath246=0 $ ] has no sensible solution for @xmath183 . the damping ratio is given by @xmath247 underlining the importance of damping in this limit . the plasmon is only well defined for a finite @xmath248 , with @xmath249 where @xmath250 sets the limit for the detectability of the plasmons , for example in the recent experiment  @xcite @xmath251 was shown to be of the order @xmath252 . ( [ eq : plasmonfrequency linear ] ) translates this into a finite momentum scale @xmath253 with the intrinsic plasmon length scale @xmath254 , given by the coulomb interaction strength times the charge decay length @xmath255 , see sec . [ sec : screening_int ] . we interpret @xmath256 as the length scale up to which charge separation due to coulomb interaction can occur and give rise to the interband plasmons , in an undoped and therefore overall neutral system .    in the opposite limit of high frequencies , the term @xmath257 spoils a simple @xmath258 behaviour of the plasmon dispersion . in this limit , we can extract the analytic form of the damping rate @xmath259 with euler s number @xmath145 , yet the plasmon dispersion can only be calculated numerically .    in the following discussion of the different excitation spectra , we will use these analytic results to check our numerics in the limits of small momenta and low and high frequencies . the non - interacting single - particle excitation spectrum is given by @xmath96 , which we plot in fig . [ fig : intrinsic spectra ] ( a )     ( a ) and @xmath180 for @xmath260 ( b ) and @xmath261 ( c ) . ( d ) shows linecuts for fixed @xmath262 with @xmath263 in black solid , red dot - dashed , blue long dashed and green short dashed lines , respectively . @xmath173 and @xmath174 . [ fig : intrinsic spectra],title="fig:",width=158 ]   ( a ) and @xmath180 for @xmath260 ( b ) and @xmath261 ( c ) . ( d ) shows linecuts for fixed @xmath262 with @xmath263 in black solid , red dot - dashed , blue long dashed and green short dashed lines , respectively . @xmath173 and @xmath174 . [ fig : intrinsic spectra],title="fig:",width=158 ]   ( a ) and @xmath180 for @xmath260 ( b ) and @xmath261 ( c ) . ( d ) shows linecuts for fixed @xmath262 with @xmath263 in black solid , red dot - dashed , blue long dashed and green short dashed lines , respectively . @xmath173 and @xmath174 . [ fig : intrinsic spectra],title="fig:",width=158 ]   ( a ) and @xmath180 for @xmath260 ( b ) and @xmath261 ( c ) . ( d ) shows linecuts for fixed @xmath262 with @xmath263 in black solid , red dot - dashed , blue long dashed and green short dashed lines , respectively . @xmath173 and @xmath174 . [ fig : intrinsic spectra],title="fig:",width=158 ]    for @xmath173 and @xmath174 . due to energy conservation , there are no excitations beneath a frequency @xmath264 . in contrast to graphene , where one observes a diverging behaviour of the polarization at @xmath264 , here @xmath96 increases continously from @xmath62 . this is due to the broken particle - hole symmetry ( @xmath265 ) which ensures that the lowest energy excitations correspond to processes exciting particles from the valence band to the proximity of the dirac point , where , however , the density of states is zero . the excitation spectrum shows a maximum for small momenta @xmath202 which lies beneath the plasmon dispersion given by the black line , perturbatively calculated from eq . ( [ eq : plasmon equation ] ) up to order @xmath266 for @xmath260 .    considering a finite coulomb interaction , the excitation spectrum is given by @xmath180 plotted in fig . [ fig : intrinsic spectra ] ( b ) for @xmath260 and ( c ) for @xmath261 . the maximum of the spectrum shifts to higher energies compared to the non - interacting one , indicating the formation of a collective excitation in the system , i.e. a plasmon . this is proven by solving the plasmon equation ( eq . ( [ eq : plasmon equation ] ) ) perturbatively up to order @xmath267 , with the dispersion plotted as a black line on top of the spectrum . additionally , the dispersion based on the expansion of @xmath239 in the limit @xmath197 for @xmath268 and @xmath148 are plotted as gray lines in fig . [ fig : intrinsic spectra](c ) . the plasmon dispersion relation starts linearly for small @xmath192 , as one would expect for a neutral system without doping . at high energies on the other hand , a free - particle behaviour could be expected , leading to the usual @xmath269 dispersion known from doped systems . although eq . ( [ eq : pi undoped expanded ] ) shows that this picture is only partly true due to the logarithmic correction of @xmath93 , fig . [ fig : intrinsic spectra ] ( c ) indicates a qualitative agreement . [ fig : intrinsic spectra ] ( d ) shows linecuts of @xmath180 for fixed @xmath262 with @xmath263 . additionally , the black vertical lines indicate the plasmon frequency for @xmath270 ( left line ) and @xmath261 ( right line ) . for @xmath271 the maximum of the interacting spectrum lies between the maximum of the non - interacting spectrum and the plasmon frequency , indicating that single - particle and collective excitations are equally strong . increasing the interaction to @xmath260 , the maximum of the interacting spectrum and the plasmon frequency almost coincide , therefore the plasmon dominates over the single - particle excitation . at very large interactions @xmath261 , the plasmon is the only relevant excitation in the system . increasing the coulomb interaction broadens the plasmon peak and reduces its height as shown in fig . [ fig : intrinsic spectra ] ( d ) . this seems contrary to the picture of a plasmon as a sharp interaction - induced charge resonance , suggesting that these interband plasmons may not be well - defined for high energies . yet this is a false conclusion . in sec . [ f - sum rule interpretation exc spectrum ] we discussed that the contribution of the resonance to the f - sum rule is the actual measure of importance of a resonance . it can be estimated by multiplying the peak height in @xmath180 by @xmath183 , while the relevant peak width is given by @xmath186 . the latter is decreasing with @xmath183 according to eq . ( [ eq : damping expansion ] ) . from this normalization of the peak we conclude that the discussed interband plasmons fulfill the interpretation as sharp interaction - induced charge resonances , with the width @xmath186 decreasing with increasing plasmon frequency , above the critical frequency @xmath272 as defined in sec . [ sec : expansion of pi ] . the f - sum rule provides a check for our numerics . in fig . [ fig : f - sum numerical undoped ]    , with @xmath273 being the difference between the numercial and analytical f - sum rule . black dots are for the non - interacting spectrum , while blue stars stand for @xmath260 and green triangles for @xmath261 . @xmath274 , @xmath175 , @xmath173 and @xmath174 . the deviations around @xmath275 stem from numerical instabilities , which are however negligibly small . [ fig : f - sum numerical undoped],width=226 ]    we plot the ratio @xmath276 , with @xmath277 where @xmath278 is the numerical calculated f - sum rule and @xmath279 the analytic one . the deviation are of the order @xmath280 , comparable to the analytical uncertainty , see sec . [ f - sum rule introduction ] , and thus negligible . the f - sum has to be the same for interacting and not interacting systems . we find a slight dependence on the interaction strength @xmath70 , which could be a numerical artifact , depending on @xmath146 , or a real @xmath70 dependence like in graphene , where spectral weight is missing for small frequencies , cf . ( 14 ) in ref . ( @xmath281 for the undoped dirac model ) . as the effect declines with increasing cutoff @xmath146 , we conclude that the rpa approximation in the bhz model misses no spectral weight compared to the full coulomb interaction , even in the undoped limit . in this section , we extend our analysis to finite doping @xmath282 , where a net charge density is present in the system . doping the system has two effects : one is the fermi blocking of interband excitations ( red arrow in fig . [ fig : ab bandstructure ] ( a ) ) for small @xmath51 and @xmath52 . the other is the appearance of intraband excitations ( green arrow in fig . [ fig : ab bandstructure ] ( a ) ) , which are absent in the intrinsic limit . again , first we study the polarization and screening properties of the system in the static limit , where we also study friedel oscillations due to the scattering on a charged impurity . then , we study the dynamical polarization function in the long - wavelength limit , where we obtain an analytical expression for the collective plasmonic modes of the system . finally we numerically compute the dynamical polarization function in the full range of momenta and frequencies in the full parameter space of the bhz model , analyzing the effect of each of the model parameters . particular emphasis is put on the coexistence of interband and intraband plasmons and on how the bhz model interpolates between the dirac and 2deg behavior .      in fig . [ fig : normalized ] ,     of the bhz model for @xmath283 at finite doping for different value of @xmath284 , normalized by the dos @xmath285 . in the inset , details on the value of @xmath286 as a function of @xmath284 are given . , width=302 ]    we present the static polarization function @xmath287 at finite doping , conveniently normalized by the density of states at the fermi level @xmath285 . this normalization stands out naturally from the long wavelength property of the polarization function @xmath288 for the bhz model at finite doping , @xmath289 has a pronounced dependence on the extrinsic parameter @xmath290 . for @xmath291 ( @xmath292 ) the fermi level falls in a region where locally the dispersion curve has predominant dirac ( 2deg ) character . in a 2deg system , the static polarization assumes the following analytic form  @xcite @xmath293 while in the dirac limit we have  @xcite @xmath294.\end{aligned}\ ] ] our calculations for the bhz model with @xmath283 correctly reproduce the dirac and 2deg limits for @xmath291 and @xmath292 , respectively . we note that with a finite @xmath16 term and nonzero @xmath284 the polarization will always have a decay behavior for @xmath295 . in the 2deg and dirac limit one finds @xmath296 for @xmath297 , coincidence due to the balancing effect of dispersion curve and overlap factor . interestingly , in the bhz model we observe instead a deviation from unity , shown in details in the inset of fig . [ fig : normalized ] , which has a maximum for @xmath298 . in the 2deg limit , @xmath299 , @xmath300 has a strong discontinuity in its first derivative at @xmath301 , while for decreasing @xmath284 this discontinuity decreases and finally vanishes in the dirac limit , where the discontinuity affects only the second derivative . we already analyzed in section  [ sec : screening_int ] the intrinsic response of a bhz system to a test charge , when no net charge density is present in the system . while the intrinsic response is realized on intrinsic scales of the model @xmath232 , the metallic response ( at finite electronic density ) is characterized by the fermi wave length @xmath302 . therefore it is convenient to express @xmath114 as a function of dimensionless units @xmath303 , due to the presence of a discontinuity at @xmath304 . the induced charge density @xmath305 is given by @xmath306 where , using the property eq . ( [ eq : lwl_pi ] ) , we have emphasized the dependence of the induced density on the dos at the fermi level , which now appears in the scaling factor @xmath307 . we note that the integral also depends on the parameters @xmath70 and @xmath284 ( and naturally on @xmath167 and @xmath151 , when finite ) .     and @xmath308 , @xmath252 and @xmath309 . data in different panels belong to systems with @xmath310 , @xmath252 , @xmath309 , @xmath311 , @xmath312 and @xmath313 . all calculation are obtained by keeping @xmath314 and @xmath315  nm@xmath316 , and varying @xmath317 , @xmath318 and @xmath319  ev  nm , for @xmath308 , @xmath252 and @xmath309 , respectively , while varying the parameter @xmath16 accordingly to @xmath284 . in the panel @xmath320 ( 2deg limit ) , the three curves with @xmath308 , @xmath252 and @xmath309 are quite close and correspond to similar @xmath321 parameter ( @xmath322 , @xmath323 and @xmath324 , respectively ) . , width=321 ]    in fig . [ fig : n_dop ] , we present the induced screening electronic radial density for the bhz model for @xmath283 due to a point - like test charge . each panel corresponds to a different value of the ratio @xmath290 , and within each panel curves differing by the dirac fine structure constant @xmath70 are presented . friedel oscillations appear of period @xmath302 , which become more defined for larger @xmath70 . we also note that density oscillations are more prominent for @xmath325 than in the dirac ( @xmath326 ) and 2deg limits ( @xmath327 ) . in the 2deg limit the @xmath70 parameter is ill defined and should be replaced by the more general parameter @xmath328 , characterizing the dielectric response of the system . the presence of friedel oscillations and their asymptotic behavior are related through the lighthill theorem  @xcite to discontinuities in the static polarization function and its derivatives ( see for example ref .   for a detailed discussion ) . a discontinuity like @xmath329 in @xmath330 , with @xmath331 the heaviside step function and @xmath332 , translates into a decay of the oscillations in @xmath333 with leading order @xmath334 . one finds @xmath335 ( @xmath336 ) for the leading order discontinuity of a 2deg ( dirac ) system , such that the first ( second ) and all higher derivatives of the static polarization function are discontinuous at @xmath337 . analyzing the friedel oscillations for the bhz model , one finds a composition of two different contributions with an asymptotic decay at large distances as @xmath234 ( 2deg contribution ) and @xmath338 ( dirac contribution ) , respectively . as a consequence , the discontinuity in the rpa polarization function of the bhz model at @xmath339 can be very well approximated by a combination of 2deg ( @xmath335 ) and dirac ( @xmath336 ) contribution . in the dirac ( 2deg ) limit , the effect of the discontinuity in the second ( first ) derivative becomes predominant and oscillations purely decay in leading order as @xmath338 ( @xmath234 ) .      at finite doping , for small momenta @xmath51 , the polarization function is governed by intraband excitations , as the interband excitations are fermi - blocked . we perform an expansion in this limit , for @xmath240 , to gain an analytical insight into the physics at finite doping and derive an analytical formula for the plasmon dispersion . in particular , intraband plasmons are expected to be the dominant excitation for small momenta , similarly to the 2deg and dirac case . we expand the polarization function up to order @xmath340 @xmath341= \pi_{44}\frac{x^{4}}{\omega^{4}}+\pi_{42}\frac{x^{4}}{\omega^{2}}+\pi_{40}x^{4}\\ + \pi_{22}\frac{x^{2 } } { \omega^{2 } } + \pi_{20}x^{2}+\mathcal{o}\left(\omega^{2}\right)\end{aligned}\ ] ] and use it to solve eq . ( [ eq : plasmon equation ] ) . we obtain a plasmon dispersion of @xmath342 with the leading coefficient @xmath343 with @xmath144 $ ] and @xmath344 . in the limit of zero mass , @xmath345 interpolates smoothly between @xmath62 for @xmath346 and @xmath347 for @xmath348 . the former case corresponds to the dirac limit , where one finds the plasmon frequency @xmath349 in the literature  @xcite , being identical to eqs . ( [ eq : plas_freq_bhz],[eq : pi doped expanded dirac limit ] ) . the latter case is the 2deg limit , where one finds the plasmon dispersion @xmath350 in the literature  @xcite , with @xmath351 the carrier density and @xmath352 . this in in agreement with eqs . ( [ eq : plas_freq_bhz],[eq : pi doped expanded 2deg limit ] ) . thus the bhz model as a function of its parameters reproduces the plasmon dispersion in the dirac and 2deg limits and interpolates between them . we note that for @xmath353 the term @xmath354 is zero and the intraband plasmon disappears . in this limit , the leading order contribution @xmath355 of the intrinsic polarization , eq . [ eq : pi undoped expanded ] , takes the place of @xmath354 . the crucial difference between the extrinsic and the intrinsic polarization is that the latter has a finite imaginary part of order @xmath355 , leading to the linear dispersion of the interband plasmons . yet for finite @xmath356 , these interband plasmons are supressed due to the fermi blockade of the interband excitations and only exist if their plasmon frequency exceeds both the chemical potential @xmath357 and the critical frequency @xmath272 as defined in sec . [ sec : expansion of pi ] , see for example fig . [ fig : im epsi bhz xf01 ] .    besides the different scaling with momenta in the limit @xmath358 , also the scaling with @xmath70 is different for the inter- and intraband plasmons , eqs . ( [ eq : plasmonfrequency linear ] ) and ( [ eq : plas_freq_bhz ] ) : linear vs. square root . this will have important consequences in the following when we will discuss how to separate the two different collective excitations . we begin the discussion of the doped spectrum by looking at the limiting results of 2deg and dirac system . from this , we then find that we can interpolate between them by changing the fermi momentum . interestingly , by considering the cases of broken particle - hole symmetry and large masses , we also find regimes which are distinct from the dirac and 2deg limit . as an example , these regimes support both inter- and intraband plasmons at parameters which are realistic for hgte qws .    in all the following plots , the boundaries of the single - particle spectrum will be indicated by faint black lines , the isolines @xmath95 by red lines . the plasmon dispersions are plotted as black curves ( full result from perturbation theory ) and gray curves ( expanded result in limit @xmath197 ) .      in the 2deg limit , only intraband excitations are possible . the polarization function has a well - known analytical form  @xcite , therefore we can easily plot the non - interacting spectrum in fig .  [ fig : pi 2deg ] ( a ) .     with @xmath359 and @xmath360 the degeneracy factor . ( b ) @xmath180 for @xmath361 , with @xmath362 the coulomb interaction . we add an artificial damping in the region of @xmath363 to make the plasmons visible . [ fig : pi 2deg],title="fig:",width=158 ]   with @xmath359 and @xmath360 the degeneracy factor . ( b ) @xmath180 for @xmath361 , with @xmath362 the coulomb interaction . we add an artificial damping in the region of @xmath363 to make the plasmons visible . [ fig : pi 2deg],title="fig:",width=158 ]    @xmath96 is peaked for @xmath364 closely to the upper boundary of the spectrum . it decays to zero instead for large momenta and frequencies like @xmath365 , if one considers a fixed ratio @xmath366 within the spe region .    the interacting spectrum is shown in fig .  [ fig : pi 2deg ] ( b ) . an intraband plasmon appears with the usual @xmath269 dispersion for @xmath358 . it absorbes all of the spectral weight in this limit , thus @xmath180 is suppressed in the spe region . for intermediate momenta , the plasmon dispersion lies in the spe region and the plasmon decays and broadens . for larger momenta and frequencies , the interacting and non - interacting spectra agree qualitatively . the dirac spectrum comprises both inter- and intraband excitations . the polarization function still has a well - known analytical expression  @xcite , of which we plot the non - interacting spectrum @xmath96 in fig . [ fig : pi graphene ] ( a ) .     with @xmath367 and @xmath360 the degeneracy factor . in the gray area , the colorscale is exceeded due to the divergency of @xmath96 . ( b ) @xmath180 for @xmath368 , with @xmath362 the coulomb interaction . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : pi graphene],title="fig:",width=158 ]   with @xmath367 and @xmath360 the degeneracy factor . in the gray area , the colorscale is exceeded due to the divergency of @xmath96 . ( b ) @xmath180 for @xmath368 , with @xmath362 the coulomb interaction . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : pi graphene],title="fig:",width=158 ]    the intraband excitations occur for higher energies @xmath369 , while for intraband excitations less energy is needed , @xmath370 . both excitation spectra touch at @xmath371 , where they diverge . only the fermi - blockade suppresses the interband transitions in @xmath96 for @xmath372 and cures the divergency , see fig . [ fig : pi graphene ] ( a ) for @xmath369 . one finds a @xmath373 decay for high frequencies .    the interacting spectrum @xmath180 is plotted in fig .  [ fig : pi graphene ] ( b ) for @xmath374 . similar to the 2deg , all of the intraband spectral weight is absorbed by a plasmon in the limit @xmath358 and the divergence at @xmath371 is cured . interestingly , for sufficient large interaction strength @xmath70 the plasmon decays in the interband spectrum . for larger momenta and frequencies , we note that the intraband polarization does not recover the non - interacting value , as it does for the 2deg , but remains much smaller . therefore single - particle intraband excitations are blocked altogether for all momenta and frequencies in this limit . the missing spectral weight goes into a charge resonance at higher frequencies in the interband spectrum  @xcite . yet , this resonance is not a solution of the plasmon equation and therefore not a plasmon  @xcite . the bandstructure of the bhz model without mass and particle - hole symmetry breaking is shown in fig . [ fig : ab bandstructure ] ( a ) .    . interband spectrum in red , intraband spectrum in blue and mixed area in purple . [ fig : ab bandstructure],title="fig:",width=158 ] . interband spectrum in red , intraband spectrum in blue and mixed area in purple . [ fig : ab bandstructure],title="fig:",width=158 ] . interband spectrum in red , intraband spectrum in blue and mixed area in purple . [ fig : ab bandstructure],title="fig:",width=158 ] . interband spectrum in red , intraband spectrum in blue and mixed area in purple . [ fig : ab bandstructure],title="fig:",width=158 ]    the interband single - particle excitations lying lowest in energy are symmetric in momentum as shown by the red arrow in fig . [ fig : ab bandstructure ] ( a ) , going from @xmath376 to @xmath377 . due to particle - hole symmetry , this leads to nesting and thus one expects these excitations to dominate the interband spectrum . interband excitations as indicated by the dashed , black arrow on the other hand , going from momentum @xmath378 to @xmath379 with @xmath380 , are suppressed due to imperfect nesting of the different sized electron and hole cones , as well as by a small overlap factor . the latter can be cured by introducing a large negative mass , as will be shown in sec . [ sec_bhz_mass ] . then these excitations have a considerable influence onto the polarization for small energies , helping with the formation of interband plasmons , following the ideas presented in sec . [ ( anti-)screening and intrinsic plasmons ] . in the pure dirac system , these processes are forbidden by helicity .    by varying the doping level we can modify the excitation spectrum of the system [ see figs . [ fig : ab bandstructure ] ( b)-(d ) ] to resemble that of a dirac system ( @xmath381 ) or of a 2deg ( @xmath382 ) , or to obtain an intermediate behavior ( @xmath383 ) . in the pictures we highlight the boundaries of the excitation spectra , with the red area corresponding to the interband spectrum and the blue area to the intraband spectrum . the overlap between the two is indicated by the purple area . the boundaries of the spectra vary from the linear graphene behavior to the @xmath384 dependence of the 2deg . in general , the mixing of linear and quadratic dispersion leads to an overlap of the inter- and intraband spectrum . this affects the visibility of the interband plasmons , which can be hidden due to strong single - particle damping . [ [ weak - doping - of - x_f0.1 ] ] weak doping of @xmath385 : + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the extrinsic ( @xmath25 ) and intrinsic ( @xmath386 ) scales of the system are separated by one order of magnitude . as the fermi surface lies in the ( almost ) linear part of the spectrum , we expect that on the @xmath25 scale we resemble graphene . the physics on the @xmath386 scale on the other hand should be more or less untouched by the doping , and the system should behave as in the intrinsic limit . we plot @xmath239 in fig . [ fig : imp bhz xf01 ] .    . the red line indicates @xmath95 . [ fig : imp bhz xf01],title="fig:",width=158 ] . the red line indicates @xmath95 . [ fig : imp bhz xf01],title="fig:",width=158 ]    comparing panel ( a ) to fig . [ fig : pi graphene ] , one finds good agreement with the dirac case . the biggest deviation is found in the peak of @xmath96 at @xmath387 , which is not symmetric as for a dirac system due to the overlap of inter- and intraband spectrum [ fig . [ fig : ab bandstructure ] ( b ) ] . the finite quadratic part in the spectrum cures the divergency formerly occuring in the dirac limit . the real part of @xmath388 is strongly negative only at the upper boundary of the intraband spectrum . this indicates that for small interactions , only one plasmon will dominate the excitation spectrum on the fermi scale .    as we are interested in the regime where both inter- and intraband plasmons are visible , we look at the interacting spectrum , given in fig .  [ fig : im epsi bhz xf01 ] by plotting @xmath180 , for a strong interaction @xmath261 .     for @xmath385 with @xmath261 on the @xmath389 scale ( a ) and the @xmath386 scale ( b ) . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : im epsi bhz xf01],title="fig:",width=158 ]   for @xmath385 with @xmath261 on the @xmath389 scale ( a ) and the @xmath386 scale ( b ) . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : im epsi bhz xf01],title="fig:",width=158 ]    on the fermi scale [ panel ( a ) ] , the intraband plasmon absorbs all spectral weight from the intraband spectrum . the dispersion agrees with the perturbative dispersion from the expansion in eq . ( [ eq : plas_freq_bhz ] ) in the limit @xmath197 , plotted as a gray curve . the green , dashed line shows the linear dispersion of the interband plasmon in the undoped limit , based on eq . ( [ eq : plasmonfrequency linear ] ) . on the fermi scale , it is not obvious that there is an interband plasmon , although the interacting polarization function develops a smeared resonance around the perturbative interband plasmon dispersion for high momenta . switching to the intrinsic scale , fig . [ fig : im epsi bhz xf01 ] ( b ) , one finds the interband plasmon , corresponding to the single peak in @xmath390 , unperturbed by doping for momenta much larger than @xmath389 . the dispersion is the same as for a plasmon in the undoped limit  @xcite . the two black lines near the peak are just the boundaries of the intraband excitation spectrum , which does not play a role here .    as in the limit of @xmath197 the interband plasmon dispersion scales linearly with @xmath70 , @xmath391 see eq . ( [ eq : plasmonfrequency linear ] ) , while the intraband plasmon frequency is proportional to @xmath392 , lowering the interaction strength will lead to an overlap of the two resonances below some critical @xmath70 . [ [ strong - doping - of - x_f3 ] ] strong doping of @xmath393 : + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    strong doping of the system significantly increases the spectral weight , as shown in fig . [ fig : f - sum rule ] ( b ) , with the increase of intraband excitations , while most of the interband excitations are fermi - blocked , leading therefore to an effective decoupling of the two bands . we expect the overall spectrum to be governed by intraband excitations and to resemble the spectrum of a 2deg , as the fermi surface lies in the ( almost ) quadratic part of the spectrum . the corresponding @xmath239 is plotted in fig . [ fig : imp bhz xf3 ] .    . the red line indicates @xmath95 . [ fig : imp bhz xf3],title="fig:",width=158 ] . the red line indicates @xmath95 . [ fig : imp bhz xf3],title="fig:",width=158 ]    the single - particle spectrum in panel ( a ) is peaked at small momenta and at energies close to the upper bound of the intraband spectrum . the interband part of the spectrum leads only to minor deviations from the 2deg case [ compare with fig .  [ fig : pi 2deg ] ( a ) ] . the real part of @xmath388 in panel ( b ) is strongly negative at the upper boundary of the intraband spectrum , indicating that only a single intraband plasmon will dominate the interacting spectrum . we additionally note that the static limit property for which the polarization is a constant @xmath394 for @xmath395 , discussed in sec . [ static limit ] , extends also to an area of finite @xmath52 .    the interacting spectrum is shown in fig .  [ fig : im epsi bhz xf3 ] by plotting @xmath180 for the interaction strength @xmath261 .     for @xmath393 and @xmath261 . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : im epsi bhz xf3],width=158 ]    even for this large coulomb interaction , we only find the intraband plasmon . this is as expected due to the combined effects of fermi blocking of interband excitations and increased spectral weight for intraband transitions . the interband plasmon lies in the large overlap of inter- and intraband spectrum , cf . [ fig : im epsi bhz xf01 ] ( b ) , and it is therefore heavily damped and not visible in the overall spectrum .    [ [ intermediate - doping - of - x_f1 ] ] intermediate doping of @xmath396 : + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    for intermediate doping levels like @xmath396 , a mixture of dirac and 2deg behaviour is expected , due to the similar importance of inter- and intraband excitations . we plot the polarization function @xmath239 in fig . [ fig : imp bhz xf1 ] .    . the red line indicates @xmath95 . [ fig : imp bhz xf1],title="fig:",width=158 ] . the red line indicates @xmath95 . [ fig : imp bhz xf1],title="fig:",width=158 ]    indeed , the single - particle spectrum in panel ( a ) looks like a combination of figs . [ fig : imp bhz xf01 ] ( a ) and [ fig : imp bhz xf3 ] ( a ) . while the the shape of the polarization resembles the one of the 2deg , the interband spectrum is now more pronounced and even dominating for @xmath397 . therefore we could expect both kinds of excitations giving rise to a plasmon mode . the real part of @xmath388 in panel ( b ) shows again just a single minimum , following the upper boundary of the intraband spectrum . the deviations from the constant behaviour @xmath394 for @xmath395 in the case of intermediate doping , see sec . [ static limit ] , are also found for finite @xmath52 .    the interacting spectrum is shown in fig .  [ fig : im epsi bhz xf1 ] by plotting @xmath180 for an interaction strength @xmath261 .     for @xmath396 with @xmath261 . we add an artificial damping in the regions of @xmath363 to make the plasmons visible . [ fig : im epsi bhz xf1],width=158 ]    it is dominated by a single resonance , lying above the intraband part of the single - particle spectrum . for small momenta , this resonance corresponds to the intraband plasmon . yet for intermediate momenta , a comparison with the interband plasmon dispersion in fig . [ fig : im epsi bhz xf01 ] ( b ) indicates that also the interband plasmon contributes to the resonance . a clear distinction between the two is then not possible anymore .    in summary , doping the system offers the possibility to change the excitations spectrum on the fermi scale from a dirac to a 2deg type . the interacting excitation spectrum is usually governed by a single intraband plasmon , while the interband plasmon is hidden in the single - particle background . only large interaction strengths offer a possibility to see both plasmons in the spectrum . in the following , we will now analyze the influence of both broken p - h symmetry and finite masses , which both offer a way to separate the two plasmons and make them visible in the total spectrum .      a broken particle - hole symmetry , @xmath398 , with small or vanishing mass is the experimental relevant case for hgte qws . it also offers the possibility of blocking the interband spe spectrum close to the minimal excitation energy @xmath387 , resulting in less damped interband plasmons  @xcite .      here , we want to use a similar effect for the intraband excitations in order to separate the inter- and intraband spectrum as well as the two plasmon modes . the broken p - h symmetry introduces an inflection point into the spectrum , @xmath399 , with momentum @xmath400 and energy @xmath401 . for @xmath402 , it lies in the hole part [ @xmath92 ] of the spectrum . with a sufficiently small fermi momentum , @xmath404 , the highest energy intraband excitations involve the dirac point for momenta on the order of the fermi momentum , see fig . [ fig : bhz spectrum finite d ] .    . both low energy interband excitations ( red arrow ) and high energy intraband excitations ( green arrow ) involve the dirac point . [ fig : bhz spectrum finite d],width=207 ]    the same is true for the lowest energy interband excitations . due to the vanishing density of states at the dirac point , both kind of excitations are suppressed , and therefore inter- and intraband spe spectrum are effectively separated in energy and momentum . this situation is shown in fig . [ fig : polarization bhz abkinf12 ] ( a )     ( a ) and @xmath93 ( b ) for @xmath405 and @xmath174 . the red line indicates @xmath95 . [ fig : polarization bhz abkinf12],title="fig:",width=158 ]   ( a ) and @xmath93 ( b ) for @xmath405 and @xmath174 . the red line indicates @xmath95 . [ fig : polarization bhz abkinf12],title="fig:",width=158 ]    and fig . [ fig : im eps bhz abkinf12 ] ( c)-(e ) , where the imaginary part of the polarization goes to zero between inter- and intraband parts of the spectrum , fully separating them . there , we choose @xmath406 , @xmath174 and @xmath173 . this blocking effect holds for small momenta up to roughly @xmath407 , indicated by the black vertical line in fig . [ fig : polarization bhz abkinf12 ] at @xmath408 . for larger momenta the high energy intraband excitation go from deep in the valence band directly to the fermi surface - the blocking effect of the dirac point is gone , see dashed arrow in fig . [ fig : bhz spectrum finite d ] .    in fig . [ fig : polarization bhz abkinf12 ] ( b ) , @xmath93 shows one major difference in comparison to the p - h symmetric case of weak doping in fig . [ fig : imp bhz xf01 ] ( b ) . at the border of intra- and interband spectrum a strong antiscreening region is formed . for sufficiently low @xmath70 a plasmon should exist there , clearly separated from the second antiscreening region at higher @xmath52 , giving rise to the possibility of observing both intra- and interband plasmons . this can be seen in fig . [ fig : im eps bhz abkinf12 ] , where we plot @xmath180     for @xmath405 , @xmath174 , @xmath173 and @xmath409 ( a ) and @xmath260 ( b ) . ( c)-(e ) show linecuts for fixed @xmath410 , @xmath411 and @xmath412 , respectively , with @xmath413 in black solid line , red dot - dashed line , blue long dashed line and green short dashed line , respectively . the black , vertical line separates the inter- and intraband spe region . [ fig : im eps bhz abkinf12],title="fig:",width=158 ]   for @xmath405 , @xmath174 , @xmath173 and @xmath409 ( a ) and @xmath260 ( b ) . ( c)-(e ) show linecuts for fixed @xmath410 , @xmath411 and @xmath412 , respectively , with @xmath413 in black solid line , red dot - dashed line , blue long dashed line and green short dashed line , respectively . the black , vertical line separates the inter- and intraband spe region . [ fig : im eps bhz abkinf12],title="fig:",width=158 ]   for @xmath405 , @xmath174 , @xmath173 and @xmath409 ( a ) and @xmath260 ( b ) . ( c)-(e ) show linecuts for fixed @xmath410 , @xmath411 and @xmath412 , respectively , with @xmath413 in black solid line , red dot - dashed line , blue long dashed line and green short dashed line , respectively . the black , vertical line separates the inter- and intraband spe region . [ fig : im eps bhz abkinf12],title="fig:",width=158 ]   for @xmath405 , @xmath174 , @xmath173 and @xmath409 ( a ) and @xmath260 ( b ) . ( c)-(e ) show linecuts for fixed @xmath410 , @xmath411 and @xmath412 , respectively , with @xmath413 in black solid line , red dot - dashed line , blue long dashed line and green short dashed line , respectively . the black , vertical line separates the inter- and intraband spe region . [ fig : im eps bhz abkinf12],title="fig:",width=158 ]   for @xmath405 , @xmath174 , @xmath173 and @xmath409 ( a ) and @xmath260 ( b ) . ( c)-(e ) show linecuts for fixed @xmath410 , @xmath411 and @xmath412 , respectively , with @xmath413 in black solid line , red dot - dashed line , blue long dashed line and green short dashed line , respectively . the black , vertical line separates the inter- and intraband spe region . [ fig : im eps bhz abkinf12],title="fig:",width=158 ]    for @xmath409 ( a ) and @xmath260 ( b ) . panel ( c ) - ( e ) show line cuts for fixed momenta @xmath414 and different interaction strengths @xmath413 .    for large interaction strength @xmath409 , the intraband plasmon decays into the interband spe spectrum , see panel ( a ) for @xmath415 and the green short dashed line in panel ( c ) . most of the spectral weight stays there also for larger momenta , as @xmath180 is close to 0 in the intraband spe region and the resonance between inter- and intraband spe spectrum is weak . the latter can be best seen in the insets of panel ( d ) and ( e ) , represented by the green short dashed line peaked slightly above [ ( d ) ] or below [ ( e ) ] the black vertical line separating intra- and interband spe region . yet even with the peak being small , it indicates the formation of a slighly damped plasmon , but with small spectral weight . the missing spectral weight is transferred to higher energies into the interband spe region . for intermediate momenta , a second plasmon branch forms , see panel ( a ) for @xmath416 and the second peak of the green short dashed line in the inset of panel ( d ) . for even higher momenta , @xmath417 , it overlaps with the forming interband plasmon leading to a broad charge resonance without clear peak , see green short dashed line in panel ( e ) for @xmath418 . the picture changes for smaller interaction strength . for @xmath260 , the intraband plasmon decays in the region between inter- and intraband spe spectrum , indicated by the strong peak of the blue long dashed line in panel ( c ) . as the single - particle excitations in this region are suppressed due to the dirac point , the plasmon leads to a high and narrow peak of @xmath180 . considering larger momenta @xmath419 , the resonance is split : one part forms an intraband plasmon in the intraband spe region , see blue long dashed line peaked slightly below the black vertical line in panels ( d ) and ( e ) . the second part stays in the interband spe region , where it enhances the spe peak [ black line in the inset of panel ( d ) ] for intermediate momenta @xmath416 . for momenta @xmath420 , an interband plasmon forms , as shown in panel ( e ) . there , the broad single - particle peak [ black line ] around @xmath421 gets reshaped into a clear peaked resonance [ blue long dashed line ] - the interband plasmon . taking the experimental parameters from sec . [ sec : experimental_parameters ] , @xmath119 and @xmath120 , one finds for the plots in fig . [ fig : im eps bhz abkinf12 ] the fermi momentum @xmath422 and chemical potential @xmath423 . the plot range is therefore @xmath424q_{0}=\left[0,0.3\right]\frac{1}{\mathrm{nm}}$ ] and @xmath425\mathrm{\frac{e_{0}}{\hbar}}=\left[0,210\right]\mathrm{thz}$ ] and thus of the right order of magnitude for experimental techniques like raman spectroscopy or electron loss spectroscopy .      both plasmonic resonances in fig . [ fig : im eps bhz abkinf12 ] ( b ) overlap for @xmath415 , before they separate for higher momenta . therefore the question arises whether one can really speak of a clear distinction between inter- and intraband plasmons for larger momenta . here , we want to study the f - sum rule and thus the spectral weight of the different excitations . the relative deviations of numerical to analytical f - sum rule are again of the order @xmath280 and thus negligible . [ fig : f - sum numerical doped ] ( a )     of interacting over non - interacting f - sum rule of the interband excitations . black dots are for @xmath271 , blue stars stand for @xmath260 and green triangles for @xmath409 . ( b ) the same for the ratio @xmath426 of the intraband excitations . for all plots : @xmath427 and @xmath428 . [ fig : f - sum numerical doped],title="fig:",width=158 ]   of interacting over non - interacting f - sum rule of the interband excitations . black dots are for @xmath271 , blue stars stand for @xmath260 and green triangles for @xmath409 . ( b ) the same for the ratio @xmath426 of the intraband excitations . for all plots : @xmath427 and @xmath428 . [ fig : f - sum numerical doped],title="fig:",width=158 ]    shows the ratio of spectral weight in the interband spe region for the interacting over the non - interacting case , @xmath429 , and panel ( b ) the same for the intraband spe region . the intraband plasmon lying between these two regions for @xmath430 is excluded . as for cutoffs @xmath431 one usually has @xmath432 , transfer of spectral weight from one region to the other can lead to quantitatively different relative changes of spectral weight in panels ( a ) and ( b ) . as a key result , we find that there is always spectral weight missing in the intraband spe region . for small momenta , @xmath430 , the weight goes into the undamped intraband plasmon [ this follows directly from the conservation of the f - sum rule for interacting and non - interacting systems ] , while at larger momenta it is transferred to higher energies into the interband spe region . yet , the increase is only about @xmath433 at @xmath420 , such that we can conclude that the plasmon between the inter- and intraband spe region is a pure intraband plasmon with a reduced spectral weight . the plasmon in the interband region is the interband plasmon we know already from the undoped system , see fig . [ fig : intrinsic spectra ] , with a slight increased spectral weight from the intraband spe region .      deviations in the thickness of a hg(cd)te qw lead to the opening of a small gap in the bandstructure , resulting in a topological trivial @xmath435 or non - trivial @xmath436 system . apart from the possible appearance of edge states , which is beyond the scope of this paper , a small mass works in opposition to the blocking effect of finite @xmath54 , as it generates a finite density of states for @xmath437 . in the following , we therefore show that the blocking effect of a finite @xmath54 is robust against the opening of small gaps .    in fig . [ fig : polarization bhz mabkinf004 ]     ( a ) and the interacting one @xmath180 ( b ) with @xmath260 . @xmath438 , @xmath439 and @xmath174 . [ fig : polarization bhz mabkinf004],title="fig:",width=158 ]   ( a ) and the interacting one @xmath180 ( b ) with @xmath260 . @xmath438 , @xmath439 and @xmath174 . [ fig : polarization bhz mabkinf004],title="fig:",width=158 ]    we plot the non - interacting and interacting spectrum for @xmath174 and a small mass @xmath440 . a comparison with figs . [ fig : polarization bhz abkinf12 ] and [ fig : im eps bhz abkinf12 ] shows that the small mass has just the effect of separating the inter- and intraband spe region additionally . thus we conclude that the idea of observing both plasmons in experiments is robust against slight deviations in the mass and therefore the thickness of the hg(cd)te qw . a finite dirac mass opens a gap in the bandstructure and changes the pseudospin , and therefore the overlap factor , in a non - trivial fashion . thus we can expect in general a quite different behavior for positive and negative mass . yet , for these differences to occur on the intrinsic scale and thus influence the interband plasmons , @xmath441 should be of the order of 1 . in the following , we study such large masses , both negative and positive , with p - h symmetry . while not experimentally relevant for hgte qws , it offers the possibility to study the effect of a topological bandstructure on the electronic excitations , including plasmons . we also note here that the dispersion of the bhz model becomes purely parabolic for the mass @xmath442 : @xmath443 . in this limit , the polarization function , eq . ( [ eq : pi_par ] ) , can be calculated analytically .      for the parameters @xmath444 and @xmath60 we plot the polarization @xmath239 in fig . [ fig : im epsi bhz mnab xf033 ] ( a ) and ( b ) .    . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath60 and @xmath444 in all plots . [ fig : im epsi bhz mnab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath60 and @xmath444 in all plots . [ fig : im epsi bhz mnab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath60 and @xmath444 in all plots . [ fig : im epsi bhz mnab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath60 and @xmath444 in all plots . [ fig : im epsi bhz mnab xf033],title="fig:",width=158 ]    the mass separates intra- and interband spe regions for @xmath445 . compared to the massless cases of @xmath446 , fig . [ fig : imp bhz xf01 ] , and @xmath447 , fig . [ fig : imp bhz xf1 ] , the interband spe spectrum is enhanced due to the combination of enhanced overlap factor and low doping , thus small fermi blockade . due to the flat bandstructure , even for @xmath448 the chemical potential is just barely above the gap . an interesting consequence of this strong interband transition can be seen in panel ( b ) , where we find two distinct areas where @xmath449 becomes negative . as a consequence , inter- and intraband plasmons will always be separated , with the intraband plasmon being confined to low energies . this stems from the fact that the electrons in the conduction band are pseudo - spin polarized , such that intraband excitations to much higher momenta and energies , where the pseudo - spin shows in the opposite direction , are not possible . this is confirmed in panels ( c ) and ( d ) , where we plot @xmath180 with @xmath261 and @xmath409 , respectively . all the spectral weight of the intraband spe region goes into the plasmon , which at least for @xmath409 follows very well the @xmath258 law . the interband spectrum is dominated by the interband plasmon , having of course a much broader peak due to damping ( finite @xmath450 $ ] ) .    the dashed , green line in the interband spectrum in fig . [ fig : im epsi bhz mnab xf033 ] indicates the energy , at which excitation processes going from momentum @xmath378 to @xmath379 with @xmath380 are possible , see black , dashed arrow in fig . [ fig : ab bandstructure ] ( a ) . usually suppressed by the overlap factor , a large negative mass enhances the overlap of these excitations to near unity for small fermi momenta . [ fig : im epsi bhz mnab xf033 ] ( b ) and ( d ) show that the interband plasmons mainly occur above this line , indicating that the described excitation process is important for the collective excitation . as the process is forbidden by helicity in the pure dirac system , it is one reason why the bhz model supports intrinsic plasmons while the dirac model does not .      for the parameters @xmath444 and @xmath451 we plot the polarization in fig . [ fig : im epsi bhz mab xf033 ] ( a ) and ( b ) .    . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath451 and @xmath444 in all plots . [ fig : im epsi bhz mab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath451 and @xmath444 in all plots . [ fig : im epsi bhz mab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath451 and @xmath444 in all plots . [ fig : im epsi bhz mab xf033],title="fig:",width=158 ] . ( c ) and ( d ) show @xmath180 with @xmath261 and @xmath409 , respectively . @xmath451 and @xmath444 in all plots . [ fig : im epsi bhz mab xf033],title="fig:",width=158 ]    compared to the negative mass , the interband spectrum is much weaker . this is a result of the lower overlap factor and the higher chemical potential ( the bandstructure is not as flat as in the ti phase ) , leading to a stronger fermi blockade . for the real part of the polarization , this has the effect that the two former distinct areas of sign reversal now almost merge . the interband excitations are so weak that the minimum @xmath449 always lies closely above the intraband spe region - indicating that it is the main source for plasmons .    in panel ( c ) for @xmath261 , one can identify both inter- and intraband plasmon . interestingly , the polarization is clearly higher in the pure interband spe region than in the mixed inter- and intraband spe spectrum , suggesting that the latter one serve as an additional damping for the interband plasmon . going to the smaller interaction strength @xmath409 in panel ( d ) , one finds just a single resonance following the upper boundary of the intraband spe spectrum . thus we conclude that the interacting spectrum for moderate interaction strength is governed by just intraband plasmons . the interband excitations are too weak to support an additional plasmon but for very high interactions - a consequence of the effective decoupling of the bands by the overlap factor . we have analyzed the dynamical and static polarization properties in random phase approximation of hg(cd)te quantum wells described by the bernevig - hughes - zhang ( bhz ) model . in the static undoped limit , due to the presence of quadratic terms in the bhz model and hence to the natural length scale @xmath452 , the induced charge density in response to a test charge has a finite spatial extent . this is in contrast to the point - like screening charge obtained with the continuous dirac model of graphene . in the doped regime , we have observed friedel oscillations with an intermediate decay behavior between the dirac ( @xmath338 ) and the 2deg ( @xmath234 ) cases . the discussion of the full dynamical polarization function has been focused on the appearance of new interband plasmons due to the interplay of dirac and schrdinger physics . in principle , we expect these plasmons to appear in multiband systems where the imaginary part of the polarization function decays faster with energy than the one in the dirac case ( @xmath373 ) , which is the case for the bhz model ( decay as @xmath453 ) . these plasmons appear already in the undoped system at experimentally relevant parameters , but it is also possible to observe them in the doped regime , where they coexists with the usual intraband plasmons . this is favored by broken particle hole symmetry in the bhz model , which allows for the presence of both a dirac point and an inflection point in the bandstructure . the behavior of these two collective modes is also influenced by the topology of the bandstructure . indeed the two plasmons tend to merge into one another in a gapped trivial insulator , while they remain distinct resonances in the topological insulator phase . we have shown that these new plasmons should appear for momenta and energies on the right order of magnitude for experimental techniques like raman spectroscopy or electron loss spectroscopy on hg(cd)te quantum wells . the wide range of parameters considered in this paper , including the regime of topological trivial and non - trivial insulators , should make our results applicable to all kinds of materials described by phenomenological models interpolating between dirac and schrdinger fermion physics .    throughout this article , we have only discussed bulk excitations of this peculiar two - dimensional system . hence , we have totally ignored the influence of edge states in the topologically non - trivial regime of the model in the presence of physical boundaries . an extension of our analysis to finite size systems might yield exciting new physics , where we expect an interplay of one - dimensional and two - dimensional collective charge excitations . we acknowledge interesting discussions with e. hankiewicz , m. polini , t. stauber and financial support by the dfg ( spp1666 and the dfg - jst research unit _ topotronics _ ) , the helmholtz foundation ( viti ) , and the enb graduate school on topological insulators .
Abstract:

Answer: under certain conditions , hg(cd)te quantum wells ( qws ) are known to realize a time - reversal symmetric , two - dimensional topological insulator phase . its low - energy excitations are well - described by the phenomenological bernevig - hughes - zhang ( bhz ) model that interpolates between schrdinger and dirac fermion physics . we study the polarization function of this model in random phase approximation ( rpa ) in the intrinsic limit and at finite doping . while the polarization properties in rpa of dirac and schrdinger particles are two comprehensively studied problems , our analysis of the bhz model bridges the gap between these two limits , shedding light on systems with intermediate properties . we gain insight into the screening properties of the system and on its characteristic plasma oscillations . interestingly , we discover two different kinds of plasmons that are related to the presence of intra- and interband excitations . observable signatures of these plasmons are carefully analyzed in a variety of distinct parameter regimes , including the experimentally relevant ones for hg(cd)te qws . we conclude that the discovered plasmons are observable by raman or electron loss spectroscopy .