Skip to main content

Full text of "Local singlets, frustration, and unconventional superconductivity in the organic charge-transfer solids"

See other formats

Local singlets, frustration, and unconventional superconductivity in 

the organic charge-transfer solids * 

R.T. Clay* 1 ' 1 S. Mazumdar b H. Li b 

a Department of Physics and Astronomy and HPC 2 Center for Computational Sciences, Mississippi State University, 

Mississippi State MS 39762 
Department of Physics, University of Arizona Tucson, AZ 85721 


We suggest that superconductivity (SC) in the organic charge transfer solids (CTS) is reached from a Bond-Charge 
Density Wave (BCDW). We discuss an effective model for the BCDW to SC transition, an attractive U extended 
Hubbard Hamiltonian with repulsive nearest neighbor interaction V. We discuss experimental consequences of the 
theory for different classes of CTS superconductors as well as related inorganic materials. 

1. Introduction 

After thirty years of experimental and theoretical 
effort, the mechanism of superconductivity (SC) in 
the organic charge-transfer solids (CTS) is still elu- 
sive. Electron-electron (e-e) interactions, electron- 
phonon (e-p) interactions, and lattice frustration 
all play key roles in the unusual insulating phases 
found in the CTS, but it is not clear how these ap- 
parently very different effects come together in the 
SC state. In this paper we outline a unified theo- 
retical approach to unconventional SC in the entire 
family of organic CTS, and argue that all three of 
the above interactions are essential for SC. Unlike 
existing theories of SC in the CTS, our work can 
potentially explain the pairing mechanism in all 
j -filled molecular superconductors, spanning from 
the quasi-one-dimensional (ID) (TMTTF) 2 X to the 
nearly isotropic two-dimensional (2D) k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 and EtMe 3 Z[Pd(dmit) 2 ] 2 . The 
theory may additionally explain SC found in inor- 
ganic j -filled superconductors. 

* Supported by DOE Grant No. DE-FG02-06ER46315 and 
NSF Grant No. DMR-0705163. 
1 email: 

SC in many CTS, for example in k-(ET) 2 X is 
proximate to antiferromagnetism (AFM). Further, 
these 2D CTS are strongly dimerized suggesting the 
possibility of describing the material by an effective 
^-filled band (one carrier per dimer). It has con- 
sequently been proposed that the i-filled band re- 
pulsive Hubbard model on an anisotropic triangular 
lattice can explain the superconducting behavior of 
the 2D CTS [1,2]. 

The i-filled band triangular lattice Hubbard 
model believed to show SC has two parameters, 
the onsite Coulomb interaction, U, and the lat- 
tice anisotropy, t'/t. The limit t' = corresponds 
to the square-lattice Hubbard model, which is an 
AFM insulator for any U > 0. With increasing t' 
the AFM state is destroyed giving a paramagnetic 
metal (PM). Mean field and variational calculations 
find a narrow d-wave SC region in between AFM 
and PM phases [1,2]. 

A necessary condition for SC mediated by antifer- 
romagnetic fluctuations is that U must enhance the 
superconducting pair-pair correlations. We have cal- 
culated exact superconducting pair-pair correlation 
functions for the i -filled band triangular Hubbard 
model, and find instead that the strength of pair- 
ing correlations always decays monotonically with 

Preprint submitted to Physica B 

9 September 2008 

(a) • = "fc=0---O=«iii* 

Fig. 1. (a) The ...1100... ID BCDW state. Double, double- 
dashed and single-dashed lines indicate strong, weak and 
weakest bonds, respectively, (b) The zigzag ladder BCDW 
state. Double, single and dashed lines again indicate strong, 
weak and weakest bonds, (c) 2D BCDW (Valence-Bond 
Solid (VBS)) state in the CTS. Filled (open) ellipses arc 
charge-rich (charge-poor) molecules. Circled pairs of sites 
constitute a single site in the effective negative- U model in 
(d), where the filled (open) circles are doubly occupied (va- 
cant) sites, respectively. 

increasing U [3] . These results show that there is no 
SC within this model. 

2. Effective model for BCDW/SC transition 

SC in CTS is proximate to not only AFM but to 
many other exotic insulating states such as charge 
ordering (CO) and spin-Peierls (SP) states. We ar- 
gue that to understand the SC in the CTS, one 
must first understand these insulating states. This 
requires returning to the full Hamiltonian with the 
correct band filling (^-filled or a carrier density of 
n = i). For the ID systems our Hamiltonian con- 
tains inter-site e-p coupling, intra-site e-p coupling, 
and Hubbard and extended Hubbard e-e interac- 
tions, H = i?i n ter + ^intra + H ee : 

Winter = ~t ^ [1 + a(Xj - Xi)} (c] )0 .Cj ;(7 + H.C.) (1) 


-ffintra = 9 ^ + + ^intra ^ vf 

i i 

H cc = *7 E n »>l n i,i + ^ X/ UiUj 

i (ij) 

c\ a creates an electron of spin a on site i, x% is the po- 
sition coordinate of molecule i, and 7ij ;Cr = c\ a Ci_ a , 
n-i = n, ; -f + rij^. Vi is the coordinate of the intra- 
molecular vibration on site i, and K int er and K int ra 
are the spring constant for inter-site and intra-site 
phonons. The density of carriers is n = |. 

In a series of works, we have investigated this 
Hamiltonian and its 2D counterparts [4-9]. In gen- 
eral two different kinds of charge ordered states are 
possible: (i) Wigner-crystal ..1010.. CO driven by 
the nearest-neighbor interaction V, and (ii) bond- 
charge density-wave (BCDW) ordering with CO 
pattern ..1100.. driven by the cooperation of e-e 
and e-p interactions. In this notation, '1' indicates 
a charge density of 0.5+5, and '0' a charge density 
of 0.5-5. For realistic U/t ~ 6 - 8 and V/t ~ 1 - 3 
the BCDW state tends to dominate [4-9]. Fig. 1(a)- 
(c) shows the BCDW state for ID, a zig-zag ladder 
lattice, and 2D. 

The BCDW in all cases consists of nearest- 
neighbor singlet bipolarons. It gains energy from 
the exchange interaction between electrons in the 
pairs and also from e-p interaction via the lattice 
distortion. Singlet formation within the BCDW 
pairs leads to a non-magnetic ground state in the 
ID SP state [5], the zigzag ladder [8], as well as the 
2D BCDW [9]. Such nonmagnetic insulating states 
coexisting with CO are found in many CTS. 

We map the nearest-neighbor pairs of singly oc- 
cupied (vacant) sites into double occupancies (va- 
cancies) as shown in Fig. 1 (d) . The Hamiltonian de- 
scribing this effective model now has a negative U 
and a positive V [10]: 

H eS = -t (4,^ + H.c.) (2) 
~ E ( C L<V + H.c.) 


- \U\ y^rti,T"i,j. + V ^ njUj + V'^^n k ni 
{ij} [hi] 

Eq. 2 describes the effective lattice of Fig. 1(d), with 
hopping t along the x and y directions (nearest- 
neighbor bonds (ij)) and t' along the x + y direction. 
Both V and V' nearest-neighbor Coulomb interac- 
tions are included. 

Using exact diagonalizations we have shown that 
the ground state of Eq. 2 is charge ordered for t' = 
with the checkerboard CO of Fig. 1(d). As t' in- 
creases there is a transition from CO to SC [10]. 
Fig. 2 shows the phase diagram as a function of t'- 






1 1 1 1 

i r 1 

1 1 


CO y 

— | \+- 1 h 


1 1 













Fig. 2. Phase diagram of the effective model Eq. 2 from ex- 
act calculations on a 4x4 cluster, (a) |f |-t' phase diagram 
(V = V = 1) (b) V-t' phase diagram (|C|=4). The transi- 
tion is continuous (thick line) for small V. 

U and t'-V. The SC pairing in the effective model 
is of the onsite singlet form. We note however that 
in terms of the original lattice, pairing is primarily 
nearest-neighbor, and a full calculation of the origi- 
nal U > Hamiltonian would likely lead to further 
pair structure such as nodes in the gap function. 
Furthermore, modifications of the simple square+f' 
lattice structure shown in Fig. 1(d) are needed to 
describe the lattices of specific CTS, for example 
the more anisotropic Bechgaard salt superconduc- 
tors. We believe however that this effective model 
can describe the BCDW/SC transition found across 
the various CTS families. 

3. Discussion: CTS materials 

Quasi-W CTS (TMTSF) 2 X: We believe that our 
theory gives the correct insight into understanding 
the SC state in all the ^-filled CTS superconduc- 
tors, including both quasi- ID and 2D CTS. Even for 
the quasi-lD materials triangular lattice structure is 
present in the interchain hopping integrals [11]. The 
addition of intrinsic dimerization to Eq. 1 acts to 
further stabilize the ..1100.. BCDW over the ..1010.. 
CO state [12]. The role of pressure is to increase in- 
terchain coupling and the effective frustration, giv- 
ing mobility to the insulating pairs in the BCDW. 
In the quasi-lD (TMTCF) 2 X (C=Se,S), we have 
shown that while the high temperature CO found in 

this series can have the CO pattern ..1010.., the CO 
pattern found in the SP state and coexisting with 
the SDW adjacent to SC is the ..1100.. BCDW [5, 6]. 

Very large upper critical fields are found in 
(TMTSF) 2 X [13] and other CTS SC [14], along with 
an upward curvature of H c2 as a function of temper- 
ature. These features are a common characteristic 
of superconductivity with short-range local pairing 
and would be expected within our theory [15]. 

k-(ET) 2 X: While many authors have suggested that 
the mechanism of SC in the 2D CTS is linked to 
the presence of the AFM state, AFM adjacent to 
SC is perhaps more unusual than common in the 
2D CTS. For example, within the k-(ET) 2 X series, 
X=Cu 2 (CN) 3 is a superconductor where the insulat- 
ing state is a spin liquid with no apparent spin order- 
ing to millikclvin temperatures [16, 17] . As discussed 
below, there are many examples of CO/SC transi- 
tions within the CTS; this suggests that the AFM 
and BCDW states are quite close in energy, and that 
under pressure the BCDW can replace AFM order. 

Precision dilatometry measurements on n- 
(ET) 2 Cu[N(CN) 2 ]Br find strong lattice effects near 
the AFM/SC transition [18]. Magneto-optical mea- 
surements on k-(ET) 2 Cu(SCN) 2 find changes in 
the vibrational spectra between SC and normal 
states, suggesting that intramolecular vibrations 
are involved in the superconducting transition [19]. 
These experiments suggest strongly the importance 
of lattice degrees of freedom in the SC transition 
even when AFM is present, in agreement with our 

Dmit salts: EtMc3Z[Pd(dmit) 2 ] 2 has a charge/bond 
ordered insulating state that has been described as 
valence-bond solid (VBS) ordering [20, 21]. We point 
out that the charge ordering found in the VBS is 
the same as in our ..1100.. BCDW, see Fig. 3 in ref- 
erence [20]. Under pressure, the VBS state becomes 

j3-(meso-DMBED T- TTF) 2 PF 6 : [3- (meso-DMBEDT- 
TTF) 2 PF 6 is another 2D CTS that has CO and 
undergoes a CO/SC transition under pressure [22]. 
The CO state has been described as a "checker- 
board" pattern, implying that the local CO is 
..1010.. driven by the V interaction. We point out 
however, that the checkerboard here is in terms of 
dimers, and the CO pattern again corresponds to 
the BCDW shown in Fig. 1(c). 


4. Relationship to other theories 

Our theory has strong similarities with RVB and 
bipolaron models of SC simultaneously. As in the 
RVB model, pairing is (partly) driven by AFM cor- 
relations also in our model. The key difference in our 
model is the filling: at density n = | because of the 
fewer number of neighbors it is much favorable en- 
ergetically to form a nearest-neighbor singlet state. 
Furthermore, at n = i e-e and e-p interactions tend 
to act cooperatively [4-9] . 

Traditional bipolaronic SC is based on over- 
screening of the e-e interaction by strong e-p in- 
teractions [15]. This leads to an unphysically large 
effective mass of the paired electrons. In contrast, 
because of the co-operation between the AFM 
and e-p interactions, no overscreening is required 
in our theory. Furthermore, it has recently been 
shown that bipolarons are especially mobile on tri- 
angular lattices [23]. Both n = i and frustration 
are hence essential ingredients of the insulator-to- 
superconductor transition in the CTS. 

5. Application to other materials 

We believe that our theory may be relevant to 
other ^-filled superconductors. Several of these ma- 
terials were known before superconducting CTS. 
LiTi 2 04 with i-filled Ti bands is one of the first 
oxide superconductors discovered with at T c of 12 
K [24]. Because of the presence of frustration in 
the spinel lattice structure, a resonating valence 
bond (RVB) ground state has been suggested [25]. 
However, e-p coupling also appears to be strong in 
this material, and it is still unclear of the relative 
importance of e-e and e-p interactions [26] . 

Very complicated CO ordering occurs in other 
^-filled spinels and is in fact a combination of or- 
bital and charge ordering. In CuIr 2 S4, CO and a 
nonmagnetic (spin singlet) ground state are found. 
This state has been explained as an orbitally driven 
Peierls state, with the pattern Ir 3 +-Ir 3+ -Ir 4+ -Ir 4+ 
[27]. This is again the ..1100.. BCDW CO pattern. 

A CO/SC transition was found under pressure in 
/3-Nao.33V 2 5 [28]. While the stoichiometry of 0- 
Nao.33V 2 05 gives filling of vanadium sites on av- 
erage, there occur here three inequivalent vanadium 
sites with unequal electron distribution. Experiment 
indicates ^-filled zig-zag chains as in the quasi-lD 
CTS [29]. The SC state in this material is extremely 
sensitive to the concentration of Na, with the SC 

transition disappearing for a tiny change in Na con- 
centration (a;=0.32 rather than x=0.33) [30]. This 
strongly suggests that the CO/SC transition here is 
a property of the exact band filling of the system, 
rather than being driven by doping where a broad 
peak around the optimal x would be expected. 


[I] B. J. Powell, R. H. McKcnzic, Phys. Rev. Lett. 94 (2005) 

[2] B. Kyung, A. M. S. Tremblay, Phys. Rev. Lett. 97 (2006) 

[3] R. T. Clay, H. Li, S. Mazumdar, preprint 

arxiv:0805.0590 (2008). 
[4] R. T. Clay, S. Mazumdar, D. K. Campbell, Phys. Rev. 

B 67 (2003) 115121. 
[5] R. T. Clay, R. P. Hardikar, S. Mazumdar, Phys. Rev. B 

76 (2007) 205118. 
[6] S. Mazumdar, S. Ramasesha, R. T. Clay, D. K. 

Campbell, Phys. Rev. Lett. 82 (1999) 1522. 
[7] S. Mazumdar, R. T. Clay, D. K. Campbell, Phys. Rev. 

B 62 (2000) 13400. 
[8] R. T. Clay, S. Mazumdar, Phys. Rev. Lett. 94 (2005) 


[9] R. T. Clay, S. Mazumdar, D. K. Campbell, J. Phys. Soc. 

Jpn. 71 (2002) 1816. 
[10] S. Mazumdar, R. T. Clay, Phys. Rev. B 77 (2008) 


[II] D. L. Pevelen, et al., Eur. Phys. J. B 19 (2001) 363. 
[12] Y. Shibata, S. Nishimoto, Y. Ohta, Phys. Rev. B 64 

(2001) 235107. 

[13] I. J. Lee, M. J. Naughton, G. M. Danner, P. M. Chaikin, 

Phys. Rev. Lett. 78 (1997) 3555. 
[14] F. Zuo, et al., Phys. Rev. B 61 (2000) 750. 
[15] R. Micnas, J. Ranninger, S. Robaszkiewicz, Rev. Mod. 

Phys. 62 (1990) 113. 
[16] Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, 

G. Saito, Phys. Rev. Lett. 91 (2003) 107001. 
[17] Kurosaki, et al., Phys. Rev. Lett. 95 (2005) 17701. 
[18] M. de Souza, et al., Phys. Rev. Lett. 99 (2007) 037003. 
[19] I. Olcjniczak, et al., Phys. Rev. B 67 (2003) 174502. 
[20] M. Tamura, A. Nakao, R. Kato, J. Phys. Soc. Jpn. 75 

(2006) 093701. 

[21] Y. Shimizu, H. Akimoto, H. Tsujii, A. Tajima, R. Kato, 

Phys. Rev. Lett. 99 (2007) 256403. 
[22] S. Kimura, et al., J. Am. Chcm. Soc. 128 (2006) 1456. 
[23] J. P. Hague, P. E. Kornilovitch, J. H. Samson, A. S. 

Alexandrov, Phys. Rev. Lett. 98 (2007) 037002. 
[24] D. C. Johnston, J. Low. Temp. Phys. 25 (1976) 145. 
[25] S. Satpathy, R. M. Martin, Phys. Rev. B 36 (1987) 7269. 
[26] C. P. Sun, et al., Phys. Rev. B 70 (2004) 054519. 
[27] D. I. Khomskii, T. Mizokawa, Phys. Rev. Lett. 94 (2005) 


[28] T. Yamauchi, Y. Ueda, N. Mori, Phys. Rev. Lett. 89 

(2002) 057002. 

[29] K. Okazaki, et al., Phys. Rev. B 69 (2004) 140506(R). 
[30] T. Yamauchi, Y. Ueda, Phys. Rev. B 77 (2008) 104529.