Absence of superconductivity in the Hubbard dimer model for -(BEDT-TTF)2X

In the most studied family of organic superconductors κ-(BEDT-TTF)2X, the BEDT-TTF molecules that make up the conducting planes are coupled as dimers. For some anions X, an antiferromagnetic insulator is found at low temperatures adjacent to superconductivity. With an average of one hole carrier per dimer, the BEDT-TTF band is effectively 2 -filled. Numerous theories have suggested that fluctuations of the magnetic order can drive superconducting pairing in these models, even as direct calculations of superconducting pairing in monomer 2 -filled band models find no superconductivity. Here we present accurate zero-temperature Density Matrix Renormalization Group (DMRG) calculations of a dimerized lattice with one hole per dimer. While we do find an antiferromagnetic state in our results, we find no evidence for superconducting pairing. This further demonstrates that magnetic fluctuations in the effective 2 -filled band approach do not drive superconductivity in these and related materials.


Introduction
Correlated-electron superconductivity (SC) continues to be one of the most challenging problems in condensed matter physics. Even as the large majority of investigators in this research area have focused on cuprates and iron-based superconductors (and more recently, ruthenates), it is now widely recognized that there exist many other unconventional or correlated-electron superconductors, although with lower superconducting critical temperature T c [1][2][3][4][5][6][7][8][9]. Superconducting organic charge-transfer (CT) solids hold a special place in this context, as SC in CT solids was discovered significantly prior to the discovery of the phenomenon in the cuprates [10]. In spite of intensive investigations over the past four decades, however, the mechanism of SC in CT solids remains a mystery, and organic SC continues to be an active research area, as illustrated by the appearance of many review articles on the subject over the past decade [11][12][13][14][15]. In the present paper we report the results of our latest calculations that have probed a key question in the field of organic SC.
The common theme among many (though not all, see below) correlated-electron superconductors is the proximity of SC to spin-density wave (SDW) or antiferromagnetism (AFM) . This property is shared by organic superconductors based on the cation TMTSF as well as some κ-phase BEDT-TTF compounds. Unlike in the cuprates, though, transition from the insulating magnetic phase to the superconducting phase is driven by the application of moderate pressure as opposed to doping. Both TMTSF and BEDT-TTF based superconductors (hereafter (TMTSF) 2 X and (BEDT-TTF) 2 X) are characterized by cation:anion ratio of 2:1, with closed-shell monovalent anions X − , the implication being that each individual cation monomer has charge + 1 2 . In both families, the cation lattice is highly dimerized, such that the overall charge on a dimer of cations is exactly +1. This has led to various effective half-filled band theories of organic CT superconductors, within which the dimer of cations is considered as a single site in a lattice described within anisotropic triangular lattice Hubbard models (we emphasize that these theories are therefore not applicable to nondimerized (BEDT-TTF) 2 X families, such as those belonging to the θ-family, which also exhibit SC; we revisit this issue later). Within meanfield and dynamic mean field theory (DMFT) calculations on the half-filled band anisotropic triangular lattice Hubbard model, pressure increases the lattice frustration, causing thereby a bandwidth-driven transition from the AFM insulator to the superconducting state (or in the case of a lattice that is already strongly frustrated, as in κ-(BEDT-TTF) 2 Cu 2 (CN) 3 , from a quantum spin liquid (QSL) to SC) [16][17][18][19][20][21][22][23][24][25][26][27][28]. Direct numerical calculations of superconducting pair-pair correlations within the 1 2 -filled Hubbard band on triangular lattices have, however, consistently found absence of SC [29][30][31][32][33]. Taken together with the very recent convincing demonstration of the absence of SC within weakly doped the square lattice two-dimensional (2D) Hubbard Hamiltonian [34], numerical calculations then cast serious doubt on the existing approximate spin fluctuation mechanisms of AFM-to-SC transition.
Conceivably, one possible origin of disagreement between the mean-field theories of organic SC and numerical calculations is that retention of the explicit dimeric structure of the Hubbard lattice sites is essential for superconducting correlations to dominate. That is, even as intra-dimer charge fluctuations are precluded within the effective 1 2 -filled band models (since that would make the effective 1 2 -filled description invalid and can even destroy the AFM in the weakly frustrated region of the Hamiltonian), the superconducting singlets should be defined such that they span over a pair of dimers, as opposed to a pair of effective single sites. Needless to say such a numerical calculation is far more involved than where dimers are replaced with single sites, as the required lattice sizes needed are twice that in the previous calculations [29,30,33]. We report precisely such a calculation here, which again finds absence of SC in such a dimer lattice.

Methods
We consider the Hubbard Hamiltonian, In Eq. 1 c † i,σ creates an electron of spin σ on molecule i and n i,σ = c † i,σ c i,σ . U is the onsite Hubbard interaction. It is important to note that this U is the Coulomb repulsion between carriers on each monomer site and is not the same as the effective dimer U d . The lattice we consider is a rectangular lattice of dimers (see Fig. 1). In this lattice there are three hopping integrals t ij : t d is the intra-dimer hopping, t the inter-dimer hopping, and t a frustrating hopping (see Fig. 1). We take the system to be 1 4 filled with an average electron density of 1 2 per monomer site. This lattice is expected to show Néel AFM order in the unfrustrated limit of small t /t. To solve for the ground state of Eq. 1 we employ the Density Matrix Renormalization Group (DMRG) method [35] using the ITensor library [36] with real-space parallelization [37]. DMRG can calculate essentially numerically exact correlation functions in quasi-onedimensional systems. However, its primary limitation is an exponential scaling in the transverse size of systems that can be studied. Here we consider lattices of dimension 4 sites in y and up to 16 sites in x (see Fig. 1), corresponding to an 8×4 lattice in terms of dimers. We choose t d =1.5 and t=0.5 such that the average hopping along the x axis is 1.0; our choice of t d ∼ 3t is comparable to the dimerization found in the CT solid superconductors [14]. We take boundary conditions that are open (periodic) in x (y). We used a DMRG bond dimension m of up to 15,000 with minimum truncation error of 10 −8 -10 −7 ; all results we show are extrapolated to zero truncation error.

Results
We measure dimer spin-spin and pair-pair correlation functions. The dimer spin operator is n d i,σ = n i 1 ,σ + n i 2 ,σ , where the sites i 1 and i 2 make up the dimer i. We calculate the dimer spin-spin correlation Fig. 2 we plot S d z (1, j) (−1) j where dimer 1 is on the first chain and dimer j is on the neighboring chain for t = 0.2 and 0.4. As shown in Fig. 2(a), for t = 0.2 and U 8 this quantity is positive for all spin correlations, indicating the presence of Néel AFM order. For t = 0.4 we do not find Néel AFM order up to at least U = 12. These results show that Néel AFM order is only present in this system when the frustration is low, and for realistic U values AFM will not be present when t /t 0.4.
To define superconducting pair-pair correlations in the effective 1 2 -filled representation, we first define an operator creating an effective particle on dimer i with equal monomer populations, as required within effective 1 2 -filled band theories, The ∆ † i,i+δ operator creates a singlet pair between nearest-neighbor (n.n.) dimers separated by δ: Note that Eq. 3 involves four different lattice sites. We consider two kinds of n.n. pairs with δ taken as the n.n. distance between dimers in the x and y directions, In the effective dimer model pairing is expected to have d-wave symmetry. Because of the small transverse dimensions of our lattice we do not define full d x 2 −y 2 pairs involving four n.n. singlets, but instead calculate only single singlet-singlet correlations. To check if the pairing has the expected d-wave sign structure we calculate two different correlation functions, P || (r) ≡ P(r) δ=x,δ =x and P ⊥ (r) ≡ P(r) δ=x,δ =y . P || (r) corresponds to correlations between two singlets both oriented along x, while P ⊥ (r) corresponds to one singlet oriented along x and one along y. For d x 2 −y 2 pairing P ⊥ (r) should be negative. In a 2D superconductor long-range order should be present in P(r) at zero temperature. Because of the quasi-one-dimensional nature of cylindrical geometry of the lattice that we can solve in DMRG long-range order in P(r) is not possible even in the limit x → ∞. If such a quasi-one-dimensional system showed a tendency to SC however, the system would behave as a Luther-Emery liquid [38] with P(r) decaying as a power law r −α with α < 1. In addition to long-range order in P(r), if SC in the system is driven by Coulomb interactions, one expects that P(r, U) > P(r, U = 0) [29].
In Fig. 3 we plot P || (r) as a function of r. In Fig. 3(a) the indices i and j in Eq. 4 are taken on the same chain, in Fig. 3(b) they are on neighboring chains, and in Fig. 3(c) they are on   Fig. 3 except that the pair-pair correlation P ⊥ (r) is for the perpendicular orientation of pairs. Note that in (d)-(f),P ⊥ is negative, with decreasing magnitude as U increases.P ⊥ in (e) is strongly affected by shorter-range correlations which lead to the unusual U dependence (see text).
second neighbor chains. In each panel r is the center-to-center distance between pairs in units of the dimer-dimer spacing. In Fig. 3(a) we also plot the function r −1 . While the distances we have on the lattice are limited, and we cannot determine if P(r) decays with r as a power law or as an exponential, we find that on each chain P(r) decays significantly faster than r −1 .
As seen in Fig. 3(a) and (c), only at very short range (r ≤ 2) we find that P(r, U) is enhanced very slightly over its value at U = 0. This same effect is seen in calculations in the monomer 1 2 -filled Hubbard model on an anisotropic triangular lattice [29,30]. As noted in [39], at r = 0 the pair-pair correlation can be written as a linear combination of charge-charge and spin-spin correlations, which are enhanced by U [39]. To measure the effect of correlations on pairing we define the average long-range pair-pair correlation [40], In Eq. 5 N p is the number of terms in the sum, and the sum is over r greater than 2 in dimerdimer units but less than a maximum distance r max . The upper cutoff is necessary because of the open boundary conditions along x in our lattice and is chosen so that the maximum pair-pair distance along the chain is 4 (one half of the lattice length of 8 dimers). Fig. 3(d)-(f) showP || as a function of U. In each case the magnitude ofP || decreases with U apart from an insignificant increase at very small U. A large decrease inP || also takes place for U 8 upon entering the AFM insulating phase. A similar decrease is seen in calculations within the monomer 1 2 -filled band [29,30]. Fig. 4(a)-(c) shows perpendicular correlations P ⊥ (r) as a function of r. Like P || , we find that P ⊥ along the chain direction decays faster than r −1 . These are negative in many cases ( Fig. 4(a) and (c)), which would be consistent with a d x 2 −y 2 type pairing symmetry. Fig. 4(d) shows a clear decrease in magnitude of P ⊥ (r) with increasing U. The U dependence of Figs. 4(e)-(f) are harder to interpret. We believe that the unusual U dependence here is due to shorter range correlations, for example as in Fig. 4(b) where P ⊥ (r) changes discontinuously at small U at r ∼ 2 [41]. If distances r < 3 are excluded, we again find a continuous decrease of P ⊥ with increasing U.  Fig. 6 show the pair-pair correlations at t = 0.6. Here AFM order is absent over the whole range of U we studied. We again find that P || (r) decreases faster than r −1 (see Fig. 5(a)). As at smaller t we again find that some pair-pair correlations are enhanced by U at short range (r ≤ 2). However, here bothP || andP ⊥ decrease continuously with U, again with the exception of one point at shorter r where there is a discontinuous change at small U (see Fig.6(b) and (e)). Summarizing our data, while we cannot show the presence of long range AFM order from our single lattice calculation, we do find strong AFM correlations for weak frustration and moderately strong U. We find no evidence that within the effective 1 2 -filled picture that the superconducting pair-pair correlations are enhanced over their value for non-interacting electrons.

Discussion
The conclusions of our work are, (i) SC is absent within the model, and (ii) the presence or absence of AFM order is not related to the tendency towards SC.
In κ-(BEDT-TTF) 2 X with X = Cu[N(CN) 2 ]Cl where AFM order is found under ambient pressure, estimates for t /t vary between 0.44 and 0.53 [14]. Similarly, no AFM order is found for X = Cu 2 (CN) 3 where 0.80 t /t 0.99 [14]. In our present results we find AFM order at t = 0.2 (t /t = 0.4) but no AFM order at t /t = 0.8, which is roughly consistent with what is observed experimentally. SC is however found in both X = Cu[N(CN) 2 ]Cl and X = Cu 2 (CN) 3 , under pressures of 0.03 and 0.6 GPa, respectively. Moreover, there is no obvious connection between the superconducting T c in κ-(BEDT-TTF) 2 X and the value of t /t [14]. For example SC is also found for t /t > 1 (t /t = 1.3 in X = CF 3 SO 3 [42]). This suggests that SC in κ-(BEDT-TTF) 2 X is not related to the proximity to an AFM or QSL phase. This is also consistent with our present results, where we find little difference in the pair-pair correlations between t = 0.2 where AFM is present and t = 0.6 where AFM is absent.
Our results on the absence of SC in the effective 1 2 -filled Hubbard dimer model is in complete agreement with our previous results on the 1 2 -filled Hubbard monomer model [29,30,33]. Independent of whether one considers the dimers of BEDT-TTF molecules as effective single sites with charge occupancy +1 (previous calculations) or as dimers with monomer charge occupancies of 1 2 each (present calculations), the superconducting pair correlations are found to be suppressed by Hubbard U continuously from U = 0. This fundamental observation was arrived at from exact diagonalization [29,33], Path Integral Renormalization Group [30] and now DMRG on completely different lattices, thereby confirming that the result is not an artifact of the choice of a specific lattice or a particular computational technique. Taken together, these results very clearly indicate that spin fluctuations are not driving organic SC. With hindsight, this conclusion should not be surprising: as already pointed out in the above, AFM occurs in only three κ-(BEDT-TTF) 2 X, X = Cu[N(CN) 2 Cl, deuterated X = Cu[N(CN) 2 Br and CF 3 SO 3 . Whether or not X = Cu 2 (CN) 3 and Ag 2 (CN) 3 are true QSLs continue to be debated, and both claims of a CO phase below 6 K [43] as well its absence [44] exist in the literature. What is perhaps very relevant in this context that transition from CO to SC is far more common in CT solids, including in the β, β , β and θ-families [14]. Even among the κ-phase materials, CO has been detected in X = Hg(SCN) 2 Cl [45][46][47] while in the related material X = Hg(SCN) 2 Br there occurs a dipolar liquid phase [48] the formation of which is related to the mechanism of CO formation [49,50]. While these latter materials are also strongly correlated, CO requires an an effective 1 4 -filled (or 3 4 -filled description of CT solids, within which intra-dimer charge dimer degrees of freedom are explicitly taken into account [14,[51][52][53][54]. The charge-charge correlations in the CO phase in these systems are different from that in Wigner crystals. The CO phase is a paired-electron crystal, with nearest neighbor spin-singlet charge-rich molecules separated by pairs of charge-poor molecules. Pressure increases frustration, causing the spin singlets to be mobile. Calculations on an anisotropic triangular lattice [40] as well as on κ-lattices [41,55] have indeed found that superconducting pair correlations are enhanced by the Hubbard U uniquely at 1 4 -filling. Additionally, electron-electron and electron-bond phonon interactions act co-operatively at this bandfilling to enhance superconducting pair correlations, even at they compete with one another at all other fillings [56]. These results are encouraging and suggest further research along this direction.