We introduce a density-based multilevel Hartree-Fock (HF) method where the electronic density is optimized in a given region of the molecule (the active region). Active molecular orbitals (MOs) are generated by a decomposition of a starting guess atomic orbital (AO) density, whereas the inactive MOs (which constitute the remainder of the density) are never generated or referenced. The MO formulation allows for a significant dimension reduction by transforming from the AO basis to the active MO basis. All interactions between the inactive and active regions of the molecule are retained, and an exponential parametrization of orbital rotations ensures that the active and inactive density matrices separately, and in sum, satisfy the symmetry, trace, and idempotency requirements. Thus, the orbital spaces stay orthogonal, and furthermore, the total density matrix represents a single Slater determinant. In each iteration, the (level-shifted) Newton equations in the active MO basis are solved to obtain the orbital transformation matrix. The approach is equivalent to variationally optimizing only a subset of the MOs of the total system. In this orbital space partitioning, no bonds are broken and no a priori orbital assignments are carried out. In the limit of including all orbitals in the active space, we obtain an MO density-based formulation of full HF.

Density-Based Multilevel Hartree-Fock Model

Koch, Henrik;
2017

Abstract

We introduce a density-based multilevel Hartree-Fock (HF) method where the electronic density is optimized in a given region of the molecule (the active region). Active molecular orbitals (MOs) are generated by a decomposition of a starting guess atomic orbital (AO) density, whereas the inactive MOs (which constitute the remainder of the density) are never generated or referenced. The MO formulation allows for a significant dimension reduction by transforming from the AO basis to the active MO basis. All interactions between the inactive and active regions of the molecule are retained, and an exponential parametrization of orbital rotations ensures that the active and inactive density matrices separately, and in sum, satisfy the symmetry, trace, and idempotency requirements. Thus, the orbital spaces stay orthogonal, and furthermore, the total density matrix represents a single Slater determinant. In each iteration, the (level-shifted) Newton equations in the active MO basis are solved to obtain the orbital transformation matrix. The approach is equivalent to variationally optimizing only a subset of the MOs of the total system. In this orbital space partitioning, no bonds are broken and no a priori orbital assignments are carried out. In the limit of including all orbitals in the active space, we obtain an MO density-based formulation of full HF.
2017
Settore CHIM/02 - Chimica Fisica
Computer Science Applications1707 Computer Vision and Pattern Recognition; Physical and Theoretical Chemistry
File in questo prodotto:
File Dimensione Formato  
acs.jctc.7b00689.pdf

Accesso chiuso

Tipologia: Published version
Licenza: Non pubblico
Dimensione 3.62 MB
Formato Adobe PDF
3.62 MB Adobe PDF   Richiedi una copia

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/11384/69634
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 32
  • ???jsp.display-item.citation.isi??? 32
  • OpenAlex ND
social impact