Probing phase coherence via density of states for strongly correlated excitons

We present the calculation of the coherent spectral functions and density of states (DOS) for excitonic systems in the frame of the three dimensional extended Falicov-Kimball model. By using gauge-invariant U(1) transformation to the usual fermions, we represent the electron operator as a fermion attached to the U(1) phase-flux tube. The emergent bosonic gauge field, related to the phase variables is crucial for the Bose-Einstein condensation (BEC) of excitons. Employing the path-integral formalism, we manipulate the bosonic and fermionic degrees of freedom to obtain the effective actions related to fermionic and bosonic sectors. Considering the normal and anomalous excitonic Green functions, we calculate the spectral functions, which have the forms of convolutions in the reciprocal space between bosonic and fermionic counterparts. For the fermionic incoherent part of the DOS we have found the strong evidence of the hybridization-gap in DOS spectra. Furthermore, considering Bogoliubov coherence mechanism, we calculate the coherent DOS spectra. For the coherent normal fermionic DOS, there is no hybridization-gap found in the system due to strong coherence effets and phase stiffness. The similar behavior is observed also for the condensate part of the anomalous excitonic DOS spectra. We show that for small values of the Coulomb interaction, fermionic DOS exhibits a Bardeen-Cooper-Schriffer (BCS) -like double-peak structure. In the BEC region of the BCS-BEC crossover, the double-peak structure disappears totally for both: coherent and incoherent DOS spectra. We discuss also, temperature dependence of DOS functions.


INTRODUCTION
The attractive Coulomb interaction between conduction band electrons and valence band holes plays a crucial role in forming the bosonic electron-hole (e-h) bound pairs 1 (excitons). The Bose-Einstein condensation of the excitons represents a very interesting subject in solid state physics. Mostly, the e-h systems are realized in photoexcited semiconducting materials and the properties of such systems strongly depend on the e-h density, temperature of the system, pressure, and other important physical quantities. There is a number of works, both experimental and theoretical, where all those effects are investigated intensively. [1][2][3][4][5][6] Low density system of excitons behaves like the usual Bose gas, and, at cryogenic temperatures, the BEC transition can be expected. 4,5 On the other hand, the high density system of bound e-h pairs behaves like the system of usual Cooper pairs in superconductors. 7 In this limit we have the Bardeen-Cooper-Schrieffer (BCS) state of e-h pairs. 8,9 Despite many experimental investigations to observe the coherent exciton condensates [10][11][12][13][14] there is not yet definitive evidence for such states. Therefore, an expected BCS-BEC crossover represents actually a fascinating problem typical to the excitonic systems. [15][16][17][18][19][20] Especially, it is interesting from the viewpoint of the difference from the similar crossover in superconductors or trapped atomic Fermi gases. [21][22][23] The transition to the e-h pair condensed phase, in the limit of weak-coupling, is related to the relative motion between electrons and holes, 15 implying the BCS regime and is in contrast to the case of strongcoupling regime, when the BEC state is related to the motion of the center of mass of excitons. 15 In the whole BCS-BEC transition region, the e-h mass difference leads to a large suppression of the BEC transition temperature, which is proved to not be same as the temperature of excitonic pair (EP) formation, and hence the excitonic insulator (EI) phase. 15 Recently, it is shown theoretically that the excitonic insulator and the excitonic condensate are not exactly the same states of the matter. 15,24,25 The author in Refs. 24 and 25 shows from general considerations that in the low density limit of the excitonic pairs, the critical temperature of excitonic condensation should be much smaller than the temperature of EP formation, in contrast to the previous treatments, [16][17][18][19][20] where the EI state is ad hoc associated with the BEC state of excitons as the same. Similarly, in Ref. 15 it is shown that the EI state is an excitonium state, where incoherent e-h bound pairs are formed and, furthermore, at the lower temperatures, the BEC of excitons appears in consequence of the reconfiguration and coherent condensation of the preformed excitonic pairs. Obviously, in the low density limit, the gas of free excitons undergoes the BEC phase transition at the very low temperatures, and in general the BEC temperature transition line is not coinciding with that of EP formation. The Bose condensation of excitonic pairs is possible only when the macroscopic phase coherence is attended by the system. 24 However, the experimental evidence of the existence of two distinct transition temperatures, for the general case of a three dimensional (3D) bulk system, is yet lacking in the literature. In fact, the question, whether a true, "coherent" BEC transition is present in the system of excitons, (in the case of high exciton density limit) is still ambiguous. The experimental proof of it, is a very cumbersome problem, because of dominant role of quantum fluctuations at low temperatures k B T E eh Ry (with E eh Ry being the binding energy of a Mott-Wannier exciton), when very large zero-point oscillations are present.
The continuing growing interest, to the problem of the coherent excitonic condensates motivates us to calculate excitonic spectral functions and density of states, as a direct probe mechanism, to compare the results with the high-resolution studies of angle-resolved photoemission spectroscopy (ARPES) [26][27][28][29] , and with spectral weight measurements on the excitonic materials at the very low temperatures.
The general strategy of our calculations is based on the effective actions method. First, we transform the initial total action of the system to a gauge invariant form, by applying the U(1) gauge transformation to the fermion operators. As a result, the electron appears in the theory like a composite object of that of the fermion with the attached U(1) phase-"flux-tube". The electron factorization in terms of two variables has an unprecedented impact on the whole theory. Then we integrate out the phase variables and we get the effective fermionic action in the theory and we derive a set of self-consistent equations for the EI state. Furthermore, we discuss shortly the results of the quantum rotor model obtained after the integration of fermionic degrees of freedom.
The path-integral formalism elaborated here permits to calculate the correlation functions in the system and, as a result, we obtain the expressions of normal (both incoherent and coherent) and anomalous excitonic spectral functions, and shapes of density of states (DOS) corresponding. A special attention is payed, when calculating the bosonic phase-stiffness DOS function, which is negative. Furthermore, it is crucial for calculation of the coherent normal and anomalous excitonic DOS functions. Namely, considering the expressions of the bosonic and fermionic DOS functions, we calculate the total, phase coherent DOS functions, as convolutions from bosonic and fermionic counterparts.
For the incoherent partial normal fermionic DOS functions at T = 0 we obtain a hybridization-gap in the excitation spectra as a direct consequence of the presence of the Hartree-type gap in the single-particle energy scales. For the case of coherent normal and anomalous excitonic DOS functions, this hybridization-gap is absent totally for all frequency modes and for all values of the Coulomb interaction. This is due of strong coherence effects in the strongly correlated fermion system at low temperatures and at low densities. For the anomalous excitonic DOS function, we found that the hybridization-gap is absent for the case of small and intermediate values of the Coulomb interaction parameter, but there is a finite constant small gap that is opening in the strong coupling regime, signaling the passage to the SC (BEC) side of the SC-SM phase transition (the BCS-BEC crossover).
The paper is organized as follows: The paper is organized as follows:in Section II we introduce the model Hamiltonian and we discuss the main calculation schemes. In Section III, we describe the EI state, with a short discussion about important energy scales in the system. In Section IV we present the analytical calculation of spectral functions and DOS functions of the system in consideration. At the end of the Section IV we show the results of numerical evaluations for calculated DOS functions.

A. The EFKM Hamiltonian
As the model for study the excitonic condensation at low temperatures we have chosen the two-band extended Falicov-Kimball model (EFKM), 16,[30][31][32][33] due to its large applicability for treatment of the electronic correlations. 38,39 The Hamiltonian of the EFKM model is Heref (r) (c(r)) is the f (c) electron creation operator at the lattice position r, the summation r, r ′ is over nearest neighbors (n.n) sites on the 3D cubic lattice. The short hand notations are introduced n(r) = n c (r)+ n f (r) andñ(r) = n c (r) − n f (r) in order to simplify the calculations. Next, t c is the hopping amplitude for c-electrons and ǫ c is the corresponding on-site energy level parameter. Similarly, t f is the hopping amplitude for f -electrons and ǫ f is the on-site energy level parameter for f -orbital.
For t c t f < 0 (t c t f > 0) we have a direct (indirect) band gap semiconductor. The on-site (local) Coulomb interaction U in the last term of the Hamiltonian in Eq.(1) plays the coupling role between two bands. As we will see later on, the strength of the local Coulomb interaction will tune the semi-metal (SM)-semiconductor (SC) transition in the system. The chemical potentialμ isμ = µ−ǭ, whereǭ = (ǫ c + ǫ f ) /2. Naturally, we adjust the chemical potentials µ f and µ c in order to maintain the number of electrons in f and c orbitals separately. Then, the equilibrium value of the chemical potential µ ≡ µ f = µ c in Eq.(1) will be determined from the half-filling condition, i.e., we suppose that n c (r) + n f (r) = 1. We will use t c = 1 as the unit of energy and we fix the band parameter values ǫ c = 0 and ǫ f = −1. The Fermi energy level is assumed to be situated at the level of the c -band, thus ǫ F = 0. For the f -band hopping amplitude t f we consider the values t f = −0.3 and t f = −0.1 corresponding to the heavy hole and light hole f -bands.
Throughout the paper, we set k B = 1, = 1 and lattice constant a = 1. Also, we keep the frequency symbol ν for fermions and ω -for bosons, throughout the paper. In the case of degenerated f and c -bands, i.e., when ǫ f = ǫ c and t f = t c , the EFKM model reduces to the standard Hubbard model. 36 The principal advantage of the EFKM, in comparison with the genuine Falicov-Kimball model (FKM), 37 is that it is taking into account the direct nearest-neighbors f electron hoppings 38,39 (t f ) and it can be shown 36 that the EI state is unstable when the pure FKM is approached. In fact, the EFKM Hamiltonian in Eq.(1) is equivalent to the asymmetric Hubbard model, if we associate to the orbitals c and f the spin variables, by replacing the fermionic Hilbert space with the pseudo-fermionic one and then linearizing the interaction term via the bosonic states (see in Ref. 17).

The partition function
The Hamiltonian in Eq.(1) is containing two separate quadratic terms and is suitable for decoupling by functional path integration method. 41,42 We employ the imaginary-time, fermionic path-integral method and we introduce the fermionic Grassmann variables 41 f (rτ ) and c(rτ ) at each site r and each time τ varying in the interval 0 ≤ τ ≤ β, where β = 1/T with T being the thermodynamic temperature. 43 The time-dependent variables f (rτ ) and c(rτ ) are satisfying the anti-periodic boundary conditions The grand canonical partition function of system of fermions written as a functional integral over the Grassmann fields is where the action in the exponent is given in the pathintegral formulation in the form Here, S B [f , f ] and S B [c, c] are fermionic Berry terms. They are defined as follows andẋ(rτ ) = ∂ τ x(rτ ) is the time-derivative .
2. Decoupling of term proportional to n 2 We will decouple the quadratic density terms in Eq.(1) using the Hubbard-Stratonovich linearisation procedure 42 and by introducing new variables V (rτ ) and ̺(rτ ), conjugated respectively to the density terms n(rτ ) andñ(rτ ). For quadratic term, proportional to n 2 (rτ ) in the exponent of the partition function in Eq.(2), we have Furthermore, we combine the exponent in Eq.(5) with the effective chemical potential term, linear in total electron density n(r) (see the second term, in Eq.(1)). Then, we decompose the variables V (rτ ) into a static and the periodic part where, from time-periodicity ofṼ (rτ ), it follows the relation β 0 dτṼ (rτ ) = 0. As a result, the integration over variables V (rτ ) becomes the integration over the scalar static variables V 0 (r) and the integration over the periodic fieldṼ (rτ ): For the periodic part in Eq.(6), we introduce U(1) phase field variables φ(rτ ) by using a Faraday-type relation 44 Thus, for the dynamic part, we transform the integration over the gauge variablesṼ (rτ ), into an integration over generic phase variables φ(rτ ) The periodicity ofṼ (rτ ) implies that φ (rβ) = φ (r0). The integration measure, in Eq.(9), over the variables φ is defined as where the notations φ i and φ f mean the initial and final paths. The path-integral in Eq.(10) could be transformed into path integration over the compact U(1) group-manifold, since the electromagnetic group U(1), governing the phase field, is compact, i.e., φ(rτ ) has the topology of a circle (S 1 ), thus we have a non-homotopic mapping of the configuration space onto the U(1) gauge group S 1 → U (1). The paths, which loop around a circle in different number of times, are in different homotopy classes and they cannot be continuously deformed into one another. All these paths can be characterized by their proper winding numbers 45 m (r). Any two paths which have different winding numbers, cannot be continuously transformed one to another, and in order to include all possible phase path contributions, we have to sum over all topologically inequivalent phase configurations described by their winding numbers. 45 Accordingly, the path-integral in Eq.(10) is transformed as where the integration measure is now In performing the integration over the phase field one should take into account that the field configurations satisfy the boundary conditions 45 Thus, integration over all phases φ(rτ ) amounts the integration over the β-periodic field ϕ(rτ ) and the summation over a set of U(1) integer winding numbers m(r).
For the scalar static part V 0 (r), we get the following functional integral The saddle-point value of V 0 (r) is given by where n is total average particle density n = n c + n f (furthermore, we will fix n as equal to 1, corresponding to the case of half-filling). Thereby, after decoupling of the quadratic term proportional to n 2 in the Hamiltonian in Eq.(1), we get a contribution to the partition function in Eq.(2) in the form where the emergent phase-only action S[ϕ] is given as and the effective chemical potential µ n , attached to the total density operator in Eq. (16), is given in the form µ n = Un 2 −μ.
3. Decoupling of term proportional toñ 2 The decoupling of the quadratic term proportional tõ n 2 (rτ ) in the exponent of the partition function in Eq. (2) is also straightforward. We obtain After combining the expression in the exponent in Eq. (18) with the similar linear term in the expression of the Hamiltonian in Eq.(1) (see the forth-term in Eq.(1)), we have The saddle-point evaluation for ̺ gives whereñ = ñ(rτ ) is the average of the particle density difference function. As a result of the decoupling, we obtain a "Zeeman"-like contribution to the partition function with the attached effective chemical potential µñ =

Linearized action with phase-field contribution
To summarize, the grand canonical partition function of the system, after of both procedures of decoupling, is where the action S[c, c,f , f, ϕ] in the exponent is given by After the Hubbard-Stratanovich linearisation, we got the total action of the system that is linear in terms of fermion densities and contains in addition a phasedependent term S [ϕ], and also the terms, proportional to the effective chemical potentials µ n and µñ.

B. The U(1) transformation
In the perspective to treat the local and non-local correlations in our excitonic system, it is important to separate the U(1) gauge degrees of freedom related to the phase sector. To this end, we perform the local gauge transformation to new fermion Grassmann variablesf (rτ ) andc(rτ ). Meanwhile, this procedure will automatically eliminates also the last imaginary term, appearing in the expression of the phase action in Eq. (17).
For the electrons of f and c-orbitals, the U(1) transformation is whereÛ(ϕ) is the U(1) transformation matrixÛ(ϕ) = I · cos ϕ(rτ ) + iσ z · sin ϕ(rτ ) with the unit matrixÎ and σ z being the Pauli matrix. The variablesx =f ,c, and we used the bosonic phase variables ϕ introduced in Eqs. (11) and (12). After transformations given in Eq. (24), we obtain the total action of the system in the U(1) gauge invariant form (for comparison, see the action in Eq.(23) before transformations) Now,t and t in Eq.(25) are respectivelyf -band andc -band fermion transfer integrals. We got in Eq.(25) also, a new, emergent, quadratic phase action S 0 [ϕ] The partition function of the system in the new variables f andc is From this form of the partition function we will generate the effective actions for fermions and for bosonic phase sector (see the general procedure presented in Fig. 1).

III. EFFECTIVE ACTION FOR FERMIONS
By following left-lowest root, presented in Fig. ??, we will integrate out the phase variables. We obtain where the effective, phase-averaged fermionic action in the exponent is The Fourier transformation of fermionic variablesf and c is where N , is the number of lattice sites, and ν n = π(2n + 1)/β are the Fermi-Matsubara frequencies 43 with n = 0, ±1, ±2, ... . Then, the effective phase-averaged fermionic action of the system in the Fourier space takes the following form Here, G −1 (k, ν n ) is the inverse of the Green function matrix where we have introduced the local excitonic order parameter ∆ = c (rτ )f (rτ ) . The single-particle Bogoliubov quasienergies Ef k (ν n ) and Ec k (ν n ) are Ef k (ν n ) = ǫf − iν n −t k and Ec k (ν n ) =ǭc − iν n − t k . Next,t k and t k are band-renormalized hopping amplitudest k = 2tg B γ k and t k = 2tg B γ k , where g B is the 3D bandwidth renormalization factor and , being the components of the lattice spacing vector d = r − r ′ with r and r ′ , being nearest neighbors site positions. For the simple cubic lattice we have d α ≡ a = 1.
The quasiparticle energiesǭf andǭc are of Hartreetype and they are defined in the theory by relation ǫx = ǫ x − µ + U nỹ + i φ(rτ ) , whereỹ means orbital, opposite tox. The EI low-temperature phase is characterized by the local excitonic order parameter ∆. Without any loss of generality, we can suppose the case of the EI state, with uniform real gap parameter ∆ =∆. The EI state develops from local on-site electron-hole correlations. The expectation value, given in the expression of local EI order parameter, could be calculated in the frame of path-integral method, as well as the fermion density averages of the respective band levels nx = x (rτ )x(rτ ) . We get a set of the coupled self-consistent equations for the EI order parameter ∆, single-particle fermion densities nx, and EI chemical potential µ Here f (ǫ) = 1/ e βǫ + 1 is the Fermi -Dirac distribution function, ξ k = −t k +ǭc +t k −ǭf is the quasiparticle dispersion and the energy parameters E + k and E − k are defined as In fact, the difference between Eqs.(34) -(36) and the Hartree-Fock results, 20 is in the presence of bandwidth renormalization factor g B . The calculation of the factor g B (r − r ′ ) could be done effectively within the self-consistent-harmonic-approximation (SCHA) method. [47][48][49] In this approximation the quantum rotor description is reduced to classical Hamiltonian one by the Feynmann-Kleinert minimization procedure. 47,48 Our results for SCHA show that the factor g B is equal identically to 1 at T = 0 as in the twodimensional (2D) case (see in Ref. 47 for details). For higher temperatures, it differs from unity, but not much. The difference between the energy parameters in Eq.(37) defines the charge transfer gap in the system The numerical solution of the system of equation  [15][16][17][18][19][20]31,32 This good correspondence is ascribed to a rather weak band renormalization at T = 0 (in our case g B = 1). At finite temperatures, the particle number fluctuations are important and the band renormalization becomes necessary, especially, when approaching from the band insulator (BI) high-temperature side. Indeed, as the numerical evaluations show, the transition temperature T EI of the e-h pair formation is not vanishing for the case of the vanishing narrow band hoppingt = 0.
The exact solutions for the chemical potential could be obtained from Eqs.(34) - (36), both at the boundary of EI transition (i.e., when ∆ (T EI , U ) = 0) and in the EI state region (∆ (T < T EI , U ) = 0). The results are discussed in Ref. (51).
Meanwhile, it is also shown in Ref. (51) that the excitonic BEC transition critical temperature T c is much smaller than the critical temperature T EI of the excitonic pair formation, in good agreement with previous theoretical predictions. 15,[23][24][25]47,52 The self-consistent numerical solutions for T EI and T c are shown in Figs. 2 fort = −0.3. It is clear that, for the case of intermediate and strong interaction limits, the coherent BEC transition critical temperature T c is smaller of about two orders of magnitude than temperature T EI . Contrary, in the very small interaction case, we have the coincidence of both transition temperatures. This result is expected also from general considerations, because at small U/t we are in the BCS limit, which means that the pairing and condensation occur simultaneously.
A similar effect is found recently in Ref. 47, considering 2D excitonic systems, where the exciton-superfluid transition critical temperature is found much smaller than the critical temperature T EI of exciton pair formation. Such a reduction of the transition temperature, due to the coherent pairing scheme, is given also in Ref. 52, where the BCS-Bose crossover is studied in the 2D attractive Hub- bard model. Especially, a HF pair formation temperature is estimated, which corresponds to the regime of the incoherent or local pairs (for a comparison see the correspondence to our T EI , and order parameter ∆, which is also local), in difference with superconducting pairing temperature, at which the coherent Cooper pairs start to be formed. The general idea used in Ref. (51) is based on the non local n.n. excitonic exchange correlation mechanism and Bogoliubov's mean-field self-consistency assumption 53 for the effective phase action obtained after fermion integration procedure in the action in Eq. (27) (this is represented by the right-lowest root in Fig. 1, given in the Section II B). The obtained Hamiltonian for the phase sector is very similar to the classical Hamiltonian one, [47][48][49]51 with an effective phase-stiffness parameter J that emerges after the fermion Wick averaging procedure. 43 The corresponding phase action is 51 (38) It appears that the non-zero value of quantity J is directly related with the pairing gap ∆, since it is shown 51 that J vanishes, when ∆ = 0. Here, we present only the final analytical result for J (for details see in Ref. (51)) where z is the number of n.n. sites on the 3D lattice. The parameters Λ 1 (k, k ′ ) and Λ 2 (k, k ′ ) entering in Eq.(38) are given by Also, J is strictly positive for all the regions of the normalized Coulomb interaction parameter U/t. In addition, it follows from the analytical form of J, that the macroscopic phase coherence in the system is characterized by an energy scale J ∼ t e · t h /(t e + t h ) for all values of the Coulomb interaction parameter U . This is related to the motion of the center of mass of e-h composed quasiparticle, 15 because t e · t h /(t e + t h ) ∼ (m e + m h ) −1 . For the strong interaction case, we are converging with the hard core Boson model, with the kinetic energy scale ∆t e · t h /U (with ∆ being the local excitonic order parameter). Thereby, it is shown in Ref. (51) that non-local correlations between the electrons and holes of different n.n. excitonic pairs, are related with the excitonic BEC condensation. Furthermore, in the frame of the quantum rotor model, the excitonic BEC transition probability function is derived and its temperature dependence is found. 51 As a systematic continuation of theoretical study given in Ref. (51), we elaborate here on the analytical forms of normal, f and c -band, Green functions (incoherent and coherent) and coherent, anomalous excitonic Green function. Also, we will discuss in details, the spectral functions and density of states (DOS) corresponding, in the next. The numerical evaluations of calculated DOS functions are shown in the Section IV D.

A. Normal and excitonic DOS functions
The spectral density functions of the system of interacting exciton gas will determine the excitonic centerof-mass distribution related to the condensation in the low-temperature limit. Therefore the calculation of these functions represents an important task. Within our theoretical approach we can access a variety of correlation functions in the system. We will concentrate now on the c and f -band normal spectral functions and also the excitonic anomalous spectral function, represented in terms of the initial operators. The c and f -band normal excitonic Green functions are the anomalous excitonic Green function is defined as After introducing the U(1) transformations, defined in Eq.(24), we will have the Green function's decomposition The phase factors, appearing in the definitions of Green functions, define in fact a charge-bosonic propagator G z (rτ, r ′ τ ′ ) as follows Then, we pass to the Fourier space representation for the U(1) -transformed Green functions and Furthermore, the Fourier transformation of the functions in Eqs.(42)-(43) will be rewritten as a convolution in the reciprocal k-space and It is worth to mention that frequency summations in Eqs. (48) and (49) are over Bose-Matsubara frequencies ω n = 2πn/β. The Fermionic Green functions will be calculated using the formalism discussed in the Sections II and III and also, functional derivation techniques. 41 Particularly, for the f and c -band Green functions Gx ,x (k, ν n ), we get In general, the experimental observation of hybridization between the valence band and conduction band could be done by examining the ARPES spectra, which measure the spectral intensities just above and below the temperature T EI of excitonic pair formation. In ARPES experiments, one observes the imaginary part of the real retarded Green function, therefore the calculation of it represents a remarkable importance. Indeed, the singleparticle density of states is related with the imaginary part of the retarded Green functions, thus we need to calculate real retarded function, which corresponds to the normal Matsubara Green functionGx ,x (k, iν n ). This could be done by the analytical continuation into the upper-half complex semi-plane (ν n > 0) of frequency modes iν ñ The single-particle DOS is defined then as where and k-summed DOS will be The summations over the wave vectors in Eq.(55) can be simplified by introducing the appropriate DOS function for the 3D cubic lattice ρ 3D (x) = 1 N k δ(x − γ k ). Then, it is not difficult to show that Here,t(x) = 2tx and t(x) = 2tx, and the energy parameters E ± (x), are continuous versions of parameters defined in Eq. (37). The density of state ρ 3D (x), for the simple cubic 3D lattice, is given by where Θ(x) is the Heaviside step function, and K(x) is the elliptic function of the first kind: K(x) = π/2 0 dt1/ 1 − x 2 sin 2 t. For the anomalous excitonic Green function, we havẽ The retarded function, which corresponds to it, is then 43 The single-particle DOS is forthcoming then as Obviously, it has more simple form than the function in Eq. (52). The k-summed DOS for excitons will be where the dimensionless parameters Λ 1,2 (ν) are given by following expressions and the functions χ i [Λ 1 (ν)] (i = 1, 2), in the denominators in the right-hand side in Eq.(61) are and Now, turning to the convolution forms for total fermionic and excitonic Green functions in Eqs. (48) and (49), we need an explicit expression for the phase-bosonic Green function G z (k, ω n ). We will calculate it in the formalism of the effective phase action given in the quantum rotor model, discussed earlier in Ref. 51, where we have derived the effective phase-only action S eff [ϕ] by integrating the fermions in the partition function in Eq. (27). In the following, we cast S eff [ϕ] into the quantum rotor representation. 51

B. The phase-stiffness DOS
To proceed, we replace the phase degrees of freedom with the complex, unimodular field z(rτ ) = e iϕ(rτ ) , which satisfies the time-periodic boundary condition z(rβ) = z(r0). The spherical constraint, imposed on a set of the unimodular variables is 1/N k |z(rτ )| 2 = 1. Now, we introduce these new variables z(rτ ) into the partition function in Eq. (27), in a way, consistent with the Faddeev-Popov ghost-field method 54 The phase-phase propagator G z (rτ, r ′ τ ′ ) will be rewritten in terms of z(rτ )-variables as follows Furthermore, the variables z(rτ ) play the role of the phase-flux attached to the fermions (see discussions in Section II). In general case, the local expression of the phase-phase correlation function in Eq.(67) is equal to unity, but, at very low temperatures (especially at T = 0), this law breaks down, because we have to consider the symmetry breaking related to the bosonic sector, thus, critically, we have fluctuation form z(rτ ) = e iϕ (rτ ) +z(rτ ), and the unimodularity constraint for z-field is violated.
Indeed, in the very low temperature limit, considering the BEC of excitons, we have the spontaneous breaking of local U(1) gauge-symmetry, related to the phase field, leading to the non-vanishing expectation value of the e iϕ (rτ ) . In order to demonstrate this, we separate the single particle states k = 0 via the Bogoliubov displacement operation (see for details in Refs.1, 51 and 55). This is so-called Bogoliubov phase coherence mechanism discussed in details in Ref. 1. Then, we write for the complex variables z(k, ω n ) where ψ 0 is the BEC transition amplitude ψ 0 = z(k, ω n ) . Next,z(k, ω n ) are the excitation part 58 (oncondensate) of effective bose-field. The Fourier transformation of the phase-phase propagator G z (rτ, r ′ τ ′ ) in Eq. (47) is We consider the expectation value z(k, ω n )z(k, ω n ) in the local limit, i.e., when d = r − r ′ = 0 and τ − τ ′ = 0 and we should draw the condensate part, by applying the transformation in Eq.(68). Hence, we have Thereby, in Eq.(70) we have defined the coherent macroscopic state for the excitonic system in the low temperature limit, and an excitonic BEC is expected in the next. The Fourier Green function G z (k, ω n ), defined in Eq. (47), could be calculated within the quantum rotor model, and we give here only the final result (for details, see in Ref. 51).
where γ −1 (ω n ) is the inverse of the Fourier transformation of the two-point phase-phase correlation function. 51 We have where Z 0 is the partition function of the non-interacting Bose sector The summations, in Eqs. (72) and (73), run over topological winding numbers m of the group U(1). The retarded bosonic Green function 58 is related to the Matsubara Green function, by the analytical continuation And the k-summed DOS for bosons reads as After non difficult algebraic manipulations and replacing the summation in Eq.(74) by integration with the help of 3D density of states ρ 3D ( where κ i (x) i = 1, 2 are given by the following relations κ 1,2 (x) = −μ ± μ 2 + 4U J (3 − x) and the stiffness parameter J is given in Eq. (38) in Section III. As it could be expected the Bosonic DOS function Eq.(76) is negative ρ z (ω) < 0. This is consistent with the general considerations of the weakly non-ideal Bose gas. 58

C. Spectral density functions and fermionic DOS
Furthermore, we separate the condensate modes {q = 0, ω n = 0} in Eqs. (48) and (49). We have and As we see, the normal and excitonic propagators are composed of two parts, one, responsible for the condensate state and the other-the on-condensate excitation part (see discussion in Ref. 58 for the case of the pure Bosegas). Note also, that first terms in the right-hand sides in Eqs.(77) and (78) consist of the condensate-transition probability function |ψ 0 | 2 , multiplied with the fermionic propagators Gx ,x (k, ν n ) and Gc ,f (k, ν n ). Now, we are ready to calculate the analytical forms of the normal excitonic spectral functions A x,x (k, ν) (x = f, c) and anomalous excitonic spectral function A c,f (k, ν) and, later on, the profiles of the respective DOS, including states of the condensate. We introduce here the spectral functions A x,x (k, ν) and A c,f (k, ν) that carries the same physical information as the correlation functions G x,x (k, ν n ) and G c,f (k, ν n ). We have The integration here, is over continuous frequencies.
Note, that G x,x (k, ν n ) and G c,f (k, ν n ) are total fermionic Green functions, including also the convolutions with bosonic parts. In the same way, we can introduce the spectral functions A z (k, ν), Ax ,x (k, ν) and Ac ,f (k, ν), associated with the charge and the pure fermionic parts (without bosonic sector). They correspond respectively to the correlation functions G z (k, ω n ), Gx ,x (k, ν n ) and Gc ,f (k, ν n ). We have the following equations for these counterparts and Using these definitions we get for the total spectral density functions A x,x (k, ν) and A c,f (k, ν) see Appendix ? ? and where n(ǫ) = 1/ e βǫ − 1 is the Bose-Einstein distribution function. The proof of these relations is given in the Appendix ??. From the spectral functions we can obtain the corresponding DOS, by summing over the reciprocal wave vectors k, hence, the total DOS are ρ x,x (ν) = 1 N k A x,x (k, ν) and ρ c,f (ν) = 1 N k A c,f (k, ν). Furthermore, using the expressions for A x,x (k, ν) and A c,f (k, ν) in Eqs.(84) and (85), we get for the total DOS functions and whereρx ,x (ν) andρc ,f (ν) are DOS, corresponding to the excitation part of the system and are given, as convolutions, in terms of continuous frequency modes (see Ap- A key feature of the results in Eqs.(86) and (87), is that we have separated DOS contributions coming from the condensate and excitation parts. We define also the total DOS function as In the Section IV D, we present the numerical evaluations of all discussed DOS functions.

D. Total DOS functions
Employing Eqs. (56), (61) and (76), we can obtain the explicite analytical expressions for DOS functions in Eqs. (86) and (87). For the normal f and c -band excitonic DOS functions, we obtain For the anomalous excitonic DOS function, we have The numerical evaluations of calculated DOS functions at T = 0 are given in Figs. 3 -7 fort = −0.3. The presence of singularities in the integration region causes that we used an adaptive 21-point integration routine combined with the Wynn ǫ-algorithm, 59 to calculate those integrals numerically. The accuracy for adaptive evaluations is achieved with an absolute error of order of 10 −4 and with a relative error of order of 10 −7 . Particularly, in Figs. 3 and Fig. 4, we have presented purely fermionic normal single-particle (incoherent) DOS ρx ,x (ω). We examine their behavior over the entire BCS-BEC crossover region (i.e., for different values of the Coulomb interaction U ). An artificial Lorentzian broadening η = 0.01 is used in numerical evaluations for the incoherent partial and total DOS functions of f and corbitals. The chemical potential values are inputed along the upper-bound µ max , where they are maximal (see in Ref. 51). The principal reason of it is that the BEC transition amplitude ψ 0 has no physical solutions along the lower-bound µ min of the chemical potential. On the other hand, the values µ max are most convenient, because they are minimalizing the Hamiltonian of the system. Furthermore, as the reference for the coherent BEC transition amplitude |ψ 0 | 2 , we considered the self-consistent calculation-results from the work in Ref. 51, where this function is calculated both analytically and numerically for different values of the Coulomb interaction parameter U and for different temperatures, including the zero temperature limit. We use here the numerical data, which are evaluated there.
We see in Fig. 3, that at the small interaction limit, when 2 U 6 (i.e., the BCS limit), the incoherent excitonic DOS exhibits a BCS-like double-peak structure (see the panels I -III in Fig. 3) and the peaks are separated with a well defined hybridization-gap. The principal reason of it is the non-vanishing Hartree-gap ∆ H = 0 in the single-particle energy-spectrum discussed above, in the Section III. The hybridization-gap is proportional to the parameter U and it is increasing with the increase of parameter U . We observe also that the peaks becomes more separated when increasing of U . In the strong interaction limit this displacement is stabilizing and we have practically constant value of the hybridization-gap, when further increasing the interaction (see the panels I -III, in Fig. 4). The results in Figs. 3 and 4 are very similar with the previous theoretical results. 15,20,56,57 Especially, they are close to those presented in Refs. 17 and 20, where the partial incoherent f and c -band normal DOS functions and total DOS is calculated using HF and SO(2)-invariant slave boson approaches.
In Fig. 5 (see the panels a -c), we have shown the coherent condensate part of the anomalous excitonic DOS given by the first term in the right-hand side in Eq.(A.4), corresponding to the fundamental state (k = 0). Different values of the Coulomb interaction are considered, including the small and strong interaction cases. Here, we observe again the double-peak fermionic structure, but there is no the hybridization-gap for the small and medium values of the Coulomb interaction parameter (see the panel-a and panel-b in Fig. 5), and we have finite number of states for all values of the frequency modes ν. This is due to the coherence effects and the presence of the coherent excitonic condensate at the fundamental mode ν = 0. For higher values of U in Fig. 5 (see the panel-c), this double-peak structure in DOS is smoothing, but, in contrast to the incoherent normal DOS behavior (see in Figs. 3 and 4), here a small Mott-gap appears at the very high values of the Coulomb interaction parameter (see the DOS curves for U = 8 and U = 9.6 in the panel-c in Fig. 5). This is due to the fact, that the very strong Coulomb interaction has a destructive role on the condensate state and, in the large-U limit of interaction, we have the destruction of the excitonic condensate, and a very small, Mott-type hybridization-gap is enhanced. In this region of the interaction, we have the coherent exciton DOS separation into two separate parts, similar to the case of the incoherent DOS functions in Figs. 3 and 4. In the positive (negative) frequency regions, the DOS spectrum in Figs. 3, 4, and 5 is slightly displacing and broadening into the direction of higher (smaller) frequencies, for both -normal (incoherent) single-particle fermionic, and coherent excitonic DOS functions and, for both, we observe also a gradual decreases in the DOS amplitudes, across the whole BCS-BEC crossover region, when increasing the Coulomb interaction parameter U . In Fig. 6, we have shown the total coherent fermionic DOS functions ρ(ν), given by the Eq.(90), for different values of the Coulomb interaction parameter U and in the limit of zero temperature. An artificial Lorentzian broadening η = 0.01 is again used during the numerical calculations. Contrary the case of the incoherent DOS in Figs. 3 and 4, the coherent total fermionic DOS in Fig. 6 shows a different, gapless behavior. Here, again, as in the case of the incoherent DOS functions, for small and intermediate interactions U , we have typical doublepeak fermionic structure in the DOS, which is smoothing and disappearing totally in the strong interaction limite. But, in difference with incoherent DOS functions, here we have not the presence of the hybridization-gap in the spectra and we have always have a finite number of states for all values of the frequency modes (see in Fig. 6). The reason of this gapless DOS behavior could be related to the strong coherence effects between two bands, which is due to the presence of the phase stiffness mechanism considered here. It is worth to mention that    another gapless-type behavior in the DOS spectrum is found recently in Ref. 60, where this effect is associated with metallic charge-density-wave phase and is driven by strong electron correlations. In Fig. 6, all figures are combined together, in the way, to see the total DOS evolution with variation of the interaction parameter U . We see clearly, how the tuning of the interaction parameter U , affects the general DOS behavior in the model, by reducing the DOS amplitudes with increasing of U and reducing also the number of states at the Fermi level ρ(ν = ǫ F ) when increasing U . Thereby, in the case of normal fermionic DOS (both incoherent and coherent), presented in Figs. 3, 4, 5 and Fig. 6, the double-peak structure disappears for U 6, signaling the appearance of the SM (BEC) limit of the transition. In Fig. 7 a comparison is given for the incoherent and coherent c -band normal DOS behaviors. Three different values of the Coulomb interaction U are chosen (see the panels -a, b and c). We see clearly in Fig. 7, how the coherence effects are reducing the incoherent DOS amplitudes and DOS spectra become broader, along frequency axis, and also, there is no gap near the Fermi level (ρ(ν ∼ ǫ F ) = 0).
In Figs. 8 and 9, we have presented the temperature dependence of the c -band normal fermionic (Fig. 8) and anomalous excitonic DOS functions (Figs. 9). They are given by the first terms in the right-hand side in Eqs.(A.4) and (92). In the left-panel in Fig. 8, we have presented the temperature dependence of the singleparticle DOS |ψ 0 | 2 ρ c,c (ν), which corresponds to the fundamental state k = 0, and for U = 6. The caset = −0.3 is considered. The corresponding values of the amplitude ψ 0 of BEC transition, are taken again from the work in Ref. 51. In the right-panel in Fig. 9, the same function is plotted for the case U = 9. As we see in Fig. 8, the hybridization-gap is still open for all values of the temperature. So the system is always an insulator.
In Fig. 9, the temperature dependence of the phasecoherent anomalous excitonic condensate part of the DOS is presented (i.e., k = 0, and without excitation part) for U = 2 and U = 6. We observe in all Figs. 8 and 9, that the temperature has a destructive effect on the DOS amplitudes.
We realize also, from Figs. 8 and 9, that the anomalous excitonic DOS functions vanish at the temperatures, that are far away from the region of the EI transition of about two orders of magnitude (compare the temperature scales in Figs. 8 and 9 with those given in Fig. 2 in the Section III). This result, could be regarded as a good proof of the theory elaborated in Ref. 51, and we can theoretically clearly state that the excitonic BEC and EI states are not the same phases of matter. We see also, in Figs. 8 and 9, that the normal single-particle DOS persists for a rather large values of temperature, than the anomalous condensate DOS, because of the presence of the hybridization-gap in the low-energy spectra.

V. FINAL REMARKS AND CONCLUSIONS
We have studied 3D system of conduction band electrons and valence band holes in the frame of the extended Falicov-Kimball model. We have implemented the path-integral formalism, in which the Coulomb interaction term is expressed in terms of U(1) quantum phase variables ϕ conjugated to the local particle number, providing a useful interpretation of the problem. In Section III, we have shown that at low temperatures, the electron-hole system becomes unstable with respect to the formation of the excitons at T = T EI , and the local gap ∆ is present in the excitation spectrum, controlled by the Coulomb interaction parameter U/t, which gives the relevant energy scales for the excitonic insulator state. Here, as a result of the spontaneous symmetry breaking, an expectation value of e iϕ = 0 appears, which is signaling of the presence of the phase coherence in the system. Furthermore, pairing and condensation are not generally the same, as it was admitted in the literature, except the weak interaction limit, when we have a BCS-like condensation of excitonic pairs. However, in the excitonic system with the strong pairing, we have the situation, where the pairs are strongly bound, but are uncorrelated one with each other, until they become phase coherent at temperatures T T c . This situation was studied in details in Ref. 51.
We have evaluated the normal and anomalous excitonic spectral functions for the f and c -band. We have determined, both analytically and numerically, density of states (DOS) spectra, governed by the pure-fermionic   part (due to the condensate modes k = 0), and the excitation spectra. We have shown that there is a usual hybridization-gap in the normal (incoherent) f and cband DOS structures. Contrary, in the case of the coherent normal DOS functions (presented in Fig. 6) this gap is lacking, and there is always a finite number of states at all frequency modes ν. We associate this result to the strong coherence effects, that are present in the system at low temperatures. In the anomalous excitonic DOS structure we have found that the hybridization-gap is absent for the weak and intermediate values of the Coulomb interaction parameter U and this is due to the presence of the coherent excitonic condensate and strong coherence effects. A very small gap is opening in the spectra in the strong interaction limit, signaling the destruction of the coherence and condensate state.
The excitonic phase coherence may be evidenced by the coherence of their light emission, which can be studied by interferometry measurements. 61 Therefore, measurements of the intensity of the line-shape of the excitons decay (by emitting photons, upon electron-hole recombi-nation) may be a powerful probe of the DOS spectra in the excitonic systems. 62 However, a final remark that should be featured is such, that for the experimental determination of the excitonic Bose-Einstein condensate, the ARPES measurements should be provided at temperatures much lower than temperatures at which the excitonic insulator state is being determined. 20 This is important for the achievement of macroscopic phase coherece between excitonic pairs, and the strong coherence effects are manifesting at the very low temperatures.
As a continuation of our studies, we would like to consider the role of the charge-density-wave-like excitations in the excitonic systems, within the EFKM model and consider the temperature effects. It is especially interesting to find out how the coherent excitonic DOS will be affected in the case of presence of such elementary excitations. This will give us a more complete picture on the excitonic phase transition scenario.  77) and (78). Then, it is easy to see that the following identity holds for the integrals of the normal x = f, c -band fermionic spectral functions.
dν A x,x (k, ν) iν n − ν = |ψ 0 | 2 · dν Ax ,x (k, ν) iν n − ν + Furthermore, for calculating the Matsubara sums over bosonic frequencies ω n , we will use the property of the Bose-Einstein distribution function n (ǫ) (see the Section IV C) of complex argument where f (ǫ) is the Fermi-Dirac distribution function. In Eq.(A.2), we used the fact that ν n are even Fermion-Matsubara frequencies ν n = π β (2n + 1). Next, we sum the Bosonic Matsubara frequencies in Eq.(A.1) amd we rewrite the equality in Eq.(A.1) in the following form The derivation of the convolution form, given in Eq.(85), is exactly the same. Summing over the wave vectors k in Eq.(A.3), we will have the forms of the normal f and c -bands density of states ρ x,x (ν) and also the excitonic DOS function ρ c,f (ν) given in Eqs.(86) and (87) in the Section IV C.

Normalization conditions for DOS
We have a composed nature of the interacting electron, given in the Section II, thus, we have various sum rules, corresponding to different counterparts of the total coherent DOS functions. As a general rule, the partial normal DOS, and anomalous excitonic DOS functions ρ x,x (ν) and ρ c,f (ν) satisfy following normalization conditions