[论文解读] Tilt-driven antiferroelectricity in PbZrO$_3$

Konstantin Shapovalov, Massimiliano Stengel|arXiv (Cornell University)|Dec 22, 2021

Ferroelectric and Piezoelectric Materials被引用 4

一句话总结

本文識別出八面體傾角、其梯度與極化之間的三線性耦合——稱為「旋轉極化耦合」——作為 PbZrO₃ 中反鐵電性的關鍵機制。透過基於第一原理的朗道連續體模型，作者表明該機制可解釋共格 $Pbam$ 反鐵電基態的穩定性，統一了先前基於彈性電致伸縮效應與三線性 $\Sigma$--$R$--$S$ 耦合的相互矛盾理論。

ABSTRACT

Antiferroelectricity is a state of matter that has so far eluded a clear-cut definition. Even in the best-known material realization, PbZrO$_3$, the physical nature of the driving force towards an antipolar order has not been settled yet. Here, by building a Landau-like continuum Hamiltonian from first-principles via an exact long-wave approach, we reconcile the existing theories in terms of a single physical mechanism. In particular, we find that a formerly overlooked trilinear coupling between tilts, tilt gradients and polarization provides a surprisingly accurate description of the energetics and structure of the antiferroelectric ground state of PbZrO$_3$. We discuss the relevance of our findings to other ferrielectric and incommensurate polar structures that were recently observed in perovskites.

研究动机与目标

• 解決 PbZrO₃ 中反鐵電性物理起源長期存在的模糊性。
• 調和兩種競爭性理論模型：基於彈性電致伸縮效應的非共格轉變與三線性 $\Sigma$--$R$--$S$ 耦合模型。
• 建立一個統一機制，解釋反極性有序與反螢狀畸變八面體傾角同時凝聚的現象。
• 透過長波極限下對第一原理計算的精確映射，驗證連續體模型。
• 將該機制推廣至其他具有調制極化與傾角結構的鈣钛礦材料。

提出的方法

• 使用來自第一原理密度泛函理論（DFT）計算的精確長波映射，構建類朗道的連續體哈密頓量。
• 將序參數識別為高對稱點處力常數矩陣的特徵向量：極化 ($\mathbf{P}$)、位移 ($\mathbf{u}$) 與八面體傾角 ($\bm{\phi}$)。
• 引入三線性耦合項 $E_{\rm tr} = K\phi_{y}^{2}(\partial P_{y}/\partial x)^{2}$，代表傾角梯度與極化之間的旋轉極化耦合。
• 推導出非共格相形成的廣義不穩定性準則 $|\gamma(\phi_{y0}W^{*})^{4}| > D_{11}G^{*}$，其中包含傾角振幅的影響。
• 執行約束 DFT 漲落計算，以抑制 $\phi_{x}$ 傾角模式，並測試旋轉極化耦合的必要性。
• 透過比較不同傾角振幅下預測的能階與結構與第一原理結果，驗證模型。

实验结果

研究问题

• RQ1 在彈性電致伸縮效應與三線性耦合的相互矛盾解釋下，PbZrO₃ 中反鐵電基態的真正物理機制為何？
• RQ2 是否存在單一機制，能統一解釋非共格漲落、反螢狀畸變傾角與反極性有序的角色？
• RQ3 此前被視為複雜集體模式的 $S$-模式，是否可透過高階耦合物理上解釋為調制傾角模式？
• RQ4 均勻八面體傾角在多大程度上影響非共格相的穩定性？此影響是否改變所提出耦合的相關性？
• RQ5 旋轉極化耦合機制能否推廣至其他鈣鈦礦材料（如 $Pnma$--80 相）中調制相的解釋？

主要发现

• 識別出三線性旋轉極化耦合 $E_{\rm tr} = K\phi_{y}^{2}(\partial P_{y}/\partial x)^{2}$ 為 PbZrO₃ 中反鐵電轉變的主要驅動力。
• 該模型定量重現了 $Pbam$ 基態結構與能量，預測的最優週期與第一原理結果相符。
• 不穩定性準則 $|\gamma(\phi_{y0}W^{*})^{4}| > D_{11}G^{*}$ 即使在包含有限傾角修正的情況下，仍能正確預測非共格性的起始。
• 約束 DFT 計算顯示，抑制 $\phi_{x}$ 傾角（進而抑制旋轉極化耦合）使能量相較於 $Ima2$ 相升高 256 meV/f.u.，證實其必要性。
• $G^{*}$ 與 $W^{*}$ 系數隨傾角振幅增加而減小，但相圖仍保持不穩定，確認該機制的穩健性。
• $S$-模式被重新詮釋為調制傾角模式，為其在 $Pbam$ 結構中的角色提供了明確的物理基礎，解決了先前的詮釋模糊性。