Fermi liquid-like behaviour of cuprates in the pseudogap phase simulated via T-dependent electron-boson spectral density

We investigated the temperature- and frequency-dependent optical scattering rates in the pseudogap phase of cuprates using model pseudogap and electron-boson spectral density (EBSD) functions. We obtained the scattering rates at various temperatures below and above a given pseudogap temperature using a generalized Allen’s (or Sharapov’s) formula, which has been used to analyse the measured optical spectra of correlated electron systems with a non-constant density of states at finite temperatures. The pseudogap and EBSD functions should be temperature dependent to simulate the Fermi liquid-like behaviour of underdoped cuprate systems observed in optical studies. Therefore, the observed Fermi liquid-like behaviour can be understood by considering the combined contribution from the T-dependent EBSD function and the T-dependent pseudogap. We also obtained the optical conductivity spectra from the optical scattering rates and analyzed them to investigate intriguing electronic properties. We expect that our results will aid in understanding the Fermi liquid-like optical response in the pseudogap phase and in revealing the microscopic pairing mechanism for superconductivity in cuprates.


I. INTRODUCTION
Copper oxide superconductors (or cuprates) have been studied intensively and most thoroughly since their discovery 1 .The temperature-doping phase diagram has become increasingly complicated 2 .Each characteristic phase region in the phase diagram is still not fully understood.The pseudogap 3 is the most studied phase among the phase regions.However, the origin of the pseudogap is not well understood.The origin of the pseudogap is crucial for understanding the superconductivity in cuprates.A study 4 claimed to have found spectroscopic evidence for Fermi liquid-like behaviour in the temperature-and frequency-dependent relaxation rate of heavily underdoped single crystals of HgBa 2 CuO 4+δ (Hg1201) with T c =67 K.This Fermi liquid-like behaviour occurs at low temperatures below the onset temperature of the pseudogap phase (T * ) and at low frequencies below the pseudogap energy (∆ pg ).
These findings can provide important information to better understand the metallic state and high-temperature superconductivity in cuprates.There are some recent intriguing optical studies on Fermi-liquid-like behavior and the strange (or bad) metal states of cuprates 5,6 .
The authors showed that the optical conductivity in the low-energy region consisted of two components, which are the Drude component and the mid-infrared absorption one.They analyzed the measured optical spectra to reveal temperature-and doping-dependencies of these components over a wide doping range, from underdoped to overdoped.
In our study, we mainly focus on simulating the measured optical spectra of a undedoped cuprate starting from a model electron-boson spectral density (EBSD) function using the generalized Allen's (or Sharapov's) formula 7,8 and the extended Drude model [9][10][11][12] .We used a well-established phenomenological approach 7,8,13,14 to simulate the Fermi liquid-like behaviour of cuprates in the pseudogap phase.This approach has been used to extract the electron-boson spectral density functions (or information on correlations) from measured optical spectra of cuprates 7,8,[13][14][15][16][17][18] .We started with a model electron-boson spectral density (EBSD) function 18 and a pseudogap model 19 for a given temperature, and obtained a corresponding optical scattering rate (or the imaginary part of the optical (or two-particle) self-energy) using a generalized Allen's (or Sharapov's) formula 8 , which can be used for analyzing a measured spectrum of a correlated electron system with a non-constant density of states, such as the pseudogap in cuprates 15 , and at a finite temperature.We then used the Kramers-Kronig relation between the imaginary and real parts of the optical energy 20 to obtain the corresponding real part of the optical self-energy.Furthermore, we obtained the optical conductivity from the complex optical self-energy using the extended Drude model formalism [9][10][11][12] .We used two different sets of the EBSD functions (one is temperature dependent, and the other temperature independent) and found that the experimentally observed optical scattering rate could be simulated with the temperature-dependent EBSD function.
The results showed that the temperature dependencies of the model EBSD function and model pseudogap are crucial in describing the observed optical response of cuprates in the pseudogap phase.Therefore, this can be a piece of important evidence that the EBSD function intrinsically contains temperature dependency, which may be critical in understanding superconductivity in cuprates.Our results will help us to understand the Fermi liquid-like optical response in the pseudogap phase and to find out the elusive microscopic pairing mechanism for superconductivity in cuprates.

II. THEORETICAL FORMALISM
Information regarding the correlations between charge carriers in a correlated electron system is encoded in the measured optical spectrum.In more detail, because the occupied and unoccupied electron states are involved in the optical spectrum, the band renormalization caused by the correlations will naturally appear in the measured optical spectrum, which contains intraband and interband transitions from filled states to empty states and is governed by the joint density of states and the dipole selection rule.The encoded correlation information can be extracted from the measured optical spectrum using a well-established approach [7][8][9][10][11][12][13][14] .The established approach consists of a couple of formalisms, such as the extended Drude formalism and the generalized Allen formalism.The correlated carriers can be described by the extended Drude formalism, whereas the free electrons (or Fermi gas) can be described by the Drude formalism.In the extended Drude formalism [9][10][11][12] , the complex optical conductivity (σ(ω)) can be written as follows: where Ω p is the plasma frequency of the itinerant electrons.Here the impurity elastic scattering rate in the Drude model is simply replaced with a frequency-dependent quantity to describe the energy transfer between the charge carriers caused by the correlations and the frequency-dependent quantity should be a complex function to maintain the causality.
The complex quantity ( Σop (ω, T )) is called the optical self-energy and carries the information of correlations between charge carriers.The real and imaginary parts of the optical self-energy form a Kramers-Kronig pair.In principle, the correlations may result in band renormalizations, and the band renormalizations can be encoded in the measured optical spectrum.Therefore, the information on correlations can be encoded in the complex optical self-energy, which is the deviation from the Drude band or, more generally, the bare band.
The bare band can be obtained using the local density approximation calculations.The optical self-energy is associated with a two-particle (or optical) process because both the filled and empty states (or electrons and holes) are involved in the optical (or absorption) process.Furthermore, the optical self-energy is closely related to the quasiparticle self-energy ( Σqp (ω)), which is associated with a single particle process because either only the filled states (or electrons) or the empty states (or holes) are involved in the process.Therefore, the optical self-energy is more complicated and can be related to the quasiparticle self-energy as, Σqp (ω) = d[ω Σop (ω)] dω 21 .In principle, the quasiparticle self-energy can be measured using angle-resolved photoemission spectroscopy (ARPES), whereas the optical self-energy can be measured using infrared/optical spectroscopy.We note that using ARPES, the quasiparticle self-energy of only the filled states (or electrons) can be obtained because of the experimental limitation; only the photoelectrons from the filled states can be measured with the technique.
A force-mediated boson exchanging model can be used to describe the correlations.In this model, the correlations can be described in terms of the EBSD function, I 2 B(Ω), where I is the coupling constant between an electron and a boson and B(Ω) is the boson spectrum.
In this case, the optical self-energy can be expressed in terms of the EBSD function.Allen originally developed an expression 7 , which can be used for the case of the normal state at T = 0 and the s-wave superconducting state.The Allen's expression has been generalized by Shulga et al. 13 for use in the case of a finite temperature, and further generalization has been achieved by Sharapov et al. for use in the case of a finite temperature and a non-constant density of states, such as the pseudogap phase in underdoped cuprates 8 .This generalized Allen formalism 8 can be written as follows: where n B (ω, T ) and f (ω, T ) are the Bose-Einstein and Fermi-Dirac distribution functions, respectively, and N (x, T ) is the normalized density of states, which is used to describe the pseudogap.The pseudogap (PG) can be modelled 22,23 as follows: where N 0 (T ) is the normalized density of states at the Fermi level and ∆ pg is the size of the pseudogap.This pseudogap model has been used for analyzing the optical spectra of underdoped cuprate systems 18,19,[22][23][24] .For this pseudogap model, the density of states loss in the pseudogap is fully recovered just above the pseudogap within 2∆ pg 19,22 .For Fermi surface arcs, the normalized density of states at the Fermi level is proportional to the temperature 22 , i.e.N 0 (T ) = T /T * , where T * is the pseudogap temperature.Because the measured optical spectrum is k-space averaged, a k-space averaged PG model is used, and the Fermi arc model is used for the T -dependent depth (1 − N 0 ) of the PG, which is closely associated with the d-wave-like PG and may also be related to the d-wave SC gap.
The EBSD functions, I 2 B(Ω, T ), of cuprates have been obtained using various spectroscopic experimental techniques, such as tunneling spectroscopy, ARPES, optical spectroscopy, Raman spectroscopy, and inelastic neutron scattering 25 .The obtained EBSD functions showed a unified temperature and doping-dependent evolution.A broad spectrum at high temperatures transforms into a peak between 30 and 60 meV and a featureless highenergy background through a spectral weight redistribution as the temperature decreases.
The overall spectral weight shifts to a higher energy region and the spectral weight redistribution weaken as the doping increases. 26,27.In this study, we used the temperature-dependent model I 2 B(Ω, T ).The model I 2 B(Ω, T ) consists of two (sharp and broad) components.For the broad component, we used an MMP mode 28 , which was introduced by Millis, Monien, and Pines (MMP) and was used to describe the antiferromagnetic fluctuations.For the sharp component, we used a sharp mode, which has been previously used to analyze the measured optical spectra of underdoped cuprates 18 .The model I 2 B(Ω) can be written as follows: where the first term is the broad MMP mode and the second one is the sharp mode.A b and ω b are the area under and peak frequency (or energy) of the broad MMP mode, respectively.
A s and ω s / 4 √ 3 are the area under and peak frequency of the sharp mode, respectively, and ω c is the cutoff frequency.The area under each mode and the peak frequency of each mode depend upon temperature.Note that I 2 B(Ω) is exactly zero at zero energy in this model.
When we obtain the optical scattering rate, we can go further and get the complex optical conductivity 20 .Because the optical scattering rate is identical to -2 times the imaginary part of the optical self-energy, i.e. 1/τ op (ω) ≡ −2Σ op 2 (ω) and the real and imaginary parts of the optical self-energy form a Kramers-Kronig pair, we can get the real part of the optical self-energy from the imaginary part using the Kramers-Kronig relation between them, i.e. 20 , where P is the principle part of the improper integral.The Kramers-Kronig relation is a Cauchy principal value integral.And then, we obtain the complex optical conductivity from the complex optical self-energy using the extended Drude model, Eq. (1).

III. RESULTS AND DISCUSSION
Fig. 1 shows our model EBSD functions (I 2 B(Ω)) at various temperatures ranging from 0 to 500 K with a 10 K increment.The EBSD function consists of two components, namely, sharp and broad modes (Eq.( 4)).Note that the sharp mode linearly shifts to a higher energy and its intensity (or area) linearly decreases as temperature decreases, whereas the broad mode linearly shifts to a higher energy and its intensity linearly increases as temperature decreases.As the temperature is reduced, the total area of the two modes gradually  that the thick purple is 0 K, the thick blue curve is 100 K, the thick green curve is 200 K, the thick orange curve is 300 K, the thick red curve is 400 K, and the thick dark brown curve is 500 K.In the inset, the temperature-dependent coupling constant, λ(T ), is shown.self-energy at temperatures ranging from 0 to 500 K with a 10 K increment.Note that the thick purple is 0 K, the thick blue curve is 100 K, the thick green curve is 200 K, the thick orange curve is 300 K, the thick red curve is 400 K, and the thick dark brown curve is 500 K.Note that the real part was obtained from the imaginary part using the Kramers-Kronig relation between them (see the text for a more detailed description).) Energy (meV)    ) Energy (meV) imaginary parts of the optical conductivity at temperatures ranging from 0 to 500 K with a 10 K increment.In the inset, the frequency-and temperature-dependent impurity scattering rates at several selected temperatures are shown.Note that the scattering rate is normalized by 1/τ imp .
The optical conductivity is obtained from the optical self-energy using the extended Drude model with Ω p of 2.0 eV and 1/τ imp of 15 meV (see the text for a more detailed description).
Figs. 6(a) and (b) show the real and imaginary parts of the optical conductivity obtained from the real and imaginary parts of the optical self-energy using the extended Drude model.
To get the optical conductivity, the plasma frequency (Ω p ) was set at 2.0 eV, and the impurity scattering rate was included in −2Σ op 2 (ω).For the case of pseudogap, the impurity scattering rate is energy-and temperature-dependent 26 as shown in the inset of Fig. 6(a When the material enters the superconducting state below the SC transition temperature (T c ), the electrons are paired and condensed for the macroscopic SC ground state, and simultaneously the SC gap is opened.In this study, we only deal with the normal state, including the pseudogap state.Therefore, our discussions in this study are valid only at temperatures above T c in the case of material systems.Also, note that the total scattering rate (−2Σ op,Total 2 (ω)) shows similar ξ 2 (≡ (hω) 2 + (pπk B T ) 2 ) dependencies as the optical scattering rate (−2Σ op 2 (ω)) (see Fig. S4 in the Supplementary Material), where h is the reduced Planck constant, k B is the Boltzmann constant, and p is an adjusting parameter, which may depend upon the material systems 4, 33,34 .We obtained the amplitude and phase angle of the optical complex conductivity and showed them in Fig. 4(c) and (d).The results seem to be consistent with those of underdoped cuprate in a published paper 5 .
We fitted the optical conductivity data at five selected temperatures with the Drude-Lorentz model up to 300 meV and showed the results in Figs.7(a-e).As we can see in the figure, the conductivity consists of two components: a Drude component and a Lorentz one.
The Lorentz one is known as the mid-infrared absorption (or the conformal tail) 5,6,35 .The results are consistent with those of underdoped cuprates 5,6 ; the quality of the fit is getting better as the temperature increases.In Fig. 7(f), we show the fitting parameters as functions of temperature.The Lorentz mode is almost temperature-independent, whereas the Drude mode exhibits a strong T -dependency, particularly the width.These T -dependent trends of fitting parameters are consistent with those of underdoped cuprates in the literature 6 .We also fitted the optical conductivity at the same selected temperatures in the case without including the PG with the same scheme as in the reference [6], i.e., using an almost Tindependent Lorentz mode and a T -dependent Drude model and showed the results in Fig. S3 in the Supplementary Material.The quality of the fits is not as good as in Fig. 7 and is getting worse as the temperature increases.Therefore, the results are not consistent with those of underdoped cuprates in the literature 6 .
Furthermore, in Figs.Note that the thick blue curve is 100 K, the thick green curve is 200 K, and the thick orange curve is 300 K.The dash-dotted lines in (b) and (d) are guides for the eyes.
optical scattering rates at temperatures between 50 and 300 K are shown for including only the pseudogap phase and reasonable pseudogap depths.When p is 2.0, all curves of the optical scattering rates at various temperatures fall into a single curve at values of ξ 2 below ∼10000 meV 2 and are linear to ξ 2 , indicating that the optical scattering rates show a Fermi liquid behaviour in the pseudogap phase.This result is very similar to the experimentally observed results of underdoped cuprate systems 4 .We could not find the reason why our data show a p-value of 2.0 rather than values less than 2 as in the published papers 4,6 yet.
It is worth noting that the optical scattering rate in this paper (1/τ op (ω) = −2Σ op 2 (ω)) is different from the optical scattering rate (h/τ (ω)) in the literature by Kumar et al. 6  We also performed the same calculations with a temperature-independent I 2 B(Ω) (see The overall temperature-dependent exhibits a strong T -dependence and the T -dependent change of the imaginary part of the optical self-energy is significantly greater than that of the measured optical scattering rate 4 .We obtained the real part from the imaginary part using the Kramers-Kroing relation between them.We also showed the total optical self-energy in Fig. S9 for the case of including the PG, which are used to get the optical conductivity.We also obtained the optical conductivity spectra from the optical self-energy using the extended Drude model and showed them in Figs.S7 and S10.The real part of the optical conductivity in the case of including the PG (Fig. S10(a)) differs from the measured conductivity 4, 15,32 .We also showed the optical scattering rate as a function of ξ 2 in Fig. S11.The curves are very widely spread vertically, even though they fall into a single line at low ξ 2 values below ∼6000 meV 2 (Fig. S11(b)) for the case of including PG.We also show the total optical scattering rate as a function of ξ 2 in Fig. S12.These quantities differ from those observed experimentally 4 .These results support the fact that the I 2 B(Ω) should be temperature dependent for simulating the experimentally observed temperate-dependent optical spectra using the procedure described in this paper.
Our results indicate that the observed Fermi liquid-like behaviour is intimately associated with both the T -dependent EBSD function and the T -dependent pseudogap.The two physical quantities directly contribute to the optical scattering rate.The EBSD function represents correlations in a correlated electron system such as cuprates and Fe-based superconductors, whereas the pseudogap reduces the density of states near the Fermi level, resulting in a partial gap near the Fermi level.Therefore, these two quantities give oppo-site contributions to the optical scattering rate.Generally, the EBSD function enhances the optical scattering due to the correlations (or scatterings) between electrons, whereas the pseudogap suppresses the optical scattering rate in the pseudogap region due to the reduction of scatterers.The EBSD function gives an up-step feature when one measures from zero energy to high energy, whereas the pseudogap gives a down-step feature when one measures from high energy above the pseudogap energy to low energy.Both give a step-like feature even their line shapes are not the same.However, it is not easy to distinguish the two contributions clearly in a measured optical spectrum.In this paper, we separately show the two features in the optical scattering rate using the approach described in this paper.We also show that the two quantities are temperature dependent and the combined contributions of the two quantities give the Fermi liquid-like behaviour observed experimentally.Therefore, the observed Fermi liquid-like behaviour can be understood by considering the combined contributions from both the correlations between electrons and the pseudogap phenomenon.There were some interesting optical studies on Fermi-liquid-like behavior and strange (or bad) metal state in the optical spectra of cuprates 5,6 .These studies can be categorized into two groups based on the methods of analysis 35 : one is a one-component analysis, and the other is a two-component analysis.These analyses are also called the Mathiessen interpretation (one-component analysis) and anti-Mathiessen interpretation (two-component analysis) 5 .In these previous works 5,6 , the authors used the two analyses to analyze the measured optical spectra in a wide doping range, from underdoped to overdoped, whereas in our study, we simulated the optical self-energy and optical conductivity of an underdoped cuprate starting from the T -dependent EBSD function using the generalized Allen's formula and extended Drude model.The optical scattering rate (the imaginary part of the optical self-energy) exhibits Fermi-liquid-like behavior (see Fig. 8), and the optical conductivity shows two components: the Drude component and the mid-infrared absorption (or conformal tail) (see Fig. 7).The T -dependent evolutions of these two components are consistent with those of measured underdoped cuprates 5,6 .Additionally, in our approach, the correlations between electrons were explicitly revealed through the EBSD function.

IV. CONCLUSION
We investigated the Fermi liquid-like behaviour observed by optical spectroscopy4 in the pseudogap phase of underdoped cuprates.We applied the reverse process of a well-established approach, which has been used to extract the EBSD function from the measured optical spectrum, to obtain the optical scattering rate in the pseudogap phase of cuprate systems.We successfully simulated the observed Fermi liquid-like behaviour using the reverse process with the temperature-dependent EBSD function and Fermi T -dependent, particularly the width.These results are consistent with previous optical studies 5,6 .Furthermore, our results indicate that the T -dependent EBSD functions extracted from the measured optical spectra of cuprates using the established approach are intrinsic and contain information regarding the glue spectrum for the Cooper pair formation.
decreases.The detailed temperature-dependent changes of the all parameters in the T -

Figure 1 :
Figure 1: (Color online) Electron-boson spectral density (EBSD) functions at temperatures ranging from 0 to 500 K with a 10 K increment.The EBSD function consists of two components.Note Fig.2(a)shows the optical scattering rates (or the imaginary parts of the optical selfenergy) at various temperatures ranging from 0 to 500 K obtained from the temperaturedependent model I 2 B(ω, T ) shown in Fig.1using Eq. (2) and Eq.(3) with N 0 (T ) = 1.0, i.e. without including the pseudogap.At each temperature, we observe a rapid increase near the peak position of I 2 B(Ω) resulting in a step-like feature.The step-like feature weakens with increasing temperature, which is consistent with the temperature-dependent intensity of the sharp mode (Fig.1).In most regions, the scattering rate increases monotonically increases as the temperature increases.Note that when the peak of I 2 B(Ω) is very sharp, the order of the scattering rate can be reversed, as seen at low temperatures above ∼300 meV.Fig.2(b)shows the corresponding real parts of the optical self-energy at various temperatures which were obtained from the imaginary parts using the Kramers-Kronig relation between the real and imaginary parts.A systematic temperature-dependent trend is shown.A broad peak is observed at the energy where the maximum slope in the scattering rate is located.The peak is shifted to higher energy as the temperature increases.However, the shape of the peak is different from that of the peak in the measured real part of the optical self-energy of underdoped cuprate systems 4,15,19 .Figs.3(a) and (b) show the real and imaginary parts of the optical conductivity, respectively, which were obtained from the real and imaginary parts of the optical self-energy and the extended Drude formula, Eq. (1).To obtain the optical conductivity, the plasma frequency (Ω p ) was set at 2.0 eV, and the impurity scattering rate (1/τ imp ) of 15 meV was

Figure 2 :
Figure 2: (Color online) Complex optical self-energy, without including the pseudogap.(a) Optical scattering rates (imaginary parts of the optical self-energy) and (b) the real parts of the optical included.The temperature-dependent real part of the optical conductivity is significantly different from that of the measured one in the paper by Mirzaei et al.4  .We also obtained the amplitude (|σ(ω)| ≡ [σ 1 (ω)] 2 + [σ 2 (ω)] 2) and phase angle (tan −1 [σ 2 (ω)/σ 1 (ω)]) of the complex optical conductivity5,31,32 and showed them in Figs.4(a) and (b).Interestingly, the results seem to be consistent with those of underdoped cuprate in a published paper5 , even though the conductivity is much different from that of a measured underdoped cuprate 4,15 , which may indicate that the pseudogap does not affect the conformal tail (or the mid-

Figure 3 :Figure 4 :
Figure 3: (Color online) Complex optical conductivity, without including the pseudogap.(a) Real and (b) imaginary parts of the optical conductivity at various temperatures ranging from 0 to 500K with a 10 K increment.The optical conductivity was obtained from the optical self-energy using the extended Drude model with the plasma frequency of 2.0 eV and 1/τ imp of 15 meV.Here, the impurity scattering rate is independent of both frequency and temperature.

Figs. 5
Figs.5(a) and (b), respectively, shows the real and imaginary parts of the optical selfenergy at various temperatures ranging from 0 to 500 K obtained from the temperaturedependent model I 2 B(Ω, T ) using Eq.(2) and Eq.(3) with N 0 (T ) ≤ 1.0, i.e. with including the pseudogap for T ≤ T * , where T * is the pseudogap temperature (in our case T * = 300 K).To include the temperature-dependent pseudogap, we used the Fermi arc model22 , i.e.N 0 (T ) = T /T * for T ≤ T * in Eq. (3) and N 0 (T ) = 1.0 for T > T * .The model pseudogap is

Figure 5 :
Figure 5: (Color online) Complex optical self-energy, including the pesudogap.(a) Optical scattering rates (or imaginary parts of the optical self-energy) (b) real parts of the optical self-energy at various temperatures ranging from 0 to 500 K with a 10 K increment.In the inset, the model pseudogap is shown for a case of N 0 = 0.5.The real part was obtained from the imaginary part using the Kramers-Kronig relation between them.

Figure 6 :
Figure 6: (Color online) Complex optical conductivity, including the pseudogap.(a) Real and (b) Fig. S2(b) in the Supplementary Material).It should be noted that the pseudogap state is a normal state with an onset temperature (T * ) below which the pseudogap state appears.

Figure 7 :
Figure 7: (Color online) (a-e) The optical conductivity data and Drude-Lorentz model fits at five selected temperatures in the PG phase.The data were fitted up to 300 meV with a Drude mode and a Lorentz mode.(f) The temperature-dependent fitting parameters are the Drude plasma frequency (Ω Dp ) and width, and the Lorentz center frequency, width, and plasma frequency (Ω Lp ).Note that the center frequency and width of the Lorentz mode are fixed.Both the Drude and Lorentz spectral weights show almost no temperature dependence.
8(a) and 8(b), we plot the optical scattering rates as functions of ξ 2 in the case of including the pseudogap for p = 1.5 and 2.0, respectively.Note that the

Fig. S5 in
Fig. S5 in the Supplementary Material).The T -independent coupling constant (λ) as a function of temperature is shown in the inset of Fig. S5.All the results are shown in Figs.S5-S12 in the Supplementary Material.The optical scattering rates (or imaginary parts of the optical self-energy) and the real parts of the optical self-energy (see Figs. and (b) and S8(a) and (b)) significantly .
arc modelled pseudogap.Therefore, the observed Fermi liquid-like behaviour can be understood by considering the combined contributions from the T -dependent EBSD function and the T -dependent pseudogap.Additionally, we obtained optical conductivity from the optical self-energy using the extended Drude model.The optical conductivity in the low-energy region shows two components: a Drude component and a Lorentz one.The Lorentz component is almost T -independent, whereas the Drude one is strongly