arXiv:0907.1836v3 [cond-mat.supr-con] 6 Sep 2010 n neatos n edt te ye funconven- of types phases. superconducting other and to ordered charge lead pair- tional and electronic phonons the interactions, and affect ing excitons significantly of still may influence however combined to phase The superconducting exciton-driven impos- form. truly it a makes for phonons sible and excitons between coupling ua ntblt oocrbfr uecnutvt a a has superconductivity arise before to struc- chance a occur re- cause to to typically tends instability excitons which tural the renor- to spectrum coupling the from phonon material sults real the any of in because malisation turned elusive, be has to phase out superconducting exciton-mediated an t omlyipsdb hnn,tu aigtewyfor superconductivity way the temperature paving thus room re- phonons, lim- by soon imposed temperature was normally all avoids its It ‘glue’ excitonic construction. the BCS that alised usual the the replacing in electrons, phonons between interaction pairing a ate hs mriguo t upeso.A a suggested was ago As decades four suppression. superconducting already its the upon to emerging interest phase cou- extra electron-phonon adding as thus well pling, as formation exciton volves ree rniinmtldichalcogenides metal transition ordered n fcag re nmtrassc stenal fer- nearly the as such materials SrTiO in roelectric order charge materials of fermion heavy bor- and the in on antiferromagnetism found of been ders has Non-BCS type used this interactions. of be pairing superconductivity may novel and mediate point, directly critical to val- quantum parameter the tuning and around physics temperature ues of the range influence wide a temperature in or- zero the at critical of parameter The disappearance der the transition. with phase associated quantum fluctuations a elec- to the order tuning upon tronic arise often super- which unconventional states, of conducting realisation the for ground fertile xio-hnnitrcin n uecnutvt borde superconductivity and interactions Exciton-phonon pcfi xml famtra nwihteexcitonic the which in material a of example specific A h D hs nsm ftelte aeil in- materials latter the of some in phase CDW The togycreae lcrnsseshv rvnt be to proven have systems electron correlated Strongly h hreodrdsae uhasprodcigpaehsr has phase superconducting a Such exc TiSe same state. the in by ordered mediated charge materials, the same the in conductivity est ae nmtrassc stetasto-ea dic transition-metal the as such materials in waves density ASnmes 74.40.Kb,71.35.Gg,63.20.Kd,71.45.Lr numbers: report PACS with tr finding the our th on compare study order and charge also region, We suppressed superconducting sugg the phase. is from superconducting emerging it single tuations results, a of these part on form Based regions. superconducting h hredniywv ntblt nTiSe in instability exciton-phonon wave combined density the elect charge that and the find formation We exciton description. of model interplay the investigate we 7 h togculn ewe atc oe n hre which charges and modes lattice between coupling strong The aeds aoaoy nvriyo abig,Madingley Cambridge, of University Laboratory, Cavendish h tog n emnl unavoidable, seemingly and strong, The . .INTRODUCTION I. 3 n h hredniywv (CDW) charge-density-wave the and 2 ohuo neclto ihCpe n ntepitn mat pristine the in and Copper with intercalation upon both , aprvnWzl alNhiWlimo,adSdat .Saxena S. Siddarth and Nahai-Williamson, Paul Wezel, van Jasper 5,6 xioscndrcl medi- directly can excitons , 7 npats such practise In . 3,4 . 2 1,2 r lorsosbefrtepiigo lcrn nits in electrons of pairing the for responsible also are oei h omto fteCWi noe TiSe undoped in CDW the of formation the in role cnl otefraino hreodri ie ythe by given is TiSe order of charge phase of signif- CDW formation contribute the both to freedom icantly of degrees phononic and fasfeigo the of softening experiments a scattering of x-ray quasi-elastic in observation pressure ambient xiosas r nidsesbeigein ntede- the in ingredient system this indispensable of an scription are also excitons h D distortion CDW the nteitraae opudCu compound superconductivity intercalated of the discoveries in recent phase the CDW backfolded context, the to of this weight measurements spectral ARPES of in transfer bands large the by ple al. eae al. et seva hrsprse rehne.I TiSe In enhanced. ei- or be exci- example, suppressed could For of ther temperature involvement mechanism. transition the pairing superconducting of electronic the be the will in effect tons the what ever n hthsbe ugse ob epnil o sta- for responsible the be to to order charge by suggested similar the mediated been phonons bilising be has and to that excitons likely one of is order combination CDW a the of pression a on.Tedtiso h prahtwrscriticality different towards very approach the are of however, criti- details quantum The potential superconducting point. the a cal covers cases which by both emerges or in dome Copper and with pressure, intercalation applying by either suppressed be h alcecetcagssg sCpe sadd sug- added, TiSe is the Copper that as gesting sign changes coefficient Hall the h upesdCWorder from CDW emerging in phases suppressed difference superconducting the a the at of hint nature to the suggested been previously have character two-band the while as- increases is constant, density remain electron hand, to the and other sumed occurs, the change on sign pressure such under no compound pure the In h uecnutn hs hc mre pnsup- upon emerges which phase superconducting The 24,25 acgndsmyas iers osuper- to rise give also may halcogenides o-hnnculn ihnasimplified a within coupling ron-phonon n npr TiSe pure in and , se htbt fteosre domes observed the of both that ested 27,28 deprmna results. experimental ed npr rprisdrcl bv the above directly properties ansport tno hnnmdsta dominate that modes phonon or iton 26 ffc fteqatmciia fluc- critical quantum the of effect e oe rvosysgetdt drive to suggested previously modes od abig B H,UK 0HE, CB3 Cambridge Road, cnl enosre o example for observed been ecently r fgetinterest. great of are , hs ieecsi aeilproperties material in differences These . ed otefraino charge of formation the to leads 8–17 ra ne rsue Here pressure. under erial igcag re nTiSe in order charge ring 18 scnb enfreapeb the by example for seen be can as , 2 2 L h atc ly nimportant an plays lattice The . ttesm iei scerthat clear is it time same the At . aesbcm oe yelectrons. by doped become layers 16,17,19–23 1 − 17 2 hnnmd soitdwith associated mode phonon ti not is It . 26 ne rsueb Kusmart- by pressure under 26,29,30 ntedpdcompound doped the In . swtesdfrexam- for witnessed as , x n as o such for cause One . TiSe priori a 2 2 h D can CDW the , yMrsnet Morosan by la how- clear 23 2 In . 2 at 2

FIG. 1: Feynman diagrams for some of the instabilities inher- ent in the Hamiltonian of eq. (1). (a) Cooper pair formation via phonon exchange. (b) Cooper pair formation arising from the exciton binding energy. (c) CDW formation promoted by electron-hole interaction via renormalised phonons. (d) FIG. 2: left A schematic representation of the band structure Cooper pair formation via renormalised phonon exchange. very close to the Fermi energy, assuming a small indirect gap between the valence and conduction bands. right Constant energy maps near EF forming electron and hole pockets in the first Brillouin zone of TiSe2 (adapted from [40]). discordance could be a change of the role played by the excitons in these compounds. We recently introduced a generic model system to study the effect of the interplay between exciton for- mation and electron-phonon coupling on the formation of CDW order16,17,31. Within this model it was shown that excitons can compete, coexist or even cooperate with phonons to drive the transition. Here we extent the pre- FIG. 3: The quasi one-dimensional double chain used as a vious model to also include superconducting pairing in- model to study the interplay of electron-phonon coupling and teractions (section II). We study the model in two differ- exciton formation. Hoppings between sites on the upper and ent regimes, and describe the effects of electron doping lower chains are described by combinations of the electronic (section III A) as well as introducing an effective pres- c and d operators, while lattice deformations are included sure (section III B). We find that although the result- through the a and b phonon operators. ing superconducting phases are BCS-like, the presence of excitons does play an essential role in renormalising the phonon properties. The cooperation of exciton for- (1), while a purely exciton-driven superconducting insta- mation and electron-phonon coupling may even lead to bility is generated by the second term (see figures 1a and a slight enhancement of the superconducting transition b). The combined effect of both of these excitations, temperature. Finally, in section IV, we consider the ef- which is thought to underlie the CDW formation in f.e. 17 fect that critical fluctuations emanating from the critical TiSe2 , may also lead to superconductivity (figures 1c point covered by the superconducting dome may have on and d). the temperature dependence of the resistivity in the nor- In the following, we will study the possible roles played mal phase. The critical exponents induced by such fluc- by excitons and phonons in the emergence of a supercon- tuations form an important experimental tool in studying 32–35 ducting phase within a specific realisation of the generic the properties of the hidden quantum critical point , structure of equation (1) which has previously been sug- and the nature of the modes driving the instability2,36–39. 16,17,31 gested to describe the CDW formation in TiSe2 . The model will be general enough to give a qualitative understanding of the influence of excitons and phonons II. A MODEL SYSTEM on the superconductivity observed in materials such as TiSe2, but it will not yield any quantitative results. The Hamiltonian governing the interplay between exci- To be specific, we only consider the two bands closest tons, phonons and conduction electrons takes the general to the Fermi energy which are known to be the Se-4p va- form: lence band located at the centre of the first Brillouin zone, and the Ti-3d conduction band forming electron pockets Hˆ =Hˆ0 + Hˆexc + Hˆe-p. (1) at the zone boundaries, as in figure 2. The Hamiltonian will be defined on the double chain lattice schematically This includes the bands of non-interacting electrons H0, shown in figure 3. The c sites represent the Se-4p orbitals, the effective Coulomb interaction between electrons and while the d† operators create electrons in the Ti-3d or- holes responsible for creating excitons Hexc, and the bitals. electron-phonon interaction He-p. A CDW instability This greatly simplified quasi one-dimensional model may arise either directly from the electron-phonon cou- has previously been argued to qualitatively capture the 16,17,31 pling (such as in a Jahn-Teller driven scenario), directly physics of the CDW transition in TiSe2 . The two- from the exciton binding energy (forming an excitonic dimensional layers making up the hexagonal TiSe2 struc- insulator), or from a combination of the two effects (for ture are separated from each other by a Van der Waal’s example via the renormalisation of the phonon spectrum gap, which allows us to neglect the effects of inter-layer by excitons). Likewise, the standard BCS type of super- coupling. Moreover, from a tight-binding fit to the band conductivity arises solely out of the last term in equation structure it becomes clear that the electronic structure 3 is dominated by a small set of orbital overlap integrals, To also include the superconducting pairing interac- which indicate that to a first approximation, the sys- tions between electrons, we follow the usual BCS pre- tem can be thought of as a network of weakly coupled scription by expanding Hˆe-p to second order in the quasi one-dimensional Ti chains and their surrounding electron-phonon coupling. The expansion terms can be Se environment31. Multiple copies of the double chain decoupled by introducing superconducting Cooper pair model of figure 3 may be put together to form a larger, densities of s-wave symmetry both within the unit cell three-dimensional structure, but because the coupling (σin = cˆi cˆi ) and between unit cells (σout = cˆi cˆi±1 ). between different double chains is weak it can be taken Noticeh thati in principle also interband pairingh termsi 17 to be a higher order effect, and we will neglect it here . of the form cˆdˆ could arise both from the electron The non-interacting part of the model can be con- phonon coupling,h i and directly from the excitonic binding structed from the tight binding fit to the band structure energy7,41–43. In the present model these pairing chan- and experimentally accessible parameters. Together they nels turn out to be repulsive, and we will neglect them form the bare potential and kinetic energy terms: from here on. The remaining pairing terms are then:

† † ∆ † † 2 ˆ ˆ ˆ ˆ ˆ ′ α † † c c H0 =¯hω aˆi aˆi + bi bi + (di di cˆi cˆi ) ˆ ˆ ˆ 2 − HMF = 2 2 dkd−k [(2 cos(ka))σin σout] i   i − mω a − − X X Xk n t ˆ† ˆ † † † d d + (d d +1 +ˆc cˆ +1 + H.c.) +ˆc cˆ (2 cos(ka))σ σ + H.c. . 2 i i i i k −k − in − out Xi   (6) ′ ˆ† ˆ† +t (di cˆi + di+1cˆi + H.c.). (2) Xi Within this mean field formulation we can investigate the effect of electron doping by adjusting the value of the The exciton binding energy is approximated as a chemical potential to yield a higher total electron density. purely local term which favours electrons and holes to In doing so, we assume that the shape of the band struc- be no further than one lattice spacing apart, while the ture around the Fermi energy will not be significantly af- electron-phonon term, encoding for Jahn-Teller effects, fected by the presence of interlayer dopants26. The effect † couples the c and d orbitals to the phonon operatorsa ˆ of pressure can be included in the model by scaling the and ˆb†, and their associated dimensionless distortions, orbital overlap integrals, in order to replicate the exper- Xˆ a = (ˆa† +ˆa)/√2mωa2 and Xˆ b = (ˆb† + ˆb)/√2mωa2: imentally observed behaviour of an increased two-band character at constant total electron density27,28. † † Hˆexc = V dˆ dˆ 1 cˆ cˆ − i i − i i Xi h   ˆ† ˆ † III. TOWARDS QUANTUM CRITICALITY + di di 1 cˆi−1cˆi−1 , (3)  − i Solving the mean field equations (5) self consistently for varying values of the exciton binding energy V and ˆ ˆ a ˆ b ˆ† † ˆ He-p = α Xi + Xi di cˆi +ˆci di the strength of the electron-hole coupling α, we find three Xi h   distinct phases at zero temperature16,17. At low V and ˆ a ˆ b ˆ† † ˆ α, the system is in its semiconducting ground state (see Xi+1 + Xi di+1cˆi +ˆci di+1 . (4) −   i figure 4). Upon increasing the exciton binding energy a first order transition is encountered at which charges These interaction terms can be decoupled on the mean from the occupied valence band transfer to the unoccu- field level by introducing parameters for the average lat- ˆ a ˆ b pied conduction band and form excitonic bound states tice distortions (u = Xi + Xi ), the electron densities h † i ˆ† ˆ with the positively charged holes left behind. This state on the different sites (ρc = cˆi cˆi and ρd = di di ), and is an (excitonic) insulator, but because we did not include the charge transfer betweenh Ti andi Se sitesh bothi within † † any direct exciton-phonon coupling in the Hamiltonian of and between unit cells (τin = cˆ dˆ and τout = cˆ dˆ ). h i i i h i i+1i equation (5), it does not break the translational symme- In terms of these, the particle number conserving part of try of the lattice. the Hamiltonian reads: Keeping a low exciton binding energy and increasing ∆ the electron-phonon coupling instead leads to a second Hˆ = cˆ† cˆ t cos(ka) +2Vρ µ order transition into a CDW ordered phase. Increasing MF  k k  − 2 d −  Xk the exciton formation at any given value of the electron- ∆ phonon coupling enhances the tendency towards charge + dˆ† dˆ t cos(ka)+ +2V (ρ 1) µ k k  2 c − −  order, and indeed there is a region in phase space where a finite exciton binding energy is required in order for a + dˆ† cˆ t′(1 + eika)+ αu(1 eika) k k − CDW to form. As V is increased even further, an ex-  ika citonic insulator will once again be formed and transla- V (τin + e τout) + H.c. (5) − tional invariance restored. The CDW observed in pristine  4

FIG. 5: The transition temperature in the mean field solution of equation (5) as a function of the exciton binding energy V/t, for a fixed value of the electron-hole coupling α/t = 0.04. The FIG. 4: The mean field solution of equation (5) for the charge − transition temperatures of both the superconducting and the transfer ρc ρd. Starting from the undistorted semiconducting charge ordered phase are enhanced by an increase in binding phase at the origin, increasing the electron-phonon coupling energy due to its renormalisation of the phonon spectrum. At α/t past a critical value induces a CDW transition, while very high values of V the charge transfer is exhausted, and an for any given value of α an increase in the exciton binding excitonic insulator is formed instead. Notice that the region energy V/t is seen to enhance the tendency towards CDW of overlap between CDW and superconducting order does not formation. At even higher values of V , the exciton formation necessarily indicate coexistence. saturates the charge transfer and a uniform excitonic insulator is formed. A. Electron Doping

The CDW observed in pristine TiSe2 under ambient pressure can be suppressed by intercalating with Copper. As the CDW transition temperature is tuned to zero, a 24,25 TiSe2 is believed to be of a hybrid exciton-phonon char- superconducting phase has been observed to emerge . acter, as in the CDW phase of figure 4. Thus excitons To describe this behaviour in the present model of equa- as well as phonons play an important role in establishing tions (5) and (6), we assume that the main consequence charge order, and both are expected to play an impor- of intercalation will be the donation of electrons from the tant role in the superconducting phase that appears as Cu intercalants to the Ti conduction band. The addition the CDW is suppressed. of Copper is known to also slightly alter the lattice pa- rameters of the TiSe2 layers, acting as an effectively neg- ative pressure24. We will assume this to be the smaller If we also include the Cooper pair densities of equation effect, and focus instead on the electron doping, which (6) in the mean field analysis, we find that a supercon- can be modelled as an effective increase in the average ducting phase may be realised for high enough values total electron density, stabilised by an enhancement of of the electron-phonon coupling. The pairing interac- the electronic chemical potential µ. tion in this phase is mediated by phonons, which in turn Upon introducing more electrons into the system de- are renormalised due to the presence of excitons. If the scribed by equation (5), the CDW transition is pushed excitonic binding energy is increased, the renormalisa- towards higher values of the electron-phonon coupling, as tion increases, and the superconducting order parameter, seen in figure 6. This initial suppression of the charge or- as well as the superconducting TC, may be slightly en- der can be easily understood in terms of the local physics, hanced, as seen figure 5. At even higher values of V there by realising that the extra electrons on the c sites block is again an instability towards the formation of an exci- the transfer of charge along the c-d bond preferred by the tonic insulator, and both superconductivity and charge lattice distortions. The transition into the excitonic insu- order are lost. Notice that in the present mean field treat- lator phase at low values of α is not affected in the same ment it is not possible to determine whether the CDW way, because in the absence of a direct exciton-phonon and superconducting phases truly coexist in the region coupling term the formation of excitons happens inde- of overlap in figure 5. In a more detailed theory the su- pendently of the lattice deformation. The order of the perconducting order may well be suppressed by the pres- transition between normal state and excitonic insulator ence of the CDW phase. Coexistence cannot be directly does change, transforming from a true first order transi- ruled out however; Kiss et al. show that in conventional tion at ρc + ρd = 1 to a smooth crossover at high dopant charge ordered materials, including the layered transi- levels. tion metal dichalcogenides, CDW and superconducting Including also the Cooper pair formation of equation order can coexist, and that the superconductivity is in (6) in the mean field analysis, we find in figure 7 the fact boosted at points in k-space connected by the CDW emergence of a dome shaped superconducting phase sur- ordering vector44. rounding the point where the CDW transition reaches 5

sition temperatures as a function of the exciton binding energy (see figure 7), it becomes clear that the excitons play an essential part in this mechanism, strengthening the role of the phonons via the renormalisation of their spectral properties, and enhancing both the CDW order and the superconducting TC . In the present model the charge order is once again cut off at very high values of V due to the formation of the spatially homogeneous excitonic insulator.

FIG. 6: The mean field phase diagram arising from equations B. Pressure (5) and (6), for various doping levels. The total charge in- creases in the direction of the arrows, and takes on the values ρc + ρd = 1.00, 1.05 and 1.10. The addition of extra electrons Superconductivity has also been observed in pure 26 pushes the formation of CDW order to higher values of the TiSe2 under pressure . Based on differences in the be- electron-phonon coupling, without affecting the onset of the haviour of the Hall coefficient and the sensitivity to mag- excitonic insulator phase. netic fields it has been argued that the nature of the superconducting phase in this case differs from that of the Copper doped compound. In the model of equations (5) and (6) we can simulate the effects of pressure by assuming that its main influence on the electronic struc- ture arises from a change in the orbital overlap inte- grals. We thus introduce the rescaled hopping param- ′ ′ eters t = t0(1 + Peff) and t = t0(1 + cPeff) with c a constant of order 1. Due to the simplified nature of the model we cannot make a quantitative comparison with the experimentally observed behaviour under pressure, FIG. 7: The transition temperature in the mean field solution but we can use the effective pressure Peff to qualitatively of equations (5) and (6), in the presence of electron doping: study its effect on the interplay between exciton forma- left as a function of the total electron density ρc + ρd, at fixed values for the exciton binding energy (V/t = 0.65) and tion and electron-phonon coupling. the electron-hole coupling (α/t = 0.032); and right as a func- Figure 8 shows the evolution of the CDW phase as a tion of the exciton binding energy V/t, at a fixed value of the function of effective pressure. It is clear that the excitonic total electron density (ρc + ρd = 1.1). insulating phase and the charge ordered phase are pushed As in figure 5, the transition temperatures of both the super- to higher values of the excitonic binding energy and the conducting and the charge ordered phase are enhanced by an electron-phonon coupling respectively. At the same time, increase in excitonic binding energy, indicating the coopera- Cooper pair formation in the regions directly adjacent tion between the exciton driven renormalisation of the phonon to the charge ordered phase allow superconductivity to properties and the phonon mediated pairing. arise as the CDW is suppressed. For any particular set of values for the parameters α and V , this implies that the typical behaviour upon tuning both temperature and zero temperature. In the present mean field treatment effective pressure is as shown in figure 9. Starting from the superconductivity extends to relatively high temper- the CDW phase, a superconducting dome again arises atures, but taking into account fluctuation effects should around the point at which the CDW transition reaches reduce the extent of the superconducting phase consid- zero temperature. As before, the superconducting phase erably. The positioning of the dome just beyond the seems to extend to high values of both temperature and point where the disappearing charge order should give effective pressure because of the mean field nature of the rise to a quantum phase transition can be interpreted present treatment, but is expected to be considerably as indicating that the superconductivity in this region is reduced upon taking into account fluctuation effects. mediated by the critical fluctuations associated with the The pair formation in the superconductor under pres- parent CDW phase. As before, the pairing “glue” thus sure is again mediated by the same cooperating exciton- consists of phonons, renormalised by the presence of ex- phonons that also stabilise the parent CDW phase. The citons. The addition of dopant electrons serves the dual excitonic presence and the associated phonon renormal- purpose of both directly augmenting the density of car- isation enhances both the charge order at ambient pres- riers available for Cooper pair formation and of driving sure, and the superconducting order at higher pressures, the system towards the point where the hybrid exciton- as can be seen in figure 9. Even though the mean field phonon modes become soft and can mediate the attrac- treatment given here is expected to largely overestimate tive interaction between the conduction electrons. By the transition temperature, the qualitative trend of in- plotting the superconducting and charge ordering tran- creasing TC as V is increased and the quantum critical 6

hole balance, which may explain the observed differences in the behaviours of the Hall coefficient. Within the present model this does not affect the pairing mechanism in the two regimes though, and the nature of the critical fluctuations is the same throughout. The superconduct- ing dome in the doped compound can then be contin- uously connected to the dome in the pure system un- der pressure, creating a ‘tunnel’ of superconductivity in the doping-pressure-temperature phase diagram, as sug- 26 gested by Kusmartseva et al. for the case of TiSe2 .

FIG. 8: The mean field phase diagram arising from equations (5) and (6), for various values of the effective pressure Peff. Pressure increases in the direction of the arrows, and takes IV. CRITICAL FLUCTUATIONS on the values Peff = 0.0, 0.5, 0.66, 0.75, 1.0 and 1.25. The application of effective pressure affects both the transition The exponent of the resistivity in the normal (non- into the CDW phase and the onset of the excitonic insulator superconducting) state of TiSe2 has been observed to phase, pushing them to higher values of the electron-phonon undergo a dip in the pressure range surrounding the crit- coupling and the exciton binding energy respectively. ical point at which the CDW order fully disappears26. Both at higher and lower pressures the low temperature resistivity scales as ρ T 3, while at the critical pressure the exponent is seen to∝ be reduced towards ρ T 2.5. In fact, the exponent may well be reduced even further,∝ as it is not clear from the data presented in reference [26] whether the reduced exponent has fully converged to its minimal value close to the critical point. Further de- tailed experiments at low temperatures close to the crit- FIG. 9: The transition temperature in the mean field solu- ical pressure will be required to settle the precise value tion of equations (5) and (6): left as a function of effective of the resistivity exponent in this regime. Regardless of pressure Peff, at fixed values for the exciton binding energy its exact depth however, the presence of a dip in the ex- (V/t0 = 0.65) and the electron-hole coupling (α/t0 = 0.035); ponent is suggestive of the influence of quantum critical and right as a function of the exciton binding energy V/t0, fluctuations, like the ones observed in nearly magnetic at a fixed value of the effective pressure Peff = 0.5 and materials, heavy fermion compounds and nearly ferro- α/t0 = 0.04. 2,32–39 Notice that the CDW transition becomes first order at high electric systems . Studying these and related crit- values of the effective pressure, and that as in figure 5, the ical exponents is an important experimental tool in un- transition temperatures of both the superconducting and the derstanding both the nature of the critical modes and the charge ordered phase are enhanced by an increase in exciton properties of the quantum critical point. binding energy, indicating that the cooperation between the To calculate the influence of the phononic and excitonic exciton driven renormalisation of the phonon properties and critical modes on the resistivity in our model system, we the phonon mediated pairing persists under pressure. assume that a description in terms of a Boltzmann trans- port equation will be applicable. Notice that this implies neglecting any fully quantum mechanical effects, like f.e. point approached, is expected to survive. weak localisation. Close to the quantum critical point, The nature of the critical fluctuations mediating super- the collision term in the Boltzmann equation is assumed conductivity in this parameter range is thus no different to be dominated by the scattering of conduction electrons from that encountered in the electron doped case. In off the bosonic critical modes associated with the (hid- both systems an exciton-enhanced phonon mechanism is den) zero temperature CDW transition. These modes are the driving force behind both the CDW and the super- expected to be the hybrid exciton-phonon excitations dis- conducting pairing. The two cases may however vary cussed in the previous sections. In fact, the main mecha- considerably in the values of their superconducting tran- nism for scattering conduction electrons, both in the case sition temperatures, which depend both on the pairing of excitons and of phonons, consists of the formation of mechanism (i.e. the effective coupling strength) and the local distortions in the crystal polarisation field. Scat- electronic density of states. The augmented carrier den- tering from such polarisation clouds is well understood sity in the doped compound may thus be expected to in terms of the Fr¨ohlich Hamiltonian45–47. Although the support the higher transition temperature. The trans- value of the coupling strength may depend on the extent port in the normal (high temperature) state is also af- of the contribution from each mode, the k-dependence fected differently in the two regimes: whereas doping ex- of the resulting matrix elements does not distinguish be- plicitly introduces a surplus of electron-like (rather than tween excitonic or phononic scatterers. We can thus treat hole-like) carriers, pressure retains the overall electron- the excitons and phonons on an equal footing, and in the 7 following we will refer to both, as well as to their combi- developing a local minimum in their dispersion at the or- nation, simply as the bosonic scattering modes. dering wavevector. Rewriting k′ k = Q + q, with q − Following the standard arguments given for example small compared to the ordering wavevector Q, we find by Ziman48, we write the stationary electron distribution ω a + bq2, so that the scattering matrix elements be- 0 0 0 ∼ as fk = fk Φk ∂fk/∂ǫk, where fk is the Fermi-Dirac come nearly constant in q in this regime. Further sim- − · distribution and Φk is a function to be determined. The plification can be achieved by rewriting the momentum resistivity arising from the linearised collision term in the integrals as dk d2k dǫ . After also introducing ∝ k Boltzmann equation is then given by: ǫk′ = ω + ǫkR, the integralR R over ǫk can be performed ex- plicitly, and the integrations over constant energy sur- 51,52,54 1 ′ 2 faces are simplified to yield : ρ dkdk (Φk Φk′ ) Pk,k′ . (7) ∝ T Z − 1 ρ dqdω q2ωn(ω) (n(ω)+1) Im χ (q + Q,ω) . Here Pk,k′ is the scattering rate of electrons with mo- ∝ T Z mentum k into a state with momentum k′. Notice that (9) this expression assumes the bosonic scatterers to be in their equilibrium distribution. The function Φk is a vari- Close to the quantum critical region, the bosonic ational function whose shape is normally well approxi- modes become dissipative, and the susceptibility can be mated by the ansatz Φk k E, with E the applied approximated in analogy to the overdamped harmonic 34 electric field48–50. Although∝ the· lack of isotropy in our oscillator , model system should in principle be reflected in the vari- ω ational function, it has been pointed out before that at −1 −1 χ(q + Q,ω) = χq+Q 1 i −1 2 . (10) very low temperatures impurity scattering will typically  − γ(χ + cq ) be strong enough to completely determine the electron Here we introduced both the static susceptibility χq+Q distribution, while the resistivity is set by the sum of the 51 and the uniform susceptibility χ. Notice that the relax- terms due to impurity and the critical mode scattering . ation rates γ are nearly momentum independent for scat- In this regime the standard ansatz for Φk thus remains tering processes at the finite wavevector q +Q, regardless applicable. of the microscopic damping mechanism involved. Insert- The scattering rate in equation (7) can be related to ing this form into equation (9), the energy integral can the dynamic susceptibility using Fermi’s golden rule and 48,49 now be evaluated. Its low temperature expansion is well the fluctuation-dissipation theorem : known34, and to first order in T the resistivity can be shown to be: ′ 2 0 0 ′ Pk,k Mk,k′ fk 1 fk′ n (ǫk ǫk) ∝ − − −1 ′ −1 2 Im χ
(k k,ǫk′ ǫk) . (8) 2 γ χ + cq /T − − ρ dq q . (11) ∝ Z 1+3 /π [γ (χ−1 +
cq2) /T ] Here n(ω) is the distribution function for the bosonic scattering mode, while Mk,k′ are the scattering matrix Depending on how closely the quantum critical point elements. Because both excitons and phonons scatter is approached, this expression can in principle give rise conduction electrons through a distortion in the local po- to two different temperature dependencies of the resis- larisation field, these matrix elements are expected to be tivity. In the regime near enough to the transition to 2 ′ 2 render the bosonic modes overdamped, but far enough well approximated by Mk k′ k k /ωk′−k, with ω , ∝ | − | to ensure that χ−1 is finite and the integrand of equation the dispersion of the bosonic scattering mode52. (11) non-singular, dimensional analysis suffices to show Away from the quantum critical region, where the sus- that the resistivity will scale as ρ T 2. At the critical ceptibility has no particular features, the expressions of pressure on the other hand, the uniform∝ susceptibility equations (7) and (8) are of the same generic form as becomes identically zero, and the integral is dominated the corresponding equations in for example the elemen- by its low-q, singular part. Introducing x =3γcq2/(πT ), tal transition metals with their partially occupied s and and cutting off the integration at a temperature indepen- d bands. In fact, it was pointed out by Wilson that in- dent point x , we then find34: terband transitions in these materials, which necessar- c ily involve finite momentum bosonic modes, will typi- dx 3 ρ (T )3/2 . (12) cally give rise to a ρ T temperature dependence of ∝ Z √x the resistivity53. The∝ normal state resistivity observed in TiSe2 far from the transition can thus be straightfor- Whether both of these regimes are in fact realised in wardly interpreted as a consequence of the scattering of TiSe2, and whether they occupy any appreciable portion conduction electrons between the hole pocket at the cen- of phase space, depends largely on the detailed proper- tre of the first Brillouin zone and the electron pockets at ties of the involved bosonic modes and their couplings its boundary. to the conduction electrons. As the transition is ap- As the transition is approached, χ(q,ω) starts to de- proached however, the softening of the boson spectrum, velop a singularity and the bosonic modes begin to soften, the increased damping and the rising singularity in the 8 susceptibility will lead to a renormalised temperature de- by imposing an altered ratio of exciton binding energy to pendence of the resistivity. In the regime precisely above electron-phonon coupling V/α in stretched thin films on the quantum critical point, the expression of equation a mismatched substrate. (12) should thus be applicable at very low temperatures. The critical fluctuations associated with the quantum The fact that the experimentally observed temperature critical point below the superconducting dome are exper- dependence does not match the ρ T 1.5 resistivity found ∝ imentally accessible through their influence on the tem- here may be due to the superconducting dome prevent- perature dependence of various thermodynamic quanti- ing measurements to be made at sufficiently low temper- ties. The observed dip in the exponent of the resistivity atures. The reported dip in resistivity may then be inter- of TiSe2 under pressure may thus be interpreted as a con- preted as the onset of a crossover from the normal state sequence of the presence of such fluctuations. We have conductivity towards the full quantum critical behaviour shown in this paper that regardless of the precise origin at even lower temperatures. of the critical fluctuations (whether they are excitons, phonons or a combination thereof), it can be expected on very general grounds that the ρ T 3 law which is a V. CONCLUSIONS consequence of the indirect coupling∝ between two bands in the normal state, will be reduced to a ρ T 1.5 depen- The superconducting phase which has been observed in dence at very low temperatures close to the∝ critical point TiSe2 to arise either upon electron-doping the material because of the critical fluctuations. Depending on the ef- or by applying hydrostatic pressure is likely to be me- fective damping rates of the modes as the critical region diated by the same bosonic modes that are responsible is approached, there may also be an additional region of for the charge density wave order in the pristine parent ρ T 2 behaviour between these two extremes. Based on compound at ambient pressure. Because the charge order these∝ findings, we suggest that the reported dip in the instability in pure TiSe2 is known to involve both exciton resistivity exponent may not have reached its absolute formation and a strong electron-phonon coupling, it may minimum yet, and should be interpreted as an indication be expected that the superconducting phase will likewise that critical fluctuations play an important part in the involve a combination of both excitonic and phononic physics of TiSe2 in the investigated pressure range, but contributions. Within the simplified model discussed in that the full quantum critical behaviour will emerge only this paper, it has been shown that the presence of ex- at lower temperatures. Measurement of related critical citons in a system already prone to lattice deformation exponents, such as those of the specific heat or suscepti- due to a strong electron-phonon coupling may lead to an bility may help to shed further light on this issue. enhancement of the CDW order. Upon suppression of that order, either by raising the effective chemical po- If the main difference between the two supercon- tential to inject more electrons or by altering the model ducting domes previously observed in TiSe2 is the de- parameters to simulate an effective pressure, a supercon- tailed balance between the effects of exciton formation ducting dome appears around the point where the CDW and electron-phonon coupling, as suggested in this pa- transition reaches zero temperature. In both cases the per, they should be connected by a continuous ‘tunnel’ of superconducting order in the pressure-doping phase superconducting order may be enhanced by enhancing 26 the excitonic binding energy, thus revealing the support- diagram . Completion of the experimental phase dia- ing role the excitons play also in this phase. In fact, gram by mapping out the pressure dependence of var- the cooperation between exciton formation and electron- iously doped compounds, as well as considering other phonon coupling seems to be ubiquitous throughout the (electron or hole) dopants may thus provide further sup- emerging phase diagram. The extent of hybridisation be- port for the hybrid exciton-phonon picture presented tween excitonic and phononic modes may be sensitive to here. The maximum critical temperature of this super- the experimental route taken, but the superconducting conducting phase will be enhanced in areas where the Cooper pairs are always held together by some combina- exciton binding energy effectively cooperates with the tion of the two modes. phonons, and where a large electronic density of states is To assess precisely which portion of the superconduct- available to utilise these cooperating bosonic modes. The amount of enhancement however, will be sensitive to the ing ‘glue’ in TiSe2 consists of excitons, more detailed the- oretical modelling (evaluating the full frequency depen- detailed material properties, and a substantial increase of dent Eliashberg equations for a more realistic model sys- TC through this mechanism will likely only be achievable tem) as well as new experimental work (focusing specif- with intricate fine tuning of all involved parameters. ically on the differences between excitons and phonons) will be required. Possible experimental signatures of the Acknowledgements influence of excitons may be sought for example by alter- The authors would like to thank the EPSRC and Jesus ing the ratio of exciton binding energy to phonon energy and Homerton Colleges of the University of Cambridge V/ω through the use of isotope effects, or more directly for financial support. 9

1 J. R. Waldram, Superconductivity of Metals and Cuprates (1984). (Taylor & Francis, 1996). 28 R. H. Friend and A. D. Yoffe, Adv. Phys. 36, 1 (1987). 2 N. D. Mathur, F. M. Grosche, S. R. Julian, I. R. Walker, 29 S. Y. Li, G. Wu, X. H. Chen, and L. Taillefer, Phys. Rev. D. M. Freye, R. K. W. Haselwimmer, and G. G. Lonzarich, Lett. 99, 107001 (2007). Nature 394, 39 (1998). 30 R. A. Jishi and H. M. Alyahyaei, Phys. Rev. B 78, 144516 3 J. F. Schooley, W. R. Hosler, E. Ambler, J. H. Becker, (2008). M. L. Cohen, and C. S. Koonce, Phys. Rev. Lett. 14, 305 31 J. van Wezel, P. Nahai-Williamson, and S. S. Saxena, Phys. (1965). Stat. Sol. b 247, 592 (2010). 4 B. Sipos, A. F. Kusmartseva, A. Akrap, H. Berger, 32 S. Sachdev, Quantum Phase Transitions. (Cambridge Uni- L. Forr´o, and E. Tutis, Nat. Mater. 7, 960 (2008). versity Press, 2001). 5 W. A. Little, Phys. Rev. 134, A1416 (1964). 33 T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 6 V. L. Ginzburg, Sov. Phys. JETP 20, 1549 (1965). netism. (Springer, Berlin, 1985). 7 O. V. Dolgov and E. G. Maksimov, Sov. Phys. Usp. 25, 34 G. G. Lonzarich, in Electron., edited by M. Springford 688 (1982). (Cambridge University Press, 1997). 8 K. Rossnagel, L. Kipp, and M. Skibowski, Phys. Rev. B 35 G. G. Lonzarich and L. Taillefer, J. Phys. C 18, 4339 65, 235101 (2002). (1985). 9 T. E. Kidd, T. Miller, M. Y. Chou, and T.-C. Chiang, 36 S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, Phys. Rev. Lett. 88, 226402 (2002). R. K. W. Haselwimmer, M. J. Steiner, E. Pugh, I. R. 10 M. H. Whangbo and E. Canadell, J. Am. Chem. Soc. 114, Walker, S. R. Julian, P. Monthoux, et al., Nature 406, 9587 (1992). 587 (2000). 11 H. P. Hughes, J. Phys. C: Solid State Phys. 10, L319 37 G. G. Lonzarich, P. Monthoux, and D. Pines, Nature 450, (1977). 1117 (2008). 12 J. A. Wilson, Phys. Stat. Sol. (b) 86, 11 (1978). 38 R. P. Smith, M. Sutherland, G. G. Lonzarich, S. S. Sax- 13 J. A. Wilson, A. S. B. Jr., F. J. D. Salvo, and J. A. Ditzen- ena, N. Kimura, S. Takashima, M. Nohara, and H. Takagi, berger, Phys. Rev. B. 18, 2866 (1978). Nature 455, 1220 (2008). 14 K. Motizuki, N. Suzuki, Y. Yoshida, and Y. Takaoka, Solid 39 S. E. Rowley, L. J. Spalek, R. P. Smith, M. P. M. Dean, State Comm. 40, 995 (1981). G. G. Lonzarich, J. F. Scott, and S. S. Saxena (2009), 15 N. Suzuki, A. Yamamoto, and K. Motizuki, J. Phys. Soc. arxiv:Cond-mat 0903.1445. Jpn. 54, 4668 (1985). 40 C. Monney, E. F. Schwier, M. G. Garnier, C. Battaglia, 16 J. van Wezel, P. Nahai-Williamson, and S. S. Saxena, Eu- N. Mariotti, C. Didiot, H. Cercellier, J. Marcus, H. Berger, roPhys. Letters 89, 47004 (2010). A. N. Titov, et al., ArXiv:cond-mat p. 0912.5283v1 (2009). 17 J. van Wezel, P. Nahai-Williamson, and S. S. Saxena, Phys. 41 C. M. Varma, S. Schmitt-Rink, and E. Abrahams, Solid Rev. B 81, 165109 (2010). State Comm. 62, 681 (1987). 18 M. Holt, P. Zschack, H. Hong, M. Y. Chou, and T.-C. Chi- 42 P. B. Littlewood, Phys. Rev. B 42, 10075 (1990). ang, Phys. Rev. Lett. 86, 3799 (2001). 43 N. Nakai, M. Ichioka, and K. Machida, J. Phys. Soc. Japan 19 J. A. Wilson, Solid State Comm. 22, 551 (1977). 71, 23 (2002). 20 A. Zunger and A. J. Freeman, Phys. Rev. B 17, 1839 44 T. Kiss, T. Yokoya, A. Chainani, S. Shin, T. Hanaguri, (1978). M. Nohara, and H. Takagi, Nat. Phys. 3, 720 (2007). 21 G. Li, W. Z. Hu, D. Qian, D. Hsieh, M. Z. Hasan, E. Mo- 45 H. Fr¨ohlich, H. Pelzer, and S. Zienau, Phil. Mag. 41, 221 rosan, R. J. Cava, and N. L. Wang, Phys. Rev. Lett. 99, (1950). 027404 (2007). 46 H. Fr¨ohlich, Adv. Phys. 3, 325 (1954). 22 C. Monney, H. Cercellier, F. Clerc, C. Battaglia, E. F. 47 C. G. Kuper and G. D. Whitfield, eds., Polarons and exci- Schwier, C. Didiot, M. G. Garnier, H. Beck, P. Aebi, tons : Scottish Universities’ Summer School, 1962 (Oliver H. Berger, et al., Phys. Rev. B 79, 045116 (2009). and Boyd, Edinburgh, 1963). 23 H. Cercellier, C. Monney, F. Clerc, C. Battaglia, L. De- 48 J. Ziman, Electrons and Phonons. (Clarendon, Oxford, spont, M. G. Garnier, H. Beck, P. Aebi, L. Patthey, 1960). H. Berger, et al., Phys. Rev. Lett. 99, 146403 (2007). 49 D. G. Dungate, Ph.D. thesis, University of Cambridge 24 E. Morosan, H. W. Zandbergen, B. S. Dennis, J. W. G. (1990), phd.17451. Bos, Y. Onose, T. Klimczuk, A. P. Ramirez, N. P. Ong, 50 R. Hlubina and T. M. Rice, Phys. Rev. B 51, 9253 (1995). and R. J. Cava, Nat. Phys. 2, 544 (2006). 51 K. Ueda, J. Phys. Soc. Jpn. 43, 1497 (1977). 25 E. Morosan, L. Li, N. P. Ong, and R. J. Cava, Phys. Rev. 52 See for exmaple: A. H. Wilson, The Theory of Metals. B 75, 104505 (2007). (Cambridge University Press, Cambridge, 1953). 26 A. F. Kusmartseva, B. Sipos, H. Berger, L. Forr´o, and 53 A. H. Wilson, Proc. R. Soc. Lond. A 167, 580 (1938). E. Tutiˇs, Phys. Rev. Lett. 103, 236401 (2009). 54 K. Ueda and T. Moriya, J. Phys. Soc. Jpn. 39, 605 (1975). 27 P. C. Klipstein and R. H. Friend, J. Phys. C 17, 2713