TD
TDA
DESCRIPTION
This method keyword requests an excited state calculation using the timedependent HartreeFock or DFT method [Bauernschmitt96a, Casida98, Stratmann98, VanCaillie99, VanCaillie00, Furche02, Scalmani06]; analytic gradients are available in Gaussian 09 [Furche02, Scalmani06]. Timedependent DFT calculations can employ the TammDancoff approximation, via the TDA keyword. TDDFTB calculations can also be performed [Trani11].
Note that the normalization criteria used is <X+YXY>=1.
Electronic circular dichroism (ECD) analysis is also performed during these calculations [Helgaker91, Bak93, Bak95, Olsen95, Hansen99, Autschbach02].
GENERAL OPTIONS
Singlets
Solve only for singlet excited states. Only effective for closedshell systems, for which it is the default.
Triplets
Solve only for triplet excited states. Only effective for closedshell systems.
5050
Solve for half triplet and half singlet states. Only effective for closedshell systems.
Root=N
Specifies the “state of interest”. The default is the first excited state (N=1).
NStates=M
Solve for M states (the default is 3). If 5050 is requested, NStates gives the number of each type of state for which to solve (i.e., the default is 3 singlets and 3 triplets).
Add=N
Read converged states off the checkpoint file and solve for an additional N states. This option implies Read as well.
Read
Reads initial guesses for the states off the checkpoint file. Note that, unlike for SCF, an initial guess for one basis set cannot be used for a different one.
EqSolv
Whether to perform equilibrium or nonequilibrium PCM solvation. NonEqSolv is the default except for excited state optimizations and when the excited state density is requested (e.g., with Density=Current or All).
IVOGuess
Force use of IVO guess. This is the default for TD HartreeFock. NoIVOGuess forces the use of canonical single excitations for guess, and it is the default for TDDFT. The HFIVOGuess option forces the use of HartreeFock IVOs for the guess, even for TDDFT.
SOS
Do sumover states polarizabilities, etc. By default, all excited states are solved for. A list of frequencies at which to do the sums is read in. Zero frequency is always done and need not be in the list.
Conver=N
Sets the convergence calculations to 10^{N} on the energy and 10^{(N2)} on the wavefunction. The default is N=4 for single points and N=6 for gradients.
ENERGY RANGE OPTIONS
An energy range can be specified for CIS and TD excitation energies using the following options to CIS, TD and TDA.
GOccSt=N
Generate initial guesses using only active occupied orbitals N and higher.
GOccEnd=N
Generate initial guesses: if N>0, use only the first N active occupied orbitals; if N<0, do not use the highest N occupieds.
GDEMin=N
Generate guesses having estimated excitation energies ≥ N/1000 eV.
DEMin=N
Converge only states having excitation energy ≥ N/1000 eV; if N=2, read threshold from input; if N<2, set the threshold to N/1000 Hartrees.
IFact=N
Specify factor by which the number of states updated during initial iterations is increased.
WhenReduce=M
Reduce to the desired number of states after iteration M.
The default for IFact is Max(4,g) where g is the order of the Abelian point group. The default for WhenReduce is 1 for TD and 2 for TDA and CIS. Larger values may be needed if there are many states in the range of interest.
AVAILABILITY
Energies and gradients using HartreeFock or a DFT method.
RELATED KEYWORDS
CIS, ZIndo, Output
EXAMPLE
Here is the key part of the output from a TD excited states calculation:
Excitation energies and oscillator strengths:
Excited State 1: SingletA2 4.0147 eV 308.83 nm f=0.0000 <S**2>=0.000
8 > 9 0.70701
This state for optimization and/or secondorder correction.
Copying the excited state density for this state as the 1particle RhoCI density.
Excited State 2: SingletB1 9.1612 eV 135.34 nm f=0.0017 <S**2>=0.000
6 > 9 0.70617
Excited State 3: SingletB2 9.5662 eV 129.61 nm f=0.1563 <S**2>=0.000
8 > 10 0.70616
The results on each state are summarized, including the spin and spatial symmetry, the excitation energy, the oscillator strength, the S^{2}, and (on the second line for each state) the largest coefficients in the CI expansion.
The ECD results appear slightly earlier in the output as follows:
1/2[<0rb>*<brxdel0> + (<0rxdelb>*<br0>)*]
Rotatory Strengths (R) in cgs (10**40 ergesucm/Gauss)
state XX YY ZZ R(length) R(au)
1 0.0000 0.0000 0.0000 0.0000 0.0000
2 0.0000 0.0000 0.0000 0.0000 0.0000
3 0.0000 0.0000 0.0000 0.0000 0.0000
1/2[<0delb>*<br0> + (<0rb>*<bdel0>)*] (Au)
state X Y Z Dip. S. Osc.(frdel)
1 0.0000 0.0000 0.0000 0.0000 0.0000
2 0.0050 0.0000 0.0000 0.0050 0.0033
3 0.0000 0.2099 0.0000 0.2099 0.1399
Last update: 31 May 2013
