# Structure and thermal property relationships in the thermomaterial di-n-butylammonium tetrafluoroborate for multipurpose cooling and cold-storage

Javier García-Ben, Juan Manuel Bermúdez-García, Richard J. C. Dixey, Ignacio Delgado-Ferreiro, Antonio Luis Llamas-Saiz, Jorge López-Beceiro, Ramón Artiaga, Alberto García-Fernández, Ute Cappel, Bruno Alonso, et al.

► To cite this version:

Javier García-Ben, Juan Manuel Bermúdez-García, Richard J. C. Dixey, Ignacio Delgado-Ferreiro, Antonio Luis Llamas-Saiz, et al. Structure and thermal property relationships in the thermomaterial di-n-butylammonium tetrafluoroborate for multipurpose cooling and cold-storage. Journal of Materials Chemistry A, 2023, 11 (41), pp.22232-22247. 10.1039/d3ta04063a. hal-04258310

HAL Id: hal-04258310

https://hal.science/hal-04258310v1

Submitted on 14 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Journal of Materials Chemistry A

Materials for energy and sustainability

Accepted Manuscript

This article can be cited before page numbers have been issued, to do this please use: J. Garcia Ben, J. M. Bermúdez-García, R. J. C. Dixey, I. Delgado-Ferreiro, A. L. Llamas-Saiz, J. López-Beceiro, R. Artiaga, A. Garcia-Fernandez, U. B. Cappel, B. Alonso, S. Castro-Garcia, A. E. Phillips, M. Sánchez-Andújar and M. A. Señarís-Rodríguez, *J. Mater. Chem. A*, 2023, DOI: 10.1039/D3TA04063A.

This is an Accepted Manuscript, which has been through the Royal Society of Chemistry peer review process and has been accepted for publication.

Accepted Manuscripts are published online shortly after acceptance, before technical editing, formatting and proof reading. Using this free service, authors can make their results available to the community, in citable form, before we publish the edited article. We will replace this Accepted Manuscript with the edited and formatted Advance Article as soon as it is available.

You can find more information about Accepted Manuscripts in the **Information for Authors**.

Please note that technical editing may introduce minor changes to the text and/or graphics, which may alter content. The journal's standard **Terms & Conditions** and the **Ethical guidelines** still apply. In no event shall the Royal Society of Chemistry be held responsible for any errors or omissions in this Accepted Manuscript or any consequences arising from the use of any information it contains.

<figure><img src="image_2.png" /></figure>
View Article Online
DOI: 10.1039/D3TA04063A

ARTICLE

# Structure and thermal properties relationships in the thermomaterial di-n-butylammonium tetrafluoroborate for multipurpose cooling and cold-storage

Javier García-Bena,‡, Juan Manuel Bermúdez-Garcíaa,‡, Richard J. C. Dixeyb, Ignacio Delgado-Ferreiroa, Antonio Luis Llamas-Saizc, Jorge López-Beceirod, Ramón Artiagad, Alberto García-Fernándezc, Ute B. Cappelc, Bruno Alonsof, Socorro Castro-Garcíaa, Anthony E. Phillipsb, Manuel Sánchez-Andújara* and María Antonia Señarís-Rodrígueza**

Nowadays around 46% of food production around the world requires refrigeration, which is generally provided either by active vapour-compression (based on refrigerants with liquid-gas transitions) or passive cold-storage (based on solid-to-liquid phase change materials, SL-PCMs). However, in order to avoid fluid losses during the transitions, new thermomaterials with solid-solid transitions are desired for both applications. In this work, we find that [DBA][BF₄] (DBA = di-n-butylammonium) is a promising thermomaterial with solid-solid phase transitions. This compound present thermal properties of great interest not only for active barocaloric refrigeration, but also for passive cold-storage, which render this as a unique multipurpose thermomaterial. The observed cold-storage capacity is very close to commercial SL-PCMs (E ~ 135 kJ kg⁻¹), while the pressure-induced thermal changes (ΔS ~ [200-270] J K⁻¹ kg⁻¹) are superior to most barocaloric materials and operating under lower pressures (p ~ [500-1000] bar). Moreover, the operating temperature range of this material is very adequate for food preservation (250-310 K), which is a great advantage over most barocalorics. Beyond the thermal properties, we perform a deep structural characterization, which reveal a progressive structural disorder of the [DBA]⁺ cations and [BF₄]⁻ anions as the origin of such thermal properties, which will help towards future rational design of enhanced thermomaterials.

Received 00th January 20xx,
Accepted 00th January 20xx

DOI: 10.1039/x0xx00000x

## 1. Introduction

There is a clear interest in finding new alternatives for our cooling systems to increase their energy efficiency and decrease the carbon emissions associated with this sector. This is especially important since around 20% of the global electric energy consumption is already devoted to refrigeration,¹ both at the industrial and domestic levels, with the aim of preserving food, medicines, achieve body comfort temperature, cooling electronics, and so on.² Unfortunately, this extensive use of refrigeration is also responsible for around 8% of global greenhouse gas emissions.³

Moreover, energy consumption in the refrigeration sector is expected to triple by the year 2050 due, among other things, to global warming. This would mean that, for the first time in history, more energy will be employed for cooling than for heating.

Food preservation is a particularly important aspect for maintaining the sustainable development and growth, since currently 46% of the total food production in the world requires refrigeration. Nevertheless, only half of it is properly refrigerated, which leads to around 13% of food losses every year.¹ Most perishable food must be preferably stored at chilled temperatures of between 272 and 287 K (-1 °C and 14 °C).4,5

To reach and/or maintain those temperatures two main technologies are normally used, namely, vapour-compression refrigeration and passive cold-storage technologies.

The first type of technology, for example commercial and domestic fridges, employs refrigerant gases in pressurization/depressurization cycles, where refrigerant gases can decrease the temperature from ambient down to sub-zero Celsius. However, those gases present risks for the user and/or the environment, since they are flammable, toxic and/or greenhouse gases. Nowadays, the Kigali Amendment to the Montreal Protocol and the F-gas regulation (EU 517/2014) is phasing out up-to 80% of fluorinated refrigeration gases (the vast majority of commercial refrigerant gases) by the year 2050. Meanwhile, hydrocarbons and other flammable alternatives are strictly controlled by the ATEX directive (2014/34/EU). In addition, gases could leak, ending up in the atmosphere and needing replacement, while solids are easier to handle.

<footnote>a. University of A Coruña, QUIMOLMAT Group, Dpt. Chemistry, Faculty of Science and Centro Interdisciplinar de Química e Biologia (CICA), Zapateira, 15071 A Coruña, Spain.
E-mail: m.andujar@udc.es, m.señaris.rodriguez@udc.es.</footnote>

<footnote>b. School of Physical and Chemical Sciences, Queen Mary University of London, London, E1 4NS UK.</footnote>

<footnote>c. Research Infrastructures Area, X-ray Unit, University of Santiago de Compostela, 15782 Santiago de Compostela, Spain.</footnote>

<footnote>d. Universidade da Coruña, Campus Industrial de Ferrol, CITENI-Grupo Proterm, Campus de Esteiro, 15403 Ferrol, Spain.</footnote>

<footnote>e. Division of Applied Physical Chemistry, Department of Chemistry, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden.</footnote>

<footnote>f. ICGM, CNRS, Université de Montpellier, ENSCM, 34095 Montpellier Cedex 5, France.</footnote>

<footnote>‡ These authors have equally contributed to this work.</footnote>

<footnote>Electronic Supplementary Information (ESI) available: [CCDC 2213758 and CCDC 2213757 contains the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/cif]. See DOI: 10.1039/x0xx00000x</footnote>
ARTICLE

Cold-storage technologies use thermally-induced phase change materials (PCMs) that can absorb heat from the surroundings and store it, mainly in the form of latent heat, so they can help to cool down when the ambient temperature increases above their phase transition.⁶ Most PCMs used for cold-storage applications are solid materials that undergo solid-to-liquid phase transitions at transition temperatures (from 20 °C down to sub-zero Celsius) depending on the cooling application.⁶ These include eutectic mixtures (mixture of salts and salt/water), hydrated salts, paraffins, fatty acids and clathrates.⁶⁻⁹ These substances exhibit very large thermal changes, although they also present important drawbacks for devices, such as: leaks of the liquid phase, incongruent melting, or mechanical strain due to large volume changes, among others.⁶

Trying to solve some of the limitations of refrigerant gases and the solid-liquid cold-storage materials, fluids-free alternatives are gaining an increasing interest. That is the case, for example of the classic thermoelectric refrigeration based on semiconductors¹⁰,¹¹ or the more emerging alternative based on solid-state PCMs (SS-PCMs). These latter thermomaterials undergo solid-to-solid (instead of solid to liquid) phase transitions either upon application of pressure (known then as barocaloric materials¹²⁻¹⁴) or upon change in ambient temperature (known as passive thermal energy storage materials¹⁵,¹⁶). These phase transition have associated large thermal changes which are of interest for cooling and heating technologies. Among other advantages, these thermomaterials with solid-to-solid phase transitions avoid any leak to the environment of hazardous and/or greenhouse fluids, avoid any incongruent melting and exhibit lower thermal changes that would result into more compact devices and less mechanical strain of such devices.

In the recent years, barocaloric materials have been demonstrated to exhibit very large thermal changes, some of them as large as $\Delta S \ge 100 \text{ J K}^{-1} \text{ kg}^{-1}$ (the same order of magnitude as refrigerant gases) under the application of low/moderate pressures of $p = 70–1000 \text{ bar}$ (while refrigerant gases normally operate at $p \le 150 \text{ bar}$).¹⁷⁻³⁴ These barocaloric materials belong to many different families of compounds, such as ammonium or phosphate salts³⁵⁻⁴³, superionic conductors⁴⁴,⁴⁵, spin crossover materials⁴⁶⁻⁵², n-alkanes⁵³, hybrid organic-inorganic materials⁵⁴⁻⁶¹, organic plastic crystals⁶²⁻⁶⁹ and polymers⁷⁰⁻⁷². Even more recently, barocaloric effects have been combined with gas adsorption/desorption processes in solid-to-solid breathing-transitions in MOFs, giving rise to larger thermal changes of $\Delta S \sim 300 \text{ J K}^{-1} \text{ kg}^{-1}$ under pressures as small as $p \le 16 \text{ bar}$.⁷³

Meanwhile, for thermal energy storage, the most studied SS-PCMs are typically organic polyalcohols (such as, pentaerythritol, neopentylglycol, trimethylol propane and pentaglycerine) with plastic crystal transitions.⁷⁴,⁷⁵ However, those materials work at room temperature and above, therefore, are not adequate for cold-storage applications.

Very recently, some of these SS-PCMs used for thermal energy storage applications at high temperature have been also tested for barocaloric applications.⁶²,⁶³,⁷⁶ However, due to the very large thermal hysteresis (generally required for thermal energy storage), the operating pressure needed to reach reversible thermal changes increases well above $p \ge 2500 \text{ bar}$, which in turn increases the energy demand of a given device.⁶²,⁶³,⁷⁶

These recent findings suggest that materials useful for passive thermal energy storage (including cold-storage) will exhibit some

important limitations for barocaloric refrigeration. For example, cold storage materials require relatively large thermal hysteresis for separating the temperatures of heat storage and heat release. Meanwhile, barocaloric materials ideally present ultralow thermal hysteresis (of only a few degrees or even smaller), in order to minimize the pressure required for driving the transition in a reversible manner.⁷⁷,⁷⁸ For instance, the plastic crystal neopentylglycol exhibits a large thermal hysteresis of around 15 °C, which makes it suitable for thermal energy storage applications.⁷⁴,⁷⁵ However, this large hysteresis makes necessary to apply a pressure of around 2500 bar in order to get reversible barocaloric effects in this neopentylglycol compound.⁶² On the contrary, the hybrid perovskite [TPrA][Mn(dca)₃] exhibits an ultralow thermal hysteresis of 0.9 °C, which give rise to large and reversible barocaloric effects under the application of only 70 bar.⁵⁴ Nevertheless, this low thermal hysteresis prevents its use for any type of thermal energy storage.

In this work, we show that it is possible to find thermomaterials that combine moderate thermal hysteresis with a large pressure responsiveness, which make them good candidates for multipurpose heating, cooling and thermal energy storage (or cold-storage in the specific case of this study). Here, we focus on the salt di-*n*-butylammonium tetrafluoroborate ([DBA][BF₄]). This compound, initially prepared by K. Syed Mohamed in 1988⁷⁹, shows two solid-to-solid phase transitions in the region of interest for cold-storage and food preservation. These transitions, which were first studied by G. Zabinska in 1988, take place at $T_{11} \sim 268 \text{ K}$ and at $T_{12} \sim 284 \text{ K}$, with an associated enthalpy change of $\Delta H_1 \sim 11.7 \text{ kJ mol}^{-1}$ and $\Delta H_2 \sim 2.7 \text{ kJ mol}^{-1}$, respectively.⁸⁰ To our knowledge, no further information about other thermal parameters (such as the thermal hysteresis or specific heat, for instance) or any structural information about the phase transitions has been reported.

For that reason, we have investigated the thermal properties of this material in order to evaluate its interest for passive cold-storage, as well as for barocaloric heating and cooling, using differential scanning calorimetry (DSC) at different pressures. Furthermore, we have performed a detailed structural characterization across the phase transitions, using single crystal X-ray diffraction (SC-XRD), variable-temperature synchrotron powder X-ray diffraction (VT-SPXRD) and solid-state nuclear magnetic resonance spectroscopy (SS-NMR), among others.

Remarkably, this material exhibits a large capacity for cold-storage of $E \sim 135 \text{ kJ kg}^{-1}$ in the range of 250-310 K, as well as very large barocaloric effects of $\Delta S \sim 200 \text{ J K}^{-1} \text{ kg}^{-1}$ at 275 K under the application of only 500 bar of pressure. We furthermore investigate the relationship between the structure and thermal properties, which will help the future design of multipurpose thermomaterials with enhanced properties.

## 2. Experimental section

### 2.1. Experimental section

Commercially available reagent grade di-*n*-butylamine, (≥99% Sigma-Aldrich) and tetrafluoroboric acid solution, HBF₄, (48% wt.% in H₂O Sigma-Aldrich) were used as starting materials without further purification.

Journal Name

Journal of Materials Chemistry A Accepted Manuscript

<figure><img src="image_3.png" /></figure>
## 2.2 Synthesis of [DBA][BF₄]

The synthesis procedure was similar to that previously reported⁸⁰ with a slight modification on the recrystallization process. In a typical experiment, a solution of the amine (4.6 mmol) in water (10 ml) was prepared. Then, HBF₄ (4.6 mmol) was slowly added to the previous solution at room temperature and cooled using an ice-water bath. The solvent was removed using a rotary evaporator. The obtained polycrystalline and colorless powder was stored inside a vacuum desiccator. The obtained material was highly hygroscopic, remaining wet after vacuum treatment. For that reason, the obtained powders were dispersed and recrystallized in p-xylene and, subsequently filtered and dried (Yield: 82%). This treatment appears to provide this material with a long-term protection (around few months) against water adsorption.

## 2.3. Compositional analysis

Elemental analysis was performed in an Elemental Analyzer ThermoQuest Flash EA 1112 using around 2 mg of sample (calc. %C: 44.27, %H: 9.29, %N: 6.45, found %C: 44.13, %H: 9.59, %N: 6.22).

## 2.4. Nuclear magnetic resonance

Liquid-state ¹H-NMR spectra were recorded for samples dissolved in CDCl₃ (98.5% ACROS Organics) using a Bruker Advance 300 MHz equipped with a dual cryoprobe in the University of A Coruña. Chemical shifts were expressed as parts per million (ppm, δ) downfield of tetramethylsilane (TMS) (see Figure S1 of ESI). The description of signals included: t = triplet and m = multiplet. All coupling constants were absolute values and were expressed in Hertz (Hz). δ 6.78 (s, 2H), 3.04 (t, J = 7.8 Hz, 4H), 1.72 (m, J = 7.5 Hz, 4H), 1.39 (m, J = 7.5 Hz, 4H), 0.94 (t, J = 7.2 Hz, 6H) ppm.

Solid-state ¹H single-pulse Magic Angle Spinning (MAS) spectra and ¹³C Cross-Polarization Magic Angle Spinning (CPMAS) NMR were performed in a Bruker NEO 750 spectrometer (B₀ = 17.6 T) using 1.3 mm zirconia rotors spun at *ν*<sub>MAS</sub> = 50 and 15 kHz respectively. Both experiments were carried out at 318 K in the University of Santiago de Compostela.

Variable temperature solid-state NMR (SS-NMR) MAS experiments were performed in a Varian-Agilent 600 WB spectrometer (B₀ = 14.1 T) using 3.2 mm zirconia rotors spun at different frequencies (*ν*<sub>MAS</sub> = 5.5, 10 or 15 kHz) (located in the CNRS campus of Montpellier). ¹H single pulse spectra were recorded using π/2 flip angle pulses of 2.5 µs and recycling delays of 5 s. ¹H spin-echo spectra were recorded using MAS synchronized interpulse delays, and π/2 and π flip angle pulses of 2.5 and 5.0 µs respectively. Chemical shifts were referenced using adamantane (δ<sub>iso</sub> = 1.8 ppm) as external reference. ¹¹B single pulse MAS spectra were recorded using π/10 flip angle pulses of 1 µs and recycling delays of 0.5 s. Chemical shifts were referenced using dry NaBH₄ (δ<sub>iso</sub> = -42.1 ppm) as external reference. The temperatures were adjusted using chilled compressed air. The differences between *T* setpoints and sample temperatures are below 5 K. Pb(NO₃)₂ was used for the calibration of the sample temperatures following a previously published procedure.⁸¹

## 2.5. Thermal analysis

Variable-temperature differential scanning calorimetry (VT-DSC) in the lower-temperature range (115-150K) at ambient pressure was performed in a Netzsch DSC 214 Polyma at KTH Royal - Institute of Technology. Meanwhile, VT-DSC in the higher-temperature range (250-310 K) at ambient pressure was carried out in using a TA-Instruments Q2000 coupled with a RCS 90 cooler. In both cases, around 10 mg of sample were tested under 10 K min⁻¹ heating/cooling rate and under nitrogen atmosphere in the Department of Naval and Industrial Engineering of the University of A Coruña.

Heat capacity at ambient pressure, Cₚ, was calculated by the classic ASTM E1269 standard test method, by which the measured heat flow of the test material is compared to that of a reference material according to⁸²,⁸³

$$Cp_m = C p_s \frac{D_m m_s}{m_m D_s} \quad (1)$$

Where Cₚₘ is the [DBA][BF₄] heat capacity, Dₘ is the vertical displacement between the baseline and the sample heat flow curve, Cₚₛ is the tabulated heat capacity of the sapphire, mₛ is the sapphire mass, mₘ is the sample mass, and Dₛ is the vertical displacement between the baseline and the measured sapphire heat flow.

Variable pressure differential scanning calorimetry (VP-DSC) was performed at various applied pressures from 1 to 1000 bar using a Setaram µDSC7 EVO microcalorimeter equipped with a Teledyne ISCO 65D pressure pump (Queen Mary University of London). Samples of mass ~ 20 mg were tested in the temperature range of 250 – 310 K, using a 1.2 K min⁻¹ heating/cooling rate under isobaric conditions. For evaluating long-term cyclability, we have performed an over 30 hours study under isothermal conditions at 286 K using 10 bar min⁻¹ continuous pressurization/depressurization ramps.

The used DSC equipment shows a temperature accuracy of ± 0.1 °C and enthalpy accuracy of ± 0.5 %.

The thermal conductivity was evaluated at ambient conditions using a C-Term Trident thermal conductivity analyser equipped with a 6 mm flexible transient plane source (TPS) sensor.

## 2.6. Variable temperature synchrotron powder X-ray diffraction (VT-SPXRD).

VT-SPXRD patterns were recorded at the I11 beamline of the Diamond Light Source Synchrotron (Harwell, UK) using a wavelength of λ = 0.825955(3) Å. The wavelength was determined by Rietveld refinement of a NIST640D silicon standard at room temperature. The sample was enclosed in a glass capillary (inner diameter φ = 0.5 mm) and kept in continuous rotation during data collection to improve powder averaging. Patterns were collected using Mythen position sensitive detector while heating the sample from 170 K to 335 K. The working temperature was set using a FMB Oxford hot-air blower. Rietveld and Le Bail refinements were carried out using the program GSAS-II.⁸⁴

## 2.7. Single crystal X-ray diffraction (SC-XRD).

Single-crystal X-ray diffraction experiments were carried out at two temperatures well above and below the observed DSC transitions. Two single-crystal data sets of the same crystal were collected at 100 K and 240 K in a Bruker D8 VENTURE Kappa

View Article Online
DOI: 10.1039/D3TA04063A

Journal of Materials Chemistry A Accepted Manuscript

<figure><img src="image_4.png" /></figure>
<table><caption>Table 1. Thermodynamic parameters of observed phase transitions by DSC data for [DBA][BF4] compound.</caption><thead><tr><th>Transition</th><th>T<sub>II(h)</sub> (K)</th><th>T<sub>II(c)</sub> (K)</th><th>ΔH<sub>(h)</sub> (kJ kg<sup>-1</sup>)</th><th>ΔH<sub>(c)</sub> (kJ kg<sup>-1</sup>)</th><th>ΔS<sub>(h)</sub> (J kg<sup>-1</sup> K<sup>-1</sup>)</th><th>ΔS<sub>(c)</sub> (J kg<sup>-1</sup> K<sup>-1</sup>)</th><th>ΔS<sub>(h)</sub> (J kg<sup>-1</sup> K<sup>-1</sup>)</th><th>ΔS<sub>(c)</sub> (J kg<sup>-1</sup> K<sup>-1</sup>)</th></tr></thead><tbody><tr><td>PI ↔ PII</td><td>282</td><td>279</td><td>12.7</td><td>13.7</td><td>45.0</td><td></td><td>49.1</td><td></td></tr><tr><td>PII ↔ PIII</td><td>269</td><td>263</td><td>58.3</td><td>56.6</td><td>216.7</td><td></td><td>215.2</td><td></td></tr><tr><td>PIII ↔ PIV</td><td>134</td><td>129</td><td>0.37</td><td>0.43</td><td>2.65</td><td></td><td>3.24</td><td></td></tr></tbody></table>

Note: (h) = heating, (c) = cooling.

X-ray diffractometer equipped with a PHOTON III detector and using monochromatic CuKα radiation (λ = 1.54178 Å) at the University of Santiago de Compostela. Several crystals among the best-looking ones in the crystallization batch were tested and one was selected for data collection. It was mounted on a MiTeGen MicroMount™ using Paratone® N (Hampton Research). The crystal temperature was changed and stabilized with a cold stream of nitrogen gas from an Oxford Cryosystem 800plus cooler. The data integration and reduction were performed using the APEX3 v2019.1-0 (Bruker AXS, 2019) software suite. The intensities collected were corrected for Lorentz and polarization effects and for absorption by semi-empirical methods based on symmetry-equivalent data using TWINABS 2012/1 of the software suite. The structures were solved by the novel dual-space algorithm implemented in SHELXT2018/2⁸⁵ program and were refined by least squares method on SHELXL2018/3⁸⁶. As the presence of twinning in all data collected was clear from visual inspection of diffraction images, all datasets were indexed using CELL NOW 2008/4. The four primary twin domains at both temperatures could be indexed, and the corresponding orientation matrices obtained, although the diffraction patterns displayed the presence of several additional domains with minor contribution. The integrations of the reflections were performed with SAINT 8.40B considering simultaneously the predictions given by the four orientation matrices corresponding to the twin domains found. The software used to treat the twinned data is included in the APEX3 suite.

The electron density maxima corresponding with the hydrogen atoms of the NH₂⁺ group of the amine cation were clearly visible in the (Fo-Fc) Fourier map. The rest of the hydrogen atoms of the DBA cation were introduced at idealized positions. All hydrogen atoms were restrained during refinement using the riding model implemented in SHELXL2018/3.

Detailed experimental crystallographic data for phases IV (T = 100 K) and III (T = 240 K) of the [DBA][BF₄] compound are included in Supporting Information, Table S2 of ESI. CCDC-2213758 (phase IV) and CCDC-2213757 (phase III) contain the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/structures.

## 3. Caloric properties

### 3.1. Cold-storage properties

To analyse the cold-storage properties of the material upon changes in the ambient temperature, we performed VT-DSC at ambient pressure. The results reveal that [DBA][BF₄] undergoes three reversible first-order phase transitions in the temperature range studied (115 - 310 K), see Figure 1 and Table 1. The first two transitions occur at T₁₁ ~ 282 K and T₁₂ ~ 269 K on heating, in agreement with a previous report.⁸⁰ These transitions exhibit a thermal hysteresis of 3 K and 6 K in the case of T₁₁ and T₁₂, respectively. Furthermore, we also found a new phase transition at T₁₃ ~ 134 K with a thermal hysteresis of 5 K, which has not been previously reported.

These transitions occur between different phases (PI ↔ PII, PII ↔ PIII, and PIII ↔ PIV) as we will further discuss in the following sections. Interestingly, we observe that the two transitions occurring at higher temperature, exhibit very large latent heats of ΔH₁ ~ 13 kJ kg⁻¹ and ΔH₂ ~ 58 kJ kg⁻¹. Even more interestingly, these two transitions are very close in temperature and appear in the range of 250-310 K, covering the temperature range adequate for food preservation and cold-storage. Furthermore, even if the thermal hysteresis is still moderate for cold-storage, both transitions occur very close in temperature and can be combined for thermal energy storage applications.

<figure><img src="image_5.png" /><figcaption>Figure 1. VT-DSC curves for the [DBA][BF4] compound in the (a) lower- and (b) higher-temperature range.</figcaption></figure>

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>Figure 2. Specific heat capacity, C_p, measurements of [DBA][BF_4] material from 250 to 295 K on heating. Note: solid red lines represent the calculated C_p for the different polymorphs in the studied temperature range, meanwhile dashed lines represent the discontinuity of the C_p across the transitions. Black dashed lines represents the transition peaks.</figcaption></figure>

In that way, the heat absorbed by the combination of both transitions will be maximized for a charging temperature of ~290 K, meanwhile the stored heat would not be fully released until the temperature decreases back to ~250 K. This combined use of both close transitions increases the separation between the full charging and discharging temperatures beyond the individual thermal hysteresis of both transitions.

With this in mind, we calculate the thermal storage capacity, E, of the material considering the latent heat of both combined transitions and the heat capacity heat throughout temperature range from 250 to 310 K, following<sup>87</sup>

$$E = m\left[\left[\int_{T_{i1}}^{T_{m1}} C_p dT\right] + \Delta H_1 + \left[\int_{T_{m1}}^{T_{m2}} C_p dT\right] + \Delta H_2 + \left[\int_{T_{m2}}^{T_{m3}} C_p dT\right]\right] + m\Delta H_3 \quad (2)$$
DOI: 10.1039/D3TA04063A

where C_p is the calculated specific heat capacities for the different polymorphs observed between 250 and 295 K (see red lines in Figure 2); m is the mass of the sample, and ΔH is the latent heat obtained for the two phase transitions by the integration of the transition's peaks, which are very close in temperature.

Accordingly, the sensible heat for the temperature range of 250-295 K can be estimated as the combination of the specific heat integration (red line) and the total enthalpy change (peak's area integration), which gives a storage capacity as large as E = 134.6 kJ kg⁻¹. If we consider only the total enthalpy change (ΔH₁+₂ ~ 71.0 kJ kg⁻¹), the observed value is very close to many cold-storage materials reported in the literature (such as polyglycol E400 (ΔH = 99.6 kJ kg⁻¹), isopropyl palmitate (ΔH = 95 kJ kg⁻¹) or K₂HPO₄·6H₂O (ΔH = 109 kJ kg⁻¹), and even to some commercially ones, such as TH31 (ΔH = 131 kJ kg⁻¹), ClimSel C7 (ΔH = 130 kJ kg⁻¹) or MPCM-30 (ΔH = 140 kJ kg⁻¹).<sup>[6]</sup>

Furthermore, according to TGA results, the [DBA][BF₄] compound also remains thermally stable up to ~450 K, similar to other commercial PCMs, see Figure S2 of ESI.

## 3.2. Barocaloric properties

To study the barocaloric properties of [DBA][BF₄], we use VP-DSC at various applied pressures from 1 to 1000 bar in the temperature range of 250 - 310 K, focusing on those transitions with larger thermal changes, T₁₁ and T₁₂ as previously discussed. Both transitions (between the PIII ↔ PII and PII ↔ PI, as we will detail in the following sections) exhibit conventional barocaloric effects, with an increase of the transition temperature upon pressurization, see Figure 3a.

<figure><img src="image_8.png" /><figcaption>Figure 3. (a) VP-DSC curves of [DBA][BF₄] in isobaric conditions from 1 bar to 1000 bar. (b) Temperature of the maximum of each peak of the VP-DSC curves on cooling and heating as a function of the applied pressure. (c) ΔS<sub>ib</sub> curves at different pressures during heating and (d) cooling.</figcaption></figure>

<figure><img src="image_9.png" /></figure>
Remarkably, this material has a particular behaviour predicted by M. V. Gorev et al.⁸⁸ where a triple-point between phases PI, PII and PIII can be observed in the temperature-pressure phase diagram (see Figure 3b). To the best of our knowledge, this is the first material reported with the merging transitions phenomenon where both transitions (PIII ↔ PII and PII ↔ PI) show direct barocaloric effect. In that way, both transitions are fully merged at 600 bar and a single transition is observed for $p \ge 600$ bar, which involves a PIII ↔ PI phase transition. This merging is due to the existence of a triple-point, where the stability of phase 2 decreases when pressurizing, until it disappears at pressures above 600 bar. This behaviour, although very scarce, already appear in a few barocaloric materials.⁸⁸

The barocaloric tunability or the barocaloric coefficient ($dT_{t1}/dp$) of each phase transition is calculated from the slope of the linear regions at Figure 3b. Due to the merging effect, we can separate three different barocaloric coefficients in heating and cooling.

Below 600 bar, there are two barocaloric coefficients of $dT_{t1}/dp \sim 5.5$ K kbar⁻¹ and $dT_{t2}/dp \sim 29$ K kbar⁻¹ for the PII ↔ PI and PIII ↔ PII phase transition, respectively. This noticeable difference between the barocaloric coefficients is responsible for the phase transitions merging into a new single-phase transition for $p \sim 600$ bar. This new PIII ↔ PI phase transition exhibits a new barocaloric coefficient of $dT_{tm}/dp \sim 22.6$ K kbar⁻¹.

The calculated latent heat and the entropy change of the transitions $T_{t1}$ and $T_{t2}$ are similar to the one previously described using VT-DSC at ambient pressure. Meanwhile, the values for the new PIII ↔ PI phase transition is approximately the sum of the values of $\Delta H$ and $\Delta S$ of the

individual phase transition, reaching values of 75.0 J g⁻¹ and 263.1 J K⁻¹ kg⁻¹ (see the structural characterization for more details).

In addition, we calculated the barocaloric effects of the [DBA][BF₄] compound in terms of isothermal entropy change, $\Delta S_{it}$, by using quasi-direct methods.⁵⁹,⁸⁹ For that purpose, we first estimated the isobaric entropy change ($\Delta S_{ib}$) as a function of temperature following

$$ \Delta S_{ib} = S(T_f,p) - S(T_0,p) = \int_{T_0}^{T_f} \frac{1}{T} \left[ C_p(T, p_{atm}) + \frac{dQ}{dT}(T, p) \right] dT \quad (3) $$

where $T_0$ and $T_f$ are the starting and end temperature chosen by eye below and above of the transition, respectively ($T_0 = 250$ K and $T_f = 310$ K), $p$ is the constant pressure of each isobaric curve, $p_{atm}$ is atmospheric pressure and $C_p$ is the heat capacity of the material for each temperature at atmospheric pressure. Here, it should be noted as previously observed in Figure 2 that $C_p$ is not a constant value, but is temperature dependent and specific of each polymorph. Therefore, the total isobaric entropy change is a contribution of the specific heat capacity of each polymorph and of the heat flow observed at each transition, as described elsewhere.⁶⁴

Figure 3c and 3d shows the obtained $\Delta S_{ib}$ versus temperature at different applied pressures on heating and cooling. For the pressure range of 1–400 bar, the $\Delta S_{ib}$ vs temperature curves exhibit a two step-like curve, in where the largest step is due the PIII ↔ PII phase transition with $|\Delta S_{ib}| \sim 211$ J K⁻¹ kg⁻¹ (at 1 bar) and a small step is related to the PII ↔ PI phase transition with $|\Delta S_{ib}| \sim 36.0$ J K⁻¹ kg⁻¹ (at 1 bar). For pressures above 400 bar, the $\Delta S_{ib}$ vs temperature curves exhibit a single and sharp increase of $\Delta S_{ib}$ due to the PIII ↔ PI phase transition with $|\Delta S_{ib}| \sim 254$ J K⁻¹ kg⁻¹ (at 1000 bar).

<figure><img src="image_10.png" /><figcaption>Figure 4. (a) Pressure-driven reversible isothermal entropy changes on applying ($0 \to p$) and removing ($p \to 0$) pressure. (b) Reversible adiabatic temperature changes on applying ($0 \to p$) and removing ($p \to 0$) pressure at different pressures. (c) Temperature span as a function of applied pressure for fixed intervals of reversible isothermal entropy change and (d) reversible adiabatic temperature changes for [DBA][BF₄].</figcaption></figure>

<figure><img src="image_11.png" /></figure>
Subsequently, the isothermal entropy change (ΔS<sub>it</sub>) is calculated by difference following

$$ \Delta S_{it} = \Delta S_{ib}(p \neq 1, T) - \Delta S_{ib}(p = 1, T) \quad (4) $$

It should be noted that, similarly to other rather flexible barocaloric materials, additional changes in isothermal entropy<sup>36</sup> were not considered since the thermal expansion coefficient might significantly decrease under pressurization,<sup>57</sup> which can lead to overestimations of the total entropy changes.

The difference between the cooling curve at p > 1 bar and the heating curve at p = 1, bar, gives the isothermally reversible entropy change (ΔS<sub>rev</sub>). Figure S3a of ESI and Figure 4a shows the ΔS<sub>it</sub> and ΔS<sub>rev</sub> curves from 200 bar to 1000 bar, respectively. Interestingly, a ΔS<sub>rev</sub> > 200 J K<sup>-1</sup> kg<sup>-1</sup> is achieved for pressures higher than 500 bar and the largest value ΔS<sub>rev</sub> = ~237.7 J K<sup>-1</sup> kg<sup>-1</sup> is obtained at 1000 bar at T= 285 K (see Figure S4 of ESI).

Furthermore, we can estimate the reversible adiabatic temperature change (ΔT<sub>rev</sub>) of the material according to

$$ |\Delta T_{rev}(S; p_{atm} \leftrightarrow p)| = |T_{C}(S; p) - T_{H}(S; p_{atm})| \quad (5) $$

where T<sub>C</sub>(S; p) is each temperature of the cooling isobaric entropy change (see Figure 4b) at the high pressure and T<sub>H</sub>(S; p<sub>atm</sub>) is each temperature of the heating isobaric entropy change (see Figure 4a) at the atmospheric pressure. As it is shown in Figure 4b, the largest ΔT<sub>rev</sub> is about ~17 K at 1000 bar, a difference that is usable in cooling devices that work in sequential compression-decompression cycles. Figure 4c and 4d shows the operating temperature range or span temperature (T<sub>span</sub>) as a function of the applied pressure where a given interval of isothermal entropy change or reversible adiabatic temperature change values can be achieved. In both cases, T<sub>span</sub> increases almost linearly with applied pressure. At p = 1000 bar, under isothermally reversible conditions, T<sub>span</sub> = 13 K with reversible ΔS<sub>it</sub> = 200 J K<sup>-1</sup> kg<sup>-1</sup>, while in the adiabatically reversible case, T<sub>span</sub> = 10 K with reversible ΔT<sub>rev</sub> = 10 K. This barocaloric performance is in a temperature range from 267 K to 290 K (-6 °C to 17 °C).

The properties of the material presented in this section are quite promising for barocaloric cooling/heating applications: Its temperature transition is slightly below room temperature, an

<figure><img src="image_12.png" /><figcaption>Figure 5. Heat flow on pressurization/depressurization cycles over time at 286 K from 1 to 1000 bar.</figcaption></figure>

excellent range for refrigeration applications. The material presents a high sensitivity to applied pressures, showing a barocaloric coefficient of 28.9 K kbar<sup>-1</sup> for PIII ↔ PII transition. This transition is higher than barocaloric coefficients reported up to date, and 22.6 K kbar<sup>-1</sup> for PIII ↔ PI transition when both merged. The material also presents other excellent barocaloric parameters such as ΔS<sub>ib</sub> > 250 J K<sup>-1</sup> kg<sup>-1</sup> and ΔT<sub>rev</sub> > 15 K with a wide operational temperature.

As noted above, the active temperature range of this material corresponds well to the needs of food preservation. Relatively few barocaloric materials have previously been reported in this temperature range.<sup>4,5</sup>

Furthermore, in order to evaluate the long-time performance, we have carried out direct calorimetry studies under several pressurization and depressurization cycles under isothermal conditions. The selected temperature was 286 K, where we can take advantage of the thermal changes of both merged transitions. Figure 5 shows that the material cyclability remains stable up-to, at least, 30 hours with constant values of ΔS ~270 J K<sup>-1</sup> kg<sup>-1</sup>, which are slightly larger than those observed by quasi-direct methods.

For comparison purposes, Figure 6a shows the transition temperature and the entropy change of [DBA][BF<sub>4</sub>] together with a selection of the best barocaloric materials reported to date. In this graphic we only include pure barocaloric materials with solid-to-solid phase transitions showing colossal thermal changes of ΔS > 100 J K<sup>-1</sup> kg<sup>-1</sup>.

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>Figure 6. (a) All barocaloric materials temperature transition (T<sub>t</sub>) and entropy change (ΔS) comparison. (b) Operating temperature range and the needed pressure of the largest entropy change materials (>100 J K<sup>-1</sup> kg<sup>-1</sup>). See more information of all materials represented in Table S2 of ESI.</figcaption></figure>
<table><caption>**Table 2.** Main structural features of the different polymorphs observed for [DBA.BF₄].</caption><thead><tr><th>Polymorph</th><th><i>T</i> (K)</th><th>Space group</th><th>[DBA]<sup>+</sup> configuration</th><th>[BF<sub>4</sub>]<sup>+</sup> configuration</th><th>DOI: 10.1039/D6TA0033A</th></tr></thead><tbody><tr><td>PIV</td><td>&lt; 134</td><td><i>Ccce</i> (orthorhombic)</td><td>Ordered zig-zag chain</td><td>Distorted tetrahedral</td><td>SC-XRD</td></tr><tr><td>PIII</td><td>134-260</td><td><i>Cmme</i> (orthorhombic)</td><td>Ordered zig-zag chain</td><td>4-fold disorder of F atoms</td><td>SC-XRD, SS-NMR</td></tr><tr><td>PII</td><td>260-280</td><td><i>Cmme</i> (orthorhombic)</td><td>Disorder of terminal -CH<sub>3</sub></td><td>Full rotation of F atoms</td><td>VT-SPXRD, SS-NMR</td></tr><tr><td>PI</td><td>&gt; 280</td><td><i>I4/mmm</i> (tetragonal)</td><td>Disorder of complete butyl chains</td><td>[BF<sub>4</sub>]<sup>+</sup> split into 2 sites</td><td>VT-SPXRD, SS-NMR</td></tr></tbody></table>

Note: SC-XRD = single crystal X-ray diffraction, (h) = heating, (c) = cooling.

We omit other materials that exhibits barocaloric effects due to solid-to-liquid phase transitions (such as n-alkanes with $\Delta S > 700$ J K⁻¹ kg⁻¹)⁹⁰ or that are combined with simultaneous adsorption/desorption thermal changes (such as the breathing-caloric MIL-53 with $\Delta S > 300$ J K⁻¹ kg⁻¹).⁷³

Meanwhile Figure 6b shows the operating temperature range of the [DBA][BF₄] together with other barocaloric materials with large entropy changes (>100 J K⁻¹ kg⁻¹) as a function of the applied pressure, where the most relevant materials are those which need less pressure to show large reversible thermal changes.

In addition, we have also explore the thermal conductivity of the material that shows values of κ ~0.2 W m⁻¹ K⁻¹ similar to the rest of non-metallic barocaloric materials.¹³ For practical refrigeration applications, even if the thermal conductivity is relatively low, it can be improved by creating composites with thermal conductive materials, such as carbon materials.¹³

## 4. Structural characterization

In order to deep into the structural dynamics and origin associated to the observed thermal changes, we perform single crystal and synchrotron powder X-ray diffraction, as well as solid-state ¹H-, ¹³C- and ¹¹B-NMR. As will be later explain in detail, the thermal changes are found to be related to a progressive increase in the structural disorder of both di-*n*-butylammonium cations and tetrafluoroborate anions.

Single crystal X-ray diffraction (SC-XRD) allowed the elucidation of the crystal structures of the low-temperature phases IV and III, which exhibit low or even negligible disorder. Nevertheless, at higher temperatures the presence of twinning and strong diffuse scattering (characteristic of plastic crystal transitions) hindered resolving the crystal structure.

<figure><img src="image_15.png" /><figcaption>**Figure 7.** Structure of [DBA][BF₄] at 100 K (PIV), where the tetrahedra correspond to [BF₄]⁻ anions, the chains to [DBA]⁺ cations and the black cage represents the unit cell.</figcaption></figure>

Therefore, the structures with higher disorder (PI, PII and PIII) were further studied by using variable temperature synchrotron powder X-ray diffraction (VT-SPXRD) and solid-state ¹H-, ¹³C- and ¹¹B-NMR. Table 2 summarizes the main features of the different polymorphs.

### 4.1. SC-XRD characterization (PIV and PIII)

**Phase IV (T < 134 K).** This phase was elucidated by single crystal X-ray diffraction. It should be noted that this PIV is here described for the first time. The polymorph crystallizes in an orthorhombic symmetry, with space group *Ccce* and lattice parameters *a* = 7.3980(5) Å, *b* = 18.4562(11) Å and *c* = 16.8535(10) Å (at T = 100 K, see more details in Table S1 of ESI and Figure 7).

The asymmetric unit of phase IV consists of one half of [DBA]⁺ and [BF₄]⁻ ions, as both moieties are placed over a 2-fold symmetry axis. In the case of the [BF₄]⁻ anions, the boron atoms are in a slightly distorted tetrahedral environment with two different B-F bond lengths (1.404(5) Å and 1.380(5) Å). Meanwhile, the [DBA]⁺ cations show a perfect zig-zag configuration of the alkylammonium chains (see Figure 7).

Concerning the molecular arrangement within the crystal structure, this compound displays a significantly distorted NaCl-type structure where the [DBA]⁺ cations are surrounded by six [BF₄]⁻ anions in a distorted octahedral environment, and vice versa (see Figure S5 of ESI).

The interactions between the anionic and cationic components are electrostatic and via H-bonds. In this latter case, the donor atoms are the N-atoms of the amine groups and the acceptor atoms are the F-atoms of [BF₄]⁻ anions (see Figure S6 of ESI). That the distances between hydrogen and fluorine atoms are around 2.19 Å and 2.30 Å; the two fluorine atoms involved in the H-bonds are those with the longest B-F bond lengths (1.405(5) Å).

**Phase III (134 K < T < 260 K).** According to SC-XRD phase III also crystallizes in an orthorhombic crystal system with space group *Cmme* and lattice parameters *a* = 7.5762(10) Å, *b* = 18.642(3) Å and *c* =

<figure><img src="image_16.png" /><figcaption>**Figure 8.** Crystal structure of [DBA][BF₄] at 240 K (PIII), where the [BF₄]⁻ anions are represented in four different disordered orientations, the chains correspond to [DBA]⁺ cations and the black cage represents the unit cell. **Note:** One [DBA]⁺ and [BF₄]⁻ species are shown as ORTEP drawing (ellipsoid probability at 50%).</figcaption></figure>

<figure><img src="image_17.png" /></figure>
8.5564(12) Å at T = 240 K (see more details in Table S1 of ESI and Figure 8).

The obtained crystal structure of this higher temperature phase III is very similar to that previously described for phase IV, although the local symmetry in the environment of both ions has increased from point group 2 (C2) to mm2 (C2v). The molecular geometry of both ions is compatible with this increased symmetry; however, while the DBA cations remain ordered, the BF4 anions are now best represented by a fourfold disorder model. The two fluoride atoms involved in N-H···F hydrogen bonds display longer B-F bond distances (1.381(14) and 1.384(13) Å) compared to the other two (1.352(13) and 1.360(15) Å) that are the acceptors of weaker C-H···F interactions. The [BF4]⁻ anions remain recognisably tetrahedral on average, rather than appearing spherical as in the case of full orientational order; this could explain the small enthalpy and entropy changes experimentally observed for the transition between phases IV and III.

## 4.2. VT-SPXRD (PII, PI) and SS-NMR (PIII, PII, PI) characterization

As previously mentioned, the progressive increase in the structural disorder hinder the structural characterization through SC-XRD and, therefore, we used VT-SPXRD for completing the structural characterization of PII and PI. Moreover, we use SS-NMR to deep into further structural details of the also disordered PIII. These analyses reveal a progressive disorder of the [DBA]⁺ cations and [BF4]⁻ anions, as represented in Figure 8, which will be the structural origin of the entropy changes observed by VT-DSC.

**VT-SPXRD characterization.** To characterize the high temperature crystalline phases and their phase transitions, VT-SPXRD experiments were carried out from 180 K to 320 K.

From the analysis of VT-SPXRD patterns (see Figure S7 of ESI), two phase transitions can be clearly observed at $T_{t1}$ ~280 K (PI ↔ PII) and $T_{t2}$ ~260 K (PII ↔ PIII), in fully agreement with VT-DSC data. These data were fitted using the Le Bail method and PIII as a starting model (Figure S8 of ESI). Subsequently, PI and PII were indexed by comparison with PIII.

Here, PII shows the same orthorhombic space group *Cmme* as PIII (see Figure S9 of ESI) with a large increase of the *b* axis (from *b* ~ 18.5 Å to *b* ~ 20.3 Å) and, in turn, a very large volume increase of ~8.4%. Meanwhile, PI shows a body-centred tetragonal symmetry with space group *I4/mmm*, where the lattice

<figure><img src="image_18.png" /><figcaption>**Figure 9.** Scheme with the molecular configuration of ions for the different phases of [DBA][BF4].</figcaption></figure>

parameters are *a* ~ 5.7 Å and *c* ~ 20.5 Å (Figure S10 of ESI). Here, the basis vectors are related by: DOI: 10.1039/D3TA04063A

$$ (a \ b \ c)_o = (a \ b \ c)_t \begin{pmatrix} 1 & 0 & 1 \\ 1 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}. $$

Using the initial Le Bail refinement, we propose structural models for PI and PII, which were generated and refined by the Rietveld method. For phase II, we used the known crystal structure obtained from SC-XRD of phase III as the initial model. The obtained crystal structure and details of the refinement are summarized in Figure S11 and Table S3 of ESI.

The most noticeable feature of the crystal structure of PII is the disorder of the terminal methyl group of the DBA cations, which is split into two different positions (Figure 9), similar to other 2D perovskites.⁵⁹,⁹¹,⁹² This disorder means that it is no longer possible for adjacent butyl chains along the *b* direction to overlap as in phase III (see Figure 10), and is thus responsible for the sudden increase in the *b* lattice parameter. In addition [BF4]⁻ anions show a more severe disorder than that previously observed in PIII.

Meanwhile, the disorder in PI affects the whole butyl chain. Moreover, the [BF4]⁻ ion sits slightly away from an inversion centre, so that the B atom is now disordered between two sites in addition to the disorder of the F atoms about the fourfold axis. Details of this refinement are given in Figure 9, Figure S12 and Table S4 of ESI.

These data reveal that the structural disorder of the [DBA]⁺ and [BF4]⁻ are the main responsible of the thermal changes observed by DSC. As it can be observed in Figure 9, the transition PI↔PII involve a lower structural disorder than PII↔PIII and PIII↔PIV, which explains the larger thermal changes experience in these latter transitions. Accordingly, we have focused on the thermal characterization of such transitions.

In the same line, the larger barocaloric coefficients observed by VT-DSC can be explained by the volume changes experience by the material.

<figure><img src="image_19.png" /><figcaption>**Figure 10.** Detail of the distance difference along the *b*-axis between PIII and PII.</figcaption></figure>

<figure><img src="image_20.png" /></figure>
<figure><img src="image_21.png" /><figcaption>Figure 11. Evolution of (a) the lattice parameters and (b) the volume of PI, PII and PIII of [DBA][BF₄] as a function of temperature. Note: blue, green and red points represent the lattice parameters and volume of the polymorps PIII, PII and PI, respectively. Note: error bars, as calculated by the software, are smaller than symbols.</figcaption></figure>

In that regard, Figure 11a and 11b show the evolution of the lattice parameters and the volume as a function of temperature as obtained from the Le Bail refinement of the SPXRD patterns.

All phases exhibit conventional thermal expansion upon heating in all cell parameters, except PII that exhibits negative thermal expansion along the c-axis.

The observed volume changes can be used to indirectly estimate the barocaloric coefficients associated with the different transitions following the Clausius-Clapeyron equation:¹²

$$dT/dp = \Delta V/\Delta S \quad (6)$$

where $dT/dp$ is the pressure dependence of the transition temperature (barocaloric coefficient), $\Delta V$ is the cell volume change and the $\Delta S$ is the entropy change as obtained by ambient pressure DSC. Using this indirect method, the estimated barocaloric coefficients for the two higher-temperature transitions are $dT_{11}/dp \sim 4.4$ K kbar⁻¹ and $dT_{12}/dp \sim 39.5$ K kbar⁻¹, which are in agreement with the experimental barocaloric coefficients obtained by VP-DSC, with values of $dT_{11}/dp \sim 5.5$ K kbar⁻¹ and $dT_{12}/dp \sim 29$ K kbar⁻¹. The significantly larger

barocaloric coefficient of $dT_{12}/dp$ can be attributed to a much larger volume change of ~8.4%.

In summary, the structural characterization confirm that the large entropy changes are related to structural disorder of the [DBA]⁺ cations and [BF₄]⁻ anions, while the large barocaloric coefficients are originated through very large volume changes of up to ~8.4% in the case of the transition PII↔PIII.

SS-NMR of ¹H, ¹³C and ¹¹B. In order to gain further insights of the dynamic disorder of the different ions, we perform ¹H and ¹³C{¹H} MAS spectra of [DBA][BF₄] at high magnetic field (17.6 T), where fast MAS show almost resolved peaks for [DBA]⁺ (Figure 12a and Figure S13 of ESI, respectively). Under these conditions and at room temperature, the crystal phase corresponds to PI. By decreasing MAS frequency and temperature important variations can be observed in the ¹H line-shapes of phases PI, PII and PIII (Figure 12b). These variations are related to a gradual broadening of the peaks when decreasing the temperature. Table 3 shows the transverse relaxation times $T_2'$, measured using a spin-echo pulse sequence. The main mechanisms of relaxation are here the dipolar couplings, notably between protons.

<figure><img src="image_22.png" /><figcaption>Figure 12. ¹H MAS NMR spectra: (a) Fast MAS spectrum of phase PI ($v_0$ = 750 MHz, $v_{MAS}$= 50 kHz); (b) MAS spectra of phase PI, PII and PIII obtained by decreasing temperature ($v_0$ = 600 MHz, $v_{MAS}$=15 kHz).</figcaption></figure>

<figure><img src="image_23.png" /></figure>
<table><caption>Table 3. ¹H spin-echo transverse relaxation times $T_2'$.</caption><thead><tr><th rowspan="2">Phase</th><th colspan="4">$T_2'$ (ms)</th></tr><tr><th>CH<sub>3</sub>- (1)</th><th>C-CH<sub>2</sub>-C (2-3)</th><th>N-CH<sub>2</sub>- (4)</th><th>-NH<sub>2</sub> (5)</th></tr></thead><tbody><tr><td>PI</td><td>0.66</td><td>0.50</td><td>0.65</td><td>0.75</td></tr><tr><td>PII</td><td>0.48</td><td></td><td>0.26</td><td>0.37</td></tr><tr><td>PIII</td><td></td><td>0.2-0.3</td><td></td><td>0.3-0.4</td></tr></tbody></table>

View Article Online

DOI: 10.1039/D3TA04063A

Note: Obtained using a spin-echo at $v_{MAS}$ = 15 kHz and T settled at: 278 K (PI), 263 K and 261 K (PII), 243 K (PIII).

The highest values obtained for PI can therefore be explained by a higher degree of motion of the [DBA]$^+$ cation allowing for a decrease in the effective dipolar couplings. Also, this cation appears slightly more mobile in phase II compared to phase III, but the differences are smaller than with phase I. This gradual increase in [DBA]$^+$ mobility from phase III to phase I is consistent with the decrease in density and the increase in alkyl chain disorder that can be extracted from X-rays data.

As VT-SPXRD analysis and the solid-state ¹H single-pulse MAS spectra showed, the PIII ↔ PII phase transition shows larger disorder/mobility in the terminal -CH₃, while the PII ↔ PI phase transition presents larger disorder/mobility in the rest of the chain. Meanwhile, for higher temperatures and pressures (after the triple point), the transitions will involve a direct disorder from PIII to PI. Here, the dibutylammonium chains get completely disordered in one step, rather than in a progressive manner with an initial disorder of the terminal -CH₃ groups.

On one hand, the PIII ↔ PII phase transition shows a very large volume change due to the disorder of the terminal -CH₃ of both chains of [DBA]$^+$ cation (see Figure 9 and Figure 12), while the PII ↔ PI

phase transition, involving the rest of the chain, is less relevant in volume changes (see Figure 9 and Figure 12). In that regard, the new transition PIII ↔ PI after the triple point, would have a volume change close to the volume change of PIII ↔ PII transition, which also explains the similarity in the barocaloric coefficients for both transitions. In the same line, the PIII ↔ PI phase transition involves order/disorder process affecting the four carbons of each chain, so the resulting values of $\Delta H$ and $\Delta S$ would correspond to the addition of the thermal changes of the individual PIII ↔ PII and the PII ↔ PI transitions, as confirmed by DSC.

Information on the thermal behaviour of [BF₄]$^-$ might be gained from ¹¹B NMR spectroscopy. Figure 13 present a series of ¹¹B spectra obtained by increasing T. Clear discontinuities in the spinning sidebands (SSBs) patterns (Figure 13a) and in the isotropic chemical shifts (Figure 13b) are observed and assigned to the transitions between phases occurring in these ranges of temperature. Fitting the SSBs allowed to estimate the parameters of the ¹¹B quadrupolar couplings (coupling constant $C_Q$ and asymmetry parameter $\mu_Q$) at each temperature.

<figure><img src="image_24.png" /><figcaption>**Figure 13.** ¹B MAS NMR spectra obtained at different temperatures ($v_0$ = 193 MHz, $v_{MAS}$= 5.5 kHz): (a) full spectral window showing the MAS satellite transitions sidebands; (b) zoom on the isotropic centre band.</figcaption></figure>

<figure><img src="image_25.png" /></figure>
ARTICLE

Journal Name

280 K). In this phase the [DBA]$^+$ cations are partially disordered on the ending methyl groups. Finally, phase I is stable for T > 280 K, with both [DBA]$^+$ cations and [BF$_4$]$^-$ anions heavily disordered.

These results highlight once again the importance of structural
disorder for enhancing barocaloric and thermal changes, which is
especially large in molecular and/or ionic compounds. Furthermore,
the results suggest that the combination of multiple phase transitions
that converge in a triple point can maximize thermal changes, which
will bring future triple-point materials under the spotlight of
barocaloric and cold-storage applications. In addition, it should be
noted that further efforts must be made in order to increase the
relatively low thermal conductivity of non-metallic thermomaterials
for practical applications. On the other hand, these results shows that,
despite apparent incompatibilities, it is possible to find solid-state
thermomaterials that can act as multipurpose compounds for active
barocaloric cooling and passive cold-storage. This open the door
towards the design of more efficient devices that could switch
between active (energy-consuming) and passive (non-energy
consuming) modes depending on the climate and/or surrounding
conditions, and/or that can enter from an active refrigeration into a
temporary cold-storage mode if an electric fail occur.

Conclusions

We have been able to obtain easily and by mild conditions the [DBA][BF₄] ionic compound, which is composed of lightweight atoms and has low density 1.25 g/cm³ (from single crystal diffraction at 100 K). This material undergoes three order-disorder transitions, exhibiting four phases in the temperature range 100-300 K, of which the lowest-temperature phase IV is newly reported here. The PIV ↔ PIII phase transition has a low entropy change ~0.4 J K⁻¹ kg⁻¹ at T₁₃ ~130 K and it is associated with partial disorder of the [BF₄]⁻ anions. On the other hand, the PIII ↔ PII phase transition exhibits a large entropy change ~205 J K⁻¹ kg⁻¹ at T₁₂ ~269 K and it is mainly due to order-disorder processes of the two-ending methyl group of the [DBA]⁺ cations. The PII ↔ PI phase transition has a moderate entropy change ~38 J K⁻¹ kg⁻¹ at T₁₁ ~282 K mainly due to disorder of all carbon atoms of the [DBA]⁺ cations. This material presents potential for cold-storage applications due to the fact that PIII ↔ PII and PII ↔ PI phase transitions show a total enthalpy change (~71.0 kJ kg⁻¹) similar to commercial materials, along with a moderate thermal hysteresis (3 K and 6 K for T₁₁ and T₁₂, respectively). Furthermore, in addition to the transition temperature (T₁₁ and T₁₂) and the large thermal changes (ΔS = 245 J K⁻¹ kg⁻¹), the [DBA][BF₄] compound exhibits a large sensitive to pressure (dT₂/dp ~29 K kbar⁻¹ and dTₘ/dp ~22.6 K kbar⁻¹) allowing the material to have excellent barocaloric properties, including a large pressure-driven reversible isothermal entropy change (ΔS<sub>rev</sub> > 200 J K⁻¹ kg⁻¹) and a large reversible adiabatic temperature change (ΔT<sub>rev</sub> ~17 K). Additionally, the material exhibits large barocaloric effects with a temperature span of 15 K, from 266 K to 281 K under the application of only 500 bar, which can be extended to 24K, from 266 K to 290 K, under the application of 1000 bar, resulting in an excellent material for barocaloric heating or refrigeration.

These results highlight that [DBA][BF₄] is not only a good candidate
as solid-solid cold-storage material, but also as a new barocaloric
material for cooling applications near 273 K, which make this new

<figure><img src="image_26.png" /><figcaption>**Figure 14.** ¹¹B quadrupolar parameters as a function of sample temperature (on heating).</figcaption></figure>

Their variations are plotted in Figure 14. The discontinuities observed for C_Q and μ_Q coincide with those observed by X-ray diffraction (see Figure 9). They are again related to the transitions between phases. For phase III, the increase in T leads to a slight and gradual increase in C_Q before a sudden increase at the transition to phase II. Such increase in C_Q might be explained by a loss in the degree of symmetry of the electric field gradient due: to a deformation of [BF_4]^⁻ and/or to a hindrance of the orientational motion of the anion. For phase II, when the T is increased, we observe a gradual decrease in C_Q and μ_Q that becomes sudden at the transition to phase I. The change in μ_Q is important (from 0.77 to 0.12 at the transition) and denotes either a strong reorganisation in the distribution of anions and cations or the appearance of new motions leading to a different time-averaged quadrupolar coupling tensor. This later argument bears consistency with the increase in motions detected by ¹H MAS NMR. In such a case, the small values of μ_Q points to an almost axial motion, compatible with the refined structure of phase I (see Figure 9). Therefore, these data further confirm a large disorder increase within the phase transitions, which would be the origin of the barocaloric properties.

In summary, all these structural studies have shown that the four
phases -and the corresponding three phase transitions- of [DBA][BF₄]
are related to an increase in the disorder of the ions (see Figure 9). At
low temperature (T < 130 K), phase IV is the most stable phase and
both cations and anions are fully ordered. In phase III, observed in the
temperature range 130 K to 260 K, the [DBA]<sup>+</sup> cations remain ordered
with a zig-zag conformation, but the [BF<sub>4</sub>]<sup>-</sup> anions are slightly
disordered. Phase II is stable in a small range of temperatures (260 -

<figure><img src="image_27.png" /></figure>

Journal of Materials Chemistry A Accepted Manuscript
compound a very versatile multipurpose material for food preservation, either through passive cold-storage or active barocaloric refrigeration, opening the door to develop new multipurpose solid thermomaterials.

Finally, this work emphasises that relating the thermal properties to the structural transitions of the [DBA][BF₄], thanks to the combination of various characterization techniques, is essential to understand the structure-properties relationships that are necessary to design new thermomaterials.

## Conflicts of interest

There are no conflicts to declare.

## Authors contributions

J.G.-B. and J.M.B.-G. contributed equally to this work and were involved in most of experimental work. S.C.-G., M.S.-A. and M.A.S.-R. conceived the project. R.J.C.D, and A.E.P. were involved in VT-SPXRD and high pressure DSC experiments. M.S.-A. and I.D.-F. were involved in VT-SPXRD experiments. S.C.-G., M.S.-A. and A.L.L.-S. were involved SC-XRD experiments. J.L.-B. and R.A. were involved in low-pressure DSC and Cp calculations. A.G.-F. and U.B.C. were involved in low-temperature DSC and VT-XRD experiments. B.A. performed solid-state NMR. All authors were involved in the results discussion and in the manuscript writing and revision. M.S.-A. and M.A.S.-R. directed and supervised this project.

## Acknowledgements

The authors thank financial support from the project PID2021-122532OB-100 funded by MCIN/AEI/10.13039/501100011033/ and by FEDER Una manera de hacer Europa, the project PDC2021-121076-100 funded by MCIN/AEI/10.13039/501100011033 and by the European Union Next GenerationEU/PRTR and the project ED431C 2022/39 funded by Xunta de Galicia. This publication is part of the grant RYC2021-033040-I, funded by MCIN/AEI/10.13039/501100011033 and from European Union «NextGenerationEU»/PRTR», and granted to J.M.B.G. This researcher is also grateful for the support received by UDC-Inditex InTalent Programme. J.G.-B. acknowledge Xunta de Galicia for Predoctoral Fellowship. A. G.-F. and U.B.C. acknowledge the financial support of the Göran Gustafsson's foundation. The R.J.C.D. and A.E.P. would like to thank the financial support of EPSRC throughout provision of the funding (EP/S03577X/1). The authors thank Diamond Light Source for the award of beamtime (CY30448-1 on I11), and Eamonn Connolly for his assistance on the experiment. The authors also thank Xavier Moya from the Department of Materials Science at the University of Cambridge, for granting access to his lab for preliminary VP-DSC measurements on still non-stable and wet samples and the use of the analytical facilities available at the Research Infrastructures Area of the University of Santiago de Compostela.

## Notes and references

The Role of Refrigeration in Worldwide nutrition (2020), 6th Informatory Note on Refrigeration and Food, https://iifiir.org/en/fridoc/the-role-of-refrigeration-in-worldwide-nutrition-2020-142029, (accessed 19 January 2023).

IIF-IIR, D. Coulomb, J. Dupont and A.- Pichard, The Role of Refrigeration in the Global Economy (2015), 29th Note on Refrigeration Technologies, https://iifiir.org/en/fridoc/the-role-of-refrigeration-in-the-global-economy-2015-138763, (accessed 13 February 2023).

IIF-IIR, V. Morlet, D. Coulomb and J. Dupont, The impact of the refrigeration sector on climate change, 35th Informatory Note on refrigeration technologies, https://iifiir.org/en/fridoc/the-impact-of-the-refrigeration-sector-on-climate-change-141135, (accessed 13 February 2023).

S. J. James and C. James, *Chilling and Freezing*, Elsevier Inc., 2014.

B. H. Lado and A. E. Yousef, *Microbes Infect.*, 2002, **4**, 433–440.

E. Oró, A. de Gracia, A. Castell, M. M. Farid and L. F. Cabeza, *Appl. Energy*, 2012, **99**, 513–533.

G. Li, Y. Hwang and R. Radermacher, *Int. J. Refrig.*, 2012, **35**, 2053–2077.

G. Li, Y. Hwang, R. Radermacher and H. H. Chun, *Energy*, 2013, **51**, 1–17.

X. Q. Zhai, X. L. Wang, T. Wang and R. Z. Wang, *Renew. Sustain. Energy Rev.*, 2013, **22**, 108–120.

W. Sun, W. Di Liu, L. Li, D. Z. Wang, L. C. Yin, M. Li, X. L. Shi, Q. Liu and Z. G. Chen, *Energy Convers. Manag.*, 2023, **290**, 117222.

W. Y. Chen, X. L. Shi, J. Zou and Z. G. Chen, *Small Methods*, 2022, **6**, 2101235.

X. Moya, S. Kar-Narayan and N. D. Mathur, *Nat. Mater.*, 2014, **13**, 439–450.

P. Lloveras and J.-L. Tamarit, *MRS Energy Sustain.*, 2021, **8**, 3–15.

D. Boldrin, *Appl. Phys. Lett.*, 2021, **118**, 170502.

A. Fallahi, G. Guldemonts, M. Tao, S. Granados-Focil and S. Van Dessel, *Appl. Therm. Eng.*, 2017, **127**, 1427–1441.

C. R. Raj, S. Suresh, R. R. Bhavsar and V. K. Singh, *J. Therm. Anal. Calorim.*, 2020, **139**, 3023–3049.

L. Mañosa, D. González-Alonso, A. Planes, E. Bonnot, M. Barrio, J. L. Tamarit, S. Aksoy and M. Acet, *Nat. Mater.*, 2010, **9**, 478–481.

E. Stern-Taulats, A. Planes, P. Lloveras, M. Barrio, J. L. Tamarit, S. Pramanick, S. Majumdar, S. Yüce, B. Emre, C. Frontera and L. Mañosa, *Acta Mater.*, 2015, **96**, 324–332.

J. Lin, P. Tong, X. Zhang, Z. Wang, Z. Zhang, B. Li, G. Zhong, J. Chen, Y. Wu, H. Lu, L. He, B. Bai, L. Ling, W. Song, Z. Zhang and Y. Sun, *Mater. Horizons*, 2020, **7**, 2690–2695.

S. Yuce, M. Barrio, B. Emre, E. Stern-Taulats, A. Planes, J. L. Tamarit, Y. Mudryk, K. A. Gschneidner, V. K. Pecharsky and L. Mañosa, *Appl. Phys. Lett.*, 2012, **101**, 071906.
ARTICLE

Journal Name

21 L. Mañosa, D. González-Alonso, A. Planes, M. Barrio, J. L. Tamarit, I. S. Titov, M. Acet, A. Bhattacharyya and S. Majumdar, *Nat. Commun.*, 2011, **2**, 1–5.

22 L. Mañosa, E. Stern-Taulats, A. Planes, P. Lloveras, M. Barrio, J. L. Tamarit, B. Emre, S. Yüce, S. Fabbrici and F. Albertini, *Phys. Status Solidi Basic Res.*, 2014, **251**, 2114–2119.

23 X. J. He, K. Xu, S. X. Wei, Y. L. Zhang, Z. Li and C. Jing, *J. Mater. Sci.*, 2017, **52**, 2915–2923.

24 Z. Wei, Y. Shen, Z. Zhang, J. Guo, B. Li, E. Liu, Z. Zhang and J. Liu, *APL Mater.*, 2020, **8**, 051101.

25 X. He, S. Wei, Y. Kang, Y. Zhang, Y. Cao, K. Xu, Z. Li and C. Jing, *Scr. Mater.*, 2018, **145**, 58–61.

26 X. He, Y. Kang, S. Wei, Y. Zhang, Y. Cao, K. Xu, Z. Li, C. Jing and Z. Li, *J. Alloys Compd.*, 2018, **741**, 821–825.

27 P. Lloveras, T. Samanta, M. Barrio, I. Dubenko, N. Ali, J. L. Tamarit and S. Stadler, *APL Mater.*, 2019, **7**, 0–9.

28 A. Aznar, A. Gràcia-Condal, A. Planes, P. Lloveras, M. Barrio, J. L. Tamarit, W. Xiong, D. Cong, C. Popescu and L. Mañosa, *Phys. Rev. Mater.*, 2019, **3**, 1–7.

29 H. Liu, Z. Li, Y. Zhang, Z. Ni, K. Xu and Y. Liu, *Scr. Mater.*, 2020, **177**, 1–5.

30 D. Matsunami, A. Fujita, K. Takenaka and M. Kano, *Nat. Mater.*, 2015, **14**, 73–78.

31 R. R. Wu, L. F. Bao, F. X. Hu, H. Wu, Q. Z. Huang, J. Wang, X. L. Dong, G. N. Li, J. R. Sun, F. R. Shen, T. Y. Zhao, X. Q. Zheng, L. C. Wang, Y. Liu, W. L. Zuo, Y. Y. Zhao, M. Zhang, X. C. Wang, C. Q. Jin, G. H. Rao, X. F. Han and B. G. Shen, *Sci. Rep.*, 2015, **5**, 1–11.

32 D. Boldrin, E. Mendive-Tapia, J. Zemen, J. B. Staunton, T. Hansen, A. Aznar, J. L. Tamarit, M. Barrio, P. Lloveras, J. Kim, X. Moya and L. F. Cohen, *Phys. Rev. X*, 2018, **8**, 41035.

33 A. Aznar, P. Lloveras, J. Y. Kim, E. Stern-Taulats, M. Barrio, J. L. Tamarit, C. F. Sánchez-Valdés, J. L. Sánchez Llamazares, N. D. Mathur and X. Moya, *Adv. Mater.*, 2019, **31**, 1–6.

34 E. Stern-Taulats, A. Planes, P. Lloveras, M. Barrio, J. L. Tamarit, S. Pramanick, S. Majumdar, C. Frontera and L. Mañosa, *Phys. Rev. B - Condens. Matter Mater. Phys.*, 2014, **89**, 1–8.

35 C. Yu, J. Huang, J. Qi, P. Liu, D. Li, T. Yang, Z. Zhang and B. Li, *APL Mater.*, 2022, **10**, 011109.

36 P. Lloveras, E. Stern-Taulats, M. Barrio, J. L. Tamarit, S. Crossley, W. Li, V. Pomjakushin, A. Planes, L. Mañosa, N. D. Mathur and X. Moya, *Nat. Commun.*, 2015, **6**, 8801.

37 Q. Ren, J. Qi, D. Yu, Z. Zhang, R. Song, W. Song, B. Yuan, T. Wang, W. Ren, Z. Zhang, X. Tong and B. Li, *Nat. Commun.*, 2022, **13**, 1–9.

38 M. V. Gorev, E. V. Bogdanov, I. N. Flerov, A. G. Kocharova and N. M. Laptash, *Phys. Solid State*, 2010, **52**, 167–175.

39 M. V. Gorev, E. V. Bogdanov, I. N. Flerov, V. N. Voronov and N. M. Laptash, *Ferroelectrics*, 2010, **397**, 76–80.

40 I. N. Flerov, A. V. Kartashev, M. V. Gorev, E. V. Bogdanov, S. V. Mel'Nikova, M. S. Molokeev, E. I. Pogoreltsev and N. M. Laptash, *J. Fluor. Chem.*, 2016, **183**, 1–9.

M. V. Gorev, E. A. Mikhaleva, I. N. Flerov and E. V. Bogdanov, *J. Alloys Compd.*, 2019, **806**, 1047–1051.

Z. Zhang, X. Jiang, T. Hattori, X. Xu, M. Li, C. Yu, Z. Zhang, D. Yu, R. Mole, S. I. Yano, J. Chen, L. He, C. W. Wang, H. Wang, B. Li and Z. Zhang, *Mater. Horizons*, 2023, **10**, 977–982.

X. Xu, F. Li, C. Niu and H. Wang, *Appl. Phys. Lett.*, 2023, **122**, 043901.

A. Aznar, P. Lloveras, M. Romanini, M. Barrio, J. L. Tamarit, C. Cazorla, D. Errandonea, N. D. Mathur, A. Planes, X. Moya and L. Mañosa, *Nat. Commun.*, 2017, **8**, 1851.

J. Min, A. K. Sagotra and C. Cazorla, *Phys. Rev. Mater.*, 2020, **4**, 1–10.

M. Romanini, Y. Wang, G. Kübra, G. Ornelas, P. Lloveras, Y. Zhang, W. Zheng, M. Barrio, A. Aznar, A. Gràcia-Condal, B. Emre, O. Atakol, C. Popescu, H. Zhang, Y. Long, L. Balicas, J. L. Tamarit, A. Planes, M. Shatruk and L. Mañosa, *Adv. Mater.*, 2021, **33**, 2008076.

S. P. Vallone, A. N. Tantillo, A. M. dos Santos, J. J. Molaison, R. Kulmaczewski, A. Chapoy, P. Ahmadi, M. A. Halcrow and K. G. Sandeman, *Adv. Mater.*, 2019, **31**, 1–7.

P. J. Von Ranke, B. P. Alho, P. H. S. Da Silva, R. M. Ribas, E. P. Nobrega, V. S. R. De Sousa, M. V. Colaço, L. F. Marques, M. S. Reis, F. M. Scaldini, L. B. L. Escobar and P. O. Ribeiro, *J. Appl. Phys.*, 2020, **127**, 165104.

J. Seo, J. D. Braun, V. M. Dev and J. A. Mason, *J. Am. Chem. Soc.*, 2022, **144**, 6493–6503.

P. J. Von Ranke, B. P. Alho, R. M. Ribas, E. P. Nobrega, A. Caldas, V. S. R. De Sousa, M. V. Colaço, L. F. Marques, D. L. Rocco and P. O. Ribeiro, *Phys. Rev. B*, 2018, **98**, 2–6.

P. J. von Ranke, B. P. Alho and P. O. Ribeiro, *J. Alloys Compd.*, 2018, **749**, 556–560.

P. J. von Ranke, B. P. Alho, E. P. Nobrega, A. Caldas, V. S. R. de Sousa, M. V. Colaço, L. F. Marques, G. M. Rocha, D. L. Rocco and P. O. Ribeiro, *J. Magn. Magn. Mater.*, 2019, **489**, 165421.

C. M. Miliante, A. M. Christmann, R. P. Soares, J. R. Bocca, C. S. Alves, A. M. G. Carvalho and A. R. Muniz, *J. Mater. Chem. A*, 2022, **10**, 8344–8355.

J. M. Bermúdez-García, M. Sánchez-Andújar, S. Castro-García, J. López-Beceiro, R. Artiaga and M. A. Señarís-Rodríguez, *Nat. Commun.*, 2017, **8**, 15715.

J. M. Bermúdez-García, S. Yáñez-Vilar, A. García-Fernández, M. Sánchez-Andújar, S. Castro-García, J. López-Beceiro, R. Artiaga, M. Dilshad, X. Moya and M. A. Señarís-Rodríguez, *J. Mater. Chem. C*, 2018, **6**, 9867–9874.

M. Szfrański, W. J. Wei, Z. M. Wang, W. Li and A. Katrusiak, *APL Mater.*, 2018, **6**, 100701.

J. Salgado-Beceiro, A. Nonato, R. X. Silva, A. García-Fernández, M. Sánchez-Andújar, S. Castro-García, E. Stern-Taulats, M. A. Señarís-Rodríguez, X. Moya and J. M. Bermúdez-García, *Mater. Adv.*, 2020, **1**, 3167–3170.

J. Salgado-Beceiro, J. M. Bermúdez-García, E. Stern-Taulats, J. García-Ben, S. Castro-García, M. Sánchez-Andújar, X. Moya and M. A. Señarís-Rodríguez,

View Article Online

Journal of Materials Chemistry A Accepted Manuscript

<figure><img src="image_28.png" /></figure>
ChemRxiv, 2021, 17–19.

59 J. Li, M. Barrio, D. J. Dunstan, R. Dixey, X. Lou, J. L. Tamarit, A. E. Phillips and P. Lloveras, *Adv. Funct. Mater.*, 2021, **31**, 1–8.

60 J. Seo, R. D. McGillicuddy, A. H. Slavney, S. Zhang, R. Ukani, A. A. Yakovenko, S. L. Zheng and J. A. Mason, *Nat. Commun.*, 2022, **13**, 1–15.

61 Y. Gao, H. Liu, F. Hu, H. Song, H. Zhang, J. Hao, X. Liu, Z. Yu, F. Shen, Y. Wang, H. Zhou, B. Wang, Z. Tian, Y. Lin, C. Zhang, Z. Yin, J. Wang, Y. Chen, Y. Li, Y. Song, Y. Shi, T. Zhao, J. Sun, Q. Huang and B. Shen, *NPG Asia Mater.*, 2022, **14**, 34.

62 P. Lloveras, A. Aznar, M. Barrio, P. Negrier, C. Popescu, A. Planes, L. Mañosa, E. Stern-Taulats, A. Avramenko, N. D. Mathur, X. Moya and J. L. Tamarit, *Nat. Commun.*, 2019, **10**, 1–7.

63 A. Aznar, P. Lloveras, M. M. Barrio, P. Negrier, A. Planes, L. L. Mañosa, N. D. Mathur, X. Moya and J. L. L. Tamarit, *J. Mater. Chem. A*, 2020, **8**, 639–647.

64 J. Li, D. Dunstan, X. Lou, A. Planes, L. Mañosa, M. Barrio, J. L. Tamarit and P. Lloveras, *J. Mater. Chem. A*, 2020, **8**, 20354–20362.

65 A. Aznar, P. Negrier, A. Planes, L. Mañosa, E. Stern-Taulats, X. Moya, M. Barrio, J. L. Tamarit and P. Lloveras, *Appl. Mater. Today*, 2021, **23**, 101023.

66 A. Salvatori, P. Negrier, A. Aznar, M. Barrio, J. L. Tamarit and P. Lloveras, 2022, **10**, 111117.

67 K. Zhang, R. Song, J. Qi, Z. Zhang, Z. Zhang, C. Yu, K. Li, Z. Zhang and B. Li, *Adv. Funct. Mater.*, 2022, **32**, 2112622.

68 F. B. Li, M. Li, X. Xu, Z. C. Yang, H. Xu, C. K. Jia, K. Li, J. He, B. Li and H. Wang, *Nat. Commun.*, 2020, **11**, 1–8.

69 F. Li, C. Niu, X. Xu, M. Li and H. Wang, *Appl. Phys. Lett.*, 2022, **121**, 223902.

70 W. Imamura, É. O. Usuda, L. S. Paixão, N. M. Bom, A. M. Gomes and A. M. G. Carvalho, *Chinese J. Polym. Sci.* (English Ed.), 2020, **38**, 999–1005.

71 S. Patel, A. Chauhan, R. Vaish and P. Thomas, *Appl. Phys. Lett.*, 2016, **108**, 072903.

72 A. M. G. Carvalho, W. Imamura, E. O. Usuda and N. M. Bom, *Eur. Polym. J.*, 2018, **99**, 212–221.

73 J. García-Ben, J. López-Beceiro, R. Artiaga, J. Salgado-Beceiro, I. Delgado-Ferreiro, Y. V. Kolen'ko, S. Castro-García, M. A. Señarís-Rodríguez, M. Sánchez-Andújar and J. M. Bermúdez-García, *Chem. Mater.*, 2022, **34**, 3323–3332.

74 J. Font, J. Muntasell, J. Navarro, J.-L. Tamarit and J. Lloveras, *Sol. Energy Mater.*, 1987, **15**, 299–310.

75 M. Barrio, J. Font, J. Muntasell, J. Navarro and J. L. Tamarit, *Sol. Energy Mater.*, 1988, **18**, 109–115.

76 B. Li, Y. Kawakita, S. Ohira-Kawamura, T. Sugahara, H. Wang, J. Wang, Y. Chen, S. I. Kawaguchi, S. Kawaguchi, K. Ohara, K. Li, D. Yu, R. Mole, T. Hattori, T. Kikuchi, S. ichiro Yano, Z. Zhang, Z. Zhang, W. Ren, S. Lin, O. Sakata, K. Nakajima and Z. Zhang, *Nature*, 2019, **567**, 506–510.

77 M. Masche, L. Ianniciello, J. Tušek and K. Engelbrecht, *Int. J. Refrig.*, 2021, **121**, 302–312.

78 T. Hess, L. M. Maier, N. Bachmann, P. Corhan, O. Schäfer, Welsen, J. Wöllenstein and K. Bartholome, *J. Am. Chem. Soc.*, 2020, **142**, 1039–1046.

79 S. K. Mohamed and D. K. Padma, *Indian J. Chem. A Inorganic, Phys. Theor. Anal. Chem.*, 1988, **27**, 759–763.

80 G. Zabinska, P. Ferloni and M. Sanesi, *Thermochim. Acta*, 1988, **137**, 39–49.

81 A. Bielecki and D. P. Burum, *J. Magn. Reson. Ser. A*, 1995, **116**, 215–220.

82 G. Ferrer, C. Barreneche, A. Solé, I. Martorell and L. F. Cabeza, *J. Energy Storage*, 2017, **11**, 1–6.

83 ASTM E1269-11, *Standard test method for determining specific heat capacity by differential scanning calorimetry*, 2018.

84 B. H. Toby and R. B. Von Dreele, *J. Appl. Crystallogr.*, 2013, **46**, 544–549.

85 G. M. Sheldrick, *Acta Crystallogr. Sect. A Found. Crystallogr.*, 2015, **71**, 3–8.

86 G. M. Sheldrick, *Acta Crystallogr. Sect. C Struct. Chem.*, 2015, **71**, 3–8.

87 S. D. Sharma and K. Sagara, *Int. J. Green Energy*, 2005, **2**, 1–56.

88 M. V. Gorev, E. V. Bogdanov and I. N. Flerov, *J. Phys. D. Appl. Phys.*, 2017, **50**, 384002.

89 X. Moya and N. D. Mathur, *Science (80-.)*, 2020, **803**, 797–803.

90 J. Lin, P. Tong, K. Zhang, K. Tao, W. Lu, X. Wang, X. Zhang, W. Song and Y. Sun, *Nat. Commun.*, 2022, **13**, 596.

91 L. Ricard, R. Cavagnat and M. Rey-Lafon, *J. Phys. Chem.*, 1985, **89**, 4887–4894.

92 Y. Abid, M. Kamoun, A. Daoud and F. Romain, *J. Raman Spectrosc.*, 1990, **21**, 709–716.

<figure><img src="image_29.png" /></figure>