Publications

Order by: publication date number of citations

• Entanglement evolution after a global quench across a conformal defect

  Luca Capizzi, Viktor Eisler

  SciPost Phys. 14, 070 (2023) · published 13 April 2023 | Toggle abstract · pdf

  We study the evolution of entanglement after a global quench in a one-dimensional quantum system with a localized impurity. For systems described by a conformal field theory, the entanglement entropy between the two regions separated by the defect grows linearly in time. Introducing the notion of boundary twist fields, we show how the slope of this growth can be related to the effective central charge that emerges in the study of ground-state entropy in the presence of the defect. On the other hand, we also consider a particular lattice realization of the quench in a free-fermion chain with a conformal defect. Starting from a gapped initial state, we obtain the slope via a quasiparticle ansatz and observe small discrepancies between the field theory and lattice results, which persist even in the limit of a vanishing gap.

  7 citations

• Replica method for eigenvalues of real Wishart product matrices

  Jacob A. Zavatone-Veth, Cengiz Pehlevan

  SciPost Phys. Core 6, 026 (2023) · published 5 April 2023 | Toggle abstract · pdf

  We show how the replica method can be used to compute the asymptotic eigenvalue spectrum of a real Wishart product matrix. For unstructured factors, this provides a compact, elementary derivation of a polynomial condition on the Stieltjes transform first proved by Müller [IEEE Trans. Inf. Theory. 48, 2086-2091 (2002)]. We then show how this computation can be extended to ensembles where the factors are drawn from matrix Gaussian distributions with general correlation structure. For both unstructured and structured ensembles, we derive polynomial conditions on the average values of the minimum and maximum eigenvalues, which in the unstructured case match the results obtained by Akemann, Ipsen, and Kieburg [Phys. Rev. E 88, 052118 (2013)] for the complex Wishart product ensemble.

  2 citations

• Codebase release 0.2 for Pychastic

  Radost Waszkiewicz, Maciej Bartczak, Kamil Kolasa, Maciej Lisicki

  SciPost Phys. Codebases 11-r0.2 (2023) · published 3 April 2023 | Toggle abstract · src

  In the last decade, Python-powered physics simulations ecosystem has been growing steadily, allowing greater interoperability, and becoming an important tool in numerical exploration of physical phenomena, particularly in soft matter systems. Driven by the need for fast and precise numerical integration in colloidal dynamics, here we formulate the problem of Brownian Dynamics (BD) in a mathematically consistent formalism of the Itō calculus, and develop a Python package to assist numerical computations. We show that, thanks to the automatic differentiation packages, the classical truncated Taylor-Itō integrators can be implemented without the burden of computing the derivatives of the coefficient functions beforehand. Furthermore, we show how to circumvent the difficulties of BD simulations such as calculations of the divergence of the mobility tensor in the diffusion equation and discontinuous trajectories encountered when working with dynamics on $S^2$ and $SO(3)$. The resulting Python package, Pychastic, is capable of performing BD simulations including hydrodynamic interactions at speeds comparable to dedicated implementations in lower-level programming languages, but with a much simpler end-user interface.

  1 citation

• Pychastic: Precise Brownian dynamics using Taylor-Itō integrators in Python

  Radost Waszkiewicz, Maciej Bartczak, Kamil Kolasa, Maciej Lisicki

  SciPost Phys. Codebases 11 (2023) · published 3 April 2023 | Toggle abstract · pdf

  In the last decade, Python-powered physics simulations ecosystem has been growing steadily, allowing greater interoperability, and becoming an important tool in numerical exploration of physical phenomena, particularly in soft matter systems. Driven by the need for fast and precise numerical integration in colloidal dynamics, here we formulate the problem of Brownian Dynamics (BD) in a mathematically consistent formalism of the Itō calculus, and develop a Python package to assist numerical computations. We show that, thanks to the automatic differentiation packages, the classical truncated Taylor-Itō integrators can be implemented without the burden of computing the derivatives of the coefficient functions beforehand. Furthermore, we show how to circumvent the difficulties of BD simulations such as calculations of the divergence of the mobility tensor in the diffusion equation and discontinuous trajectories encountered when working with dynamics on $S^2$ and $SO(3)$. The resulting Python package, Pychastic, is capable of performing BD simulations including hydrodynamic interactions at speeds comparable to dedicated implementations in lower-level programming languages, but with a much simpler end-user interface.

  3 citations

• Macroscopic fluctuation theory and current fluctuations in active lattice gases

  Tal Agranov, Sunghan Ro, Yariv Kafri, Vivien Lecomte

  SciPost Phys. 14, 045 (2023) · published 22 March 2023 | Toggle abstract · pdf

  We study the current large deviations for a lattice model of interacting active particles displaying a motility-induced phase separation (MIPS). To do this, we first derive the exact fluctuating hydrodynamics of the model in the large system limit. On top of the usual Gaussian noise terms the theory also presents Poissonian noise terms, that we fully account for. We find a dynamical phase transition between flat density profiles and sharply phase-separated traveling waves, and we derive the associated phase diagram together with the large deviation function for all phases, including the one displaying MIPS. We show how the results can be obtained using methods similar to those of equilibrium phase separation, in spite of the nonequilibrium nature of the problem.

  20 citations

• Learning a restricted Boltzmann machine using biased Monte Carlo sampling

  Nicolas Béreux, Aurélien Decelle, Cyril Furtlehner, Beatriz Seoane

  SciPost Phys. 14, 032 (2023) · published 14 March 2023 | Toggle abstract · pdf

  Restricted Boltzmann Machines are simple and powerful generative models that can encode any complex dataset. Despite all their advantages, in practice the trainings are often unstable and it is difficult to assess their quality because the dynamics are affected by extremely slow time dependencies. This situation becomes critical when dealing with low-dimensional clustered datasets, where the time required to sample ergodically the trained models becomes computationally prohibitive. In this work, we show that this divergence of Monte Carlo mixing times is related to a phenomenon of phase coexistence, similar to that which occurs in physics near a first-order phase transition. We show that sampling the equilibrium distribution using the Markov chain Monte Carlo method can be dramatically accelerated when using biased sampling techniques, in particular the Tethered Monte Carlo (TMC) method. This sampling technique efficiently solves the problem of evaluating the quality of a given trained model and generating new samples in a reasonable amount of time. Moreover, we show that this sampling technique can also be used to improve the computation of the log-likelihood gradient during training, leading to dramatic improvements in training RBMs with artificial clustered datasets. On real low-dimensional datasets, this new training method fits RBM models with significantly faster relaxation dynamics than those obtained with standard PCD recipes. We also show that TMC sampling can be used to recover the free-energy profile of the RBM. This proves to be extremely useful to compute the probability distribution of a given model and to improve the generation of new decorrelated samples in slow PCD-trained models. The main limitations of this method are, first, the restriction to effective low-dimensional datasets and, second, the fact that the Tethered MC method breaks the possibility of performing parallel alternative Monte Carlo updates, which limits the size of the systems we can consider in practice.

  9 citations

• Multi species asymmetric simple exclusion process with impurity activated flips

  Amit Kumar Chatterjee, Hisao Hayakawa

  SciPost Phys. 14, 016 (2023) · published 8 February 2023 | Toggle abstract · pdf

  We obtain an exact matrix product steady state for a class of multi species asymmetric simple exclusion process with impurities, under periodic boundary condition. Alongside the usual hopping dynamics, an additional flip dynamics is activated only in the presence of impurities. Although the microscopic dynamics renders the system to be non-ergodic, exact analytical results for observables are obtained in steady states for a specific class of initial configurations. Interesting physical features including negative differential mobility and transition of correlations from negative to positive with changing vacancy density, have been observed. We discuss plausible connections of this exactly solvable model with multi lane asymmetric simple exclusion processes as well as enzymatic chemical reactions.

  3 citations

• Time-reparametrization invariances, multithermalization and the Parisi scheme

  Jorge Kurchan

  SciPost Phys. Core 6, 001 (2023) · published 17 January 2023 | Toggle abstract · pdf

  The Parisi scheme for equilibrium and the corresponding slow dynamics with multithermalization - same temperature common to all observables, different temperatures only possible at widely separated timescales – imply one another. Consistency requires that two systems brought into infinitesimal coupling be able to rearrange their timescales in order that all their temperatures match: this time reorganisation is only possible because the systems have a set of time-reparametrization invariances, that are thus seen to be an essential component of the scenario.

  3 citations

• Electrostatic solution of massless quenches in Luttinger liquids

  Paola Ruggiero, Pasquale Calabrese, Thierry Giamarchi, Laura Foini

  SciPost Phys. 13, 111 (2022) · published 18 November 2022 | Toggle abstract · pdf

  The study of non-equilibrium dynamics of many-body systems after a quantum quench received a considerable boost and a deep theoretical understanding from the path integral formulation in imaginary time. However, the celebrated problem of a quench in the Luttinger parameter of a one dimensional quantum critical system (massless quench) has so far only been solved in the real-time Heisenberg picture. In order to bridge this theoretical gap and to understand on the same ground massive and massless quenches, we study the problem of a gaussian field characterized by a coupling parameter $K$ within a strip and a different one $K_0$ in the remaining two semi-infinite planes. We give a fully analytical solution using the electrostatic analogy with the problem of a dielectric material within a strip surrounded by an infinite medium of different dielectric constant, and exploiting the method of charge images. After analytic continuation, this solution allows us to obtain all the correlation functions after the quench within a path integral approach in imaginary time, thus recovering and generalizing the results in real time. Furthermore, this imaginary-time approach establishes a remarkable connection between the quench and the famous problem of the conductivity of a Tomonaga-Luttinger liquid coupled to two semi-infinite leads: the two are in fact related by a rotation of the spacetime coordinates.

  6 citations

• Navigating through the O(N) archipelago

  Benoit Sirois

  SciPost Phys. 13, 081 (2022) · published 4 October 2022 | Toggle abstract · pdf

  A novel method for finding allowed regions in the space of CFT-data, coined navigator method, was recently proposed in arXiv:2104.09518. Its efficacy was demonstrated in the simplest example possible, i.e. that of the mixed-correlator study of the 3D Ising Model. In this paper, we would like to show that the navigator method may also be applied to the study of the family of $d$-dimensional $O(N)$ models. We will aim to follow these models in the $(d,N)$ plane. We will see that the "sailing" from island to island can be understood in the context of the navigator as a parametric optimization problem, and we will exploit this fact to implement a simple and effective path-following algorithm. By sailing with the navigator through the $(d,N)$ plane, we will provide estimates of the scaling dimensions $(\Delta_{\phi},\Delta_{s},\Delta_{t})$ in the entire range $(d,N) \in [3,4] \times [1,3]$. We will show that to our level of precision, we cannot see the non-unitary nature of the $O(N)$ models due to the fractional values of $d$ or $N$ in this range. We will also study the limit $N \xrightarrow[]{} 1$, and see that we cannot find any solution to the unitary mixed-correlator crossing equations below $N=1$.

  6 citations