Abstract

Realizing scalable quantum photonic networks requires precise, dynamically tunable control over light–matter interactions at subwavelength scales. This capability is currently bottlenecked by the static nature of conventional planar metamaterials, which limits real-time quantum state routing. This work introduces a dynamically reconfigurable planar metamaterial platform featuring an atomically thin graphene monolayer to achieve deterministic, continuous tuning of quantum emitter coupling environments. By exploiting transition boundary conditions in rigorous full-wave finite element simulations using COMSOL Multiphysics Software, we map the evolution of highly localized near-field resonant modes, identifying a dominant localized surface plasmon resonance at 4.3 THz. We demonstrate that actively modulating the graphene chemical potential (Ef) from 0.05 eV to 0.6 eV provides external control over the Local Density of Optical States (LDOS). This control enables a transition from weak uncoupled states to a near-perfect critically coupled light-matter absorption regime (98%). Crucially, the spatial field distribution transitions from an edge-confined mode to a symmetric multi-polar profile at 5 THz, culminating in a complete spatial mode inversion at 6 THz. This framework provides the missing design rules for actively steering complex quantum states, establishing a concrete pathway for fully integrated, dynamically programmable quantum photonic circuits.

1. Introduction

The generation, routing, and control of complex quantum states—such as single-photon states, squeezed light, and entangled photons—are foundational to quantum computing, secure quantum key distribution, and high-precision quantum metrology by ¹. Traditional quantum optical setups rely on bulky crystals or high-quality-factor conventional cavities to enhance light-matter interactions observed by ². However, these approaches face severe challenges regarding scalability, physical footprint, and rigid operating bandwidths that hinder real-time quantum state routing ³.

Planar metamaterials have emerged as a disruptive paradigm to overcome these limitations. By engineering subwavelength-scaled, artificially patterned structures (meta-atoms), researchers can manipulate the phase, amplitude, and quantum vacuum fluctuations of light at will like the case of ⁴. When combined with two-dimensional (2D) materials like graphene, these planar platforms gain crucial dynamic tunability and high field confinement, making them ideal for manipulating light-matter interactions at the single-quantum limit as seen by ⁵.

Despite these advances, a critical bottleneck remains: current tunable metamaterial designs lack a scalable, continuous, and deterministic mechanism to dynamically control the coupling efficiency of quantum states ⁶. Most existing models rely on static geometries that cannot adapt to real-time quantum state routing or varying emitter frequencies ⁷. To address this gap, this study introduces a comprehensive, dynamically reconfigurable planar metamaterial framework optimized for deterministic quantum emitter coupling and light-matter interaction control.

Rather than treating the platform as a static system, we deploy a transition boundary condition to explicitly map how a 1 nm graphene monolayer continuously modulates the near-field quantum environment. Specifically, we reveal how actively tuning the graphene chemical potential from 0.05 to 0.6 eV serves as a direct external dial to control localized surface plasmon resonances near 4.3 THz, governing the local density of optical states (LDOS) ⁸. By evaluating the precise correlation between unit cell geometry, material response, and localized field enhancement, this work establishes the practical design rules required to build integrated, dynamically programmable quantum photonic circuits.

2. Theoretical and Numerical Methods

The computational and theoretical workflow engineered to model the quantum-cooperative planar metamaterial bridges geometric design, material characterization, boundary approximation, and full-wave electromagnetic solutions. The complete schematic sequence from unit cell layout to frequency domain quantum-field environment calculation is systematically detailed in Figure 1.

The design begins by defining the subwavelength unit cell geometry (Figure 1a) to optimize localized light-matter interactions. The meta-atom comprises a patterned concentric split-ring-like resonator configuration. This structure is patterned atop a dielectric spacer layer, which rests on a continuous substrate backbone. The structural stack consists of three functional regimes (Figure 1b): an active monolayer graphene sheet, a dielectric Silicon Dioxide) spacer layer, and a Silicon metal substrate base. Each domain is mathematically is defined by assigning its constitutive parameters: frequency-dependent complex permittivity (), permeability and (electrical conductivity (σ) which dictates the electronic backing of the quantum control interface.

Massive mesh counts caused by modeling an atom-thin layer of 1 nm in a macroscopic mesh are avoided by modeling graphene as an ultra-thin conductive sheet (Figure 1c). We implemented a Transition Boundary Condition (TBC)) ⁹, governed by the surface current density in equation (1):

(1)

where is the surface current density vector, represents the tangential component of the localized electric field of the interface.

The complex surface conductivity of the monolayer graphene sheet (σ(ω)) governs the light-matter coupling parameters and is modeled using the comprehensive Kubo formulation by ¹⁰. The complex conductivity is dictated by the summation of intraband electron-phonon scattering and interband electron-electron transitions ¹¹, given by:

σ(ω)= σ_intra(ω)+ σ_inter(ω)(2)

The intraband scattering contribution is mathematically evaluated using:

(3)

while the interband electron transition contribution is approximated via:

(4)

where e is the is the elementary electron charge, is the reduced plank constant is the Boltzmann constant, is the angular operating frequency, T is the active operating temperature (300 K), is the quantum periodic time (1x 10^-13 s) and is the dynamically tuned chemical potential.

To rigorously solve the frequency-domain vector wave equation that dictates the spatial mapping of light-matter interactions, precise boundary conditions, waveguide ports, and localized mesh refinement strategies were enforced. Perfect electric conductor (PEC) was assigned along parallel walls perpendicular to the electric field vector by enforcing:

(5)

Similarly, perfect magnetic conductor (PMC) was enforced on the orthogonal pair of lateral boundaries to establish a perfect magnetic wall ^{12, 13}:

(6)

The vector wave equation is solved via ¹⁴:

(7)

where is the electric field vector, is the free-space wave number, and relative permittivity of the material.

Port 1 (Incident Port) launches a normally incident wave along the negative z-direction, while Port 2 (Reference Port) captures transmitted power. Scattering parameters are gathered to compute macroscopic optical coefficients (R, T, and A).

To link these macroscopic results to quantum state modulation, the Local Density of Optical States (LDOS), ρ(r, ω), which determines the emission decay rate of a quantum emitter via the Purcell effect as seen by ¹⁵, is evaluated directly from the green's tensor of the classical electric field enhancement:

(8)

where is the free-space LDOS, and Gsc is the scattered green's tensor, proportional to our calculated local electric field norm |E|. A global 3D unstructured tetrahedral mesh discretizes the volume (Figure 2).

The full-wave numerical modeling of the graphene-based planar metamaterial is governed by the specific geometric, thermodynamic, and material constants detailed in Table

3. Results and Discussion

We demonstrated the active control of light-matter interactions and quantum emitter environments by evaluating the spectral absorptance (A) across a frequency range from 0.1 to 6.0 THz as a function of graphene's Fermi energy E_f. As illustrated in Figure 3, modulating E_f from 0.05 eV to 0.6 eV reveals a highly tunable behavior transitioning from weak light-matter coupling to near-perfect critical coupling, which agrees with ¹⁶.

At low levels, graphene behaves like a lossy dielectric with low carrier concentration. For E_f = 0.05 eV, light-matter coupling remains weak, and absorptance stays below 13% across the entire band. At E_f = 0.1 eV, a broad resonance begins to appear near 2.2 THz, peaking at 23%. At E_f = 0.2 eV, a stable resonant profile forms with a peak absorptance of 52% at 2.85 THz. At E_f = 0.5 eV, the system reaches a near-perfect absorption peak of approximately 98% at 4.3 THz due to the critical coupling condition, matching trends in hybrid architectures.

A continuous blueshift of the resonance peak frequency is observed as E_f increases. By varying the gate voltage to tune E_f, the resonance frequency shifts dynamically across a wide bandwidth from 2.2 THz to 4.3 THz. This features a dynamic mechanism to tune the metamaterial to match specific quantum emitter emission frequencies.

function of incident frequency (THz) for varying graphene Fermi energy levels

(E_f = 0.05 eV to 0.6 eV).

At low levels, graphene behaves like a lossy dielectric with low carrier concentration. For Ef = 0.05 eV, absorptance remains below 13% across the entire band. At Ef = 0.1eV, a broad resonance begins to appear near 2.2 THz peaking at 23%. At Ef = 0.2 eV, a resonant profile form with a peak absorptance of 52% at 2.85 THz. At Ef = 0.5 eV, the system reaches a near-perfect absorption peak of approximately 98% at 4.3 THz due to the critical coupling condition which is similar to what was gotten by ¹⁶.

A continuous blueshift of the resonance peak frequency is observed as E_f increases. By varying the gate voltage to tune E_f, the resonance frequency shifts dynamically across a wide bandwidth from 2.2 THz to 4.3THz. In summary, when E_f increases, intra-band transitions dominate the surface conductivity.

We investigated how the metamaterial enhances high-precision light-matter interactions; the global integrated electric field norm was tracked as a function of frequency. As shown in Figure 4, the electric field amplitude exhibits a clear resonant peak.

The curve rises steadily from near-zero up to a maximum localized value at approximately 3.3 THz. This peak marks the global localized surface plasmon resonance of the structural array, which agrees with ¹³. Beyond 3.3 THz, the field profile encounters a slight dip and settles into a plateau around 4.5 to 6.0 THz, mapping the system's transition into higher-order multi-polar regimes. This pronounced spectral field enhancement provides the physical mechanism required to increase the local density of optical states (LDOS), delivering the high-precision coupling needed to manipulate and route localized quantum states.

The curve rises steadily from near-zero up to a maximum localized value at approximately 3.3 THz. This peak marks the global localized surface plasmon resonance of the structural array which is in agreement with ¹⁷. Beyond 3.3 THz, the field profile encounters a slight dip and settles into a plateau around 4.5 to 6 THz, mapping the system's transition into higher-order multi-polar regimes similar to results by ¹⁸. This pronounced spectral field enhancement provides the physical mechanism required to increase the local density of optical states (LDOS), delivering the high-precision coupling needed to manipulate quantum emitters as seen by ¹⁹.

To unveil the microscopic origins of this quantum state control, 2D spatial profiles of the surface electric field norm were mapped at different frequencies across the band.

At 1 THz, the field maximum is limited to 4.97 x 10⁶ V/m, representing a weak light-matter interaction regime. At 3 THz, aligning with the strong resonance slope, the field is amplified significantly to 14.5 x 10⁶ V/m, localizing intensely along the right structural boundary and establishing a spatial gradient for asymmetric quantum state interactions.

At 1 THz as shown in Figure 5a, the field maximum is limited to 4.97 x 10⁶ V/m. At 3 THz in Figure 5c, aligning with the strong resonance slope, the field is amplified significantly to 14.5 x 10⁶ V/m, localizing intensely along the right structural boundary.

The peak field intensity reaches its absolute maximum of 1.72 x 10⁷ V/m at 4 THz, matching the near-perfect absorption threshold under high gate voltage conditions. This maximum spatial field confinement creates a deep subwavelength hot-zone that drastically enhances the Purcell factor, enabling strong coupling between the metamaterial and local quantum states. The same trend is also observed in the research by ²⁰.

At 5 THz, the edge-localized gradient transforms into a central oval hotspot (x = 0 um) peaking at 1.03 x 10⁷ V/m, indicating a higher-order symmetric mode similar to results by Jahangiri et al. At 6 THz, the field layout flips completely to the left boundary (x = -2 um) at 1.24 x 10⁷ V/m due to an out-of-phase scattering regime. This complete spatial mode inversion demonstrates that a quantum emitter's decay pathway can be dynamically routed to different spatial channels depending on the excitation frequency ²¹; as shown in Figure 6:

A comparative summary of the near-field performance parameters across the simulated spectrum is organized in Table 3.

4. Conclusion

This work has presented a robust computational modeling framework to engineer strong light-matter interactions using a graphene-integrated planar metamaterial. By employing a transition boundary condition and solving the frequency-domain vector wave equation under exact PEC and PMC conditions, we successfully evaluated both macroscopic optical properties and localized quantum field environments. The platform demonstrates near-perfect absorption of 98% and strong near-field confinement of 1.72 x 10⁷ V/m with a gate-tunable resonant blueshift from 2.2 THz to 4.3 THz. These dynamic spatial and spectral field transformations alter the local density of optical states (LDOS), offering an effective route for chip-scale quantum state generation, quantum emission control, and high-precision quantum optical circuits. Future work should focus on the fabrication and experimental characterization of the proposed planar metamaterial structures via THz time-domain spectroscopy (THz-TDS) to validate these design rules in physical quantum hardware environments.

References