Record high-Tc and large practical utilization level of electric polarization in metal-free molecular antiferroelectric solid solutions.

Metal-free antiferroelectric materials are holding a promise for energy storage application, owing to their unique merits of wearability, environmental friendliness, and structure tunability. Despite receiving great interests, metal-free antiferroelectrics are quite limited and it is a challenge to acquire new soft antiferroelectric candidates. Here, we have successfully exploited binary CMBrxI1-x and CMBrxCl1-x solid solution as single crystals (0 ≤ x ≤ 1, where CM is cyclohexylmethylammonium). A molecule-level modification can effectively enhance Curie temperature. Emphatically, the binary CM-chloride salt shows the highest antiferroelectric-to-paraelectric Curie temperature of ~453 K among the known molecular antiferroelectrics. Its characteristic double electrical hysteresis loops provide a large electric polarization up to ~11.4 μC/cm2, which endows notable energy storage behaviors. To our best knowledge, this work provides an effective solid-solution methodology to the targeted design of new metal-free antiferroelectric candidates toward biocompatible energy storage devices.

Introduction

Electroactive materials that enable to store and release electrical energy have been widely used as basic ingredients of energy-storage systems, electromagnetic devices, and high-power microwaves[1-4]. Among them, antiferroelectric (AFE) materials, featured by the antipolar array of electric dipoles inside the neighboring sublattices[5,6], have long been a fascinating branch for energy storage. Particularly, such initially antiparallel dipoles can be switched into parallel alignment under high electric fields, namely, AFE-to-ferroelectric (FE) structural transformation[7,8]. Compared with linear dielectrics and ferroelectrics, this distinctive electric-induced phase transition in AFEs leads to characteristic double P–E hysteresis loops, thereby allowing to store a larger amount of electric energy[9]. The current mainstream of AFE materials for energy storage is still dominated by inorganic oxides containing heavy metals[10,11], such as PbZrO3, (Ba, Na)TiO3, and AgNbO3, of which antiferroelectricity stems from ionic displacement and their optimized performances highly depend on chemical doping to these archetypes. To work toward a sustainable and green society, nevertheless, metal-free molecular AFEs are booming as an alternative family. For example, purely organic squaric acid (SQA) shows a comparable storage efficiency with inorganic oxides[12], and a giant AFE-related electrostriction was observed in 2-trifluoromethylnaphthimidazole[13]. However, the majority of existing molecular AFEs suffer from a quite low Curie temperature (Tc), such as NH4H2PO4 (Tc = 147 K)[14], 5,5′-dimethyl-2,2′-bipyridinium chloranilate (Tc = 318 K), and SQA (Tc = 373 K)[7], indicating the low operating temperature condition (<100 °C) that is not suitable for practical applications[15-17]. Therefore, it is of significance and urgent to exploit new metal-free molecular AFE candidates with high Tc and large spontaneous polarization (Ps).

Solid solution methodology has been substantially devoted to optimizing and modifying the physical properties of ferroic-ordered materials, including both inorganic oxides and inorganic-organic hybrid ferroelectrics[18-20]. For instance, AgNbO3 crystals exhibit quite weak room-temperature ferroelectricity, while the metal-doping of La3+ and Sm3+ results in strong AFE orders in AgNbO3-based ceramics that facilitate high-efficiency energy storage. Recently, solid solutions of inorganic-organic hybrid ferroelectrics, (TMFM)(TMCM)1–CdCl3 (0 ≤ x ≤ 1), were successfully achieved through cation modifications, which display a huge enhancement of piezoelectric coefficients up to 1540 pC/N near the morphotropic phase boundaries. For purely metal-free systems, dynamic orderings of molecular moieties usually dominate the potential energy barriers that are crucial to modulating Tc for phase transitions[21-23]. Nevertheless, the formation of metal-free solid solutions is quite challenging, since the mixing of different molecules in a single homogeneous phase closely involve with their spatial orientation and configurational distortion. This probably accounts for the lacking of high-Tc molecular AFEs achieved by the solid solution method. In this context, we here attempt to construct new metal-free AFEs through the solid solution methodology, of which the countered anions serve as variables to form the multi-composition solid solutions.

Here, we have successfully obtained new metal-free molecular AFEs in a series of binary CMBrI1− and CMBrCl1− solid solutions (0 ≤ x ≤ 1, where CM is cyclohexylmethylammonium) as single-phase homogeneous crystals. A subtle modification of halogen-bonded dimers raises phase transition energy barriers that result in the enhanced Tc. Among them, the CM-chloride crystal (CMC, x = 0 of CMBrCl1−) exhibits an AFE-to-paraelectric (PE) phase transition at ~453 K, falling in the range of the record-high values among molecular AFEs. Particularly, a large practical utilization level of Ps ~ 11.4 μC/cm2 in CMC, as deduced from the double P–E hysteresis loops, reveals its potentials for energy storage application. This is the first time to achieve such a record high-Tc and large practical utilization level of Ps in the molecular AFE systems. This result sheds light on exploiting metal-free AFE candidates for future biocompatible energy storage.

Results and discussion

Molecule-level design of metal-free AFEs

The family of CMBrI1− and CMBrCl1− solid solutions adopt a hydrogen-bonding dimer motif, resembling the binary stator–rotator system in AFE-to-PE phase transition (Fig. 1a). Organic CM+ cation is the dynamically active moiety with order–disorder genius like “rotator”, while halogen anions linked to cations via N-H···X hydrogen bonds stabilize the structure. The cooperativity of dipolar components in the adjacent crystal lattices afford the power source to create long-range AFE orders, mainly determined by the cationic reorientation. A hint is the stator–rotator interactions play a crucial role to increase energy barriers, as well as the realization of high-Tc AFE orders. Notably, this binary solid-solution family has a delicate compositional dependence on intermolecular H-bonds, as shown by Hirshfeld surface, which results in the enhanced Tc. Taking CMB as a prototype[24], the height of energy barriers can be successfully regulated by I/Br/Cl halogen modification (Fig. 1b, c), which leads to a wide tunable rang of Tc from 324 K of CMI to 453 K of CMC.

Fig. 1

Diagram for the molecule-level design of metal-free molecular AFEs in a series of CMBrxI1− and CMBrxCl1− (0 ≤ x ≤ 1) solid solutions.

a Basic structures for solid solutions. b The relationship between Curie temperatures and D···A distances of hydrogen bonds, where donor (D) is N of organic CM+ cation and acceptor (A) is halogen (X). The D···A distances can be obtained from the crystal structures; the averaged electronegativity of the solid solution is estimated from the ratio x based on the known Cl/Br/I. c Tc for the solid solutions is dramatically enhanced through the gradual substitution, the background refers to Hirshfeld surface to quantify intermolecular interactions.

Temperature-dependent crystal structure transition

Colorless crystals of CMBrI1− and CMBrCl1− (0 ≤ x ≤ 1) were synthesized from aqueous solutions with the stoichiometric ratios (Supplementary Fig. 1). Results of X-ray photoelectron spectroscopy and Powder X-ray diffraction of Cl/Br/I elements coincide with the theoretical calculations (Supplementary Figs. 2, 3). Structure analyses reveal that both CMC and CMI, isostructural to our reported CMB analog, crystallize in orthogonal system with the point group mmm at room temperature, thus enabling the compatibility of solid solutions. Figure 2a shows the basic unit consisting of H-bonding dimers, resembling the famous stator–rotator systems[25-29], in which the orientations of H-bonding dimers are arranged in the opposite directions (Fig. 2b), which eliminate bulk polarization of unit cell and demonstrate possible AFE order. Moreover, the N-H···X distances between stator and rotator are dependent on the halogen sizes. The N···Cl distance of CMC ≈ 3.148 Å is much shorter than N···Br ≈ 3.300 Å and N···I ≈ 3.529 Å. Such different interactions related to ion radius and electronegativity cause a significant change of Tc. Take CMC as an example, CM+ cations become highly disordered at high-temperature paraelectric phase accompanied by a halving of the volume (Fig. 2c–f). Dynamic ordering and reorientation of C-N bonds and six-member rings of two adjacent bipartite arrays require much more energy to overcome barriers, corresponding to the higher Tc.

Fig. 2

Crystal structures of binary solid solutions.

a, b Schematic diagram of CMBrI1− and CMBrCl1− binary solid solutions. The organic cations afford a driving force to trigger phase transitions. c Projection of AFE structure of CMC at LTP. d The cell edges of both the AFEP and PEP were drawn to show volume multiplication of AFE. e, f Structure packing of CMC at the paraelectric HTP P4/nmm, showing the complete disordering of organic CM cations.

The CMBrI1− and CMBrCl1− solid solutions have an isostructural high-temperature phase (HTP) with the nonpolar space group P4/nmm (point group 4/mmm), in which organic CM+ cations are rotationally disordered and adopt 4-fold rotation axis (C4) symmetry (Supplementary Figs. 4, 5). The C-N bonds of cations occupy eight equivalent sites that obey the basic requirement of crystallographic symmetry. Structure comparison reveals that the order–disordering accounts for their phase transitions and the vanishing of bulk polarization. The halogen anions of these solid solutions exhibit relative atomic displacements in FEP, thus coinciding with the variation of crystal symmetry[24]. Hence, the origin of AFE-FE-PE phase transitions can ascribe to the collaboration between organic cationic ordering and anionic displacement.

Phase transitions and related properties

To verify the occurrence of structural phase transitions, we analyzed differential scanning calorimetry (DSC) on the typical compounds in the CMBrI1− and CMBrCl1− solid solutions (Fig. 3a and Supplementary Fig. 6). Sequences of sharp endo/exothermic peaks in the heating/cooling modes indicate that their Tc values can be efficiently regulated between 324 K (CMI) and 453 K (CMC), which depend on the compositions of anionic halogen. It is intriguing to note that the CMB and CMBr0.8Cl0.2 solid solutions undergo two phase transitions, of which the intermediate phase is determined as a polar space group of P4mm. Furthermore, strong λ-peak anomalies are clearly observed in the temperature-dependent dielectric measurement (Fig. 3b). With the composition changing from I to Br to Cl, the Tc for observing sharp dielectric anomalies is also increasing, in agreement with DSC and structure analyses (Fig. 3c). Especially, the monocomposition CMC exhibits both the highest Tc and strongest dielectric response, including low dielectric loss (<0.01, at room temperature) and large dielectric constant (εʹ) variation (Fig. 3d). This should be reminiscent of its significant AFE order. Subsequently, we further determine the phase structures as a function of temperature and the halogen composition. As depicted in Fig. 3e and Supplementary Fig. 7, CMC has the record-high Tc among the metal-free molecular AFEs at 453 K, e.g. 5,5′-dimethyl-2,2′-bipyridinium chloranilate (~318 K), SQA (~373 K), and catches up with some oxide counterparts, such as Pb0.97La0.02(Zr0.50Sn0.43Ti0.07)O3 (~433 K) and Pb0.93Ba0.04La0.02Zr0.95Ti0.05O3 (~ 450 K)[12,30,31]. As far as we know, this is the first time to report metal-free molecular solid solutions with such a widely adjustable Tc, and CMC has the record-high Tc which is well below the melting point (Fig. 3f), which probably knocks the door toward practical application.

Fig. 3

CM content-induced phase transition in solid solutions.

a DSC curves of CMBrI1− and CMBrCl1− solid solutions. b Temperature variation of εʹ upon cooling. c Schematic diagram of temperature dependence versus cell lattice parameters obtained directly from the X-ray single-crystal diffraction, as depicted by the short dash lines. d Temperature dependence of εʹ and dielectric loss measured at different frequencies for CMC. e Phase diagram and the point groups of CMBrI1− and CMBrCl1− solid solutions. The dotted boundary is roughly estimated due to the difficulty of growing homogeneous crystals. f TG/DTA curves of CMC with heating rate 15 °C/min.

Antiferroelectric and energy storage behaviors

Generally, the typical P–E hysteresis loop is the most convincing evidence to reveal bulk AFE and FE orders, which are microscopically manifested as dipoles flipping under external electric field. Therefore, we prepared electrodes by observing the optical axis and ferroelastic domains under a polarized light microscope to observe the characteristic P–E hysteresis loops from the current density versus electric field (J–E) curves for CMBrI1− and CMBrCl1− solid solutions (Fig. 4a–c and Supplementary Fig. 9). Figure 4d depicts the selected members show obvious current peaks in the J–E curves at different temperatures. To emphasize, an additional current peak emerges in CMBr0.2Cl0.8, suggesting the coexistence of both AFE and FE orders that resemble many other inorganic solid solutions[32,33]. For CMC, two pairs of current peaks can also be clearly observed just below Tc, which act as the direct proof of AFE properties. From the hysteresis loops in Fig. 4e, the AFE Ps of CMC is estimated as 11.4 μC/cm2, much higher than other molecular AFEs, such as 2-trifluoromethylbenzimidazole (~5.6 μC/cm2), (5,5′-dimethyl-2,2′-bipyridine) (chloranilic acid) (~4.9 μC/cm2), and hybrid perovskite (isopentylammonium)2CsPb2Br7 (~6 μC/cm2)[12,34,35]. Besides, the coercive electric field (Ec) that defines the strength to reverse dipolar moments displays temperature-dependent behavior; for instance, the forward switching field (EA, AFE-to-FEP transition) increases from 41 kV/cm at 443 K to 55 kV/cm at 433 K (Supplementary Fig. 8), and even higher Ec is require to switch electric dipoles. Such large practical utilization level of Ps and high Ec provide great potential for CMC toward energy storage applications.

Fig. 4

AFE properties of CMBrI1− and CMBrCl1− solid solutions.

a Variable-temperature domains of CMC. b, c Preparation of electrodes according to the relationship between optical axis and the direction of electric polarization. d J–E curves collected at different temperatures. e Variable-temperature double P–E hysteresis loops of CMC.

Above studies reveal that CMC is a metal-free molecular AFE with the record high-Tc and large practical utilization level of Ps, which are the essential criteria for energy storage application. Figure 5a presents the temperature range for retaining antiferroelectricity in several molecule-based AFE materials, of which CMC has a quite wide working temperature range. Furthermore, energy storage characters of CMC, including total energy density (Wst), recoverable energy density (Wre) and storage efficiency (ƞ), are obtained from the integration of P–E curves (Fig. 5b). The calculated Wre and Wst are 0.28 and 0.54 J/cm3, respectively; the corresponding ƞ reaches up to 52% at 433 K, which is beyond that of CMBr0.8Cl0.2 (~20%, at 387 K). This efficiency is also comparable to some AFE oxides, e.g. AgNbO3 (~40%), Hf0.3Zr0.7O2 (~51%), Pb0.97La0.02(Zr0.50Sn0.39Ti0.11)O3 (~50%), 0.95NBT-0.05BH (~48.2%), etc[17,36-38]. Given structural tunability and environment-friendly merits, these metal-free AFE solid solutions hold great application potentials in flexible and wearable devices.

Fig. 5

AFE and energy storage related characteristics.

a AFEP Curie temperature (Tc) of CMC compared with some reported molecule-based AFE compounds. The columnar with twill background present the temperature range for experimentally measuring the antiferroelectric properties (for details, see Table S9). b Schematic illustration of the calculation of Wre (green) and Wloss (yellow) from the double P–E hysteresis loop.

To summarize, we have presented a family of metal-free molecular AFEs in CMBrI1− and CMBrCl1− (0 ≤ x ≤ 1) solid solutions by a precise molecule-level modification of halogen compositions. The CM-chloride salt displays the highest Tc of molecular AFEs (~453 K) and large practical utilization level of Ps up to 11.4 μC/cm2. Such AFE attributes are unprecedented for the molecular AFEs, which allow notable energy storage performances in a relatively broad temperature range. This work provides guidance on the design of new metal-free AFE candidates, and might knock the door toward biocompatible, environmentally-friendly and flexible device applications.

Methods

Synthesis of CMX (X = Cl, Br, I)

The series of binary CMBrI1− and CMBrCl1− (0 ≤ x ≤ 1) solid solutions as single crystals were obtained via the temperature-cooling method. Cyclohexanemethylamine (1.5 ml, Aladdin, 99.9%) was dissolved in HX solution (20 ml, Aladdin, 48 wt% in water) by heating under strong stirring to colorless solution. Subsequently, the colorless sheet crystals were obtained by the temperature-cooling technology with a speed of ~1 K/day and dried for subsequent use.

Synthesis of CMBrCl1− and CMBrI1− (0 < x < 1)

Stoichiometric CMB (x × 5 mmol) and CMC/CMI ((1−x) × 5 mmol) were added into gamma-Butyrolactone/water (10 ml) to get a clear solution. After slow cooling at a rate of 1 K/day, the colorless sheet crystals can be obtained. The distribution of C, N, Cl, Br and I in the whole selected area reveals that the elements are uniformly dispersed on the single crystals. (Supplementary Figs. 10–12).

Single crystal structure determination

Single crystal X-ray diffraction data of CMBrCl1− and CMBrI1− were collected on a Bruker D8 using the Mo Kα radiation at different temperatures. The direct method solved all the crystal structures and then refined them by the full-matrix least-squares refinements on F2 using SHELXLTL software package. The non-hydrogen atoms were refined anisotropically based on all reflections with I > 2σ (I).

Thermal analysis

Differential scanning calorimetry (DSC) measurements were carried out using the NETZSCH DSC 200 F3 with a heating and cooling rate of 20 K/min under the nitrogen atmosphere. Thermogravimetric analysis (TGA) was measured using a STA449C Thermal Analyser ranging from room temperature to 900 °C with a heating rate of 15 °C/min.

Dielectric and ferroelectric measurements

The dielectric constants (εʹ) and dielectric loss were measured using the two-probe AC impedance method with an Impedance Analyzer (TH2828A). Moreover, the variable temperature current density versus electric field (J–E) curves and the polarization versus electric field (P–E) hysteresis loops were recorded with a ferroelectric analyzer (Radiant Precision Premier II).

Calculation of energy storage properties

Total energy density Wst, recoverable energy density Wre and energy density efficiency ƞ can be estimated by the following equations:ƞ = Wst/Wre, and Wloss = Wst−Wre, where E is the applied electric field, Pmax and Pr represent the maximum and remnant polarization, respectively. These energy storage characteristics can be obtained from the integration of P–E curves.