Publications of Tomás Ramos

Whether long-range interactions allow for a form of causality in non-relativistic quantum models remains an open question with far-reaching implications for the propagation of information and thermalization processes. Here, we study the out-of-equilibrium dynamics of the one-dimensional transverse Ising model with algebraic long-range exchange coupling. Using a state of the art tensor-network approach, complemented by analytic calculations and considering various observables, we show that a weak form of causality emerges, characterized by non-universal dynamical exponents. While the local spin and spin correlation causal edges are sub-ballistic, the causal region has a rich internal structure, which, depending on the observable, displays ballistic or super-ballistic features. In contrast, the causal region of entanglement entropy is featureless and its edge is always ballistic, irrespective of the interaction range. Our results shed light on the propagation of information in long-range interacting lattice models and pave the way to future experiments, which are discussed.

Topologically-ordered phases of matter, although stable against local perturbations, are usually restricted to relatively small regions in phase diagrams. Their preparation requires thus a precise fine tunning of the system’s parameters, a very challenging task in most experimental setups. In this work, we investigate a model of spinless fermions interacting with dynamical $\mathbb{Z}_2$ gauge fields on a cross-linked ladder, and show evidence of topological order throughout the full parameter space. In particular, we show how a magnetic flux is spontaneously generated through the ladder due to an Aharonov-Bohm instability, giving rise to topological order even in the absence of a plaquette term. Moreover, the latter coexists here with a symmetry-protected topological phase in the matter sector, that displays fractionalised gauge-matter edge states, and intertwines with it by a flux-threading phenomenon. Finally, we unveil the robustness of these features through a gauge frustration mechanism, akin to geometric frustration in spin liquids, allowing topological order to survive to arbitrarily large quantum fluctuations. In particular, we show how, at finite chemical potential, topological solitons are created in the gauge field configuration, which bound to fermions forming $\mathbb{Z}_2$ deconfined quasi-particles. The simplicity of the model makes it an ideal candidate where 2D gauge theory phenomena, as well as exotic topological effects, can be investigated using cold-atom quantum simulators.

Lattice gauge theories, which originated from particle physics in the context of Quantum Chromodynamics (QCD), provide an important intellectual stimulus to further develop quantum information technologies. While one long-term goal is the reliable quantum simulation of currently intractable aspects of QCD itself, lattice gauge theories also play an important role in condensed matter physics and in quantum information science. In this way, lattice gauge theories provide both motivation and a framework for interdisciplinary research towards the development of special purpose digital and analog quantum simulators, and ultimately of scalable universal quantum computers. In this manuscript, recent results and new tools from a quantum science approach to study lattice gauge theories are reviewed. Two new complementary approaches are discussed: first, tensor network methods are presented – a classical simulation approach – applied to the study of lattice gauge theories together with some results on Abelian and non-Abelian lattice gauge theories. Then, recent proposals for the implementation of lattice gauge theory quantum simulators in different quantum hardware are reported, e.g., trapped ions, Rydberg atoms, and superconducting circuits. Finally, the first proof-of-principle trapped ions experimental quantum simulations of the Schwinger model are reviewed.

We consider the time evolution of the gaps of the entanglement spectrum for a block of consecutive sites in finite transverse field Ising chains after sudden quenches of the magnetic field. We provide numerical evidence that, whenever we quench at or across the quantum critical point, the time evolution of the ratios of these gaps allows to obtain universal information. They encode the low-lying gaps of the conformal spectrum of the Ising boundary conformal field theory describing the spatial bipartition within the imaginary time path integral approach to global quenches at the quantum critical point.

We benchmark the accuracy of a variational quantum eigensolver based on a finite-depth quantum circuit encoding ground state of local Hamiltonians. We show that in gapped phases, the accuracy improves exponentially with the depth of the circuit. When trying to encode the ground state of conformally invariant Hamiltonians, we observe two regimes. A finite-depth regime, where the accuracy improves slowly with the number of layers, and a finite-size regime where it improves again exponentially. The cross-over between the two regimes happens at a critical number of layers whose value increases linearly with the size of the system. We discuss the implication of these observations in the context of comparing different variational ansatz and their effectiveness in describing critical ground states.

We excite the vacuum of a relativistic theory of bosons coupled to a $U(1)$ gauge field in 1+1 dimensions (bosonic Schwinger model) out of equilibrium by creating a spatially separated particle-antiparticle pair connected by a string of electric field. During the evolution, we observe a strong confinement of bosons witnessed by the bending of their light cone, reminiscent of what has been observed for the Ising model [Nat. Phys. 13, 246 (2017)]. As a consequence, for the time scales we are able to simulate, the system evades thermalization and generates exotic asymptotic states. These states are made of two disjoint regions, an external deconfined region that seems to thermalize, and an inner core that reveals an area-law saturation of the entanglement entropy.

Tensor network (TN), a young mathematical tool of high vitality and great potential, has been undergoing extremely rapid developments in the last two decades, gaining tremendous success in condensed matter physics, atomic physics, quantum information science, statistical physics, and so on. In this lecture notes, we focus on the contraction algorithms of TN as well as some of the applications to the simulations of quantum many-body systems. Starting from basic concepts and definitions, we first explain the relations between TN and physical problems, including the TN representations of classical partition functions, quantum many-body states (by matrix product state, tree TN, and projected entangled pair state), time evolution simulations, etc. These problems, which are challenging to solve, can be transformed to TN contraction problems. We present then several paradigm algorithms based on the ideas of the numerical renormalization group and/or boundary states, including density matrix renormalization group, time-evolving block decimation, coarse-graining/corner tensor renormalization group, and several distinguished variational algorithms. Finally, we revisit the TN approaches from the perspective of multi-linear algebra (also known as tensor algebra or tensor decompositions) and quantum simulation. Despite the apparent differences in the ideas and strategies of different TN algorithms, we aim at revealing the underlying relations and resemblances in order to present a systematic picture to understand the TN contraction approaches.

We study the non-equilibrium dynamics of a 1D Bose-Hubbard model in a gradient potential and a superlattice, beginning from a deep Mott insulator regime with an average filling of one particle per site. Studying a quench that is near resonance to tunnelling of the particles over two lattice sites, we show how a spin model emerges consisting of two coupled Ising chains that are coupled by interaction terms in a staggered geometry. We compare and contrast the behavior in this case with that in a previously studied case where the resonant tunnelling was over a single site. Using optimized tensor network techniques to calculate finite temperature behavior of the model, as well as finite size scaling for the ground state, we conclude that the universality class of the phase transition for the coupled chains is that of a tricritical Ising point. We also investigate the out-of-equilibrium dynamics after the quench in the vicinity of the resonance and compare dynamics with recent experiments realized without the superlattice geometry. This model is directly realizable in current experiments, and reflects a new general way to realize spin models with ultracold atoms in optical lattices.

We show that atoms in tilted optical superlattices provide a platform for exploring coupled spin chains of forms that are not present in other systems. In particular, using a period-2 superlattice in 1D, we show that coupled Ising spin chains with XZ and ZZ spin coupling terms can be engineered. We use optimized tensor network techniques to explore the criticality and non-equilibrium dynamics in these models, finding a tricritical Ising point in regimes that are accessible in current experiments. These setups are ideal for studying low-entropy physics, as initial entropy is “frozen-out” in realizing the spin models, and provide an example of the complex critical behaviour that can arise from interaction-projected models.

We propose and analyze a mechanism for rectification of spin transport through a small junction between two spin baths or leads. For interacting baths we show that transport is conditioned on the spacial asymmetry of the quantum junction mediating the transport, and attribute this behavior to a gapped spectral structure of the lead-system-lead configuration. For non-interacting leads a minimal quantum model that allows for spin rectification requires an interface of only two interacting two-level systems. We obtain approximate results with a weak-coupling Born-master-equation in excellent agreement with matrix-product-state calculations that are extrapolated in time by mimicking absorbing boundary conditions. These results should be observable in controlled spin systems realized with cold atoms, trapped ions, or in electrons in quantum dot arrays.

The fact that the computational cost of simulating a many-body quantum system on a computer increases with the amount of entanglement has been considered as the major bottleneck for simulating its out-of-equilibrium dynamics. Some aspects of the dynamics are, nevertheless, robust under appropriately devised approximations. Here we present a possible algorithm that allows to systematically approximate the equilibration value of local operators after a quantum quench. At the core of our proposal there is the idea to transform entanglement between distant parts of the system into mixture, and at the same time preserving the local reduced density matrices of the system. We benchmark the resulting algorithm by studying quenches of quadratic Fermionic Hamiltonians.

We show how to accurately study 2D quantum critical phenomena using infinite projected entangled-pair states (iPEPS). We identify the presence of a finite correlation length in the optimal iPEPS approximation to Lorentz-invariant critical states which we use to perform a finite correlation-length scaling (FCLS) analysis to determine critical exponents. This is analogous to the one-dimensional (1D) finite entanglement scaling with infinite matrix product states. We provide arguments why this approach is also valid in 2D by identifying a class of states that despite obeying the area law of entanglement seems hard to describe with iPEPS. We apply these ideas to interacting spinless fermions on a honeycomb lattice and obtain critical exponents which are in agreement with Quantum Monte Carlo results. Furthermore, we introduce a new scheme to locate the critical point without the need of computing higher order moments of the order parameter. Finally, we also show how to obtain an improved estimate of the order parameter in gapless systems, with the 2D Heisenberg model as an example.

We study the spreading of information in a wide class of quantum systems, with variable-range interactions. We show that, after a quench, it generally features a double structure, whose scaling laws are related to a set of universal microscopic exponents that we determine. When the system supports excitations with a finite maximum velocity, the spreading shows a twofold ballistic behavior. While the correlation edge spreads with a velocity equal to twice the maximum group velocity, the dominant correlation maxima propagate with a different velocity that we derive. When the maximum group velocity diverges, as realizable with long-range interactions, the correlation edge features a slower-than-ballistic motion. The motion of the maxima is, instead, either faster-than-ballistic, for gapless systems, or ballistic, for gapped systems. The phenomenology that we unveil here provides a unified framework, which encompasses existing experimental observations with ultracold atoms and ions. It also paves the way to simple extensions of those experiments to observe the structures we describe in their full generality.

Fundamental forces of Nature are described by field theories, also known as gauge theories, based on a local gauge invariance. The simplest of them is quantum electrodynamics (QED), which is an example of an Abelian gauge theory. Such theories describe the dynamics of massless photons and their coupling to matter. However, in two spatial dimension (2D) they are known to exhibit gapped phases at low temperature. In the realm of quantum spin systems, it remains a subject of considerable debate if their low energy physics can be described by emergent gauge degrees of freedom. Here we present a class of simple two-dimensional models that admit a low energy description in terms of an Abelian gauge theory. We find rich phase diagrams for these models comprising exotic deconfined phases and gapless phases – a rare example for 2D Abelian gauge theories. The counter-intuitive presence of gapless phases in 2D results from the emergence of additional symmetry in the models. Moreover, we propose schemes to realize our model with current experiments using ultracold bosonic atoms in optical lattices.

We introduce a theoretical scheme for the analog quantum simulation of long-range XYZ models using current trapped-ion technology. In order to achieve fully-tunable Heisenberg-type interactions, our proposal requires a state-dependent dipole force along a single vibrational axis, together with a combination of standard resonant and detuned carrier drivings. We discuss how this quantum simulator could explore the effect of long-range interactions on the phase diagram by combining an adiabatic protocol with the quasi-periodic drivings and test the validity of our scheme numerically. At the isotropic Heisenberg point, we show that the long-range Hamiltonian can be mapped onto a non-linear sigma model with a topological term that is responsible for its low-energy properties, and we benchmark our predictions with Matrix-Product-State numerical simulations.