Publications

Order by: publication date number of citations

• Full counting statistics for interacting trapped fermions

  Naftali R. Smith, Pierre Le Doussal, Satya N. Majumdar, Grégory Schehr

  SciPost Phys. 11, 110 (2021) · published 22 December 2021 | Toggle abstract · pdf

  We study $N$ spinless fermions in their ground state confined by an external potential in one dimension with long range interactions of the general Calogero-Sutherland type. For some choices of the potential this system maps to standard random matrix ensembles for general values of the Dyson index $\beta$. In the fermion model $\beta$ controls the strength of the interaction, $\beta=2$ corresponding to the noninteracting case. We study the quantum fluctuations of the number of fermions ${\cal N}_{\cal D}$ in a domain $\cal{D}$ of macroscopic size in the bulk of the Fermi gas. We predict that for general $\beta$ the variance of ${\cal N}_{\cal D}$ grows as $A_{\beta} \log N + B_{\beta}$ for $N \gg 1$ and we obtain a formula for $A_\beta$ and $B_\beta$. This is based on an explicit calculation for $\beta\in\left\{ 1,2,4\right\} $ and on a conjecture that we formulate for general $\beta$. This conjecture further allows us to obtain a universal formula for the higher cumulants of ${\cal N}_{\cal D}$. Our results for the variance in the microscopic regime are found to be consistent with the predictions of the Luttinger liquid theory with parameter $K = 2/\beta$, and allow to go beyond. In addition we present families of interacting fermion models in one dimension which, in their ground states, can be mapped onto random matrix models. We obtain the mean fermion density for these models for general interaction parameter $\beta$. In some cases the fermion density exhibits interesting transitions, for example we obtain a noninteracting fermion formulation of the Gross-Witten-Wadia model.

  17 citations

• Engineering spectral properties of non-interacting lattice Hamiltonians

  Ali G. Moghaddam, Dmitry Chernyavsky, Corentin Morice, Jasper van Wezel, Jeroen van den Brink

  SciPost Phys. 11, 109 (2021) · published 22 December 2021 | Toggle abstract · pdf

  We investigate the spectral properties of one-dimensional lattices with position-dependent hopping amplitudes and on-site potentials that are smooth bounded functions of the position. We find an exact integral form for the density of states (DOS) in the limit of an infinite number of sites, which we derive using a mixed Bloch-Wannier basis consisting of piecewise Wannier functions. Next, we provide an exact solution for the inverse problem of constructing the position-dependence of hopping in a lattice model yielding a given DOS. We confirm analytic results by comparing them to numerics obtained by exact diagonalization for various incarnations of position-dependent hoppings and on-site potentials. Finally, we generalize the DOS integral form to multi-orbital tight-binding models with longer-range hoppings and in higher dimensions.

  4 citations

• Conserved momenta of ferromagnetic solitons through the prism of differential geometry

  Xingjian Di, Oleg Tchernyshyov

  SciPost Phys. 11, 108 (2021) · published 21 December 2021 | Toggle abstract · pdf

  The relation between symmetries and conservation laws for solitons in a ferromagnet is complicated by the presence of gyroscopic (precessional) forces, whose description in the Lagrangian framework involves a background gauge field. This makes canonical momenta gauge-dependent and requires a careful application of Noether's theorem. We show that Cartan's theory of differential forms is a natural language for this task. We use it to derive conserved momenta of the Belavin--Polyakov skyrmion, whose symmetries include translation, global spin rotation, and dilation.

  3 citations

• Entanglement dynamics in Rule 54: Exact results and quasiparticle picture

  Katja Klobas, Bruno Bertini

  SciPost Phys. 11, 107 (2021) · published 20 December 2021 | Toggle abstract · pdf

  We study the entanglement dynamics generated by quantum quenches in the quantum cellular automaton Rule 54. We consider the evolution from a recently introduced class of solvable initial states. States in this class relax (locally) to a one-parameter family of Gibbs states and the thermalisation dynamics of local observables can be characterised exactly by means of an evolution in space. Here we show that the latter approach also gives access to the entanglement dynamics and derive exact formulas describing the asymptotic linear growth of all Rényi entropies in the thermodynamic limit and their eventual saturation for finite subsystems. While in the case of von Neumann entropy we recover exactly the predictions of the quasiparticle picture, we find no physically meaningful quasiparticle description for other Rényi entropies. Our results apply to both homogeneous and inhomogeneous quenches.

  30 citations

• Exact relaxation to Gibbs and non-equilibrium steady states in the quantum cellular automaton Rule 54

  Katja Klobas, Bruno Bertini

  SciPost Phys. 11, 106 (2021) · published 20 December 2021 | Toggle abstract · pdf

  We study the out-of-equilibrium dynamics of the quantum cellular automaton Rule 54 using a time-channel approach. We exhibit a family of (non-equilibrium) product states for which we are able to describe exactly the full relaxation dynamics. We use this to prove that finite subsystems relax to a one-parameter family of Gibbs states. We also consider inhomogeneous quenches. Specifically, we show that when the two halves of the system are prepared in two different solvable states, finite subsystems at finite distance from the centre eventually relax to the non-equilibrium steady state (NESS) predicted by generalised hydrodynamics. To the best of our knowledge, this is the first exact description of the relaxation to a NESS in an interacting system and, therefore, the first independent confirmation of generalised hydrodynamics for an inhomogeneous quench.

  27 citations

• Generalized eigenproblem without fermion doubling for Dirac fermions on a lattice

  M. J. Pacholski, G. Lemut, J. Tworzydło, C. W. J. Beenakker

  SciPost Phys. 11, 105 (2021) · published 14 December 2021 | Toggle abstract · pdf

  The spatial discretization of the single-cone Dirac Hamiltonian on the surface of a topological insulator or superconductor needs a special "staggered" grid, to avoid the appearance of a spurious second cone in the Brillouin zone. We adapt the Stacey discretization from lattice gauge theory to produce a generalized eigenvalue problem, of the form ${\mathcal H}\psi=E {\mathcal P}\psi$, with Hermitian tight-binding operators ${\mathcal H}$, ${\mathcal P}$, a locally conserved particle current, and preserved chiral and symplectic symmetries. This permits the study of the spectral statistics of Dirac fermions in each of the four symmetry classes A, AII, AIII, and D.

  10 citations

• On computing non-equilibrium dynamics following a quench

  Neil J. Robinson, Albertus J. J. M. de Klerk, Jean-Sébastien Caux

  SciPost Phys. 11, 104 (2021) · published 9 December 2021 | Toggle abstract · pdf

  Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.

  6 citations

• Dual applications of Chebyshev polynomials method: Efficiently finding thousands of central eigenvalues for many-spin systems

  Haoyu Guan, Wenxian Zhang

  SciPost Phys. 11, 103 (2021) · published 8 December 2021 | Toggle abstract · pdf

  Computation of a large group of interior eigenvalues at the middle spectrum is an important problem for quantum many-body systems, where the level statistics provides characteristic signatures of quantum chaos. We propose an exact numerical method, dual applications of Chebyshev polynomials (DACP), to simultaneously find thousands of central eigenvalues, where the level space decreases exponentially with the system size. To disentangle the near-degenerate problem, we employ twice the Chebyshev polynomials, to construct an exponential semicircle filter as a preconditioning step and to generate a large set of proper basis states in the desired subspace. Numerical calculations on Ising spin chain and spin glass shards confirm the correctness and efficiency of DACP. As numerical results demonstrate, DACP is $30$ times faster than the state-of-the-art shift-invert method for the Ising spin chain while $8$ times faster for the spin glass shards. In contrast to the shift-invert method, the computation time of DACP is only weakly influenced by the required number of eigenvalues, which renders it a powerful tool for large scale eigenvalues computations. Moreover, the consumed memory also remains a small constant (5.6 GB) for spin-$1/2$ systems consisting of up to $20$ spins, making it desirable for parallel computing.

  1 citation

• Rotating 5D black holes: Interactions and deformations near extremality

  Alejandra Castro, Juan F. Pedraza, Chiara Toldo, Evita Verheijden

  SciPost Phys. 11, 102 (2021) · published 7 December 2021 | Toggle abstract · pdf

  We study a two-dimensional theory of gravity coupled to matter that is relevant to describe holographic properties of black holes with two equal angular momenta in five dimensions (with or without cosmological constant). We focus on the near-horizon geometry of the near-extremal black hole, where the effective theory reduces to Jackiw-Teitelboim (JT) gravity coupled to a massive scalar field. We compute the corrections to correlation functions due to cubic interactions present in this theory. A novel feature is that these corrections do not have a definite sign: for AdS$_5$ black holes the sign depends on the mass of the extremal solution. We discuss possible interpretations of these corrections from a gravitational and holographic perspective. We also quantify the imprint of the JT sector on the UV region, i.e. how these degrees of freedom, characteristic for the near-horizon region, influence the asymptotically far region of the black hole. This gives an interesting insight on how to interpret the IR modes in the context of their UV completion, which depends on the environment that contains the black hole.

  15 citations

• Emergent fractal phase in energy stratified random models

  Anton G. Kutlin, Ivan M. Khaymovich

  SciPost Phys. 11, 101 (2021) · published 6 December 2021 | Toggle abstract · pdf

  We study the effects of partial correlations in kinetic hopping terms of long-range disordered random matrix models on their localization properties. We consider a set of models interpolating between fully-localized Richardson's model and the celebrated Rosenzweig-Porter model (with implemented translation-invariant symmetry). In order to do this, we propose the energy-stratified spectral structure of the hopping term allowing one to decrease the range of correlations gradually. We show both analytically and numerically that any deviation from the completely correlated case leads to the emergent non-ergodic delocalization in the system unlike the predictions of localization of cooperative shielding. In order to describe the models with correlated kinetic terms, we develop the generalization of the Dyson Brownian motion and cavity approaches basing on stochastic matrix process with independent rank-one matrix increments and examine its applicability to the above set of models.

  18 citations