The Ancient Romans’ Route to Charge Density Waves in Cuprates

: An account is given of the main steps that led the research group in Rome, to which the author belongs, to the formulation of the charge-density-wave scenario for high- T c superconducting cuprates. The early ﬁnding of the generic tendency of strongly correlated electron systems with short range interactions to undergo electron phase separation was subsequently contrasted with the homogenizing effect of the long-range Coulomb interaction. The two effects can ﬁnd a compromise in the formation of incommensurate charge density waves. These charge density waves are inherently dynamical and are overdamped as a consequence of the possibility to decay in electron-hole pairs, yet tend to maintain a (quantum) critical character, which is mirrored in their marked momentum and frequency dependence and in their strong variation with temperature and doping. These dynamical incommensurate charge density waves act as mediators of pairing lading to high- T c superconductivity, and provide the scattering mechanism that produces the observed violation of the Fermi-liquid paradigm in the metallic phase. corrections due to processes with electron-hole pairs in the virtual intermediate state, contributing to Raman and optical spectra.


Introduction
High-T c superconducting cuprates represent one of the most intriguing challenges of condensed matter physics [1]. Several ingredients and experimental facts contribute to create a halo of mystery around these systems. The parent compounds are antiferromagnetic insulators and all the families share the same perovskite structure based on copper-oxygen (CuO 2 ) planes intercalated with rare-earth slabs. Upon substitution with of the rare earth atoms with heterovalent dopants, charge carriers (holes or electrons, in the following, for simplicity, I discuss the hole-doped case) are introduced in the copper-oxygen planes and the system becomes a metal, although the metallic properties are anisotropic, the resistance along the copper-oxygen planes being smaller that the resistance in the direction perpendicular to them. Moreover, the resistance displays a linear temperature dependence over a very wide temperature range (except for very heavily doped samples), that is hardly reconcilable with the paradigmatic Fermi-liquid behaviour of normal metals. Superconductivity, of d-wave symmetry, occurs below a dome-shaped doping-dependent critical line T c (p), that reaches its maximum at optimal doping p opt ⇡ 0.16 and declines with further doping, in the so-called overdoped region, where the metallic phase above T c is well described by a Fermi liquid theory. In the underdoped region, p < p opt , and in a temperature range T c (p) < T < T ⇤ (p), where T ⇤ (p) decreases with p and merges with T c near optimal doping, various phenomena apparently connected with a suppression of the density of states at the Fermi level are detected (the so-called pseudogap phenomena), that indicate an even more pronounced violation of the Fermi liquid paradigm. Along the years, the research group in Rome, to which the author belongs, (the ancient romans named in the title of this piece of work) developed a scenario that tried to identify the crucial ingredients to account for the above complicated phenomenology, in the perspective of a modification of the Fermi liquid properties while the doping is reduced, starting from the overdoped region. My aim here is to retrace the ancient romans' route to a theory for the anomalies of the metallic state of cuprates and their possible connections with superconductivity, that relies on charge density waves [2][3][4][5].
It must be emphasized that, over more than three decades, other groups in Rome have been working in the field of high-T c cuprates. Their role and relevance is witnessed by the many references cited in this piece of work. It is likewise evident that many scenarios have been developed for cuprates along all those years. The focus of this piece of work is the development of the scenario developed by the research group in Rome, to which the author belongs, that is based on charge density waves, but it may be interesting to compare this scenario at least with similar scenarios. The scenario that shares more similarities with the one presented in this piece of work relies on nearly antiferromagnetic critical spin fluctuations [6][7][8][9][10], remnants of the proximity to the antiferromagnetic phase. These fluctuations have characteristic wave vectors substantially larger than those of the charge density waves, and pointing along the (1, 1) direction of the Brillouin zone, so a careful analysis often allows to separate their contribution to various observed quantities from the contribution of charge density waves with characteristic wave vectors along the (1, 0) and (0, 1) directions of the Brillouin zone.
This paper is organised as follows: Section 2 is devoted to illustrate the main ideas behind supporting the possible occurrence of electron phase separation in cuprates; in Section 3 the homogeneising role of Coulomb interaction is taken into account, turning electron phase separation into incommensurate charge density waves; Sections 4-6 deal with the long search of evanescent fingerprints of the charge density waves in various experimental probes, like angle resolved photoemission spectroscopy, transport, optics, Raman spectroscopy, and other indirect evidences; in Section 7 the role of resonant X-ray scattering in providing direct evidence of charge density waves is discussed; final remarks are found in Section 8.

The Electron Phase Separation Era
The fact that the undoped parent compounds are antiferromagnetic insulators, was immediately recognised as the consequence of the strong correlation U in the state with two holes residing on the same copper d orbital [11]. The configuration with two holes on neighbouring sites is much less costly and the energy is further lowered if the spins of the two holes are antiparallel, because virtual hopping processes are possible in this configuration, with an energy gain of the order J ⇠ t 2 /U, where t is the hopping amplitude. This is the famous superexchange mechanism. Then, in the early days of high-T c superconductivity, a great effort was devoted to the investigation of the properties of various singleand multi-band models for strongly correlated electrons on the CuO 2 planes.
It became rapidly evident that a property common to all those models was a marked tendency toward electron phase separation [12][13][14][15][16]. This can be understood as the consequence of the suppression of the main mechanism that leads to a uniform electron density in normal metals, namely the kinetic energy gain when the electrons move freely across the system. Strong correlations severely reduce this energy gain so that residual attractive forces of various origin can bring the system on the verge of an instability that favours electron phase separation in hole-rich (more metallic) and hole-poor (more insulating) regions.
One possible mechanism stabilising the system against electron phase separation was found to be superconductivity [12,17], pointing to the possibility that the same attractive forces responsible for electron phase separation in the particle-hole channel, may lead to Cooper pairing in the particle-particle channel.

The Charge Density Wave Era
In its crude simplicity, the electron phase separation scenario discussed in the previous section misses a crucial ingredient, namely the charged character of doped holes. It is rather clear that long-range Coulomb repulsion makes charge segregation energetically very costly, so that genuine electron phase separation can only occur in systems where the countercharges that guarantee the overall charge neutrality of the system are sufficiently mobile [18]. However, the system may reach a compromise between these two opposed tendencies and produce charge segregation on short length scales. This is the so-called frustrated electron phase separation mechanism [2,3], that leads to the formation of incommensurate charge density waves, for which an increasing experimental evidence was being gathered [19][20][21]. The electron phase separation mechanism for charge density waves was shaped as an instability of the Fermi liquid by the research group in Rome, to which the author belongs [22][23][24][25][26].
The idea may be schematised as follows: in a correlated Fermi liquid one can identify three main contributions to the electron-electron interaction, V ee = V LR + U SR + L PS , where V LR is the long-range Coulomb repulsion, U SR is a short-range residual repulsion that embodies the strongly-correlated nature of the parent undoped insulator, and L PS is the residual attraction leading to electron phase separation if long-range Coulomb forces are neglected. In a Fermi liquid, screening processes are active and the electron-electron interaction is modified by the polarisation associated with virtual electron-hole pairs, leading to a screened interaction where P is the electron-hole polarisation bubble in the random phase approximation (also called Lindhard function). This is a complex function that includes damping effects resulting from the decay process involving electron-hole pairs [see below, Equation (2), and the text therby]. Within this picture, instabilities are seen as zeroes in the denominator of Equation (1). Since U SR is minimum at small wave vectors q, while V LR diverges for |q| ! 0 and decreases at larger q, if the residual attraction L PS has no particular structure as a function of q, then the denominator of Equation (1) may vanish at |q| = 0 along a certain line T PS (p) in the temperature-vs.-doping phase diagram, marking the electron phase separation instability of the Fermi liquid. This is simpler to achieve if the electron bandwidth is suppressed by strong correlation effects, thereby enhancing the polarisation bubble P. However, electron phase separation can be achieved only if the term V LR is neglected. In the presence of V LR , the instability at |q| = 0 is strictly impeded by the divergence of the unscreened Coulomb interaction at small wave vector, but an instability may yet be achieved along a line T 0 CDW (p) < T PS (p), marking a finite wave vector instability of the Fermi liquid, at q = q c (actually, the instability occurs at a whole star of wave vectors that are equivalent under point group symmetries of the lattice). Since this wave vector is the result of the form of V ee and P, as functions of q, the instability wave vector q c is unrelated to the periodicity of the lattice: the compromise between the tendency toward electron phase separation and Coulomb repulsion is realised as an incommensurate charge density wave.
When Equation (1) is expanded for small frequency and wave vectors close to q c , then the effective interaction takes the singular form where g 2 is a constant setting the overall strength of the interaction, in the field-theoretic language, is called the propagator of the collective excitations mediating the effective interaction e V ee (in our case, charge density waves), m 0 µ T T 0 CDW (p) vanishes along the instability line T CDW (p), marking the divergence of the correlation length x 0 = p n/m 0 , n is a constant whose order of magnitude is set by an electron energy scale (the Fermi energy # F ), and G 1 ⌘ W sets a high energy cutoff, that depends on various energy scales. For instance, if Equation (1) is derived within the Hubbard-Holstein model including the coupling to a phonon with frequency w 0 [3], then W ⇠ w 2 0 /# F < w 0 . Apart from an overall multiplicative constant, this is the same form of the dressed density-density dynamical correlation function where the bare term is c 0 rr ⌘ P, and the result of the right-hand side is the resummation of the diagrams in Figure 1a. This is exactly the form of a dynamical Ornstein-Zernike correlator, the term iw representing the Landau damping due to charge density waves decaying in electron-hole pairs in the Fermi liquid state, giving an alternative but equivalent interpretation of the charge density wave propagator D in Equation (3). The details of the band structure of cuprates, determining P, and of the effective interaction V ee , fix the direction of q c along the CuO bonds consistent with all the evidence existing at the time when the theory was developed, as well as with further evidence gathered in the subsequent years. Of course, as we stated above, the magnitude of q c is not related to a periodicity of the lattice, therefore the resulting charge density wave is incommensurate.
The above scenario is still rather qualitative. As a matter of fact, the loosely layered structure of cuprates makes fluctuations play an important role: the interaction between dynamical charge density wave [which can be seen as the result of corrections to Equation (1) beyond the random phase approximation] suppresses the order in a temperature region T CDW (p) < T < T 0 CDW (p), so that charge density waves are expected to stay dynamical and fluctuating until the real instability line is reached [4]. This can be seen as the correction to the m 0 term in Equation (2), where u is a dimensionless coupling for dynamical charge density waves (see Figure 1b, where the interaction between two ingoing and two outgoing charge density waves (broken lines) is mediated by virtually excited electrons in the medium (solid line); at low energy this interaction is a constant u) and the correlation length x = p n/m diverges at T = T CDW (p) < T 0 CDW (p), marking the true instability line [4,5,27]. Further corrections dressing the vertices in the diagram of Figure 1b are expected to be small for incommensurate charge density waves [28]. The resulting T CDW (p) instability line may be dome shaped, and completely masked underneath the superconducting dome T c (p) [5]. This result is obtained, e.g., in a model where the dynamical charge density wave in the underdoped region have a different nature, the modulus of the complex charge density wave order parameter being well formed (contrary to the case of the optimally and overdoped region, where it fluctuates) while its phase fluctuates.
Of course, if a coupling to the lattice degrees of freedom is included, the maximum of T CDW (p) will likely adjust to match a nearby nesting condition between the (doping-dependent) period of the charge density wave and the periodicity of the lattice, which in the doping region of interest may occur at p ⇡ 1/8.
The resulting theoretical phase diagram looks like the one sketched in Figure 2.

Looking for Fingerprints: Angle Resolved Photoemission Spectroscopy
Once the possibility for a charge density wave as an instability of the Fermi liquid was assessed, the question of detecting fingerprints of this state naturally arose. Charge density waves in cuprates involve rather tiny modulations of the charge density profile, presumably because strong correlations reduce the charge density wave amplitude [29], and have a dynamical character over most of the phase diagram above T c (p) [4,5,27], being therefore rather hard to observe directly (of course, one could find traces of the electron charge modulation in the accompanying lattice deformation, but here my focus is the direct observation of charge density waves), unless favourable conditions like nesting or pinning occur, and charge density waves are seen, usually in the form of charge stripes, with a highly anharmonic profile. In the early times of the detection of charge density waves in cuprates, this elusive phenomenon had to be studied mostly by means of indirect probes, apt to reveal modifications of the electron properties caused by the interaction between electrons and charge density waves.
The first fingerprints of charge density waves were looked for in angle resolved photoemission spectra [30]. The idea is to calculate the electron self-energy correction due to the interaction with dynamical charge density waves (see Figure 1c) and try to recover the various feature observed in angle resolved photoemission spectra by adjusting the parameters that characterise the dynamical charge density waves within reasonable ranges of variation, as well as entering values extracted from other experiments, see, e.g., Section 5.
It must be pointed out that angle resolved photoemission spectra of cuprates (as well as many other spectral porperties) are also affected by nearly antiferromagnetic spin density fluctuations [30,40], with their characteristic wave vector q AFM c ⇡ (p, p). The specificity of charge density waves is that their characteristic wave vector q c is along the CuO bonds, in the (1, 0) and (0, 1) directions of the Brillouin zone, and contribute with features that cannot be associated to scattering processes involving a mediator with a larger characteristic wave vector, q AFM c , pointing along the diagonals of the Brillouin zone. Furthermore, charge density waves introduce lower energy scales that are appropriate to describe the kinks (while waterfalls have both contributions from charge and spin fluctuations). As it was recalled in Section 3, if Equation (1) is derived within the Hubbard-Holstein model including the coupling to a phonon with frequency w 0 [3], then the characteristic energy scale W of a charge density wave turns out to be W ⇠ w 2 0 /# F . The role of charge density waves in the formation of the pseudogap was also investigated. The idea is that the interaction in Equation (2) is apt to gap electron states near the M points of the Brillouin zone [42][43][44]. The pseudogap temperature T ⇤ is interpreted in terms of a reduction of the density of states due to incipient charge density waves and, at lower temperature to the possible formation of incoherent superconducting pairs [45].

Looking for Fingerprints: Transport, Optics, Raman Spectroscopy
Of course, the theory provides only estimates of the order of magnitude of the various parameters characterising the charge density waves, see Equations (2) and (4), so other independent determinations of the parameters g, n and W = 1/G, and an experimental determination of the evolution of m with temperature and doping became highly desirable.
A major breakthrough came from the idea that nearly critical dynamical charge density waves would affect the Raman response with a channel associated to the possible emission of two virtual charge density waves in the intermediate state [46] (see Figure 1d). It was soon realised that this contribution had the right form and the right symmetry properties to explain the anomalous peak observed in Raman experiment [47].
A thorough analysis of Raman spectra began [46,[48][49][50][51][52][53], that enjoyed feedbacks from, and provided further inputs to, the analysis of angle resolved photoemission spectra, discussed in Section 4, allowed to extract the first reliable experimental estimates of the characteristic parameters of the charge density waves. It was confirmed that n ⇡ 1 eV is an electron energy scale, while W ⇡ 50 meV is a much lower energy scale. Interestingly enough, the analysis of Raman spectra confirmed the fact that charge density waves have a very weak tendency to order and stay dynamical an fluctuating over most of the accessible phase diagram above T c (p), confirming their elusive character. The parameter m, that should vanish at the charge density wave phase transition, is found to be a linear function of the temperature at sufficiently high temperature, but is always found to saturate to a small, but finite value (⇠ 1-10 meV) at low temperature. Extrapolation of the linear part gives an estimate of the temperature at which the charge density waves would tend to order. This temperature scale is negative (i.e., unattainable) in overdoped samples, nearly vanishing at optimal doping, thus marking the presence of a missed quantum critical point, and finite, but smaller than or comparable to T c (p) in moderately underdoped samples.
A subsequent deeper analysis of Raman spectra [53,54] led to the complete characterisation of the glue function a 2 F that should mediate pairing between electrons. This glue function was found to be mostly due to charge density waves in optimally an overdoped samples, whereas upon underdoping, the role of antiferromagnetic spin fluctuation increases, as it is expected, since the system is now closer to the antiferromagnetic phase. It is interesting that the strength of the two channel becomes equal around optimal doping. These spin fluctuations do contribute to the Raman spectra, but can never account for the anomalous peak via a process like the one depicted in Figure 1d, because they never become nearly critical over all the doping region where the anomalous peak is observed.
Subleading contributions to Raman spectra in channels where the intermediate state always contains an electron-hole pair (see Figure 1e), as well as contributions to optical spectra and to electron transport have provided further evidence of the presence of charge density waves [48][49][50]. These evidences being weaker than those gained by the analysis of Raman and angle resolved photoemission spectra, had mostly the role of testing the consistency of the parameters extracted thereby. It must be emphasized that also spin fluctuations do contribute to Raman spectra with subleading corrections, and this is precisely the reason why the full analysis of Raman spectra allowed to disentangle the spin and charge contribution to the glue function a 2 F.
For instance, the most important contribution to Raman spectra (see Figure 1d) does not contribute to optical spectra, because of the neutral character of the dynamical charge density waves in the intermediate state, unless the conditions occur (e.g., in the presence of strong disorder) to prevent this term from vanishing [55,56]. Likewise, self-energy and vertex corrections (see Figure 1e) to the optical conductivity (and to transport) crucially depend on the cutoff W as a consequence of the strong constraints imposed by charge conservation: an exact cancellation of self-energy and vertex corrections occurs, that is enforced by a Ward identity, if no retardation effect due to degrees of freedom different from the electrons themselves is present in the effective interaction V ee . If one includes retardation effects in L PS , e.g., due to phonons, then the charge density wave contribution to optical spectra (and transport) does not vanish [57].

Other Fingerprints
Yet other indirect fingerprints of the charge density waves are found in other contexts. For instance, the interplay of charge and spin degrees of freedom in the heavily underdoped region, near the antiferromagnetic state, can give rise to charge segments [58] that can be seen as the precursor (the melt [59]) of stripy charge density waves and may account for the properties of the nematic state reported in this region of the phase diagram [60][61][62][63]. Dynamical charge density waves may coexist with (or evolve into) nematic fluctuations affecting the Raman spectra [64].
An even more important context where fingerprints of the charge density waves may be detected is in the connection and interplay with superconductivity.
Since the early days of the electron phase separation era discussed in Section 2, it was found that superconductivity might stabilise the system against phase separation, indicating that the same attractive forces responsible for the instability in the particle-hole channel, could instead produce Cooper pairing in the particle-particle channel [12,17].
When the charge density wave era discussed in Section 3 began, it was likewise clear that dynamical charge density wave described by Equation (4) could mediate Cooper pairing in the particle-particle channel [65], contribute to the properties of the pseudogap state [42][43][44], explain the possible occurrence of a time-reversal-symmetry-broken superconducting state [66,67], and be possibly related to a cutoff of the damped Cooper pair fluctuations contributing to the paraconductivity above T c [68], within the Aslamazov-Larkin theory.
Even more interestingly, in a recent work, the possibility that in the underdoped region superconductivity and charge density waves might be different breaking of a common symmetry, led to the proposal that the occurrence of filamentary superconductivity might be promoted by disorder in a region of the phase diagram where charge order would be present in the absence of disorder. This phenomenon is related to the fact that, in order to match with the different phases of the charge order parameter in different domains induced by disorder, the order parameter must rotate to the superconducting state at the domain boundary, thus producing topologically protected filaments [69].
The fact that superconductivity and charge order compete is witnessed by the fact that once superconductivity is suppressed by a strong enough magnetic field, static charge order emerges below a dome-shaped curve T CDW (p) [70,71] very similar to the one reported in Figure 2, that was hidden underneath the superconducting dome T c (p), at zero magnetic field. Interestingly enough, superconductivity seems to be more resilient around the two endpoints where the curve T CDW (p) meets the axis T = 0 [72]. This seems to indicate that quantum fluctuations of the competing phase act indeed as mediators of Cooper pairing in the particle-particle channel, like it happens in several heavy fermion superconductors. If this were the case, the two superconducting domes persisting around the endpoints of the curve T CDW (p) at strong magnetic field are two-different superconducting phases, in which pairing is mediated by different quantum critical fluctuations. Upon reducing the magnetic field, the two dome expand and eventually coalesce in what appears as a single superconducting phase at zero magnetic field. The two plateaux observed, e.g., at T c ⇡ 60 K and T c ⇡ 90 K in YBCO should be the remnant of the two different superconducting phases. It must be emphasized that spin fluctuations will certainly contribute to Cooper pairing, especially in the underdoped region (see the discussion on the glue function in Section 5). However, since the quantum critical point for spin fluctuations is outside the superconducting dome and not underneath it, the quantum critical fluctuations boosting T c must be related with the quantum critical point(s) underneath the superconducting dome, associated to charge density waves.
It is interesting to discuss here the results of another group in Rome, led by A. Bianconi. Following the appearance of photo-induced X-ray diffraction satellites in cuprates at low temperature below 150 K, it has been shown the predicted self-organization of anisotropic polarons into striped puddles [73][74][75]. The data analysis has shown that photoinduced polarons have an anisotropic shape with a finite size of about 1 nm; it has been inferred that they are in the intermediate coupling regime and coexist with itinerant carriers in the doping range between 0.06 and 0.21 holes per Cu site. From these experimental results, it was deduced that the intermediate polarons show internal fluctuations which generate an anisotropic polaron-polaron van Der Waals attraction giving a fluid-like behavior exhibiting a complex phase separation. They have proposed a phase diagram formally similar to the gas-to-liquid-liquid transition in supercooled water, where the anisotropy interaction generates a novel intermediate liquid phase between the standard gas and liquid phases [76][77][78].
The intermediate phase has been identified with the 1/8 polaronic generalized Wigner charge density wave phase. In this scenario, it was proposed that the superconducting phase is determined by the complex filamentary network of interface pathways which could be mapped by an hyperbolic geometry [79].

The Resonant X-Ray Scattering Era
An exciting new era begun when resonant X-ray scattering experiment were gradually refined to reach the accuracy and sensitivity needed to directly probe charge density waves [80][81][82][83]. In a rather short lapse of time, what that far had only been seen in few experiments and under peculiar conditions, became ubiquitous and charge density waves were now finally seen, beyond any shadow of doubt, in all family of cuprates. They may even re-enter in the overdoped region, in those families of cuprates where the van Hove singularity in the electron spectrum is known to play a role [84], and may indeed enhance the polarisation bubble P in Equation (1) if the doping brings the van Hove singularity sufficiently close to the Fermi level, promoting the reappearance of charge density waves to the right of p 0 c in Figure 2. The most important fact, here, is that resonant X-ray scattering processes may excite a single charge density wave in the intermediate channel, so that the absorption spectrum is directly proportional to the dressed density-density dynamical correlation function in Equation (4). Then, for the first time, experimental data may give direct access to the parameters characterizing the charge density waves After the phase diagram sketched in Figure 2 had been published [5], purely based on theoretical calculations, with the parameters characterising the charge density waves extracted from the various indirect evidences discussed in Sections 4-6, a thorough analysis was carried out in YBCO and NBCO at various dopings and temperatures. The analysis not only confirmed the theoretical phase diagram, but also provided values for n and W, and allowed to trace the evolution of m with varying the temperature, in very good agreement with all previous indirect estimates [85].
Surprisingly, beyond the dynamical charge density wave, whose intensity decreases with the temperature and becomes unobservable at some T dyn CDW (p) < T ⇤ (p), very robust and nearly temperature independent charge density fluctuations, characterised by correlation lengths much shorter that those of the charge density waves, were found to pervade the phase diagram up to the highest temperature reached in the experiment. It is not clear at the moment, whether the two kind of excitations coexist or are spatially separated (giving rise to an inhomogeneous landscape, see, e.g., [86]), but it is quite evident that charge density waves compete with (and are weakened by) superconductivity, whereas the much broader charge density fluctuations are apparently ataraxic, hardly disturbed by anything, included superconductivity, with whom they seem to peacefully coexist.
It might be speculated that the broad charge density fluctuations may be remnant of the re-entrant charge density waves due to the proximity to a van Hove singularity, whenever the singularity is not close enough to the Fermi level and only produces broad resonances in the dressed density-density dynamical correlation function in Equation (4). This hypothesis is presently under scrutiny.

Concluding Remarks
Despite the fact that the pace of convergence between an old theory and the experiment significantly accelerated in the resonant X-ray scattering era discussed in Section 7, many major challenges await ahead.
Open questions regard crucial aspects, like: i) the fate of charge density waves in the heavily underdoped region of the phase diagram and their interplay with the spin degrees of freedom when approaching the antiferromagnetic phase (see Reference [58]); ii) the role of charge density waves in the marginal Fermi liquid behaviour observed in the anomalous metallic phase above T c around optimal doping, in the formation of the pseudogap, and in the reconstruction of the Fermi surface once (nearly) static order sets in; iii) the role of disorder and inhomogeneity, of the charged dopants residing away, but not too far, from the CuO 2 planes, of commensuration effects when charge density waves are coupled to the lattice, and of lattice strain; iv) the interplay between charge density waves and superconductivity, both as static and as fluctuating orders, including the possibility that quantum critical fluctuations of the charge density wave static order act as mediators of Cooper pairing in the particle-particle channel.
Very recently [87], the marginal Fermi liquid behaviour seems to have found an explanation based on the ubiquitous presence of the broad and nearly temperature independent charge fluctuations observed in resonant X-ray scattering experiments. How this new scenario is accommodated in the attempt to explain marginal-Fermi-liquid-like features of the Raman spectra by means of the narrower dynamical charge density waves, is a further question to be clarified.
As a final remark, although it must be admitted that often different scientist move along very different routes, the various pathways do cross several times. After a while, one always realises that Science is indeed a collective endeavour, and the contributions of other researchers, that appeared obscure at a first glance, become clearer and clearer as time goes by, making every theory, included the one highlighted in this piece of work, the result of ideas that were piecemeal triggered by all other ideas, all experimental findings. Funding: S.C. acknowledge financial support from Sapienza Università di Roma, Project n. RM11715C642E8370.