Publications

Order by: publication date number of citations

• Probing black hole microstate evolution with networks and random walks

  Anthony M. Charles, Daniel R. Mayerson

  SciPost Phys. 8, 077 (2020) · published 14 May 2020 | Toggle abstract · pdf

  We model black hole microstates and quantum tunneling transitions between them with networks and simulate their time evolution using well-established tools in network theory. In particular, we consider two models based on Bena-Warner three-charge multi-centered microstates and one model based on the D1-D5 system; we use network theory methods to determine how many centers (or D1-D5 string strands) we expect to see in a typical late-time state. We find three distinct possible phases in parameter space for the late-time behaviour of these networks, which we call ergodic, trapped, and amplified, depending on the relative importance and connectedness of microstates. We analyze in detail how these different phases of late-time behavior are related to the underlying physics of the black hole microstates. Our results indicate that the expected properties of microstates at late times cannot always be determined simply by entropic arguments; typicality is instead a highly non-trivial, emergent property of the full Hilbert space of microstates.

  1 citation

• The $D_{03}(2380)$ dibaryon resonance excitation in the $pd\to pd\pi\pi$ reaction.

  Nurbek Tursunbayev, Yuriy Uzikov

  SciPost Phys. Proc. 3, 056 (2020) · published 27 February 2020 | Toggle abstract · pdf

  The recent ANKE@COSY data on the differential cross section of reaction $pd \to pd\pi\pi$ demonstrate a peak at invariant mass of the final $d\pi\pi$ system $2.38~GeV$, that corresponds to the isocalar $J^P = 3^+$ dibaryon $D_{03}$, and also enhancement in the distribution over the invariant mass of two final pions. The two-resonance model involving the $t$-channel $\sigma$-meson exchange between the proton and deuteron in the subprocess $pd\to pD_{03}$ and the sequential decays $D_{03}\to D_{12}+\pi$ and $D_{12}\to d +\pi$ was applied to describe the shapes of these distributions with the lowest orbital angular momenta in the corresponding vertices. A possible role of higher orbital momenta in those vertices is studied here.

  1 citation

• A somewhat random walk through nuclear and particle physics

  Thomas D. Cohen, Nicholas R. Poniatowski

  SciPost Phys. Lect. Notes 34 (2021) · published 20 September 2021 | Toggle abstract · pdf

  These notes are an outgrowth of an advanced undergraduate course taught at the University of Maryland, College Park. They are intended as an introduction to various aspects of particle and nuclear physics with an emphasis on the role of symmetry. The basic philosophy is to introduce many of the fundamental ideas in nuclear and particle physics using relatively sophisticated mathematical tools -- but to do so in as a simplified a context to explain the underlying ideas. Thus, for example, the Higgs mechanism is discussed in terms of an Abelian Higgs model. The emphasis is largely, but not entirely theoretical in orientation. The goal is for readers to develop an understanding of many of the underlying issues in a relatively sophisticated way.

  1 citation

• Scattering theory of higher order topological phases

  R. Johanna Zijderveld, Isidora Araya Day, Anton R. Akhmerov

  SciPost Phys. 19, 058 (2025) · published 27 August 2025 | Toggle abstract · pdf

  The surface states of intrinsic higher order topological phases are protected by the spatial symmetries of a finite sample. This property makes the existing scattering theory of topological invariants inapplicable: the scattering geometry is either incompatible with the symmetry or does not probe the bulk topology. We resolve this obstacle by using a symmetric scattering geometry that probes transport from the inside to the outside of the sample. We demonstrate that the intrinsic higher order topology is captured by the flux dependence of the reflection matrix. Our finding follows from identifying the spectral flow of a flux line as a signature of higher order topology. We show how this scattering approach applies to several examples of higher order topological insulators and superconductors. Our theory provides an alternative approach for proving bulk-edge correspondence in intrinsic higher order topological phases, especially in the presence of disorder.

  1 citation

• Variance reduction via simultaneous importance sampling and control variates techniques using vegas

  Prasanth Shyamsundar, Jacob L. Scott, Stephen Mrenna, Konstantin T. Matchev, Kyoungchul Kong

  SciPost Phys. Codebases 28 (2024) · published 4 March 2024 | Toggle abstract · pdf

  Monte Carlo (MC) integration is an important calculational technique in the physical sciences. Practical considerations require that the calculations are performed as accurately as possible for a given set of computational resources. To improve the accuracy of MC integration, a number of useful variance reduction algorithms have been developed, including importance sampling and control variates. In this work, we demonstrate how these two methods can be applied simultaneously, thus combining their benefits. We provide a python wrapper, named CoVVVR, which implements our approach in the VEGAS program. The improvements are quantified with several benchmark examples from the literature.

  1 citation

• Quantum robustness of the toric code in a parallel field on the honeycomb and triangular lattice

  Viktor Kott, Matthias Mühlhauser, Jan Alexander Koziol, Kai Phillip Schmidt

  SciPost Phys. 17, 053 (2024) · published 15 August 2024 | Toggle abstract · pdf

  We investigate the quantum robustness of the topological order in the toric code on the honeycomb lattice in the presence of a uniform parallel field. For a field in $z$-direction, the low-energy physics is in the flux-free sector and can be mapped to the transverse-field Ising model on the honeycomb lattice. One finds a second-order quantum phase transition in the 3D Ising$^\star$ universality class for both signs of the field. The same is true for a postive field in $x$-direction where an analogue mapping in the charge-free sector yields a ferromagnetic transverse-field Ising model on the triangular lattice and the phase transition is still 3D Ising$^\star$. In contrast, for negative $x$-field, the charge-free sector is mapped to the highly frustrated antiferromagnetic transverse-field Ising model on the triangular lattice which is known to host a quantum phase transition in the 3D XY$^\star$ universality class. Further, the charge-free sector does not always contain the low-energy physics for negative $x$-fields and a first-order phase transition to the polarized phase in the charge-full sector takes place at larger negative field values. We quantify the location of this transition by comparing quantum Monte Carlo simulations and high-field series expansions. The full extension of the topological phase in the presence of $x$- and $z$-fields is determined by perturbative linked-cluster expansions using a full graph decomposition. Extrapolating the high-order series of the charge and the flux gap allows to estimate critical exponents of the gap closing. This analysis indicates that the topological order breaks down by critical lines of 3D Ising$^\star$ and 3D XY$^\star$ type with interesting potential multi-critical crossing points. All findings for the toric code on the honeycomb lattice can be transferred exactly to the toric code on the triangular lattice.

  1 citation

• Axions in the early universe

  Francesco D'Eramo

  SciPost Phys. Proc. 12, 005 (2023) · published 3 July 2023 | Toggle abstract · pdf

  The Peccei-Quinn solution to the strong CP problem provides a motivated framework rich in cosmological consequence. Thermal axion production is unavoidable if there is a thermal bath at early times. Scattering and decay processes of bath particles can dump relativistic axions in the primordial plasma, and they can leave observable signatures in cosmological observables probing both the early and the late universe if produced with a significant abundance. We present recent and significant improvements for the calculation of the axion production rate for different scenarios and apply these results to predict the abundance of produced axions. Finally, we provide updated cosmological bounds on the QCD axion mass.

  1 citation

• Estimating time in quantum chaotic systems and black holes

  Haifeng Tang, Shreya Vardhan, Jinzhao Wang

  SciPost Phys. 19, 095 (2025) · published 13 October 2025 | Toggle abstract · pdf

  We characterize new universal features of the dynamics of chaotic quantum many-body systems, by considering a hypothetical task of "time estimation". Most macroscopic observables in a chaotic system equilibrate to nearly constant late-time values. Intuitively, it should become increasingly difficult to estimate the precise value of time by making measurements on the state. We use a quantity called the Fisher information from quantum metrology to quantify the minimum uncertainty in estimating time. Due to unitarity, the uncertainty in the time estimate does not grow with time if we have access to optimal measurements on the full system. Restricting the measurements to act on a small subsystem or to have low computational complexity leads to results expected from equilibration, where the time uncertainty becomes large at late times. With optimal measurements on a subsystem larger than half of the system, we regain the ability to estimate the time very precisely, even at late times. Hawking's calculation for the reduced density matrix of the black hole radiation in semiclassical gravity contradicts our general predictions for unitary quantum chaotic systems. Hawking's state always has a large uncertainty for attempts to estimate the time using the radiation, whereas our general results imply that the uncertainty should become small after the Page time. This gives a new version of the black hole information loss paradox in terms of the time estimation task. By restricting to simple measurements on the radiation, the time uncertainty becomes large. This indicates from a new perspective that the observations of computationally bounded agents are consistent with the semiclassical effective description of gravity.

  1 citation

• Variational solution of the superconducting Anderson impurity model and the band-edge singularity phenomena

  Teodor Iličin, Rok Žitko

  SciPost Phys. 19, 006 (2025) · published 2 July 2025 | Toggle abstract · pdf

  We propose a set of variational wavefunctions for the sub-gap spin-doublet and spin-singlet eigenstates of the particle-hole symmetric superconducting Anderson impurity model. The wavefunctions include up to two Bogoliubov quasiparticles in the continuum which is necessary to correctly capture the weak-coupling asymptotics in all parameter regimes. The eigenvalue problems reduce to solving transcendental equations. We investigate how the lowest singlet state evolves with increasing charge repulsion $U$, transitioning from a proximitized state (a superposition of empty and doubly occupied impurity orbitals, corresponding to an Andreev bound state) to a local moment that is Kondo screened by Bogoliubov quasiparticles (Yu-Shiba-Rusinov state). This change occurs for $U = 2\Delta$, where $\Delta$ is the BCS gap. At this point, the band-edge effects make the eigenenergy scale in a singular way as $\Gamma^{2/3}$, where $\Gamma$ is the hybridization strength. Away from this special point, regular $\Gamma$-linear behavior is recovered, but only for $\Gamma ≲ (U/2-\Delta)^2/\Delta$. The singular behavior thus extends over a broad range of parameters, including those relevant for some quantum devices in current use. The singular state is an equal-superposition state with maximal fluctuations between the local impurity charge configurations. Accurately capturing the band-edge singularity requires a continuum model, and it cannot be correctly described by discrete (truncated) models such as the zero-bandwidth approximation or the superconducting atomic limit. We determine the region of parameter space where the second spin-singlet state exists: in addition to the whole $U<2\Delta$ ABS region, it also includes a small part of the $U>2\Delta$ YSR region for finite values of $\Gamma$, as long as some ABS wavefunction component is admixed.

  1 citation

• Weak vs. strong breaking of integrability in interacting scalar quantum field theories

  Bence Fitos, Gábor Takács

  SciPost Phys. 15, 137 (2023) · published 5 October 2023 | Toggle abstract · pdf

  The recently proposed classification of integrability-breaking perturbations according to their strength is studied in the context of quantum field theories. Using random matrix methods to diagnose the resulting quantum chaotic behaviour, we investigate the $\phi^4$ and $\phi^6$ interactions of a massive scalar, by considering the crossover between Poissonian and Wigner-Dyson distributions in systems truncated to a finite-dimensional Hilbert space. We find that a naive extension of the scaling of crossover coupling with the volume observed in spin chains does not give satisfactory results for quantum field theory. Instead, we demonstrate that considering the scaling of the crossover coupling with the number of particles yields robust signatures, and is able to distinguish between the strengths of integrability breaking in the $\phi^4$ and $\phi^6$ quantum field theories.

  1 citation