Publications

Order by: publication date number of citations

• Boosting determinant quantum Monte Carlo with submatrix updates: Unveiling the phase diagram of the 3D Hubbard model

  Fanjie Sun, Xiao Yan Xu

  SciPost Phys. 18, 055 (2025) · published 17 February 2025 | Toggle abstract · pdf

  Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally developed for Hirsch-Fye QMC with onsite interactions at finite temperatures, their comprehensive application in DQMC has remained unexplored despite noted algorithmic similarities. We present the first comprehensive application of submatrix updates in DQMC, significantly extending beyond the original scope by enabling simulations with extended interactions and at zero temperature. Building upon conventional fast updates and delay updates, our generalized implementation achieves an order-of-magnitude improvement in computational efficiency, enabling simulations of the half-filled Hubbard model on lattices up to 8,000 sites - a scale previously challenging with standard DQMC implementations. This enhanced computational capability allows us to accurately determine the finite-temperature phase diagram of the 3D Hubbard model at half-filling. Our findings not only shed light on the phase transitions within these complex systems but also pave the way for more effective simulations of strongly correlated electrons, potentially guiding experimental efforts in cold atom simulations of the 3D Hubbard model.

  1 citation

• Describing hadronization via histories and observables for Monte-Carlo event reweighting

  Christian Bierlich, Phil Ilten, Tony Menzo, Stephen Mrenna, Manuel Szewc, Michael K. Wilkinson, Ahmed Youssef, Jure Zupan

  SciPost Phys. 18, 054 (2025) · published 17 February 2025 | Toggle abstract · pdf

  We introduce a novel method for extracting a fragmentation model directly from experimental data without requiring an explicit parametric form, called Histories and Observables for Monte-Carlo Event Reweighting (HOMER), consisting of three steps: the training of a classifier between simulation and data, the inference of single fragmentation weights, and the calculation of the weight for the full hadronization chain. We illustrate the use of HOMER on a simplified hadronization problem, a $q\bar{q}$ string fragmenting into pions, and extract a modified Lund string fragmentation function $f(z)$. We then demonstrate the use of HOMER on three types of experimental data: (i) binned distributions of high-level observables, (ii) unbinned event-by-event distributions of these observables, and (iii) full particle cloud information. After demonstrating that $f(z)$ can be extracted from data (the inverse of hadronization), we also show that, at least in this limited setup, the fidelity of the extracted $f(z)$ suffers only limited loss when moving from (i) to (ii) to (iii). Public code is available at https://gitlab.com/uchep/mlhad.

  2 citations

• Quenching a Fermi superfluid across the BEC-BCS crossover

  Moritz Breyer, Daniel Eberz, Andreas Kell, Michael Köhl

  SciPost Phys. 18, 053 (2025) · published 14 February 2025 | Toggle abstract · pdf

  We study the response of a Fermi superfluid to a rapid change of the interaction strength. For a broad range of quench parameters, the order parameter exhibits damped oscillations, however, we also identify quench regimes in which these oscillations are absent. By comparing our data to a numerical model we find that the damping time constants are explainable by a dephasing in a local density approximation, however, the oscillation frequencies, while being comparable to the superconducting order parameter, display in detail significant differences to the integrable BCS (Bardeen-Cooper-Schrieffer) model.

• Necklace Ansatz for strongly repulsive spin mixtures on a ring

  Gianni Aupetit-Diallo, Giovanni Pecci, Artem Volosniev, Mathias Albert, Anna Minguzzi, Patrizia Vignolo

  SciPost Phys. Core 8, 022 (2025) · published 14 February 2025 | Toggle abstract · pdf

  We propose an alternative to the Bethe Ansatz method for repulsive strongly-interacting fermionic (or bosonic) mixtures on a ring. Starting from the knowledge of the solution for single-component non-interacting fermions (or strongly-interacting bosons), we explicitly impose periodic condition on the amplitudes of the spin configurations. This reduces drastically the number of independent complex amplitudes that we determine by constrained diagonalization of an effective Hamiltonian. This procedure allows us to obtain a complete basis for the exact low-energy many-body solutions for mixtures with a large number of particles, both for $SU(\kappa)$ and symmetry-breaking systems.

• Fast Perfekt: Regression-based refinement of fast simulation

  Moritz Wolf, Lars O. Stietz, Patrick L. S. Connor, Peter Schleper, Samuel Bein

  SciPost Phys. Core 8, 021 (2025) · published 14 February 2025 | Toggle abstract · pdf

  The availability of precise and accurate simulation is a limiting factor for interpreting and forecasting data in many fields of science and engineering. Often, one or more distinct simulation software applications are developed, each with a relative advantage in accuracy or speed. The quality of insights extracted from the data stands to increase if the accuracy of faster, more economical simulation could be improved to parity or near parity with more resource-intensive but accurate simulation. We present Fast Perfekt, a machine-learned regression to refine the output of fast simulation that employs residual neural networks. A deterministic morphing model is trained using a unique schedule that makes use of the ensemble loss function MMD, with the option of an additional pair-based loss function such as the MSE. We explore this methodology in the context of an abstract analytical model and in terms of a realistic particle physics application featuring jet properties in hadron collisions at the CERN Large Hadron Collider. The refinement makes maximum use of existing domain knowledge, and introduces minimal computational overhead to production.

  1 citation

• Higher-matter and Landau-Ginzburg theory of higher-group symmetries

  Ruizhi Liu, Ran Luo, Yi-Nan Wang

  SciPost Phys. 18, 052 (2025) · published 13 February 2025 | Toggle abstract · pdf

  Higher-matter is defined by higher-representation of a symmetry algebra, such as the $p$-form symmetries, higher-group symmetries or higher-categorical symmetries. In this paper, we focus on the cases of higher-group symmetries, which are formulated in terms of the strictification of weak higher-groups. We systematically investigate higher-matter charged under 2-group symmetries, defined by automorphism 2-representations. Furthermore, we construct a Lagrangian formulation of such higher-matter fields coupled to 2-group gauge fields in the path space of the spacetime manifold. We interpret such model as the Landau-Ginzburg theory for 2-group symmetries, and discuss the spontaneous symmetry breaking (SSB) of 2-group symmetries under this framework. Examples of discrete and continuous 2-groups are discussed. Interestingly, we find that a non-split 2-group symmetry can admit an SSB to a split 2-group symmetry, where the Postnikov class is trivialized. We also briefly discuss the strictification of weak 3-groups, weak 3-group gauge fields and 3-representations in special cases.

  2 citations

• Critical behavior of a phase transition in the dynamics of interacting populations

  Thibaut Arnoulx de Pirey, Guy Bunin

  SciPost Phys. 18, 051 (2025) · published 13 February 2025 | Toggle abstract · pdf

  Many-variable differential equations with random coefficients provide powerful models for the dynamics of many interacting species in ecology. These models are known to exhibit a dynamical phase transition from a phase where population sizes reach a fixed point, to a phase where they fluctuate indefinitely. Here we provide a theory for the critical behavior close to the phase transition. We show that timescales diverge at the transition and that temporal fluctuations grow continuously upon crossing it. We further show the existence of three different universality classes, with different sets of critical exponents, highlighting the importance of the migration rate coupling the system to its surroundings.

• Possible causes of false general relativity violations in gravitational wave observations

  Anuradha Gupta, K. G. Arun, Enrico Barausse, Laura Bernard, Emanuele Berti, Sajad A. Bhat, Alessandra Buonanno, Vitor Cardoso, Shun Yin Cheung, Teagan A. Clarke, Sayantani Datta, Arnab Dhani, Jose María Ezquiaga, Ish Gupta, Nir Guttman, Tanja Hinderer, Qian Hu, Justin Janquart, Nathan K. Johnson-McDaniel, Rahul Kashyap, N. V. Krishnendu, Paul D. Lasky, Andrew Lundgren, Elisa Maggio, Parthapratim Mahapatra, Andrea Maselli, Purnima Narayan, Alex B. Nielsen, Laura K. Nuttall, Paolo Pani, Lachlan Passenger, Ethan Payne, Lorenzo Pompili, Luca Reali, Pankaj Saini, Anuradha Samajdar, Shubhanshu Tiwari, Hui Tong, Chris Van Den Broeck, Kent Yagi, Huan Yang, Nicolás Yunes, B. S. Sathyaprakash

  SciPost Phys. Comm. Rep. 5 (2025) · published 13 February 2025 | Toggle abstract · pdf

  General relativity (GR) has proven to be a highly successful theory of gravity since its inception. The theory has thrivingly passed numerous experimental tests, predominantly in weak gravity, low relative speeds, and linear regimes, but also in the strong-field and very low-speed regimes with binary pulsars. Observable gravitational waves (GWs) originate from regions of spacetime where gravity is extremely strong, making them a unique tool for testing GR, in previously inaccessible regions of large curvature, relativistic speeds, and strong gravity. Since their first detection, GWs have been extensively used to test GR, but no deviations have been found so far. Given GR's tremendous success in explaining current astronomical observations and laboratory experiments, accepting any deviation from it requires a very high level of statistical confidence and consistency of the deviation across GW sources. In this paper, we compile a comprehensive list of potential causes that can lead to a false identification of a GR violation in standard tests of GR on data from current and future ground-based GW detectors. These causes include detector noise, signal overlaps, gaps in the data, detector calibration, source model inaccuracy, missing physics in the source and in the underlying environment model, source misidentification, and mismodeling of the astrophysical population. We also provide a rough estimate of when each of these causes will become important for tests of GR for different detector sensitivities. We argue that each of these causes should be thoroughly investigated, quantified, and ruled out before claiming a GR violation in GW observations.

  25 citations

• Pymablock: An algorithm and a package for quasi-degenerate perturbation theory

  Isidora Araya Day, Sebastian Miles, Hugo K. Kerstens, Daniel Varjas, Anton R. Akhmerov

  SciPost Phys. Codebases 50 (2025) · published 12 February 2025 | Toggle abstract · pdf

  A common technique in the study of complex quantum-mechanical systems is to reduce the number of degrees of freedom in the Hamiltonian by using quasi-degenerate perturbation theory. While the Schrieffer-Wolff transformation achieves this and constructs an effective Hamiltonian, its scaling is suboptimal, it is limited to two subspaces, and implementing it efficiently is both challenging and error-prone. We introduce an algorithm for constructing an equivalent effective Hamiltonian as well as a Python package, Pymablock, that implements it. Our algorithm combines an optimal asymptotic scaling and the ability to handle any number of subspaces with a range of other improvements. The package supports numerical and analytical calculations of any order and it is designed to be interoperable with any other packages for specifying the Hamiltonian. We demonstrate how the package handles constructing a k.p model, analyses a superconducting qubit, and computes the low-energy spectrum of a large tight-binding model. We also compare its performance with reference calculations and demonstrate its efficiency.

  4 citations

• Wilson loops and spherical branes

  Davide Astesiano, Pieter Bomans, Fridrik Freyr Gautason, Valentina Giangreco M. Puletti, Alexia Nix

  SciPost Phys. 18, 050 (2025) · published 12 February 2025 | Toggle abstract · pdf

  We study $\frac{1}{2}$-BPS Wilson loop operators in maximally supersymmetric Yang-Mills theory on $d$-dimensional spheres. Their vacuum expectation values can be computed at large $N$ through supersymmetric localisation. The holographic duals are given by back-reacted spherical D-branes. For $d≠ 4$, the resulting theories are non-conformal and correspondingly, the dual geometries do not possess an asymptotic AdS region. The main aim of this work is to compute the holographic Wilson loops by evaluating the partition function of a probe fundamental string and M2-brane in the dual geometry, focusing on the next-to-leading order. Along the way, we highlight a variety of issues related to the presence of a non-constant dilaton. In particular, the structure of the divergences of the one-loop partition functions takes a non-universal form in contrast to examples available in the literature. We devise a general framework to treat the divergences, successfully match the sub-leading scaling with $\lambda$ and $N$, and provide a first step towards obtaining the numerical prefactor.

  4 citations