[論文レビュー] Third and fourth density and acoustic virial coefficients of neon from first-principles calculations

Robert Hellmann, Giovanni Garberoglio|arXiv (Cornell University)|Jan 29, 2026

Quantum, superfluid, helium dynamics被引用数 0

ひとこと要約

本研究は、相対論的およびポストブリオン–オッペンハイマー効果を含む高レベルのab initio二体、三体、四体ポテンシャルを用いたパス積分モンテカルロ法により、ネオンの三重および四重密度 virial係数と音響 virial係数を第一原理から算出する。

ABSTRACT

The third and fourth density and acoustic virial coefficients of neon were determined at temperatures between 10 and 5000 K from first principles employing the path-integral Monte Carlo (PIMC) approach. For these calculations, we used the pair potential of Hellmann $ extit{et al.}$ [J. Chem. Phys. 154, 164304 (2021)], which is based on supermolecular $ extit{ab initio}$ calculations with basis sets of up to octuple-zeta quality and levels of theory up to coupled cluster with single, double, triple, quadruple, and perturbative pentuple excitations [CCSDTQ(P)]. The potential also accounts for relativistic, retardation, and post-Born$-$Oppenheimer effects and is provided with reliable uncertainty estimates. To incorporate nonadditive interactions, we developed a nonadditive three-body potential based on extensive supermolecular CCSD(T), CCSDT, and CCSDT(Q) calculations with basis sets of up to sextuple-zeta quality. This potential also accounts for relativistic effects. The very small nonadditive four-body contributions to the fourth virial coefficients were considered using a relatively simple nonadditive four-body potential based on supermolecular CCSD(T) calculations. We calculated the third and fourth density and third acoustic virial coefficients directly by PIMC and the fourth acoustic virial coefficient indirectly using thermodynamic relations between the density and acoustic virial coefficients. The uncertainties of the pair potential and those estimated for our nonadditive three-body potential were rigorously propagated in the PIMC calculations into uncertainties for the virial coefficients. These uncertainties are distinctly smaller than those of almost all of the corresponding experimental virial coefficient data.

研究の動機と目的

• ガス計量と多体相互作用の基礎的検証のためのネオンの熱物理モデリングの正確性を動機づける。
• ネオンのab initio 二体、三体、四体ポテンシャルを高精度で開発・検証する。
• 広範な温度範囲（10–5000 K）にわたり、三重・四重密度 virial係数および三重・四重音響 virial係数を計算する。
• 電子構造の不確かさを virial係数へ伝搬し、実験データと比較する。

提案手法

• relativistic、遅延効果、およびポスト Born–Oppenheimer 効果を含む CCSDTQ(P) から構築された Hellmann らのab initio 二体ポテンシャルを使用。
• relativistic補正を含む広範な CCSD(T)、CCSDT、CCSDT(Q) 計算に基づく非加法的三体ポテンシャルを構築。
• CCSD(T) 計算からの簡略化ポテンシャルを用いて非常に小さな非加法的四体寄与を組み込む。
• 三重および四重密度 virial係数をパス積分モンテカルロ法で直接計算し、熱力学的関係から第四音響 virial係数を導出。
• ab initio の非加法ポテンシャル（三体・四体）に解析形を適合させ、評価と不確かさ伝搬を効率化。
• 非加法ポテンシャルの評価および±不確かさ変動（V123の±シフト ΔV123 および ΔV123±）を扱う Fortran 90 ルーチンを提供。

実験結果

リサーチクエスチョン

• RQ1 10–5000 K におけるネオンの三重・四重密度 virial係数とは何か。
• RQ2 同じ温度範囲に対するネオンの三重・四重音響 virial係数とは何か。
• RQ3 非加法的三体および四体相互作用はこれらの高次 virial係数にどの程度寄与し、不確かさはどうなるか。
• RQ4 第一原理の結果はネオンの virial係数の実験データとどのように比較されるか。
• RQ5 ab initio ポテンシャルの不確かさを virial係数へ厳密に伝搬できるか。

主な発見

• 研究は10–5000 Kにわたるネオンの三重・四重密度 virial係数を第一原理から算出して提供する。
• 第三音響 virial係数は直接決定され、第四音響 virial係数は熱力学的関係から間接的に得られる。
• 非加法的三体ポテンシャルは第三・第四 virial係数に有意に寄与する一方、四体項は非常に小さく、簡略化されたポテンシャルで扱われる。
• 二体ポテンシャルおよび非加法的三体ポテンシャルの不確かさは virial係数へ伝搬され、他の多くの実験データより小さい。
• 解析的非加法ポテンシャルはab initio の結果を良好な忠実度で再現し、±シフト（ΔV123 および ΔV123±）による不確かさ伝搬を可能にする。
• 比較から、第一原理の結果がネオンの理論データおよび実験データと整合するか、あるいは改善するかが示される。