Search

We report a bond-valence method (BVM) parameterization framework that captures density functional theory (DFT)-computed relative stabilities using the BVM global instability index (GII). We benchmarked our framework against a dataset of 188 experimentally observed ABO

The catalytic partial oxidation of methane is achieved at low temperatures (90% at ∼35% methane conversion at approximately 6 turnovers. Under our catalytic methane oxidation reaction conditions, MeTFA is stable against overoxidation, which explains the likely high selectivity at conversions >15%. Using combined experimental studies and DFT calculations, a mechanism involving soluble and molecular Mn species in the catalytic cycle is proposed. The proposed reaction pathway involves initial activation of Mnᴵᴵ by dioxygen, cleavage of a methane C–H bond by a Mnᴵᴵᴵ hydroxo intermediate, rebound of the methyl radical to generate MeTFA, and finally regeneration of the starting MnII complex. Also, this process is shown to be applicable to the oxidation of ethane, favoring the mono-oxidized product ethyl trifluoroacetate (EtTFA) and reaching ∼46% conversion.

We use quinone reduction potentials to showcase molecularly tunable CO2 binding strengths for a family of electrochemically generated quinoid molecular sorbents. In our study, we find that CO2 binding stoichiometry varies with CO2 concentration.

To understand the reaction mechanism and kinetics of N2 reduction reaction and competing hydrogen evolution reaction on the dual-atom catalyst as a function of applied potential by applying the Grand canonical potential kinetics.

A highly active heterogenized molecular CO2 reduction catalyst on a conductive carbon support is investigated to identify if its improved catalytic activity can be attributed to strong electronic interactions between catalyst and support. The molecular structure and electronic character of a [Re+1(tBu-bpy)(CO)3Cl] (tBu-bpy = 4,4′-tert-butyl-2,2′-bipyridine) catalyst deposited on multiwalled carbon nanotubes are characterized using Re L3-edge x-ray absorption spectroscopy under electrochemical conditions and compared to the homogeneous catalyst. The Re oxidation state is characterized from the near-edge absorption region, while structural changes of the catalyst are assessed from the extended x-ray absorption fine structure under reducing conditions. Chloride ligand dissociation and a Re-centered reduction are both observed under applied reducing potential. The results confirm weak coupling of [Re(tBu-bpy)(CO)3Cl] with the support, since the supported catalyst exhibits the same oxidation changes as the homogeneous case. However, these results do not preclude strong interactions between a reduced catalyst intermediate and the support, preliminarily investigated here using quantum mechanical calculations. Thus, our results suggest that complicated linkage schemes and strong electronic interactions with the initial catalyst species are not required to improve the activity of heterogenized molecular catalysts.

We report the synthesis, characterization, and electrocatalytic water oxidation activity of two cobalt complexes, (6-FP)Co(NO₃)₂ (1) (6-FP = 8,8′-(1,2-phenylene)diquinoline) and (5-FP)Co(NO₃)₂ (2) (5-FP = 1,2-bis(N-7-azaindolyl)benzene), containing “capping arene” bidentate ligands with nitrogen atom donors. The cobalt complexes 1 and 2 were supported on ordered mesoporous carbon (OMC) by π–π stacking, resulting in heterogenized cobalt materials 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, respectively, and studied for electrocatalytic water oxidation. We find that 6-FP-Co-OMC-1 exhibits an overpotential of 355 mV for a current density of 10 mA cm⁻² and a turnover frequency (TOF) of ∼0.53 s–1 at an overpotential of 400 mV at pH 14. 6-FP-Co-OMC-1 exhibits activity that is ∼1.6 times that of 5-FP-Co-OMC-2, which gives a TOF of 0.32 s⁻¹ at 400 mV overpotential. The structural stability of the single-atom Co site was demonstrated for 6-FP-Co-OMC-1 using X-ray absorption spectroscopy for the molecular complex supported on OMC, but slow degradation in catalyst activity can be attributed to eventual formation of Co oxide clusters. DFT computations of electrocatalytic water oxidation using the molecular complexes as models provide a description of the catalytic mechanism. These studies reveal that the mechanism for O–O bond formation involves an intermediate Coᴵⱽ oxo complex that undergoes an intramolecular reductive O–O coupling to form a Coᴵᴵ–OOH species. Further, the calculations predict that the molecular 6-FP-Co structure is more active for electrocatalytic water oxidation than 5-FP-Co, which is consistent with experimental studies of 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, highlighting the possibility that the ligand structure influences the catalytic activity of the supported molecular catalysts.

Support-induced strain engineering is a powerful strategy to modulate the electronic structure of two-dimensional materials. However, controlling strain of planar molecules such as metallophthalocyanines and metalloporphyrins is technically challenging due to their sub–2 nm lateral size. In addition, the effect of strain on molecular properties remains poorly understood. Starting with cobalt phthalocyanine (CoPc), a model CO2 reduction reaction (CO2RR) catalyst, we show that carbon nanotubes (CNTs) are ideal substrates for inducing optimum properties through molecular curvature. Using a tandem-flow electrolyzer with monodispersed CoPc on single-walled CNTs (CoPc/SWCNT) as the catalyst, we achieve a methanol partial current density of >90 mA cm-2 with a selectivity of >60%. CoPc on wide multi-walled CNTs (MWCNTs) leads to only 16.6% selectivity. We report X-ray spectroscopic characterizations to unravel the distinct local coordinations and electronic structures induced by the strong molecule-support interactions. These results agree with our Grand Canonical Density Functional Theory that calculates the energetics as a function of applied potential. We find that SWCNTs induce curvature in CoPc, which improves *CO binding to enable subsequent formation of methanol, while wide MWCNTs favor CO desorption. Thus, we demonstrate that the SWCNT-induced molecular strain increases methanol formation. We also show that induced strain can accelerate the oxygen reduction reaction and CO2RR for other catalysts. Our results show the important role of SWCNTs beyond catalyst dispersion and electron conduction.

The electronic structure and local coordination of binary (Mo6T8, T = S, Se, Te) and ternary Chevrel Phases (MxMo6T8, M = transition metal) are investigated for a range of metal intercalant and chalcogen compositions through XANES and EXAFS analyses.

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO

We show how to improve the yield of MeX from CH4 activation catalysts from 12% to 90% through the use of “capping arene” ligands. Four (FP)RhIII(Me)(TFA)2 {FP = “capping arene” ligands, including 8...

We report a combined experimental and computational study focused on the mechanism of oxidative conversion of benzene and ethylene to styrene using [(η2-C2H4)2Rh(μ-OAc)]2 as the catalyst precursor ...

Selective reduction of CO2 to formate represents an ongoing challenge in photoelectrocatalysis. To provide mechanistic insights, we investigate the kinetics of hydride transfer (HT) from a series o...

Electrochemical routes provide an attractive alternative to the Haber-Bosch process for cheaper and more efficient ammonia (NH3) synthesis from N2 while avoiding the onerous environmental impact of the Haber-Bosch process. We prototype a strategy based on a eutectic mixture of phosphate molten salt. Using quantum-mechanics (QM)-based reactive molecular dynamics, we demonstrate that lithium nitride (Li3N) produced from the reduction of nitrogen gas (N2) by a lithium electrode can react with the phosphate molten salt to form ammonia. We extract reaction kinetics of the various steps from QM to identify conditions with favorable reaction rates for N2 reduction by a porous lithium electrode to form Li3N followed by protonation from phosphate molten salt (Li2HPO4-LiH2PO4 mixture) to selectively form NH3.

The design of multinary solid-state material systems that undergo reversible phase changes via changes in temperature and pressure provides a potential means of safely storing hydrogen. However, fully mapping the stabilities of known or newly targeted compounds relative to competing phases at reaction conditions has previously required many stringent experiments or computationally demanding calculations of each compound's change in Gibbs energy with respect to temperature, G(T). In this work, we have extended the approach of constructing chemical potential phase diagrams based on ΔGf(T) to enable the analysis of phase stability at non-zero temperatures. We first performed density functional theory calculations to compute the formation enthalpies of binary, ternary, and quaternary compounds within several compositional spaces of current interest for solid-state hydrogen storage. Temperature effects on solid compound stability were then accounted for using our recently introduced machine learned descriptor for the temperature-dependent contribution Gδ(T) to the Gibbs energy G(T). From these Gibbs energies, we evaluated each compound's stability relative to competing compounds over a wide range of conditions and show using chemical potential and composition phase diagrams that the predicted stable phases and H2 release reactions are consistent with experimental observations. This demonstrates that our approach rapidly computes the thermochemistry of hydrogen release reactions for compounds at sufficiently high accuracy relative to experiment to provide a powerful framework for analyzing hydrogen storage materials. This framework based on G(T) enables the accelerated discovery of active materials for a variety of technologies that rely on solid-state reactions involving these materials.

Herein, we detail an approach to accelerate the computational screening of materials for properties dictated by the kinetics of solid-state diffusion through reliably and rapidly identifying upper and lower bounds to the transition state (TS) energy through our proposed modified single iteration synchronous-transit (MSIST) approach. While this sacrifices providing detailed information of the explicit TS structure, it requires only 30% of the force evaluations of a full nudged elastic band (NEB) TS search and reduces the computational demand to compute estimated diffusion barriers by ∼70% on average. In all 53 cases in which we explicitly compared our results to those of an NEB calculation, the upper and lower bounds identified using this approach bracketed the TS energy calculated with explicit NEB calculations. We use the applications of diffusion of Na

State-of-the-art high temperature oxide melt solution calorimetry and density functional theory were employed to produce the first systematic study of thermodynamic stability in a series of binary ...

In this work, we implement a facile microwave-assisted synthesis method to yield three binary Chevrel-Phase chalcogenides (Mo6X8; X = S, Se, Te) and investigate the effect of increasing chalcogen e...

Radical photopolymerization (RPP) has grown into a multibillion-dollar technology for reduced energy consumption and waste with increased productivity. However, radical-mediated polymerization ceas...

The iron aluminate spinel hercynite (FeAl2O4) is a promising redox material for solar thermochemical hydrogen production (STCH). Although it has a high H2 production capacity, the kinetics of its o...

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs2BB'Cl6) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (τ) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs2BB'Cl6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy

Materials with crystal structures containing tetrahedral motifs are preferable for optoelectronic applications because they often have direct band gaps and low electron effective masses. However, crystal structures of manganese chalcogenides typically contain octahedral motifs, such as in rock salt (RS) MnS and MnSe materials. Here, we experimentally show that MnS1−xSex alloys with tetrahedrally bonded wurtzite (WZ) structure can form between MnSe and MnS parent compounds with octahedral RS structures, at S-rich compositions (x < 0.4) and low synthesis temperatures (∼300 °C). The calculated mixing enthalpies of MnS1−xSex alloys in RS and WZ structures cannot explain this experimental observation, so we hypothesize that WZ stabilization may be related to smaller structure density and lower surface energy compared with RS. The resulting WZ MnS1−xSex alloys have 3.0–3.2 eV optical absorption onset and lower electrical conductivity (

This work examines the relative reactivities of ReI and MnI tricarbonyl pyridine-2,6-bis-N-heterocyclic carbene pincers M(CO)3CNCBnX (M = Re, Mn and X = Cl and Br) towards catalysis for the electrochemical conversion of CO2 to CO. Unlike prior well-studied group VII catalysts, Mn(CO)3CNCBnX is extraordinarily active, while the new Re(CO)3CNCBnX complex surprisingly does not exhibit catalytic response. DFT calculations shed light on this puzzling behavior and show that the redox-active pyridine-2,6-bis-N-heterocyclic carbene ligand facilitates the reduction of the ground-state complexes; however, the extent of electronic delocalization in the reduced intermediates differs in the degree of metalloradical character. The highly-active Mn(CO)3CNCBnX complex proceeds through an intermediate with nucleophilic metalloradical character in which 66% of the unpaired electron spin resides on Mn. In contrast, Re(CO)3CNCBnX reduction proceeds through an intermediate with less metalloradical character in which only 38% of the unpaired spin is localized on Re with the remainder delocalized over the ligand. The energetic penalty of the electron delocalization of an electron on the ligand affects the M-CO bond strengths and related kinetic barriers. We discuss these observations in the context of turnover-enabling effects in CO2 reductions mediated by group VII NHC pincer molecular electrocatalysts.

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

Ammonia is an essential chemical owing to its importance in fertilizer production and other industrial applications. Electrocatalytic nitrate reduction to ammonia (NO₃RR) holds great promise for low-temperature ammonia production while simultaneously addressing nitrate-based environmental concerns. To provide the mechanistic understanding needed to design an effective electrocatalyst, we systematically investigated the catalytic performance of metal-based dual-atom catalysts (DACs) anchored on two-dimensional (2D) expanded phthalocyanine (Pc) for NO₃RR. We found that NO3RR can efficiently produce ammonia on Cr₂-Pc, V₂-Pc, Ti₂-Pc, and Mn₂-Pc surfaces with low limiting potentials of − 0.02, − 0.25, − 0.34, and − 0.41 VRHE, respectively. Moreover, using the free energy difference of *NO₃⁻ and *H as a descriptor, we found that the hydrogen evolution reaction is significantly suppressed on the DAC surface due to an ensemble effect in which the two metal atoms cooperate to selectively form ammonia. We performed high-throughput screening to develop an efficient metal-based DAC for NO₃⁻ reduction, followed by a mechanistic study to elucidate the NO₃RR pathway on the DAC. This work provides design information for advancing sustainable ammonia synthesis under ambient conditions.

This work maps the thermodynamics of electrochemically generated C-nucleophiles for reactive capture of CO

A combined experimental and computational study of the reactivities of seven commonly used Michael acceptors paired with two thiols within the framework of photobase-catalyzed thiol-Michael reactions is reported. The rate coefficients of the propagation (k(P)), reverse propagation (k(−P)), chain-transfer (k(CT)), and overall reaction (k(overall)) were experimentally determined and compared with the well-accepted electrophilicity parameters of Mayr and Parr, and DFT-calculated energetics. Both Mayr’s and Parr’s electrophilicity parameters predict the reactivities of these structurally varying vinyl functional groups well, covering a range of overall reaction rate coefficients from 0.5 to 6.2 s(−1). To gain insight into the individual steps, the relative energies have been calculated using DFT for each of the stationary points along this step-growth reaction between ethanethiol and the seven alkenes. The free energies of the individual steps reveal the underlying factors that control the reaction barriers for propagation and chain transfer. Both the propagation and chain transfer steps are under kinetic control. These results serve as a useful guide for Michael acceptor selection to design and predict thiol-Michael-based materials with appropriate kinetic and material properties

Conversion plateaus rapidly in radical photopolymerizations (RPPs) following discontinuation of irradiation due to rapid termination of reactive radicals, which restricts the wider use of RPPs in applications that involve nonuniform light access including those with attenuated light transmission or irregular surfaces. Based on our recent report of a radical dark-curing photoinitiator (DCPI) that continues polymerization beyond the cessation of irradiation by enabling latent redox initiation with photo-released amine in the presence of a suitable oxidant, we developed a new DCPI with an absorption spectrum that extends well into the visible range. Our design process involved a series of computational investigations of candidate molecules, including a systematic study of substituents and their position-dependent effects on absorption characteristics, electronic transitions, and the photochemical mechanism and its associated energetics. Our quantum chemical computations identified the target compound 5,7-dimethoxy-6-bromo-3-aroylcoumarin-DMPT/BPh(4) and predicted that it would facilitate the dark-curing mechanism by concurrent photo-radical generation and photo-induced release of an efficient redox reductant under visible irradiation. This reductant-tethered chromophore was then synthesized and optically characterized with UV–vis spectroscopy that revealed its strong visible-light absorption with a molar absorptivity of 5710 M(−1) cm(−1) at 405 nm and 50 M(−1) cm(−1) at 455 nm. We then demonstrated extensive dark-curing of >35% additional conversion over 25 min following brief activation of the shelf-stable one-part system by irradiation with a 455 nm LED that was ceased at 20% conversion. In contrast, shuttering irradiation of the control formulation at that same point resulted in immediate cessation of conversion, which plateaued at 20%. We determined a remarkable initiator efficiency of 2.82 that results from the additional redox-generated radicals with a 77% photo-reductant generation quantum yield. The combination of superior photo- and dark-curing efficiencies of this new visible DCPI is expected to open new application opportunities in RPP, especially those involving resins that are highly light attenuating, surfaces that possess irregular features that produce uneven irradiance, and production lines where continued dark-curing downstream of the light activation step enhances line efficiencies.