Publications

We introduce a methodology to estimate nonstabilizerness or “magic,” a key resource for quantum complexity, with neural quantum states (NQS). Our framework relies on two schemes based on Monte Carlo sampling to quantify nonstabilizerness via stabilizer Rényi entropy (SRE) in arbitrary variational wave functions. When combined with NQS, this approach is effective for systems with strong correlations and in dimensions larger than one, unlike tensor network methods. First, we study the magic content in an ensemble of random NQS, demonstrating that neural network parametrizations of the wave function capture finite nonstabilizerness besides large entanglement. Second, we investigate the nonstabilizerness in the ground state of the J1-J2 Heisenberg model. In one dimension, we find that the SRE vanishes at the Majumdar-Ghosh point J2 = J1/2, consistent with a stabilizer ground state. In two dimensions, a dip in the SRE is observed near maximum frustration around J2/J1 ≈ 0.6, suggesting a valence bond solid between the two antiferromagnetic phases.

We propose quantum simulation experiments of the Kondo impurity problem using cold alkaline-earth(-like) atoms (AEAs) in a combination of optical lattice and optical tweezer potentials. Within an ab initio model for atomic interactions in the optical potentials, we analyze hallmark signatures of the Kondo effect in a variety of observables accessible in cold-atom quantum simulators. We identify additional terms not part of the textbook Kondo problem, mostly ignored in previous works and giving rise to a direct competition between spin and charge correlations, strongly suppressing Kondo physics. We show that the Kondo effect can be restored by locally adjusting the chemical potential on the impurity site, and we identify realistic parameter regimes and preparation protocols suited to current experiments with AEA arrays. Our work paves the way for quantum simulations of the Kondo problem and offers insights into Kondo physics in unconventional regimes.

Exponentially small spectral gaps are known to be the crucial bottleneck for traditional Quantum Annealing (QA) based on interpolating between two Hamiltonians, a simple driving term and the complex problem to be solved, with a linear schedule in time. One of the simplest models exhibiting exponentially small spectral gaps is a ferromagnetic Ising ring with a single antiferromagnetic bond introducing frustration. Previous studies of this model have explored continuous-time diabatic QA, where optimized non-adiabatic annealing schedules provided good solutions, avoiding exponentially large annealing times. In our work, we move to a digital framework of Variational Quantum Algorithms, and present two main results: (1) we show that the model is digitally controllable with a scaling of resources that grows quadratically with the system size, achieving the exact solution using the Quantum Approximate Optimization Algorithm; (2) We combine a technique of quantum control—the Chopped RAndom Basis method—and digitized quantum annealing to construct smooth digital schedules yielding optimal solutions with very high accuracy.

The spectral properties of the Heisenberg spin-1/2 chain with random fields are analyzed in light of recent works on the renormalization-group flow of the Anderson model in infinite dimension. We reconstruct the β function of the order parameter from the numerical data, and observe that it may not admit a one-parameter scaling form and a simple Wilson-Fisher fixed point. Rather, it appears to be more compatible with a two parameter, Berezinskii–Kosterlitz-Thouless-like flow with a line of fixed points (the many-body localized phase) terminating at the localization transition critical point. We argue that this renormalization group framework provides a more coherent and intuitive explanation of numerical data, up to the system sizes available with the present technology.

We study monitored quantum dynamics of infinite-range interacting bosonic systems in the thermodynamic limit. We show that under semiclassical assumptions, the quantum fluctuations along single monitored trajectories adopt a deterministic limit for both quantum-jump and state-diffusion unravelings, and they can be exactly solved. In particular, the hierarchical structure of the equations of motion explains the coincidence of entanglement criticalities and dissipative phase transitions found in previous finite-size numerical studies. We illustrate the findings on a Bose-Hubbard dimer and a collective spin system.

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. This study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water-methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

We analyze two simple model planar molecules: an ionic molecule with D3 symmetry and a covalent molecule with D6 symmetry. Both symmetries allow the existence of chiral molecular orbitals and normal modes that are coupled to each other in a Jahn-Teller manner, invariant under U(1) symmetry with the generator a pseudo-angular momentum. In the ionic molecule, the chiral mode possesses an electric dipole but lacks physical angular momentum, whereas in the covalent molecule, the situation is reversed. In spite of that, we show that in both cases, the chiral modes can be excited by a circularly polarized light and are subsequently able to induce rotational motion of the entire molecule.

Recent work has shown that the entanglement of finite-temperature eigenstates in chaotic quantum many-body local Hamiltonians can be accurately described by an ensemble of random states with an internal U(1) symmetry. We build upon this result to investigate the universal symmetry-breaking properties of such eigenstates. As a probe of symmetry breaking, we employ the entanglement asymmetry, a quantum information observable that quantifies the extent to which symmetry is broken in a subsystem. This measure enables us to explore the finer structure of finite-temperature eigenstates in terms of the U(1)-symmetric random state ensemble; in particular, the relation between the Hamiltonian and the effective conserved charge in the ensemble. Our analysis is supported by analytical calculations for the symmetric random states, as well as exact numerical results for the Mixed-Field Ising spin-1/2 chain, a paradigmatic model of quantum chaoticity.

In quantum resource theory, one is often interested in identifying which states serve as the best resources for particular quantum tasks. If a relative comparison between quantum states can be made, this gives rise to a partial order, where states are ordered according to their suitability to act as a resource. In the literature, various different partial orders for a variety of quantum resources have been proposed. In discrete variable systems, vector majorization ofWigner functions in discrete phase space provides a natural partial order between quantum states. In the continuous variable case, a natural counterpart would be continuous majorization of Wigner functions in quantum phase space. Indeed, this concept was recently proposed and explored (mostly restricting to the single-mode case) by Van Herstraeten et al. [Quantum 7, 1021 (2023)]. In this work, we develop the theory of continuous majorization in the general N-mode case. In addition, we propose extensions to include states with finite Wigner negativity. For the special case of the convex hull of N-mode Gaussian states, we prove a conjecture made by Van Herstraeten, Jabbour, and Cerf.We also prove a phase space counterpart of Uhlmann's theorem of majorization.

We discuss the dependence of the critical properties of the Anderson model on the dimension d in the language of β-function and renormalization group recently introduced in Vanoni et al. [C. Vanoni et al., Proc. Natl. Acad. Sci. U.S.A. 121, e2401955121 (2024)] in the context of Anderson transition on random regular graphs. We show how in the delocalized region, including the transition point, the one-parameter scaling part of the β-function for the fractal dimension D1 evolves smoothly from its d = 2 form, in which β2 ≤ 0, to its β∞ ≥ 0 form, which is represented by the random regular graph (RRG) result. We show how the ε = d − 2 expansion and the 1/d expansion around the RRG result can be reconciled and how the initial part of a renormalization group trajectory governed by the irrelevant exponent y depends on dimensionality. We also show how the irrelevant exponent emerges out of the high-gradient terms of expansion in the nonlinear sigma model and put forward a conjecture about a lower bound for the fractal dimension. The framework introduced here may serve as a basis for investigations of disordered many-body systems and of more general nonequilibrium quantum systems.

In this paper, we investigate the quench dynamics of the negativity and fermionic negativity Hamiltonians in free fermionic systems. We do this by generalizing a recently developed quasiparticle picture for the entanglement Hamiltonians to tripartite geometries. We obtain analytic expressions for these quantities, which are then extensively checked against previous results and numerics. In particular, we find that the standard negativity Hamiltonian contains both non-local hopping terms and four-fermion interactions, whereas the fermionic version is purely quadratic. However, despite their marked difference, we show that the logarithmic negativity obtained from either is identical in the ballistic scaling limit, as are their symmetry resolution.

Quantum many-body scarred systems exhibit atypical dynamical behavior, evading thermalization and featuring periodic state revivals. In this Letter, we investigate the impact of projective measurements on the dynamics in the scar subspace for the paradigmatic PXP model, revealing that they can either disrupt or enhance the revivals. Local measurements performed at random times rapidly erase the system's memory of its initial conditions, leading to fast steady state relaxation. In contrast, a periodic monitoring amplifies recurrences and preserves the coherent dynamics over extended timescales. We identify a measurement-induced phase resynchronization, countering the natural dephasing of quantum scars, as the key mechanism underlying this phenomenon.

Accurate and reliable algorithms to solve complex impurity problems are instrumental to a routine use of quantum embedding methods for material discovery. In this context, we employ an efficient selected configuration-interaction impurity solver to investigate the role of bath discretization—specifically, bath size and parametrization—in Hamiltonian-based cluster dynamical mean-field theory (CDMFT) for the one- and two-orbital Hubbard models. We consider two- and four-site clusters for the single-orbital model and a two-site cluster for the two-orbital model. Our results demonstrate that, for small bath sizes, the choice of parametrization can significantly influence the solution, highlighting the importance of systematic convergence checks. Comparing different bath parametrizations not only reveals the robustness of a given solution but can also provide insights into the nature of different solutions and potential instabilities of the paramagnetic state. We present an extensive analysis of the zero-temperature Mott transition of the paramagnetic half-filled single-band Hubbard model, benchmarking our findings against previous literature. We find that, for the single-band model the dependence on parametrization is weak for the largest bath sizes accessible with ASCI, while a tendency towards a nematic solution can be seen when the bath size is small. Building on this, we extend our study to the multiband regime, where we present a systematic analysis at zero temperature for two orbitals and a two-site cluster. This setup allows us to assess the effect of nearest-neighbor dynamical correlations on the multi-orbital Mott transition. In this case, some quantitative dependence on the parametrization is retained for the two-orbital model, for instance in the value of the critical interaction for a Mott transition.

We address weak-coupling frictional sliding with phononic dissipation by means of analytic many-body techniques. Our model consists of a particle (the “slider”) moving through a two-or three-dimensional crystal and interacting weakly with its atoms, and therefore, exciting phonons. By means of linear response theory, we obtain explicit expressions for the friction force slowing down the slider as a function of its speed, and compare them to the friction obtained by simulations, demonstrating a remarkable accord.

We study the quantum transport generated by the bipartite entanglement in two-dimensional conformal field theory at finite density with the U(1) × U(1) symmetry associated to the conservation of the electric charge and of the helicity. The bipartition given by an interval is considered, either on the line or on the circle. The continuity equations and the corresponding conserved quantities for the modular flows of the currents and of the energy-momentum tensor are derived. We investigate the mean values of the associated currents and their quantum fluctuations in the finite density representation, which describe the properties of the modular quantum transport. The modular analogues of the Johnson- Nyquist law and of the fluctuation-dissipation relation are found, which encode the thermal nature of the modular evolution.

The Mpemba effect, in which a hotter system can equilibrate faster than a cooler one, has long been a subject of fascination in classical physics. In the past few years, notable theoretical and experimental progress has been made in understanding its occurrence in both classical and quantum systems. In this Perspective, we provide a concise overview of recent work and open questions on the Mpemba effect in quantum systems, with a focus on both open and isolated dynamics, which give rise to distinct manifestations of this anomalous non-equilibrium phenomenon. We discuss key theoretical frameworks, highlight experimental observations and explore the fundamental mechanisms that give rise to anomalous relaxation behaviours. Particular attention is given to the role of quantum fluctuations, integrability and symmetry in shaping equilibration pathways.

The essence of the Mpemba effect is that nonequilibrium systems may relax faster the further they are from their equilibrium configuration. In the quantum realm, this phenomenon arises in closed systems dynamics and is witnessed by features of symmetry and entanglement. Here, we study the quantum Mpemba effect in charge-preserving random circuits, combining extensive numerical simulations and analytical arguments. We show that the more asymmetric certain classes of initial states (tilted ferromagnets) are, the faster they restore symmetry and reach the grand-canonical ensemble. Conversely, other classes of states (tilted antiferromagnets) do not show the Mpemba effect. We provide a simple and general mechanism underlying the effect, based on the spreading of nonconserved operators in terms of conserved densities. Grounded only in locality, unitarity, and symmetry, our analysis clarifies the emergence of Mpemba physics in chaotic quantum systems.

We propose powering a quantum clock with the nonthermal resources offered by the stationary state of an integrable quantum spin chain, driven out of equilibrium by a sudden quench. We describe the dynamics of the clock as a “biased random walk.” The bias level is directly related to the negativity of the steady-state response function. Using experimentally relevant examples of quantum spin chains, we suggest that crossing a phase transition point is crucial for the clock’s operation. The coupling takes place through a global observable and, in this case, the battery lifespan is found to be extensive in its size.

Quantum secure direct communication (QSDC) is an evolving quantum communication framework based on transmitting secure information directly through a quantum channel, without relying on key-based encryption such as in quantum key distribution (QKD). Optical QSDC protocols, utilizing discrete and continuous variable encodings, show great promise for future technological applications. We present the first table-top continuous-variable QSDC proof of principle, analyzing its implementation and comparing the use of coherent against squeezed light sources. A simple beam-splitter attack is analyzed by using Wyner wiretap channel theory. Our study illustrates the advantage of squeezed states over coherent ones for enhanced security and reliable communication in lossy and noisy channels. Our practical implementation, utilizing mature telecom components, could foster secure quantum metropolitan networks compatible with advanced multiplexing systems.

In the collective photon emission from atomic clouds both the atomic transition frequency and the decay rate are modified compared to a single isolated atom, leading to the effects of superradiance and subradiance. In this article, we analyze the properties of the Euclidean random matrix associated with the radiative dynamics of a cold atomic cloud, previously investigated in the contexts of photon localization and Dicke super- and subradiance.We present evidence of a phase transition, surprisingly controlled by the cooperativeness parameter, rather than the spatial density or the diagonal disorder. The numerical results corroborate the occurrence of such a phase transition at a critical value of the cooperativeness parameter, above which the lower edge of the spectrum vanishes, exhibiting a macroscopic accumulation of eigenvalues. Independent evaluations based on the two phenomena provide the same value for the critical cooperativeness parameter.