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.     

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     

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     

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.     

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.     

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.     

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.     

Theoretical Auger electron spectra of polymers by density functional theory calculations using model dimers

We propose a new approach for analysis of Auger electron spectra (AES) of polymers by density functional theory (DFT) calculations with the Slater's transition-state concept. Simulated AES and X-ray photoelectron spectra (XPS) of four polymers [(CH2CH2)n (PE), (CH2CH(CH3))n (PP), (CH2CH(OCH3))n (PVME), and (CH2CH(COCH3))n (PVMK)] by DFT calculations using model dimers are in a good accordance with the experimental ones. The experimental AES of the polymers can be classified in each range of 1s-2p2p, 1s-2s2p, and 1s-2s2s transitions for C KVV and O KVV spectra, and in individual contributions of the functional groups from the theoretical analysis.     

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.     

Surface-enhanced Raman scattering of 4,4′-bipyridine on silver by density functional theory calculations

Surface-enhanced Raman scattering (SERS) of 4,4′-bipyridine (BPy) on silver foil substrate was measured using the 488, 514.5, and 1064 nm excitation lines. Density functional theory (DFT) methods were used to calculate the structure and vibrational spectra of Ag–BPy, Ag₃–BPy and Ag₄–BPy complexes with B3LYP/6-31++G(d,p)(C,H,N)/Lanl2dz(Ag) basis set. The Raman bands of BPy were assigned on the basis of the calculation of potential energy distribution. The calculated spectra of Ag–BPy and Ag₄–BPy complexes were much closer to the experimental results of BPy adsorbed on silver surface than that of Ag₃–BPy complexes. The vibrational frequencies that are sensitive to the planar or non-planar structure of BPy and to the dihedral angle of two pyridyl rings were discussed. The DFT results showed that the angles between two pyridyl rings for Ag–BPy and Ag₄–BPy were skewed by about 38.44° and 37.1°, respectively. The energy gaps of the HOMO and LUMO from DFT were 415–912 nm for BPy–Ag complexes. The relative intensities of SERS bands changed with different excitation laser lines. Thus, a chemical enhancement mechanism should play an important role in the SERS of BPy on silver substrate.     

Basis set dependence of phosphate frequencies in density functional theory calculations

The addition of extravalence, polarization and diffuse functions, were studied in order to conclude how they affect the P? O stretching frequencies of several biological relevant phosphate molecules. The results show that the polarization and the diffuse functions have opposite effects on the frequencies: the polarization functions downshift while the diffuse functions upshift the frequencies. The effect of the valence functions was more difficult to interpret. The effect of the conductor‐like screening model (CPCM)‐continuum model was also studied. The results show that the CPCM‐continuum model has a substantial effect on the frequencies for these small molecules. The continuum model's efficiency is mainly due to its effect on the geometries and not on the frequencies. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Basis set convergence of explicitly correlated double-hybrid density functional theory calculations

The basis set convergence of explicitly correlated double-hybrid density functional theory (DFT) is investigated using the B2GP-PLYP functional. As reference values, we use basis set limit B2GP-PLYP-F12 reaction energies extrapolated from the aug(')-cc-pV(Q+d)Z and aug(')-cc-pV(5+d)Z basis sets. Explicitly correlated double-hybrid DFT calculations converge significantly faster to the basis set limit than conventional calculations done with basis sets saturated up to the same angular momentum (typically, one "gains" one angular momentum in the explicitly correlated calculations). In explicitly correlated F12 calculations the VnZ-F12 basis sets converge faster than the orbital A(')VnZ basis sets. Furthermore, basis set convergence of the MP2-F12 component is apparently faster than that of the underlying Kohn-Sham calculation. Therefore, the most cost-effective approach consists of combining the MP2-F12 correlation energy from a comparatively small basis set such as VDZ-F12 with a DFT energy from a larger basis set such as aug(')-cc-pV(T+d)Z.     

Kinetic energies to analyze the experimental auger electron spectra by density functional theory calculations

In the Auger electron spectra (AES) simulations, we define theoretical modified kinetic energies of AES in the density functional theory (DFT) calculations. The modified kinetic energies correspond to two final-state holes at the ground state and at the transition-state in DFT calculations, respectively. This method is applied to simulate Auger electron spectra (AES) of 2nd periodic atom (Li, Be, B, C, N, O, F)-involving substances (LiF, beryllium, boron, graphite, GaN, SiO₂, PTFE) by deMon DFT calculations using the model molecules of the unit cell. Experimental KVV (valence band electrons can fill K-shell core holes or be emitted during KVV-type transitions) AES of the (Li, O) atoms in the substances agree considerably well with simulation of AES obtained with the maximum kinetic energies of the atoms, while, for AES of LiF, and PTFE substance, the experimental F KVV AES is almost in accordance with the spectra from the transitionstate kinetic energy calculations.     

A characterization study on 2,6-dimethyl-4-nitropyridine N-oxide by density functional theory calculations

This study deals with the identification of a title compound, 2,6-dimethyl-4-nitropyridine N-oxide by means of theoretical calculations. The optimized molecular structures, vibrational frequencies, corresponding vibrational assignments, thermodynamic properties and atomic charges of the title compound in the ground state were evaluated using density functional theory (DFT) with the standard B3LYP/6-311G(d,p) method and basis set combination for the first time. Theoretical vibrational spectra were interpreted with the aid of normal coordinate analysis based on scaled density functional force field. The results show that the optimized geometric parameters (bond lengths and bond angles) and vibrational frequencies were observed to be in good agreement with the available experimental results. Based on the results of comparison between experimental results and theoretical data, the chosen calculation level is powerful approach for understanding the molecular structures and vibrational spectra of the 2,6-dimethyl-4-nitropyridine N-oxide. Moreover, we not only simulated frontier molecular orbitals (FMO) and molecular electrostatic potential (MEP) but also determined the transition state and energy band gap. Based on the investigations, the title compound is found to be useful to bond metallically and interact intermolecularly. Infrared intensities and Raman activities were also reported.     

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.     

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.     

Periodic Hartree-Fock and density functional theory calculations for Li-doped polyacetylene chains

We have performed periodic restricted Hartree-Fock/6-31G** and B3LYP6-31G** density functional theory calculations on Li-doped trans-polyacetylene at various dopant concentrations, using C(2m)H(2m)Li2 unit cells (m = 7-14). Except for maintaining P1 rod symmetry the geometry was completely optimized for both uniform and nonuniform doping structures. In addition to geometry we obtain atomic charges, along with soliton formation and dopant binding energies, as well as band structures and densities of states. A thorough analysis of the band structure and density of states, as a function of dopant concentration, is presented. We also characterize the complex nature of the binding interaction between Li and the polyacetylene chain.     

Condensed phase ionic polarizabilities from plane wave density functional theory calculations

A method is presented to allow the calculation of the dipole polarizabilities of ions and molecules in a condensed-phase coordination environment. These values will be useful for understanding the optical properties of materials and for developing simulation potentials which incorporate polarization effects. The reported values are derived from plane wave density functional theory calculations, though the method itself will apply to first-principles calculations on periodic systems more generally. After reporting results of test calculations on atoms to validate the procedure, values for the polarizabilities of the oxide ion and various cations in a range of materials are reported and compared with experimental information as well as previous theoretical results.