Publications

We develop a systematic approach to compute the subsystem trace distances and relative entropies for subsystem reduced density matrices associated to excited states in different symmetry sectors of a 1+1 dimensional conformal field theory having an internal U(1) symmetry. We provide analytic expressions for the charged moments corresponding to the resolution of both relative entropies and distances for general integer n. For the relative entropies, these formulas are manageable and the analytic continuation to n = 1 can be worked out in most of the cases. Conversely, for the distances the corresponding charged moments become soon untreatable as n increases. A remarkable result is that relative entropies and distances are the same for all symmetry sectors, i.e. they satisfy entanglement equipartition, like the entropies. Moreover, we exploit the OPE expansion of composite twist fields, to provide very general results when the subsystem is a single interval much smaller than the total system. We focus on the massless compact boson and our results are tested against exact numerical calculations in the XX spin chain.

We study the critical properties of the 3d O(2) universality class in bounded domains through Monte Carlo simulations of the clock model. We use an improved version of the latter, chosen to minimize finite-size corrections at criticality, with 8 orientations of the spins and in the presence of vacancies. The domain chosen for the simulations is the slab configuration with fixed spins at the boundaries. We obtain the universal critical magnetization profile and two-point correlations, which favorably compare with the predictions of the critical geometry approach based on the Yamabe equation. The main result is that the correlations, once the dimensionful contributions are factored out with the critical magnetization profile, are shown to depend only on the distance between the points computed using a metric found solving the corresponding fractional Yamabe equation. The quantitative comparison with the corresponding results for the Ising model at criticality is shown and discussed. Moreover, from the magnetization profiles the critical exponent η is extracted and found to be in reasonable agreement with up-to-date results.

We consider the problem of the decomposition of the Rényi entanglement entropies in theories with a non-abelian symmetry by doing a thorough analysis of Wess-Zumino-Witten (WZW) models. We first consider SU(2)k as a case study and then generalise to an arbitrary non-abelian Lie group. We find that at leading order in the subsystem size L the entanglement is equally distributed among the different sectors labelled by the irreducible representation of the associated algebra. We also identify the leading term that breaks this equipartition: it does not depend on L but only on the dimension of the representation. Moreover, a log log L contribution to the Rényi entropies exhibits a universal prefactor equal to half the dimension of the Lie group.

In the conformal field theories given by the Ising and Dirac models, when the system is in the ground state, the moments of the reduced density matrix of two disjoint intervals and of its partial transpose have been written as partition functions on higher genus Riemann surfaces with symmetry. We show that these partition functions can be expressed as the grand canonical partition functions of the two-dimensional two component classical Coulomb gas on certain circular lattices at specific values of the coupling constant.

Ultracold fermionic atoms are proposed to be used in 1D optical lattices to quantum simulate the electronic transport in quantum cascade laser (QCL) structures. The competition between the coherent tunneling among (and within) the wells and the dissipative decay at the basis of lasing is discussed. In order to validate the proposed simulation scheme, such competition is quantitatively addressed in a simplified 1D model. The existence of optimal relationships between the model parameters is shown, maximizing the particle current, the population inversion (or their product), and the stimulated emission rate. This substantiates the concept of emulating the QCL operation mechanisms in cold-atom optical lattice simulators, laying the groundwork for addressing open questions, such as the impact of electron–electron scattering and the origin of transport-induced noise, in the design of new-generation QCLs.

We study the impact of quenched disorder on the dynamics of locally constrained quantum spin chains, that describe 1D arrays of Rydberg atoms in both the frozen (Ising-type) and dressed (XY-type) regime. Performing large-scale numerical experiments, we observe no trace of many-body localization even at large disorder. Analyzing the role of quenched disorder terms in constrained systems we show that they act in two, distinct and competing ways: as an on-site disorder term for the basic excitations of the system, and as an interaction between excitations. The two contributions are of the same order, and as they compete (one towards localization, the other against it), one does never enter a truly strong disorder, weak interaction limit, where many-body localization occurs. Such a mechanism is further clarified in the case of XY-type constrained models: there, a term which would represent a bona fide local quenched disorder term acting on the excitations of the clean model must be written as a series of nonlocal terms in the unconstrained variables. Our observations provide a simple picture to interpret the role of quenched disorder that could be immediately extended to other constrained models or quenched gauge theories.

We revisit early suggestions to observe spin-charge separation (SCS) in cold-atom settings in the time domain by studying one-dimensional repulsive Fermi gases in a harmonic potential, where pulse perturbations are initially created at the center of the trap. We analyze the subsequent evolution using generalized hydrodynamics (GHD), which provide an exact description, at large space-time scales, for arbitrary temperature T, particle density, and interactions. At T=0 and vanishing magnetic field, we find that, after a nontrivial transient regime, spin and charge dynamically decouple up to perturbatively small corrections which we quantify. In this limit, our results can be understood based on a simple phase-space hydrodynamic picture. At finite temperature, we solve numerically the GHD equations, showing that for low T>0 effects of SCS survive and characterize explicitly the value of T for which the two distinguishable excitations melt.

We study the reduced dynamics of open quantum spin chains of arbitrary length N with nearest-neighbor XX interactions, immersed within an external constant magnetic field along the z direction, the end spins of which are weakly coupled to heat baths at different temperatures, via energy-preserving couplings. We find the analytic expression of the unique stationary state of the master equation obtained in the so-called global approach based on the spectralization of the full-chain Hamiltonian. Hinging upon the explicit stationary state, we reveal the presence of sink and source terms in the spin-flow continuity equation and compare their behavior with that of the stationary heat flow. Moreover, we also obtain analytic expressions for the steady-state two-spin reduced density matrices and for their concurrence. We then set up an algorithm suited to compute the stationary bipartite entanglement along the chain and to study its dependence on the Hamiltonian parameters and on the bath temperatures.

The study of the entanglement dynamics plays a fundamental role in understanding the behaviour of many-body quantum systems out of equilibrium. In the presence of a globally conserved charge, further insights are provided by the knowledge of the resolution of entanglement in the various symmetry sectors. Here, we carry on the program we initiated in Parez et al (2021 Phys. Rev. B 103 L041104), for the study of the time evolution of the symmetry resolved entanglement in free fermion systems. We complete and extend our derivations also by defining and quantifying a symmetry resolved mutual information. The entanglement entropies display a time delay that depends on the charge sector that we characterise exactly. Both entanglement entropies and mutual information show effective equipartition in the scaling limit of large time and subsystem size. Furthermore, we argue that the behaviour of the charged entropies can be quantitatively understood in the framework of the quasiparticle picture for the spreading of entanglement, and hence we expect that a proper adaptation of our results should apply to a large class of integrable systems. We also find that the number entropy grows logarithmically with time before saturating to a value proportional to the logarithm of the subsystem size.

Out-of-time-order correlators (OTOCs) have become established as a tool to characterise quantum information dynamics and thermalization in interacting quantum many-body systems. It was recently argued that the expected exponential growth of the OTOC is connected to the existence of correlations beyond those encoded in the standard Eigenstate Thermalization Hypothesis (ETH). We show explicitly, by an extensive numerical analysis of the statistics of operator matrix elements in conjunction with a detailed study of OTOC dynamics, that the OTOC is indeed a precise tool to explore the fine details of the ETH. In particular, while short-time dynamics is dominated by correlations, the long-time saturation behavior gives clear indications of an operator-dependent energy scale ωGOE associated to the emergence of an effective Gaussian random matrix theory. We provide an estimation of the finite-size scaling of ωGOE for the general class of observables composed of sums of local operators in the infinite-temperature regime and found linear behavior for the models considered.

We study the quantum quench in two coupled Tomonaga-Luttinger Liquids (TLLs), from the off-critical to the critical regime, relying on the conformal field theory approach and the known solutions for single TLLs. We consider a squeezed form of the initial state, whose low energy limit is fixed in a way to describe a massive and a massless mode, and we encode the non-equilibrium dynamics in a proper rescaling of the time. In this way, we compute several correlation functions, which at leading order factorize into multipoint functions evaluated at different times for the two modes. Depending on the observable, the contribution from the massive or from the massless mode can be the dominant one, giving rise to exponential or power-law decay in time, respectively. Our results find a direct application in all the quench problems where, in the scaling limit, there are two independent massless fields: these include the Hubbard model, the Gaudin-Yang gas, and tunnel-coupled tubes in cold atoms experiments.

Atomtronics deals with matter-wave circuits of ultracold atoms manipulated through magnetic or laser-generated guides with different shapes and intensities. In this way, new types of quantum networks can be constructed in which coherent fluids are controlled with the know-how developed in the atomic and molecular physics community. In particular, quantum devices with enhanced precision, control, and flexibility of their operating conditions can be accessed. Concomitantly, new quantum simulators and emulators harnessing on the coherent current flows can also be developed. Here, the authors survey the landscape of atomtronics-enabled quantum technology and draw a roadmap for the field in the near future. The authors review some of the latest progress achieved in matter-wave circuits' design and atom-chips. Atomtronic networks are deployed as promising platforms for probing many-body physics with a new angle and a new twist. The latter can be done at the level of both equilibrium and nonequilibrium situations. Numerous relevant problems in mesoscopic physics, such as persistent currents and quantum transport in circuits of fermionic or bosonic atoms, are studied through a new lens. The authors summarize some of the atomtronics quantum devices and sensors. Finally, the authors discuss alkali-earth and Rydberg atoms as potential platforms for the realization of atomtronic circuits with special features.

Quantum many-body systems are characterized by patterns of correlations defining highly nontrivial manifolds when interpreted as data structures. Physical properties of phases and phase transitions are typically retrieved via correlation functions, that are related to observable response functions. Recent experiments have demonstrated capabilities to fully characterize quantum many-body systems via wave-function snapshots, opening new possibilities to analyze quantum phenomena. Here, we introduce a method to data mine the correlation structure of quantum partition functions via their path integral (or equivalently, stochastic series expansion) manifold. We characterize path-integral manifolds generated via state-of-the-art quantum Monte Carlo methods utilizing the intrinsic dimension (ID) and the variance of distances between nearest-neighbor (NN) configurations: the former is related to data-set complexity, while the latter is able to diagnose connectivity features of points in configuration space. We show how these properties feature universal patterns in the vicinity of quantum criticality, that reveal how data structures simplify systematically at quantum phase transitions. This is further reflected by the fact that both ID and variance of NN distances exhibit universal scaling behavior in the vicinity of second-order and Berezinskii-Kosterlitz-Thouless critical points. Finally, we show how non-Abelian symmetries dramatically influence quantum data sets, due to the nature of (noncommuting) conserved charges in the quantum case. Complementary to neural-network representations, our approach represents a first elementary step towards a systematic characterization of path-integral manifolds before any dimensional reduction is taken, that is informative about universal behavior and complexity, and can find immediate application to both experiments and Monte Carlo simulations.

We study the heat statistics of a multilevel N-dimensional quantum system monitored by a sequence of projective measurements. The late-time, asymptotic properties of the heat characteristic function are analyzed in the thermodynamic limit of a high, ideally infinite, number M of measurements (M→∞). In this context, the conditions allowing for an infinite-temperature thermalization (ITT), induced by the repeated monitoring of the quantum system, are discussed. We show that ITT is identified by the fixed point of a symmetric random matrix that models the stochastic process originated by the sequence of measurements. Such fixed point is independent on the nonequilibrium evolution of the system and its initial state. Exceptions to ITT, which we refer to as partial thermalization, take place when the observable of the intermediate measurements is commuting (or quasicommuting) with the Hamiltonian of the quantum system or when the time interval between measurements is smaller or comparable with the system energy scale (quantum Zeno regime). Results on the limit of infinite-dimensional Hilbert spaces (N→∞), describing continuous systems with a discrete spectrum, are also presented. We show that the order of the limits M→∞ and N→∞ matters: When N is fixed and M diverges, then ITT occurs. In the opposite case, the system becomes classical, so that the measurements are no longer effective in changing the state of the system. A nontrivial result is obtained fixing M/N2 where instead partial ITT occurs. Finally, an example of partial thermalization applicable to rotating two-dimensional gases is presented.

A number of experimental platforms relevant for quantum simulations, ranging from arrays of superconducting circuits hosting correlated states of light to ultracold atoms in optical lattices in the presence of controlled dissipative processes. Their theoretical understanding is hampered by the exponential scaling of their Hilbert space and by their intrinsic nonequilibrium nature, limiting the applicability of many traditional approaches. In this work, we extend the nonequilibrium bosonic dynamical mean-field theory (DMFT) to Markovian open quantum systems. Within DMFT, a Lindblad master equation describing a lattice of dissipative bosonic particles is mapped onto an impurity problem describing a single site embedded in its Markovian environment and coupled to a self-consistent field and to a non-Markovian bath, where the latter accounts for fluctuations beyond Gutzwiller mean-field theory due to the finite lattice connectivity. We develop a nonperturbative approach to solve this bosonic impurity problem, which dresses the impurity, featuring Markovian dissipative channels, with the non-Markovian bath, in a self-consistent scheme based on a resummation of noncrossing diagrams. As a first application of our approach, we address the steady state of a driven-dissipative Bose-Hubbard model with two-body losses and incoherent pump. We show that DMFT captures hopping-induced dissipative processes, completely missed in Gutzwiller mean-field theory, which crucially determine the properties of the normal phase, including the redistribution of steady-state populations, the suppression of local gain, and the emergence of a stationary quantum-Zeno regime. We argue that these processes compete with coherent hopping to determine the phase transition toward a nonequilibrium superfluid, leading to a strong renormalization of the phase boundary at finite connectivity. We show that this transition occurs as a finite-frequency instability, leading to an oscillating-in-time order parameter, that we connect with a quantum many-body synchronization transition of an array of quantum van der Pol oscillators.

Inflation solves several cosmological problems at the classical and quantum level, with a strong agreement between the theoretical predictions of well-motivated inflationary models and observations. In this Letter, we study the corrections induced by dynamical collapse models, which phenomenologically solve the quantum measurement problem, to the power spectrum of the comoving curvature perturbation during inflation and the radiation-dominated era. We find that the corrections are strongly negligible for the reference values of the collapse parameters.

Dirac fermions play a central role in the study of topological phases, for they can generate a variety of exotic states, such as Weyl semimetals and topological insulators. The control and manipulation of Dirac fermions constitute a fundamental step toward the realization of novel concepts of electronic devices and quantum computation. By means of Angle-Resolved Photo-Emission Spectroscopy (ARPES) experiments and ab initio simulations, here, we show that Dirac states can be effectively tuned by doping a transition metal sulfide, BaNiS2, through Co/Ni substitution. The symmetry and chemical characteristics of this material, combined with the modification of the charge-transfer gap of BaCo1−xNixS2 across its phase diagram, lead to the formation of Dirac lines, whose position in k-space can be displaced along the Γ − M symmetry direction and their form reshaped. Not only does the doping x tailor the location and shape of the Dirac bands, but it also controls the metal-insulator transition in the same compound, making BaCo1−xNixS2 a model system to functionalize Dirac materials by varying the strength of electron correlations.

We report V51 NMR and inelastic neutron scattering (INS) measurements on a quasi-1D antiferromagnet BaCo2V2O8 under transverse field along the [010] direction. The scaling behavior of the spin-lattice relaxation rate above the Néel temperatures unveils a 1D quantum critical point (QCP) at Hc1D≈4.7 T, which is masked by the 3D magnetic order. With the aid of accurate analytical analysis and numerical calculations, we show that the zone center INS spectrum at Hc1D is precisely described by the pattern of the 1D quantum Ising model in a magnetic field, a class of universality described in terms of the exceptional E8 Lie algebra. These excitations are nondiffusive over a certain field range when the system is away from the 1D QCP. Our results provide an unambiguous experimental realization of the massive E8 phase in the compound, and open a new experimental route for exploring the dynamics of quantum integrable systems as well as physics beyond integrability.

Shooting glass beads across the inside of a satellite could probe the limits of quantum wave behaviour. Here’s how. [Figure not available: see fulltext.]

We study the connection between the magnetization profiles of models described by a scalar field with marginal interaction term in a bounded domain and the solutions of the so-called Yamabe problem in the same domain, which amounts to finding a metric having constant curvature. Taking the slab as a reference domain, we first study the magnetization profiles at the upper critical dimensions d=3, 4, 6 for different (scale-invariant) boundary conditions. By studying the saddle-point equations for the magnetization, we find general formulas in terms of Weierstrass elliptic functions, extending exact results known in literature and finding ones for the case of percolation. The zeros and poles of the Weierstrass elliptic solutions can be put in direct connection with the boundary conditions. We then show that, for any dimension d, the magnetization profiles are solution of the corresponding integer Yamabe equation at the same d and with the same boundary conditions. The magnetization profiles in the specific case of the four-dimensional Ising model with fixed boundary conditions are compared with Monte Carlo simulations, finding good agreement. These results explicitly confirm at the upper critical dimension recent results presented in Gori and Trombettoni [J. Stat. Mech: Theory Exp. (2020) 0632101742-546810.1088/1742-5468/ab7f32].