Sakharovの3条件
5. 電弱バリオン数生成
5.1 電弱バリオン数生成の概要
電弱理論やその拡張に基づく −→ 実験で検証可能 (逆に言えば、制限がきつい)
✬
✫
✩
✪
(1) バリオン数非保存過程 = スファレロン過程 (∆(B + L) ̸= 0) 平衡ならwash out −→ 生成直後にdecouple
(2) CP対称性の破れ
KM位相では不十分[後述]
=⇒ 標準理論の拡張 — SUSY-SM, extra Higgs, · · ·
(3) 非平衡状態
t¯EW= 10GeV < t¯(sym)sph = 103GeV ≪ H(T)−1= 1014GeV at T ≃ 100GeV
−→ 宇宙膨張は無視できる
=⇒ 電弱相転移が、相境界の形成・成長を伴う一次転移
高温(対称)相
∆B 6= 0
∆B = 0
非対称相 vC vco
v
wbroken phase
symmetric phase
z v
∆B 6= 0
∆B = 0 TR
→L
s TR
→R s
,
RR
→L s
ψR
bubble wall (= Higgs + gauge config.)とのCPを破る相互作用 — B-conserving
⇓
ψLとψRの反射率の差 + bubbl wallの運動
⇓
chiral charge (QL ̸= QR) が対称相に流れ込む flux ∝ (QL − QR) × (RsR→L − RsL→R)
⇓
スファレロン過程にバイアス µB ̸= 0
⇓
バリオン数生成 n˙B = −µB
T Γ(s)sph =⇒ 生成後直ちに非対称相で凍結
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 3/29
解くべき問題
⋆ 電弱相転移 Electroweak Phase Transition 静的
{ 相転移の次数
相転移の強さ vC/TC — sphaleron過程のdecoupling
動的 一次相転移のダイナミクス
{ bubble nucleation rate wall velocity
⋆ CP対称性の破れ — EWPTでのCP violation
←− bubble wall profile = gauge-Higgs系の古典解
⋆ 非平衡過程の取扱 — Leptogenesisのように空間的に一様でない 主にmacroscopicな拡散方程式を使って評価されている
初めの2つはモデルの選択も含むが、何れも場の理論の基本的問題と関係
標準理論でEWBGは可能か?
◃ mh > 114GeV (LEPII) =⇒ 電弱相転移が二次またはcross over [後述]
◃ CP対称性の破れはKM位相だけ
W
U mj U†
q q
mi ̸= mj
⇒ dispersionにO(αW)だけのCP violation effect
[Farrar and Shaposhnikov, Phys. Rev. D50 ('94)]
Weak int.の1-loop効果なので、かなり小さいが、さらに、
• QCD correction (short range) → decoherence [Gavela, et al., Nucl. Phys. B430 ('94)]
• bubble wallとの多重散乱 [Huet and Sather, Phys. Rev. D51 ('95)]
−→ ¯¯¯nB s
¯¯¯ < 10−26
非平衡状態の実現とCP対称性の破れの2つの条件を満たすことが出来ない。
−→ 1st-order EWPT & CP violtion の2条件を満たすようなモデルは?
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 5/29
5.2 電弱相転移
T = 100GeV程度では, t¯EW ≃ 10GeV−1≪H(T)−1 ≃ 1014GeV−1
→ 平衡系の統計力学が使える — (相転移の静的な性質)
相転移の記述
磁性体(Landauの現象論) 電弱理論
秩序変数 自発磁化 (M) ⟨Φ(x)⟩ = v
自由エネルギー F(M;T) = a(T)M2 + b(T)M4 Effective potential Veff(v;T) 計算法 例えば、スピン模型の平均場近似 有限温度の場の理論
場の理論の有効ポテンシャル (T = 0) ⇔ ⟨Φ⟩ = v という拘束条件付きの最小エネルギー
−→ Veff(v) が最小となる v が真空期待値 [Coleman, Secret Symmetry]
有効ポテンシャルが必要な場合
= 古典論の対称性が、量子補正により自発的に破れる
[Coleman and Weinberg, Phys. Rev. D7 ('73)]
ゼロ温度での計算法
Veff(v) = −Γ[φ(x) = v]/
∫
d4x とも表される
ここで、Γ[φ] = Effective Action = generating functional of 1PI vertex functions
摂動論的には、order parameter v をbackgroundとするpropagatorとvertexを使って vacuum-to-vacuum amp.を計算 systematicな計算法: [Jackiw, Phys. Rev. D9 ('74)]
Veff(v) = Vtree(v) − i
2Tr log(
iD−1(v))
+ 2-loop contrib.
= Vtree(v) − i 2
∫ d4k
(2π)4 log (
k2 − m(v)2)
+ · · · + 2-loop contrib.
有限温度では、Effective potential = 自由エネルギー密度 Veff(v;T) = − 1
V T log Z = − 1
V T log Tr [
e−H/T ]
⟨Φ⟩=v → 最小にする v ⇒ 相 at T
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 7/29
相転移の次数
T > T > 0c T > T > 0c Veff(v;T) – Veff(0;T)
v v
v0
v0
vC
v0
v(T)
T
v(T)
v0
vC
TC T
TC
2nd order PT 1st order PT
標準理論
order parameter:
⟨Φ⟩ = 1
√2
( 0 φ
)
... 1st order EWPT
⇕ vC ≡ lim
T↑TC φ(T) ̸= 0
Effective potential Veff(v;T) の計算
• Dolan and Jackiw, Phys. Rev. D9 ('74)
• Landsman and van Weert, Phys. Rep. 145 ('87)
• Kapusta, Finite-Temperature Field Theory (Cambridge Univ. P.)
Imaginary-time formalism — Feynman graphはT = 0と同じ。T = 0部分との分離が面倒。
path integral: Tr(e−H/T) = N(T)
∫
pbc
[dϕ] exp (
−
∫ 1/T 0
d4xE LE )
◃ Euclidean ActionをEuclidean time [0,1/T] で定義
◃ 場の境界条件は、分配函数の Tr のために、
{ ϕ(0,x) = ϕ(1/T,x) boson
ψ(0,x) = −ψ(1/T, x) fermion
−→ Fourier mode k0 = iωn with ωn =
{ 2nπT (boson)
(2n + 1)πT (fermion)
Real-time formalism (Thermofield dynamics) ここでは使わない。
Trを表現するために、1つの場について1対の場 (ϕ(x),ϕ(x))˜ が必要。
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 9/29
摂動計算 T = 0でのloop積分で、時間積分をMatsubara frequency ωn の和に置き換える Veff(v) −→ Veff(v;T) by
∫ d4k
(4π)4 −→ i β
∑∞ n=−∞
∫ d3k
(2π)3(· · ·)
¯¯¯¯
¯k0=iωn
例えば、1-loop effective potentail に現れる積分は i
β
∑
n
∫ d3k
(2π)3 log(k2 − m2) =
∫ d4k
(4π)4 log(k2 − m2)±2i β
∫ d3k
(2π)3 log (
1 ∓ e−β
√k2+m2
)
=
∫ d4k
(4π)4 log(k2 − m2)±iT4 π2
∫ ∞
0
dx x2 log (
1 ∓ e−
√x2+(m/T)2)
↑
T = 0-contribution — T = 0のcountertermで繰り込み可 一般的に、T = 0のcountertermで全てのループ積分を有限にできる。
標準理論 — 1-loop の摂動論 (W, Z, top quarkのloop)
Veff(φ;T) = −1
2µ2φ2 + λ
4φ4 + 2Bv02φ2 + Bφ4 [
log
(φ2 v02
)
− 3 2
]
+ V¯(φ;T)
where B = 3
64π2v04(2m4W + m4Z − 4m4t), V¯ (φ;T) = T4
2π2 [6IB(aW) + 3IB(aZ) − 6IF(at)], (aA = mA(φ)/T) IB,F(a) ≡
∫ ∞
0
dx x2 log (
1 ∓ e−
√x2+a2
) .
high-temperature expansion [m/T << 1] γE = 0.5772· · ·
IB(a) = −π4
45+π2
12a2−π
6(a2)3/2 − a4 16 log
√a2
4π − a4 16
(
γE − 3 4
)
+ O(a6) IF(a) = 7π4
360−π2
24a2 − a4 16 log
√a2
π − a4 16
(
γE − 3 4
)
+ O(a6)
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 11/29
T > mW, mZ, mtとして展開を適用すると、
Veff(φ;T) ≃ D(T2 − T02)φ2−E Tφ3 + λT 4 φ4
where
D = 1
8v02(2m2W + m2Z + 2m2t), E = 1
4πv03(2m3W + m3Z) ∼ 10−2 λT = λ − 3
16π2v04 (
2m4W log m2W
αBT2 + m4Z log m2Z
αBT2 − 4m4t log m2t αFT2
)
T02 = 1
2D(µ2 − 4Bv02), logαF(B) = 2 log (4)π − 2γE TCで、φ = 0と縮退した極小がφCに存在: φC = 2E TC
λTC Γ(br)sph < H(TC) ⇐⇒ φC
TC >∼ 1 =⇒ λ に上限 [mH = √
2λv0]
−→ mh <∼ 46GeV =⇒ 標準理論はexcluded! — mSMh > 114GeV by LEPII
Sphaleron decoupling condition
broken phase での バリオン数変化率 [Arnold and McLerran, Phys. Rev. D36 ('87)]
1 B
dB
dt ≃ −13Nf 128π2
ω−
α3W κNtrNrot e−Esph/T これがH(T) = 1.66√
g∗T2/mPより小さいとおくと, Esph ≡ 4πv
g2 Eとすれば, v
T > 0.050 E
[
43.17 + log(κNtrNrot) + log
( ω− mW
)
− 2 log
( T 100GeV
)]
E = 2.00, NtrNrot = 80.13, ω−2 = 2.3m2W ( at λ/g22 = 1) を採用し, κ = 1, T = 100GeV とすると,
v
T > 0.025 × (43.17 + 4.38 + 0.416) = 1.20
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 13/29
⋆ Monte Carlo simulations [標準理論]
effective fermion mass : mf(T) ∼ O(T) ← |ωn| = |(2n + 1)πT| ≥ πT ... bosonsだけでsimulation
格子場の理論
{ scalar fields: ϕ(x)→ 格子点 (site) gauge fields: Uµ(x)→リンク
Z =
∫
[dϕ dUµ]exp{−SE[ϕ, Uµ]}
• 3-dim. SU(2) system with a Higgs doublet and a triplet time-component of Uµ
[Laine & Rummukainen, hep-lat/9809045]
• 4-dim. SU(2) system with a Higgs doublet [Csikor, hep-lat/9910354]
EWPT is first order for mh < 66.5 ± 1.4GeV Both the simulations found end-point of EWPT at
mh =
{ 72.3 ± 0.7 GeV
72.1 ± 1.4 GeV ⇒ no PT (cross-over) in the MSM !
5.3 理論の拡張 — 電弱バリオン数生成のために 標準理論での困難 =
{ 電弱一次相転移
十分なCP対称性の破れ
◃ EWPTが一次転移となるには、
boson loop からの Veff(v;T) への寄与 ∼ −T(
m(v)2)3/2
Higgs bosonと相互作用するbosonで, m(v)2 ∼ g2v2 (for v ∼ 0)となるもの 例) 2HDMのHiggs場, SUSY-SMのsfermion
m(v)2 = m20 + g2v2 (m20 ≪ g2v02)
◃ 新たなCP対称性の破れ
⋆ scalar self-interactionのcomplex parameters λ6,7 in 2HDM, µB, A in SUSY-SM
⋆ complex Majorama mass gaugino, Higgsino soft masses in SUSY-SM
⋆ スカラー場の期待値がcomplex — 2つ以上の複素スカラー場の相対位相
EWBGに必要なのは、TCで形成されるbubble wall近傍でのCP violation
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 15/29
標準理論を拡張した模型の例
• MSSM (Minimal Supersymmetric Standard Model)
[標準理論 +{1 Higgs doublet] + それらのsuperpartners
Supersymmetric Yukawa int. (∈ superpotential) gauge anomaly cancellation
• 2HDM (two-Higgs-doublet Model)
[標準理論 + 1 Higgs doublet], 即ち、2つのHiggs doublet Φ1, Φ2 Yukawa int.にどちら or 両方のHiggs doubletを含むか,と
Higgs potential により幾つかのvariation
• NMSSM (Next-to-Minimal Supersymmetric Standard Model) MSSM + 1 Singlet Superfield
新しいタイプの強い一次相転移 [KF, Tao and Toyoda, Prog. Theor. Phys. 114 ('05)]
以下では、MSSMにおける電弱相転移を紹介
Higgs potential (tree level)
MSSM:
superpotential W = fAB(e) HdLAEB + fAB(d) HdQADB − fAB(u)HuQAUB − µHdHu V0 = m21Φ†dΦd + m22Φ†uΦu − (
m23ΦdΦu + h.c.)
+ g22 + g12 8
(
Φ†dΦd − Φ†uΦu )2
+ g22 2
¯¯¯Φ†dΦu¯¯¯2 soft-SUSY-br. terms D-term potential
2HDM: 最も一般的な、gauge-inv. renormalizable potential V0 = m21Φ†1Φ1 + m22Φ†2Φ2 + (m23Φ†1Φ2 + h.c.)
+ 1
2λ1(Φ†1Φ1)2 + 1
2λ2(Φ†2Φ2)2 + λ3(Φ†1Φ1)(Φ†2Φ2) − λ4(Φ†1Φ2)(Φ†2Φ1)
−
{1
2λ5(Φ†1Φ2)2 + [
λ6(Φ†1Φ1) + λ7(Φ†2Φ2)]
(Φ†1Φ2) + h.c.
}
2HDM→MSSM by
Φ˜1 → Φd, Φ2 → Φu, λ1 = λ2 = −λ3 = g22 + g12
4 , λ4 = −g22
4 , λ5,6,7 = 0
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 17/29
Higgs particles in the MSSM (tree level)
Φd,Φu: 4つの複素成分 (実8成分) − 3つのNG modes = 5 = 3 (neutral) + 2 (H±) Higgs potential V0の最小とHiggs mass (tree level)
V0 = min at Φd = 1
√2
(v0cosβ 0
)
, Φu = 1
√2
( 0 v0sinβ
)
但し、m21 = m23 tanβ − 1
2m2Z cos(2β), m22 = m23cotβ + 1
2m2Z cos(2β) Φd,uの位相の採り方により、m23を常にreal positiveにできる
V0のHiggs fieldsによる2階微分 → mass matrix of the Higgs particles 中性Higgs
m2A = m23
sinβ cosβ CP-odd
m2h,H = 1 2
[
m2Z + m2A ∓ √
(m2Z + m2A)2 − 4m2Zm2A cos2(2β) ]
CP-even 荷電Higgs m2H± = m2W + m2A
Higgs potentialのϕ4項の係数 ∼ g21, g12 −→ 軽いHiggs粒子 実際、 mh ≤ min{
m2Z, m2A}
, mH ≥ max{
m2Z, m2A}
−→ LEPII の結果と矛盾
第3世代のquark, squarkのloop correctionが重要
[Okada, Yamaguchi and Yanagida, Prog. Theor. Phys. 85 ('91)]
new upper bound on mh: m2h ≤ m2Z cos2(2β) + 3 2π2
m4t v02 log
(m2t˜+ m2t m2t
)
最近は、light-Higgs scenario も注目されている。
→ h-H mixingのために、ZZh-couplingが小さくなる
e
+
e
−
Z
h Z
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 19/29
MSSMのHiggs massに対する制限
1 10
0 20 40 60 80 100 120 140 1
10
mh (GeV/c2)
tanβ
Excluded by LEP
Theoretically Inaccessible mh-max
(b)
99.7%CL 95%CL
mt= 169.3,174.3,179.3, 183.0GeV
mh-max benchmark scenario
[PDG: W.-M. Yao et al., Journal of Physics G 33, 1 (2006)]
mA (GeV/c2)
80 100 120 140 160 180 200
0 20 40 60 80 100
no mixing
LEP 2
no mixing
mmaxh
mmaxh
CDF
CDF and D0
MSSM Higgs Searches Preliminary
D0
tanβ
allowed region for the lightest neutral Higgs boson allowed region for the pseudoscalar Higgs boson
⋆ MSSM の EWPT
order parameter: ⟨Φd⟩ = 1
√2
(v1 0
)
, ⟨Φu⟩ = 1
√2
( 0 v2 + iv3
)
= eiθ
√2
( 0 vu
)
Veff(v;T)を v = (v1, v2, v3)の函数として最小点を求める。(at each T) 特殊な場合には1次元問題に帰着
• symmetric under Φ1 ↔ Φ2
• m1 = m2 = m3 = 0 (tan4β = λ1/λ2) [KF, Kakuto, Takenaga, Prog. Theor. Phys. 91]
Higgs bosonの質量とEWPの関係
tree-levelの時と同様に
Higgs mass2 matrix: M2ij ≃
⟨∂2Veff(T = 0)
∂ϕi∂ϕj
⟩
= Γ(2)1PI(p = 0)
−→
{ 質量固有値 mHi
対角化する行列 → gZZHi
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 21/29
1-loop levelの計算 [Carena, Ellis, Pilaftsis and Wagner, Nucl. Phys. B586 ('00)]
Veff(v;T = 0) = V0(v) + 1 64π2
∑
a
ca( ¯m2a)2 (
log m¯ 2a
M2 − 3 2
)
¯
m2a: field (Φd,Φu)-dependent mass2, ca: 統計因子 ループを回る粒子: a = t, t, b,˜ ˜b, W, Z, · · ·
field parametrization — VEV + fluctuation Φd =
( √1
2(vd + hd + iad) ϕ−d
)
, Φu = eiθ
( ϕ+u
√1
2(vu + hu + iau) )
VEV = Veff のglobal min. −→ vd = v0 cosβ, vu = v0sinβ, θ をinputにする。−→ soft mass 0 = 1
vd
⟨∂Veff
∂hd
⟩
= m21 − Re(m23eiθ)tanβ + 1
2m2Zcos(2β) + · · · , 0 = 1
vu
⟨∂Veff
∂hu
⟩
= m22 − Re(m23eiθ)cotβ − 1
2m2Zcos(2β) + · · ·, 0 = 1
vu
⟨∂Veff
∂ad
⟩
= 1 vd
⟨∂Veff
∂au
⟩
= Im(m23eiθ) + · · ·.
ここで ⟨· · ·⟩は、vacuumで評価した値 — fluctuationで微分してそれを0と置く
neutral Higgs boson と charged Higgs boson�mass: [NG modeを除去した後]
M2 =
⟨∂2Veff
∂h2
d
⟩ ⟨
∂2Veff
∂hd∂hu
⟩ 1
cosβ
⟨ ∂2Veff
∂hd∂au
⟩
⟨ ∂2Veff
∂hd∂hu
⟩ ⟨
∂2Veff
∂h2u
⟩ 1
sinβ
⟨ ∂2Veff
∂hu∂ad
⟩
1 cosβ
⟨ ∂2Veff
∂hd∂au
⟩ 1 sinβ
⟨ ∂2Veff
∂hu∂ad
⟩ 1
sinβcosβ
⟨ ∂2Veff
∂ad∂au
⟩
m2H± = 1
sinβ cosβ
⟨ ∂2Veff
∂ϕ+d∂ϕ−u
⟩
= 1
sinβ cosβRe(m23eiθ) + m2W + · · ·
−→ input m2H± → Re(m23eiθ) mass eigenstates Hi
hd hu a
= O
H1 H2 H3
, OTM2O = diag(m2H
1, m2H
2, m2H
3)
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 23/29
gauge and Yukawa interactions Lgauge ∼ g2mW gV V Hi
(
Wµ+W−µ + ZµZµ 2 cos2θW
)
Hi + g2
2 cos θW gZHiHjZµ (
Hi↔∂µHj )
LY ∼ −g2mb 2mW
¯b(gbbhS
i + iγ5gbbhP
i)bHi
corrections to the couplings [標準理論: gV V H = 1, gZHH = 0, gbbH = 1]
gV V Hi = O1i cosβ + O2isinβ gZHiHj = 1
2 [(O3iO1j − O3jO1i) sinβ + (O3iO2j − O3jO2i) cos β]
gbbHS
i = O1i 1
cosβ, gbbHP
i = −O3i tanβ, gbbH2
i = (
gbbHS
i
)2 + (
gbbHP
i
)2
v0 = 246GeV,tanβ, mH± と、loopを通して効くパラメータ (µ, Aq,scalar soft mass,· · ·) をinputとして、これらの量と m21, m22, m23 を計算
−→ Higgs mass vs EWPT
EWPTが一次転移になるparameter
light stop scenario [de Carlos, Espinosa, Nucl. Phys. B503 ('97)]
stop mass2 matrix:
M2t˜ =
m2t˜
L+ (g12
24 − g822)
(vu2 − vd2) + y2t2vu2 √yt
2 (µvd + A(v2 − iv3))
∗ m2t˜
R − g612(vu2 − vd2)+ y22tvu2
m2t˜
L = 0 or m2t˜
R = 0 =⇒ smaller eigenvalue: m2t˜
1 ∼ O(v2) ... 高温展開
V¯t˜(v;T) ⇒ − T
6π(m2t˜
1)3/2 ∼ Tv3
−→ stronger 1st order PT
effective for larger yt — smaller tanβ mt = √1
2 ytv0sinβ
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 25/29
��� (CP-conserving case) [KF, Prog. Theor. Phys. 101 ('99)]
tan β = 6, mh = 82.3GeV, mA = 118GeV, mt˜
1 = 168GeV TC = 93.4GeV, vC = 129GeV
0 50 100 150
0 5 10 15 20 25
0.0x100 2.0x105 4.0x105 6.0x105 8.0x105 1.0x106 1.2x106 1.4x106 1.6x106 1.8x106 2.0x106
v1 v2
0.6 0.8 1 1.2 1.4
6 6.1 6.2 6.3 6.4 6.5
0 10 20 30 40 50 60 70 80
60 80 100 120 140 160 180 200
0 10 20 30 40 50 60 70 80
(GeV) (GeV)
GeV
m~t
R m~t
R
v / TC C tanβ
mh mA m~t1
Veff(v1, v2, v3 = 0;TC) mtR-dependence (tanβ = 6)
⋆ Lattice MC studies
• 3d reduced model [Laine et al. hep-lat/9809045]
strong 1st order for mt˜1 <∼ mt and mh ≤ 110GeV
• 4d model [Csikor, et al. hep-lat/0001087]
with SU(3), SU(2) gauge bosons, 2 Higgs doublets, stops, sbottoms At,b = 0, tan β ≃ 6 −→ errorの範囲内で摂動論と一致
v /T =1C C mh
mA = 500 GeV
vC/TC > 1
below the steeper lines
⇓
max. mh = 103 ± 4 GeV for mt˜L ≃ 560 GeV
—ůOHFŔRZHDN%DU\RJHQHVLϳ — 27/29
light-stop scenarioでもEWPTが強い一次相転移は軽いHiggsの領域だけ
— CP-conserving MSSM
vc/T c > 1
LHCでHiggs bosonが見つかったとき、
{ mh > 110GeV ⇒ EWBG in the MSSM
X
mh > 135GeV ⇒ MSSM