第3章：カイラルフォノンの実験的検出

ビデオ講義

このビデオは以下のテキストと同じ内容をカバーしています。お好みの学習形式をお選びください。

学習目標

• 円偏光ラマン分光法の実験セットアップと測定原理を理解する
• カイラルフォノンのラマン選択則を導出・適用できる
• σ⁺/σ⁺とσ⁺/σ⁻配置の物理的意味を説明する
• フォノン吸収の円二色性とフォノン角運動量の関係を理解する
• ヘリシティ分解フォトルミネッセンスによるバレー偏極検出の原理を学ぶ
• 超高速分光法によるカイラルフォノンのダイナミクス測定を習得する
• 非弾性中性子・X線散乱の可能性と課題を評価する
• 主要な実験結果（WSe₂、α-石英）を解釈できる
• Pythonで偏光依存性をシミュレーションできる

3.1 円偏光ラマン分光法

3.1.1 実験セットアップ

円偏光ラマン分光法は、カイラルフォノンを検出する最も強力で広く使われている実験技術です。この手法の鍵は、入射光と散乱光の両方の偏光状態を制御することにあります。

graph LR
    A[レーザー] --> B[直線偏光子]
    B --> C[QWP1: σ⁺/σ⁻選択]
    C --> D[対物レンズ]
    D --> E[試料]
    E --> F[対物レンズ]
    F --> G[QWP2: 分析]
    G --> H[分析用偏光子]
    H --> I[分光器]
    I --> J[CCD検出器]

3.1.2 1/4波長板の動作原理

1/4波長板（QWP）は、直線偏光を円偏光に変換する複屈折光学素子です。その動作は以下のように理解できます：

直線偏光から円偏光へ

QWPの速軸（fast axis）に対して45°で入射する直線偏光を考えます。電場は2つの成分に分解されます：

E_入射 = E₀(cos 45° x̂ + sin 45° ŷ) = (E₀/√2)(x̂ + ŷ)

QWP内部では、2つの偏光成分が異なる速度で伝播し、位相差δ = π/2（90°）が生じます：

E_出射 = (E₀/√2)(x̂ + e^(iπ/2) ŷ) = (E₀/√2)(x̂ + i ŷ)

これは右円偏光（σ⁺）を表します。QWPを90°回転させると左円偏光（σ⁻）が得られます：

E_出射 = (E₀/√2)(x̂ - i ŷ)  → σ⁻

3.1.3 σ⁺/σ⁺とσ⁺/σ⁻配置

カイラルフォノンの測定では、2つの標準的な配置が使用されます：

3.1.4 カイラルフォノンのラマン選択則

カイラルフォノンのラマン散乱選択則は、角運動量保存則から導出されます。光子とフォノンの角運動量を考えます：

角運動量収支

L_入射光 + L_フォノン = L_散乱光 + L_電子励起

バレー励起を無視すると（または基底状態での測定）：

L_入射光 - L_散乱光 = L_フォノン

各配置での選択則

3.1.5 円偏光度の測定

円偏光度（degree of circular polarization, DCP）は、カイラルフォノンの強さを定量化する重要な指標です：

DCP = (I_σ⁺/σ⁺ - I_σ⁺/σ⁻) / (I_σ⁺/σ⁺ + I_σ⁺/σ⁻)

ここで、I_σ⁺/σ⁺とI_σ⁺/σ⁻は、それぞれの配置でのラマン散乱強度です。

DCPの物理的解釈

• DCP = +1：完全な右円偏光（純粋な右カイラルフォノン）
• DCP = -1：完全な左円偏光（純粋な左カイラルフォノン）
• DCP = 0：非偏光またはアキラル
• |DCP| < 1：部分的なカイラリティ（混合状態）

3.2 フォノン吸収の円二色性

3.2.1 赤外円二色性（IR-CD）

赤外円二色性は、物質が左円偏光と右円偏光の赤外光を異なる程度に吸収する現象です。カイラルフォノンを持つ材料では、フォノン共鳴周波数でIR-CDが観測されます。

測定原理

円二色性シグナルは、2つの円偏光状態での吸収係数の差として定義されます：

ΔA(ω) = A_σ⁺(ω) - A_σ⁻(ω)

または規格化された形式：

CD(ω) = (A_σ⁺(ω) - A_σ⁻(ω)) / (A_σ⁺(ω) + A_σ⁻(ω))

3.2.2 フォノン角運動量との関連

IR-CDとフォノン角運動量の関係は、電磁場と格子の相互作用ハミルトニアンから導出できます。カイラルフォノンモードνでは：

ΔA_ν ∝ Im[ε_xy(ω_ν)] ∝ L_ph,ν

ここで、ε_xy(ω)は誘電率テンソルの非対角成分で、フォノン角運動量L_ph,νに比例します。

3.2.3 IR-CDとラマン散乱の相補性

3.3 ヘリシティ分解フォトルミネッセンス

3.3.1 バレー偏極検出の原理

2次元TMDにおいて、ヘリシティ分解フォトルミネッセンス（HR-PL）は、バレー物理とカイラルフォノンの結合を調べる強力な手法です。

測定手順

1. バレー選択励起：σ⁺（またはσ⁻）円偏光でK（またはK'）谷を選択的に励起
2. フォノン支援過程：励起子-フォノン相互作用が谷間散乱を引き起こす可能性
3. 偏光分解検出：σ⁺とσ⁻成分のPL強度を分離測定
4. バレー偏極度計算：P_v = (I_σ⁺ - I_σ⁻)/(I_σ⁺ + I_σ⁻)

3.3.2 フォノン支援過程

カイラルフォノンは、谷間散乱において重要な役割を果たします：

|K, σ⁺⟩ + カイラルフォノン(L_ph = -2ℏ) → |K', σ⁻⟩

この過程は角運動量保存を満たします：

L_初期 = +ℏ (K谷) + (-2ℏ) (フォノン) = -ℏ (K'谷) = L_最終

3.4 超高速分光法

3.4.1 カイラルフォノンの時間分解検出

超高速分光法は、カイラルフォノンのダイナミクスを実時間で観測できる強力な技術です。典型的なセットアップはポンプ-プローブ配置を使用します。

実験配置

1. ポンプパルス：フェムト秒（10⁻¹⁵秒）の円偏光パルスで試料を励起
2. 遅延制御：機械的遅延ステージでプローブパルスのタイミングを制御
3. プローブパルス：弱い円偏光パルスで時間発展を観測
4. ヘテロダイン/ホモダイン検出：微弱な信号変化を高感度検出

3.4.2 コヒーレントフォノン生成

コヒーレントフォノンは、パルス励起によって生成される位相が揃ったフォノン振動です。カイラルフォノンの場合、円偏光ポンプが重要です：

誘導ラマン散乱機構（ISRS）

強いポンプパルスと弱いプローブパルスが仮想中間状態を経由してフォノンを生成：

ポンプ(σ⁺) + ポンプ(σ⁺) → カイラルフォノン(L = 0)
ポンプ(σ⁺) + ポンプ(σ⁻) → カイラルフォノン(L = ±2ℏ)

3.4.3 ポンプ-プローブ技術

差分反射率の時間発展は以下のように表されます：

ΔR/R(t, ω_probe) = A₀ · cos(ω_ph t + φ₀) · exp(-t/τ_ph) + 非振動成分

ここで：

• A₀：振幅（電子-フォノンカップリングに比例）
• ω_ph：フォノン角周波数
• φ₀：初期位相
• τ_ph：減衰時間

3.5 非弾性中性子・X線散乱

3.5.1 偏極中性子散乱

非弾性中性子散乱（INS）は、フォノン分散を測定する古典的な手法ですが、カイラルフォノンの検出には課題があります。

原理

中性子がフォノンと相互作用してエネルギーと運動量を交換：

ℏω = E_入射 - E_散乱
ℏq = k_入射 - k_散乱

偏極中性子の利点

スピン偏極した中性子ビームを使用すると、磁気相互作用を通じてフォノンの角運動量成分に感度を持つ可能性があります：

• 中性子スピン：s = 1/2ℏ
• フォノン角運動量：L_ph = ±2ℏ
• 角運動量移行：ΔL = L_ph - Δs

3.5.2 共鳴非弾性X線散乱（RIXS）

RIXSは、円偏光X線を使用してフォノン励起を測定する新興技術です：

利点

• 高いエネルギー分解能（~100 meV → 数meV）
• 小さな試料サイズで測定可能（μmオーダー）
• 元素選択性（特定の原子のフォノンに焦点）
• 運動量分解能

カイラルフォノン検出への応用

円偏光軟X線RIXSは、原理的にカイラルフォノンの角運動量を検出できます：

入射X線(σ⁺) → 内殻励起 → フォノン放出 → 散乱X線(σ⁺ or σ⁻)

しかし、実用化にはX線光学系の改良と検出器の高感度化が必要です。

3.6 主要な実験結果

3.6.1 WSe₂実験（2018–2019）

単層WSe₂におけるカイラルフォノンの観測は2018年の報告（Zhu et al., Science）に端を発し、 後続研究（例：Chen et al., 2019, Nat. Phys.）でバレー–フォノン結合と角運動量の関係が精査されました。

実験手法

• 試料：CVD成長単層WSe₂、SiO₂/Si基板上
• レーザー：532 nm、出力 < 1 mW（過熱防止）
• 温度：77 K（液体窒素冷却）
• 配置：σ⁺/σ⁺、σ⁺/σ⁻、σ⁻/σ⁺、σ⁻/σ⁻の4配置で測定

主要な発見

1. E'モード（~250 cm⁻¹）のカイラリティ：
   • σ⁺/σ⁺配置：強いピーク
   • σ⁺/σ⁻配置：ほぼゼロ
   • σ⁻/σ⁺配置：ほぼゼロ
   • σ⁻/σ⁻配置：強いピーク
   • → 完璧な谷選択性（DCP ≈ 1）
2. A₁モード（~250 cm⁻¹）の非カイラル性：
   • すべての配置で同じ強度
   • DCP ≈ 0
3. 温度依存性：
   • 室温ではE'モードのDCPが減少（~0.3）
   • フォノン-フォノン散乱による谷間混合

3.6.2 α-石英実験（Zhu et al. 2018）

Zhuらは、3次元カイラル結晶α-石英でフォノン角運動量を光学的に検出しました。

実験アプローチ

• 試料：高品質α-石英単結晶、c軸方向測定
• 技術：フーリエ変換赤外（FTIR）円二色性分光法
• 温度範囲：10 K - 300 K
• 周波数範囲：400 - 1400 cm⁻¹

主要な結果

1. 強いIR-CDシグナル：
   • ~1100 cm⁻¹（Si-O伸縮モード）で最大
   • CD ≈ 5%（α-石英の螺旋対称性を反映）
2. フォノン角運動量の定量化：
   • 第一原理計算と比較
   • L_ph ≈ 0.3ℏ（格子単位胞あたり）
3. エナンチオマー依存性：
   • 左旋性石英：正のCD
   • 右旋性石英：負のCD
   • → 構造カイラリティとフォノンカイラリティの一対一対応

3.7 Pythonコード：ラマン偏光依存性シミュレーション

3.7.1 カイラルフォノンのラマン強度計算

import numpy as np
import matplotlib.pyplot as plt

# ラマン選択則に基づくカイラルフォノンのシミュレーション

class ChiralPhononRaman:
    """カイラルフォノンのラマン散乱強度を計算するクラス"""

    def __init__(self, omega_phonon, gamma, chirality='right'):
        """
        パラメータ:
        -----------
        omega_phonon : float
            フォノン周波数 (cm⁻¹)
        gamma : float
            減衰定数 (cm⁻¹)
        chirality : str
            'right' (L = +2ℏ) または 'left' (L = -2ℏ)
        """
        self.omega = omega_phonon
        self.gamma = gamma
        self.chirality = chirality

    def raman_tensor(self):
        """
        カイラルフォノンのラマンテンソルを返す

        E'モード（二重縮退）の場合:
        右カイラル: R ∝ [[1, i], [-i, 1]]
        左カイラル: R ∝ [[1, -i], [i, 1]]
        """
        if self.chirality == 'right':
            # 右カイラル (K谷)
            R = np.array([[1, 1j], [-1j, 1]], dtype=complex)
        else:
            # 左カイラル (K'谷)
            R = np.array([[1, -1j], [1j, 1]], dtype=complex)
        return R

    def intensity(self, omega_laser, polarization_in, polarization_out):
        """
        ラマン散乱強度を計算

        パラメータ:
        -----------
        omega_laser : float
            レーザー周波数 (cm⁻¹)
        polarization_in : str
            入射偏光 ('sigma+', 'sigma-', 'linear_x', 'linear_y')
        polarization_out : str
            散乱偏光 ('sigma+', 'sigma-', 'linear_x', 'linear_y')

        戻り値:
        -------
        intensity : float
            ラマン散乱強度
        """
        # 偏光ベクトルの定義
        pol_vectors = {
            'sigma+': np.array([1, 1j]) / np.sqrt(2),
            'sigma-': np.array([1, -1j]) / np.sqrt(2),
            'linear_x': np.array([1, 0]),
            'linear_y': np.array([0, 1])
        }

        e_in = pol_vectors[polarization_in]
        e_out = pol_vectors[polarization_out].conj()

        # ラマンテンソル
        R = self.raman_tensor()

        # ラマン散乱振幅: e_out · R · e_in
        amplitude = np.dot(e_out, np.dot(R, e_in))

        # ローレンツ型共鳴
        omega_stokes = omega_laser - self.omega
        resonance = 1 / ((omega_stokes - self.omega)**2 + self.gamma**2)

        # 強度 ∝ |振幅|² × 共鳴因子
        intensity = np.abs(amplitude)**2 * resonance

        return intensity

    def spectrum(self, omega_laser, omega_range, polarization_in, polarization_out):
        """
        ラマンスペクトルを計算（ストークスシフトの関数）

        パラメータ:
        -----------
        omega_laser : float
            レーザー周波数 (cm⁻¹)
        omega_range : array
            ラマンシフト範囲 (cm⁻¹)
        polarization_in, polarization_out : str
            偏光配置

        戻り値:
        -------
        spectrum : array
            ラマン強度スペクトル
        """
        spectrum = np.zeros_like(omega_range)

        for i, omega_shift in enumerate(omega_range):
            # ローレンツ型ピーク
            spectrum[i] = 1 / ((omega_shift - self.omega)**2 + self.gamma**2)

        # 偏光依存性の振幅因子
        pol_vectors = {
            'sigma+': np.array([1, 1j]) / np.sqrt(2),
            'sigma-': np.array([1, -1j]) / np.sqrt(2),
            'linear_x': np.array([1, 0]),
            'linear_y': np.array([0, 1])
        }

        e_in = pol_vectors[polarization_in]
        e_out = pol_vectors[polarization_out].conj()
        R = self.raman_tensor()

        amplitude = np.dot(e_out, np.dot(R, e_in))
        intensity_factor = np.abs(amplitude)**2

        return spectrum * intensity_factor


# シミュレーション実行
# WSe₂のE'モードをモデル化
omega_E_prime = 250  # cm⁻¹
gamma = 2  # cm⁻¹
omega_laser = 18797  # 532 nm → cm⁻¹

# ラマンシフト範囲
omega_shift = np.linspace(200, 300, 500)

# 右カイラルフォノン（K谷）
phonon_K = ChiralPhononRaman(omega_E_prime, gamma, chirality='right')

# 左カイラルフォノン（K'谷）
phonon_Kp = ChiralPhononRaman(omega_E_prime, gamma, chirality='left')

# 4つの偏光配置でスペクトルを計算
configs = [
    ('sigma+', 'sigma+', 'σ⁺/σ⁺'),
    ('sigma+', 'sigma-', 'σ⁺/σ⁻'),
    ('sigma-', 'sigma+', 'σ⁻/σ⁺'),
    ('sigma-', 'sigma-', 'σ⁻/σ⁻')
]

fig, axes = plt.subplots(2, 2, figsize=(14, 10))
axes = axes.flatten()

for i, (pol_in, pol_out, label) in enumerate(configs):
    # K谷（右カイラル）とK'谷（左カイラル）の寄与
    spec_K = phonon_K.spectrum(omega_laser, omega_shift, pol_in, pol_out)
    spec_Kp = phonon_Kp.spectrum(omega_laser, omega_shift, pol_in, pol_out)

    # プロット
    axes[i].plot(omega_shift, spec_K, 'b-', linewidth=2, label='K谷 (右カイラル)')
    axes[i].plot(omega_shift, spec_Kp, 'r--', linewidth=2, label="K'谷 (左カイラル)")
    axes[i].fill_between(omega_shift, 0, spec_K, alpha=0.3, color='blue')
    axes[i].fill_between(omega_shift, 0, spec_Kp, alpha=0.3, color='red')

    axes[i].set_xlabel('ラマンシフト (cm⁻¹)', fontsize=11)
    axes[i].set_ylabel('ラマン強度 (任意単位)', fontsize=11)
    axes[i].set_title(f'配置: {label}', fontsize=12, fontweight='bold')
    axes[i].legend(fontsize=10)
    axes[i].grid(True, alpha=0.3)
    axes[i].set_xlim(200, 300)

plt.tight_layout()
plt.savefig('chiral_phonon_raman_polarization.png', dpi=300, bbox_inches='tight')
plt.show()

# 円偏光度（DCP）の計算
I_pp = phonon_K.spectrum(omega_laser, np.array([omega_E_prime]), 'sigma+', 'sigma+')[0]
I_pm = phonon_K.spectrum(omega_laser, np.array([omega_E_prime]), 'sigma+', 'sigma-')[0]
DCP_K = (I_pp - I_pm) / (I_pp + I_pm)

I_mm = phonon_Kp.spectrum(omega_laser, np.array([omega_E_prime]), 'sigma-', 'sigma-')[0]
I_mp = phonon_Kp.spectrum(omega_laser, np.array([omega_E_prime]), 'sigma-', 'sigma+')[0]
DCP_Kp = (I_mm - I_mp) / (I_mm + I_mp)

print(f"\n円偏光度（DCP）:")
print(f"K谷（右カイラル）: DCP = {DCP_K:.3f}")
print(f"K'谷（左カイラル）: DCP = {DCP_Kp:.3f}")
print(f"\n理想的なカイラルフォノンではDCP ≈ ±1を示します")

3.7.2 円二色性スペクトルのシミュレーション

import numpy as np
import matplotlib.pyplot as plt

# フォノン円二色性（CD）スペクトルのシミュレーション

def phonon_cd_spectrum(omega, omega_phonon, gamma, L_phonon, oscillator_strength):
    """
    フォノンモードの円二色性スペクトルを計算

    パラメータ:
    -----------
    omega : array
        周波数範囲 (cm⁻¹)
    omega_phonon : float
        フォノン共鳴周波数 (cm⁻¹)
    gamma : float
        減衰定数 (cm⁻¹)
    L_phonon : float
        フォノン角運動量 (ℏ単位)
    oscillator_strength : float
        振動子強度

    戻り値:
    -------
    CD : array
        円二色性シグナル ΔA = A_σ⁺ - A_σ⁻
    """
    # ローレンツ型吸収（虚部）
    absorption = gamma / ((omega - omega_phonon)**2 + gamma**2)

    # CD ∝ L_phonon × 吸収の微分
    # 簡略化：CD ∝ L_phonon × Im[ε_xy]
    CD = L_phonon * oscillator_strength * absorption

    return CD


# α-石英のパラメータ（Zhu et al. 2018に基づく）
phonon_modes = [
    {'omega': 400, 'gamma': 5, 'L': 0.1, 'f': 0.3},   # 弱いカイラルモード
    {'omega': 700, 'gamma': 8, 'L': 0.2, 'f': 0.5},   # 中程度
    {'omega': 1100, 'gamma': 10, 'L': 0.3, 'f': 1.0}, # 強いSi-O伸縮
]

omega_range = np.linspace(300, 1400, 2000)

# CDスペクトル計算
CD_total = np.zeros_like(omega_range)

for mode in phonon_modes:
    CD_mode = phonon_cd_spectrum(
        omega_range,
        mode['omega'],
        mode['gamma'],
        mode['L'],
        mode['f']
    )
    CD_total += CD_mode

# プロット
fig, axes = plt.subplots(2, 1, figsize=(12, 10))

# 上パネル：個別モード
axes[0].axhline(y=0, color='k', linestyle='-', linewidth=0.8)
for i, mode in enumerate(phonon_modes):
    CD_mode = phonon_cd_spectrum(
        omega_range,
        mode['omega'],
        mode['gamma'],
        mode['L'],
        mode['f']
    )
    axes[0].plot(omega_range, CD_mode, linewidth=2,
                label=f"ω = {mode['omega']} cm⁻¹, L = {mode['L']}ℏ")
    axes[0].fill_between(omega_range, 0, CD_mode, alpha=0.3)

axes[0].set_xlabel('周波数 (cm⁻¹)', fontsize=12)
axes[0].set_ylabel('円二色性 ΔA (任意単位)', fontsize=12)
axes[0].set_title('個別フォノンモードのCD寄与', fontsize=14, fontweight='bold')
axes[0].legend(fontsize=10)
axes[0].grid(True, alpha=0.3)
axes[0].set_xlim(300, 1400)

# 下パネル：全CDスペクトル
axes[1].axhline(y=0, color='k', linestyle='-', linewidth=0.8)
axes[1].plot(omega_range, CD_total, 'b-', linewidth=2.5, label='全CDスペクトル')
axes[1].fill_between(omega_range, 0, CD_total, alpha=0.4, color='blue')

# 強度の最大値を示す
max_idx = np.argmax(np.abs(CD_total))
axes[1].plot(omega_range[max_idx], CD_total[max_idx], 'ro', markersize=10,
            label=f'最大CD: {omega_range[max_idx]:.0f} cm⁻¹')

axes[1].set_xlabel('周波数 (cm⁻¹)', fontsize=12)
axes[1].set_ylabel('円二色性 ΔA (任意単位)', fontsize=12)
axes[1].set_title('α-石英の全フォノンCDスペクトル（シミュレーション）', fontsize=14, fontweight='bold')
axes[1].legend(fontsize=11)
axes[1].grid(True, alpha=0.3)
axes[1].set_xlim(300, 1400)

plt.tight_layout()
plt.savefig('phonon_circular_dichroism.png', dpi=300, bbox_inches='tight')
plt.show()

# CDの温度依存性シミュレーション
temperatures = [10, 100, 200, 300]  # K
fig2, ax2 = plt.subplots(figsize=(12, 7))

for T in temperatures:
    # 温度によるブロードニング: γ(T) = γ₀(1 + 2n_B(T))
    # ボース分布: n_B = 1/(exp(ℏω/k_BT) - 1)
    k_B = 0.695  # cm⁻¹/K（ボルツマン定数）

    CD_T = np.zeros_like(omega_range)
    for mode in phonon_modes:
        # 温度依存ブロードニング
        n_B = 1 / (np.exp(mode['omega'] / (k_B * T)) - 1) if T > 0 else 0
        gamma_T = mode['gamma'] * (1 + 2 * n_B)

        CD_mode = phonon_cd_spectrum(
            omega_range,
            mode['omega'],
            gamma_T,
            mode['L'],
            mode['f']
        )
        CD_T += CD_mode

    ax2.plot(omega_range, CD_T, linewidth=2, label=f'T = {T} K')

ax2.axhline(y=0, color='k', linestyle='-', linewidth=0.8)
ax2.set_xlabel('周波数 (cm⁻¹)', fontsize=12)
ax2.set_ylabel('円二色性 ΔA (任意単位)', fontsize=12)
ax2.set_title('フォノンCDスペクトルの温度依存性', fontsize=14, fontweight='bold')
ax2.legend(fontsize=11, loc='upper right')
ax2.grid(True, alpha=0.3)
ax2.set_xlim(300, 1400)

plt.tight_layout()
plt.savefig('phonon_cd_temperature.png', dpi=300, bbox_inches='tight')
plt.show()

print(f"\nフォノン円二色性の主要特性:")
print(f"最大CDピーク: {omega_range[max_idx]:.1f} cm⁻¹")
print(f"最大CD強度: {CD_total[max_idx]:.3f} (任意単位)")
print(f"\n温度上昇に伴い:")
print(f"- ピークがブロードニング（フォノン-フォノン散乱）")
print(f"- CD強度が減少（熱的乱れ）")

まとめ

演習問題

ΔR/R(t) = 0.05 × cos(250 × 2πc × t) × exp(-t/2 ps)