Showing 1–50 of 150 results for author: Berges, J

Particle production and hadronization temperature of the massive Schwinger model

Authors: Laura Batini, Lara Kuhn, Jürgen Berges, Stefan Floerchinger

Abstract: We study the pair production, string breaking, and hadronization of a receding electron-positron pair using the bosonized version of the massive Schwinger model in quantum electrodynamics in 1+1 space-time dimensions. Specifically, we study the dynamics of the electric field in Bjorken coordinates by splitting it into a coherent field and its Gaussian fluctuations. We find that the electric field… ▽ More We study the pair production, string breaking, and hadronization of a receding electron-positron pair using the bosonized version of the massive Schwinger model in quantum electrodynamics in 1+1 space-time dimensions. Specifically, we study the dynamics of the electric field in Bjorken coordinates by splitting it into a coherent field and its Gaussian fluctuations. We find that the electric field shows damped oscillations, reflecting pair production. Interestingly, the computation of the asymptotic total particle density per rapidity interval for large masses can be fitted using a Boltzmann factor, where the temperature can be related to the hadronization temperature in QCD. Lastly, we discuss the possibility of an analog quantum simulation of the massive Schwinger model using ultracold atoms, explicitly matching the potential of the Schwinger model to the effective potential for the relative phase of two linearly coupled Bose-Einstein condensates. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Coupling renormalization flow in the strongly interacting regime of an asymptotically free quantum field theory in four dimensions

Authors: Jürgen Berges, Razvan Gurau, Hannes Keppler, Thimo Preis

Abstract: We consider a scalar quantum field theory with global $O(N)^3$ symmetry in four Euclidean dimensions and solve it numerically in closed form in the large-N limit. For imaginary tetrahedral coupling the theory is asymptotically free, with stable and real quantum effective action. We demonstrate the dynamical build-up of a strong interaction as the correlation length increases in a regime where the… ▽ More We consider a scalar quantum field theory with global $O(N)^3$ symmetry in four Euclidean dimensions and solve it numerically in closed form in the large-N limit. For imaginary tetrahedral coupling the theory is asymptotically free, with stable and real quantum effective action. We demonstrate the dynamical build-up of a strong interaction as the correlation length increases in a regime where the coupling renormalization flow remains well-defined in the infrared. This is in contrast to perturbative results of asymptotically free theories, which predict that the coupling becomes ill-defined in the infrared, like in quantum chromodynamics. These properties make the model an important laboratory for the study of strong-coupling phenomena in quantum field theory from first principles. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 5+4 pages, 4+2 figures

What ultracold atoms tell us about the real-time dynamics of QCD in extreme conditions

Authors: J. Berges

Abstract: I review developments of how compact table-top setups with ultracold atoms can help us to understand the more complex real-time dynamics of QCD probed in heavy-ion collision experiments. I review developments of how compact table-top setups with ultracold atoms can help us to understand the more complex real-time dynamics of QCD probed in heavy-ion collision experiments. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Proceedings for Quark Matter 2023 plenary talk

Real-time dynamics of false vacuum decay

Authors: Laura Batini, Aleksandr Chatrchyan, Jürgen Berges

Abstract: We investigate false vacuum decay of a relativistic scalar field initialized in the metastable minimum of an asymmetric double-well potential. The transition to the true ground state is a well-defined initial-value problem in real time, which can be formulated in nonequilibrium quantum field theory on a closed time path. We employ the non-perturbative framework of the two-particle irreducible (2PI… ▽ More We investigate false vacuum decay of a relativistic scalar field initialized in the metastable minimum of an asymmetric double-well potential. The transition to the true ground state is a well-defined initial-value problem in real time, which can be formulated in nonequilibrium quantum field theory on a closed time path. We employ the non-perturbative framework of the two-particle irreducible (2PI) quantum effective action at next-to-leading order in a large-N expansion. We also compare to classical-statistical field theory simulations on a lattice in the high-temperature regime. By this, we demonstrate that the real-time decay rates are comparable to those obtained from the conventional Euclidean (bounce) approach. In general, we find that the decay rates are time dependent. For a more comprehensive description of the dynamics, we extract a time-dependent effective potential, which becomes convex during the nonequilibrium transition process. By solving the quantum evolution equations for the one- and two-point correlation functions for vacuum initial conditions, we demonstrate that quantum corrections can lead to transitions that are not captured by classical-statistical approximations. △ Less

Submitted 7 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 18 pages, 12 figures. v2: journal version

Report number: DESY-23-158, NORDITA 2023-082

Journal ref: Phys. Rev. D 109 (2024) 2, 023502

Unconventional charge-density-wave gap in monolayer NbS$_2$

Authors: Timo Knispel, Jan Berges, Arne Schobert, Erik G. C. P. van Loon, Wouter Jolie, Tim Wehling, Thomas Michely, Jeison Fischer

Abstract: Using scanning tunneling microscopy and spectroscopy, for a monolayer of transition metal dichalcogenide H-NbS$_2$ grown by molecular beam epitaxy on graphene, we provide unambiguous evidence for a charge density wave (CDW) with a 3$\times$3 superstructure, which is not present in bulk NbS$_2$. Local spectroscopy displays a pronounced gap of the order of 20 meV at the Fermi level. Within the gap l… ▽ More Using scanning tunneling microscopy and spectroscopy, for a monolayer of transition metal dichalcogenide H-NbS$_2$ grown by molecular beam epitaxy on graphene, we provide unambiguous evidence for a charge density wave (CDW) with a 3$\times$3 superstructure, which is not present in bulk NbS$_2$. Local spectroscopy displays a pronounced gap of the order of 20 meV at the Fermi level. Within the gap low energy features are present. The gap structure with its low energy features is at variance with the expectation for a gap opening in the electronic band structure due to a CDW. Instead, comparison with \it{ab initio} calculations indicates that the observed gap structure must be attributed to combined electron-phonon quasiparticles. The phonons in question are the elusive amplitude and phase collective modes of the CDW transition. Our findings advance the understanding of CDW mechanisms in two dimensional materials and their spectroscopic signatures. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Order parameters for gauge invariant condensation far from equilibrium

Authors: Jürgen Berges, Kirill Boguslavski, Lillian de Bruin, Tara Butler, Jan M. Pawlowski

Abstract: Nuclear collisions at sufficiently high energies are expected to produce far-from-equilibrium matter with a high density of gluons at early times. We show gauge condensation, which occurs as a consequence of the large density of gluons. To identify this condensation phenomenon, we construct two local gauge-invariant observables that carry the macroscopic zero mode of the gauge condensate. The firs… ▽ More Nuclear collisions at sufficiently high energies are expected to produce far-from-equilibrium matter with a high density of gluons at early times. We show gauge condensation, which occurs as a consequence of the large density of gluons. To identify this condensation phenomenon, we construct two local gauge-invariant observables that carry the macroscopic zero mode of the gauge condensate. The first order parameter for gauge condensation investigated here is the correlator of the spatial Polyakov loop. We also consider, for the first time, the correlator of the gauge invariant scalar field, associated to the exponent of the Polyakov loop. Using real-time lattice simulations of classical-statistical $SU(2)$ gauge theory, we find gauge condensation on a system-size dependent time scale $t_{\text{cond}} \sim L^{1/ζ}$ with a universal scaling exponent $ζ$. Furthermore, we suggest an effective theory formulation describing the dynamics using one of the order parameters identified. The formation of a condensate at early times may have intriguing implications for the early stages in heavy ion collisions. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 11 pages, 9 figures, 1 table

Probing universal dynamics with topological data analysis in a gluonic plasma

Authors: Daniel Spitz, Kirill Boguslavski, Jürgen Berges

Abstract: We study nonequilibrium dynamics of SU(2) lattice gauge theory in Minkowski space-time in a classical-statistical regime, where characteristic gluon occupancies are much larger than unity. In this strongly correlated system far from equilibrium, the correlations of energy and topological densities show self-similar behavior related to a turbulent cascade towards higher momentum scales. We employ p… ▽ More We study nonequilibrium dynamics of SU(2) lattice gauge theory in Minkowski space-time in a classical-statistical regime, where characteristic gluon occupancies are much larger than unity. In this strongly correlated system far from equilibrium, the correlations of energy and topological densities show self-similar behavior related to a turbulent cascade towards higher momentum scales. We employ persistent homology to infer topological features of the gluonic plasma via a hierarchy of simplicial and cubical complexes. All topological observables under investigation are also manifestly gauge invariant and are shown to exhibit self-similar evolution, which relate the spatial and temporal properties of the plasma in terms of universal scaling exponents and functions. The findings may help to understand the early stages of heavy-ion collisions in the limit of high energies, and our methods can also facilitate the topological analysis of other complex systems such as encountered in experiments with ultracold quantum gases. △ Less

Submitted 5 December, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 17 pages, 11 figures

Journal ref: Phys. Rev. D 108, 056016 (2023)

Ab initio electron-lattice downfolding: potential energy landscapes, anharmonicity, and molecular dynamics in charge density wave materials

Authors: Arne Schobert, Jan Berges, Erik G. C. P. van Loon, Michael A. Sentef, Sergey Brener, Mariana Rossi, Tim O. Wehling

Abstract: The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the si… ▽ More The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the simulation of electronic properties and nuclear motion - in particular molecular dynamics (MD) simulations. Three different downfolding strategies based on constraining, unscreening, and combinations thereof are benchmarked against full density functional calculations for selected charge density wave (CDW) systems, namely 1H-TaS$_2$, 1T-TiSe$_2$, 1H-NbS$_2$, and a one-dimensional carbon chain. We find that the downfolded models can reproduce potential energy surfaces on supercells accurately and facilitate computational speedup in MD simulations by about five orders of magnitude in comparison to purely ab initio calculations. For monolayer 1H-TaS$_2$ we report classical replica exchange and quantum path integral MD simulations, revealing the impact of thermal and quantum fluctuations on the CDW transition. △ Less

Submitted 16 January, 2024; v1 submitted 13 March, 2023; originally announced March 2023.

Journal ref: SciPost Phys. 16, 046 (2024)

Asymptotic freedom in a strongly interacting scalar quantum field theory in four Euclidean dimensions

Authors: Jürgen Berges, Razvan Gurau, Thimo Preis

Abstract: We show that scalar quantum field theory in four Euclidean dimensions with global $O(N)^3$ symmetry and imaginary tetrahedral coupling is asymptotically free and bounded from below in the large-N limit. While the Hamiltonian is non-Hermitian, the full quantum effective action for the large-N theory only depends on the square of that coupling which is real. A perturbative analysis uncovers that the… ▽ More We show that scalar quantum field theory in four Euclidean dimensions with global $O(N)^3$ symmetry and imaginary tetrahedral coupling is asymptotically free and bounded from below in the large-N limit. While the Hamiltonian is non-Hermitian, the full quantum effective action for the large-N theory only depends on the square of that coupling which is real. A perturbative analysis uncovers that the renormalization group flow of the quartic couplings connects a Gaussian ultraviolet fixed point to a strongly interacting theory in the infrared. This realizes a renormalizable field theory which exhibits non-trivial dynamics, such as direct scattering, while still being analytically tractable also non-perturbatively. Our findings open up a way to address outstanding problems in strongly coupled theories from first principles. △ Less

Submitted 28 July, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 5 pages, 2 figures; published version

Journal ref: Phys. Rev. D 108, 016019 (2023)

Phonon Self-Energy Corrections: To Screen, or Not to Screen

Authors: Jan Berges, Nina Girotto, Tim Wehling, Nicola Marzari, Samuel Poncé

Abstract: First-principles calculations of phonons are often based on the adiabatic approximation and on Brillouin-zone samplings that might not always be sufficient to capture the subtleties of Kohn anomalies. These shortcomings can be addressed through corrections to the phonon self-energy arising from the low-energy electrons. The exact self-energy involves a product of a bare and a screened electron-pho… ▽ More First-principles calculations of phonons are often based on the adiabatic approximation and on Brillouin-zone samplings that might not always be sufficient to capture the subtleties of Kohn anomalies. These shortcomings can be addressed through corrections to the phonon self-energy arising from the low-energy electrons. The exact self-energy involves a product of a bare and a screened electron-phonon vertex [Rev. Mod. Phys. 89, 015003 (2017)]; still, calculations often employ two adiabatically screened vertices, which have been proposed as a reliable approximation for self-energy differences [Phys. Rev. B 82, 165111 (2010)]. We assess the accuracy of both approaches in estimating the phonon spectral functions of model Hamiltonians and the adiabatic low-temperature phonon dispersions of monolayer TaS$_2$ and doped MoS$_2$. We find that the approximate method yields excellent corrections at low computational cost, due to its designed error cancellation to first order, while using a bare vertex could in principle improve these results but is challenging in practice. We offer an alternative strategy based on downfolding to partially screened phonons and interactions [Phys. Rev. B 92, 245108 (2015)]. This is a natural scheme to include electron-electron interactions and tackle phonons in strongly correlated materials and the frequency dependence of the electron-phonon vertex. △ Less

Submitted 18 October, 2023; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: Published version, 26 pages, 12 figures

Journal ref: Phys. Rev. X 13, 041009 (2023)

Stable and unstable perturbations in universal scaling phenomena far from equilibrium

Authors: Thimo Preis, Michal P. Heller, Jürgen Berges

Abstract: We study the dynamics of perturbations around nonthermal fixed points associated to universal scaling phenomena in quantum many-body systems far from equilibrium. For an N-component scalar quantum field theory in 3+1 space-time dimensions, we determine the stability scaling exponents using a self-consistent large-N expansion to next-to-leading order. Our analysis reveals the presence of both stabl… ▽ More We study the dynamics of perturbations around nonthermal fixed points associated to universal scaling phenomena in quantum many-body systems far from equilibrium. For an N-component scalar quantum field theory in 3+1 space-time dimensions, we determine the stability scaling exponents using a self-consistent large-N expansion to next-to-leading order. Our analysis reveals the presence of both stable and unstable perturbations, the latter leading to quasi-exponential deviations from the fixed point in the infrared. We identify a tower of far-from-equilibrium quasi-particle states and their dispersion relations by computing the spectral function. With the help of linear response theory, we demonstrate that unstable dynamics arises from a competition between elastic scattering processes among the quasi-particle states. What ultimately renders the fixed point dynamically attractive is the phenomenon of a scaling instability, which is the universal scaling of the unstable regime towards the infrared due to a self-similar quasi-particle cascade. Our results provide ab initio understanding of emergent stability properties in self-organized scaling phenomena. △ Less

Submitted 18 December, 2022; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 6 pages, 5 figures, PRL version (minor text changes)

Journal ref: Phys. Rev. Lett. 130, 031602 (2023)

Condensation and thermalization of an easy-plane ferromagnet in a spinor Bose gas

Authors: Maximilian Prüfer, Daniel Spitz, Stefan Lannig, Helmut Strobel, Jürgen Berges, Markus K. Oberthaler

Abstract: The extensive control of spin makes spintronics a promising candidate for future scalable quantum devices. For the generation of spin-superfluid systems, a detailed understanding of the build-up of coherence and relaxation is necessary. However, to determine the relevant parameters for robust coherence properties and faithfully witnessing thermalization, the direct access to space- and time-resolv… ▽ More The extensive control of spin makes spintronics a promising candidate for future scalable quantum devices. For the generation of spin-superfluid systems, a detailed understanding of the build-up of coherence and relaxation is necessary. However, to determine the relevant parameters for robust coherence properties and faithfully witnessing thermalization, the direct access to space- and time-resolved spin observables is needed. Here, we study the thermalization of an easy-plane ferromagnet employing a homogeneous one-dimensional spinor Bose gas. Building on the pristine control of preparation and readout we demonstrate the dynamic emergence of long-range coherence for the spin field and verify spin-superfluidity by experimentally testing Landau's criterion. We reveal the structure of the emergent quasi-particles: one 'massive'(Higgs) mode, and two 'massless' (Goldstone) modes - a consequence of explicit and spontaneous symmetry breaking, respectively. Our experiments allow for the first time to observe the thermalization of an easy-plane ferromagnetic Bose gas; we find agreement for the relevant momentum-resolved observables with a thermal prediction obtained from an underlying microscopic model within the Bogoliubov approximation. Our methods and results pave the way towards a quantitative understanding of condensation dynamics in large magnetic spin systems and the study of the role of entanglement and topological excitations for its thermalization. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Equal-time approach to real-time dynamics of quantum fields

Authors: Robert Ott, Torsten V. Zache, Jürgen Berges

Abstract: We employ the equal-time formulation of quantum field theory to derive effective kinetic theories, first for a weakly coupled non-relativistic Bose gas, and then for a strongly correlated system of self-interacting N-component fields. Our results provide the link between state-of-the-art measurements of equal-time effective actions using quantum simulator platforms, as employed in Refs. [1, 2], an… ▽ More We employ the equal-time formulation of quantum field theory to derive effective kinetic theories, first for a weakly coupled non-relativistic Bose gas, and then for a strongly correlated system of self-interacting N-component fields. Our results provide the link between state-of-the-art measurements of equal-time effective actions using quantum simulator platforms, as employed in Refs. [1, 2], and observables underlying effective kinetic or hydrodynamic descriptions. New non-perturbative approximation schemes can be developed and certified this way, where the a priori time-local formulation of the equal-time effective action has crucial advantages over the conventional closed-time-path approach which is non-local in time. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 18 pages, 3 figures

Non-equilibrium dynamics of fluctuations in an ultra-cold atomic mixture

Authors: Apoorva Hegde, Robert Ott, Andy Xia, Valentin Kasper, Jürgen Berges, Fred Jendrzejewski

Abstract: We investigate an ultra-cold mixture of Bose gases interacting via spin-changing collisions by studying the dynamics of spin fluctuations. The experimental implementation employs $^{23}$Na and $^{7}$Li atoms, which are prepared out of equilibrium across a wide range of initial conditions. We identify three regimes in the dynamics of the system for different initial states: a long-lived metastable… ▽ More We investigate an ultra-cold mixture of Bose gases interacting via spin-changing collisions by studying the dynamics of spin fluctuations. The experimental implementation employs $^{23}$Na and $^{7}$Li atoms, which are prepared out of equilibrium across a wide range of initial conditions. We identify three regimes in the dynamics of the system for different initial states: a long-lived metastable regime, an instability range with strong growth of fluctuations, and a fast relaxing regime approaching thermal equilibrium. Theoretical modelling of the data allows us to reconstruct effective potentials which characterize the different dynamical regimes of the system. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 9 pages, 5 figures

Stability analysis of non-thermal fixed points in longitudinally expanding kinetic theory

Authors: Aleksandr N. Mikheev, Aleksas Mazeliauskas, Jürgen Berges

Abstract: We use the Hamiltonian formulation of kinetic theory to perform a stability analysis of non-thermal fixed points in a non-Abelian plasma. We construct a perturbative expansion of the Fokker-Planck collision kernel in an adiabatic approximation and show that the (next-to-)leading order solutions reproduce the known non-thermal fixed point scaling exponents. Working at next-to-leading order, we deri… ▽ More We use the Hamiltonian formulation of kinetic theory to perform a stability analysis of non-thermal fixed points in a non-Abelian plasma. We construct a perturbative expansion of the Fokker-Planck collision kernel in an adiabatic approximation and show that the (next-to-)leading order solutions reproduce the known non-thermal fixed point scaling exponents. Working at next-to-leading order, we derive the stability equations for scaling exponents and find the relaxation rate to the non-thermal fixed point. This approach provides the basis for an understanding of the prescaling phenomena observed in QCD kinetic theory and non-relativistic Bose gas systems. △ Less

Submitted 29 June, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 9 pages, 2 figures. Minor changes, published version

Report number: CERN-TH-2022-030

Journal ref: Phys. Rev. D 105, 116025 (2022)

Engineering Interlayer Hybridization in Energy Space via Dipolar Overlayers

Authors: Bin Shao, Xiao Jiang, Jan Berges, Sheng Meng, Bing Huang

Abstract: The interlayer hybridization (IH) of van der Waals (vdW) materials is thought to be mostly associated with the unignorable interlayer overlaps of wavefunctions ($t$) in real space. Here, we develop a more fundamental understanding of IH by introducing a new physical quantity, the IH admixture ratio $α$. Consequently, an exotic strategy of IH engineering in energy space can be proposed, i.e., inste… ▽ More The interlayer hybridization (IH) of van der Waals (vdW) materials is thought to be mostly associated with the unignorable interlayer overlaps of wavefunctions ($t$) in real space. Here, we develop a more fundamental understanding of IH by introducing a new physical quantity, the IH admixture ratio $α$. Consequently, an exotic strategy of IH engineering in energy space can be proposed, i.e., instead of changing t as commonly used, $α$ can be effectively tuned in energy space by changing the onsite energy difference ($2Δ$) between neighboring-layer states. In practice, this is feasible via reshaping the electrostatic potential of the surface by deposing a dipolar overlayer, e.g., crystalline ice. Our first-principles calculations unveil that IH engineering via adjusting $2Δ$ can greatly tune interlayer optical transitions in transition-metal dichalcogenide bilayers, switch different types of Dirac surface states in Bi$_2$Se$_3$ thin films, and control magnetic phase transition of charge density waves in 1H/1T-TaS$_2$ bilayers, opening new opportunities to govern the fundamental optoelectronic, topological, and magnetic properties of vdW systems beyond the traditional interlayer-distance or twisting engineering. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 12 pages, 4 figures

Journal ref: Chin. Phys. Lett. 40, 087303 (2023)

Doping fingerprints of spin and lattice fluctuations in moiré superlattice systems

Authors: Niklas Witt, José M. Pizarro, Jan Berges, Takuya Nomoto, Ryotaro Arita, Tim O. Wehling

Abstract: Twisted Van der Waals systems offer the unprecedented possibility to tune different states of correlated quantum matter with an external non-invasive electrostatic doping. The nature of the superconducting order presents a recurring open question in this context. In this work, we quantitatively assess the case of spin-fluctuation-mediated pairing for $Γ$-valley twisted transition metal dichalcogen… ▽ More Twisted Van der Waals systems offer the unprecedented possibility to tune different states of correlated quantum matter with an external non-invasive electrostatic doping. The nature of the superconducting order presents a recurring open question in this context. In this work, we quantitatively assess the case of spin-fluctuation-mediated pairing for $Γ$-valley twisted transition metal dichalcogenide homobilayers. We self-consistently and dynamically calculate the doping dependent superconducting transition temperature $T_{\mathrm{c}}$ revealing a superconducting dome with a maximal $T_{\mathrm{c}}\approx 0.1-1$ K depending on twist angle. We compare our results with conventional phonon-mediated superconductivity and identify clear fingerprints in the doping dependence of $T_{\mathrm{c}}$, which allow experiments to distinguish between different pairing mechanisms. △ Less

Submitted 18 March, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 9 pages, 3 figures + SM with 18 pages, 11 figures

Thermalization dynamics of a gauge theory on a quantum simulator

Authors: Zhao-Yu Zhou, Guo-Xian Su, Jad C. Halimeh, Robert Ott, Hui Sun, Philipp Hauke, Bing Yang, Zhen-Sheng Yuan, Jürgen Berges, Jian-Wei Pan

Abstract: Gauge theories form the foundation of modern physics, with applications ranging from elementary particle physics and early-universe cosmology to condensed matter systems. We perform quantum simulations of the unitary dynamics of a U(1) symmetric gauge field theory and demonstrate emergent irreversible behavior. The highly constrained gauge theory dynamics is encoded in a one-dimensional Bose--Hubb… ▽ More Gauge theories form the foundation of modern physics, with applications ranging from elementary particle physics and early-universe cosmology to condensed matter systems. We perform quantum simulations of the unitary dynamics of a U(1) symmetric gauge field theory and demonstrate emergent irreversible behavior. The highly constrained gauge theory dynamics is encoded in a one-dimensional Bose--Hubbard simulator, which couples fermionic matter fields through dynamical gauge fields. We investigate global quantum quenches and the equilibration to a steady state well approximated by a thermal ensemble. Our work may enable the investigation of elusive phenomena, such as Schwinger pair production and string-breaking, and paves the way for simulating more complex higher-dimensional gauge theories on quantum synthetic matter devices. △ Less

Submitted 15 July, 2022; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: Accepted version, 6+6 pages, 4+5 figures

Journal ref: Science 377, 311-314 (2022)

Dimensional crossover for universal scaling far from equilibrium

Authors: Lasse Gresista, Torsten V. Zache, Jürgen Berges

Abstract: We perform a dynamical finite-size scaling analysis of a nonequilibrium Bose gas which is confined in the transverse plane. Varying the transverse size, we establish a dimensional crossover for universal scaling properties far from equilibrium. Our results suggest that some aspects of the dynamical universal behavior of anisotropic systems can be classified in terms of fractional spatial dimension… ▽ More We perform a dynamical finite-size scaling analysis of a nonequilibrium Bose gas which is confined in the transverse plane. Varying the transverse size, we establish a dimensional crossover for universal scaling properties far from equilibrium. Our results suggest that some aspects of the dynamical universal behavior of anisotropic systems can be classified in terms of fractional spatial dimensions. We discuss our findings in view of recent experimental results with quasi one-dimensional setups of quenched ultracold quantum gases. △ Less

Submitted 25 July, 2021; originally announced July 2021.

Comments: 9 pages, 7 figures

Downfolding the Su-Schrieffer-Heeger model

Authors: Arne Schobert, Jan Berges, Tim Wehling, Erik van Loon

Abstract: Charge-density waves are responsible for symmetry-breaking displacements of atoms and concomitant changes in the electronic structure. Linear response theories, in particular density-functional perturbation theory, provide a way to study the effect of displacements on both the total energy and the electronic structure based on a single ab initio calculation. In downfolding approaches, the electron… ▽ More Charge-density waves are responsible for symmetry-breaking displacements of atoms and concomitant changes in the electronic structure. Linear response theories, in particular density-functional perturbation theory, provide a way to study the effect of displacements on both the total energy and the electronic structure based on a single ab initio calculation. In downfolding approaches, the electronic system is reduced to a smaller number of bands, allowing for the incorporation of additional correlation and environmental effects on these bands. However, the physical contents of this downfolded model and its potential limitations are not always obvious. Here, we study the potential-energy landscape and electronic structure of the Su-Schrieffer-Heeger (SSH) model, where all relevant quantities can be evaluated analytically. We compare the exact results at arbitrary displacement with diagrammatic perturbation theory both in the full model and in a downfolded effective single-band model, which gives an instructive insight into the properties of downfolding. An exact reconstruction of the potential-energy landscape is possible in a downfolded model, which requires a dynamical electron-biphonon interaction. The dispersion of the bands upon atomic displacement is also found correctly, where the downfolded model by construction only captures spectral weight in the target space. In the SSH model, the electron-phonon coupling mechanism involves exclusively hybridization between the low- and high-energy bands and this limits the computational efficiency gain of downfolded models. △ Less

Submitted 17 August, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: Submission to SciPost. Python file to generate the plots included in the arxiv source

Journal ref: SciPost Phys. 11, 079 (2021)

Collisional strong-field QED kinetic equations from first principles

Authors: Gregor Fauth, Jürgen Berges, Antonino Di Piazza

Abstract: Starting from nonequilibrium quantum field theory on a closed time path, we derive kinetic equations for the strong-field regime of quantum electrodynamics (QED) using a systematic expansion in the gauge coupling $e$. The strong field regime is characterized by a large photon field of order $\mathcal{O}(1/e)$, which is relevant for the description of, e.g., intense laser fields, the initial stages… ▽ More Starting from nonequilibrium quantum field theory on a closed time path, we derive kinetic equations for the strong-field regime of quantum electrodynamics (QED) using a systematic expansion in the gauge coupling $e$. The strong field regime is characterized by a large photon field of order $\mathcal{O}(1/e)$, which is relevant for the description of, e.g., intense laser fields, the initial stages of off-central heavy ion collisions, and condensed matter systems with net fermion number. The strong field enters the dynamical equations via both quantum Vlasov and collision terms, which we derive to order $\mathcal{O}(e^2)$. The kinetic equations feature generalized scattering amplitudes that have their own equation of motion in terms of the fermion spectral function. The description includes single photon emission, electron-positron pair photoproduction, vacuum (Schwinger) pair production, their inverse processes, medium effects and contributions from the field, which are not restricted to the so-called locally-constant crossed field approximation. This extends known kinetic equations commonly used in strong-field QED of intense laser fields. In particular, we derive an expression for the asymptotic fermion pair number that includes leading-order collisions and remains valid for strongly inhomogeneous fields. For the purpose of analytically highlighting limiting cases, we also consider plane-wave fields for which it is shown how to recover Furry-picture scattering amplitudes by further assuming negligible occupations. Known on-shell descriptions are recovered in the case of simply peaked ultrarelativistic fermion occupations. Collisional strong-field equations are necessary to describe the dynamics to thermal equilibrium starting from strong-field initial conditions. △ Less

Submitted 10 August, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 45 pages, 8 figures; journal version (previous version only implemented changes of first proof), minor changes

Journal ref: Phys. Rev. D 104, 036007 (2021)

Downfolding approaches to electron-ion coupling: Constrained density-functional perturbation theory for molecules

Authors: Erik G. C. P. van Loon, Jan Berges, Tim O. Wehling

Abstract: Constrained electronic-structure theories enable the construction of effective low-energy models consisting of partially dressed particles. However, the interpretation and physical content of these theories is not straightforward. Here, we carefully explore the properties of downfolding theories for electron-ion problems, in particular constrained density-functional perturbation theory (cDFPT). We… ▽ More Constrained electronic-structure theories enable the construction of effective low-energy models consisting of partially dressed particles. However, the interpretation and physical content of these theories is not straightforward. Here, we carefully explore the properties of downfolding theories for electron-ion problems, in particular constrained density-functional perturbation theory (cDFPT). We show that the dipole selection rules determine whether the partially dressed phonons satisfy Goldstone's theorem, and we prove that electronic screening always lowers the phonon frequencies. We illustrate the theory with cDFPT calculations for minimal example systems: the nitrogen and benzene molecule as well as graphene. △ Less

Submitted 10 May, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Comments: Revised manuscript with new appendix and ancillary scripts and data

Journal ref: Phys. Rev. B 103, 205103 (2021)

A full gap above the Fermi level: the charge density wave of monolayer VS2

Authors: Camiel van Efferen, Jan Berges, Joshua Hall, Erik van Loon, Stefan Kraus, Arne Schobert, Tobias Wekking, Felix Huttmann, Eline Plaar, Nico Rothenbach, Katharina Ollefs, Lucas Machado Arruda, Nick Brookes, Gunnar Schoenhoff, Kurt Kummer, Heiko Wende, Tim Wehling, Thomas Michely

Abstract: In the standard model of charge density wave (CDW) transitions, the displacement along a single phonon mode lowers the total electronic energy by creating a gap at the Fermi level, making the CDW a metal--insulator transition. Here, using scanning tunneling microscopy and spectroscopy and ab initio calculations, we show that VS$_2$ realizes a CDW which stands out of this standard model. There is a… ▽ More In the standard model of charge density wave (CDW) transitions, the displacement along a single phonon mode lowers the total electronic energy by creating a gap at the Fermi level, making the CDW a metal--insulator transition. Here, using scanning tunneling microscopy and spectroscopy and ab initio calculations, we show that VS$_2$ realizes a CDW which stands out of this standard model. There is a full CDW gap residing in the unoccupied states of monolayer VS$_2$. At the Fermi level, the CDW induces a topological metal-metal (Lifshitz) transition. Non-linear coupling of transverse and longitudinal phonons is essential for the formation of the CDW and the full gap above the Fermi level. Additionally, x-ray magnetic circular dichroism reveals the absence of net magnetization in this phase, pointing to coexisting charge and spin density waves in the ground state. △ Less

Submitted 8 September, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Scalable cold-atom quantum simulator for two-dimensional QED

Authors: Robert Ott, Torsten V. Zache, Fred Jendrzejewski, Jürgen Berges

Abstract: We propose a scalable analog quantum simulator for quantum electrodynamics (QED) in two spatial dimensions. The setup for the U(1) lattice gauge field theory employs inter-species spin-changing collisions in an ultra-cold atomic mixture trapped in an optical lattice. Building on the previous one-dimensional implementation scheme of Ref. [1] we engineer spatial plaquette terms for magnetic fields,… ▽ More We propose a scalable analog quantum simulator for quantum electrodynamics (QED) in two spatial dimensions. The setup for the U(1) lattice gauge field theory employs inter-species spin-changing collisions in an ultra-cold atomic mixture trapped in an optical lattice. Building on the previous one-dimensional implementation scheme of Ref. [1] we engineer spatial plaquette terms for magnetic fields, thus solving a major obstacle towards experimental realizations of realistic gauge theories in higher dimensions. We apply our approach to pure gauge theory and discuss how the phenomenon of confinement of electric charges can be described by the quantum simulator. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 127, 130504 (2021)

Analog cosmological reheating in an ultracold Bose gas

Authors: Aleksandr Chatrchyan, Kevin T. Geier, Markus K. Oberthaler, Jürgen Berges, Philipp Hauke

Abstract: Cosmological reheating describes the transition of the post-inflationary universe to a hot and thermal state. In order to shed light on the underlying dynamics of this process, we propose to quantum-simulate the reheating-like dynamics of a generic cosmological single-field model in an ultracold Bose gas. In our setup, the excitations on top of an atomic Bose-Einstein condensate play the role of t… ▽ More Cosmological reheating describes the transition of the post-inflationary universe to a hot and thermal state. In order to shed light on the underlying dynamics of this process, we propose to quantum-simulate the reheating-like dynamics of a generic cosmological single-field model in an ultracold Bose gas. In our setup, the excitations on top of an atomic Bose-Einstein condensate play the role of the particles produced by the decaying inflaton field after inflation. Expanding spacetime as well as the background oscillating inflaton field are mimicked in the non-relativistic limit by a time dependence of the atomic interactions, which can be tuned experimentally via Feshbach resonances. As we illustrate by means of classical-statistical simulations for the case of two spatial dimensions, the dynamics of the atomic system exhibits the characteristic stages of far-from-equilibrium reheating, including the amplification of fluctuations via parametric instabilities and the subsequent turbulent transport of energy towards higher momenta. The transport is governed by a non-thermal fixed point showing universal self-similar time evolution as well as a transient regime of prescaling with time-dependent scaling exponents. While the classical-statistical simulations can capture only the earlier stages of the dynamics for weak couplings, the proposed experiment has the potential of exploring the evolution up to late times even beyond the weak coupling regime. △ Less

Submitted 6 August, 2021; v1 submitted 5 August, 2020; originally announced August 2020.

Comments: 18+12 pages, 6+4 figures. Version published in PRA. Changes in this version: improved presentation of the analogy to cosmological reheating, moved some technical paragraphs of Secs. II and III to the appendix. A video of the preheating dynamics can be downloaded from the section "ancillary files" on the article's arXiv page

Journal ref: Phys. Rev. A 104, 023302 (2021)

QCD thermalization: Ab initio approaches and interdisciplinary connections

Authors: Jürgen Berges, Michal P. Heller, Aleksas Mazeliauskas, Raju Venugopalan

Abstract: Heavy-ion collisions at BNL's Relativistic Heavy Ion Collider and CERN's Large Hadron Collider provide strong evidence for the formation of a quark-gluon plasma, with temperatures extracted from relativistic viscous hydrodynamic simulations shown to be well above the transition temperature from hadron matter. How the strongly correlated quark-gluon matter forms in a heavy-ion collision, its proper… ▽ More Heavy-ion collisions at BNL's Relativistic Heavy Ion Collider and CERN's Large Hadron Collider provide strong evidence for the formation of a quark-gluon plasma, with temperatures extracted from relativistic viscous hydrodynamic simulations shown to be well above the transition temperature from hadron matter. How the strongly correlated quark-gluon matter forms in a heavy-ion collision, its properties off-equilibrium, and the thermalization process in the plasma, are outstanding problems in QCD. We review here the theoretical progress in this field in weak coupling QCD effective field theories and in strong coupling holographic approaches based on gauge-gravity duality. We outline the interdisciplinary connections of different stages of the thermalization process to non-equilibrium dynamics in other systems across energy scales ranging from inflationary cosmology, to strong field QED, to ultracold atomic gases, with emphasis on the universal dynamics of non-thermal and of hydrodynamic attractors. We survey measurements in heavy-ion collisions that are sensitive to the early non-equilibrium stages of the collision and discuss the potential for future measurements. We summarize the current state-of-the art in thermalization studies and identify promising avenues for further progress. △ Less

Submitted 17 August, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 73 pages, 32 figures, prepared for Reviews of Modern Physics; version 2: small improvements and additions, submitted version; version 3: accepted version; title changed, differs in style and format from published version in Reviews of Modern Physics

Report number: CERN-TH-2020-080

Journal ref: Rev. Mod. Phys. 93, 35003 (2021)

Thermalization and dynamical spectral properties in the quark-meson model

Authors: Linda Shen, Jürgen Berges, Jan M. Pawlowski, Alexander Rothkopf

Abstract: We investigate the nonequilibrium evolution of the quark-meson model using two-particle irreducible effective action techniques. Our numerical simulations, which include the full dynamics of the order parameter of chiral symmetry, show how the model thermalizes into different regions of its phase diagram. In particular, by studying quark and meson spectral functions, we shed light on the real-time… ▽ More We investigate the nonequilibrium evolution of the quark-meson model using two-particle irreducible effective action techniques. Our numerical simulations, which include the full dynamics of the order parameter of chiral symmetry, show how the model thermalizes into different regions of its phase diagram. In particular, by studying quark and meson spectral functions, we shed light on the real-time dynamics approaching the crossover transition, revealing e.g. the emergence of light effective fermionic degrees of freedom in the infrared. At late times in the evolution, the fluctuation-dissipation relation emerges naturally among both meson and quark degrees of freedom, confirming that the simulation successfully reaches thermal equilibrium. △ Less

Submitted 2 August, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Journal ref: Phys. Rev. D 102, 016012 (2020)

Finding self-similar behavior in quantum many-body dynamics via persistent homology

Authors: Daniel Spitz, Jürgen Berges, Markus K. Oberthaler, Anna Wienhard

Abstract: Inspired by topological data analysis techniques, we introduce persistent homology observables and apply them in a geometric analysis of the dynamics of quantum field theories. As a prototype application, we consider data from a classical-statistical simulation of a two-dimensional Bose gas far from equilibrium. We discover a continuous spectrum of dynamical scaling exponents, which provides a ref… ▽ More Inspired by topological data analysis techniques, we introduce persistent homology observables and apply them in a geometric analysis of the dynamics of quantum field theories. As a prototype application, we consider data from a classical-statistical simulation of a two-dimensional Bose gas far from equilibrium. We discover a continuous spectrum of dynamical scaling exponents, which provides a refined classification of nonequilibrium self-similar phenomena. A possible explanation of the underlying processes is provided in terms of mixing strong wave turbulence and anomalous vortex kinetics components in point clouds. We find that the persistent homology scaling exponents are inherently linked to the geometry of the system, as the derivation of a packing relation reveals. The approach opens new ways of analyzing quantum many-body dynamics in terms of robust topological structures beyond standard field theoretic techniques. △ Less

Submitted 22 September, 2021; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 21 pages, 18 figures

Journal ref: SciPost Phys. 11, 060 (2021)

Spectral, statistical and vertex functions in scalar quantum field theory far from equilibrium

Authors: Linda Shen, Jürgen Berges

Abstract: We compute the far-from-equilibrium dynamics of relativistic scalar quantum fields in 3+1 space-time dimensions starting from over-occupied initial conditions. We determine universal scaling exponents and functions for two-point correlators and the four-vertex in a self-similar regime in time and space or momenta. The scaling form of the momentum-dependent four-vertex exhibits a dramatic fall-off… ▽ More We compute the far-from-equilibrium dynamics of relativistic scalar quantum fields in 3+1 space-time dimensions starting from over-occupied initial conditions. We determine universal scaling exponents and functions for two-point correlators and the four-vertex in a self-similar regime in time and space or momenta. The scaling form of the momentum-dependent four-vertex exhibits a dramatic fall-off towards low momenta. Comparing spectral functions (commutators) and statistical correlations (anti-commutators) of field operators allows us to detect strong violations of the fluctuation-dissipation relation in this non-perturbative infrared regime. Based on a self-consistent expansion in the number of field components to next-to-leading order, a wide range of interaction strengths is analyzed and compared to weak-coupling estimates in effective kinetic theory and classical-statistical field theory. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Comments: 15 pages, 12 figures

Journal ref: Phys. Rev. D 101, 056009 (2020)

Ab initio phonon self-energies and fluctuation diagnostics of phonon anomalies: Lattice instabilities from Dirac pseudospin physics in transition metal dichalcogenides

Authors: Jan Berges, Erik G. C. P. van Loon, Arne Schobert, Malte Rösner, Tim O. Wehling

Abstract: We present an ab initio approach for the calculation of phonon self-energies and their fluctuation diagnostics, which allows us to identify the electronic processes behind phonon anomalies. Application to the transition-metal-dichalcogenide monolayer 1H-TaS$_2$ reveals that coupling between the longitudinal-acoustic phonons and the electrons from an isolated low-energy metallic band is entirely re… ▽ More We present an ab initio approach for the calculation of phonon self-energies and their fluctuation diagnostics, which allows us to identify the electronic processes behind phonon anomalies. Application to the transition-metal-dichalcogenide monolayer 1H-TaS$_2$ reveals that coupling between the longitudinal-acoustic phonons and the electrons from an isolated low-energy metallic band is entirely responsible for phonon anomalies such as the mode softening and associated charge-density waves observed in this material. Our analysis allows us to distinguish between different mode-softening mechanisms including matrix-element effects, Fermi-surface nesting, and Van Hove scenarios. We find that matrix-element effects originating from a peculiar type of Dirac pseudospin textures control the charge-density-wave physics in 1H-TaS$_2$ and similar transition metal dichalcogenides. △ Less

Submitted 3 November, 2021; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: Revised manuscript with ancillary scripts and data, 15 pages, 9 figures

Journal ref: Phys. Rev. B 101, 155107 (2020)

Extracting the field theory description of a quantum many-body system from experimental data

Authors: Torsten V. Zache, Thomas Schweigler, Sebastian Erne, Jörg Schmiedmayer, Jürgen Berges

Abstract: Quantum field theory is a powerful tool to describe the relevant physics governing complex quantum many-body systems. Here we develop a general pathway to extract the irreducible building blocks of quantum field theoretical descriptions and its parameters purely from experimental data. This is accomplished by extracting the one-particle irreducible (1PI) vertices from which one can construct all o… ▽ More Quantum field theory is a powerful tool to describe the relevant physics governing complex quantum many-body systems. Here we develop a general pathway to extract the irreducible building blocks of quantum field theoretical descriptions and its parameters purely from experimental data. This is accomplished by extracting the one-particle irreducible (1PI) vertices from which one can construct all observables. To match the capabilities of experimental techniques used in quantum simulation experiments, our approach employs a formulation of quantum field theory based on equal-time correlation functions only. We illustrate our procedure by applying it to the quantum sine-Gordon model in thermal equilibrium. The theoretical foundations are illustrated by estimating the irreducible vertices at equal times both analytically and using numerical simulations. We then demonstrate explicitly how to extract these quantities from an experiment where we quantum simulate the sine-Gordon model by two tunnel-coupled superfluids. We extract the full two-point function and the interaction vertex (four-point function) and their variation with momentum, encoding the `running' of the couplings. The measured 1PI vertices are compared to the theoretical estimates, verifying our procedure. Our work opens new ways of addressing fundamental questions in quantum field theory, which are relevant in high-energy and condensed matter physics, and in taking quantum phenomena from fundamental science to practical technology. △ Less

Submitted 11 April, 2020; v1 submitted 27 September, 2019; originally announced September 2019.

Comments: 18 pages, 11 figures. Updated Phys. Rev. X version with minor changes

Journal ref: Phys. Rev. X 10, 011020 (2020)

Probing the evolution of heavy-ion collisions using direct photon interferometry

Authors: Oscar Garcia-Montero, Nicole Löher, Aleksas Mazeliauskas, Jürgen Berges, Klaus Reygers

Abstract: We investigate the measurement of Hanbury Brown-Twiss (HBT) photon correlations as an experimental tool to discriminate different sources of photon enhancement, which are proposed to simultaneously reproduce the direct photon yield and the azimuthal anisotropy measured in nuclear collisions at RHIC and the LHC. To showcase this, we consider two different scenarios in which we enhance the yields fr… ▽ More We investigate the measurement of Hanbury Brown-Twiss (HBT) photon correlations as an experimental tool to discriminate different sources of photon enhancement, which are proposed to simultaneously reproduce the direct photon yield and the azimuthal anisotropy measured in nuclear collisions at RHIC and the LHC. To showcase this, we consider two different scenarios in which we enhance the yields from standard hydrodynamical simulations. In the first, additional photons are produced from the early pre-equilibrium stage computed from the \textit{bottom-up} thermalization scenario. In the second, the thermal rates are enhanced close to the pseudo-critical temperature $T_c\approx 155\,\text{MeV}$ using a phenomenological ansatz. We compute the correlators for relative momenta $q_o, \,q_s$ and $q_l$ for different transverse pair momenta, $K_\perp$, and find that the longitudinal correlation is the most sensitive to different photon sources. Our results also demonstrate that including anisotropic pre-equilibrium rates enhances non-Gaussianities in the correlators, which can be quantified using the kurtosis of the correlators. Finally, we study the feasibility of measuring a direct photon HBT signal in the upcoming high-luminosity LHC runs. Considering only statistical uncertainties, we find that with the projected $\sim 10^{10}$ heavy ion events a measurement of the HBT correlations for $K_\perp<1\, \text{GeV}$ is statistically significant. △ Less

Submitted 30 October, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. C 102, 024915 (2020)

Realizing a scalable building block of a U(1) gauge theory with cold atomic mixtures

Authors: Alexander Mil, Torsten V. Zache, Apoorva Hegde, Andy Xia, Rohit P. Bhatt, Markus K. Oberthaler, Philipp Hauke, Jürgen Berges, Fred Jendrzejewski

Abstract: In the fundamental laws of physics, gauge fields mediate the interaction between charged particles. An example is quantum electrodynamics -- the theory of electrons interacting with the electromagnetic field -- based on U(1) gauge symmetry. Solving such gauge theories is in general a hard problem for classical computational techniques. While quantum computers suggest a way forward, it is difficult… ▽ More In the fundamental laws of physics, gauge fields mediate the interaction between charged particles. An example is quantum electrodynamics -- the theory of electrons interacting with the electromagnetic field -- based on U(1) gauge symmetry. Solving such gauge theories is in general a hard problem for classical computational techniques. While quantum computers suggest a way forward, it is difficult to build large-scale digital quantum devices required for complex simulations. Here, we propose a fully scalable analog quantum simulator of a U(1) gauge theory in one spatial dimension. To engineer the local gauge symmetry, we employ inter-species spin-changing collisions in an atomic mixture. We demonstrate the experimental realization of the elementary building block as a key step towards a platform for large-scale quantum simulations of continuous gauge theories. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Gauge-invariant condensation in the nonequilibrium quark-gluon plasma

Authors: Jürgen Berges, Kirill Boguslavski, Mark Mace, Jan M. Pawlowski

Abstract: The large density of gluons, which is present shortly after a nuclear collision at very high energies, can lead to the formation of a condensate. We identify a gauge-invariant order parameter for condensation based on elementary non-perturbative excitations of the plasma, which are described by spatial Wilson loops. Using real-time lattice simulations, we demonstrate that a self-similar transport… ▽ More The large density of gluons, which is present shortly after a nuclear collision at very high energies, can lead to the formation of a condensate. We identify a gauge-invariant order parameter for condensation based on elementary non-perturbative excitations of the plasma, which are described by spatial Wilson loops. Using real-time lattice simulations, we demonstrate that a self-similar transport process towards low momenta builds up a macroscopic zero mode. Our findings reveal intriguing similarities to recent discoveries of condensation phenomena out of equilibrium in table-top experiments with ultracold Bose gases. △ Less

Submitted 27 July, 2020; v1 submitted 13 September, 2019; originally announced September 2019.

Comments: 5 pages, 3 figures; v2: journal version, add. section comparing to thermal case and extended discussion

Journal ref: Phys. Rev. D 102, 034014 (2020)

Experimental extraction of the quantum effective action for a non-equilibrium many-body system

Authors: Maximilian Prüfer, Torsten V. Zache, Philipp Kunkel, Stefan Lannig, Alexis Bonnin, Helmut Strobel, Jürgen Berges, Markus K. Oberthaler

Abstract: Far-from-equilibrium situations are ubiquitous in nature. They are responsible for a wealth of phenomena, which are not simple extensions of near-equilibrium properties, ranging from fluid flows turning turbulent to the highly organized forms of life. On the fundamental level, quantum fluctuations or entanglement lead to novel forms of complex dynamical behaviour in many-body systems for which a d… ▽ More Far-from-equilibrium situations are ubiquitous in nature. They are responsible for a wealth of phenomena, which are not simple extensions of near-equilibrium properties, ranging from fluid flows turning turbulent to the highly organized forms of life. On the fundamental level, quantum fluctuations or entanglement lead to novel forms of complex dynamical behaviour in many-body systems for which a description as emergent phenomena can be found within the framework of quantum field theory. A central quantity in these efforts, containing all information about the measurable physical properties, is the quantum effective action. Though the problem of non-equilibrium quantum dynamics can be exactly formulated in terms of the quantum effective action, the solution is in general beyond capabilities of classical computers. In this work, we present a strategy to determine the non-equilibrium quantum effective action using analog quantum simulators, and demonstrate our method experimentally with a quasi one-dimensional spinor Bose gas out of equilibrium. Building on spatially resolved snapshots of the spin degree of freedom, we infer the quantum effective action up to fourth order in an expansion in one-particle irreducible correlation functions at equal times. We uncover a strong suppression of the irreducible four-vertex emerging at low momenta, which solves the problem of dynamics in the highly occupied regime far from equilibrium where perturbative descriptions fail. Similar behaviour in this non-pertubative regime has been proposed in the context of early-universe cosmology. Our work constitutes a new realm of large-scale analog quantum computing, where the high level of control of synthetic quantum systems provides the means for the solution of long-standing theoretical problems in high-energy and condensed matter physics with an experimental approach. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Journal ref: Nature Physics 16, 1012-1016 (2020)

Non-cancellation of the parity anomaly in the strong-field regime of QED$_{2+1}$

Authors: Robert Ott, Torsten V. Zache, Niklas Mueller, Jürgen Berges

Abstract: Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilso… ▽ More Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in $2+1$ space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 6 pages, 4 figures

Foamy Dark Matter from Monodromies

Authors: Jürgen Berges, Aleksandr Chatrchyan, Joerg Jaeckel

Abstract: We investigate the dynamics of axion-like particle (ALP) dark matter where the field range is enlarged by a monodromy. The monodromy potential allows sufficient production of dark matter also at larger couplings to the Standard Model particles. The potential typically features a number of "wiggles" that lead to a rapid growth of fluctuations. Using classical-statistical field theory simulations we… ▽ More We investigate the dynamics of axion-like particle (ALP) dark matter where the field range is enlarged by a monodromy. The monodromy potential allows sufficient production of dark matter also at larger couplings to the Standard Model particles. The potential typically features a number of "wiggles" that lead to a rapid growth of fluctuations. Using classical-statistical field theory simulations we go beyond the linear regime and treat the system in the non-linear and even non-perturbative regime. For sufficiently strong wiggles the initially homogeneous field is completely converted into fluctuations. The fluctuations correspond to dark matter particles with a non-vanishing velocity and we consider the corresponding restrictions from structure formation as well as the effects on today's dark matter density. Since all the dark matter is made up from these strong fluctuations, the dark matter density features large, $\mathcal{O}(1)$ fluctuations at scales $\lesssim 10^{6}\,{\rm km}\sqrt{{\rm eV}/m_a}$. △ Less

Submitted 8 November, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 26 pages, 8 figures; v2: extended discussion about the field ending up in the true vs false minimum, extended discussion about the field becoming inhomogeneous vs staying homogeneous, Fig. 8 updated, clarifications and references added, published version

Entanglement and thermalization

Authors: Jürgen Berges, Stefan Floerchinger, Raju Venugopalan

Abstract: In a quantum field theory, apparent thermalization can be a consequence of entanglement as opposed to scatterings. We discuss here how this can help to explain open puzzles such as the success of thermal models in electron-positron collisions. It turns out that an expanding relativistic string described by the Schwinger model (which also underlies the Lund model) has at early times an entanglement… ▽ More In a quantum field theory, apparent thermalization can be a consequence of entanglement as opposed to scatterings. We discuss here how this can help to explain open puzzles such as the success of thermal models in electron-positron collisions. It turns out that an expanding relativistic string described by the Schwinger model (which also underlies the Lund model) has at early times an entanglement entropy that is extensive in rapidity. At these early times, the reduced density operator is of thermal form, with an entanglement temperature $T_τ=\hbar/(2πk_Bτ)$, even in the absence of any scatterings. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: Contribution to the proceedings of Quark Matter 2018, to appear in Nuclear Physics A

Schwinger pair production and string breaking in non-Abelian gauge theory from real-time lattice improved Hamiltonians

Authors: D. Spitz, J. Berges

Abstract: Far-from-equilibrium dynamics of SU(2) gauge theory with Wilson fermions is studied in 1+1 space-time dimensions using a real-time lattice approach. Lattice improved Hamiltonians are shown to be very efficient in simulating Schwinger pair creation and emergent phenomena such as plasma oscillations. As a consequence, significantly smaller lattices can be employed to approach continuum physics in th… ▽ More Far-from-equilibrium dynamics of SU(2) gauge theory with Wilson fermions is studied in 1+1 space-time dimensions using a real-time lattice approach. Lattice improved Hamiltonians are shown to be very efficient in simulating Schwinger pair creation and emergent phenomena such as plasma oscillations. As a consequence, significantly smaller lattices can be employed to approach continuum physics in the infinite-volume limit as compared to unimproved implementations. This allows us to compute also higher-order correlation functions including four fermion fields, which give unprecedented insights into the real-time dynamics of the fragmentation process of strings between fermions and antifermions. △ Less

Submitted 4 March, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 19 pages, 12 figures

Journal ref: Phys. Rev. D 99, 036020 (2019)

Breaking the Fluctuation-Dissipation Relation by Universal Transport Processes

Authors: A. Piñeiro Orioli, J. Berges

Abstract: Universal phenomena far from equilibrium exhibit additional independent scaling exponents and functions as compared to thermal universal behavior. For the example of an ultracold Bose gas we simulate nonequilibrium transport processes in a universal scaling regime and show how they lead to the breaking of the fluctuation-dissipation relation. As a consequence, the scaling of spectral functions (co… ▽ More Universal phenomena far from equilibrium exhibit additional independent scaling exponents and functions as compared to thermal universal behavior. For the example of an ultracold Bose gas we simulate nonequilibrium transport processes in a universal scaling regime and show how they lead to the breaking of the fluctuation-dissipation relation. As a consequence, the scaling of spectral functions (commutators) and statistical correlations (anticommutators) between different points in time and space become linearly independent with distinct dynamic scaling exponents. As a macroscopic signature of this phenomenon we identify a transport peak in the statistical two-point correlator, which is absent in the spectral function showing the quasiparticle peaks of the Bose gas. △ Less

Submitted 3 May, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 122, 150401 (2019)

Prescaling and far-from-equilibrium hydrodynamics in the quark-gluon plasma

Authors: Aleksas Mazeliauskas, Jürgen Berges

Abstract: Prescaling is a far-from-equilibrium phenomenon which describes the rapid establishment of a universal scaling form of distributions much before the universal values of their scaling exponents are realized. We consider the example of the spatio-temporal evolution of the quark-gluon plasma explored in heavy-ion collisions at sufficiently high energies. Solving QCD kinetic theory with elastic and in… ▽ More Prescaling is a far-from-equilibrium phenomenon which describes the rapid establishment of a universal scaling form of distributions much before the universal values of their scaling exponents are realized. We consider the example of the spatio-temporal evolution of the quark-gluon plasma explored in heavy-ion collisions at sufficiently high energies. Solving QCD kinetic theory with elastic and inelastic processes, we demonstrate that the gluon and quark distributions very quickly adapt a self-similar scaling form, which is independent of initial condition details and system parameters. The dynamics in the prescaling regime is then fully encoded in a few time-dependent scaling exponents, whose slow evolution gives rise to far-from-equilibrium hydrodynamic behavior. △ Less

Submitted 26 April, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 6 pages, 6 figures, v2 rewrote conclusion, other minor changes, published version

Journal ref: Phys. Rev. Lett. 122, 122301 (2019)

Dynamical topological transitions in the massive Schwinger model with a θ-term

Authors: T. V. Zache, N. Mueller, J. T. Schneider, F. Jendrzejewski, J. Berges, P. Hauke

Abstract: Aiming at a better understanding of anomalous and topological effects in gauge theories out-of-equilibrium, we study the real-time dynamics of a prototype model for CP-violation, the massive Schwinger model with a $θ$-term. We identify dynamical quantum phase transitions between different topological sectors that appear after sufficiently strong quenches of the $θ$-parameter. Moreover, we establis… ▽ More Aiming at a better understanding of anomalous and topological effects in gauge theories out-of-equilibrium, we study the real-time dynamics of a prototype model for CP-violation, the massive Schwinger model with a $θ$-term. We identify dynamical quantum phase transitions between different topological sectors that appear after sufficiently strong quenches of the $θ$-parameter. Moreover, we establish a general dynamical topological order parameter, which can be accessed through fermion two-point correlators and, importantly, which can be applied for interacting theories. Enabled by this result, we show that the topological transitions persist beyond the weak-coupling regime. Finally, these effects can be observed with table-top experiments based on existing cold-atom, superconducting-qubit, and trapped-ion technology. Our work, thus, presents a significant step towards quantum simulating topological and anomalous real-time phenomena relevant to nuclear and high-energy physics. △ Less

Submitted 23 August, 2018; originally announced August 2018.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 122, 050403 (2019)

Nonequilibrium quark production in the expanding QCD plasma

Authors: Naoto Tanji, Juergen Berges

Abstract: We present real-time lattice simulation results for nonequilibrium quark production from an over-occupied gluon plasma in longitudinally expanding geometry. The quark number density per unit transverse area and rapidity shows almost linear growth in time, and its growth rate appears to be consistent with a simple kinetic theory estimate involving only two-to-two scattering processes in small-angle… ▽ More We present real-time lattice simulation results for nonequilibrium quark production from an over-occupied gluon plasma in longitudinally expanding geometry. The quark number density per unit transverse area and rapidity shows almost linear growth in time, and its growth rate appears to be consistent with a simple kinetic theory estimate involving only two-to-two scattering processes in small-angle approximation. We also find that quarks produced at early times satisfy a nonequilibrium scaling law. △ Less

Submitted 13 July, 2018; originally announced July 2018.

Comments: Quark Matter 2018 proceedings

Non-equilibrium Quantum Spin Dynamics from 2PI Functional Integral Techniques in the Schwinger Boson Representation

Authors: A. Schuckert, A. Piñeiro Orioli, J. Berges

Abstract: We present a non-equilibrium quantum field theory approach to the initial-state dynamics of spin models based on two-particle irreducible (2PI) functional integral techniques. It employs a mapping of spins to Schwinger bosons for arbitrary spin interactions and spin lengths. At next-to-leading order (NLO) in an expansion in the number of field components, a wide range of non-perturbative dynamical… ▽ More We present a non-equilibrium quantum field theory approach to the initial-state dynamics of spin models based on two-particle irreducible (2PI) functional integral techniques. It employs a mapping of spins to Schwinger bosons for arbitrary spin interactions and spin lengths. At next-to-leading order (NLO) in an expansion in the number of field components, a wide range of non-perturbative dynamical phenomena are shown to be captured, including relaxation of magnetizations in a 3D long-range interacting system with quenched disorder, different relaxation behaviour on both sides of a quantum phase transition and the crossover from relaxation to arrest of dynamics in a disordered spin chain previously shown to exhibit many-body-localization. Where applicable, we employ alternative state-of-the-art techniques and find rather good agreement with our 2PI NLO results. As our method can handle large system sizes and converges relatively quickly to its thermodynamic limit, it opens the possibility to study those phenomena in higher dimensions in regimes in which no other efficient methods exist. Furthermore, the approach to classical dynamics can be investigated as the spin length is increased. △ Less

Submitted 11 December, 2018; v1 submitted 6 June, 2018; originally announced June 2018.

Journal ref: Phys. Rev. B 98, 224304 (2018)

Observation of universal dynamics in an isolated one-dimensional Bose gas far from equilibrium

Authors: S. Erne, R. Buecker, T. Gasenzer, J. Berges, J. Schmiedmayer

Abstract: We provide experimental evidence of universal dynamics far from equilibrium during the relaxation of an isolated one-dimensional Bose gas. Following a rapid cooling quench, the system exhibits universal scaling in time and space, associated with the approach of a non-thermal fixed point. The time evolution within the scaling period is described by a single universal function and scaling exponent,… ▽ More We provide experimental evidence of universal dynamics far from equilibrium during the relaxation of an isolated one-dimensional Bose gas. Following a rapid cooling quench, the system exhibits universal scaling in time and space, associated with the approach of a non-thermal fixed point. The time evolution within the scaling period is described by a single universal function and scaling exponent, independent of the specifics of the initial state. Our results provide a quantum simulation in a regime, where to date no theoretical predictions are available. This constitutes a crucial step in the verification of universality far from equilibrium. If successful, this may lead to a comprehensive classification of systems based on their universal properties far from equilibrium, relevant for a large variety of systems at different scales. △ Less

Submitted 30 May, 2018; originally announced May 2018.

Comments: 12 pages, 9 figures

Journal ref: Nature 563, 225 (2018)

Observation of universal dynamics in a spinor Bose gas far from equilibrium

Authors: Maximilian Prüfer, Philipp Kunkel, Helmut Strobel, Stefan Lannig, Daniel Linnemann, Christian-Marcel Schmied, Jürgen Berges, Thomas Gasenzer, Markus K. Oberthaler

Abstract: The dynamics of quantum systems far from equilibrium represents one of the most challenging problems in theoretical many-body physics. While the evolution is in general intractable in all its details, relevant observables can become insensitive to microscopic system parameters and initial conditions. This is the basis of the phenomenon of universality. Far from equilibrium, universality is identif… ▽ More The dynamics of quantum systems far from equilibrium represents one of the most challenging problems in theoretical many-body physics. While the evolution is in general intractable in all its details, relevant observables can become insensitive to microscopic system parameters and initial conditions. This is the basis of the phenomenon of universality. Far from equilibrium, universality is identified through the scaling of the spatiotemporal evolution of the system, captured by universal exponents and functions. Theoretically, this has been studied in examples as different as the reheating process in inflationary universe cosmology, the dynamics of nuclear collision experiments described by quantum chromodynamics, or the post-quench dynamics in dilute quantum gases in non-relativistic quantum field theory. Here we observe the emergence of universal dynamics by evaluating spatially resolved spin correlations in a quasi one-dimensional spinor Bose-Einstein condensate. For long evolution times we extract the scaling properties from the spatial correlations of the spin excitations. From this we find the dynamics to be governed by transport of an emergent conserved quantity towards low momentum scales. Our results establish an important class of non-stationary systems whose dynamics is encoded in time independent scaling exponents and functions signaling the existence of non-thermal fixed points. We confirm that the non-thermal scaling phenomenon involves no fine-tuning, by preparing different initial conditions and observing the same scaling behaviour. Our analog quantum simulation approach provides the basis to reveal the underlying mechanisms and characteristics of non-thermal universality classes. One may use this universality to learn, from experiments with ultra-cold gases, about fundamental aspects of dynamics studied in cosmology and quantum chromodynamics. △ Less

Submitted 6 December, 2018; v1 submitted 30 May, 2018; originally announced May 2018.

Comments: 10 pages, 4 figures, 3 supplementary figures

Journal ref: Nature 563, 217-220 (2018)

Electronic structure of single layer 1T-NbSe$_2$: interplay of lattice distortions, non-local exchange, and Mott-Hubbard correlations

Authors: E. Kamil, J. Berges, G. Schönhoff, M. Rösner, M. Schüler, G. Sangiovanni, T. O. Wehling

Abstract: Using ab-initio calculations we reveal the nature of the insulating phase recently found experimentally in monolayer 1T-NbSe$_2$. We find soft phonon modes in a large parts of the Brillouin zone indicating the strong-coupling nature of a charge-density-wave instability. Structural relaxation of a $\sqrt{13}\times\sqrt{13}$ supercell reveals a Star-of-David reconstruction with an energy gain of 60… ▽ More Using ab-initio calculations we reveal the nature of the insulating phase recently found experimentally in monolayer 1T-NbSe$_2$. We find soft phonon modes in a large parts of the Brillouin zone indicating the strong-coupling nature of a charge-density-wave instability. Structural relaxation of a $\sqrt{13}\times\sqrt{13}$ supercell reveals a Star-of-David reconstruction with an energy gain of 60 meV per primitive unit cell. The band structure of the distorted phase exhibits a half-filled flat band which is associated with orbitals that are delocalized over several atoms in each Star of David. By including many-body corrections through a combined GW, hybrid-functional, and DMFT treatment, we find the flat band to split into narrow Hubbard bands. The lowest energy excitation across the gap turns out to be between itinerant Se-$p$ states and the upper Hubbard band, determining the system to be a charge-transfer insulator. Combined hybrid-functional and GW calculations show that long-range interactions shift the Se-$p$ states to lower energies. Thus, a delicate interplay of local and long-range correlations determines the gap nature and its size in this distorted phase of the monolayer 1T-NbSe$_2$. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 10 pages, 11 figures

Journal ref: E Kamil et al 2018 J. Phys.: Condens. Matter 30 325601

Plasmonic Superconductivity in Layered Materials

Authors: M. Rösner, R. E. Groenewald, G. Schönhoff, J. Berges, S. Haas, T. O. Wehling

Abstract: Plasmonic excitations behave fundamentally different in layered materials in comparison to bulk systems. They form gapless modes, which in turn couple at low energies to the electrons. Thereby they can strongly influence superconducting instabilities. Here, we show how these excitations can be controlled from the outside via changes in the dielectric environment or in the doping level, which allow… ▽ More Plasmonic excitations behave fundamentally different in layered materials in comparison to bulk systems. They form gapless modes, which in turn couple at low energies to the electrons. Thereby they can strongly influence superconducting instabilities. Here, we show how these excitations can be controlled from the outside via changes in the dielectric environment or in the doping level, which allows for external tuning of the superconducting transition temperature. By solving the gap equation for an effective system, we find that the plasmonic influence can both strongly enhance or reduce the transition temperature, depending on the details of the plasmon-phonon interplay. We formulate simple experimental guidelines to find plasmon- induced elevated transition temperatures in layered materials. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 13 pages, 10 figures

Quantum simulation of lattice gauge theories using Wilson fermions

Authors: T. V. Zache, F. Hebenstreit, F. Jendrzejewski, M. K. Oberthaler, J. Berges, P. Hauke

Abstract: Quantum simulators have the exciting prospect of giving access to real-time dynamics of lattice gauge theories, in particular in regimes that are difficult to compute on classical computers. Future progress towards scalable quantum simulation of lattice gauge theories, however, hinges crucially on the efficient use of experimental resources. As we argue in this work, due to the fundamental non-uni… ▽ More Quantum simulators have the exciting prospect of giving access to real-time dynamics of lattice gauge theories, in particular in regimes that are difficult to compute on classical computers. Future progress towards scalable quantum simulation of lattice gauge theories, however, hinges crucially on the efficient use of experimental resources. As we argue in this work, due to the fundamental non-uniqueness of discretizing the relativistic Dirac Hamiltonian, the lattice representation of gauge theories allows for an optimization that up to now has been left unexplored. We exemplify our discussion with lattice quantum electrodynamics in two-dimensional space-time, where we show that the formulation through Wilson fermions provides several advantages over the previously considered staggered fermions. Notably, it enables a strongly simplified optical lattice setup and it reduces the number of degrees of freedom required to simulate dynamical gauge fields. Exploiting the optimal representation, we propose an experiment based on a mixture of ultracold atoms trapped in a tilted optical lattice. Using numerical benchmark simulations, we demonstrate that a state-of-the-art quantum simulator may access the Schwinger mechanism and map out its non-perturbative onset. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 19 pages, 11 figures

Journal ref: Quantum Sci. Technol. 3 034010 (2018)

Dynamics of entanglement in expanding quantum fields

Authors: Jürgen Berges, Stefan Floerchinger, Raju Venugopalan

Abstract: We develop a novel real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in… ▽ More We develop a novel real-time approach to computing the entanglement between spatial regions for Gaussian states in quantum field theory. The entanglement entropy is characterized in terms of local correlation functions on space-like Cauchy hypersurfaces. The framework is applied to explore an expanding light cone geometry in the particular case of the Schwinger model for quantum electrodynamics in 1+1 space-time dimensions. We observe that the entanglement entropy becomes extensive in rapidity at early times and that the corresponding local reduced density matrix is a thermal density matrix for excitations around a coherent field with a time dependent temperature. Since the Schwinger model successfully describes many features of multiparticle production in $e^+ e^-$ collisions, our results provide an attractive explanation in this framework for the apparent thermal nature of multiparticle production even in the absence of significant final state scattering. △ Less

Submitted 26 December, 2017; originally announced December 2017.

Comments: 39 pages, 4 figures