Piezoaxial coupling for strain-selected ferroaxial domain control

技術要約：歪み選択的なフェロアキシャル・ドメイン制御のためのピエゾアキシャル結合

問題提起
フェロアキシャル秩序は、秩序パラメータが軸性ベクトル（axial vector）として変換され、空間反転（$P$）および時間反転（$T$）の両方に対して偶であるという特性を持つフェロイック状態である。この対称性は、フェロ電気的（$P$-奇）および強磁性的（$T$-奇）秩序とは異なる。その結果、標準的な電場や磁場は、フェロアキシャル・ドメインの符号を制御するための普遍的な線形共役場として機能することができない。なぜなら、それらは必要な $P$-偶、$T$-偶の軸性ベクトルの対称性を共有していないからである。非平衡アキシャル場を介してスイッチングを誘起する、円偏光テラヘルツパルスを用いた動的な経路は確立されているが、平衡状態においてフェロアキシャル・ドメインにバイアスをかけることができる、静的で均一な共役場はこれまで未解明であった。歪み－フェロアキシャル結合の実験的兆候（例：1T-TiSe$_2$におけるもの）は現れているものの、すべての点群に対する特定の歪み由来の場を特定する一般的な対称性の枠組みや、結果として生じるエネルギー分裂の第一原理計算による検証は確立されていない。

手法
著者らは、選択されたフェロアキシャル秩序パラメータ（$A$）と同じ既約表現で変換される、静的な均一歪み多項式の対称性に基づく階層構造を開発した。

1. 対称性の定式化： 著者らは、自由エネルギーにおけるピエゾアキシャル結合項 $F_A = -h(\boldsymbol{\varepsilon})A$ を定義した。ここで、$h(\boldsymbol{\varepsilon})$ は歪みに由来するアキシャル場である。著者らは、非パイロアキシャル（non-pyroaxial）親点群に対して、最低次の許容される歪み多項式を系統的に導出し、結晶族（斜方晶、正方晶、三方晶、六方晶、および立方晶）ごとに分類した。
2. 第一原理計算による検証： 理論的枠組みを検証するため、著者らは三方晶フェロアキシャル化合物である Na$_2$BaMg(PO$_4$)$_2$ に対して密度汎関数理論（DFT）計算を行った。この物質では、秩序パラメータは三方軸周りの PO$_4$ 四面体の集団的な回転に対応している。
3. 計算プロトコル： 以下の3つの相補的な計算を用いている。
   • クランプ座標ドメイン分裂： 固定歪力下で、あらかじめ構築された反対方向の回転符号を持つフェロアキシャル状態間のエネルギー差を計算する。
   • ダブルウェル・ポテンシャル・スキャン： 回転座標に沿ったエネルギープロファイルをマッピングし、ポテンシャルの傾きを観察する。
   • 固定歪緩和： パラアキシャル（$\phi=0$）構造から固定歪の下で原子位置を緩和させ、歪み場のみによって特定のドメイン・ベイシンが選択されるかどうかを決定する。

主要な貢献と結果

• 歪み場の対称性階層： 本論文は、主要な次数のピエゾアキシャル場が親点群に依存することを確立した：
  • 斜方晶： 歪みに対して線形（せん断場）。
  • 正方晶： 歪みに対して二次。
  • 三方晶および六方晶： 歪みに対して三次。
  • 立方晶： 2つのクラスに分かれる。「Cubic-I」群（23, $\bar{m}3$）は、選択された軸に対して線形フル歪み場を許容する一方、「Cubic-II」群（432, $\bar{4}3m$, $m\bar{3}m$）は線形場を禁止し、選択された軸（[001], [110], または [111]）に応じて二次または三次結合を必要とする。
• 三方晶の特異性： 三方晶系（例：Na$_2$BaMg(PO$_4$)$_2$）において、主要な底面内場は三次である：$h(\boldsymbol{\varepsilon}) \propto Y(3X^2 - Y^2)$、ここで $X$ と $Y$ は偏差歪成分である。符号付き歪振幅 $\varepsilon_u$ と主軸角 $\theta$ を用いると、これは $h \propto \varepsilon_u^3 \sin(6\theta)$ となる。
• 第一原理による検証：
  • 角度および振幅依存性： 計算されたドメイン分裂（$\Delta E_{dom}$）は、予測された $\sin(6\theta)$ の角度依存性と、三次スケーリング $\Delta E_{dom} \propto \varepsilon_u^3$ を厳密に追従する。
  • ポテンシャルの傾き： エネルギースキャンにより、歪み場が $\phi$ に関して奇関数である項を導入し、ダブルウェル・ポテンシャルにバイアスを与え、2つのドメイン間の縮退を解くことが確認された。
  • 歪み選択的緩和： パラアキシャル構造から開始した固定歪緩和は、印加された歪みの符号（$\varepsilon_u$）または歪み軸の向き（$\theta$）のみに基づいて、反対のフェロアキシャル最小値へと進化する。
• 実験的実現可能性： 著者らは、$\sim 100$ 単位格子からなる相関核について、歪み $|\varepsilon_u| = 0.020$ における歪み誘起バイアスが転移温度（$T_c \approx 540$ K）付近で熱力学的に有意になることを推定しており、これは歪み場冷却実験の実行可能性を示唆している。

意義
本論文は、静的な均一歪みが、フェロアキシャル秩序に対する対称性的に許容された共役場であることを確立した。また、様々な結晶系においてフェロアキシャル・ドメインを制御するために必要な特定の歪み多項式を特定するための、一般的な対称性の枠組みを提供している。歪みが平衡状態においてフェロアキシャル・ドメインをエネルギー的に選択できることを示すことで、本研究は、磁場冷却に類似した「歪み場冷却」プロトコルを提案している。これは、非平衡な動的場に依存することなく、フェロアキシャル・ドメインのエンジニアリングおよび制御への道を開くものである。この枠組みは、単純な構造転移を超えて、電子的および有限波長ベクトルの軸性状態（例：CDW状態）にも拡張可能であるが、それは歪み場が誘導された一様なアキシャル・コンポジットに結合すると理解される場合に限られる。著者らは、本フレームワークは一般的であるが、最も直接的な実験的実現は、ここで研究されたグラセライト型リン酸塩のような、適切な構造的フェロアキシャル転移において期待されると述べている。