高エネルギー重イオン衝突反応における流体揺らぎの影響
上智大学
坂井あづみ
Heavy Ion Pub
2021年6月4日

目次
1. はじめに2. 高エネルギー重イオン衝突反応における測定量
3. 統合的動的模型
4. 結果
5. まとめQCD物質の物性研究

温度

ハドロン相
カラーの閉じ込め
重初イ 高期オ エ宇ン ネ宙
衝 ル
突 ギ
実 ー験

QGP相
臨界点

ハドロン相
バリオン化学ポテンシャル 𝝁𝑩

予想されるQCD相図

クォークグルーオンプラズマ（QGP）相カラーの閉じ込めから解放➢ 初期宇宙➢ 中性子星内部➢ 高エネルギー重イオン衝突実験
QGPの物性研究高エネルギー重イオン衝突反応の衝突過程
衝突前

衝突直後

ハドロン化

測定

http://wl33.web.rice.edu/research.html

QGP

光速近くまで加速
ローレンツ収縮すり抜け
カラーフラックス
チューブ生成QGPは数fm/𝑐でハドロン化

測定量：
ハドロンの種類や運動量高エネルギー重イオン衝突反応の衝突過程
衝突前

衝突直後ハドロン化

測定

http://wl33.web.rice.edu/research.html

QGP

光速近くまで加速ローレンツ収縮

すり抜け
カラーフラックス
チューブ生成QGPは数fm/𝑐でハドロン化

測定量：
ハドロンの種類や運動量

QGPの物性？

高エネルギー重イオン衝突反応の衝突過程
衝突前

衝突直後

ハドロン化

測定

http://wl33.web.rice.edu/research.html

QGP

光速近くまで加速ローレンツ収縮

統合的動的模型

すり抜け
カラーフラックス
チューブ生成初期条件

QGPは数fm/𝑐でハドロン化

流体模型

測定量：
ハドロンの種類や運動量

ハドロンカスケード

測定量高エネルギー重イオン衝突反応における揺らぎ
衝突前

衝突直後

ハドロン化

測定

http://wl33.web.rice.edu/research.html

QGP

統合的動的模型

初期条件

流体模型

初期揺らぎ

流体揺らぎ

ハドロンカスケード

測定量

初期揺らぎ
エントロピー密度

エントロピー密度

初期揺らぎ

衝突事象毎に核子分布を考慮
衝突軸

中心衝突
滑らかな分布

衝突軸

中心衝突ホットスポットが存在

初期揺らぎ⇒フローの理解
B. Alver, G. Roland, Phys. Rev. C 81, 054905 (2010)

衝突軸方向の初期揺らぎ
紐破砕をモデル化
粒子生成の分布を考慮
𝜂𝑠
𝜂𝑠

衝突軸方向の初期揺らぎ
紐破砕をモデル化
粒子生成の分布を考慮
𝜂𝑠
𝜂𝑠

カラーフラックスチューブ
衝突軸方向に滑らか
T.Hirano et al., Phys.Lett.B 636, 299 (2006)

粒子分布による揺らぎ
衝突軸方向に凹凸
T. Sjöstrand et al., Comput. Phys. Commun. 191, 159 (2015)
M. Okai et al., Phys. Rev. C 95, 054914 (2017)

縦初期揺らぎ⇒衝突軸方向の相関の理解

流体揺らぎ
エントロピーマクロ：
エントロピー最大＝熱平衡状態

散逸

ミクロ：
エントロピー最大の周辺で揺らぐ

揺らぎ

状態

揺動散逸関係⇒揺らぎと散逸は同時に存在

流体揺らぎ⇒QGPの輸送的性質の解明

流体揺らぎの影響
揺動散逸関係流体揺らぎ 𝛿𝜋 2 = 4𝑇𝜂/Δ𝑡Δ𝑥 3K. Murase talk at heavy ion café (2019)

𝜂: ずれ粘性
𝑇: 温度Δ𝑡, Δ𝑥: 典型的なスケール

コップの水𝜂 ≅ 10−3 Pa sec
𝑇 ≅ 300 K
Δ𝑥 ≅ 10−3 mΔ𝑡 ≅ 10−1 sec

重イオン衝突反応で生成されるQGP𝜂 ≅ 0.1𝑠 ≅ 0.2 fm−3
𝑇 ≅ 300 MeV
Δ𝑥 ≅ 1 fm
Δ𝑡 ≅ 1 fm

𝛿𝜋 ≅ 4 × 10−8 Pa
𝑃 ≅ 105 Pa

𝛿𝜋 ≅ 2 × 102 MeV/fm3
𝑃 ≅ 4 × 103 MeV/fm3

𝛿𝜋/𝑃 ≅ 4 × 10−13𝛿𝜋/𝑃 ≅ 𝟎. 𝟎𝟓

流体揺らぎは無視できない！

研究目的
QGPの物性の探究
本研究の目的:
流体揺らぎ、初期揺らぎが高エネルギー重イオン衝突実験における測定量に与える影響を解析
QGPの輸送的性質の理解

目次
1. はじめに2. 高エネルギー重イオン衝突反応における測定量– 多重度と中心度– フローの方位角異方性– 因子化比3. 統合的動的模型
4. 結果
5. まとめ多重度と中心度
多重度：
衝突によって生成されたハドロン数多重度
生成エントロピー
衝突反応領域

大大
円形多重度
生成エントロピー
衝突反応領域

小小
楕円形

多重度 最小・・・

最大

中心度 100%

・・・

0%

ALICE Collaboration, Phys. Rev. C 88 (2013) 044909

方位角異方性
衝突軸∞𝑑𝑁∝ 1 + 2 ෍ 𝑣𝑛 cos 𝑛 𝜙 − Ψ𝑛𝑑𝜙

𝜙: 方位角Ψ𝑛 : 事象平面角𝑣𝑛 : フロー係数

𝑛=1𝜙

Ψ2𝑑N𝑑𝜙0
https://www.zmescience.com/scienc
e/physics/quark-gluon-plasma-lhc26052011/

cos 2 𝜙 − Ψ22𝜋𝜙

J. Y. Ollitrault, Phys. Rev. D 46, 229 (1992).

方位角異方性∞𝑑𝑁∝ 1 + 2 ෍ 𝑣𝑛 cos 𝑛 𝜙 − Ψ𝑛𝑑𝜙

衝突軸𝑛=1Ψ2𝜙

𝜙: 方位角Ψ𝑛 : 事象平面角𝑣𝑛 : フロー係数

Ψ2𝑑N𝑑𝜙cos 2 𝜙 − Ψ2

大きな𝑣2(楕円型フロー)

02𝜋𝜙

圧力勾配大→フローを生成
流体模型で𝑣2 の記述が可能→粘性の小さいQGP流体の発見

ラピデティー方向のフローの相関
擬ラピディティー𝜃𝜂𝑝 = −ln tan2

𝜂𝑝1

𝜂𝑝1

𝜂𝑝2
𝜂𝑝𝑣𝑛 𝜂𝑝1 ≈ 𝑣𝑛 𝜂𝑝2
Ψ𝑛 𝜂𝑝1 ≈ Ψ𝑛 𝜂𝑝2
ラピディティー方向の相関あり

𝜂𝑝2
𝜂𝑝𝑣𝑛 𝜂𝑝1 ≠ 𝑣𝑛 𝜂𝑝2
Ψ𝑛 𝜂𝑝1 ≠ Ψ𝑛 𝜂𝑝2
ラピディティー方向の相関喪失⇒事象毎の初期揺らぎ・流体揺らぎの影響？

フローの因子化比−𝜂𝑝𝑎

𝜂𝑝𝑎

𝜂𝑝𝑏 ：基準𝜙 −𝑎𝜂𝑝𝜙𝑎𝜙𝑏相関?
𝑟𝑛 𝜂𝑝𝑎 , 𝜂𝑝𝑏 =𝑉𝑛Δ −𝜂𝑝𝑎 , 𝜂𝑝𝑏
𝑉𝑛Δ 𝜂𝑝𝑎 , 𝜂𝑝𝑏

𝑟𝑛 𝜂𝑝𝑎 , 𝜂𝑝𝑏 ~1
相関あり𝜂𝑝𝑎𝑎 , 𝜂 𝑏 = ⟪cos 𝑛Δ𝜙 ⟫𝑉𝜂, 𝑛Δ 𝑝 𝑝,Δ𝜙 = 𝜙 𝑎 − 𝜙 𝑏𝑟𝑛 𝜂𝑝𝑎 , 𝜂𝑝𝑏 < 1
相関喪失

因子化比の減少
相関の喪失の要因⇒初期揺らぎ？
⇒流体揺らぎ？

CMS Collaboration, Phys. Rev. C 92, 034911 (2015)
F.G. Gardim et al., Phys. Rev . C 87, 031901 (2013)

目次
1. はじめに2. 高エネルギー重イオン衝突反応における測定量
3. 統合的動的模型
4. 結果5. まとめ高エネルギー重イオン衝突反応の時空発展
衝突前

衝突直後

測定

ハドロン化

http://wl33.web.rice.edu/research.html

QGP

光速近くまで加速
ローレンツ収縮すり抜け
カラーフラックス
チューブ生成QGPは数fm/𝑐でハドロン化

測定量：
ハドロンの種類や運動量ハドロンの分布からQGPの性質を調べる方法？
⇒統合的動的模型

統合的動的模型
時間

4. ハドロンガス

T. Hirano et al., Prog. Part. Nucl. Phys. 70,108 (2013)

ハドロンカスケード (JAM) Y. Nara et al., Phys. Rev. C61, 024901 (2000)

3. 粒子化Cooper and G. Frye,
クーパー・フライ公式 (𝑇sw = 155 MeV) F.Phys.
Rev. D10, 186 (1974)

2. (3+1)次元QGP流体

衝突軸

流体揺らぎ
(3+1)次元揺動流体模型
EoS: s95p-v1.1 (格子QCD+共鳴ハドロンガス)
ずれ粘性: 𝜂Τ𝑠 = 1/4𝜋
緩和時間： 𝜏𝜋 = 3/4𝜋𝑇

1. 初期状態
縦: PYTHIA x 修正BGK

K. Murase, Ph. D thesis, The Univ. of Tokyo (2015)

縦初期揺らぎT. Hirano et al., Phys.Lett.B 636, 299 (2006)
M. Okai et al., Phys. Rev. C 95, 054914 (2017)

横平面: MCグラウバー模型

流体方程式
エネルギー運動量保存則𝜇𝜈

𝜕𝜇 𝑇 = 0𝑇𝜇𝜈 = (𝑒 + 𝑃)𝑢𝜇 𝑢𝜈 − 𝑃𝑔𝜇𝜈 + 𝜋 𝜇𝜈𝑇𝜇𝜈 ： エネルギー運動量テンソル𝑒: エネルギー密度,𝜋 𝜇𝜈 : 応力テンソル
𝑢𝜇 : 4元流速（Landau系）, 𝑃: 圧力

状態方程式 𝑃 = 𝑃(𝑒)格子QCD + 共鳴ハドロンガス
P. Huovinen and P. Petreczky, Nucl. Phys. A837, 26 (2010)

状態方程式

状態方程式 𝑃 = 𝑃(𝑒)格子QCD + 共鳴ハドロンガス
P. Huovinen and P. Petreczky, Nucl. Phys. A837, 26 (2010)25

流体揺らぎ
応力テンソルに対する構成方程式𝜇𝜐𝜏𝜋 ∆𝛼𝛽 𝑢𝜆𝜕𝜆 𝜋𝛼𝛽4+ 1 + 𝜏𝜋 𝜕𝜆 𝑢 𝜆 𝜋 𝜇𝜈 − 2𝜂∆𝜇𝜐 𝛼𝛽 𝜕 𝛼 𝜋 𝛽 = 𝛿𝜋 𝜇𝜈
3

緩和時間を含む項
𝜂: ずれ粘性𝑢𝜇 : 4元流速
𝜏𝜋 : 緩和時間∆𝜇𝜈 𝛼𝛽 : 射影演算子
1 𝜇 𝜈1 𝜇𝜈𝜇 ν𝜇𝜈∆ 𝛼𝛽 = Δ 𝛼 Δ 𝛽 + Δ 𝛽 Δ 𝛼 − Δ Δ𝛼𝛽23𝜇𝜈𝜇𝜈𝜇Δ = 𝑔 − 𝑢 𝑢𝜈

流体揺らぎ

構成方程式と揺動散逸関係4𝜏𝜋 ∆𝜇𝜐 𝛼𝛽 𝑢 𝜆 𝜕𝜆 𝜋 𝛼𝛽 + 1 + 𝜏𝜋 𝜕𝜆 𝑢 𝜆 𝜋 𝜇𝜈 − 2𝜂∆𝜇𝜐 𝛼𝛽 𝜕 𝛼 𝜋 𝛽 = 𝛿𝜋 𝜇𝜈3エントロピー増加

エントロピー減少バランス

揺動散逸関係 = 熱平衡状態の安定性条件
𝛿𝜋𝜇𝜈 (𝑥)𝛿𝜋 𝛼𝛽 (𝑥 ′ ) = 4𝑇𝜂Δ𝜇𝜈𝛼𝛽 𝛿 (4) 𝑥 − 𝑥 ′
𝛿

(4)

𝑥−𝑥′

1
⇒
Δ𝑡

13𝑒−2′𝒙−𝒙4𝜆2

4𝜋𝜆2 2𝜆: ガウシアン幅 揺らぎの強さをコントロール

初期条件 修正BGK模型𝑠0 𝜏 = 𝜏0 , 𝑟⊥ , 𝜂s𝜌part 𝜂𝑠 , 𝑟⊥1−𝛼∝𝜌part 𝜂𝑠 , 𝑟⊥ + 𝛼𝜌coll 𝑟⊥
2𝑌𝑏 − 𝜂𝑠
𝑌𝑏 + 𝜂𝑠=𝜌part,A 𝑟⊥ +𝜌part,𝐵 𝑟⊥
𝑌𝑏𝑌𝑏𝛼: パラメータ𝜌part : 衝突に関与した核子数密度
𝜌coll : 衝突回数密度

𝜏0 = 0.6 fm
𝑌𝑏 ：ビームラピディティー
T.Hirano et al., Phys.Lett.B 636, 299 (2006).

初期条件 修正BGK模型
𝜌part 𝜂𝑠 , 𝑟⊥

𝑌𝑏 − 𝜂𝑠
𝑌𝑏 + 𝜂𝑠=𝜌part,A 𝑟⊥ +𝜌part,𝐵 𝑟⊥
𝑌𝑏
𝑌𝑏Rapidity trapezoid を取り入れた形

𝜌part,A + 𝜌part,B2

𝜌part,A𝜌part,B𝑌−𝑌beam

A. Adil, M. Gyulassy, Phys. Rev. C 72, 034907 (2005)

0𝑌beamPYTHIA x 修正BGK模型

イベントジェネレーター
PYTHIA𝜂𝑠

𝜂𝑠

紐破砕をモデル化
粒子生成の分布を考慮

𝑥𝜂𝑠

T.Hirano et al., Phys.Lett.B 636, 299 (2006)𝑥𝜂𝑠T. Sjöstrand et al., Comput. Phys. Commun. 191, 159 (2015)
M. Okai et al., Phys. Rev. C 95, 054914 (2017)

PYTHIA x 修正BGK模型PYTHIA

修正BGK

p-p衝突における縦方向の揺らぎ

𝑁part スケーリング
𝑁A + 𝑁𝐵
2𝑑𝑁 pp𝑑𝑌

PYTHIA

𝑁A

𝑁B
𝑌−𝑌beam0𝑌beam

𝑌
T. Sjöstrand et al., Comput. Phys. Commun.
191, 159 (2015)

M. Okai et al., Phys. Rev. C 95, 054914 (2017)

PYTHIA x 修正BGK模型
MC-グラウバー模型
衝突に関与した核子数: 𝑁A (𝒙⊥ ), 𝑁B (𝒙⊥ )衝突回数:𝑁coll (𝒙⊥ )PYTHIA
生成ハドロン 𝑁coll × pp collisions
𝑁A

𝑁B

棄却化• 低 𝑝𝑇 : 𝑁part にスケール• 高 𝑝𝑇 : 𝑁coll にスケール−𝑌beam  ラピデティー → 𝑌beam
S.J.Brodsky, J.F.Gunion and J.H.Kuhn, Phys.Rev.Lett.39, 1120 (1977)
T.Hirano et al., Phys.Lett.B 636, 299 (2006)

M. Okai et al., Phys. Rev. C 95, 054914 (2017)

初期条件
𝑠0 𝜏0 , 𝜂s , 𝑥⊥𝐾= ෍
𝜏0𝑖12𝜋𝜎𝜂

1exp −22 2𝜋𝜎⊥2𝑖𝑥−𝑥+2𝜎⊥ 22𝑖𝑦−𝑦−2𝑖𝜂s − 𝜂s2𝜎𝜂 2

規格化
𝐾 = 4.8 for 𝑠𝑁𝑁 = 2.76 TeV

ガウシアン幅𝜎⊥ = 0.3 fm𝜎𝜂 = 0.3

核子の位置

初期時刻
𝜏0 = 0.6 [fm]𝑥 𝑖 , 𝑦 𝑖 , 𝜂s 𝑖
MC Glauber + PYTHIA + 修正 BGK

33

初期条件 エントロピー密度分布の例
修正BGK模型PYTHIA x 修正BGK模型

𝑥𝑥
𝜂𝑠

𝜂𝑠

縦初期揺らぎ無縦初期揺らぎ有34

流体模型 エントロピー密度分布の時間発展
粘性流体模型

揺動流体模型𝑥𝑥
𝜂𝑠

𝜂𝑠

縦初期揺らぎ無𝜂𝑠 方向に滑らか𝜂𝑠 方向に凹凸

流体模型 エントロピー密度分布の時間発展
粘性流体模型

揺動流体模型𝑥𝑥
𝜂𝑠

𝜂𝑠

縦初期揺らぎ有大きな違いはない？→定量的な評価が必要

目次
1. はじめに2. 高エネルギー重イオン衝突反応における測定量
3. 統合的動的模型
4. 結果– 生成粒子数の中心度依存性– 因子化比𝑟𝑛 𝜂𝑝𝑎 , 𝜂𝑝𝑏– ルジャンドル係数
5. まとめ生成粒子数の中心度依存性
縦初期揺らぎ無Pb+Pb, 𝑠NN = 2.76 TeV

縦初期揺らぎ有Pb+Pb, 𝑠NN = 2.76 TeV

ALICE Collaboration,
Phys.Rev.Lett.106, 032301 (2011)

実験の生成粒子数を再現するように初期条件のパラメータを調整
同じパラメータセットで因子化比を評価

縦初期揺らぎ無𝑝𝑇 スペクトル

縦初期揺らぎ有低横運動量：実験結果をよく再現
高横運動量：流体模型が実験より小⇒JET等の効果
揺動流体模型＞完全・粘性流体

39

𝑣2 の𝑝𝑇 依存性
縦初期揺らぎ無完全流体
実験データより大きい
粘性・揺動流体(𝜂Τ𝑠 = 1Τ4𝜋)→ 𝑝𝑇 ~1.5GeV 以下でALICEデータを再現
ALICE Collaboration,
Phys. Rev. Lett. 116 (2016) 132302

𝑣2 の𝑝𝑇 依存性
縦初期揺らぎ有粘性・揺動流体(𝜂Τ𝑠 = 1Τ4𝜋)→ 𝑝𝑇 ~1 GeV 以下でALICEデータを再現

ALICE Collaboration,
Phys. Rev. Lett. 116 (2016) 132302

𝑎 𝑏因子化比 𝑟2 𝜂𝑝 , 𝜂𝑝
縦初期揺らぎ無𝑠NN = 2.76 TeV

3.0 < 𝜂𝑝𝑏 < 4.01 ≈ 粘性流体 > CMS data ≈ 揺動流体
粘性流体揺らぎ

CMS, Phys. Rev. C 92, 034911 (2015).

相関を喪失させない
因子化比の減少 相関の部分的喪失

𝑎 𝑏因子化比 𝑟2 𝜂𝑝 , 𝜂𝑝
縦初期揺らぎ有𝑠NN = 2.76 TeV

3.0 < 𝜂𝑝𝑏 < 4.01 > 粘性流体 > CMS data ≈ 揺動流体
縦初期揺らぎCMS, Phys. Rev. C 92, 034911 (2015).

因子化比の減少 相関の部分的喪失

流体揺らぎ + 縦初期揺らぎ

因子化比のラピディティ依存性の再現

因子化比 𝑟2 𝜂𝑝𝑎 , 𝜂𝑝𝑏 の中心度依存性
縦初期揺らぎ無Pb+Pb, 𝑠NN = 2.76 TeV

縦初期揺らぎ有2.0 < 𝜂𝑝𝑎 < 2.5, 3.0 < 𝜂𝑝𝑏 < 4.0

Pb+Pb, 𝑠NN = 2.76 TeV

流体揺らぎ + 縦初期揺らぎ

因子化比 𝑟2 (𝜂𝑝𝑎 , 𝜂𝑝𝑏 )の中心度依存性を再現

因子化比の理解

Jiangyong Jia, Peng Huo, Phys. Rev. C 90, 034905 (2014)𝑣𝑛 𝜂𝑝𝑎 ≠ 𝑣𝑛 𝜂𝑝−𝑎→ 𝑟𝑛 減少Ψ𝑛 𝜂𝑝𝑎 ≠ Ψ𝑛 𝜂𝑝−𝑎→ 𝑟𝑛 減少

フロー係数の相関 + 事象平面角の相関を評価

𝜂𝑝𝑎CMS, Phys. Rev. C 92, 034911 (2015)

線形を仮定
Ψ𝑛 𝜂𝑝𝑎 = −Ψ𝑛 𝜂𝑝−𝑎

• フロー係数・事象平面角のどちらの相関が喪失しているのか？
• 線形だけでなく高次の構造？

ルジャンドル係数

Jiangyong Jia, Peng Huo
Phys. Rev. C 90, 034905 (2014)∞𝑣2 𝜂𝑝= ෍ 𝑎2𝑘 𝑃𝑘
𝑘=0𝜂𝑝,
𝜂𝑝max∞

Ψ2 𝜂𝑝= ෍ 𝑏2𝑘 𝑃𝑘
𝑘=0𝜂𝑝,
𝜂𝑝max

𝜂𝑝max = 2.5

𝑃𝑘 ： ルジャンドル多項式
1𝑃1 𝑥 ＝𝑥, 𝑃2 𝑥 = 2 3𝑥 2 − 1

ルジャンドル係数 𝑎2𝑘 , 𝑏2𝑘⇒ 𝑣2 𝜂𝑝 とΨ2 𝜂𝑝 の相関の喪失を独立に評価⇒ 𝑣2 𝜂𝑝 とΨ2 𝜂𝑝 の非線形構造も評価可能

縦初期揺らぎなし

ルジャンドル係数
楕円型フロー 𝑣2
𝐴12 =(𝑎12 )2２次の事象平面角 Ψ2

𝐵21 =(𝑏21 )2

𝑎2𝑘 , 𝑏2𝑘 : ルジャンドル係数
揺動流体 ＞ 粘性流体

Pb+Pb, 𝑠NN = 2.76 TeV流体揺らぎ⇒楕円型フロー 𝑣2⇒ ２次の事象平面角 Ψ2両方の相関の部分的喪失

ルジャンドル係数
縦初期揺らぎなしΨ2 𝜂𝑝

Ψ2 𝜂𝑝′

Ψ2 𝜂𝑝′′

𝐵22 ≠ 0×線形（Ψ𝑛 ∝ 𝜂𝑝 )⇒ルジャンドル係数：高次の構造を評価可能

𝜂𝑝縦初期揺らぎのルジャンドル係数への影響？

縦初期揺らぎ有 ＞ 縦初期揺らぎ無縦初期揺らぎのルジャンドル係数への影響

縦初期揺らぎ有 ＞ 縦初期揺らぎ無縦初期揺らぎのルジャンドル係数への影響
Ψ2 𝜂𝑝

Ψ2 𝜂𝑝′

Ψ2 𝜂𝑝′′

𝐵22 ≠ 0縦初期揺らぎ有≃縦初期揺らぎ無𝜂𝑝まとめと今後の展望◆統合的動的模型• 縦初期揺らぎ• 流体揺らぎ◆因子化比• 流体揺らぎ⇒相関をより喪失• 縦初期揺らぎ⇒相関をより喪失• 流体揺らぎ 𝜆 = 2.0 fm + 縦初期揺らぎ ⇒ 𝑟2 𝜂𝑝𝑎 , 𝜂𝑝𝑏 の中心度依存性を再現• 流体揺らぎと縦初期揺らぎの重要性まとめと今後の展望◆ ルジャンドル係数• フロー係数の相関と事象平面角の相関を独立に評価する方法を提案• 流体揺らぎはフロー係数の相関と事象平面角の両方の相関をより喪失• 縦初期揺らぎはフロー係数の相関をより喪失
高エネルギー重イオン衝突反応における流体揺らぎの影響を初めて定量的に評価

今後の展望• 流体揺らぎ＋縦初期揺らぎを取り入れた模型で輸送係数を評価