Publications

Order by: publication date number of citations

• RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles

  Nicholas P. Bailey, Trond S. Ingebrigtsen, Jesper Schmidt Hansen, Arno A. Veldhorst, Lasse Bøhling, Claire A. Lemarchand, Andreas E. Olsen, Andreas K. Bacher, Lorenzo Costigliola, Ulf R. Pedersen, Heine Larsen, Jeppe C. Dyre, Thomas B. Schrøder

  SciPost Phys. 3, 038 (2017) · published 14 December 2017 | Toggle abstract · pdf

  RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPU's). RUMD addresses the challenge of utilizing the many-core nature of modern GPU hardware when simulating small to medium system sizes (roughly from a few thousand up to hundred thousand particles). It has a performance that is comparable to other GPU-MD codes at large system sizes and substantially better at smaller sizes.RUMD is open-source and consists of a library written in C++ and the CUDA extension to C, an easy-to-use Python interface, and a set of tools for set-up and post-simulation data analysis. The paper describes RUMD's main features, optimizations and performance benchmarks.

  64 citations

• Whispers from the dark side: Confronting light new physics with NANOGrav data

  Wolfram Ratzinger, Pedro Schwaller

  SciPost Phys. 10, 047 (2021) · published 23 February 2021 | Toggle abstract · pdf

  The NANOGrav collaboration has recently observed first evidence of a grav- itational wave background (GWB) in pulsar timing data. Here we explore the possibility that this GWB is due to new physics, and show that the signal can be well fit also with peaked spectra like the ones expected from phase transitions (PTs) or from the dynamics of axion like particles (ALPs) in the early universe. We find that a good fit to the data is obtained for a very strong PT at temperatures around 1 MeV to 10 MeV. For the ALP explanation the best fit is obtained for a decay constant of F ≈ 5 × 10^17 GeV and an axion mass of 2 × 10^−13 eV. We also illustrate the ability of PTAs to constrain the parameter space of these models, and obtain limits which are already comparable to other cosmological bounds.

  64 citations

• Emergence of curved light-cones in a class of inhomogeneous Luttinger liquids

  Jérôme Dubail, Jean-Marie Stéphan, Pasquale Calabrese

  SciPost Phys. 3, 019 (2017) · published 6 September 2017 | Toggle abstract · pdf

  The light-cone spreading of entanglement and correlation is a fundamental and ubiquitous feature of homogeneous extended quantum systems. Here we point out that a class of inhomogenous Luttinger liquids (those with a uniform Luttinger parameter $K$) at low energy display the universal phenomenon of curved light cones: gapless excitations propagate along the geodesics of the metric $ds^2=dx^2+v(x)^2 d\tau^2$, with $v(x)$ being the calculable spatial dependent velocity induced by the inhomogeneity. We confirm our findings with explicit analytic and numerical calculations both in- and out-of-equilibrium for a Tonks-Girardeau gas in a harmonic potential and in lattice systems with artificially tuned hamiltonian density.

  64 citations

• Reconstructing the gluon

  Anton K. Cyrol, Jan M. Pawlowski, Alexander Rothkopf, Nicolas Wink

  SciPost Phys. 5, 065 (2018) · published 20 December 2018 | Toggle abstract · pdf

  We reconstruct the gluon spectral function in Landau gauge QCD from numerical data for the gluon propagator. The reconstruction relies on two novel ingredients: Firstly we derive analytically the low frequency asymptotics of the spectral function. Secondly we construct a functional basis from a careful consideration of the analytic properties of the gluon propagator in Landau gauge. This allows us to reliably capture the non-perturbative regime of the gluon spectrum. We also compare different reconstruction methods and discuss the respective systematic errors.

  63 citations

• Exotic $\mathbb{Z}_N$ symmetries, duality, and fractons in 3+1-dimensional quantum field theory

  Nathan Seiberg, Shu-Heng Shao

  SciPost Phys. 10, 003 (2021) · published 8 January 2021 | Toggle abstract · pdf

  Following our earlier analyses of nonstandard continuum quantum field theories, we study here gapped systems in 3+1 dimensions, which exhibit fractonic behavior. In particular, we present three dual field theory descriptions of the low-energy physics of the X-cube model. A key aspect of our constructions is the use of discontinuous fields in the continuum field theory. Spacetime is continuous, but the fields are not.

  62 citations

• Effective universality in quantum gravity

  Astrid Eichhorn, Peter Labus, Jan M. Pawlowski, Manuel Reichert

  SciPost Phys. 5, 031 (2018) · published 4 October 2018 | Toggle abstract · pdf

  We investigate the asymptotic safety scenario for a scalar-gravity system. This system contains two avatars of the dynamical Newton coupling, a gravitational self-coupling and a scalar-graviton coupling. We uncover an effective universality for the dynamical Newton coupling on the quantum level: its momentum-dependent avatars are in remarkable quantitative agreement in the scaling regime of the UV fixed point. For the background Newton coupling, this effective universality is not present, but qualitative agreement remains.

  62 citations

• Symmetry decomposition of negativity of massless free fermions

  Sara Murciano, Riccarda Bonsignori, Pasquale Calabrese

  SciPost Phys. 10, 111 (2021) · published 19 May 2021 | Toggle abstract · pdf

  We consider the problem of symmetry decomposition of the entanglement negativity in free fermionic systems. Rather than performing the standard partial transpose, we use the partial time-reversal transformation which naturally encodes the fermionic statistics. The negativity admits a resolution in terms of the charge imbalance between the two subsystems. We introduce a normalised version of the imbalance resolved negativity which has the advantage to be an entanglement proxy for each symmetry sector, but may diverge in the limit of pure states for some sectors. Our main focus is then the resolution of the negativity for a free Dirac field at finite temperature and size. We consider both bipartite and tripartite geometries and exploit conformal field theory to derive universal results for the charge imbalance resolved negativity. To this end, we use a geometrical construction in terms of an Aharonov-Bohm-like flux inserted in the Riemann surface defining the entanglement. We interestingly find that the entanglement negativity is always equally distributed among the different imbalance sectors at leading order. Our analytical findings are tested against exact numerical calculations for free fermions on a lattice.

  61 citations

• Entanglement evolution and generalised hydrodynamics: interacting integrable systems

  Vincenzo Alba, Bruno Bertini, Maurizio Fagotti

  SciPost Phys. 7, 005 (2019) · published 8 July 2019 | Toggle abstract · pdf

  We investigate the dynamics of bipartite entanglement after the sudden junction of two leads in interacting integrable models. By combining the quasiparticle picture for the entanglement spreading with Generalised Hydrodynamics we derive an analytical prediction for the dynamics of the entanglement entropy between a finite subsystem and the rest. We find that the entanglement rate between the two leads depends only on the physics at the interface and differs from the rate of exchange of thermodynamic entropy. This contrasts with the behaviour in free or homogeneous interacting integrable systems, where the two rates coincide.

  60 citations

• The ALF (Algorithms for Lattice Fermions) project release 1.0. Documentation for the auxiliary field quantum Monte Carlo code

  Martin Bercx, Florian Goth, Johannes S. Hofmann, Fakher F. Assaad

  SciPost Phys. 3, 013 (2017) · published 16 August 2017 · licensed under CC BY-SA (4.0) | Toggle abstract · pdf

  The Algorithms for Lattice Fermions package provides a general code for the finite temperature auxiliary field quantum Monte Carlo algorithm. The code is engineered to be able to simulate any model that can be written in terms of sums of single-body operators, of squares of single-body operators and single-body operators coupled to an Ising field with given dynamics. We provide predefined types that allow the user to specify the model, the Bravais lattice as well as equal time and time displaced observables. The code supports an MPI implementation. Examples such as the Hubbard model on the honeycomb lattice and the Hubbard model on the square lattice coupled to a transverse Ising field are provided and discussed in the documentation. We furthermore discuss how to use the package to implement the Kondo lattice model and the $SU(N)$-Hubbard-Heisenberg model. One can download the code from our Git instance at https://alf.physik.uni-wuerzburg.de and sign in to file issues.

  59 citations

• Quantum quenches to the attractive one-dimensional Bose gas: exact results

  Lorenzo Piroli, Pasquale Calabrese, Fabian H. L. Essler

  SciPost Phys. 1, 001 (2016) · published 14 September 2016 | Toggle abstract · pdf

  We study quantum quenches to the one-dimensional Bose gas with attractive interactions in the case when the initial state is an ideal one-dimensional Bose condensate. We focus on properties of the stationary state reached at late times after the quench. This displays a finite density of multi-particle bound states, whose rapidity distribution is determined exactly by means of the quench action method. We discuss the relevance of the multi-particle bound states for the physical properties of the system, computing in particular the stationary value of the local pair correlation function $g_2$.

  59 citations