Search Results


9.11 Ab Initio Molecular Dynamics with Complex Absorbing Potentials

9.11.2 Finding Electronic Resonance States of Temporary Anions

(May 7, 2024)

We limit our discussion here to the HF method, but the general principles outlined also applies to higher level methods.

  1. 1.

    The correct description of the electronic state of temporary anions usually requires the use of extra diffuse functions in the basis set. You may thus consider having enough extra diffuse functions in your basis set.

  2. 2.

    It is best you use the core Hamiltonian as guess for the CAP-HF SCF procedure (i.e., set SCF_GUESS = CORE).

  3. 3.

    For CAP calculations, Q-Chem first solves the SCF problem using the old SCF drivers in Q-Chem; the solution is then used as the starting point for the CAP part of the calculation. This preliminary SCF solution usually leads to wrong CAP-HF solutions. You may avoid this issue by setting SKIP_OLD_SCFMAN = TRUE in $rem. See Example 9.11.3.

  4. 4.

    Make sure the CAP contribution to the real and imaginary parts of the SCF solution found is negligible. You can get these CAP contributions by asking for an energy decomposition of the complex energy by setting CS_SCF_FINAL_PRINT = 1 in $rem. An alternative is to look at how close the imaginary (and real) parts of the CAP-HF energy and the corrected CAP-HF energy are; the ratio of the imaginary parts should be close to unity. See Example 9.11.3.

For CAP-HF calculations, set the $rem variable CS_STRICT = TRUE to print out the correct properties (e.g., Mulliken charges and multipole moments) of the solution. See Example 9.11.3.