We present a new coupled Hartree–Fock(HF)/Kohn–Sham DFT perturbation method that accounts for the effect of enlarging the basis set in electronic structure calculations. In contrast with previous approaches, our dual basis set treatment yields not only a correction for the total energy but also for the orbital eigenvalues and density. The zeroth order solution is obtained from the projection of the small basis set coefficients. Diagonalization of the full Fock matrix in the large basis set is avoided. In this first paper of a series, we develop the theoretical foundations of our approach for molecules, including the coupled-perturbed equations through second order and the energy expressions through fourth order—as our method complies with Wigner’s 2n + 1 rule. The first-order perturbation equation turns out to be uncoupled, and odd-order terms in the energy expansion vanish. In calculations on simple molecules, our method recovers over 93% (84%) of the missing DFT(HF) energy when going from the cc-pVDZ to the aug-cc-pVDZ basis, and over about 95% in all cases if an energy extrapolation formula is used. Mulliken charges, the orbital eigenvalue spectrum, and HOMO–LUMO gaps of the large basis are well reproduced. Charge density maps show that the differences between the perturbatively corrected density and the reference nearly vanish through second-order.

Coupled Perturbation Theory Approach to Dual Basis Sets for Molecules and Solids. 1. General Theory and Application to Molecules

Maschio, Lorenzo;Kirtman, Bernard
2019-01-01

Abstract

We present a new coupled Hartree–Fock(HF)/Kohn–Sham DFT perturbation method that accounts for the effect of enlarging the basis set in electronic structure calculations. In contrast with previous approaches, our dual basis set treatment yields not only a correction for the total energy but also for the orbital eigenvalues and density. The zeroth order solution is obtained from the projection of the small basis set coefficients. Diagonalization of the full Fock matrix in the large basis set is avoided. In this first paper of a series, we develop the theoretical foundations of our approach for molecules, including the coupled-perturbed equations through second order and the energy expressions through fourth order—as our method complies with Wigner’s 2n + 1 rule. The first-order perturbation equation turns out to be uncoupled, and odd-order terms in the energy expansion vanish. In calculations on simple molecules, our method recovers over 93% (84%) of the missing DFT(HF) energy when going from the cc-pVDZ to the aug-cc-pVDZ basis, and over about 95% in all cases if an energy extrapolation formula is used. Mulliken charges, the orbital eigenvalue spectrum, and HOMO–LUMO gaps of the large basis are well reproduced. Charge density maps show that the differences between the perturbatively corrected density and the reference nearly vanish through second-order.
2019
16
1
340
353
Maschio, Lorenzo; Kirtman, Bernard
File in questo prodotto:
File Dimensione Formato  
dual_basis_mol_revision.pdf

Open Access dal 04/12/2020

Descrizione: Articolo in versione post-print
Tipo di file: POSTPRINT (VERSIONE FINALE DELL’AUTORE)
Dimensione 1.93 MB
Formato Adobe PDF
1.93 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2318/1725527
Citazioni
  • ???jsp.display-item.citation.pmc??? 1
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 7
social impact