筑波大学 計算科学研究センター 令和二年度 年次報告書

3 各研究部門の報告
I. 素粒子物理研究部門
1.メ
メンバー
教授

藏増

嘉伸、金谷

和至（共同研究員、数理物質系特命教授）

准教授

石塚

成人、谷口

裕介、吉江

友照、山﨑

助教

大野

浩史

研究員

浮田

尚哉、新谷

栄悟、吉村

友佑

学生

大学院生

剛（共同研究員、数理物質系）

8 名、学類生 2 名

概要
2.概
当部門では、数理物質系との密接な連携のもと、格子 QCD の大型シミュレーション研究を推進している。2016 年秋から JCAHPC(最先端共同 HPC 基盤施設:筑波大学と東京大学両機関の教職員が中心となり設計するスーパー コンピュータシステムを設置し、最先端の大規模高性能計算基盤を構築・運営するための組織)において Oakforest-PACS(略称「OFP」:ピーク演算性能 25PFLOPS の超並列クラスタ計算機、HPC(High Performance Computing)向けとしては「富岳」に次ぐ性能を有する)が稼働を開始した。昨年度に引き続き、本年度も筑波大学を中心とした PACS Collaboration に基づく共同研究体制のもと、OFP を用いた大型プロジェクト研究を推進した。これと並行して、有限温度・有限密度 QCD の研究、 テンソルネットワーク(TN)
形式に基づく格子ゲージ理論・スピンモデルの研究、標準理論を超える物理の探求など、活
発な研究活動を行った。さらに、格子 QCD 配位やその他のデータを共有する為のデータグリ
ッド ILDG/JLDG の構築・整備を推進した。
国内の計算科学全体の動向として、2015 年度で終了した HPCI 戦略プログラムの後継として、「ポスト「京」で重点的に取り組むべき社会的・科学的課題」に関するアプリケーション
開発・研究開発が始まった。9 つの重点課題が設定され、9 番目の課題である「宇宙の基本法
則と進化の解明」が素粒子物理・原子核物理・宇宙物理分野が対象とする基礎科学的研究課
題である。その活動は、http://www.jicfus.jp/jp に詳しい。また、重点課題と並行して、2016 年
度から 4 つの萌芽的課題が設定され、1 番目の課題である「基礎科学のフロンティア − 極限への挑戦」は基礎科学における分野横断的な研究課題であり、本グループも分担機関として参加している。重点課題・萌芽的課題ともに 2019 年度末で終了した。ポスト「京」の名称は正式に「富岳」と決定し、2020 年度から重点課題・萌芽的課題の成果創出フェーズである成果創出加速プログラムが 3 年間実施される。また、それと並行して 2020 年度は「富岳」の早期利用課題の募集が行われ、2021 年度以降には一般利用が開始される予定である。

- 14 ― 14 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

なお、2020 年度は新型コロナウィルス COVID-19 の世界的な感染拡大により、ほぼすべての国際会議や国内学会・研究会がオンライン開催あるいは中止となった。また、日々の研究活動においてもオンラインでの議論・会議が強く推奨されるなど、研究遂行に大きな制約が
生じた。
3.研
研究成果
【1】

PACS Collaboratio n による Oakfore st-PACS を用いた大規模シミュレーション（藏増、石塚、谷口、山崎、吉江、浮田、新谷）
2016 年秋に JCAHPC において Oakforest-PACS(OFP)が導入され、稼働を開始した。OFP はピーク演算性能が 25PFlops であり、HPC 向けとしては「富岳」に次ぐ性能を有する。昨年度に引き続き、本年度も PACS Collaboration に基づく共同研究体制のもと、OFP を用いて物理点における 2+1 フレーバーQCD の大規模シミュレーションを推進した。また、「富岳」の早期利用課題制度を利用して「富岳」向けのプログラムの最適化を行い、「富岳」を用いた本格計算の準備が整った。
過去 30 年以上にわたり、格子 QCD は主にハドロン単体の諸性質解明を目指して来た。現在の世界的な状況においては、2 つの大きな問題点が存在する。まず、物理点直上でのシミュレーションが可能になったことは事実だが、実際には物理点のみで物理量の評価を行えるほどの精度を得るレベルには至っていない。次に、現在の格子 QCD シミュレーションにおける物理量計算は“テーラーメイド”であると評されている。これは、目的とする物理量計算に応じて、適当と思われる物理パラメータ(クォーク質量や空間体積など)を選んでシミュレーシ
ョンすることを意味している。この場合、例えば、同じゲージ配位を用いた計算であっても、ある物理量に対しては良く実験値と合うが、他の物理量に関しては実験値を再現しないということが起こりうる。OFP を用いたプロジェクトでは、複数の格子間隔において物理点直上で(10fm)3 超の大空間体積を持つシミュレーションを行うことによって、上記 2 つの課題を克服
した計算を実現する。
ゲージ配位は異なる格子間隔 3 点(a=0.085fm,0.065fm,0.045fm)において生成し、系統誤差となる格子間隔依存性を取り除くために連続極限(a→0)を取る。2019 年度までに(格子間隔,
格子サイズ)=(0.085fm,1284)、(0.065fm,1604)のゲージ配位生成は完了しており、現在は(格子間隔,格子サイズ)=(0.045fm,2564)での物理点シミュレーションに注力している。また、並行してこれまで生成した格子サイズ 1284 と 1604 のゲージ配位を用いて、以下に説明するような物理量計算を行っている。

【2】

現実的クォーク質量を用いた K 中間子セミレプトニック崩壊形状因子計算（藏増、石塚、谷口、山崎、吉江、浮田、新谷）- 15 ― 15 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

賀数研究員と山崎は現実クォーク質量での K 中間子セミレプトニック崩壊の動的 2+1 フレーバー大規模格子 QCD シミュレーションを行うことで、CKM 行列要素の一つである Vus の決定を行った。この物理量はクォークの世代間混合を表す行列である CKM 行列の行列要素のうち、アップクォークとストレンジクォークの混合の度合いを表す量である。この行列は標準理論においてユニタリー性を持つので、ユニタリー性の確認を行うことで標準理論を超える物理の検証を行うことができる。
図 1 に示すように、本研究で得られた|Vus|(赤丸)は、これまでの多くの計算結果(Kl3
Nf=2+1+1, Nf=2+1)と異なり、標準理論から予測される値(灰色帯)と無矛盾な結果になった(研究論文 B-1)。2020 年度には、この計算よりも格子間隔を小さくした計算を実行し、以前の結果と一致する中間結果(赤四角)を得た。この結果には一部の系統誤差の見積りが含まれていないため、誤差が以前の結果よりも小さくなっている。
これらの結果は K 中間子レプトニック崩壊から決定される|Vus|(青丸、緑星)とも一致しており、標準模型を超える物理のシグナルに対しては否定的な結論を示唆している。しかし、近年の標準模型の予測では図中の灰色帯を超えた大きな|Vus|が報告(水色帯)されており、標準理論を超える物理の検証を行うために、さらに精密な計算結果を得る必要がある。今後は、有限格子間隔起因の系統誤差を取り除くために、異なる格子間隔で得られた結果を用いて連続極限への外挿を実行する計画である。また、さらに小さな格子間隔のゲージ配位が今年度中に生成される予定であるため、そのゲージ配位を用いた計算も予定している。

図 1: |Vus|の比較。灰色帯と水色帯は CKM 行列のユニタリー性から求まる標準理論の予測。

- 16 ― 16 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書【3】

格子 QCD を用いた核子構造研究（藏増、山崎、新谷）陽子と中性子(核子)はクォークの束縛状態であり、その構造を詳細に調べるためには、強い相互作用の第一原理計算である格子 QCD を用いた計算が必要である。これまでに格子 QCDを用いて、核子構造に関係する核子形状因子研究が行われてきたが、非常に良い精度で測定されている実験値を再現できていなかった。
藏増、山崎は、広島大学石川健一准教授、東北大学佐々木勝一准教授、理研計算科学研究センター青木保道チームリーダー、新谷栄悟研究員、東北大学大学院生塚本夏基氏、辻竜太
朗氏とともに、PACS Collaboration において、現実的クォーク質量直上で核子形状因子計算を行った。2020 年度は格子間隔 a=0.065fm での形状因子計算を進めるとともに、核子スカラ
ー電荷と核子テンソル電荷(研究論文 B-2)に加え、核子構造関数に関係する物理量に対する繰り込み計算を実行した。
また、現実的クォーク質量直上での形状因子計算結果(研究論文 A-2)で明らかになった荷電半径計算に含まれる系統誤差の原因を探るため、異なる格子間隔を用いた計算による有限格子間隔系統誤差の調査を行っている。また、これまでとは異なる観点から計算を実行する
ため、新しい荷電半径計算方法の開発も行っている。

【4】

格子 QCD によるクォークを自由度とした原子核の直接構成（藏増、山崎）
藏増、山崎は宇川名誉教授との共同研究により、2010 年世界で初めて格子 QCD によるヘリ
ウム原子核の構成に成功し、そののち 2 核子系の束縛状態である 重陽子の構成にも成功した。これらの計算は、計算コストを抑えるためにクェンチ近似かつ重いクォーク質量を用いた試験的なものであった。その後、広島大学石川健一准教授を共同研究者に加え、真空偏極効果を取り入れた 2+1 フレーバーQCD シミュレーションを行い、試験的計算より現実世界に近い状況でのヘリウム原子核および 2 核子系の束縛エネルギー計算に成功した。この計算はπ中間子質量 0.5GeV と 0.3GeV のクォーク質量を用いたものであり、物理点(π中間子質量0.14GeV に相当)よりも重い質量を用いていた。この成果を踏まえ、京コンピュータで生成された 964 格子サイズのゲージ配位を用いた現実に近いπ中間子質量 0.146GeV での軽原子核束
縛エネルギー計算を進めている。この計算は統計誤差を抑えることが非常に難しく、ヘリウ
ム原子核については有意な結果は得られていないが、重陽子については現状で実験値から予測された値を再現する結果が見え始めている。
一方、現実よりも重いπ中間子を用いた他グループの計算からは束縛状態を観測できないという問題が報告されている。この計算には 2 体核子散乱を考慮した計算になってはいるが、束縛エネルギー計算で問題となるのは核子励起状態散乱の寄与と考えられる。今後、この状態を考慮した計算方法を開発し、束縛エネルギーに含まれる、それら状態の系統誤差を見積もる予定である。
- 17 ― 17 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書【5】

格子 QCD による相互作用範囲内の Bethe-Salpeter 波動関数を用いた散乱振幅（山崎）
散乱長は、粒子間の相互作用を特徴づける基本的な物理量である。これまで散乱長は、相互作用範囲外の Bethe-Salpeter 波動関数から導出された有限体積法により計算されてきた。
京都大学滑川特任助教と山崎は 2 パイ中間子系において、相互作用範囲内の Bethe-Salpeter
波動関数を用いて散乱振幅を決定した。得られた散乱振幅から散乱長を計算し、既存の手法で求められた結果との一致を確認した。加えて、半オフシェル散乱振幅計算に世界で初めて成功した。さらに散乱振幅のパイ中間子質量依存性を調べ、演算子依存性や虚時間依存性な
どの散乱振幅の性質や、相互作用領域外に設定した境界上の波動関数から散乱位相差が得られることを議論した(研究論文 B-3)。

【6】

改良 Wilson クォークを用いた格子 QCD シミュレーションによる有限温度・有限密度 QCDの研究（金谷、谷口）
ビッグバン直後の宇宙初期に実現したと考えられるクォーク・グルオン・プラズマ状態から通常のハドロン物質への相転移前後のクォーク物質の様々な熱力学的性質は、初期宇宙の物質進化や物質創成メカニズムの解明において重要である。これは本質的に非摂動的な問題であり、格子 QCD による QCD 第一原理からの大規模シミュレーションが不可欠である。谷口、
金谷らは、改良 Wilson 型格子クォークを用いた有限温度・有限密度 QCD のシミュレーション
研究を行なっている。2020 年度も、QCD の相構造やクォーク物質の熱力学的諸性質の研究を推進した。クォーク物質の相構造や熱力学特性を引き出すために、有限温度・有限密度格子
QCD を研究し、またそのための技術開発を進めた。

[6-1] Gradient flow に基づく SFtX 法を用いた有限温度 QCD の研究
有限温度・有限密度 QCD のシミュレーション研究の多くは、計算量が少ないスタガード型
格子クォークを用いて行われているが、連続極限で QCD を再現 することが証明されていないという本質的問題を孕んでいる。我々は、理論的基礎が確立している Wilson 型格子クォーク
を用いて QCD 相転移近傍の温度でクォーク物質がどのような熱力学特性を示すかの大規模シ
ミュレーション研究を推進している。Wilson 型クォークは、連続極限の正しさが保証されている反面、有限の格子上ではカイラル対称性を陽に壊しているため、カイラル対称性に関わる物理量に関して格子化誤差が大きく、それを取り除いて物理量を計算するために膨大な計算資源が要求されるという困難があった。また、並進対称性に伴う保存カレントとして定義されるエネルギー運動量テンソルは系の力学特性を調べる上で基本的な観測量だが(例えば、対角成分はエネルギー密度や圧力などの情報を含み、2 点相関関数から様々な粘性率が導かれる)、格子上では連続な並進対称性が陽に壊されているため、従来の方法では、5 種類の演- 18 ― 18 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

算子の非自明なくりこみと混合を非摂動論的に決定するという、複雑なくりこみ操作を行わなければ意味のある評価が出来なかった。
谷口、金谷らは、九州大学鈴木博教授、大阪大学北澤正清助教、新潟大学江尻信司准教授、
広島大学梅田貴士准教授らと、Gradient flow(勾配流)に基づ いて鈴木博らにより開発された SFtX 法(small flow-time expansion method)を応用して、これらの課題を克服した大きなブレークスルーを目指している。Gradient flow とは、仮想的な時間パラメータ t(flow-time)を導入して、作用の勾配で与えられる発展方程式(フロー方程式)により場の量を変形させる
理論的手法である。この発展方程式は一種の拡散方程式になっており、t>0 までフローさせた
結果は、元の場の量を√8t の物理的領域で平準化(smear)させたものと解釈することができる。
さらに、フローさせた場で作る演算子が紫外発散も同一点特異性も持たないという目覚ましい特性を持っていることが、Lüscher と Weisz により証明された。
SFtX 法は、Gradient flow のこの有限性を活用して、連続極限のくりこまれた物理量に対応する量を格子上で評価する一般的な計算方法である。連続理論で何らかの物理量を非摂動論的に評価しようとすると、通常は、その物理量を格子理論で定義し、格子上で評価された値を連続極限まで外挿(a→0)するが、数値的なくりこみに加えて、格子上で重要な対称性が壊されている場合にはそれによる不要な演算子との混合を数値的に除去する必要があり、十分な精度を出すためにはしばしば重い計算となる。SFtX 法では、フローさせた演算子が有限であることを利用して、対応する演算子を格子上で計算することにより、くりこみ操作や混合の除去無しに直接評価する。ただし、フローさせた演算子は求める物理量そのものでは無いので、格子の結果を、連続極限(a→0)とフロー時間ゼロ極限(t→0)に 2 重外挿する。
SFtX はどんな物理量にも使うことができるので、格子化で並進対称性が陽に壊されるためにこれまで扱いが難しかったエネルギー運動量テンソルの格子計算に、最初に応用された。
我々は、SFtX 法が並進対称性の破れだけでなく、Wilson 型クォークによるカイラルの破れの困難にも有効であることに着目し、2+1 フレーバーの動的なクォークを含む現実的 QCD に
SFtX 法を応用した一連の研究を推進している。
我々は、動的クォークを含む QCD への応用の第一段階として、2016-2017 年度に u,d クォークが現実よりやや重い場合の Nf=2+1 QCD の研究を実行した。エネルギー運動量テンソルの対角成分から計算した状態方程式が、従来の方法による結果を良く再現することを示し、さらに、カイラル感受率の disconnected 部分がクロスオーバー温度でピークを示すことをWilson 型クォー クとして初めて示した。また、位相感受率を、グルーオンを用いた定義式と、それを、連続理論のカイラル関係式を用いてクォークを用いて表し直した評価式の両方で計算し、有限の格子間隔でも両者が極めてよく一致することを示した。通常の計算方法でスタガード型クォークを用いた研究では、ここでシミュレーションした程度の格子間隔では- 19 ― 19 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

両者は 2 桁も違っており、SFtX 法が物理的に信頼できる結果を得る上で極めて有用であることをあらわしている。
この研究を発展させ、現実のクォーク質量(物理点)での研究、u,d クォークが重い場合に格子間隔を変えたシミュレーション、及び、
SFtX 法の更なる改良に向けた研究を推進している。
SFtX 法により、物理量の観測に関しては計算時間の大きな削減ができたが、物理点や格子間隔が細かい格子の配位生成には膨大な計算が要求され、様々な計算機資源を動員して大規模シミュレーションを系統的に遂行している。2020 年度には、以下の研究成果が得られた。

SFtX 法の改良
SFtX 法では、格子上で測った t>0 でのフローさせた演算子の結果を、連続極限(a→0)とフ
ロー時間ゼロ極限(t→0)に 2 重外挿する。ここで、t/a2〜0 領域は格子化誤差が大きいので、それを避けながら 2 重外挿する必要がある。従って、格子が粗い、もしくは格子化誤差が大きい演算子の場合には、 外挿に必要な線形領域が不明確になる可能性がある。これらの場合に SFtX 法を応用するためには、なんらかの改善が要求される。研究論文 A-3 では、マッチン
グ係数におけるくりこみスケールと摂動高次項の効果を研究した。
マッチング係数とは、フローさせた演算子と求めるくりこまれた物理量とを結びつける係数である。SFtX 法では、物理量を、格子上で測った t>0 でのフローさせた演算子にマッチン
グ係数を掛けたものを組み合わせてあらわし、その 2 重外挿を実行する。マッチング係数は、
QCD を含む漸近自由な理論では、t が大きくなければ摂動計算することができる。SFtX 法では、摂動であらかじめわかっている t 依存性をマッチング係数に取り込むことで、t→0 外挿が滑らかになるように改良している。実際、我々の u,d クォークが重い Nf=2+1 QCD の研究(格子間隔 0.07fm)では、マッチング係数のおかげで、t→0 外挿を数値的に可能とする線形領域が実現した。他方、現実のクォーク質量での Nf=2+1 QCD の研究(格子間隔 0.09fm)や、u,d ク
ォークが重い Nf=2+1 QCD で格子間隔を粗くした研究において、いくつかの物理量に関して、
t→0 外挿を行うための線形領域が十分明確ではなくなるという現象を経験した。
マッチング係数を摂動計算するときのくりこみスケール μ は、摂動展開が破綻しない限
り、どんな値を取っても良く、t→0 極限の最終的な物理量は μ のとりかたに依存しないはずである。通常 μ としては、フローさせた演算子の自然なスケールの一つである、フローによる smearing スケール μd=1/(√8t)が採用されてきた。上に述べたように、μd と同程度ならば、原理的にどんな μ をとっても良いはずである。しかし実際は、摂動展開を有限次で切っているので、摂動展開の質は、μ をどう取るかで影響を受けるる。特に、SFtX 法で t→0 外
挿を実行する上で、どこまで大きな t を利用できるかは実用上重要である。μd=1/(√8t)の形から分かるように、t を大きくすると赤外領域に近づき、マッチング係数の摂動計算が破綻する。
- 20 ― 20 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

研究論文 A-3 では、Harlander らが最近の摂動の高次計算で提案した μ0=1/√(2eγEt)(γEはオイラー数)の効果を調べた。数値的には μ0≃1.5μd なので、μ0 スケールは、摂動領域をより大きな t まで拡大し、t が大きな領域で t→0 外挿のシグナルを改善してくれる可能性がある。
図 2 に、u,d クォークが現実より重い Nf=2+1 QCD におけるエントロピー密度(左図)と ud ク
ォークのカイラル感受率(右図)の t→0 外挿を示す。 緑が μd スケール、青が μ0 スケールの結果である。t→0 極限の結果はどちらもコンシステントだが、μ0 スケールにより線形シグナ
ルが改善されていることがわかる。現実のクォーク質量での Nf=2+1 QCD での試験研究で、μdスケールでは t→0 外挿が不定な場合でも、μ0 スケールにより線形シグナルが劇的に改善され、SFtX 法の安定性と適用範囲が大きく改善されることを確認した。
研究論文 A-3 では、マッチング係数における 2 ループ項の効果や、2 ループ計算で使われた運動方程式の影響も研究したが、我々の研究している Nf=2+1 QCD では SFtX 法の改良には
大きくは寄与しないらしいことがわかったので、ここでは省略する。

図 2: マッチング係数におけるくりこみスケールの効果。ud クォークが重い Nf=2+1 QCD でのエントロピー密度(左図)と ud クォークのカイラル感受率(右図)。横軸はフロー時間で、直線(点線、ダッシュ線)は、μ0 スケールを用いた t→0 線形外挿(非線形外挿、linear+log 外挿)で、下の矢印がフィット領域。t〜0 近傍のシンボルは、右から、線形外挿、非線形外挿、
linear+log 外挿によるエントロピー密度(左図)と ud クォークのカイラル感受率(右図)の外挿結果(研究論文 A-3)。

クエンチ QCD の潜熱と SFtX 法における 2 重外挿の順序依存性の研究
研究論文 A-5 では、SFtX 法を応用して、クエンチ QCD(SU(3) Yang-Mills 理論)の有限温度
1 次相転移点における潜熱と圧力ギャップ∆p を研究した。1 次相転移点では、共存する 2 相の動的平衡から∆p=0 となるはずで、その数値的確認は結果の信頼性の良い指標とされてきた。
- 21 ― 21 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

また、クエンチ QCD では様々な格子間隔を系統的に計算できるので、上で議論した、a→0 とt→0 の 2 重外挿における極限操作の順番の問題を試験した。SFtX 法における a→0 と t→0 の
2 重外挿では、格子化誤差が大きい t/a2〜0 領域のデータを取り除いて外挿する必要がある。
そして、この t/a2〜0 領域の格子化誤差が正しく取り除かれているならば、2 つの外挿の順序によらず同じ結果が得られると期待される。
研究論文 A-5 では、3 つの格子間隔(Nt=8,12,16 に対応する a=1/(NtTc))、2 つの空間体積(Ns/Nt=6,8 に対応する V=(aNs)3=(Ns/Nt)3/Tc3)でクエンチ QCD シミュレーションを実行し、我々が開発した再重み付け法とヒストグラム法を組み合わせた手法を使って臨界点 T=Tc に調整した。そして、臨界点直上で、SFtX 法を使って臨界点における潜熱と∆p を測定した。
図 3 の左図に、μ0 スケールを使った潜熱と圧力ギャップ∆p の結果を示す。この図から、(a)a→0 と t→0 の極限操作の順番を変えても結果が一致すること、および、(b)∆p がゼロとコンシステントであることが確認される。これにより、連続極限における潜熱を精密で信頼性が高く測定することに成功した。同時に、図 3 の右図に示したように、潜熱の結果に空間
体積依存性確認されるので、より大きな空間体積での追試が望まれる。

SFtX 法を用いた物理点有限温度(2+1)-flavor QCD の研究
これらの研究を発展させ、クォーク質量を現実のクォーク質量(物理点)に調整した Nf=2+1
QCD の熱力学研究を推進している。ゼロ温度配位として PACS-CS Collaboration が生成したものを活用し、それと同じ格子間隔 a≃0.09fm で有限温度シミュレーションを実行している。
物理点の配位生成には大きな計算が要求されるので、様々な計算機資源を動員して系統的な大規模シミュレーションを系統的に遂行している。この物理点シミュレーションでも、上で紹介した SFtX 法の改良が重要な役割を果たすと期待される。
これまでに得た様々な物理量の挙動から、T≃122–146MeV が相転移温度近傍の臨界領域にあることが示唆されるが、まだ統計が十分ではなく、明確な結論を得るためには T≤146MeV の低温領域で更なる統計が必要である。我々の Wilson 型クォークによる相転移温度は、スタガー
ド型クォークによる先行結果より低い可能性があるが、相転移温度の精密な評価は実験的・現象論的にも重要である。現在そのための低温領域の配位生成を集中的に進めており、その中間結果を各種会議で報告した(国際会議発表 3, 国内学会・研究会発表 4,10,11,13,14,15)。

その他
SFtX 法でエネルギー運動量テンソルそのものが計算できるようになったことにより、その非対角要素や二点相関関数も直接評価可能となった。そこから、輸送係数など、クォーク物
質の様々な熱力学特性を引き出す試みを進めている(国内学会・研究会発表 2)。

- 22 ― 22 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

また、SFtX 法を 2 点関数に拡大して、カイラル感受率の connected 部分と U(1)A 感受率の評価を行った。u,d クォークが重い Nf=2+1 QCD では、カイラル感受率の disconnected 部分と違って、connected 部分や全カイラル感受率は QCD 相転移にそれほど敏感では無いことがわかった。これは u,d クォークが重いことによる影響と考えられるので、u,d クォークを物理点まで軽くした Nf=2+1 QCD でも検証する必要がある。また、U(1)A 感受率の結果から、QCD 相
転移温度では U(1)A 対称性がまだ回復していないことを示唆する中間結果を得た(国内学会・研究会発表 5,17)。

[6-2] SFtX 法によるゼロ温度カイラル物理量の研究
SFtX 法により、通常の方法では格子化に伴う対称性の破れのために評価が難しい物理量も直接非摂動論的に評価できる可能性が拓かれた。有限温度の物理量だけでなく、ゼロ温度物
理量にも、格子による非摂動論的評価が望まれる物理量は多く存在する。
谷口、鈴木遊、金谷は、九州大学鈴木博教授と共同で、SFtX 法のゼロ温度物理量への応用
として、B 中間子バッグパラメータ BK を SFtX 法で計算するために必要なマッチング係数の摂
CP 対称性の破れを K0-K0bar 振動実験から抽出する上で鍵となる量で、動計算を行った。
BK は、その非摂動論的評価は極めて重要であるが、格子上のカイラル対称性の破れのために通常の方法では複雑なくりこみと混合の解消が要求される。我々は、SFtX 法をクォーク 4 点関数の評価に拡張して、BK などの評価に必要なマッチング係数を計算した(研究論文 A-4)。

【7】

Wilson 型クォークを用いた 2+1 フレーバー有限温度 QCD における臨界終線の研究（藏増、大野）有限温度 QCD の相構造はクォークの質量やフレーバー数に応じて多様に変化すると予想されており、特に二次相転移となる臨界終点、あるいはこれが連なった臨界終線の位置を決定することが非常に重要である。実際、特に 3 フレーバーの場合について臨界終点の位置を決定するため、格子 QCD に基づく様々な計算が行われてきた。しかしながら、これまでの研究
から臨界終点の位置は格子化誤差の影響を強く受けることが分かってきており、異なるフェ
ルミオン作用を用いた計算結果の間で一貫性のある結果が未だに得られていないなど多くの課題が存在し、より信頼できる連続極限の結果を得ることが求められている。
大野と藏増は、金沢大の武田真滋准教授及び理化学研究所の中村宜文研究員とともに、2+1フレーバーQCD の臨界終点の位置を、O(a)改良された Wilson 型クォーク作用を用いた格子 QCD
シミュレーションにより調べた。本年度はまず、時間方向格子サイズ(NT)を 6 に固定し、以前の研究で調べた SU(3)対称点の近傍からより広いパラメータ領域へと計算を拡張した(研究論文 B-5)。その結果、図 4 に示すような臨界終線を得た。ただし NT=8 及び 10 の結果は、SU(3)- 23 ― 23 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

対称点における NT 依存性をもとに NT=6 の結果を外挿したものである。現在、より正確な結果を得るため、NT=8 の計算を遂行中である。

図 4: 2+1 フレーバーQCD における臨界終線。CEP は Critical End Point(臨界終点)を表す。

【8】

Highly improved staggered quark を用いた有限温度・密度 QCD の研究（大野）有限温度・密度 QCD の相構造について調べることは、強い相互作用する物質の高温・高密度媒質中での性質を理解する上で非常に重要であり、QCD 臨界点の探索等を目的として理論・実験両方の側面から様々研究が行われている。
大野は、独国 Bielefeld 大の Frithjof Karsch 教授を中心とする HotQCD Collaboration に参加
し、highly improved staggered quark 作用を用いた(2+1)-flavor 格子 QCD シミュレーションにより、有限温度・密度 QCD の相構造に関する研究を行った(研究論文 A-6)。本研究ではクォー
ク化学ポテンシャルに対する Taylor 展開を用いて有限密度の効果を取り入れる手法を採用し、
化学ポテンシャルの 8 次のオーダーまで展開係数の計算を行った。これによりバリオン数揺
らぎについて 6 次のキュムラントまで計算することが可能となった。また、この結果を STAR
実験で得られた陽子数揺らぎについてのキュムラントの結果と比較したところ、3 次及び 4
次のキュムラント(歪度及び尖度)に関して、実験結果は若干大きいもののほぼ我々の計算結果と同様の振る舞いをすることが分かった(図 5)。このことから、クォーク・グルーオン・プ
ラズマが再びハドロンに冷え固まる温度(freeze-out 温度)は QCD のクロスオーバーが起こる- 24 ― 24 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

温度(pseudo-critical 温度)の近傍にあり、恐らくより低いところにあると予想された。さらに、
5 次及び 6 次のキュムラントについても同様の比較を行ったが、我々の計算結果から実験結果と同様の振る舞いを得ることはできなかった。その理由としては、Taylor 展開のより高次の効果を取り入れる必要があることが考えられる。
2 [GeV]: 2 0 0s 1/
NN62.4 54.4

RB
nm (Tpc)

1.2

39

27

dashed lines:
joint t toSTAR data for
R P31 , R P42

1
0.8
0.6

NNLO, R B
31 (Tpc)RB
42 (Tpc)
STAR 2020: R p31
R p42

0.4
0.2

RB
12 (Tpc)

STAR preliminary: open symbols
0
0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

図 5: バリオン数揺らぎのキュムラントの格子 QCD 計算結果と陽子数揺らぎのキュムラント
の STAR 実験結果の比較。
【9】

クォーコニウムスペクトル関数の研究（大野）チャームやボトムクォークといった重いクォークとその反クォークの束縛状態であるクォーコニウムのスペクトル関数は、高温媒質中でのクォーコニウムの振る舞いや重クォーク輸送に関する情報をすべて含んでおり、理論的にその性質を調べることは、重イオン衝突実験の結果を理解する上で非常に重要である。しかしながら、格子 QCD に基づく第一原理計算で直接得られる量は相関関数であり、相関関数からスペクトル関数を計算することは ill-posedな問題で、解くことが非常に困難であることが知られている。従って、より信頼できるスペクトル関数を計算するために様々な方法が試みられている。
大野は、中国華中師範大の Heng-Tong Ding 教授、独国 Bielefeld 大の Olaf Kacz- marek 博士
らとともに、クエンチ近似を用いた大規模な格子 QCD シミュレーションにより、連続極限におけるクォーコニウム相関関数を計算し、これを摂動論的モデルにフィットすることでクォーコニウムスペクトル関数を計算した。これまでの計算ではまず輸送ピークの存在しない擬
スカラーチャネルのスペクトル関数を計算を行い、そしてこれを輸送ピークを持つベクター
チャネルのスペクトル関数へと拡張した。今年度はさらにベクターチャネルの相関関数に含まれる輸送の凡その寄与を見積り、チャームクォークに対する抗力係数がボトムクォークのそれよりも大きいという予想を得た(研究論文 B-7)。
- 25 ― 25 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書【10】

テンソルネットワーク形式に基づく格子ゲージ理論の研究（藏増、吉村）
格子 QCD 計算では、近年の計算機能力の向上や新規アルゴリズムの開発・改良の結果、自然界の u、d、s クォーク質量上でのシミュレーションや、更には軽原子核の束縛エネルギー計算までもが可能となりつつある。その一方で、 解決すべき長年の課題がそのまま残されていることも事実である。最も重要な課題は、フェルミオン系を扱う際の負符号問題および複素作用を持つ系のシミュレーションである。これらは、軽いクォークのダイナミクス、Strong
CP 問題、有限密度 QCD、格子 SUSY の研究において避けて通れない問題である。 われわれは、近年物性物理分野で提案されたテンソルネットワーク形式に基づく分配関数の数値計算手法を格子ゲージ理論へ応用し、モンテカルロ法に起因する負符号問題および複素作用問題を解決し、これまでの格子 QCD 計算が成し得なかった新たな物理研究の開拓を目指している。なお、本研究課題は、ポスト「京」で重点的に取り組むべき社会的・科学的課題における
4 つの萌芽的課題のうち、1 番目の課題である「基礎科学のフロンティア − 極限への挑戦」に含まれており、本グループも分担機関として参加し、テンソルネットワーク法の素粒子物理学への応用に取り組んできた。萌芽的課題は 2019 年度末で終了したが、2020 年度以降は科学研究費基盤(A)を新たな資金として研究の発展を図っている。
2014 年、藏増と理研計算科学研究機構(現理研計算科学研究センター)の清水特別研究員は、テンソル繰り込み群をグラスマン数も扱えるように拡張し(グラスマンテンソル繰り込み群)、
世界で初めてフェルミオン入りのゲージ理論への応用に成功した。具体的には、グラスマンテンソル繰り込み群を用いて、θ項が有る場合と無い場合の 1 フレーバーの 2 次元格子
Schwinger モデル(2 次元格子 QED)における相構造を調べた(論文発表済)。この研究により、グラスマンテンソル繰り込み群が、現在の格子 QCD 計算が抱える負符号問題や複素作用問題を解決していることを示すことに成功した。今後は、最終目標である 4 次元 QCD への応用に向け、(i)非可換ゲージ理論への拡張、(ii)高次元モデルへの応用、(iii)物理量計算のための手法
開発、(iv)興味深い低次元素粒子論モデルへの応用、という 4 つの課題に取り組む。
2020 年度における課題(i)~(iv)のうち、特に課題(ii)において大きな進展があった。一般的に、テンソル繰り込み群はモデルの次元が上がるにつれて計算コストが増大する。そのため、これまでテンソル繰り込み群の主な応用例は 2 次元モデルに限られており、4 次元モデルへの適用例は、2019 年にわれわれが行った 4 次元イジングモデルの計算のみであった(論文発表済)。われわれは、まず Anisotropic TRG(ATRG)コードを並列化し、4 次元有限密度複素スカラー理論における Silver Blaze 現象の確認に成功した。図 6(左)は粒子数密度⟨n⟩の化学ポテンシ
ャル μ 依存性を 4 次元体積を変化させながらプロットしたものである。⟨n⟩はある化学ポテン
シャル μc の値までゼロのままであるが、μc≈0.65 を越えると突然増加し始めることがわかる(Silver Blaze 現象)。なお、平均場近似では、μMF≈0.7 と予想されている。このモデルは複素作
用問題の典型例であり、テンソル繰り込み群が複素作用問題・符号問題を解決していること- 26 ― 26 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書がわかる(研究論文 A-8)。このモデルは、複素 Langevin 法や Lefscetz-thimble 法でも計算されているが、いずれも O(104)程度の格子体積でしか成功しておらず、ゼロ温度・熱力学極限とは程遠い。4 次元有限密度複素スカラー理論の計算に成功した後、われわれは並列版 ATRG
コードをグラスマン数化し、4 次元有限密度 Nambu−Jona-Lasinio モデルの極低温・高密度領域におけるカイラル相転移を調べた。図 6(右)はオーダーパラメータであるカイラル凝縮⟨χχ⟩の化学ポテンシャル μ 依存性を表している。μ≈3.0 において⟨χχ⟩の有限の値からゼロへの不連続な変化が確認できるため、カイラル対称性の回復が一次相転移で起きていることがわかる(研究論文 A-9)。

図 6: (左)4 次元有限密度複素スカラー理論における粒子数密度⟨n⟩の化学ポテンシャル μ 依存性。(右)4 次元有限密度 Nambu−Jona-Lasinio モデルにおけるカイラル凝縮⟨χχ⟩の化学ポテンシ
ャル μ 依存性。m はフェルミオンの質量を表す。
【11】

素粒子標準模型を超えた理論の探索（山崎）ウォーキングテクニカラー模型は素粒子標準模型を超えた理論の有力な候補の一つである。
この模型は、強結合ゲージ理論のダイナミクスにより、素粒子標準模型では手で与えられていた電弱対称性の自発的破れの起源を説明できる可能性がある。しかし、この模型を構築するために必要な強結合ゲージ理論には、近似的共形対称性を持つなど、特殊な条件が課されている。山崎は名古屋大学山脇幸一名誉教授や理研計算科学研究センター 青木保道チームリーダーらと共に、LatKMI Collaboration において、格子ゲージ理論を用いた数値計算から、そのような条件を満たすゲージ理論が存在するかの探索を行っている。
これまでの 4、8、12 フレーバーSU(3)ゲージ理論の研究から、8 フレーバー理論がそれら条件を満たす可能性があることを示唆した。2020 年度はフレーバー1 重項スカラー中間子及びフレーバー1 重項擬スカラー中間子の質量について、4、8、12 フレーバー理論の比較及び、8フレーバー理論から得られる S パラメータ計算を行った。

- 27 ― 27 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書【12】

有限サイズスケーリングによるコンフォーマル対称性の検証の研究（吉江）コンフォーマル対称性は、素粒子標準模型を超える理論の構築に重要な働きをする、と考えられており、場の理論の非摂動的定式化である格子場理論とその数値シミュレーションに
よって、多フレーバ QCD などのコンフォーマル理論の候補が調べられている。数値シミュレーションは有限サイズの格子で行われるので、コンフォーマル対称性は厳密には成り立たず、
格子サイズの有限性に起因する相構造が現れる。この『有限の IR cutoff を持つコンフォーマ
ル理論の相構造』の研究は、この 10 年ほどの間に大きく進展し、繰り込み群の帰結である『格子サイズに対する相関関数のスケーリング則』を調べて理論が(体積無限大で)コンフォーマ
ル対称性を持つ(赤外固定点が存在すること)か否かを判別できる様になってきた。
これらの一連の研究は、素粒子物理にとって興味深い、多フレーバ QCD に対して行われたものであるが、同様な議論が、別の(遙かに簡単な)系でも成り立つことが示されれば、上述の論理構造に穴がないことの傍証となる。この目的で、吉江と三好は、2 次元 Ising 場理論のコンフォーマル対称性を有限体積スケーリンに基づき調べた。2 次元 Ising 模型は、相関距離無限大の critical surface 上では、コンフォーマル対称性を持つ場の理論が実現している、と考えられている。実際、隣接相互作用のみを持つ Ising 模型は Onsager の critical point でコンフォーマル対称性を持つことが解析的に示されている。我々は、まず、critical point 上でいくつかの格子サイズでの相関関数を数値的に求め、期待されるスケーリング則が成り立っていることを確かめた。最近接と次に近接するスピン間の相互作用を含む Ising 模型で critical point の知られている幾つかを選んで、同様な計算を行い、近接作用の場合と同様のスケーリング則が成り立っていることを示したが、critical surface 上の赤外固定点の位置に関する情報は、統計不足の為、得られなかった。これらの結果は、三好の修士論文に纏めている。Critical Surface上の有限サイズスケーリングの方法は、繰り込み群の beta 関数(繰り込み群のフロー)や gamma関数を決める方法を含んで おり、コンフォーマル理論の研究手段として有用なものになり得る。その方向での研究を計画中である。

【13】

格子 QCD 研究用データグリッド ILDG/JLDG の運用（吉江、大野）JLDG(Japan Lattice Data Grid)は、計算素粒子物理研究者の為のデータグリッドであり、国内7 拠点を SINET5 VPN で結んだネットワークにより、拠点のファイルサーバを統合した単一のファイルシステムを運用している。JLDG は、研究者が研究グループ内で日々のデータを共有する、複数のスパコンでデータ生成・解析を分担する、基礎データを ILDG(International Data
Grid)の枠組みで世界の研究者に公開する、等に利用されている。本学計算科学研究センター
が代表を務める、拠点代表やユーザー代表から構成される JLDG 管理者グループが運用の一切を行っている。2008 年に実運用を開始して以来 12 年経過した現在、国内の複数の大きな- 28 ― 28 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

研究グループが研究インフラとして使用している。JLDG は実用システムとして、数年前から一定の完成の域に達しており、今年度もシステムの増強・安定運用を主眼に活動を行った。
日常のメンテナンス・ユーザ対応以外の主な活動は、以下の通りである。
(a) 新拠点設置:『理化学研究所計算科学研究センター(RCCS)』拠点の立ち上げ完了。HPCI 最大の計算資源である富岳のクライアントで JLDG ファイルシステムをマウントし、データ転送を並列に行える様になった。
(b) 筑波大 CCS 拠点のスーパーコンピュータ Cygnus に JLDG クライアント機能を導入する計画が完了した。JLDG ファイルシステムを Cygnus フロントエンドでマウント利用可能である。
(c) 格子 QCD アンサンブルの DOI 登録の枠組みが完成し、PACS-CS Collaboration 提供の 6 アンサンブルに DOI が付与された。
(d) 管理用サーバの OS を EOL を迎えた CentOS6 から。CentOS7/8 へ更新した。
(e) 筑波大キャンパスネットワークのセキュリティゾーン化の一環として、センターより提供された『外部レベル相当ネットワーク』の一つのサブネットへの機器移転を行った。その際、センターの運用するファイアウォールシステムに登録するファイアウォールルールを検討し、システムのセキュリティレベル向上を図った。
教育
4.教[修士論文]
1. 井元 航希「テンソル繰り込み群を用いた O(3)非線形シグマモデルの解析」
2. 三好 拓己
「有限体積スケーリングによる 2 次元イジング場理論のコンフォーマル対称性の検証」[博士論文]
1. 馬場 惇「Measurement of Chiral Susceptibility in Lattice QCD with Wilson Fermion using
Gradient Flow」
5.受受賞、外部資金、知的財産権等受賞
1.

日本物理学会学生優秀発表賞(素粒子論領域)
秋山

進一郎「テンソル繰り込み群による 3+1 次元有限密度 Nambu–Jona-Lasinio 模型の研究」- 29 ― 29 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

2020 年 9 月
外部資金
1.

藏増嘉伸（代表）、科学研究費補助金・基盤研究(A)、令和 2 年度採択、「テンソルネットワーク法による計算物理学の新展開」、10,300 千円

2.

金谷和至（代表）、科学研究費補助金・基盤研究(C)、平成 31 年度採択、「グラジエントフローによるクォーク・ハドロン物質の熱力学特性」、1,200 千円

3.

金谷和至（分担）、科学研究費補助金・基盤研究(B)、令和 2 年度採択、「有限温度
QCD における物理量の決定へ向けて」、100 千円

4.

谷口裕介（代表）、科学研究費補助金・基盤研究(C)、平成 30 年度採択、「勾配流法を用いたクォーク・グルーオンプラズマの物性的研究」、700 千円

5.

山崎剛（代表）、科学研究費補助金・基盤研究(B)、平成 31 年度採択、「強い相互作
用の第一原理計算による軽ハドロン形状因子の総合理解」、3,400 千円

6.

山崎剛（代表）、二国間交流事業(共同研究・セミナー)、令和 2 年度採択、「真空および媒質中での格子 QCD によるハドロン測定の新しい可能性」、2,375 千円、7.

浮田尚哉（代表）、科学研究費補助金・基盤研究(C)、令和 2 年度採択、「機械学習を取り入れた格子 QCD による超精密物理量測定のための計算手法の開発と実践」、
1,200 千円

知的財産権

6.研
研究業績
A) 研究論文A) 査読付き論文
1.

PACS Collaboration: J. Kakazu, K.-I. Ishikawa, N. Ishizuka, Y. Kuramashi, Y. Nakamura, Y.
Namekawa, Y. Taniguchi, N. Ukita, T. Yamazaki, and T. Yoshié, “Kl3 form factors at thephysical point on (10.9fm)3 volume”, Phys. Rev. D101 (2020) no.9, ref.094504.

2.

PACS Collaboration: E. Shintani, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, and T. Yamazaki,“Erratum: Nucleon form factors and root- mean-square radii on a (10.8fm)4 lattice at thephysical point”, Phys. Rev. D102 (2020) no.1, ref.019902.

3.

WHOT-QCD Collaboration: Yusuke Taniguchi, Shinji Ejiri, Kazuyuki Kanaya, Masakiyo
Kitazawa, Hiroshi Suzuki, and Takashi Umeda, “Nf=2+1 QCD thermodynamics with gradient
flow using two-loop matching coefficients”, Phys. Rev. D102 (2020) no.1, ref.014510.
- 30 ― 30 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

4.

Asobu Suzuki, Yusuke Taniguchi, Hiroshi Suzuki, and Kazuyuki Kanaya, “Four quark
operators for kaon bag parameter with gradient flow”, Phys. Rev. D102 (2020) no.3,
ref.034508

5.

Mizuki Shirogane, Shinji Ejiri, Ryo Iwami, Kazuyuki Kanaya, Masakiyo Kitazawa, Hiroshi
Suzuki, Yusuke Taniguchi, and Takashi Umeda, “Latent heat and pressure gap at the firstorder deconfining phase transition of SU(3) Yang-Mills theory using the small flow-time
expansion method”, Prog. Theor. Exp. Phys. 2021 (2021) ref.013B08.

6.

HotQCD Collaboration: A. Bazavov, D. Bollweg, H.-T. Ding, P. Enns, J. Goswami, P. Hegde,O. Kaczmarek, F. Karsch, R. Larsen, Swagato Mukherjee, H. Ohno, P. Petreczky, C. Schmidt,
S. Sharma, and P. Steinbrecher, “Skewness, kurtosis, and the fifth and sixth order cumulantsof net baryon-number distributions from lattice QCD confront high-statistics STAR data”,
Phys. Rev. D101 (2020) no.7, ref.074502.

7.

Y. Kuramashi and Y. Yoshimura, “Tensor renormalization group study of two-dimensional
U(1) lattice gauge theory with a θ term”, JHEP 2004 (2020) ref.089.

8.

S. Akiyama, D. Kadoh, Y. Kuramashi, T. Yamashita, and Y. Yoshimura, “Tensor
renormalization group approach to four-dimensional complex φ4 theory at finite density”,
JHEP 2009 (2020) 177.

9.

S. Akiyama, Y. Kuramashi, T. Yamashita, and Y. Yoshimura, “Restoration of chiral symmetry
in cold and dense Nambu-Jona-Lasinio model with tensor renormalization group”, JHEP 2101
(2021) 121.

B) 査読無し論文
1.

PACS Collaboration: J. Kakazu, K.-I. Ishikawa, N. Ishizuka, Y. Kuramashi, Y. Nakamura, Y.
Namekawa, Y. Taniguchi, N. Ukita, T. Yamazaki, and T. Yoshié, “Kl3 form factors in Nf=2+1
QCD at physical point on large volume”, PoS(LATTICE2019) (2020) ref.186, pp.1-7.

2.

PACS Collaboration: Natsuki Tsukamoto, Yasumichi Aoki, Ken-ichi Ishikawa, Yoshinobu
Kuramashi, Shoichi Sasaki, and Takeshi Yamazaki, “Nucleon isovector couplings from 2+1
flavor lattice QCD at the physical point”, PoS(LATTICE2019) (2020) ref.132, pp.1-7.

3.

Takeshi Yamazaki and Yusuke Namekawa, “Two-pion scattering amplitude from BetheSalpeter wave function at the interaction boundary”, PoS(LATTICE2019) (2020) ref.032,
pp.1-7.

4.

PACS Collaboration: Y. Aoki, Y. Kuramashi, E. Shintani, and N. Tsukamoto, “Proton decay
matrix elements with physical quark masses”, PoS(LATTICE2019) (2020) ref.141, pp.1-7.

- 31 ― 31 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

5.

Y. Nakamura, Y. Kuramashi, H. Ohno, and S. Takeda, “Critical endpoint in the continuum
limit and critical endline at NT=6 of the finite temperature phase transition of QCD with clover
fermions”, PoS(LATTICE2019) (2020) ref.053, pp.1-7.

6.

HotQCD Collaboration: Heng-Tong Ding, Prasad Hegde, Olaf Kaczmarek, Frithjof Karsch,
Anirban Lahiri, Sheng-Tai Li, Swagato Mukherjee, Hiroshi Ohno, Peter Petreczky, Chris- tianSchmidt, and Patrick Steinbrecher, “The chiral phase transition temperature in (2+1)-flavor
QCD”, 18th International Conference on Hadron Spectroscopy and Structure, pp.672-677.

7.

A.-L. Lorenz, H.-T. Ding, O. Kaczmarek, H. Ohno, H. Sandmeyer, and H.- T. Shu, “Thermal
modifications of quarkonia and heavy quark diffusion from a comparison of continuumextrapolated lattice results to perturbative QCD”, PoS(LATTICE2019) (2020) ref.207, pp.17.

8.

S. Akiyama, Y. Kuramashi, T. Yamashita, and Y. Yoshimura, “Phase transition of fourdimensional Ising model with tensor network scheme”, PoS(LATTICE2019) (2020) ref.138,
pp.1-7.

B ) 国際会議発表
A) 招待講演1.

H. Ohno「In meedium properties of heavy quarks on the lattice」, Extreme Nonequilibrium
QCD (Online, Oct. 5-9, 2020).

B) 一般講演
1.

R. Tsuji, Y. Aoki, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, E. Shintani, and T. Yamazaki「Nucleon structure at physical point from 2+1 flavor Lattice QCD」, Asia-Pacific Symposium
for Lattice Field Theory (APLAT2020) (Online, Aug. 4-7, 2020).

2.

S. Akiyama, D. Kadoh, Y. Kuramashi, T. Yamashita, and Y. Yoshimura 「 Tensor
renormalization group approach to four-dimensional complex φ4 theory at finite density」,
Asia-Pacific Symposium for Lattice Field Theory (APLAT2020) (Online, Aug. 4-7, 2020).

3.

Kazuyuki Kanaya「Thermodynamic observables in (2+1)-flavor QCD applying the gradientflow method」, CCS 12th international symposium 2020 on Discovery, Fusion, Creation of
New Knowledge by Multidisciplinary Computational Sciences (CCS, Tsukuba, Japan (online),
Oct. 6, 2020).

4.

R. Tsuji, Y. Aoki, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, E. Shintani, and T. Yamazaki「Nucleon structure with 2+1 flavor lattice QCD at physical quark mass」, International
School for Strangeness Nuclear Physics (SNP school 2020) (J- PARC, Tokai, Dec. 2-5, 2020).
- 32 ― 32 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

5.

Shininchiro Akiyama「Restoration of chiral symmetry in cold and dense Nambu-Jona-Lasiniomodel with tensor renormalization group」, KEK Theory Workshop 2020 (Online, Dec. 1518, 2020).

6.

Shininchiro Akiyama「Restoration of chiral symmetry in cold and dense Nambu-Jona-Lasiniomodel with tensor renormalization group」, YITP workshop on Quantum Computing forquantum field theories (Online, Jan. 5-8, 2021).

C ) 国内学会・研究会発表
A ) 招待講演1.

辻竜太朗, 青木保道, 石川健一, 藏増嘉伸, 佐々木勝一, 新谷栄吾, 山崎剛「格子場の理論に基づく核子構造研究」, 日本のスピン物理学の展望 (島根 (online),
2021 年 2 月 23-24 日).

B ) その他の発表
1.

金谷和至「Thermodynamics of 2+1 flavor QCD with the SFtX method based on the gradientflow」, 理化学研究所計算科学研究センター(R-CCS)セミナー (理化学研究所, 神戸市,
兵庫県 (online), 2020 年 5 月 20 日).

2.

谷口裕介, 金谷和至 (Y. Taniguchi, K. Kanaya, A. Baba, S. Ejiri, M. Ki- tazawa, H. Suzuki,
T. Umeda)「物理的なクォーク質量におけるエネルギー運動量テンソルの研究」, 学際大規模情報基盤共同利用・共同研究拠点(JHPCN)第 12 回拠点シンポジウム (THEGRAND HALL, 品川区, 東京都 (online), 2020 年 7 月 9 日).

3.

鈴木博 (Hiroshi Suzuki, Kazuyuki Kanaya, Yusuke Taniguchi, Shinji Ejiri, Takashi Umeda,
Masakiyo Kitazawa)「有限温度量子色力学のダイナミクス」, 学際大規模情報基盤共
同利用・共同研究拠点(JHPCN)第 12 回拠点シンポジウム(THE GRAND HALL, 品川
区, 東京都 (online), 2020 年 7 月 9 日).

4.

金谷和至「Improvement of the SFtX method based on the gradient flow in the study of finite
temperature Nf=2+1 QCD」, 熱場の量子論とその応用 2020(TFQT 2020) (KEK, Tsukuba,
Japan (online), 2020 年 8 月 24 日-26 日).

5.

馬場惇「SFtX 法を用いたカイラル感受率の測定」, 熱場の量子論とその応用
2020(TFQT 2020) (KEK, Tsukuba, Japan (online), 2020 年 8 月 24 日-26 日).

6.

秋山進一郎「テンソル繰り込み群による 4 次元有限密度複素スカラー場の理論の研究」, 熱場の量子論とその応用 2020(TFQT 2020) (KEK, Tsukuba, Japan (online), 2020
年 8 月 24 日-26 日).

- 33 ― 33 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

7.

秋山進一郎「テンソル繰り込み群による 4 次元有限密度複素スカラー場の理論の研究」, 素粒子若手オンライン研究会 2020 (Online, 2020 年 8 月 26 日-28 日).

8.

辻竜太朗, 青木保道, 石川健一, 藏増嘉伸, 佐々木勝一, 新谷栄吾, 山崎剛「物理点2+1 フレーバー格子 QCD による核子構造研究」, 日本物理学会 2020 年秋季大会 (筑波大学, つくば市, 茨城県 (online), 2020 年 9 月 14 日-17 日).

9.

馬場惇, 梅田貴士, 江尻信司, 金谷和至, 北沢正清, 鈴木遊, 鈴木博, 谷口裕介,「SFtX 法を用いたカイラル感受率の測定」, 日本物理学会 2020 年秋季大会 (筑波大学, つくば市, 茨城県 (online), 2020 年 9 月 14 日-17 日).

10. 金谷和至, 梅田貴士, 江尻信司, 北沢正清, 鈴木博, 谷口裕介「Improvement of the
SFtX method based on the gradient flow in the study of finite temperature Nf=2+1 QCD」,
日本物理学会 2020 年秋季大会 (筑波大学, つくば市, 茨城県 (online), 2020 年 9 月 14日-17 日).
11. 鈴木博, 梅田貴士, 江尻信司, 金谷和至, 北沢正清, 鈴木遊, 谷口裕介, 馬場 惇
「Small Flow time eXpansion(SFtX)法による 2+1 フレーバーQCD の熱力学」, 日本物理学会 2020 年秋季大会 (筑波大学, つくば市, 茨城県 (online), 2020 年 9 月 14 日-17日).
12. 秋山進一郎, 藏増嘉伸, 山下巧, 吉村友佑「テンソル繰り込み群による 3+1 次元有限
密度 Nambu–Jona-Lasinio 模型の研究」, 本物理学会 2020 年秋季大会 (筑波大学, つくば市, 茨城県 (online), 2020 年 9 月 14 日-17 日).
13. 谷口裕介, 金谷和至「勾配流法を用いたクォーク・グルオンプラズマの物性的研究」,
第 7 回「京」を中核とする HPCI システム利用研究課題成果報告会 (THE GRAND
HALL, 品川区, 東京都 (online), 2020 年 10 月 30 日).
14. 鈴木博, 金谷和至, 谷口裕介, 江尻信司, 梅田貴士, 北澤正清, 馬場惇「ウィルソンフ
ェルミオンに基づいた有限温度量子色力学の研究」, 第 7 回「京」を中核とする HPCI
システム利用研究課題成果報告会 (THE GRAND HALL, 品川区, 東京都 (online),
2020 年 10 月 30 日).
15. 金谷和至「Gradient flow に基づく SFtX 法による物理点 QCD の熱力学特性の研究」,
大阪大学サイバーメディアセンター2020 年度公募型利用制度成果報告会 (大阪大学
吹田キャンパス, 茨木市, 大阪府 (online), 2021 年 3 月 10 日).
16. 辻竜太朗, 青木保道, 石川健一, 藏増嘉伸, 佐々木勝一, 新谷栄吾, 山崎剛「核子内クォークの担う運動量及びヘリシティ割合の物理点格子 QCD 計算」, 日本物理学会第76 回年次大会 (Online, 2021 年 3 月 12 日-15 日).

- 34 ― 34 ―

筑波大学 計算科学研究センター 令和二年度 年次報告書

17. 馬場惇, 梅田貴士, 江尻信司, 金谷和至, 北沢正清, 鈴木遊, 鈴木博, 谷口裕介「SFtX
法を用いた Wilson fermion の下でのカイラル感受率の測定」, 日本物理学会第 76 回
年次大会 (Online, 2021 年 3 月 12 日-15 日).
18. 江尻信司, 白銀瑞樹, 石見涼, 金谷和至, 北沢正清, 鈴木博, 谷口裕介, 梅田貴士
「Small flow-time expansion 法による SU(3)格子ゲージ理論の一次相転移点における潜熱と圧力差」, 日本物理学会第 76 回年次大会 (Online, 2021 年 3 月 12 日-15).
19. 秋山進一郎, 藏増嘉伸, 吉村友佑「4 次元φ4 理論のテンソル繰り込み群による解析」,
日本物理学会第 76 回年次大会 (Online, 2021 年 3 月 12 日-15 日).
D ) 著書、解説記事等
7.異異分野間連携・産学官連携・国際連携・国際活動等
異分野間連携（センター内外）
1. 計算基礎科学連携拠点
http://www.jicfus.jp/jp/
産学官連携
国際連携・国際活動
1. International Lattice Data Grid (ILDG)
http://ildg.sasr.edu.au/Plone
2. Japan Lattice Data Grid (JLDG)
http://www.jldg.org/jldg/, http://ws.jldg.org/QCDArchive/index.jsp
8.シ
シンポジウム、研究会、スクール等の開催実績
なし
管理・運営
9.管
1.

藏増嘉伸

運営委員会委員、運営協議会委員

2.
3.

吉江友照
吉江友照

共同研究運用委員会委員
藏増嘉伸 計算機システム運用委員会委員

10.社
社会貢献・国際貢献
11.そその他
海外長期滞在、フィールドワークなど
該当なし- 35 ― 35 ―