Structure and information measures of few-electron systems under a spherically symmetric Gaussian potential within a density functional approach

技術的要約：ガウス型球対称ポテンシャル下における少電子系の構造と情報量

問題提起
量子ドット（QD）などのナノ構造における量子閉じ込めは、原子の電子および光学的特性を自由な状態と比較して著しく変化させる。これらの系を記述するために、様々なモデルポテンシャル（調和ポテンシャル、ウッドス・サクソン・ポテンシャルなど）が提案されてきたが、ガウス型ポテンシャル（GP）は、核反応、分子プロセス、およびフラーレンケージ内に閉じ込められた原子の有効な相互作用をモデル化するのに有用な、ソフトで連続的な境界を提供する。GPにおける二電子系のエネルギーに関する既存の文献はあるものの、多電子系、励起状態、および情報理論的尺度に関する包括的な研究は著しく不足している。本研究は、ガウス型球対称ポテンシャル内に閉じ込められたH、He様イオン（$Z=2-18$）、Li、およびBeの構造的および情報的特性を調査することで、これらの空白を埋めるものである。

手法
本研究では、任意の原子状態に適用可能な時間独立な密度汎関数理論（DFT）の枠組みを採用している。アプローチの中核は、以下の径方向のコーン・シャム（KS）方程式を解くことである：
$$\hat{H}(\vec{r})\psi_i(\vec{r}) = \epsilon_i(\vec{r})\psi_i(\vec{r})$$
ここで、有効ポテンシャルは、電子・核間の引力、古典的なクーロン斥力、交換相関（XC）ポテンシャル、および閉じ込めガウスポテンシャル $V_{conf}(r) = V_0 [1 - \exp(-r^2/R^2)]$ を含む。

主な手法構成要素は以下の通りである：

• 交換ポテンシャル： 物理的に動機付けられた、仕事関数に基づく交換ポテンシャルを使用し、交換エネルギーを電子とその関連するフェルミ・クーロン空孔電荷密度間の相互作用として扱う。
• 相関汎関数： 電子相関は、局所的なパラメータ化されたウィグナー汎関数と、非線形で勾配およびラプラシアンに依存するリー・ヤン・パー（LYP）汎関数の2つの異なる汎関数を用いて分析される。
• 数値解法： KS方程式を解くために、一般化擬スペクトル（GPS）法を採用している。この基底関数を用いない手法は、ディリクレ境界条件（$\psi(0) = \psi(r_c) = 0$）を満たすために、非一様で最適に分布された空間格子を利用しており、固有関数および固有値における高い精度を保証している。
• 情報量： 本研究では、$r$空間の密度およびそのフーリエ変換から導出される、位置空間のシャノンエントロピー（$S_r$）、運動量空間のシャノンエントロピー（$S_p$）、および位置空間のフィッシャー情報量（$I_r$）を算出する。相関による寄与は、XC汎関数と交換のみ（XX）の結果を比較することで推定される。

主な結果
調査範囲は、H、He、Li、Beの基底状態、およびHeの低励起状態（$1s2s$ $^1S$ および $^3S$）をカバーしている。閉じ込めパラメータ（深さ $V_0$ および幅 $R$）および原子番号（$Z$）の変化に関する主要な知見は以下の通りである：

1. エネルギー論：

   • エネルギーは、閉じ込め幅 $R$ が増加するにつれて単調に減少し、$R \to \infty$ において自由原子の極限に近づく。
   • $R \to 0$ のとき、運動エネルギーの増加によりエネルギーは正となる。
   • 固定された $V_0$ に対して、エネルギーが負から正へと遷移する点は、Heでは $R \approx 5$ であるが、LiおよびBeでははるかに小さい $R \approx 1$ である。
   • $Z$ が増加するにつれ、核引力（$-Z/r$）がガウス閉じ込めを支配する。その結果、高$Z$イオン（$Z > 8$）に対するGPの影響は減少し、系は中程度の閉じ込め条件下でも自由原子のように振る舞う。
   • 相関エネルギー（$|E_{corr}|$）は、中間領域（$V_0=5$ のとき $R \sim 0.5-0.6$ 付近）において極大を示し、これはGPと核ポテンシャルの競合が最も激しい領域である。

2. 情報理論的尺度：

   • シャノンエントロピー（$S_r, S_p$）： $S_r$ は中間領域の $R$ で極小を示し、最大の空間的局在化を示す一方、$S_p$ は対応する極大（運動量の非局在化）を示す。$Z$ が増加するにつれ、$S_r$ は減少し（局在化の強化）、$S_p$ は増加する。
   • フィッシャー情報量（$I_r$）： $I_r$ は中間領域の $R$ で極大に達し、電子密度の勾配の鋭さを反映する。このピークは、$Z$ が増加するにつれてより低い $R$ 値へとシフトする。
   • 相関効果： エネルギーおよび情報量への相関による寄与（$E_{corr}, S_{corr}, I_{corr}$）の大きさは、一般に低$Z$系で大きく、$Z$ が増加するにつれて減少する。ウィグナー汎関数は、一般にLYP汎関数よりも大きな相関の大きさを与える。

3. フィッシャー・シャノン平面およびコリンの予想：

   • フィッシャー・シャノン平面（$J_r$ 対 $I_r$）は、$Z$ の増加に伴う電子密度の漸進的な局在化と圧縮を明らかにしており、これは相対的な電子相関効果の弱体化を示している。
   • $R$ および $V_0$ の関数としての相関シャノンエントロピー（$S_{corr}$）と相関エネルギー（$E_{corr}$）の関係の分析は、非線形なループ状の特徴（コリンの予想）を明らかにした。このループ構造は、より高い $V_0$ においてより顕著であり、強固な閉じ込めと自由原子の振る舞いの間の系の遷移を反映している。

意義および主張
著者らは、本研究が、ガウス型ポテンシャル内に閉じ込められた原子およびイオンの基底状態および励起状態のエネルギー、ならびに情報エントロピー（$S_r, S_p, I_r$）を包括的に評価した初めての事例であることを主張している。これは、二電子系を超えてLi、Be、および励起He状態を含むものである。

本研究は、ドットのパラメータ（$V_0$ および $R$）を調整することで、エネルギーおよび情報量を操作できることを示している。結果は、閉じ込めが低$Z$系では特性を著しく変化させる一方で、高$Z$イオンでは核引力が支配的となり、閉じ込め効果が最小限になることを強調している。$S_{corr}$ 対 $E_{corr}$ 平面におけるループ状の特徴の特定は、従来のエネルギーに基づく記述を超えた、電子相関の多次元的な特徴付けを提供するものである。著者らは、これらの知見が、他のモデル閉じ込めポテンシャルを用いた将来の研究や、複雑性およびもつれ（エンタングルメント）の尺度への研究の拡張を動機付ける可能性があると示唆している。