[論文レビュー] A Semidefinite Program Solver for the Conformal Bootstrap

David Simmons–Duffin|arXiv (Cornell University)|Feb 6, 2015

Markov Chains and Monte Carlo Methods参考文献 66被引用数 57

ひとこと要約

この論文では、自己相互作用スピン1/2のイジング模型に由来する共形ブートストラップにおける多項式行列計画問題（PMP）に特化した、オープンソースで並列化された任意精度の半定値計画法ソルバ―であるSDPBを紹介する。ブロック構造を活用し、内点法を最適化することで、3次元イジングCFTの臨界指数の計算を桁違いに高速化し、Δσ = 0.518151(6)、Δε = 1.41264(6)という新たな高精度の境界を達成した。

ABSTRACT

We introduce SDPB: an open-source, parallelized, arbitrary-precision semidefinite program solver, designed for the conformal bootstrap. SDPB significantly outperforms less specialized solvers and should enable many new computations. As an example application, we compute a new rigorous high-precision bound on operator dimensions in the 3d Ising CFT, $Δ_σ=0.518151(6)$, $Δ_ε=1.41264(6)$.

研究の動機と目的

• スピンを持つ複数の自己相関関数を含む共形ブートストラップ問題に特化した、高性能な半定値計画法ソルバ―の開発。
• 一般用途のソルバ―（例：SDPA-GMP）では、1回の3次元イジングCFT最適化に数週間を要するなど、著しい計算制限に直面する問題を克服すること。
• 保存則を伴う電流、ストレステンソル、混合自己相関関数を含む複雑なCFTにおける、新たな高精度ブートストラップ計算を可能にすること。
• アルゴリズムの特化と並列化が、数値的ブートストラップ研究における劇的な性能向上をもたらす可能性を示すこと。
• 今後の共形場理論および数値最適化分野の研究を促進するため、オープンソースで拡張可能なツール（SDPB）を提供すること。

提案手法

• 共形ブートストラップ問題を、多項式行列制約を持つ半定値計画問題（SDP）の構造的クラスである多項式行列計画問題（PMP）として定式化する。
• 自由変数yと行列変数Yを分離する変数分割法を用いて、PMPを標準的なSDPに変換する。
• ロバストな収束性と高速な超線形収束を実現するため、Mehrotraの予測子・修正子加速を用いた原双対内点法を実装する。
• PMPにおけるスチェル補行列のブロック構造を活用し、ブロックごとのコレスキー分解と行列・ベクトル演算を用いて、その計算を効率化する。
• PMP構造に特化した線形方程式系の解法を、コレスキー分解とLU分解を組み合わせて実装し、メモリと計算コストの最小化を図る。
• 大規模かつ高精度なブートストラップ問題を処理できるように、任意精度算術と複数CPUコアへの並列化を実装する。

実験結果

リサーチクエスチョン

• RQ1 カスタム半定値計画法ソルバ―は、一般用途のソルバーよりも共形ブートストラップ問題の解決において顕著な性能向上を達成できるか？
• RQ2 CFTブートストラップ計算に生じる大規模な多項式行列計画問題を解くために、最も効果的なアルゴリズム的最適化は何か？
• RQ3 並列化と任意精度算術は、3次元イジングCFTにおける臨界指数の境界の精度と速度をどの程度向上できるか？
• RQ4 高スピンオブザーバブルや複数の自己相関関数を含む場合、共形ブートストラップの計算複雑度はどのように変化するか？
• RQ5 SDPBのような特化型ソルバ―は、保存則を伴う電流や混合自己相関関数を含む、従来では不可能だった新たなブートストラップ計算を可能にするか？

主な発見

• SDPBにより、1回の3次元イジングCFTの妥当性計算の実行時間が、SDPA-GMPで約2週間かかっていたのを、1～3CPU時間（16コアマシンでは4～12分）に短縮した。
• 3次元イジングCFTにおいて、SDPBは新たな高精度境界 Δσ = 0.518151(6) および Δε = 1.41264(6) を計算し、精度の世界記録を更新した。
• SDPA-GMPに比べ100倍～1000倍の高速化を達成し、SDPBではΛ=19で1イテレーションあたり約2×10⁸回の乗算で済ませるのに対し、SDPA-GMPは約1×10¹¹回にのぼった。
• 性能向上の背景には、スチェル補行列のブロック単位での効率的計算と、PMP構造に最適化された線形代数ルーチンの実装がある。
• SDPBの任意精度算術により、収束の信頼性と高精度な境界の達成が可能となり、厳密なブートストラップ結果を求める上で不可欠である。
• ソルバ―はオープンソースであり、https://github.com/davidsd/sdpb で公開されており、再現性とコミュニティによる拡張が可能である。