Intervention Strategies for Polarization Switching in Hybrid Improper Ferroelectrics

TL;DR

The paper tackles polarization switching in hybrid improper ferroelectrics (HIFs) by fusing causal-discovery methods with first-principles simulations to identify material descriptors that govern the switching barrier $E_s$. Using LiNGAM-based causal discovery and physics-constrained interventions, it reveals that tolerance factor $\tau$, A-site radii mismatch, epitaxial strain, and octahedral rotation/tilt causally modulate $E_s$, and it uncovers three rotation–tilt pathways, including a cooperative mechanism where both distortions rise yet the barrier drops under strain. The authors validate these causal mechanisms with DFT and NEB calculations across a broad set of A-site layered double perovskites and superlattices, showing that substrate-induced strain (e.g., NdScO$_3$, NdGaO$_3$) can selectively realize the predicted distortion patterns to reduce $E_s$ by tens of percent. The study provides actionable, mechanism-driven design rules linking composition, structure, and substrate choice to low-barrier HIFs, and demonstrates a general causal-reasoning workflow that can guide rational design in complex oxides beyond the specific systems studied.

Abstract

The potential of hybrid improper ferroelectrics (HIFs) in electronic and spintronic devices hinges on their ability to switch polarization. Although the coupling between octahedral rotation and tilt is well established, the factors that govern switching barriers remain elusive. In this study, we explore this area to demonstrate the critical role of causal reasoning in uncovering the mechanisms to control the ferroelectric switching barrier in HIFs. By combining causal discovery, causal interventions, and first-principles simulations, we identify tolerance factor, A-site cation radii mismatch, epitaxial strain, and octahedral rotation/tilt as key parameters and quantify how their interplay directly influences switching barrier. Three key insights emerge from our work: (a) the analysis identifies the structural descriptors controlling polarization reversal across a broad family of A-site-layered double perovskites and superlattices, (b) it uncovers non-trivial, material-specific rotation-tilt mechanisms, including a counterintuitive cooperative pathway where both rotation and tilt change while lowering the barrier, an effect mostly inaccessible to conventional Landau or first-principles-based approaches and (c) it maps these material-specific mechanisms to experimentally realizable parameters, showing that epitaxial strain from orthorhombic substrates (e.g., NdScO$_3$, NdGaO$_3$) selectively tunes octahedral distortions to achieve barrier reduction across varied compositions. These results establish actionable, materials-by-design principles linking composition, structure, and strain to polarization switching, while highlighting the potential of causal reasoning to guide intelligent, mechanism-driven strategies for engineering complex functional oxides.

Figures (6)

• Figure 1: Outline of how causal interventions can elucidate the fundamental atomistic mechanisms governing ferroelectric switching in hybrid improper ferroelectric oxides.
• Figure 2: Key steps in integrating causal modeling with physics simulations to derive mechanisms for reducing switching barriers. The approach employs algorithms for causal modeling to address design tasks effectively, followed by causal interventions that incorporate physics-based constraints to gain insights into the underlying processes guiding the phenomena.
• Figure 3: Schematic representation of the (a) high symmetry phase, key structural modes such as in-phase rotation, tilt in superlattices (upper panel), double perovskite oxides (middle panel) and low symmetry phase (bottom panel). (b) The structural distortions corresponding to polarization switching via out-of-phase rotation of the BO$_6$ octahedra. (c) Illustration of the compositional space showcasing the chemical diversity of cations, as considered in the study.
• Figure 4: (a) Directed acyclic graph illustrating the cause-effect relations between various features and the target. (b) Two-dimensional contour representation of the changes in rotation and tilt that lead to changes in the switching barrier, as estimated by the causal interventions conducted across all materials with $\tau$$>$0.8 (Group I in our study).All materials exhibit negative changes, indicating a reduction in the E$_\text{s}$. Upon intervention, three distinct response pathways are identified: (1) concurrent and complementary variations in the rotation and tilt angles, (2) an increase in $\theta_{t}$ with $\theta_{r}$ remaining approximately constant, and (3) a decrease in $\theta_{r}$ with $\theta_{t}$ remaining approximately constant. These pathways correspond to distinct distortion mechanisms that reduce $E_{\text{s}}$. The materials following each pathway are denoted in green, blue, and red, respectively.
• Figure 5: Selected Group I materials ($\tau > 0.8$) for which complementary changes in rotation and tilt angles reduce the switching barrier. The required distortions are identified via causal interventions, while DFT calculations provide validation of these material-specific pathways.