In the last decade there has been an intense activity aimed at the quantum simulation of interacting many-body systems using cold atoms [1, 2]. The idea of quantum simulations traces back to Feynman [3], who argued that the ideal setting to study quantum systems would be a quantum experimental setup rather then a classical one - the latter one being fundamentally limited due to its hardware classical structure. This is a particularly important problem given the intrinsic complexity of interacting many-body problems, and the difficulties that arise when tackling them with numerical simulations - two paradigmatic examples being the sign(s) problem affecting Monte Carlo simulations of fermionic systems, and the real time dynamics in more than one spatial dimension. Ultracold atoms offer a very powerful setting for quantum simulations. Atoms can be trapped in tailored optical and magnetic potentials, also controlling their dimensionality. The inter-atomic interactions can be tuned by external knobs, such as Feshbach resonances. This gives a large freedom on model building and, with suitable mappings, they allow the implementation of desired target models. This allowed an impressive exploitation of quantum simulators on the context of condensed matter physics. The simulation of high-energy physics is an important line of research in this field and it is less direct. In particular it requires the implementation of symmetries like Lorentz and gauge invariance which are not immediately available in a cold atomic setting. Gauge fields are ubiquitous in physics ranging from condensed matter [4–6] and quantum computation [7, 8] to particle physics [9], an archetypical example being Quantum Chromodynamics (QCD) [10,11], the theory of strong nuclear forces. Currently open problems in QCD, providing a long-term goal of cold atomic simulations, include confinement/deconfinement and the structure of color superconducting phases at finite chemical potential [12]. Even though QCD is a very complicated theory (to simulate or study), it is possible to envision a path through implementation of simpler models. Furthermore, it is also expected that interesting physics is found on such “intermediate models” which may deserve attention irrespectively of the QCD study. A very relevant model in this regard is the Schwinger model (Quantum Electrodynamics in 1+1 dimensions) [13]. This theory exhibits features of QCD, such as confinement [14], and is at the very same time amenable to both theoretical studies and simpler experimental schemes. This model was the target of the first experimental realization of a gauge theory with a quantum simulator [15]. The work on this Thesis is, in part, motivated by the study of toy models which put in evidence certain aspects that can be found in QCD. Such toy models provide also intermediate steps in the path towards more complex simulations. The two main aspects of QCD which are addressed here are symmetry-locking and confinement. The other main motivation for this study is to develop a systematic framework, through dimensional mismatch, for theoretical understanding and quantum simulations of long-range theories using gauge theories. The model used to study symmetry-locking consists of a four-fermion mixture [16]. It has the basic ingredients to exhibit a non-Abelian symmetry-locked phase: the full Hamiltonian has an SU (2) × SU (2) (global) symmetry which can break to a smaller SU (2) group. Such phase is found in a extensive region of the phase diagram by using a mean-field approach and a strong coupling expansion. A possible realization of such system is provided by an Ytterbium mixture. Even without tuning interactions, it is shown that such mixture falls inside the the locked-symmetry phase pointing towards a possible realization in current day experiments. The models with dimensional mismatch investigated here have fermions in a lower dimensionality d + 1 and gauge fields in higher dimensionality D + 1. They serve two purposes: establish mappings to non-local theories by integration of fields [17] and the study of confinement [18]. In the particular case of d = 1 and D = 2 it is found that some general non-local terms can be obtained on the Lagrangian [17]. This is found in the form of power-law expansions of the Laplacian mediating either kinetic terms (for bosons) or interactions (for fermions). The fact that such expansions are not completely general is not surprising since constraints do exist, preventing unphysical features like breaking of unitarity. The non-local terms obtained are physically acceptable, in this regard, since they are derived from unitary theories. The above mapping is done exactly. In certain cases it is shown that it is possible to construct an effective long-range Hamiltonian in a perturbative expansion. In particular it is shown how this is done for non-relativistic fermions (in d = 1) and 3 + 1 gauge fields. These results are relevant in the context of state of the art experiments which implement models with long-range interactions and where theoretical results are less abundant than for the case of local theories. The above mappings establish a direct relation with local theories which allow theoretical insight onto these systems. Examples of this would consist on the application of Mermin- Wagner-Hohenberg theorem [19, 20] and Lieb-Robinson bounds [21], on the propagation of quantum correlations, to non-local models. In addition they can also provide a path towards implementation of tunable long-range interactions with cold atoms. Furthermore, in terms of quantum simulations, they are in between the full higher dimensional system and the full lower dimensional one. Such property is attractive from the point of view of a gradual increase of complexity for quantum simulations of gauge theories. Another interesting property of these models is that they allow the study of confinement beyond the simpler case of the Schwinger model. The extra dimensions are enough to attribute dynamics to the gauge field, which are no longer completely fixed by the Gauss law. The phases of the Schwinger model are shown to be robust under variation of the dimension of the gauge fields [18]. Both the screened phase, of the massless case, and the confined phase, of the massive one, are found for gauge fields in 2+1 and 3+1 dimensions. Such results are also obtained in the Schwinger- Thirring model. This shows that these phases are very robust and raises interesting questions about the nature of confinement. Robustness under Thirring interactions are relevant because it shows that errors on the experimental implementation will not spoil the phase. Even more interesting is the case of gauge fields in higher dimensions since confinement in the Schwinger model is intuitively atributed to the dimensionality of the gauge fields (creating linear potentials between particles). This Thesis is organized as follows. In Chapter 1, some essential background regarding quantum simulations of gauge theories is provided. It gives both a brief introduction to cold atomic physics and lattice gauge theories. In Chapter 2, it is presented an overview over proposals of quantum simulators of gauge potentials and gauge fields. At the end of this Chapter, in Section 2.3 it is briefly presented ongoing work on a realization of the Schwinger model that we term Half Link Schwinger model. There, it is argued, some of the generators of the gauge symmetry on the lattice can be neglected without comprimising gauge invariance. In Chapter 3, the results regarding the phase diagram of the fourfermion mixture, exhibiting symmetry-locking, is presented. In Chapter 4, the path towards controlling non-local kinetic terms and interactions is provided, after a general introduction to the formalism of dimensional mismatch. At the end the construction of effective Hamiltonians is described. Finally, Chapter 5 concerns the study of confinement and the robustness of it for 1+1 fermions. The first part regards the Schwinger-Thirring with the presence of a -term while, in the second part, models with dimensional mismatch are considered. The thesis ends with conclusions and perspectives of future work based on the results presented here.
In the last decade there has been an intense activity aimed at the quantum simulation of interacting many-body systems using cold atoms [1, 2]. The idea of quantum simulations traces back to Feynman [3], who argued that the ideal setting to study quantum systems would be a quantum experimental setup rather then a classical one - the latter one being fundamentally limited due to its hardware classical structure. This is a particularly important problem given the intrinsic complexity of interacting many-body problems, and the difficulties that arise when tackling them with numerical simulations - two paradigmatic examples being the sign(s) problem affecting Monte Carlo simulations of fermionic systems, and the real time dynamics in more than one spatial dimension. Ultracold atoms offer a very powerful setting for quantum simulations. Atoms can be trapped in tailored optical and magnetic potentials, also controlling their dimensionality. The inter-atomic interactions can be tuned by external knobs, such as Feshbach resonances. This gives a large freedom on model building and, with suitable mappings, they allow the implementation of desired target models. This allowed an impressive exploitation of quantum simulators on the context of condensed matter physics. The simulation of high-energy physics is an important line of research in this field and it is less direct. In particular it requires the implementation of symmetries like Lorentz and gauge invariance which are not immediately available in a cold atomic setting. Gauge fields are ubiquitous in physics ranging from condensed matter [4–6] and quantum computation [7, 8] to particle physics [9], an archetypical example being Quantum Chromodynamics (QCD) [10,11], the theory of strong nuclear forces. Currently open problems in QCD, providing a long-term goal of cold atomic simulations, include confinement/deconfinement and the structure of color superconducting phases at finite chemical potential [12]. Even though QCD is a very complicated theory (to simulate or study), it is possible to envision a path through implementation of simpler models. Furthermore, it is also expected that interesting physics is found on such “intermediate models” which may deserve attention irrespectively of the QCD study. A very relevant model in this regard is the Schwinger model (Quantum Electrodynamics in 1+1 dimensions) [13]. This theory exhibits features of QCD, such as confinement [14], and is at the very same time amenable to both theoretical studies and simpler experimental schemes. This model was the target of the first experimental realization of a gauge theory with a quantum simulator [15]. The work on this Thesis is, in part, motivated by the study of toy models which put in evidence certain aspects that can be found in QCD. Such toy models provide also intermediate steps in the path towards more complex simulations. The two main aspects of QCD which are addressed here are symmetry-locking and confinement. The other main motivation for this study is to develop a systematic framework, through dimensional mismatch, for theoretical understanding and quantum simulations of long-range theories using gauge theories. The model used to study symmetry-locking consists of a four-fermion mixture [16]. It has the basic ingredients to exhibit a non-Abelian symmetry-locked phase: the full Hamiltonian has an SU (2) × SU (2) (global) symmetry which can break to a smaller SU (2) group. Such phase is found in a extensive region of the phase diagram by using a mean-field approach and a strong coupling expansion. A possible realization of such system is provided by an Ytterbium mixture. Even without tuning interactions, it is shown that such mixture falls inside the the locked-symmetry phase pointing towards a possible realization in current day experiments. The models with dimensional mismatch investigated here have fermions in a lower dimensionality d + 1 and gauge fields in higher dimensionality D + 1. They serve two purposes: establish mappings to non-local theories by integration of fields [17] and the study of confinement [18]. In the particular case of d = 1 and D = 2 it is found that some general non-local terms can be obtained on the Lagrangian [17]. This is found in the form of power-law expansions of the Laplacian mediating either kinetic terms (for bosons) or interactions (for fermions). The fact that such expansions are not completely general is not surprising since constraints do exist, preventing unphysical features like breaking of unitarity. The non-local terms obtained are physically acceptable, in this regard, since they are derived from unitary theories. The above mapping is done exactly. In certain cases it is shown that it is possible to construct an effective long-range Hamiltonian in a perturbative expansion. In particular it is shown how this is done for non-relativistic fermions (in d = 1) and 3 + 1 gauge fields. These results are relevant in the context of state of the art experiments which implement models with long-range interactions and where theoretical results are less abundant than for the case of local theories. The above mappings establish a direct relation with local theories which allow theoretical insight onto these systems. Examples of this would consist on the application of Mermin- Wagner-Hohenberg theorem [19, 20] and Lieb-Robinson bounds [21], on the propagation of quantum correlations, to non-local models. In addition they can also provide a path towards implementation of tunable long-range interactions with cold atoms. Furthermore, in terms of quantum simulations, they are in between the full higher dimensional system and the full lower dimensional one. Such property is attractive from the point of view of a gradual increase of complexity for quantum simulations of gauge theories. Another interesting property of these models is that they allow the study of confinement beyond the simpler case of the Schwinger model. The extra dimensions are enough to attribute dynamics to the gauge field, which are no longer completely fixed by the Gauss law. The phases of the Schwinger model are shown to be robust under variation of the dimension of the gauge fields [18]. Both the screened phase, of the massless case, and the confined phase, of the massive one, are found for gauge fields in 2+1 and 3+1 dimensions. Such results are also obtained in the Schwinger- Thirring model. This shows that these phases are very robust and raises interesting questions about the nature of confinement. Robustness under Thirring interactions are relevant because it shows that errors on the experimental implementation will not spoil the phase. Even more interesting is the case of gauge fields in higher dimensions since confinement in the Schwinger model is intuitively atributed to the dimensionality of the gauge fields (creating linear potentials between particles). This Thesis is organized as follows. In Chapter 1, some essential background regarding quantum simulations of gauge theories is provided. It gives both a brief introduction to cold atomic physics and lattice gauge theories. In Chapter 2, it is presented an overview over proposals of quantum simulators of gauge potentials and gauge fields. At the end of this Chapter, in Section 2.3 it is briefly presented ongoing work on a realization of the Schwinger model that we term Half Link Schwinger model. There, it is argued, some of the generators of the gauge symmetry on the lattice can be neglected without comprimising gauge invariance. In Chapter 3, the results regarding the phase diagram of the fourfermion mixture, exhibiting symmetry-locking, is presented. In Chapter 4, the path towards controlling non-local kinetic terms and interactions is provided, after a general introduction to the formalism of dimensional mismatch. At the end the construction of effective Hamiltonians is described. Finally, Chapter 5 concerns the study of confinement and the robustness of it for 1+1 fermions. The first part regards the Schwinger-Thirring with the presence of a -term while, in the second part, models with dimensional mismatch are considered. The thesis ends with conclusions and perspectives of future work based on the results presented here.
Gauge theories and topological phases play a fundamental role in different areas of physics. The first ones are at the basis of the Standard Model in the field of particle physics, describing the electroweak and strong interactions through a non-Abelian gauge theory. In condensed matter and statistical physics, gauge theories arise as low-energy effective descriptions of strongly correlated phenomena, such as quantum spin liquids and quantum Hall effect. In this realm there is a strong connection with topological phases and order, as emergent gauge fermions and bosons often describe collective excitations of new exotic states of matter of spin models. The discretization on a lattice is a possible way of dealing with the strongly coupled nature of these theories, due to the fact that this formulation at finite volume provides natural regularizing cut-offs, i.e. the lattice size and spacing. This allows for the investigation of different non-perturbative properties both numerically and analytically. Despite the success of these methods, there are various aspects which remain intractable due to the sign or complex action problems, like the out of equilibrium real time evolution or the analysis of quantum chromodynamics with finite chemical potential. In this respect, quantum simulators of many-body systems to simulate high energy physics arise as promising alternatives to face these problems in the near future. They are quantum systems that can be controlled and used to simulate more complicated systems, whose properties could not be analysed with classical computational, experimental or theoretical tools. In the last decades, there has been a huge development in the fields of quantum optics and atomic physics, allowing for the realization highly precise and controllable platforms by means of trapped ions, superconducting circuits, Rydberg atoms and ultracold atoms in optical lattices. This work of Thesis is part of these fields and has different purposes, all of them connected with the study of gauge theories and topological phases. Firstly, we want to develop a reformulation of lattice gauge theories in terms of gauge invariant fields, in a way to deal solely with physical variables directly in the action. Among several possible advantages, a crucial point is that this can be particularly helpful for the construction of consistent approximation schemes, such as mean-field theories, in order to understand and capture some of the physical features of the theory and make the mean-field approximation consistent with the Elitzur theorem, stating that a local gauge symmetry can not be spontaneously broken. As a second point, we face the problem of simulating higher dimensional gauge theories with ultracold atoms. One of the challenges in more than one dimension is indeed the realization of plaquette interaction terms in the Hamiltonian of lattice gauge theories: these can be engineered through four correlated hoppings in perturbation theory, or by means of constrained hoppings in the dual formulation. In this respect, our main target is to set up an ultracold atomic platform generating the plaquette term in two dimensions using only two correlated hoppings, and protecting gauge invariance through angular momentum conservation. The last point we address in the Thesis is related to the analysis of particles moving in static background gauge potentials, i.e. when the considered gauge field has no dynamical term in the action of the system. Indeed, the dynamics of quantum particles in the presence of static gauge fields gives rise to intriguing physical phenomena. In particular, using ultracold atom setups, the realization of artificial gauge potentials can be used to investigate the physics of topological semimetals, such as Weyl or Dirac type. In the last ten years, a lot of attention has been paid to their characterization, due to the appearence of clear theoretical predictions and very well-controlled experimental techniques. We investigate here the relation between topological phase transitions and van Hove singularities, i.e. the discontinuities in the energy derivative of the density of states, in three-dimensional gapless systems. In such materials, topological phase transitions can be defined by changes of the topological invariants of the Fermi sheets, happening at specific singular points. Moreover, Fermi surface singularities result in the presence of the so-called van Hove points. We then present a general argument to relate topological phase transitions and van Hove singularities, and show observable consequences that are related to the transport properties of the system. We exemplify our argument in Weyl systems by analyzing the three-dimensional Hofstadter model for various commensurate fluxes, which offers the opportunity to consider different kinds of Weyl metals and to understand the features of their density of states.
Can high-energy physics be simulated by low-energy, non-relativistic, many-body systems such as ultracold atoms? Such ultracold atomic systems lack the type of symmetries and dynamical properties of high energy physics models: in particular, they manifest neither local gauge invariance nor Lorentz invariance, which are crucial properties of the quantum field theories which are the building blocks of the standard model of elementary particles. However, it turns out, surprisingly, that there are ways to configure an atomic system to manifest both local gauge invariance and Lorentz invariance. In particular, local gauge invariance can arise either as an effective low-energy symmetry, or as an exact symmetry, following from the conservation laws in atomic interactions. Hence, one could hope that such quantum simulators may lead to a new type of (table-top) experiments which will be used to study various QCD (quantum chromodynamics) phenomena, such as the confinement of dynamical quarks, phase transitions and other effects, which are inaccessible using the currently known computational methods. In this report, we review the Hamiltonian formulation of lattice gauge theories, and then describe our recent progress in constructing the quantum simulation of Abelian and non-Abelian lattice gauge theories in 1 + 1 and 2 + 1 dimensions using ultracold atoms in optical lattices.
Gauge theories are ubiquitous in physics. Many intriguing phenomena in condensed matter physics owe to the action of the electromagnetic field, which is an Abelian gauge theory. The numerical treatment of many-body systems is inherently complex due to the exponentially growing size of the Hilbert space. While in one dimension an area law guarantees that numerical methods on classical computers can deal with strongly correlated systems, in higher dimensions the quantum simulation comes as the panacea for the many-body problem. The present thesis comprises the elaboration of experimentally feasible methods for the quantum simulation of dynamical Abelian gauge fields with ultra-cold gases of neutral atoms and the theoretical analysis of the related model Hamiltonians. As neutral atoms do not interact with external vector potentials like charged particles would do, the gauge fields have to be artificially engineered. The elements of a gauge theory that need to be replicated on a quantum simulator vary depending on the subject of investigation. The key ingredient at the root of many condensed matter phenomena, from the quantum Hall effect to superconductivity and chiral topological insulators, is the Berry phase. Whilst artificial static gauge fields have been widely explored, much remains to do regarding the realization of artificial dynamical gauge fields. In Chapter 3 we present a method based on the amplitude modulation of a one-dimensional optical lattice, which allows for an unprecedented degree of control over a wide range of parameters. The method also comprises the generation of a density-dependent complex phase, fundamental to the creation of anyonic pseudo-particles. The anyons are amenable of observation through interferometric measurement, realizable with the same experimental set-up. With regard to gauge theories, the Berry phase is just the visible tip of the iceberg. Below the waterline, there is more to consider in order to comprehensively reproduce a gauge theory, like the electric and magnetic fields in quantum electrodynamics. Moreover, a full account for the inherent symmetry is crucial to investigate phenomena proper of non-Abelian gauge theories in the context of high-energy physics, such as confinement. For this collection of topics, one can turn to lattice gauge theories. In Chapter 5, we consider a class of lattice gauge theories particularly suitable for quantum simulation, the Quantum Link Model. The study of the Abelian U(1) Quantum Link Model on a ladder geometry reveals a highly non-trivial phase diagram, featuring a symmetry-protected topological phase. In both Chapters, innovative solutions for the experimental realization of the model Hamiltonians are designed and proposed. To gain numerical access to the ground-state properties and the dynamics of the systems investigated we make use of state-of-the-art numerical methods based on Tensor Networks. The elements of the numerical analysis carried out throughout this thesis are presented in Chapter 6. In the last part we offer an outlook on research perspectives related to the topics discussed in the thesis.
Predicting the quantum dynamics of charged matter interacting with dynamical gauge fields poses an outstanding challenge in theoretical physics. Lacking a generally applicable computational method, quantum simulators offer a promising alternative. \n \nIn this thesis, we contribute to the quantum simulation of high-energy physics, focusing on the platform of ultracold atoms. Using Wilson fermions, we propose to improve implementations of lattice gauge theories based on mixtures of cold atoms in optical lattices. Numerical benchmarks indicate that this makes the realization of Schwinger pair production feasible with current technology. Our proposal is modular and an elementary building is demonstrated experimentally. We further identify dynamical topological transitions, which we discovered in the massive Schwinger model, as a suitable target for quantum simulators. Defining a gauge-invariant order parameter, these transitions are shown to persist beyond weak coupling. \n \nIn the second part of this thesis, we develop a framework for analyzing quantum simulators in terms of experimentally accessible irreducible correlation functions at equal times. We verify this approach numerically for the sine-Gordon model in thermal equilibrium, quantum simulated by two tunnel-coupled superfluids. Finally, we apply our analysis to the non-equilibrium dynamics of a spinor Bose gas, revealing suppressed effective interactions in a strongly-correlated infrared regime.
Can high-energy physics be simulated by low-energy, non-relativistic, many-body systems such as ultracold atoms? Such ultracold atomic systems lack the type of symmetries and dynamical properties of high energy physics models: in particular, they manifest neither local gauge invariance nor Lorentz invariance, which are crucial properties of the quantum field theories which are the building blocks of the standard model of elementary particles. However, it turns out, surprisingly, that there are ways to configure an atomic system to manifest both local gauge invariance and Lorentz invariance. In particular, local gauge invariance can arise either as an effective low-energy symmetry, or as an exact symmetry, following from the conservation laws in atomic interactions. Hence, one could hope that such quantum simulators may lead to a new type of (table-top) experiments which will be used to study various QCD (quantum chromodynamics) phenomena, such as the confinement of dynamical quarks, phase transitions and other effects, which are inaccessible using the currently known computational methods. In this report, we review the Hamiltonian formulation of lattice gauge theories, and then describe our recent progress in constructing the quantum simulation of Abelian and non-Abelian lattice gauge theories in 1 + 1 and 2 + 1 dimensions using ultracold atoms in optical lattices.
Simulating quantum-mechanical systems is a very challenging issue due to the high computational costs. The main reason for this is related to the dimension of the Hilbert space that grows exponentially with the number of degrees of freedom, making an exact implementation impossible even for the most powerful today’s supercomputers. \nOn this ground it appears clear the need of a new simulation method, i.e. quantum simulation. The general idea is very simple: using a controllable quantum system, called quantum simulator, to emulate and to analyze another quantum system that usually results less controllable or accessible. Quantum simulation is currently a growing and multidisciplinary physical area that involves theoretical and experimental research. From the theoretical side, the main focus is on the study of simulation-models that create a mapping between the target physics and the simulator. \nIn the first two parts of this thesis, we explore these ideas by studying the ground-state properties and the real-time dynamics of a class of Zn lattice gauge theories in 1+1 dimensions, in which the gauge fields are coupled to fermionic matter. These models can be considered as quantum simulation-models of lattice 1+1 Quantum Electrodynamics (QED), which is the simplest gauge theory that shows non-trivial phenomena, like confinement, that are also observed in more complicated gauge theories, such as 3+1 Quantum Chromodynamics (QCD). \nIn the third part of this thesis, we unveil an interesting interplay of symmetry and topology by demonstrating the existence of symmetry-protected topological (SPT) phases in the presence of gauge interactions. More specifically, we introduce an alternative discretization of the Schwinger model leading to a new type of topological-QED. We use bosonization and Density-Matrix Renormalization Group (DMRG) techniques to analyze the phase diagram of the model, opening an interesting route to study topological phases of matter in the context of gauge theories.