arxiv: v1(astro-ph.co)

最も一般的なスカラーテンソル理論における

宇宙論的密度揺らぎのバイスペクトル

広島大学  修士2年   宅嶌 祐一郎

       共同研究者 照喜名 歩 

       山本 一博 

１、導入

２、基礎方程式とモードカップリング

３、密度揺らぎとバイスペクトルの表式

４、具体的なモデルを用いた考察

５、まとめと今後の課題

arXiv:1311.0281v1(astro-ph.CO)

イントロダクション

一般相対性理論（宇宙項模型）

修正重力理論

•   重力理論

ラグランジアンに の関数  を導入

•    重力理論

高階微分を含むスカラー場  の導入

•

グラビトンに質量を与える     

宇宙項 を導入

宇宙の加速膨張を説明する理論

Galileon

Massive Gravity

etc...

f(R)

f(R)

R

φ(t)

R

µν

−

1

2

Rg

µν

+ Λg

µν

=

8πG

c

4

T

µν

Λ

重力理論の特徴を調べて、検証に役立てる

最も一般的なスカラーテンソル理論

場の方程式が２階の微分方程式になるという条件の下で、

最も一般的な理論

Action

S =

�

d

4

x

√

_{−g (L}

_GG

+

_L

_m

)

L

GG

= K(φ, X)

− G

3

(φ, X)

�φ + G

4

(φ, X)R + G

4X

�

(

�φ)

2

_{− (∇}

_µ

_∇

_ν

φ)

2

�

+G

₅

(φ, X)G

_µν

_∇

µ

_∇

ν

φ

₋

1

6

G

5X

�

(

�φ)

3

_{− 3�φ(∇}

_µ

_∇

_ν

φ)

2

+ 2(

_∇

_µ

_∇

_ν

φ)

3

�

本研究では、宇宙論バックグラウンドの下で、      機構が働く 可能性がある模型を考える（          ）

Vainshtein

R.Kimura,T.Kobayashi,K.Yamamoto 2012 φ : scalar field, X _{≡ −}1 2g µν ∇µφ∇νφ, �φ ≡ ∇µ∇µφ, GiX ≡ ∂Gi ∂X

宇宙論的密度揺らぎのバイスペクトル

宇宙論的密度揺らぎ

δ(t, x)

_≡

ρ(t, x)

− ¯

ρ(t)

¯

ρ(t)

ρ(t, x) :

物質の質量密度 の空間平均

ρ(t, x)

バイスペクトル

�δ(t, k

1

)δ(t, k

2

)δ(t, k

3

)

� ≡ (2π)

3

δ

(3)

(k

₁

+ k

₂

+ k

₃

)B(t, k

₁

, k

₂

, k

₃

)

B(t, k

₁

, k

₂

, k

₃

)

δ(t, x)

実空間 フーリエ変換 δ(t, x) = 波数空間 1 (2π)3 � d3kδ(t, k)eik·x

δ(t, x)

δ(t, k)

¯

ρ(t) :

初期 ガウス分布に従う 現在 重力や流体の非線形性

δ :

δ :

ガウス分布からずれる

B = 0

B

�= 0

基礎方程式(5つ)

重力場の方程式１ 計量 スカラー場の摂動 −k2 (_FTΨ(t, k) _{− GT} Φ(t, k) _{− A1}Q(t, k)) = B1 2a2_H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Q(t, k1)Q(t, k2) + B3 a2_H2 1 (2π)3 � dk₁dk₂δ(3)(k₁ + k_{2 − k)} �k₁2k₂2 _{− (k1 · k2})2� Q(t, k₁)Φ(t, k₂) ds2 = _{−(1 + 2Φ)dt}2 + a2(1 _{− 2Ψ)dx}2 φ _{→ φ(t) + δφ(t, x), Q ≡ H} δφ ˙ φ FT = 2 � G4 − X � ¨ φG5X + G5φ �� , GT = 2 � G4 − 2XG4X − X � H ˙φG5X − G5φ �� , A₀ = . . . , .. .

H, φ, X, K, G

i もしくはそれらの微分の関数 方程式の係数

Real space

� ∇2Q�2 _{− (∂i}∂_jQ)2

重力場の方程式２ スカラー場の方程式 × � −k12k22k32 + 3k12(k2 · k3)2 − 2(k1 · k2)(k2 · k3)(k3 · k1) � Q(t, k₁)Q(t, k₂)Q(t, k₃) − B3 a2_H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Q(t, k1)Ψ(t, k2) − C1 3a4_H4 1 (2π)6 � dk1dk2dk3δ(3)(k1 + k2 + k3 − k) −k2(_GT Ψ(t, k) + A₂Q(t, k)) ₋ a 2 2 ρmδ(t, k) = ₋ B2 2a2H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Q(t, k1)Q(t, k2) −k2(A0Q(t, k) − A1Ψ(t, k) − A2Φ(t, k)) = ₋ B0 a2_H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Q(t, k1)Q(t, k2) + B1 a2_H2 1 (2π)3 � dk₁dk₂δ(3)(k₁ + k_{2 − k)} �k₁2k₂2 _{− (k1 · k2})2� Ψ(t, k₁)Q(t, k₂) + B2 a2H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Φ(t, k1)Q(t, k2) + B3 a2_H2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � k₁2k₂2 _{− (k}1 · k2)2 � Φ(t, k1)Ψ(t, k2) +_O(QQQ) ３次のモードカップリングは バイスペクトルのLeading order には効かない

連続の方程式 オイラー方程式

モードカップリングを与える関数は２種類

β(k

₁

, k

₂

) = α(k

₁

, k

₂

)

_{− γ(k}

₁

, k

₂

)

２次のオーダーで、方程式に現れるモードカップリングを与える関数は   の３種類α, β, γ 連続の方程式 オイラー方程式 重力場の方程式 スカラー場の方程式 α(k₁, k₂) = 1 + k1 · k2(k 2 1 + k22) 2k₁2k₂2 , β(k1, k2) = (k1 · k2)|k1 + k2|2 2k₁2k₂2 , γ(k1, k2) = 1 − (k1 · k2)2 k₁2k₂2 1 H ∂δ(t, k) ∂t + θ(t, k) = − 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � 1 + k1 · k2 k₂2 � δ(t, k1)θ(t, k2), 1 H ∂θ(t, k) ∂t + � 2 + H˙ H2 � θ(t, k) ₋ k 2 a2_H2 Φ(t, k) = ₋1 2 1 (2π)3 � dk1dk2δ(3)(k1 + k2 − k) � (k_{1 · k2})_|k1 + k_2|2 k₁2k₂2 � θ(t, k1)θ(t, k2).

θ

_{≡ ∇u/aH}

u : velocity field

宇宙論的密度揺らぎ 1次の場合

重力場の方程式(2つ)、

スカラー場の方程式

密度揺らぎの微分方程式

Y (t, k) = � n=1 Yn(t, k)

1次

Φ₁(t, k), Ψ₁(t, k), Q₁(t, k) が δ₁(t, k) で表される

連続の方程式

θ1(t, k) が δ1(t, k) で表される

δ

₁

(t, k) = D

₊

(t)δ

_L

(k)

δL(k) : 成長モード D+(t) : 線形密度揺らぎ

オイラー方程式

＋

L(t) = (A0FT − A 2 1)ρm 2(A0GT + 2A1A2GT + A2₂FT )

∂

2

δ

1

(t, k)

∂t

2

+ 2H

∂δ

1

(t, k)

∂t

+ L(t)δ

1

(t, k) = 0

初期はガウシアン Y = δ, θ, Φ, Ψ, Q

宇宙論的密度揺らぎ 2次の場合

２次の場合、１次の解を用いて解くと、（初期はガウシアンを仮定）

δ

_2K

(t, k) =

�

dk

₁

dk

₂

δ

(3)

(k

₁

+ k

₂

_{− k)F}

₂

(t, k

₁

, k

₂

)δ

_L

(k

₁

)δ

_L

(k

₂

)

F

₂

(t, k

₁

, k

₂

) = κ(t)α(k

₁

, k

₂

)

₋

2

7

λ(t)γ(k

1

, k

2

)

宇宙論的密度揺らぎの２次のカーネル 減衰モード           などの複雑な関数 D₋(t) : 成長モード D₊(t) : f = d ln D+ d ln a : growth rate

２次の解を特徴づけるパラメータは２つか？

κ(t) = 1 D₊2 (t) � t 0 D₋(t)D+(t�) − D+(t)D₋(t�) W [D+(t�), D₋(t�)] D₊2 (t�)�2f2H2 _{− L(t}�)� λ(t) = 7 2D₊2 (t) � t 0 D₋(t)D₊(t�) _{− D+}(t)D₋(t�) W [D+(t�), D₋(t�)] D₊2 (t�)�f2H2 _{− Nγ}(t�)�

δ

₂

(t, k) = D

₊2

(t)δ

_2K

(t, k)

N_γ(t) : H, φ, X, K, G_i

パラメータに関する微分方程式

κの微分方程式（    に比例する項の方程式）

 モードカップリングを与える関数が2つ

＋それらが互いに1次独立

2次の密度揺らぎの微分方程式

λの微分方程式（    に比例する項の方程式）

¨

κ + (4f + 2)H ˙κ + (2f

2

H

2

_{− L)κ = 2f}

2

H

2

_{− L}

¨

λ + (4f + 2)H ˙λ + (2f

2

H

2

_{− L)λ =}

7

2

(f

2

_H

2

− N

γ

)

となるので、パラメータは1つになる

κ = 1

α(k

₁

, k

₂

)

γ(k

₁

, k

₂

)

Bispectrumの表式

バイスペクトルをパワースペクトルの積の和で規格化 縦軸に          、 横軸に  と の間の角度 のグラフを用いて、 パラメータ依存性を調べる Reduced Bispectrum k₁ k₂

k

k

₃

Reduced Bispectrum

�δL(k1)δL(k2)� ≡ (2π)3δ(3)(k1 + k2)P11(k1) F2(t, k1, k2) = α(k1, k2) − 2 7λ(t)γ(k1, k2)

B(t, k

₁

, k

₂

, k

₃

) = 2D

4

(t)(F

₂

(t, k

₁

, k

₂

)P

₁₁

(k

₁

)P

₁₁

(k

₂

) + 2cyc.)

θ

₁₂

k

2 θ₁₂

δ

(3)

(k

₁

+ k

₂

+ k

₃

)

Q

₁₂₃

(t, k

₁

, k

₂

, θ

₁₂

)

_≡

B(t, k

1

, k

2

, k

3

)

P (t, k

₁

)P (t, k

₂

) + P (t, k

₂

)P (t, k

₃

) + P (t, k

₃

)P (t, k

₁

)

Reduced Bispectrum

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.5 1.0 1.5 2.0 Θ₁₂�Π Q123 k1 � k2� 0.1hMpc�1 0.0 0.2 0.4 0.6 0.8 1.0 1.0 1.5 2.0 2.5 Θ₁₂�Π Q123 k1 � 5 � k2� 0.5hMpc�1

のとき、全体的に下がり、

のとき、全体的に上がる

λ > 1

λ < 1

バイスペクトルは のみで決まる

_λ

: λ = 1 : λ = 0.8 : λ = 1.2

k

₁

= k

₂

= 0.1hM pc

−1

k

₁

= 5

_{× k}

₂

= 0.5hM pc

−1 F2(t, k1, k2) = α(k1, k2) − 2 7λ(t)γ(k1, k2) α(k1, k2) > 0, γ(k1, k2) > 0, P11(k) > 0

Reduced Bispectrum

     は    と  の関数

θ

₁₂

が同じならおおよそ同じ形（特徴）になる

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 Θ₁₂�Π Q123 k1 � k2� 0.01hMpc�1 0.0 0.2 0.4 0.6 0.8 1.0 0.8 1.0 1.2 1.4 Θ₁₂�Π Q123 k1 � 5 � k2� 0.05hMpc�1 : λ = 1 : λ = 0.8 : λ = 1.2

_k

₁

_{= 5}

_{× k}

₂

_{= 0.05hM pc}

−1

k

₁

= k

₂

= 0.01hM pc

−1

k

₂

/k

₁

k

₂

/k

₁

α(k

_i

, k

_j

), γ(k

_i

, k

_j

)

Kinetic Gravity Braiding Model

ラグランジアン

アトラクター解

3 ˙

φHG

_3X

= 1

0.0 0.2 0.4 0.6 0.8 1.0 0.98 0.99 1.00 1.01 1.02 1.03 a Λ� a� : n = 1 : n = 2 : n = 5 : ΛCDM

2.5％程度ずれる

λ(a = 1) 1.019 0.994 1.003 1.011 （       ）_{R.Kimura,K.Yamamoto 2011}

•ゴーストなし

•加速膨張解

r

_c

: crossover scale

M

_pl

=

√

8πG

K = _{−X, G}3 = Mpl � r_c2 M_pl2 X �n , G4 = M_pl2 2 , G5 = 0

Reduced Bispectrumの比較

0.0 0.2 0.4 0.6 0.8 1.0 �0.015 �0.010 �0.005 0.000 Θ₁₂�Π Q 123 �Q 123 � � 1 k1 � k2� 0.1hMpc�1 0.0 0.2 0.4 0.6 0.8 1.0 �0.015 �0.010 �0.005 0.000 Θ₁₂�Π Q 123 �Q 123 � � 1 k1� 5�k2 � 0.5hMpc�1

2％程度のずれる

: n = 1 : n = 2 : n = 5

縦軸

Q

₁₂₃

Q

_123Λ

− 1 :

KGBモデルと宇宙項模型との差

まとめ

•2次ではモードカップリングを与える関数は2種類

パラメータは1つ

•KGB modelでは  は２.5％ほどのずれた。

今後の課題

•なぜ になるのか?

•3次のオーダーの研究

•他の具体的なモデルでのパラメータ計算

•最も一般的なスカラーテンソル理論における

宇宙論的密度揺らぎのバイスペクトルを計算した。

κ(t) = 1

λ(t)

λ(t)

α(k

₁

, k

₂

), γ(k

₁

, k

₂

)