Efficient Transition State Optimization of Periodic Structures through Automated Relaxed Potential Energy Surface Scans

Ligand-modified nanoparticle surfaces influence CO electroreduction selectivity

Improving the kinetics and selectivity of CO2/CO electroreduction to valuable multi-carbon products is a challenge for science and is a requirement for practical relevance. Here we develop a thiol-modified surface ligand strategy that promotes electrochemical CO-to-acetate. We explore a picture wherein nucleophilic interaction between the lone pairs of sulfur and the empty orbitals of reaction intermediates contributes to making the acetate pathway more energetically accessible. Density functional theory calculations and Raman spectroscopy suggest a mechanism where the nucleophilic interaction increases the sp2 hybridization of CO(ad), facilitating the rate-determining step, CO* to (CHO)*. We find that the ligands stabilize the (HOOC-CH2)* intermediate, a key intermediate in the acetate pathway. In-situ Raman spectroscopy shows shifts in C-O, Cu-C, and C-S vibrational frequencies that agree with a picture of surface ligand-intermediate interactions. A Faradaic efficiency of 70% is obtained on optimized thiol-capped Cu catalysts, with onset potentials 100 mV lower than in the case of reference Cu catalysts.

Cascade electrocatalysis via AgCu single-atom alloy and Ag nanoparticles in CO2 electroreduction toward multicarbon products

Electrocatalytic CO2 reduction into value-added multicarbon products offers a means to close the anthropogenic carbon cycle using renewable electricity. However, the unsatisfactory catalytic selectivity for multicarbon products severely hinders the practical application of this technology. In this paper, we report a cascade AgCu single-atom and nanoparticle electrocatalyst, in which Ag nanoparticles produce CO and AgCu single-atom alloys promote C-C coupling kinetics. As a result, a Faradaic efficiency (FE) of 94 ± 4% toward multicarbon products is achieved with the as-prepared AgCu single-atom and nanoparticle catalyst under ~720 mA cm-2 working current density at -0.65 V in a flow cell with alkaline electrolyte. Density functional theory calculations further demonstrate that the high multicarbon product selectivity results from cooperation between AgCu single-atom alloys and Ag nanoparticles, wherein the Ag single-atom doping of Cu nanoparticles increases the adsorption energy of *CO on Cu sites due to the asymmetric bonding of the Cu atom to the adjacent Ag atom with a compressive strain.

Structure-dependence and metal-dependence on atomically dispersed Ir catalysts for efficient n-butane dehydrogenation

Single-site pincer-ligated iridium complexes exhibit the ability for C-H activation in homogeneous catalysis. However, instability and difficulty in catalyst recycling are inherent disadvantages of the homogeneous catalyst, limiting its development. Here, we report an atomically dispersed Ir catalyst as the bridge between homogeneous and heterogeneous catalysis, which displays an outstanding catalytic performance for n-butane dehydrogenation, with a remarkable n-butane reaction rate (8.8 mol·g_Ir^-1·h^-1) and high butene selectivity (95.6%) at low temperature (450 °C). Significantly, we correlate the BDH activity with the Ir species from nanoscale to sub-nanoscale, to reveal the nature of structure-dependence of catalyst. Moreover, we compare Ir single atoms with Pt single atoms and Pd single atoms for in-depth understanding the nature of metal-dependence at the atomic level. From experimental and theoretical calculations results, the isolated Ir site is suitable for both reactant adsorption/activation and product desorption. Its remarkable dehydrogenation capacity and moderate adsorption behavior are the key to the outstanding catalytic activity and selectivity.

Anharmonic Correction to Free Energy Barriers from DFT-Based Molecular Dynamics Using Constrained Thermodynamic Integration

For the calculation of anharmonic contributions to free energy barriers, constrained thermodynamic λ-path integration (λ-TI) from a harmonic reference force field to density functional theory is presented as an alternative to the established Blue Moon ensemble method (ξ-TI), in which free energy gradients along the reaction coordinate ξ are integrated. With good agreement in all cases, the λ-TI method is benchmarked against the ξ-TI method for several reactions, including the internal CH₃ group rotation in ethane, a nucleophilic substitution of CH₃Cl, a retro-Diels-Alder reaction, and a proton transfer in zeolite H-SSZ-13. An advantage of λ-TI is that one can use virtually any reference state to compute anharmonic contributions to reaction free energies or free energy barriers. This is particularly relevant for catalysis, where it is now possible to compute anharmonic corrections to the free energy of a transition state relative to any reference, for example, the most stable state of the active site and the reactants in the gas phase. This is in contrast to ξ-TI, where free energy barriers can only be computed relative to an initial state with all reactants coadsorbed. Finally, the Bennett acceptance ratio method combined with λ-TI is demonstrated to reduce the number of required integration grid points with tolerable accuracy, favoring thus λ-TI over ξ-TI in terms of computational efficiency.

A computational investigation of the decomposition of acetic acid in H-SSZ-13 and its role in the initiation of the MTO process

The zeolite-catalyzed reaction of acetic acid is important in the direct utilization of biomass and also plays a role in the reactivity of oxygenates in the methanol-to-olefins (MTO) process.

Fully-exposed Pt-Fe cluster for efficient preferential oxidation of CO towards hydrogen purification

Hydrogen is increasingly being discussed as clean energy for the goal of net-zero carbon emissions, applied in the proton-exchange-membrane fuel cells (PEMFC). The preferential oxidation of CO (PROX) in hydrogen is a promising solution for hydrogen purification to avoid catalysts from being poisoned by the trace amount of CO in hydrogen-rich fuel gas. Here, we report the fabrication of a novel bimetallic Pt-Fe catalyst with ultralow metal loading, in which fully-exposed Pt clusters bonded with neighbor atomically dispersed Fe atoms on the defective graphene surface. The fully-exposed PtFe cluster catalyst could achieve complete elimination of CO through PROX reaction and almost 100% CO selectivity, while maintaining good stability for a long period. It has the mass-specific activity of 6.19 (molCO)*(gPt)−1*h−1 at room temperature, which surpasses those reported in literatures. The exhaustive experimental results and theoretical calculations reveal that the construction of fully-exposed bimetallic Pt-Fe cluster catalysts with maximized atomic efficiency and abundant interfacial sites could facilitate oxygen activation on unsaturated Fe species and CO adsorption on electron-rich Pt clusters to hence the probability of CO oxidation, leading to excellent reactivity in practical applications.

Theoretical investigation of the olefin cycle in H-SSZ-13 for the ethanol-to-olefins process using ab initio calculations and kinetic modeling

The formation of the hydrocarbon pool (HCP) in the ethanol-to-olefins (ETO) process catalyzed by H-SSZ-13 is studied in a kinetic model with ab initio computed reaction barriers.

A new mechanistic proposal for the aromatic cycle of the MTO process based on a computational investigation for H-SSZ-13

The paring mechanism of the aromatic cycle of the hydrocarbon pool is reinvestigated based on the heptamethylbenzenium cation adsorbed within H-SSZ-13 using quantum chemical calculations.

Against the NEER principle: the third type of photochromism for GFP chromophore derivatives

Photochromism is the heart of photochromic fluorescent proteins. Excited-state proton transfer (ESPT) is the major cause of photochromism for the green fluorescent protein (GFP) and Z-E photoisomerization through τ-torsion is the major cause of photochromism for Dronpa (a GFP mutant). In this article, s-E-1 opens a third type of photochromism for GFP chromophore derivatives, which involves light-driven φ-torsion with no τ-torsion, followed by excited-state intramolecular proton transfer (ESIPT), and is gated by environmental polarity. Since s-E-1 does not follow Z-E photoisomerization through τ-torsion but undergoes light-driven φ-torsion, which involves equilibration of the excited-state rotamers, it is clearly against the NEER (Non-Equilibration of Excited-state Rotamers) principle.

Tuning the selectivity of catalytic nitriles hydrogenation by structure regulation in atomically dispersed Pd catalysts

The product selectivity in catalytic hydrogenation of nitriles is strongly correlated with the structure of the catalyst. In this work, two types of atomically dispersed Pd species stabilized on the defect-rich nanodiamond-graphene (ND@G) hybrid support: single Pd atoms (Pd₁/ND@G) and fully exposed Pd clusters with average three Pd atoms (Pd_n/ND@G), were fabricated. The two catalysts show distinct difference in the catalytic transfer hydrogenation of nitriles. The Pd₁/ND@G catalyst preferentially generates secondary amines (Turnover frequency (TOF@333 K 709 h⁻¹, selectivity >98%), while the Pd_n/ND@G catalyst exhibits high selectivity towards primary amines (TOF@313 K 543 h⁻¹, selectivity >98%) under mild reaction conditions. Detailed characterizations and density functional theory (DFT) calculations show that the structure of atomically dispersed Pd catalysts governs the dissociative adsorption pattern of H₂ and also the hydrogenation pathway of the benzylideneimine (BI) intermediate, resulting in different product selectivity over Pd₁/ND@G and Pd_n/ND@G, respectively. The structure-performance relationship established over atomically dispersed Pd catalysts provides valuable insights for designing catalysts with tunable selectivity.

The selective hydrogenation of nitriles to prepare corresponding amines is highly desired in chemistry industry. Here, the authors selectively obtained secondary amines and primary amines over two types of atomically dispersed Pd catalysts supported on the nanodiamond-graphene hybrid support.

Electron penetration triggering interface activity of Pt-graphene for CO oxidation at room temperature

Achieving CO oxidation at room temperature is significant for gas purification but still challenging nowadays. Pt promoted by 3d transition metals (TMs) is a promising candidate for this reaction, but TMs are prone to be deeply oxidized in an oxygen-rich atmosphere, leading to low activity. Herein we report a unique structure design of graphene-isolated Pt from CoNi nanoparticles (PtǀCoNi) for efficiently catalytic CO oxidation in an oxygen-rich atmosphere. CoNi alloy is protected by ultrathin graphene shell from oxidation and therefore modulates the electronic property of Pt-graphene interface via electron penetration effect. This catalyst can achieve near 100% CO conversion at room temperature, while there are limited conversions over Pt/C and Pt/CoNiO_x catalysts. Experiments and theoretical calculations indicate that CO will saturate Pt sites, but O₂ can adsorb at the Pt-graphene interface without competing with CO, which facilitate the O₂ activation and the subsequent surface reaction. This graphene-isolated system is distinct from the classical metal-metal oxide interface for catalysis, and it provides a new thought for the design of heterogeneous catalysts.

Achieving CO oxidation at room temperature is significant for gas purification but remains challenging to perform. Here, the authors report design of graphene-isolated Pt from cobalt-nickel nanoparticles for efficiently catalytic CO oxidation in an oxygen-rich atmosphere.

Regulating coordination number in atomically dispersed Pt species on defect-rich graphene for n-butane dehydrogenation reaction

Metal nanoparticle (NP), cluster and isolated metal atom (or single atom, SA) exhibit different catalytic performance in heterogeneous catalysis originating from their distinct nanostructures. To maximize atom efficiency and boost activity for catalysis, the construction of structure-performance relationship provides an effective way at the atomic level. Here, we successfully fabricate fully exposed Pt3 clusters on the defective nanodiamond@graphene (ND@G) by the assistance of atomically dispersed Sn promoters, and correlated the n-butane direct dehydrogenation (DDH) activity with the average coordination number (CN) of Pt-Pt bond in Pt NP, Pt3 cluster and Pt SA for fundamentally understanding structure (especially the sub-nano structure) effects on n-butane DDH reaction at the atomic level. The as-prepared fully exposed Pt3 cluster catalyst shows higher conversion (35.4%) and remarkable alkene selectivity (99.0%) for n-butane direct DDH reaction at 450 °C, compared to typical Pt NP and Pt SA catalysts supported on ND@G. Density functional theory calculation (DFT) reveal that the fully exposed Pt3 clusters possess favorable dehydrogenation activation barrier of n-butane and reasonable desorption barrier of butene in the DDH reaction.

Effect of Impurities on the Initiation of the Methanol-to-Olefins Process: Kinetic Modeling Based on Ab Initio Rate Constants

Abstract

The relevance of a selection of organic impurities for the initiation of the MTO process was quantified in a kinetic model comprising 107 elementary steps with ab initio computed reaction barriers (MP2:DFT). This model includes a representative part of the autocatalytic olefin cycle as well as a direct initiation mechanism starting from methanol through CO-mediated direct C–C bond formation. We find that the effect of different impurities on the olefin evolution varies with the type of impurity and their partial pressures. The reactivity of the considered impurities for initiating the olefin cycle increases in the order formaldehyde < di-methoxy methane < CO < methyl acetate < ethanol < ethene < propene. In our kinetic model, already extremely low quantities of impurities such as ethanol lead to faster initiation than through direct C–C bond formation which only matters in complete absence of impurities.

Graphic Abstract

Theoretical investigation of the side-chain mechanism of the MTO process over H-SSZ-13 using DFT and ab initio calculations

The side-chain mechanism of the methanol-to-olefins process over the H-SSZ-13 acidic zeolite was investigated using periodic density functional theory with corrections from highly accurate ab intio calculations on large cluster models.

The transformation of cuboctahedral to icosahedral nanoparticles: atomic structure and dynamics

The rearrangement of transition metal nanoparticles from cuboctahedral to icosahedral structures is studied for up to 923 atoms. The atomic structure and temperature dependence of the transition are investigated with a well-defined collective variable. This collective variable describes the folding of the square fcc(100) facets into two triangular facets through a linear combination of the diagonals of all fcc(100) facets of all shells of the particle. Activation barriers are determined through harmonic transition state theory and constrained molecular dynamics simulations based on force field potentials. These calculations predict an activation entropy larger than 1 meV K^-1, leading to strongly temperature dependent activation barriers. Density functional theory calculations were additionally performed both as single point calculations and as full optimizations. Cu, Ag, Au and Ni clusters show low barriers for concerted, symmetric transition up to the 309-atomic clusters. In contrast, for Pd, Pt, Rh and Ir higher barriers are required, already for the 147-atomic clusters. With increasing barriers, an asymmetric but still concerted rearrangement becomes energetically more favorable than the fully symmetric transformation. The material-dependence of the transition can be correlated with the melting point of the bulk metals.

How Accurately Do Approximate Density Functionals Predict Trends in Acidic Zeolite Catalysis?

Density functional theory (DFT) is increasingly used for computational screening procedures with the aim of finding new catalysts. To achieve this, it is critical that relative differences between materials are predicted with high accuracy. How DFT at the generalized gradient approximation (GGA) level performs in this respect is investigated here for catalytic reactions employing acidic zeotypes using highly accurate DLPNO-CCSD(T) calculations as the reference. This is studied for 65 reaction energies and 130 reaction barriers related to zeolite catalysis. Our results obtained for the PBE-D3 and BEEF-vdW functionals show that while these functionals are prone to large errors, they predict trends occurring from one catalyst to another with an accuracy of about 5 kJ/mol, strongly supporting the widespread use of DFT calculations for the computational screening and design of new catalytic materials.

Describing Polytopal Rearrangements of Fluxional Molecules with Curvilinear Coordinates Derived from Normal Vibrational Modes: A Conceptual Extension of Cremer-Pople Puckering Coordinates

In this work a new curvilinear coordinate system is presented for the comprehensive description of polytopal rearrangements of N-coordinate compounds (N = 4-7) and systems containing an N-coordinate subunit. It is based on normal vibrational modes and a natural extension of the Cremer-Pople puckering coordinates ( J. Am. Chem. Soc. 1975, 97, 1354) together with the Zou-Izotov-Cremer deformation coordinates ( J. Phys. Chem. A 2011, 115, 8731) for ring structures to N-coordinate systems. We demonstrate that the new curvilinear coordinates are ideal reaction coordinates describing fluxional rearrangement pathways by revisiting the Berry pseudorotation and the lever mechanism in sulfur tetrafluoride, the Berry pseudorotation and two Muetterties' mechanisms in pentavalent compounds, the chimeric pseudorotation in iodine pentafluoride, Bailar and Ray-Dutt twists in hexacoordinate tris-chelates as well as the Bartell mechanism in iodine heptafluoride. The results of our study reveal that this dedicated curvilinear coordinate system can be applied to most coordination compounds opening new ways for the systematic modeling of fluxional processes.

S1/S0 Potential Energy Surfaces Experience Different Types of Restricted Rotation: Restricted Z/ E Photoisomerization and E/ Z Thermoisomerization by an Out-of-Plane Benzyl Group or In-Plane m-Pyridinium Group?

Any method that can enhance the fluorescence of fluorophores is highly desirable. Fluorescence enhancement accomplished by restricted Z/ E photoisomerization through intramolecular steric hindrance or relatively high bond order of a C═C double bond in a S₁ excited state has rarely been studied. In this article, we used green fluorescent protein (GFP) chromophore analogues as a model to get new physical insights into the restricted Z/ E photoisomerization and E/ Z thermoisomerization phenomena. We found that the S₁ and S₀ potential energy surfaces (PESs) of the GFP chromophore analogues experience two dramatically different types of restricted rotation, and 2b can be a representative example. In its S₁ PES, it is not the intramolecular steric hindrance between the out-of-plane benzyl group and the in-plane m-pyridinium group but the relatively high bond order of the I-bond in the S₁ excited state of 2b that makes it have a higher barrier for the Z/ E photoisomerization, a smaller Z/ E photoisomerization quantum yield, and a higher fluorescence quantum yield. In its S₀ PES, it is not the reduced bond order of the I-bond in the S₀ ground state of 2b but the intramolecular steric hindrance between the out-of-plane benzyl group and the in-plane m-pyridinium group that makes it have an extra higher barrier for E/ Z thermoisomerization and a much smaller E/ Z thermoisomerization rate constant.

Olefin methylation and cracking reactions in H-SSZ-13 investigated with ab initio and DFT calculations

The olefin cycle of the methanol-to-olefins process is investigated for the zeolite H-SSZ-13 using periodic, van-der-Waals corrected DFT calculations, together with MP2 corrections derived from cluster models, which are essential for accurate barriers.