Publications

Order by: publication date number of citations

• Introducing iFluid: a numerical framework for solving hydrodynamical equations in integrable models

  Frederik S. Møller, Jörg Schmiedmayer

  SciPost Phys. 8, 041 (2020) · published 13 March 2020 | Toggle abstract · pdf

  We present an open-source Matlab framework, titled iFluid, for simulating the dynamics of integrable models using the theory of generalized hydrodynamics (GHD). The framework provides an intuitive interface, enabling users to define and solve problems in a few lines of code. Moreover, iFluid can be extended to encompass any integrable model, and the algorithms for solving the GHD equations can be fully customized. We demonstrate how to use iFluid by solving the dynamics of three distinct systems: (i) The quantum Newton's cradle of the Lieb-Liniger model, (ii) a gradual field release in the XXZ-chain, and (iii) a partitioning protocol in the relativistic sinh-Gordon model.

  32 citations

• Universal bounds on transport in holographic systems with broken translations

  Matteo Baggioli, Wei-Jia Li

  SciPost Phys. 9, 007 (2020) · published 15 July 2020 | Toggle abstract · pdf

  We study the presence of universal bounds on transport in homogeneous holographic models with broken translations. We verify numerically that, in holographic systems with momentum dissipation, the viscosity to entropy bound might be violated but the shear diffusion constant remains bounded by below. This confirms the idea that $\eta/s$ loses its privileged role in non-relativistic systems and that, in order to find more universal bounds, one should rather look at diffusion constants. We strengthen this idea by showing that, in presence of spontaneously broken translations, the Goldstone diffusion constant satisfies a universal lower bound in terms of the Planckian relaxation time and the butterfly velocity. Additionally, all the diffusive processes in the model satisfy an upper bound, imposed by causality, which is given in terms of the thermalization time -- the imaginary part of the first non-hydrodynamic mode in the spectrum -- and the speed of longitudinal sound. Finally, we discuss the existence of a bound on the speed of sound in holographic conformal solids and we show that the conformal value acts as a lower (and not upper) bound on the speed of longitudinal phonons. Nevertheless, we show that the stiffness $\partial p/\partial \epsilon$ is still bounded by above by its conformal value. This suggests that the bounds conjectured in the past have to be considered on the stiffness of the system, related to its equation of state, and not on the propagation speed of sound.

  32 citations

• NetKet 3: Machine Learning Toolbox for Many-Body Quantum Systems

  Filippo Vicentini, Damian Hofmann, Attila Szabó, Dian Wu, Christopher Roth, Clemens Giuliani, Gabriel Pescia, Jannes Nys, Vladimir Vargas-Calderón, Nikita Astrakhantsev, Giuseppe Carleo

  SciPost Phys. Codebases 7 (2022) · published 24 August 2022 | Toggle abstract · pdf

  We introduce version 3 of NetKet, the machine learning toolbox for many-body quantum physics. NetKet is built around neural-network quantum states and provides efficient algorithms for their evaluation and optimization. This new version is built on top of JAX, a differentiable programming and accelerated linear algebra framework for the Python programming language. The most significant new feature is the possibility to define arbitrary neural network ansätze in pure Python code using the concise notation of machine-learning frameworks, which allows for just-in-time compilation as well as the implicit generation of gradients thanks to automatic differentiation. NetKet 3 also comes with support for GPU and TPU accelerators, advanced support for discrete symmetry groups, chunking to scale up to thousands of degrees of freedom, drivers for quantum dynamics applications, and improved modularity, allowing users to use only parts of the toolbox as a foundation for their own code.

  32 citations

• Carroll Expansion of General Relativity

  Dennis Hansen, Niels A. Obers, Gerben Oling, Benjamin T. Søgaard

  SciPost Phys. 13, 055 (2022) · published 8 September 2022 | Toggle abstract · pdf

  We study the small speed of light expansion of general relativity, utilizing the modern perspective on non-Lorentzian geometry. This is an expansion around the ultra-local Carroll limit, in which light cones close up. To this end, we first rewrite the Einstein-Hilbert action in pre-ultra-local variables, which is closely related to the 3+1 decomposition of general relativity. At leading order in the expansion, these pre-ultra-local variables yield Carroll geometry and the resulting action describes the electric Carroll limit of general relativity. We also obtain the next-to-leading order action in terms of Carroll geometry and next-to-leading order geometric fields. The leading order theory yields constraint and evolution equations, and we can solve the evolution analytically. We furthermore construct a Carroll version of Bowen-York initial data, which has associated conserved boundary linear and angular momentum charges. The notion of mass is not present at leading order and only enters at next-to-leading order. This is illustrated by considering a particular truncation of the next-to-leading order action, corresponding to the magnetic Carroll limit, where we find a solution that describes the Carroll limit of a Schwarzschild black hole. Finally, we comment on how a cosmological constant can be incorporated in our analysis.

  32 citations

• Quench dynamics of the Ising field theory in a magnetic field

  Kristóf Hódsági, Márton Kormos, Gábor Takács

  SciPost Phys. 5, 027 (2018) · published 25 September 2018 | Toggle abstract · pdf

  We numerically simulate the time evolution of the Ising field theory after quenches starting from the $E_8$ integrable model using the Truncated Conformal Space Approach. The results are compared with two different analytic predictions based on form factor expansions in the pre-quench and post-quench basis, respectively. Our results clarify the domain of validity of these expansions and suggest directions for further improvement. We show for quenches in the $E_8$ model that the initial state is not of the integrable pair state form. We also construct quench overlap functions and show that their high-energy asymptotics are markedly different from those constructed before in the sinh/sine-Gordon theory, and argue that this is related to properties of the ultraviolet fixed point.

  32 citations

• Dissipation-induced topological insulators: A no-go theorem and a recipe

  Moshe Goldstein

  SciPost Phys. 7, 067 (2019) · published 27 November 2019 | Toggle abstract · pdf

  Nonequilibrium conditions are traditionally seen as detrimental to the appearance of quantum-coherent many-body phenomena, and much effort is often devoted to their elimination. Recently this approach has changed: It has been realized that driven-dissipative dynamics could be used as a resource. By proper engineering of the reservoirs and their couplings to a system, one may drive the system towards desired quantum-correlated steady states, even in the absence of internal Hamiltonian dynamics. An intriguing category of equilibrium many-particle phases are those which are distinguished by topology rather than by symmetry. A natural question thus arises: which of these topological states can be achieved as the result of dissipative Lindblad-type (Markovian) evolution? Beside its fundamental importance, it may offer novel routes to the realization of topologically-nontrivial states in quantum simulators, especially ultracold atomic gases. Here I give a general answer for Gaussian states and quadratic Lindblad evolution, mostly concentrating on 2D Chern insulator states. I prove a no-go theorem stating that a finite-range Lindbladian cannot induce finite-rate exponential decay towards a unique topological pure state above 1D. I construct a recipe for creating such state by exponentially-local dynamics, or a mixed state arbitrarily close to the desired pure one via finite-range dynamics. I also address the cold-atom realization, classification, and detection of these states. Extensions to other types of topological insulators and superconductors are also discussed.

  32 citations

• Anomaly matching in the symmetry broken phase: Domain walls, CPT, and the Smith isomorphism

  Itamar Hason, Zohar Komargodski, Ryan Thorngren

  SciPost Phys. 8, 062 (2020) · published 17 April 2020 | Toggle abstract · pdf

  Symmetries in Quantum Field Theory may have 't Hooft anomalies. If the symmetry is unbroken in the vacuum, the anomaly implies a nontrivial low-energy limit, such as gapless modes or a topological field theory. If the symmetry is spontaneously broken, for the continuous case, the anomaly implies low-energy theorems about certain couplings of the Goldstone modes. Here we study the case of spontaneously broken discrete symmetries, such as Z/2 and T. Symmetry breaking leads to domain walls, and the physics of the domain walls is constrained by the anomaly. We investigate how the physics of the domain walls leads to a matching of the original discrete anomaly. We analyze the symmetry structure on the domain wall, which requires a careful analysis of some properties of the unbreakable CPT symmetry. We demonstrate the general results on some examples and we explain in detail the mod 4 periodic structure that arises in the Z/2 and T case. This gives a physical interpretation for the Smith isomorphism, which we also extend to more general abelian groups. We show that via symmetry breaking and the analysis of the physics on the wall, the computations of certain discrete anomalies are greatly simplified. Using these results we perform new consistency checks on the infrared phases of 2+1 dimensional QCD.

  32 citations

• Riemannian optimization of isometric tensor networks

  Markus Hauru, Maarten Van Damme, Jutho Haegeman

  SciPost Phys. 10, 040 (2021) · published 17 February 2021 | Toggle abstract · pdf

  Several tensor networks are built of isometric tensors, i.e. tensors satisfying $W^\dagger W = \mathrm{I}$. Prominent examples include matrix product states (MPS) in canonical form, the multiscale entanglement renormalization ansatz (MERA), and quantum circuits in general, such as those needed in state preparation and quantum variational eigensolvers. We show how gradient-based optimization methods on Riemannian manifolds can be used to optimize tensor networks of isometries to represent e.g. ground states of 1D quantum Hamiltonians. We discuss the geometry of Grassmann and Stiefel manifolds, the Riemannian manifolds of isometric tensors, and review how state-of-the-art optimization methods like nonlinear conjugate gradient and quasi-Newton algorithms can be implemented in this context. We apply these methods in the context of infinite MPS and MERA, and show benchmark results in which they outperform the best previously-known optimization methods, which are tailor-made for those specific variational classes. We also provide open-source implementations of our algorithms.

  31 citations

• Domain wall problem in the quantum XXZ chain and semiclassical behavior close to the isotropic point

  Grégoire Misguich, Nicolas Pavloff, Vincent Pasquier

  SciPost Phys. 7, 025 (2019) · published 29 August 2019 | Toggle abstract · pdf

  We study the dynamics of a spin-1/2 XXZ chain which is initially prepared in a domain-wall state. We compare the results of time-dependent Density Matrix Renormalization Group simulations with those of an effective description in terms of a classical anisotropic Landau-Lifshitz (LL) equation. Numerous quantities are analyzed: magnetization (x, y and z components), energy density, energy current, but also some spin-spin correlation functions or entanglement entropy in the quantum chain. Without any adjustable parameter a quantitative agreement is observed between the quantum and the LL problems in the long time limit, when the models are close to the isotropic point. This is explained as a consequence of energy conservation. At the isotropic point the mapping between the LL equation and the nonlinear Schr\"odinger equation is used to construct a variational solution capturing several aspects of the problem.

  31 citations

• Investigation of the Néel phase of the frustrated Heisenberg antiferromagnet by differentiable symmetric tensor networks

  Juraj Hasik, Didier Poilblanc, Federico Becca

  SciPost Phys. 10, 012 (2021) · published 20 January 2021 | Toggle abstract · pdf

  The recent progress in the optimization of two-dimensional tensor networks [H.-J. Liao, J.-G. Liu, L. Wang, and T. Xiang, Phys. Rev. X 9, 031041 (2019)] based on automatic differentiation opened the way towards precise and fast optimization of such states and, in particular, infinite projected entangled-pair states (iPEPS) that constitute a generic-purpose Ansatz for lattice problems governed by local Hamiltonians. In this work, we perform an extensive study of a paradigmatic model of frustrated magnetism, the J 1 − J 2 Heisenberg an- tiferromagnet on the square lattice. By using advances in both optimization and subsequent data analysis, through finite correlation-length scaling, we re- port accurate estimations of the magnetization curve in the Néel phase for J 2 /J 1 ≤ 0.45. The unrestricted iPEPS simulations reveal an U (1) symmetric structure, which we identify and impose on tensors, resulting in a clean and consistent picture of antiferromagnetic order vanishing at the phase transition with a quantum paramagnet at J 2 /J 1 ≈ 0.46(1). The present methodology can be extended beyond this model to study generic order-to-disorder transitions in magnetic systems.

  31 citations