The local orbital energy and density functional theory

From the viewpoint of density functional theory, an expression is derived which improves the average energy of a trial density. Applications to atoms and molecules are made using wave function methods and are based on properties of the variance, which is defined as $ (\overline {\varepsilon ^2 } - (\overline \varepsilon)^2)^{1/2} $, where ? is the local orbital energy. Calculated results for both Hartree-Fock and correlated wave functions are quite encouraging.     

The magnetic and electronic structure of vanadyl pyrophosphate from density functional theory

We have studied the magnetic structure of the high symmetry vanadyl pyrophosphate ((VO)(2)P(2)O(7), VOPO), focusing on the spin exchange couplings, using density functional theory (B3LYP) with the full three-dimensional periodicity. VOPO involves four distinct spin couplings: two larger couplings exist along the chain direction (a-axis), which we predict to be antiferromagnetic, J(OPO) = -156.8 K and J(O) = -68.6 K, and two weaker couplings appear along the c (between two layers) and b directions (between two chains in the same layer), which we calculate to be ferromagnetic, J(layer) = 19.2 K and J(chain) = 2.8 K. Based on the local density of states and the response of spin couplings to varying the cell parameter a, we found that J(OPO) originates from a super-exchange interaction through the bridging -O-P-O- unit. In contrast, J(O) results from a direct overlap of 3d(x(2)-y(2)) orbitals on two vanadium atoms in the same V(2)O(8) motif, making it very sensitive to structural fluctuations. Based on the variations in V-O bond length as a function of strain along a, we found that the V-O bonds of V-(OPO)(2)-V are covalent and rigid, whereas the bonds of V-(O)(2)-V are fragile and dative. These distinctions suggest that compression along the a-axis would have a dramatic impact on J(O), changing the magnetic structure and spin gap of VOPO. This result also suggests that assuming J(O) to be a constant over the range of 2-300 K whilst fitting couplings to the experimental magnetic susceptibility is an invalid method. Regarding its role as a catalyst, the bonding pattern suggests that O(2) can penetrate beyond the top layers of the VOPO surface, converting multiple V atoms from the +4 to +5 oxidation state, which seems crucial to explain the deep oxidation of n-butane to maleic anhydride.     

Time dependent density functional theory study of the near-edge x-ray absorption fine structure of benzene in gas phase and on metal surfaces

The near-edge x-ray absorption fine structure of benzene in the gas phase and adsorbed on the Au(111) and Pt(111) surfaces is studied with time dependent density functional theory. Excitation energies computed with hybrid exchange-correlation functionals are too low compared to experiment. However, after applying a constant shift the spectra are in good agreement with experiment. For benzene on the Au(111) surface, two bands arising from excitation to the e(2u)(pi(*)) and b(2g)(pi(*)) orbitals of benzene are observed for photon incidence parallel to the surface. On Pt(111) surface, a broader band arises from excitation to benzene orbitals that are mixed with the surface and have both sigma(*)(Pt-C) and pi(*) characters.     

The hydrogen bond in ice probed by soft x-ray spectroscopy and density functional theory

We combine photoelectron and x-ray absorption spectroscopy with density functional theory to derive a molecular orbital picture of the hydrogen bond in ice. We find that the hydrogen bond involves donation and back-donation of charge between the oxygen lone pair and the O-H antibonding orbitals on neighboring molecules. Together with internal s-p rehybridization this minimizes the repulsive charge overlap of the connecting oxygen and hydrogen atoms, which is essential for a strong attractive electrostatic interaction. Our joint experimental and theoretical results demonstrate that an electrostatic model based on only charge induction from the surrounding medium fails to properly describe the internal charge redistributions upon hydrogen bonding.     

C-C bond unsaturation degree in monosubstituted ferrocenes for molecular electronics investigated by a combined near-edge x-ray absorption fine structure, x-ray photoemission spectroscopy, and density functional theory approach

We present the results of an experimental and theoretical investigation of monosubstituted ethyl-, vinyl-, and ethynyl-ferrocene (EtFC, VFC, and EFC) free molecules, obtained by means of synchrotron-radiation based C 1s photoabsorption (NEXAFS) and photoemission (C 1s XPS) spectroscopies, and density functional theory (DFT) calculations. Such a combined study is aimed at elucidating the role played by the C-C bond unsaturation degree of the substituent on the electronic structure of the ferrocene derivatives. Such substituents are required for molecular chemical anchoring onto relevant surfaces when ferrocenes are used for molecular electronics hybrid devices. The high resolution C 1s NEXAFS spectra exhibit distinctive features that depend on the degree of unsaturation of the hydrocarbon substituent. The theoretical approach to consider the NEXAFS spectrum made of three parts allowed to disentangle the specific contribution of the substituent group to the experimental spectrum as a function of its unsaturation degree. C 1s IEs were derived from the experimental data analysis based on the DFT calculated IE values for the different carbon atoms of the substituent and cyclopentadienyl (Cp) rings. Distinctive trends of chemical shifts were observed for the substituent carbon atoms and the substituted atom of the Cp ring along the series of ferrocenes. The calculated IE pattern was rationalized in terms of initial and final state effects influencing the IE value, with special regard to the different mechanism of electron conjugation between the Cp ring and the substituent, namely the σ/π hyperconjugation in EtFC and the π-conjugation in VFC and EFC.     

Vibronic structure of the permanganate absorption spectrum from time-dependent density functional calculations

The UV absorption spectrum of the permanganate anion is a prototype transition-metal complex spectrum. Despite this being a simple d0 Td system, for which a beautiful spectrum with detailed vibrational structure has been available since 1967, the assignment of the second and third bands is still very controversial. The issue can be resolved only by an elucidation of the intricate vibronic structure of the spectrum. We investigate the vibronic coupling by means of linear-response time-dependent density functional calculations. By means of a diabatizing scheme that employs the transition densities obtained in the TDDFT calculations in many geometries around Re, we construct a Taylor series expansion in the normal coordinates of a diabatic potential energy matrix, coupling 24 excited states. The simulated vibronic structure is in good agreement with the experimental absorption spectrum after the adjustment of some of the calculated vertical excitation energies. The peculiar blurred vibronic structure of the second band, which is a very distinctive feature of the experimental spectrum, is fully reproduced in the calculations. It is caused by the double-well shape of the adiabatic energy surface along the Jahn-Teller active e mode of the allowed 1E state arising from the second 1T2 state, which exhibits a Jahn-Teller splitting into 1B2 and 1E states. We trace the double-well shape to an avoided crossing between two diabatic states with different orbital-excitation character. The crossing can be explained at the molecular orbital level from the Jahn-Teller splitting of the set of 7t2{3d(xy), 3d(xz), 3d(yz)} orbitals (the LUMO + 1), to which the excitations characterizing the diabatic states take place. In contrast to its character in the two well regions, at Re the 2(1)T2 state is not predominantly an excitation to the LUMO + 1, but has more HOMO - 1 --> LUMO (2e = {3d(x2-y2), 3d(z2)}) character. The changing character of the 2(1)T2 - 1E state along the e mode implies that the assignment of the experimental bands to single orbital transitions is too simplistic intrinsically. This spectrum, and notably the blurring of the vibronic structure in the second band, can be understood only from the extensive configurational mixing and vibronic coupling between the excited states. This solves the long-standing assignment problem of these bands.     

A high-pressure mass spectrometric and density functional theory investigation of the thermochemical properties and structure of protonated dimers and trimers of glycine

A new modification of pulsed-ionization high-pressure mass spectrometry (PHPMS) has been used to perform equilibrium thermochemical studies for relatively nonvolatile biomolecules such as amino acids. Binding enthalpy and entropy changes have been measured for proton-bound clusters of glycine, which are in good agreement with both theoretical (DFT) results of this work and a previous blackbody infrared dissociation experiment. Experimental data indicate that a number of conformers of the proton-bound dimer of glycine may coexist in the explored temperature range (360-460 K). Several new, conceptually different isomers (two of them zwitterionic) have been found by DFT calculations, one of which is 7 kJ mol(-1) lower in energy than the structure previously reported to be the energy minimum.     

A combined density functional and x-ray diffraction study of Pt nanoparticle structure

The structure of 1.7 nm Pt nanoparticles is investigated using x-ray diffraction (XRD) measurements and density functional theory (DFT) calculations. Two types of particles are compared, those made by solution chemistry which are capped either by thiol or amine ligands, and dendrimer encapsulated particles (DENs) which do not have capping ligands. All particles were dried before analyzing their structure. Pair distribution function (PDF) data from XRD measurements show that the ligand-capped particles are more disordered than the DENs. To determine the structure of the particles and the nature of the ligand-induced disorder, we use a hybrid reverse Monte Carlo approach. A weighted average of the calculated binding energy of the particles and a goodness-of-fit parameter to the PDF data is taken as the object function, which is minimized to determine the optimal structure. A scan over different weights gives the set of pareto optimal structures, which show how well simultaneous agreement can be reached to both experiment and theory. Using an embedded atom potential to sample configuration space and DFT to refine the optimal structures, we show that the DEN structure is most consistent with a face centered cubic lattice of truncated octahedral shape. The disorder induced by the capping ligands is consistent with surface relaxation of the particle rather than disorder of the crystal structure.     

Molecular density functional theory of solvation: from polar solvents to water

A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics.     

The structure of the second-order reduced density matrix in density matrix functional theory and its construction from formal criteria

Some formal requirements for the second-order reduced density matrix are discussed in the context of density matrix functional theory. They serve as a basis for the ad hoc construction of the second-order reduced density matrix in terms of the first-order reduced density matrix and lead to implicit functionals where the occupation numbers of the natural orbitals are obtained as diagonal elements of an idempotent matrix the elements of which represent the variational parameters to be optimized. The numerical results obtained from a first realization of such an implicit density matrix functional give excellent agreement with the results of full configuration interaction calculations for four-electron systems like LiH and Be. Results for H2O taken as an example for a somewhat larger molecule are numerically less satisfactory but still give reasonable occupation numbers of the natural orbitals and indicate the capability of density matrix functional theory to cope with static electron correlation.     

Combined X-ray absorption spectroscopy and density functional theory examination of ferrocene-labeled peptides

A combination of soft X-ray absorption spectroscopy (XAS) measurements and StoBe density functional theory (DFT) calculations has been used to study the electronic structures of the ferrocene-labeled peptides Fc-Pro(n)-OBz (n = 1-4). Excellent agreement between the measured and the simulated data is observed in all cases, and the origin of all major spectral features was assigned. The breaking of the degeneracy of the ferrocene 3e(2u)-like unoccupied molecular orbital under the influence of a substituent attached to a Cp ring was observed experimentally. The influence of the bonding environment on the O 1s and N 1s XAS spectra was examined. A corrected assignment of one of the major features in the Fe 2p XAS spectra of ferrocene is proposed and supported by the DFT simulations, as well as the measured spectra.     

Thomas–Fermi theory from a perturbative treatment of the hydrogenic limit energy density functional: A possible link to local density functional theory

A perturbation expansion which connects the hydrogenic limit energy density functional to the Thomas–Fermi functional is discussed. This perturbation series, where the Coulomb energy density functional is treated as the perturbation to the hydrogenic limit functional, is, in fact, the q = (N/Z) expansion of Thomas–Fermi theory. A truncated form of the first-order correction to the functional provides further insight into the model which treats the ground state energy as a local functional of the electron density.     

Prediction of iron K-edge absorption spectra using time-dependent density functional theory

Iron K-edge X-ray absorption pre-edge features have been calculated using a time-dependent density functional approach. The influence of functional, solvation, and relativistic effects on the calculated energies and intensities has been examined by correlation of the calculated parameters to experimental data on a series of 10 iron model complexes, which span a range of high-spin and low-spin ferrous and ferric complexes in O(h) to T(d) geometries. Both quadrupole and dipole contributions to the spectra have been calculated. We find that good agreement between theory and experiment is obtained by using the BP86 functional with the CP(PPP) basis set on the Fe and TZVP one of the remaining atoms. Inclusion of solvation yields a small improvement in the calculated energies. However, the inclusion of scalar relativistic effects did not yield any improved correlation with experiment. The use of these methods to uniquely assign individual spectral transitions and to examine experimental contributions to backbonding is discussed.     

Exact exchange-only density functional theory by means of local scaling transformations

The density-constrained variation of the kinetic energy of a non-interacting system carried out within the framework of local-scaling transformations is employed in the present work for the purpose of determining exchange-only Kohri Sham orbitals for the beryllium atom as well as the Kohn-Sham exchange-only correlation potential. The starting basis functions are those of Clementi-Roetti and Raffenetti-types. We have also performed an optimization via a densitydriven method. The resulting exchange-only potentials for the examples studied are almost indistinguishable from the Talman-Shadwick potential for the beryllium atom.     

Hybrid density functional theory band structure engineering in hematite

We present a hybrid density functional theory (DFT) study of doping effects in α-Fe(2)O(3), hematite. Standard DFT underestimates the band gap by roughly 75% and incorrectly identifies hematite as a Mott-Hubbard insulator. Hybrid DFT accurately predicts the proper structural, magnetic, and electronic properties of hematite and, unlike the DFT+U method, does not contain d-electron specific empirical parameters. We find that using a screened functional that smoothly transitions from 12% exact exchange at short ranges to standard DFT at long range accurately reproduces the experimental band gap and other material properties. We then show that the antiferromagnetic symmetry in the pure α-Fe(2)O(3) crystal is broken by all dopants and that the ligand field theory correctly predicts local magnetic moments on the dopants. We characterize the resulting band gaps for hematite doped by transition metals and the p-block post-transition metals. The specific case of Pd doping is investigated in order to correlate calculated doping energies and optical properties with experimentally observed photocatalytic behavior.     

Energies of ions in water and nanopores within density functional theory

Accurate calculations of electrostatic potentials and treatment of substrate polarizability are critical for predicting the permeation of ions inside water-filled nanopores. The ab initio molecular dynamics method, based on density functional theory (DFT), accounts for the polarizability of materials, water, and solutes, and it should be the method of choice for predicting accurate electrostatic energies of ions. In practice, DFT coupled with the use of periodic boundary conditions in a charged system leads to large energy shifts. Results obtained using different DFT packages may vary because of the way pseudopotentials and long-range electrostatics are implemented. Using maximally localized Wannier functions, we apply robust corrections that yield relatively unambiguous ion energies in select molecular and aqueous systems and inside carbon nanotubes. Large binding energies are predicted for ions in metallic carbon nanotube arrays, while Na+ and Cl- energies are found to exhibit asymmetry in water that is smaller than but comparable with those computed using nonpolarizable water force fields.     

Magnetizability tensors from auxiliary density functional theory

The working equations for the calculation of the magnetizability tensor in the framework of auxiliary density functional theory with gauge including atomic orbitals (ADFT-GIAO) are derived. Unlike in the corresponding conventional density functional theory implementations the numerical integration of the GIAOs is avoided in ADFT-GIAO. Our validation shows that this simplification has no effect on the accuracy of the methodology. As a result, a reliable and efficient implementation for the calculation of magnetizabilities of systems with more than 1000 atoms and 14 000 basis functions is presented.