manual   quickstart   instguide   update   basis

Next: 33.3 Parameters on LAPLACE Up: 33 LOCAL METHODS FOR Previous: 33.1 Local CC2 and   Contents   Index   PDF

33.2 Options for EOM

see also section 25.1.

In the case of LT-DF-LCC2/ADC(2) multistate calculations are possible, and it is recommended to calculate more states as needed.

The parameters on the EOM card:

EOM,-n.1, key1=value1, key2=value2,...

where n.1 is the last state of interest, e.g., with EOM,-5.1 the four lowest excited states will be calculated. The following keywords key are possible:

SINGLET=ising
If set to 1, singlet states will be calculated (default).
TRIPLET=itrip
If set to 1, triplet states will be calculated (not implemented for ADC(2)).
EXFILE=record.ifil
Record for converged CIS eigenvectors (default 6100.2).
SAVE=record.ifil
Record for save of restart information.
SAVET=record.ifil
Record for restart information for triplet states (if calculated together with singlet states).
START=record.ifil
Record for restart of previous calculation.
STARTT=record.ifil
Record for restart of previous calculation (if calculated together with singlet states).
NSRCH4ST=nst
In the first nst iterations in the Davidson diagonalisation the excited state domains are determined for each basis vector in the Davidson subspace ("search for states") (default 7).
THRLCH=thrlch
threshold for the Davidson procedure. If smaller than zero, the Davidson procedure is skipped and DIIS is started directly instead (possible only for restart, SAVE and START have to be identical).
THRLCD=thrlcd
threshold for DIIS.
ADC2=adc2
if $=1$, do ADC(2) calculation instead of CC2.
LTGS=ltgs
if $=1$, use LT-DF-LMP2 for the ground state.

Default local approximations are defined according to procedure described in Ref. [3] (Laplace domains).

INTFRAC=fracint
Rough criterion for specifying eom-domains from laplace-transformed integrals (default 0.8).
INTEXC=excint
Criterion for specifying important orbitals from laplace-transformed integrals (default 0.999).
REALFRAC=fracreal
Exact criterion for specifying eom-domains from laplace-transformed integrals (default 0.98).
FULLFRAC=fracfull
Check for all orbital domains (complete sum over all orbitals) (default 0.95).

To switch to local aproximations calculated according to Ref. [6] (Boughton-Pulay procedure for excited states), set INTFRAC to zero.

Occupied orbital pair lists are calculated from important orbitals (Refs. [1,5]).

EWEAKPAIR=ewpair
Distance criterion for excited state orbital pairs definition (default 5).
ALLSTRONG=allstr
Different possibilities for excited state orbital pairs:
0: $[ij]\leq$ ewpair, $[im]\leq$ ewpair, $[mn]\leq$ ewpair;
1: $\forall [ij]$, $\forall [im]$, $[mn]\leq$ ewpair;
2: $\forall [ij]$, $[im]\leq$ ewpair, $[mn]\leq$ ewpair (default);
where $i$, $j$ are important orbitals and $m$, $n$ non-important. More detailed see Ref. [1].

One can try to improve the convergence of iterative procedures by changing following parameters

PRECOND=prec
Type of preconditioner in LT-DF-LCC2/ADC(2) and DF-CIS:
0: canonical orbital energies;
3: solving linear equations including diagonal part of $(H-F)^{\rm CIS}$ with MINRES (default).
HSCAL=hscl
Scaling factor for diagonal part of $(H-F)^{\rm CIS}$ matrix in the case of linear equations preconditioner (default 0.7).

Properties are also activated on the EOM card:

PROPES=prop
States for which properties (dipole moments) should be calculated, e.g., PROPES=-3.1+5.1-8.1+15.1 $\Longrightarrow$ 2.1 3.1 5.1 6.1 7.1 8.1 15.1
TRANES=tran
States for which transition moments should be calculated, syntax like for properties. Is not implemented for ADC(2) and for triplet states.
DENSAVE=rec.ifil
Record, where densities calculated for states prop can be saved (for printing, see 33.4)



Next: 33.3 Parameters on LAPLACE Up: 33 LOCAL METHODS FOR Previous: 33.1 Local CC2 and   Contents   Index   PDF

manual   quickstart   instguide   update   basis

molpro@molpro.net 2013-05-24