Search NASA⌕ Search

SEARCH · Search NASA

Results for “microkinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational Investigation of the Catalytic Hydrodeoxygenation of Propanoic Acid over a Cu(111) Surface

Cu-based alloy catalysts have recently been investigated experimentally for the hydrodeoxygenation (HDO) of biomass-derived organic acids. Here, the HDO of propanoic acid (PAc) has been studied over Cu(111) by mean-field microkinetic modeling based on parameters obtained from first-principles calculations. Models were developed for the gas- and liquid-phase HDO in condensed water and 1,4-dioxane. In agreement with experimental observations, the gas-phase PAc conversion rate is low at 573 K and increases in liquid water by 1 order of magnitude. In all reaction environments, the decarboxylation mechanism is dominant at low hydrogen partial pressures less than 0.1 bar, and the C–COO bond dissociation is the rate-controlling elementary step. This observation contrasts with the rate-controlling step identified over most group VIII metal surfaces, which is the C–OH bond dissociation in the decarbonylation mechanism. At high hydrogen (H 2 ) partial pressures greater than 10 bar, the HDO of PAc produces propionaldehyde that can readsorb and further react through decarbonylation to produce C 2 alkane products, which is conceptually different from the low H 2 partial pressure scenario. At high H 2 partial pressures, the initial hydrogenation at the carbonyl carbon of PAc becomes the rate-controlling elementary step.

Catalysts↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Insights into the Hydrogen Evolution Reaction on 2D Transition-Metal Dichalcogenides

Understanding hydrogen evolution reaction (HER) behaviors over two-dimensional transition-metal dichalcogenides (2D-TMDs) is critical for the development of nonprecious HER electrocatalysts with better activity. Here, in this work, by combining density functional theory calculations with microkinetic modeling, we thoroughly investigated the HER mechanism on 2D-TMDs. We find an important dependence of simulated cell size on the calculated hydrogen adsorption energy and the activation barrier for MoS 2 . Distinct from previous “H migration” mechanisms proposed for the Heyrovsky reaction, the rate-determining step for MoS 2 , we propose that the Mo site only serves as the stabilized transition state rather than H adsorption. In comparison to transition-metal electrocatalysts, we find that the activation barrier of the Heyrovsky reaction on 2D-TMDs scales with the hydrogen adsorption energy exactly as for transition metals except that all activation energies are displaced upward by ca. 0.4 eV. This higher Heyrovsky activation barrier is responsible for the substantially lower activity of 2D-TMDs. We further show that this higher activation barrier stems from the more positively charged adsorbed hydrogen on the chalcogenides interacting repulsively with the incoming proton. Based on these insights, we discuss potential strategies for the design of nonprecious HER catalysts with activity comparable to Pt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Point Defects Control Guest Molecule Diffusion in the 1D Pores of Zn(tbip)

Molecular diffusion plays a critical role in metal-organic frameworks (MOFs) within the application of kinetic chemical separations. We carefully study in this work the unexpected role of point defects for short-chain alkanes diffusing in Zn(tbip), an MOF with rigid one-dimensional (1D) channels. Inside a defect-free Zn(tbip), guest molecules are expected to follow single-file diffusion along 1D channels. It has been found previously by Heinke et al. that these parallel 1D channels are connected for molecular diffusion. Our density functional theory (DFT) calculations suggest that linker vacancy defects could arise under experimentally relevant conditions by removing a pair of linkers. Further climbing-image nudged elastic band (cNEB) DFT calculations indicate that hopping of short-chain alkanes between adjacent 1D channels over defect windows can occur at moderate temperatures. In addition to providing insights into connected adjacent 1D channels in Zn(tbip), Heinke et al. also inferred that most 1D pores are blocked from a microkinetic model. Additionally, we explored the influence of hydrolyzed linker created by the formation of linker removal inside 1D pores. Our DFT calculations show that the linkers can effectively block the pores and the linker diffusion in 1D channels is slow. Our results, for the first time, offer a mechanistic explanation of the unexpected molecular diffusion behavior in this MOF with 1D channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Unraveling 1-Hexene Hydrogenation over Dilute Pd-in-Au Alloys

Dilute Pd-in-Au alloys are valuable materials for selectively hydrogenating alkynes and isomerizing alkenes. By diluting Pd in a Au host, the selectivity toward semihydrogenated alkene isomers can be significantly enhanced and the unfavorable overhydrogenation to alkanes is suppressed. However, a detailed mechanistic study on the origin of the enhanced alkene selectivity over dilute alloy catalysts is still missing. Here, we combine experiment and theory to unravel the reaction mechanism, identifying rate-limiting and selectivity-controlling steps in 1-hexene hydrogenation over dilute Pd-in-Au catalysts. Using isotope-exchange hydrogenation experiments, we show that 1-hexene and hydrogen over a bimetallic Pd4Au96 in silica catalyst preferentially form 1-hexene isomers, (trans and cis) 2- and 3-hexene and only a small amounts of hexane. The reaction is consistent with a Horiuti–Polanyi mechanism, similar to a monometallic Pd nanoparticle catalyst. Computation of the free-energy profiles for 1-hexene hydrogenation and isomerization over a single Pd atom in a Au surface using first principles calculations indicated that the isomerization of 1-hexene to 2-hexene is energetically favorable due to the relatively large barrier for H2 dissociation preventing hydrogenation to n-hexane. Microkinetic modeling established that H2 dissociation on the single-atom Pd sites and H spillover from these sites onto the Au host are rate-limiting and key in steering the selectivity of dilute Pd-in-Au alloys toward the hexene isomers. Furthermore, the mechanistic insights from this study contribute to the rational design of optimized dilute alloy catalysts for selective alkene isomerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modifying Interface Solvation and Oxygen Reduction Electrocatalysis with Hydrophobic Species

Ionic liquid interlayers improve the oxygen reduction reaction (ORR) kinetics on bulk and nanostructured catalysts for both Pt and alloyed-Pt materials. Despite the demonstrated performance enhancement at the half-cell and membrane electrode assembly level, the mechanism of the improvement is not fully understood. In this work, we combine single-crystal experiments with microkinetic modeling to uncover the origin of the ORR kinetic improvement on Pt(111) in the presence of ionic liquids and hydrophobic cations. With the incorporation of a modified Frumkin isotherm, our model accurately simulates the disorder–order transition observed in hydroxyl and bisulfate adsorption on Pt(111) under acidic conditions. Voltametric analysis shows that ionic liquids impact solvation to break so-called scaling relations between the adsorption strength of OH ad and O ad , but these effects have little impact on ORR activity. Instead, destabilized OH ad reduces the overall hydroxyl (spectator) coverage, resulting in higher availability of active sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Understanding the Unique Reactivity of Cu for Electrochemical CO 2 Reduction with a 3-Site Model

Cu-based catalysts for the electrochemical reduction of CO 2 and CO exhibit a perplexingly unique reactivity toward multicarbon based products compared to other studied electrocatalysts. Here, in this study, we use insights gained from a recent phenomenological 3-site microkinetic model and grand-canonical density functional theory calculations to clarify the importance of an underemphasized aspect critical to Cu’s unique reactivity: a population of so-called “reservoir” sites. Using model Cu surface motifs, we discuss how these types can be represented by undercoordinated structural defects like step edges and grain boundaries which form a network of highly anisotropic migration channels. These pathways are found to be amenable for feeding *CO over time to reactive sites like Cu adatoms more active toward C–C coupling. These results highlight an often overlooked aspect of catalyst optimization: reservoir site engineering, which exploits surface mobility and presents an equally important avenue for electrocatalyst engineering oversimply maximizing active site densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2

Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.

08 HYDROGEN↗
DOE OSTI2020

Deconvolution of Water-Splitting on the Triple-Conducting Ruddlesden–Popper-Phase Anode for Protonic Ceramic Electrolysis Cells

Triple-conducting materials have been proved to improve the performance of popular protonic ceramic electrolysis cells. However, partially because of the complexity of the water splitting reaction involving three charge carriers, that is, oxygen (O 2– ), proton (H + ), and electron (e – ), the triple-conducting reaction mechanism was not clear, and the reaction conducting pathways have seldom been addressed. In this study, the triple conducting Ruddlesden–Popper phase Pr 1.75 Ba 0.25 NiO 4+δ as an anode on the BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3–δ electrolyte was fabricated and its electroresponses were characterized by electrochemical impedance spectroscopy with various atmospheres and temperatures. The impedance spectra are deconvoluted by means of the distribution of the relaxation time method. The surface exchange rate and chemical diffusivity of H + and O 2– are characterized by electrical conductivity relaxation. The physical locations of electrochemical processes are also identified by atomic layer deposition with a surface inhibitor. A microkinetics model is proposed toward conductivities, triple-conducting pathways, reactant dependency, surface exchange and bulk diffusion capabilities, and other relevant properties. Lastly, the rate-limiting steps and suggestions for further improvement of electrode performance are presented.

25 ENERGY STORAGE↗
DOE OSTI2020

Plasma-Catalytic Ammonia Synthesis beyond the Equilibrium Limit

In this work, we explore the consequences of nonthermal plasma-activation on product yields in catalytic ammonia synthesis, a reaction that is equilibrium-limited at elevated temperatures. We employ a minimal microkinetic model that incorporates the influence of plasma-activation on N 2 dissociation rates to predict NH 3 yields into and across the equilibrium-limited regime. NH 3 yields are predicted to exceed bulk thermodynamic equilibrium limits on materials that are thermal-rate-limited by N 2 dissociation. In all cases, yields revert to bulk equilibrium at temperatures at which thermal reaction rates exceed plasma-activated ones. Beyond-equilibrium NH 3 yields are observed in a packed bed dielectric barrier discharge reactor and exhibit sensitivity to catalytic material choice in a way consistent with model predictions. The approach and results highlight the opportunity to exploit synergies between nonthermal plasmas and catalysts to affect transformations at conditions inaccessible through thermal routes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Near-Surface Imaging of the Multicomponent Gas Phase above a Silver Catalyst during Partial Oxidation of Methanol

Fundamental chemistry in heterogeneous catalysis is increasingly explored using operando techniques in order to address the pressure gap between ultrahigh vacuum studies and practical operating pressures. Because most operando experiments focus on the surface and surface-bound species, there is a knowledge gap of the near-surface gas phase and the fundamental information the properties of this region convey about catalytic mechanisms. In this work, we demonstrate in situ visualization and measurement of gas-phase species and temperature distributions in operando catalysis experiments using complementary near-surface optical and mass spectrometry techniques. The partial oxidation of methanol over a silver catalyst demonstrates the value of these diagnostic techniques at 600 Torr (800 mbar) pressure and temperatures from 150 to 410 °C. Planar laser-induced fluorescence provides two-dimensional images of the formaldehyde product distribution that show the development of the boundary layer above the catalyst under different flow conditions. Raman scattering imaging provides measurements of a wide range of major species, such as methanol, oxygen, nitrogen, formaldehyde, and water vapor. Near-surface molecular beam mass spectrometry enables simultaneous detection of all species using a gas sampling probe. Detection of gas-phase free radicals, such as CH 3 and CH 3 O, and of minor products, such as acetaldehyde, dimethyl ether, and methyl formate, provides insights into catalytic mechanisms of the partial oxidation of methanol. The combination of these techniques provides a detailed picture of the coupling between the gas phase and surface in heterogeneous catalysis and enables parametric studies under different operating conditions, which will enhance our ability to constrain microkinetic models of heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Direct Comparison of the Activity and Selectivity of Rh 1 Cu and Ni 1 Cu Single-Atom Alloy Sites for Ethanol Decomposition

Ethanol is an important source of clean hydrogen, acetaldehyde, acetic acid, acetate esters, and light hydrocarbons. Controlling the divergent reaction pathways to these products requires understanding how different active sites influence the elementary steps involved. Herein, we present a combined surface science, theory, and nanoparticle catalysis study demonstrating how two single-atom dopants (Rh and Ni) in a Cu host can distinctively alter the selectivity of alcohol conversion. Specifically, our model studies reveal that ethanol reacts on Ni 1 Cu single-atom alloys to selectively produce acetaldehyde, whereas methane and CO are also formed on Rh 1 Cu single-atom alloys. Interestingly, these different reactivities are in contrast to the behavior of the pure metals as Ni(111) and Rh(111) surfaces favor methane/CO and surface carbon/CO, respectively. DFT calculations of reaction pathways and simulated product desorption based on microkinetic analyses explain these reactivity differences, demonstrating that C–C cleavage leading to methane formation has a lower barrier on Rh single-atom sites. To test the catalytic relevance of these fundamental results we synthesized and characterized supported Ni 1 Cu and Rh 1 Cu single-atom alloy nanoparticles with dopant:Cu ratios of 1:200. Flow reactor results revealed that both Ni and Rh increased ethanol conversion over Cu and that Ni 1 Cu catalysts were >99.9% selective to acetaldehyde, while Rh 1 Cu also produced 0.6%–2.6% of equimolar methane and CO between 433 and 493 K, demonstrating that C–C bond cleavage is enabled by isolated Rh sites. Furthermore, these catalytic results bridge the pressure and materials gaps, and together, this study provides insights into how different isolated dopant sites promote different catalytic pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Active Sites in the Dealuminated Beta Zeolite-Supported Cobalt Catalyst for Non-Oxidative Ethane Dehydrogenation

Dispersed metal species in siliceous zeolites have been actively studied for non-oxidative dehydrogenation of ethane (NDE). Fundamental insights into the dynamics of metal species in zeolites under reaction conditions have rarely been explored. Herein, we report an atomic level understanding of the dynamics and activity of cobalt (Co) sites in dealuminated Beta zeolite (DeAl-BEA) for NDE during induction and reaction conditions with extensive characterization techniques such as diffuse reflectance UV–vis, solid state nuclear magnetic resonance and X-ray photoelectron, X-ray diffraction along with in situ Fourier transform infrared and X-ray absorption spectroscopy. For a catalyst with 0.5 mass % Co loading, tetrahedral Co 2+ mononuclear sites, di-coordinated to the zeolite framework and with two silanol groups in vicinity (i.e., (≡SiO) 2 Co(HO–Si≡) 2 ), form upon exposure to hydrogen during induction and persist through the NDE reaction. Increasing the Co loading to 3.0 mass % yielded Co sites with similar electronic and coordination structures but slightly elongated Co–O bonds. Upon cooling to room temperature, the Co sites persisted in the same coordination environment, though the disappearance of a feature in the Co K-edge near-edge region revealed changes in the active site’s electronic structure coinciding with modest shifts in bond lengths. The electronic structure and activity of (≡SiO) 2 Co(HO–Si≡) 2 sites were studied comparatively to a few other hypothetical Co 2+ coordination structures, using electronic structure calculations and microkinetic simulations. The simulations showed that NDE is controlled by β-hydride elimination following C–H bond activation and that Co-sites possessing flexibility because of neighboring silanol defects are more active. Interestingly, dinuclear Co–O–Co sites (i.e., (≡SiO)Co(HO–Si≡) 2 –O–(HO–Si≡) 2 Co(≡SiO)) were more active than the mononuclear (≡SiO) 2 Co(HO–Si≡) 2 sites because of favorable hydrogen bonding with the vicinal silanol groups. In conclusion, the present study bridges the gap between the knowledge acquired by ex-situ characterizations and the active sites under the reaction conditions in alkane dehydrogenation chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling the Structural Sensitivity of Metal Catalysts in Ethylene Hydroformylation: Insights from Theory and Experiments

Here, in this study, we combined experimental and theoretical methods to investigate the structural sensitivity of metal catalysts in the ethylene hydroformylation reaction. Among Rh, Pt, Ir, Ni, Au, Ag, Pd, and Cu catalysts studied using experimental and theoretical methods, Rh showed the highest selectivity toward the C-C coupling product from CO and C 2 H 5 (i.e., C 2 H 5 CHO). The results from DFT and microkinetic simulations revealed that the activation energy barrier for C-C coupling is lowest on the Rh nanocluster, which explains the experimentally observed highest C 2 H 5 CHO selectivity on the Rh catalyst. Furthermore, DFT results demonstrated that the sites located on the flat surfaces of nanoparticles primarily promote the hydrogenation reaction, leading to the formation of undesired C 2 H 6 . In contrast, undercoordinated edge and corner sites of the nanocluster promote the C-C coupling reaction. Thus, our results illustrate that the selectivity toward C 3 oxygenates in ethylene hydroformylation reaction can be steered by tuning the size of Rh nanoparticles (the best-performing catalyst) to optimize the active (edge and corner) sites that preferentially promote the C-C coupling reaction.

58 GEOSCIENCES↗
DOE OSTI2025

Effect of Electrolyte Ions on Iridium Oxide-Based Water Oxidation Catalysis

The oxygen evolution reaction often limits the efficiency of renewable fuel syntheses due to its sluggish reaction kinetics. Of the factors that have been studied to improve this important reaction is how the choice of the electrolyte may alter the reaction kinetics. Despite its importance, sys-tematic studies of this effect have been relatively rare. Herein, we report an effort toward correcting this deficiency by investigating the effect of nitrate on water oxidation catalyzed by IrO x . The results show that nitrate can suppress the reaction, resulting in a decrease of the rate and an increase in the Tafel slope. The effect was found to be consistent with a microkinetic model incorporating competitive adsorption between reac-tion intermediates and nitrate, suggesting that the reaction mechanism was unaffected by anion identity. Moreover, this blocking effect exhibited dependence on the cations, following a trend Li + ≈ Na + ≈ K + > Cs + > TEA + . The results are expected to find broad applications in electrocatalysis.

Tafel slope analysis↗
DOE OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
DOE OSTI2025

Enhancing the Optically Detected Magnetic Resonance Signal of Organic Molecular Qubits

In quantum information science and sensing, electron spins are often purified into a specific polarization through an optical-spin interface, a process known as optically detected magnetic resonance (ODMR). Diamond-NV centers and transition metals are both excellent platforms for these so-called color centers, while metal-free molecular analogues are also gaining popularity for their extended polarization lifetimes, milder environmental impacts, and reduced costs. In our earlier attempt at designing such organic high-spin π-diradicals, we proposed to spin-polarize by shelving triplet M S = ±1 populations as singlets. This was recently verified by experiments albeit with low ODMR contrasts of <1% at temperatures above 5 K. In this work, we propose to improve the ODMR signal by moving singlet populations back into the triplet M S = 0 sublevel, designing a true carbon-based molecular analogue to the NV center. Our proposal is based upon transition-orbital and group-theoretical analyses of beyond-nearest-neighbor spin–orbit couplings, which are further confirmed by ab initio calculations of a realistic trityl-based radical dimer. Microkinetic analyses point toward high ODMR contrasts of around 30% under experimentally feasible conditions, a stark improvement from previous works. Finally, in our quest toward ground-state optically addressable molecular spin qubits, we exemplify how our symmetry-based design avoids Zeeman-induced singlet–triplet mixings, setting the scene for realizing electron spin qubit gates.

Group theory↗