# Magnetic Ordering in the Frustrated Heisenberg Chain System Cupric Chloride, CuCl

###### Abstract

We report a detailed examination the magnetic structure of anhydrous cupric chloride CuCl carried out by powder neutron diffraction, magnetic susceptibility and specific heat measurements on polycrystalline and single crystal samples as well as an evaluation of the spin exchange interactions by first principles density functional theory (DFT) calculations. Anhydrous CuCl shows one dimensional antiferromagnetic behavior and long range antiferromagnetic ordering below a Nel temperature of 23.9 K. Neutron powder and single crystal diffraction reveal that, below 23.9 K, CuCl undergoes a phase transition into an incommensurate magnetic structure (propagation vector (1,0.2257,0.5) with a spin-spiral propagating along and the moments confined in the crystallographic plane. Our DFT calculations show that the spin-spiral results from competing ferromagnetic nearest neighbor and antiferromagnetic next-nearest neighbor spin-exchange interaction along the spin chains. Implications for possible multiferroic behavior of CuCl are discussed.

###### pacs:

75.40.Cx;75.40.Cx;75.25.+z;75.50.Ee## I Introduction

Cupric chloride CuCl has a crystal structure in which layers of composition CuCl stack along the -direction with van der Waals interactions between them. Each CuCl layer consists of CuCl ribbons made up of edge-sharing CuCl square planes such that the axial positions of each CuCl square are capped by the Cl atoms of its adjacent ribbons with long Cu…Cl distances. Thus, each Cu (d, = 1/2) ion is contained in an axially elongated CuCl octahedron, as expected for a Jahn-Teller active Cu ion. The magnetic susceptibility measurements exhibit a broad maximum, characteristic of a short-range ordering predicted for low-dimensional magnetic systems.DeHaas1931 ; Starr1940 ; Stout1962 ; Barraclough1964 Historically, CuCl was among the very first magnetic systems whose bulk magnetic properties were analyzed in terms of a linear spin chain model. Heat capacity and NQR measurements indicated a Nel temperature 23.9 K for CuCl. To the best of our knowledge the magnetic structure of CuCl has remained unresolved until now.

CuCl is a chemically simple one-dimensional (1D) AFM spin 1/2 quantum chain system. However, by analogy to the CuO ribbon chains found in the magnetic oxides such as LiCuVO and LiCuO, the next-nearest-neighbor (NNN) spin exchange coupling along the chain may be essential in addition to the nearest-neighbor (NN) spin exchange in understanding the magnetic properties of CuCl. For these CuO chains, made up of edge-sharing CuO square planes, the NN spin exchange is ferromagnetic (FM), while the NNN spin exchange is strongly AFM, and the resulting geometric spin frustration gives rise to an incommensurate spiral spin ordering and ferroelectric polarization. Similarly, CuCl might undergo a spiral spin ordering within each Cu chain and a long-range Nel ordering as do LiCuVO and LiCuO. In the present work we examine the magnetic structure of CuCl in detail by powder neutron diffraction, magnetic susceptibility and specific heat measurements as well as by evaluating its spin exchange interactions on the basis of first principles density functional DFT calculations.

## Ii Experimental

Single crystals of CuCl were grown in quartz glass ampoules by the Bridgman technique using commercially available powder (Alfa Aesar, ultra dry, purity 99.995%). Special care was taken in the design of the ampoules to support the elevated vapor pressure of chlorine above the melting point of CuCl. Well-crystallized needle-shape crystals, obtained with the needle axis oriented along the crystallographic axis were found to be twinned in the -plane. Due to the moisture sensitivity of CuCl all handling of the samples was done in an Ar filled glove box. Temperature and field dependent magnetizations (0 T 7 T) were measured using a superconducting quantum interference device magnetometer (MPMS, Quantum Design, 6325 Lusk Boulevard, San Diego, CA.). Magnetizations up to 50 T were measured in pulsed magnetic fields at the Laboratoire National des Champs Magntiques Pulss (Toulouse, France) at temperatures of 1.4 and 4.2 K with the field aligned along the crystal axis. The heat capacities of crystals ( 20 mg) were measured using a commercial Physical Property Measurements System calorimeter (Quantum Design, 6325 Lusk Boulevard, San Diego, CA.) employing the relaxation method in magnetic fields up to 9 T aligned parallel and perpendicular to the crystal axis. To thermally anchor the crystals to the heat capacity platform, a minute amount of Apiezon N grease was used. The heat capacity of the platform and the grease was determined in a separate run and subtracted from the total heat capacity. The heat capacity of powder samples ( 300 mg) was measured in a home-built fully automatic Nernst-type adiabatic calorimeter similar to that described in detail in Refs. MBThesis, ; Gmelin, . The powder samples were encapsulated in Duran glass flasks under 900 mbar He atmosphere to enable rapid thermal equilibration. The flasks were attached with a minute amount of Apiezon N vacuum grease to the sapphire platform which carries an deposited thin-film stainless steel heater and a calibrated Cernox CX-1050-SD thermometer (Lake Shore Cryotronics, Inc. 575 McCorkle Blvd. Westerville OH 43082). The heat capacities of the sapphire sample platform, the glass flask and the Apiezon vacuum grease were determined in separate runs and subtracted from the total heat capacity. EPR spectra were measured with a Bruker ER040XK microwave X-band spectrometer and a Bruker BE25 magnet equipped with a BH15 field controller calibrated against DPPH. A stepper-motor-controlled goniometer was used to rotate the crystals around an axis perpendicular to the magnetic field plane. The crystals were oriented optically with a precision of . 100 kHz field modulation was used to record the first derivative of the signal intensity. The spectra were fitted to a resonance curve with Lorentzian lineshape (zero dispersion) by varying the resonance position as well as the linewidth, intensity and background parameters.

Powder neutron diffraction patterns were collected on the high-intensity two-axis diffractometer D20 (Institut Laue-Langevin, Grenoble) at a wavelength of 2.4 Å and 1.88 Å in the temperature range 1.8 K 50K. D20 Single crystal diffraction was carried out on the four-circle diffractometer D10 (Institut Laue-Langevin, Grenoble) at a wavelength of 2.36 Å on a crystal of 50 mm D10

## Iii Experimental Results

### iii.1 Crystal Structure

At room temperature, anhydrous cupric chloride CuCl was reported to crystallize with a monoclinic structure (space group 2/) with 2 formula units per unit cell.Wells1947 ; Burns1993 The crystal structure of CuCl consists of CuCl slabs parallel to the -plane that are made up of CuCl chains running along the crystallographic axis. Such slabs are interconnected along the axis via van der Waals contacts (cf. Figure 1).

To confirm the crystal structure details of CuCl and to search for structural phase transitions at low temperatures, we collected a series of neutron diffraction patterns on the diffractometer D20 at temperatures 2 K 50K. Patterns with sufficient intensity but improved resolution were collected using neutrons with a wavelength of = 1.88 Å (Ge monochromator) suitable for Rietveld profile refinements. These were done with the FullProf program packageFullProf . The patterns were refined assuming the space group 2/ with the Cu atoms in the Wyckoff position (2a) (0,0,0) and the Cl atoms in the position (4i) (,0,) using a slightly asymmetric pseudo-Voigt profile shape to model the Bragg reflections. Within experimental error the atom coordinates and for the Cl atoms are independent of the temperature and amount to =0.5065(2) and =0.2365(3). Figure 2 displays the experimental pattern collected at 42.9 K and the results of the Rietveld refinement (=1.88 Å ). Apart from changes in the lattice parameters, the structure remains unchanged with respect to the room temperature crystal structure of CuCl described by Wells and Burns et al..Wells1947 ; Burns1993 The atom position parameters show no noticeable difference to those of the room temperature structure. The structure parameters at 2.5 and 42.9 K are summarized in Table 1.

(K) | (Å ) | (Å ) | (Å ) | () | method |
---|---|---|---|---|---|

2.55 | 6.7986(3) | 3.29418(12) | 6.7718(3) | 122.45(2) | neut |

42.9 | 6.8010(3) | 3.29437(12) | 6.7768(3) | 122.47(2) | neut |

293 | 6.8973(4) | 3.2961(2) | 6.8160(4) | 122.239(3) | x-ray |

293 | 6.9038(9) | 3.2995(4) | 6.824(1) | 122.197(8) | x-ray (Ref. Burns1993, ) |

The temperature dependence of the lattice parameters is displayed in Figure 3. The lattice parameters show a monotonic decrease with temperature. Below 15 K they level off and remain constant to the lowest temperatures. Near 24 K, especially and , exhibit noticeable changes in the slope, which indicate a magnetoelastic response of the lattice to the onset of long-range AFM ordering, as is evidenced by our magnetization and heat capacity experiments discussed in the following paragraphs.

Single-crystal x-ray diffraction on crystals of CuCl were carried out using a laboratory based diffractometer (PDS II, Stoe & Cie., Darmstadt, Germany). These measurements confirmed the crystal structure determined from neutron powder diffraction and showed the crystals to be twinned in the ac- plane. The twinning matrix T relating the Miller indices of the twinned individua was determined to

(1) |

The twinning reverses the -axis and rotates the the axes in the reciprocal space, such that twin-related axes enclose an angle of 180 - 57.

## Iv Magnetic Properties

### iv.1 Magnetic Susceptibility Measurements

Figure 4 displays the magnetic susceptibility of a polycrystalline sample of CuCl. The susceptibility is characterized by a broad maximum centered at 70 K and a Curie-Weiss law (Eq. (5) behavior at high temperatures (see the detailed discussion of the susceptibility measurement on a crystal below) with a negative Curie-Weiss temperature of -60 K and a powder -factor (see Eq. (5) of 2.13.

The negative Curie-Weiss temperature proves predominant AFM exchange interactions, and the broad maximum reveals AFM short range ordering indicative of a low-dimensional magnetic behavior. A kink-like change of the slope near 24 K is due to the onset of long-range AFM ordering. A -type anomaly is clearly revealed by the temperature derivative of the quantity (’Fisher’s heat capacity’, cf. Ref. Fisher1962, ) shown in the inset of Figure 4 from which the Néel temperature is determined to be 23.9(3) K. coincides with the anomalies observed in the temperature-dependence of the lattice parameters (Figure 3). When fitting a Curie-Weiss law to the high temperature data of the powder susceptibility it is found that the fitted parameters (the slope of the reciprocal susceptibility i.e. the effective magnetic moment, and the intercept) are strongly correlated and depend essentially on temperature-independent contributions to the susceptibility (diamagnetic and Van Vleck susceptibilities).

In order to obtain a reliable fit of the magnetic susceptibility of CuCl over a wide temperature range and as well as the temperature independent diamagnetic and paramagnetic corrections to the susceptibility we determined the magnetic susceptibility on a crystal (=1T ) up to 600 K. At high temperatures the magnetic susceptibility follows a Curie-Weiss law as demonstrated by the plot of the reciprocal susceptibility in Figure 5. Deviations from the Curie-Weiss law are seen already at temperature below 175 K.

At high temperatures, besides the spin susceptibility, , temperature independent diamagnetic orbital contributions, , from the core electronic shells and the Van Vleck paramagnetism, , from the open shells of the Cu ions become important. can be estimated from the increments using Selwood’s tabulated values as -1210 cm/mole and -2610 cm/mole for Cu and Cl, respectively, giving a total temperature independent diamagnetic susceptibility for CuCl of -6410 cm/mole.Selwood1956 An estimate of the Van Vleck contribution, can be obtained from (cf. Ref Lueken1999, ).

(2) |

where is Avogadro’s constant and = 10 is the energy separation of the orbital states in a cubic crystal field, which amounts to approximately 2510 cm. With / 0.375 cm K/mole one estimates a Van Vleck susceptibility of Cu in a cubic crystal field of the order of +4310 cm/mole. In CuCl due to the Jahn-Teller distortion the symmetry is reduced and the Van Vleck contribution depends on the actual matrix elements and the energy separation of the orbital states to the state. For a detailed discussion see for example Ref. Takigawa1989, . Since these energy splittings are not known at present, we used an alternative method to determine the temperature independent contributions : In order to perform this analysis, we assume that the magnetic susceptibility of CuCl can be described by a sum of the temperature dependent spin susceptibility and a temperature independent contribution

(3) |

where is composed by the sum of the diamagnetic and the Van Vleck contribution according to

(4) |

At sufficiently high temperatures, , the spin susceptibility (T) can be approximated by a Curie-Weiss law according to

(5) |

When fitting the high temperature susceptibility to the sum of equations (4) and (5) it is found that the fitting parameters, , and , are strongly correlated, and meaningful results cannot be obtained without fixing one of them. In order to reduce the number of parameters we therefore determined the g-factor, independently by an EPR experiment on a crystal (see below) to = 2.050(1). The EPR experiment also revealed that our CuCl crystals consist of mirror twins in the crystallographic plane with the mirror twins having the same axis. By using = 2.05 and fitting Eq. (3) to the (uncorrected) susceptibility data for (see Figure 5), we determine the following parameters:

### iv.2 Magnetization Measurements

Laboratory based magnetization experiments at 1.8 K with evidenced a spin-flop transition near 4 T. Measurements in a pulsed magnetic field at 1.4 K up to 50 T confirmed this result. Evidence for further magnetic phase transitions at fields larger than 4 T could could not be found. Figure 6 displays the magnetizations measured in the pulsed field experiment. Saturation is not seen up to 50 T. The 1.4 K magnetization trace shows a slight upward deviation from nonlinearity while the 4.2 K traces within experimental error show a linear dependence of the magnetization over the full field range.

Magnetic susceptibility experiments in various magnetic fields reveal that depends only weakly on the field. The combination of the magnetization and the magnetic susceptibility experiments allows us to construct a magnetic phase diagram for as displayed in Figure 7 with a bicritical point at 5.4(1) T and 24.2(1) K where the paramagnetic, spin flop and Néel phases coexist.

### iv.3 Electron Paramagnetic Resonance Measurements

In order to determine the -factors along the principal axes of the -tensor we carried out EPR measurements on small crystals of CuCl at room temperature at a microwave frequency of 9.46 GHz (X-band) (resonance field 0.33 T). In a monoclinic system, one principal axis of the -tensor is fixed by symmetry and oriented parallel to the twofold axis, i.e. the crystallographic axis for CuCl. The remaining two principal axes of the -tensor lie within the plane. To obtain the -factor along , a crystal was oriented with the -axis within the magnetic field plane. Rotation was carried out around an arbitrary axis perpendicular to . A single resonance line is observed. Figure 8 shows the angular dependence of the -factor obtained from fits of a field derivative of a Lorentzian to the resonance line. The angular dependence of the factor shows the typical oscillating behavior. The minima occur when the field is oriented parallel to the axis. The -factor with field oriented parallel to the -axis amounts to

The linewidths of the resonance lines exhibit angular dependence and they vary between 20 mT and 40 mT due to small anisotropic exchange components.MBThesis Details will be discussed in a forthcoming paper. Rekretbp

If the magnetic field was oriented in the plane of the crystal (rotation around -axis), we observed an EPR spectrum that consists of two overlapping resonance lines, indicating that the crystals were twinned. The resonance fields and the linewidths of each line and the intensity were obtained by fitting to the spectra the sum of two Lorentzian absorption lines (see Figure 9).

The -factors of the two lines are shown in Figure 10. One observes the typical oscillating angular dependence of the -factors on the angle between magnetic field the crystal orientation. Both lines show the same oscillating angular dependence of the -factor, however with respect to each other shifted by 57. The extremal values are identical within experimental errors. The intensities for the two crystallites are about the same indicating that the investigated twinned crystallites have about the same volume. At angles where both twin-individua have the same resonance field and -factors, it is difficult to deconvolute the two resonance line unambiguously. This leads to the rather high error level in the fitted parameters especially at angles where the two -factor curves intersect.

The principal -factors, and were determined by comparing the results with the general angular dependence of the -factor for an axial system (Ref. Abragam1970, )

(6) |

where is the angle between the external field and the direction of the principal axis . is a phase factor adjusted individually for each individual crystal.

With = 2.050(1), these values result in an average -factor of

(7) |

This valueis in good agreement with the -factor obtained from the powder susceptibility measurement (see. Figure 4).

and are almost identical indicating that the local symmetry for Cu in CuCl is very close to axial with =2.255 and 2.05. Here, we have identified the large with the -factor along the quasi-fourfold symmetry axis, perpendicular to the CuCl squares in the ribbons and and with the -factors for the magnetic field aligned within the planes of the ribbon.

This finding reflects the local symmetry, which is very close to tetragonal given the elongated CuCl octahedron containing each Cu ion. The four Cl atoms in the equatorial plane are all at the same distance of 2.253 Å and the apical Cl atoms at a distance of 2.911 Å ( = 42.9 K); the normal of the equatorial plane and the direction to the apical Cl atoms encompass an angle of only 3.2.

The deviation of the -factors of the ground state of a 3 configuration in octahedral field with tetragonal splitting from the -factor of the free electron is determined by the spin-orbit coupling parameter = -/2, where is the one-electron spin-orbit coupling parameter and =1/2. For Cu amounts to 829 cm (Ref. Abragam1970, ) and the deviations from the free electron -factor are given by

where and are the energy separations between the state and the and the electronic states, respectively. Abragam1970 Using our experimental results, and amount to 2610 cm and 3010 cm, respectively, consistent with the estimate we made from the magnitude of the temperature independent Van Vleck contribution.

The principal axes of the two crystallites are shifted by 180 - 57 with respect to each other, where is the monoclinic angle (cf. Table 1). This finding indicates that CuCl crystals generally are twinned with the twins related to each other by a mirror symmetry operation with the mirror symmetry plane in the plane making the crystals to appear pseudo-hexagonal. This finding becomes essential for the assignment of the propagation vector of the magnetic structure from elastic neutron diffraction on CuCl crystals.

## V Thermal Properties

### v.1 Heat Capacity Measurements

The heat capacity of anhydrous CuCl has been determined before by Stout and Chisholm (Ref. Stout1962, ) and by Billerey et al. (Ref. Billerey1982, ). Figure 11 shows our heat capacities of a CuCl crystal and a powder sample. Our data agree very well with the data published by Stout et al. and those of Billerey and collaborators. The data of Stout and Chisholm are shown in Figure 11 for comparison. The heat capacities of our powder sample and the crystal are very similar, and they exhibit -shaped anomalies due to the onset of long-range AFM ordering. The anomalies are centered at 23.8 K (crystal: 23.87 K; powder 23.82 K) in good agreement with the values reported by Stout et al. (23.910.1 K) and Billerey et al. (23.90.1 K). The maximum value in the -anomaly is somewhat larger than the value reported by Stout et al. and Billerey et al. possibly indicating that our samples structurally are better ordered and contain less structural defects. The -anomaly does not change its shape nor shifts the position when a magnetic field of up to 9 T is applied perpendicular to the crystal axis. Sufficiently below the anomaly where the lattice contribution can be approximated by a Debye- power law, the total heat capacities follow well the same power law = + indicating magnon contributions from a three-dimensional antiferromagnet.

The entropy removed by long-range ordering signaled by the anomaly at 23.8 K amounts to 0.60(5) J/molK equivalent to 10% of ln2 as expected for a = spin system. The remaining 90% of the entropy apparently are removed by short-range ordering above .

The total heat capacity, , of CuCl contains the spin (magnetic) and contributions from the phonon system (lattice contribution) according to

(9) |

In order to separate the magnetic contributions, , from the total heat capacity of CuCl, one has to determine the lattice heat contributions, . The lattice contributions can be estimated by fitting a power series expansion to the high temperature regime, where magnetic contributions are negligible, and extending the series to low temperatures. A significantly improved estimate of the lattice heat capacity is usually achieved by measuring the heat capacities of isostructural nonmagnetic compounds. The lattice reference systems appropriate for CuCl would be MgCl and CdCl. Both compounds crystallize with a layered crystal structure similar to that of CuCl.Partin1991 However, the MCl (M=Mg, Cd) octahedra of MgCl and CdCl do not exhibit a Jahn-Teller distortion, and consequently the two-valent anions are coordinated regularly by six Cl atoms at the same distance of 2.64 Å and 2.96 Å for CdCl and MgCl, respectively.Partin1991 .

Billerey et al. determined the heat capacities of MgCl in the temperature range 2 K 100 K and constructed the lattice contribution to the heat capacity of CuCl. We re-measured the heat capacities of a polycrystalline sample and of crystalline pieces of high purity (Alfa Aesar 99.99%) MgCl in an extended temperature range. We found that our data for MgCl obtained on two independent samples with two different experimental setups, especially, in the temperature regime relevant for the magnetic part below 100 K deviate markedly from Billerey’s data. Our values being up to 10% smaller then Billerey’s results.

Additionally we determined the heat capacity of a small polycrystalline beads, typical diameter 2mm) of CdCl (Alfa Aesar 99.998%) The heat capacity of CdCl has been measured before by Itskevich et al.(Ref. Itskevich1956, ) in a very limited temperature range. Takeda et al. (Ref. Takeda1984, ) also determined the heat capacity of CdCl in order to subtract the lattice contribution to the heat capacity of magnetic VCl, but no data are available from the work of Takeda et al.. Only the low-temperature cubic term is quoted. Figure 12 shows the heat capacities of MgCl and CdCl as determined in this work in comparison with literature data.

Figure 12 displays the heat capacities of CuCl and of MgCl and CdCl for comparison. As expected, the heat capacity of MgCl is significantly smaller, that of CdCl markedly greater than that of CuCl. A detailed discussion how to construct an adequate lattice reference for CuCl and the magnetic contributions to the heat capacity of CuCl will be given in Section VIII.

## Vi Neutron Diffraction

### vi.1 Elastic Magnetic Scattering of a Powder Sample and a Twinned Crystal

In order to search for magnetic Bragg reflections of the long-range ordered AFM state we examined in detail the low-angular parts of the powder diffraction patterns. A comparison of two diffraction patterns collected in D20’s highest intensity mode at a wavelength of =2.4 Å at 2 K and 30K reveals additional magnetic scattering up to 2 40 (Figure 13(a)). Most pronounced is a reflection 23 at the low angle shoulder of the nuclear 0 0 1 reflection. Additional, though significantly less intense, magnetic scattering intensity is visible in the 2 K pattern at 34. Both features are not present any more in the pattern collected at 30K and thus can be clearly identified as magnetic Bragg reflections.

Towards low scattering angles, the 30K pattern lies clearly above the 2 K pattern indicating significant paramagnetic scattering in the 30 K pattern. In the difference pattern the magnetic Bragg reflection near 19 is slightly deformed owing to a small shift of the 0 0 1 nuclear reflection arising from the contraction of the lattice between 30 K and 2 K.

In order to obtain the temperature dependence of this first magnetic Bragg reflection, we collected a series of patterns between 2 K and 45 K in the D20’s higher resolution mode using a wavelength of =1.888 ÅṪhe difference between the pattern collected at the respective temperatures and the pattern at 28.9 K was fitted to a Gaussian the intensity of which is displayed versus temperature in Figure 14.

The temperature dependence of the integrated intensities can be well fitted by a critical power law according to Eq. (10) with the critical exponent

(10) |

A fit of the integrated intensities converges to a Nel temperature =24.7(1) K consistent with the magnetic susceptibility and the heat capacity and a critical exponent =0.16(1) (see Figure 14(a)).

Intensities of a set of magnetic Bragg reflections were collected on a twinned crystal. Figure 14(b,c) display the intensity of a strong magnetic reflection collected on a twinned crystal on ILL’s four-circle diffractometer D10.MBThesis ; Rekretbp Initially, this reflection was indexed as 0 0.776 0.5. In view of the magnetic structure refinement of the powder data (see below) we must conclude that this reflection originates from the related twin individuum in the crystal. Using the inverse of the twin matrix eq. (1) the indexes transform to 0 -0.776 -0.5 i.e. 1 1 0 - and 1 1 1 - , where is the magnetic propagation vector refined from the powder data below. A fit of the integrated intensities of this reflection (solid line in Figure 14(b,c)) gives =23.9(1) and =0.23(1) consistent with the powder diffraction results.

Although the investigated temperature regime is limited, and corrections for critical scattering near have not been applied, the refined critical exponents obtained from the powder and crystal are significantly lower than those characterizing the well known three-dimensional universality classes (Ising, XY and Heisenberg) for which ranges between 0.325 and 0.365.Guillou1977 The refined critical exponent 0.2 for CuCl, however, comes close to the predictions of Kawamura for the universality class of of the SO(3) or the Z S symmetry for systems with chiral order. For such systems ranges between 0.22 and 0.25.Kawamura ; Collins1989 ; Gaulin1994 Figure 14(c) displays the intensity of a magnetic reflection collected on a twinned crystal on ILL’s four-circle diffractometer D10 together with a power law with a critical exponent = 0.24(2)

The difference between the 2 K and the 30 K powder diffraction patterns (lower panel of Figure 13) reveals an additional weak triple of magnetic Bragg reflections between 32 and 37 with a peak shape very similar to those of neighboring nuclear reflections. Additionally, the difference patterns allows one to identify magnetic scattering around a scattering angle of 70.

All magnetic Bragg reflections can be indexed with a propagation vector (1, 0.225,0.5) indicating a doubling of the magnetic cell along the c axis and an incommensurate AFM ordering along b (see Figure 15).

The magnetic structure of CuCl was refined using the D20 2 K and the difference powder diffraction data sets using the program Fullprof. In order to obtain reliable parameters for the nuclear structure (nuclear scale factor, profile parameters), we first refined the 2 K pattern assuming nuclear scattering only. Since magnetic scattering is very small this refinement converges well and provides the lattice and profile parameters and the scale factor. Subsequently a difference pattern obtained by subtracting the 30 K pattern from the 2 K pattern was used to refine the magnetic structure on the basis of the propagation vector and the nuclear scale factor, which was kept fixed in the refinements. The background of the difference pattern was chosen manually and subtracted. The refinement based on 18 magnetic reflections was successful assuming the nuclear space group and a helix with the moments confined to the plane. Several tests were carried with the moment directions pointing in a more general direction out of the plane but did not lead to significantly improved fits.

Assuming this model for the magnetic structure the ordered moment at 2 K converged to

The propagation vector at 2 K was refined to

implying an angle of 81.2(3) between neighboring moments along the chain. Figure 16 displays the magnetic structure of CuCl at 2 K.

## Vii Density functional characterization of the spin exchange interactions and the spin-spiral state

To discuss the magnetic structure of CuCl, we consider the intra-chain and the inter-chain spin exchange interactions defined in Figure 17. and are the NN and NNN intra-chain spin interactions, respectively. and are the inter-chain spin exchange interactions along the -direction, and is that along the a-direction. The structural parameters associated with the spin exchange paths - are summarized in Table 2.

path | Cu…Cu (Å ) | Cl…Cl (Å ) | Cu-Cl…Cl ( |
---|---|---|---|

3.299 | - | - | |

6.599 | 3.299 (2) | 136.8 | |

6.824 | 3.729 (2) | 133.1 | |

9.493 | 4.979 | 174.1 | |

7.608 | 3.808 | 131.7 |

To estimate the spin exchange parameters - , we first determine the total energies of the six ordered spin states of CuCl presented in Figure 18 on the basis of first principles density functional theory (DFT) electronic band structure calculations. Our DFT calculations employed the Vienna ab initio simulation package (Ref. Kresse1993, ; Kresse1996a, ; Kresse1996b, ) with the projected augmented-wave method, the generalized gradient approximation (GGA) for the exchange and correlation functional, (Ref. Perdew1996, ) the plane-wave cut-off energy of 330 eV, and the sampling of the irreducible Brillouin zone with 36 -points. To take into consideration of the strong electron correlation associated the Cu 3d state, we performed GGA plus onsite repulsion (GGA+U) calculations (Ref. Dudarev1998, ) with = 3, 5 and 7 eV on Cu. The relative energies of the six ordered spin states obtained from our GGA+U calculations are summarized in Table 3.

state | = 3eV | = 5eV | = 7 eV |
---|---|---|---|

FM | 0 | 0 | 0 |

AFM1 | -11.23 | -8.62 | -6.52 |

AFM2 | 21.68 | 20.06 | 17.43 |

AFM3 | 13.77 | 14.09 | 12.95 |

AFM4 | -31.01 | -22.12 | -15.93 |

AFM5 | -39.72 | -28.85 | -21.05 |

The energies of the six ordered spin states can also be written in terms of the spin Hamiltonian,

(11) |

(1) where (= - ) is the spin exchange parameter for the spin exchange interaction between the spin sites and , while and are the spin angular momentum operators at the spin sites and , respectively. By applying the energy expressions obtained for spin dimers with N unpaired spins per spin site (in the present case, N = 1), (Ref. Dai2001, ; Dai2003, ) the total spin exchange energies per formula unit of the six ordered spin states (see Figure 18) are written as

state | = 3eV | = 5eV | = 7 eV |
---|---|---|---|

23.3 | 21.4 | 18.4 | |

-41.4 | -31.9 | -24.5 | |

-9.1 | -7.0 | -5.3 | |

-0.2 | -0.1 | -0.1 | |

-0.8 | -0.7 | -0.5 |

By mapping the relative energies of the six ordered spin states given in terms of the spin exchange parameters to the corresponding relative energies obtained from the GGA+U calculations, we obtain the values of the spin exchange parameters - summarized in Table 4.

For all values of employed, the two strongest spin exchange interactions are the intra-chain interactions and . As found for the magnetic oxides with CuO ribbon chains, (Ref. Xiang2007a, ; Xiang2007b, ) is ferromagnetic while is antiferromagnetic with larger in magnitude than . As a consequence, the intra-chain spin exchange interactions are geometrically frustrated. Except for , the inter-chain spin exchange interactions are negligible; is antiferromagnetic and is weaker than by a factor of approximately four. Because of the intra-chain spin frustration, the magnetic ground state of each CuCl chain is expected to be a spin-spiral state with a certain incommensurate repeat vector along the chain direction, i.e., the -direction. The inter-chain spin arrangement is expected to be antiferromagnetic along the -direction due to the spin exchange . Our calculations of the energy E (0, , 1/2) as a function of the modulation wave vector = (0, , 1/2) shows a minimum at 0.22.

The magnetic state with the spin-spiral arrangement removes inversion symmetry and hence should induce ferroelectric polarization. To confirm this prediction, we carried out GGA+U plus spin-orbit coupling (SOC) (Ref. Kunes2001, ) calculations for CuCl that has the commensurate spin-spiral arrangement with = (0, 0.25, 0) to reduce the computational task. Indeed, our calculations of the ferroelectric polarization using the Berry phase method (Ref. King1993, ; Resta1994, ) for the spin-spiral state of CuCl determined by the GGA+U+SOC calculations lead to nonzero ferroelectric polarizations; with the spin-spiral plane parallel () and perpendicular ( ) to the CuCl ribbon plane, the -, - and - components of the polarizations are calculated to be = (-44.5, 0, 71.6) and = ( 5.8, 0, 28.3) in units of C/m.

## Viii Experiment versus Theoretical Results

In the following we will compare the implications of the theoretical results for the magnetic susceptibility and the thermal properties. The calculations indicate that the intra-chain exchange is dominant and that the magnetic properties of CuCl have to be described as those of a =1/2 frustrated Heisenberg chain with FM nearest and AFM NN interaction. According to the calculations, the ratio / is negative and is expected to be of the order of -0.5 to -0.75. Additionally, there is small but noticeable inter-chain interaction being dominant along the axis that leads to cooperative AFM ordering and the AFM arrangement along the direction with the doubling of the magnetic unit cell along . Magnetic and thermal properties of =1/2 zig-zag Heisenberg chains have been calculated in great detail by exact diagonalization by Heidrich-Meisner and coworkers, and numerical tables are available for a wide range of the ratio /. Heidrich2006 ; Heidrichweb These works comprise also the magnetic susceptibility of a =1/2 Heisenberg chain with NN AFM exchange only, i.e., =0 as calculated before by Klümper et al. and parametrized by a Pad approximant by Johnston et al..Kluemper2000 ; Johnston2000

### viii.1 Magnetic Susceptibility

Correcting the magnetic susceptibility results with the sum of the temperature independent diamagnetic and Van Vleck corrections to the susceptibility () gives the pure spin susceptibility of CuCl. Using the -factor along as determined by the EPR measurements, =2.050(1), enables us to compare the magnetic susceptibility of a CuCl crystal measured with magnetic field along the -axis with model calculations. The only remaining adjustable parameters are the intra-chain exchange constants and . The effect of the finite inter-chain exchange on the paramagnetic susceptibilities calculated for the zig-zag chains has been taken into account by employing a mean field approach, which gives as inter-chain interaction corrected susceptibility as(Ref. Carlin1986, )

(13) |

The leading intra-chain exchange is provided by (see Table 3 and Figure 14) and amounts to about 20% of the next-nearest exchange along the chain. Including intra-chain exchange according to Eq. 13 renormalizes the calculated susceptibilities by about 8% near the maximum at 73 K. At room temperature the decrease amounts to 4%.

Figure 19 shows the comparison of the spin susceptibility of CuCl as determined in this work with the inter-chain exchange corrected susceptibility of a Heisenberg chain with NN and NNN exchange interaction and , respectively.Heidrich2006 ; Heidrichweb

According to this comparison the magnetic susceptibility of CuCl is best reproduced by an inter-chain exchange corrected susceptibility of a Heisenberg chain with NNN exchange interaction =-117 K and a ratio -0.25, i.e., ferromagnetic nearest neighbor exchange as predicted by the DFT calculations. However, the ratio is markedly smaller than predicted by the DFT calculations (see Table 3). Including inter-chain exchange interaction is essential to reproduce the measured susceptibility correctly.

### viii.2 Magnetic Heat Capacity

The magnetic contributions to the heat capacity of CuCl has been evaluated before by Stout and Chisholm and subsequently by Billerey and coworkers. Stout1962 ; Billerey1982 Our measurements confirm (see Figure 11(a)) their peak-anomaly at 23 K indicating long-range AFM ordering. Our data agree very well with those of the previous works, and the anomaly is somewhat sharper in our experiment possibly due to the fact that polycrystalline samples have been used before.Stout1962 ; Billerey1982 Stout and Chishom and Billerey et al. used the heat capacities of MnCl and MgCl, respectively, in order to construct an estimate of the lattice heat capacity of CuCl. We redetermined the heat capacities of MgCl and found marked deviations of our data from those of Billerey. We considered also CdCl as another possible system to estimate the phonon contributions to the heat capacity of CuCl and measured the heat capacity of CdCl over a temperature range much wider than covered by the experiments of Itskevich et al.(see Figure 12(a)). MgCl and CdCl are diamagnetic and have a crystal structure containing MgCl and CdCl layers, which closely resemble those in the crystal structure of CuCl. However, the MCl (M=Mg, Cd) octahedra are not axially elongated as they are in the CuO octahedra. In addition, owing to the different atomic masses of Cu and Mg, the phonon spectrum of CuCl and MgCl and Cd may be expected to be different, especially at low energies. Although MnCl is magnetic, its AFM ordering occurs at very low temperature, 1.81 K and 1.96 K, and magnetic contributions to its heat capacity are negligible above 10K.Murray1955 ; Chisholm1962 Above 10K, MnCl may therefore provide a reasonable approximation to the lattice heat capacity of CuCl. In order to compensate for the difference in the cation masses, we employed the method of corresponding states and assume that the phonon spectra of MgCl, CdCl and MnCl are very similar to that of CuCl and the different phonon energies can be compensated for by modifying the temperatures.Boo1977 Compensation of the different cation mass was done by stretching or compressing the temperature scale by a constant factor. This factor was chosen such that at sufficiently high temperatures where magnetic contributions to the heat capacity vanish the heat capacities of the lattice references match with those of CuCl (see Figure 11). We found factors of 0.88, 0.93 and 1.18 to be appropriate for MgCl, MnCl and CdCl, respectively.

A comparison of the heat capacities of CuCl and MCl (M=Mg, Mn, Cd) reveals significant short-range ordering magnetic contributions to the heat capacity of CuCl up to 100 K, as already concluded by Stout et al. and Billerey Stout1962 ; Billerey1982 .

In order to construct a lattice heat capacity reference for CuCl from the heat capacities of MgCl, MnCl and CdCl, we tried various combinations of the heat capacities of MgCl, MnCl and CdCl with temperatures scaled by the factors given above and subtracted these from the total heat capacity of CuCl. The resulting magnetic heat capacity of CuCl, / was subsequently integrated and the total magnetic entropy was compared with ln2, the entropy expected for a =1/2 magnetic system as a crosscheck.

Due to the very large deviations of the heat capacity of MgCl from that of CuCl, especially in the temperature regime below 50K where magnetic contributions are essential we discarded MgCl and achieved an adequate lattice heat capacity reference by averaging for each temperature the heat capacities of MnCl and CdCl.

The magnetic heat capacity of CuCl obtained after subtracting the phonon contributions is shown in Figure 20(b) and (c). In addition to the -type anomaly at 23 K the magnetic contributions to the heat capacity are characterized by a broad anomaly due to short-range AFM ordering with maximum at 35 K (in and a maximum value of 0.11 J/molK). Our data are close to those of Billerey et al. but deviate markedly from the data published by Stout et al., who found a significantly higher value for the short-range ordering maximum. The total entropy of the magnetic heat capacity amounts to 0.66(2), very close to ln2 expected for the entropy of a spin system with =1/2. Figure 20(c) shows a comparison of the magnetic heat capacity with various model calculations, e.g., the heat capacity of an AFM Heisenberg chain with NN only and with NN and NNN interactions with exchange parameters and the ratio / similar to those discussed for the magnetic susceptibility.Johnston2000 ; Heidrich2006 Both chain models deviate significantly from the experimental data. Surprisingly good agreement of the broad magnetic heat capacity anomaly is obtained for a simple model assuming a spin =1/2 dimer according to

(14) |

with the heat capacity discussed in detail in Ref. Johnston2000, for an exchange constant -117 K. The data are also well reproduced by an anisotropic three-dimensional Heisenberg model on a cubic lattice with in-plane exchange and interplane coupling .Sengupta2003 With -70K and maximum value and position of the short-range ordering anomaly are rather well reproduced. The onset of three-dimensional ordering showing up as a well separated peak well below the short-range ordering maximum appears close to the cooperative ordering temperature observed for CuCl.

## Ix Discussion

Complex compounds containing Cu ions with low-dimensional magnetic behavior are legion and have been widely investigated.DeJongh1974 However, the detailed magnetic behavior of anhydrous Copper(II)-chloride, CuCl, itself, however, has largely remained a mystery until now. From the early susceptibility studies by deHaas and Gorter and Starr et al. it had become evident that CuCl shows unusual AFM properties with a broad maximum in the magnetic susceptibility.DeHaas1931 ; Starr1940 The specific heat of CuCl measured by Stout and Chisholm proved long-range AFM ordering at around 23 K and the excess heat capacity extending well above was analyzed in terms of linear chain behavior. Because of a lack of better theoretical foundation this analysis was done using an Ising-type approach.Stout1962 A neutron diffraction study in order to determine the long-range ordered groundstate has not been carried out so far.

Our study presents a detailed confirmation of the early susceptibility study carried out on polycrystalline samples and twinned single-crystals of anhydrous CuCl. Susceptibility and heat capacity confirm the one-dimensional AFM properties of CuCl and the Nel temperature 23 K as already found by Stout and Chisholm.Stout1962 This rather high Nel temperature as compared to the maximum temperature of the short-range ordering susceptibility maximum occurring at about 70K indicates substantial inter-chain exchange interaction and underlines the conclusion by deJongh and Miedema that CuCl is a rather poor example of a chain structure. DeJongh1974

CuCl consists of ribbons of edge-sharing CuCl square planes with Cu ions. In such systems a spin helicoidal groundstate can set in along the chain direction at low temperature due to the competition between NN ferromagnetic and NNN antiferromagnetic interactions along the Cu chains. The first clear example proven to show such a behavior was LiCuVO that crystallizes with an inverse spinel structure and contains CuO ribbon chains. Helicoidal magnetic ordering with the helix propagating along the chain direction was found by a single-crystal neutron diffraction study.Gibson2004 The frustration scenario with ferromagnetic NN and larger AFM NNN exchange was proven by an inelastic neutron scattering study.Enderle2005 Special interest in the magnetic properties of LiCuVO arose lately from the observation of multiferroic behavior found by Naito et al. and Schrettle et al..Naito2007 ; Yasui2008 ; Schrettle2008 LiCuVO shows a dielectric polarization below the Nel temperature and the polarization can favorably be switched by applying a magnetic field. The origin of this dielectric polarization is due to spin-orbit coupling on the Cu sites, however, the asymmetric charge density distribution necessary for the electric polarization occurs around the O ions.Xiang2007a

LiCuO is another Cu system with CuO ribbon chains exhibiting helicoidal spin order. Masuda2004 ; Masuda2005 ; Gippius2004 ; Drechsler2005 ; Papagno2006 ; Mihaly2006 ; Xiang2007a In LiCuO, in contrast to LiCuVO, the magnitude of the FM NN (-113 meV) exchange is larger than the AFM NNN exchange (71 meV). Ferroelectricity in LiCuO has been shown by Park and coworkers.Park2007

Another related system isostructural with LiCuO, which lately has attracted special attention, is NaCuO. A helicoidal magnetic ground state has been found by neutron powder diffraction and by NMR experiments.Capogna2005 ; Choi2006 ; Drechsler2006 ; Gippius2008 A review of recent systems with helicoidal ground states has been compiled by Drechsler and collaborators.Drechsler2007

All helicoidal systems investigated so far contain oxygen ions to provide superexchange between the Cu moments. Clear proof of the helicoidal ground state by magnetic neutron scattering is available only for a few of them. Evidence for systems containing anions other than O has not been reported yet.

Anhydrous CuCl is therefore the first example for a halide system for which a helicoidal ground state has been proven unequivocally. Our electronic structure calculations give clear theoretical evidence for a scenario of competing FM NN and AFM NNN exchange interactions. Our experimental results for the bulk magnetic properties and the comparison with their predictions for a =1/2 AFM Heisenberg chain with competing FM NN and AFM NNN exchange support the theoretical predictions. Clear proof for the helicoidal AFM ground state comes from the neutron diffraction work on powder and single crystals. We find a magnetic structure indicating a =1/2 helix in the CuCl ribbon chains. The component of the magnetic propagation vector along the -axis amounts to 0.2257(6) and is in best agreement with the result of the electronic structure calculations. The spin-spiral as refined from diffraction data to is found to lie in the plane, consistent with the spin-flop observed in the magnetization data when the magnetic field is oriented along the axis. The refined magnetic structure indicates a magnetic unit cell which is close to being quadrupled as compared to the nuclear cell with neighboring moments enclosing an angle of 81 which leads to a significant reduction of NN exchange energy. With respect to the helix this scenario is analogous to that in LiCuV (propagation vector (0,0.468,0) (Ref. Propvector, ) and NaCuO (propagation vector (0.5,0.227,0). Capogna2005

The magnetic susceptibility data indicate the importance of inter-chain interaction, which according to the electronic structure calculation amounts to 20% - 30% of the intra-chain interaction. Interchain interaction decreases the temperature for the short-range order maximum and provides a better agreement of experimental data and theoretical prediction (Figure 19).

Particularly surprising and striking are the deviations of the short-range ordering contributions to the magnetic heat capacity for CuCl and the heat capacity expected for a =1/2 quantum chain with FM NN and AFM NNN exchange interactions (see Fig. 20(c)). The short-range order contributions match significantly better with a =1/2 dimer singlet-triplet excitation scenario with AFM exchange -117 K corresponding to the exchange found for the NNN interactions along the ribbons. We believe that this coincidence is due to the competing intra-chain competing interactions and the additional AFM inter-chain coupling. The magnetic configuration energy measures the short-range pair correlation function according to (cf. e.g. Ref. Fisher1962, )

(15) |

where the subscript 1 denotes a spin adjacent to 0. Thus and the magnetic heat capacity in a frustrated chain system with coupling to neigboring chains may primarily sense magnetic correlations in a short-ranged cluster in which the configuration energy is determined by the leading AFM exchange to NNN spins.

In summary, anhydrous CuCl shows 1D AFM behavior and long range AFM ordering below a Nel temperature of 23.9 K, below which CuCl adopts an incommensurate magnetic structure (1,0.2257,0.5) with a spin-spiral propagating along and the moments confined in the crystallographic plane. The spin-spiral results from competing FM NN and AFM NNN spin exchange interactions along . Anhydrous CuCl is the first halide quantum system containing CuCl ribbon chains for which a helicoidal magnetic ground state is realized, and is expected to be ferroelectric below 23.9 K.

###### Acknowledgements.

We thank E. Brücher, S. Höhn, S. Lacher and G. Siegle for experimental assistance. The work at North Carolina State University was supported by the Office of Basic Energy Sciences, Division of Materials Sciences, U. S. Department of Energy, under Grant DE-FG02-86ER45259.## References

- (1) W. J. de Haas, C. J. and Gorter, Commun. Leiden 215a (1931).
- (2) C. Starr, F. Bitter, A. F. Kaufmann Phys. Rev 58, 977 (1940).
- (3) J.W. Stout and R.C. Chisholm, J. Chem. Phys. 26, 979 (1962).
- (4) C.G. Barraclough and C.F. Ng, Trans. Faraday Soc. 60, 836 (1964).
- (5) M. G. Banks, Ph.D. thesis, Loughborough University, Loughborough U.K. (2007).
- (6) E. Gmelin, Thermochimica Acta 29, 1 (1979); E. Gmelin and P. Rödhammer, J. Phys. E (Sci. Instrum.) 14, 223 (1981); E. Gmelin and K. Ripka, Cryogenics 21, 117 (1981).
- (7) http://www.ill.eu/d20/home/
- (8) http://www.ill.eu/d10/
- (9) A.F. Wells, J. Chem. Soc. 1947, 1670 (1947).
- (10) P.C. Burns and F.C. Hawthorne, American Mineralogist 78, 187 (1993).
- (11) J. Rodríguez-Carvajal, Physica (Amsterdam) 192B, 55 (1993).
- (12) M. E. Fisher, Philos. Mag. 7, 1731 (1962).
- (13) P. W. Selwood, Magnetochemistry 2nd ed. (Interscience, New York, 1956), p. 78.
- (14) H. Lueken, Magnetochemie (Teubner, Leipzig, 1999).
- (15) M. Takigawa, P. C. Hammel, R. H. Heffner, Z. Fisk, J. L. Smith, and R. B. Schwarz, Phys. Rev. 39, 300 (1989).
- (16) A. Abragam and B. Bleaney, Electron Paramagnetic Resonance of Transition Ions, (Oxford University Press, 1970).
- (17) M. B. Banks and R. K. Kremer, unpublished results.
- (18) D. Billerey, C. Terrier, R. Mainard, and N. Criet, J. Mag. Mag. Mater., 30, 55 (1982).
- (19) D. E. Partin and M. O’Keeffe, J. Solid State Chem. 95, 176 (1991).
- (20) E. S. Itskevich and P. G. Strelkov, Soviet Phys. Doklady 1, 98 (1956).
- (21) K. Takeda, K. Ubukoshi, T. Haseda, and K. Hirakawa, J. Phys. Soc. Jpn. 53, 1480 (1984).
- (22) J. C. LeGuillou and J. Zinn-Jusin, Phys. Rev. Lett. 39, 95 (1977).
- (23) H. Kawamura, J. Phys. Soc. Jpn. 54, 3220 (1985); H. Kawamura, J. Phys. Soc. Jpn. 55, 2095 (1986).
- (24) M. F. Collins, Magnetic Critical Scattering, (Oxford University Press 1989).
- (25) B. D. Gaulin, in: Magnetic Systems with Competing Interactions Ed. H. T. Diep (World Scientific Publishing Co. 1994).
- (26) G. Kresse, J. Hafner, Phys. Rev. B 62, 558 (1993).
- (27) G. Kresse, J. Furthmüller, Comput. Mater. Sci. 6, 15 (1996).
- (28) G. Kresse, J. Furthmüller, Phys. Rev. B 54, 11169 (1996).
- (29) J. P. Perdew, S. Burke, M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
- (30) S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, A. P. Sutton, Phys. Rev. B 57 1505 (1998).
- (31) D. Dai, M.-H. Whangbo, J. Chem. Phys. 114 2887 (2001).
- (32) D. Dai, M.-H. Whangbo, J. Chem. Phys. 118, 29 (2003).
- (33) H. J. Xiang and M.-H. Whangbo, Phys. Rev. Lett. 99, 257203 (2007).
- (34) H. J. Xiang, C. Lee and M.-H. Whangbo, Phys. Rev. B 76, 220411(R) (2007).
- (35) K. Kune, P. Novk, M. Divi and P. M. Oppeneer, Phys. Rev. B 63, 205111 (2001).
- (36) R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 (1993).
- (37) R. Resta, Rev. Mod. Phys. 66, 899 (1994).
- (38) F. Heidrich-Meisner, A. Honecker, T. Vekua, Phys. Rev. B 74 020403(R) (2006).
- (39) For detailed numerical tables see http://www.theorie.physik.uni-goettingen.de/ honecker/j1j2-td/
- (40) R. L. Carlin, Magnetochemistry, (Springer Berlin-Verlag 1986).
- (41) A. Klümper and D. C. Johnston, Phys. Rev. Lett 84, 4701 (2000).
- (42) D. C. Johnston. R. K. Kremer, M. Troyer, X. Wang, A. Klümper, S. L. Bud’ko, A. F. Panchula, and P. C. Canfield, Phys. Rev. B 61, 9558 (2000).
- (43) W. O. Boo and J. W. Stout, J. Phys. E 10, 485 (1977).
- (44) R. C. Chisholm and J. W. Stout, J. Chem. Phys. 36, 972 (1962).
- (45) R. B. Murray, Phys. Rev. 100, 1071 (1955).
- (46) P. Sengupta, A. W. Sandvik, and R. R. P. Singh, Phys. Rev B 68, 094423 (2003).
- (47) L. J. de Jongh and A. R. Miedema, Adv. Phys. 23, 1 (1974).
- (48) B. J. Gibson, R. K. Kremer, A. V. Prokofiev, W. Assmus, and G. J. McIntyre, Physica B 350, e253 (2004).
- (49) M. Enderle, C. Mukherjee, B. Fåk, R. K. Kremer, J.-M. Broto, H. Rosner, S.-L. Drechsler, J. Richter, J. Malek, A. Prokofiev, W. Assmus, S. Pujol, J.-L. Raggazzoni, H. Rakoto, M. Rheinstädter, and H. M. Rnnow, Europhys. Lett. 70, 237 (2005).
- (50) Y. Naito, K. Sato, Y. Yasui, Y. Kobayashi, and M. Sato, J. Phys. Soc. Jpn. 76, 023708 (2007).
- (51) F. Schrettle, S. Krohns, P. Lunkenheimer, J. Hemberger, N. Büttgen, H.-A. Krug von Nidda, A. V. Prokofiev, and A. Loidl, Phys. Rev. B 77, 144101 (2008).
- (52) Y. Yasui, Y. Naito, K. Sato, T. Moyoshi, M. Sato, and K. Kakurai, J. Phys. Soc. Jpn. 77, 023712 (2008).
- (53) T. Masuda, A. Zheludev, A. Bush, M. Markina, and A. Vasiliev, Phys. Rev. Lett. 92, 177201 (2004).
- (54) T. Masuda, A. Zheludev, B. Roessli, A. Bush, M. Markina, and A. Vasiliev, Phys. Rev. B 72, 014405 (2005).
- (55) A. A. Gippius, E. N. Morozova, A. S. Moskvin, A. V. Zalessky, A. A. Bush, M. Baenitz, H. Rosner, and S.-L. Drechsler, Phys. Rev. B 70, 020406(R) (2004).
- (56) S.-L. Drechsler, J. Málek, J. Richter A. S. Moskvin, Phys. Rev. Lett. 94, 039705 (2005).
- (57) M. Papagno, D. Pacilé, G. Caimi, H. Berger, L. Degiorgi, and M. Grioni, Phys. Rev. B 73, 115120 (2006).
- (58) L. Mihály, B. Dóra, A. Ványolos, H. Berger, and L. Forró, Phys. Rev. Lett. 97, 067206 (2006).
- (59) T. Masuda, A. Zheludev, A. Bush, M. Markina, and A. Vasiliev, Phys. Rev. Lett. 92, 177201 (2004).
- (60) S. Park, Y. J. Choi, C. L. Zhang, and S-W. Cheong, Phys. Rev. Lett. 98, 057601 (2007).
- (61) L. Capogna, M. Mayr, P. Horsch, M. Raichle, R. K. Kremer, M. Sofin, A. Maljuk, M. Jansen, and B. Keimer, Phys. Rev. B 71, 140402(R) (2005).
- (62) K.-Y. Choi, V. P. Gnezdilov, P. Lemmens, L. Capogna, M. R. Johnson, M. Sofin, A. Maljuk, M. Jansen, and B. Keimer, Phys. Rev. B 73, 094409 (2006).
- (63) S.-L. Drechsler, J. Richter, A. A. Gippius, A. Vasiliev, A. A. Bush, A. S. Moskvin, J. M lek, Yu. Prots, W. Schnelle and H. Rosner, Europhys. Lett. 73, 83 (2006).
- (64) A. A. Gippius, A. S. Moskvin, and S. L. Drechsler, Phys. Rev. B 77, 180403(R) (2008).
- (65) S.-L. Drechsler, J. Richter, R. Kuzian, J. Mlek, N. Tristan, B. Büchner, A.S. Moskvin, A.A. Gippius, A. Vasiliev, O. Volkova, A. Prokofiev, H. Rakotoh, J.-M. Broto, W. Schnelle, M. Schmitt, A. Ormeci, C. Loison, and H. Rosner, J. Magn. Magn. Mater. 316, 306 (2007).
- (66) Note that in LiCuV the nuclear cell contains two Cu atoms along due to the VO-tetrahedra alternating in their orientation. Accordingly, the propagation vector (0,0.468,0) indicates an approximate doubling of the magnetic cell which again contains four =1/2 entities.