Non-Markovian dynamics of a qubit on the zigzag edge of photonic graphene

AI summary

Overview

Research area: Quantum optics and condensed-matter-inspired photonic lattices — specifically non-Markovian open quantum systems, where a single quantum emitter (qubit) interacts with a structured photonic bath.

Technical level: Advanced. The paper is a fully analytical treatment built on tight-binding lattice Green's functions, resolvent operator methods, complete elliptic integrals, and Riemann-sheet analytic continuation. It assumes familiarity with self-energy formalism and complex contour integration.

Scope (one sentence): An exact analytical study of how a qubit coupled to the zigzag edge of a photonic honeycomb lattice decays, showing that the edge's zero-energy flat band replaces the slow branch-cut tails of bulk graphene with damped Rabi oscillations.

What This Paper Is About

When a qubit is placed inside a 2D photonic lattice patterned after graphene, it does not decay exponentially: the Dirac points produce long-lived non-exponential tails, and the qubit can even stay partially trapped. Nearly all of this physics has been studied for a qubit embedded in the bulk. This paper asks what changes when the qubit sits on the zigzag edge instead — a boundary that hosts a zero-energy flat band localized on one sublattice. The goal is to derive the qubit's population dynamics exactly and classify every spectral contribution that shapes it.

Key Contributions

1. A closed-form self-energy for the edge geometry. The authors derive an exact expression for 𝚺(z) = g²⟨vac|a₀₀(z − H_B)⁻¹a₀₀†|vac⟩, written in terms of complete elliptic integrals K̃(z) and Ẽ(z) with parameters α(z), β(z), C(z) and t(z) = (z² − 3)/2, t₀ = −3/2. The expression separates cleanly into a bulk term and a flat-band term Ω²/z, with Ω = g·sqrt(sqrt(3)/π − 1/3).

2. A three-way decomposition of the dynamics. By closing the integration contour in the lower half-plane, the qubit amplitude is split exactly into real bound states (r_BS), unstable complex poles (r_UP), and branch-cut detours c_BC(t). Each contribution is characterized with explicit existence conditions and relative weights.

3. Identification of the flat band as an effective cavity. In the regime |Re{z}| ≪ Ω, the self-energy reduces to 𝚺(z) ≃ Ω²/z, so the pole equation becomes z − Δ − Ω²/z = 0 with roots z± = (Δ ± sqrt(Δ² + 4Ω²))/2. The edge flat band therefore behaves like a single-mode cavity coupled to the qubit.

4. A new long-time scaling law and its suppression. The Dirac-point branch cut produces c_BC(t) ≃ (g/J)²(sqrt(3)/π)(Ωt)⁻⁴, a t⁻⁴ decay that contrasts sharply with the 1/log(t) scaling of bulk graphene — and which the authors show is practically unobservable because it is masked by Rabi oscillations or fractional decay.

Main Findings

• The flat band suppresses branch-cut tails. Unlike a qubit in bulk graphene, the zero-energy flat band connecting the Dirac points removes the long-time branch-cut-induced tails and instead drives damped Rabi oscillations.

• Bound-state threshold: Real bound states exist outside the band (|ω_BS| > 3J). Because 𝚺(E) is monotonically increasing, exactly one bound state appears only when Δ > J[3 − 0.297(g/J)²] or Δ < −J[3 − 0.297(g/J)²]. The resulting fractional decay leaves the population stabilized at |r_BS|².

• Two unstable poles near Δ = 0: The flat band produces two complex poles shifted into the lower half-plane by the bulk contribution, giving damped Rabi oscillations at a frequency proportional to the effective vacuum Rabi splitting sqrt(Δ² + 4Ω²).

• Population trapping for large detuning: For |Δ| ≫ Ω, one pole approaches zero, preserving population trapping independently of system size.

• Band structure and edge modes: Bulk modes ω±(k,q) = ±sqrt(J² + 𝒥²(k) + 2J𝒥(k)cos q) with 𝒥(k) = 2J cos(k/2) fill a continuous band [−3J, 3J] with Dirac cones at (k,q) = (±2π/3, π). For 2π/3 < |k| ≤ π, zero-energy topological edge modes appear, exponentially localized with length λ_k⁻¹ = ln|J/𝒥(k)|, forming a flat band (ω_FB = 0) strictly on the A sublattice.

• Three dynamical regimes illustrated at g = 0.2J: Δ = 0 (two unstable poles relevant, blue curve), Δ = 1.5J (single unstable pole relevant, red curve), and Δ = 3J (bound state appears, green curve).

• The t⁻⁴ tail is invisible in practice. It is fully masked by Rabi oscillations when Δ < Ω and by fractional decay when Δ > Ω.

• Branch-cut handling: The elliptic integrals are analytically continued beyond the first Riemann sheet to avoid branch cuts at z ∈ {−3J, −J, 0, J, 3J}.

Methodology in Plain English

The authors model the photonic bath as a tight-binding honeycomb lattice where photons hop between nearest-neighbor cavities at rate J. The zigzag edge is placed at n = 0, with open boundaries perpendicular to the edge (n ∈ ℕ) and periodic boundaries along it (m ∈ ℤ). A single qubit of bare frequency Δ couples with strength g to one edge cavity, a₀₀, on the A sublattice, under the rotating-wave approximation.

Rather than simulating the dynamics, they compute the qubit's self-energy exactly — the object that encodes how the bath feeds back on the qubit. This is possible because the lattice Green's function for this geometry can be summed in closed form, yielding elliptic integrals. The qubit amplitude is then obtained from a resolvent integral, c_e(t) = −(1/2πi)∫dE G_e(E + i0⁺)e^{−iEt} with G_e(z) = [z − Δ − 𝚺(z)]⁻¹. Deforming the contour into the lower half-plane splits the result into pole contributions (bound states and unstable poles) plus a branch-cut integral. Expanding the analytically continued Green's function near E = 0 gives the asymptotic long-time tail. The whole analysis is analytical; no numerical diagonalization of the full many-body problem is required.

Why This Matters

Impact on research: The paper establishes lattice boundary conditions as a practical design knob for quantum dynamics. It provides an exact, closed-form benchmark for edge-coupled emitter problems and shows that the same bath can produce qualitatively different physics — Rabi oscillations versus power-law tails — depending purely on whether the emitter sits in the bulk or on the edge.

Potential applications (as motivated by the paper's discussion):

• Circuit QED hardware: Transmon qubits coupled to high-impedance cavity arrays, where the authors note edge-coupling is far easier to implement than bulk-coupling.
• Photonic lattice quantum simulation: Engineered 2D photonic baths that reproduce graphene-like band structures with controlled boundaries.
• Cold-atom quantum simulators: Cited as an alternative platform for realizing graphene-like baths.
• Quantum state storage: Fractional decay and population trapping offer a route to preserving excitation in a chosen emitter, with the trapping behavior shown to be independent of system size for large detuning.

Industry relevance: Superconducting quantum computing and integrated photonics both depend on understanding and controlling decoherence in structured environments. A predictive analytical model for edge-coupled emitters is directly useful for designing cavity-array architectures and for engineering controlled qubit-qubit interactions mediated by a shared flat band.

Future Directions

• Experimental realization: The paper discusses potential platforms only briefly; a concrete demonstration in transmon/high-impedance cQED or photonic lattices remains open.
• Beyond the single-qubit case: How two or more emitters interact through the edge flat band — and whether the effective-cavity picture extends to long-range qubit-qubit coupling — is not addressed.
• Other edge terminations and geometries: Only the zigzag edge is treated; armchair edges and finite ribbons are not analyzed.
• Robustness: The effects of disorder, finite lattice size, and temperature on the bound-state threshold and Rabi frequency are not reported.

Target Audience

Graduate students and researchers in quantum optics, condensed-matter-inspired photonics, and circuit QED who work on non-Markovian open quantum systems, structured photonic baths, or emitter-lattice coupling. Readers looking for experimental data, numerical benchmarks, or device fabrication details will not find them here — the contribution is purely analytical theory.

Authors’ abstract

We present an analytical study of the dynamics of a qubit coupled to the zigzag edge of a photonic honeycomb lattice. Leveraging a resolvent operator approach and exact lattice Green's functions, we map the spectrum of qubit-photon bound states and complex poles to uncover the long-time scaling laws of the qubit population. Unlike a qubit coupled to the bulk of graphene, the presence of a zero-energy flat band connecting the Dirac points suppresses long-time branch-cut-induced tails, instead driving damped Rabi oscillations. Potential experimental platforms are briefly discussed.