Publications

Order by: publication date number of citations

• Disorder-free spin glass transitions and jamming in exactly solvable mean-field models

  Hajime Yoshino

  SciPost Phys. 4, 040 (2018) · published 26 June 2018 | Toggle abstract · pdf

  We construct and analyze a family of $M$-component vectorial spin systems which exhibit glass transitions and jamming within supercooled paramagnetic states without quenched disorder. Our system is defined on lattices with connectivity $c=\alpha M$ and becomes exactly solvable in the limit of large number of components $M \to \infty$. We consider generic $p$-body interactions between the vectorial Ising/continuous spins with linear/non-linear potentials. The existence of self-generated randomness is demonstrated by showing that the random energy model is recovered from a $M$-component ferromagnetic $p$-spin Ising model in $M \to \infty$ and $p \to \infty$ limit. In our systems the quenched disorder, if present, and the self-generated disorder act additively. Our theory provides a unified mean-field theoretical framework for glass transitions of rotational degree of freedoms such as orientation of molecules in glass forming liquids, color angles in continuous coloring of graphs and vector spins of geometrically frustrated magnets. The rotational glass transitions accompany various types of replica symmetry breaking. In the case of repulsive hardcore interactions in the spin space, continuous the criticality of the jamming or SAT/UNSTAT transition becomes the same as that of hardspheres.

  15 citations

• Generating dense packings of hard spheres by soft interaction design

  Thibaud Maimbourg, Mauro Sellitto, Guilhem Semerjian, Francesco Zamponi

  SciPost Phys. 4, 039 (2018) · published 26 June 2018 | Toggle abstract · pdf

  Packing spheres efficiently in large dimension $d$ is a particularly difficult optimization problem. In this paper we add an isotropic interaction potential to the pure hard-core repulsion, and show that one can tune it in order to maximize a lower bound on packing density. Our results suggest that exponentially many (in the number of particles) distinct disordered sphere packings can be effectively constructed by this method, up to a packing fraction close to $7\, d\, 2^{-d}$. The latter is determined by solving the inverse problem of maximizing the dynamical glass transition over the space of the interaction potentials. Our method crucially exploits a recent exact formulation of the thermodynamics and the dynamics of simple liquids in infinite dimension.

  7 citations

• The Inhomogeneous Gaussian Free Field, with application to ground state correlations of trapped 1d Bose gases

  Yannis Brun, Jérôme Dubail

  SciPost Phys. 4, 037 (2018) · published 25 June 2018 | Toggle abstract · pdf

  Motivated by the calculation of correlation functions in inhomogeneous one-dimensional (1d) quantum systems, the 2d Inhomogeneous Gaussian Free Field (IGFF) is studied and solved. The IGFF is defined in a domain $\Omega \subset \mathbb{R}^2$ equipped with a conformal class of metrics $[{\rm g}]$ and with a real positive coupling constant $K: \Omega \rightarrow \mathbb{R}_{>0}$ by the action $\mathcal{S}[h] = \frac{1}{8\pi } \int_\Omega \frac{\sqrt{{\rm g}} d^2 {\rm x}}{K({\rm x})} \, {\rm g}^{i j} (\partial_i h)(\partial_j h)$. All correlations functions of the IGFF are expressible in terms of the Green's functions of generalized Poisson operators that are familiar from 2d electrostatics in media with spatially varying dielectric constants. This formalism is then applied to the study of ground state correlations of the Lieb-Liniger gas trapped in an external potential $V(x)$. Relations with previous works on inhomogeneous Luttinger liquids are discussed. The main innovation here is in the identification of local observables $\hat{O} (x)$ in the microscopic model with their field theory counterparts $\partial_x h, e^{i h(x)}, e^{-i h(x)}$, etc., which involve non-universal coefficients that themselves depend on position --- a fact that, to the best of our knowledge, was overlooked in previous works on correlation functions of inhomogeneous Luttinger liquids ---, and that can be calculated thanks to Bethe Ansatz form factors formulae available for the homogeneous Lieb-Liniger model. Combining those position-dependent coefficients with the correlation functions of the IGFF, ground state correlation functions of the trapped gas are obtained. Numerical checks from DMRG are provided for density-density correlations and for the one-particle density matrix, showing excellent agreement.

  44 citations

• Lecture notes on Diagrammatic Monte Carlo for the Fröhlich polaron

  Jonas Greitemann, Lode Pollet

  SciPost Phys. Lect. Notes 2 (2018) · published 11 June 2018 | Toggle abstract · pdf

  These notes are intended as a detailed discussion on how to implement the diagrammatic Monte Carlo method for a physical system which is technically simple and where it works extremely well, namely the Fröhlich polaron problem. Sampling schemes for the Green function as well as the self-energy in the bare and skeleton (bold) expansion are disclosed in full detail. We discuss the Monte Carlo updates, possible implementations in terms of common data structures, as well as techniques on how to perform the Fourier transforms for functions with discontinuities. Control over the variety of parameters, especially in the bold scheme, is demonstrated. Sample codes are made available online along with extensive documentation. Towards the end, we discuss various extensions of the method and their applications. After working through these notes, the reader will be well equipped to explore the richness of the diagrammatic Monte Carlo method for quantum many-body systems.

  9 citations

• Liu-Nagel phase diagrams in infinite dimension

  Giulio Biroli, Pierfrancesco Urbani

  SciPost Phys. 4, 020 (2018) · published 27 April 2018 | Toggle abstract · pdf

  We study Harmonic Soft Spheres as a model of thermal structural glasses in the limit of infinite dimensions. We show that cooling, compressing and shearing a glass lead to a Gardner transition and, hence, to a marginally stable amorphous solid as found for Hard Spheres systems. A general outcome of our results is that a reduced stability of the glass favors the appearance of the Gardner transition. Therefore using strong perturbations, e.g. shear and compression, on standard glasses or using weak perturbations on weakly stable glasses, e.g. the ones prepared close to the jamming point, are the generic ways to induce a Gardner transition. The formalism that we discuss allows to study general perturbations, including strain deformations that are important to study soft glassy rheology at the mean field level.

  23 citations

• Time evolution of effective central charge and signatures of RG irreversibility after a quantum quench

  Axel Cortes Cubero

  SciPost Phys. 4, 016 (2018) · published 27 March 2018 | Toggle abstract · pdf

  At thermal equilibrium, the concept of effective central charge for massive deformations of two-dimensional conformal field theories (CFT) is well understood, and can be defined by comparing the partition function of the massive model to that of a CFT. This temperature-dependent effective charge interpolates monotonically between the central charge values corresponding to the IR and UV fixed points at low and high temperatures, respectively. We propose a non-equilibrium, time-dependent generalization of the effective central charge for integrable models after a quantum quench, $c_{\rm eff}(t)$, obtained by comparing the return amplitude to that of a CFT quench. We study this proposal for a large mass quench of a free boson, where the charge is seen to interpolate between $c_{\rm eff}=0$ at $t=0$, and $c_{\rm eff}\sim 1$ at $t\to\infty$, as is expected. We use our effective charge to define an "Ising to Tricritical Ising" quench protocol, where the charge evolves from $c_{\rm eff}=1/2$ at $t=0$, to $c_{\rm eff}=7/10$ at $t\to\infty$, the corresponding values of the first two unitary minimal CFT models. We then argue that the inverse "Tricritical Ising to Ising" quench is impossible with our methods. These conclusions can be generalized for quenches between any two adjacent unitary minimal CFT models. We finally study a large mass quench into the "staircase model" (sinh-Gordon with a particular complex coupling). At short times after the quench, the effective central charge increases in a discrete "staircase" structure, where the values of the charge at the steps can be computed in terms of the central charges of unitary minimal CFT models. When the initial state is a pure state, one always finds that $c_{\rm eff}(t\to\infty)\geq c_{\rm eff}(t=0)$, though $c_{\rm eff}(t)$, generally oscillates at finite times. We explore how this constraint may be related to RG flow irreversibility.

• Fluctuations of observables for free fermions in a harmonic trap at finite temperature

  Aurélien Grabsch, Satya N. Majumdar, Grégory Schehr, Christophe Texier

  SciPost Phys. 4, 014 (2018) · published 24 March 2018 | Toggle abstract · pdf

  We study a system of 1D noninteracting spinless fermions in a confining trap at finite temperature. We first derive a useful and general relation for the fluctuations of the occupation numbers valid for arbitrary confining trap, as well as for both canonical and grand canonical ensembles. Using this relation, we obtain compact expressions, in the case of the harmonic trap, for the variance of certain observables of the form of sums of a function of the fermions' positions, $\mathcal{L}=\sum_n h(x_n)$. Such observables are also called linear statistics of the positions. As anticipated, we demonstrate explicitly that these fluctuations do depend on the ensemble in the thermodynamic limit, as opposed to averaged quantities, which are ensemble independent. We have applied our general formalism to compute the fluctuations of the number of fermions $\mathcal{N}_+$ on the positive axis at finite temperature. Our analytical results are compared to numerical simulations. We discuss the universality of the results with respect to the nature of the confinement.

  11 citations

• (Non) equilibrium dynamics: a (broken) symmetry of the Keldysh generating functional

  Camille Aron, Giulio Biroli, Leticia F. Cugliandolo

  SciPost Phys. 4, 008 (2018) · published 31 January 2018 | Toggle abstract · pdf

  We unveil the universal (model-independent) symmetry satisfied by Schwinger-Keldysh quantum field theories whenever they describe equilibrium dynamics. This is made possible by a generalization of the Schwinger-Keldysh path-integral formalism in which the physical time can be re-parametrized to arbitrary contours in the complex plane. Strong relations between correlation functions, such as the fluctuation-dissipation theorems, are derived as immediate consequences of this symmetry of equilibrium. In this view, quantum non-equilibrium dynamics -- e.g. when driving with a time-dependent potential -- are seen as symmetry-breaking processes. The symmetry-breaking terms of the action are identified as a measure of irreversibility, or entropy creation, defined at the level of a single quantum trajectory. Moreover, they are shown to obey quantum fluctuation theorems. These results extend stochastic thermodynamics to the quantum realm.

  23 citations

• Drude Weight for the Lieb-Liniger Bose Gas

  Benjamin Doyon, Herbert Spohn

  SciPost Phys. 3, 039 (2017) · published 15 December 2017 | Toggle abstract · pdf

  Based on the method of hydrodynamic projections we derive a concise formula for the Drude weight of the repulsive Lieb-Liniger $\delta$-Bose gas. Our formula contains only quantities which are obtainable from the thermodynamic Bethe ansatz. The Drude weight is an infinite-dimensional matrix, or bilinear functional: it is bilinear in the currents, and each current may refer to a general linear combination of the conserved charges of the model. As a by-product we obtain the dynamical two-point correlation functions involving charge and current densities at small wavelengths and long times, and in addition the scaled covariance matrix of charge transfer. We expect that our formulas extend to other integrable quantum models.

  138 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.

  69 citations