[論文レビュー] Lattice QCD with light Wilson quarks

Martin Lüscher|ArXiv.org|Sep 27, 2005

Quantum Chromodynamics and Particle Interactions参考文献 26被引用数 10

ひとこと要約

本論文は、微小なクォーク質量および細かい格子定数（a ≤ 0.1 fm）における軽いウィルソンクォークを伴う格子QCDシミュレーションを効率的に行える、ドメイン分割プリコンディショニング付きハイブリッドモンテカルロ（HMC）アルゴリズムを提示する。この手法は、従来のアルゴリズムと比較して、自己相関時間とコスト係数（ν ≈ 13–53）が著しく短縮されており、大規模な格子（例：96×48³）においてパイオン質量が200 MeVまで低下した状態のシミュレーションが、一般のPCクラスタ上で可能になる。

ABSTRACT

Wilson's formulation of lattice QCD is attractive for many reasons, but perhaps mainly because of its simplicity and conceptual clarity. Numerical simulations of the Wilson theory (and of its improved versions) tend to be extremely demanding, however, to the extent that they rapidly become impractical at small quark masses. Recent advances in simulation algorithms now allow such simulations to be pushed to significantly smaller masses without having to compromise in other ways. Contact with chiral perturbation theory can thus be made and many physics questions can be addressed that previously appeared to be inaccessible.

研究の動機と目的

• 軽いクォークを伴う格子QCDシミュレーションにおける深刻なクリティカル・スローワンを克服すること。特に、微小な質量および細かい格子定数の条件下で。
• 従来のアルゴリズムが実用的でなくなる「クリティカル領域」におけるウィルソンフェルミオン形式の数値的シミュレーションを可能にすること。
• さまざまな格子サイズおよびクォーク質量にわたって良好な性能を維持する、スケーラブルで効率的なシミュレーションアルゴリズムの開発と検証すること。
• パイオン質量を微小に保ちながら格子の不正確さを制御できるようにすることで、手頃なPCクラスタ上でチャーミカル摂動論と接続できるようにすること。
• 標準的なPCクラスタ上でのO(a)補正付き作用素および複数の動的クォークを含む大規模シミュレーションの実現可能性を示すこと。

提案手法

• 空間的格子を重複する部分領域に分割することで、格子ディラック作用素をプリコンディショニングするドメイン分割アプローチを採用する。
• HMCアルゴリズムにおけるディラック方程式の解法を高速化するために、シュバルツ交替法手順をプリコンディショナーとして適用する。
• 特に微小なクォーク質量におけるHMCの効率を向上させるために、イーブン・オッズプリコンディショニングと2つの擬似フェルミオン場を組み合わせる。
• MPI（メッセージパッシングインターフェース）を用いてPCクラスタ上に実装し、並列処理を活用するとともに通信オーバーヘッドを低減する。
• 非摂動的O(a)補正として、調整されたcsw係数を用いたSheikholeslami–Wohlert項を統合し、離散化誤差を低減する。
• アルゴリズムの効率とスケーリング特性を定量化するために、プラケットの統合自己相関時間とコスト係数νを測定する。

実験結果

リサーチクエスチョン

• RQ1 ドメイン分割プリコンディショニングは、軽いウィルソンクォークを伴う格子QCDシミュレーションの計算コストを著しく低減できるか？
• RQ2 新しいHMCアルゴリズムは、特にクリティカル領域においてクォーク質量および格子サイズに伴う良好なスケーリングを維持するか？
• RQ3 パイオン質量が小さい（例：mπ ≈ 200 MeV）および格子定数が細かい（a ≤ 0.1 fm）条件下でも、効率的なサンプリングが可能か？
• RQ4 自己相関とコストの観点から、ドメイン分割HMCは、イーブン・オッズプリコンディショニングを施した標準HMCと比較して優れた性能を示すか？
• RQ5 O(a)補正を組み込んだ場合にアルゴリズムは頑健であるか？また、ストレンジ海クォークの導入も可能か？

主な発見

• ドメイン分割HMCアルゴリズムにより、さまざまな格子サイズおよびクォーク質量においてコスト係数νが13～53の範囲にまで低下しており、従来のシミュレーション（例：mπ ≈ 255 MeVでν ≈ 29）と比較して顕著に低減されている。
• アルゴリズムはクォーク質量および格子サイズに伴う良好なスケーリングを示しており、微小な質量や大規模な格子でもνの著しい増加は見られない。
• a ≈ 0.06–0.08 fm、パイオン質量が200 MeVまで低下した64×32³および96×48³格子上でのシミュレーションが、一般のPCクラスタ上で現在可能である。
• コスト係数νは、さまざまな格子サイズおよびクォーク質量で安定しており、頑健性と最小限のクリティカル・スローダウンを示している。
• O(a)補正を組み込んだ場合でもアルゴリズムは良好に動作し、Sheikholeslami–Wohlert項を含んでも顕著な性能劣化は観察されない。
• 32ノードクラスタ上では1日あたり1000件を超える有効なゲージ配置が受理されており、標準的なハードウェア上での実用的効率の高さが示されている。