搜索

x

留言板

姓名
邮箱
手机号码
标题
留言内容
验证码

引用本文:
Citation:

郭常青, 杨乐陶, 王静, 黄厚兵
cstr: 32037.14.aps.74.20250334

Phase-field simulation of domain switching in ferroelectric trilayer films under bending-induced strain gradient

GUO Changqing, YANG Letao, WANG Jing, HUANG Houbing
cstr: 32037.14.aps.74.20250334
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 柔性铁电材料在可穿戴电子领域具有重要应用前景, 然而其动态弯曲过程中应变梯度与极化翻转的力-电耦合物理机制仍缺乏系统性的研究. 本研究基于相场模拟, 系统地探讨了(SrTiO3)10/(PbTiO3)10/(SrTiO3)10三层异质膜在U型和N型弯曲下的畴结构演化及其宏观电学响应. 研究表明, 通过改变弯曲变形方向可以产生方向相反的挠曲电场, 导致电滞回线发生相应方向的偏移. 此外, 弯曲应变和应变梯度可驱动极性涡旋态与单畴态之间的拓扑相变, 其中界面静电能、弹性约束及梯度能的协同作用对拓扑结构的稳定性起关键作用. 本研究揭示了弯曲变形通过力-电耦合效应实现畴构型与电响应的定向调控机制, 为高密度柔性存储器和能量收集器件的跨尺度设计奠定了理论基础, 并进一步拓展了拓扑态工程在柔性电子领域的应用前景.
    Flexible ferroelectric materials possess considerable potentials for wearable electronics and bio-inspired devices, yet their mechano-electric coupling mechanisms under dynamic bending conditions remain incompletely understood. In his work, the effects of bending deformation on domain structures and macroscopic ferroelectric responses in (SrTiO3)10/(PbTiO3)10/(SrTiO3)10 flexible ferroelectric trilayer films are systematically investigated using phase-field simulations. By constructing computational models for upward-concave (U-shaped) and downward-concave (N-shaped) bending configurations, the strain distribution and its regulation mechanism on polarization patterns under different curvature radii are analyzed. The results reveal distinct strain gradients across bending modes: U-shaped bending induces compressive strain in the upper layer and tensile strain in the lower layer, generating a negative out-of-plane strain gradient. Conversely, N-shaped bending reverses this strain distribution. Such inhomogeneous strains drive significant polarization reconfiguration within the PTO layer. At a moderate curvature (large R), the system retains stable vortex-antivortex pairs. Reducing bending radius (smaller R) promotes divergent topological transitions—U-shaped bending facilitates vortex pair transformation into zigzag-like domains, while N-shaped bending drives vortex-to-out-of-plane c-domain evolution. Notably, bending-induced strain gradients impose transverse flexoelectric fields that markedly change trilayer hysteresis loops. U-shaped bending introduces a negative flexoelectric field, shifting loops rightward with maximum polarization (Pmax) decreasing. In contrast, N-shaped bending generates a positive field, enhancing Pmax via leftward loop shifting. The polarization switching analysis under electric field further demonstrates bending-mediated control of domain evolution pathway and reversal dynamics. These findings not only elucidate profound bending effects on flexible ferroelectrics’ domain architectures and functional properties but also provide theoretical guidance for designing strain-programmable ferroelectric memories, adaptive sensors, and neuromorphic electronics.
      通信作者: 黄厚兵, hbhuang@bit.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 92463306, 52372100, 52472119)和北京市自然科学基金(批准号: 2242057)资助的课题.
      Corresponding author: HUANG Houbing, hbhuang@bit.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 92463306, 52372100, 52472119) and the Natural Science Foundation of Beijing, China (Grant No. 2242057).
    [1]

    [2]

    [3]

    [4]

    [5]

    [6]

    [7]

    [8]

    [9]

    [10]

    [11]

    [12]

    [13]

    [14]

    [15]

    [16]

    [17]

    [18]

    [19]

    [20]

    [21]

    [22]

    [23]

    [24]

    [25]

    [26]

    [27]

    [28]

    [29]

    [30]

    [31]

    [32]

    [33]

    [34]

    [35]

    [36]

    [37]

    [38]

  • 变量 数值 变量 数值
    PTO $ {\alpha _1}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} m {\cdot} {{\mathrm{C}}^{ - 2}}) $ $ - 1.706 $ $ {{{Q}}_{11}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0.089
    $ {\alpha _{11}}/({10^7}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^5} {\cdot} {{\mathrm{C}}^{ - 4}}) $ $ - 7.3 $ $ {{{Q}}_{12}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ –0.026
    $ {\alpha _{12}}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^5} {\cdot} {{\mathrm{C}}^{ - 4}}) $ $ 7.5 $ $ {{{Q}}_{44}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0.0675
    $ {\alpha _{111}}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^9} {\cdot} {{\mathrm{C}}^{ - 6}}) $ $ 2.6 $ $ {{{G}}_{11}}/({10^{ - 10}}{\text{ }}{\mathrm{N}} {\cdot} {{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ $ 1.44 $
    $ {\alpha _{112}}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^9} {\cdot} {{\mathrm{C}}^{ - 6}}) $ $ 6.1 $ $ {{{G}}_{12}}/({\mathrm{N}} {\cdot} {{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0
    $ {\alpha _{123}}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^9} {\cdot} {{\mathrm{C}}^{ - 6}}) $ $ - 3.7 $ $ G_{44},G_{44}'/(10^{-11}\text{ }\mathrm{N}{\cdot}\mathrm{m}^4{\cdot}\mathrm{C}^{-2}) $ $ 7.2 $
    $ {{c}_{11}}/({10^{11}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 2.3 $ $ {{f}_{11}}/{\mathrm{V}} $ 1.6
    $ {{c}_{12}}/({10^{11}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 1 $ $ {{f}_{12}}/{\mathrm{V}} $ –0.8
    $ {{c}_{44}}/({10^{10}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 7 $ $ {{f}_{44}}/{\mathrm{V}} $ 0.15
    STO $ {\alpha _1}/({10^8}{\text{ }}{\mathrm{J}} {\cdot} {\mathrm{m}} {\cdot} {{\mathrm{C}}^{ - 2}}) $ $ 2.017 $ $ {{{Q}}_{44}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0.00957
    $ {\alpha _{11}}/({10^9}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^5} {\cdot} {{\mathrm{C}}^{ - 4}}) $ $ 1.7 $ $ {{{G}}_{11}}/({10^{ - 10}}{\text{ }}{\mathrm{N}} {\cdot} {{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ $ 1.44 $
    $ {\alpha _{12}}/({10^9}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^5} {\cdot} {{\mathrm{C}}^{ - 4}}) $ $ 4.45 $ $ {{{G}}_{12}}/({\mathrm{N}} {\cdot} {{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0
    $ {{c}_{11}}/({10^{11}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 3.3 $ $ {{{G}}_{44}}, {{G}}_{44}'/({10^{ - 11}}{\text{ }}{\mathrm{N}} {\cdot} {{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ $ 7.2 $
    $ {{c}_{12}}/({10^{11}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 1 $ $ {{f}_{11}}/{\mathrm{V}} $ –3.21
    $ {{c}_{44}}/({10^{11}}{\text{ }}{\mathrm{J}} {\cdot} {{\mathrm{m}}^{ - 3}}) $ $ 1.25 $ $ {{f}_{12}}/{\mathrm{V}} $ 1.47
    $ {{{Q}}_{11}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ 0.0457 $ {{f}_{44}}/{\mathrm{V}} $ 1.07
    $ {{{Q}}_{12}}/({{\mathrm{m}}^4} {\cdot} {{\mathrm{C}}^{ - 2}}) $ –0.0135 εr (PTO/STO) 20
    下载: 导出CSV

    U型弯曲-
    R/nm
    εxx,z/
    (106 m–1)
    Ec/
    (kV·cm–1)
    Pmax/
    (μC·cm–2)
    Pr/
    (μC·cm–2)
    未弯曲 0 0 43.90 0
    1200 –0.61 2.78 43.72 –0.29
    937 –0.77 2.78 43.64 –0.31
    600 –1.14 5.56 43.52 –0.46
    400 –1.82 11.11 43.30 –0.57
    300 –2.42 13.89 43.05 –0.64
    240 –3.03 16.67 42.86 –0.69
    下载: 导出CSV

    N型弯曲-
    R/nm
    εxx,z/
    (106 m–1)
    Ec/
    (kV·cm–1)
    Pmax/
    (μC·cm–2)
    Pr/
    (μC·cm–2)
    未弯曲 0 0 43.90 0
    1200 0.61 –2.78 44.13 0.21
    600 1.14 –5.56 44.33 0.43
    400 1.82 –8.33 44.53 0.67
    300 2.42 –13.89 44.78 0.93
    240 3.03 –16.67 44.88 1.18
    200 3.61 –19.44 45.15 1.42
    下载: 导出CSV
    Baidu
  • [1]

    [2]

    [3]

    [4]

    [5]

    [6]

    [7]

    [8]

    [9]

    [10]

    [11]

    [12]

    [13]

    [14]

    [15]

    [16]

    [17]

    [18]

    [19]

    [20]

    [21]

    [22]

    [23]

    [24]

    [25]

    [26]

    [27]

    [28]

    [29]

    [30]

    [31]

    [32]

    [33]

    [34]

    [35]

    [36]

    [37]

    [38]

  • [1] 解晓洁, 孙俊松, 秦吉红, 郭怀明. 弯曲应变下六角晶格量子反铁磁体的赝朗道能级. 必威体育下载 , 2024, 73(2): 020202. doi: 10.7498/aps.73.20231231
    [2] 刘续希, 高士森, 喇永孝, 玉栋梁, 柳文波. Zr-2.5Sn合金高温腐蚀过程的相场模拟. 必威体育下载 , 2024, 73(14): 148201. doi: 10.7498/aps.73.20240393
    [3] 廖宇轩, 申文龙, 吴学志, 喇永孝, 柳文波. 陶瓷型复合燃料烧结过程的相场模拟研究. 必威体育下载 , 2024, 73(21): 210201. doi: 10.7498/aps.73.20241112
    [4] 王继光, 李珑玲, 邱嘉图, 陈许敏, 曹东兴. 钙钛矿超晶格材料界面二维电子气的调控. 必威体育下载 , 2023, 72(17): 176801. doi: 10.7498/aps.72.20230573
    [5] 夏文强, 赵彦, 刘振智, 鲁晓刚. 应变诱发四方相小角度对称倾侧晶界位错反应的晶体相场模拟. 必威体育下载 , 2022, 71(9): 096102. doi: 10.7498/aps.71.20212278
    [6] 姜彦博, 柳文波, 孙志鹏, 喇永孝, 恽迪. 外加应力作用下 UO2 中空洞演化过程的相场模拟. 必威体育下载 , 2022, 71(2): 026103. doi: 10.7498/aps.71.20211440
    [7] 杨朝曦, 柳文波, 张璁雨, 贺新福, 孙正阳, 贾丽霞, 师田田, 恽迪. Fe-Cr合金晶界偏析及辐照加速晶界偏析的相场模拟. 必威体育下载 , 2021, 70(11): 116101. doi: 10.7498/aps.70.20201840
    [8] 刘迪, 王静, 王俊升, 黄厚兵. 相场模拟应变调控PbZr(1–x)TixO3薄膜微观畴结构和宏观铁电性能. 必威体育下载 , 2020, 69(12): 127801. doi: 10.7498/aps.69.20200310
    [9] 张军, 陈文雄, 郑成武, 李殿中. γ-α相变中不同晶界特征下铁素体生长形貌的相场模拟. 必威体育下载 , 2017, 66(7): 070701. doi: 10.7498/aps.66.070701
    [10] 宋宗根, 邓科, 何兆剑, 赵鹤平. 高对称型声子晶体自准直弯曲及分束. 必威体育下载 , 2016, 65(9): 094301. doi: 10.7498/aps.65.094301
    [11] 高英俊, 全四龙, 邓芊芊, 罗志荣, 黄创高, 林葵. 剪切应变下刃型位错的滑移机理的晶体相场模拟. 必威体育下载 , 2015, 64(10): 106104. doi: 10.7498/aps.64.106104
    [12] 段培培, 邢辉, 陈志, 郝冠华, 王碧涵, 金克新. 镁基合金自由枝晶生长的相场模拟研究. 必威体育下载 , 2015, 64(6): 060201. doi: 10.7498/aps.64.060201
    [13] 曹晔, 裴庸惟, 童峥嵘. 仅用一根局部微结构长周期光纤光栅实现温度与弯曲曲率的同时测量. 必威体育下载 , 2014, 63(2): 024206. doi: 10.7498/aps.63.024206
    [14] 冯小勤, 贾建明, 陈贵宾. 弯曲BN纳米片的电子性质及其调制. 必威体育下载 , 2014, 63(3): 037101. doi: 10.7498/aps.63.037101
    [15] 杜立飞, 张蓉, 邢辉, 张利民, 张洋, 刘林. 横向限制下凝固微观组织演化的相场法模拟. 必威体育下载 , 2013, 62(10): 106401. doi: 10.7498/aps.62.106401
    [16] 潘诗琰, 朱鸣芳. 双边扩散枝晶生长的定量相场模型. 必威体育下载 , 2012, 61(22): 228102. doi: 10.7498/aps.61.228102
    [17] 魏承炀, 李赛毅. 温度梯度对晶粒生长行为影响的相场模拟. 必威体育下载 , 2011, 60(10): 100701. doi: 10.7498/aps.60.100701
    [18] 邵建立, 何安民, 段素青, 王裴, 秦承森. 单轴应变驱动铁bcc—hcp相转变的微观模拟. 必威体育下载 , 2010, 59(7): 4888-4894. doi: 10.7498/aps.59.4888
    [19] 王刚, 徐东生, 杨锐. Ti-6Al-4V合金中片层组织形成的相场模拟. 必威体育下载 , 2009, 58(13): 343-S348. doi: 10.7498/aps.58.343
    [20] 吴自玉, 汪克林. 弯曲时空的半整数自旋场方程. 必威体育下载 , 1985, 34(5): 588-593. doi: 10.7498/aps.34.588
计量
  • 文章访问数:  459
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2025-03-13
  • 修回日期:  2025-04-14
  • 上网日期:  2025-04-24
  • 刊出日期:  2025-06-20

返回文章
返回
Baidu
map