Search

The electric field gradient (EFG) at the nucleus, the generalized Sternheimer shielding constants, and the EFG hyperpolarizabilities of a set of reference molecules are computed using analytic density-functional theory (DFT) quadratic response theory. At the three-parameter Becke--parameter-Lee--Yang--Parr (B3LYP) level, DFT underestimates correlation effects compared with accurate coupled-cluster and multiconfigurational self-consistent field estimates. For the prediction of EFG properties of hydrogen nuclei and electron-rich atoms such as halides, B3LYP even provides less reliable results than Hartree--Fock theory.

We present an analysis of the dispersion interaction energy and forces in density-functional theory from the point of view of the adiabatic connection between the Kohn–Sham non-interacting and fully interacting systems. Accurate coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] densities are computed for the helium dimer and used to construct the exchange-correlation potential of Kohn–Sham theory, showing agreement with earlier results presented for the Hartree–Fock–Kohn–Sham method [M. Allen and D. J. Tozer, J. Chem. Phys. 117, 11113 (2002)]. The accuracy of the methodology utilized to determine these solutions is checked by calculation of the Hellmann–-Feynman forces based on the Kohn–-Sham densities, which are compared with analytic CCSD(T) forces. To ensure that this comparison is valid in a finite atomic-orbital basis set, we employ floating Gaussian basis functions throughout and all results are counterpoise corrected. The subtle charge-rearrangement effects associated with the dispersion interaction are highlighted as the origin of a large part of the dispersion force. To recover the exchange-correlation components of the interaction energy, adiabatic connections are constructed for the supermolecular system and for its constituent atoms; subtraction of the resulting adiabatic-connection curves followed by integration over the interaction strength recovers the exchange-correlation contribution relevant to the density-functional description of the dispersion interaction. The results emphasize the long-ranged, dynamically correlated nature of the dispersion interaction between closed-shell species. An alternative adiabatic-connection path is also explored, where the electronic interactions are introduced in a manner that emphasizes the range of the electronic interactions, highlighting their purely long-ranged nature, consistent with the success of range-separated hybrid approaches in this context.

Using a recently implemented technique for the calculation of the adiabatic connection (AC) of density functional theory (DFT) based on Lieb maximization with respect to the external potential, the AC is studied for atoms and molecules containing up to ten electrons: the helium isoelectronic series, the hydrogen molecule, the beryllium isoelectronic series, the neon atom, and the water molecule. The calculation of AC curves by Lieb maximization at various levels of electronic-structure theory is discussed. For each system, the AC curve is calculated using Hartree–Fock (HF) theory, second-order Møller–Plesset (MP2) theory, coupled-cluster singles-and-doubles (CCSD) theory, and coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] theory, expanding the molecular orbitals and the effective external potential in large Gaussian basis sets. The HF AC curve includes a small correlation-energy contribution in the context of DFT, arising from orbital relaxation as the electron-electron interaction is switched on under the constraint that the wave function is always a single determinant. The MP2 and CCSD AC curves recover the bulk of the dynamical correlation energy and their shapes can be understood in terms of a simple energy model constructed from a consideration of the doubles-energy expression at different interaction strengths. Differentiation of this energy expression with respect to the interaction strength leads to a simple two-parameter doubles model (AC-D) for the AC integrand (and hence the correlation energy of DFT) as a function of the interaction strength. The structure of the triples-energy contribution is considered in a similar fashion, leading to a quadratic model for the triples correction to the AC curve (AC-T). From a consideration of the structure of a two-level configuration-interaction (CI) energy expression of the hydrogen molecule, a simple two-parameter CI model (AC-CI) is proposed to account for the effects of static correlation on the AC. When parametrized in terms of the same input data, the AC-CI model offers improved performance over the corresponding AC-D model, which is shown to be the lowest-order contribution to the AC-CI model. The utility of the accurately calculated AC curves for the analysis of standard density functionals is demonstrated for the BLYP exchange-correlation functional and the interaction-strength-interpolation (ISI) model AC integrand. From the results of this analysis, we investigate the performance of our proposed two-parameter AC-D and AC-CI models when a simple density functional for the AC at infinite interaction strength is employed in place of information at the fully interacting point. The resulting two-parameter correlation functionals offer a qualitatively correct behavior of the AC integrand with much improved accuracy over previous attempts. The AC integrands in the present work are recommended as a basis for further work, generating functionals that avoid spurious error cancellations between exchange and correlation energies and give good accuracy for the range of densities and types of correlation contained in the systems studied here.

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.

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field theories for the calculation of frequency-dependent molecular response properties and excitation energies is presented, based on a nonredundant exponential parametrization of the one-electron density matrix in the atomic-orbital basis, avoiding the use of canonical orbitals. The response equations are solved iteratively, by an atomic-orbital subspace method equivalent to that of molecular-orbital theory. Important features of the subspace method are the use of paired trial vectors (to preserve the algebraic structure of the response equations), a nondiagonal preconditioner (for rapid convergence), and the generation of good initial guesses (for robust solution). As a result, the performance of the iterative method is the same as in canonical molecular-orbital theory, with five to ten iterations needed for convergence. As in traditional direct Hartree-Fock and Kohn-Sham theories, the calculations are dominated by the construction of the effective Fock/Kohn-Sham matrix, once in each iteration. Linear complexity is achieved by using sparse-matrix algebra, as illustrated in calculations of excitation energies and frequency-dependent polarizabilities of polyalanine peptides containing up to 1400 atoms.

The relationship between the random-phase-approximation (RPA) correlation energy and the continuous algebraic Riccati equation is examined and the importance of a stabilizing solution is emphasized. The criterion to distinguish this from non-stabilizing solutions can be used to ensure that physical, smooth potential energy surfaces are obtained. An implementation of analytic RPA molecular gradients is presented using the Lagrangian technique. Illustrative calculations indicate that RPA with Hartree-Fock reference orbitals delivers an accuracy similar to that of second-order Mo̸ller-Plesset perturbation theory.