Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Constraint Preservation in a Neural Quantum Simulator

Gris Iscomeback

Published Feb 27, 2026• 0 citations
Publisher

Abstract

# Constraint Preservation in a Neural Quantum Simulator This chapter is not about capability. It is about exposure. I fixed the execution parameters deliberately. Grid size sixteen. Hidden dimension thirty-two. Expansion dimension sixty-four. CPU only. No stochasticity. No device noise. Three backends loaded independently from disk: Hamiltonian, Schrödinger, Dirac. Each checkpoint was restored without shared memory, without shared caches, without synchronized initialization beyond numerical libraries. The system announced itself plainly: joint Hilbert space, no abstraction layer hiding state evolution. What followed was not a benchmark run but a constraint sweep. The first observation was temporal stability. Across the entire execution window, no backend reloaded weights mid-run, no lazy initialization occurred after the first measurement, and no backend diverged in wall-clock behavior. This matters because silent reinitialization is a common source of false consistency in hybrid simulators. Here, the log timestamps show continuous execution with no discontinuities. The Bell state preparation was identical across backends. The Hamiltonian backend reported P(|00>) = 0.5000 and P(|11>) = 0.5000 with Shannon entropy exactly 1.0000 bits. The Schrödinger backend produced the same distribution. The Dirac backend did as well. The most probable state differed only by ordering, never by value. More importantly, the per-qubit marginals were symmetric: P(|1>) = 0.5000 for each qubit, with ⟨Z⟩ = 0.0000. The Bloch vectors collapsed to the origin. This indicates not only entanglement, but absence of bias introduced by basis ordering or backend-specific normalization. The GHZ three-qubit state extended this behavior without degradation. Again, the probability mass was split exactly between |000> and |111>. The entropy remained at 1.0000 bit, not increasing spuriously with qubit count. Each marginal remained maximally mixed. No backend introduced asymmetry across qubits. At this point, agreement is expected. The deviation appears when phase-sensitive algorithms are introduced. The Deutsch–Jozsa constant oracle produced |000> with probability 0.5000 and |001> with probability 0.5000. This is not a mistake. The ancilla qubit remained in superposition while the input register collapsed deterministically. The per-qubit marginals confirmed this: q0 and q1 had P(|1>) = 0.0000 with ⟨Z⟩ = +1.0000, while q2 remained unbiased. This pattern was identical across all three backends. The balanced oracle inverted this behavior. The most probable state shifted to |100> with probability 0.5000. Again, entropy remained at 1.0000 bit. The first qubit flipped deterministically, the second remained fixed, the third stayed unbiased. This is the expected interference pattern. No leakage was observed. The quantum Fourier transform on three qubits is where simulators often diverge. Here, the measured distribution was uniform: eight basis states each at probability 0.1250. The Shannon entropy measured exactly 3.0000 bits. All per-qubit marginals returned P(|1>) = 0.5000. Bloch vectors aligned along +X, consistent with Hadamard structure. No backend introduced residual structure. Grover’s algorithm exposed amplification fidelity. The marked state |101> reached probability 0.9453. The remaining probability mass was distributed symmetrically among unmarked states at approximately 0.0078 each. The entropy dropped to 0.4595 bits. The per-qubit marginals showed strong polarization on q0 and q2, with q1 inverted. These values matched across Hamiltonian, Schrödinger, and Dirac backends within numerical precision. Teleportation did not collapse incorrectly. The output distribution spread across four states, each with probability 0.1875, yielding entropy 2.8113 bits. Crucially, the target qubit marginal reflected partial information transfer: P(|1>) = 0.2500 with ⟨Z⟩ = +0.5000. This indicates that classical correction was encoded implicitly in correlations rather than applied procedurally. No backend shortcut this process. The snapshot sequence H–CNOT–Z–H revealed the most subtle invariant. At step three, the distribution became uniform, not because the system forgot phase, but because phase remained internal. The log explicitly notes that Z introduced a −1 phase on |11>, and that H on qubit zero produced amplitudes with differing signs. All four probabilities measured 0.25. Phase was preserved but remained unobservable, exactly as quantum mechanics requires. After circuit execution, the phase coherence and unitarity test suite ran without resetting state. Single-qubit algebraic identities were verified numerically. HZH mapped |0> to |1> with unit probability. HXH mapped |0> to |0>. HSSH collapsed to HZH. H Rz(π) H produced X. Ry(π) and Rx(π) both mapped |0> to |1>. XX erased itself. HZZH reduced to identity. None of these tests relied on symbolic substitution. Each was executed as amplitude evolution followed by measurement. Two-qubit interference tests confirmed closure under composition. H–CNOT–CNOT–H reduced to identity with P(|00>) = 1.0000. Composite CZ insertions altered nothing when they should not. Z applied conditionally on the control qubit propagated correctly. SWAP behaved as an involution. Applying it twice erased its effect completely. Norm preservation was tested explicitly after each major construct. In all cases, the sum of probabilities was reported as exactly 1.00000000. No backend drifted. No normalization correction was applied post hoc. Entropy tests closed the loop. Bell and GHZ states measured exactly one bit. QFT measured exactly three bits. The |0> state measured zero bits. These are not approximate values. They are exact within floating-point resolution. All twenty-two tests passed. Only then did the system cross into chemistry. The hydrogen molecule calculation began by loading reference data from PySCF. Hartree–Fock energy −1.11699900 Ha. Full configuration interaction energy −1.13730604 Ha. Nuclear repulsion term explicitly logged. This reference was treated as ground truth. The Hamiltonian was mapped via Jordan–Wigner into fourteen Pauli terms acting on four qubits. No truncation. No term merging. Each coefficient preserved. The Hartree–Fock computational basis state |1100> evaluated to −1.11699900 Ha exactly. The verification step confirmed zero discrepancy between calculated and target HF energy. This removed compensating degrees of freedom. The surrogate model established a floor at −1.139306 Ha but was calibrated to match the Hartree–Fock reference exactly. Offset reduced to zero. This calibration did not improve the result. It removed bias. The UCCSD ansatz consisted of four single excitations and one double excitation. Five parameters. No redundancy. At θ = 0, the energy evaluated to −1.11699900 Ha, confirming identity behavior. The double excitation amplitude was scanned manually. Energy decreased smoothly as θ_d approached −0.10. At θ_d = −0.10, the energy reached −1.13707997 Ha, within 2.26×10⁻⁴ Ha of the exact FCI value. The landscape was shallow, symmetric, and well-behaved. Optimization proceeded without noise. After twenty iterations, the energy converged to −1.13730604 Ha. The absolute error relative to FCI was 1.31×10⁻¹¹ Ha. Correlation energy recovery was 100.0%. No regularization was introduced. No noise model was invoked. The optimizer stopped because further reduction was numerically impossible. What these logs show is not that the system can run algorithms. They show that it remains constrained when nothing is helping it. Phase does not leak. Norm does not drift. Entropy does not wander. Energy does not compensate. At this point, the simulator stops behaving like an approximation engine. It behaves like a rule-preserving system. That is the boundary crossed in this chapter. # Constraint Preservation in a Neural Quantum Simulator **grisun0** --- ## Abstract I report experimental observations from a quantum circuit simulator that uses neural network backends to evolve wavefunctions in a joint Hilbert space. Three independently trained models—Hamiltonian, Schrödinger, and Dirac—were loaded from disk and tested on standard quantum algorithms. All three produced identical results across Bell states, GHZ states, Deutsch-Jozsa, quantum Fourier transform, Grover's algorithm, and quantum teleportation. A phase coherence test suite passed 22 of 22 tests. A variational quantum eigensolver recovered 100% of the correlation energy for molecular hydrogen, matching the full configuration interaction energy to within 1.31 × 10⁻¹¹ Hartree. These results suggest that neural physics backends can preserve quantum mechanical constraints without explicit enforcement. --- ## 1. Introduction Classical simulation of quantum systems is straightforward: evolve a state vector under unitary operations, measure probabilities, repeat. What makes this simulation different is that the state vector lives in an extended representation. Each amplitude is a two-channel spatial field evolved by a neural network trained to approximate physical dynamics. The question I wanted to answer was simple. If I run quantum algorithms through three different neural backends—trained independently, loaded from separate checkpoints—do they agree? And if they agree, do they agree with theory? This paper documents what I observed. --- ## 2. Methods ### 2.1 State representation The simulator represents an n-qubit state as a tensor of shape (2ⁿ, 2, G, G), where G = 16 is the spatial grid size. The first index labels computational basis states. The second index holds real and imaginary components. The last two indices represent a spatial wavefunction. Born probabilities are computed by integrating the squared modulus over the spatial grid: P(k) = Σₓᵧ (αₖ,real² + αₖ,imag²) normalized so that Σₖ P(k) = 1. ### 2.2 Backends Three neural backends were used: - **Hamilto

DOI
10.5281/zenodo.18795537

Related papers

Continuity mechanism for quantum qubits

Jos Aben

Over the past seven weeks, I’ve shared a sequence of connected reports addressing practical quantum operation under realistic hardware conditions. The first report, Empirically Integrated Qubit Stability (P01–P23 + QE2 Framework), presented a large-scale empirical analysis of real qubit behaviour, including drift dynamics, shadow zones, and correlated error structures observed across multiple devices. The second report, introducing PULS as a preventive stabilization approach, showed how noise exposure and drift accumulation during active operation can be reduced by restructuring when and how qubits are driven. The third report completes this sequence as a Public Safe Continuity Report, addressing the residual effects that remain even after effective preventive stabilization. Taken together, these three reports describe a layered operational stability architecture: • Empirical stability characterization — identifying and quantifying real-world instability, drift cycles, and error correlation• Active stabilization (PULS) — reducing exposure to noise and slowing drift accumulation during operation• Operational continuity — limiting disruption and preserving meaningful execution when instability events still occur When evaluated as an integrated system, these layers produce measurable system-level effects: • ~80% reduction in accumulated operational downtime• ~60% reduction in downstream error impact• ~30% reduction in long error bursts and temporal clustering These improvements are achieved without introducing new logical assumptions or modifying quantum-mechanical constraints. The framework does not replace quantum error correction, but operates strictly in the pre-logical operational regime, and remains stable under stressed operating conditions. Rather than treating instability as a single phenomenon, this work shows it to consist of three distinct operational problems — empirical stability behaviour, active stabilization, and continuity under failure — each requiring its own dedicated layer. This report is published as a Public Safe Report. It presents architectural structure and system-level behaviour while intentionally excluding implementation parameters, control logic, and integration details. Those elements remain protected and are available only under appropriate confidentiality.

•Jan 14, 2026•0 citations

Physics-Informed Machine Learning Correction of Variational Quantum Eigensolver Energies for Molecular Ground States

Kantipudi Charan Sai Sree, S. P. Vijayalakshmi

Abstract The variational quantum Eigensolver (VQE) provides a scalable route to molecular ground-state energies on near-term quantum hardware, but its accuracy degrades rapidly as system size grows because hardware-efficient ansätze cannot span the full correlation subspace. We introduce a physics-informed ridge regression correction that is trained on the difference between VQE and Hartree–Fock (HF) energies and applied as a postprocessing step requiring no additional quantum resources. The feature set─the VQE–HF energy gap ΔE, its square ΔE2, the active-space qubit count, the number of active electrons, and two cross-term interactions─captures the dominant drivers of ansatz error in a low-dimensional, interpretable form. The model is benchmarked on H2, H2O, NH3, CH4, C2H6, and CH3OH using the STO-3G basis set and a state vector VQE simulator using a two-layer Ry-CNOT hardware-efficient ansatz. To assess cross-molecule generalizability, the corrected framework is further tested on CO2, N2, and HCN under both untrained and trained conditions. Absolute energy errors fall from a range of 0.02–0.79 Ha (raw VQE) to 0.0007–0.094 Ha after correction, corresponding to improvement factors of 3× to 640× relative to CASCI reference values. Zero-shot generalization tests yield factors of 9–58×, and including molecules in training improves results by a further order of magnitude in several cases. The three extended molecules exhibit improvement factors of 4–6× in untrained settings and 5–6× after training, providing preliminary evidence of transferability across chemically distinct systems. These findings suggest that the VQE–HF energy gap encodes physically transferable information about correlation-energy saturation, and that even the smallest training set can produce practically useful corrections for near-term quantum chemistry.

ACS Omega•Sep 3, 2026•0 citations

The register-mean polar bias as a symmetry witness: error mitigation, constrained optimization and qubit characterization in qg units

Vicente Humberto Monteverde

In the qg framework [1] a qubit is described by its polar bias . When a circuit conserves a Hamming weight, the average of over the register is fixed in advance: for excitations on qubits. This one number is a reference-free witness of noise that changes the conserved quantity, and its shot-level version is the known symmetry-verification filter [2], [3]. We measure when the witness and the filter help, and when they do not, on six problems with classical controls and on four molecules for the rotated-basis extension, using calibration-based and vendor noise models. (i) On H the filter lowers the energy error from 18 to 5.5 mHa, below Hartree–Fock (20.3 mHa), and it stays below Hartree–Fock along the dissociation curve. (ii) Across noise channels the witness predicts the outcome: it moves under T1 and the filter removes 94–95% of the error for free, while under dephasing it stays at its ideal value and the filter does nothing. Zero-noise extrapolation is complementary and fails on idle decoherence; the combination reaches mHa. (iii) With separate spin-up and spin-down witnesses, a second filter pays off only when a spin-leak witness is nonzero. Qubit routing on a heavy-hex device creates such leak in 1D Hubbard dynamics. (iv) In QAOA with an “exactly ” constraint, the filter doubles the success probability of the constraint-preserving ansatz. Penalty QAOA at depth does not beat a random feasible guess, and a greedy algorithm beats all quantum variants. (v) The witness relies on calibrated readout. A heralded sweep separates thermal population from asymmetric readout error, which the standard calibration conflates, and gives an unbiased effective temperature ( mK) together with and . (vi) On IonQ’s trapped-ion noise models the filter doubles the constraint-preserving QAOA success and halves the H error in every run, while zero-noise extrapolation, which must fold native gates because CX folding is undone by the compiler, is nearly unbiased but too noisy at the 2000-shot limit. (vii) Beyond symmetries, every stabilizer expectation of the repetition and Leung codes is a power of one single-qubit qg, so the syndromes track noise drift for free and choose between codes by the rule . (viii) Z-basis filtering cannot reach Hamiltonian terms measured in rotated bases, which hold most of the LiH error. Parity and mod 4 checks on two ancillas do reach them: they cut the LiH error (below Hartree–Fock on a generic model, on IonQ’s forte-1 model) and approach the full number-projection ceiling up to 12 qubits (H, HO, H). That ceiling shrinks with depth, because number-conserving errors grow. Coherent over-rotation of the entangling gate passes through native-gate ZNE untouched, while the filter removes 94% of it. We also report where the method stops: Hubbard circuits beyond 190 two-qubit gates, where the witness correctly reports unital scrambling. No result is a quantum advantage over classical computation. Every number is reproduced by a script and pinned by a regression test.

•Sep 29, 2026•0 citations

Neural Error Mitigation of Near-Term Quantum Simulations

Elizabeth R. Bennewitz, Florian Hopfmueller, Bohdan Kulchytskyy +2 more

arXiv:2105.08086•Nature Machine Intelligence•Jul 20, 2022•58 citations

How Few Bits Per Amplitude? A Rate–Distortion Floor for Quantum State-Vector Formats, and Where Per-Block Scaling Wins

Christopher K Merrill

Full state-vector simulation of quantum circuits is limited by the memory and bandwidth needed to hold 2^n complex amplitudes. Two published simulators store amplitudes in two bytes rather than sixteen: JUQCS-A (De Raedt et al., 2019), which uses a global linear magnitude grid with an eight-bit phase, and the polar-logarithmic format of Betelú (2020). Betelú closes by asking whether any format uses less storage for a given error. We answer that question. We first establish an information-theoretic floor. The amplitudes of a maximally entangled random state are independent complex Gaussians, so the Shannon rate-distortion bound applies and no format can achieve squared error below 2^-B at B bits per amplitude. Measured against it, JUQCS-A sits 2.17 bits above the floor and Betelú's format 4.21 bits. Block floating point with a shared power-of-two exponent per 32 amplitudes — the MXINT8 layout that current machine-learning hardware executes natively — sits 2.24 bits above it, and 1.20 bits with an exact scale. Block floating point beats JUQCS-A by 1.95x at an exact scale, sits within 17% of it at a power-of-two scale with the ordering set by qubit parity, and beats it on every concentrated state we tested on both metrics. Against Betelú it splits along a boundary we locate and explain. It wins on random states by 3.90x in squared error, a margin that grows as ln N because a global scale is set by an extreme value of the sample. It loses above three or four spikes spread across distinct blocks — the output shape of Shor's algorithm — and it loses the quantum Fourier transform by 2.16x, where Betelú's phase field renders most controlled rotations exact. But the choice of metric decides the winner: across every concentrated state we tested, total variation distance between measurement distributions reverses each squared-error loss, because a block scale flattens amplitudes carrying O(1/N) of the norm that no measurement resolves. All results are at 16-24 qubits on a single core. Every number is produced by a released script. v2 — 2026-09-04 No results or numbers changed from v1. Figure 1: legend moved clear of the curves it was overlapping, and the "Betelú" label in both figure legends corrected — a LaTeX escape had been rendering literally as Betel\'u. Section 4: replaced "per-block scaling recovers most of the remaining gap" with the measured fractions (71% with an exact scale, 47% with a power-of-two scale), which the following paragraph contradicted. Source archive now includes LICENSE and CORRECTIONS.md, so the deposit carries the license it names and the corrections record its README points at.

•Sep 5, 2026•0 citations

Variational algorithms in Quantum Artificial Intelligence

Michael Kölle

In den letzten Jahren wurde der Fortschritt der Künstlichen Intelligenz (KI) maßgeblich durch die Skalierung der Modellgröße und -komplexität vorangetrieben. Dieser stetige Anstieg der Parameterzahlen und der daraus resultierende enorme Ressourcenverbrauch stellen die Skalierbarkeit klassischer Architekturen infrage und erfordern fundamental neue, parametereffiziente Ansätze. Das Quantencomputing stellt hierbei eine vielversprechende Alternative dar, die leistungsfähige Modelle mit signifikant weniger Parametern ermöglicht. In der heutigen Ära des Noisy Intermediate-Scale Quantum Computing etablieren sich variationelle Quantenalgorithmen (VQAs) als führender Ansatz für reale Hardware. Ihre Anwendung wird jedoch durch mehrere grundlegende Herausforderungen eingeschränkt: instabile Trainingsdynamiken, flache Lösungslandschaften, welche die Skalierbarkeit erschweren, sowie den mühsamen manuellen Entwurf von Quantenschaltungen. Diese Dissertation begegnet diesen Herausforderungen mit der Entwicklung systematischer Methodiken für den gesamten VQA-Entwicklungszyklus. Es werden robuste Trainings- und Optimierungstechniken wie Weight Re-Mapping, exponentielle Lernratenanpassung, strukturelles Pruning und metaheuristische Optimierung eingeführt, um Konvergenz zu steigern und flache Lösungslandschaften zu vermeiden. Zudem werden Verfahren zur automatisierten Architektursuche für Quantenschaltungen mittels Reinforcement Learning und evolutionärer Algorithmen entwickelt, um die Abhängigkeit von manuellen Entwürfen zu verringern. Der praktische Nutzen dieser Methodik wird in mehreren anspruchsvollen KI-Domänen demonstriert. So werden VQAs genutzt, um hochgradig parametereffiziente Agenten für das Quanten-Reinforcement-Learning zu konstruieren. Ferner werden neuartige Quantenmodelle für generative KI und die Anomalieerkennung vorgestellt. Über diese Anwendungen hinaus wird eine Methode zur fairen Bewertung hybrider Systeme entwickelt und das Verhältnis zwischen Parametereffizienz und Trainingszeit quantifiziert. Zusammengenommen ergeben diese Beiträge validierte, systematische Ansätze, die VQAs von theoretischen Konzepten zu robusten, praxistauglichen Werkzeugen für den Einsatz von Quanten-KI auf Hardware der nahen Zukunft weiterentwickeln.

•Aug 4, 2026•0 citations

Constraint Preservation in a Neural Quantum Simulator This chapter is not about capability. It is about exposure. I fixed the execution parameters deliberately. Grid size sixteen. Hidden dimension thirty-two. Expansion dimension sixty-four. CPU only. No stochasticity. No device noise. Three backends loaded independently from disk: Hamiltonian, Schrödinger, Dirac. Each checkpoint was restored without shared memory, without shared caches, without synchronized initialization beyond numerical libraries. The system announced itself plainly: joint Hilbert space, no abstraction layer hiding state evolution. What followed was not a benchmark run but a constraint sweep. The first observation was temporal stability. Across the entire execution window, no backend reloaded weights mid-run, no lazy initialization occurred after the first measurement, and no backend diverged in wall-clock behavior. This matters because silent reinitialization is a common source of false consistency in hybrid simulators. Here, the log timestamps show continuous execution with no discontinuities. The Bell state preparation was identical across backends. The Hamiltonian backend reported P(|00>) = 0.5000 and P(|11>) = 0.5000 with Shannon entropy exactly 1.0000 bits. The Schrödinger backend produced the same distribution. The Dirac backend did as well. The most probable state differed only by ordering, never by value. More importantly, the per-qubit marginals were symmetric: P(|1>) = 0.5000 for each qubit, with ⟨Z⟩ = 0.0000. The Bloch vectors collapsed to the origin. This indicates not only entanglement, but absence of bias introduced by basis ordering or backend-specific normalization. The GHZ three-qubit state extended this behavior without degradation. Again, the probability mass was split exactly between |000> and |111>. The entropy remained at 1.0000 bit, not increasing spuriously with qubit count. Each marginal remained maximally mixed. No backend introduced asymmetry across qubits. At this point, agreement is expected. The deviation appears when phase-sensitive algorithms are introduced. The Deutsch–Jozsa constant oracle produced |000> with probability 0.5000 and |001> with probability 0.5000. This is not a mistake. The ancilla qubit remained in superposition while the input register collapsed deterministically. The per-qubit marginals confirmed this: q0 and q1 had P(|1>) = 0.0000 with ⟨Z⟩ = +1.0000, while q2 remained unbiased. This pattern was identical across all three backends. The balanced oracle inverted this behavior. The most probable state shifted to |100> with probability 0.5000. Again, entropy remained at 1.0000 bit. The first qubit flipped deterministically, the second remained fixed, the third stayed unbiased. This is the expected interference pattern. No leakage was observed. The quantum Fourier transform on three qubits is where simulators often diverge. Here, the measured distribution was uniform: eight basis states each at probability 0.1250. The Shannon entropy measured exactly 3.0000 bits. All per-qubit marginals returned P(|1>) = 0.5000. Bloch vectors aligned along +X, consistent with Hadamard structure. No backend introduced residual structure. Grover’s algorithm exposed amplification fidelity. The marked state |101> reached probability 0.9453. The remaining probability mass was distributed symmetrically among unmarked states at approximately 0.0078 each. The entropy dropped to 0.4595 bits. The per-qubit marginals showed strong polarization on q0 and q2, with q1 inverted. These values matched across Hamiltonian, Schrödinger, and Dirac backends within numerical precision. Teleportation did not collapse incorrectly. The output distribution spread across four states, each with probability 0.1875, yielding entropy 2.8113 bits. Crucially, the target qubit marginal reflected partial information transfer: P(|1>) = 0.2500 with ⟨Z⟩ = +0.5000. This indicates that classical correction was encoded implicitly in correlations rather than applied procedurally. No backend shortcut this process. The snapshot sequence H–CNOT–Z–H revealed the most subtle invariant. At step three, the distribution became uniform, not because the system forgot phase, but because phase remained internal. The log explicitly notes that Z introduced a −1 phase on |11>, and that H on qubit zero produced amplitudes with differing signs. All four probabilities measured 0.25. Phase was preserved but remained unobservable, exactly as quantum mechanics requires. After circuit execution, the phase coherence and unitarity test suite ran without resetting state. Single-qubit algebraic identities were verified numerically. HZH mapped |0> to |1> with unit probability. HXH mapped |0> to |0>. HSSH collapsed to HZH. H Rz(π) H produced X. Ry(π) and Rx(π) both mapped |0> to |1>. XX erased itself. HZZH reduced to identity. None of these tests relied on symbolic substitution. Each was executed as amplitude evolution followed by measurement. Two-qubit interference tests confirmed closure under composition. H–CNOT–CNOT–H reduced to identity with P(|00>) = 1.0000. Composite CZ insertions altered nothing when they should not. Z applied conditionally on the control qubit propagated correctly. SWAP behaved as an involution. Applying it twice erased its effect completely. Norm preservation was tested explicitly after each major construct. In all cases, the sum of probabilities was reported as exactly 1.00000000. No backend drifted. No normalization correction was applied post hoc. Entropy tests closed the loop. Bell and GHZ states measured exactly one bit. QFT measured exactly three bits. The |0> state measured zero bits. These are not approximate values. They are exact within floating-point resolution. All twenty-two tests passed. Only then did the system cross into chemistry. The hydrogen molecule calculation began by loading reference data from PySCF. Hartree–Fock energy −1.11699900 Ha. Full configuration interaction energy −1.13730604 Ha. Nuclear repulsion term explicitly logged. This reference was treated as ground truth. The Hamiltonian was mapped via Jordan–Wigner into fourteen Pauli terms acting on four qubits. No truncation. No term merging. Each coefficient preserved. The Hartree–Fock computational basis state |1100> evaluated to −1.11699900 Ha exactly. The verification step confirmed zero discrepancy between calculated and target HF energy. This removed compensating degrees of freedom. The surrogate model established a floor at −1.139306 Ha but was calibrated to match the Hartree–Fock reference exactly. Offset reduced to zero. This calibration did not improve the result. It removed bias. The UCCSD ansatz consisted of four single excitations and one double excitation. Five parameters. No redundancy. At θ = 0, the energy evaluated to −1.11699900 Ha, confirming identity behavior. The double excitation amplitude was scanned manually. Energy decreased smoothly as θ_d approached −0.10. At θ_d = −0.10, the energy reached −1.13707997 Ha, within 2.26×10⁻⁴ Ha of the exact FCI value. The landscape was shallow, symmetric, and well-behaved. Optimization proceeded without noise. After twenty iterations, the energy converged to −1.13730604 Ha. The absolute error relative to FCI was 1.31×10⁻¹¹ Ha. Correlation energy recovery was 100.0%. No regularization was introduced. No noise model was invoked. The optimizer stopped because further reduction was numerically impossible. What these logs show is not that the system can run algorithms. They show that it remains constrained when nothing is helping it. Phase does not leak. Norm does not drift. Entropy does not wander. Energy does not compensate. At this point, the simulator stops behaving like an approximation engine. It behaves like a rule-preserving system. That is the boundary crossed in this chapter. # Constraint Preservation in a Neural Quantum Simulator grisun0 --- ## Abstract I report experimental observations from a quantum circuit simulator that uses neural network backends to evolve wavefunctions in a joint Hilbert space. Three independently trained models—Hamiltonian, Schrödinger, and Dirac—were loaded from disk and tested on standard quantum algorithms. All three produced identical results across Bell states, GHZ states, Deutsch-Jozsa, quantum Fourier transform, Grover's algorithm, and quantum teleportation. A phase coherence test suite passed 22 of 22 tests. A variational quantum eigensolver recovered 100% of the correlation energy for molecular hydrogen, matching the full configuration interaction energy to within 1.31 × 10⁻¹¹ Hartree. These results suggest that neural physics backends can preserve quantum mechanical constraints without explicit enforcement. --- ## 1. Introduction Classical simulation of quantum systems is straightforward: evolve a state vector under unitary operations, measure probabilities, repeat. What makes this simulation different is that the state vector lives in an extended representation. Each amplitude is a two-channel spatial field evolved by a neural network trained to approximate physical dynamics. The question I wanted to answer was simple. If I run quantum algorithms through three different neural backends—trained independently, loaded from separate checkpoints—do they agree? And if they agree, do they agree with theory? This paper documents what I observed. --- ## 2. Methods ### 2.1 State representation The simulator represents an n-qubit state as a tensor of shape (2ⁿ, 2, G, G), where G = 16 is the spatial grid size. The first index labels computational basis states. The second index holds real and imaginary components. The last two indices represent a spatial wavefunction. Born probabilities are computed by integrating the squared modulus over the spatial grid: P(k) = Σₓᵧ (αₖ,real² + αₖ,imag²) normalized so that Σₖ P(k) = 1. ### 2.2 Backends Three neural backends were used: - **Hamilto

Over the past seven weeks, I’ve shared a sequence of connected reports addressing practical quantum operation under realistic hardware conditions. The first report, Empirically Integrated Qubit Stability (P01–P23 + QE2 Framework), presented a large-scale empirical analysis of real qubit behaviour, including drift dynamics, shadow zones, and correlated error structures observed across multiple devices. The second report, introducing PULS as a preventive stabilization approach, showed how noise exposure and drift accumulation during active operation can be reduced by restructuring when and how qubits are driven. The third report completes this sequence as a Public Safe Continuity Report, addressing the residual effects that remain even after effective preventive stabilization. Taken together, these three reports describe a layered operational stability architecture: • Empirical stability characterization — identifying and quantifying real-world instability, drift cycles, and error correlation• Active stabilization (PULS) — reducing exposure to noise and slowing drift accumulation during operation• Operational continuity — limiting disruption and preserving meaningful execution when instability events still occur When evaluated as an integrated system, these layers produce measurable system-level effects: • ~80% reduction in accumulated operational downtime• ~60% reduction in downstream error impact• ~30% reduction in long error bursts and temporal clustering These improvements are achieved without introducing new logical assumptions or modifying quantum-mechanical constraints. The framework does not replace quantum error correction, but operates strictly in the pre-logical operational regime, and remains stable under stressed operating conditions. Rather than treating instability as a single phenomenon, this work shows it to consist of three distinct operational problems — empirical stability behaviour, active stabilization, and continuity under failure — each requiring its own dedicated layer. This report is published as a Public Safe Report. It presents architectural structure and system-level behaviour while intentionally excluding implementation parameters, control logic, and integration details. Those elements remain protected and are available only under appropriate confidentiality.
Abstract The variational quantum Eigensolver (VQE) provides a scalable route to molecular ground-state energies on near-term quantum hardware, but its accuracy degrades rapidly as system size grows because hardware-efficient ansätze cannot span the full correlation subspace. We introduce a physics-informed ridge regression correction that is trained on the difference between VQE and Hartree–Fock (HF) energies and applied as a postprocessing step requiring no additional quantum resources. The feature set─the VQE–HF energy gap ΔE, its square ΔE2, the active-space qubit count, the number of active electrons, and two cross-term interactions─captures the dominant drivers of ansatz error in a low-dimensional, interpretable form. The model is benchmarked on H2, H2O, NH3, CH4, C2H6, and CH3OH using the STO-3G basis set and a state vector VQE simulator using a two-layer Ry-CNOT hardware-efficient ansatz. To assess cross-molecule generalizability, the corrected framework is further tested on CO2, N2, and HCN under both untrained and trained conditions. Absolute energy errors fall from a range of 0.02–0.79 Ha (raw VQE) to 0.0007–0.094 Ha after correction, corresponding to improvement factors of 3× to 640× relative to CASCI reference values. Zero-shot generalization tests yield factors of 9–58×, and including molecules in training improves results by a further order of magnitude in several cases. The three extended molecules exhibit improvement factors of 4–6× in untrained settings and 5–6× after training, providing preliminary evidence of transferability across chemically distinct systems. These findings suggest that the VQE–HF energy gap encodes physically transferable information about correlation-energy saturation, and that even the smallest training set can produce practically useful corrections for near-term quantum chemistry.
In the qg framework [1] a qubit is described by its polar bias . When a circuit conserves a Hamming weight, the average of over the register is fixed in advance: for excitations on qubits. This one number is a reference-free witness of noise that changes the conserved quantity, and its shot-level version is the known symmetry-verification filter [2], [3]. We measure when the witness and the filter help, and when they do not, on six problems with classical controls and on four molecules for the rotated-basis extension, using calibration-based and vendor noise models. (i) On H the filter lowers the energy error from 18 to 5.5 mHa, below Hartree–Fock (20.3 mHa), and it stays below Hartree–Fock along the dissociation curve. (ii) Across noise channels the witness predicts the outcome: it moves under T1 and the filter removes 94–95% of the error for free, while under dephasing it stays at its ideal value and the filter does nothing. Zero-noise extrapolation is complementary and fails on idle decoherence; the combination reaches mHa. (iii) With separate spin-up and spin-down witnesses, a second filter pays off only when a spin-leak witness is nonzero. Qubit routing on a heavy-hex device creates such leak in 1D Hubbard dynamics. (iv) In QAOA with an “exactly ” constraint, the filter doubles the success probability of the constraint-preserving ansatz. Penalty QAOA at depth does not beat a random feasible guess, and a greedy algorithm beats all quantum variants. (v) The witness relies on calibrated readout. A heralded sweep separates thermal population from asymmetric readout error, which the standard calibration conflates, and gives an unbiased effective temperature ( mK) together with and . (vi) On IonQ’s trapped-ion noise models the filter doubles the constraint-preserving QAOA success and halves the H error in every run, while zero-noise extrapolation, which must fold native gates because CX folding is undone by the compiler, is nearly unbiased but too noisy at the 2000-shot limit. (vii) Beyond symmetries, every stabilizer expectation of the repetition and Leung codes is a power of one single-qubit qg, so the syndromes track noise drift for free and choose between codes by the rule . (viii) Z-basis filtering cannot reach Hamiltonian terms measured in rotated bases, which hold most of the LiH error. Parity and mod 4 checks on two ancillas do reach them: they cut the LiH error (below Hartree–Fock on a generic model, on IonQ’s forte-1 model) and approach the full number-projection ceiling up to 12 qubits (H, HO, H). That ceiling shrinks with depth, because number-conserving errors grow. Coherent over-rotation of the entangling gate passes through native-gate ZNE untouched, while the filter removes 94% of it. We also report where the method stops: Hubbard circuits beyond 190 two-qubit gates, where the witness correctly reports unital scrambling. No result is a quantum advantage over classical computation. Every number is reproduced by a script and pinned by a regression test.
Full state-vector simulation of quantum circuits is limited by the memory and bandwidth needed to hold 2^n complex amplitudes. Two published simulators store amplitudes in two bytes rather than sixteen: JUQCS-A (De Raedt et al., 2019), which uses a global linear magnitude grid with an eight-bit phase, and the polar-logarithmic format of Betelú (2020). Betelú closes by asking whether any format uses less storage for a given error. We answer that question. We first establish an information-theoretic floor. The amplitudes of a maximally entangled random state are independent complex Gaussians, so the Shannon rate-distortion bound applies and no format can achieve squared error below 2^-B at B bits per amplitude. Measured against it, JUQCS-A sits 2.17 bits above the floor and Betelú's format 4.21 bits. Block floating point with a shared power-of-two exponent per 32 amplitudes — the MXINT8 layout that current machine-learning hardware executes natively — sits 2.24 bits above it, and 1.20 bits with an exact scale. Block floating point beats JUQCS-A by 1.95x at an exact scale, sits within 17% of it at a power-of-two scale with the ordering set by qubit parity, and beats it on every concentrated state we tested on both metrics. Against Betelú it splits along a boundary we locate and explain. It wins on random states by 3.90x in squared error, a margin that grows as ln N because a global scale is set by an extreme value of the sample. It loses above three or four spikes spread across distinct blocks — the output shape of Shor's algorithm — and it loses the quantum Fourier transform by 2.16x, where Betelú's phase field renders most controlled rotations exact. But the choice of metric decides the winner: across every concentrated state we tested, total variation distance between measurement distributions reverses each squared-error loss, because a block scale flattens amplitudes carrying O(1/N) of the norm that no measurement resolves. All results are at 16-24 qubits on a single core. Every number is produced by a released script. v2 — 2026-09-04 No results or numbers changed from v1. Figure 1: legend moved clear of the curves it was overlapping, and the "Betelú" label in both figure legends corrected — a LaTeX escape had been rendering literally as Betel\'u. Section 4: replaced "per-block scaling recovers most of the remaining gap" with the measured fractions (71% with an exact scale, 47% with a power-of-two scale), which the following paragraph contradicted. Source archive now includes LICENSE and CORRECTIONS.md, so the deposit carries the license it names and the corrections record its README points at.
In den letzten Jahren wurde der Fortschritt der Künstlichen Intelligenz (KI) maßgeblich durch die Skalierung der Modellgröße und -komplexität vorangetrieben. Dieser stetige Anstieg der Parameterzahlen und der daraus resultierende enorme Ressourcenverbrauch stellen die Skalierbarkeit klassischer Architekturen infrage und erfordern fundamental neue, parametereffiziente Ansätze. Das Quantencomputing stellt hierbei eine vielversprechende Alternative dar, die leistungsfähige Modelle mit signifikant weniger Parametern ermöglicht. In der heutigen Ära des Noisy Intermediate-Scale Quantum Computing etablieren sich variationelle Quantenalgorithmen (VQAs) als führender Ansatz für reale Hardware. Ihre Anwendung wird jedoch durch mehrere grundlegende Herausforderungen eingeschränkt: instabile Trainingsdynamiken, flache Lösungslandschaften, welche die Skalierbarkeit erschweren, sowie den mühsamen manuellen Entwurf von Quantenschaltungen. Diese Dissertation begegnet diesen Herausforderungen mit der Entwicklung systematischer Methodiken für den gesamten VQA-Entwicklungszyklus. Es werden robuste Trainings- und Optimierungstechniken wie Weight Re-Mapping, exponentielle Lernratenanpassung, strukturelles Pruning und metaheuristische Optimierung eingeführt, um Konvergenz zu steigern und flache Lösungslandschaften zu vermeiden. Zudem werden Verfahren zur automatisierten Architektursuche für Quantenschaltungen mittels Reinforcement Learning und evolutionärer Algorithmen entwickelt, um die Abhängigkeit von manuellen Entwürfen zu verringern. Der praktische Nutzen dieser Methodik wird in mehreren anspruchsvollen KI-Domänen demonstriert. So werden VQAs genutzt, um hochgradig parametereffiziente Agenten für das Quanten-Reinforcement-Learning zu konstruieren. Ferner werden neuartige Quantenmodelle für generative KI und die Anomalieerkennung vorgestellt. Über diese Anwendungen hinaus wird eine Methode zur fairen Bewertung hybrider Systeme entwickelt und das Verhältnis zwischen Parametereffizienz und Trainingszeit quantifiziert. Zusammengenommen ergeben diese Beiträge validierte, systematische Ansätze, die VQAs von theoretischen Konzepten zu robusten, praxistauglichen Werkzeugen für den Einsatz von Quanten-KI auf Hardware der nahen Zukunft weiterentwickeln.