# Structural Characterisation and Chemical Stability of Commercial Fibrous Carbons in Molten Lithium Salts

Blagoj Karakashov, Vanessa Fierro, Sandrine Mathieu, Philippe Gadonneix, Ghouti Medjahdi, Alain Celzard

► To cite this version:

Blagoj Karakashov, Vanessa Fierro, Sandrine Mathieu, Philippe Gadonneix, Ghouti Medjahdi, et al. Structural Characterisation and Chemical Stability of Commercial Fibrous Carbons in Molten Lithium Salts. Materials, 2019, 12 (24), pp.4232. 10.3390/ma12244232. hal-04962803

HAL Id: hal-04962803

https://hal.univ-lorraine.fr/hal-04962803v1

Submitted on 23 Feb 2025

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 Attribution 4.0 International License
Article

# Structural Characterisation and Chemical Stability of Commercial Fibrous Carbons in Molten Lithium Salts

Blagoj Karakashov <fnref n="1" />,*, Vanessa Fierro <fnref n="1" />,*, Sandrine Mathieu <fnref n="2" />, Philippe Gadonneix <fnref n="1" />, Ghouti Medjahdi <fnref n="2" /> and Alain Celzard <fnref n="1" />,*

Received: 19 November 2019; Accepted: 13 December 2019; Published: 17 December 2019

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

check for
updates

**Abstract:** The growing trend towards sustainable energy production, while intermittent, can meet all the criteria of energy demand through the use and development of high-performance thermal energy storage (TES). In this context, high-temperature hybrid TES systems, based upon the combination of fibrous carbon hosts and peritectic phase change materials (PCMs), are seen as promising solutions. One of the main conditions for the operational viability of hybrid TES is the chemical inertness between the components of the system. Thus, the chemical stability and compatibility of several commercial carbon felts (CFs) and molten lithium salts are discussed in the present study. Commercial CFs were characterised by elemental analysis, X-ray diffraction (XRD) and Raman spectroscopy before being tested in molten lithium salts: LiOH, LiBr, and the LiOH/LiBr peritectic mixture defined as our PCM of interest. The chemical stability was evaluated by gravimetry, gas adsorption and scanning electron microscopy (SEM). Among the studied CFs, the materials with the highest carbon purity and the most graphitic structure showed improved stability in contact with molten lithium salts, even under the most severe test conditions (750 °C). The application of the Arrhenius law allowed calculating the activation energy (in the range of 116 to 165 kJ mol⁻<fnref n="1" />), and estimating the potential stability of CFs at actual application temperatures. These results confirmed the applicability of CFs as porous hosts for stabilising peritectic PCMs based on molten lithium salts.

**Keywords:** carbon fibres; carbon felts; characterisation; chemical stability; high-temperature; lithium salts; peritectic compounds; thermal energy storage

## 1. Introduction

Population growth and the economic and industrial development of many countries around the world require and support the search for constant and stable sources of energy [1]. More importantly, global concerns about climate change or pollution make it crucial to develop new technologies and sustainable materials that improve yields, especially for storing renewable energy in its various forms. In recent decades, important scientific work has been done to meet the pressing demand for alternative technologies for the collection and storage of sustainable energies, such as solar energy [2–4]. Thermal Energy Storage (TES) is therefore considered as one of the main solutions to meet storage expectations [5]. The latter system should play an important role in the development of clean energy storage technologies, in the improvement of global energy management and in the incorporation of intermittent solar energy, or other sources of thermal energy, thus releasing a large excess of energy that is even more sustainable [6,7].

<footnote n="1">¹ Institut Jean Lamour, Université de Lorraine, CNRS, IJL, F-88000 Epinal, France; blagoj.karakashov@univ-lorraine.fr (B.K.); vanessa.fierro@univ-lorraine.fr (V.F.); philippe.gadonneix@univ-lorraine.fr (P.G.)</footnote>

<footnote n="2">² Institut Jean Lamour, Université de Lorraine, CNRS, IJL, F-54000 Nancy, France; sandrine.mathieu@univ-lorraine.fr (S.M.); ghouti.medjahdi@univ-lorraine.fr (G.M.)</footnote>

<footnote>* Correspondence: alain.celzard@univ-lorraine.fr; Tel.: +33-372-74-96-14</footnote>
To date, there are three types of TES systems, classified as sensible heat, latent heat and thermal chemical storage systems, or a combination of these [8]. In TES, latent heat, and more particularly phase-change materials (PCMs), play a key role in the implementation of this type of energy storage [9]. The storage via PCMs takes advantage of endothermic and exothermic isothermal phase-change processes (e.g., solid to liquid and vice-versa), exploiting increased energy densities involved at much higher temperatures than in sensible heat storage systems [10]. Thus, recent scientific advances make PCMs extremely attractive and effective for large-scale energy storage, with a small environmental footprint and a competitive price compared to other solutions [11].

Peritectic compounds (PCs) have been recently recognised as a different subcategory of PCMs, because they store energy, not only by melting/solidification, but also by a liquid/solid chemical reaction [11,12]. Indeed, PCs can be used as a one-pot storage solution, in which the chemical reaction is driven by temperature and carried out at atmospheric pressure, both in charge and discharge cycles, considered advantageous compared to gas-solid thermal chemical reaction storage systems [12]. The breakthrough of these promising PCs opens a new route for the development of novel ultra-compact systems, through the joint use of sensible, latent and thermo-chemical reaction, with high potential to provide economically viable TES [12].

However, previous studies have shown that the best performance of such PCMs and the improvement of their stability and efficiency can only be achieved by using host materials to confine them [13]. Thus, the use of carbonaceous materials, in the form of carbon fibres, carbon nanotubes, expanded graphite, graphene, alone or mixed together or with other materials, for the preparation of hybrid materials, is one of the most explored methods for improving PCMs [14]. Many properties of carbon-based materials, such as high chemical stability, low density and lower cost compared to different metals and ceramics, are indeed favourable for their use as supports of PCMs [15]. Furthermore, the use of fibrous carbons as porous hosts for PCMs, such as non-woven carbon felts (CFs), is considered a suitable solution in terms of cost, technological readiness, ability to be modified, and properties to achieve the requested TES performance [16].

The chemical inertness between the carbon host and the storage material is a critical factor in maintaining complete and reversible peritectic reactions, i.e., in the long-term, the effective operation of PC systems. Previously, the chemical stability of carbon fibres was considered a proven fact for their use as hosts of low-temperature organic PCMs [17–24]. However, future use of CFs in hybrid applications requires information on the compatibility and stability of composites under high-temperature working conditions and in the presence of anhydrous inorganic PCs. Some of the recently reported PCs [12] and other PCMs could potentially be corrosive in certain circumstances, depending in particular on the increase in the temperature of use, the presence of impurities (moisture, oxygen), and the occurrence of chemical reactions with carbon host materials [10,25,26]. Nowadays, the chemical compatibility of TES supports made of ceramics, inert asbestos or metals with molten inorganic PCMs, such as hydroxides, fluorides, carbonates, sulphates, chlorides, or their mixtures, is extensively studied, and the results have reported the technical challenges for developing long-lasting TES systems [27–29]. On the other hand, the chemical stability of carbon host candidates for high-temperature TES, via PCMs, and more specifically, via PCs, is not well examined, and therefore in-depth studies are required to test their chemical stability.

Consequently, the main purpose of this paper is the study of the materials' properties in terms of selection criteria for chemically-stable fibrous carbons, as well as the search for the probable reasons for the occurrence of unwanted reactions between the investigated carbon host and the PC of interest here. This paper therefore proposes the necessary examination of the elemental, structural and textural properties of CFs and their constituents, i.e., the carbon fibres of ex-polyacrylonitrile (PAN) or ex-Rayon origins, and studies of their high-temperature chemical stability in molten lithium salts, LiOH and LiBr, and their peritectic phase, Li₄(OH)₃Br. Research in the literature has revealed that there is no similar systematic study on the effect of carbon fibre properties on their chemical stability in contact with these salts in the molten state. Li₄(OH)₃Br is a unique inorganic anhydrous binary system of stoichiometric
composition, recently considered as a promising peritectic PCM intermediate for storing the heat produced by concentrated solar power plants or other thermal energy sources [30]. The last part of the study deals with the kinetics of the chemical reactions, using the Arrhenius law, thus considering the influence of temperature upon the chemical reaction rate. Fitting the results from CFs' chemical reactivity data in molten LiOH, the main purpose was to obtain information on the apparent activation energy of the dominant chemical reaction, thus making it possible to estimate the burn-off of the carbon in the temperature range of the application.

## 2. Materials and Methods

### 2.1. Materials

In the present work, 12 commercial carbon felts (CFs), received from five different suppliers, were sorted into groups for an easier evaluation of their characteristics (see Table 1). In the context of this work focusing on particular CFs, Karakashov et al. (2019) [31] recalled how they are produced, and how the origin of the materials and their manufacturing history can affect their final structural properties, and their physical properties in general.

<table><caption>**Table 1.** Sorting of the carbon felts (CFs) investigated here into groups, depending on their main characteristics.</caption><thead><tr><th>Final Heat Treatment</th><th>Precursor</th><th>Commercial Name</th><th>Used Sample Code</th></tr></thead><tbody><tr><td>Carbonised</td><td rowspan="7">Rayon</td><td>Carbon (Rayon) felt CeraMaterials</td><td>CR1</td></tr><tr><td>Graphitised</td><td>Graphite (Rayon) felt CeraMaterials</td><td>GR1</td></tr><tr><td>Carbonised</td><td>RSF1 Beijing Great Wall Co.</td><td>CR2</td></tr><tr><td>Graphitised</td><td>RSF2 Beijing Great Wall Co.</td><td>GR2</td></tr><tr><td>Carbonised</td><td>SIGRATHERM® KFA5</td><td>CR3</td></tr><tr><td>Graphitised</td><td>SIGRATHERM® GFA10</td><td>GR3</td></tr><tr><td>Graphitised</td><td>GF2 Schunk</td><td>GR4</td></tr><tr><td>Carbonised</td><td rowspan="5">PAN</td><td>PX 35 ZOLTEK™</td><td>CP5</td></tr><tr><td>Graphitised</td><td>GFE-1 CeraMaterials</td><td>GP6</td></tr><tr><td>Carbonised</td><td>BESF Beijing Great Wall Co.</td><td>CP7</td></tr><tr><td>Carbonised</td><td>Carbon (PAN) felt CeraMaterials</td><td>CP8</td></tr><tr><td>Graphitised</td><td>Graphite (PAN) felt CeraMaterials</td><td>GP8</td></tr></tbody></table>

A preliminary de-sizing [32,33] step was necessary to ensure the purity (in terms of carbon content) of the commercial CFs, thus avoiding the risk of additional chemical contamination and unwanted reactions with molten inorganic salts. Therefore, slow pyrolysis (heating rate of 3 °C min⁻¹ up to 900 °C followed by a dwell time of 2 h, under a flow of pure N₂) has always been performed after receiving the CFs and before doing any further materials' analysis. The set point of 900 °C was chosen, since the CFs were already industrially heat-treated at higher temperatures (1000–2400 °C), and thus was considered high enough for the carbonisation of any kind of polymer coating.

### 2.2. Chemical Reactivity Tests

Prior to the chemical reactivity tests, the inorganic powdered salts LiOH (anhydrous, ≥ 98%, Alfa Aesar, Kandel, Germany, melting temperature 462 °C) and LiBr (anhydrous, ≥ 99.8%, Alfa Aesar, city, country, melting temperature 550 °C) were always handled inside a glove bag (Glas-Col TM 108D X-27-27H, Terre Haute, IN, USA) under inert atmosphere because of the hygroscopic behaviour of these salts, and heated up to 120 °C to eliminate any unwanted water content.

To test the chemical stability of CFs in the molten peritectic compound (PC), an additional synthesis step was carried out by mixing LiOH and LiBr at their stoichiometric weight ratios (45 wt % LiOH/ 55 wt % LiBr), then synthesising the PC according to the protocol of Achchaq et al. (2019) [30].
This additional synthesis allows the formation of Li₄(OH)₃Br at its very precise stoichiometric composition, which is required for testing its effect on the investigated CFs.

Indeed, if a simple physical melting had been performed prior to the chemical reactivity studies, different experimental results could possibly be obtained because of the different melting temperatures of LiOH and LiBr.

With regard to the chemical reactivity test (shown schematically in Figure 1), the CFs with predefined dimensions were first placed in a nickel crucible (50 mL, VWR), then covered with calculated and fixed weight ratios of LiOH/CF (40:1), LiBr/CF (80:1) or Li₄(OH)₃Br/CF (40:1). Separate tests were performed with LiOH and LiBr to observe the effect of the different natures of these inorganic salts [34]. By doing so, a conclusion can be drawn for the final chemical stability of the CFs in molten Li₄(OH)₃Br. The other important experimental variable was the heat-treatment temperature, between 500 and 750 °C. The available information presented in the literature [35] is therefore used to establish the limits of the test and the occurrence of chemical reactions, significantly controlled by the temperature, thus clearly distinguishing the different CFs. In addition, to ensure complete impregnation, an additional smaller nickel crucible (10 mL, VWR) was placed on the lightweight CF to keep it still at the bottom and in contact with the molten LiOH, LiBr or Li₄(OH)₃Br. In this context, Achchaq et al. (2018) [36] used a simple infiltration protocol to confirm the complete infiltration process of the nonwoven CFs with the molten PC, and the corresponding process was also successfully performed. Then, the prepared test samples were transferred to a horizontal stainless steel tubular reactor (Carbolite Gero furnace, Neuhausen, Germany) constantly under an N₂ flow rate of 100 mL min⁻¹. The reactivity tests consisted of heating the sample from room temperature to different plateau temperatures (500–750 °C), with a heating rate of 5 °C min⁻¹ and a dwell time of 4h, before free cooling to ambient temperature. The solid salt block recovered with the CF immobilised therein was dissolved in a beaker filled with 2 L deionized water at 80 °C overnight. The liberated CF was transferred into a funnel with installed filter paper (Carl Roth GmbH+Co., Rotilabo®—round filter, type 114A, Ø 70 mm, Karlsruhe, Germany) and washed with 1 mol L⁻¹ HCl to neutral pH, and then dried at 105 °C overnight.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Chemical reactivity test protocol of carbon felts (CFs) in molten Li-based salts.</figcaption></figure>

The examined CFs were labelled according to the scheme presented in Figure 1. Whether carbonised (i.e., manufactured at a final production temperature of less than 1600 °C) or graphitised (i.e., manufactured at a final production temperature above 2000 °C), the first part of the CF code consisted of a capital letter G or C, respectively. The latter was followed by another capital letter, relating to the carbon fibre precursor (R for Rayon and P for polyacrylonitrile (PAN)), but also the number linking the sample with the commercial supplier code (see Table 1). After performing the chemical reactivity experiments, the presented code labels were followed by a low-case letter relating to the molten salts
tested as chemical reagents (a for LiOH, b for Li₄(OH)₃Br, or c for LiBr), ending with the temperature used. Thus, for instance, CR1a500 refers to the carbonised CF from the Rayon precursor, obtained from CeraMaterials, after a reactivity test in molten LiOH carried out at 500 °C.

## 2.3. Characterisation

### 2.3.1. Elemental Analysis

The ultimate composition of CFs was determined using a Vario EL Cube (Elementar, Langenselbold, Germany) elemental analyser. At least two samples of each CF were analysed to check the repeatability of the results. The weight fractions of C, H, N and S were determined from the gases obtained by burning a piece of sample in an excess of dioxygen at a temperature of about 1700 °C, using as the combustion catalyst a tin foil in which the samples were wrapped. Gaseous CO₂ and H₂O were directly obtained, whereas N₂ and SO₂ were produced after reduction of NOₓ and SOₓ on hot copper chips. These four gases were separated and quantified by a chromatographic column. As the amount of oxygen present in the sample could not be determined in the same run, additional samples were burnt, and the resultant CO₂ and H₂O were reduced into CO + H₂, and the quantification of CO led to the oxygen fraction in the materials. The same analysis was also used to observe changes in the elemental content of de-sized CFs with respect to as-received (pristine) CFs.

### 2.3.2. X-ray Diffraction (XRD)

High-resolution X-ray diffraction (XRD) was used to get the CFs' diffraction patterns, and was performed with a Bruker D8 Advance diffractometer (Billerica, MA, USA), at 50 kV and 45 mA, using a Mo anode (λ = 0.7093 Å) and a Lynxeye detector. The diffraction patterns were recorded in the 2θ range between 8 and 45° (with a scan step of 0.009°) at room temperature.

From the XRD patterns of the examined CFs' materials, parameters such as interlayer spacing (d₀₀₂) and crystallite size perpendicular to the carbon layers (L_c) were calculated from the position and the width of the (002) band, respectively. The Bragg equation [37]:

$$d_{002} = \frac{\lambda}{2 \sin \theta} \qquad (1)$$

where λ is the wavelength of the used radiation (λ = 0.7093 Å) and θ is the Bragg angle of the (002) reflection, was thus used to obtain the value of d₀₀₂, whereas L_c was calculated from the conventional Scherrer equation [38]:

$$L_c = \frac{K\lambda}{B\cos\theta} \qquad (2)$$

where K is a dimensionless shape constant with a present value of 0.9, and B is the full width at half maximum (FWHM) of the (002) reflection.

### 2.3.3. Raman Spectroscopy

The Raman spectra of selected carbonised and graphitised CFs were studied by performing Raman scattering with a Horiba Scientific XploRa Raman spectrometer (Kyoto, Japan), carried out on samples without additional preparation. The investigation was performed for both the outer surface and cross-section of carbon fibres (see Figure S1). The dispersion of the Raman-scattered light was done by a 1200 lines per mm holographic grating and detected by a CCD camera. A sample of silicon (521 ± 2 cm⁻¹) was used as a reference material. The Raman spectra were recorded with a laser at a wavelength of 532 nm, circularly polarised and filtered at 10% of its nominal power. Such filtration led to an incident power too low to overheat and damage the samples, and thus to influence the results. Once the measurement conditions were established, each sample spectrum was obtained by accumulating two spectra recorded in the range of 800 to 3500 cm⁻¹. The deconvolution of the bands was carried out with Horiba's LabSpec 6 software (Horiba France SAS, Longjumeau, France).
Typical Raman spectra of partially-graphitised or graphitised carbon materials comprise two main parts. The first order comprises $D_1$ band (around $1350 \text{ cm}^{-1}$), G band (around $1600 \text{ cm}^{-1}$) and sometimes $D_2$ band ($1620 \text{ cm}^{-1}$), and the second-order part mainly comprises the $S_1$ (or $2D_1$) band (around $2700 \text{ cm}^{-1}$) [39]. The G band corresponds to the C–C elongation mode in the aromatic rings, and therefore corresponds to the vibrations of the $sp^2$ carbons [40].

In contrast, the D bands observed in the first order part have been attributed to in-plane defects and to heteroatoms. The $S_1$ band is as an overtone of the $D_1$ band, and is related to the stacking order along the c-axis [41]. From the $D_1$ and G bands, their full width at half maximum (FWHM) values and intensity ratio ($I_{D1}/I_G$) are calculated, which gives information on the extent of structural defects and the graphitisation degree in a series of carbon materials.

### 2.3.4. Scanning Electron Microscopy

The carbon fibre microstructure was then evaluated by scanning electron microscopy (SEM) (Philips XL30 SFEG, Amsterdam, The Netherlands, equipped with field-enhanced thermionic emission gun operating at 4 kV) for CFs before and after chemical reaction with molten salts. The SEM micrographs of different samples were used to determine the alterations of the carbon fibres relative to the corresponding initial materials. Measurements of the mean fibre diameter (assuming circularity) were carried out for at least 10 fibres at five different locations in selected CFs, leading to a standard error of less than 10%, and thereby ensuring representative average values.

### 2.3.5. Textural Properties of CFs

The textural properties of the materials and their modifications were measured by gas adsorption analysis. All the samples were first degassed under vacuum for at least 48 h at 110 °C. Next, a manometric method was carried out by injecting a purified gas into a cell of known volume, already containing a known amount of adsorbent material. The nanoporous structure of the pristine and chemically-reacted CFs were derived from nitrogen ($N_2$) and carbon dioxide ($CO_2$) adsorption at –196 °C and 0 °C, respectively, performed with automatic adsorption analysers ASAP 2020 and ASAP 2420 (Micromeritics, Norcross, GA, USA), respectively.

Using Microactive® and SAIEUS® software (version 3.0, Norcross, GA, USA), the following parameters were calculated by application of the 2D non-local density functional theory for the heterogeneous surface of pore walls (2D Non-Localized Density Functional Theory (NLDFT)-HS) to $N_2$ and $CO_2$ isotherms: specific surface area, $S_{NLDFT}$ (m² g⁻¹); total pore volume, $V_{tot}$, 2D-NLDFT HS (cm³ g⁻¹); micropore volume, $V_μ$, 2D NLDFT-HS (cm³ g⁻¹); and mesopore volume, $V_{meso}$ (cm³ g⁻¹) = $V_{tot}$, 2DNLDFT-HS – $V_μ$, 2D NLDFT-HS.

The 2D NLDFT-HS model, rather than the traditional BET model, was used here because it improves the modelling by taking simultaneously into account the adsorption results of nitrogen ($N_2$) and carbon dioxide ($CO_2$), and considering carbon pore walls with heterogeneous surfaces. Additionally, this modern method combines all pore widths to give a more accurate pore size distribution, thus better representing carbonaceous materials with a poorly developed porous structure and a low surface area [42].

## 3. Results and Discussion

### 3.1. CFs Elemental Composition

In addition to differences in manufacturing processes and precursors [31], CFs also differ in terms of carbon and heteroelement contents. The performance of the analysis showed that the elemental composition of commercial CFs is highly repeatable, with maximal standard deviations of 0.45 wt % for C and O, or 0.01 wt % for the other heteroatoms. As no sulphur was found, this element was omitted from the discussion of the results, given in Table 2 in units of wt %. As expected, C was by far the most abundant element in as-received commercial materials, with minimum values of 76 wt % and
95 wt % for ex-Rayon CFs and ex-PAN CFs, respectively. The carbon amount was further increased when the CFs derived from a given precursor were manufactured at a higher temperature. Oxygen was the second most abundant element, with a maximum value of 18.4 wt % for the ex-Rayon CFs and a value as low as 0.6 wt % for the ex-PAN CFs. A meaningful reason for the presence of O and N is that the regenerated cellulose (Rayon) fibrils are rich in O, and that N can also be originally present as a production impurity, whereas PAN contains large numbers of C≡N functionalities [43–45]. Compared to the C content, the amounts of H were found to be low and variable from one material to another, but decreasing again with the increase in the final heat treatment temperature. In general, pyrolysis in inert atmosphere controls the increase of the carbon content through several complex elimination processes of CH₄, H₂, N₂, HCN, H₂O, CO, CO₂, NH₃, tars and other volatiles [46–48].

<table><caption>**Table 2.** Elemental analysis results (in wt %) of as-received and de-sized CFs.</caption><thead><tr><th rowspan="2">Sample Code</th><th rowspan="2"></th><th colspan="4">As-Received CFs</th><th colspan="4">De-Sized CFs</th></tr><tr><th>C (%)</th><th>N (%)</th><th>O (%)</th><th>H (%)</th><th>C (%)</th><th>N (%)</th><th>O (%)</th><th>H (%)</th></tr></thead><tbody><tr><td rowspan="6">Ex-Rayon CFs</td><td>CR1</td><td>75.95</td><td>3.46</td><td>18.37</td><td>2.14</td><td>77.66</td><td>3.18</td><td>16.63</td><td>2.45</td></tr><tr><td>CR2</td><td>77.55</td><td>3.17</td><td>17.15</td><td>2.07</td><td>78.45</td><td>2.77</td><td>16.40</td><td>2.33</td></tr><tr><td>CR3</td><td>90.94</td><td>0.50</td><td>7.64</td><td>0.92</td><td>89.16</td><td>0.12</td><td>9.27</td><td>1.46</td></tr><tr><td>GR1</td><td>91.27</td><td>0.35</td><td>7.46</td><td>0.92</td><td>97.83</td><td>0.13</td><td>1.92</td><td>0.12</td></tr><tr><td>GR2</td><td>98.70</td><td>0.43</td><td>0.78</td><td>0.09</td><td>99.01</td><td>0.10</td><td>0.81</td><td>0.08</td></tr><tr><td>GR4</td><td>99.53</td><td>0.01</td><td>0.39</td><td>0.07</td><td>99.66</td><td>0.03</td><td>0.25</td><td>0.05</td></tr><tr><td rowspan="4">Ex-PAN CFs</td><td>GR3</td><td>99.66</td><td>0.04</td><td>0.24</td><td>0.07</td><td>99.72</td><td>0.01</td><td>0.22</td><td>0.05</td></tr><tr><td>CP5</td><td>96.82</td><td>2.82</td><td>0.31</td><td>0.05</td><td>97.01</td><td>2.52</td><td>0.41</td><td>0.06</td></tr><tr><td>CP8</td><td>95.07</td><td>4.09</td><td>0.65</td><td>0.19</td><td>95.19</td><td>3.91</td><td>0.60</td><td>0.23</td></tr><tr><td>CP7</td><td>97.89</td><td>1.59</td><td>0.46</td><td>0.06</td><td>97.85</td><td>1.71</td><td>0.35</td><td>0.09</td></tr><tr><td rowspan="2"></td><td>GP6</td><td>99.46</td><td>0.03</td><td>0.44</td><td>0.06</td><td>99.58</td><td>0.02</td><td>0.32</td><td>0.08</td></tr><tr><td>GP8</td><td>99.44</td><td>0.04</td><td>0.46</td><td>0.06</td><td>99.68</td><td>0.04</td><td>0.22</td><td>0.06</td></tr></tbody></table>

A complementary elemental analysis allowed comparing as-received CFs with de-sized CFs. In the case of all ex-PAN CFs, no significant changes were observed in terms of elemental composition. However, in some cases of carbonised, ex-Rayon CFs (CR1 and CR2), weight losses (averaged on a dry basis) of 7% and 8%, respectively, were observed after the de-sizing treatment. This could produce a slight decrease of O content and an increase of C content. Another exception is the graphitised, ex-Rayon GR1, for which an increase in C content and a decrease in O content were seen, without clear net weight loss. These observations of changes in the elemental composition and weight losses evidence the sizing [49] or any other post-production modification of some CFs, and confirm the need to proceed to an additional de-sizing step before performing all characterisation analyses and chemical reactivity tests.

Van Krevelen diagrams of ex-Rayon and ex-PAN CFs confirm that both O/C and H/C ratios decrease, and that the degree of carbonisation increases with the final treatment temperature (see Figure 2a,b, respectively). In general, most ex-PAN CFs present lower O/C and H/C ratios than ex-Rayon CFs, due to the different elemental composition of the precursors. It can be noted that the y-axis is used to represent the H/C ratios, whereas the x-axis describes the variation of the O/C ratio, with significant differences in scales. This indicates that the elimination of O is favoured, but that a higher number of H atoms is still present in the final structure of the materials. This is clearly visible for the graphitised CFs. Thus, the thermal annealing of CFs leads to a sharp drop of ratios for most of the corresponding materials (such as CP8 → GP8; CR1 → GR1; CR2 → GR2; and CR3 → GR3), also suggesting the improvement of the structural and chemical stability of the carbon by decreasing the heteroatom content. Remarkably, it is evident that different precursors or fabrication processes produce CFs having different elemental compositions, while all were produced in an inert atmosphere, and with similar final heat-treatment temperatures.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Van Krevelen diagrams of: (a) ex-Rayon and (b) ex-polyacrylonitrile (PAN) CFs. The dashed lines connect CFs before (empty symbols) and after (full symbols) graphitisation of the CFs.</figcaption></figure>

## 3.2. CFs Nanostructure

### 3.2.1. XRD Analysis

The XRD patterns of the 12 commercial materials are given in Figure 3a,c for ex-Rayon CFs and ex-PAN CFs, respectively, and these are presented as a function of the gradual increase in their structural order. Figure 3b for ex-Rayon and Figure 3d for ex-PAN CFs relate the characteristics of the (002) reflection: position, width and intensity, to the extents of interlayer spacing (d<sub>002</sub>) and of graphitic layers' stacking (L<sub>c</sub>) using Equations (1) and (2), respectively.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** X-ray diffraction (XRD) patterns of: (a) ex-Rayon CFs, and (c) ex-PAN CFs, presented as a function of the gradual increase in their structural order. The arrows connect corresponding CFs, and the patterns are shifted with respect to each other by a constant number of counts per second (CPS, arbitrary units). Corresponding results of interlayer spacing (d<sub>002</sub>) and nanocrystallite thickness (L<sub>c</sub>), plotted against increasing structural order for: (b) ex-Rayon CFs, and (d) ex-PAN CFs.</figcaption></figure>
The (002) reflections of the ex-Rayon, carbonised CFs appear very broad and not completely
defined, especially for the materials CR1, CR2 and CR3, because of the large number of highly
disordered carbon regions. Their low structural order results from semi-crystalline conformations of
the precursor, which have been found to be connected with disordered transition regions containing
heteroatoms, thus limiting rearrangement during pyrolysis processes [50,51]. As for the corresponding
GR1, GR2 and GR3, an absence or a partial improvement of the structural order was observed after
the graphitisation step (above 2000 °C). The higher heat-treatment temperature clearly defined these
materials as non-graphitising (GR1) or partially-graphitising (GR2, GR3). The splitting of the (002)
reflection suggests a structural heterogeneity of the partially-graphitising material (clearly seen for
GR2), due to the presence of graphitic nanocrystallites of different sizes and different interlayer
spacing [52]. The graphitised GF4 gave narrower reflections slightly shifted to a higher Bragg angle.
The performance of pyrolysis at low heating rates and the use of impregnants or organosilicon
compounds (stabilisers) [53,54] can explain the structural improvements with the highest evaluated
structural order. In general, the XRD patterns of the graphitised, ex-Rayon CFs also suggest the
absence of a hot stretching step at the late carbonisation stage in the manufacturing process of
the carbon fibre [53], crucial for the development or improvement of the carbon structural order.
The diffraction patterns obtained from ex-PAN CFs prepared at temperatures below 1200 °C exhibited
a better ordered carbon structure than that of carbonised, ex-Rayon CFs. Graphitised, ex-PAN CFs are
materials with improved turbostratic structure, resulting in sharper and narrower (002) reflections [45]
compared to carbonised CFs. However, no tri-periodic stacking formation of graphitic structure is
observed in ex-PAN or ex-Rayon CFs, according to the low definition or the absence of (hkl) reflections.
Moreover, the discussion of the structural properties of CFs is omitted with respect to two-dimensional
XRD reflections (hk0) because of their poor appearance and the possibility of important estimation
errors [55,56].

The samples having a very poor structural arrangement produce broad (002) diffraction bands
of asymmetrical appearance (CR3, GR1 and, to a certain extent, GR2), the deconvolution of which
did not result in accurate fits. When they could be clearly deconvoluted, the (002) reflections were
used to calculate the values of structural parameters d<sub>002</sub> and L<sub>c</sub>. The values of d<sub>002</sub> and L<sub>c</sub> decrease
and increase, respectively, for materials having a higher structural organisation and alignments of
graphite nanocrystallite layers. In general, d<sub>002</sub> values are observed in the range from 3.5 to 3.68 Å for
ex-Rayon CFs, and from 3.5 to 3.6 Å for ex-PAN CFs, i.e., these are much higher than those of graphite,
3.35 Å [57]. Therefore, the low structural order and the absence of tri-periodic graphitic stacking have
been confirmed for these fibrous carbons, assumed to be of general use because they are produced
without advanced methods of structural improvement.

Comparing elemental compositions and structural properties was then performed to obtain
additional cross-property relationships. For this purpose, Figure 4 shows the O/C and H/C atomic
ratios as a function of the interlayer spacing (d<sub>002</sub>) for ex-Rayon and ex-PAN CFs.

There is a reasonable correlation between the decrease of H and O contents with respect to that
of C, and the increase of structural order. In general, d<sub>002</sub> decreases when reducing the amount of
heteroatoms in the carbon lattice, which allows additional arrangement within the graphitic layers.
These observations confirm that the structural information derived from the (002) reflection also reflects
the degree of carbon purity in the CFs, indicated previously for other carbons [58,59]. The ex-Rayon
CFs have slightly scattered results, but still follow the same trend as ex-PAN CFs. The corresponding
pairs of CFs from the two precursors, i.e., CR2 vs. GR2 and CP8 vs. GR8, show a very good correlation
between the elemental composition and the properties of the graphitic nanocrystallites. Moreover,
for both ex-PAN CFs and ex-Rayon CFs, the impact of the structural properties of the precursor and
the preheating treatment steps is clearly observed for carbonised CP5 and for graphitised GR2 and
GR3, respectively. Although the carbonised CP5 sample shows comparable or higher carbon purity
and structural properties than many graphitised CFs, the latter samples have a high purity, but only a
partially or completely non-graphitic structure. As a result, a moderate structural improvement is
often observed for biosourced carbon fibres, even after having been subjected to high-temperature treatments [52,53,60].

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Elemental composition as a function of interlayer spacing (d<sub>002</sub>) of: (a) ex-PAN CFs, and (b) ex-Rayon CFs. Empty symbols represent carbonised CFs and full symbols represent graphitised CFs. The difference in colour corresponds to O/C (blue) and H/C (black) atomic ratios.</figcaption></figure>

### 3.2.2. Raman Spectroscopy Analysis

After carrying out the overall structural XRD analysis, Raman spectroscopy was performed for local structural studies of selected CFs. Random zones on different carbon fibres were examined from the outer surface or cross-sections of the fibres (presented in Figure S1). Figure 5 shows the Raman spectra of the studied CFs, shifted with respect to each other for easier visualisation of the similarities or differences. Moreover, the Raman spectra obtained from the same CFs at the cross-section (dashed lines) and the outer fibre surface (solid lines) are completely superimposed, suggesting structural homogeneity. A first-order part is seen (in the range from 1000 to 1800 cm⁻¹) presenting an intense and broad D band (around 1350 cm⁻¹) and a slightly narrower G band (around 1590 cm⁻¹). A second-order part (in the range from 2200 to 3400 cm⁻¹) is also observed, with various degrees of structuration, depending on the investigated CFs. Therefore, clear differences in spectra are observed when comparing the corresponding CFs, with the only difference being in the final manufacturing treatment temperature, as seen for CR2 and GR2 and for CP8 and GP8 (linked by arrows).

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Raman spectra of: (a) ex-Rayon CFs and (b) ex-PAN CFs, where the arrows connect the corresponding CFs. The dashed or solid lines represent spectra of outer fibre surface or cross-section of carbon fibres, respectively. For easier viewing, an identical intensity shift has been applied.</figcaption></figure>
From the deconvolution of the Raman spectra (see Figures S2 and S3), the features of carbonised
CFs are considered very typical of highly disordered carbon fibres [60–63]. Fabricated with additional
heat treatment, the graphitised CFs have narrower first-order bands, with a difference in band intensity
and a much deeper valley between D₁ and G bands.

The bands D₁, D₃ and D₄ correspond to in-plane defects, out-of-plane defects and impurities
in the sp² graphitic network, respectively [40], whereas the G band solely represents the stretching
vibration in the aromatic layers of graphite-like carbon [59]. The constant presence of the G band in the
examined CFs explains the clear existence of sp²-hybridised carbon, and the intensity and width of this
band are both influenced by the three dimensional (3D) evolution of the graphitic nanocrystallites. It is
also found that the D₂ band, observed only for graphitised CFs, is related to surface graphene sheet
vibrations, and can therefore be used as a possible indicator of the surface-to-volume ratio of crystallites
in nanocrystalline, graphitic carbons [39]. Changes in these bands result from the improvement of the
structural order of the materials, accompanied by the decrease of D₄ and D₃ and the appearance of
the D₂ band. The presence of D₁ in all CFs samples indicates sp³-hybridised carbon atoms, as such
or more probably related to the presence of oxygen, hydrogen and/or nitrogen in the CFs, in good
correlation with the results of elemental analysis. In addition to the structural disorder, the D₁ band
also comes from the presence of small graphitic nanocrystallites with a low number of graphitic layers,
confirmed by the results of XRD.

The same kind of conclusion is obtained from the estimated full width at half maximum (FWHM) values of the deconvoluted D₁ and G bands, ranging from 169 to 48 cm⁻¹ and from 79 to 45 cm⁻¹, respectively (see Figure 6). The corresponding pairs of CFs (from the two precursors) show a very good correlation of the FWHM of the D₁ and G bands with the increase of structural order, comparing CR2 and CP8 with GR2 and GR8, respectively. However, even the lowest FWHM (GR4, at 45 cm⁻¹), is much higher than that of highly-oriented graphite, about 15–16 cm⁻¹ [40]. Therefore, the analysis indicates a low degree of local graphitic order in the examined CFs, which confirms the validity of the d₀₀₂ and L_c values from the XRD analysis.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Full width at half maximum (FWHM) of deconvoluted D₁ and G bands of ex-Rayon CFs (left) and ex-PAN (right) CFs, where the arrows connect the corresponding CFs. Empty symbols represent carbonised CFs and full symbols represent graphitised CFs, the colour differences corresponding to D₁ (black) and G (red).</figcaption></figure>

When the band deconvolution is performed for the second order of the Raman spectra, a clear
appearance of the intense S<sub>1</sub> band (attributed to the first overtone of D<sub>1</sub> at around 2700 cm<sup>-1</sup>) has been
found only for graphitised CFs. The observed intensity increase of the S<sub>1</sub> band (as in GR4) should
correspond to the lateral growth of the graphitic nanocrystallites, and to the improvement of the
stacking order of the graphitic layers [41,60]. However, it is known that the S<sub>1</sub> band splits into two
bands when the structure acquires a tri-periodic stacking order of the carbon layers [39,64], which is obviously not the case for any of the studied materials, in agreement with the XRD analysis.

The intensity of the D₁ band (I<sub>D1</sub>) is higher than that of the G band (I<sub>G</sub>) for the examined CFs, including the ones (GR2, GR4, GP6 and GP8) manufactured with a higher final heat-treatment temperature, above 2000 °C. For the corresponding CFs pairs (see Figure S4), this I<sub>D1</sub>/I<sub>G</sub> ratio increases gradually, which is an additional indication that these CFs are in the carbonisation regime, whereas in the graphitisation regime, one should expect the contrary [38]. Again, it is likely that only lateral graphite nanocrystallite growth occurs with a supposed lack of the tri-periodic stacking of carbon layers.

### 3.3. High-Temperature CFs Chemical Stability in Molten Li-Based Salts

After having presented the results of elemental composition and structural properties, the chemical stability of these commercial CFs is considered below. Many experimental variables, such as precursor, manufacturing technique or temperature, have been observed to influence the structural and textural properties of the CFs, and therefore should affect their chemical stability in molten Li-based salts. Knowledge of the chemical stability of the CFs is needed to identify the best materials to serve as host structures for PCMs for the intended TES application. It should be emphasised that, as discussed in Section 2.2, the chemical reactivity tests were repeated at least twice for each of the many CFs examined, showing a very low standard error of less than 2% on the carbon mass loss (defined as burn-off, B-O (%) = 100 – percent of carbon disappeared by initial weight of CF). Finding a high reproducibility of the reactivity tests makes it possible to compare the CF samples tested in molten LiOH, LiBr and Li₄(OH)₃Br, and to correlate the textural and morphological modifications studied.

#### 3.3.1. Chemical Stability of CFs in LiOH with Regard to Their Structural Properties

The chemical stability tests revealed that losses of carbon mass (i.e., B-O) always occurred, and increased with temperature (values indicated in Tables S1 and S2). This is obvious because of the considerable increase in the reaction rate and the existence of surface functional groups, considered as catalysts for the following reaction [35]:

$$2\text{C} + 6\text{LiOH} \rightarrow 2\text{Li} + 3\text{H}_2 + 2\text{Li}_2\text{CO}_3 \quad \text{(Reaction 1)}$$

Clearly, for both ex-PAN and ex-Rayon materials, the reactivity of LiOH is reduced as carbon purity increases, and the CFs become very stable even at the highest temperature tested (see Figure 7a,b). This confirms that the rate of reaction (1) and its activation energy can be modified when testing the reaction of LiOH with carbonaceous materials having different graphitisation stages. In general, by comparing only carbonised materials obtained from different precursors, the ex-Rayon CFs show a higher B-O than the ex-PAN CFs, in good agreement with their higher number of heteroatoms. Moreover, the differences in reactivity are clearly observed between the corresponding CFs (CR1a → GR1a, CR2a → GR2a, CR3a → GR3a or CP8a → GP8a). In this context, the presence of O, N and H in the CFs is associated with the existence of surface functional groups (–C–H; –C–OH; –C–O–C; –C=O; C=N; –COOH; etc.), the number of which decreases with subsequent heat treatments at a higher temperature [33,43,45,65–69]. Thus, the B-O results correlate well with the Raman analysis, where the carbonised CFs spectra showed broad and intense D₁, D₄ and D₃ bands, indicating a significant number of heteroatoms and aliphatic side chains.

In addition to the discussed correlation with their elemental composition, the graphitic improvement of the graphitic order of the CFs is also compared to the results of the chemical stability tests. As seen in Figure 7c,d, the values of B-O are generally correlated with the interlayer spacing (d<sub>002</sub>). The plotted results clearly show the downward trend in the CFs’ chemical reactivity in molten LiOH as the structural order increases, but this correlation should be considered with caution. Indeed, some CFs had a similar interlayer spacing (d<sub>002</sub>), such as in GR3a ↔ GR2a and GP8a ↔ GP6a, but with different B-O values. Coming back to Figure 4, a clear indication of this behaviour can be inferred from the
different amounts of O and H in the mentioned CFs. Thus, even if a similar structural development is observed, especially for GP8a ↔ GP6a (from both XRD and Raman analyses), the higher number of heteroatoms in the latter clearly distinguished this material as being more reactive with LiOH. A slight difference from the downward trend is also observed for CR2a and CR1a, again in good correlation with the elemental composition. Thus, the current correlations show that LiOH is more reactive (higher B-O at a lower test temperature) in the presence of low-purity CFs. In most cases, this was less dependent on the graphitic order achieved during the last manufacturing process (carbonisation or graphitisation), taking into account CFs from different suppliers. However, the graphitised CFs with improved structural order and higher crystallite size clearly appear as better candidates in the TES application, which should further improve the thermal management properties of the PCM-CF composite materials [70,71].

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Burn-off (B-O) (%) as a function of reaction temperature for: (a) ex-Rayon CFs and (b) ex-PAN CFs. Empty and full symbols represent carbonised CFs and graphitised CFs, respectively; the dashed or solid lines are only guides for the eye. Three dimensional (3D) bar representation of changes in burn-off (%) versus test temperature and interlayer spacing (d<sub>002</sub>) for: (c) ex-Rayon CFs and (d) ex-PAN CFs, tested for chemical reactivity in molten LiOH.</figcaption></figure>

### 3.3.2. Chemical Stability of CFs in LiOH with Regard to Their Textural Properties

A higher structural order also prevents the development of porosity in carbon fibres, as presented by the Franklin carbon models [72] in Figure 8. The imperfections of the graphitic network are reduced by the development of the graphitic nanocrystallites, thus inhibiting the microporosity.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Franklins' structural models of different carbon nanostructures and the presence of nanopores.</figcaption></figure>

Therefore, since the chemical stability of CFs has been discussed with respect to their elemental and structural properties, the same correlations have been sought with their nanotexture properties. The textural analysis carried out was based on N₂ and CO₂ adsorption isotherms on the initial (de-sized) CFs and those reacted in molten Li-based salts (the detailed results of textural analysis are available in Section S2 of the Supplementary Material). A large variety of nanotextures has indeed been observed for CFs with different structural properties, as shown by the adsorption isotherms of Figures S5–S7. Then, the application of the 2D NLDFT-HS model was chosen as quite advantageous in terms of the improved accuracy of the calculated textural properties. Moreover, the use of an up-to-date calculation model is very important when adsorption equilibrium is difficult to reach during the nanotexture analysis of carbons with heterogeneous and very narrow porosity [73].

The 2D NLDFT-HS surface areas obtained from the CFs reacted in molten LiOH were found to vary between 9.6 and 1034 m² g⁻¹. Although LiOH showed greater reactivity (given the high B-O values) with carbonised CFs with respect to graphitised ones, the same upward trend was not observed from the results of the 2D NLDFT-HS surface area. As seen in Figure 9a,b,d, CR1a, GR1a and GR2a show a monotonous increase of surface area, at a test temperature of 500 or 600 °C. In addition, CR2a (Figure 9c) is an exception to the monotonic trend with an initial increase (from 907 to 986 m² g⁻¹ at 500 °C) followed by a decrease in surface area at 600 °C (814 m² g⁻¹). By further increasing the test temperature, the carbonised CR1a and CR2a were completely consumed, and textural analysis was not possible. Whereas, their corresponding graphitised GR1a and GR2a again showed a slight increase of surface area when tested at higher temperatures, 650 and 700 °C, respectively.

These results of the low development of porosity, or even loss of porosity, together with the high values of B-O, are observed due to the high LiOH/CF molar ratio of 20/1, leading to an intense chemical reaction with the carbon, whose reaction rate increases with temperature. Moreover, a high chemical reactivity in molten LiOH has been reported elsewhere for different carbon materials, even when using a much lower LiOH/C molar ratio of 1.5/1 [74]. Another important feature is the presence and role of amorphous carbon that disorderly crosslinks graphitic nanocrystallites in the carbon fibre structure, leading to the presence of accessible nanopores and heteroatoms. When comparing the carbonised CFs to the graphitised CFs, the fraction of disordered phase is larger, and the very accessible atoms are more likely to perform chemical reactions than the carbon atoms of the graphitic nanocrystallites. The opposite was evident when analysing GR3a and GR4a CFs, as these materials have a significantly higher structural organisation and carbon purity. Even at higher temperatures, 700 and 750 °C, LiOH reacts to a much lesser extent with these CFs, resulting in negligibly low surface areas (see Figure 10a,b).

These results confirm the information of high structural order, deduced from the Raman spectroscopy (as for GR4). In addition, the results are in good correlation with the lower oxygen content, whose presence has been shown to facilitate the fixation of alkaline ions, thus catalysing carbon gasification reactions [75].

The same trends of surface area were observed when analysing the textural properties of the ex-PAN CFs. Likewise, the carbonised CFs (represented by CP8a) showed a higher surface area than the corresponding graphitised GP8a, after being heat-treated in molten LiOH at an identical test temperature. However, textural analysis was performed on only a few reacted ex-PAN CFs, due to the poorly developed porosity and high experimental error risk. The obtained results are presented in Table S2.
<figure><img src="image_11.png" /><figcaption>**Figure 9.** 2D Non-Localized Density Functional Theory (NLDFT)-HS textural properties: specific surface area (scatter symbols, where the dotted or dashed lines are only guides for the eye), micro and mesopore volumes (stacked columns), and pore size distributions (3D plots) of corresponding CFs: (a) CR1 and (b) GR1, (c) CR2 and (d) GR2, as a function of reaction temperature in molten LiOH or $Li_4(OH)_3Bi$.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 10.** 2D NLDFT-HS textural properties: specific surface area (scatter symbols, where the dotted or dashed lines are only guides for the eye), micro- and mesopore volumes (stacked columns), and pore size distributions (3D plots) of graphitised CFs: (a) GR3 and (b) GR4.</figcaption></figure>

In addition, the total pore volume and pore size distributions (PSDs) obtained with the 2D NLDFT-HS model further confirm the observed trends for the specific surface areas. Figures 9 and 10 also show that very different porosities are obtained, in terms of micro-/mesoporous volume, depending on the initial CFs properties and their reactivity in molten LiOH at the temperatures tested. In general, for all CFs examined, the increase in micropore volume occurred only by widening the initial micropores, but also because of the appearance of medium-sized micropores (0.7–2.0 nm) and mesopores, in good correlation with the PSDs. At a lower temperature of 500 °C, the PSDs were slightly modified by widening the medium-sized micropores of the carbonised CFs, due to the slow chemical reaction (see the PSDs in Figure 9a,c). As for the corresponding graphitised CFs (Figure 9b,d), the initial PSD was obviously different due to structural differences, leading to narrower pores. Thus, GR1 had both ultramicropores (< 0.7 nm) and medium-sized micropores, the latter extending their range to mesopores, whereas GR2 was purely microporous, in good agreement with the structural information. When testing the latter two CFs at 500 °C, it was observed that the chemical reaction in LiOH resulted in a widening of the medium-sized micropores in GR1a500, whereas medium-sized micropores and mesopores formed in GR2a500. In both CFs, this clearly influenced the increase in mesopore volume (from 9.0% to 20.5% and from 0.1% to 20.4%, respectively). Increasing the test temperatures led to an increase in the total pore volume in most of the CFs studied. For poorly ordered carbonised CFs (CR1 and CR2), this was mainly due to the significant opening of the microporosity initially observed and the increase in mesoporosity (e.g. from 0.30 to 0.86 cm³ g⁻¹ for CR1a600, and from 0.23 to 0.59 cm³ g⁻¹ for CR2a600 at 600 °C). The PSDs of the aforementioned CFs had merged regions of micro- and mesopores, which even caused a clear decrease in the specific surface area of CR2a600. With regard to the corresponding graphitised CFs (GR1 and GR2, respectively) of improved structural order, only the medium-sized microporosity was involved in the development of a higher micro-/mesopore volume without clear
evolution of the initial ultramicropores (see Figure 9b,d). Finally, the examined textural results of GR3 and GR4 with higher carbon purity and improved structural order confirm the important effect of the initial properties of the materials on the significant reduction of their reactivity in molten LiOH (see Figure 10a,b).

It is difficult to make an in-depth comparison between the results presented here and those obtained for different carbon materials, because of the different variables used, such as carbon precursor and experimental conditions. In general, the reviewed bibliography is consistent with our results on the chemical reactivity of CFs in molten LiOH and the observed loss of carbon weight, but with poor development of surface area and microporosity [74,76–83]. The main reason is always indicated as being the strong reaction of LiOH under the experimental conditions used and the high ratio LiOH/carbon material. The small size of the alkaline atom should also have promoted high reactivity, allowing easier access to the initial material porosity and further pore widening [84]. This resulted in a higher carbon mass loss in carbonised CFs compared to graphitised CFs. No metallic Li intercalation was observed in any of the CFs examined, as evidenced by XRD and supported by the textural analysis, still showing poor porosity development. This can be explained by the fact that Li cannot be converted to metal vapour (in the temperature range studied, up to 750 °C), then easily diffused and intercalated between the carbon layers [75,84–86]. This is in a good agreement with the similar chemical stability of CFs with a different structural order, but with high carbon purity (see Figure 4).

### 3.3.3. Chemical Stability of CFs in LiOH with Regard to Their Morphology

SEM micrographs of the carbon fibre morphology of de-sized materials before and after reaction in molten LiOH are shown in Figures 11 and 12 for ex-Rayon CFs and ex-PAN CFs, respectively. As detailed in a former study from our group [31], the morphology of the outer surface and the cross-section of de-sized carbon fibres prior to reaction varies with the precursor and the manufacturing technology, since these commercial materials come from different suppliers.

Moreover, the only observations common to these carbon fibres are grooves and striations parallel to the axis of ex-Rayon fibres and of ex-PAN fibres, respectively, because of the manufacturing process of the fibres' precursor. The microscopic impurities visible on the surface of most carbon fibres are probably by-products of de-sizing. Ex-Rayon CFs had pores between the constituent fibrils (as seen for GR1), clearly observed from the inhomogeneous cross-section of the carbon fibres. Whereas, ex-PAN CFs showed an obvious absence of these pores both on the outer surface of the fibres or in the cross-section. Therefore, further investigations were conducted to determine how the carbon fibre structure was modified during surface contact with LiOH, either uniformly or in a defined pattern, due to the various characteristics of the different CFs.

Most carbon fibres constituting the CFs retain their original global structure after the reactivity tests. However, some of them have been destroyed due to a severe chemical reaction with LiOH (as GR1a700), in good agreement with the results of B-O at high temperatures. From carbonised ex-Rayon CFs, we observed that the prolonged chemical reaction with LiOH deepened the initial grooves and created new cracks, clearly triggering the longitudinal division of the carbon fibres (as shown by the carbonised CR1a600).

The same is not clearly observed in the other carbonised CR2a600 and CR3a600 when they reacted at the same temperature, but the destruction of the fibres occurred after the increase in the test temperature (CR3a650) (see Figure S8). This confirms the use of different materials and methods of manufacturing ex-Rayon carbon fibres and their influence on the structure of the final materials, even with similar final heat-treatment temperatures. The influence of the materials' manufacturing history was also seen when comparing the graphitised ex-Rayon CFs. The sample GR1a700 was almost completely oxidised in the reaction with LiOH, and only indeterminate particles were observed as residues. On the contrary, in the case of samples GR3a750 and GR4a750, even after activation at the highest temperature tested, only pores appear on the initial surface, caused by the etching of disordered carbon zones. Moreover, the core of these latter samples remains unchanged in the reacted carbon fibres, again confirming their
higher purity and structural order. Thus, the larger pits observed in GR3a750 relative to GR4a750 can be attributed to the higher fraction of disordered domains, originating from heterogeneities in cellulosic fibres [48]. In addition, the chemical stability of GR1, GR3 and GR4 CFs in LiOH must also be the result of their manufacturing differences from non-, partially- or well-stabilised Rayon fibres, respectively [53]. Therefore, SEM microscopy data confirm structural properties' results, and ex-Rayon carbon fibres with less disordered domains should always be preferred when high chemical stability of CFs in molten LiOH is required.

<figure><img src="image_13.png" /><figcaption>**Figure 11.** Scanning electron microscopy (SEM) micrograph of de-sized ex-Rayon CFs before (**left**), and after reaction with either LiOH (**middle**) or Li₄(OH)₃Br (**right**) at the highest temperature tested.</figcaption></figure>

For ex-PAN CFs, no swelling or longitudinal peeling of carbon fibres was observed, as was the case for some ex-Rayon CFs or for stabilised PAN fibres chemically activated with different alkali hydroxides at much lower weight ratios and reaction times [87]. Here, only surface etching and smoothening were evidenced, which clearly indicates that LiOH alters the fibre morphology differently due to differences in chemical reactivity. As presented elsewhere [49,56], it was also observed that subsequent carbonisation and graphitisation treatments had an influence on the structural improvement of ex-PAN carbon fibres, which controlled morphological changes during air oxidation and chemical activation in
KOH or NaOH. Here, a clear effect of the chemical reactivity in molten LiOH was only observed at higher temperature, above 650 °C. Thus, the studied carbonised ex-PAN CFs (CP5a750 and CP7a750) experienced surface smoothening and a shrinkage of the carbon fibre diameter of 2 µm on average (see Figure 12). On the other hand, the outer surface of the graphitised ex-PAN carbon fibres showed modified grooves and striations, elongated parallel to the fibre axis direction, following that observed in the initial materials (see Figure S9). Again, an overall decrease in fibre diameter of 0.5 µm was observed, but much lower compared to carbonised CFs (by comparing the corresponding CP8a650 and GP8a650). The additional increase of test temperature widened the submicron surface pits, and a larger decrease in fibre diameter was observed due to the increased chemical reaction kinetics (as for GP8a750 and GP6a750).

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Same as Figure 11, but for ex-PAN CFs.</figcaption></figure>

The chemical stability of the investigated CFs, regardless of their originating precursor (ex-Rayon or ex-PAN), is directly related to the stability and persistence of their graphitic structure. Again, from the SEM analysis, the chemical stability of the CFs is identified when comparing materials of similar crystallites size ($L_c$) and interlayer spacing ($d_{002}$), but also similar carbon purity. Both GR4a750 and GF8a750 show only slight surface changes after the chemical reaction in molten LiOH, without damaging the overall structure of the constituent carbon fibres.

### 3.3.4. Chemical Stability of CFs in LiBr or $Li_4(OH)_3Br$ Compared to that in LiOH

The study of the chemical stability of CFs in molten LiBr or $Li_4(OH)_3Br$ was then carried out in the same way as for LiOH, and compared to it. The results of this section reveal that the nature of the inorganic lithium salts plays a major role in the chemical reactivity of CFs. Thus, in all examined samples (Table S1), a lower chemical reactivity was observed in molten $Li_4(OH)_3Br$ with respect to LiOH (visible graphically for the corresponding carbonised and graphitised CFs in Figure 13a with values of B-O 50% lower, or even less. As previously observed for LiOH, it was again found that the chemical reactivity in molten $Li_4(OH)_3Br$ is higher for ex-Rayon CFs than for ex-PAN CFs with a similar or lower heat-treatment temperature. Indeed, the main properties that regulate the chemical stability of CFs in molten $Li_4(OH)_3Br$ are found to be the purity of carbon and the structural order, in the same way as for LiOH.
<figure><img src="image_15.png" /><figcaption>**Figure 13.** (a) Burn-off (%) as a function of reaction temperature of ex-Rayon CFs; empty and full symbols represent CFs reacted in LiOH, and crossed or half-full symbols present CFs reacted in Li₄(OH)₃Br; the dashed or solid lines are only guides for the eye. (b) Upper figure: Theoretical phase diagram of a LiOH/LiBr binary system, with the corresponding peritectic transition (Tₚ) and liquidus temperatures (Tₘ); Bottom figure: Gravimetric energy density dependence of the LiBr/(LiBr + LiOH) ratio as seen from the calculated (yellow) line (after [30]).</figcaption></figure>

In addition, to study the role of the second component of the peritectic compound on the chemical reactivity of CFs, tests were carried out in molten LiBr only for the two CFs (carbonised CR1 and CR2) which proved to be the most reactive with LiOH. The chemical reactivity tests were performed at the highest temperature experienced (750 °C), all other conditions remaining unchanged.

By repeating the same experiment more than twice, no significant difference in weight before and after the performed treatment was observed. Thus, LiBr alone does not affect the carbon structure of the materials, and its chemical reactivity with the CFs is considered negligible. The chemical inertness of CFs in molten LiBr can be explained by the very high melting temperature of LiBr (about 550 °C) [34]. Moreover, the low reactivity of LiBr-KBr eutectic salts with carbons (or their initial precursors) was previously used to produce 3D carbon structures from biosourced materials under similar heat treatment conditions to those described herein [88]. Therefore, the obtained results indicate that the undesirable chemical reactivity of CFs with Li₄(OH)₃Br could only occur because of reactions between the carbon fibres and the hydroxide. The chemical reaction should start even at lower temperature, compared to LiOH alone, since the melting of Li₄(OH)₃Br (PC) begins with the formation of the pre-peritectic phase (EC, around 304 °C, see Figure 13b) [30]. The reaction should be further facilitated when the residual solid LiOH is melted at a higher temperature (about 380 °C) in stoichiometric proportion with LiBr, as shown in the theoretical phase diagram of Figure 13b.

The textural properties of the CFs also show less structural changes when submitted to molten Li₄(OH)₃Br, compared to molten LiOH (see Figures 9 and 10). Hence, as with the LiOH experiments, the surface area of CFs tested in molten Li₄(OH)₃Br tends to increase with the final test temperature, and is again highly dependent on the CFs structural properties. Although the surface area of the CFs has increased as was usually observed with LiOH, it generally comes from different micro/mesopore volume distributions. Therefore, when the CFs are tested in molten Li₄(OH)₃Br, the textures remain essentially microporous, whereas those of materials reacted with LiOH are much more mesoporous.
This is clearly visible in Figure 9a,c for samples CR1 and CR2 tested at 600 °C. Thus, the strong reaction in molten LiOH destroys most of the micropores, merging them with the developed high mesoporosity, which is observed in molten Li₄(OH)₃Br to a much lesser extent. When tests are carried out on CFs with higher carbon purity and structural order (GR3 and GR4), the porosity is less modified in Li₄(OH)₃Br with respect to LiOH (see Figure 10). In agreement, the PSDs are more modified for all the materials tested in molten LiOH compared to those in Li₄(OH)₃Br (see also Figures 9 and 10). The results with Li₄(OH)₃Br also suggest that only initial micropore widening seems to occur at lower test temperatures, but the development of larger mesopores is moderate, and observed only at higher test temperatures (at 600 °C, for carbonised CFs, and above for graphitised CFs).

The SEM micrographs of the tested CFs confirm that the modifications of the carbon fibres depend largely upon the presence and nature of the inorganic salt. Thus, all ex-Rayon CFs tested in molten Li₄(OH)₃Br retain the initial fibre structure, and only surface etching and pitting are observed, compared with the severe destruction of the fibres during the reaction with LiOH (see Figure 11). Uniform, smooth and flat surface carbon fibres without the initial surface striations, were observed after reaction of the carbonised ex-PAN CFs in molten Li₄(OH)₃Br (see Figure 12). Again, the lower B-O values observed suggest a smaller decrease in fibre diameter compared to the decrease seen in LiOH.

In general, the chemical reactions of CFs in the investigated molten binary peritectic compound, Li₄(OH)₃Br, should be considered similar to those in LiOH, but the presence of the LiBr phase only inhibits the severity of the redox reactions with the hydroxide. The reaction of the hydroxide with the active carbon sites should be inhibited only because of the presence of LiBr, no chemical reaction being observed between it and the carbon. The present LiBr phase does not completely prevent the reactions, but delays them differently depending on the structural properties of the CFs and the test temperatures. This is evident when one sees that the PSDs of the same CF tested in LiOH at a given temperature on the one hand, and in Li₄(OH)₃Br at a temperature of 50 to 100°C higher, are the same (see Figures 9 and 10). In particular, the consumption of LiOH during the reaction with the CFs shifts the molar ratio of lithium salts in the binary system from its stoichiometric composition (pure Li₄(OH)₃Br phase, presented by the blue dashed line in Figure 13b) to a higher concentration of LiBr phase. The continuous change in ratio and enrichment in the LiBr phase, as well as the production of Li₂CO₃ (from reaction (1)) must introduce multiple variations of the LiOH–LiBr binary system, ultimately leading to changes in the stoichiometry at the nanoscale.

As suggested by the adsorption analysis, the LiBr phase might therefore be trapped in the microporosity after the complete consumption of the LiOH phase, and the subsequent reactions of the latter would only cause a widening of the mesopores in the materials tested in Li₄(OH)₃Br.

Depending on the chemical stability of the CFs in the LiOH phase within Li₄(OH)₃Br, the relevant stoichiometric ratio for the latter PC can be modified to a different degree, thus jeopardising the high initial energy density (see Figure 13b). Therefore, choosing the most inert CFs should inhibit the chemical reaction in molten LiOH at the application temperature, thus ensuring one of the major advantages of using unmodified Li₄(OH)₃Br in its composition: its high energy-storage potential.

Our study also shows a beneficial reduction in the chemical reactivity of CFs with the LiOH phase in Li₄(OH)₃Br, when compared to LiOH alone. Moreover, the reduction in chemical reactivity can lead to an adjustable textural modification of some CFs with a poorly developed graphitic structure (such as CR1, GR1, CR2 and GR2). The use of such reactive phase/inert phase binary systems can play a dual role of pore modification (through carbon activation) and micropore protection (through the presence of the inert phase) resulting in well-developed micro/mesoporous carbon materials. Indeed, the results obtained here show a potential and beneficial use of the binary inorganic salt system for other fields of research and energy/environmental applications, such as the “salt-templated” carbonisation and synthesis of carbons of tuned porosity [88–90].
### 3.3.5. Estimation of CFs Chemical Stability at the TES Application Temperature

The last part of this study was performed to determine the activation energy of the reactions occurring with CFs in molten LiOH. The main idea of this part is to obtain more information about the results of B-O and to use them to estimate the B-O of the materials at the application temperature (between 300 and 400 °C). Therefore, it is expected that the apparent activation energy will change depending on the properties of the CFs and the chemical reactions that occur in the fixed experimental tests. As discussed in Section 3.3.1, the dominant chemical reaction (1) was expected within the limits of the experiments performed. However, it was previously clear that for some carbons of low structural order and low purity, the carbon/metal hydroxide reaction (1) can be followed by consecutive carbon oxidation reactions within the applied experimental limits [77]:

$$C + Li_2CO_3 \rightarrow Li_2O + 2CO \quad (Reaction 2)$$

Although $Li_2CO_3$ has a high melting temperature of about 730 °C [91], earlier decomposition can be observed at lower temperatures and even under a constant flow of inert gas, especially in contact with carbons of low graphitisation stage and bearing a high number of surface groups [92]. Therefore, reaction (2) is expected due to the decrease of the activation energy of $Li_2CO_3$ as a product of reaction (1), and its subsequent reaction with unreacted carbon. However, even if it occurs at about 730 °C or less, the rate of reaction (2) is initially low and quite negligible below 800 °C in the presence of a less reactive carbon material containing not more than 10% of foreign elements [93]. In our results, this is seen by the higher B-O and the dramatic change in the textural properties of CFs, with clear destruction of microporosity (as for CR1a600 and CR2a600 seen in Figure 9a,c), both generally seen above 600 °C. However, for all the CFs studied (under fixed test parameters, i.e., by changing only the material), an exponential increase of the B-O is evident by performing the experiment at a higher final temperature. These results clearly show that the kinetics of the chemical reactions that occur depend on the temperature, which makes it possible to use the Arrhenius law to estimate the apparent activation energy. The Arrhenius equation reads as follows:

$$k = Ae^{-E_a/RT} \quad (3)$$

where k is the reaction rate, T is the absolute temperature, A is a pre-exponential factor, $E_a$ is the apparent activation energy of the reaction, and R is the universal gas constant (8.314 J mol⁻¹ K⁻¹). As in other studies dealing with carbons oxidation reactions [94], in which the reaction rate constant was simply calculated as the inverse of the test temperature, we implicitly assume that first-order kinetics applies to most of the chemical reactivity process.

Although the present redox reactions can be expected to occur through an elementary mechanism, a combination of several mechanisms cannot be ruled out from the results obtained. However, previous studies of carbon reactions with oxygen or other catalytic agents [93,95] have shown that the pre-exponential factor and the apparent activation energy change in the same way, even when several reactions occur. This gives the possibility of using the Arrhenius law to obtain a consistent set of comparative values when considering different carbon materials. Therefore, the Arrhenius law is used to calculate the corresponding activation energy, thus determining the B-O of the CFs adjusted to a different temperature, by perfectly applying the linearised form of Equation (3), as follows:

$$\ln(k) = \frac{-E_a}{R} \left(\frac{1}{T}\right) + \ln(A) \quad (4)$$

The application of Equation (3) led determination coefficients (R²) higher than 0.919 (see Figure S10). It is found that the apparent activation energy is in the range of 116 to 165 kJ mol⁻¹, from carbonised to graphitised CFs, respectively. The carbonised CFs, manufactured at temperatures below 1200 °C, with a highly disordered structure, high porosity (surface area) and active surface (such as surface heteroatoms),
were more reactive and more prone to oxidation, even at lower temperatures. For thermodynamic reasons, the existence of an active carbon surface induces a higher number of reactive centres that decrease the activation energy of the reaction [95]. The relationship between all the CFs' characteristics controlling the chemical reaction justifies the reason for which different slopes were found when fitting the experimental results. Therefore, the apparent activation energy is lower for the carbonised CFs than for the graphitised CFs, with an improved structure and a lower number of surface defects and higher amount of heteroatom content. Thus, it is found that the estimation of activation energy values corresponds to the differences of CFs structure. This is observed when the activation energy values are plotted as a function of the interlayer spacing (d<sub>002</sub>) or the materials' elemental composition (see Figure 14a,b, respectively). The results follow a consistent order, with a minor scattering of the observed trends due to possible experimental or subsequent calculation deviations.

<figure><img src="image_16.png" /><figcaption>**Figure 14.** Correlation trends between estimated apparent activation energy and CFs properties, as a function of (a) interlayer spacing (d<sub>002</sub>) or (b) materials elemental composition. Carbonised or graphitised CFs are represented by empty or solid symbols, respectively.</figcaption></figure>

Furthermore, based on the performed bibliographic research, no activation energy values are reported for similar carbon materials tested under the current test conditions and in molten LiOH. However, the estimated values of apparent activation energy are in good agreement with others for carbon gasification in the presence of Li salts (such as Li<sub>2</sub>CO<sub>3</sub> and LiOH) under inert atmosphere or in the presence of oxygen, carbon monoxide/dioxide, or in steam [93,96–98]. In addition, Cuesta et al. (1993) [95] observed that ex-PAN high-purity carbon fibres react with oxygen at similar apparent activation energy values of between 117 and 134 kJ mol<sup>-1</sup>.

Thus, the slightly lower activation energy of high-purity carbon fibres, compared to the present results, was evident because it is known that the reaction of oxygen with carbon is quite exothermic, which decreases the activation energy, unlike CO<sub>2</sub>, steam or other chemical activating agents [99]. Once again, in all the studies reviewed, it was considered that the structural characteristics of the carbons play an important role in the trends in chemical reactivity. Therefore, the calculated parameters related to the chemical reactions that occur in the CFs' stability tests seem to correlate well with the bibliography. Hence, the subsequent use of the derived experimental results is considered valid for estimating the level of B-O, due to an undesirable chemical reaction between the CF hosts and LiOH at the application temperature (presented in Figure 15).
<figure><img src="image_17.png" /><figcaption>**Figure 15.** Estimated trends in B-O from the Equation (4), for different CFs, in the application temperatures of interest. Carbonised or graphitised CFs are presented with dashed (ex-Rayon)/dotted (ex-PAN) or solid trend lines, respectively.</figcaption></figure>

The set of CFs analysed is extremely diverse, so that materials with different structural properties allow meaningful estimates of reactivity at the application temperature (see Figure 15). The estimated trends of B-O are in agreement with the experimental results, distinguishing the graphitised CFs (i.e., with an improved purity and structure, e.g. ex-Rayon, GR4, and ex-PAN, GP6 and GP8), from the rest. In addition, we must consider that these estimates are based on the results of B-O only with LiOH and that a larger decrease of reactivity is expected with $Li_4(OH)_3Br$, as discussed in Section 3.3.4. Therefore, only those candidates that are chemically most stable and found to be suitable will be considered for the preparation and testing of future CF-PC hybrid materials.

## 4. Conclusions

In this work, we discussed how the structural and textural parameters of commercial fibrous carbons influence their chemical stability in Li- based inorganic salts and promote their use as carbon hosts in TES application. We focused on commercial CFs as mass-produced materials using well-known techniques, considered here as cost-/performance-effective solutions for hosting phase-change materials. The application of different characterisation techniques, including elemental analysis, XRD, Raman spectroscopy, gas adsorption and SEM, allowed us to evaluate clearly the range of properties of the examined CFs, which affected their stability during the chemical reaction tests.

For both ex-PAN and ex-Rayon materials, as carbon purity and structural order increase, the presence of nanopores is reduced, and CFs become chemically stable and fully retain their fibre morphology when in contact with molten LiOH, even at the highest temperature tested. The chemical stability tests confirmed the effect of temperature by increasing the B-O due to the increase in the chemical reaction rate.

Thus, the main reasons for the observed carbon loss are indicated as being the high LiOH/CF ratio, and the reaction of LiOH with the surface of carbon fibres under the used experimental conditions.

The effect of the salts' nature was evaluated by performing additional CFs' chemical stability tests in LiBr or $Li_4(OH)_3Br$. The results show the complete CFs' inertness in molten LiBr, and increased stability in molten $Li_4(OH)_3Br$, compared to the higher chemical reactivity in molten LiOH alone. In general, the chemical reactions of CFs in $Li_4(OH)_3Br$ are considered similar to those in LiOH alone, but the presence of an inert LiBr phase plays a beneficial role in inhibiting the severity of redox reactions with the hydroxide. Finally, it was confirmed that the experimental results of chemical stability tests are in excellent agreement with the values of activation energy estimated by use of the Arrhenius law, distinguishing again graphitised CFs, with improved purity and structure, from the rest.
Finally, the criteria for selecting potential candidates for the TES application based on the peritectic compound Li₄(OH)₃Br hosted in fibrous carbon matrices, i.e., CFs with high chemical stability in contact with Li salts, may be obtained from the information presented in this study. In addition, the identification of the critical properties of CFs is a useful guide to further improve materials, but also the use of the CFs and Li₄(OH)₃Br in different applications.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/1996-1944/12/24/4232/s1, Section S1: Raman spectroscopy analysis; Figure S1: Raman laser spot position (yellow midpoint), for achieving Raman spectra for both cross-section (left) and surface (right micrographs) of ex-PAN (CP8: top row) and ex-Rayon (GR4: bottom row) carbon fibres; Figure S2: Deconvoluted Raman spectra of ex-Rayon CFs, where the scatter curves indicate the measured intensities and the solid smooth lines represent the fits; Figure S3: Same as Figure S2 but for ex-PAN CFs; Figure S4: Intensity ratio between deconvoluted D₁ and G bands of corresponding: (a) ex-Rayon and (b) ex-PAN CFs, with an arrow indicating the increase of structural order. Empty and full symbols stand for carbonised CFs and graphitised CFs, respectively; Section S2: Textural analysis; Figure S5: (a,b) N₂ adsorption isotherms at -196 °C, and (c,d) CO₂ adsorption isotherms at 0 °C for de-sized CFs before reaction in molten Li-based salts; Figure S6: N₂ adsorption isotherms of corresponding carbonised and graphitised CFs before and after chemical reaction with LiOH (black symbols) and with Li₄(OH)₃Br (green symbols). The arrow is shown for CFs reacted at a same temperature but in different molten salts; Figure S7: Same as Figure S6 but for graphitised ex-Rayon CFs (of highest structural order); Table S1: Textural parameters from both N₂ and CO₂ isotherms calculated by 2D NLDFT-HS model of ex-Rayon CFs examined before and after chemical reaction in molten LiOH (normal style) or Li₄(OH)₃Br (underlined style). The bold style distinguishes graphitised CFs from carbonised CFs; Table S2: Same as Table S1 but for ex-PAN CFs; Section S3: Scanning electron microscopy; Figure S8: SEM micrographs of de-sized ex-Rayon CFs before (left) and after (middle and right) reaction with LiOH at different temperatures; Figure S9: Same as Figure S8 but for de-sized ex-PAN CFs; Section S4: Estimation of CFs chemical stability at the TES application temperature; Figure S10: Application of Equation (4) to the experimentally derived B-O data of most CFs shown in Figure 7a,b). Carbonised or graphitised CFs are represented by empty or solid symbols and dashed or solid linear fits, respectively.

**Author Contributions:** Conceptualisation and methodology, B.K., V.F. and A.C.; bibliographic investigation, B.K.; elemental analysis, P.G.; XRD, G.M.; SEM, S.M.; other experiments, B.K., V.F. and A.C.; data analysis, all; data curation, B.K.; funding acquisition and project administration, V.F. and A.C.; writing—original draft preparation, B.K.; writing—review and editing, V.F. and A.C.; supervision, V.F. and A.C.

**Funding:** This work was partly supported by a grant overseen by the French National Research Agency (Pc2TES ANR-16-CE06-0012-01). The authors also gratefully acknowledge the financial support of European FEDER funds (TALiSMAN project, 2019-000214).

**Acknowledgments:** The authors are thankful to all suppliers having kindly delivered free samples for our research: Sean McFeely from CeraMaterials (UK), Zhao (Johnson) Xi from Beijing Great Wall Co. Ltd (China), Alain Pichon from SGL CARBON S.A.S. (France), Francis Guérard and Elena Ogneva from Schunk Electrographite (France), and István Fáykiss and Kinga Kincses from Zoltek (Hungary).

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Nejat, P.; Jomehzadeh, F.; Taheri, M.M.; Gohari, M.; Majid, M.Z.A. A global review of energy consumption, CO₂ emissions and policy in the residential sector (with an overview of the top ten CO₂ emitting countries). *Renew. Sustain. Energy Rev.* **2015**, *43*, 843–862. [CrossRef]

2. Dincer, I.; Acar, C. Chapter 1.1—Potential Energy Solutions for Better Sustainability. In *Exergetic, Energetic and Environmental Dimensions*; Dincer, I., Colpan, C.O., Kizilkan, O., Eds.; Academic Press: New York, NY, USA, 2018; pp. 3–37. ISBN 978-0-12-813734-5.

3. Gielen, D.; Boshell, F.; Saygin, D.; Bazilian, M.D.; Wagner, N.; Gorini, R. The role of renewable energy in the global energy transformation. *Energy Strategy Rev.* **2019**, *24*, 38–50. [CrossRef]

4. *Global Market Outlook for Solar Power 2018–2022*; SolarPower: Brussels, Belgium, 2018.

5. Calderón, A.; Barreneche, C.; Hernández-Valle, K.; Galindo, E.; Segarra, M.; Fernández, A.I. Where is Thermal Energy Storage (TES) research going?—A bibliometric analysis. *Sol. Energy* **2019**. [CrossRef]

6. Palomba, V.; Frazzica, A. Recent advancements in sorption technology for solar thermal energy storage applications. *Sol. Energy* **2018**. [CrossRef]
7. Pelay, U.; Luo, L.; Fan, Y.; Stitou, D.; Rood, M. Thermal energy storage systems for concentrated solar power plants. *Renew. Sustain. Energy Rev.* **2017**, *79*, 82–100. [CrossRef]

8. Ibrahim, N.I.; Al-Sulaiman, F.A.; Rahman, S.; Yilbas, B.S.; Sahin, A.Z. Heat transfer enhancement of phase change materials for thermal energy storage applications: A critical review. *Renew. Sustain. Energy Rev.* **2017**, *74*, 26–50. [CrossRef]

9. Liu, M.; Steven Tay, N.H.; Bell, S.; Belusko, M.; Jacob, R.; Will, G.; Saman, W.; Bruno, F. Review on concentrating solar power plants and new developments in high temperature thermal energy storage technologies. *Renew. Sustain. Energy Rev.* **2016**, *53*, 1411–1432. [CrossRef]

10. Zhou, C.; Wu, S. Medium-and high-temperature latent heat thermal energy storage: Material database, system review, and corrosivity assessment. *Int. J. Energy Res.* **2019**, *43*, 621–661. [CrossRef]

11. Pandey, A.K.; Hossain, M.S.; Tyagi, V.V.; Abd Rahim, N.; Selvaraj, J.A.L.; Sari, A. Novel approaches and recent developments on potential applications of phase change materials in solar energy. *Renew. Sustain. Energy Rev.* **2018**, *82*, 281–323. [CrossRef]

12. Achchaq, F.; Barrio, E.P.D. A proposition of peritectic structures as candidates for thermal energy storage. *Energy Procedia* **2017**, *139*, 346–351. [CrossRef]

13. Achchaq, F.; Palomo del Barrio, E.; Renauld, A.; Ben-Khemis, S. Characterization of Li₂K(OH)₃ as material for thermal energy storage at high temperature. In Proceedings of the 13th International Conference on Energy Storage, Beijing, China, 19–21 May 2015; Volume 2, pp. 1–8.

14. Badenhorst, H. A review of the application of carbon materials in solar thermal energy storage. *Sol. Energy* **2019**, *192*, 35–68. [CrossRef]

15. Nomura, T.; Tabuchi, K.; Zhu, C.; Sheng, N.; Wang, S.; Akiyama, T. High thermal conductivity phase change composite with percolating carbon fiber network. *Appl. Energy* **2015**, *154*, 678–685. [CrossRef]

16. Jiang, Z.; Ouyang, T.; Yang, Y.; Chen, L.; Fan, X.; Chen, Y.; Li, W.; Fei, Y. Thermal conductivity enhancement of phase change materials with form-stable carbon bonded carbon fiber network. *Mater. Des.* **2018**, *143*, 177–184. [CrossRef]

17. Fukai, J.; Kanou, M.; Kodama, Y.; Miyatake, O. Thermal conductivity enhancement of energy storage media using carbon fibers. *Energy Convers. Manag.* **2000**, *41*, 1543–1556. [CrossRef]

18. Hamada, Y.; Otsu, W.; Fukai, J.; Morozumi, Y.; Miyatake, O. Anisotropic heat transfer in composites based on high-thermal conductive carbon fibers. *Energy* **2005**, *30*, 221–233. [CrossRef]

19. Frusteri, F.; Leonardi, V.; Vasta, S.; Restuccia, G. Thermal conductivity measurement of a PCM based storage system containing carbon fibers. *Appl. Therm. Eng.* **2005**, *25*, 1623–1633. [CrossRef]

20. Karaipekli, A.; Sari, A.; Kaygusuz, K. Thermal conductivity improvement of stearic acid using expanded graphite and carbon fiber for energy storage applications. *Renew. Energy* **2007**, *32*, 2201–2210. [CrossRef]

21. Wang, J.; Xie, H.; Xin, Z.; Li, Y.; Yin, C. Investigation on thermal properties of heat storage composites containing carbon fibers. *J. Appl. Phys.* **2011**, *110*, 094302. [CrossRef]

22. Samimi, F.; Babapoor, A.; Azizi, M.; Karimi, G. Thermal management of a Li-ion battery using carbon fiber-PCM composites. *Appl. Therm. Eng.* **2015**, *82*, 281–290. [CrossRef]

23. Huang, X.; Alva, G.; Liu, L.; Fang, G. Microstructure and thermal properties of cetyl alcohol/high density polyethylene composite phase change materials with carbon fiber as shape-stabilized thermal storage materials. *Appl. Energy* **2017**, *200*, 19–27. [CrossRef]

24. Zhang, Q.; Luo, Z.; Guo, Q.; Wu, G. Preparation and thermal properties of short carbon fibers/erythritol phase change materials. *Energy Convers. Manag.* **2017**, *136*, 220–228. [CrossRef]

25. Ye, F.; Ge, Z.; Ding, Y.; Yang, J. Multi-walled carbon nanotubes added to Na₂CO₃/MgO composites for thermal energy storage. *Particuology* **2014**, *15*, 56–60. [CrossRef]

26. Singh, D.; Kim, T.; Zhao, W.; Yu, W.; France, D.M. Development of graphite foam infiltrated with MgCl₂ for a latent heat based thermal energy storage (LHTES) system. *Renew. Energy* **2016**, *94*, 660–667. [CrossRef]

27. Guillot, S.; Faik, A.; Rakhmatullin, A.; Lambert, J.; Veron, E.; Echegut, P.; Bessada, C.; Calvet, N.; Py, X. Corrosion effects between molten salts and thermal storage material for concentrated solar power plants. *Appl. Energy* **2012**, *94*, 174–181. [CrossRef]
28. Calvet, N.; Gomez, J.C.; Faik, A.; Roddatis, V.V.; Meffre, A.; Glatzmaier, G.C.; Doppiu, S.; Py, X. Compatibility of a post-industrial ceramic with nitrate molten salts for use as filler material in a thermocline storage system. *Appl. Energy* **2013**, *109*, 387–393. [CrossRef]

29. Vasu, A.; Hagos, F.Y.; Noor, M.M.; Mamat, R.; Azmi, W.H.; Abdullah, A.A.; Ibrahim, T.K. Corrosion effect of phase change materials in solar thermal energy storage application. *Renew. Sustain. Energy Rev.* **2017**, *76*, 19–33. [CrossRef]

30. Achchaq, F.; Palomo del Barrio, E.; Lebraud, E.; Péchev, S.; Toutain, J. Development of a new LiBr/LiOH-based alloy for thermal energy storage. *J. Phys. Chem. Solids* **2019**, *131*, 173–179. [CrossRef]

31. Karakashov, B.; Toutain, J.; Achchaq, F.; Legros, P.; Fierro, V.; Celzard, A. Permeability of fibrous carbon materials. *J. Mater. Sci.* **2019**, *54*, 13537–13556. [CrossRef]

32. Naito, K.; Yang, J.-M.; Xu, Y.; Kagawa, Y. Enhancing the thermal conductivity of polyacrylonitrile- and pitch-based carbon fibers by grafting carbon nanotubes on them. *Carbon* **2010**, *6*, 1849–1857. [CrossRef]

33. Szabó, L.; Imanishi, S.; Tetsuo, F.; Hirose, D.; Ueda, H.; Tsukegi, T.; Ninomiya, K.; Takahashi, K. Lignin as a Functional Green Coating on Carbon Fiber Surface to Improve Interfacial Adhesion in Carbon Fiber Reinforced Polymers. *Materials* **2019**, *12*, 159. [CrossRef]

34. Wietelmann, U.; Bauer, R.J. Lithium and Lithium Compounds. *Ullmanns Encycl. Ind. Chem.* **2000**, *21*, 339–363. [CrossRef]

35. Linares-Solano, A.; Lozano-Castello, D.; Lillo-Rodenas, M.; Cazorla-Amoros, D. Carbon Activation by Alkaline Hydroxides. In *Chemistry and Physics of Carbon*; Radovic, L., Ed.; CRC Press: New York, NY, USA, 2007; pp. 1–62.

36. Achchaq, F.; Risueño, E.; Mahroug, I.; Legros, P.; Lebraud, E.; Karakashov, B.; Palomo del Barrio, E.; Celzard, A.; Fierro, V.; Toutain, J. Development of a Carbon Felt/Salt-Based Hybrid Material for Thermal Energy Storage Applications. *J. Energy Power Eng.* **2018**, *12*. [CrossRef]

37. Pope, C.G. X-Ray Diffraction and the Bragg Equation. *J. Chem. Educ.* **1997**, *74*, 129. [CrossRef]

38. Celzard, A.; Pasc, A.; Schaefer, S.; Mandel, K.; Ballweg, T.; Li, S.; Medjahdi, G.; Nicolas, V.; Fierro, V. Floating hollow carbon spheres for improved solar evaporation. *Carbon* **2019**, *146*, 232–247. [CrossRef]

39. Pimenta, M.A.; Dresselhaus, G.; Dresselhaus, M.S.; Cançado, L.G.; Jorio, A.; Saito, R. Studying disorder in graphite-based systems by Raman spectroscopy. *Phys. Chem. Chem. Phys.* **2007**, *9*, 1276–1290. [CrossRef]

40. Sadezky, A.; Muckenhuber, H.; Grothe, H.; Niessner, R.; Pöschl, U. Raman microspectroscopy of soot and related carbonaceous materials: Spectral analysis and structural information. *Carbon* **2005**, *43*, 1731–1742. [CrossRef]

41. Manciu, F.S.; Oh, Y.; Barath, A.; Rusheen, A.E.; Kouzani, A.Z.; Hodges, D.; Guerrero, J.; Tomshine, J.; Lee, K.H.; Bennet, K.E. Analysis of Carbon-Based Microelectrodes for Neurochemical Sensing. *Materials* **2019**, *12*, 3186. [CrossRef]

42. Chen, J. (Ed.) *Activated Carbon Fiber and Textiles*, 1st ed.; Woodhead Publishing: Cambridge, UK, 2016.

43. Chiang, Y.-C.; Lee, C.-Y.; Lee, H.-C. Surface chemistry of polyacrylonitrile- and rayon-based activated carbon fibers after post-heat treatment. *Mater. Chem. Phys.* **2007**, *101*, 199–210. [CrossRef]

44. Rahaman, M.S.A.; Ismail, A.F.; Mustafa, A. A review of heat treatment on polyacrylonitrile fiber. *Polym. Degrad. Stab.* **2007**, *92*, 1421–1432. [CrossRef]

45. Nunna, S.; Maghe, M.; Rana, R.; Varley, R.J.; Knorr, D.B.; Sands, J.M.; Creighton, C.; Henderson, L.C.; Naebe, M. Time Dependent Structure and Property Evolution in Fibres during Continuous Carbon Fibre Manufacturing. *Materials* **2019**, *12*, 1069. [CrossRef]

46. Tang, M.M.; Bacon, R. Carbonization of cellulose fibers—I. Low temperature pyrolysis. *Carbon* **1964**, *2*, 211–220. [CrossRef]

47. Park, S.-J.; Kim, B.-J. Carbon Fibers and Their Composites. In *Carbon Fibers*; Springer Series in Materials Science; Springer: Dordrecht, The Netherlands, 2015; pp. 275–317. ISBN 978-94-017-9477-0.

48. Al Aiti, M.; Jehnichen, D.; Fischer, D.; Brünig, H.; Heinrich, G. On the morphology and structure formation of carbon fibers from polymer precursor systems. *Prog. Mater. Sci.* **2018**, *98*, 477–551. [CrossRef]

49. Li, C.; Xian, G. Experimental and Modeling Study of the Evolution of Mechanical Properties of PAN-Based Carbon Fibers at Elevated Temperatures. *Materials* **2019**, *12*, 724. [CrossRef] [PubMed]
50. Fink, H.-P.; Hofmann, D.; Philipp, B. Some aspects of lateral chain order in cellulosics from X-ray scattering. *Cellulose* **1995**, *2*, 51–70. [CrossRef]

51. Wada, M.; Okano, T.; Sugiyama, J. Allomorphs of native crystalline cellulose I evaluated by two equatorial spacings. *J. Wood Sci.* **2001**, *47*, 124–128. [CrossRef]

52. Kim, D.-Y.; Nishiyama, Y.; Wada, M.; Kuga, S. Graphitization of highly crystalline cellulose. *Carbon* **2001**, *39*, 1051–1056. [CrossRef]

53. Dumanli, A.G.; Windle, A.H. Carbon fibres from cellulosic precursors: A review. *J. Mater. Sci.* **2012**, *47*, 4236–4250. [CrossRef]

54. Byrne, N.; Chen, J.; Fox, B. Enhancing the carbon yield of cellulose based carbon fibres with ionic liquid impregnates. *J. Mater. Chem. A* **2014**, *2*, 15758–15762. [CrossRef]

55. Carrott, P.J.M.; Nabais, J.M.V.; Ribeiro Carrott, M.M.L.; Pajares, J.A. Preparation of activated carbon fibres from acrylic textile fibres. *Carbon* **2001**, *39*, 1543–1555. [CrossRef]

56. Maciá-Agulló, J.A.; Moore, B.C.; Cazorla-Amorós, D.; Linares-Solano, A. Influence of carbon fibres crystallinities on their chemical activation by KOH and NaOH. *Microporous Mesoporous Mater.* **2007**, *101*, 397–405. [CrossRef]

57. Dobiášová, L.; Starý, V.; Glogar, P.; Valvoda, V. Analysis of carbon fibers and carbon composites by asymmetric X-ray diffraction technique. *Carbon* **1999**, *37*, 421–425. [CrossRef]

58. Lu, L.; Sahajwalla, V.; Kong, C.; Harris, D. Quantitative X-ray diffraction analysis and its application to various coals. *Carbon* **2001**, *39*, 1821–1833. [CrossRef]

59. Sonibare, O.O.; Haeger, T.; Foley, S.F. Structural characterization of Nigerian coals by X-ray diffraction, Raman and FTIR spectroscopy. *Energy* **2010**, *35*, 5347–5353. [CrossRef]

60. Kong, K.; Deng, L.; Kinloch, I.A.; Young, R.J.; Eichhorn, S.J. Production of carbon fibres from a pyrolysed and graphitised liquid crystalline cellulose fibre precursor. *J. Mater. Sci.* **2012**, *47*, 5402–5410. [CrossRef]

61. Cuesta, A.; Dhamelincourt, P.; Laureyns, J.; Martínez-Alonso, A.; Tascón, J.M.D. Raman microprobe studies on carbon materials. *Carbon* **1994**, *32*, 1523–1532. [CrossRef]

62. Zickler, G.A.; Smarsly, B.; Gierlinger, N.; Peterlik, H.; Paris, O. A reconsideration of the relationship between the crystallite size $L_a$ of carbons determined by X-ray diffraction and Raman spectroscopy. *Carbon* **2006**, *44*, 3239–3246. [CrossRef]

63. Zabihi, O.; Shafei, S.; Fakhrhoseini, S.M.; Ahmadi, M.; Ajdari Nazarloo, H.; Stanger, R.; Anh Tran, Q.; Lucas, J.; Wall, T.; Naebe, M. Low-Cost Carbon Fibre Derived from Sustainable Coal Tar Pitch and Polyacrylonitrile: Fabrication and Characterisation. *Materials* **2019**, *12*, 1281. [CrossRef]

64. Dresselhaus, M.S.; Dresselhaus, G.; Sugihara, K.; Spain, I.L.; Goldberg, H.A. (Eds.) *Graphite Fibers and Filaments*; Springer Series in Materials Science; Springer: Berlin/Heidelberg, Germany, 1988.

65. Zhong, S.; Padeste, C.; Kazacos, M.; Skyllas-Kazacos, M. Comparison of the physical, chemical and electrochemical properties of rayon-and polyacrylonitrile-based graphite felt electrodes. *J. Power Sources* **1993**, *45*, 29–41. [CrossRef]

66. Ko, T.-H.; Liao, Y.-K.; Liu, C.-H. Effects of graphitization of PAN-based carbon fiber cloth on its use as gas diffusion layers in proton exchange membrane fuel cells. *New Carbon Mater.* **2007**, *22*, 97–101. [CrossRef]

67. Dai, Z.; Shi, F.; Zhang, B.; Li, M.; Zhang, Z. Effect of sizing on carbon fiber surface properties and fibers/epoxy interfacial adhesion. *Appl. Surf. Sci.* **2011**, *257*, 6980–6985. [CrossRef]

68. He, Z.; Shi, L.; Shen, J.; He, Z.; Liu, S. Effects of nitrogen doping on the electrochemical performance of graphite felts for vanadium redox flow batteries. *Int. J. Energy Res.* **2015**, *39*, 709–716. [CrossRef]

69. Qian, X.; Zhang, Y.G.; Wang, X.F.; Heng, Y.J.; Zhi, J.H. Effect of carbon fiber surface functionality on the moisture absorption behavior of carbon fiber/epoxy resin composites. *Surf. Interface Anal.* **2016**, *48*, 1271–1277. [CrossRef]

70. Pradere, C.; Batsale, J.C.; Goyhénèche, J.M.; Pailler, R.; Dilhaire, S. Thermal properties of carbon fibers at very high temperature. *Carbon* **2009**, *47*, 737–743. [CrossRef]

71. Qiu, L.; Zheng, X.H.; Zhu, J.; Su, G.P.; Tang, D.W. The effect of grain size on the lattice thermal conductivity of an individual polyacrylonitrile-based carbon fiber. *Carbon* **2013**, *51*, 265–273. [CrossRef]

72. Harris, P.J.F. Rosalind Franklin’s work on coal, carbon, and graphite. *Interdiscip. Sci. Rev.* **2001**, *26*, 204–210. [CrossRef]
73. Jagiello, J.; Olivier, J.P. 2D-NLDFT adsorption models for carbon slit-shaped pores with surface energetical heterogeneity and geometrical corrugation. *Carbon* **2013**, *55*, 70–80. [CrossRef]

74. Girón, R.P.; Gil, R.R.; Suárez-Ruiz, I.; Fuente, E.; Ruiz, B. Adsorbents/catalysts from forest biomass fly ash. Influence of alkaline activating agent. *Microporous Mesoporous Mater.* **2015**, *209*, 45–53. [CrossRef]

75. Mikova, N.M.; Chesnokov, N.V.; Kuznetsov, B.N. Study of High Porous Carbons Prepared by the Alkaline Activation of Anthracites. *J. Sib. Fed. Univ. Chem.* **2009**, *1*, 3–10.

76. Mora, E.; Blanco, C.; Pajares, J.A.; Santamaría, R.; Menéndez, R. Chemical activation of carbon mesophase pitches. *J. Colloid Interface Sci.* **2006**, *298*, 341–347. [CrossRef]

77. Sutcu, H.; Dural, A. Effect of Hydroxides on Carbonization of Bituminous Coal. *Coal Prep.* **2006**, *26*, 201–208. [CrossRef]

78. Aweed, K.A. Preparation of Activated Carbon from Wastes Tires by Minerals Hydroxide (LiOH, NaOH and KOH). *Tikrit J. Pure Sci.* **2008**, *13*, 71–74.

79. Tamarkina, Y.V.; Shendrik, T.G.; Kucherenko, V.A.; Khabarova, T.V. Conversion of Alexandriya brown coal into microporous carbons under alkali activation. *Chemistry* **2012**, *5*, 24–36.

80. Chowdhury, Z.Z.; Hamid, S.B.A.; Das, R.; Hasan, M.R.; Zain, S.M.; Khalid, K.; Uddin, M.N. Preparation of Carbonaceous Adsorbents from Lignocellulosic Biomass and Their Use in Removal of Contaminants from Aqueous Solution. *BioResources* **2013**, *8*, 6523–6555. [CrossRef]

81. Cabasso, I.; Yuan, Y. *Nano Structured Activated Carbon for Hydrogen Storage*; Project Final Technical Report (May 2, 2005-Dec. 31, 2012); Research Foundation of State University of New York, Polymer Research Institute, State University of New York-esf: New York, NY, USA, 2013. [CrossRef]

82. Kucherenko, V.A.; Tamarkina, Y.V.; Raenko, G.F.; Chernyshova, M.I. Thermolysis of brown coal in the presence of alkali metal hydroxides. *Solid Fuel Chem.* **2017**, *51*, 147–154. [CrossRef]

83. Krstić, S.S.; Kragović, M.M.; Dodevski, V.M.; Marinković, A.D.; Kaluđerović, B.V.; Žerjav, G.; Pintar, A.; Pagnacco, M.C.; Stojmenović, M.D. Influence of temperature and different hydroxides on properties of activated carbon prepared from saccharose. Characterization, thermal degradation kinetic and dyes removal from water solutions. *Sci. Sinter.* **2018**, *50*, 255–273. [CrossRef]

84. Zhou, J.; Li, Z.; Xing, W.; Shen, H.; Bi, X.; Zhu, T.; Qiu, Z.; Zhuo, S. A New Approach to Tuning Carbon Ultramicropore Size at Sub-Angstrom Level for Maximizing Specific Capacitance and CO₂ Uptake. *Adv. Funct. Mater.* **2016**, *26*, 7955–7964. [CrossRef]

85. Celzard, A.; Fierro, V.; Marêché, J.F.; Furdin, G. Advanced Preparative Strategies for Activated Carbons Designed for the Adsorptive Storage of Hydrogen. *Adsorpt. Sci. Technol.* **2007**, *25*, 129–142. [CrossRef]

86. Tamarkina, Y.V.; Maslova, L.A.; Khabarova, T.V.; Kucherenko, V.A. Adsorption properties of carbon materials produced by thermolysis of brown coal in the presence of alkali metal hydroxides. *Russ. J. Appl. Chem.* **2008**, *81*, 1167–1170. [CrossRef]

87. Lee, Y.J.; Kim, J.H.; Kim, J.S.; Lee, D.B.; Lee, J.C.; Chung, Y.J.; Lim, Y.S. Fabrication of Activated Carbon Fibers from Stabilized PAN-Based Fibers by KOH. *Mater. Sci. Forum* **2004**, *449*, 217–220. [CrossRef]

88. Vilian, A.T.E.; Song, J.Y.; Lee, Y.S.; Hwang, S.-K.; Kim, H.J.; Jun, Y.-S.; Huh, Y.S.; Han, Y.-K. Salt-templated three-dimensional porous carbon for electrochemical determination of gallic acid. *Biosens. Bioelectron.* **2018**, *117*, 597–604. [CrossRef]

89. Silvestre-Albero, A.; Silvestre-Albero, J.; Martinez-Escandell, M.; Molina-Sabio, M.; Kovacs, A.; Rodriguez-Reinoso, F. Novel synthesis of a micro-mesoporous nitrogen-doped nanostructured carbon from polyaniline. *Microporous Mesoporous Mater.* **2015**, *218*, 199–205. [CrossRef]

90. Pampel, J.; Fellinger, T.-P. Opening of Bottleneck Pores for the Improvement of Nitrogen Doped Carbon Electrocatalysts. *Adv. Energy Mater.* **2016**, *6*. [CrossRef]

91. Schiemann, M.; Bergthorson, J.; Fischer, P.; Scherer, V.; Taroata, D.; Schmid, G. A review on lithium combustion. *Appl. Energy* **2016**, *162*, 948–965. [CrossRef]

92. Epstein, M.; Yogev, A.; Yao, C.; Berman, A. Carbothermal reduction of alkali hydroxides using concentrated solar energy. *Energy* **2001**, *26*, 441–455. [CrossRef]

93. McKee, D.W.; Chatterji, D. The catalytic behavior of alkali metal carbonates and oxides in graphite oxidation reactions. *Carbon* **1975**, *13*, 381–390. [CrossRef]

94. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. *Anal. Chem.* **1957**, *29*, 1702–1706. [CrossRef]
95. Cuesta, A.; Martinez-Alonso, A.; Tascon, J.M.D. Correlation between Arrhenius kinetic parameters in the reaction of different carbon materials with oxygen. *Energy Fuels* **1993**, *7*, 1141–1145. [CrossRef]

96. McKee, D.W.; Chatterji, D. The catalyzed reaction of graphite with water vapor. *Carbon* **1978**, *16*, 53–57. [CrossRef]

97. Kapteijn, F.; Abbel, G.; Moulijn, J.A. CO₂ gasification of carbon catalysed by alkali metals: Reactivity and mechanism. *Fuel* **1984**, *63*, 1036–1042. [CrossRef]

98. Sheth, A.; Agrawal, P.; Yeboah, Y.D. *Catalytic Gasification of Coal using Eutectic Salt Mixtures*; Federal Energy Technology Center: Morgantown, WV, USA; Pittsburgh, PA, USA, 1998. [CrossRef]

99. Sevilla, M.; Mokaya, R. Energy storage applications of activated carbons: Supercapacitors and hydrogen storage. *Energy Environ. Sci.* **2014**, *7*, 1250–1280. [CrossRef]

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

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).