Publications

Order by: publication date number of citations

• Thermal pure matrix product state in two dimensions: Tracking thermal equilibrium from paramagnet down to the Kitaev honeycomb spin liquid state

  Matthias Gohlke, Atsushi Iwaki, Chisa Hotta

  SciPost Phys. 15, 206 (2023) · published 24 November 2023 | Toggle abstract · pdf

  We present the first successful application of the matrix product state (MPS) representing a thermal quantum pure state (TPQ) in equilibrium in two spatial dimensions over almost the entire temperature range. We use the Kitaev honeycomb model as a prominent example hosting a quantum spin liquid (QSL) ground state to target the two specific-heat peaks previously solved nearly exactly using the free Majorana fermionic description. Starting from the high-temperature random state, our TPQ-MPS framework on a cylinder precisely reproduces these peaks, showing that the quantum many-body description based on spins can still capture the emergent itinerant Majorana fermions in a ${\mathbb Z}_2$ gauge field. The truncation process efficiently discards the high-energy states, eventually reaching the long-range entangled topological state approaching the exact ground state for a given finite size cluster. An advantage of TPQ-MPS over exact diagonalization or purification-based methods is its lowered numerical cost coming from a reduced effective Hilbert space even at finite temperature.

  1 citation

• Topological interpretation of color exchange invariants: hexagonal lattice on a torus

  Olivier Cépas, Peter M. Akhmetiev

  SciPost Phys. 10, 042 (2021) · published 18 February 2021 | Toggle abstract · pdf

  We explain a correspondence between some invariants in the dynamics of color exchange in a 2d coloring problem, which are polynomials of winding numbers, and linking numbers in 3d. One invariant is visualized as linking of lines on a special surface with Arf-Kervaire invariant one, and is interpreted as resulting from an obstruction to transform the surface into its chiral image with special continuous deformations. We also consider additional constraints on the dynamics and see how the surface is modified.

  1 citation

• Improved Hilbert space exploration algorithms for finite temperature calculations

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

  SciPost Phys. Core 6, 039 (2023) · published 23 May 2023 | Toggle abstract · pdf

  Computing correlation functions in strongly-interacting quantum systems is one of the most important challenges of modern condensed matter theory, due to their importance in the description of many physical observables. Simultaneously, this challenge is one of the most difficult to address, due to the inapplicability of traditional perturbative methods or the few-body limitations of numerical approaches. For special cases, where the model is integrable, methods based on the Bethe Ansatz have succeeded in computing the spectrum and given us analytical expressions for the matrix elements of physically important operators. However, leveraging these results to compute correlation functions generally requires the numerical evaluation of summations over eigenstates. To perform these summations efficiently, Hilbert space exploration algorithms have been developed which has resulted most notably in the ABACUS library. While this performs quite well for correlations on ground states or low-entropy states, the case of high entropy states (most importantly at finite temperatures or after a quantum quench) is more difficult, and leaves room for improvement. In this work, we develop a new Hilbert space exploration algorithm for the Lieb-Liniger model, specially tailored to optimize the computational order on finite-entropy states for correlations of density-related operators.

  1 citation

• Probabilistic deconstruction of a theory of gravity, part I: Flat space

  Sunok Josephine Suh

  SciPost Phys. 15, 174 (2023) · published 19 October 2023 | Toggle abstract · pdf

  We define and analyze a stochastic process in anti-de Sitter Jackiw-Teitelboim gravity, induced by the quantum dynamics of the boundary and whose random variable takes values in $AdS_2$. With the boundary in a thermal state and for appropriate parameters, we take the asymptotic limit of the quantum process at short time scales and flat space, and show associated classical joint distributions have the Markov property. We find that Einstein's equations of the theory, sans the cosmological constant term, arise in the semi-classical limit of the quantum evolution of probability under the asymptotic process. In particular, in flat Jackiw-Teitelboim gravity, the area of compactified space solved for by Einstein's equations can be identified as a probability density evolving under the Markovian process.

  1 citation

• Effective description of non-equilibrium currents in cold magnetized plasma

  Nabil Iqbal

  SciPost Phys. 12, 078 (2022) · published 28 February 2022 | Toggle abstract · pdf

  The dynamics of cold strongly magnetized plasma -- traditionally the domain of force-free electrodynamics -- has recently been reformulated in terms of symmetries and effective field theory, where the degrees of freedom are the momentum and magnetic flux carried by a fluid of cold strings. In physical applications where the electron mass can be neglected one might expect the presence of extra light charged modes -- electrons in the lowest Landau level -- propagating parallel to the magnetic field lines. We construct an effective description of such electric charges, describing their interaction with plasma degrees of freedom in terms of a new collective mode that can be thought of as a bosonization of the electric charge density along each field line. In this framework QED phenomena such as charged pair production and the axial anomaly are described at the classical level. Formally, our construction corresponds to gauging a particular part of the higher form symmetry associated with magnetic flux conservation. We study some simple applications of our effective theory, showing that the scattering of magnetosonic modes generically creates particles and that the rotating Michel monopole is now surrounded by a cloud of electric charge.

  1 citation

• Codebase release r1.1 for SmoQyDEAC.jl

  James Neuhaus, Nathan S. Nichols, Debshikha Banerjee, Benjamin Cohen-Stead, Thomas A. Maier, Adrian Del Maestro, Steven Johnston

  SciPost Phys. Codebases 39-r1.1 (2024) · published 12 November 2024 | Toggle abstract · src

  We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.

  1 citation

• Logarithmic correlation functions for critical dense polymers on the cylinder

  Alexi Morin-Duchesne, Jesper Lykke Jacobsen

  SciPost Phys. 7, 040 (2019) · published 30 September 2019 | Toggle abstract · pdf

  We compute lattice correlation functions for the model of critical dense polymers on a semi-infinite cylinder of perimeter $n$. In the lattice loop model, contractible loops have a vanishing fugacity whereas non-contractible loops have a fugacity $\alpha \in (0,\infty)$. These correlators are defined as ratios $Z(x)/Z_0$ of partition functions, where $Z_0$ is a reference partition function wherein only simple half-arcs are attached to the boundary of the cylinder. For $Z(x)$, the boundary of the cylinder is also decorated with simple half-arcs, but it also has two special positions $1$ and $x$ where the boundary condition is different. We investigate two such kinds of boundary conditions: (i) there is a single node at each of these points where a long arc is attached, and (ii) there are pairs of adjacent nodes at these points where two long arcs are attached. We find explicit expressions for these correlators for finite $n$ using the representation of the enlarged periodic Temperley-Lieb algebra in the XX spin chain. The resulting asymptotics as $n\to \infty$ are expressed as simple integrals that depend on the scaling parameter $\tau = \frac {x-1} n \in (0,1)$. For small $\tau$, the leading behaviours are proportional to $\tau^{1/4}$, $\tau^{1/4}\log \tau$, $\log \tau$ and $\log^2 \tau$. We interpret the lattice results in terms of ratios of conformal correlation functions. We assume that the corresponding boundary changing fields are highest weight states in irreducible, Kac or staggered Virasoro modules, with central charge $c=-2$ and conformal dimensions $\Delta = -\frac18$ or $\Delta = 0$. With these assumptions, we obtain differential equations of order two and three satisfied by the conformal correlation functions, solve these equations in terms of hypergeometric functions, and find a perfect agreement with the lattice results. We use the lattice results to compute structure constants and ratios thereof which appear in the operator product expansions of the boundary condition changing fields. The fusion of these fields is found to be non-abelian.

  1 citation

• Nonequilibrium quasiparticle distribution in superconducting resonators: Effect of pair-breaking photons

  Paul B. Fischer, Gianluigi Catelani

  SciPost Phys. 17, 070 (2024) · published 2 September 2024 | Toggle abstract · pdf

  Many superconducting devices rely on the finite gap in the excitation spectrum of a superconductor: thanks to this gap, at temperatures much smaller than the critical one the number of excitations (quasiparticles) that can impact the device's behavior is exponentially small. Nevertheless, experiments at low temperature usually find a finite, non-negligible density of quasiparticles whose origin has been attributed to various non-equilibrium phenomena. Here, we investigate the role of photons with energy exceeding the pair-breaking threshold $2\Delta$ as a possible source for these quasiparticles in superconducting resonators. Modeling the interacting system of quasiparticles, phonons, sub-gap and pair-breaking photons using a kinetic equation approach, we find analytical expressions for the quasiparticles' density and their energy distribution. Applying our theory to measurements of quality factor as function of temperature and for various readout powers, we find they could be explained by assuming a small number of photons above the pair-breaking threshold. We also show that frequency shift data can give evidence of quasiparticle heating.

  1 citation

• Development of the detectors for the DeeMe experiment

  Natsuki Teshima on behalf of the DeeMe Collaboration

  SciPost Phys. Proc. 1, 051 (2019) · published 22 February 2019 | Toggle abstract · pdf

  In the DeeMe experiment, approximately $70\ \mathrm{GHz/mm^{2}}$ prompt-charged particles will hit the multi-wire proportional chambers (MWPCs) before signal electrons come. To avoid the space charge effect, we developed fast HV-switching MWPCs in order to control the gas gain dynamically. All four of the MWPCs were manufactured. Last year, the circuit of readout amplifiers is slightly modified to further improve the efficiency of the detector, and we also started to investigate other possible choices of the gas mixture. In this article, the development of the detectors and results of performance tests will be presented.

  1 citation

• Renormalization-group improvement in hadronic $τ$ decays in 2018

  Diogo Boito, Pere Masjuan, Fabio Oliani

  SciPost Phys. Proc. 1, 049 (2019) · published 22 February 2019 | Toggle abstract · pdf

  One of the main sources of theoretical uncertainty in the extraction of the strong coupling from hadronic tau decays stems from the renormalization group improvement of the series. Perturbative series in QCD are divergent but are (most likely) asymptotic expansions. One needs knowledge about higher orders to be able to choose the optimal renormalization-scale setting procedure. Here, we discuss the use of Pad\'e approximants as a model-independent and robust method to extract information about the higher-order terms. We show that in hadronic \boldmath $\tau$ decays the fixed-order expansion, known as fixed-order perturbation theory (FOPT), is the most reliable mainstream method to set the scale. This fully corroborates previous conclusions based on the available knowledge about the leading renormalon singularities of the perturbative series.

  1 citation