Switching of the Chiral Magnetic Domains in the Hybrid Molecular/Inorganic Multiferroic (ND4)2[FeCl5(D2O)].

(ND4)2[FeCl5(D2O)] represents a promising example of the hybrid molecular/inorganic approach to create materials with strong magneto-electric coupling. Neutron spherical polarimetry, which is directly sensitive to the absolute magnetic configuration and domain population, has been used in this work to unambiguously prove the multiferroicity of this material. We demonstrate that the application of an electric field upon cooling results in the stabilization of a single-cycloidal magnetic domain below 6.9 K, while poling in the opposite electric field direction produces the full population of the domain with opposite magnetic chirality. We prove the complete switchability of the magnetic domains at low temperature by the applied electric field, which constitutes a direct proof of the strong magnetoelectric coupling. Additionally, we refine the magnetic structure of the ordered ground state, deducing the underlying magnetic space group consistent with the direction of the ferroelectric polarization, and we provide evidence of a collinear amplitude-modulated state with magnetic moments along the a-axis in the temperature region between 6.9 and 7.2 K.

Introduction

Unlike “type I” multiferroics, where electric and magnetic orders coexist but are weakly coupled, the magneto-electric (ME) coupling is very strong in “type II” (or “improper”) multiferroics, where the magnetic ordering breaks spatial inversion symmetry and induces electric polarization[1-5]. This direct coupling leads to remarkable possibilities of manipulation of the electric order by magnetic fields and of the magnetic order by electric fields, which are at the origin of the significant attention devoted to this class of materials. The number of such spin-driven multiferroics is, however, still scarce, and mainly limited to transition metal oxides. Most of the examples involve non-collinear spiral spin arrangements, like in RMnO3 (R = rare earth)[6,7], MnWO4[8,9], CoCr2O4[10], CuFeO2[11] and RbFe(MoO4)2[12], although collinear or almost collinear magnetic structures can also give rise to spin-driven multiferroics, like in HoMnO3[13] or in the RMn2O5 family[14]. The difficult control of the delicate balance between exchange couplings in these systems impedes a rational design of these interesting materials. This fact highlights the interest for innovative approaches for obtaining magnetically driven multiferroic compounds. Important efforts have been made in this sense following metal-organic routes, which allow the combination of building blocks with the desired properties in a bottom-up approach[15-18]. However, the cross coupling between magnetic and electric orders is absent or very weak in this type of systems[19,20]. As an alternative, hybrid approaches inspired by the architectures of the successful organic-inorganic photovoltaic materials have been proposed[21]. For this research, microscopic techniques that give direct proof of the multiferroicity and give insight on the underlying magneto-electric coupling mechanism are crucial, since a good understanding of these mechanisms may help to design new improper multiferroic materials with improved functional properties.

The family A2[FeX5(H2O)][22], where A stands for an alkali metal or ammonium ion and X for a halide ion, represents a new route to obtain materials with strong ME coupling. The occurrence of ME coupling in the alkali-based compounds K2[FeCl5(H2O)], Rb2[FeCl5(H2O)] and Cs2[FeCl5(H2O)] has been described in a thorough study of their macroscopic physical properties[23], and understood from the analysis of their crystal and magnetic structures[24,25]. While the previous compounds are linear magneto-electric materials where electric polarization is only observed when a magnetic field is applied (no spontaneous electric polarization is observed in the ground state), another member of the family, (NH4)2[FeCl5(H2O)][21], has been proven to be a veritable spin-driven multiferroic, i.e. the appearance of an electric polarization is spontaneous and directly related to the magnetic ordering in the ground state. We recently determined the mechanisms of multiferroicity in the deuterated form of this material in different regions of its rich magnetic field-temperature (B-T) phase diagram from a detailed determination of its crystal and magnetic structures by neutron diffraction[26,27]. In absence of magnetic field, this compounds becomes magnetically ordered at 7.2 K (TN) and ferroelectric below 6.9 K (TFE). The magnetic structure in the ground state is cycloidal spiral propagating along the c-axis and with magnetic moments mainly contained in the ac-plane. A ferroelectric polarization, primarily directed along the a-axis, develops in this phase, which is compatible with the spin current mechanism[26]. In a scenario reminiscent of that observed in canonical compounds like TbMnO3 or MnWO4[6-9], the symmetry of the magnetic structure in the ground state is consistent with two successive transitions (as indicated by two anomalies observed in the low temperature region of the heat capacity), suggesting that two order parameters condense in turn at TFE and TN. This would correspond to a transition on warming from a cycloidal (T < TFE) to a collinear antiferromagnetic structure (TFE < T < TN), the latter being compatible with the absence of electric polarization.

In the present work, we prove the possibility of controlling the “chiral” magnetic domains of this material by an applied electric field (E-field). We use neutron spherical polarimetry to determine the absolute magnetic configuration and domain population of this system under different external E-fields. The cycloidal magnetic structure of the ground state implies two “chiral” domains with opposite rotation of the cycloids (Fig. 1), which directly correspond to electric domains with opposite polarities. Therefore, by applying an E-field, we can manipulate the population of the “chiral” magnetic domains, and we show that the system is completely switchable, which represents a direct evidence of its multiferroic character. In addition, we use neutron spherical polarimetry data to refine the previously proposed magnetic structure of the ground state, obtaining a more accurate magnetic structure model. Finally, we confirm the existence of a collinear sinusoidal magnetic phase with magnetic moments along the a-axis in the temperature region between 6.9 and 7.2 K[28].

Figure 1

Schematic representation of the magnetic structure of the two magnetic domains, with the direction of the resulting electric polarization.

Results

Spherical neutron polarimetry considerations

Spherical neutron polarimetry[29] is an ideally suited technique to study complex magnetic structures. This is a zero magnetic field technique which allows to access the direction and phase of the magnetic scattering, unreachable by unpolarized methods. In fact, neutron polarization analysis[30], and in particular neutron spherical polarimetry[31,32], have demonstrated to be particularly adapted for studies very similar to the one presented here. The experimental setup allows setting the polarization of the incident neutron beam in an arbitrary direction and analyzing the polarization of the scattered beam in another spatial direction[33]. In the present case, with a magnetic propagation vector k = (0, 0, k) [k = 0.23], the magnetic and nuclear Bragg peaks never superimpose, and therefore the polarization of the scattered beam for each Bragg peak can be calculated, for any incident polarization, merely from the magnetic structure factor M(Q) (a complex vector)[34,35].. The most convenient setup for these calculations is to use the Blume reference frame[34], with the X-axis parallel to the scattering vector Q, the Z-axis perpendicular to the scattering plane (vertical in our case), and the Y-axis completing the right-hand set. Magnetic neutron diffraction is only sensitive to the component of M(Q) perpendicular to Q, M⊥(Q), the so-called magnetic interaction vector, therefore lying in the YZ plane: M⊥(Q) = (0, M⊥y, M⊥z). The scattering of a fully polarized incident neutron beam can be expressed by the polarization matrix P, with matrix elements Pij (i, j = X, Y, Z), representing the polarization of the scattered beam in the direction j, for an incident beam polarized in the direction i (see Supporting Information for details). In the present case, for a pure magnetic Bragg peak probed at Q (for simplicity, we will drop the argument Q in the expressions and consider a perfectly polarized neutron beam):

The off-diagonal components, Pyx, Pzx, are the so-called ‘chiral’ terms, and can only be non-zero for non-collinear structures. Their value not only depends on the magnetic structure, but also on the populations of the possible magnetic domains, hence the particular sensitivity of spherical neutron polarimetry to the domain structure. In particular, in the case of cycloidal structures as the one of the ground state of the present compound, the rotation of the cycloid is opposite for each possible “chiral” domain, and thus the signs of the ‘chiral’ terms are reversed for the same magnetic reflection probed at Q. For a single domain, the ‘chiral’ terms will have a finite value, with the sign indicating the handedness (clockwise or counter-clockwise) of the “chiral” magnetic domain. Finally, the chirality of a single domain is only observable when the projection of the cycloidal envelope along Q is non-zero. The sense of rotation of the cycloid is coupled with the direction of the ferroelectric polarization through the spin current mechanism[26], and therefore the application of an external E-field can modify the populations of the described magnetic domains. If both domains are equally populated, the ‘chiral’ terms average out and P = P = 0 (as would be the case for a collinear structure), but if the domain populations are unbalanced, which can be realized by cooling in an applied E-field in the direction of the ferroelectric polarization, finite values are measured, which allow the determination of the relative fraction of cycloidal domains.

For collinear arrangements, like the proposed at TFE < T < TN, off-diagonal terms are zero, but the diagonal terms contain information about the orientation of the magnetic moments that is often undistinguishable by unpolarized neutron diffraction techniques. In the general case, the final values of both the diagonal and off-diagonal terms for each magnetic reflection further depend on the details of the magnetic structure and the exact orientation of the moments with respect to Q. In helical structures, the diagonal Pyy and Pzz components contain information about the ellipticity (for example, for a reflection with the magnetic moment nearly perpendicular to Q, these terms are close to zero when the envelope is circular), while the off-diagonal terms Pyz and Pzy are sensitive to the inclination of the plane of rotation of the magnetic moments. Therefore, the measurement of polarization matrices in a series of magnetic reflections allows for a particularly precise refinement of the magnetic structure (especially if combined with a set of magnetic Bragg peaks intensities obtained by unpolarized neutron diffraction), together with the determination of the relative fraction of spiral domains.

Collinear phase

The spherical neutron polarimetry results in the TFE < T < TN phase are displayed in Fig. 2-Top – right (measurements at T = 7 K after zero-E-field cooling). As expected, the polarization matrix elements of the (0 1 k) reflection have non-zero values only for the diagonal terms, Pii, consistent with the magnetic structure being collinear, and thereby sinusoidal, since it remains incommensurate[26]. The external electric field does not affect the neutron polarization matrix in this non-ferroelectric phase (the only domains in this phase are the indistinguishable anti-phase domains, corresponding with a π/2 shift on the sinusoidal structure, and not coupled to the electric polarization). The symmetry lowering from the paramagnetic phase to the ground-state cycloidal phase indicates that two order parameters condense successively at TFE and TN. Depending on what is the first one that condenses, two sinusoidal structures are possible in the TFE 28]. Figure 2-Bottom displays the temperature dependence of selected terms of the polarization matrix, showing how the value of the off-diagonal Pyx term becomes zero when entering in the 6.9 K < T < 7.2 K phase, as corresponds to its collinear character, while the diagonal terms evolve to absolute values close to one, with signs consistent with a collinear sinusoidal magnetic structure with the magnetic moments directed along a.

Figure 2

Neutron spherical polarization matrix elements for the (0 1 0.23) magnetic reflection of (ND4)2[FeCl5(D2O)]. (Top) Observed (solid color columns and error bars) and calculated (black rectangles) neutron spherical polarization matrix elements, Pij, for the (0 1 0.23) magnetic reflection of a (ND4)2[FeCl5(D2O)] crystal oriented with a* parallel to the Z-axis. Left: Results corresponding to the cycloidal phase (T = 4 K) for the sample cooled in zero electric field (orange) and under a negative (purple) and positive (red) electric field of 25 kV cm−1 applied along a*. The calculated values are the results of the joint fit of neutron spherical polarimetry and unpolarized neutron diffraction intensities. Right: Experimental results (orange bars) corresponding to the sinusoidal phase (T = 7 K) for the sample cooled at zero-field, compared with the calculated values for a sinusoidal magnetic structure with moments parallel to c and a directions (green and black rectangles, respectively). Inset: Scheme of the experimental geometry showing the crystal (red/green/blue arrows) and Blume (gray arrows) reference frames, together with the cycloid rotation plane with one particular magnetic moment (both in magenta).(Bottom) Temperature dependence of observed neutron spherical polarization matrix elements, Pij, for the (0 1 0.23) magnetic reflection in zero-E-field cooling and negative E-field cooling conditions.

Cycloidal phase

We examined in detail the magnetic structure and the domain population in the cycloidal phase at T < TFE by measuring polarization matrices for a series of magnetic reflections in both zero-E-field cooled and E-field cooled conditions (positive and negative fields along the a* direction). The results corresponding to measurements of the (0 1 0.23) magnetic reflection at 4 K are shown in Fig. 2-Top – left, while Fig. 2-Bottom shows the temperature dependence of selected polarization matrix terms. Further results for other magnetic reflections, entirely consistent with the ones in Fig. 2, are provided in the Supplemental Material. The results of the full polarization matrix measurement for the single (0 1 0.23) magnetic reflection immediately evidence:

(i) The finite values of the ‘chiral’ components, Pyx (and Pzx), confirming the non-collinear character of the magnetic structure in this phase and the E-field induced unbalance of domain populations.

(ii) The change of sign of these terms when the sample is cooled in either negative or positive E-field (while the rest of the terms remain the same), which demonstrates the reversal of the domain populations by a reversal of the E-field.

A joint fit of data of three full polarization matrices and the intensities of 127 independent magnetic reflections obtained by unpolarized neutron diffraction (see Figures S2 and S3 in Supplemental Material) yields a refined model of the magnetic structure in the ground state consisting on a cycloid with an ellipticity of 0.97(2), and an angle between the rotation plane of the magnetic moments and the ac-plane of 2.0(2)° (the full set of parameters of the refinement is given in Table 1). This model slightly improves the fit with respect to the previously reported in ref.[26]. The four Fe(III) atoms in the primitive unit-cell, labeled as Fe(1), Fe(2), Fe(3) and Fe(4), 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) respectively, were grouped in two orbits where the magnetic moments were constrained to be equal and the envelopes of the spin cycloid to be circular, in order to avoid over parametrization of the fit. In the present case, the last constraint is lifted, allowing an elliptical envelope. Concerning the symmetry of the system, once it becomes ferroelectric below 6.9 K, the non-polar P1121/a space group is not strictly correct, although the structural changes responsible of this type of electric polarization are usually too subtle to be observable. However, the magnetic symmetry should account for the polar character of the system. The magnetic space group should be a polar subgroup of the parent space group P1121/a1′. The highest polar symmetry subgroup is P11211′(00γ)000s, but this magnetic space group belongs to the 21′ point group, which only allows ferroelectric polarization along the c-axis and consequently is not compatible with the spontaneous electric polarization observed from the pyroelectric measurements[21]. Therefore, the only magnetic super-space group deriving from the parent group compatible with the observed polarization is P11′(αβγ)0s, which allows ferroelectric polarization in whatever direction, and should correspond to the symmetry of our system. In this way, the breaking of symmetry from the parent space group explains the physical properties of the system. The fitting of the domain population gives values of 0.03/0.97(2) for negative E-field cooling and 0.97/0.03(3) for positive E-field cooling, that is, a nearly full population of one or the other domain depending on the sign of the applied E-field. A small magnetic domain unbalance of 0.39/0.61(2) is observed when the sample is cooled without applied E-field, with values of Pyx and Pzx different from zero.

Table 1

Refined Fourier components of the magnetic moment, S1 = 1/2(R + iI) and S2 = 1/2(R + iI)exp(−2πiϕ) for the two independent Fe atoms (Ri and Ii are the components along the crystallographic axes) from the joint fit of neutron spherical polarimetry and unpolarized neutron diffraction intensities.

Rx	Ry	Rz
−3.96(6)	0.14(8)	0
Ix	Iy	Iz
−0.01(1)	−0.36(8)	−3.83(6)

Agreement factors: R = 12.52% for the integrated intensities (slightly better than in ref.[26]) and χ2 = 7.20 for the polarization matrix elements. Both independent magnetic atoms are constrained to have equal values of the components of the magnetic moment, with a refined phase difference of ϕ = 0.495(5) between both.

Switching of magnetic domains

Once established the coupling between the magnetic domain population and the E-field applied on cooling, we demonstrate how we can switch the magnetic domains by E-field reversal at constant temperature, which constitutes an unambiguous and direct proof of the multiferroicity on this material. Starting from a negative E-field-cooled state, we have followed one of the chiral terms, Pyx, of the polarization matrix of the (0 1 0.23) magnetic reflection–directly related with the magnetic domain population–as a function of a variable E-field (Fig. 3). The magnetic domains are strongly pinned at the lowest temperatures, with the polarization matrix term P, changing by only ca. 25% for applied positive E-fields up to 25 kV cm−1 (corresponding to a domain unbalance of ca. 0.10/0.90). Therefore, we performed the switching measurements at 6.5 K, sufficiently close to the ferroelectric transition. We can clearly observe a complete hysteresis loop, with the population of the magnetic domains fully and reversibly switched by the applied E-field, that is, we have direct control over the magnetic domain population by E-field.

Figure 3

Hysteresis loop measured on the off-diagonal neutron polarization matrix element Pyx of the (0 1 0.23) reflection of a (ND4)2[FeCl5(D2O)] crystal. The sample was initially cooled down to 4 K under a negative E-field of 17.5 kV cm−1 and then warmed to 6.5 K, temperature that was kept constant for the subsequent measurements as a function of a variable E-field. Inset: Domain populations (purple and red bars refer to negative and positive chiral magnetic domains, respectively) corresponding to the increasing E-field branch of the hysteresis loop.

In summary, by means of spherical neutron polarimetry experiments, combined with unpolarized neutron diffraction data, we have determined in detail the magnetic structures of (ND4)2[FeCl5(D2O)] hybrid compound in both its cycloidal (T < 6.9 K) and sinusoidal (6.9 K < T < 7.2 K) phases. In the temperature region between 6.9 and 7.2 K we have determined a collinear sinusoidal magnetic structure with magnetic moments along the a-axis, consistent with the symmetry lowering scheme. At temperatures below 6.9 K, we have obtained an improved model for the magnetic structure of the ground state, determining the underlying magnetic space group, which allows deducing a symmetry for the system compatible with the observed ferroelectric polarization. We are able to stabilize a nearly full population of either chiral magnetic domain by cooling in an external electric field. Furthermore, we can tune the domain populations by varying the electric field, reaching a full population reversal by switching the field in a complete hysteresis loop, thus providing a direct proof at the microscopic level of the multiferroicity of (ND4)2[FeCl5(D2O)], which represents a promising example of the hybrid molecular/inorganic approach to materials with strong magneto-electric coupling.

Methods

Sample preparation

A large prism-shaped single crystal of (ND4)2[FeCl5(D2O)] (P1121/a, with a ≈ 13.50 Å, b ≈ 9.96 Å and c ≈ 6.90 Å and γ ≈ 90.1° in the low temperature phase) of dimensions ca. 2 × 3 × 4 mm along the crystallographic a-, b- and c- directions was obtained by the seeded growth technique as described in our previous work[26].

Spherical neutron polarimetry

Experiments were conducted on the hot neutron diffractometer D3 at the Institut Laue-Langevin (ILL, Grenoble, France), equipped with the CryoPAD device[33], operating at λ = 0.84 Å. The sample was mounted onto a sample stick allowing applying up to 10 kV high voltage by a potential difference between two parallel horizontal aluminum plates. The crystal was fixed to the lower plate by silver epoxy with the a*-axis in the vertical direction (see Figure S1), and the upper plate positioned at ca. 1 mm from the sample surface (total distance between electrodes, 3.15 mm, sample thickness, ca. 2 mm). The aluminum sample chamber was indium-sealed and evacuated, and installed into the CryoPAD device at D3 instrument. Measurements were carried out after either zero-field cooling or electric-field cooling through the multiferroic transition. The crystal alignment with an accuracy better than 0.5 degrees was obtained by a previous orientation of the crystal using the neutron Laue diffractometer Orient Express[36]. We have probed several magnetic reflections, (0 1 k), (0 1 − k) and (0 0 1 − k), selected considering their magnetic signal and geometrical conditions to have the propagation vector in the scattering plane and a significant component of the magnetic moment perpendicular to Q (see inset of Fig. 2 top).

Unpolarized neutron diffraction

Data were collected at the hot-neutrons four-circle diffractometer D9 at ILL (λ = 0.832 Å). The crystal was placed into a closed-cycle cryostat and 127 independent magnetic reflections were recorded at 2 K (together with a large set of nuclear reflections that was used to accurately set the scale factor).

Data analysis

We used the Mag2Pol program[37] for analyzing and fitting the spherical neutron polarimetry data jointly with unpolarized neutron diffraction intensities. The calculated polarization matrix components were corrected for the imperfect incident neutron beam polarization and the efficiency of the neutron spin-filter in the scattered beam. The spin-filter efficiency was monitored by measuring the polarization term Pzz on the purely nuclear Bragg peak (0 2 0).

Data availability

The data supporting the findings of this study are available from the corresponding authors on request. Neutron scattering data obtained at Institut Laue-Langevin are available at: https://doi.ill.fr/10.5291/ILL-DATA. 5-51-503.

Supporting Information