Charge order affects most of the electronic properties but is believed not to alter the spin arrangement since the magnetic susceptibility remains unchanged. We present electron-spin-resonance experiments on quasi-one-dimensional salts ( PF6, AsF6, and SbF6), which reveal that the magnetic properties are modified below when electronic ferroelectricity sets in. The coupling of anions and organic molecules rotates the g-tensor out of the molecular plane creating magnetically nonequivalent sites on neighboring chains at domain walls. Due to anisotropic Zeeman interaction a novel magnetic interaction mechanism in the charge-ordered state is observed as a doubling of the rotational periodicity of .

1. Introduction

Charge disproportionation is a common ordering mechanism in correlated electron systems that is intensively investigated in transition-metal oxides as well as organic compounds [1, 2] because it is recognized to influence charge, lattice, and magnetic degrees of freedom. For static charge order (CO) the electrons are (partially) localized resulting in an insulating state, most evident in manganites or at the anomaly observed in cuprates [3]. Fluctuating stripes, however, may be relevant for high-temperature superconductivity [4], but also in organic systems evidence accumulates that charge fluctuations mediate superconductivity [57].

As far as the magnetic properties are concerned, the situation is even more puzzling. CO can induce ferroelectricity when the inversion symmetry is broken, as suggested (and also disputed) for some mixed-valence oxides [9, 10], magnetite [1115], manganites [16, 17], or nickelates [18, 19]. For perovskites, however, there seems to be a mutual exclusion of magnetism and ferroelectricity [2022], albeit the field of multiferroics drew enormous attention over the last years [2325]. About ten years ago, electronic ferroelectricity was discovered in quasi-one-dimensional organic compounds (where TMTTF stands for tetramethyltetrathiafulvalene and denotes a monovalent anion such as PF6, AsF6, SbF6, Br) around 50–160 K [26]. X-ray investigations could not identify any superstructure, on this account it is known as a “structureless” transition [27].

In these strongly correlated electron systems, nearest-neighbor Coulomb repulsion causes CO upon cooling below characterized by a charge gap in the transport properties [28]. From the magnetic point of view, these are antiferromagnetic chains with some ordering at low temperatures [29, 30]. Surprisingly, by now no indications of the CO transition could be observed in the magnetic properties, such as spin susceptibility (cf. Figure 1). Here we, report comprehensive electron spin resonance (ESR) investigations that reveal a novel magnetic interaction mechanism caused by charge disproportionation; we observe first evidence that charge order breaks the symmetry of the magnetic degree of freedom in these quantum spin chains.

2. Results and Analysis

Single crystals of with octahedral anions , AsF6, and SbF6, are grown electrochemically. The planar TMTTF molecules stack in a slight zig-zag fashion normal to the molecular plane; this -axis is the best-conducting direction. The stacks are arranged next to each other with some minor interaction in direction and separated in the direction by the anions. Since one electron is transferred to each anion, the organic molecules are left with half a hole, on average. Triggered by the stoichiometry, the TMTTF molecules form dimers along the -axis; a fact that enhances the charge localization observed as a broad minimum in resistivity at K [3135]. In Figures 1(a)–1(c), is plotted for the three compounds under consideration [28].

In addition, basically all TMTTF salts develop charge imbalance for , resulting in charge-rich and charge-poor molecules as first proven by nuclear magnetic resonance [36, 37]. The charge disproportionation can be best visualized by the splitting of the intramolecular vibration plotted in Figures 1(d)–1(f) [3840]. The results are summarized in Table 1; for further details, see [41].

Charge order causes an additional gap to open in the density of states at , leading to a kink in the electrical resistivity indicated by arrows in Figures 1(a)–1(c); as the temperature decreases it evolves in a mean-field fashion to K, for the example of [28].

Although the effect of charge ordering on the spin degree of freedom has attracted considerable attention from the theoretical side [42], it is not completely clear how the spin and orbital magnetic moments behave in the charge-order state and in the metallic state under the influence of charge fluctuations. Naively, one would expect that a phase transition, which has such a severe effect on the electric transport, evidences also in the spin susceptibility, albeit spin-charge separation is an issue for one-dimensional conductors [34, 35]. Despite numerous charge-ordered systems investigated over the years, no sign of the charge-order was revealed in the magnetic behavior. This is demonstrated in Figures 1(a)–1(c) which show for various TMTTF salts [29, 30]. The obtained coupling constants are listed in Table 1.

We have performed comprehensive multifrequency ESR investigations (X-, Q- and W-band) in the temperature range from 4 to 300 K at different orientations of the magnetic field and by rotating around the molecular axes (cf. Figure 6). In Figure 2 the linewidth and the -shift (from the free-electron value ) are plotted as a function of temperature of . The -factor is determined by the electronic wavefunction on the TMTTF molecule, the anisotropy arises from spin-orbit coupling of the conduction electrons [43, 44]; in other words the arrangement of the g-tensor can be solely deduced from the molecular symmetry. At room temperature the -factor and linewidth are maximum when the external magnetic field is oriented along the molecular axis , while and are smallest when measured normal to the TMTTF molecules (-axis which is along the stacks). No indications of the charge-order transition are observed in the -shift; is basically temperature independent except a small reduction of the in-plane anisotropy upon cooling. The intriguing observation is the clear anomaly in the linewidth at which can also be seen in [45]: the linear decrease in becomes weaker indicating an additional interaction channel to develop for . Basically, the same behavior at is observed for , while for the effect is smaller due to stronger coupling and less charge imbalance , which coincide with the much lower transition temperature and almost vanishing charge gap , summarized in Table 1.

In order to gain more insight how charge order affects the spins, we performed detailed angular-dependent measurements of and at selected temperatures. While the rotational symmetry of the remains unchanged with temperature, exhibits a very unusual behavior as demonstrated in Figure 3. For the -factor and the linewidth are dominated by spin-phonon scattering with an angular dependence [29] in the plane, for instance. Most remarkably, the periodicity doubles as the charge-ordered state is entered. In the CO state the orientational dependence of the linewidth can be fitted by the sum of (1) and (2): where is the excess linewidth. The enhancement of the linewidth observed in the diagonal direction becomes more pronounced as the frequency increases, as depicted in Figure 4. The quadratic behavior strongly indicates anisotropic Zeeman interaction [46]. The same period doubling in can be observed in our X-band experiments on , although the effect is much weaker.

In the case of staggered fields, antisymmetric Dzyaloshinskii-Moriya interaction also causes a frequency dependence of the linewidth; however, leading to a temperature dependence [4749] in contrast to our findings. Note also that in those references Heisenberg spin chains with significant intrachain Dzyaloshinskii-Moriya interaction were considered. In the present case, the alternating contributions result from interaction between the chains, for which the iostropic exchange is already small. Thus the antisymmetric Dzyaloshinskii-Moriya exchange can be neglected.

From similar observations on inorganic linear-chain compounds, Pilawa derived a contribution of the anisotropic Zeeman interaction to the linewidth [50] where is the exchange interaction between two inequivalent magnetic sites on the neighboring chains; is the difference between their g-tensors. In our case the interchain exchange constant is K [51]. Note, their Hamiltonian differs by a factor of 2 compared to the one used in [50]. With  kOe for the W-band, we obtain and for and , respectively. In the case of we performed only X-band measurements which in general yield results similar to the Q-band; we estimate . The fact that we do not observe two distinct features but only a single Lorentzian line implies that the hopping rate between the nonequivalent sites is much larger than the splitting (strong-coupling limit). The difference in Larmor frequencies in the W-band GHz corresponds to only 27 mK, which is much smaller than .

The magnetically inequivalent sites have to be located on adjacent stacks; intrachain interaction is so strong that it would suppress the anisotropic Zeeman effect at these frequencies. From the fact that the maximum of occurs at and that is not zero, we conclude that the g-tensors of the two sites do not coincide with the principle magnetic axes any more but are rotated around the -axis by an angle in opposite directions as depicted in Figure 5; at temperature well below we obtain for , and for , respectively Note that this effect is distinct from the gradual rotation of the g-tensor around the -axis due to thermal contraction, reported by Furukawa et al. [44]. That rotation causes a reduction of the g-value anisotropy in the plane as the temperature is reduced (Figure 2).

3. Discussion

The crucial question now is why CO triggers a rotation of the g-tensors. ESR experiments by Furukawa et al. [44] revealed that the Coulomb potential of the anions influences the electronic wavefunctions on the TMTTF molecules. Depending on whether the adjacent molecules are charge rich or charge poor, this constitutes nonequivalent charge and spin distributions as sketched in Figure 5. The imbalance of the spin density causes a rotation of the g-tensor by an angle of . Our findings are in accord with the suggestion by Riera and Poilblanc [52] that charge order in TMTTF salts is a cooperative effect between the Coulomb interaction and the electronic coupling of the TMTTF stacks to the anions The charge disproportionation also causes an imbalance in the Coulomb attraction on the anions, which become slightly distorted in bc plane in accord with thermal expansion results by de Souza et al. [53].

The occurrence of spontaneous dipole moments [26] and second harmonic generation [55] evidences the loss of inversion symmetry (with one center located at the anions, for instance) for . Our vibrational spectra (Figure 1) show the presence of charge-rich and charge-poor molecules, which, however, can be arranged on the dimerized stacks in two possible ways: +− +− +− or −+ −+ −+. This immediately implies that the rotation of the g-tensor can occur in both directions . Thus at domains are formed with inequivalent magnetic sites on neighboring stacks at the boundaries, as depicted in Figure 5. Without interaction two distinct ESR signals are expected below , but the strong coupling leads to a significant broadening, that might get resolved by going to even higher frequencies.

Structural investigations could not identify any superstructure or change in symmetry, however, a strong increase in mosaicity at , as shown in Figure 7 [54]. We conclude that the ordering occurs on a local scale with nanometer-sized CO domains. Similar, but larger ferroelectric domains of electronic origin have been directly observed in the two-dimensional organic compound - below the CO transition [56].

As the averaged g-tensor does not abruptly change at , the charge ordering can be pictured as follows: For random charge disproportionation exists on the TMTTF stacks giving rise to an averaged -factor as observed, but without any anisotropic Zeeman effect, which is suppressed due to narrowing by large intrachain exchange. Below , the charge imbalance becomes cooperative within each stack; there is some weak interaction of adjacent chains. At the boundary of two domains defined by the direction of charge and spin distortion, translational symmetry is lost. This gives rise to two inequivalent chains that yield the same averaged g-tensor along the -Bindestrich and -Bindestrich axes as above . In the diagonal directions, however, the anisotropic Zeeman effect now becomes visible in our ESR data because the signal is narrowed by the much weaker interchain exchange.

4. Conclusions

In summary, the charge order observed in the linear-chain compounds by the splitting of the molecular vibrations opens a gap in the charge excitations. It also affects the magnetic behavior via mutual interaction with the anions. For the inversion symmetry is broken; charge disproportionation leads to a different electronic distortion on the organic molecules that eventually rotates the g-tensors, with two possible directions . From our multifrequency ESR measurements we conclude that in the charge-ordered regime two inequivalent magnetic TMTTF chains coexist, which produce a doubling in the angular periodicity of the linewidth as well as the characteristic quadratic frequency dependence. Thus, our experiments show that charge order not only causes ferroelectricity, but can also break the symmetry of the magnetic degree of freedom in quantum spin chains.


A. Structure

The Fabre salts are charge transfer salts consisting of stacks of the planar organic molecules TMTTF (which stands for tetramethyltetrathiafulvalene) along the -axis that are separated in direction by monovalent anions . The overlap of the orbitals in direction makes them model systems of one-dimensional metals. In -direction the distance of the stacks is comparable to the van der Waals radii. Basically, no coupling between the molecular stacks exists in the direction. In the case of the Bechgaard salts the TMTSF molecule contains selenium instead of sulphur with more extended orbitals leading to better metallic conduction and even superconductivity [32, 5759].

All compounds of the TMTTF family are isostructural (triclinic space group, symmetry with ) with the TMTTF stacked normal to the molecular plane in direction, slightly displaced in direction to optimize the orbital overlap of the selenium atoms along . In spite of the stoichiometry and the corresponding dimerization, the separation of the TMTTF molecules along the stacks only differs by approximately 3% at room temperature and even less at low temperature [28, 60], implying almost equal overlap integrals along the -axis. Due to the triclinic symmetry, commonly denotes the projection of the -axis perpendicular to , and is normal to the plane. Here we introduce a Cartesian coordinate system linked to the TMTTF molecule, with points normal to the molecular plane, is the extended molecular axis, and is perpendicular to both (cf. inset of Figure 2). In a good approximation is parallel to the stacking axis , and is almost , as depicted in Figure 6.

X-ray investigations of single crystals performed at different temperatures down to K rule out a structural transition at the charge-order temperature [54], confirming previous attempts by Coulon and collaborators [27, 60, 61]. Structural data do not unambiguously reveal the loss of inversion symmetry. Their analysis, however, yields a marginal better fit with no symmetry restrictions; but this also holds (to a minor degree) at room temperature. Charge order does not constitute a superstructure; nevertheless, the spots broaden significantly upon cooling with a distinct change at . This is a measure of the mosaicity of the sample: in Figure 7 the temperature dependence of the mosaicity is plotted normalized to the room-temperature value. We interpret this as the development of domains with the anions slightly distorted in one or the other direction. The distortion, but also the domain size changes with temperature leading to a maximum around . Typical for an effect of disorder, the phenomenon depends on cooling rate and slightly vary from crystal to crystal.

B. Experimental Details

Single crystals of with , and were grown by electrochemical methods in an H-type glass cell at room temperature [62]. A constant voltage of 1.5 V was applied between platinum electrodes with an area of approximately . The current through the solution was between 1 and 2 A. After several weeks we were able to harvest needle-shaped single crystals of several millimeters in length and less than a millimeter in width.

The electrical transport was characterized by dc resistivity measured as a function of temperature in all three directions [28]. We evaporated small gold contacts onto the natural crystal surface and attached thin gold wires by carbon paste. Along the long -axis of the crystals and also for the direction four-point measurements could be performed, while for the direction two contacts were applied on opposite sides of the crystal. The samples were attached to a sapphire plate in order to ensure good thermal contact and slowly cooled down to helium temperatures.

Measurements of the magnetic susceptibility were performed using a Quantum Design SQUID (superconducting quantum interference device) magnetometer between 1.8 and 380 K. A large number of single crystals were glued parallel to each other inside a plastic straw using vacuum grease, and the magnetic field was applied parallel to the -axis. The background signal of the sample holder and the vacuum grease was measured separately and subtracted in order to obtain the intrinsic magnetization of the sample. To estimate the spin susceptibility, we subtracted a temperature independent diamagnetic contribution of the core electrons /mole from the original data of all specimen under consideration.

The infrared reflectivity spectra polarized along the direction were recorded by utilizing a Bruker IFS 66v/s and a Bruker Vertex 80 Fourier transform spectrometer extended with a IR microscope HYPERION which is purged with nitrogen to suppress any disturbing influence of the water absorption bands. The reflectivity spectra were measured in a frequency range from until 8000 cm−1 at temperatures from 10 K up to 290 K. Here we focus on the perpendicular direction, that is . To receive the real optical conductivity, we performed a Kramers-Kronig analysis with a constant extrapolation for the low-frequency range and a and decrease for the high-frequency part. With a linear interpolation between the mode frequency in fully neutral (1628 ) and fully ionized (1548 ) TMTTF molecules [63], we can quantitatively determine the charge imbalance [], and the molecular ionicity.

The electron spin resonance (ESR) spectra were carried out in a continuous wave X-band spectrometer Bruker ESP 300 at 9.5 GHz and a W-band spectrometer Bruker Elexsys 680 at 95 GHz at Stuttgart University, and a Q-band spectrometer Bruker Elexsys 500 at 34 GHz at Augsburg University. Complementary X-band measurements were performed at the Hochfeld-Magnetlabor Dresden. The temperature dependence of the ESR properties was measured down to K by utilizing continuous-flow helium cryostats [29, 30]. The typical size of the single crystals used in the X-band ESR measurements was , while very small samples of less than were used in the Q- and W-band ESR measurements in order to avoid any sample-size effect that can lead to a broadening of the linewidth [41].

C. Magnetic Characterization

The large thermal expansion of the organic compounds has strong effects on the temperature dependence of the spin susceptibility. To compare the experimental results (usually obtained at ) with the theoretical predictions (in general calculated for ), the spin susceptibility at constant pressure has to be transformed in the spin susceptibility at constant volume . In the case of the temperature dependence of was estimated by Wzietek et al. [64] from NMR and X-ray measurements under pressure. We assumed that the substitution of sulfur for selenium and the exchange of the inorganic anions has no considerable influence on the thermal expansion, thus we took the ratio for different temperatures to rescale our susceptibility data of . At high temperatures the spin susceptibility at constant volume of resembles the well-known behavior of a spin-1/2 Heisenberg chain with antiferromagnetic coupling as depicted in Figure 1 by the dashed blue lines. The thermodynamic and magnetic properties of such a system were studied by Bonner and Fisher [65]. The advanced model of Eggert, Affleck, and Takahashi [8] using the Bethe ansatz differs appreciably only at low temperatures (). For K the ESR intensity at constant volume can be modeled using , 410, and 400 K for , , and , respectively.

The temperature dependence of the ESR results along the three crystal axes of , , and are shown in Figure 9. The -shift for all three crystal axes is very small. At elevated temperatures, the -values show no significant temperature dependence but a distinct anisotropy with a negative value of along the chain direction and positive values perpendicular to the chain direction, as seen in Figure 8. This behavior is found in all TMTTF and TMTSF salts as presented in [29, 30]. In all cases the small anisotropy of between the - and the -axes increases with temperature. The -value is largest for and closest to the free-electron value parallel to the stacks. As the temperature is lowered, the linewidth decreases almost linearly with in all directions. In the direction is slightly larger; the most narrow lines are always observed along the chains.

At low temperatures and undergo spin-Peierls transitions at and 13 K, respectively. The formation of states with singlet-paired holes due to the tetramerization causes a decreasing susceptibility with an exponential drop of the intensity down to lowest temperatures, but it does not completely vanish (Figure 1). This phase transition also leads to a broadening of the ESR linewidth and a small change in the -values [29]. , on the other hand, orders antiferromagnetically at K, leading to a decrease in with decreasing temperature before it eventually turns up again. This phase transition is accompanied by a decrease in the -shift for the and axes; along the -axis the -factor becomes larger below K. The low-temperature ordered states are out of the scope of the present paper and discussed elsewhere [29, 66].

From the fact that the maximum of occurs at the diagonal direction of the plane and that the -values of neighboring molecules differ by , we conclude that the g-tensors of the spins corresponding to the two different magnetic sites do not coincide with the principle magnetic axes. They are rotated around the -axis by in opposite directions as depicted in Figure 6. As shown in Figure 5, the two different sites contribute to the observed -value according to where and correspond to the g-tensor of type 1 and type 2 spins, respectively. The principal axes of the two kinds of spins span the angle . and are the maximum and the minimum of the -values for each kind of spins at the angle , respectively. For symmetry reasons the -values have to be identical at and , while at the difference of the -values is obtained from the anisotropic Zeeman effect. We obtain for , , and ; for , , and ; and , , and , for the example of .

In Figure 10 the linewidth of . , and is plotted as a function of angle when the magnetic field is rotated around the -axis. While at room temperature (upper frames) the angular dependence follows (1), strongly increases in the diagonal direction () as (lower frames). The effect becomes stronger when going from to and to , corresponding to the increase in charge disproportionation and transition temperature .


The authors thank S. Brown, A. Jánossy, J.-P. Pouget, and E. Rose for useful discussions. The crystals were grown by G. Untereiner, H.J. Kümmerer, H.-A. Krug von Nidda, F. Lissner, and M. Ozerov assisted during some experiments. The work is supported by the Deutsche Forschungsgemeinschaft (DFG). T. Knoblauch acknowledges a fellowship of the Carl-Zeiss-Stiftung.