[論文レビュー] Computational General Relativity in the Wolfram Language using Gravitas I: Symbolic and Analytic Computation

Jonathan Gorard|arXiv (Cornell University)|Aug 15, 2023

Algebraic and Geometric AnalysisMathematics被引用数 3

ひとこと要約

本論文は、Wolfram言語で記述されたオープンソースの計算フレームワーク「Gravitas」を紹介する。このフレームワークは、記号的・解析的一般相対性理論を実装し、記号的テンソル計算、任意の座標系、幾何構造をシームレスに統合可能である。広範な時空幾何構造および物質場においてアインシュタイン方程式やアインシュタイン＝マクスウェル方程式の自動解法を可能とし、数値的時間発展のための新規なハイパーグラフベースの適応的メッシュ細分化機構を備える。これは、統合的記号的・数値的相対性理論計算分野における顕著な前進を示している。

ABSTRACT

We introduce a new, open-source computational general relativity framework for the Wolfram Language called Gravitas, which boasts a number of novel and distinctive features as compared to the many pre-existing computational and numerical relativity frameworks currently available within the open-source community. These include, but are not limited to: seamless integration of its powerful symbolic and numerical subsystems, and, by extension, seamless transition between analytic/continuous representations and numerical/discrete representations of arbitrary spacetime geometries; highly modular, general and extensible representations of spacetime geometries, spacetime topologies, gauge conditions, coordinate systems, matter fields, evolution equations and initial data; ability to set up and run complex numerical relativity simulations, and to perform 2D and 3D visualizations, symbolic computations and numerical analysis (including the extraction of gravitational wave signals) on the resulting data, all from within a single notebook environment; and a totally-unstructured adaptive refinement scheme based on hypergraph rewriting, allowing for exceedingly efficient discretization and numerical evolution of Cauchy initial data for a wide range of challenging computational problems involving strong relativistic field dynamics. In this first in a series of two articles covering the framework, we focus on the design and capabilities of Gravitas's symbolic subsystem, including its general and flexible handling of arbitrary geometries parametrized by arbitrary curvilinear coordinate systems (along with an in-built library of standard metrics and coordinate conditions), as well as its various high-level tensor calculus and differential geometry features. We proceed to show how this subsystem can be used to solve the Einstein field equations both analytically and numerically.

研究の動機と目的

• 任意の座標系および時空幾何構造において、アインシュタインの場の方程式を解析的にかつ数値的に解くという長年の課題に取り組む。
• 既存のフレームワークの限界を克服し、1つのノートブック環境内で記号的計算と適応的数値的相対性理論を統合する。
• 任意のエネルギー・物質分布および電磁場に対して、アインシュタイン方程式およびアインシュタイン＝マクスウェル方程式の自動的で一般化可能な解法を可能にする。
• 高度な相対論的理論や今後の数学的相対性理論および量子重力の発展を支援できる、非常に拡張性の高いモジュラーなフレームワークを開発する。
• 適応的メッシュ細分化と強力な可視化機能を備えた次世代の数値的相対性理論シミュレーションの基盤を築く。

提案手法

• 任意の曲線座標系および一般時空幾何構造を扱えるWolfram言語における記号的テンソル計算サブシステムを実装する。
• 標準的な計量、ゲージ条件、運動量テンソル、エネルギー条件の組み込みライブラリを統合し、場の方程式の自動設定を可能にする。
• Wolfram言語の記号的エンジンを活用して、リーマン、リッチ、アインシュタインテンソルの計算を含む高レベルなテンソル演算を実行する。
• ハイパーグラフ書き換えを用いて、完全に非構造的で適応的な時空離散化の細分化を実現する。
• 記号的表現と数値的表現の間でのシームレスな移行を可能にするモジュラーなアーキテクチャを設計する。
• 専用の記号的および数値的ソルバーを介して、相対論的電磁気学および物質場モデルの統合を支援する。

実験結果

リサーチクエスチョン

• RQ1 記号的および数値的一般相対性理論を、任意の幾何構造において統合的に解くことのできる、拡張性のある計算フレームワークとして実現するにはどうすればよいか？
• RQ2 ハイパーグラフベースの適応的細分化スキームは、従来の構造的メッシュ手法に比べて、強力な場の相対的力学をシミュレートする際にどの程度優れているか？
• RQ3 Wolfram言語における記号的サブシステムは、複雑な物質および場の配置に対して、アインシュタイン＝マクスウェル方程式の導出と解法を自動化できるか？
• RQ4 このフレームワークのモジュラー設計は、アインシュタイン＝カルタン理論やスカラー・テンソル理論といった高度な理論への拡張をどの程度可能にするか？
• RQ5 記号的テンソル計算と適応的ハイパーグラフベースの離散化を組み合わせることで、どのような新しい数値的相対性理論の能力が生じるか？

主な発見

• Gravitasは、標準的な計量および運動量テンソルの組み込みサポートを備えた、任意の座標系においてアインシュタイン方程式およびアインシュタイン＝マクスウェル方程式の完全な自動記号的計算を可能にする。
• このフレームワークは、時空の記号的表現と数値的表現との間でシームレスな移行を実現し、初期データから重力波抽出までをエンドツーエンドで計算可能にする。
• ハイパーグラフベースの適応的細分化スキームにより、カウチ初期データの効率的かつ柔軟な離散化が可能となり、複雑で非自明な時空トポロジーや座標構造をサポートする。
• 記号的サブシステムは、リーマン、リッチ、アインシュタインテンソルの計算をネイティブに処理でき、今後の拡張でフレーム場やスピン接続といった高度な構造もサポートする。
• このフレームワークは、2次元および3次元の可視化、重力波信号抽出を含む、数値的相対性理論シミュレーションの自動設定と実行を支援する。
• 記号的および数値的サブシステムの統合により、複雑な相対論的問題を解く際の人的作業を大幅に削減できる高レベルで一般化可能なワークフローが実現される。