---
title: 'Stacking Ferroelectricity: Interlayer Registry Effects'
url: https://www.emergentmind.com/topics/stacking-ferroelectricity
type: topic
---

# Stacking Ferroelectricity: Interlayer Registry Effects

Stacking ferroelectricity is a form of ferroelectric order in which spontaneous polarization is created by **interlayer registry** rather than by the intrinsic polar instability of an isolated monolayer or by the conventional bulk displacive mechanism. In this framework, relative translation, rotation, reversed stacking, twist, or more general heterodeformation can remove the symmetries that enforce vanishing dipole moment, so that bilayers, multilayers, moiré superlattices, and even trilayer sandwiches assembled from nonpolar constituents become polar and electrically switchable [2210.16542], [2309.00728], [2509.00384]. The subject now spans interfacial sliding ferroelectrics in hexagonal boron nitride and rhombohedral transition-metal dichalcogenides, mixed-stacking graphene, vertical trilayer ferroelectrics, stacking-controlled multiferroics, and quasi-one-dimensional van der Waals crystals whose macroscopic polarity is fixed by the stacking of polar subunits [2010.06600], [2108.07659], [2305.10896], [2307.11267], [2209.13070].

## 1. Definition and distinguishing features

The general bilayer formulation treats the second layer as a stacked image of the first,
\[
S'=\hat t_z \hat O S,
\]
where \(\hat t_z\) is an out-of-plane translation and \(\hat O=\{O\mid \mathbf T_O\}\) is the essential stacking operator with rotational part \(O\) and in-plane translation \(\mathbf T_O\) [2210.16542]. In that sense, stacking ferroelectricity is a symmetry problem posed on the *assembled bilayer or multilayer*, not on the monolayer alone. The relevant order parameter is therefore a stacking coordinate—most commonly a relative sliding field, but in other realizations a layer flip, a relative rotation, or a vertical off-centering of a sandwiched layer.

This distinguishes stacking ferroelectricity from conventional displacive ferroelectricity. In perovskite-like systems, polarization is associated with ionic off-centering within a three-dimensional unit cell. In stacking ferroelectrics, by contrast, polarization is “rooted in the interlayer charge compensation and the resultant dipole formation between adjacent layers,” and ferroelectric behavior can be created “even from non-ferroelectric monolayers by breaking the inversion symmetry through specific stacking configurations” [2509.00384]. The underlying monolayers may be nonpolar, centrosymmetric, or only weakly polar; what matters is whether the assembled stack retains or destroys the symmetries that forbid a net dipole.

Sliding ferroelectricity is the most prominent realization of this idea, but it is not the whole subject. The general bilayer-stacking-ferroelectricity framework includes pure translation, rotation-related stacking, mirrored or reversed bilayers, twisted bilayers, and combined in-plane and out-of-plane polar states [2210.16542]. A distinct but related branch is **vertical stacking ferroelectricity**, in which an A/B/A trilayer acquires bistable out-of-plane polarization through the vertical off-centering of the middle layer rather than through in-plane sliding [2307.11267].

## 2. Symmetry principles and microscopic theory

A systematic group-theory analysis over all **80 layer groups** classifies monolayers into 15 IP, 13 OP, 4 CP, and 48 NP layer groups and then determines which bilayer stackings can become ferroelectric [2210.16542]. The central symmetry statement is that layer-exchanging operations \(R^{-}\), such as inversion, \(m_z\), \(C_{2\alpha}\), or \(S_{nz}\), forbid the corresponding polarization component if they survive in the bilayer. A particularly important rule is that **pure sliding cannot break monolayer inversion for a centrosymmetric monolayer**; therefore sliding ferroelectricity is impossible in that case unless the stacking operation includes a nontrivial rotational part [2210.16542]. This is why general bilayer stacking ferroelectricity is broader than sliding ferroelectricity.

Beyond symmetry selection rules, a quantum-geometric description has been developed for out-of-plane stacking ferroelectricity in van der Waals bilayers. In that approach, the bilayer at fixed in-plane momentum is mapped to the two-cell limit of an SSH/Rice–Mele chain, so that the vertical polarization is controlled by a Berry phase,
\[
P = \frac{e}{2\pi} \oint \langle u_{-}(k)|-i\partial_k|u_{-}(k)\rangle\, dk,
\]
with robust polarity when the broken-sublattice-symmetry scale and the asymmetric interlayer hopping are comparable, summarized by the criterion \(\Delta_{AB}\sim g\) [2309.00728]. This framework unifies AB-stacked honeycomb bilayers, 3R bilayer TMDs, and \(T_d\) bilayer TMDs, and it explains why SiC-type bilayers are strongly polar, 3R-TMDs are intermediate, and \(T_d\)-TMDs are weak and cancellation-prone [2309.00728].

A complementary microscopic route appears in mixed-stacking graphene. There the polarization is computed directly from layer charge redistribution,
\[
P_z=ed \sum_{L=1}^{n+m+3} L\, n_L,
\]
using a self-consistent full hybrid \(\mathbf{k}\cdot\mathbf{p}\)-tight-binding Slonczewski–Weiss–McClure Hamiltonian [2305.10896]. Structural asymmetry alone is not sufficient: a reduced electron-hole-symmetric model still gives \(P_z=0\), so a nonzero equilibrium dipole requires both broken inversion and electron-hole-symmetry-breaking terms such as \(v_4\), \(\Delta'\), \(\gamma_2\), and \(\gamma_5\) [2305.10896]. Screening is essential; in ABCB tetralayer graphene, the unscreened \(P_z^u \approx 0.50\, e/\mu\text{m}\) is reduced to a screened \(P_z \approx -0.07\, e/\mu\text{m}\), and the sign can even reverse [2305.10896].

These theories converge on a common picture. Stacking ferroelectricity is produced when stacking lowers symmetry so that layer-resolved charge, orbital hybridization, or effective Wannier-center positions become inequivalent between the two sides of the bilayer. The resulting dipole may be predominantly electronic, mixed ionic-electronic, or coupled to a structural relaxation of the stacking coordinate, depending on the material class.

## 3. Bilayer and interfacial realizations

Parallel-stacked bilayer h-BN is the canonical sliding-ferroelectric realization. In experimentally assembled \(0^\circ\) bilayers, the ferroelectric states are AB and BA, while the natural bulk-like \(180^\circ\) AA′ stacking is nonpolar [2010.06600]. Transport measurements with an adjacent graphene sensor yield
\[
P_{2D} = 2.25(0.37)\times 10^{-12}\ \mathrm{C\,m^{-1}},
\]
equivalent to
\[
P_{3D} = 0.68\ \mu \mathrm{C\,cm^{-2}},
\]
in good agreement with the Berry-phase value \(P_{2D,\mathrm{theory}} = 2.08\times 10^{-12}\ \mathrm{C\,m^{-1}}\) [2010.06600]. The switching fields extracted from dual-gate maps are approximately \(-0.06\ \mathrm{V\,nm^{-1}}\) for BA \(\rightarrow\) AB and \(0.16\ \mathrm{V\,nm^{-1}}\) for AB \(\rightarrow\) BA, while room-temperature writing is achieved with \(V_B=+1.8\ \mathrm{V}\) \((+0.19\ \mathrm{V\,nm^{-1}})\) and \(V_B=-1.2\ \mathrm{V}\) \((-0.13\ \mathrm{V\,nm^{-1}})\); retention remains essentially unchanged after 14 and 31 days at \(V_B=0\) [2010.06600].

A direct interfacial view of the same mechanism was obtained in naturally grown h-BN flakes stacked in a metastable parallel orientation. Kelvin probe force microscopy resolves alternating AB and BA domains with
\[
\Delta V_{\mathrm{KP}} = 100 \pm 10~\text{mV},
\]
corresponding to
\[
P_z \approx 0.66~\text{Debye}/\text{nm}^2,
\]
and reversible switching occurs by biased-tip-driven domain-wall motion coupled to lateral sliding [2010.05182]. The required relative translation is \(1.44~\text{\AA}\), the domain-wall width is about \(10~\text{nm}\), and domain-wall motion is observed for fields exceeding approximately \(0.3~\text{V/nm}\) [2010.05182]. The microscopic origin was traced to a subtle interplay between charge redistribution and a tiny out-of-plane displacement \(\Delta d \sim 2\times10^{-3}~\text{\AA}\) [2010.05182].

Rhombohedral-stacked bilayer TMDs extend the same mechanism to semiconducting \(MX_2\) bilayers. In bilayer WSe\(_2\), MoSe\(_2\), WS\(_2\), and MoS\(_2\), parallel stacking produces two local minima, MX and XM, with opposite out-of-plane dipoles generated by asymmetric interlayer hybridization and charge redistribution [2108.07659]. Nearly parallel stacked bilayers reconstruct into triangular MX/XM moiré ferroelectric domains visualized by PFM, and electric-field-induced domain-wall motion is directly observed [2108.07659]. In graphene-sensor devices, the built-in interlayer potential is about \(55\)–\(57\) mV across all four TMDs, while a WSe\(_2\) device exhibits a coercive field of about \(0.30\ \text{V nm}^{-1}\) [2108.07659].

A growth-engineered route was demonstrated in bilayer MoS\(_2\) homoepitaxy with co-existing 3R polytypic domains. Under low-Mo, low-temperature, sulfur-rich two-step CVD, fan-shaped bilayer patterns form with alternating AA, AB, and BA domains separated by partial dislocations [2205.12118]. Since AB and BA correspond to opposite out-of-plane interlayer ferroelectric polarizations, their coexistence provides a structural basis for room-temperature switching. Back-gated FETs with these 3R polytypic domains show repeatable counterclockwise hysteresis, and the memory window exceeds that of compact-shaped 3R bilayer control devices [2205.12118].

## 4. Multilayer, moiré, vertical, and elemental variants

Graphene broadened stacking ferroelectricity beyond binary compounds. Mixed-stacking few-layer graphene is predicted to host **weak ferroelectricity** when the stacking sequence breaks inversion and \(z\to -z\) symmetry, as in ABCB and ABAC tetralayers or in \(n\)ABA\(m\) films with an asymmetrically placed twin boundary [2305.10896]. For ABCB tetralayer graphene, the self-consistent screened polarization is \(P_z \approx -0.07\, e/\mu\text{m}\), and in reconstructed small-angle structures the largest local polarization occurs at domain walls and especially near domain-wall intersections [2305.10896]. A distinct across-layer mechanism in pure multilayer graphene appears for \(N>3\): tetralayer ABAC has \(P_z=0.21\ \mathrm{pC/m}\), five-layer ABABC and ABACB have \(P_z=0.17\) and \(0.32\ \mathrm{pC/m}\), and six-layer polar states span approximately \(0.05\) to \(1.0\ \mathrm{pC/m}\); the polar and nonpolar states are nearly degenerate, within \(0.3\)–\(0.4\ \mathrm{meV/unitcell}\), and the switching barriers are below \(5\ \mathrm{meV/unitcell}\) [2305.07286].

The same subject has now been accessed experimentally without moiré alignment. Dual-gated, non-aligned ABCB tetralayer graphene encapsulated in hBN exhibits pronounced hysteresis under both top and bottom gate modulation, with
\[
\Delta n_t = 0.63 \times 10^{12}\ \mathrm{cm}^{-2},\qquad
\Delta n_b = 0.7 \times 10^{12}\ \mathrm{cm}^{-2},
\]
and an inferred polarization
\[
P_{2D}=0.04\ \mu \mathrm{C}\,\mathrm{cm}^{-2},
\]
persisting up to room temperature [2504.07935]. The proposed microscopic origin is reversible layer-polarized charge reordering driven by gate-induced transitions between ABCB and BCBA-related stackings, without requiring a moiré superlattice [2504.07935].

A separate extension is **vertical stacking ferroelectricity** in A/B/A trilayers. In this model, the two outer layers are fixed symmetrically while the middle layer is free to move vertically, so the energy can be written as
\[
E=F\left(A,B,d_{\text{device}},d_{AB}\right),
\]
with a single-well centered state for \(d_{\text{device}} \le 2d_{\mathrm{vdw}}\) and a double well of two symmetry-related off-center states for \(d_{\text{device}} \gg 2d_{\mathrm{vdw}}\) [2307.11267]. First-principles calculations give \(P_s=1.83\ \mathrm{pC/m}\) and a \(10\ \mathrm{meV/atom}\) switching barrier for \(h\)-BN/\(h\)-BN/\(h\)-BN, and \(P_s=1.35\ \mathrm{pC/m}\) with an \(8\ \mathrm{meV/atom}\) barrier for \(h\)-BN/graphene/\(h\)-BN; in both cases the crossover to polar stability occurs around \(d_{\text{device}}=7\ \text{\AA}\) [2307.11267].

Large-heterodeformation h-BN shows that stacking ferroelectricity is not confined to the near-AA moiré limit. Bicrystallography and a DFT-informed BFIM continuum model demonstrate ferroelectricity not only in AA-vicinal systems but also in configurations vicinal to the large-angle \(\Sigma 7\) parent at \(21.786789^\circ\) [2510.01419]. In the \(\Sigma 7\)-vicinal case, the two degenerate GSFE minima carry opposite out-of-plane polarization with
\[
|P_z| = 1.35\times 10^{-4}\,\text{C/m}^2,
\]
and the relevant interface dislocations have Burgers vector magnitude \(\approx 0.55\,\text{\AA}\), much smaller than in small-twist AA-vicinal h-BN [2510.01419].

Stacking engineering also reaches elemental topological ferroelectrics and polar metals. In buckled-honeycomb Bi films, 2BL ABCA is nonpolar and centrosymmetric, whereas 2BL ABAC is polar with layer group \(p3m1\), and CABA is its inversion partner with opposite polarization; the two-bilayer Bi film with a polar stacking sequence is identified as an elemental topological ferroelectric, while 3BL and 4BL polar structures are elemental polar metals with topological nontrivial edge states [2309.14609].

## 5. Domains, switching kinetics, and field control

Because the order parameter is a stacking coordinate, switching commonly proceeds through **domain-wall motion** rather than coherent uniform reversal. A first-principles and machine-learning study of bilayer h-BN showed that AB and BA carry
\[
P_z \approx \pm 1.46 \times 10^{-12}\ \mathrm{C/m},
\]
while the saddle-point state along the sliding path has vanishing \(P_z\) and a large in-plane polarization of \(1.38 \times 10^{-1}\ \mathrm{C/m}\) [2212.14203]. The domain walls are unusually wide, with widths of \(9.7\), \(17.6\), \(32.1\), and \(40.7\) nm for \(0^\circ\), \(30^\circ\), \(60^\circ\), and \(90^\circ\) walls, respectively, described by
\[
w = u_0 \sqrt{\frac{2\lambda_{1D}}{\Delta}},
\]
where \(\Delta\) is the switching barrier between stacking-polarization states and \(\lambda_{1D}\) is the effective in-plane elastic stiffness [2212.14203]. Coherent monodomain reversal requires \(E_v^c \approx 1.41\ \mathrm{V/nm}\) at 100 K and \(1.32\ \mathrm{V/nm}\) at 300 K, but domain-wall motion lowers the threshold to \(0.026\ \mathrm{V/nm}\) for a \(0^\circ\) wall and \(0.11\ \mathrm{V/nm}\) for a \(90^\circ\) wall; the predicted wall velocity is about \(6000\ \mathrm{m/s}\), enabling switching in about \(9\)–\(15\) ps [2212.14203]. In ideal twisted h-BN, the same analysis yields a field-responsive but nearly reversible state termed **super-paraelectric** [2212.14203].

Continuum theory for moiré superlattices clarifies this distinction. In 3R MoS\(_2\), local stackings carry local spontaneous polarization \(p_0(\mathbf r)\), and the relaxed moiré divides into AB and BA domains separated by narrow walls because of the competition between stacking and elastic energies [2108.10838]. An applied field enters the local energy density through
\[
V_{\text{tot}} = V_{\text{elastic}} + V_{\text{stack}} + V_{\text{elec}},
\]
with a linear \(-Ep_0(\mathbf r)\) term that splits AB and BA and a quadratic dielectric term that reduces the stacking energy and softens the domains [2108.10838]. The resulting net polarization of the moiré superlattice arises from unequal domain areas rather than from flipping the sign of local polarization inside each domain. The authors therefore argue that ideal moiré superlattices governed by this mechanism are generally **not truly ferroelectric** in the strict zero-field-switchable sense, even though they contain local polar domains and exhibit field-tunable average polarization [2108.10838]. This is a central terminology dispute in the field.

Photoexcitation adds another control channel. In 3R bilayer MoS\(_2\), the dark-state polarization is \(P_z=0.82\ \mathrm{pC/m}\), but constrained-DFT calculations show a non-monotonic photoinduced evolution: \(P_z=0.47\ \mathrm{pC/m}\) at \(n_{\rm ph}=0.01\ e/{\rm u.c.}\), \(0.03\ \mathrm{pC/m}\) at \(0.2\ e/{\rm u.c.}\), and \(0.4\ \mathrm{pC/m}\) at \(0.8\ e/{\rm u.c.}\), while the sliding barrier remains in the \(17\)–\(21\ \mathrm{meV/u.c.}\) range [2410.00776]. At \(0.9\)–\(1.0\ e/{\rm u.c.}\), the original \(P3m1\) structure becomes unstable and reconstructs into a \(2\times2\) \(Cm\) phase with \(P_z\) enhanced to about \(\pm 4.5\) to \(\pm 5.4\ \mathrm{pC/m}\) [2410.00776]. This underscores that stacking ferroelectricity is often electronically soft because its origin is interfacial charge transfer rather than a large ionic displacement.

## 6. Multiferroic extensions, non-bilayer analogues, and open problems

Stacking ferroelectricity has increasingly merged with multiferroicity. In bilayer ScI\(_2\), aligned AB and BA are polar with opposite out-of-plane polarization
\[
P = 0.18 \times 10^{-12}\ \text{C/m},
\]
while AA is nonpolar because it preserves \(M_z\); AB and BA are also the lowest-energy aligned stackings, and sliding between them requires an \(8.87\ \text{meV/f.u.}\) barrier [2510.16379]. These same stackings are interlayer ferromagnetic and, with SOC and out-of-plane spin orientation, exhibit spontaneous valley polarization with splitting around \(101.0\ \text{meV}\) [2510.16379]. In reversed bilayer PtBr\(_3\), AC′ stacking yields a polar antiferromagnet with
\[
P_{\rm OOP}=2.4\ \text{pC/m},\qquad P_{\rm IP}=9.7\ \text{pC/m},
\]
and a CI-NEB switching barrier of \(18\ \text{meV/f.u.}\); the ferroelectric polarization controls the spin splitting and the sign of the Kerr signal, enabling electrical writing and magneto-optical readout without net ferromagnetism [2411.13182]. More generally, bilayer stacking-ferroelectricity theory predicts that even centrosymmetric 2D ferromagnets can become multiferroic after appropriate stacking, and that in some cases the out-of-plane polarization is interlocked with an in-plane component, allowing deterministic control of \(P_z\) by an in-plane electric field [2210.16542].

The stacking principle is not confined to two-dimensional bilayers. In quasi-one-dimensional van der Waals oxyhalides NbOX\(_3\), the local Nb off-centering creates polar double chains, and different ways of stacking those double chains generate ferroelectric or antiferroelectric bulk phases with meV-scale energy differences [2209.13070]. The reported spontaneous polarizations of the ferroelectric phases range from \(14.2\) to \(27.2\ \mu\text{C/cm}^2\), and \(d_{33}\) reaches \(60.1\ \text{pC/N}\) in monoclinic NbOI\(_3\) [2209.13070]. This suggests that stacking ferroelectricity is best understood as a broader **registry-controlled polarity** principle, not solely as a van der Waals bilayer phenomenon.

Several limitations remain. Many proposed systems are still intrinsically theoretical, often based on semilocal DFT with empirical vdW corrections, and frequently omit finite-temperature fluctuations, disorder, realistic contacts, substrate effects, and depolarization screening [2307.11267]. Even in experimentally established systems, the microscopic switching path is often inferred from transport or PFM rather than directly imaged under bias, and quantitative memory metrics such as coercive field, endurance, or retention remain incomplete in several material families [2205.12118], [2504.07935]. The terminology of “ferroelectricity” in moiré systems also remains contested, because domain-redistribution-induced net polarization can mimic ferroelectric response without satisfying the strict criterion of a switchable spontaneous polarization at zero field [2108.10838].

Taken together, the literature defines stacking ferroelectricity as a distinct ferroic paradigm in which **interlayer geometry is itself the order parameter**. It can be realized by lateral sliding between AB and BA, by across-layer symmetry breaking in multilayer graphene, by vertical off-centering in A/B/A trilayers, by bicrystallographic domain formation in twisted or strained bilayers, or by stacking polar subunits into ferroelectric and antiferroelectric manifolds. The unifying statement is that polarization is created, reversed, or annihilated by changing how layers or chains are registered in space.

Source: https://www.emergentmind.com/topics/stacking-ferroelectricity