physics-based device simulation - tu dresden · 2015-12-16 · pros & cons 4 i realistic device...
TRANSCRIPT
1
Physics-based Device Simulation
C. Jungemann
Institut für Theoretische ElektrotechnikRWTH Aachen University
52056 Aachen
2OutlineI MotivationI Pros & ConsI Semiclassical transportI ResultsI Conclusions
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
3MotivationSimulation hierarchy
I Quantum transportI Semiclassical transportI TCADI Compact models
Speed
Acc
urac
y
Com
plexity
Why semiclassical transport simulation (Boltzmann equation)?
I Complex band structuresI Detailed scattering modelsI Only material-dependent parametersI DC, AC and noise analysisI No adjustable parameters in the RF noise model
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
4Pros & ConsI Realistic device structureI Microscopic, accurate physics (Boltzmann equation)I Detailed physicsI Consistent modeling of transport and noiseI Numerically challengingI Very slow (a CPU day per bias point for a 2D device)I Limited device detail (2D, only few parasitics)I Only as good as the models
5
Semiclassical transport
6Semiclassical transportPath of an electron in a semiconductor
0 100 200 300 400
50
100
200
300
X [nm]
Y [
nm
]
Newton’s equations of motion
d~pdt
= ~F ,d~rdt
= ~v
Two first order ODEsTwo initial conditions (position, momentum)
6Semiclassical transportPath of an electron in a semiconductor
0 100 200 300 400
50
100
200
300
X [nm]
Y [
nm
]
Newton’s equations of motion
d~pdt
= ~F ,d~rdt
= ~v
Two first order ODEsTwo initial conditions (position, momentum)
6Semiclassical transportPath of an electron in a semiconductor
0 100 200 300 400
50
100
200
300
X [nm]
Y [
nm
]
Newton’s equations of motion
d~pdt
= ~F ,d~rdt
= ~v
Two first order ODEsTwo initial conditions (position, momentum)
6Semiclassical transportPath of an electron in a semiconductor
0 100 200 300 400
50
100
200
300
X [nm]
Y [
nm
]
Newton’s equations of motion
d~pdt
= ~F ,d~rdt
= ~v
Two first order ODEsTwo initial conditions (position, momentum)
6Semiclassical transportPath of an electron in a semiconductor
0 100 200 300 400
50
100
200
300
X [nm]
Y [
nm
]
Newton’s equations of motion
d~pdt
= ~F ,d~rdt
= ~v
Two first order ODEsTwo initial conditions (position, momentum)
7Semiclassical TransportParticle ensemble described by a distribution function (~p = ~~k )
f (~r , ~k , t)
Number of particles in a volume of phase space
dN =2
(2π)3 f (~r , ~k , t)d3rd3k
Equilibrium (Boltzmann distribution)
f0 = c exp
(−ε(~r , ~k)
kT
)
7Semiclassical TransportParticle ensemble described by a distribution function (~p = ~~k )
f (~r , ~k , t)
Number of particles in a volume of phase space
dN =2
(2π)3 f (~r , ~k , t)d3rd3k
Equilibrium (Boltzmann distribution)
f0 = c exp
(−ε(~r , ~k)
kT
)
7Semiclassical TransportParticle ensemble described by a distribution function (~p = ~~k )
f (~r , ~k , t)
Number of particles in a volume of phase space
dN =2
(2π)3 f (~r , ~k , t)d3rd3k
Equilibrium (Boltzmann distribution)
f0 = c exp
(−ε(~r , ~k)
kT
)
8Semiclassical TransportBoltzmann Transport Equation
∂f(~r , ~k , t
)∂t
− q~~E(~r) · ∇~k f
(~r , ~k , t
)+ ~v
(~k)· ∇~r f
(~r , ~k , t
)=
Ωsys
(2π)3
∫W(~r , ~k |~k ′
)f(~r , ~k ′, t
)−W
(~r , ~k ′|~k
)f(~r , ~k , t
)d3k ′
Integro-differential equation
Direct discretization of the 6D phase space is not possible
Spherical Harmonics Expansion solver (SHE) for the BTE(Bologna group, Maryland group et al.)
8Semiclassical TransportBoltzmann Transport Equation
∂f(~r , ~k , t
)∂t
− q~~E(~r) · ∇~k f
(~r , ~k , t
)+ ~v
(~k)· ∇~r f
(~r , ~k , t
)=
Ωsys
(2π)3
∫W(~r , ~k |~k ′
)f(~r , ~k ′, t
)−W
(~r , ~k ′|~k
)f(~r , ~k , t
)d3k ′
Integro-differential equation
Direct discretization of the 6D phase space is not possible
Spherical Harmonics Expansion solver (SHE) for the BTE(Bologna group, Maryland group et al.)
8Semiclassical TransportBoltzmann Transport Equation
∂f(~r , ~k , t
)∂t
− q~~E(~r) · ∇~k f
(~r , ~k , t
)+ ~v
(~k)· ∇~r f
(~r , ~k , t
)=
Ωsys
(2π)3
∫W(~r , ~k |~k ′
)f(~r , ~k ′, t
)−W
(~r , ~k ′|~k
)f(~r , ~k , t
)d3k ′
Integro-differential equation
Direct discretization of the 6D phase space is not possible
Spherical Harmonics Expansion solver (SHE) for the BTE(Bologna group, Maryland group et al.)
8Semiclassical TransportBoltzmann Transport Equation
∂f(~r , ~k , t
)∂t
− q~~E(~r) · ∇~k f
(~r , ~k , t
)+ ~v
(~k)· ∇~r f
(~r , ~k , t
)=
Ωsys
(2π)3
∫W(~r , ~k |~k ′
)f(~r , ~k ′, t
)−W
(~r , ~k ′|~k
)f(~r , ~k , t
)d3k ′
Integro-differential equation
Direct discretization of the 6D phase space is not possible
Spherical Harmonics Expansion solver (SHE) for the BTE(Bologna group, Maryland group et al.)
9Spherical Harmonics ExpansionSpherical coordinates for k -space
(kx , ky , kz) → (k , ϑ, ϕ) → (ε, ϑ, ϕ)
Dependence on angles is expanded with spherical harmonics Yl,m(ϑ, ϕ)
Y0,0 =1√4π
Y1,−1 =
√3
4πsinϑ sinϕ
Y1,0 =
√3
4πcosϑ
Y1,1 =
√3
4πsinϑ cosϕ
∮Yl,m(ϑ, ϕ)Yl ′,m′(ϑ, ϕ)dΩ = δl,l ′δm,m′ , dΩ = sinϑdϑdϕ
9Spherical Harmonics ExpansionSpherical coordinates for k -space
(kx , ky , kz) → (k , ϑ, ϕ) → (ε, ϑ, ϕ)
Dependence on angles is expanded with spherical harmonics Yl,m(ϑ, ϕ)
Y0,0 =1√4π
Y1,−1 =
√3
4πsinϑ sinϕ
Y1,0 =
√3
4πcosϑ
Y1,1 =
√3
4πsinϑ cosϕ
∮Yl,m(ϑ, ϕ)Yl ′,m′(ϑ, ϕ)dΩ = δl,l ′δm,m′ , dΩ = sinϑdϑdϕ
9Spherical Harmonics ExpansionSpherical coordinates for k -space
(kx , ky , kz) → (k , ϑ, ϕ) → (ε, ϑ, ϕ)
Dependence on angles is expanded with spherical harmonics Yl,m(ϑ, ϕ)
Y0,0 =1√4π
Y1,−1 =
√3
4πsinϑ sinϕ
Y1,0 =
√3
4πcosϑ
Y1,1 =
√3
4πsinϑ cosϕ
∮Yl,m(ϑ, ϕ)Yl ′,m′(ϑ, ϕ)dΩ = δl,l ′δm,m′ , dΩ = sinϑdϑdϕ
9Spherical Harmonics ExpansionSpherical coordinates for k -space
(kx , ky , kz) → (k , ϑ, ϕ) → (ε, ϑ, ϕ)
Dependence on angles is expanded with spherical harmonics Yl,m(ϑ, ϕ)
Y0,0 =1√4π
Y1,−1 =
√3
4πsinϑ sinϕ
Y1,0 =
√3
4πcosϑ
Y1,1 =
√3
4πsinϑ cosϕ
∮Yl,m(ϑ, ϕ)Yl ′,m′(ϑ, ϕ)dΩ = δl,l ′δm,m′ , dΩ = sinϑdϑdϕ
9Spherical Harmonics ExpansionSpherical coordinates for k -space
(kx , ky , kz) → (k , ϑ, ϕ) → (ε, ϑ, ϕ)
Dependence on angles is expanded with spherical harmonics Yl,m(ϑ, ϕ)
Y0,0 =1√4π
Y1,−1 =
√3
4πsinϑ sinϕ
Y1,0 =
√3
4πcosϑ
Y1,1 =
√3
4πsinϑ cosϕ
∮Yl,m(ϑ, ϕ)Yl ′,m′(ϑ, ϕ)dΩ = δl,l ′δm,m′ , dΩ = sinϑdϑdϕ
10Semiclassical transportSpherical harmonics expansion of the distribution function
gl,m(~r , ε, t) =1
(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ)f (~r , ~k , t)d3k
g(~r , ~k , t) = g(~r , ε, ϑ, ϕ, t) ≈lmax∑l=0
m=l∑m=−l
gl,m(~r , ε, t)Yl,m(ϑ, ϕ)
Electron density
n(~r , t) =2
(2π)3
∫f (~r , ~k , t)d3k = 2
√4π∫ ∞
0g0,0(~r , ε, t)dε
Electron current density (spherical bands)
~j(~r , t) = 2
√4π3
∫ ∞0
v(ε)
g1,1(~r , ε, t)g1,−1(~r , ε, t)g1,0(~r , ε, t)
dε
10Semiclassical transportSpherical harmonics expansion of the distribution function
gl,m(~r , ε, t) =1
(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ)f (~r , ~k , t)d3k
g(~r , ~k , t) = g(~r , ε, ϑ, ϕ, t) ≈lmax∑l=0
m=l∑m=−l
gl,m(~r , ε, t)Yl,m(ϑ, ϕ)
Electron density
n(~r , t) =2
(2π)3
∫f (~r , ~k , t)d3k = 2
√4π∫ ∞
0g0,0(~r , ε, t)dε
Electron current density (spherical bands)
~j(~r , t) = 2
√4π3
∫ ∞0
v(ε)
g1,1(~r , ε, t)g1,−1(~r , ε, t)g1,0(~r , ε, t)
dε
10Semiclassical transportSpherical harmonics expansion of the distribution function
gl,m(~r , ε, t) =1
(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ)f (~r , ~k , t)d3k
g(~r , ~k , t) = g(~r , ε, ϑ, ϕ, t) ≈lmax∑l=0
m=l∑m=−l
gl,m(~r , ε, t)Yl,m(ϑ, ϕ)
Electron density
n(~r , t) =2
(2π)3
∫f (~r , ~k , t)d3k = 2
√4π∫ ∞
0g0,0(~r , ε, t)dε
Electron current density (spherical bands)
~j(~r , t) = 2
√4π3
∫ ∞0
v(ε)
g1,1(~r , ε, t)g1,−1(~r , ε, t)g1,0(~r , ε, t)
dε
10Semiclassical transportSpherical harmonics expansion of the distribution function
gl,m(~r , ε, t) =1
(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ)f (~r , ~k , t)d3k
g(~r , ~k , t) = g(~r , ε, ϑ, ϕ, t) ≈lmax∑l=0
m=l∑m=−l
gl,m(~r , ε, t)Yl,m(ϑ, ϕ)
Electron density
n(~r , t) =2
(2π)3
∫f (~r , ~k , t)d3k = 2
√4π∫ ∞
0g0,0(~r , ε, t)dε
Electron current density (spherical bands)
~j(~r , t) = 2
√4π3
∫ ∞0
v(ε)
g1,1(~r , ε, t)g1,−1(~r , ε, t)g1,0(~r , ε, t)
dε
11Semiclassical transportSpherical harmonics expansion of the Boltzmann equation
1(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ) BEd3k
yields balance equations over energy/real space
∂gl,m
∂t− q~E ·
∂~jl,m∂ε
+∇~r ·~jl,m − Γl,m = Wl,mg
with
~jl,m(~r , ε, t) = v(ε)∑l ′,m′
~al,m,l ′,m′gl ′.m′(~r , ε, t)
The balance equations are coupled
11Semiclassical transportSpherical harmonics expansion of the Boltzmann equation
1(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ) BEd3k
yields balance equations over energy/real space
∂gl,m
∂t− q~E ·
∂~jl,m∂ε
+∇~r ·~jl,m − Γl,m = Wl,mg
with
~jl,m(~r , ε, t) = v(ε)∑l ′,m′
~al,m,l ′,m′gl ′.m′(~r , ε, t)
The balance equations are coupled
11Semiclassical transportSpherical harmonics expansion of the Boltzmann equation
1(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ) BEd3k
yields balance equations over energy/real space
∂gl,m
∂t− q~E ·
∂~jl,m∂ε
+∇~r ·~jl,m − Γl,m = Wl,mg
with
~jl,m(~r , ε, t) = v(ε)∑l ′,m′
~al,m,l ′,m′gl ′.m′(~r , ε, t)
The balance equations are coupled
11Semiclassical transportSpherical harmonics expansion of the Boltzmann equation
1(2π)3
∫δ(ε− ε(~k)
)Yl,m(ϑ, ϕ) BEd3k
yields balance equations over energy/real space
∂gl,m
∂t− q~E ·
∂~jl,m∂ε
+∇~r ·~jl,m − Γl,m = Wl,mg
with
~jl,m(~r , ε, t) = v(ε)∑l ′,m′
~al,m,l ′,m′gl ′.m′(~r , ε, t)
The balance equations are coupled
12Semiclassical transportH-transform
5eV
25eV
0eV x
H
H(~r , ~k
)= ε
(~k)− qΨ
(~r)⇒ g′(~r ,H, t) = g(~r , ε, t)
Cancels the energy derivative for the stationary state
H
HHHHHH
−q~E ·∂~j ′l,m∂ε
+∇~r · ~j ′l,m − Γ′l,m = Wl,mg′
Captures ballistic transport
12Semiclassical transportH-transform
5eV
25eV
0eV x
H
H(~r , ~k
)= ε
(~k)− qΨ
(~r)⇒ g′(~r ,H, t) = g(~r , ε, t)
Cancels the energy derivative for the stationary state
H
HHHHHH
−q~E ·∂~j ′l,m∂ε
+∇~r · ~j ′l,m − Γ′l,m = Wl,mg′
Captures ballistic transport
12Semiclassical transportH-transform
5eV
25eV
0eV x
H
H(~r , ~k
)= ε
(~k)− qΨ
(~r)⇒ g′(~r ,H, t) = g(~r , ε, t)
Cancels the energy derivative for the stationary state
H
HHHHHH
−q~E ·∂~j ′l,m∂ε
+∇~r · ~j ′l,m − Γ′l,m = Wl,mg′
Captures ballistic transport
12Semiclassical transportH-transform
5eV
25eV
0eV x
H
H(~r , ~k
)= ε
(~k)− qΨ
(~r)⇒ g′(~r ,H, t) = g(~r , ε, t)
Cancels the energy derivative for the stationary state
H
HHHHHH
−q~E ·∂~j ′l,m∂ε
+∇~r · ~j ′l,m − Γ′l,m = Wl,mg′
Captures ballistic transport
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
13Semiclassical transportNumerics
I H-transformI Maximum Entropy Dissipation SchemeI Dimensional splittingI Box Integration (exact particle number conservation)I First order expansion has property MI Newton-Raphson method to solve BE and PE self-consistently
14Semiclassical transport
Complex band structure effects
First conduction band of Si Density of states of Si
0 1 2 3 4 5 6
0
1
2
3·1028
Energy [eV]
DO
S[1/
eVcm
3 ]
Calculated by the nonlocal empirical pseudopotential method
Only DOS and v2 are required for SHE
14Semiclassical transport
Complex band structure effects
First conduction band of Si Density of states of Si
0 1 2 3 4 5 6
0
1
2
3·1028
Energy [eV]
DO
S[1/
eVcm
3 ]
Calculated by the nonlocal empirical pseudopotential method
Only DOS and v2 are required for SHE
15Semiclassical transportElectron-phonon scattering rate of silicon at 300K
0 1 2 3 4 5 6
0
1
2
3·1014
Energy [eV]
Sca
tterin
gra
te[1/s
]
All parameters determined by theory or basic experiments
16Semiclassical transportElectron mobility
0 100 200 300 400 500
103
104
105
106
Temperature [K]
Low
-fiel
dm
obili
ty[c
m2 /
Vs]
BECanali
Electron velocity at 300Kin < 111 > direction
10−1 100 101 102105
106
107
Electric field [kV/cm]D
riftv
eloc
ity[c
m/
s]
BECanali
Good agreement between simulation and experiment for silicon
17Semiclassical transport
Impact ionization
Ionization coefficient
2 4 6 8 10100
101
102
103
104
105
Inverse electric field [cm/kV]
IIco
effic
ient
[1/c
m]
BEExp.
Quantum yield
1 1.5 2 2.5 3
10−8
10−6
10−4
10−2
100
Energy [eV]
Qua
ntum
yiel
dBE
DiMaria
Good agreement between simulation and experiment for silicon
18
Results
19
Bipolar transistor
20Bipolar transistor
Doping profile Collector current
SHE can handle high doping concentrations without problems
SHE can handle small currents without problems
20Bipolar transistor
Doping profile Collector current
SHE can handle high doping concentrations without problems
SHE can handle small currents without problems
21Bipolar transistor
VBE = 0.55V , VCE = 0.5VEnergy distribution function
VCE = 0.5VElectron density
SHE can handle huge variations in the density without problems
22Bipolar transistor
Dependence on the maximum order of SHEVBE = 0.85V , VCE = 0.5V
Error in current Electron velocity
Transport in nanometric devices requires at least 5th order SHE
23
N+NN+ structure
24N+NN+ structure
25N+NN+ structure
26N+NN+ structure
27N+NN+ structure
28N+NN+ structure
29
-0.30 -0.20 -0.10 0.000 0.10 0.20 0.30
0.00
0.05
0.10
0.15
0.20
0.25
Silicon
Bottom oxide
GS D
Top oxide
Partially depleted SOI NMOSFET
30PDSOI NMOSFET
N+
N+
VSource
VBulk
VDrain
VGate
ISub
Impactionization
event
n−MOSFET
P
−4
−2
0
2
4
6
ε (
eV
)
L Γ X
k
[111] [100]
Every impact ionization event generates a secondary electron/hole pair
31PDSOI NMOSFET
0 0.2 0.4 0.6 0.8 10
5
10
15
20
25
Drain voltage [V]
Dra
incu
rren
t[A
/m]
with IIw/o II
Kink effect due to impact ionization (II) (Vgate = 1.0V )
Exact stationary solutions
31PDSOI NMOSFET
0 0.2 0.4 0.6 0.8 10
5
10
15
20
25
Drain voltage [V]
Dra
incu
rren
t[A
/m]
with IIw/o II
Kink effect due to impact ionization (II) (Vgate = 1.0V )
Exact stationary solutions
32PDSOI NMOSFET
0 0.2 0.4 0.6 0.8 10
5
10
15
20
25
Drain voltage [V]
Cur
rent
[A/m
]ID with IIID w/o II
10−16
10−15
10−14
10−13
10−12
10−11
10−10
II-curr.
About 17 orders of magnitude difference in currents at kink
33
pn Junction
34pn Junction
Simulation of avalanche breakdown
35 36 37 38 39
10−17
10−14
10−11
10−8
10−5
Reverse bias [V]
Cur
rent
[A/µ
m2 ]
35 36 37 38 39
0
1
2
3
·10−5
Reverse bias [V]
Cur
rent
[A/µ
m2 ]
Extremely steep increase of current with voltage (20µV/dec)
34pn Junction
Simulation of avalanche breakdown
35 36 37 38 39
10−17
10−14
10−11
10−8
10−5
Reverse bias [V]
Cur
rent
[A/µ
m2 ]
35 36 37 38 39
0
1
2
3
·10−5
Reverse bias [V]
Cur
rent
[A/µ
m2 ]
Extremely steep increase of current with voltage (20µV/dec)
35pn JunctionTransient DD simulation of switching a silicon pn-diode from 0 to -39V
100 101 102 10310−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
Time [ps]
Cur
rent
[A/µ
m2 ]
Avalanche breakdown is slow!
35pn JunctionTransient DD simulation of switching a silicon pn-diode from 0 to -39V
100 101 102 10310−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
Time [ps]
Cur
rent
[A/µ
m2 ]
Avalanche breakdown is slow!
36pn Junction
Simulation for 37.172V and 10pA/µm2
0
10
20
30
40
Pote
ntia
l[V
]0 1 2 3 4 510−14
10−5
104
1013
1022
Position [µm]
IIge
nera
tion
rate
[a.u
.]
elec.holesPot.
0 1 2 3 4 5101
105
109
1013
1017
Position [µm]
Den
sity
[cm−
3 ]
NANDelec.hole
II generation increases number of particles in the space charge region
37pn Junction
Increase of current by a factor of 1000 (dashed lines)
0 1 2 3 4 5101
105
109
1013
1017
Position [µm]
Den
sity
[cm−
3 ]
0 1 2 3 4 510−16
10−7
102
1011
1020
Position [µm]
IIge
nera
tion
rate
[µm−
1 C−
1 ]Particle density in the space charge region increases 1000 times
Relevant II generation rate scales with current
Hot electron/hole effects scale with current
37pn Junction
Increase of current by a factor of 1000 (dashed lines)
0 1 2 3 4 5101
105
109
1013
1017
Position [µm]
Den
sity
[cm−
3 ]
0 1 2 3 4 510−16
10−7
102
1011
1020
Position [µm]
IIge
nera
tion
rate
[µm−
1 C−
1 ]Particle density in the space charge region increases 1000 times
Relevant II generation rate scales with current
Hot electron/hole effects scale with current
37pn Junction
Increase of current by a factor of 1000 (dashed lines)
0 1 2 3 4 5101
105
109
1013
1017
Position [µm]
Den
sity
[cm−
3 ]
0 1 2 3 4 510−16
10−7
102
1011
1020
Position [µm]
IIge
nera
tion
rate
[µm−
1 C−
1 ]Particle density in the space charge region increases 1000 times
Relevant II generation rate scales with current
Hot electron/hole effects scale with current
37pn Junction
Increase of current by a factor of 1000 (dashed lines)
0 1 2 3 4 5101
105
109
1013
1017
Position [µm]
Den
sity
[cm−
3 ]
0 1 2 3 4 510−16
10−7
102
1011
1020
Position [µm]
IIge
nera
tion
rate
[µm−
1 C−
1 ]Particle density in the space charge region increases 1000 times
Relevant II generation rate scales with current
Hot electron/hole effects scale with current
38pn Junction
Increase of current by a factor of 1000 (dashed lines)
Scaled electron/hole distribution function in the middle of the pn junction
0 1 2 3 4 5 6 710−4
10−1
102
105
108
Energy [eV]
Dis
trib
utio
nfu
nctio
n[µ
m−
1 eV−
1 A−
1 ]
Hot electron/hole effects scale with current
38pn Junction
Increase of current by a factor of 1000 (dashed lines)
Scaled electron/hole distribution function in the middle of the pn junction
0 1 2 3 4 5 6 710−4
10−1
102
105
108
Energy [eV]
Dis
trib
utio
nfu
nctio
n[µ
m−
1 eV−
1 A−
1 ]
Hot electron/hole effects scale with current
39
Vertical power transistor
40Power MOSFETAvalanche breakdown of a power transistor
BE simulation for 35.83V and 10pA/µm2
Electrostatic Potential [V]
-0.6 6.8 14 22 36
e-Imp-Ioniz [cm^-3*s^-1]
2.6e+12 2.6e+16 2.6e+22
h-Imp-Ioniz [cm^-3*s^-1]
5.2e+13 5.2e+17 5.2e+23
45 · 106 variables, 2 CPU hours
41Power MOSFET
BE simulation for 35.83V and 10pA/µm2
Electrostatic Potential [V]
-0.6 6.8 14 22 36
0.3 0.4 0.5 0.6 0.7 0.8 0.9 110−2
101
104
107
1010
1013
Normalized path-length
Hig
h-en
ergy
part
icle
dens
ity[c
m−
3 ]1eV2eV3eV
Injection of electrons and holes into the oxide can be calculated
42
Structure, DC and AC behaviorof the npn THz SiGe HBT
43Comparison with experiments
Comparison for an ST-Microelectronics SiGe HBTat 10GHz and VCE = 1.2V
Cutoff frequency
10−1 100 101 102
101
102
Collector current [mA/µm2]
Cut
offf
requ
ency
[GH
z]
BESchröter
Minimum noise figure
100 1010
1
2
3
4
Collector current [mA/µm2]
Noi
sefig
ure
[dB
]
BESchröter
Good agreement between simulation and measurement
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
44Device structure
x
y0.0
10.0
20.0
30.0
90.0
100.00.0 10.0 20.0 30.0 40.0 50.0
OxideEmitter
Base
Collector
Emitter-window width: 2 × 25nm
I 149 nodes in x-directionI 20 nodes in y -directionI 5th order SHE expansion (21
spherical harmonics)I Up to 1.42 · 108 unknown
variables (spacing in energyis 5.17meV and number ofenergy grid points dependson bias)
I Over 200 GB of memoryrequired (for AC)
I Up to 105 CPU-seconds for aDC point
45Doping and germanium profile
0
5
10
15
20
25
30
Ge
cont
ent[
%]
0.01 0.02 0.03 0.04·10−2
1016
1017
1018
1019
1020
1021
x [µm]
Dop
ing
[cm−
3 ]
NANDGe
Base thickness: 7nm
46Input characteristic
101
102
103
Cur
rent
gain
0.4 0.5 0.6 0.7 0.8 0.9 1
10−8
10−6
10−4
10−2
100
102
Base/emitter bias [V]
Cur
rent
[mA/µ
m2 ]
ICIBIC/IB
VCE = 1.0VNo SRH recombination
47Output characteristic
10−3
10−2
10−1
100
101
Bas
ecu
rren
t[µ
A/µ
m2 ]
0 0.5 1 1.5 2 2.5 30
20
40
60
Collector/emitter bias [V ]
Col
lect
orcu
rren
t[µ
A/µ
m2 ]
IC|IB|
VBE = 0.7V
BVCE0 = 1.4V
48Gains
108 109 1010 1011 1012 1013
0
10
20
30
40
50
Frequency [Hz]
Gai
n[d
B]
|HCB|2UMSG/MAG
108 109 1010 1011 1012 1013
0
10
20
30
40
50
Frequency [Hz]
Gai
n[d
B]
|HCB|2UMSG/MAG
VBE = 0.88V, VCE = 1.0V, IC = 18.2mA/µm2
49Cutoff and max. oscil. frequencies
10−1 100 101 102
101
102
103
Collector current [mA/µm2]
Freq
uenc
y[G
Hz]
fT QSfT [email protected] [email protected]
VCE = 1.0V , fmax by extrapolation of Upeak fT = 1.2THz, peak fmax = 1.45THz
49Cutoff and max. oscil. frequencies
10−1 100 101 102
101
102
103
Collector current [mA/µm2]
Freq
uenc
y[G
Hz]
fT QSfT [email protected] [email protected]
VCE = 1.0V , fmax by extrapolation of Upeak fT = 1.2THz, peak fmax = 1.45THz
50Local contribution to the transit timefT =
12πτ
, τ =
∫dτdx
dx
5 10 15 20 25 30 35
−20
−10
0
10
20Base
x [nm]
dτ/dx
[fs/
nm]
IC = 54.6mA/µm2
IC = 19.6mA/µm2
VCE = 1.0V
51Minimum noise figure
10−1 100 101 102
0
2
4
6
Collector current [mA/µm2]
Min
imum
nois
efig
ure
[dB
] NFmin @ 1GHzNFmin @ 100GHzNFmin @ 300GHz
VCE = 1.0V
52Electrothermal Transport
High power densities lead to self-heating
Heating is non-local!
52Electrothermal Transport
High power densities lead to self-heating
Heating is non-local!
53Electrothermal Transport
Dirichlet boundary conditions
Scattering between optical and acoustic phonons goes both ways
53Electrothermal Transport
Dirichlet boundary conditions
Scattering between optical and acoustic phonons goes both ways
53Electrothermal Transport
Dirichlet boundary conditions
Scattering between optical and acoustic phonons goes both ways
54Electrothermal Transport
Junction temperature determines the current
Self-heating can lead to thermal runaway
54Electrothermal Transport
Junction temperature determines the current
Self-heating can lead to thermal runaway
54Electrothermal Transport
Junction temperature determines the current
Self-heating can lead to thermal runaway
55Electrothermal Transport
2D electrothermal simulation with SHE
Exact thermal simulation requires 3D simulation including metalization
55Electrothermal Transport
2D electrothermal simulation with SHE
Exact thermal simulation requires 3D simulation including metalization
56
Conclusions
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex
57ConclusionsI Semiclassical simulations give detailed physics at the price of high
CPU times and a 2D real spaceI Various effects can be included at a microscopic levelI DC, AC and noise analysisI Device degradation, self-heating, mechanical strain, ...I Transport at the nanoscale is very complex