welcome: please sign in

Revision 10 as of 2018-05-24 18:28:01

Clear message
location: TD-DFT with SOC: open-shell systems

TD-DFT with SOC: open-shell systems

Contents

  1. TD-DFT with SOC: open-shell systems
    1. Open-shell systems with SOC: HI+

Examples from the paper "Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin–orbit splittings of light and heavy elements"", see http://pubs.rsc.org/-/content/articlehtml/2017/cp/c6cp07588f

Open-shell systems with SOC: HI+

$COMPASS
Title
 HI+ ground state fine structure splittings
Basis
 ano-rcc-vtzp
uncontracted
Geometry
 H  0.0   0.0 0.0
 I 1.62 0.0 0.0
END geometry
$END

$XUANYUAN
scalar
heff
3
soint
hsoc
2
$END

# Ground state for Si=1/2 with a sigma^+ excited state as reference to avoid symmetry breaking, the true ground state 2Pi will be computed in TDDFT as 'excited state' via de-excitations
$SCF
ROKS
DFT
BHHLYP
charge
1
spin
2
Alpha
13       2       6       6
Beta
12       2       6       6
THRESHCONV
1.d-12 1.d-10
$END

$TRAINT
utddft
orbi
hforb
$END

# Excited state calculations for Sf=Si=1/2 via spin-adapted TD-DFT (X-TDA)
$TDDFT
IMETHOD
 2
ISF
 0
ITDA
 1
IDIAG
 1
istore
 1
iexit
 200
crit_e
1.d-8
icorrect
1
itest
1
itrans
1
$END

# Excited state calculations for Sf=Si+1=3/2 via spin-flip TD-DFT (flip-up TDA with ALDA0 kernel)
$TDDFT
IMETHOD
 2
ISF
 1
ITDA
 1
IDIAG
 1
istore
 2
iexit
 200
crit_e
1.d-8
ialda
 2
$END

# TD-DFT SOC-SI step
$TDDFT
isoc
2
nfiles
2
ifgs
1
$END

 *** List of SOC-SI results ***

  No.      ExEnergies            Dominant Excitations         Esf        dE      Eex(eV)     (cm^-1)

    1      -3.9295 eV    50.0%  Spin: |So,1>    1-th   B2   -3.5859   -0.3436    0.0000         0.00
    2      -3.9295 eV    50.0%  Spin: |So,2>    1-th   B2   -3.5859   -0.3436    0.0000         0.00
    3      -3.3110 eV    49.4%  Spin: |So,2>    1-th   B1   -3.5859    0.2750    0.6186      4989.14
    4      -3.3110 eV    49.4%  Spin: |So,1>    1-th   B1   -3.5859    0.2750    0.6186      4989.14
    5      -0.1435 eV    97.9%  Spin: |Gs,2>    0-th   A1    0.0000   -0.1435    3.7860     30535.97
    6      -0.1435 eV    97.9%  Spin: |Gs,1>    0-th   A1    0.0000   -0.1435    3.7860     30535.97
    7       4.6229 eV    49.9%  Spin: |So,2>    2-th   B2    4.9533   -0.3304    8.5524     68980.01
    8       4.6229 eV    49.9%  Spin: |So,1>    2-th   B2    4.9533   -0.3304    8.5524     68980.01
    9       5.2520 eV    49.1%  Spin: |So,1>    2-th   B1    4.9533    0.2987    9.1816     74054.18
   10       5.2520 eV    49.1%  Spin: |So,2>    2-th   B1    4.9533    0.2987    9.1816     74054.18
...

Small difference with the reported value 5005.5cm^-1 is due to the no. of scalar states included or (ialda != 2 used here). For heavy elements like Iodine, convergence with respect to this number should be tested carefully.