Cloud parameter estimation for interacting BEC after time-of-flight

技術要約：飛行時間（Time-of-Flight）後の相互作用するBECにおけるクラウド・パラメータ推定

問題提起
極低温原子物理学において、飛行時間（TOF）展開後の吸収イメージングは、ボース・アインシュタイン凝縮体（BEC）から温度（$T$）、原子数（$N$）、および凝縮率といったシステム・パラメータを抽出するための標準的な手法である。純粋な熱的クラウド（ボース増幅を伴う理想気体としてモデル化される）および純粋なBEC（グロス・ピタエフスキー方程式およびトーマス・フェルミ近似によってモデル化される）の自由膨張については解析的なモデルが存在するが、これらのモデルは部分的に凝縮したクラウドに対しては機能しない。具体的には、標準的な解析的手法は、熱的成分と凝縮成分が互いにデカップリングしていることを前提としている。しかし、実際には、高密度の凝縮原子が熱的原子に対して斥力的な平均場力を及ぼし、熱的クラウドの膨張プロファイルを変化させている。相互作用する部分凝縮クラウドの自由膨張に関する解析解は存在しないため、研究者は通常、非相互作用のフィッティングに依存しており、これが系統的な不確実性を導入している。

手法
著者らは、部分凝縮クラウドの膨張をモデル化するために、弾道的近似に基づいた数値シミュレーションを用いている。その手法は以下の通りである：

1. 平均場ポテンシャル： シミュレーションには、トラップ内におけるBECと熱的原子の間の斥力相互作用が組み込まれている。熱的原子が経験するポテンシャルは、$U_0$ を相互作用エネルギー、$n_0$ をBEC密度として、$2U_0n_0$ という項を含むように修正されている。また、熱的クラウドの密度が低いため、BEC自体は熱的クラウドの影響を受けないものとして扱われている。
2. 数値積分： 非相互作用のクラウドに対する解析解とは異なり、有効ポテンシャル $V_{eff}$ にBEC密度項が含まれる場合、熱的密度分布 $n_{th}(r, t)$ の積分は解析的に解くことができない。著者らは、初期グリッド（トラップ内のクラウドを表す）と膨張グリッド（TOFプロファイルを表す）を総和をとることにより、この分布を数値的に計算している。各膨張グリッドの点について、必要な運動量を備えた特定の初期グリッドの点から原子が由来する確率を算出している。
3. 誤差推定プロトコル： 相互作用を無視することによって導入される誤差を定量化するために、既知の入力パラメータ（$N_{s}, T_{s}$）を用いたシミュレーションを用いて合成TOF画像を作成している。これらの合成画像に対し、セクション2.3で述べられている標準的なボース増幅分布（式7）を用いたフィッティングを行う。シミュレーションによる入力値とフィッティング値の差から、熱的原子数（$\Delta N_{th}$）および温度（$\Delta T$）の相対誤差を算出する。
4. 実験的検証： シミュレーションは、$^{87}$Rb系からの実験データに対するフィッティング関数としても使用されている。実験画像のシミュレーションによるフィッティング結果を、標準的なボース増幅関数および理論的なセミ・イデアル・モデルを用いたフィッティングと比較している。

主要な結果
本研究は、相互作用を無視することが、長い飛行時間においてもパラメータ推定において重大かつ無視できない誤差をもたらすことを明らかにしている：

• プロファイルの変形： BECからの斥力により、短時間のTOFでは熱的クラウドの密度中心に「ディップ（窪み）」が生じる。クラウドが膨張するにつれて、熱的原子がこの穴を埋めていき、最終的には、長時間のTOFにおいて標準的なボース増幅の予測よりもピークの鋭い分布へと至る。
• パラメータ誤差：
  • 原子数： 非相互作用モデルによるフィッティングは、一般に、短時間のTOFにおいて熱的原子数を過大評価する（中心のディップが低密度として解釈されるため）ことがあり、長時間のTOFでは過小評価する場合がある。誤差は短時間のTOFで最大約64%に達し、長時間のTOFでは約-5%に収束する。
  • 温度： 温度誤差は一般に $\pm 10\%$ 以内である。誤差の符号は、TOFの期間とトラップの形状に依存する。例えば、球状トラップにおける短時間のTOFでは温度は過大評価され、長時間のTOFでは過小評価される。
• 幾何学的形状への依存性： 誤差の大きさと符号は、トラップの幾何学的形状（シガー型、パンケーキ型、または球状）およびアスペクト比に大きく依存する。球状トラップは温度推定において最適（温度誤差を最小化）であるが、原子数の推定においては最悪の結果となる。逆に、非対称なトラップ（シガー型やパンケーキ型）は、原子数の決定には適しているが、より大きな温度誤差を導入する。
• 実験との一致： 実験データに適用した場合、シミュレーションから導出された分布を用いたフィッティングは、標準的なボース増幅フィットよりも、セミ・イデアル・モデルとよく一致する凝縮率を与える。標準的なフィットは、温度を過小評価し、BEC分率をわずかに過大評価する傾向がある。

意義と主張
本論文は、凝縮原子と熱的原子の間の斥力相互作用が、膨張中の熱的分布を根本的に再形成するため、標準的な解析的フィットでは高精度な測定には不十分であることを主張している。

• 系統的誤差の定量化： 著者らは、現在の実験手法に内在する誤差を系統的に定量化しており、これらの誤差は、BECが熱的原子の運動量分布を変化させるため、長い膨張時間においても持続することを実証している。
• 改良されたフィッティングモデル： 本研究の主要な貢献は、数値シミュレーションをフィッティング関数として用いることで、より正確なパラメータが得られることを示した点にある。このアプローチは、理想気体の仮定と、複雑な相互作用を持つクラウドという現実との間の溝を埋めるものである。
• トレードオフ： 本研究は、簡略化されたモデルを使用する場合、トラップの幾何学的形状とTOFの時間の選択によって、温度の正確な決定と原子数の正確な決定のどちらを優先するかが決まるというトレードオフを浮き彫りにしている。
• 限界： 著者らは、自身のシミュレーションが弾道的膨張を仮定していること（TOF中の相互作用を無視していること）を控えめに述べている。このモデルが標準的なフィットよりもセミ・イデアル・モデルとの一致を改善していることは認めつつも、飛行中の相互作用を考慮しないことによる小さな不一致が残る可能性があり、それには根本的に異なる理論的アプローチが必要であることを認めている。

結論として、本論文は、部分凝縮クラウドを含む実験において、非相互作用の解析的フィットに依存することは系統的なバイアスを導入することを確立している。提案された数値的フィッティング手法は、相互作用する膨張問題に対する完全な解析解を必要とせずに、クラウド・パラメータをより高精度に抽出するための道筋を提供するものである。