Publications

Order by: publication date number of citations

• Fast counting with tensor networks

  Stefanos Kourtis, Claudio Chamon, Eduardo R. Mucciolo, Andrei E. Ruckenstein

  SciPost Phys. 7, 060 (2019) · published 12 November 2019 | Toggle abstract · pdf

  We introduce tensor network contraction algorithms for counting satisfying assignments of constraint satisfaction problems (#CSPs). We represent each arbitrary #CSP formula as a tensor network, whose full contraction yields the number of satisfying assignments of that formula, and use graph theoretical methods to determine favorable orders of contraction. We employ our heuristics for the solution of #P-hard counting boolean satisfiability (#SAT) problems, namely monotone #1-in-3SAT and #Cubic-Vertex-Cover, and find that they outperform state-of-the-art solvers by a significant margin.

  24 citations

• Navigator function for the conformal bootstrap

  Marten Reehorst, Slava Rychkov, David Simmons-Duffin, Benoit Sirois, Ning Su, Balt van Rees

  SciPost Phys. 11, 072 (2021) · published 28 September 2021 | Toggle abstract · pdf

  Current numerical conformal bootstrap techniques carve out islands in theory space by repeatedly checking whether points are allowed or excluded. We propose a new method for searching theory space that replaces the binary information “allowed”/“excluded” with a continuous “navigator” function that is negative in the allowed region and positive in the excluded region. Such a navigator function allows one to efficiently explore high-dimensional parameter spaces and smoothly sail towards any islands they may contain. The specific functions we introduce have several attractive features: they are well-defined in large regions of parameter space, can be computed with standard methods, and evaluation of their gradient is immediate due to an SDP gradient formula that we provide. The latter property allows for the use of efficient quasi-Newton optimization methods, which we illustrate by navigating towards the 3d Ising island.

  24 citations

• Euler-scale dynamical correlations in integrable systems with fluid motion

  Frederik S. Møller, Gabriele Perfetto, Benjamin Doyon, Jörg Schmiedmayer

  SciPost Phys. Core 3, 016 (2020) · published 22 December 2020 | Toggle abstract · pdf

  We devise an iterative scheme for numerically calculating dynamical two-point correlation functions in integrable many-body systems, in the Eulerian scaling limit. Expressions for these were originally derived in Ref. [1] by combining the fluctuation-dissipation principle with generalized hydrodynamics. Crucially, the scheme is able to address non-stationary, inhomogeneous situations, when motion occurs at the Euler-scale of hydrodynamics. In such situations, in interacting systems, the simple correlations due to fluid modes propagating with the flow receive subtle corrections, which we test. Using our scheme, we study the spreading of correlations in several integrable models from inhomogeneous initial states. For the classical hard-rod model we compare our results with Monte-Carlo simulations and observe excellent agreement at long time-scales, thus providing the first demonstration of validity for the expressions derived in Ref. [1]. We also observe the onset of the Euler-scale limit for the dynamical correlations.

  23 citations

• Bosonic and fermionic Gaussian states from Kähler structures

  Lucas Hackl, Eugenio Bianchi

  SciPost Phys. Core 4, 025 (2021) · published 22 September 2021 | Toggle abstract · pdf

  We show that bosonic and fermionic Gaussian states (also known as "squeezed coherent states") can be uniquely characterized by their linear complex structure $J$ which is a linear map on the classical phase space. This extends conventional Gaussian methods based on covariance matrices and provides a unified framework to treat bosons and fermions simultaneously. Pure Gaussian states can be identified with the triple $(G,\Omega,J)$ of compatible K\"ahler structures, consisting of a positive definite metric $G$, a symplectic form $\Omega$ and a linear complex structure $J$ with $J^2=-1\!\!1$. Mixed Gaussian states can also be identified with such a triple, but with $J^2\neq -1\!\!1$. We apply these methods to show how computations involving Gaussian states can be reduced to algebraic operations of these objects, leading to many known and some unknown identities. We apply these methods to the study of (A) entanglement and complexity, (B) dynamics of stable systems, (C) dynamics of driven systems. From this, we compile a comprehensive list of mathematical structures and formulas to compare bosonic and fermionic Gaussian states side-by-side.

  23 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

• Mean-field model of interacting quasilocalized excitations in glasses

  Corrado Rainone, Pierfrancesco Urbani, Francesco Zamponi, Edan Lerner, and Eran Bouchbinder

  SciPost Phys. Core 4, 008 (2021) · published 16 April 2021 | Toggle abstract · pdf

  Structural glasses feature quasilocalized excitations whose frequencies $\omega$ follow a universal density of states ${\cal D}(\omega)\!\sim\!\omega^4$. Yet, the underlying physics behind this universality is not yet fully understood. Here we study a mean-field model of quasilocalized excitations in glasses, viewed as groups of particles embedded inside an elastic medium and described collectively as anharmonic oscillators. The oscillators, whose harmonic stiffness is taken from a rather featureless probability distribution (of upper cutoff $\kappa_0$) in the absence of interactions, interact among themselves through random couplings (characterized by strength $J$) and with the surrounding elastic medium (an interaction characterized by a constant force $h$). We first show that the model gives rise to a gapless density of states ${\cal D}(\omega)\!=\!A_{\rm g}\,\omega^4$ for a broad range of model parameters, expressed in terms of the strength of stabilizing anharmonicity, which plays a decisive role in the model. Then --- using scaling theory and numerical simulations --- we provide a complete understanding of the non-universal prefactor $A_{\rm g}(h,J,\kappa_0)$, of the oscillators' interaction-induced mean square displacement and of an emerging characteristic frequency, all in terms of properly identified dimensionless quantities. In particular, we show that $A_{\rm g}(h,J,\kappa_0)$ is a non-monotonic function of $J$ for a fixed $h$, varying predominantly exponentially with $-(\kappa_0 h^{2/3}\!/J^2)$ in the weak interactions (small $J$) regime --- reminiscent of recent observations in computer glasses --- and predominantly decays as a power-law for larger $J$, in a regime where $h$ plays no role. We discuss the physical interpretation of the model and its possible relations to available observations in structural glasses, along with delineating some future research directions.

  22 citations

• Large fluctuations of the KPZ equation in a half-space

  Alexandre Krajenbrink, Pierre Le Doussal

  SciPost Phys. 5, 032 (2018) · published 12 October 2018 | Toggle abstract · pdf

  We investigate the short-time regime of the KPZ equation in $1+1$ dimensions and develop a unifying method to obtain the height distribution in this regime, valid whenever an exact solution exists in the form of a Fredholm Pfaffian or determinant. These include the droplet and stationary initial conditions in full space, previously obtained by a different method. The novel results concern the droplet initial condition in a half space for several Neumann boundary conditions: hard wall, symmetric, and critical. In all cases, the height probability distribution takes the large deviation form $P(H,t) \sim \exp( - \Phi(H)/\sqrt{t})$ for small time. We obtain the rate function $\Phi(H)$ analytically for the above cases. It has a Gaussian form in the center with asymmetric tails, $|H|^{5/2}$ on the negative side, and $H^{3/2}$ on the positive side. The amplitude of the left tail for the half-space is found to be half the one of the full space. As in the full space case, we find that these left tails remain valid at all times. In addition, we present here (i) a new Fredholm Pfaffian formula for the solution of the hard wall boundary condition and (ii) two Fredholm determinant representations for the solutions of the hard wall and the symmetric boundary respectively.

  22 citations

• A classical density functional from machine learning and a convolutional neural network

  Shang-Chun Lin, Martin Oettel

  SciPost Phys. 6, 025 (2019) · published 26 February 2019 | Toggle abstract · pdf

  We use machine learning methods to approximate a classical density functional. As a study case, we choose the model problem of a Lennard Jones fluid in one dimension where there is no exact solution available and training data sets must be obtained from simulations. After separating the excess free energy functional into a "repulsive" and an "attractive" part, machine learning finds a functional in weighted density form for the attractive part. The density profile at a hard wall shows good agreement for thermodynamic conditions beyond the training set conditions. This also holds for the equation of state if it is evaluated near the training temperature. We discuss the applicability to problems in higher dimensions.

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

  20 citations

• Exact out-of-equilibrium steady states in the semiclassical limit of the interacting Bose gas

  Giuseppe Del Vecchio Del Vecchio, Alvise Bastianello, Andrea De Luca, Giuseppe Mussardo

  SciPost Phys. 9, 002 (2020) · published 6 July 2020 | Toggle abstract · pdf

  We study the out-of-equilibrium properties of a classical integrable non-relativistic theory, with a time evolution initially prepared with a finite energy density in the thermodynamic limit. The theory considered here is the Non-Linear Schrodinger equation which describes the dynamics of the one-dimensional interacting Bose gas in the regime of high occupation numbers. The main emphasis is on the determination of the late-time Generalised Gibbs Ensemble (GGE), which can be efficiently semi-numerically computed on arbitrary initial states, completely solving the famous quench problem in the classical regime. We take advantage of known results in the quantum model and the semiclassical limit to achieve new exact results for the momenta of the density operator on arbitrary GGEs, which we successfully compare with ab-initio numerical simulations. Furthermore, we determine the whole probability distribution of the density operator (full counting statistics), whose exact expression is still out of reach in the quantum model.

  18 citations