Ferroelectric Single-Molecule Magnet with Toroidal Magnetic Moments.

Materials that coexist magnetic and electric properties on the molecular scale in single-molecule magnets (SMMs) with peculiar quantum behaviors have promise in molecular electronics and spintronics. Nevertheless, such molecular materials are limited in potentials because their magnetic signal cannot be transformed into an electrical signal through magnetoresistance or Hall effects for their high insulativity. The discovery of an entirely new material, ferroelectric SMMs (FE SMMs) is reported. This FE SMM also shows single-molecule magnetic behaviors, toroidal magnetic moments, and room-temperature ferroelectricity. The toroidal moment is formed by a vortex distribution of magnetic dipoles in triangular Dy3 clusters. The analysis of ac magnetic susceptibility reveals the coexistence of three distinct magnetic relaxation processes at low temperatures. The ferroelectricity is introduced by incorporating polar alcohol molecules in the structure, which is confirmed by the X-ray diffraction and optical second harmonic generation (SHG) measurements. Moreover, the dielectric measurements reveal a ferroelectric-to-ferroelectric phase transition around 150 K due to the symmetry change from Pc to Pna21 . The coexistence of toroidal moment and ferroelectricity along with quantum magnetism in the rare-earth single-molecule magnets yields a unique class of multiferroics.

Introduction

Single‐molecule magnets (SMMs) are discrete molecular species with quantum tunneling behaviors of magnetism.[ ] Many SMMs with different structures and magnetic ions have been designed and synthesized in the past two decades.[ , , , ] Recently, SMMs with fantastic toroidal moments have attracted much attention owing to their special spin arrangement and abundant quantum magnetism,[ ] where the toroidal moment is characterized by a vortex distribution of magnetic dipoles.[ ] The toroidal moment is regarded as an electromagnetic moment in addition to traditional electric polarization and magnetization,[ ] and the concept of multiferroics has been extended to include the ferro‐toroidal order.[ ] A lot of intriguing properties based on toroidal moments have been observed in recent years.[ , , , , ]

Although the magnetic properties of SMMs have been extensively studied, their electric properties have been less concerned because they are generally insulating dielectrics without magnetoresistance or Hall effects. This feature makes SMMs unfavorable for applications in electronic and spintronic devices.[ , ] However, when ferroelectricity is introduced into SMMs, their functionality would be greatly expanded in terms of multiferroicity and magnetoelectric effects. In fact, there has been growing interest in achieving ferroelectricity in organic molecular materials and metal‐organic frameworks beyond traditional ferroelectric oxides.[ , , ] Nevertheless, ferroelectricity in SMMs has not been convincingly achieved. In this letter, we report the discovery of ferroelectricity above room temperature in a peculiar SMM with toroidal magnetic moments.

In this work, we constructed and characterized the solid evidences of SMM and ferroelectricity in a ferroelectric SMM,[ ] [Dy3(HL)(H2L)(NO3)4]·C2H5OH (H4L = N,N,N’,N’‐tetrakis(2‐hydroxyethyl)‐ethylene‐diamine), which posesses toroidal moment formed by a vortex distribution of magnetic dipoles in triangular Dy3 clusters. Interestingly, we find three‐step magnetic relaxation processes coexisting in the same temperature range. A modified Debye model with three relaxation times was introduced to study the multiple dynamic mechanisms as well. The ferroelectricity with the Curie temperatures T C ≈ 470 K was evidenced by the SHG signal that is only allowed in non‐centrosymmetric structures. Moreover, the dielectric anomaly around 160–170 K suggests a ferroelectric‐to‐ferroelectric phase transition originated from the structural transformation from Pc to Pna21.

Results and Discussion

X‐ray diffraction studies reveal that the SMM crystallizes in the polar space group Pna21 at 298 K (Figure ). There are three crystallography‐independent Dy3+ ions and one alcohol molecule in the polar Pna21 space group. One eight‐coordinated Dy3+ ion (Dy1) is coordinated by eight oxygen atoms (DyO8) from two nitrates and two ligands, while the other two nine‐coordinated Dy3+ ions (Dy2 and Dy3) are surrounded by two nitrogen atoms and seven oxygen atoms (DyN2O7) provided by one nitrate and two ligands. At 30 K, the space group changes to another polar space group Pc. The main difference between the crystal structures of Pna21 and Pc is the Z‐value and the ordering of solvent alcohol molecule. As a result, there are two types of molecular structure in the polar Pc space group, as shown in Figure 1b. Four different orientation triangular molecules and four ethanol molecules surround the 21 screw axis in the Pna21 phase but two different orientation triangular molecules and two ethanol molecules in the Pc phase which induced different molecular packing. The local coordination symmetry of each Dy3+ center was analyzed by continuous shape measure [ ] calculations on the DyO8 and DyN2O7 sites. At both 298 and 30 K, the calculations reveal that Dy1 is situated in a slightly distorted D 4d geometry while Dy2 and Dy3 are situated within a C 4v coordination environment with different deviation values.

Figure 1

Crystalline structure and toroidal magnetic moments. a) Crystal structure and space group at 298 and 30 K. b) Toroidal magnetic moments in two configurations of Dy3 clusters at low temperatures. The arrows represent the calculated orientations of the local magnetic easy axes on Dy3+ ions. Color scheme: green Dy, red O, grey C, and blue N. The hydrogen atoms and the isopropanol molecules are omitted for clarity.

Based on the determined geometries of the Dy3 complex at 30 K, Complete‐active‐space self‐consistent field (CASSCF) calculations on individual Dy3+ fragments were carried out. The energy levels, tensors, and the predominant m values of the lowest eight Kramers doublets (KDs) of individual Dy3+ fragments were obtained. POLY_ANISO program[ , , , ] was used to fit the Dy3+‐Dy3+ exchange interactions and the intermolecular interactions through comparison of the computed and measured magnetic susceptibilities. The computational details are described in the Supporting Information. The results suggest that the Dy3+‐Dy3+ interactions in the Dy3 clusters are antiferromagnetic. The magnetic axes on Dy3+ ions in the two molecular structures are indicated in Figure 1b. The included angles between magnetic easy axes of Dy3+ ions are all larger than 97o. Therefore, the local magnetization vectors of Dy3+ ions can be approximately regarded as a triangle. Subsequently, a toroidal moment is formed based on the Dy3 magnetic triangles.

A vanishing small magnetic susceptibility at low temperatures would be expected for the state with toroidal magnetic moments, as reported in previous studies.[ ] Figure shows the temperature dependence of dc magnetization of the Dy3 complex. The magnetization exhibits a maximum of around 4.5 K and decreases toward zero at low temperatures. This feature is consistent with the anticipated behavior of a SMM with toroidal magnetic moments.

Figure 2

Magnetic relaxation analyses. a) Temperature dependence of magnetization measured in 1000 Oe. b–d) Imaginary component of ac magnetic susceptibility as a function of frequency in different temperature ranges. e) The Cole‐Cole plots at 6–12 K show three distinct magnetic relaxation processes, named faster, medium, and slower, respectively. f) The relaxation time ln(τ) as a function of the reciprocal of temperature. The different relaxation processes are marked in different symbols. The solid lines are the best fitting results using the sum of three modified Debye functions.

The ac magnetic susceptibility of the Dy3 complex exhibits a strongly frequency‐dependent behavior in the temperature range from 2 to 20 K (see Figure S1, Supporting Information), which is another characteristic of SMMs. Especially, the imaginary component (χ“) shows an unusual feature of multiple peaks, indicating multiple magnetic relaxation processes. In the low‐temperature region (2–4.5 K), χ” shows two separate peaks, corresponding to a slow (low frequency) and fast (high frequency) relaxation process, respectively, as shown in Figure 2b. While the slow relaxation shifts to higher frequencies with increasing temperature, the fast relaxation process is almost temperature‐independent, which suggests that the fast relaxation is related to quantum tunneling of magnetization (QTM).[ ] In the intermediate temperature region (6–12 K), three distinguishable relaxation processes, named as the faster, medium, and slower processes, coexist (Figure 2c). Both the slower and medium relaxations shift to higher frequencies as temperature increases. In the high‐temperature region, the medium relaxation process gradually shifts to high frequencies and becomes undetectable while the slow relaxation process remains (Figure 2d). The strong shift with the temperature of the slower and medium relaxation processes indicates that they are related to thermally activated processes.

Interestingly, the experimental data between 6 and 12 K can be nicely simulated and depicted as the Cole–Cole plots[ ] which provide a more indicative view of the evolution and coexistence of multiple relaxation processes, as shown in Figure 2e. The relaxation time (τ) for each process can be extracted by fitting the χ"(ν) curves to three relaxation processes using the sum of three modified Debye functions,[ ] where χ S = χ S1 + χ S2 + χ S3 represents the sum of the adiabatic susceptibilities of the three and two relaxing species; Δχ is the difference between the adiabatic susceptibility (χ S,) and the isothermal susceptibility (χ T,) of each magnetic phase; ω = 2πν and the parameter α is the quantifying the width of the τ distribution. The small α values (α < 0.30) indicate that each relaxation phase has a very narrow distribution of relaxation times. And the remainder two temperature ranges of 2 to 4.5 K and 13 to 20 K are fitted by the sum of two modified Debye functions with χ S = χ S1 + χ S2.

The correlation between the relaxation time (τ) and temperature (T) can be obtained from a plot of ln(τ) versus ln(T) to give n values for the three relaxation processes (n is the parameter that represents the relation between the relaxation time and temperature in the equation τ = T –n). These results indicate that the medium and slow processes are dominated by the Raman process and Orbach relaxation pathways, while QTM dictates the faster process (Figure 2f and Figure S2, Supporting Information).

The corresponding magnetization blocking barriers of individual Dy3+ fragments are shown in Figure , where the transversal magnetic moments in the ground KDs of individual Dy3+ fragments are all smaller than 10–2 µ B, and thus the quantum tunneling of magnetization (QTM) in their ground KDs could be suppressed at low temperature. For individual Dy3+ fragments of and , the transversal magnetic moments in their first excited KDs are all smaller than 10–1 µ B, and thus the relaxation can probably proceed through higher excited KDs. For individual Dy3+ fragments of and , the transversal magnetic moments in the second excited KDs are all larger than 10–1 µ B, therefore, allowing a fast QTM in their second excited KDs. Hence, the calculated energy barriers of individual Dy3+ fragments of them according to the scheme of Figure 3 are 450.0, 316.8, 507.0, 437.1, 326.2, and 512.2 cm–1, respectively. Due to the unfavorable effects of anharmonic phonons, Raman magnetic relaxation, QTM, etc on the energy barrier, the experimental effective U eff are much smaller than the calculated energy gaps between the lowest three KDs of individual Dy3+ fragments.

Figure 3

Magnetization blocking barriers of individual Dy3+ fragments. The thick black lines represent the KDs as a function of their magnetic moment along the magnetic axis. The green lines correspond to diagonal quantum tunneling of magnetization (QTM); the blue line represents off‐diagonal relaxation process. The numbers at each arrow stand for the mean absolute value of the corresponding matrix element of transition magnetic moment.

The presence of multiple relaxation channels in SMMs is not unusual, but they usually play a dominating role in different temperature regimes. The coexistence of three relaxation processes at the same temperature has never been observed. Chibotaru et al. concluded that reorientation of the toroidal magnetization requires consecutive transitions through three exchanging excited KDs of Dy3 with the nonmagnetic ground state, whose excitation energies represent the barrier of blockage of this magnetization.[ ] Three exchange excited doublets for each compound are close to each other in Table S7, Supporting Information. Thus, we deduced the coexistence of three relaxation processes at the same temperature may arise from the three near‐degenerate exchange excited KDs.

The polar space groups at both high and low temperatures imply that this toroidal SMM could be ferroelectric. We employed the SHG technique to testify intrinsic ferroelectricity in the sample (Figure ).[ , , ] The relationship between the intensity of the SHG signal from a noncentrosymmetric crystal and its nonlinear polarization is I (2ɷ)/|P (2ɷ)|2. As shown in Figure 4a, the intensity of the reflection SHG signal of the sample presents a nonlinear variation with incident fundamental optical power, which can be well defined by the quadratic equation, confirming the second nonlinear response of the nanocrystals. In addition, the SHG signals obtained at high temperatures are almost zero. In contrast, a nonzero SHG signal is visible below ≈470 K, which corroborates the emergence of electric polarization and confirms intrinsic ferroelectricity with a Curie temperature above room temperature.

Figure 4

Second nonlinear response. a) Quadratic dependence of the SHG signal as a function of the laser power. As the incident light power increases, the SHG signal increases following the fitting curve y = A*x2 with a high R 2 (Adjusted R Square) value, confirming the second nonlinear response of the nanocrystals. Inset is the Schematic diagram of SHG measurement. The red line is the fitting curve to quadratic dependence. The inset shows the scheme of the SHG measurement. b) SHG intensity as a function of temperature. A ferroelectric phase transition is evidenced by a rapid increase of the intensity around 470 K.

Figure shows the dielectric permittivity (ε r) and loss (tanδ) as a function of temperature between 80 and 300 K. The sudden rise of ε r as well as the clear peak in tanδ indicate that there is another phase transition around 150 K. Moreover, the SHG intensity abruptly decreases to a minimum but nonzero value of around 150 K (see Figure S6, Supporting Information). The sharp decline in the SHG signal is likely due to the space group change from Pna21 to Pc. Since both space groups are polar, this transition around 150 K corresponds to a ferroelectric–ferroelectric phase transition.

Figure 5

Temperature dependence of a, b) the dielectric permittivity and c,d) loss tangent of the SMM. The dielectric anomaly indicates a ferroelectric‐to‐ferroelectric phase transition around 150 K at 1 Hz.

Conclusion

In summary, our study verifies a unique Dy3 SMM with toroidal magnetic moments and ferroelectricity, where the local magnetization vectors of Dy3+ ions form a triangle and the ferroelectricity is introduced by polar alcohol molecules. Interestingly, three distinct magnetic relaxation processes, corresponding to the Raman, Orbach, and QTM relaxation pathways, respectively, can coexist in the same temperature range. Hence, a modified Debye model with three relaxation times was introduced to study the multiple dynamic mechanisms. The ferroelectricity with the Curie temperatures T C ≈ 470 K was evidenced by the SHG signal that is only allowed in noncentrosymmetric structures. Moreover, the dielectric anomaly around 150 K (1 Hz) suggests a ferroelectric‐to‐ferroelectric phase transition originated from the symmetry change from Pc to Pna21. A significantly enhanced ME coupling would be expected in polar SMMs. This could pave the way for the design of molecular multiferroics and magnetoelectric materials using ferroelectric SMMs in the near future. The successful incorporation of ferroelectricity into SMMs not only broadens the family of multiferroics but also opens more opportunities for the electrical control of quantum magnetism of SMMs.

Experimental Section

The SMM in this study was a trinuclear (Dy3) complex with a formula of [Dy3(HL)(H2L)(NO3)4]·C2H5OH. The samples were synthesized by a solution method.[ ] Dy(NO3)3·6H2O (0.5 mmol, 228.3 mg) was slowly added to a 20 mL isopropanol solution containing H4L (1.0 mmol, 273 mg) and LiOH (0.4 mmol, 16.8 mg). The mixture was moved to a vial after 2 h of stirring. Colorless crystals were obtained under 100 ˚C solvothermal condition three days later in a 46% yield (based on Dy). Single‐crystal X‐ray diffraction was performed on an Agilent SuperNova diffractometer (298 K) and a Bruker D8 Venture diffractometer (30 K) equipped with graphite monochromated Mo‐Kα radiation (λ = 0.71073 Å). The structures were solved by direct methods and refined by the full‐matrix least‐squares method on F2 with anisotropic thermal parameters for all nonhydrogen atoms using the SHELXL program. Hydrogen atoms were located geometrically and refined isotropically (see Table S1, Supporting Information). Powder X‐ray diffraction patterns were measured on a Rigaku Ultima IV diffractometer using Cu‐Kα radiation (shown in Figure S7, Supporting Information). Magnetic properties were measured on powder samples using a Quantum Design superconducting quantum interference device magnetometer. The dielectric permittivity was measured on a palate made of pressed powder samples with an Aglient 4980A LCR meter in a cryogen‐free superconducting magnet system (Oxford Instruments, Teslatron PT). The optical second harmonic generation (SHG) measurements were performed by using a mode‐lock femtosecond (fs) Ti:sapphire oscillator (Tsunami 3941‐X1BB, Spectra‐Physics), which generates the incident fundamental laser with a central wavelength at 800 nm, a pulse duration of ≈100 fs, and a repetition rate of 82 MHz. The energy of the incident light was attenuated to 40 mW before being focused. A 400 nm filter was used in the reflective light path to ensure that only SH photons arrived at the photomultiplier tube. A single‐photon counting technique was conducted to count the SH photons, indicating the intensity of SHG signals generated from the samples. CASSCF calculations on individual Dy3+ fragments of trinuclear complex based on X‐ray determined geometries at 30 K have been carried out with MOLCAS 8.4 and SINGLE_ANISO programs.

Conflict of Interest

The authors declare no conflict of interest.

Supporting Information

Click here for additional data file.