## Abstract

Unconventional superconductivity in molecular conductors is observed at the border of metal-insulator transitions in correlated electrons under the influence of geometrical frustration. The symmetry as well as the mechanism of the superconductivity (SC) is highly controversial. To address this issue, we theoretically explore the electronic properties of carrier-doped molecular Mott system *κ*-(BEDT-TTF)_{2}X. We find significant electron-hole doping asymmetry in the phase diagram where antiferromagnetic (AF) spin order, different patterns of charge order, and SC compete with each other. Hole-doping stabilizes AF phase and promotes SC with *d*_{xy}-wave symmetry, which has similarities with high-*T*_{c} cuprates. In contrast, in the electron-doped side, geometrical frustration destabilizes the AF phase and the enhanced charge correlation induces another SC with extended-*s* + \(d_{x^2 - y^2}\)wave symmetry. Our results disclose the mechanism of each phase appearing in filling-control molecular Mott systems, and elucidate how physics of different strongly-correlated electrons are connected, namely, molecular conductors and high-*T*_{c} cuprates.

## Introduction

Understanding the intimate correlation among metal-insulator (MI) transition, magnetism, and superconductivity (SC) is one of the most challenging issues in modern condensed matter physics. The most well-studied example is the high-*T*_{c} cuprates, where SC is observed when mobile carriers are doped into the parent antiferromagnetic (AF) Mott insulators^{1,2,3}. A general understanding there, supported by various experiments and theories, is that strong AF spin fluctuation mediates the *d*-wave SC that appears through the filling-control Mott MI transition generating mobile charge carriers. However, can we export this mechanism to other strongly correlated materials? To address this question, it is crucial to make a comparison among different classes of materials. In this respect, heavy fermion compounds and molecular conductors provide such opportunities^{4,5,6,7,8}.

The family of quasi two-dimensional molecular conductors *κ*-(ET)_{2}*X* (ET = BEDT-TTF, and *X* takes different monovalent anions^{9}) is in fact compared often with the cuprates^{10}. They indeed have common factors: simple quasi-two-dimensional electronic structure to begin with in the non-interacting limit and the Mott MI transition and SC closely related with each other. However, there are important differences: First, in *κ*-(ET)_{2}*X*, SC appears through the bandwidth-control Mott transition; the carrier density is usually unchanged but the pressure (either physically or chemically) is the controlling factor. Although the variation of the carrier density is necessary for direct comparisons, it has not been realized in *κ*-(ET)_{2}*X* for a long time due to experimental difficulties. Second, while the cuprates are basically governed by the physics nearby 1/2-filling, *κ*-(ET)_{2}*X* is a 3/4-filled system. The similarity enters when the so-called dimer approximation is applied in the latter^{11}, resulting in the effective 1/2-filled system (dimer model). Although the dimer model has been extensively studied using various theoretical methods^{12,13,14,15,16,17,18,19,20,21,22,23,24,25}, the validity of the dimer approximation itself is recently reexamined^{26,27,28,29,30,31,32}. Especially, the importance of the intradimer charge degree of freedom^{33,34} and intersite Coulomb interactions^{27,29}, which are discarded in the dimer approximation, has attracted much attention because of recent experimental suggestions^{35,36,37,38}. Third, in *κ*-(ET)_{2}*X*, SC is observed not only next to the AF insulators but also to nonmagnetic (candidate of gapless spin-liquid) insulators^{8,39,40,41}. The strong influence of geometrical frustration owing to the anisotropic triangular arrangement of dimers is present in this family.

Recently, carrier doping has been realized either chemically in *κ*-(ET)_{4}Hg_{3−δ}*Y*_{8} (*Y* = Br or Cl)^{42,43,44,45} or in *κ*-(ET)_{2}Cu[N(CN)_{2}]Cl (*κ*-Cl) by using electric-double-layer transistor (EDLT) technique^{46,47}, revealing intriguing phenomena such as a dome-shaped SC region, anomalous metallic behaviors, and significant electron-hole doping asymmetry, which are all reminiscent of the high-*T*_{c} cuprates. Therefore, *κ*-type ET systems can now provide a unique playground of both filling-control and bandwidth-control Mott transitions with SC phases nearby, for which a unified theoretical understanding is highly desired.

In this paper, we theoretically study the ground-state properties of *κ*-(ET)_{2}*X* varying the carrier number from 3/4-filling, in order to elucidate the electronic phases appearing near the Mott transition in this system, especially SC, and to investigate their stabilities beyond mean field treatments. The intradimer charge degree of freedom and intersite Coulomb interactions are explicitly considered. We find that the ground-state phase diagram shows significant electron-hole asymmetry in the stability of AF phase and in terms of competing two types of SC. While in the hole-doped side *d*_{xy}-wave SC is favored by the AF spin fluctuation as in the high-*T*_{c} cuprates, the electron doping highlights the geometrical frustration and the charge degree of freedom, which are unique in *κ*-(ET)_{2}*X*, stabilizing extended-*s* +\(d_{x^2 - y^2}\)-wave SC. The electron-hole doping asymmetry, including the symmetry of SC, is attributed to the degree of frustration that is controlled by carrier doping. This is a conceptually new perspective to *κ*-(ET)_{2}*X*, which can also be applied to other frustrated systems in general. Our results, beyond the usual description based on the 1/2-filled dimer model, thus provide new understanding of how physics of molecular conductors and high-*T*_{c} cuprates are distinct.

## Results

### Model derivation and framework

The electronic properties of molecular conductors are modeled by a simple model where the molecules are replaced by lattice sites^{48}. They are described by the extended Hubbard model (EHM)^{11,49}, a textbook model for studying correlated electrons. The Hamiltonian is given as

where \(c_{i\sigma }^\dagger\) (*c*_{iσ}) is a creation (annihilation) operator of electron at molecular site *i* with spin *σ*(=↑,↓), \(n_{i\sigma } = c_{i\sigma }^\dagger c_{i\sigma }\), and *n*_{i} = *n*_{i↑} + *n*_{i↓}. *U* and *V*_{ij} are on-site and intersite Coulomb repulsions, respectively. <*i*, *j*> denotes pairs of neighboring molecules in the *κ*-type geometry, labeled by *b*_{1}, *b*_{2}, *p*, and *q*, as shown in Fig. 1a.

The tight-binding parameters *t*_{ij} are set for the deuterated *κ*-(ET)_{2}Cu[N(CN)_{2}]Br (*κ*-Br), which locates very close to the MI transition^{50}, and are adopted from a first-principles band calculation as \((t_{b_1},t_{b_2},t_p,t_q) = (196,65,105, - 39)\,{\mathrm{meV}} = (1.0,0.332,0.536, - 0.199)\,t_{b_1}\)^{51}. We set the largest hopping integral \(t_{b_1}\) as the unit of energy. The unit cell is a rectangle with *R*_{x} × 2*R*_{y} and **δ**_{±} = (*δ*_{x}, ±*δ*_{y}) are vectors connecting the centers of molecules facing each other in a dimer (see Fig. 1a). Here, we set (*R*_{x}, *R*_{y},* δ*_{x},* δ*_{y}) = (1.0, 0.7, 0.3, 0.3)*R*_{x} with *R*_{x} as a unit of length^{29}. The non-interacting band structure is shown in Fig. 1b. Among the four energy bands, the upper two bands (bands 1 and 2) contribute to form the Fermi surface (FS). For the undoped case, the electron density per molecular site is *n* = 3/2 (3/4-filling) and it corresponds to the hole density *n*_{hole} = 2 − *n* = 1/2. In this study, we change *n*_{hole} from 1/3 to 2/3 to investigate the doping dependence of the system. The corresponding Fermi energies are indicated in Fig. 1b by dotted lines.

The effect of Coulomb interactions is treated using a variational Monte Carlo (VMC) method^{52,53,54}. The trial wave function considered here is a Gutzwiller–Jastrow type, \(\left| \Psi \right\rangle = P_{{\mathrm{J}}_{\mathrm{c}}}P_{{\mathrm{J}}_{\mathrm{s}}}\left| {\mathrm{\Phi }} \right\rangle\). |Φ> is a one-body part constructed by diagonalizing the one-body Hamiltonian, and \(P_{{\mathrm{J}}_{\mathrm{c}}}\) and \(P_{{\mathrm{J}}_{\mathrm{s}}}\) are charge and spin Jastrow factors, respectively. The explicit form of them are described in Methods. In the following, we show results for 1152 molecular sites (corresponding to *L* = 24, see in Methods), which is large enough to avoid finite size effects.

### Ground-state phase diagram

Figure 2a shows the ground-state phase diagram. The hole density *n*_{hole} = 2 − *n* and the on-site Coulomb interaction \(U{\mathrm{/}}t_{b_1}\) are varied as parameters, while the ratio between *U* and the largest intersite Coulomb interaction \(V_{b_1}\) is fixed at \(V_{b_1}{\mathrm{/}}U = 0.50\). The other intersite Coulomb interactions are set as \((V_{b_2},V_p,V_q) = (0.56,0.66,0.58)\,V_{b_1}\), assuming the 1/*r*-dependence. At *n*_{hole} = 1/2 (undoped case)^{29}, a first-order phase transition occurs, with increasing \(U/t_{b_1}\), from a paramagnetic metal (PM) to a dimer-type AF (DAF) phase in which the spins between dimers order in a staggered way as shown in Fig. 2b. This transition corresponds to the Mott MI transition. As \(U/t_{b_1}\) increases further, there appears a polar charge-ordered (PCO) phase breaking the inversion symmetry^{31,34} with AF spin order, which can avoid the energy loss of \(V_{b_1}\), \(V_{b_2}\), and *V*_{p}, at the expense of the energy loss of *V*_{q} as shown schematically in Fig. 2c. The DAF and PCO phases are insulating at *n*_{hole} = 1/2. They have also been found in previous studies for the 3/4-filled Hubbard models^{11,31,49} and the effective strong coupling models^{33,34}, and are stabilized in the relevant parameter regions for *κ*-(ET)_{2}*X*. Experimentally, the DAF phase is widely observed in *κ*-(ET)_{2}*X* as AF dimer-Mott insulator and the PCO phase is proposed to be related to the dielectric anomaly observed in *κ*-(ET)_{2}Cu_{2}(CN)_{3} (*κ*-CN)^{36} and the insulating phase in *κ*-(ET)_{2}Hg(SCN)_{2}Cl^{55}. Note that SC is a metastable state for \(U/t_{b_1}\) = 7–11.5, lying on each side of the DAF-PCO boundary.

Away from *n*_{hole} = 1/2, significant doping asymmetry is observed and several different phases appear. For the hole-doped side (*n*_{hole} > 1/2), while the PCO phase is rapidly suppressed, the DAF phase is enhanced to a smaller \(U/t_{b_1}\) region toward *n*_{hole} = 2/3. Note that in these phases the system becomes metallic once the doping is finite. Furthermore, a 3-fold charge-ordered (3-fold CO-1) phase appears for larger \(U/t_{b_1}\). The 3-fold CO-1 phase shows charge disproportionation and magnetic order as shown in Fig. 2d; hole-rich sites form a two-dimensional network with AF spin order. This phase is insulating at *n*_{hole} = 2/3 (along the brown line in Fig. 2a) since the electron density fits the commensurability, and metallic for other hole densities because the excess holes can move through the ordered holes.

For the electron-doped side (*n*_{hole} < 1/2), the situation is much different. The PCO and DAF (both become metallic) phases are rapidly suppressed and another CO (3-fold CO-2) phase and a SC phase appear. The pattern of the charge disproportionation is opposite to that in the 3-fold CO-1 phase (hole rich ↔ hole poor) as shown in Fig. 1e; this configuration can fully avoid the intersite Coulomb interactions. Similar to the 3-fold CO-1 phase, the 3-fold CO-2 phase is insulating at *n*_{hole} = 1/3 (along the blue line in Fig. 2a) and metallic for other hole densities. Note that the 3-fold CO-2 phase is stabilized also for *n*_{hole} = 1/2 when \(V_{b_1}/U\) is larger (≥0.55)^{29,31} and is smoothly connected in the parameter space. While the 3-fold CO-1 phase is accompanied by the magnetic order, 3-fold CO-2 is nonmagnetic. We have tried several magnetic ordering patterns that coexist with 3-fold CO-2. However, none of them are stabilized because the CO pattern in the hole-rich sites forms a triangular-like structure and the spin degree of freedom is fully frustrated, and furthermore the distances between hole-rich sites are much longer than the original intermolecular bonds and therefore the effective magnetic exchange couplings are quite small. For smaller \(U/t_{b_1}\), the SC phase is realized by doping, located between the DAF/3-fold CO-2 and PM phases. The symmetry of the SC is the extended-*s* + \(d_{x^2 - y^2}\)-wave type, same with the one shown in our previous study at *n*_{hole} = 1/2^{29}. Details are discussed later.

Although SC does not appear as the ground state in the hole-doped side, we find finite superconducting condensation energy in the phase diagram. Figure 2f shows the region where the condensation energy is finite, ignoring other ordered phases by setting Weiss fields to be zero in |Φ〉. It is possible that the hidden SC phase appears if the DAF phase is destabilized by, e.g., disorder effect associated with doping or phase separation. Therefore, it is worthwhile to study the most favored SC phase even if it is a metastable state. While the *d*_{xy}-wave SC is dominant for most of the hole-doped side, the extended-*s* + \(d_{x^2 - y^2}\)-wave SC is stabilized for the electron-doped side. Namely, the symmetry of SC changes with carrier doping. Note that the charge correlation is greatly enhanced toward regions indicated by gray shade in Fig. 2f. In these regions, the mobility of holes are strongly restricted due to the strong intersite Coulomb interactions, and stable VMC simulations are difficult unless additional Weiss fields that induce long-range CO are introduced in |Φ〉.

### Fermi surface and spin structure factor

The electron-hole doping asymmetry is closely related to the shape of the FS and the interdimer magnetic fluctuations. Figure 3a–c show the non-interacting FS for *n*_{hole} = 2/3, 1/2, and 1/3. As *n*_{hole} increases from 1/2 (hole doping), the FS shifts toward the right and left edges of the first Brillouin zone (M-X line in Fig. 3b). Since the energy gap of the DAF order opens along the Brillouin zone edge^{11}, the DAF order becomes more favored for hole doping. For *n*_{hole} = 2/3, the FS almost touches the Brillouin zone edge, and there the DAF region extends down to \(U{\mathrm{/}}t_{b_1}\sim 2.5\), as shown in Fig. 2a. Note that the Fermi energy is located in the vicinity of van Hove singularity at *n*_{hole} = 2/3 as shown in Fig. 1b. Around this hole density, anomalous behavior such as pseudogap phenomena is naively expected^{46}. In clear contrast, the FS departs from the M-X line for electron doping, consistent with the tendency of the DAF order being rapidly suppressed for *n*_{hole} < 1/2.

Next, Fig. 3d–f show the interdimer spin structure factor defined as,

for *n*_{hole} = 2/3, 1/2, and 1/3. Here, *N*_{dim}(=*L*^{2}) is the total number of dimers and \(M_l^{{\mathrm{dim}}} = (n_{2l - 1 \uparrow } + n_{2l \uparrow }) - (n_{2l - 1 \downarrow } + n_{2l \downarrow })\) is the total spin density within *l*-th dimer formed by molecular sites 2*l* − 1 and 2*l* with the central position **r**_{l}^{56}. Since the dimer centers form the anisotropic triangular lattice^{11}, the corresponding first Brillouin zone is the anisotropic hexagon. As shown in Fig. 3d, *S*^{dim}(**q**) for *n*_{hole} = 2/3 peaks around (0, ±*π*/*R*_{y}), which corresponds to the DAF spin configuration, suggesting that the AF spin fluctuation is enhanced by the Coulomb interactions and thus the DAF order is favored. On the other hand, for *n*_{hole} = 1/2, the peaks appear around six vertices of the Brillouin zone, as shown in Fig. 3e. This implies that the spin structure becomes more triangular-lattice like (frustrated) and the AF spin fluctuation is suppressed as compared with that at *n*_{hole} = 2/3. For *n*_{hole} = 1/3, the peak structures almost diminish (see Fig. 3f), and the DAF order is not stabilized around this hole density.

### Superconducting gap functions

The above mentioned electron-hole asymmetry in the spin and the charge degrees of freedoms are the keys to understand the competition between the two types of SC. The main contribution of the gap function for the *d*_{xy}-wave SC is given as

where *α*(=1, 2) denotes a band index, and \({\mathrm{\Delta }}_i^\alpha\) is the pairing with the *i*-th neighbor dimers in the real space and treated as a variational parameter. We optimize the real space pairing up to 22nd neighbor dimers and find that the overall feature of the gap function is determined within the fourth neighbor, i.e., \({\mathrm{\Delta }}_m^\alpha\) for *m* ≤ 4. The term contaning \({\mathrm{\Delta }}_1^\alpha\) in Eq. (3) gives nodes in the horizontal (along *k*_{x}-axis) and vertical (along *k*_{y}-axis) directions. This is because the two diagonal pairings (orange and blue bars in Fig. 4a) have different sign, giving a *d*_{xy}-type contribution. Therefore, this gap symmetry is referred to as *d*_{xy}-wave^{26}. Note that the terms corresponding to real space pairings parallel to horizontal (along *x*-axis) and vertical (along *y*-axis) directions vanish since they are along the nodal directions. This SC phase has been discussed in analogy with that of high-*T*_{c} cuprates, ascribing the diagonal directions in *κ*-type structure to an approximate square lattice^{10}. As in the case for high-*T*_{c} cuprates, the sign of the gap function on the FS changes four times (see Fig. 4c).

On the other hand, the main contribution of gap function for the extended-*s*+\(d_{x^2 - y^2}\)-wave SC is given as

The first term contaning \({\mathrm{\Delta }}_1^\alpha\) in Eq. (4) (blue bars in Fig. 4b) does not change sign within the first Brilloin zone and becomes zero only along the zone boundary; this term gives an extended-*s*-like contribution. Furthermore, \({\mathrm{\Delta }}_1^\alpha\) changes sign between different bands, namely, \(sgn{\mathrm{\Delta }}_1^1 = - sgn{\mathrm{\Delta }}_1^2\). In this respect, the pairing symmetry can also be referred to as *s*_{±}, similar to that of iron-based SC^{57,58}. The second and third terms containing \({\mathrm{\Delta }}_2^\alpha\) and \({\mathrm{\Delta }}_3^\alpha\) in Eq. (4) (orange and yellow bars in Fig. 4b, respectively) give nodes in the diagonal direction; these terms give a \(d_{x^2 - y^2}\)-like contribution. The gap symmetry is thus referred to as an extended-*s* + \(d_{x^2 - y^2}\)-wave^{26,28,32,59}. The sign changes of the gap function on the FS is shown in Fig. 4d.

The competition between the two types of SC has been already discussed for the undoped case. First, in the effective 1/2-filled dimer model, many early studies inferred the *d*_{xy}-wave type^{12,13,14,16,17,18}. However, its stability over the AF phase has recently been doubted in several numerical studies including the VMC approach^{20,22}. Kuroki et al.^{26} were the first to point out the importance of treating the 3/4-filled model, namely, considering the intradimer charge degree of freedom. In the 3/4-filled Hubbard model, the two types of SC compete and either of them is favored depending on the parameters, especially the degree of dimerization^{26,28,32}. It is only recently that the importance of the intersite Coulomb interaction *V*_{ij} on SC was pointed out^{27,29}. In fact, our previous VMC study for the undoped case found that, although both symmetries are enhanced by the intersite Coulomb interactions, the extended-*s* + \(d_{x^2 - y^2}\)-wave is slightly more favored^{29}.

As shown in Fig. 2f, our calculations find that this competition is released when the mobile carriers are doped into the system. For the hole-doped side, the AF spin fluctuation toward the DAF order is enhanced due to the shape of the FS and the Coulomb interactions, as already seen in Fig. 3, and consistent with the recent calculations based on the dimer model^{46}. Similarly to the high-*T*_{c} cuprates, the AF spin fluctuation mediated SC is then developed with the strong singlet correlation along diagonal bonds shown in Fig. 4a, resulting in the *d*_{xy}-wave symmetry. On the other hand, the AF spin fluctuation is suppressed with electron doping as seen in Fig. 3f and the spin singlet correlation along horizontal and vertical bonds compete with that along diagonal bonds. The extended-*s* + \(d_{x^2 - y^2}\)-wave is eventually favored since all bonds can contribute to the singlet pairing. The carrier doping deforms the shape of the FS and modifies the AF spin fluctuation, thus inducing the change of the symmetry of SC. Although the electron-hole asymmetry shown in Fig. 2a appears similar to that of high-*T*_{c} cuprates at a glance, they are much different especially for the SC phase. In both cases, the van Hove singularity appears only in the hole-doped side, causing the electron-hole asymmetry. However, different physics are delicately involved in *κ*-(ET)_{2}*X*, as we have shown so far, and even the change of the symmetry of SC occurs by carrier doping. This is due to the unique geometrical frustration inherent in the triangular-like lattice structure of *κ*-(ET)_{2}*X*. Furthermore, this asymmetry is expected to be robust for *κ*-(ET)_{2}*X* in general because the van Hove singularity is always located in the hole-doped side for the realistic parameter set of these conductors.

We can show more directly how the spin and charge correlations are correlated to the stability of the SC phases. The interdimer charge structure factor is defined as

where \(n_l^{{\mathrm{dim}}} = (n_{2l - 1 \uparrow } + n_{2l \uparrow }) + (n_{2l - 1 \downarrow } + n_{2l \downarrow })\) is the total charge density within *l*-th dimer formed by molecular sites 2*l* − 1 and 2*l* with the central position **r**_{l}^{56}. Figure 5a–d show *S*^{dim}(**q**) and *N*^{dim}(**q**) for *n*_{hole} = 640/1152 = 0.556 and 544/1152 = 0.472, where the *d*_{xy}-wave and the extended-*s* + \(d_{x^2 - y^2}\)-wave SC are stabilized, respectively. For *n*_{hole} = 0.556, *S*^{dim}(**q**) peaks around (0, ±*π*/*R*_{y}) that are favorable for the *d*_{xy}-wave SC, while for *n*_{hole} = 0.472, *S*^{dim}(**q**) shows frustrated spin structure that are favorable for the extended-*s* + \(d_{x^2 - y^2}\)-wave SC. These are consistent with the mechanism of SC described above. On the other hand, *N*^{dim}(**q**) peaks around (0, ±3*π*/2*R*_{y}) for both *n*_{hole} = 0.556 and 0.472. This is because they are located near the instability of 3-fold CO-1 and 2, respectively, and the corresponding wave vectors are the same.

The superconducting condensation energy Δ*E* and spin/charge correlations are contrasting between the two SC phases. Figure 5e, f show the \(U/t_{b_1}\) dependence of *N*_{3-fold}, *S*_{peak}, and Δ*E*. *N*_{3-fold} = *N*^{dim}(0, ±3*π*/2*R*_{y}) and *S*_{peak} = *S*^{dim}(0, ±*π*/*R*_{y}) for *n*_{hole} = 0.556. The peak position in *S*^{dim}(**q**) changes along the upper and lower edges of the Brillouin zone for *n*_{hole} = 0.472 and we take the maximum value for *S*_{peak}. For *n*_{hole} = 0.556, Δ*E* does enhance, but despite the more rapid increase of *N*_{3-fold} toward 3-fold CO-1 instability, it rather follows *S*_{peak}. This is consistent with the usual AF spin fluctuation picture discussed in high-*T*_{c} cuprates. On the other hand, for *n*_{hole} = 0.472, Δ*E* is greatly enhanced following the rapid increase of *N*_{3-fold}, which suggests the close correlation between the extended-*s* + \(d_{x^2 - y^2}\)-wave SC and the 3-fold CO-2 type charge fluctuation.

## Discussion

Let us note that the intersite Coulomb interactions *V*_{ij} are indispensable to the stability of SC, not only for the extended-*s* + \(d_{x^2 - y^2}\)-wave SC phase but also for the metastable *d*_{xy}-wave SC phase. As we have shown in the previous work^{29}, no long-range ordered phases are stabilized in the absence of *V*_{ij} for the undoped condition (or unphysically large *U* is necessary). This is also the case in the doped condition studied here. This indicates that in *κ*-(ET)_{2}X, the charge degree of freedom is still active even with large *U* and the cooperation between *U* and *V*_{ij} induces various phases such as the DAF, PCO, 3-fold COs, and SC. Especially, the extended-*s* + \(d_{x^2 - y^2}\)-wave SC is enhanced toward both polar and 3-fold type CO instabilities. Recent experiment on photoinduced phase transition in *κ*-Br suggests that the polar charge oscillation is enhanced near the superconducting transition^{60}, consistent with our picture that the SC is enhanced toward the CO instability.

Experimentally, the symmetry of the SC for the undoped case is still controversial^{38,61,62,63,64,65,66,67}. This is consistent with our result suggesting that the competition between the two types of SC is most pronounced near the undoped *n*_{hole} = 1/2. Slight modification of parameters may alter the stability of the two. Now in *κ*-Cl, both electron and hole doping are realized using the EDLT^{46,47}. The filling-control MI transition and the emergence of SC have been confirmed with significant electron-hole doping asymmetry. Our result suggests that the different SC appears by carrier doping; the *d*_{xy}-wave SC on the hole-doped side and the extended-*s* + \(d_{x^2 - y^2}\)-wave SC on the electron-doped side. The former is due to the strong AF spin fluctuation, similar to high-*T*_{c} cuprates, and the latter is the consequence of frustrated spin structure and enhanced charge correlation under the geometrical frustration characteristic of the *κ*-type ET compounds. Experiments for the chemically doped *κ*-(ET)_{4}Hg_{3−δ}Y_{8} suggest similarities with the high-*T*_{c} cuprates, especially for the non-Fermi-liquid behaviors above the SC transition temperature^{42,43,44}. This is overall consistent with our results because this system is hole-doped. However, the tight-binding parameters for this system are rather closer to the isotropic triangular-lattice like arrangement of dimers, where the AF phase is expected to be unfavored because of the stronger geometrical frustration. Indeed, the temperature dependence of the spin susceptibility is similar to that of *κ*-CN^{45}, which does not show any magnetic long-range orders. It is difficult to fully explain the character of *κ*-(ET)_{4}Hg_{3−δ}Y_{8} within the present study and the analysis with appropriate tight-binding parameters is necessary for further understanding.

The nonmagnetic (candidate of gapless spin-liquid) insulator and neighboring SC observed in *κ*-CN are also an intriguing phenomena for the undoped case. Although the tight-binding parameters and the resulting electronic structure of *κ*-CN is different from the present study, the intradimer charge degree of freedom and the intersite Coulomb interactions should play crucial roles as pointed out previously^{29}. Indeed, the anomalous dielectric response^{36} and Raman spectroscopy^{37} indicate that the intradimer charge degree of freedom is active in *κ*-CN. The stability of spin-liquid and SC phase in the 1/2-filled Hubbard model on the anisotropic triangular lattice, which is the approximate dimer model of *κ*-CN, is still controversial despite the long and extensive studies^{12,13,14,15,17,19,20,21,22,23,24,25}. The analysis for the 3/4-filled EHM employed in this study will be an alternative way to investigate this issue and it is left for future studies.

## Methods

### Details of the VMC method

Here, we show the details of the VMC method. Trial wave function is a Gutzwiller-Jastrow type, \(\left| {\mathrm{\Psi }} \right\rangle = P_{{\mathrm{J}}_{\mathrm{c}}}P_{{\mathrm{J}}_{\mathrm{s}}}\left| {\mathrm{\Phi }} \right\rangle\). |Φ> is a one-body part constructed by diagonalizing the one-body Hamiltonian including the off-diagonal elements {*D*}, {*M*}, and {Δ} to treat long-range orders of charge, spin, and SC, respectively. The renormalized hopping integrals (\(\tilde t_{b_1},\tilde t_{b_2},\tilde t_p,\tilde t_q\)) are also included in |Φ> as variational parameters, where \(\tilde t_{b_1} = t_{b_1}\) is fixed as a unit. \(P_{{\mathrm{J}}_{\mathrm{c}}} = {\mathrm{exp}}\,[ - \mathop {\sum}\nolimits_{i,j} \,v_{ij}^{\mathrm{c}}n_in_j]\) and \(P_{{\mathrm{J}}_{\mathrm{s}}} = {\mathrm{exp}}\,[ - \mathop {\sum}\nolimits_{i,j} \,v_{ij}^{\mathrm{s}}s_i^zs_j^z]\) are charge and spin Jastrow factors which control long-range charge and spin correlations, respectively. Here, \(s_i^z = n_{i \uparrow } - n_{i \downarrow }\), and we assume \(v_{ij}^{\mathrm{c}} = v^{\mathrm{c}}(|r_i - r_j|)\) and \(v_{ij}^{\mathrm{s}} = v^{\mathrm{s}}(|r_i - r_j|)\), where *r*_{i} is the position of molecular site *i*. The variational parameters in |Ψ> are \(\tilde t_{b_2}\), \(\tilde t_p\), \(\tilde t_q\), {*D*}, {*M*}, \(\{ \Delta \}\), \(\{ v_{ij}^{\mathrm{c}}\}\), and \(\{ v_{ij}^{\mathrm{s}}\}\), and they are simultaneously optimized using the stochastic reconfiguration method^{68}. The total number of molecular sites is 4 × *L* × *L*/2 = 2*L*^{2} and varied from *L* = 12 to *L* = 24 with antiperiodic boundary conditions in both *x* and *y* directions of the primitive lattice vectors (see Fig. 1a).

The one-body part |Φ> for the PM, DAF, and PCO states can be obtained by diagonalizing the one-body Hamiltonian,

with the hopping matrix elements given as

where \(c_{{\bf{k}}m\sigma }^\dagger\) (*c*_{kmσ}) is a creation (annihilation) operator of electron at molecule *m*(=1–4), as indicated in Fig. 1a, with momentum *k* and spin *σ*(=↑,↓), and *s*_{σ} = 1(−1) for *σ* = ↑(↓). \(M_m^z\) is a variational parameter which induces the staggered AF long-range order aligned to the *z* direction for molecule *m*; \(M_1^z = M_2^z = M_3^z = M_4^z = 0\) for the PM state, \(M_1^z = M_2^z = M_3^z = M_4^z \ne 0\) for the DAF state, and \(M_1^z = M_3^z \ne M_2^z = M_4^z\) for the PCO state. \(a_{{\bf{k}}\alpha \sigma }^\dagger\) (*a*_{kασ}) in Eq. (7) is a creation (annihilation) operator of quasiparticle in band *α* with momentum *k* and spin *σ*(=↑,↓), obtained by diagonalizing Eq. (6), and \(\tilde E_\alpha (\mathbf{k})\) is the corresponding quasipaticle energy.

For the 3-fold CO-1 and 2 states, the unit cell is three times larger than the original one and contains 12 molecules, as shown in Fig. 2d, e. Therefore, |Φ> is obtained by diagonalizing a 12 × 12 matrix. The corresponding Weiss field is \(D_{m\sigma }\mathop {\sum}\nolimits_{{\bf{k}}\sigma } \,c_{{\bf{k}}m\sigma }^\dagger c_{{\bf{k}}m\sigma }\), which induces charge disproportionation and spin ordering within the unit cell through the variational parameters *D*_{mσ} for *m* = 1, 2, …, 12.

Finally, |Φ> for SC is obtained by diagonalizing the BCS-type mean-field Hamiltonian,

where Δ^{α} denotes gap function for band *α* and \(\xi _{\alpha} = \tilde E_{\alpha} ({\mathbf{k}}) - {\tilde {\mu}}\) is a quasiparticle energy measured from the renormalized chemical potential \(\tilde \mu\). \({\mathrm{\Delta }}^\alpha = {\mathrm{\Delta }}^\alpha (\{ {\mathrm{\Delta }}_i^\alpha \} )\) is constructed from the real space pairing \({\mathrm{\Delta }}_i^\alpha\) up to the 22nd neighbor dimers (*i* = 1, 2, …, 22). Since the band 3 and 4 are located much below the Fermi energy (see Fig. 1b), their contribution to the pairing is negligible and therefore we set Δ^{3} = Δ^{4} = 0.

## Data availability

The data that support the findings of this study are available from the corresponding author on reasonable request.

## Code availability

The code that support the findings of this study are available from the corresponding author on reasonable request.

## References

- 1.
Bednorz, J. G. & Muller, K. A. Possible high-

*T*_{c}superconductivity in the Ba-La-Cu-O system.*Z. Phys. B***64**, 189 (1986). - 2.
Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator transitions.

*Rev. Mod. Phys.***70**, 1039 (1998). - 3.
Uchida, S.

*High Temperature Superconductivity, The Role to Higher Critical Temperature*(Springer, Berlin, 2015). - 4.
Uemura, Y. J. et al. Basic similarities among cuprate, bismuthate, organic, Chevrel-phase, and heavy-fermion superconductors shown by penetration-depth measurements.

*Phys. Rev. Lett.***66**, 2665 (1991). - 5.
Sigrist, M. & Ueda, K. Phenomenological theory of unconventional superconductivity.

*Rev. Mod. Phys.***63**, 239 (1991). - 6.
Matsuda, Y., Izawa, K. & Vekhter, I. Nodal structure of unconventional superconductors probed by angle resolved thermal transport measurements.

*J. Phys. Condens. Matter***18**, R705 (2006). - 7.
Taillefer, L. Scattering and pairing in cuprate superconductors.

*Annu. Rev. Condens. Matter Phys.***1**, 51 (2010). - 8.
Ardavan, A. et al. Recent topics of organic superconductors.

*J. Phys. Soc. Jpn.***81**, 011004 (2012). - 9.
Lebed, A.

*The Physics of Organic Superconductors and Conductors*(Springer, New York, 2008). - 10.
McKenzie, R. H. Similarities between organic and cuprate superconductors.

*Science***278**, 820 (1997). - 11.
Kino, H. & Fukuyama, H. Phase diagram of two-dimensional organic conductors: (BEDT-TTF)

_{2}X.*J. Phys. Soc. Jpn.***65**, 2158 (1996). - 12.
Kino, H. & Kontani, H. Phase diagram of superconductivity on the anisotropic triangular lattice Hubbard model: An effective model of

*κ*-(BEDT-TTF) salts.*J. Phys. Soc. Jpn.***67**, 3691 (1998). - 13.
Kondo, H. & Moriya, T. Spin fluctuation-induced superconductivity in organic compounds.

*J. Phys. Soc. Jpn.***67**, 3695 (1998). - 14.
Schmalian, J. Pairing due to spin fluctuations in layered organic superconductors.

*Phys. Rev. Lett.***81**, 4232 (1998). - 15.
Morita, H., Watanabe, S. & Imada, M. Nonmagnetic insulating states near the Mott transitions on lattices with geometrical frustration and implications for

*κ*-(ET)_{2}Cu_{2}(CN)_{3}.*J. Phys. Soc. Jpn.***71**, 2109 (2002). - 16.
Liu, J., Schmalian, J. & Trivedi, N. Pairing and superconductivity driven by strong quasiparticle renormalization in two-dimensional organic charge transfer salts.

*Phys. Rev. Lett.***94**, 127003 (2005). - 17.
Kyung, B. & Tremblay, A.-M. S. Mott transition, antiferromagnetism, and d-Wave superconductivity in two-dimensional organic conductors.

*Phys. Rev. Lett.***97**, 046402 (2006). - 18.
Sahebsara, P. & Sénéchal, D. Antiferromagnetism and superconductivity in layered organic conductors: variational cluster approach.

*Phys. Rev. Lett.***97**, 257004 (2006). - 19.
Koretsune, T., Motome, Y. & Furusaki, A. Exact diagonalization study of Mott transition in the Hubbard model on an anisotropic triangular lattice.

*J. Phys. Soc. Jpn.***76**, 074719 (2007). - 20.
Watanabe, T., Yokoyama, H., Tanaka, Y. & Inoue, J. Predominant magnetic states in the Hubbard model on anisotropic triangular lattices.

*Phys. Rev. B***77**, 214505 (2008). - 21.
Shinaoka, H., Misawa, T., Nakamura, K. & Imada, M. Mott transition and phase diagram of

*κ*-(BEDT-TTF)_{2}Cu(NCS)_{2}studied by two-dimensional model derived from ab initio method.*J. Phys. Soc. Jpn.***81**, 034701 (2012). - 22.
Dayal, S., Clay, R. T. & Mazumdar, S. Absence of longrange superconducting correlations in the frustrated half-filled-band Hubbard model.

*Phys. Rev. B***85**, 165141 (2012). - 23.
Tocchio, L. F., Feldner, H., Becca, F., Valentí, R. & Gros, C. Spin-liquid versus spiral-order phases in the anisotropic triangular lattice.

*Phys. Rev. B***87**, 035143 (2013). - 24.
Laubach, M., Thomale, R., Platt, C., Hanke, W. & Li, G. Phase diagram of the Hubbard model on the anisotropic triangular lattice.

*Phys. Rev. B***91**, 245125 (2015). - 25.
Shirakawa, T., Tohyama, T., Kokalj, J., Sota, S. & Yunoki, S. Ground-state phase diagram of the triangular lattice Hubbard model by the density-matrix renormalization group method.

*Phys. Rev. B***96**, 205130 (2017). - 26.
Kuroki, K., Kimura, T., Arita, R., Tanaka, Y. & Matsuda, Y. d

_{x}^{2}_{−y}^{2}- versus d_{xy}- like pairings in organic superconductors*κ*-(BEDT-TTF)_{2}X.*Phys. Rev. B***65**, 100516(R) (2002). - 27.
Sekine, A., Nasu, J. & Ishihara, S. Polar charge fluctuation and superconductivity in organic conductors.

*Phys. Rev. B***87**, 085133 (2013). - 28.
Guterding, D., Altmeyer, M., Jeschke, H. O. & Valentí, R. Near-degeneracy of extended s + d

_{x}^{2}_{−}_{y}^{2}and d_{xy}order parameters in quasi-two-dimensional organic superconductors.*Phys. Rev. B***94**, 024515 (2016). - 29.
Watanabe, H., Seo, H. & Yunoki, S. Phase competition and superconductivity in

*κ*-(BEDT-TTF)_{2}X: Importance of intermolecular Coulomb interactions.*J. Phys. Soc. Jpn.***86**, 033703 (2017). - 30.
Powell, B. J., Kenny, E. P. & Merino, J. Dynamical reduction of the dimensionality of exchange interactions and the “spin-Liquid” phase of

*κ*-(BEDT-TTF)_{2}X.*Phys. Rev. Lett.***119**, 087204 (2017). - 31.
Kaneko, R., Tocchio, L. F., Valentí, R. & Becca, F. Charge orders in organic charge-transfer salts.

*New J. Phys.***19**, 103033 (2017). - 32.
Zantout, K., Altmeyer, M., Backes, S. & Valentí, R. Superconductivity in correlated BEDT-TTF molecular conductors: critical temperatures and gap symmetries.

*Phys. Rev. B***97**, 014530 (2018). - 33.
Hotta, C. Quantum electric dipoles in spin-liquid dimer Mott insulator

*κ*-ET_{2}Cu_{2}(CN)_{3}.*Phys. Rev. B***82**, 241104(R) (2010). - 34.
Naka, M. & Ishihara, S. Electronic ferroelectricity in a dimer Mott insulator.

*J. Phys. Soc. Jpn.***79**, 063707 (2010). - 35.
Manna, R. S., de Souza, M., Brühl, A., Schlueter, J. A. & Lang, M. Lattice effects and entropy release at the low-temperature phase transition in the spin-liquid candidate

*κ*-(BEDT-TTF)_{2}Cu_{2}(CN)_{3}.*Phys. Rev. Lett.***104**, 016403 (2010). - 36.
Abdel-Jawad, M. et al. Anomalous dielectric response in the dimer Mott insulator

*κ*-(BEDT-TTF)_{2}Cu_{2}(CN)_{3}.*Phys. Rev. B***82**, 125119 (2010). - 37.
Yakushi, K., Yamamoto, K., Yamamoto, T., Saito, Y. & Kawamoto, A. Raman spectroscopy study of charge fluctuation in the spin-liquid candidate

*κ*-(BEDT-TTF)_{2}Cu_{2}(CN)_{3}.*J. Phys. Soc. Jpn.***84**, 084711 (2015). - 38.
Guterding, D. et al. Evidence for eight-node mixed-symmetry superconductivity in a correlated organic metal.

*Phys. Rev. Lett.***116**, 237001 (2016). - 39.
Kagawa, F., Miyagawa, K. & Kanoda, K. Unconventional critical behaviour in a quasi-two-dimensional organic conductor.

*Nature***436**, 534 (2005). - 40.
Kurosaki, Y., Shimizu, Y., Miyagawa, K., Kanoda, K. & Saito, G. Mott transition from a spin liquid to a fermi liquid in the spin-frustrated organic conductor

*κ*-(ET)_{2}Cu_{2}(CN)_{3}.*Phys. Rev. Lett.***95**, 177001 (2005). - 41.
Kanoda, K. Metal-Insulator transition in

*κ*-(ET)_{2}X and (DCNQI)_{2}M: two contrasting manifestation of electron correlation.*J. Phys. Soc. Jpn.***75**, 051007 (2006). - 42.
Naito, A. et al. Anomalous enhancement of electronic heat capacity in the organic conductors

*κ*-(BEDT-TTF)_{4}Hg_{3−δ}X_{8}(X = Br, Cl).*Phys. Rev. B***71**, 054514 (2005). - 43.
Taniguchi, H. et al. Anomalous pressure dependence of superconductivity in layered organic conductor,

*κ*-(BEDT-TTF)_{4}Hg_{2.89}Br_{8}.*J. Phys. Soc. Jpn.***76**, 113709 (2007). - 44.
Oike, H., Miyagawa, K., Taniguchi, H. & Kanoda, K. Pressure-induced Mott transition in an organic superconductor with a finite doping level.

*Phys. Rev. Lett.***114**, 067002 (2015). - 45.
Oike, H. et al. Anomalous metallic behaviour in the doped spin liquid candidate

*κ*-(BEDT-TTF)_{4}Hg_{2.89}Br_{8}.*Nat. Commun.***8**, 756 (2017). - 46.
Kawasugi, Y. et al. Electron-hole doping asymmetry of Fermi surface reconstructed in a simple Mott insulator.

*Nat. Commun.***7**, 12356 (2016). - 47.
Sato, Y., Kawasugi, Y., Suda, M., Yamamoto, H. M. & Kato, R. Critical behavior in doping-driven metal-insulator transition on single-crystalline organic Mott-FET.

*Nano Lett.***17**, 708 (2017). - 48.
Seo, H., Hotta, C. & Fukuyama, H. Toward systematic understanding of diversity of electronic properties in low-dimensional molecular solids.

*Chem. Rev.***104**, 5005 (2004). - 49.
Seo, H. Charge ordering in organic ET compounds.

*J. Phys. Soc. Jpn.***69**, 805 (2000). - 50.
Miyagawa, K., Kawamoto, A. & Kanoda, K. Proximity of pseudogapped superconductor and commensurate antiferromagnet in a quasi-two-dimensional organic system.

*Phys. Rev. Lett.***89**, 017003 (2002). - 51.
Koretsune, T. & Hotta, C. Evaluating model parameters of the

*κ*- and*β*′-type Mott insulating organic solids.*Phys. Rev. B***89**, 045102 (2014). - 52.
McMillan, W. L. Ground state of liquid He

^{4}.*Phys. Rev.***138**, A442 (1965). - 53.
Ceperley, D., Chester, G. V. & Kalos, M. H. Monte Carlo simulation of a many-fermion study.

*Phys. Rev. B***16**, 3081 (1977). - 54.
Yokoyama, H. & Shiba, H. Variational Monte-Carlo studies of Hubbard Model. I.

*J. Phys. Soc. Jpn.***56**, 1490 (1987). - 55.
Drichko, N. et al. Metallic state and charge-order metal-insulator transition in the quasi-two-dimensional conductor

*κ*-(BEDT-TTF)_{2}Hg(SCN)_{2}Cl.*Phys. Rev. B***89**, 075133 (2014). - 56.
Yoshimi, K., Seo, H., Ishibashi, S. & Brown, S. E. Tuning the magnetic dimensionality by charge ordering in the molecular TMTTF salts.

*Phys. Rev. Lett.***108**, 096402 (2012). - 57.
Mazin, I. I., Singh, D. J., Johannes, M. D. & Du, M. H. Unconventional superconductivity with a sign reversal in the order parameter of LaFeAsO

_{1−x}F_{x}.*Phys. Rev. Lett.***101**, 057003 (2008). - 58.
Kuroki, K. et al. Unconventional pairing originating from the disconnected fermi surfaces of superconducting LaFeAsO

_{1−x}F_{x}.*Phys. Rev. Lett.***101**, 087004 (2008). - 59.
Powell, B. J. & McKenzie, R. H. Symmetry of the superconducting order parameter in frustrated systems determined by the spatial anisotropy of spin correlations.

*Phys. Rev. Lett.***98**, 027005 (2007). - 60.
Kawakami, Y. et al. Nonlinear charge oscillation driven by a single-cycle light field in an organic superconductor.

*Nat. Photonics***12**, 474 (2018). - 61.
Elsinger, H. et al.

*κ*-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Br: a fully gapped strong-coupling superconductor.*Phys. Rev. Lett.***84**, 6098 (2000). - 62.
Izawa, K., Yamaguchi, H., Sasaki, T. & Matsuda, Y. Superconducting gap structure of

*κ*-(BEDT-TTF)_{2}Cu(NCS)_{2}probed by thermal conductivity tensor.*Phys. Rev. Lett.***88**, 027002 (2001). - 63.
Taylor, O. J., Carrington, A. & Schlueter, J. A. Specific-heat measurements of the gap structure of the organic superconductors

*κ*-(ET)_{2}Cu[N(CN)_{2}]Br and*κ*-(ET)_{2}Cu(NCS)_{2}.*Phys. Rev. Lett.***99**, 057001 (2007). - 64.
Ichimura, K., Takami, M. & Nomura, K. Direct observation of d-wave superconducting gap in

*κ*-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Br with scanning tunneling microscopy.*J. Phys. Soc. Jpn.***77**, 114707 (2008). - 65.
Malone, L., Taylor, O. J., Schlueter, J. A. & Carrington, A. Location of gap nodes in the organic superconductors

*κ*-(ET)_{2}Cu(NCS)_{2}and*κ*-(ET)_{2}Cu[N(CN)_{2}]Br determined by magnetocalorimetry.*Phys. Rev. B***82**, 014522 (2010). - 66.
Milbradt, S. et al. In-plane supefluid density and microwave conductivity of the organic superconductor

*κ*-(BEDT-TTF)_{2}Cu[N(CN)_{2}]Br: Evidence for d-wave pairing and resilient quasiparticles.*Phys. Rev. B***88**, 064501 (2013). - 67.
Oka, Y. et al. Tunneling spectroscopy in organic superconductor

*κ*-(BEDT-TTF-d[3,3])_{2}Cu[N(CN)_{2}]Br.*J. Phys. Soc. Jpn.***84**, 064713 (2015). - 68.
Sorella, S. Generalized Lanczos algorithm for variational quantum Monte Carlo.

*Phys. Rev. B***64**, 024512 (2001).

## Acknowledgements

The authors thank R. Kato, Y. Kawasugi, H. Itoh, S. Iwai, Y. Kawakami, and M. Naka for useful discussions. The computation has been done using the facilities of the Supercomputer Center, Institute for Solid State Physics, University of Tokyo. This work has been supported by JSPS KAKENHI (Grant Nos 26800198, 26400377, 16H02393, and 18H01183) and Waseda University Grant for Special Research Projects (Project Number: 2018B-352).

## Author information

### Affiliations

### Contributions

H.W. performed the VMC simulations and prepared the figures. Results were analyzed and the paper was written by H.W., H.S., and S.Y.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Peer review information:** *Nature Communications* thanks Kazushi Kanoda and other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Publisher’s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Watanabe, H., Seo, H. & Yunoki, S. Mechanism of superconductivity and electron-hole doping asymmetry in *κ*-type molecular conductors.
*Nat Commun* **10, **3167 (2019). https://doi.org/10.1038/s41467-019-11022-1

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.