2025-11-23T14:37:16.213090

Heitler effect and resonance fluorescence in quantum magnonics

IlbuÄa, Dobrovitski, Blanter

We consider a coupled system of a qubit and a magnon mode in which the qubit is weakly driven. We demonstrate that the spectral steady-state responses of both the qubit and the magnon show, in addition to two sidebands split by the coupling, also a peak at the driving frequency with virtually zero linewidth. This phenomenon, which persists at both strong and weak coupling, is an analog of the Heitler effect in atomic physics, and shows the path towards building of coherent magnon sources.

academic

ハイトラー効果と量子マグノニクスにおける共鳴蛍光

基本情報

論文ID: 2501.01254
タイトル: Heitler effect and resonance fluorescence in quantum magnonics
著者: Enes Ilbuğa, V. V. Dobrovitski, Ya. M. Blanter
分類: cond-mat.mes-hall quant-ph
発表日: 2025年1月3日（arXivプレプリント）
論文リンク: https://arxiv.org/abs/2501.01254
著者所属: デルフト工科大学カブリナノサイエンス研究所、デルフト工科大学QuTech

要旨

本研究は、弱駆動を受ける量子ビットとマグノン模式が結合したシステムを考察している。量子ビットとマグノンのスペクトル定常応答は、結合分裂による2つのサイドバンドを示すだけでなく、駆動周波数においてほぼゼロ線幅のピークを示すことが明らかになった。この現象は強結合と弱結合の両方の場合に存在し、原子物理学におけるハイトラー効果の類似物であり、コヒーレントマグノン源の構築への道を示している。

研究背景と動機

問題背景

量子マグノニクスの発展需要：マグノニクスは磁激励（マグノン）の性質を研究する学問として、古典的および量子情報伝送において巨大な可能性を有している。機能化を実現するには、単一マグノンの制御された励起、操作、および検出が必要である。
材料の制限：現在主に使用されているイットリウム鉄ガーネット（YIG）材料は最高の磁気品質因子（10⁴を超える可能性がある）を有しているが、GHz周波数でのマグノン寿命は約100 nsであり、量子ビットと空洞の寿命よりはるかに短く、量子技術応用の重要な制限要因となっている。
既存方法の不足：従来の空洞マグノニクス方法は強結合を実現しているが、マグノン寿命は依然として制限要因であり、材料パラメータの改善の見通しは限定的である。

研究動機

本論文は、マグノニクスにおけるハイトラー効果の類似物を実現することにより、ほぼゼロ線幅の長寿命マグノンビームを生成し、材料寿命の制限を克服し、量子情報処理に新しい解決策を提供することを目指している。

核心的貢献

理論的発見：量子マグノニクスにおいてハイトラー効果の類似現象を初めて発見・分析し、弱駆動下の量子ビット-マグノン結合システムが駆動周波数においてゼロ線幅のスペクトルピークを生成することを示した。
普遍性の検証：この効果が強結合と弱結合の両方の場合に存在することを証明し、現象の適用範囲を拡張した。
コヒーレントマグノン源スキーム：高度に安定した、位相と周波数が固定された単一マグノン源を構築するための理論的基礎を提案した。
量子相関の確立：量子ビットとマグノン模式間に固定位相の長寿命安定相関を確立し、量子情報タスクのためのリソースを提供した。

方法の詳細

システムモデル

本研究は線形結合した量子ビット-マグノンシステムを考察し、回転座標系でのハミルトニアンは以下の通りである：

$\hat{H} = \frac{1}{2}\Delta_q\hat{\sigma}_z + \Delta_m\hat{m}^\dagger\hat{m} + g(\hat{\sigma}_+\hat{m} + \hat{\sigma}_-\hat{m}^\dagger) + \Omega(\hat{\sigma}_+ + \hat{\sigma}_-)$

ここで：

$\Delta_q = \omega_q - \omega_L$ 、 $\Delta_m = \omega_m - \omega_L$ はそれぞれ量子ビットとマグノンの周波数の駆動周波数に対する失調
$g$ は相互作用強度
$\Omega$ は駆動強度

マスター方程式記述

システムの進化は散逸を含むリンドブラッド主方程式で記述される：

$\frac{d}{dt}\hat{\rho}(t) = i[\hat{\rho}(t), \hat{H}] + \frac{\gamma}{2}(2\hat{\sigma}_-\hat{\rho}(t)\hat{\sigma}_+ - \{\hat{\sigma}_+\hat{\sigma}_-, \hat{\rho}(t)\}) + \frac{\kappa}{2}(2\hat{m}\hat{\rho}(t)\hat{m}^\dagger - \{\hat{m}^\dagger\hat{m}, \hat{\rho}(t)\})$

ここで $\gamma$ と $\kappa$ はそれぞれ量子ビットとマグノンの自発減衰率である。

スペクトル応答関数

システムのスペクトル定常応答は相関関数のフーリエ変換により定義される：

$S_q(\omega) = 2\text{Re}\int_0^\infty d\tau e^{-i(\omega-\omega_L)\tau}\langle\hat{\sigma}_+(\tau)\hat{\sigma}_-(0)\rangle_{ss}$

$S_m(\omega) = 2\text{Re}\int_0^\infty d\tau e^{-i(\omega-\omega_L)\tau}\langle\hat{m}^\dagger(\tau)\hat{m}(0)\rangle_{ss}$

弱駆動近似

実験設定

パラメータ範囲

結合強度： $g = 0.5$ から $g = 100$ （任意単位）
駆動強度： $\Omega = 0.01$ から $\Omega = 1$ （弱駆動条件）
減衰率： $\gamma = 5$ から $\gamma = 15$ 、 $\kappa = 0$ から $\kappa = 50$

数値方法

時間進化：主方程式の運動方程式系の数値積分
スペクトル計算：離散フーリエ変換を用いたスペクトル分布の計算
定常状態解：行列反転により定常状態演算子期待値を取得

分析方法

強結合極限： $g \gg \kappa, \gamma$
弱結合極限： $g \ll \kappa, \gamma$
共鳴条件： $\omega_q = \omega_m = \omega_L \equiv \omega_0$

実験結果

主要な発見

1. デルタ関数ピークの出現

量子ビットとマグノンのスペクトルの両方で観察される：

サイドバンドピーク： $\omega_0 \pm g$ に位置する有限線幅ローレンツピーク
弾性散乱ピーク：駆動周波数 $\omega_0$ に位置するデルタ関数ピーク、線幅はほぼゼロ

2. 強結合の場合

$g = 100, \gamma = 15, \kappa = 0$ の場合：

量子ビットとマグノンの散乱振幅はほぼ同一
デルタピークの高さは積分時間に線形に増加し、デルタ関数特性を確認

3. 弱結合の場合

$g = 0.5, \gamma = 10, \kappa = 50$ の場合：

マグノンは振幅は小さいもののデルタピークを示す
量子ビット散乱振幅はマグノン減衰に無関係

解析結果

弱駆動極限では、定常状態演算子期待値の解析表現は以下の通りである：

量子ビット： $|\langle\hat{\sigma}_+\rangle_{ss}|^2 = \frac{\Omega^2(\Delta^2 + \frac{1}{4}\kappa^2)}{g^4 + \Delta^4 + \frac{1}{4}\gamma^2\Delta^2 + \frac{1}{4}\kappa^2\Delta^2 - 2\Delta^2g^2 + \frac{1}{16}\kappa^2\gamma^2 + \frac{1}{2}g^2\gamma\kappa}$

マグノン： $|\langle\hat{m}^\dagger\rangle_{ss}|^2 = \frac{\Omega^2g^2}{g^4 + \Delta^4 + \frac{1}{4}\gamma^2\Delta^2 + \frac{1}{4}\kappa^2\Delta^2 - 2\Delta^2g^2 + \frac{1}{16}\kappa^2\gamma^2 + \frac{1}{2}g^2\gamma\kappa}$