Next Article in Journal
Optimizing Human–AI Collaboration in Chemistry: A Case Study on Enhancing Generative AI Responses through Prompt Engineering
Previous Article in Journal
Efficient Two-Step Synthesis of Novel Pyrimido[4,5-d] Pyrimidines with Potent Neuroprotective, Antioxidant, and Aβ Anti-Aggregation Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermodynamic Exercises for the Kinetically Controlled Hydrogenation of Carvone

by
Artemiy A. Samarov
1,
Sergey V. Vostrikov
2,
Aleksandr P. Glotov
3 and
Sergey P. Verevkin
4,5,*
1
Department of Chemical Thermodynamics and Kinetics, Saint Petersburg State University, Peterhof, 198504 Saint Petersburg, Russia
2
Chemical-Technological Department, Samara State Technical University, 443100 Samara, Russia
3
Department of Physical and Colloid Chemistry, Gubkin Russian State University of Oil and Gas, 65 Leninsky Prospekt, 119991 Moscow, Russia
4
Competence Centre CALOR at the Department Life, Light & Matter of Faculty of Interdisciplinary Research at University of Rostock, 18059 Rostock, Germany
5
Department of Physical Chemistry, Kazan Federal University, 420008 Kazan, Russia
*
Author to whom correspondence should be addressed.
Chemistry 2024, 6(4), 706-722; https://doi.org/10.3390/chemistry6040042
Submission received: 26 June 2024 / Revised: 3 August 2024 / Accepted: 8 August 2024 / Published: 10 August 2024
(This article belongs to the Section Theoretical and Computational Chemistry)

Abstract

:
Carvone belongs to the chemical family of terpenoids and is the main component of various plant oils. Carvone and its hydrogenated products are used in the flavouring and food industries. A quantitative thermodynamic analysis of the general network of carvone hydrogenation reactions was performed based on the thermochemical properties of the starting carvone and all possible intermediates and end products. The enthalpies of vaporisation, enthalpies of formation, entropies and heat capacities of the reactants were determined by complementary measurements and a combination of empirical, theoretical and quantum chemical methods. The energetics and entropy change in the hydrogenation and isomerisation reactions that take place during the conversion of carvone were derived, and the Gibbs energies of the reactions were estimated. It was shown that negative Gibbs energies are recorded for all reactions that may occur during the hydrogenation of carvone, although these differ significantly in magnitude. This means that all these reactions are thermodynamically feasible in a wide range from ambient temperature to elevated temperatures. Therefore, all these reactions definitely take place under kinetic and not thermodynamic control. Nevertheless, the numerical Gibbs energy values can help to establish the chemoselectivity of catalysts used to convert carvone to either carvacarol or to dihydro- and terahydrocarvone, either in carvotanacetone or carveol.

Graphical Abstract

1. Introduction

Renewable raw materials have received increasing attention as alternatives to petrochemical raw materials with a view to modern sustainable industrial development. The fragrance and flavour industry is based on terpenes. Carvone belongs to this chemical family and occurs naturally as the main component of caraway, dill, gingergrass, and spearmint oils. The various applications of carvone in the flavour and food industries as a starting material for pharmaceutical compounds, antimicrobial agents, etc. justify research into increasing the global production of this monoterpene [1]. Carvone is obtained at low cost by the steam distillation of these oils or recovered from citrus-derived limonene [2]. The production of carvone is also possible with specialised microorganisms using biotechnological techniques; e.g., the biosynthesis of (R)-(˗)-carvone from glucose was possible with the help of genetically modified Escherichia coli [3].
The upgrading of biomass can improve the sustainability of processes, but it is still at an early stage of development. For example, the catalytic hydrogenation of carvone produces valuable chemicals (Figure 1 and Figure 2) for the food, pharmaceutical, and agricultural industries.
However, the general good selectivity of these reactions is challenging, and most of the catalysts used are not selective, leading to product mixtures [4,5]. To address this problem, a thermodynamic analysis of the reactions shown in Figure 1 and Figure 2 is required. Carvone is an interesting substrate for the study of thermodynamic stability and catalytic activity in stereo- and chemoselective hydrogenation because it contains an asymmetric centre and three types of double bonds that can be hydrogenated: carbonyl C=O, a conjugated endo-double bond C=C (incorporated within the six-membered ring) and an isolated exo-double bond C=C outside the six-membered ring.
As shown in Figure 1, there are four possible directions in the hydrogenation of carvone (reactions 1,4,6) and isomerisation reaction 12. Each of these directions determines a cascade of subsequent reactions that lead to products of different chemical classes. The focus of this study is to quantify the energetics of the reactions shown in Figure 1 and Figure 2 to understand which of the hydrogenation directions is preferable from a thermodynamic perspective. These preferences are usually helpful in selecting specific catalysts that will drive the process in the desired direction.
The thermodynamic analysis of chemical reactions usually involves a complex of basic thermodynamic properties such as enthalpies, entropies, heat capacities, etc. of the reactants shown in Figure 1 and Figure 2. The thermochemical experiments to measure these properties require large quantities of high-purity (99.9%) samples. Since the carvone derivatives have very similar physico-chemical properties, the separation and purification of the individual compounds is a very challenging and demanding task. In order to possibly reduce the experimental effort, a hybrid thermodynamic method for evaluating the feasibility of chemical reactions is outlined in this paper. This hybrid approach combines the experimental and theoretical methods in an appropriate balance that leads to a realistic thermodynamic analysis of the feasibility and energetics of the network of chemical reactions depicted in Figure 1 and Figure 2. The thermodynamic workflow of this hybrid approach is as follows.

2. Thermodynamic Workflow

A chemical reaction is generally feasible if the sign of the standard Gibbs energy of the reaction, r G m o , is negative (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
where r H m o is the standard molar reaction enthalpy, r S m o is the change in the standard molar entropy of a chemical reaction and K p is the gas-phase thermodynamic equilibrium constant. If the reaction takes place in the liquid phase, K p in Equation (1) is replaced by K a , which represents the liquid-phase thermodynamic equilibrium constant. The latter Equation (1) is known in textbooks as the Gibbs–Helmholtz equation. The K p and K a values are directly related to the yield of the desired reaction products. The equilibrium constants that are ≥1 promise the already acceptable yields at a particular temperature. The negative and large r G m o values are an indication of a large thermodynamic driving force for the chemical reactions of interest. As a consequence, the reaction to the product is favorable and proceeds under thermodynamic control. However, this does not exclude the need to search for a suitable catalyst and temperature, as kinetics is also important to achieve high reaction rates in order to reach equilibrium.
There are two contributions to the Gibbs–Helmholtz equation, the enthalpic, r H m o , and entropic, r S m o . In the first contribution, the standard molar reaction enthalpy can be calculated for each chemical reaction according to Hess’s Law:
r H m o ( T ) = Σ [ f H m o ( T ,   products ) ] Σ [ f H m o ( T ,   educts ) ]
using the standard molar enthalpies of formation, f H m o (T), of reactants. Most chemical reactions in industry take place either in the gas or liquid phase, so the corresponding enthalpies of formation must be used in Equation (2). These values are linked using the following equation:
f H m o ( liq ,   T ) = f H m o ( g ,   T ) l g H m o ( T )
where f H m o (liq, T) is the liquid-phase enthalpy of formation, f H m o (g, T) is the gas-phase enthalpy of formation and l g H m o (T) is the standard molar enthalpy of vaporisation. The standard molar enthalpies of formation and enthalpies of vaporisation can either be of experimental origin or obtained using various empirical methods, e.g., from group-additivity rules or structure–property correlations.
The second contribution to the Gibbs–Helmholtz equation, the change in the standard molar entropy of the reaction, is derived as follows from the standard molar entropies, S m o (T), of the reactants:
r S m o ( T ) = Σ [ S m o ( T ,   products ) ] Σ [ S m o ( T ,   educts ) ]
The liquid and gas phase entropies are linked using the following equation:
S m o ( liq ,   T ) = S m o ( g ,   T ) l g S m o ( T )
where S m o (liq, T) is the liquid-phase standard molar entropy at temperature T, which can be measured using cryogenic adiabatic calorimetry or estimated using empirical methods; S m o (g, T) is the gas-phase entropy, which is usually calculated using statistical thermodynamics or a suitable quantum chemical (QC) method; and l g S m o (T) is the standard molar entropy of vaporisation, which can be derived from the slope of the temperature dependence of the experimental vapour pressures as follows:
l g S m o T = l g H m o / T + R × l n p i / p o
where p i is the absolute vapour pressure of a pure compound i at a given temperature T and p o = 0.1 MPa is the standard pressure.
It is common practice to carry out thermodynamic calculations at the reference temperature T = 298 K. However, the typical temperature range for carvone hydrogenation reactions is between 293 and 523 K [4,5]. The variation in the reaction enthalpies and entropies with temperature is given by the thermodynamic relations:
r H m o ( T ) = r H m o ( 298   K ) + r C p , m o × ( T 298   K )
r S m o ( T ) = r S m o ( 298   K ) + r C p , m o × ln ( T / 298   K )
where r C p , m o is the change in standard molar heat capacity for the reaction, which is calculated from the standard molar heat capacities, C p , m o (liq or g), of educts and products. The latter values can be easily measured using differential scanning calorimetry (DSC) or estimated using empirical methods.
The final input properties of the thermodynamic workflow according to Equations (2)–(6) are f H m o (liq, 298 K), l g H m o (298 K), and S m o (liq), for all the reactants shown in Figure 1 and Figure 2. Since most of these compounds are not commercially available and their synthesis and purification is not trivial, only a reasonable combination of experimental, empirical and quantum chemical methods can contribute to a reliable thermodynamic analysis and optimisation of hydrogenation reactions. The following steps leading to the solution of the Gibbs–Helmholtz equation for each of the fifteen reactions of the carvone hydrogenation reactions were performed in this work.
-
Step I: the high-level QC calculations are first performed in order to obtain f H m o (g) and S m o (g) values.
-
Step II: vapour pressure–temperature dependencies for the molecules involved in the reaction network are measured or taken from the literature. All possible experimental and empirical data and methods are used to evaluate and validate the l g H m o and l g S m o values.
-
Step III: the f H m o (liq) and S m o (liq) values in the liquid phase are derived from the results of the first and second steps. The enthalpies of formation of the liquid phase are validated (e.g., with the help of combustion experiments in this work). The r H m o (liq) and r S m o (liq) values are derived according to Hess’s Law and inserted into Equation (1) to derive the r G m o (liq) of a desired reaction, and an analysis of all three quantities is performed.
This workflow helps to follow the combination of experimental, empirical and quantum chemical methods used in this work to evaluate and compare the feasibility and energetics of the reaction network depicted in Figure 1 and Figure 2.

3. Theoretical and Experimental Methods

The Gaussian 16 series software [6] was used for quantum chemical calculations under the “rigid rotator-harmonic oscillator” assumption. The total H298 enthalpies of the most stable conformers were calculated using the G4 method [7]. The H298 values were finally converted to the theoretical f H m o (g, 298 K) values and discussed. The S m o (g, 298 K)QC values were calculated according to (Equation (1)) using H298 and G298 from the output file.
The transpiration method was used to measure the absolute vapour pressures of (-)-carvone. The vapour pressures temperature dependence was used to derive the enthalpy and entropy of vaporisation of (-)-carvone. The necessary details can be found in ESI.
The standard molar energy of combustion of (-)-carvone was measured with a home-made isoperibolic calorimeter with a static bomb and a stirred water bath. The standard molar enthalpy of formation of (-)-carvone in the liquid phase, f H m o (liq), was obtained from the specific energy of combustion. The sample of (-)-carvone used in this work was of commercial origin (see Table S1).

4. Results and Discussion

4.1. Step I: Theoretical Gas-Phase Enthalpies of Formation from Quantum Chemistry

The most stable conformer for the reactants shown in Figure 1 and Figure 2 was localised using a computer code CREST (conformer-rotamer ensemble sampling tool) [8] and optimised with the B3LYP/6-31g(d,p) method [9]. The structures of the most stable conformers of carvone are shown in Figure 3 as an example.
The structures of the most stable conformers of carvone derivatives are shown in Figure S1. The general atomisation (AT) reaction
CmHnOk = m × C + n × H + k × O
was used to convert the total H298 values available directly from the output file to the standard molar enthalpies of formation f H m o (g, 298 K)AT. However, the AT method must be corrected using a linear correlation between the f H m o (g, 298 K)AT values and the experimental enthalpies of formation, f H m o (g, 298 K)exp, of a number of molecules containing similar functional groups (as the reactants in Figure 1 and Figure 2) with reliable experimental data (see Table 1, column 3).
We deliberately included the aliphatic and aromatic alcohols (phenols) as well as the cyclic aliphatic alcohols and ketones in the correlation. However, all these structures lie on the same line with a very high correlation coefficient R2 = 0.9997 as follows:
f H m o ( g ,   G 4 ) corr / kJ · mol 1 = 0.9818 × f H m o ( g ,   AT ) 4.3   with   R 2 = 0.9997
This equation was used to obtain the “corrected” G4 results for the carvone derivatives (see Figure 1 and Figure 2). The summary of the QC calculations is given in Table 2.
Table 2. Compilation of the G4 calculated gas-phase enthalpies of formation for the carvone derivatives (at T = 298 K, = 0.1 MPa, in kJ·mol−1).
Table 2. Compilation of the G4 calculated gas-phase enthalpies of formation for the carvone derivatives (at T = 298 K, = 0.1 MPa, in kJ·mol−1).
CASCompound f H m o (g, AT) a f H m o (g, G4)corr b l g H m o c f H m o (liq)theor d
6485-40-1(-)-carvone−122.1−124.261.7−185.9
99-48-9trans-carveol−168.7−169.972.6−242.5
5524-05-0trans-(±)-dihydrocarvone−221.1−221.457.2−278.6
3792-53-8cis-(±)-dihydrocarvone−220.0−220.356.8−277.1
499-71-8carvotanacetone−238.4−238.459.3−297.7
499-70-7(±)-carvomenthone−337.1−335.358−393.3
cis-(±)-tetrahydrocarvone−335.1−333.358.4−391.7
499-75-2carvacrol−210.8−211.372.5−283.8
499-69-4carvomenthol−397.0−394.169.9−464.0
619-01-2trans-dihydrocarveol−281.8−281.071.9−352.9
22567-21-1cis-dihydrocarveol−274.0−273.370.1−343.4
536-30-1carvotanacetol−280.5−279.773.1−352.8
a Calculated using the atomisation reactions Equation (9) with the expanded uncertainties of ±3.5 kJ·mol−1 [7]. b The result obtained from the atomisation reaction Equation (9) was corrected with Equation (10). c From Table 3. d Difference columns 4 and 5 in this table.
Table 3. Compilation of the standard molar enthalpies of vaporisation, l g H m o , of the carvone derivatives (in kJ·mol−1).
Table 3. Compilation of the standard molar enthalpies of vaporisation, l g H m o , of the carvone derivatives (in kJ·mol−1).
Compounds, CAST-Range/
K
l g H m o /
Tav
l g H m o / a
298 K
(˗)-carvone, 6485-40-1279.6–333.061.4 ± 0.261.7 ± 0.3 b
320–50452.2 ± 0.560.2 ± 1.7
cis-(±)-dihydrocarvone, 3792-53-8322–49549.7 ± 0.856.8 ± 1.6
trans-(±)-dihydrocarvone, 5524-05-0343–49549.4 ± 0.657.2 ± 1.7
cis-(±)-tetrahydrocarvone335–49650.6 ± 0.558.4 ± 1.6
trans-(±)-tetrahydrocarvone (carvomenthone), 499-70-7335–49450.0 ± 0.558.0 ± 1.7
carvotanacetone, 499-71-8317–50152.4 ± 0.759.3 ± 1.5
carvomenthol, 499-69-4354–49153.7 ± 0.569.9 ± 3.3
(±)-cis-carveol, 1197-06-4338–50457.0 ± 1.271.0 ± 3.0
(±)-trans-carveol, 99-48-9314–50158.5 ± 0.872.6 ± 2.9
trans-dihydrocarveol, 619-01-2350–49855.0 ± 0.671.9 ± 3.4
cis-dihydrocarveol, 22567-21-1353–49855.4 ± 1.370.1 ± 3.6
carvotanacetol, 536-30-1337–49559.6 ± 1.373.1 ± 3.0
carvacrol, 499-75-2 72.5 ± 0.1 c
2-methyl-2-cyclohexen-1-one, 1121-18-2325–45244.3 ± 0.348.2 ± 0.8 d
2-methylcyclohex-2-en-1-ol, 20461-30-7298–44154.0 ± 1.060.2 ± 1.3 d
a Vapour pressures from Table S2 were treated using Equations (11) and (12) using the heat capacity differences from Table S4 to calculate the enthalpies of vaporisation at 298 K. Uncertainties of the 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 K. Uncertainties in the temperature adjustment of vaporisation enthalpies to the reference temperature T = 298 K are estimated to account with 20% to the total adjustment. b From Table S3. c From Ref. [13]. d Used as model compounds for the calculations according to the “centrepiece” approach (see Section 4.2.3).
The “corrected” f H m o (g, G4)corr values for carvone derivatives (Table 2, column 4) obtained in this step are now ready for thermochemical calculations in further steps.

4.2. Step II: Vaporisation Thermodynamics

4.2.1. Absolute Vapour Pressures

The vapour pressures of carvone derivatives are of interest for the appropriate design of the technology, the separation of the products and for safety reasons. Unfortunately, the systematic vapour pressure–temperature dependencies for carvone derivatives are practically non-existent in the literature (except for carvacrol [13]). However, the boiling points (BPs) of these compounds at various reduced pressures were collected from the literature (see a compilation in Table S2), and these data were evaluated to obtain vapour pressure–temperature dependencies for each compound and derive the corresponding enthalpies of vaporisation. It is important to note that these boiling points were recorded during the distillation of the reaction mixtures after synthesis. Therefore, the temperatures were given in the range of a few degrees, and the pressures were measured with uncalibrated manometers. In our experience, even such rough data generally show a reasonable trend of vapour pressure–temperature dependencies [14]. To prove this general observation for the carvone derivatives, we collected the available BP for (-)-carvone in the range 320–504 K (see Table S2), and to verify these data, the vapour pressure–temperature dependence in the range 279.6–333.0 K for (-)-carvone was measured by the transpiration method (see Table S3).
In Figure S2, the experimental vapour pressures for (-)-carvone were compared with the boiling points found in the literature for this compound at reduced pressures (see Table S2). As can be seen from this figure, the results from both data sets generally show the same trend. Once the applicability of the BP method for carvone derivatives was established, the boiling points at reduced pressure for each compound were collected (see Table S2) and involved in the thermodynamic analysis of the carvone hydrogenation reaction network.
The vapour pressures collected in Tables S2 and S3 were fit by the following equation:
R · ln ( p i / p r e f ) = a + b T + l g C p , m o · l n T T 0
where a and b are adjustable parameters, T0 applied in Equation (11) is the arbitrary temperature (was chosen to be T0 = 298 K in this work), R = 8.314462 J·K−1·mol−1 is the molar gas constant, and the reference pressure p r e f = 1 Pa ,   c r , l g C p , m o = C p , m o (g) − C p , m o (liq) is the difference between the molar heat capacities of the gaseous C p , m o (g) and the liquid phase C p , m o (liq), respectively. The l g C p , m o values used in Equation (11) are given in Table S4. The vapour pressures and fitting parameters of Equation (11) are given in Tables S2 and S3.

4.2.2. Vaporisation Enthalpies

The standard molar enthalpies of vaporisation of carvone derivatives at temperatures T were derived from the temperature dependence of the vapour pressures, which was approximated by Equation (11) using the following equation:
l g H m o T = b + l g C p , m o × T
where b is one of the adjustable parameters of Equation (11). The vaporisation enthalpies at the reference temperature T = 298 K, l g H m o (298 K), of (carvone derivatives are summarised in Table 3.
As can be seen from this table, for (-)-carvone, our transpiration result l g H m o (298 K) = (61.7 ± 0.3) kJ·mol−1 agrees very well with the vaporisation enthalpy, l g H m o (298 K) = (60.2 ± 1.7) kJ·mol−1 derived in this work from the boiling points collected for this compound in Table S2.

4.2.3. Validation of Enthalpies of Vaporisation Using the “Centrepiece” Approach

Before the enthalpies of vaporisation of the carvone derivatives summarised in Table 3 are used for further thermochemical calculations, it is useful to validate these results with empirical methods based on structure–property correlations. The crucial advantage of empirical methods is that they are developed from reliable experimental data and have internal consistency. Therefore, the new results should be integrated into the network of already known data. If the new results are consistent with the already known trends, they can be validated for the thermochemical calculations. If this is not the case, the reason for the specific behaviour should be uncovered and justified.
A group addition (GA) method is one of the best manifestation of structure–property correlation. The conventional methodology of GA methods is to divide the experimental enthalpies of the molecules into the smallest possible groups and calculate numerical contributions for them from the matrix of well-established data. The prediction is then made by building a framework of the desired model molecule from the corresponding number and type of these contributions. Admittedly, the conventional GA methods have difficulties with the cyclic molecules. To overcome these difficulties, a series of unique “ring” corrections are proposed to improve the prediction. To address this shortcoming, the GA method can be modified by taking the specific ring as the centrepiece as the starting point for constructing the frame of the desired model molecule when the small additional contributions are attached to the “centrepiece”.
The general idea of this approach is to first select a potentially large “centrepiece” molecule that has a reliable enthalpy and can generally mimic the structure of other molecules of interest [15]. Then, the required substituents (or groups) are added to the “centrepiece” molecule, leading to the construction of the desired molecule. In this work, cyclohexane derivatives (cyclohexanone and cyclohexanol derivatives) can serve as a suitable “centrepiece” for the prediction of the enthalpy of vaporisation of carvone derivatives. Two types of substituents are specific to the carvone derivatives: the iso-propenyl and the iso-propyl substituents for the exchange of the H atom in the cyclohexane ring, as shown in Figure 4.
The idea of the “centrepiece” approach for carvone derivatives is illustrated in Figure 5.
To construct the carvone, the contribution for the substituent H→iPr(en) = 10.8 kJ·mol−1 (see Figure 4) should be added to the enthalpy of vaporisation of 2-methyl-2-cyclohexen-1-one l g H m o (298 K) = 48.2 ± 0.8 kJ·mol−1 (see Table S2), which is used as the “centrepiece”. The resulting value l g H m o (298 K) = 59.0 ± 1.4 kJ·mol−1 (see Figure 5) has acceptable agreement with the value l g H m o (298 K) = 61.7 ± 0.3 kJ·mol−1 (see Table 3) measured with transpiration, as well as with the value l g H m o (298 K) = 60.2 ± 1.7 kJ·mol−1 (see Table S2) derived from BP in this work.
To construct the trans-tetrahydrocarvone, the contribution for the substituent H→iPr = 10.9 kJ·mol−1 (see Figure 4) should be added to the enthalpy of vaporisation of 2-methyl-cyclohexanone l g H m o (298 K) = 46.5 ± 0.2 kJ·mol−1 (see Table S5), which is used as the “centrepiece”. The resulting value l g H m o (298 K) = 57.4 ± 1.0 kJ·mol−1 (see Figure 5) is in agreement with the value l g H m o (298 K) = 58.0 ± 1.7 kJ·mol−1 (see Table 4) derived from BP in this work.
Similarly, the enthalpies of vaporisation for other carvone derivatives were calculated, as shown in Figures S3–S5. The “empirical” enthalpies of vaporisation derived using the “centrepiece” approach are compared in Table S6 and agree well with results obtained by other methods. This validation using the structure–property correlations makes it possible to recommend the l g H m o (298 K) values (Table 3) for further thermochemical calculations.

4.2.4. Entropies of Vaporisation

The change in the liquid-phase standard molar entropy of a chemical reaction, r S m o (liq), as the entropic contribution to the Gibbs–Helmholtz equation is required for further thermochemical calculations. The standard molar entropies, S m o (liq, 298 K), of reactants (see Figure 1 and Figure 2) in the liquid phase are necessary as input for the calculation of the r S m o (liq, 298 K) values according to Equation (4). The standard molar entropies, S m o (liq, 298 K), of the carvone derivatives were obtained using Equation (5). The gas-phase entropies were calculated using the G4 method. The standard molar entropies of vaporisation, l g S m o (298 K), were calculated according to Equation (6). The compilation of the standard molar entropies of vaporisation, l g S m o (298 K), and the absolute standard molar entropies, S m o (g or liq, 298 K), of carvone derivatives is shown in Table 4.
The resulting liquid-phase standard molar entropies, S m o (liq), required to estimate the entropic contribution to the Gibbs–Helmholtz equation (see Equation (1)) are given in Table 4, column 5.

4.3. Step III: Thermodynamic Analysis of the Reaction Network for the Hydrogenation of Carvone

4.3.1. Liquid-Phase Enthalpies of Formation and Liquid-Phase Reaction Enthalpies

The “corrected” quantum–chemical gas-phase enthalpies of formation, f H m o (g, 298 K)G4, of the carvone derivatives were calculated using the G4 method (Table 2, column 4). They were recalculated to the liquid-phase enthalpies of formation, f H m o (liq, 298 K), via the enthalpies of vaporisation, l g H m o (298 K), which are shown in Table 3. The f H m o (liq, 298 K) values were calculated according to Equation (3), and the results are listed in Table 2 (last column).
The enthalpies of formation of the carvone derivatives in the liquid phase determined in this way must also be validated before they can be used for the thermodynamic analysis of the reaction network shown in Figure 1 and Figure 2. The theoretical result for (-)-carvone f H m o (liq)theor = −185.9 ± 3.5 kJ·mol−1 (Table 2) was compared with the experimental result, which was directly measured in this work using high-precision combustion calorimetry as follows.
The standard specific energy of combustion c u o (liq) = −38222.5 ± 4.3 J·g−1 of the (-)-carvone was determined from five experiments. The results of the combustion experiments are summarised in Table S8. The c u o (liq) value was used to calculate the experimental standard molar enthalpy of combustion, c H m o (liq) = −5479.1 ± 1.9 kJ·mol−1, which refers to the reaction
C10H14O (liq) + 13 O2 (g) = 10 CO2 (g) + 7 H2O (liq)
The experimental standard molar enthalpy of formation in the liquid state f H m o (liq) = −186.8 ± 2.3 kJ·mol−1 for (-)-carvone is in very good agreement with the theoretical result f H m o (liq)theor = −185.9 ± 3.5 kJ·mol−1 (Table 2). Such a good agreement supports the results derived for the carvone derivatives in Steps I and II.
The enthalpies of formation in the liquid phase have now been validated, and the reliable r H m o (liq) values for reactions of carvone hydrogenation (see Figure 1 and Figure 2) were derived according to Equation (2). The resulting reaction enthalpies in the liquid phase at the reference temperature, r H m o (liq, 298 K), are shown in Table 5, column 2.
It is interesting to compare the energetics of the hydrogenation of the double bonds (DB): the endo-, the exo-, and the carbonyl C=O, all three of which are present in carvone.
The enthalpy of Reaction 1, r H m o (liq, 298 K) = −92.7 kJ·mol−1 (Table 5), represents the energetics of endo-DB hydrogenation. For comparison, the hydrogenation enthalpy of cyclohexene (−117.9 kJ·mol−1, calculated from data given in Table S5) which also has the endo-DB in the six-member ring, is significantly more exothermic, reflecting the effect of the conjugation of double bonds in the carvone.
In contrast, the enthalpy of Reaction 4, r H m o (liq, 298 K) = −111.8 kJ·mol−1 (Table 5), which represents the energetics of exo-DB hydrogenation, is indistinguishable from the hydrogenation of exo-DB of α-methyl-styrene (−111.6 kJ·mol−1, calculated from data given in Table S5), demonstrating the absence of conjugation of exo-DB with the double bonds on the carvone ring.
The enthalpy of Reaction 6, r H m o (liq, 298 K) = −56.9 kJ·mol−1 (Table 5), represents the energetics of the carbonyl DB hydrogenation. For comparison, the hydrogenation enthalpy of cyclohexanone (−77.0 kJ·mol−1, Table 5) which has a similar shape, is more exothermic. It is likely that the energetics of Reaction 6 are influenced by the conjugation of all three double bonds in carvone and by steric interactions of the substituents on the cyclohexane ring. This conclusion is supported by the fact that the enthalpy of hydrogenation of tetrahydrocarvone in carvomenthol, r H m o (liq, 298 K) = −70.7 kJ·mol−1 (R3, Table 5), lies energetically between the values of R3 and cyclohexanone (R20). Moreover, the specific effect of the carbonyl double bond in the carvone on the hydrogenation energetics becomes clear when one considers the quite comparable (at the level of −110 kJ·mol−1) hydrogenation enthalpies of reactions R7, R8, R9 and R10, in which the hydroxyl is present on the cyclohexane ring instead of the carbonyl group.
The only endothermic reaction in the network shown in Figure 1 is the dehydrogenation of carvotanacetone to carvacrol, r H m o (liq, 298 K) = 13.9 kJ·mol−1 (R11, Table 5). However, this very moderate energetics indicates that this dehydrogenation reaction is not far from equilibrium conditions even at low temperatures.
The most surprising aspect of the network shown in Figure 1 is that the isomerisation of carvone to carvacrol, r H m o (liq, 298 K) = −97.9 kJ·mol−1 (R12, Table 5), is strongly exothermic, which is not typical for this type of reaction. The very strong energetics indicates that the isomerisation reaction competes with the hydrogenation reactions even at low temperatures.
The reaction enthalpies of the cis-trans isomerisation of the carvone derivatives (see Figure 2) as a side reaction of carvone hydrogenation reactions are generally not exceeding 10 kJ·mol−1 (see R13, R14 and R15 in Table 5). The low energetics is also typical for this type of isomerisation reaction.
The network of hydrogenation reactions shown in Figure 1 illustrates the gradual stepwise hydrogenation of carvone. However, it is also of interest to compare the energetics of the partial hydrogenation of endo-DB and exo-DB in carvone (see Figure 6), r H m o (liq, 298 K) = −207.4 kJ·mol−1 (R16, Table 5) and the complete hydrogenation of all three bonds in carvone, r H m o (liq, 298 K) = −278.1 kJ·mol−1 (R17, Table 5).
It is obvious that the reaction enthalpy of R17 is greater than the reaction enthalpy of R16 exactly by the contribution to the energetics of C=O bond hydrogenation (see R3 in Table 5). Since different numbers of hydrogen molecules are involved R16 and R17, the reaction enthalpies are usually related to the amount of hydrogen consumed (as kJ·mol−1/H2 unit) in order to be able to correctly compare the energetics of hydrogenation reactions, taking the stoichiometry into account. In these new units, the hydrogenation enthalpy for R16, r H m o (liq, 298 K) = −103.7 kJ·mol−1/H2 and r H m o (liq, 298 K) = −92.4 kJ·mol−1/H2 for R17 were recalculated for comparison. Indeed, the two results do not differ significantly, but both hydrogenation enthalpies of the aliphatic carvone are significantly higher than those of its aromatic isomer carvacrol: r H m o (liq, 298 K) = −180.2/3 = −60.1 kJ·mol−1/H2 for R18 (see Table 5). This significant difference between the energetics of the hydrogenation of aromatic and aliphatic unsaturated systems indicates that the reaction is facilitated for the aliphatic unsaturated molecules.

4.3.2. Liquid-Phase Reaction Entropies

The reaction entropies, r S m o (liq, 298 K), which were calculated for R1 to R20 according to Equation (4), are summarised in Table 5, column 3. Admittedly, the reaction entropies (in relation to the amount of hydrogen consumed) do not differ significantly. As can be seen from Table 5, the r S m o (liq, 298 K) values are quite similar for all hydrogenation reactions of endo- or exo-DB (see R1, R2, R4, R5, R7 to R10) at the 120 J·mol−1·K−1 level. The reactions R3, R4 as well as R19 (hydrogenation of acetone) and R20 (hydrogenation of cyclohexanone), in which a C=O bond is involved, are somewhat higher, but they are also very close to the 150 J·mol−1·K−1 level. As expected, the reaction entropies, r S m o (liq, 298 K), for the isomerisation reactions R12 to R15 are very small and do not exceed 10 J·mol−1·K−1 (see Table 5).
In this way, both the enthalpic (Section 4.3.1) and entropic (Section 4.3.2) contributions to the Gibbs–Helmholtz equation for the network reactions of carvone hydrogenation were derived (see Table 5) according to Equation (1), and the feasibility of these chemical reactions in terms of the standard molar Gibbs energies can be performed.

4.3.3. Gibbs Energies and Thermodynamic Analysis of the Hydrogenation of Carvone

From the experimental studies on the hydrogenation of carvone in the presence of various catalysts, the following order of decreasing reactivity was determined empirically: exo-DB > endo DB > C=O [17]. This sequence corresponds to the preferred formation of unsaturated and saturated ketones (see R1, R2 and R4, R5 in Figure 1) as the main products in the hydrogenation of carvone and the low possibility of forming the hydroxyl derivatives of carvone (see R3 and R6 to R10). The r G m o (liq, 298 K) values calculated for reactions 1 to 15 fundamentally support and quantify (see Table 5) these empirical observations. Indeed, the r G m o (liq, 298 K) = −75.3 kJ·mol−1 for the hydrogenation of the exo-DB of carvone to carvotaneacetone (R4) is noticeably greater than the r G m o (liq, 298 K) = −58.2 kJ·mol−1 for the hydrogenation of the endo-DB of carvone to dihydrocarvone (R1). At the same time, the r G m o (liq, 298 K) = −7.0 kJ·mol−1 for the hydrogenation of the C=O from carvone to carveol (R4) is drastically lower than that for endo- and exo-DB. Thus, the results of the order of thermodynamic stability correspond exactly to the empirically determined order of decreasing reactivity.
The same conclusion can be drawn for the hydrogenation of the remaining hydrogen bonds in the second step: the r G m o (liq, 298 K) = −76.8 kJ·mol−1 for the hydrogenation of the exo-DB of dihydrocarvone to tetrahydrocarvone (R2) is greater compared to r G m o (liq, 298 K) = −58.7 kJ·mol−1 for the hydrogenation of the endo-DB of carvotanacetone to tetrahydrocarvone (R5). Also, the r G m o (liq, 298 K) = −23.4 kJ·mol−1 for the hydrogenation of the C=O from tetrahydrocarvone to carvomenthol (R3).
Interestingly, although the formation of the hydroxyl derivatives of carvone is unfavourable, their thermodynamic stabilities are quite similar: the r G m o (liq, 298 K) values of ≈ −75 kJ·mol−1 are the same for the hydrogenation of the exo-DB in dihydrocarveol (R8) and carveol (R9) and for the hydrogenation of the endo-DB in carvotanacetol (R10). However, the same trend has already been discussed for the reaction enthalpies of the hydroxyl derivatives of carvone (see Section 4.3.1).
The r G m o (liq, 298 K) values of the cis-trans isomerisation reactions R13, R14, and R15 are only slightly negative (see Table 5, column 5), which indicates their general feasibility at 298 K. Most strikingly, the r G m o (liq, 298 K) = −96.9 kJ·mol−1 for the isomerisation of carvone to carvacrol (R12) is by far the largest in the reaction network shown in Figure 1. This proves that the highly selective direct conversion of carvone to carvacrol is possible provided a suitable catalyst is found that does not promote the partial hydrogenation reactions. The additional formation of carvacrol, as the end product of the network reactions in Figure 1, can be achieved by the dehydrogenation of carvotanacetone to carvacrol according to R11. The r G m o (liq, 298 K) = −21.6 kJ·mol−1 for this reaction is large enough to make a sufficient contribution.
As part of the thermodynamic analysis of chemical reactions, it is usually discussed whether the process is enthalpically or entropically driven. This assessment is usually made by comparing the size of the two contributions to the Gibbs energy according to Equation (1) as shown in Table S9. Interestingly, the entropy contributions ( T × r S m o ≈ 35 kJ·mol−1) to the Gibbs energy of all hydrogenation reactions in Figure 1 are approximately constant and amount to only 1/3 compared to the enthalpic term (see Table S9). Consequently, from a thermodynamic point of view, these reactions equilibrium are mainly determined by the enthalpic term.
Basically, the thermodynamic calculations are initially carried out at the reference temperature T = 298 K. However, the typical temperature range for the reaction network is between 293 and 500 K. Using the standard molar isobaric heat capacities C p , m o of the reactions’ participants (see Table S4), the Gibbs energies and equilibrium constants of reactions 1 to 15 at 400 K and 500 K were calculated using Equations (7) and (8). The results are shown in Table 5 and help to understand whether increasing the temperature shifts the equilibrium toward an increase or decrease in yield. As the temperature increases, the equilibrium constants for R1 decrease considerably. In contrast, as the temperature increases, the equilibrium constants for R11 increase. The equilibrium of R12, R14 and R15 is practically independent of the temperature. Only in the cis-trans isomerisation of dihydrocarvone does the equilibrium constant increase with increasing temperature.

4.3.4. Hydrogenation of Carvone: Catalytic Aspects

Quantitative thermodynamic analysis has shown that negative Gibbs energies are recorded for all reactions that may occur during the hydrogenation of carvone (Figure 1 and Figure 2), although these differ significantly in magnitude. This means that all these reactions are thermodynamically feasible in a wide range from ambient temperature to elevated temperatures. Therefore, all these reactions definitely take place under kinetic and not thermodynamic control. Nevertheless, the numerical r G m o (liq, 298 K) values can help to establish the chemoselectivity of catalysts used to convert carvone to either carvacarol or to dihydro- and terahydrocarvone, either in carvotanacetone or carveol.
For example, Benavente et al. [5] demonstrated the exclusive formation of carvacrol with the complete carvone conversion via Pd/Al2O3 and Pd/C. In the same study, it was shown that Pd/CeO2 promotes the formation of carveol due to the activation of C=O bonds at certain oxygen vacancies of the catalyst.
The selective hydrogenation of carvone to dihydrocarvone was carried out at 373 K with an Au/TiO2 catalyst [4]. The gold catalyst showed high activity as well as stereo- and chemoselectivity in the conjugated endo-DB hydrogenation with a predominant formation of trans-dihydrocarvone. The overall selectivity to dihydrocarvone 62% was achieved with an almost complete carvone conversion (90%) after 13 h, whereby the ratio of trans- to cis-dihydrocarvone was about 1.8 (see reaction R13 in Figure 1). According to the calculations with the Gibbs energy of R13 given in Table 5, a ratio of trans- to cis-dihydrocarvone of 3.7 is expected when the system reaches equilibrium. This means that extending the reaction time or doubling the amount of catalyst could increase the carvone conversion and the yield of dihydrocarvones.

5. Conclusions

The thermodynamic analysis of the network of carvone hydrogenation reactions shows that even if the hydrogenations reactions run under kinetic control, the quantitative knowledge of thermodynamic characteristics of reactions is essential to assess the level of equilibrium constants at a given temperature in order to manipulate the amount of catalyst, time, and temperature and achieve sufficient yields of industrially important products of carvone hydrogenation.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/chemistry6040042/s1, Figure S1: Structures of the most stable conformers of carvone derivatives as calculated using the G4 method; Figure S2: Vapour pressure - temperature dependence for (-)-carvone; Figure S3: Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of dihydrocarvone using the 2-methyl-cyclohexanone as the “centerpiece” (left). Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of carvonacetol using the 2-methyl-cyclohex-2-en-1-ol as the “centerpiece” (right). All values in kJ·mol−1. The experimental values of vaporisation enthalpies of the “centerpieces” are given in Tables S2 and S5; Figure S4: Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of dihydrocarveol using the 2-methyl-cyclohexanol as the “centerpiece” (left). Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of carvomenthol using the 2-methyl-cyclohexanol as the “centerpiece” (right). All values in kJ·mol−1. The experimental values of vaporisation enthalpies of the “centerpieces” are given in Table S5; Figure S5: Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of carveol using the 2-methyl-2-cyclohex-2-en-1-ol as the “centerpiece” (left). Calculation the enthalpy of vaporisation, l g H m o (298.15 K), of tcarvotanacetone using the 2-methyl-2-cyclohexen-1-one as the “centerpiece” (right). All values in kJ·mol−1. The experimental values of vaporisation enthalpies of the “centerpieces” are given in Tables S2 and S5; Figure S6: Reactions 19 and 20 (see Table 5) for hydrogenation of acetone and cyclohexanone; Table S1: Provenance and purity of the materials; Table S2: 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 S3: Absolute vapour pressures p, and standard molar vaporisation enthalpies and entropies determined using the transpiration method; Table S4: Compilation of data on molar heat capacities C p , m o (liq or cr) and heat capacity differences l , cr g C p , m o (in J·K−1·mol−1) at T = 298.15 K; Table S5: The enthalpies of vaporisation, l g H m o , and enthalpies of formation, f H m o , of the auxiliary reference compounds at T = 298.15 K ( = 0.1 MPa) in kJ·mol−1; Table S6: Compilation of the standard molar enthalpies of vaporisation. l g H m o . of the carvone derivatives using the “centerpiece” approach (in kJ·mol−1); Table S7: 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 S8: Compilation of combustion results at T = 298.15 K ( = 0.1 MPa) obtained for (-)-carvone; Table S9: Calculation of the liquid phase thermodynamic properties for reactions R1 to R20 (see Table 5) [10,13,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44].

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, A.P.G., S.V.V.; resources, A.P.G., 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, A.P.G., S.P.V. All authors have read and agreed to the published version of the manuscript.

Funding

AAS acknowledges gratefully 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). APG thanks GZ FSZE-2022-0002. SPV acknowledges financial support from German Science Foundation in the frame of SPP 1807 “Control of London Dispersion Interactions in Molecular Chemistry”, grant VE 265-9/2. The work is carried out in accordance with the Strategic Academic Leadership Program “Priority 2030” of the Kazan Federal University of the Government of the Russian Federation.

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. Jiang, Y.; Wang, X.; Cao, Q.; Dong, L.; Guan, J.; Mu, X. Sustainable Production of Bulk Chemicals; Xian, M., Ed.; Springer: Dordrecht, The Netherlands, 2016; pp. 19–49. [Google Scholar]
  2. Bicas, J.L.; Dionísio, A.P.; Pastore, G.M. Bio-Oxidation of Terpenes: An Approach for the Flavor Industry. Chem. Rev. 2009, 109, 4518–4531. [Google Scholar] [CrossRef]
  3. Ascue Avalos, G.A.; Toogood, H.S.; Tait, S.; Messiha, H.L.; Scrutton, N.S. From Bugs to Bioplastics: Total (+)-Dihydrocarvide Biosynthesis by Engineered Escherichia coli. ChemBioChem 2019, 20, 785–792. [Google Scholar] [CrossRef]
  4. Demidova, Y.S.; Suslov, E.V.; Simakova, O.A.; Simakova, I.L.; Volcho, K.P.; Salakhutdinov, N.F.; Murzin, D.Y. Selective Carvone Hydrogenation to Dihydrocarvone over Titania Supported Gold Catalyst. Catal. Today 2015, 241, 189–194. [Google Scholar] [CrossRef]
  5. Benavente, P.; Cárdenas-Lizana, F.; Keane, M.A. Selective Production of Carvacrol from Carvone over Supported Pd Catalysts. Catal. Commun. 2017, 96, 37–40. [Google Scholar] [CrossRef]
  6. 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]
  7. Curtiss, L.A.; Redfern, P.C.; Raghavachari, K. Gaussian-4 Theory. J. Chem. Phys. 2007, 126, 084108. [Google Scholar] [CrossRef] [PubMed]
  8. 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]
  9. 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]
  10. 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]
  11. Verevkin, S.P. Weaving a Web of Reliable Thermochemistry around Lignin Building Blocks: Phenol, Benzaldehyde, and Anisole. J. Therm. Anal. Calorim. 2022, 147, 6073–6085. [Google Scholar] [CrossRef]
  12. Verevkin, S.P. Thermochemistry of Phenols: Buttress Effects in Sterically Hindered Phenols. J. Chem. Thermodyn. 1999, 31, 1397–1416. [Google Scholar] [CrossRef]
  13. Štejfa, V.; Fulem, M.; Růžička, K. Thermodynamic Study of Selected Aromatic Monoterpenoids. J. Mol. Liq. 2023, 380, 121724. [Google Scholar] [CrossRef]
  14. Samarov, A.A.; Verevkin, S.P. Comprehensive Thermodynamic Study of Alkyl-Biphenyls as a Promising Liquid Organic Hydrogen Carriers. J. Chem. Thermodyn. 2022, 174, 106872. [Google Scholar] [CrossRef]
  15. Verevkin, S.P.; Andreeva, I.V.; Zherikova, K.V.; Pimerzin, A.A. Prediction of Thermodynamic Properties: Centerpiece Approach—How Do We Avoid Confusion and Get Reliable Results? J. Therm. Anal. Calorim. 2022, 147, 8525–8534. [Google Scholar] [CrossRef]
  16. Gallis, H.E.; van Miltenburg, J.C.; Oonk, H.A.J.; van der Eerden, J.P. Mixtures of D- and l-Carvone. III. Thermodynamic Properties of l-Carvone. Thermochim. Acta 1996, 286, 307–319. [Google Scholar] [CrossRef]
  17. Barbaro, P.; Liguori, F. Heterogenized Homogeneous Catalysts for Fine Chemicals Production: Materials and Processes; Springer: Dordrecht, The Netherlands, 2010. [Google Scholar]
  18. SciFinder. Chemical Abstracts Service. Available online: http://scifinder.cas.org/ (accessed on 25 June 2024).
  19. Reaxys. Available online: https://www.reaxys.com/ (accessed on 25 June 2024).
  20. Guidechem. Chemical Network. Available online: https://www.guidechem.com (accessed on 25 June 2024).
  21. Kulikov, D.; Verevkin, S.P.; Heintz, A. Determination of Vapor Pressures and Vaporization Enthalpies of the Aliphatic Branched C 5 and C 6 Alcohols. J. Chem. Eng. Data 2001, 46, 1593–1600. [Google Scholar] [CrossRef]
  22. 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]
  23. Verevkin, S.P.; Sazonova, A.Y.; Emel’yanenko, V.N.; Zaitsau, D.H.; Varfolomeev, M.A.; Solomonov, B.N.; Zherikova, K.V. Thermochemistry of Halogen-Substituted Methylbenzenes. J. Chem. Eng. Data 2015, 60, 89–103. [Google Scholar] [CrossRef]
  24. Emel’yanenko, V.N.; Verevkin, S.P. Benchmark Thermodynamic Properties of 1,3-Propanediol: Comprehensive Experimental and Theoretical Study. J. Chem. Thermodyn. 2015, 85, 111–119. [Google Scholar] [CrossRef]
  25. Acree, W.; Chickos, J.S. Phase Transition Enthalpy Measurements of Organic and Organometallic Compounds and Ionic Liquids. Sublimation, Vaporization, and Fusion Enthalpies from 1880 to 2015. Part 2. C11-C192. J. Phys. Chem. Ref. Data 2017, 46, 013104. [Google Scholar] [CrossRef]
  26. Chickos, J.S.; Hosseini, S.; Hesse, D.G.; Liebman, J.F. Heat Capacity Corrections to a Standard State: A Comparison of New and Some Literature Methods for Organic Liquids and Solids. Struct. Chem. 1993, 4, 271–278. [Google Scholar] [CrossRef]
  27. Clarke, E.C.W.; Glew, D.N. Evaluation of Thermodynamic Functions from Equilibrium Constants. Trans. Faraday Soc. 1966, 62, 539. [Google Scholar] [CrossRef]
  28. Ambrose, D.; Ghiassee, N. Vapour Pressures and Critical Temperatures and Critical Pressures of C5 and C6 Cyclic Alcohols and Ketones. J. Chem. Thermodyn. 1987, 19, 903–909. [Google Scholar] [CrossRef]
  29. Aucejo, A.; Monton, J.B.; Munoz, R.; Sanchotello, M. Isobaric Vapor-Liquid Equilibrium Data for the Cyclohexanone + N-Methylacetamide System. J. Chem. Eng. Data 1993, 38, 160–162. [Google Scholar] [CrossRef]
  30. Palczewska-Tulińska, M.; Oracz, P. Vapor Pressures of Hexanal, 2-Methylcyclohexanone, and 2-Cyclohexen-1-One. J. Chem. Eng. Data 2006, 51, 639–641. [Google Scholar] [CrossRef]
  31. Cruz Burguet, M.; Monton, J.B.; Sanchotello, M.; Vazquez, M.I. Vapor-Liquid Equilibria of Cyclohexanone + Cyclohexanol and Cyclohexanone + 2-Methylcyclohexanone Systems at 4.00 and 26.66 KPa. J. Chem. Eng. Data 1993, 38, 328–331. [Google Scholar] [CrossRef]
  32. Aucejo, A.; Cruz Burguet, M.; Munoz, R.; Vazquez, M.I. Isobaric Vapor-Liquid Equilibria for Cyclohexanone + 3-Methylcyclohexanone or 4-Methylcyclohexanone Systems at 4.00 and 26.66 KPa. J. Chem. Eng. Data 1993, 38, 379–382. [Google Scholar] [CrossRef]
  33. Herz, W.; Bloch, W. Physikalisch-Chemische Untersuchungen an Verbindungen Der Cyklohexanreihe. Z. Phys. Chem. 1924, 110, 23–39. [Google Scholar] [CrossRef]
  34. Vilas-Boas, S.M.; Pokorný, V.; Štejfa, V.; Ferreira, O.; Pinho, S.P.; Růžička, K.; Fulem, M. Vapor Pressure and Thermophysical Properties of Eugenol and (+)-Carvone. Fluid Phase Equilib. 2019, 499, 112248. [Google Scholar] [CrossRef]
  35. Majer, V.; Svoboda, V. Enthalpies of Vaporization of Organic Compounds: A Critical Review and Data Compilation; Blackwell Scientific Publications: Oxford, UK, 1985. [Google Scholar]
  36. Goodwin, S.R.; Newsham, D.M.T. Vapor-Liquid Equilibriums in Mixtures of o-Methylcyclohexanol and o-Methylcyclohexyl Acetate from 50-296 Mm Hg. J. Chem. Eng. Data 1975, 20, 178–179. [Google Scholar] [CrossRef]
  37. 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]
  38. Verevkin, S.P.; Schick, C. Substituent Effects on the Benzene Ring. Determination of the Intramolecular Interactions of Substituents in Tert-Alkyl-Substituted Catechols from Thermochemical Measurements. J. Chem. Eng. Data 2000, 45, 946–952. [Google Scholar] [CrossRef]
  39. 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]
  40. 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]
  41. Verevkin, S.P. Thermochemistry of Phenols: Quantification of the Ortho-, Para-, and Meta-Interactions in Tert-Alkyl Substituted Phenols. J. Chem. Thermodyn. 1999, 31, 559–585. [Google Scholar] [CrossRef]
  42. Bardyshev, I.I.; Kozlov, N.G.; Vyalimyaé, T.K.; Pekhk, T.I. Synthesis and Study of the Structure of New N-Substituted 2-Methyl-5-(1-Methylethyl) Cyclohexylamines. Chem. Nat. Compd. 1980, 16, 397–402. [Google Scholar] [CrossRef]
  43. 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 Corporation: New York, NY, USA, 1989. [Google Scholar]
  44. 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]
Figure 1. General network of carvone hydrogenation reactions.
Figure 1. General network of carvone hydrogenation reactions.
Chemistry 06 00042 g001
Figure 2. Cis-trans isomerisation of the carvone derivatives as a concomitant of carvone hydrogenation reactions.
Figure 2. Cis-trans isomerisation of the carvone derivatives as a concomitant of carvone hydrogenation reactions.
Chemistry 06 00042 g002
Figure 3. Structures of the most stable conformers of carvone (CAS 99-49-0) as calculated using the G4 method.
Figure 3. Structures of the most stable conformers of carvone (CAS 99-49-0) as calculated using the G4 method.
Chemistry 06 00042 g003
Figure 4. Calculating of l g H m o (298 K). Development of the contribution H→iPr(en) for the exchange of the H-atom in the cyclohexane ring for the iso-propenyl substituent (left). Development of the contribution H→iPr for the exchange of the H-atom in the cyclohexane ring for the iso-propyl substituent (right). All values in kJ·mol−1.
Figure 4. Calculating of l g H m o (298 K). Development of the contribution H→iPr(en) for the exchange of the H-atom in the cyclohexane ring for the iso-propenyl substituent (left). Development of the contribution H→iPr for the exchange of the H-atom in the cyclohexane ring for the iso-propyl substituent (right). All values in kJ·mol−1.
Chemistry 06 00042 g004
Figure 5. Calculation of the enthalpy of vaporisation, l g H m o (298 K), of carvone using the 2-methyl-2-cyclohexen-1-one as the “centrepiece” (left). Calculation the enthalpy of vaporisation, l g H m o (298 K), of tetrahydrocarvone using the 2-methyl-cyclohexanone as the “centrepiece” (right). All values in kJ·mol−1. The experimental values of vaporisation enthalpies of the “centrepieces” are given in Tables S2 and S5.
Figure 5. Calculation of the enthalpy of vaporisation, l g H m o (298 K), of carvone using the 2-methyl-2-cyclohexen-1-one as the “centrepiece” (left). Calculation the enthalpy of vaporisation, l g H m o (298 K), of tetrahydrocarvone using the 2-methyl-cyclohexanone as the “centrepiece” (right). All values in kJ·mol−1. The experimental values of vaporisation enthalpies of the “centrepieces” are given in Tables S2 and S5.
Chemistry 06 00042 g005
Figure 6. Reactions of the partial and complete hydrogenation of double bonds in carvone and carvacrol (reactions 16–18 in Table 5).
Figure 6. Reactions of the partial and complete hydrogenation of double bonds in carvone and carvacrol (reactions 16–18 in Table 5).
Chemistry 06 00042 g006
Table 1. Correlation of the G4 calculated and experimental gas-phase enthalpies of formation f H m o (g) at T = 298 K and = 0.1 MPa (in kJ·mol−1).
Table 1. Correlation of the G4 calculated and experimental gas-phase enthalpies of formation f H m o (g) at T = 298 K and = 0.1 MPa (in kJ·mol−1).
CompoundAT aExp bATcorr cΔ d
cyclohexanol−289.8−290.0 ± 2.1 [10]−288.8−1.2
butanol−274.4−275.0 ± 0.4 [10]−273.9−1.1
pentanol−295.7−294.7 ± 0.5 [10]−294.90.2
hexanol−317.0−315.8 ± 0.6 [10]−315.80.0
heptanol−338.4−336.4 ± 1.0 [10]−336.70.3
octanol−359.7−356.5 ± 1.5 [10]−357.71.2
cylopentanone−192.9−192.1 ± 1.8 [10]−193.91.8
cyclohexanone−228.4−226.1 ± 2.1 [10]−228.72.6
phenol−90.8−92.5 ± 1.2 [11]−93.61.1
2-methylphenol−127.0−128.4 ± 0.9 [10]−129.10.7
4-methylphenol−122.0−125.8 ± 1.5 [10]−124.2−1.6
2,3-dimethylphenol−156.6−157.2 ± 1.1 [10]−158.21.0
2,5-dimethylphenol−159.5−161.6 ± 1.0 [10]−161.1−0.5
2,6-dimethylphenol−161.1−161.8 ± 0.5 [10]−162.60.8
3,4-dimethylphenol−156.2−156.6 ± 0.6 [10]−157.81.2
3,5-dimethylphenol−157.7−161.6 ± 0.7 [10]−159.3−2.2
2,3,6-trimethylphenol−190.3−192.9 ± 1.4 [10]−191.3−1.6
2,4,6-trimethylphenol−192.0−192.4 ± 1.3 [12]−193.00.6
carvacrol (5-isopropyl-2-methylphenol)−210.8 −211.4
a Calculated according to atomisation reaction in Equation (9). b Experimental values from the literature. c Calculated according to Equation (10). d Difference between columns 3 and 4.
Table 4. 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 carvone derivatives (all values at T = 298 K in J·mol−1·K−1).
Table 4. 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 carvone derivatives (all values at T = 298 K in J·mol−1·K−1).
CASCompound l g S m o  a S m o (g) b S m o (liq) c
6485-40-1(-)-carvone130.8 ± 0.7441310.2 d
99-48-9trans-carveol159.5 ± 2.8434274.5
3792-53-8cis-(±)-dihydrocarvone122.8 ± 2.8439316.2
5524-05-0trans-(±)-dihydrocarvone123.0 ± 2.0448325.0
499-71-8carvotanacetone126.6 ± 2.5445318.4
499-70-7trans-(±)-tetrahydrocarvone (carvomenthone)124.8 ± 1.8450325.2
cis-(±)-tetrahydrocarvone125.9 ± 1.9453327.1
619-01-2trans-dihydrocarveol159.8 ± 2.0451291.2
22567-21-1cis-dihydrocarveol155.4 ± 4.5449293.6
536-30-1carvotanacetol163.2 ± 4.5450286.8
499-75-2carvacrol157.1 ± 0.5 [13]455 [13]297.9
a From the vapour pressure measurements compiled in Tables S2 and S3. b Calculated with the G4 method [7]. c Calculated according to Equation (3) using entries from columns 3 and 4 from this table. d For comparison, the experimental value S m o (liq) =314.3 ± 0.6 J·mol−1·K−1 [16] was measured for (-)-carvone with adiabatic calorimetry.
Table 5. Calculation of the liquid-phase thermodynamic properties for reactions R1 to R20 (see Figure 1, Figure 2 and Figure S6).
Table 5. Calculation of the liquid-phase thermodynamic properties for reactions R1 to R20 (see Figure 1, Figure 2 and Figure S6).
R r H m o  a
(298 K)
r S m o  b
(298 K)
r C p , m o  c
(298 K)
r G m o  d
(298 K)
r G m o  e
(400 K)
r G m o  f
(500 K)
lnKa g
(298 K)
lnKa g
(400 K)
lnKa g
(500 K)
kJ·mol−1J·mol−1·K−1J·mol−1·K−1kJ·mol−1kJ·mol−1kJ·mol−1
partial hydrogenation
1−92.7−115.9−29.0−58.2−45.9−33.123.518.513.4
2−114.7−130.5−19.0−75.8−62.2−48.430.625.119.5
3−70.7−158.6−0.7−23.4−7.38.69.52.9−3.5
4−111.8−122.5−21.6−75.3−62.5−49.330.425.219.9
5−95.6−123.9−26.4−58.7−45.7−32.123.718.413.0
6−56.9−166.4−3.3−7.010.026.82.8−4.0−10.8
7−110.4−114.0−3.9−76.4−64.8−53.230.826.121.5
8−111.1−124.6−41.5−74.0−60.7−46.429.924.518.7
9−110.3−118.4−19.0−75.0−62.7−50.030.325.320.2
10−111.2−120.2−26.4−75.4−62.8−49.630.425.320.0
dehydrogenation
1113.9119.267.2−21.6−34.7−49.68.714.020.0
isomerisation
12−97.9−3.345.6−96.9−97.2−98.939.139.239.9
13−1.58.8-−4.1−5.0−5.91.72.02.4
14−1.6−1.9-−1.0−0.8−0.70.40.30.3
15−9.5−2.4-−8.8−8.5−8.33.53.43.4
full hydrogenation
16−207.4−246.4−48.0−133.9−108.2−81.454.143.732.9
17−278.1−404.9−48.7−157.4−115.4−72.863.546.629.4
18−180.2−401.6−94.3−60.5−18.226.124.47.4−10.5
simple ketones h
19−70.0−150.56.9−25.1−9.94.810.14.0−2.0
20−77.0−155.87.6−30.6−14.80.412.36.0−0.2
a Calculated according to Equation (2) from the standard molar enthalpies of formation of the reactants from Table 2. b Calculated according to Equation (4) from the standard molar entropies of the reactants from Table 4. c Difference in the standard molar heat capacities of reactants calculated from the values in Table S4. d Calculated according to Equation (1) from the results given in columns 2 and 3 and referenced to 298 K. e Calculated according to Equation (1) from the results given in columns 2 and 3 and adjusted to 400 K according to Kirchhoff’s Law (see text). f Calculated according to Equation (1) from the results given in column 2 and 3 and adjusted to 500 K according to Kirchhoff´s Law (see text). g The logarithm of the thermodynamic equilibrium constant Ka, calculated according to Equation (1) for the corresponding temperature. h The reactions are shown in Figure S6.
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

Samarov, A.A.; Vostrikov, S.V.; Glotov, A.P.; Verevkin, S.P. Thermodynamic Exercises for the Kinetically Controlled Hydrogenation of Carvone. Chemistry 2024, 6, 706-722. https://doi.org/10.3390/chemistry6040042

AMA Style

Samarov AA, Vostrikov SV, Glotov AP, Verevkin SP. Thermodynamic Exercises for the Kinetically Controlled Hydrogenation of Carvone. Chemistry. 2024; 6(4):706-722. https://doi.org/10.3390/chemistry6040042

Chicago/Turabian Style

Samarov, Artemiy A., Sergey V. Vostrikov, Aleksandr P. Glotov, and Sergey P. Verevkin. 2024. "Thermodynamic Exercises for the Kinetically Controlled Hydrogenation of Carvone" Chemistry 6, no. 4: 706-722. https://doi.org/10.3390/chemistry6040042

APA Style

Samarov, A. A., Vostrikov, S. V., Glotov, A. P., & Verevkin, S. P. (2024). Thermodynamic Exercises for the Kinetically Controlled Hydrogenation of Carvone. Chemistry, 6(4), 706-722. https://doi.org/10.3390/chemistry6040042

Article Metrics

Back to TopTop