-
Ferroelectric material possesses spontaneous polarization at room temperature, which can be switched by an external electric field. The diverse domain structures within ferroelectric materials, consisting of polarizations in various directions, often significantly affect their physical properties and practical applications. Numerical simulations can aid in comprehending and validating the complex domains observed in experiments. They can also provide guidance for controlling such structures. One popular method for finding dipole configurations is to create an energy model and employ it in Monte-Carlo simulations to find dipole configuration. However, since these simulations usually reaches the ground state of the system (the state with the lowest global energy), they often miss the dipole configurations of interest, such as topological domain structures, which are usually metastable. Here, in order to simulate complex domain, we introduce Brown's equation, which is originally used for micromagnetic simulation, into the large-scale simulation of ferroelectric materials. Using the effective Hamiltonian as the energy model, we derive the Brown's equations with respect to the electric dipoles in ferroelectric materials, and invesitgate perovskites such as ${\rm{BaTiO_3}}$ bulk,${\rm{PbTiO_3}}$ bulk, and${\rm{SrTiO_3}}$ /${\rm{PbTiO_3}}$ /${\rm{SrTiO_3}}$ sandwiched structures. We demonstrate the reliability and feasibility of Brown's equation in ferroelectrics through the simulation of${\rm{BaTiO_3}}$ bulk and${\rm{PbTiO_3}}$ bulk, which are consistent with experiments. Then, using Brown's equation derived in our work, we obtain various domain structures in${\rm{SrTiO_3}}$ /${\rm{PbTiO_3}}$ /${\rm{SrTiO_3}}$ sandwiched structures, including periodic stripe domains and vortex domains. The simulation results are compared with related exprimental results.-
Keywords:
- topological domains /
- Brown’s equation /
- perovskites
[1] Seidel J, Ramesh R, Scott J, Catalan G 2012 Rev. Mod. Phys. 84 119Google Scholar
[2] Nataf G F, Guennou M, Gregg J M, Meier D, Hlinka J, Salje E K H, Kreisel J 2020 Nat. Rev. Phys. 2 634Google Scholar
[3] Scott J 2016 Ferroelectrics 503 117Google Scholar
[4] 谭丛兵, 钟向丽, 王金斌 2020 69 127702Google Scholar
Tan C B, Zhong X L, Wang J B 2020 Acta Phys. Sin. 69 127702Google Scholar
[5] Lu X M, Huang F Z, Zhu J S 2020 Acta Phys. Sin. 69Google Scholar
[6] Tian G, Yang W, Chen D, Fan Z, Hou Z, Alexe M, Gao X 2019 Natl. Sci. Rev. 6 684Google Scholar
[7] Das S, Tang Y L, Hong Z, Gonçalves M A P, McCarter M R, Klewe C, Nguyen K X, Gómez-Ortiz F, Shafer P, Arenholz E, Stoica V A, Hsu S L, Wang B, Ophus C, Liu J F, Nelson C T, Saremi S, Prasad B, Mei A B, Schlom D G, Íñiguez J, García-Fernández P, Muller D A, Chen L Q, Junquera J, Martin L W, Ramesh R 2019 Nature 568 368Google Scholar
[8] Huang J, Tan P, Wang F, Li B 2022 Crystals 12 786Google Scholar
[9] Das S, Ghosh A, McCarter M R, Hsu S L, Tang Y L, Damodaran A R, Ramesh R, Martin L W 2018 APL Mater. 6 100901Google Scholar
[10] Seidel J 2016 Topological Structures in Ferroic Materials, volume 228 of Springer Series in Materials Science (Cham: Springer International Publishing)
[11] Lu L, Nahas Y, Liu M, Du H, Jiang Z, Ren S, Wang D, Jin L, Prokhorenko S, Jia C L, Bellaiche L 2018 Phys. Rev. Lett. 120 177601Google Scholar
[12] Marton P, Rychetsky I, Hlinka J 2010 Phys. Rev. B 81 144125Google Scholar
[13] Zhong W, Vanderbilt D, Rabe K M 1995 Phys. Rev. B 52 6301Google Scholar
[14] Chen L Q 2008 J. Am. Ceram. Soc. 91 1835Google Scholar
[15] Liu J, Chen W, Wang B, Zheng Y 2014 Materials 7 6502Google Scholar
[16] Jiang Z, Xu B, Li F, Wang D, Jia C L 2015 Phys. Rev. B 91 014105Google Scholar
[17] Wojdeł J C, Íñiguez J 2014 Phys. Rev. Lett. 112 247603Google Scholar
[18] Zhong W, Vanderbilt D, Rabe K M 1994 Phys. Rev. Lett. 73 1861Google Scholar
[19] Liu Y, Tang Y L, Zhu Y L, Wang W Y, Ma X L 2016 Adv. Mater. Interfaces 3 1600342Google Scholar
[20] Brown W F 1941 Phys. Rev. 60 139Google Scholar
[21] Brown W F 1940 Phys. Rev. 58 736Google Scholar
[22] Aharoni A, et al. 2000 Introduction to the Theory of Ferromagnetism (Vol. 109) (Clarendon Press)
[23] García-Palacios J L, Lázaro F J 1998 Phys. Rev. B 58 14937Google Scholar
[24] Kumar D, Adeyeye A O 2017 J. Phys. D: Appl. Phys. 50 343001Google Scholar
[25] Zhu B, Lo C C H, Lee S J, Jiles D C 2001 J. Appl. Phys. 89 7009Google Scholar
[26] Shen K, Tatara G, Wu M W 2011 Phys. Rev. B 83 085203Google Scholar
[27] Li B L, Liu X P, Fang F, Zhu J L, Liu J M 2006 Phys. Rev. B 73 014107Google Scholar
[28] Wang D, Liu J, Zhang J, Raza S, Chen X, Jia C L 2019 Comput. Mater. Sci. 162 314Google Scholar
[29] Hong Z, Damodaran A R, Xue F, Hsu S L, Britson J, Yadav A K, Nelson C T, Wang J J, Scott J F, Martin L W, et al. 2017 Nano Lett. 17 2246Google Scholar
[30] Chen L Q, Zhao Y 2022 Prog. Mater. Sci. 124 100868Google Scholar
[31] Petković M, Petković L, Herceg ff 2010 J. Comput. Appl. Math. 233 1755Google Scholar
[32] Sugiura H, Hasegawa T 2019 J. Comput. Appl. Math. 358 136Google Scholar
[33] 黄云清, 舒适, 陈艳萍 2009 数值计算方法 (北京: 科学出版社)
Huang Y Q, Shu S, Chen Y P 2009 Numerical Computation Method (Beijing: Science Press) (in Chinese)
[34] Nishimatsu T, Iwamoto M, Kawazoe Y, Waghmare U V 2010 Phys. Rev. B 82 134106Google Scholar
[35] Nishimatsu T, Aoyagi K, Kiguchi T, J Konno T, Kawazoe Y, Funakubo H, Kumar A, V Waghmare U 2012 J. Phys. Soc. Jpn. 81 124702Google Scholar
[36] Nishimatsu T, Grünebohm A, Waghmare U V, Kubo M 2016 J. Phys. Soc. Jpn. 85 114714Google Scholar
[37] Kwei G H, Lawson A C, Billinge S J L, Cheong S W 1993 J. Phys. Chem. 97 2368Google Scholar
[38] Bersuker I B 1966 Phys. Lett. 20 589Google Scholar
[39] Ravel B, Stern E A, Vedrinskii R I, Kraizman V 1998 Ferroelectrics 206 407Google Scholar
[40] Liu J, Jin L, Jiang Z, Liu L, Himanen L, Wei J, Zhang N, Wang D, Jia C L 2018 J. Chem. Phys. 149 244122Google Scholar
[41] Tinte S, Stachiotti M, Sepliarsky M, Migoni R, Rodriguez C 1999 J. Phys.: Condens. Matter 11 9679Google Scholar
[42] Zhang F, Zhang J, Jing H, Li Z, Wang D, Jia C L 2022 Phys. Rev. B 105 024106Google Scholar
[43] Tadmor E, Waghmare U, Smith G, Kaxiras E 2002 Acta Mater. 50 2989Google Scholar
[44] Sani A, Hanfland M, Levy D 2002 J. Solid State Chem. 167 446Google Scholar
[45] Abid A Y, Sun Y, Hou X, Tan C, Zhong X, Zhu R, Chen H, Qu K, Li Y, Wu M, et al. 2021 Nat. Commun. 12 1Google Scholar
[46] Lemee N, Infante I C, Hubault C, Boulle A, Blanc N, Boudet N, Demange V, Karkut M G 2015 ACS Appl. Mater. Interfaces 7 19906Google Scholar
[47] Pereiro M, Yudin D, Chico J, Etz C, Eriksson O, Bergman A 2014 Nat. Commun. 5 1Google Scholar
[48] Sichuga D, Bellaiche L 2011 Phys. Rev. Lett. 106 196102Google Scholar
[49] Kim J H, Lange F F 1999 J. Mater. Res. 14 1626Google Scholar
[50] Zhang J, Wang Y J, Liu J, Xu J, Wang D, Wang L, Ma X L, Jia C L, Bellaiche L 2020 Phys. Rev. B 101 060103Google Scholar
-
图 1 (a)—(c)是BTO系统模拟得到的三种偶极子极化构型, 宏观极化分别沿
$ \langle111\rangle $ ,$ \langle110\rangle $ ,$ \langle100\rangle $ 方向; (d)计算得到的所有样本中偶极子构型能量和宏观极化方向的关系, 横轴表示不同次计算获得的结果, 纵轴表示每个偶极子的平均能量(meV); (e)说明了图(a)—(c)中箭头颜色表示的偶极子极化方向, 不同颜色的箭头表明偶极子极化方向处于不同的象限. 下节PTO块体的构型也沿用了图(e)表示的颜色方向Figure 1. (a)–(c) Three types of dipole polarization configurations in BTO simulation, where the macroscopic polarizations are along
$ \langle111\rangle $ ,$ \langle 110 \rangle $ ,$ \langle100\rangle $ , respectively; (d) energy (per dipole) with respect to the calculated configurations obtained from simulations; (e) dipole polarization direction indicating the colors used by the arrows in panels (a)–(c). The colors shown in panel (e) are also used for PTO bulk in the next section图 3 (a) STO
$ _{6} $ /PTO$ _{8} $ /STO$ _{6} $ 体系模拟结果, 红色箭头和蓝色箭头分别表示极化沿$ \left[001\right] $ 正方向和负方向; (b) 图(a)中沿$ \left[100\right] $ 方向第28层偶极子平面投影图, 即图(a)黑色纵切面, 黄色圆表示涡旋畴结构Figure 3. (a) Simulation results of STO
$ _{6} $ /PTO$ _{8} $ /STO$ _{6} $ system, where red arrows and blue arrows indicate that the polarization is along the positive and negative directions of$ \left[001\right] $ respectively; (b) projection of dipoles in the 28th$ (100) $ plane along the$ \left[100\right] $ direction (the black arrows in panel (a)). The yellow circles show the vortex domain structures图 4 STO5/PTO20/STO5体系模拟结果 (a) STO/PTO/STO偶极子极化构型图, 为使图像更加清晰, 沿
$ \left[010\right] $ 方向每隔四层画一层, 红色箭头和蓝色箭头分别表示极化沿$ \left[001\right] $ 正方向和负方向; (b)图(a)中沿$ \left[010\right] $ 方向的第25层偶极子的纵切面投影图; (c)图(b)中黄色区域的放大图, 包括其沿$ \left[010\right] $ 方向的投影图和三维结构, 该涡旋畴沿$ \left[010\right] $ 方向呈柱状; (d)位于图(a)中沿$ [010] $ 方向的第13层偶极子纵切面底部; 在图(a)中, 图(b)—(d)表示结构的所在平面都用黑色箭头标注Figure 4. Simulation results of STO
$ _{5} $ /PTO$ _{20} $ / STO$ _{5} $ system: (a) STO/PTO/STO dipole polarization configuration, drawing every 5th layers along the$ \left[010\right] $ direction in order to make the figure clearer; (b) the projected dipole pololarization of the 25th layer along the$ \left[010\right] $ direction; (c) detailed view of the yellow area in panel (b), the vortex domain is cylindrical along the$ \left[010\right] $ direction; (d) the dipole configuration of the 13th layer along the$ [010] $ direction. The planes shown in panels (b)–(d) are marked by black arrows in panel (a)图 5 STO/PTO/STO模拟结果 (a) STO/PTO/STO偶极子极化构型示例(和图4(a)是同一个三维图), 沿[001]方向每隔三层画出一层中的偶极子, 同时画出PTO和STO交界层; (b) PTO中间层, 即沿
$ \left[001\right] $ 方向第15层的投影图; (c) PTO底部和STO交界层沿$ \left[001\right] $ 方向的投影图Figure 5. (a) Dipole polarization configuration of STO/PTO/STO system, which is the configuration as Fig. 4(a), drawing every 4th layer along the
$ \left[001\right] $ direction in order to make the figure clearer, the interface layer between PTO and STO is also shown; (b) projection of the middle PTO layer along the$ \left[001\right] $ direction; (c) projection of the bottom interface PTO layer along$[001]$ 表 1 图1(a)—(c)所示构型中的近似沿
$ \langle 111\rangle $ 方向的偶极子数目Table 1. The number of unit cell dipoles that are approximately along
$ \langle 111\rangle $ direction for the configurations in Figs. 1(a)–(c), respectively宏观极化方向 $ \left\langle 111\right\rangle $ $ \left\langle 110\right\rangle $ $ \left\langle 100\right\rangle $ 近似沿$ \left\langle 111\right\rangle $方向的偶极子极化个数 1728 1536 1348 近似沿$ \left\langle 111\right\rangle $方向的偶极子极化占比 100% 89% 78% -
[1] Seidel J, Ramesh R, Scott J, Catalan G 2012 Rev. Mod. Phys. 84 119Google Scholar
[2] Nataf G F, Guennou M, Gregg J M, Meier D, Hlinka J, Salje E K H, Kreisel J 2020 Nat. Rev. Phys. 2 634Google Scholar
[3] Scott J 2016 Ferroelectrics 503 117Google Scholar
[4] 谭丛兵, 钟向丽, 王金斌 2020 69 127702Google Scholar
Tan C B, Zhong X L, Wang J B 2020 Acta Phys. Sin. 69 127702Google Scholar
[5] Lu X M, Huang F Z, Zhu J S 2020 Acta Phys. Sin. 69Google Scholar
[6] Tian G, Yang W, Chen D, Fan Z, Hou Z, Alexe M, Gao X 2019 Natl. Sci. Rev. 6 684Google Scholar
[7] Das S, Tang Y L, Hong Z, Gonçalves M A P, McCarter M R, Klewe C, Nguyen K X, Gómez-Ortiz F, Shafer P, Arenholz E, Stoica V A, Hsu S L, Wang B, Ophus C, Liu J F, Nelson C T, Saremi S, Prasad B, Mei A B, Schlom D G, Íñiguez J, García-Fernández P, Muller D A, Chen L Q, Junquera J, Martin L W, Ramesh R 2019 Nature 568 368Google Scholar
[8] Huang J, Tan P, Wang F, Li B 2022 Crystals 12 786Google Scholar
[9] Das S, Ghosh A, McCarter M R, Hsu S L, Tang Y L, Damodaran A R, Ramesh R, Martin L W 2018 APL Mater. 6 100901Google Scholar
[10] Seidel J 2016 Topological Structures in Ferroic Materials, volume 228 of Springer Series in Materials Science (Cham: Springer International Publishing)
[11] Lu L, Nahas Y, Liu M, Du H, Jiang Z, Ren S, Wang D, Jin L, Prokhorenko S, Jia C L, Bellaiche L 2018 Phys. Rev. Lett. 120 177601Google Scholar
[12] Marton P, Rychetsky I, Hlinka J 2010 Phys. Rev. B 81 144125Google Scholar
[13] Zhong W, Vanderbilt D, Rabe K M 1995 Phys. Rev. B 52 6301Google Scholar
[14] Chen L Q 2008 J. Am. Ceram. Soc. 91 1835Google Scholar
[15] Liu J, Chen W, Wang B, Zheng Y 2014 Materials 7 6502Google Scholar
[16] Jiang Z, Xu B, Li F, Wang D, Jia C L 2015 Phys. Rev. B 91 014105Google Scholar
[17] Wojdeł J C, Íñiguez J 2014 Phys. Rev. Lett. 112 247603Google Scholar
[18] Zhong W, Vanderbilt D, Rabe K M 1994 Phys. Rev. Lett. 73 1861Google Scholar
[19] Liu Y, Tang Y L, Zhu Y L, Wang W Y, Ma X L 2016 Adv. Mater. Interfaces 3 1600342Google Scholar
[20] Brown W F 1941 Phys. Rev. 60 139Google Scholar
[21] Brown W F 1940 Phys. Rev. 58 736Google Scholar
[22] Aharoni A, et al. 2000 Introduction to the Theory of Ferromagnetism (Vol. 109) (Clarendon Press)
[23] García-Palacios J L, Lázaro F J 1998 Phys. Rev. B 58 14937Google Scholar
[24] Kumar D, Adeyeye A O 2017 J. Phys. D: Appl. Phys. 50 343001Google Scholar
[25] Zhu B, Lo C C H, Lee S J, Jiles D C 2001 J. Appl. Phys. 89 7009Google Scholar
[26] Shen K, Tatara G, Wu M W 2011 Phys. Rev. B 83 085203Google Scholar
[27] Li B L, Liu X P, Fang F, Zhu J L, Liu J M 2006 Phys. Rev. B 73 014107Google Scholar
[28] Wang D, Liu J, Zhang J, Raza S, Chen X, Jia C L 2019 Comput. Mater. Sci. 162 314Google Scholar
[29] Hong Z, Damodaran A R, Xue F, Hsu S L, Britson J, Yadav A K, Nelson C T, Wang J J, Scott J F, Martin L W, et al. 2017 Nano Lett. 17 2246Google Scholar
[30] Chen L Q, Zhao Y 2022 Prog. Mater. Sci. 124 100868Google Scholar
[31] Petković M, Petković L, Herceg ff 2010 J. Comput. Appl. Math. 233 1755Google Scholar
[32] Sugiura H, Hasegawa T 2019 J. Comput. Appl. Math. 358 136Google Scholar
[33] 黄云清, 舒适, 陈艳萍 2009 数值计算方法 (北京: 科学出版社)
Huang Y Q, Shu S, Chen Y P 2009 Numerical Computation Method (Beijing: Science Press) (in Chinese)
[34] Nishimatsu T, Iwamoto M, Kawazoe Y, Waghmare U V 2010 Phys. Rev. B 82 134106Google Scholar
[35] Nishimatsu T, Aoyagi K, Kiguchi T, J Konno T, Kawazoe Y, Funakubo H, Kumar A, V Waghmare U 2012 J. Phys. Soc. Jpn. 81 124702Google Scholar
[36] Nishimatsu T, Grünebohm A, Waghmare U V, Kubo M 2016 J. Phys. Soc. Jpn. 85 114714Google Scholar
[37] Kwei G H, Lawson A C, Billinge S J L, Cheong S W 1993 J. Phys. Chem. 97 2368Google Scholar
[38] Bersuker I B 1966 Phys. Lett. 20 589Google Scholar
[39] Ravel B, Stern E A, Vedrinskii R I, Kraizman V 1998 Ferroelectrics 206 407Google Scholar
[40] Liu J, Jin L, Jiang Z, Liu L, Himanen L, Wei J, Zhang N, Wang D, Jia C L 2018 J. Chem. Phys. 149 244122Google Scholar
[41] Tinte S, Stachiotti M, Sepliarsky M, Migoni R, Rodriguez C 1999 J. Phys.: Condens. Matter 11 9679Google Scholar
[42] Zhang F, Zhang J, Jing H, Li Z, Wang D, Jia C L 2022 Phys. Rev. B 105 024106Google Scholar
[43] Tadmor E, Waghmare U, Smith G, Kaxiras E 2002 Acta Mater. 50 2989Google Scholar
[44] Sani A, Hanfland M, Levy D 2002 J. Solid State Chem. 167 446Google Scholar
[45] Abid A Y, Sun Y, Hou X, Tan C, Zhong X, Zhu R, Chen H, Qu K, Li Y, Wu M, et al. 2021 Nat. Commun. 12 1Google Scholar
[46] Lemee N, Infante I C, Hubault C, Boulle A, Blanc N, Boudet N, Demange V, Karkut M G 2015 ACS Appl. Mater. Interfaces 7 19906Google Scholar
[47] Pereiro M, Yudin D, Chico J, Etz C, Eriksson O, Bergman A 2014 Nat. Commun. 5 1Google Scholar
[48] Sichuga D, Bellaiche L 2011 Phys. Rev. Lett. 106 196102Google Scholar
[49] Kim J H, Lange F F 1999 J. Mater. Res. 14 1626Google Scholar
[50] Zhang J, Wang Y J, Liu J, Xu J, Wang D, Wang L, Ma X L, Jia C L, Bellaiche L 2020 Phys. Rev. B 101 060103Google Scholar
Catalog
Metrics
- Abstract views: 4109
- PDF Downloads: 116
- Cited By: 0