Ferromagnetism out of charge fluctuation of strongly correlated electrons in κ-(BEDT-TTF)2Hg(SCN)2Br

We perform magnetic susceptibility and magnetic torque measurements on the organic κ-(BEDT-TTF)2Hg(SCN)2Br, which is recently suggested to host an exotic quantum dipole-liquid in its low-temperature insulating phase. Below the metal-insulator (MI) transition temperature, the magnetic susceptibility follows a Curie–Weiss law with a positive Curie–Weiss temperature, and a particular M∝H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M\propto \sqrt{H}$$\end{document} curve is observed. The emergent ferromagnetically interacting spins amount to about 1/6 of the full spin moment of localized charges. Taking account of the possible inhomogeneous quasi-charge-order that forms a dipole-liquid, we construct a model of antiferromagnetically interacting spin chains in two adjacent charge-ordered domains, which are coupled via fluctuating charges on a Mott-dimer at the boundary. We find that the charge fluctuations can draw a weak ferromagnetic moment out of the spin singlet domains.


INTRODUCTION
Typical phase transitions in condensed matter accompany either universal critical singularities or the competitions between two different orderings. The former is easily converted to the latter when additional degrees of freedom become relevant. In reality, there often appear intermediate situations where the interplay of several degrees of freedom affects the nature of the growth of correlations and low-lying excitations. In such cases, the phase transitions at low temperatures can be easily masked, and the order parameters suffer intrinsic inhomogeneities. Historical examples are the dynamically disordered charge stripes in high-T c superconducting cuprates 1,2 , and the orbital disorders that trigger the colossal magnetoresistance in manganites 3 . The complexity of dealing with multiple correlated degrees of freedom such as charge, spin, orbital, and lattice often makes it difficult to pin down their dominant mechanism.
Organic κ-(BEDT-TTF) 2 X materials ( Fig. 1a, b) become an ideal playground to study such an issue in a simpler setup. These materials form quarter-filled two-dimensional strongly correlated electronic systems, where molecular dimer (BEDT-TTF) 2 connected by a large transfer integral t d (see Fig. 1b) serves as a lattice site of a Mott insulator by hosting one charge per dimer. In these Mott insulators, the spin-1/2 interact antiferromagnetically as shown Fig. 1c and a quantum spin liquid phase is observed in κ-(BEDT-TTF) 2 Cu 2 (CN) 3 (abbreviated as κ-CN) as well as a typical antiferromagnetism in κ-(BEDT-TTF) 2 Cu[N(CN) 2 ]Cl (κ-Cu-Cl) 4,5 . If the degree of dimerization, namely the ratio of t d to other interdimer transfer integrals (t B , t p , and t q ), is weakened, the charges would no longer stay on the dimer-orbital, but rather localize to one side of the dimerized molecules to gain the intermolecular Coulomb interaction V ij , resulting in a charge-ordered phase. A charge degree of freedom enclosed in the dimer is interpreted as quantum electric dipole 6,7 , which is detected by anomalous frequency-dependence of dielectricity in many materials including κ-CN 7 , κ-Cu-Cl 8 , and EtMe 3 Sb[Pd(dmit) 2 ] 2 9,10 . In that context, the dimer Mott and charge-ordered phases are interpreted as para and ferroelectricity 6,11 , separated by a typical Ising type secondorder phase transition (see the phase diagram in Fig. 1d). However, when the universal criticality of dipoles couples to magnetism or lattice degrees of freedom, this transition can be masked and some inhomogeneous phases may emerge 12 . Indeed, the subtleties of the transition are recently disclosed by the fresh members of this family, κ-(BEDT-TTF) 2 Hg(SCN) 2 Br (κ-Hg-Br) and κ-(BEDT-TTF) 2 Hg(SCN) 2 Cl (κ-Hg-Cl), which have a relatively weak dimerization 13 and fill the empty region of materials parameter space. In contrast to a simple Mott insulator, which shows a crossover from the high-temperature metallic regime 4 , these compounds show an abrupt increase of resistivity at the metalinsulator (MI) transition 14 . Raman spectroscopy reveals a distinct charge order in κ-Hg-Cl in the temperature range 15-30 K 15 whereas κ-Hg-Br does not show any sign of regular charge ordering down to the lowest temperature 16 . The absence of magnetic order in κ-Hg-Br is also shown by the specific heat measurements down to 100 mK 16 . A picture of "quantum dipole liquid" is provided as an interpretation to the latter intriguing phase 16 , possibly consisting of dynamical charge-ordered domains enclosing electric dipole moments maximally amounting to 0.1e per dimer.
We report the experimental evidence of intrinsic ferromagnetic exchange interactions emerging in the clean bulk crystal of κ-Hg-Br, indicated by a positive Curie-Weiss temperature of Θ CW = 16 K, where about 1/6 of the full spins of localized charges contribute. We find that the M-H curve at low temperature follows M / ffiffiffi ffi H p , showing a very rapid onset with a small field. Although the square-root onset of the M-H curve is well-known for a gapped quasi-one-dimensional quantum magnet near the critical field 17,18 , it is qualitatively different from the present gapless M / ffiffiffi ffi H p that continues up to a large field. It does not resemble any of the M-H profile of the magnetism known so far, such as the H-linear antiferromagnetic magnetization or the paramagnetic Brillouin curve. Such robust ferromagnetic Curie-Weiss law just above the antiferromagnetic singlet ground state can be scarcely found in nature, except for those originating from magnetic impurities or a spin glass, both of which are excluded in the present case by the lack of remnant filed or hysteresis in magnetization. Since no existing theory on bulk magnetism both for the localized spins or itinerant electrons can be applied, we construct a synergetic quantum-spin model that includes the effect of charge fluctuation. The starting point is the low-temperature inhomogeneous state of charges that appear by masking the phase transition in Fig. 1d. We take account of the already existing idea of a short charge correlation length and the robustly remaining charge fluctuation at the simplest level 16,19 . The model represents spins on two charge-ordered domains which couple to dimer-spins carried by fluctuating charges at the domain boundary, and successfully shows how ferromagnetic behavior can originate from the charge fluctuation. The theory thus explains the properties disclosed by the magnetic susceptibility and torque measurements.  Fig. 2b) gives a positive Θ CW~1 6 K with the Curie constant C = 0.060 emu K mol −1 . This positive Θ CW provides strong evidence of a ferromagnetic interaction between spins. The C value shows that~1/6 of the total spins contribute to the Curie-Weiss paramagnetism with the ferromagnetic interaction. The 1/6-concentration is intrinsic to the ferromagnetic behavior, as χ M (T) does not depend on the measured field strength below 5 T in this temperature range (see Supplementary Fig. 2 in Supplemental Information (SI)). Similar χ M (T) was observed previously 20 . However, our data shows the ferromagnetic Θ CW more clearly in a wider temperature range (see Section A in SI for a comparison). To the best of our knowledge, this compound is the first to show such ferromagnetic behavior in a family of organic Mott insulators κ-(BEDT-TTF) 2 X and X[Pd (dmit) 2 ] 2 . If one interprets this χ as the one from the ferromagnetic Heisenberg chain 21 , the effective ferromagnetic coupling constant is evaluated as J F~ΘCW /Θ = 53 K, with Θ = 0.3036.

Magnetization measurements
The particular ferromagnetic behavior is also found in the field dependence of the magnetization M. The linear M-H curve at high temperatures becomes nonlinear at 2.0 K as shown in Fig. 2c (see also Supplementary Fig. 2 in SI for the temperature dependence). Below 2 T, M increases faster than the Brillouin curve (the dashed line in Fig. 2c) 13 ). c Schematic illustration of the dimer Mott insulator and charge order. Charges on Mott-dimers (red circles) carry spin-1/2, and the exchange interactions between them J AF / t 2 B (vertical) and J AF / ðt p À t q Þ 2 (diagonal directions) form an antiferromagnetic triangular lattice (red lines). In the charge-ordered state, the charges localized on one side of the dimer (green circles) form an antiferromagnetic quasi-one-dimensional spin-1/2 chain of J AF / t 2 q (green line). Inset: the Mott-dimer state is the linear combination of charge located on the left/right molecule, supported by large t d . The charge order keeps the charge on one side of the dimerized pair to avoid the inter-site (inter-dimer) Coulomb interaction V ij = V B , V p , V q (indices follow those of t ij ). d Schematic phase diagram (Hotta, C., Yoshida, T. & Harada, K. unpublished) of the present system for low energy effective model of charges proposed in ref. 19 , where (V p − V q ) and V B account for the Coulomb interactions between charges on different dimers in diagonal and vertical directions, respectively. When V ij /2t d ≳ 0.5-1, the charge order is realized. According to the first principlesbased evaluation 32  single crystal. Figure 3a shows the field dependence of the magnetic torque obtained from a fixed-angle high-field torque measurement at 0.12 and 1.7 K. As shown in Fig. 3a, the magnetic torque shows the field dependence of H 3/2 (the dashed line in Fig.  3a). Given the form of the magnetic torque M × H, the field dependence of H 3/2 shows M / ffiffiffi ffi H p . Note that the free impurity spins are not responsible for this magnetization, since otherwise, the saturation should take place at~1 T for 0.12 K, which is not observed in our data.
Throughout the whole sweep of H, the M-H curve shows neither a remnant field nor a hysteresis (Fig. 2c). The absence of hysteresis is further confirmed down to the lowest fields by our magnetic torque measurements (Fig. 3b). Our data excludes the spin-glass-based weak ferromagnetism picture presented in the previous study 20 , because both a ferromagnetic state 22 and a spin glass state 23 is known to exhibit clear hystereses in the torque  Low-field data of the magnetic torque at 5, 10 and 20 K. The field angle was fixed at 45°from the caxis to the a-axis. Field sweep is given as (1) 0 ! þH max , (2) þH max ! ÀH max , and (3) ÀH max ! 0. measurements from a remnant field and a frozen moment, respectively. In ref. 20 , the authors analyzed the magnetic susceptibility consistent with ours (see Section A in SI) using a narrower fitting range and concluded the antiferromagnetic Curie-Weiss temperature. They speculated that such behavior comes from the spin-glass and the ferromagnetic polarons that may emerge in the presence of charge carriers. However, the studied material has no such carriers, and our torque measurement excludes the spin-glass picture. Their picture on weak ferromagnetism is qualitatively different from ours.
To investigate the magnetic state below T MI in detail, we measured the magnetic torque curves by rotating the magnetic field in the a-c plane (Fig. 4a), where θ denotes the angle between the field and the a-axis (see Fig. 4d). The magnetic torque signal, τ mag ¼ τ 2θ sin 2ðθ À ϕ 2θ Þ, is obtained after subtracting the sin θ component that comes from the gravity of the sample mass (see Section C in SI for details). These magnetic torque curve measurements allow us to detect the magnitude of the magnetic anisotropy, which is proportional to the amplitude of the 2θ component divided by H 2 (Fig. 4b), and the direction of the magnetic principle axis by the phase ϕ 2θ (Fig. 4c). In the metallic T > T MI phase, ϕ 2θ stays at around −20°, which is close to the angle between the long axis of BEDT-TTF molecules and the a-axis (see Fig. 1a), showing that the magnetic anisotropy comes from the spins on the BEDT-TTF dimers [24][25][26] .
At T MI , ϕ 2θ shows a sharp jump which is followed by a rapid shift of ϕ 2θ toward zero, while at the same time τ 2θ stays nearly temperature independent in contrast to the increase of χ M (T). These contrasting temperature dependencies indicate that the magnetic principle axis varies concomitantly with the decrease of the magnetic anisotropy below T MI . Since Raman 16 and IR 14 measurements observed no change of the phonon spectrum below T MI , the change of ϕ 2θ cannot be attributed to the rotation of the BEDT-TTF molecules. Therefore, this ϕ 2θ shift is given by an emergence of a magnetic easy axis parallel to the a axis caused by the ferromagnetic interaction appearing below T MI . A similar but much smaller phase shift has been observed in κ-CN 26 , which may be ascribed to an additional moment from valence bond defects 27 . We further find a characteristic temperature T *~2 4 K. Below T * , ϕ 2θ drops, τ 2θ /H 2 increases, and both τ 2θ /H 2 and ϕ 2θ starts to depend on the field strength. The increase of the magnetic anisotropy particularly developing below T * is consistent with the anisotropy of χ M (T) observed in the previous measurement 20 , supporting the magnetic origin of these temperature changes. This temperature-dependent change is larger for lower fields; as we saw in Fig. 3b the torque data at |H| ≲ 0.5 T changes its sign below 20 K. Another bump-like feature in ϕ 2θ is observed around 7 K, implying a further change of the magnetic state. These features might be related to the changes of the relaxation times observed in NMR measurements done at higher fields 28 . We thus observed a distinct change in the magnetic property already starting below T MI via two torque parameters. Further magnetic torque measurements performed in a dilution refrigerator reveal no change in the magnetic torque below 2 K (Fig. 3a), showing saturation of the temperature dependence.

DISCUSSION
Our magnetic measurements on κ-Hg-Br disclose an unconventional magnetic state, which to our best knowledge has never been observed in the other family members of κ-(BEDT-TTF) 2 Fig. 4 Development of magnetic structures in lowering temperature observed by the magnetic torque measurements. a Magnetic torque curves at different temperatures, normalized by H 2 to compare the ones for different fields (see Section C in SI for details). The data below (above) 40 K was measured at 1 T (3 T). b, c Temperature dependence of the torque amplitude τ 2θ divided by H 2 (b) and the phase shift ϕ 2θ (c) of sin 2θ component of the torque curve. The magnetic susceptibility shown in Fig. 2a is also plotted on the right axis as the dashed line in (b). d Schematic drawing of magnetic torque measurement by a resistive cantilever.
Let us first start by elucidating the way the charges are localized at T < T MI . Most of the previously known κ-(BEDT-TTF) 2 X become a dimer Mott insulator depicted schematically in Fig. 1c. In a Mott phase, the dominant magnetic interactions between the spins carried by the localized charge are always antiferromagnetic as they originate from the kinetic exchange as, J AF ¼ 4t 2 ij =ðU À V ij Þ where t ij and V ij are the transfer integrals and intermolecular Coulomb interaction along the exchange bond, and U is the onmolecular Coulomb interaction. Then, the antiferromagnetic order of κ-Cu-Cl and quantum spin liquid nature of κ-CN are roughly understood by the square-like and triangular lattice geometry of J AF which amounts to 500 29 and 250 K 30 , respectively. These values of J AF are evaluated consistently from both the experiments as well as first-principles calculation where one may consider the actual Coulomb interaction to be 30-40% smaller than the ab initio based ones (see Section D in SI). Therefore, the positive Θ CW observed in κ-Hg-Br cannot be explained by the magnetism of a dimer Mott insulator.
In fact, the abrupt increase of resistivity just below T MI in both κ-Hg-Br and κ-Hg-Cl is different from the crossover behavior usually observed in dimer Mott materials 4 , signaling some sort of translational symmetry breaking of charge distribution. However, the Raman spectroscopy measurements indicate the absence of static charge ordering in κ-Hg-Br 16 . A scenario compatible with all these findings is the dynamical and inhomogeneous charge distribution in between the dimer Mott and charge-ordered state. The intra-dimer transfer integral from the first-principles calculation on κ-Hg-Br is t d~1 20 meV (Valentí, R. & Jeschke, H. private communications and ref. 13 ), much weaker than the typical valuẽ 200 meV of the κ-salts 31 , and thus a quasi-charge-order by the inter-dimer Coulomb interactions is a reasonable expectation.
In Fig. 1d, we locate κ-Hg-Br and κ-Hg-Cl according to the firstprinciples-based evaluations 32 . In the uniform charge-ordered case possibly realized in κ-Hg-Cl, J AF forms long quasi-one-dimensional (1D) chains (see the green lines in Fig.  1c). Here the vertical stripe charge configuration is possibly favored for the Coulomb-interaction-strength of V q < V p , V B of the material 32 .
When the static and bulk charge order is no longer stabilized in κ-Hg-Br, these chains shall break up into short fragments separated by Mott-dimers, as shown in Fig. 5a. The way to construct the domain is not really random; we assume that the chain length N roughly corresponds to the correlation length of charges, and a Mott-dimer is inserted between the chains running in the t q direction, while in reality sometimes there will be connections with the dimers through t p in the other directions. Inside the 1D fragment, the spins interact along the t q -bonds via J AF~1 70-300 K, which will give t d /J AF~4 -8 (see Section D in SI). The charge on a Mott-dimer fluctuates, with fluctuation parameters defined by t d values. During these fluctuations, the charge (and relevant spin) occupies either the left or right molecule on the dimer and interacts with S = 1/2 (green circle) at the adjacent left/right end of the chain via J AF . For the charge configuration shown in Fig. 5a, J AF shown in green and red bonds have the same amplitude.
To elucidate how these quantum fluctuations modify the dominant antiferromagnetism, we construct a synergetic model consisting of two open chains with N L and N R spin-1/2's and a single electron with S = 1/2 (which we call dimer-spin) as shown in Fig. 5d. The Hamiltonian is given as n L=R ¼ c y L=R c L=R , and S d is the dimer spin. This model cuts out the locally interacting manifold of spins shown in Fig. 5a. Such charge configurations behind the model are expected for κ-Hg-Br at temperatures less than T MI , where we find no indications of longrange order of both charges and spins. Experimentally, the charges in an ordered domain show intra-dimer fluctuation altogether via 1.3 THz 15 . In Eq. 1 we have J AF /k B~1 70 K of the same energy scale, which however originates from the inter-dimer charge fluctuation of t q~4 65 K, larger by an order of magnitude. The former is regarded as nearly "static" and are not explicitly seen in our model. The dimers having the t d -term fluctuate alone and do not contribute to the fluctuations of the chains with 1.3 THz frequency. The details of an electronic state below T * are not really known, but the present model does not contradict with the experimental reports given so far.
The model is solved numerically by combining the exact diagonalization calculations (see Section D in SI). Since the total-S z of Eq. 1 is a conserved quantity, we analyze the model by dividing the Hilbert space into different S z -sectors and evaluating the lowest energy levels for each sector. Figure 5c shows 〈S L ⋅ S R 〉 between spins on the left and right chains, S γ = ∑ j∈γ S j for several different series of N γ and system length N. One finds a dominant ferromagnetic correlation (〈S L ⋅ S R 〉 > 0) for large portions of the two lowest excited states. Representative spatial distribution of spin moments for a slightly polarized state is shown in Fig. 5d; the contribution from the constituents of the wave function with dimer-spin on the right and left are separately drawn. The leftupper panel is a typical spin distribution with twofold periodic Friedel oscillation generated by the two open edges of the chain 33 . The dimer-spin hops back and forth, and mixes quantum mechanically with the spins on closer edges of the chains and suppresses their moments. The moments are redistributed throughout the chains and are accumulated densely on the further edges from the center. They point in the same orientation mediated by fluctuating spins closer to the dimer-spin (top panel of Fig. 5d). This interplay of t d and J AF generates robust quantum ferromagnetism (Fig. 5e), which is insensitive to the value of t d and the choice of N γ (see Section E in SI for details).
The theory explains the square-root behavior of the M-H curve at low temperature in κ-Hg-Br. Our calculations show that the ground state of Eq. 1 with even N is always nonmagnetic and has a dominant antiferromagnetic correlation. Let us consider exciting a magnetic moment by applying a magnetic field. Suppose that for an isolated chain with fixed N L and N R , the lowest eigenenergy of Eq. 1 for each S z sector is given as E(S z ). In an applied field H, the system acquires a finite magnetization S z that gives the minimum of energy E H ¼ min Sz ðEðS z Þ À S z HÞ. The "magnetization curve" at finite N L , N R is given as such that H(S z ) = ΔE/ΔS z , where ΔE is the energy difference E(S z ) − E(S z − 1) for ΔS z = 1. As mentioned earlier, the magnetism of the short-range charge-ordered phase shall be described by the assemblage of small magnetic subsystems, interacting with each other, connected with more than two neighboring subsystems. Since the information on the distribution of the chain length is missing, and since the calculation is dealing with only two interacting segments, the direct comparison of the theory and experiments may seem difficult. However, we find that an unbiased comparison is possible as shown in Fig. 5e, where we plotted the magnetization density S z /N as a function of ΔE/N of the subsystems with various different N L and N R . Here, since all the data form a universal square-root curve regardless of the chain length N, it can be interpreted as a stochastic magnetization curve against magnetic field H/J AF . As found in the logarithmic plot, the functional form, ΔE=N $ C ffiffiffiffiffiffiffiffiffi ffi S z =N p , always holds regardless of the length of the chains, while the constant C may depend on the ratio of N L and N R . The universal square-root behavior insensitive to N means that the energy is determined locally. Accordingly, if we consider a bulk assemblage of segments of chains connected by Mott-dimers, their energy shall be an extensive quantity, i.e., the summation of local energy gains. Therefore, we consider this functional form to be intrinsic.
The experimental data is plotted together in Fig. 5e where we add the magnetic torque data (data in Fig. 3a divided by H) into the field dependence of the magnetization data by SQUID (the data in Fig. 2c as it is) so that M estimated from τ 2θ /H coincides to M of the SQUID data at 5 T. The horizontal axis of the experimental data is determined by the value of J AF , and is illustrated for the two parameter choices of J AF /k B = 170 and 300 K discussed in Section D in SI. The shaded region represents the vertical range over which the absolute value of the magnetic moment may vary if the distribution of chain lengths has a large variance, and hence one may state that the theory shows good qualitative agreement with the experiment for any comparable choice of J AF .
An extrapolation of the experimental data in Fig. 5e shows that the magnetization reaches μ B /6 at about 20-30 T. Therefore, approximately, the field strength of 20-30 T which is comparable to T * = 24 K, gives the energy scale to excite the 1/6 magnetic moment from the nonmagnetic ground state. At the temperature range T * < T < T MI , such 1/6 moment is thermally excited and contributes to the ferromagnetism; the ferromagnetically coherent orientation of the moment would contribute to the phase shift of ϕ 2θ . There, M-H curve no longer has a square root, because the low energy magnetic excitations are smeared out. The energy scale of μ B H~0.1J AF to have the μ B /6 moment is consistent with 20-30 T.
Also, the preserved SU(2) symmetry in Eq. 1 matches with the restored isotropy in the magnetic torque at T < T MI . Notice that this ferromagnetic phase is not a long-range order but a correlation because of the one-dimensionality, as can also be suggested from the lack of hysteresis. Below T * the nonmagnetic ground state component becomes dominant. From Raman spectroscopy measurements, the static charge ordering is excluded, whereas the broad peak in ν 2 mode is still compatible with quasi-chargeordered domains with a variant charge disproportionation maximally amounting to ±0.1e, which are coherently fluctuating together inside the domain with a frequency estimated as 1.3 THz 14, 16 .
Once the temperature falls below T * , the dipole (charge) degrees of freedom become correlated over a length scale ξ whose T-dependence saturates and fluctuates slowly together at a corresponding timescale. These fluctuations can safely be integrated out (see ref. 6 ), leading to the effective model of Eq. 1 for spins with antiferromagnetic interactions on chain segments of average length ξ.
While evaluating the precise character of the charge distribution is beyond the scope of any theory currently available, in Section F in SI we provide a phenomenological treatment performed by assuming a functional form for ξ(T) that is valid throughout the critical regime. Within this approach, we show that χ manifests a Curie-Weiss-like behavior that reflects the ferromagnetic correlations between thermally excited spins at temperatures T > T * . Further experimental information concerning the functional form of ξ(T) is required to verify this type of treatment.
One may suspect that the spin models with extrinsic impurities can also explain the phenomena. Although the possibility of spin glass is experimentally excluded, its quantum analog, the random singlet phase may share a similar feature with the present magnetism 34 ; most of the spin moments form a singlet and the remaining spins may contribute to the magnetism. However, for such a state to happen one needs a large amount of static randomness in the distribution of J AF = J(1 ± Δ) that amounts to Δ ≳ 0.6 34 , which cannot happen in the present system.
Ferromagnetism is elusive; For molecular-based materials with only light elements, few ferromagnetic compounds are known, e.g., p-NPNN 35 2 37 , whose exchange interactions originate from the higher-order Goodenough-Kanamori rule. Other mechanisms of bulk ferromagnetism known so far are the Nagaoka ferromagnetism 38 , flat-band 39 , and double exchange or multi-orbitals Hubbard models [40][41][42] , which are applied to metals. The present finding should thus be the first proposal of generating a robust ferromagnetic exchange from the inhomogeneous charge distribution forming dominant spin singlet (paramagnet) by the leading antiferromagnetism interactions. Naively, our ferromagnetism can be viewed as a local doubleexchange; a single charge hops back and forth inside the dimer, and since it interacts with the spin chains on both sides, it is favorable to have both points in the same direction to maximally gain the antiferromagnetic exchange interaction. The difference from the double exchange is that the kinetics of charge is local, keeping the insulating character, and the weakly coupled onedimensional fragments of spin chains do not allow the development of coherent ferromagnetic long-range order.
By designing a three-dimensional critical phase in between the charge order and dimer Mott insulator a two-dimensional ferromagnet might be available, in which case the ferromagnetic long-range order is allowed at a finite temperature.

METHODS Synthesis
Single crystals of κ-(BEDT-TTF) 2 Hg(SCN) 2 Br were grown by electrochemical oxidation of the BEDT-TTF solution 16 . In this synthesis, many single crystals with a variety of the size were obtained. We used a batch of the smaller crystals (powder samples, 3.73 mg) for the SQUID measurements and picked up larger ones for the magnetic torque measurements.

Measurement setups
The magnetic susceptibility of powder samples was measured by SQUID for 2.0-300 K. The magnetic torque measurements were carried out for one single crystal with dimensions 0.75 mm 3 × 0.57 mm 3 × 0.13 mm 3 attached to a piezo-resistive cantilever by a tiny amount of grease (see Supplementary Fig.  3a in SI) by using a variable temperature insert (1.6-200 K) and a dilution refrigerator (0.09-2.5 K). The magnetic torque (τ mag = M × H) is measured by the change of the resistance of the piezo-resistive cantilever. The temperature dependence of the sensitivity of the piezo-resistive cantilever is calibrated by the sin θ component in the torque curve by the gravity of the sample mass as described in Section C in SI. The standard deviations of all the data shown in the figures are smaller than the symbol size.

DATA AVAILABILITY
All the data that support the findings of this study are available from the corresponding authors (M.Y. and N.D.) upon reasonable request.