welcome: please sign in
location: Diff for "tddft"
Differences between revisions 58 and 59
Revision 58 as of 2014-11-28 08:30:44
Size: 4299
Editor: 162
Comment:
Revision 59 as of 2014-11-28 08:37:24
Size: 4217
Editor: 162
Comment:
Deletions are marked like this. Additions are marked like this.
Line 13: Line 13:
Line 17: Line 16:
1. [[ Closed-shell Systems : R-TD-DFT]] 1. [[Closed-shell Systems : R-TD-DFT]]
Line 19: Line 18:
2. [[ Open-shell Systems : U-TD-DFT and spin-adapted TD-DFT for spin-conserving excitations ]] 2. [[Open-shell Systems : U-TD-DFT and spin-adapted TD-DFT for spin-conserving excitations]]
Line 21: Line 20:
3. [[ Open-shell Systems : Spin-flip TD-DFT for spin-flip excitations ]] 3. [[Open-shell Systems : Spin-flip TD-DFT for spin-flip excitations]]
Line 23: Line 22:
4. [[ TD-DFT with SOC ]] 4. [[TD-DFT with SOC]]
Line 26: Line 25:
Line 28: Line 26:
Line 30: Line 27:
Line 32: Line 28:
Line 34: Line 29:
Line 36: Line 30:
Line 38: Line 31:
Line 40: Line 32:
Line 42: Line 33:
Line 44: Line 34:
Line 46: Line 35:
Line 48: Line 36:
Line 50: Line 37:
Line 52: Line 38:
Line 54: Line 39:
Line 56: Line 40:
Line 58: Line 41:
Line 60: Line 42:
Line 62: Line 43:
Line 64: Line 44:
Line 66: Line 45:
Line 69: Line 47:
Line 71: Line 48:
Line 81: Line 57:
Line 83: Line 58:
Line 85: Line 59:
Line 101: Line 74:
Line 125: Line 99:
isf  isf
Line 144: Line 118:
}}}
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:
Line 145: Line 121:
}}}  * The string "0 0 0" always treat as the ground state.
 * For other states, three numbers "n m n" represent "ith-tddft", "symmetry" and "ithstate" respectively. Therefore, the first matrix element "0 0 0 2 1 1" means SOC matrix of <S0|HSOC|T1>. The third matrix element "1 1 1 2 1 1" means SOC matrix <S1|HSOC|T1>. Here, the first "1" in bra state "1 1 1" means the state from first TDDFT calculation. The second and third "1" in the bra state "1 1 1" means this state has spatial symmetry "1" and is the first excited state.
Line 147: Line 124:
In this input, 7 means seven of SOC matrices will be calculate. '''If the number <0, then ALL possible HSOC mat will be calculated !'''
Line 149: Line 125:
Here, it is very tricky to specify states. First, "0 0 0" always treat as the ground state. Second: For other states, three numbers "n m n" represent "ith-tddft", "symmetry" and "ithstate" respectively. Therefore, the first matrix element "0 0 0 2 1 1" means SOC matrix of <S0|HSOC|T1>. The third matrix element "1 1 1 2 1 1" means SOC matrix <S1|HSOC|T1>. Here, the first "1" in bra state "1 1 1" means the state from first TDDFT calculation. The second and third "1" in the bra state "1 1 1" means this state has spatial symmetry "1" and is the first excited state.
Line 155: Line 130:
Line 161: Line 135:
Line 176: Line 151:
Here, "irso" is set to 1 to enable transition dipole moment calculation. Then, "imatros" 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 ground state and five low-lying states are printed.  Here, "irso" is set to 1 to enable transition dipole moment calculation. Then, "imatros" 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 ground state and five low-lying states are printed.
Line 179: Line 154:
Line 181: Line 155:
Line 183: Line 156:
Line 185: Line 157:
By default, idiag=0 uses full diagonalization (preferred for small model space).  By default, idiag=0 uses full diagonalization (preferred for small model space).
Line 190: Line 162:
=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.
=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.
Line 193: Line 164:
==  Stability analysis ==
== Stability analysis ==
Line 196: Line 166:
Line 198: Line 167:
Line 200: Line 168:
Line 202: Line 169:
Line 204: Line 170:
Line 206: Line 171:

tddft: time-dependent density functional theory

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

General keywords

imethod

isf

itda

idrpa

ispa

ialda

thrdab

itest

icorrect

itrans

iro

icv

ioo

iksf

iact

elw

eup

idiag

ndiag

aokxc

States specification

iext

next

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
iext
 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:

  • The string "0 0 0" always treat as the ground state.
  • For other states, three numbers "n m n" represent "ith-tddft", "symmetry" and "ithstate" respectively. Therefore, the first matrix element "0 0 0 2 1 1" means SOC matrix of <S0|HSOC|T1>. The third matrix element "1 1 1 2 1 1" means SOC matrix <S1|HSOC|T1>. Here, the first "1" in bra state "1 1 1" means the state from first TDDFT calculation. The second and third "1" in the bra state "1 1 1" means this state has spatial symmetry "1" and is the first excited state.

irsf

Unused.

irso

imatrsf

Unused.

imatrso

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

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

Here, "irso" is set to 1 to enable transition dipole moment calculation. Then, "imatros" 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 ground state and five low-lying states are printed.

insf

inso

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

tddft (last edited 2022-10-28 13:20:58 by bsuo)