N2 Bond Research Articles

Conformational transformations of N-[arylsulfonylimino(methyl)methyl]-1,4-benzoquinonemonoimines

In solutions of N-[arylsulfonylimino(methyl)methyl]-1,4-benzoquinonemonoimines, several Z,E-isomerization processes occur simultaneously: a fast inversion isomerization process (topomerization) relative to the C=N1 quinoneimine bond; a slow inversion isomerization process relative to the C=N2 exocyclic bond; and a process of inhibited rotation around the =N1–С= bond that connects two imine fragments. In the 1H NMR spectra, isomerization with respect to the C=N2 exocyclic bond was only detected. The theoretical values of the isomerization barriers for the C=N1 and C=N2 bonds, as well as for the inhibited rotation around the =N1–C= bond were determined, along with the experimental Z,E-isomerization barriers for to the C=N2 bond. It was established that in the 1H NMR spectra for N-[arylsulfonylimino(methyl)methyl]-1,4-benzoquinonemonoimines, only Z, E-isomerization relative to the C=N2 bond can be observed.

Deciphering the Nitrogen Activation Mechanisms on Group VIII Single Atoms at MoS2.

The activation of nitrogen (N2) is vital for sustainable ammonia production and nitrogen fixation technologies. This study employs density functional theory (DFT) to investigate the nitrogen activation and reduction capabilities of Group VIII single-atom catalysts anchored on MoS2. Among these, osmium anchored on MoS2 (Os@MoS2) emerged as the most promising catalyst, exhibiting the highest N2 activation and the lowest nitrogen reduction reaction (NRR) overpotential (0.624 V). A pronounced "electron drift" effect was observed for Os@MoS2, leading to significant charge redistribution that weakens the N ≡ N triple bond, facilitating its activation. The N-N dissociation energy barrier at the *N-NH2 intermediate was calculated to be only 0.82 eV, confirming Os@MoS2's superior catalytic efficiency. Detailed analyses, including electrostatic potential maps, electron localization functions, spin density, and charge transfer, revealed the pivotal role of orbital interactions in driving N2 activation. Interestingly, the trends in adsorbed N2 bond energies and NRR overpotentials showed a consistent diagonal pattern across the Group VIII catalysts, emphasizing the importance of electronic and geometric factors. This work offers valuable insights into nitrogen activation mechanisms and provides a framework for designing efficient catalysts, highlighting Os@MoS2's potential in sustainable ammonia synthesis.

Final E5 to E8 Steps in the Nitrogenase Mechanism for Nitrogen Fixation.

Nitrogenase converts nitrogen in the air to ammonia. It is often regarded as the second most important enzyme in nature after photosystem II. The mechanism for how nitrogenase is able to perform the difficult task of cleaving the strong bond in N2 is debated. It is known that for every electron that is donated to N2, two ATP are hydrolyzed. In the experimentally suggested mechanism, the activation occurs after four reductions of the ground state, but there is no suggestion for how the enzyme uses the hydrolysis energy to perform catalysis. In the theoretical mechanism, it is suggested that hydrolysis is used to reduce the electron donor. In previous papers, the steps leading to the activation of N2 in the so-called E4 state has been investigated, using both the experimental and theoretical mechanism, showing that only the theoretical one leads to agreement with EPR observations for E4. In the present paper, the four steps following E4, leading to the release of two ammonia molecules, are described using the same methodology as used in the previous studies.

Activity of CeO2(111) supported Mnx(x = 1–6) for electrochemical N2 reduction reaction: Insights from density functional theory

It is challenging to realize an efficient nitrogen reduction reaction (NRR) under mild conditions, which suffers from low ammonia yield and low Faraday efficiency due to the extremely stable NN triple bond of N2 as well as competitive hydrogen evolution reaction. In this work, the NRR reactivity of Mnx(x= 1–6) clusters supported on CeO2(111) (Mnx(x = 1–6)/CeO2(111)) was systematically investigated using density functional theory. A volcanic relationship between the limiting potential of NRR on Mnx(x = 1–6)/CeO2(111) and the atom number of Mnx was found. Mn3/CeO2(111) shows the highest activity for NRR with a limiting potential of -0.36, -0.55 and -0.53 V along distal, alternating and enzymatic reaction pathway, respectively. Its high activity is attributed to the triangular geometry and optimal average number of electrons every Mn transferred from Mn3 to CeO2(111), which leads to the strong N2 activation and the stabilization of nitrogen-containing intermediates. Also, Mn3/CeO2(111) exhibits a high NRR selectivity by hindering H adsorption and a high thermal stability at both 298 and 773 K, suggesting its promising potential as effective NRR catalyst. This work provides new insights into the rational design of single cluster catalysts.

Sulfide release and rebinding in the mechanism for nitrogenase.

Nitrogenases are the only enzymes that activate the strong triple bond in N2. The mechanism for the activation has been very difficult to determine in spite of decades of work. In previous modeling studies it has been suggested that the mechanism for nitrogen activation starts out by four pre-activation steps (A0-A4) before catalysis. That suggestion led to excellent agreement with experimental Elecrtron Paramagnetic Resonance (EPR) observations in the step where N2 becomes protonated (E4). An important part of the pre-activation is that a sulfide is released. In the present paper, the details of the pre-activation are modeled, including the release of the sulfide. Several possible transition states for the release have been obtained. An A4(E0) state is reached which is very similar to the E4 state. For completeness, the steps going back from A4(E0) to A0 after catalysis are also modeled, including the insertion of a sulfide.

Exploring dynamics in single atom catalyst research: A comprehensive DFT-kMC study of nitrogen reduction reaction with focus on TM aggregation

Transition metal (TM)-based single atom catalysts (SACs) have emerged as a promising solution for the electrochemical nitrogen reduction reaction (NRR) due to their unique d-orbitals and low coordination, which facilitate efficient N2 bond weakening and increased reactivity. However, conventional theoretical approaches fail to incorporate the kinetics of TM aggregation, a crucial factor in SAC systems that affects catalyst stability and activity over time. Addressing this, we combined kinetic evaluations with thermodynamic and electronic analyses on 216 SAC candidates, using boron carbon nitride (BCN) as a substrate. Our findings revealed that Cr anchored to BCN exhibits a high turnover frequency (3.1 × 10−6 s−1) under mild conditions (300 K, 1 bar), attributed to its minimal TM aggregation over time. This research not only fills the gap in kinetic research within the SAC field but also proposes a SAC candidate distinguished by its superior kinetics, thermodynamic, and electronic properties.

Theoretical Insights into N-Glycoside Bond Cleavage of 5-Carboxycytosine by Thymine DNA Glycosylase: A QM/MM Study.

Thymine DNA glycosylase (TDG)-mediated excision of 5-formylcytosine and 5-carboxylcytosine (5-caC) is a critical step in active DNA demethylation. Herein, we employed a combined quantum mechanics/molecular mechanics approach to investigate the reaction mechanism of TDG-catalyzed N-glycosidic bond cleavage of 5-caC. The calculated results show that TDG-catalyzed 5-caC excision follows a concerted (SN2) mechanism in which glycosidic bond dissociation is coupled with nucleophile attack. Protonation of the 5-caC anion contributes to the cleavage of the N-glycoside bond, in which the N3-protonated zwitterion and imino tautomers are more favorable than carboxyl-protonated amino tautomers. This is consistent with the experimental data. Furthermore, our results reveal that the configuration rearrangement process of the protonated 5-caC would lower the stability of the N-glycoside bond and substantially reduce the barrier height for the subsequent C1'-N1 bond cleavage. This should be attributed to the smaller electrostatic repulsion between the leaving base and the negative phosphate group as a result of the structural rearrangement.

Iron-catalyzed azidation of cyclobutanol by C[sbnd]C bond cleavage

An iron-catalyzed ring-opening azidation of cyclobutanols through CC bond cleavage is disclosed for the first time. A range of tertiary cyclobutanols are tolerated providing γ-carbonyl-containing secondary alkyl azides. The method provides a straightforward approach for introducing azide group into aliphatic chain of organic molecules. A plausible mechanism involving radical-mediated sequential CC bond cleavage and C−N3 bond formation is proposed.

Nitrogen-Rich Molybdenum Nitride Synthesized in a Crucible under Air.

The triple bond in N2 is significantly stronger than the double bond in O2, meaning that synthesizing nitrogen-rich nitrides typically requires activated nitrogen precursors, such as ammonia, plasma-cracked atomic nitrogen, or high-pressure N2. Here, we report a synthesis of nitrogen-rich nitrides under ambient pressure and atmosphere. Using Na2MoO4 and dicyandiamide precursors, we synthesized nitrogen-rich γ-Mo2N3 in an alumina crucible under an ambient atmosphere, heated in a box furnace between 500 and 600 °C. Byproducts of this metathesis reaction include volatile gases and solid Na(OCN), which can be washed away with water. X-ray diffraction and neutron diffraction showed Mo2N3 with a rock salt structure having cation vacancies, with no oxygen incorporation, in contrast to the more common nitrogen-poor rock salt Mo2N with anion vacancies. Moreover, an increase in temperature to 700 °C resulted in molybdenum oxynitride, Mo0.84N0.72O0.27. This work illustrates the potential for dicyandiamide as an ambient-temperature metathesis precursor for an increased effective nitrogen chemical potential under ambient conditions. The classical experimental setting often used for solid-state oxide synthesis, therefore, has the potential to expand the nitride chemistry.

Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH)4H]- Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods.

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.

Bond dissociation energy of N2 measured by state-to-state resolved threshold fragment yield spectra.

The precise determination of the bond dissociation energy of N2 is crucial for thermochemistry database and theoretical calculations. However, there has been ongoing debate regarding its exact value. In this study, we used the velocity map imaging method combined with an extreme ultraviolet laser to measure the threshold fragment yield (TFY) spectra of N2 in the N(2D) + N(2D) photodissociation channels. By integrating the signals within a small circular area on the fragment velocity map images, we were able to obtain TFY spectra at nine different dissociation thresholds. These spectra are rotational state-resolved for the N2(J″) molecules and spin-orbit state-resolved for the dissociation channels involving N(2D) fragments. By employing the Wigner threshold law to simulate the TFY spectra and conducting statistical analysis on the comprehensive dataset, we determined the N2 bond dissociation energy to be 78 691.09 ± 0.15cm-1. This work now places N2 among the few diatomic molecules with bond dissociation energies measured at sub-wavenumber precision.

H2 formation from the E2-E4 states of nitrogenase.

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 ×102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).

Experimental and kinetic modeling study on the low-temperature decomposition and autoignition of 2-Azido-N,N-dimethylethanamine: A promising green mono- and bi-propellant

2-Azido-N,N-Dimethylethanamine (DMAZ) is a promising candidate for mono- and bi- propellant. However, the fundamental gas-phase combustion experiments have not been reported, and its chemical kinetic mechanism is not well understood. Therefore, the ignition delay times (IDTs) of DMAZ were measured utilizing a rapid compression machine and a shock tube, covering a wide temperature range of 570 – 960 K, at 10 and 20 bar with varying equivalence ratios. DMAZ was surprisingly found to undergo decomposition at temperatures as low as 600 K, leading to a pressure rise within the chamber. The low-temperature decomposition characteristics of DMAZ were systematically investigated under various fuel concentrations and pressures. A kinetic model of DMAZ was developed, incorporating quantum chemistry calculations for the thermodynamic data of new species and the rate constants of H-atom abstractions. The newly measured IDTs and characteristic decomposition times (CDTs) were further adopted in the model validation. Results show that DMAZ mainly decomposes through N–N2 bond fissions, which are also the major reaction pathways during autoignition. O2 addition to the radicals derived from the decomposition products and subsequent reactions contribute to most of the low-temperature reactivity in DMAZ oxidation. The current kinetic model can reasonably predict the IDTs and CDTs, as well as their dependencies on pressure, equivalence ratio, and fuel concentration.

Contrasting potential for biological N2 fixation at three polluted central European Sphagnum peat bogs: combining the 15N2-tracer and natural-abundance isotope approaches

Abstract. Availability of reactive nitrogen (Nr) is a key control on carbon (C) sequestration in wetlands. To complement the metabolic demands of Sphagnum in pristine rain-fed bogs, diazotrophs supply additional Nr via biological nitrogen fixation (BNF). As breaking the triple bond of atmospheric N2 is energy-intensive, it is reasonable to assume that increasing inputs of pollutant Nr will lead to BNF downregulation. However, recent studies have also documented measurable BNF rates in Sphagnum-dominated bogs in polluted regions, indicating the adaptation of N2 fixers to changing N deposition. Our aim was to quantify BNF in high-elevation peatlands located in industrialized central Europe. A 15N2-tracer experiment was combined with a natural-abundance N-isotope study at three Sphagnum-dominated peat bogs in the northern Czech Republic in an attempt to assess the roles of individual BNF drivers. High short-term BNF rates (8.2 ± 4.6 g N m2 d−1) were observed at Malé mechové jezírko, which receives ∼ 17 kg Nr ha−1 yr−1. The remaining two peat bogs, whose recent atmospheric Nr inputs differed from Malé mechové jezírko by only 1–2 kg ha−1 yr−1 (Uhlír̆ská and Brumiště), showed zero BNF. The following parameters were investigated to elucidate the BNF difference: the NH4+-N / NO3--N ratio, temperature, wetness, Sphagnum species, organic-N availability, possible P limitation, possible molybdenum (Mo) limitation, SO42- deposition, and pH. At Malé mechové jezírko and Uhlír̆ská, the same moss species (S. girgensohnii) was used for the 15N2 experiment; therefore, the host identity could not explain the difference in BNF at these sites. Temperature and moisture were also identical in all incubations and could not explain the between-site differences in BNF. The N : P stoichiometry in peat and bog water indicated that Brumiště may have lacked BNF due to P limitation, whereas non-detectable BNF at Uhlír̆ská may have been related to the 70-fold higher SO42- concentration in bog water. Across the sites, the mean natural-abundance δ15N values increased in the following order: atmospheric deposition (−5.3 ± 0.3 ‰) < Sphagnum (−4.3 ± 0.1 ‰) < bog water (−3.9 ± 0.4 ‰) < atmospheric N2 (0.0 ‰). Only at Brumiště was N in Sphagnum significantly isotopically heavier than in atmospheric deposition, possibly indicating a longer-term BNF effect. Collectively, our data highlight spatial heterogeneity in BNF rates under high Nr inputs as well as the importance of environmental parameters other than atmospheric Nr pollution in regulating BNF.

Observation of Continuum State Dissociation Enables the Determination of N2 Bond Dissociation Energy to Spectroscopic Accuracy.

Nitrogen (N) is one of the most fundamental elements of life. Precise determination of the bond dissociation energy (BDE) of the corresponding homonuclear diatomic molecule N2 is not only important for calculating the enthalpies of formation for any N-containing molecules but also provides the best benchmark for evaluating theoretical computational methods. Thus, it has attracted many experimental and theoretical studies, while controversies still exist. Here, we report the observation of continuum state dissociation of N2 into the channel N(2D5/2,3/2)+N(2D5/2,3/2) for the first time by using the vacuum ultraviolet (VUV)-pump-VUV-probe time-sliced velocity-mapped imaging setup. The quantum-state-resolved images enable the direct visualization of the dissociation onsets corresponding to each of the correlated spin-orbit fine-structure channels within a few tenths of wavenumber. The BDEs of 14N2 and 15N2 are directly determined to be 78691.8 ± 0.3 cm-1 and 78731.5 ± 0.3 cm-1, respectively, which should represent the most accurate BDE of N2 thus far.

Protonation of Homocitrate and the E1 State of Fe-Nitrogenase Studied by QM/MM Calculations.

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological life. There are three isozymes of nitrogenase, differing in the composition of the active site, viz., Mo, V, and Fe-nitrogenase. Recently, the first crystal structure of Fe-nitrogenase was presented. We have performed the first combined quantum mechanical and molecular mechanical (QM/MM) study of Fe-nitrogenase. We show with QM/MM and quantum-refinement calculations that the homocitrate ligand is most likely protonated on the alcohol oxygen in the resting E0 state. The most stable broken-symmetry (BS) states are the same as for Mo-nitrogenase, i.e., the three Noodleman BS7-type states (with a surplus of β spin on the eighth Fe ion), which maximize the number of nearby antiferromagnetically coupled Fe-Fe pairs. For the E1 state, we find that protonation of the S2B μ2 belt sulfide ion is most favorable, 14-117 kJ/mol more stable than structures with a Fe-bound hydride ion (the best has a hydride ion on the Fe2 ion) calculated with four different density-functional theory methods. This is similar to what was found for Mo-nitrogenase, but it does not explain the recent EPR observation that the E1 state of Fe-nitrogenase should contain a photolyzable hydride ion. For the E1 state, many BS states are close in energy, and the preferred BS state differs depending on the position of the extra proton and which density functional is used.

Mössbauer and Nuclear Resonance Vibrational Spectroscopy Studies of Iron Species Involved in N-N Bond Cleavage.

Diketiminate-supported iron complexes are capable of cleaving the strong triple bond of N2 to give a tetra-iron complex with two nitrides (Rodriguez et al., Science, 2011, 334, 780-783). The mechanism of this reaction has been difficult to determine, but a transient green species was observed during the reaction that corresponds to a potential intermediate. Here, we describe studies aiming to identify the characteristics of this intermediate, using a range of spectroscopic techniques, including Mössbauer spectroscopy, electronic absorption spectroscopy, Raman spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, and nuclear resonance vibrational spectroscopy (NRVS) complemented by density functional theory (DFT) calculations. We successfully elucidated the nature of the starting iron(II) species and the bis(nitride) species in THF solution, and in each case, THF breaks up the multiiron species. Various observations on the green intermediate species indicate that it has one N2 per two Fe atoms, has THF associated with it, and has NRVS features indicative of bridging N2. Computational models with a formally diiron(0)-N2 core are most consistent with the accumulated data, and on this basis, a mechanism for N2 splitting is suggested. This work shows the power of combining NRVS, Mössbauer, NMR, and vibrational spectroscopies with computations for revealing the nature of transient iron species during N2 cleavage.

Molecular structure of 3,4-dicyanofuroxan studied by gas electron diffraction. Application of DFT theory and coupled cluster computations

Equilibrium molecular structure of 3,4-dicyanofuroxan has been determined for the first time by means of gas electron diffraction (GED) and quantum chemical (QC) calculations. The GED values have been compared to those obtained from QC computations at different levels of theory. It has been discovered that the value of the O1−N2 bond is very sensitive to the QC method applied. The best fit to the experimental data has performed CCSD(T)/cc-pVTZ, whereas DFT methods applied in this work have been found to be less accurate. Peculiar features of 3,4-dicyanofuroxan molecular structure have been revealed by applying NBO, AIM and NCI analyses.

D-Orbital Electron Delocalization Realized by Axial Fe4 C Atomic Clusters Delivers High-Performance Fe-N-C Catalysts for Oxygen Reduction Reaction.

Fe-N-C catalyst for oxygen reduction reaction(ORR) has been considered as the most promising nonprecious metal catalyst due to its comparable catalytic performance to Pt in proton exchange membrane fuel cells (PEMFCs). The active centers of Fe-pyrrolic N4 have been proven to be extremely active for ORR. However, forming a stable Fe-pyrrolic N4 structure is a huge challenge. Here, a Cyan-Fe-N-C catalyst with Fe-pyrrolic N4 as the intrinsic active center is constructed with the help of axial Fe4 C atomic clusters, which shows a half-wave potential of up to 0.836V (vs. RHE) in the acid environment. More remarkably, it delivers a high power density of 870 and 478mW cm-2 at 1.0bar in H2 -O2 and H2 -Air fuel cells, respectively. According to theoretical calculation and in situ spectroscopy, the axial Fe4 C can provide strong electronic perturbation to Fe-N4 active centers, leading to the d-orbital electron delocalization of Fe and forming the Fe-pyrrolic N4 bond with high charge distribution, which stabilizes the Fe-pyrrolic N4 structure and optimizes the OH* adsorption during the catalytic process. This work proposes a new strategy to adjust the electronic structure of single-atom catalysts based on the strong interaction between single atoms and atomic clusters.

MXene-supported transition metal single-atom catalysts for nitrogen dissociation

Industrial ammonia production follows the Haber-Bosch process, whose rate-limiting step is the dissociation of the nitrogen molecule (N2), hence requiring suitable catalysts to break its triple bond. MXenes, a class of two-dimensional transition metal carbides and nitrides, have been proposed as very efficient catalysts for N2 dissociation. Here, by employing density functional theory-based calculations, we assess whether the deposition of one atom of a transition metal element (TM) on the Ti2C MXene surface further improves the catalytic potential of the MXene, serving as a single-atom catalyst. The results show that, for 21 of the 30 TMs considered, N2 can exothermically bind to the TM adatom, this bonding being favourable with respect to adsorption on the pristine Ti2C MXene surface for TMs of groups 3 to 6 of the Periodic Table. All the 21 TMs that successfully bind to N2 effectively reduce the N2 dissociation energy barrier when compared to the bulk Ti2C MXene by 18 to 84 %. Our results strongly indicate that doping the Ti2C MXene with atoms of transition metal elements significantly reduces the energy required to break the triple bond in N2, which may impact the nitrogen-to-ammonia process.