Improved exchange-correlation potential for polarizability and dissociation in density functional theory

The authors propose a novel approach to the problem of polarizabilities and dissociation in electric fields from the static limit of the Vignale-Kohn (VK) functional. The response to the purely scalar part of the VK response potential is considered. This potential has ground-state properties that notably improve over the full VK response density and over usual (semi-)local functionals. The correct qualitative behavior of our potentials means that it is expected to work well for polarizabilities in cases such as the H(2) chain, and it will also correctly dissociate open-shell fragments in a field.     

Role of the exchange-correlation potential in ab initio electron transport calculations

The effect of the exchange-correlation potential in ab initio electron transport calculations is investigated by constructing optimized effective potentials using different energy functionals or the electron density from second-order perturbation theory. The authors calculate electron transmission through two atomic chain systems, one with charge transfer and one without. Dramatic effects are caused by two factors: changes in the energy gap and the self-interaction error. The error in conductance caused by the former is about one order of magnitude while that caused by the latter ranges from several times to two orders of magnitude, depending on the coupling strength and charge transfer. The implications for accurate quantum transport calculations are discussed.     

An ab initio and density functional theory study of radical-clock reactions

Density functional theory (DFT) and ab initio (CBS-RAD) calculations have been used to investigate a series of "radical clock" reactions. The calculated activation energies suggest that the barriers for these radical rearrangements are determined almost exclusively by the enthalpy effect with no evidence of significant polar effects. The ring-closure reactions to cyclopentylmethyl radical derivatives and the ring opening of cyclopropylmethyl radicals give different correlations between the calculated heat of reaction and barrier, but the two types of reaction are internally consistent.     

Vibrational spectroscopy investigation using ab initio and density functional theory on p-anisaldehyde

The FTIR and FT Raman spectra of p-anisaldehyde has been recorded in the regions 4,000-400 and 3,500-100 cm(-1), respectively. The optimized geometry, frequency and intensity of the vibrational bands of p-anisaldehyde were obtained by ab initio and DFT levels of theory with complete relaxation in the potential energy surface using 6-31G(d,p) basis set. A complete vibrational assignment aided by the theoretical harmonic frequency analysis has been proposed. The harmonic vibrational frequencies calculated have been compared with experimental FTIR and FT Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically constructed bar type spectrograms.     

Time-dependent density functional theory based on a noncollinear formulation of the exchange-correlation potential

In this study we have introduced a formulation of time-dependent density functional theory (TDDFT) based on a noncollinear exchange-correlation potential. This formulation is a generalization of conventional TDDFT. The form of this formulation is exactly the same as that of the conventional TDDFT for the excitation energies of transitions that do not involve spin flips. In addition, this noncollinear TDDFT formulation allows for spin-flip transitions. This feature makes it possible to resolve more fully excited state spin multiplets, while for closed-shell systems, the spin-flip transitions will result in singlet-triplet excitations and this excitation energy calculated from this formulation of TDDFT is exactly the same as that from ordinary TDDFT. This formulation is applied to the dissociation of H(2) in its (1)Sigma(g) (+) ground state and (1)Sigma(u) (+) and (3)Sigma(u) (-) excited states with (3)Sigma(u) (-) (M(s)=+1) as the reference state and the multiplets splitting of some atoms.     

Electron affinities of metals computed by density functional theory and ab initio methods

An extensive computational study of the meal electron affinity was performed using the ab initio and density functional theory (DFT) methods. HF, MP2, MP3, MP4, QCISD, and QCISD(T) was used as computational methods, while the hybrid, local, and nonlocal DFT methods with the LYP, P86, PW91, and VWN correlation functionals were used. Two basis sets, one small and applicable to almost all metals (LanL2DZ) and one large [6-311 + + G(3df, 3 pd)] used only for small metals, were employed. The computed results were compared with the experimental data and the capabilities of the DFT methods to perform this study were discussed. © 1997 John Wiley & Sons, Inc.     

Carbohydrate-protein recognition probed by density functional theory and ab initio calculations including dispersive interactions

Carbohydrate-protein recognition has been studied by electronic structure calculations of complexes of fucose and glucose with toluene, p-hydroxytoluene and 3-methylindole, the latter aromatic molecules being analogues of phenylalanine, tyrosine and tryptophan, respectively. We use mainly a density functional theory model with empirical corrections for the dispersion interactions (DFT-D), this method being validated by comparison with a limited number of high level ab initio calculations. We have calculated both binding energies of the complexes as well as their harmonic vibrational frequencies and proton NMR chemical shifts. We find a range of minimum energy structures in which the aromatic group can bind to either of the two faces of the carbohydrate, the binding being dominated by a combination of OH-pi and CH-pi dispersive interactions. For the fucose-toluene and alpha-methyl glucose-toluene complexes, the most stable structures involve OH-pi interactions, which are reflected in a red shift of the corresponding O-H stretching frequency, in good quantitative agreement with experimental data. For those structures where CH-pi interactions are found we predict a corresponding blue shift in the C-H frequency, which parallels the predicted proton NMR shift. We find that the interactions involving 3-methylindole are somewhat greater than those for toluene and p-hydroxytoluene.     

Metahybrid density functional theory and correlated ab initio studies on microhydrated adenine-thymine base pairs

Monohydrated and dihydrated adenine-thymine base pairs are characterized using metahybrid density functional theory and correlated ab initio approaches. The motivation of this work is twofold. First, the high-level geometries and interaction energies computed for different complexes serve as a reference for the testing of recently developed density functional theory (DFT) with respect to its ability to correctly describe the balance between the electrostatic and the dispersion contributions that bind these complexes. Second, these studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for a better understanding of the influence of the solvent on the stability of the structure of DNA duplexes.     

Differences between ab initio and density functional electron densities

We analyzed the energy contributions and the spatial distribution differences of several electron densities of atoms and small molecules. The results show the insensitivity of local spin density correlation functionals in respect to differences in the electron densities. On the other hand, significant changes in one-electron and two-electron energy contributions are observed, although both compensate each other. The projection of the differences between these electron densities, referred to as the Hartree-Fock density, shows a qualitative resemblance between multideterminantal and Kohn-Sham wavefunctions. Finally, a comparative analysis of the optimized conformational parameters obtained using several methods shows that the inclusion of the correlation energy in SCF or in post-SCF procedures gives similar results and that the exchange potential is more important than is the correlation potential to improve these conformational parameters. © 1997 John Wiley & Sons, Inc.     

The structure of 1-chlorosilatrane: An ab initio molecular orbital and a density functional theory study

The molecular geometries of the 1-chloro-, 1-fluoro-, 1-methyl-, and 1-hydrogenosilatranes were fully optimized by the restricted Hartree-Fock (HF) method supplemented with 3-21G, 3-21G(d), 6-31G(d), and CEP-31G(d) basis sets; by MP2 calculations using 6-31G(d) and CEP-31G(d) basis sets; and by GGA-DFT calculations using 6-31G(d5) basis set with the aim of locating the positions of the local minima on the energy hypersurface. The HF/6-31G(d) calculations predict long (>254 pm) and the MP2/CEP calculations predicted short (∼225 pm) equilibrium Si(SINGLE BOND)N distances. The present GGA-DFT calculations reproduce the available gas phase experimental Si(SINGLE BOND)N distances correctly. The solid phase experimental results predict that the Si(SINGLE BOND)N distance is shorter in 1-chlorosilatrane than in 1-fluorosilatrane. In this respect the HF results show a strong basis set dependence, the MP2/CEP results contradict the experiment, and the GGA-DFT results in electrolytic medium agree with the experiment. The latter calculations predict that 1-chlorosilatrane is more polarizable than 1-fluorosilatrane and also support a general Si(SINGLE BOND)N distance shortening trend for silatranes during the transition from gas phase to polar liquid or solid phase. The calculations predict that the ethoxy links of the silatrane skeleton are flexible. Consequently, it is difficult to measure experimentally the related bond lengths and bond and torsion angles. This is the probable origin of the surprisingly large differences for the experimental structural parameters. On the basis of experimental analogies, ab initio calculations, and density functional theory (DFT) calculations, a gas phase equilibrium (r_e) geometry is predicted for 1-chlorosilatrane. The semiempirical methods predict a so-called exo minimum (at above 310 pm Si(SINGLE BOND)N distance); however, the ab initio and GGA-DFT calculations suggest that this form is nonexistent. The GGA-DFT geometry optima were characterized by frequency analysis. © 1996 by John Wiley & Sons, Inc.     

In cation interactions with some organics: ab initio molecular orbital and density functional theory

Ab initio molecular orbital and density functional theory studies were undertaken to investigate the structural and energetic characteristics of complexes of In⁺ with several different organic molecules for the first time. HF, MP2, QCISD, and CCD levels of theory in ab initio MO as well as B3LYP, B3PW91 hybrid functionals in density functional theory were used. A valence TZ+P basis set with relativistic effective core potentials was used for the In atom while the 6-311++G(3d, 2p) basis set was utilized for all other atoms. Both closed-shell (H₂O, CH₄, CH₃OH, and C₆H₆) and open-shell (CH₃ and C₂H₃) molecules were considered for complexation with In⁺. In⁺ affinities of 21.5, 24.8, 28.6, 18.4, and 23.0 kcal/mol were obtained with the B3PW91 hybrid functional for H₂O, CH₃OH, C₆H₆, CH₃, and C₂H₃, respectively. The large values for the calculated affinities indicate the validity of our recent experimental detection of In⁺ ion attachment to some organic molecules.     

Halide adsorption on single-crystal silver substrates: dynamic simulations and ab initio density functional theory

We investigate the static and dynamic behaviors of a Br adlayer electrochemically deposited onto single-crystal Ag(100) using an off-lattice model of the adlayer. Unlike previous studies using a lattice-gas model, the off-lattice model allows adparticles to be located at any position within a two-dimensional approximation to the substrate. Interactions with the substrate are approximated by a corrugation potential. Using density functional theory (DFT) to calculate surface binding energies, a sinusoidal approximation to the corrugation potential is constructed. A variety of techniques, including Monte Carlo and Langevin simulations, are used to study the behavior of the adlayer. The lateral root-mean-square (rms) deviation of the adparticles from the binding sites is presented along with equilibrium coverage isotherms, and the thermally activated Arrhenius barrier-hopping model used in previous dynamic Monte Carlo simulations is tested.     

The performance of time-dependent density functional theory based on a noncollinear exchange-correlation potential in the calculations of excitation energies

In the present work we have studied the accuracy of excitation energies calculated from spin-flip transitions with a formulation of time-dependent density functional theory based on a noncollinear exchange-correlation potential proposed in a previous study. We compared the doublet-doublet excitation energies from spin-flip transitions and ordinary transitions, calculated the multiplets splitting of some atoms, the singlet-triplet gaps of some diradicals, the energies of excited quartet states with a doublet ground state. In addition, we attempted to calculate transition energies with excited states as reference. We compared the triplet excitation energies and singlet-triplet separations of the excited state from spin-flip and ordinary transitions. As an application, we show that using excited quartet state as reference can help us fully resolve excited states spin multiplets. In total the obtained excitation energies calculated from spin-flip transitions agree quite well with other theoretical results or experimental data.     

Molecular structure and vibrational spectra of phenobaraitone by density functional theory and ab initio hartree-Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra at RHF/6-31++G** and B3LYP/6-31++G** levels for phenobarbitone (C₁₂H₁₂N₂O₃) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. Theoretical vibrational spectra of the title compound were interpreted by means of potential energy distributions (PEDs) using MOLVIB program. A detailed interpretation of the infrared spectra of the title compound is reported. On the basis of the agreement between the calculated and observed results, the assignments of fundamental vibrational modes of phenobarbitone were examined and some assignments were proposed. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title compound have been constructed.     

An ab initio and density functional theory study of the structure and bonding of sulfur ylides

Sulfur ylides are useful synthetic intermediates that are formed from the interaction between singlet carbenes and sulfur-containing molecules. Partial double-bond character frequently has been proposed as a key contributor to the stability of sulfur ylides. Calculations at the B3LYP, MP2, and CCSD(T) levels of theory employing various basis sets have been performed on the sulfur ylides H(2)S-CH(2) and (CH(3))(2)S-CH(2) in order to investigate the structure and bonding of these systems. The following general properties of sulfur ylides were observed from the computational studies: C-S bond distances that are close in length to that of a typical C-S double bond, high charge transfer from the sulfide to the carbene, and large torsional rotation barriers. Analysis of the sulfur ylide charge distribution indicates that the unusually short C-S bond distance can be attributed in part to the electrostatic attraction between highly oppositely charged carbon and sulfur atoms. Furthermore, n --> sigma* stabilization arising from donation of electron density from the carbon lone pair orbital into S-H or S-C antibonding orbitals leads to larger than expected torsional barriers. Finally, natural resonance theory analysis indicates that the bond order of the sulfur ylides H(2)S-CH(2) and (CH(3))(2)S-CH(2) is 1.4-1.5, intermediate between a single and double bond.     

The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations

A new algorithm for density-functional-theory-based ab initio molecular dynamics simulations is presented. The Kohn–Sham orbitals are expanded in Gaussian-type functions and an augmented-plane-wave-type approach is used to represent the electronic density. This extends previous work of ours where the density was expanded only in plane waves. We describe the total density in a smooth extended part which we represent in plane waves as in our previous work and parts localised close to the nuclei which are expanded in Gaussians. Using this representation of the charge we show how the localised and extended part can be treated separately, achieving a computational cost for the calculation of the Kohn–Sham matrix that scales with the system size N as O(NlogN). Furthermore, we are able to reduce drastically the size of the plane-wave basis. In addition, we introduce a multiple-cutoff method that improves considerably the performance of this approach. Finally, we demonstrate with a series of numerical examples the accuracy and efficiency of the new algorithm, both for electronic structure calculations and for ab initio molecular dynamics simulations. Received: 15 December 1998 /Accepted: 18 February 1999 /Published online: 14 July 1999     

Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations

Ab initio MP2/6-31G* interaction energies were calculated for more than 80 geometries of stacked cytosine dimer. Diffuse polarization functions were used to properly cover the dispersion energy. The results of ab initio calculations were compared with those obtained from three electrostatic empirical potential models, constructed as the sum of a Lennard-Jones potential (covering dispersion and repulsion contributions) and the electrostatic term. Point charges and point multipoles of the electrostatic term were also obtained at the MP2/6-31G* level of theory. The point charge MEP model (atomic charges derived from molecular electrostatic potential) satisfactorily reproduced the ab initio data. Addition of π-charges localized below and above the cytosine plane did not affect the calculated energies. The model employing the distributed multipole analysis gave worse agreement with the ab initio data than the MEP approach. The MP2 MEP charges were also derived using larger sets of atomic orbitals: cc-pVDZ, 6-311 + G(2d, p), and aug-cc-pVDZ. Differences between interaction energies calculated using these three sets of point charges and the MP2/6-31G* charges were smaller than 0.8 kcal/mol. The correlated ab initio calculations were also compared with the density functional theory (DFT) method. DFT calculations well reproduced the electrostatic part of interaction energy. They also covered some nonelectrostatic short-range effects which were not reproduced by the empirical potentials. The DFT method does not include the dispersion energy. This energy, approximated by an empirical term, was therefore added to the DFT interaction energy. The resulting interaction energy exhibited an artifact secondary minimum for a 3.9-4.0 vertical separation of bases. This defect is inherent in the DFT functionals, because it is not observed for the Hartree-Fock + dispersion interaction energy.© 1996 John Wiley & Sons, Inc.     

The C-N rotation barrier of the lithium enolate of acetamide: an ab initio and density functional theory investigation

Semiempirical (PM3), ab initio (HF/6-31+G(d) and MP2/6-31+G(d)), and density functional (pBP/DN) calculations are used to investigate the rotation barrier of the carbon-nitrogen bond in a simple enolate anion: lithium acetamide, 1. For comparison, the amidate anion 2, vinylamine 3, and a simulated dimer 4 were also calculated. In all systems, the barrier to rotation was found to be less than 10 kcal x mol(-1) in agreement with experiment. The correlated calculations show the barrier to be lowest for the anion 2. The results show conjugation effects in 1 and 2 comparable to that in vinylamine 3 and imply that polarization effects are more important than charge transfer in amine conjugation.     

Computation of bond dissociation energy for sulfides and disulfides with ab initio and density functional theory methods

The geometries and S-H, S-S, and S-C bond dissociation energies for hydrogen sulfide, hydrogen disulfide, methanethiol, dimethyl disulfide, and dimethyl disulfide were calculated with both ab initio (ROHF and MP2), hybrid (BHandH, BHandHLYP, Becke3LYP and Becke3P86), and nonlocal (BLYP and BP86) density functional theory (DFT) methods. In all studies the 6–31 + G(d) basis set is used. The computed results are compared to the experimentally obtained values, targeting the selection of a suitable ab initio or DFT method for the study of these systems. © 1997 John Wiley & Sons, Inc.