Biomimetic Photoexcited Cobaloxime Catalysis in Organic Synthesis.

Drawing inspiration from nature has long been a cornerstone of chemical innovation, with natural systems offering a wealth of untapped potential for discovery. In this minireview, we delve into the burgeoning field of cobaloxime catalysis in organic synthesis, which mimics the catalytic activity of the natural organometallic alkylcobalamine enzymes. Our focus lies on elucidating the latest advancements in this area, as well as delineating the primary mechanistic pathways at play. By describing, and comparing these mechanisms, we provide a comprehensive overview of the current state-of-the-art, while also shedding light on the key unresolved challenges that await further exploration.

Concluding remarks: Sustainable nitrogen activation - are we there yet?

The activation of dinitrogen as a fundamental step in reactions to produce nitrogen compounds, including ammonia and nitrates, has a cornerstone role in chemistry. Bringing together research from disparate fields where this can be achieved sustainably, this Faraday Discussion seeks to build connections between approaches that can stimulate further advances. In this paper we set out to provide an overview of these different approaches and their commonalities. We explore experimental aspects including the positive role of increasing nitrogen pressure in some fields, as well as offering perspectives on when 15N2 experiments might, and might not, be necessary. Deconstructing the nitrogen reduction reaction, we attempt to provide a common framework of energetic scales within which all of the different approaches and their components can be understood. On sustainability, we argue that although green ammonia produced from a green-H2-fed Haber-Bosch process seems to fit the bill, there remain many real-world contexts in which other, sustainable, approaches to this vital reaction are urgently needed.

Stabilisation of dianion dimers trapped inside cyanostar macrocycles.

Interanionic H-bonds (IAHBs) are unfavourable interactions in the gas phase becoming favoured when anions are in solution. Dianion dimers are also susceptible to being trapped inside the cavities of cyanostar (CS) macrocycles, and thus, the formation of 2 : 2 anion : cyanostar aggregates is mainly supported by three kinds of interactions: IAHBs between the dianions, π-π stacking between the confronted cyanostars, and the presence of an intricate network of multiple C(sp2)HO H-bonds between cyanostar ligands and the anionic moieties. An analysis of the interaction energies supported by NBO reveals a slight cooperative effect of the CSs on the IAHB stabilisation.

A Conceptual DFT Study of Phosphonate Dimers: Dianions Supported by H-Bonds.

A conceptual DFT study of the dissociation of anionic and neutral phosphonate dimers has been carried out. In addition, the dianion complexes have been studied in the presence of two solvents, water and tetrahydrofuran. The dissociation of the dianion complexes in the gas phase and in solution present a maximum along the reaction coordinate that is not present in the neutral-neutral and anion-neutral complexes. The principal chemical descriptors (chemical potential, reaction electronic flux, hardness, and global electrophilicity index) do not show changes in their trends along the dissociation profiles even when there is an energy maximum in the case of the anion-anion complexes.

Hydrogen Evolution in [NiFe] Hydrogenases: A Case of Heterolytic Approach between Proton and Hydride.

The mechanism for Hydrogen Evolution Reaction (HER) in [NiFe] hydrogenase enzymes distinguishes them from inorganic catalysts. The first H+/e- pair injected to the active site of the hydrogenases transforms into hydride, while the second H+/e- pair injection leads to the formation of the H-/H+ pair both binding to the active site. The two opposite charged hydrogens heterolytically approach each other in order to form dihydrogen (H2), which is enhanced by the Coulomb force. Two previously proposed reaction routes for this process have been examined by Conceptual Density Functional Theory (DFT) in this work. One presents better agreement with experimental spectra, while the other is thermodynamically more favorable. Both paths suggest that the approach and the charge transfer between the proton and hydride are motivated by the stabilization of the electronic activity and the electrophilicity of Ni. After the heterolytic approach of the proton and hydride moieties, the two hydrogen atoms attach to the Ni ion and combine homolytically.

Quantifying electronic similarities between NHC-gold(i) complexes and their isolobal imidazolium precursors.

A series of NHC-gold(i) (NHC = N-heterocyclic carbene) complexes has been studied by DFT calculations, enabling comparison of electronic and NMR behaviour with related protonated and free NHC molecules. Based on calculations, the NMR resonances of the carbenic C2 carbon atom in [Au(NHC)(Cl)] and [NHC(H)][Cl] exhibit increased shielding when compared to the free N-heterocyclic carbenes by an average of 46.6 ± 2.2 and 73.7 ± 4.3 ppm, respectively. A similar trend is observed when analysing the paramagnetic term of the magnetic shielding tensor. Although gold(i) and proton are considered isolobal fragments, imidazolium compounds lack π-backdonation due to the energetic unavailability of d-orbitals in H+. We propose that NHC-gold(i) complexes exhibit important π-backdonation irrespective of the relative amount of σ-donation between the NHC and gold(i)-X (X = anionic ligand) moieties in Au-NHC complexes. Interestingly, a correlation exists between the calculated shielding for gold (197Au) and the π-donation and π-backdonation contributions. We describe that this correlation also exists when analysing the σ-backdonation term, a property generally ignored yet representing a significant energetic contribution to the stability of the C2-Au bond.

Hydrogen bonding effect between active site and protein environment on catalysis performance in H2-producing [NiFe] hydrogenases.

The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H+/e- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S- in cysteine residues shows a larger influence on the energy profile than that via the CN- ligands on Fe.

Hydrogenation of CO2 -Derived Carbonates and Polycarbonates to Methanol and Diols by Metal-Ligand Cooperative Manganese Catalysis.

The first base-metal-catalysed hydrogenation of CO2 -derived carbonates to alcohols is presented. The reaction proceeds under mild conditions in the presence of a well-defined manganese complex with a loading as low as 0.25 mol %. The non-precious-metal homogenous catalytic system provides an indirect route for the conversion of CO2 into methanol with the co-production of value-added (vicinal) diols in yields of up to 99 %. Experimental and computational studies indicate a metal-ligand cooperative catalysis mechanism.

Single-Site Molybdenum on Solid Support Materials for Catalytic Hydrogenation of N2 -into-NH3.

Very stable in operando and low-loaded atomic molybdenum on solid-support materials have been prepared and tested to be catalytically active for N2 -into-NH3 hydrogenation. Ammonia synthesis is reported at atmospheric pressure and 400 °C with NH3 rates of approximately 1.3×103  µmol h-1 gMo -1 using a well-defined Mo-hydride grafted on silica (SiO2-700 ). DFT modelling on the reaction mechanism suggests that N2 spontaneously binds on monopodal [(≡Si-O-)MoH3 ]. Based on calculations, the fourth hydrogenation step involving the release of the first NH3 molecule represents the rate-limiting step of the whole reaction. The inclusion of cobalt co-catalyst and an alkali caesium additive impregnated on a mesoporous SBA-15 support increases the formation of NH3 with rates of circa 3.5×103  µmol h-1 gMo -1 under similar operating conditions and maximum yield of 29×103  µmol h-1 gMo -1 when the pressure is increased to 30 atm.

Feasibility of N2 Binding and Reduction to Ammonia on Fe-Deposited MoS2 2D Sheets: A DFT Study.

Based on the structure of the nitrogenase FeMo cofactor (FeMoco), it is reported that Fe deposited on MoS2 2D sheets exhibits high selectivity towards the spontaneous fixation of N2 against chemisorption of CO2 and H2 O. DFT predictions also indicate the ability of this material to convert N2 into NH3 with a maximum energy input of 1.02 eV as an activation barrier for the first proton-electron pair transfer.

Experimental and Computational Study of an Unexpected Iron-Catalyzed Carboetherification by Cooperative Metal and Ligand Substrate Interaction and Proton Shuttling.

An iron-catalyzed cycloisomerization of allenols to deoxygenated pyranose glycals has been developed. Combined experimental and computational studies show that the iron complex exhibits a dual catalytic role in that the non-innocent cyclopentadienone ligand acts as proton shuttle by initial hydrogen abstraction from the alcohol and by facilitating protonation and deprotonation events in the isomerization and demetalation steps. Molecular orbital analysis provides insight into the unexpected and selective formation of the 3,4-dihydro-2H-pyran.

Why is a proton transformed into a hydride by [NiFe] hydrogenases? An intrinsic reactivity analysis based on conceptual DFT.

The hydrogen evolution reaction (HER) catalysed by [NiFe] hydrogenases entails a series of chemical events involving great mechanistic interest. In an attempt to understand and delve into the question about 'Why does nature work in that way?', an in-depth intrinsic reactivity analysis based on conceptual DFT has been carried out focusing on the so-called to step, i.e. our work tries to answer how and why the proton attached to the reactive sulphur atom from one of the exo-cyclic cysteine residues is transformed into a bridging hydride to be shared between the Ni/Fe metals in the active site of [NiFe] hydrogenases, which involves not only H migration, but also a change of the charge state on Ni from Ni(i) to Ni(iii). Our DFT results suggest that the transformation is motivated by spontaneous rearrangements of the electron density, and stabilisation comes from the decrease of both electronic activity and electrophilicity index from Ni.

A DFT study of planar vs. corrugated graphene-like carbon nitride (g-C3N4) and its role in the catalytic performance of CO2 conversion.

Graphene-like carbon nitride (g-C3N4), a metal-free 2D material that is of interest as a CO2 reduction catalyst, is stabilised by corrugation in order to minimise the electronic repulsions experienced by the N lone pairs located in their structural holes. This conformational change not only stabilises the Fermi level in comparison with the totally planar structure, but also increases the potential depth of the π-holes, representing the active sites where the catalytic CO2 conversion takes place. Finally, as a result of corrugation, our DFT-D3 calculations indicate that the reaction Gibbs free energy for the first H(+)/e(-) addition decreases by 0.49 eV with respect to the totally planar case, suggesting that corrugation not only involves the material's stabilisation but also enhances the catalytic performance for the selective production of CO/CH3OH.

Understanding the Aldo-Enediolate Tautomerism of Glycolaldehyde in Basic Aqueous Solutions.

The biochemically important interconversion process between aldoses and ketoses is assumed to take place via 1,2-enediol or 1,2-enediolate intermediates, but such intermediates have never been isolated. The current work was undertaken in an attempt to detect the presence of the 1,2-enediol structure of glycolaldehyde in alkaline medium, actually a 1,2-enediolate, and to try to clarify the scarce data existing about both the formation of deprotonated enediol and the aldo-enediolate equilibrium. The Raman spectra of neutral and basic solutions were recorded as a function of time for eleven days. Several bands associated with the presence of the enediolate were observed in alkaline medium. Glycolaldehyde exists as three different structures in aqueous solution at neutral pH, that is, hydrated aldehydes, aldehydes and dimers, with a respective ratio of approximately 4:0.25:1. Additionally, the formation of Z-enediolate forms takes place at basic pH, together with an increase in the concentration of aldehyde species, such as 2-oxoethan-1-olate, and a decrease in the concentrations of the hydrated aldehyde and dimeric forms. The theoretical ratio of ≈1.5:1 for aldehyde:Z-enediolate reproduces the experimental Raman spectrum in basic medium, with an additional contribution of the previously mentioned ratio between the hydrated aldehyde and dimeric forms. Finally, Raman spectroscopy allowed us to monitor the enolization of this carbohydrate model and conclude that aldo-enediol tautomerism-formally aldo-enediolate-happens when a suitable amount of basic species is added.

Chalcogen bonds in complexes of SOXY (X, Y = F, Cl) with nitrogen bases.

SOF2, SOFCl, and SOCl2 were each paired with a series of N bases. The potential energy surface of the binary complexes were characterized by MP2 calculations with double and triple-ξ basis sets, extrapolated to complete sets. The most stable configurations contained a S···N chalcogen bond with interaction energies as high as 6.8 kcal/mol. These structures are stabilized by a Nlp → σ*(S-Z) electron transfer (Z = O, F, Cl), complemented by Coulombic attraction of N to the σ-hole opposite the Z atom. N···S-F and N···S-Cl chalcogen bonds are stronger than N···S═O interactions. Formation of each chalcogen bond elongates all of the internal covalent bonds within SOXY, especially the S-Cl bond. Halogen-bonded (N···Cl-S) complexes were also observed, but these are more weakly bound, by less than 3 kcal/mol.

Complexes containing CO2 and SO2. Mixed dimers, trimers and tetramers.

Mixed dimers, trimers and tetramers composed of SO2 and CO2 molecules are examined by ab initio calculations to identify all minimum energy structures. While AIM formalism leads to the idea of a pair of C···O bonds in the most stable heterodimer, bound by some 2 kcal mol(-1), NBO analysis describes the bonding in terms of charge transfer from O lone pairs of SO2 to the CO π* antibonding orbitals. The second minimum on the surface, just slightly less stable, is described by AIM as containing a single O···O chalcogen bond. The NBO picture is that of two transfers in opposite directions: one from a SO2 O lone pair to a π* antibond of CO2, supplemented by CO2 Olp → π*(SO). Decomposition of the interaction energies points to electrostatic attraction and dispersion as the dominant attractive components, in roughly equal measure. The various heterotrimers and tetramers generally retain the dimer structure as a starting point. Cyclic oligomers are favored over linear geometries, with a preference for complexes containing larger numbers of SO2 molecules.

Strongly bound noncovalent (SO3)n:H2CO complexes (n = 1, 2).

The potential energy surfaces (PES) for the SO3:H2CO and (SO3)2:H2CO complexes were thoroughly examined at the MP2/aug-cc-pVDZ computational level. Heterodimers and trimers are held together primarily by SO chalcogen bonds, supplemented by weaker CHO and/or OC bonds. The nature of the interactions is probed by a variety of means, including electrostatic potentials, AIM, NBO, energy decomposition, and electron density redistribution maps. The most stable dimer is strongly bound, with an interaction energy exceeding 10 kcal mol(-1). Trimers adopt the geometry of the most stable dimer, with an added SO3 molecule situated so as to interact with both of the original molecules. The trimers are strongly bound, with total interaction energies of more than 20 kcal mol(-1). Most such trimers show positive cooperativity, with shorter SO distances, and three-body interaction energies of nearly 3 kcal mol(-1).

Substituent Effects in the Noncovalent Bonding of SO2 to Molecules Containing a Carbonyl Group. The Dominating Role of the Chalcogen Bond.

The SO2 molecule is paired with a number of carbonyl-containing molecules, and the properties of the resulting complexes are calculated by high-level ab initio theory. The global minimum of each pair is held together primarily by a S···O chalcogen bond wherein the lone pairs of the carbonyl O transfer charge to the π* antibonding SO orbital, supplemented by smaller contributions from weak CH···O H-bonds. The binding energies vary between 4.2 and 8.6 kcal/mol, competitive with even some of the stronger noncovalent forces such as H-bonds and halogen bonds. The geometrical arrangement places the carbonyl O atom above the plane of the SO2 molecule, consistent with the disposition of the molecular electrostatic potentials of the two monomers. This S···O bond differs from the more commonly observed chalcogen bond in both geometry and origin. Substituents exert their influence via inductive effects that change the availability of the carbonyl O lone pairs as well as the intensity of the negative electrostatic potential surrounding this atom.

Complexation of n SO2 molecules (n = 1, 2, 3) with formaldehyde and thioformaldehyde.

Ab initio and density functional theory calculations are used to examine complexes formed between H2CO and H2CS with 1, 2, and 3 molecules of SO2. The nature of the interactions is probed by a variety of means, including electrostatic potentials, natural bond orbital, atoms in molecules, energy decomposition, and electron density redistribution maps. The dimers are relatively strongly bound, with interaction energies exceeding 5 kcal/mol. The structures are cyclic, containing both a O/S⋯S chalcogen bond and a CH⋯O H-bond. Addition of a second SO2 molecule leads to a variety of heterotrimer structures, most of which resemble the original dimer, where the second SO2 molecule engages in a chalcogen bond with the first SO2, and a C⋯O attraction with the H2CX. Some cooperativity is apparent in the trimers and tetramers, with an attractive three-body interaction energy and shortened intermolecular distances.

An exploration of the ozone dimer potential energy surface.

The (O3)2 dimer potential energy surface is thoroughly explored at the ab initio CCSD(T) computational level. Five minima are characterized with binding energies between 0.35 and 2.24 kcal/mol. The most stable may be characterized as slipped parallel, with the two O3 monomers situated in parallel planes. Partitioning of the interaction energy points to dispersion and exchange as the prime contributors to the stability, with varying contributions from electrostatic energy, which is repulsive in one case. Atoms in Molecules analysis of the wavefunction presents specific O⋯O bonding interactions, whose number is related to the overall stability of each dimer. All internal vibrational frequencies are shifted to the red by dimerization, particularly the antisymmetric stretching mode whose shift is as high as 111 cm(-1). In addition to the five minima, 11 higher-order stationary points are identified.