EFISH: Electric-Field Induced Second-Harmonic Generation

• EFISH is a nonlinear optical process where a static electric field breaks inversion symmetry in centrosymmetric media, enabling effective second-order responses.
• It leverages mixing via χ(3) and applied dc fields to induce a χ(2)-like behavior, as demonstrated in silicon waveguides, bilayer MoS₂, and other platforms.
• Applications include electrically tunable SHG devices, spatial electric field diagnostics, and integrated nonlinear photonics through optimized field and mode engineering.

Searching arXiv for recent and foundational EFISH papers to ground the article. Electric-field-induced second-harmonic generation (EFISH) is the activation, enhancement, or modulation of second-harmonic generation (SHG) by a static or low-frequency electric field. In the standard picture, the applied field breaks inversion symmetry and converts a third-order nonlinear response into an effective second-order one, allowing a centrosymmetric medium with vanishing bulk electric-dipole $\chi^{(2)}$ to radiate at $2\omega$. EFISH therefore provides an electrically tunable route to SHG in platforms that would otherwise forbid bulk second-order conversion, including silicon, bilayer transition-metal dichalcogenides, gases, ferroelectrics, and several classes of nanophotonic structures (Fan et al., 12 Jul 2025, Timurdogan et al., 2016, Klein et al., 2017).

1. Constitutive description and weak-field scaling

Ordinary SHG is described by the second-order polarization

$$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$

or, equivalently, by

$$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$

In a centrosymmetric crystal, inversion symmetry forces the bulk electric-dipole $\boldsymbol{\chi}^{(2)}$ to vanish, so bulk SHG is forbidden in the electric-dipole approximation (Klein et al., 2017, Fan et al., 12 Jul 2025).

EFISH bypasses that restriction by mixing the optical field with a dc field $\mathbf{E}_{\mathrm{DC}}$ through $\boldsymbol{\chi}^{(3)}$. A standard form is

$$\mathbf{P}_{2\omega} = \left( \boldsymbol{\chi}^{(2)} + 3\,\boldsymbol{\chi}^{(3)}:\mathbf{E}_{\mathrm{DC}} \right) :\mathbf{E}_\omega\otimes\mathbf{E}_\omega,$$

which reduces, for a centrosymmetric medium, to an induced effective second-order response $$\boldsymbol{\chi}^{(2)}_{\mathrm{eff}} \propto \boldsymbol{\chi}^{(3)}:\mathbf{E}_{\mathrm{DC}}$$ (Fan et al., 12 Jul 2025). In silicon waveguides this relation is written explicitly as

$$\chi_{xxxx}^{(2)} = 3\chi_{xxxx}^{(3)} E_{\mathrm{DC}},$$

emphasizing that the field-induced $2\omega$0 is proportional to the large intrinsic $2\omega$1 of silicon (Timurdogan et al., 2016).

In the perturbative weak-field regime, the SHG intensity follows the standard EFISH scaling

$2\omega$2

This quadratic dependence is observed, for example, in 2H bilayer MoS$2\omega$3, where the second-harmonic intensity varies approximately as $2\omega$4 or $2\omega$5, and in gas-phase EFISH diagnostics, where the second-harmonic intensity is proportional to the square of the external electric field (Klein et al., 2017, Cui et al., 2019). Later work, however, shows that this scaling is not universal outside the perturbative regime.

2. Symmetry breaking and microscopic mechanisms

The defining operation in EFISH is inversion-symmetry breaking by an electric field. In a classical description of centrosymmetric media, a dc field adds a linear term to an initially even potential, shifts the equilibrium position, and produces an effective cubic term in the expansion around the shifted minimum; in the Simplified Bond Hyperpolarizability Model, that $2\omega$6 term “immediately leads to second harmonic generation” (Alejo-Molina et al., 2015). In tensor language, the fourth-rank $2\omega$7 contracted with the dc-field direction becomes an effective third-rank EFISH tensor (Alejo-Molina et al., 2015).

In layered materials, the same principle is realized through field-controlled structural asymmetry. A 2H-stacked MoS$2\omega$8 bilayer is inversion symmetric because the two layers are rotated by $2\omega$9 relative to each other. A perpendicular dc field breaks that inversion symmetry and turns SHG on, or strongly enhances it, in the bilayer. The resulting response is not purely geometric: the observed spectral dependence reflects bandstructure and wave-function admixture, and interlayer hybridization is essential. When the two layers are artificially decoupled, the SHG signal is suppressed by at least three orders of magnitude (Klein et al., 2017).

Correlated systems provide a distinct microscopic setting. In the Kitaev honeycomb model, even-order nonlinear susceptibilities vanish unless a reflection-related symmetry $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$0 is broken. A dc electric field renormalizes the Kitaev couplings as $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$1, with $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$2, and thereby enables second-harmonic generation. In that context the low-temperature response is governed by Majorana quasiparticles, while thermally excited visons modify the spectrum at higher temperature through randomization of the $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$3 gauge background (Krupnitska et al., 2023). This shows that EFISH need not be limited to structural asymmetry; it can also probe symmetry breaking in strongly correlated quantum matter.

3. Material platforms and field-engineering strategies

EFISH has been implemented with externally applied fields, built-in depletion fields, space-charge fields, and photoinduced internal fields. Representative realizations span semiconductor waveguides, van der Waals crystals, thin films, and gases (Fan et al., 12 Jul 2025).

Platform	Field engineering	Reported behavior
2H bilayer MoS$$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$4 (Klein et al., 2017)	Microcapacitor geometry; perpendicular field up to $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$5	SHG minimal near $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$6; up to 60-fold enhancement
Silicon ridge waveguide (Timurdogan et al., 2016)	Lateral p-i-n junctions; periodically structured dc field for QPM	$$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$7 at $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$8; $$P^{(2)}_{\alpha}(2\omega)=\varepsilon_0\chi^{(2)}_{\alpha\beta\gamma}(-2\omega;\omega,\omega)E_{\beta}(\omega)E_{\gamma}(\omega),$$9
Bi$$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$0Se$$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$1 thin films (Glinka et al., 2015)	Photoinduced capacitor-type dc field and plasmon-associated dc field	Overall $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$2 SHG enhancement; Lorentz-shaped resonance near $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$3 thickness
Ambient-air corona gap (Cui et al., 2019)	External discharge field in coaxial cylindrical electrodes	Spatial electric-field distribution reconstructed from EFISH signal

These implementations illustrate three recurrent design variables emphasized in the review literature: $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$4 engineering, $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$5 engineering, and optical-mode engineering (Fan et al., 12 Jul 2025). External electrodes provide the most direct route, as in microcapacitors, p-i-n waveguides, and interdigitated structures. Built-in or dynamically generated fields appear in depletion regions, space-charge layers, and photoexcited thin films. In gases and plasmas, the target field itself becomes the object of EFISH measurement rather than merely a control parameter.

4. Resonant, phase-matched, and nonperturbative regimes

Although the weak-field formalism is simple, the observed EFISH response is often highly spectrally structured. In bilayer MoS$$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$6, tuning the pump over $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$7 revealed strong dependence of the SHG response on energy. The largest enhancement appears below the C resonance: at $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$8, detuned by $$\mathbf{P}_{2\omega}=\boldsymbol{\chi}^{(2)}:\mathbf{E}_\omega\otimes\mathbf{E}_\omega.$$9 from $\boldsymbol{\chi}^{(2)}$0, the nonlinear signal increases by about 60-fold under strong field. The response is therefore governed not only by symmetry breaking in the abstract, but by the orbital composition and field sensitivity of the relevant Bloch states (Klein et al., 2017).

Integrated photonics adds a phase-matching constraint. In silicon ridge waveguides, the pump and second-harmonic propagation constants do not naturally satisfy $\boldsymbol{\chi}^{(2)}$1. Quasi-phase matching is achieved by periodically structuring the p-i-n junctions with period $\boldsymbol{\chi}^{(2)}$2, chosen to satisfy $\boldsymbol{\chi}^{(2)}$3. Under reverse bias up to $\boldsymbol{\chi}^{(2)}$4, SHG was observed at multiple pump wavelengths from $\boldsymbol{\chi}^{(2)}$5 to $\boldsymbol{\chi}^{(2)}$6, with maximum efficiency $\boldsymbol{\chi}^{(2)}$7 at $\boldsymbol{\chi}^{(2)}$8 in a $\boldsymbol{\chi}^{(2)}$9 long waveguide (Timurdogan et al., 2016).

Several recent systems depart explicitly from the textbook perturbative picture. In time-varying amorphous-Si metasurfaces, the dc field does more than induce a weak $\mathbf{E}_{\mathrm{DC}}$0: it participates in resonantly enhanced parametric oscillation of a time-dependent resonance. The authors state that, “Despite some similarities to conventional EFISH… the mechanism in SHG enhancement is inherently different.” On resonance, the SHG enhancement factor reaches $\mathbf{E}_{\mathrm{DC}}$1, the modulation depth is $\mathbf{E}_{\mathrm{DC}}$2, and the dc-field dependence becomes super-quadratic (Guo et al., 2020). In an inversion-symmetric MoS$\mathbf{E}_{\mathrm{DC}}$3 homobilayer under resonant excitonic excitation, SHG grows linearly with static field only in the weak-field limit; for sufficiently intense optical excitation, the SHG increases superlinearly with static field amplitude, while THG can also increase, due to an interplay of static and transient Stark shifts, exciton ionization, Wannier–Stark localization, off-resonant Rabi oscillations, and modified interference between optical nonlinearities induced by intraband acceleration (Zuo et al., 7 Nov 2025).

A distinct spatiotemporal generalization appears in diamond, where an applied field induces a moving effective $\mathbf{E}_{\mathrm{DC}}$4 through $\mathbf{E}_{\mathrm{DC}}$5, and spatiotemporal quasi-phase matching predicts second-harmonic output at $\mathbf{E}_{\mathrm{DC}}$6 rather than exactly $\mathbf{E}_{\mathrm{DC}}$7 (Ghalandari, 2015). This suggests that EFISH can be embedded in moving-modulation and nonstationary-matching frameworks.

5. Relation to ordinary SHG, current-induced SHG, and field-controlled analogues

EFISH is distinct from ordinary SHG in noncentrosymmetric media because the latter requires no dc field. It is also distinct from current-induced SHG, in which symmetry breaking occurs in momentum space rather than real space. The review literature formulates the contrast explicitly: EFISH arises from a static electrostatic field, whereas current-induced SHG arises from an asymmetric nonequilibrium carrier distribution (Fan et al., 12 Jul 2025).

The GaAs work on current-induced SHG makes this distinction quantitative. There the induced second-order susceptibility is proportional to current density,

$\mathbf{E}_{\mathrm{DC}}$8

and reversing the current reverses the sign of the SHG signal. In the doped GaAs device studied there, the current contribution was estimated to be more than an order of magnitude larger than the corresponding electric-field-induced contribution expected from conventional EFISH-like mixing (Ruzicka et al., 2011).

A second distinction concerns mechanisms that are only EFISH-like. In III-V wire waveguides aligned with a crystallographic axis, SHG can be enabled by longitudinal electric-field components of guided modes because the zinc-blende tensor permits only mixed-component products containing all three Cartesian directions. That paper explicitly notes an analogy to EFISH “in the broader sense,” but also states that the mechanism is not the usual external-field EFISH where a static $\mathbf{E}_{\mathrm{DC}}$9 mixes with $\boldsymbol{\chi}^{(3)}$0 (Poulvellarie et al., 2020).

Time-resolved ferroelectrics introduce a further interpretive caution. For a low-frequency pulsed electric field $\boldsymbol{\chi}^{(3)}$1, the polarization modulation satisfies approximately $\boldsymbol{\chi}^{(3)}$2, but the SH intensity modulation obeys

$\boldsymbol{\chi}^{(3)}$3

Accordingly, $\boldsymbol{\chi}^{(3)}$4 can depend on both the field and its time derivative, and can even be in antiphase with the pulsed field and with $\boldsymbol{\chi}^{(3)}$5 (Ono, 10 Jul 2025). EFISH measurements of transient order parameters therefore require amplitude-level and phase-sensitive interpretation rather than the simplistic identification of SH intensity with polarization amplitude.

6. Applications, diagnostics, and open problems

EFISH has been proposed and demonstrated as a route to electrically switchable nonlinear devices, compact optical modulators, tunable frequency-doubling elements, integrated nonlinear photonics, and electrically controlled on/off SHG devices (Klein et al., 2017). In silicon photonics, the same p-i-n waveguide platform supports both phase-only modulation through the dc Kerr effect and second-harmonic generation through EFISHG, with CMOS-compatible fabrication and operation spanning near- to mid-infrared wavelengths (Timurdogan et al., 2016).

The concept also extends beyond scalar SHG tuning to symmetry-selective control. In the van der Waals antiferromagnet MnPS$\boldsymbol{\chi}^{(3)}$6, an in-plane electric field transforms the magnetic point group to its unitary subgroup and tunes interference among intrinsic electric-dipole, magnetic-dipole, and field-induced electric-dipole SHG channels. The result is electrically tunable SHG domain contrast and nonreciprocity over a broad spectral range (Wang et al., 2024). This is not merely a bulk $\boldsymbol{\chi}^{(3)}$7 background, but a field-controlled reweighting of symmetry-allowed amplitudes.

As a diagnostic, EFISH is notable for being species-independent, non-resonant, and non-contact. In ambient-air negative corona discharges, it enables spatially resolved field mapping and shows that the electric field at the conductor surface is roughly proportional to corona current density with a negative constant of proportionality. The study concludes that Kaptzov’s assumption is valid only when the discharge current approaches zero or is small (Cui et al., 2019). In discharge measurements more generally, EFISH is a coherent line-of-sight probe rather than a purely local one. The “inverse EFISH problem” arises because the measured second-harmonic power depends on an axial integral containing phase mismatch and the Gouy phase. An operator-learning model termed Decoder-DeepONet was introduced to reconstruct electric-field profiles from EFISH signal profiles; it was trained on 703,682 profile pairs sampled at 109 points and found to remain useful even with incomplete inputs, with Integrated Gradients indicating that practical sampling within about $\boldsymbol{\chi}^{(3)}$8 HWHM of the focus captures the key information (Yang et al., 29 Nov 2025).

The principal open issues are likewise well defined. Review work emphasizes the difficulty of integrating strong, well-controlled static fields into nanophotonic devices, separating EFISH from current-induced SHG, interface SHG, hot-carrier effects, and quantum-confined Stark shifts, mitigating screening and carrier saturation, and balancing resonant enhancement against spectral bandwidth and interaction length. The broader implication is that progress in EFISH depends on simultaneous optimization of material nonlinearity, static-field profile, and optical-mode structure rather than on any one of these factors in isolation (Fan et al., 12 Jul 2025).