---
title: Site-Nematic Insulator in Correlated Materials
url: https://www.emergentmind.com/topics/site-nematic-insulator
type: topic
---

# Site-Nematic Insulator in Correlated Materials

A site-nematic insulator is an interaction-driven insulating state in which rotational symmetry is spontaneously broken through inequivalent sites, sublattices, or closely related intra-unit-cell degrees of freedom. In the most explicit usage, developed for spinless fermions on the checkerboard lattice at half filling, it denotes a strong-coupling, time-reversal-preserving, \(C_4\)-breaking insulator with a charge-density imbalance between the two sublattices [2507.15668]. In other materials classes, closely related language is used for insulating states with local site inequivalence tied to magneto-structural coupling, or for gapped moiré phases with broken \(C_3\) or \(C_{3z}\) symmetry and anisotropic real-space charge patterns [2010.00968; 2008.05041; 2101.04711]. This suggests that the term has both a narrow model-specific meaning and a broader descriptive role for insulating phases with locally broken rotational symmetry.

## 1. Conceptual definition and scope

In the checkerboard-lattice literature, the site-nematic insulator is defined as a charge-ordered insulating state with unequal occupation of the two sublattices \(A\) and \(B\) inside each unit cell, but without the additional spatial modulation characteristic of a stripe phase [1805.01116]. The associated order parameter is the sublattice imbalance,
\[
\frac{\langle n_{i,A}-n_{i,B}\rangle}{2},
\]
which is finite in the site-nematic insulator and vanishes in the symmetric semimetal [1805.01116]. In the newer constrained-path quantum Monte Carlo formulation, the same phase is diagnosed through the staggered sublattice density correlation
\[
\delta_{\mathrm{SNI}} = \frac{1}{N_s} \sum_{\mathbf r} \left\langle \left(c^\dagger_{\mathbf r,1}c_{\mathbf r,1}-c^\dagger_{\mathbf r,2}c_{\mathbf r,2}\right) \left(c^\dagger_{\mathbf r_0,1}c_{\mathbf r_0,1}-c^\dagger_{\mathbf r_0,2}c_{\mathbf r_0,2}\right) \right\rangle,
\]
whose finite long-range value signals site-nematic order [2507.15668].

In insulating transition-metal compounds, the term is used more locally. In La\(_2\)O\(_3\)Fe\(_2\)Se\(_2\), the relevant phenomenon is the emergence of two inequivalent La sites below \(T_{\rm N}\simeq 90\) K, with different quadrupole couplings and likely domain segregation, together with low-frequency fluctuations in the paramagnetic phase that are interpreted as nematic-like [2010.00968]. In the cuprate context, the closely related order is often intra-unit-cell nematicity on oxygen sites,
\[
\eta \equiv \frac{n_x-n_y}{n_x+n_y},
\]
rather than a Cu-site charge imbalance, and it appears as part of a hierarchy including spin stripes, charge stripes, and isotropic phases [1406.2711].

The expression should not be extended indiscriminately. In Cu\(_x\)Bi\(_2\)Se\(_3\), the phase above \(T_c\) is a vestigial nematic phase generated by superconducting fluctuations. It has \(\langle \mathbf{\Delta}\rangle = 0\) but \(\langle \hat q_{\alpha\beta}\rangle \neq 0\), so it is non-superconducting but not an insulator in the localized-electron sense [1712.07523].

## 2. Checkerboard-lattice realization

A canonical realization is the half-filled spinless-fermion checkerboard lattice. One formulation studies
\[
H = -\sum_{ij} (t_{ij} c^\dagger_i c_j + h.c.) + V_1 \sum_{\langle ij\rangle} n_i n_j + V_2 \sum_{\langle\langle ij\rangle\rangle} n_i n_j,
\]
with \(t=1\), \(t' = 0.5\), and \(t''=-0.5\), so that the noninteracting system has a quadratic band touching at half filling [1805.01116]. A related formulation fixes
\[
t_1=1,\qquad t_2=0.5,
\]
and uses nearest-neighbor repulsion \(V\) as the control parameter; in momentum space, the noninteracting model has a quadratic band touching at \(\mathbf k=(\pi,\pi)\), protected by time-reversal symmetry and a fourfold lattice rotational symmetry \(C_4\) [2507.15668].

At strong coupling, the interactions favor localized charge configurations that spontaneously break \(C_4\). In the two-interaction classical limit, the competing insulating energies are
\[
E_{\rm nematic}=V_2,\qquad E_{\rm stripe}=V_1/2,\qquad E_{\rm phase\text{-}separation}=V_1+V_2,
\]
so the site-nematic insulator is one of the strong-coupling endpoints of the phase diagram [1805.01116]. In the constrained-path quantum Monte Carlo study with nearest-neighbor repulsion only, the half-filled system displays three interaction-induced phases: a quantum anomalous Hall insulator at weak coupling, a bond-nematic Dirac semimetal at intermediate coupling, and a site-nematic insulator at strong coupling [2507.15668].

The symmetry content is sharp. The site-nematic insulator breaks \(C_4\) lattice rotation symmetry, preserves time-reversal symmetry, and exhibits an imbalance of electron density between sublattices 1 and 2 [2507.15668]. It differs from the bond-nematic Dirac semimetal because both phases break the same \(C_4\) symmetry, but the bond-nematic state remains gapless with two Dirac cones, whereas in the site-nematic insulator the Dirac cones annihilate and a full gap opens [2507.15668]. It also differs from the stripe insulator, which has additional translation-symmetry breaking; in the stripe state the neighboring-site difference
\[
\frac{\langle n_{i,A}-n_{i+\hat{x},A}\rangle}{2}
\]
is finite, while it is zero in the site-nematic insulator [1805.01116].

The numerical evidence is correspondingly specific. On finite torus geometries, the bond-nematic Dirac semimetal to site-nematic insulator transition occurs near \(V \approx 1.56\text{ to }1.7\), depending on size and boundary conditions, and is identified as first order [2507.15668]. Finite-size scaling of the site-nematic order parameter gives
\[
\delta_{\mathrm{SNI}}(1/L=0)=0.1421(2)
\]
at \(V=1.80\), supporting long-range order in the thermodynamic limit [2507.15668]. At the same time, weak-coupling analytical work rejects a weak-coupling site-nematic instability in this model: the nematic semimetal proposed for \(|V_2| \gg V_1\) at weak coupling is absent, and the quantum anomalous Hall state is the only weak-coupling instability of the spinless quadratic-band-touching semimetal [1805.01116].

## 3. Mott and charge-transfer settings

A different, experimentally anchored realization appears in the Mott insulator La\(_2\)O\(_3\)Fe\(_2\)Se\(_2\). \(^{139}\)La NQR shows a single La site in the paramagnetic phase and two inequivalent La sites, La1 and La2, in the antiferromagnetically ordered phase below \(T_{\rm N}\simeq 90\) K [2010.00968]. The two sites have different quadrupole couplings,
\[
\nu_Q^{(1)} \neq \nu_Q^{(2)},
\]
implying different local electric-field gradients and therefore distinct lattice and/or charge configurations [2010.00968]. Their comparable intensity supports an intrinsic bulk origin rather than impurity physics, and the authors suggest distorted and undistorted domains, or possibly a charge-disproportionated phase-separated state [2010.00968]. Within the paper’s interpretation, this is the microscopic signature of a site-nematic Mott insulator: local \(C_4\) symmetry is broken and atomic sites become inequivalent without requiring a large global orthorhombic distortion [2010.00968].

The temperature dependence links this site inequivalence to magnetism. The quadrupole frequency of La2 follows the magnetic order-parameter-like evolution below \(T_{\rm N}\), which the authors interpret as evidence that the magnetic order parameter drives the structural distortion [2010.00968]. Dynamically, the transverse relaxation rate \(1/T_2\) begins increasing around \(150\) K, has a clear maximum around \(108\) K, diverges at \(T_{\rm N}\), and reflects fluctuations with characteristic times in the tens of microseconds, corresponding to MHz-scale frequencies far below the Heisenberg exchange scale [2010.00968]. In the fast-fluctuation regime,
\[
\frac{1}{T_2} \simeq \gamma^2 \langle \Delta h_Z^2 \rangle \tau_c,
\qquad
\tau_c \propto e^{E_A/T},
\qquad
E_A \approx 990 \pm 90\ \text{K},
\]
which the paper identifies as strongly reminiscent of nematic fluctuations in frustrated spin systems [2010.00968].

The cuprate literature develops the same theme more microscopically, though not always as a standalone insulating phase. In a three-band Emery model, quadrupolar fluctuations in the oxygen \(p\)-orbitals inevitably generate a biquadratic coupling between Cu spins, and this coupling favors local stripe-like magnetic fluctuations that enhance the \(B_{1g}\) nematic susceptibility [1703.02210]. The resulting susceptibility,
\[
\chi_{\rm nem}(T) = \int_0^{1/T}d\tau\sum_i \left\langle \mathcal{T}_\tau\,\varphi_i(\tau)\varphi_0(0)\right\rangle,
\]
is amplified by the oxygen-fluctuation-induced biquadratic exchange, but its peak remains controlled by the magnetic scale \(J\), not by the biquadratic coupling itself [1703.02210]. This mechanism explains enhanced nematic fluctuations near an antiferromagnetic Mott insulator, rather than a spontaneous long-range site-nematic phase by itself [1703.02210].

A related strong-coupling three-band model for hole-doped charge-transfer insulators finds that doped oxygen holes frustrate the Cu antiferromagnetic background through Kondo-type Cu–O exchange, stabilizing a hierarchy of spin stripes, charge stripes, intra-unit-cell nematic order, and isotropic phases [1406.2711]. The order parameter
\[
\eta \equiv \frac{n_x-n_y}{n_x+n_y}
\]
measures oxygen-site inequivalence within a CuO\(_2\) unit cell [1406.2711]. The work supports a robust intra-unit-cell nematic sector tied to charge-transfer insulating physics, while also emphasizing that the clearest explicitly insulating states are stripe insulators rather than an isolated pure nematic insulator across generic parameters [1406.2711].

## 4. Moiré graphene nematic insulators

In moiré graphene systems, “nematic insulator” refers to gapped correlated states that spontaneously break the rotational symmetry of the moiré lattice. Near charge neutrality in twisted bilayer graphene, a generalized Bistritzer–MacDonald continuum model with corrugation effects and symmetry-classified contact interactions yields two classes of quartic interactions: one favoring gap opening and one favoring nematic ordering [2008.05041]. The combined group-theory and renormalization-group analysis concludes that proximity to the first magic angle favors a layer-polarized, gapped state with a spatial modulation of interlayer correlations, termed a nematic insulator [2008.05041]. In the paper’s language, the state is insulating because a charge gap opens, and nematic because \(C_3\) is broken by stripe-like modulation of interlayer density or interlayer correlations [2008.05041].

At half filling of the flat bands in twisted bilayer-monolayer graphene and twisted double-bilayer graphene under finite displacement fields, unrestricted Hartree–Fock calculations find spin-polarized, \(C_{3z}\)-broken insulator states with zero total Chern number [2101.04711]. The dominant order parameters are of the form
\[
\tau _{0,z}s_{0,z}\sigma _k \quad (k=x,y),
\]
and the resulting real-space charge density becomes stripe-like and nematic [2101.04711]. If only the dominant intra-valley Coulomb interaction is included, the spin-polarized insulator is quasi-degenerate with a valley-polarized state, but adding atomic on-site Hubbard interactions lowers the energy of the spin-polarized state by about \(0.4\!\sim\!0.6\) meV per electron, making it the unique ground state [2101.04711]. In some displacement-field regimes these nematic insulators are also quantum valley Hall states, with zero total Chern number but nontrivial valley Chern numbers [2101.04711].

A nearby semimetallic counterpart also exists. Self-consistent Hartree–Fock calculations for magic-angle twisted bilayer graphene at charge neutrality, restricted to valley-\(U(1)\)-preserving orders, find three classes of competitive states: \(C_2\mathcal T\)-breaking insulators, spin- or valley-polarized insulators, and \(C_3\)-breaking semimetals whose gaplessness is protected by the topology of the moiré flat bands [1905.07409]. Weak strains that break \(C_3\) stabilize both the nematic semimetal and, somewhat unexpectedly, the \(C_2\mathcal T\)-breaking insulators [1905.07409]. This places the moiré nematic insulator in a landscape where gapped and gapless rotational-symmetry-broken states can be very close in energy.

## 5. Diagnostics and empirical signatures

The identification of site-nematic insulating behavior depends strongly on the microscopic platform, but the core diagnostic is always a static or fluctuating anisotropy tied to site, sublattice, or local orbital inequivalence.

| Platform | Primary diagnostic | Interpretation |
|---|---|---|
| Checkerboard lattice | Finite \(\delta_{\mathrm{SNI}}\); full gap | Sublattice density imbalance with broken \(C_4\) |
| La\(_2\)O\(_3\)Fe\(_2\)Se\(_2\) | One La site above \(T_{\rm N}\), two below; slow \(1/T_2\) fluctuations | Local site inequivalence and magnetically induced distortion |
| Twisted graphene | Broken \(C_3\) or \(C_{3z}\); stripe-like/local filling anisotropy | Gapped nematic moiré state |

In the checkerboard model, the order parameter itself is a density-density correlation of sublattice imbalance, and finite-size scaling supplies a direct thermodynamic test. The constrained-path quantum Monte Carlo results show long-range site-nematic order in the strong-coupling regime, while the simultaneous disappearance of bond-nematic order tracks the transition from the bond-nematic Dirac semimetal into the gapped site-nematic state [2507.15668].

In La\(_2\)O\(_3\)Fe\(_2\)Se\(_2\), the evidence is spectroscopic and dynamical. Two inequivalent quadrupole frequencies below \(T_{\rm N}\) imply distinct local electric-field gradients, while the activated slowing down seen in \(1/T_2\) reveals low-frequency fluctuations already in the paramagnetic phase [2010.00968]. The coexistence of local inequivalence and slow fluctuations is central to the paper’s site-nematic interpretation [2010.00968].

In moiré graphene, the observables are spatial and topological. The charge-neutral twisted-bilayer analysis predicts a gapped, layer-polarized state with broken \(C_3\) and spatial modulation of interlayer correlations [2008.05041]. In twisted multilayer graphene, the broken-\(C_{3z}\) insulator produces stripe-like charge density and can carry nontrivial valley topology even when the total Chern number is zero [2101.04711]. In the magic-angle bilayer Hartree–Fock study, the local filling fraction
\[
\nu(\mathbf{r}) = 8\left(\frac{\rho_{\rm LB}(\mathbf{r})}{\rho_{\rm LB}(\mathbf{r})+\rho_{\rm UB}(\mathbf{r})}-\frac12\right)\in[-4,4]
\]
displays clear \(C_3\)-breaking patterns compatible with STM observations near charge neutrality [1905.07409].

## 6. Relation to neighboring phases and common confusions

The site-nematic insulator is often confused with several related but distinct phases. The first distinction is with the stripe insulator. In the checkerboard model, the site-nematic insulator has finite sublattice imbalance but zero stripe order, whereas the stripe insulator has both sublattice imbalance and translation-symmetry breaking along one direction [1805.01116]. The second distinction is with the bond-nematic Dirac semimetal: both break \(C_4\), but the bond-nematic state remains gapless, while the site-nematic state is fully gapped after Dirac-cone annihilation [2507.15668].

A different confusion concerns spin-nematic insulators. In the spin-1 bilinear-biquadratic model on the triangular lattice, the spin-nematic insulator is a ferro-quadrupolar state with no conventional magnetic dipole order,
\[
\langle \mathbf S \rangle = 0
\]
in the zero-field limit, but with quadrupolar order and hidden spin symmetry breaking [2409.08246]. This is not a charge-ordered site-nematic insulator; its order parameter is quadrupolar in spin space, and its experimental probe in the cited work is spin pumping from a ferromagnetic insulator into a spin-nematic insulator [2409.08246].

The vestigial nematic phase proposed for Cu\(_x\)Bi\(_2\)Se\(_3\) is likewise distinct. There the nematic order is composite and fluctuation-induced, derived from a two-component superconducting order parameter,
\[
\mathbf{\Delta}=\big(\Delta^{(x)},\Delta^{(y)}\big)^T,
\]
with \(T_{\rm nem}>T_c\) and no superconducting long-range order above \(T_c\) [1712.07523]. The phase is “insulating” only in the sense that it is not superconducting; it is not an electronically localized insulator [1712.07523].

Finally, nematic order in Dirac systems need not imply an insulating phase. On the surface of a three-dimensional topological insulator, time-reversal-invariant nematic order spontaneously breaks rotational symmetry but distorts rather than generically gaps the Dirac cone [1702.07364]. At the Dirac point, the paper finds a first-order transition between isotropic and nematic Dirac semimetals, while at finite doping it finds a second-order transition between isotropic and nematic helical Fermi liquids [1702.07364]. This suggests that “nematic” and “insulating” should be treated as separate descriptors unless the underlying work explicitly establishes both.

Across the literature, the most restrictive meaning of site-nematic insulator remains the checkerboard-lattice phase with sublattice density imbalance and broken \(C_4\) symmetry. Broader usages extend the idea to Mott, charge-transfer, and moiré systems whenever an insulating state acquires local site inequivalence or symmetry-breaking anisotropy in real space.

Source: https://www.emergentmind.com/topics/site-nematic-insulator