MEL Ginzburg–Landau Framework for Superconductivity

• The MEL Ginzburg–Landau framework is a theoretical model that integrates momentum‐dependent charge modulations and superconducting order parameter coupling to predict superconductivity in metals.
• It establishes a unified phase classification of metals into MEL-enhanced superconductors, conventional BCS superconductors, and non-superconducting metals using quantitative instability criteria.
• The framework extends conventional GL theory by incorporating nonlocal gradient expansions and higher-order coupling terms to enhance predictions for inhomogeneous phases and interface phenomena.

The Modulated Electron Lattice (MEL) Ginzburg–Landau (GL) framework provides a unified theoretical approach to understanding the emergence and suppression of superconductivity in metals, extending and generalizing the conventional Bardeen–Cooper–Schrieffer (BCS) formalism by incorporating momentum-dependent charge modulations and their coupling to the superconducting order parameter. The framework resolves longstanding selection problems in metallic superconductivity through a predictive criterion based on the spectral properties of the MEL sector. This article systematically develops the core concepts, mathematical structure, classification scheme, and experimental implications of the MEL–GL formalism.

1. Mathematical Structure of the MEL–GL Functional

The MEL–GL theory employs two coarse-grained fields: $\psi(\mathbf{r})$, a complex superconducting order parameter, and $\rho_{\mathrm{MEL}}(\mathbf{r})$, a real scalar field representing slow charge modulations near reciprocal lattice vectors. The coupled free energy functional is

$$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$

where coefficients $\alpha_s(T)$, $\beta_s$, $K_s$ correspond to standard Ginzburg–Landau parameters for superconductivity, and $\alpha(-i\nabla)$ is a Fourier-space operator encapsulating the momentum-dependent stiffness of the MEL mode. $\beta_\rho > 0$ stabilizes the MEL amplitude. The coupling terms ($\gamma_1$, $\gamma_2$) encode symmetry-allowed interactions, with $\rho_{\mathrm{MEL}}(\mathbf{r})$0 relevant for finite-$\rho_{\mathrm{MEL}}(\mathbf{r})$1 MEL modes and $\rho_{\mathrm{MEL}}(\mathbf{r})$2 for homogeneous backgrounds.

2. Momentum-Dependent MEL Stiffness and Instability Criterion

The distinguishing feature of MEL–GL is the explicit momentum dependence of the MEL stiffness kernel,

$\rho_{\mathrm{MEL}}(\mathbf{r})$3

with $\rho_{\mathrm{MEL}}(\mathbf{r})$4 the bare term, $\rho_{\mathrm{MEL}}(\mathbf{r})$5 the analytic part, and renormalizations from electronic charge susceptibility $\rho_{\mathrm{MEL}}(\mathbf{r})$6 and phonon propagator $\rho_{\mathrm{MEL}}(\mathbf{r})$7, weighted by respective couplings $\rho_{\mathrm{MEL}}(\mathbf{r})$8 and $\rho_{\mathrm{MEL}}(\mathbf{r})$9. Softening in $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$0—signaled by a minimum $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$1 at wavevector $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$2—induces a condensation of the MEL field, either at finite $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$3 (modulated phase) or at $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$4 (homogeneous case corresponding to conventional BCS).

The transition to superconductivity in a MEL-enhanced metal requires the existence of a “MEL enhancement window,” characterized by a soft MEL fluctuation mode ($$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$5) and sufficient coupling to $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$6. The effective SC quadratic coefficient is renormalized as

$$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$7

where large static or fluctuation expectation values in $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$8 can drive $$F[\psi,\rho_{\rm MEL}] = \int d^3r\,\left\{ \underbrace{\alpha_s(T)\,|\psi|^2 + \tfrac{\beta_s}{2}\,|\psi|^4 + K_s\,|\nabla\psi|^2}_{\text{Superconducting sector}} + \underbrace{\tfrac12\,\rho_{\rm MEL}\,\alpha(-i\nabla)\,\rho_{\rm MEL} + \tfrac{\beta_\rho}{4}\,\rho_{\rm MEL}^4}_{\text{MEL sector}} + \underbrace{\gamma_1\,\rho_{\rm MEL}\,|\psi|^2 + \gamma_2\,\rho_{\rm MEL}^2\,|\psi|^2}_{\text{Coupling sector}} \right\}$$9 negative, inducing superconductivity in a system that would otherwise remain normal.

3. Classification of Metallic Elements

The MEL–GL formalism naturally partitions metals into three universal classes:

Class	Condition and Mechanism	Examples
I: MEL-enhanced SC	$\alpha_s(T)$0, $\alpha_s(T)$1; superconductivity from large $\alpha_s(T)$2	CDW-prone metals (transition-metal dichalcogenides), some elements under pressure
II: BCS superconductors	$\alpha_s(T)$3, $\alpha_s(T)$4; uniform MEL background stabilizes standard pairing	Al, Sn, Pb, Nb
III: Non-superconducting metals	$\alpha_s(T)$5 for all $\alpha_s(T)$6; MEL amplitude suppressed, no enhancement	Noble metals Cu, Ag, Au

In Class I, the mechanism is enhancement through finite-$\alpha_s(T)$7 MEL fluctuations, often correlated with charge-density-wave tendencies. Class II recovers standard BCS theory as a homogeneous limit within MEL–GL, while Class III metals remain strictly normal due to the absence of soft MEL modes. First-principles electron–phonon coupling parameter $\alpha_s(T)$8 (Allen–Dynes) is correlated: $\alpha_s(T)$9–$\beta_s$0 for Class II, $\beta_s$1 for Class III (Kim et al., 20 Jan 2026).

4. Relation to and Extension of Conventional Ginzburg–Landau Theory

When the minimum of $\beta_s$2 is at $\beta_s$3, MEL condensation yields a homogeneous background: $\beta_s$4 and the effective SC free energy reduces to the standard GL form with renormalized quadratic coefficient,

$\beta_s$5

recovering BCS expressions for coherence length $\beta_s$6 and penetration depth $\beta_s$7 as special cases (Kim et al., 20 Jan 2026). The MEL–GL framework, therefore, generalizes GL theory to spatially modulated electron–lattice states and bridges homogeneous and modulated superconducting regimes.

5. Improved Ginzburg–Landau Construction for Inhomogeneous Phases

Advanced MEL–GL functionals integrate nonlocal amplitude-resummed free energy and systematic gradient expansions. The Improved GL (IGL) formalism (Mannarelli, 2018) employs:

• A moving-average (nonlocal) amplitude $\beta_s$8 with smoothing scale $\beta_s$9 to resum all homogeneous terms;
• Gradient expansion $K_s$0 or explicit sum over high-frequency modes to account for short-wavelength inhomogeneities.

This construction allows for controlled variational evaluation of modulated ansätze—stripes, hexagonal lattices, bubbles—and minimization over amplitudes and wavevectors. The recipe applies directly to charge-density waves, Wigner crystals, and other MEL-type orders.

6. Phenomenology and Experimental Implications

In high-$K_s$1 cuprates, MEL–GL describes short-range modulations coupled to $K_s$2-wave condensates (Kim et al., 3 Dec 2025). Preferred wavevector $K_s$3 r.l.u. is set by electronic/phononic susceptibility profiles. The interplay produces an enhancement window in doping, temperature, correlation length, and disorder where the superfluid stiffness is increased (up to $K_s$4), as shown by Monte Carlo simulations.

Key experimentally testable predictions include:

• FT-LDOS peak sharpening at $K_s$5 concurrent with SC phase coherence below $K_s$6,
• Positive gap-amplitude/MEL correlation in scanning tunneling spectroscopy,
• Disorder-driven percolation thresholds for global coherence,
• Vortex pinning energy landscapes modulated by local MEL amplitude.

In conventional systems, the absence of superconductivity in noble metals is traced to strictly positive $K_s$7 for all $K_s$8; thus, calculating $K_s$9 via DFPT and RPA yields a direct material-selective prediction for the possibility of superconductivity (Kim et al., 20 Jan 2026).

7. Numerical Techniques and Implementation

The hybrid Schrödinger–GL approach (Stojewski et al., 2023) utilizes the mathematical correspondence between the two formalisms to seed GL relaxation algorithms with ground-state Schrödinger wavefunctions, incrementally introducing the nonlinearity to converge to self-consistent modulated solutions. The spatial patterns of $\alpha(-i\nabla)$0 and $\alpha(-i\nabla)$1 encode arbitrary MEL geometries, with iterative schemes enabling the study of interfaces, patterned structures, and time-dependent dynamics.

8. Couplings and Stability of Modulated Phases

Extended MEL–GL models incorporate higher-order gradient and biquadratic couplings to additional long-range order parameters (structural, magnetic, or ferroelectric), as formulated analytically in (Morozovska et al., 2022). The general energy density includes local ($\alpha(-i\nabla)$2) and gradient couplings ($\alpha(-i\nabla)$3, $\alpha(-i\nabla)$4, $\alpha(-i\nabla)$5), enabling explicit calculation of modulation period, amplitude, and stability boundaries for intertwined phases. Closed-form criteria for onset of incommensurate states and analytic phase diagrams guide the interpretation of experimental data in charge–ordered quantum materials.

---

The Modulated Electron Lattice Ginzburg–Landau framework thus provides a rigorous, material-selective, and predictive theoretical foundation for spatially modulated superconductivity, resolving both microscopic selection problems and guiding experimental investigations of intertwined electronic orders (Kim et al., 20 Jan 2026, Mannarelli, 2018, Kim et al., 3 Dec 2025, Stojewski et al., 2023, Morozovska et al., 2022).