Next Article in Journal
Muscular Oxygen Saturation and Hemoglobin Concentration during Freediving: A Case Study
Previous Article in Journal
Evaluation of Biochemical and Oxidative Stress Markers in the Early Stages of Rheumatoid Arthritis in a Comparative Study of Two Different Therapeutic Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Does the Oxygen Functionality Really Improve the Thermodynamics of Reversible Hydrogen Storage with Liquid Organic Hydrogen Carriers?

by
Sergey P. Verevkin
1,2,*,
Artemiy A. Samarov
3 and
Sergey V. Vostrikov
4
1
Competence Centre CALOR, Faculty of Interdisciplinary Research, University of Rostock, 18059 Rostock, Germany
2
Department of Physical Chemistry, Kazan Federal University, 420008 Kazan, Russia
3
Department of Chemical Thermodynamics and Kinetics, Saint Petersburg State University, Peterhof, 198504 Saint Petersburg, Russia
4
Chemical-Technological Department, Samara State Technical University, 443100 Samara, Russia
*
Author to whom correspondence should be addressed.
Oxygen 2024, 4(3), 266-284; https://doi.org/10.3390/oxygen4030015
Submission received: 30 May 2024 / Revised: 25 June 2024 / Accepted: 27 June 2024 / Published: 2 July 2024

Abstract

:
Liquid organic hydrogen carriers (LOHCs) are aromatic molecules that are being considered for the safe storage and release of hydrogen. The thermodynamic properties of a range of aromatic ethers were investigated using various experimental and theoretical methods to assess their suitability as LOHC materials. The absolute vapour pressures were measured for benzyl phenyl ether, dibenzyl ether and 2-methoxynaphthalene using the transpiration method. The standard molar enthalpies and entropies of vaporisation/sublimation were derived from the temperature dependence of the vapour pressures. The combustion energies of benzyl phenyl ether and dibenzyl ether were measured using high-precision combustion calorimetry, and their standard molar enthalpies of formation were derived from these data. High-level quantum chemical calculations were used to calculate the standard molar enthalpies of formation in the gas phase for benzyl phenyl ether, dibenzyl ether and 2-methoxynaphthalene. The latter values agreed very well with the experimental results obtained in this work. The thermodynamic properties of the hydrogenation/dehydrogenation reactions in liquid phase in LOHC systems based on methoxy–benzene, diphenyl ether, benzyl phenyl ether, dibenzyl ether and 2-methoxynaphthalene were derived and compared with the data for similarly structured hydrogen carriers based on benzene, diphenylmethane, 1,2-diphenylethane, 1,3-diphenylpropane and naphthalene. The influence of the oxygen functionality on the thermodynamic properties of the hydrogenation/dehydrogenation reactions was evaluated.

1. Introduction

Hydrogen has proven to be a promising alternative to conventional fossil fuels and has the potential to solve the environmental and energy problems of the 21st century. As a clean and versatile energy carrier, hydrogen is the key to a sustainable and low-carbon future. However, the efficient and safe storage and transport of hydrogen remains a huge challenge. Liquid organic hydrogen carriers (LOHCs) have recently attracted much attention as a breakthrough solution to this challenge, revolutionising the way we handle and use hydrogen. LOHC is a fuel-like organic aromatic molecule that is loaded with hydrogen in a catalytic hydrogenation reaction [1]. The perhydrogenated product can be stored for an unlimited time and transported to the customer, and the hydrogen can be released on demand. From a chemical point of view, the LOHC system, therefore, consists of a pair of organic molecules, the hydrogen-lean (HL) form of the system (e.g., benzene, naphthalene, etc.) and the hydrogen-rich (HR) form of the system (e.g., cyclohexane, decalin, etc.). An important aspect of all LOHC technologies is that the existing infrastructure of traditional fossil-fuel-based energy distribution technologies can generally be reused, which facilitates the establishment of the hydrogen economy.
Among the multitude of promising LOHCs, aromatic hydrocarbons are the most attractive from an economic point of view, although the dehydrogenation of these compounds is both kinetically and thermodynamically unfavourable. The high temperature of the methylcyclohexane/toluene system dehydrogenation process has led to a closer examination of other suitable polycyclic aromatic systems. Among the polycyclic compounds considered in the literature were decalin [2,3] bicyclohexyl, tercyclohexane [4], perhydroanthracene, perhydropyrene, perhydrofluorene, and perhydrofluoranthene [5], as well as hydrogenated forms of chrysene, picene, coronene and even hexabenzocoronene [6]. The gravimetric capacity, e.g., for the decalin/naphthalene system, is 7.2% wt. and 7.4% wt. for the perhydrocoronene/coronene system [7]. Of the compounds listed, decalin is the only one that is liquid under ambient conditions, while the other compounds are solids. As a rule, the aromatic counterparts of LOHC systems have high melting points, e.g., 148 °C for pyrene and 438 °C for coronene. To overcome this obstacle, it is necessary to use toluene to obtain a liquid, which reduces process safety. As an alternative, Pez et al. [6] suggested the use of compositions of two or more components that can form liquids or low-melting eutectic mixtures. The introduction of substituents, e.g., n-alkyl, alkoxyl or ether groups, into the ring structure of a polycyclic molecule also allows the melting point to be lowered, but the amount of hydrogen released decreases. According to Cooper et al. [8], compounds with high hydrogenation enthalpies generally require high dehydrogenation temperatures, and compounds with a large number of rings are characterised by relatively low hydrogenation enthalpies. Experimental data on the dehydrogenation of pentacene (73.2 kJ/mol H2) and coronene (57.7 kJ/mol H2) support this statement; however, for these compounds, it is shown that pentacene and coronene are able to release only 0.15 wt.% and 1 wt.%, respectively [9]. It has also been shown that π-conjugated aromatic molecules containing five-membered rings are more efficient hydrogen substrates because they have a lower enthalpy of dehydrogenation than the corresponding conjugated system in a six-membered ring. Decalin is an attractive carrier because it is liquid under normal conditions and has a high hydrogen capacity (7.2 wt%), with an energy density of 2.2 kWh/L. However, if toluene is used as a solvent, the capacity and energy density of the decalin/naphthalene system drop to 3.8 wt. % and 1.1 kWh/L, respectively [10]. In addition, the evaporation losses during storage are low due to the high boiling points, but the high temperature of the dehydrogenation reaction can lead to deactivation of the catalyst. A common disadvantage of the aromatic systems considered is the low selectivity in the dehydrogenation reaction phase caused by hydrocracking and hydrogenolysis reactions [7,10].
LOHC systems with oxygen functionality have been shown to offer very attractive hydrogen release performance at low temperatures and very attractive properties for potential technical applications [11]. For example, it has been shown that benzophenone can be selectively and completely hydrogenated to dicyclohexylmethanol over Ru/Al2O3 catalysts at 50 bar hydrogen pressure and temperatures in the range of 90 to 180 °C. Hydrogen can be released from dicyclohexylmethanol with Pt-based catalysts to recover benzophenon at a temperature of 250 °C. This reversible process is characterised by an excellent hydrogen storage capacity of 7.2 mass% [11]. This successful example has motivated further investigations of oxygen-containing LOHC systems. The focus of this work is the question of whether the introduction of oxygen functionality into aromatic molecules generally leads to an improvement in the thermodynamic properties of hydrogen storage compared to similarly shaped aromatic molecules without oxygen functionality.
It is known that the hydrogenation of aromatics is a thermodynamically favourable reaction. However, reversible dehydrogenation is thermodynamically unfavourable, and only the optimisation of the reaction conditions helps to shift the equilibrium towards a sufficient release of hydrogen. Understanding the general regularities between the structures of LOHC counterparts and their thermodynamic properties is essential for screening molecules suitable for practical applications. Typical dehydrogenation reactions of hydrocarbon-based LOHCs and their analogues with oxygen functionality are shown in Figure 1.
In this work, we carefully collected the available thermodynamic data for the HR and HL counterparts of the LOHC systems shown in Figure 1 and for related molecules. These data were evaluated by our own complementary measurements, structure–property correlations and high-level quantum chemical calculations and finally recommended for chemical engineering calculations.
The general feasibility of chemical reactions is determined by the sign and size of the standard Gibbs energy of the reaction,   r G m o , as follows (e.g., for a reaction in the ideal gas state):
r G m o = r H m o T × r S m o = R T × ln K p
with r H m o = the standard molar reaction enthalpy, r S m o is the change in the standard molar entropy of a chemical reaction, and K p = the gas-phase thermodynamic equilibrium constant. The Equation (1) is known as the Gibbs–Helmholtz equation, and it is broadly used in chemical engineering for the assessment of the K p -values, which are directly related to the yield of the desired products of the reaction. The negative and large r G m o values provide evidence of a large thermodynamic driving force for the chemical reactions of interest. The standard Gibbs energy of reaction,   r G m o , is not directly measurable; therefore, it is calculated from the measurable enthalpic and entropic contributions involved in Equation (1). The standard molar reaction enthalpies, r H m o , can be measured directly using reaction calorimetry or derived from temperature dependencies of equilibrium constants. However, most frequently the r H m o values are calculated from the standard molar enthalpies of formation, f H m o , of reaction participants. The change in the standard molar entropy of reactions, r S m o , are generally derived from absolute entropies of reaction participants.
Given the countless molecules suitable for hydrogen storage, it is hardly possible to obtain the required thermodynamic properties for the numerous hydrogen-rich and hydrogen-poor counterparts of LOHC systems. The only way to collect sufficient data for thermodynamic analysis and optimisation of hydrogen storage is a clever combination of experimental and theoretical methods, including quantum chemical (QC) methods and empirical correlations of structural properties. In particular, the development of modern QC methods enables significant accumulation thermodynamics for the screening of promising LOHC systems since molar enthalpies of gas-phase formation calculated with QC high-level methods, f H m o (g), nowadays reach a ‘chemical accuracy’ of ±(2–4) kJ·mol−1 for organic molecules with up to twenty heavy atoms in the structure [12].
The absolute gas-phase entropies, S m o (g), for molecules with moderate flexibility could also be calculated with a reasonable accuracy of ±(5–10) J·mol−1·K−1 [13]. However, it should be noted that these achievements refer to the thermodynamic data of a single molecule in the ideal gas state. Progress in the development of QC methods for the calculation of properties in the liquid or crystalline state is still insufficient for the modelling of processes related to hydrogen storage.
Admittedly, hydrogen storage with the LOHC systems based on aromatic compounds is conveniently carried out in the liquid phase under a hydrogen pressure of up to 30 bar [1]. For such experimental conditions, the corresponding reaction enthalpies,   r H m o , required for Equation (1) are calculated from the liquid-phase standard molar enthalpies of formation of reactants, f H m o (liq), using Hess’s Law. The liquid phase r S m o values are also calculated from the corresponding liquid-phase entropies S m o (liq) of the reaction participants. Consequently, the equilibrium constant of the liquid phase, K a , appears in Equation (1) instead of the thermodynamic equilibrium constant of the gas phase, K p .
A valuable way of overcoming the limitations of QC methods and making them usable for the liquid phase is given by Equation (2):
f H m o ( liq ,   298.15   K ) = f H m o ( g ,   298.15   K ) QC l g H m o ( 298.15   K )
where f H m o (g, 298.15 K)QC is the theoretical gas-phase enthalpy of formation calculated by any suitable high-level QC method, and l g H m o (298.15 K) is the standard molar enthalpy of vaporisation. There are numerous experimental methods to accurately measure l g H m o (298.15 K) values, and they can be relatively reliably estimated using various group additivity rules and structure–property correlations [14].
Similarly, the entropic contribution, S m o (liq, 298.15 K), to the Gibbs–Helmholtz equation (Equation (1)) can be derived by combining the QC and experimental methods as follows:
S m o ( liq ,   298.15   K ) = S m o ( g ,   298.15   K ) l g S m o ( 298.15   K )
where S m o (g, 298.15 K)QC is the theoretical gas-phase entropy calculated by any suitable high-level QC method and l g S m o (298.15 K) is the standard molar entropy of vaporisation, which can be derived from the slope of the temperature dependence of the experimental vapor pressures.
Such a combination of experimental and theoretical methods can greatly facilitate the screening of the promising LOHC and a reasonable workflow could be proposed to comprehensively determine the thermodynamic quantities essential for reversible hydrogen storage using Equation (1):
-
Step I: The typical representatives of the promising LOHC systems are studied using the conventional experimental methods (the high-precision combustion calorimetry for f H m o (liq/cr, 298.15 K) values and the transpiration method for l g H m o (298.15 K) and l g S m o (298.15 K) in this work)
-
Step II: The high-level QC calculations are performed to derive f H m o (g, 298.15 K)QC values and f S m o (g, 298.15 K)QC values and the possible experimental and empirical methods are involved to validate these theoretical results.
-
Step III: The liquid-phase f H m o (liq, 298.15 K) and r S m o (liq, 298.15 K) values are derived and used to calculate r G m o (liq, T). In addition, the equilibrium constants are derived according to Equation (1), and a thermodynamic analysis of the LOHC system is carried out.
These steps represent the general roadmap of the present work, which aims to evaluate the impact and practical potential of introducing oxygen functionality in LOHC systems based on aromatic ethers compared to similarly structured aromatics.

2. Theoretical and Experimental Methods

Provenance and purities of the commercial sample of of benzyl phenyl ether, dibenzyl ether and 2-methoxy-napthalene are shown in Table S1. A high-precision self-made calorimeter equipped with a static bomb [15] was used for measurements of standard molar combustion energies of benzyl phenyl ether and dibenzyl ether. A transpiration method [16] was used to measure the vapor pressures over the solid and liquid samples of benzyl phenyl ether, dibenzyl ether, and 2-methoxy-napthalene. The temperature dependence of the absolute vapour pressures for each compound was used to derive the enthalpies and entropies of vaporisation of the aromatic ethers.
The melting temperature and enthalpy of fusion of benzyl phenyl ether were measured by DSC [17]. The main details of the experimental techniques are contained in the Supplementary Materials.

Quantum-Chemical Calculations

The conformational analysis was performed with a computer code called CREST (conformer-rotamer ensemble sampling tool). The software package uses a metadynamics simulation-based screening procedure to generate molecular conformations in the gas phase. A history-dependent bias potential is applied, where the collective variables for the metadynamics are the previous structures on the potential energy surface, expressed as the root-mean-square deviation of the atoms between them, which are calculated according to the quaternion algorithm. Conformations are generated at the GFNn-xTB level in the iMTD-GC workflow. More details are provided in the Ref. [18]. The obtained structures were optimised with the B3LYP/6–31 g(d,p) method [19]. It was found that the energetics of the conformers were quite close to each other within a few kJ·mol−1. The most stable conformers for the aromatic compounds and their perhydrogenated products are shown in Figure S2.
The Gaussian 16 series software [20] was used for QC calculations. The total H298-enthalpies and the most stable conformers were calculated using the G3MP2 method [21]. The G3MP2 method was chosen as the optimal approach, combining high accuracy of energy calculation with reasonable computational cost. This method has been systematically used in our laboratory for calculations of C-, H-, N- and O-containing molecules, and in our experience, G3MP2 provides the theoretical enthalpies of formation in very good agreement with the reliable experimental values. The H298 values were finally converted to the theoretical f H m o (g, 298.15 K)QC values and discussed. The well-established assumption “rigid rotator-harmonic oscillator” was used for the quantum-chemical calculations. The S m o (g, 298.15 K)QC values were calculated according to (Equation 1) using H298 and G298 from the output file. The general details of the calculation procedure have been published elsewhere [22].

3. Results and Discussion

3.1. Step I: Evaluation of Experimental and Empirical Thermochemical Results

3.1.1. Enthalpies of Formation from Combustion Calorimetry

The standard specific energies of combustion c u o (cr) of benzyl phenyl ether and c u o (liq) of dibenzyl ether were determined from five experiments for each compound. These c u o (cr/liq) values were used to calculate the experimental standard molar enthalpies of combustion, c H m o , which refer to the reactions:
benzyl phenyl ether  C13H12O (cr) + 15.5 O2 (g) = 13 CO2 (g) + 6 H2O (liq)
dibenzyl ether   C14H14O (liq) + 17 O2 (g) = 14 CO2 (g) + 7 H2O (liq)
The experimental standard molar enthalpies of combustion, c H m o , were used to derive the standard molar enthalpies of formation in the crystal or liquid state f H m o (cr/liq) for both compounds. The results of combustion experiments are given in Table 1.
The standard molar enthalpies of formation, f H m o (cr/liq), (see Table 1) were calculated based on the c H m o values of the reaction according to Equations (4) and (5) using Hess’s Law and the standard molar enthalpies of formation of H2O(liq) and CO2(g) assigned by CODATA [24]. The total uncertainties of c H m o and f H m o values have been calculated according to the guidelines presented by Hubbard et al. [25] and Olofsson [23]. The uncertainties of combustion energies, c u o (cr/liq), are expressed as the standard deviation of the mean. According to the thermochemical practice, the uncertainties assigned to the f H m o (cr/liq) values are twice the overall standard deviations and include the uncertainties of the calibration, the combustion energies of the auxiliary materials, and the uncertainties of the enthalpies of formation of the reaction products H2O and CO2. The combustion experiments with benzyl phenyl ether and dibenzyl ether were performed for the first time.

3.1.2. Absolute Vapour Pressures and Vaporisation Enthalpies

The primary vapour pressure–temperature dependencies for benzyl phenyl ether, dibenzyl ether and 2-methoxy-naphthalene measured in this work (see Table S3), and those from the literature (see Table S4) were fitted uniformly using a three-parametric equation (see ESI for details). From these data, the standard molar enthalpies of vaporisation at the respective temperatures T were derived. The final results at the reference temperature T = 298.15 K, l g H m o (298.15 K), are summarised and compared in Table 2.
Only a few vapour pressure temperature dependencies of aromatic and aliphatic ethers relevant to this work can be found in the literature [26,27,28]. Furthermore, the enthalpies of vaporisation were never derived from these results. The available data for these compounds were collected and treated consistently in this work (see details in ESI) to calculate the l g H m o (298.15 K) values in the same way as our own results from the transpiration method.
The results show that no systematic vapour pressure measurements were available for most of the hydrogenated aromatic ethers. However, the experimental boiling points (BP) of the desired aliphatic cyclic ethers at different reduced pressures can be found in the literature (see a compilation in Table S4) and may be helpful in obtaining the missing thermodynamic information. These boiling points usually originate from the distillation of reaction mixtures after synthesis and are generally only used to identify the compounds. In such experiments, the pressures are usually measured with uncalibrated manometers, and the temperatures are given in the range of a few degrees. However, our previous work (e.g., on methylbiphenyls [33]) has shown that, in general, reasonable trends can be derived even from such rough data. This general conclusion also applies to the aromatic ethers. As shown in Table 2, the l g H m o (298.15 K) values for benzyl phenyl ether and dibenzyl ether derived from BP at various reduced pressures agree with the enthalpies of vaporisation measured by conventional methods within the experimental uncertainties.
Therefore, we carefully collected and evaluated the available thermodynamic data and BP for all reactions R1 to R10 in Figure 1 from the literature and used these data to calculate their l g H m o (298.15 K) and l g S m o (298.15 K) values (see Tables S3 and S4), which are required for the thermodynamic analysis of the dehydrogenation reactions shown in Figure 1.
As shown in Figure 1, the thermodynamics of the dehydrogenation of aromatic ethers are compared with those of similarly structured aromatic hydrocarbons and their fully dehydrogenated aliphatic counterparts. While reliable experimental data from the literature are available for the pairs benzene and cyclohexane, diphenyl–methane and dicyclohexyl–methane, naphthalene, and decalin (see Table 2), the thermodynamic data for the pairs 1,2-diphenyl-ethane and 1,2-dicyclohexyl-ethane, 1,3-diphenyl-propane, and 1,3-dicyclohexyl-propane were evaluated in this work (see Table 2) and these data were recommended for thermochemical calculations.

3.1.3. Validation of Vaporisation Enthalpies by Structure–Property Correlations

Structure–property correlations are the backbone of physical–organic chemistry, as the new experimental results are integrated into the network of already known and evaluated data. Provided that the new results are consistent with the trends already known from the available results, they can be recommended for the thermochemical calculations. If this is not the case, the unusual structural effects could be anticipated, or the measurements should be repeated to avoid confounding. In this context, before using l g H m o (298.15 K) values derived from single experiments using conventional methods or from BP, it is useful to validate the results by structure–property correlations, e.g., with relationships between the enthalpy of vaporisation and normal boiling points (Tb). For example, in our recent work [34], a robust l g H m o (298.15 K)−Tb relationship was established for a large series of alkyl-cyclohexanes:
l g H m o ( 298.15   K ) = 0.1972 × T b 37.6   with   R 2 = 0.989
We used this correlation to assess the l g H m o (298.15 K)-values for 1,2-dicyclohexyl-ethane and 1,3-dicyclohexyl-propane, and the results (see Table 2, indicated as Tb) agreed very well with the vaporisation enthalpies derived from the approximation of BP selected in Table S4.
Group additivity (GA) methods are another valuable way of correlating structural properties. These methods are generally based on the idea of the Lego box, where a property of a molecule is assembled from small building blocks with well-established numerical contributions. Comprehensive systems of group contributions (or increments) have been developed for most major classes of organic compounds [14]. The set of group contributions used in this work for the prediction of thermodynamic properties is listed in Table S8. However, it is well known that the conventional GA methods have difficulties with cyclic molecules. To improve the quality of the prediction, numerous and unique ‘ring’ corrections are proposed [14,35]. However, for the molecules investigated in this work, it might be more practical to use a phenyl fragment (C6H5) and a cyclohexyl fragment (C6H11-) developed as shown in Figure S1 instead of the ‘ring’ corrections. These fragments were used to calculate the thermodynamic properties of aromatic and aliphatic ethers and hydrocarbons and their hydrogenated counterparts. The examples of GA calculations for the LOHC pairs dibenzyl ether/dicyclohexylmethyl ether and 1,3-diphenylpropane/1,3-dicyclohexylpropane are shown in Figure 2.
The results of the GA calculations are summarised in Table 2. For all aromatic and aliphatic compounds, the enthalpies of vaporisation l g H m o (298.15 K) evaluated in this way are in good agreement with results obtained using other methods. To gain more confidence in the enthalpies of vaporisation, the weighted average values were calculated for each species considered in this work (see Table 2) and used to calculate the enthalpies of formation of the liquid phase in Section 3.3

3.2. Step II: Quantum-Chemical Calculations of the Gas-Phase Enthalpies of Formation

The aromatic and aliphatic ethers are flexible molecules and are represented in the gas phase by a mixture of conformers. The energies E0 and the total enthalpies H298 of the most stable conformers were finally calculated using the method G3MP2 implemented in the software Gaussian 16. The resulting H298 enthalpies were converted to the standard molar enthalpies of formation based on conventional atomisation (AT) procedure [22] (e.g., for benzyl phenyl ether):
C13H12O → 13 × C + 12 × H + 1 × O
The results of the quantum chemical G3MP2 calculations for the hydrogen-lean and hydrogen-rich counterparts of the LOHC systems are summarised in Table 3 (column 3).
Admittedly, the AT reaction to convert the resulting H298 enthalpies to the standard molar enthalpies of formation is the least biased choice compared to other isodesmic, isogyric, homodesmotic, etc., reactions [39]. The main advantage of the AT reaction is that the enthalpies of the atoms are precisely known, compared to the less precise enthalpies of the formation of the various molecules involved in the construction of the aforementioned reactions. However, it is also known [36,39] that the QC results calculated with the G3MP2 and the AT reactions,   f H m o (g, AT), show a small but noticeable deviation from the experimental results (see Table 3). Nevertheless, a simple correction derived from a set of similarly structured molecules with reliable experimental enthalpies of formation, f H m o (g)exp, (see Table 3, Tables S10 and S11) can reconcile QC and experimental results if the f H m o (g, AT) values are correlated with the experimental values. The following equation,
f H m o ( g ) QC / kJ · mol 1 = 1.0194   ×   f H m o ( g , AT ) G 3 MP 2 + 12.3    with R 2 = 0.9996
was used to obtain the “corrected” QC results for the aromatic ethers and their hydrogenated products. The correction for 1,2-diphenyl ethane and 1,3-diphenyl-propane was made with the following equation:
f H m o ( g ) QC / kJ · mol 1 = 0.9576   ×   f H m o ( g , AT ) G 3 MP 2 + 18.5    with R 2 = 0.9979
which was derived from the data for similarly structured alkyl-biphenyls (see Table S10). For their hydrogenated products, the following equation was applied:
f H m o ( g ) QC / kJ · mol 1 = 0.970   ×   f H m o ( g , AT ) G 3 MP 2 3.2    with R 2 = 0.9998
derived from the data for similarly structured alkyl-cyclohexanes (see Table S11). It is obvious that the coefficients of Equations (8) to (10) are slightly different because the elementary compositions of the molecules in the sets are different, as well as the structures of the molecules included in the sets (e.g., aromatic molecules in Table S10 and aliphatic molecules in Table S11). The enthalpies of formation calculated with the corrected atomisation reactions are given in Table 3, column 5, and used for further thermochemical calculations (see Section 3.3).
The QC results for benzyl phenyl ether, dibenzyl ether and 2-methoxy-naphthalene were validated with the experimental values evaluated in Table 4.

3.3. Step III: Thermodynamic Analysis of the LOHC Systems

3.3.1. Comparison of the Energetics of Dehydrogenation Reactions with and without Oxygen Functionality

The standard molar enthalpies of chemical reactions were calculated from the standard molar enthalpies of formation, f H m o , of the reactants according to Hess´s Law:
r H m o ( T ) = [ f H m o ( T ,   products ) ] [ f H m o ( T ,   educts ) ]
Unfortunately, f H m o values for some cyclic aliphatic molecules of the hydrogen-proximal LOHC counterparts of the R1–R10 reactants are missing in the literature. However, they can be reliably calculated at a high level using QC methods, as shown in our recent work for the series of alkyl-cyclohexanes [34]. Therefore, we calculated the missing gas-phase formation enthalpies of aromatic ethers, aromatic hydrocarbons and their perhydrogenated counterparts with G3MP2 (see Table 3, column 5) and converted them into the liquid-phase formation enthalpies by subtracting the enthalpies of vaporisation as shown in Equation (2). The required l g H m o (298.15 K) contributions in Equation (2) for aromatic and hydrogenated LOHC counterparts were already evaluated in Table 2. The f H m o (liq, 298.15 K) values, derived according to Equation (2), are given in Table 5 and can now be used as input to Hess´s Law (Equation (11)).
The results of the calculation of the liquid-phase reaction enthalpies, r H m o (liq), of the dehydrogenation reactions R1 to R10 are shown in Table 6.
As can be seen from this table, the reactions R1 to R10 are strongly endothermic from 196.8 to 409.8 kJ·mol−1. However, it should be noted that different numbers of hydrogen molecules are involved in these reactions. Therefore, the enthalpy of the reaction is usually related to the amount of hydrogen released as a kJ·mol−1/H2 unit in order to be able to correctly compare the energetics of hydrogenation reactions, taking the stoichiometry into account. As can be seen from Table 6, the dehydrogenation enthalpies in these units are generally at the level of ≈65 kJ·mol−1/H2 for all pairs. In general, the lower the reaction enthalpy, the better and more effective the dehydrogenation process is from a thermodynamic point of view. According to the results shown in Table 6, the influence of oxygen functionality is favourable for the cyclohexane/benzene, methoxy–cyclohexane/methoxy–benzene, and 2-methoxy-decalin/2-methoxy-naphthalene pairs. For other LOHC pairs, the influence is quite small, and even the dehydrogenation of oxygenated cycloalkanes is somewhat less thermodynamically favoured compared to the aliphatic hydrocarbons.

3.3.2. Comparison of the Reaction Entropies of Dehydrogenation of Lohc Based on Ethers and Hydrocarbons

The liquid-phase standard molar entropies, S m o (liq, 298.15 K), were calculated according to Equation (3) using the gas-phase standard molar entropies S m o (g, 298.15 K) from G3MP2 calculations and the standard molar entropy of vaporisation, l g S m o (298.15 K), derived (details are given in ESI) from the intercept of the temperature dependencies of the experimental vapor pressures. The results for l g S m o (298.15 K), S m o (g, 298.15 K), and final entropies, S m o (liq, 298.15 K), are shown in Table 7.
r S m o ( T ) = [ S m o ( T ,   products ) ] [ S m o ( T ,   educts ) ]
The resulting total liquid-phase reaction entropies, r S m o (liq, 298.15 K), of the dehydrogenation reactions R1 to R10 are given in Table 8, column 3.
As can be seen from Table 8, the reaction entropies of R1 to R10 are between 357 and 757 J·mol−1·K−1. For an appropriate comparison, however, they should also be related to the amount of hydrogen released, whereby the stoichiometry must be taken into account. As shown in Table 8, the dehydrogenation entropies in all LOHC pairs are barely distinguishable at the level of ≈120 J·mol−1·K−1/H2, and this observation is important for chemical engineering calculations with the LOHC systems. As also shown in Table 8, the entropy contributions ( T × r S m o ) to the Gibbs energy of the reactions R1 to R10 very close for all reactions on the level of ≈36 J·mol−1·K−1. From a thermodynamic point of view, the equilibrium position for these reactions is, therefore, mainly determined by the enthalpic term. Thus, both the required enthalpic ( f H m o (liq, 298.15 K) and entropic r S m o (liq, 298.15 K) contributions to the Gibbs–Helmholtz equation (Equation (1)) have now been developed and are ready for the comparison of the Gibbs energies of the dehydrogenation reactions R1 to R10 in the following section.

3.3.3. Comparison of the Gibbs Energies of Dehydrogenation Reactions R1 to R10

The total liquid-phase Gibbs energies and equilibrium constants of dehydrogenation reactions R1 to R10 calculated according to Equation (1) are given in Table 9.
As shown in Table 9, the r G m o values for all pairs of LOHC molecules are positive and strongly temperature-dependent. The oxygen functionality significantly decreases the r G m o values for pairs methoxy–cyclohexane/methoxy–benzene and 2-methoxy-decalin/2-methoxy-naphthalene. For the system dicyclohexyl/diphenyl ether, the introduction of oxygen instead of the methylene group does not change the energetics practically. For systems cyclohexylmethoxy–cyclohexane/benzyl phenyl ether and dicyclohexylmethyl ether/dibenzylether, the oxygen functionality even increases the r G m o values (compared to the similarly shaped hydrocarbons), making the thermodynamics of these LOHC systems even less favourable.
The influence of the oxygen functionality becomes clearer when comparing the equilibrium constants Ka (e.g., at 500 K), as the size of the equilibrium constant is directly related to the amount of the compound with a low hydrogen content in the reaction mixture of the dehydrogenation reaction. For example, for the methoxy–cyclohexane/methoxy–benzene pair, the constant Ka = 3.0 × 10−2 (see Table 9, last column) is approximately an order of magnitude greater than the constant Ka = 6.2 × 10−3 (see Table 9, last column) for the cyclohexane/benzene pair.
For the 2-methoxy-decalin/2-methoxy-naphthalene pair, the constant Ka = 1.4 × 10−1 (see Table 9, last column) is also an order of magnitude greater than the constant Ka = 2.3 × 10−2 (see Table 9, last column) for the decalin/naphthalene pair. For the dicyclohexyl/diphenyl ether system, the constant Ka = 4.5 × 10−3 (see Table 9, last column) is about the same as the constant Ka = 4.9 × 10−3 (see Table 9, last column) for the dicyclohexyl–methane/diphenyl–methane system and the introduction of oxygen instead of the methylene group practically does not change the energetics. In the cyclohexylmethoxy–cyclohexane/benzylphenyl ether and dicyclohexylmethyl ether/dibenzyl ether systems, the oxygen functionality reduces the Ka values compared to the corresponding hydrocarbon-based LOHC systems.

3.3.4. Thermodynamic Analysis of the Reversible Hydrogenation/Dehydrogenation Process in the Gas Phase

Admittedly, the hydrogen release from the hydrogen-rich LOHC counterparts is a thermodynamically unfavourable process [1]. However, a theoretical optimal dehydrogenation temperature can be evaluated based on quantum chemical calculations of the Gibbs free energy according to Equation (1). If the system approximates the equilibrium, the free reaction energy approaches zero. A theoretically optimal dehydrogenation temperature can be determined on the basis of quantum chemical calculations of the Gibbs free energy according to Equation (1). When the system approaches equilibrium, the free energy of the reaction approaches zero. When the r G m o = 0, the enthalpy term r H m o and the entropy term T × r S m o in Equation (1) become equal. This occurs at a certain “equilibrium temperature”, T e q , which can imply the reversal of thermodynamic feasibility from dehydrogenation to hydrogenation tendency:
T e q = r H m o r S m o
This temperature can be considered a valuable indicator for selecting the optimum temperature for a particular LOHC system. This T e q for the dehydrogenation step should not be too high, but sufficient to achieve good conversion and selectivity. For the hydrogenation step, the T e q should not be too low in order to achieve acceptable reaction rates. The results of the T e q calculations according to Equation (13) for the reactions R1 to R10 are compared in Table 10.
As can be seen from Table 10, the T e q values for dehydrogenation in LOHC systems based on aliphatic cyclic hydrocarbons vary between 530 and 569 K. The T e q values for dehydrogenation in aliphatic cyclic ethers are generally somewhat lower and vary between 525 and 554 K. The lowest equilibrium temperature is observed for trans-2-methoxydecaline (525 K). However, the positive influence of the oxygen function is only significant for methoxy–cyclohexane and trans-2-methoxy-decaline. For other pairs, the T e q values are slightly higher for oxygen-containing compounds (see Table 10).

4. Conclusions

The thermodynamic properties of a series of aromatic ethers have been studied using different experimental and theoretical methods in order to evaluate their applicability as LOHC materials. The consistent sets of phase transition enthalpies and entropies, enthalpies of formation, and absolute entropies were obtained using vapour pressure measurements, combustion calorimetry, and empirical and quantum chemical calculations. These results were recommended for chemical engineering calculations to optimise the storage and release of hydrogen using the reversible hydrogenation/dehydrogenation reactions. The thermodynamic properties of the liquid-phase hydrogenation/dehydrogenation reactions in LOHC systems based on methoxy–benzene, diphenyl ether, benzyl phenyl ether, dibenzyl ether, and 2-methoxy-naphthalene were derived and compared with the data for the similarly structured hydrogen carriers based on benzene, diphenyl-methane, 1,2-diphenyl-ethane, 1,3-diphenyl-propane and naphthalene. The effects of oxygen functionality on the thermodynamic properties of hydrogenation/dehydrogenation reactions were evaluated. It was found that for methoxy–benzene and 2-methoxy-naphthalene, a reduction in enthalpies and Gibbs-free energies for the dehydrogenation reactions can be observed, allowing more favourable reaction conditions to be achieved. For systems based on diphenyl ether, the thermodynamic properties of the reaction are very similar to those of the diphenyl–methane system. In the systems based on benzyl phenyl ether and dibenzyl ether, the introduction of oxygen functionality is not advantageous compared to LOHC systems based on 1,2-diphenylethane and 1,3-diphenylpropane. These new findings contribute to a better understanding of the reaction thermodynamic properties of LOHC systems and reveal the structural features of aromatic and aliphatic cyclic compounds that are essential for their technical applications. Further investigation of the thermodynamic properties of hydrogenation/dehydrogenation reactions for LOHC systems containing one or two oxygen functions in the aliphatic ring should complete this interesting structure–property evaluation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/oxygen4030015/s1, Figure S1: The development of the contributions C6H5- (left) and C6H11- (right) for calculating the vaporization enthalpies (or gas-phase enthalpies of formation) of aromatic and cyclic aliphatic hydrocarbons using the GA method. All values in kJ·mol−1. The numerical values of the contributions are given in Table S8. Results of calculations are shown in Table S9; Figure S2: The most stable conformers for the aromatic compounds and their perhydrogenated products; Table S1: Provenance and purity of the materials; Table S2: Auxiliary quantities: formula, density ρ(293 K), massic heat capacity cp(298.15 K), and expansion coefficients (δVT)p of the materials used in the present study; Table S3: Absolute vapour pressures p, and standard molar vaporisation enthalpies and entropies determined using the transpiration method; Table S4: The vapour pressures p, and standard vaporisation enthalpies and entropies obtained by the approximation of boiling points at different pressures available in the literature; Table S5: Enthalpies of fusion ( cr l H m o ), vaporization ( l g H m o ) and sublimation ( cr g H m o ) at melting points (Tfus) and at 298.15 K of benzyl phenyl ether and 2-methoxy-naphthalene; Table S6: Compilation of data on molar heat capacities C p , m o (liq) and heat capacity differences l g C p , m o (in J·K−1·mol−1) at T = 298.15 K; Table S7: Coefficients of the Clarke and Glew equation (S6) for 1,2-diphenyl-ethane from Osborn et al; Table S8: Group-additivity values Γi for calculation of enthalpies of vaporization, l g H m o , and enthalpies of formation, f H m o (g), of alkanes, ethers, and aromatics at 298.15 K in kJ mol−1; Table S9: Comparison of experimental thermochemical data with those calculated by group additivity at T = 298.15 K ( = 0.1 MPa) for ethers and hydrocarbons (in kJ·mol−1); Table S10: Calculation of the G3MP2 corrected gas-phase enthalpies at T = 298.15 K (in kJ·mol−1); Table S11: Comparison of experimental and G3MP2 calculated gas-phase enthalpies of formation for alkyl-cyclohexanes (at T = 298.15 K, = 0.1 MPa, in kJ·mol−1); Table S12: Thermochemical data at T = 298.15 K ( = 0.1 MPa) for auxiluary compounds (in kJ·mol−1); Table S13: Liquid phase thermodynamic properties of the hydrogenation of the hydrogen-lean LOHC counterparts.

Author Contributions

Conceptualization, S.P.V.; methodology, A.A.S., S.V.V. and S.P.V.; software, A.A.S.; validation, A.A.S., S.V.V. and S.P.V.; formal analysis, A.A.S. and S.P.V.; investigation, S.V.V.; resources, S.V.V.; data curation, S.V.V. and S.P.V.; writing—original draft preparation, S.P.V. and A.A.S.; writing—review and editing, S.P.V. and A.A.S.; visualization, A.A.S. and S.V.V.; supervision, S.P.V.; project administration, S.P.V.; funding acquisition, S.P.V. All authors have read and agreed to the published version of the manuscript.

Funding

SPV gratefully acknowledges the financial support of the German Research Foundation within the framework of SPP 1807 “Control of London Dispersion Interactions in Molecular Chemistry”, grant VE 265-9/2. This paper has been supported by the Kazan Federal University Strategic Academic Leadership Program (“PRIORITY-2030”). AAS gratefully acknowledges the Committee on Science and Higher Education of the Government of St. Petersburg. The work was supported by the Ministry of Science and Higher Education of the Russian Federation (theme No. FSSE-2024-0021) as part of the state task of the Samara State Technical University (creation of new youth laboratories).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Preuster, P.; Papp, C.; Wasserscheid, P. Liquid Organic Hydrogen Carriers (LOHCs): Toward a Hydrogen-Free Hydrogen Economy. Acc. Chem. Res. 2017, 50, 74–85. [Google Scholar] [CrossRef] [PubMed]
  2. Hodoshima, S.; Takaiwa, S.; Shono, A.; Satoh, K.; Saito, Y. Hydrogen Storage by Decalin/Naphthalene Pair and Hydrogen Supply to Fuel Cells by Use of Superheated Liquid-Film-Type Catalysis. Appl. Catal. A Gen. 2005, 283, 235–242. [Google Scholar] [CrossRef]
  3. Sebastián, D.; Bordejé, E.G.; Calvillo, L.; Lázaro, M.J.; Moliner, R. Hydrogen Storage by Decalin Dehydrogenation/Naphthalene Hydrogenation Pair over Platinum Catalysts Supported on Activated Carbon. Int. J. Hydrogen Energy 2008, 33, 1329–1334. [Google Scholar] [CrossRef]
  4. Sung, J.; Choo, K.; Kim, T.; Tarasov, A.; Tkachenko, O.; Kustov, L. A New Hydrogen Storage System Based on Efficient Reversible Catalytic Hydrogenation/Dehydrogenation of Terphenyl. Int. J. Hydrogen Energy 2008, 33, 2721–2728. [Google Scholar] [CrossRef]
  5. Kustov, L.M.; Kalenchuk, A.N.; Bogdan, V.I. Systems for Accumulation, Storage and Release of Hydrogen. Russ. Chem. Rev. 2020, 89, 897–916. [Google Scholar] [CrossRef]
  6. Pez, G.P.; Scott, A.R.; Cooper, A.C.; Cheng, H.; Wilchelm, F.C.; Abdourazak, A.H. Hydrogen storage by reversible hydrogenation of pi-conjugated substrates. US patent US7351395B1, 2008. Available online: https://patents.google.com/patent/US7351395B1/en (accessed on 26 June 2024).
  7. Bourane, A.; Elanany, M.; Pham, T.V.; Katikaneni, S.P. An Overview of Organic Liquid Phase Hydrogen Carriers. Int. J. Hydrogen Energy 2016, 41, 23075–23091. [Google Scholar] [CrossRef]
  8. Cooper, A.C.; Campbell, K.M.; Pez, G.P. An Integrated Hydrogen Storage and Delivery Approach Using Organic Liquid-Phase Carriers. In Proceedings of the 16th World Hydrogen Energy Conference 2006, Lyon, France, 13–16 June 2006; Volume 3, pp. 2164–2175. [Google Scholar]
  9. Products, A.; Cooper, A.C.; Manager, P. Design and Development of New Carbon-Based Sorbent Systems for an Effective Containment of Hydrogen Final Report Revised April 2012 Hydrogen Storage by Reversible Hydrogenation of Solid- and Liquid-Phase Hydrogen Carriers; Air Products and Chemicals, Inc.: Allentown, PA, USA, 2012. [Google Scholar]
  10. Crabtree, R.H. Hydrogen Storage in Liquid Organic Heterocycles. Energy Environ. Sci. 2008, 1, 134. [Google Scholar] [CrossRef]
  11. Zakgeym, D.; Hofmann, J.D.; Maurer, L.A.; Auer, F.; Müller, K.; Wolf, M.; Wasserscheid, P. Better through Oxygen Functionality? The Benzophenone/Dicyclohexylmethanol LOHC-System. Sustain. Energy Fuels 2023, 7, 1213–1222. [Google Scholar] [CrossRef]
  12. Emel’yanenko, V.N.; Zaitseva, K.V.; Agapito, F.; Martinho Simões, J.A.; Verevkin, S.P. Benchmark Thermodynamic Properties of Methylanisoles: Experimental and Theoretical Study. J. Chem. Thermodyn. 2015, 85, 155–162. [Google Scholar] [CrossRef]
  13. Verevkin, S.P.; Zaitsau, D.H.; Emel’yanenko, V.N.; Zhabina, A.A. Thermodynamic Properties of Glycerol: Experimental and Theoretical Study. Fluid Phase Equilib. 2015, 397, 87–94. [Google Scholar] [CrossRef]
  14. Verevkin, S.P.; Emel’yanenko, V.N.; Diky, V.; Muzny, C.D.; Chirico, R.D.; Frenkel, M. New Group-Contribution Approach to Thermochemical Properties of Organic Compounds: Hydrocarbons and Oxygen-Containing Compounds. J. Phys. Chem. Ref. Data 2013, 42, 033102. [Google Scholar] [CrossRef]
  15. Emel’yanenko, V.N.; Verevkin, S.P.; Heintz, A. The Gaseous Enthalpy of Formation of the Ionic Liquid 1-Butyl-3-Methylimidazolium Dicyanamide from Combustion Calorimetry, Vapor Pressure Measurements, and Ab Initio Calculations. J. Am. Chem. Soc. 2007, 129, 3930–3937. [Google Scholar] [CrossRef] [PubMed]
  16. Verevkin, S.P.; Emel’yanenko, V.N. Transpiration Method: Vapor Pressures and Enthalpies of Vaporization of Some Low-Boiling Esters. Fluid Phase Equilib. 2008, 266, 64–75. [Google Scholar] [CrossRef]
  17. Emel’yanenko, V.N.; Zaitsau, D.H.; Shoifet, E.; Meurer, F.; Verevkin, S.P.; Schick, C.; Held, C. Benchmark Thermochemistry for Biologically Relevant Adenine and Cytosine. A Combined Experimental and Theoretical Study. J. Phys. Chem. A 2015, 119, 9680–9691. [Google Scholar] [CrossRef] [PubMed]
  18. Pracht, P.; Bohle, F.; Grimme, S. Automated Exploration of the Low-Energy Chemical Space with Fast Quantum Chemical Methods. Phys. Chem. Chem. Phys. 2020, 22, 7169–7192. [Google Scholar] [CrossRef] [PubMed]
  19. Petersson, G.A.; Bennett, A.; Tensfeldt, T.G.; Al-Laham, M.A.; Shirley, W.A.; Mantzaris, J. A Complete Basis Set Model Chemistry. I. The Total Energies of Closed-shell Atoms and Hydrides of the First-row Elements. J. Chem. Phys. 1988, 89, 2193–2218. [Google Scholar] [CrossRef]
  20. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  21. Curtiss, L.A.; Raghavachari, K.; Redfern, P.C.; Pople, J.A. Gaussian-3 Theory Using Scaled Energies. J. Chem. Phys. 2000, 112, 1125–1132. [Google Scholar] [CrossRef]
  22. Verevkin, S.P.; Emel’yanenko, V.N.; Notario, R.; Roux, M.V.; Chickos, J.S.; Liebman, J.F. Rediscovering the Wheel. Thermochemical Analysis of Energetics of the Aromatic Diazines. J. Phys. Chem. Lett. 2012, 3, 3454–3459. [Google Scholar] [CrossRef] [PubMed]
  23. Olofsson, G. Assignment of Uncertainties. In Combustion Calorimetry: Experimental Chemical Thermodynamics; Sunner, S., Månsson, M., Eds.; Pergamon: New York, NY, USA, 1979; pp. 137–161. [Google Scholar]
  24. Cox, J.D.; Wagman, D.D.; Medvedev, V.A. CODATA Key Values for Thermodynamics: Final Report of the CODATA Task Group on Key Values for Thermodynamics; CODATA Series on Thermodynamic; Hemisphere Publishing Corp.: New York, NY, USA, 1989. [Google Scholar]
  25. Hubbard, W.N.; Scott, D.W.; Waddington, G. Standard States and Corrections for Combustions in a Bomb at Constant Volume. In Experimental Thermochemistry; Rossini, F.D., Ed.; Interscience Publishers: New York, NY, USA, 1956; pp. 75–128. [Google Scholar]
  26. Stull, D.R. Vapor Pressure of Pure Substances. Organic and Inorganic Compounds. Ind. Eng. Chem. 1947, 39, 517–540. [Google Scholar] [CrossRef]
  27. Stephenson, R.M.; Malanowski, S. Handbook of the Thermodynamics of Organic Compounds; Springer: Dordrecht, The Netherlands, 1987; ISBN 978-94-010-7923-5. [Google Scholar]
  28. Barontini, F.; Cozzani, V. Thermogravimetry as a Screening Tool for the Estimation of the Vapor Pressures of Pure Compounds. J. Therm. Anal. Calorim. 2007, 89, 309–314. [Google Scholar] [CrossRef]
  29. Emel’yanenko, V.N.; Zaitsau, D.H.; Pimerzin, A.A.; Verevkin, S.P. Benchmark Properties of Diphenyl Oxide as a Potential Liquid Organic Hydrogen Carrier: Evaluation of Thermochemical Data with Complementary Experimental and Computational Methods. J. Chem. Thermodyn. 2018, 125, 149–158. [Google Scholar] [CrossRef]
  30. Osborn, A.G.; Scott, D.W. Vapor Pressures of 17 Miscellaneous Organic Compounds. J. Chem. Thermodyn. 1980, 12, 429–438. [Google Scholar] [CrossRef]
  31. Bolmatenkov, D.N.; Yagofarov, M.I.; Notfullin, A.A.; Solomonov, B.N. Calculation of the Vaporization Enthalpies of Alkylaromatic Hydrocarbons as a Function of Temperature from Their Molecular Structure. Fluid Phase Equilib. 2022, 554, 113303. [Google Scholar] [CrossRef]
  32. Verevkin, S.P.; Samarov, A.A.; Vostrikov, S. V Thermodynamics of Reversible Hydrogen Storage: Does Alkoxy-Substitution of Naphthalene Yield Functional Advantages for LOHC Systems? J. Chem. Phys. 2024, 160, 154709. [Google Scholar] [CrossRef] [PubMed]
  33. Samarov, A.A.; Verevkin, S.P. Hydrogen Storage Technologies: Methyl-Substituted Biphenyls as an Auspicious Alternative to Conventional Liquid Organic Hydrogen Carriers (LOHC). J. Chem. Thermodyn. 2022, 165, 106648. [Google Scholar] [CrossRef]
  34. Verevkin, S.P.; Samarov, A.A.; Vostrikov, S.V.; Wasserscheid, P.; Müller, K. Comprehensive Thermodynamic Study of Alkyl-Cyclohexanes as Liquid Organic Hydrogen Carriers Motifs. Hydrogen 2023, 4, 42–59. [Google Scholar] [CrossRef]
  35. Benson, S.W. Thermochemical Kinetics, 2nd ed.; John Wiley & Sons: New York, NY, USA; London, UK; Sydney, Australia; Toronto, ON, Canada, 1976; pp. 1–320. [Google Scholar]
  36. Verevkin, S.P.; Emel’yanenko, V.N.; Pimerzin, A.A.; Vishnevskaya, E.E. Thermodynamic Analysis of Strain in the Five-Membered Oxygen and Nitrogen Heterocyclic Compounds. J. Phys. Chem. A 2011, 115, 1992–2004. [Google Scholar] [CrossRef]
  37. Steele, W.V.; Chirico, R.D.; Knipmeyer, S.E.; Nguyen, A.; Smith, N.K.; Tasker, I.R. Thermodynamic Properties and Ideal-Gas Enthalpies of Formation for Cyclohexene, Phthalan (2,5-Dihydrobenzo-3,4-Furan), Isoxazole, Octylamine, Dioctylamine, Trioctylamine, Phenyl Isocyanate, and 1,4,5,6-Tetrahydropyrimidine. J. Chem. Eng. Data 1996, 41, 1269–1284. [Google Scholar] [CrossRef]
  38. Chirico, R.; Archer, D.; Hossenlopp, I.; Nguyen, A.; Steele, W.; Gammon, B. The Thermodynamic Properties of Chroman and Isochroman. J. Chem. Thermodyn. 1990, 22, 665–682. [Google Scholar] [CrossRef]
  39. Wheeler, S.E.; Houk, K.N.; Schleyer, P.v.R.; Allen, W.D. A Hierarchy of Homodesmotic Reactions for Thermochemistry. J. Am. Chem. Soc. 2009, 131, 2547–2560. [Google Scholar] [CrossRef]
  40. Verevkin, S.P.; Vostrikov, S.V.; Leinweber, A.; Wasserscheid, P.; Müller, K. Thermochemical Properties of Benzyltoluenes and Their Hydro- and Perhydro-Derivatives as Key Components of a Liquid Organic Hydrogen Carrier System. Fuel 2023, 335, 126618. [Google Scholar] [CrossRef]
  41. Roux, M.V.; Temprado, M.; Chickos, J.S.; Nagano, Y. Critically Evaluated Thermochemical Properties of Polycyclic Aromatic Hydrocarbons. J. Phys. Chem. Ref. Data 2008, 37, 1855–1996. [Google Scholar] [CrossRef]
  42. Pedley, J.B.; Naylor, R.D.; Kirby, S.P. Thermochemical Data of Organic Compounds; Chapman and Hall: New York, NY, USA, 1986; pp. 1–792. [Google Scholar]
  43. Oliver, G.D.; Eaton, M.; Huffman, H.M. The Heat Capacity, Heat of Fusion and Entropy of Benzene 1. J. Am. Chem. Soc. 1948, 70, 1502–1505. [Google Scholar] [CrossRef]
  44. Ambrose, D.; Ellender, J.; Sprake, C.H.; Townsend, R. Thermodynamic Properties of Organic Oxygen Compounds XLIII. Vapour Pressures of Some Ethers. J. Chem. Thermodyn. 1976, 8, 165–178. [Google Scholar] [CrossRef]
  45. Chirico, R.D.; Steele, W.V. Thermodynamic Properties of Diphenylmethane. J. Chem. Eng. Data 2005, 50, 1052–1059. [Google Scholar] [CrossRef]
  46. Chirico, R.D.; Knipmeyer, S.E.; Nguyen, A.; Steele, W.V. The Thermodynamic Properties to the Temperature 700 K of Naphthalene and of 2,7-Dimethylnaphthalene. J. Chem. Thermodyn. 1993, 25, 1461–1494. [Google Scholar] [CrossRef]
  47. Ruehrwein, R.A.; Huffman, H.M. Thermal Data. XVII. The Heat Capacity, Entropy, and Free Energy of Formation of Cyclohexane. A New Method of Heat Transfer in Low Temperature Calorimetry. J. Am. Chem. Soc. 1943, 65, 1620–1625. [Google Scholar] [CrossRef]
  48. Verevkin, S.P.; Samarov, A.A.; Vostrikov, S.V. Thermodynamics of Reversible Hydrogen Storage: Are Methoxy-Substituted Biphenyls Better through Oxygen Functionality? Hydrogen 2023, 4, 862–880. [Google Scholar] [CrossRef]
  49. McCullough, J.P.; Finke, H.L.; Messerly, J.F.; Todd, S.S.; Kincheloe, T.C.; Waddington, G. The Low-Temperature Thermodynamic Properties of Naphthalene, l-Methylnaphthalene, 2-Methylnaphthalene, 1,2,3,4-Tetrahydronaphthalene,Trans-Decahydronaphthalene and Cis-Decahydronaphthalene. J. Phys. Chem. 1957, 61, 1105–1116. [Google Scholar] [CrossRef]
Figure 1. The dehydrogenation reactions of liquid organic hydrogen carriers.
Figure 1. The dehydrogenation reactions of liquid organic hydrogen carriers.
Oxygen 04 00015 g001
Figure 2. Calculating the enthalpy of vaporisation, l g H m o (298.15 K) or f H m o (g, 298.15 K) of for LOHC pairs dibenzyl ether/dicyclohexylmethyl ether and 1,3-diphenyl-propane/1,3-dicyclohexyl-propane. The numerical values of contributions required for calculations are given in Table S8.
Figure 2. Calculating the enthalpy of vaporisation, l g H m o (298.15 K) or f H m o (g, 298.15 K) of for LOHC pairs dibenzyl ether/dicyclohexylmethyl ether and 1,3-diphenyl-propane/1,3-dicyclohexyl-propane. The numerical values of contributions required for calculations are given in Table S8.
Oxygen 04 00015 g002
Table 1. Compilation of combustion results at T = 298.15 K ( = 0.1 MPa) obtained for benzyl phenyl ether and dibenzyl ether a.
Table 1. Compilation of combustion results at T = 298.15 K ( = 0.1 MPa) obtained for benzyl phenyl ether and dibenzyl ether a.
Compound c u o (cr/liq) c H m o (cr/liq) f H m o (cr/liq)exp
J·g−1kJ·mol−1kJ·mol−1
benzyl phenyl ether (cr)−36,691.3 ± 5.1−6765.9 ± 2.4−64.7 ± 3.0
dibenzyl ether (liq)−37,589.6 ± 3.8−7459.9 ± 2.3−50.0 ± 2.9
a Uncertainties related to the combustion experiments were estimated according to the procedure recommended by Olofsson [23]. The auxiliary information for the combustion experiments is given in Table S2.
Table 2. Compilation of the standard molar enthalpies of vaporisation of aromatic and aliphatic ethers and hydrocarbons (in kJ·mol−1).
Table 2. Compilation of the standard molar enthalpies of vaporisation of aromatic and aliphatic ethers and hydrocarbons (in kJ·mol−1).
CompoundsMethod aT-Range/
K
l g H m o
Tav
l g H m o
298.15 K b
benzyl phenyl ethern/a368.6–560.259.2 ± 1.071.3 ± 2.6[26]
946-80-5T318.3–363.270.6 ± 0.373.8 ± 0.4Table S3
BP358–56162.8 ± 0.674.7 ± 2.4Table S4
73.7 ± 0.5 caverage
GA298.1575.8 ± 1.5this work
dibenzyl ethern/a275–41750.1 ± 1.0(54.2 ± 1.3)[27]
103-50-4n/a413–56157.4 ± 1.073.8 ± 3.4[27]
TGA392–42569.4 ± 1.078.9 ± 2.1[28]
T309.4 ± 358.175.6 ± 0.478.6 ± 0.6Table S3
BP365–57167.0 ± 0.479.5 ± 2.5Table S4
78.5 ± 0.5average
GA298.1578.2 ± 1.5this work
dicyclohexyl etherGA68.3 ± 1.5[29]
4645-15-2BP348–51757.4 ± 1.168.1 ± 2.4Table S4
68.2 ± 1.3 caverage
1,2-diphenylethaneIP333.2–413.264.0 ± 0.169.2 ± 1.0[30]
103-29-7BP338–55859.6 ± 0.568.8 ± 1.9Table S4
69.1 ± 0.9 caverage
GA298.1571.2 ± 1.5this work
1,2-dicyclohexylethaneBP348–54756.7 ± 0.868.2 ± 2.4Table S4
3321-50-4GA298.1567.2 ± 1.5this work
Tb298.1567.9 ± 2.0Equation (6)
67.6 ± 1.1 caverage
cyclohexylmethoxy–cyclohexaneGA298.1569.7 ± 1.5this work
118161-77-6
1,3-diphenylpropanen/a342–57760.0 ± 1.071.4 ± 2.5[27]
1081-75-0SC298.1574.2 ± 2.0[31]
BP342–57663.4 ± 0.773.8 ± 2.2Table S4
73.3 ± 1.3average
GA298.1575.7 ± 1.5this work
1,3-dicyclohexylpropaneBP382–56558.1 ± 0.572.9 ± 3.0Table S4
3178-24-3GA298.1571.8 ± 1.5this work
Tb298.1571.3 ± 2.0Equation (6)
71.8 ± 1.1 caverage
dicyclohexylmethyl etherBP356–57459.2 ± 2.173.9 ± 3.6Table S4
14315-63-0GA298.1572.8 ± 1.5this work
73.0 ± 1.4 caverage
2-methoxy-naphthaleneBP365–54959.8 ± 1.168.9 ± 2.1[32]
93-04-9T348.2–378.265.4 ± 0.869.1 ± 0.9[32]
T348.7–372.466.1 ± 0.569.6 ± 0.6Table S3
PhT69.2 ± 0.8Table S5
69.4 ± 0.4 caverage
2-methoxy-decalin (liq)GA58.4 ± 1.0[32]
55473-38-6Tb58.8 ± 1.0[32]
BP349–50949.6 ± 1.758.7 ± 2.5Table S4
58.6 ± 0.7 caverage
a Methods: n/a = method is not available; T = transpiration method; SC = solution calorimetry; Tb = from correlation with the normal boiling temperatures; BP = estimated from boiling points at different pressures (see Table S4); PhT = from the consistency of phase transitions (see Table S5); GA = calculated according to the group additivity (see text). b Vapour pressures available in the literature were treated using Equations (S2) and (S3) with the help of heat capacity differences from Tables S6 and S7 to calculate the enthalpies of vaporisation at 298.15 K. Uncertainties of the sublimation/vaporisation enthalpies U( l g H m o ) are the expanded uncertainties (0.95 level of confidence). They include uncertainties from the fitting equation and uncertainties from temperature adjustment to T = 298.15 K. Uncertainties in the temperature adjustment of vaporisation enthalpies to the reference temperature T = 298.15 K are estimated to account for 20% of the total adjustment. c Weighted mean value (uncertainties were taken as the weighting factor). Values given in bold are recommended for further thermochemical calculations.
Table 3. Results of quantum-chemical calculations with the G3MP2 method for aromatic compounds and their hydrogenated products at T = 298.15 K ( = 0.1 MPa). All data are in kJ·mol−1.
Table 3. Results of quantum-chemical calculations with the G3MP2 method for aromatic compounds and their hydrogenated products at T = 298.15 K ( = 0.1 MPa). All data are in kJ·mol−1.
CASCompound f H m o (g, AT)G3MP2 a f H m o (g)exp b f H m o (g)QC cΔ d
496-16-2dihydrobenzofuran−57.5−47.0 ± 1.4 [36]−46.3−0.7
496-14-0dihydroisobenzofuran−36.9−30.1 ± 1.0 [37]−25.3−4.8
493-05-0isochromane−78.2−63.1 ± 1.1 [38]−67.44.3
2216-69-51-methoxynaphthalene−19.0−6.7 ± 1.8 [32]−7.10.4
93-04-92-methoxynaphthalene−20.3−6.5 ± 1.4 [32]−8.41.9
493-08-3dihydrobenzopyran−96.2−82.4 ± 1.2 [38]−85.83.4
21720-89-81-methoxydecalin−333.6−328.7 ± 3.5 [32]−327.8−0.9
55473-38-62-methoxydecalin−334.6−329.5 ± 3.5 [32]−328.8−0.7
101-84-8diphenyl ether38.650.9 ± 1.4 [29]51.6−0.7
4645-15-2dicyclohexyl ether−361.4-−356.1
946-80-5phenyl benzyl ether20.730.6 ± 3.1 [Table 4]33.4−2.8
118161-77-6cyclohexylmethoxy–cyclohexane−377.0-−372.0
103-50-4dibenzyl ether15.028.5 ± 3.0 [Table 4]27.60.9
14315-63-0dicyclohexylmethyl ether−391.5-−386.8
a Calculated by G3MP2 and atomisation reaction with the expanded uncertainties of ±4.1 kJ·mol−1 [21]. b Uncertainties in this table are expressed as two times the standard deviation. c Calculated according to Equation (8). d Difference between column 3 and column 4.
Table 4. Comparison of experimental and theoretical thermochemical data at T = 298.15 K ( = 0.1 MPa) for aromatic ethers (in kJ·mol−1) a.
Table 4. Comparison of experimental and theoretical thermochemical data at T = 298.15 K ( = 0.1 MPa) for aromatic ethers (in kJ·mol−1) a.
Compound f H m o (cr/liq)exp b c r / l g H m o  c c r l H m o  d f H m o (g)exp e f H m o (g)QC f
benzyl phenyl ether (cr)−64.7 ± 3.073.7 ± 0.521.6 ± 0.430.6 ± 3.133.4 ± 4.1
dibenzyl ether (liq)−50.0 ± 2.978.5 ± 0.528.5 ± 3.027.6 ± 4.1
2-methoxy-naphthalene (cr)−96.6 ± 1.3 [32]90.1 ± 0.5 [32]−6.5 ± 1.4−8.4 ± 4.1
a Uncertainties correspond to expanded uncertainties of the mean (0.95 level of confidence). b From Table 1. c The evaluated values from Table 2. d From Table S5. e The sum of columns 2 to 4. f Calculated with the G3MP2 method using the corrected atomisation reactions (see Table 3, column 5).
Table 5. Compilation of thermochemical data at T = 298.15 K ( = 0.1 MPa) for LOHC counterparts (in kJ·mol−1) a.
Table 5. Compilation of thermochemical data at T = 298.15 K ( = 0.1 MPa) for LOHC counterparts (in kJ·mol−1) a.
Compound; CAS f H m o (g) b l g H m o  c f H m o (liq) d
diphenyl–methane; 101-81-5164.7 ± 0.9 [40]68.1 ± 0.3 [40]96.6 ± 0.8 [40]
dicyclohexyl–methane; 3178-23-2 −242.4 ± 2.0 [34]64.7 ± 1.0 [34]−307.1 ± 2.2 [34]
diphenyl ether, 101-84-8 50.9 ± 1.4 [29]66.7 ± 0.2 [29]−15.8 ± 1.4 [29]
dicyclohexyl ether; 4645-15-2−356.1 ± 4.1 [Table 3]68.2 ± 1.3 −424.3 ± 4.3
1,2-diphenyl-ethane; 103-29-7133.7 ± 4.1 [Table S10]69.1 ± 0.9 64.6 ± 4.2
1,2-dicyclohexyl-ethane; 3321-50-4−261.0 ± 4.1 [Table S11]67.6 ± 1.1 −328.6 ± 4.2
benzyl phenyl ether; 946-80-530.6 ± 3.1 [Table 4]73.7 ± 0.5 −43.1 ± 3.1
cyclohexylmethoxy–cyclohexane; 118161-77-6−372.0 ± 4.1 [Table 3]69.7 ± 1.5 −441.7 ± 4.4
1,3-diphenyl-propane; 1081-75-0114.7 ± 4.1 [Table S10]73.3 ± 1.3 41.4 ± 4.3
1,3-dicyclohexyl-propane; 3178-24-3−281.6 ± 4.1 [Table S11]71.8 ± 1.1 −353.4 ± 4.2
dibenzyl ether; 103-50-4−50.0 ± 2.9 [Table 4]78.5 ± 0.5 −128.5 ± 2.9
dicyclohexylmethyl ether; 14315-63-0−386.8 ± 4.1 [Table 3]73.0 ± 1.4 −459.8 ± 4.3
naphthalene; 91-20-3 150.6 ± 1.5 [41]55.4 ± 1.4 [41]94.9 ± 1.6 [41]
trans-decalin; 493-02-7 −182.0 ± 1.0 [42]48.6 ± 0.2 [42]−230.6 ± 1.0 [42]
2-methoxy-naphthalene; 93-04-9−6.5 ± 1.4 [32]69.4 ± 0.4 −75.9 ± 1.5
trans-2-methoxy-decalin; 55473-38-6 −329.5 ± 3.6 [32]58.6 ± 0.7 −388.1 ± 3.6
a Uncertainties correspond to expanded uncertainties of the mean (0.95 level of confidence). b Calculated with the G3MP2 method using the corrected atomisation reactions (see Table 3, Tables S10 and S11) or taken from the original literature. c From Table 2 or from original literature. d Calculated as the difference between columns 2 and 3 in this table.
Table 6. Calculation of the liquid-phase reaction enthalpies, r H m o (liq), of the dehydrogenation reactions R1 to R10 in Figure 1, at T = 298.15 K ( = 0.1 MPa, in kJ·mol−1).
Table 6. Calculation of the liquid-phase reaction enthalpies, r H m o (liq), of the dehydrogenation reactions R1 to R10 in Figure 1, at T = 298.15 K ( = 0.1 MPa, in kJ·mol−1).
Compound f H m o (liq)HL a f H m o (liq)HR b r H m o (liq) c r H m o (liq)/H2 d
benzene49.0 ± 0.9 [Table S12]−156.4 ± 0.8 [Table S12]205.468.5
methoxy–benzene−116.9 ± 0.7 [Table S12]−313.7 ± 3.5 [Table S12]196.865.6
diphenyl–methane96.6 ± 0.8−307.1 ± 2.2403.767.3
diphenyl ether−15.8 ± 1.4−424.3 ± 4.3408.568.1
1,2-diphenyl-ethane64.6 ± 4.2−328.6 ± 4.2393.265.5
benzyl phenyl ether−43.1 ± 3.1−441.7 ± 4.4398.666.4
1,3-diphenyl-propane41.4 ± 4.3−353.4 ± 4.2394.865.8
dibenzyl ether−50.0 ± 2.9−459.8 ± 4.3409.868.3
naphthalene94.9 ± 1.6−230.6 ± 1.0325.565.1
2-methoxy-naphthalene−75.9 ± 1.5−388.1 ± 3.6312.262.4
a Liquid-phase enthalpies of formation of of aromatic ethers and hydrocarbons (hydrogen-lean (HL) counterparts of the LOHC system) from Table 5. b Liquid-phase enthalpies of formation of perhydrogenated aromatic ethers and hydrocarbons (hydrogen-rich (HR) counterparts of the LOHC system) from Table 5. c Calculated according to Hess´s Law applied to the reactions R1 to R10 shown in Figure 1. d Reaction enthalpy per mole H2.
Table 7. Compilation of the standard molar entropies of vaporisation, l g S m o , and the absolute standard molar entropies, S m o (g or liq), of reactants of reactions R1–R10 (all values at T = 298.15 K in J·mol−1·K−1).
Table 7. Compilation of the standard molar entropies of vaporisation, l g S m o , and the absolute standard molar entropies, S m o (g or liq), of reactants of reactions R1–R10 (all values at T = 298.15 K in J·mol−1·K−1).
Compound l g S m o  a S m o (g) b S m o (liq) c
benzene173.3 [43]
methoxy–benzene111.8 ± 0.1 [44]346.6234.8
diphenyl–methane140.5± 0.5 [45]440.3299.8
diphenyl ether134.5± 4.3 [29]439.8305.3
1,2-diphenyl-ethane134.9 ± 1.8466.4331.5
benzyl-phenyl ether142.2 ± 1.1456.9314.7
1,3-diphenyl-propane140.9 ± 2.5498.9358.0
dibenzyl ether151.8 ± 1.3502.4350.6
naphthalene217.6 [46]
2-methoxy-naphthalene137.9 ± 1.6401.5263.6
cyclohexane204.4 [47]
methoxy–cyclohexane106.7 ± 3.2 [48]376.8270.1
dicyclohexyl-methane136.3 ± 2.8 [34]470.5334.2
1,2-dicyclohexyl-ethane137.1 ± 2.6510.5373.4
1,3-dicayclohexyl-propane144.1 ± 1.7544.4400.3
dicyclohexyl ether142.9 ± 3.9475.1332.2
cyclohexylmethoxy–cyclohexane145.0 ± 5.0 d511.1366.1
dicyclohexylmethyl ether147.9 ± 7.1539.2391.3
trans-decalin264.9 [49]
2-methoxy-decalin123.8 ± 5.7446.3322.5
a From the results, the vapour pressure measurements are compiled in Tables S3 and S4. b Calculated with the G3MP2 method. c Calculated according to Equation (14) using entries from columns 2 and 3 from this table. d Assessed. The changes in the standard molar entropies of the dehydrogenation reactions R1 to R10 were calculated according to Equation (12) using the standard molar entropies of the reactants from Table 7, as follows.
Table 8. Liquid-phase thermodynamic properties of the dehydrogenation of LOHC systems.
Table 8. Liquid-phase thermodynamic properties of the dehydrogenation of LOHC systems.
Reactants r H m o  a
(298 K)
r S m o  b
(298 K)
T × r S m o
(298 K)
r G m o  c
(298 K)
r G m o  d
(400 K)
r G m o  e
(500 K)
kJ·mol−1J·mol−1·K−1kJ·mol−1kJ·mol−1kJ·mol−1kJ·mol−1
benzene205.4360.9107.697.860.121.1
cyclohexane68.5120.335.932.620.07.0
methoxy–benzene197.7356.7106.491.323.914.6
methoxy–cyclohexane65.9118.935.530.418.04.9
diphenyl–methane403.7749.7223.5180.2102.222.1
dicyclohexyl–methane67.3125.037.330.017.03.7
diphenyl–ether408.5757.2225.8182.7103.822.5
dicyclohexyl ether68.1126.237.630.517.33.8
1,2-diphenyl-ethane393.2742.2221.3171.994.514.5
1,2-dicyclohexyl-ethane65.5123.736.928.715.72.4
benzyl phenyl ether398.6732.7218.5180.2103.624.3
Cyclohexylmethoxy–cyclohexane66.4122.136.430.017.34.1
1,3-diphenyl-propane394.8741.8221.2173.696.216.3
1,3-dicyclohexyl-propane65.8123.636.928.916.02.7
dibenzyl ether409.8743.4221.6188.2110.530.2
dicyclohexylmethyl ether68.3123.936.931.418.45.0
naphthalene325.5606.1180.7144.881.415.7
trans-decalin65.1121.236.129.016.33.1
2-methoxy-naphthalene312.2594.5177.3134.972.88.2
trans-2-methoxy-decalin62.4118.935.527.014.61.6
a Calculated according to Equation (11) from standard molar enthalpies of formation of reactants from Table 6. Values highlighted in bold are related to the amount of hydrogen released and given in kJ·mol−1/H2 units. More details are shown in Table S13. b Calculated according to Equation (12) from standard molar entropies of reactants from Table 7. Values highlighted in bold are related to the amount of hydrogen released and given in J·mol−1·K−1/H2 units. c Calculated according to Equation (1) from results given in columns 2 and 3 and referenced to 298 K. d The r G m o (298.15 K) value from column 4 was adjusted to 400 K, as shown in ESI. e The r G m o (298.15 K) value from column 4 was adjusted to 500 K, as shown in ESI.
Table 9. The total Gibbs free reaction enthalpies (in kJ·mol−1) and equilibrium constants of dehydrogenation of aliphatic cyclic hydrocarbons and ethers in the liquid phase, calculated according to Equation (1).
Table 9. The total Gibbs free reaction enthalpies (in kJ·mol−1) and equilibrium constants of dehydrogenation of aliphatic cyclic hydrocarbons and ethers in the liquid phase, calculated according to Equation (1).
LOHC System r G m o
(298 K)
r G m o
(400 K)
r G m o
(500 K)
Ka
(298 K)
Ka
(400 K)
Ka
(500 K)
cyclohexane97.860.121.17.2 × 10−181.4 × 10−86.2 × 10−3
methoxy–cyclohexane91.323.914.69.9 × 10−177.6 × 10−43.0 × 10−2
dicyclohexyl–methane180.2102.222.12.6 × 10−324.5 × 10−144.9 × 10−3
dicyclohexyl ether182.7103.822.59.4 × 10−332.8 × 10−144.5 × 10−3
1,2-dicyclohexyl-ethane171.994.514.57.4 × 10−314.6 × 10−133.1 × 10−2
Cyclohexylmethoxy–cyclohexane180.2103.624.32.6 × 10−323.0 × 10−142.9 × 10−3
1,3-dicyclohexyl-propane173.696.216.33.7 × 10−312.7 × 10−132.0 × 10−2
dicyclohexylmethyl ether188.2110.530.21.0 × 10−333.7 × 10−157.0 × 10−4
decalin144.881.415.74.1 × 10−262.3 × 10−112.3 × 10−2
2-methoxy-decalin134.972.88.22.3 × 10−243.1 × 10−101.4 × 10−1
Table 10. Calculation of the equilibrium temperature, T e q , for dehydrogenation/hydrogenation reactions with aliphatic cyclic hydrocarbons and ethers in the liquid phase calculated using the G3MP2 method ( = 0.1 MPa).
Table 10. Calculation of the equilibrium temperature, T e q , for dehydrogenation/hydrogenation reactions with aliphatic cyclic hydrocarbons and ethers in the liquid phase calculated using the G3MP2 method ( = 0.1 MPa).
Reactants r H m o  a r S m o  b T e q  c
kJ·mol−1J·mol−1·K−1K
cyclohexane = benzene + 3 × H2205.4360.9569
methoxy–cyclohexane = methoxy–benzene + 3 × H2197.7356.7554
Dicyclohexyl–methane = diphenyl–methane + 6 × H2403.7749.7538
dicyclohexyl ether = diphenyl ether + 6 × H2408.5757.2539
1,2-dicyclohexyl-ethane = 1,2-diphenyl-ethane + 6 × H2393.2742.2530
cyclohexylmethoxy–cyclohexane = benzyl phenyl ether + 6 × H2398.6732.7544
1,3-dicyclohexyl-propane = 1,3-diphenyl-propane + 6 × H2394.8741.8532
dibenzyl ether + 6 × H2 = dicyclohexylmethyl ether409.8743.4551
trans-decalin = naphthalene + 5 × H2325.5606.1537
trans-2-methoxy-decalin = 2-methoxy-naphthalene + 5 × H2312.2594.5525
a The total reaction enthalpies from Table 9. b The total reaction entropies from Table 9. c The equilibrium temperature calculated according to Equation (13).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Verevkin, S.P.; Samarov, A.A.; Vostrikov, S.V. Does the Oxygen Functionality Really Improve the Thermodynamics of Reversible Hydrogen Storage with Liquid Organic Hydrogen Carriers? Oxygen 2024, 4, 266-284. https://doi.org/10.3390/oxygen4030015

AMA Style

Verevkin SP, Samarov AA, Vostrikov SV. Does the Oxygen Functionality Really Improve the Thermodynamics of Reversible Hydrogen Storage with Liquid Organic Hydrogen Carriers? Oxygen. 2024; 4(3):266-284. https://doi.org/10.3390/oxygen4030015

Chicago/Turabian Style

Verevkin, Sergey P., Artemiy A. Samarov, and Sergey V. Vostrikov. 2024. "Does the Oxygen Functionality Really Improve the Thermodynamics of Reversible Hydrogen Storage with Liquid Organic Hydrogen Carriers?" Oxygen 4, no. 3: 266-284. https://doi.org/10.3390/oxygen4030015

APA Style

Verevkin, S. P., Samarov, A. A., & Vostrikov, S. V. (2024). Does the Oxygen Functionality Really Improve the Thermodynamics of Reversible Hydrogen Storage with Liquid Organic Hydrogen Carriers? Oxygen, 4(3), 266-284. https://doi.org/10.3390/oxygen4030015

Article Metrics

Back to TopTop