Next Article in Journal
Analysis of Waveguides on Lithium Niobate Thin Films
Next Article in Special Issue
Raman Scattering as a Probe of the Magnetic State of BEDT-TTF Based Mott Insulators
Previous Article in Journal
Comparative Study on the Sand Bioconsolidation through Calcium Carbonate Precipitation by Sporosarcina pasteurii and Bacillus subtilis
Previous Article in Special Issue
Spatially Nonuniform Superconductivity in Quasi-Two-Dimensional Organic Charge-Transfer Salts
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Crystals
Volume 8
Issue 5
10.3390/cryst8050190
crystals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrodynamics in Organic Dimer Insulators Close to Mott Critical Point

by
Marko Pinterić
1,2,
David Rivas Góngora
1,
Željko Rapljenović
1,
Tomislav Ivek
1,
Matija Čulo
1,
Bojana Korin-Hamzić
1,
Ognjen Milat
1,
Branko Gumhalter
1,
Predrag Lazić
3,
Miriam Sanz Alonso
4,
Weiwu Li
4,
Andrej Pustogow
4,
Guilherme Gorgen Lesseux
4,
Martin Dressel
4,* and
Silvia Tomić
1,*
1
Institut za fiziku, Bijenička 46, HR-10000 Zagreb, Croatia
2
Faculty of Civil Engineering, University of Maribor, Smetanova 17, SI-2000 Maribor, Slovenia
3
Rudjer Bošković Institute, Bijenička 56, HR-10000 Zagreb, Croatia
4
1. Physikalisches Institut, Universität Stuttgart, Pfaffenwaldring 57, D-70550 Stuttgart, Germany
*
Authors to whom correspondence should be addressed.
Crystals 2018, 8(5), 190; https://doi.org/10.3390/cryst8050190
Submission received: 26 March 2018 / Revised: 18 April 2018 / Accepted: 18 April 2018 / Published: 27 April 2018
(This article belongs to the Special Issue Advances in Organic Conductors and Superconductors)
Browse Figures
Versions Notes

Abstract

:
Organic layered charge-transfer salts κ -(BEDT-TTF) 2 X form highly frustrated lattices of molecular dimers in which strong correlations give rise to Mott insulating states situated close to the metal-to-insulator phase boundary. The salts κ -(BEDT-TTF) 2 Cu 2 (CN) 3 and κ -(BEDT-TTF) 2 Ag 2 (CN) 3 have been considered as prime candidates for a quantum spin liquid, while κ -(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl has been suggested as a prototypical charge-order-driven antiferromagnet. In this paper, we summarize and discuss several key results, including some not reported previously, obtained in search to clarify the competition of these two ground states. The origin of anomalous dielectric response found at low temperatures in all three salts is also discussed. We conclude by pointing out the relevant new insights into the role of frustration and random disorder in the suppression of magnetic ordering and formation of the spin liquid state.

1. Introduction

Layered organic solids κ -(BEDT-TTF) 2 X (X = Cu[N(CN) 2 ]Cl, Cu 2 (CN) 3 and Ag 2 (CN) 3 ) are systems in which strong correlations lead to Mott insulating states situated close to the metal-to-insulator phase boundary [1,2]. These systems are charge-transfer salts, i.e., they consist of conducting BEDT-TTF layers separated by non-conducting anion layers (Figure 1). Within BEDT-TTF layers, the triangular lattices of molecular dimers with one hole per dimer are formed with different degrees of frustration (Figure 2) [3,4,5,6]. Moderate pressure in the range of a few hundred bars up to 10 kbar is sufficient to recover metallicity and superconducting ground state. At ambient pressure in these salts two different ground states can be established. Canted antiferromagnetic, i.e., weak ferromagnetic, state is formed in κ -(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl (denoted as κ -Cl) [7,8,9]. On the other hand, quantum spin liquid (QSL) develops in the two latter salts κ -(BEDT-TTF) 2 X , X = Cu 2 (CN) 3 and Ag 2 (CN) 3 (denoted as κ -CuCN and κ -AgCN, respectively) [10,11,12]. QSL is commonly considered as a highly correlated fluctuating quantum spin state expected to appear in geometrically frustrated systems when quantum fluctuations are strong enough to preclude spin ordering down to zero temperature [13,14]. Historically, QSL was predicted theoretically as a resonating valence bond state to be the ground state of a perfect triangular lattice S = 1 / 2 spin antiferromagnet instead of the conventional long-range antiferromagnetic (AF) order [15]. Experimentally, however, few real materials with triangular lattice, no spontaneously broken symmetry, and emergent fractional excitations have been discovered. Among them stand out the most recent discoveries of two novel QSL types. First is a high-temperature QSL in the layered transition metal dichalcogenide 1T-TaS 2 [16] as a long-awaited prediction by Anderson. Second is a quantum liquid of magnetic and electric dipoles developed in a hydrogen-bonded molecular compound κ -H 3 (Cat-EDT-TTF) 2 [17].
Common to all three κ -(BEDT-TTF) 2 X organic charge-transfer salts is an anomalous dielectric response [8,18,19,20,21,22,23], the origin of which remained elusive due to the missing experimental evidence of any sizeable electric dipoles on molecular dimers [19,23,24,25] proposed theoretically [26,27,28,29,30,31,32,33,34]. However, relying on the theory and based on the observation of Curie-like dielectric peak, hysteresis and time-dependent phenomena, a dipolar order has been proposed to set in κ -Cl acting as a driving force for the magnetic ground state [18]. Related open issue in this salt is that, instead of Curie-like peak, in the majority of measured single crystals, a relaxor response has been detected, suggesting that the dielectric response rather arises from the motion of charge domain walls between insulating ferromagnetic domains in the background impurity potential [8,19]. In two salts κ -CuCN and κ -AgCN with quantum spin liquid ground state, dielectric response is consistently found to be of relaxor type in all measured single crystals. The suppression of magnetic state has been mostly associated with highly frustrated molecular dimer lattices, but it has been argued that an additional mechanism, such as an exotic spin-dipolar coupling, is needed to induce quantum fluctuations which are strong enough to suppress magnetic order [26,27,28,29,30,31,32,33,34]. The reason for that is twofold. First is the proximity of Mott insulator-to-metal phase boundary where the charge degrees of freedom start to delocalize and thus a spin-only description cannot be sufficient. Second is the inability of frustrated triangular lattices to destroy the long-range magnetic order on their own [35,36,37]. The role of random disorder in this competition may also be important, as indicated by an experimental study on κ -Cl showing that X-rays suppress magnetic order and allow a spin liquid to form [38]. Finally, we have recently suggested that the non-conducting anions which contain frustration and disorder in the organic spin liquids may be a critical ingredient needed for the prevalence of quantum spin liquid over magnetic ordering [39].
In this paper, we address the issue of anomalous dielectric response and inherent randomness, and explore its origin and role in the competition between the AF and QSL ground state formation.

2. Antiferromagnet with Ferroelectric Character κ -(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl

At ambient pressure, κ -Cl system behaves as an insulator in dc limit in the whole temperature range, but, in the infrared range, a large temperature-dependent charge gap (Mott gap) develops only at temperatures lower than 50 K [19,40]. The system was reported to show zero resistance and weak superconducting diamagnetic response at pressures as small as 100 bar, and a bulk superconductivity (SC) below 12.5 K at 300 bar [41,42,43]. On cooling, dc resistivity increases only weakly down to about 100 K. This behavior is observed for majority of measured samples; the behavior of three single crystals is displayed in Figure 3, Inset in Left panel.
On further cooling below 100 K, the dc resistivity increases steeply, indicating an insulator-to- insulator phase transition, right in the temperature range where the Mott gap develops. Interestingly, the rate by which dc resistivity rises seems to be much larger perpendicular to the dimer planes (i.e., along the b-axis) than within them [18,44]. Overall, dc resistivity data accumulated by us and the other authors however indicates that the phase transition temperature and the width of the transition extracted from the logarithmic resistivity derivative differ for different single crystals. Figure 3 shows distinct behaviors of the dc resistivity along the b-axis for three single crystals of κ -Cl. The resistivity derivative of S2, as well as of the sample in Ref. [19] (see Figure 3 in [19]), shows a broad maximum centered at about 30 K, indicative of the formation of a low-temperature phase at short-range scales. Conversely, samples S1 and S3 present a peak that is twice as narrow, higher, and positioned at a slightly lower temperature of 27 K. However, the full width at half maximum is ten times larger than what is typical for the phase transition developed at the long-range scales, indicating a certain amount of disorder even in the crystals S1 and S3.
The nature of phase transition, based on dielectric data displaying a peak in the same temperature range (Figure 4, Right), was proposed to be ferroelectric and ascribed to the ordering of electric dipoles on BEDT-TTF dimers [18]. Moreover, it was argued that the dipolar ordering releases the frustration of dimers’ spins on the triangular lattice and in this way drives the magnetic ordering which sets in at same temperatures. The exciting proposal of multiferroicity awaits for its final confirmation due to missing evidence for inversion symmetry breaking [19,25]. Thus, the nature of phase transitions detected in dc resistivity and dielectric response, and the origin of disorder in κ -Cl are still a mystery. We come back to this issue at the end of this section. An additional broad feature, situated at about 50 K, is detected only for sample S3 (Figure 3, Right), right in the temperature range where the Mott gap opens.
Further, the most intriguing behavior is found for sample S2. Upon reducing temperature the resistivity rises with a much smaller rate, reaches maximum at about 12 K and falls down rapidly. Note that 12 K is the temperature below which bulk SC sets in at 300 bar. Notably, a similar drop of resistivity measured in dc limit, as well as at 35 GHz, has been reported before independently by Sushko et al. [45] and one of us [44], and attributed to SC islands growing within the insulating sample bulk. The suggestion has been supported by magnetization measurements indicating that the low temperature state of κ -Cl below critical pressure of 300 bar is an inhomogeneous mixture of SC and insulating magnetic phases [45]. Indeed, for sample S2, we do find evidence for the presence of SC domains at T < T SC ≈ 12 K. These domains do not percolate down to 5 K since the resistivity remains finite (Figure 3). At these low temperatures only, a non-linear decrease of the low-frequency conductivity is found when increasing ac voltage (not shown) or applying magnetic field (Figure 5, Left). At zero magnetic field, there is a rise in conductivity below about 12 K, as observed in dc resistivity measurements. Increasing the magnetic field shifts the SC phase transition to lower temperatures. For magnetic field of 5 T the SC phase transition is suppressed down to 5 K.
Next, we address the dielectric function measured as a function of temperature and frequency along the E ∥ b i.e., perpendicular to the BEDT-TTF dimer layers. Figure 4 presents the real part of dielectric function ε ′ of three representative samples of κ -Cl. Common to all of them is that below 40 K, ε ′ features a pronounced peak, the high-temperature shoulders of which align to follow the same curve. Conversely, the peak’s amplitudes and the positions in temperature strongly differ between them. The former seem to correlate with the magnitude of dc conductivity, but might also be influenced by different levels of disorder present in the three samples. The ε ′ peaks of sample S1 and S3 are situated at about 25 K, showing a weak relaxor and a Curie-like behavior i.e., nearly frequency independent position of the peak, and reach similar values ≤100 for frequencies ≥200 Hz. A remarkably distinct result is observed for sample S2: the ε ′ peak is shifted to lower temperatures and centered at T c ≈ 9 K for all applied frequencies attaining a value as high as 1800 at 200 Hz. Notably, T c seems to coincide with T SC . Applied magnetic field reduces the onset of conductivity rise (Figure 5, Left) and concomitantly smooths the low-temperatures sides of the ε ′ peak and lowers its amplitude, the effect showing up only at temperatures below T SC ≈ 12 K (Figure 5, Right). Finally, the susceptibility behavior observed for sample S2 shows an increase, albeit rather broad, below about 22 K (Figure 6), as expected for the intra-layer antiferromagnetic ordering documented previously [7,8,9].
Importantly, for sample S2 both the magnetic phase transition and the insulator-to-insulator phase transition seem to be developed only at short-length scales, however they appear within the same temperature range as commonly observed in κ -Cl samples. On the other hand, the dielectric peak if taken as an indication of a ferroelectric-like ordering reveals itself as a feature separated from the phase transitions visible in dc transport and in susceptibility. Moreover, it is shifted to lower temperatures and appears concomitantly with the superconducting phase transition indicating the coexistence of ferroelectric-like order and SC. Thus, despite ferroelectricity being expected only in insulating systems, here it seems to coexist with dilute metallicity, indicating that ferroelectricity, together with magnetism, shares a common boundary with SC in the phase diagram. However, note that, while magnetic and SC are competing phases, the observed data suggest that ferroelectricity and SC support each other. Remarkably, Clay and Mazumdar recently developed a valence-bond theory of SC within which the SC state is reached by the destabilization of the molecular dimer order called paired-electron crystal (PEC) [31,32,46]. According to their theory, the low-temperature phase in sample S2 might be considered as consisting of segregated magnetic, PEC and SC phases. This outstanding result in organics bears similarities to the superconducting high- T c cuprates, where the proximity of a charge-ordered state to the AF and the SC state has been established.
It is challenging to reconcile this exciting proposal with the missing experimental evidence of structural inversion symmetry breaking needed for formation of ferroelectricity. Namely, X-ray diffraction measurements have shown that the space group of κ -Cl is P n m a implying that all eight BEDT-TTF molecules and all four dimers in the unit cell are crystallographically equivalent [47]. Interestingly, NMR data reported in [9] seem to indicate four nonequivalent spin dimers sites present already at room temperature, but remarkably the nonequivalence disappeared once the magnetic state was established. The experimentally evidenced absence of any considerable electric dipoles on BEDT-TTF dimer sites [25] has promoted an alternative interpretation for the relaxor-like dielectric response observed in the most number of κ -Cl samples. This scenario, backed up by the generally observed tendency to phase segregation, invokes charged domain-walls between insulating ferromagnetic domains in the disordered background of the canted AF i.e., weak ferromagnetic phase [8,19,22].
Here, we present vibrational spectroscopy data performed on sample S2 of κ -Cl. Figure 7 displays the most charge-sensitive intramolecular vibrational mode ν 27 centered at 1469 cm − 1 at 5 K , which involves the anti-symmetric C=C stretching of the rings within BEDT-TTF molecules, as a function of temperature [48,49,50]. The side band at 1473 cm − 1 is present at all temperatures due to inequivalent lattice sites [25] and becomes visible at low temperatures as the modes become more narrow. The small satellite at 1478 cm − 1 involves ethylene motions and is not related to the ν 27 mode [50]. A possible charge inequality among the donor molecules can be quantitatively estimated from the well-known linear frequency-charge relation concluded by Yamamoto et al. [49]. The charge disproportionation can be directly determined from the frequency separation of the respective ν 27 peaks via 2 δ ρ = δ ν 27 / ( 140 cm − 1 / e ) . In the case of κ -Cl, however, no splitting of the mode occurs on lowering the temperature down to 5 K, providing clear evidence that, even in sample S2 of κ -Cl, no sizeable static charge redistribution exceeding the limit of 2 δ ρ ≈ ± 0.01 e takes place. The observation is in full accord with previous reports on κ -Cl [19,25]. Given the estimated extremely small upper bound of charge imbalance and fast rate of charge oscillations within the dimer [25], such static or fluctuating electric dipoles associated with BEDT-TTF dimers can hardly explain the low-frequency dielectric response, in particular the one showing a sizeable Curie-like anomaly, as observed for sample S2 of κ -Cl. Importantly, in sample S2 the magnetic and dielectric anomaly, setting in at separated temperatures of 25 K and 9 K, seem to suggest that it is a ferroelectricity which is driven by the magnetism, and not vice versa as proposed in [18].
Finally, origin of disorder in single crystals of κ -Cl is a very intriguing issue. The presence of disorder is suggested by the broadness of phase transition anomalies seen in dc transport and in susceptibility, as well as relaxor-like dielectric anomalies [51]. However, here the possibility that disorder originates in the anion layer, as suggested in the case of κ -CuCN and κ -AgCN [21,23,24], must be discarded since all cyanide (CN) groups positioned between the copper (Cu) atoms are ordered in a zigzag line along the a-axis so that the anion layer consists of one-dimensional chains (Figure 8). Alternatively, disorder may arise in the organic dimer layers due to the formation of segregated phases of similar energy [46]. This possibility should be explored further.

3. Quantum Spin Liquids in Disordered κ -(BEDT-TTF) 2 Cu 2 (CN) 3 and κ -(BEDT-TTF) 2 Ag 2 (CN) 3

κ -CuCN and κ -AgCN are Mott–Anderson insulators in which combined effects of strong Coulomb interaction and randomness have been well documented, experimentally and theoretically. dc charge transport within dimer planes happens via two-dimensional variable-range hopping in a wide low-temperature range. At elevated temperatures, it crosses-over to nearest-neighbor hopping [21,23,52,53]. Anomalously broad anisotropic dielectric response exhibits the relaxor-like ferroelectric behavior: with lower frequencies the magnitude of peak increases and its position shifts to lower temperature indicating a gradual freezing out [20,21,23]. Figure 9 displays the temperature-dependent real part of the dielectric function at selected frequencies along the a * direction where dimer and anion layers alternate. Notably, the out-of-plane response is one to three orders of magnitude weaker in strength, which correlates with the anisotropy of dc conductivity. Extensive study on different single crystals from different syntheses have revealed strong sample dependence of the dc transport and the dielectric response parameters [21]. Dielectric peak becomes smaller for larger disorder and the difference becomes more pronounced upon reducing temperatures. Similar effects due to disorder have been found previously in the dielectric response of the charge-density wave and ferroelectric charge ordered states in quasi-one-dimensional inorganic and organic compounds [22,54,55]. Disorder also affects the width of dielectric spectrum the broadening of which decreases implying a lower degree of cooperativity, the effect typically occurring in glassy systems [56]. Cooperativity appears well established within the dimer planes, while it develops gradually upon decreasing temperature along the out-of-plane direction [21,23]. In addition, for κ -CuCN, sample-dependent dielectric and lattice anomalies were also detected in microwave [57] and thermal expansion measurements [58], respectively.
Dielectric response observed in κ -CuCN and κ -AgCN thus bears typical fingerprints of cooperative motion and glassy freezing in relaxor ferroelectrics [59]. Relaxor phenomenology assumes the existence of polar domains of low-symmetry and anticipates basically two mechanisms of dielectric response. It can be either due to the thermally activated reorientation of dipole moments within the domain, or to the motion of the interphase boundaries of the polar regions.
To examine if dipole moments emerge upon lowering the temperature in κ -CuCN and κ -AgCN, vibrational spectroscopy was employed to probe the most charge-sensitive intramolecular vibrational mode ν 27 , which involves the anti-symmetric ring C=C stretching of the BEDT-TTF molecule [23,25]. In Figure 10 and Figure 11, the ν 27 is shown for various temperatures in κ -CuCN and κ -AgCN, respectively. In contrast to κ -Cl, the side band at 1473 cm − 1 is absent in different lattice symmetry; the satellite slightly below 1480 cm − 1 is similar for all compounds discussed here. At 300 K the ν 27 shows as a very strong peak centered at 1460 cm − 1 , as expected for half a hole per BEDT-TTF molecule. The mode is rather broad indicating the presence of fast charge fluctuations. Upon cooling no splitting of the mode occurs down to 10 K, providing clear evidence that no sizeable static charge redistribution takes place in both κ -CuCN and κ -AgCN. However, there is a slight difference. In κ -CuCN, only one Fano mode was sufficient to describe the ν 27 allowing only to conclude that there is no static charge imbalance larger than 2 δ ρ ≈ ± 0.01 e (Figure 10). For κ -AgCN, two Fano functions were required to fit the spectra, revealing the two maxima which indicate that two unequal crystallographic sites exist in the unit cell of κ -AgCN (Figure 11) [23]. With a linear shift of 140 cm − 1 per unit charge, the charge imbalance 2 δ ρ is calculated by 2 δ ρ = δ ν 27 / ( 140 cm − 1 / e ) , where δ ν 27 is the difference in frequency positions between two maxima associated with two non-equal molecular or dimer sites. The two maxima are separated by 6 cm − 1 which would correspond to an average charge imbalance of 2 δ ρ ≈ ± 0.05 e .
Given that vibrational spectroscopy detects only a tiny charge imbalance and fast rate of oscillations, such dipoles cannot be invoked to explain the dielectric response at low frequencies. These fluctuating dipoles may show up in the microwave [57] and terahertz response [60]; however, the latter has been shown to involve coupled anion–dimer vibrations instead [24]. Alternatively, the collective motion of interphase boundaries of polar regions has to be considered as the pertinent mechanism in these systems. Therefore, the clue to understand relaxor properties of dielectric response lies in identifying the randomness and domain structure in κ -CuCN and κ -AgCN. The structure of these salts is commonly solved within P 2 1 / c symmetry implying that all four molecules and two dimers within the unit cell are equivalent [4,6,61]. This means that any charge disproportionation between molecular and/or dimer sites is excluded.
However, randomness is revealed by X-ray diffraction data already at 300 K which show that, in addition to ordered CN groups (chain CN) between Cu/Ag atoms in the chains along the b-axis, there are CN groups (bridging CN) located at the inversion centers. These bridging CN groups connect Cu/Ag atoms along the c-axis so that the anion network is formed in two-dimensions (Figure 12). The triangular coordination of Cu and Ag implies frustration since each Cu/Ag atom can be linked either to two N and one C atom, or to one N atom and two C atoms introducing intrinsic disorder. This disorder is extended via hydrogen bonds connecting CN and ethylene groups onto the molecular dimers where it alters the charge distribution, and thus can cause anti-phase boundaries separating frustration-limited domains [21]. The perturbation of BEDT-TTF layers is expected to be stronger in κ -CuCN, than in κ -AgCN, since in the former all contacts between bridging CN and the terminal ethylene groups of the BEDT-TTF molecule are shorter than the sum of van der Waals radii [23]. Note that the disorder level in κ -CuCN and κ -AgCN represents an intrinsic property of these charge-transfer salts.
To get an insight as how such a disorder affects the electronic ground state, ab initio density functional theory (DFT) calculations have been carried out on κ -CuCN and κ -AgCN [23,24]. Since two CN bridging groups in the unit cell occupy one of the two equally probable orientations, four different crystallographic configurations were constructed, based on experimental X-ray data, for κ -CuCN and κ -AgCN. In each configuration, the entire system, with atomic positions and unit cell parameters, was allowed to relax resulting in electronic states with reduced symmetry as compared to P 2 1 / c and quasi-degenerate in energy; all states lie only 5–25 meV above the relaxed state with minimum energy. Reduced symmetry implies that non-equivalent crystallographic sites (BEDT-TTF molecules and/or dimers) exist in the unit cell. In κ -CuCN, the space group of all four states was reduced to P 1 , while, in κ -AgCN, the energy states were associated with the configurations of higher symmetry P 2 1 and P c . These results indicate a higher level of disorder in κ -CuCN than in κ -AgCN, which explains why it was possible to resolve two Fano modes in the ν 27 mode in the κ -AgCN system.
Importantly, the quasi-degenerate electronic ground state implies a random domain structure in the real space of κ -CuCN and κ -AgCN at all temperatures. It thus supports the consideration of relaxor dielectric response explained above: low-energy excitations are caused by charged domain walls rather than by static or fluctuating electric dipoles. Upon cooling, screening gets gradually weaker due to the reduced number of charge carriers. This allows domains to be more “visible”, the number of domain walls increases, and concomitantly the motion of domain walls become increasingly correlated. The detailed description of how and why the dielectric response evolve in temperature can be found in Refs. [21,23]. In closing, we note that the dielectric response scenario presented in this section, was recently backed up by the theoretical approach developed to explain anomalous dielectric response in strongly correlated molecular solids [62]. Most importantly, in-depth X-ray diffraction measurements of κ -CuCN in the wide temperature range and refinements, complemented by DFT calculations, have demonstrated that two nonequivalent dimers exist in the unit cell due to the reduced symmetry of the whole crystal [63].

4. Van de Waals Interactions in κ -(BEDT-TTF) 2 X ; X = Cu[N(CN) 2 ]Cl, Cu 2 (CN) 3 and Ag 2 (CN) 3

While there exists a general consent in the community of organic crystals on the importance of vdW interaction in their formation, its importance has never been clearly evidenced and quantified in the organic materials κ -(BEDT-TTF) 2 X with X = Cu[N(CN) 2 ]Cl, Cu 2 (CN) 3 and Ag 2 (CN) 3 . These three kappa-BEDT-TTF materials, as well as other BEDT-TTF compounds are not simple molecular crystals, rather they should be classified as complex chemically bound charge transfer or cation–anion bulk solids. The previous DFT studies on different organic charge-transfer salts have neither included nor discussed the van der Waals (vdW) interactions at all [64,65,66,67,68,69,70]. The only exception is found in our work [23,24,39].
In the present DFT calculations, the experimental structure obtained from X-ray diffraction measurements at 100 K, 150 K and 100 K was used for the atomic coordinates of κ -CuCN, κ -AgCN and κ -Cl, respectively [23,24,47]. Two types of DFT functionals have been used: the standard semilocal PBE functional and nonlocal vdW-DF functional. The structure was relaxed to allow for the change of volume and shape of the unit cell together with the atomic positions. The parameters of the resulting relaxed state structures with minimum energy are displayed in Table 1 ( κ -CuCN), Table 2 ( κ -AgCN) and Table 3 ( κ -Cl).
The inadequacy of DFT calculations without the vdW interactions is found in κ -CuCN and κ -AgCN for which the absence of vdW interaction produces deviations of the computed unit cell parameters from the experimental ones at the level of 10% and 14%, respectively (Table 1 and Table 2). Hence, only the DFT calculations with the non-local vdW-DF functional can clarify the impact of disorder on the electronic ground state (see Section 3). Intriguingly enough, the vdW forces are almost irrelevant in κ -Cl (Table 3) thereby indicating that the latter structure is dominantly determined by the one-dimensional chains of anion Cu[N(CN) 2 ]Cl chemically bonded to the organic dimers.

5. Materials and Methods

All measurements were performed on high-quality single crystals grown by J. A. Schlueter, Argonne ( κ -(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl); K. Miyagawa and K. Kanoda, Tokyo University ( κ -(BEDT-TTF) 2 Cu 2 (CN) 3 ); and T. Hiramatsu, Y. Yoshida and G. Saito, Meijo University ( κ -(BEDT-TTF) 2 Ag 2 (CN) 3 ). Sample S1 of κ -Cl was of cube-like shape with an area of 0.4 × 0.5 mm 2 , and a thickness of 0.6 mm. Sample S2 of κ -Cl was plate-like with an a c -plane area of 0.8 × 0.9 mm 2 . and a thickness of 0.4 mm. Contacts for dc and dielectric measurements along the out-of-plane direction were made by applying conductive carbon paint directly to the surface of the sample [21,23]. DC resistivity was measured by a standard four contact technique. Temperature-dependent dielectric measurements were performed from room temperature down to 4.2 K (crystals S1 and S2 of κ -Cl) and 1.8 K (crystal of κ -CuCN) in the frequency range 40 Hz–10 MHz using Agilent 4294A and Keysight E4980AL-102 precision impedance analyzers. Care was taken to subtract the background capacitance of the set-up and sample holder, as well as to exclude any possible extrinsic effects due to sample preparation. Magnetic-field-dependent dielectric measurements of sample S2 of κ -Cl were performed in the magnetic field up to 5 T aligned along the out-of-plane crystallographic axis. Vibrational spectroscopy measurements were made on the thin side of the crystals by an IR microscope attached to a Bruker Vertex 80v Fourier-transform spectrometer with 1 cm − 1 resolution. Magnetic susceptibility of sample S2 of κ -Cl was performed by using a superconducting quantum interference device (SQUID) magnetometer and the data were corrected for the contribution of sample holder.
First principles calculations of the electronic structure was carried out by using the self-consistently implemented nonlocal van der Waals density functional (vdW-DF) [71,72,73] for correlation with optB88 for exchange [74]. The expansion in the plane waves was done with a cutoff energy of 700 eV. The Brillouin zone was sampled by 1 × 2 × 2 Monkhorst–Pack choice of k-points [75]. For the atomic coordinates, the experimental structure obtained by X-ray diffraction measurements was used. The structure was relaxed allowing for change of the volume and the shape of the unit cell together with atomic positions within, until the forces on atoms dropped below 1 meV/Å.
To ensure the chemical composition and exclude presence of any bromine in sample S2 of κ -Cl, we examined the crystal symmetry and composition [76]. The X-ray data were collected at 300 K with graphite monochromated Mo-K α radiation by ω -scans on an Oxford Diffraction KM4 XCALIBUR2 CCD diffractometer. Data acquisition and reduction was performed using the CrysAlisPro software package. The structure was solved by direct methods using SHELXS and the refinement was performed using SHELXL [77]. Energy Dispersive X-ray Spectra (EDXS) of sample S2 of κ -Cl, were acquired using Bruker detector within Tescan Vega III Easyprobe Scanning Electron Microscope (SEM).
Figure 13 shows the result of structural refinements of X-ray data within models containing either Cl (Figure 13, Left) or Br (Figure 13, Right). All atoms are drawn with 50% probability ellipsoids of thermal displacement. Importantly, the size of ellipsoids is convincingly small for Cl, while their size is unrealistically inflated if the presence of Br is assumed. In addition, the distance of 2.22 Å between Cu and neighboring Cl atom (represented by lines) is close to 2.26 Å which is typical for Cu-Cl bond, thus confirming the presence of Cl in the structure. Finally, the presence of Br larger than a few percent can be rejected since the corresponding bond in the case of Br would be longer and amount to 2.42 Å. Therefore, overall results of the structural refinements strongly indicate that the measured single crystal is κ -Cl and not κ -Br.
Figure 14 shows EDXS recorder at 300 K. Characteristic lines for C, N, S and Cu were observed with intensities that correspond to the κ -Cl material. On the other hand, no significant intensity of Br lines, expected if trace amount of Br is present, were observed.

6. Conclusions

Mott insulators κ -CuCN and κ -AgCN with the quantum spin liquid ground state exhibit high geometrical frustration in the organic dimers and in the non-conducting anions. Frustration in the anions arises as a combined effect due to the triangular coordination of copper/silver atoms and the bridging cyanide groups sitting in the inversion centers. The resulting disorder is mapped via anion-cation coupling onto the organic dimers. It gives rise to charge redistribution and a ground state which consists of several quasi-degenerate electronic states with reduced symmetry.
Remarkably, QSL is replaced by the antiferromagnetic ground state in κ -Cl where disorder due to inherent randomness in the anions and related symmetry reduction are absent. This result strongly suggests that the inherent randomness may be of a crucial importance for the prevalence of QSL over the antiferromagnetic ordering. In addition, van der Waals forces are found to contribute to inter-dimer interactions and the ground state structure in κ -CuCN and κ -AgCN, whereas their role appears insignificant in κ -Cl where chemical bonding dominantly determines the inter-dimer interactions besides the cation–anion coupling.
Anomalous dielectric peak shows up as an ubiquitous property independent on the nature of the ground state in the spin sector. It is now well established that, in κ -CuCN and κ -AgCN, this arises due to the cooperative motion of charged domain walls. On the other hand, in κ -Cl, some mysteries remain unresolved. The charged domain wall motion might also be in action in κ -Cl if the tendency to phase segregation is invoked due to canted spins in the antiferromagnetic state. Further experimental efforts are vital in search for structural inversion symmetry breaking in order to clarify the ferroelectricity in κ -Cl. The possibility that ferroelectricity is proximate to the magnetic and superconducting phases is most exciting and deserves more attention in future studies. Experiments such as X-ray diffraction and dielectric spectroscopy under pressure, on irradiated samples and in magnetic field are envisaged as tools to elucidate respective roles of correlations, randomness and charge-spin coupling in the phase diagram.

Author Contributions

DC and dielectric measurements and the data analysis were performed by Marko Pinteric, David Rivas Gongora, and Zeljko Rapljenovic, with help of Tomislav Ivek and Matija Culo. Dielectric measurements in magnetic field were designed by Tomislav Ivek and Matija Culo. Structural data analysis was performed by Ognjen Milat. DFT calculations were performed by Predrag Lazic. Vibrational experiments and data analysis were performed by Miriam Sanz Alonso, Andrej Pustogow, and Weiwu Li. Susceptibility experiments and data analysis were performed by Guilherme Gorgen Lesseux. The interpretation of the data and the text of the manuscript were drafted by Silvia Tomic and Tomislav Ivek. All authors contributed to the discussion and the final manuscript.

Acknowledgments

X-ray diffraction measurements by V. Stilinović and SEM measurements by J. Macan are gratefully acknowledged. We would like to thank V. Dobrosavljević, E. Canadell, R. T. Clay, P. Foury, V. Ilakovac, S. Mazumdar and J-P. Pouget for many enlightening discussions. This work has been supported by the Croatian Science Foundation project IP-2013-11-1011 and by the Deutsche Forschungsgemeinschaft (DFG).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kanoda, K.; Kato, R. Mott Physics in Organic Conductors with Triangular Lattices. Annu. Rev. Condens. Matter Phys. 2011, 2, 167. [Google Scholar] [CrossRef]
  2. Powell, B.J.; McKenzie, R.H. Quantum frustration in organic Mott insulators: From spin liquids to unconventional superconductors. Rep. Prog. Phys. 2011, 74, 056501. [Google Scholar] [CrossRef]
  3. Kandpal, H.C.; Opahle, I.; Zhang, Y.-Z.; Jeschke, H.O.; Valentí, R. Revision of Model Parameters for κ-Type Charge Transfer Salts: An Ab Initio Study. Phys. Rev. Lett. 2009, 103, 067004. [Google Scholar] [CrossRef] [PubMed]
  4. Jeschke, H.O.; de Souza, M.; Valentí, R.; Manna, R.S.; Lang, M.; Schlueter, J.A. Temperature dependence of structural and electronic properties of the spin-liquid candidate κ-(BEDT-TTF)2Cu2(CN)3. Phys. Rev. B 2012, 85, 035125. [Google Scholar] [CrossRef]
  5. Koretsune, T.; Hotta, C. Evaluating model parameters of the κ and β′-type Mott insulating organic solids. Phys. Rev. B 2014, 89, 045102. [Google Scholar] [CrossRef]
  6. Hiramatsu, T.; Yoshida, Y.; Saito, G.; Otsuka, A.; Yamochi, H.; Maesato, M.; Shimizu, Y.; Ito, H.; Nakamura, Y.; Kishida, H.; et al. Design and preparation of a new quantum spin liquid candidate, κ-(ET)2Ag2(CN)3, having a nearby superconductivity. Bull. Chem. Soc. Jpn. 2017, 90, 1073–1082. [Google Scholar] [CrossRef]
  7. Miyagawa, K.; Kawamoto, A.; Nakazawa, Y.; Kanoda, K. Antiferromagnetic Ordering and Spin Structure in the Organic Conductor, κ-(BEDT-TTF)2Cu[N(CN)2]Cl. Phys. Rev. Lett. 1995, 75, 1174. [Google Scholar] [CrossRef] [PubMed]
  8. Pinterić, M.; Miljak, M.; Biškup, N.; Milat, O.; Aviani, I.; Tomić, S.; Schweitzer, D.; Strunz, W.; Heinen, I. Magnetic anisotropy and low frequency dielectric response of weak ferromagnetic phase in κ-(BEDT-TTF)2Cu[N(CN)2]Cl. Eur. Phys. J. B 1999, 11, 217. [Google Scholar]
  9. Smith, D.F.; De Soto, S.M.; Slichter, C.P. Dzialoshinskii-Moriya interaction in the organic superconductor κ-(BEDT-TTF)2Cu[N(CN)2]Cl. Phys. Rev. B 2003, 68, 024512. [Google Scholar] [CrossRef]
  10. Shimizu, Y.; Miyagawa, K.; Kanoda, K.; Maesato, M.; Saito, G. Spin liquid state in an organic Mott insulator with a triangular lattice. Phys. Rev. Lett. 2003, 91, 107001. [Google Scholar] [CrossRef] [PubMed]
  11. Shimizu, Y.; Hiramatsu, T.; Maesato, M.; Otsuka, A.; Yamochi, H.; Ono, A.; Itoh, M.; Yoshida, M.; Takigawa, M.; Yoshida, Y.; et al. Pressure-tuned exchange coupling of a quantum spin liquid in the molecular triangular lattice κ-(BEDT-TTF)2Ag2(CN)3. Phys. Rev. Lett. 2016, 117, 107203. [Google Scholar] [CrossRef] [PubMed]
  12. Zhou, Y.; Kanoda, K.; Ng, T.-K. Quantum spin liquid states. Rev. Mod. Phys. 2017, 89, 025003. [Google Scholar] [CrossRef]
  13. Balents, L. Spin liquids in frustrated magnets. Nature 2010, 464, 199. [Google Scholar] [CrossRef] [PubMed]
  14. Savary, L.; Balents, L. Quantum spin liquids: A review. Rep. Prog. Phys. 2017, 80, 016502. [Google Scholar] [CrossRef] [PubMed]
  15. Anderson, W.P. Resonating valence bonds: A new kind of insulator? Mat. Res. Bull. 1973, 8, 153–160. [Google Scholar] [CrossRef]
  16. Klanjšek, M.; Zorko, A.; Žitko, R.; Mravlje, J.; Jagličić, Z.; Kumar Biswas, P.; Prelovšek, P.; Mihailovic, D.; Arčon, D. A high-temperature quantum spin liquid with polaron spins. Nat. Phys. 2017, 13, 1130–1134. [Google Scholar] [CrossRef]
  17. Shimozawa, M.; Hashimoto, K.; Ueda, A.; Suzuki, Y.; Sugii, K.; Yamada, S.; Imai, Y.; Kobayashi, R.; Itoh, K.; Iguchi, S.; et al. Quantum-disordered state of magnetic and electric dipoles in an organic Mott system. Nat. Commun. 2017, 8, 1821. [Google Scholar] [CrossRef] [PubMed]
  18. Lunkenheimer, P.; Müller, J.; Krohns, S.; Schrettle, F.; Loidl, A.; Hartmann, B.; Rommel, R.; de Souza, M.; Hotta, C.; Schlueter, J.A.; et al. Multiferroicity in an organic charge-transfer salt that is suggestive of electric-dipole-driven magnetism. Nat. Mater. 2012, 11, 755. [Google Scholar] [CrossRef] [PubMed]
  19. Tomić, S.; Pinterić, M.; Ivek, T.; Sedlmeier, K.; Beyer, R.; Wu, D.; Schlueter, J.A.; Schweitzer, D.; Dressel, M. Magnetic ordering and charge dynamics in κ-(BEDT-TTF)2Cu[N(CN)2]Cl. J. Phys. Condens. Matter 2013, 25, 436004. [Google Scholar] [CrossRef] [PubMed]
  20. Abdel-Jawad, M.; Terasaki, I.; Sasaki, T.; Yoneyama, N.; Kobayashi, N.; Uesu, Y.; Hotta, C. Anomalous dielectric response in the dimer Mott insulator κ-(BEDT-TTF)2Cu2(CN)3. Phys. Rev. B 2010, 82, 125119. [Google Scholar] [CrossRef]
  21. Pinterić, M.; Čulo, M.; Milat, O.; Basletić, M.; Korin-Hamzić, B.; Tafra, E.; Hamzić, A.; Ivek, T.; Peterseim, T.; Miyagawa, K.; et al. Anisotropic charge dynamics in the spin-Liquid candidate κ-(BEDT-TTF)2Cu2(CN)3. Phys. Rev. B 2014, 90, 195139. [Google Scholar] [CrossRef]
  22. Tomić, S.; Dressel, M. Ferroelectricity in molecular solids: A review of electrodynamic properties. Rep. Prog. Phys. 2015, 78, 096501. [Google Scholar] [CrossRef] [PubMed]
  23. Pinterić, M.; Lazić, P.; Pustogow, A.; Ivek, T.; Kuveždić, M.; Milat, O.; Gumhalter, B.; Basletić, M.; Čulo, M.; Korin-Hamzić, B.; et al. Anions effects on the electronic structure and electrodynamic properties of the Mott insulator κ-(BEDT-TTF)2Ag2(CN)3. Phys. Rev. B 2016, 94, 161105(R). [Google Scholar] [CrossRef]
  24. Dressel, M.; Lazić, P.; Pustogow, A.; Zhukova, E.; Gorshunov, B.; Schlueter, J.A.; Milat, O.; Gumhalter, B.; Tomić, S. Lattice vibrations of the charge-transfer salt κ-(BEDT-TTF)2Cu2(CN)3: Comprehensive explanation of the electrodynamic response in a spin-liquid compound. Phys. Rev. B 2016, 93, 081201(R). [Google Scholar] [CrossRef]
  25. Sedlmeier, K.; Elsässer, S.; Neubauer, D.; Beyer, R.; Wu, D.; Ivek, T.; Tomić, S.; Schlueter, J.A.; Dressel, M. Absence of charge order in the dimerized κ-phase BEDT-TTF salts. Phys. Rev. B 2012, 86, 245103. [Google Scholar] [CrossRef]
  26. Hotta, C. Quantum electric dipoles in spin-liqiod dimer Mott insulator κ-(ET)2Cu2(CN)3. Phys. Rev. B 2010, 82, 241104. [Google Scholar] [CrossRef]
  27. Hotta, C. Theories on Frustrated Electrons in Two-Dimensional Organic Solids. Crystals 2012, 2, 1155–1200. [Google Scholar] [CrossRef]
  28. Naka, M.; Ishihara, S. Electronic ferroelectricity in a dimer mott insulator. J. Phys. Soc. Jpn. 2010, 79, 063707. [Google Scholar] [CrossRef]
  29. Naka, M.; Ishihara, S. Collective charge excitation in a dimer Mott insulating system. J. Phys. Soc. Jpn. 2013, 82, 023701. [Google Scholar] [CrossRef]
  30. Naka, M.; Ishihara, S. Quantum melting of magnetic order in an organic dimer-Mott insulating system. Phys. Rev. B 2016, 93, 195114. [Google Scholar] [CrossRef]
  31. Li, H.; Clay, R.T.; Mazumdar, S. The paired-electron crystal in the two-dimensional frustrated quarter-filled band. J. Phys. Condens. Matter 2010, 22, 272201. [Google Scholar] [CrossRef] [PubMed]
  32. Dayal, S.; Clay, R.T.; Li, H.; Mazumdar, S. Paired electron crystal: Order from frustration in the quarter-filled band. Phys. Rev. B 2011, 83, 245106. [Google Scholar] [CrossRef]
  33. Gomi, H.; Ikenaga, M.; Hiragi, Y.; Segawa, D.; Takahashi, A.; Inagaki, T.J.; Aihara, M. Ferroelectric states induced by dimer lattice disorder in dimer Mott insulators. Phys. Rev. B 2013, 87, 195126. [Google Scholar] [CrossRef]
  34. Gomi, H.; Inagaki, T.J.; Takahashi, A. Ferroelectric charge order enhanced by magnetic frustration in dimer Mott insulators. Phys. Rev. B 2016, 93, 035105. [Google Scholar] [CrossRef]
  35. Huse, D.A.; Elser, V. Simple Variational Wave Functions for Two-Dimensional Heisenberg Spin-1/2 Antiferromagnets. Phys. Rev. Lett. 1988, 60, 2351–2354. [Google Scholar] [CrossRef] [PubMed]
  36. Bernu, B.; Lhuillier, C.; Pierre, L. Signature of Néel order in exact spectra of quantum antiferromagnets on finite lattices. Phys. Rev. Lett. 1992, 69, 2590–2593. [Google Scholar] [CrossRef] [PubMed]
  37. Capriotti, L.; Trumper, A.E.; Sorella, S. Long-range Néel order in the triangular Heisenberg model. Phys. Rev. Lett. 1999, 82, 3899–3902. [Google Scholar] [CrossRef]
  38. Furukawa, T.; Miyagawa, K.; Itou, T.; Ito, M.; Taniguchi, H.; Saito, M.; Iguchi, S.; Sasaki, T.; Kanoda, K. Quantum spin liquid emerging from antiferromagnetic order by introducing disorder. Phys. Rev. Lett. 2015, 115, 077001. [Google Scholar] [CrossRef] [PubMed]
  39. Lazić, P.; Pinterić, M.; Rivas Góngora, D.; Pustogow, A.; Treptow, K.; Ivek, T.; Milat, O.; Gumhalter, B.; Došlić, N.; Dressel, M.; et al. Importance of van der Waals interactions and cation-anion coupling in an organic quantum spin liquid. arXiv, 2017; arXiv:1710.01942. [Google Scholar]
  40. Kornelsen, K.; Eldridge, J.E.; Wang, H.H.; Charlier, H.A.; Williams, J.M. Infrared study of the metal-insulator transition in the organic conductor κ-(BEDT-TTF)2Cu[N(CN)2]Cl. Solid State Commun. 1992, 4, 343–349. [Google Scholar] [CrossRef]
  41. Williams, J.M.; Kini, A.M.; Wang, H.H.; Carlson, K.D.; Geiser, U.; Montgomery, L.K.; Pyrka, G.J.; Watkins, D.M.; Kommers, J.M.; Boryschuk, S.J.; et al. From Semiconductor-Semiconductor Transition (42 K) to the Highest-Tc Organic Superconductor, κ-(BEDT-TTF)2Cu[N(CN)2]Cl (Tc = 12.5 K). Inorg. Chem. 1990, 29, 3272–3274. [Google Scholar] [CrossRef]
  42. Sushko, Y.V.; Andres, K. Superconducting Meissner effect under hydrostatic pressure in the ambient-pressure semiconductor κ-(BEDT-TTF)2Cu[N(CN)2]Cl, where BEDT-TTF is bis(ethylenedithio)tetrathiafulvalene. Phys. Rev. B 1993, 47, 330–333. [Google Scholar] [CrossRef]
  43. Lefebvre, S.; Wzietek, P.; Brown, S.; Bourbonnais, C.; Jérome, D.; Mézière, C.; Fourmigué, M.; Batail, P. Mott Transition, Antiferromagnetism, and Unconventional Superconductivity in Layered Organic Superconductors. Phys. Rev. Lett. 2000, 85, 5420–5423. [Google Scholar] [CrossRef] [PubMed]
  44. Dressel, M.; Gruener, G.; Carlson, K.D.; Wang, H.H.; Williams, J.M. Studies of the Microwave Resistivity of κ-(BEDT-TTF)2Cu[N(CN)2]Cl. Synth. Metals 1995, 70, 927–928. [Google Scholar] [CrossRef]
  45. Sushko, Y.V.; Ito, H.; Ishiguro, T.; Horiuchi, S.; Saito, G. Reentrant superconductivity in κ-(BEDT-TTF)2Cu[N(CN)2]Cl and its pressure phase diagram. Solid State Commun. 1993, 87, 997–1000. [Google Scholar] [CrossRef]
  46. Clay, R.T.; Mazumdar, S. From charge- and spin-ordering to superconductivity in the organic charge-transfer solids. arXiv, 2018; arXiv:1802.01551. [Google Scholar]
  47. Hiramatsu, T.; Yoshida, Y.; Saito, G.; Otsuka, A.; Yamochi, H.; Maesato, M.; Shimizu, Y.; Ito, H.; Kishida, H. Quantum spin liquid: Design of a quantum spin liquid next to a superconducting state based on a dimer-type ET Mott insulator. J. Mater. Chem. C 2015, 3, 1378. [Google Scholar] [CrossRef]
  48. Dressel, M.; Drichko, N. Optical Properties of Two-Dimensional Organic Conductors: Signatures of Charge Ordering and Correlation Effects. Chem. Rev 2004, 104, 5869. [Google Scholar] [CrossRef] [PubMed]
  49. Yamamoto, T.; Uruichi, M.; Yamamoto, K.; Yakushi, K.; Kawamoto, A.; Taniguchi, H. Examination of the Charge-Sensitive Vibrational Modes in Bis(ethylenedithio)tetrathiafulvalene. J. Phys. Chem. B 2005, 109, 15226. [Google Scholar] [CrossRef] [PubMed]
  50. Girlando, A. Charge Sensitive Vibrations and Electron-Molecular Vibration Coupling in Bis(ethylenedithio)- tetrathiafulvalene (BEDT-TTF). J. Phys. Chem. C 2011, 115, 19371. [Google Scholar] [CrossRef]
  51. Lang, M.; Lunkenheimer, P.; Müller, J.; Loidl, A.; Hartmann, B.; Hoang, N.H.; Gati, E.; Schubert, H.; Schlueter, J.A. Multiferroicity in the Mott insulating charge-transfer salt κ-(BEDT-TTF)2Cu[N(CN)2]Cl. IEEE Trans. Magn. 2014, 50, 2700107. [Google Scholar]
  52. Kawamoto, A.; Honma, Y.; Kumagai, K.I. Electron localization in the strongly correlated organic system κ-(BEDT-TTF)2 probed with nuclear magnetic resonance 13C-NMR. Phys. Rev. B 2004, 70, 060510(R). [Google Scholar] [CrossRef]
  53. Čulo, M.; Tafra, E.; Basletić, M.; Tomić, S.; Hamzić, A.; Korin-Hamzić, B.; Dressel, M.; Schlueter, J.A. Two-dimensional variable range hopping in the spin-liquid candidate κ-(BEDT-TTF)2Cu2(CN)3. Phys. B 2015, 460, 208–210. [Google Scholar] [CrossRef]
  54. Nad, F.; Monceau, P.; Kaboub, L.; Fabre, J.M. Divergence of the relaxation time in the vicinity of the ferroelectric charge-ordered phase transition in (TMTTF)2AsF6. Europhys. Lett. 2006, 73, 567–573. [Google Scholar] [CrossRef]
  55. Cava, R.J.; Schneemeyer, L.F.; Fleming, R.M.; Littlewood, P.B.; Rietman, E.A. Effects of impurities on the dielectric response of the charge-density wave in K0.3MoO3. Phys. Rev. B 1985, 32, 4088–4096. [Google Scholar] [CrossRef]
  56. Ngai, K.L. Relaxation and Diffusion in Complex Systems; Springer, New York, NY, USA, 2011; ISBN 978-1-4419-7649-9. [Google Scholar]
  57. Poirier, M.; Parent, S.; Côté, A.; Miyagawa, K.; Kanoda, K.; Shimizu, Y. Magnetodielectric effects and spin-charge coupling in the spin-liquid candidate κ-(BEDT-TTF)2Cu2(CN)3. Phys. Rev. B 2012, 85, 134444. [Google Scholar] [CrossRef]
  58. Manna, R.S.; Hartmann, S.; Gati, E.; Schlueter, J.A.; de Souza, M.; Lang, M. Low-Temperature Lattice Effects in the Spin-Liquid Candidate κ-(BEDT-TTF)2Cu2(CN)3. Crystals 2018, 8, 87. [Google Scholar] [CrossRef]
  59. Cross, L.E. Relaxor ferroelectrics. Ferroelectrics 1987, 76, 241–267. [Google Scholar] [CrossRef]
  60. Itoh, K.; Itoh, H.; Naka, M.; Saito, S.; Hosako, I.; Yoneyama, N.; Ishihara, S.; Sasaki, T.; Iwai, S. Collective Excitation of an Electric Dipole on a Molecular Dimer in an Organic Dimer-Mott Insulator. Phys. Rev. Lett. 2013, 110, 106401. [Google Scholar] [CrossRef] [PubMed]
  61. Geiser, U.; Wang, H.H.; Carlson, K.D.; Williams, J.M.; Charlier, H.A.; Heindl, J.E.; Yaconi, G.A.; Love, B.J.; Lathrop, M.W.; Schirber, J.E.; et al. Superconductivity at 2.8 K and 1.5 kbar in κ-(BEDT-TTF)2Cu2(CN)3: First Organic Superconductor Containing a Polymeric Copper Cyanide Anion. Inorg. Chem. 1991, 30, 2586. [Google Scholar] [CrossRef]
  62. Fukuyama, H.; Kishine, J.; Ogata, M. Energy Landscape of Charge Excitations in the Boundary Region between Dimer–Mott and Charge Ordered States in Molecular Solids. J. Phys. Soc. Jpn. 2017, 86, 123706. [Google Scholar] [CrossRef]
  63. Foury-Leylekian, P.; Ilakovac, V.; Balédent, V.; Fertey, P.; Arakcheeva, A.; Milat, O.; Petermann, D.; Guillier, G.; Miyagawa, K.; Kanoda, K.; et al. κ-(BEDT-TTF)2Cu2(CN)3 spin liquid: Beyond the average structure. Crystals 2018, 8, 158. [Google Scholar] [CrossRef]
  64. Nakamura, K.; Yoshimoto, Y.; Imada, M. Ab initio two-dimensional multiband low-energy models of EtMe3Sb[Pd(dmit)2]2 and κ-(BEDT-TTF)2Cu2(CN)3 with comparisons to single-band models. Phys. Rev. B 2012, 86, 205117. [Google Scholar] [CrossRef]
  65. Scriven, E.P.; Powell, B.J. Geometrical Frustration in the Spin Liquid β′-EtMe3Sb[Pd(dmit)2]2 and the Valence-Bond Solid Me3EtP[Pd(dmit)2]2. Phys. Rev. Lett. 2012, 109, 097206. [Google Scholar] [CrossRef] [PubMed]
  66. Tsumuraya, T.; Seo, H.; Tsuchiizu, M.; Kato, R.; Miyazaki, T. Cation Dependence of the Electronic States in Molecular Triangular Lattice System β′-X[Pd(dmit)2]2. J. Phys. Soc. Jpn. 2013, 82, 033709. [Google Scholar] [CrossRef]
  67. Jacko, A.C.; Tocchio, L.F.; Jeschke, H.O.; Valenti, R. Importance of anisotropy in the spin-liquid candidate EtMe3Sb[Pd(dmit)2]2. Phys. Rev. B 2013, 88, 155139. [Google Scholar] [CrossRef]
  68. Alemany, P.; Pouget, J.-P.; Canadell, E. Essential role of anions in the charge ordering transition of α-(BEDT-TTF)2I3. Phys. Rev. B 2012, 85, 195118. [Google Scholar] [CrossRef]
  69. Alemany, P.; Pouget, J.-P.; Canadell, E. Electronic structure and anion ordering in (TMTSF)2ClO4 and (TMTSF)2NO3: A first-principles study. Phys. Rev. B 2014, 89, 155124. [Google Scholar] [CrossRef]
  70. Alemany, P.; Pouget, J.-P.; Canadell, E. Structural and electronic control of the metal to insulator transition and local orderings in the θ-(BEDT-TTF)2X organic conductors. J. Phys. Condens. Matter 2015, 27, 465702. [Google Scholar] [CrossRef] [PubMed]
  71. Dion, M.; Rydber, H.; Schröder, E.; Langreth, D.C.; Lundqvist, B.I. Van derWaals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. [Google Scholar] [CrossRef] [PubMed]
  72. Román-Pérez, G.; Soler, J.M. Efficient Implementation of a van derWaals Density Functional: Application to Double-Wall Carbon Nanotubes. Phys. Rev. Lett. 2009, 103, 096102. [Google Scholar] [CrossRef] [PubMed]
  73. Klimeš, J.; Bowler, D.R.; Michaelides, A. Van der Waals density functionals applied to solids. Phys. Rev. B 2011, 83, 195131. [Google Scholar] [CrossRef]
  74. Mittendorfer, F.; Garhofer, A.; Redinger, J.; Klimeš, J.; Harl, J.; Kresse, G. Graphene on Ni(111): Strong interaction and weak adsorption. Phys. Rev. B 2011, 84, 201401. [Google Scholar] [CrossRef]
  75. Monkhorst, H.J.; Pack, J.D. Special points for Brillonin-zone integrations. Phys. Rev. B 1976, 13, 5188. [Google Scholar] [CrossRef]
  76. Yasin, S.; Dumm, M.; Salameh, B.; Batail, P.; Meziere, C.; Dressel, M. Transport studies at the Mott transition of the two-dimensional organic metal κ-(BEDT-TTF)2Cu[N(CN)2]Brx Cl1−x. Eur. Phys. J. B 2011, 79, 383–390. [Google Scholar] [CrossRef] [Green Version]
  77. Sheldrick, G.M. A short history of SHELX. Acta Cryst. A 2008, 64, 112–122. [Google Scholar] [CrossRef] [PubMed]
Crystals 08 00190 g001 550
Figure 1. The crystal structure of the three layered organic charge-transfer salts. The lines mark the unit cell. Carbon, sulfur and hydrogen atoms of the BEDT-TTF molecule are colored in dark gray, yellow and light gray, respectively. In the anion network, chlorine, cooper, silver, carbon and nitrogen are colored in green, red, pink, dark grey and violet, respectively. (Left) Unit cell of κ -Cl. The space group is P n m a . In each of two organic layers, related by mirror symmetry, all four BEDT-TTF molecules are equivalent. (Middle,Right) Extended unit cell of κ -CuCN and κ -AgCN, respectively. The space group is commonly solved in P 2 1 / c in which all four BEDT-TTF molecules are equivalent.
Figure 1. The crystal structure of the three layered organic charge-transfer salts. The lines mark the unit cell. Carbon, sulfur and hydrogen atoms of the BEDT-TTF molecule are colored in dark gray, yellow and light gray, respectively. In the anion network, chlorine, cooper, silver, carbon and nitrogen are colored in green, red, pink, dark grey and violet, respectively. (Left) Unit cell of κ -Cl. The space group is P n m a . In each of two organic layers, related by mirror symmetry, all four BEDT-TTF molecules are equivalent. (Middle,Right) Extended unit cell of κ -CuCN and κ -AgCN, respectively. The space group is commonly solved in P 2 1 / c in which all four BEDT-TTF molecules are equivalent.
Crystals 08 00190 g001
Crystals 08 00190 g002 550
Figure 2. View of BEDT-TTF dimers arranged in anisotropic triangles in the three layered organic charge-transfer salts: (Left) κ -Cl, projection of one of the layers is shown along the [110] direction; (Middle) κ -CuCN, projected along the [100] direction; and (Right) κ -AgCN, projected along the [10-1] direction. Neighboring dimers are rotated by 90 ∘ with respect to each other. The interdimer transfer integrals are denoted by t and t ′ , and the intradimer transfer integral by t d . The ratio t ′ / t measures the degree of frustration. t ′ / t ≈ 0.5 indicates medium frustration for κ -Cl, while t ′ / t ≈ 0.85 indicates high frustration for κ -CuCN and κ -AgCN.
Figure 2. View of BEDT-TTF dimers arranged in anisotropic triangles in the three layered organic charge-transfer salts: (Left) κ -Cl, projection of one of the layers is shown along the [110] direction; (Middle) κ -CuCN, projected along the [100] direction; and (Right) κ -AgCN, projected along the [10-1] direction. Neighboring dimers are rotated by 90 ∘ with respect to each other. The interdimer transfer integrals are denoted by t and t ′ , and the intradimer transfer integral by t d . The ratio t ′ / t measures the degree of frustration. t ′ / t ≈ 0.5 indicates medium frustration for κ -Cl, while t ′ / t ≈ 0.85 indicates high frustration for κ -CuCN and κ -AgCN.
Crystals 08 00190 g002
Crystals 08 00190 g003 550
Figure 3. (Left) Distinct behaviors of the resistivity along the b axis below 100 K for three single crystals of κ -Cl: S1 (red), S2 (blue) and S3 (black). Data for S3 are after [18]. The inset displays the resistivity behavior in the whole temperature range demonstrating a weak activated behavior similar for all three samples. (Right) Logarithmic resistivity derivative versus temperature.
Figure 3. (Left) Distinct behaviors of the resistivity along the b axis below 100 K for three single crystals of κ -Cl: S1 (red), S2 (blue) and S3 (black). Data for S3 are after [18]. The inset displays the resistivity behavior in the whole temperature range demonstrating a weak activated behavior similar for all three samples. (Right) Logarithmic resistivity derivative versus temperature.
Crystals 08 00190 g003
Crystals 08 00190 g004 550
Figure 4. Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ b for three representative single crystals of κ -Cl: S1 (Left); S2 (Middle); and S3, after [18] (Right panel) The full lines in (Right) are guides for the eye.
Figure 4. Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ b for three representative single crystals of κ -Cl: S1 (Left); S2 (Middle); and S3, after [18] (Right panel) The full lines in (Right) are guides for the eye.
Crystals 08 00190 g004
Crystals 08 00190 g005 550
Figure 5. Magnetic field influence on sample S2 of κ -Cl. (Left) Conductivity σ measured at 500 Hz for E ∥ b versus temperature measured at different magnetic field strengths ( H ∥ b ). At H = 0 T there is a rise in conductivity below 12 K, while the application of magnetic field shifts the onset of conductivity rise to lower temperatures, thus indicating a superconducting phase transition at T SC ≈ 12 K. (Right) Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ b , and measured at different magnetic field strengths ( H ∥ b ). The effect of the magnetic field can be seen only below 12 K. Increasing magnetic field lowers the peak amplitude.
Figure 5. Magnetic field influence on sample S2 of κ -Cl. (Left) Conductivity σ measured at 500 Hz for E ∥ b versus temperature measured at different magnetic field strengths ( H ∥ b ). At H = 0 T there is a rise in conductivity below 12 K, while the application of magnetic field shifts the onset of conductivity rise to lower temperatures, thus indicating a superconducting phase transition at T SC ≈ 12 K. (Right) Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ b , and measured at different magnetic field strengths ( H ∥ b ). The effect of the magnetic field can be seen only below 12 K. Increasing magnetic field lowers the peak amplitude.
Crystals 08 00190 g005
Crystals 08 00190 g006 550
Figure 6. Susceptibility χ of sample S2 of κ -Cl as a function of temperature measured in a magnetic field of 500 Oe applied H ∥ a c -plane.
Figure 6. Susceptibility χ of sample S2 of κ -Cl as a function of temperature measured in a magnetic field of 500 Oe applied H ∥ a c -plane.
Crystals 08 00190 g006
Crystals 08 00190 g007 550
Figure 7. (Left) Temperature evolution of the intramolecular vibration ν 27 in sample S2 of κ -Cl. (Right) Temperature dependence of the resonance frequency and damping obtained by Fano function fit of the mode.
Figure 7. (Left) Temperature evolution of the intramolecular vibration ν 27 in sample S2 of κ -Cl. (Right) Temperature dependence of the resonance frequency and damping obtained by Fano function fit of the mode.
Crystals 08 00190 g007
Crystals 08 00190 g008 550
Figure 8. View of the anion layer of κ -Cl in the a c plane projected along the b-axis. Copper is colored in red, and chlorine in green. Carbon and nitrogen are colored in blue and orange, respectively. The unit cell is marked as a rectangle.
Figure 8. View of the anion layer of κ -Cl in the a c plane projected along the b-axis. Copper is colored in red, and chlorine in green. Carbon and nitrogen are colored in blue and orange, respectively. The unit cell is marked as a rectangle.
Crystals 08 00190 g008
Crystals 08 00190 g009 550
Figure 9. Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ a * for single crystals of: κ -CuCN (Left); and κ -AgCN (Right).
Figure 9. Real part of the dielectric function ε ′ as a function of temperature with the ac electric field applied E ∥ a * for single crystals of: κ -CuCN (Left); and κ -AgCN (Right).
Crystals 08 00190 g009
Crystals 08 00190 g010 550
Figure 10. (Left) Temperature evolution of the intramolecular vibration ν 27 of κ -CuCN. (Right) Temperature dependence of the resonance frequency and damping obtained by Fano function fit of the mode.
Figure 10. (Left) Temperature evolution of the intramolecular vibration ν 27 of κ -CuCN. (Right) Temperature dependence of the resonance frequency and damping obtained by Fano function fit of the mode.
Crystals 08 00190 g010
Crystals 08 00190 g011 550
Figure 11. (Left) Temperature evolution of the intramolecular vibration ν 27 of κ -AgCN. The two Fano functions required to fit the spectra indicate two crystallographically distinct sites. (Right) Temperature dependence of the resonance frequency and damping obtained by the fit using two Fano functions.
Figure 11. (Left) Temperature evolution of the intramolecular vibration ν 27 of κ -AgCN. The two Fano functions required to fit the spectra indicate two crystallographically distinct sites. (Right) Temperature dependence of the resonance frequency and damping obtained by the fit using two Fano functions.
Crystals 08 00190 g011
Crystals 08 00190 g012 550
Figure 12. View of the anion network of κ -CuCN (Left) and κ -AgCN (Right) in the b c plane projected along the a * -axis. Cooper and silver are colored in red and violet, respectively. Carbon and nitrogen of ordered CN − groups are colored in blue and orange, while they are colored in black for CN − groups located at inversion centers. The unit cell is marked as a rectangle.
Figure 12. View of the anion network of κ -CuCN (Left) and κ -AgCN (Right) in the b c plane projected along the a * -axis. Cooper and silver are colored in red and violet, respectively. Carbon and nitrogen of ordered CN − groups are colored in blue and orange, while they are colored in black for CN − groups located at inversion centers. The unit cell is marked as a rectangle.
Crystals 08 00190 g012
Crystals 08 00190 g013 550
Figure 13. Perspective view of zigzag polymeric anion chains as obtained by structural refinements of X-ray data assuming either: κ -Cl (Left); or κ -Br (Right). View shows the ac plane. Copper, carbon and nitrogen are shown in red, blue and orange, respectively. Chlorine and bromine ligands on each copper atom are colored in dark and light green, respectively. All atoms are drawn with 50% probability ellipsoids of anisotropic thermal displacement parameters. Unit cell is marked as a rectangle.
Figure 13. Perspective view of zigzag polymeric anion chains as obtained by structural refinements of X-ray data assuming either: κ -Cl (Left); or κ -Br (Right). View shows the ac plane. Copper, carbon and nitrogen are shown in red, blue and orange, respectively. Chlorine and bromine ligands on each copper atom are colored in dark and light green, respectively. All atoms are drawn with 50% probability ellipsoids of anisotropic thermal displacement parameters. Unit cell is marked as a rectangle.
Crystals 08 00190 g013
Crystals 08 00190 g014 550
Figure 14. Energy Dispersive X-ray Spectra (EDXS) recorded at 300 K. Only characteristic lines for C, N, S and Cu, and in particular for ClK α at 2.59 keV and for ClK β at 2.83 keV were observed with intensities that correspond to the chemical formula of κ -Cl. On the other hand, no significant intensity for lines BrK α 1 and BrK α 2 at 11.92 keV, for BrK β at 13.29 keV, and for BrL α at 1.49 keV were observed, which would be expected in the presence of Br.
Figure 14. Energy Dispersive X-ray Spectra (EDXS) recorded at 300 K. Only characteristic lines for C, N, S and Cu, and in particular for ClK α at 2.59 keV and for ClK β at 2.83 keV were observed with intensities that correspond to the chemical formula of κ -Cl. On the other hand, no significant intensity for lines BrK α 1 and BrK α 2 at 11.92 keV, for BrK β at 13.29 keV, and for BrL α at 1.49 keV were observed, which would be expected in the presence of Br.
Crystals 08 00190 g014
Table
Table 1. Unit cell parameters of κ -CuCN obtained from X-ray diffraction measurements at 100 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Table 1. Unit cell parameters of κ -CuCN obtained from X-ray diffraction measurements at 100 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Unit Cell ParametersExpCalc:PBECalc:vdW-DF
a15.9644 Å16.9706 Å (+6.3%)16.1809 Å (+1.4%)
b8.5618 Å8.75745 Å (+2.3%)8.62586 Å (+0.8%)
c13.2662 Å13.6497 Å (+2.9%)13.3523 Å (+0.7%)
α 90.000 ∘ 89.9968 ∘ 90.0037 ∘
β 114.067 ∘ 116.409 ∘ 114.088 ∘
γ 90.000 ∘ 90.0037 ∘ 89.9994 ∘
V1655.65 Å 3 1816.90 Å 3 (+9.7%)1701.35 Å 3 (+2.8%)
Table
Table 2. Unit cell parameters of κ -AgCN obtained from X-ray diffraction measurements at 150 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Table 2. Unit cell parameters of κ -AgCN obtained from X-ray diffraction measurements at 150 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Unit Cell ParametersExpCalc:PBECalc:vdW-DF
a14.969900  Å15.6693 Å (+4.7%)15.1161 Å (+0.98%)
b8.656500 Å9.04653 Å (+4.5%)8.72513 Å (+0.79%)
c13.216900  Å13.8238 Å (+4.6%)13.3487 Å (+1.0%)
α 90.000000 ∘ 90.0000 ∘ 90.0000 ∘
β 91.389000 ∘ 94.6823 ∘ 91.3952 ∘
γ 90.000000 ∘ 90.0000 ∘ 90.0000 ∘
V1712.23 Å 3 1953.03 Å 3 (+14.1%)1760.03 Å 3 (+2.8%)
Table
Table 3. Unit cell parameters of κ -Cl obtained from X-ray diffraction measurements at 100 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Table 3. Unit cell parameters of κ -Cl obtained from X-ray diffraction measurements at 100 K (left column), ab initio calculations based solely on PBE functional (central column) and on vdW-DF functional (right column). Relative deviations from experimental values are given in parentheses.
Unit Cell ParametersExpCalc:PBECalc:vdW-DF
a12.885200 Å12.9769 Å (+0.7%)12.8334 Å (−0.4%)
b29.575899 Å29.6637 Å (+0.3%)29.5323 Å (−0.2%)
c8.416100 Å8.45028 Å (+0.4%)8.23316 Å (−2.2%)
α 90.000000 ∘ 90.0000 ∘ 90.0000 ∘
β 90.000000 ∘ 90.0000 ∘ 90.0000 ∘
γ 90.000000 ∘ 90.0000 ∘ 90.0000 ∘
V3207.303 Å 3 3252.875 Å 3 (+1.4%)3120.366 Å 3 ( − 2 . 7 % )

Share and Cite

MDPI and ACS Style

Pinterić, M.; Rivas Góngora, D.; Rapljenović, Ž.; Ivek, T.; Čulo, M.; Korin-Hamzić, B.; Milat, O.; Gumhalter, B.; Lazić, P.; Sanz Alonso, M.; et al. Electrodynamics in Organic Dimer Insulators Close to Mott Critical Point. Crystals 2018, 8, 190. https://doi.org/10.3390/cryst8050190

AMA Style

Pinterić M, Rivas Góngora D, Rapljenović Ž, Ivek T, Čulo M, Korin-Hamzić B, Milat O, Gumhalter B, Lazić P, Sanz Alonso M, et al. Electrodynamics in Organic Dimer Insulators Close to Mott Critical Point. Crystals. 2018; 8(5):190. https://doi.org/10.3390/cryst8050190

Chicago/Turabian Style

Pinterić, Marko, David Rivas Góngora, Željko Rapljenović, Tomislav Ivek, Matija Čulo, Bojana Korin-Hamzić, Ognjen Milat, Branko Gumhalter, Predrag Lazić, Miriam Sanz Alonso, and et al. 2018. "Electrodynamics in Organic Dimer Insulators Close to Mott Critical Point" Crystals 8, no. 5: 190. https://doi.org/10.3390/cryst8050190

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop