Using classical Monte Carlo simulations, we study a simple statistical mechanical model of relevance to the emergence of polarization from local displacements on the square and cubic lattices. Our model contains two key ingredients: a Kitaev-like orientation-dependent interaction between nearest neighbors and a steric term that acts between next-nearest neighbors. Taken by themselves, each of these two ingredients is incapable of driving long-range symmetry breaking, despite the presence of a broad feature in the corresponding heat-capacity functions. Instead, each component results in a “hidden” transition on cooling to a manifold of degenerate states; the two manifolds are different in the sense that they reflect distinct types of local order. Remarkably, their intersection, i.e., the ground state when both interaction terms are included in the Hamiltonian, supports a spontaneous polarization. In this way, our study demonstrates how local-order mechanisms might be combined to break global inversion symmetry in a manner conceptually similar to that operating in the “hybrid” improper ferroelectrics. We discuss the relevance of our analysis to the emergence of spontaneous polarization in well-studied ferroelectrics such as ${\mathrm{BaTiO}}_3$ and ${\mathrm{KNbO}}_3$.

• Received 15 December 2017
• Revised 21 March 2018

DOI:https://doi.org/10.1103/PhysRevB.97.134106