Time-dependent density-functional theory calculations of triplet-triplet absorption

We present density-functional theory calculations of triplet-triplet absorption by three different approaches based on time-dependent density-functional theory (DFT): unrestricted DFT linear response, open-shell restricted DFT linear response applied to the triplet state, and quadratic response with triplet excitations applied to the ground state. Comparison is also made with corresponding results obtained by Hartree-Fock and multiconfiguration self-consistent-field response theory. Two main conclusions concerning triplet-triplet transitions are drawn in this study: First, the very good agreement between unrestricted and restricted DFT results indicates that spin contamination of the triplet state is not a serious problem when computing triplet-triplet spectra of common organic molecules. Second, DFT response calculations of triplet-triplet transitions can be affected by triplet instability problems, especially for the combination of DFT quadratic response with functionals containing fractional exact Hartree-Fock exchange.     

Energy gradients with respect to atomic positions and cell parameters for the Kohn-Sham density-functional theory at the Gamma point

The application of theoretical methods based on density-functional theory is known to provide atomic and cell parameters in very good agreement with experimental values. Recently, construction of the exact Hartree-Fock exchange gradients with respect to atomic positions and cell parameters within the Gamma-point approximation has been introduced. In this article, the formalism is extended to the evaluation of analytical Gamma-point density-functional atomic and cell gradients. The infinite Coulomb summation is solved with an effective periodic summation of multipole tensors. While the evaluation of Coulomb and exchange-correlation gradients with respect to atomic positions are similar to those in the gas phase limit, the gradients with respect to cell parameters needs to be treated with some care. The derivative of the periodic multipole interaction tensor needs to be carefully handled in both direct and reciprocal space and the exchange-correlation energy derivative leads to a surface term that has its origin in derivatives of the integration limits that depend on the cell. As an illustration, the analytical gradients have been used in conjunction with the QUICCA algorithm to optimize one-dimensional and three-dimensional periodic systems at the density-functional theory and hybrid Hartree-Fock/density-functional theory levels. We also report the full relaxation of forsterite supercells at the B3LYP level of theory.     

Efficient hybrid density functional calculations in solids: assessment of the Heyd-Scuseria-Ernzerhof screened Coulomb hybrid functional

The present work introduces an efficient screening technique to take advantage of the fast spatial decay of the short range Hartree-Fock (HF) exchange used in the Heyd-Scuseria-Ernzerhof (HSE) screened Coulomb hybrid density functional. The screened HF exchange decay properties and screening efficiency are compared with traditional hybrid functional calculations on solids. The HSE functional is then assessed using 21 metallic, semiconducting, and insulating solids. The examined properties include lattice constants, bulk moduli, and band gaps. The results obtained with HSE exhibit significantly smaller errors than pure density functional theory (DFT) calculations. For structural properties, the errors produced by HSE are up to 50% smaller than the errors of the local density approximation, PBE, and TPSS functionals used for comparison. When predicting band gaps of semiconductors, we found smaller errors with HSE, resulting in a mean absolute error of 0.2 eV (1.3 eV error for all pure DFT functionals). In addition, we present timing results which show the computational time requirements of HSE to be only a factor of 2-4 higher than pure DFT functionals. These results make HSE an attractive choice for calculations of all types of solids.     

Ultrasoft pseudopotentials in time-dependent density-functional theory

We describe an efficient formulation allowing the use of ultrasoft pseudopotentials (USPPs) in plane wave based time-dependent density-functional theory. The practical steps required to implement USPP functionality within real time propagation schemes and linear-response schemes based on Lanczos algorithms are provided. The functioning of the methodology is demonstrated by calculations of the optical absorption spectra of the fullerene C(60), using both real time propagation and the Lanczos/linear-response approaches. Comparisons between the rates of convergence of the optical spectra with the number of applications of the Hamiltonian required in calculations with ultrasoft pseudopotentials and norm-conserving pseudopotentials show clearly the benefits provided by the use of USPP.     

Multiconfiguration optimized effective potential method for a density-functional treatment of static correlation

An approach to treat static correlation within a density-functional framework is presented. To that end, a multiconfiguration optimized effective potential (MCOEP) method is derived. In contrast to standard multiconfiguration self-consistent field (MCSCF) methods and previous combinations of MCSCF procedures with density-functional theory, the MCOEP method yields well-defined physically meaningful orbital and eigenvalue spectra. In addition to the electronic ground state also excited electronic states can be described. The MCOEP method is implemented invoking the localized Hartree-Fock approximation, leading to a multiconfiguration localized Hartree-Fock approach. Applications of the new method to the dissociation of the hydrogen molecule and the isomerization of ethene and cyclobutadiene show that it is capable of describing situations that are characterized by strong static correlation.     

Multipole-based integral estimates for the rigorous description of distance dependence in two-electron integrals

We derive multipole-based integral estimates (MBIE) as rigorous and tight upper bounds to four-center two-electron integrals in order to account for the 1/R distance decay between the charge distributions, which is missing in the Schwarz screening commonly used in ab initio methods. Our screening criteria are valid for all angular momenta and can be formulated for any order of multipoles. We have found the expansion limited to dipoles to be sufficiently tight for estimating the integrals in Hartree-Fock and density-functional theories, while the screening effort is negligible. For, e.g., a DNA fragment with 1052 atoms and 10,674 basis functions (6-31G*) the exchange part is faster by a factor of 2.1 as compared to the Schwarz screening both within our linear exchange scheme, whereas a smaller factor of 1.3 is gained for the Coulomb part within the continuous fast multipole method. Most importantly, our new MBIE screening is perfectly suited to exploit the strong distance decay of electron-correlation effects of at least 1/R4 in atomic-orbital-based formulations of correlation methods.     

Excitonic effects in a time-dependent density functional theory

Excited state properties of one-dimensional molecular materials are dominated by many-body interactions resulting in strongly bound confined excitons. These effects cannot be neglected or treated as a small perturbation and should be appropriately accounted for by electronic structure methodologies. We use adiabatic time-dependent density functional theory to investigate the electronic structure of one-dimensional organic semiconductors, conjugated polymers. Various commonly used functionals are applied to calculate the lowest singlet and triplet state energies and oscillator strengths of the poly(phenylenevinylene) and ladder-type (poly)(para-phenylene) oligomers. Local density approximations and gradient-corrected functionals cannot describe bound excitonic states due to lack of an effective attractive Coulomb interaction between photoexcited electrons and holes. In contrast, hybrid density functionals, which include long-range nonlocal and nonadiabatic corrections in a form of a fraction of Hartree-Fock exchange, are able to reproduce the excitonic effects. The resulting finite exciton sizes are strongly dependent on the amount of the orbital exchange included in the functional.     

The planar-to-tubular structural transition in boron clusters from optical absorption

The optical response of the lowest-energy isomers of the B20 family is calculated using time-dependent density-functional theory within a real-space, real-time scheme. Significant differences are found among the absorption spectra of the clusters studied. We show that these differences can be easily related to changes in the overall geometry. Optical spectroscopy is thus an efficient tool to characterize the planar-to-tubular structural transition, known to be present in these boron-based systems.     

Open-shell localized Hartree-Fock method based on the generalized adiabatic connection Kohn-Sham formalism for a self-consistent treatment of excited states

An effective exact-exchange Kohn-Sham approach for the treatment of excited electronic states, the generalized adiabatic connection open-shell localized Hartree-Fock (GAC-OSLHF) method is presented. The GAC-OSLHF method is based on the generalized adiabatic connection Kohn-Sham formalism and therefore capable of treating excited electronic states, which are not the energetically lowest of their symmetry. The method is self-interaction free and allows for a fully self-consistent computation of excited valence as well as Rydberg states. Results for atoms and small- and medium-size molecules are presented and compared to restricted open-shell Hartree-Fock (ROHF) and time-dependent density-functional results as well as to experimental data. While GAC-OSLHF and ROHF results are quite close to each other, the GAC-OSLHF method shows a much better convergence behavior. Moreover, the GAC-OSLHF method as a Kohn-Sham method, in contrast to the ROHF approach, represents a framework which allows also for a treatment of correlation besides an exchange by appropriate functionals. In contrast to the common time-dependent density-functional methods, the GAC-OSLHF approach is capable of treating doubly or multiply excited states and can be easily applied to molecules with an open-shell ground state. On the nodal planes of the energetically highest occupied orbital, the local multiplicative GAC-OSLHF exchange potential asymptotically approaches a different, i.e., nonzero, value than in other regions, an asymptotic behavior which is known from exact Kohn-Sham exchange potentials of ground states of molecules.     

On the role of the nonlocal Hartree-Fock exchange in ab initio quantum transport: H2 in Pt nanocontacts revisited

We propose a practical way to overcome the ubiquitous problem of the overestimation of the zero-bias and zero-temperature conductance, which is associated with the use of local approximations to the exchange-correlation functional in density-functional theory when applied to quantum transport. This is done through partial substitution of the local exchange term in the functional by the nonlocal Hartree-Fock exchange. As a nontrivial example of this effect we revisit the smallest molecular bridge studied so far: a H(2) molecule placed in between Pt nanocontacts. When applied to this system the value of the conductance diminishes as compared to the local-exchange-only value, which is in close agreement with the results predicted from time-dependent current-density-functional theory. Our results issue a warning message on recent claims of perfect transparency of a H(2) molecule in Pt nanocontacts.     

Study of electronic and spectroscopic properties on a newly synthesized red fluorescent material

The ground state (S(0)) and the lowest singlet excited state (S(1)) of a newly synthesized red fluorescent material, 2-[3-(2-{4-[(2-Hydroxy-ethyl)-methyl-amino]-phenyl}-vinyl)-5,5-dimethyl-cyclohex-2-enylidene]-malononitrile (A31), are investigated. The S(0) and S(1) geometries are optimized at the ab initio Hartree-Fock and the singles configuration interaction (CIS) levels of theory, respectively. The CIS and semiempirical Zerner's Intermediate Neglect of Differential Overlap (ZINDO) methods provide the results for the absorption (S(0)-->S(1)) and emission (S(1)-->S(0)) transition energies. The Stokes shifts calculated at the CIS and ZINDO levels of theory are obtained. The absorption spectra in various solvents are calculated using the time-dependent density-functional theory method in combination with the polarized continuum model, which are in very good agreement with our experimental measurements. The solvent effects are discussed.     

Assessment of dressed time-dependent density-functional theory for the low-lying valence states of 28 organic chromophores

Almost all time-dependent density-functional theory (TDDFT) calculations of excited states make use of the adiabatic approximation, which implies a frequency-independent exchange-correlation kernel that limits applications to one-hole/one-particle states. To remedy this problem, Maitra et al. [N.T. Maitra, F. Zhang, R.J. Cave, K. Burke, Double excitations within time-dependent density functional theory linear response theory, J. Chem. Phys. 120 (2004) 5932 ] proposed dressed TDDFT (D-TDDFT), which includes explicit two-hole/two-particle states by adding a frequency-dependent term to adiabatic TDDFT. This paper offers the first extensive test of D-TDDFT, and its ability to represent excitation energies in a general fashion. We present D-TDDFT excited states for 28 chromophores and compare them with the benchmark results of Schreiber et al. [M. Schreiber, M.R. Silva-Junior, S.P.A. Sauer, W. Thiel, Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3, J. Chem. Phys. 128 (2008) 134110]. We find the choice of functional used for the A-TDDFT step to be critical for positioning the 1h1p states with respect to the 2h2p states. We observe that D-TDDFT without HF exchange increases the error in excitations already underestimated by A-TDDFT. This problem is largely remedied by implementation of D-TDDFT including Hartree-Fock exchange.     

Hartree-Fock and Kohn-Sham time-dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling

We present a second-quantization based atomic-orbital method for the computation of time-dependent response functions within Hartree-Fock and Kohn-Sham density-functional theories. The method is suited for linear scaling. Illustrative results are presented for excitation energies, one- and two-photon transition moments, polarizabilities, and hyperpolarizabilities for hexagonal BN sheets with up to 180 atoms.     

Cubic response functions in time-dependent density functional theory

We present density-functional theory for time-dependent response functions up to and including cubic response. The working expressions are derived from an explicit exponential parametrization of the density operator and the Ehrenfest principle, alternatively, the quasienergy ansatz. While the theory retains the adiabatic approximation, implying that the time-dependency of the functional is obtained only implicitly-through the time dependence of the density itself rather than through the form of the exchange-correlation functionals-it generalizes previous time-dependent implementations in that arbitrary functionals can be chosen for the perturbed densities (energy derivatives or response functions). In particular, general density functionals beyond the local density approximation can be applied, such as hybrid functionals with exchange correlation at the generalized-gradient approximation level and fractional exact Hartree-Fock exchange. With our implementation the response of the density can always be obtained using the stated density functional, or optionally different functionals can be applied for the unperturbed and perturbed densities, even different functionals for different response order. As illustration we explore the use of various combinations of functionals for applications of nonlinear optical hyperpolarizabilities of a few centrosymmetric systems; molecular nitrogen, benzene, and the C(60) fullerene. Considering that vibrational, solvent, and local field factors effects are left out, we find in general that very good experimental agreement can be obtained for the second dynamic hyperpolarizability of these systems. It is shown that a treatment of the response of the density beyond the local density approximation gives a significant effect. The use of different functional combinations are motivated and discussed, and it is concluded that the choice of higher order kernels can be of similar importance as the choice of the potential which governs the Kohn-Sham orbitals.     

A computational study of some electric and magnetic properties of gaseous BF3 and BCl3

We present the results of an extended computational study of the electric and magnetic properties connected to Cotton-Mouton birefringences, on the trifluoro- and trichloroborides in the gas phase. The electric dipole polarizabilities, magnetizabilities, quadrupole moments, and higher-order hypersusceptibilities--expressed as quadratic and cubic frequency-dependent response functions--are computed within Hartree-Fock, density-functional, and coupled-cluster response theories employing singly and doubly augmented correlation-consistent basis sets and London orbitals in the magnetic property calculations. The results, which illustrate the capability of time-dependent density-functional theory for electron-rich systems, are compared with available experimental data. Revised values of both experimentally derived quadrupole moment of BF3, 2.72 +/- 0.15 a.u., and magnetizability anisotropy of BCl3, -0.45 +/- 0.09 a.u., both obtained in birefringence experiments that neglect the effects of higher-order hypersusceptibilities, are presented. In the theoretical limit the traceless quadrupole moments of BF3 and BCl3 are determined to be 3.00 +/- 0.01 and 0.71 +/- 0.01 a.u., respectively.     

Efficient exact exchange approximations in density-functional theory

Two approaches to approximate the Slater potential component of local exact exchange of density-functional theory are investigated. The first approach employs density fitting of the electrostatic potential integrals over two occupied orbitals and the other approach approximates the "exact" Slater potential with the potential derived from the Becke-Roussel [Phys. Rev. A. 39, 3761 (1989)] model of the exchange hole. In both cases significant time savings can be achieved for larger systems compared to the calculation of the numerical Slater potential. It is then analyzed how well the orbitals obtained from the various total exchange potentials reproduce Hartree-Fock energies and molecular properties. A large range of atoms and small molecules has been utilized, including the three DNA bases adenine, thymine, and cytosine.     

Communication: Monte Carlo calculation of the exchange energy

In recent generalized Kohn-Sham (GKS) schemes for density functional theory (DFT) Hartree-Fock type exchange is important. In plane waves and grid approaches the high cost of exchange energy calculations makes these GKS considerably more expensive than Kohn-Sham DFT calculations. We develop a stochastic approach for speeding up the calculation of exchange for large systems. We show that stochastic error per particle does not grow and can even decrease with system size (at a given number of iterations). We discuss several alternative approaches and explain how these ideas can be included in the GKS framework.     

Perturbational calculations of parity-violating effects in nuclear-magnetic-resonance parameters

We investigate the effects of the parity-violating electroweak interaction in the spectral parameters of nuclear magnetic resonance. Perturbational theory of parity-violating effects in the nuclear magnetic shielding is presented to the order of G(F)alpha, and in the indirect spin-spin coupling, to the order of G(F)alpha3. These leading-order parity-violating corrections are evaluated using analytical linear-response theory methods based on Hartree-Fock and density-functional theory reference states. Parity-violating contributions to spin-spin couplings are evaluated for the first time at the first-principles level. Calculations are carried out for two chiral halomethanes, bromochlorofluoromethane and bromofluoroiodomethane.