[论文解读] Deep Potential and Deep Wannier models for KH2PO4&KD2PO4

Bingjia Yang, Pinchen Xie|arXiv (Cornell University)|Apr 11, 2024

Molecular Sensors and Ion DetectionChemistry被引用 3

一句话总结

本研究采用有效质量的路径积分分子动力学（PIMD）方法模拟KDP和DKDP晶体中的核量子效应，结果表明质子穿隧效应会形成量子偶极缺陷，从而降低KDP中的极化强度和熵。氘代作用抑制了这些效应，使体系呈现更经典的、稳定的铁电相，其相变温度更高，热力学性质更优。

ABSTRACT

DPGEN_param.json: Parameters for the active learning agent DPGEN (for training SCAN model) KDP_DP_Train.json & KDP_DW_Train.json: Training parameters for the SCAN Deep Potential and Deep Wannier models KDP_DP_Script.sh & KDP_DW_Script.sh: Training scripts for the SCAN Deep Potential and Deep Wannier models KDP_DP_Model.pb & KDP_DW_Model.pb: The productive SCAN Deep Potential and Deep Wannier models. Tested with DeePMD-Kit v2.0.3 DPGEN_param_vdwDF.json: Parameters for the active learning agent DPGEN (for training vdw-DF model) KDP_DP_vdwDF_Train.json: Training parameters for the vdw-DF Deep Potential model. KDP_DP_vdwDF_Model.pb: The vdw-DF Deep Potential model.

研究动机与目标

• 理解KDP和DKDP晶体中同位素效应的微观起源，特别是质子化与氘代形式之间居里温度（Tc）的显著差异。
• 研究核量子效应（NQEs），尤其是质子穿隧效应，如何影响氢键动力学、极化强度及氢键结合铁电体的相变行为。
• 开发一种计算精确的模型，采用路径积分分子动力学（PIMD）结合有效质量方法，重现KDP和DKDP中的实验Tc值与非中心位移量。
• 阐明量子离域在形成偶极缺陷中的作用，这些缺陷会降低KDP中自发极化强度与剩余熵，而DKDP则表现出更经典的特性。
• 建立亚稳态相中H/D非中心位移与实验测得Tc之间的关联，验证量子模拟在预测铁电性质方面的有效性。

提出的方法

• 采用全原子路径积分分子动力学（PIMD）模拟，将核量子效应（NQEs）纳入KDP和DKDP体系，将质子和氘核视为具有有限零点能的量子粒子。
• 采用有效质量方法（m* = μm）通过调节弹簧常数μ来调控H/D离子的量子离域程度，实现对实验Tc和非中心位移Δxy的标定。
• 利用实验数据对PIMD模拟进行标定：KDP的Tc为123 K，DKDP的Tc为229 K，以及通过中子衍射测得的亚稳态相中Δxy值。
• 计算H/D离子在xy平面内的位置概率密度函数n(δxy)，以分析量子离域性及Tc以上区域氢键的双模特征。
• 分析Tc与Δxy(T→Tc+)之间的相关性，以验证模型的预测能力，并识别控制相变的关键几何参数。
• 将经典极限（无限质量）PIMD与结果对比，确认DKDP比KDP表现出更明显的经典行为，其Tc ≈ 250 K，与实验值一致。

Figure 1: (a) Schematic representation [ 21 ] of the Fdd2 equilibrium structure of KDP. The z-axis is parallel to the $c$ crystallographic axis. (b) Sketch of a ferroelectric KDP unit. K (shared), H (shared), O, and P ions are represented by purple, light blue, red, and yellow spheres. The grey arro

实验结果

研究问题

• RQ1 核量子效应，尤其是质子穿隧效应，如何影响KDP相比DKDP的稳定性与热力学性质？
• RQ2 KDP与DKDP之间居里温度（Tc）存在超过100 K的显著同位素效应，其微观起源是什么？
• RQ3 KDP中质子的量子离域在多大程度上导致偶极缺陷，从而降低自发极化强度与剩余熵？
• RQ4 采用有效质量的路径积分分子动力学能否定量重现KDP和DKDP中的实验Tc值与H/D非中心位移？
• RQ5 质子概率分布n(δxy)的双模特征如何与有序-无序机制关联，并解释相变过程中无长程结构畸变的现象？

主要发现

• KDP中质子穿隧效应产生量子偶极缺陷——即不贡献净极化的局域化、非中心质子构型，导致自发极化强度和剩余熵降低。
• 在DKDP中，氘代作用抑制了质子穿隧效应，使n(δxy)分布更局域化，Tc提高至229 K，高于KDP的123 K，与经典行为一致。
• 采用有效质量m* = 2.5 Da（H）和m* = 32 Da（D）的PIMD模拟成功重现了KDP和DKDP的实验Tc值与Δxy(T→Tc+)值。
• 模拟中观察到Tc与Δxy(T→Tc+)之间存在线性相关性，与实验数据高度吻合，证实H/D非中心位移是控制相变的关键几何参数。
• PIMD的经典极限（无限质量）得到Tc ≈ 250 K，表明DKDP的相变行为主要为经典性质，而KDP则表现出强烈的量子力学特征。
• KDP中的量子离域导致在低有效质量下n(δxy)为单峰分布，但在较高m*时出现双峰分布，与中子衍射数据一致，证实了无序但极化的氢键构型的存在。

Figure 2: (a) The longitudinal spontaneous polarization $|\mathcal{P}_{z}|$ as a function of $T$ . The experimental data are for pure KDP and KD 1.96 H 0.04 PO 4 [ 1 ] . (b) Predicted and experimental [ 18 ] lattice constants $a$ and $b$ as functions of $T$ . (c) The joint PDFs $n(d_{x},0,d_{z})$ (u