Transport and optical response
- f molecular junctions
UCSD July 20-21, 2009
Michael Galperin University of California at San Diego
UCSD July 20-21, 2009 – p.1
Transport and optical response of molecular junctions UCSD July - - PowerPoint PPT Presentation
Transport and optical response of molecular junctions UCSD July 20-21, 2009 Michael Galperin University of California at San Diego UCSD July 20-21, 2009 p.1 Introduction UCSD July 20-21, 2009 p.2 Introduction Timescale BO Energy
UCSD July 20-21, 2009
Michael Galperin University of California at San Diego
UCSD July 20-21, 2009 – p.1
UCSD July 20-21, 2009 – p.2
E
L R
M = L+ R
UCSD July 20-21, 2009 – p.2
Science 319, 1056 (2008).
IETS RIETS e-e interaction Noise Non-linear conductance Heating Light-matter interaction
UCSD July 20-21, 2009 – p.2
Nano Lett. 9, 758 (2009)
UCSD July 20-21, 2009 – p.2
J.Zhang et al. Nano Lett. 7, 2101 (2007)
UCSD July 20-21, 2009 – p.3
S.W.Wu et al. Phys. Rev. B 77, 205430 (2008)
UCSD July 20-21, 2009 – p.3
S.Nie and S.R.Emory. Science 275, 1102 (1997)
UCSD July 20-21, 2009 – p.3
D.R.Ward et al. Nano Lett. 8, 919 (2008)
UCSD July 20-21, 2009 – p.3
Z.Ioffe et al. Nature Nanotechnology 3, 727 (2008)
UCSD July 20-21, 2009 – p.3
Nano Lett. 9, 758 (2009); J. Chem. Phys. 130, 144109 (2009)
UCSD July 20-21, 2009 – p.4
L R
UCSD July 20-21, 2009 – p.4
L R
|2> |1> e e e
UCSD July 20-21, 2009 – p.4
EF |2> |1> e e
UCSD July 20-21, 2009 – p.4
|2> |1> e ph |2> |1> e ph
UCSD July 20-21, 2009 – p.4
mˆ
kˆ
αˆ
km
kˆ
kk′
kˆ
2ˆ
α
2ˆ
mk
km′
2 1
MK,mm′ = −iδmm′ΓMK,m/2
MK,mm′(E) = iδmm′fK(E)ΓMK,m
MK,mm′(E) = −iδmm′[1 − fK(E)]ΓMK,m
UCSD July 20-21, 2009 – p.4
k=k′∈K
kk′
2 1
UCSD July 20-21, 2009 – p.4
Projections Σ<
NK,mm(E) =
dω 2π BNK(ω, µK) G<
¯ m ¯ m(E + ω)
Σ>
NK,mm(E) =
dω 2π BNK(ω, µK) G>
¯ m ¯ m(E − ω)
with BNK(ω, µK) = 2π
kk′
×δ(E − εk)δ(E + ω − εk′)fK(E)[1 − fK(E + ω)] ≡ 2π
2 ρe−h
K (ω)
UCSD July 20-21, 2009 – p.4
NK = iBNK
NK = −iBNK
UCSD July 20-21, 2009 – p.4
α
α
2 1
UCSD July 20-21, 2009 – p.4
P(E) = α
α
22(E + ωα)
11(E − ωα)
P(E) = α
α
22(E + ωα)
11(E − ωα)
UCSD July 20-21, 2009 – p.4
P
P
α
α
UCSD July 20-21, 2009 – p.4
−∞
<(E) G>(E) − ΣB >(E) G<(E)]
UCSD July 20-21, 2009 – p.5
General expression Iabs(ω0) = +∞
−∞
dE 2π
P0,22(E) G> 22(E) − Σ> P0,22(E) G< 22(E)
ε1 ≪ µL,R ≪ ε2 (low bias) coupling to the photon field is weak Γ1,2 ≪ ε21, |ε1,2 − EF| Iabs(ω0) =
(ε2 − ω0 − ε1)2 + (Γ/2)2 × ΓM,1ΓM,2 Γ1Γ2
UCSD July 20-21, 2009 – p.6
1 2 3 4
0 (eV)
4 8 12 ( 10
+9)
Iabs (photons/s) =1.5 V =1.2 V =1.1 V =1.0 V =0.5 V
ε21 = 2 eV γP = 10−6 eV T = 300 K BNL = BNR = 0.1 eV ΓML/R,1 = 0.01 eV ΓML/R,2 = 0.2 eV
UCSD July 20-21, 2009 – p.6
Radiation field in resonance with the molecular optical transition Molecules with strong charge-transfer transitions
DMEANS (4-Dimethylamino-4’-nitrostilbene) 7D (ground) → 31D (first excited singlet) all-trans Retinal in Poly-methyl methacrylate films 6.6D → 19.8D (1Bu electronic state) 40Å CdSe nanocrystals 0D → 32D (first excited state)
If optical charge transfer is parallel to the wire axis
UCSD July 20-21, 2009 – p.7
−∞
ML(E) G>(E) − Σ> ML(E) G<(E)]
Simplified version (ω0 ∼ ε21, Φ = 0, Γ1,2 ≪ ε21) Isd = |V (P) |2
(ε2 − ω0 − ε1)2 + (Γ/2)2 ΓML,1ΓMR,2 − ΓML,2ΓMR,1 Γ1Γ2 The yield of the effect Yc = Isd Iabs
= ΓML,1ΓMR,2 − ΓML,2ΓMR,1 ΓM,1ΓM,2
UCSD July 20-21, 2009 – p.7
1 2 3 4
0 (eV)
1 3 5 ( 10
Isd (A)
Φ = 0 ε21 = 2 eV γP = 10−6 eV T = 300 K BNL = BNR = 0.1 eV ΓML/R,1 = 0.2 eV ΓML,2 = 0.02 eV ΓMR,2 = 0.3 eV V (P) = 10−3 eV
UCSD July 20-21, 2009 – p.7
0.0 0.5 1.0
(V)
5 10 15 20 ( 10-7)
Isd (A)
0=1.5 eV 0=1.7 eV 0=1.9 eV
Φ = 0 ε21 = 2 eV γP = 10−6 eV T = 300 K BNL = BNR = 0.1 eV ΓML,1/MR,2 = 0.2 eV ΓML,2/MR,1 = 0.02 eV V (P) = 0.02 eV
UCSD July 20-21, 2009 – p.7
UCSD July 20-21, 2009 – p.8
Frequency resolved spectrum I′
em(ω)
= ρP(ω) +∞
−∞
dE 2π
P,11(E, ω) G> 11(E) − Σ> P,11(E, ω) G< 11(E)
Itot
em
= ∞ dω I′
em(ω)
= +∞
−∞
dE 2π
P,11(E) G> 11(E) − Σ> P,11(E) G< 11(E)
When coupling to radiation field is weak and Γ1,2 ≪ ε21 I′
em(ω)
= γP(ω)
−∞
dE 2π fL(E + ω)ΓML,2 + fR(E + ω)ΓMR,2 (E + ω − ε2)2 + (Γ2/2)2 × [1 − fL(E)]ΓML,1 + [1 − fR(E)]ΓMR,1 (E − ε1)2 + (Γ1/2)2
em
= γP(ε21)
UCSD July 20-21, 2009 – p.8
When in addition µL ≫ ε2 > ε1 ≫ µR Itot
em = γP
Γ1Γ2 Also in this case Isd = 1
ΓML,mΓMR,m Γm + BN + γP
Γ1Γ2 So that the yield Yem = Itot
em
Isd = γP
ΓMR,2 ΓMR,1Γ1 + ΓML,1 ΓML,2Γ2 + BN + γP
UCSD July 20-21, 2009 – p.8
L R
1 2 3
L R
1 2 3
UCSD July 20-21, 2009 – p.8
5 10 15 ( 10
7)
Iem
tot (s
1 2 ( 10
Isd (A)
BNL=BNR=0.1 eV 5 10 15 ( 10
7)
Iem
tot (s
1 2 ( 10
Isd (A)
BNL=BNR=1 eV 1 2 3 4 5 6
(V)
2 6 10 ( 10
Iem
tot/Isd BNL=BNR=1.0 eV BNL=BNR=0.1 eV
T = 300 K ε21 = 2 eV ΓMK,m = 0.1 eV γP = 10−6 eV BNL = BNR = 0.1 eV
UCSD July 20-21, 2009 – p.8
1 2 3 4
(eV)
1 3 5 7 ( 10
9)
Iem’( )
=5 V =3 V =2 V =1 V
UCSD July 20-21, 2009 – p.8
1 2 3 4 5 6
(V)
1 2 ( 108)
dIem
tot/d (s
MK 3
BNK 10
1 2 3 4
(eV)
1 3 5 7 ( 109)
Iem’( )
=3 V
UCSD July 20-21, 2009 – p.8
1 2 3 4 5 6
(V)
4 8 12 16 ( 107)
Iem
tot (s
=0.99 =0.66 =0.50
η = ΦL/Φ = ΓMR,m/Γm = BNR/BN η → 1 no emission (either LUMO is empty
Fluorescence in STM Φ should fall at the molecule-substrate interface
Spacers reduce energy losses into substrate (BN) But enable light emission at the molecule
UCSD July 20-21, 2009 – p.8
MG, A.Nitzan, M.A.Ratner, PRL 96, 166803 (2006)
UCSD July 20-21, 2009 – p.9
L R
|2> |1> e e e
Fluxes considered electronic current through the molecule energy flow between the molecule and electron-hole excitations in the leads
UCSD July 20-21, 2009 – p.9
EF |2> |1> e e
Fluxes considered electronic current through the molecule energy flow between the molecule and electron-hole excitations in the leads
UCSD July 20-21, 2009 – p.9
mˆ
kˆ
km
kˆ
kk′
kˆ
2ˆ
Isd = IL
sd + Ie−h sd
IL
sd
= e
−∞
dE 2π
Γ(ML)
m
Gr
mm(E)Γ(MR) m
Ga
mm(E)
× [fL(E) − fR(E)] Ie−h
sd
= e B
2
1
Γ1 − n(MR)
1
2
1
Γ1 − n(ML)
1
m
≪ ε21 is assumed in the last
UCSD July 20-21, 2009 – p.9
sd = e
m
m
sd
2
1
1
2
L R
UCSD July 20-21, 2009 – p.9
L R
UCSD July 20-21, 2009 – p.9
L R
UCSD July 20-21, 2009 – p.9
DMEANS (4-Dimethylamino-4’-nitrostilbene) 7D (ground) → 31D (first excited singlet) all-trans Retinal in Poly-methyl methacrylate films 6.6D → 19.8D (1Bu electronic state) 40Å CdSe nanocrystals 0D → 32D (first excited state)
UCSD July 20-21, 2009 – p.9
1 ( 10
Isd (A)
e-h Landauer Isd
1,2 (ML/R)=0.1 eV
2 4 6 e (eV)
1 ( 10
Isd (A)
2 (MR)=0.01 eV 2 (ML)=0.19 eV 1 (ML/R)=0.1 eV
T = 300 K ε1 = 0 eV ε2 = 2 eV Γ(M)
1,2 = 0.2 eV
UCSD July 20-21, 2009 – p.9
m
m
m
UCSD July 20-21, 2009 – p.9
2 4 6 8 10 12 14 R (Angstrom) 10
10
10
10
10-7 10-6 10
Isd (A)
e-h Landauer Isd =3 V
A(ML/R)
1
= 0.27 eV A(ML)
2
= 0.52 eV A(MR)
2
= 0.027 eV α(MK)
m
= 1 Å−1 β(K) = 0.01 eV
UCSD July 20-21, 2009 – p.9
D.R.Ward et al. Nano Lett. 8, 919 (2008)
UCSD July 20-21, 2009 – p.10
Nano Lett. 9, 758 (2009); J. Chem. Phys. 130, 144109 (2009)
UCSD July 20-21, 2009 – p.11
ˆ H0 =
εm ˆ d†
m ˆ
dm + ωvˆ b†
vˆ
bv +
εkˆ c†
kˆ
ck +
ωβˆ b†
βˆ
bβ +
ναˆ a†
αˆ
aα ˆ V (e−v) =
V (e−v)
m
ˆ Qv ˆ d†
m ˆ
dm ˆ V (et) =
km ˆ
c†
k ˆ
dm + V (et)
mk ˆ
d†
mˆ
ck
V (v−b) =
U (v−b)
β
ˆ Qv ˆ Qβ ˆ V (e−h) =
k1k2
ˆ d†
1 ˆ
d2ˆ c†
k1ˆ
ck2 + H.c.
α
1 ˆ
∗
α
α ˆ
2 ˆ
kmˆ
k ˆ
mk ˆ
mˆ
α)
UCSD July 20-21, 2009 – p.11
Introducing excitation operators ˆ D ≡ ˆ d†
1 ˆ
d2 ˆ Dmk ≡ ˆ d†
mˆ
ck m = 1, 2 k ∈ {L, R} and after small polaron transformation ˆ V (e−p) =
a†
α ˆ
Oα + ˆ O†
αˆ
aα
Oα =
∗
U (e−p)
α
ˆ D ˆ X +
1k ˆ
D1k ˆ X†
1 + V α k2 ˆ
Dk2 ˆ X2 + V α
k1 ˆ
Dk1 ˆ X1 + V α
2k ˆ
D2k ˆ X†
2
O(M)
α
+ ˆ O(1)
α + ˆ
O(2)
α + ˆ
O(3)
α + ˆ
O(4)
α
UCSD July 20-21, 2009 – p.11
α(t)ˆ
−∞
α (t, t′) G> α (t′, t) + G> α (t, t′) F < α (t′, t)
α (t, t′) G< α (t′, t) − G< α (t, t′) F > α (t′, t)]
α(τ ′) >
α(τ ′) >
UCSD July 20-21, 2009 – p.12
−∞
f (t′ − t) G< f (t − t′)
UCSD July 20-21, 2009 – p.12
f (t − t′)
−∞
f (t′ − t) Fi(τ1, τ2)
f(t′) ˆ
i(τ1) ˆ
O)
UCSD July 20-21, 2009 – p.12
i (t1 − t2) = −ie−iνi(t1−t2)
f (t′ − t) = −ie−iνf(t′−t)
UCSD July 20-21, 2009 – p.12
i→f =
−∞
−∞
−∞
f(t′) ˆ
i(t1) >
t1 t2 t’ t
UCSD July 20-21, 2009 – p.12
i→f =
−∞
t
t′
f(t′) ˆ
i(t1) ˆ
t t’ t2 t1
UCSD July 20-21, 2009 – p.12
i→f
−∞
−∞
t′
f(t′) ˆ
i(t1) >
t1 t t2 t2 t t1 t’
UCSD July 20-21, 2009 – p.12
1 2 3
V (V)
1 2 3 10
I (a.u.)
8 10 12 10
i
{ f} (a.u.)
¯ J
4
4
2
UCSD July 20-21, 2009 – p.13
1 2 3
V (V)
1 2 3 10-14
i i v (a.u.)
anti-Stokes Stokes
¯ J
UCSD July 20-21, 2009 – p.13
1
( f- i)/
v
6 12 10-14
i f (a.u.)
V=0 V V=2.5 V
¯ J
UCSD July 20-21, 2009 – p.13
1
( f- i)/
v
5 10 15 10-14
i f (a.u.) (e-h)=0.004 eV (e-h)=0.010 eV
¯ J
UCSD July 20-21, 2009 – p.13
1 2
( i- )/
v
6 10 10-13
i
{ f} (a.u.)
V=0 V V=2.5 V
¯ J
UCSD July 20-21, 2009 – p.13
0.0 0.01 0.02
(e-h) (eV)
0.0
/
0.0 0.01 0.02
m (eV)
/
0.001
¯ J(intR)
νi→νi−ωv ¯
Jνi→νi−ωv ¯ J(intR)
νi→νi−ωv ¯
Jνi→νi−ωv
UCSD July 20-21, 2009 – p.13
1 2 3
V (V)
500 1500
T (K) TS-aS Tv
UCSD July 20-21, 2009 – p.13
1.0 1.01 1.02
TS-aS/Tv
0.005 0.01 0.015 0.02
m (eV)
0.0 0.01 0.02
(e-h) (eV)
0.0
/
0.0 0.01 0.02
m (eV)
/
0.001
¯ J(intR)
νi→νi−ωv ¯
Jνi→νi−ωv ¯ J(intR)
νi→νi−ωv ¯
Jνi→νi−ωv
UCSD July 20-21, 2009 – p.13
J.R.Lombardi et al. JCP 84, 4174 (1986)
UCSD July 20-21, 2009 – p.14
Origin of the (metal-to-molecule) peak dE1 2π dE2 2π
S(K)<
2
(E1) G>
2 (E2)
νi + E1 + ωvvin − E2 − ωvv′ + iΓ(e−h)/2 At µL = µR = EF and for T → 0 integral on E1 yields ln
D where D is leads half-bandwidth. This gives a peak at νi = E2 − EF − ωv(vin − v′)
UCSD July 20-21, 2009 – p.14
1.0 1.1 1.2
i (eV)
10 20 10-20
i i- v (a.u.)
20
¯ J
UCSD July 20-21, 2009 – p.14
0.0 0.1 0.2
m (eV)
4 8 12 10-21
i=1eV
molecule met-to-mol 0.01 0.02 0.03
m (eV)
4 8 10-23
i=0.5eV
i i- v
¯ J
UCSD July 20-21, 2009 – p.14
UCSD July 20-21, 2009 – p.15
UCSD July 20-21, 2009 – p.15