Papers

The gate and qubit requirements of quantum computations of electronic structure have been extensively studied. However, the quantum resources present in electronic ground states, as measured by entanglement and magic, remain less well understood. We study the relationship between correlation in electronic structure Hamiltonians and magic as measured by the 2-stabilizer Renyi entropy (2-SRE). Perturbative calculations show that the 2-SRE of a given state is proportional to its overlap with a reference stabilizer state. In the context of quantum chemistry, this links the magic of electronic structure ground states to their Hartree-Fock weight, an established measure of electronic correlation. We then show that the 2-SRE of post-Hartree-Fock ground states is proportional to the correlation energy they recover. We explore this connection through the contextual subspace (CS) method. We present a theoretical framework showing that the CS method can be used to monotonically vary the magic of approximate CS ground states, and we prove that the correlation energy recovered by the CS ground states is proportional to the magic present in the approximate ground state. We present simulation results using 190 molecular species under Jordan-Wigner encoding at a range of bond lengths. The linear relationships between magic and correlation are robust across the Hamiltonians in our dataset, but break down at bond lengths beyond the Coulson-Fischer point, where Hartree-Fock fails to capture key physical features of the true ground state wavefunction. By establishing linear relationships for both correlation energy and Hartree-Fock reference weight with the 2-SRE, we conclude that for weakly- and moderately-correlated electronic structure Hamiltonians, the correlation is directly represented by 2-SRE, and thus by the magic.

Two-qubit-gate error dominates the failure budget of near-term quantum circuits, so the decomposition chosen for each Toffoli (CCX) gate should minimize hardware two-qubit infidelity, not gate count. The cheapest decompositions - relative-phase and approximate Toffolis - are only correct in context: their residual phase or bounded error must be cancelled or absorbed downstream. We present the first compiler pass that selects a per-Toffoli decomposition to minimize a two-qubit-infidelity error budget. It admits each context-dependent decomposition only when an exact, instance-specific equivalence check certifies its validity in that circuit context, coupling an error-budget objective with per-instance verification and closing the gap between context-aware-but-unverified and verified-but-context-free optimizers. The central result is a safety one: pattern-matched relative-phase substitution is silently incorrect. Our verifier flags 66 library rewrites of a deployed open optimizer as non-equivalent without a context check, and count-greedy substitution silently corrupts 6 of 12 benchmark circuits; the verification gate certifies 0 errors while still applying every valid decomposition. The two-qubit-gate reduction is real but workload-dependent: up to 39.5% fewer two-qubit gates and 36.7% lower infidelity over exact-only on a compute/uncompute-heavy suite (approx. 39%/35% versus Qiskit opt-3 and tket), and 15.6% aggregate on a larger 12-24-qubit suite, with decision-diagram checking certifying every substitution past the exhaustive-verification limit. At current superconducting and trapped-ion error rates, the certified substitutions lower estimated circuit infidelity by 36-43%, and on a quantum state-resetting circuit, the pass removes 48.8% of the native two-qubit gates, every substitution verified.

In Ref. 1, a method was given for self-consistently generating sum-of-squares decompositions of quartic fermionic Hamiltonians. Perturbation theory was used to generate a useful choice of cubic operators in this sum-of-squares. On a range of model problems, this method, which is only a fragment of degree-six sum-of-squares, was able to outperform the full degree-four sum-of-squares in both speed and accuracy. Unfortunately for applications, many problems in chemistry have strong density-density interaction terms, as well as moderately strong density-dependent hopping and spin-spin interaction terms, limiting the power of the perturbative choice of the cubic operators. Here we propose a method for generating these decompositions in the presence of these strong interaction terms, hopefully extending the range of applicability of this method.

We propose PIFLUX, a single-molecule localization scheme combining deep-subwavelength plasmonic illumination with widefield detection. Interference between counter-propagating gap plasmons and a normally incident optical field generates an illumination pattern whose position can be tuned through the plasmon phase while preserving its spatial period. A Cramér-Rao analysis shows PIFLUX reaches few-nanometer precision matching MINFLUX while doubling that of SIMFLUX over a micrometer field of view, and a maximum-likelihood estimator confirms this on a synthetic nuclear pore complex.

Recently, Doosti et al. introduced the problem of distributed quantum state verification, where $m$ distributed nodes are given a copy of an unknown state $\rho$, and can send limited one way communication to a central node, who has a complete description of a known state $\sigma$. They ask how many distributed nodes $m$ are required, before the central node can succeed at distinguishing whether $\rho=\sigma$ or $\|\rho-\sigma\|_1\geq\varepsilon$ with high probability. In the setting where only quantum communication is allowed, Doosti et al. exhibit conditional lower bounds in both the public and private-coin settings, and a matching upper bound in the public-coin setting. We extend these results, and show unconditional lower bounds for when both classical and quantum communication are permitted. We show the public-coin lower bound is tight by giving an algorithm with a matching upper bound. We also show an almost tight upper bound in the private-coin setting when only quantum communication is permitted.

The circular Unruh effect is the excitation of a detector moving along a planar circular trajectory within an electromagnetic vacuum. We demonstrate that the magnetic dipole transitions in an atom, acting as the detector, dominate the electric dipole transitions. Our analysis of both free-space and cavity schemes shows that the sensitivity to the circular Unruh effect can be maximized by balancing the minimization of mode volume against the resulting decrease in mode density. Moreover, we propose a novel measurement scheme that uses the atom's multilevel structure to suppress the spontaneous emission rate, thereby enabling the experimental detection of the circular Unruh effect.

Electron and positron can form a meta-stable bound state called positronium that decays via pair annihilation. We show how polarization-dependent Compton scattering can be used to verify that the two annihilation photons in the spin-zero case (para-positronium) are emitted in a maximally entangled Bell state. Our theoretical approach based on two-photon density matrices connects concepts from relativistic quantum electrodynamics and quantum information theory.

Following the holographic proposal that identifies the growth rate of Krylov complexity with the proper radial momentum of an infalling massive probe, we study Krylov complexity in Lifshitz and hyperscaling-violating backgrounds. For pure Lifshitz geometries, we find quadratic complexity growth for all values of the dynamical exponent. For hyperscaling-violating backgrounds, we show that the hyperscaling-violating exponent directly controls the late-time growth exponent, with a special limiting case exhibiting oscillatory behavior with a logarithmic envelope. We also extend our analysis to top-down string constructions, confirming the robustness of the scaling behaviors obtained in the bottom-up approach. Our analysis establishes that the momentum-Krylov correspondence extends naturally to non-relativistic holographic settings and remains well-defined despite the causal pathologies of Lifshitz spacetimes.

Merging optomechanics and polaritonics opens stimulating perspectives like the giant enhancement of optomechanical interaction and the enrichment of optomechanics with effective nonlinear photons. The experimental implementation of these concepts has however remained elusive. Here we report on the resonant optical control of polaritonic optomechanical resonators constituted of semiconductor disks embedding quantum wells. Whispering gallery photons and quantum well excitons strongly couple, leading to the emergence of polaritons that couple to the mechanical vibrations of the disk. We perform resonant optomechanical frequency response experiments on these resonators, modeled introducing a minimal set of constitutive equations, from which we extract the polariton-modified optomechanical coupling $g_0$ and the polariton nonlinearities. We observe a boost of $g_0$ by more than a decade compared to bare photons, reaching to a record $g_0$ for whispering gallery resonators of $22$ MHz, and analyze experimentally and theoretically its evolution as function of the polariton's composition. We also measure a clear hierarchy of three polaritonic nonlinearities, again analyzed as function of polariton composition, establishing a bridge between past unconciliated reports in polaritonics. Grounded on experimental and theoretical foundations, resonant polaritonic optomechanics is set ready for an optomechanical exploration of quantum fluids of polaritons.

We propose a quantum collocation framework for approximating solutions of one-dimensional linear and nonlinear boundary value problems. The method formulates the search for admissible solutions as a residual-based quantum search over a discretized ansatz space, where candidate solutions are evaluated through residual conditions imposed at collocation points. A residual-threshold oracle is constructed that acts jointly on spatial and parameter registers. This joint oracle structure leads to amplification dynamics that decompose into a coherent superposition of spatially conditioned amplitude-amplification processes rather than a single global amplification mechanism. We derive the corresponding amplification geometry and show that the success probability is governed by a weighted combination of spatially dependent amplification angles. Furthermore, we prove that the reversible residual oracle can be implemented with gate complexity polynomial in the logarithm of the number of collocation points, while retaining the quadratic search acceleration associated with amplitude amplification in the parameter space. We analyze how the spatially dependent oracle structure influences the amplification dynamics and corresponding success probabilities. Furthermore, we investigate how discretization, ansatz expressivity, oracle tolerance, and finite-precision effects influence both approximation quality and amplification behavior. Numerical experiments validate the theoretical predictions and illustrate the resulting search dynamics across different discretization and precision regimes.

Exceptional points, where eigenvalues and eigenvectors coalesce, are a defining feature of non-Hermitian systems and have been extensively observed in photonic, atomic, and condensed matter systems. However, they have received little attention in quantum chromodynamics (QCD), which is the fundamental theory of quarks, gluons, and hadrons. We propose that imaginary magnetic fields provide a simple realization of non-Hermitian dynamics in hadronic systems. Based on two theoretical approaches, a hadronic effective Lagrangian and a constituent quark model, we compute mass spectra of neutral mesons and find exceptional points separating the real-spectrum and complex-eigenvalue regimes. In small fields, the real spectrum exhibits level attraction between hadronic states, whereas in larger fields, hadrons are deconfined, which is a signature of a field-induced inverted potential. Our findings open a new avenue for studying QCD dynamics in non-Hermitian environments.

Gravity-induced entanglement has been proposed as a method for testing the non-classical nature of gravity via tabletop experiments. While most existing proposals are restricted to the Newtonian limit, the frame dragging effect offers access to genuinely post-Newtonian features of the gravitational interaction and remains comparatively less explored. Here, we study gravity-induced entanglement generated by frame dragging in an interferometric setting and compute the entanglement phase between the rotational degrees of freedom of a source mass and the paths of a particle in two complementary ways: (i) via Schrödinger evolution with a quantized Lense-Thirring Hamiltonian in the large angular momentum limit, and (ii) via the on-shell action of linearized quantum gravity within the stationary phase approximation. Both approaches yield the same entanglement phase, consistent with the proper time difference between the interferometer arms. The path integral derivation further reveals how gravitational retardation modifies the entanglement phase, thereby making the local, relativistically causal linearized-gravity description explicit. Under the standard locality/mediator assumptions used in existing arguments, the resulting entanglement would witness non-classicality of the gravitational interaction.

Periodically driven quantum systems offer remarkable flexibility in tailoring effective Hamiltonians and synthetic band structures. However, such driving also induces heating and dynamical instabilities that limit the coherence and lifetime of many-body states. Here, we demonstrate that these instabilities can be suppressed by employing multi-tone driving schemes. Using a Bose-Einstein condensate of cesium atoms in an optical lattice, we experimentally explore two approaches: pulsed driving composed of odd harmonics and two-tone driving with tunable amplitude and relative phase. We show that both methods allow independent control of the effective tunneling amplitude and Peierls phase factor, while significantly reducing phonon excitation and the resulting rapid decay of the condensate. Numerical simulations and theoretical modeling based on Bogoliubov-de Gennes equations confirm the suppression of unstable modes under optimized driving conditions. Our results establish multifrequency drives as powerful tools for stabilizing driven many-body systems and pave the way toward robust Floquet engineering with interactions.

The motion of a particle in vacuum near macroscopic bodies gives rise to a Magnus-like contribution to the nonequilibrium Casimir-Polder force. This effect originates from the interplay between particle dynamics and material-modified electromagnetic quantum fluctuations, inducing in the particle a direction-dependent angular momentum coupled to the electromagnetic field spin. The resulting drift force is proportional to the cross product of the particle's angular and translational velocities, revealing a rotational transport component in the nonequilibrium Casimir-Polder interaction. Our results establish a striking connection between quantum fluctuations-induced forces and the classical Magnus effect in fluid dynamics.

The core step in quantum simulations is typically matrix vector multiplication $\phi = \Hmat \psi$. Executing this step is limited by memory requirement to store the Hamiltonian. We present a memory-scalable, hardware-adaptive matrix-free framework for applying large operators on vectors without materializing the full matrix on a single accelerator. The operator is represented through a block-procedural interface: blocks may be generated, loaded, cached, distributed, or applied directly only when their action is needed. For quantum simulation, it provides the core kernel for quantum operations. An adaptive planner selects block size, cache strategy, GPU grouping, row distribution, and task parallelization from memory and workload estimates. We describe analytic, measured, and learned planning strategies that choose between procedural generation, partial caching, full caching, and row-distributed caching. The method removes the requirement that the full dense matrix fit in the accelerator memory. This shifts large simulations from a fixed memory barrier to a tunable balance between block generation, cache reuse, data movement, parallel scheduling, and numerical accuracy.

We study scattering in noninteracting multi-terminal quantum wire junctions and show that junctions with dihedral symmetry can exhibit exact zero-reflection points for $N \ge 4$ terminals. By analyzing the scattering matrix, we identify these reflectionless points in the $(E,t')$ parameter space, where $E$ is the incident particle energy and $t'$ is the junction hopping amplitude. These points exhibit an even-odd dependence on $N$ and converge asymptotically to a common limiting value in the large-$N$ limit. We show that the reflectionless points are characterized by an integer winding number associated with the phase of the reflection amplitude, providing a topological description for their stability against weak on-site disorder. We also consider junctions with broken time-reversal symmetry and find that a magnetic flux can induce additional reflectionless points, including for the $N = 3$ case. For a four-terminal junction threaded by a $\pi$-flux, we identify a unique parameter regime in which the reflection amplitude vanishes over the entire energy band. Finally, we discuss experimental signatures through the behavior of Friedel oscillations and examine the stability of these reflectionless points in the presence of weak interactions.

Three ingredients of the elementary Landau-Zener problem determine the familiar expression $a_{LZ}\equiv\exp\left[-\pi/(2\epsilon)\right]$ for the asymptotic value of the probability amplitude for remaining in the initial level: (i) A wave whose phase is determined by the product of a contour integral over a simple pole at the origin of the complex plane and the inverse of twice the scaled chirp parameter $\epsilon$. (ii) An asymptotic limit of the associated path connecting the points $\pm 1$ along the real axis and circumventing the pole in the upper half-plane, and (iii) a half-circle in the lower half plane enclosing together with the asymptotic path the pole. The Cauchy theorem immediately provides us with the value $\ii\pi$ of the asymptotic contour, and thus with $a_{LZ}$. Our analysis demonstrates not only that $a_{LZ}$ is the consequence of a logarithmic phase singularity but also explains why the Markov approximation also leads to $a_{LZ}$.

We introduce a squeezed-state pulse-position modulation (S-PPM) format, where the empty slots are squeezed vacuum states and the pulse slot is a displaced squeezed state. Based on this property, we propose an inverse-squeezing conditional pulse-nulling (IS-CPN) receiver. In the ideal case, inverse squeezing maps S-PPM into an equivalent coherent-state PPM signal with a large pulse energy, leading to a closed-form expression for the receiver error probability. We further analyze IS-CPN under common phase diffusion using a finite-path MAP formulation with phase-averaged likelihoods. Numerical results show that IS-CPN outperforms conventional CPN under the same energy constraint and remains advantageous under phase noise and finite photon-number resolution. These results demonstrate that combining squeezed-state modulation with inverse-squeezing conditional nulling can improve photon-efficient optical communication.

As Quantum Machine Learning (QML) transitions toward practical implementation, the field faces a critical architectural bottleneck that challenges the fundamental assumptions of classical statistical learning theory. In classical deep learning, increasing model capacity typically risks overfitting. However, this study advances a counter-intuitive paradigm: unstructured contemporary QML architectures suffer from a profound state of quantum underfitting, driven by the "expressivity-trainability paradox." We demonstrate that the vast Hilbert space capacity of Parameterized Quantum Circuits (PQCs)-traditionally chased as the source of quantum advantage is the direct mathematical cause of Barren Plateaus (BPs), where gradient landscapes become exponentially flat. By synthesizing recent breakthroughs in Dynamical Lie Algebras (DLAs) and Geometric QML, we establish a comprehensive framework linking the algebraic dimension of circuit generators to their optimization dynamics. Furthermore, we empirically validate this framework on a non-linear binary classification task, illuminating a uniquely quantum manifestation of the bias-variance tradeoff: while unstructured architectures achieve near-perfect training accuracy via unscalable parameterization (quantum overfitting), embedding group-theoretic geometric priors acts as a structural regularizer. By restricting the DLA growth to a polynomial regime, our symmetry-preserving approach sacrifices raw memorization capacity to guarantee scalable, gradient-rich training landscapes, offering a robust roadmap for "Trainability-by-Design" in scalable quantum neural networks.

Practical hardware limitations, including finite transmitter and receiver filter lengths as well as the finite resolution of digital-to-analog and analog-to-digital converters, lead to mode mismatch and degrade the performance of continuous-variable quantum key distribution systems. To address this, we develop a machine learning-based end-to-end optimization framework that jointly optimizes transmitter pulse shaping and receiver matched filtering. The approach employs reinforcement learning under realistic hardware constraints, including a limited number of filter taps, finite digital-to-analog and analog-to-digital converter resolution, analog low-pass filtering, and the optimal mean photon number. By mitigating mode mismatch and accounting for implementation constraints, the proposed method improves overall system performance. Simulation results demonstrate enhanced secure key rates compared to conventional approaches, demonstrating the effectiveness of the proposed framework.

Interference is commonly regarded as the most direct manifestation of the superposition principle. This association is natural for the linear Schrödinger equation, where coherent alternatives combine at the level of probability amplitudes. However, the situation becomes less transparent when nonlinear couplings are present, or when the field is only partially coherent. In this work, we argue that a more robust organizing principle is provided by the local flow generated by phase variations. In this sense, phase-induced flow acts as a unifying mechanism for interference-like dynamics in nonlinear and partially coherent Schrödinger systems. The discussion is developed from a hydrodynamic, or Bohmian, perspective, understood here as a practical probing tool rather than as an additional ontology. Three representative situations are considered: interfering Bose–Einstein condensates described by the Gross–Pitaevskii equation, nonlinear Schrödinger dynamics obtained by modifying the quantum-potential contribution, and partially coherent Airy beams described through their cross-spectral density. Although these systems differ in physical origin and mathematical implementation, they share a common dynamical structure: density-related observables are shaped by velocity fields determined by phase, or ensemble-phase, information. From this viewpoint, interference-like traits, localization, self-acceleration and coherence loss can be interpreted in terms of the preservation, deformation or breaking of the symmetries displayed by the underlying flow. This provides a compact way of connecting interference, nonlinear dynamics, classicality, coherence loss, and structured-light propagation within a single trajectory-based framework.

A prevailing intuition holds that quantum protocols using only free states confer no operational advantage. This intuition is contradicted by free-state discrimination gaps in which restricted measurements fail to optimally distinguish even orthogonal free states. Known instances include nonlocality without entanglement and, more recently, nonstabilizerness without magic. We trace these examples to a single convex-geometric mechanism: whenever the set of free measurements is closed, convex, and strict subset the set of all measurements, and the free states is a convex set with an interior, a gap-witnessing ensemble can be drawn entirely from the free states. The resulting gap is operationally rigid: no finite-dimensional assistance – catalyst or quantum memory – can asymptotically improve the discrimination rate beyond the single-shot restricted limit. By contrast, non-free ensembles admit memory-assisted attacks that fully erase the gap, exposing a sharp operational asymmetry between free and resource-carrying ensembles.

Motivated by the temporal dynamics identified in the $(1+1)$ Carroll-Schrödinger theory, we derive a post-Carrollian Schrödinger dynamics in flat Klein space with signature $(2,2)$. Starting from the tachyonic Klein-Gordon equation in double-polar coordinates and removing a spacelike carrier, the spatial radius behaves as an effective evolution parameter, whereas the temporal two-plane $(t_1,t_2)$ serves as the equal-radius configuration space. The additional time direction supplies an $SO(2)$ temporal angular momentum $J$, produces temporal vortex sectors, and gives the centrifugal contribution to the post-Carrollian momentum $P_{\mathrm{PC}}=E_{\tau}^{2}/(2M_{\mathrm{eff}})+J^{2}/(2M_{\mathrm{eff}}\tau^{2})$ in the Hamilton-Jacobi limit. We determine the regular Bessel modes, Gaussian packets, oscillator spectrum, radial $SU(1,1)$ tower, equal-$r$ continuity equation, $\mathfrak{sch}(2)$ symmetry algebra, radial-ordered propagator, and the metaplectic organization of the quadratic sectors. Effective flat connections on the temporal configuration plane give Aharonov-Bohm, Landau, and Fock-Darwin analogues, while the two-body relative sector admits anyonic boundary conditions on the punctured temporal plane. As a curved extension, we derive a branch-dependent carrier reduction and apply it to an illustrative $SO(2,1)$-symmetric Kleinian Schwarzschild exterior, where the Kleinian gravitational source produces a lensing-type angular deviation on the temporal plane.

We introduce a unified, computable measure of multipartite entanglement based on the spectral properties of an entanglement graph and its associated entanglement matrix. This framework quantifies quantum correlations among arbitrary subsystems and partitions of a composite system. We prove that the resulting spectral entanglement measure satisfies the fundamental requirements of entanglement measures. Furthermore, we derive a generic multipartite monogamy relation that extends residual entanglement beyond qubit systems and introduces spectral residual entanglement for arbitrary multipartite states.

We introduce the average confidence width $\Delta_{a} x=\int_0^1\Delta_{c} x (\theta)\,\mathrm{d}\theta$: the confidence width $\Delta_{c} x(\theta)$—the measure of the smallest region of position space carrying probability $\theta$—averaged over all confidence levels. It is the first moment of the decreasing rearrangement of the position density, an $L^1$ (mean-absolute-deviation) measure of localization, and $\Delta_{a} x\,\Delta_{a} p$ is dilation invariant. We study the uncertainty relation $\Delta_{a} x\,\Delta_{a} p\ge c^\ast\hbar$ for pure and mixed states. A mean–entropy argument with the Bialynicki-Birula–Mycielski relation gives $c^\ast\ge\pi/e$, while the ground state $\psi_0$ of the Fourier-invariant operator $|x|+|p|$, with $E_0=1.1040744$, gives $c^\ast\le E_0^2\approx1.2190<4/\pi$: Gaussian states, which saturate the Heisenberg–Kennard and entropic relations, are not optimal here. We prove that $\psi_0$ is nondegenerate, nodeless, even, self-dual, and symmetric decreasing, with algebraic tails, and that the sharp hybrid relations $\Delta_{a} x \cdot 2 \langle |p-p_0| \rangle \ge E_0^2\hbar$ and $2\langle|x-x_0| \rangle \cdot\Delta_{a} p\ge E_0^2\hbar$ hold for all states, so $\Delta_{a} x\,\Delta_{a} p\ge E_0^2\hbar$ whenever either density is symmetric unimodal; the second variation of $\Delta_{a} x\,\Delta_{a} p$ at $\psi_0$ is positive definite apart from the symmetry directions. Extensive numerical minimization, including coherent states, never finds a product below $E_0^2\hbar$ and consistently returns $\psi_0$, supporting the conjecture $c^\ast=E_0^2$. In two dimensions the relation is sharp and saturated by Gaussians, identifying the dimension as the origin of the non-Gaussian optimizer in one dimension. Finally, we discuss other complementary pairs and the relation of the width to interferometric predictability and visibility.

In both classical and quantum error correction, heralded erasures are known to be easier to tolerate than unheralded general stochastic errors. Whilst an established benefit of loss-dominant quantum architectures such as photonic qubits, this fact has received renewed interest, with a pivot towards reconstructing other architectures to be erasure-dominant, such as dual-rail transmons. This work investigates exploiting these 'erasure qubits' in the near term by using postselection as a technique for error mitigation, wherein circuit shots detecting any erased qubits are discarded from the computational ensemble and repeated. Firstly, we outline a numerical framework for representing circuit-level erasure noise and present 'erado', an open-source library capable of simulating erasure noise and postselection. Secondly, we investigate the effects of both erasure noise and noise in the erasure checks themselves on the quantum Fourier transform (QFT), in the additional presence of gate depolarising noise. A worked example is provided of postselection fully mitigating against the erasure channel for erasure check error rates less than 3.0%. We also show how a postselected dual-rail system can surpass a fundamental noise floor at the kiloquop scale where a comparable single-rail system cannot, justifying this approach in the NISQ regime before (and, perhaps, combined with) the practical arrival of QEC.

We introduce a programmable optical parametric amplifier (OPA) synthesizer that, under a heralded photon-number-resolving framework, generates high-fidelity cubic phase states and amplifies Schrodinger cat states. By systematically exploring both the catalytic configuration, where the idler input and output contain the same number of photons ($m=n$), and non-catalytic configurations ($m\neq n$), we discover two qualitatively different functionalities. First, with a coherent-state signal input, our protocol generates cubic phase states with fidelity exceeding 0.99 across a broad range of $(m,n)$ configurations. Second, using a Schrödinger cat state as the signal input, the same framework amplifies the cat state: an input cat with amplitude $\alpha_{\mathrm{in}}\le 1$ is transformed into an output squeezed cat with $\alpha_{\mathrm{out}}\ge 2$ while maintaining fidelity above 0.99. The catalytic configuration preserves the input parity and restores the idler state, whereas non-catalytic configurations enable parity-flipping amplification with higher success rates. Moreover, the amplified output can serve as a seed for subsequent amplification rounds, offering a self-seeding pathway to progressively larger cat states. Our protocol requires only moderate-gain OPA operation and low-order photon-number-resolving detection, providing a flexible and experimentally accessible platform for cubic phase state preparation and amplified squeezed cat state generation.

We introduce a hybrid approach utilising a quantum machine learning surrogate model to approximate the non-linear collision dynamics of the LBM. It effectively offloads the non-unitary operations that challenge pure quantum solvers. The expressivity of the surrogate is built on the ability of parameterised quantum circuits to implement partial Fourier series, with data re-uploading extending the spectrum of representable frequencies. Unlike previous approaches with a fixed relaxation parameter, the surrogate recovers the complete Bhatnagar-Gross-Krook (BGK) collision dynamics across the full physically admissible range of relaxation without retraining. We reassess the relevance of standard variational quantum circuit (VQC) metrics, including expressibility, entanglement, and effective dimension, by relating them directly to task-specific surrogate performance and identifying the key architectural parameters that determine approximation accuracy. The proposed surrogate is validated against the classical BGK collision operator using established benchmark problems, including the Taylor-Green vortex for evaluating energy dissipation and the double shear layer for assessing shear-driven instabilities and nonlinear flow evolution. Our results demonstrate that the hybrid model achieves high accuracy and generalisability while closely replicating classical solutions. These findings suggest that hybrid quantum-classical strategies offer a practical path toward realising the potential of quantum computing in fluid engineering.

Quantum Markov models have been successfully used to accurately model various physical quantum systems in fields such as quantum optics, optomechanics and superconducting circuits and they provide the basis for (measurement-based) quantum feedback control. However, the quantum Markov assumption is a strong one and it is not expected to hold for general quantum systems of interest. The projection operator approach is one approach that has been developed to model non-Markovian quantum systems by considering its embedding in a larger Markovian quantum system, but mainly in the context of quantum master equations for the dynamics of the unmonitored reduced quantum state of a quantum system. This approach was recently adapted for continuously measured non-Markovian quantum systems, which enables open-loop control but did not yet consider the presence of feedback of the stochastic measurement record, deriving non-Markovian SDEs for the evolution of the projected state of the Markovian embedding. This paper generalizes these stochastic equations to the setting of stochastic feedback based on the continuous-measurement record and shows that the equations take the same form but that previously deterministic terms become stochastic ones which depend on the measurement record, as would be intuitively expected. The stochastic equations are obtained for a generalized class of measurements that includes continuous (possibly adaptive) homodyne and photon counting measurements.

This work presents a full-wave Green's function framework for modeling collective and coherent single-photon emission from multiple quantum emitters embedded in complex electromagnetic structures. Starting from a transverse modal completeness relation of modified Langevin noise formalism, we derive a closed set of coupled equations for population dynamics and frequency-resolved field amplitudes in the single-excitation regime. Since the electromagnetic reservoir is not traced out at the level of the dynamical amplitudes, the emitted single-photon dynamics can be modeled within the same closed set of equations without Markovian approximation in open and dissipative environments. We demonstrate that finite-bandwidth truncation of the spectral density leads to systematic deviations in coherent dispersive interactions, even when dissipative rates appear converged. To restore causal consistency, we introduce a counter-term compensation scheme that restores the missing dispersive contributions without modifying the retained non-Markovian memory kernel. To validate the scheme and demonstrate the practicality of the proposed framework, we present numerical examples ranging from benchmark configurations to a three-dimensional dispersive ring-resonator structure via finite element method. These capabilities provide a practical route for rigorously incorporating full-wave electromagnetic simulations into non-Markovian multi-emitter quantum electrodynamics, enabling predictive modeling of collective emission, coherent energy exchange, and single-photon radiation in realistic open structures.