![Page 1: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/1.jpg)
235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana
COMPUTATIONAL SPECTROSCOPY
Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry University of Calgary,Alberta, Canada T2N 1N4
Time-dependent density functional theory as a practical tool in the study of MCD and CD spectra of transition metal complexes. Implementations and applications.
New Orleans National Meeting
Organizers: H. B. Schlegel and Krishnan RaghavachariPreciding:H. B. SchlegelWednesday April 9 10:20 -11:00 am Morial Convention Center -- Rm. 342,
![Page 2: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/2.jpg)
ADF• Solves Kohn-Sham equations• Properties
– NMR, EFG, EPR, Raman, IR, UV/Vis, NLO, CD, MCD…– Potential energy surfaces (transition states, geometry
optimization)• Environment effects
– QM/MM, COSMO• Relativistic effects
– Scalar relativistic effects, spin-orbit coupling– Transition and heavy metal compounds
• Uses Slater functions
![Page 3: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/3.jpg)
hv
Cl
C
C
C
C
C
Si Zr
C
C
C
C
Cl
C
Inorganic SpectroscopyInorganic Spectroscopy
![Page 4: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/4.jpg)
Basic Time Dependent Density Functionl TheoryBasic Time Dependent Density Functionl Theory
Basic Equation :Basic Equation :
€
Aia, jb = (εa − εi)0δ ijδab +∂F ia
∂Pjb
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟0
€
Bia,bj =∂F ia
∂Pbj
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟0
Definition of A and B Matrices :Definition of A and B Matrices :
M.E.CasidaM.E.Casida
Gross,E.K.; Kohn W.Gross,E.K.; Kohn W.
€
ΩF (λ ) =Wλ2F (λ )
€
Ω=−S−1/2(A + B)S−1/2
€
S−1/ 2 = (A − B)1/ 2
Where :Where :
![Page 5: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/5.jpg)
Basic Time Dependent Density Functionl TheoryBasic Time Dependent Density Functionl Theory
€
Wλ = ΔE o,λTransition Energy :Transition Energy :
Basic Equation :Basic Equation :M.E.CasidaM.E.Casida
Gross,E.K.; Kohn W.Gross,E.K.; Kohn W.
€
ΩF (λ ) = Wλ2F (λ )
Electric Transition Dipole Moment :Electric Transition Dipole Moment :
€
Aα ˆ M Jλ =1
WJ
μ iaFia(Jλ ) (εa − εi)
ia
∑
€
μ ia = − ir r a
€
Jλ ˆ L Aα = WJ liaFia( Jλ ) 1
(ε a −ε i )ia
∑
Magnetic Transition Dipole Moment :Magnetic Transition Dipole Moment :
€
l jb = −iμB jr r ×
r ∇ b
![Page 6: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/6.jpg)
A
C
B
A
B
C
Absorption Spectra and TD-DFTAbsorption Spectra and TD-DFT
Transition Energy :Transition Energy :
€
fλ =2
3μ iaFia
(Jλ ) (εa − εi)ia
∑ ⎡
⎣ ⎢
⎤
⎦ ⎥⋅ μ jbF jb
(Jλ ) (εb − ε j )jb
∑∫ ⎡
⎣ ⎢ ⎢
⎤
⎦ ⎥ ⎥
€
Wλ = ΔE0,λ
![Page 7: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/7.jpg)
hv
Cl
C
C
C
C
C
Si Zr
C
C
C
C
Cl
C
Inorganic SpectroscopyInorganic Spectroscopy
N
MN N
N
H
![Page 8: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/8.jpg)
![Page 9: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/9.jpg)
Why MCD and MOR ? Why MCD and MOR ?
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
More information about each excited stateMore information about each excited state
In absorption spectroscopy onlypositive (often overlapping) bandsIn absorption spectroscopy onlypositive (often overlapping) bands
![Page 10: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/10.jpg)
Why MCD ? Why MCD ?
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
In MCD bands of different shapesMore information about each excited stateIn MCD bands of different shapesMore information about each excited state
![Page 11: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/11.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
€ €
Origin of MCD ? Origin of MCD ?
€
AJ = γoω(NAαg − N Jλ j )
Nαλgj∑ Aαg ˆ M Jλ j
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟× ρ Aαg.Jλ j (ω)
€
AJ = γoω(NAαg − N Jλ j )
Nαλgj∑ Aαg ˆ M Jλ j
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟× ρ Aαg.Jλ j (ω)
Electric dipole operator:Electric dipole operator:
€
ˆ M = ˆ m ii∑ = − (xii
∑ r e xi
+ yi
r e yi
+ yi
r e yi
)
€
ˆ M = ˆ m ii∑ = − (xii
∑ r e xi
+ yi
r e yi
+ yi
r e yi
)
Absorbance in dipole approximation.Absorbance in dipole approximation.
Aα
Jλ
![Page 12: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/12.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
Electric dipole operatorFor circular polarizedLight:
Electric dipole operatorFor circular polarizedLight:
€
ˆ M ± = ˆ m ±,ii∑
€
ˆ M ± = ˆ m ±,ii∑€
ΔAJ
ω= γo
(NAαg − N Jλj )
Nαλgj∑ Aαg ˆ M − Jλ j
2− Aαg ˆ M + Jλ j
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟× ρ Aαg.Jλ j (ω)
Difference in absorbance of left and right circular polarized lightDifference in absorbance of left and right circular polarized light
Circular Polarized LightCircular Polarized Light
€
ˆ m − =1
2(x
r e x − iy
r e y )
€
ˆ m + =1
2(x
r e x + iy
r e y )
Origin of MCD ? Origin of MCD ?
![Page 13: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/13.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
€
ΔAJ'
ω=
A−,J'
ω−
A+,J'
ω
γo
(NAαg − N Jλj )
Nαλgj∑ Aαg ˆ M − Jλ j
2− Aαg ˆ M + Jλ j
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟.ρ Aαg.Jλ j (ω)
⎛
⎝ ⎜
⎞
⎠ ⎟
'
The difference in absorption of left and right circularly polarized light in the presence of a magnetic field as a function of photon energy
The difference in absorption of left and right circularly polarized light in the presence of a magnetic field as a function of photon energy
Origin of MCD ? Origin of MCD ?
€
ρAαg.Jλj (ω) ≈ fJ (ω −ωJλ ) =1
πWJ
e−
ωJλ −ω
WJ
⎛
⎝ ⎜
⎞
⎠ ⎟
2
![Page 14: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/14.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
€
ΔAJ'
hω=
γo
3 γ
3
∑ (NAα − N Jλ )
Nαλ∑ Aα ˆ M −
γ Jλ2
− Aα ˆ M +γ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟ ∂
∂Bγ
fJ (ω −ωJλ )( )oB
+γo
3 γ
3
∑ (NAα − N Jλ )
Nαλ∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
fJ (ω −ωJ )B
+γo
3
∂
∂B γ
3
∑ (NAα − N J )
N
⎛
⎝ ⎜
⎞
⎠ ⎟o
α∑ Aα ˆ M −
γ Jλ2
− Aα ˆ M +γ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟ fJ (ω −ωJ )B
€
ΔAJ'
hω=
γo
3 γ
3
∑ (NAα − N Jλ )
Nαλ∑ Aα ˆ M −
γ Jλ2
− Aα ˆ M +γ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟ ∂
∂Bγ
fJ (ω −ωJλ )( )oB
+γo
3 γ
3
∑ (NAα − N Jλ )
Nαλ∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
fJ (ω −ωJ )B
+γo
3
∂
∂B γ
3
∑ (NAα − N J )
N
⎛
⎝ ⎜
⎞
⎠ ⎟o
α∑ Aα ˆ M −
γ Jλ2
− Aα ˆ M +γ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟ fJ (ω −ωJ )B
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B +γoC J fJ (ω −ωJ )B
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B +γoC J fJ (ω −ωJ )B
Origin of MCD ? Origin of MCD ?
![Page 15: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/15.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
The MCD disprsion The MCD disprsion
€
ΔA '
ω= γoB -A J
∂B J fJ (ω −ωJ )B
∂ω+(B J +
C J
kT)B J fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
ΔA '
ω= γoB -A J
∂B J fJ (ω −ωJ )B
∂ω+(B J +
C J
kT)B J fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
MCD dispersion over all bands MCD dispersion over all bands
P.J.Stephens. Ph.D. Thesis 1964
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
€
γo
3 γ
3
∑ (NAα − N Jλ )
Nλ∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
fJ (ω −ωJ )B
€
ΔA '
ω= γoB -A J
∂fJ (ω −ωJ )B
∂ω+(B J +
C J
kT) fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
ΔA '
ω= γoB -A J
∂fJ (ω −ωJ )B
∂ω+(B J +
C J
kT) fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
![Page 16: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/16.jpg)
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
The MCD disprsion The MCD disprsion
P.J.Stephens. Ph.D. Thesis 1964
A
B
C(T)
A
![Page 17: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/17.jpg)
Origin of B-term in Magnetic Circular DichroismOrigin of B-term in Magnetic Circular Dichroism
The MCD dispersionThe MCD dispersion
MCD dispersion due to B-term MCD dispersion due to B-term P.J.Stephens. Ph.D. Thesis 1964
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
€
ΔAJ'
ω= −γoA J
∂fJ (ω −ωJ )
∂ωB +γ0B J fJ (ω −ωJ )B
+γoC J fJ (ω −ωJ )B
€
γ0B J fJ (ω −ωJ )B =
γo
3 γ
3
∑ (NAα − N Jλ )
Nλ∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
fJ (ω −ωJ )B
ThusThus
€
B J =1
3 γ
3
∑ (NAα − N Jλ )
Nλ∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
=1
3 γ
3
∑λ
∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
![Page 18: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/18.jpg)
Origin of B-TermOrigin of B-Term
The B term The B term
O
X+iaY
Y-iaX
M- M+
B>0
ΔA-A+
A-
B>0
O
X
Y
M- M+
B=0
-A+
A-ΔA
B=0
€
ΔA '
ω= γoB -A J
∂fJ (ω −ωJ )B
∂ω+(B J +
C J
kT) fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
ΔA '
ω= γoB -A J
∂fJ (ω −ωJ )B
∂ω+(B J +
C J
kT) fJ (ω −ωJ )B
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
B J =1
3 γ
3
∑λ
∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
![Page 19: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/19.jpg)
Expression for the B-TermExpression for the B-Term
The B term The B term
€
B J =1
3 γ
3
∑λ
∑ ∂
∂Bγ
Aα ˆ M −γ Jλ
2− Aα ˆ M +
γ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
Or by using the identity Or by using the identity
€
t
∑ A ˆ M t J2
− A ˆ M t J2
= i ε rst
r,s ,t
∑ A ˆ M r J J ˆ M s A = i ε rst
r ,s,t
∑ α rs (ωL )
€
Here ε rst is the three - dimensional Levi - Civita symbol
€
Here ε rst is the three - dimensional Levi - Civita symbol
We thus haveWe thus have
€
B J =i
3ε stu
∂α st (ω)
∂Bu
⎛
⎝ ⎜
⎞
⎠ ⎟ω=ωL( )
![Page 20: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/20.jpg)
The Calculation of the B-termThe Calculation of the B-term
The B term : practical calculations The B term : practical calculations
We have:We have:
TD-DFT calculationsTD-DFT calculations
€
€
B J =i
3ε stu
∂α st (ω)
∂Bu
⎛
⎝ ⎜
⎞
⎠ ⎟
s,t ,u
∑ω=ωL
Where:Where:
€
α st (ω)ω= ω L= M s[XL (ω) + YL (ω)][M t (XL (ω) + YL (ω)]
€
ωC 0
0 −C
⎛
⎝ ⎜
⎞
⎠ ⎟X
Y
⎛
⎝ ⎜
⎞
⎠ ⎟=
A B
B* A*
⎛
⎝ ⎜
⎞
⎠ ⎟X
Y
⎛
⎝ ⎜
⎞
⎠ ⎟
Early work:
J.Michl, J.Am.Chem.Soc. 100,6801 (1978)
Early work:
J.Michl, J.Am.Chem.Soc. 100,6801 (1978)
![Page 21: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/21.jpg)
The B term : practical calculations The B term : practical calculations
€
where A, B C are defined by
€
where A, B C are defined by
€
Aaiσ ,bjτ =δστ δabδij
εbτ −ε jτ
nbτ − n jτ
− Kaiσ ,bjτ
€
Kaiσ , jbτ = dr dr 'ϕ aσ* (r)∫∫ ϕ iσ (r)
1
r − r'ϕ bτ (r' )ϕ jτ
* (r' )
+ dr dr 'ϕ aσ* (r)∫∫ fXC (r,r',ω)ϕ bτ (r' )ϕ jτ
* (r' )
€
Baiσ ,bjτ = Kaiσ , jbτ
€
Caiσ ,bjτ =δστ δabδij
1
nbτ − n jτ
The Calculation of the B-termThe Calculation of the B-term
![Page 22: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/22.jpg)
The B term : practical calculations The B term : practical calculations
We have:We have:
TD-DFT calculationsTD-DFT calculations
Solve:
€
€
B J = −i
3ε rst
∂α st (ω)
∂Br
⎛
⎝ ⎜
⎞
⎠ ⎟
r,s ,t
∑ω=ωL
Where:Where:
€
ω(0) C 0
0 −C
⎛
⎝ ⎜
⎞
⎠ ⎟X (0)
Y (0)
⎛
⎝ ⎜
⎞
⎠ ⎟=
A(0) B(0)
B(0) A(0)
⎛
⎝ ⎜
⎞
⎠ ⎟X (0)
Y (0)
⎛
⎝ ⎜
⎞
⎠ ⎟
€
BAJ = −2i
3ε stu
stu
∑ M s (XJ(1) t +YJ
(1)u )M t (XJ(0) +YJ
(0) )
The Calculation of the B-termThe Calculation of the B-term
![Page 23: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/23.jpg)
The B term : practical calculations The B term : practical calculations
By differentiation ofBy differentiation of
€
ω(0) C 0
0 −C
⎛
⎝ ⎜
⎞
⎠ ⎟X (0)
Y (0)
⎛
⎝ ⎜
⎞
⎠ ⎟=
A(0) B(0)
B(0) A(0)
⎛
⎝ ⎜
⎞
⎠ ⎟X (0)
Y (0)
⎛
⎝ ⎜
⎞
⎠ ⎟
€
BAJ = −2i
3εαβγ
αβγ
∑ M β (XJ(1)α +YJ
(1)α )M λ (XJ(0) +YJ
(0) )
€
Implementation - X(1), Y(1)( )
€
The equation that we use for evaluating (X(1), Y(1) ) is
€
A(0) B(0)
B(0)* A(0)*
⎛
⎝ ⎜
⎞
⎠ ⎟−ωI
(0) −I 0
0 I
⎛
⎝ ⎜
⎞
⎠ ⎟
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟X I
(1)
YI(1)
⎛
⎝ ⎜
⎞
⎠ ⎟=
ωI(1) −I 0
0 I
⎛
⎝ ⎜
⎞
⎠ ⎟−
A(1) B(1)
B(1)* A(1)*
⎛
⎝ ⎜
⎞
⎠ ⎟
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟X I
(0)
YI(0)
⎛
⎝ ⎜
⎞
⎠ ⎟
The Calculation of the B-termThe Calculation of the B-term
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 24: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/24.jpg)
€
Evaluation of - X(1),Y(1)( )
€
Evaluation of - X(1),Y(1)( )
€
Introducing the unitary transformation U =1 1
1 -1
⎛
⎝ ⎜
⎞
⎠ ⎟
€
UA(0) B(0)
B(0)* A(0)*
⎛
⎝ ⎜
⎞
⎠ ⎟−ωI
(0) −I 0
0 I
⎛
⎝ ⎜
⎞
⎠ ⎟
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟U
+UX I
(1)
YI(1)
⎛
⎝ ⎜
⎞
⎠ ⎟=U ωI
(1) −I 0
0 I
⎛
⎝ ⎜
⎞
⎠ ⎟−
A(1) B(1)
B(1)* A(1)*
⎛
⎝ ⎜
⎞
⎠ ⎟
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟U
+UX I
(0)
YI(0)
⎛
⎝ ⎜
⎞
⎠ ⎟
€
ω I(0)I −Ω[ ]Z I
(1) = −ωI(0)S1/2(A(1) + B(1) )S−1/2FI
(0)
−ωI(0)S−1/2(A(1) − B(1) )S1/2FI
(0)
Here:Here:
€
Z I(1) = ωI
(0) S1/2(X I(1) +YI
(1) )
AffordsAffords
The Calculation of the B-termThe Calculation of the B-term
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 25: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/25.jpg)
The B term : practical calculations The B term : practical calculations
€
ω I(0)I −Ω[ ]Z I
(1) = −ωI(0)S1/2(A(1) + B(1) )S−1/2FI
(0)
−ωI(0)S−1/2(A(1) − B(1) )S1/2FI
(0)
€
An Expression for Kai,bj(1)
€
An Expression for Kai,bj(1)
€
We need φp(1). A well known expresson exists that is particularly simple because we have an
imaginary perturbation
€ €
€
φp(1) = Uqp
(1)φq(0) Uqp
(1) =H pq
(1)
εq(0) −ε p
(0)q≠ p
∑
€
φp(1) = Uqp
(1)φq(0) Uqp
(1) =H pq
(1)
εq(0) −ε p
(0)q≠ p
∑
€
Where H(1) is the Hamiltonian describing the perurbation
Thus
€
Kai,bj(1) = U pa
(1)*K pi,bj(0)
p≠a
∑ + U pi(1)*Kap,bj
(0) + U pb(1)*Kai,pj
(0)
p≠b
∑p≠i
∑ + U pb(1)*Kap,bp
(0)
p≠ j
∑
€
Kai,bj(1) = U pa
(1)*K pi,bj(0)
p≠a
∑ + U pi(1)*Kap,bj
(0) + U pb(1)*Kai,pj
(0)
p≠b
∑p≠i
∑ + U pb(1)*Kap,bp
(0)
p≠ j
∑
The Calculation of the B-termThe Calculation of the B-term
Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 26: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/26.jpg)
The B term : Direct method The B term : Direct method
We must solve We must solve
€
AX = b€
ω I(0)I −Ω[ ]Z I
(1) = −ωI(0)S1/2(A(1) + B(1) )S−1/2FI
(0)
−ωI(0)S−1/2(A(1) − B(1) )S1/2FI
(0)
€
Our equation has the form
€
Our equation has the form
€
With A a known matrix, b a known vector and X the unknown vector to be determined. This
equation can be solved easily if we have A−1. There are two problems however
€
With A a known matrix, b a known vector and X the unknown vector to be determined. This
equation can be solved easily if we have A−1. There are two problems however
€
(a) The matrix A =ωI(0)I −Ω. This matrix has no inverse because ωI
(0)I is an eigenvalue of Ω
€
(b) The matrix A is extremely large and we don' t want to try and invert it directly.
€
To avoid this problem we :
€
To avoid this problem we :
€
(i) Solve the equations iteratively by expanding the solution in a Krylov
subspace(the space b,Ab, A2b,...Aib in the ith iteration)
€
(ii) Project out from the Krylov supspaces any contribution from FI(0)
The Calculation of the B-term by Direct MethodThe Calculation of the B-term by Direct Method
Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 27: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/27.jpg)
The B term : Direct Method The B term : Direct Method
We must solve We must solve
€
Ax = b
€
(i) Can be used in conjunction with an unperturbed
TDDFT calculation that yields only a few solutions F(0).
€
(ii)Degree of convergence is known
ProsPros
ConsCons
€
(i) The iterative procedure is often slowly convergent.
We are attempting to improve convergence by adding
the unperturbed TDDFT solutions FJ(0), J ≠ I to Krylov subspace
€
BAJ = −2i
3εαβγ
αβγ
∑ M β S−1/2
ωJ(0)
(Z J(1)α )M λ (XJ
(0) +YJ(0) )
€
ω I(0)I −Ω[ ]Z I
(1) = −ωI(0)S1/2(A(1) + B(1) )S−1/2FI
(0)
−ωI(0)S−1/2(A(1) − B(1) )S1/2FI
(0)
The Calculation of the B-term by Direct MethodThe Calculation of the B-term by Direct Method
Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 28: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/28.jpg)
The B term : Sum Over State The B term : Sum Over State
€
ω I(0)I −Ω[ ]Z I
(1) = −ωI(0)S1/2(A(1) + B(1) )S−1/2FI
(0)
−ωI(0)S−1/2(A(1) − B(1) )S1/2FI
(0)
€
Z (1) by Sum - Over - State
€
Z (1) by Sum - Over - State
€
Z I(1) = CJIFJ
(0)
J≠I
∑
€
Z I(1) = CJIFJ
(0)
J≠I
∑
€
Substitute into first order equation and multiply by FJ(0) from left affords
€
FJ(0)+ ωI
(0)I −Ω[ ]( (CJIFJ(0)
J≠I
∑ ) = −ωI(0)FJ
(0)+S1/2 (A(1) + B(1) )S−1/2FI(0)
−ωI(0)FJ
(0)+S−1/2(A(1) − B(1) )S1/2FI(0)
OrOr
€
CJI =−ωI
(0)(FJ(0)+S1/2(A(1) + B(1) )S−1/2FI
(0) + FJ(0)+S−1/2(A(1) − B(1) )S1/2FI
(0)
ωI(0) −ωJ
(0)
Writing Z(1) in terms of the complete set F(0) affordsWriting Z(1) in terms of the complete set F(0) affords
The Calculation of the B-term by Sum-over-State MethodThe Calculation of the B-term by Sum-over-State Method
Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 29: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/29.jpg)
The B term : Sum Over State The B term : Sum Over State
€
Z (1) by Sum - Over - State : Z I(1) = CJI
J≠I
∑ FJ(0)
€
CJI =−ωI
(0)(FJ(0)+S1/2(A(1) + B(1) )S−1/2FI
(0)
ωI(0) −ωJ
(0)+
+FJ(0)+S−1/2(A(1) − B(1) )S1/2FI
(0)
ωI(0) −ωJ
(0)
€
BAJ = −2i
3εαβγ
αβγ
∑ M β S−1/2
ωJ(0)
(Z J(1)α )M λ (XJ
(0) +YJ(0) )
€
BAJ = −2i
3εαβγ
αβγ
∑ M β S−1/2
ωJ(0)
(Z J(1)α )M λ (XJ
(0) +YJ(0) )
ProsPros
€
Interpretation easy in terms of contributions
from different excited statesConsCons
€
May need to calculate many FJ(0) in unperturbed
TDDFT and convergence of summation is unknown
The Calculation of the B-term by Sum-over-State MethodThe Calculation of the B-term by Sum-over-State Method
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 30: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/30.jpg)
Convergence of SOS-method for Ethylene Convergence of SOS-method for Ethylene
Comparison of Sum-over-State and Direct Method for B-termsComparison of Sum-over-State and Direct Method for B-terms
€
π → π *
€
π → 3s
€
π → 3s
€
π → π *
Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 31: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/31.jpg)
Comparison of Direct Method for B-terms with ExperimentComparison of Direct Method for B-terms with Experiment
Exp: H.-P.Klein, R.T. Oakley, J.Michl Inorg.Chem. 25,3194 (1986)
Exp: H.-P.Klein, R.T. Oakley, J.Michl Inorg.Chem. 25,3194 (1986) Seth+Ziegler JCP,2008,in pressSeth+Ziegler JCP,2008,in press
![Page 32: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/32.jpg)
Comparison of Direct Method for B-terms with Experiment and other MethodsComparison of Direct Method for B-terms with Experiment and other Methods
Exp: H.-P.Klein, R.T. Oakley, J.Michl Inorg.Chem. 25,3194 (1986)
T.Kjœrgaard, B.Jansik, P.Jørgensen,S.Coriani, J.Michl, J.Phys.Chem. A 111,11278 (2007))
![Page 33: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/33.jpg)
TD-DFT calculations of B-term. TD-DFT calculations of B-term.
Furan
Thiophene
Selenophen
Tellurophen
S
Se
Te
OW. Hieringer, S. J. A. van Gisbergen, and E. J. BaerendsJ. Phys. Chem. A 2002, 106, 10380
X
1a2
1b1
2b1
11A1 --> 21A1
11A1 --> 11B2
X
X
X 1b2 11A1 --> 11B1
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 34: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/34.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
W. Hieringer, S. J. A. van Gisbergen, and E. J. BaerendsJ. Phys. Chem. A 2002, 106, 10380
aVarsanyi, G.; Nyulazi, L.; Vezpremi, T.; Narisawa, T. J. Chem.Soc., Perkin Trans. 2 1982, 761.
![Page 35: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/35.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
Norden, B.; Hansson, R.; Pedersen, P. B.; Thulstrup, E. W. Chem.Phys. 1978, 33, 355.
![Page 36: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/36.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
Norden, B.; Hansson, R.; Pedersen, P. B.; Thulstrup, E. W. Chem.Phys. 1978, 33, 355.
![Page 37: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/37.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
Furan
O
6.0 6.2
0.20
-5.05 0.0 3.37
.13
€
1a2 → 3b1
€
2b1 → 3b1
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 38: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/38.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
Thiophene
S
5.5 5.7
450 6-477
.04
5.9
.13
€
1a2 → 3b1
€
2b1 → 3b1
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 39: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/39.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
Selonophene
Se
5.1 5.3
0.22
59.1 -3 -101
.07
5.5
€
1a2 → 3b1
€
2b1 → 3b1
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 40: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/40.jpg)
TD-DFT calculations of B-term. Sum-over-state formulationTD-DFT calculations of B-term. Sum-over-state formulation
TellurophenX
1a2
1b1
2b1
11A1 --> 21A1
11A1 --> 11B2
X
X
X 1b2 11A1 --> 11B1
Te
4.4 4.8
0.64
-5.1 -28.012.8
5.24.4 4.8
0.64
-5.1 -28.012.8
5.2
Seth+Ziegler JCP,2008,134108Seth+Ziegler JCP,2008,134108
![Page 41: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/41.jpg)
MOR and MCDMOR and MCD
Consider a planar polarized light traveling a distance l through a media of randomly oriented molecules along the direction ofa constant magnetic field with strength B.
Consider a planar polarized light traveling a distance l through a media of randomly oriented molecules along the direction ofa constant magnetic field with strength B.
€
€
α =V (ω)Bl
Here V( ) is called the Verdet constantHere V( ) is called the Verdet constant
BE E
α
l
For such a system the plane of polarization will rotate by an angle given byFor such a system the plane of polarization will rotate by an angle given by
€
Vsos (ω) ≡ −μ 0cN
3
hω2B J
WJ2 − hω2
J
∑
€
Vsos (ω) ≡ −μ 0cN
3
hω2B J
WJ2 − hω2
J
∑
aM. Krykunov, A. Banerjee, T. Ziegler,J. Autschbach J. Chem. Phys. 2005, 122, 075105,
A.Banerjee,J.Autschbach,T.Ziegler Int.J.Quant.Chem.2006,101,572
![Page 42: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/42.jpg)
MOR and MCDMOR and MCD
€
€
Vsos (ω) ≡ −μ 0cN
3
hω2B J
ωJ2 −ω2
J
∑
€
Vsos (ω) ≡ −μ 0cN
3
hω2B J
ωJ2 −ω2
J
∑€
VRe s (ω)€
VRe s (ω) =ωμ ocN
12Im[
∂
∂Bx
{α yz (ω) −α zy(ω}]ω≠ωJ
+ωμ ocN
12Im[
∂
∂By
{α zx (ω) −α xz (ω}]ω≠ωJ
+ωμ ocN
12[Im
∂
∂Bz
{α xy(ω) −α yz (ω)}]ω≠ωJ
ωμ ocN
12Imε stu
∂
∂Bu
(α st (ω)) Bu =0( )
€
VRe s (ω) =ωμ ocN
12Im[
∂
∂Bx
{α yz (ω) −α zy(ω}]ω≠ωJ
+ωμ ocN
12Im[
∂
∂By
{α zx (ω) −α xz (ω}]ω≠ωJ
+ωμ ocN
12[Im
∂
∂Bz
{α xy(ω) −α yz (ω)}]ω≠ωJ
ωμ ocN
12Imε stu
∂
∂Bu
(α st (ω)) Bu =0( )
€
∂αst (ω)
∂Bu
⎡
⎣ ⎢
⎤
⎦ ⎥ω≠ωJ
=
−Im [∂
∂Bu
Δρ (t )(ω,r r )]ω≠ωJ∫ xsd
r r
€
∂αst (ω)
∂Bu
⎡
⎣ ⎢
⎤
⎦ ⎥ω≠ωJ
=
−Im [∂
∂Bu
Δρ (t )(ω,r r )]ω≠ωJ∫ xsd
r r
€
B J
€
∂ααβ (ω)
∂Bγ
⎡
⎣ ⎢
⎤
⎦ ⎥
ω=ωJ
=
−Im [∂
∂Bλ
Δρ (β )(ω,r r )]ω=ωJ∫ xα d
r r
€
∂ααβ (ω)
∂Bγ
⎡
⎣ ⎢
⎤
⎦ ⎥
ω=ωJ
=
−Im [∂
∂Bλ
Δρ (β )(ω,r r )]ω=ωJ∫ xα d
r r
€
B J =
−i
3ε rst
∂α st (ω)
∂Br
⎛
⎝ ⎜
⎞
⎠ ⎟
r,s ,t
∑ω=ωL
M. Krykunov, A. Banerjee, T. Ziegler,J. Autschbach J. Chem. Phys. 2005, 122, 075105,
![Page 43: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/43.jpg)
MOR and MCDMOR and MCD
€
Vsos (ω) = −γω2B J
ωJ2 −(ω)2
J
∑
€
Vsos (ω) = −γω2B J
ωJ2 −(ω)2
J
∑The expressionThe expression
Allows us to calculate the MOR parameter V( ) from the MCDparameters BJ after summing over all statesAllows us to calculate the MOR parameter V( ) from the MCDparameters BJ after summing over all states
aM. Krykunov, A. Banerjee, T. Ziegler,J. Autschbach J. Chem. Phys. 2005, 122, 075105,
![Page 44: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/44.jpg)
MOR and MCDMOR and MCD
€
V (ω) = −γω2B J
ωJ2 −(ω)2
J
∑
€
V (ω) = −γω2B J
ωJ2 −(ω)2
J
∑The expressionThe expression
Vres(
€
Vdamp(ω)€
The expression for V(ω) diverges for ω = ωJ
€
The expression for V(ω) diverges for ω = ωJ
We need a TD-DFT formulation in which damping includedWe need a TD-DFT formulation in which damping included
![Page 45: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/45.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation without damping TD-DFT formulation without damping
€
We solve the equation
€
We solve the equation
€
ˆ h ks (v r )+V ext (
v r , t) − i
∂
∂t
⎡ ⎣ ⎢
⎤ ⎦ ⎥φk (r)× exp[−iε kt] = 0
To obtain the solutionTo obtain the solution
€
φk' (
v r , t) = C j (t)φ j
j≠i
∑ (v r )× exp[−iε jt]
€
Δρ ( y)(ω,r r ) = [X(ω)ai +Y (ω)ai ]φaφi
a
vir
∑i
occ
∑
From which we obtain density change in frequency domainFrom which we obtain density change in frequency domain
With:With:
€
(X(ω)+Y (ω)) = 2S−1/2[ω2 −Ω]−1S−1/2V (ω)
![Page 46: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/46.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
We solve the equation
€
We solve the equation
€
ˆ h ks (v r )+V ext (
v r , t) − i
∂
∂t
⎡ ⎣ ⎢
⎤ ⎦ ⎥φk (r)× exp[−iε kt]exp[−Γt] = 0
To obtain finite lifetime solutionsTo obtain finite lifetime solutions
€
φk' (
v r , t) = C j (t)φ j
j≠i
∑ (v r )× exp[−iε jt]exp[−λ t]
€
Δρ ( y)(ω,r r ) = [X(ω)ai +Y (ω)ai ]φaφi
a
vir
∑i
occ
∑
From which we obtain density change in frequency domainFrom which we obtain density change in frequency domain
With:
€
(X(ω)+Y (ω)) = 2S−1/2[(ω + iγ )2 −Ω]S−1/2V (ω)L.Jensen; J.Autchbach; G.C.Schatz J.Chem.Phys.2005,122,224115
![Page 47: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/47.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
Vresdm (ω) =Vres
R,dm (ω)+ iVresI ,dm
HereHere
€
VresR,dm (ω) =Vsos
R,dm (ω) ≡J
∑ 2ω2 (ωJ2 −ω2 )B J
(ωJ2 −ω2 )2 + 4ω2γ 2
€
VresI ,dm (ω) =Vsos
I ,dm (ω) ≡J
∑ 4ω3γB J
(ωJ2 −ω2 )2 + 4ω2γ 2andand
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,127,244107
![Page 48: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/48.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
VresR,dm (ω) =Vsos
R,dm (ω) ≡ γ0
J
∑ 2ω2(ωJ2 −ω2 )B J
(ωJ2 −ω2 )2 + 4ω2γ 2
oror
€
VsosR,dm (ω) =Vsos
udm (ω) fd (ω)
€
VsosR,dm (ω) =Vsos
udm (ω) fd (ω)
€
fd (ω) =
2(ωJ2 −ω2 )2
(ωJ2 −ω2 )2 + 4ω2γ 2
HereHere
€
Vsosudm (ω)
€
VsosR,dm (ω)
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,127,244107
![Page 49: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/49.jpg)
MOR and MCDMOR and MCD
TD-DFT formulation with damping TD-DFT formulation with damping
€
VresI ,dm (ω) =Vsos
I ,dm (ω) ≡J
∑ 4ω3γB J
(ωJ2 −ω2 )2 + 4ω2γ 2
€
VresI ,dm (ω) /ω = γ0 ωB J
J
∑ fJ,B (ω) = Δε MCD (ω)oror
€
fJ,B (ω) =4ω γ
(ωJ2 −ω2 )2 + 4ω2γ 2
€
fJ,B (ω)
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,127,244107
![Page 50: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/50.jpg)
MOR and MCDMOR and MCD
€
VresI ,dm (ω) /ω = γ0 ωB J
J
n
∑ fJ,B (ω) = Δε MCD (ω)
€
We can obtain B J (j =1,n) from a least square fit of
€
D =i=1
m
∑ VresI ,dm (ωi ) /ωi −γ0 ωiB J
J
n
∑ fJ,B (ωi ) ⎛
⎝ ⎜
⎞
⎠ ⎟
2
For m>n
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,127,244107
![Page 51: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/51.jpg)
MOR and MCDMOR and MCD
€
Δε
Furan Thiophene
Selenophen Tellurophene
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,127,244107
![Page 52: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/52.jpg)
MOR and MCDMOR and MCD
Transitions Energy eV ImV( ) SOS Furan 1b1‡ 2b1 6.00 -7.80 -5.05 1a2‡ 2b1 6.33 3.08 3.37 Thiophene 1b1‡ 2b1 5.56 420 450 1a2‡ 2b1 5.60 -453 -477 Selenophene 1b1‡ 2b1 5.23 48.5 59.1 1a2‡ 2b1 5.21 -98.9 -101 Tellurophene 1b1‡ 2b1 4.42 7.75 12.8 1a2‡ 2b1 5.06 -30.6 -27.8
Calculated B-term constants from
X
1a2
1b1
2b1
11A1 --> 21A1
11A1 --> 11B2
X
X
X 1b2 11A1 --> 11B1
€
VresI ,dm (ω)
![Page 53: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/53.jpg)
A-term of MCDA-term of MCD
€
ΔAJ'
hω=
γo
1
3γ
∑ (NA − N Jλ )
NJλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅
∂
∂Bγ
f (ω −ωJλ )( )0B
+γo
1
3γ
∑ (NA − N Jλ )
N
∂
∂BγJλ
∑ A ˆ M γ− Jλ
2− A ˆ M γ
+ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
⋅ f (ω −ωJλ )( )B
+γo
1
3γ
∑ ∂
∂Bγ
(NA − N Jλ )
N
⎛
⎝ ⎜
⎞
⎠ ⎟0
Jλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅ f (ω −ωJλ )( )B
€
ΔAJ'
hω=
γo
1
3γ
∑ (NA − N Jλ )
NJλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅
∂
∂Bγ
f (ω −ωJλ )( )0B
+γo
1
3γ
∑ (NA − N Jλ )
N
∂
∂BγJλ
∑ A ˆ M γ− Jλ
2− A ˆ M γ
+ Jλ2 ⎛
⎝ ⎜ ⎞
⎠ ⎟o
⋅ f (ω −ωJλ )( )B
+γo
1
3γ
∑ ∂
∂Bγ
(NA − N Jλ )
N
⎛
⎝ ⎜
⎞
⎠ ⎟0
Jλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅ f (ω −ωJλ )( )B
€
ΔAJ'
ω= −γoA J
∂f (ω −ωJλ )
∂ω+γ0B J f (ω −ωJλ )B +γoC J f (ω −ωJλ )B
€
ΔAJ'
ω= −γoA J
∂f (ω −ωJλ )
∂ω+γ0B J f (ω −ωJλ )B +γoC J f (ω −ωJλ )B
Origin of A-term Origin of A-term
M.Seth,T.Ziegler,J.Chem.Phys. 2004,120,10943
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
![Page 54: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/54.jpg)
The A-term of Magnetic Circular Dichroism (MCD) SpectroscopyThe A-term of Magnetic Circular Dichroism (MCD) Spectroscopy
The A term The A term
€
γo -A J
∂f (ω −ωJ )
∂ω
⎡ ⎣ ⎢
⎤ ⎦ ⎥J
∑ B
= γo
1
3γ
∑ (NA − N Jλ )
NJλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅
∂
∂Bγ
f (ω −ωJλ )( )0B
= γo
1
3γ
∑ (NA − N Jλ )
NJλ∑ A ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟⋅
∂ωJλ
∂Bγ
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟0
∂f (ω −ωJ )
∂ωB
ThusThus
€
A J = −λ
∑ 1
3γ
∑ (NA − N Jλ )
NA ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟0
⋅∂ωJλ
∂Bγ
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟
M.Seth,T.Ziegler,J.Chem.Phys. 2004,120,10943
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
![Page 55: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/55.jpg)
The A-term of Magnetic Circular Dichroism (MCD) SpectroscopyThe A-term of Magnetic Circular Dichroism (MCD) Spectroscopy
The A term The A term
€
A J = −λ
∑ 1
3γ
∑ (NA − N Jλ )
NA ˆ M γ
− Jλ2
− A ˆ M γ+ Jλ
2 ⎛ ⎝ ⎜ ⎞
⎠ ⎟0
⋅∂ωJλ
∂Bγ
⎛
⎝ ⎜ ⎜
⎞
⎠ ⎟ ⎟
HereHere
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
ThusThus
![Page 56: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/56.jpg)
The A term The A term
€
A−'
hω−
A+'
hω=
ΔA'
hω= γβB -A J
∂ρA .J (hω)
∂(hω)+ (B J +
C J
kT)ρA .J (hω)
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
A−'
hω−
A+'
hω=
ΔA'
hω= γβB -A J
∂ρA .J (hω)
∂(hω)+ (B J +
C J
kT)ρA .J (hω)
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
-A+
A-ΔA
B=0
RCP LCP
1S
1P
B=0
-A+
A-
ΔA-
B>0
RCP LCP
O
-1
1
0
B>0
€
A ˆ M Jλ =r
∑ M rS−1/ 2FLλ(0)v
e r
€
A ˆ M Jλ =r
∑ M rS−1/ 2FLλ(0)v
e r
The A-term of Magnetic Circular Dichroism (MCD) SpectroscopyThe A-term of Magnetic Circular Dichroism (MCD) Spectroscopy
€
ˆ L = ˆ l ij=1
n
∑ = −ir r i
j=1
n
∑ ×r
∇ i
![Page 57: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/57.jpg)
Applications:A/D
Se42+ Te4
2+
Fe(CN)64-
Ni(CN)42- C6Cl6
C6H3Br3
Oh
D4hD4h
D6h
D3h
Exp: 0.72Calc: 0.63
Exp: 0.60Calc: 0.55Exp: 0.40 Calc: 0.48
Exp:-0.66Calc:-0.72Exp:-0.50Calc:-0.80
M.Seth,T.Ziegler,J.Chem.Phys. 2004,120,10943
![Page 58: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/58.jpg)
Positive A-term
Negative A-termNegative
B-termPositive B-term
Different MCD-terms
Absorption band
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
3t2
2e
t1
2t2
Metal
Ligand
3t2
2e
t1
2t2
Metal
Ligand
![Page 59: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/59.jpg)
MCD-terms for Oxyanions MCD-terms for Oxyanions
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
![Page 60: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/60.jpg)
Exp.Theor
MCD-terms for Thioanions MCD-terms for Thioanions
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. Submitted
![Page 61: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/61.jpg)
N
MN N
N
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
MCD spectra of Porphyrins containing Mg,Ni and ZnMCD spectra of Porphyrins containing Mg,Ni and Zn
![Page 62: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/62.jpg)
2e1.g
2a2.u
1a2.u
1a1.u
1b2.u1e1.g
ZnP
1b1.g
Orbital level diagram for ZnPOrbital level diagram for ZnP
2a2u
1a1u
1b2u
E.J. Baerends , G. Ricciardi , A. Rosa , S.J.A van GisbergenJ.Phys.Chem. A2001,105,3311
E.J. Baerends , G. Ricciardi , A. Rosa , S.J.A van GisbergenCoord.Chem.Rev. 2002,230,5
2eg2
2eg1
![Page 63: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/63.jpg)
Exc. Energ. (eV) Complex Symmetry
exp. calc. Composition % h
f Assig n.
2a 2u -> 2e g 52.10 1E u
2.03 c, 2.21 d,
2.23 e , 2.18 f
2.28 1a 1u -> 2e g 46.63
0.001 Q
1b 2u -> 2e g 68.44
1a 1u -> 2e 1g 17.54 2E u 2.95 c, 3.09 d,
3.18 e
, 3.13 f 3.25
2a 2u -> 2e g 10.05 0.496
1b 2u -> 2e g 29.88 2a 2u -> 2e g 29.31 1a 1u -> 2e g 27.13
ZnP
3E u
3.32
1a 2u -> 2e g 10.30
0.943
Bg
Experimental Spectrum for ZnPExperimental Spectrum for ZnP
E.J. Baerends , G. Ricciardi , A. Rosa , S.J.A van GisbergenCoord.Chem.Rev. 2002,230,5
E.J. Baerends , G. Ricciardi , A. Rosa , S.J.A van GisbergenCoord.Chem.Rev. 2002,230,5
1Eu C1Φ(2a2u → 2eg)+C2Φ(1a1u --> 2eg)Gouterman State
2Eu
3Eu
C2Φ(2a2u → 2eg)-C1Φ(1a1u --> 2eg)Conjugated Gouterman State
Φ(1b2u → 2eg)
1A1g Ground State
1Eu C1Φ(2a2u → 2eg)+C2Φ(1a1u --> 2eg)Gouterman State
2Eu
3Eu
C2Φ(2a2u → 2eg)-C1Φ(1a1u --> 2eg)Conjugated Gouterman State
Φ(1b2u → 2eg)
1A1g Ground State
2e1.g
2a2.u
1a2.u
1a1.u
1b2.u1e1.g
ZnP
1b1.g
![Page 64: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/64.jpg)
21Eu 31Eu
5 10-2
Experimental Spectrum for ZnPExperimental Spectrum for ZnP
L.Edwards,D.H.Dolphin,M.Goutermn J.Mol.Spectrosc 35(1970)90
E.J. Baerends , G. Ricciardi , A. Rosa , S.J.A van GisbergenCoord.Chem.Rev. 2002,230,5
1Eu C1Φ(2a2u → 2eg)+C2Φ(1a1u --> 2eg)Gouterman State
2Eu
3Eu
C2Φ(2a2u → 2eg)+C1Φ(1a1u --> 2eg)Conjugated Gouterman State
Φ(1b2u → 2eg)
1A1g Ground State
1Eu C1Φ(2a2u → 2eg)+C2Φ(1a1u --> 2eg)Gouterman State
2Eu
3Eu
C2Φ(2a2u → 2eg)+C1Φ(1a1u --> 2eg)Conjugated Gouterman State
Φ(1b2u → 2eg)
1A1g Ground State
€
D(1Eu ) = C1 2a2u y 2egy +C2 1a1u y 2egx[ ]2
1
22.92 −
1
23.25
⎡ ⎣ ⎢
⎤ ⎦ ⎥
2
= 2.27x10−2
€
D(3Eu ) = C1 2a2u y 2egy +C2 1a1u y 2egx[ ]2
1
22.92 +
1
23.25
⎡ ⎣ ⎢
⎤ ⎦ ⎥
2
= 9.51
![Page 65: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/65.jpg)
Simulated Spectrum for ZnP with A-term onlySimulated Spectrum for ZnP with A-term only
ZnP Exp
A-only
Comp l ex Sy m m e try h
A
h
A/D
1E u 0.05 5.49
2E u - 3.37 - 1.62
ZnP
3E u - 0.57 - 0.15
1Eu
2Eu+3Eu
Q
Q
S
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem. Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem. Inorg. Chem. 2007,46, 9111-9125.
![Page 66: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/66.jpg)
Influence of ring distortion on MCD spectrum of ZnPInfluence of ring distortion on MCD spectrum of ZnP
€
B(nB1) = Im−2
3
nB1ˆ L z nB2 A1
ˆ M x nB1 nB2ˆ M y A1
ΔWn
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
€
B(nB2 ) = Im2
3
nB1ˆ L z nB2 A1
ˆ M x nB1 nB2ˆ M y A1
ΔWn
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
€
B (nB1)
€
B (nB2)
€
B (nB1) + B (nB2)
€
A (nEu)
N
N
N
N
M
ω
nB1
nB2
C2vD4h
nEu
N
N
N
N
M
![Page 67: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/67.jpg)
Influence of ring distortion on MCD spectrum of ZnPInfluence of ring distortion on MCD spectrum of ZnP
2.00 2.50 3.00 3.50E(eV)
ZnPx10
Normalized Intensity-0.5
0.50.0
Dist C2V
2.00 2.50 3.00 3.50E(eV)
ZnPx10
Normalized Intensity-0.5
0.50.0
D4h
€
B(nB1) = Im−2
3
nB1ˆ L z nB2 A1
ˆ M x nB1 nB2ˆ M y A1
ΔWn
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
€
B(nB2 ) = Im2
3
nB1ˆ L z nB2 A1
ˆ M x nB1 nB2ˆ M y A1
ΔWn
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
N
N
N
N
M
ω
nB1
nB2
C2vD4h
nEu
N
N
N
N
M
![Page 68: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/68.jpg)
Simulated Spectrum for ZnP with B-term onlySimulated Spectrum for ZnP with B-term only
Exp.
€
B(nEu ) = Im−4
3 p≠n
∑nEux
ˆ L z pEuy A1gˆ M x nEux pEuy
ˆ M y A1g
W ( pE1uy )− W (nE1uy )
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
1Eu 2Eu3Eu
B-terms
![Page 69: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/69.jpg)
Simulated Spectrum forZnP with A+B-term onlySimulated Spectrum forZnP with A+B-term only
E (eV) E (eV)
x 100
ZnP
-0.50
0.00
0.50
1.00
Normalized Intensity
2.00 2.50 3.00 3.50
Exp.
€
B (3Eu ) =
Im−4
3 p≠n
∑3Eux
ˆ L z 2Euy A1gˆ M x 3Eux 2Euy
ˆ M y A1g
W (2E1uy ) −W (3E1uy )
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
€
B (2Eu ) =
Im−4
3 p≠n
∑2Eux
ˆ L z 3Euy A1gˆ M x 2Eux 3Euy
ˆ M y A1g
W (3E1uy ) −W (2E1uy )
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
€
B (3Eu ) =
Im−4
3 p≠n
∑2Eux
ˆ L z 3Euy A1gˆ M x 3Eux 2Euy
ˆ M y A1g
W (2E1uy ) −W (3E1uy )
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪
= −B (2Eu )
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
![Page 70: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/70.jpg)
Simulated Spectrum for MgP and NiP with A+B-termSimulated Spectrum for MgP and NiP with A+B-term
2eg 2eg 2eg
2a2u2a2u2a2u
1a1u
1a2u 1a2u 1a2u
1a1u1a1u
1b2u1b2u1b2u
dxy
1eg 1eg
1egdxz, dyz
1eu
dxz, dyz
dx2-y2
dz2
MgP NiP ZnP-7.00
-9.50
-12.00
E(eV)
1b1g
2eg 2eg 2eg
2a2u2a2u2a2u
1a1u
1a2u 1a2u 1a2u
1a1u1a1u
1b2u1b2u1b2u
dxy
1eg 1eg
1egdxz, dyz
1eu
dxz, dyz
dx2-y2
dz2
MgP NiP ZnP-7.00
-9.50
-12.00
E(eV)
1b1g
2.0 2.5 3.0 3.5
E(eV)
(a) MgP
0.0
0.5
0.5
x100
2.0 2.5 3.0 3.5
E(eV)
(a) MgP
0.0
0.5
0.5
x100
2.0 2.5 3.0 3.5
E(eV)
0.0
0.5
0.5
(b) NiP
x100
2.0 2.5 3.0 3.5
E(eV)
0.0
0.5
0.5
(b) NiP
x100
1Eu
2Eu3Eu
![Page 71: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/71.jpg)
Substituted Porphyrins Substituted Porphyrins
N
M
m
β
N N
N
N
MN N
N
MTPP
N
MN N
N
MOEPtetraphenylporphyrin octaethylporphyrin
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
![Page 72: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/72.jpg)
Excited States for Substituted Porphyrins Excited States for Substituted Porphyrins
2.00 2.50 3.00 3.50E(eV)
Normalized Intensity-0.5
0.5
0.0
NiTPP
N
NiN N
N
€
A (1Eu )
€
B (2Eu )
€
B (3Eu )
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
![Page 73: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/73.jpg)
Excited States for Substituted Porphyrins Excited States for Substituted Porphyrins
2.00 2.50 3.00 3.50E(eV)
Normalized Intensity-0.5
0.5
0.0
ZnTPP
x10
N
ZnN N
N
€
A (1Eu )
€
B (2Eu )
€
B (3Eu )
€
A (1Eu )
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
![Page 74: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/74.jpg)
Tetraazaporphyrins and PhthalocyaninesTetraazaporphyrins and Phthalocyanines
N
M
m
β
N N
N
![Page 75: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/75.jpg)
Tetraazaporphyrins and PhthalocyaninesTetraazaporphyrins and Phthalocyanines
N
M
m
β
N N
N
N
M
N N
N N
NN N
MTAPtetraazaporphyrin
![Page 76: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/76.jpg)
Tetraazaporphyrins and PhthalocyaninesTetraazaporphyrins and Phthalocyanines
N
M
m
β
N N
N
N
M
N N
N N
NN N
MTAPtetraazaporphyrin
N
M
N N
N N
NN N
MPcphthalocyanine
Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125. Alejandro Gonzales, Mike Seth, Tom Ziegler Inorg.Chem.
Inorg. Chem. 2007,46, 9111-9125.
![Page 77: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/77.jpg)
Magnetic Circular Dichroism (MCD) SpectroscopyMagnetic Circular Dichroism (MCD) Spectroscopy
The C term The C term
€
A−'
hω−
A+'
hω=
ΔA'
hω= γβB -A J
∂ρA .J (hω)
∂(hω)+ (B J +
C J
kT)ρA .J (hω)
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
€
A−'
hω−
A+'
hω=
ΔA'
hω= γβB -A J
∂ρA .J (hω)
∂(hω)+ (B J +
C J
kT)ρA .J (hω)
⎡
⎣ ⎢ ⎤
⎦ ⎥J∑
1P
1S
M- M+
B=0
-A+
A-ΔA
B=0
1P+
1S
M- M+
B>0
1P-
ΔA-A+
A-
B>0
If
€
NP+− NP+
N tot
≈EP+
− EP+
3kT
€
EP+− EP+
<< kT
€
C = −i
3 AAα '
αa 'λ
∑ ˆ L Aα ⋅ Aα ˆ M Jλ × Jλ ˆ M Aα ' ⎛ ⎝ ⎜
⎞ ⎠ ⎟
€
C = −i
3 AAα '
αa 'λ
∑ ˆ L Aα ⋅ Aα ˆ M Jλ × Jλ ˆ M Aα ' ⎛ ⎝ ⎜
⎞ ⎠ ⎟
![Page 78: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/78.jpg)
Electron configuration t1u6t2u
6t1u6t2g
5
Seth,Ziegler,Autschbach,Ziegler JCP, 2005,09412
![Page 79: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/79.jpg)
Application: Fe(CN)63-Application: Fe(CN)63-
Electron configuration t1u6t2u
6t1u6t2g
5
Excitations are ligand-metal charge transfer. C term of a transition to a T1u state is positive andto a T2u state is negative.
Transition Exp. Calc.
1 1.21/0.61 0.86
2 -0.68 -0.86
3 0.56 0.86
Seth,Ziegler,Autschbach,Ziegler JCP, 2005,09412
![Page 80: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/80.jpg)
More Applications
RuCl63- [Fe(CN)5SCN]3-
MnPc
Exp. Calc.
0.58 0.84
-0.60 -0.84
Exp Calc
7.5 7.3
6.9 7.3
-6.9 -7.3
6.3 7.3
-3.1 -7.3
2.2 7.3
Exp. Calc.
0.03 0.90
0.23 0.90
Seth,Ziegler,Autschbach,Ziegler JCP, 2005,09412
![Page 81: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/81.jpg)
Limitations of Traditional TD-DFTLimitations of Traditional TD-DFT
ix iy
What do we do with adegenerate ground statethat can not be representedby a single Slater determinant ?
What do we do with adegenerate ground statethat can not be representedby a single Slater determinant ?
Degenerate Ground StateDegenerate Ground State
a
ix iy
a
ix iy
What are thefundamentalequations ?
What are thefundamentalequations ?
How do we calculateexcitationenergies
How do we calculateexcitationenergies
![Page 82: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/82.jpg)
Transformed Reference with an Intermediate ConfigurationKohn Sham (TRICKS) TDDFT
Solution:Solution:
TRICKS of the Trade: Calculating the Excitation Energies of Molecules with Degenerate Ground States using TD-DFTTRICKS of the Trade: Calculating the Excitation Energies of Molecules with Degenerate Ground States using TD-DFT
Idea:Idea:Avoid problems with a degenerate ground state by taking an excited state that is nondegenerate as the (Transformed) Reference Intermediate Configuration.
A. I. Krylov, Acc. Chem. Res. 2006, 39, 83-91
![Page 83: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/83.jpg)
Example 2:d1 transition metal complexes of Td symmetry,d-d transition
Example 2:d1 transition metal complexes of Td symmetry,d-d transition
Application of the TRIC methodApplication of the TRIC method
€
VCl4
![Page 84: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/84.jpg)
Result 2:d1 transition metal complexes of Td symmetry,d-d transition
Result 2:d1 transition metal complexes of Td symmetry,d-d transition
Application of the TRIC methodApplication of the TRIC method
€
VCl4
![Page 85: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/85.jpg)
Example 3:d1 transition metal complexes of Td symmetry,charge transfer
Example 3:d1 transition metal complexes of Td symmetry,charge transfer
Application of the TRIC methodApplication of the TRIC method
€
VCl4
![Page 86: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/86.jpg)
Result 3:d1 transition metal complexes of Td symmetry,charge transfer
Result 3:d1 transition metal complexes of Td symmetry,charge transfer
Application of the TRIC methodApplication of the TRIC method
€
VCl4
![Page 87: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/87.jpg)
Conclusion• Method for calculating the MCD A term (and dipole strength D)
within TD-DFT is outlined. Procedure for calculating C/D more straightforward.
• Implemented into the Amsterdam Density Functional Theory (ADF) program
• Applications to a range of small molecules
• Further information can be found in M. Seth, T Ziegler, A Banerjee, J. Autschbach, S.J.A. van Gisbergen E. J. Baerends, J. Chem. Phys. 120,10942, 2004 and M. Seth, T. Ziegler, J. Autschbach, J. Chem. Phys. accepted for publication.
![Page 88: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/88.jpg)
TD-DFT/MCD
Dr. Mike SethDr.Jochen Autschbach
Alejandro Gonzalez Peralta
Dr. Mykhaylo Krykunov
Fan Wang
Hristina Zhekova
![Page 89: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/89.jpg)
PRF
Mitsui
![Page 90: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/90.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation without damping TD-DFT formulation without damping
€
We solve the equation
€
We solve the equation
€
ˆ h ks (v r )+V ext (
v r , t) − i
∂
∂t
⎡ ⎣ ⎢
⎤ ⎦ ⎥φk (r)× exp[−iε kt] = 0
To obtain the solutionTo obtain the solution
€
φk' (
v r , t) = C j (t)φ j
j≠i
∑ (v r )× exp[−iε jt]
€
Δρ ( y)(ω,r r ) = [X(ω)ai +Y (ω)ai ]φaφi
a
vir
∑i
occ
∑
From which we obtain density change in frequency domainFrom which we obtain density change in frequency domain
With:With:
€
(X(ω)+Y (ω)) = 2S−1/2[ω2 −Ω]−1S−1/2V (ω)
![Page 91: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/91.jpg)
MOR and MCDMOR and MCD
€
Vsos (ω) = −γhω2B J
WJ2 −(hω)2
J
∑
€
Vsos (ω) = −γhω2B J
WJ2 −(hω)2
J
∑The expressionThe expression
Allows us to calculate the MOR parameter V( ) from the MCDparameters BJ after summing over all statesAllows us to calculate the MOR parameter V( ) from the MCDparameters BJ after summing over all states
aM. Krykunov, A. Banerjee, T. Ziegler,J. Autschbach J. Chem. Phys. 2005, 122, 075105,
![Page 92: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/92.jpg)
MOR and MCDMOR and MCD
€
V (ω) = −γhω 2B J
WJ2 − (hω)2
J
∑
€
V (ω) = −γhω 2B J
WJ2 − (hω)2
J
∑The expressionThe expression
Vres(
€
Vdamp(ω)
€
The expression for V(ω) diverges for hω = WJ
€
The expression for V(ω) diverges for hω = WJ
We need a TD-DFT formulation in which damping includedWe need a TD-DFT formulation in which damping included
![Page 93: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/93.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
We solve the equation
€
We solve the equation
€
ˆ h ks (v r )+V ext (
v r , t) − i
∂
∂t
⎡ ⎣ ⎢
⎤ ⎦ ⎥φk (r)× exp[−iε kt]exp[−Γt] = 0
To obtain finite lifetime solutionsTo obtain finite lifetime solutions
€
φk' (
v r , t) = C j (t)φ j
j≠i
∑ (v r )× exp[−iε jt]exp[−λ t]
€
Δρ ( y)(ω,r r ) = [X(ω)ai +Y (ω)ai ]φaφi
a
vir
∑i
occ
∑
From which we obtain density change in frequency domainFrom which we obtain density change in frequency domain
With:
€
(X(ω)+Y (ω)) = 2S−1/2[(ω + iγ )2 −Ω]S−1/2V (ω)L.Jensen; J.Autchbach; G.C.Schatz J.Chem.Phys.2005,122,224115
![Page 94: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/94.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
Vresdm (ω) =Vres
R,dm (ω)+ iVresI ,dm
HereHere
€
VresR,dm (ω) =Vsos
R,dm (ω) ≡J
∑ 2ω2 (ωJ2 −ω2 )B J
(ωJ2 −ω2 )2 + 4ω2γ 2
€
VresI ,dm (ω) =Vsos
I ,dm (ω) ≡J
∑ 4ω3γB J
(ωJ2 −ω2 )2 + 4ω2γ 2andand
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,submitted
![Page 95: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/95.jpg)
MOR and MCD`MOR and MCD`
TD-DFT formulation with damping TD-DFT formulation with damping
€
VresR,dm (ω) =Vsos
R,dm (ω) ≡ γ0
J
∑ 2ω2(ωJ2 −ω2 )B J
(ωJ2 −ω2 )2 + 4ω2γ 2
oror
€
VsosR,dm (ω) =Vsos
udm (ω) fd (ω)
€
VsosR,dm (ω) =Vsos
udm (ω) fd (ω)
€
fd (ω) =
2(ωJ2 −ω2 )2
(ωJ2 −ω2 )2 + 4ω2γ 2
HereHere
€
Vsosudm (ω)
€
VsosR,dm (ω)
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,submitted
![Page 96: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/96.jpg)
MOR and MCDMOR and MCD
TD-DFT formulation with damping TD-DFT formulation with damping
€
VresI ,dm (ω) =Vsos
I ,dm (ω) ≡J
∑ 4ω3γB J
(ωJ2 −ω2 )2 + 4ω2γ 2
€
VresI ,dm (ω) /ω = γ0 ωB J
J
∑ fJ,B (ω) = Δε MCD (ω)oror
€
fJ,B (ω) =4ω γ
(ωJ2 −ω2 )2 + 4ω2γ 2
€
fJ,B (ω)
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,submitted
![Page 97: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/97.jpg)
MOR and MCDMOR and MCD
M.Krykunov,M.Seth,T.Ziegler,J.Autschbach J.Chem.Phys. 2007,submitted
€
VresI ,dm (ω) /ω = γ0 ωB J
J
n
∑ fJ,B (ω) = Δε MCD (ω)
€
We can obtain B J (j =1,n) from a least square fit of
€
D =i=1
m
∑ VresI ,dm (ωi ) /ωi −γ0 ωiB J
J
n
∑ fJ,B (ωi ) ⎛
⎝ ⎜
⎞
⎠ ⎟
2
For m>n
![Page 98: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/98.jpg)
MCD spectra of Porphyrins containing Mg,Ni and ZnMCD spectra of Porphyrins containing Mg,Ni and Zn
N
MN N
N
N
MN N
N
MTPP
tetraphenylporphyrin
MOEP
tetraphenylporphyrin
N
M
N N
N N
NN N
MTAP
N
M
N N
N N
NN N
tetraazaporphyrin MPcphthalocyanine
porphyrin
MP
N
M
m
β
N N
N
![Page 99: 235nd ACS National Meeting April 5-10, 2008 New Orleans, Louisiana COMPUTATIONAL SPECTROSCOPY Mike Seth, Mykhaylo Krykunov, Tom Ziegler Department of Chemistry](https://reader030.vdocument.in/reader030/viewer/2022032522/56649d6c5503460f94a4c0f6/html5/thumbnails/99.jpg)
Example 2:Double ExcitationsExample 2:Double Excitations
Application of the TRIC methodApplication of the TRIC method
Seth,M., Ziegler,T. , J. Chem. Phys., 2005,123, 144105,
Seth,M. ; Ziegler,T. J. Chem. Phys. 2006, 124, 144105
Groundstate
Excited state
TRIC state