Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni,Pd)-Based Transition-Metal Catalysts

The performance of 23 density functionals, including one LDA, four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid meta-GGAs, and ten double-hybrid functionals, was investigated for the computation of activation energies of various covalent main-group single bonds by four catalysts: Pd, PdCl⁻, PdCl₂, and Ni (all in the singlet state). A reactant complex, the barrier, and reaction energy were considered, leading to 164 energy data points for statistical analysis. Extended Gaussian AO basis sets were used in all calculations. The best functional for the complete benchmark set relative to estimated CCSD(T)/CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal mol⁻¹ followed by PW6B95-D3, the double hybrid PWPB95-D3, and B3LYP-D3 (1.9 kcal mol⁻¹ each). The other tested hybrid meta-GGAs perform less well (M06-HF: 7.0 kcal mol⁻¹; M06-2X: 6.3 kcal mol⁻¹; M06: 4.9 kcal mol⁻¹) for the investigated reactions. In the Ni case, some double hybrids show larger errors due to partial breakdown of the perturbative treatment for the correlation energy in cases with difficult electronic structures (partial multi-reference character). Only double hybrids either with very low amounts of perturbative correlation (e.g., PBE0-DH) or that use the opposite-spin correlation component only (e.g., PWPB95) seem to be more robust. We also investigated the effect of the D3 dispersion correction. While the barriers are not affected by this correction, significant and mostly positive results were observed for reaction energies. Furthermore, six very recently proposed double-hybrid functionals were analyzed regarding the influence of the amount of Fock exchange as well as the type of perturbative correlation treatment. According to these results, double hybrids with <50–60 % of exact exchange and ∼30 % perturbative correlation perform best.     

Hartree-Fock and Density Functional Theory Studies on the Molecular Recognition of the Cyclodextrin

Cyclodextrins(CD),cyclicoligomersofa-D-glucoseconnectedthroughglycosidical,4bonds,areseductivemoleculesappealingtoresearchersinpureacademicfieldsandappliedtechnologies'.ModelstudiesonCDinclusioncomplexationofferimportantinsightsintoenZyme-substrateillteractionsZandhenceattractgreatattention.Althoughmanyexperimentalapproachesareavailable,molecularmodelingprovidesanimportantalternativewayinstudyingtheCDchemistry'.Duetotheirlargesize,mosttheoreticalstudiesonCDchosemolecularmechanicsormolecula…     

乙炔基自由基C2H与氧气反应的密度泛函理论研究

应用量子化学从头算和密度泛函理论(DFT)对C₂H自由基和O₂的反应进行了研究.在B3LYP/6-311G^**水平上优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,并计算出它们的振动频率和零点振动能(ZPVE).各物种的总能量由CCSD(T)/6-311G^**//B3LYP/6-311G^**给出,并对能量进行了零点能校正.计算结果表明,反应物中自由基C₂H中的边端C进攻O₂形成了中间体1 (HCCOO),中间体1是一个加合产物.由中间体1经过不同的反应通道可以生成不同的产物P₁ (HCO+CO), P₂ (HCCO+O), P₃(CO₂+CH), P₄ (C₂O+OH)和P₅ (2CO+H).反应通道之间存在着竞争机制.其中P₁, P₂是主要产物,其次还有一定比例的P5生成,而产物P₃, P₄的生成几率较低.各条反应通道化学反应热的计算与实验吻合较好.     

Exchange coupling constants using Density Functional Theory: long-range corrected functionals

The presence of the self-interaction error in common exchange correlation functionals is one of the main drawbacks of calculations using methods based on Density Functional Theory. The recently developed asymptotically corrected functionals are an alternative way to solve the erroneous asymptotic behavior that may be at the origin of this error. Some of these functionals have been tested in the calculation of the exchange coupling constants of a number of transition metal dinuclear complexes. Correction of the self-interaction error is especially important for calculations of this kind. This error leads to the inclusion of some static correlation terms that overestimate the J values. This study shows that some of the long-range corrected functionals tested give promising results.     

Ethanol Oxidation Reaction Mechanism on Gold Nanowires from Density Functional Theory

Thin gold nanowires (NWs) are materials that could be used as support in different chemical reactions. Using density functional theory (DFT) it was shown that NWs that form linear atomic chains (LACs) are suitable for stimulating chemical reactions. To this end, the oxidation reaction of ethanol supported on the LACs of Au−NWs was investigated. Two types of LACs were used for the study, one pure and the other with an oxygen impurity. The results showed that the oxygen atom in the LAC fulfills important functions throughout the reaction pathway. Before the chemical reaction, it was observed that the LAC with impurity gains structural stability, that is, the oxygen acts as an anchor for the gold atoms in the LAC. In addition, the LAC was shown to be sensitive to disturbances in its vicinity, which modifies its nucleophilic character. During the chemical reaction, the oxidation of ethanol occurs through two different reaction paths and in two stages, both producing acetaldehyde (CH₃CHO). The different reaction pathways are a consequence of the presence of oxygen in the LAC (oxygen conditions the formation of reaction intermediates). In addition, the oxygen in the LAC also modifies the kinetic behavior in both reaction stages. It was observed that, by introducing an oxygen impurity in the LAC, the activation energy barriers decrease ∼69 % and ∼97 % in the first and second reaction stages, respectively.     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Density Functional Theory Study of Water Diffusion and Clustering on Pd（111）

1 INTRODUCTION The interaction of water molecules with metal sur- faces plays a vital role in a number of important pro- cesses, such as corrosion, heterogeneous catalysis, electrochemical processes in aqueous solutions, hydrogen production, etc.[1] The structure and pro- perties of water adsorbed on well-defined metal sur- faces have been the subject of numerous experi- mental and theoretical investigations. There have been a number of experimental studies of water on metal surfaces throu…     

Ab initio static and molecular dynamics study of 4-styrylpyridine.

We report an in‐depth theoretical study of 4‐styrylpyridine in its singlet S₀ ground state. The geometries and the relative stabilities of the trans and cis isomers were investigated within density functional theory (DFT) as well as within Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and coupled cluster (CC) theories. The DFT calculations were performed using the B3LYP and PBE functionals, with basis sets of different qualities, and gave results that are very consistent with each other. The molecular structure is thus predicted to be planar at the energy minimum, which is associated with the trans conformation, and to become markedly twisted at the minimum of higher energy, which is associated with the cis conformation. The results of the calculations performed with the post‐HF methods approach those obtained with the DFT methods, provided that the level of treatment of the electronic correlation is high enough and that sufficiently flexible basis sets are used. Calculations carried out within DFT also allowed the determination of the geometry and the energy of the molecule at the biradicaloid transition state associated with the thermal cis?trans isomerization and at the transition states associated with the enantiomerization of the cis isomer and with the rotations of the pyridinyl and phenyl groups in the trans and cis isomers. Car–Parrinello molecular dynamics simulations were also performed at 50, 150, and 300 K using the PBE functional. The studies allowed us to evidence the highly flexible nature of the molecule in both conformations. In particular, the trans isomer was found to exist mainly in a nonplanar form at finite temperatures, while the rotation of the pyridinyl ring in the cis isomer was incidentally observed to take place within ≈1 ps during the simulation carried out at 150 K on this isomer.     

The Stability of α‐Hydroperoxyalkyl Radicals

High‐level ab initio and Born–Oppenheimer molecular dynamic calculations have been carried out on a series of hydroperoxyalkyl (α‐QOOH) radicals with the aim of investigating the stability and unimolecular decomposition mechanism into QO+OH of these species. Dissociation was shown to take place through rotation of the C?O(OH) bond rather than through elongation of the CO?OH bond. Through the C?O(OH) rotation, the unpaired electron of the radical overlaps with the electron density on the O?OH bond, and from this overlap the C=O π bond forms and the O?OH bond breaks spontaneously. The CH₂OOH, CH(CH₃)OOH, CH(OH)OOH, and α‐hydroperoxycycloheptadienyl radical were found to decompose spontaneously, but the CH(CHO)OOH has a decomposition energy barrier of 5.95 kcal mol^?1 owing to its steric and electronic features. The systems studied in this work provide the first insights into how structural and electronic effects govern the stabilizing influence on elusive α‐QOOH radicals.     

Relativistic Density Functional NMR Tensors Analyzed with Spin-free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the ¹⁵N shielding in a heavy metal complex.     

Time‐Dependent Density Functional Theory Molecular Dynamics Simulations of Liquid Water Radiolysis

The early stages of the Coulomb explosion of a doubly ionized water molecule immersed in liquid water are investigated with time‐dependent density functional theory molecular dynamics (TD–DFT MD) simulations. Our aim is to verify that the double ionization of one target water molecule leads to the formation of atomic oxygen as a direct consequence of the Coulomb explosion of the molecule. To that end, we used TD–DFT MD simulations in which effective molecular orbitals are propagated in time. These molecular orbitals are constructed as a unitary transformation of maximally localized Wannier orbitals, and the ionization process was obtained by removing two electrons from the molecular orbitals with symmetry 1B₁, 3A₁, 1B₂ and 2A₁ in turn. We show that the doubly charged H₂O²⁺ molecule explodes into its three atomic fragments in less than 4 fs, which leads to the formation of one isolated oxygen atom whatever the ionized molecular orbital. This process is followed by the ultrafast transfer of an electron to the ionized molecule in the first femtosecond. A faster dissociation pattern can be observed when the electrons are removed from the molecular orbitals of the innermost shell. A Bader analysis of the charges carried by the molecules during the dissociation trajectories is also reported.     

Band‐Gap Manipulations of Monolayer Graphene by Phenyl Radical Adsorptions: A Density Functional Theory Study

Phenyl radical (Ph^.) adsorption on monolayer graphene sheets is used to investigate the band‐gap manipulation of graphene through density functional theory. Adsorption of a single Ph^. on graphene breaks the aromatic π‐bond and generates an unpaired electron, which is delocalized to the ortho or para position. Adsorption of a second radical at the ortho or para position saturates the radical by electron pairing and results in semiconducting graphene. Adsorption of a second radical at the ortho position (ortho–ortho pairing) is found to be more favorable than adsorption at the para position (ortho–para pairing), and the ortho–ortho pairing has stronger effects on band‐gap opening compared with ortho–para pairing. Adsorption of even numbers of Ph^. on graphene by ortho–ortho and ortho–para pairings, in general, increases the band gap. Our study shows promise of band‐gap manipulation in monolayer graphene by Ph^. adsorption, leading to potential wider applications of graphene.     

Understanding the Role of Water in Aqueous Ruthenium‐Catalyzed Transfer Hydrogenation of Ketones.

We report an accurate computational study of the role of water in transfer hydrogenation of formaldehyde with a ruthenium‐based catalyst using a water‐specific model. Our results suggest that the reaction mechanism in aqueous solution is significantly different from that in the gas phase or in methanol solution. Previous theoretical studies have shown a concerted hydride and proton transfer in the gas phase (M. Yamakawa, H. Ito, R. Noyori, J. Am. Chem. Soc. 2000 , 122, 1466–1478;J.‐W. Handgraaf, J. N. H. Reek, E. J. Meijer, Organometallics 2003 , 22, 3150–3157; D. A. Alonso, P. Brandt, S. J. M. Nordin, P. G. Andersson, J. Am. Chem. Soc. 1999 , 121, 9580–9588; D. G. I. Petra, J. N. H. Reek, J.‐W. Handgraaf, E. J. Meijer, P. Dierkes, P. C. J. Kamer, J. Brussee, H. E. Schoemaker, P. W. N. M. van Leeuwen, Chem. Eur. J. 2000 , 6, 2818–2829), whereas a delayed, solvent‐mediated proton transfer has been observed in methanol solution (J.‐W. Handgraaf, E. J. Meijer, J. Am. Chem. Soc. 2007 , 129, 3099–3103). In aqueous solution, a concerted transition state is observed, as in the previous studies. However, only the hydride is transferred at that point, whereas the proton is transferred later by a water molecule instead of the catalyst.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

New Insights into the Jahn–Teller Effect through ab initio Quantum-Mechanical/Molecular-Mechanical Molecular Dynamics Simulations of CuII in Water

The Cu^II hydration shell structure has been studied by means of classical molecular dynamics (MD) simulations including three-body corrections and hybrid quantum-mechanical/molecular-mechanical (QM/MM) molecular dynamics (MD) simulations at the Hartree–Fock level. The copper(II ) ion is found to be six-fold coordinated and [Cu(H₂O)₆]²⁺ exhibits a distorted octahedral structure. The QM/MM MD approach reproduces correctly the experimentally observed Jahn–Teller effect but exhibits faster inversions (<200 fs) and a more complex behaviour than expected from experimental data. The dynamic Jahn–Teller effect causes the high lability of [Cu(H₂O)₆]²⁺ with a ligand-exchange rate constant some orders of magnitude higher than its neighbouring ions Ni^II and Zn^II. Nevertheless, no first-shell water exchange occurred during a 30-ps simulation. The structure of the hydrated ion is discussed in terms of radial distribution functions, coordination numbers, and various angular distributions and the dynamical properties as librational and vibrational motions and reorientational times were evaluated, which lead to detailed information about the first hydration shell. Second-shell water-exchange processes could be observed within the simulation time scale and yielded a mean ligand residence time of ≈20 ps.