Linear magnetoelastic coupling and magnetic phase diagrams of the buckled-kagomé antiferromagnet [Formula: see text].

Single crystals of Cu[Formula: see text]Bi(SeO[Formula: see text])[Formula: see text]O[Formula: see text]Cl were investigated using high-resolution capacitance dilatometry in magnetic fields up to 15 T. Pronounced magnetoelastic coupling is found upon evolution of long-range antiferromagnetic order at [Formula: see text] [Formula: see text] K. Grüneisen analysis reveals moderate effects of uniaxial pressure on [Formula: see text], of 1.8(4) K/GPa, [Formula: see text] K/GPa and 0.33(10) K/GPa for [Formula: see text], b, and c, respectively. Below 22 K Grüneisen scaling fails which implies the presence of competing interactions. The structural phase transition at [Formula: see text] [Formula: see text] K is much more sensitive to uniaxial pressure than [Formula: see text], with strong effects of up to 27(3) K/GPa ([Formula: see text]). Magnetostriction and magnetization measurements reveal a linear magnetoelastic coupling for [Formula: see text] below [Formula: see text], as well as a mixed phase behavior above the tricritical point around 0.4 T. An analysis of the critical behavior in zero-field points to three-dimensional (3D) Ising-like magnetic ordering. In addition, the magnetic phase diagrams for fields along the main crystalline axes are reported.

Introduction

Frustrated magnetism has been a highly active research field in the past three decades[1]. Besides simple corner-sharing triangular lattice geometries, the Kagomé, pyrochlore, and hyperkagome lattices offer potential playgrounds of strongly geometrically frustrated magnetism[2-4]. Among them, the ideal Kagomé Heisenberg antiferromagnet is the most prominent realization of a system in which macroscopic ground state degeneracies, i.e., a (quantum) spin liquid state, are expected to prevent the evolution of any long-range order[5]. Other phenomena of frustrated magnetism include fractionalized magnetization plateaus, chiral and helical spin arrangements, as well as spin glass, spin nematic, and spin ice behaviors[1]. A number of geometrically frustrated systems such as [6], [7], and [8,9] were found to additionally exhibit multiferroic behavior. From a technological perspective multiferroics – materials combining more than one ferroic property such as ferromagnetism, ferroelectricity and ferroelasticity[10,11] – are especially sought-after for the control of magnetism via electric fields, with the promise of substantially lower energy consumption than manipulating magnetic states via magnetic fields[12]. Potential applications range from ultra-low power logic-memory[13] to radio- and high-frequency devices, including electric field-tunable radio-frequency/microwave signal processing, magnetic field sensors, magnetoelectric random access memory (MERAM)[14] and voltage-tunable magnetoresistance[15,16].

CuBi(SeO)OCl, eponym of the francisite family[17-22], crystallizes in a special buckled realization of the Kagomé lattice[23] and exhibits multiferroic properties at low temperatures[24,25]. Cu ions ( = 1/2) situated at two different crystal sites, Cu1 and Cu2, are the magnetic centers forming the Kagomé lattice in the ab plane. Both Cu1 and Cu2 ions are found in a square planar coordination with Cu-O bond lengths of 1.933 Å to 1.978 Å. These plaquettes around Cu1 and Cu2 sites are non-parallel with respect to each other. Along the c axis Cu ions are connected by long Bi–O bonds with a bond length of about 2.8 Å. Below   K CuBi(SeO)OCl (in short: CBSCl) develops an A-type antiferromagnetic order along the c axis of spins aligned ferromagnetically (FM) in the ab plane. Cu2 spins are aligned strictly (anti)parallel with the c axis whereas Cu1 spins are canted towards the b axis by 59(4)[25]. This peculiar magnetic order arises from the competition of large FM nearest neighbor exchange interactions ( and , see Fig. S2[26]]) on the order of  K to  K and AFM next-nearest neighbor interactions () on the order of 60 K, in conjunction with small exchange couplings along the c direction ( and ,  K), Dzyaloshinskii-Moriya interactions and a symmetric anisotropic exchange component[25-27]. Furthermore, a linear magnetoelectric coupling has been observed below [24,25]. At higher temperatures around 120 K a structural phase transition from an orthorhombic Pmmn to an orthorhombic and nonpolar, possibly antiferroelectric (AFE), Pcmn space group occurs upon cooling[25,28,29].

In this paper we investigate the interplay of the lattice and spin degrees of freedom in CBSCl with high-resolution thermal expansion and magnetostriction as well as magnetization measurements. From these measurements we construct the complete magnetic phase diagram for the first time, and quantify the magnetoelastic coupling by the effects of uniaxial pressure on the phase boundaries. Furthermore, magnetostriction and magnetization measurements reveal a linear magnetoelastic coupling below . Finally, a critical scaling analysis suggests that three-dimensional Ising-like magnetic ordering occurs at , whereas magnetic correlations above are constrained to the ab plane.

Results and discussion

Thermal expansion

Zero-field

Thermal expansion measurements in zero field show pronounced anomalies at   K and   K (Fig. 1). these anomalies provide evidence for the presence of significant magnetoelastic coupling in CBSCl. In the thermal expansion coefficients the anomalies have -like shapes, signaling continuous phase transitions at both and . Notably, the transition at shows a jump superimposing the -like behavior for the a and c axis, but not for the b axis. The length changes observed in CBSCl (Fig. 1b) are highly anisotropic. While the a axis shrinks above upon warming up to , the c axis strongly expands in this temperature regime. Expansion along the b axis, in contrast, is much smaller than for the other two axes. The volume expansion is positive in the whole measured temperature regime (see supplement, Fig. S3).

Figure 1

Zero-field thermal expansion: (a) Thermal expansion coefficients and (b) relative length changes in zero-field for the crystallographic a, b, and c axis of CBSCl. Insets show a magnification of the low temperature regime around with zero y-values indicated by horizontal dashed lines. Vertical dotted lines in (b) mark and .

Grüneisen analysis

The magnetoelastic coupling in CBSCl can be quantified by means of the uniaxial pressure dependence of . This pressure dependence may be derived from the magnetic contributions to thermal expansion and specific heat using the Ehrenfest relation which is valid for continuous phase transitions[30,31]: is the (magnetic) Grüneisen ratio and the molar volume, which for CBSCl is  m/mol[32]. The magnetic contributions to the thermal expansion, , and specific heat, , in Fig. 2 were obtained by subtracting a phononic background (see supplement, Fig. S4) estimated by a sum of Debye and Einstein functions (for more details see supplement). In order to assess the validity of a constant in Eq. (1), in Fig. 2 the ordinates are scaled for an optimal overlap of and around . Above and down to about 22 K (24 K for c) the data scale well within the experimental errors which implies a single dominating energy scale in this temperature regime[30]. Below 22 K, however, a difference in the behavior of and is clearly visible for all axes as well as the volume. This difference indicates competing interactions arising below . A discussion of possible sources for this behavior is given in the supplemental material.

Figure 2

Grüneisen Scaling: Comparison of magnetic contributions to the specific heat (red circles, left axes) and thermal expansion (empty black circles, right axes) along the (a) a axis, (b) b axis, (c) c axis, and (d) for the volume expansion.

The magnetic Grüneisen ratio in Eq. (1) is obtained from the ratio of the scaled ordinates in Fig. 2. This results in  mol/J,  mol/J, and  mol/J, which yields uniaxial pressure dependencies at of 1.8(4) K/GPa (),  K/GPa () and 0.33(10) K/GPa (). Long-range antiferromagnetic (AFM) order is thus strengthened by pressure and c whereas influences the competition of interactions in a way that long-range order is suppressed. From a qualitative point of view, increase of upon in-plane pressure which likely increases the lattice spacing along the c direction suggests that the perpendicular exchange couplings () are not crucially driving long-range AFM order. We also speculate that uniaxial pressure () changes the Cu1-O-Cu1 bond angle further towards (away from) 90, thereby increasing (decreasing) ferromagnetic exchange while becomes smaller for . Our results, i.e., and hence indicate a crucial role of for stabilizing A-type AFM order.

The analysis of the structural phase transition yields much larger values for the uniaxial pressure dependencies in comparison with . In-plane pressure strongly suppresses the structural phase transition ( K/GPa,  K/GPa) whereas pressure applied out-of-plane enhances even more strongly (/GPa). The Grüneisen ratios corresponding to these values are derived in a different manner than for , by area-conserving – and for entropy-conserving – interpolations of the jump heights of the respective quantities at (see supplement, Fig. S5). Table 1 summarizes the pressure dependence of the transition temperatures in CBSCl.

Table 1

Pressure dependencies of the transition temperatures calculated from the magnetic Grüneisen ratios at , as well as from the jumps and at , for the three main crystallographic axes. The hydrostatic pressure dependence is calculated as the sum of all three uniaxial ones and denoted as “volume”.

	γ_i,mag	∂T_N/∂p_i	γ_i	∂T_S/∂p_i
	(10^-7 mol/J)	(K/GPa)	(10^-7 mol/J)	(K/GPa)
a axis	5.0	1.8±0.4	–11.9	–19±6
b axis	–1.8	–0.62±0.15	–3.5	–5.7±1.1
c axis	0.93	0.33±0.10	16.9	27±3
volume	4.0	1.4±0.3	–	2.3±1.0

We find that both and are enhanced by hydrostatic pressure. Previous studies investigated the effects of chemical doping on the physical properties of CBSCl by substituting the lone-pair Bi site by Y[18] or lanthanide elements[19,21], exchanging the halide ion for Br or I[32], or Te doping the Se site[24]. Hydrostatic pressure applied to CBSCl was shown to enhance by about 1 K/GPa[24] which is in line with our results. For the pressure dependence of no literature data are available. The structural phase transition is absent in the compounds where Cl is exchanged for Br and I due to their larger ionic radii, and was not reported for any of the other doped versions of francisite.

Thermal Expansion for

The c axis is the easy-axis in CBSCl and shows a metamagnetic transition in magnetic fields , rendering it of special interest. Application of a magnetic field suppresses the Néel transition to lower temperatures (Fig. 3b,c). Up to a field of 0.5 T an increase in the thermal expansion coefficient signals field-induced enhancement of the pressure dependence of the entropy changes (). At even higher fields the character of the phase transition changes from a continuous to a discontinuous transition. This change is evidenced by a broadening of the -like feature to a plateau-like behavior between two jumps (up and down) in as the temperature is increased. A magnetic field of 1 T fully suppresses AFM order while a ferromagnetic phase appears. Notably, the exact same behavior described for at low fields is also visible in the Fisher specific heat[33], , a derivative of the static magnetic susceptibility (Fig. 3a). This behavior and especially the intermediate phase emerging in an applied field will be discussed in more detail below.

Figure 3

Fisher’s specific heat and thermal expansion up to  T: (a) Fisher’s specific heat calculated from the static magnetic susceptibility, (b) thermal expansion coefficient , and (c) relative length changes, for 0 T  1 T. The inset in (a) presents a magnification of the low temperature region. Vertical bars in (c) mark the relative length changes from 0 T to 0.8 T (olive green) and 1 T (orange), respectively, obtained from magnetostriction measurements.

Magnetostriction for

Application of magnetic fields drives the system into a ferromagnetic (FM) phase as indicated by sharp jumps in the magnetization (Fig. 4c,d). The metamagnetic transition is characterized by sharp jumps in the macroscopic length as seen in the magnetostriction measurements along the c axis (Fig. 4a,b). In both magnetization and length changes, a broad hysteresis region becomes visible upon ramping down the applied field at 2 K (0.86 T to 0.77 T), indicating the first-order nature of the transition. The FM state features the magnetization of Cu thereby indicating that the spins are not fully aligned but presumably canted. At 10 K and above a linear increase is visible in between a regime where at low fields and a constant value above the phase transition. The hysteresis region shrinks as the temperature is increased towards and does not extend over the whole regime of linear increase.

Figure 4

Linear magnetoelastic coupling: Magnetostriction and magnetization. Comparison of relative length changes (a, b) and isothermal magnetization (c, d) for . Note that (b) and (d) show a magnification of the transition region and only data for selected temperatures. Hysteresis is marked by colored areas. Magnetostriction down-sweeps in (b) are shifted by 6 mT to correct for the remanent field of the magnet.

Mixed-phase for

The linear increase and hysteresis observed in points to a mixed phase between the low-field AFM ordered phase and the field-induced FM state. In this mixed intermediate phase both AFM and FM domains are present and the ratio FM:AFM linearly increases with the applied magnetic field, as the demagnetizing field is overcome and AFM regions align with the field. Such mixed-phase behavior is common for metamagnets and has also been observed in the brother compound CuBi(SeO)OBr[34]. There, the mixed-phase was found to exhibit broadband absorption with excitations extending over at least ten decades of frequency[34]. A similar behavior is naturally also expected in CuBi(SeO)OCl.

Linear magnetoelastic coupling

On top of the mixed-phase behavior, a remarkable direct proportionality between the magnetization and the relative length changes , i.e., a linear magnetoelastic coupling  [/Cu], is observed (Fig. 4) at . A comparison of the jump heights at the metamagnetic transition, normalized to the jump heights at 2 K, illustrates this observation even more clearly (Fig. 5). The temperature evolution of both and are proportional to each other and obeys an order parameter-like behavior.[35], linear magnetostriction was first observed in by Borovik-Romanov and co-workers[36]. Among many others, antiferromagnetic examples include hematite (-FeO)[37], dysprosium orthoferrite[35] (), and more recently the multiferroic [38]. In all of these compounds linear magnetostriction is closely related to the presence of antiferromagnetic domains. The magnetic point group of CBSCl, , is not among the point groups in which a linear magnetostriction is expected to be dominant[39], which does not mean, however, that it may not, as observed, be dominant.

Figure 5

Linear magnetoelastic coupling: Jump sizes. Comparison of the jumps in magnetization (left ordinate) and relative length changes (right ordinate) at the metamagnetic transition normalized by their value at 2 K.

In addition to the linear magnetoelastic coupling, it was argued that the point group is also compatible with a linear magnetoelectric coupling[25], and indeed, a linear magnetoelectric coupling was shown experimentally[24].

In-plane magnetic correlations far above

Thermal expansion measurements at 15 T reveal a large temperature regime above where significant in-plane magnetostriction is present (see supplement, Fig. S7). Notably, magnetostriction along the c axis is negligible above , whereas for the a and b axis there are pronounced magnetic field effects up to above 80 K. For the a axis the length changes are about one magnitude larger than for the b axis. Qualitatively, application of a magnetic field acts in the same way as the evolution of long-range order below : the a axis shrinks whereas the b axis elongates. The strong difference between magnetostriction in the ab plane and along the c axis evidences strong in-plane magnetic correlations above , in line with the large and competing FM and AFM exchange couplings on the order of 60 K to 70 K and the small inter-plane couplings on the order of up to 2 K[26,27]. Negligible magnetostriction along the out-of-plane direction, signals that the exchange interactions J are nearly independent of pressure .

Phase diagrams

The low-temperature magnetic phase diagrams for CBSCl as measured by thermal expansion (c axis) and magnetization (all axes) are presented in Fig. 6. Corresponding magnetization measurements are shown in the supplement (Fig. S9, S10, and S11). The saturation field of  T is extrapolated from the data for , whereas  T and  T were obtained by measurements at 2 K for and , respectively. For the c axis, a change in behavior from a continuous phase transition at zero- and low fields to a discontinuous transition above 0.4 T signals the presence of a tricritical point around 0.4 T. Moreover, a large magnetic field hysteresis at the metamagnetic transition is visible at low temperatures, which decreases as the temperature is increased.

Figure 6

Magnetic Phase Diagrams: (a) Low-temperature magnetic phase diagram for . The hatched area displays the hysteresis region visible in isothermal magnetization and magnetostriction curves. Where hysteresis and the mixed phase overlap, white lines indicate phase boundaries obtained from up-sweep data. (b) Low-temperature magnetic phase diagram for (black boundary) and (white boundary) constructed from magnetization measurements. The phase boundary ( T) is obtained by scaling ().

Quantitative analysis of the phase boundaries for

By measuring the jump-like discontinuities in the magnetization and length at the metamagnetic phase boundary , the uniaxial pressure dependencies of the critical field and the critical temperature can be derived from the Clausius-Clapeyron equation and related thermodynamic equations[30]where and are the observed jumps in the length and the magnetic moment. is calculated from two polynomial fits in different temperature regimes to the phase boundary such that the changes in entropy, , can be calculated from Eq. (3). The pressure dependence of the critical field is small and positive, on the order of 80(6) mT/GPa at all temperatures. Derived changes of the transition temperature at the metamagentic (AFM to FM) phase transition under pressure are large at low temperatures ( K/GPa at 2 K) and decrease strongly as the temperature is increased ( K/GPa at 24 K). Correspondingly, the entropy related to the phase transition increases from 30 mJ/(mol K) at 2 K to 360(130) mJ/(mol K) at 24 K. This analysis shows that uniaxial pressure changes the exchange couplings in such a way that the AFM phase is stabilized with respect to both the FM and the paramagnetic (PM) phase. An overview of calculated values is given in the supplement in Table S2.

Critical scaling

In the vicinity of a critical point the specific heat is expected to behave aswhere is the reduced temperature[40]. As seen in the Grüneisen scaling in Fig. 2, the magnetic contributions to the specific heat and the thermal expansion around can be scaled to each other. Therefore, we can safely assume that a critical scaling of close to , i.e., above the shoulder-like feature at low temperatures ( for and for ), with the expression from Eq. (5) is allowed. Although we use the magnetic contributions, and , for the fitting, there is an uncertainty on the phonon background correction. Therefore, and to account for possible further contributions on top of the critical behavior, we initially used the canonical expression (similar to the one in Ref.[41])where “” (“−”) denotes fitting parameters for (). It turned out, however, that the data can be fitted very well even when setting and .

Depending on the critical exponent and the ratio of the amplitudes, , the critical behavior can be categorized by one of many universality classes. The most well-known universality classes for a d-dimensional lattice and an order parameter of dimensionality D are the 3D Heisenberg model (, ), the 3D XY model (, ), the 3D Ising model (, ) and the 2D Ising model (, ). Calculations for these different models predict , (3D Heisenberg), , (3D XY), , (3D Ising)[42] , as well as (2D Ising)[43].

The fits to the critical region around , roughly around , for the linear thermal expansion coefficients, the volume expansion coefficient and the specific heat are shown in Fig. 7.S3). The two main results which stand out from the fits are that (1) the critical exponents are all positive and lie between and 0.23(4), and (2) the ratio of the amplitudes for and , , lies in the range from to 0.46(3) for the thermal expansion data, and for the specific heat. Considering these two main results, values around 0.11 and ratios close to 0.5 – but far away from 1 or even 1.5 – enables us to clearly categorize the transition around as a transition with a one-dimensional order parameter on a three-dimensional lattice, i.e., a 3D Ising-type transition. This result is in line with the strong anisotropy visible in the magnetization and the phase diagrams of CBSCl as well as the strict (uniaxial) alignment of the Cu2 spins along the c axis.

Figure 7

Critical scaling analysis near . Semi-logarithmic plots of the magnetic contributions to (a–c) the thermal expansion coefficients , , (d) the magnetic volume expansion coefficient and (e) the specific heat versus reduced temperature , for  K. Black squares (red circles) mark data for (). Red and black solid lines are fits to the data according to Eq. (5), vertical dashed lines indicate the fitting window. Values for with no error given were fixed for the fitting.

Previous calculations successfully described many of the properties of (X = Cl,Br) by an effective 3D spin model, i.e., by Heisenberg spins combined with Dzyaloshinskii-Moriya (DM) interactions[26] and symmetric anisotropic exchange interactions[25]. Also, neutron diffraction measurements suggested that the direction of the Cu spins is constrained to the bc plane[25]. For the brother compound CuBi(SeO)OBr (CBSBr), on the other hand, a crossover from a 2D XY () to a 3D () character was suggested from neutron diffraction measurements near  = 27.4 K[44]. Deviations from the critical behavior are also observed in our data below , roughly for . Together with the observation of dispersionless magnon modes along the c direction at 2 K[25] this may point to 2D correlations evolving below . However, whether or not the low-temperature behavior is caused by a 3D-to-2D crossover can not be concluded unambiguously from our data and necessitates further neutron or NMR studies. Therefore, while a crossover to a spatial 2D behavior as reported for CBSBr cannot be neither concluded nor excluded from our data on CBSCl, our results clearly point to a 3D Ising-like behavior in the critical region around .

Conclusions

The buckled-kagomé antiferromagnet CuBi(SeO)OCl was investigated by thermal expansion, magnetostriction, magnetization and specific heat measurements. Its highly anisotropic lattice changes in temperature and in magnetic field reveal two phase transitions at  K and  K. The low-temperature and low-field AFM phase for , which experiences a field-driven metamagnetic transition to an FM phase, exhibits linear magnetoelastic coupling and a sizable mixed phase between the AFM and field-induced FM phases. Uniaxial pressure stabilizes the AFM phase at the expense of the surrounding FM and paramagnetic (PM) phases. is not affected by magnetic fields but strongly suppressed by uniaxial pressure along the in-plane directions, whereas strongly enhances it by about 27 K/GPa. The critical behavior in the vicinity of is in line with calculations for the 3D Ising model. Furthermore, while in-plane magnetic correlations extend to temperatures far above the c axis shows no significant magnetostriction and effects of a magnetic field on the thermal expansion above .

The francisite CuBi(SeO)OCl is an exciting compound with strong physical effects from both lattice and spin degrees of freedom at different temperatures. Our study provides additional evidence that the low-temperature phase of CuBi(SeO)OCl presents a versatile playground for investigating multiferroic effects, combining both a linear magnetoelastic and linear magnetoelectric coupling.

Methods

Crystal growth

Single crystals of CuBi(SeO)OCl were grown by the chemical vapor transport method as reported in Ref. [28]. Different single crystals of sizes  mm ( mg),  mm ( mg) and  mm ( mg) were used for our studies.

High-resolution dilatometry

High-resolution capacitance dilatometry measurements were performed in two three-terminal high-resolution capacitance dilatometers from Kuechler Innovative Measurement Technology with a home-built setup placed inside a variable temperature insert of an Oxford magnet system[45,46]. Linear thermal expansion coefficients   were derived for all crystallographic axes, i.e., , in a temperature range from 2 K to 200 K, with magnetic fields up to 15 T applied along the measurement direction. Field-induced length changes were measured at various fixed temperatures between 1.7 K and 200 K and the magnetostriction coefficients   were derived.

Magnetization

Magnetization measurements were performed in a Magnetic Properties Measurement System (MPMS3, Quantum Design) up to 7 T and in a Physical Property Measurement System (PPMS-14, Quantum Design) in fields up to 14 T. A rotatable sample holder was used in the MPMS3 for measurements perpendicular to the c axis in order to determine the a and b axis for further magnetization measurements. Laue XRD measurements confirmed the orientation and quality of the crystals (see supplement, Fig. S1).

Specific heat

Specific heat measurements were performed on a PPMS calorimeter using a relaxation method, on a sample of  mg.

Below is the link to the electronic supplementary material.

Supplementary Information.