Magnetically-induced ferroelectricity in the (ND4)2[FeCl5(D2O)] molecular compound
Jose Alberto Rodriguez-Velamazan, Oscar Fabelo, Angel Millan, Javier Campo, Roger D. Johnson, Laurent Chapon
11 Magnetically-induced ferroelectricity in the (ND ) [FeCl (D O)] molecular compound.
José Alberto Rodríguez-Velamazán,
Óscar Fabelo, Ángel Millán, Javier Campo, Roger D. Johnson, Laurent Chapon. Instituto de Ciencia de Materiales de Aragón (ICMA), CSIC – Universidad de Zaragoza, 50009 Zaragoza, Spain. Institut Laue-Langevin, 38042 Grenoble Cedex 9, France. Univ. Oxford, Dept Phys, Clarendon Lab, Oxford OX1 3PU, England
The number of magnetoelectric multiferroic materials reported to date is scarce, as magnetic structures that break inversion symmetry and induce an improper ferroelectric polarization typically arise through subtle competition between different magnetic interactions. The (NH ) [FeCl (H O)] compound is a rare case where such improper ferroelectricity has been observed in a molecular material. We have used single crystal and powder neutron diffraction to obtain detailed solutions for the crystal and magnetic structures of (NH ) [FeCl (H O)], from which we determined the mechanism of multiferroicity. From the crystal structure analysis, we observed an order-disorder phase transition related to the ordering of the ammonium counterion. We have determined the magnetic structure below T N , at 2K and zero magnetic field, which corresponds to a cycloidal spin arrangement with magnetic moments contained in the ac -plane, propagating parallel to the c -axis. The observed ferroelectricity can be explained, from the obtained magnetic structure, via the inverse Dzyaloshinskii-Moriya mechanism. This work is in honor of Prof. Fernando Palacio on the occasion of his retirement
Introduction
Multiferroic materials, compounds presenting at least two ferroic- orders simultaneously, have attracted considerable interest due to the possible cross control of physical properties. In particular, magnetoelectric multiferroics where an electrical polarization is induced by a magnetic transition, the so-called type-II multiferroics, have been extensively studied for their fundamental properties as well as for their potential applications in memories, sensors, transducers, etc.[1] The complex magnetic ordered state in these improper (or pseudo-proper) ferroelectric materials [2] is stabilized by competitions between different exchange couplings which break inversion symmetry and induce a ferroelectric polarization. The most studied multiferroics of this type are the “cycloidal” magnets such as RMnO [3] MnWO [4] or CoCr O [5], where ferroelectricity emerges as a consequence of a cycloidal modulation of the magnetic moments. The coupling mechanism between magnetization and ferroelectricity has been explained by the so-called inverse Dzyaloshinskii-Moriya (DM) effect [6] or spin current mechanism, that predict that the electric polarization, P , is proportional to r ij × ( S i × S j ), where r ij is the vector connecting the nearest spins, S i and S j , and ( S i × S j ) is the so-called spin-chirality vector.[7] One of the approaches to the design of materials combining different properties is the association of different building blocks that carry the different properties of interest. When a ferroelectric material is wanted, the use of ammonium molecules substituting monovalent counterions (i.e. Na, K, Rb or Cs) in inorganic networks has been proved as a successful approach in order to induce a switchable ferroelectric transition. This is the case of peroxychromates of general formula M CrO ,[8] where the substitution of the alkali metals by an ammonium molecule produces a system with a well-defined ferroelectric phase transition.[9] A similar approach has been used in metal-organic compounds, where the combination of two different networks with different properties has been shown to be an excellent strategy to design new multiferroic materials.[10] Nevertheless, in all these examples the existence of ferroelectric or antiferroelectric transitions is due to the occurrence of an order−disorder phase transition where the hydrogen bonds between the guest molecule and the host-framework play a central role. No interplay between the magnetic and electric orders has been observed in these cases, and therefore improper ferroelectricity has never been described for this type of compounds. However, a recently published work focused on ammonium pentacloroaquaferrate(III), with formula (NH ) [FeCl (H O)], has characterized this molecular compound as a new multiferroic material, where magnetic and electric properties are connected.[11] The heat capacity of the deuterated sample, of formula (ND ) [FeCl (D O)] (see Figure S1), shows a behavior equivalent to that of the hydrogenated sample reported by Ackermann et al ,[11] presenting two closely spaced peaks in the low temperature region at ca. 6.9 and 7.2 K.[11,12] At high temperature ( ca
79 K), this material reveals a third signal in the heat capacity, associated to a suspected structural phase transition, although the previous X-ray diffraction measurements above and below this transition did not show any significant changes.[11] The single crystal susceptibility measurements indicate an XY anisotropy with the ac -plane being the magnetic easy plane. The spontaneous electric polarization obtained integrating the pyroelectric current shows that the system becomes ferroelectric below T FE ~ 6.9K at zero magnetic field, corresponding with the low temperature cusp in the heat capacity. The polarization has two components, the main one along the a -axis with a value of 3 µ Cm −2 and a secondary one, with a value ten times lower, along the b -axis. Thus, the spontaneous electric polarization lies in the ab- plane. In this paper, we investigate the temperature dependence of the crystal structure using X-ray and neutron diffraction measurements (on hydrogenated and deuterated samples respectively). Below T N , (at 2 K) we have determined unambiguously the magnetic structure at zero magnetic field. A simple model based on a set of four exchange coupling interactions has been proposed to explain the magnetic ground state. We present also the temperature evolution of some selected magnetic reflections in order to evaluate the possible spin-reorientation above 6.9 K. Finally, we propose the possible multiferroicity mechanism in (ND ) [FeCl (D O)] compound.
Experimental methods
Materials
Deuterated reagents and solvents used in the synthesis were purchased from commercial sources and used without further purification. The synthesis of non-deuterated samples has been carried out with the same protocol described below but with hydrogenated reagents.
Synthesis of (ND ) [FeCl (D O)]
Single crystals of (ND ) [FeCl (D O)] of suitable size for neutron scattering were prepared by the seeded growth technique, by slow evaporation of saturated solutions at 50 ºC. The seeds were formed by cooling of hot saturated solutions. The resulting compound was directly checked by single crystal X-ray diffraction and the results were in agreement with the unit cell parameters and the crystal structure previously reported. [11-12]
X-ray Single Crystal Refinements
Data collections of (NH ) [FeCl ·H O] were carried out at room temperature (HT) and 50 K (LT) in an Agilent SuperNova X-ray m -focus source equipped with a helium cryostream, using Mo-K α radiation ( λ = 0.71073 Å). The reflections were indexed, integrated and scaled using CrysAlis Pro program.[13] The structures of and were solved by direct methods using the SHELXS97 program at the space group Pnma (No. 62). All non-hydrogen atoms were refined anisotropically by full-matrix least-squares technique based on F using SHELXL97.[14] The hydrogen atoms were positioned geometrically and refined using the difference electron density map applying DFIX soft-constrains, giving rise to very close models for both temperatures; therefore no clear signal of the structural phase transition was observed. The only observed change comparing the HT with the LT measurements is the cell volume compression, which is around 2 %, value which is compatible with the thermal expansion effect. The final geometrical calculations and the graphical manipulations were carried out with PARST97,[15] PLATON[16] and DIAMOND[17] programs. Neutron Diffraction Measurements.
Single-crystal neutron diffraction data were collected at the hot-neutrons four-circle diffractometer D9 at Institut Laue Langevin (ILL, Grenoble, France) with a wavelength of 0.8322(1) Å obtained with a copper (113) monochromator operating in transmission geometry. D9 diffractometer is equipped with a 2D detector of 6 x 6 cm (32 x 32 pixels) that allow us a reciprocal space survey. The crystal was mounted onto specific aluminum sample holders that produce a low background and sealed into a closed-cycle cryostat. In order to determine the crystal structure above and below the nuclear phase transition (79 K) data collections consisting in a combination of omega- and omega-2theta-scans of each individual reflection where carried out at 100 and 2 K, yielding 4900 and 2710 independent reflections respectively. The cell parameters where obtained at both temperatures and the crystal structure refined at 100 K and accurately determined below the phase transition. Due to the occurrence of a magnetic propagation vector different of k = (0,0,0), nuclear and magnetic reflections are not overlapped below the magnetic phase transition (T N ca 7.25 K) and therefore both the magnetic and nuclear phases can be measured independently. The magnetic structure was determined at 2 K from 127 independent reflections. Additionally, the evolution of two different magnetic reflections [(0,0,-1- k z ) and (0,-3, k z )] was followed in the temperature range from 2 to 12 K. The main crystallographic data are summarized in Table 1. The program RACER[18] was used to integrate the omega- and omega-2theta-scans and to correct them for the Lorentz factor. The crystal attenuation corrections were performed with a prism model using DATAP program,[19] with an estimated total neutron absorption coefficient of 0.057 cm -1 .[20] Additional high-resolution neutron powder diffraction measurements at D2B instrument (ILL) were performed at 45 K in order to verify the structural transformation at 79 K. The sample was contained in a Ø 6 mm cylindrical vanadium-can inside an Orange Cryostat. The neutron diffraction pattern was collected using 1.5642 Å wavelength. The refinements of the nuclear structures were performed using the programs SHELX[14] and FullProf Suite.[21] The magnetic structure was treated with the program FullProf using the integrated intensities obtained from the single-crystal measurements. The models for the magnetic structure were deduced from the symmetry analysis techniques implemented in the program BasIreps included in the FullProf Suite.[21] The magnetic moments were localized on the Fe(III) atoms and, in the final refinement, the real and imaginary parts of the Fourier coefficients of the magnetic moments were constrained to have the same magnitude with directions described by spherical angles. The nuclear and magnetic contributions to the diffraction intensities were treated as two separate patterns, with the magnetic phase described in space group P -1 (No. 2), using the magnetic form-factor curve of Fe(III) and the scale factor obtained from the refinement of the nuclear phase at the same temperature. Results and discussion
Neutron Studies Crystal structure at 100 K
The title compound is built up from one [FeCl (D O)] ion and two ND counterions, and the crystal structure is held together by an extensive network of H-bonds in concurrence with the ionic interaction (see Figure 1). The KPI packing index gives a percent of filled space of 81.5,[22] therefore there is no accessible space in the network for solvents. The crystallographic properties of the two phases (above and below the structural phase transition), including space groups, lattice parameters as well as some details of the data refinement, are shown in Table 1. Formula [ND ] [Fe III Cl (D O)]
Empirical
Formula Cl D FeN M r (g·mol -1 ) Temperature (K)
Inst. l (Å) D9, 0.83220
D9, 0.83220
Crystal system
Orthorhombic
Monoclinic
Space group (No.)
Pnma (62) P / a (14) Crystal size (mm) a (Å) b (Å) c (Å) a (°) b (°) g (°) V (Å ) Z 4 4 r c (g·cm -3 ) Meas. Reflections/ (R int ) (0.0270) (0.0234) Indep. ref. [I > 2 s (I)] Parameters/restraints.
174 / 0
Hydrogen treatment
Refall
Refall
Goodness of fit
Final R indices [I > 2 s (I)]: R / wR R indices (all data): R / wR *The crystallographic details of the data collections at RT and 50 K using single crystal X-ray diffraction have been provided in the Supplementary Information (see text and Table S1). Table 1.
Experimental parameters and main structural crystallographic data for the studied compounds.
Table 2.
Main structural variations between high and low temperature phases obtained from single crystal neutron diffraction refinement.
Fe(1)-O(1W) 2.095(2) Fe(1)-O(1w) 2.095(2) Fe(1)-Cl(1) 2.3155(15) Fe(1)-Cl(1) 2.391(2) Fe(1)-Cl(2) 2.3926(7) Fe(1)-Cl(2) 2.400(2) Fe(1)-Cl(2) a Angles (°) Angles (°)
Cl(1)-Fe(1)-Cl(2) 94.40(1) Cl(1)-Fe(1)-Cl(2) 170.73(1) Cl(1)-Fe(1)-Cl(2) a a a a a a Symmetry codes:(a) x , 0.5- y , z The Fe(III) atom is located on a mirror plane perpendicular to the b -axis, and presents a slightly distorted octahedral environment. The shortest Fe-Cl bond distances correspond to Cl(1), which is coordinated in trans - conformation with respect to the water molecule (see Table 2). The interbond angles in the anion are slightly different of those of an ideal octahedron. These deviations are related with the strength of the different hydrogen bonds present in this compound. A clear example is the coordination water molecule hydrogen bond [O(1w)-D···Cl(2)], with a D···Cl distance of 2.217(2) Å. As result of this interaction the Fe-Cl(2) bond distances are slightly elongated. The influence of the H-bonds due to the ND counterion is more difficult to quantify due to the extension of the H-bond network in which the ND molecule is involved. A detailed list with all the possible H-bonds can be consulted in Table 3. The refinement of the ND molecule at high temperature was made by several successive attempts. Our first model considered the ND molecule ordered. An initial model was obtained using soft constraints, with the ND molecule geometrically fixed and only small deviations of this ideal geometry allowed. A similar protocol was also applied using rigid bodies in FullProf suite [21] but in that case the ND molecule forms an undeformable rigid body with a perfect tetrahedral symmetry. In order to define a realistic model, the N-D bond length was fixed at 1.025 Å and the D-N-D angle was fixed to 109°, values in agreement with other ammonium-containing compounds. [23] The values of the Fe–Cl, Fe-O and O-D bond lengths were left variable. Unfortunately all our attempts to refine this model produce unstable solutions. A more complex model considering the ND molecule disordered was subsequently tested. For the starting point of this model, the D atoms of the ND molecule were set in 8 different sites (Wyckoff position 8 d ) with an occupancy of 0.5 that is equivalent to two disordered ND units with half occupancy. The first set of refinements was made using the Simulated Annealing procedure included on the FullProf suite.[21] The refinement produced a solution where the 8 D atoms around the N can be seen as two superimposed tetrahedrons (see Figure 1). This model was the starting point for the latter refinement using the SHELX program.[14] The final refinement was made without restraints. The values of the N-D bond distances as well as the D-N-D angles are in agreement with those observed in other ammonium-containing compounds.[23] Figure 1.
Views of the 100K crystal packing along the c -axis (top) and b -axis (bottom) of (ND ) [FeCl (D O)] compound, where the O-D···Cl and the N-D···Cl hydrogen bonds have been represented as blue and green dashed lines, respectively. For the sake of clarity only one conformation of the disordered ND molecule has been drawn in polyhedron representation. Table 3.
Relevant hydrogen bonds for 100 K ( ) and 2 K ( ) measurements obtained from single crystal neutron diffraction refinement.
D-H···A * H···A / Å D-H···A / ° 100K O(1W)-D(1W)···Cl(2) b c d e d a f g h i j g f h k l * D and A stand for donor and acceptor, respectively. Symmetry codes: (a) x , 0.5- y , z ; (b) 2- x , -0.5+ y , - z ; (c) x , 0.5- y , 1+ z ; (d) 2- x , -0.5+ y , 1- z ; (e) 1.5- x , -0.5+ y , 0.5+ z ; (f) 1.5- x , - y , 0.5+ z ; (g) x , y , 1+ z ; (h) 2- x , 1- y , - z ; (i) 2- x , - y , - z ; (j) 2- x , - y , 1- z ; (k) x , y , -1+ z ; (l) 1.5- x , 1- y , -0.5+ z The deuterium atoms of the disordered ND act as donors in an extensive hydrogen-bond network where each deuterium is involved in a single H-bond, with distances ranging from 2.329(4) to 2.540(6) Å. The shortest H-bonds are those involving the coordinated water molecule and both trans -Cl(2) atoms, which connect the Iron(III) atoms forming a zig-zag chain along the b -axis. The shortest Fe···Fe distance connects iron atoms in the Fe(1)-O(1w)-D···Cl(2)-Fe(1) zig-zag chain running along the b -axis, with a distance between metallic centers of 6.457(1) Å, while the shortest interchain distance connects the Fe atoms within the ac -plane, with a distance of 6.8125(14) Å. Crystal structure determination at low temperature
There have been several previous attempts to refine the crystal structure of this compound at low temperature (below the structural phase transition);[11,12,24] nevertheless the reported structures above and below the phase transition were mainly the same, and neither changes in the space group nor a significant displacement on the atomic positions were observed. These observations have been also verified in this work with X-ray single crystal diffraction data collections above and below the nuclear phase transition. The best model fitting the X-ray data is always obtained with the structure refined on the space group
Pnma , independently of the temperature (in our case RT and 45 K, see Supplementary Information for details). These results suggest that the hydrogen atoms present in the structure are responsible for the phase transition, and therefore X-ray diffraction would not be the best probe to analyze it due to the low scattering power of X-rays by hydrogen. In order to verify our hypothesis, neutron diffraction measurements of the deuterated compound have been carried out on a powder sample (at 45 K) and on a single crystal (at 100 and 2K). Assuming that the possible phase transition is a second order (“continuous”) phase transition, we determined the list of possible subgroups of the high temperature space group (
Pnma ). The only possible space groups are
Pna (No. 33), Pmn (No. 31), Pmc (No. 26), P (No. 19), P / c (No. 14), P / m (No. 11), Pc (No. 7), Pm (No. 6), P (No. 4), P -1 (No. 2) and P (No. 1), where the polar ones are written in bold. Given that the structural transition (ca 79 K) is not associated with the onset of a ferroelectric state,[11] the only possible space groups immediately below this phase transition are in principle the non-polar ones. Once the system becomes ferroelectric below 6.9 K, these non-polar groups are no more strictly correct, but the structural changes responsible of this type of electric polarization are usually too subtle to be observable,[4] so we assumed the same space group at 45 K and 2 K. All possible space groups were tested by trial and error and the only subgroup that refines properly the experimental data was found to be P / c (No. 14). However, in order to facilitate comparisons between the high- (100K) and low-temperature phases the refinement was carried out in the P / a space group, which is a non-standard setting of P / c (No. 14). The reflection splitting due to the subtle monoclinic distortion can be clearly observed in the high resolution neutron powder diffraction data at 45K (see Figure S2). The results of the refinement of the single crystal data show a crystal structure at low temperature which is built up from one and two crystallographically independent [FeCl5·H2O] and ND units, respectively (see Figure 2). As previously mentioned, the low temperature phase crystalizes in P / a space group, which is a maximal non-isomorphic subgroup of Pnma with index 2. Therefore the number of crystallographically independent atoms at low temperature is doubled with respect to the high temperature phase. The phase transition involves the loss of a n -glide parallel to the a -axis and a mirror plane which is contained in the ac -plane, together with a slight variation in the cell parameters and the cell volume (see Table 1). The KPI packing index gives a percent of filled space of 77.2, a value which is slightly lower than that observed at high temperature.[22] However, as occurs in the Figure 2.
Views of the 2K crystal packing along the c -axis (left) and b -axis (right) of (ND ) [FeCl (D O)] compound, where the O-D···Cl and the N-D···Cl hydrogen bonds have been represented as blue and green dashed lines, respectively high-temperature phase, there is no accessible space in the network for solvents. At low temperature, the Fe(III) atoms are located on a general position due to the loss of the mirror plane perpendicular to the b -axis. The Fe(III) environment is slightly more distorted compared with the high temperature structure, but the octahedral environment is held. The largest variations are observed in the long bond distances involving Cl(1), Cl(2) and Cl(5), where a significant elongation of the Fe-Cl distance has been identified (see Table 2). The interbond angles in the anion are very similar to those observed at high temperature (see Table 2). The most remarkable difference between high- and low-temperature phases concerns the ND counterions which are well ordered below the phase transition (ca 79 K). An order-disorder transition is therefore proposed as the origin of this structural phase transition. The influence of the H-bond network in this transition has been taken into account by analyzing the bond-distances and the bond-orientation. At low-temperature all the bond-distances are slightly shorter than those observed at high-temperature. The bond-orientations are very similar, and in both cases the donor atoms interact with only one acceptor, so the structural transition can be seen as the splitting of the previously disordered ammonium into two independent positions, triggered by an reduction of the ND···Cl bond distances between the ND molecules and the [FeCl5·H2O] units (see Figure 3). A detailed list with all the possible H-bonds can be consulted in Table 3. The Fe···Fe distance along the zig-zag chain running along the b -axis is slightly longer than the observed at high temperature, with a value of 6.466(2) Å, while the interchain distance, which connects the Fe atoms in the ac -plane, is Fe···Fe 6.8061(13) Å, slightly longer compared with the high temperature one. Magnetic Structure
The magnetic structure of (ND ) [FeCl (D O)] was obtained from single-crystal neutron diffraction data taken at 2 K. is. A survey of reciprocal space was carried out with a series of Q -scans in the first Brillouin zone showing the appearance of superlattice reflections indexed by a propagation vector parallel to the c -axis, k = (0,0, k z ), with k z = 0.2288(4). This corresponds to a period of 30.2 Å in real space. In the following section, the magnetic structure is described by labeling the four Fe(III) atoms in the primitive unit-cell as Fe(1), Fe(2), Fe(3) and Fe(4) respectively with crystallographic coordinates (0.388, 0.249, 0.313), (0.119, 0.751, 0.813), (0.619, 0.751, 0.687) and (0.881, 0.249, 0.187). The magnetic moment [ m l ( j )] for atoms at positions Fe(j) (j=1 to 4) in the unit-cell l , can be calculated by means of the Fourier expansion: (cid:1) (cid:2) (j)= Re (Sj)cos{2π·[ k · R (cid:2) (cid:3) (cid:4)(cid:5)(cid:6)(cid:7)(cid:8)(cid:9) (cid:3) Im (Sj)sin{2π·[ k · R (cid:2) (cid:3) (cid:4)(cid:5)(cid:6)(cid:7)(cid:8)(cid:9) where R l , is the position vector of the unit-cell l with respect to the origin (cid:10) (cid:2) (cid:11) (cid:12) (cid:13) a (cid:3) (cid:12) (cid:14) b (cid:3) (cid:12) (cid:15) c , where l i are integers, j (j) are relative phases (in fractions of 2 π ) and Re (S j ) and Im (S j ) are the real and imaginary parts of the Fourier vectors for each site j . In order to determine the possible magnetic structures compatible with the crystallographic space group ( P / a ) and with the k = (0,0, k z ) propagation vector, we have used representational analysis as described initially by Bertaut.[25] The decomposition of the magnetic representation ( G) in the group of the wave-vector (little group) involves two one-dimensional irreducible representations ( G and G ). In this symmetry group, only Fe(1)/Fe(2) on one hand, and Fe(3)/Fe(4) on the other, are related by symmetry elements, so that the orbit splits into two sets. For each of them, G=3G +3G , generating three sets of basis vectors for each representation (Table 4). Magnetic structures compatible with a single representation (magnetic space group P a g )00 s ) have been tested for all combinations of the real and imaginary parts of the Fourier vectors S(j) (i.e. describing all collinear and non- Figure 3.
View of the hydrogen bond interactions between the [FeCl (D O)] unit and the adjacent molecules obtained from the neutron diffraction crystal structure at 100 (top) and 2 K (bottom). The O-D···Cl and the N-D···Cl hydrogen bonds have been represented as blue and green dashed lines, respectively. A detailed list of relevant hydrogen bond distances and angles can be consulted in Table 3. The asymmetric units are given in ORTEP representation with atoms at 50 % of probability. Ammonium hydrogen atoms in the HT phase, have been drawn on lime and orange colors in order to denote the two different conformations. collinear magnetic arrangements). In all cases, the refinements lead to poor agreement factors (over 25 % R F2 ) and therefore these solutions can be excluded. In particular, the amplitude-modulated models systematically overestimate the Fe(III) magnetic moment, with amplitude maxima of 5.23(3) m B , not physically meaningful. It deserves to be noted that none of these magnetic models would generate a macroscopic ferroelectric polarization, as observed experimentally. In order to fit successfully the experimental data, a combination of magnetic modes belonging to the G and G irreducible representations is required. The best solution was obtained with a cycloid model with moments mainly contained in the ac -plane and propagating along the c -direction (see Figure 4). This corresponds to the admixture in phase quadrature of a mode along the a -axis belonging to G and a mode along the c -axis belonging to G ; it should be noted that the v and q component of the G and G irreducible representations permit a small component along the b -axis. This solution has been confirmed by global optimization using a simulated annealing procedure.[21] In the final refinements, the real and imaginary parts of the Fourier components S(j) were constrained to be equal, i.e. (cid:16) (cid:17)(cid:18) = | Re(Sj) | = | Im(Sj) | , which corresponds to a spin cycloid with a circular envelope. The magnetic moments rotate within a -that includes the propagation vector and which forms an angle of 4.4(1)° with the ac -plane. A refinement using an elliptical envelope instead (Re(Sj) ≠ Im(Sj)), does not improve significantly the agreement factors: the R F2 reaches a value of 12.6 % and the difference between real and imaginary parts of the Fourier components is ca 0.3 m B . The best fit between the observed and calculated intensities was obtained with phase angles of zero for Fe(1) and Fe(3) and k z /2 for Fe(2) and Fe(4); therefore Fe(1) and Fe(3) as well as Fe(2) and Fe(4) are strictly antiparallel, while the angle between the magnetic moments of Fe(1)-Fe(4) and Fe(2)-Fe(3) is ca 41.5° (see Figure 4). A full list of structural and magnetic parameters is given in Table 5, as well as the result of the magnetic refinement at 2K (see also Figure S3). The refined value of the Fe(III) magnetic moment, 3.805(2) m B , is in agreement with those previously reported for the family of general formula A [FeX (H O)], where A stands for an alkali metal or ammonium ion and X for a halide ion.[12] This value, slightly below the expected for a Fe(III) ion, indicates a non-negligible spin delocalization from the Fe(III) ions to the coordinated chloride atoms.[26]
Table 4.
Magnetic moments of the content of a primitive cell deduced for the two possible irreducible representations ( G and G ) and for the combination of both, for the magnetic sites: Fe(1) = (0.388, 0.249, 0.313), Fe(2) = (0.119, 0.751, 0.813), Fe(3) = (0.619, 0.751, 0.687) and Fe(4) = (0.881, 0.249, 0.187). G G m Fe(1) (u,v,w) (u,v,w) m Fe(2) e (-u,-v,w) e (u,v,-w) m Fe(3) (p,q,r) (p,q,r) m Fe(4) -e (-p,-q,r) -e (p, q,-r) The symmetry operators relating the positions are: (1) = x , y , z ; (2) = - x +1/2,- y +1, z -1/2; e = exp{-i·2 p · 0.1136}, the common phase can be set to zero on the single representation, while for combined representations the phase difference between the basis vectors of the respective representations is π /2. The admixture of two magnetic modes belonging to different irreducible representations lowers the point symmetry to m ( ab -mirror plane) in perfect agreement with the emergence of an electrical polarization ( P ) in the ab -plane. The presence of cycloids rotating in the ac -plane (with a small tilt of 4.4° along b ) suggests that P should be preferentially directed along the a -axis, through the spin current mechanism induced via the inverse Dzyaloshinskii-Moriya interaction. The small tilt of the cycloid rotation plane towards b , is consistent with the weak electric polarization component along the b -axis previously observed.[11] Figure 4.
View along the b -axis of the superposition of the nuclear and magnetic structures for the combination of G G
2 irreducible representations. For the sake of clarity only [FeCl5(D2O)] units, have been represented. It should be noted that this combination of irreducible representations produces the same chirality for every spin chain running along the c -axis. Table 5.
Magnetic parameters of (NH ) [FeCl (H O)] obtained from the refinement of the magnetic reflections recorded at 2 K. a (Å) b (Å) 9.9578(5) c (Å) 6.9049(4) a (°) 90.00 b (°) 90.00 g (°) 90.109(4) k (0,0, k z ) k z M Fe ( m B ) 3.805(2) j (1)(2 p ) 0.0 j (2)(2 p ) a Anisotropic f Re (°) f Im (°) -92.89(6) R exp (%) 4.25 q Re (°) q I m (°) 175.6(1) R F (%) 12.8 f and q are the spherical angles of the Fourier components of the magnetic moment (real and imaginary). a Fixed values from the structural refinement at 2K.
The symmetry of the magnetic structure is consistent with either a first-order transition or two successive transitions, the latter being consistent with the presence of two anomalies observed in the low temperature region of the heat capacity (ca 6.9 and 7.2 K). This scenario is also in agreement with experiments conducted by Mösssbauer spectroscopy [27] and muon spin relaxation,[28] both pointing to the occurrence of two distinct magnetic structures. This sequence of transitions suggests that two order parameters directed respectively along the a - or c - crystallographic axis (and belonging to G and G ), condense in turn at T FE and T N . This corresponds to a transition on warming from cycloidal (T Temperature evolution of the (0,0,-1- k z ) and (0,-3, k z ) magnetic reflections, represented as blue and red circles, respectively. The change of curvature in the (0,0,-1.23) reflection could be associated to a magnetic phase transition from cycloid to collinear magnetic structure. The solid lines are guides to the eye. The proposed magnetic structure at zero magnetic field can be explained by simple considerations regarding the exchange integrals, discussed using the notation already employed in ref 16. All magnetic exchange couplings are mediated by super-super exchange interactions, either through Cl-Cl or Cl-water pathways, and bound to be antiferromagnetic in nature. From previous DFT calculation in the K [FeCl (H O)] compound,[29] J that couples Fe(1)-Fe(3) and Fe(2)-Fe(4) ions, must be the strongest interaction due to the short halogen-oxygen together with the enhancement of the magnetic interaction due to the well oriented hydrogen bond. Indeed our magnetic configuration shows that the spins of these pairs of ions are strictly antiparallel. For the sake of simplicity, we assume J =1 for the rest of the discussion. Figure 6. View along b - and c -axis of the four magnetic interactions using to model the magnetic ground state at 2K (see main text). The different exchange couplings have been represented in color dashed lines. For the sake of clarity only the Fe(III) ions have represented. The next most relevant topological feature is the presence of buckled anisotropic triangular layers lying in the bc -plane, involving the J and J interactions (see Figure 6). This minimal set of parameters ( J , J , J ) allows to obtain a magnetic ground state which reproduces the experimental wave-vector ( k z =0.23), since for classical Heisenberg spins, the analytic solution of the anisotropic triangular lattice ( J / J ) is a spiral structure with k z = 2·acos[ J /(2· J )]. In the present case, one can derive J » . It shall be noted that exact diagonalization of the first-ordered ground state shows that the minimum energy for any set { J , J , J , J } is either at the Gamma point or along the L symmetry line (0,0, k z ) in CDML notation.[30] The magnetic structure determined experimentally shows that the relative angles between spins on neighboring sites are slightly different from those obtained using { J , J , J } alone. A more realistic solution matching closely the observed angles requires the addition of a weak J coupling. By including J however, the J / J ratio must be renormalized in order to preserve the experimental wave-vector. For example, for J = 0.1, and J / J » J / J ratio compared to that reported for the K [FeCl (H O)] compound seems compatible with the shorter Cl(5)···O(1w) and Cl(3)···Cl(4) distances mediating the J coupling. The set of exchange couplings proposed here may not be a unique solution since some values were arbitrarily fixed, but their relative magnitude ( J > J > J > J ) are consistent with the observed structure. Ab-initio calculations and/or measurements of magnetic excitations are required to provide quantitative values of these exchange couplings. Conclusion With the ultimate goal of establishing the mechanism at the origin of multiferroicity in the title compound, we have carried out crystal and magnetic structure studies as a function of the temperature. In previous works, the transition identified at 79 K by heat capacity measurements was hypothetically attributed to a structural phase transition, but all efforts to determine the crystal structure below this temperature led to solutions indistinguishable from the high temperature structure. In this work, we have taken advantage of the contrast provided by neutron diffraction for the precise determination of the positions of light atoms (in this case mainly the hydrogen/deuterium atoms of the ammonium molecules) to characterize this structural phase transition, which is ultimately explained by a blocking of the ammonium counterions. The space group analysis at low temperature has shown that the title compound crystalizes in the P / a space group. The centrosymmetric character of the space group precludes ferroelectric polarization, as in fact is observed in the previously reported pyroelectric current measurements just below the structural phase transition.[11] The compound becomes magnetically ordered at 7.2 K (T N ) and becomes ferroelectric below 6.9 K (T FE ), with a spontaneous electric polarization presenting two different components: a strong one along the a -axis and a weak one along the b -axis, with a value ten times lower. Our magnetic neutron diffraction unambiguously shows that at 2K, an incommensurate cycloidal spin structure is stabilized at zero magnetic field. The magnetic moments are contained mostly in the ac -plane (ca 4° tilt along b ) and propagate along the c -direction. The symmetry lowering at 2K confirms that the ferroelectricity is a direct consequence of the complex magnetic ordering and suggests that the inverse Dzyaloshinskii-Moriya mechanism is at play in the magnetoelectric coupling. This is in sharp contrast with all other compounds of the series (A [FeX (H O)] , where A stands for an alkali metal and X is an halogen atom) presenting commensurate collinear antiferromagnetic structures and not ferroelectric. This highlights that subtle chemical modifications produce remarkable changes in the physical behavior of this class of compounds. Specifically, the substitution of the alkali-metal for an ammonium molecule modifies the exchange coupling interactions giving rise to a frustrated magnetic structure. We have provided a minimalistic model based on a set of four isotropic exchange interactions which explains well the observed magnetic state. Strong competition within anisotropic triangular layers seems to be responsible for the observed behavior. It is likely that external parameters such as magnetic field or pressure influence drastically the spin state and potentially lift the frustration, as recently suggested from magnetoelectric measurements under external magnetic field. Further neutron diffraction experiments under external magnetic field are needed in order to elucidate the different magnetic models and determine the mechanism of magnetoelectric coupling in the different regions of the phase diagram. Notes and references Acknowledgements Partial funding for this work is provided by the Ministerio Español de Ciencia e Innovación through projects MAT2010-16981, MAT2011-27233-C02-02. JARV acknowledges CSIC for a JAEdoc contract. We are grateful to the Institut Laue Langevin for the neutron beam-time allocated. The authors are especially grateful to Dr. C. Ritter (ILL) for the access to the D2B diffractometer for the high resolution low temperature neutron measurement. Corresponding Authors * OF [email protected] and LC [email protected] Supporting Information Single crystal X-ray crystallographic data at RT and 50 K, detail of the heat capacity curves of (ND ) [FeCl (D O)] compound, high-resolution neutron powder diffraction pattern at 45 K and Rietveld refinement using two different space groups, squared magnetic structure factors observed versus calculated for D9 instrument at 2K. Abbreviations , Compound High Temperature; , Compound Low Temperature; KPI is the acronym of Kitaigorodskii Packing Index; CDML is the acronym of Cracknell, A. P., Davies, B.L., Miller, S. C. and Love, W. F. References [1] Fiebig, M. J. Phys. , D 38 , R123. [2] Toledano P. Phys. Rev. B 2009, 79, 094416 [3] a) Goto, T.; Kimura, T.; Lawes, G.; Ramirez, A. P.; Tokura, Y. Phys. Rev. Lett. , , 257201. b) Van Aken, B.B.; Palstra, T.T.M.; Filippetti, A.; Spaldin, N.A. Nature Materials , , 164. c) Sergienko, I. A.; Dagotto, E. Phys. Rev. B , , 094434. d) Kimura, T.; Lawes, G.; Goto, T.; Tokura, Y.; Ramirez, A. P. Phys. Rev. B , , 224425. [4] a) Taniguchi, K.; Abe, N.; Takenobu, T.; Iwasa, Y.; Arima, T. Phys. Rev. Lett. , , 097203. b) Arkenbout, A. H.; Palstra, T. T. M.; Siegrist, T.; Kimura, T. Phys. Rev. B , , 184431. c) Heyer, O.; Hollmann, N.; Klassen, I.; Jodlauk, S.; Bohatý, L.; Becker, P.; Mydosh, J. A.; Lorenz, T.; Khomskii, D. J. Phys.: Condens. Matter , , L471. [5] a) Yamasaki, Y.; Miyasaka, S.; Kaneko, Y.; He, J.-P.; Arima, T.; Tokura, Y. Phys. Rev. Lett. , , 207204. b) Tomiyasu, K.; Fukunaga, J.; Suzuki, H. Phys. Rev. B , , 214434. [6] a) Dzyaloshinskii, I. J. Phys. Chem. Solids , , 241. b) Moriya, T. Phys. Rev. , , 91. c) Sergienko A.; Dagotto, E. Phys. Rev. B , , 094434. [7] Katsura, H.; Nagaosa, N.; Balatsky, A.V. Phys. Rev. Lett. , , 057205. [8] a) Harter, A; Cage, B; Nguyen, P; Abboud, K. A.; Dalal, N. S. Polyhedron , , 2350. b) Cage, B.; Dalal, N. S. Chem. Mater. , J. Am. Chem. Soc. , , 3792. b) Samantaray, R.; Clark, R. J.; Choi, E. S.; Dalal, N. S. J. Am. Chem. Soc. , CrystEngComm , , , 3558. b) Liu, B.; Shang, R.; Hu, K.-L.; Wang, Z.-M.; Gao S. Inorg. Chem. , , 13363. c) Xu, G.-C.; Ma, X.-M.; Zhang, L.; Wang, Z.-M.; Gao S. J. Am. Chem. Soc. , , 9588. d) Xu, G.-C.; Zhang, W.; Ma, X.-M.; Chen, Y.-H.; Zhang, L.; Cai, H.;-L.; Wang, Z.-M.; Xiong, R.-G.; Gao S. J. Am. Chem. Soc . , , 14948. e) Cañadillas-Delgado, L.; Fabelo, O.; Rodríguez-Velamazán, J. A.; Lemée-Cailleau, M.-H.; Mason, S. A.; Pardo, E.; Lloret, F.; Zhao, J.-P.; Bu, X.-H.; Simonet, V.; Colin, C. V.; Rodríguez-Carvajal J. J. Am. Chem. Soc. , , 19772. f) Jain, P.; Ramachandran, V.; Clark, R. J.; Zhou, H. D.; Toby, B. H.; Dalal, N. S.; Kroto, H. W.; Cheetham A. K. J. Am. Chem. Soc . , , 13625. [11] Ackermann, M.; Brüning, D.; Lorenz, T.; Becker, P.; Bohatý L. New Journal of Physics , , 123001. [12] Carlin, R. L.; Palacio, F. Coord. Chem. Rev. , , 141. [13] Agilent. CrysAlis PRO Agilent Technologies Ltd, Yarnton, England, . [14] Sheldrick, G. M. Acta Cryst. , A64 , 112. [15] Nardelli, M. J. Appl. Crystallogr . , , 659. [16] Spek, A. L. J. Appl. Cryst . , , 7. [17] DIAMOND 2.1d, Crystal Impact GbR, CRYSTAL IMPACT, K. Brandenburg & H. Putz GbR, Postfach 1251, D-53002 Bonn, Germany, . [18] Lehmann, M. S.; Larsen, F. K. Acta Cryst. , A30 , 580. [19] Coppens, P. Crystallographic Computing , Ed. F. R. Ahmed, Munksgaard International Booksellers and Publishers Ltd.: Copenhagen , 255. [20] Howard, J.A.K.; Johnson, O.; Schultz, A.J.; Stringer, A.M. J. Appl. Cryst. , , 120. [21] Rodríguez-Carvajal, J . Physica B , Molecular Crystals and Molecules , New-York, Academic Press, . [23] a) Rubín, J.; Palacios, E.; Bartolomé, J.; Rodríguez-Carvajal J. J. Phys.: Condens. Matter , Chem. Phys. , , 5161. [25] Bertaut, E. F.; in: Magnetism, Eds. Rado, G. T. and Shul H.; Academic: NewYork , Vol.III , Chap.4 [26] Luzón, J.; Campo, J.; Palacio, F.; McIntyre, G. J.; Millán, A; Phys. Rev. B , 78, 054414. [27] a) Partiti, C.S.M.; Rechenberg, H.R.; Sanchez, J.P. J. Phys. C: Solid State Phys. , , 5825. b) Calage, Y.; Dormann, J.L.; Moron M.C.; Palacio, F. Hyp. Interact. , , 483. c) Brown, S.R.; Hall, I. J. Phys.: Cond. Matter , , 9191.. [28] Attendoborough, M.; Hall, I.; Nikolov, O.; Brown, S.R.; Hyp. Interact. , , 435. [29] Campo, J.; Luzón, J.; Palacio, F.; McIntyre, G. J.; Millán, A.; Wildes, A. R. Phys. Rev. B , , 054415. [30] Cracknell, A. P.; Davies, B.L.; Miller, S. C.; and Love, W. F. ( ). Kronecker Product Tables. Vol. 1. General Introduction and Tables of Irreducible Representations of Space Groups. New York: IFI/Plenum. upporting information for: Magnetically-induced ferroelectricity in the (ND ) [FeCl ·D O] molecular compound. José Alberto Rodríguez-Velamazán, Óscar Fabelo, Ángel Millán, Javier Campo, Roger Johnson, Laurent Chapon. Instituto de Ciencia de Materiales de Aragón (ICMA), CSIC – Universidad de Zaragoza, 50009 Zaragoza, Spain. Institut Laue-Langevin, 38042 Grenoble Cedex 9, France. Univ Oxford, Dept Phys, Clarendon Lab, Oxford OX1 3PU, England Table S1. Experimental parameters and main structural crystallographic data for the studied compounds determined from single crystal X-ray diffraction. Formula [NH ] [Fe III Cl (H O)] Empirical Formula Cl H FeN Mr (g·mol-1) 287.20 Temperature (K) 293(2) 50(2) l (Å) 0.71073 0.71073 Crystal system Orthorombic Orthorombic Space group (No.) Pnma (62) Pnma (62) Crystal size (mm) 0.08 × 0.08 × 0.06 0.08 × 0.08 × 0.06 a (Å) 13.728(5) 13.5088(6) b (Å) 9.934(5) 9.9408(4) c (Å) 7.040(5) 6.9070(4) a (º) 90.00 90.00 b (º) 90.00 90.00 g (º) V (Å ) 960.0(9) 927.53(8) Z 4 4 r c (g·cm -3 ) 1.987 2.057 Meas. Reflections/ (Rint) 26820/(0.0442) 16496/(0.0931) Indep. ref. [I > 2 s (I)] 943 923 Parameters/ restraints. 69/ 4* 69 / 4* Hydrogen treatment Refall Refall Goodness of fit 1.139 1.132 Final R indices [I > 2 s (I)]: R / wR R / wR Figure S1. Details of the heat capacity curves of (ND ) [FeCl (D O)] compound in the low temperature range. The two different phase transitions observed (long range magnetic ordering at ca 7.4 K and ferroelectric ordering at ca 7.1K ) are in agreement with those previously characterized for the non-deuterated compound. Figure S2. Neutron diffraction pattern of (ND ) [FeCl (D O)] compound collected at 45K using D2B high resolution instrument with l = 1.5942 Å. The Rietveld refinements were done in the space group Pnma (up) and in P / a (bottom), the insets show the splitting due to the nuclear phase transition from orthorhombic to monoclinic system. The experimental data have represented as blue circles, the calculated curve as solid red line and the difference between them as solid green line. The green vertical lines represent the Bragg positions for the selected space group. Figure S3.