Publications

Increasing the probability of quantum tunneling between two states, while keeping constant the resources of the underlying physical system, is a task of key importance in several physical contexts and platforms, including ultracold atoms confined by double-well potentials and superconducting qubits. We propose a novel ancillary assisted protocol showing that when a quantum system—such as a qubit—is coupled to an ancilla, one can learn the optimal ancillary component and its coupling, to increase the tunneling probability. As a case study, we consider a quantum system that, due to the presence of an energy detuning between two modes, cannot transfer by tunneling the particles from one mode to the other. However, it does it through a learned coupling with an ancillary system characterized by a detuning not smaller than the one of the primary system. We provide several illustrative examples for the paradigmatic case of a two-mode system and a two-mode ancilla in the presence of interacting particles. This reduces to a qubit coupled to an ancillary qubit in the case of one particle in the system and one in the ancilla. Our proposal provides an effective method to increase the tunneling probability in all those physical situations where no direct improvement of the system parameters, such as tunneling coefficient or energy detuning, is either possible or resource efficient. Finally, we also argue that the proposed strategy is not hampered by weak coupling to noisy environments.

Measurement-induced phases exhibit unconventional dynamics as emergent collective phenomena, yet their behavior in tailored interacting systems – crucial for quantum technologies – remains less understood. We develop a systematic toolbox to analyze monitored dynamics in long-range interacting systems, relevant to platforms like trapped ions and Rydberg atoms. Our method extends spin-wave theory to general dynamical generators at the quantum trajectory level, enabling access to a broader class of states than approaches based on density matrices. This allows efficient simulation of large-scale interacting spins and captures nonlinear dynamical features such as entanglement and trajectory correlations. We showcase the versatility of our framework by exploring entanglement phase transitions in a monitored spin system with power-law interactions in one and two dimensions, where the entanglement scaling changes from logarithm to volume law as the interaction range shortens, and by dwelling on how our method mitigates experimental post-selection challenges in detecting monitored quantum phases.

Foundation models are highly versatile neural-network architectures capable of processing different data types, such as text and images, and generalizing across various tasks like classification and generation. Inspired by this success, we propose Foundation Neural-Network Quantum States (FNQS) as an integrated paradigm for studying quantum many-body systems. FNQS leverage key principles of foundation models to define variational wave functions based on a single, versatile architecture that processes multimodal inputs, including spin configurations and Hamiltonian physical couplings. Unlike specialized architectures tailored for individual Hamiltonians, FNQS can generalize to physical Hamiltonians beyond those encountered during training, offering a unified framework adaptable to various quantum systems and tasks. FNQS enable the efficient estimation of quantities that are traditionally challenging or computationally intensive to calculate using conventional methods, particularly disorder-averaged observables. Furthermore, the fidelity susceptibility can be easily obtained to uncover quantum phase transitions without prior knowledge of order parameters. These pretrained models can be efficiently fine-tuned for specific quantum systems. The architectures trained in this paper are publicly available at https://huggingface.co/nqs-models, along with examples for implementing these neural networks in NetKet.

Measuring fluctuations in matter’s low-energy excitations is the key to unveiling the nature of the non-equilibrium response of materials. A promising outlook in this respect is offered by spectroscopic methods that address matter fluctuations by exploiting the statistical nature of light-matter interactions with weak few-photon probes. Here we report the first implementation of ultrafast phase randomized tomography, combining pump-probe experiments with quantum optical state tomography, to measure the ultrafast non-equilibrium dynamics in complex materials. Our approach utilizes a time-resolved multimode heterodyne detection scheme with phase-randomized coherent ultrashort laser pulses, overcoming the limitations of phase-stable configurations and enabling a robust reconstruction of the statistical distribution of phase-averaged optical observables. This methodology is validated by measuring the coherent phonon response in α-quartz. By tracking the dynamics of the shot-noise limited photon number distribution of few-photon probes with ultrafast resolution, our results set an upper limit to the non-classical features of phononic state in α-quartz and provide a pathway to access non-equilibrium quantum fluctuations in more complex quantum materials.

This summary of the second Terrestrial Very-Long-Baseline Atom Interferometry (TVLBAI) Workshop provides a comprehensive overview of our meeting held in London in April 2024 (Second Terrestrial Very-Long-Baseline Atom Interferometry Workshop, Imperial College, April 2024), building on the initial discussions during the inaugural workshop held at CERN in March 2023 (First Terrestrial Very-Long-Baseline Atom Interferometry Workshop, CERN, March 2023). Like the summary of the first workshop (Abend et al. in AVS Quantum Sci. 6:024701, 2024), this document records a critical milestone for the international atom interferometry community. It documents our concerted efforts to evaluate progress, address emerging challenges, and refine strategic directions for future large-scale atom interferometry projects. Our commitment to collaboration is manifested by the integration of diverse expertise and the coordination of international resources, all aimed at advancing the frontiers of atom interferometry physics and technology, as set out in a Memorandum of Understanding signed by over 50 institutions (Memorandum of Understanding for the Terrestrial Very Long Baseline Atom Interferometer Study).

We investigate a hole-doped Mott insulator in a slab geometry using the dynamical cluster approximation. We show that the enhancement of the correlation strength at the surface results in the remarkable evolution of the layer-projected Fermi surface, which exhibits holelike pockets in the superficial layers, but gradually evolves into a single electronlike surface in the innermost layers. We further analyze the behavior of the Friedel oscillations induced by the surface and identify distinct signatures of the Fermi surface reconstruction as function of hole doping. In addition, we introduce a computationally tractable quantity that diagnoses the same Fermi surface variation by the concurrent breakdown of Luttinger's theorem. Both the latter quantity and the Friedel oscillations serve as reliable indicators of the change in Fermi surface topology, without the need for any periodization in momentum space.

In most noninteracting quantum systems, the scaling theory of localization predicts one-parameter scaling flow in both ergodic and localized regimes. A corresponding scaling theory of many-body ergodicity breaking is still missing. Here, we introduce a scaling theory of ergodicity breaking in interacting systems, in which the divergent relaxation time follows from the Fermi golden rule, and the observable fluctuations in proximity of the ergodicity breaking critical point are described by the recently introduced fading ergodicity scenario. We argue that, in general, the one-parameter scaling is insufficient, and we show that the scaling theory predicts the critical exponent ν=1 at the ergodicity breaking critical point. Our theoretical framework may serve as a building block for two-parameter scaling theories of many-body systems.

We study the behavior of magic as a bipartite correlation in the quantum Ising chain across its quantum phase transition and at finite temperature. To quantify the magic of partitions rigorously, we formulate a hybrid scheme that combines stochastic sampling of reduced density matrices via quantum Monte Carlo, with state-of-the-art estimators for the robustness of magic - a bona fide measure of magic for mixed states. This allows us to compute the mutual robustness of magic for partitions up to eight sites, embedded into a much larger system. We show how mutual robustness is directly related to critical behaviors: at the critical point, it displays a power-law decay as a function of the distance between partitions, whose exponent is related to the partition size. Once finite temperature is included, mutual magic retains its low temperature value up to an effective critical temperature, whose dependence on size is also algebraic. This suggests that magic, differently from entanglement, does not necessarily undergo a sudden death.

In quantum mechanics, the probability distribution function and full counting statistics play a fundamental role in characterizing the fluctuations of quantum observables, as they encode the complete information about these fluctuations. In this Letter, we measure these two quantities in a trapped-ion quantum simulator for the transverse and longitudinal magnetization within a subsystem. We utilize the toolbox of classical shadows to postprocess the measurements performed in random bases. The measurement scheme efficiently allows access to the full counting statistics and probability distribution function of all possible operators on desired choices of subsystems of an extended quantum system.

In this paper we discuss and revisit the finite temperature extension of the renormalization group (RG) treatment of T = 0 field theories, focusing as a case study on the ϕ⁴ model.We first discuss the extension of RG equations of the very same model from T = 0 to finite T in the usual way by resorting to sums on the Matsubara frequencies and fixing the physical temperature parameter T. We show that this approach, although useful for a variety of applications, may lead to the disappearance of the critical points as extracted from the RG flow. Since the identification of fixed points is key in the study of classical and quantum phase transitions, we propose a modification of the usual finite-temperature RG approach by relating the temperature parameter to the running RG scale, T ≡ kT = τk, where kT is the running cutoff for thermal and k is for the quantum fluctuations. Once this dimensionless temperature τ is introduced, we investigate the consequences on the thermal RG approach for the ϕ⁴ model and construct its phase diagram. Finally, we formulate requirements for the phase diagram of the ϕ⁴ theory based on known properties of the quantum and classical phase diagrams of the Ising model.

We study the static and dynamical properties of a harmonically confined Rouse polymer coupled to a fluctuating correlated medium, which affect each other reciprocally during their stochastic evolution. The medium is modeled by a scalar Gaussian field which can feature modes with slow relaxation and long-range spatial correlations. We show that these modes affect the long-time behavior of the average position of the center of mass of the polymer, which, after a displacement, turns out to relax algebraically towards its equilibrium value. This is a manifestation of the non-Markovian nature of the effective evolution of the position of the center of mass, once the degrees of freedom of the medium have been integrated out. In contrast, we show that the coupling to the medium speeds up the relaxation of higher Rouse modes. We further characterize the typical size of the polymer as a function of its polymerization degree and of the correlation length of the medium, particularly when the system is driven out of equilibrium via the application of a constant external driving force. Finally, we study the response of a linear polymer to a tensile force acting on its terminal monomers.

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.

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.

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 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.