[論文レビュー] Local Discontinuous Galerkin for the Functional Renormalisation Group

Friederike Ihssen, Jan M. Pawlowski+2 more|arXiv (Cornell University)|Jul 25, 2022

Advanced Numerical Methods in Computational Mathematics被引用数 5

ひとこと要約

本稿では、O(N)-模型における局所的ポテンシャル近似（LPA）における関数的ランダム化群（fRG）フロー方程式を解くための局所不連続ガラーキン（LDG）有限要素法を提案する。従来の拡散に弱いスキームに代えてLDGを導入することで、特に一次相転移付近において優れた数値的安定性を達成し、場の空間におけるショック形成の検出が信頼できるものとなった。実装はオープンソースであり、DUNE PDEフレームワークに統合されている。

ABSTRACT

We apply the Local Discontinuous Galerkin discretisation to flow equations of the O(N)-model in the Local Potential Approximation. The improved stability is directly observed by solving the flow equation for various $N$ and space-time dimensions $d$. A particular focus of this work is the numerical discretisation and its implementation. The code is publicly available, and is explained in detail here. It is realised as a module within the high performance PDE framework DUNE.

研究の動機と目的

• fRGフロー方程式における数値的不安定性、特にショックが形成される一次相転移付近の問題を解決する。
• 過剰な数値的拡散を引き起こす既存の不連続ガラーキンスキームの限界を克服する。
• ショック形成や競合する秩序項効果を含む非自明な場のダイナミクスを捉える、強固で高精度なfRG用数値フレームワークを開発する。
• 複雑なポテンシャル構造を有するスカラー場理論における相境界や臨界的挙動を信頼性高く計算可能にする。
• fRG研究分野における広範な採用を支援するため、DUNE PDEフレームワークに統合された公開可能でモジュラーな実装を提供する。

提案手法

• d=0–4次元におけるO(N)-模型のfRGフロー方程式に由来する2階偏微分方程式に、局所不連続ガラーキン（LDG）離散化を適用する。
• 2階微分を補助変数を用いて扱う混合有限要素式を採用することで、標準的なDG法よりも安定性を向上させる。
• 保存則を保ち、人工的拡散を低減するためのフラックス再構築および数値フラックス関数を実装する。
• 得られる方程式系の高性能かつ並列処理可能な解法のために、DUNE PDEフレームワークを活用する。
• 流れの進化中に急峻な特徴を解像できるように、適応的ステップサイズを用いた時間積分スキームを導入する。
• ゼロ次元QFTおよび非自明な初期ポテンシャル（セクシックおよびオクティック項を含む）を有するO(N)-模型に対して、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ショック形成、特に一次相転移付近において、LDG法がfRGフロー方程式を安定化できるか。
• RQ2 数値的拡散および収束性の観点から、LDG法は従来の直接不連続ガラーキン（DDG）スキームと比べてどのように異なるか。
• RQ3 高次結合項（例：λ₆, λ₈）は、有限Nにおけるマスのσモードにおけるショック発展にどのような役割を果たすか。
• RQ4 本手法は、パイオンおよびシグマ質量の流れにおける非解析的特徴やキック型構造を信頼性高く解像できるか。
• RQ5 LDGアプローチは、他のfRG截断や場の理論へ一般化可能か。

主な発見

• LDG法は、従来のDDGスキームと比較して数値的拡散を顕著に低減し、fRGフロー方程式における安定性および収束性が向上した。
• 初期ポテンシャルにλ₆ = -0.5およびλ₈ = 0.25を用いることで、マスのσモードにショックが観測され、有限NにおけるO(N)-模型におけるショック形成が確認された。
• パイオン質量mπ²のキック構造が、シグマ質量mσ²においてショックに発展するという、モード間の非線形結合を示した。
• 本手法は、凸性回復ダイナミクスを的確に捉え、流れ全体を通してk² + (u + 2ρv)の正値性を維持し、数値的安定性を確保した。
• 極めて小さなRG時間ステップ（≤10⁻¹⁰）でさえも、数値的に安定であり、相転移付近の微細構造の解像が可能となった。
• オープンソースコードはhttps://github.com/satfra/dune-FRGDGにて公開されており、再利用および拡張を可能にする詳細なドキュメンテーションを併記している。