An approximate density-functional method using the Harris-Foulkes functional

We present a method which uses the results of a molecular Kohn-Sham calculation at a reference geometry to approximate the energy at many different geometries. The Kohn-Sham electron density of the reference geometry is decomposed into atomic fragments, which move with the nuclei to approximate the density at a new geometry and the energy is evaluated with the Harris-Foulkes functional. Preliminary results for a biological quantum-mechanics/molecular-mechanics trajectory are promising: the errors of reference-geometry Harris-Foulkes (compared to full self-consistent Kohn-Sham) for the PBE exchange-correlation functional have the same magnitude as the difference between the energies of PBE and BLYP.     

Time-dependent density-functional theory in the projector augmented-wave method

We present the implementation of the time-dependent density-functional theory both in linear-response and in time-propagation formalisms using the projector augmented-wave method in real-space grids. The two technically very different methods are compared in the linear-response regime where we found perfect agreement in the calculated photoabsorption spectra. We discuss the strengths and weaknesses of the two methods as well as their convergence properties. We demonstrate different applications of the methods by calculating excitation energies and excited state Born-Oppenheimer potential surfaces for a set of atoms and molecules with the linear-response method and by calculating nonlinear emission spectra using the time-propagation method.     

Time-dependent density-functional theory/localized density matrix method for dynamic hyperpolarizability

Time-dependent density-functional theory/localized density matrix method (TDDFT/LDM) was developed to calculate the excited state energy, absorption spectrum and dynamic polarizability. In the present work we generalize it to calculate the dynamic hyperpolarizabilities in both time and frequency domains. We show that in the frequency domain the 2n+1 rule can be derived readily and the dynamic hyperpolarizabilities are thus calculated efficiently. Although the time-domain TDDFT/LDM is time consuming, its implementation is straightforward because the evaluation of the derivatives of exchange-correlation potential with respect to electron density is avoided. Moreover, the time-domain method can be used to simulate higher order response which is very difficult to be calculated with the frequency-domain method.     

Calculation of optical rotation with time-periodic magnetic-field-dependent basis functions in approximate time-dependent density-functional theory

We report the implementation of a method for the calculation of optical rotation. This method is based on the time-dependent density-functional theory and utilizes time-periodic magnetic-field-dependent basis functions. The calculations are based on a density fit. It is demonstrated that additional terms in the analytical expression appearing from derivatives of the approximated Coulomb potential are necessary to provide the gauge-origin independence of the results within a given numerical accuracy. Contributions from these terms also restore the symmetry between the electric and magnetic perturbations in the optical rotation tensor.     

The trust-region self-consistent field method in Kohn-Sham density-functional theory

The trust-region self-consistent field (TRSCF) method is extended to the optimization of the Kohn-Sham energy. In the TRSCF method, both the Roothaan-Hall step and the density-subspace minimization step are replaced by trust-region optimizations of local approximations to the Kohn-Sham energy, leading to a controlled, monotonic convergence towards the optimized energy. Previously the TRSCF method has been developed for optimization of the Hartree-Fock energy, which is a simple quadratic function in the density matrix. However, since the Kohn-Sham energy is a nonquadratic function of the density matrix, the local energy functions must be generalized for use with the Kohn-Sham model. Such a generalization, which contains the Hartree-Fock model as a special case, is presented here. For comparison, a rederivation of the popular direct inversion in the iterative subspace (DIIS) algorithm is performed, demonstrating that the DIIS method may be viewed as a quasi-Newton method, explaining its fast local convergence. In the global region the convergence behavior of DIIS is less predictable. The related energy DIIS technique is also discussed and shown to be inappropriate for the optimization of the Kohn-Sham energy.     

Fluctuation-dissipation theorem density-functional theory

Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)-that is, the first such results with fxc not equal to 0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.     

A subsystem density-functional theory approach for the quantum chemical treatment of proteins

We present an extension of the frozen-density embedding (FDE) scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] that can be applied to subsystems connected by covalent bonds, as well as a practical implementation of such an extended FDE scheme. We show how the proposed scheme can be employed for quantum chemical calculations of proteins by treating each constituting amino acid as a separate subsystem. To assess the accuracy of the extended FDE scheme, we present calculations for several dipeptides and for the protein ubiquitin.     

Time-dependent density-functional theory for open electronic systems

Time-dependent density-functional theory(TDDFT)has been successfully applied to predict excited-state properties of isolated and periodic systems.However,it cannot address a system coupled to an environment or whose number of electrons is not conserved.To tackle these problems,TDDFT needs to be extended to accommodate open systems.This paper provides a comprehensive account of the recent developments of TDDFT for open systems(TDDFT-OS),including both theoretical and practical aspects.The practicality and accuracy of a latest TDDFT-OS method is demonstrated with two numerical examples:the time-dependent electron transport through a series of quasi-one-dimensional atomic chains,and the real-time electronic dynamics on a two-dimensional graphene surface.The advancement of TDDFT-OS may lead to promising applications in various fields of chemistry,including energy conversion and heterogeneous catalysis.     

Real-time linear response for time-dependent density-functional theory

We present a linear-response approach for time-dependent density-functional theories using time-adiabatic functionals. The resulting theory can be performed both in the time and in the frequency domain. The derivation considers an impulsive perturbation after which the Kohn-Sham orbitals develop in time autonomously. The equation describing the evolution is not strictly linear in the wave function representation. Only after going into a symplectic real-spinor representation does the linearity make itself explicit. For performing the numerical integration of the resulting equations, yielding the linear response in time, we develop a modified Chebyshev expansion approach. The frequency domain is easily accessible as well by changing the coefficients of the Chebyshev polynomial, yielding the expansion of a formal symplectic Green's operator.     

A multiconfigurational hybrid density-functional theory

We propose a multiconfigurational hybrid density-functional theory which rigorously combines a multiconfiguration self-consistent-field calculation with a density-functional approximation based on a linear decomposition of the electron-electron interaction. This gives a straightforward extension of the usual hybrid approximations by essentially adding a fraction λ of exact static correlation in addition to the fraction λ of exact exchange. Test calculations on the cycloaddition reactions of ozone with ethylene or acetylene and the dissociation of diatomic molecules with the Perdew-Burke-Ernzerhof and Becke-Lee-Yang-Parr density functionals show that a good value of λ is 0.25, as in the usual hybrid approximations. The results suggest that the proposed multiconfigurational hybrid approximations can improve over usual density-functional calculations for situations with strong static correlation effects.     

Generator coordinate method in time-dependent density-functional theory: memory made simple

The generator coordinate (GC) method is a variational approach to the quantum many-body problem in which interacting many-body wave functions are constructed as superpositions of (generally nonorthogonal) eigenstates of auxiliary Hamiltonians containing a deformation parameter. This paper presents a time-dependent extension of the GC method as a new approach to improve existing approximations of the exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT). The time-dependent GC method is shown to be a conceptually and computationally simple tool to build memory effects into any existing adiabatic XC potential. As an illustration, the method is applied to driven parametric oscillations of two interacting electrons in a harmonic potential (Hooke's atom). It is demonstrated that a proper choice of time-dependent generator coordinates in conjunction with the adiabatic local-density approximation reproduces the exact linear and nonlinear two-electron dynamics quite accurately, including features associated with double excitations that cannot be captured by TDDFT in the adiabatic approximation.     

Double-hybrid density-functional theory made rigorous

We provide a rigorous derivation of a class of double-hybrid approximations, combining Hartree-Fock exchange and second-order M?ller-Plesset correlation with a semilocal exchange-correlation density functional. These double-hybrid approximations contain only one empirical parameter and use a density-scaled correlation energy functional. Neglecting density scaling leads to a one-parameter version of the standard double-hybrid approximations. We assess the performance of these double-hybrid schemes on representative test sets of atomization energies and reaction barrier heights, and we compare to other hybrid approximations, including range-separated hybrids. Our best one-parameter double-hybrid approximation, called 1DH-BLYP, roughly reproduces the two parameters of the standard B2-PLYP or B2GP-PLYP double-hybrid approximations, which shows that these methods are not only empirically close to an optimum for general chemical applications but are also theoretically supported.     

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.     

Inverse density-functional theory as an interpretive tool for measuring colloid-surface interactions in dense systems

Recent advances in optical microscopy, such as total internal reflection and confocal scanning laser techniques, now permit the direct three-dimensional tracking of large numbers of colloidal particles both near and far from interfaces. A novel application of this technology, currently being developed by one of the authors under the name of diffusing colloidal probe microscopy (DCPM), is to use colloidal particles as probes of the energetic characteristics of a surface. A major theoretical challenge in implementing DCPM is to obtain the potential energy of a single particle in the external field created by the surface, from the measured particle trajectories in a dense colloidal system. In this paper we develop an approach based on an inversion of density-functional theory (DFT), where we calculate the single-particle-surface potential from the experimentally measured equilibrium density profile in a nondilute colloidal fluid. The underlying DFT formulation is based on the recent work of Zhou and Ruckenstein [Zhou and Ruckenstein, J. Chem. Phys. 112, 8079 (2000)]. For model hard-sphere and Lennard-Jones systems, using Monte Carlo simulation to provide the "experimental" density profiles, we found that the inversion procedure reproduces the true particle-surface-potential energy to an accuracy within typical DCPM experimental limitations (approximately 0.1 kT) at low to moderate colloidal densities. The choice of DFT closures also significantly affects the accuracy.     

A dual-level approach to density-functional theory

An efficient approximate scheme for density-functional theory (DFT) calculations, which eliminates the time-consuming self-consistent-field (SCF) procedure, is proposed using a dual-level DFT approach. In this approach, dual levels of basis sets and exchange-correlation functionals are adopted. The dual-level DFT approach is based on the idea that the total electron density in the ground state can be represented in terms of the density evaluated using the low-quality basis set and the low-cost exchange-correlation functional. Since the SCF procedure is avoided in the total energy evaluation, the dual-level DFT approach drastically reduces the computational cost. The applications of several dual-level DFT calculations to molecular systems show that our approach is more efficient than the self-consistent DFT approach with a moderate accuracy.     

A dual-level state-specific time-dependent density-functional theory

A highly efficient new algorithm for time-dependent density-functional theory (TDDFT) calculations is presented. In this algorithm, a dual-level approach to speed up DFT calculations (Nakajima and Hirao, J Chem Phys 2006, 124, 184108) is combined with a state-specific (SS) algorithm for TDDFT (Chiba et al., Chem Phys Lett 2006, 420, 391). The dual-level SS-TDDFT algorithm was applied to excitation energy calculations of typical small molecules, the Q bands of the chlorophyll A molecule, the charge-transfer energy of the zincbacteriochlorin-bacteriochlorin model system, and the lowest-lying excitation of the circumcoronene molecule. As a result, it was found that the dual-level SS-TDDFT gave correct excitation energies with errors of 0.2-0.3 eV from the standard TDDFT approach, with much lower CPU times for various types of excitation energies of large-scale molecules.     

On the theory underlying the Car-Parrinello method and the role of the fictitious mass parameter

The theory underlying the Car-Parrinello extended-Lagrangian approach to ab initio molecular dynamics (CPMD) is reviewed and reexamined using "heavy" ice as a test system. It is emphasized that the adiabatic decoupling in CPMD is not a decoupling of electronic orbitals from the ions but only a decoupling of a subset of the orbital vibrational modes from the rest of the necessarily coupled system of orbitals and ions. Recent work [J. Chem. Phys. 116, 14 (2002)] has pointed out that, due to the orbital-ion coupling that remains once adiabatic decoupling has been achieved, a large value of the fictitious mass mu can lead to systematic errors in the computed forces in CPMD. These errors are further investigated in the present work with a focus on those parts of these errors that are not corrected simply by rescaling the masses of the ions. It is suggested that any comparison of the efficiencies of Born-Oppenheimer molecular dynamics (BOMD) and CPMD should be performed at a similar level of accuracy. If accuracy is judged according to the average magnitude of the systematic errors in the computed forces, the efficiency of BOMD compares more favorably to that of CPMD than previous comparisons have suggested.