Linear response across interaction regimes in two-dimensional ferromagnets

技術要約：二次元強磁性体における相互作用レジームを横断する線形応答

問題提起
原子層レベルの厚さを持つ二次元（2D）強磁性体の最近の発見は、交換相互作用の運動量依存性に起因するユニークな輸送現象を浮き彫りにした。これらの系では、マグノン・マグノン散乱の頂点が、入射運動量のスカラー積に比例する（$T_{ex} \propto -\mathbf{k} \cdot \mathbf{p}$）。その結果、低運動量のマグノンは弱く相互作用し、高運動量のマグノンは強く散乱するため、集団的な流体力学的挙動を示す可能性がある。このような系の輸送特性を理解するには、弾道的（ballistic）領域と流体力学的（hydrodynamic）領域を等しく扱うことができる微視的な枠組みが必要である。自然な枠組みは量子ボルツマン方程式（QBE）であるが、エネルギーおよび運動量保存則によって制約された多次元衝突積分を解くことは極めて困難である。既存のアプローチは、制御不能な近似（例：単一緩和時間近似）に頼るか、計算負荷の高い力まかせの離散化を行うことが多く、温度や磁場の全範囲にわたって、完全な周波数および運動量分解能を持つ線形応答 $\chi(\mathbf{q}, \omega)$ を正確に計算することは不可能であった。

手法
著者らは、平衡ボース・アインシュタイン分布関数をガウス関数の和として表現する（ガウス混合モデル、GMM）ことにより、線形化されたQBEを解く数値的手法を開発した。

1. ガウス表現： 平衡分布 $n^{(0)}_p$ を、少数のガウス関数の和（$n^{(0)}_p \approx \sum_s a_s e^{-p^2/(2\gamma_s^2)}$）として近似する。この表現は、約10個のガウス関数を用いて精度$10^{-10}$オーダーを実現しており、これにより多次元衝突積分を解析的に扱いやすくしている。
2. 線形化動力学理論： 微弱な外部ポテンシャル摂動に対する線形応答を定式化する。平衡からの偏差を角調和関数（angular harmonics）で展開することにより、2Dの積分微分方程式を、結合された1次元積分方程式の階層へと簡約化する。
3. 数値実装： 積分方程式を径方向の運動量グリッド上で離散化し、角調和関数で打ち切りを行うことで、問題を疎な線形代数系へと変換する。この手法により、広範な温度および磁場の範囲にわたって、衝突カーネルの行列要素を効率的に計算することが可能となる。
4. CrCl$_3$への適用： 本フレームワークを、放物線型の分散と交換駆動型の相互作用頂点を持つ2Dボゾン系としてモデル化した単層CrCl$_3$に適用する。このモデルは、特定の格子構造、交換結合（$J \approx 10.9$ K）、および異方性を考慮に入れている。

主な結果

1. 弾道から流体力学へのクロスオーバー： 本研究は、マグノン輸送における温度駆動のクロスオーバーを解明した。低温では、マグノンガスは希薄であり、弱く相互作用する低運動量モードが支配的であるため、弾道的輸送を示す。温度が上昇すると、高運動量モードが熱的に占有されるようになる。これらのモードは強く相互作用し、音波モード（sound mode）の伝播を特徴とする集団的な流体力学的レジムへとシステムを駆動する。
2. 完全な線形応答： 著者らは、完全な密度密度応答関数 $\chi(\mathbf{q}, \omega)$ を計算した。流体力学的極限において、応答は線形分散を持つ音波モードに対応する鋭いピークを示す。解析の結果、この音波モードは等温（isothermal）レジムではなく、等エントロピー（isentropic）レジムにおいて伝播することが明らかになった。
3. 量子ビットのデフェージング・シグネチャ： マグノンガスの磁気ノイズをプローブする窒素空孔（NV）センターのスピン量子ビットのデフェージング率（$1/T_2$）を計算した。結果は非単調な温度依存性を予測している：
   • 低温域： 温度とともにノイズが増加する（弾道的レジム）。
   • 高温域： 温度とともにノイズが減少する（流体力学的レジム）。
     この逆転現象は、交換駆動型マグノン粘性（$\mu_{ex} \sim 1/T$）の特異なスケーリングが、ノイズの公式における明示的な温度因子を相殺することに起因する。
4. 接触相互作用との比較： 著者らは、交換駆動型の流体と接触相互作用を持つボゾン流体を対比させている。接触流体も流体力学的な音波モードを示すが、高運動量の感受率は急速に減衰する。そのため、量子センサー（高運動量をフィルタリングするもの）は、接触流体においては非相互作用的な弾道的テールをプローブすることになり、結果として交換流体で見られるような異常な減少ではなく、温度に対する単調な増加を示す。

意義と主張
本論文は、量子ボルツマン動力学によって記述可能な、相互作用する2次元量子系における運動量および周波数分解能を持つ線形応答を計算するための一般的なフレームワークを確立したと主張している。

• 異常ノイズの解決： 本手法は、単層CrCl$_3$上のNVセンター実験で近年観察された異常な磁気ノイズ（温度とともに減少するノイズ）に対して、微視的な解釈を与えるものである。著者らは、これが交換駆動型マグノン流体力学の直接的な指紋であると主張している。
• 計算の実行可能性： GMM表現によって衝突積分のボトルネックを克服することで、制御不能な近似や計算不可能な力まかせの手法に頼ることなく、線形化されたQBEの正確な解出を可能にした。
• 汎用性： このフレームワークは強磁性体に限定されず、懸架されたグラフェンの屈曲フォノンや半導体中のボゾン・フェルミ混合系など、放物線型の分散とソフトな散乱位相空間を持つあらゆる2D系に適用可能である。著者らは、本手法がマルチモード系に拡張可能であり、将来的には非放物線型バンドにも一般化できる可能性があると述べている。

著者らは、予測されるデフェージング率はサンプルとプローブの距離（$d_{NV}^{-4}$）に敏感であることや、簡略化された理論が擬ギャップやマグノン質量の温度依存性を無視していることを挙げ、具体的な実験との定量的な比較については控えめな立場をとっている。しかし、定性的な一致と非単調な傾向の説明は、交換流体力学の図式を検証するものとして堅牢であると提示されている。