Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set

We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.     

Time-dependent density functional theory scheme for efficient calculations of dynamic (hyper)polarizabilities

The authors present an efficient perturbative method to obtain both static and dynamic polarizabilities and hyperpolarizabilities of complex electronic systems. This approach is based on the solution of a frequency-dependent Sternheimer equation, within the formalism of time-dependent density functional theory, and allows the calculation of the response both in resonance and out of resonance. Furthermore, the excellent scaling with the number of atoms opens the way to the investigation of response properties of very large molecular systems. To demonstrate the capabilities of this method, they implemented it in a real-space (basis-set-free) code and applied it to benchmark molecules, namely, CO, H2O, and para-nitroaniline. Their results are in agreement with experimental and previous theoretical studies and fully validate their approach.     

Lithium adsorption on graphite from density functional theory calculations

The structural, energetic, and electronic properties of the Li/graphite system are studied through density functional theory (DFT) calculations using both the local spin density approximation (LSDA), and the gradient-corrected Perdew-Burke-Ernzerhof (PBE) approximation to the exchange-correlation energy. The calculations were performed using plane waves basis, and the electron-core interactions are described using pseudopotentials. We consider a disperse phase of the adsorbate comprising one Li atom for each 16 graphite surface cells, in a slab geometry. The close contact between the Li nucleus and the graphene plane results in a relatively large binding energy (larger than 1.1 eV). A detailed analysis of the electronic charge distribution, density difference distribution, and band structures indicates that one valence electron is entirely transferred from the atom to the surface, which gives rise to a strong interaction between the resulting lithium ion and the cloud of pi electrons in the substrate. We show that it is possible to explain the differences in the binding of Li, Na, and K adatoms on graphite considering the properties of the corresponding cation/aromatic complexes.     

Embedding wave function theory in density functional theory

We present a framework for embedding a highly accurate coupled-cluster calculation within a larger density functional calculation. We use a perturbative buffer to help insulate the coupled-cluster region from the rest of the system. Regions are defined, not in real space, but in Hilbert space, though connection between the two can be made by spatial localization of single-particle orbitals. Relations between our embedding approach and some similar techniques are discussed. We present results for small sample systems for which we can extract essentially exact results, demonstrating that our approach seems to work quite well and is generally more reliable than some of the related approaches due to the introduction of additional interaction terms.     

Benchmark density functional theory calculations for nanoscale conductance

We present a set of benchmark calculations for the Kohn-Sham elastic transmission function of five representative single-molecule junctions. The transmission functions are calculated using two different density functional theory methods, namely an ultrasoft pseudopotential plane-wave code in combination with maximally localized Wannier functions and the norm-conserving pseudopotential code SIESTA which applies an atomic orbital basis set. All calculations have been converged with respect to the supercell size and the number of k|| points in the surface plane. For all systems we find that the SIESTA transmission functions converge toward the plane-wave result as the SIESTA basis is enlarged. Overall, we find that an atomic basis with double zeta and polarization is sufficient (and in some cases, even necessary) to ensure quantitative agreement with the plane-wave calculation. We observe a systematic downshift of the SIESTA transmission functions relative to the plane-wave results. The effect diminishes as the atomic orbital basis is enlarged; however, the convergence can be rather slow.     

Recent applications of density functional theory calculations to biomolecules

The purpose of this overview is to highlight the broad scope and utility of current applications of density functional theory (DFT) methods for the study of the properties and reactions of biomolecules. This is illustrated using examples selected from research carried out within our research group and in collaboration with others. The examples include the hyperfine coupling constants of amino acid radicals, the use of an amino acid as a chiral catalyst for the formation of carbon–carbon bonds in the aldol reaction, hydrogen-bond mediated catalysis of an aminolysis reaction, radiation-induced protein–DNA cross-links, and the mechanism by which an antitumor drug cleaves DNA. We demonstrate that DFT-based methods can be applied successfully to a broad range of problems that remain beyond the scope of conventional electron-correlation methods. Furthermore, we show that contemporary computational quantum chemistry complements experiment in the study of biological systems. Received: 19 December 2001 / Accepted: 8 April 2002 / Published online: 4 July 2002     

Studies of iridium nanoparticles using density functional theory calculations

The energetics and the electronic and magnetic properties of iridium nanoparticles in the range of 2-64 atoms were investigated using density functional theory calculations. A variety of different geometric configurations were studied, including planar, three-dimensional, nanowire, and single-walled nanotube. The binding energy per atom increases with size and dimensionality from 2.53 eV/atom for the iridium dimer to 6.09 eV/atom for the 64-atom cluster. The most stable geometry is planar until four atoms are reached and three-dimensional thereafter. The simple cubic structure is the most stable geometric building block until a strikingly large 48-atom cluster, when the most stable geometry transitions to face-centered cubic, as found in the bulk metal. The strong preference for cubic structure among small clusters demonstrates their rigidity. This result indicates that iridium nanoparticles intrinsically do not favor the coalescence process. Nanowires formed from linear atomic chains of up to 4-atom rings were studied, and the wires formed from 4-atom rings were extremely stable. Single-walled nanotubes were also studied. These nanotubes were formed by stacking 5- and 6-atom rings to form a tube. The ring stacking with each atom directly above the previous atom is more stable than if the alternate rings are rotated.     

Methylation of zebularine investigated using density functional theory calculations

Deoxyribonucleic acid (DNA) methylation is an epigenetic phenomenon, which adds methyl groups into DNA. This study reveals methylation of a nucleoside antibiotic drug 1‐(β‐D ‐ribofuranosyl)‐2‐pyrimidinone (zebularine or zeb) with respect to its methylated analog, 1‐(β‐D ‐ribofuranosyl)‐5‐methyl‐2‐pyrimidinone (d5) using density functional theory calculations in valence electronic space. Very similar infrared spectra suggest that zeb and d5 do not differ by types of the chemical bonds, but distinctly different Raman spectra of the nucleoside pair reveal that the impact caused by methylation of zeb can be significant. Further valence orbital‐based information details on valence electronic structural changes caused by methylation of zebularine. Frontier orbitals in momentum space and position space of the molecules respond differently to methylation. Based on the additional methyl electron density concentration in d5, orbitals affected by the methyl moiety are classified into primary and secondary contributors. Primary methyl contributions include MO8 (57a), MO18 (47a), and MO37 (28a) of d5, which concentrates on methyl and the base moieties, suggest certain connection to their Frontier orbitals. The primary and secondary methyl affected orbitals provide useful information on chemical bonding mechanism of the methylation in zebularine. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Performance of density functional theory in computing nonresonant vibrational (hyper)polarizabilities

A set of exchange‐correlation functionals, including BLYP, PBE0, B3LYP, BHandHLYP, CAM‐B3LYP, LC‐BLYP, and HSE, has been used to determine static and dynamic nonresonant (nuclear relaxation) vibrational (hyper)polarizabilities for a series of all‐trans polymethineimine (PMI) oligomers containing up to eight monomer units. These functionals are assessed against reference values obtained using the Møller–Plesset second‐order perturbation theory (MP2) and CCSD methods. For the smallest oligomer, CCSD(T) calculations confirm the choice of MP2 and CCSD as appropriate for assessing the density functionals. By and large, CAM‐B3LYP is the most successful, because it is best for the nuclear relaxation contribution to the static linear polarizability, intensity‐dependent refractive index second hyperpolarizability, static second hyperpolarizability, and is close to the best for the electro‐optical Pockels effect first hyperpolarizability. However, none of the functionals perform satisfactorily for all the vibrational (hyper)polarizabilities studied. In fact, in the case of electric field‐induced second harmonic generation all of them, as well as the Hartree–Fock approximation, yield the wrong sign. We have also found that the Pople 6–31+G(d) basis set is unreliable for computing nuclear relaxation (hyper)polarizabilities of PMI oligomers due to the spurious prediction of a nonplanar equilibrium geometry. © 2013 Wiley Periodicals, Inc.     

Atomic multiplet energies from density functional calculations

Atomic multiplet term energies for dⁿ configurations have been estimated within density functional theory (DFT) exploiting symmetry to the largest possible extent. The electrostatic two-electron integrals, as well as term energies, are expressed in function of only three non-redundant single determinants (NRSDs), each of them being obtained from density functional calculations. The influence of correlation effects described with a gradient-corrected functional (GGA) is examined and discussed. Comparison with experimental data shows the reliability of this symmetry-based density functional approach.     

Prediction of cellulose nanotube models through density functional theory calculations

We report the generation of a nano-scale tubular structure of cellulose molecules (CelNT), through density functional theory (DFT) calculations. When a cellulose III_I (1 0 0) chain sheet model is optimized by DFT calculations, the sheet models spontaneously roll into tubes. The oligomers arrange in a right-handed, four-fold helix with one-quarter chain staggering, oriented with parallel polarity similar to the original crystal structure. Based on a one-quarter chain staggering relationship, six large CelNT models, consisting of 16 cellulose chains with DP = 80, are constructed by combinations of two types of chain polarities and three types of symmetry operations to generate a circular arrangement of molecular chains. All six CelNT models are examined by molecular dynamics (MD) calculations in chloroform. While four CelNT models retain a tubular form throughout MD calculations, the remaining two deform. 3D-RISM theory model is used to estimate the solvation free energies of the four CelNT models. The results suggest that the CelNT model with a chain arrangement of parallel polarity and right-handed helical symmetry forms the most stable tube structure.     

Spectral analysis of acetylcholine halides by density functional theory calculations

The optimized molecular structures, vibrational frequencies and ¹H and ¹³C NMR chemical shifts of acetylcholine halides (F, Cl, and Br) have been investigated using density functional theory (B3LYP) method with 6-311G(d) basis set. The comparison of their experimental and calculated IR, R and NMR spectra of the compounds has indicated that the spectra of three optimized minimum energy conformers can simultaneously exist in one experimental spectrum. Thus, it was concluded that the compounds simultaneously exist in three conformations in the ground state. The calculated optimized geometric parameters (bond lengths and bond angles), vibrational frequencies and NMR chemical shifts for the minimum energy conformers were seen to be in a good agreement with the corresponding experimental data. All the assignments of the theoretical frequencies were performed by potential energy distributions using VEDA 4 program.     

Time-dependent density functional theory calculations of the photoabsorption of fluorinated alkanes

Time-dependent density functional theory (TD-DFT) calculations of the transition energies and oscillator strengths of fluorinated alkanes have been performed. The TD-DFT method with the non-local B3LYP potential yields transition energies for the methanes, which are smaller by about 10% as compared to the experimental values. An empirical linear correlation was found between the calculated and experimental transition energies both at the B3LYP/DZ+Ryd(C, F) and B3LYP/cc-pVTZ+Ryd(C, F, H) levels for a total of 19 transitions of the fluorinated methanes with linear correlation coefficients of 0.987 for the former and 0.988 for the latter. This empirical correlation for fluorinated methane molecules is found to agree well with the previously obtained empirical correlations between calculated and experimental values for non-fluorinated molecules. The results show that a single empirical-correlation relationship can be used for both non-fluorinated and fluorinated molecules to predict transition energies. This linear relationship is then used to predict the photoabsorption spectra of ethane, propane, butane, and partially and fully fluorinated derivatives. A key result of these calculations is the dominance of Rydberg transitions in the spectral region of interest.     

Matrix infrared spectra and density functional theory calculations of molybdenum hydrides

Laser-ablated Mo atoms react with H2 upon condensation in excess argon, neon, and hydrogen. The molybdenum hydrides MoH, MoH2, MoH4, and MoH6 are identified by isotopic substitution (H2, D2, HD, H2 + D2) and by comparison with vibrational frequencies calculated by density functional theory. The MoH2 molecule is bent, MoH4 is tetrahedral, and MoH6 appears to have the distorted trigonal prism structure.     

Projected gradient algorithms for Hartree-Fock and density matrix functional theory calculations

We present projected gradient algorithms designed for optimizing various functionals defined on the set of N-representable one-electron reduced density matrices. We show that projected gradient algorithms are efficient in minimizing the Hartree-Fock or the Muller-Buijse-Baerends functional. On the other hand, they converge very slowly when applied to the recently proposed BBk (k=1,2,3) functionals [O. Gritsenko et al., J. Chem. Phys. 122, 204102 (2005)]. This is due to the fact that the BBk functionals are not proper functionals of the density matrix.     

Robust acceleration of self consistent field calculations for density functional theory

We show that the type 2 Broyden secant method is a robust general purpose mixer for self consistent field problems in density functional theory. The Broyden method gives reliable convergence for a large class of problems and parameter choices. We directly mix the approximation of the electronic density to provide a basis independent mixing scheme. In particular, we show that a single set of parameters can be chosen that give good results for a large range of problems. We also introduce a spin transformation to simplify treatment of spin polarized problems. The spin transformation allows us to treat these systems with the same formalism as regular fixed point iterations.