- Search
- Download PDF

(May 16, 2021)

If we have a wave function $|\mathrm{\Psi}\u27e9={\displaystyle \sum _{i}}{C}_{i}|{D}_{i}\u27e9$ (where $|{D}_{i}\u27e9$ are Slater determinants with coefficients ${C}_{i}$) as an eigenstate of the Hamiltonian, then

$${C}_{i}=\frac{{\sum}_{j\ne i}{H}_{ij}{C}_{j}}{({H}_{ii}-E)},$$ | (6.75) |

where ${H}_{ij}=\u27e8{D}_{i}|\mathbf{H}|{D}_{j}\u27e9$ is the Hamiltonian matrix element between determinants $i$ and $j$, and $E$ is the energy of the eigenstate $|\mathrm{\Psi}\u27e9$. This exact relationship can be generalized to a metric to predict the expected weight of a determinant $|{D}_{i}\u27e9$ in a CI expansion, by how it connects to other determinants in an approximate trial wave function. This metric is also used in Epstein-Nesbet Perturbation theory as coefficients for the determinants in the first order wave function.

In the ASCI method, all determinants $|{D}_{i}\u27e9$ that are single or double excitations away from the most important determinants (as ranked by magnitude of coefficients) in the trial wave function $|{\psi}_{k}\u27e9$ are assigned an estimated importance ${A}_{i}$ given as

$${A}_{i}=\frac{{\sum}_{|{D}_{j}\u27e9\in |{\psi}_{k}\u27e9}{H}_{ij}{C}_{j}}{({H}_{ii}-{E}_{k})},$$ | (6.76) |

where ${E}_{k}$ is the energy of the trial wave function $|{\psi}_{k}\u27e9$. The search and selection is only done in the space spanned by determinants connected to the top $c$ determinants in $|{\psi}_{k}\u27e9$ because unimportant determinants are unlikely to be the sole generator for a top ranked determinant, and this pruning of the search space greatly accelerates the algorithm. The top $t$ determinants (as ranked by magnitude of ${A}_{i}$) connected to $|{\psi}_{k}\u27e9$ are used to form the new wave function $|{\psi}^{k+1}\u27e9$ by exact diagonalization within that Hilbert subspace.

Once several cycles of ASCI has been completed, the wave function will contain all (or
very nearly all) of the largest weight determinants in the FCI wave function and the
remaining determinants not included should be of small weight. The effect of these many
small remaining determinants are estimated by second order Epstein-Nesbet perturbation
theory (PT2).
^{
298
}
Phys. Rev.

(1926),
28,
pp. 695.
Link
^{,}
^{
790
}
Proc. Roy. Soc. Ser. A

(1955),
230,
pp. 312.
Link
This final PT2 correction gives extremely
accurate results, often within a kcal/mol of the absolute FCI energies even when only a
tiny fraction of the Hilbert space is included in the ASCI wave function.
^{
1107
}
J. Chem. Phys.

(2016),
145,
pp. 044112.
Link
An extrapolation of the variational energy against the PT2 correction (to the FCI limit
of of zero PT2 correction) can also be carried out to generate more accurate estimates,
and predict a metric for error in the final estimate. Indeed, it has been shown
that linear or quadratic fits are quite accurate for extrapolation of SCI energies
against the PT2 correction.
^{
399
}
J. Chem. Theory Comput.

(2019),
15,
pp. 5370.
Link
^{,}
^{
17
}
J. Chem. Phys.

(2017),
147,
pp. 164111.
Link
^{,}
^{
697
}
J. Chem. Theory Comput.

(2018),
14,
pp. 4360.
Link
We observe essentially linear behavior in the case of ASCI.

ASCI may be used as the full-CI solver for a CASSCF calculation, permitting the extension
of CASSCF to active spaces of $\approx 50$ electrons in $\approx 50$ orbitals. The resulting
method is termed ASCI-SCF
^{
654
}
J. Chem. Theory Comput.

(2020),
Link
. See section 6.19 for details on
CASSCF job control.