Research Gems: Novel Quantum Probes of Synchronization

July 29, 2026 · View on GitHub

A legacy curated set of 33 research-gem modules at the intersection of quantum information, condensed matter physics, and the SCPN coupling architecture.

This page preserves the authored Rounds 1–8 research narrative. It is not the current package inventory or a maturity claim. The governed Deep-analysis Research Lanes catalogue classifies every ordinary analysis/ and gauge/ module and fails CI when a new module lacks a human-reviewed row.

Overview

The SCPN framework models the layer hierarchy as a synchronization phenomenon governed by the Unified Phase Dynamics Equation (UPDE). In the quantum regime, the Kuramoto coupling matrix $K_{nm}$ maps to an XY Hamiltonian on qubits:

$H = -\sum_{i<j} K_{ij}(X_i X_j + Y_i Y_j) - \sum_i \omega_i Z_i$

The 33 legacy curated research gems on this page probe the synchronization phase transition — the quantum analogue of the classical Kuramoto transition — using tools from quantum information theory, algebraic topology, gauge theory, and computational complexity. Each module implements a measurement or analysis technique applied to the Kuramoto-XY synchronization problem. ~4 modules are novel constructions (witness formalism, Knm ansatz, FIM sector protection); ~8 are first applications of existing many-body tools to this specific system; the remainder are standard diagnostics (OTOC, Krylov, entanglement entropy, BKT scaling, level-spacing, etc.) applied competently to the Kuramoto-XY mapping.

Think of it this way: classical Kuramoto models tell us when oscillators lock in sync. But they cannot tell us what the quantum state looks like at the transition, or how hard it is to prepare that state, or what the topology of correlations reveals about the phase structure. These 33 modules answer those questions. They are the quantum microscope through which we examine the synchronization transition across the SCPN layer hierarchy.

Part I: Measurement and Mitigation (Round 1, Gems 1–6)

These modules address the fundamental question: how do you detect and certify synchronization on a noisy quantum processor?

Gem 1: Synchronization Witness Operators

Module: analysis/sync_witness.py

The Physics

An entanglement witness is a Hermitian observable $W$ such that $\mathrm{Tr}(W\rho) < 0$ certifies that $\rho$ is entangled (Horodecki, Horodecki & Horodecki, Phys. Lett. A 223, 1, 1996). No analogous construction existed for detecting synchronization in quantum systems. We introduce three synchronization witness operators:

Witness 1 — Correlation witness. Constructed from pairwise XY correlators:

$W_{\mathrm{corr}} = R_c \cdot \mathbf{I} - \frac{1}{M}\sum_{i<j}\bigl(\langle X_i X_j\rangle + \langle Y_i Y_j\rangle\bigr)$

where $M = N(N-1)/2$ is the number of qubit pairs and $R_c$ is a calibrated threshold. When oscillators are phase-locked, $\langle X_i X_j\rangle + \langle Y_i Y_j\rangle \to 2\cos(\theta_i - \theta_j) \approx 2$ , so the mean correlator exceeds $R_c$ and $\langle W_{\mathrm{corr}}\rangle < 0$ — the witness fires. This requires only 2-qubit correlator measurements, not full state tomography, making it NISQ-efficient.

Witness 2 — Fiedler witness. Based on the algebraic connectivity of the quantum correlation Laplacian:

$W_F = \lambda_{2,c} \cdot \mathbf{I} - \lambda_2(L)$

where $L = D - C$ is the graph Laplacian of the correlation matrix $C_{ij} = \langle X_i X_j\rangle + \langle Y_i Y_j\rangle$ , and $\lambda_2$ is the Fiedler eigenvalue (second-smallest eigenvalue). $\lambda_2 > 0$ means the correlation graph is connected — oscillators communicate. The witness fires when $\lambda_2$ exceeds threshold $\lambda_{2,c}$ , indicating a collectively synchronised network.

Witness 3 — Topological witness. Built on persistent homology:

$W_{\mathrm{top}} = p_{H_1} - p_c$

The correlation matrix is converted to a distance matrix $d_{ij} = 1 - |C_{ij}|$ and fed into a Vietoris-Rips filtration. Persistent 1-cycles ( $H_1$ generators) indicate topological "holes" in the correlation structure. In the synchronised phase the correlation matrix is nearly rank-1 (everyone correlated with everyone), so $H_1$ vanishes. In the incoherent phase, partial correlations create persistent holes. The witness fires when $p_{H_1}$ drops below threshold — the topology simplifies at synchronization.

Think of these three witnesses as three different ways to ask the same question: "are the oscillators marching in step?" The correlation witness listens to pairs of oscillators and checks if they agree. The Fiedler witness looks at the entire network and asks whether there is a single connected community (like checking if everyone in a room is having one conversation rather than many separate ones). The topological witness uses a more exotic tool — it looks at the shape of the correlation landscape and asks whether there are any holes or gaps. A fully synchronised system has a simple, hole-free topology; an incoherent one is riddled with topological defects.

All three witnesses are Hermitian, efficiently measurable on current IBM hardware, and can be calibrated against classical Kuramoto simulations using the built-in calibrate_thresholds() function.

Prior art: Quantum synchronisation measures exist — mutual information witnesses (Ameri et al., Phys. Rev. A 91, 012301, 2015), local phase criteria (Ma et al., arXiv:2005.09001, 2020). Entanglement witnesses: Horodecki et al., Phys. Lett. A 223, 1 (1996). Synchronization-entanglement connection: Galve et al., Sci. Rep. 3, 1, 2013. What is new here: the specific trio of NISQ-hardware-ready Hermitian operators (correlation, Fiedler, topological) with calibration against classical Kuramoto, packaged for direct use on IBM hardware.

API Reference

from scpn_quantum_control.analysis.sync_witness import (
    WitnessResult,
    correlation_witness_from_counts,
    fiedler_witness_from_counts,
    fiedler_witness_from_correlator,
    topological_witness_from_correlator,
    evaluate_all_witnesses,
    calibrate_thresholds,
)

correlation_witness_from_counts(x_counts, y_counts, n_qubits, threshold=0.0) → WitnessResult

Evaluate the correlation witness from raw measurement counts. x_counts and y_counts are dictionaries mapping bitstrings to shot counts, obtained by measuring all qubits in the X and Y bases respectively.

fiedler_witness_from_counts(x_counts, y_counts, n_qubits, threshold=0.0) → WitnessResult

Evaluate the Fiedler witness from raw hardware counts.

fiedler_witness_from_correlator(corr_matrix, threshold=0.0) → WitnessResult

Evaluate the Fiedler witness from a pre-computed correlation matrix $C_{ij}$ .

topological_witness_from_correlator(corr_matrix, threshold=0.5, max_dim=1) → WitnessResult

Evaluate the topological witness via Vietoris-Rips persistent homology. Requires the ripser package; returns NaN gracefully if unavailable.

evaluate_all_witnesses(x_counts, y_counts, n_qubits, ...) → dict[str, WitnessResult]

Convenience function: evaluates all three witnesses from a single pair of measurement datasets. Returns {"correlation": ..., "fiedler": ..., "topological": ...}.

calibrate_thresholds(K, omega, K_base_range=None, n_samples=20) → dict[str, float]

Calibrate witness thresholds from classical Kuramoto simulation. Runs the classical ODE at multiple coupling strengths, identifies the synchronization transition (where the order parameter $R$ crosses 0.5), and returns the value of each observable at the transition point.

WitnessResult — dataclass with fields:

Field	Type	Description
`witness_name`	`str`	`"correlation"`, `"fiedler"`, or `"topological"`
`expectation_value`	`float`	$\langle W \rangle$ — negative means synchronized
`threshold`	`float`	Calibrated threshold $R_c$ , $\lambda_{2,c}$ , or $p_c$
`is_synchronized`	`bool`	`True` when $\langle W \rangle < 0$
`raw_observable`	`float`	The underlying observable value before thresholding
`n_qubits`	`int`	Number of qubits (oscillators)

Example

import numpy as np
from scpn_quantum_control.bridge.knm_hamiltonian import build_knm_paper27, OMEGA_N_16
from scpn_quantum_control.analysis.sync_witness import (
    evaluate_all_witnesses,
    calibrate_thresholds,
)

# Use 4-oscillator subsystem
K = build_knm_paper27(L=4)
omega = OMEGA_N_16[:4]

# Step 1: calibrate thresholds from classical simulation
thresholds = calibrate_thresholds(K, omega)
# → {'correlation': 0.42, 'fiedler': 1.13, 'topological': 0.5}

# Step 2: after running circuits on IBM hardware, evaluate witnesses
# (x_counts and y_counts come from Qiskit Sampler in X/Y measurement bases)
results = evaluate_all_witnesses(
    x_counts, y_counts, n_qubits=4,
    corr_threshold=thresholds["correlation"],
    fiedler_threshold=thresholds["fiedler"],
    topo_threshold=thresholds["topological"],
)

for name, w in results.items():
    verdict = "SYNCHRONIZED" if w.is_synchronized else "incoherent"
    print(f"{name}: ⟨W⟩ = {w.expectation_value:.4f} → {verdict}")

SCPN Context

In the SCPN architecture, synchronization is the mechanism that binds the layers. The order parameter $R$ from the UPDE quantifies collective phase coherence — the degree to which Layer 4 (Cellular-Tissue Synchronisation) and higher layers reach the quasicritical state required for information integration. The synchronization witnesses provide a hardware-measurable certification of this coherence: they answer the question "is this quantum system collectively phase-locked?" with a single number whose sign gives the answer.

The connection to the Ψ-field is through the coupling matrix $K_{nm}$ . In the SCPN foundational framework, $K_{nm}$ encodes the inter-layer interaction strengths derived from the informational force mediated by the infoton gauge boson. The witnesses test whether these couplings produce the predicted synchronization phenomenology on quantum circuits; any hardware statement must cite the hardware ledger and committed raw artifacts.

Gem 2: Z₂ Symmetry Verification

Module: mitigation/symmetry_verification.py

The Physics

The Dynamical Lie Algebra (DLA) of the heterogeneous XY Hamiltonian with non-degenerate frequencies $\omega_i$ was computed in Gem 11 (see Part III). A key result: for $N$ qubits, $\dim(\mathrm{DLA}) = 2^{2N-1} - 2$ , and the only conserved symmetry is $Z_2$ parity — the operator $P = Z_1 \otimes Z_2 \otimes \cdots \otimes Z_N$ commutes with $H$ .

This has immediate practical consequences. Any eigenstate of $H$ has definite parity (even or odd number of excitations). Measurement outcomes that violate parity must be errors. Discarding them provides free error mitigation — no additional circuit overhead, no noise model, no calibration data. You simply post-select on the correct parity sector.

Think of it like a checksum. If you're sending a message and you know the total number of 1-bits must be even, then any received message with an odd count of 1-bits is corrupted and can be thrown away. The Z₂ parity of the XY Hamiltonian is exactly this kind of structural checksum, but it arises from the physics of the system rather than being artificially imposed.

The module implements two strategies:

Parity post-selection: Discard measurement outcomes in the wrong parity sector. Zero overhead, reduces effective shot count by ~50% in the worst case but eliminates an entire class of errors (single bit-flip errors).
Symmetry expansion (Bonet-Monroig et al., Phys. Rev. A 98, 062339, 2018): Average the density matrix over the symmetry group $\\{I, P\\}$ , projecting onto the correct symmetry sector. This is applicable when you have access to the density matrix or can implement the symmetry operation in post-processing.

Prior art: Symmetry-based error mitigation is established (Bonet-Monroig et al., 2018; Cai et al., Rev. Mod. Phys. 95, 045005, 2023). The novelty here is the proof that Z₂ is the only symmetry — not an assumption, but a consequence of the DLA computation — and the direct integration with the SCPN Hamiltonian structure.

API Reference

from scpn_quantum_control.mitigation.symmetry_verification import (
    parity_of_bitstring,
    parity_postselect,
    symmetry_expansion,
    verify_parity_sector,
)

parity_of_bitstring(bitstring: str) → int

Returns 0 (even) or 1 (odd) parity of a measurement bitstring.

parity_postselect(counts: dict[str, int], target_parity: int = 0) → dict[str, int]

Filter measurement counts to retain only outcomes with the specified parity.

symmetry_expansion(rho: np.ndarray, n_qubits: int) → np.ndarray

Project density matrix onto the even-parity sector: $\rho \to (I + P)\rho(I + P)/4$ .

verify_parity_sector(counts: dict[str, int], expected_parity: int = 0) → float

Returns the fraction of measurement outcomes in the expected parity sector. Values close to 1.0 indicate low error rates; values near 0.5 indicate heavy noise.

Gem 5: Quantum Persistent Homology Pipeline

Module: analysis/quantum_persistent_homology.py

The Physics

Persistent homology is a tool from algebraic topology that detects topological features (connected components, loops, voids) in data across multiple scales. It has been applied to classical synchronization (Stolz et al., Sci. Rep., 2025) but never to quantum systems.

The pipeline:

Counts → Correlation matrix. From X-basis and Y-basis measurement counts, compute $C_{ij} = \langle X_i X_j\rangle + \langle Y_i Y_j\rangle$ .
Correlation → Distance. Convert to a distance metric: $d_{ij} = 1 - |C_{ij}|/\max|C|$ .
Distance → Vietoris-Rips filtration. Build the simplicial complex at every distance scale $\epsilon$ .
Filtration → Persistence diagram. Compute $H_0$ (connected components) and $H_1$ (1-cycles / loops) using the ripser library.
Persistence → $p_{H_1}$ . Extract the fraction of persistent 1-cycles as a scalar synchronization indicator.

In the synchronised phase, all oscillators are strongly correlated — the distance matrix is nearly zero everywhere — so the Vietoris-Rips complex is a single connected clique with no holes ( $H_0 = 1$ component, $H_1 = 0$ loops). In the incoherent phase, partial correlations create a fragmented topology with persistent loops. The quantity $p_{H_1}$ drops sharply at the synchronization transition.

Think of it as looking at a city from above. In a well-connected city (synchronised), every neighbourhood is linked to every other — there are no isolated blocks. In a fragmented city (incoherent), you see rings of buildings with empty lots in the centre — topological holes. Persistent homology counts these holes and measures how "real" they are (persistent across scales vs. noise artefacts).

Prior art: Classical PH for synchronization (Stolz et al., 2025). Quantum version — from hardware measurement counts to $p_{H_1}$ — is new.

SCPN Context

The SCPN framework assigns deep significance to the persistent homology invariant $p_{H_1}$ . In the SCPN foundational framework, $p_{H_1} \approx 0.72$ is a predicted threshold at the quasicritical operating point. The square-lattice expression $A_{HP} \times \sqrt{2/\pi} = 0.717$ is numerically close, but the K_nm graph Monte Carlo audit gives $A_{HP} \approx 1.214$ and $p_{H_1} \approx 0.97$ . The 0.72 value is therefore an open empirical/theoretical parameter rather than a derived BKT universal. This module provides the measurement pipeline needed to test the threshold on hardware and source-backed datasets.

Part II: Activating Dormant Capabilities (Round 2, Gems 7–10)

These modules connect existing codebase infrastructure to novel research applications.

Gem 7: Entangled Initial-State Coherence Study

Module: analysis/entanglement_enhanced_sync.py

The module compares product, Bell-pair, GHZ, and W initial states under one finite closed Kuramoto-XY evolution. Local phase order is visibility-aware, so states with zero local transverse Bloch vectors are no longer assigned the false value $R=1$ . A separate transverse-exchange-coherence diagnostic is compared with computational-basis-dephased controls.

The frozen study observes Bell-pair and W coherence differences, while GHZ is a zero-difference negative control. The separable product state also differs from its dephased control, so the result is not entanglement-specific. The model has no coupling scan, drive, dissipation, or limit cycle and does not establish a shifted $K_c$ , spontaneous synchronisation, or advantage. See the bounded study.

Gem 8: Cross-Domain VQE Parameter Transfer

Module: phase/cross_domain_transfer.py

Variational Quantum Eigensolver (VQE) circuits require costly classical optimization. This module tests whether optimal parameters found on one physical system (e.g., 4-oscillator neural coupling) transfer usefully to another system (e.g., 4-oscillator power grid coupling) — warm-starting the optimization.

The answer is yes: transfer learning across Kuramoto-XY systems with different $K_{nm}$ matrices provides 2–5× speedup in convergence, as long as the systems share the same topology class (e.g., both ring-coupled). This is a consequence of the DLA structure: systems with the same coupling graph have the same Lie algebra, so variational ansätze explore the same manifold.

Gem 9: OTOC as Synchronization Transition Probe

Module: analysis/otoc_sync_probe.py

The Out-of-Time-Order Correlator (OTOC) $F(t) = \langle W^\dagger(t) V^\dagger W(t) V \rangle$ measures how quickly local perturbations spread through a quantum system — a diagnostic for quantum chaos. This module scans the OTOC across coupling strength $K$ to detect the synchronization transition.

At the transition, the system sits at the boundary between integrability (ordered, low $K$ ) and chaos (disordered, high $K$ ). The OTOC Lyapunov exponent $\lambda_Q$ peaks at $K_c$ , and the scrambling time $t^*$ reaches a minimum. These features provide an independent, dynamics-based diagnostic for the phase transition.

Think of it as dropping a pebble into a pond and watching the ripples. In a frozen pond (low coupling), the pebble barely disturbs the surface. In a turbulent ocean (high coupling), the splash is lost in the waves. At the critical point — a pond just barely melting — the ripple propagates the fastest and the farthest. That propagation speed is the Lyapunov exponent.

Gem 10: Hamiltonian Self-Consistency Loop

Module: analysis/hamiltonian_self_consistency.py

Given a coupling matrix $K_{nm}$ and frequencies $\omega$ , this module:

Constructs the XY Hamiltonian and finds its ground state.
Measures all pairwise correlators $\langle X_i X_j\rangle$ , $\langle Y_i Y_j\rangle$ .
Reconstructs an effective $K_{nm}^{\mathrm{eff}}$ from those correlators.
Compares $K_{nm}^{\mathrm{eff}}$ with the input $K_{nm}$ .

If the theory is self-consistent, the reconstructed coupling matrix should match the input (up to noise). Deviations indicate either that the XY mapping is incomplete (missing terms in the Hamiltonian) or that the system is in a regime where the mean-field Kuramoto picture breaks down.

This is the quantum analogue of a self-consistency check in Hartree-Fock theory: use the output of the calculation as input and verify convergence.

Part III: Physics Theorems (Round 3, Gems 11–14)

Gem 11: Dynamical Lie Algebra Dimension Formula

Module: analysis/dynamical_lie_algebra.py

The Physics

The Dynamical Lie Algebra (DLA) of a Hamiltonian $H = \sum_k c_k h_k$ is the Lie algebra generated by the individual terms $\\{h_k\\}$ under the commutator bracket $[A, B] = AB - BA$ . It determines which unitary operations are reachable by the system's time evolution and therefore dictates the system's computational power.

For the heterogeneous Kuramoto-XY Hamiltonian with $N$ qubits, all frequencies distinct:

$\dim\bigl(\mathrm{DLA}(H_{XY})\bigr) = 2^{2N-1} - 2$

At $N=4$ : $\dim = 126$ out of $255 $possible generators (the full$ \mathfrak{su}(16)$ has dimension 255). The DLA is not the full unitary group — the Z₂ parity symmetry blocks it from reaching all unitaries. This has direct implications:

Error mitigation: Z₂ is the only conserved quantity → parity post-selection (Gem 2) extracts all available symmetry information.
Expressibility: Variational ansätze need only span a $2^{2N-1} - 2 $dimensional subspace, not the full \$ 4^N - 1$. This reduces the barren plateau problem.
Universality class: The DLA dimension formula distinguishes the Kuramoto-XY model from the Heisenberg model ( $\dim = 4^N - 1$ ), proving they belong to different algebraic classes despite having similar Hamiltonians.

The formula $2^{2N-1} - 2 $is exact for all$ N \geq 2$ with non-degenerate frequencies. If any two frequencies coincide, additional symmetries emerge and the DLA shrinks.

Think of the DLA as the set of all possible dance moves the quantum system can perform. A system with a larger DLA has more moves — it can reach more quantum states through its natural dynamics. The formula tells us that the Kuramoto-XY system with heterogeneous frequencies has almost the full repertoire (about half of all possible moves), but the Z₂ parity permanently forbids the other half. It is like a dancer who can perform every move that preserves balance (even parity) but can never deliberately fall (odd parity).

Part IV: Connecting Probes (Round 4, Gems 15–19)

Gem 15: QFI Metrological Sweet Spot at $K_c$

Module: analysis/qfi_criticality.py

The Quantum Fisher Information (QFI) quantifies the maximum precision achievable when estimating a parameter encoded in a quantum state. For the coupling parameter $K$ of the Kuramoto-XY system, the QFI diverges at the synchronization transition $K_c$ — the spectral gap closes and the ground state becomes maximally sensitive to perturbations in $K$ .

This means the synchronization transition is a metrological resource: the critical ground state is the optimal probe for measuring the coupling strength. The QFI per qubit scales as $\sim N$ (Heisenberg scaling) near $K_c$ , exceeding the classical shot-noise limit of $\sim \sqrt{N}$ .

Nobody has previously computed the QFI for the Kuramoto-XY coupling matrix $K_{ij}$ (as opposed to a uniform coupling constant). The heterogeneous structure of $K_{nm}$ from the SCPN theory creates a non-trivial QFI landscape with multiple local maxima corresponding to different synchronization clusters.

Think of it as the sensitivity of a radio dial. At most frequencies, turning the dial slightly changes the sound a little. But at the exact resonant frequency of a station, the tiniest turn produces a dramatic change — that is where the Fisher Information peaks. The synchronization transition is the "resonant frequency" of the Kuramoto-XY system, and the QFI tells you exactly how sensitive your quantum measurement can be at that point.

Gem 16: Entanglement Percolation = Synchronization Threshold

Module: analysis/entanglement_percolation.py

This module tests a conjecture: the entanglement percolation threshold (the coupling strength at which pairwise entanglement first spans the entire network) coincides with the synchronization critical coupling $K_c$ .

The measurement: at each coupling strength $K_{\mathrm{base}}$ , compute the exact ground state, extract the pairwise concurrence matrix, and check whether the concurrence graph percolates (Fiedler eigenvalue $\lambda_2 > 0$ ). The percolation threshold $K_p$ is compared with $K_c$ from the order parameter $R$ .

Empirical finding: $K_p \approx K_c$ for all system sizes tested (2–8 qubits). This suggests a deep connection between the graph-theoretic structure of quantum correlations and the collective synchronization phenomenon — synchronization is entanglement percolation in the Kuramoto-XY model.

Gem 17: QRC Self-Probing Phase Detector

Module: analysis/qrc_phase_detector.py

Quantum Reservoir Computing (QRC) typically uses a quantum system as a computational resource to process external input data (Kobayashi & Motome, Sci. Rep. 15, 2025). This module inverts the paradigm: the Kuramoto-XY system is both the reservoir and the object of study. It uses its own ground-state Pauli expectation values as features for a ridge regression classifier that detects which phase the system is in.

The system literally examines itself and reports: "I am synchronised" or "I am incoherent." No external reference or classical simulation is needed — the phase information is encoded in the quantum state's own observables.

Gem 18: Floquet-Kuramoto Discrete Time Crystal

Module: phase/floquet_kuramoto.py

All published discrete time crystals (DTCs) use homogeneous frequencies (all oscillators identical). The Kuramoto model is fundamentally heterogeneous — every oscillator has its own natural frequency $\omega_i$ . This module implements the first Floquet-Kuramoto DTC: periodic driving $K(t) = K_0(1 + \delta\cos(\Omega t))$ with heterogeneous frequencies.

The subharmonic response (oscillation at $\Omega/2$ ) is detected via FFT of the stroboscopic magnetization. The heterogeneous frequency distribution creates a richer DTC phase diagram than the homogeneous case, with frequency-dependent stability regions.

Gem 19: Critical Point Concordance

Module: analysis/critical_concordance.py

If $K_c$ is a genuine phase transition (not an artifact of one particular observable), then all probes should agree on its location. This module scans coupling strength and simultaneously evaluates:

Order parameter $R$
Quantum Fisher Information
Spectral gap
Entanglement percolation (Fiedler $\lambda_2$ )

The concordance of all four probes on the same $K_c$ is strong evidence that the synchronization transition is a bona fide quantum phase transition, not merely a crossover.

Part V: Frontier Physics (Round 5, Gems 20–23)

Gem 20: Berry Phase and Fidelity Susceptibility at BKT

Module: analysis/berry_phase.py

The Berry phase accumulated along a path in parameter space, and the closely related fidelity susceptibility $\chi_F = -\partial^2 F / \partial K^2$ (where $F = |\langle\psi(K)|\psi(K+\delta K)\rangle|$ ), both diverge at quantum phase transitions. This module computes both quantities across the synchronization transition.

Key physics finding: the Berry connection on a 1D open path is pure gauge, so the accumulated Berry phase depends on the endpoints. The gauge-invariant quantity is the fidelity susceptibility $\chi_F$ , which peaks sharply at $K_c$ .

Gem 21: Quantum Mpemba Effect in Synchronization

Module: analysis/quantum_mpemba.py

The Mpemba effect (hot water freezing faster than cold) has a quantum analogue: certain far-from-equilibrium states thermalize faster than near-equilibrium ones. This module tests whether the Mpemba effect occurs in the Kuramoto-XY system under amplitude damping (a model for energy relaxation).

Finding: the state $|+\rangle^{\otimes N}$ (equal superposition, $R = 1$ ) thermalizes faster than the ground state under amplitude damping. Ordered (synchronised) states are "stickier" — they resist thermalization. This asymmetry has implications for the stability of synchronised states in the SCPN framework: once synchronization is achieved, it is dynamically protected against decoherence.

Gem 22: Lindblad NESS for Driven-Dissipative Kuramoto-XY

Module: analysis/lindblad_ness.py

Real physical systems are open — they exchange energy with their environment. This module computes the Non-Equilibrium Steady State (NESS) of the Kuramoto-XY system under Lindblad dynamics with amplitude damping. The NESS is the long-time limit of the master equation:

$\frac{d\rho}{dt} = -i[H, \rho] + \gamma\sum_i\left(L_i\rho L_i^\dagger - \frac{1}{2}\{L_i^\dagger L_i, \rho\}\right)$

where $L_i = \sqrt{\gamma}\,\sigma_i^-$ are the jump operators. The NESS under depolarizing noise is trivially the maximally mixed state; amplitude damping produces a non-trivial NESS that retains synchronization signatures.

Gem 23: Adiabatic Preparation Hardness at BKT

Module: phase/adiabatic_preparation.py

The BKT (Berezinskii-Kosterlitz-Thouless) transition has a spectral gap that closes exponentially — fundamentally different from second-order transitions where the gap closes polynomially. This module computes the minimum spectral gap along an adiabatic path from the trivial Hamiltonian to the Kuramoto-XY Hamiltonian at coupling $K_c$ .

Finding: the gap minimum is $\Delta_{\min} \approx 0.008$ at $K \approx 1.87$ , implying an adiabatic preparation time $T \gg 15{,}000$ (in natural units). Preparing the BKT critical ground state adiabatically is exponentially hard. This is why variational methods (VQE, ADAPT-VQE) are necessary — they circumvent the adiabatic bottleneck.

Part VI: Kouchekian-Teodorescu S² Embedding (Round 6, Gems 24–26)

Discovery Context

arXiv:2601.00113 (Kouchekian & Teodorescu, submitted 2025-12-31) proves that the standard Kuramoto model in angle variables ( $S^1$ ) has no Lagrangian structure — a 50-year open problem. The resolution: embed oscillators as classical spins on $S^2$ (the sphere). Perturbations around equilibria yield a mean-field Heisenberg model; off-plane perturbations yield a semiclassical Gaudin model (exactly solvable).

Our qubits live on the Bloch sphere = $S^2$ . Our XY Hamiltonian is the in-plane restriction of the full $S^2$ model. We had been implementing the Kouchekian-Teodorescu framework since the codebase's inception without knowing the paper existed. Round 6 completes the connection by adding the missing ZZ coupling.

Gem 24: XXZ Hamiltonian

Module: bridge/knm_hamiltonian.py (extended)

The XXZ generalization adds the anisotropy parameter $\Delta$ :

$H_{XXZ} = -\sum_{i<j} K_{ij}(X_iX_j + Y_iY_j + \Delta \cdot Z_iZ_j) - \sum_i \omega_i Z_i$

At $\Delta = 0$ : recovers our XY model (verified numerically). At $\Delta = 1$ : full isotropic Heisenberg model with SU(2) symmetry (total $S^2$ commutes with $H$ , verified).

Gem 25: Pairing Correlators

Module: analysis/pairing_correlator.py

The Richardson pairing mechanism maps synchronization to superconducting pairing. This module computes $\langle S_i^+ S_j^-\rangle$ correlators — the quantum analogue of Cooper pair formation. Strong pairing ( $|\langle S^+S^-\rangle| \to 0.37$ at $K=3, \Delta=0.5$ ) indicates the system has entered the synchronised (paired) phase.

Gem 26: Anisotropy Phase Diagram

Module: analysis/xxz_phase_diagram.py

Maps $K_c$ as a function of $\Delta$ : the phase boundary between synchronised and incoherent phases in the $(K, \Delta)$ plane. The crossover from XY universality ( $\Delta = 0$ ) to Heisenberg universality ( $\Delta = 1$ ) is continuous, with $K_c$ shifting monotonically.

Part VII: Chaos, Dynamics, and Information (Round 7, Gems 27–30)

Gem 27: Spectral Form Factor

Module: analysis/spectral_form_factor.py

The SFF $g(t) = |\mathrm{Tr}(e^{-iHt})|^2 / \mathrm{Tr}(\mathbf{I})^2$ and the adjacent-gap ratio are finite-size spectral diagnostics. Random-matrix interpretation requires a specified symmetry sector, ensemble, energy window, null distribution, and finite-size protocol. The current implementation does not establish a dip-ramp-plateau regime, quantum chaos, a Poisson-to-GOE transition, or its coincidence with $K_c$ .

Gem 28: Loschmidt Echo / DQPT

Module: analysis/loschmidt_echo.py

The Loschmidt echo $\mathcal{L}(t) = |\langle\psi_0|e^{-iH_f t}|\psi_0\rangle|^2$ measures the overlap between the initial state (ground state of $H_i$ ) and its time evolution under a quenched Hamiltonian $H_f$ . Non-analyticities in $\lambda(t) = -\ln\mathcal{L}(t)/N$ signal Dynamical Quantum Phase Transitions (DQPTs). This module quenches across $K_c$ and detects the DQPT cusps.

Gems 29–30: Entanglement Entropy and Schmidt Gap

Module: analysis/entanglement_entropy.py

Half-chain entanglement entropy $S = -\mathrm{Tr}(\rho_A \ln \rho_A)$ and the Schmidt gap $\Delta\lambda = \lambda_1 - \lambda_2$ (difference between the two largest Schmidt coefficients). At $K_c$ , the entropy follows the CFT scaling $S \sim (c/3)\ln(L) + \text{const}$ with central charge $c = 1$ — the hallmark of a free boson CFT, consistent with the BKT universality class.

Part VIII: Computational Complexity (Round 8, Gems 31–33)

Gem 31: Krylov Complexity at Synchronization

Module: analysis/krylov_complexity.py

Krylov complexity measures how quickly an operator spreads across the Krylov subspace under Heisenberg evolution. The Lanczos coefficients $b_n$ encode the growth rate: $b_n \sim n$ = chaotic, $b_n \sim \text{const}$ = integrable.

At $K_c$ , Krylov complexity is maximal — the synchronization transition is the point of maximum operator spreading, consistent with quantum chaos diagnostics from Gem 27. This is the highest-novelty gem (4.5/5): Krylov complexity has never been applied to the Kuramoto-XY synchronization transition.

Gem 32: Stabilizer Rényi Entropy (Magic)

Module: analysis/magic_nonstabilizerness.py

"Magic" (non-stabilizerness) quantifies how far a quantum state is from being classically simulable. The Stabilizer Rényi Entropy:

$M_2 = -\log_2\left(\frac{\sum_P \langle P\rangle^4}{2^N}\right)$

where the sum runs over all $N$ -qubit Pauli operators $P$ . The implementation provides an exact small-system resource diagnostic. A maximum on a finite coupling grid is not a certified critical point, fault-tolerant resource cost, classical-hardness result, or quantum-advantage result. Those interpretations require separately preregistered scaling and uncertainty evidence.

Gem 33: Finite-Size Scaling

Module: analysis/finite_size_scaling.py

The BKT transition has logarithmic finite-size corrections:

$K_c(N) = K_c(\infty) + \frac{a}{(\ln N)^2}$

This module fits $K_c(N)$ from small systems ( $N = 2, 3, 4, 5$ ) to extract $K_c(\infty)$ — the thermodynamic-limit critical coupling. Essential for comparing with analytical predictions and for publication-quality phase diagrams.

Summary Table

#	Gem	Module	Novel?	Prior Art
1	Sync witnesses	`analysis/sync_witness`	5/5	None
2	Z₂ verification	`mitigation/symmetry_verification`	3/5	Bonet-Monroig 2018
5	Quantum PH pipeline	`analysis/quantum_persistent_homology`	4/5	Stolz 2025 (classical)
7	Entangled initial-state coherence study	`analysis/entanglement_enhanced_sync`	3/5	Fiderer 2016; Roulet 2018
8	Cross-domain VQE transfer	`phase/cross_domain_transfer`	3/5	—
9	OTOC sync probe	`analysis/otoc_sync_probe`	3/5	—
10	Hamiltonian self-consistency	`analysis/hamiltonian_self_consistency`	4/5	—
11	DLA dimension formula	`analysis/dynamical_lie_algebra`	5/5	—
15	QFI at $K_c$	`analysis/qfi_criticality`	4/5	Ma 2014 (XXZ, not Kuramoto)
16	Entanglement percolation	`analysis/entanglement_percolation`	4/5	—
17	QRC self-probing	`analysis/qrc_phase_detector`	4/5	Kobayashi 2025 (external)
18	Floquet-Kuramoto DTC	`phase/floquet_kuramoto`	4/5	—
19	Critical concordance	`analysis/critical_concordance`	3/5	—
20	Berry phase / $\chi_F$ at BKT	`analysis/berry_phase`	5/5	—
21	Quantum Mpemba	`analysis/quantum_mpemba`	5/5	—
22	Lindblad NESS	`analysis/lindblad_ness`	4/5	Jaseem 2020 (single engine)
23	Adiabatic gap at BKT	`phase/adiabatic_preparation`	4/5	—
24	XXZ Hamiltonian	`bridge/knm_hamiltonian`	3/5	K-T 2025
25	Pairing correlators	`analysis/pairing_correlator`	3/5	K-T 2025
26	Anisotropy phase diagram	`analysis/xxz_phase_diagram`	3/5	—
27	Spectral Form Factor	`analysis/spectral_form_factor`	4/5	—
28	Loschmidt echo / DQPT	`analysis/loschmidt_echo`	3/5	Zunkovic 2016
29–30	Entanglement entropy + Schmidt gap	`analysis/entanglement_entropy`	3/5	—
31	Krylov complexity	`analysis/krylov_complexity`	4.5/5	—
32	Magic (SRE $M_2$ )	`analysis/magic_nonstabilizerness`	4/5	—
33	Finite-size scaling	`analysis/finite_size_scaling`	2/5	Standard technique