1 2 LDA Local Density Approximation 2 LDA 1 LDA LDA N N N H = N [ 2 j + V ion (r j ) ] + 1 e 2 2 r j r k j j k (3) V ion V ion (r) = I Z I e 2 r

Size: px
Start display at page:

Download "1 2 LDA Local Density Approximation 2 LDA 1 LDA LDA N N N H = N [ 2 j + V ion (r j ) ] + 1 e 2 2 r j r k j j k (3) V ion V ion (r) = I Z I e 2 r"

Transcription

1 11 March [ { } ] 3 1/3 2 + V ion (r) + V H (r) 3α 4π ρ σ(r) ϕ iσ (r) = ε iσ ϕ iσ (r) (1) KS Kohn-Sham [ 2 + V ion (r) + V H (r) + V σ xc(r) ] ϕ iσ (r) = ε iσ ϕ iσ (r) (2)

2 1 2 LDA Local Density Approximation 2 LDA 1 LDA LDA N N N H = N [ 2 j + V ion (r j ) ] + 1 e 2 2 r j r k j j k (3) V ion V ion (r) = I Z I e 2 r j R I (4) Z I N Ψ HΨ = EΨ (5) e Ψ N r j Ψ = Ψ (r 1,, r N ) (6) 2

3 Ψ (r 1... r N ) = ϕ 1 (r 1 ),, ϕ N (r N ) (7) 7 N 3 [ 2 + V ion (r) + V H (r) ] σ =σ ϕ iσ (r) j ˆF jσ,iσ(r)ϕ jσ (r) = ε iσ ϕ iσ (r) (8) ˆF jσ,iσ(r) r 8 i It is now perfectly clear that a single configuration (Slater determinant) wave function must inevitably lead to a poor energy Coulson (1960) N N Ψ N K = {ϕ 1 ϕ N } Ψ N = K c K Ψ(K) (9) 3 3

4 CI Configuration Interaction CI CI α 2/3 α Xα Hohenberg-Kohn N Hohenberg-Kohn HK [1] N v 1 ( ) E G ρ(r) ρ(r) v 2 ( ) E G [ρ] N ρ (r) ρ (r) ρ(r) E G [ρ(r)] < E G [ρ (r)] (10) Hohenberg-Kohn 4

5 ψ ~ e ikr E = h2 k 2 2m ρ ~ const 1: E v [Ψ] Ψ Hohenberg-Kohn HK Ψ Ψ N 7 N N N 1 Ψ ρ HK ρ Ψ 1 Kohn 5

6 1 ψ k = e ikr k k k 2 ρ(r) = ψ(r) 2 k HK ψ ρ v ρ v ψ (11) HK ρ ψ v Hohenberg-Kohn ρ v v ρ v ψ ρ (12) 12 v ψ ρ 12 v ψ HK Hohenberg-Kohn 12 ρ v (13) 6

7 ρ ρ ψ v v 2: ρ v ρ v v ρ Hohenberg-Kohn v ρ E G [ρ] E G [ρ] = F HK [ρ] + ρ(r)v(r)dr (14) Hohenberg- Kohn F HK [ρ] F HK [ρ] N Ψ F HK [ρ] = Ψ ˆT + ˆV ee Ψ F HK [ρ] ρ v ρ v Löwdin A. J. Coleman, (15)

8 Q ρ ψ 3: Levy ρ v n Hohenberg-Kohn 3.2 ρ n ρ HK Hohenberg Kohn Sham ρ v ρ v Levy n [3] 3 v ρ ρ v F HK [ρ] v ρ ψ ρ 8

9 HK ψ v ψ ρ Levy Q[ρ] = min Ψ ρ ˆT + ˆVee Ψρ Ψ ρ Q[ρ] HK 16 ρ ρ Ψ ρ ˆT + ˆV ee Ψ ρ Q[ρ] Q[ρ] ρvdr E[ρ] ρ v Q[ρ] 15 F HK [ρ] v 16 ρ v HK 3 ρ (16) 3.3 Thoms-Fermi HK N 15 N E[ρ] = T [ρ] + ρ(r)v(r)dr + U ee [ρ] (17) U ee [ρ] U ee [ρ] = 1 ρ(r)v H (r)dr + E xc [ρ] (18) 2 V H ρ(r ) V H (r) = r r dr (19) U ee Thoms-Fermi T [ρ] = C TF dr [ρ(r)] 5/3 (20) 9

10 5 3 C TFρ(r) 2/3 + v(r) + ρ(r ) r r dr = µ (21) µ N 21 ρ(r) v(r) Thoms-Fermi Thoms-Fermi Kohn-Sham [2] ρ(r) N {ϕ i (r)} ρ(r) = N ϕ i (r) 2 (22) i 22 2 T T s [ρ] = N ϕi 2 ϕ i i (23) 23 N E xc T s [ρ] T E xc 2 v xc (r) = δe xc[ρ] δρ 2 ρ(r) Ψ 22 ρ(r) (24) 10

11 ρ N N 2 ε i E = N ε i 1 2 i=1 ρ(r)v H (r)dr + E xc [ρ] ρ(r)v xc (r)dr (25) N 2 22 ϕ i (r) 22 {ϕ i (r)} ρ(r) {ϕ i (r)} 4.2 E xc [ρ] E xc [ρ] E xc [ρ] = Exc HOM (ρ) E xc [ρ] LDA ɛ xc (ρ(r)) [] () 11

12 LDA LDA LDA [9] LDA E xc [ρ] E xc [ρ] LDA LDA E xc [ρ] It is advisable to stop at the simple LDA W. Kohn (1984) LDA Exc HOM (ρ) Exc HOM (ρ) LDA LDA LDA 4.3 KS KS KS 12

13 i i ε i I i i E(n 1,, n i,, n N ) E(n 1,, n i 1,, n N ) I i = E(, n i, ) E(, n i 1, ) (26) I (1) = I N I (2) N 0 E(N) = N I (i) (27) N 27 E = i ε i I i = ε i 27 I (i) = ɛ N+1 i 25 I i = ε i I (i) = ɛ N+1 i N i=1 E(, n i, ) E(, n i 1, ) = ε i (28) 28 E n i = ε i (29) 13

14 Janak [4] E i E(, n i, ) n i E(, n i,, n N ) E(n 1,, n i 1, ) ε i (, n i 0.5, ) (30) KS ε N 0 < n N 1 n N [5] E(N) E(N 1) = ε N (31) E(N) = N ε i (i) (32) i=1 N ɛ i = ɛ i (N)

15 1: HF DFT Ψ(x 1,, x N ) ρ(r) Ψ ρ ρ Ψ Ψ Ψ H Ψ E[ρ] LDA HF(S) KS [ 2 + V ion (r) + V H (r) + V Xα (r) ] [ 2 + V ion (r) + V H (r) + V xc (r) ] ϕ i (r) = ε i ϕ i (r) ϕ i (r) = ε i ϕ i (r) CI CI N LDA N ρ(r) ρ (r) ρ (r) ρ (r) ρ (r) 2 2 ρ αβ (r) LDA LSD 15

16 [1] P. Hohenberg and W. Kohn, Phys. Rev. 136 B864 (1964). [2] W. Kohn and L. Sham, Phys. Rev. 140 A1133 (1965). [3] M. Levy, Proc. Natl. Acad. Sci. USA (1979). [4] J. F. Janak, Phys. Rev. B (1978). [5] J. P. Perdew, R. G. Parr, M. Levy, and J. L. Balduz, Jr., Phys. Rev. Lett (1982). [6] O. Gunnarsson, M. Jonson, and B. I. Lundqvist, Phys. Rev. B 20, 3136 (1979). [7] S. Lundqvist and N. H. March, eds., Theory of the Inhomogeneous Electron Gas (Plenum, New York, 1983). [8] J. Callaway and N. H. March, Solid State Physics 38 (Academic, New York, 1984) p [9] R. O. Jones and O. Gunnarsson, Rev. Mod. Phys. 61, 689 (1989) [10] S. B. Trickey ed., Adv. in Quantum Chemistry 21, (Academic, San Diego, 1989) [11] R. G. Parr and W. Yang, Density-Functional Theory of Atoms and Molecules (Oxford, New York, 1989). [12] C. Fiolhais, F. Nogueira, M. Marques, eds., A Primer in Density Functional Theory (Springer, Berlin, 2003). 16

1: Sheldon L. Glashow (Ouroboros) [1] 1 v(r) u(r, r ) ( e 2 / r r ) H 2 [2] H = ( dr ψ σ + (r) 1 2 ) σ 2m r 2 + v(r) µ ψ σ (r) + 1 dr dr ψ σ + (r)ψ +

1: Sheldon L. Glashow (Ouroboros) [1] 1 v(r) u(r, r ) ( e 2 / r r ) H 2 [2] H = ( dr ψ σ + (r) 1 2 ) σ 2m r 2 + v(r) µ ψ σ (r) + 1 dr dr ψ σ + (r)ψ + 1 1.1 21 11 22 10 33 cm 10 29 cm 60 6 8 10 12 cm 1cm 1 1.2 2 1 1 1: Sheldon L. Glashow (Ouroboros) [1] 1 v(r) u(r, r ) ( e 2 / r r ) H 2 [2] H = ( dr ψ σ + (r) 1 2 ) σ 2m r 2 + v(r) µ ψ σ (r) + 1 dr dr

More information

講 座 熱電研究のための第一原理計算入門 第1回 密度汎関数法による第一原理バンド計算 桂 1 はじめに ゆかり 東京大学 2 密度汎関数理論 第一原理 first-principles バンド計算とは 結晶構造 Schrödinger 方程式は 量子力学を司る基本方程式で 以外の経験的パラメータや

講 座 熱電研究のための第一原理計算入門 第1回 密度汎関数法による第一原理バンド計算 桂 1 はじめに ゆかり 東京大学 2 密度汎関数理論 第一原理 first-principles バンド計算とは 結晶構造 Schrödinger 方程式は 量子力学を司る基本方程式で 以外の経験的パラメータや 講 座 熱電研究のための第一原理計算入門 第1回 密度汎関数法による第一原理バンド計算 桂 1 はじめに ゆかり 東京大学 2 密度汎関数理論 第一原理 first-principles バンド計算とは 結晶構造 Schrödinger 方程式は 量子力学を司る基本方程式で 以外の経験的パラメータや任意パラメータを使わず 基 ある 定常状態において電子 i の状態を定義する波動 本的な物理方程式のみを用いて行う電子状態計算であ

More information

和佐田P indd

和佐田P indd 2000 B3LYP/6-31G Gaussian 98 03 B3LYP/6-31G* Gaussian STO-3G RHF Gaussian RHF/STO-3G B3LYP RHF 6-31G* STO-3G Schrödinger Schrödinger s p d Schrödinger Schrödinger Hohenberg-Kohn Kohn-Sham Kohn-Sham [1-3]

More information

TQFT_yokota

TQFT_yokota , TY, Naito, Phys. Rev. B 99, 115106 (2019),, 2019 9 2 1 (DFT) (DFT)? HΨ(x 1,, x N ) = EΨ(x 1,, x N ) N DFT! Hohenberg, Kohn, PR (1964) Kohn, Sham, PRA (1965) (EDF) E[ρ] = F[ρ] + dxv(x)ρ(x) δe[ρ] δρ(x)

More information

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising

,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising ,, Andrej Gendiar (Density Matrix Renormalization Group, DMRG) 1 10 S.R. White [1, 2] 2 DMRG ( ) [3, 2] DMRG Baxter [4, 5] 2 Ising 2 1 Ising 1 1 Ising Model 1 Ising 1 Ising Model N Ising (σ i = ±1) (Free

More information

研修コーナー

研修コーナー l l l l l l l l l l l α α β l µ l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l

More information

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) * * 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *1 2004 1 1 ( ) ( ) 1.1 140 MeV 1.2 ( ) ( ) 1.3 2.6 10 8 s 7.6 10 17 s? Λ 2.5 10 10 s 6 10 24 s 1.4 ( m

More information

i

i 009 I 1 8 5 i 0 1 0.1..................................... 1 0.................................................. 1 0.3................................. 0.4........................................... 3

More information

( ) ) AGD 2) 7) 1

( ) ) AGD 2) 7) 1 ( 9 5 6 ) ) AGD ) 7) S. ψ (r, t) ψ(r, t) (r, t) Ĥ ψ(r, t) = e iĥt/ħ ψ(r, )e iĥt/ħ ˆn(r, t) = ψ (r, t)ψ(r, t) () : ψ(r, t)ψ (r, t) ψ (r, t)ψ(r, t) = δ(r r ) () ψ(r, t)ψ(r, t) ψ(r, t)ψ(r, t) = (3) ψ (r,

More information

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H 199 1 1 199 1 1. Vx) m e V cos x π x π Vx) = x < π, x > π V i) x = Vx) V 1 x /)) n n d f dξ ξ d f dξ + n f = H n ξ) ii) H n ξ) = 1) n expξ ) dn dξ n exp ξ )) H n ξ)h m ξ) exp ξ )dξ = π n n!δ n,m x = Vx)

More information

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' =

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' = y x = α + β + ε =,, ε V( ε) = E( ε ) = σ α $ $ β w ( 0) σ = w σ σ y α x ε = + β + w w w w ε / w ( w y x α β ) = α$ $ W = yw βwxw $β = W ( W) ( W)( W) w x x w x x y y = = x W y W x y x y xw = y W = w w

More information

Aharonov-Bohm(AB) S 0 1/ 2 1/ 2 S t = 1/ 2 1/2 1/2 1/, (12.1) 2 1/2 1/2 *1 AB ( ) 0 e iθ AB S AB = e iθ, AB 0 θ 2π ϕ = e ϕ (ϕ ) ϕ

Aharonov-Bohm(AB) S 0 1/ 2 1/ 2 S t = 1/ 2 1/2 1/2 1/, (12.1) 2 1/2 1/2 *1 AB ( ) 0 e iθ AB S AB = e iθ, AB 0 θ 2π ϕ = e ϕ (ϕ ) ϕ 1 13 6 8 3.6.3 - Aharonov-BohmAB) S 1/ 1/ S t = 1/ 1/ 1/ 1/, 1.1) 1/ 1/ *1 AB ) e iθ AB S AB = e iθ, AB θ π ϕ = e ϕ ϕ ) ϕ 1.) S S ) e iθ S w = e iθ 1.3) θ θ AB??) S t = 4 sin θ 1 + e iθ AB e iθ AB + e

More information

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1)

φ 4 Minimal subtraction scheme 2-loop ε 2008 (University of Tokyo) (Atsuo Kuniba) version 21/Apr/ Formulas Γ( n + ɛ) = ( 1)n (1 n! ɛ + ψ(n + 1) φ 4 Minimal subtraction scheme 2-loop ε 28 University of Tokyo Atsuo Kuniba version 2/Apr/28 Formulas Γ n + ɛ = n n! ɛ + ψn + + Oɛ n =,, 2, ψn + = + 2 + + γ, 2 n ψ = γ =.5772... Euler const, log + ax x

More information

SPring-8_seminar_

SPring-8_seminar_ X 21 SPring-8 XAFS 2016 (= ) X PC cluster Synchrotron TEM-EELS XAFS / EELS HΨ k = E k Ψ k XANES/ELNES DFT ( + ) () WIEN2k, Elk, OLCAO () CASTEP, QUANTUM ESPRESSO FEFF, GNXAS, etc. Bethe-Salpeter (BSE)

More information

N cos s s cos ψ e e e e 3 3 e e 3 e 3 e

N cos s s cos ψ e e e e 3 3 e e 3 e 3 e 3 3 5 5 5 3 3 7 5 33 5 33 9 5 8 > e > f U f U u u > u ue u e u ue u ue u e u e u u e u u e u N cos s s cos ψ e e e e 3 3 e e 3 e 3 e 3 > A A > A E A f A A f A [ ] f A A e > > A e[ ] > f A E A < < f ; >

More information

( ) 1 1.1? ( ) ( ) ( ) 1.1(a) T m ( ) 1.1(a) T g ( ) T g T g 500 74% ( ) T K ( 1.1(b) 15 T g T g 10 13 T g T g T g [ ] A ( ) exp (1.1) T T 0 Vogel-Fulcher T 0 T 0 T K T K Ortho-Terphenil (OTP) SiO 2 (1.1)

More information

2 G(k) e ikx = (ik) n x n n! n=0 (k ) ( ) X n = ( i) n n k n G(k) k=0 F (k) ln G(k) = ln e ikx n κ n F (k) = F (k) (ik) n n= n! κ n κ n = ( i) n n k n

2 G(k) e ikx = (ik) n x n n! n=0 (k ) ( ) X n = ( i) n n k n G(k) k=0 F (k) ln G(k) = ln e ikx n κ n F (k) = F (k) (ik) n n= n! κ n κ n = ( i) n n k n . X {x, x 2, x 3,... x n } X X {, 2, 3, 4, 5, 6} X x i P i. 0 P i 2. n P i = 3. P (i ω) = i ω P i P 3 {x, x 2, x 3,... x n } ω P i = 6 X f(x) f(x) X n n f(x i )P i n x n i P i X n 2 G(k) e ikx = (ik) n

More information

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2 2005 9/8-11 2 2.2 ( 2-5) γ ( ) γ cos θ 2πr πρhr 2 g h = 2γ cos θ ρgr (2.1) γ = ρgrh (2.2) 2 cos θ θ cos θ = 1 (2.2) γ = 1 ρgrh (2.) 2 2. p p ρgh p ( ) p p = p ρgh (2.) h p p = 2γ r 1 1 (Berry,1975) 2-6

More information

( ) ) ) ) 5) 1 J = σe 2 6) ) 9) 1955 Statistical-Mechanical Theory of Irreversible Processes )

( ) ) ) ) 5) 1 J = σe 2 6) ) 9) 1955 Statistical-Mechanical Theory of Irreversible Processes ) ( 3 7 4 ) 2 2 ) 8 2 954 2) 955 3) 5) J = σe 2 6) 955 7) 9) 955 Statistical-Mechanical Theory of Irreversible Processes 957 ) 3 4 2 A B H (t) = Ae iωt B(t) = B(ω)e iωt B(ω) = [ Φ R (ω) Φ R () ] iω Φ R (t)

More information

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n 003...............................3 Debye................. 3.4................ 3 3 3 3. Larmor Cyclotron... 3 3................ 4 3.3.......... 4 3.3............ 4 3.3...... 4 3.3.3............ 5 3.4.........

More information

DVIOUT-fujin

DVIOUT-fujin 2005 Limit Distribution of Quantum Walks and Weyl Equation 2006 3 2 1 2 2 4 2.1...................... 4 2.2......................... 5 2.3..................... 6 3 8 3.1........... 8 3.2..........................

More information

(a) (b) (c) (d) 1: (a) (b) (c) (d) (a) (b) (c) 2: (a) (b) (c) 1(b) [1 10] 1 degree k n(k) walk path 4

(a) (b) (c) (d) 1: (a) (b) (c) (d) (a) (b) (c) 2: (a) (b) (c) 1(b) [1 10] 1 degree k n(k) walk path 4 1 vertex edge 1(a) 1(b) 1(c) 1(d) 2 (a) (b) (c) (d) 1: (a) (b) (c) (d) 1 2 6 1 2 6 1 2 6 3 5 3 5 3 5 4 4 (a) (b) (c) 2: (a) (b) (c) 1(b) [1 10] 1 degree k n(k) walk path 4 1: Zachary [11] [12] [13] World-Wide

More information

液晶の物理1:連続体理論(弾性,粘性)

液晶の物理1:連続体理論(弾性,粘性) The Physics of Liquid Crystals P. G. de Gennes and J. Prost (Oxford University Press, 1993) Liquid crystals are beautiful and mysterious; I am fond of them for both reasons. My hope is that some readers

More information

a L = Ψ éiγ c pa qaa mc ù êë ( - )- úû Ψ 1 Ψ 4 γ a a 0, 1,, 3 {γ a, γ b } η ab æi O ö æo ö β, σ = ço I α = è - ø çèσ O ø γ 0 x iβ γ i x iβα i

a L = Ψ éiγ c pa qaa mc ù êë ( - )- úû Ψ 1 Ψ 4 γ a a 0, 1,, 3 {γ a, γ b } η ab æi O ö æo ö β, σ = ço I α = è - ø çèσ O ø γ 0 x iβ γ i x iβα i 解説 4 matsuo.mamoru jaea.go.jp 4 eizi imr.tohoku.ac.jp 4 maekawa.sadamichi jaea.go.jp i ii iii i Gd Tb Dy g khz Pt ii iii Keywords vierbein 3 dreibein 4 vielbein torsion JST-ERATO 1 017 1. 1..1 a L = Ψ

More information

2 2 1?? 2 1 1, 2 1, 2 1, 2, 3,... 1, 2 1, 3? , 2 2, 3? k, l m, n k, l m, n kn > ml...? 2 m, n n m

2 2 1?? 2 1 1, 2 1, 2 1, 2, 3,... 1, 2 1, 3? , 2 2, 3? k, l m, n k, l m, n kn > ml...? 2 m, n n m 2009 IA I 22, 23, 24, 25, 26, 27 4 21 1 1 2 1! 4, 5 1? 50 1 2 1 1 2 1 4 2 2 2 1?? 2 1 1, 2 1, 2 1, 2, 3,... 1, 2 1, 3? 2 1 3 1 2 1 1, 2 2, 3? 2 1 3 2 3 2 k, l m, n k, l m, n kn > ml...? 2 m, n n m 3 2

More information

4/15 No.

4/15 No. 4/15 No. 1 4/15 No. 4/15 No. 3 Particle of mass m moving in a potential V(r) V(r) m i ψ t = m ψ(r,t)+v(r)ψ(r,t) ψ(r,t) = ϕ(r)e iωt ψ(r,t) Wave function steady state m ϕ(r)+v(r)ϕ(r) = εϕ(r) Eigenvalue problem

More information

80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t = 0 i r 0 t(> 0) j r 0 + r < δ(r 0 x i (0))δ(r 0 + r x j (t)) > (4.2) r r 0 G i j (r, t) dr 0

80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t = 0 i r 0 t(> 0) j r 0 + r < δ(r 0 x i (0))δ(r 0 + r x j (t)) > (4.2) r r 0 G i j (r, t) dr 0 79 4 4.1 4.1.1 x i (t) x j (t) O O r 0 + r r r 0 x i (0) r 0 x i (0) 4.1 L. van. Hove 1954 space-time correlation function V N 4.1 ρ 0 = N/V i t 80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t

More information

橡超弦理論はブラックホールの謎を解けるか?

橡超弦理論はブラックホールの謎を解けるか? 1999 3 (Can String Theory Solve the Puzzles of Black Holes?) 305-0801 1-1 [email protected] D-brane 1 Schwarzschild 60 80 2 [1] 1 1 1 2 2 [2] 25 2.2 2 2.1 [7,8] Schwarzschild 2GM/c 2 Schwarzschild

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 3.............................. 3.............................. 4.3 4................... 5.4........................ 6.5........................ 8.6...........................7

More information

Kroneher Levi-Civita 1 i = j δ i j = i j 1 if i jk is an even permutation of 1,2,3. ε i jk = 1 if i jk is an odd permutation of 1,2,3. otherwise. 3 4

Kroneher Levi-Civita 1 i = j δ i j = i j 1 if i jk is an even permutation of 1,2,3. ε i jk = 1 if i jk is an odd permutation of 1,2,3. otherwise. 3 4 [2642 ] Yuji Chinone 1 1-1 ρ t + j = 1 1-1 V S ds ds Eq.1 ρ t + j dv = ρ t dv = t V V V ρdv = Q t Q V jdv = j ds V ds V I Q t + j ds = ; S S [ Q t ] + I = Eq.1 2 2 Kroneher Levi-Civita 1 i = j δ i j =

More information

1 911 9001030 9:00 A B C D E F G H I J K L M 1A0900 1B0900 1C0900 1D0900 1E0900 1F0900 1G0900 1H0900 1I0900 1J0900 1K0900 1L0900 1M0900 9:15 1A0915 1B0915 1C0915 1D0915 1E0915 1F0915 1G0915 1H0915 1I0915

More information

128 3 II S 1, S 2 Φ 1, Φ 2 Φ 1 = { B( r) n( r)}ds S 1 Φ 2 = { B( r) n( r)}ds (3.3) S 2 S S 1 +S 2 { B( r) n( r)}ds = 0 (3.4) S 1, S 2 { B( r) n( r)}ds

128 3 II S 1, S 2 Φ 1, Φ 2 Φ 1 = { B( r) n( r)}ds S 1 Φ 2 = { B( r) n( r)}ds (3.3) S 2 S S 1 +S 2 { B( r) n( r)}ds = 0 (3.4) S 1, S 2 { B( r) n( r)}ds 127 3 II 3.1 3.1.1 Φ(t) ϕ em = dφ dt (3.1) B( r) Φ = { B( r) n( r)}ds (3.2) S S n( r) Φ 128 3 II S 1, S 2 Φ 1, Φ 2 Φ 1 = { B( r) n( r)}ds S 1 Φ 2 = { B( r) n( r)}ds (3.3) S 2 S S 1 +S 2 { B( r) n( r)}ds

More information

第5章 偏微分方程式の境界値問題

第5章 偏微分方程式の境界値問題 October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ

More information

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100 positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) 0.5 1.5MeV : thermalization 10 100 m psec 100psec nsec E total = 2mc 2 + E e + + E e Ee+ Ee-c mc

More information

Einstein 1905 Lorentz Maxwell c E p E 2 (pc) 2 = m 2 c 4 (7.1) m E ( ) E p µ =(p 0,p 1,p 2,p 3 )=(p 0, p )= c, p (7.2) x µ =(x 0,x 1,x 2,x

Einstein 1905 Lorentz Maxwell c E p E 2 (pc) 2 = m 2 c 4 (7.1) m E ( ) E p µ =(p 0,p 1,p 2,p 3 )=(p 0, p )= c, p (7.2) x µ =(x 0,x 1,x 2,x 7 7.1 7.1.1 Einstein 1905 Lorentz Maxwell c E p E 2 (pc) 2 = m 2 c 4 (7.1) m E ( ) E p µ =(p 0,p 1,p 2,p 3 )=(p 0, p )= c, p (7.2) x µ =(x 0,x 1,x 2,x 3 )=(x 0, x )=(ct, x ) (7.3) E/c ct K = E mc 2 (7.4)

More information

Microsoft Word - 九大物理集中講義テキスト'13.doc

Microsoft Word - 九大物理集中講義テキスト'13.doc [email protected] 23 2 23 . DFT DFT Hohenberg Kohn DFT DFT 999 DFT DFT DFT 964 J.. ercus [,2] DFT ercus Hyper-netted chan HC [2,3] ercus DFT DFT 2 DFT. DFT + ercus 2. DFT 2. DFT DFT. 2. DFT 3. ercus 2 2 3 3 ercus

More information

MS#sugaku(ver.2).dvi

MS#sugaku(ver.2).dvi 1 1,,,.,,,,,.,.,,,.,, Computer-Aided Design).,,, Boltzmann,, [1]., Anderson, []., Anderson, Schrödinger [3],[4]., nm,,.,,,.,, Schrödinger.,, [5],[6].,,.,,, 1 .,, -. -, -., -,,,. Wigner-Boltzmann Schrödinger,

More information

スライド 1

スライド 1 Matsuura Laboratory SiC SiC 13 2004 10 21 22 H-SiC ( C-SiC HOY Matsuura Laboratory n E C E D ( E F E T Matsuura Laboratory Matsuura Laboratory DLTS Osaka Electro-Communication University Unoped n 3C-SiC

More information

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T NHK 204 2 0 203 2 24 ( ) 7 00 7 50 203 2 25 ( ) 7 00 7 50 203 2 26 ( ) 7 00 7 50 203 2 27 ( ) 7 00 7 50 I. ( ν R n 2 ) m 2 n m, R = e 2 8πε 0 hca B =.09737 0 7 m ( ν = ) λ a B = 4πε 0ħ 2 m e e 2 = 5.2977

More information

006 11 8 0 3 1 5 1.1..................... 5 1......................... 6 1.3.................... 6 1.4.................. 8 1.5................... 8 1.6................... 10 1.6.1......................

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 1 19 3 19.1................... 3 19.............................. 4 19.3............................... 6 19.4.............................. 8 19.5.............................

More information

18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α

18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α 18 I ( ) (1) I-1,I-2,I-3 (2) (3) I-1 ( ) (100 ) θ ϕ θ ϕ m m l l θ ϕ θ ϕ 2 g (1) (2) 0 (3) θ ϕ (4) (3) θ(t) = A 1 cos(ω 1 t + α 1 ) + A 2 cos(ω 2 t + α 2 ), ϕ(t) = B 1 cos(ω 1 t + α 1 ) + B 2 cos(ω 2 t

More information

II ( ) (7/31) II ( [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re

II ( ) (7/31) II (  [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re II 29 7 29-7-27 ( ) (7/31) II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I Euler Navier

More information

Bethe-Bloch Bethe-Bloch (stopping range) Bethe-Bloch FNAL (Fermi National Accelerator Laboratory) - (SciBooNE ) SciBooNE Bethe-Bloch FNAL - (SciBooNE

Bethe-Bloch Bethe-Bloch (stopping range) Bethe-Bloch FNAL (Fermi National Accelerator Laboratory) - (SciBooNE ) SciBooNE Bethe-Bloch FNAL - (SciBooNE 21 2 27 Bethe-Bloch Bethe-Bloch (stopping range) Bethe-Bloch FNAL (Fermi National Accelerator Laboratory) - (SciBooNE ) SciBooNE Bethe-Bloch FNAL - (SciBooNE ) Bethe-Bloch 1 0.1..............................

More information

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 2 版 1 刷発行時のものです. 医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987

More information

2.5 (Gauss) (flux) v(r)( ) S n S v n v n (1) v n S = v n S = v S, n S S. n n S v S v Minoru TANAKA (Osaka Univ.) I(2012), Sec p. 1/30

2.5 (Gauss) (flux) v(r)( ) S n S v n v n (1) v n S = v n S = v S, n S S. n n S v S v Minoru TANAKA (Osaka Univ.) I(2012), Sec p. 1/30 2.5 (Gauss) 2.5.1 (flux) v(r)( ) n v n v n (1) v n = v n = v, n. n n v v I(2012), ec. 2. 5 p. 1/30 i (2) lim v(r i ) i = v(r) d. i 0 i (flux) I(2012), ec. 2. 5 p. 2/30 2.5.2 ( ) ( ) q 1 r 2 E 2 q r 1 E

More information

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 =

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 = 3 3.1 3.1.1 kg m s J = kg m 2 s 2 MeV MeV [1] 1MeV=1 6 ev = 1.62 176 462 (63) 1 13 J (3.1) [1] 1MeV/c 2 =1.782 661 731 (7) 1 3 kg (3.2) c =1 MeV (atomic mass unit) 12 C u = 1 12 M(12 C) (3.3) 41 42 3 u

More information

Hanbury-Brown Twiss (ver. 2.0) van Cittert - Zernike mutual coherence

Hanbury-Brown Twiss (ver. 2.0) van Cittert - Zernike mutual coherence Hanbury-Brown Twiss (ver. 2.) 25 4 4 1 2 2 2 2.1 van Cittert - Zernike..................................... 2 2.2 mutual coherence................................. 4 3 Hanbury-Brown Twiss ( ) 5 3.1............................................

More information

2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h)

2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h) 1 16 10 5 1 2 2.1 a a a 1 1 1 2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h) 4 2 3 4 2 5 2.4 x y (x,y) l a x = l cot h cos a, (3) y = l cot h sin a (4) h a

More information

QMI_10.dvi

QMI_10.dvi ... black body radiation black body black body radiation Gustav Kirchhoff 859 895 W. Wien O.R. Lummer cavity radiation ν ν +dν f T (ν) f T (ν)dν = 8πν2 c 3 kt dν (Rayleigh Jeans) (.) f T (ν) spectral energy

More information

1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l

1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l 1 1 ϕ ϕ ϕ S F F = ϕ (1) S 1: F 1 1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l : l r δr θ πrδr δf (1) (5) δf = ϕ πrδr

More information

meiji_resume_1.PDF

meiji_resume_1.PDF β β β (q 1,q,..., q n ; p 1, p,..., p n ) H(q 1,q,..., q n ; p 1, p,..., p n ) Hψ = εψ ε k = k +1/ ε k = k(k 1) (x, y, z; p x, p y, p z ) (r; p r ), (θ; p θ ), (ϕ; p ϕ ) ε k = 1/ k p i dq i E total = E

More information

: 2005 ( ρ t +dv j =0 r m m r = e E( r +e r B( r T 208 T = d E j 207 ρ t = = = e t δ( r r (t e r r δ( r r (t e r ( r δ( r r (t dv j =

: 2005 ( ρ t +dv j =0 r m m r = e E( r +e r B( r T 208 T = d E j 207 ρ t = = = e t δ( r r (t e r r δ( r r (t e r ( r δ( r r (t dv j = 72 Maxwell. Maxwell e r ( =,,N Maxwell rot E + B t = 0 rot H D t = j dv D = ρ dv B = 0 D = ɛ 0 E H = μ 0 B ρ( r = j( r = N e δ( r r = N e r δ( r r = : 2005 ( 2006.8.22 73 207 ρ t +dv j =0 r m m r = e E(

More information