Periodic Lattice Distortion (PLD)

• PLD is a modulation of atomic positions where each nucleus is displaced periodically, defining structural supercells in charge density wave materials.
• Techniques like HAADF-STEM, STM, and reciprocal-space diffraction enable sub-picometer resolution mapping of PLD fields and defect structures.
• Coupled electron-phonon interactions, Fermi surface nesting, and excitonic effects drive PLD, linking electronic and elastic free energies in various systems.

Searching arXiv for recent and foundational papers on periodic lattice distortion to ground the article in published work. Periodic lattice distortion (PLD) denotes a static or quasi-static modulation of atomic positions in which each nucleus is displaced from its high-symmetry site by a periodic field. In CDW materials, PLD is the structural counterpart of the charge modulation: conduction electrons condense at wavevector(s) $\mathbf{q}$, and the lattice follows by displacing its nuclei so as to lower the combined electronic and elastic free energy. A general representation is

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$

or, in a complex three-$q$ form,

$$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$

Commensurate PLDs use $\mathbf{q}_j$ that are simple rational fractions of reciprocal-lattice vectors, whereas nearly commensurate and incommensurate states deviate slightly in magnitude or direction. Although the concept emerged most prominently in CDW studies, the same language also describes spin-Peierls distortions, charge-ordered oxides, and torsional relaxation fields in twisted two-dimensional moiré systems (Hovden et al., 2016, Hildebrand et al., 2017, Sung et al., 2022).

1. Mathematical form and symmetry content

A PLD is a displacement field rather than a density field. In the simplest one-mode description,

$$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$

while for two coexisting modulations one may write

$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$

Here $\mathbf{q}_i$ fixes the periodicity, $A_i$ is the amplitude vector, and $\phi_i(\mathbf{r})$ carries local information about defects, dislocations, and domain walls. In BSCMO, the amplitude vectors are transverse to $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$0; in twisted bilayers, the analogous field is written as $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$1 with $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$2 (Savitzky et al., 2017, Sung et al., 2022).

The symmetry content of a PLD is often more informative than its amplitude alone. In the commensurate $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$3 CDW of $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$4-TiSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$5, the three principal wavevectors are $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$6, $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$7, and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$8 in reciprocal-lattice units, with dominant in-plane polarizations and an antiphase relation between adjacent layers,

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$9

which restores inversion symmetry and forbids net handedness along $q$0 (Hildebrand et al., 2017). In phenomenological form, the coupled electronic and lattice sectors may be expressed through

$q$1

or, alternatively, through a charge-lattice functional

$q$2

both of which encode the lock-in between electronic order and ionic displacement (Payne et al., 2020, Hildebrand et al., 2017).

2. Canonical realizations in layered chalcogenides

Layered chalcogenides provide the standard exemplars of PLD physics because their electronic phase transitions are frequently inseparable from real-space superstructures. In $q$3-TaS$q$4, the low-temperature commensurate phase contains three principal CDW/PLD wavevectors,

$q$5

forming the “Star-of-David” $q$6 supercell. Each vector lies roughly $q$7 from the nearest $q$8–M direction and has magnitude $q$9, corresponding to a $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$0 periodicity. In the nearly commensurate phase, reached for $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$1, each $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$2 rotates back toward $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$3 and contracts by $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$4 to $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$5, breaking the modulation into domains of width $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$6–$$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$7 separated by discommensurations (Hovden et al., 2016).

In $$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$8-TiSe$$\mathbf{u}(\mathbf{r})=\sum_{j=1}^3\left[A_j\hat e_j e^{i\mathbf{q}_j\cdot\mathbf{r}}+\text{c.c.}\right].$$9, the PLD is fully commensurate and triples the periodicity in all three directions through a $\mathbf{q}_j$0 superstructure. Experiment and DFT relaxation following Di Salvo et al. find that the dominant displacements lie in-plane; Se atoms in the upper Se layer displace by $\mathbf{q}_j$1 at $\mathbf{q}_j$2, while direct STM places any vertical displacement below resolution, $\mathbf{q}_j$3, consistent with DFT estimates $\mathbf{q}_j$4. The three branches have equal amplitude, and adjacent TiSe$\mathbf{q}_j$5 sandwiches are locked in antiphase (Hildebrand et al., 2017).

Under compression, SnSe$\mathbf{q}_j$6 develops a pressure-induced commensurate PLD characterized by a $\mathbf{q}_j$7 supercell at

$\mathbf{q}_j$8

Above $\mathbf{q}_j$9, single-crystal XRD reveals new Bragg peaks indexed as $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$0, Raman spectroscopy shows new low-frequency modes at $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$1 and $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$2, and transport measurements display a kink in $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$3 at $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$4 only for $$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$5 (Ying et al., 2018).

System	Modulation vector / supercell	Quantitative signature
$$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$6-TaS$$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$7	$$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$8 “Star-of-David”	$$\mathbf{u}(\mathbf{r})=A\cos(\mathbf{q}\cdot \mathbf{r}+\phi),$$9, $$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$0 period, $$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$1–$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$2
$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$3-TiSe$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$4	$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$5	Se displacement $$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$6, $$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$7, antiphase stacking
SnSe$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$8 under pressure	$$\Delta_i(\mathbf{r})=A_i\sin(\mathbf{q}_i\cdot \mathbf{r}+\phi_i).$$9 at $\mathbf{q}_i$0	onset above $\mathbf{q}_i$1, Raman modes at $\mathbf{q}_i$2 and $\mathbf{q}_i$3

These examples show that PLD may be nearly commensurate, strictly commensurate, or pressure induced, and that the relevant length scale may range from picometer displacements to nanometer supercells (Hovden et al., 2016, Hildebrand et al., 2017, Ying et al., 2018).

3. Microscopic driving mechanisms

The standard starting point is a coupled electron-phonon Hamiltonian,

$\mathbf{q}_i$4

with a static lattice distortion $\mathbf{q}_i$5 below the transition and an electronic order parameter $\mathbf{q}_i$6. This framework permits both the canonical Peierls interpretation and a broader class of strong-coupling scenarios in which lattice, charge, and sometimes excitonic degrees of freedom are co-primary (Payne et al., 2020).

In SnSe$\mathbf{q}_i$7, density functional theory attributes the pressure-induced PLD to the combined effect of strong Fermi surface nesting and electron-phonon coupling at $\mathbf{q}_i$8. The static nesting function $\mathbf{q}_i$9 develops a pronounced peak at precisely that wavevector under compression, while the corresponding DFPT phonon branch softens to $A_i$0 and the mode-resolved EPC constant becomes anomalously large, of order unity. This is explicitly contrasted with canonical CDW-bearing TMDs, where significant Fermi surface nesting is not usually involved (Ying et al., 2018).

In $A_i$1-TiSe$A_i$2, a different microscopic route has been formulated. An exciton-condensate mean-field model with order parameter

$A_i$3

using $A_i$4 and $A_i$5, generates a condensate-induced force on the lattice through the electron-phonon term. Minimization of the phonon and exciton-phonon energies yields a static normal-coordinate shift and, for the Ti transverse mode at the $A_i$6 point, a calculated displacement amplitude $A_i$7 at $A_i$8, compared with the neutron-diffraction value $A_i$9 at $\phi_i(\mathbf{r})$0. The calculation reproduces the order of magnitude of the experimental PLD amplitude without invoking an independent Jahn-Teller instability (Monney et al., 2010).

The monolayer limit of $\phi_i(\mathbf{r})$1-TiSe$\phi_i(\mathbf{r})$2 reinforces the mixed electronic-lattice character. DFT finds a Kohn-type soft mode of $\phi_i(\mathbf{r})$3 symmetry at the $\phi_i(\mathbf{r})$4 point in the undistorted $\phi_i(\mathbf{r})$5 structure, a relaxed $\phi_i(\mathbf{r})$6 PLD with in-plane displacements $\phi_i(\mathbf{r})$7 and $\phi_i(\mathbf{r})$8, and a total-energy lowering

$\phi_i(\mathbf{r})$9

in GGA, or $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$00 in LDA. The relaxed phase opens a direct gap of $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$01 in GGA and shows a “Mexican-hat” valence-band profile. The paper identifies these results as evidence that lattice dynamics and electron-electron interactions conspire to stabilize the CDW/PLD state (Singh et al., 2017).

In $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$02-NbSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$03, time-resolved optical spectroscopy provides direct evidence against a purely electronic ordering picture. At equilibrium, CDW gap opening and PLD formation occur simultaneously at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$04, but under moderate photoexcitation the electronic order is destroyed while the PLD survives for tens of picoseconds. The same material is also discussed in terms of finite-$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$05 phonon softening and strong anharmonic phonon-phonon interactions rather than an adequate low-energy nesting instability (Payne et al., 2020). Independent STM work on $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$06-TaS$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$07, TaSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$08, and NbSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$09 likewise reports dominant lattice contributions rather than the electronic modulation expected from a weak-coupling Peierls transition, with measured PLD amplitudes of $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$10 in $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$11-TaS$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$12, $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$13 in $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$14-NbSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$15, and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$16 in $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$17-TaSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$18 (Dai et al., 2013).

4. Imaging and quantifying PLD

The recent literature treats PLD as a directly measurable vector field rather than an inferred secondary order parameter. In exfoliated $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$19-TaS$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$20, atomic-resolution HAADF-STEM at room temperature and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$21 was performed at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$22 with a semi-angle of $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$23, with plan-view data on the $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$24 zone axis and cross-section data on $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$25. Ta atomic columns were fitted to six-parameter Gaussians and centroided with sub-$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$26 precision; FFTs of both raw images and centroid-extracted positions were then used to isolate genuine displacement peaks. In the commensurate phase at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$27, the resulting Ta displacements show clear $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$28 periodic modulations with amplitudes $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$29–$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$30 (Hovden et al., 2016).

A related STEM strategy was developed for the charge-ordered manganite Bi$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$31Sr$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$32Ca$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$33MnO$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$34. Aberration-corrected HAADF-STEM at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$35 with a $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$36 convergence beam and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$37–$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$38 collection range was combined with rigid and non-rigid registration of $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$39–$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$40 frames and two-dimensional Gaussian fits of atomic columns to achieve $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$41 precision. The PLD field was extracted by damping selected FT satellites, inverse transforming to construct a reference lattice, and subtracting the reference positions from the original coordinates. This yielded peak amplitudes $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$42 and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$43 (Savitzky et al., 2017).

STM accesses PLD through several distinct observables. In $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$44-TiSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$45, intentionally introduced interstitial Ti atoms in the van der Waals gap serve as local markers of the PLD: their asymmetric triangular contrast at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$46, together with Tersoff-Hamann simulations, identifies which Se sites are displaced and which correspond to CDW maxima. At a monoatomic step of height $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$47, defects on the two terraces reveal opposite handedness and directly image the $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$48 phase shift between adjacent sandwiches (Hildebrand et al., 2017). In $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$49-TaS$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$50, TaSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$51, and NbSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$52, constant-current topographs at opposite biases are decomposed as

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$53

thereby separating the symmetric PLD component from the antisymmetric electronic modulation (Dai et al., 2013). In Bi$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$54Sr$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$55CaCu$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$56O$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$57, affine correction and the Lawler-Fujita drift-correction algorithm applied at coarse-graining lengths $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$58 and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$59 yield a bond-length map with $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$60 precision over a $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$61 field of view (Du et al., 2023).

Reciprocal-space diffraction remains central, but in several systems it now serves as a quantitative PLD metrology rather than only a superstructure detector. In twisted bilayer graphene and related moiré systems, a transverse torsional PLD produces superlattice peaks whose $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$62th-order intensity scales as

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$63

allowing the harmonic amplitudes to be extracted directly from selected-area electron diffraction with picometer sensitivity (Sung et al., 2022).

5. Domains, defects, stacking, and non-equilibrium dynamics

PLD fields are rarely uniform over macroscopic length scales. In $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$64-TaS$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$65, cryogenic STEM resolves one–Ta–Ta–bond-length stacking faults in which the trigonal $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$66 registry switches to hexagonal $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$67 over a narrow $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$68 boundary. The lateral shift is $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$69, and cross-sectional geometric phase analysis finds $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$70 local compressive strain at the boundary. Yet FFTs on either side continue to show nearly commensurate PLD peaks, so the electronic CDW persists across the stacking fault. The same study reports that NC PLDs exist inside both the stacking domains and their boundaries (Hovden et al., 2016).

In BSCMO, the combined displacement field $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$71 reveals locally unidirectional stripe domains as small as $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$72, despite an apparently bidirectional Fourier pattern on larger length scales. A $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$73 domain boundary can host an edge-type PLD dislocation in which one wavefront terminates abruptly, the Burgers vector equals one PLD wavelength, the phase winds by $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$74 around the core, and the local amplitude collapses in an inlet extending from the boundary to the defect. Shear strain extracted from the phase field,

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$75

maximizes where $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$76 is minimal and nascent $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$77 appears (Savitzky et al., 2017).

Ultrafast experiments show that PLD can persist even when the associated electronic order is transiently suppressed. In $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$78-NbSe$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$79, a $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$80 probe detects a sign change in the slow differential-reflectivity component $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$81 at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$82 and a divergent relaxation time $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$83 on approaching $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$84 from below. At $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$85, the fluence dependence shows a critical window between $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$86 and $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$87 in which the electronic CDW gap collapses while the PLD remains, yielding a non-equilibrium state with $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$88 but $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$89 for tens of picoseconds (Payne et al., 2020).

In K$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$90MoO$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$91, femtosecond x-ray diffraction tracks the PLD coordinate $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$92 in an effective quartic potential,

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$93

with equation of motion

$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$94

where $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$95. A two-pulse optical scheme with $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$96 produces a second revival of the PLD superlattice intensity at $$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$97 by re-exciting the electronic population and further suppressing the damping. This directly demonstrates optical control over PLD coherence in a photoinduced high-symmetry potential (Neugebauer et al., 2019).

6. Beyond canonical CDWs: cuprates, spin-Peierls systems, moiré materials, and controversies

PLD is not restricted to transition-metal dichalcogenides. In Bi$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$98Sr$$\mathbf{u}(\mathbf{r})=\sum_j A_j\cos(\mathbf{q}_j\cdot \mathbf{r}+\phi_j),$$99CaCu$q$00O$q$01, the bond-length map $q$02 has mean $q$03 and a full width at half maximum $q$04, corresponding to $q$05 disorder. Fourier analysis of $q$06 and $q$07 yields peaks at $q$08, indicating an approximately four-unit-cell modulation. Group-theoretic projection onto the $q$09 representations identifies PLD components in the $q$10 and $q$11 channels, and the same wavevector appears in the d-symmetry CDW measured via the spectroscopic ratio $q$12. The paper interprets the $q$13-sector PLDs as a locally frozen version of the bond-stretching soft phonon at $q$14 (Du et al., 2023).

A spin-Peierls realization is proposed for CeRu$q$15Al$q$16. There, the relevant distortion involves Al(1)–Al(4) atoms moving toward their unique nearest-neighbor Ce atoms, while Ce, Ru, and Al(5) remain immobile. The displacements are written as $q$17 with magnitudes of order $q$18, and the resulting three-dimensional order is characterized by

$q$19

This pattern removes the base-centering translation, is compatible with the appearance of forbidden reflections such as $q$20, and accounts for the two-site splitting of the $q$21Al-NQR lines for Al(1)–Al(4) below $q$22 while leaving Al(5) unsplit (Hanzawa, 2010).

Moiré systems introduce an additional structural category: torsional PLD. In twisted bilayer graphene, the amplitude of the single-harmonic torsional order parameter obeys

$q$23

with $q$24 and $q$25. The torsional PLD begins to form below $q$26, reaches $q$27 around the magic angle of $q$28, and has a geometric upper bound of $q$29. At very low twist angles, additional harmonics arise, the second, third, and fourth becoming $q$30 at approximately $q$31, $q$32, and $q$33, respectively, thereby sharpening Bernal domains and narrowing solitonic boundaries to a minimal width of $q$34. Similar torsional distortions are reported in twisted WS$q$35, CrI$q$36, and WSe$q$37/MoSe$q$38 (Sung et al., 2022).

Several controversies in the PLD literature concern causality and symmetry. In the canonical 2H dichalcogenides, STM separation of symmetric and antisymmetric topographic channels finds that the PLD contribution dominates over the purely electronic modulation, challenging a weak-coupling Peierls interpretation (Dai et al., 2013). In $q$39-NbSe$q$40, nonequilibrium measurements go further by exhibiting a state in which electronic order is melted but the lattice distortion remains (Payne et al., 2020). In $q$41-TiSe$q$42, direct STM imaging with interstitial Ti markers shows inversion-symmetric achiral stacking and explicitly contradicts claims of helical CDW stacking and associated chiral order; FFT spot intensities on adjacent terraces are identical to within $q$43 (Hildebrand et al., 2017). Taken together, these results place PLD not as a passive consequence of electronic order in every material, but as an experimentally resolvable degree of freedom whose symmetry, amplitude, defects, and dynamics can determine the character of the ordered phase itself.