Terahertz oscillation of $180^{\circ}$ domain walls in ferroelectric membranes

技術要約：強誘電体メンブレンにおける180ºドメイン壁のテラヘルツ振動

問題提起
強誘電体におけるドメイン壁（DW）の集団励起は、ナノエレクトロニクスやオプトエレクトロニクスデバイスへの応用が期待されるものの、根本的な理解が不十分な現象である。動的フェーズフィールドモデリング（DPFM）は、バルク結晶や薄膜における分極応答を予測するために成功裏に適用されてきたが、マルチドメイン強誘電体メンブレンにおける分極波と音響波の結合ダイナミクスについては、比較的未開拓の領域である。特に、一軸歪みを加えたイジング型の180ºドメイン壁の高周波ダイナミクスに関して、とりわけ非従来型のスライディングモードや、キャビティ強化された音響波との相互作用に関する理解が不足している。

手法
著者らは、一軸歪みを加えたバリウムチタン酸塩（BaTiO$_3$）メンブレンに含まれるイジング型180ºドメイン壁の、結合された分極・歪みダイナミクスを調査するために、動的フェーズフィールドシミュレーション（DPFM）を用いている。

• シミュレーション・フレームワーク： 本モデルは、ランダウ・ギンツブルグ・デボーン理論によって支配される分極の時間発展方程式を、機械的変位に関する弾性力学方程式と結合させて解くものである。このフレームワークには、脱分極電界や放射電界を含む、分極、歪み、および電磁波（EM波）の間の相互作用が含まれている。
• システム構成： シミュレーションでは、歪み軸（$x_1$）方向に周期境界条件を、横方向の軸（$x_2, x_3$）方向に自由表面を持つ、BaTiO$_3$メンブレンを表す3次元離散化システムを利用する。システムには、300 Kにおいて0.5%から1.1%の範囲の一軸歪みが加えられる。
• 励起および解析： 研究には主に2つのアプローチが含まれる。(1) ドメインパターンの自発的形成を観察し、時空間分極データの2次元フーリエ変換から固有の分散関係を抽出すること、および (2) 平衡状態にあるストライプドメインに対して広帯域の単一サイクルTHz電界パルスを印加し、特定の固有モードを選択的に励起・再構成することである。
• 解析的支援： シミュレーション結果は、横方向および縦方向の分極波の分散関係、結合された音響・光学フォノン分散、およびバンド折り畳みされた横方向音響フォノンに関する解析的な計算を用いて解釈される。さらに、簡略化された1次元数値固有モード解析をメンブレンのケースに拡張し、標準的なDWモードを特定している。

主要な貢献および結果
本研究は、テラヘルツ（THz）領域におけるいくつかの明確な集団モードを特定し、特徴付けを行っている：

1. 固有分極波： 単一ドメインおよびマルチドメイン領域において、シミュレーションは固有の光学フォノンモード、すなわち、4–6 THzの範囲にある横方向光学（A1(TO)）モード（ソフトモード）および1–4 THzの範囲にある縦方向光学（E(LO)）モードの存在を裏付けている。これらのモードの分散関係は、自由エネルギー景観の曲率に基づく解析的な予測と一致している。

2. ドメイン壁（DW）固有モード： 著者らは3つの異なるDWモードを特定した：

   • DWブリージングモード： ドメイン壁の中心は固定されたまま、ドメイン壁の厚さが振動するモード。このモードは2.94–3.01 THz付近の周波数で観察される。
   • 高次DWモード： 3次固有モードに類似した非対称モードであり、3.35–3.49 THz付近で観察される。
   • 非従来型のDWスライディングモード： ドメイン壁がその法線方向に沿って滑る、新規なサブTHzモード（0.17–0.90 THz）。以前に報告されているゼロ周波数のゴールドストーン・スライディングモードや表面電荷駆動のスライディングモードとは異なり、このモードは非ゼロの共鳴周波数を示す。

3. 非従来型スライディングモードのメカニズム： 本研究は、この非従来型スライディングモードの復元力が、表面電荷ではなくバルク分極束縛電荷（$\rho_b = -\nabla \cdot \Delta \mathbf{P}$）によって駆動されていることを明らかにしている。空間的に変化する縦方向フォノン成分（$\Delta P_1$）が、歪み軸に沿ったバルク脱分極電界（$E_1^d$）を生成し、これが横方向のスライディング運動に対する復元力として作用する。このメカニズムは、表面脱分極電界に依存する従来のモデルとは明確に異なる。

4. 歪みによる変調： 固有およびDWモードの両方の周波数は、一軸歪みによって調整可能であることが示されている。

   • A1(TO)およびブリージングモードの周波数は、自由エネルギー景観の曲率が増大するため、一般に歪みの増加とともに上昇する。
   • 非従来型スライディングモードの周波数は、低歪み（0.5%–0.8%）では比較的安定しているが、高歪み（0.9%–1.1%）では大幅に上昇する。この上昇は、ドメイン壁の幅が狭まることでバルク束縛電荷が集中し、脱分極電界が強化されることに起因する。
   • ドメイン壁の有効質量密度は、歪みの増加に伴い減少すると推定されており、これは壁の厚さの減少と一致している。

5. キャビティ強化共鳴： メンブレンは分極波と音響波の両方のキャビティとして機能する。本研究は、圧電結合により、E(LO)モードがバンド折り畳みされた横方向音響（TA）フォノンによって共鳴励起され得ることを示しており、これはマグノン・フォノン結合に類似した現象である。しかし、現在のシステム構成においては、A1(TO)モードおよびDWモードはこの共鳴結合を示さない。

意義および主張
本論文は、強誘電体ドメイン壁の高周波ダイナミクスに関する新たな知見を提供し、特に、バルク電荷によって駆動される、サブTHzの共鳴周波数を持つ非従来型のスライディングモードの存在を強調することを主張している。著者らは、これらの知見が以下の機会を示唆していると述べている：

• 歪み制御： THz領域におけるフォノン-DW共鳴の、歪みによる調整可能な制御の実現。
• デバイス応用： ドメイン設計された強誘電体ナノメンブレンに基づく、再構成可能なTHzおよび光学デバイスの実現の可能性。
• 手法の拡張： 本理論的枠組みおよび数値的手法は、他のドメイン構成（例：90ºまたは109ºの壁、混合ブロホ・ネール・アイザック特性）や、トポロジカルに非自明な極性テクスチャにおける集団モードの発見へと拡張可能である。
• ハイブリッドシステム： コヒーレントな分極波と音響波の間のキャビティ強化された強結合を介した、ハイブリッド状態の創出の可能性であり、これは量子通信やセンシングへの応用にも利益をもたらす可能性がある。

著者らは、シミュレーションでは特定の周期的なドメインパターンを用いているものの、予測されたDWモードの高度な局在性は、これらが完全な周期性を必要とせず、180ºストライプドメインを含む様々な強誘電体メンブレンにおいて励起可能であることを示唆していると強調している。