welcome: please sign in

Revision 82 as of 2016-03-30 11:25:56

Clear message
location: tddft

tddft: time-dependent density functional theory

Contents

  1. tddft: time-dependent density functional theory
    1. Quick guides
    2. General keywords
      1. imethod
      2. isf
      3. itda
      4. idrpa
      5. ispa
      6. ialda
      7. thrdab
      8. itest
      9. icorrect
      10. itrans
      11. iro
      12. icv
      13. ioo
      14. iksf
      15. iact
      16. elw
      17. eup
      18. idiag
      19. ndiag
      20. aokxc
    3. Convergence threshold
      1. crit_e
      2. crit_vec
    4. States specification
      1. iexit
      2. nexit
    5. Save eigenvectors
      1. istore
      2. lefteig
    6. output eigenvector control
      1. nprt
      2. cthrd
    7. TD-DFT/SOC and Property evaluation
      1. nfiles
      2. isoc
      3. ifgs
      4. imatsoc
      5. imatrsf
      6. imatrso
      7. imatnso
      8. imatnsf
      9. idiag
      10. iact
    8. Stability analysis
      1. isab
      2. isave
    9. memory control
      1. memjkop
    10. Others
      1. isgn
      2. ivo
    11. Modified Davidson algorithm
      1. Eneshift
    12. AO-TDDFT and AO-TDA

Time dependent DFT/HF calculation. Support Full TDDFT, TDA and RPA.

Quick guides

The following examples give the minimal inputs for starting TD-DFT calculations.

1. Closed-shell Systems : R-TD-DFT

2. Open-shell Systems : U-TD-DFT and spin-adapted TD-DFT for spin-conserving excitations

3. Open-shell Systems : Spin-flip TD-DFT for spin-flip excitations

4. TD-DFT with SOC

5. TD-DFT with SOC: Kramers pairs

6. Excitation analyze based on molecular fragments

7. TD-DFT with spin-flip calculations

General keywords

imethod

  imethod 1, R-TDDFT, start from RKS
  imethod 2, U-TDDFT, start from UKS or ROKS
  imethod 3, X-TDDFT, start from ROKS

isf

  Spin flip TDDFT. 
  isf 0, do not flip
  isf 1, spin flip up
  isf -1, spin flip down

itda

  itda 0, TDDFT, do not use TDA
  itda 1, TDA

idrpa

ispa

ialda

thrdab

itest

icorrect

itrans

iro

icv

ioo

iksf

iact

iact = 1: define active space based on energy [elw,eup]

elw

lower bound in eV (not in au!).

eup

upper bound in eV.

idiag

ndiag

aokxc

Convergence threshold

crit_e

crit_vec

States specification

iexit

The number of calculated excited states of each irreducible representation for a specific point group

nexit

Same as above

Save eigenvectors

istore

Integer: specify the file no. to store TDDFT information

lefteig

By default, in TD-DFT the left eigenvector X-Y is also stored.

output eigenvector control

nprt

cthrd

TD-DFT/SOC and Property evaluation

nfiles

No. of TD-DFT calculations to be loaded.

isoc

=1, Only work for closed-shell case (NOT recommended!)

=2, General SOC state interaction

=3, just print SOC matrix elements between two spin-free states (without diagonalization Hsoc).

ifgs

=0, default for not including ground state (GS) in SOC treatment; =1, include GS.

imatsoc

Define SOC matrices need to be calculated. Input format looks like 

...
#SCF calculation for the ground state S0. It is a singlet.
$scf
spin
 0
...
$end

#First TDDFT, singlets S0-S9.
$tddft
imethod
 1
isf
 0
iext
 10
....
$end

#Second TDDFT, triplet T1-T10
$tddft
imethod
 1
isf
 1
iexit
 10
$end

$tddft
....
imatsoc
  7
0 0 0 2 1 1
0 0 0 2 1 2
1 1 1 2 1 1
1 1 1 2 1 2
1 1 2 2 1 1
1 1 2 2 1 2
2 1 1 2 1 1
2 1 1 2 1 2
$end

In this input, 7 means seven of SOC matrices will be calculate (If the number <0, then ALL possible HSOC mat will be printed !). Here, it is very tricky to specify states:

imatrsf

Transition dipole between Spin-free states. The input is similar to imatsoc (but currently selected printing is not implemented). Simply use -1 to print all of them.

imatrso

Define transition dipole moment need to be printed between two SOC-included states. Input format looks like(notice we omit other input in TDDFT module) 

$TDDFT
...
imatrso
5
1 1
1 2
1 3
1 4
1 5
...
$END

Then, "imatrso" is specified to define transition dipole moments need to be printed. The number "5" require transition dipoles between 5-pairs of states to be print. The following 5 lines define which pairs will be printed. Here, we require transition dipoles between the first state and five states are printed.

imatnso

imatnsf

idiag

By default, idiag=0 uses full diagonalization (preferred for small model space).

If idiag=1, then TD-DFT/SOC can use Davidson's algorithm also, along with a specification for the no. of states by iexit.

iact

=1, allows to use active space specification for the projected active-orbital SOC Hamiltonian (P*HSOC*P), eup can be specified in (eV) to give a cut off to define active physically interested excited states.

Stability analysis

isab

isave

memory control

memjkop

Others

isgn

ivo

Modified Davidson algorithm

Eneshift

Specify an energy window. States with excitation energies close to input value will be calculate. The energy unit is eV. 

$TDDFT
Eneshift
 9.0
...
$End

AO-TDDFT and AO-TDA

AO-TDDFT supports R-TDDFT, U-TDDFT, R-TDDFT-SF+1. The possible combinations are
     imethod=1, itda=0, isf=0
     imethod=1, itda=0, isf=1
     imethod=2, itda=0, isf=0
AO-TDA  supports R-TDA,U-TDA, R-TDA-SF+1, R-TDA-SF-1, U-TDA-SF3. The possible combinations are
     imethod=1, itda=1, isf=0
     imethod=1, itda=1, isf=1
     imethod=2, itda=1, isf=0
     imethod=2, itda=1, isf=-1
     imethod=2, itda=1, isf=3