Massive Phason in Condensed Matter

• Massive phason is defined as a collective mode in incommensurate structures that acquires a finite gap or damping, deviating from the conventional gapless acoustic behavior.
• Mechanisms such as unscreened long-range Coulomb interactions, impurity pinning, and overdamped diffusion underpin its altered optical, transport, and thermodynamic properties.
• Experimental observations in CDW insulators and moiré superlattices confirm massive phasons through finite-frequency measurements, anomalous damping, and diffusive dynamics.

A massive phason is a phason whose long-wavelength response is no longer that of a gapless acoustic Goldstone mode. Across the literature, the term is used in several related but not identical senses. In incommensurate charge density waves (CDWs), it denotes the phase collective mode of the order parameter after it acquires a finite optical frequency through long-range Coulomb interaction or a pinning gap through disorder. In quasicrystals and other incommensurate structures, the same phrase may refer more broadly to a phason with finite pinning, wavevector-independent damping, overdamped diffusion, or another departure from symmetry-protected acoustic propagation. In moiré systems, “massive phason” is often used more loosely for a soft, low-energy, damped relative-layer sliding mode of the moiré superlattice rather than for a strictly gapped mode at $\mathbf q=0$ (Kim et al., 2022, Jiang et al., 2023, Baggioli et al., 2020, Birkbeck et al., 2024).

1. Conceptual basis and terminology

A phason is a collective degree of freedom associated with phase or internal-coordinate shifts of an incommensurate or quasiperiodic structure. In the standard CDW notation,

$$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$

the amplitudon is the fluctuation of $|\Psi|$ and the phason is the fluctuation of $\phi$ (Kim et al., 2022). In an ideal incommensurate CDW, the phason is expected to be a Goldstone mode with zero gap at $\mathbf q=0$, whereas in quasicrystals the analogous phason variable describes internal rearrangements or shifts in superspace rather than ordinary rigid motion in physical space (Baggioli et al., 2020).

The phrase “massive phason” therefore depends on what destroys or modifies the massless limit. The mechanisms explicitly represented in the literature considered here include unscreened long-range Coulomb interaction, impurity pinning, wavevector-independent damping, and explicit reconstruction or adhesion effects in moiré structures (Kim et al., 2022, Fukuyama et al., 2020, Jiang et al., 2023, Boschi et al., 3 Nov 2025).

Context	Phason meaning	“Massive” meaning
Incommensurate CDW insulator	Phase mode of the CDW order parameter	Finite optical frequency at $\mathbf q=0$ from long-range Coulomb interaction
Pinned Peierls/CDW	Sliding phase mode of the condensate	Finite pinning scale or gap from disorder
Incommensurate or quasiperiodic structure	Collective internal shift mode	Overdamped, diffusive, pinned, or non-acoustic response
Moiré superlattice	Relative-layer sliding mode of the moiré pattern	Usually a soft damped low-energy collective mode, not necessarily a true $\mathbf q=0$ gap

This diversity of usage is central to the subject. A strict encyclopedic treatment must therefore distinguish between a genuine optical gap, a pinning gap, and a relaxational or overdamped scale.

2. Lee–Fukuyama massive phason in charge-density-wave insulators

The most literal usage of the term occurs in the Lee–Fukuyama scenario for incommensurate CDWs. The conventional expectation is that the longitudinal phason remains gapless or nearly so, perhaps with only a small pinning frequency from disorder. Lee and Fukuyama instead argued that if long-range Coulomb interactions are poorly screened, the longitudinal phason is pushed up to a finite frequency at $\mathbf q=0$, producing a fully gapped spectrum and an optical, IR-active phason (Kim et al., 2022).

In the LRA framework, the characteristic frequency is

$$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$

For $(\mathrm{TaSe}_4)_2\mathrm{I}$, the values $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$0 and $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$1 imply

$$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$2

which closely matches the observed mode at $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$3 (Kim et al., 2022).

The material $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$4 is a quasi-one-dimensional CDW insulator with $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$5, a single-particle CDW gap $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$6, and unusually high low-temperature resistivity, implying poor screening of long-range Coulomb interactions (Kim et al., 2022). Time-domain THz emission spectroscopy revealed coherent, narrow-band THz radiation with a sharp line at $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$7, horizontally polarized along the quasi-1D chain direction. The coherent oscillation persisted for more than $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$8, was strongest at $$\Psi(\mathbf{r},t) = |\Psi| e^{i\phi(\mathbf{r},t)} ,$$9, decreased with increasing temperature, and dropped sharply near $|\Psi|$0, approximately $|\Psi|$1 (Kim et al., 2022).

The identification of this mode as a massive phason rests on several convergent observations: the frequency agrees with the Lee–Fukuyama estimate, the polarization follows the CDW chain direction, the intensity is strongest where Coulomb screening is weakest, and the temperature dependence matches prior RPA expectations for a massive phason (Kim et al., 2022). This established the first direct experimental evidence for the massive “Lee-Fukuyama” phason.

A later multimodal ultrafast study on the same material used time-resolved tunneling, time-resolved point-contact, and optical pump-probe reflectance to detect a principal mode at $|\Psi|$2 with the temperature dependence expected for the massive phason, and also a second mode at $|\Psi|$3 with a robust $|\Psi|$4 frequency ratio (Bae et al., 15 Jul 2025). That work interpreted the lower-frequency mode as a daughter phason generated by parametric splitting of the $|\Psi|$5 massive phason into two massless phasons, and framed the mass generation through a Higgs-mechanism analogy tied to unscreened long-range Coulomb interaction (Bae et al., 15 Jul 2025).

3. Pinning, disorder, and transport interpretations in Peierls systems

A second strict meaning of massive phason arises from impurity pinning in a Peierls or CDW condensate. In the clean limit, the phason is an acoustic sliding mode with propagator

$|\Psi|$6

where

$|\Psi|$7

This is the gapless regime associated with ideal sliding conductivity (Fukuyama et al., 2020).

Impurity coupling explicitly breaks the continuous sliding symmetry and modifies the propagator to

$|\Psi|$8

Here $|\Psi|$9 is the pinning-induced restoring term and $\phi$0 is the damping term (Fukuyama et al., 2020). At $\phi$1, the finite constant $\phi$2 acts as the analog of a mass term: the phason is no longer gapless, the zero-frequency divergence in the conductivity is cut off, and the system becomes dielectric or insulating at $\phi$3 (Fukuyama et al., 2020).

This framework also produces a transport meaning of massive phason that differs from the Lee–Fukuyama optical case. In the clean limit, the conductivity takes the Fröhlich form

$\phi$4

while in the pinned regime the static conductivity vanishes at zero temperature (Fukuyama et al., 2020). At finite temperature, thermally activated phason motion still contributes to transport. Within Kubo and Luttinger linear response, the paper derived a phason-drag contribution to thermoelectricity and argued that the Seebeck coefficient can acquire a large magnitude scaling roughly like the square root of resistivity in the regime $\phi$5, with sign opposite to the electronic contribution in the absence of a Peierls gap (Fukuyama et al., 2020). In this transport context, “massive phason” means a pinned phase mode that no longer supports free dc sliding but remains a low-energy collective transport channel.

4. Overdamped, diffusive, and effectively massive phasons

In a broader incommensurate-structure literature, a phason is often called massive when it is no longer a symmetry-protected acoustic mode, even if the dominant effect is damping rather than a literal gap. A generic starting point is

$\phi$6

with the important feature that neither $\phi$7 nor $\phi$8 is required to vanish at $\phi$9 for a phason. Defining $\mathbf q=0$0 and $\mathbf q=0$1, one reaches the overdamped model

$\mathbf q=0$2

in the limit $\mathbf q=0$3 (Jiang et al., 2023).

This description yields a crossover from overdamped diffusion at small $\mathbf q=0$4,

$\mathbf q=0$5

to underdamped propagation at larger $\mathbf q=0$6, with crossover scales

$\mathbf q=0$7

The corresponding spectral function produces three low-frequency density-of-states regimes: Debye-like $\mathbf q=0$8, anomalous $\mathbf q=0$9, and a strong-damping plateau $\mathbf q=0$0 (Jiang et al., 2023). The heat capacity then crosses from $\mathbf q=0$1 to $\mathbf q=0$2 and finally to linear $\mathbf q=0$3, the last described as a glass-like regime (Jiang et al., 2023). A hypothetical Eliashberg analysis further found that phason damping enhances the effective coupling $\mathbf q=0$4 but produces a non-monotonic $\mathbf q=0$5 with a maximum near the underdamped-to-overdamped crossover (Jiang et al., 2023).

A closely related finite-temperature effective field theory for quasicrystals, formulated with Schwinger–Keldysh methods, derived the same telegraph-like structure,

$\mathbf q=0$6

and emphasized that the phason is diffusive at long wavelengths,

$\mathbf q=0$7

because internal phason shifts are symmetries with no associated physical Noether current (Baggioli et al., 2020). Under explicit pinning, the dispersion becomes

$\mathbf q=0$8

and the relaxation rate satisfies

$\mathbf q=0$9

Here the “massive” character is explicitly not a conventional relativistic mass, but rather the presence of a relaxation or pinning scale that removes ordinary sound-like behavior (Baggioli et al., 2020).

Continuum mechanics of quasicrystals gives yet another variant. There the phason is an internal field $\mathbf q=0$0 describing local atomic flips or phase rearrangements, distinct from the displacement field. A diffusive linear model yields a well-posed phason-relaxation dynamics, while a nonlinear gyroscopic term

$\mathbf q=0$1

introduces rotational inertia-like coupling between macroscopic vorticity and phason rate (Bisconti et al., 2015). This suggests a further interpretation of “massive phason” as inertial phason behavior beyond simple diffusion, although the paper does not introduce a separate standard mass term for $\mathbf q=0$2 (Bisconti et al., 2015).

5. Moiré superlattices and the loose usage of “massive phason”

In moiré materials, phasons are the lowest-energy collective modes associated with relative sliding of one layer with respect to the other. Their physical meaning is that of a global or slowly varying shift of the moiré pattern generated by relative layer displacement. Because the moiré pattern acts as an amplifier, tiny atomic displacements can produce large moiré-scale motion (Maity et al., 2023, Birkbeck et al., 2024).

In twisted bilayer graphene, a cryogenic quantum twisting microscope measured a low-energy mode whose coupling to electrons increases as the twist angle decreases. The paper identified this with a layer-antisymmetric moiré phason that modulates the interlayer tunnelling,

$\mathbf q=0$3

and argued that its coupling scales as

$\mathbf q=0$4

in the small-$\mathbf q=0$5 limit for the moiré zone-boundary momentum selected by the experiment (Birkbeck et al., 2024). The same source notes an important subtlety: if momentum is held fixed and $\mathbf q=0$6, the coupling vanishes because a uniform $\mathbf q=0$7 shift merely translates the moiré pattern without changing the electronic energies for incommensurate twist angles (Birkbeck et al., 2024). In this literature, “massive” therefore does not usually mean a strict $\mathbf q=0$8 optical gap.

The transport literature makes this looser usage explicit. For moiré superlattices with interlayer friction, the phason susceptibility is

$\mathbf q=0$9

with $\mathbf q=0$0 (Ochoa et al., 2023). Friction transfers spectral weight from coherent peaks into a broad low-energy diffusive peak, and the resulting electron–phason scattering yields resistivity that remains linear in temperature down to

$\mathbf q=0$1

where $\mathbf q=0$2 (Ochoa et al., 2023). The same mechanism produces a metallic-like linear-in-$\mathbf q=0$3 specific heat of mechanical origin (Ochoa et al., 2023). This is an overdamped-phason picture rather than a strict gapped-mode picture.

Experimentally, minimally twisted bilayer graphene with $\mathbf q=0$4 exhibits a quadratic low-temperature resistivity crossing over to linear-in-$\mathbf q=0$5 above about $\mathbf q=0$6 K, with trends compatible with scattering mediated by longitudinal phasons that dominate over conventional acoustic phonons (Boschi et al., 3 Nov 2025). The phason spectral function used there,

$\mathbf q=0$7

treats the relevant modes as soft but damped, and the paper states that what one may loosely think of as a “massive phason” is not necessarily a true mass gap in the symmetry sense, but a soft mode shifted away from a clean acoustic response by pinning, damping, and reconstruction (Boschi et al., 3 Nov 2025).

A related moiré bilayer study on twisted MoSe$\mathbf q=0$8/WSe$\mathbf q=0$9 found ultra-soft phason modes with energy about $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$0 at $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$1, an almost rigid moiré translation by several nanometers at $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$2 K, and carrier “surfing” speeds $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$3 for $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$4 and $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$5 for $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$6 (Maity et al., 2023). That paper explicitly states that it does not derive a standalone mass formula; its relevance to “massive phason” lies instead in showing that phasons can be dynamic, transport-active collective modes with finite velocity and finite-temperature coherence (Maity et al., 2023).

6. Non-examples, adjacent notions, and scope of the term

Several important phason phenomena are not examples of a massive phason in the strict sense. In icosahedral quasicrystal nucleation, phasons appear as hidden degrees of freedom in the high-dimensional projection framework: uniform phason shifts change the real-space symmetry of an IQC without changing its diffraction pattern or bulk free energy, and they select distinct critical nuclei and minimum-energy paths. That work explicitly states that it does not discuss a “massive phason” in the sense of a gapped phason mode or an effective mass term (Cui et al., 15 Feb 2026). This is a structural and kinetic use of phason shifts, not a dynamical massive-mode theory.

In modulated martensitic Ni$$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$7MnGa, synchrotron x-ray powder diffraction revealed higher-order satellite reflections and phason broadening of satellite peaks, modeled through a fourth-rank covariant strain tensor in the superspace formalism (Singh et al., 2014). The observed broadening was tied to fluctuations in the phase of the incommensurate modulation and favored a soft-mode phonon origin over the adaptive phase model (Singh et al., 2014). This is evidence for phason disorder or phason strain, not for a massive phason in the Lee–Fukuyama or pinned-CDW sense.

Phason engineering in metamaterials likewise uses “phason” differently. In a reconfigurable 2D acoustic crystal, the phason is a tunable coordinate $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$8 controlling a quasiperiodic pattern-generation orbit and giving access to effective 4D quantum Hall physics with $$\omega_{\mathrm{LF}}^2 = \frac{3}{2}\lambda \omega_Q^2 .$$9 and phason-resolved boundary Weyl singularities (Cheng et al., 2020). Here the phason functions as an engineered topological parameter rather than as a massive collective excitation.

Cold-atom phasonic spectroscopy provides another adjacent meaning. In a 1D bichromatic optical lattice, moving only the secondary lattice drives the phasonic coordinate dynamically and produces a nonperturbative high-harmonic plateau with multiphoton resonances up to the $(\mathrm{TaSe}_4)_2\mathrm{I}$0th order (Rajagopal et al., 2019). The experiment probes phasonic dynamics directly, but it does not measure a phason mass in the field-theoretic sense (Rajagopal et al., 2019).

The literature therefore uses “massive phason” with different degrees of strictness. In its narrowest and most specific form, it denotes a CDW phase mode that acquires a finite optical or pinning frequency at $(\mathrm{TaSe}_4)_2\mathrm{I}$1. In a broader incommensurate and moiré literature, it may denote any phason whose low-energy behavior is no longer that of an ordinary massless acoustic mode because of pinning, dissipation, friction, reconstruction, or explicit symmetry breaking. The distinction between these usages is essential for interpreting both experiments and theory.