Papers

There has been intensive research on increasing the utility and performance of Parameterized Quantum Circuits (PQCs) in the past couple of years. Owing to this research, there are now several inductive biases available to a quantum algorithms researchers to design a good circuit for their chosen task. In this paper, we focus on the problem of finding performant initial parameters for a given PQC. Different from previous research that focuses on finding the right distribution, we focus on finding the hyperparameters for any given distribution. To that end we introduce an evolutionary-search based algorithm that finds optimal hyperparameter given a PQC and quantum task. Our empirical results indicate that our algorithm consistently leads to selection of performant initial parameters tuned specifically to the ansatz and the quantum task leading to faster convergence and performance. More importantly, our algorithm does not negatively affect the barren plateau phenomenon. In other words, the initial parameters suggested by algorithm do not worsen the gradient variance scaling for a given initializing distribution.

Fast qubit reset is essential for qubit reuse in the noisy intermediate-scale quantum computing era, yet it conflicts with the weak decoherence required for high-fidelity computation. We solve the time-optimal reset problem for a frequency-tunable qubit coupled to a structural environment under realistic spectral and control constraints. The optimal strategy consists of a switch–restore–switch sequence, where the qubit is moved from a low-decoherence computational configuration to a high-decoherence restoring configuration and then returned for reuse. For superconducting qubits in four representative environments, this strategy reduces the reset time from typically $\gtrsim\SI{100}{\nano\second}$ to $\SI{20}{\nano\second}$, about $40\%$ of a typical two-qubit gate time, while achieving a reset precision of $10^{-5}$. Our results identify environmental spectral structure as a practical resource for rapid, high-fidelity qubit reset and provide a design principle for qubit reuse on qubit-limited processors.

In continuously monitored quantum systems, the feedback protocol of García-Pintos, Liu, and Gorshkov reshapes the arrow of time: a Hamiltonian $H_{\mathrm{meas}} = r A / \tau$ applied with gain $X$ tilts the distribution of measurement trajectories, with $X < -2$ producing statistically time-reversed outcomes. Why this specific Hamiltonian achieves reversal, and how the mechanism relates to score-based diffusion models in machine learning, has remained unexplained. We compute the functional derivative of the log path probability of the quantum trajectory distribution directly in density-matrix space. Combining Girsanov's theorem applied to the measurement record, Fréchet differentiation on the Banach space of trace-class operators, and Kähler geometry on the pure-state projective manifold, we prove that $\delta \log P_F / \delta \rho = r A / \tau = H_{\mathrm{meas}}$. The García-Pintos feedback Hamiltonian is the score function of the quantum trajectory distribution – exactly the object Anderson's reverse-time diffusion theorem requires for trajectory reversal. The identification extends to multi-qubit systems with independent measurement channels, where the score is a sum of local operators. Two consequences follow. First, the feedback gain $X$ generates a continuous one-parameter family of path measures (for feedback-active Hamiltonians with $[H, A] \neq 0$), with $X = -2$ recovering the backward process in leading-order linearization – a structure absent from classical diffusion, where reversal is binary. Second, the score identification enables machine learning (ML) score estimation methods – denoising score matching, sliced score matching – to replace the analytic formula when its idealizations (unit efficiency, zero delay, Gaussian noise) fail in real experiments.

We study the dynamics of $N$ photons in a Fock state, initially located inside one cavity, and coupled by an optical fiber to a second cavity. The entanglement of the photons is monitored by projective measurements, repeated with a fixed time step. This approach is applied to the formation of a photonic N00N state. We calculate the probability of the transition of $N$ photons from the left to the right cavity and the probability of the return of $N$ photons to the left cavity under repeated projective measurements. The entanglement is analyzed for the N00N state by its fidelity and its phase sensitivity, while for the entanglement between the states in the two cavities the entanglement entropy is calculated. In addition, we study the monitored evolution of photons in a single cavity, which are coupled to a single qubit, using the Jaynes-Cummings model. Photon entanglement is analyzed in terms of the entanglement entropy. In all these cases we find that entanglement is sensitive to the details of monitoring protocol, which can be used to control photon entanglement for specific applications.

We study the magnetic-field dependence of the interactions between two alkaline-earth(-like) Rydberg atoms, ${}^{88}$Sr and ${}^{174}$Yb. Considering the pair of Rydberg states $|ns,{}^3S_1,m_J\rangle$ and $|(n+1)s,{}^3S_1,m_J\rangle$, we show that the effective Hamiltonian takes the form of an XXZ-type quantum spin model, as in the alkali-atom case [M. Kunimi and T. Tomita, Phys. Rev. A 112, L051301 (2025)]. We find that the behavior of the anisotropy parameter for ${}^{174}$Yb at zero magnetic field is significantly different from that for other atomic species. This behavior arises from the interplay of strong spin-orbit coupling and the resulting multichannel redistribution of Förster defects in ${}^{174}$Yb. We systematically calculate the interaction parameters of the XXZ model in the presence of a magnetic field and show that they can be tuned by the field. As applications to quantum many-body problems, we investigate one-dimensional systems in the large-anisotropy regime and show that the folded XXZ model can be realized in ${}^{174}$Yb systems without fine-tuning of the field. We also investigate two-dimensional square-lattice systems and show that a supersolid phase can emerge in the ground state at the mean-field level.

In a previous manuscript we addressed the possibility of generating a reflection in a region of spacetime inside a null surface under electromagnetic gauge transformations. In this manuscript we will deal with the possibility of full inversions in a region of spacetime inside a null surface under electromagnetic gauge transformations. Since new tetrads whose construction depends on electromagnetic gauge were introduced it has been proved that the physical kinematic nature of spacetime can be altered as proved previously. In this case we will see how to create a field configuration that will reverse the flow of time. We will prove that we can turn a future directed timelike vector into a past directed timelike vector by physical means. Experiments devised in order to achieve time reversal or full inversion have been never discussed before.

We quantify the excitonic coupling in the homodimer of the dimeric Venus fluorescent protein using a quantum-classical hybrid workflow. Because the anionic chromophore carries appreciable double-excitation character, we obtain its site energy from domain-based local-pair-natural-orbital similarity-transformed equation-of-motion coupled-cluster theory (DLPNO-STEOM-CCSD), which places the in-protein bright $\pi\pi^{*}$ state at $523.90~\mathrm{nm}$, in good agreement with the experimental absorption maximum near $515~\mathrm{nm}$ and validated against an explicit triples-corrected EOM-CCSD(fT) benchmark; single-excitation time-dependent density functional theory (TDDFT) blueshifts this state. Applying a transition-density coupling (TDC) formalism to the STEOM transition density and simulating the VenusA206 tandem-dimer construct (two $\beta$-barrels covalently joined by an inter-domain linker), the covalent tether holds the A206 interface docked throughout $1~\mathrm{ns}$ of room-temperature dynamics. We obtain a thermally averaged coupling of $J = 32.8 \pm 1.6~\mathrm{cm^{-1}}$ across the $1000$-frame trajectory (Davydov splitting $2|J| = 65.6 \pm 3.1~\mathrm{cm^{-1}}$), representing a $19\%$ near-field enhancement over the ensemble point-dipole estimate ($J_{\mathrm{PDA}} = 27.6 \pm 1.3~\mathrm{cm^{-1}}$) at the $24.69 \pm 0.32~\text{\AA}$ chromophore centroid separation. Furthermore, we argue that a separation of timescales resolves the apparent theoretical tension between intermediate experimental excitonic couplings and the highly decoherent biological environment. While it has been hypothesised that the fluorescent protein $\beta$-barrel scaffold sustains coupling by shielding the chromophore from thermal fluctuations, we emphasise that the separation of timescales applies irrespective of the exact degree of environmental decoherence and dissipation.

We characterize where, in the noise parameter space of the uniform dephasing–depolarizing channel $\mathcal N_\gamma^\delta=\mathcal E_\delta\!\circ\!\mathcal D_\gamma$, a deterministic, nonlinear, target-informed denoising heuristic yields its largest fidelity gain. The procedure pulls the off-diagonal magnitudes of the noisy state toward those of a known target, with efficiency set by a SWAP-test-purified catalyst; it is not a quantum channel, and target access is an explicit classical resource, so the results describe a benchmark procedure, not blind error correction. Our main tool is the covariant purification map $\mathcal P_\mathrm{cov}(\rho)=(\rho+\rho^2)/(1+\mathrm{Tr}\,\rho^2)$, an exact closed form for one SWAP-test purification round (a rederivation of symmetrization purification: Barenco et al., Cirac–Ekert–Macchiavello) that reduces the catalyst to a scalar eigenvalue iteration. With it we derive the $d\to\infty$ fidelity-gain peak location on Haar-random pure states (Theorem~3): $\gamma_{\rm peak}(d)\to\gamma^\star(r,\delta)$, with $\gamma^\star(2,0.1)=0.4725$ and limiting magnitude $0.2262$. Bootstrap-quantified sweeps to $d=256$ are consistent with both limits. The peak is resource- and protocol-dependent: a catalyst-only reference moves it from $\approx0.50$ to $\approx0.34$, and one purification round instead of two to $\approx0.39$. Bell and uniform $d=4$ states admit unique-peak theorems for the $r=0$, $\delta=0$ protocol member. All results are reproducible from the open-source organic-qc-bench package with seed~42.

Bell inequalities play a central role in certifying quantum correlations and underpin protocols such as device-independent quantum key distribution. However, enumerating all Bell inequalities for a given scenario remains intractable beyond the simplest cases, as it requires solving a computationally hard facet enumeration problem on the associated Bell polytope. We propose the Adjacency Sampling method, which builds on the Adjacency Decomposition method but sacrifices completeness for speed. On previously solved Bell polytopes, the method reproduces every known class of inequalities. For scenarios where no complete enumeration exists, it greatly exceeds existing partial results: in $\mathcal{L}_{3,3,3,3}$ we obtain over $1.29 \times 10^8$ classes, more than 25 times the previous count; in $\mathcal{L}_{4,5,2,2}$ we nearly triple the known list to 49\,358 classes; and for $\mathcal{L}_{4,6,2,2}$ we report over 4.3 million classes.

We present a general and efficient approach to compute phase-resolved multidimensional spectra of anharmonic molecular polaritons, based on a semiclassical evolution of the molecular Hamiltonian and cavity field in the large-$\mathcal{N}$ limit of many molecules coupled to a confined photonic mode. By systematically expanding the response in both amplitudes and phases of the input fields, our method enables a transparent and computationally simple construction of phase-cycled two-dimensional single- and double-quantum polariton spectra from the underlying nonlinear signal components. Here, phase cycling acts as an analogue of phase matching with oblique pulses, allowing for the isolation of the contributing nonlinear pathways in Liouville space. We specialize to vibrational polaritons and benchmark the method through direct comparison with experimentally measured single-quantum spectra, providing an explanation for the longstanding puzzle of the polariton bleach effect observed at short waiting times. Further, we show how the imprint of various types of anharmonicities on the double-excitation manifold can be directly probed and analyzed through double-quantum coherence spectroscopy. Taken together, our results establish a practical and powerful framework for the modeling and interpretation of nonlinear spectroscopic experiments on strongly coupled light-matter platforms and for guiding the design of cavity-enhanced molecular platforms.

We introduce a qubit- and gate-efficient higher-order unconstrained binary optimization (HUBO) encoding for graph partitioning problems requiring label-count minimization. This widely applicable class of problems includes minimum graph coloring, minimum k-cut, and community detection. To the best of our knowledge, this is the first work to address the optimization versions of these problems in a quantum setting, rather than only their decision counterparts. Our construction encodes each vertex variable in a number of bits logarithmic in the number of available labels—the information-theoretic minimum—and employs a novel lexicographic penalty system that implicitly minimizes partition count without requiring dedicated indicator variables. We derive provably sufficient conditions on all penalty coefficients, including those arising from Rosenberg quadratization, guaranteeing feasibility and optimality of the lowest-energy solution. Analogous conditions are derived for a one-hot encoding to enable controlled comparison. We also show that our encoding reduces two-qubit gate count per layer of the quantum approximate optimization algorithm (QAOA). Benchmarking on a quantum annealer demonstrates that our logarithmic encoding greatly improves solution quality and time-to-solution for minimum graph coloring relative to one-hot encoding, with greater advantage as problem size increases.

Estimating local observables in Gibbs states is a central problem in quantum simulation. While this task is BQP-complete at asymptotically low temperatures, the possibility of quantum advantage at constant temperature remains open. The Sachdev-Ye-Kitaev (SYK) model is a natural candidate: at any constant temperature, its Gibbs states have polynomial quantum circuit complexity and are not described by Gaussian states. Rigorous analyses of the SYK model are difficult due to the failure of known techniques using random matrix theory, cluster expansions, and rigorous formulations of the quantum path integral and replica trick. Despite this, we give a rigorous proof of a quasipolynomial-time classical algorithm that estimates SYK local thermal expectations at sufficiently high constant temperature. Our result introduces a new Wick-pair cluster expansion that we expect to be broadly useful for disordered quantum many-body systems.

We investigate how laser-driven, cooperative dipole-dipole interactions in weakly trapped atomic arrays give rise to self-organized configurations. Starting from an analytically tractable two-emitter system, we identify the possible steady-state spatial arrangements accessible to the atoms. We then extend this analysis to larger ensembles in both linear and ring geometries. In linear chains, we demonstrate the emergence of topologically nontrivial dimerized configurations across a range of initial interatomic spacings. In ring geometries, we find that the system undergoes self-organized contraction and expansion, enabling access to length scales below those set by the trapping lattice. Our results demonstrate that collective light-matter interactions in free space can spontaneously generate modified ordered geometries, even when the emitters are initially separated by distances larger than their transition wavelength.

Accurate modeling of ion-molecule reaction networks is essential for understanding the chemical evolution of planetary ionospheres, particularly for giant planets where proton-transfer chains drive atmospheric composition. However, predicting reaction rates in these ultracold environments remains a challenge due to the non-trivial interplay between vibrational dynamics and quantum tunneling. In this work we present a chaos-diagnostic framework that integrates multireference electronic structure theory, Adiabatic Gauge Potentials (AGP), and Random Matrix Theory (RMT) to characterize the microscopic dynamics of proton transport. Using the formation of H+3 and the proton-bound cluster H+5 as representative model systems relevant to Jovian atmospheres, we demonstrate that the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity. This diagnostic approach offers a generalizable, data-driven metric for identifying vibrationally gated pathways in complex astrochemical networks, providing a theoretical basis for refining kinetic models of planetary and interstellar plasmas

We investigate the effects of coupling a quantum-magnetic cavity field to molecules. Our high-precision auxiliary-field quantum Monte Carlo calculations capture the effect of the cavity field in the presence of electron correlations, and their interplay and competition. In H$_2$, we find that a strong enough cavity coupling makes the original bound ground state metastable, along with inverting the singlet-triplet gap. In ring molecules (e.g., H$_n$), the magnetic cavity coupling stabilizes symmetric geometries. As a consequence, open-shell rings such as H$_4$, H$_8$, or C$_4$H$_4$, which would undergo Jahn-Teller distortions outside of the cavity, obtain exotic spin or ring-current polarized, antiaromatic ground states. These effects are enhanced by increasing the molecule concentration inside the cavity. Our results suggest cavity quantum electrodynamics beyond the long-wavelength approximation as a promising avenue for cavity-altered chemistry.

We explore the application of the Variational Quantum Eigensolver (VQE) to investigate the ground state properties, particularly the entanglement entropy, of the Transverse Field Ising Model (TFIM) in one, two, and three dimensions, considering systems of up to 27 spins. By benchmarking VQE results against exact diagonalization and analyzing the entanglement properties across different system sizes, we assess the algorithm's effectiveness in capturing faithful ground state. Using results of TFIM, we also investigate how VQE's expressivity and optimization influence the simulation of highly entangled quantum states. We employ different ansätze: the hardware-efficient EfficientSU2 from Qiskit, the physics-inspired Hamiltonian Variational ansätz (HVA) and HVA with symmetry breaking, and benchmark their performance using energy variance, entanglement entropy, spin correlations, and magnetization. We further discuss the implications for scaling these methods to larger quantum systems.

This work investigates the emergent thermalization regimes in a chaotic Tavis-Cummings (TC) model and their implications in quantum spectroscopy. While the TC model is a cornerstone of cavity quantum electrodynamics, traditional treatments often overlook many-body effects that arise in the thermodynamic limit. We utilize the Eigenstate Thermalization Hypothesis to demonstrate that a non-integrable excitonic Hamiltonian within the material manifold drives local thermalization. By tuning the polariton splitting $g$, we observe two dynamical regimes: a thermalizing regime at low interactions driven by quantum chaos and ergodicity, and a non-thermalizing regime at high interactions where strong coupling suppresses ergodicity. We further show that these regimes have direct implications on output photon statistics, specifically influencing the correlation times $\tau_c$ of the cavity population and the second-order correlation function $g^{(2)}(t+\tau)$. We propose that entangled-biphoton spectroscopy serves as an ideal experimental platform to probe these effects and to allow the characterization of the underlying many-body exciton-coupling disorder $\sigma$ through coincidence measurements of the output. Taken together, these results exploit a naturally occurring many-body phenomenon to bridge theoretical predictions with experimental observables.

The Berry phase encodes the geometry of a closed Hamiltonian path, complementing the dynamical phase determined by the accumulated energy. We uncover an adiabatic error-cancellation principle that arises from this geometric character and has no counterpart for energy-based quantities. To extract the Berry phase while canceling the dynamical phase, we combine finite-runtime evolutions generated by $\pm H$ along the loop. This construction exactly cancels the leading $O(T^{-1})$ and all higher odd-order nonoscillatory phase errors. Richardson extrapolation further reduces the leading residual error to an oscillatory contribution whose amplitude is set by the endpoints of the loop alone. Beyond this deterministic cancellation, runtime randomization suppresses the remaining oscillatory contribution, reducing the bias after $r$ Richardson levels to $O(T^{-(2r+2)})$ for fixed $r$. The same average also suppresses the bias from an imperfectly prepared guiding state, so that the admissible preparation infidelity no longer depends on the target accuracy \(\varepsilon_B\). Combining these principles, we obtain a randomized Hadamard-test algorithm for Berry phase estimation that reduces the coherent evolution time to \(O(\varepsilon_B^{-1/(2r+2)})\) while retaining the standard \(O(\varepsilon_B^{-2})\) sample complexity. This turns the geometric origin of the cancellation into an algorithmic advantage and makes Berry phase estimation a promising task for early fault-tolerant quantum computers.

Quantum resources are increasingly integrated into high-performance computing (HPC) and cloud environments, but quantum high-performance computing (QHPC) software stacks remain isolated, often proprietary, full-stack solutions lacking common interfaces across runtime, resource management, orchestration, and execution layers. This paper analyzes nine production QHPC stacks and identifies common design patterns and emerging requirements, covering deployment models, application interaction patterns, SDK support, and readiness for fault-tolerant operation. The survey exposes consistent needs in runtime abstraction, resource management, interconnect semantics, and observability. Based on these findings, we propose the open quantum-HPC software ecosystem ( openQSE) reference architecture as a first step toward unifying the state-of-the-practice. openQSE defines a set of layer boundaries that allow different implementations to interoperate while preserving deployment flexibility, and is structured to support both current noisy intermediate-scale quantum (NISQ) workloads and future fault-tolerant quantum computing (FTQC) systems without changes to upper-layer application interfaces.

We construct a quantum low-density parity-check code family from a length-$512$ Calderbank–Shor–Steane base matrix pair. The base pair is permutation-equivalent to the known SPC(3) product CSS code, and the present affine-coset description gives a direct proof that both Tanner graphs are $(3,8)$-regular with girth $8$. The base code has parameters $[[512,174,8]]$. We then apply circulant permutation matrix (CPM) lifts. The main decoding experiment uses the CPM-lifted code with lift factor $P=32$, which has parameters $[[16384,4142,18\le d\le32]]$, under the code-capacity depolarizing model. Complete enumeration excludes nontrivial logical representatives of weight at most $16$ on both sides, giving $d_X,d_Z\ge18$. Explicit weight-$32$ non-stabilizer logical representatives give the upper bounds $d_X,d_Z\le32$. A belief-propagation decoder with post-processing achieved frame error rate about $10^{-8}$ at $p=0.085$; an independently observed logical residual of weight $40$ is consistent with the sharper structural bound.

We study the ground-state phase diagram of a spin-1 Heisenberg chain with staggered long-range (LR) interactions decaying as $\propto r^{-\alpha}$ using a quantum Monte Carlo approach based on the split-spin representation. This formulation enables efficient large-scale simulations by mapping the spin-1 model onto spin-$1/2$ degrees of freedom with local projection constraints. We resolve the continuous quantum phase transition between the gapped Haldane phase at large $\alpha$ (short-range regime) and a gapless antiferromagnetically ordered Néel phase at small $\alpha$ (LR regime), where the continuous SU(2) symmetry is broken. From finite-size scaling and crossing point analyses, we determine the critical point to be at $\alpha_c = 2.49(1)$ and extract the associated critical exponents, which indicate unconventional criticality. In particular, the transition is found to be nonconformal, characterized by a dynamical exponent $z \neq 1$. We further analyze the scaling of entanglement entropy and bipartite fluctuations across the transition, and determine the corresponding universal scalings in both phases and at criticality.

Measurement-induced phase transitions (MIPTs) in monitored quantum systems are typically diagnosed using entanglement-based measures. Here, we develop a complementary thermodynamic perspective based on the arrow of time (AoT), which arises from the intrinsic irreversibility of the quantum measurements driving these transitions. We study the AoT - defined as the logarithmic ratio of forward and backward trajectory probabilities - across a family of models exhibiting MIPTs. We find that, like entanglement entropy, the AoT is a nonlinear functional of the averaged density matrix; however, in contrast to entanglement, it is associated with a local operator. To determine whether the AoT exhibits critical behavior, we formulate and exactly solve a model of a random quantum circuit with non-projective measurements. This allows us to analytically demonstrate that the AoT displays nonanalytic behavior and identify its critical exponent. Our results establish the AoT as a novel diagnostic for phase transitions in monitored quantum systems.

We present a method for learning quantum hardware noise from a measurement distribution of a single device experiment. Each noise channel is represented by automatically differentiable Kraus operators obtained from a Stinespring-based parameterization that is completely positive and trace preserving by construction, and circuits are simulated with a matrix product density operator forward model. Independent channels are attached to each native gate type, to each nearest-neighbor crosstalk interaction, and to state preparation and measurement, and all channels are optimized end-to-end against a distance between the simulated and observed measurement distributions. On ibm_fez, a Heron-generation superconducting processor, training on a ripple-carry adder circuit reproduces the device output distribution, and the same learned parameters, applied without retraining, also track the device distribution of an unrelated multiplier circuit, indicating that the method captures intrinsic device characteristics rather than overfitting to the training circuit. A systematic evaluation across a range of benchmark circuits confirms that this generalization is consistent. We further use the learned model to perform an offline feasibility assessment of the quantum approximate optimization algorithm with an error detection scheme, demonstrating the kind of noise-aware prediction the framework is designed to enable.

When a quantum electronic device is coupled to an electrical resonator, admittance changes of the quantum subsystem may be detected. The effective reactance may include capacitive and inductive terms that incorporate geometric, quantum, and tunneling components; while the effective resistance may be composed of Sisyphus and Hermes terms linked to relaxation and decoherence, respectively. Such reflectometry is usually studied when all characteristic times of the quantum system are much shorter than the resonator's period, in which case only stationary quantum states are probed. We present a rigorous description of a driven-dissipative qudit-resonator system. Our approach demonstrates how to strictly introduce quantum and tunneling capacitances as well as Hermes and Sisyphus resistances, and how these values are modified when the dynamics of the subsystems becomes mutually dependent. We present the cases of a Cooper-pair box, a single-Cooper-pair transistor, a double quantum dot, and a single-electron box. Our approach can be applied to describe any quantum system coupled to any classical resonator.

We develop a path integral representation for the dynamics of quantum systems with a finite-dimensional Hilbert space, formulated entirely within a discrete phase space. Starting from the discrete Wigner function on ${\mathbb{Z}_d} \times {\mathbb{Z}_d}$ ($d$ an odd prime) and the associated Weyl transform built from generalized displacement operators, we derive an exact kernel that propagates the discrete Wigner function in time and, by iterating its composition law through a short-time approximation, obtain a sum-over-paths expression weighted by a discrete phase-space action – the natural finite-dimensional counterpart of Marinov's functional. For Hamiltonians linear in the phase-space coordinates and at times strictly commensurate with the lattice, the fluctuation sum collapses to a deterministic shift, realizing the discrete analog of classical Hamiltonian flow. Applying the formalism to one and to two interacting qutrits ($d=3$), we show that the full entanglement dynamics – captured by a closed-form linear entropy valid for all times – requires the coherent contribution of all fluctuation sectors; the boundary-term (mean-field) sector alone fails to reproduce it. For a non-stabilizer Hamiltonian, where the short-time kernel is only approximate, the time-sliced path integral converges to the exact dynamics, including the dynamical generation of Wigner negativity. We discuss implications for the semiclassical simulation of many-body spin systems and for the characterization of non-classicality through Wigner negativity.

The eigenstates of two test-masses (where each test-mass is placed inside of a harmonic trap) separated by a distance, can get entangled where gravity acts as the mediator of entanglement and it has been argued in \href{https://doi.org/10.48550/arXiv.2511.07348}{arXiv:2511.07348 [quant-ph]} that this entanglement of masses cannot be generated without the underlying quantum nature of gravity. In this work, we consider two non-relativistic Bose-Einstein condensates (formed inside of harmonic trap potentials with identical trapping frequencies) separated by a distance. We take a linearized quantum gravity model and investigate the generation of entanglement while gravitons serve as the mediator of entanglement. The entanglement is generated between the phonon modes of the two condensates, and we observe that for very low separation distance, the entanglement generated is significantly higher than that observed for the quantum gravity induced entanglement of masses or QGEM protocol; however, the fall of entanglement is faster than the two-particle case for two separated Bose-Einstein condensates. We observe that when the number of particles in the condensate is increased, the degree of entanglement for a smaller separation distance becomes substantially higher compared to the case discussed in \href{https://doi.org/10.1103/PhysRevD.105.106028}{Phys. Rev. D 105 (2022) 106028}, which allows for a more robust experimental proposal using this quantum gravity induced entanglement of phonons or QGEP protocol.

The main claim of Baumeler et al. (Proc. R. Soc. A 481, 20240601; arXiv:2206.08467) is that "functional (deterministic) no-signalling resources can be stronger than probabilistic ones" in a certain nonlocal game on a Bell scenario with countably many parties. We disagree and argue that (i) under standard definitions, deterministic no-signalling resources are always probabilistic no-signalling resources; (ii) the deterministic strategy considered in Baumeler et al. can be promoted to a genuinely probabilistic strategy with similar properties; and (iii) a key step in a derivation in Baumeler et al., claimed to hold for all no-signalling strategies, implicitly assumes measurability, leaving a gap in the argument. We propose measurability assumptions which we conjecture would make this derivation rigorous. Taken together, the phenomenon highlighted by Baumeler et al. is best understood as a difference between measurable and non-measurable no-signalling resources.

Hydrogen-like ions cannot usefully bound the gyrogravitational ratio $\chi$: in ${}^{95}$Mo$^{41+}$ gravity shifts the hyperfine constants by fractions of $10^{-39}$–$10^{-38}$ with the tensor structure of the hyperfine interaction, centroids cancel it, and the $2p$ lines are 1.3 eV wide. Isotopes with spin can, however, enter the King plots that search for new bosons: for the cadmium clock line, the second-order hyperfine shift of ${}^{111}$Cd is $-0.96$~MHz, and a specific difference with ${}^{113}$Cd reduces it to about $2$~Hz after calculable anomaly and third-order corrections.

Sampling problems are promising candidates for demonstrating quantum advantage, and one approach known as quantum-enhanced Markov chain Monte Carlo [Layden, D. et al., Nature 619, 282-287 (2023)] uses quantum samples as a proposal distribution to accelerate convergence to a target distribution. On the other hand, many practical problems are large-scale and constrained, making it difficult to construct efficient proposal distributions in classical methods and slowing down MCMC mixing. In this work, we propose a divide-and-conquer neural network surrogate framework for quantum sampling to accelerate MCMC under fixed Hamming weight constraints. Our method divides the interaction graph for an Ising problem into subgraphs, generates samples using QAOA for those subproblems with an XY mixer, and trains neural network surrogates conditioned on the Hamming weight to provide proposal distributions for each subset while preserving the constraint. In numerical experiments of Boltzmann sampling on 3-regular graphs, our method consistently accelerated mixing as the system size $N$ increased, with average improvements in the autocorrelation decay rate constant by speedup factors of about $20.3$ and $7.6$ over classical pair-flip methods based on nearest-neighbor and non-nearest-neighbor exchanges, respectively. We also applied the method to an MNIST feature mask optimization problem with $N=784$, obtaining faster energy convergence and a $2.03\%$ higher classification accuracy. These results show that our method enables efficient and scalable MCMC and can outperform classical methods for practical applications on NISQ devices.

While the $^{229}$Th nuclear isomer has recently been observed and laser-excited, converting optical nuclear manipulation into a chip-scale solid-state frequency standard remains an open challenge. Here, we present a nanophotonic platform to realize an all-solid-state nuclear clock based on the low-energy isomeric transition of $^{229}$Th embedded in high-$Q$ fluoride photonic resonators. By coupling ensembles of thorium nuclei to confined optical modes, we show that resonant field build-up in the cavity can substantially enhance the nuclear excitation rate, enabling optical interrogation at practical laser intensities. We model the nuclei-photon interaction dynamics and outline a technological roadmap toward addressing this challenge, including resonator fabrication in fluoride crystals, thorium implantation, nuclear excitation with integrated lasers, and on-chip detection of vacuum-ultraviolet photons. As an initial proof of concept, we implant a crystalline fluoride whispering-gallery-mode resonator with $^{229}$Th and assess the impact of implantation-induced damage on resonator performance. Our platform leverages recent advances in materials integration and nanophotonics to chart a realistic route toward compact and scalable nuclear frequency standards.