Temperature effects in excitonic condensation driven by the lattice distortion
Abstract
The stability of the excitonic condensation at low temperature driven by a coupling of electrons to vibrational degrees of freedom in semimetal two-dimensional electronic system is discussed. In the framework of the unrestricted Hartree-Fock approximation, we derive a set of equations to determine both the excitonic condensate order parameter and lattice displacement self-consistently. By lowering temperature we find out a semimetal-insulator transition in the system if the coupling is large enough. The insulating state typifies an excitonic condensation accompanied by a finite lattice distortion. Increasing temperature, both excitonic condensate order parameter and the lattice distortion decrease and then disappear in the same manner. Microscopic analysis in momentum space strongly specifies that the excitonic condensate driven by the lattice distortion favours the BCS type.
pacs:
71.45.Lr, 71.35.Lk, 63.20.kk, 71.30.+h, 71.28.+dI Introduction
Electron-hole pairing or excitonic condensation recently has become one of the most attractive subjects in many-particle physics Snoke and Kavoulakis (2014); Combescot and Shiau (2015). At sufficiently low temperature, high density excitons may condense and the system stabilizes in an insulating state with respect to the spontaneous formation of a new macroscopic phase-coherent quantum state called an excitonic insulator (EI). Investigating the EI state has been intensively focused in the literature but mainly on the purely electronic manners Batista (2002); Farkašovský (2008); Phan et al. (2010); Zenker et al. (2012). In doing so the coupling of electrons or excitons to the phonon was completely neglected.
However, recently, several experiments have opened an issue that the lattice distortion is non-negligible to anticipate the EI phase in the quasi-two dimensional transition metal dichalgogenide 1-TiSe. In this material, the charge density wave (CDW) has been manifested to accompany with the weak periodic lattice distortion Di Salvo et al. (1976). At low temperature, photoemission signatures indicate that the exciton condensate strongly associates to the CDW state Monney et al. (2010). By measuring the thermodynamic properties of TmSeTe, Wachter and co-workers proposed a strong phonon-exciton couple in the excitonic condensate phase, forming exciton-polaron quasiparticles Wachter (2001); Wachter et al. (2004). In their studies, the heat conductivity shows a divergence for . This anomaly is analogous to the observation for He II, a typical superfluidity below K where the heat is being carried by phonons Tilley and Tilley (1990). In TmSeTe, at an appropriate pressure, an excitonic bound state of a hole at the -point and a electron at the -point can be created. An interaction of electron/hole with a phonon thus needs to assist that transfer. Wachter (2001) At sufficiently low temperature, those excitons condense to a superfluid forming the EI state Wachter (2001); Wachter et al. (2004). Without any doubt, lattice distortion or phonon effects are extremely important in this kind of material, particularly, in studying the EI state. The electron/hole-phonon interaction therefore seems to be non-negligible and needs to be considered thoroughly.
On the theoretical side, the lattice distortion causing the EI state has been studied intensively however only for the ground state, i.e., at zero temperature Kaneko et al. (2013); Zenker et al. (2013); Phan et al. (2013). In general, as a kind of superfluidity, the EI state possibly occurs at finite temperature. At high temperature it might be deformed by thermal fluctuations. Studying the influence of temperature on the EI state therefore is an important effort. The phase diagram specifying the EI region in the temperature-pressure plane has been experimentally measured, which emphasizes that the excitons condense only at low temperature with intermediate pressure around the semimetal-semiconductor transition Neuenschwander and Wachter (1990); Wachter et al. (2004). Without the electron/hole-phonon coupling, the temperature-pressure phase diagram of the EI state has been considered theoretically in the extended Falicov-Kimball model Ihle et al. (2008); Zenker et al. (2010a, b).
By considering the tight-binding formalism, the temperature effects in excitonic condensate exerting a force on the lattice generating periodic ionic displacements in 1-TiSe have been discussed Monney et al. (2011). Of course, here, the mean-field form of the EI order parameter depending on temperature has been assumed from the beginning. Moreover, the photoemission temperature dependence measured for -TiSe has shown that it fits quite well with the mean-field form for temperatures below the EI transition temperature Monney et al. (2010). In the present work, we intend to develop the unrestricted Hartree-Fock approximation, a kind of the mean-field approach but allowing for decoupling with respect to the excitonic order parameter, to a two-dimensional two-band - electron model with a coupling to the phonon degrees of freedom. This coupling results in a ‘hybridization’ between and electrons. In the one-dimensional case, this model has been studied. It shows that at zero temperature the EI exists with respect to a lattice displacement happening only if the electron-phonon coupling is larger than a critical value Phan et al. (2013).
The paper is organized as follows. In Sec. 2 we introduce the two-band - electronic model with a coupling to the phonon. The theoretical approach is outlined in Sec. 3, where the unrestricted Hartree-Fock approximation has been developed to the specific model mentioned in Sec. 2. Sec. 4 presents and discusses in detail numerical results. Our main conclusions can be found in Sec. 5.
Ii Two-band electron-phonon interaction model
The system with two spinless electronic bands involving electron-phonon coupling is modeled by the following Hamiltonian
(1) | |||||
where () and () correspond to creation (annihilation) operators of and spinless electrons carrying momentum . The kinetic energy of a phonon with dispersionless energy is given in the third term in Eq. (1) where and are phonon creation and annihilation operators at momentum . Here, the electronic excitation energies are given by
(2) |
where represents the local part of the () electron excitation, the next term , with
(3) |
accounts for the nearest-neighbor hopping in a 2D lattice, and is the chemical potential. The last term in Eq. (1) addresses a local electron-phonon interaction (with coupling constant ), written in -space.
The Hamiltonian given in Eq. (1) is identical to the Holstein model in which by the canonical transformations, the effective electron-hole attraction can be performed Fröhlich (1952). Electron-hole bound states or excitons might therefore exist due to the coupling with the lattice - a CDW state. Like Cooper pairs in superconductors, at low temperature these excitons would prefer to form a superfluid state Combescot and Shiau (2015).
To mimic the situation of TmSeTe, where the quasilocalized state has its maximum at the -point and the strongly dispersive state has its minimum at the point, we choose and , whereas is chosen as the unit of energy. It characterizes an indirect coupling assisting by the transfer phonon. At sufficiently low temperature, the bound pairs with finite momentum might condense, indicated by a nonzero value of the order parameter, i.e.,
(4) |
where in two dimensions. The bound pair or exciton in this case is considered to have nonzero center-mass momentum which is distance between and points. That bound pairs raise the CDW state accompanied by a phonon Wachter and Bucher (2013). In the following, we also consider
(5) |
as the excitonic condensation order parameter.
Iii Unrestricted Hartree-Fock approximation
As a kind of Hartree-Fock approximation, the unrestricted Hartree-Fock approximation allows decoupling with respect to the off-diagonal expectation values Schneider and Czycholl (2008), such as the excitonic order parameter, , in our case. In this respect, we introduce the fluctuation operator for an arbitrary operator , and write the electron-phonon interaction operator in Eq. (1) as
(6) | |||
Assuming that the fluctuations are small, the first term on the right hand side in Eq. (III) can be eliminated. In this case, the Hamiltonian in Eq. (1) reduces to the so-called unrestricted Hartree-Fock Hamiltonian which can be separated into two parts
(7) |
where the electronic part reads
(8) |
and
(9) |
is the phononic one. Note here that the additional constant has been neglected in Eq. (7). In Eq. (8) reads
(10) |
indicating the hybridization between the and electrons, on one hand it expresses the effective bound state of electron-hole pairs, on the other hand it indicates the CDW instability like in the Holstein model in the frozen phonon approximation. The factor in Eq. (9) is given by
(11) |
From the expressions in Eqs. (10-11) and the Hamiltonians in Eqs. (8) and (9) we can easily realized that both and are mutually dependent. Moreover, since , the field contribution has been replaced by . Therefore a finite lattice displacement approximating to would give rise to the formation of a CDW state connected to a doubling of the lattice unit cell.
In order to diagonalize the unrestricted Hartree-Fock Hamiltonian written in Eq. (7), firstly, we define a new phonon operator
(12) |
to diagonalize the phononic part in Eq. (9). Meanwhile, the electronic part can be diagonalized itself by using a Bogoliubov transformation, where the new quasi-particle fermionic operators read
(13) | ||||
(14) |
Here, the prefactors and are chosen to satisfy . Then finally, we are led to a completely diagonalized Hamiltonian
(15) |
where the electronic quasiparticle energies read
(16) |
and the prefactors addressed in Eqs. (13-14) are given by
(17) | ||||
(18) |
with
(19) |
The quadratic form of Eq. (15) allows to compute all expectation values formed with , resulting in
(20) | |||
(21) | |||
(22) |
Here is the Fermi-Dirac distribution function and is the inverse of the temperature.
To calculate in Eq. (10) we note here that and from Eq. (12) one obtains
(23) |
Then finally we obtain the amplitude of the lattice displacement in the EI state for a single center-mass momentum
(24) |
Let us also consider the electronic one-particle spectral functions. For the -electron, its spectral function can be found following the definition
(25) |
From the diagonal form of the effective Hamiltonian in Eq. (15) one derives
(26) |
In the same way, we can construct the spectral function of the -electron, that reads
(27) |
Note that according to the spectral function we can evaluate the corresponding density of state by taking a summation over all momenta in the first Brillouin zone.
Iv Numerical results
In this section, we present numerical results to discuss the influence of temperature on the excitonic condensate of the two-band electrons involving the electron-phonon interaction. For the two dimensional system consisting of lattice sites, the numerical results are obtained by solving self-consistently Eqs. (10), (11), (22), and (23) starting from some guess values for and . In what follows, we fix and consider a half-filled band case, i.e.,
(28) |
where
(29) |
and
(30) |
are respectively the densities of valance and conduction electrons. The chemical potential therefore has to be adjusted in such a way that Eq. (28) is satisfied. The difference between the on-site energies of a -electron and an -electron is .
As a first step, we discuss the ground state of the electron-phonon systems modeled by the Hamiltonian in Eq. (1). Figure 1 shows us a dependence of the EI order parameter and the lattice displacement on the electron-phonon coupling at zero temperature. We see that and are intimately related. For a given value of phonon frequency , and appear to be nonzero only if the exciton-phonon coupling is larger than a critical value . This increases when increasing the phonon frequency. This feature has been addressed in the one-dimensional case in which a linear relation between a square critical coupling and the phonon frequency has been confirmed both in numerical and analytical calculations Phan et al. (2013). In the case of , the EI/CDW state exhibits a finite lattice distortion. In the contrast, i.e., with , the system settles in a semi-metal state with an undistorted lattice structure. This semimetal EI/CDW transition is Kosterlitz-Thouless type, typically observed in general two-dimensional systems Kosterlitz and Thouless (1973). Note here that in the one-dimensional case, the model written in Eq. (1) also consists the Kosterlitz-Thouless transition form of the semimetal EI transition Phan et al. (2013).
Figure 1 indicates that our situation can be split into two different regimes characterized by the phonon frequency . In the adiabatic regime (), for instance , the lattice displacement is larger than the EI order parameter, whereas in the opposite case, i.e., in the anti-adiabatic regime () the lattice displacement is smaller and goes more slowly than the EI order parameter. For that reason, the temperature effect on the EI/CDW semimetal transition needs to be discussed separately in both regions. In the following we focus on two typical values of , and .
In Figure 2(a) the EI order parameter and the lattice displacement are shown as functions of temperature for some values of at . Obviously, the EI state stabilizes at low temperature. That is analogous to the case of superconductivity. The EI order parameter decreases if temperature is increased. It completely disappears at a critical temperature, . here therefore stands for an EI/CDW transition temperature. The EI transition temperature is weakened by lowering the coupling of the electron-phonon interaction. This once more reminds us of the BCS theory in which the condensation of Cooper pairs is driven by the electron-phonon interaction. The BCS type condensation of electron-hole pairs in the present problem will be discussed in detail in the remained figures. The temperature dependence of shown in Figure 2(a) here fits quite well with the recent experimental observation in the quasi-two dimensional -TiSe system Monney et al. (2010). That good description of the mean-field approach strongly indicates the macroscopic condensation of coherent excitons below the critical temperature. Above the critical temperature, the strong fluctuations of incoherent electron-hole pairs would be significant Monney et al. (2010)
In Figure 2(a) we show also the lattice displacement evaluated for the same set of parameters. As in the ground state (cf. Figure 1), at finite temperature we still recognize the analogous behaviour of the lattice displacement and the EI order parameter. They simultaneously disappear if the temperature is larger than the EI/CDW transition temperature. Below the transition temperature, the lattice distortion occurs according to the stability of the exciton condensation. This temperature dependence of the displacement in the effect of the electron-phonon interaction has been discussed in the tight-binding approximation where the mean-field like temperature dependence of the EI order parameter has been assumed Monney et al. (2011). In our task both lattice displacement and EI order parameter are calculated self-consistently, i.e., no manual form of them needs to be chosen at the beginning. The temperature dependence of the lattice displacement agrees qualitatively with extracted data from neutron diffraction experiments at low temperatures below Di Salvo et al. (1976). That once more confirms the validity of the mean-field approach at low temperature to study the lattice distortion induced by the EI/CDW state. Similarly to Figure 2(a), Figure 2(b) displays the temperature dependence of the EI order parameter and the lattice displacement for some values but in the anti-adiabatic case. In this situation the lattice displacement is always smaller than the EI order parameter in the whole range of temperature. In what below, we focus on the adiabatic regime to discuss the effects of temperature on the microscopic properties in our system. In the anti-adiabatic case we find (not shown here) the same physical scenarios.
Figure 3 shows the renormalized quasiparticle energy bands and along the diagonal direction of the first 2D Brillouin zone, i.e., , for different temperatures at and . Note that we have chosen the semi-metallic situation meaning that in the non-interacting case, both - and -bands overlap. The Fermi surface in this case is large and both types of quasiparticles participate to form the Fermi surface. At low temperature, it shows us that a sufficiently large electron-phonon coupling triggers a gap opening at the Fermi level. The opening of the gap indicates the bound state of the electron-hole pairs. This feature makes manifest the experimental observation in TmSeTe that at large pressure, two - and -bands overlap. Due to the phonon scattering, -holes couple to -electrons to form excitons. At sufficiently low temperature, these excitons condense Wachter (2001). This once again reminds of a similar relevance to the BCS theory of the superconductivity where Cooper pairs are formed. The width of the gap is proportional to the EI order parameter, which therefore decreases when enhancing temperature (cf. Fig. 2). Below the EI/CDW transition temperature electron-hole pairs are formed. Whereas, above the EI/CDW transition temperature, large thermal fluctuations destroy the bound electron-hole pairs and the gap disappears. For a strictly 2D system, of course, the critical temperature for exciton condensation would be zero, but the superfluid properties should survive for temperatures smaller than the Kosterlitz-Thouless transition temperature Zhang and Joglekar (2008); Min et al. (2008); Phan and Fehske (2012).
Next, in Figure 4 we address the momentum dependence of the order parameter in the first Brillouin zone for different temperatures at the same set of parameters given in Figure 3. Note that besides , also characterizes the EI order parameter. For low temperatures and close to the Fermi momentum, is strongly peaked, otherwise is a rather smooth function of . As a matter of course, increasing temperature, the amplitude of the peaks goes down and it becomes insignificant if the temperature is higher than the EI transition temperature. At low temperature, Fermi surface apparently plays an important role to form the electron-hole bound state, indicating the typical BCS-type of the EI stability Phan et al. (2010, 2011).
To understand more about the temperature effects in the electron-hole bound state we now present the wave-vector and frequency resolved single-particle spectral functions associated with the photoemission or inverse photoemission (injection) of and electrons. Figure 5 shows the variation of (left) and (right) for different temperatures as discussed in the last two figures following the high-symmetry diagonal direction in the Brillouin zone.
At zero temperature (see top panels) we find the gap feature opened at the Fermi level (see also Figure 3). In this case, - and -electron states strongly hybridize close to the Fermi energy with large spectral weight transferring. Electron-hole pairs are therefore created and then condense in BCS type. The single-particle excitation gap is suppressed by increasing temperature (see middle panels) and completely disappears when the temperature is larger than the critical value (see bottom panels). The gapless excitation at high temperature indicates that the bounding state of electron-hole pairs is broken. At high temperature, the system settles in the electron hole plasma state.
To analyze the dynamical properties of the electron-hole systems in the full range of temperature we discuss below a scenario of density of state (DOS). In Figure 6 intensity plots are shown for the DOS of the -band (left panel) and -band (right panel) as functions of temperature at and . At high temperature the systems clearly show metallic state with nonzero DOS at the Fermi energy. Lowering temperature, a correlation induced “hybridization” gap opens, indicating long-range order of the non-vanishing polarization and also the lattice displacement. A strong enhancement of the DOS at the upper valence and lower conduction band edges due to the -mixing state reminds us of a BCS-type structure evolving from a semimetallic state with a large Fermi surface above Ihle et al. (2008).
The microscopic investigations of the system depending on temperature are in qualitatively good agreement with experimental observations Cercellier et al. (2007); Monney et al. (2010). Our results indicate that the excitons in the system are condensed in a BCS-like manner and give rise to a CDW state, characterized by the signatures of the EI order parameter and the lattice displacement as discussed before.
V Summary
To summarize, in this paper we have developed the unrestricted Hartree-Fock approximation adapting to the two-band electronic-phonon interaction model to address the influence of temperature on the excitonic condensate and lattice displacement in quasi two-dimensional systems. With an expectation that the electron-hole pair bound state might be formed at low temperature we have derived self-consistent equations permitting us to determine the EI order parameter and the lattice displacement. Numerical results show us that, in the ground state, the EI stability and the lattice displacement are intimately related at large electron-phonon interaction. Below a critical value of the electron-phonon coupling strength the system settles in the metallic state. That happens for both adiabatic and anti-adiabatic limitations. Extending the calculation over a wide temperature range, we find a critical temperature at which the bound state of electron-hole pairs start to be completely destroyed by thermal fluctuations. By inspecting in further detail the momentum dependence of quasi-particle bands, the EI order parameter and photoemission spectral functions, we once more confirm the BCS type condensation of the electron-hole pairs at low temperature. In particular, the photoemission spectral function reveals a pronounced back-folding of the spectral signature in the EI state that directly makes manifest the transition forward a CDW state. Analyzing the temperature dependence of and electron DOSs also exhibits the striking BCS type structure of the EI state evolving from a semimetallic state with a large Fermi surface above . In this way our work has pointed out the prominent role played by the lattice degrees of freedom establishing a charge-density wave in semimetallic systems with weak (indirect) band overlap and in mixed-valent semi-conductors with band gaps comparable to the exciton binding energy, such as quasi-two-dimensional 1-TiSe. The good agreement with experimental observations of our results strongly indicates that the unrestricted Hartree-Fock approach is applicable to consider excitonic condensation inducing the charge density wave state with lattice distortion. Our work therefore gives strong support for exciton condensation as a purely phonontic mechanism responsible for the CDW phase in materials. Combining with the electronic mechanism (due to Coulomb interactions) to consider competition of the EI/CDW state and the lattice distortion will be left to the future.
Acknowledgements.
The authors would like to thank Prof. H. Fehske for valuable discussions. This research is funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.01-2014.05.References
- D. Snoke and G. M. Kavoulakis, Rep. Prog. Phys. 77, 116501 (2014).
- M. Combescot and S.-Y. Shiau, Excitons and Cooper Pairs: Two Composite Bosons in Many-Body Physics (Oxford University Press, New York, 2015).
- C. D. Batista, Phys. Rev. Lett. 89, 166403 (2002).
- P. Farkašovský, Phys. Rev. B 77, 155130 (2008).
- V.-N. Phan, K. W. Becker, and H. Fehske, Phys. Rev. B 81, 205117 (2010).
- B. Zenker, D. Ihle, F. X. Bronold, and H. Fehske, Phys. Rev. B 85, 121102R (2012).
- F. J. Di Salvo, D. E. Moncton, and J. V. Waszczak, Phys. Rev. B 14, 4321 (1976).
- C. Monney, E. F. Schwier, M. G. Garnier, N. Mariotti, C. Didiot, H. Beck, P.Aebi, C. Cercellier, J. Marcus, C. Battaglia, et al., Phys. Rev. B 81, 155104 (2010).
- P. Wachter, Solid State Commun. 118, 645 (2001).
- P. Wachter, B. Bucher, and J. Malar, Phys. Rev. B 69, 094502 (2004).
- D. R. Tilley and J. Tilley, in Superfluidity and Superconductivity, edited by D. F. Brewer (Adam Hilger, Bristol and New York, 1990), p. 35.
- T. Kaneko, T. Toriyama, T. Konishi, and Y. Ohta, Phys. Rev. B 87, 035121 (2013), ibid, 87, 199902(E) (2013).
- B. Zenker, H. Fehske, H. Beck, C. Monney, and A. R. Bishop, Phys. Rev. B 88, 075138 (2013).
- V.-N. Phan, K. W. Becker, and H. Fehske, Phys. Rev. B 88, 205123 (2013).
- J. Neuenschwander and P. Wachter, Phys. Rev. B 41, 12693 (1990).
- D. Ihle, M. Pfafferott, E. Burovski, F. X. Bronold, and H. Fehske, Phys. Rev. B 78, 193103 (2008).
- B. Zenker, H. Fehske, and C. D. Batista, Phys. Rev. B 82, 165110 (2010a).
- B. Zenker, D. Ihle, F. X. Bronold, and H. Fehske, Phys. Rev. B 81, 115122 (2010b).
- C. Monney, C. Battaglia, H. Cercellier, P. Aebi, and H. Beck, Phys. Rev. Lett. 106, 106404 (2011).
- H. Fröhlich, Proc. R. Soc. London A 215, 291 (1952).
- P. Wachter and B. Bucher, Physica B 408, 51 (2013).
- C. Schneider and G. Czycholl, Eur. Phys. J. B 64, 43 (2008).
- J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973).
- C.-H. Zhang and Y. N. Joglekar, Phys. Rev. B 77, 233405 (2008).
- H. Min, R. Bistritzer, J. J. Su, and A. H. MacDonald, Phys. Rev. B 78, 121401(R) (2008).
- N. V. Phan and H. Fehske, New J. Phys. 14, 075007 (2012).
- N. V. Phan, H. Fehske, and K. W. Becker, Europhys. Lett. 95, 17006 (2011).
- H. Cercellier, C. Monney, F. Clerc, C. Battaglia, L. Despont, M. G. Garnier, H. Beck, P. Aebi, L. Patthey, H. Berger, et al., Phys. Rev. Lett. 99, 146403 (2007).