Publications

Order by: publication date number of citations

• A new derivation of the relationship between diffusion coefficient and entropy in classical Brownian motion by the ensemble method

  Yi Liao, Xiao-Bo Gong

  SciPost Phys. Core 4, 015 (2021) · published 2 June 2021 | Toggle abstract · pdf

  The diffusion coefficient--a measure of dissipation, and the entropy--a measure of fluctuation are found to be intimately correlated in many physical systems. Unlike the fluctuation dissipation theorem in linear response theory, the correlation is often strongly non-linear. To understand this complex dependence, we consider the classical Brownian diffusion in this work. Under certain rational assumption, i.e. in the bi-component fluid mixture, the mass of the Brownian particle $M$ is far greater than that of the bath molecule $m$, we can adopt the weakly couple limit. Only considering the first-order approximation of the mass ratio $m/M$, we obtain a linear motion equation in the reference frame of the observer as a Brownian particle. Based on this equivalent equation, we get the Hamiltonian at equilibrium. Finally, using canonical ensemble method, we define a new entropy that is similar to the Kolmogorov-Sinai entropy. Further, we present an analytic expression of the relationship between the diffusion coefficient $D$ and the entropy $S$ in the thermal equilibrium, that is to say, $D =\frac{\hbar}{eM} \exp{[S/(k_Bd)]}$, where $d$ is the dimension of the space, $k_B$ the Boltzmann constant, $\hbar $ the reduced Planck constant and $e$ the Euler number. This kind of scaling relation has been well-known and well-tested since the similar one for single component is firstly derived by Rosenfeld with the expansion of volume ratio.

  1 citation

• Self-congruent point in critical matrix product states: An effective field theory for finite-entanglement scaling

  Jan T. Schneider, Atsushi Ueda, Yifan Liu, Andreas M. Läuchli, Masaki Oshikawa, Luca Tagliacozzo

  SciPost Phys. 18, 142 (2025) · published 29 April 2025 | Toggle abstract · pdf

  We set up an effective field theory formulation for the renormalization flow of matrix product states (MPS) with finite bond dimension, focusing on systems exhibiting finite-entanglement scaling close to a conformally invariant critical fixed point. We show that the finite MPS bond dimension $\chi$ is equivalent to introducing a perturbation by a relevant operator to the fixed-point Hamiltonian. The fingerprint of this mechanism is encoded in the $\chi$-independent universal transfer matrix's gap ratios, which are distinct from those predicted by the unperturbed Conformal Field Theory (CFT). This phenomenon defines a renormalization group self-congruent point, where the relevant coupling constant ceases to flow due to a balance of two effects; When increasing $\chi$, the infrared scale, set by the correlation length $\xi(\chi)$, increases, while the strength of the perturbation at the lattice scale decreases. The presence of a self-congruent point does not alter the validity of the finite-entanglement scaling hypothesis, since the self-congruent point is located at a finite distance from the critical fixed point, well inside the scaling regime of the CFT. We corroborate this framework with numerical evidence from the exact solution of the Ising model and density matrix renormalization group (DMRG) simulations of an effective lattice model.

  1 citation

• Spectral properties, localization transition and multifractal eigenvectors of the Laplacian on heterogeneous networks

  Jeferson D. da Silva, Diego Tapias, Peter Sollich, Fernando L. Metz

  SciPost Phys. 18, 047 (2025) · published 7 February 2025 | Toggle abstract · pdf

  We study the spectral properties and eigenvector statistics of the Laplacian on highly-connected networks with random coupling strengths and a gamma distribution of rescaled degrees. The spectral density, the distribution of the local density of states, the singularity spectrum and the multifractal exponents of this model exhibit a rich behaviour as a function of the first two moments of the coupling strengths and the variance of the rescaled degrees. In the case of random coupling strengths, the spectral density diverges within the bulk of the spectrum when degree fluctuations are strong enough. The emergence of this singular behaviour marks a transition from non-ergodic delocalized states to localized eigenvectors that exhibit pronounced multifractal scaling. For constant coupling strengths, the bulk of the spectrum is characterized by a regular spectral density. In this case, the corresponding eigenvectors display localization properties reminiscent of the critical point of the Anderson localization transition on random graphs.

  1 citation

• Codebase release 1.1 for TASEPy

  Luca Ciandrini, Richmond L. Crisostomo, Juraj Szavits-Nossan

  SciPost Phys. Codebases 22-r1.1 (2023) · published 13 December 2023 | Toggle abstract · src

  The totally asymmetric simple exclusion process (TASEP) is a paradigmatic lattice model for one-dimensional particle transport subject to excluded-volume interactions. Solving the inhomogeneous TASEP in which particles' hopping rate vary across the lattice is a long-standing problem. In recent years, a power series approximation (PSA) has been developed to tackle this problem, however no computer algorithm currently exists that implements this approximation. This paper addresses this issue by providing a Python-based package TASEPy that finds the steady state solution of the inhomogeneous TASEP for any set of hopping rates using the PSA truncated at a user-defined order.

  1 citation

• Lightcone Hamiltonian for Ising field theory I: $T < T_c$

  Andrew Liam Fitzpatrick, Emanuel Katz, Yuan Xin

  SciPost Phys. 18, 179 (2025) · published 5 June 2025 | Toggle abstract · pdf

  We study 2d Ising Field Theory (IFT) in the low-temperature phase in lightcone quantization, and show that integrating out zero modes generates a very compact form for the effective lightcone interaction that depends on the finite volume vacuum expectation value of the $\sigma$ operator. This form is most naturally understood in a conformal basis for the lightcone Hilbert space. We further verify that this simple form reproduces to high accuracy results for the spectra, the $c$-function, and the form-factors from integrability methods for the magnetic deformation of IFT. For generic non-integrable values of parameters we also compute the above observables and compare our numeric results to those of equal-time truncation. In particular, we report on new measurements of various bound-state form-factors as well as the stress-tensor spectral density. We find that the stress tensor spectral density provides additional evidence that certain resonances of IFT are surprisingly narrow, even at generic strong coupling. Explicit example code for constructing the effective Hamiltonian is included in an appendix.

  1 citation

• Exact first-order effect of interactions on the ground-state energy of harmonically-confined fermions

  Pierre Le Doussal, Naftali R. Smith, Nathan Argaman

  SciPost Phys. 17, 038 (2024) · published 7 August 2024 | Toggle abstract · pdf

  We consider a system of $N$ spinless fermions, interacting with each other via a power-law interaction $\epsilon/r^n$, and trapped in an external harmonic potential $V(r) = r^2/2$, in $d=1,2,3$ dimensions. For any $0 < n < d+2$, we obtain the ground-state energy $E_N$ of the system perturbatively in $\epsilon$, $E_{N}=E_{N}^{(0)}+\epsilon E_{N}^{(1)}+O(\epsilon^{2})$. We calculate $E_{N}^{(1)}$ exactly, assuming that $N$ is such that the "outer shell" is filled. For the case of $n=1$ (corresponding to a Coulomb interaction for $d=3$), we extract the $N \gg 1$ behavior of $E_{N}^{(1)}$, focusing on the corrections to the exchange term with respect to the leading-order term that is predicted from the local density approximation applied to the Thomas-Fermi approximate density distribution. The leading correction contains a logarithmic divergence, and is of particular importance in the context of density functional theory. We also study the effect of the interactions on the fermions' spatial density. Finally, we find that our result for $E_{N}^{(1)}$ significantly simplifies in the case where $n$ is even.

  1 citation

• Unifying description of the damping regimes of a stochastic particle in a periodic potential

  Antonio Piscitelli, Massimo Pica Ciamarra

  SciPost Phys. 3, 018 (2017) · published 4 September 2017 | Toggle abstract · pdf

  We analyze the classical problem of the stochastic dynamics of a particle confined in a periodic potential, through the so called Il'in and Khasminskii model, with a novel semi-analytical approach. Our approach gives access to the transient and the asymptotic dynamics in all damping regimes, which are difficult to investigate in the usual Brownian model. We show that the crossover from the overdamped to the underdamped regime is associated with the loss of a typical time scale and of a typical length scale, as signaled by the divergence of the probability distribution of a certain dynamical event. In the underdamped regime, normal diffusion coexists with a non Gaussian displacement probability distribution for a long transient, as recently observed in a variety of different systems. We rationalize the microscopic physical processes leading to the non-Gaussian behavior, as well as the timescale to recover the Gaussian statistics. The theoretical results are supported by numerical calculations and are compared to those obtained for the Brownian model.

  1 citation

• A minimal tensor network beyond free fermions

  Carolin Wille, Maksimilian Usoltcev, Jens Eisert, Alexander Altland

  SciPost Phys. 18, 196 (2025) · published 19 June 2025 | Toggle abstract · pdf

  This work proposes a minimal model extending the duality between classical statistical spin systems and fermionic systems beyond the case of free fermions. A Jordan-Wigner transformation applied to a two-dimensional tensor network maps the partition sum of a classical statistical mechanics model to a Grassmann variable integral, structurally similar to the path integral for interacting fermions in two dimensions. The resulting model is simple, featuring only two parameters: one governing spin-spin interaction (dual to effective hopping strength in the fermionic picture), the other measuring the deviation from the free fermion limit. Nevertheless, it exhibits a rich phase diagram, partially stabilized by elements of topology, and featuring three phases meeting at a multicritical point. Besides the interpretation as a spin and fermionic system, the model is closely related to loop gas and vertex models and can be interpreted as a parity-preserving (non-unitary) circuit. Its minimal construction makes it an ideal reference system for studying non-linearities in tensor networks and deriving results by means of duality.

  1 citation

• 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

• Codebase release 0.2 for JFSD

  Kim William Torre, Raoul D. Schram, Joost de Graaf

  SciPost Phys. Codebases 56-r0.2 (2025) · published 30 June 2025 | Toggle abstract · src

  Stokesian Dynamics (SD) is a powerful computational framework for simulating the motion of particles in a viscous Newtonian fluid under Stokes-flow conditions. Traditional SD implementations can be computationally expensive as they rely on the inversion of large mobility matrices to determine hydrodynamic interactions. Recently, however, the simulation of thermalized systems with large numbers of particles has become feasible [Fiore and Swan, J. Fluid Mech. 878, 544 (2019)]. Their "fast Stokesian dynamics" (FSD) method leverages a saddle-point formulation to ensure overall scaling of the algorithm that is linear in the number of particles $\mathcal{O}(N)$; performance relies on dedicated graphics-processing-unit computing. Here, we present a different route toward implementing FSD, which instead leverages the Just-in-Time (JIT) compilation capabilities of Google JAX. We refer to this implementation as JFSD and perform benchmarks on it to verify that it has the right scaling and is sufficiently fast by the standards of modern computational physics. In addition, we provide a series of physical test cases that help ensure accuracy and robustness, as the code undergoes further development. Thus, JFSD is ready to facilitate the study of hydrodynamic effects in particle suspensions across the domains of soft, active, and granular matter.

  1 citation