Semiclassical spin-bath calculation of the nitrogen-isotope effect on NV$^-$ ensemble coherence in diamond

技術要約：窒素同位体効果によるNV⁻アンサンブル・コヒーレンスに関する半古典的スピンバス計算

問題提起
ダイヤモンド中の窒素空孔（NV⁻）アンサンブルは固体状態磁気センサとして機能し、その感度は、広帯域DCセンシングにおけるアンサンブル・デフェージング時間（$T_2^*$）および、ハーンエコー・コヒーレンス時間（$T_2$）によって決定される。窒素含有量が高いダイヤモンドでは、両方の時間は窒素置換体（P1）電子スピンのバスによって制限され、窒素濃度$[N]$に反比例してスケールする。天然の同位体存在比を持つダイヤモンド（$^{14}\text{N}$、核スピン $I=1$ が支配的）では、$T_2/T_2^*$ の比は約16であると観察されているが、最近の研究において、意図的に $^{15}\text{N}$ ドープ（$I=1/2$）されたアンサンブル磁気センサを用いた際、シングル量子慣習において、著しく低い $8.05 \pm 0.18$ という比が報告された。ここでの中心的な問いは、窒素核の同位体そのものがこのコヒーレンスの減少を説明できるかということであり、具体的には、核スピンの多重性がP1電子バスのダイナミクスを変化させるかどうかを調査することである。

手法
著者らは、Jahn-Teller（JT）分解されたハイパーファイン・テンソル（Parkら[5]より導出）を用いた、Bauchら[2]によって確立されたフレームワークを拡張した半古典的スピンバスモデルを採用している。このアプローチは、意図的に多体問題ではないものとして設計されており、将来のクラスター展開法（CCE）処理のベンチマークとなるよう、透明性の高い、パラメータフリーの計数論に基づいている。

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

• モデルの設定: ダイヤモンド格子内の原点にNV⁻センターを配置し、確率 $[N]$ に基づいてP1センターを置換サイトにランダムに割り当てる。
• ハイパーファイン分解: モデルはP1センターのJT歪み軸を組み込んでいる。有効ハイパーファイン分裂（$A_{\text{eff}}$）は、JT軸と磁場のなす角度に依存し、オン軸およびオフ軸の配向に対して異なる値をとる。
• フリップフロップ・ダイナミクス: P1–P1のスピン交換（フリップフロップ）は、二準位問題として扱われる。二つのP1センター間のハイパーファイン・エネルギー差（デチューニング $\delta$）がゼロのとき、ポピュレーション転送率は鋭くピークを示す。
• 共鳴チャネルの計数: バスの相関時間 $\tau_c$ は、ハイパーファイン的に縮退しているP1ペアの割合（$f_{\text{res}}$）によって支配される。著者らは、核スピンの多重性（$2I+1$）およびJT縮退に基づいて、この割合を解析的に計算している：
  • $^{14}\text{N}$ ($I=1$) の場合：$f_{\text{res}} = 1/4$。
  • $^{15}\text{N}$ ($I=1/2$) の場合：$f_{\text{res}} = 5/16$。
• シミュレーション: モデルは、窒素濃度範囲 $[N] = 0.1\text{--}100$ ppm におけるアンサンブル・コヒーレンス時間をシミュレートする。これは、双極子場強度（$\Delta_{\text{single}}$）と相関時間（$\tau_{c,\text{single}}$）の分布にわたって単一NVの減衰を平均化することで、アンサンブルの $T_2^*$ および $T_2$ を抽出する。

主な貢献と結果

1. $T_2^*$ の同位体独立性: シミュレーションは、$T_2^*$ が両方の同位体で同一であることを確認している。これは、$T_2^*$ がP1バスの準静的な双極子場（式3）によってのみ決定され、それがP1の位置に依存し、核スピンの状態には依存しないためである。シミュレーションされた傾き（$11.21 \pm 0.08$ $\mu\text{s}\cdot\text{ppm}$）は、約17%以内で実験値と一致している。
2. $T_2$ の同位体依存性: モデルは、あらゆる濃度において $T_2$ が $^{14}\text{N}$ よりも $^{15}\text{N}$ の方が短くなると予測している。$^{15}\text{N}$ における共鳴フリップフロップ・チャネルの増加（5/16 対 1/4）は、バスの相関時間を速める。
3. 定量的比率: モデルは、濃度に依存しない比率 $T_2(^{14}\text{N})/T_2(^{15}\text{N}) = 1.118 \pm 0.011$ を算出する。この値は、単位（1）に対して $11\sigma$ 上回っており、$^{15}\text{N}$ 濃縮が固定された窒素濃度において $T_2$ を本質的に短縮することを裏付けている。
4. 実験との乖離: モデルは、同位体効果の方向性と濃度独立性は再現しているものの、その大きさについては過小評価している。実験的な比率（Ref. [2] と Ref. [9] を組み合わせて導出）は $2.07 \pm 0.25$ であり、シミュレーション値の $1.118$よりも$3.8\sigma$ 大きい。モデルは、対数スケールにおける観察されたコントラストの約17%しか説明できていない。

意義と主張
本論文は、コヒーレンス減少に対する同位体寄与の下限を確立したと主張している。主な意義は、窒素核の同位体が、純粋に組合せ論的なメカニズム（P1のハイパーファイン縮退ペアの割合）を通じて、$T_2$（ただし $T_2^*$ は除く）を制御していることを示す点にある。

著者らは、以下の3つの具体的な知見を強調している：

• 構造的定数: 同位体比は、サンプル固有のパラメータではなく、核スピンの多重性とJT縮退に由来する、真の材料定数である。
• メカニズムの区別: 希薄な核バスにおいては、より高い核スピンが通常コヒーレンスを低下させるが、P1バスは逆の挙動を示す。ここでは、核は電子のフリップフロップのデチューナーとして機能する。すなわち、より大きなスピン多重性（$^{14}\text{N}$）は遷移周波数をより多くの値に分散させ、フリップフロップを抑制してバスを遅らせる。一方で、$^{15}\text{N}$ はより多くの共鳴チャネルを許容し、バスを加速させる。
• 多体物理学のベンチマーク: 半古典的な計数論（1.12）と実験的な観察（2.07）の間の乖離は、多体効果が重要であることを示唆している。著者らは、近接縮退チャネルや単純な速度方程式を超えた相関進化を考慮する、完全なCCE処理が、このギャップを埋めるために必要とされる次のステップであると断じている。

結論として、半古典的モデルは同位体効果の方向性と濃度独立性を予測することには成功しているが、観察されたコントラストの全容を解明するには、平均場的な確率論的取り扱いを超えた物理、具体的には相関のある多体ダイナミクスの導入が必要である。