Search

We present a comprehensive numerical investigation of the cluster Truncated Wigner Approximation (cTWA) applied to quench dynamics in bond-disordered Heisenberg spin chains with power-law interactions. We find that cTWA yields highly accurate results over a wide parameter range. However, its accuracy hinges on a suitable choice of clusters. By using a clustering strategy inspired by the strong disorder renormalisation group (SDRG)/real-space renormalization group (RSRG), clusters of two spins are sufficient to obtain essentially exact results in the regime of strong disorder. Surprisingly, even for rather weak disorder, e.g.\ in the presence of very long-range interactions, this choice of clustering outperforms a naive choice of clusters of consecutive spins. Additionally, we develop a discrete sampling scheme for the initial Wigner function, as an alternative to the originally introduced scheme based on Gaussian approximations. This sampling scheme puts cTWA on the same conceptional footing as regular dTWA for single spins and yields some reduction in the Monte Carlo shot noise compared to the Gaussian scheme., Comment: 16 pages, 13 figures, submitted to Phys. Rev. B

We investigate the information extractable from measurement distributions of two non-commuting spin observables in a multi-well spin-1 Bose-Einstein condensate. We provide a variety of analytic and numerical evidence that suitably chosen classical entropies and classical mutual informations thereof contain the typical feature of quantum entropies known in quantum field theories, that is, the area law, even in the non-Gaussian regime and for a non-zero temperature. Towards a feasible experimental implementation, we estimate entropic quantities from a finite number of samples without any additional assumptions on the underlying quantum state using k-nearest neighbor estimators., Comment: 25 pages; comments welcome!

The scaling of local quantum entropies is of utmost interest for characterizing quantum fields, many-body systems, and gravity. Despite their importance, theoretically and experimentally accessing quantum entropies is challenging as they are nonlinear functionals of the underlying quantum state. Here, we show that suitably chosen classical entropies capture the very same features as their quantum analogs for an experimentally relevant setting. We describe the post-quench dynamics of a multi-well spin-1 Bose-Einstein condensate from an initial product state via measurement distributions of spin observables and estimate the corresponding entropies using the asymptotically unbiased k-nearest neighbor method. We observe the dynamical build-up of quantum correlations signaled by an area law, as well as local thermalization revealed by a transition to a volume law, both in regimes characterized by non-Gaussian distributions. We emphasize that all relevant features can be observed at small sample numbers without assuming a specific functional form of the distributions, rendering our method directly applicable to a large variety of models and experimental platforms., Comment: 5+2 pages; comments welcome!

Time reversal in a macroscopic system is contradicting daily experience. It is practically impossible to restore a shattered cup to its original state by just time reversing the microscopic dynamics that led to its breakage. Yet, with the precise control capabilities provided by modern quantum technology, the unitary evolution of a quantum system can be reversed in time. Here, we implement a time-reversal protocol in a dipolar interacting, isolated many-body spin system represented by Rydberg states in an atomic gas. By changing the states encoding the spin, we flip the sign of the interaction Hamiltonian, and demonstrate the reversal of the relaxation dynamics of the magnetization by letting a demagnetized many-body state evolve back-in-time into a magnetized state. We elucidate the role of atomic motion using the concept of a Loschmidt echo. Finally, by combining the approach with Floquet engineering, we demonstrate time reversal for a large family of spin models with different symmetries. Our method of state transfer is applicable across a wide range of quantum simulation platforms and has applications far beyond quantum many-body physics, reaching from quantum-enhanced sensing to quantum information scrambling.

We investigate a disorder-free quantum Ising chain subject to a time-periodic drive that rotates each spin by an angle $\pi(1-\epsilon_i)$. In case all spins experience the same deviation $\epsilon$ and the system is initialized in a fully polarized state, the dynamics is known to be time-crystalline: the magnetization of the system exhibits period-doubled oscillations for timescales that grow exponentially with the length of the chain. In this work, we study the effect of a deviation $\epsilon$ that differs between spins. We find that reducing $\epsilon$ for a single spin drastically enhances the lifetime of spatio-temporal order, suggesting the name ``metronome" spin. Employing perturbative arguments in an average Hamiltonian picture, we explain this observation for initial states with macroscopic bulk magnetization. Furthermore, in the case of random bitstring initial states, we report the enhancement of the lifetime of a topological edge mode, which can also be understood in the same picture. Finally, we discuss an altered geometry in which the metronome spin is not directly part of the chain, affecting the dynamics in different ways in the two scenarios considered. Our findings unveil the intricate dynamics that emerge in Floquet systems under the influence of a spatially varying drive, thereby uncovering new avenues for Floquet engineering.

Quantum technologies rely heavily on accurate control and reliable readout of quantum systems. Current experiments are limited by numerous sources of noise that can only be partially captured by simple analytical models and additional characterization of the noise sources is required. We test the ability of readout error mitigation to correct realistic noise found in systems composed of quantum two-level objects (qubits). To probe the limit of such methods, we designed a beyond-classical readout error mitigation protocol based on quantum state tomography (QST), which estimates the density matrix of a quantum system, and quantum detector tomography (QDT), which characterizes the measurement procedure. By treating readout error mitigation in the context of state tomography the method becomes largely readout mode-, architecture-, noise source-, and quantum state-independent. We implement this method on a superconducting qubit and evaluate the increase in reconstruction fidelity for QST. We characterize the performance of the method by varying important noise sources, such as suboptimal readout signal amplification, insufficient resonator photon population, off-resonant qubit drive, and effectively shortened $T_1$ and $T_2$ decay times. As a result, we identified noise sources for which readout error mitigation worked well, and observed decreases in readout infidelity by a factor of up to 30., Comment: 18 pages, 12 figures, corrections, title change and added additional information

We explore a supervised machine learning approach to estimate the entanglement entropy of multi-qubit systems from few experimental samples. We put a particular focus on estimating both aleatoric and epistemic uncertainty of the network's estimate and benchmark against the best known conventional estimation algorithms. For states that are contained in the training distribution, we observe convergence in a regime of sample sizes in which the baseline method fails to give correct estimates, while extrapolation only seems possible for regions close to the training regime. As a further application of our method, highly relevant for quantum simulation experiments, we estimate the quantum mutual information for non-unitary evolution by training our model on different noise strengths., Comment: 6 pages, 4 figures

This paper introduces assignment flows for density matrices as state spaces for representing and analyzing data associated with vertices of an underlying weighted graph. Determining an assignment flow by geometric integration of the defining dynamical system causes an interaction of the non-commuting states across the graph, and the assignment of a pure (rank-one) state to each vertex after convergence. Adopting the Riemannian Bogoliubov-Kubo-Mori metric from information geometry leads to closed-form local expressions which can be computed efficiently and implemented in a fine-grained parallel manner. Restriction to the submanifold of commuting density matrices recovers the assignment flows for categorial probability distributions, which merely assign labels from a finite set to each data point. As shown for these flows in our prior work, the novel class of quantum state assignment flows can also be characterized as Riemannian gradient flows with respect to a non-local non-convex potential, after proper reparametrization and under mild conditions on the underlying weight function. This weight function generates the parameters of the layers of a neural network, corresponding to and generated by each step of the geometric integration scheme. Numerical results indicates and illustrate the potential of the novel approach for data representation and analysis, including the representation of correlations of data across the graph by entanglement and tensorization.

Quantum entanglement has been identified as a crucial concept underlying many intriguing phenomena in condensed matter systems, such as topological phases or many-body localization. Recently, instead of considering mere quantifiers of entanglement like entanglement entropy, the study of entanglement structure in terms of the entanglement spectrum has shifted into the focus, leading to new insights into fractional quantum Hall states and topological insulators, among others. What remains a challenge is the experimental detection of such fine-grained properties of quantum systems. The development of protocols for detecting features of the entanglement spectrum in cold-atom systems, which are one of the leading platforms for quantum simulation, is thus highly desirable and will open up new avenues for experimentally exploring quantum many-body physics. Here, we present a method to bound the width of the entanglement spectrum, or entanglement dimension, of cold atoms in lattice geometries, requiring only measurements in two experimentally accessible bases and utilizing ballistic time-of-flight (TOF) expansion. Building on previous proposals for entanglement certification for photon pairs, we first consider entanglement between two atoms of different atomic species and later generalize to higher numbers of atoms per species and multispecies configurations showing multipartite high-dimensional entanglement. Through numerical simulations, we show that our method is robust against typical experimental noise effects and thus will enable high-dimensional entanglement certification in systems of up to eight atoms using currently available experimental techniques., Comment: 20 pages + 8 pages Appendix

Neural quantum states are a new family of variational ans\"atze for quantum-many body wave functions with advantageous properties in the notoriously challenging case of two spatial dimensions. Since their introduction a wide variety of different network architectures has been employed to study paradigmatic models in quantum many-body physics with a particular focus on quantum spin models. Nonetheless, many questions remain about the effect that the choice of architecture has on the performance on a given task. In this work, we present a unified comparison of a selection of popular network architectures and symmetrization schemes employed for ground state searches of prototypical spin Hamiltonians, namely the two-dimensional transverse-field Ising model and the J1-J2 model. In the presence of a non-trivial sign structure of the ground states, we find that the details of symmetrization crucially influence the performance. We describe this effect in detail and discuss its consequences, especially for autoregressive models, as their direct sampling procedure is not compatible with the symmetrization procedure that we found to be optimal., Comment: 7 + 3 pages, 3 figures; Fixed references and a typo in table 1

We explore the advantages of a class of entanglement criteria for continuous variable systems based on the Husimi $Q$-distribution in scenarios with sparse experimental data. The generality of these criteria allows optimizing them for a given entangled state and experimental setting. We consider the scenario of coarse grained measurements, or finite detector resolution, where the values of the Husimi $Q$-distribution are only known on a grid of points in phase space, and show how the entanglement criteria can be adapted to this case. Further, we examine the scenario where experimental measurements amount to drawing independent samples from the Husimi distribution. Here, we customize our entanglement criteria to maximize the statistical significance of the detection for a given finite number of samples. In both scenarios optimization leads to clear improvements enlarging the class of detected states and the signal-to-noise ratio of the detection, respectively., Comment: Has been merged with arXiv:2211.17165

We prove a general class of continuous variable entanglement criteria based on the Husimi $Q$-distribution, which represents a quantum state in canonical phase space, by employing a theorem by Lieb and Solovej. We discuss their generality, which roots in the possibility to optimize over the set of concave functions, from the perspective of continuous majorization theory and show that with this approach families of entropic as well as second moment criteria follow as special cases. All derived criteria are compared to corresponding marginal based criteria and the strength of the phase space approach is demonstrated for a family of prototypical example states where only our criteria flag entanglement. Further, we explore their optimization prospects in two experimentally relevant scenarios characterized by sparse data: finite detector resolution and finite statistics. In both scenarios optimization leads to clear improvements enlarging the class of detected states and the signal-to-noise ratio of the detection, respectively., Comment: 23 pages, 11 figures, published version

We present a general class of entanglement criteria for continuous variable systems. Our criteria are based on the Husimi $Q$-distribution and allow for optimization over the set of all concave functions rendering them extremely general and versatile. We show that several entropic criteria and second moment criteria are obtained as special cases. Our criteria reveal entanglement of families of states undetected by any commonly used criteria and provide clear advantages under typical experimental constraints such as finite detector resolution and measurement statistics., Comment: 5 pages, 4 figures, published version

An important step towards a comprehensive understanding of far-from-equilibrium dynamics of quantum many-body systems is the identification of unifying features that are independent of microscopic details of the system. We experimentally observe such robust features in the magnetization relaxation dynamics of disordered Heisenberg XX-, XXZ- and Ising Hamiltonians. We realize these Heisenberg spin models with tunable anisotropy parameter and power-law interactions in an ensemble of Rydberg atoms by encoding the spin in suitable Rydberg state combinations. We consistently observe stretched-exponential relaxation of magnetization for all considered spin models, collapsing onto a single curve after appropriate rescaling of time. This robust short-time relaxation behavior is explained by a perturbative treatment that exploits the strong disorder in pairwise couplings, which leads to a description in terms of approximately independent pairs of spins. In numerical simulations of small systems, we show that these pairs of spins constitute approximate local integrals of motion, which remain at least partially conserved on a timescale exceeding the duration of the relaxation dynamics of the magnetization.

Strongly interacting quantum systems subject to quenched disorder exhibit intriguing phenomena such as glassiness and many-body localization. Theoretical studies have mainly focused on disorder in the form of random potentials, while many experimental realizations naturally feature disorder in the interparticle interactions. Inspired by cold Rydberg gases, where such disorder can be engineered using the dipole blockade effect,we study a Heisenberg XXZ spin model where the disorder is exclusively due to random spin-spin couplings, arising from power-law interactions between randomly positioned spins. Using established spectral and eigenstate properties and entanglement entropy, we show that this system exhibits a localization crossover and identify strongly interacting pairs as emergent local conserved quantities in the system, leading to an intuitive physical picture consistent with our numerical results., Comment: 10 pages, 6 main figures, 1 supplementary figure, close to published version

Understanding how closed quantum systems dynamically approach thermal equilibrium presents a major unresolved problem in statistical physics. Generically, non-integrable quantum systems are expected to thermalize as they comply with the Eigenstate Thermalization Hypothesis. However, in the presence of strong disorder, the dynamics can possibly slow down to a degree that systems fail to thermalize on experimentally accessible timescales, as in spin glasses or many-body localized systems. In general, particularly in long-range interacting quantum systems, the specific nature of the disorder necessary for the emergence of a prethermal, metastable state--distinctly separating the timescales of initial relaxation and subsequent slow thermalization--remains an open question. We study an ensemble of Heisenberg spins with a tunable distribution of random coupling strengths realized by a Rydberg quantum simulator. We observe a drastic change in the late-time magnetization when increasing disorder strength. The data is well described by models based on pairs of strongly interacting spins, which are treated as thermal for weak disorder and isolated for strong disorder. Our results indicate a crossover into a pair-localized prethermal regime in a closed quantum system of thousands of spins in the critical case where the exponent of the power law interaction matches the spatial dimension.

The accurate numerical solution of partial differential equations is a central task in numerical analysis allowing to model a wide range of natural phenomena by employing specialized solvers depending on the scenario of application. Here, we develop a variational approach for solving partial differential equations governing the evolution of high dimensional probability distributions. Our approach naturally works on the unbounded continuous domain and encodes the full probability density function through its variational parameters, which are adapted dynamically during the evolution to optimally reflect the dynamics of the density. For the considered benchmark cases we observe excellent agreement with numerical solutions as well as analytical solutions in regimes inaccessible to traditional computational approaches., Comment: 6 + 3 pages, 4 figures

Recent research has demonstrated the usefulness of neural networks as variational ansatz functions for quantum many-body states. However, high-dimensional sampling spaces and transient autocorrelations confront these approaches with a challenging computational bottleneck. Compared to conventional neural networks, physical-model devices offer a fast, efficient and inherently parallel substrate capable of related forms of Markov chain Monte Carlo sampling. Here, we demonstrate the ability of a neuromorphic chip to represent the ground states of quantum spin models by variational energy minimization. We develop a training algorithm and apply it to the transverse field Ising model, showing good performance at moderate system sizes ($N\leq 10$). A systematic hyperparameter study shows that scalability to larger system sizes mainly depends on sample quality, which is limited by temporal parameter variations on the analog neuromorphic chip. Our work thus provides an important step towards harnessing the capabilities of neuromorphic hardware for tackling the curse of dimensionality in quantum many-body problems., Comment: 14 pages, 7 figures

Modern day quantum simulators can prepare a wide variety of quantum states but the accurate estimation of observables from tomographic measurement data often poses a challenge. We tackle this problem by developing a quantum state tomography scheme which relies on approximating the probability distribution over the outcomes of an informationally complete measurement in a variational manifold represented by a convolutional neural network. We show an excellent representability of prototypical ground- and steady states with this ansatz using a number of variational parameters that scales polynomially in system size. This compressed representation allows us to reconstruct states with high classical fidelities outperforming standard methods such as maximum likelihood estimation. Furthermore, it achieves a reduction of the estimation error of observables by up to an order of magnitude compared to their direct estimation from experimental data., Comment: 11 pages, 6 figures

We numerically study out-of-equilibrium dynamics in a family of Heisenberg models with $1/r^6$ power-law interactions and positional disorder. Using the semi-classical discrete truncated Wigner approximation (dTWA) method, we investigate the time evolution of the magnetization and ensemble-averaged single-spin purity for a strongly disordered system after initializing the system in an out-of-equilibrium state. We find that both quantities display robust glassy behavior for almost any value of the anisotropy parameter of the Heisenberg Hamiltonian. Furthermore, a systematic analysis allows us to quantitatively show that, for all the scenarios considered, the stretch power lies close to the one analytically obtained in the Ising limit. This indicates that glassy relaxation behavior occurs widely in disordered quantum spin systems, independent of the particular symmetries and integrability of the Hamiltonian.

Path integrals with complex actions are encountered for many physical systems ranging from spin- or mass-imbalanced atomic gases and graphene to quantum chromo-dynamics at finite density to the non-equilibrium evolution of quantum systems. Many computational approaches have been developed for tackling the sign problem emerging for complex actions. Among these, complex Langevin dynamics has the appeal of general applicability. One of its key challenges is the potential convergence of the dynamics to unphysical fixed points. The statistical sampling process at such a fixed point is not based on the physical action and hence leads to wrong predictions. Moreover, its unphysical nature is hard to detect due to the implicit nature of the process. In the present work we set up a general approach based on a Markov chain Monte Carlo scheme in an extended state space. In this approach we derive an explicit real sampling process for generalized complex Langevin dynamics. Subject to a set of constraints, this sampling process is the physical one. These constraints originate from the detailed-balance equations satisfied by the Monte Carlo scheme. This allows us to re-derive complex Langevin dynamics from a new perspective and establishes a framework for the explicit construction of new sampling schemes for complex actions., Comment: 29 pages, 9 figures

We derive entropic inseparability criteria for the phase space representation of quantum states. In contrast to criteria involving differential entropies of marginal phase space distributions, our criteria are based on a joint distribution known as the Husimi Q-distribution. This distribution is experimentally accessible in cold atoms, circuit QED architectures and photonic systems and bears practical advantages compared to the detection of marginals. We exemplify the strengths of our entropic approach by considering several classes of non-Gaussian states where second-order criteria fail. We show that our criteria certify entanglement in previously undetectable regions highlighting the strength of using the Husimi Q-distribution for entanglement detection., Comment: 7 pages, 2 figures

We study excitation transport in a two-dimensional system of randomly assembled spins with power-law hopping in two dimensions. This model can be realized in cold atom quantum simulators with Rydberg atoms. In these experiments, due to the Rydberg blockade effect, the degree of disorder in the system is effectively tunable by varying the spin density. We study dynamics and eigenstate properties of the model as a function of disorder strength and system size and discuss potential limitations for experiments. At strong disorder we predominantly observe localized eigenstates with power-law tails. In this regime the spectral and eigenstate properties can be understood in a perturbative picture of states localized on small clusters of spins. As the disorder strength is weakened eigenstates become increasingly delocalized and a set of seemingly multifractal states appears in the low-energy tail of the spectrum. A detailed study of the system-size scaling of the eigenstate properties indicates that in the infinite-size limit all states eventually become localized. We discuss the feasibility of observing localization effects experimentally in the spatial spreading of an initially localized excitation and identify limited system sizes and finite decoherence rates as major challenges. Our study paves the way towards an experimental observation of localization effects in Rydberg spin systems with tunable disorder., Comment: v2 close to journal version

A prerequisite for the comprehensive understanding of many-body quantum systems is a characterization in terms of their entanglement structure. The experimental detection of entanglement in spatially extended many-body systems describable by quantum fields still presents a major challenge. We develop a general scheme for certifying entanglement and demonstrate it by revealing entanglement between distinct subsystems of a spinor Bose-Einstein condensate. Our scheme builds on the spatially resolved simultaneous detection of the quantum field in two conjugate observables which allows the experimental confirmation of quantum correlations between local as well as non-local partitions of the system. The detection of squeezing in Bogoliubov modes in a multi-mode setting illustrates its potential to boost the capabilities of quantum simulations to study entanglement in spatially extended many-body systems., Comment: 11 pages, 4 figures, 3 supplementary figures

We study the out-of-equilibrium dynamics in the quantum Ising model with power-law interactions and positional disorder. For arbitrary dimension $d$ and interaction range $\alpha \geq d$ we analytically find a stretched exponential decay of the global magnetization and ensemble-averaged single-spin purity with a stretch-power $\beta = d/\alpha$ in the thermodynamic limit. Numerically, we confirm that glassy behavior persists for finite system sizes and sufficiently strong disorder. We identify dephasing between disordered coherent pairs as the main mechanism leading to a relaxation of global magnetization, whereas genuine many-body interactions lead to a loss of single-spin purity which signifies the build-up of entanglement. The emergence of glassy dynamics in the quantum Ising model extends prior findings in classical and open quantum systems, where the stretched exponential law is explained by a scale-invariant distribution of time scales, to both integrable and non-integrable quantum systems.

We develop a variational approach to simulating the dynamics of open quantum many-body systems using deep autoregressive neural networks. The parameters of a compressed representation of a mixed quantum state are adapted dynamically according to the Lindblad master equation by employing a time-dependent variational principle. We illustrate our approach by solving the dissipative quantum Heisenberg model in one and two dimensions for up to 40 spins and by applying it to the simulation of confinement dynamics in the presence of dissipation., Comment: 7 + 5 pages, 3 figures

We experimentally investigate the nonlinear transmission spectrum of coherent light fields propagating through a Rydberg EIT medium with strong atomic interactions. In contrast to previous investigations, which have largely focused on resonant control fields, we explore here the full two-dimensional spectral response of the Rydberg gas. Our measurements confirm previously observed spectral features for a vanishing control-field detuning, but also reveal significant differences on two-photon resonance. In particular, we find qualitative deficiencies of mean-field models and rate-equation simulations as well as a third-order nonlinear susceptibility that accounts for pair-wise interaction effects at low probe-field intensities in describing the nonlinear probe-field response under EIT conditions. Our results suggest that a more complete understanding of Rydberg-EIT and emerging photon interactions requires to go beyond existing simplified models as well as few-photon theories.

We study experimentally accessible lower bounds on entanglement measures based on entropic uncertainty relations. Experimentally quantifying entanglement is highly desired for applications of quantum simulation experiments to fundamental questions, e.g., in quantum statistical mechanics and condensed-matter physics. At the same time it poses a significant challenge because the evaluation of entanglement measures typically requires the full reconstruction of the quantum state, which is extremely costly in terms of measurement statistics. We derive an improved entanglement bound for bipartite systems, which requires measuring joint probability distributions in only two different measurement settings per subsystem, and demonstrate its power by applying it to currently operational experimental setups for quantum simulation with cold atoms. Examining the tightness of the derived entanglement bound, we find that the set of pure states for which our relation is tight is strongly restricted. We show that for measurements in mutually unbiased bases the only pure states that saturate the bound are maximally entangled states on a subspace of the bipartite Hilbert space (this includes product states). We further show that our relation can also be employed for entanglement detection using generalized measurements, i.e., when not all measurement outcomes can be resolved individually by the detector. In addition, the impact of local conserved quantities on the detectable entanglement is discussed., Comment: 21 pages, 8 figures; v2: accepted version; see also the companion paper arXiv:2101.08847 [Phys. Rev. Lett. 126, 190503 (2021)]

Entanglement is not only the resource that fuels many quantum technologies but also plays a key role for some of the most profound open questions of fundamental physics. Experiments controlling quantum systems at the single quantum level may shed light on these puzzles. However, measuring, or even bounding, entanglement experimentally has proven to be an outstanding challenge, especially when the prepared quantum states are mixed. We use entropic uncertainty relations for bipartite systems to derive measurable lower bounds on distillable entanglement. We showcase these bounds by applying them to physical models realizable in cold-atom experiments. The derived entanglement bounds rely on measurements in only two different bases and are generically applicable to any quantum simulation platform., Comment: 6 pages, 2 figures; v2: accepted version; see also the companion paper arXiv:2101.08848 [Phys. Rev. A 103, 052412 (2021)]

The approximation of quantum states with artificial neural networks has gained a lot of attention during the last years. Meanwhile, analog neuromorphic chips, inspired by structural and dynamical properties of the biological brain, show a high energy efficiency in running artificial neural-network architectures for the profit of generative applications. This encourages employing such hardware systems as platforms for simulations of quantum systems. Here we report on the realization of a prototype using the latest spike-based BrainScaleS hardware allowing us to represent few-qubit maximally entangled quantum states with high fidelities. Bell correlations of pure and mixed two-qubit states are well captured by the analog hardware, demonstrating an important building block for simulating quantum systems with spiking neuromorphic chips., Comment: 9+13 pages, 4+2 figures; Submission to SciPost

We study the performance of efficient quantum state tomography methods based on neural network quantum states using measured data from a two-photon experiment. Machine learning inspired variational methods provide a promising route towards scalable state characterization for quantum simulators. While the power of these methods has been demonstrated on synthetic data, applications to real experimental data remain scarce. We benchmark and compare several such approaches by applying them to measured data from an experiment producing two-qubit entangled states. We find that in the presence of experimental imperfections and noise, confining the variational manifold to physical states, i.e. to positive semi-definite density matrices, greatly improves the quality of the reconstructed states but renders the learning procedure more demanding. Including additional, possibly unjustified, constraints, such as assuming pure states, facilitates learning, but also biases the estimator., Comment: v3: published version

We study the out-of-equilibrium dynamics of a dilute, lattice-confined Bose-Fermi mixture initialized in a highly excited state consisting of boson-fermion pairs (doublons) occupying single lattice sites. This system represents a paradigmatic case for studying relaxation dynamics in strongly correlated systems, and provides a versatile platform for studying thermalization and localization phenomena. We provide analytical expressions for the short-time decay of isolated doublons and small doublon clusters due to the competition between tunneling and interparticle interactions. We also discuss a mechanism for long-time decay that crucially depends on the quantum statistics of the particles constituting the doublon, namely, the conversion of pairs of neighboring doublons into an unpaired fermion and a site with a fermion and two bosons. Building on these insights, we develop a cluster expansion method to describe the dynamics in extended systems and compare it to numerically exact matrix product state simulations in one dimension. Finally, we discuss how our predictions can be observed in experiments with ultracold heteronuclear molecules., Comment: 6+6 pages, 3+2 figures, published version

Due to the complexity of the space of quantum many-body states the computation of expectation values by statistical sampling is, in general, a hard task. Neural network representations of such quantum states which can be physically implemented by neuromorphic hardware could enable efficient sampling. A scheme is proposed which leverages this capability to speed up sampling from so-called neural quantum states encoded by a restricted Boltzmann machine. Due to the complex network parameters a direct hardware implementation is not feasible. We overcome this problem by considering a phase reweighting scheme for sampling expectation values of observables. Applying our method to a set of paradigmatic entangled quantum states we find that, in general, the phase-reweighted sampling is subject to a form of sign problem, which renders the sampling computationally costly. The use of neuromorphic chips could allow reducing computation times and thereby extend the range of tractable system sizes., Comment: 15 pages, 7 figures

We study the dynamics of a three-mode bosonic system with mode-changing interactions. For large mode occupations the short-time dynamics is well described by classical mean-field equations allowing us to study chaotic dynamics in the classical system and its signatures in the corresponding quantum dynamics. By introducing a symmetry-breaking term we tune the classical dynamics from integrable to strongly chaotic which we demonstrate by calculating Poincar\'e sections and Lyapunov exponents. The corresponding quantum system features level statistics that change from Poissonian in the integrable to Wigner-Dyson in the chaotic case. We investigate the behavior of out-of-time-ordered correlators (OTOCs), specifically the squared commutator, for initial states located in regular and chaotic regions of the classical mixed phase space and find marked differences between the two cases. The short-time behavior is well captured by semi-classical truncated Wigner simulations directly relating these features to properties of the underlying classical mean field dynamics. We discuss a possible experimental realization of this model system in a Bose-Einstein condensate of rubidium atoms which allows reversing the sign of the Hamiltonian required for measuring OTOCs experimentally., Comment: 13 pages, 6 figures, v2: Sec. IV.C revised and extended. Appendix on truncated Wigner simulations added. v3: added new method for calculation of classical Lyapunov exponents, discussion of OTOCs updated, wording changes

We augment the information extractable from a single absorption image of a spinor Bose-Einstein condensate by coupling to initially empty auxiliary hyperfine states. Performing unitary transformations in both, the original and auxiliary hyperfine manifold, enables the simultaneous measurement of multiple spin-1 observables. We apply this scheme to an elongated atomic cloud of $ ^{87} $Rb to simultaneously read out three orthogonal spin directions and with that directly access the spatial spin structure. The readout even allows the extraction of quantum correlations which we demonstrate by detecting spin nematic squeezing without state tomography., Comment: 11 pages, 4 figures, 5 supplementary figures

The semi-classical discrete truncated Wigner approximation (dTWA) has recently been proposed as a simulation method for spin-$1/2$ systems. While it appears to provide a powerful approach which shows promising results in higher dimensions and for systems with long-range interactions, its performance is still not well understood in general. Here we perform a systematic benchmark on the one-dimensional transverse-field Ising model and point to limitations of the approximation arising after sudden quenches into the quantum critical regime. Our procedure allows to identify the limitations of the semi-classical simulations and with that to determine the regimes and questions where quantum simulators can provide information which is inaccessible to semi-classics., Comment: 9 pages, 10 figures

The near-critical unitary dynamics of quantum Ising spin chains in transversal and longitudinal magnetic fields is studied using an artificial neural network representation of the wave function. A focus is set on strong spatial correlations which build up in the system following a quench into the vicinity of the quantum critical point. We compare correlations observed following reinforcement learning of the network states with analytical solutions in integrable cases and tDMRG simulations, as well as with predictions from a semi-classical discrete Truncated Wigner analysis. While the semi-classical approach excells mainly at short times and for small transverse fields, the neural-network representation provides accurate results for a much wider range of parameters. Where long-range spin-spin correlations build up in the long-time dynamics we find qualitative agreement with exact results while quantitative deviations are of similar size as for the semi-classically predicted correlations, and slow convergence is observed when increasing the number of hidden neurons., Comment: 11 pages, 8 figures

A key resource for distributed quantum-enhanced protocols is entanglement between spatially separated modes. Yet, the robust generation and detection of nonlocal entanglement between spatially separated regions of an ultracold atomic system remains a challenge. Here, we use spin mixing in a tightly confined Bose-Einstein condensate to generate an entangled state of indistinguishable particles in a single spatial mode. We show experimentally that this local entanglement can be spatially distributed by self-similar expansion of the atomic cloud. Spatially resolved spin read-out is used to reveal a particularly strong form of quantum correlations known as Einstein-Podolsky-Rosen steering between distinct parts of the expanded cloud. Based on the strength of Einstein-Podolsky-Rosen steering we construct a witness, which testifies up to genuine five-partite entanglement., Comment: 27 pages, 4 figures, 6 supplementary figures

Out-of-time-order correlations (OTOCs) characterize the scrambling, or delocalization, of quantum information over all the degrees of freedom of a system and thus have been proposed as a proxy for chaos in quantum systems. Recent experimental progress in measuring OTOCs calls for a more thorough understanding of how these quantities characterize complex quantum systems, most importantly in terms of the buildup of entanglement. Although a connection between OTOCs and entanglement entropy has been derived, the latter only quantifies entanglement in pure systems and is hard to access experimentally. In this work, we formally demonstrate that the multiple-quantum coherence spectra, a specific family of OTOCs well known in NMR, can be used as an entanglement witness and as a direct probe of multiparticle entanglement. Our results open a path to experimentally testing the fascinating idea that entanglement is the underlying glue that links thermodynamics, statistical mechanics, and quantum gravity., Comment: 4 pages + references and 8 pages of supplementary materials. v2: as published

Controllable arrays of ions and ultra-cold atoms can simulate complex many-body phenomena and may provide insights into unsolved problems in modern science. To this end, experimentally feasible protocols for quantifying the buildup of quantum correlations and coherence are needed, as performing full state tomography does not scale favorably with the number of particles. Here we develop and experimentally demonstrate such a protocol, which uses time reversal of the many-body dynamics to measure out-of-time-order correlation functions (OTOCs) in a long-range Ising spin quantum simulator with more than 100 ions in a Penning trap. By measuring a family of OTOCs as a function of a tunable parameter we obtain fine-grained information about the state of the system encoded in the multiple quantum coherence spectrum, extract the quantum state purity, and demonstrate the buildup of up to 8-body correlations. Future applications of this protocol could enable studies of many-body localization, quantum phase transitions, and tests of the holographic duality between quantum and gravitational systems., Comment: main text: 7 pages, 4 figures; supplement: 9 pages, 4 figures

We propose a scheme for preparing an ensemble of atoms in a maximally entangled W state by exploiting the Rydberg blockade effect. The success of our protocol is indicated by the detection of an ion, which thus serves as a herald for successful entangled state preparation. The fidelity of the preparation scheme is independent of the number of atoms in the ensemble. Therefore, a small cloud of atoms in a dipole trap randomly loaded from a background gas can be reliably prepared in a maximally entangled state despite of atom number fluctuations., Comment: 7 pages, 3 figures

Understanding and controlling the electronic as well as ro-vibrational motion and, thus, the entire chemical dynamics in molecules is the ultimate goal of ultrafast laser and imaging science. In photochemistry, laser-induced dissociation has become a valuable tool for modification and control of reaction pathways and kinetics. Here, we present a pump-probe study of the dissociation dynamics of H$_2^+$ using ultrashort extreme-ultraviolet (XUV) and near-infrared (IR) laser pulses. The reaction kinematics can be controlled by varying the pump-probe delay. We demonstrate that the nuclear motion through the transition state can be reduced to isolated pairs of initial vibrational states. The dynamics is well reproduced by intuitive semi-classical trajectories on a time-dependent potential curve. From this most fundamental scenario we gain insight in the underlying mechanisms which can be applied as design principles for molecular quantum control, particularly for ultrafast reactions involving protons., Comment: 7 pages and 3 figures. Supplementary materials: 6 pages and 3 figures

In an ensemble of laser-driven atoms involving strongly interacting Rydberg states, the excitation probability is usually strongly suppressed. In contrast, here we identify a regime in which the steady-state Rydberg excited fraction is enhanced by the interaction. This effect is associated with the build-up of many-body coherences, induced by coherent multi-photon excitations between collective states. The excitation enhancement should be observable under currently-existing experimental conditions, and may serve as a direct probe for the presence of coherent multi-photon dynamics involving collective quantum states., Comment: 7 pages, 4 figures

We study the impact of dephasing on the excitation dynamics of a cloud of ultracold two-level Rydberg atoms for both resonant and off-resonant laser excitation, using the wave function Monte Carlo (MCWF) technique. We find that while for resonant laser driving, dephasing mainly leads to an increase of the Rydberg population and a decrease of the Mandel Q parameter, at off-resonant driving strong dephasing toggles between direct excitation of pairs of atoms and subsequent excitation of single atoms, respectively. These two excitation mechanisms can be directly quantified via the pair correlation function, which shows strong suppression of the two-photon resonance peak for strong dephasing. Consequently, qualitatively different dynamics arise in the excitation statistics for weak and strong dephasing in off-resonant excitation. Our findings show that time-resolved excitation number measurements can serve as a powerful tool to identify the dominating process in the system's excitation dynamics., Comment: 10 pages, 10 figures

Light propagation through an ensemble of ultra-cold Rydberg atoms in electromagnetically induced transparency (EIT) configuration is studied. In strongly interacting Rydberg EIT media, non-linear optical effects lead to a non-trivial dependence of the degree of probe beam attenuation on the medium density and on its initial intensity. We develop a Monte Carlo rate equation model that self-consistently includes the effect of the probe beam attenuation to investigate the steady state of the Rydberg medium driven by two laser fields. We compare our results to recent experimental data and to results of other state-of-the-art models for light propagation in Rydberg EIT-media. We find that for low probe field intensities, our results match the experimental data best if a density-dependent dephasing rate is included in the model. At higher probe intensities, our model deviates from other theoretical approaches, as it predicts a spectral asymmetry together with line broadening. These are likely due to off-resonant excitation channels, which however have not been observed in recent experiments. Atomic motion and coupling to additional Rydberg levels are discussed as possible origins for these deviations., Comment: 10 pages, 8 figures

We present a theoretical study of the impact of an electric field combined with a magnetic field on the rotational dynamics of open shell diatomic molecules. Within the rigid rotor approximation, we solve the time-independent Schr\"odinger equation including the fine-structure interactions and the \Lambda-doubling effects. We consider three sets of molecule specific parameters and several field regimes and investigate the interplay between the different interactions identifying the dominant one. The possibility of inducing couplings between the spin and rotational degrees of freedom is demonstrated., Comment: 11 pages, 16 figures

Interfacing light and matter at the quantum level is at the heart of modern atomic and optical physics and enables new quantum technologies involving the manipulation of single photons and atoms. A prototypical atom-light interface is electromagnetically induced transparency, in which quantum interference gives rise to hybrid states of photons and atoms called dark-state polaritons. We have observed individual dark-state polaritons as they propagate through an ultracold atomic gas involving Rydberg states. Strong long-range interactions between Rydberg atoms give rise to an effective interaction blockade for dark-state polaritons, which results in large optical nonlinearities and modified polariton number statistics. The observed statistical fluctuations drop well below the quantum noise limit indicating that photon correlations modified by the strong interactions have a significant back-action on the Rydberg atom statistics., Comment: 7 pages, 4 figures