Comparative study of BSSE correction methods at DFT and MP2 levels of theory

A comparative study of intermolecular potential energy curves is performed on the complexes H₂O(SINGLE BOND)HF, H₂O(SINGLE BOND)H₂O, H₂O(SINGLE BOND)H₂S, and H₂S(SINGLE BOND)H₂S using nine different basis sets at the MP2 and DFT (BLYP and B3LYP) levels of theory. The basis set superposition error is corrected by means of the counterpoise scheme and based on the “chemical Hamiltonian approach.” The counterpoise and CHA-corrected DFT curves are generally close to each other. Using small basis sets, the B3LYP functional cannot be favored against the BLYP one because the BLYP results sometimes get closer to the MP2 values than those of B3LYP. From the results—including the available literature data—we suggest that one has to use at least polarized-valence triple-zeta-quality basis sets (TZV, 6-311G) for the investigation of hydrogen-bonded complexes. Special attention must be paid to the physical nature of the binding. If the dispersion forces become significant DFT methods are not able to describe the interaction. Proper correction for the basis set superposition error is found to be mandatory in all cases. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 575–584, 1998     

Density functional study of chlorine–oxygen compounds related to the ClO self-reaction

Geometrical parameters, vibrational frequencies, relative stabilities, and dissociation energies of the three stable Cl₂O₂ isomers and the OClO and ClOO radicals were investigated by density functional theory (DFT). The present analysis shows that DFT using hybrid functionals is capable of describing these systems to at least the same degree of accuracy as ab initio methods. The average absolute bond-length deviation of ClClO₂, ClOOCl, and ClO₂ from experimental results is 0.024/0.027 Å, with a maximum deviation for the dichlorine peroxide O(SINGLE BOND)O bond equal to 0.072/0.063 Å, for the B3PW91 and B3LYP functionals, respectively. The average absolute bond-angle deviation for the hybrid functionals is 0.8°. Harmonic vibrational frequencies calculated with DFT give for all Cl(SINGLE BOND)O compounds good agreement with experiments. The dissociation energies of ClOOCl, OClO, and ClOO were found to be in good agreement with experiments, the average error being less than 1.2 kcal/mol. The two isomers chloryl chloride (ClClO₂) and dichlorine peroxide (ClOOCl) were found to be approximately 9 kcal/mol more stable than the chlorine chlorite (ClOClO) isomer. The ClOO isomer is predicted to be 3.0 kcal/mol more stable than OClO, in accordance with the experimental value of 4 kcal/mol. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 203–217, 1998     

Accurate Structure and Bonding Description of the Transition Metal‐Disulfur Monoxide Complexes [(PMe3)2M(S2O)] (M = Ni,Pd, Pt): Grimme Dispersion Corrected DFT Study

Geometry and bonding energy analysis of M–S₂O bonds in the metal‐disulfur monoxide complexes [(PMe₃)₂M(S₂O)] of nickel, palladium, and platinum were investigated at DFT, DFT‐D3, and DFT‐D3(BJ) methods using three different functionals (BP86, PBE, and TPSS). The TPSS/DFT‐D3(BJ) yields better geometry, while the BP86 geometry is least accurate for studied complexes. The geometry of platinum complex optimized at TPSS/DFT‐D3(BJ) level is in excellent agreement with the available experimental values. The M–S bonds are shorter than the M–S(O) bonds. The Mayer bond orders suggest the presence of M–S and M–S(O) single bonds. Both the M–S and M–S(O) bond lengths vary with the density functionals as TPSS‐D3(BJ) < TPSS < PBE < BP86. The Hirshfeld charge distribution indicates that the overall charge flows from metal fragment to [S₂O]. The Ni–S₂O bond has greater degree of covalent character than the ionic. The contribution of dispersion interactions is large in computing accurate bond dissociation energies between the interacting fragments. The BDEs are largest for the functional TPSS and smallest for the functional BP86. The DFT‐D3 dispersion corrections to the BDEs between the metal fragments [(PMe₃)₂M] and ligand fragment [(S₂O)] for the TPSS functional are in the range 7.1–7.3 kcal · mol^–1, which are smaller than the corresponding DFT‐D3(BJ) dispersion corrections (9.4–10.6 kcal · mol^–1).     

A theoretical study on the structure and hygroscopicity of ammonium dinitramide

Density functional theory (DFT) and the dispersion corrected DFT have been used to investigate the hygroscopicity of ammonium dinitramide (ADN). Calculation results show that the gaseous ADN has a strong hydrogen bond. But the ionic pair structure NH₄ ⁺ · N(NO₂)^? is stabilized upon the addition of water molecules. Natural bond orbital calculations suggest that the intra- and intermolecular orbital interactions LP(O) → σ*(N–H) or LP(O) → σ*(O–H) make the system stabilized as a whole. En energy decomposition analysis reveals that the interactions between ADN and H₂O are dominated by the electrostatic and orbital interactions. The formation reactions become more spontaneous with the increasing number of water molecules but can be weakened by the growing temperature from 200 to 400 K. Moreover, the molecular dynamic method is applied to explore a more realistic cluster model to study the interactions between ADN and H₂O.     

The reliability of DFT methods to predict electronic structures and minimum energy crossing point for [FeIVO](OH)2 models: A comparison study with MCQDPT method

In order to test the reliability of DFT methods for calculating electronic structures of [Fe^IVO] system, detailed calculations of [Fe^IVO](OH)₂ models were performed for several low‐energy states using multiconfiguration quasidegenerate perturbation theory (MCQDPT) as well as DFT‐based methods. The minimum energy crossing points (MECP) of ⁵A₁/⁵B₂ and ³B₂/⁵B₂ were investigated based on Lagrange–Newton approach. The results show that M06 functional produce energy gaps close to those of MCQDPT results. Another topic in this article is that the electron configurations of [Fe^IVO](OH)₂ models strongly depend on the type of surface ligand used, and the two lowest states of these can facile transition each other by the MECP. The practicability of M06 method in locating the MECP is validated by the results of MCQDPT which demonstrate the two‐state reactivity (TSR) can be studied with proper DFT method. These inspections provide the basis for further TSR study of larger [Fe^IVO] system. © 2014 Wiley Periodicals, Inc.     

UO22＋•nH2O和PuO22＋•nH2O (n＝2, 4, 5, 6)密度泛函理论研究

首先用密度泛函理论(DFT)方法研究了铀酰和钚酰离子的几何与电子结构, 计算结果与实验基本符合, 表明DFT方法也能用于含铀和钚重原子的化合物计算. 然后对铀酰和钚酰水合离子的几何构型、Mulliken集居数分布以及铀酰(钚酰)与配体水分子的结合能进行计算, 计算结果表明UO₂^2＋•5H₂O和PuO₂^2＋•5H₂O分别为铀酰和钚酰系列水合离子中最稳定的配合物.     

Unexpected Vulnerability of DPEphos to C−O Activation in the Presence of Nucleophilic Metal Hydrides

C−O bond activation of DPEphos occurs upon mild heating in the presence of [Ru(NHC)₂(PPh₃)₂H₂] (NHC=N-heterocyclic carbene) to form phosphinophenolate products. When NHC=IEt₂Me₂, C−O activation is accompanied by C−N activation of an NHC ligand to yield a coordinated N-phosphino-functionalised carbene. DFT calculations define a nucleophilic mechanism in which a hydride ligand attacks the aryl carbon of the DPEphos C−O bond. This is promoted by the strongly donating NHC ligands which render a trans dihydride intermediate featuring highly nucleophilic hydride ligands accessible. C−O bond activation also occurs upon heating cis-[Ru(DPEphos)₂H₂]. DFT calculations suggest this reaction is promoted by the steric encumbrance associated with two bulky DPEphos ligands. Our observations that facile degradation of the DPEphos ligand via C−O bond activation is possible under relatively mild reaction conditions has potential ramifications for the use of this ligand in high-temperature catalysis.     

Substituting CF2 for O4′ in Components of Nucleic Acids: Towards Systems with Reduced Propensity to Form Abasic Lesions

Intrinsic structural features and energetics of nucleotides containing variously fluorinated sugars as potential building blocks of DNA duplexes and quadruplexes are explored systematically using the modern methods of density functional theory (DFT) and quantum chemical topology (QCT). Our results suggest that fluorination at the 2′‐β or 2′‐α,β positions somewhat stabilizes in vacuo the AI relative to the BI conformations. In contrast, substitution of the CF₂ group for the O4′ atom (O4′‐CF₂ modification) leads to a preference of the BI relative to AI DNA‐like conformers. All the studied modifications result in a noticeable increase in the stability of the glycosidic bond [estimated by the relaxed force constants (RFC) approach], with particularly encouraging results for the O4′‐CF₂ derivative. Consequently, the O4′‐CF₂ modified systems are suggested and explored as promising scaffolds for the development of duplex and quadruplex structures with reduced propensity to form abasic lesions and to undergo DNA damage.     

Optimized structures and vibration frequencies of the ether–water complex: a DFT and FTIR study

The infrared spectrum of ether was studied using Fourier transform infrared spectroscopy in conjunction with the density functional theory (DFT). The optimized structures and vibrational frequencies of the ether·(H₂O) _n (n = 1–3) complexes were obtained at B3LYP/6-31G(d) theory levels. Compared to those of free-form ether, the C–O stretching vibrational frequencies of the ether–water complexes are found to shift to red by up to 39 cm^?1 with an increase in the C–O length of 0.016 Å. Meanwhile, the frequency of the O–H stretching modes of water in the complexes appears significantly redshifted to a varying degree. The DFT calculations suggest that these shifts are caused by the hydrogen bonding between ether and water.     

Erratum: Separability between valence and conduction bands in transition metal clusters

Density Functional Theory (DFT) and direct ab initio molecular dynamics (MD) calculations were applied to the hydrogen molecule trapped in a water cluster composed of 12 water molecules (H₂O)₁₂. The static DFT calculation showed that the H₂ molecule is trapped in the center of mass of (H₂O)₁₂. The vibrational frequency of the H–H stretching mode of the H₂ molecule trapped in the water cluster was blueshifted from that in vacuo. On the other hand, the vibrational frequency of H₂ in water‐hydrogen 1:1 complex (H₂O–H₂) was redshifted. A direct ab initio MD calculation of H₂(H₂O)₁₂ at 50 K indicated that the H₂ molecule is rotated freely around the center of mass of the water cluster. The origin of the spectral shift of H₂ in water ice is discussed on the basis of the theoretical results. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Nitrous Oxide as a Hydrogen Acceptor for the Dehydrogenative Coupling of Alcohols

The oxidation of alcohols with N₂O as the hydrogen acceptor was achieved with low catalyst loadings of a rhodium complex that features a cooperative bis(olefin)amido ligand under mild conditions. Two different methods enable the formation of either the corresponding carboxylic acid or the ester. N₂ and water are the only by‐products. Mechanistic studies supported by DFT calculations suggest that the oxygen atom of N₂O is transferred to the metal center by insertion into the Rh?H bond of a rhodium amino hydride species, generating a rhodium hydroxy complex as a key intermediate.     

H2O和OH在UO(100)表面吸附的密度泛函研究

运用密度泛函理论中的广义梯度近似(GCA)的PW91方法结合周期性平板模型,研究了H2O分子和OH在UO(100)表面上的吸附.通过对不同吸附位的吸附能和几何结构参数的计算和比较发现:水分子在UO(100)表面的吸附为化学吸附,水分子平面与UO(100)表面夹角为15°的吸附构型最稳定,吸附能最大,近89 kJ·mol-1.对H2O吸附前后的态密度分析可知,H2O通过其O原子的P轨道与底物U原子的d轨道作用.本文还进一步探讨H2O在UO(100)表面的解离机理.     

Synthesis of an Iron(IV) Aqua–Oxido Complex Using Ozone as an Oxidant

The iron(IV) oxido complex [(tmc)Fe=O(OTf)]OTf with the macrocyclic ligand 1,4,8,11‐tetramethyl‐1,4,8,11‐tetraazacyclo‐tetradecane (tmc) has been synthesized using ozone as an oxidant. By adding water to this compound the complex [(H₂O)(tmc)Fe=O)](OTf)₂ could be prepared. This complex is important in regard to a better understanding of the reactivity of Fe^IV oxido complexes. Mössbauer measurements using the solid compound showed an isomer shift of δ=0.19 mm s^?1 and a quadrupole splitting ΔE_Q=1.38 mm s^?1, confirming the high‐valent Fe^IV state. DFT calculations were performed and led to an assignment of triplet spin multiplicity. Crystallographic characterization of [(H₂O)(tmc)Fe=O)](OTf)₂ as well as of starting materials [(tmc)Fe(CH₃CN)](OTf)₂ and [(tmc)Fe(OTf)]OTf together with previous results strongly suggest that [(H₂O)(tmc)Fe=O)](OTf)₂ was formed similar to the oxido–hydroxido tautomerism analogous to heme systems.     

The Critical Role of Reductive Steps in the Nickel-Catalyzed Hydrogenolysis and Hydrolysis of Aryl Ether C−O Bonds

The hydrogenolysis of the aromatic C−O bond in aryl ethers catalyzed by Ni was studied in decalin and water. Observations of a significant kinetic isotope effect (k_H/k_D=5.7) for the reactions of diphenyl ether under H₂ and D₂ atmosphere and a positive dependence of the rate on H₂ chemical potential in decalin indicate that addition of H to the aromatic ring is involved in the rate-limiting step. All kinetic evidence points to the fact that H addition occurs concerted with C−O bond scission. DFT calculations also suggest a route consistent with these observations involving hydrogen atom addition to the ipso position of the phenyl ring concerted with C−O scission. Hydrogenolysis initiated by H addition in water is more selective (ca. 75 %) than reactions in decalin (ca. 30 %).     

The Critical Role of Reductive Steps in the Nickel‐Catalyzed Hydrogenolysis and Hydrolysis of Aryl Ether C−O Bonds

The hydrogenolysis of the aromatic C?O bond in aryl ethers catalyzed by Ni was studied in decalin and water. Observations of a significant kinetic isotope effect (k_H/k_D=5.7) for the reactions of diphenyl ether under H₂ and D₂ atmosphere and a positive dependence of the rate on H₂ chemical potential in decalin indicate that addition of H to the aromatic ring is involved in the rate‐limiting step. All kinetic evidence points to the fact that H addition occurs concerted with C?O bond scission. DFT calculations also suggest a route consistent with these observations involving hydrogen atom addition to the ipso position of the phenyl ring concerted with C?O scission. Hydrogenolysis initiated by H addition in water is more selective (ca. 75 %) than reactions in decalin (ca. 30 %).     

The synthesis,molecular, crystal,and electronic structures of [ReBr2(O)(OCH3)(PPh3)2]

[ReBr₂(O)(OCH₃)(PPh₃)₂] has been obtained in the reaction of [ReBr₃O(PPh₃)₂] or [ReBr₂(η²-N₂COPh-N′,O)(PPh₃)₂] with an excess of methanol. [ReBr₂O(OMe)(PPh₃)₂] crystallizes in the triclinic space group P-1. The complex was characterized by infrared, UV-Vis, and ¹H NMR spectra. The electronic structure of the obtained compound has been calculated using the DFT/TD–DFT method.     

Ligand‐Exchange Processes on Solvated Zinc Cations: Water Exchange on [Zn(H2O)4(L)]2+⋅2 H2O (L=Heterocyclic Ligand)

The water‐exchange mechanisms of [Zn(H₂O)₄(L)]²⁺?2 H₂O (L=imidazole, pyrazole, 1,2,4‐triazole, pyridine, 4‐cyanopyridine, 4‐aminopyridine, 2‐azaphosphole, 2‐azafuran, 2‐azathiophene, and 2‐azaselenophene) have been investigated by DFT calculations (RB3LYP/6‐311+G**). The results support limiting associative reaction pathways that involve the formation of six‐coordinate intermediates [Zn(H₂O)₅(L)]²⁺?H₂O. The basicity of the coordinated heterocyclic ligands shows a good correlation with the activation barriers, structural parameters, and stability of the transition and intermediate states.     

Synthesis, Raman, X-ray diffraction, and density functional studies of antimony(III) heterotetracycles displaying intramolecular transannular interactions O → Sb

A set of compounds of general formula [{S(C₆H₃S)₂O}SbHal] [Hal = Cl (1), Br (2), I (3)] has been synthesized and studied by Raman and NMR spectroscopy as well as quantum chemical DFT calculations. X-ray diffraction studies of compound 2 revealed that the oxygen atom participates as donor and the antimony atom plays the role of acceptor, adopting a Ψ-distorted trigonal bi-pyramidal geometry, where the eight-membered central ring displays a boat–boat conformation. Furthermore, a series of DFT calculations was performed on compounds 1–3 as well as calculations on the non-synthesized heterotetracyclic systems [{S(C₆H₃S)₂O}SbF] (4a) and the cation [{S(C₆H₃S)₂O}Sb]⁺ (5a). The theoretic study at DFT level indicates as the electronegativity increases at exocyclic substituent along the set of compounds, the interaction is stronger. Moreover, the topological analysis of electronic density showed the existence of critical points along the O → Sb direction which prompted us to suggest the existence of a dative bond.     

Ca2+‐Induced Oxygen Generation by FeO42− at pH 9–10

Although FeO₄^2? (ferrate(IV)) is a very strong oxidant that readily oxidizes water in acidic medium, at pH 9–10 it is relatively stable (<2 % decomposition after 1 h at 298 K). However, FeO₄^2? is readily activated by Ca²⁺ at pH 9–10 to generate O₂. The reaction has the following rate law: d[O₂]/dt=k_Ca[Ca²⁺][FeO₄^2?]². ¹⁸O‐labeling experiments show that both O atoms in O₂ come from FeO₄^2?. These results together with DFT calculations suggest that the function of Ca²⁺ is to facilitate O–O coupling between two FeO₄^2‐ions by bridging them together. Similar activating effects are also observed with Mg²⁺ and Sr²⁺.     

DFT Simulation of Electron Spectra for Auger Electron and Photoelectron Spectra of Lithium Compounds

(Li, O, F)-Auger electron, and X-ray photoelectron spectra (AES, VXPS) of solid lithium compounds (Li metal, LiCl, LiF, Li₂O) are simulated by deMon density functional theory (DFT) calculations using the model molecules of the unit cell. Calculated valence XPS, core-electron binding energies (CEBE)s, and Li-, O-, and F-KVV AES for the substances correspond considerably well to experimental results. For the calculation of VXPS, the observed spectra of Li₂O pellet with chemisorbed CO₂ almost show agreement with simulation curve of the valence XPS according to the model for the 1/1 ratio of Li₂O/Li₂CO₃. In the case of AES calculation, we analyze the experimental AES with our modified Auger electron kinetic energy calculation method which corresponds to the two final-state holes at the ground state and at the transition-state in DFT calculation by removing 1 and 2 electrons, respectively. Experimental KVV AES of the Li atom, and (O, F) KVV AES of (Li₂O and LiF) in the substances almost agree well to the AES calculated with maximum kinetic energies at the ground state, and at the transition-state, respectively.