Thermodynamics of Reversible Hydrogen Storage: Are Methoxy-Substituted Biphenyls Better through Oxygen Functionality?

: The reversible hydrogenation/dehydrogenation of aromatic molecules, known as liquid organic hydrogen carriers, is considered as an attractive option for the safe storage and release of elemental hydrogen. The recently reported efﬁcient synthetic routes to obtain methoxy-biphenyls in high yield make them promising candidates for hydrogen storage. In this work, a series of methoxy-substituted biphenyls and their structural parent compounds were studied. The absolute vapour pressures were measured using the transpiration method and the enthalpies of vaporisation/sublimation were determined. We applied a step-by-step procedure including structure–property correlations and quantum chemical calculations to evaluate the quality of thermochemical data on the enthalpies of phase transitions and enthalpies of formation of the studied methoxy compounds. The data sets on thermodynamic properties were evaluated and recommended for calculations in chemical engineering. A thermodynamic analysis of chemical reactions based on methoxy-biphenyls in the context of hydrogen storage was carried out and the energetics of these reactions were compared with the energetics of reactions of common LOHCs. The inﬂuence of the position of the methoxy groups in the rings on the enthalpies of the reactions relevant for hydrogen storage was discussed.


Introduction
The sustainable development of technologies using renewable energy sources such as wind, solar, biomass, or geothermal instead of traditional fossil energy sources is the obvious trend of the last decades.Currently, one of the most promising concepts for the alternative energy sources is based on the conversion of the redundant electrical energy obtained from the sun, wind, etc., into elemental hydrogen to be used as an energy carrier.Conventional compression or liquefaction of hydrogen is challenging and not safe enough.In contrast, the LOHCs (or chemical hydrides) have the advantage of storing a sufficient amount of hydrogen at ambient temperatures and releasing it through reversible catalytic dehydrogenation processes over several cycles.
Generally, all compounds containing the double bonds in the molecule can be considered as a potential LOHC.However, if hydrogenation of the double bonds is the thermodynamically favourable reaction, reverse dehydrogenation is thermodynamically very unfavourable and requires high temperatures to shift the equilibrium toward sufficient hydrogen release.The dehydrogenation temperature is directly related to the chemical structures of the hydrogen-lean and hydrogen-rich counterparts of the LOHC system.Understanding the general regularities between the structures of the LOHC counterparts and their thermodynamic properties is the only way to find the optimal compounds for practical applications.The present work extends our previous efforts [1] to establish these general regularities via thermochemical studies of LOHC systems based on methoxy-substituted biphenyls, as shown in Figure 1.
Hydrogen 2023, 4, FOR PEER REVIEW 2 unfavourable and requires high temperatures to shift the equilibrium toward sufficient hydrogen release.The dehydrogenation temperature is directly related to the chemical structures of the hydrogen-lean and hydrogen-rich counterparts of the LOHC system.Understanding the general regularities between the structures of the LOHC counterparts and their thermodynamic properties is the only way to find the optimal compounds for practical applications.The present work extends our previous efforts [1] to establish these general regularities via thermochemical studies of LOHC systems based on methoxy-substituted biphenyls, as shown in Figure 1.The recently reported efficient synthetic routes to obtain these materials in high yield [2][3][4] make them promising candidates for hydrogen storage.The practically relevant LOHC systems and chemical reactions of the reversible hydrogen storage are given in Figure 2. The perhydrogenated methoxy-bicyclohexanes formed in the Reaction 1 may be involved in the transfer hydrogenation (TH) reactions (see Reaction 2).The transfer hydrogenation reactions can be considered as an attractive alternative to the reversible catalytic hydrogenation/dehydrogenation systems, since these reactions occur under mild conditions without expensive high-pressure reactors [5].The TH reactions are very suitable for feeding fuel cells, e.g., for the reduction of acetone to the corresponding iso-propanol [6], using perhydro-dibenzyltoluene and platinum-on-silica (Pt/SiO2) catalyst [7].Oxygenated LOHC systems have been shown to release hydrogen efficiently, especially at low temperatures, and have a high storage capacity, which makes them attractive for technical applications [8].Consequently, alkoxy-naphthalenes/alkoxy-bicyclohexanes are also expected to be efficient as LOHC systems.

Methodology
The hydrogenation/dehydrogenation reactions proceed under pressure, and the reactants are generally liquid under the reaction conditions (with the exception of hydrogen).The hydrogenation/dehydrogenation reaction enthalpies, ∆  m o , for any liquid-phase The recently reported efficient synthetic routes to obtain these materials in high yield [2][3][4] make them promising candidates for hydrogen storage.The practically relevant LOHC systems and chemical reactions of the reversible hydrogen storage are given in Figure 2.
Hydrogen 2023, 4, FOR PEER REVIEW 2 unfavourable and requires high temperatures to shift the equilibrium toward sufficient hydrogen release.The dehydrogenation temperature is directly related to the chemical structures of the hydrogen-lean and hydrogen-rich counterparts of the LOHC system.Understanding the general regularities between the structures of the LOHC counterparts and their thermodynamic properties is the only way to find the optimal compounds for practical applications.The present work extends our previous efforts [1] to establish these general regularities via thermochemical studies of LOHC systems based on methoxy-substituted biphenyls, as shown in Figure 1.The recently reported efficient synthetic routes to obtain these materials in high yield [2][3][4] make them promising candidates for hydrogen storage.The practically relevant LOHC systems and chemical reactions of the reversible hydrogen storage are given in Figure 2. The perhydrogenated methoxy-bicyclohexanes formed in the Reaction 1 may be involved in the transfer hydrogenation (TH) reactions (see Reaction 2).The transfer hydrogenation reactions can be considered as an attractive alternative to the reversible catalytic hydrogenation/dehydrogenation systems, since these reactions occur under mild conditions without expensive high-pressure reactors [5].The TH reactions are very suitable for feeding fuel cells, e.g., for the reduction of acetone to the corresponding iso-propanol [6], using perhydro-dibenzyltoluene and platinum-on-silica (Pt/SiO2) catalyst [7].Oxygenated LOHC systems have been shown to release hydrogen efficiently, especially at low temperatures, and have a high storage capacity, which makes them attractive for technical applications [8].Consequently, alkoxy-naphthalenes/alkoxy-bicyclohexanes are also expected to be efficient as LOHC systems.

Methodology
The hydrogenation/dehydrogenation reactions proceed under pressure, and the reactants are generally liquid under the reaction conditions (with the exception of hydrogen).The hydrogenation/dehydrogenation reaction enthalpies, ∆  m o , for any liquid-phase The perhydrogenated methoxy-bicyclohexanes formed in the Reaction 1 may be involved in the transfer hydrogenation (TH) reactions (see Reaction 2).The transfer hydrogenation reactions can be considered as an attractive alternative to the reversible catalytic hydrogenation/dehydrogenation systems, since these reactions occur under mild conditions without expensive high-pressure reactors [5].The TH reactions are very suitable for feeding fuel cells, e.g., for the reduction of acetone to the corresponding iso-propanol [6], using perhydro-dibenzyltoluene and platinum-on-silica (Pt/SiO 2 ) catalyst [7].Oxygenated LOHC systems have been shown to release hydrogen efficiently, especially at low temperatures, and have a high storage capacity, which makes them attractive for technical applications [8].Consequently, alkoxy-naphthalenes/alkoxy-bicyclohexanes are also expected to be efficient as LOHC systems.

Methodology
The hydrogenation/dehydrogenation reactions proceed under pressure, and the reactants are generally liquid under the reaction conditions (with the exception of hydrogen).The hydrogenation/dehydrogenation reaction enthalpies, ∆ r H o m , for any liquid-phase hydrogen storage LOHC system are usually calculated from the standard molar liquid phase enthalpies of the formation of reactants, ∆ f H o m (liq), using Hess's Law.With increasing interest in reversible hydrogen storage, it is impractical to measure ∆ f H o m (liq) values for the numerous LOHC systems that hold promise for the future.Therefore, a judicious combination of experimental, empirical, and quantum chemical (QC) methods seems to be the only way forward in the thermodynamic analysis and selection of LOHC systems of industrial interest.In particular, the use of modern QC methods enables a significant acceleration in the screening of promising LOHC systems, since standard molar gas phase formation enthalpies, ∆ f H o m (liq), calculated using high-level methods nowadays reach a "chemical accuracy" of ±4-5 kJ•mol −1 [9].
The quantum chemical gas-phase enthalpies of formation and the practically important liquid phase enthalpies of formation are related via the textbook equation where  1) and a thermodynamic analysis of the LOHC system is performed.This roadmap helps to follow the logic of the combination of experimental and theoretical methods used in this work to evaluate the practical potential of the methoxybiphenyl/methoxy-bicyclohexyl LOHC systems.

Materials
The samples of methoxy-biphenyls used in this work were commercially available (see Table S1).They were additionally purified before the thermochemical measurements.Purities better than a 0.9998 mass fraction were determined using gas-liquid chromatography on the capillary column.Details on purification and purity determination are given in ESI.

Experimental and Theoretical Methods
The absolute vapour pressures over methoxy-biphenyls were measured using the transpiration method.A nitrogen stream was passed through the saturator filled with the sample at various constant temperatures (±0.1 K).To eliminate hydraulic resistance and create sufficient contact area between the sample and the nitrogen stream, the sample was spread over the small glass beads.After saturation equilibrium was reached, the material transported by the gas stream was collected in a cold trap and its amount determined using GC.The temperature dependences of the absolute vapour pressures for each compound were used to calculate thermodynamic functions of vaporisation/sublimation of methoxybiphenyls.A concise description and the necessary details are given in ESI.
The melting temperatures and enthalpies of fusion for 4-methoxy-biphenyl and 4,4dimethoxy-biphenyl were measured via DSC.The necessary experimental details are given in ESI.
The Gaussian 16 series [10] software was used for quantum chemical (QC) calculations.The H 298 enthalpies of the most stable conformers were calculated using the G4 method [11] and they were finally converted to the theoretical ∆ f H o m (g, 298.15 K) values and discussed.

Step I: Experimental Absolute Vapour Pressures
The vapour pressure-temperature dependencies for methoxy-biphenyls were measured using the transpiration method and fitted using the following equation: where R = 8.314462 J•K −1 •mol −1 is the molar gas constant, the reference pressure p ref = 1 Pa, a and b are adjustable parameters, the arbitrary temperature T 0 applied in Equation ( 2) was chosen to be T 0 = 298.2) are compiled in Table S2.The vapour pressures measured in this work and fitting parameters of Equation ( 2) are collected in Table 1.
The systematic vapour pressure measurements for 2-and 4-methoxy-biphenyls were performed for the first time.It turned out that the single experimental boiling temperatures at various reduced pressures found in the literature (see Table S3) can be used to support these new results.These boiling points are usually measured during purification of the synthesised product and are used to characterise the compound.As a rule, the pressures are determined using non-calibrated manometers and the boiling temperatures are usually monitored in the range of a few degrees.Recently [12], however, we have shown that even from such rough data, reasonable trends can generally be derived and used for comparison if no systematic vapour pressures are available.The boiling temperatures compiled for methoxy-biphenyls in Table S3 were measured at the elevated temperature and cannot be directly compared to the "low temperature" transpiration results.However, the ∆ g l H o m (298.15K) values derived from these literature results can be fairly compared to the results of this work, as shown in Section 4.2.

Step I: Experimental Standard Molar Enthalpies of Sublimation/Vaporisation from Vapour Pressure Measurements
The molar vaporisation/sublimation enthalpies of methoxy-biphenyls were derived from the temperature dependences of the vapour pressures, approximated by Equation (3) using the following equation: where b is one of the adjustable parameters of Equation (3).The vaporisation/sublimation enthalpies ∆ g cr,l H o m (298.15K) of the methoxy-biphenyls were calculated with help of Equations ( 2) and (3) with the ∆ g cr,l C o p,m values from Table S2.Details of the calculation procedures can be found in ESI.As mentioned in the introduction, in this work the auxiliary sets of the perhydrogenated methoxy-biphenyls (methoxy-bicyclohexanes) and structurally similar molecules of methoxy-iso-propylbenzenes and their perhydrogenated products (methoxy-iso-propyl-cyclohexanes) were studied using empirical methods.The boiling points available in the literature were collected for each compound and also treated using Equations ( 2) and (3) to calculate the vaporisation enthalpies at 298.15 K.The uniformly treated results are summarized for comparison in Table 2.
The standard molar vaporisation enthalpy for 2-methoxy-biphenyl as well as the standard molar vaporisation/sublimation enthalpies for 4-methoxy-biphenyl were derived for the first time.As can be seen from Table 2, the enthalpy of vaporisation calculated from the boiling points (BP) for 4-methoxybiphenyl, ∆ g l H o m (298.15K) = (76.5 ± 5.0) kJ•mol −1 , agrees within the experimental uncertainties with the more accurate experimental result, ∆ g l H o m (298.15K) = (79.9± 0.9) kJ•mol −1 , measured using the transpiration method.Also, the sufficient consistency of solid-liquid, liquid-gas, and solid-gas phase transitions for 4-methoxybiphenyl is shown in Table 3.  a Saturation temperature measured using standard uncertainty (u(T) = 0.1 K). b Mass of transferred sample condensed at T = 273 K. c Volume of nitrogen (u(V) = 0.005 dm 3 ) used to transfer m (u(m) = 0.0001 g) of the sample.Uncertainties are given as standard uncertainties.d T a is the temperature of the soap bubble meter used for the measurement of the gas flow.e Vapour pressure is given at temperature T, calculated from the m and the residual vapour pressure at the condensation temperature calculated using an iteration procedure.a Methods: T = transpiration method; J x = from correlation with Kovats indices; PhT = from consistency of phase transitions (see Table 3); CP = calculated according to the "centerpiece" approach (see text); T b = from correlation with the normal boiling temperatures; BP = from boiling points at various pressures (see Table S3).b Weighted mean value (uncertainties were taken as the weighting factor).Values given in bold are recommended for further thermochemical calculations.Indeed, the sublimation, vaporisation, and fusion enthalpies are related by Equation ( 4): when all enthalpies are adjusted to the common temperature T = 298.15K as explained in ESI.The 4-methoxybiphenyl is solid at room temperatures.The solid-liquid, liquid-gas, and solid-gas phase transitions for this compound are summarised in  2), demonstrating the consistency of the phase transition data for this compound.Similarly, the sublimation data for methyl 4,4 -dimethoxybiphenyl were calculated (see calculations in Table 3 and the evaluated value in Table 2).Unfortunately, for all other entries in Table 2, there are no experimental data for comparison.Therefore, any type of validation of ∆ g l H o m (298.15K) values is required for the compounds of interest listed in Table 2.

Step I: Validation of Vaporisation Enthalpies Using Structure-Property Correlations
Structure-property correlations are one of the basic principles of physical chemistry.This type of correlation helps to establish consistency in the new experimental results with the network of physicochemical data already available for similarly shaped molecules.

Correlation with the Kovats Indices J x
The Kovats retention indices, J x , are values that are readily measurable using gas chromatography (GC) [13].The ∆ g l H o m (298.15K) values in structural parent series of molecules usually correlate linearly with the Kovats indices.The Kovats indices of 3-and 4-substituted methoxy-biphenyls, methoxy-benzenes, and alkyl phenyl ethers on non-polar column Ultra 1 (see Table S4) were taken from the literature [14].
The collected S4 show a very good linear correlation with the corresponding J x values: Table S4 (column 6) shows that the differences between experimental values and those calculated according to Equation (5) values are mostly below 1.5 kJ•mol −1 .The "empirical" enthalpies of vaporisation derived from Equation (5) (see Table S4, column 5) are referred to as J x values in Table 2 and agree well with results obtained using other methods.
It has been found that the "ortho"-methoxy-substituted molecules (see Table S5) are in a separate line in this type of correlation: However, such behaviour is typical for the methoxy-substituted benzenes [15].
Hydrogen 2023, 4 Table S5 (column 6) shows that the disagreement between values calculated according to Equation ( 6) and the experimental values is well below 1.0 kJ•mol −1 .The "empirical" enthalpies of vaporisation derived from Equation (6) (see Table S5, column 5) also agree well with results obtained using other methods, as shown in Table 2.
Similarly, we correlated the ∆ g l H o m (298.15K) values of di-alkyl ethers with the corresponding J x values (see Tables S6 and S7) to obtain the auxiliary value of the enthalpy of vaporisation of methoxy-cyclohexane, which is needed for the thermochemical calculations in this work.The value ∆ g l H o m (298.15K) = 42.8 ± 1.0 kJ•mol −1 obtained in this way is in excellent agreement with the results obtained using other methods, as shown in Table 2.

Correlation Vaporisation Enthalpies of the Parent Structures
It is obvious that methoxy-biphenyls and the hydroxy-biphenyls (or phenyl-phenols) are structurally related.Therefore, a correlation of the experimental ∆ g l H o m (298.15K) values for the series of methoxy-biphenyls with the experimental vaporisation enthalpies for the series of hydroxy-biphenyls was obtained.The vaporisation enthalpies for the latter series were well established in our earlier work (see Tables S8-S10).The summary of experimental data for this correlation is given in Table 4.  referred to as Jx values in Table 2 and agree well with results obtain ods.
It has been found that the "ortho"-methoxy-substituted molec in a separate line in this type of correlation: However, such behaviour is typical for the methoxy-substituted be Table S5 (column 6) shows that the disagreement between val ing to Equation ( 6) and the experimental values is well below 1.0 kJ• enthalpies of vaporisation derived from Equation (6) (see Table S5 well with results obtained using other methods, as shown in Table Similarly, we correlated the ∆ l g  m o (298.15K) values of di-alkyl sponding Jx values (see Tables S6 and S7) to obtain the auxiliary va vaporisation of methoxy-cyclohexane, which is needed for the the tions in this work.The value ∆ l g  m o (298.15K) = 42.8 ± 1.0 kJ•mol −1 o in excellent agreement with the results obtained using other metho 2.

Correlation Vaporisation Enthalpies of the Parent Structures
It is obvious that methoxy-biphenyls and the hydroxy-bipheny are structurally related.Therefore, a correlation of the experimenta ues for the series of methoxy-biphenyls with the experimental vapo the series of hydroxy-biphenyls was obtained.The vaporisation en series were well established in our earlier work (see Tables S8-S10 perimental data for this correlation is given in Table 4.  2 and agree well with results obtained using other ods. It has been found that the "ortho"-methoxy-substituted molecules (see However, such behaviour is typical for the methoxy-substituted benzenes [15].Table S5 (column 6) shows that the disagreement between values calculated a ing to Equation ( 6) and the experimental values is well below 1.0 kJ•mol −1 .The "emp enthalpies of vaporisation derived from Equation (6) (see Table S5, column 5) also well with results obtained using other methods, as shown in Table 2.
Similarly, we correlated the ∆ l g  m o (298.15K) values of di-alkyl ethers with the sponding Jx values (see Tables S6 and S7) to obtain the auxiliary value of the entha vaporisation of methoxy-cyclohexane, which is needed for the thermochemical c tions in this work.The value ∆ l g  m o (298.15K) = 42.8 ± 1.0 kJ•mol −1 obtained in this in excellent agreement with the results obtained using other methods, as shown in 2.

Correlation Vaporisation Enthalpies of the Parent Structures
It is obvious that methoxy-biphenyls and the hydroxy-biphenyls (or phenyl-ph are structurally related.Therefore, a correlation of the experimental ∆  (298.15ues for the series of methoxy-biphenyls with the experimental vaporisation enthalp the series of hydroxy-biphenyls was obtained.The vaporisation enthalpies for the series were well established in our earlier work (see Tables S8-S10).The summary perimental data for this correlation is given in Table 4. a Experimental data evaluated in Tables S8 and S9.b Experimental data evaluated in Table 2.
As is apparent from Table 4, the vaporisation enthalpies of 2-substituted bip are drastically lower compared to those of 3-and 4-substituted biphenyls.The vap tion enthalpies of 3-and 4-isomers are indistinguishable within their experim a Experimental data evaluated in Tables S8 and S9.b Experimental data evaluated in Table 2.
As is apparent from Table 4, the vaporisation enthalpies of 2-substituted biphenyls are drastically lower compared to those of 3-and 4-substituted biphenyls.The vaporisation enthalpies of 3-and 4-isomers are indistinguishable within their experimental uncertainties.Such behaviour is typical for ortho-, meta-, and para-substituted benzenes [15].The correlation of the vaporisation enthalpies, ∆ uncertainties.Such behaviour is typical for ortho-, meta-, and para-substituted benzenes [15].The correlation of the vaporisation enthalpies, ∆ l g  m o (298.15K), of hydroxy-biphenyls and methoxy-biphenyls exhibits a practically straight line (see Figure 3, left), which can be taken as an indicator of the internal consistency of the ∆ l g  m o (298.15K) values in both series.Structurally, the series of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes is quite comparable with the series of hydroxy-biphenyls and methoxy-biphenyls shown in Figure S1.The correlation of the ∆ l g  m o (298.15K) values of hydroxy-isopropylbenzenes and methoxy-iso-propylbenzenes (see Table 5) also exhibits a practically straight line (see Figure S1), and it is considered as proof for the consistency of evaluated Structurally, the series of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes is quite comparable with the series of hydroxy-biphenyls and methoxy-biphenyls shown in Figure S1.The correlation of the ∆ g l H o m (298.15K) values of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes (see Table 5) also exhibits a practically straight line (see Figure S1), and it is considered as proof for the consistency of evaluated vaporisation enthalpies of methoxy-iso-propylbenzenes with the reliable set of data for hydroxy-isopropylbenzenes (see Tables S8-S10).Structurally, the series of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes is quite comparable with the series of hydroxy-biphenyls and methoxy-biphenyls shown in Figure S1.The correlation of the ∆ l g  m o (298.15K) values of hydroxy-isopropylbenzenes and methoxy-iso-propylbenzenes (see Table 5) also exhibits a practically straight line (see Figure S1), and it is considered as proof for the consistency of evaluated vaporisation enthalpies of methoxy-iso-propylbenzenes with the reliable set of data for hydroxy-iso-propylbenzenes (see Tables S8-S10).a Experimental data evaluated in Tables S8-S10.b Experimental data evaluated in Table 2.

hydroxy-iso-propylbenzenes
An interesting energetic aspect in structurally related behaviour is revealed by comparing the quantitative differences in in ∆ l g  m o (298.15K) values: between para-and orthomethoxy biphenyls and between para-and ortho-methoxy iso-propylbenzenes and the corresponding differences between para-and ortho-hydroxy biphenyls and for the paraand ortho-hydroxy iso-propylbenzenes, as shown in Figure 4. Structurally, the series of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes is quite comparable with the series of hydroxy-biphenyls and methoxy-biphenyls shown in Figure S1.The correlation of the ∆ l g  m o (298.15K) values of hydroxy-isopropylbenzenes and methoxy-iso-propylbenzenes (see Table 5) also exhibits a practically straight line (see Figure S1), and it is considered as proof for the consistency of evaluated vaporisation enthalpies of methoxy-iso-propylbenzenes with the reliable set of data for hydroxy-iso-propylbenzenes (see Tables S8-S10).
An interesting energetic aspect in structurally related behaviour is revealed by comparing the quantitative differences in in ∆ l g  m o (298.15K) values: between para-and orthomethoxy biphenyls and between para-and ortho-methoxy iso-propylbenzenes and the corresponding differences between para-and ortho-hydroxy biphenyls and for the paraand ortho-hydroxy iso-propylbenzenes, as shown in Figure 4. a Experimental data evaluated in Tables S8-S10.b Experimental data evaluated in Table 2.
An interesting energetic aspect in structurally related behaviour is revealed by comparing the quantitative differences in in ∆ g l H o m (298.15K) values: between para-and orthomethoxy biphenyls and between para-and ortho-methoxy iso-propylbenzenes and the corresponding differences between para-and ortho-hydroxy biphenyls and for the paraand ortho-hydroxy iso-propylbenzenes, as shown in Figure 4.  between para-and ortho-methoxy-biphenyls and between para-and ortho-methoxy-iso-propylbenzenes (left) and the corresponding differences between para-and ortho-hydroxy-biphenyls and between para-and ortho-hydroxy-iso-propylbenzenes (right).Numerical data are given in Tables S8 and S10.All values are in kJ•mol −1 .
It can be seen from Figure 4 that the differences in ∆ l g  m o (298.15K) values between para-and ortho-methoxy-iso-propylbenzenes, Δ = −7.1 kJ•mol −1 , and the corresponding differences between para-and ortho-hydroxy-iso-propylbenzenes, Δ = −5.9kJ•mol −1 , are quite comparable in magnitude.It can be concluded that ortho-methoxy and ortho-hydroxy substituents have a similar effect on reducing the enthalpies of vaporisation compared to the corresponding para-methoxy and para-hydroxy substituents on isopropylbenzene.Similar behaviour was expected for the interplay of the methoxy-and hy- m (298.15K)-values: between para-and ortho-methoxy-biphenyls and between para-and ortho-methoxy-iso-propylbenzenes (left) and the corresponding differences between para-and ortho-hydroxy-biphenyls and between para-and orthohydroxy-iso-propylbenzenes (right).Numerical data are given in Tables S8 and S10.All values are in kJ•mol −1 .
It can be seen from Figure 4 that the differences in ∆ g l H o m (298.15K) values between para-and ortho-methoxy-iso-propylbenzenes, ∆ = −7.1 kJ•mol −1 , and the corresponding differences between para-and ortho-hydroxy-iso-propylbenzenes, ∆ = −5.9kJ•mol −1 , are quite comparable in magnitude.It can be concluded that ortho-methoxy and ortho-hydroxy substituents have a similar effect on reducing the enthalpies of vaporisation compared to the corresponding para-methoxy and para-hydroxy substituents on iso-propylbenzene.Similar behaviour was expected for the interplay of the methoxy-and hydroxy-substituents on the biphenyl, considering that iso-propyl and phenyl are not spatially significantly different.To our surprise, the differences in ∆ g l H o m (298.15K) values between para-and ortho-methoxy-biphenyls, ∆ = −21.1 kJ•mol −1 , and the corresponding differences between para-and ortho-hydroxy-biphenyls, ∆ = −12.9kJ•mol −1 , are significantly different.How can this discrepancy be explained?Most likely, the reason is the intra-molecular hydrogen bonding specific to 2-phenylphenol.Indeed, it is known that in the liquid phase of such molecules the intermolecular hydrogen bond is predominant.However, upon leaving the liquid phase, the hydroxyl of the 2-phenylphenol immediately forms an intra-molecular H bond to the phenyl ring, which is in close proximity.The energy of this bond reduces the total energy required to evaporate 1 mole of 2-phenylphenol.In contrast, no intra-or inter-hydrogen bonding is possible in 2-methoxybiphenyl, and the very low enthalpy of vaporisation of this compound is due to the very loose packing of the sterically curled molecules.This plausible explanation for the energetic differences observed between the similarly shaped molecules is an indirect validation of the experimental and empirical vaporisation enthalpies of the methoxy-biphenyls evaluated in Table 2.

Correlation with Normal Boiling Points T b
One of the quantitative methods of structure-property correlations is the correlation of vaporisation enthalpies with normal boiling points.The enthalpies of vaporisation of methoxy-cyclohexane and 1-methoxy-2-iso-propyl-cyclohexane are of interest for this work and their boiling temperatures are known.This type of correlation, therefore, makes it possible to obtain the desired property in an independent manner.The vaporisation enthalpies, ∆ g l H o m (298.15K), and normal boiling points, T b , of di-alkyl ethers, methoxycyclohexane, and 1-methoxy-2-iso-propyl-cyclohexane (see Table S11) correlated as follows: The "empirical" vaporisation enthalpy ∆ g l H o m (298.15K) = 42.9 ± 0.6 kJ•mol −1 of methoxycyclohexane derived from Equation (7) (see Table S11, column 5) agrees well with the values obtained for this compound using other methods, as shown in Table 2.This good agreement can again be considered as evidence of good consistency of the experimental data included in the correlation and provides confidence the value ∆ g l H o m (298.15K) = 50.7 ± 0.6 kJ•mol −1 calculated for 1-methoxy-2-iso-propyl-cyclohexane according to Equation (7).As is finally shown in Table 2, various structure-property correlations have definitely demonstrated the internal consistency of the datasets of methoxy-biphenyls, methoxy-iso-propylbenzenes, and methoxy-cyclohexane, evaluated in this table.

Step I: Calculation of Vaporisation Enthalpies Using the "Centerpiece" Group-Additivity Concept
The methoxy-di-cyclohexanes are involved in the hydrogenation/dehydrogenation reactions shown in Figure 2 as the hydrogen-rich counterpart of the LOHC system methoxybiphenyl/methoxy-di-cyclohexane.Unfortunately, the experimental thermochemical data for these species are completely lacking.Numerous group additivity methods have been developed to predict both of the energetic properties, including enthalpies of formation and enthalpies vaporisation [16].The idea behind conventional group additivity (GA) methods is to divide the structures of molecules into small groups and to attribute these groups to well-defined numerical contributions for the corresponding property of interest.For the prediction, a frame of the desired molecule is collected from the corresponding number and type of these contributions.Recently, we modified the GA method and proposed a socalled "centerpiece" (CP) approach [17].In this approach, a potentially large "centerpiece" molecule with a reliable enthalpy is selected.This "centerpiece" molecule should generally mimic the structure of other molecules of interest.Then, the required increments or groups (see Figure S2) are added to the "centerpiece" molecule, resulting in the frame of the desired molecule.
For example, bi-cyclohexane can serve as a perfect "centerpiece" for predicting the thermodynamic properties of methoxy-bi-cyclohexanes.The enthalpy of vaporisation (and also the enthalpy of formation) for this compound are well established (see Table S10).The idea is illustrated in Figure 5.To construct the 4-methoxy-bi-cyclohexane the contribution H(C6H)→OCH3 = 9.9 kJ•mol −1 (see Figure S2) should be added to the vaporisation enthalpy of bi-cyclohexane ∆ l g  m o (298.15K) = 59.4 ± 0.3 kJ•mol −1 (see Table S10).The resulting value ∆ l g  m o (298.15K) = 69.3± 1.5 kJ•mol −1 (see Table 2) can be used for thermochemical calculations according to the reactions shown in Figure 2. To construct the 4,4′di-methoxy-bi-cyclohexane the two contributions H(C6H)→OCH3 should be added to the vaporisation enthalpy of bi-cyclohexane as shown in Figure 5.
The construction of 2-methoxy-bi-cyclohexane is somewhat more demanding.For this molecule, an additional contribution for the steric interaction of methoxy-substituent and the adjacent cyclohexane ring is required (see Figure 6).Figure 6.Calculating the ortho-interaction between methoxy-group and cyclohexane ring required for calculating the enthalpy of vaporisation of 2-methoxy-bi-cyclohexane using the bi-cyclohexane as the "centerpiece".
The latter contribution was derived as follows.Since thermodynamic data for 2-methoxy-bi-cyclohexane are completely lacking, we used the 2-methoxy-iso-propyl-cyclohexane to model the required steric interaction, assuming that the spatial characteristics of the iso-propyl and cyclohexyl substituents are essentially similar.The calculation was based on the empirical value for 2-methoxy-iso-propyl-cyclohexane, ∆ l g  m o (298.15K) = 50.7 ± 0.6 kJ•mol −1 (see Table 2).The "strain-free" reference value 53.9 ± 1.5 kJ•mol −1 for the To construct the 4-methoxy-bi-cyclohexane the contribution H(C 6 H)→OCH 3 = 9.9 kJ•mol −1 (see Figure S2) should be added to the vaporisation enthalpy of bi-cyclohexane ∆  S10).The resulting value ∆ g l H o m (298.15K) = 69.3± 1.5 kJ•mol −1 (see Table 2) can be used for thermochemical calculations according to the reactions shown in Figure 2. To construct the 4,4 di-methoxy-bi-cyclohexane the two contributions H(C 6 H)→OCH 3 should be added to the vaporisation enthalpy of bi-cyclohexane as shown in Figure 5.
The construction of 2-methoxy-bi-cyclohexane is somewhat more demanding.For this molecule, an additional contribution for the steric interaction of methoxy-substituent and the adjacent cyclohexane ring is required (see Figure 6).To construct the 4-methoxy-bi-cyclohexane the contribution H(C6H)→OCH3 = 9.9 kJ•mol −1 (see Figure S2) should be added to the vaporisation enthalpy of bi-cyclohexane ∆ l g  m o (298.15K) = 59.4 ± 0.3 kJ•mol −1 (see Table S10).The resulting value ∆ l g  m o (298.15K) = 69.3± 1.5 kJ•mol −1 (see Table 2) can be used for thermochemical calculations according to the reactions shown in Figure 2. To construct the 4,4′di-methoxy-bi-cyclohexane the two contributions H(C6H)→OCH3 should be added to the vaporisation enthalpy of bi-cyclohexane as shown in Figure 5.
The construction of 2-methoxy-bi-cyclohexane is somewhat more demanding.For this molecule, an additional contribution for the steric interaction of methoxy-substituent and the adjacent cyclohexane ring is required (see Figure 6).Figure 6.Calculating the ortho-interaction between methoxy-group and cyclohexane ring required for calculating the enthalpy of vaporisation of 2-methoxy-bi-cyclohexane using the bi-cyclohexane as the "centerpiece".
The latter contribution was derived as follows.Since thermodynamic data for 2-methoxy-bi-cyclohexane are completely lacking, we used the 2-methoxy-iso-propyl-cyclohexane to model the required steric interaction, assuming that the spatial characteristics of the iso-propyl and cyclohexyl substituents are essentially similar.The calculation was based on the empirical value for 2-methoxy-iso-propyl-cyclohexane, ∆ l g  m o (298.15K) = 50.7 ± 0.6 kJ•mol −1 (see Table 2).The "strain-free" reference value 53.9 ± 1.5 kJ•mol −1 for the vaporisation enthalpy of this compound was estimated using the summation enthalpy of Figure 6.Calculating the ortho-interaction between methoxy-group and cyclohexane ring required for calculating the enthalpy of vaporisation of 2-methoxy-bi-cyclohexane using the bi-cyclohexane as the "centerpiece".
The latter contribution was derived as follows.Since thermodynamic data for 2-methoxybi-cyclohexane are completely lacking, we used the 2-methoxy-iso-propyl-cyclohexane to model the required steric interaction, assuming that the spatial characteristics of the iso-propyl and cyclohexyl substituents are essentially similar.The calculation was based on the empirical value for 2-methoxy-iso-propyl-cyclohexane, ∆ g l H o m (298.15K) = 50.7 ± 0.6 kJ•mol −1 (see Table 2).The "strain-free" reference value 53.9 ± 1.5 kJ•mol −1 for the vaporisation enthalpy of this compound was estimated using the summation enthalpy of vaporisation of methoxycyclohexane (Table 2) with the contribution H(C 6 H)→iso-Pr given in Figure S2.The difference ∆ = 50.7 − 53.9 = −3.2kJ•mol −1 can be considered as a reasonable estimate for the decrease in vaporisation enthalpy due to steric repulsions of the substituents in both 2-methoxy-iso-propylcyclohexane and in the 2-methoxy-bi-cyclohexane.With this contribution, the vaporisation enthalpy of 2-methoxy-bi-cyclohexane ∆ g l H o m (298.15K) = 66.1 ± 1.5 kJ•mol −1 (see Table 2) was assessed (see Table 2).The numerical values for the contributions that can be involved in the additive calculations of the hydrogen rich counterparts of the LOHC systems are given in ESI (see Figure S2 and Table S12).
Finally, with the comprehensive evaluation and validation of the ∆  2, the first step is completed.These results can now be combined with the ∆ f H o m (g) values, which are calculated in the following second step.

4.5.
Step II: Gas-Phase Enthalpies of Formation from Quantum Chemical Calculations The most stable conformers for methoxy-biphenyls, methoxy-iso-propylbenzenes, and their perhydrogenated products (see Table S13) were found after conformational analysis performed using a computer code CREST [18].The final optimisation was performed with the B3LYP/6-31g(d,p) method [19].The G4 method was used to calculate the energies E 0 and the enthalpies H 298 of the most stable conformers.The standard molar enthalpies of formation were derived from the resulting H 298 enthalpies using the atomisation (AT) procedure [20] (e.g., for mono-methoxy-biphenyls) as well as using the well-balanced reactions (WBR) shown in Figure 7.
Hydrogen 2023, 4, FOR PEER REVIEW 13 values for the contributions that can be involved in the additive calculations of the hydrogen rich counterparts of the LOHC systems are given in ESI (see Figure S2 and Table S12).
Finally, with the comprehensive evaluation and validation of the ∆ l g  m o (298.15K) values performed in Sections 4.3 and 4.4, the consistent set of vaporisation enthalpies was established.With the results evaluated in Table 2, the first step is completed.These results can now be combined with the ∆  m o (g) values, which are calculated in the following second step.

Step II: Gas-Phase Enthalpies of Formation from Quantum Chemical Calculations
The most stable conformers for methoxy-biphenyls, methoxy-iso-propylbenzenes, and their perhydrogenated products (see Table S13) were found after conformational analysis performed using a computer code CREST [18].The final optimisation was performed with the B3LYP/6-31g(d,p) method [19].The G4 method was used to calculate the energies E0 and the enthalpies H298 of the most stable conformers.The standard molar enthalpies of formation were derived from the resulting H298 enthalpies using the atomisation (AT) procedure [20] (e.g., for mono-methoxy-biphenyls) as well as using the well-balanced reactions (WBR) shown in Figure 7.The theoretical gas-phase enthalpies of formation were calculated according to Hess's Law applied to reactions shown in Figure 7 using reliable experimental gas-phase enthalpies of formation ∆  m o (g, 298.15 K) of reactions participants (see Table S10).The results of the G4 quantum chemical calculations are summarised in Table 6.
Table 6.Comparison of gas-phase enthalpies of formation: calculated using G4, the centrepiece ap- The theoretical gas-phase enthalpies of formation were calculated according to Hess's Law applied to reactions shown in Figure 7 using reliable experimental gas-phase enthalpies of formation ∆ f H o m (g, 298.15 K) of reactions participants (see Table S10).The results of the G4 quantum chemical calculations are summarised in Table 6.Table 6.Comparison of gas-phase enthalpies of formation: calculated using G4, the centrepiece approach, and group additivity at T = 298.15K (p • = 0.1 MPa) for methoxy-substituted biphenyls, iso-propyl-benzenes, and cyclohexanes (in kJ•mol −1 ) a .a The expanded uncertainties of the G4 method were assessed to be ±3.5 kJ•mol −1 [11].b Calculated in this work according to the general atomisation reactions according to Equation ( 8).c Calculated in this work according to well-balanced reactions shown in Figure 7. d Calculated using the "centerpiece" approach as shown in Figure S2.

Compound
e Calculated using group-additivity contributions from Table S12.
As shown in Table 6, the enthalpies of formation calculated using the atomization reaction are systematically more negative than those obtained using the well-balanced reactions.This phenomenon is already known from the literature [20] and, admittedly, the WBR results are preferable for the further thermochemical calculations (see Table 6, column 3).Unfortunately, there are no experimental data on the methoxy-substituted compounds listed in Table 6.Therefore, validation of these theoretical gas-phase enthalpies of formation is required.Such validation can be performed using either the "centrepiece" (CP) approach or conventional group additivity (GA).The contributions required for the application of the CP approach are shown in Figure S2.The contributions required for the application of the GA method are given in Table S12.The results of the empirical calculations using both methods are shown in columns 4 and 5 of Table 6.Admittedly, both empirical methods are not applicable for molecules where the substituents are located in close proximity on the ring (ortho-substituted species in Table 6).However, there is remarkable agreement between the ∆ f H o m (g) WBR and ∆ f H o m (g) CP results for the meta-and para-substituted biphenyls, benzenes, bicyclohexanes, and cyclohexanes (see comparison of columns 3 and 4 in Table 6).The agreement with the GA method is less satisfactory, especially for the series of methoxy-substituted bicyclohexanes (see comparison of columns 3 and 5 in Table 6).The better predictive performance of the CP approach motivates further development and refinement of this method.In any case, both empirical methods adequately support the ∆ f H o m (g) WBR results from quantum chemical calculations.Following, these theoretical values ∆ f H o m (g, 298.15 K) WBR (see Table 6, column 3) can now be safely combined with the vaporisation enthalpies evaluated in step I (see Table 2) to determine the theoretical liquid-phase enthalpies of formation, ∆ f H o m (liq, 298.15 K) theor , according to Equation (1).The latter values are essential for calculations of energetics of hydrogenation and transfer hydrogenation reactions in the third step.

Step III: Energetics of Hydrogenation and Transfer Hydrogenation Reactions Based on Methoxy-Biphenyls
The reaction enthalpies in the liquid phase, ∆ r H o m (liq), are an important parameter for optimizing the conditions for hydrogen storage using LOHCs.They are usually derived from the liquid-phase enthalpies of formation of the reactants according to Hess's Law.The ∆ f H o m (liq, 298.15 K) values for the participants of the reactions considered in Figure 2 are not known.However, with the enthalpies of vaporisation evaluated in Table 2 and the theoretical for the gas-phase enthalpies of formation, ∆ f H o m (g, 298.15 K) theor , the desired ∆ f H o m (liq, 298.15 K) theor values were calculated (see Table 7, column 4) and used for the thermodynamic analysis.a Uncertainties correspond to expanded uncertainties of the mean (0.95 level of confidence).b Calculated in this work according to the well-balanced reactions (see Figure 7).The uncertainties of the G4 method were assessed to be ±3.5 kJ•mol −1 [11].c The evaluated values from Table 2. d The difference in columns 2 and 3. e Calculated via the G4 method using the atomisation reaction.
The enthalpies of reactions hydrogenation/dehydrogenation and transferhydrogenation considered in Figure 2 were derived according to Hess's Law using the standard molar enthalpies of formation of the reactants summarized in Table 7, column 4. The compilation of data required for these calculations and the resulting ∆ r H o m (liq) values are shown in Table 8. a Liquid-phase enthalpies of formation of methoxy-biphenyls (hydrogen-lean (HL) counterparts of the LOHC system) from Table 7. b Liquid-phase enthalpies of formation of perhydrogenated methoxy-biphenyls (or methoxycyclohexanes) from (hydrogen-rich (HR) counterparts of the LOHC system) Table 7. c Calculated according to the Hess's Law applied to the Reaction (1) shown in Figure 2. d Reaction enthalpy per mole H 2 .e From Table S10.
As shown in Table 8, the hydrogenation of the methoxy-biphenyls is highly exothermic from −371 to −402 kJ•mol −1 (see Table 8, column 4), and this energy must be taken into account when optimising the design of chemical reactors.It is convenient to relate the enthalpy of reaction to the amount of hydrogen released (kJ•mol −1 / H2 ).This allows a comparison of the reaction enthalpies of LOHC systems with different stoichiometries (see Table 8, column 5).In these units the reaction enthalpies of the methoxy-biphenyls (at −61.7 to −67.1 kJ•mol −1 / H2 ) are very similar to those of methoxy-benzene (−65.6 kJ•mol −1 / H2 ) or commercially available thermofluids benzyltoluenes and dibenzyltoluenes (−63.5 and −65.4 kJ•mol −1 / H2 ), respectively [21]).
The enthalpies of TH reactions according to Reactions 2 (see Figure 2) were derived using the standard molar enthalpies of formation of the reactants summarized in in Ta-ble 7, column 4. The compilation of data required for these calculations and the resulting ∆ r H o m (liq) values are shown in Table 9. a Liquid-phase enthalpy of formation of acetone from Table S10.b Liquid-phase enthalpies of formation of iso-propanol from Table S10.c Liquid-phase enthalpies of formation of perhydrogenated methoxy-biphenyls (or methoxy-cyclohexanes) from (hydrogen-rich (HR) counterparts of the LOHC system) Table 7. d Liquid-phase enthalpies of formation of methoxy-biphenyls (hydrogen-lean (HL) counterparts of the LOHC system) from Table 7. e Calculated according to the Hess's Law applied to the Reactions 2 in Figure 2, using liquid-phase enthalpies of formation of acetone and iso-propyl alcohol from Table S10.
As it shown in Table 9, column 6, the TH reactions based on the methoxy-cyclohexanes are significantly less exothermic compared to direct hydrogenation.Such an energetic advantage also speaks to the mild reaction conditions in TH compared to conventional dehydrogenation/hydrogenation processes.Therefore, our ongoing studies on TH reactions with different types of hydrogen acceptors and donors could open up new possibilities for hydrogen storage and utilisation in the future, especially if the donors and acceptors could come from renewable sources.

Thermodynamic Analysis of the Reversible Hydrogenation/Dehydrogenation Process in the Gas Phase
Admittedly, the dehydrogenation of the hydrogen-rich LOHC counterpart is a thermodynamically unfavourable process [22].Therefore, proper selection of optimal dehydrogenation temperatures is a direct way to increase hydrogen yield.In our recent work [23], we proposed an approach to evaluate the optimal dehydrogenation temperatures based on quantum chemical calculations of the Gibbs free energy.According to the Gibbs free energy equation, at low temperatures the enthalpy term predominates over the entropy term, so that hydrogenation is possible.When ∆ r G o m = 0, the enthalpy term ∆ r H o m and the entropy term T × ∆ r S o m in Equation ( 9) become equal.This occurs at a certain "equilibrium temperature", T eq , which can imply the reversal of thermodynamic feasibility from hydrogenation to dehydrogenation tendency: This temperature should be a valuable indicator for the selection and development of an effective LOHC system.This T eq should not be too low for the hydrogenation step to achieve acceptable reaction rates; the T eq for the dehydrogenation step should not be too high but sufficient to achieve good conversion and selectivity.Results for T eq calculations according to Equation (10) for methoxy-cyclohexanes are compared in Table 10.
As can be seen in Table 10, the T eq values for dehydrogenation in methoxy-substituted LOHC systems vary between 516 and 539 K, depending on the position of the CH 3 O group on the aliphatic rings.The lowest equilibrium temperature is observed for 4,4'-dimethoxybicyclohexane (516 K).However, except for 2-methoxy-bicyclohexane (527 K), the T eq values for all methoxy-substituted species are very close (see Table 10).
Table 10.Thermodynamic analysis of the dehydrogenation of methoxy-cyclohexanes in the gas phase calculated using the G4 method (at 298.15K and p • = 0.1 MPa).It is interesting to note that the dehydrogenation of 2-methoxy-bicyclohexane occurs with the highest T eq among the methoxy-compounds considered in Table 10, and at the same time the dehydrogenation enthalpy for this compound (of 67.1 kJ•mol −1 / H2 , see Table 10) is also the highest.This observation is important for the proper selection of LOHC systems.Obviously, the ortho-substituted methoxy-substituted aromatic compounds are less suitable for hydrogen storage via reversible hydrogenation/dehydrogenation.In contrast, hydrogen storage via transfer hydrogenation with 2-methoxybicyclohexane is more advantageous by a factor of 1.5 compared to the meta-and para-isomers, since the reaction enthalpy of −17.6 kJ•mol −1 is the lowest in this series (see Table 9).These new findings contribute to a better understanding of the reaction thermodynamic properties of LOHC systems and reveal the structural features of aromatic compounds that are essential for their technical applications.

Conclusions
Methoxy-substituted biphenyls have been studied thermodynamically regarding their applicability as LOHC materials.The consistent sets of phase transition enthalpies and formation enthalpies were obtained using vapour pressure determinations, thermal behaviour studies, empirical, and quantum chemical calculations.These results were recommended for chemical engineering calculations to optimize reversible reactions for the storage and release of hydrogen.The reaction enthalpies of hydrogenation of methoxybiphenyls in the liquid phase were calculated and compared with other types of hydrogen carriers.Generally, a reduction in enthalpy of reaction for the dehydrogenation reaction is desired to achieve more beneficial reaction conditions.The hydrogenation enthalpies were found to be of the same order of magnitude as for the common aromatic compounds.The results show that substitution with methoxy groups on the ring can have a positive effect on the reaction thermodynamics.However, the effect is dependent on the position of the substitution.The meta-and para-substitution of the aromatic ring increases the positive effect of the introduction of the methoxy group.The "equilibrium temperatures" of the meta-and para-methoxy-substituted compounds are significantly lower compared to those of the ortho-substituted isomers.The energetics of the TH reaction based on the methoxycyclohexanes was also calculated and proved to be significantly less exothermic compared to direct hydrogenation.The results of this study may contribute to a better understanding of the reaction thermodynamics of potential LOHC materials and their chemical structure.However, from a technical point of view, it is important to consider that the appending of an ether functionality to a LOHC compound introduces some risk of ether cleavage during the catalytic hydrogenation step that loads the LOHC compound with hydrogen.This issue needs to be addressed in catalyst development studies.The thermodynamic analysis presented here shows that such an approach could be very rewarding to provide promising LOHC systems for efficient hydrogen release in mild temperature conditions.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/hydrogen4040052/s1,Table S1 S13: Structures of the most stable conformers and the G4 calculated gas-phase enthalpies of formation ∆ f H o m (g) G4 at T = 298.15K (p • = 0.1 MPa) for methoxybicyclohexyls (in kJ•mol −1 ); Figure S1: Correlation of vaporisation enthalpies, ∆ g l H o m (298.15K), of hydroxy-iso-propylbenzenes and methoxy-iso-propylbenzenes; Figure S2: The development of the contributions for calculating the vaporisation enthalpies (or gas-phase enthalpies of formation) of aromatic and cyclic aliphatic hydrocarbons using the "centerpiece" approach.The auxiliary experimental and computational details necessary to understand the main points of the paper are given in Refs. .

(Figure 2 .
Figure 2. Chemical reactions of the reversible hydrogen storage using methoxy-biphenyls as the liquid organic hydrogen carriers (LOHCs).

(Figure 2 .
Figure 2. Chemical reactions of the reversible hydrogen storage using methoxy-biphenyls as the liquid organic hydrogen carriers (LOHCs).

Figure 2 .
Figure 2. Chemical reactions of the reversible hydrogen storage using methoxy-biphenyls as the liquid organic hydrogen carriers (LOHCs).
g l H o m (298.15K), of hydroxy-biphenyls and methoxy-biphenyls exhibits a practically straight line (see Figure 3, left), which can be taken as an indicator of the internal consistency of the ∆ g l H o m (298.15K) values in both series.Hydrogen 2023, 4, FOR PEER REVIEW 9

Figure 3 .
Figure 3. Correlation of vaporisation enthalpies, ∆ l g  m o (298.15K), of hydroxy-biphenyls and methoxy-biphenyls.All values are in kJ•mol −1 with R = position of the substituent on the benzene ring.

Figure 3 .
Figure 3. Correlation of vaporisation enthalpies, ∆ g l H o m (298.15K), of hydroxy-biphenyls and methoxybiphenyls.All values are in kJ•mol −1 with R = position of the substituent on the benzene ring.

Figure 3 .
Figure 3. Correlation of vaporisation enthalpies, ∆ l g  m o (298.15K), of hydroxy-biphenyls and methoxy-biphenyls.All values are in kJ•mol −1 with R = position of the substituent on the benzene ring.

Figure 3 .
Figure 3. Correlation of vaporisation enthalpies, ∆ l g  m o (298.15K), of hydroxy-biphenyls and methoxy-biphenyls.All values are in kJ•mol −1 with R = position of the substituent on the benzene ring.

Figure 4 .
Figure 4. Comparison of the quantitative differences in ∆ lg  m o (298.15K)-values: between para-and ortho-methoxy-biphenyls and between para-and ortho-methoxy-iso-propylbenzenes (left) and the corresponding differences between para-and ortho-hydroxy-biphenyls and between para-and ortho-hydroxy-iso-propylbenzenes (right).Numerical data are given in TablesS8 and S10.All values are in kJ•mol −1 .

Figure 4 .
Figure 4. Comparison of the quantitative differences in ∆ g l H om (298.15K)-values: between para-and ortho-methoxy-biphenyls and between para-and ortho-methoxy-iso-propylbenzenes (left) and the corresponding differences between para-and ortho-hydroxy-biphenyls and between para-and orthohydroxy-iso-propylbenzenes (right).Numerical data are given in TablesS8 and S10.All values are in kJ•mol −1 .

Figure 5 .
Figure 5. Calculating the enthalpy of vaporisation of methoxy-bi-cyclohexanes using the bicyclohexane as the "centerpiece".
values performed in Sections 4.3 and 4.4, the consistent set of vaporisation enthalpies was established.With the results evaluated in Table
15 K) is the standard molar enthalpy of vaporisation.In our experience, ∆ g l H o m (298.15K) values can be relatively easily measured using various methods or estimated using various group additivity rules and structure-property correlations.Under the latter prerequisite, a reasonable roadmap can be proposed to comprehensively obtain the thermodynamic quantities essential for reversible hydrogen storage.
-Step I: the possible experimental and empirical methods are involved to obtain and validate the ∆ g l H o m (298.15K) values; -Step II: the high-level quantum chemical calculations are performed to compute ∆ f H o m (g, 298.15 K) values and the possible experimental and empirical methods are involved to validate these results; -Step III: the liquid-phase ∆ f H o m (liq, 298.15 K) values are derived according to Equation (

Table 1 .
Vapour pressures p and standard molar functions of sublimation/vaporisation measured in this work using the transpiration method.
The ∆ l cr H o m at T fus values were adjusted to 298.15 K as shown in ESI.Uncertainties of ∆ l cr H o m (298 K) were estimated with 30% of the total adjustment.c Enthalpies of vaporisation evaluated in Table 2. d Sum of columns 4 and 5 in this table.
a Uncertainties are presented as expanded uncertainties (0.95 level of confidence with k = 2).b

Table S7 :
: Provenance and purity of the materials; TableS2: Compilation of data on molar heat capacities C o p,m (liq) and heat capacity differences ∆ The vapour pressures p, standard molar vaporisation enthalpies, and entropies obtained using the approximation of boiling points at different pressures available in the literature; TableS4.Correlation of vaporisation enthalpies, ∆ g l H o m (298.15K), of alkoxy-substituted aromatics with their Kovats indices (J x ); TableS5.Correlation of vaporisation enthalpies, ∆ g l H o m (298.15K), of alkoxy-substituted aromatics with their Kovats indices (J x ); TableS6: Correlation of vaporisation enthalpies, ∆ Correlation of the normal boiling temperatures, T b , of dialkyl ethers with their J x (retention indices); TableS8: Enthalpies of fusion (∆ l melting points (T fus ) and 298 K (in kJ•mol −1 ); TableS9: Correlation of vaporisation enthalpies, ∆ g l H o m (298.15K), of substituted phenols with their Kovats indices (J x ); TableS10: Thermochemical data at T = 298.15K(p• = 0.1 MPa) for reference compounds (in kJ•mol −1 ); TableS11: Correlation of the normal boiling temperatures, T b , of dialkyl ethers with the standard molar enthalpies of vaporisation, ∆ g l H o m (298.15K); Table S12: Group-additivity values Γ i for calculation of enthalpies of vaporisation, ∆ g l H o m , and enthalpies of formation, ∆ f H o m (g), of alkanes, ethers, and aromatics at 298.15 K in kJ mol −1 ; Table g l C o p,m (in J.K −1 •mol −1) at T = 298.15K;TableS3: