---
title: 'Sliding Ferroelectrics: Interlayer Charge Switching'
url: https://www.emergentmind.com/topics/sliding-ferroelectrics
type: topic
---

# Sliding Ferroelectrics: Interlayer Charge Switching

Sliding ferroelectrics are layered, usually van der Waals, systems in which a spontaneous out-of-plane polarization is generated by stacking-dependent interfacial charge transfer and reversed by in-plane relative sliding of the constituent layers rather than by the conventional intralayer ion-displacement mechanism. The class now spans bilayer and multilayer inorganic materials such as \(h\)-BN, rhombohedral \(MX_2\), MBi\(_2\)Te\(_4\), CuF\(_2\), and Fe\(_5\)GeTe\(_2\), as well as an amphidynamic coordination polymer, \((15\text{-Crown-5})\)Cd\(_3\)Cl\(_6\), in which direct macroscopic \(P\)–\(E\) hysteresis was reported [2505.09084] [2211.05277] [2506.09317] [2302.03414].

## 1. Structural origin and polarization formalism

The defining structural motif is a stack of nonpolar or weakly polar monolayers whose relative registry removes inversion or mirror symmetry and creates a vertical dipole. In bilayer MBi\(_2\)Te\(_4\) \((M=\mathrm{Ge},\mathrm{Sn},\mathrm{Pb})\), the centrosymmetric AA′ bilayer is unstable and spontaneously slides into two equivalent symmetry-broken minima, AB′\(_1\) and AB′\(_2\), carrying \(P_\downarrow\) and \(P_\uparrow\), respectively; the two states are reversibly connected by in-plane sliding along either \([100]\) or \([110]\) [2506.09317]. In 3R-bilayer MoS\(_2\), two energy-degenerate stackings, XM and MX, are related by a shear translation \(\Delta x \approx 1.23\,\text{\AA}\), and the interface dipole reverses when that lateral displacement is executed [2401.16150]. In Janus TMD bilayers \(XMY\), the sliding-induced interlayer polarization coexists with an intrinsic monolayer electric field set by the electronegativity difference between the two chalcogens, so the monolayer polarity itself becomes a control parameter for the bilayer ferroelectric response [2505.22207].

The polarization is typically evaluated with the modern Berry-phase formalism. In MBi\(_2\)Te\(_4\), the polarization is written as
$$
\mathbf{P}=\frac{e}{\Omega}\sum_{n}^{\rm occ}\int_{\rm BZ}\frac{d\mathbf{k}}{(2\pi)^3}\,\langle u_{n\mathbf{k}}|\,i\nabla_{\mathbf{k}}\,|u_{n\mathbf{k}}\rangle,
$$
and, for the sliding path from AA′ to a polar AB′ state,
$$
P_z=\frac{e}{2\pi A_{\rm cell}}\,\Delta\varphi.
$$
The same systems admit an equivalent charge-density picture in which
$$
\Delta\rho(z)=\rho_{\rm AB'_i}(z)-\rho_{\rm AA'}(z), \qquad
P_z=\int dz\, z\,\Delta\rho(z),
$$
making explicit that the dipole originates from layer-resolved charge redistribution across the vdW gap [2506.09317]. Closely related Berry-phase and charge-integration forms are used in bilayer MoS\(_2\), Janus TMDs, and \(\beta\)-ZrI\(_2\) [2410.00776] [2505.22207] [2510.03738].

Representative computed values show that the scale of polarization is system dependent rather than universal. For MBi\(_2\)Te\(_4\), \(|P_z|=0.571\), \(0.615\), and \(0.623\) pC/m for GeBi\(_2\)Te\(_4\), SnBi\(_2\)Te\(_4\), and PbBi\(_2\)Te\(_4\), respectively [2506.09317]. Bilayer CuF\(_2\) reaches \(P_z=\pm 1.23\) pC/m [2603.10907], whereas \(\beta\)-ZrI\(_2\) is reported as \(P_s=0.30\,\mu\text{C}/\text{cm}^2\) [2510.03738]. This spread suggests that sliding ferroelectricity is best viewed as a mechanism class rather than a single materials family with one characteristic polarization scale.

## 2. Symmetry, energy landscapes, and switching coordinates

The symmetry reduction induced by sliding controls which polarization components are allowed and how electric fields couple to the sliding coordinate. In rhombohedral \(MX_2\) bilayers, the two ferroelectric ground states M and N have point group \(C_{3v}\), which allows only \(P_z\), while the midpoint A along the sliding path has \(C_{2v}\), which forbids \(P_z\) but allows an in-plane \(P_y\). For a generic intermediate registry of point group \(C_s\), both \(P_z\) and \(P_y\) are permitted, while \(P_x\) remains forbidden for sliding along \(y\). For MoS\(_2\), the calculated evolution is \(P_z(\lambda\approx 0)=+0.77\) pC/m, \(P_z(\lambda\approx 1/2)=0\), \(P_z(\lambda\approx 1)=-0.77\) pC/m, with \(P_y(\lambda\approx 1/2)\approx 2.5\) pC/m [2308.15241]. The existence of a large in-plane dipole at the intermediate state is a recurring feature of sliding pathways.

In \(h\)-BN bilayers, the coupling between a vertical electric field and lateral motion is governed by the off-diagonal Born effective charge tensor. At AB or BA, \(Z^*_{xz}=Z^*_{yz}=0\) by \(C_3\) symmetry, but along the AB\(\rightarrow\)SP\(\rightarrow\)BA path the off-diagonal components become sharply nonzero and peak near the saddle point, so an out-of-plane field \(E_z\) generates in-plane forces \(F_x=Z^*_{xz}E_z\) and \(F_y=Z^*_{yz}E_z\) that drive the required sliding [2505.09084]. This is a symmetry-mediated electromechanical coupling rather than the direct longitudinal force familiar from oxide ferroelectrics.

The sliding potential is commonly represented by a double well. For BN, a fitted one-dimensional form is
$$
U_0(t)=e_4 t^4 - 2e_4 t_{\rm BA}^2 t^2 + e_0,
$$
with \(t_{\rm BA}\simeq -0.725\,\text{\AA}\), \(e_4=27.017\,\text{meV}/\text{\AA}^4\), and \(e_0=7.41\,\text{meV}\) [2407.15081]. Because the saddle state also carries a finite in-plane dipole, an inclined electric field lowers the switching threshold more efficiently than a purely vertical field. The critical field is
$$
E_{t,c}(\theta)=\frac{\Delta U_0}{P_\perp \cos\theta + P_\parallel \sin\theta},
$$
and, for BN, it is minimized at \(\theta \simeq 24^\circ\) below the normal, giving a reduction of over \(70\%\) relative to purely vertical switching [2407.15081]. This result is presented as universal for sliding ferroelectric systems whose monolayers belong to the symmetry group \(p\bar{6}m2\).

## 3. Domain walls, nonuniform switching, and superlubric motion

A central controversy concerns whether polarization reversal proceeds by homogeneous rigid translation of an entire atomic layer or by domain-wall-mediated collective motion. For Bernal \(h\)-BN bilayers, the synchronized long-distance translation picture has been explicitly challenged: it is described as a “spinodal decomposition-like homogeneous switching process” that violates Neumann’s principle and is unlikely to occur because, in a pure AB or BA domain, \(C_3\) symmetry forces the unit-cell-averaged in-plane force under \(E_z\) to vanish [2502.01007]. In the same system, any perturbation breaking the in-plane symmetry is reported to trigger avalanche-like switching, so the practical coercive field is governed by fluctuations and symmetry breaking rather than by a conventional intrinsic maximum of the restoring-force curve [2505.09084].

The alternative picture centers on broad domain walls. In \(h\)-BN, a “wave-like” wall of width \(\sim 10\) nm is described by an order parameter
$$
\lambda(x)=2u_y(x)/u_0-1,
$$
with free-energy density
$$
f(\lambda,\partial_x\lambda)=\tfrac12 A(\partial_x\lambda)^2+\tfrac14 B(1-\lambda^2)^2,
$$
equilibrium width \(w_0=2\sqrt{A/B}\), and profile \(\lambda(x)=\tanh(2x/w_0)\) [2502.01007]. The key claim is that coherent propagation of this wide interface has an almost zero activation barrier because the local energy increase at the leading edge is balanced by the local decrease at the trailing edge. In the same material, the off-diagonal Born charges peak near the \(P=0\) region inside the wall, producing a nonuniform driving force and a “wriggling” wall motion under \(E_z\) [2505.09084].

The resulting dynamics are superlubric or soliton-like rather than creep-like. Finite-field MD for \(h\)-BN gives \(v \approx 3\times 10^3\) m/s at \(T=300\) K and \(E_3=0.3\) V/nm for the \(\Sigma_0\) wall, with extrapolated velocities up to \(\approx 4\times 10^3\) m/s and an anomalous cooling-promoted increase of velocity [2502.01007]. In bilayer 3R-MoS\(_2\), the domain wall width is \(\xi\approx 100\,\text{\AA}\), the characteristic speed is \(V_c \simeq 3.0\times 10^3\) m/s, and the translational coordinate obeys \(M\ddot X=qE\), yielding uniformly accelerated motion under field and constant coasting velocity after field removal [2502.02137]. Machine-learning molecular dynamics on ferroelectric MoS\(_2\) moiré superlattices sharpen this distinction further: the rigid-sliding barrier is \(\approx 3.8\) meV/atom, whereas the domain-wall-mediated barrier is \(<0.05\) meV/atom, sufficient to explain thermally driven interlayer drift of order \(1\) m/s at 300 K; sulfur vacancies of about \(0.1\%\) already induce a sliding-to-pinning transition [2604.20277].

## 4. Thickness dependence and multistate ferroelectricity

Layer number introduces an additional discrete degree of freedom absent in ordinary bilayer descriptions. In multilayer 3R-MoS\(_2\), each MoS\(_2\)/MoS\(_2\) interface is modeled by an interfacial dipole \(P_i\), with free energy
$$
F[\{P_i\};E]=\sum_{i=1}^{N-1}\alpha_i P_i^2+\sum_{i=1}^{N-2}\Gamma_{i,i+1}P_iP_{i+1}-E\sum_{i=1}^{N-1}P_i.
$$
This construction yields \(N\) distinct net-polarization values in an \(N\)-layer stack and explains the anomalous intermediate polarization states observed in dual-gate 3R-MoS\(_2\) devices [2211.05277]. In trilayers, the reported switching path consists of two sequential interfacial flips, each with a barrier of \(15.4\) meV per unit cell, separated by an intermediate valley of about \(0.2\) meV per cell. Bilayer, tetralayer, and pentalayer barriers remain in the narrow range \(15.0\)–\(16.1\) meV per unit cell, consistent with predominantly local interfacial energetics [2211.05277].

A distinct multilayer realization appears in altermagnetic CuF\(_2\). Bilayer CuF\(_2\) has two ferroelectric minima, FE-I and FE-II, at sliding displacements \(t_x=\pm a/6\), separated by a nonpolar intermediate state at \(t_x=0\), with \(\Delta E_{\rm IM-FE}=3.83\) meV/atom and \(P_z=\pm 1.23\) pC/m [2603.10907]. In quadrilayer CuF\(_2\), four inequivalent ferroelectric states are identified:
\[
(+2,+2): P_z=+4.22~\text{pC/m},\quad
(+1,-1): P_z=+1.45~\text{pC/m},
\]
\[
(-1,+1): P_z=-1.45~\text{pC/m},\quad
(-2,-2): P_z=-4.22~\text{pC/m}.
\]
Neighboring-state switching requires sliding one bilayer by \(a/6\) with a barrier of about \(3.8\) meV/atom [2603.10907]. This suggests a general multilayer design rule: once interfacial dipoles become individually addressable, sliding ferroelectrics can support multilevel nonvolatile states without changing composition.

## 5. Coupling to topology, magnetism, and spin/layer transport

Sliding ferroelectricity has become a practical route to electrically tunable topological phases. In bilayer MBi\(_2\)Te\(_4\), sliding between oppositely polarized states induces reversible band inversion accompanied by topological phase transitions. Without SOC, the AB′\(_1\) bilayers are ordinary insulators; with SOC, a Bi-\(p\)/Te-\(p\) inversion closes and reopens the gap, producing indirect gaps of 31 meV in GeBi\(_2\)Te\(_4\), 36 meV in SnBi\(_2\)Te\(_4\), and 35 meV in PbBi\(_2\)Te\(_4\). The nontrivial phase is confirmed both by the WCC-derived \(Z_2=1\) invariant and by gapless helical edge states spanning the bulk gap [2506.09317]. The corresponding low-energy description,
$$
H(\mathbf{k})=M\sigma_z + A k_x \sigma_x s_z - A k_y \sigma_y,
$$
encodes the central switching effect: sliding changes the interlayer potential and hybridization such that \(M\to -M\), toggling band inversion [2506.09317]. A related lattice-model and first-principles program on bilayer ScI\(_2\) predicts a sliding-induced sequence of second-order TI, spin-hybrid-order TI, quantum spin Hall, quantum anomalous Hall, and trivial phases, with distinct anomalous Nernst signatures [2506.01402].

Magnetic and metallic extensions show that polarization reversal can be coupled directly to spin splitting and transport. In bilayer Fe\(_5\)GeTe\(_2\), the nonpolar AFM AA′ phase converts, by interlayer sliding, into FE AB or BA ferrimagnetic metallic states with \(\pm P_s \approx 0.075\) pC/m, sliding barrier \(77.2\) meV per unit cell, and net magnetization \(\mp 0.19\,\mu_B\) per cell [2508.07947]. Broken in-plane mirror symmetry lifts the Kramers degeneracy, generating \(k\)-dependent spin splitting of \(\sim 50\)–\(100\) meV near the Fermi crossings, while FE switching reverses polarization, spin splitting, and magnetization simultaneously. The same system exhibits a linear magnetoelectric coefficient \(\alpha \approx 7.1\times 10^{-14}\,\mathrm{G\cdot cm^2/V}\), anomalous Hall conductivity \(\sigma_{xy}(E_F)\approx \pm 0.49\,e^2/h\), and anomalous Nernst conductivity \(\alpha_{xy}^{\rm ANE}(100\,\mathrm K)\approx \pm 0.36\,k_B e/h\) [2508.07947].

In altermagnetic CuF\(_2\), the ferroelectric polarization is directly coupled to layer-locked \(d\)-wave spin splitting. The data are summarized by
$$
\Delta E_{\rm spin}(k)\approx \kappa\,P_z\,[\cos k_x-\cos k_y],
$$
so reversing \(P_z\) reverses the sign of the spin splitting and the layer localization of the spin-polarized carriers [2603.10907]. More generally, rhombohedral \(MX_2\) sliding ferroelectrics preserve a large Zeeman-type splitting throughout the entire sliding path; reported \(K\)-valley values include 148 meV in MoS\(_2\), 184 meV in MoSe\(_2\), 417 meV in WS\(_2\), and 453 meV in WSe\(_2\) [2308.15241]. This indicates that the sliding coordinate can reconfigure polarization while leaving strong SOC-driven spin polarization intact.

## 6. Experimental manifestations, optical responses, and device implementations

Direct experimental access to sliding ferroelectricity has been hindered by weak polarization and leakage in many 2D systems, but an important exception is the amphidynamic crystal \((15\text{-Crown-5})\)Cd\(_3\)Cl\(_6\). Its large band gap \(>4\) eV permits direct macroscopic ferroelectric measurements, with \(P\)–\(E\) loops showing \(P_s \simeq 0.30\)–\(0.40\,\mu\text{C}/\text{cm}^2\) and \(E_c\) rising from \(\sim 40\) kV/cm at 293 K to \(\sim 60\) kV/cm at 273 K [2302.03414]. The same material exhibits SHG intensity reaching \(6\times\) quartz at 300 K, pyroelectric polarization \(P_{\rm pyro}\simeq 0.35\,\mu\text{C}/\text{cm}^2\), and switchable \(180^\circ\) PFM domains under \(\pm 220\) V tip bias [2302.03414]. Because the net polarization decomposes into geometric and sliding contributions,
$$
P_{\rm total}=P_{\rm geo}+P_{\rm slide},
$$
it also demonstrates that sliding ferroelectricity can coexist with other polar mechanisms in a single lattice [2302.03414].

Optical and electronic control channels are now diverse. In 3R-bilayer MoS\(_2\), photoexcitation tunes the sliding polarization over a large range at fixed sliding coordinate, with
$$
P_z(S,n_{\rm ph})\simeq P_{\max}(n_{\rm ph})\sin(2\pi S),
$$
\(P_{\max}(0)=0.82\) pC/m, quenching to \(\approx 0.03\) pC/m at \(n_{\rm ph}=0.2\) e/u.c. and recovering to \(\approx 0.4\) pC/m at \(n_{\rm ph}=0.8\) e/u.c.; the switching barrier changes only modestly, from \(\approx 17\) meV/u.c. in the dark to \(21\) meV/u.c. at \(n_{\rm ph}=0.4\) and \(19\) meV/u.c. at \(n_{\rm ph}=0.8\) [2410.00776]. A separate nonlinear-optical line predicts photo-injected spin-plus-orbital magnetization of order \(0.1\)–\(1\,\mu_B\) in ZrI\(_2\), WTe\(_2\), and MoS\(_2\), with the sign reversing when the sliding dipole reverses [2201.04367]. In \(\beta\)-ZrI\(_2\), the linear electro-optic response is unusually electronic rather than ionic, with \(r_{33}^{\eta}\approx 28.7\) pm/V and \(r_{33}^{\eta,\mathrm{ion}}\lesssim 0.5\) pm/V; under biaxial strain, \(r_{33}^{\eta}\) varies nearly linearly from 16.1 to 43.2 pm/V and follows \(r_{33}^{\eta}\,[\mathrm{pm/V}]\approx 65.5-193.6\,E_g[\mathrm{eV}]\) [2510.03738]. In WTe\(_2\), exciton condensation has been proposed to stabilize the sliding-ferroelectric state beyond single-particle DFT, increasing the barrier from \(\simeq 0.3\) meV to \(\simeq 4\)–\(6\) meV and yielding excitonic gaps of \(\simeq 3\)–\(7\) meV [2510.01465].

Device demonstrations establish that the phenomenon is technologically operative rather than merely model-level. Bilayer-MoS\(_2\) sliding ferroelectric memories exhibit a memory window \(>8\) V, conductance ratio above \(10^6\), retention time \(>10\) years, and programming endurance \(>10^4\) cycles at room temperature; flexible devices retain their performance after \(10^3\) bending cycles and support 32 conductance levels, Hebbian plasticity, and 97.81% MNIST accuracy when mapped into a five-layer CNN [2401.16150]. The sliding ferroelectric resonant tunnel junction extends the concept to two-terminal transport, combining atomically thin polar barriers, superlubric shear-soliton switching, and momentum-conserving resonant tunneling to realize writing voltages below \(0.5\) V, reading biases under \(0.1\) V, current densities exceeding \(50\) nA \(\mu\)m\(^{-2}\), a room-temperature ON/OFF ratio \(>7\), and switching energies below \(1\) fJ within the SLAP architecture [2603.00817].

Sliding ferroelectrics therefore occupy a distinct place within ferroic condensed matter. Their order parameter is set by interlayer registry, their switching can be domain-wall-mediated and superlubric rather than nucleation-limited, and their polarization readily hybridizes with topological, magnetic, optical, and resonant-transport phenomena. A plausible implication is that the most consequential future developments will come from co-designing these couplings—rather than optimizing polarization alone—because the literature already shows nonvolatile electrical control of \(Z_2\) topology, anomalous Hall transport, altermagnetic spin splitting, electro-optic response, and multistate memory in the same sliding-based framework [2506.09317] [2508.07947] [2603.10907] [2510.03738].

Source: https://www.emergentmind.com/topics/sliding-ferroelectrics