when matter is strongly coupled to light , the interaction can not simply be thought of in terms of absorption and emission processes .
instead we must consider the eigenstates of the fully coupled matter - light system .
the paradigmatic example of this is the existence of exciton polaritons , hybrid matter - light particles formed by the strong interaction between excitons and photons  @xcite .
matter - light coupling can be engineered by confining light in optical cavities , so as to modify the density of states and the coupling to matter . for weak coupling , or a bad cavity , cavity losses
are fast so one can eliminate virtual processes where photons are in the cavity .
this gives fermi s golden rule , but with the cavity density of states modifying the emission rate , as first discussed by purcell  @xcite .
when coupling is strong , first order perturbation theory ( i.e.  fermi s golden rule ) fails , as there can instead be coherent emission and reabsorption of photons before light leaks out of the cavity  @xcite .
a natural context in which strong matter - light coupling arises is between organic molecules and light in semiconductor microcavities . because of the existence of conjugated @xmath0 bonds in organic molecules
, electronic transitions can acquire large dipole moments  @xcite , leading to very strong coupling to light .
when such molecules are placed in optical microcavities this leads to huge polariton splittings  @xcite .
these scales allow such experiments to be performed at room temperature , whereas for many inorganic materials , cryogenic temperatures are required .
the polariton splitting is due to a collective phenomenon : the electronic transitions of many molecules couple to radiation , and as such the polariton splitting grows as the square root of the molecule density .
in contrast , in weak coupling , the rate at which one molecule emits is independent of whether or not any other molecules are present  @xcite .
much of the recent work on organic microcavity polaritons ( see , e.g.ref .
@xcite for a recent review ) has been focused on condensation and lasing  @xcite , involving a strongly pumped system , and the appearance of macroscopic quantum coherence .
there has however also been significant recent work on the effects of matter - light coupling in the vacuum state , i.e.  without strong pumping .
such work aims to understand how the physical and chemical properties of organic molecules are affected by strong coupling to electromagnetic modes .
examples of this include modifying the rates of photochemical reactions  @xcite , or modifying the transport properties of organic semiconductors  @xcite .
more recently , there has also been experimental  @xcite and theoretical  @xcite work on coupling the vibrational state of organic molecules to infra - red radiation , leading to molecular optomechanics .
theoretical work  @xcite has also studied how strong matter - light coupling to electronic states can suppress the effects of disorder and vibronic features in the polariton spectrum .
of particular interest for the present paper is a recent work from the ebbesen group  @xcite , in which the optical spectra of strongly - coupled organic microcavities were studied by varying molecular concentration and temperature , and paying particular attention to the relative weights of the resonant features in the absorption spectra : the two polariton peaks , and a third peak at the bare excitonic energy  @xcite .
our aim in this manuscript is to examine the behavior of such strongly - coupled organic microcavities starting from various microscopic models , allowing quantitative predictions of the extent to which a self - consistent adaptation of the molecular state , driven by coupling with light , may occur .    to understand the variation of the optical spectra with both concentration and temperature the models which we consider all contain a variable degree of coupling to light .
this is because a molecule that has strictly zero coupling to light is not visible in the absorption spectrum , while molecules with a small but non - zero coupling will lead to absorption at the bare molecular energy . in order that the coupling to light can vary self consistently ( in response to the rabi splitting ) , it must depend on some adaptable feature of the state or environment of the molecule .
i.e.  there must be some physical property that can vary , which determines the strength of matter - light coupling .
we refer to this concept hereafter as `` self - consistent molecular adaptation '' .
we refer to this process as `` self - consistent '' because the effective matter - light coupling depends on ( some aspect of ) the molecular state , and the molecular state is modified because of how its energy depends on the matter - light coupling . in the first part of this manuscript we investigate in detail two candidates that could lead to self - consistent adaptation : rotational and vibrational degrees of freedom , we also consider an extension of these models to generic ( classical ) aspects of the molecules physical or chemical state . while we do find a temperature dependence of the optical spectra , the involved energy scales turn out to be incompatible with the observation of ref .  @xcite . in the final part of the present work ,
we examine how our model of vibrationally dressed polaritons naturally predicts an effect whose energy scale is of the right magnitude to explain the data . this effect does not involve the renormalization of the coupling strength , but instead involves the effect of vibrational replicas and their coupling to the excitonic transition on the optical spectra .
our results thus show a first example of the rich , and presently poorly understood behavior that can stem from the interplay of strong matter - light coupling with strong coupling to vibrational / conformational modes of the molecules .
we start by noting that the existence of a peak at the bare energy of the exciton , brought forward as evidence of novel physics in ref .
@xcite , is not unexpected ; such a `` residual excitonic peak '' has been seen in many cases , for example ref .
@xcite discussed theoretically , and demonstrated experimentally , the appearance of such a feature in a gaas / algaas heterostructure , containing quantum wells inside a dbr microcavity .
while such a peak comes from the spectral weight of the exciton line , it is important to note that this peak _ can not _ be viewed simply as excitons which do not couple to light : if they did not couple , they would not be visible in the absorption or transmission spectrum .
the origin of the peak can be understood physically as coming from the subradiant excitonic states due to inhomogeneous broadening . in a disordered system
, the coupling to the photon mode picks out a specific superradiant state , which forms the polaritons , while the other states  orthogonal to this superradiant state  remain at the bare exciton energies .
however , because of energetic disorder , the superradiant state is not an energy eigenstate , and a residual coupling between the superradiant and subradiant states exists , so that the spectral weight of the subradiant states is visible in the optical spectrum  @xcite .
a simple analytical treatment shows that the weight of this residual excitonic peak does decrease as the matter - light coupling increases .
the existence of this peak is thus consistent with the behavior seen in  @xcite .
however , on its own this explanation can not account for the temperature dependence observed , as the residual excitonic peak should be unaffected by temperature , unless @xmath1 approaches optical energies , of the order of @xmath2 ( for comparison 300k@xmath3mev ) .
one of the main goals of this paper is to address how this temperature dependence may occur .
as will become clear in the following , to discuss the microscopic theory of such effects , it will be crucial to consider the physics of ultrastrong matter - light coupling  @xcite , and the breakdown of the rotating wave approximation ( rwa ) .
this requires retaining `` counter - rotating '' terms in the matter - light coupling hamiltonian .
these terms , which involve simultaneous creation of pairs of excitations , are typically considered to be non - resonant and so are often neglected . however ,
if the matter - light coupling is a significant fraction of the bare exciton and photon energies , then these terms have a non - negligible impact .
such behavior has been seen in both inorganic  @xcite and organic  @xcite systems , with a current record of a coupling strength @xmath4 of the bare oscillator frequency  @xcite .
our focus in this paper is on the more typical regime where such counter - rotating terms can not be neglected , but remain sufficiently small to be treated perturbatively .
the rest of the paper is structured in two main sections . in section  [ sec : rotational - freedom ] we consider how temperature dependence can arise due to self - consistent adaptation of the rotational and vibrational degrees of freedom of the molecules , via a mechanism very similar to that proposed in ref .
@xcite . for the orientational degree of freedom ,
we consider both free molecules , and molecules with randomly pinned orientations as appropriate in a polymer matrix .
we will see that in such systems we do predict a temperature dependence of the residual excitonic peak .
however , while this effect could potentially be observed in other experimental realizations , the energy scales ( temperatures ) required and the scaling with molecular concentration are not compatible with the experimental observations reported in ref .
@xcite .    in section  [ sec : vibrational - freedom ] we instead consider a different effect , arising from the interplay of vibrational modes with the matter - light coupling , which is able to reproduce similar behavior to that observed in experiments .
specifically we find that vibrational excitations dress the residual excitonic peak in a strongly temperature dependent manner .
moreover , the form of the vibrational dressed spectrum shows that the spectral feature at the exciton energy can have a more complex interpretation than that previously considered  @xcite .    a brief but self - contained account of the main theoretical methods used throughout this paper is given in the appendices .
in this section we consider whether self - consistent molecular adaptation can enhance matter - light coupling by renormalizing the bare matter - light coupling strength .
we consider models in which the effective matter - light coupling strength of a given molecule depends on the configuration of that molecule , such as its orientation , or its vibrational state .
we then ask how this same matter - light coupling modifies the energy landscape for the auxiliary parameters describing the configuration .
this leads to the idea of self consistency  if strong coupling leads to a reduction of the ground state energy , the energy landscape is deformed so as to favor auxiliary parameters for which the effective matter - light coupling is as large as possible .
our aim is to derive this from a microscopic model , and so quantify this effect . in the following we consider two potential scenarios involving adaptation of either orientational or vibrational degrees of freedom .
if such a self - consistent enhancement of matter - light coupling occurs , then this can lead to a temperature dependent effective coupling , and thus to a temperature dependence of the residual excitonic peak .
we show that such an effect exists , but that its strength is relatively weak , and that the relevant energy scale shows no collective enhancement , i.e.  the presence of @xmath5 molecules does not lead to a @xmath5 enhancement of this energy scale , because it must compete with the extensive entropy gain from orientational or vibrational disorder . as such , while increasing the molecular concentration will increase the polariton splitting , it has little effect on the self - consistent orientation . changing the bare oscillator strength of the molecules
does however affect both the polariton splitting and the self - consistent molecular adaptation energy scale .    before introducing any auxiliary variables , the basic hamiltonian which we consider is an extended dicke model , including diamagnetic terms : @xmath6,\end{gathered}\ ] ] where the field @xmath7 is written in terms of the bosonic creation and annihilation operators @xmath8 describing photon modes labeled by their in - plane momentum @xmath9 , and energy @xmath10 .
the pauli matrices @xmath11 describe the electronic state of the molecules . for completeness
, we therefore also included the diamagnetic terms , arising from the @xmath12 term of the minimal coupling hamiltonian  @xcite .
furthermore , in order to consider varying the cavity mode volume while respecting the thomas - reiche - kuhn sum rule  @xcite we use @xmath13 where @xmath14 is the oscillator strength of the given molecule .
the coefficient @xmath15 depends on the electric field strength of a single photon , and the properties of the effective charges that respond to the field .
if we assume an uniform distribution of molecules , momentum conservation can be used to write the diamagnetic term in eq .
( [ eq:1 ] ) as @xmath16 .
this term can then be removed by a bogoliubov transformation @xmath17 , yielding the effective hamiltonian : @xmath18 where the on - site hamiltonian is given by @xmath19 and the renormalized parameters @xmath20 are : @xmath21 with @xmath5 the number of molecules in the mode volume .    in the dipole approximation , for a molecule with a single mobile electron @xmath22 with @xmath23 the mode volume and @xmath24 quantifying the electronic response of a single electron in terms of the vacuum permittivity @xmath25 and
its reduced mass @xmath26 . for molecules involving many conjugated bonds ,
the coefficient @xmath27 is replaced by a sum over all mobile charges .
it is important to note that changing the cavity length changes both the mode volume @xmath23 ( and hence both @xmath15 and @xmath28 ) and the spectrum of photon modes @xmath29 . in the following numerical results we will fix the values of the polariton splitting @xmath30 , exciton energy @xmath31 , and coupling strength @xmath32 , and use these to determine @xmath33 .    before considering the role of vibrational and orientational degrees of freedom ,
we consider how the existence of a temperature dependent matter - light coupling strength , @xmath34 would be seen in the absorption spectrum .
for this purpose , it is sufficient to consider the absorption spectrum of eq .
( [ eq:2 ] ) , and its dependence on @xmath35 .
appendix  [ sec : absorpt - transm - refl ] summarizes how the absorption spectrum can be calculated , including the counter - rotating terms in the matter - light coupling  @xcite . in performing these calculations , as noted above ,
it is necessary to include disorder in order to see the residual excitonic feature .
we will consider disorder in the exciton energies , denoted by the energy distribution @xmath36 for simplicity we consider only disorder in the energies and we ignore the subleading effect of the exciton energy distribution on the coupling strength @xmath28 , by using the averaged value @xmath37 . as discussed in the appendix
, we consider the quantity @xmath38 $ ] , in terms of the retarded green s @xmath39 , which is proportional to the absorption spectrum for a good cavity .
the green s function has the form : @xmath40 where the excitonic self energy for eq .
( [ eq:2 ] ) , can be written as @xmath41 with : @xmath42 here @xmath43 is the inverse temperature and @xmath44 is the homogeneous linewidth of the excitons . in the following
we will present results both for @xmath45 and for small but non - zero @xmath44 as indicated in the figure captions .
the spectrum is shown in fig .  [ fig : abs_simple ] , focusing on the residual excitonic peak , to show its dependence on the polariton splitting . as noted above ,
such a feature has been observed and commented on several times before , e.g.  @xcite . when the polariton splitting @xmath46 is increased , the splitting between the lower and upper polaritons increases , and the exciton spectral weight of the feature at the exciton energy decreases .
the asymmetry between the shift of the lower and upper polaritons arises due to the diamagnetic term @xmath47 renormalizing the photon energy . if the coupling to light is weak , so that the rwa is valid , this asymmetry vanishes . in fig .
[ fig : abs_simple ] ( b ) we show how the spectrum is modified by including the effects of cavity losses and non - radiative excitonic decay as described in appendix  [ sec : absorpt - transm - refl ] .
the main effect these processes have is to broaden and thus reduce the height of the polariton peaks , there is also additional broadening to the central peak .    because the only temperature dependence of eq .
( [ eq:3],[eq:4 ] ) is via the combination @xmath48 , the spectrum is temperature independent while @xmath49 ( with @xmath50 of the order of the exciton bare energy ) .
therefore , as noted in the introduction , the experimentally observed temperature dependence can not occur from this mechanism alone , unless the effective value of @xmath28 is made temperature dependent via its dependence on configuration .
we now go on to consider if this effect is plausible
.     absorption spectrum , @xmath51 vs polariton splitting @xmath52 for the system with no auxiliary degrees of freedom . plotted for @xmath53ev ( i.e.  @xmath54k ) , @xmath55ev , and a truncated gaussian distribution @xmath56 with @xmath57ev , @xmath58ev .
panel ( a ) shows the results with a perfect cavity while ( b ) includes the effects of cavity losses at rate @xmath59ev and excitonic non - radiative decay at rate @xmath60ev also the width of the energy distribution is larger , @xmath61ev as discussed in appendix  [ sec : absorpt - transm - refl].,title="fig:",width=307 ]    absorption spectrum , @xmath51 vs polariton splitting @xmath52 for the system with no auxiliary degrees of freedom .
plotted for @xmath53ev ( i.e.  @xmath54k ) , @xmath55ev , and a truncated gaussian distribution @xmath56 with @xmath57ev , @xmath58ev .
panel ( a ) shows the results with a perfect cavity while ( b ) includes the effects of cavity losses at rate @xmath59ev and excitonic non - radiative decay at rate @xmath60ev also the width of the energy distribution is larger , @xmath61ev as discussed in appendix  [ sec : absorpt - transm - refl].,title="fig:",width=307 ]    to consider molecular adaptation , the hamiltonian in eq .
( [ eq:2 ] ) will be modified to include either orientational or vibrational degrees of freedom .
these are illustrated in fig .
[ fig : cartoon ] .
we next introduce these modifications and then discuss how they may be treated .
the hamiltonian including an orientational degree of freedom takes the same form as eq .   with the on - site part @xmath62 where the angle @xmath63 parametrizes the orientation of the dipole moment of the @xmath64th molecule with respect to the polarization of the cavity electric field , thus reducing the oscillator strength .
the term @xmath65 represents the bare dependence of the hamiltonian on orientation .
this term allows one to model pinning of the orientation @xmath63 . for simplicity
we consider only a classical orientational degree of freedom ; the corresponding quantum theory would require us to also include a rotational kinetic energy term , and diagonalize the resulting hamiltonian .
the effective matter - light coupling strength , which depends on the distribution of angles @xmath63 adopted by the molecule , can be written as @xmath66 , where the double angle brackets represent both an ensemble and thermal average .    to consider vibrational degrees of freedom , we again start from the transformed hamiltonian , eq .
( [ eq:2 ] ) , and now consider the following modification : @xmath67 here @xmath68 describe the @xmath69th harmonic vibrational mode of molecule @xmath64 , with the mode having frequency @xmath70 and its coupling to the electronic state being parametrized by the huang - rhys parameter @xmath71 . in this case , defining the effective oscillator strength is more involved : the effective oscillator strength depends on the matrix element describing the overlap between the vibrational states in the ground and excited state manifold .
we return to this point in later sections .    for both the orientational and vibrational degrees of freedom ,
our aim is to find how the matter - light coupling is self - consistently modified by these auxiliary degrees of freedom : i.e.  how the presence of matter - light coupling modifies the distribution of orientational or vibrational states , and how this in turn affects the effective matter - light coupling strength .
we consider a case without any strong pumping , and with a temperature such that @xmath72 , which is typically satisfied even at room temperature for organic polaritons . as such , the origin of the `` self - consistent '' dependence of configuration on the matter - light coupling coupling arises due to the existence of the counter - rotating terms in the original hamiltonian : if these terms were neglected then the energy in the ground state sector can be trivially found as the ground state would correspond to the empty state , and its energy would therefore not involve the matter - light coupling strength at all  @xcite .
the presence of counter - rotating terms means that the ground state sector also involves an admixture of all even parity sectors , and the degree of admixture depends on the effective matter - light coupling .    as discussed below , while exact solutions are possible in some limiting cases of the orientational problem , these are not generally possible at finite temperature , nor for the vibrational problem .
this is because thermal or quantum fluctuations of the auxiliary degrees of freedom break translational symmetry , preventing simple exact diagonalization . as such , we proceed using the schrieffer - wolff formalism  @xcite , which allows us to consider perturbatively the effects of these counter - rotating terms , and how they modify the energy landscape seen by the auxiliary orientational or vibrational degrees of freedom . for completeness ,
appendix  [ sec : schr - wolff - transf ] provides a brief summary of the schrieffer - wolff formalism .
the essential point is to separate @xmath73 , where @xmath74 are the terms treated perturbatively . at leading order
this gives an effective hamiltonian : @xmath75 ,    \quad    { \hat{g } } : \     [ { \hat{g } } , \hat{h}_0 ] \equiv i \hat{h}_1.\ ] ] taking the counter - rotating terms as @xmath74 , the perturbation theory is controlled by the small parameter @xmath76 , which is indeed small for the physical parameters we consider   can still be comparable to @xmath77 , as is indeed the case for the parameters we consider . ] .
we are thus in a regime where the counter - rotating terms can not be ignored , but where they can be included perturbatively . in the following sections we apply this approach in turn to the orientational and vibrational degrees of freedom , and
see how the effective matter - light coupling can be derived self - consistently .
as discussed above , we consider first the classical orientational degrees of freedom @xmath63 , subject to a pinning potential @xmath65 .
if all molecules are identical , @xmath78 , then one can find the zero temperature ground state by choosing @xmath79 .
for the ground state , all that is required is to find the quantum ground state of eq .
( [ eq:6 ] ) as a function of @xmath80 and then minimize over @xmath80 .
since eq .
( [ eq:6 ] ) is translationally invariant in the case @xmath79 , it is possible to find the exact ground state in the bosonic approximation by fourier transforming .
the bosonic approximation assumes the occupation of each excited molecule is small , a result that is valid unless @xmath81  @xcite .
however , at finite temperature , even for identical molecules , it is crucial to allow independent fluctuations of each @xmath63 ; assuming @xmath79 massively underestimates the entropy at finite temperature . at zero temperature
, one may compare the exact solution to the schrieffer - wolff perturbative expansion used below , and one finds that these indeed match to leading order .    for this rotational case , the form of @xmath82 required in eq .
( [ eq:8 ] ) can be found trivially , and the resulting hamiltonian can most conveniently be written as @xmath83 , with the molecular hamiltonian in the form : @xmath84 where @xmath85 is the bare molecular hamiltonian , including the rwa coupling to light and including the pinning term @xmath65 .    in order to consider the thermal distribution of @xmath63 , we must specify the orientational potential @xmath65 .
we consider a form @xmath86 which tries to pin the molecules at angle @xmath87 , relative to the cavity electric field , with strength @xmath88 .
we may thus consider both the free orientation case , @xmath89 , and the pinned case simultaneously . in
what follows we will assume that the pinning angles have a uniform distribution such as would be found in a polymer matrix .
this treatment , however , ignores effects which would be important in systems such as organic crystals in which the pinning angle has a fixed direction .
the energy landscape for each angle @xmath63 , given the pinning angle @xmath87 is then : @xmath90 where for simplicity we have assumed that each molecule has the same values of @xmath91 and the only disorder is in the pinning angles . in the following we will define the quantity : @xmath92 the quantity @xmath93 characterizes the self - consistent energy favoring alignment of molecules .
replacing the summation by an integral , and inserting the explicit forms of @xmath94 and @xmath15 written above we have that @xmath95 where @xmath96 is again the cavity length , @xmath27 the combination defined following eq .
( [ eq:2 ] ) , and @xmath97 is a cutoff reflecting the breakdown of the dipole approximation . to evaluate such integrals
it is useful to write the dispersion in the form @xmath98 which allows us to find the exact result @xmath99    a notable feature of eq .
( [ eq : k0exact ] ) , as anticipated above , is that this quantity _ does not simply increase _ as one increases the polariton splitting by varying the density of emitters , @xmath100 .
there is a sub - leading dependence on the molecule density , via the renormalization @xmath101 given in eq .
( [ eq:14 ] ) .
however this effect is only significant in the deep strong coupling limit  @xcite .
physically this lack of scaling with @xmath5 is because this `` molecular adaptation energy '' depends on the shift seen for _ each _ molecule . inserting typical experimental values for an organic system @xmath102 , @xmath103ev , @xmath104 nm , @xmath105 , @xmath106 nm , @xmath107ev@xmath108nm@xmath109 which correspond to the values extracted from ref .
@xcite , one finds that @xmath110 which is much smaller than @xmath111 at room temperature .
as noted earlier , the effective matter - light coupling strength is given by @xmath112 at zero temperature , this corresponds to minimizing the energy , leading to @xmath113 .
however , since the value of @xmath93 given above is such that @xmath114 , it is crucial to consider finite temperatures .
the smallness of the ratio @xmath115 will also allow us to make further perturbative expansions in the following .
defining @xmath116 , this ratio can be calculated as @xmath117,\ ] ] where the partition function is @xmath118.\ ] ] as noted above @xmath119 , and so we may taylor expand in this small parameter to get a closed form for @xmath120 . assuming that the pinning angle distribution @xmath121 is uniform we obtain the simple result for the coupling @xmath122\ ] ] with @xmath123 the imaginary bessel function @xmath124 at this point we have made no assumption about the pinning strength @xmath88 , and the value of the effective coupling is controlled by the combination @xmath125 .
if the pinning is strong @xmath126 then we find the asymptotic form @xmath127 this no longer depends on temperature as the strong pinning limit means entropy become unimportant .
thus , in the @xmath128 limit the coupling takes on the isotropic value of @xmath129 , corresponding to the uniform distribution of angles @xmath130 .
the effective coupling increases as the pinning @xmath88 decreases . in the limit of vanishing pinning @xmath131 , the imaginary bessel function @xmath132 vanishes and
so we have : @xmath133 both eq .
( [ eq:9 ] ) and eq .
( [ eq:10 ] ) indicate that as long as @xmath134 , the modification and temperature dependence of @xmath135 is very small .
we note that while in the organic systems which we focus on here @xmath93 is relatively small , in systems which have very small mode volumes this parameter could be engineered to be much larger and hence the coupling strength renormalization would be much more pronounced .
i.e. , since there is no scaling with @xmath5 , the crucial feature to see a strong renormalization is to minimize the mode volume in absolute terms , and not the mode volume _ per molecule_. this suggests evanescently confined radiation modes in plasmonic  @xcite or phonon polariton  @xcite systems may be a promising venue to explore this physics .      in the vibrational case , even without disorder , an exact solution of eq .
( [ eq:7 ] ) via fourier transformation is no longer possible , because the vibrational degrees of freedom are modeled as quantum degrees of freedom with their own quantum dynamics . as such
, they break translational invariance  i.e.  localized vibrational excitations can scatter between different polariton momentum states .
thus , once again we must use the schrieffer - wolff formalism . if we start from eq .
( [ eq:7 ] ) , solving the equation @xmath136=i\hat{h}_1 $ ] is now more challenging than for the rotational case , as @xmath137 involves terms that couple the electronic state to the vibrational quantum state . the equation can however be solved in the form of a power series , @xmath138 where the operators @xmath139 are defined by the recursion relation @xmath140       \\- \hat{o}_{j-1 }      \sqrt{s_m}({\hat{b}^{}}_m+{\hat{b}^\dagger}_m ) \biggr\},\end{gathered}\ ] ] with the base case @xmath141 .
this expansion then allows one to write out the effective hamiltonian in the same form as above , but with the molecular hamiltonian , @xmath142    \sum_{{\mathbf{k}},j } \frac{\tilde{g}_{{\mathbf{k}}}^2     } { ( \tilde{\omega}_{{\mathbf{k } } } + \epsilon_n)^{j+1 } }     ( \hat{o}_{j } + \hat{o}_{j}^{\dagger}),\end{aligned}\ ] ] with @xmath85 the bare molecular hamiltonian , including the rwa coupling to light , and vibrational terms of eq .
( [ eq:7 ] ) .    this expression can be considered as a multinomial power series in the quantities @xmath143 for each vibrational mode @xmath69 . for typical parameters ,
such quantities are small , and so we may truncate at first order , i.e.  keep terms up to @xmath144 . beyond @xmath144 ,
the expression becomes considerably more complicated , as cross terms between different vibrational modes appear .
up to @xmath144 , we find an effective molecular hamiltonian which we write out in full ( neglecting constant terms ) : @xmath145       \sqrt{s_m }      ( b_m+{\hat{b}^\dagger}_m )      \right).\end{gathered}\ ] ] the @xmath146 term gave an energy shift of two - level systems with the same form @xmath93 found previously .
the @xmath144 term is on the last line , and describes a shift to the vibrational modes .
the coefficients for these terms are defined by a generalization of that used in eq . for the rotational case above , @xmath147 where we may find an analytic form for the resulting integral @xmath148
we may once again note that none of the terms @xmath149 are proportional to the number of molecules @xmath5 : vacuum - state molecular adaptation is not collectively enhanced .    from the form of eq .
( [ eq:12 ] ) it is clear that the term in @xmath150 describes a reduction in the offset between the vibrational ground state in the two electronic states , as the virtual excitations admix the excited electronic state configuration into the ground state .
this can be viewed as a reduction of the huang - rhys parameter , @xmath151 $ ] .
this point is illustrated in fig .
[ fig : cartoon ] .
as noted earlier , the dependence of @xmath120 on the vibrational degrees of freedom is more complicated : one must calculate the overlap between the vibrational states of the ground and excited electronic state manifolds .
the reason that the vibrational states differ in these manifolds is the existence of the terms @xmath152 , which correspond to displacement of the vibrational coordinate dependent on the electronic state of the molecule . as such
, a reduction of @xmath153 would lead to an enhanced matter - light coupling .
however , using the same typical experimental values as before we find that this dimensionless shift has a value of approximately @xmath154 .
thus , once again one may conclude that the characteristic scale of any vibrational molecular adaptation ( determined by @xmath150 ) is negligibly small .
the discussion so far has focused on two specific microscopic mechanisms which might have led to self - consistent adaptation of the molecules so as to enhance their coupling to light . here
we note those aspects of the above results which can be easily generalized to other microscopic mechanisms .
examples of such other potential mechanisms include solvation of the molecule , molecular configuration , and charge transfer state . in some cases
, these will require detailed modeling of the specific process , particularly for degrees of freedom with energy scales larger than temperature , where quantum effects become important , as in the example of vibrational modes above .
however , the basic idea can be illustrated in the simplest case , where some aspect of configuration can be parametrized by a classical variable .    in the case where classical parametrization applies
, the description is very similar to the discussion of orientational degrees of freedom : one considers a classical variable @xmath155 which parametrizes some aspect the state for molecule @xmath64 .
associated with this will be an energy function @xmath156 , and , for self - consistent adaptation to occur , the matter - light coupling must depend on this variable as @xmath157 . the same perturbative analysis as discussed above will then lead to the effective energy function : @xmath158 .
this in turn means that the self - consistent effective matter - light coupling will take the form : @xmath159.\ ] ] since this involves the same sum @xmath93 as defined in eq .
( [ eq : defk0 ] ) , the same absence of scaling with number of molecules occurs i.e.  it is the single - molecular , rather than collective , coupling which is important .
moreover , since @xmath134 , then : @xmath160 where @xmath161 indicates thermal averaging with the bare energy function @xmath156 , i.e.  any such vacuum - state self - consistent adaptation of molecules is suppressed by the small parameter @xmath162 .
in the previous section we have seen that while self - consistent molecular adaptation due to matter - light coupling is possible , neither the strength nor the dependence upon molecular concentration are consistent with the results reported in ref .
@xcite . in this section
we show how directly calculating the absorption spectrum from the model of vibrationally dressed polaritons introduced in the previous section can lead to a rather different mechanism which could however be responsible for the features observed .
this mechanism arises from the combination of vibrational dressing of the spectrum and the effects of disorder .
we focus in particular on the peak in the spectrum near the bare excitonic resonance .
we will see that the shape of this feature depends strongly on the vibrational state of the molecules .
this section is divided into two subsections , in section  [ sec : self - energy - incl ] we first discuss how vibrational excitations should be included in the calculation of the disordered polariton spectrum .
this follows the method outlined in appendix  [ sec : schr - wolff - transf ] , and so all that is required is to calculate the excitonic self energy .
we then discuss the resulting form of the spectrum in section  [ sec : evol - absorpt - spectr ] .      as discussed in appendix  [
sec : absorpt - transm - refl ] , the absorption , emission and transmission spectra can all be found in terms of the photon retarded green s function . in this approach , all the properties of the molecules ( i.e.  inhomogeneous broadening , vibrational dressing etc . )
are incorporated via the excitonic self energy @xmath163 .
we must therefore calculate this quantity , defined in eq .
( [ eq : self - energy],[eq:13 ] ) , for the vibrationally dressed hamiltonian , eq .
( [ eq:7 ] ) . to do this , it is useful to label the eigenstates as @xmath164 , corresponding to the vibrational eigenstates in the excited electronic state manifold and the ground electronic state manifold .
we denote the energies of these states as @xmath165 , and @xmath166 and we introduce the overlap matrix element : @xmath167 .
this matrix element describes the extent to which the state with @xmath168 vibrational excitations in the electronic ground manifold overlaps with the state in the excited electronic manifold with @xmath169 excitations . in order to find these eigenvalues and eigenfunctions we must diagonalize the vibrational problem , @xmath170 in the electronic ground and excited states . as discussed in sec .
[ sec : vibr - degr - freed ] , the renormalization of these parameters due to virtual pair creation is very small , and so we neglect it in the following .    in terms of these quantities we may write the self energy , including inhomogeneous broadening of excitonic energies in the form : @xmath171 } {   \nu+i\gamma + ( e_q^{\downarrow}-e_p^{\uparrow } ) } ,   \label{eq : rwa - self - energy}\end{aligned}\ ] ] where @xmath36 is again the distribution of exciton energies , as in sec .
[ sec : rotational - freedom ] , and we have again ignored the subleading effects of the exciton energy distribution on the coupling strength , by using @xmath172 . it should be noted that the energies @xmath173 depend ( linearly ) on the energy @xmath50 , but that the matrix elements @xmath174 are not dependent on this energy scale .    in the following , we model the inhomogeneous broadening of excitons by using the distribution @xmath175 , i.e.  a truncated gaussian distribution of excitonic energies .
the truncation has little effect , but is formally required , as negative energy states can not physically exist , and cause problems in the ultrastrong coupling limit  @xcite .
the exciton spectrum in the absence of vibrational dressing was shown previously in fig .
[ fig : abs_simple ] . before exploring the effect of vibrational modes , it is helpful to summarize how the residual exciton peak arises .
mathematically , the excitonic feature in the polaritonic spectrum can be understood directly from the form of eq .
( [ eq:3 ] ) and the definition of the absorption spectrum .
the peak at the exciton frequency occurs because the imaginary part of the retarded green s function has a numerator involving the imaginary part of the self energy , @xmath176 $ ] , and this self energy has a peak at the excitonic energy .
however , the weight of the residual excitonic peak reduces as the matter - light coupling increases , because the excitonic self energy also appears , squared , in the denominator .
it is also important to note that this residual excitonic peak does not appear in the transmission spectrum : since @xmath177 , there is no term in the numerator of @xmath178 arising from the exciton self energy .
thus , the residual exictonic peak is a feature of the absorption ( and reflection ) , but not the transmission spectrum .     for a system with a single vibrational mode with @xmath179ev .
the insets show the equivalent absorption spectra for bare molecules .
the other parameters are as in fig .
[ fig : abs_simple ] .
panel ( a ) shows the spectrum for a perfect cavity while panel ( b ) includes the effects of cavity losses and non - radiative excitonic decay exactly as in fig .
[ fig : abs_simple](b).,title="fig:",width=307 ]   for a system with a single vibrational mode with @xmath179ev .
the insets show the equivalent absorption spectra for bare molecules .
the other parameters are as in fig .
[ fig : abs_simple ] .
panel ( a ) shows the spectrum for a perfect cavity while panel ( b ) includes the effects of cavity losses and non - radiative excitonic decay exactly as in fig .
[ fig : abs_simple](b).,title="fig:",width=307 ]    figure  [ fig : abs_vs_g ] shows an equivalent set of spectra to fig .  [ fig : abs_simple ] , but with vibrational dressing . for comparison in the insets
we show the absorption spectra of a bare molecule @xmath176 $ ] , i.e.  the spectra without a cavity . panel ( a ) shows the simple case where cavity losses are ignored , while panel ( b ) shows the more experimentally relevant case where these effects are included ( discussed in appendix  [ sec : absorpt - transm - refl ] ) . in the following we use the notation @xmath180 which denotes transitions in the molecule from the state with @xmath169 vibrational excitations in the electronic ground state to the state with @xmath168 vibrational excitations in the electronic excited state .
the presence of the vibrational modes causes dramatic changes to the residual excitonic peak in the polariton spectrum not seen in the bare excitonic absorption . for the bare molecule the spectral weight associated with transition with the `` zero phonon line '' , i.e.  the transition denoted @xmath181 in the notation introduced above
is completely dominant , only a small amount of weight visible in the sideband which corresponds to @xmath182 transition .
when the molecule is placed inside a cavity , the vibrational sidebands corresponding to the @xmath183 and @xmath182 transitions become much more prominent .
physically this is because the spectral weight that had been associated with the @xmath181 transition in the bare molecular spectrum has been moved into the polariton peaks of the spectrum .
the mathematical form of the green s function makes clear that formation of the polariton spectral feature is predominantly at the expense of whatever feature dominates the excitonic emission spectrum ; here this is the @xmath181 feature .
the excitonic feature corresponds to the `` left - over '' spectral weight associated with the subradiant states .
thus , the vibrational sidebands have been `` excavated '' by removing the dominant feature from the zero - phonon line . including the the effects of cavity losses and non - radiative excitonic decays , as can be seen in fig .
[ fig : abs_vs_g ]  ( b ) , washes out this complex sideband structure of the central peaks and the spectrum as a function of coupling strength looks very similar to that obtained without coupling to vibrational modes , as in fig .
[ fig : abs_simple ] .
however as we discuss below , there are still effects which can be observed which are a direct consequence of the vibrational structure .    , ( b ) coupling to a single vibrational mode with @xmath184ev .
the inset shows the bare molecular spectrum .
( c ) including the effects of cavity losses and non - radiative excitonic decay .
all panels are for @xmath185ev , and other parameters as in previous figures.,title="fig:",width=307 ] , ( b ) coupling to a single vibrational mode with @xmath184ev .
the inset shows the bare molecular spectrum .
( c ) including the effects of cavity losses and non - radiative excitonic decay .
all panels are for @xmath185ev , and other parameters as in previous figures.,title="fig:",width=307 ] , ( b ) coupling to a single vibrational mode with @xmath184ev .
the inset shows the bare molecular spectrum .
( c ) including the effects of cavity losses and non - radiative excitonic decay .
all panels are for @xmath185ev , and other parameters as in previous figures.,title="fig:",width=307 ]    fig .
[ fig : abs_vs_t ] illustrates the evolution of spectra with temperature .
panel ( a ) shows that without coupling to vibrational modes there is no notable temperature dependence . in the presence of the vibrational dressing
, a strong temperature dependence appears . at higher temperatures ,
there is a greater thermal occupation of the vibrational modes hence the spectral weight under the vibrational peaks rises . while this has a small effect in the bare molecular spectrum ( where the @xmath181 transition dwarfs all other features ) ,
see inset in ( b ) , it is very pronounced in the polariton spectrum and is even visible in the presence of large cavity losses as in fig .
[ fig : abs_vs_t](c ) .    .
the main panel shows the optical spectrum of the strongly coupled system , the inset shows the spectrum of the bare excitons . in this case , we include two vibrational modes : @xmath186ev , and @xmath1874ev .
other parameters as in fig .
[ fig : abs_vs_t].,width=307 ]    the figures so far have shown results where disorder is relatively small , and so vibronic replicas can be clearly observed for the good cavity , but merge for the bad cavity limit . to show that small disorder is not required for the strong temperature dependence to occur fig .  [
fig : abs_large_disorder ] shows the effects of large disorder ( i.e. inhomogeneous broadening ) . just as seen for the homogeneous broadening in fig .
[ fig : abs_vs_t](c ) , a temperature dependence of the residual excitonic peak is still visible . in this figure
, we have also included a more complicated vibrational spectrum , involving two vibrational modes .
such a system will exhibit behavior similar to that of a single mode but with a large vibrational coupling  , then the effective huang - rhys parameter will be @xmath188 .
in this paper we have presented two microscopic models which could in principle describe self - consistent molecular adaptation so as to maximize the vacuum - state coupling to light . in both cases ,
the crucial feature of the model is the counter - rotating terms in the matter - light coupling .
these allow virtual fluctuations in the ground state , that lower the ground state energy depending on the configuration of the molecules .
this energy gain is the only energy gain that can be relevant in the linear response regime  i.e.  the question of whether the excited states would have lower energy is not of relevance while the system is only weakly pumped .
we found that while such a mechanism for molecular adaptation does exist , it does not show any collective enhancement , in contrast to the polariton splitting , and does not therefore lead to significant molecular adaptation , even when the polariton splitting @xmath189 .
the appearance of a residual exciton peak in the polariton spectrum would be affected by any such self - consistent molecular adaptation if its scale were sufficient . however , for relevant parameters , such effects are dwarfed by a far more dramatic effect , of vibrational dressing of the residual exciton peak .
this leads to a pronounced temperature dependence of the feature near the exciton energy in the absorption spectrum .
while the molecular adaptation energy scale is not collectively enhanced , the basic underlying physics could potentially be relevant in single molecule strong coupling , e.g.  with plasmonic resonances  @xcite .
in such cases , rather than having many molecules in a large mode volume , the mode volume is reduced so that strong , or even ultra - strong coupling occurs at the single molecule level , so that both the polariton energy and the molecular adaptation energy become large .
in such a case , one may hope to see either reorientation , or renormalization of the huang - rhys parameter due to strong coupling .
another intriguing direction for future research is to consider how the physics discussed in this manuscript interacts with the physics of polariton condensation and lasing  @xcite .
polariton condensation has been seen in both inorganic  @xcite and organic  @xcite systems .
in addition , condensation of photons has been seen for weakly coupled systems of organic molecules  @xcite .
theoretical work  @xcite has begun to address some of the peculiarities of the organic polariton system , including effects of disorder and of vibrational modes .
however , features as seen in this paper , resulting from the interplay of these may lead to further exotic behavior in the high density condensed phase .
we are grateful for comments from t.  ebbesen on an earlier version of this paper .
jk acknowledges helpful discussions with david lidzey and brendon lovett .
jac acknowledges support from epsrc .
jk and pgk acknowledge financial support from epsrc program `` topnes '' ( ep / i031014/1 ) .
jk acknowledges support from the leverhulme trust ( iaf-2014 - 025 ) .
sdl is royal society research fellow .
sdl acknowledges financial support from epsrc grant ep / m003183/1 .
pgk acknowledges support from epsrc grant ep / m010910/1 .
jk , pgk and sdl acknowledge support from the british council for the meeting which initiated this work .    [
[ note - added ] ] note added : + + + + + + + + + + +    during the final preparation of this manuscript , another paper  @xcite appeared , also reporting the fact that ground - state bond length depends on the single - molecule coupling @xmath190 , not the collective coupling @xmath46 .
in this appendix we summarize the calculation of the absorption , transmission and reflection spectra .
some subtleties arise because we wish to calculate the spectrum of a model with ultrastrong coupling , i.e.  without making the rotating wave approximation .
such results were first calculated by  @xcite , here we present a synopsis of these results , as well as a `` dictionary '' to translate the results of that paper into the language of green s functions .
we begin by defining the retarded green s function  @xcite .
because we consider both co- and counter - rotating terms , we must consider both normal and anomalous green s functions .
i.e.  we must include number non - conserving terms which appear for ultrastrong coupling , and thus we consider a matrix green s function : @xmath191 \rangle      &      \langle [ \psi^{\dagger}_{-{\mathbf{k}}}(t ) , \psi^\dagger_{{\mathbf{k}}}(t^\prime ) ] \rangle      \\
\langle [ \psi^{}_{{\mathbf{k}}}(t ) , \psi^{}_{-{\mathbf{k}}}(t^\prime ) ] \rangle      &      \langle [ \psi^{\dagger}_{-{\mathbf{k}}}(t ) , \psi^{}_{-{\mathbf{k}}}(t^\prime ) ] \rangle    \end{pmatrix}\ ] ]    in terms of the bogoliubov transformed operators , i.e.  the operators appearing in eq .
( [ eq:2 ] ) the inverse green s function takes the form @xmath192^{-1 } =     \begin{pmatrix }      \nu + i \tilde{\kappa}(\nu ) - \tilde{\omega}_{{\mathbf{k } } } + \tilde{\sigma}^{}_{{\mathbf{k}},xx}(\nu )       &      + i \tilde{\kappa}(\nu ) + \tilde{\sigma}_{{\mathbf{k}},xx}(\nu )       \\      - i \tilde{\kappa}^\ast(-\nu ) + \tilde{\sigma}^\ast_{{\mathbf{k}},xx}(-\nu )       &      -\nu - i \tilde{\kappa}^\ast(-\nu ) - \tilde{\omega}_{{\mathbf{k } } } + \tilde{\sigma}^\ast_{{\mathbf{k}},xx}(-\nu )     \end{pmatrix},\ ] ] where @xmath193 is the self energy for a photon of in - plane momentum @xmath9 , arising from the excitonic response ( discussed further below ) , and @xmath194 is the loss rate .    here , following  @xcite we have used a frequency dependent complex loss rate @xmath194 .
frequency dependence is required for physical consistency in the case of ultrastrong coupling  markovian loss and ultrastrong coupling would predict a perpetual light source . frequency dependent loss requires , via the kramers - kronig relation , a corresponding lamb shift , which is incorporated into the imaginary part of @xmath194 .
both @xmath195 and @xmath196 are written for the bogoliubov transformed operators , and so both these terms incorporate a pre - factor @xmath197 to account for the bogoliubov transformation of the combination @xmath198 .
the specific self energy required , @xmath193 , corresponds to correlation functions of the @xmath199 excitonic operators . in the absence of strong ( i.e.  beyond rwa ) excitonic damping ( see  @xcite for the more general case ) , this self energy can however be related to results in the rotating wave approximation by : @xmath200^{\star } \right ) ,    \label{eq : self - energy}\ ] ] where the expression @xmath201 is the `` standard '' self energy that would appear in the rotating wave approximation , depending on the correlation of @xmath202 operators .
these can most easily be found by analytic continuation from imaginary time to real time , starting from the matsubara self energy , @xmath203 and replacing the matsubara frequency by @xmath204 .
the sum over @xmath169 appearing here is over all states of the excitonic system , and @xmath205 is the partition function . for the `` vacuum '' state we consider  i.e.  in the absence of strong pumping
 these are the bare exciton states , including the quantum states of any auxiliary degrees of freedom .    using the input - output formalism  @xcite adapted to the ultrastrong coupling regime  @xcite
, one can write a frequency dependent scattering matrix relating input and output fields at the left and right sides of the cavity , @xmath206 where @xmath207 with @xmath208 are the real parts of the loss rates arising from the left and right mirrors and the quantity @xmath39 relates to the matrix retarded green s function as @xmath209 where @xmath210 .
this structure means the bogoliubov transformation corresponds to @xmath211 and the bogoliubov transformed green s function takes the form : @xmath212 one can then find the transmission @xmath178 , reflection @xmath213 and absorption @xmath214 coefficients by considering the modulus square of various coefficients .
clearly @xmath215 is independent of which direction light is incident from , while the absorption coefficient @xmath216 takes the form @xmath217       + \kappa(\nu ) \left\vert g^r_{{\mathbf{k}},xx}(\nu)\right\vert ^2 \right],\ ] ] with @xmath218 .
the prefactor in this expression shows the obvious dependence on the transmissivity of the input mirrors .    in order to separate mirror transmissivity dependent features from the `` intrinsic '' properties of the ultrastrong coupling we will consider below the two quantities @xmath219 as being proportional to the transmission , and @xmath38 $ ] as controlling the absorption in the limit of a good cavity , i.e.  @xmath220 .
the quantity @xmath221 differs from the full absorption as it neglects interference effects at the input mirror .    ) , plotted for @xmath222 , @xmath223 , @xmath224ev , @xmath225ev , @xmath226ev and various cavity loss rates , @xmath227 .
other parameters are as in fig .
[ fig : abs_large_disorder ] . as discussed in the text , in order that the polaritonic peaks have a finite width , one must include the effects of excitonic absorption or non - mirror cavity losses .
for this figure , we include an excitonic linewidth @xmath60ev.,width=307 ]    for comparison to the results in ref .
@xcite where silver mirrors were used , we present in figures  [ fig : abs_simple ] , [ fig : abs_vs_g ] , and [ fig : abs_vs_t ] the effects of large cavity linewidth . figure  [ fig : true - absorption ] also shows how the full absorption spectrum given by eq .
( [ eq:15 ] ) evolves with varying linewidth . in calculating these spectra for large linewidth
, an issue arises regarding the form of the absorption spectrum : the equations written above assume that the only photon loss is due to escape through the mirrors .
this means that the @xmath228 appearing explicitly in eq .
( [ eq:15 ] ) ( describing interference effects from the mirror ) is the same as the @xmath228 appearing in the denominator of the photon green s function , eq .
( [ eq:11 ] ) , and one may check that for any frequency where @xmath229 $ ] is small , this causes a near cancellation between the two contributions . for the gaussian exciton density of states @xmath36 used in this paper , this cancellation almost completely suppresses the polariton peaks .
such a cancellation in the absorption spectrum can be expected on physical grounds : if the only loss channel for photons is the mirrors , then there is no absorption .
all photons that enter eventually leave . in a real device
there are other photon loss sources ( absorbers , scattering by surface roughness ) .
similarly , in a real device , the excitons have a non - zero rate of non - radiative decay .
the effect of this is included by retaining a non - zero value of @xmath44 in the denominator of the self - energy , eq .  .
this leads to lorentzian tails of @xmath229 $ ] , giving a finite weight to the polariton peak in the absorption spectrum ; such an effect was included by  @xcite .
we follow this approach in plotting fig .
[ fig : true - absorption ] .
in section  [ sec : rotational - freedom ] we make use of the schrieffer - wolff approximation ; for completeness we provide here a brief explanation of this formalism .
the approach is based on dividing the hamiltonian into two parts , @xmath230 , where the term @xmath231 takes one between different `` sectors '' . in our case
, these sectors correspond to different numbers of polaritons  i.e.  @xmath232 is the `` counter - rotating '' part of the hamiltonian which simultaneous creates a photon and excites a molecule .    the aim of the schrieffer - wolff formulation is to make a unitary transformation @xmath233 such that the transformed hamiltonian no longer has any coupling between sectors .
physically , this corresponds to eliminating the effect of virtual pair creation and destruction , and deriving how such virtual processes renormalize the hamiltonian within a given sector .
if the hamiltonian @xmath74 can be treated perturbatively by replacing @xmath234 with @xmath235 a small parameter , then one can consider a series solution @xmath236 . in order to make the first - order terms in @xmath235 vanish
, one must choose @xmath237=i\hat{h}_1 $ ] .
this then leads ( setting @xmath238 ) to the expression : @xmath239 +   \text{h.o.t}\end{aligned}\ ] ] where the higher order terms involve all @xmath240 .
stopping at leading order gives the expression in eq .
( [ eq:8 ] ) , corresponding to the leading order effects of virtual pair creation and annihilation . to solve @xmath241=i\hat{h}_1 $ ] in practice is straightforward if one knows the eigenspectrum of @xmath242 , which then allows one to write @xmath243    58ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty http://prola.aps.org/abstract/pr/v112/i5/p1555_1 [ * * ,   ( ) ] http://adsabs.harvard.edu/abs/1958jetp....6..785p [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physreva.44.669 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.69.3314 [ * * ,   ( ) ] @noop _ _  ( ,  , ) @noop _ _  ( ,  ,  ) @noop _ _  ( , ,  ) link:\doibase    10.1038/25692 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.82.3316 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.106.196405 [ * * , ( ) ] link:\doibase    10.1103/physrevlett.95.036401 [ * * ,   ( ) ] link:\doibase 10.1103/physrev.93.99 [ * * , ( ) ] link:\doibase 10.1007/978 - 3 - 662 - 45082 - 6 [ _ _ ] ,  edited by  ,  nano - optics and nanophotonics  ( ,  ,  ) link:\doibase 10.1103/physrevlett.101.116401 [ * * ,   ( ) ] link:\doibase 10.1038/nphoton.2010.86 [ * * ,   ( ) ] link:\doibase    10.1038/nmat3825 [ * * ,   ( ) ] link:\doibase 10.1038/nmat3874 [ * * ,   ( ) ] link:\doibase    10.1002/ange.201107033 [ * * ,   ( ) ] link:\doibase    10.1038/nmat4392 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.114.196403 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.114.196402 [ * * ,   ( ) ] link:\doibase    10.1038/ncomms6981 [ * * ,   ( ) ] link:\doibase    10.1038/nnano.2015.264 [ * * ,   ( ) ] link:\doibase 10.1088/1367 - 2630/17/5/053040 [ * * ,   ( ) ] @noop `` , ''   ( ) ,   link:\doibase 10.1063/1.4919348 [ * * ,   ( ) ] link:\doibase 10.1002/anie.201301861 [ * * ,   ( ) ] link:\doibase    10.1002/cphc.201200734 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.53.2711 [ * * , ( ) ] link:\doibase 10.1103/physrevb.64.235101 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.72.115320 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.72.115303 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.92.125433 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.79.201303 [ * * ,   ( ) ] link:\doibase 10.1021/ph500266d [ * * ,   ( ) ] @noop * * ( ) link:\doibase    10.1103/physrevb.90.205309 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.35.432 [ * * ,   ( ) ] link:\doibase 10.1103/physreva.74.033811 [ * * ,   ( ) ] link:\doibase 10.1103/physrev.149.491 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.105.263603 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.112.016401 [ * * ,   ( ) ] link:\doibase 10.1021/acs.nanolett.5b01204 [ * * ,   ( ) ] link:\doibase 10.1021/nl401590 g [ * * ,   ( ) ] link:\doibase 10.1103/physreva.80.053810 [ * * ,   ( ) ] http://stacks.iop.org/0034-4885/78/i=1/a=013901 [ * * ,   ( ) ] link:\doibase 10.1103/revmodphys.85.299 [ * * , ( ) ] link:\doibase 10.1038/nature05131 [ * * ,   ( ) ] link:\doibase    10.1126/science.1140990 [ * * ,   ( ) ] link:\doibase 10.1038/nature09567 [ * * ,   ( ) ] link:\doibase 10.1002/pssa.200304067 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.74.165320 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.77.155325 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.79.035325 [ * * , ( ) ] link:\doibase 10.1103/physrevb.80.235313 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.80.235314 [ * * , ( ) ] link:\doibase 10.1103/physrevb.88.075321 [ * * ,   ( ) ] link:\doibase 10.1209/0295 - 5075/105/47009 [ * * ,   ( ) ] link:\doibase 10.1103/physrevx.5.041022 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) link:\doibase 10.1103/physreva.30.1386 [ * * ,   ( ) ]