Publications from 2024

The Lattice Gauge Theory Hilbert space is divided into gauge-invariant sectors selected by the background charges. Such a projector can be directly embedded in a tensor network ansatz for gauge-invariant states as originally discussed in [Phys. Rev. B 83, 115127 (2011)] and in [Phys. Rev. X 4, 041024 (2014)] in the context of PEPS. The original ansatz is based on sparse tensors, though parts of them are not explicitly symmetric, and thus their actual implementation in numerical simulations has been hindered by the complexity of developing ad hoc libraries. Here we provide a new PEPS tensor network formulation of gauge-invariant theories purely based on symmetric elementary tensors. The new formulation can be implemented in numerical simulation using available state-of-the-art tensor network libraries but also holds interest from a purely theoretical perspective since it requires embedding the original gauge theory with gauge symmetry G into an enlarged globally symmetric theory with symmetry GxG. By revisiting the original ansatz in the modern landscape of i) duality transformations between gauge and spin systems, ii) finite depth quantum circuits followed by measurements that allow generating topologically ordered states, and iii) Clifford enhanced tensor networks, we show that such a new formulation provides a novel duality transformation between lattice gauge theories and specific sectors of globally invariant systems.

Floquet topological phases emerge when systems are periodically driven out-of-equilibrium. They gained attention due to their external control, which allows to simulate a wide variety of static systems by just tuning the external field in the high frequency regime. However, it was soon clear that their relevance goes beyond that, as for lower frequencies, anomalous phases without a static counterpart are present and the bulk-to-boundary correspondence can fail. In this work we discuss the important role of resonances in Floquet phases. For that, we present a method to find analytical solutions when the frequency of the drive matches the band gap, extending the well-known high frequency analysis of Floquet systems. With this formalism, we show that the topology of Floquet phases with resonances can be accurately captured in analytical terms. We also find a bulk-to-boundary correspondence between the number of edge states in finite systems and a set of topological invariants in different frames of reference, which crucially do not explicitly involve the micromotion. To illustrate our results, we periodically drive a SSH chain and a $\pi$-flux lattice, showing that our findings remain valid in various two-band systems and in different dimensions. In addition, we notice that the competition between rotating and counter-rotating terms must be carefully treated when the undriven system is a semi-metal. To conclude, we discuss the implications to experimental setups, including the direct detection of anomalous topological phases and the measurement of their invariants.

This work explores the representation of univariate and multivariate functions as matrix product states (MPS), also known as quantized tensor-trains (QTT). It proposes an algorithm that employs iterative Chebyshev expansions and Clenshaw evaluations to represent analytic and highly differentiable functions as MPS Chebyshev interpolants. It demonstrates rapid convergence for highly-differentiable functions, aligning with theoretical predictions, and generalizes efficiently to multidimensional scenarios. The performance of the algorithm is compared with that of tensor cross-interpolation (TCI) and multiscale interpolative constructions through a comprehensive comparative study. When function evaluation is inexpensive or when the function is not analytical, TCI is generally more efficient for function loading. However, the proposed method shows competitive performance, outperforming TCI in certain multivariate scenarios. Moreover, it shows advantageous scaling rates and generalizes to a wider range of tasks by providing a framework for function composition in MPS, which is useful for non-linear problems and many-body statistical physics.

The Quantum Approximate Optimization Algorithm (QAOA) is an algorithm originally proposed to find approximate solutions to Combinatorial Optimization problems on quantum computers. However, the algorithm has also attracted interest for sampling purposes since it was theoretically demonstrated under reasonable complexity assumptions that one layer of the algorithm already engineers a probability distribution beyond what can be simulated by classical computers. In this regard, a recent study has shown as well that, in universal Ising models, this global probability distribution resembles pure but thermal-like distributions at a temperature that depends on internal correlations of the spin model. In this work, through an interferometric interpretation of the algorithm, we extend the theoretical derivation of the amplitudes of the eigenstates, and the Boltzmann distributions generated by single-layer QAOA. We also review the implications that this behavior has from both a practical and fundamental perspective.

In the out-of-equilibrium evolution induced by a quench, fast degrees of freedom generate long-range entanglement that is hard to encode with standard tensor networks. However, local observables only sense such long-range correlations through their contribution to the reduced local state as a mixture. We present a tensor network method that identifies such long-range entanglement and efficiently transforms it into mixture, much easier to represent. In this way, we obtain an effective description of the time-evolved state as a density matrix that captures the long-time behavior of local operators with finite computational resources.

To address the outstanding task of detecting entanglement in large quantum systems, entanglement witnesses have emerged, addressing the separable nature of a state. Yet optimizing witnesses, or accessing them experimentally, often remains a challenge. We here introduce a family of entanglement witnesses for open quantum systems, based on the electric field — its quadratures and the total fluorescence. More general than spin-squeezing inequalities, it can detect new classes of entangled states, as changing the direction for far-field observation opens up a continuous family of witnesses, without the need for a state tomography. Their efficiency is demonstrated by detecting, from almost any direction, the entanglement of collective single-photon states, such as long-lived states generated by cooperative spontaneous emission. Able to detect entanglement in large quantum systems, these electric-field-based witnesses can be used on any set of emitters described by the Pauli group, such as atomic systems (cold atoms and trapped ions), giant atoms, color centers, and superconducting qubits.

The generation of entanglement across different nodes in distributed quantum architectures plays a pivotal role for different applications. In particular, deterministic, robust, and fast protocols that prepare genuine multipartite entangled states are highly desirable. In this article, we propose a fractional quantum state transfer, in which the excitation of an emitter is partially transmitted through the quantum communication channel and then absorbed at a spatially separated node. This protocol is based on wavepacket shaping allowing for a fast deterministic generation of Bell states among two quantum registers and $W$ states for a general setting of $N$ qubits, either in a sequential or simultaneous fashion, depending on the topology of the network. By means of detailed numerical simulations, we show that genuine multipartite entangled states can be faithfully prepared within current experimental platforms and discuss the role of the main decoherence sources, qubit dephasing and relaxation, depending on the network topology.

Two non-identical quantum emitters, when placed within a cavity and coherently excited at the two-photon resonance, can reach stationary states of nearly maximal entanglement. In Vivas-Via\~na et al., we introduce a frequency-resolved Purcell effect stabilizing entangled $W$ states among strongly interacting quantum emitters embedded in a cavity. Here, we delve deeper into a specific configuration with a particularly rich phenomenology: two interacting quantum emitters under coherent excitation at the two-photon resonance. This scenario yields two resonant cavity frequencies where the combination of two-photon driving and Purcell-enhanced decay stabilizes the system into the sub- and superradiant states, respectively. By considering the case of non-degenerate emitters and exploring the parameter space of the system, we show that this mechanism is merely one among a complex family of phenomena that can generate both stationary and metastable entanglement when driving the emitters at the two-photon resonance. We provide a global perspective of this landscape of mechanisms and contribute analytical characterizations and insights into these phenomena, establishing connections with previous reports in the literature and discussing how some of these effects can be optically detected.

Anomalous topological phases, where edge states coexist with topologically trivial Chern bands, can only appear in periodically driven lattices. When the driving is smooth and continuous, the bulk-edge correspondence is guaranteed by the existence of a bulk invariant known as the winding number. However, in lattices subject to periodic discrete-step walks the existence of edge states does not only depend on bulk invariants but also on the boundary. This is a consequence of the absence of an intrinsic time-dependence or micromotion in discrete-step walks. We report the observation of edge states and a simultaneous measurement of the bulk invariants in anomalous topological phases in a two-dimensional synthetic photonic lattice subject to discrete-step walks. The lattice is implemented using time multiplexing of light pulses in two coupled fibre rings, in which one of the dimensions displays real space dynamics and the other one is parametric. The presence of edge states is inherent to the periodic driving and depends on the properties of the boundary in the implemented two-band model with zero Chern number. We provide a suitable expression for the topological invariants whose calculation does not rely on micromotion dynamics.

We introduce a driven-dissipative Bose-Hubbard chain describing coupled lossy photonic modes, in which time-reversal symmetry is broken by a coherent drive with a uniform phase gradient. We investigate this model by means of a Gaussian variational ansatz and numerically prove that the steady-state solution is stabilized by an inhomogeneous profile of the driving amplitude, which damps out boundary effects. Our calculations unveil a non-equilibrium phase diagram showing low- and high-density phases for photons separated by a phase coexistence region in which the system exhibits the phenomenon of topological amplification and is characterized by a finite non-Hermitian winding number. Our work shows the emergence of non-Hermitian topological phases in an interacting model that can be naturally implemented with superconducting circuits.

We propose to use the quantum states of an electron trapped on the inner surface of a graphene nanotorus to realize as a new kind of physical quantum bit, which can be used to encode quantum information. Fundamental tasks for quantum information processing, such as the qubit initialization and the implementation of arbitrary single qubit gates, can then be performed using external magnetic and electric fields. We also analyze the robustness of the device again systematic errors, which can be suppressed by a suitable choice of the external control fields. These findings open new prospects for the development an alternative platform for quantum computing, the scalability of which remains to be determined.

We study the spin-glass transition in several Ising models of relevance for quantum annealers. We extract the spin-glass critical temperature by extrapolating the pseudo-critical properties obtained with Replica-Exchange Monte-Carlo for finite-size systems. We find a spin-glass phase for some random lattices (random-regular and small-world graphs) in good agreement with previous results. However, our results for the quasi-two-dimensional graphs implemented in the D-Wave annealers (Chimera, Zephyr, and Pegasus) indicate only a zero-temperature spin-glass state, as their pseudo-critical temperature drifts towards smaller values. This implies that the asymptotic runtime to find the low-energy configuration of those graphs is likely to be polynomial in system size, nevertheless, this scaling may only be reached for very large system sizes — much larger than existing annealers — as we observe an abrupt increase in the computational cost of the simulations around the pseudo-critical temperatures. Thus, two-dimensional systems with local crossings can display enough complexity to make unfeasible the search with classical methods of low-energy configurations.

In this manuscript, we explore the intersection of QML and TN in the context of the one-dimensional ANNNI model with a transverse field. The study aims to concretely connect QML and TN by combining them in various stages of algorithm construction, focusing on phase diagram reconstruction for the ANNNI model, with supervised and unsupervised techniques. The model’s significance lies in its representation of quantum fluctuations and frustrated exchange interactions, making it a paradigm for studying magnetic ordering, frustration, and the presence of a floating phase. It concludes with discussions of the results, including insights from increased system sizes and considerations for future work, such as addressing limitations in QCNN and exploring more realistic implementations of QC.

We report a driven-dissipative mechanism to generate stationary entangled $W$ states among strongly-interacting quantum emitters placed within a cavity. Driving the ensemble into the highest energy state — whether coherently or incoherently — enables a subsequent cavity-enhanced decay into an entangled steady state consisting of a single de-excitation shared coherently among all emitters, i.e., a $W$ state, well known for its robustness against qubit loss. The non-harmonic energy structure of the interacting ensemble allows this transition to be resonantly selected by the cavity, while quenching subsequent off-resonant decays. Evidence of this purely dissipative mechanism should be observable in state-of-the-art cavity QED systems in the solid-state, enabling new prospects for the scalable stabilization of quantum states in dissipative quantum platforms.

We investigate a quasi-two-dimensional system consisting of two species of alkali atoms confined in a specific optical lattice potential [Phys. Rev. A 95, 053608 (2017)]. In the low-energy regime, this system is governed by a unique $\mathbb{Z}_N$ gauge theory, where field theory arguments have suggested that it may exhibit two exotic gapless deconfined phases, namely a dipolar liquid phase and a Bose liquid phase, along with two gapped (confined and deconfined) phases. We address these predictions numerically by using large-scale density matrix renormalization group simulations. Our findings provide conclusive evidence for the existence of a gapless Bose liquid phase for $N \geq 7$. We demonstrate that this gapless phase shares the same critical properties as one-dimensional critical phases, resembling weakly coupled chains of Luttinger liquids. In the range of ladder and cylinder geometries and $N$ considered, the gapless dipolar phase predicted theoretically is still elusive and its characterization will probably require a full two-dimensional treatment.

Poor man’s Majorana Bound States (MBS) arise in minimal Kitaev chains when the parameters are fine-tuned to a sweet spot. We consider an interacting two-site Kitaev chain coupled to a single-mode cavity and show that the sweet spot condition can be controlled with the cavity frequency and the hopping between sites. Furthermore, we demonstrate that photon-mediated effective interactions can be used to screen intrinsic interactions, improving the original quality of the MBS. We describe experimental signatures in the cavity transmission to detect their presence and quality. Our work proposes a new way to tune poor man’s MBS in a quantum dot array coupled to a cavity.

We introduce a variational algorithm based on Matrix Product States that is trained by minimizing a generalized free energy defined using Tsallis entropy instead of the standard Gibbs entropy. As a result, our model can generate the probability distributions associated with generalized statistical mechanics. The resulting model can be efficiently trained, since the resulting free energy and its gradient can be calculated exactly through tensor network contractions, as opposed to standard methods which require estimating the Gibbs entropy by sampling. We devise a variational annealing scheme by ramping up the inverse temperature, which allows us to train the model while avoiding getting trapped in local minima. We show the validity of our approach in Ising spin-glass problems by comparing it to exact numerical results and quasi-exact analytical approximations. Our work opens up new possibilities for studying generalized statistical physics and solving combinatorial optimization problems with tensor networks.

We propose a novel experimental protocol to measure generalized temporal entropies in many-body quantum systems. Our approach involves using local operators as probes to characterize the out-of-equilibrium dynamics induced by a geometric double quench on a replicated system. Such protocol mimics the path-integral on the corresponding Riemann surface encoding generalized temporal entanglement. We present the results of tensor network simulations of one-dimensional systems which validate the protocol and demonstrate the experimental feasibility of measuring generalized temporal entropies, and we outline the experimental requirements for implementing these quenches using state-of-the-art quantum simulators. Therefore, our results provide a physical interpretation of the meaning of generalized temporal entropies. Furthermore, they reveal that the dynamics induced on two replicas of the Ising model in a transverse field differ qualitatively from the ones of its non-integrable extension, suggesting that generalized temporal entropies can be used as a tool for identifying different dynamical classes in quantum systems.

Noisy Intermediate-Scale Quantum (NISQ) devices lack error correction, limiting scalability for quantum algorithms. In this context, digital-analog quantum computing (DAQC) offers a more resilient alternative quantum computing paradigm that outperforms digital quantum computation by combining the flexibility of single-qubit gates with the robustness of analog simulations. This work explores the impact of noise on both digital and DAQC paradigms and demonstrates DAQC’s effectiveness in error mitigation. We compare the quantum Fourier transform and quantum phase estimation algorithms under a wide range of single and two-qubit noise sources in superconducting processors. DAQC consistently surpasses digital approaches in fidelity, particularly as processor size increases. Moreover, zero-noise extrapolation further enhances DAQC by mitigating decoherence and intrinsic errors, achieving fidelities above 0.95 for 8 qubits, and reducing computation errors to the order of $10^{-3}$. These results establish DAQC as a viable alternative for quantum computing in the NISQ era.

Using topology optimization, we inverse-design nanophotonic cavities enabling the preparation of pure states of pairs and triples of quantum emitters. Our devices involve moderate values of the dielectric constant, operate under continuous laser driving, and yield fidelities to the target (Bell and W) states approaching unity for distant qubits (several natural wavelengths apart). In the fidelity optimization procedure, our algorithm generates entanglement by maximizing the dissipative coupling between the emitters, which allows the formation of multipartite pure steady states in the driven-dissipative dynamics of the system. Our findings open the way towards the efficient and fast preparation of multiqubit quantum states with engineered features, with potential applications for nonclassical light generation, quantum simulation, and quantum sensing.

Recent experimental work has demonstrated the ability to achieve reconfigurable photon localization in lossy photonic lattices by continuously driving them with lasers strategically positioned at specific locations. This localization results from the perfect, destructive interference of light emitted from different positions and, because of that, occurs only at very specific frequencies. Here, we examine this localization regime in the presence of standard optical Kerr non-linearities, such as those found in polaritonic lattices, and show that they stabilize driven-dissipative localization across frequency ranges significantly broader than those observed in the linear regime. Moreover, we demonstrate that, contrary to intuition, in most siutations this driven-dissipative localization does not enhance non-linear effects like optical bistabilities, due to a concurrent reduction in overall intensities. Nevertheless, we are able to identify certain parameter regions where non-linear enhancement is achieved, corresponding to situations where emission from different spots constructively interferes.

Interference underpins some of the most practical and impactful properties of both the classical and quantum worlds. In this work we experimentally investigate a new formalism to describe interference effects, based on collective states which have enhanced or suppressed coupling to a two-level system. We employ a single trapped ion, whose electronic state is coupled to two of the ion’s motional modes in order to simulate a multi-mode light-matter interaction. We observe the emergence of phononic bright and dark states for both a single phonon and a superposition of coherent states and demonstrate that a view of interference which is based solely on their decomposition in the collective basis is able to intuitively describe their coupling to a single atom. This work also marks the first time that multi-mode bright and dark states have been formed with the bounded motion of a single trapped ion and we highlight the potential of the methods discussed here for use in quantum information processing.

Photonic quantum metrology harnesses quantum states of light, such as NOON or Twin-Fock states, to measure unknown parameters beyond classical precision limits. Current protocols suffer from two severe limitations that preclude their scalability: the exponential decrease in fidelities (or probabilities) when generating states with large photon numbers due to gate errors, and the increased sensitivity of such states to noise. Here, we develop a deterministic protocol combining quantum optical non-linearities and variational quantum algorithms that provides a substantial improvement on both fronts. First, we show how the variational protocol can generate metrologically-relevant states with a small number of operations which does not significantly depend on photon-number, resulting in exponential improvements in fidelities when gate errors are considered. Second, we show that such states offer a better robustness to noise compared to other states in the literature. Since our protocol harnesses interactions already appearing in state-of-the-art setups, such as cavity QED, we expect that it will lead to more scalable photonic quantum metrology in the near future.

We introduce a counter-diabatic approach for deriving Hamiltonians modeling superchargable quantum batteries (QBs). A necessary requirement for the supercharging process is the existence of multipartite interactions among the cells of the battery. Remarkably, this condition may be insufficient no matter the number of multipartite terms in the Hamiltonian. We analytically illustrate this kind of insufficiency through a model of QB based on the adiabatic version for the Grover search problem. On the other hand, we provide QB supercharging with just a mild number of global connections in the system. To this aim, we consider a spin-$1/2$ chain with $n$ sites in the presence of Ising multipartite interactions. We then show that, by considering the validity of the adiabatic approximation and by adding $n$ terms of $(n-1)$-site interactions, we can achieve a Hamiltonian exhibiting maximum QB power, with respect to a normalized evolution time, growing quadratically with $n$. Therefore, supercharging can be achieved by $O(n)$ terms of multipartite connections. The time constraint required by the adiabatic approximation can be surpassed by considering a counter-diabatic expansion in terms of the gauge potential for the original Hamiltonian, with a limited $O(n)$ many-body interaction terms assured via a Floquet approach for the counter-diabatic implementation.

Quantum obesity (QO) is new function used to quantify quantum correlations beyond entanglement, which also works as a witness for entanglement. Thanks to its analyticity for arbitrary state of bipartite systems, it represents an advantage with respect to other quantum correlations, like quantum discord for example. In this work we show that QO is a fundamental quantity to observe signature of quantum phase transitions. We also describe a mechanism based on local filtering operations able to intensify the critical behavior of the QO near to the transition point. To this end, we introduce a theorem stating how QO changes under local quantum operations and classical communications. This work opens perspective for the characterization of new phenomena in quantum critical systems through the analytically computable pairwise QO.

We propose an efficient, scalable, and deterministic scheme to generate multiple indistinguishable photons over independent channels, on demand. Our design relies on multiple single-photon sources, each coupled to a waveguide, and all of them interact with a common cavity mode. The cavity synchronizes and triggers the simultaneous emission of one photon by each source, which are collected by the waveguides. For a state-of-the-art circuit QED implementation, this scheme supports the creation of single photons with purity, indistinguishability, and efficiency of $99\%$ at rates of $\sim $MHz. We also discuss conditions to produce up to 100 photons simultaneously with generation rates of hundreds of kHz. This is orders of magnitude more efficient than previous demultiplexed sources for boson sampling and enables the realization of deterministic multi-photon sources and scalable quantum information processing with photons.

The presence of gauge symmetry in 1+1D is known to be redundant, since it does not imply the existence of dynamical gauge bosons. As a consequence, in the continuum, the Abelian-Higgs model, the theory of bosonic matter interacting with photons, just possesses a single phase, as the higher dimensional Higgs and Coulomb phases are connected via non-perturbative effects. However, recent research published in [Phys. Rev. Lett. 128, 090601 (2022)] has revealed an unexpected phase transition when the system is discretized on the lattice. This transition is described by a conformal field theory with a central charge of $c=3/2$. In this paper, we aim to characterize the two components of this $c=3/2$ theory — namely the free Majorana fermionic and bosonic parts — through equilibrium and out-of-equilibrium spectral analyses.

We propose a novel quantum battery realized with a few interacting particles in a three-well system with different on-site energies, which could be realized with ultracold atom platforms. We prepare the initial state in the lowest energy well and charge the battery using a Spatial Adiabatic Passage (SAP)-based protocol, enabling the population of a higher energy well. We examine the charging under varying interaction strengths and reveal that the consideration of collective charging results in an intriguing oscillatory behavior of the final charge for finite interactions, through diabatic evolution. Our findings open a new avenue for building stable and controllable quantum batteries.

We study the computational complexity of simulating the time-dependent expectation value of a local operator in a one-dimensional quantum system by using temporal matrix product states. We argue that such cost is intimately related to that of encoding temporal transition matrices and their partial traces. In particular, we show that we can upper-bound the rank of these reduced transition matrices by the one of the Heisenberg evolution of local operators, thus making connection between two apparently different quantities, the temporal entanglement and the local operator entanglement. As a result, whenever the local operator entanglement grows slower than linearly in time, we show that computing time-dependent expectation values of local operators using temporal matrix product states is likely advantageous with respect to computing the same quantities using standard matrix product states techniques.

Quantum systems are open, continually exchanging energy and information with the surrounding environment. This interaction leads to decoherence and decay of quantum states. In complex systems, formed by many particles, decay can become correlated and enhanced. A fundamental question then arises: what is the maximal decay rate of a large quantum system, and how does it scale with its size? In this work, we address these issues by reformulating the problem into finding the ground state energy of a generic spin Hamiltonian. Inspired by recent work in Hamiltonian complexity theory, we establish rigorous and general upper and lower bounds on the maximal decay rate. These bounds are universal, as they hold for a broad class of Markovian many-body quantum systems. For many physically-relevant systems, the bounds are asymptotically tight, resulting in exact scaling laws with system size. Specifically, for large atomic arrays in free space, these scalings depend only on the arrays’ dimensionality and are insensitive to details at short length-scales. The scaling laws set fundamental limits on the decay rates of all quantum states, shed light on the behavior of generic driven-dissipative systems, and may ultimately constrain the scalability of quantum processors and simulators based on atom arrays.