Papers

Although entangled state vectors cannot be fully described in terms of variables localized in space and time, any given entanglement experiment can be built from basic quantum circuit components with well-defined locations. We analyze such quantum circuits and present evidence that the local weak values comprise a covariant tensor associated with each individual qubit. Even if the state is massively entangled, these tensors do not evolve or collapse when other qubits are measured or pass through distant circuit elements. They can therefore be viewed from different reference frames without contradiction. Furthermore, their evolution through any circuit always obeys covariant dynamical rules. Weak values are subject to both past and future constraints, so the covariant quantities can only be determined by considering the entire circuit "all-at-once", as in action principles, incorporating the future measurement basis to avoid the standard no-go theorems. Because these results hold for a set of universal quantum gates, this work lends support to the claim that any quantum circuit can be assigned a realistic, lower-level description compatible with our understanding of classical spacetime.

Quantum evolution is generally expected to drive a quantum many-body system toward equilibrium. This expectation is often justified by the Hilbert-space ergodicity of generic quantum dynamics, namely, the idea that pure-state evolution explores Hilbert space uniformly up to physical constraints. Such a statement can be made rigorous by requiring the associated state ensemble to form the Haar-random ensemble, or its more structured generalization, the Scrooge ensemble. In this Letter, we report the emergence of exact Hilbert-space ergodicity in a continuously monitored quantum many-body system. For any target density matrix $\sigma$, we construct a continuously monitored system for which we rigorously prove that the Scrooge ensemble of $\sigma$ is the unique late-time equilibrium distribution of quantum trajectories. Remarkably, this requires only that the jump operators in the monitoring form a deformed unitary 1-design, a seemingly much weaker condition than full ergodicity. We numerically demonstrate our predictions by simulating continuously monitored systems whose equilibrium states are thermal states. Our results establish a rigorous mechanism for the emergence of Hilbert-space ergodicity and provide a practical route for its investigation on quantum devices.

We present a theoretical study of the evolution of polarization entanglement and quantum coherence in 511 keV photon pairs produced by para-positronium decay during successive Compton scattering events. We start with a maximally entangled Bell state and employ the generalized Stokes-Mueller formalism to derive the two-photon density matrix following single-, double-, and triple-Compton scattering, explicitly considering both polar and azimuthal scattering geometries. Using this framework, we quantify the degradation of quantum correlations through concurrence (as a measure of entanglement) and the $l_1$-norm (as a measure of coherence). Our results demonstrate that entanglement is highly sensitive to the scattering geometry and disappears near right-angle scattering, while quantum coherence remains finite even in regimes where entanglement vanishes completely. These findings provide a unified description of polarization-dependent decoherence in annihilation photon pairs and clarify the distinct roles of entanglement and coherence in realistic two-photon interactions. These results are relevant for quantum-enhanced positron emission tomography and highlight the persistence of quantum resources in scattering-dominated media.

Applications of the proposed primitive: Quantum Primitive for Output-Hiding Function Sharing are discussed for environments in which the parties may have strategic considerations over the generated function value, and therefore may not be mutually trusted. We show that in environments in which the measurement may not be reliably controlled by either party, the primitive permits the parties to generate private, unbiased coins.

Numerical simulation of quantum computing hardware and open quantum systems governed by the Lindblad equation is challenging due to the high dimensionality of the density matrix and the need to preserve fundamental physical properties. In our previous work, we developed an arbitrary-order, low-rank, completely positive and trace preserving (CPTP) method for the Lindblad equation with time-dependent Hamiltonians by nested Picard iteration (NPI). In this work, we develop Gregory NPI schemes, which are CPTP schemes constructed by Gregory-type quadrature on equispaced nodes. The methods, which are of order up to nine, substantially reduce the computational cost compared to our previously proposed NPI schemes with Gaussian quadrature rules, while retaining high-order accuracy and structure preservation. We analyze the stability of the resulting scheme for a physics-based test equation. Numerical experiments verify the convergence of the method and demonstrate the effectiveness of the low-rank approximation. We study the performance of a previously constructed CNOT gate for both closed and open quantum systems.

Applications of the proposed primitive: Quantum Primitive for Output-Hiding Function Sharing, are discussed for secure quantum communications and computation protocols. In particular, QKD applications provide notable enhanced security and efficiency properties relative to status-quo protocols. Additionally, we provide examples when parties may wish to encode joint functions, or decisions, which remain information-theoretically hidden from external parties and those internal to the quantum system, without additional private keys, hidden randomness, or classical communication. In particular these applications may be useful in domains such as; financial transactions, joint signaling or coordination decisions, and navigation systems, among others.

We study entanglement transfer in a minimal model of two qubits that are coupled with one another through an antiferromagnetic Heisenberg exchange ($\tilde{J}$), and where one of them is additionally coupled to a fermionic environment through another antiferromagnetic Heisenberg exchange ($J_{K}$). By tuning the coupling ratio $J_{K}/\tilde{J}$, the system undergoes a quantum phase transition at $T=0$, accompanied by a redistribution of entanglement from the $d'-d$ qubit-subsystem to the environment. Remarkably, the resulting physics exhibits properties that bear analogy with a quantum-information theoretic perspective of the physics of black hole thermodynamics. Carefully selected bipartite and tripartite mutual information measures displays behaviour analogous to the dynamical evolution of black hole entropy, Hawking entropy, and the Page curve expected during the process of evaporation. An effective temperature scale is obtained from the variation of the ground state energy with respect to changes in the bipartite mutual information between the subsystem and the bath. A steady growth of this temperature with the coupling ratio resembles that of the Hawking temperature with the inverse mass of the black hole. Concomitantly, the emergence of non-Fermi liquid behaviour observed near the quantum critical point and in the strong-coupling phase resembles strange-metal-like physics expected near the event horizon from a holographic duality perspective. Our results establish the minimal model as a platform for studying entanglement transfer and information scrambling within a fully unitary quantum framework, and offer new insights into a resolution of the black hole information paradox.

Efficient manipulation of Pauli-algebraic objects is a key bottleneck in the classical emulation and benchmarking of quantum algorithms for chemistry and many-body physics. This bottleneck appears in Hamiltonian construction, variational ansatz preparation, expectation-value and gradient evaluation, and real-time propagation, all of which require repeated Pauli-algebra operations. Here, we present a high-performance Pauli-algebra framework tailored to quantum many-body and quantum-chemical simulations. The framework combines compact binary symplectic encoding, canonical coefficient reduction, and grouped sparse operator representations that exploit shared bit-flip patterns among Pauli strings. The resulting Julia/C++ implementation accelerates Pauli multiplication, Hamiltonian construction, and operator–state multiplication in sparse and symmetry-adapted many-electron spaces. Benchmarks demonstrate efficient Hamiltonian construction, large-active-space VQE and ADAPT-VQE calculations, and real-time variational dynamics on modern multicore CPU and GPU architectures. These results show that structure-aware Pauli-algebra engines provide a scalable classical backend for developing and benchmarking quantum algorithms in quantum chemistry and many-body simulation.

Surface acoustic wave (SAW) phonon coupling with nitrogen-vacancy (NV) center spins in diamond offers a promising platform for on-chip quantum phononic manipulations. Although an ensemble of NV centers coupled to a common SAW phonon mode enables superradiance and collective quantum control, achieving a tunable superradiant phase transition remains challenging. Here, we show that optically driving NV centers level transitions enhances the effective spin-phonon coupling, triggering a SAW phonon superradiant phase transition in the weak-coupling regime. We also demonstrate that above a critical threshold, the driving light rapidly switches on the phonon superradiance–a dynamic effect that persists in finite-number NV ensembles. Our results provide a controllable route to coherent phonon-NV spin manipulation in solid state quantum devices.

Moment closure is a central problem in reduced descriptions of stochastic, kinetic, and quantum dynamics, where equations for low-order observables are coupled to an unresolved hierarchy of higher-order moments. Existing closures usually impose a prescribed form on the distribution or directly truncate the hierarchy, which can become inaccurate or unstable for strongly non-Gaussian states. Here we introduce unit-circle moment closure, which recasts the problem as analytic continuation. Raw moments are mapped to bounded unit-circle moments, whose unresolved tail is reconstructed by a Takagi-Prony procedure from the effective pole structure of a mapped generating function. The resulting continuation yields stable higher-order moments without assuming a fixed distributional ansatz. Illustrative static and dynamical examples demonstrate accurate reconstruction of non-Gaussian distributions and stable evolution of moment hierarchies. Our approach provides a general perspective for moment closure based on analytic structure rather than direct truncation.

Rydberg states of atoms in vacuum are now well recognized as a resource for quantum technologies. Donors in semiconductors also display analogous states, which have been proposed for similar applications. While they benefit from permanent locations in their host crystals, electron-lattice coupling leads to much shorter excited-state lifetimes than for neutral atoms in vacuum. Here we provide a quantitative description of donor-phonon kinetics, creating a basis for engineering donor systems in realistic material stacks for quantum devices. Our theory incorporates both form factors for the Rydberg states, which given their large extents in real space provide strong selectivity in momentum space, and tabulated deformation potentials for all six phonon branches throughout the Brillouin zone. By confronting this framework with carefully controlled time-resolved free electron laser measurements, we show that the widely quoted position of the silicon conduction band minimum, $k_0$, is inconsistent with observed donor relaxation rates and that quantitative agreement is obtained for a value further from the X-point than commonly assumed. This stringent experiment-theory comparison establishes donor relaxation as a precision metrology for conduction band parameters and scattering processes in silicon, with consequences spanning from quantum devices to classical electronics.

Near-optimal solutions to Maximum Independent Set on dense graphs sit in local optima that greedy correction and maximality extension cannot escape. We encode near-optimal seeds as a uniform quantum superposition on ancilla qubits and evolve them under an excitation-preserving variational ansatz that holds the search inside the feasible Hamming-weight subspace. A preprocessing stage of spectral reordering and distance-based sparsification, together with history-guided post-processing of the sampled bitstrings, takes the method to 200 nodes. The ansatz entangles the seed branches and the bond dimension of the simulated state grows with circuit depth, which places deeper circuits outside the reach of exact matrix product state simulation at the bond dimensions available to us. Larger instances therefore need quantum hardware. Measured on the data register alone, the superposition behaves as a classical mixture over the seeds, so coherence between the branches has to be created and then looked for. We do this with a CRZ phase layer and post-selection on the ancilla qubits, which brings the branches into interference and exposes the cross terms. The idea is to see if interference between near-optimal seeds widens the range of independent sets the circuit returns. Standard VQE with this pipeline recovers the certified optimum for instances up to 125 nodes, and these run on ibm_marrakesh with parameters transferred from noiseless simulation. The ancilla construction is introduced for the sizes past that point. On a 180-node hard instance the superposition recovers the certified MIS. Five 200-node instances are solved to the certified optimum, and on the 400-node brock400-1 benchmark the method reaches size 25 against a certified optimum of 27. These are the largest hard general-graph instances we know of where a gate-based variational algorithm optimises the full circuit directly.

We reveal that engineering electrode-coupling configurations can fundamentally reshape coherent transport phenomena in quasiperiodic quantum systems. Leveraging nonequilibrium Green's function theory, we systematically analyze charge and heat transport, as well as current fluctuations, in a magnetic-flux-threaded quasiperiodic Su-Schrieffer-Heeger ring with both symmetric and asymmetric multi-site reservoir couplings. Contrary to the conventional expectation that optimal transport is achieved near the homogeneous-hopping limit, our results reveal that multi-site lead coupling fundamentally reshapes the transport landscape, extending the regime of enhanced transport deep into the topological phase. Strikingly, asymmetric source-drain coupling induces a disorder-assisted conducting phase where quasiperiodic modulation enhances, rather than suppresses, charge and energy transport. Magnetic flux exerts a dual influence: it activates additional interference-mediated transmission channels that amplify transport while simultaneously suppressing the disorder-induced re-entrant conducting regime. Furthermore, we uncover a flux-driven migration of the optimal transport window with increasing disorder strength, shifting from the topological regime toward the trivial-hopping regime. This behavior highlights the intricate interplay among quasiperiodicity, dimerization, magnetic-flux-induced quantum interference, and the geometry of the system-reservoir coupling. Collectively, our findings position coupling engineering as a powerful paradigm for the rational control of nonequilibrium transport in quasiperiodic materials and chart a route toward quantum device configurations in which transport characteristics can be precisely tuned via the interplay of disorder, topology, and quantum interference.

Spin-photon interfaces based on solid-state defects are key building blocks for scalable quantum networks and hybrid quantum platforms. Optimizing light-matter coupling in these systems requires precise knowledge of the optical transition dipole polarization, yet for many promising quantum emitters this quantity is hard to determine and therefore remains poorly characterized. Here, we develop a framework for reconstructing electric transition dipole polarization in spin-1/2 solid-state defects directly from ensemble spectroscopy. The approach combines the response of photoluminescence spectra to magnetic field, optical polarization, and strain. Applied to erbium ions in silicon, a particularly challenging system containing multiple crystallographic subsites, the framework identifies strain-induced shifts as the origin of asymmetric ensemble spectra and enables simultaneous determination of the optical dipole polarization and strain-orbital coupling tensor. The resulting model predicts how cavity-ion coupling depends on crystallographic orientation and magnetic-field direction, which we verify using single erbium ions coupled to a nanophotonic cavity. Together, these results establish a broadly applicable route for extracting microscopic properties of solid-state quantum emitters from ensemble spectroscopy and for engineering optimized spin-photon and spin-phonon interfaces.

We have performed first-principles hybrid density functional theory calculations to study the formation and properties of two-dimensional electron or hole gases (2DEG or 2DHG) at s-Si/SiGe interfaces. For small Ge concentrations $x < 0.25$, we find a novel type of band alignment with no offset in the conduction bands, implying that a 2DEG cannot be formed, though a 2DHG can. In contrast, for $x > 0.25$ the band alignment suggests that either a 2DEG or 2DHG can be formed. The electronic band structure features two nearly degenerate 2DEG states at the bottom of the conduction bands, and two 2DHG states at the top of the valence band. These states can be accessed by appropriate doping and gating. Charge density plots of these states show that they feature carriers confined to the near vicinity (2–3 atomic layers) of the interface. Calculated effective masses are anisotropic, being markedly so for the 2DHG states, and in excellent agreement with experiment. This property can be exploited to create a 1D carrier gas. Our results are especially important for s-Si/SiGe-based semiconducting spin qubits for quantum computing applications.

Bipartite Gaussian boson sampling (BipartiteGBS) produces output probabilities governed by squared permanents of submatrices of arbitrary complex matrices, matching the nonsymmetric structure of directed graphs. Most GBS-based graph algorithms, however, rely on symmetric hafnian structure and are formulated for undirected problems. Here we propose a BipartiteGBS-based framework for directed-graph heuristic optimization. We introduce Max-Perm as a canonical optimization task for BipartiteGBS and derive a closed-form sampling enhancement factor relative to uniform classical sampling in this idealized setting. We then use permanent-biased BipartiteGBS samples to guide a genetic algorithm for the celebrated directed Hamiltonian cycle problem. Numerical experiments on Erdős–Rényi random directed graphs show that the resulting BipartiteGBS-enhanced algorithms improve success rates over a standard genetic algorithm and yield longer valid paths when no Hamiltonian cycle is found, while ablation tests indicate that BipartiteGBS-guided initialization is the dominant contributor. These results show how permanent-based photonic sampling can provide useful algorithmic guidance for asymmetric combinatorial search.

We investigate the generation of optical torque in coherently prepared multi-level atomic media driven by a vector vortex beam composed of two counter-rotating components carrying opposite orbital angular momenta, $+l\hbar$ and $-l\hbar$. We consider a three-level $\Lambda$ configuration and a four-level tripod configuration. Using a perturbative steady-state solution of the optical Bloch equations, we obtain analytical expressions for both linear and nonlinear contributions to the optical torque. The results show that the torque is strongly controlled by atomic coherence, including the initial population imbalance and the relative phase between the vortex components. Nonvanishing torque can arise even when the two components have equal amplitudes, due to coherence-induced asymmetry in the atomic response. In the tripod configuration, the presence of a strong control field leads to electromagnetically induced transparency, which suppresses the torque near resonance and shifts the dominant response to finite detunings. These results establish a route for controlling light-induced rotational dynamics in atomic media using vector vortex fields, with potential applications in coherent optical manipulation and angular-momentum-based control in quantum systems.

We propose a dynamical approach to identifying ground-state quantum phases through short-time dissipative cooling. Rather than determining the phase by preparing highly accurate approximations to ground states, we prepare a representative state of a candidate phase and monitor the early-time response of phase-sensitive observables under cooling dynamics tailored to the target Hamiltonian. For a class of phase-decision problems in which the relevant observables can be inferred from the low-energy manifold, and with jump operators implementable using only short-time Hamiltonian simulation, the dissipative evolution rapidly suppresses high-energy components and drives the system into a low-energy manifold whose observables already reveal the underlying ground-state phase, well before mixing to the steady state. We demonstrate this strategy for the frustrated $J_1$–$J_2$ Heisenberg chain, the Kitaev honeycomb model, and the XXZ chain, including Berezinskii–Kosterlitz–Thouless and topological phase transitions. In particular, coarse filter resolutions and short evolution times suffice to recover phase-sensitive quantities such as the Luttinger parameter and topological diagnostics. We further provide theoretical justification that cooling dynamics with such jump operators can rigorously prepare low-energy manifolds for free-fermionic and free-bosonic systems, and investigate this mechanism for interacting fermionic systems. Our results suggest that phase decision is a plausible target for future utility-scale studies on early fault-tolerant quantum devices.

Estimating potential future exposure (PFE) for path-dependent derivatives, such as FX Target Redemption Forwards (TARFs), represents a formidable computational challenge due to the demand of nested Monte Carlo simulations. We present a hybrid quantum-classical framework that leverages Iterative Quantum Amplitude Estimation (IQAE) to address this via a reduced-order counterparty credit risk model. Our methodology maps the non-linear TARF payoff – including cumulative gains and knock-out features – into a quantum circuit via a two-step formulation, whereby a first-step percentile is computed classically and then used to condition quantum evaluation of subsequent exposure. We employ discretisation of the FX process and a linearised additive approximation of dynamics to enable implementation on current quantum platforms. Developed via the Classiq platform and validated on NVIDIA CUDA-Q and Amazon Braket SV1, our approach achieves relative errors of 1%-8% against classical benchmarks at the 97.5% and 99% confidence levels. While discretisation constraints and approximate monotonicity assumption may introduce bias and limit recovery of the full exposure distribution, our framework offers a tractable testbed for quantum acceleration. Scaling analysis suggests that $\sim$300 logical qubits could enable full 52-week exposure estimation, reducing sample complexity for tail-risk estimation via amplitude estimation at the cost of increased circuit depth.

The entanglement of Gaussian continuous-variable (CV) states is fully determined by the state's second moments. In contrast, some entangled non-Gaussian states evade every second-moment criterion, and non-Gaussian entanglement detection remains an experimental challenge. The $p_3$-PPT criterion detects entanglement using moments of the partial transpose of the density matrix. This criterion was recently extended to CV systems using photon-number-resolving detectors and multi-copy interferometry; here we introduce a single-copy homodyne protocol that detects bipartite CV entanglement via the same criterion. Unbiased U-statistic estimators for the partial-transpose moments $p_2$ and $p_3$ are constructed directly from randomized homodyne data and used to evaluate the $p_3$-PPT entanglement witnesses: a linear one for detection, and a quadratic one whose violation yields a dimension-free lower bound on the entanglement negativity. The protocol estimates $p_2$ and $p_3$ up to additive error $\varepsilon$ at Fock cutoff $N$ from $O((N+1)^{14/3}/\varepsilon^2)$ measurements at fixed confidence. We demonstrate the protocol on six families of Gaussian and non-Gaussian states, reaching $95\%$ empirical one-sided detection probability from $\sim 10^3$ to $10^4$ homodyne measurements for states with $\bar{n} \approx 2$, placing non-Gaussian entanglement detection within reach of current homodyne experiments.

Planar Germanium has emerged as a promising platform to build spin-based large scale quantum computers. By exploiting the anisotropic hyperfine interaction of holes in Ge, qubits with long T2* have been recently realized. While the performance of single qubits is still more or less limited by 73Ge nuclear spin fluctuations, the site-to-site variation of qubit sweet spot becomes obstacles to maintaining high fidelity of each qubit across the whole wafer. To achieve high performance Ge-based quantum circuit, it is therefore essential to eliminate the origin source of hyperfine noise. In its Silicon counterparts, reduction of 29Si abundance enables exceptional high-fidelity operation. In contrast, hole qubits based on isotopically purified Ge have not been demonstrated. Here, we report the synthesis of high quality 2-dimensional hole gas (2DHG) with enriched 70GeH4 precursor. Due to the suppression of nonzero spin nucleus, the qubits' T2* on the sweet spot is moderately extended beyond 20 us, surpassing the previous best reported Ge hole qubits. More importantly, the qubits' T2* off the sweet spot is enhanced to above 3 us, enabling single qubit gate fidelity exceeding 99.9% in both operating regimes. Hahn-echo spectroscopy further resolves a finite-frequency nuclear-noise channel that is distinct from the conventional Larmor-linked hyperfine response. We associate this channel with quadrupole-modified dynamics of residual 73Ge nuclei sampling local electric-field gradients near the Ge/SiGe interface. Its field scaling and angle-dependent visibility are consistent with a qubit-visible quadrupolar nuclear-noise component transduced through the anisotropic hyperfine interaction of Ge holes. These results establish isotopically purified planar Ge as a high-coherence scalable platform for hole spin qubits and provide a spectroscopic probe of interfacial quadrupolar nuclear dynamics.

We consider the encoding of graph problems as Quadratic Unconstrained Binary Optimization (QUBO) problems, which are solvable by either quantum or classical annealers. Yet, the class of problems encoded as QUBO problems has not previously included nowhere-zero flows. Nowhere-zero flows are related to Tutte's $5$-flow conjecture and appear in many contexts in graph theory. We provide an encoding of nowhere-zero flows as a QUBO Hamiltonian and prove the correctness of the construction. Our construction yields a Hamiltonian $H_{\mathrm{mod},k}$ whose ground state has zero energy if and only if the graph $G$ has a nowhere-zero $\mathbb Z_k$-flow. By Tutte's equivalence theorem, zero ground energy is equivalent to $\varphi(G)\le k$, and the zero-energy degeneracy is given by the flow polynomial $F(G;k)$. In particular, when the ground-state energy is zero, this is also the ground-state degeneracy. The construction uses one-hot variables to represent the edge flow residues modulo $k$ and auxiliary variables to represent the per-vertex modular quotient. We prove that the correctness of the construction is independent of the choice of orientation, root vertex, and positive penalty weights. We verify the construction on $59$ examples of graphs and values of $k$ that include both yes-instances and no-instances. We exhaustively sweep orientations and root choices on selected robustness instances and test a finite suite of positive penalty weights. The resulting Hamiltonian is implemented using the dimod.BinaryQuadraticModel class, which is compatible with the D-Wave Ocean SDK. Quantum-hardware runs and claims about potential speedup using these devices are left to follow-up work.

Fermionic topological orders can host 't Hooft anomalies with no bosonic counterpart. We identify a new sixteen-fold family of (2+1)D fermionic topological orders, forming a fermionic analogue of Kitaev's sixteen-fold way. This family is distinguished by the mod 16 't Hooft anomaly of a $\mathbb{Z}_2$ one-form symmetry, generated in each theory by a single nontrivial $\mathbb{Z}_2$ anyon. This intrinsically fermionic anomaly permits anyon spins that are forbidden in bosonic phases; the simplest new example is an Abelian fermionic topological order containing a single $\mathbb{Z}_2$ Abelian anyon of spin 1/8. Each theory can be realized as the gapped boundary of a (3+1)D fermionic symmetry-protected topological (SPT) phase protected by the $\mathbb{Z}_2$ one-form symmetry, which acquires a $\mathbb{Z}_{16}$ classification once the spacetime spin structure is twisted by the one-form symmetry. We realize these phases microscopically via lattice models built from Walker-Wang models coupled to local fermions.

Temporal non-analyticities in the rate function of the Loschmidt echo manifests a class of dynamical quantum phase transitions (DQPTs) that has emerged as a powerful framework for understanding far-from-equilibrium many-body dynamics. While such DQPT has been extensively studied theoretically in spin-boson systems such as the Dicke model, their experimental observation remains elusive. In particular, the dynamics of DQPT in asymmetric spin subspaces and under the influence of spin dissipation are largely unexplored. Here, we report an experimental study of temporal DQPT in a generalized Dicke model using a trapped-ion quantum simulator. By coupling a linear chain of $\rm{^{40}Ca^{+}}$ ions to a collective center-of-mass motional mode, we probe the quench dynamics starting from both symmetric and asymmetric initial states. We extract the rate function and identify temporal turn-around points that are in quantitative agreement with theoretical predictions. Additionally, we investigate the impact of spin dissipation on these dynamics. Our results establish an experimental platform for probing complex many-body out-of-equilibrium phenomena and advance the development of hybrid oscillator-spin quantum simulators.

Understanding how macroscopic thermodynamic behavior emerges from microscopic quantum dynamics remains an open problem. While equilibration of quantum observables is well established, thermodynamics also relies on variables not directly associated with linear operators, but which are defined instead as functions of expectation values. Whether and how such derived quantities inherit equilibration properties is an open question. Here, we establish that any continuously differentiable function of equilibrating expectation values also equilibrates. We apply this result to a bipartite isolated system, showing that the entropy and conjugate variables of each subsystem – defined through Jaynes' maximum entropy principle – equilibrate. Moreover, with the assumption that their equilibrium properties depend solely on local conserved quantities, we show the dynamical maximization of the total entropy, enforcing equality of conjugate variables across subsystems. These results provide a direct dynamical justification for entropy maximization and the emergence of thermodynamic equilibrium conditions, showing that fundamental principles of thermodynamics follow from the unitary evolution of quantum systems.

Parameterized quantum circuits (PQCs) provide a flexible substrate for hybrid quantum machine learning (QML), but their practical value on Noisy Intermediate-Scale Quantum (NISQ) devices remains an empirical question, especially because training depth and scale can introduce optimization challenges such as barren plateaus. Here we study how the number and topology of two-qubit entangling gates in the feature-map stage influence a fixed hybrid QNN workflow for classifying strong versus weak epitope-receptor binding in Porcine Reproductive and Respiratory Syndrome (PRRS) vaccine design. The dataset consists of docking-derived binding affinities for N=80 9-mer epitopes, labeled as Strong or Weak binding, and partitioned into training, validation, and test subsets using a 40:30:30 split. We compare a classical CNN benchmark with a hybrid Embedding-QNN architecture under four feature-map configurations: a non-entangling Z feature map, an all-to-all high-entanglement ZZ feature map, and two interleaved nearest-neighbour entanglement patterns of low and high depth. Among the configurations tested, the high-entanglement ZZ feature map is seen to provide the strongest evidence of reduced training-set overfit, with a lower training area under the accuracy curve (AUAC) and the highest test/training AUAC ratio, while preserving competitive test-set accuracy. These results do not establish a general QML advantage, but they suggest that feature-map entanglement topology is a meaningful design variable for sparse biological screening tasks and warrants further evaluation with additional metrics, larger datasets, and noise-aware or hardware-based experiments.

The controlled creation of specific quantum states is a highly challenging field of research that is also in high demand with applications in various quantum technologies. Due to its difficulty, it has been historically impossible to use quantum states as a rendering target for complex scenes and visualizations, especially in the Noisy Intermediate-Scale Quantum (NISQ) era where operations are noisy and only a few qubits are available. Even so, in this paper we propose a new, quantum compatible method for general purpose rendering by extending Composable Function Systems (CFSs) to quantum architectures. We discuss limitations in the classical implementation of CFSs and the physical maps necessary for adapting the classical method to quantum by showing the generation of object primitives and transformations of such objects, including duplications and smears, some of which are topologically non-trivial. Our results reveal the possibility of a speed-up for this, specific method on quantum hardware and our method has been used to create the first video rendered on quantum architectures.

We propose a new Brownian bridge construction for our newly devised Coherent States Path Integral Monte Carlo algorithm. We apply it to the numerically exact calculation of the thermodynamic properties of the Helium fluid on a plane at low non zero temperature. We find very good agreement with the conventional plane waves path integral Monte Carlo results.

We formulate a continuity equation for the Schrödinger equation in the complex space. We define a complex momentum by normalizing the complex current by the particle density. This momentum is a quantum analog of the classical, kinematic momentum analytically continued into the complex plane. The kinematic momentum and the gradient of the wavefunction's phase each represent a fluid-like flow in the complex plane; the phase-gradient flow is incompressible. The zeros of the wavefunction give rise to simple poles in the momentum. The poles manifest as irrotational vortexes in the phase-gradient flow, while critical points of the wavefunction present as rigid body-like rotational flows of the kinematic momentum. A discrete nature of elementary excitations comes about inherently because the quantity of the poles is automatically integer. An exact quantization condition is subsequently formulated, which reduces to the Bohr-Sommerfeld condition in the semiclassical limit. We establish a priori that the Bohr-Sommerfeld condition must be exact for the Harmonic Oscillator. We show that the kinetic energy is a sum of contributions of the average value and fluctuations, respectively, of the kinematic momentum. The zero-point vibrations within bound states are solely due to the fluctuations of the momentum and manifest as rigid-body flows at infinity. The momentum poles – and hence the wavefunction's zeros – can be viewed as emergent, consistent with the remarkable property of quantum entanglement exhibited by standing wave solutions of the Schrödinger equation.

We introduce the dimension-restricted prepare and broadcast (PAB) scenario, which generalizes standard prepare-and-measure frameworks. Here, the system prepared by a sender undergoes a broadcasting transformation before being locally measured by multiple receivers. We develop a hierarchy of classical, quantum, and nonsignalling models describing this scenario, characterize their corresponding correlation sets, and derive new families of Bell-like inequalities together with linear and semidefinite programming methods for their certification. First, assuming shared randomness, we prove that the hierarchy collapses into a single set whenever we consider only one measurement per party. Then, considering multiple possible measurements, we show that PAB scenarios allow the activation of nonclassicality, revealing genuinely nonclassical features in resources that admit classical descriptions in standard prepare-and-measure or Bell settings.