SciPost: SciPost Physics: issue detail

SciPost Physics

SciPost Physics Vol. 19 issue 4 (October 2025)

Previous issue | Vol. 19 issue 3

Vol. 19 issue 5 | Next issue

Holographic tensor network for double-scaled SYK

Kazumi Okuyama

SciPost Phys. 19, 083 (2025) · published 2 October 2025 | · pdf

We construct a holographic tensor network for the double-scaled SYK model (DSSYK). The moment of the transfer matrix of DSSYK can be mapped to the matrix product state (MPS) of a spin chain. By adding the height direction as a holographic direction, we recast the MPS for DSSYK into the holographic tensor network whose building block is a 4-index tensor with the bond dimension three.
Cryptographic tests of the python’s lunch conjecture

Alex May, Sabrina Pasterski, Chris Waddell, Michelle Xu

SciPost Phys. 19, 084 (2025) · published 2 October 2025 | · pdf

In the AdS/CFT correspondence, a subregion of the CFT allows for the recovery of a corresponding subregion of the bulk known as its entanglement wedge. In some cases, an entanglement wedge contains a locally but not globally minimal surface homologous to the CFT subregion, in which case it is said to contain a python's lunch. It has been proposed that python's lunch geometries should be modelled by tensor networks that feature projective operations where the wedge narrows. This model leads to the python's lunch (PL) conjecture, which asserts that reconstructing information from past the locally minimal surface is computationally difficult. In this work, we use cryptographic tools related to a primitive known as the Conditional Disclosure of Secrets (CDS) to develop consequences of the projective tensor network model that can be checked directly in AdS/CFT. We argue from the tensor network picture that the mutual information between appropriate CFT subregions is lower bounded linearly by an area difference associated with the geometry of the lunch. Recalling that the mutual information is also computed by bulk extremal surfaces, this gives a checkable geometrical consequence of the tensor network model. We prove weakened versions of this geometrical statement in asymptotically AdS$_{2+1}$ spacetimes satisfying the null energy condition, and confirm it in some example geometries, supporting the tensor network model and by proxy the PL conjecture.
Efficient mutual magic and magic capacity with matrix product states

Poetri Sonya Tarabunga, Tobias Haug

SciPost Phys. 19, 085 (2025) · published 6 October 2025 | · 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.
Long-time divergences in the nonlinear response of gapped one-dimensional many-particle systems

Michele Fava, Sarang Gopalakrishnan, Romain Vasseur, Siddharth A. Parameswaran, Fabian H. L. Essler

SciPost Phys. 19, 086 (2025) · published 6 October 2025 | · pdf

We consider one dimensional many-particle systems that exhibit kinematically protected single-particle excitations over their ground states. We show that momentum and time-resolved 4-point functions of operators that create such excitations diverge linearly in particular time differences. This behaviour can be understood by means of a simple semiclassical analysis based on the kinematics and scattering of wave packets of quasiparticles. We verify that our wave packet analysis correctly predicts the long-time limit of the four-point function in the transverse field Ising model through a form factor expansion. We present evidence in favour of the same behaviour in integrable quantum field theories. In addition, we extend our discussion to experimental protocols where two times of the four-point function coincide, e.g. 2D coherent spectroscopy and pump-probe experiments. Finally, focusing on the Ising model, we discuss subleading corrections that grow as the square root of time differences. We show that the subleading corrections can be correctly accounted for by the same semiclassical analysis, but also taking into account wave packet spreading.
How negative can null energy be in large N CFTs?

Jackson R. Fliss, Ben Freivogel, Eleni-Alexandra Kontou, Diego Pardo Santos

SciPost Phys. 19, 087 (2025) · published 7 October 2025 | · pdf

Smeared null energy has been shown to be bounded from below for free minimally coupled quantum field theories. This is not the case for conformally coupled free bosonic theories where states of unbounded null energy can be constructed by increasing the particle number. Little is known for interacting conformal field theories (CFTs) in dimensions larger than two. In this work we consider states that are superpositions of scalar primary operators or the stress-energy tensor itself in large $N$ CFTs. Within the large $N$ approximation we present arguments that the negative smeared null energy of such states scales at worst as the central charge of the theory, $C_T$. This provides evidence for a general bound for CFTs in $d$-dimensions proportional to the central charge.
Boundary layers, transport and universal distribution in boundary driven active systems

Pritha Dolai, Arghya Das

SciPost Phys. 19, 088 (2025) · published 8 October 2025 | · pdf

We discuss analytical results for a run-and-tumble particle (RTP) in one dimension in presence of boundary reservoirs. It exhibits "kinetic boundary layers", nonmonotonous distribution, current without density gradient, diffusion facilitated current reversal and optimisation on tuning dynamical parameters, and a new transport effect in the steady state. The spatial and internal degrees of freedom together possess a symmetry, using which we find the eigenspectrum for large systems. The eigenvalues are arranged in two bands which can mix in certain conditions resulting in a crossover in the relaxation. The late time distribution for large systems is obtained analytically; it retains a strong and often dominant "active" contribution in the bulk rendering an effective passive-like description inadequate. A nontrivial "Milne length" also emerges in the dynamics. Finally, a novel universality is proposed in the absorbing boundary problem for dynamics with short-range colored noise. Active processes driven by active reservoirs may thus provide a common physical ground for diverse and new nonequilibrium phenomena.
Duality-Symmetry Enhancement in Maxwell Theory

Shani Meynet, Daniele Migliorati, Raffaele Savelli, Michele Tortora

SciPost Phys. 19, 089 (2025) · published 8 October 2025 | · pdf

Free Maxwell theory on general four-manifolds may, under certain conditions on the background geometry, exhibit holomorphic factorization in its partition function. We show that when this occurs, new discrete symmetries emerge at orbifold points of the conformal manifold. These symmetries, which act only on a sublattice of flux configurations, are not associated with standard dualities, yet they may carry 't Hooft anomalies, potentially causing the partition function to vanish even in the absence of apparent pathologies. We further explore their non-invertible extensions and argue that their anomalies can account for zeros of the partition function at smooth points in the moduli space.
The open XYZ spin 1/2 chain: Separation of variables and scalar products for boundary fields related by a constraint

Giuliano Niccoli, Veronique Terras

SciPost Phys. 19, 090 (2025) · published 8 October 2025 | · pdf

We consider the open XYZ spin chain with boundary fields. We solve the model by the new Separation of Variables approach introduced in [J. M. Maillet and G. Niccoli, J. Stat. Mech.: Theory Exp. 094020 (2019)]. In this framework, the transfer matrix eigenstates are obtained as a particular sub-class of the class of so-called separate states. We consider the problem of computing scalar products of such separate states. As usual, they can be represented as determinants with rows labelled by the inhomogeneity parameters of the model. We notably focus on the special case in which the boundary parameters parametrising the two boundary fields satisfy one constraint, hence enabling for the description of part of the transfer matrix spectrum and eigenstates in terms of some elliptic polynomial $Q$-solution of a usual $TQ$-equation. In this case, we show how to transform the aforementioned determinant for the scalar product into some more convenient form for the consideration of the homogeneous and thermodynamic limits: as in the open XXX or XXZ cases, our result can be expressed as some generalisation of the so-called Slavnov determinant.
Worldsheet fermion correlators, modular tensors and higher genus integration kernels

Eric D'Hoker, Oliver Schlotterer

SciPost Phys. 19, 091 (2025) · published 9 October 2025 | · pdf

The cyclic product of an arbitrary number of Szegö kernels for even spin structure $\delta$ on a compact higher-genus Riemann surface $\Sigma$ may be decomposed via a descent procedure which systematically separates the dependence on the points $z_i ∈ \Sigma$ from the dependence on the spin structure $\delta$. In this paper, we prove two different, but complementary, descent procedures to achieve this decomposition. In the first procedure, the dependence on the points $z_i ∈ \Sigma$ is expressed via the meromorphic multiple-valued Enriquez kernels of e-print 1112.0864 while the dependence on $\delta$ resides in multiplets of functions that are independent of $z_i$, locally holomorphic in the moduli of $\Sigma$ and generally do not have simple modular transformation properties. The $\delta$-dependent constants are expressed as multiple convolution integrals over homology cycles of $\Sigma$, thereby generalizing a similar representation of the individual Enriquez kernels. In the second procedure, which was proposed without proof in e-print 2308.05044, the dependence on $z_i$ is expressed in terms the single-valued, modular invariant, but non-meromorphic DHS kernels introduced in e-print 2306.08644 while the dependence on $\delta$ resides in modular tensors that are independent of $z_i$ and are generally non-holomorphic in the moduli of $\Sigma$. Although the individual building blocks of these decompositions have markedly different properties, we show that the combinatorial structure of the two decompositions is virtually identical, thereby extending the striking correspondence observed earlier between the roles played by Enriquez and DHS kernels. Both decompositions are further generalized to the case of linear chain products of Szegö kernels.
Athermal creep deformation of ultrastable amorphous solids

Pinaki Chaudhuri, Ludovic Berthier, Misaki Ozawa

SciPost Phys. 19, 092 (2025) · published 10 October 2025 | · pdf

We numerically investigate the athermal creep deformation of amorphous materials having a wide range of stability. The imposed shear stress serves as the control parameter, allowing us to examine the time-dependent transient response through both the macroscopic strain and microscopic observables. Least stable samples exhibit monotonicity in the transient strain rate versus time, while more stable samples display a pronounced non-monotonic S-shaped curve, corresponding to failure by sharp shear band formation. We identify a diverging timescale associated with the fluidization process and extract the corresponding critical exponents. Our results are compared with predictions from existing scaling theories relevant to soft matter systems. The numerical findings for stable, brittle-like materials represent a challenge for theoretical descriptions. We monitor the microscopic initiation of shear bands during creep responses. Our study encompasses creep deformation across a variety of materials ranging from ductile soft matter to brittle metallic and oxide glasses, all within the same numerical framework.
Competition of light- and phonon-dressing in microwave-dressed Bose polarons

G. M. Koutentakis, S. I. Mistakidis, F. Grusdt, H. R. Sadeghpour, P. Schmelcher

SciPost Phys. 19, 093 (2025) · published 13 October 2025 | · pdf

We theoretically investigate the stationary properties of a spin-1/2 impurity immersed in a one-dimensional confined Bose gas. In particular, we consider coherently coupled spin states with an external field, where only one spin component interacts with the bath, enabling light dressing of the impurity and spin-dependent bath-impurity interactions. Through detailed comparisons with ab-initio many-body simulations, we demonstrate that the composite system is accurately described by a simplified effective Hamiltonian. The latter builds upon previously developed effective potential approaches in the absence of light dressing. It can be used to extract the impurity energy, residue, effective mass, and anharmonicity induced by the phononic dressing. Light-dressing is shown to increase the polaron residue, undressing the impurity from phononic excitations because of strong spin coupling. For strong repulsions, previously shown to trigger dynamical Bose polaron decay (a phenomenon called temporal orthogonality catastrophe), it is explained that strong light-dressing stabilizes a repulsive polaron-dressed state. Our results establish the effective Hamiltonian framework as a powerful tool for exploring strongly interacting polaronic systems and corroborating forthcoming experimental realizations.
Symmetries, conservation laws and entanglement in non-Hermitian fermionic lattices

Rafael D. Soares, Youenn Le Gal, Chun Y. Leung, Dganit Meidan, Alessandro Romito, Marco Schirò

SciPost Phys. 19, 094 (2025) · published 13 October 2025 | · pdf

Non-Hermitian quantum many-body systems feature steady-state entanglement transitions driven by the competition between unitary dynamics and dissipation. In this work, we reveal the fundamental role of conservation laws in shaping this competition. Focusing on translation-invariant non-interacting fermionic models with U(1) symmetry, we present a theoretical framework to understand the structure of the steady-state of these models and their entanglement content based on two ingredients: the nature of the spectrum of the non-Hermitian Hamiltonian and the constraints imposed on the steady-state single-particle occupation by the conserved quantities. These emerge from an interplay between Hamiltonian symmetries and initial state, due to the non-linearity of measurement back-action. For models with complex energy spectrum, we show that the steady state is obtained by filling single-particle right eigenstates with the largest imaginary part of the eigenvalue. As a result, one can have partially filled or fully filled bands in the steady-state, leading to an entanglement entropy undergoing a filling-driven transition between critical sub-volume scaling and area-law, similar to ground-state problems. Conversely, when the spectrum is fully real, we provide evidence that local observables can be captured using a diagonal ensemble, and the entanglement entropy exhibits a volume-law scaling independently on the initial state, akin to unitary dynamics. We illustrate these principles in the Hatano-Nelson model with periodic boundary conditions and the non-Hermitian Su-Schrieffer-Heeger model, uncovering a rich interplay between the single-particle spectrum and conservation laws in determining the steady-state structure and the entanglement transitions. These conclusions are supported by exact analytical calculations and numerical calculations relying on the Faber polynomial method.
Estimating time in quantum chaotic systems and black holes

Haifeng Tang, Shreya Vardhan, Jinzhao Wang

SciPost Phys. 19, 095 (2025) · published 13 October 2025 | · 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.
The double EFT expansion in quantum gravity

José Calderón-Infante, Alberto Castellano, Alvaro Herráez

SciPost Phys. 19, 096 (2025) · published 14 October 2025 | · pdf

In this work, we aim to characterize the structure of higher-derivative corrections within low-energy Effective Field Theories (EFTs) arising from a UV-complete theory of quantum gravity. To this end, we use string theory as a laboratory and argue that such EFTs should exhibit a double EFT expansion involving higher-curvature operators. The field-theoretic expansion is governed by the mass of the lightest (tower of) new degrees of freedom, as expected from standard field theory considerations. Conversely, the quantum-gravitational expansion is suppressed relative to the Einstein-Hilbert term by the quantum gravity cutoff, $\Lambda_{\text{QG}}$, above which no local gravitational EFT description remains valid. This structure becomes manifest in the so-called asymptotic regime, where a hierarchy between the Planck scale and $\Lambda_{\text{QG}}$ emerges, the latter identified herein as the species scale. Most notably, we demonstrate the features of the double EFT expansion through an amplitudes-based approach in (toroidal compactifications of) ten-dimensional Type IIA string theory, and via a detailed analysis of the supersymmetric black hole entropy in 4d $\mathcal{N}=2$ supergravities derived from Type II Calabi–Yau compactifications. We provide further evidence for our proposal across various string theory setups, including Calabi–Yau compactifications of M/F-theory and Type II string theory. Finally, we explore the implications of this framework for the Wilson coefficients of the aforementioned higher-curvature operators, revealing potentially significant constraints in the asymptotic regime and highlighting a remarkable interplay with recent results from the S-matrix bootstrap program.
Free growth under tension

Chenyun Yao, Jens Elgeti

SciPost Phys. 19, 097 (2025) · published 15 October 2025 | · pdf

Ever since the ground breaking work of Trepat et al. in 2009, we know that cell colonies growing on a substrate can be under tensile mechanical stress. The origin of tension has so far been attributed to cellular motility forces being oriented outward of the colony. Works in the field mainly revolve around how this orientation of the forces can be explained, ranging from velocity alignment to self-sorting due to self-propulsion. In this work, we demonstrate that tension in growing colonies can also be explained without cellular motility forces! Using a combination of well established tissue growth simulation technique and analytical modelling, we show how tension can arise as a consequence of simple mechanics of growing tissues. Combining these models with a minimalistic motility model shows how colonies can expand while under even larger tension. Furthermore, our results and analytical models provide novel analysis procedures to identify the underlying mechanics.
Generalized charges, part II: Non-invertible symmetries and the symmetry TFT

Lakshya Bhardwaj, Sakura Schäfer-Nameki

SciPost Phys. 19, 098 (2025) · published 15 October 2025 | · pdf

Consider a $d$-dimensional quantum field theory (QFT) $\mathfrak{T}$, with a generalized symmetry $\mathcal{S}$, which may or may not be invertible. We study the action of $\mathcal{S}$ on generalized or $q$-charges, i.e. $q$-dimensional operators. The main result of this paper is that $q$-charges are characterized in terms of the topological defects of the Symmetry Topological Field Theory (SymTFT) of $\mathcal{S}$, also known as the "Sandwich Construction". The SymTFT is a $(d+1)$-dimensional topological field theory, which encodes the symmetry $\mathcal{S}$ and the physical theory in terms of its boundary conditions. Our proposal applies quite generally to any finite symmetry $\mathcal{S}$, including non-invertible, categorical symmetries. Mathematically, the topological defects of the SymTFT form the Drinfeld Center of the symmetry category $\mathcal{S}$. Applied to invertible symmetries, we recover the result of Part I of this series of papers. After providing general arguments for the identification of $q$-charges with the topological defects of the SymTFT, we develop this program in detail for QFTs in 2d (for general fusion category symmetries) and 3d (for fusion 2-category symmetries).
Statistical physics of the polarised IKKT matrix model

Sean A. Hartnoll, Jun Liu

SciPost Phys. 19, 099 (2025) · published 16 October 2025 | · pdf

The polarised IKKT matrix model is the worldpoint theory of $N$ D-instantons in a background three-form flux of magnitude $\Omega$, and promises to be a highly tractable model of holography. The matrix integral can be viewed as a statistical physics partition function with inverse temperature $\Omega^4$. At large $\Omega$ the model is dominated by a matrix configuration corresponding to a 'polarised' spherical D1-brane. We show that at a critical value of $\Omega^2 N$ the model undergoes a first order phase transition, corresponding to tunneling into a collection of well-separated D-instantons. These instantons are the remnant of a competing saddle in the high $\Omega$ phase corresponding to spherical $(p,q)$ fivebranes. We use a combination of numerical and analytical arguments to capture the different regimes of the model.
Critical dynamics and cyclic memory retrieval in non-reciprocal Hopfield networks

Shuyue Xue, Mohammad Maghrebi, George I. Mias, Carlo Piermarocchi

SciPost Phys. 19, 100 (2025) · published 16 October 2025 | · pdf

We study Hopfield networks with non-reciprocal coupling inducing switches between memory patterns. Dynamical phase transitions occur between phases of no memory retrieval, retrieval of multiple point-attractors, and limit-cycles. The limit cycle phase is bounded by a Hopf bifurcation line and a fold bifurcation line. Autocorrelation scales as $\tilde{C}(\tau/N^\zeta)$, with $\zeta = 1/2$ on the Hopf line and $\zeta = 1/3$ on the fold line. Perturbations of strength $F$ on the Hopf line exhibit response times scaling as $|F|^{-2/3}$, while they induce switches in a controlled way within times scaling as $|F|^{-1/2}$ in the fold line. A Master Equation approach numerically verifies the critical behavior predicted analytically. We discuss how these networks could model biological processes near a critical threshold of cyclic instability evolving through multi-step transitions.
Proposal for an autonomous quantum heat engine

Miika Rasola, Vasilii Vadimov, Tuomas Uusnäkki, Mikko Möttönen

SciPost Phys. 19, 101 (2025) · published 17 October 2025 | · pdf

We propose and theoretically analyse a superconducting electric circuit which can be used to experimentally realize an autonomous quantum heat engine. Using a quasiclassical, non-Markovian theoretical model, we demonstrate that coherent microwave power generation can emerge solely from the heat flow through the circuit determined by non-linear circuit quantum electrodynamics. The predicted energy generation rate is sufficiently high for experimental observation with contemporary techniques, rendering this work a significant step toward the first experimental realization of an autonomous quantum heat engine based on Otto cycles.
Quantum chaos, randomness and universal scaling of entanglement in various Krylov spaces

Hai-Long Shi, Augusto Smerzi, Luca Pezzè

SciPost Phys. 19, 102 (2025) · published 17 October 2025 | · pdf

Multipartite entanglement is a crucial resource for advancing quantum technologies, with considerable research efforts directed toward achieving its rapid and scalable generation. In this work, we derive an analytical expression for the time-averaged quantum Fisher information (QFI), enabling the detection of scalable multipartite entanglement dynamically generated by all quantum chaotic systems governed by Dyson's ensembles. Our approach integrates concepts of randomness and quantum chaos, demonstrating that the QFI is universally determined by the structure and dimension of the Krylov space that confines the chaotic dynamics. In particular, the QFI ranges from $N^2/3$ for $N$ qubits in the permutation-symmetric subspace (e.g. for chaotic kicked top models with long-range interactions), to $N$ when the dynamics extend over the full Hilbert space with or without bit reversal symmetry or parity symmetry (e.g. in chaotic models with short-range Ising-like interactions). In the former case, the QFI reveals multipartite entanglement among $N/3$ qubits and highlights the power of chaotic collective spin systems in generating scalable multipartite entanglement. Interestingly this result can be related to isotropic substructures in the Wigner distribution of chaotic states and demonstrates the efficacy of quantum chaos for Heisenberg-scaling quantum metrology. Finally, our general expression for the QFI agrees with that obtained for random states and, differently from out-of-time-order-correlators, it can also distinguish chaotic from integrable unstable spin dynamics.
Rademacher expansion of modular integrals

Marco Maria Baccianti, Jeevan Chandra, Lorenz Eberhardt, Thomas Hartman, Sebastian Mizera

SciPost Phys. 19, 103 (2025) · published 20 October 2025 | · pdf

We develop a method to evaluate integrals of non-holomorphic modular functions over the fundamental domain of the torus with modular parameter $\tau$ analytically. It proceeds in two steps: first the integral is transformed to a Lorentzian contour by the same strategy that leads to the Lorentzian inversion formula in CFT, and then we apply a two-dimensional version of the Rademacher expansion. This computes the integral in terms of an expansion sensitive to the singular behaviour of the integrand near all the Lorentzian cusps $\tau \to i ∞$, $\bar{\tau} \to x ∈ \mathbb{Q}$. We apply this technique to a variety of examples such as the evaluation of string one-loop partition functions, where it leads to the first analytic formula for the cosmological constants of the bosonic string and the $\mathrm{SO}(16) × \mathrm{SO}(16)$ string.
Post-hoc reweighting of hadron production in the Lund string model

Benoît Assi, Christan Bierlich, Philip Ilten, Tony Menzo, Stephen Mrenna, Manuel Szewc, Michael K. Wilkinson, Ahmed Youssef, Jure Zupan

SciPost Phys. 19, 104 (2025) · published 20 October 2025 | · pdf

We present a method for reweighting flavor selection in the Lund string fragmentation model. This is the process of calculating and applying event weights enabling fast and exact variation of hadronization parameters on pre-generated event samples. The procedure is post hoc, requiring only a small amount of additional information stored per event, and allowing for efficient estimation of hadronization uncertainties without repeated simulation. Weight expressions are derived from the hadronization algorithm itself, and validated against direct simulation for a wide range of observables and parameter shifts. The hadronization algorithm can be viewed as a hierarchical Markov process with stochastic rejections, a structure common to many complex simulations outside of high-energy physics. This perspective makes the method modular, extensible, and potentially transferable to other domains. We demonstrate the approach in Pythia, including both coverage considerations and timing benefits. For the purpose of this paper, our goal is to develop and demonstrate the the formalism, and we therefore exclude several model variations for baryon production (popcorn model, junction production) needed for proton collisions. These will be the topic of a future paper.
Higgs phases and boundary criticality

Kristian Tyn Kai Chung, Rafael Flores-Calderón, Rafael C. Torres, Pedro Ribeiro, Sergej Moroz, Paul McClarty

SciPost Phys. 19, 105 (2025) · published 21 October 2025 | · pdf

Motivated by recent work connecting Higgs phases to symmetry protected topological (SPT) phases, we investigate the interplay of gauge redundancy and global symmetry in lattice gauge theories with Higgs fields in the presence of a boundary. The core conceptual point is that a global symmetry associated to a Higgs field, which is pure-gauge in a closed system, acts physically at the boundary under boundary conditions which allow electric flux to escape the system. We demonstrate in both Abelian and non-Abelian models that this symmetry is spontaneously broken in the Higgs regime, implying the presence of gapless edge modes. Starting with the U(1) Abelian Higgs model in 4D, we demonstrate a boundary phase transition in the 3D XY universality class separating the bulk Higgs and confining regimes. Varying the boundary coupling while preserving the symmetries shifts the location of the boundary phase transition. We then consider non-Abelian gauge theories with fundamental and group-valued Higgs matter, and identify the analogous non-Abelian global symmetry acting on the boundary generated by the total color charge. For $\mathrm{SU}{N}$ gauge theory with fundamental Higgs matter we argue for a boundary phase transition in the $\mathrm{O}{2N}$ universality class, verified numerically for $N=2,3$. For group-valued Higgs matter, the boundary theory is a principal chiral model exhibiting chiral symmetry breaking. We further demonstrate this mechanism in theories with higher-form Higgs fields. We show how the higher-form matter symmetry acts at the boundary and can spontaneously break, exhibiting a boundary confinement-deconfinement transition. We also study the electric-magnetic dual theory, demonstrating a dual magnetic defect condensation transition at the boundary. We discuss some implications and extensions of these findings and what they may imply for the relation between Higgs and SPT phases.
The 3d $A$-model and generalised symmetries, part I: Bosonic Chern-Simons theories

Cyril Closset, Elias Furrer, Adam Keyes, Osama Khlaif

SciPost Phys. 19, 106 (2025) · published 21 October 2025 | · pdf

The 3d $A$-model is a two-dimensional approach to the computation of supersymmetric observables of three-dimensional $\mathcal{N}=2$ supersymmetric gauge theories. In principle, it allows us to compute half-BPS partition functions on any compact Seifert three-manifold (as well as of expectation values of half-BPS lines thereon), but previous results focussed on the case where the gauge group $\widetilde G$ is a product of simply-connected and/or unitary gauge groups. We are interested in the more general case of a compact gauge group $G=\widetilde G/\Gamma$, which is obtained from the $\widetilde G$ theory by gauging a discrete one-form symmetry. In this paper, we discuss in detail the case of pure $\mathcal{N}=2$ Chern-Simons theories (without matter) for simple groups $G$. When $G=\widetilde G$ is simply-connected, we demonstrate the exact matching between the supersymmetric approach in terms of Seifert fibering operators and the 3d TQFT approach based on topological surgery in the infrared Chern-Simons theory $\widetilde G_k$, including through the identification of subtle counterterms that relate the two approaches. We then extend this discussion to the case where the Chern-Simons theory $G_k$ can be obtained from $\widetilde G_k$ by the condensation of abelian anyons which are bosonic. Along the way, we revisit the 3d $A$-model formalism by emphasising its 2d TQFT underpinning.
Sequential adiabatic generation of chiral topological states

Xie Chen, Michael Hermele, David T. Stephen

SciPost Phys. 19, 107 (2025) · published 23 October 2025 | · pdf

In previous work, it was shown that non-trivial gapped states can be generated from a product state using a sequential quantum circuit. Explicit circuit constructions were given for a variety of gapped states at exactly solvable fixed points. In this paper, we show that a similar generation procedure can be established for chiral topological states as well, despite the fact that they lack a zero-correlation-length exactly solvable form. Instead of sequentially applying local unitary gates, we sequentially evolve the Hamiltonian by changing local terms in one subregion and then the next. The Hamiltonian remains gapped throughout the process, giving rise to an adiabatic evolution mapping the ground state from a product state to a chiral topological state. We demonstrate such a sequential adiabatic generation process for free fermion chiral states like the Chern Insulator and the $p+ip$ superconductor. Moreover, we show that coupling a quantum state to a discrete gauge group can be achieved through a sequential quantum circuit, thereby generating interacting chiral topological states from the free fermion ones.
A Lorentz-equivariant transformer for all of the LHC

Johann Brehmer, Víctor Bresó, Pim de Haan, Tilman Plehn, Huilin Qu, Jonas Spinner, Jesse Thaler

SciPost Phys. 19, 108 (2025) · published 23 October 2025 | · pdf

We show that the Lorentz-Equivariant Geometric Algebra Transformer (L-GATr) yields state-of-the-art performance for a wide range of machine learning tasks at the Large Hadron Collider. L-GATr represents data in a geometric algebra over space-time and is equivariant under Lorentz transformations. The underlying architecture is a versatile and scalable transformer, which is able to break symmetries if needed. We demonstrate the power of L-GATr for amplitude regression and jet classification, and then benchmark it as the first Lorentz-equivariant generative network. For all three LHC tasks, we find significant improvements over previous architectures.
Precise quantum-geometric electronic properties from first principles

José Luís Martins, Carlos Loia Reis, Ivo Souza

SciPost Phys. 19, 109 (2025) · published 24 October 2025 | · pdf

The calculation of quantum-geometric properties of Bloch electrons – Berry curvature, quantum metric, orbital magnetic moment and effective mass – was implemented in a pseudopotential plane-wave code. The starting point was the first derivative of the periodic part of the wavefunction $\psi_{k}(r)$ with respect to wavevector $k$. This was evaluated with perturbation theory by solving a Sternheimer equation. Comparison of effective masses obtained from perturbation theory for silicon and gallium arsenide with carefully-converged numerical second derivatives of band energies confirmed the high precision of the method. Calculations of quantum-geometric quantities for gapped graphene were performed by adding a bespoke symmetry-breaking potential to first-principles graphene. As the two bands near the opened gap are reasonably isolated, the results could be compared with those obtained from an analytical two-band model, allowing to assess the strengths and limitations of such widely-used models. The final application was trigonal tellurium, where some quantum-geometric quantities flip sign with chirality.
Noninvertible symmetry and topological holography for modulated SPT in one dimension

Jintae Kim, Yizhi You, Jung Hoon Han

SciPost Phys. 19, 110 (2025) · published 24 October 2025 | · pdf

We examine noninvertible symmetry (NIS) in one-dimensional (1D) symmetry-protected topological (SPT) phases protected by dipolar and exponential-charge symmetries, which are two key examples of modulated SPT (MSPT). To set the stage, we first study NIS in the $\mathbb{Z}_N × \mathbb{Z}_N$ cluster model, extending previous work on the $\mathbb{Z}_2 × \mathbb{Z}_2$ case. For each symmetry type (charge, dipole, exponential), we explicitly construct the noninvertible Kramers-Wannier (KW) and Kennedy-Tasaki (KT) transformations, revealing dual models with spontaneous symmetry breaking (SSB). The resulting symmetry group structure of the SSB model is rich enough that it allows the identification of other SSB models with the same symmetry. Using these alternative SSB models and KT duality, we generate novel MSPT phases distinct from those associated with the standard decorated domain wall picture, and confirm their distinctiveness by projective symmetry analyses at their interfaces. Additionally, we establish a topological-holographic correspondence by identifying the 2D bulk theories-two coupled layers of toric codes (charge), anisotropic dipolar toric codes (dipole), and exponentially modulated toric codes (exponential)-whose boundaries host the respective 1D MSPT phases.
Mutual multilinearity of nonequilibrium network currents

Sara Dal Cengio, Pedro E. Harunari, Vivien Lecomte, Matteo Polettini

SciPost Phys. 19, 111 (2025) · published 27 October 2025 | · pdf

Continuous-time Markov chains have been successful in modelling systems across numerous fields, with currents being fundamental entities that describe the flows of energy, particles, individuals, chemical species, information, or other quantities. They apply to systems described by agents transitioning between vertices along the edges of a network (at some rate in each direction). It has recently been shown by the authors that, at stationarity, a hidden linearity exists between currents that flow along edges: if one controls the current of a specific "input" edge (by tuning transition rates along it), any other current is a linear-affine function of the input current [Phys. Rev. Lett. 133, 047401 (2024)]. In this paper, we extend this result to the situation where one controls the currents of several edges, and prove that other currents are in linear-affine relation with the input ones. Two proofs with distinct insights are provided: the first relies on Kirchhoff's current law and reduces the input set inductively through graph analysis, while the second utilizes the resolvent approach via a Laplace transform in time. We obtain explicit expressions for the current-to-current susceptibilities, which allow one to map current dependencies through the network. We also verify from our expression that Kirchhoff's current law is recovered as a limiting case of our mutual linearity. Last, we uncover that susceptibilities can be obtained from fluctuations when the reference system is originally at equilibrium.
Pseudospectral implementation of the Einstein-Maxwell system

Jorge Expósito Patiño, Hannes Robert Rüter, David Hilditch

SciPost Phys. 19, 112 (2025) · published 27 October 2025 | · pdf

Electromagnetism plays an important role in a variety of applications in gravity that we wish to investigate. To that end, in this work, we present an implementation of the Maxwell equations within the adaptive-mesh pseudospectral numerical relativity code BAMPS. We perform a thorough analysis of the evolution equations as a first order symmetric hyperbolic system of PDEs. This includes both the construction of the characteristic variables for use in our penalty boundary communication scheme, as well as radiation controlling, constraint preserving outer boundary conditions which, for the first time in a numerical context, are shown to be boundary-stable. After choosing a formulation of the Maxwell constraints that we may solve for initial data, we move on to show a suite of numerical tests. Our simulations, both within the Cowling approximation, and in full non-linear evolution, demonstrate rapid convergence of error with resolution, as well as consistency with known quasinormal decay rates on the Kerr background. Finally we evolve the electrovacuum equations of motion with strong data, a good representation of typical critical collapse runs.
Hasse diagrams for gapless SPT and SSB phases with non-invertible symmetries

Lakshya Bhardwaj, Daniel Pajer, Sakura Schäfer-Nameki, Alison Warman

SciPost Phys. 19, 113 (2025) · published 28 October 2025 | · pdf

We discuss (1+1)d gapless phases with non-invertible global symmetries, also referred to as categorical symmetries. This includes gapless phases showing properties analogous to gapped symmetry protected topological (SPT) phases, known as gapless SPT (or gSPT) phases; and gapless phases showing properties analogous to gapped spontaneous symmetry broken (SSB) phases, that we refer to as gapless SSB (or gSSB) phases. We fit these gapless phases, along with gapped SPT and SSB phases, into a phase diagram describing possible deformations connecting them. This phase diagram is partially ordered and defines a so-called Hasse diagram. Based on these deformations, we identify gapless phases exhibiting symmetry protected criticality, that we refer to as intrinsically gapless SPT (igSPT) and intrinsically gapless SSB (igSSB) phases. This includes the first examples of igSPT and igSSB phases with non-invertible symmetries. Central to this analysis is the Symmetry Topological Field Theory (SymTFT), where each phase corresponds to a condensable algebra in the Drinfeld center of the symmetry category. On a mathematical note, gSPT phases are classified by functors between fusion categories, generalizing the fact that gapped SPT phases are classified by fiber functors; and gSSB phases are classified by functors from fusion to multi-fusion categories. Finally, our framework can be applied to understand gauging of trivially acting non-invertible symmetries, including possible patterns of decomposition arising due to such gaugings.
Superstring amplitudes meet surfaceology

Qu Cao, Jin Dong, Song He, Fan Zhu

SciPost Phys. 19, 114 (2025) · published 28 October 2025 | · pdf

We reformulate tree-level amplitudes in open superstring theory (type-I) in terms of stringy Tr$(\phi^3)$ amplitudes with various kinematical shifts in the "curve-integral" formulation: while the bosonic-string amplitude with $n$ pairs of "scaffolding" scalars comes from a particularly simple shift of the Tr$(\phi^3)$ one (corresponding to $n$ length-$2$ cycles), the analogous superstring amplitude requires "correction" terms given by bosonic-string amplitudes with longer, even-length "cycles", which are also Tr$(\phi^3)$ ones at shifted kinematics dictated by the cycles; in total it is expressed as a sum of $(2n-3)!!$ shifted amplitudes originated from the expansion of a reduced Pfaffian. Upon taking $n$ scaffolding residues, this leads to a new formula of the $n$-gluon superstring amplitude, which is manifestly symmetric in $n-1$ legs, as a gauge-invariant combination of mixed bosonic string amplitudes with gluons and scalars, which come from length-$2$ cycles and longer ones respectively (the total sum is associated with the expansion a $n× n$ symmetrical determinant); the corresponding prefactors are nested commutators of $2n$-gon kinematical variables, which nicely become traces of field-strengths for those legs corresponding to scalars in the mixed amplitudes. These interesting linear combinations of bosonic string amplitudes must guarantee the cancellation of tachyon poles and $F^3$ vertices etc., and they give new relations between the superstring amplitude and its bosonic-string building blocks to all orders in the $\alpha'$ expansion (the first order gives a new formula for gluon amplitudes with a single $F^3$ insertion in terms of Yang-Mills-scalar amplitudes). We provide both the worldsheet and "curve-integral" derivations, and discuss applications to heterotic and type II cases.
The complex Liouville string: The gravitational path integral

Scott Collier, Lorenz Eberhardt, Beatrix Mühlmann

SciPost Phys. 19, 115 (2025) · published 28 October 2025 | · pdf

We give a rigorous definition of sine dilaton gravity in terms of the worldsheet theory of the complex Liouville string [S. Collier et al., Phys. Rev. Lett. 134, 251602 (2025)]. The latter has a known exact solution that we leverage to explore the gravitational path integral of sine dilaton gravity – a quantum deformation of dS JT gravity that admits both AdS$_2$ and dS$_2$ vacua. We uncover that the gravitational path integral receives contributions from new saddles describing transitions between vacua in a third-quantized picture. We also discuss the sphere and disk partition function in this context and contrast our findings with other recent work on this theory.
Symmetry operators and gravity

Ibrahima Bah, Patrick Jefferson, Konstantinos Roumpedakis, Thomas Waddleton

SciPost Phys. 19, 116 (2025) · published 29 October 2025 | · pdf

We argue that topological operators for continuous symmetries written in terms of currents need regularization, which effectively gives them a small but finite width. The regulated operator is a finite tension object which fluctuates. In the zero-width limit these fluctuations freeze, recovering the properties of a topological operator. When gravity is turned on, the zero-width limit becomes ill-defined, thereby prohibiting the existence of topological operators.
The effect of dephasing and spin-lattice relaxation during the switching processes in quantum antiferromagnets

Asliddin Khudoyberdiev, Götz S. Uhrig

SciPost Phys. 19, 117 (2025) · published 29 October 2025 | · pdf

The control of antiferromagnetic order can pave the way to large storage capacity as well as fast manipulation of stored data. Here achieving a steady-state of sublattice magnetization after switching is crucial to prevent loss of stored data. The present theoretical approach aims to obtain instantaneous stable states of the order after reorienting the Néel vector in open quantum antiferromagnets using time-dependent Schwinger boson mean-field theory. The Lindblad formalism is employed to couple the system to the environment. The quantum theoretical approach comprises differences in the effects of dephasing, originating from destructive interference of different wave vectors, and spin-lattice relaxation. We show that the spin-lattice relaxation results in an exponentially fast convergence to the steady-state after full ultrafast switching.
Parallel assembly of neutral atom arrays with an SLM using linear phase interpolation

Ivo H. A. Knottnerus, Yu Chih Tseng (曾于誌), Alexander Urech, Robert J. C. Spreeuw, Florian Schreck

SciPost Phys. 19, 118 (2025) · published 31 October 2025 | · pdf

We present fast parallel rearrangement of single atoms in optical tweezers into arbitrary geometries by updating holograms displayed by an ultra fast spatial light modulator. Using linear interpolation of the tweezer position and the optical phase between the start and end arrays, we can calculate and display holograms every few ms, limited by technology. To show the versatility of our method, we sort the same atomic sample into multiple geometries with success probabilities of $0.996(2)$ per rearrangement cycle. This makes the method a useful tool for rearranging large atom arrays for quantum computation and quantum simulation.