2025-07-11 | | Total: 50
Violation of local realism via Bell inequality - a profound and counterintuitive manifestation of quantum theory that conflicts with the prediction of local realism - is viewed to be intimately linked with quantum entanglement. Experimental demonstrations of such a phenomenon using quantum entangled states are among the landmark experiments of modern physics and paved the way for quantum technology. Here we report the violation of the Bell inequality that cannot be described by quantum entanglement in the system but arises from quantum indistinguishability by path identity, shown by the multi-photon frustrated interference. By analyzing the measurement of four-photon frustrated interference within the standard Bell-test formalism, we find a violation of Bell inequality by more than four standard deviations. Our work establishes a connection between quantum correlation and quantum indistinguishability, providing insights into the fundamental origin of the counterintuitive characteristics observed in quantum physics.
Hybrid solvers for combinatorial optimization problems combine the advantages of classical and quantum computing to overcome difficult computational challenges. Although their theoretical performance seems promising, their practical applicability is challenging due to the lack of a technological stack that can seamlessly integrate quantum solutions with existing classical optimization frameworks. We tackle this challenge by introducing the ProvideQ toolbox, a software tool that enables users to easily adapt and configure hybrid solvers via Meta-Solver strategies. A Meta-Solver strategy implements decomposition techniques, which splits problems into classical and quantum subroutines. The ProvideQ toolbox enables the interactive creation of such decompositions via a Meta-Solver configuration tool. It combines well-established classical optimization techniques with quantum circuits that are seamlessly executable on multiple backends. This paper introduces the technical details of the ProvideQ toolbox, explains its architecture, and demonstrates possible applications for several real-world use cases. Our proof of concept shows that Meta-Solver strategies already enable the application of quantum subroutines today, however, more sophisticated hardware is required to make their performance competitive.
Computing the excitation spectra of quantum many-body systems on noisy quantum devices is a promising avenue to demonstrate the practical utility of current quantum processors, especially as we move toward the ``megaquop'' regime. For this task, here we introduce a \textit{tangent-space excitation ansatz} for quantum circuits, motivated by the quasi-particle picture of many-body systems and the structural similarity between quantum circuits and classical tensor networks. Increasing the circuit depth by one layer to construct the tangent space around the variational optimum of a parametrized quantum circuit, we show that a large number of low-energy states can be accurately captured. We demonstrate this ansatz using various models in both one and two spatial dimensions, including the challenging kagome Heisenberg antiferromagnet. We further provide evidence that this approach, implementable using Hadamard test, is stable in the presence of measurement noise and scalable to large system size.
We show that deep neural networks, including transformers and RNNs, pretrained as usual on next-token prediction, intrinsically discover and represent beliefs over 'quantum' and 'post-quantum' low-dimensional generative models of their training data -- as if performing iterative Bayesian updates over the latent state of this world model during inference as they observe more context. Notably, neural nets easily find these representation whereas there is no finite classical circuit that would do the job. The corresponding geometric relationships among neural activations induced by different input sequences are found to be largely independent of neural-network architecture. Each point in this geometry corresponds to a history-induced probability density over all possible futures, and the relative displacement of these points reflects the difference in mechanism and magnitude for how these distinct pasts affect the future.
A three-dimensional magneto-optical trap (MOT), as an efficient method of producing cold atoms from room-temperature atomic vapor, has been widely used to develop atomic sensors. Various compact MOTs using a single laser beam have been reported, simplifying apparatuses and leading to miniaturized devices. Here, we propose single-beam MOTs based on back-to-back pyramidal and conical mirrors. In such back-to-back mirrors, a MOT trapping volume is formed by an incident laser beam, a retroreflected beam, and multiple reflections from the mirror surfaces. We present the design of back-to-back mirrors and a series of compact MOT configurations, with the potential of increasing access to the MOT and simultaneously creating multiple MOTs. We demonstrate a MOT in a back-to-back conical mirror, loading 10 million rubidium-87 atoms from background vapor and cooling the atoms to 7 {\mu}K using polarization gradients. Single-beam MOTs based on back-to-back mirrors will contribute to building compact and scalable cold-atom-based sensors.
The experimental verification of multipartite entangled states is essential for advancing quantum information processing. Entanglement witnesses (EWs) provide a widely used and experimentally accessible approach for detecting genuinely multipartite entangled states. In this work, we theoretically derive the entanglement witness for the four-qubit Dicke state and experimentally evaluate it on two distinct IBM 127-qubit Quantum Processing Units (QPUs), namely ibm\_sherbrook and ibm\_brisbane. A negative expectation value of the witness operator serves as a sufficient condition for confirming genuine multipartite entanglement. We report the maximum (negative) values of the witness achieved on these QPUs as −0.178±0.009 and -0.169 \pm 0.002, corresponding to two different state preparation protocols. Additionally, we theoretically investigate the effect of various noise channels on the genuine entanglement of a four-qubit Dicke state using the Qiskit Aer simulator. We show the behavior of the EW constructed under the assumption of Markovian and non-Markovian amplitude damping and depolarizing noises, bit-phase flip noise, and readout errors. We also investigate the effect of varying thermal relaxation time on the EW, depicting a bound on the T_1 time required for successful generation of a Dicke State on a superconducting QPU.
Semiclassical approximations for various representations of a quantum state are constructed on a single (Lagrangian) surface in phase space, but it is not available for chaotic systems. An analogous evolution surface underlies semiclassical representations of the evolution operator, albeit in a doubled phase space. It is here shown that, corresponding to the Fourier transform on a unitary operator, represented as a Green function or spectral Wigner function, a Legendre transform generates a resolvent surface as the classical basis for semiclassical representations of the resolvent operator in double phase space, independently of the integrable or chaotic nature of the system. This surface coincides with derivatives of action functions (or generating functions) depending on the choice of appropriate coordinates and its growth departs from the energy shell following trajectories in double phase space. In an initial study of the resolvent surface based on its caustics, its complex nature is revealed to be analogous to a multidimensional sponge. Resummation of the trace of the resolvent in terms of linear combinations of periodic orbits, known as pseudo orbits or composite orbits, provides a cutoff to the semiclassical sum at the Heisenberg time. It is here shown that the corresponding actions for higher times can be approximately included within true secondary periodic orbits, in which multiple windings of short periodic orbits are joined by heteroclinic orbits into larger circuits.
The Sachdev-Ye-Kitaev (SYK) model describes a strongly correlated quantum system that shows a strong signature of quantum chaos. Due to its chaotic nature, the simulation of real-time dynamics becomes quickly intractable by means of classical numerics, and thus, quantum simulation is deemed to be an attractive alternative. Nevertheless, quantum simulations of the SYK model on noisy quantum processors are severely limited by the complexity of its Hamiltonian. In this work, we simulate the real-time dynamics of a sparsified version of the SYK model with 24 Majorana fermions on a trapped-ion quantum processor. We adopt a randomized quantum algorithm, TETRIS, and develop an error mitigation technique tailored to the algorithm. Leveraging the hardware's high-fidelity quantum operations and all-to-all connectivity of the qubits, we successfully calculate the Loschmidt amplitude for sufficiently long times so that its decay is observed. Based on the experimental and further numerical results, we assess the future possibility of larger-scale simulations of the SYK model by estimating the required quantum resources. Moreover, we present a scalable mirror-circuit benchmark based on the randomized SYK Hamiltonian and the TETRIS algorithm, which we argue provides a better estimate of the decay of fidelity for local observables than standard mirror-circuits.
The variational quantum eigensolver (VQE) is a promising method for simulating molecular systems on near-term quantum computers. This approach employs energy estimation; however, other relevant molecular properties can be extracted from the one-particle reduced density matrix (1-RDM) generated by VQE. The accuracy of these properties strongly depends on the reliability and convergence of the 1-RDM, which is not guaranteed by energy-only optimization. Thus, we investigate the effect of optimizing the 1-RDM within VQE to improve the accuracy of both the energy and molecular properties. A two-step algorithm was implemented that optimizes the energy and 1-RDM by incorporating a penalty term in the cost function to enforce the convergence of the 1-RDM. The first step focuses on energy minimization, while in the second step, a weighted penalty is added to the cost function to promote simultaneous improvement of the energy and 1-RDM. This approach was tested and validated for the k-UpCCGSD and GateFabric ansätzes with active spaces (4,4) and (2,2), respectively. k-UpCCGSD produces energies close to CISD, so optimizing 1-RDM has little effect on the energy but significantly improves electronic properties such as electron density, dipole moments, and atomic charges. GateFabric initially shows higher energy deviations from CISD, but optimizing 1-RDM substantially improves both the energy accuracy and the quality of 1-RDM. These results demonstrate that simultaneous optimization of energy and 1-RDM is an effective strategy to improve the accuracy of energies and molecular properties in variational quantum algorithms.
The Anderson model in one dimension is a quantum particle on a discrete chain of sites with nearest-neighbor hopping and random on-site potentials. It is a progenitor of many further models of disordered systems, and it has spurred numerous developments in diverse branches of physics. The literature is silent, however, on practical analytic tools for computing the density-of-states of this model when the distribution of the on-site potentials is arbitrary. Here, supersymmetry-based techniques are employed to give an explicit linear integral equation whose solutions control the density-of-states. The output of this analytic procedure is in perfect agreement with numerical sampling. By Thouless formula, these results immediately provide analytic control over the localization length.
Critical phenomena at finite temperature underpin a broad range of physical systems, yet their study remains challenging due to computational bottlenecks near phase transitions. Quantum annealers have attracted significant interest as a potential tool for accessing finite temperature criticality beyond classical reach, but their utility in precisely resolving criticality has remained limited by noise, hardware constraints, and thermal fluctuations. Here we overcome these challenges, showing that careful calibration and embedding allow quantum annealers to capture the full finite-temperature critical behavior of the paradigmatic two-dimensional Ising ferromagnet. By tuning the energy scale of the system and mitigating device asymmetries, we sample effective Boltzmann distributions and extract both the critical temperature and the associated critical exponents. Our approach opens the study of equilibrium and non-equilibrium critical phenomena in a broad class of systems at finite temperature.
A flagship application of quantum computers is the simulation of other quantum systems, including quantum field theories. In this article, we show how quantum computers can be employed to naturally calculate Feynman diagrams and their interferences in Quantum Chromodynamics (QCD). We simulate the colour parts of the interactions directly on the quantum computer, while the kinematic parts are for now pre-computed classically. For processes where some of the external particles are identical, we find the first hints of a potential quantum advantage. We validate our techniques using simulated quantum computers. Furthermore, for toy examples we also demonstrate our algorithms on a 56-qubit trapped-ion quantum computer. The work constitutes a further key step towards a full quantum simulation of generic perturbative QCD processes.
The Sachdev-Ye-Kitaev (SYK) model delivers an analytically tractable framework for exotic strongly correlated phases where conventional paradigms like Landau's Fermi liquid theory collapse. This review offers a pedagogical introduction to the SYK physics, highlighting its unique capacity to model \textit{strange metals} -- systems exhibiting linear-in-temperature resistivity, Planckian dissipation, and quasiparticle breakdown. We systematically construct both Majorana and complex fermion variants, transforming them into training grounds for modern many-body physics techniques, for instance, (1) large-N formulations via disorder averaging and replica symmetry, (2) Schwinger-Dyson and Kadanoff-Baym equations, (3) imaginary time Matsubara formulation, (4) real-time dynamics via Keldysh formalism, and the associated (5) non-perturbative Keldysh contour deformations. These tools lay the foundation for equilibrium thermodynamics, quantum chaos, quench dynamics, and transport in the thermodynamic limit -- all within a solvable, chaotic quantum system. Designed as a self-contained resource, the review bridges advanced technical machinery to physical insights, with computational implementations provided. Though principally treating the SYK model as a condensed matter laboratory, we also highlight its profound connection to quantum gravity -- woven throughout this work -- underscoring how this solvable chaotic fermionic model serves as a lens onto black hole thermodynamics and holographic duality.
We investigate the properties of the vacuum state for the Proca field in the geometry of two parallel plates on background of (D+1)-dimensional Minkowski spacetime. The two-point functions for the vector potential and the field tensor are evaluated for higher-dimensional generalizations of the perfect magnetic conductor (PMC) and perfect electric conductor (PEC) boundary conditions. Explicit expressions are provided for the vacuum expectation values (VEVs) of the electric and magnetic field squares, field condensate, and for the VEV of the energy-momentum tensor. In the zero-mass limit the VEVs of the electric and magnetic field squares and the condensate reduce to the corresponding expressions for a massless vector field. The same is the case for the VEV of the energy-momentum tensor in the problem with PEC conditions. However, for PMC conditions the zero-mass limit for the vacuum energy-momentum tensor differs from the corresponding VEV for a massless field. This difference in the zero-mass limits is related to the different influences of the boundary conditions on the longitudinal polarization mode of a massive vector field. The PMC conditions constrain all the polarization modes including the longitudinal mode, whereas PEC conditions do not influence the longitudinal mode. The vacuum energy-momentum tensor is diagonal. The normal stress is uniformly distributed in the region between the plates and vanishes in the remaining regions. The corresponding Casimir forces are attractive for both boundary conditions.
We introduce a generalized \textit{Probabilistic Approximate Optimization Algorithm (PAOA)}, a classical variational Monte Carlo framework that extends and formalizes prior work by Weitz \textit{et al.}~\cite{Combes_2023}, enabling parameterized and fast sampling on present-day Ising machines and probabilistic computers. PAOA operates by iteratively modifying the couplings of a network of binary stochastic units, guided by cost evaluations from independent samples. We establish a direct correspondence between derivative-free updates and the gradient of the full 2^N \times 2^N Markov flow, showing that PAOA admits a principled variational formulation. Simulated annealing emerges as a limiting case under constrained parameterizations, and we implement this regime on an FPGA-based probabilistic computer with on-chip annealing to solve large 3D spin-glass problems. Benchmarking PAOA against QAOA on the canonical 26-spin Sherrington-Kirkpatrick model with matched parameters reveals superior performance for PAOA. We show that PAOA naturally extends simulated annealing by optimizing multiple temperature profiles, leading to improved performance over SA on heavy-tailed problems such as SK-Lévy.
We investigated the quantum effects of a frustrated Ising model on a two-dimensional square lattice using the Quantum Approximate Optimization Algorithm (QAOA). While strong spin frustration is known to induce quantum fluctuations at low temperatures, previous classical approaches restricted to binary (up or down) spin configurations have been insufficient to fully capture the quantum contributions of frustration. In this study, we introduced a quantitative metric to evaluate the quantum effects arising from frustration and employed QAOA to differentiate between classical and quantum regimes. Notably, we found that in the weakly frustrated region, QAOA measurements rarely capture first excited states, as they are energetically well separated from the ground state. In contrast, near the quantum phase transition point, excited states appear more frequently in QAOA measurements, highlighting the increased role of quantum fluctuations.
We report on an optical magnetometer enhanced by vacuum-squeezed light, employing an Mx magnetometer based on ^{87}Rb vapor in a micrometer-scale cell (~100 \mum). Using the well-established polarization self-rotation effect in a room-temperature ^{87}Rb vapor cell, we achieve -3 dB of vacuum squeezing within the noise spectral window of 100 Hz to several MHz, corresponding to 3.5 dB squeezing when accounting for optical losses. Leveraging this level of squeezing, we demonstrate a magnetic field sensitivity of approx. 1 pT/\sqrt{Hz}. The combination of vacuum-squeezed light and micrometer-scale vapor cells paves the way for compact, low-power-consumption atomic sensors with enhanced performance.
We develop a general framework based on the functional derivative to extract nonlinear dynamical response functions from the temporal evolution of physical quantities, without explicitly computing multipoint correlation functions. We validate our approach by calculating the second- and third-order optical responses in the Rice-Mele model and further apply it to a many-body interacting system using a tensor network method. This framework is broadly applicable to any method that can compute real-time dynamics, offering a powerful and versatile tool for investigating nonlinear responses in dynamical systems.
Recently, a gauge theory of unified gravity [Rep. Prog. Phys. 88, 057802 (2025)] has been developed to extend the Standard Model to include gravity. Here we present unified gravity using the ordinary four-vector and tensor field notation of the Standard Model. The main goal of the present work is to extend the original Minkowski spacetime formulation of the theory to account for graviton-graviton interaction. This is a necessary extension for problems involving interactions between gravitational fields, for example, in the propagation of gravitational waves in external gravitational potentials. The 4xU(1) gauge invariance of unified gravity is preserved in this extension.
Sensing with undetected photons has become a vibrant, application-driven research domain with a special focus on the mid-infrared (mid-IR) wavelength region. Since the mid-IR contains spectral bands with highly specific and strong molecular absorbance signatures, often referred to as fingerprints, a multitude of different samples and their compositions can be detected and quantified spectroscopically. Enhancing this inherently sample alteration-free spectroscopic method with imaging capabilities leads to a powerful technique for environmental monitoring and biomedical applications that enables automated diagnostics while omitting time-consuming and non-reversible labeling steps. To evade the shortcomings of state-of-the-art instruments for mid-IR hyperspectral microscopy related to cost, complexity, power-consumption, and performance, which are associated with technological challenges for mid-IR cameras and low-noise and broadband mid-IR sources, here, we construct a proof-of-concept nonlinear interferometer in a wide-field imaging arrangement combined with high-resolution spectral acquisition by pixelwise quantum Fourier transform infrared spectroscopy. For the broadband range of 2300-3100 cm^{-1}, covering the important CH-stretch band, we perform hyperspectral imaging that simultaneously resolves 3500 spatial modes, each with a spectral resolution of 10 cm^{-1} leveraging in total around 10^5 spatio-spectrally entangled photon modes. Our image acquisition uses a commercial, megapixel sCMOS camera, while a medium-power and compact c.w. pump laser is the only necessary light source.
We present a general theorem demonstrating non-perturbatively the decrease of the "effective central charge" c_{\text{eff}}=(d c/dR)|_{R=1} under renormalization group (RG) flow in the R\rightarrow1 replica limit of a R-copy 2D conformal field theory (CFT) action S_{*} perturbed by a replica interaction of the form -\mathbb{S}=-\sum_{a=1}^{R}S_{*}^{(a)}+\Delta\int d^2 x \sum_{\substack{a,b=1\\ a\neq b}}^{R}\varphi^{(a)}(x)\varphi^{(b)}(x). Here \varphi is a scaling field belonging to the CFT with action S_* and the coupling \Delta is relevant in the RG sense. We show that the infrared value of c_{\text{eff}} is always \textit{less} than the central charge c of the unperturbed CFT S_{*}. We refer to this result as the "c-effective theorem". As an application of this theorem, we consider replica field theories in the limit of R \to 1 replicas of the form above, shown by Nahum and Jacobsen [arXiv:2504.01264] to describe 2D classical monitored systems, where measurements introduce a form of quenched randomness via Bayes' theorem. Lastly, we discuss a possible relationship of our theorem with the effective central charge c_{\text{eff}}^{(R\rightarrow0)}=(dc/dR)|_{R=0} for the above replica action in the different R\rightarrow0 replica limit, which is of relevance to systems with generic uncorrelated impurity-type quenched disorder, as opposed to measurements.
The Hepp-Lieb-Dicke model is ubiquitous in cavity quantum electrodynamics, describing spin-cavity coupling which does not conserve excitation number. Coupling the closed spin-cavity system to an environment realizes the open Dicke model, and by tuning the structure of the environment or the system-environment coupling, interesting spin-only models can be engineered. In this work, we focus on a variation of the multi-species open Dicke model which realizes mediated nonreciprocal interactions between the spin species and, consequently, an interesting dynamical limit-cycle phase. In particular, we improve upon adiabatic elimination and, instead, employ a Redfield master equation in order to describe the effective dynamics of the spin-only system. We assess this approach at the mean-field level, comparing it both to adiabatic elimination and the full spin-cavity model, and find that the predictions are sensitive to the presence of single-particle incoherent decay. Additionally, we clarify the symmetries of the model and explore the dynamical limit-cycle phase in the case of explicit parity-time-symmetry breaking, finding a region of phase coexistence terminating at an codimension-two exceptional point. Lastly, we go beyond mean-field theory by exact numerical diagonalization of the master equation, appealing to permutation symmetry in order to increase the size of accessible systems. We find signatures of phase transitions even for small system sizes.
We give an introduction to the theory of quantum error correction using stabilizer codes that is geared towards the working computer scientists and mathematicians with an interest in exploring this area. To this end, we begin with an introduction to basic quantum computation for the uninitiated. We then construct several examples of simple error correcting codes without reference to the underlying mathematical formalism in order to develop the readers intuition for the structure of a generic code. With this in hand, we then discuss the more general theory of stabilizer codes and provide the necessary level of mathematical detail for the non-mathematician. Finally, we give a brief look at the elegant homological algebra formulation for topological codes. As a bonus, we give implementations of the codes we mention using OpenQASM, and we address the more recent approaches to decoding using neural networks. We do not attempt to give a complete overview of the entire field, but provide the reader with the level of detail needed to continue in this direction.
Monitored quantum systems, where unitary dynamics compete with continuous measurements, exhibit dynamical transitions as the measurement rate is varied. These reflect abrupt changes in the structure of the evolving wavefunction, captured by complementary complexity diagnostics that include and go beyond entanglement aspects. Here, we investigate how monitoring affects magic state resources, the nonstabilizerness, of Gaussian fermionic systems. Using scalable Majorana sampling techniques, we track the evolution of stabilizer Rényi entropies in large systems under projective measurements. While the leading extensive (volume-law) scaling of magic remains robust across all measurement rates, we uncover a sharp transition in the subleading logarithmic corrections. This measurement-induced complexity transition, invisible to standard entanglement probes, highlights the power of magic-based diagnostics in revealing hidden features of monitored many-body dynamics.
The recent topological classification of non-Hermitian `Hamiltonians' is usually interpreted in terms of pure quantum states that decay or grow with time. However, many-body systems with loss and gain are typically better described by mixed-state open quantum dynamics, which only correspond to pure-state non-Hermitian dynamics upon a postselection of measurement outcomes. Since postselection becomes exponentially costly with particle number, we here investigate to what extent the most important example of non-Hermitian topology can survive without it: the non-Hermitian skin effect and its relationship to a bulk winding number in one spatial dimension. After defining the winding number of the Lindbladian superoperator for a quadratic fermion system, we systematically relate it to the winding number of the associated postselected non-Hermitian Hamiltonian. We prove that the two winding numbers are equal (opposite) in the absence of gain (loss), and provide a physical explanation for this relationship. When both loss and gain are present, the Lindbladian winding number typically remains quantized and non-zero, though it can change sign at a phase transition separating the loss and gain-dominated regimes. This transition, which leads to a reversal of the Lindbladian skin effect localization, is rendered invisible by postselection. We also identify a case where removing postselection induces a skin effect from otherwise topologically trivial non-Hermitian dynamics.