First‐principles study of molecular hydrogen dissociation on doped Al12X (X = B,Al, C,Si, P,Mg, and Ca) clusters

Inspired by the concept of superatom via substitutionally doping an Al₁₃ magic cluster, we investigated the H₂ molecule dissociation on the doped icosahedral Al₁₂X (X = B, Al, C, Si, P, Mg, and Ca) clusters by means of density functional theory. The computed reaction energies and activation barriers show that the concept of superatom is still valid for the catalysis behavior of doped metal clusters. The hydrogen dissociation behavior on metal clusters characterized by the activation barrier and reaction energy can be tuned by controllable doping. Thus, doped Al₁₂X clusters might serve as highly efficient and low‐cost catalysts for hydrogen dissociation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Type‐I Dyotropic Reactions: Understanding Trends in Barriers

To understand the factors that control the activation barrier of type‐I 1,2‐dyotropic reactions (X‐EH₂‐CH₂‐X*→X*‐EH₂‐CH₂‐X, with E=C and Si, X=H, CH₃, SiH₃, F to I) and trends therein as a function of the migrating groups X, we have explored ten archetypal model reactions of this class using relativistic density functional theory (DFT) at ZORA‐OLYP/TZ2P. The main trends in reactivity are rationalized using the activation strain model of chemical reactivity, which had to be extended from bimolecular to unimolecular reactions. Thus, the above type‐I dyotropic reactions can be conceived as a relative rotation of the CH₂CH₂ and [X???X] fragments in X‐CH₂‐CH₂‐X. The picture that emerges from these analyses is that reduced C? X bonding in the transition state is the origin of the reaction barrier. Also the trends in reactivity on variation of X can be understood in terms of how sensitive the C? X interaction is towards adopting the transition‐state geometry. A valence bond analysis complements the analyses and confirms the picture emerging from the activation strain model.     

Ultrafast dynamics in noble metal clusters: The role of internal vibrational redistribution

We present a theoretical study of the ultrafast dynamics in noble metal clusters interacting with molecular oxygen which is of fundamental importance for the understanding and design of cluster-based heterogenous nanocatalysts. We demonstrate that intrinsic dynamical properties can significantly promote the reactivity of small noble metal clusters towards O₂. This concept is illustrated by performing collision simulations between and clusters and O₂ in the framework of the ab initio molecular dynamics (MD) using density functional theory (DFT). We show that different nature and efficiency of the internal vibrational energy redistribution (IVR) during the collisions with O₂ are responsible for considerably different sticking probabilities of O₂ to silver and gold clusters, respectively. In the case of , resonant IVR between the cluster and the O₂ subunit activates the O–O bond and promotes the subsequent oxidation reaction. In contrast, in the case of fast dissipative IVR on the time scale of 1 ps leads to a higher sticking probability for O₂ but the O–O bond is very rapidly deactivated and cannot participate in further oxidation processes. These findings allow us to introduce the nature of IVR as a criterion for promoting the reactivity of noble metal clusters. Such different behaviour of silver and gold clusters colliding with O₂ originates from difference in relativistic effects which are considerably more pronounced in the case of gold clusters causing more directional rigid bonding in contrast to silver clusters with more s-metallic floppy character. Moreover, we demonstrate that breaking of O–O bond can be induced in by a selective excitation of the O–O bond with an ultrashort pulse in the infrared spectral range. This opens the perspective to control the action of nanocatalysts by employing shaped laser pulses and thus bridges the fields of femtochemistry and cluster nanocatalysis.     

Computational Study of the Interactions between Benzene and Crystalline Ice Ih: Ground and Excited States

Ground‐state geometries of benzene on crystalline ice cluster model surfaces (I_h) are investigated. It is found that the binding energies of benzene‐bound ice complexes are sensitive to the dangling features of the binding sites. We used time‐dependent DFT to study the UV spectroscopy of benzene, ice clusters, and benzene–ice complexes, by employing the M06‐2X functional. It is observed that the size of the ice cluster and the dangling features have minor effects on the UV spectral characteristics. Benzene‐mediated electronic excitations of water towards longer wavelengths (above 170 nm) are noted in benzene‐bound ice clusters, where the cross‐section of photon absorption by water is negligible, in good agreement with recent experimental results (Thrower et al., J. Vac. Sci. Technol. A, 2008, 26 , 919–924). The intensities of peaks associated with water excitations in benzene–ice complexes are found to be higher than in isolated ice clusters. The π→π* electronic transition of benzene in benzene–ice complexes undergoes a small redshift compared with the isolated benzene molecule, and this holds for all benzene‐bound ice complexes.     

DFT calculations on the catalytic oxidation of CO over Si-doped (6,0) boron nitride nanotubes

The oxidation of carbon monoxide (CO) is important for a series of technological and environmental applications. In this work, the catalytic oxidation of CO on Si-doped (6,0) boron nitride nanotubes (BNNTs) is investigated by using density functional theory calculations. Reaction barriers and corresponding thermodynamic parameters were calculated using the M06-2X, B3LYP and wB97XD density functionals with 6-31G* basis set. Our results indicate that a vacancy defect in BNNT strongly stabilizes the Si adatom and makes it more positively charged. This charging enhances the adsorption of reaction gases (O₂ and CO) and results in the change of the electronic structure properties of the tube. The calculated barrier of the reaction CO + O₂ → CO₂ + O_ads on Si-doped BNNTs following the Langmuir–Hinshelwood is lower than that on the traditional noble metal catalysts. The second step of the oxidation would be the Eley–Rideal reaction (CO + O_ads → CO₂) with an energy barrier of about 1.8 and 10.1 kcal/mol at M06-2X/6-31G* level. This suggests that the CO oxidation catalyzed by the Si-doped BNNTs is likely to occur at the room temperature. The results also demonstrate that the activation energies and thermodynamic quantities calculated by M06-2X, B3LYP and wB97XD functionals are consistent with each other.     

O2 dissociative adsorption on the Cu‐, Ag‐, and W‐doped Al(111) surfaces from DFT computation

The O₂ adsorption and dissociation on M‐doped (M = Cu, Ag, W) Al(111) surface were studied by density functional theory. The adsorption energy of adsorbate, the average binding energy and surface energy of Al surface, and the doping energy of doping atom were calculated. All the doped atoms can be stably combined with Al atoms, while being slightly embedded in the surface to a certain depth. The TOP‐type surfaces are the most stable doped surfaces for O₂ adsorption, which is related to the orbital hybridization between the adsorbate and the surface atoms, the electronegativity, and the orbital energy level of the doping atoms. Moreover, the O atoms and doping atoms contribute significantly to the density of states (DOS), especially the O‐p orbital electrons and the d orbital electrons of doping atoms. The degree of O₂ dissociation is related to the doping atoms on Al surfaces, and the doping atoms actually resist the dissociation of O₂. W atoms have the best resistance effect on the O₂ dissociation as compared with Cu and Ag atoms, especially W‐1NN surface, which has both large barrier energy and reaction energy.     

A Theoretical Study on the Mechanism of C2H4 Oxidation over a Neutral V3O8 Cluster

Density functional theory (DFT) calculations are used to investigate the reaction mechanism of V₃O₈+C₂H₄. The reaction of V₃O₈ with C₂H₄ produces V₃O₇CH₂+HCHO or V₃O₇+CH₂OCH₂ overall barrierlessly at room temperature, whereas formation of hydrogen‐transfer products V₃O₇+CH₃CHO is subject to a tiny overall free energy barrier (0.03 eV), although the formation of the last‐named pair of products is thermodynamically more favorable than that of the first two. These DFT results are in agreement with recent experimental observations. The (O_b)₂V(O_tO_t)^. (b=bridging, t=terminal) moiety containing the oxygen radical in V₃O₈ is the active site in the reaction with C₂H₄. Similarities and differences between the reactivities of (O_b)₂V(O_tO_t)^. in V₃O₈ and the small VO₃ cluster [(O_t)₂VO_t^.] are discussed. Moreover, the effect of the support on the reactivity of the (O_b)₂V(O_tO_t)^. active site is evaluated by investigating the reactivity of the cluster VX₂O₈, which is obtained by replacing the V atoms in the (O_b)₃VO_t support moieties of V₃O₈ with X atoms (X=P, As, Sb, Nb, Ta, Si, and Ti). Support X atoms with different electronegativities influence the oxidative reactivity of the (O_b)₂V(O_tO_t)^. active site through changing the net charge of the active site. These theoretical predictions of the mechanism of V₃O₈+C₂H₄ and the effect of the support on the active site may be helpful for understanding the reactivity and selectivity of reactive O^. species over condensed‐phase catalysts.     

Bis(sulfonylimide)ruthenium(VI) Porphyrins: X‐ray Crystal Structure and Mechanism of C?H Bond Amination by Density Functional Theory Calculations

The X‐ray crystal structure of [Ru^VI(NMs)₂(tmp)] (Ms=SO₂‐ p‐MeOC₆H₄; tmp=5,10,15,20‐tetramesitylporphyrinato(2?)), a metal sulfonylimide complex that can undergo alkene aziridination and C? H bond amination reactions, shows a Ru?N distance of 1.79(3) Å and Ru‐N‐S angle of 162.5(3)°. Density functional theory (DFT) calculations on the electronic structures of [Ru^VI(NMs)₂(tmp)] and model complex [Ru^VI(NMs)₂(por⁰)] (por⁰=unsubstituted porphyrinato(2?)) using the M06L functional gave results in agreement with experimental observations. For the amination of ethylbenzene by the singlet ground state of [Ru^VI(NMs)₂(por⁰)], DFT calculations using the M06L functional revealed an effectively concerted pathway involving rate‐limiting hydrogen atom abstraction without a distinct radical rebound step. The substituent effect on the amination reactivity of ethylbenzene by [Ru^VI(NX)₂(por⁰)] (X=SO₂‐p‐YC₆H₄ with Y=MeO, Me, H, Cl, NO₂) was examined. Electron‐withdrawing Y groups lower the energy of the LUMOs of [Ru^VI(NX)₂(por⁰)], thus facilitating their interaction with the low‐lying HOMO of the ethylbenzene C? H bond and hence increasing the reactivity of [Ru^VI(NX)₂(por⁰)]. DFT calculations on the amination/aziridination reactions of [Ru^VI(NSO₂C₆H₅)₂(por⁰)] with pent‐4‐enal, an aldehyde substrate bearing acyl, homoallylic, and allylic C? H bonds and a C?C bond, revealed a lower reaction barrier for the amination of the acyl C? H bond than for both the amination of the other C? H bonds and aziridination of the C?C bond in this substrate.     

O?O bond homolysis in hydrogen peroxide

O O bond homolysis in hydrogen peroxide (H₂O₂) has been studied using theoretical methods of four conceptually different types: hybrid DFT (B3LYP, M06‐2X), double‐hybrid DFT (B2‐PLYP), coupled‐cluster (CCSD(T)), and multiconfigurational (CASPT2). In addition, the effects of basis set size have also been analyzed. For all of these methods, the O O bond homolysis in hydrogen peroxide has been found to proceed through hydrogen bonded radical pair complexes. Reaction barriers for collapse of the radical pairs to hydrogen peroxide are minute, leading to an overall very flat potential energy surface. However, hydrogen bonding energies in the radical pair complex expressed as the energy difference to two separate hydroxyl radicals are sizeable and exceed 10 kJ/mol for all theoretical methods considered in this study. © 2017 Wiley Periodicals, Inc.     

Pt原子在γ-Al2O3(001)表面的吸附及迁移

采用密度泛函理论(DFT)中广义梯度近似(GGA)方法, 对Pt原子与γ-Al₂O₃(001)面的相互作用及迁移性能进行了研究. 分析了各种可能吸附位及吸附构型的松弛和变形现象, 吸附能和迁移能垒的计算结果表明: Pt团簇能够稳定吸附在该表面. Pt原子在表面O位的吸附能明显较高, 这主要是由Pt向基底O原子转移了电子所致. 电荷布居分析表明, Pt原子显电正性, Pt和Al原子之间存在排斥作用, 导致与Al原子产生较弱相互作用. 计算的平均吸附能大小依赖于Pt团簇的大小和形状, 总体趋势是随着Pt原子数增多, 吸附能降低. Pt原子在γ-Al₂O₃(001)表面迁移过程所需克服的迁移能垒最高值为0.51 eV. 随着吸附的Pt原子数增多,更倾向于形成Pt团簇. 因此, Pt原子在γ-Al₂O₃(001)表面的吸附演变不可能形成光滑、均匀平铺的吸附构型, 而在一定条件下容易出现团聚.     

Assessment of dispersion‐improved exchange‐correlation functionals for the simulation of CO2 binding by alcoholamines

In this study, 12 bound complexes were selected to construct a database for testing 15 dispersion‐improved exchange‐correlation (XC) functionals, including hybrid generalized gradient approximation (GGA), modified using the Grimme's pairwise strategy, and double hybrid XC functionals, for specifically characterizing the CO₂ binding by alcoholamines. Bound complexes were selected based on the characteristics of their hydrogen bonds, dispersion, and electrostatic (particularly between the positive charge of CO₂ and the lone pair of N of alcoholamines) interactions. The extrapolated binding energy from the aug‐cc‐pVTZ (ATZ) to aug‐cc‐pVQZ (AQZ) basis set at the CCSD(T)/CBS(MP2+DZ) level was used as the reference for the XC functional comparison. M06‐2X produced the optimal agreement if the optimized geometries at MP2/ATZ level were chosen for all the test bound complexes. However, M06‐L, ωB97X, and ωB97, and were preferred if the corresponding density functional theory (DFT) optimized geometries were adapted for the benchmark. Simple bimolecular reaction between CO₂ and monoethanolamine simulated using polarizable continuum solvation model confirmed that ωB97, ωB97X, and ωB97XD qualitatively reproduced the energetics of MP2 level. The inconsistent performance of the tested XC functionals, observed when using MP2 or DFT optimized geometries, raised concerns regarding using the single‐point ab initio correction combined with DFT optimized geometry, particularly for determining the nucleophilic attack by alcoholamines to CO₂. © 2014 Wiley Periodicals, Inc.     

Reducing the Cost and Preserving the Reactivity in Noble‐Metal‐Based Catalysts: Oxidation of CO by Pt and Al–Pt Alloy Clusters Supported on Graphene

The oxidation mechanisms of CO to CO₂ on graphene‐supported Pt and Pt‐Al alloy clusters are elucidated by reactive dynamical simulations. The general mechanism evidenced is a Langmuir–Hinshelwood (LH) pathway in which O₂ is adsorbed on the cluster prior to the CO oxidation. The adsorbed O₂ dissociates into two atomic oxygen atoms thus promoting the CO oxidation. Auxiliary simulations on alloy clusters in which other metals (Al, Co, Cr, Cu, Fe, Ni) replace a Pt atom have pointed to the aluminum doped cluster as a special case. In the nanoalloy, the reaction mechanism for CO oxidation is still a LH pathway with an activation barrier sufficiently low to be overcome at room temperature, thus preserving the catalyst efficiency. This provides a generalizable strategy for the design of efficient, yet sustainable, Pt‐based catalysts at reduced cost.     

Computational chemistry and molecular simulations of phosphoric acid

Quantum mechanics (QM) calculations, molecular dynamics (MD) simulations using the condensed‐phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field, and the atom‐centered density matrix propagation (ADMP) approach have been used to investigate properties of phosphoric acid (PA). QM using B3LYP/6‐31++G(d,p) density functional theory were used to calculate gas‐phase proton affinities and interaction energies of PA and its derivatives. Detailed single coordinate driving, followed by quadratic synchronous transit optimization was used to determine energy barriers for different proton transfer (PT) pathways. Determined energy barrier heights in ascending order are (unit: kJ/mol): H₃O⁺→H₃PO₄ (0); H₄P₂O₇→H₃PO₄ (2.61); H₃PO₄→H₂PO (5.31); H₄PO→H₃PO₄ (～7.33); H₃PO₄→H₄P₂O₇/H₃PO₄→H₃PO₄ (15.99); H₄P₂O₇→H₂O (28.61); H₃PO₄→H₂O (47.14). The COMPASS force field was used to study condensed‐phase properties of PA. Good agreement between experimental data and MD results including density, radial distribution functions, and self‐diffusion coefficient at different temperatures provides validation of the COMPASS force field for PA. Finally, preliminary ADMP studies on a cluster of three PA molecules shows that the ADMP approach can reasonably describe the PT and self‐dissociation processes in PA. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Electronic structure and stability of binary and ternary aluminum‐bismuth‐nitrogen nanoclusters

The electronic structure and stability in binary and ternary aluminum‐bismuth‐nitrogen nanoclusters up to six atoms are studied using density functional theory (DFT). The lowest energy geometries were obtained by sampling the geometrical space with a Monte Carlo method and geometry optimizations, at DFT level, with M06L functional. The clusters stability is analyzed using formation and fragmentation energies. Our results show that a high concentration of nitrogen presents a tendency to form nitrogen clusters. highest occupied molecular orbital‐lowest unoccupied molecular orbital gaps show the well‐known oscillation as the number of atoms is increased. Bonding between Al, Bi, and N has mainly a π character. Bismuth and aluminum atoms tend to promote high multiplicity states in small clusters. These new binary and ternary materials provide a potential new field in optoelectronics and high energetic material compounds. © 2014 Wiley Periodicals, Inc.