雑誌名
筑波大学計算科学研究センター 平成24年度 年次報
告書
平成
24
年度 年次報告書
Center for Computational Sciences,
University of Tsukuba
1
目
1挾 「」 挾 敲策 改善目標挾普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普挾「挾 1責1普 敲策挾普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普普挾「挾 1責「 改善目標 法人評価 証評価 部評価 び 査･ 部 査 指摘
2
1
成
23
度
点施策
改善目標
1-1.
点施策
第 期中期計 関わ 大学全体 実行計 い 敲策
22-1 各研究 や研究 研究戦略 企 組織 実 強 当 組織 学問
特色 生 長期的展望 立 礎研究 学 融 的 研究 計 的 推逭
24 敲策 各系 研究 い 系等 特色 生 研究戦略
び研究連携策 推逭
32-1 共 利用 共 研究 戦略経費等 支援 各 第 期中期
目標 中期計 期間中 目標 計 期間中 評価 実敲 実敲状況 検証 行い い研究 実現
24 敲策 各共 利用 共 研究 国 研究 向 体制構築
係 活動 実敲
研究
挾学 計算 学 計算 礎 計算
能力 び機能 飛躍的 向 計算技術開 計算 学 革新 組 計算技術開拓 端学 計算 学教育研究 実 い 研究開 実敲
挾共 利用 共 研究 端学 計算 学共 研究 活動 特 経費 端学 計算 学 開拓 推逭 展開 業 学 共 研 究 実敲 学 計算 学 共 研究 実敲 計算 学 研 究 推逭
挾KEK 国立 文 遀営 計算 礎 学連携 い HPCI 戦略 計算 礎 学 実敲 国 計算資源 効利用
HPCI 参 国 計算 学 遀営 献
向 検討 逭 国 計算 計 献
挾 全体 研究推逭 業 大型 等 中心 計算 学
3
素粒子物理研究部門 挾 京 格子QCD 物理 間相 作用 決 計算 向 改良や 間相 作用 計算方法 研究 行う 限 密 格子QCD 研究 行う
原子 物理研究部門 挾 初 体 銀河 巨大 形
大規模 輻射流体力学 体 究 多
系 子 大規模計算 元素 関わ 原子 応や
物質 相 作用 関 理解 展
子物性研究部門 挾 第一原理 子論等 構造体 持 特 的物性 解明 世
技術 創 等 研究 行い 参 い 産業界
献 計算 学 知見 得 生 等 連携 生物学
資 知見 第一原理 子論 獲得
生 学研究部門 挾 生 学 要 課 い 計算機 理論
的解明 GPU 用い 子力学計算手法 開 行い 生 学研究 い 必
要 膨大 計算 応 う 子逭 生
物群 近縁 新 真 微生物 索 網羅的 現遺伝子解析 行う
地球環境研究部門 挾 大気大 環 NICAM 用い 極 気 研 究 推逭 極海海 激浿 関係 解明 領域気象 WRF 建物解像
LES 開 都 気候 研究 実敲
性能計算 研究部門 HA-PACS 並列言語XMP び 演算 速向
張 密結 演算 速機構 実 評価 逭 域 散
Gfarm HPCI実遀用 通 評価 改良 行う 並列数値 一層
大規模 速 演算 速機構向 逭
計算情報学研究部門 挾 大規模計算 関わ 技術 敁備 共 実
旓間実世界 利活用 い 応用的側面 技術 深 人 型
精 解析 複 現実感型 通 行え う
計算 効性 検討
4
国 連携 大学 び 米国 研究所 連携
逭 韓国KISTI 連携 逭 連携 計算機 国 共
研究 活用 う 手 確立
教育
計算 学 研究 実敲
学生 数 増や 30 計算 学 英語 講義 計算 学 関 大学院
共通 目 実敲 国 逭
1-2
改善目標
法人評価
認証評価
部評価及び
査･
部
査の指摘
対す
組等
部評価 い 指摘 計算機開 端 計算 学 推逭 学 計算 学 推逭体制 い 特 経費 び HPCI 戦略 等 確
実 推逭 あ 課 T2K後 開 調遉 い
5
2
成
24
度
実績報告
2-1.
成
24
度計
並び
改善目標
記載さ
い
項
い
の
成状況
等
挾 計算 学 計算機 学 働 学 計算 学 強 特 経
費 計算技術開拓 端学 計算 学教育研究
実 実敲 い い 開 演算 速機構(GPU)間 直接結
い 設計 実装 了 産品 調遉 終え 着実 逭行 い 東 大 理研 会津大 斌立 将来 HPCI あ
方 関 調査研究 い 国 い 検討 逭
目指 計算 学 推逭 い T2K 後
い 東京大学情報 柏 共 期
設置 遀用 組織 端共 HPC 敲設 設置 遀営
い 締結
挾 共 利用 共 研究 共 研究 主要計算設備 活用
学 共 実敲 48課 択 実敲 38件 旅費
支援 1件 会支援 設置20周 念 第3回 学
計算 学 新 知 見 統 創出 9 開催 長期的
礎研究 含 各研究 行う 課 い 学 共
実敲 研究 着実 逭 い
国 的 連携 米国 研究所 欧
大学 力関係 逭 い 人 育 い 計算 学英語
推逭や計算 学 海 支
援 行い 国 的 活動 人 育 組 い
敲策 遉 状況 研究
1 24 2 遀用 始 HA-PACS い 当 研究
供 10 学 共 利用 組 入 用い
共 研究 開始 一 要 研究 目 あ 演算 速
機構(GPU)間 直接結 PEACH2 い 設計 実装 了 産品 調遉 終
え い 組 TCA部HA-PACS 向
調遉作業 逭 い い 昨 設置 端計算 学推逭
6
2 共 利用 共 研究 共 研究 主要計算設備 あ
T2k-tsukuba FIRST 10 HA-PACS え 利用
学 共 実敲 学 計算 学 研究 推逭 当 い 48
課 択 実敲 38 件 旅費支援 件 会支援 行
設置20周 念 第3回 学 計算 学 新 知 見
統 創出 9 開催
3 KEK 国立 文 締結 計算 礎 学連携 遀営 中心 計算
礎 学 HCPI 戦略 推逭 HPCI 参
学長 理 長 え 組織遀営 あ い HPCI
構築 い 域 担当 引 散
Gfarm 用い 国 計算資源 効利用 HPCI 共 敁備
行 7 中心 東 大 理研 会津大 斌立
将来 HPCI あ 方 関 調査研究 択
国 い 検討 逭 い
4 全体 研究推逭 業 大型 等 中心 計算 学 学
的研究 遂行 実現 各部門 遉 状況 あ
素粒子物理研究部門 京 配 生 改良
了 配 生 格的 始 間相 作用
関 計算方法 改良 使 研究 逭 京
移植や京 作業 終了 限密
格子QCD 法 reweighting法 質 い
領域 相構造 解析 了 挾 挾
原子 物理研究部門 初 体 銀河 巨大 形
大規模 輻射流体力学 体 究 多 系
子 大規模計算 元素 関わ 原子 応や 物
質 相 作用 関 理解 展 挾 挾
子物性研究部門 第一原理 子論等 構造体 特 的物性 解明 行う 旓 産業界 献 知見 獲得 目指 研究 行 結
物性等 礎物理学 い 大 知見 得
7
い ワ TPEC 参
い 産業界 大 献 計算 学的知見 得 世 ワ
設計指針 獲得 挾 挾
生 学研究部門 生 学 要 課 解決 計算機
理論的解明 GPU 用い 子力学計算手法 開 行 特 生
連携 生物学 い L型 酸遃 起源 や
系学惑斴 能性 等 子論的知見 得
子逭 生物群 含 類*Roombia* sp.
生物群 含 Dysnectes brevis 載種NY0171株 網羅的 現遺伝子 行
い 使用 大規模 子系統解析 行う 新 生物種 系統
的 置 確 挾 挾
地球環境研究部門 2008 6 極海 生 極 気 水 解像 14km
(gl-9) NICAM 用い 再現実験 行い 極 気 遉 維持 関 細 析
結 査 付 国 会議 表 極海 海 浿少 伴い 極温 増幅
温 場 第2EOF 示 一方 第1EOF 極
振動 示 明 国 掲載 第 回国
極研究 (ISAR-3) 表 都 気象LES用 街路樹 開 街路
樹 中心断面 温 再現実験 行 領域気象 WRF-UCM 力学的
行い IPCC SRES A2 大都 圏 夏季気候 将来予
測 実敲
性能計算 研究部門 並列言語XMP 演算 速装置向 張 HA-PACS 評
価 共 GPU CPU 全 並列 無駄 利用 機構
組 実装 評価 密結 演算 速機構
開 終え 予備性能評価 来方法 比 通信遅延 大幅 浿
確 研究 終 向 産 終え 計算 学
研究 共 物理学 気象 等 HA-PACS
移植 速 逭 Gfarm 関 HPCI実遀用 支援 性能 検
出 改善 全国規模 遀用 推逭 並列数値 い FFT
改善 京 HPC Challenge 再 世界第一
部 空間法 改善 QCD 始 適用
8
計算情報学研究部門 中心 情報統 技術 深
共 GPU 活用 大規模計算 関わ 知識 見やXML
Web 等 技術 関 研究 推逭 GEOGid
利活用やX線 文 検索等 応用的 研究開 遂行 計算
斵像 像 い 解析 効率 行え 方法論 構築
特 複 現実感型 通 人 型 精 解析法
多種 人物 表情 元的 情報や指示 共
研究 計算 効性 検証
国
挾 3 7～8斌 研究所 ワ 方 行い 連携 状
況 い 行 共 研究 い 情報交換 研究 力
全般 MoU 締結 見 大学 筑波 行
う予 あ 方 状況 調敁 降 見 あ
韓国KISTI 連携 関 方 組織改 調敁中 あ 計算機
国 共 研究 海 利用 関 敁備 行い 手 遀用
い 当 海 利用件数 件 あ
教育
挾 計算 学 関 念 学生
数 増や 新 開設 生 環境研究
情報研究 計算 学 博士課程 学生
入 計算 学 籍 学生 い
RA 用 海 支援 い 2人
支援 計算 学英語 え
応 英語 開始 英語 計算 学 教 育 実
改善目標 遉 状況
挾 T2K 後 い 東京大学情報 柏
共 期 設置 遀用 計 い 組織
端共 HPC 敲設 東大 書 び覚書 準備
締結 行 当 敲設 組織敁備 行い 設計 調遉準備
体的 遀営 逭 い 期 設計 研究開 準備
9
波大学 敁備 い
期 設置 間 予算 26 増
強 端共 HPC 敲設 設置 期 設置
遀用体制 検討 い 現状 T2K い 当初
見直 26 2 遀用 等 調遉
移行 準備 行う
特色あ 組 実敲状況
(1) 共 JLDG 遀用 散 Gfarm 研究開
主 計 算 素 粒 子 物 理 学 研 究 者 共 JLDG(Japan
Lattice Data Grid) 遀営 い 現 主 HPCI戦略 支援
遀営 い 組 十 渡 中心 組
あ 昨 端 敁備 支援 現 7 機関 20 結ぶ
総 2.6PB 供 い
支え あ 散 Gfarm
開 行 HPCI 実遀用 供 い
(2) 部門間 連携研究 推逭
当 学諸 計算機 学 情報 学 連携 働 学 計算
学 中心的 研究活動 行 い 素粒子 超 速計算
地球環境 計算知能 体的 研究課 い ワ 設置 期開催 行い 共 研究 逭 い
計算 学 目指 特 経費 逭 部資金 獲得 い
組 い
2-2.
自己評価
課題
1 自己評価
特 経費 計算技術開拓 端学 計算 学教育研究
実 終的 目標 密結 並列演算 速機構 持
構築 向 着実 逭捗 共 利用 共 研究 学 共 利
用 い 着実 遀営 い 期 い 東
京大学情報 端共 HPC 敲設 設置 柏
共 期 設置 遀用 計 あ い 締結 来
計 2015 遀用 計 稼
10
全国的観 現 HPCI 戦略 計算 礎
学 戦略機関 推逭 HPCI 遀営 中
心 HPCI 構築 大 献 い 将来 HPCI
あ 方 関 調査研究 開始 国
開 献 い 計算 学 や計算 学
英語 人 育 国 連携 逭 い い
2 課
生 学や 学 研究体制 敁備
評価 査 指摘 研究 実績 積
研究領域 中心 推逭 生 学や 学 業 連携 活 推逭
方策 検討 必要 あ 研究者 確保や 力体制
う い 検討 逭
務支援体制 実
実 T2k-tsukuba え HA-PACS 遀用開始 世
戦略 実敲 伴い 研究企 務 総務 報等 面
検討 処理 務 大幅 増大 特 世 戦略
い 務 務組織 立 遀営 い 共 利用 共
研究 形 強 業 支援 計算機 遀営 人員 用
改善 い 依然 務体制 人員 足 改善 必要
あ
省電力 組 予算面 課
震災後 電力 情 変 電気料金 騰 予想 い 電気料金
予算 大 部 占 部 あ 騰 場 処 検討 必要
あ 部資金 間接経費 配 方式 全学的 研究 敁備 中 研
11
3.
各研究部門の報告
I.
素粒子物理研究部門
1.
教授 青木 慎也 計算 学研究 金谷 和 共
研究員
准教授 挾 挾 挾 挾 石 人 藏増 嘉伸 照 石 理 英克
講師 谷 裕
研究員 挾 挾 挾 挾 浮 尚哉 々木 健志 滑 裕
2. 概要
挾 当部門 格子QCD 大型 研究 活 研究活
動 行わ 格子場 理論 研究者 大部 2006 7 計算 学研
究 期並列計算機 PACS-CS 入 契機 新 立
研究 PACS-CS Collaboration 参 い PACS-CS
Collaboration 計算 学研究 PACS-CSやT2K-Tsukuba 主要 計
算機資源 QCD 関 近似 い物理的予言 行う 目的 3種
類 up down strange 軽い 物理的質 物理 い 動的
扱うNf=2+1 QCD 大規模 逭 更 展的課
up down 間 質 差や電磁相 作用 入 Nf=1+1+1 QCD 研究や
格子QCD He原子 研究 推逭 PACS-CS Collaboration 2011
9 PACS-CS機 伴 解散 組
研究課 HPCI戦略 後述 研究開 課 展的
引 い 並行 速器研究機構 KEK
共 用い 研究 推逭 い JLQCD Collaboration
や 子間 や 間相 作用 研究 組 い HAL QCD
Collaboration 限温 限密 QCD 研究 推逭 い WHOT-QCD
Collaboration 活 研究 行 格子QCD配 や
共 ILDG/JLDG 構築 敁備 推逭
挾 一方 国 計算 学全体 動向 京 中 革新的
構築 主 High
Performance Computing Infrastructure HPCI 戦略 文部 学省
推逭 い HPCI戦略 戦略 青木
統括 任者 物質 起源 構造 択 2011
格的 活動 始 2012 秋 一般利用 開始 京
中心 研究活動 活 い い活動 容
12
戦略 実敲 機関 青木 長 計算 礎
学連携 あ 活動 http://www.jicfus.jp/jp い
3. 研究
挾PACS-CS Collaboration 活動 全員
挾 当 17 3ヶ 計 特 教育研究経費 形
開 製作 逭 超並列 計算機PACS-CS 計算 数2560
演算性能14.3Tflops 18 7 稼働 開始 PACS-CS
Collaboration PACS-CS 主要 計算設備 格子QCD 研究 行う 目
的 筑波大学物理学域 中心 組織 い 目標 3種類
up down strange 軽い 物理的質 物理 い 動的 扱
う 行い QCD 関 近似 い物理的予言 行う あ
18 19 物理 向 up-down 質 段 的 軽
20 PACS-CS 目標 あ 物理
組 開始 21 reweighting法 用い 物理 直
後 1+1+1 QCD+QED 開
組 up down 質 差 電磁相 作用 reweighting法
入 24 論文 表
挾 PACS-CS 23 9 遀用 終了 24 2 後
機 密結 並列演算 速機構実証 HA-PACS 計算 数268 GPU
部 演算性能713Tflops CPU部 演算性能89Tflops 稼働 開始
1 1+1+1 QCD+QED
挾 来 格子QCD計算 的理 up down 質
人 的 等 2+1 電磁相 作用 効 無視 い
1+1+1 QCD+QED 自然界 再現 up
down strange 質 犉立 扱い 電磁相 作用 効
旓 評価 目指 い 電磁相 作用 びup down 質
差 reweighting法 入 変形block BiCGSTAB
開 行う reweighting法 効率実行 能 論文1
13
質 strange 質 格子間 決 4 物理
π+ 質 K0 質 K+ 質 Ω 質 用
up down 質 差 直接 的 評価 能 (論文2)
1 電磁相 作用 up down 質 差 入 K0 K+
伝 関数 比 旓間 関数 あ 赤線 K0-K+質 差 実験値 K0 K+ 均質 1% 満 期 傾 あ 計算結 誤差
範 実験値 再現 い t=0 近傍 励起状態 寄 実験値
い 計算 PACS-CS機 超並列 計算機T2K-Tsukuba 計
算 数648 演算性能94 Tflops 20 6 稼働開始 利用 い
1 K0 K+ 伝 関数 比 赤線 実験 期 傾
2 研究
挾 物理 直 極 現実 近い 含
検証 び予言 え 論文3 4 含
実験的 精 調 い 方 含
Ξcc(3520) 実験的 見 報告 い 見 報告 い 実験
あ 実験 否 的 結 報告 い
Ξcc 確 い い 計算 1 含
質 実験値 再現 Ξcc 質 値 求
3603(15)(16) MeV 得 値 報告 い 実験値3520 MeV 意 異
0 5 10 15 20 25
t
0.92 0.94 0.96 0.98 1.00 1.02 1.04 1.06 1.08 1.10
14
研究 既 実験報告値3520 MeV 誤 あ 真 値 100 MeV程
い 示唆 Ξcc え び 含 見
質 予言値 2% 精 決
挾HPCI戦略 5 研究開 課 青木 石 藏増 石
谷 浮 々木 滑
挾 5 物質 起源 構造 戦略目標 始 歴
敱 素粒子 元素 斴 銀河形 物質 起源 構造 複数
層 繋 計算 学的手法 統一的 理解 あ 目標 目指 4
研究開 課 設 う 一 あ 格子QCD 物理
間相 作用 決 任者 藏増 務 い 課 目指 格
子QCD計算 微細 鍵 新 い展開 あ 微細
電磁相 作用やup down 質 差 入 大規模
実現 意味 QCD あ 質 精
決 能 方 QCD 用い 子 複数作
子間 効相 作用 調 更 子 束縛状態 あ 原子
直接構 行う 意味 前者 青木 初 理研 中心 HAL QCD
Collaboration 組 い あ 後者 藏増 中心
推逭 い 後述
1 格子QCD 物理 間相 作用 決
挾 24 9 28斌 京 共用 開始 降 配 生
格計算 開始 964 格子 0.1 fm 格子間 用 DDHMC
Domain-Decomposed Hybrid Monte Carlo 用い 2+1
mu=md≠ms QCD 実行 い HMC 人
的 入 擬旓間 子動力学法 組 入 法 一種 あ
格子QCD い 用 い DDHMC法 4
元旓空間格子 領域 割 擬旓間 方程式 層的 積
効率 い 熱 衡状態 遉 確
後 5 子動力学旓間 擬旓間 配 保 配 用
15
強い相 作用 崩壊 い 質 計算結 実験値 比較
あ 既 1%程 統計誤差 到遉 い わ
質 mu=md≠ms 格子間 決 3 物理 π 質
mπ K 質 mK Ω 質 mΩ 用 い m
π mK 更 数%程 実験値 必要 あ 見
補 後reweighting 法 用い 行う
2 質 実験値 比較
挾 研究開 課 作用 smeared link 用い Wilson-clover
作用 用い い 谷 clover 改良係数cSW 非摂動論的 決
い 論文5
挾格子QCD 自 原子 直接構 藏増
藏増 山崎 屋大 2010 世界 初 格子QCD
原子 構 後 2 子系 束縛状態 あ 陽子 構
計算 抑え 近似 物理的 値 い
質 行 計算 あ 真空偏極効
入 2+1 格子QCD計算 張 行 結 近似
2+1 QCD い 原子 陽子 束縛 確
論文4 質 依然 物理的 値 m
π
=0.51 GeV相当 あ 3 4He原子 ΔEL(4He) 空間体積
依 性 表 い 4He原子 底状態 自 4
-0.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14
K N 16
個 子 質 差ΔEL(4He)=EL(4He)−4mN 義 限 空間体積 L3
4個 子 士 散乱 効 ΔEL(4He) 含 束縛
出 空間体積無限大極限 1/L3→0 必要 斴
束縛 実験結 表 々 計算結 程 値 再現 い
わ 3 陽子 ΔEL(3S1)=E(3S1)−2mN 空間体積
依 性 表 い 陽子 場 束縛 計算結 実験結 5
倍程 大 い 々 計算 実験的 束縛 い い 一 1S0
い 束縛状態 確 結 質 物理的
値 い い 推察 い 後 質 更 軽 現実
値 近 い 要 あ
3 L3
空間格子 4He原子 陽子
ΔEL 空間体積依 性 四角 空間体積無限大極限(1/L3→0) 値 斴
束縛 実験結 白抜 近似 計算結 表
挾格子QCD 間相 作用 研究 青木 石 々木
挾 子間 働 力 力 中遠距 引力 近距 強い斥力
実験的 知 い 力 性質 特 斥力芯 近距 斥力
理論的 素粒子原子 物理 大問 あ 青木 石
々木 理 学研究所 初 HAL QCD Collaboration 結 子
間 波動関数 子間 出 いう方法 応用 様々 粒子
間 格子QCD 数値 計算 一連 研
0 110-5
210-5 310-5 410-5 1/L3 -0.024 -0.020 -0.016 -0.012 -0.008 -0.004 0.000 experiment Nf=2+1 m=0.51GeV Nf=0 m=0.80GeV
EL(3S1) [GeV]
0 110-5
210-5 310-5 410-5 1/L3 -0.08 -0.06 -0.04 -0.02 0 experiment Nf=2+1 m=0.51GeV Nf=0 m=0.80GeV
17
究 論文5 あ 論文6 7 斥力芯 起源
演算子積展開 摂動論的QCD 繰 群 いう解析的 手法 研究 格子QCD
計算結 比較 い 論文8 HAL QCD 用い 方法 非弾性散乱
張 方法 案 あ 石 々木 び青木 学生
あ 山 2012 研究 紹
1 力 称 び 称LS力
挾 懸案 あ flavor SU(3) 極限 力 一般形 い
子系 Okubo-Marshak 方法 張 求
Octet-baryon 体系 flavor規 表現 表現 現 粒子 交換
称 物(s) 称 物(a) 類 称LS力
子系 唯一 遊い あ 様 flavor SU(3) 破 状況 一般形 求
結 off-diagonal部 子系 排
い 形 生 逆構 い
扱う い あえ 初 無視 計算 逭
あ 計算 実行 新型 contraction code 開 行
25 早期 段 結 出 現実的 計算 い い
2 2+1 flavor QCD 力
挾 22 23 開 旓間依 型 Schrödinger-like eq. 方
法 4 関数 ground state saturation 必要 い t
間 計算 方法 あ PACS-CS 配
用い 2+1 flavor QCD 力 統計 方法 適用
再解析 行 長距 one pion exchange極限 考慮 AV18型
関数形 用い 中心力 一 力 旓 行
実験的 周期境界 響 う 試 滑
使 求 相差 引力的 絶 値 実験値
遥 束縛状態 い い m
π=700 570 411 MeV 変
領域 実験値 向 引力 強 傾向 見 い い
18
4 2+1 flavor 格子QCD 計算 力 spin-singlet 中心力
結 m
π=570 MeV 使
1S0
散乱 相差 m
π=411 570
700 MeV
3 格子QCD 得 力 少数 子多体系 研究
挾 物理 格子QCD 一般 力 出 備え
Nambu-Bethe-Salpeter(NBS)波動関数 得 一般 力 用
い 散乱 相 物理 出 各種方法 確立 び軽い原子 構造
適用 準備 逭 体的 SU(3) 称 得
力 偶 1S0 び3S1−3D1状態 働 中心力 び 力 用い 4He原子 少数多体問 精密 解 束縛状態
確 5 束縛状態 結 4.4 MeV 陽子間 力
旓 考慮 4体問 解い 場 5.1 MeV 力 考慮
4体問 解い 場 あ 現実 4He 比 非常 い 結 既
行わ い 力 繰 中心力 4体問 解い 近似計算 性
的 良 似 結 あ 回用い 力 π中間子 応
質 469 MeV あ 質 現実 得 力
力 現実 場 比 非常 弱い 考え 4体問 解い 得 4He 底状態 全軌遈角遀動 L=2 状態 期 値 わ
1.3% 力 い 1.4 % 力無 あ 現象論
的 決 現実的 力 い 計算 4He 底状態 L=2状態 期 値
10%程 あ 知 い 物理 力 計算
極 要 あ 示 い
-500 0 500 1000 1500 2000 2500 3000 3500 4000
0 0.5 1 1.5 2 2.5 VC ( 1S 0) [ MeV ] r [fm] -40 -20 0 20 40 60
19
5 SU(3) 称 格子QCD計算 得 力 用
い 得 4He 底状態 子間距 r [fm] 相関関数 陽子間
考慮 い
4 結 手法 用い 間相 作用 研究
挾 S 含 間相 作用 研究 間 近距 斥力
起源 明 要 あ 々木 HAL QCD Collaboration
出方法 結 Schrödinger 方程式 張
含 2体系 行列 系統的 調
間 近距 斥力 消失 伴い現 期 い H-dibaryon 状態 着目
S=−2 間 用い 状態 関 研究 行 挾 計算
PACS-CS Collaboration 生 2+1 配 用い
質 700 MeV 570 MeV 410 MeV 応 3 異 質
い 結 解析 行 Esb1 2 3 ぶ SU(3)
称性 破 大 い 注意 6 準備段 あ
Esb1 2 3 配 得 計算 ΛΛ NΞ 散
乱 相差 示 い 質 軽 SU(3) 破 大
束縛状態 あ H-dibaryon 共鳴状態 移 変わ い 様子 確
出来 後 散乱 相差 誤差 評価 H-diabaryon 様子 細
調 計 い 系 い 解析 束縛状態 無
20
6 Esb1 2 中 3 ΛΛ NΞ 散乱 相差
5 Ω-Ω 間相 作用 研究
Ω 模型 類 10 唯一強い相
作用 あ S=−6 あ Ω-Ω間
HAL QCD Collaboration 力 出方法 用い 計算
10 粒子間 計算 回 初 あ 回
模型 計算 強い引力 あ 予想 J=0 状態 計算
計算 CP-PACS Collaboration 生 2+1 現実 い
質 質 875 MeV 応 配 用い 結 中心 斥力芯
あ 中距 引力 あ 束縛状態 断言 い 強い引力
あ 示唆 い 7
7 J=0 あ Ω-Ω間 V [MeV] 軸 間 距 r [fm]
挾 限温 限密 QCD 研究 青木 金谷
挾 金谷 青木 理 学研究所初 主任研究員 新潟大学 尻准教授 島大学梅
准教授 共 研究 Wilson型 限温 密 Ｑ 研究 引
推逭 WHOT-QCD Collaboration 主要 論文
21
1 固 格子間 T-integral法 状態方程式 研究
挾 20 開 T-integral 法 固 格子間 様々 温
一 格子 実行 方法 あ 特 固 格子
間 既 温 配 活用 あ 計算旓間 大
幅 抑え 精 い 限温 計算 遂行 能性 拓い い 20
近似 試験研究 方法 効性 確 21 現
実的 Nf=2+1 逭 い 第一段 Nf=2+1 QCD
限温 実行 u d 質 現実 い あ
型 初 Nf=2+1 状態方程式 計算 研
究 CP-PACS+JLQCD Nf=2+1 QCD 温 研究結
公開 い 温 配 利用 計算 大幅 削浿
8 論文11 固 格子間 生 Nf=2+1 QCD 配
び 20 試験研究 生 近似 配 使 固 格子間
静的 自 計算 固 格子間 法 自
研究 い 利 あ 示 9 論文12 現 u d
質 物理 PACS-CS 温 配 利用 物理
限温 向 準備 逭 い
8 固 格子間 型 用い 限温 Nf=2+1
QCD 研究 論文11 感 率 U
d 現実 T〜200 MeVあ 限温 温 相当
22
状態方程式 p 物質 力 ε 密 ε−3p
表
9 固 格子間 静的 自 研究 論文12
Nf=2+1 QCD QQ-bar 静的 自
温 T=0 自 測
限温 自 個 相関 評価 固 格子間
数 調敁無 自 距 温 依
い 直接示 出来 端 矢 個 期 値
期 自 遠距 極限 漸近値 遮蔽質 mD(T)
数依 性 縦軸 mD(T)/T 軸 T/Tc Tc 相転移温 Nf=2+1
190 MeV 仮
5 法 QCD 限温 限密 相転移 研究
挾 QGP 限温 限密 相転移 性質や 数 うえ 観測
や 数 義 効 直感的 情報 含 い
方 系 依 性 調 強力 手法 あ reweighting法
系 作用 現 物理 関 必 計算 々 物理
reweighting法 組 わ 相転移 数 容断
解析方法 法 開 論文10 13 初 方法
い領域 QCD 相構造 効
振 舞い 相構造 研究 純 理論 限温 相転移
動的 効 変わ 臨界 置 一般 学
1 2 3
0 0.2 0.4 0.6 0.8 1
F
1(r
, T ) [ G eV ]
r [fm]
Nf = 2 + 1
V(r,T = 0)
T = 232 MeV
T = 279 MeV
T = 348 MeV
T = 464 MeV
T = 696 MeV
1 2 3 4 5 6
1 1.5 2 2.5 3 3.5 4
T / Tc
23
場 計算 限密 相 う振 舞う 評価 結
臨界 置 関 極 い 示 臨界 置
質 学 関数 計算 10
軽い領域 QCD 限密 相構造 解明 向 研究 逭 い
軽 要 役割
部 部 応 行列式 関 効
研究 Nf=2 QCD 場 相 無視 phase-quenched QCD
実行 相 効 再 付 法 評価
学 大 ゆ 転移 転移 変わ
示唆 結 得 10
10 法 QCD 限温 相構造 研究 領域
Nf=2+1 QCD 臨界面 u d 学 μud
s い い場 結 κ 質
逆数 比例 κs=0 ms=∞ κud=0 mu=md=∞ 応 臨界
面 向 う側 限温 QCD転移 相転移 手前 い
軽いNf=2 QCD 相構造 研究 効 微
物理 P F 自 関係 関数
示 相転移
挾 限密 QCD 藏増
挾 限密 QCD研究 非 学 領域 格子QCD計算
必要 学 入 符 問 引 起 う
う 困 克服 藏増 武 金沢大 符 問 解
析的研究 行 体的 行列式 巻 付 展開法 用い 表
0.001 0.001
0.002 0.002
0 200 400 600 800 1000 1200 1400 1.0
1.2 1.4 1.6 1.8 2.0
F
P
0.144139, T 4.0
24
符 問 源 あ 複素 相 解析的 表現 質 展開法 適用
相 限値 格子 空間体積 比例 旓間方向 伸張 指数
関数的 抑制 見出 々 格子 用い 計算
解析的結 検証 論文 纏 論文14 後
数値的 4 QCD 限 解析 行
4 QCD あ 質 領域 い 一 相転移 起 期
示 否 々 計算手法 良い
11 旓間方向 格子 4 固 空間格子 V=63〜103 変
場 様々 物理 感 率 V依 性 表 い
依 性 線形 あ 一 相転移 結論 世界
初 限密 QCD 限 解析 応用 あ 効
性 実証 意義 大 い 現 結 論文 纏 い
11 様々 物理 感 率 空間格子 V 依 性
挾 限温 U(1)A 称性 谷
挾 谷 2個 含 Nf=2 QCD い SU(2)搬SU(2) 称
性 破 い U(1)A 称性 限温 回 能性 い
青木 深谷(大阪大)両氏 共 示 体的 SU(2)搬SU(2) 称
25
式 U(1)A 称性 Ward-Takahashi恒等式 調査 U(1)A 称
性 熱力学的極限 回 示 限温 い 回
称性 種類 限温 相転移 数 強い 響 え 回 結
Nf=2 QCD 相転移 数 一 あ 強 示唆 い 論文15
挾Conformal理論 研究
挾 Conformal理論 素粒子標準模型 越え 理論 候補 あ 理論
関連 興味深 近 多 研究 い 岩崎筑波大 誉教授 中山
Caltech教授 石 島大准教授 共 研究 格子数値
conformal理論 研究 行い SU(N) 理論 Nf個 結 系
Nf conformal window あ 場 理論 dynamics 関 新 い知見
得 infrared cutoff 持 "conformal理論" 1) 質 大
閉 相 "conformal相" 2) 相間 転移
転移 あ 3)"conformal相" 伝 関数G(t) 閉 相 旓間t
関 指数関数型decay form c exp(−mt) 異 指数補 あ 湯 型
decay form c’ exp(−m’t)/t あ 示 論文16
挾JLDG 改良 遀用 浮
挾 国 格子QCD研究者 Japan Lattice Data Grid JLDG
利便性 用性 向 建部筑波大准教授 笠筑波大准教授 び
JLDG 担当者 改良等 1) 東京大学 屋大学
立 2) gfarm 更新 3) 筑波大HA-PACS FE WG群
JLDG 機能 組 4) KEK A FE JLDG
機能 組 5) 管理 遀用 保 支援業務 部委 開始 HPCI
共用 JLDG 連携 構築 技術検討 設計
4. 教育
1 挾 学 論文
1. 斎藤華 博士論文 Finite density QCD in the heavy quark region 質
26
2. 崔 敦 士論文 Phase diagram at finite temperature and chemical potential in strong coupling limit of lattice QCD 限温 び 学
格子QCD 相 強結 極限
2 挾 非常勤講師 中講義
1普
青木慎也
東京大学大学院総
文
研究
客員教授
「01「
『
「01『
」
医
間
週
挾「01「
格子
理論
礎
応用
関
講義
行
挾
「普
青木慎也
中講義
格子
理論入門
立教大学理学部
「01「
盤
責
重
斌挾
5. 賞 部資金 知的 産 等
1 挾 賞
1. 青木慎也 初 哲男 石 理 格子 子色力学 力 研究 第23回
賞 2012 10 17斌 茨城県 学技術振興 団
2. 青木慎也 初 哲男 石 理 格子 子色力学 力 出 2012
仁 念賞 2012 12 6斌 仁 念 団
2 挾 部資金
1. 青木慎也 表 性能汎用計算機 利用 業費補 金 23 択
HPCI戦略 5 物質 起源 構造 492,939,785
2. 青木慎也 表 学研究費補 金 新学術領域研究 20 択 素
融 計算 学 い 層的物質構造 解明 5,900千
3. 金谷和 表 学研究費補 金 研究(B)(一般) 21 択 物
理的 質 限温 限密 Ｑ 格子研究 3,200千
4. 金谷和 担 学研究費補 金 研究(A) 23 択 湯
朝 坂 念 料 析 斌 素粒子物理学者 系譜100 1,000千
27
6. 藏増嘉伸 表 学研究費補 金 研究(A) 22 択 強い相
作用 織 物質形態 QCD 統一的研究 11,300千
7. 藏増嘉伸 担 学研究費補 金 新学術領域研究 20 択
子色力学 真空構造 力学 5,100千
8. 藏増嘉伸 担 戦略的創造研究推逭 業(CREST) 23 択
応 層 超並列固 値解析 開
1,428千
9. 照 担 学研究費補 金 新学術領域研究 20 択
断 計算機 10,000千
10. 谷 裕 表 学研究費補 金 研究(C) 22 択 格子
非摂動論的 繰 実践 650千
11. 浮 尚哉 表 学研究費補 金 若手研究(B) 24 択 物理
格子QCD+QED数値計算 実現 陽子荷電半 計算 1,200千
12. 滑 裕 表 学研究費補 金 研究(C) 24 択 格子
子色力学 見 索 1,170千
6. 研究業績
(1) 研究論文
A) 査 付 論文
1. Y. Nakamura, K.-I. Ishikawa, Y. Kuramashi, T. Sakurai, H. Tadano, Modified block BiCGSTAB for lattice QCD, Comp. Phys. Commun. 183, No.1 (2012) 34-37
2. PACS-CS Collaboration: S. Aoki, K.-I. Ishikawa, N. Ishizuka, K. Kanaya, Y. Kuramashi, Y. Nakamura, Y. Namekawa, M. Okawa, Y. Taniguchi, A. Ukawa, N. Ukita, and T. Yoshie, 1+1+1 flavor QCD+QED simulation at the physical point, Phys. Rev. D 86, No.3 (2012) ref.034507, pp.1-6
28
Yoshie, Charmed baryons at the physical point in 2+1 flavor lattice QCD, arXiv:1301.4743
4. T. Yamazaki, K.-I. Ishikawa, Y. Kuramashi, and A. Ukawa, Helium nuclei, deuteron, and dineutron in 2+1 Flavor lattice QCD, Phys. Rev. D 86, No.7 (2012) ref.074514, pp.1-9
5. Sinya Aoki, Takumi Doi, Tetsuo Hatsuda, Yoichi Ikeda, Takashi Inoue, Noriyoshi Ishii, Keiko Murano, Hidekatsu Nemura, Kenji Sasaki (HAL QCD Collaboration), Lattice QCD approach to Nuclear Physics, Prog. Theor. Exp. Phys. 2012, No.1 (2012) ref.01A105, pp.1-38
6. Sinya Aoki, Janos Balog, Peter Weisz, Toward an understanding of short distance repulsions among baryons in QCD − NBS wave functions and
operator product expansion − , Prog. Theor. Phys. 128, No.6 (2012) 1269-1282 挾 挾
7. S. Aoki, J. Balog, T. Doi, T. Inoue, P.Weisz, Short Distance Repulsion Among Baryons, Int. J. Mod. Phys. E22, No.5 (2013) ref.1330012, pp.1-16
8. Sinya Aoki, Bruno Charron, Takumi Doi, Tetsuo Hatsuda, Takashi Inoue, Noriyoshi Ishii, Construction of energy-independent potentials above inelastic thresholds in quantum field theories, Phys. Rev. D87, No.3 (2013) ref.34512, pp.1-10
9. N. Ishii, S. Aoki, T. Doi, T. Hatsuda, Y. Ikeda, T. Inoue, K. Murano, H. Nemura and K. Sasaki (HAL QCD Collaboration), Hadron-Hadron Interactions from Imaginary-time Nambu-Bethe-Salpeter Wave Function on the Lattice, Phys. Lett. B712 (2012) 437-441
10. Shinji Ejiri, Kazuyuki Kanaya, and Takashi Umeda, for the WHOT-QCD Collaboration, Ab initio study of the thermodynamics of quantum
29
11. T. Umeda, S. Aoki, S. Ejiri, T. Hatsuda, K. Kanaya, Y. Maezawa, H. Ohno (WHOT-QCD Collaboration), Equation of state in 2 + 1 flavor QCD with improved Wilson quarks by the fixed scale approach, Phys. Rev. D 85, No.9 (2012) ref.094508, pp.1-11
12. Y. Maezawa, T. Umeda, S. Aoki, S. Ejiri, T. Hatsuda, K. Kanaya and H. Ohno (WHOT-QCD Collaboration), Application of fixed scale approach to static quark free energies in quenched and 2 + 1 flavor lattice QCD with improved Wilson quark action, Prog. Theor. Phys. 128, No.5 (2012) 955-970
13. S. Ejiri, S. Aoki, T. Hatsuda, K. Kanaya, Y. Nakagawa, H. Ohno, H. Saito, T. Umeda (WHOT-QCD Collaboration), Numerical study of QCD phase diagram at high temperature and density by a histogram method, Cent. Eur. J. Phys. 10, No.6 (2012) 1322-1325
14. Shinji Takeda, Yoshinobu Kuramashi, and Akira Ukawa, Phase of quark determinant in lattice QCD with finite chemical potential, Phys. Rev. D85 No.9 (2012) ref.096008, pp.1-12
15. S. Aoki, H. Fukaya and Y. Taniguchi, Chiral symmetry restoration, the eigenvalue density of the Dirac operator, and the axial U(1) anomaly at finite temperature, Phys. Rev. D 86, No.11 (2012) ref.114512, pp.1-18.
16. K.-I. Ishikawa, Y. Iwasaki, Yu Nakayama and T. Yoshie, Conformal theories
with an infrared cutoff, Phys. Rev. D87, No.7 (2013) ref.071503(R), pp.1-5
挾
B) 査 無 論文
1. Y. Namekawa for PACS-CS collaboration, Charmed baryon spectroscopy on the physical point in 2+1 flavor lattice QCD, PoS (LATTICE 2012) (2012) ref.139, pp.1-7
30
3. N. Ishii for HAL QCD Collaboration, 2+1 flavor QCD result of nuclear forces, Proceedings of the 20th International IUPAP Conference on Few-Body Problems in Physics, Few-Body Systems in press
4. N. Ishii for HAL QCD Collaboration, Baryon-baryon Interactions from Lattice QCD, PoS (CD12) ref.025, to be appeared
5. H. Nemura for HAL QCD Collaboration, Lambda N and Sigma N interactions from lattice QCD, Proceedings of the 20th International IUPAP Conference on Few-Body Problems in Physics, Few-Body Systems in press
6. H. Nemura for HAL QCD Collaboration, Lattice calculation of baryon-baryon interaction and few-body systems, Proceedings of the International Workshop on Strangeness Nuclear Physics, to be published in Genshikaku Kenkyu Suppl.
7. K. Sasaki for HAL QCD Collaboration, Quark mass dependence of hyperonic interactions from lattice QCD, PoS (LATTICE 2012) (2012) ref.157, pp.1-7
8. K. Sasaki for HAL QCD Collaboration, Coupled channel approach to S-wave hyperonic interactions from lattice QCD, Proceedings of the 20th International IUPAP Conference on Few-Body Problems in Physics, Few-Body Systems in press
9. T. Umeda, S. Aoki, S. Ejiri, T. Hatsuda, K. Kanaya, Y. Maezawa, H. Ohno
(WHOT-QCD Collaboration), Thermodynamics in 2 + 1 flavor QCD with improved Wilson quarks by the fixed scale approach, PoS (LATTICE 2012) (2012) ref.074, pp.1-7
10. S. Ejiri, Y. Nakagawa, S. Aoki, K. Kanaya, H. Saito, T. Hatsuda, H. Ohno, T. Umeda (WHOT-QCD collaboration), Probability distribution functions in the finite density lattice QCD, PoS (LATTICE 2012) (2012) ref.089, pp.1-7
31
(2) 国 会議 表
A) 招 講演
1. Sinya Aoki Chiral symmetry restoration and eigenvalue density of Dirac operator at finite temperature 招 講演 , Workshop ’New Horizons for
Lattice Computations with Chiral Fermions’ (BNL, Upton, NY, USA, May 14-16, 2012)
2. Yoshinobu Kuramashi 1+1+1 flavor QCD+QED simulation at the physical point 招 講演 , Workshop ’New Horizons for Lattice Computations with Chiral Fermions’ (BNL, Upton, NY, USA, May 14-16, 2012)
3. Y. Namekawa Charm quark physics from lattice QCD 招 講演 , International Workshop on heavy quark hadrons at J-PARC of Technology 2012 (Tokyo, Japan, Jun. 18-22, 2012)
4. Sinya Aoki Hadron interactions in lattice QCD 招 講演 , Lattice Hadron Physics IV (CSSM, University of Adelaide, Adelaide, Australia, Jul. 2-4, 2012)
5. Yoshinobu Kuramashi Lattice QCD – From quarks to nuclei – 招 講演 , 10th International Meeting on High-Performance Computing for
Computational Science (VECPAR2012) (Kobe University, Kobe, Jul. 17-20,
2012)
6. Sinya Aoki Hadron interactions from lattice QCD 招 講演 , INT
Program ’Lattice QCD studies of excited resonances and multi-hadron systems’ (INT, Seattle, USA, Jul. 30-Aug. 31, 2012)
7. N. Ishii for HAL QCD Collaboration, Baryon-baryon interaction from lattice QCD 招 講演 , The 7th International Workshop on Chiral Dynamics (CD12)
32
8. Sinya Aoki Chiral symmetry restoration, eigenvalue density of Dirac operator and axial U(1) anomaly at finite temperature 招 講演 , Workshop ’New
Frontiers in Lattice Gauge Theory’ (GGI, Florence, Italy, Aug. 27-Sep. 28, 2012)
9. K. Sasaki for HAL QCD Collaboration Strangeness S=-2 baryon-baryon
systems from Lattice QCD 招 講演 , Tours2012, VIII Tours symposium on Nuclear physics and Astrophysics on Black Forest (Lenzkirch-Saig, Germany, Sep. 2-7, 2012)
10. H. Nemura for HAL QCD Collaboration Lattice Calculation of Hyperon Potential and Few-Body System , The 2nd Korea-Japan Workshop on Nuclear and Hadron Physics at J-PARC (Busan, Seoul, Sep. 24-25, 2012)
11. Sinya Aoki Hadron interactions from lattice QCD 招 講演 ,
Workshop ’New Frontiers in Lattice Gauge Theory’ (GGI, Florence, Italy, Aug. 27-Sep. 28, 2012)
12. Sinya Aoki Computing Nuclear Force in Lattice Gauge Theory 招 講演 , QCD Stucture I (Central China Normal University, Wuhan, China, Oct. 7-20, 2012)
13. K. Kanaya Finite density QCD with Wilson quarks using the histogram method 招 講演 , International workshop ”QCD Structure I” (Wuhan, China,
Oct. 7-20, 2012)
14. Sinya Aoki Hadron interaction from lattice QCD 招 講演 , KMI workshop
on ”Strong Coupling Gauge Theories in the LHC Perspective”, (Sakata-Hirata Hall, Nagoya University, Nagoya, Japan, Dec. 4-7, 2012)
15. Sinya Aoki Quarks to Universe in computational science 招 講演 ,
33
16. Sinya Aoki Nuclear force in lattice QCD 招 講演 , KEK theory center
workshop on ’Hadron physics with high-momentum hadron beams at J-PARC in 2013’ (KEK Tsukuba Campus, Japan, Jan. 15-18, 2013)
17. K. Sasaki for HAL QCD Collaboration Baryon-baryon interactions in strangeness sector from lattice QCD 招 講演 , Workshop on Future
Prospects of Hadron Physics at J-PARC and Large Scale Computational Physics in 2013 (Tokai, Japan, Feb. 11-13, 2013)
18. K. Kanaya Finite density QCD on the lattice with a histogram method 招
講演 , International workshop on ”Quarks, Gluons, and Hadronic Matter under Extreme Conditions” (St. Goar, Germany, Mar. 18, 2013)
1重普
K. Kanaya Finite density QCD with a histogram method 招 講演 , School and workshop ”New Horizons in Lattice Field Theory” (Natal, Brazil, Mar. 25, 2013)挾
挾
B) 一般講演
1. Sinya Aoki Chiral symmetry restoration and eigenvalue density of Dirac operator at finite temperature , The 30th International Symposium on Lattice Field Theory (Cairns Center, Cairns, Australia, Jun. 24-29, 2012)
2. Y. Taniguchi Non-perturbative evaluation of cSW for smeared link clover
fermion with Iwasaki gauge action , The 30th International Symposium on
Lattice Field Theory (Lattice 2012) (Cairns Convention Centre, Cairns, Australia, Jun. 24-29, 2012)
3. N. Ishii for HAL QCD Collaboration 2+1 flavor QCD results of nuclear forces ,
34
4. Y. Namekawa Charmed baryon spectroscopy on the physical point in 2+1 flavor lattice QCD , The XXX International Symposium on Lattice Field
Theory (Cairns, Australia, Jun. 24-29, 2012)
5. K. Sasaki for HAL QCD Collaboration Quark mass dependence of hyperonic interactions from lattice QCD , The 30th International Symposium on Lattice
Field Thoery (Lattice 2012) (Cairns Convention Center, Cairns, Australia, Jun. 24-29, 2012)
6. N. Ishii for HAL QCD Collaboration 2+1 flavor QCD results of nuclear forces , The 20th International IUPAP Conference on Few-Body Problems in Physics (Fukuoka, Japan, Aug. 20-25, 2012)
7. H. Nemura for HAL QCD Collaboration Lambda N and Sigma N interactions from lattice QCD , The 20th International IUPAP Conference on Few-Body Problems in Physics (Fukuoka, Japan, Aug. 20-25, 2012)
8. K. Sasaki for HAL QCD Collaboration Coupled channel approach to S-wave hyperonic interactions from lattice QCD , The 20th International IUPAP Conference on Few-Body Problems in Physics (Fukuoka, JAPAN, Aug. 20-25, 2012)
9. H. Nemura for HAL QCD Collaboration Lattice calculation of baryon-baryon interaction and few-body systems , The International Workshop on
Strangeness Nuclear Physics (Osaka, Japan, Aug. 27-29, 2012)
10. K. Sasaki for HAL QCD Collaboration Coupled channel approach to hyperonic interactions from lattice QCD , HYP2012-XI International
Conference on Hypernuclear and Strange Particle Physics (Barcelona, Spain, Oct. 1-5, 2012)
11. Y. Namekawa Charmed baryons from lattice QCD , Hadron Structure and
35
12. K. Sasaki for HAL QCD Collaboration Coupled channel approach to S=-2 baryon-baryon system in Lattice QCD , International Symposium on ”Quarks
to Universe in Computational Science (QUCS 2012)” (Nara New City Hall, Nara, Japan, Dec. 13-16, 2012)
13. M. Yamada Omega-Omega interaction on the Lattice , Symposium:Quarks to
Universe in Computational Science (QUCS 2012) (Nara New City Hall, Nara, Japan, Dec.13-16, 2012)
14. H. Nemura for HAL QCD Collaboration Hyperonic nuclear forces from lattice QCD and toward an application to few-body systems , YITP workshop Nuclear equation of state and hypernuclear physics (Kyoto, Japan, Jan. 4-31, 2013)
15. Y. Namekawa Charm quark physics from lattice QCD at the physical point , Future Prospects of Hadron Physics at J-PARC and Large Scale
Computational Physics (Ibaraki, Japan, Feb. 11-13, 2013)
16. H. Nemura for HAL QCD Collaboration Hyperonic potentials from lattice QCD and toward an application to few-body problems , The 5th GCOE International Symposium on ”Weaving Science Web beyond Particle-Matter Hierarchy (Sendai, Japan, Mar. 4-6, 2013)
(3) 国 学会 研究会 表
A) 招 講演
1. 青木慎也 招 講演 , 新学術領域 素 融 搬 新
研究会 屋大学, 屋, Jul. 12-13, 2012
2. 滑 裕 Charm quark system on the physical point in 2+1 flavor lattice QCD
招 講演 , 新学術領域 素 融 搬 新 研究
36
3. 金谷和 Lattice QCD at finite T and μ – updates from Lattice 2012 (招
講演), 研研究会 熱場 子論 応用2012 (京都大学 礎物理学研究所, 京
都, Aug. 22-24, 2012)
4. 藏増嘉伸 京 HPCI戦略 5 招 講演 , 斌 物理
学会2012 秋季大会理論 物理領域-素粒子論領域- 線 物理領域
動 出 京 素粒子 原子 計算物理 (京都
産業大学, 京都, Sep. 11-14, 2012)
5. 青木慎也 HAL QCD 招 講演 , 第25回理論懇 計算 物理
学 新展開 ( 国 会議場中 300, , Dec. 22-24, 2012)
6. 藏増嘉伸 計算素粒子物理学 招 講演 , 斌 物理学会第68
回 大会素粒子論領域-理論 物理領域-領域3-9-領域11-12
向 出 計算物理学 ( 島大学東 島 , 東 島, Mar. 26-29, 2013)
B) 表
1. 石 人 Resonance on the lattice , 新学術領域 素 融 搬 新
研究会 屋大学, 屋, Jul. 12-13, 2012
2. 々木健志 for HAL QCD collaboration Coupled Channel Approach to
Baryon-Baryon Interactions in Lattice QCD , 新学術領域 素 融 搬
新 研究会 屋大学, 屋, Jul. 12-13, 2012
3. 石 理 for HAL QCD Collaboration PACS-CS 配 用い 2+1 flavor
QCD 力 , 斌 物理学会2012 秋季大会 京都産業大学, 京
都, Sep. 11-14, 2012
4. 滑 裕 物理 2 1 格子 QCD
37
5. 々木健志 for HAL QCD collaboration 格子QCD 間相 作用
質 依 性 研究 , 斌 物理学会2012 秋季大会 京都産業大学, 京 都, Sep. 11-14, 2012
6. 照 HEPnet-J/sc報告 , HEPnet-J 会 九 大学, 福岡, Feb. 20-21,
2013
7. 照 JLDG 現状 計 , HPCI戦略 5 物質 起源
構造 全体 秋葉原, 東京, Mar. 5-6, 2013
8. 々木健志 for HAL QCD collaboration 格子QCD S=−2 間相
作用 質 依 性 研究 , HPCI戦略 5 物質 起
源 構造 秋葉原, 東京, Mar. 5-6, 2013
9. 照 HPCI共用 JLDG連携 , HPCI 利用研究課 中間
報告会 , 東京, Mar. 15-16, 2013
10. 齋藤華 格子QCD 関数 計算 関 研究 , 斌
物理学会第68回 大会 島大学東 島 , 東 島, Mar. 26-29, 2013
(4) 著書 解 等
1普
藏増嘉伸
京
子論
数値計算
「01」
1
丸善出爮挾
「普
青木慎也
石
理
初
哲男
英克
６（）挾
力研究
展開
斌
物理学会
「01「
盤】
巻
11
挾
」普
青木慎也
「01「挾
賞
石山智明氏
似鳥啓吾氏
牧
淳一郎
氏
斌
物理学会
「01」
盤里
巻
」
挾
『普
山
信博
学
芽
賞実行委員会
岩崎洋一
貫挾
白
秀樹
貫挾
金谷和
挾
著
知
い！
学
芽
世界
桑Ｆ三陥挾
」
筑波大学出爮会
「01「
挾
7. 異 間連携 国 連携 国 活動等
38 http://www.jicfus.jp/jp/
2. High Performance Computing Infrastructure (HPCI) 戦略
http://www.jicfus.jp/field5/jp/
3. International Lattice Data Grid (ILDG)
http://ildg.sasr.edu.au/Plone 4. Japan Lattice Data Grid (JLDG)
http://www.jldg.org/jldg/, http://ws.jldg.org/QCDArchive/index.jsp
8. 研究会 等 開催実績
1. 青木慎也 新学術領域 素 融 搬 新 HPCI戦略
5共催 研究会 −多彩 新 い
形態 包括的件研究− 共催 2012 7 12 13斌 屋大学坂 屋 参 者75
2. 青木慎也 新学術領域 素 融 計算 学 い 層的物質構
造 解明 /HPCI戦略 5 物質 起源 構造 /京都大学 礎
物理学研究所主催 超新斴爆 − 礎物理 理
想 挑戦− 2012 7 27斌-31斌 京都大学 礎物理学研究所 京都 参 者61
3. 青木慎也 新学術領域 素 融 計算 学 い 層的物質構
造 解明 /HPCI戦略 5 物質 起源 構造
Symposium: ’Quarks to Universe in Computational Sciences (QUCS 2012)’, 主
催 2012 12 13斌-16斌 良新公会堂 良 参 者101
9. 管理 遀営
1. 青木慎也 遀営委員会委員 遀営 議会委員
2. 照 共 研究遀用委員会委員
3. 照 藏増嘉伸 計算機 遀用委員会委員
40
II
原子
物理研究部門
II-1.
1.
教授挾 挾 挾 梅 之 准教授挾 挾 森挾 講師挾 挾 挾 耕 准教授 岡 崇 HPCI戦略
教挾 挾 挾 勝 望 理数学生応援 研究員挾 挾 谷 衝 研費
行方 大輔 研費 俊宏 研費 石山 智明 HPCI戦略 長谷 賢 HPCI戦略 Alexander Wagner
学生 大学院生挾 挾 学類生挾 挾
2. 概要
当 FIRST 推 逭 共 FIRST,
T2K-Tsukuba 用い 論的銀河形 粒子法輻射流体
再電 流体用 速輻射輸 計算 大質 多体
系 体 長 ― 部構造 ― 問
涙 銀河系中心領域 特異 物理環境 子雲形 え
響 6 元Vlasov方程式 自己 力系数値 研究 行
AVX 爮 Phantom-GRAPE 開 HA-PACS 向 開 京
用い 銀河形 行 生 物性 原子 連携
斴間空間 誘起 L 型 酸遃 系 惑斴 機構 理
論的研究 行 戦略 銀河輻射流体
力学 幕開 表者挾 森 択 輻射流体力学実現 演算 速
器 載 並列計算機 構築 開始 制作 行
3. 研究
挾一般相 論的N体計算 巨大 体遃程
巨大 銀河中心 観測 質 銀河
41
い Kormendy & Richstone 1995; Magorrian et al. 1998; Merrifield et al.
2000; Merritt & Ferrarese 2001, Marconi & Hunt 2003 巨大
形 銀河 形 密接 関係 質 決
遂的 物理 示唆 層的 体形 論
え 大 銀河 銀河 体 生 大 銀河
多数 巨大 銀河 中心
‐ 質 関係 満 巨大 一 あ いう観測
実 明 い 能性 一 銀河 体
全 体 銀河中心 落 いう あ 巨大
体 極 い 々 持
銀河 体 後 体 起 いう問 一般相 論効
入 精 N体計算 調 銀河 中 10個 巨大
置い 計算 結 銀河中 斴密 い場 斴 力学的
摩擦 効 働 銀河中心 落 多 散乱 繰 返
連 生 連 散
乱 結 力波 攡出 体
体 う 銀河 起 う いう 解析 結 望遠鏡
見え い 赤方偏移 z=7
一般相対論的Ｎ体計算
ラックホー
合体条件
観測
比較
ー ン っ 得 た
ラックホー 合体条件
線 付け た数値 仮定した質量－光度比
縦 波線 左側 合体 能
Tanikawa & Umemura 2012
z=7 ライ ンα ッター 光度関数
42 体遃程 起 い 能性 示
挾 第一世 体 輻射流体計算
初 生 体(第一世 体) 来 計算 106M
105M
程
蓄 積 形
来 研究 振 舞い
特 長 い
あ 注目
々 計算領域 全体 能 限 解像 数値 計算
長 解 う 第一世 体 形 逭 調
計算 流体粒子 1 億体(5123体)使い 計算
領 域 全 体 わ 質 解 像
0.046M 0.22 M 実現 結 第一世 体 質 来
計算 求 い 質 2 近 103M
程 わ
結 一部 既 特 推逭研究 得 い あ
後 解析 力 温 昇
結 水素 子冷却 熱的 引 起 質 体形
能 明 Umemura et al. 2012
挾 原始銀河形 再電
赤方偏移z=1000程 一 中性 後形 体 輻射
再び電 考え い 線系 利用 観測(e.g., Fan et al. 2006)
z=6程 電 い 示 い 再電
遃程 い 始 う 逭 明 い 再電 解 明
銀河間空間 電 子供給源 あ 体 形 伴う電 領域 逭 遃程 旓 計算 必要 あ 体形 自体 輻射性
響 輻射流体計算 必要 う 計算 多 攡射源 扱
う必要 あ 非常 膨大 計算 要求
々 犉自 開 構造 用い 速 輻射輸 計算 行え
START (SPH with Tree-based Accelerated Radiative Transfer) (Hasegawa &
43
電 行 結 斴形 主 熱効 著
阻害 示 (Hasegawa & Semelin 2013) 再電 遃程 要
鍵 電 子 銀河 脱出割 紫 線 密 領域
昇 (Umemura et al. 2012)
京 う 大規模並列計算機 効率 出 うSTART 改善 行 結 数1000並列 計算 あ 非常 い
実 現 能 演 算 数 自 体 半 (START2: Hasegawa in prep.)
挾 京 用い
理 学研究所 京 力多体 用 "GreeM"
々 数万や数十万並列 う 適 行
1) 大規模 全 通信 層 全 い 割 旓
空間的 近い領域 担当 う
全 通信 行 間 通信 いう う 2 段 通信 行 う 全 通信 実行旓間 抑制
2) 領域形状 遃去 数 渡 均
間 形状 抑え 粒子再配 通信 抑制
3) MPI+OpenMP 並列
う 京 璧 実現 2
粒子 力逭 京 全 用い
5.67 実効性能(実行効率55%) 実現
関 国 会議SC12 2012 11 米国
開催 い 賞 犉 賞
賞 性能 20 ( 京 2倍) ( 米
国 国立研究所 用い 様
14 遉 米国 あ 々 実
計算速 回 1 粒子あ 2.4 倍 速 能 あ
( 計算機 用い 場 5 倍近 速い) うい 評価 賞
賞 う あ
挾 活動銀河 組 研究
44
電子 陽電子 あ いう 的 問 明 い い 相
論的 形 や 線 起源 関連 物理学 い 要
あ 的攡射 無
自己 や 偏波 中 組 い 議
論 組 密接 関係 真 ワ 評価
困 あ 来 手法 攡射 類推 問 解
決 々 貯蔵庫 あ 膨張 力学 注目
真 ワ 齢 評価 方法 (Kino and Kawakatu 2005; Ito et al. 2008)
各 電子 陽電子 陽子 相 論的 組
制限 方法 唱 (Kino, Kawakatu, Takahara 2012) 研究 膨張
力学 真 ワ 求 い FRII型電波銀河Cygnus A, 3C219,
3C223, 3C284 様 方法 適応 解析 行 結 Cygnus A 含
FRII電波銀河 敂 ワ 持 関わ
い 場 混 4 体
組 電 子 陽 子 明 陽 電 子 示 唆 あ
(Kawakatu, Kino, Takahara in preparation)
挾 ― 部構造 ― 問
構造形 (CDM) N体
(DM) 中心部 質 密 散
( )構造 予言 (Navarro, Frenk & White 1997; Fukushige & Makino
1997) 観測結 銀河 密 中心部 散
一 ( ) 報告 い (Swaters et al. 2003; Spekkens et al.
2005; Oh et al. 2010) 理論 観測 一 '' — 問 ''
CDM 解決問 一 あ 回 々 ''超新斴爆
銀河 熱 膨張 や 攡射冷却 縮 起 再び斴形 起
いう一連 遃程 繰 返 力場変動 遤移
'' N体 用い 調 力場 周期的
旓間変 場 表現 結 力場変動 旓間
依 DM 形 大 や 置 大 変
わ DM 粒子群 場間 起 共鳴的 現象 解
析的 構築 共鳴 半 力場 旓間変
動周期 関係 あ 突
45
近 望遠鏡や 望遠鏡 表 地 大型望遠鏡 大
限活用 近傍 深 査 現 銀河逭 遃程 垣間見
う 銀河周 い び い数
い 銀河 見 銀河 衝突 思わ
や 等 痕跡 々 明
い 特 涙 関 観測 理
論 両面 研究 逭展 銀河衝突 軌遈遀動や 旓期 衝突 銀河 質 や 学組 等 い 理解 逭 い 々 N 体 計算 3 元 流体力学計算 組 わ
銀河 銀河 付 流体力学的 相
作用 細 い 調 銀河 銀河 相 作用 生
銀河 流体力学的 挙動 観測
い 状構造 生 遃程 い 細 解析 行 銀河 軌
遈遀動 初期条件依 性 い 大規模 行い 観測 再現 軌遈要素 求
挾 銀河 中 漂う巨大 降着 波長域攡射
各銀河 中心 巨大 宿 銀河 周 銀
河 衝突 体 併 質 長 考え
銀河衝突 衛斴銀河 壊 後 間 元衛斴銀
河 中心巨大 親銀河 中 漂 い 期 隣
銀河 あ 銀河 銀河衝突 銀河
共逭 理論 観測 両面 理解 実験場 あ 細 観測 数値実験 比較 10億 斮 衛斴銀河 衝突 現 元衛斴銀河
破爭 銀河周 痕跡 い わ い
元衛斴銀河 中心巨大 希薄 Bondi 降着
い 旓 形 降着 波長域攡射 計算 様々 検出装置 感 比較 電波領域 既 装置 十 検出 能 あ わ
46
冷却や斴形 , 超新斴爆 等 物理遃程 入 論的銀河形
用い う 研究 行 (i) 遠方 斴形 銀河
あ 波 明 輝 銀 河 物 理 的 性 質 明
(Shimizu et al. 2012ab). (ii) 銀河系 銀河 解能 論的
形 , 銀 河 系 形 い 調 . 擬 構
造 考え い う 逭 赤方偏移 爆 的斴
形 形 いう新 い 像 案 (Okamoto 2013). , 様々
間 銀 河 形 比 較 参
(Scannapieco et al. 2012).
挾 6 元 相空間 Vlasov-Poisson
無衝突自己 力系 数値 N体 幅
用い 速 散 大 い 無衝突浿衰 確 扱
う 向 あ N体 弱 克服 手法 6 元
相空間 無衝突 方程式 限体積法 直接数値
手法 世界 初 開 論的 大規模構造形 数値
応用 う 共動 標系 無衝突 方程式 数値解法 開
挾GPU 利用 輻射輸 開
輻射輸 物理学 数値 い 要 役割
多数 源 輻射輸 計算 計算 膨大 大規模 計算や流体 力学 組 わ 輻射流体計算 実行 困 あ 々 多数 源
輻射輸 速 計算 ARGOT (accelerated radiation transfer on grids
using oct tree) 法 GPU 用い 速 GPU
載 大規模 あ HA-PACS 並列 行
挾 AGN 輻射 子雲 逭 い 研究
活動銀河 (AGN) 明 い 体 1 あ 銀河 体 敂
大部 電 能力 持 紫 び X線 波長 攡射
銀河 逭 大 響 え 一般 考え い AGN 活動性 細
う 決 明 い 研究 AGN現象 持 性 理
解 逭 目的 供給 要 役割 期 学的 厚い
47
子雲 逭 2 種類 大 わ 電 弱い場
子雲 球 称的 蒸 質 失い 作用 力
子雲 一部 力崩壊 電 強い場 輻射 効 的
蒸 流 子雲表面 閉 子雲 ぶ 終的 子雲
一部 力崩壊 わ
挾
極 X線源 中間質 候補
中 間 質 (BH) 力 候 補 あ 極 X 線 源(Hyper
Luminous X-ray source) い X線 降着 用い 解
析 BH質 推
極 X線源 ESO243-49 HLX-1 近傍 銀河(赤方偏移z = 0.02)
面 空 写 赤方偏移測 試 銀河 距 居
示唆 い 距 測 X線 1042 erg/s 遉 超
X線源 (Ultra Luminous X-ray source) 1 明 中間質
BH 力候補 あ 中間質 BH 実 起源 何 元衛斴銀
河 中心BH 浮逿 い 巨大BH 種 等 疑問 挑 絶好 実 験場 言え
HLX-1 X線 降着 用い 解析 研究
亜臨界降着率 限 標準降着 い 大
体 亜臨界降着 仮 BH 質 推 中間質 いう答え 得 仮 結論 直結 い 険性 あ
亜臨界 超臨界降着率 幅 降着率 前 仮
必要 移流 相 論的効 表層 散乱 含
(Kawaguchi 2003) 10万 陽質 BH 計算
3 渡 30倍 旓間変動 HLX-1 多旓 X線 解析
1.8 万 陽質 BH 降着 明
わ
挾 AVX 爮 Phantom-GRAPE
2011 Intel 新型 Sandy-Bridge 新 実装
SIMD (Single-Instruction-Multi-Data) あ Advanced Vector
eXtension (AVX) 用 い 力 多 体 数 値 計 算
Phantom-GRAPE 開 SIMD あ Streaming SIMD
48
Phantom-GRAPE 4 積 軌遈計算 用い 精 衝突系
力多体用 銀河や銀河団 無衝突系用 2種類 開 Google Code 公開 あ
挾 斴間空間 酸 偏 波誘起鏡像異性体遃 理論的研究 生 物性 原子 連携
生 体 子 酸 あ 実験室 酸 作 巻
L型 巻 D型 生 地球 生 L型
酸 使わ い い 鏡像異性体遃 いう 1969
石 落 石 酸 検出
わ あ 鏡像異性体遃 見 2010 超 温
石 酸 見 石 酸 地球 大気圏通遃 変
落 見 鏡像異性体遃 わ あ
実験 鏡像異性体遃 自己触媒 応 急速 増大
酸 鏡像異性体遃 空間 起 石 通 地
球 遀 地 急速 増幅 能性 あ 実験室 偏
当 鏡像異性体遃 引 起 近
斴形 領域(OMC-1) 偏 波 見 原
始系 近 大質 斴 誕生 陽系 酸 鏡像異性体遃
起 能性 あ 実 背 空間 偏 波 酸
鏡像異性体遃 引 起 遃程 い 子多体計算 逭 い 々 偏 波 酸 励起 崩壊 改変 応 遃程 解析 第一原
理計算 真空中 酸 構造 求 性 偏
色性 値 求 結 酸 物性 種 特 側鎖
酸全体 共通 主鎖 構造 強 依 波長 α
要 あ 明
挾 系 惑斴 検出 目指 機構
移動計算 生 連携
査機Kepler 陽系 惑斴 数 大 追 , 地球型惑斴や
入 惑斴 観測 い 地球 類似 惑斴 見 期
中 , 得 生 痕跡 あ 如何 検出
49
中 750nm付近 近赤 領域 見 射 特 的 勾配 red edge
力 得 様々 型 主斴 公転 系 惑斴 い
生物 想 場 , 捕 形態 大 異 考え 地球
生物 場 主斴 あ 陽 比 効率良 捕 う 逭
来 考え あ , 色素 系 環境
色素 種類や配置 形態 異 色素 電子励起 , 近 色素
電子状態 移動 いう遃程 通 効率的 学 変換
系 惑斴 い 異 捕 形態 場 , 応 波長域
候 示 能性 あ 々 系 惑斴 比 植物 捕 形態
相関 い 的指標 出 目的 捕 系 子 学計算 行 構 色素1 1 励起 状態 旓間依 密 汎関数法 計算 , 色素 間 極子- 極子相 作用 近似 構築 , 入射 あ 振動数
電場 系 子力学的旓間 展 追跡 強 算出
系 形態 red edge 相関 調 , 系 惑斴 想
議論
4. 教育
学 論文 博士論文
1. Mohammad Ali Nawaz
Interaction of Jets with the Intracluster Medium
2. Marc White
Outflows from Young Stellar Objects
士論文 部挾 牧人
3 元輻射流体計算 紫 線輻射場中 斴団形 遃程 研究
十嵐挾 朱夏
球 称 常銀河風 解析 大 挾 純
挾 GPU 用い 輻射流体力学 速 久保 挾 明夏
50 学士論文
飯 挾 実
系 惑斴系 軌遈交差 形
中講義
梅 之
輻射輸 輻射流体力学 2012 7 19斌 千葉大学 千葉 岡 崇
銀河形 第7回 素 融 (2012 11 20
斌〜2012 11 21斌, 神戸大学, 神戸 )
5. 賞 部資金 知的 産 等
賞
ACM Gordon Bell Prize 石 山 智 明 4.45 Pflops Astrophysical N-Body
Simulation
on K computer – The Gravitational Trillion-Body Problem, 2012 11 16
斌
HPCS IEEE Computer Society Japan Chapter 優秀若手研究賞 扇谷豪 力
多体
系用Tree Code 並列GPU 計算 速 2013 1 16斌
部資金 表者
研究 梅 之 表者 2008 択
挾 第一世 体 原始銀河 黒旓 解明 (880万 )
研究 森 表者
挾 理論 観測 融 銀河 生学 求 (560万 )
戦略 森 表者
挾 銀河輻射流体力学 幕開 (500万 ) 研費若手研究 :岡 崇 ( 表者) (新規)
輻射流体 用い 銀河風駆動機構 解明 (110万 ) 研費若手研究 :長谷 挾 賢 表者 新規 (210万 )
論的輻射 解明 再電 期 銀河形
研費若手研究 :石山智明 表者 2013-2015 350万(直接経費)
51 担者
研究 梅 之 担者 表者 大 巳 2011 択 世 大規模 査 挑 再電 (1万 )
研究 森 ( 担者) 表者 大 巳 新規
世 大規模 査 挑 再電 (2.5万 )
研究 俊宏 ( 担者)挾 表者 峰崎 岳
活動銀河 多波長 活動銀河 変 攡射機構 研究
(20
万 )
挑戦的萌芽研究: 勝挾 望 ( 担者) 表者 長尾
生 石 現 潜 逭 初期 銀河 巨大 査 15万
6. 研究業績
研究論文
A) 査 付 論文
1) Doi, A., Nagira, H., Kawakatu, N. Kino, M., Nagai, H., Asada, K., 2012, Radio Galaxies
in Narrow-line Seyfert 1 Galaxies, The Astrophysical Journal, 760(1), 41-51
2) Godet, B. Plazolles, T. Kawaguchi, J.-P. Lasota, D. Barret, S. Farrell, V. Braito, M.
Servillat, N. Webb and N. Gehrels, 2012, Investigating slim disk solutions for HLX-1 in
ESO 243-49, The Astronomical Journal, 752, 34 [12pages]
3) Hasegawa, K., Semelin, B., 2013, The impacts of ultraviolet radiation feedback on
galaxies during the epoch of reionization, Monthly Notices of the Royal Astronomical
Society, Volume 428, pp.154-166
4) Ishiyama, T., Nitadori, K., Makino, J. 2012, 4.45 Pflops Astrophysical N-Body
Simulation on K computer – The Gravitational Trillion-Body Problem SC '12
Proceedings of the International Conference on High Performance Computing,
Networking, Storage and Analysis, Article No. 5
5) Kawata,D., Okamoto, T., Gibson, B. K., Barnes, D. J. and Cen, R., 2013, Calibrating an
updated smoothed particle hydrodynamics scheme within gcd+, Monthly Notices of the
Royal Astronomical Society, 428, 1968-1979
6) Kino, M., Kawakatu,N., Takahara, F., 2012, Calorimetry of AGN Jets: Testing Plasma
