Electrochemical hydronation of model compounds
Benzoic acid electrochemical hydrogenation tests were carried out using a H-type electrochemical cell (Fig. 1). At an initial BA concentration of 10 mM, 94.70% conversion with 100% selectivity to cyclohexane carboxylic acid (CCA) was achieved within 4 hours at 9% Faradaic efficiency (Fig. 2a-b). This high selectivity was consistently observed across all benzoic acid experiments in this study. Increasing the initial concentration from 10 to 20 mM resulted in a decreased conversion rate and increased FE (75.4% and 14%, respectively). The initial BA concentration significantly impacted both conversion rate and Faradaic efficiency (FE). Increasing the concentration from 10 to 20 mM enhanced both parameters, likely due to increased collision frequency between H+ ions and organic molecules. However, further concentration increases to 30 and 50 mM led to decreased conversion rates, despite slightly improved FEs. This suggests a mass transport limitation at higher concentrations, where the availability of electrons becomes insufficient to sustain the reaction rate. Based on this finding, the optimal initial BA concentration for this study was determined to be 20 mM and applied to other experiments across this study.
Temperature significantly influenced ECH performance. To evaluate the effect of temperature on the process, and determine the optimal condition, we tested temperatures of 40, 50, 55, 60 and 80°C (Fig. 2c-d), maintaining the same current density at 100.0 mA cm− 2 and setting the initial concentration of BA to 20 mM. At lower temperatures (40°C), the system's conductivity was reduced, leading to insufficient activation energy for bond cleavage 11,13,24 and consequently lower conversion (49.8%) and Faradaic efficiency (9%). Increasing the temperature to 50–60°C enhanced both conversion and FE, with optimal performance observed at 55°C. However, further heating to 80°C promoted hydrogen desorption from the cathode, resulting in decreased ECH activity (78.74% conversion, 14% FE).
The influence of current density on ECH was investigated using five current values (0.1, 0.2, 0.4, and 0.6 A), corresponding to current densities of 12.5, 25.0, 50.0, 100.0, and 150.0 mA cm⁻², respectively (based on a 4.0 cm² cathode area). The initial BA concentration was maintained at 20 mM, and the temperature at 55°C. Preliminary tests at 0.05 A indicated insufficient proton generation on the cathode surface, limiting reactant collisions and hindering process efficiency.
Increasing the current density from 12.5 to 25.0 mA cm− 2 significantly improved both conversion rates and Faradaic efficiency. At a current density of 50.0 mA cm− 2, the conversion rate increased by 13.55%, but the FE reduced significantly from 53–31%. At higher current densities of 100 and 150 mA cm− 2, higher overpotentials are required (0.340 and 0.416 V vs Ag/AgCl, respectively), increasing HER rate and resulted in the lowest FE values observed in the experiment, at 14% and 10%, respectively. This decline in Faradaic efficiency at higher overpotentials can be attributed to several factors. First, higher overpotentials increase the likelihood of side reactions, which in this process is hydrogen evolution reaction (HER), these compete with the desired reaction and consume part of the applied current without contributing to the target product formation, reducing FE. Additionally, the rapid generation of hydrogen gas at elevated current densities can lead to the accumulation of gas bubbles on the cathode surface, which impedes effective mass transport by blocking active sites, thereby reducing the available surface area for the reaction.11
To assess cathode stability, reproducibility, and process sustainability, the optimized conditions (temperature = 55°C, current = 0.1 A, initial BA concentration = 20 mM) were replicated over five consecutive cycles. The results (Fig. 3) indicated a consistent performance of the PtRu/ACC catalyst over the five cycles of the experiment. Despite the aggressive reaction environment, there is no sign of cathode degradation, which was confirmed by scanning electron microscopy observation (Fig. 9) in addition to assessing its consistency in performance via GC-MS. It is worth noting that the selectivity for CCA remains 100% throughout all cycles.
Phenol and guaiacol were subjected to ECH under the optimized conditions determined for benzoic acid. Phenol was observed to be the most effective under these conditions, achieving 83.17% conversion and 60% Faradaic efficiency with high selectivity (99%) towards cyclohexanol (Fig. 4a-b). In contrast, guaiacol demonstrated lower activity (68.59% conversion, 59% FE) (Fig. 4c-d) with a product distribution of 50% cyclohexanol, 46% 2-methoxycyclohexanol, and 4% 1-methoxycyclohexane, indicating a preference for demethoxylation over dehydration. Notably, complete aromatic ring hydrogenation was observed for all three model compounds, suggesting a common reaction pathway.
While higher conversion rates and FE for benzoic acid ECH have been reported in literature using smaller H-type cells (25 mL vs. 90 mL) 21, direct comparison is challenging due to differences in cell geometry. Smaller cells often exhibit lower resistance and require lower overpotentials, which in turn affects FE. Scaling up to larger H-type cells necessitates increased reactant quantities, higher currents, and overpotentials, potentially impacting efficiency and selectivity. To address these limitations, continuous flow reactors may offer advantages by mitigating mass transfer issues and enabling better process control.
Electrochemical hydrogenation of model compound mixtures
Understanding the effects of having more than one model compound during ECH is essential to progress bio-oil upgrading, which are inherently mixed compounds. To investigate potential interactions (synergistic and antagonistic effects) of each compound in the mixture, a factorial experimental design approach was applied as shown in Table 1, where a “+” signs indicate the compound is present, and a “-” signals the opposite. Notably, the applied overpotential showed minimal variance (0.003V), remaining between 0.205 and 0.208 V vs Ag/AgCl, regardless of the composition or concentration of the model compounds. This indicates that the conductivity of the system is more significantly influenced by the electrolyte and temperature rather than the model compounds, likely due to their low concentration.
Table 1
Results comparison for the ECH of benzoic acid (BA), phenol (P) and guaiacol (G) mixtures. Reaction parameters: current density = 25.0 mA cm-2, temperature = 55°C, electrolyte = H2SO4 1M.
Compound | Entry | Initial concentration (mM) | Conversion (%) | FE (%) |
BA | P | G | | |
+ - - + + - + + + - + | - + - + - + + + - + + | - - + - + + + - + + + | (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) | 20 20 20 40 40 40 60 20 20 20 20 | 77.35 88.75 68.59 64.19 57.27 47.80 42.59 80.02 63.01 74.78 33.64 | 56 64 59 74 66 55 73 58 46 54 24 |
It was observed that increasing the number of unsaturated molecules generally decreased conversion rates, with the exception of mixtures containing all three compounds (Entries 7 and 11). Conversely, FE increased in runs 4, 5 and 7. Figure 2a shows a clear correlation between decreasing conversion rates and increasing initial concentration.
Notably, in entry 4 a significant drop in conversion from 83.05% (average of BA and P) to 64.19% in the mixture was observed, accompanied by a 13% increase in FE. This suggests a possible synergistic interaction between these two compounds, which aligns with previous studies on phenol-containing mixtures (phenol + furfural 25, and phenol + benzaldehyde26,27), where hydrogen-bonded complex formation enhanced the hydrogenation of the co-reactant.
Phenol exhibited the highest conversion rates among the model compounds. In the BA + P mixture (Entry 4), phenol conversion reached 88.57% compared to 39.81% for benzoic acid. Similarly, phenol was preferentially converted in the G + P mixture (Entry 6) at 74.16% compared to 21.43% for guaiacol. This trend persisted in the ternary mixture (Entry 7), with phenol conversion at 75.79% exceeding those of benzoic acid (18.94%) and guaiacol (33.04%). Phenol also demonstrated the highest overall conversion efficiency as a standalone compound.
As a standalone compound, guaiacol consistently exhibited lower conversion rates compared to phenol and benzoic acid. This reduced reactivity is attributed to its complex molecular structure containing both methoxy and phenolic groups, resulting in higher bond energies20. Interestingly, in mixtures, the hydrogenation of guaiacol appears to be favoured when compared to BA. Mixing benzoic acid with guaiacol (entry 5), led to 50% and 65% conversion of each compound, respectively. This value is quite similar to pure guaiacol tests (69%, entry 3), despite the higher initial concentration. This trend was also observed for entry 7 (mixture of all three compounds). All these results are backed by the LSV (linear sweep voltammetry) analysis (see Fig. 8), where the hydrogenation of phenol is favoured over guaiacol, which is, in turn, favoured over benzoic acid.
The ECH of mixtures also led to different products, including, although in small quantities (< 2.0 area %), the full hydro-dehydroxylation of single aromatic compounds into cyclohexane. Guaiacol was also converted into phenol and cyclohexanone when mixed with BA, which complicates the calculations for the conversion rates of phenol. Where, if this behaviour is duplicated in other mixtures, it becomes increasingly difficult to precisely calculate to which extent the present molecules are the ones added into the system as a model compound or produced via ECH. Therefore, the conversion and FE of said processes could be higher than the reported values. The proposed hydrogenation path for guaiacol and phenol in mixtures is shown in Fig. 5.
Tests were carried out to observe the impact of the initial concentration on mixtures. These followed the same ratio between compounds, but the initial total concentration of compounds for ECH was set to 20 mM. For entries with two compounds, each compound contributed 10 mM, and for entries with three compounds, each compound contributed 6.7 mM.
Setting the initial concentration to 20 mM led to similar trends to the ones observed previously. The BA + P mixture (entry 8) remains the best result for a mixture, presenting a considerably high 69.70% conversion for BA and 90.34% for phenol, very similar results to entries 1 and 2, despite the reduction in the initial concentration for each compound (from 20 to 10 mM), once again remarking to possible synergistic effects between these two compounds. Interestingly, although a slight decrease in activity is to be expected, this did not occur for phenol. Phenol also remained the most active species for ECH, achieving 91.44% conversion, compared to 58.12% for guaiacol (entry 10) and 61.60% compared to 24.33% and 26.11% for BA and guaiacol (entry 11), respectively.
Although the mixtures with guaiacol still presented inferior results, it is worth noting that P + G at 20 mM (entry 10) is far superior to P + G at 40 mM (entry 6), possibly due to guaiacol’s resistance to ECH. That is, increasing guaiacol’s concentration meant a greater negative interference in the process than phenol’s positive one, leading to a lesser result. No phenol was detected in entry 9 (BA + G), and there was also no cyclohexane found in any of the experiments. This suggests that the concentration of these compounds may have been too low to detect, or that a larger quantity of reactant was needed to achieve the full conversion.
Factorial design analysis
Two models were developed to study the effects of mixture composition at low and high concentrations. In both cases, the initial model included main effects and interaction terms: BA, P, G, BA + P, BA + G, and P + G, with conversion and FE as responses.
For conversion at higher concentrations (Fig. 6a), significant contributions were observed from BA, BA × P, BA × G, and P × G, while P and G were excluded from the final model. The r2 value for conversion was 0.7267, indicating that the model explains a moderate proportion of the variability in the response. For FE (Fig. 6b), all terms significantly contributed, with an r2 value of 0.9290, indicating a strong model fit.
At lower concentrations, G was excluded from the final model for conversion, while BA and G were removed for FE. The r2 value for conversion was 0.9608 (Fig. 6d), and for FE, r2 was 0.9325 (Fig. 6e), both demonstrating a high correlation with the predicted model.
The conversion profiler (Figs. 6c and 6f) shows that increasing the concentration of the mixture generally leads to higher FE, with one exception. Phenol had the greatest positive impact on both conversion and FE, regardless of concentration. Conversely, Guaiacol decreased both responses in 3 out of 4 cases but led to a slight increase in FE at higher concentrations. These findings reaffirm the observed effects of these two compounds in the ECH process.
Residual plots for both models' responses showed that the residuals appeared random, with no discernible patterns or trends, confirming the model's assumptions. Additionally, no residuals exceeded control limits, indicating a good model fit. The plots of actual versus predicted responses showed a high correlation, confirming the model's predictive capability.
Density Functional Theory (DFT) calculations
DFT calculations were used to explain further the behaviour of benzoic acid on the cathode surface. The material was built with a Pt-Ru atomic ratio of 1:1 herein denoted PtRu(111). BA was placed both with a parallel placement and at a 45° angle into the surface and the structure was relaxed to find the adsorption energy. Figure 7 shows the interaction between them.
The lower energy value calculated for the parallel placement gives a good interpretation of the 100% selectivity of BA into CCA, whereas there is no bond between the carboxylic acid part of the molecule and the surface. The atom's placement shows that this group moves in the “z” direction, bending the molecule and moving away from the surface. That is, there is the formation of π-σ interactions between the aromatic ring and the metal surface, rather than the carboxylic group. There was also a difference in Ead depending on the initial position of the compound, where placing the aromatic ring closer to ruthenium atoms (Fig. 7b) led to a more stable system (Ead2 = -1.12 eV) than closer to platinum atoms (Ead3 = -0.17 eV, Fig. 7c), signifying it as the preferential site.
Interestingly, Fig. 7b shows that a defect is created in the first layer of the cathode’s surface. To accommodate the BA molecule, the two ruthenium atoms which bonded with C1, C3, C4 and C6 moved too far from each other, severing the bond. To determine whether this was a one-time event, the atoms were placed closer to each other until a bond was formed, but the relaxation process would always lead to this defect. It is worth mentioning that from the second layer onwards, this did not occur.
It is widely acknowledged that the Langmuir-Hinshelwood (L-H) adsorption mechanism is followed by organic compounds containing a benzene ring during ECH.21 The process involves H+ being quickly adsorbed into the surface, due to the synergistic effects between Pt and Ru where, in combination with an e−, produces Hads (Vomer step). The next step is the adsorption of BA in a position parallel to the catalyst surface, to form BAads. Subsequently, Hads is transferred into BAads until the full hydrogenation of the aromatic ring is accomplished, followed by the desorption of the product (cyclohexane carboxylic acid). The proposed mechanism for the ECH of benzoic acid over PtRu/ACC is shown in Fig. 8.
The hydrogenation happens in a specific order, following the most stable state of the molecule. Calculations on the one-by-one hydrogenation (Supplementary Fig. 1) of BA show that the first carbon to receive a H+ is C1 (meta), followed by the order C2-C3-C4-C5-C6. Due to the carboxylic acid`s electronegativity, a negatively charged dipole (δ−) is formed which is close to C6, giving a good justification for why it is the least stable position for the first hydrogen (H1), and why it is the last to be hydrogenated. The positive dipole on C6 (ortho), which borders C1 and C5 (meta) and the distance from C3 to the -COOH group set them as the major candidates for H1, at a relatively small 0.056 and 0.0147 eV difference from C1 to C5 and C3, respectively. The same trend follows for the rest of the process, where the neighbouring carbon to the one that receives the hydrogen is always the least stable position, until the ring is fully stabilised.
Material characterisation
Figure 9 shows the Scanning Electron Microscopy (SEM) secondary electron image for PtRu/ACC before use and after 5 cycles of the optimised BA electrochemical hydrogenation process. In both samples, the metal deposition is homogeneous throughout the surface of the carbon cloth. Elemental mapping (by energy dispersive X-ray spectroscopy) was carried out and revealed an even distribution of both Pt (in teal) and Ru (in green) across the surface; this even distribution is beneficial to optimise the reaction. No degradation of the cathode is apparent after 5 reruns, the concentration and distribution of metals remain unaltered.
LSV analysis shows that PtRu/ACC exhibits − 0.203 V vs Ag/AgCl at 10 mA cm− 2, while it is less negative when adding the compound, reaching − 0.176 V, -0.159V and − 0.149 V vs Ag/AgCl for BA, G and P, respectively (Fig. 10). This measurement is important for the well-known activity for hydrogen evolution reaction (HER) at this current density. The voltage change (∆E) of ∆E1 = 0.027V, ∆E2 = 0.044V and ∆E3 = 0.054V shows that the ECH of all the model compounds is prior to HER, the competing reaction. When adding P to BA, the potential went from − 0.176V to -0.160V vs Ag/AgCl, a difference of 0.16V, whereas when adding G to BA, this difference was 0.08V. Adding P to G increased the potential from − 0.159 to -0.156V, a 0.03V change (see Figure S2-4).