405. Second and Third Harmonic Generation in Topological Insulator-Based van der Waals Metamaterials

Introduction

High-order harmonic generation (HHG) is the nonlinear up-conversion of an electromagnetic signal into integer multiples of the fundamental frequency. In solids, HHG has become an important route to access spectral regions that are difficult to reach with conventional laser sources, particularly in the terahertz (THz) domain where compact, tunable, high-power emitters remain limited. The harmonic spectrum encodes the symmetries, band topology, and carrier dynamics of the material, making HHG not only a source technology but also a probe of condensed-matter physics.

The simplest nonlinear processes are second harmonic generation (SHG) and third harmonic generation (THG). SHG is forbidden in centrosymmetric media within the electric-dipole approximation, while THG is generally allowed. This symmetry dependence becomes especially powerful in topological materials, where bulk inversion symmetry, surface symmetry breaking, and strong spin-orbit coupling coexist. Topological insulators (TIs) are a particularly rich platform because they host an insulating bulk with inverted band order and metallic surface states protected by time-reversal symmetry. These surface states can strongly alter nonlinear optical response, including harmonic generation.

Recent work on TI-based van der Waals metamaterials has shown that embedding Bi2Se3 or (InxBi(1-x))2Se3/Bi2Se3 heterostructures into resonant split-ring resonator (SRR) arrays can dramatically enhance local THz fields and enable harmonic up-conversion in the 6.4–9.7 THz range. In these structures, the interplay of field confinement, metamaterial resonance, and topological electronic states allows both even and odd harmonics to emerge. The observed harmonics are not simply a consequence of strong driving; they are a direct manifestation of symmetry selection rules in the bulk and at the surface.

Symmetry and nonlinear optics in solids

In nonlinear optics, the polarization response of a medium is expanded as
P = ε0(χ(1)E + χ(2)E^2 + χ(3)E^3 + ...),
where χ(2) governs SHG and χ(3) governs THG. The existence and magnitude of these susceptibilities depend on crystal symmetry, electronic structure, and excitation frequency. In a centrosymmetric crystal, the second-order susceptibility vanishes in the electric-dipole approximation because inversion symmetry requires P(E) = -P(-E), eliminating even-order terms. Consequently, SHG is usually absent in the bulk of centrosymmetric materials, though it may arise at surfaces, interfaces, or from higher multipoles.

THG, by contrast, is allowed in centrosymmetric systems because χ(3) is inversion-even. In many materials, THG is therefore the dominant harmonic process under intense optical driving. However, the amplitude, phase, and spectral content of THG still depend strongly on band structure, carrier scattering, and resonant enhancement.

Topological insulators complicate this picture in a useful way. Their bulk often retains inversion symmetry, but the surface necessarily breaks it. The surface electronic states are also spin-momentum locked, meaning that the carrier dynamics are constrained by the topological band structure. As a result, a TI can exhibit odd harmonics from both bulk and surface, while even harmonics can emerge specifically from symmetry-broken surface contributions or from interfaces in heterostructures. This makes TIs ideal systems for disentangling bulk and surface nonlinearities.

Topological insulators and van der Waals heterostructures

Bi2Se3 is one of the most studied three-dimensional topological insulators. It has a narrow band gap in the bulk and a single Dirac cone on the surface. The surface state is protected against nonmagnetic disorder, but it remains sensitive to symmetry breaking, strain, thickness, and interface engineering. In thin films and van der Waals stacks, the relative contributions of top and bottom surfaces, bulk carriers, and interface states can be tuned systematically.

The alloy (InxBi(1-x))2Se3 provides an additional degree of control. By varying the indium concentration x, one can modify the band inversion strength, bulk carrier density, and the degree of symmetry breaking. When combined with Bi2Se3 in a van der Waals heterostructure, the resulting electronic structure can support hybridized topological states and modified optical selection rules. Such heterostructures are particularly attractive for nonlinear optics because they preserve sharp interfaces and allow stacking without the constraints of lattice matching typical of epitaxial growth.

In the context of HHG, these materials offer three important advantages: strong spin-orbit coupling, tunable surface/bulk balance, and compatibility with metamaterial field enhancement. The van der Waals nature of the heterostructure reduces interfacial degradation and enables controlled thicknesses on the order of a few nanometers, where surface contributions become comparatively large. This is essential for observing weak harmonic signals in the THz regime.

Metamaterial field enhancement with split-ring resonators

The central experimental challenge in observing SHG and THG in TIs at THz frequencies is the need for sufficiently strong driving fields. Harmonic generation is a nonlinear process, so the output scales nonlinearly with the incident field amplitude. In practice, the field strengths required for measurable signals can be difficult to achieve without damaging the sample or introducing excessive thermal effects.

Split-ring resonators provide a solution. SRRs are subwavelength metallic resonant structures that concentrate incident THz radiation into nanoscale gaps, generating intense localized electric fields. Single split-ring resonators primarily support a fundamental LC resonance, while double split-ring resonators can offer additional modes and stronger near-field confinement. When a TI thin film or heterostructure is embedded in or placed near the resonator gap, the local field experienced by the material can be enhanced by orders of magnitude relative to the incident field.

This local enhancement lowers the effective threshold for nonlinear conversion and can selectively amplify surface-sensitive contributions. Because the field is highly nonuniform, the metamaterial also breaks spatial symmetry, further relaxing selection rules and enabling processes that may be weak or forbidden in the bare material. In this sense, the metamaterial is not merely a passive antenna; it is an active symmetry-engineering platform.

Second harmonic generation in TI metamaterials

SHG in topological insulators is particularly interesting because it provides a direct signature of inversion symmetry breaking. In a bulk centrosymmetric TI, SHG should be suppressed, but at the surface the inversion symmetry is broken by construction. The surface nonlinear polarization can therefore radiate at 2ω even if the bulk remains silent. In a thin film, the relative contribution of the top and bottom surfaces may depend on film thickness, substrate asymmetry, and the presence of capping layers or heterostructure interfaces.

In Bi2Se3-based van der Waals metamaterials, SHG can be enhanced by resonant coupling to the SRR mode. The resonator concentrates the fundamental field at the TI surface, where the nonlinear response is strongest. Because the surface state is topological and spin-momentum locked, the nonlinear current can exhibit unusual angular dependence and polarization selection rules. This is distinct from conventional semiconductors, where SHG is typically governed by crystal point group symmetry alone.

The observation of SHG in the 6.4–9.7 THz regime is especially significant because it demonstrates frequency doubling in a spectral region where many materials are weakly nonlinear and where conventional sources are scarce. For a 3.2–4.85 THz pump, the second harmonic lies in the 6.4–9.7 THz band. Such conversion opens access to higher-frequency THz spectroscopy, imaging, and sensing applications. It also provides a testbed for studying how surface topological states contribute to even-order nonlinearities.

Third harmonic generation and bulk nonlinear response

THG is expected in centrosymmetric materials, so it serves as a natural benchmark for nonlinear conversion in TIs. In Bi2Se3, the bulk contribution to THG can arise from intraband acceleration, interband polarization, and higher-order current nonlinearities. Under strong THz driving, carriers are driven far from equilibrium, and their motion within the nonparabolic band structure leads to anharmonic current components at 3ω.

The presence of a topological surface state modifies this behavior in several ways. First, the surface state can contribute additional THG through its linear Dirac dispersion and spin-momentum locking. Second, because the surface lacks inversion symmetry, it can generate both even and odd harmonics, but odd harmonics remain allowed in both bulk and surface channels. This means THG often contains mixed contributions that are difficult to separate experimentally unless the symmetry and geometry are carefully controlled.

In TI-based metamaterials, THG can be further enhanced by the resonant response of the SRR. The resonator not only boosts the local field but can also shape the phase and polarization of the emitted harmonic. The result is a strong, spectrally narrow THG signal that reflects both the nonlinear susceptibility of the TI and the electromagnetic mode structure of the metamaterial.

Even and odd harmonics from bulk and surface states

A key feature of topological insulator nonlinear optics is the coexistence of even and odd harmonics. Odd harmonics are compatible with both centrosymmetric bulk and surface contributions, while even harmonics generally point to inversion symmetry breaking. In a TI, this symmetry breaking can originate from the topological surface itself, the interface with a substrate, strain gradients, or the asymmetry introduced by heterostructure layering.

The experimental observation of both SHG and THG in Bi2Se3 and (InxBi(1-x))2Se3/Bi2Se3 metamaterials therefore indicates that the harmonic emission is not purely bulk-driven. Instead, the data support a picture in which the bulk contributes primarily to odd-order processes, while the surface and interface states enable even-order emission. This distinction is crucial because it provides a route to separate topological surface response from bulk background, a longstanding challenge in TI spectroscopy.

The coexistence of even and odd harmonics also suggests that the nonlinear response is sensitive to the microscopic symmetry of the electronic states. In particular, the surface Dirac fermions may respond differently to the driving field depending on whether the field couples predominantly to in-plane or out-of-plane motion. The SRR geometry can therefore be used to tailor the balance between surface and bulk nonlinear currents.

Physical mechanisms of harmonic generation

Several microscopic mechanisms can contribute to HHG in TIs. In the intraband picture, the electric field accelerates carriers within a nonlinear band dispersion, producing anharmonic currents. In the interband picture, the field drives coherent polarization between valence and conduction states, and the time-dependent polarization radiates at harmonic frequencies. In real materials, both mechanisms may coexist, especially under strong THz excitation.

For TIs, the surface Dirac cone introduces a nearly linear dispersion that can support strong nonlinear currents when driven far from equilibrium. The spin texture may suppress certain scattering channels and modify dephasing, potentially increasing harmonic coherence. Meanwhile, the bulk inverted bands and strong spin-orbit coupling can enhance interband matrix elements. In heterostructures, band bending and charge transfer at the interface may further alter the local carrier density and thus the nonlinear response.

The metamaterial resonance adds another layer of complexity. The SRR can produce a highly localized quasi-static field that interacts with the TI over a subwavelength region. This means the harmonic signal may be dominated by nanoscale hotspots rather than the entire illuminated area. As a result, the measured conversion efficiency reflects not only intrinsic material nonlinearities but also the electromagnetic environment.

Spectral window and experimental significance

The 6.4–9.7 THz frequency range is technologically important because it lies in a region where sources and detectors are limited, yet many molecular and solid-state excitations are active. Generating harmonics in this band can enable compact THz spectroscopy for chemical identification, low-energy vibrational modes, and ultrafast carrier dynamics. It also bridges a gap between microwave electronics and mid-infrared photonics.

The fact that the harmonic output can be tuned between even and odd orders is especially valuable. SHG provides a direct probe of surface and interface symmetry, while THG offers a robust signal linked to bulk and surface nonlinearities. By comparing the relative strengths of the second and third harmonics across different TI thicknesses, alloy compositions, and resonator geometries, one can infer the role of topological surface states in the nonlinear optical response.

This capability is important because traditional probes of topological surface states, such as angle-resolved photoemission spectroscopy, require ultrahigh vacuum and are not easily integrated into devices. Harmonic generation, in contrast, is an all-optical technique that can be performed in situ and potentially in operating devices. It therefore offers a path toward nonlinear topological photonics and THz device engineering.

Outlook and applications

The demonstration of SHG and THG in TI-based van der Waals metamaterials suggests several future directions. First, one can optimize the SRR geometry to maximize field confinement and tailor the polarization dependence of the emitted harmonics. Second, one can engineer heterostructures with controlled asymmetry to selectively enhance even-order processes. Third, one can explore gating, temperature, and magnetic perturbations to modulate the relative weights of bulk and surface contributions.

From an applications perspective, these systems could lead to compact THz frequency multipliers, harmonic mixers, and spectroscopic sources. Because the nonlinear response is tied to topological surface states, such devices may also offer robustness against certain forms of disorder. Beyond technology, the ability to generate and detect harmonics in TIs provides a new window into nonequilibrium quantum dynamics, including carrier acceleration, scattering, and topological protection under strong fields.

More broadly, TI-based metamaterials exemplify the convergence of three research areas: topological condensed matter, van der Waals heterostructures, and nonlinear photonics. By combining symmetry-engineered resonators with quantum materials, one can access nonlinear optical processes that are otherwise weak or hidden. The resulting harmonic emission is not merely an optical output; it is a diagnostic of how topology, symmetry, and strong-field electrodynamics cooperate in solids.

Conclusion

Second and third harmonic generation in topological insulator-based van der Waals metamaterials represents a compelling platform for both fundamental physics and THz photonic applications. The bulk of a TI, with its centrosymmetric inverted band structure, naturally favors odd harmonics such as THG, while the topological surface and interfaces break inversion symmetry and enable SHG. When these materials are embedded in split-ring resonator arrays, the local THz field is strongly amplified, allowing harmonic up-conversion into the 6.4–9.7 THz range.

The observed even and odd harmonics arise from a combination of bulk centro-symmetry, symmetry breaking at the surface, and electromagnetic enhancement by the metamaterial. This makes the system uniquely suited to disentangle bulk and surface nonlinearities in a regime where direct experimental access has been limited. As a result, TI-based van der Waals metamaterials offer not only a route to efficient THz frequency conversion, but also a powerful probe of topological electronic states under intense optical driving.