next up previous contents index
Next: Excited State Geometry Optimizations Up: How to Perform Previous: Stability Analysis   Contents   Index

Vertical Excitation and CD Spectra

The calculation of excited states within the TDHF(RPA)/TDDFT approach is enabled by
$scfinstab rpas

for closed-shell singlet excitations,
$scfinstab rpat

for closed-shell triplet excitations, and
$scfinstab urpa

for excitations out of spin-unrestricted reference states.
If it is intended to use the TDA instead, specify
$scfinstab ciss

for closed-shell singlet excitations,
$scfinstab cist

for closed-shell triplet excitations,
$scfinstab ucis

for excitations out of spin-unrestricted reference states, and
$scfinstab spinflip

for spin-flip (z -component of the total spin changes by ±1 ) excitations out of spin-unrestricted reference states. For details concerning the theory see ref. [89]. In practice, this functionality can be used for the calculation of triplet-singlet, quartet-doublet, ... excitations (see ref. [90] also for further information about the implementation). It is only available within the TDA in combination with LDA functionals and the HF exchange. It is strongly recommended to increase $escfiterlimit.

Next, the IRREPs of the excitations need to be defined, which is again accomplished using $soes. For example, to calculate the 17 lowest excitations in IRREP b1g, the 23 lowest excitations in IRREP eu, and all excitations in IRREP t2g, use

b1g  17
eu   23
t2g  all
and run escf.

Note that $soes specifies the IRREP of the excitation vector which is not necessarily identical to the IRREP of the excited state(s) involved. In general, the IRREP(s) of the excitation(s) from the ground to an excited state is given by the direct product of the IRREPs of the tow states. For example, to calculate the first A2 state in a C2v -symmetric molecule with a B2 (open-shell) ground state, it is necessary to specify

b1   1

The number of excitations that have to be calculated in order to cover a certain spectral range is often difficult to determine in advance. The total number of excitations within each IRREP as provided by the define ex menu may give some hint. A good strategy is to start with a smaller number of excitations and, if necessary, perform a second escf run on a larger number of states using the already converged excitation vectors as input.

To compute absorption and CD spectra, it is often sufficient to include optically allowed transitions only. This leads to substantial reduction of computational effort for molecules with higher symmetry. For example, in the UV-VIS spectrum of an Oh symmetric molecule, only t1u excitations are optically allowed. The IRREPs of the electric and magnetic dipole moments as well as of the electric quadrupole moment are displayed automatically in the define ex menu.

If a large number of states is to be calculated, it is highly recommended to provide extra memory by specifying

$rpacor m
the integer m being the core memory size in megabytes (default is 20). The larger m, the more vectors can be processed simultaneously without re-calculation of integrals. As a rule of thumb, m should be ca. 90% of the available main memory. If RI-J is used ($ridft), it is recommended to set $ricore to a small value and $rpacor to a large value if the number of states is large, and vice versa if it is small.

By specifying

$spectrum unit
$cdspectrum unit
a list of excitation energies and oscillator and/or rotatory strengths of the optically allowed transitions is written onto file spectrum and/or cdspectrum. As above, unit specifies the energy unit and may be ev, nm, 1/cm, or a.u. (default). The files spectrum and cdspectrum may conveniently be used for further processing, e.g., using a plotting program such as Gnuplot.

By specifying

inclusion of the current-density response for MGGA calculations is disabled. Note that the results of calculations using this flag will no longer be gauge-invariant and will differ from results obtained with the standard gauge-invariant implementation.

next up previous contents index
Next: Excited State Geometry Optimizations Up: How to Perform Previous: Stability Analysis   Contents   Index