Comparative simulations of methanol steam reforming on PdZn alloy using kinetic Monte Carlo and mean-field microkinetic model.

Methanol steam reforming (MSR) is an attractive route for producing clean energy hydrogen. PdZn alloys are extensively studied as potential MSR catalysts for their stability and high CO2 selectivity. Here, we investigated the reaction mechanism using density functional calculations, mean-field microkinetic modeling (MF-MKM), and kinetic Monte Carlo (kMC) simulations. To overcome the over-underestimation of CO2 selectivity by log-kMC, an ads-kMC algorithm is proposed in which the adsorption/desorption rate constants were reduced under certain requirements and the diffusion process was treated by redistributing surface species each time an event occured. The simulations show that the dominant pathway to CO2 at low temperatures is CH3OH → CH3O → CH2O → H2COOH → H2COO → HCOO → CO2. The ads-kMC predicted OH coverage is 2-3 times that of MF-MKM, while they produce similar coverage for other species. Analyses indicate that surface OH promotes the dehydrogenation of CH3OH, CH3O, and H2COOH significantly and plays a key role in the MSR process. The dissociation of water/methanol is the most important rate-limiting/rate-inhibiting step. The CO2 selectivity obtained by the two methods is close to each other and consistent with the experimental trend with temperature. Generally, the ads-kMC results agree with the MF-MKM ones, supporting the previous finding that kMC and MF-MKM predict similar results if the diffusion is very fast and adsorbate interactions are neglected. The present study sheds light on the MSR process on PdZn alloys, and the proposed scheme to overcome the stiff problems in kMC simulations is worthy of being extended to other systems.

CO2 Hydrogenation to CH3OH on Metal-Doped TiO2(110): Mechanisms, Strain Effect and a New Thermodynamic-Kinetic Relation.

Surface strain and linear thermodynamic-kinetic relation are interesting topics in catalysis. Development of low temperature methanol catalysts of high activity and selectivity is of particularly importance for conversion of CO2 to methanol. In the present paper CO2 hydrogenation to methanol on Znx@TiO2(110) (x=0-2) was explored using density functional calculations and microkinetic simulations. The reaction mechanisms on the three model systems were determined and it is shown that Zn2@TiO2(110) is the most active. The most favorable pathway on Zn2@TiO2(110) is identified and CO2+H to HCOO is found to be the rate-controlling step. It is demonstrated that there is a linear relation (named AEB relation) between the adsorption energies of the initial states and the barriers for the controlling step on the 18 systems studied. Calculations on strained surfaces show that the AEB relation exists within ±1 % strain. Sr2@TiO2(110) and -1 % strained CaZn and ZnCu doped TiO2(110) are potential good low temperature catalysts and deserve experimental testing.

Mechanistic Insight into Acceptorless Dehydrogenation of Methanol to Syngas Catalyzed by MACHO-Type Ruthenium and Manganese Complexes: A DFT Study.

The acceptorless dehydrogenation of methanol to produce carbon monoxide (CO) and dihydrogen (H2) mediated by MACHO-type 1-Ru and 1-Mn complexes was theoretically investigated via density functional theory calculations. The 1-Ru-catalyzed process involves the formation of active species 4-Ru through a methanol-bridged H2 release pathway. Methanol dehydrogenation by 4-Ru yields formaldehyde and 1-Ru, followed by H2 release to regenerate 4-Ru (rate-determining step, ΔG‡ = 32.5 kcal/mol). Formaldehyde further reacts with methanol via nucleophilic attack of the MeO- ligand in the Ru complex (ΔG‡ = 9.6 kcal/mol), which is more favorable than the traditional methanol-to-formaldehyde nucleophilic attack (ΔG‡ = 33.8 kcal/mol) due to the higher nucleophilicity of MeO-. CO is ultimately produced through the methyl formate decarbonylation reaction. Accelerated H2 release in the early reaction stage compared to CO results from the initial methanol dehydrogenation and condensation of formaldehyde with methanol. In contrast, CO generation occurs later via methyl formate decarbonylation. The 1-Mn-catalyzed reaction has reduced efficiency compared to 1-Ru for the higher Gibbs energy barrier (ΔG‡ = 34.1 kcal/mol) of the rate-determining step. Excess NaOtBu promotes the reaction of CO and methanol, forming methyl formate, significantly reducing the CO/H2 ratio as the catalyst amount decreases. These findings deepen our understanding of the methanol-to-syngas transformation and can drive progress in this field.

Metal-free catalytic hydroboration of imine with pinacolborane using a pincer-type phosphorus compound: mechanistic insight and improvement of the reaction.

A mechanistic study of the catalytic hydroboration of imine using a pincer-type phosphorus compound 1NP was performed through the combination of DFT and DLPNO-CCSD(T) calculations. The reaction proceeds through a phosphorus-ligand cooperative catalytic cycle, where the phosphorus center and triamide ligand work in a synergistic manner. First, the pinB-H bond activation by 1NP occurs through the cooperative functions of the phosphorus center and the triamide ligand, leading to a phosphorus-hydride intermediate 2NP. This is the rate-determining step, with the Gibbs energy barrier and Gibbs reaction energy of 25.3 and -17.0 kcal mol-1, respectively. Subsequently, the hydroboration of phenylmethanimine takes place through a concerted transition state through the cooperative function of the phosphorus center and the triamide ligand. It leads to the final hydroborated product 4 with the regeneration of 1NP. Our computational results reveal that the experimentally isolated intermediate 3NP is a resting state of the reaction. It is formed through the B-N bond activation of 4 by 1NP, rather than via the insertion of the CN double bond of phenylmethanimine into the P-H bond of 2NP. However, this side reaction can be suppressed by utilizing a planar phosphorus compound AcrDipp-1NP as the catalyst, which features steric-demanding substituents on the chelated N atom of the ligand.

Amalgamation of DNAzymes and Nanozymes in a Coronazyme.

Artificial enzymes such as nanozymes and DNAzymes are economical and stable alternatives to natural enzymes. By coating Au nanoparticles (AuNPs) with a DNA corona (AuNP@DNA), we amalgamated nanozymes and DNAzymes into a new artificial enzyme with catalytic efficiency 5 times higher than AuNP nanozymes, 10 times higher than other nanozymes, and significantly greater than most of the DNAzymes on the same oxidation reaction. The AuNP@DNA demonstrates excellent specificity as its reactivity on a reduction reaction does not change with respect to pristine AuNP. Single-molecule fluorescence and force spectroscopies and density functional theory (DFT) simulations indicate a long-range oxidation reaction initiated by radical production on the AuNP surface, followed by radical transport to the DNA corona, where the binding and turnover of substrates take place. The AuNP@DNA is named coronazyme because of its natural enzyme mimicking capability through the well-orchestrated structures and synergetic functions. By incorporating different nanocores and corona materials beyond DNAs, we anticipate that the coronazymes represent generic enzyme mimics to carry out versatile reactions in harsh environments.

A new adsorption energy-barrier relation and its application to CO2 hydrogenation to methanol over In2O3-supported metal catalysts.

Herein, we report a new adsorption energy-barrier relation, the adsorbate-dependent barrier scaling (ADBS) relation, with which the catalytic activity of In2O3-supported metal catalysts for CO2 hydrogenation to methanol is predicted. It is shown that Cu, Ga, NiPt and NiPd alloys exhibit high catalytic activity for CO2 hydrogenation to methanol.

Strain effect on adsorption and reactions of AHx (A = C, N, O, x ≤ 3) on In2O3(110), TiO2(110), and ZrO2(101) surfaces.

More and more attention has been paid to strain-based regulation of catalytic activity. To guide regulation of catalytic performance via strain engineering, adsorption and reactions of AHx (A = C, N, O, x ≤ 3) were investigated on uniformly strained In2O3 (110), rutile TiO2 (110), and tetragonal ZrO2 (101) from -2% to 4%. The results show that adsorption energies vary linearly with strain; expansive strain enhances the adsorption of most adsorbates. Unlike the adsorbate scaling relations that are central atom dependent, the adsorbate scaling relations on strained surfaces are central atom independent. C-H/O-H bonds are elongated/shortened with expansive strain, and adsorption energies of CHx generally change more than those of OHx and NHx, which can be rationalized with effective medium theory and pertinent bond energies. Thermodynamically, In2O3(110)/ZrO2(101) is the most active/inactive. The estimated variation of rate constants at 300 K from 0% to 2% strain based on the Brønsted-Evans-Polanyi relationship demonstrates great strain regulation potential of catalytic performance on these oxide surfaces. Finally, it is demonstrated that strain tends to facilitate the reactions whose sum of the stoichiometric number is positive, which can be used as a rule to guide strain engineering for heterogeneous catalysis.

DFT and microkinetic study of acetylene transformation on Pd(111), M(111) and PdM(111) surfaces (M = Cu, Ag, Au).

Density functional calculations and microkinetic simulations were performed on the transformation network of acetylene on Pd(111), M(111) and PdM(111) (M = Cu, Ag, Au) surfaces. It is demonstrated that the adsorption energies on alloy surfaces linearly correlate with the values on the pure metal surfaces. A good linear relationship between the co-adsorption energies of initial states and transition states is revealed with which the barriers of most elementary steps in the reaction network were estimated. To shed light on the transformation of acetylene, microkinetic simulations were conducted on the network. The results show that CHCH and H are dominant species on the surfaces and CCH, CCH2 and CCH3 are the main intermediates. Analysis indicates that introduction of coinage metals into Pd reduces the activity, but promotes the selectivity by lowering the barrier of CHCH2 → CH2CH2. The present work provides a comprehensive overview of acetylene transformation on palladium, coinage metals and their alloy surfaces. The linear relationship of adsorption energies between the component metal and alloy surfaces and usage of the TSS relationship to evaluate barriers for microkinetic simulations are worthy of being further studied and extended to other systems.

Theoretical Studies on Bimetallic Salen Complexes Catalyzed Epoxide Hydration: Effects of Metal Centers, Substrates, and Ligands.

To provide guiding information for developing efficient and stable catalysts for epoxide hydration, we investigated the mechanism of propylene oxide (PO) to 1,2-propylene glycol (PG) using density functional theory (DFT) calculations. The mechanism was identified to follow the cooperative bimetallic mechanism in which a metal-salen complex activated H2O attacks the middle carbon atom of a metal-salen complex activated PO from the oxygen side of three-membered ring. Analyses reveal that the distortion energy correlates linearly with the barrier, and the hydrogen bonding between H2O and PO increases from reaction precursors to transition states. A nice linear relationship exists between the ratio of square root of ionic potential to the square of the distance from the metal ion spherical surface to the oxygen atom center of PO. It is demonstrated that the substrates with larger polarizability tend to have lower hydration barriers and the influence of ligands is less than that of metal centers and substrates. Modifying metal ions is the first choice for developing metal-salen catalysts, and metal ions with more formal charges and larger radius are expected to exhibit high activity. These findings shed lights on the mechanism and provide guiding information for developing efficient metal-salen catalysts for epoxide hydration.

Rationalization of Nonlinear Adsorption Energy-Strain Relations and Brønsted-Evans-Polanyi and Transition State Scaling Relationships under Strain.

Scaling relations play a vital role in high-throughput screening of catalytic materials, and more and more attention is being paid to strain-based regulation of catalytic performance. Here we investigated the variation of several energetics, including adsorption energies in the initial state, transition state, and final state, reaction energies, and energy barriers with strain, by studying CO, BH, NH, CH, and NO adsorption and dissociation on M(111) (M = Cu, Ag, Ni, Pd, or Pt) surfaces. We show that energy barriers, reaction energies, and adsorption energies can vary either linearly or nonlinearly (quadratically) with strain. Systems with stronger adsorbate-substrate interaction and weaker atom-atom interaction in substrates are more likely to exhibit nonlinear relations. The well-known Brønsted-Evans-Polanyi relationships and transition state scaling relationships under strain were also examined, and both of them can be nonlinear under strain, in principle. The observed nonlinear relationships were satisfactorily rationalized with the equations derived from Mechanics of Solids.

Exploring the mechanism and counterion activity regulation in the CoIII(salen)-catalyzed hydration of propylene oxide.

CoIII(salen)-X (X = Cl-, OAc-, and OTs-) mediated hydration of propylene oxide (PO) to propylene glycol has been investigated in detail using density functional theory (DFT) calculations. Two kinds of reaction mechanisms, the concerted and stepwise pathways, were scrutinized. For the eight concerted routes, the cooperative bimetallic route in which the middle carbon atom is attacked by the nucleophilic oxygen atom (route VI-m) was calculated to be the most favorable, and among the three catalysts examined H2O-CoIII-OTs was found to be the most active, due to the strong hydrogen bonding between the nucleophilic H2O and the ring oxygen atom in the epoxides as well as the extra π-π stacking interaction. For the stepwise mechanism which consists of the formation of H2O-CoIII-OH, the ring-opening of PO and propylene glycol formation, our studies reveal that different H2O-CoIII-Xs behave kinetically very similarly in the course of propylene glycol formation, but show a notable difference in the rate of H2O-CoIII-OH formation with Cl- > OAc- > OTs-. The rate ordering with which we rationalize the experimental phenomena well is disclosed to be consistent with the nucleophilicity of the counterions by molecular electrostatic potential, condensed Fukui function and condensed local softness. We show that the nucleophilicity of the counterion determines the favorable mechanism that PO hydration follows.

Rapid room temperature synthesis of red iridium(iii) complexes containing a four-membered Ir-S-C-S chelating ring for highly efficient OLEDs with EQE over 30.

Three red cyclometalated iridium(iii) complexes (4tfmpq)2Ir(dipdtc), (4tfmpq)2Ir(dpdtc) and (4tfmpq)2Ir(Czdtc) (4tfmpq = 4-(4-(trifluoromethyl)phenyl)quinazoline, dipdtc = N,N-diisopropyl dithiocarbamate, dpdtc = N,N-diphenyl dithiocarbamate, and Czdtc = N-carbazolyl dithiocarbamate) containing the unique four-membered Ir-S-C-S backbone ring were synthesized in five minutes at room temperature with good yields, and the Gibbs free energy calculation results indicate that all reactions are exothermic and thermodynamically favorable processes. The emission colors (λ peak = 641-611 nm), photoluminescence quantum efficiencies (Φ P = 58.3-93.0%) and bipolar properties can be effectively regulated by introducing different electron-donating substituents into the dithiocarbamate ancillary ligands. Employing these emitters, organic light emitting diodes (OLEDs) with double emissive layers exhibit excellent performances with a maximum brightness over 60 000 cd m-2, a maximum current efficiency of 40.68 cd A-1, a maximum external quantum efficiency (EQEmax) of 30.54%, and an EQE of 26.79% at the practical luminance of 1000 cd m-2. These results demonstrate that Ir(iii) complexes with sulfur-containing ligands can be rapidly synthesized at room temperature, which is key to the production of metal luminescent materials for large-scale application in highly efficient OLEDs.

Density functional theory studies on the structures and vibrational spectroscopic characteristics of nickel, copper and zinc naphthalocyanines.

Molecular geometries and electronic structure calculations and spectroscopic assignments for metallonaphthalocyanines NiNc, CuNc and ZnNc are performed on the optimized geometries at B3LYP/6-31G* level. The order of the bond lengths of the NM bonds is computed to be NiNc < CuNc < ZnNc. The Mulliken atomic charges of the central M vary in the same order as the bond lengths. The metal dependent frequencies of IR-active and Raman-active vibrational modes decrease in the order of NiNc > CuNc > ZnNc, which is in reverse sequence of the N-M bond length. The strongest Raman lines predicted at 1570 (NiNc), 1546 (CuNc) and 1525 (ZnNc) cm-1 are highly sensitive to the metal ion. Comparison of the calculated properties of MNcs and MPcs (metallphthalocyanines) compounds reveals some interesting and meaningful results due to extension of the conjugated systems.

Correction: Rapid room temperature synthesis of red iridium(iii) complexes containing a four-membered Ir-S-C-S chelating ring for highly efficient OLEDs with EQE over 30.

[This corrects the article DOI: 10.1039/C8SC05605F.].

The region-specific segregation and catalytic activity of gold-silver nanoparticles.

A coordination-dependent interatomic interaction is developed for Ag atoms and applied to simulate Au-Ag alloy nanoparticles ranging from 2.6 nm to 4.5 nm with different composition ratios at different temperatures. Using the lattice Monte Carlo method, we found that Ag segregates preferentially on the edges of nanoparticles. The region-specific segregation is revealed by our study and it indicates that depending on the size and temperature, Ag may segregate on the edges only, whereas its segregation on flat surfaces is negligible. Assuming that equimolar Au-Ag edges maximize the catalytic synergistic effect, we predicted the size and composition dependent optimal reaction temperatures for the Au-Ag catalyzed CO oxidation, which agree with the available experimental temperatures very well. Our study demonstrates that the region-specific segregation not only complements the details of surface segregation but is also very important for the utilization and tuning of faceted bimetallic nanoparticles.

Computational insight into the cooperative role of non-covalent interactions in the aza-Henry reaction catalyzed by quinine derivatives: mechanism and enantioselectivity.

Density functional theory (DFT) calculations were performed to elucidate the mechanism and the origin of the high enantioselectivity of the aza-Henry reaction of isatin-derived N-Boc ketimine catalyzed by a quinine-derived catalyst (QN). The C-C bond formation step is found to be both the rate-determining and the stereo-controlled step. The results revealed the important role of the phenolic OH group in pre-organizing the complex of nitromethane and QN and stabilizing the in situ-generated nitronate and protonated QN. Three possible activation modes for C-C bond formation involving different coordination patterns of catalyst and substrates were studied, and it was found that both the ion pair-hydrogen bonding mode and the Brønsted acid-hydrogen bonding mode are viable, with the latter slightly preferred for the real catalytic system. The calculated enantiomeric excess (ee) favouring the S enantiomer is in good agreement with the experimental result. The high reactivity and enantioselectivity can be ascribed to the cooperative role of the multiple non-covalent interactions, including classical and non-classical H bonding as well as anionπ interactions. These results also highlight the importance of the inclusion of dispersion correction for achieving a reasonable agreement between theory and experiment for the current reaction.

The Lattice Kinetic Monte Carlo Simulation of Atomic Diffusion and Structural Transition for Gold.

For the kinetic simulation of metal nanoparticles, we developed a self-consistent coordination-averaged energies for Au atoms based on energy properties of gold bulk phases. The energy barrier of the atom pairing change is proposed and holds for the microscopic reversibility principle. By applying the lattice kinetic Monte Carlo simulation on gold films, we found that the atomic diffusion of Au on the Au(111) surface undergoes a late transition state with an energy barrier of about 0.2 eV and a prefactor between 40~50 Å(2)/ps. This study also investigates the structural transition from spherical to faceted gold nanoparticles upon heating. The temperatures of structural transition are in agreement with the experimental melting temperatures of gold nanoparticles with diameters ranging from 2 nm to 8 nm.

Density Functional and Kinetic Monte Carlo Study of Cu-Catalyzed Cross-Dehydrogenative Coupling Reaction of Thiazoles with THF.

Cu-catalyzed cross-dehydrogenative coupling (CDC) reaction of thiazoles with THF has been studied with the density functional theory method and kinetic Monte Carlo (kMC) simulations. Our results show that the previously proposed concerted metalation-deprotonation mechanism is unfavorable. On the basis of the DFT calculation and kMC simulation results, a new mechanism is proposed. In the favorable mechanism, the Cu(II) catalyst first combines with the thiazoles, forming an organocopper species that then binds to the THF radical. The rate-limiting step, C-C bond formation, is realized through an intramolecular structural rearrangement. The Cu catalyst works as a matchmaker to render the C-C bond formation. Kinetic Monte Carlo simulations demonstrate that one should be careful with the conclusions drawn simply from the calculated barriers.

Where does methanol lose hydrogen to trigger steam reforming? A revisit of methanol dehydrogenation on the PdZn alloy model obtained from kinetic Monte Carlo simulations.

Pd/ZnO is a promising catalyst studied for methanol steam reforming (MSR) and the 1 : 1 PdZn alloy is demonstrated to be the active component. It is believed that MSR starts from methanol dehydrogenation to methoxy. Previous studies of methanol dehydrogenation on the ideal PdZn(111) surface show that methanol adsorbs weakly on the PdZn(111) surface and it is hard for methanol to transform into methoxy because of the high dehydrogenation barrier, indicating that the catalyst model is not appropriate for investigating the first step of MSR. Using the model derived from our recent kinetic Monte Carlo simulations, we examined the process CH3OH → CH3O → CH2O → CHO → CO. Compared with the ideal model, methanol adsorbs much more strongly and the barrier from CH3OH → CH3O is much lower on the kMC model. On the other hand, the C-H bond breaking of CH3O, CH2O and CHO becomes harder. We show that co-adsorbed water is important for refreshing the active sites. The present study shows that the first MSR step most likely takes place on three-fold hollow sites formed by Zn atoms, and the inhomogeneity of the PdZn alloy may exert significant influences on reactions.

Density functional study of methanol decomposition on clean and O or OH adsorbed PdZn(111).

Methanol is the future and clean fuel, and its chemistry on metal surfaces has received much attention. In this paper we explore methanol dissociation on the clean and O or OH covered PdZn(111) that mimics Pd∕ZnO catalyst studied as a promising catalyst for methanol steam reforming, using density functional theory at PW91 level and slab model. Our study demonstrates that unlike the situation on Pd (111), methanol preferentially undergoes the O-H bond scission on the PdZn (111). The presence of O and OH species hinders the C-H bond dissociation, but significantly reduces the O-H bond-breaking barrier. The present results indicate that in the course of methanol steam reforming, methanol first loses the hydrogen atom of the hydroxyl group, forming methoxy. This step is greatly enhanced when there are O and∕or OH species (i.e., after water dissociation happens). Analyses reveal that CH2O is formed mainly from CH3O, not from CH2OH.