Microwave-free vector magnetometry and crystal orientation determination with Nitrogen-Vacancy centers using Bayesian inference

技術要約：ベイズ推論を用いた窒素空孔中心によるマイクロ波フリーのベクトル磁気計測および結晶配向決定

問題提起
ダイヤモンド中の窒素空孔（NV）中心は、磁気計測、温度計測、電場検出のための高感度量子センサとして確立されている。従来のNV磁気計測は通常、スピン状態を操作するためにマイクロ波（MW）場を必要とする光検出磁気共鳴（ODMR）に依存している。この手法は堅牢ではあるものの、MWの導入は、環境への摂動、発熱、および、特に極低温条件下での統合を困難にするガルバニック接続の必要性といった、実用的な制限をもたらす。

マイクロ波を用いない代替手法も存在するが、例えば基底状態の準位反交差や、磁場によるスピン副準位の混合を利用するものなどである。しかし、これらの手法は、外部磁場を特定の結晶軸に厳密に整列させるか、あるいはフォトルミネセンス（PL）パターンの高解能な光学計測を必要とすることが多い。さらに、ベクトル磁気計測における決定的なボトルネックは、NV中心の正確な結晶配向を知る必要があることである。従来の校正は、制御された磁場下でのODMR分光によって行われるが、マイクロ波フリーの配向決定法は、撮像精度や双極子モデリングによって制限されることが多い。NV-NV間のクロスリラクゼーションを利用した最近の研究は、感度を示しているものの、事前に結晶学的な知識を必要とし、かつ共鳴を明確に割り当てるために精密な磁場整列を必要とするプロトコルによって制約を受けてきた。

手法
著者らは、2次元フォトルミネセンス（PL）マップを解析するためにベイズ推論を活用した、マイクロ波フリーのベクトル磁気計測および結晶配向決定のための一般的なフレームワークを提案している。この手法の核となるのは、クロスリラクゼーション現象である。すなわち、不等価なNV配向の間のスピン遷移が縮退するとき、フリップフロップ過程が強化され、PL強度の低下として観測される。

1. 物理モデル： システムは、固定された軸（実験室の$z$軸）に沿って印加される調整可能なバイアス磁場（$B_{bias}$）と、静的な外部磁場（$B_{ext}$）を利用する。ダイヤモンド試料は、この軸の周りで回転する（または、横方向の磁場が回転する）。試料フレームにおける全磁場は、回転角$\phi$とバイアス強度によって決定される。
2. 共鳴条件： クロスリラクゼーションは、磁場空間における特定の平面（例：$B^s_x = \pm B^s_y$のような磁場成分間の線形関係で定義される面）上で発生する。これらの条件は、PL信号における鋭いディップとして現れる。モデルは、NVの四面体対称性（$O$群）に固有の離散的な縮退を考慮しており、複数の物理的に異なる配向が同一のPLパターンを示すことを考慮に入れている。
3. 解析的順モデル： 著者らは、ハミルトニアンを繰り返し対角化する代わりに、期待されるPL信号の閉形式の解析的表現を導出した。このモデルは、共鳴条件を、外部磁場ベクトル、試料の配向行列、線幅、およびコントラストの重みによってパラメータ化されたローレンツ型のディップとして組み込んでいる。
4. ベイズ推論： 磁場ベクトルと結晶配向をPLマップから抽出するという逆問題は、ベイズ推定を用いて解決される。
   • 事後分布： 本手法は、未知のパラメータ（磁場成分、配向角）を表す$\theta$と、実験的なPLデータ$y$との間の完全な事後確率分布$P(\theta|y)$を計算する。
   • 縮退の処理： 単一の解を強制するピークフィッティング法とは異なり、ベイズアプローチは離散的な縮退を、事後分布における明確に異なるピークとして自然に捉え、不確実性と曖昧さを定量化する。
   • 対称性の削減： 計算効率を確保するため、NVの四面体対称性を利用して、配向のパラメータ空間を最小限の、対称的に不等価な領域に制限している。これにより、探索空間を24個の等価な構成から管理可能なサブセットへと削減している。

主な貢献

• アライメントフリーのフレームワーク： 本研究は、結晶配向の事前知識や、外部磁場を結晶軸に精密に整列させる必要のないプロトコルを導入している。磁場ベクトルとNV配向の両方が、データから同時に抽出される。
• 解析的効率性： クロスリラクゼーション共鳴の解析的モデルを導出することで、著者らはハミルトニアンの繰り返しの対角化という計算上のボトルネックを回避し、高速かつスケーラブルなベイズ推論を可能にした。
• 縮退への堅牢性： 本フレームワークは、NV系の離散的な対称性を明示的に扱い、単に最適な適合を選択したり破棄したりするのではなく、反転に関連する配向などの複数の有効な解を示す完全な事後分布を提供する。
• 実験的実証： 著者らは、広視野蛍光顕微鏡とダイヤモンドエンセルブルを用いて、本手法を実験的に検証した。彼らは、未知の整列を持つダイヤモンド結晶の配向を再構成することに成功し、その後、校正されたこの配向を用いて未知の外部磁場ベクトルを再構成した。

結果

• 配向決定： プルーフ・オブ・コンセプト実験において、チームは実験室フレームに対するダイヤモンド結晶の配向を決定した。実験的なPLマップのベイズ解析により、反転に関連する配向に対応する2つの高確率領域が得られ、これはシステムの対称性と一致していた。再構成されたパラメータ（$\alpha, \beta, \zeta$）は、実験的なPLパターンを高精度に再現した。
• ベクトル磁気計測： 先に決定された配向を用いて、著者らは未知の外部磁場を再構成した。磁場成分（$b_z, b_\perp, \phi_0$）の周辺事後分布は、よく局在化したピークを示した。再構成された値は、独立したODMR測定の結果と整合していた。
• 対称性解析： シミュレーションにより、回転軸が高対称な結晶方向（例：[100], [110], [111]）に一致する場合、PLマップは回転および鏡映対称性を示し、磁場再構成における曖昧さを生じることが示された。しかし、任意の回転軸を選択して対称性を破ることにより、一意の再構成が可能となる。

意義
本論文は、加熱や統合の制約により従来のODMR技術が困難な実用的なセンシングアプリケーションに適した、コンパクトでアライメントフリーなNV磁気計への一般的な経路を確立している。クロスリラクゼーションの特徴に対してベイズ推論を活用することで、本手法は、厳密な磁場整列や結晶学的な事前知識の必要性を排除し、任意の磁場構成におけるベクトル磁気計測および配向決定のための堅牢なソリューションを提供する。本研究は、不確実性を自然に定量化し、NV対称性に固有の離刻的な縮退を扱うフレームワークを提供することで、NVベースのセンシングの範囲を広げるものである。