|
Size: 4063
Comment:
|
Size: 4217
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: |
=== ifgs === |
|
| Line 89: | Line 60: |
| No. of TD-DFT calculations to be loaded. | |
| Line 95: | Line 67: |
| =3, just print SOC matrix elements between two spin-free states | =3, just print SOC matrix elements between two spin-free states (without diagonalization Hsoc). |
| Line 97: | Line 69: |
| === irsf === Unused. === irso === === insf === === inso === |
=== ifgs === =0, default for not including ground state (GS) in SOC treatment; =1, include GS. |
| Line 108: | Line 74: |
| Line 132: | Line 99: |
| isf | isf |
| Line 152: | Line 119: |
| In this input, 7 means seven of SOC matrices will be calculate. 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. | 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 154: | 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. === irsf === Unused. === irso === |
|
| Line 159: | Line 135: |
| Line 174: | 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 176: | Line 153: |
| === insf === === inso === |
|
| Line 177: | Line 156: |
| Line 179: | 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 184: | 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 187: | Line 164: |
| == Stability analysis == |
== Stability analysis == |
| Line 190: | Line 166: |
| Line 192: | Line 167: |
| Line 194: | Line 168: |
| Line 196: | Line 169: |
| Line 198: | Line 170: |
| Line 200: | Line 171: |
tddft: time-dependent density functional theory
Contents
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
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