• 検索結果がありません。

Expanding Polyhedral Universe in Regge Calculus

N/A
N/A
Protected

Academic year: 2021

シェア "Expanding Polyhedral Universe in Regge Calculus"

Copied!
111
0
0

読み込み中.... (全文を見る)

全文

(1)

24 June 2017

Expanding Polyhedral Universe in Regge Calculus

津田 廉

茨城大学 素粒子論研究室

共同研究者:藤原高徳(茨城大学)

核物理

×

物性セミナー

@

千葉工業大学 新習志野キャンパス

(2)

Contents

1. 単体分割

2. Regge Calculus における Hilbert 作用

3. Regge Calculus における Einstein 方程式 4. Expanding Polyhedral Universe

5. Summary

2

(3)

1. 単体分割

3

(4)

2 次元曲面の三角分割

滑らかな曲面

4

(5)

2 次元曲面の三角分割

三角形で近似

5

(6)

曲率 · · · 三角形の頂点に定義

平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)

6

(7)

曲率 · · · 三角形の頂点に定義

平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)

6

(8)

曲率 · · · 三角形の頂点に定義

過不足なし

平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)

不足 過剰

6

(9)

n- 単体 · · · n + 1 個の頂点を持つ n 次元凸面体

0- 単体

1- 単体

2- 単体

3- 単体

  .. .

 点

.. .

線分 三角形 四面体

  .. .

7

(10)

単体の構造

=

2- 単体

× 3 × 3

=

3- 単体

× 4 × 6 × 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

(11)

単体の構造

2- 単体 1- 単体

× 3 0- 単体

× 3

=

,

3- 単体

× 4 × 6 × 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

(12)

単体の構造

2- 単体 1- 単体

× 3 0- 単体

× 3

, ,

,

3- 単体 2- 単体

× 4 1- 単体

× 6 0- 単体

× 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

(13)

単体の構造

Volume : A

b

A

h

名称 : Boundary Hinge: h

n- 単体 ∋ { (n 1) - 単体 , (n 2) - 単体 , · · · , 0- 単体 }

×

n+1

C

1

×

n+1

C

2

×

n+1

C

n

9

(14)

Delaunay 格子と Voronoi ポリゴン

単体分割によって作られた格子 · · · Delaunay 格子

1. 各単体の外心を求める

2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン

10

(15)

Delaunay 格子と Voronoi ポリゴン

単体分割によって作られた格子 · · · Delaunay 格子

1. 各単体の外心を求める

2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン

10

(16)

Delaunay 格子と Voronoi ポリゴン

単体分割によって作られた格子 · · · Delaunay 格子

1. 各単体の外心を求める

2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン

10

(17)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

11

(18)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

11

(19)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

11

(20)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

11

(21)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

11

(22)

Hinge Voronoi ポリゴンと Deficit Angle

上から見た図 開いた図 Hinge: h

Voronoi Polygon: h

Area: A

h

Deficit Angle: ε

h

h h ε h : 1 1 1 対応

11

(23)

Regge Calculus とは

(T. Regge, Il Nouvo Cim. 19, 558 (1961))

多様体の Dynamics を単体分割で表現するには

1. 各単体の形状を変化させる

2. 単体の配置を変化させる

12

(24)

Regge Calculus とは

(T. Regge, Il Nouvo Cim. 19, 558 (1961))

多様体の Dynamics を単体分割で表現するには

1. 各単体の形状を変化させる

−→ Regge Calculus 2. 単体の配置を変化させる

12

(25)

Regge Calculus とは

(T. Regge, Il Nouvo Cim. 19, 558 (1961))

多様体の Dynamics を単体分割で表現するには

1. 各単体の形状を変化させる

−→ Regge Calculus 2. 単体の配置を変化させる

−→ Dynamical Triangulation

12

(26)

Regge Calculus とは

(T. Regge, Il Nouvo Cim. 19, 558 (1961))

多様体の Dynamics を単体分割で表現するには

1. 各単体の形状を変化させる

−→ Regge Calculus 2. 単体の配置を変化させる

−→ Dynamical Triangulation

12

(27)

2. Regge Calculus における Hilbert 作用

(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))

Ricci Scalar Curvature

Proper Volume Element

Hilbert Action

13

(28)

Ricci Scalar Curvature の変形

Dimension

R = g µρ g νσ R µνρσ

Gaussian Curvature

14

(29)

Ricci Scalar Curvature の変形

R = D (D 1) K

Dimension

R = g µρ g νσ R µνρσ

Gaussian Curvature

14

(30)

Ricci Scalar Curvature の変形

R = D (D 1) K

Dimension

R = g µρ g νσ R µνρσ

Gaussian Curvature

14

(31)

Ricci Scalar Curvature の変形

R = D (D 1) K

Dimension

R = g µρ g νσ R µνρσ

Gaussian Curvature

14

(32)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

r

15

(33)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

Vi

15

(34)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

Vf

Vi

15

(35)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

Vf

Vi

15

(36)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

r

Vf

Vi

15

(37)

Gaussian Curvature

K = lim

Area 0

Angle that Vector is Rotated

Area Circumnavigated = 1 r 2

r

Vf

Vi

15

(38)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

16

(39)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

16

(40)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

16

(41)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

16

(42)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

V

f

16

(43)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

V

f

ε

h

16

(44)

Gaussian Curvature

Vi Vf

Area Circumnavigated

= A h

V

i

V

f

ε

h

Angle that Vector

is Rotated

= ε h

16

(45)

Ricci Scalar Curvature in Regge Calculus

K Regge = ε h A h

(

−→ 1 r 2

)

R Regge = D (D 1) ε h A h

17

(46)

2. Regge Calculus における Hilbert 作用

(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))

Ricci Scalar Curvature

Proper Volume Element

Hilbert Action

18

(47)

複体 · · · 有限個の単体からなる立体

= + + · · · +

19

(48)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

S

Nh1

20

(49)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

h = { H

0

, H

1

, · · · , H

D2

}

S

Nh1

20

(50)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

h = { H

0

, H

1

, · · · , H

D2

}

S

1

S

Nh1

20

(51)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

h = { H

0

, H

1

, · · · , H

D2

}

S

0

S

1

S

2

S

Nh1

20

(52)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

h = { H

0

, H

1

, · · · , H

D2

}

S

0

S

1

S

2

S

Nh1

{ H

0

, H

1

, · · · , H

D2

, S

0

, S

1

}

20

(53)

複体の構造

N

h

個の単体 {S

0

, S

1

, · · · , S

Nh1

} から成る D 次元複体

S

Nh1

S

1

S

0

h = { H

0

, H

1

, · · · , H

D2

}

S

0

S

1

S

2

S

Nh1

{ H

0

, H

1

, · · · , H

D2

, S

0

, S

1

}

S i = { H 0 , H 1 , · · · , H D 2 , S i , S i+1 mod N

h

}

20

(54)

Voronoi ポリゴンの分解

C

0

C

1

C

Nh1

: Circumcenter of S

i

C

i

S

0

S

1

S

Nh1

h

= { C

0

, C

1

, · · · , C

Nh1

}

21

(55)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

: Circumcenter of S

i

: Circumcenter of Hinge C

i

O

h

S

0

S

1

S

Nh1

h

h

= { C

0

, C

1

, · · · , C

Nh1

}

21

(56)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

h

21

(57)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

h

21

(58)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

= O

h

C

0

C

1

+

+ +

+ · · · + h

21

(59)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

= O

h

C

0

C

1

+

+ +

+ · · · + h

h

0

21

(60)

Voronoi ポリゴンの分解

O

h

C

0

C

1

C

Nh1

= O

h

C

0

C

1

+

+ +

+ · · · + h

h

0

h

i

= { O

h

, C

i

, C

i+1 mod Nh

}

A

h

i

: Area of h

i

A

h

=

N

h1 i=0

A

h

i

21

(61)

Hinge の分解

Oh

H0

H2

H1

h

h

2

H3

22

(62)

Hinge の分解

Oh

H0

H2

H1

H3

h

h

2

22

(63)

Hinge の分解

=

Oh Oh

H0

H2

H1

H3

H1

H2

H3

Oh

H0

H2

Oh

H0

H1

Oh

H0

H2

H1

H3

H3

h

h

2

22

(64)

Hinge の分解

=

Oh Oh

H0

H2

H1

H3

H1

H2

H3

Oh

H0

H2

Oh

H0

H1

Oh

H0

H2

H1

H3

H3

h h

3

h

2

22

(65)

Hinge の分解

=

Oh Oh

H0

H2

H1

H3

H1

H2

H3

Oh

H0

H2

Oh

H0

H1

Oh

H0

H2

H1

H3

H3

h h

3

h

0

h

1

h

2

22

(66)

Hinge の分解

=

Oh Oh

H0

H2

H1

H3

H1

H2

H3

Oh

H0

H2

Oh

H0

H1

Oh

H0

H2

H1

H3

H3

h h

3

h

0

h

1

h

2

h

i

= { O

h

, H

0

, · · · , H

i1

,

, H

i+1

, · · · , H

D2

}

A

hi

: Volume of h

i

A

h

=

D

2 i=0

A

hi

22

(67)

Proper Volume Element の定義

C1

C2 C0

23

(68)

Proper Volume Element の定義

C1

C2 C0

h

23

(69)

Proper Volume Element の定義

C1

C2 C0

h

23

(70)

Proper Volume Element の定義

C1

C2 C0

23

(71)

Proper Volume Element の定義

C1

C2 C0

V

Volume: ∆V

proper

23

(72)

Proper Volume Element の定義

C1

C2 C0

V

Volume: ∆V

proper

23

(73)

Proper Volume Element の定義

C1

C2 C0

=

C1 C0

C1

C2

+ +

C2 C3

V

Volume: ∆V

proper

+ · · ·

23

(74)

Proper Volume Element の定義

C1

C2 C0

=

C1 C0

C1

C2

+ +

C2 C3

V

Volume: ∆V

proper

+ · · ·

V

0

V

1

V

2

V

0

V

1

V

2

23

(75)

Proper Volume Element の定義

C1

C2 C0

=

C1 C0

C1

C2

+ +

C2 C3

V

Volume: ∆V

proper

+ · · ·

V

0

V

1

V

2

V

0

V

1

V

2

V

i

= { H

0

, H

1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

}

∆V

proper

=

N

h1 i=0

V

i

23

(76)

V i の分解

V

i

= { H

0

, H

1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

}

= h

−→D2

j=0

h

j = Dj=02{Oh,H0,···,Hj1,,Hj+1,···,HD2}

24

(77)

V i の分解

V

i

= { H

0

, H

1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

}

= h

−→D2

j=0

h

j = Dj=02{Oh,H0,···,Hj1,,Hj+1,···,HD2}

24

(78)

V i の分解

V

i

= { H

0

, H

1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

}

= h

−→D2

j=0

h

j = Dj=02{Oh,H0,···,Hj1,,Hj+1,···,HD2}

24

(79)

V i の分解

V

i

= { H

0

, H

1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

}

= h

−→D2

j=0

h

j = Dj=02{Oh,H0,···,Hj1,,Hj+1,···,HD2}

V

ij

≡ { O

h

, H

0

, · · · , H

j1

,

, H

j+1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

} Vol. (

V

ij

) = V

ij

V

i

=

D

2 j=0

V

ij

24

(80)

V i j の構造

V

ij

≡ { O

h

, H

0

, · · · , H

j1

,

, H

j+1

, · · · , H

D2

, C

i

, C

i+1 mod Nh

} 頂点 D + 1 · · · V

ij

D- 単体

Σ

D2

H1

H0 Oh

H2

C

i+1 mod Nh

C

i

Hj+1

Hj1

h

ij

−−−→ O

h

C

i

⊥ −−−→ O

h

H

j

25

(81)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

26

(82)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

26

(83)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

26

(84)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

26

(85)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

26

(86)

V i j の値

V

ij

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

∧−−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 1 D!

( −−−→ O

h

H

0

∧ · · · ∧ −−−−−→ O

h

H

j1

∧ −−−−−→ O

h

H

j+1

∧ · · · ∧ −−−−−−→ O

h

H

D2

)

= (D 2)!A

hj

× ( −−−→ O

h

C

i

∧ −−−−−−−−−−→ O

h

C

i+1 mod Nh

)

= 2A

hi

= 2

D (D 1) A

hj

A

hi

26

(87)

Proper Volume Element in Regge Calculus

∆V

proper

=

N

h

1 i=0

D 2

j =0

V i

j

=

N

h

1 i=0

D 2

j =0

2

D (D 1) A h

j

A h

i

= 2

D (D 1) A h A h

27

(88)

2. Regge Calculus における Hilbert 作用

(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))

Ricci Scalar Curvature

Proper Volume Element

Hilbert Action

28

(89)

Hilbert Action in Regge Calculus

S Regge = 1 16π

hinges

∆V proper R Regge

= 1 16π

hinges

( 2

D (D 1) A

h

A

h

) (

D (D 1) ε

h

A

h

)

= 1 8π

hinges

A

h

ε

h

29

(90)

3. Regge Calculus における Einstein 方程式

30

(91)

Regge Analogue of Metric

計量テンソル · · · (大雑把に言うと)時空の形を決めるもの

n- 単体 · · ·

n(n+1)2

本の辺の長さ全てを指定すると,形が一意的に決定する.

※一般の立体ではこうはいかない

2 次元の例:三角形 · · · 辺の長さを決めると形が決まる.

:四角形 · · · 辺の長さが等しい菱形が存在する.

基本変数 : g

µν

−→ l

p

: Delaunay 格子に含まれる辺のうちの p 番目のものの 長さ

31

(92)

Hilbert Action の変分

(L. Schl¨afli, Quart. J. Pure Appl. Math. 2, 269 (1858))

∂S

Regge

∂l

p

= 1 8π

hinges

∂l

p

(A

h

ε

h

) = 1 8π

 ∑

hinges

∂A

h

∂l

p

ε

h

+ ∑

hinges

A

h

∂ε

h

∂l

p

= 0 ( Schl¨ afli Identity)

真空の Einstein 方程式

0 = 1 8π

hinges

∂A

h

∂l

p

ε

h

※ Schl¨ afli Identity は 表面項が Gauss の定理によって落とされるのに似てい る.しかし両者の関係はあまり分かっていない.

32

参照

関連したドキュメント

From (3.2) and (3.3) we see that to get the bound for large deviations in the statement of Theorem 3.1 it suffices to obtain a large deviation bound for the continuous function ϕ k

These results will then be used to study Sobolev spaces on Lie manifolds with true boundary, which in turn, will be used to study weighted Sobolev spaces on polyhedral domains. The

Goal of this joint work: Under certain conditions, we prove ( ∗ ) directly [i.e., without applying the theory of noncritical Belyi maps] to compute the constant “C(d, ϵ)”

The Mixmaster (Bianchi IX) model of the early history of the universe is neatly explained in this picture by postulating that the reverse Wick rotation follows a hyperbolic

(Recent result: Yes, but consistent quantum gravity is delicate.) Early universe cosmology: Observations of cosmic microwave background, maybe even earlier stages with

We develop applications of the results obtained and some other techniques in variational analysis to generalized differential calculus involving normal cones to nonsmooth and

We develop applications of the results obtained and some other techniques in variational analysis to generalized differential calculus involving normal cones to nonsmooth and

In Section 3, we give applications of the RMR property and related results to first- order calculus rules for generalized normals, coderivatives, and subgradients of sets,