Publications

Order by: publication date number of citations

• Efficient mutual magic and magic capacity with matrix product states

  Poetri Sonya Tarabunga, Tobias Haug

  SciPost Phys. 19, 085 (2025) · published 6 October 2025 | Toggle abstract · pdf

  Stabilizer Rényi entropies (SREs) probe the non-stabilizerness (or "magic") of many-body systems and quantum computers. Here, we introduce the mutual von-Neumann SRE and magic capacity, which can be efficiently computed in time $O(N\chi^3)$ for matrix product states (MPSs) of bond dimension $\chi$. We find that mutual SRE characterizes the critical point of ground states of the transverse-field Ising model, independently of the chosen local basis. Then, we relate the magic capacity to the anti-flatness of the Pauli spectrum, which quantifies the complexity of computing SREs. The magic capacity characterizes transitions in the ground state of the Heisenberg and Ising model, randomness of Clifford+T circuits, and distinguishes typical and atypical states. Finally, we make progress on numerical techniques: we design two improved Monte-Carlo algorithms to compute the mutual $2$-SRE, overcoming limitations of previous approaches based on local update. We also give improved statevector simulation methods for Bell sampling and SREs with $O(8^{N/2})$ time and $O(2^N)$ memory, which we demonstrate for $24$ qubits. Our work uncovers improved approaches to study the complexity of quantum many-body systems.

  3 citations

• Floquet engineering of non-equilibrium superradiance

  Lukas Broers, Ludwig Mathey

  SciPost Phys. 14, 018 (2023) · published 13 February 2023 | Toggle abstract · pdf

  We demonstrate the emergence of a non-equilibrium superradiant phase in the dissipative Rabi-Dicke model. This phase is characterized by a photonic steady state that oscillates with a frequency close to the cavity frequency, in contrast to the constant photonic steady state of the equilibrium superradiant phase in the Dicke model. We relate this superradiant phase to the population inversion of Floquet states by introducing a Schwinger representation of the driven two-level systems in the cavity. This inversion is depleted near Floquet energies that are resonant with the cavity frequency to sustain a coherent light-field. In particular, our model applies to solids within a two-band approximation, in which the electrons act as Schwinger fermions. We propose to use this Floquet-assisted superradiant phase to obtain controllable optical gain for a laser-like operation.

  3 citations

• Codebase release 0.1 for jVMC

  Markus Schmitt, Moritz Reh

  SciPost Phys. Codebases 2-r0.1 (2022) · published 23 August 2022 | Toggle abstract · src

  The introduction of Neural Quantum States (NQS) has recently given a new twist to variational Monte Carlo (VMC). The ability to systematically reduce the bias of the wave function ansatz renders the approach widely applicable. However, performant implementations are crucial to reach the numerical state of the art. Here, we present a Python codebase that supports arbitrary NQS architectures and model Hamiltonians. Additionally leveraging automatic differentiation, just-in-time compilation to accelerators, and distributed computing, it is designed to facilitate the composition of efficient NQS algorithms.

  3 citations

• Transitions in Xenes between excitonic, topological and trivial insulator phases: Influence of screening, band dispersion and external electric field

  Olivia Pulci, Paola Gori, Davide Grassano, Marco D'Alessandro, Friedhelm Bechstedt

  SciPost Phys. 15, 025 (2023) · published 24 July 2023 | Toggle abstract · pdf

  Using a variational approach, the binding energies $E_b$ of the lowest bound excitons in Xenes under varying electric field are investigated. The internal exciton motion is described both by Dirac electron dispersion and in effective-mass approximation, while the screened electron-hole attraction is modeled by a Rytova-Keldysh potential with a 2D electronic polarizability $\alpha_{2D}$. The most important parameters as spin-orbit-induced gap $E_g$, Fermi velocity $v_F$ and $\alpha_{2D}$ are taken from ab initio density functional theory calculations. In addition, $\alpha_{2D}$ is approximated in two different ways. The relation of $E_b$ and $E_g$ is ruled by the screening. The existence of an excitonic insulator phase with $E_b>E_g$ sensitively depends on the chosen $\alpha_{2D}$. The values of $E_g$ and $\alpha_{2D}$ are strongly modified by a vertical external electric bias $U$, which defines a transition from the topological into a trivial insulator at $U=E_g/2$, with the exception of plumbene. Within the Dirac approximation, but also within the effective mass description of the kinetic energy, the treatment of screening dominates the appearance or non-appearance of an excitonic insulator phase. Gating does not change the results: the prediction done at zero electric field is confirmed when a vertical electric field is applied. Finally, Many-Body perturbation theory approaches based on the Green's function method, applied to stanene, confirm the absence of an excitonic insulator phase, thus validating our results obtained by ab initio modeling of $\alpha_{2D}$.

  3 citations

• Physical and unphysical regimes of self-consistent many-body perturbation theory

  Kris Van Houcke, Evgeny Kozik, Riccardo Rossi, Youjin Deng, Félix Werner

  SciPost Phys. 16, 133 (2024) · published 27 May 2024 | Toggle abstract · pdf

  In the standard framework of self-consistent many-body perturbation theory, the skeleton series for the self-energy is truncated at a finite order N and plugged into the Dyson equation, which is then solved for the propagator $G_N$. We consider two examples of fermionic models, the Hubbard atom at half filling and its zero space-time dimensional simplified version. First, we show that $G_N$ converges when $N\to∞$ to a limit $G_∞\,$, which coincides with the exact physical propagator $G_{exact}$ at small enough coupling, while $G_∞ ≠ G_{exact}$ at strong coupling. This follows from the findings of [Phys. Rev. Lett. 114, 156402 (2015)] and an additional subtle mathematical mechanism elucidated here. Second, we demonstrate that it is possible to discriminate between the $G_∞=G_{exact}$ and $G_∞≠G_{exact}$ regimes thanks to a criterion which does not require the knowledge of $G_{exact}$, as proposed in [Phys. Rev. B 93, 161102 (2016)].

  2 citations

• Charge order and antiferromagnetism in twisted bilayer graphene from the variational cluster approximation

  B. Pahlevanzadeh, P. Sahebsara, D. Sénéchal

  SciPost Phys. 13, 040 (2022) · published 30 August 2022 | Toggle abstract · pdf

  We study the possibility of charge order at quarter filling and antiferromagnetism at half-filling in a tight-binding model of magic angle twisted bilayer graphene. We build on the model proposed by Kang and Vafek, relevant to a twist angle of $1.30^\circ$, and add on-site and extended density-density interactions. Applying the variational cluster approximation with an exact-diagonalization impurity solver, we find that the system is indeed a correlated (Mott) insulator at fillings $\frac14$, $\frac12$ and $\frac34$. At quarter filling, we check that the most probable charge orders do not arise, for all values of the interaction tested. At half-filling, antiferromagnetism only arises if the local repulsion $U$ is sufficiently large compared to the extended interactions, beyond what is expected from the simplest model of extended interactions.

  2 citations

• 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.

  2 citations

• Codebase release 1.0 for SpinParser

  Finn Lasse Buessen

  SciPost Phys. Codebases 5-r1.0 (2022) · published 24 August 2022 | Toggle abstract · src

  We present the SpinParser open-source software [ https://github.com/fbuessen/SpinParser ]. The software is designed to perform pseudofermion functional renormalization group (pf-FRG) calculations for frustrated quantum magnets in two and three spatial dimensions. It aims to make such calculations readily accessible without the need to write specialized program code; instead, custom lattice graphs and microscopic spin models can be defined as plain-text input files. Underlying symmetries of the model are automatically analyzed and exploited by the numerical core written in C++ in order to optimize the performance across large-scale shared memory and/or distributed memory computing platforms.

  2 citations

• Parametrization and thermodynamic scaling of pair correlation functions for the fractional quantum Hall effect

  Jørgen Fulsebakke, Mikael Fremling, Niall Moran, Joost K. Slingerland

  SciPost Phys. 14, 149 (2023) · published 9 June 2023 | Toggle abstract · pdf

  The calculation of pair correlations and density profiles of quasiholes are routine steps in the study of proposed fractional quantum Hall states. Nevertheless, the field has not adopted a standard way to present the results of such calculations in an easily reproducible form. We develop a polynomial expansion that allows for easy quantitative comparison between different candidate wavefunctions, as well as reliable scaling of correlation and quasihole profiles to the thermodynamic limit. We start from the well-known expansion introduced by Girvin [PRB, 30 (1984)] (see also [Girvin, MacDonald and Platzman, PRB, 33 (1986)]), which is physically appealing but, as we demonstrate, numerically unstable. We orthogonalize their basis set to obtain a new basis of modified Jacobi polynomials, whose coefficients can be stably calculated. We then apply our expansion to extract pair correlation expansion coefficients and quasihole profiles in the thermodynamic limit for a wide range of fractional quantum Hall wavefunctions. These include the Laughlin series, composite fermion states with both reverse and direct flux attachment, the Moore-Read Pfaffian state, and BS hierarchy states. The expansion procedure works for both abelian and non-abelian quasiholes, even when the density at the core is not zero. We find that the expansion coefficients for all quantum Hall states considered can be fit remarkably well using a cosine oscillation with exponentially decaying amplitude. The frequency and the decay length are related in an intuitive, but not elementary way to the filling fraction. Different states at the same filling fraction can have distinct values for these parameters. Finally, we also use our scaled correlation functions to calculate estimates for the magneto-roton gaps of the various states.

  2 citations

• Analysis of radiation damage in silicon charge-coupled devices used for dark matter searches

  Steven J. Lee

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

  Nuclear recoils in crystal detectors generate radiation damage in the form of crystal defects that can be measured in scientific-grade CCDs as local hot spots of leakage current stimulated by temperature increases in the devices. In this proceeding, we use a neutron source to generate defects in DAMIC-M CCDs, and using increases in leakage current at different temperatures, we demonstrate a procedure for identifying crystal defects in the CCDs of the DAMIC-M experiment. This is the first time that individual defects generated from nuclear recoils have been studied. This technique could be used to distinguish nuclear recoils from electron recoils in some energy ranges, which would improve the ability of CCD detectors to search for weakly interacting dark matter.

  2 citations