DFT Studies on Thermal Stabilities,Electronic Structures,and ~（13）C Chemical Shifts of C_（24）O_2 Based on Fullerene C_（24）（D_6）

Quantum chemical calculations on some possible equilibrium geometries of C24O2 isomers derived from C24 (D6) and C24O have been performed using density functional theory (DFT) method. The geometric and electronic structures as well as the relative energies and thermal stabilities of various C24O2 isomers at the ground state have been calculated at the B3LYP/6-31G(d) level of theory. And the 1,4,2,5-C24O2 isomer was found to be the most stable geometry where two oxygen atoms were added to the longest carbon-carbon bonds in the same pentagon from a thermodynamic point of view. Based on the optimized neutral geometries, the vertical ionization potential and vertical electron affinity have been obtained. Meanwhile, the vibrational frequencies, IR spectrum, and 13C chemical shifts of various C24O2 isomers have been calculated and analyzed.     

(BN)n富勒烯和单层BN纳米管的图形理论方法研究

应用图形理论方法对(BN)12等富勒烯和单层BN纳米管的能级分布及其稳定性进行研究,发现(BN)n比同构型的(C2)n稳定,且与用DFT方法计算的结果一致.计算结果表明,采用图形理论方法是一种很有意义的解释和预测BN纳米材料的结构和性质的定性研究方法.     

五员碳环、氮环和磷环比较的量子化学研究

用密度泛函(DFT)理论中的B3LYP和从头算(abinitio)理论中UHF,在6-31G基组水平上,对五员碳环、氮环和磷环进行几何优化计算,由所得结果讨论了分子的成键情况,对两类环(一类是C5H5、N5、P5;另一类是C5H5-、N5-和P5-)的相对稳定性分别作比较。结果表明:有机五员碳环比无机五员氮环、磷环稳定,而五员氮环与五员磷环相比,磷环更稳定。     

Theoretical investigation of the tautomerism of isoorotic acid in gaseous and aqueous phases

In the present investigation, the tautomeric and conformational equilibrium of isoorotic acid have been studied using Møller–Plesset second‐order (MP2) and density functional theory (DFT) methods in the gas phase and aqueous solution (ε = 78.5) using the IPCM model. The relative energies of these tautomers have been calculated at the two levels of theory using 6‐311++G** basis set. Energetics and relative stabilities of the tautomers were compared and analyzed in both the gaseous and aqueous phases. The results indicate that the diketo tautomer (iso) is the most stable form in the gas phase and water. The carboxylic substitution in the uracil ring does not alter its relative stability order of the tautomers. The proton affinity of the oxygen atoms and the deprotonation enthalpy of the NH bonds of isoorotic acid have been compared with recent data of uracil. The relative stability of both syn‐ and anti‐conformations was investigated and the syn form was found to be more stable by 17.65 kcal/mol. It was determined in ab initio calculations that an electron can attach to isoorotic acid, forming a stable anion better than uracil. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Evaluation of <S2> in density functional theory

The evaluation of in density functional theory (DFT) is considered. Wang et al. [J. Chem. Phys. 102, 3477 (1995)] have derived an approximate, local density expression for and in the present study their formula is evaluated using densities from unrestricted Hartree-Fock (UHF) and a range of DFT exchange-correlation functionals. The results are compared with those obtained by evaluating the conventional UHF expression using the Kohn-Sham orbitals, which is appropriate for the noninteracting system. A generalized gradient approximation for is then proposed and investigated.     

Structure and Stability of TaON Polymorphs

The stabilities and electronic structures of several polymorphs of tantalum oxynitride, TaON, were studied quantum‐chemically at density‐functional level. Results obtained by complementary quantum‐chemical techniques with wave‐functions either expanded in atom‐centered functions or in plane waves were compared. Close agreement was obtained for the relative stabilities of the baddeleyite, anatase, rutile, and fluorite phases of TaON. The effect of anion distribution on the structural parameters and the lattice energies of the anatase and rutile polymorphs was investigated. The calculated band structure of the polymorphs is compared with available experimental information.     

Structural properties and energetics of Li2FeSiO4 polymorphs and their delithiated products from first-principles

Structural properties, thermodynamic stability and delithiation process for Li(2)FeSiO(4) polymorphs are investigated by using density functional theory (DFT) within the DFT + U framework. Three Li(2)FeSiO(4) polymorphs crystallizing in space group Pmn2(1), P2(1)/n, and Pmnb have been considered. The investigations demonstrate that the strong Si-O bonds remain almost unchanged during the lithiation-delithiation process for all the polymorphs, which contribute significantly to the structural stability. On the other hand, the differences in local environment around FeO(4) tetrahedra will be translated into varying degrees of distortion, which shows a significant influence on the structural stability and average voltages. The average voltages obtained here are in good agreement with the experimental values. Furthermore, the possibility of extracting more than one lithium ions per formula unit from Li(2)FeSiO(4) of P2(1)/n is also discussed.     

Relative energy of the high-(5T2g) and low-(1A1g) spin states of the ferrous complexes [Fe(L)(NHS4)]: CASPT2 versus density functional theory

High-level ab initio calculations using multiconfigurational perturbation theory [complete active space with second-order perturbation theory (CASPT2)] were performed on the transition energy between the lowest high-spin (corresponding to (5T2g) in Oh) and low-spin (corresponding to 1A1g in Oh) states in the series of six-coordinated Fe(II) molecules [Fe(L)(NHS4)], where NHS4 is 2,2'-bis(2-mercaptophenylthio)diethylamine dianion and L=NH3, N2H4, PMe3, CO, and NO+. The results are compared to (previous and presently obtained) results from density functional theory (DFT) calculations with four functionals, which were already shown previously by Casida and co-workers [Fouqueau et al., J. Chem. Phys. 120, 9473 (2004); Ganzenmuller et al., ibid. 122, 234321 (2005); Fouqueau et al., ibid. 122, 044110 (2005); Lawson Daku et al., ChemPhysChem 6, 1393 (2005)] to perform well for the spin-pairing problem in these and other Fe(II) complexes, i.e., OLYP, PBE0, B3LYP, and B3LYP*. Very extended basis sets were used both for the DFT and CASPT2 calculations and were shown to be necessary to obtain quantitative results with both types of method. This work presents a sequel to a previous DFT/CASPT2 study of the same property in the complexes [Fe(H2O)6]2+, [Fe(NH3)6]2+, and [Fe(bpy)3]2+ [Pierloot et al., J. Chem. Phys. 125, 124303 (2006)]. The latter work was extended with new results obtained with larger basis sets and including the OLYP functional. For all considered complexes, the CASPT2 method predicts the correct ground state spin multiplicity. Since experimental data for the actual quintet-singlet (free) energy differences are not available, the performance of the different DFT functionals was judged based on the comparison between the DFT and CASPT2 results. From this, it was concluded that the generalized gradient OLYP functional performs remarkably well for the present series of ferrous compounds, whereas the success of the three hybrid functionals varies from case to case.     

Theoretical study of I2O

We present high-level density functional calculations (DFT) on the unknown I₂O molecule. The results suggest that the compound may be sufficiently stable for detection and synthesis. Our results also suggest that the DFT method is a reliable and computationally cheap alternative to G2, for estimating thermodynamic properties. The trends in relative stabilities within the HOX and X₂O series are discussed (X=halogen). © 1998 John Wiley & Sons, Inc. Heteroatom Chem 9:383–385, 1998     

Acceleration of catalyst discovery with easy,fast, and reproducible computational alchemy

The expense of quantum chemistry calculations significantly hinders the search for novel catalysts. Here, we provide a tutorial for using an easy and highly cost‐efficient calculation scheme, called alchemical perturbation density functional theory (APDFT), for rapid predictions of binding energies of reaction intermediates and reaction barrier heights based on the Kohn‐Sham density functional theory (DFT) reference data. We outline standard procedures used in computational catalysis applications, explain how computational alchemy calculations can be carried out for those applications, and then present benchmarking studies of binding energy and barrier height predictions. Using a single OH binding energy on the Pt(111) surface as a reference case, we use computational alchemy to predict binding energies of 32 variations of this system with a mean unsigned error of less than 0.05 eV relative to single‐point DFT calculations. Using a single nudged elastic band calculation for CH₄ dehydrogenation on Pt(111) as a reference case, we generate 32 new pathways with barrier heights having mean unsigned errors of less than 0.3 eV relative to single‐point DFT calculations. Notably, this easy APDFT scheme brings no appreciable computational cost once reference calculations are performed, and this shows that simple applications of computational alchemy can significantly impact DFT‐driven explorations for catalysts. To accelerate computational catalysis discovery and ensure computational reproducibility, we also include Python modules that allow users to perform their own computational alchemy calculations.     

Computational study on energetic properties of nitro derivatives of furan substituted azoles

Density functional theory has been used to investigate geometries, heats of formation (HOFs), C-NO₂ bond dissociation energies (BDEs), and relative energetic properties of nitro derivatives of azole substituted furan. HOFs for a series of molecules were calculated by using density functional theory (DFT) and Møller–Plesset (MP2) methods. The density is predicted using crystal packing calculations; all the designed compounds show density above 1.71 g/cm³. The calculated detonation velocities and detonation pressures indicate that the nitro group is very helpful for enhancing the detonation performance for the designed compounds. Thermal stabilities have been evaluated from the bond dissociation energies. Charge on the nitro group was used to assess the impact sensitivity in this study. According to the results of the calculations, tri- and tetra-nitro substituted derivatives reveal high performance with better thermal stability.     

Structural and spectral properties of 3-substitutedphenyl-1,5-diphenylformazans: A quantum chemical study

The structural and optical properties of 3-substitutedphenyl-1,5-diphenylformazans are studied by quantum chemical methods. The density functional theory (DFT) is employed to optimize the ground state geometries of formazans substituted with different electron donating and withdrawing groups in both gas and solvent phases. The absorption spectra of formazan derivatives are calculated using time dependent density functional theory (TD-DFT). The polarizable continuum model (PCM) calculations of 3-substitutedphenyl-1,5-diphenylformazans are performed for bulk solvent effects. The geometrical parameters, vibrational frequencies, and relative stabilities of isomers of 3-substitutedphenyl-1,5-diphenylformazans are studied. The results obtained by TD-DFT calculations reveal that the substitution of electron withdrawing and donating substituents affects the absorption spectra of 3-substitutedphenyl-1,5-diphenylformazans. The calculated maximum absorption wavelengths (λ_max) are highly consistent with the experimental values as found from UV-vis spectra.     

Ab initio base-pairing energies of an oxidized thymine product, 5-formyluracil, with standard DNA bases at the BSSE-free DFT and MP2 theory levels

Oxidation of the thymine methyl group produces two stable products, non-mutagenic 5-hydroxymethyluracil and highly mutagenic 5-formyluracil. We have calculated the interaction energy of base-pair formation involving 5-formyluracil bound to the natural DNA bases adenine (A), cytosine (C), guanine (G), and thymine (T), and discuss the effects of the 5-formyl group with respect to similar base-pairs containing uracil, 5-hydroxyuracil, thymine (5-methyluracil), and 5-hydroxycytosine. The interaction geometries and energies were calculated four ways: (a) using density functional theory (DFT) without basis set super-position error (BSSE) corrections, (b) using DFT with BSSE correction of geometries and energies, (c) using M?ller-Plesset second order perturbation theory (MP2) without BSSE correction, and (d) using MP2 with BSSE geometry and energy correction. All calculations used the 6-311G(d,p) basis set. Notably, we find that the A:5-formyluracil base-pair is more stable than the precursor A:T base-pair. The relative order of base-pair stabilities is A:5-Fo-U > G:5-Fo-U > C:5-Fo-U > T:5-Fo-U.     

Derivation of force field parameters for SnO2-H2O surface systems from plane-wave density functional theory calculations

Plane-wave density functional theory (DFT-PW) calculations were performed on bulk SnO2 (cassiterite) and the (100), (110), (001), and (101) surfaces with and without H2O present. A classical interatomic force field has been developed to describe bulk SnO2 and SnO2-H2O surface interactions. Periodic density functional theory calculations using the program VASP (Kresse et al., 1996) and molecular cluster calculations using Gaussian 03 (Frisch et al., 2003) were used to derive the parametrization of the force field. The program GULP (Gale, 1997) was used to optimize parameters to reproduce experimental and ab initio results. The experimental crystal structure and elastic constants of SnO2 are reproduced reasonably well with the force field. Furthermore, surface atom relaxations and structures of adsorbed H2O molecules agree well between the ab initio and force field predictions. H2O addition above that required to form a monolayer results in consistent structures between the DFT-PW and classical force field results as well.     

Two- and three-body interactions among nanoparticles in a polymer melt

We perform direct three-dimensional density functional theory (DFT) calculations of two- and three-body interactions in polymer nanocomposites. The nanoparticles are modeled as hard spheres, immersed in a hard-sphere homopolymer melt of freely jointed chains. The two-particle potential of mean force obtained from the DFT is in near quantitative agreement with the potential of mean force obtained from self-consistent polymer reference interaction site model theory. Three-body interactions among three nanoparticles are found to be significant, such that it is not possible to describe these systems with a polymer-mediated two-body interaction calculated from the potential of mean force.