*QUADRATIC RESPONSE¶
This section gives directives for the calculation of quadratic response functions [Saue2002a].
General control statements¶
.PRINT¶
Print level.
Default:
.PRINT
0
Definition of the quadratic response function¶
.DIPLEN¶
Specification of dipole operators for A, B, and C (see One-electron operators for details).
.A OPERATOR¶
Specification of the A operator (see One-electron operators for details).
.B OPERATOR¶
Specification of the B operator (see One-electron operators for details).
.C OPERATOR¶
Specification of the C operator (see One-electron operators for details).
.B FREQ¶
Specify frequencies of operator B.
Example: 3 different frequencies.
.B FREQ
3
0.001
0.002
0.01
Default: Static case.
.B FREQ
1
0.0
.C FREQ¶
Specify frequencies of operator C (see .B FREQ).
.ALLCMB¶
Evaluate all nonzero quadratic response functions and thereby disregarding analysis of overall permutational symmetry.
Default: Evaluate only unique, nonzero, response functions.
Excited state properties¶
This page describes unreleased functionality. The keywords may not be available in your version of DIRAC.
First order properties of excited states can be computed from the quadratic response function.
.EXCPRP¶
Give the number of “left” and “right” states in each boson symmetry.
Example:
.EXCPRP
5 5 5 5
0 0 0 0
Compute the excited state expectation values |\langle i|\hat{A}|i\rangle|, where i goes from 1 to 5 in each symmetry (four symmetries in this case). The zeros can be substituted with positive integers to generate transition state moments |\langle i|\hat{A}|j\rangle|.
Control variational parameters¶
.SKIPEE¶
Exclude all rotations between occupied positive-energy and virtual positive-energy orbitals.
.SKIPEP¶
Exclude all rotations between occupied positive-energy and virtual negative-energy orbitals.
Control reduced equations¶
.MAXITR¶
Maximum number of iterations.
Default:
.MAXITR
30
.MAXRED¶
Maximum dimension of matrix in reduced system.
Default:
.MAXRED
100
.THRESH¶
Threshold for convergence of reduced system.
Default:
.THRESH
1.0D-5
Control integral contributions¶
The user is encouraged to experiment with these options since they may have an important effect on run time.
.INTFLG¶
Specify what two-electron integrals to include (default: .INTFLG under **HAMILTONIAN).
.CNVINT¶
Set threshold for convergence before adding SL and SS integrals to SCF-iterations.
2 (real) Arguments:
.CNVINT
CNVXQR(1) CNVXQR(2)
Default: Very large numbers.
.ITRINT¶
Set the number of iterations before adding SL and SS integrals to SCF-iterations.
Default:
.ITRINT
1 1
Control trial vectors¶
.XQRNRM¶
Normalize trial vectors. Using normalized trial vectors will reduce efficiency of screening.
Default: Use un-normalized vectors.
Advanced/debug flags¶
.NOPREC¶
No preconditioning of initial trial vectors.
Default: Preconditioning of trial vectors.
.RESFAC¶
New trial vector will be generated only for variational parameter classes whose residual has a norm that is larger than a fraction 1/RESFAC of the maximum norm.
Default:
.RESFAC
1000.0