development of fast methods for electronic structure ......hal id: pastel-00001655 submitted on 8...
TRANSCRIPT
![Page 1: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/1.jpg)
HAL Id: pastel-00001655https://pastel.archives-ouvertes.fr/pastel-00001655
Submitted on 8 Jun 2006
HAL is a multi-disciplinary open accessarchive for the deposit and dissemination of sci-entific research documents, whether they are pub-lished or not. The documents may come fromteaching and research institutions in France orabroad, or from public or private research centers.
L’archive ouverte pluridisciplinaire HAL, estdestinée au dépôt et à la diffusion de documentsscientifiques de niveau recherche, publiés ou non,émanant des établissements d’enseignement et derecherche français ou étrangers, des laboratoirespublics ou privés.
Development of fast methods for electronic structurecalculationsMaxime Barrault
To cite this version:Maxime Barrault. Development of fast methods for electronic structure calculations. domain_other.Ecole des Ponts ParisTech, 2005. English. pastel-00001655
![Page 2: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/2.jpg)
![Page 3: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/3.jpg)
![Page 4: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/4.jpg)
P
❯❱❲ ❳❨❩❬❭ ❪
❫ ❴❯ ❭ ❵
![Page 5: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/5.jpg)
![Page 6: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/6.jpg)
![Page 7: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/7.jpg)
❲
![Page 8: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/8.jpg)
❲ ❴ ❴❴ ❫❲ ❫ ❬ ❲ ❴❴ ❳❴ P ❬ ❳ ❩❬❭ ❪ ❴P ❬ ❩ ❲
❩ ❴ ❲ P ❴ ❲ ❴ ❨ ❬ P ❲ ❴ ❲❲ ❲ ❱❲ ❳❨ ❲ ❨
❲ P ❭ ❵ ❯ ❯ ❴ ❲ ❴ ❲ ❲ ❴ ❴ ❴ ❴ ❯ ❲
❲ ❨❴ ❩❬ ❫❨ ❳ ❫P❴ ❨ ❨ ❵❴ ❲ P ❱❲ ❳❨ P ❴
❲ ❲ ❨ ❫ ❳ ❬ ❯ ❯ ❭ ❵ ❱ ❲❵ ❵ ❴ ❭ ❱ ❭ ❨ ❭❬❯❲ ❭ ❴ ❵ ❳ ❩ ❨ ❵ P
![Page 9: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/9.jpg)
❲
❩ ❯ ❨❲ ❵ ❩ ❫ ❳ ❫ ❴ ❳❬ ❳ ❴ ❯ ❫ ❫❨❴ ❫ ❭ ❭ ❴ ❴ ❫ ❲
❴ P P ❴ ❲ ❴❯
![Page 10: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/10.jpg)
❲
❴❯ ❨❯ ❴❨❬❴
❯ ❯ ❴ P ❴❬ ❴ ❴ ❬ ❬ ❯ ❯
❯ ❲ ❯ ❴ ❴ ❲ ❨❯ ❴ ❨❯ ❴ ❫ ❴ ❲ ❴ ❨❯
❯ ❴ ❬❴ ❯ ❫ P ❴ ❨❯
H+2
H2
❯ ❬ ❨ ❯ ❨❯ ❫ ❴ ❯
❭ ❴ ❬ ❴ ❴❯ ❲P
❲ ❴ ❴ ❴ P ❴ P P ❨
❴ ❴ ❴ ❴ ❴ P ❴ P ❬❴❴ P ❴ P ❴
❴ ❴ P ❴ ❲ ❴❴ ❴ P ❴ P ❴ ❴ ❲ ❴ ❨ ❴ ❴ P ❴ ❨ P ❴ P ❲ P ❴ ❴ ❴ ❴ ❴ ❴ P ❨
❴ ❲ ❴ P ❴ ❴ ❴ ❬ ❨ ❴❲ ❴ ❨ H+
2
H2
P ❴ P ❴ ❴❴ ❴ P ❴ P ❴ ❨ ❴ ❬
![Page 11: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/11.jpg)
❲
❴ P ❨ P ❴ ❴ ❨ P ❴ P
❨ ❴ ❴ ❴ ❴ ❴ ❴ P P❴ ❴
![Page 12: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/12.jpg)
![Page 13: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/13.jpg)
![Page 14: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/14.jpg)
❭ ❯
❳❴ ❯
❫ ❯ ❳❴ ❯ P
❭ ❴ ❲N
❩ ❭ ❴ ❴ ❭ ❨❯
❪ ❨
❭ ❬ ❴ ❴❯
❭ ❯ ❴ ❬ ❴❯ ❭ ❨
❭ ❴ ❬
❯ ❨ ❴ ❯ P ❯ ❯
❭ ❴ ❭ ❳❯ ❲
❭ ❨❯ P ❬
❴ ❴❨
❴
Tz
❳❴ ❨❯
![Page 15: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/15.jpg)
❲
❭ ❴ ❭ ❴ ❨ ❫ P ❴ ❴
❫ P P ❫ P ❴ ❴ ❴ ❴
P ❴ ❴ P ❴ P ❴ P ❴ ❴ ❴ ❴
❲ ❲
❯ ❲ ❲
❴ ❨❯
❴ ❴ ❫ ❴ ❴ ❫❳❳
S 6= INb
❭ ❴
❭ ❯ ❴ ❳❲ ❪ ❯ ❯
❴❬❴ ❫ ❯
❴
![Page 16: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/16.jpg)
❫ µ
❫
uµ
❯ ❲ ❯
❭ H+
2
❯
❭ H2
❯
❭ ❬ ❭ ❴ ❭ ❨❯ ❭ ❴
![Page 17: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/17.jpg)
❲
![Page 18: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/18.jpg)
⋆
[
]❳ ❵ ❬ ❫
[
]❳ ❵ ❯ ❭ ❵
❩ P ❬ ❴❨ ❩ P ❴❨
[
]❳ ❵ ❱ ❳❨
❬ ❳❴ ❬
⋆ ❯
[
] ❯ ❴ ❯ ❭ ❬ ❳
[
] ❳❳ P ❨ ❬
[
] ❳ ❫ ❳❴ ❲ ❳ ❬ ❩
[
] ❩ ❳ ❬❳ ❭❨ ❫
![Page 19: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/19.jpg)
❲
![Page 20: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/20.jpg)
❭ ❴ ❴❯ ❴❨❬❴ ❬ ❨❯
N
M❨ ❴
(zk)k=1,M
(mk)k=1,M
T ❯
❯ ❨ ❲ ❯ ❨
(xk)k=1,M
❭ ❴ ❬ ❯
NT
P ❴❯ ❲ ❲ ❴ ❲ ❲
❫ ❯ ❲ ❨❯
N
me = 1, e = 1, ~ = 1,1
4πǫ0= 1.
❭ k
l M(k, l)
MS(k)
❨ kIk
kδij
❨ ❭ P ❲ ❴
❴ ❴
❭ P ❴ ❨❯ ❬
P Ψ(t, .)
❲ ❲ ❬ H L2(IRd)
❲d = 3(M+N)
❲ ❴ ❯ ❭ H ❨❯ P ❯ ❨
![Page 21: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/21.jpg)
❨ ❲ ❨❬ P ❭❨ ❲ ❲ ❨ ❴ P ❯
inf〈Ψe|He|ψe〉, Ψe ∈ He, ‖Ψe‖ = 1
He
❴ He
❬ L2(IR3N)
P ❯ ❨ ❬ ❨ ❴❬❴
NP φii=1,N
P ❯ ❯ P
infE(φ1, . . . , φN), φi ∈ H1(IR
3),
∫
IR3
φiφj = δij , ∀ 1 ≤ i, j ≤ N
❭ P E ❨xk ❭ ❯ P
❨ ∂E/∂xkk=1,M ❲ ❨ ❬ ❴
∀ 1 ≤ k ≤M, mkd2xk
dt2(t) = − ∂E
∂xk− ∂Wnuc
∂xk, Wnuc =
∑
1≤i<j≤M
zizj
|xi − xj |.
❭ P ❨
❯ P ❯ ❲ ❬❬ ❬❭ ❯ ❴
φi
❲ ❲ ❲ ❴ ❨❯ ❨❯ ❴ ❭ ❲
φi
![Page 22: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/22.jpg)
❭ ❴ ❯ P χµµ=1,Nb
❴❴φi
P
∀ 1 ≤ i ≤ N, φi =
Nb∑
µ=1
Cµiχµ.
Pχµ
C
❭ ❲ P ❯ ❲
infEd(CC
t), C ∈ M(Nb, N), CtSC = IN
❯ ❲
∣∣∣∣∣∣∣
ERHFd (D) = Tr(hD) +
(J(D)D
)− 1
2
(K(D)D
),
EKSd (D) = 2Tr(hD) + 2
(J(D)D
)+ Exc(D).❲ ∀ 1 ≤ µ, ν ≤ Nb,
Sµν =
∫
IR3
χµχν ,
hµν = hLµν + hV
µν =1
2
∫
IR3
∇χµ∇χν +
∫
IR3
Vxχµχν , Vx = −M∑
k=1
zk
|x− xk|,
J(X)µν =
Nb∑
κ,λ=1
(µν|κλ)Xκλ,
K(X)µν =
Nb∑
κ,λ=1
(µλ|νκ)Xκλ,
∀ 1 ≤ κ, λ ≤ Nb, (µν|κλ) =
∫
IR3
∫
IR3
χµ(x)χν(x)χκ(x′)χλ(x
′)
|x− x′| dx dx′.
❭ S
h
Vx
❨ Exc(D)
P ❲ P ❭❫
diteExc(D) =
∫
IR3
ρ(x)εLDAxc
(ρ(x)
)dx
❲ρ(x) = 2
Nb∑
µ,ν=1
Dµνχµ(x)χν(x).
![Page 23: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/23.jpg)
❭ (µν|κλ)
❭ O(N4
b )
❭ Ed
D = CCt ❲ ❬ ❴❬❴ ❴
φi ❴
P ❲ ❲
infEd(D), D ∈ MS(Nb), DSD = D,
(DS) = N
❭ ❬❭ ❲ ❲❬
❴ C
❭
F (D)C = SCE, E =❫
(ǫ1, . . . , ǫN ),CtSC = IN ,D = CCt.
❭ F (D)
❯
FRHF (D) = h + 2J(D) −K(D), FKS(D) = h + 2J(D) + Fxc(D)
❲∀ 1 ≤ µ, ν ≤ Nb,
(Fxc(D)
)
µν=
1
2
∫
IR3
µxc(ρ)χµχν
❲
µxc(ρ) = ∂Exc/∂ρ
ρ
P
❨❯ ❴❨ ❴❯ ❴❬
❲ ❲ ❴ ❲ ❴❨❬❴ P ❲ ❨ ❨❯ ❴ ❲ ❲ ❴
❴❨ ❯ ❨❯
![Page 24: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/24.jpg)
❨❯ ❯ Vx
❬ P
V ps ❴
∀ ψ(r, θ, ϕ) V ps.ψ(r, θ, ϕ) = V psloc(r)ψ(r, θ, ϕ)
+
lmax∑
l=0
l∑
m=−l
(∫
S2
ψYml dθdϕ
)V ps
l (r)Yml (θ, ϕ)
(Ym
l )lm ❴ ❴ ❭ P
V psloc
V ps
l rc
❨ ❴❯ P
φi
V psloc
V psl
[0, rc]
❭
hV hps
∀ 1 ≤ µ, ν ≤ Nb, hpsµν =
∫
IR3
χµ
(V ps.χν
)=(hps
loc
)µν
+
lmax∑
l=1
(hps
l
)µν.
❭ ❬ ❬ ❴ ❴ ❬ ❴
❫ P ❯ P
❯ ❲ ❯ P
❭ ❯ ❬ ❭ ❭ ❴
φi
❬ ❲❬P Nb = O(N)
❯ P ❴ ❲ ❲
S
![Page 25: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/25.jpg)
❭ ❴
(S = INb)
❨❯ ❯ ❴❨ ❭ ❲ P ❭ ❯ P
φi
❯ ❨
Nb ≃ 102N
❭ ❴ ❨❯
❭ ❨❯ ❲❴ P ❨❯ ❫ ❯ ❴ ❯
❭ ❲ ❵❨ ❨
(µν|κλ) ❯
❭ ❴
F (D)❲
J
ρ
2Jµν =
∫
IR3
(ρ ⋆
1
|x|)χµχν .
❭ J
❲ O(Nα
b
) ❲α ≤ 3
S = INb
hL
❯ ❴ ❴❨ ❴❬❴ ❴
hV J Fxc
❲
Vx
ρ ⋆ 1
|x|
µxc(ρ)
P ❴ ❲ O
(Nblog(Nb)
) ❭
F (D) O(N2
b
) S
Fxc
❭ hL P ❬
❴ W = Vx + ρ ⋆
1
|x|W
−W = 4π(ρx − ρ)
ρx
![Page 26: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/26.jpg)
❴ ❭ P O(Nblog(Nb)
) ❴ ❬
∫
IR3
Wχµχν = hVµν + 2Jµν(D).
F (D)
S
❴ ❴ ❨ ❴❬❴ ❭
F (D)
S O(Nblog(Nb)
) ❴ ❴❬ ❭❫ ❲ ❲
Nb
❭ Vx
P P❲ ❴ P P
φi ❯ ❯ ❨
❴ ❲ ❬ ❨ ❬ ❯ ❭
hpsloc
❯
Fxc
V loc
ps
J
P
ρ ⋆ 1|x|
❲
Vx
❬ ❭
hpsl
❭ P ❬❵❨ ❯ ❨ ❯
Nb(Nb + 1)/2
Nb
❭ ❯ ❬ ❲❴ P
hpsloc
hps
l
V psloc
V ps
l
❴
❯ V ps
loc
V ps
l hps P
χµ
❬ ❬ ❲
Nb
F (D)
❲❴
hpsl
O(N3b )
❫❲ ❲ O
(N2
b log(Nb))
F (D)
P
N ≪ Nb
❴ ❲ ❲ P O(Nblog(Nb))
![Page 27: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/27.jpg)
❭ P ❲ ❲ ❯ ❲ P ❴ Ed
❴ ❲ ❲
❫ ❯ Dn
Dn
Dn+1
F (Dn)
❫ ❴
Dn+1 ❴❴
P ❲ ❴ ❯ ❬
D⋆ ❬❯ ❲ F (D⋆)
ǫi
N ❲
F (D⋆)
Dn+1
Dn = CCt
C P N
❲ F (Dn)
❴
Dn+1 = D ❲ ❴ ❴
❴ ❫ ❴ ❲ ❲ ❴ ❴ ❲ ❴
Dn
F (Dn)
❭ ❬ ❯ ❴ ❴
❯ ❲ ❭ ❯ ❴ O(N3
b
) ❯
F (D)S
F (D)
S
O(N2
b ) ❭
χµ
Nb = O(N)
❯ O(N3)
N ≪ Nb
❴ ❲ O(NbN2)
O(N3)
❴
F (D) ❲
O(N2) O
(Nblog(Nb)
) ❯ ❴ ❴
![Page 28: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/28.jpg)
N
❭ P P ❨
❯
∂Ed
∂x=(∂h
∂xD
)+
1
2
(∂G∂x
(D)D
)− (
DE∂S
∂x
)
❲DE = CECt = S−1F (D)D
G(D) = 2J(D) −K(D)
❭ ∂h∂x
∂G∂x
∂S∂x
P ❲ ❭ ❬ ∂hps
∂x
P P
χµ
x∂G
∂x= ∂S
∂x= 0
N ❴
❲ ❯ ❲ Nb
inf
(FCCt), C ∈ M(Nb, N), CtSC = IN
.
❭ ❴
N
F
O(N3) P
❯ D⋆ = CsolC
tsol
Csol P
❭ N
F
❬ P ❲
inf
(FD), D ∈ MS(Nb), DSD = D,
(SD) = N,
D⋆ ❲ ❴
D⋆
D⋆ ❫
γ = (ǫN+1 − ǫN)/(ǫNb− ǫ1)
❲ ǫF = (ǫN + ǫN+1)/2
(ǫi)i ❲
F
❭ ❴
❴❨ ❨❯ ❲ ❫ ❲ P
![Page 29: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/29.jpg)
ψw
i
P
ψi
❴❨
|r − r′| → +∞ D
D(r, r′) = 2N∑
i=1
ψi(r)ψi(r′) = 2
N∑
i=1
ψwi (r)ψw
i (r′) = 2
Nb∑
µ,ν=1
D⋆µνχµ(r)χν(r
′).
❫ ❲ ❨❯ ❲
❭ ❬ γ 6= 0
D(r, r′) ≃ exp(− γα|r − r′|
).
α = 1/2
α = 1
❬ ❭ ❨❯
γ = 0
❨ ❨❯ ❬ ❲ ❭
ψwi
❨❯ γ = O(1)
N
❬ P ψw
i
❯ χµ
C⋆ D⋆ = C⋆(C⋆)t ❫
P ❯ D⋆
D
χµ
❯ ❯ O(N)
C⋆
D⋆ ❫ ❲ ❲ ❴ ❲
N
C⋆ D⋆ ❴
❯
C⋆ ❲ D⋆ ❲ P
C
D
❲
N
❫ ❴❨❴❯ ❯ C⋆
D⋆ ❲ ❯
C⋆
![Page 30: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/30.jpg)
N
❴ ❲❴
❴
FS
C
D
❬ ❴
FS ❲
N ❲ P ❲❬
❲ ❴ ❨❯ ❴ ❲
γ ❨❯
❴ ❲N
❫ ❲
S = INb
❨ ❴ ❯ ❨ ❴
❴ ❬ ❲ ❯ ❲ ❴ ❭
D⋆ ❬ ❴ ❲
❫ ❴ P D
P ❲ ❲ ❯ ❫ GC
❲ GD
C⋆ ❲ D⋆ ❲ GC
❲ GD
P
ψwi
❲ D
❭ ❴ ❲ P ❯ ❲ P
inf
ΩC(C), C ∈ M(Nb, N), inf
ΩD(D), D ∈ MS(Nb)
.
❵ ❲ ❲ PΩC
ΩD
P D
❴ ❫❳❳ P
C
❴ ❳ ❴ ❫❳❳
ΩD(D) =(
(F − ǫF INb)(3D2 − 2D3)
) ❴ ❳
ΩC(C) =(
(F − ǫF INb)C(2IN − CtC)Ct
)
D⋆ ΩD
❫
C⋆ ΩC
![Page 31: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/31.jpg)
❴ ❯
∣∣∣∣∣∣
inf
ΩC(C), C ∈ M(Nb, N), C ∈ GC
,
inf
ΩD(D), D ∈ MS(Nb), D ∈ GD
.
❴ ❬ ❲ ❴❴ ❴ ❲❬ ❲ GC,D
❭ ❴ ❬ ❭ P P ❨ ❲❯ ❨
❴ ❲ ❲
−∞ ❴ ❫❳❳ ❭S = INb
D0 = (1/2)INb
❲ ❲
❲ ❲ ❯ ❭ ❯ ❨❯ ❨ ❴ ❴❴ ❲❯ ❴ ❳
❴ ❳ ❲ ❲ GC
GC
❯ ❲ ❲
ΩC
❴ C
❴ ❲
(CCt) = N ❲ ❴
D
Tr(D) = N
❭ ❴
ρ ❲ ❯
❴ ❲ ❯ ❴
❯
❭ ❴ D
D = q(F ) ∀ q ∣∣∣∣ q(ǫi) = 1 ∀ 1 ≤ i ≤ Nq(ǫi) = 0 ∀ N + 1 ≤ i ≤ Nb
.
❴ ❲
![Page 32: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/32.jpg)
N
❭ ❴ q P
❨ ❴❨❴❲(Tn)n
❲ ❴ ❳
q P
pn ❲p : x 7−→ 3x2 − 2x3.❫ ❨ ❴❨❴❲❬❩ P
❴ Pq ❴ ❫❳❳
❴ D0
P
F
❭ ❴ ❴ P
D0 = F = αF + β
αβ
ǫ1ǫNb
ǫF
DnDn+1
Dn+1 = Dn + cnTn(F )
(cn)n
❴❬❨❴❲ P ❲
[−1, 1]
❳ Dn+1 = 3D2
n − 2D3n. ❴
Dn
❲ GD
❴ ❲
❴ ❲ ❲ ❬ ❳ P ❨ ❲ ❴ ❲ ❲ ❲ GD
❭ ❴ ❴ ❬
P ❨
D⋆ ❴ C⋆
❭ ❴ ❴ ❬P ❴ ❬ ❬ ❲❯ ❲ ❨ P ❲ GC,D ❴
C⋆ D⋆
❴ P
ǫF ❴ ❲ ǫF
❬
ǫ1ǫNb
❭ ❯ ❬ ❲ O(N)
❴ ❲ F
P ❴ ❲❴ǫF
![Page 33: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/33.jpg)
❴ ❴❲
ǫF
S 6= INb
❯ S = INb
❴❨
S−1 ❴ ❭ ❴
F
F = S−1F
❯ ❲
SF = F
F
F ❴
S−1 P
S−1 ❯
SD⋆ ❯
D⋆
❭ ❴ ❲ ❲ ❯ ❲❬ ❬
S−1 S−1F
P ❫❳❳ P ❫❳❳
S−1F P
❫❳❳ ❳ ❴❲
S 6= INb
❭ ❴ ❬ ❫ ❴ P ❭ ❴ P ❲
C❭ ❴ D
❲❯ P ❴
C P ❯ ❬
❲ P ❫ ❴ ❨❯ ❯
C⋆ D⋆
❬ ❴ P ❲
CCt
iCj = δij
D2 = D
❲ C⋆
D⋆ ❲
C⋆ ❲
![Page 34: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/34.jpg)
❵
C⋆
❭ ❴
❨❯
nit
❲ nit ≃ γ−α ❲
α❬
P ❯ ❨❯ γ ≃ 0
❫
C⋆ ❬ ❨❯
D⋆ ❴ ❲
N ❭ ❴
❨❯ P
❭ ❴ ❨❬
❯ ❨ ❨❯ ❨
10−15 ❲N
103
T
10−12 ❲ ❲
104 10−6 ❨ ❬ ❬
❭ ❲ ❴❯ ❲ ❬
❴
❭ ❬ ❬
❯ ❴❨ ❴ ❭ ❴ ❴❨❬❴ P ❬ ❲ ❲ ❨❯ ❯❲ ❴ ❨ ❬ ❴ ❲ ❬ ❭ ❴ ❫
[
] ❲ ❴
P ❴ ❬ ❳ ❨❯ ❴ ❴
![Page 35: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/35.jpg)
❬ ❴ ❬ ❴ ❲ P ❴
❭ ❯
❯ ❲ F
S ❴ ❨ ❭ ❲ ❲
❲ ❨❯ ❴ P P ❭ ❴ ❯ P ❲ ❲❬ ❴
❲ ❲ ❴ ❬
C ❨❯ ❴ ❬
❴ ❴
C⋆ ❲N
❯
C GC
inf
(FCCt), C ∈ M(Nb, N), CtSC = IN , C ∈ GC
.
❴ ❴ ❲ ❬ ❴
P
F ❲
S 6= INb
❬
❴ ❴ P ❴
❲ ❨ ❯
FS ❴ ❲ ❯
❲ ❫❳❳ ❯ ❴ P ❫❳❳ ❲ P ❲❯
![Page 36: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/36.jpg)
❵
❴ ❫❳❳
❴ ❫❳❳ ❴ ❲ ❲❴ ❴ ❨❯ ❴❨
C⋆ ❨❯
F ❴
P ❨ ❬
❴ [
]
❴ ❯ P ❴ ❲ ❲ ❴ ❯ ❵ ❲ P ❨ ❴ ❴ P ❴ ❲
❫ ❴ ❴
❨❯ ❲ ❲ ❪ ❫ ❨ ❵ ❫ P ❫ ❬ ❯
[
] [
] [
] [
] ❯
γ = O(1)
F
❲ ❴
F ❲ ❫
❴ ❳ ❫❳❳ P ❴
D⋆ij
D⋆
Fij 6= 0 ❲
N
❭ P ❭ ❴ P ❴ ❫❳❳
D0 ❨
F
❨ ❲
ux
Cx
Dx
❲ Pφi
❲ ❯ ❯
(Px)
❯ ❯ ❨ x ❴
![Page 37: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/37.jpg)
x❲
P ❴❨ ux
❲
ux
❯ P ❫ ❴ ❴ ❴ ❯
ux
❬ Ps(ux) ❨
❭ ❴ ❲ ❯
(Px)m
Pxk
R
uxk
B P P
uxk
❲ ❯
xl
❴uxl
uxl
P
uxl=
m∑
k=1
αkuxk.
uxl
❯ (Px)
m
αk(Px)
(Px)
R
xk
(Px)
(Px)
m P ❯
(Px) ❯
xl
❨ ❴ P
xk
❬ ❯ ❬ ❴ ❴
uxl
B ❴
uxl
❭ ❲
❭ ❯ B ❭ ❲ ❯ ❴ P ❲ ❭ ❨ ❲ ❨❯ ❨❯ B
m ❲❴
❴ ux
❲ ❭
![Page 38: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/38.jpg)
❵
❴ ❨ B ❫ ❯ B
❭ ❴ ❲
❯ (Px)
❴ ❴ ❴
uxux
❲ ❴ ❴ ❯ ❯
(Px)❲
x P
(Px)
−u+ xu = f
Ω,
u = 0
∂Ω.
❭ ❯ ❴z ❭ ❬
❬❭ ❯ ❲ ❯ ❬
−1
2φi + Vxφi +
(ρ ⋆
1
|x|)φi −
N∑
j=1
(φiφj ⋆
1
|x|)φj = ǫiφi i = 1, N,
ρ(x) = 2N∑
i=1
φ2i (x), Vx =
z
|x− x| ,∫
IR3
φiφj = δij .
P
φi
❯
x
❯ ❲ N
Pφi)
N(N + 1)/2
❯ ❴ ❴ ❴ ❴ ❯ P
[
] ❯ ❴
❭ ❯ ❯ ❴ ❲ ❯ ❲ ❯
❴φi
Bi
P m
P(φi)µi
k
∀ 1 ≤ i ≤ N, φi =
mi∑
k=1
αi,k(φi)µik.
![Page 39: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/39.jpg)
❴ ❴ ❬ O(N)
❯ ❯ O(N)
αi, k)
(φk
1, . . . , φkN)k
P m N
❬(φk
i )i=1,N
❴ ❴❴
φi
P
∀ 1 ≤ i ≤ N, φi =m∑
k=1
αkφki .
αk
❴ ❴φi
❬❴
N ❫
m
❬
mi
❴ ❨❯
❴ ❯ ❯
N❲ ❯ P
P P ❲ ❲ ❴ux
❴ ux ❫ ❴ ❲❬
❴ ❯ ❲ ❳ P ❴ P ❲
N ❲
φi
❨xk
❭ ❴
φi
❭
χi
❨ ❯ ❴ P Tx
φi
ψi = Txφi ❨ ❯
ψi
❴ ❴
ψi
❭ χi
❨ ❴
CxDx
❴ ❲ P Tx
❯ ❯
IR3 ❲❴ P ❴ ❯ ❴ ❲ ❬ P
φi
❨ P❬ Tx
❲
![Page 40: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/40.jpg)
❵
❴ ❯ ❯ ❨❯ ❨❯ ❴
H+2
H2
❯
N = 1 ❭ ❯
µ ❨ ❴❨❯
P ❴❨❴❯ ❨ ❨ Pφµ
χi
P1
❭ Tµ
❯ ❯ P
ψµ = Tµφµ
❴
❫ ❯ ❴ ❴ ❫ ❨❯ ❲ P ❴ ❨❯
![Page 41: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/41.jpg)
![Page 42: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/42.jpg)
❭ ❬ ❴❨ ❴ ❴ ❭ ❲ ❴ ❴❯ ❴ ❲ ❲ ❴❨❬❴ P ❲ ❨ ❨❯ ❴ ❲ ❲ ❴
❫ ❴ ❬ ❲ P ❴❨ ❯ ❨❯ ❬ ❨ ❲
❭ ❯ ❴ ❴ ❴ ❨❯ ❴ ❨❯ P ❭ ❬ ❲ ❲ ❨❯ ❬ ❬ P ❲ ❲ ❲ ❴❨❬❴ ❬❲ ❴ ❬ ❬ P
![Page 43: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/43.jpg)
❲ ❲ ❴ ❴ ❭ ❬ ❲
❵ ❴ ❯ ❴❨❬❴ ❬ ❲ ❨❯ ❬ ❲ ❬ P ❲ ❲ P ❬ ❭ ❬ ❲❬P ❴ ❴ ❴ ❴ ❴❨ ❯ ❯ ❨❯ ❴❨ ❴❴
❫ ❯ ❲ ❫ ❯ ❴ ❵ ❴ ❴ ❯ ❴❲ ❴
❬ ❲ ❭ ❯ ❲ ❫ ❬ ❯ ❲ ❬ ❲ ❴ ❬ ❲ ❴ ❯ ❲ ❲ P ❲ ❬ ❫
(Yml )l≥0,−l≤m≤l
❴ ❴
❴ ❬
❨❯ ❭ ❲ P ❬ ❨❯
❯ ❨ ❴
z
x ∈ IR3 N
Vex
❭ P
![Page 44: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/44.jpg)
❨❯ ❲
(P)
∥∥∥∥∥∥∥
❭(ψi, ǫi)i
N
∀ 1 ≤ i ≤ N, −1
2ψi +
(Vnuc + Vex
)ψi + KRHF,KS(Ψ).ψi = ǫiψi,
∀ 1 ≤ i, j ≤ N,(ψi, ψj
)= δij ,
❲ǫ1 ≤ . . . ≤ ǫi ≤ . . . ≤ ǫN .❭
Vnuc
∀ x ∈ IR3 Vnuc(x) = − z
|x− x| , Ψ
P P(ψi)i=1,N
∀ ψ ∈ L2(IR3),(KRHF (Ψ).ψ
)(x) = J(ρΨ)(x)ψ(x) − 1
2
(K(τΨ).ψ
)(x),
∀ ψ ∈ L2(IR3),(KKS(Ψ).ψ
)(x) = J(ρΨ)(x)ψ(x) + µxc(ρΨ)(x)ψ(x),
❲
∀ x ∈ IR3, J(ρΨ)(x) =
∫
IR3
ρΨ(y)
| x− y | dy,
∀ x ∈ IR3,(K(τΨ).ψ
)(x) =
∫
IR3
τΨ(x, y)
|x− y| ψ(y) dy,
∀ x ∈ IR3, ρΨ(x) =N∑
i=1
|ψi(x)|2,
∀ x, y ∈ IR3, τΨ(x, y) =
N∑
i=1
ψi(x)ψ∗i (y).
µxc(ρ)
❴❬ ❬
ρ
Nc ❲
Nv ❬
❲ ❲ ❲
1 ≤ i ≤ Nc, ψci = ψi,
1 ≤ i ≤ Nv, ψvi = ψi+Nc
, ❲
1 ≤ i ≤ Nc, ǫci = ǫi
1 ≤ i ≤ Nv, ǫvi = ǫi+Nc
❲
❲ ❨❯ ❫ K KRHF,KS P P KKS KRHF P
![Page 45: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/45.jpg)
❯ P
(Pr) ❬
❴ P ❭ ❬❭
(Pr) ❯ ❲
(Pr)
∥∥∥∥∥∥∥
❭(ψr
i , ǫri )i
N
∀ 1 ≤ i ≤ N, −1
2ψr
i + Vnucψri + K(Ψr).ψ
ri = ǫriψ
ri ,
∀ 1 ≤ i, j ≤ N,(ψr
i , ψrj
)= δij .
❲Ψr = (ψr
i )1≤i≤N
ǫr1 ≤ . . . ≤ ǫri ≤ . . . ≤ ǫrN .❫ P
∀ 1 ≤ i ≤ Nc, ψr,ci = ψr
i ,∀ 1 ≤ i ≤ Nv, ψr,v
i = ψri+Nc
,∀ 1 ≤ i ≤ Nc, ǫr,ci = ǫri ,∀ 1 ≤ i ≤ Nv, ǫr,vi = ǫri+Nc
.
❭ ❬
∀ 1 ≤ i ≤ Nc, ∀ x ∈ IR, ψc
i (x) ≃ ψr,ci (x),
∀ 1 ≤ i ≤ Nc, ǫci ≃ ǫr,ci ,
∀ 1 ≤ i ≤ Nv, ∀ x ∈ Sc, ψv
i (x) ≃ ψr,vi (x),
Sc
ψr,ci ❴❯ ❨
rc ❨
❭ ❯
(P)
❭(ψi, ǫi)i
N
(2.1a)
∀ 1 ≤ i ≤ Nc, −1
2ψc
i +(Vnuc + Vex
)ψc
i + K(Ψc).ψci + K(Ψv).ψc
i = ǫciψci , (2.1b)
∀ 1 ≤ i ≤ Nv, −1
2ψv
i +(Vnuc + Vex
)ψv
i + K(Ψc).ψvi + K(Ψv).ψv
i = ǫviψvi , (2.1c)
∀ 1 ≤ i, j ≤ Nc,(ψc
i , ψcj
)= δij, (2.1d)
∀ 1 ≤ i, j ≤ Nv,(ψv
i , ψvj
)= δij , (2.1e)
∀ 1 ≤ i ≤ Nv, ∀ 1 ≤ j ≤ Nc, (ψvi , ψ
cj) = 0. (2.1f)
![Page 46: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/46.jpg)
❲Ψc = (ψc
i )1≤i≤Nc
Ψv = (ψv
i )1≤i≤Nv
❭
(2.1a) − (2.1f) K P
N
❲ KRHF ❲❴ ❲ KKS P ❴❨❴❯ ❲ ❭ ❴❨❴❯ ❲❯ P
(ψvi , ǫ
vi )i
❯ ❨❯
❭(ψv
i , ǫvi )i
Nv
∀ 1 ≤ i ≤ Nv, −1
2ψv
i +(Vnuc + Vex
)ψv
i + K(Ψc).ψvi + K(Ψv).ψv
i = ǫviψvi ,
∀ 1 ≤ i, j ≤ Nv,(ψv
i , ψvj
)= δij ,
∀ 1 ≤ i ≤ Nv, ∀ 1 ≤ j ≤ Nc, (ψvi , ψ
cj) = 0,
❲ǫv1 ≤ . . . ≤ ǫvi ≤ . . . ≤ ǫvNv
.
❨❯ ❬
Nv ❨
Vex
V tot =
Vnuc + K(Ψc) ❫ ❴
(2.1f)
❲ ❲ ❯ ❨ ❬ ❬
V ps ❴ ❯ (P)
N =
Nc +Nv
❯ (Pps)
Nv
(Pps)
∥∥∥∥∥∥∥
❭(φps
i , ǫpsi )i
Nv
∀ 1 ≤ i ≤ Nv, −1
2φps
i +(Vnuc + Vex
)φps
i + V psφpsi + K(Φps).φps
i = ǫpsi φ
psi ,
∀ 1 ≤ i, j ≤ Nv,(φps
i , φpsj
)= δij.
❭ (Pps)
(P)
Nv < N ❬
φpsi
❲ ❯
ψvi
❲ ❲ ❲ ❴❨❬❴ ❲ P ❴ ❯❴ ❴❴
V ps
∀ 1 ≤ i ≤ Nv, ∀ x /∈ Sc φpsi (x) = ψv
i (x),
∀ 1 ≤ i ≤ Nv, ǫpsi = ǫvi .
❫ ❴❴ ❬
φpsi
Sc
(Pps)
![Page 47: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/47.jpg)
❭ ❬ V ps
ψi
Sc ❲ ❲ ❲ ❬
V ps = V psr
V psr
❬ ❯
V ps ❬ ❯ P(Pr)
P P
(Pr) ❯
(Ppsr )
(Ppsr )
∥∥∥∥∥∥∥
❭(φr,ps
i , ǫr,psi )i
Nv
∀ 1 ≤ i ≤ Nv, −1
2φr,ps
i + Vnucφr,psi + V ps
r φr,psi + K(Φps
r ).φr,psi = ǫr,ps
i φr,psi ,
∀ 1 ≤ i, j ≤ Nv,(φr,ps
i , φr,psj
)= δij,❲
ǫr,ps1 ≤ . . . ≤ ǫr,ps
j ≤ . . . ≤ ǫr,psNv,
∀ 1 ≤ i ≤ Nv, ∀ x /∈ Sc, φr,psi (x) = ψr,v
i (x),
∀ 1 ≤ i ≤ Nv, ǫr,psi = ǫr,vi ,
❲ ❬
φr,psi
Sc
❬ V ps ❬
V ps
r
❨❯ P P ❬ ❬V ps
r
❫ ❯ ❲ ❬ ❴
(2.1f r)
(2.1f) ❯
(Pr)
∀ 1 ≤ i ≤ Nv, ∀ 1 ≤ j ≤ Nc,(ψr,v
i , ψr,cj
)= 0. (2.1f r)
ψr,c
i
(Pps
r )❲
V psr = K(Ψc
r) ❴
(2.1f r) ❴
V psr
P P (φr,ps
i , ǫr,psi )
❲ ❲ ❲ −1
2 + Vnuc + V ps
r + K(Φpsr )
❭ ❯ ❴ P❬ ❲
ψr,ci
❭ ❬ φr,ps
i (Pps
r )
ψr,vi
IR3 ❭
(2.1f r)
❲ ❬ φr,ps
i
![Page 48: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/48.jpg)
❫ ❴ P
ψri
P −1
2+ Vnuc + V ps
r + K(Ψvr)
❴❴ ❬
φr,psi
(Pps
r )
ψr,vi
❴ Sc Sc
ψr,vi
❭❬
φr,psi
Sc ❲ ψr,v
i
(2.1f r) ❲ ❬
φr,psi
❬
V psr
ψv
i
Sc
ψr,vi
Sc
(2.1f) ❨❯
❲ ❬V ps ❵ ❴
❲ ❴
❭ ❳❯ ❯ ❬
❯ ❴ ❭ P ❯ ❴❬❴ ❴❨
❫ ❴ K(Ψcr)
V pm
r
P ❲
V pmr =
Nc
r−
na∑
j=1
Aje−αjr2
r.
Vnuc + V pm
r
V add
r
V addr =
Nc∑
i=1
ǫ|ψr,ci 〉〈ψr,c
i |,❲
ǫ
ǫ+ ǫr,c1 ≫ ǫr,vNv.
❯ ❲ ǫr,ci
❴
(2.1f r) ❭ ❯
(Aj, αj)j=1,na
na
P ❲
ψr,vi
❲ǫr,vi
(Pps
r )❲
V psr = V pm
r + V addr .
(11)
(12)
❲ ❬ φr,ps
i
ψr,v
i
❬ ❯
(Ppsr )
❲ P
![Page 49: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/49.jpg)
❭ ❬
φr,psi
❴ ❳❬❯ ❯
ψr,vi
Sc ❬
V psr
P ❬ φr,ps
i
(Pps
r ) Sc ❲ ψr,v
i
❴ Sc
❭ ❯
(Pr) ❨❯ ❨ ❴ ❬
❴❴ V ps
r
P ❨
V psr (r, r′) = Vloc
(|r|)δ(r − r′)
+
lmax∑
l=0
l∑
m=−l
∣∣Yml
⟩(Vl(∣∣r|)− Vloc
(|r|))δ(|r| − |r′|
)⟨Ym
l
∣∣,
❭ lmax
l0 +1
l0
❲ ❲ ❭ P ❯ ❲ ❴ ❫ ❲ ❴ ❪ ❲ ❯ ❴ P
(13)
❬ φr,ps
i
❴ ❬ ❨
ECP❭ ❯ ❴
V psr
❭ ❯❬ P ❲
ψr,vi
❴ Sc ❲
ǫr,vi
(Pps
r ) ❴ ❴
❳❯ ❬ V ps
r
P ❲ P ❲ P
❭ ❴ ❬ φr,ps
i ❭ ❯ P ❲ ψr,v
i❴ Sc ❬ φr,ps
i
Sc ❬
V psr
❲ ❯ ❬
φr,psi
❲
ǫr,vi
❲ ❬ ❲
(φr,psi , ǫr,vi )
(Ppsr )
❲ (Pps
r )
![Page 50: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/50.jpg)
❴ ❴ ❬ φr,ps
i
P ❲ ψr,v
i
❴ Sc ECP
❨ ❲ ❴ ❭ P ❬
V psr
❬ P ❲ ❬ ❭ ❬
V psr
❴ ❬ ❯ ❬
ECP❴ ❭
ECP❴ ❲
❯ ❬ ❯ ❫ ❭ ❴❬ ❲ ❯ ❫ ❫ ❬ ❴ V ps
r
❫ ❴ ❬
ECP
(i)
ψr
i
❲
(Pr)
(ii) ❴ P
Vloc
Vl
0 ≤ l ≤ lmax
(iii)
❯ P (φr,ps
i , ǫr,psi ) ❯
(Ppsr )
❲V ps
r
(|r|, |r′|
)= Vloc
(|r|)δ(r − r′)
+lmax∑
l=0
l∑
m=−l
∣∣Yml
⟩(Vl
(|r|)− Vloc
(|r|))δ(|r| − |r′|
)⟨Ym
l
∣∣.❲
φr,psi = ψr,v
i
,
φr,psi
,
ǫr,psi = ǫr,vi .
❭ (iii)
❨ (iv)
(Pps)
❲V ps
V ps
r
❲ P
(Ppsr )
❯ (Pps
r )
ECP
(i)
ψr
i
❲
(Pr)
![Page 51: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/51.jpg)
(ii) ❴❴ ❴
φr,ps
i
❴ ψr,v
i
,
φr,psi
.
(iii) ❲ ❯ ❲
V psr
❬❬ (a)
❴l
i(l) ❲
l
n
❲ (b)
0 ≤ l ≤ lmax
Vl
(|r|)
=ǫr,vi(l)φ
r,psi(l) + 1
2φr,ps
i(l) − Vnucφr,psi(l) −K(Φps
r ).φr,psi(l)
φr,psi(l)
,
K KKS KRHF ❭ P |r|
r ❯
❬ φr,v
i(l)
P Vl
(|r|)
(c)
Vloc = Vlmax.
(d)
V psr
(|r|, |r′|
)= Vloc
(|r|)δ(r − r′)
+lmax∑
l=0
l∑
m=−l
∣∣Yml
⟩(Vl
(|r|)− Vloc
(|r|))δ(|r| − |r′|
)⟨Ym
l
∣∣.
(iv)
(Pps)❲
V ps
V psr
❲ P
(Ppsr )
(Pps
r ) ❲
(φr,ps
i , ǫr,vi )
(φr,psi , ǫr,ps
i )
(Pps
r )
❭ ❬ ❬
V psr
❯ ❯
(Pps)❯
(χµ)i=1,Nb
Nb(Nb + 1)/2 P
∀ 1 ≤ µ ≤ ν ≤ Nb
∫
r
(∫
θ,ϕ
Yml χµ
)( ∫
θ,ϕ
Yml χν
)(Vl − Vloc
)(r) dr.
❲ ❲❯ ❯ ❬❵❨
![Page 52: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/52.jpg)
P ❬ ❬ V ps
r
lmax∑
l=0
l∑
m=−l
∣∣φr,pslm (Vl − Vloc)
⟩⟨φr,ps
lm
(Vl − Vloc
)∣∣⟨φr,ps
lm
∣∣Vl − Vloc
∣∣φr,pslm
⟩
❴ Pφr,ps
lm
Pφr,ps
i Yml
φr,ps
i
lm
Nb
P
∀ 1 ≤ µ ≤ Nb,
∫
IR3
φr,pslm (Vl − Vloc)χµ.
❭ ❯ P ❬ ❲ ❬ P ❲ P
φr,psi
❬❲ ❴
Vloc ❬
Vl
❲
Vlmax
❨ rc
❯ ❲ ❬ ❲
❯ ❨❯
M(Ak)k=1,M
(zk)1≤k≤M
N
❬ ❨❯
(xk)k=1,M
❨ ❭ P ❬ ❨❯ ❲
(P)
∥∥∥∥∥∥∥
❭(ψj , ǫj)j
N
∀ 1 ≤ j ≤ N, −1
2ψj + Vnucψj + KRHF,KS(Ψ).ψj = ǫjψj ,
∀ 1 ≤ i, j ≤ N, (ψi, ψj) = δij ,
Ψ
P P(ψj)1≤j≤N
ǫ1 ≤ . . . ≤ ǫj ≤ . . . ≤ ǫN .
❭ Vnuc
∀ x ∈ IR3, Vnuc(x) =
M∑
k=1
Vk(x) = −M∑
k=1
zk
|x− xk|,
∀ x ∈ IR3, ρΨ(x) =N∑
j=1
|ψj(x)|2,
![Page 53: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/53.jpg)
∀ x, y ∈ IR3, τΨ(x, y) =
N∑
j=1
ψj(x)ψ∗j (y).
❫ ψj
❨ ❬P ❨❯ ❲
ψj
❬ ❲ P ❯ ❨❯ P
µ
ψr1s
1s
rc
ψr
1s
❲ µ
2rc
ψ1
ψ2
❯
(P)
0 0−µ/2
µ/2
−µ/2 µ/2
ψ1 ψ2
Pψ1
ψ2
❲
x ∈ IR3
ψ1(x) ≃(ψr
1s(x+ µ/2) + ψr1s(x− µ/2)
)/√
2
ψ2(x) ≃(ψr
1s(x+ µ/2) − ψr1s(x− µ/2)
)/√
2.
φ1 = (ψ1 + ψ2)/
√2
φ2 = (ψ1 − ψ2)/
√2
P
x ∈ IR3 φ1(x) ≃ ψr
1s(x+ µ/2)φ2(x) ≃ ψr
1s(x− µ/2).
❯ ❨❯ ❴
N c
k
❨
rc,k
❯ P
(Pr,k)
Ak
❬
V psr,k
![Page 54: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/54.jpg)
00 −µ/2−µ/2 µ/2
φ1 φ2
µ/2
Pφ1
φ2
❴ ❯ (P)
N
❯ (Pps)
Nv =
N −M∑
k=1
N ck
(Pps)
∥∥∥∥∥∥∥
❭(φps
j , ǫpsj )j
Nv
∀ 1 ≤ j ≤ Nv, −1
2φps
j + Vnucφpsj + V psφps
j + K(Φps).φpsj = ǫps
j φpsj ,
∀ 1 ≤ i, j ≤ Nv, (φpsi , φ
psj ) = δij ,
❲ǫps1 ≤ . . . ≤ ǫps
j ≤ . . . ≤ ǫpsNv
V ps =M∑
k=1
V psr,k.
❭ P ❴ Ak
❴
∀ 1 ≤ j ≤ Nv, ∀ x /∈ Sc φpsj (x) ≃ ψv
j+Nc(x),
∀ 1 ≤ j ≤ Nv, ǫpsj ≃ ǫvj+Nc
❲ Sc =
⋃Mk=1 Sc
k
❭ P ❴
❲ ❲ ❬ ❬❲ ❬ ❴ P ❯ ❲ P ❬ ❨❯ P ❲ P ❬
![Page 55: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/55.jpg)
❴ ❬ ❬ ❲ ❲ ❴ Ak
❲ ❲ ❲
❲ ❲
P
❭ (i− iii)
❯ P ❲ ❲ ❴ P
(i− iii) ❴
P P ❲ ❭ ❴ ❬ ❴ ❲ ❲ ❲ P ❬ ❨❯ ❲ ❲ ❴ ❴❨❬❴ ❨❯
SC
❬ ❲ ❴ ❫ ❯
LC ❬
❲ ❴ ❬ ❭ ❬
SC
P ❲ ❲ ❴ ❬ O(N7)
LC
SC
❲
❯ ❬
LC
(i − iii) ❲❴
❲ 3d
4s
f
5f
6s
❯ ❨❯ ❯ ❲ ❴
❯ ❴ ❲
(Pps)
P ❲ ❯ ❲ Pps P V ps ❬
φr,psik
❴ ❬
![Page 56: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/56.jpg)
V ps ❭
❬
V ps
❭ ❴ ❴ ❯ ❯ ❬ ❫ ❯
❯ ❬ ❨
rkc
P
Ak
❬
φpsik
❨ ❬ ❯
❲❴ P ❴ O(N7)
❲ ❬SC
❲ ❬ ❴ ❯ ❴ ❲ ❨
LC ❯ ❲
P ❲ ❬ ❲ ❴
LC
❭ ❴❨ ❴❨
❨❯ ❲
❭ ❬ ❬ ❯ ❲ ❯ ❬ ❨ P ❲❴ ❴ ❨ ❴❴ ❯ ❴ ❴ ❨❬❯
❬ ❯ ❬ P ❲ P ❬ ❴ ❲ ❲
![Page 57: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/57.jpg)
❭ ❬ ❯ ❯ ❬ P ❯ ❯ ❨ ❴ ❴ ❭ ❴ P P ❴ ❬ ❲ ❵ ❴ ❴ ❬ ❴ ❬❴ ❲ ❴ ❫ ❴❨❴❯ ❬ ❴ ❲ P ❴ ❴
❫ ❯ ❴❬❴ P
Wρ(x) = −M∑
k=1
zk
|x− xk|+
(ρ ⋆
1
|x|
)(x) + µxc[ρ](x)
1/|x|
xk ❨ ❬
No ≥ N
❴❬❴ ψi1≤i≤No
Hρψi = ǫiψi ǫ1 < ǫ2 ≤ ǫ3 ≤ · · ·
(ψi, ψj)L2 = δij
ρ(x) =
+∞∑
i=1
fi|ψi(x)|2
∣∣∣∣∣∣
fi = 1ǫi < εF
0 ≤ fi ≤ 1ǫi = εF
fi = 0ǫi > εF
+∞∑
i=1
fi = N
No
fi
N
No
N −1+nεF
nεF
εF
No
Hρ
0
![Page 58: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/58.jpg)
Hρ = −1
2∆ +W
xk
❭ ❴ ❲
T
L2(IR3)
I + T ❲ ψi1≤i≤N0
ψi
1≤i≤N0
Hρψi = ǫiSψi ǫ1 < ǫ2 ≤ ǫ3 ≤ · · ·
(ψi, Sψj)L2 = δij
ρ(x) =
N∑
i=1
fi
([(I + T )ψi
](x))2
∣∣∣∣∣∣
fi = 1ǫi < εF
0 ≤ fi ≤ 1ǫi = εF
fi = 0ǫi > εF
+∞∑
i=1
fi = N
Hρ = (I + T T )Hρ(I + T )
S = (I + T T )(I + T )
P
ψi = (I + T )ψi.
P ❲ ❬ ❭ ❴
T
(I + T )
❲ ❯ ❯
ψi❯ P
❴ P ❯ P ❲
❴ T
❴❬❴ ❲ ❲ ❯ ❫
Tz
T
P z
❨
Tz
❴
∀u ∈ L2(IR3), ∀x ∈ IR3 \Brzc(0), (Tzu)(x) = 0,
rzc
❨ P ❬
Br(x)
IR3 x
![Page 59: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/59.jpg)
❨r
(I+Tz) P
Brz
c(0)
P
Tz
❴
❨❯ M
❨ ❴z1, · · · , zM
x1, · · · , xM
T P
T =
M∑
k=1
τxkTzk
τxkTzk
Tzk
xk ❬❬
L2(IR3)
∀ u ∈ L2(IR3), ∀ x ∈ IR3,((τxkTzk
)u)(x) =
(Tzk
(u(· + xk)
))(x− xk).
1 ≤ k < l ≤M, Brzkc
(xk) ∩ Brzlc(xl) = ∅
❬❲ ❲ P ❲ P
T
❯ ❲ ❴ Tzk
❲❴ ❲ ❲ ❨❯ ❴ ❯
Tz
❯ P ❲❬❬❲
Tz
N = z
Tz
P❲
rzc
❨ ❴❯ ❴❯ ❲
φz,µ1≤µ≤Np
Np
P L2(IR3)
1 ≤ µ ≤ Np
φz,µ(x) =
µ∑
ν=1
αµνφ0z,ν(x)
![Page 60: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/60.jpg)
αµν
φ0
z,ν
1≤ν≤Np
Np
❴❬❴
Hzρφ
0z,ν = ǫ0z,νφ
0z,ν ǫ0z,1 < ǫ0z,2 ≤ ǫ0z,3 ≤ · · ·
(φ0z,µ, φ
0z,ν)L2 = δµν
ρ(x) =
+∞∑
ν=1
fν |φ0z,ν(x)|2
∣∣∣∣∣∣
fν = 1ǫ0z,ν < ε0
z,F
0 ≤ fν ≤ 1ǫ0z,ν = ε0
z,F
fν = 0ǫ0z,ν > ε0
z,F
+∞∑
ν=1
fν = N
❲Hz
ρ = −1
2∆ − z
|x| +
(ρ ⋆
1
|x|
)+ µxc[ρ];
φz,ν
1≤ν≤Np
P L2(IR3)
❯ ❲❬
∀ 1 ≤ ν ≤ Np, ∀ x ∈ IR3 \Brzc(0), φz,ν(x) = φz,µ(x);
pz,ν1≤ν≤Np
P L2(IR3)
φz,ν
1≤ν≤Np
❲
1 ≤ µ, ν ≤ Np, (pz,µ, φz,ν)L2 = δµν ;
Tz
❲ ❲ ❬
∀ u ∈ L2(IR3), Tzu =
Np∑
ν=1
(pz,ν , u)L2 (φz,ν − φz,ν).
1 ≤ ν ≤ Np
(I + Tz)φz,ν = φz,ν.
![Page 61: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/61.jpg)
ψi
1≤i≤N
Hzρ ψi = ǫiSψi ǫ1 < ǫ2 ≤ ǫ3 ≤ · · ·
(ψi, Sψj)L2 = δij
ρ(x) =
+∞∑
i=1
fi
([(I + Tz)ψi
](x))2
,
∣∣∣∣∣∣
fi = 1ǫi < εz,F
0 ≤ fi ≤ 1ǫi = εz,F
fi = 0ǫi > εz,F
+∞∑
i=1
fi = N
❲Hz
ρ = (I + T Tz )Hz
ρ(I + Tz)
S = (I + T Tz )(I + Tz),
P φz,ν
1≤ν≤N
ψi(x) =
i∑
ν=1
βiνφz,ν(x)
∀ 1 ≤ i ≤ ∞, ǫi = ǫz,i. P
[β] ❲ ❬
P[α]
❭ ❯ P
ψi
1≤i≤N
❯ P ❯ ❴ φz,ν
1≤ν≤N
❭ ❴
Tz
P ❲ ❴ ❲❬ ❯ P P ❬
Tz
P P
T ❨❯ ❲❬
P ❬❲ ❲
![Page 62: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/62.jpg)
1❯
φ0z,ν .
❯ ❴❬❴ ❨ ❴ ❫ P
φ0z,ν(x) =
φ0z,nνlν
(r)
rYmν
lν(θ, ϕ)
❨ ❴ ❬
k❬❯
ρk
❨ ❴
Hzρk
❨ ❴ P❬ P
φk+1n,l,m(x) =
φk+1n,l (r)
rYm
l (θ, ϕ).
(n, l)
fn,l,m
❬
m
ρk+1(x) =∑
n,l,m
fnlm|φk+1n,l,m(x)|2
=∑
n,l
fnl
∣∣∣∣∣φk+1
n,l (r)
r
∣∣∣∣∣
2 ∑
−l≤m≤l
|Yml (θ, ϕ)|2
=∑
n,l
(2l + 1) fnl
|φk+1n,l (r)|24πr2
❨ ❴ ❲
ρSCF (x) =+∞∑
ν=1
wnν lν
|φ0z,nνlν (r)|24πr2
.
2❯
φ0z,ν.
❴ P
rzc
❨ ❴❯ P ❯
k : IR+ → IR+ C1 ❲
k(0) = 1k′(0) = 0
,
k(rz
c ) = 0k′(rz
c ) = 0,
∀ r ≥ rzc , k(r) = 0.
❭ Pk
∀ 0 ≤ r ≤ rzc , k(r) =
[sin(π r/rz
c )
(π r/rzc )
]2
;
![Page 63: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/63.jpg)
V0
❴❴ P
φ0z,ν
1≤ν≤Np
P
φ0z,ν(x) =
φ0z,nνlν
(r)
rYmν
lν(θ, ϕ)
❯
(HPS
ρ − ǫ0z,ν
)φ0
z,ν = Cν k(r) φ0z,ν
φ0z,ν(r
zc ) = φ0
z,ν(rzc )
∂φ0z,ν
∂r(rz
c ) =∂φ0
z,ν
∂r(rz
c )
ρ(x) = ρ(r) =
+∞∑
ν=1
wnν lν
|φ0z,nνlν
(r)|24πr2
.
❲HPS
ρ = −1
2∆ + veff
ρ
veff
ρ = vloc +
((ρ+ ρρ) ⋆
1
|x|
)+ µxc(ρ)
vloc(x) = V0 k(|x|)ρρ(x) = Q00
ρ g00(x)
g00(x) =k(|x|)
4π∫ +∞
0r2 k(r) dr
( ∫IR3
g00 =
∫
Brzc(0)
g00 = 1)
Q00
ρ = −z +
∫
Brzc(0)
(ρSCF − ρ).
wnν lν
❲ ❲ ❯ ❴❬❴
φ0z,ν
P
k ❴
Brzc(0)
❴❯
∀ x ∈ IR3 \Brzc(0),
(( ρ|Brz
c(0) + ρρ) ⋆
1
|x|
)(x) =
−z +∫
Brzc(0)ρSCF
|x|
= − z
|x| +
(ρSCF |Brz
c(0) ⋆
1
|x|
)(x),
![Page 64: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/64.jpg)
❲ IR3 \Brz
c(0)
(−1
2∆ − z
|x| +
((ρSCF |Brz
c(0) + ρ|IR3\Brz
c(0)
)⋆
1
|x|
)(x) + µxc[ρ](x)
)φ0
z,ν = ǫ0z,νφ0z,ν
φ0z,ν(r
zc ) = φ0
z,ν(rzc )
∂φ0z,ν
∂r(rz
c ) =∂φ0
z,ν
∂r(rz
c )
ρ(x) = ρ(r) =+∞∑
ν=1
wnν lν
|φ0z,nνlν
(r)|24πr2
.
P ❴❬ ❨ ❭❫ ❲
φ0z,ν
IR3 \
Brzc(0)
P ❲
Brzc(0)
❯ ❴ ❲
uν(r) = rφz,nνlν(r) ❯
−1
2u′′ν(r) +
lν(lν + 1)
r2uν(r) + veff
ρ (r)uν(r) − ǫ0z,νuν(r) = Cν k(r) uν(r)
uν(0) = 0
uν(rzc ) = rz
cφ0z,nνlν
(rzc )
u′ν(rzc ) = rz
c
dφ0z,nνlν
dr(rz
c ) + φ0z,nνlν (r
zc )
ρ(r) =1
4π
+∞∑
ν=1
wnν lν |uν(r)|2.
−1
2u′′(r) +W (r)u(r) = Cν k(r) u(r)
u(0) = 0
u(rzc ) = rz
cφ0z,nν lν
(rzc )
u′(rzc ) = rz
c
dφ0z,nνlν
dr(rz
c ) + φ0z,nνlν (r
zc )
![Page 65: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/65.jpg)
P u(r) = rz
cφ0z,nν lν(r
zc ) v(r)/v(r
zc )
v ❯ ❲
−12v′′(r) +W (r)v(r) = λ k(r) v(r)
v(0) = 0
v′(rzc ) =
(1
φ0z,nνlν
(rzc )
dφ0z,nνlν
dr(rz
c ) +1
rzc
)v(rz
c )
❯ n ∈ IN
❨ n
]0, rz
c [ ❭ ❲
❴ uν
❴❲
l
ǫ0z,ν
Pφ0
z,ν
lν = l
❲ Puν
nν
❬❯ ❲ǫ0z,ν
n− 1
]0, rzc [
3
❯ p0
z,ν .
p0z,ν(x) =
k(|x|) φ0z,ν(x)
〈φ0z,ν|k|φ0
z,ν〉
p0z,ν(x) =
k(r) φ0z,nνlν
(r)
r∫ +∞
04πs2 k(s) |φ0
z,nνlν(s)|2 ds
Yml (θ, ϕ)
P ❲ (HPS
ρ − ǫ0z,ν
)φ0
z,ν = p0z,ν 〈φ0
z,ν|HPS
ρ − ǫ0z,ν |φ0z,ν〉
〈φ0z,ν|p0
z,ν〉 = 1.
❭ Pp0
z,ν
Brz
c(0)
4
❯
φz,1 = φ0z,1, φz,1 = φ0
z,1, pz,1 = p0z,1, ❴ P
ν ≥ 2 P
Pφz,ν
pz,ν
❲
1 ≤ µ, ν ≤ Np, (pz,µ, φz,ν)L2 = δµν .
![Page 66: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/66.jpg)
❴ ❨ ❬❴
pz,µ = Fz,µ
(p0
z,µ −µ−1∑
ν=1
(p0z,µ, φz,ν)L2 pz,ν
)
pz,µ
❴ φz,ν
1 ≤ ν < µ
φz,µ = Fz,µ
(φ0
z,µ −µ−1∑
ν=1
(pz,ν, φ0z,µ)L2 φz,ν
)
φz,µ
❴ pz,ν
1 ≤ ν < µ
φz,µ = Fz,µ
(φ0
z,µ −µ−1∑
ν=1
(pz,ν, φ0z,µ)L2 φz,ν
),
❲
∀ x ∈ IR3 \Brzc(0), φz,ν(x) = φz,ν(x).
❭ Fz,µ
❲
(pz,µ, φz,µ)L2 = 1.
µ = 2
Fz,µ =(1 − (p1, φ
02)L2(p0
2, φ1)L2
)−1/2
P P
φz,ν
φz,ν
pz,ν
❲ ❴ P ❲
l
P
pz,ν
Brz
c(0)
P
φz,ν
φz,ν
❴ Brz
c(0)
P
Brzc(0)
IR3 ❬
Tz
❭
Tz
I+Tz
❲ P ❲ ❴ P ❲ ❲ P
Np
P(φz,ν − φz,ν)
![Page 67: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/67.jpg)
[(pz,µ, φz,ν)]
❲
I+Tz
P 1 ❲
[(pz,µ, φz,ν)]
(φz,ν − φz,ν)
❫ ❯ ❴
φz,ν
❯ P ❯
φz,ν
Brz
c(0)
❫ ❬ ❲❬P ❴ ❲ ❯
k(r) V0
❫ ❬ ❲ ❯ ❯ ❯
❨❯ ❲ P
Tz
P τxkTzk P
∀ u ∈ L2(IR3), ∀ x ∈ IR3, ((τxkTzk
)u) (x) =
Np,k∑
ν=1
(pk,ν(·−xk), u)L2
(φk,ν(x− xk) − φk,ν(x− xk)
).
Nk = Np,k
φkν = φzk,ν(· − xk), φk
ν = φzk,ν(· − xk), pkν = pzk,ν(· − xk), φk
ν = φkν − φk
ν ,
Jku =
Np,k∑
ν=1
(pkν , u)L2 φk
ν , Jku =
Np,k∑
ν=1
(pkν , u)L2 φk
ν ,
Tku = τxkTzk
u = Jku−Jku,
T =M∑
k=1
(Jk − Jk) =M∑
k=1
Tk.
❭ Jk
Jk
Tk
❲
JTk =
Nk∑
ν=1
(φkν , ·) pk
ν, JTk =
Nk∑
ν=1
(φkν , ·) pk
ν,
T Tk =
Nk∑
ν=1
(φkν , ·) pk
ν.
∀ 1 ≤ k ≤M, ∀ x ∈ IR3 \Brzkc
(xk),(Jku)(x) =
(Jku)(x).
![Page 68: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/68.jpg)
P ❲ ❴❨❴❯ Br
zkc
(xk)
φkνP
L2(Brzkc
(xk)) ❲ P
❨❴❯ ∀ 1 ≤ k ≤M, ∀ x ∈ Brzkc,(Jku)
(x) = u(x). ❴❨❴❯ ❭
ψi = (I + T )ψi = ψi +
M∑
k=1
(Jkψi − Jkψi
).
P ❴❨❴❯ ❯
ρ P
❴❬❴ψi
❵ ❴
ρ(x) =
N∑
i=1
|ψi(x)|2 = ρ(x) + ρ1(x) − ρ1(x) + ρinc(x)
❲ρ(x) =
N∑
i=1
|ψi(x)|2
ρ1(x) =
M∑
k=1
ρk1(x), ρ1(x) =
M∑
k=1
ρk1(x), ρinc(x) =
M∑
k=1
ρkinc(x),
❲ρk
1(x) =N∑
i=1
|(Jkψi)(x)|2,
ρk1(x) =
N∑
i=1
|(Jkψi)(x)|2,
ρkinc(x) = 2
N∑
i=1
[ψi(x) (Jkψi)(x) − (Jkψi)(x) (Jkψi)(x) − ψi(x) (Jkψi)(x) + |(Jkψi)(x)|2
].
ρk
inc
Br
zkc
(xk)
❴❨❴❯ ❲ρk
inc = 0
Brzkc
(xk)
❲ ❴❨❴❯ P
![Page 69: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/69.jpg)
A
❬ L2(IR3)
〈ψ,Aψ〉 = 〈ψ, Aψ〉
A ❬
A = A+
M∑
k=1
[T T
k A+ ATk + T Tk ATk
].
A
P P
A = A +M∑
k=1
[JT
k AJk − JTk AJk
]
+
M∑
k=1
[AJk + JT
k A− JTk AJk − JT
k AJk − AJk − JTk A + 2JT
k AJk
].
A
P ❨(Aφ)(x) = v(x)φ(x)
v(x) ❲ P ❴❨❴❯
❬
A = A+
M∑
k=1
[JT
k AJk − JTk AJk
]
= A+
M∑
k=1
[Nk∑
µ,ν=1
|pkµ〉(〈φk
µ, Aφkν〉 − 〈φk
µ, Aφkν〉)
〈pkν |]
;
❲ ❲ ❲ P P
∂Brkc(xk)
❯ ❲ P P
ψi
❫ ❨❯ ❬
❭ ❯ ❴ ❬ ❲ ❴
![Page 70: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/70.jpg)
❨❯ ❬ ❬ ❬ ❴ ❲ ❯ ❲ ❯ ❲ P P
ψi
❭ ❴
❴❬❴ ❬ P ❵❨ ❲ P ❴
vloc ❯
❲ V0
❲ ❴ ❴ ❯
Tz
❨❯ ❫ ❲ ❴ ❨
rzc
Pk(r)
V0 ❬
❴ ❯ ❴ ❴ ❴❬ ❲ ❴ ❴❯ ❫ ❲ ❴
φz,ν
φz,ν
pz,ν
❴ P
![Page 71: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/71.jpg)
![Page 72: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/72.jpg)
❴ [
]
❴ ❲ ❲ ❴ ❯
![Page 73: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/73.jpg)
❳ ❵1 ❯1 2 ❭ ❵1
1
2
❴ ❨ ❴ P ❴
P ❨ P❬❨ ❨ ❴ ❴❨
❨ P
N ❴❨
❲ P ❴ ❬ ❵❬ ❴ P P ❴ P ❴ ❴ ❫ ❴ ❲ ❬ ❴
IR3N ❴ P ❴ P ❲ PN
P ❨ P P ❲ P
N ❴ P
❴ ❨ ❴❬ ❴ ❴ ❫
![Page 74: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/74.jpg)
❳❬ ❴ ❴ ❴ ❴ P ❴ ❴ ❴❲ ❲ ❴ ❴ ❨ ❴ ❴ ❴ P P ❴ P ❨ ❲
H
S ❲
Nb × Nb
❨ Nb × Nb
❨ ❲ ❴Nb > N
D⋆P ❴
Hci = ǫiSci, ǫ1 ≤ . . . ≤ ǫN ≤ ǫN+1 ≤ . . . ≤ ǫNb,
ctiScj = δij ,
D⋆ =N∑
i=1
cicti.
❭ ❴ P ❴ ❴ ❴ ❨ ❴
N ❴ ❴ P
❴ ❴S
❴ ❨ D⋆
❴ ❴ ❴ ❨ ❴
N❲
❴ ❴ N
❲ PH
❴❴ ❴ ❴ ❲ ❴ ❴ ❬ ❴
N❲
ci ❲ ❴
❴ P ❴ S ❴ ❲
χi1≤i≤Nb
P ❴N
❬ ❲P ❴ H
❬ ❬ P ❴ ❴❬❴ ❴❲
Hij =1
2
∫
IR3
∇χi · ∇χj +
∫
IR3
V χiχj
❴V
❬ ❴ S
❴ ❲ ❴ ❴ χi1≤i≤Nb
Sij =
∫
IR3
χiχj.
❴ P ❴ ❭ ❴ ❴ ❲❨ ❴ P❬ χi
❨ ❴ P❬ ❴ ❴ ❴ ❬
❴ ❴ ❭ ❴ ❴ ❲ ❴ P ❴ P
Nb
N
❴ ❴ ❨Nb ∼ 2N
![Page 75: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/75.jpg)
❴ ❴ P ❲ P ❴ ❬ ❲ ❴❴ ❴
H
S
P ❨ P ❴❨ ❴ ❨ ❴ ❴ ❴ ❴ P
H
S P ❴ P
102 ❴ ❴ P ❲ ❴❴ ❴ ❴ ❬
Nb
❴ ❴N
❨Nb ∼ 100N
❴ S
❴ ❨ ❴
H P ❴ ❴ ❴ P
❴ ❴ D⋆
❴ ❲ci
❴ ❲ ❴ ❴ P ❴ P ❴
D⋆
P D⋆ = C⋆C
t⋆
❴
C⋆
❴
inf
(HCCt
), C ∈ MNb,N(IR), CtSC = IN
.
❴ ❴ ❨ P (HCCt
) ❲ ❴ ❨ P(
CtHC) Mk,l ❴ ❲ P ❴
k×l ❴ ❴ ❨ P
C⋆
C⋆UP ❨
❴N ×N
U ❲ ❴ ❴ ❴
N❬
❴ ❲ PH
❨ ❴ ❴(N + 1)
❬❴ ❴ D∗
❨ P ❴ ❴ P ❴ C⋆
P ❴ ❴ ❬❭ P ❴ ❨ P P P ❨ ❴ P P ❴ P
C⋆ = (c1| · · · |cN)
❴ ❴ D⋆
❲ ❴ ❲ ❴
C⋆
❴D⋆
❨ ❴ N
❲ P
H ❴ ❴ ❨ ❴ ❴
NP
❨ ❴103
P ❴ ❴ P ❴ ❨ ❴ ❬
❨N3 P ❴ ❨
N ❴
❴❨ ❴ ❲ ❴ P ❴ ❲❨ ❴ ❴ ❴ ❴ ❴ ❴ ❴ ❴ ❴ ❴ P ❴ P ❴ H
❴ ❲ P ❴ ❴
![Page 76: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/76.jpg)
❴ ❴ P ❴ ❴ ❨ P ❲ ❨ ❴ ❴ ❴ ❨ ❴ ❴ ❴ P ❴ ❴ P ❴ ❲ ❴ ❴ ❲
❴❨ P D⋆
P❨ ❴ P ❴ ❴ ❨ P ❴ ❲ ❴ ❴ ❲ ❨ ❴ ❴ ❨ ❴
N❲ ❴ ❲ P ❴ N
❬❲ ❵ ❴ ❨ ❨ ❴ ❴ P P ❴
H
S ❴
D⋆
❴ H
S
❴ ❴ P ❨ ❴ P ❬
N ❴
❲ P P ❴ ❨ P ❴❬❴ ❴ P ❴❴ ❴ P ❴ ❨
❴ ❴ D⋆
P ❴ ❴ P ❴ ❲
γ =ǫN+1 − ǫNǫNb
− ǫ1.
P ❴ P ❴ ❴ ❲ ❨ ❲ ❴❨ ❴ ❴ ❴ P ❨ ❴ P ❴ P❨
❬ ❴ P ❴❴ P ❴❴ ❴ P ❴ ❴ ❨ ❴
γ ❲❨ ❴
D ❨ ❴
❨ ❴ ❴ P❨ ❲ P ❴ ❴ ❴ ❨ ❴ ❴ ❴P ❲ ❴ ❴ ❴ P P ❴
❵P ❴ ❴ P ❴ ❴ P P P ❴
❴ ❲ ❴ ❴
H P❬❨
D⋆
P P ❴ ❴ ❨
![Page 77: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/77.jpg)
P ❲ ❴ ❲❨ ❴ ❲ ❨ ❴ P ❴ ❴ ❴ ❬
H ❴
❲ S
❴ ❴ ❨ ❲❨❴ ❴ ❴ ❴ P ❴❲ ❴ ❨ ❴ ❨ ❨ ❴ P P ❴ P ❴ ❬ ❨ ❴ ❴ ❲ ❵ P ❲ ❨ ❴❲ P ❴ ❴ P ❲ P ❲ ❴ P P
H
❨ ❲ P ❲ ❴ ❴ ❴
H + δH
S + δS
H
S
❴ P P ❴ ❴ ❴ ❴ ❨ ❴ ❴ ❴❨ ❴ P ❴ P❨
χi
H
bN
bN H =
0
0
H
H
![Page 78: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/78.jpg)
S
N
H
q ≪ Nb
Nb/q
q
n = 2q
C⋆
Ci
H
D
C
D = Dt D2 = D
D) N
H p H =
0
0
n
N = (p+1) n/2 b
H 1
H
0
0
n
N = (p+1) n/2 b
D =
D
q
n = 2q
![Page 79: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/79.jpg)
n = 2q
S = INb
S 6= INb
n = 2q + nbs
2nbs− 1
S
S = INb
χi1≤i≤Nb
S 6= INb
inf
(HCCt
), C ∈ MNb,N(IR), CtC = IN
.
(HCCt
) C
CtC = IN
p∑
i=1
(HiCiC
ti
), Ci ∈ Mn,mi(IR), mi ∈ IN, Ct
iCi = Imi∀ 1 ≤ i ≤ p,
CtiTCi+1 = 0 ∀ 1 ≤ i ≤ p− 1,
p∑
i=1
mi = N
.
T ∈ Mn,n(IR)
Tkl =
1
k − l = q
0
Hi ∈ Mn,n(IR)
H
C 1
0
0
C p
C 1
H p
H 1
C i H i Tr
t
0
0
C p
Trt
C i Σ =
p
i=1
C 1 C 1
0
0
C p
C i C i t
C i+1
0
0
C p
=
t
0
0
ti C T
![Page 80: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/80.jpg)
bN = (p+1) n/2
C 1
C =
0
0 n
C p
m
N = m + ... + m
m
1 p
p
1
C
N(N+1)2
CtC = IN
Nb
∑pi=1
mi(mi+1)2
Ct
iCi = Imi
∑p−1
i=1 mimi+1
Ct
iTCi+1 = 0
n
N
H
HiCi = CiEi + T tCi−1Λi−1,i + TCi+1Λti,i+1 1 ≤ i ≤ p,
CtiCi = Imi
1 ≤ i ≤ p,Ct
iTCi+1 = 0 1 ≤ i ≤ p− 1,
C0 = Cp+1 = 0.
(Ei)1≤i≤p
(Λi,i+1)1≤i≤p−1
CtiCi = Imi
Ct
iTCi+1 = 0
mi × mi
Ei
Λi,i+1
mi×mi+1
L(Ci , Ei , Λi,i+1
)=
p∑
i=1
Tr(HiCiC
ti
)+
p∑
i=1
Tr( (Ct
iCi − Imi
)Ei
)
+
p−1∑
i=1
Tr(Ct
iTCi+1Λti,i+1
).
![Page 81: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/81.jpg)
(mi)1≤i≤p
inf〈H1Z1, Z1〉 + 〈H2Z2, Z2〉, Zi ∈ IRNb , 〈Zi, Zi〉 = 1, 〈Z1, Z2〉 = 0
.
〈·, ·〉
IRNb
(Z0
1 , Z02)
(Zk
1 , Zk2 )k∈IN
(Zk1 , Z
k2 )
Zk
1 =〈H1Z1, Z1〉, Z1 ∈ IRNb , 〈Z1, Z1〉 = 1 〈Z1, Z
k2 〉 = 0
,
Zk2 =
〈H2Z2, Z2〉, Z2 ∈ IRNb , 〈Z2, Z2〉 = 1 〈Zk1 , Z2〉 = 0
;
α∗ =〈H1Z1, Z1〉 + 〈H2Z2, Z2〉, α ∈ IR
Z1 =Zk
1 + αZk2√
1 + α2, Z2 =
−αZk1 + Zk
2√1 + α2
,
Zk+1
1 =Zk
1 + α∗Zk2√
1 + (α∗)2, Zk+1
2 =−α∗Zk
1 + Zk2√
1 + (α∗)2.
k
Z2 = Zk
2
Z1
Zk1
Z1 = Zk
1
Z2
Zk
2
〈Z1, Z2〉 = 0
![Page 82: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/82.jpg)
〈H1Z1, Z1〉 + 〈H2Z2, Z2〉
Zk
1
Zk
2
Zk1
Zk
2
Nb
2Nb
2
H1 = H2 = H
H
α
H
H
ǫ2 − ǫ1 < ǫ3 − ǫ2 ǫ3 − ǫ2
p
(Ci)1≤i≤p
E((Ci)1≤i≤p
)=
p∑
i=1
(HiCiC
ti
),
U0 = Up = 0.
ǫ
k
(mki )1≤i≤p
(Cki )1≤i≤p
Ck
i ∈ Mn,mki (IR)
[Ck
i ]tCki = Imk
i
[Ck
i ]tTCki+1 = 0
(mk+1i )1≤i≤p
(Ck+1
i )1≤i≤p
•
![Page 83: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/83.jpg)
i
H2i+ǫ
V k2i+ǫ =
x ∈ IRn,
[Ck
2i+ǫ−1
]tTx = 0, xtTCk
2i+ǫ+1 = 0,
P k
2i+ǫH2i+ǫPk2i+ǫ
P k
2i+ǫ
V k2i+ǫ
n−mk
2i+ǫ−1 −mk2i+ǫ+1
λk
2i+ǫ,1 ≤ λk2i+ǫ,2 ≤ · · ·
xk2i+ǫ,j
T
Cki−1
Ck
i+1
(λk2i+ǫ,j)i,j
∑
i
m2i+ǫ
i
#2i + ǫ
λk2i+ǫ,j
mk
2i+ǫ
i
mk2i+ǫ
xk
2i+ǫ,j
n× mk
2i+ǫ
C
k
2i+ǫ
i
H2i+ǫ+1
V k2i+ǫ+1 =
x ∈ IRn,
[C
k
2i+ǫ
]tTx = 0, xtTC
k
2i+ǫ+2 = 0
λk
2i+ǫ+1,1 ≤ λk2i+ǫ+1,2 ≤ · · ·
xk
2i+ǫ+1,j
T
Ck
2i+ǫ
C
k
2i+ǫ+2
(λk2i+ǫ+1,j)i,j, (λ
k2i+ǫ,j)i,j
N
l
#l
λkl,j
(mk+1l )1≤l≤p
Ckl =
[xk
l,1| · · · |xkl,mk+1
l
]
ǫ
1 − ǫ
• U∗ =
f(U), U = (Ui)i, ∀1 ≤ i ≤ p− 1 Ui ∈ Mmi+1,mi(IR)
,
f(U) = E((
Ci(U)(Ci(U)tCi(U)
)− 1
2
)
i
),
Ci(U) = Ck
i + TCki+1Ui
([Ck
i ]tTT tCki
)− T tCk
i−1Uti−1
([Ck
i ]tT tTCki
).
1 ≤ i ≤ p
Ck+1i = Ci
(U∗) (
Ci
(U∗)tCi
(U∗))−1/2
.
[Ck+1
i
]tTCk+1
i+1 = 0
T 2 = 0
![Page 84: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/84.jpg)
C2i
C2i+1
U = (Ui)i
α
(Ci)1≤i≤p
r
(Gl)1≤l≤r
(G2l+1)
(G2l)
N
81 32 4 5 6 7 9 10
G1
G2
G3
p = 10
r = 3
ǫ = 1
(mk2i+1, C
k
2i+1)i
∑i
(H2i+1C2i+1C
t2i+1
), C2i+1 ∈ Mn,m2i+1(IR), Ct
2i+1C2i+1 = Im2i+1,
[Ck2i]
tTC2i+1 = 0, Ct2i+1TC
k2i+2 = 0,
m2i+1 ∈ IN,∑
i
m2i+1 =∑
i
mk2i+1
.
![Page 85: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/85.jpg)
Ck
2i
Ck
2i+1 p
p/2
n
p∑
i=1
(HiCiC
ti
), Ci ∈ Mn,mi(IR), Ct
iCi = Imi, mi ∈ IN,
∑
i
mi = N
[Ck
2j−1]tTC2j = 0, [C2j ]
tT [C2j+1]k = 0,
0 ≤ m2j+1 ≤ mk2j+1, C2j+1 ⊂ C
k
2j+1
,
C2j+1 ⊂ C
k
2j+1
C2j+1
C
k
2j+1
ǫ = 1
ǫ
1 − ǫ
mi
T
∀ 1 ≤ i ≤ p− 1,∥∥[Ck
i ]tTCki+1
∥∥ ≤ ǫL,
ǫL > 0
m1 = m2 = m
m1
m2
(Ck+11 , Ck+1
2 ) = (Ck1 , C
k2 )
CtTCk
2 = 0
ǫ = 1
[Ck
1 ]tTC = 0
ǫ = 1
![Page 86: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/86.jpg)
U∗ = 0
k
k
(C1, C2) = (Ck1 , C
k2 )
f(U) =(
J1(U)C1(U)tH1C1(U))
+(
J2(U)C2(U)tH2C2(U))
Ji(U) =
(Ci(U)tCi(U)
)−1 i = 1, 2
(J1(U)
)−1
= Im +(Ct
1TTtC1
)U t(Ct
2TtTC2
)U(Ct
1TTtC1
),
(J2(U)
)−1
= Im +(Ct
2TtTC2
)U(Ct
1TTtC1
)U t(Ct
2TtTC2
),
∇J1(0) = ∇J2(0) = 0
U
m
1 ≤ i, j ≤ m
1
2
∂f
∂Uij
(0) =([ ∂C1
∂Uij
(0)
]t
H1C1
)+
([ ∂C2
∂Uij
(0)
]t
H2C2
)
=
((Ct
1TTtC1
)Ct
1H1TC2
)
ji
−(Ct
1TH2C2
(Ct
2TtTC2
))
ji
=
((Ct
1TTtC1
)(Λ1 − Λ2)
(Ct
2TtTC2
))
ji
,
Λ1
Λ2
H1C1 = C1E1 + TC2Λt1,
H2C2 = C2E2 + T tC1Λ2.
U∗ = 0
∀ 1 ≤ i, j ≤ m∂f
∂Uij(0) = 0,
Λ1 = Λ2 (
Ct1TT
tC1
) (Ct
2TtTC2
)
n ≫ 2m
(C1, C2)
n
2m
∀ 1 ≤ i ≤ p, Ci(U) = Cki + TCk
i+1Ui
([Ck
i ]tTT tCki
)− T tCk
i−1Uti−1
([Ck
i ]tT tTCki
)
![Page 87: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/87.jpg)
Ck
i
Ck
i
Cki
Ck
i
Cki
n
m
(Ck+1
i )i
nm
P1
P2
P3
3
2
nm
3
P1
P2
P3
nm
P1
P2
H
S
nm
H
S
D⋆
![Page 88: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/88.jpg)
P1 P2 P3
S
H
H 10−12 10−12 10−10
D 10−11 10−11 10−7
m1 = 67 m1 = 105 m1 = 136 mp = 67 mp = 106 mp = 137 mi = 56 mi = 84 mi = 104
nm
nm
H
S
H
S
D⋆
E0
nm
γ
H
nm
P1
P2
P3
mi
Ci
I1
C
I2
Ci
mi
Hi
Si
H
S
H
S
I1 I3
I2
![Page 89: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/89.jpg)
eE =|E − E0|
|E0|
L∞
e∞ = sup(i,j) |Hij |≤ε
∣∣∣Dij − [D⋆]ij
∣∣∣ ,
ε = 10−10
H
ε
D
Tr(AD)
A
H
(i, j)
|Hij |
mi
S1
S2
S3
S4
P1
nm = 801
![Page 90: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/90.jpg)
Nb = 8050
N = 5622 P2 P3
E0 = −27663.484
p = 100
99
2
S1
S2
S3
S3 I2
I3
S3
S4
I3
S2
I2
0 500 1000 1500 2000 2500 3000 3500 4000 450010
−11
10−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−1
CPU Time in seconds
En
erg
y e
rro
r
S1S2S3S4
0 500 1000 1500 2000 2500 3000 3500 4000 450010
−4
10−3
10−2
10−1
100
CPU Time in seconds
De
nsity e
rro
r
S1S2S3S4
I1
S4
P2 P3
I2
I3
![Page 91: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/91.jpg)
0 500 1000 1500 2000 2500 3000 3500 4000 4500 500010
−12
10−10
10−8
10−6
10−4
10−2
100
CPU Time in seconds
En
erg
y e
rro
r
S1S2S3S4
0 500 1000 1500 2000 2500 3000 3500 4000 4500 500010
−4
10−3
10−2
10−1
100
CPU Time in seconds
De
nsity e
rro
r
S1S2S3S4
I2
S
QR
N3b
N2b
Nb
D = CCt
C
![Page 92: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/92.jpg)
N
Nb
P1
P2
P3
nm
Nb
103105
I3
Nb
104
Nb
Nb
Nb
P3
103
P1
4001
Nb = 40050
P2
2404
Nb = 24080
P3
208
Nb = 854
I3
![Page 93: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/93.jpg)
102
103
104
105
106
100
101
102
103
104
105
106
107
108
109
Nb
CP
U T
ime
in
se
co
nd
s
MDD LAPACKDMM
102
103
104
105
106
103
104
105
106
107
108
109
Nb
Me
mo
ry r
eq
uire
me
nt
in K
byte
s
MDD LAPACKDMM
P1
102
103
104
105
106
101
102
103
104
105
106
107
108
Nb
CP
U T
ime
in
se
co
nd
s
MDD LAPACKDMM
102
103
104
105
106
103
104
105
106
107
108
109
Nb
Me
mo
ry r
eq
uire
me
nt
in K
byte
sMDD LAPACKDMM
P2
102
103
104
105
101
102
103
104
105
106
107
Nb
CP
U T
ime
in
se
co
nd
s
MDD LAPACK
102
103
104
105
103
104
105
106
107
108
Nb
Me
mo
ry r
eq
uire
me
nt
in K
byte
s
MDD LAPACK
P3
![Page 94: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/94.jpg)
‖Dn −Dn−1‖ ≥ ‖Dn−1 −Dn−2‖ ‖Dn −Dn−1‖ ≤ ǫa
ǫa
ǫa = 10−4 ǫa = 10−3
P1
P2
P3
P2
I3
nm
nm
![Page 95: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/95.jpg)
0 2 4 6 8 10 12 14 16 18
x 104
10−5
10−4
10−3
10−2
10−1
100
CPU Time in seconds
De
nsi
ty e
rro
r
MDD DMM MDD+DMM
P1
0 2 4 6 8 10 12 14
x 104
10−5
10−4
10−3
10−2
10−1
100
CPU Time in seconds
De
nsi
ty e
rro
r
MDD DMM MDD+DMM
P2
0 0.5 1 1.5 2 2.5 3 3.5
x 104
10−4
10−3
10−2
10−1
100
CPU Time in seconds
De
nsi
ty e
rro
r
MDD DMM MDD+DMM
P3
P1 4001
P2
2404
P3 208
![Page 96: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/96.jpg)
C
P3
ǫa
![Page 97: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/97.jpg)
0 2 4 6 8 10 12 14
x 104
10−14
10−12
10−10
10−8
10−6
10−4
10−2
100
CPU Time in seconds
En
erg
y a
nd
de
nsity e
rro
rs
Energy error Density error
P1
0 5 10 15
x 104
10−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−1
100
CPU Time in seconds
En
erg
y a
nd
de
nsity e
rro
rs
Energy error Density error
P2
1 2 3 4 5 6 7 8 9 10 11
x 104
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−1
100
CPU Time in seconds
En
erg
y a
nd
de
nsi
ty e
rro
rs
Energy error Density error
P3
P1
20001
Nb = 200050
P2
12004
Nb = 120080 P3
5214
Nb = 36526
![Page 98: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/98.jpg)
N = 2
X
Y
C
inf(HX,X) + (HY, Y ), (X,X) = 1, (Y, Y ) = 1, (X, Y ) = 0
.
(Xn, Y n)
(Xn+1, Y n+1)
(I)
∣∣∣∣X =
(HX,X), (X,X) = 1, (X, Y n) = 0
,
Y =
(HY, Y ), (Y, Y ) = 1, (X, Y ) = 0.
(II)
∣∣∣∣∣∣∣∣∣∣∣∣∣
(Xn+1, Y n+1) =(Xα∗ , Yα∗
),
α∗ = infα∈IR
(HXα, Xα) + (HYα, Yα)
,
(Xα, Yα) =
(X + αY√
1 + α2,Y − αX√
1 + α2
).
(II)
∀ α ∈ IR, (HXα, Xα) + (HYα, Yα) = (HX,X) + (HY, Y ).
(ci, ǫi)
H
ǫi
Dci
=cic
ti
(Dn
X , DnY , D
n)n
∀ n ∈ IN, (DnX , D
nY , D
n) =(Xn(Xn)t, Y n(Y n)t, Dn
X +DnY
).
HXn+1 = λn+1Xn+1 + νn+1Y n,HY n+1 = δn+1Xn+1 + µn+1Y n+1,
(Xn+1, Y n) = (Xn+1, Y n+1) = 0,(Xn+1, Xn+1) = (Y n+1, Y n+1) = 1.
![Page 99: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/99.jpg)
λn+1 = (HXn+1, Xn+1),νn+1 = (HXn+1, Y n),δn+1 = (HY n+1, Xn+1),µn+1 = (HY n+1, Y n+1).
(λn)n, (νn)n, (δn)n, (µn)n
(λn)n
(µn)n
λ
µ
(λn)n
(µn)n,
ǫ1
(HXn+1, Xn+1
)≤(HXn, Xn
)
λn+1 ≤ λn.
(λn)n
λ
(µn)n
µ
λ ≤ µ.
X0
Y 0
(δn)n
(νn)n
ν
νn ≥ 0
δn ≥ 0
n
Xn Y n
(δn)n
(νn)n
H
δn+1 = (HY n+1, Xn+1),= (Y n+1, HXn+1),= λn+1(Y n+1, Xn+1) + νn+1(Y n+1, Y n),= νn+1(Y n+1, Y n).
νn+1 = δn(Xn+1, Xn).
Xn
Yn
|(Xn+1, Xn)| ≤ 1 |(Y n+1, Y n)| ≤ 1.
![Page 100: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/100.jpg)
. . . ≤ δn+1 ≤ νn+1 ≤ δn ≤ νn ≤ . . .
(δn)n
(νn)n
ν
ν = 0
(Dn
X)n
(DnY )n
(DnX , λ
n)n(Dn
Y , µn)n
(Dcα, ǫα)
(Dcβ
, ǫβ)α < β
H
∀ 1 ≤ k ≤ Nb, (ck, HX
n+1) = (Hck, Xn+1),
= ǫk(ck, Xn+1).
∀ 1 ≤ k ≤ Nb, (ǫk − λn+1)(ck, X
n+1) = νn+1(ck, Yn).
|(ck, Y n)| ≤ 1
ν∀ 1 ≤ k ≤ Nb, lim
n→∞(ǫk − λ)(ck, X
n+1) = 0.
λ
H
(ci)i
IRNb
(Xn)n
∀ k 6= α, limn→∞
(ck, Xn+1) = 0,
(DnX)n
Dcα
H
µ
ǫβH
(Dn
Y )n
Dcβ
α 6= β
α < β
(Dα, Dβ)
(Dc1, Dc2) i0
1, . . . , Nb
ei0
ei0 = ci0 − (ci0 , X
n+1)Xn+1.
(Xn+1, ei0) = 0
µn+1 ≤(H
ei0
‖ei0‖,ei0
‖ei0‖
),
≤ 1
‖ei0‖2
(ǫi0 − 2ǫi0(ci0, X
n+1)2 + λn+1(ci0, Xn+1)2
).
![Page 101: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/101.jpg)
∀ i0 6= α, β ǫα < ǫβ ≤ ǫi0 .
α
β
1
i0 = 1
α
1
β 6= 2
i0 = 2
ν 6= 0
((Xn, Y n
))n
(
Xp(n), Y p(n))
n
limn→∞
Xp(n) = Xp, limn→∞
Y p(n)) = Yp.
DXp
DYp
Dp
(DXp, DYp
, Dp) =(Xp(Xp)
t, Yp(Yp)t, DXp
+DYp
).
(Xp, Yp)
(cαp, cβp
)H
ǫαp
< ǫβp
p(n)
δp(n)+1 = νp(n)+1(Y p(n)+1, Y p(n)
).
ν 6= 0
limn→∞
(Y p(n)+1, Y p(n)
)= 1.
∣∣Y p(n)+1−Y p(n)∣∣ = 2
(1−(Y p(n)+1, Y p(n)
))
limn→∞
∣∣Y p(n)+1, Y p(n)∣∣ = 0.
Xp(n)
limn→∞
Xp(n)+1 = Xp, limn→∞
Y p(n)+1 = Yp.
![Page 102: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/102.jpg)
p(n)(Xp, Yp)
(S)
HXp = λXp + νYp,HYp = νXp + µYp.
(Xp, Yp)
H
(cαp, αp)
(cβp
, βp)
H
H
ǫαp=
(λ+ µ−
√(µ− λ)2 + ν2
)/2,
ǫβp=
(λ+ µ+
√(µ− λ)2 + ν2
)/2,
ǫαp+ ǫβp
= λ+ µ.
ν 6= 0
(DXp
, DYp)
(Dcαp
, Dcβp)
(cαp, cβp
)
Dp
p
(α, β)
(αp, βp) α < β (DXp
, DYp)
p
λµ
(ǫαp
, ǫβp)
p
Dp
(Xp, Yp)
p
p
DXp+Dyp
= Dcα+Dcβ
.
Xp
Yp
θp ∈]0, 2π[
cos(θp)sin(θp) 6= 0
Xp = cos(θp)cα + sin(θp)cβ,Yp = −sin(θp)cα + cos(θp)cβ .
H
(S)
sin(θp)cos(θp) =ν
ǫβ − ǫα,
sin2(θp) =λ− ǫαǫβ − ǫα
=ǫβ − µ
ǫβ − ǫα,
cos2(θp) =ǫβ − λ
ǫβ − ǫα=
µ− ǫαǫβ − ǫα
.
![Page 103: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/103.jpg)
DXp
DYp
∣∣∣∣∣∣
DXp= cos2(θp)Dcα
+ cos(θp)sin(θp)(cβctα + cαc
tβ) + sin2(θp)Dcβ
DYp= sin2(θp)Dcα
− cos(θp)sin(θp)(cβctα + cαc
tβ) + cos2(θp)Dcβ
.
(DXp
, DYp)
p
α
1 β
2
ǫ3 − ǫ2ǫ3 − ǫ1
>1
2.
i0
1, . . . , Nb
ei0
ei0 = ci0 −(ci0 , X
p(n)+1)Xp(n)+1.
(Xp(n)+1, ei0
)= 0
µp(n)+1 ≤(H
ei0
‖ei0‖,ei0
‖ei0‖
),
≤ 1
‖ei0‖2
(ǫi0 − 2ǫi0(ci0 , X
p(n)+1)2 + λp(n)+1(ci0 , Xp(n)+1)2
).
∀ i0 6= α, β limn→∞
(Xp(n), ci0
)= lim
n→∞
(Xp(n)+1, ci0
)=(Xp, ci0
)= 0.
∀ i0 6= α, β λ ≤ µ ≤ ǫi0 .
ǫα < λ ≤ µ < ǫβ .
α
1
ǫα < λ < ǫ1,
α
1
β
2
λ+ µ = ǫ1 + ǫβ ,
≤ 2ǫ2.
![Page 104: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/104.jpg)
2ǫ2 < ǫ1 + ǫβ .
H
β
2
2ǫ2 < ǫ1 + ǫ3 < ǫ1 + ǫk
k ≥ 3
ǫ3 − ǫ2ǫ3 − ǫ1
>1
2.
H
3 ǫ1, ǫ2, ǫ3
ǫ3−ǫ2 < ǫ2−ǫ1
(X0, Y0)
(X0, Y0) =((e1 + e3)/
√2, (e1 − e3)/
√2),
(λ0, µ0
)=(X t
0HX0, Yt0HY0
)=((ǫ1 + ǫ3)/2, (ǫ1 + ǫ3)/2
).
ǫiλ0
µ0 ǫ2
n ∈ IN(Xn, Yn) = (±X0,±Y0)
(Xn, Yn)
(e1, e3)
H
(Dn)n
(cα, cβ)
H
α = 1α < β
β
2
ǫ3 − ǫ2ǫ3 − ǫ1
>1
2.
![Page 105: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/105.jpg)
![Page 106: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/106.jpg)
[
] [
][
] [
]
![Page 107: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/107.jpg)
S
INb
D⋆
Fci = ǫici, ǫ1 ≤ . . . ≤ ǫN ≤ ǫN+1 ≤ . . . ≤ ǫNb,
cticj = δij ,
D⋆ =N∑
i=1
cicti.
F
−1 = ǫ1 ≤ . . . ≤ ǫN < 0 < ǫN+1 ≤ . . . ≤ ǫNb≤ 1.
ǫ = ǫN+1 = −ǫN ≃ 0
ǫF
0
γF
ǫN+1 − ǫNǫNb
− ǫ1=
2ǫ
ǫNb− ǫ1
≃ 0
D⋆D⋆ = H(F )
H
[−1, 1] H = 1I[−1,0]
D = q(F )
q ∣∣∣∣ q(ǫi) = 1 ∀ 1 ≤ i ≤ N
q(ǫi) = 0 ∀ N + 1 ≤ i ≤ Nb.
H
ǫi
H [−1,−ǫ]∪[ǫ, 1]
H
ǫ
q
ǫ
γ
D⋆
γ
γ ≃ 0
D⋆
q(F )
O(N3)
Nb = O(N)
![Page 108: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/108.jpg)
F
0
2d
(ci, ǫi)i=N−d+1,N+d
F
(ci, ǫi) 2d
F−1
ǫF = 0
F−1
d
F
F = F −N+d∑
i=N−d+1
ǫicicti,
γF =
ǫN+d+1 − ǫN−d
2ǫγF
F
Dq
D⋆
Dq = q(F ) +N+d∑
i=N−d+1
1Iǫi<0 cicti.
F
F
nm
N
A F−1 Id
(., .)
IRNb
p
(uj, λj)i=j,p
A
O(Nb)
F
![Page 109: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/109.jpg)
(u1, λ1)
A
• v1
‖v1‖ = 1
• i = 1,
(vi, θi)⋄ v = Avi,⋄ θi = vt
iv,⋄ vi+1 = v/‖v‖,
•
(v1, u1) 6= 0(vi, θi)
(u1, λ1)
p > 1
• V1 ∈ M(Nb, p)
• i = 1,
Vi
⋄ V = AVi,⋄
V V = QR
⋄ Vi+1 = Q,•
∀ 1 ≤ j ≤ p, V t1uj 6= 0IRp
Vi
uj
V
QR
p = Nb
QR
ui
A
• V1 ∈ M(Nb, p)
P1
• i = 1,
⋄ V = Pi(A)Vi
⋄
V V = Vi+1R
⋄ Hi+1 = V t
i+1AVi+1
⋄ p
(zj
i , µji )j=1,p
Hi+1
⋄
Pi+1
µj
i
•
![Page 110: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/110.jpg)
V1
(zj
i , µji )j
(uj, λj)j
Pi
Hi
Vi
• v1
‖v1‖ = 1
• m = 1
Km(A, v1) n
⋄ Vm
Km(A, v1)
⋄ Hm = V t
mAVm
⋄ (zj
m, µjm)j=1,p
Hm
• Km(A, v1)
A
m
v1
Km(A, v1) =
(v1, Av1, ..., Am−1v1).
V
H
Vm
Hm
• v1
‖v1‖ = 1
• 1 ≤ j ≤ m− 1
⋄ w = Avj
⋄ i = 1, j
(hm)ij = (w, vi)
w = w − (hm)ijvi
⋄ ⋄ (hm)j+1,j = ‖w‖⋄ vj+1 =
w
(hm)j+1,j
•
∀ m (hm)m+1,m 6= 0
∣∣∣∣∣∣
AVm = Vm+1Hm+1,m,= VmHm + (hm)m+1,mvm+1e
tm,
V tmAVm = Hm.
![Page 111: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/111.jpg)
(ei)
Hm
(hm)ij
Hm+1,m
Hm
(hm)m+1,m
Vm
V t
mVm = Im
0
Hm
(hm)m+1,m
Hm+1,m =
A
Hm
(hm)m+1,m = 0 ⇐⇒ AKm(A, v1) = Km(A, v1).
Km(A, v1)
(zj
m, µjm)
(uj, λj)
(v1, uj) 6= 0
p
Vm
m
m
vnew1
(zj , µj)
Hp+1,p
p
p
vnew1
• v1
m
p+m
H(1)p+1,p
p
v1
• k = 1
p
(uj , λj)
vt1uj 6= 0
![Page 112: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/112.jpg)
⋄ H
(k)p+m
m
H
(k)p+1,p
⋄ (zj
k, µjk)j=1, p+m
A
H
(k)p+m
⋄
vk+1
(zj
k, µjk)
⋄
H(k+1)p+1,p
p
vk+1
•
H
(+1)+1 (νi
k)i=1,m
µjk
vk+1
vk+1 = γ
m∏
i=1
(A− νikId)vk
γ
Hk+1
p+1,p
m
QR
p
vk+1
QR
H
Q
R H = QR
H
Q
QR
Q
QR
QR
k
k
jhj+1,j
(hj)j+1,j
AVp+m = Vp+mHp+m + hp+m+1,p+mvp+m+1etp+m.
ν1 QR
Hp+m − ν1Id = Q1R1
(A− ν1Id)Vp+m − Vp+m(Hp+m − ν1Id) = hp+m+1,p+mvp+m+1etp+m
(A− ν1Id)Vp+m − Vp+mQ1R1 = hp+m+1,p+mvp+m+1etp+m.
Q1
AV(1)p+m − V
(1)p+mH
(1)p+m = hp+m+1,p+mvp+m+1e
tp+mQ1
![Page 113: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/113.jpg)
H
(1)p+m = R1Q1 + ν1Id = Qt
1Hp+mQ1
V
(1)p+m = Vp+mQ1
Q1
V(1)p+m
Q1
H(1)p+m
hp+m+1,p+metp+mQ1 = βet
p+m−1 +
p+m+1∑
i=p+m
αieti.
m
νi
µj
AV(m)p+m − V
(m)p+mH
(m)p+m = hp+m+1,p+mvp+m+1e
tp+mQ
H(m)p+m = QtHp+mQ
V
(m)p+m = Vp+mQ
Q = Q1Q2 . . . Qm
Qi
QR
hp+m+1,p+metp+mQ = βet
p +
p+m+1∑
i=p+1
αieti.
Vp
p
V(m)p+m
Qi Vtp Vp = Ip
AVp = VpH(m)p + h
(m)p+1,pv
(m)p+1e
tp + βvp+m+1e
tp,
= VpHp + hp+1,pvp+1etp,
Hp = H
(m)p
vp+1 =
h(m)p+1,pv
(m)p+1 + βvp+m+1
‖h(m)p+1,pv
(m)p+1 + βvp+m+1‖
=h
(m)p+1,pv
(m)p+1 + βvp+m+1
hp+1,p
vp+1
Vp
Q
Vp+m+1
p
v1
Vp
F
q(F )
Dq
q(F )
Dq = q(F ) +
N+d∑
i=N−d+1
1Iǫi<0 cicti
F
Dq = q(F −
N+d∑
i=N−d+1
ǫicicti
)+
N+d∑
i=N−d+1
1Iǫi<0 cicti.
![Page 114: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/114.jpg)
F
F
k ∈ IN
(F −
N+d∑
i=N−d+1
ǫicicti
)k
= F k −N+d∑
i=N−d+1
ǫki cicti.
Dq = q(F ) −N+d∑
i=N−d+1
q(ǫi)cicti +
N+d∑
i=N−d+1
1Iǫi<0 cicti,
Dq = q(F ) +
N+d∑
i=N−d+1
(1Iǫi<0 − q(ǫi)
)cic
ti.
D⋆
q(F )
q
1
ci
Dq
(Dq)ij
Fij 6= 0
C(d,N)
2d
N
D⋆ Nb = O(N)
d⋆(N) C(d,N)
l
l O(N)
F
Nb = O(N)
Hl,l
F
F−1 Nb = O(N)
l − 2d
QR
O(l2)
Hl,l
l = O(p)
p
p
2d
F
O(dN + d3)
![Page 115: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/115.jpg)
q(F )
O(NM2)
F
O(n2
mNM)
F
d⋆
MdN
d
γd
F
ed = sup
|x|≥γd/2
|H − q|
ǫγd
M
q
sup|x|≥γd/2
|H − qM | < ed
γd
γd = O( dN
).
Eex = (FD⋆)Eap = (FDq)
∣∣∣Eex − Eap
Eex
∣∣∣ ≤ ǫ 2ed(Nb − 2d)
γd(N − d)≤ ǫ.
ed
M = O(− log2(ed)γ
−1d
).
C(d,N)
C(d,N) = O(N3
d2+Nd+ d3
).
d⋆ = O(N3/5) C = O(N9/5)
d⋆ = O(N1/3) C = O(N4/3)
![Page 116: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/116.jpg)
q
Ω(F,D) = (F (3D2 − 2D3))
F
F
D0 = p0(F )
D0
F
F = F −N+d∑
i=N−d+1
ǫicicti,
D0 = D0 −N+d∑
i=N−d+1
p0(ǫi)cicti,
(D0, F )
Dn
Dn
2d
(pi
n)i=N−d+1,N+d
Dn = Dn −N+d∑
i=N−d+1
pincic
ti.
Dn
pn
Gn
Ω(F , D)
Dn−1
Gn =(∇DΩ(F , D)
)(Dn−1
),
= 6F Dn−1
(1 − Dn−1
)
Ω(F,D)
Gn = 6FDn−1
(1 −Dn−1
)+
N+d∑
i=N−d+1
gincic
ti,
=(∇DΩ(F,D)
)(Dn−1
)+
N+d∑
i=N−d+1
gincic
ti
![Page 117: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/117.jpg)
N − d+ 1 ≤ i ≤ N + d
gi
n = −6ǫipin−1(1 − pi
n−1)
Xn
Xn = −Gn + βCGn Xn−1,
= Xn +N+d∑
i=N−d+1
xincic
ti
N − d+ 1 ≤ i ≤ N + d
xi
n = −gin + βCG
n xin−1
αn
αn =
loc
Ω(F , Dn−1 + αXn
), α ∈ IR
,
=
loc
Ω(F,Dn−1 + αXn
)+
N+d∑
i=N−d+1
ǫi
(3f 2
n,i(α) − 2f 3n,i(α)
), α ∈ IR
N − d+ 1 ≤ i ≤ N + d
fn,i(α) = −pi
n−1 + αxin
Dn
Dn = Dn−1 + αnXn,
= Dn −N+d∑
i=N−d+1
pincic
ti.
Dn = Dn−1+αnXn
N−d+1 ≤ i ≤ N+d
pi
n = pin−1−αnx
in
Dq
D⋆
Dq = Dn +
N+d∑
i=N−d+1
1Iǫi<0 cicti,
Dq = Dn +
N+d∑
i=N−d+1
(1Iǫi<0 − pi
n
)cic
ti.
ci
pn
S = INb
![Page 118: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/118.jpg)
F
F
F =
α β 0 · · · 0
β
0
0 β
0 · · · 0 β α
, Nb = 3N, (α, β) =
(2,−1
).
γNb= 1/Nb
Nb
ǫa
ǫD < ǫa
ǫD = maxi,j |∆Dij|
(∆D)ij = 2D⋆
ij − (Dq)ij
|D⋆ij| + |(Dq)ij|
Fij 6= 0,
= 0
N
d⋆(N)
d
d⋆
C1
nop
C2
tCPU
Nb103
104 ǫa = 10−3
Nb
n∗op
t∗CPU
C1
C2
O(N2)
O(N2.4)
Nb
O(N9/5)
![Page 119: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/119.jpg)
102
103
104
100
101
102
103
104
105
N
t CP
U*
,no
p*
Résultats obtenus avec FOE
C1
C2
102
103
104
102
103
104
105
N
t CP
U*
,no
p*
Résultats obtenus avec DMM
C1
C2
Nb
104 4.104
ǫa = 10−2
200
n∗
op
t∗CPU
C1
C2
O(N1.3)
O(N1.6)
O(N4/3)
10%
20%
N
d⋆
![Page 120: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/120.jpg)
S 6= INB
103
104
105
100
101
102
103
104
N
t CP
U*
,no
p*
Résultats obtenus avec FOE
C1
C2
103
104
105
102
103
104
105
N
t CP
U*
,no
p*
Résultats obtenus avec DMM
C1
C2
S 6= INb
S 6= INb
Fci = ǫiSci, ǫ1 ≤ . . . ≤ ǫN ≤ ǫN+1 ≤ . . . ≤ ǫNb,
ctiScj = δij ,
D⋆ =N∑
i=1
cicti.
Dq
F =
Nb∑
i=1
ǫiScictiS,
q(F ) 6=
Nb∑
i=1
p(ǫi)cicti
q(S−1F ) = q
( Nb∑
i=1
ǫicictiS
)=
Nb∑
i=1
q(ǫi)cictiS,
D⋆ = q
(S−1F
)S−1
q
∣∣∣∣q(ǫi) = 1
1 ≤ i ≤ N
q(ǫi) = 0
N + 1 ≤ i ≤ Nb.
![Page 121: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/121.jpg)
(ci, ǫi)i=N−d+1,N+d
2d
F
F
F = F −N+d∑
i=N−d+1
ǫiScictiS.
Dq
D⋆
Dq = q(S−1F )S−1 +N+d∑
i=N−d+1
1Iǫi<0 cicti.
S = INb
Dq
q(S−1F )S−1
1
S
ci
q(S−1F ) = q(S−1F ) −N+d∑
i=N−d+1
q(ǫi)cictiS.
Dq
Dq = q(S−1F )S−1 −N+d∑
i=N−d+1
q(ǫi)cicti +
N+d∑
i=N−d+1
1Iǫi<0 cicti,
Dq = q(S−1F )S−1 +
N+d∑
i=N−d+1
(1Iǫi<0 − q(ǫi)
)cic
ti,
= Dq +N+d∑
i=N−d+1
(1Iǫi<0 − q(ǫi)
)cic
ti.
Dq
DqS = q(S−1F )
S−1F
S(S−1F ) = F
S−1F
q
Dq D⋆
Dq
E = (FDq) = (S−1FDqS)
Ω(F,D) = (3DSD − 2DSDSD)
D0
D0SF = FSD0
F−1S
O(Nb)
F
S
![Page 122: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/122.jpg)
S 6= INB
d⋆
F
SCF
S 6= INb
S−1
![Page 123: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/123.jpg)
![Page 124: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/124.jpg)
[
]
![Page 125: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/125.jpg)
(Pµ)
nh
L
(µj)j=1,L
nh
(Pµ)
µj
uµj
(Pµj
)
m
uµk
m
(µk)k=1,m
m≪ nh
uµj
(Pµ)
(Pµ)
uµ
uµ
m
uµk
(Pµ)
m
(αµ,k)1≤k≤m
uµ
uµ =
m∑
k=1
αµ,kuµk
uµ
m
µk
L2 H1 uµj
uµ
uµj
uµ1
(Pµ1
) B = uµ1
j = 2, L
uµj
ε
uµj
uµj
ε
uµj
B = B ∪ uµj
(Pµ)
![Page 126: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/126.jpg)
µ
∣∣∣∣−u+ µu = f
Ω,
u = 0
∂Ω.
Ω
IR3
f ∈ L2(Ω)
(χi)i=1,nh nh
µ
Pµ
∣∣∣∣∣∣∣∣
∫
Ω
∇uµ.∇χi dΩ + µ
∫
Ω
uµχi dΩ =
∫
Ω
fχi dΩ ∀ 1 ≤ i ≤ nh,
uµ =
nh∑
i=1
Uµ,iχi.
Uµ
(A + µM)Uµ = F
Aij =
∫
Ω
∇χi.∇χj dΩ,
Mij =
∫
Ω
χiχj dΩ,
Fi =
∫
Ω
fχi dΩ.
O(nγh)
γ
13
AM
nh
m
(uµk)k=1,m
(Pµ)
µ
uµk
(Pµ)
∣∣∣∣∣∣∣∣
∫
Ω
∇uµ.∇uµk+ µ
∫
Ω
uµuµk=
∫
Ω
fuµk∀ 1 ≤ k ≤ m,
uµ =m∑
k=1
αµ,kuµk.
αµ
(A+ µM)αµ = F
![Page 127: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/127.jpg)
1 ≤ k, l ≤ m
Akl =
∫
Ω
∇uµk.∇uµl
dΩ,
Mkl =
∫
Ω
uµkuµl
dΩ,
Fk =
∫
Ω
fuµkdΩ.
V
M(nh, m)
Uµk
A = V tAV
M = V tMV
F = V tF
O(mγnηh)
uµk
µ
nh
A
M
AM
O(m3) (≪ O(nγ
h)
m≪ nh
(P1
µ)
(P2µ)
(P1µ)
∣∣∣∣−u+ g(µ, x)u = f
Ω,
u = 0
∂Ω.
g(µ, x) 6= f(µ)h(x)
IR
µ
(P2µ)
∣∣∣∣−u+ g(µ, u) = f
Ω,
u = 0
∂Ω.
g
IR
(P1
µ)
µ
g
Gµ
∀ 1 ≤ k, l ≤ m, (Gµ)kl =
∫
Ω
g(µ, x)uµkuµl.
gGµ 6= f(µ)M
µ
Gµ
O(m2nηh)
η = 1, 2
M
nh
(γ, η) = (1, 2)
A
(γ, η) = (2, 1)
A
![Page 128: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/128.jpg)
(P2
µ)
u
uk+1
µ = ukµ + δk
uµ
δkuµ
∣∣∣∣∣∣∣∣∣∣∣∣
∫
Ω
∇δkuµ.∇χi dΩ +
∫
Ω
g′(µ, ukµ)δ
kuµχi dΩ =
∫
Ω
fχi −∫
Ω
∇ukµ.∇χi dΩ
−∫
Ω
g(µ, ukµ)χi dΩ ∀ 1 ≤ i ≤ nh,
δkuµ
=
nh∑
i=1
Uµ,iχi.
Gu
G′u
∀ 1 ≤ k ≤ m, (Gu)k =
∫
Ω
g(µ, uµ)uµk,
∀ 1 ≤ k, l ≤ m, (G′u)kl =
∫
Ω
g′(µ, uµ)uµkuµl.
µ
nh
N
M
(zk)k=1,M
µ = (xk)k=1,M
nh
B1
nh = O(N)
B2
nh >> N
Ω
d
![Page 129: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/129.jpg)
L
L
O(nβh)
β
1
3
N
104
L
105
m
m≪ nh
N
Ψ = (ψi)i=1,N
EKS(Ψ) =
N∑
i=1
∫
Ω
|∇ψi|2 +
∫
Ω
ρVµ +1
2
∫
Ω
∫
Ω
ρ(x)ρ(x′)
|x− x′| dxdx′ + Exc(ρ)
∀ 1 ≤ i ≤ j ≤ N,
∫
Ω
ψiψj = δij
ρ(x) = 2N∑
i=1
|ψi(x)|2,
Vµ(x) = −M∑
k=1
zk
|x− xk|.
Exc
![Page 130: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/130.jpg)
∣∣∣∣∣∣∣∣∣∣∣∣∣∣
∀ 1 ≤ i ≤ N, −1
2ψi + Vµψi +
(ρ ⋆
1
|x|)ψi + vxc(ρ)ψi =
N∑
j=1
λijψj ,
ρ = 2N∑
i=1
|ψi|2,
∀ 1 ≤ i ≤ j ≤ N,
∫
Ω
ψiψj = δij.
vxc(ρ)
Exc(ρ)
ρ
φ
−φ = 4πρ
IRd,
u IRd
∫
IRd
(ρ ⋆
1
|x|)(x)u(x) dx =
∫
IRd
∫
IRd
ρ(x′)u(x)
|x− x′| dx′dx =
∫
IRd
φ(x)u(x) dx.
Ω
−φ = 4πρ
Ω,φ = 0
∂Ω.
Ω
u
Ω ∫
Ω
∫
Ω
ρ(x′)u(x)
|x− x′| dx′dx ≃∫
Ω
φ(x)u(x) dx.
O(n2
h)
O(nh) φ
O(nh)
∣∣∣∣∣∣∣∣∣∣∣∣∣∣
∀ 1 ≤ i ≤ N, −1
2ψi + Vµψi + φψi + vxc(ρ)ψi =
N∑
j=1
λijψj ,
ρ = 2N∑
i=1
|ψi|2, −φ = 4πρ
Ω, φ = 0
∂Ω,
∀ 1 ≤ i ≤ j ≤ N,
∫
Ω
ψiψj = δij.
![Page 131: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/131.jpg)
ψi
ψi
ψi
L2
µ
Vµ
µ
ψi µ
ψi
xk
B1 ηk
µ
B2
ηk
µ
µ
µ
B1
P1
B1
µ
Ω
µ
N
N
O(N)
![Page 132: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/132.jpg)
(P1
µ)
(P2µ)
Bu
B u
Bu
supµ∈Sµ
infα∈IRm
∥∥∥∥u(µ) −m∑
k=1
αku(µk)
∥∥∥∥∞
≤ ǫ
Sµ
µ
µ
ǫ
µ
uµ
(P1
µ)
∣∣∣∣−u+ g(µ, x)u = f
Ω,
u = 0
∂Ω.
Bg =(
g(µk, x))
k
µk
supµ∈Sµ
infγ∈IRK
∥∥∥∥g(µ, x) −K∑
k=1
γkg(µk, x)
∥∥∥∥∞
≤ ǫ.
g(µ, x) =
K∑
k=1
γk(µ)g(µk, x).
γ(µ) ∈ IRK
∀ 1 ≤ i ≤ K, g(µ, xi) =K∑
k=1
γk(µ)g(µk, xi)
(xi)i=1,K
Ω
[
]
r1(x) = g(µ1, x),
x1 = argmax
|r1(x)|
B111 = g(µ1, x1),
2 ≤ k ≤ K
![Page 133: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/133.jpg)
Bk−1δ = [g(µk, x1), . . . , g(µk, xk−1)]
t, δ ∈ IRk−1 rk(x) = g(µk, x) −
k−1∑
j=1
δjg(µj, x) xk = argmax
|rk(x)|
Bk
ij = g(µi, xj), ∀ 1 ≤ i, j ≤ k, Bk ∈ IRk×k.
γ(µ)
BKγ(µ) = [g(µ, x1), . . . , g(µ, xK)]t.
µ BK
xk
AF
(Gk)k=1,K
Gkij =
∫
Ω
g(µk, x)uµiuµj
dΩ.
O(Km2nh)
µ uµ
(A+
K∑
k=1
γk(µ)Gk
)α = F .
α
O(m3 + K2 + Km2)
µ Bu
Bg
uµ
uµ
(P2
µ)
∣∣∣∣−u+ g(µ, u) = f
Ω,
u = 0
∂Ω.
g(µ, uµ) =K∑
k=1
γk(α)g(µk, uµk)
![Page 134: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/134.jpg)
Bg =(g(µk, uµk
))
k
µk
γ
∀ 1 ≤ i ≤ K, g
(µ,
m∑
j=1
αjuµj
)(xi
)=
K∑
k=1
γk(α)g(µk, uµk)(xi)
(xi)i=1,K
γ(α)
BKγ(α) = G(α)
∀ 1 ≤ i ≤ K, G(α)i = g
(µ,
m∑
j=1
αjuµj
)(xi).
µ BK
xk
A
F G
Gik =
∫
Ω
g(µk, u(µk)
)uµi
dΩ.
O(Kmnh)
µ
∣∣∣∣∣Aα + Gγ(α) = F ,
γ(α) =(BK)−1
G(α).
αn
n αn+1 = αn + δα
(A + G
(BK)−1
G′(αn))δα = F − Aαn − Gγ(αn)
∀ 1 ≤ k ≤ K, ∀ 1 ≤ i ≤ m, G′ki(αn) =
∂g
∂u
(µ,
m∑
j=1
αjuµj
)(xk)uµi
(xi).
α
O(m3 +K2 +Km2)
m
K
![Page 135: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/135.jpg)
ψi
Ω
[0, a]
∫Ω
|∇u|2 dΩ − µ
∫
Ω
u4 dΩ,
∫
Ω
u2 dΩ = 1, u ∈ H10 (Ω)
.
(uµ, λµ)
−uµ − 2µu3µ = λµ uµ
Ω,∫
Ω
u2µ dΩ = 1,
uµ ∈ H10 (Ω).
g
g(u) = u3 Bu
Bg ǫ
(xk)k=1,K
Bg BK A M
G
µ
(α, λµ)
Aα− 2µGγ(α) = λµMα,
αtMα = 1.
γ(α)
(5.16) g g(u) = u4 α α
![Page 136: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/136.jpg)
a = 20
200
µ
[0, a]
Sµ
nh
nh
200
ǫ = 10−5
Bu
Bg
Bu = 0; 3.21; 7.44; 14.86; 17.89; 19.90, (Bu
)= 6,
Bg = 0; 2.91; 7.09; 9.55; 11.76; 15.58; 18.39; 19.40; 20, (Bg
)= 9.
Bg (xk)k=1,9
(xk)k =3.8; 6.8; 7.7; 9.2; 10; 11.3; 13.8; 14.6; 17.5
.
µ
eλ eL
eH
eλ = |λµ − λµ|, eL = ‖uµ − uµ‖L2(Ω), eH = |uµ − uµ|H1(Ω).
H1 uµ
L2 uµ
λµ
µ
(λµ, uµ)
(λµ, uµ)
µ
10
m
3
6
K
3
9
eλ eL
eH
(m,K)
µ = 10
m K
m
K
m
K
Km
mK
Bg
Bu
µ = 10
(uref , λref) nh = 2560
nh
(λµ, uµ)
(λµ, umu) eref
λ (.) erefL (.) eref
H (.)
eref
λ (.) = |λref − . |, erefL = ‖uref − . ‖L2(Ω), e
refH = |uref − . |H1(Ω). (λµ, uµ)
(λµ, uµ)
(λµ, uµ)
nh ∈ [80, 2560]
(λµ, uµ)
nh = 80
(λµ, uµ)
nh ∈ [80, 2560]
(λµ, uµ)
nh = 80
nh = 80 (λµ, uµ)
(λµ, uµ)
![Page 137: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/137.jpg)
H+2
Vµ
H+
2
µ (a, b, c)
c > µ/2
a >> µ
Ω = [−a−µ/2, a+µ/2]× [0, a] uµ
1
2
∫
Ω
|∇u|2 dΩ −∫
Ω
Vµu2 dΩ,
∫
Ω
u2 dΩ = 1, u ∈ H10 (Ω)
∀ (x, y) ∈ Ω, Vµ(x, y) =1√
(x+ µ/2)2 + y2+
1√(x− µ/2)2 + y2
.
uµ
(−µ/2, 0)
(µ/2, 0)
H2
fµ
∣∣∣∣∣∣∣∣
f1,µ(x, y) = (x+ µ/2 − b− c, y) (x, y) ∈ [−a− µ/2,−µ/2 + c] × [0, a],
f2,µ(x, y) =( 2b
µ− 2cx, y)
(x, y) ∈ [−µ/2 + c, µ/2 − c] × [0, a],
f3,µ(x, y) = (x− µ/2 + b+ c, y) (x, y) ∈ [µ/2 − c, a+ µ/2] × [0, a]. Ω = fµ(Ω) = [−a− b− c, a+ b+ c]
µ
Ω µ Ω
f−1
µ
uµ
uµ = inf
L(u, µ) − V(u, µ), C(u, µ) = 1, u ∈ H1
0 (Ω)
L(u, µ) =1
2
∫
Ω1
|∇u|2 dΩ1 +1
2
∫
Ω3
|∇u|2 dΩ3
+1
2
∫
Ω2
2b
µ− 2c
( ∂u∂x2
)2
+µ− 2c
2b
( ∂u∂y2
)2
dΩ2,
C(u, µ) =
∫
Ω1
u2 dΩ1 +µ− 2c
2b
∫
Ω2
u2 dΩ2 +
∫
Ω3
u2 dΩ3
![Page 138: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/138.jpg)
H+
2
V(u, µ) =
∫
Ω1
g1(µ, x1, y1)u2 dΩ1 +
∫
Ω3
g3(µ, x3, y3)u2 dΩ3
+µ− 2c
2b
∫
Ω2
g2(µ, x2, y2)u2 dΩ2
g1(µ, x1, y1) =1√
(x1 + b+ c)2 + y21
+1√
(x1 − µ+ b+ c)2 + y21
,
g2(µ, x2, y2) =1√(
µ−2c2bx2 + µ
2
)2
+ y22
+1√(
µ−2c2bx2 − µ
2
)2
+ y22
,
g3(µ, x3, y3) =1√
(x3 + µ− b− c)2 + y23
+1√
(x3 − b− c)2 + y23
.
(uµ, λµ)
Aµuµ = λµMµuµ
utµMµuµ = 1
Aµ = A1 +2b
µ− 2cAx
2 +µ− 2c
2bAy
2 + A3 + Gµ1 +
µ− 2c
2bGµ
2 + Gµ3 ,
Mµ = M1 +µ− 2c
2bM2 + M3.
(A.i)i
(Mi)i
µ
∀ 1 ≤ k ≤ 3,(Gµ
k
)
ij=
∫
Ωk
gk(µ)ηiηj dΩk,
ηj
Fµ
Fµ =∂λµ
∂µ=ut
µAµuµ
∂µ= ut
µKuµ
K = − 2b
(µ− 2c)2Ax
2 +1
2bAy
2 −1
2bGµ
2 + G′µ,1 +
µ− 2c
2bG′
µ,2 + G′µ,3 −
λµ
2bM2,
∀ 1 ≤ k ≤ 3,(G′
µ,k
)ij
=
∫
Ωk
g′k(µ)ηiηj dΩk.
![Page 139: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/139.jpg)
µ gk
Sgk
= (µkl )l=1, Kk
3
(xkl , y
kl )l=1, Kk
µ
(γk)k=1, 3
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ l ≤ Kk, gk(µ, xkl , y
kl ) =
Kk∑
p=1
γkpgk(µ
kp, x
kl , y
kl ).
Vµ
uµ =m∑
i=1
αiuµi α
Aµα = λµMµα
αtMµα = 1
Aµ = A1 +2b
µ− 2cAx
2 +µ− 2c
2bAy
2 + A3 + Gµ1 − µ− 2c
2bGµ
2 − Gµ3 ,
Mµ = M1 +µ− 2c
2bM2 + M3.
(A.i)i (Mi)i
(Gµ
i )i
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ i, j ≤ m,
∣∣∣∣∣
(A.
k
)ij
= utµiA.
kuµj(Mk
)ij
= utµiMkuµj
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ i, j ≤ m,(Gµ
k
)
ij=
Kk∑
l=1
γkl u
tµi
(G
µkl
k
)uµj
.
Fµ =∂λµ
∂µ=∂αtAµα
∂µ= αtKα
K = − 2b
(µ− 2c)2Ax
2 +1
2bAy
2 −1
2bGµ
2 −3∑
k=1
Kk∑
l=1
∂γkl
∂µG
µkl
k − λµ
2bM2.
![Page 140: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/140.jpg)
H2
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ l ≤ Kk,∂γk
l
∂µ=
((Bk)−1dgk
)
l,
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ l, p ≤ Kk, (Bk)lp = gk(µkp, x
kl , y
kl ),
∀ 1 ≤ k ≤ 3, ∀ 1 ≤ l ≤ Kk,(dgk
)l
=∂gk
∂µ(µ, xk
l , ykl ).
a = 20
100
µ
[1, 5]
Sµ
Bu
Bkg
gk
ǫ = 10−5 K
Kk |F (µ)−FN (µ)| |E(µ)−EN(µ)| ‖u(µ)−uN(µ)‖H(Ω)
eF = |Fµ − Fµ| eE = eλ
eH
µ m = 8
K
6
11
eE eF
eu
K
λµ
uµ
(λµ, Fµ, uµ)
K
K
eE eF
eu
µ = 2
µ = 2.6
uµ
λµ K
µ
eE eF
eu
4
11 Bu
µ
[1, 5]
H2
H2
![Page 141: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/141.jpg)
µ
(ax, ay, b, c)
c > µ/2
ax, ay >> µ
Ω =
[−ax − µ/2, ax + µ/2] × [0, ay] uµ
inf
1
2
∫
Ω
|∇u|2 dΩ −∫
Ω
Vµu2 +
∫
Ω
∫
Ω
u2(r)u2(r′)
|r − r′| dΩ dΩ,
∫
Ω
u2 dΩ = 1, u ∈ H10 (Ω)
.
∀ u ∈ L3(IR3),
∫
IR3
∫
IR3
u2(r)u2(r′)
|r − r′| dIR3dIR3 =
∫
IR3
φu2 dIR3
φ ∈ L2(IR3)
−φ = 4πu2 L2(IR3).
Ω
∫
Ω
∫
Ω
u2(r)u2(r′)
|r − r′| dΩ dΩ ≃∫
Ω
φu2 dΩ
−φ = 2u2 Ω,
φ = 0
∂Ω.
Ω a
φ ∂Ω φ
Eh
Ω
(uµ, φµ)
∀v ∈ Eh,1
2
∫
Ω
∇uµ.∇v dΩ −∫
Ω
Vµuµv dΩ + 2
∫
Ω
φµuµv dΩ = λ
∫
Ω
uµv dΩ,
∀w ∈ Eh,
∫
Ω
∇φµ.∇w dΩ = 2
∫
Ω
u2µw dΩ,
∫
Ω
u2µ dΩ = 1.
![Page 142: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/142.jpg)
H2
fµ
(λµ, uµ, φµ
Eh
Ω
λµ, uµ, φµ)
(Pµ)
∀v ∈ Eh, F (uµ, φµ, v) = 0,∀w ∈ Eh, G(φµ, uµ, w) = 0,
C(uµ) = 1.
F (u, φ, v) =1
2
(∫
Ω1
∇u.∇v dΩ1 +
∫
Ω3
∇u.∇v dΩ3
)
+1
2
∫
Ω2
2b
µ− 2c
∂u
∂x.∂v
∂x+µ− 2c
2b
∂u
∂y.∂v
∂ydΩ2
)
−∫
Ω1
g1(µ, x1, y1)uv dΩ1 −µ− 2c
2b
∫
Ω2
g2(µ, x2, y2)uv dΩ2
−∫
Ω3
g3(µ, x3, y3)uv dΩ3
+ 2(∫
Ω1
φuv dΩ1 +µ− 2c
2b
∫
Ω2
φuv dΩ2 +
∫
Ω3
φuv dΩ3
)
− λ(∫
Ω1
uv dΩ1 +µ− 2c
2b
∫
Ω2
uv dΩ2 +
∫
Ω3
uv dΩ3
),
G(φ, u, w) =
∫
Ω1
∇φ.∇w dΩ1 +
∫
Ω3
∇φ.∇w dΩ3
+
∫
Ω2
2b
µ− 2c
∂φ
∂x.∂w
∂x+µ− 2c
2b
∂φ
∂y.∂w
∂ydΩ2
− 2(∫
Ω1
u2w dΩ1 +µ− 2c
2b
∫
Ω2
u2w dΩ2 +
∫
Ω3
u2w dΩ3
),
C(u) =
∫
Ω1
u2 dΩ1 +µ− 2c
2b
∫
Ω2
u2 dΩ2 +
∫
Ω3
u2 dΩ3.
Vµ
H+
2
uµ Bφ
φµ
µ uµ
φµ
uµ =mu∑
i=1
αiuµi,
φµ =
mφ∑
j=1
βjφµj.
![Page 143: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/143.jpg)
αβ
(v, w) = (uµi, φµj
)
1 ≤ i ≤ mu
1 ≤ j ≤ mφ
ax = 40 ay = 8
100
µ
[1, 10]
Sµ
ǫ = 10−5 Bu Bφ
Bkg
gk
µ
g1
g2
m
K Bu
Bφ Bkg
eE = |E(µ) − EN(µ)|
eF = |F (µ)−FN(µ)| eλ = |λ(µ)−λN(µ)| eu = ‖u(µ)−uN(µ)‖H(Ω)
|Eµ − Eµ| |Fµ − Fµ| |λµ − λµ| ‖uµ − uµ‖H1(Ω)
(m,K)
eE eF
eλ eu
µ
K
(λµ, uµ, φµ) m
K
m
(λµ, umu, φmu)
K
m
K
[6, 9]
m
5, 7 Bu
Bφ
Bkg
φ
(mu, mφ)K
(λµ, uµ, φµ)
(λµ, uµ, φµ)
µ
(λµ, uµ, φµ)
2000
(λµ, uµ, φµ)
φ
2000 Ω
2, 3
![Page 144: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/144.jpg)
fµ
IR3
O(N2)
![Page 145: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/145.jpg)
![Page 146: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/146.jpg)
![Page 147: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/147.jpg)
![Page 148: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/148.jpg)
![Page 149: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/149.jpg)
![Page 150: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/150.jpg)
![Page 151: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/151.jpg)
−µ/2 + c
a
µ/2 a + µ/2
a
−b b b + c a + b + c
0
0
Ω2Ω1
Ω1 Ω2
Ω3
Ω3
−µ/2 µ/2 − c
−a − b − c −b − c
−a − µ/2
(x2, y2) = f−12,µ
(x2, y2)(x1, y1) = f−11,µ
(x1, y1) (x3, y3) = f−13,µ
(x3, y3)
f−1
i,µ
H+
2
H2
![Page 152: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/152.jpg)
K = 6 K = 7
K = 8 K = 9
K = 10 K = 11
eE eF
eu
m = 8
K
![Page 153: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/153.jpg)
K = 6 K = 7
K = 8 K = 9
K = 10 K = 11
eE eF eλ
eu
m = 5
K
![Page 154: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/154.jpg)
K = 6 K = 7
K = 8 K = 9
K = 10 K = 11
eE eF eλ
eu
m = 7
K
![Page 155: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/155.jpg)
K = 6 K = 7
K = 8 K = 9
K = 10 K = 11
eE eF eλ
eu
m = 9
K
![Page 156: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/156.jpg)
0 2 4 6 8 10 12 14 16 18 2010
−14
10−12
10−10
10−8
10−6
10−4
10−2
µ
Err
eur
eλe
L
eH
(m,K) = (6, 9)
0 2 4 6 8 10 12 14 16 18 20−1
0
1
2
3
4
5
6
7
8
9
µ
Tem
ps d
e ca
lcul
FEMBRD
![Page 157: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/157.jpg)
3
4
5
6
7
8
9
33.5
44.5
55.5
6
10−10
10−8
10−6
10−4
10−2
K
m
eλ
3
4
5
6
7
8
9
3
3.5
4
4.5
5
5.5
6
10−6
10−5
10−4
10−3
10−2
Km
eL
34
56
78
9
33.5
44.5
55.5
6
10−4
10−3
10−2
10−1
100
Km
eH
(m,K)
µ = 10
![Page 158: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/158.jpg)
0 500 1000 1500 2000 2500 300010
−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
nh
ere
fλ
FEMBRD
0 500 1000 1500 2000 2500 30000
20
40
60
80
100
120
nh
Te
mp
s d
e c
alc
ul
FEMBRD
0 500 1000 1500 2000 2500 300010
−10
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−1
nh
ere
fL
FEMBRD
0 500 1000 1500 2000 2500 300010
−10
10−8
10−6
10−4
10−2
100
102
nh
ere
fH
FEMBRD
µ = 10
![Page 159: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/159.jpg)
uµ
H2
µ
![Page 160: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/160.jpg)
eE eF
eu
K
M
µ = 2
µ = 2.6
eE eF
eu
(m,K) = (8, 8)
![Page 161: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/161.jpg)
u
φ φ g1
(V r1µ )
g2 (V r2
µ )
g1
g2
![Page 162: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/162.jpg)
uµ
uµ
µ
![Page 163: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/163.jpg)
![Page 164: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/164.jpg)
Tz
rzc
k(r) V0
H2
Be2
Be2
![Page 165: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/165.jpg)
H2
S
F
![Page 166: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/166.jpg)
SF
γ
O(Nα
) α < 3
χµ
O(N)
ψl
i
O(N) ψd
i
ψd
i
χµ
χl
µ
ψli
ψd
i
ρd χd
µ
ψdi
χdµ
Cd ψd
i
χd
µ
ρ
ψli
ψd
i
![Page 167: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/167.jpg)
F d Sd = Id
Cd C F
H+
2
H2
(
Nb = O(N))
φi N
N
ρ
![Page 168: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/168.jpg)
![Page 169: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/169.jpg)
![Page 170: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/170.jpg)
![Page 171: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/171.jpg)
O2
![Page 172: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/172.jpg)
![Page 173: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/173.jpg)
![Page 174: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/174.jpg)
![Page 175: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/175.jpg)
/
![Page 176: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/176.jpg)
O(N)
12
![Page 177: Development of fast methods for electronic structure ......HAL Id: pastel-00001655 Submitted on 8 Jun 2006 HAL is a multi-disciplinary open access archive for the deposit](https://reader035.vdocument.in/reader035/viewer/2022071016/5fcea5507b21ab4e557652d3/html5/thumbnails/177.jpg)