Second-Order Nonlinear Electrical Response

• Second-order NLER is the quadratic response to an applied electric field, characterized by the first hyperpolarizability tensor at the molecular scale and second-order conductivity in bulk media.
• It manifests in phenomena such as second-harmonic generation, nonlinear Hall effects, and electro-optic modulation through mechanisms like symmetry breaking, structural reconfiguration, and gradient-induced responses.
• Recent research highlights its tunability in hyperbolic media, poled dielectrics, and twisted graphene systems, providing practical routes for advanced photonic and electronic applications.

Searching arXiv for the cited papers and related NLER work to ground the article in current literature. Second-order nonlinear electrical response (NLER) is the component of an electrical or optical response that is quadratic in the applied field. At the molecular level it appears through the first hyperpolarizability tensor $\beta_{ijk}$ in $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$, while in macroscopic media it is encoded in a second-order conductivity or susceptibility tensor, as in $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$ or $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$ (Srivastava, 2019). Although the conventional electric-dipole $\chi^{(2)}$ vanishes in centrosymmetric media for spatially uniform fields, the broader NLER landscape includes gradient-induced bulk response in centrosymmetric hyperbolic media, extrinsic skew-scattering and side-jump mechanisms in inversion-broken conductors, quantum-geometric contributions governed by Berry curvature dipoles and quantum metric dipoles, and formal many-body constructions based on Matsubara response theory (Narimanov, 22 May 2026).

1. Definition and constitutive structure

Second-order NLER is the part of the response described by a term quadratic in the driving electric field. In nonlinear optics and response theory, the dipole moment of a molecule in an external field is expanded as

$$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$

with $\beta_{ijk}$ the first hyperpolarizability tensor and therefore the molecular second-order response (Srivastava, 2019). In condensed-matter transport, the analogous expansion is

$$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$

or equivalently $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$, depending on notation (Duan et al., 2022).

This quadratic response underlies second-harmonic generation, sum-frequency generation, difference-frequency generation, the nonlinear Hall effect, bilinear magnetoelectric response, and the Pockels electro-optic effect (Srivastava, 2019). In integrated photonics and nonlinear optics, the relevant macroscopic constitutive law is often written as

$$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$

while in transport the observable is a second-order conductivity tensor entering $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$0 (Narimanov, 22 May 2026).

A recurring symmetry statement is that inversion symmetry forbids the ordinary electric-dipole second-order response for spatially uniform fields. In a centrosymmetric molecule at zero field, $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$1; in a centrosymmetric bulk medium, $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$2 vanishes in the dipole approximation (Srivastava, 2019). The data block, however, also documents several mechanisms by which a measurable second-order response survives or is engineered despite nominal centrosymmetry: externally induced symmetry breaking in molecules, electrical poling in amorphous dielectrics, and gradient-induced response in hyperbolic media (Zhang et al., 2022).

2. Symmetry, inversion breaking, and routes to finite response

In conventional nonlinear optics, inversion symmetry eliminates the uniform-field electric-dipole term because under $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$3, both $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$4 and $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$5 change sign, so a term proportional to $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$6 is incompatible with centrosymmetry. This is the standard reason why a bulk $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$7 is absent in centrosymmetric crystals and why inversion breaking is a prerequisite for ordinary second-order response (Narimanov, 22 May 2026).

One route to finite NLER is direct symmetry breaking by an external field. Hexalithiobenzene, $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$8, is planar and highly symmetric at zero field, with $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$9 D and $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$0 a.u., so its second-order response is effectively negligible. Under an oriented external electric field applied along the molecular $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$1-axis, the symmetry is broken, a finite dipole develops, and the first mean hyperpolarizability rises to $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$2 a.u. at $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$3 a.u. (Srivastava, 2019). This is a field-induced activation of second-order response through controlled symmetry lowering.

A second route is structural reconfiguration. Stoichiometric $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$4 is amorphous and statistically centrosymmetric, so its bulk $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$5 vanishes. Electrical poling at high temperature, with an applied DC field of approximately $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$6, aligns Si–N bonds or local polar units and freezes in a noncentrosymmetric configuration. The resulting poled material exhibits a bulk $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$7, inferred from electro-optic measurements to be approximately $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$8 (Zhang et al., 2022). In that case the symmetry class is modeled as analogous to $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \cdots$$9, and the induced electro-optic response is consistent with tensor components $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$0 and $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$1.

A third route does not require breaking structural centrosymmetry at all. In centrosymmetric hyperbolic media, the gradient term $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$2 is inversion-even, because both $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$3 and $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$4 pick up a minus sign under inversion. Consequently, a second-order polarization proportional to $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$5 is symmetry-allowed even in a centrosymmetric bulk medium (Narimanov, 22 May 2026). This mechanism is central to recent work on hyperbolic media and is distinct from ordinary bulk $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$6 in noncentrosymmetric crystals.

3. Microscopic mechanisms and quantum-geometric origins

A central distinction in the literature is between intrinsic and extrinsic mechanisms. Intrinsic mechanisms are tied to equilibrium band geometry or linear response tensors. Extrinsic mechanisms arise from disorder-mediated dynamics such as skew scattering and side jump.

In inversion-broken but time-reversal-symmetric conductors, the intrinsic nonlinear Hall effect is closely related to the Berry curvature dipole. The semiclassical framework summarized in the data block defines a Berry curvature dipole tensor

$P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$7

and the second-harmonic Hall current at lowest order in the field is controlled by this dipole (Zhang et al., 2021). In the Matsubara framework, the corresponding second-order Hall conductivity appears through derivatives of an interband Berry curvature object $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$8, yielding $P_i^{(2)} = \varepsilon_0 \chi^{(2)}_{ijk}E_jE_k$9, $\chi^{(2)}$0, and the first-order SHG conductivity $\chi^{(2)}$1 in terms of Berry-curvature-dipole structures (Zhang et al., 9 Jul 2025).

The quantum metric provides a distinct intrinsic channel. In the Matsubara theory of nonlinear Ohmic electromagnetic response, the order-$\chi^{(2)}$2 Ohmic SHG conductivity is

$\chi^{(2)}$3

so the intrinsic nonlinear Ohmic response is governed by the fully symmetrized normalized quantum metric dipole (Zhang et al., 10 May 2026). In a related non-Hermitian framework, the narrow-wavepacket second-order DC conductivity contains an intrinsic, scattering-time-independent term

$\chi^{(2)}$4

so the real part of the non-Hermitian quantum metric controls an intrinsic nonlinear conductivity even in open systems with a spectral line gap (Chen et al., 15 Sep 2025).

Extrinsic mechanisms dominate in several moiré materials. In twisted bilayer graphene aligned with hBN, the $\chi^{(2)}$5 point-group symmetry forbids the intrinsic Berry-curvature-dipole contribution, so the observed giant nonlinear Hall effect is extrinsic. The measured second-order Hall conductivity reaches $\chi^{(2)}$6, and the effect is attributed to skew scattering from static Coulomb impurities at low temperature together with phonon skew scattering at elevated temperature (Duan et al., 2022). Twisted double bilayer graphene shows an even larger extrinsic second-order conductivity, $\chi^{(2)}$7, near a mid-band van Hove singularity, again tied to extrinsic side-jump and skew-scattering channels (Ahmed et al., 8 Jul 2025).

Strongly correlated systems add another microscopic layer. In a noncentrosymmetric Kondo lattice with Rashba-type spin-orbit coupling, ferromagnetism is required for a finite second-order conductivity, and the response becomes finite only perpendicular to the ferromagnetic magnetization. In the low-density approximation, the second-order conductivity matrix scales as

$\chi^{(2)}$8

so for $\chi^{(2)}$9 only the $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$0-direction nonlinear response survives (Shinada et al., 2021).

4. Centrosymmetric hyperbolic media and gradient-induced bulk response

The recent hyperbolic-media mechanism is unusual because it produces a bulk second-order response in a centrosymmetric material without invoking crystal-potential anharmonicity. The starting point is a more complete expansion of the polarization,

$$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$1

which is explicitly gradient-dependent (Narimanov, 22 May 2026). The tensor $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$2 is determined entirely by linear susceptibilities through

$$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$3

with $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$4 the atomic number density and $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$5 the electron charge (Narimanov, 22 May 2026). The paper explicitly states that this mechanism does not rely on any anharmonicity of the crystalline potential and is entirely governed by the linear response of the medium.

Hyperbolic media enable this response because their dielectric tensor has opposite signs along different principal directions, for example

$$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$6

with $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$7 and $$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$8 of opposite sign. The resulting TM isofrequency surface,

$$\mu_i = \mu_i^{(0)} + \sum_j \alpha_{ij}E_j + \frac{1}{2}\sum_{jk}\beta_{ijk}E_jE_k + \frac{1}{6}\sum_{jkl}\gamma_{ijkl}E_jE_kE_l + \cdots,$$9

is a hyperboloid, and large-$\beta_{ijk}$0 modes are propagating rather than evanescent (Narimanov, 22 May 2026). Consequently, the internal field can vary on deeply subwavelength scales, so $\beta_{ijk}$1 is no longer limited by $\beta_{ijk}$2. This activates the otherwise negligible $\beta_{ijk}$3 nonlinearity.

For a structured hyperbolic slab, the far field sees an effective bulk second-order susceptibility defined by

$\beta_{ijk}$4

For the dominant $\beta_{ijk}$5 component, the explicit result is

$\beta_{ijk}$6

with defect size $\beta_{ijk}$7, spacing $\beta_{ijk}$8, and slab thickness $\beta_{ijk}$9 (Narimanov, 22 May 2026). Using $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$0, $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$1, $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$2, the estimated magnitude is

$$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$3

which is of the same order as ADP and KDP values quoted in the paper (Narimanov, 22 May 2026). This supports the claim that centrosymmetric hyperbolic media such as hBN can exhibit bulk second-harmonic, sum-frequency, and difference-frequency generation efficiencies comparable to established nonlinear crystals.

A plausible implication is that this mechanism broadens the meaning of “bulk” second-order response. The effective $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$4 is not a bare crystal parameter but an emergent property of the hyperbolic medium, subwavelength couplers, and the resulting field distribution (Narimanov, 22 May 2026).

5. Materials platforms and representative magnitudes

The data block documents several experimentally or theoretically important material classes in which second-order NLER is large, tunable, or unusually structured.

Platform	Mechanism emphasized	Representative magnitude
hBN hyperbolic medium	Gradient-induced effective bulk $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$5 in centrosymmetric medium	$$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$6 (Narimanov, 22 May 2026)
Poled $$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$7	Electrical poling and bond alignment	$$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$8 (Zhang et al., 2022)
$$j_\alpha = \sigma^{(1)}_{\alpha\beta}E_\beta + \chi^{(2)}_{\alpha\beta\gamma}E_\beta E_\gamma + \dots,$$9 under OEEF	Field-induced symmetry breaking and gap reduction	$$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$0 a.u. at $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$1 a.u. (Srivastava, 2019)
TBG/hBN	Extrinsic skew scattering	$$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$2 (Duan et al., 2022)
tDBLG	Extrinsic NLER near vHSs	$$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$3 (Ahmed et al., 8 Jul 2025)
BLG	Lifshitz-transition-sensitive NLER	$$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$4 exceeds $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$5 at 3 K (Ahmed et al., 8 Jul 2025)

In poled $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$6, the induced $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$7 manifests as a high-speed Pockels response with $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$8 up to approximately $$j_i = \sigma^{(1)}_{ij}E_j + \sigma^{(2)}_{ijk}E_jE_k + \dots$$9 and a 3 dB electro-optic bandwidth of at least $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$0, in contrast to the $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$1 carrier-dominated response of non-poled devices (Zhang et al., 2022). In $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$2, the same external field that enhances $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$3 also reduces the HOMO–LUMO gap from $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$4 eV to $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$5 eV, consistent with the paper’s two-level argument $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$6 (Srivastava, 2019).

In bilayer graphene, NLER becomes a probe of Fermi-surface topology. The second-order conductivity changes sign near Lifshitz transitions and remains diagnostic even at $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$7 K (Ahmed et al., 8 Jul 2025). In twisted double bilayer graphene, a “cascade of singularities” appears as sharp sign reversals and extrema in $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$8 near multiple van Hove singularities across several moiré bands (Ahmed et al., 8 Jul 2025). These observations suggest that second-order conductivity is exceptionally sensitive to Fermi-surface reconstructions.

6. Many-body, Matsubara, and coarse-grained response theory

A general many-body formulation of second-order NLER requires a causal three-point response function. The imaginary-time formalism paper proves that fully causal $$P_i^{(2)}(2\omega)=\varepsilon_0\sum_{jk}\chi^{(2)}_{ijk}(2\omega;\omega,\omega)E_j(\omega)E_k(\omega),$$9-th order response functions can be obtained by analytic continuation of Matsubara $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$00-point functions, with the main theorem

$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$01

implemented at the level of the Lehmann representation (Sinha et al., 26 Jun 2025). For $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$02, this means the physical causal second-order conductivity or susceptibility can be computed from an imaginary-time three-point correlator and analytically continued to $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$03.

The same work provides the time-domain Kubo structure through nested commutators,

$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$04

and derives a high-frequency sum rule for $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$05-th harmonic generation,

$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$06

which for $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$07 constrains the asymptotic second-harmonic response (Sinha et al., 26 Jun 2025). This provides a rigorous route from Matsubara diagrammatics to measurable second-order NLER in interacting and disordered systems.

At the level of metamaterial homogenization, second-order NLER is richer than a purely electric $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$08. In a bi-anisotropic medium one may define

$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$09

to include electric, magnetic, and magneto-electric second-order couplings. The total number of effective second-order susceptibilities is $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$10, and a nonlinear transfer-matrix retrieval scheme can recover them from a sufficient set of sum-frequency-generation measurements or simulations (Larouche et al., 2017). This is a formal generalization of NLER from conventional dielectric polarization to full effective electromagnetic response.

A different formal generalization arises in coarse-grained response theory. For equilibrium stochastic systems, the second-order response of an observable can be written as

$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$11

where $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$12 is the time-antisymmetric entropy-production part and $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$13 is the time-symmetric dynamical-activity part (Müller et al., 2020). The related path-space formulation emphasizes that linear response depends only on entropy production, while second-order response explicitly probes frenetic, dynamical details of the system (Basu et al., 2014). This suggests that in coarse-grained electrical systems, second-order response can reveal kinetic structure invisible in linear transport.

7. Experimental signatures, selection rules, and open directions

A standard experimental signature of second-order NLER under AC drive $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$14 is a voltage or current at $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$15. In twisted bilayer graphene, the transverse voltage $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$16 scales linearly with $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$17 and does not change sign when current direction is reversed, consistent with a second-order signal (Duan et al., 2022). In bilayer graphene and tDBLG, the second-harmonic longitudinal and transverse voltages are extracted to obtain $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$18 and $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$19, with sign reversals near Lifshitz transitions or van Hove singularities (Ahmed et al., 8 Jul 2025).

Selection rules depend sharply on time-reversal and inversion symmetry. In the semiclassical harmonic-order framework, with time-reversal symmetry preserved and inversion symmetry broken, even harmonics of the charge current are purely transverse while odd harmonics are longitudinal (Zhang et al., 2021). With both time-reversal and inversion preserved, even harmonics vanish. With both broken, all harmonic orders of charge and spin currents can in principle appear (Zhang et al., 2021). These results organize when second-order NLER is allowed and whether it should be Hall-like or longitudinal.

The recent literature also separates Hall-like and Ohmic second-order responses. In the Matsubara theory of nonlinear Ohmic electromagnetic response, the order-$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$20 Ohmic conductivity vanishes for both SHG and the bilinear magnetoelectric effect, while the intrinsic order-$$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$21 contribution survives and is governed by the normalized quantum metric dipole (Zhang et al., 10 May 2026). In the broader quantum theory of nonlinear electromagnetic response, Berry-curvature-dipole terms dominate time-reversal-symmetric nonlinear Hall effects, while quantum-metric-dipole terms govern intrinsic nonlinear Hall effects in time-reversal-breaking systems (Zhang et al., 9 Jul 2025).

Several open directions are explicit in the source material. The hyperbolic-media work notes that once an effective $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$22 is established, conventional phase matching can be used to enhance output (Narimanov, 22 May 2026). The poled $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$23 work points to higher poling fields and lower-loss geometries as routes toward larger $$\mu_i = \mu_i^{(0)} + \alpha_{ij}E_j + \frac{1}{2}\beta_{ijk}E_jE_k + \cdots$$24 (Zhang et al., 2022). The moiré-material studies suggest that NLER is a reliable tool for locating Lifshitz transitions and Fermi-surface reconstructions under time-reversal-symmetric conditions (Ahmed et al., 8 Jul 2025). A plausible implication is that second-order conductivity will continue to function both as a spectroscopy of electronic quantum geometry and as a device figure of merit for rectification, electro-optic modulation, and frequency conversion.