Molecular Electrochemical Catalysis of CO-to-Formaldehyde Conversion with a Cobalt Complex.

Formox, a highly energy-intensive process, currently serves as the primary source of formaldehyde (HCHO), for which there is a crucial and steadily growing chemical demand. The alternative electrochemical production of HCHO from C1 carbon sources such as CO2 and CO is still in its early stages, with even the few identified cases lacking mechanistic rationalization. In this study, we demonstrate that cobalt phthalocyanine (CoPc) immobilized on multiwalled carbon nanotubes (MW-CNTs) constitutes an excellent electrocatalytic system for producing HCHO with productivity through the direct reduction of CO, the two-electron reduction product of CO2. By carefully adjusting both the pH and the applied potential, we identified conditions that enable the production of HCHO with a partial current density of 0.64 mA cm-2 (17.5% Faradaic efficiency, FE) and a total FE of 61.2% for the liquid products (formaldehyde and methanol). A reduction mechanism is proposed.

Aliphatic Ketone Claisen Rearrangement: Troubleshooting the Transetherification Step by Identifying a Stable Acid Catalyst.

After optimization for retention of catalytic activity, 4-chlorobenzoic acid emerged as the optimal catalyst for the aliphatic ketone Claisen rearrangement. The optimal catalyst enables a one-pot, metal-free, catalytic protocol from allylic alcohols to γ,δ-unsaturated ketones. The optimized process tolerates a range of substrates, including substituents with acid-labile protecting groups. Reaction monitoring and DFT studies of the aliphatic ketone Claisen process agree that the ultimate rearrangement step typically has the highest activation barrier.

Exploration of Vitamin B6-Based Redox-Active Pyridinium Salts Towards the Application in Aqueous Organic Flow Batteries.

Pyridoxal hydrochloride, a vitamin B6 vitamer, was synthetically converted to a series of diverse redox-active benzoyl pyridinium salts. Cyclic voltammetry studies demonstrated redox reversibility under basic conditions, and two of the most promising salts were subjected to laboratory-scale flow battery tests involving galvanostatic cycling at 10 mM in 0.1 M NaOH. In these tests, the battery was charged completely, corresponding to the transfer of two electrons to the electrolyte, but no discharge was observed. Both CV analysis and electrochemical simulations confirmed that the redox wave observed in the experimental voltammograms corresponds to a two-electron process. To explain the irreversibility in the battery tests, we conducted bulk electrolysis with the benzoyl pyridinium salts, affording the corresponding benzylic secondary alcohols. Computational studies suggest that the reduction proceeds in three consecutive steps: first electron transfer (ET), then proton-coupled electron transfer (PCET) and finally proton transfer (PT) to give the secondary alcohol. 1H NMR deuterium exchange studies indicated that the last PT step is not reversible in 0.1 M NaOH, rendering the entire redox process irreversible. The apparent reversibility observed in CV at the basic media likely arises from the slow rate of the PT step at the timescale of the measurement.

Ion-Specific Effects on Ion and Polyelectrolyte Solvation.

Ion-specific effects on aqueous solvation of monovalent counter ions, Na + ${^+ }$ , K + ${^+ }$ , Cl - ${^- }$ , and Br - ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl - ${^- }$ and Br - ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

DFT Mechanistic Investigation into Ni(II)-Catalyzed Hydroxylation of Benzene to Phenol by H2O2.

Introduction of oxygen into aromatic C-H bonds is intriguing from both fundamental and practical perspectives. Although the 3d metal-catalyzed hydroxylation of arenes by H2O2 has been developed by several prominent researchers, a definitive mechanism for these crucial transformations remains elusive. Herein, density functional theory calculations were used to shed light on the mechanism of the established hydroxylation reaction of benzene with H2O2, catalyzed by [NiII(tepa)]2+ (tepa = tris[2-(pyridin-2-yl)ethyl]amine). Dinickel(III) bis(µ-oxo) species have been proposed as the key intermediate responsible for the benzene hydroxylation reaction. Our findings indicate that while the dinickel dioxygen species can be generated as a stable structure, it cannot serve as an active catalyst in this transformation. The calculations allowed us to unveil an unprecedented mechanism composed of six main steps as follows: (i) deprotonation of coordinated H2O2, (ii) oxidative addition, (iii) water elimination, (iv) benzene addition, (v) ketone generation, and (vi) tautomerization and regeneration of the active catalyst. Addition of benzene to oxygen, which occurs via a radical mechanism, turns out to be the rate-determining step in the overall reaction. This study demonstrates the critical role of Ni-oxyl species in such transformations, highlighting how the unpaired spin density value on oxygen and positive charges on the Ni-O• complex affect the activation barrier for benzene addition.

Ab initio molecular dynamics investigation of the Pt(111)-water interface structure in an alkaline environment with high surface OH-coverages.

In this study, we investigate the structure of the Pt(111)-water interface in an alkaline environment with large OH coverages of 1/3, 2/3 and 1 monolayer using a large well-equilibrated system. We observe that the OH coverage influences both the orientational distribution of the water molecules and their density, with more structure associated with higher coverage. At the same time, there is evidence of a highly dynamic hydrogen bond network on the lower coverage systems with substantial exchange of water between the surface and the solvent. In addition to OH and H2O species, which are preferentially located at the top sites, the 1/3 and 2/3 monolayer surfaces also contain O atoms, which are relatively stable and prefer the hollow sites. In contrast, the 1 monolayer surface shows none of these dynamics, and is unlikely to be active. The dynamic coexistence of O, OH and H2O on Pt(111) electrodes in alkaline conditions necessitates the investigation of several possible reaction paths for processess like ORR and water splitting. Finally, the exchange processes observed between the solvent and the interface underscore the need to explicitly include liquid water in simulations of systems similar to Pt(111).

Modelling the growth reaction pathways of zincone ALD/MLD hybrid thin films: a DFT study.

ALD/MLD hybrid thin films can be fabricated by combining atomic layer deposition (ALD) and molecular layer deposition (MLD). Even though this deposition method has been extensively used experimentally, the computational work required to acquire the reaction paths during the thin film deposition process is still in dire demand. We investigated hybrid thin films consisting of diethyl zinc and either 4-aminophenol or hydroquinone using both gas-phase and surface reactions to gain extensive knowledge of the complex phenomena occurring during the process of hybrid thin film deposition. We used density functional theory (DFT) to obtain the activation energies of these kinetic-dependent deposition processes. Different processes of ethyl ligand removal as ethane were discovered, and we found that the hydroxyl group of 4-aminophenol was more reactive than the amino group in the migration of hydrogen to an ethyl ligand within a complicated branching reaction chain.

Universal Trends between Acid Dissociation Constants in Protic and Aprotic Solvents.

pKa values in non-aqueous solvents are of critical importance in many areas of chemistry. Our knowledge is, despite their relevance, still limited to the most fundamental properties and few pKa values in the most common solvents. Taking advantage of a recently introduced computationally efficient procedure we computed the pKa values of 182 compounds in 21 solvents. This data set is used to establish for the first time universal trends between all solvents. Our computations indicate, that the total charge of the molecule and the charge of the acidic group combined with the Kamlet-Taft solvatochromic parameters are sufficient to predict pKa values with at least semi- quantitative accuracy. We find, that neutral acids such as alcohols are strongly affected by the solvent properties. This is contrasted by cationic acids like ammonium ions whose pKa is often almost completely independent from the choice of solvent.

Hydrogen adsorption trends on two metal-doped Ni2P surfaces for optimal catalyst design.

In this study, we looked at the hydrogen evolution reaction on the doubly doped Ni3P2 terminated Ni2P surface. Two Ni atoms in the first three layers of the Ni2P surface model were exchanged with two transition metal atoms. We limited our investigation to combinations of Al, Co, and Fe based on their individual effectiveness as Ni2P dopants in our previous computational studies. The DFT calculated hydrogen adsorption free energy was employed as a predictor of the materials' catalytic HER activity. Our results indicate that the combination of Co and Fe dopants most improves the catalytic activity of the surface through the creation of multiple novel and active catalytic sites.

From absolute potentials to a generalized computational standard hydrogen electrode for aqueous and non-aqueous solvents.

We describe a simple and efficient procedure to compute a conversion factor for the absolute potential of the standard hydrogen electrode in water to any other solvent. In contrast to earlier methods our procedure only requires the pKa of an arbitrary acid in water and few simple quantum chemical calculations as input. Thus, it is not affected adversely by experimental shortcomings related to measurements in non-aqueous solvents. By combining this conversion factor with the absolute potential in water, the absolute potential in the solvent of interest is obtained. Based on this procedure a new generalized computational standard hydrogen electrode for the computation of electron transfer and proton-coupled electron transfer potentials in non-aqueous solvents and ionic liquids is developed. This enables for the first time the reliable prediction of redox potentials in any solvent. The method is tested through calculation of absolute potentials in 36 solvents. Using the Kamlet-Taft linear solvation energy model we find that the relative absolute potentials consistently increase with decreasing polarisability and decreasing hydrogen bonding ability. For protic solvents good agreement with literature is observed while significant deviations are found for aprotic solvents. The obtained conversion factors are independent of the quantum chemical method, while minor differences are observed between solvation models. This does, however, not affect the global trends.

Revisiting the Volmer-Heyrovský mechanism of hydrogen evolution on a nitrogen doped carbon nanotube: constrained molecular dynamics versus the nudged elastic band method.

Density functional theory (DFT) based computational electrochemistry has the potential to serve as a tool with predictive power in the rational development and screening of electrocatalysts for renewable energy technologies. It is, however, of paramount importance that simulations are conducted rigorously at a level of theory that is sufficiently accurate in order to obtain physicochemically sensible results. Herein, we present a comparative study of the performance of the static climbing image nudged elastic band method (CI-NEB) vs. DFT based constrained molecular dynamics simulations with thermodynamic integration in estimating activation and reaction (free) energies of the Volmer-Heyrovský mechanism on a nitrogen doped carbon nanotube. Due to cancellation of errors within the CI-NEB calculations, static and dynamic activation barriers are observed to be surprisingly similar, while a substantial decrease in reaction energies is seen upon incorporation of solvent dynamics. This finding is attributed to two competing effects; (1) solvent reorganization that stabilizes the transition and, in particular, the product states with respect to the reactant state and (2) destabilizing entropic contributions due to solvent fluctuations. Our results highlight the importance of explicitly sampling the interfacial solvent dynamics when studying hydrogen evolution at solid-liquid interfaces.

Method for the accurate prediction of electron transfer potentials using an effective absolute potential.

A protocol for the accurate computation of electron transfer (ET) potentials from ab initio and density functional theory (DFT) calculations is described. The method relies on experimental pKa values, which can be measured accurately, to compute a computational setup dependent effective absolute potential. The effective absolute potentials calculated using this protocol display strong variations between the different computational setups and deviate in several cases significantly from the "generally accepted" value of 4.28 V. The most accurate estimate, obtained from CCSD(T)/aug-ccpvqz, indicates an absolute potential of 4.14 V for the normal hydrogen electrode (nhe) in water. Using the effective absolute potential in combination with CCSD(T) and a moderately sized basis, we are able to predict ET potentials accurately for a test set of small organic molecules (σ = 0.13 V). Similarly we find the effective absolute potential method to perform equally good or better for all considered DFT functionals compared to using one of the literature values for the absolute potential. For, M06-2X, which comprises the most accurate DFT method, standard deviation of 0.18 V is obtained. This improved performance is a result of using the most appropriate effective absolute potential for a given method.

Computational exploration of Fe55@C240-catalyzed Fischer-Tropsch synthesis.

Single-shell carbon-encapsulated iron nanoparticles (SCEINs), Fe@C, have been shown to be charge-transfer complexes that can act as effective catalysts in the hydrogen and oxygen evolution reactions. A new generation of Fe-based catalysts for Fischer-Tropsch synthesis (FTS) which resembles SCEINs, e.g. single carbide nanoparticles encapsulated in carbon shells, has demonstrated enhanced activity and stability for FTS as compared to other carbon-supported Fe-based FTS. Thus the catalytic ability of SCEINs for the water splitting reactions and the Fe-based FTS catalyst with SCEINs-like features stimulated our exploration of SCEINs-catalyzed FTS. We performed ab initio DFT calculations using a realistic SCEINs model system Fe55@C240 to investigate for the first time the adsorption of the main reactants in FTS (CO and H/H2) and further to evaluate the catalytic ability of Fe55@C240 by reproducing the key steps of the well-known Fe-based FTS mechanisms (carbide, enol and CO insertion). Our calculations revealed: (i) a determinant role of Fe in enhancing CO adsorption (ii) strong cooperative effects with the adsorbates that stabilize the binding (iii) a less favourable two-sites reaction on Fe55@C240 due to preferential positions of the reactants farther to each other which prevent enol and carbide FTS mechanisms. We propose therefore a possible CO insertion path for hydrocarbon growth on Fe55@C240.

Hydrogen adsorption trends on Al-doped Ni2P surfaces for optimal catalyst design.

Nanoparticles of nickel phosphide are promising materials to replace the currently used rare Pt-group metals at cathode-side electrodes in devices for electrochemical hydrogen production. Chemical modification by doping can be used to fine-tune the electrocatalytic activity, but this path requires theoretical, atomic-level support which has not been widely available for Ni-P. We present a density functional theory analysis of Al-doped Ni2P surfaces to identify structural motifs that could contribute to the improved behavior of the catalyst. Based on the formation energies of substitutionally Al-doped Ni sublattices, we find doping to take place preferably at the topmost layers. The Ni-Ni bridge and the P-top sites are the optimal ones in terms of hydrogen bonding energies. The Ni-Ni bridge site is not present on pristine surfaces but is a consequence of Al doping and provides a candidate to explain the experimentally observed high activities in doped Ni-P nanoparticles. Similar structural motifs can be recommended to be engineered for other Ni-P structures for improved electrocatalytic activity.

Hydrogen adsorption trends on various metal-doped Ni2P surfaces for optimal catalyst design.

In this study, we looked at the hydrogen evolution reaction on Mg-, Mo-, Fe-, Co-, V-, and Cu-doped Ni3P2 and Ni3P2 + P terminated Ni2P surfaces. The DFT calculated hydrogen adsorption free energy was employed as a predictor of the materials' catalytic HER activity. Our results indicate that doping can substantially improve the catalytic activity of the Ni3P2 terminated surface. In contrast, the Ni3P2 + P terminated one seems to be catalytically active irrespective of the type of doping, including in the absence of doping. Based on our doping energy and adsorption free energy calculations, the most promising dopants are iron and cobalt, whereas copper is less likely to function well as a doping element.

Diabatic model for electrochemical hydrogen evolution based on constrained DFT configuration interaction.

The accuracy of density functional theory (DFT) based kinetic models for electrocatalysis is diminished by spurious electron delocalization effects, which manifest as uncertainties in the predicted values of reaction and activation energies. In this work, we present a constrained DFT (CDFT) approach to alleviate overdelocalization effects in the Volmer-Heyrovsky mechanism of the hydrogen evolution reaction (HER). This method is applied a posteriori to configurations sampled along a reaction path to correct their relative stabilities. Concretely, the first step of this approach involves describing the reaction in terms of a set of diabatic states that are constructed by imposing suitable density constraints on the system. Refined reaction energy profiles are then recovered by performing a configuration interaction (CDFT-CI) calculation within the basis spanned by the diabatic states. After a careful validation of the proposed method, we examined HER catalysis on open-ended carbon nanotubes and discovered that CDFT-CI increased activation energies and decreased reaction energies relative to DFT predictions. We believe that a similar approach could also be adopted to treat overdelocalization effects in other electrocatalytic proton-coupled electron transfer reactions, e.g., in the oxygen reduction reaction.

Hydrogen adsorption on MoS2-surfaces: a DFT study on preferential sites and the effect of sulfur and hydrogen coverage.

We report a comprehensive computational study of the intricate structure-property relationships governing the hydrogen adsorption trends on MoS2 edges with varying S- and H-coverages, as well as provide insights into the role of individual adsorption sites. Additionally, the effect of single- and dual S-vacancies in the basal plane on the adsorption energetics is assessed, likewise with an emphasis on the H-coverage dependency. The employed edge/site-selective approach reveals significant variations in the adsorption free energies, ranging between ∼±1.0 eV for the different edges-types and S-saturations, including differences of even as much as ∼1.2 eV between sites on the same edge. The incrementally increasing hydrogen coverage is seen to mainly weaken the adsorption, but intriguingly for certain configurations a stabilizing effect is also observed. The strengthened binding is seen to be coupled with significant surface restructuring, most notably the splitting of terminal S2-dimers. Our work links the energetics of hydrogen adsorption on 2H-MoS2 to both static and dynamic geometrical features and quantifies the observed trends as a function of H-coverage, thus illustrating the complex structure/activity relationships of the MoS2 catalyst. The results of this systematical study aims to serve as guidance for experimentalists by suggesting feasible edge/S-coverage combinations, the synthesis of which would potentially yield the most optimally performing HER-catalysts.

Charge distribution and Fermi level in bimetallic nanoparticles.

Upon metal-metal contact, a transfer of electrons will occur between the metals until the Fermi levels in both phases are equal, resulting in a net charge difference across the metal-metal interface. Here, we have examined this contact electrification in bimetallic model systems composed of mixed Au-Ag nanoparticles containing ca. 600 atoms using density functional theory calculations. We present a new model to explain this charge transfer by considering the bimetallic system as a nanocapacitor with a potential difference equal to the work function difference, and with most of the transferred charge located directly at the contact interface. Identical results were obtained by considering surface contacts as well as by employing a continuum model, confirming that this model is general and can be applied to any multimetallic structure regardless of geometry or size (going from nano- to macroscale). Furthermore, the equilibrium Fermi level was found to be strongly dependent on the surface coverage of different metals, enabling the construction of scaling relations. We believe that the charge transfer due to Fermi level equilibration has a profound effect on the catalytic, electrocatalytic and other properties of bimetallic particles. Additionally, bimetallic nanoparticles are expected to have very interesting self-assembly for large superstructures due to the surface charge anisotropy between the two metals.

First principles study of the atomic layer deposition of alumina by TMA-H2O-process.

Atomic layer deposition (ALD) is a coating technology used to produce highly uniform thin films. Aluminiumoxide, Al2O3, is mainly deposited using trimethylaluminium (TMA) and water as precursors and is the most studied ALD-process to date. However, only few theoretical studies have been reported in the literature. The surface reaction mechanisms and energetics previously reported focus on a gibbsite-like surface model but a more realistic description of the surface can be achieved when the hydroxylation of the surface is taken into account using dissociatively adsorbed water molecules. The adsorbed water changes the structure of the surface and reaction energetics change considerably when compared to previously studied surface model. Here we have studied the TMA-H2O process using density functional theory on a hydroxylated alumina surface and reproduced the previous results for comparison. Mechanisms and energetics during both the TMA and the subsequent water pulse are presented. TMA is found to adsorb exothermically onto the surface. The reaction barriers for the ligand-exchange reactions between the TMA and the surface hydroxyl groups were found to be much lower compared to previously presented results. TMA dissociation on the surface is predicted to saturate at monomethylaluminium. Barriers for proton diffusion between surface sites are observed to be low. TMA adsorption was also found to be cooperative with the formation of methyl bridges between the adsorbants. The water pulse was studied using single water molecules reacting with the DMA and MMA surface species. Barriers for these reactions were found to reasonable in the process conditions. However, stabilizing interactions amongst water molecules were found to lower the reaction barriers and the dynamical nature of water is predicted to be of importance. It is expected that these calculations can only set an upper limit for the barriers during the water pulse.

CO oxidation catalyzed by neutral and anionic Cu20 clusters: relationship between charge and activity.

Reactions of CO and O2 on neutral and anionic Cu20 clusters have been investigated by spin-polarized density functional theory. Three reaction mechanisms of CO oxidation are explored: reactions with atomic oxygen (dissociated O2) as well as reactions with molecular oxygen, including Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) mechanisms. The adsorption energies, reaction pathways, and reaction barriers for CO oxidation are calculated systematically. The anionic Cu20(-) cluster can adsorb CO and O2 more strongly than the neutral counterpart due to the superatomic shell closing of 20 valence electrons which leaves one electron above the band gap. The activation of O2 molecule upon adsorption is crucial to determine the rate of CO oxidation. The CO oxidation proceeds efficiently on both Cu20 and Cu20(-) clusters, when O2 is pre-adsorbed dissociatively. The ER mechanism has a lower reaction barrier than the LH mechanism on the neutral Cu20. In general, CO oxidation occurs more readily on the anionic Cu20(-) (effective reaction barriers 0.1-0.3 eV) than on the neutral Cu20 cluster (0.3-0.5 eV). Moreover, Cu20(-) exhibits enhanced binding for CO2. From the analysis of the reverse direction of CO oxidation, it is observed that the transition of CO2 to CO + O can occur on the Cu20(-) cluster, which demonstrates that Cu clusters may serve as good catalyst for CO2 chemistry.