Di‑Alkyl Adipates as New Phase Change Material for Low Temperature Energy Storage

This work is a contribution to the thermal characterization of a selected binary sys-tem of two di-n -alkyl adipates that can be used as phase change material for thermal energy storage at low temperatures. The construction of the solid–liquid phase diagram using differential scanning calorimetry (DSC), complemented with Raman Spectroscopy studies for the system composed by diethyl and dibutyl adipates is presented. The solidus and liquidus equilibrium temperatures were determined by DSC for the pure components and 30 binary mixtures at selected molar compositions were used to construct the corresponding solid–liquid phase diagram. The binary system of diethyl and dibutyl adipates presents eutectic behaviour at low temperatures. The eutectic temperature was found at 240.46 K, and the eutectic composition was determined to occur at the molar fraction x dibutyl = 0.46. Additionally, the system shows a polymorphic transition, characteristic of dibutyl adipate, occurring at ca. 238 K, confirmed by optical microscopy. To the best of our knowledge, no reference to the phase diagram of the present system could be found in the literature. Raman spectroscopy was essential to complement the construction of the phase equilibrium diagram, enabling the identification of the solid and liquid phases of the system. Finally, the liquidus curve of the phase diagram was also successfully predicted using a suitable fitting equation, being the root mean square deviation of the data from the correlation equal to 0.54 K. In addition, this fitting operation enabled a correct prediction of the eutectic composition of the system.


Introduction
The increasing worldwide need for energy requires urgent and effective solutions.It is widely recognized that the supply of renewable energies is principally intermittent, which often leads to a disturbing gap between energy availability and its consumption.Therefore, energy storage becomes a very significant tool for increasing the efficiency in the use of the produced renewable energy.Energy Storage is thus considered to be a fundamental means to enable the transition to alternative cleaner and more efficient use of energy.As a consequence, it has received worldwide attention and increasing study interest [1].
Thermal energy storage (TES) is one of the different forms of energy storage and its use has been increasing in recent decades and one growing technology makes use of the latent heat storage, which relies on the energy involved in the phase transition, often between the solid and liquid phases of a material.Latent heat storage (LHS) became attractive as a result of some specific characteristics such as the nearly isothermal nature of the phase-change process, compactness and the high latent heat inherently involved [2][3][4][5][6][7][8][9].
Phase change materials (PCMs), are materials specifically used for latent heat storage applications [10,11].These materials are characterized for being able to absorb and release heat from the surroundings by changing phase, thus enabling the storage of a substantial amount of thermal energy [10,12].In fact, the use of PCMs for latent heat storage provides a high energy storage density due to the large amount of heat exchanged during the phase-change process.Furthermore, this exchange occurs with a small variation of the temperature of the PCM, as the temperature must remain practically constant while the two phases are in equilibrium.Additionally, the volume of the PCM does not change significantly during the solid-liquid phase-change, which is an interesting characteristic for their industrial applications; one advantage is a more efficient energy transport to a heat transfer fluid (HTF).Consequently, latent heat storage using PCMs is one of the most effective systems to store thermal energy [4], having attracted significant attention from several authors [11][12][13][14][15][16][17].Additionally, the use of these materials for latent heat storage has taken a very important role, especially for active cooling systems, in particular, eutectic systems that can be customized for each application, allowing significant energy savings, being relevant to the economic and environmental sustainability [18].
Phase equilibrium studies aimed to obtain phase diagrams which are important to characterize the behaviour of PCM systems, including the desired melting temperatures, revealing their aptitude for an industrial application [19].Differential scanning calorimetry (DSC) is one of the most effective and used techniques to determine whether binary mixtures form solid solutions or eutectic systems [1].In fact, it is an important tool to construct solid-liquid phase diagrams for binary organic and inorganic systems [20].
In recent years, alkanes and paraffins have been under massive research because of their high latent heat enthalpy and chemical stability.However, in general, the high cost of paraffins when compared to other organic compounds restricts their use in a wide range of applications [21].To solve this problem, various non-paraffinic-based compounds, such as esters, fatty alcohols and fatty acids have also been developed as PCMs due to their advantageous characteristics such as their availability, low cost, and renewability.Despite the numerous studies available in the literature with organic PCMs, they are not focused on low temperature energy storage applications, and this can be attributed to the lack of information concerning essential characteristics as the fusion temperature range and their thermophysical properties, such as the thermal conductivity.Thus, the search for new suitable organic PCMs for cold TES systems is becoming crucial to fulfil the needs for the desired applications [19,21].Regarding this, n-alkyl adipates, and other esters, have been recently started to be studied and proved to be viable PCMs for TES applications at low temperatures [21][22][23].Additionally, they also present suitable properties for industrial applications such as low cost, low toxicity and chemical stability, and are considered to be environmentally friendly [24,25].
This work aims to the characterization of a selected system of two di-n-alkyladipates that can be used for energy storage at low temperatures, starting by the construction of solid-liquid phase diagrams using DSC and Raman spectroscopy studies.In addition, to the best of our knowledge, no research studies exploring binary system of diethyl and dibutyl adipates, for cool thermal storage applications, has been reported in the literature.

Materials
Diethyl adipate was supplied by Sigma Aldrich, Lot number STBD3299V, with a purity > 99.9 %, according to the Certificate of Analysis.Dibutyl adipate was supplied by TCI Chemicals, Lot number SG2JN-ML, with a purity of 99.9 %, according to the corresponding Certificate of Analysis.
The water content of the materials used in this work were measured using a Karl-Fisher 831 KF Coulometer from Metrohm.The values of water content were 617.4 ppm for diethyl adipate and 6.8 ppm for dibutyl adipate.

Differential Scanning Calorimetry (DSC)
The calorimetric measurements were performed with a 2920 MDSC system from TA Instruments Inc.The sample masses of 3.0-10.0mg were sealed in air inside aluminium pans, weighed with a precision of ± 0.1 µg in a Mettler UMT2 ultramicro balance and analysed at a scanning rating of β = 5 K•min -1 .Helium (Air Liquide N55), at a flow rate of 30 cm 3 min −1 was used as purging gas.The baseline was corrected by scanning the temperature range of the experiments with an empty pan.
The temperature and heat flow scales of the instrument were calibrated at different heating rates, based on the onsets of the fusion peaks of several standards.Details of the calibration are described and presented elsewhere [26].The software used was TA Universal Analysis 2000 v4.0.

Hot Stage Microscopy (HSM)
Polarized optical microscopy observations were performed on an Olympus BX51 optical microscope.The temperature changes and stabilization were set by a Linkam LTS360 liquid nitrogen cooled cryostage and was measured with a Pt resistance thermometer.The liquid samples were first placed on the glass plate, covered by a second glass and cooled down at 5 K•min −1 .After solidification, the microstructure of the sample was monitored with an Olympus C5060 wide zoom camera for picture and/or movie record.Images were recorded with 250× magnification in the temperature range from 183.15 K to 293.15 K at a heating rate of 5 K•min −1 .

Raman Spectroscopy
The Raman spectra were obtained using a micro-Raman Horiba LabRam HR Evolution system, with excitation provided by a 532 nm solid-state laser, using an approximate power of 50 mW at the sample.A 50× long working distance objective was used, with a laser spot diameter of around 1 µm in the sample.The equipment was calibrated using a Si crystal (reference band at 520.5 cm -1 ).The final spectra were the average of 5 spectra accumulations obtained with an acquisition-time of 5 s and resolution of 0.5 cm -1 .A hot stage THMS 600 (Linkam Scientific Instruments), controlled by a T95-PE Linkpad controlling unit, was used for the temperature-variation Raman experiments.

Differential Scanning Calorimetry (DSC)
DSC heating curves of the pure compounds, diethyl adipate and dibutyl adipate and some of the most significant binary mixtures are presented in Fig. 1.The relevant temperatures and enthalpies of fusion were obtained using the average of two corresponding values for each cycle and are given in Table 1; the supplementary material contains the remaining DSC curves in Fig. S1, and the corresponding temperature and enthalpy of fusion values in Table S1.
The curves shown in Fig. 1 evidence two different types of solid-liquid transition.In the case of pure compounds and the eutectic mixture, only a single peak corresponding to the fusion is detectable while for the mixtures with molar fraction x dibutyl = 0.80 and x dibutyl = 0.20 it is possible to identify two distinct peaks; also, for the mixture with x dibutyl = 0.20 the peak assigned to the polymorphism of dibutyl adipate is detectable.Nonetheless, the temperature of the polymorphic transition is always the same and does not affect the solid-liquid phase transition for any of the studied mixtures, including the eutectic mixture which is the most interesting one in terms of TES application.These results enable to rationalise the different shape of the curves according to the molar fraction of diethyl adipate.
Concerning the pure components, it was observed that diethyl adipate melts at T max = 253.55K with an enthalpy of fusion, Δ fus H = 102.8J•g −1 and dibutyl adipate melts at 251.17 K with Δ fus H = 153.1 J•g −1 .A eutectic mixture is found with molar fraction x dibutyl = 0.46 with a eutectic temperature T max = 240.46K with an enthalpy of fusion, Δ fus H = 131.7 J•g −1 .Furthermore, it has been observed that dibutyl adipate presents polymorphism, being the polymorph  These preliminary findings obtained from isolated DSC experiments, have been complemented by Hot Stage Microscopy and needed confirmation using another sort of instrumentation and studies, like low temperature Raman spectroscopy.

Hot Stage Microscopy
Hot stage microscopy enabled to verify the results obtained using DSC.By employing light polarization, it was possible to observe not only the solid-liquid phase transition, but also the polymorphic transition.In Fig. 2, it is possible to visualize two different cases: (a) the polymorphism of dibutyl adipate for the eutectic sample and, (b) the phase transition for the x dibutyl = 0.10 mixture.
In the first case, Fig. 2a, the mixture x dibutyl = 0.1 reveals that a partial phase transition with the appearance of some drops of liquid first occurs, and that this mixture of solid and liquid phases continues to evolve until there is only the liquid phase.In Fig. 2b, the eutectic mixture is shown, where not only it is possible to see the polymorphic transition of dibutyl adipate by the change of form, with concomitant colour modification of the solid, but also the sudden phase transition similar to a pure compound, a typical signature of a eutectic transition.

Raman Spectroscopy
Raman spectra of pure diethyl adipate and dibutyl adipate for both solid and liquid phases, of the eutectic mixture and of one binary mixture having x dibutyl = 0.20 for different temperatures are presented in Figs. 3, 4, 5 and 6, respectively; being the remaining Raman spectra demonstrated in Figs.S2, S3 and S4 in the supplementary information.Additionally, the Raman equipment enabled to obtain images of the system during the experiments, which is particularly significant to monitor the system in real time and corroborate visually the results obtained using DSC.These For the x dibutyl = 0.20 mixture at 193.15 K, the marker bands of diethyl adipate (e.g., 863 cm −1 ) and dibutyl adipate polymorph 2 (e.g., 843 cm −1 and 932 cm −1 ) are clearly identifiable.At 223.15 K, the dibutyl adipate polymorph 2 disappears (as identified by the circle in Fig. 6) and a new band at 937 cm −1 appears, which might indicate the evolution of dibutyl adipate polymorph 2 into polymorph 1.At 238.15 K, the marker bands of dibutyl adipate polymorph 1 (e.g., 295 cm −1 , 537 cm −1 , 846 cm −1 and 942 cm −1 ) are identified, meaning that the polymorph (2 → 1) transition is fully completed at this temperature.At 243.15 K, the dibutyl adipate polymorph 1 marker bands are no longer present, corresponding to the complete melting of this compound, and the marker bands of diethyl adipate start to broaden, corresponding to the beginning of the melting process of this solid, which is completed at 248.15 K.

Solid-Liquid Binary Phase Diagram
The DSC thermograms and the complementary studies carried out using Hot Stage Microscopy and, especially, Low Temperature Raman Spectroscopy results have led to the construction of a solid-liquid phase diagram for the system (dibutyl adipate + diethyl adipate) as shown in Fig. 9.This phase diagram was obtained using the experimental results for 30 different molar fractions of the binary mixture and for the two pure components.It is possible to identify two The freezing-point depression curve can, in principle, be used to characterise the liquidus line of a solid-liquid binary phase diagram of a non-isomorphous eutectic system [19].It can be described by Eq. 1, where x is the molar fraction, T is the liquidus temperature for molar fraction x , T Mi and L Mi are the temperature and enthalpy of melting of the pure compound i, respectively, and R is the gas constant.This equation results from assuming that: the two components of an ideal liquid solution do not form solid solutions with each other; the ideal solution is in equilibrium with one or the other of the pure solid phases and that L is independent of the temperature [27].This same equation has been previously used by other authors, e.g., Campanella et al. [20], Cortesão et al. [28], and Han et al. [29], for that purpose.
In the present work, the liquidus line was obtained by a fitting procedure, making use of all the experimental points obtained in the DSC studies.These 32 experimental points cover the whole range of compositions, including the fusion points for the pure components, and 30 binary mixtures.For that purpose, it was decided that an equal statistical weight should be attributed to all the data, including the fusion temperatures of the pure components.Therefore, Eq. 2 has been applied to define the liquidus line of the present system, as shown in Fig. 9. Equation 2 is derived from Eq. 1, letting both L Mi and T Mi lose their precise original physical significance and be included in the fitting parameters, a and b, and can be estimated by Eqs. 3 and 4, for comparison.
The liquidus line obtained using the fitting Eq. 2 shows an excellent agreement with the 32 experimental data points obtained by differential scanning calorimetry; the coordinates of the eutectic point being well predicted by the intersection of Eq. 2, when applied to each side of the diagram, as depicted in Fig. 9.
The experimental solidus and liquidus temperature points were obtained from the maximum of the first and second endothermic peaks in the DSC curves, respectively.Figure 10 shows the deviations of the thermal analysis experimental data from the correlation Eq. 2, using the parameters in Table 2.The maximum absolute deviation of all data are smaller than ± 1.5 K.
(1)  Additionally, the eutectic composition and temperature were calculated by the intersection of the two fitted liquidus lines for each side of the binary phase diagram.According to the fitting Eq. 2, the eutectic composition, x dibutyl , is 0.44 and the eutectic temperature is T fus = 239.78K, which corresponds to an absolute deviation of 0.68 K compared to the experimental result.The root mean square deviation of the experimental data points to the fitting equation for the liquidus line, hereby designated as σ, was determined as 0.54 K. Therefore, the uncertainty is well within a 0.95 confidence level.

Conclusions
Organic PCMs have been under massive research due to their high latent heat enthalpy and chemical stability.In particular non-paraffin-based compounds have, recently, shown good potential as PCMs for TES applications supported by the availability, low cost and renewability.Although, these studies are not focused on low temperature energy storage applications and, thus, the search for new suitable organic PCMs for cold TES systems is becoming crucial to fulfil the needs for the desired applications.
In the present article, a binary system of two di-n-alkyl adipates, in particular, diethyl and dibutyl adipates, is characterized in order to be used for low temperature energy storage applications.DSC tests performed with 30 different mole fractions of the binary mixture, complemented with studies achieved by Raman spectroscopy with temperature-variation allowed the detailed construction of the solid-liquid phase diagram for this system.Additionally, the liquidus curve of the phase diagram was also successfully described using a suitable fitting equation.The studied binary system reveals eutectic behaviour at low temperatures, being the eutectic temperature found at 240.46 K with a molar fraction x dibutyl = 0.46 and an enthalpy of fusion Δ fus H = 131.7 J•g −1 .The characteristic polymorphic transition of dibutyl adipate is also observed in the binary mixtures at around 238 K with an enthalpy value of Δ trans H = 2.50 J•g −1 .The optical image of the polymorphic transition, under polarized light, was also identified by HSM.
In summary, the present studied binary system of di-n-alkyl adipates proved to be a good candidate as PCM for TES applications at low temperatures.As far as the authors are aware, no low temperature TES research studies for the present binary system, has been reported in the literature.
Future work will be focused on the study of the thermophysical properties of this binary system, e.g., the thermal conductivity, which is imperative for the practical application as PCM.

Fig. 1
Fig. 1 DSC heating curves of pure diethyl, dibutyl, and of selected binary mixtures, with dibutyl molar fraction x dibutyl

Fig. 2
Fig. 2 Hot stage microscopy results for the x dibutyl = 0.1 mixture (a) and for the eutectic mixture (b) in temperature region between 193.15 K and 293.15 K, with 250× magnification

Fig. 9
Fig. 9 Binary solid-liquid phase diagram of the system of diethyl and dibutyl adipates; •• liquidus fitting curve according to the Eq. 2, Δ experimental liquidus line data points, ‫-•-‬ solidus line, ○ experimental solidus line data points, ▬ ▬ polymorphic line, □ experimental polymorphic transition temperature data points

Table 1
DSC data of pure diethyl adipate, pure dibutyl adipate and three of their binary mixtures polymorph 2 at ca. 238.15K on heating, with an enthalpy around Δ trans H = 2.50 J•g −1 .

Table 2
3 Absolute deviations of the experimental data obtained with thermal analysis, from the fitting Eq. 2, as a function of dibutyl adipate molar fraction Fitting parameters a and b of Eq. 2 and the correspondent relative root mean square deviation σ R for the experimental liquidus data points using DSC

Table 3
Comparison of the results obtained by fitting for L Mi and T Mi and the experimental results of Δ fus H and T fus for pure diethyl and dibutyl adipates According to Eqs. 3 and 4, it is possible to estimate the values for L Mi and T Mi from the fitting parameters a and b.In Table3the values of L Mi and T Mi obtained by the fitting of Eq. 2 are compared to the experimental values of Δ fus H and T fus for pure diethyl and dibutyl adipates.