24 June 2017
Expanding Polyhedral Universe in Regge Calculus
津田 廉
茨城大学 素粒子論研究室
共同研究者:藤原高徳(茨城大学)
核物理
×
物性セミナー@
千葉工業大学 新習志野キャンパスContents
1. 単体分割
2. Regge Calculus における Hilbert 作用
3. Regge Calculus における Einstein 方程式 4. Expanding Polyhedral Universe
5. Summary
2
1. 単体分割
3
2 次元曲面の三角分割
滑らかな曲面
4
2 次元曲面の三角分割
三角形で近似
5
曲率 · · · 三角形の頂点に定義
平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)
6
曲率 · · · 三角形の頂点に定義
平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)
6
曲率 · · · 三角形の頂点に定義
過不足なし
平面(曲率ゼロ) 球面(正曲率) 鞍状面(負曲率)
不足 過剰
6
n- 単体 · · · n + 1 個の頂点を持つ n 次元凸面体
0- 単体
1- 単体
2- 単体
3- 単体
.. .
点
.. .
線分 三角形 四面体
.. .
7
単体の構造
=
2- 単体
× 3 × 3
=
3- 単体
× 4 × 6 × 4
8
単体の構造
∋
2- 単体 1- 単体
× 3 0- 単体
× 3
=
,
3- 単体
× 4 × 6 × 4
8
単体の構造
∋
2- 単体 1- 単体
× 3 0- 単体
× 3
∋
, ,
,
3- 単体 2- 単体
× 4 1- 単体
× 6 0- 単体
× 4
8
単体の構造
Volume : A
bA
h名称 : Boundary Hinge: h
n- 単体 ∋ { (n − 1) - 単体 , (n − 2) - 単体 , · · · , 0- 単体 }
×
n+1C
1×
n+1C
2×
n+1C
n9
Delaunay 格子と Voronoi ポリゴン
単体分割によって作られた格子 · · · Delaunay 格子
1. 各単体の外心を求める
2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン
10
Delaunay 格子と Voronoi ポリゴン
単体分割によって作られた格子 · · · Delaunay 格子
1. 各単体の外心を求める
2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン
10
Delaunay 格子と Voronoi ポリゴン
単体分割によって作られた格子 · · · Delaunay 格子
1. 各単体の外心を求める
2. 隣り合う外心を結ぶ · · · Voronoi ポリゴン
10
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
h11
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
h11
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
h11
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
h11
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
h11
Hinge と Voronoi ポリゴンと Deficit Angle
上から見た図 開いた図 Hinge: h
Voronoi Polygon: h
∗Area: A
∗hDeficit Angle: ε
hh と h ∗ と ε h : 1 対 1 対 1 対応
11
Regge Calculus とは (T. Regge, Il Nouvo Cim. 19, 558 (1961))
多様体の Dynamics を単体分割で表現するには
1. 各単体の形状を変化させる
2. 単体の配置を変化させる
12
Regge Calculus とは (T. Regge, Il Nouvo Cim. 19, 558 (1961))
多様体の Dynamics を単体分割で表現するには
1. 各単体の形状を変化させる
−→ Regge Calculus 2. 単体の配置を変化させる
12
Regge Calculus とは (T. Regge, Il Nouvo Cim. 19, 558 (1961))
多様体の Dynamics を単体分割で表現するには
1. 各単体の形状を変化させる
−→ Regge Calculus 2. 単体の配置を変化させる
−→ Dynamical Triangulation
12
Regge Calculus とは (T. Regge, Il Nouvo Cim. 19, 558 (1961))
多様体の Dynamics を単体分割で表現するには
1. 各単体の形状を変化させる
−→ Regge Calculus 2. 単体の配置を変化させる
−→ Dynamical Triangulation
12
2. Regge Calculus における Hilbert 作用
(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))
• Ricci Scalar Curvature
• Proper Volume Element
• Hilbert Action
13
Ricci Scalar Curvature の変形
Dimension
R = g µρ g νσ R µνρσ
Gaussian Curvature
14
Ricci Scalar Curvature の変形
R = D (D − 1) K
Dimension
R = g µρ g νσ R µνρσ
Gaussian Curvature
14
Ricci Scalar Curvature の変形
R = D (D − 1) K
Dimension
R = g µρ g νσ R µνρσ
Gaussian Curvature
14
Ricci Scalar Curvature の変形
R = D (D − 1) K
Dimension
R = g µρ g νσ R µνρσ
Gaussian Curvature
14
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
r
15
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
Vi
15
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
Vf
Vi
15
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
Vf
Vi
15
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
r
Vf
Vi
15
Gaussian Curvature
K = lim
Area → 0
Angle that Vector is Rotated
Area Circumnavigated = 1 r 2
r
Vf
Vi
15
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
i16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
i16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
i16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
iV
f16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
iV
fε
h16
Gaussian Curvature
Vi Vf
Area Circumnavigated
= A ∗ h
V
iV
fε
hAngle that Vector
is Rotated
= ε h
16
Ricci Scalar Curvature in Regge Calculus
K Regge = ε h A ∗ h
(
−→ 1 r 2
)
R Regge = D (D − 1) ε h A ∗ h
17
2. Regge Calculus における Hilbert 作用
(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))
• Ricci Scalar Curvature
• Proper Volume Element
• Hilbert Action
18
複体 · · · 有限個の単体からなる立体
= + + · · · +
19
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0S
Nh−120
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0h = { H
0, H
1, · · · , H
D−2}
S
Nh−120
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0h = { H
0, H
1, · · · , H
D−2}
S
1S
Nh−120
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0h = { H
0, H
1, · · · , H
D−2}
S
0S
1S
2S
Nh−120
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0h = { H
0, H
1, · · · , H
D−2}
S
0S
1S
2S
Nh−1{ H
0, H
1, · · · , H
D−2, S
0, S
1}
20
複体の構造
N
h個の単体 {S
0, S
1, · · · , S
Nh−1} から成る D 次元複体
S
Nh−1S
1S
0h = { H
0, H
1, · · · , H
D−2}
S
0S
1S
2S
Nh−1{ H
0, H
1, · · · , H
D−2, S
0, S
1}
S i = { H 0 , H 1 , · · · , H D − 2 , S i , S i+1 mod N
h}
20
Voronoi ポリゴンの分解
C
0C
1C
Nh−1: Circumcenter of S
iC
iS
0S
1S
Nh−1h
∗= { C
0, C
1, · · · , C
Nh−1}
21
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1: Circumcenter of S
i: Circumcenter of Hinge C
iO
hS
0S
1S
Nh−1h
h
∗= { C
0, C
1, · · · , C
Nh−1}
21
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1h
∗21
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1h
∗21
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1= O
hC
0C
1+
+ +
+ · · · + h
∗21
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1= O
hC
0C
1+
+ +
+ · · · + h
∗h
∗021
Voronoi ポリゴンの分解
O
hC
0C
1C
Nh−1= O
hC
0C
1+
+ +
+ · · · + h
∗h
∗0h
∗i= { O
h, C
i, C
i+1 mod Nh}
A
∗hi
: Area of h
∗iA
∗h=
N
∑
h−1 i=0A
∗hi
21
Hinge の分解
Oh
H0
H2
H1
h
h
2H3
22
Hinge の分解
Oh
H0
H2
H1
H3
h
h
222
Hinge の分解
=
Oh Oh
H0
H2
H1
H3
H1
H2
H3
Oh
H0
H2
Oh
H0
H1
Oh
H0
H2
H1
H3
H3
h
h
222
Hinge の分解
=
Oh Oh
H0
H2
H1
H3
H1
H2
H3
Oh
H0
H2
Oh
H0
H1
Oh
H0
H2
H1
H3
H3
h h
3h
222
Hinge の分解
=
Oh Oh
H0
H2
H1
H3
H1
H2
H3
Oh
H0
H2
Oh
H0
H1
Oh
H0
H2
H1
H3
H3
h h
3h
0h
1h
222
Hinge の分解
=
Oh Oh
H0
H2
H1
H3
H1
H2
H3
Oh
H0
H2
Oh
H0
H1
Oh
H0
H2
H1
H3
H3
h h
3h
0h
1h
2h
i= { O
h, H
0, · · · , H
i−1,
∧, H
i+1, · · · , H
D−2}
A
hi: Volume of h
iA
h=
D
∑
−2 i=0A
hi22
Proper Volume Element の定義
C1
C2 C0
23
Proper Volume Element の定義
C1
C2 C0
h
23
Proper Volume Element の定義
C1
C2 C0
h
23
Proper Volume Element の定義
C1
C2 C0
23
Proper Volume Element の定義
C1
C2 C0
V
Volume: ∆V
proper23
Proper Volume Element の定義
C1
C2 C0
V
Volume: ∆V
proper23
Proper Volume Element の定義
C1
C2 C0
=
C1 C0
C1
C2
+ +
C2 C3
V
Volume: ∆V
proper+ · · ·
23
Proper Volume Element の定義
C1
C2 C0
=
C1 C0
C1
C2
+ +
C2 C3
V
Volume: ∆V
proper+ · · ·
V
0V
1V
2V
0V
1V
223
Proper Volume Element の定義
C1
C2 C0
=
C1 C0
C1
C2
+ +
C2 C3
V
Volume: ∆V
proper+ · · ·
V
0V
1V
2V
0V
1V
2V
i= { H
0, H
1, · · · , H
D−2, C
i, C
i+1 mod Nh}
∆V
proper=
N
∑
h−1 i=0V
i23
V i の分解
V
i= { H
0, H
1, · · · , H
D−2, C
i, C
i+1 mod Nh}
= h
−→ ∑D−2j=0
h
j = ∑Dj=0−2{Oh,H0,···,Hj−1,∧,Hj+1,···,HD−2}24
V i の分解
V
i= { H
0, H
1, · · · , H
D−2, C
i, C
i+1 mod Nh}
= h
−→ ∑D−2j=0
h
j = ∑Dj=0−2{Oh,H0,···,Hj−1,∧,Hj+1,···,HD−2}24
V i の分解
V
i= { H
0, H
1, · · · , H
D−2, C
i, C
i+1 mod Nh}
= h
−→ ∑D−2j=0
h
j = ∑Dj=0−2{Oh,H0,···,Hj−1,∧,Hj+1,···,HD−2}24
V i の分解
V
i= { H
0, H
1, · · · , H
D−2, C
i, C
i+1 mod Nh}
= h
−→ ∑D−2j=0
h
j = ∑Dj=0−2{Oh,H0,···,Hj−1,∧,Hj+1,···,HD−2}V
ij≡ { O
h, H
0, · · · , H
j−1,
∧, H
j+1, · · · , H
D−2, C
i, C
i+1 mod Nh} Vol. (
V
ij) = V
ijV
i=
D
∑
−2 j=0V
ij24
V i j の構造
V
ij≡ { O
h, H
0, · · · , H
j−1,
∧, H
j+1, · · · , H
D−2, C
i, C
i+1 mod Nh} 頂点 D + 1 個 · · · V
ijは D- 単体
Σ
D−2H1
H0 Oh
H2
C
i+1 mod NhC
iHj+1
Hj−1
h
ij−−−→ O
hC
i⊥ −−−→ O
hH
j25
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi26
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi26
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi26
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi26
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi26
V i j の値
V
ij= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2∧−−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 1 D!
( −−−→ O
hH
0∧ · · · ∧ −−−−−→ O
hH
j−1∧ −−−−−→ O
hH
j+1∧ · · · ∧ −−−−−−→ O
hH
D−2)
= (D − 2)!A
hj× ( −−−→ O
hC
i∧ −−−−−−−−−−→ O
hC
i+1 mod Nh)
= 2A
∗hi= 2
D (D − 1) A
hjA
∗hi26
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
jA ∗ h
i
= 2
D (D − 1) A h A ∗ h
27
2. Regge Calculus における Hilbert 作用
(W. A. Miller, Class. Quantum Crav. 14, 199 (1997))
• Ricci Scalar Curvature
• Proper Volume Element
• Hilbert Action
28
Hilbert Action in Regge Calculus
S Regge = 1 16π
∑
hinges
∆V proper R Regge
= 1 16π
∑
hinges
( 2
D (D − 1) A
hA
∗h) (
D (D − 1) ε
hA
∗h)
= 1 8π
∑
hinges
A
hε
h29
3. Regge Calculus における Einstein 方程式
30
Regge Analogue of Metric
計量テンソル · · · (大雑把に言うと)時空の形を決めるもの
n- 単体 · · ·
n(n+1)2本の辺の長さ全てを指定すると,形が一意的に決定する.
※一般の立体ではこうはいかない
2 次元の例:三角形 · · · 辺の長さを決めると形が決まる.
:四角形 · · · 辺の長さが等しい菱形が存在する.
基本変数 : g
µν−→ l
p: Delaunay 格子に含まれる辺のうちの p 番目のものの 長さ
31
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