Solid-angle based nearest-neighbor algorithm adapted for systems with low coordination number

原子や分子という目に見えない世界において、粒子の配列の仕方がその物質の性質を決定づける。ある物質が硬いダイヤモンドなのか、滑らかな潤滑剤なのか、あるいは流動的な液体であるかは、その構成要素の局所的な近傍関係に完全に依存している。これらの材料を理解するために、科学者たちはまず、「隣人とは誰か？」という、一見単純に見える問いに答えなければならない。粒子の密集した群衆の中では、どの粒子が接触しており、どの粒子が単に通り過ぎているだけなのかを判別することは必ずしも容易ではない。この区別は極めて重要である。なぜなら、粒子が持つ即時的な隣人の数（配位数として知られる）が、システム全体の構造を決定するからである。数十年にわたり、研究者たちはこれらの目に見えない境界線を描くために数学的なツールに頼ってきたが、群衆がまばらであったり、配列が特殊であったりする場合、これらのツールは苦戦し、遠くの粒子を親しい友人と誤認してしまうことがよくあった。

ユトレヒト大学とパリ・サクレー大学の研究チームは、この特定の問題を解決するための洗練された手法を開発した。彼らは、粒子の周囲の視界が仲間によってどの程度遮られているかをチェックすることで隣人を判定する、「立体角に基づく最近接アルゴリズム（solid-angle-based nearest-neighbor algorithm）」と呼ばれる既存の技術に着目した。この手法は、高密度で混沌とした系ではうまく機能するが、粒子が少なくまばらな、開いた構造を持つ格子状の構造においては、系統的な誤差を生じる傾向がある。このような低密度環境では、元のアルゴリズムは「手を広げすぎて」しまい、第1層の隣人に属するはずの粒子を、あたかも最初の円の一部であるかのようにカウントしてしまう。研究者たちは、計算に新しい調整可能な設定を追加することなく、この過剰カウントを修正するための幾何学的な補正を導入した。彼らがmSANNと呼ぶこの修正されたアプローチは、ハニカムパターンからダイヤモンド結晶に至るまで、複雑な構造において正しい隣人の数を特定することに成功し、微視的な世界のより正確な地図を提供している。

隣人を特定する上での核心的な課題は、「接触」を構成するものについての単一の普遍的な定義が存在しないことにある。完全な結晶においては答えは明白だが、実際の材料では熱エネルギーによって粒子が小刻みに震え（ジグル）、層の境界を曖昧にする。伝統的な手法は、多くの場合、固定された距離のカットオフ（閾値）に依存し、粒子の周囲に円を描いて、その中に入っている全員をカウントする。しかし、これは材料内で密度が変化する場合に失敗する。別の一般的な手法は、空間の幾何学的な分割を用い、各粒子の周囲を独自のセル（細胞）に分割する。これは任意の距離制限を回避できる一方で、微細な振動に敏感であり、ハニカム格子のような、各粒子がわずか3つの隣人しか持たない低配位構造において、遠くの粒子を誤って含んでしまうことがある。立体角法は、固定された距離の設定を必要としない、堅牢な代替案として設計された。これは、中心となる粒子の周りに球体を想定し、潜在的な隣人が占める角度の空間を計算することで機能する。アルゴリズムは、隣人たちが球体全体を埋め尽くすまで境界を広げていく。これは高密度の系では見事に機能するが、開いた構造においては、状況の幾何学がアルゴリズムを欺いてしまう。

研究者たちは、開いた構造において、元の手法が実質的に大きすぎる円を描いていることを発見した。3つの最近接隣人に囲まれた三角形の中心に粒子が位置している場面を想像してほしい。中心の粒子を囲む空間を満たすために、アルゴリズムはその三角形の角（頂点）を通る半径を計算する。その過程で、アルゴリズムは意図せずして、その三角形のすぐ外側、つまり構造の次の層に位置する粒子まで含んでしまう。これは、アルゴリズムが、隣人を単に「含む」ための円ではなく、隣人を「通過する」円によって空間を満たそうとするために起こる。この幾何学的な見落としにより、一貫した過大評価が生じ、第1層の隣人と第2層の隣人を混同してしまうのである。

これを修正するために、著者らは、隣人を通過する円と、それらが形成する形状の内側に収まる円との関係に基づいた、単純な幾何学的調整を提案した。彼らは、元の手法が外側の円を使用している一方で、これらの疎な構造に対してより正確なアプローチは、内側と外側の限界の間のどこかに位置する半径を使用することであると気づいた。彼らは、計算された半径をスケールダウンさせる修正を導入し、遠くの粒子を除外するのに十分なほど境界を縮小させつつ、実際の材料に見られる熱的な揺らぎも許容できるようにした。この調整は純粋に幾何学的なものであり、新しいパラメータやチューニングを必要とせず、元の手法の簡潔さを維持している。それは、開いた格子においてアルゴリズムが手を広げすぎるのを防ぐフィルターとして機能すると同時に、実在の材料に存在する熱的ノイズに対しても柔軟に対応できる。

チームは、新しい手法であるmSANNを、元のアルゴリズムおよび伝統的な幾何学的分割法を用いて、幅広いシミュレーション・システムに対してテストした。二次元のハニカム格子および正方格子のシミュレーションにおいて、元の手法は頻繁に隣人を誤認し、正しい数である3または4の代わりに、6または5とカウントすることが多かった。しかし、修正された手法は、すべての粒子に対して正確な配位数を一貫して特定し、理論的な構造と一致する鮮明で明確な分布を生み出した。ダイヤモンドやグラファイト構造を含む、低い配位数を持つ三次元のテストにおいても、元の手法は第1層と第2層を区別するのに再び苦戦した。mSANNの補正はこの問題を解決し、単純立方格子や体心立方格子を含む、テストされたすべての結晶タイプにおいて正しい隣人の数を特定した。

研究者たちはまた、異なる形状や隣人の数が混在する、クアジクリスタル（準結晶）のような、より複雑で無秩序なシステムについても調査した。これらの不均質な環境において、元の立体角法は、正方形の隙間の対角線方向に偽の接続を作り出し、事実上、別々の領域を結合させてしまうことがあった。修正されたアルゴリズムは、これらの偽のリンクを回避し、構造の真のトポロジーを保持した。さらに、結晶と無秩序な領域の境界のような、異なる相が共存するシステムにおいて、新手法は界面を横断して一貫した隣人の特定を提供したが、他の手法は顕著な不一致を示した。この堅牢性は、局所的な環境が地点ごとに大きく異なる、完全に秩序だったわけではない材料の研究において、この修正が特に価値があることを示唆している。

精度を超えて、研究者たちは計算の速度についても懸念していた。なぜなら、隣人の特定は多くの大規模シミュレーションにおける基本的なステップだからである。彼らは、現代のコンピューティングパワーを活用するようにアルゴリズムを実装し、並列処理を用いて計算を効率的に処理した。ベンチマークによれば、1,000個未満の粒子を持つ小さなシステムでは、伝統的な幾何学的手法が依然として最も速い選択肢である。しかし、システムサイズが大きくなるにつれ、修正されたアルゴリズムは大幅に高速化し、数百万個の粒子を含む非常に大きなシステムでは、伝統的な手法よりも2倍近く速い速度を実現した。この効率性と、低密度環境における改善された精度の組み合わせにより、この新手法は複雑な材料を分析するための強力なツールとなっている。

この研究は、問題の根底にある幾何学を注意深く見ることが、物理的な世界のモデリング方法にいかに大きな改善をもたらすかを示している。元の手法における「隣人」の定義が開いた構造において寛容すぎたことを認識することで、研究者たちは、数学的にエレガントでありながら実用的な効果も高い補正を導入することができた。修正されたアルゴリズムは、単に特定の誤りを修正するだけでなく、結晶の硬い格子から、無秩序な相の変動する配置に至るまで、物質の局所構造をマッピングするためのより信頼できる方法を提供するものである。原子スケールでの材料の挙動を研究する科学者にとって、その間の隙間に惑わされることなく、隣人の数を正確に数えることができるツールを持つことは、私たちの世界を構成する材料の性質を理解するための重要な一歩となる。