Polar morphologies from first principles: PbTiO${}_3$ films on SrTiO${}_3$ substrates and the $p(2\times \Lambda )$ surface reconstruction

Abstract

Low dimensional structures comprised of ferroelectric (FE) PbTiO${}_3$ (PTO) and quantum paraelectric SrTiO${}_3$ (STO) are hosts to complex polarization textures such as polar waves, flux-closure domains and polar skyrmion phases. Density functional theory (DFT) simulations can provide insight into this order, but, are limited by the computational effort needed to simulate the thousands of required atoms. To relieve this issue, we use the novel multi-site support function (MSSF) method within DFT to reduce the solution time for the electronic groundstate whilst preserving high accuracy. Using MSSFs, we simulate thin PTO films on STO substrates with system sizes $>2000$ atoms. In the ultrathin limit, the polar wave texture with cylindrical chiral bubbles emerges as an intermediate phase between full flux closure domains and in-plane polarization. This is driven by an internal bias field born of the compositionally broken inversion symmetry in the [001] direction. Since the exact nature of this bias field depends sensitively on the film boundary conditions, this informs a new principle of design for manipulating chiral order on the nanoscale through the careful choice of substrate, surface termination or use of overlayers. Antiferrodistortive (AFD) order locally interacts with these polar textures giving rise to strong FE/AFD coupling at the PbO terminated surface driving a $p(2\times \Lambda )$ surface reconstruction. This offers another pathway for the local control of ferroelectricity.

13 pages, 7 figures

References in corpus (6)

• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Large scale and linear scaling DFT with the CONQUEST code
• Efficient Calculations with Multisite Local Orbitals in the large-scale DFT Code CONQUEST
• Highly accurate local basis sets for large-scale DFT calculations in CONQUEST
• First principles soft mode lattice dynamics of PbZr${}_{0.5}$Ti${}_{0.5}$O${}_3$ and shortcomings of the virtual crystal approximation
• The pseudoatomic orbital basis: electronic accuracy and soft-mode distortions in ABO${}_3$ perovskites

Cited by in corpus (6)

• Roadmap on Electronic Structure Codes in the Exascale Era
• Two-dimensional ferroelectricity induced by octahedral rotation distortion in perovskite oxides
• Giant voltage amplification from incipient ferroelectric states
• Assessing the nature of nanoscale ferroelectric domain walls in lead titanate multilayers
• Large-scale DFT methods for calculations of materials with complex structures
• Origin of ferroelectric domain wall alignment with surface trenches in ultrathin films