35.7 Choosing the ansatz and the level of approximation
The Ansatz can be chosen using the ANSATZ option and/or by options on the command line.
- Ansatz 3A;
- Ansatz 3A with EBC approximation
- Ansatz 3B
- Ansatz 3B with EBC approximation
- Ansatz 3C
- Ansatz 3C with EBC approximation
The ansatz can be further detailed by appending options in parenthesis, e.g.
These options can be one of
- Use diagonal ansatz
- Use diagonal ansatz with fixed amplitudes (orbital invariant)
- Use diagonal ansatz with fixed amplitudes and canonical orbitals
- Use diagonal ansatz and assume X-matrix to be diagonal (only for Ansatz 3A)
- Use GBC approximation (only for 3B, default in 3A)
- Use EBC approximation (same as *)
- Use HX1 approximation (only for 3B).
- Use HY1 approximation (only for 3B and 3C).
- Use HY2 approximation (only for 3B and 3C).
- Default hybrid approximation. HX1 and HY2 approximation in 3B, HY2 in 3C.
- Neglect Z terms (only for 3B and 3C).
- Neglect X terms (only for 3A and 3B).
Several options separated by commas can be given. For instance
uses the diagonal ansatz 3C(D) with fixed coefficients and the hybrid (HY1) approximation.
Alternatively or in addition, the following options can be given on the command line:
- Use diagonal ansatz.
- Use diagonal ansatz and assume X-matrix to be diagonal.
- Use GBC approximation (only for 3B, default in 3A).
- Use EBC approximation (same as *).
- Use HYn approximation.
- Use HX1 approximation.
- Neglect Z terms (only 3B and 3C).
- Neglect X terms (only 3A and 3B).
- Use diagonal ansatz with fixed coefficient approximation (orbital invariant).
- Use diagonal ansatz with fixed coefficient approximation. Evaluate only first order energy expression,
not the Hylleraas functional. Very fast but less accurate and reliable!
- Use diagonal ansatz with fixed coefficient approximation and canonical orbitals. A non-iterative method
is used to evaluate the energy. This is equivalent to FIX=1,CANONICAL=1 and is most efficient.
- As FIXCAN=1, but equations are solved iteratively (test purpose only).
- Use CABS (default). If CABS=0 is given, CABS is disabled. However, if RI_BASIS=OPTRI,
the orbital and OPTRI basis sets are automatically merged, and then exactly the same results as with CABS=1 are obtained.
- Construct CABS basis from orthogonal MOs and ABS basis rather than AO and RI basis.
- Threshold for smallest eigenvalue of S in auxiliary ABS (only used with ORTHO_CABS=1; default=THRCABS).
- Threshold for smallest eigenvalue of S in CABS (default 1.d-8).
- Relative CABS threshold (default 1.d-9). The actual threshold is
max(thrcabs,eigmax*thrcabs_rel, where eigmax is the largest eigenvalue of the overlap matrix.
- Print parameter. PRINT=1 give information about all computed integrals
and the iterations.
- Can be used to obtain extended debug print.
- Use a most efficient pair-specific fully iterative method (default).
- Use simple fully iterative method.
- Use pair specific iterative method (more expensive).
- Use pair specific non iterative method (most expensive, only with canonical orbitals).
- Use canonical orbitals and full domains.
- Use full domains and localized orbitals (unless CANONICAL=1 is given).
- Scale triples energy as explained in section 35.10.
- If set to one, include CABS singles correction (default=1)
- If set to one, include CABS singles correction for core orbitals (default=0)
- For open-shell systems: If 1 (default,recommended), include all occupied valence orbital
pairs for in , independent of spin (as
described in J. Chem. Phys. 130, 054104 (2009), section II.E). If 0, use only pairs
mn where the spins of and , and and are equal.
implies a canonical 3C calculation with diagonal ansatz 3C, using fixed coefficient and hynrid approximations. The
combination of the options fix=1 and canonical=1 implies a non-iterative calculation of the energy
and is recommended. The above is equivalent to all of the following:
Note that the HF convergence threshold should be rather strict to obtain
accurate results (use ACCU,14 in the HF).
Numerous further options are for specialist use only and not described here. See explicit.registry for a full list.