Publications

Order by: publication date number of citations

• 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 | Toggle abstract · 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.

  1 citation

• EERAD3 version 2: QCD corrections in hadronic colour-singlet decays

  Benjamin Campillo Aveleira, Aude Gehrmann-De Ridder, Thomas Gehrmann, Nigel Glover, Gudrun Heinrich, Christian T. Preuss

  SciPost Phys. Codebases 59 (2025) · published 22 October 2025 | Toggle abstract · pdf

  We present a major update of the publicly available EERAD3 package to calculate perturbative corrections in the strong coupling in hadronic Higgs and $Z$-boson decays. We describe the theoretical framework underlying the numerical implementation and provide a guide to the usage of the program.

  1 citation

• Codebase release 1.1 for GAPS

  Michael H. Seymour, Siddharth Sule

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

  The Single Instruction, Multiple Thread (SIMT) paradigm of GPU programming does not support the branching nature of a parton shower algorithm by definition. However, modern GPUs are designed to schedule threads with diverging processes independently, allowing them to handle such branches. With regular thread synchronisation and careful treatment of the individual steps, one can simulate a parton shower on a GPU. We present a Sudakov veto algorithm designed to simulate parton branching on multiple events in parallel. We also release a CUDA C++ program that generates matrix elements, showers partons and computes jet rates and event shapes for LEP at 91.2 GeV on a GPU. To benchmark its performance, we also provide a near-identical C++ program designed to simulate events serially on a CPU. While the consequences of branching are not absent, we demonstrate that a GPU can provide the throughput of a many-core CPU. As an example, we show that the time taken to shower $10^6$ events on one NVIDIA TESLA V100 GPU is equivalent to that of 295 Intel Xeon E5-2620 CPU cores.

  1 citation

• An algorithm to parallelise parton showers on a GPU

  Michael H. Seymour, Siddharth Sule

  SciPost Phys. Codebases 33 (2024) · published 12 August 2024 | Toggle abstract · pdf

  The Single Instruction, Multiple Thread (SIMT) paradigm of GPU programming does not support the branching nature of a parton shower algorithm by definition. However, modern GPUs are designed to schedule threads with diverging processes independently, allowing them to handle such branches. With regular thread synchronisation and careful treatment of the individual steps, one can simulate a parton shower on a GPU. We present a Sudakov veto algorithm designed to simulate parton branching on multiple events in parallel. We also release a CUDA C++ program that generates matrix elements, showers partons and computes jet rates and event shapes for LEP at 91.2 GeV on a GPU. To benchmark its performance, we also provide a near-identical C++ program designed to simulate events serially on a CPU. While the consequences of branching are not absent, we demonstrate that a GPU can provide the throughput of a many-core CPU. As an example, we show that the time taken to shower $10^6$ events on one NVIDIA TESLA V100 GPU is equivalent to that of 295 Intel Xeon E5-2620 CPU cores.

  1 citation

• TorchGDM: A GPU-accelerated Python toolkit for multi-scale electromagnetic scattering with automatic differentiation

  Sofia Ponomareva, Adelin Patoux, Clément Majorel, Antoine Azéma, Aurélien Cuche, Christian Girard, Arnaud Arbouet, Peter R. Wiecha

  SciPost Phys. Codebases 60 (2025) · published 22 October 2025 | Toggle abstract · pdf

  We present "torchGDM", a numerical framework for nano-optical simulations based on the Green's Dyadic Method (GDM). This toolkit combines a hybrid approach, allowing for both fully discretized nano-structures and structures approximated by sets of effective electric and magnetic dipoles. It supports simulations in three dimensions and for infinitely long, two-dimensional structures. This capability is particularly suited for multi-scale modeling, enabling accurate near-field calculations within or around a discretized structure embedded in a complex environment of scatterers represented by effective models. Importantly, torchGDM is entirely implemented in PyTorch, a well-optimized and GPU-enabled automatic differentiation framework. This allows for the efficient calculation of exact derivatives of any simulated observable with respect to various inputs, including positions, wavelengths or permittivity, but also intermediate parameters like Green's tensor components, which can be interesting for physics informed deep learning applications. We anticipate that this toolkit will be valuable for applications merging nano-photonics and machine learning, as well as for solving nano-photonic optimization and inverse problems, such as the global design and characterization of metasurfaces, where optical interactions between structures are critical.

  1 citation

• On the Virasoro fusion and modular kernels at any irrational central charge

  Julien Roussillon

  SciPost Phys. 17, 138 (2024) · published 20 November 2024 | Toggle abstract · pdf

  We propose a series representation for the Virasoro fusion and modular kernels at any irrational central charge. Two distinct, yet closely related formulas are needed for the cases $c∈ \mathbb{C} \backslash (-∞,1]$ and $c <1$. Our proposal for $c <1$ agrees numerically with the fusion transformation of the four-point spherical conformal blocks, whereas our proposal for $c∈ \mathbb{C} \backslash (-∞,1]$ agrees numerically with Ponsot and Teschner's integral formula for the fusion kernel. The case of the modular kernel is studied as a special case of the fusion kernel.

  1 citation

• Spin-lattice couplings and effect of displacements on magnetic interactions of a skyrmion system PdFe/Ir(111)

  Banasree Sadhukhan, Anders Bergman, Johan Hellsvik, Patrik Thunström, Anna Delin

  SciPost Phys. 18, 064 (2025) · published 21 February 2025 | Toggle abstract · pdf

  PdFe/Ir(111) has attracted tremendous attention for next-generation spintronics devices due to existence of magnetic skyrmions with the external magnetic field. Our density functional theoretical calculations in combination with spin dynamics simulation suggest that the spin spiral phase in fcc stacked PdFe/Ir(111) flips into the skyrmion lattice phase around B$_{ext} \sim$ 6 T. This leads to the microscopic understanding of the thermodynamic and kinetic behaviours affected by the intrinsic spin-lattice couplings (SLCs) in this skyrmion material for magneto-mechanical properties. Here we calculate fully relativistic SLC parameters from first principle computations and investigate the effect of SLC on dynamical magnetic interactions in skyrmion multilayers PdFe/Ir(111). The exchange interactions arising from next nearest-neighbors (NN) in this material are highly frustrated and responsible for enhancing skyrmion stability. We report the larger spin-lattice effect on both dynamical Heisenberg exchanges and Dzyaloshinskii-Moriya interactions for next NN compared to NN which is in contrast with recently observed spin-lattice effect in bulk bcc Fe and CrI$_3$ monolayer. Based on our analysis, we find that the effective measures of SLCs in fcc (hcp) stacking of PdFe/Ir(111) are $\sim 2.71 ( \sim 2.36)$ and $\sim 14.71 ( \sim21.89)$ times stronger for NN and next NN respectively, compared to bcc Fe. The linear regime of displacement for SLC parameters is $≤$ 0.02 Å which is 0.72% of the lattice constant for PdFe/Ir(111). The microscopic understanding of SLCs provided by our current study could help in designing spintronic devices based on thermodynamic properties of skyrmion multilayers.

  1 citation

• Dynamics of systems with varying number of particles: From Liouville equations to general master equations for open systems

  Mauricio J. del Razo, Luigi Delle Site

  SciPost Phys. 18, 001 (2025) · published 6 January 2025 | Toggle abstract · pdf

  A varying number of particles is one of the most relevant characteristics of systems of interest in nature and technology, ranging from the exchange of energy and matter with the surrounding environment to the change of particle number through internal dynamics such as reactions. The physico-mathematical modeling of these systems is extremely challenging, with the major difficulty being the time dependence of the number of degrees of freedom and the additional constraint that the increment or reduction of the number and species of particles must not violate basic physical laws. Theoretical models, in such a case, represent the key tool for the design of computational strategies for numerical studies that deliver trustful results. In this manuscript, we review complementary physico-mathematical approaches of varying number of particles inspired by rather different specific numerical goals. As a result of the analysis on the underlying common structure of these models, we propose a unifying master equation for general dynamical systems with varying number of particles. This equation embeds all the previous models and can potentially model a much larger range of complex systems, ranging from molecular to social agent-based dynamics.

  1 citation

• Integrable models on Rydberg atom chains

  Luke Corcoran, Marius de Leeuw, Balázs Pozsgay

  SciPost Phys. 18, 139 (2025) · published 25 April 2025 | Toggle abstract · pdf

  We initiate a systematic study of integrable models for spin chains with constrained Hilbert spaces; we focus on spin-1/2 chains with the Rydberg constraint. We extend earlier results for medium-range spin chains to the constrained Hilbert space, and formulate an integrability condition. This enables us to construct new integrable models with fixed interaction ranges. We classify all time- and space-reflection symmetric integrable Rydberg-constrained Hamiltonians of range 3 and 4. At range 3, we find a single family of integrable Hamiltonians: the so-called RSOS quantum chains, which are related to the well-known RSOS models of Andrews, Baxter, and Forrester. At range 4 we find two families of models, the first of which is the constrained XXZ model. We also find a new family of models depending on a single coupling $z$. We provide evidence of two critical points related to the golden ratio $\phi$, at $z=\phi^{-1/2}$ and $z=\phi^{3/2}$. We also perform a partial classification of integrable Hamiltonians for range 5.

  1 citation

• Codebase release 0.2 for JFSD

  Kim William Torre, Raoul D. Schram, Joost de Graaf

  SciPost Phys. Codebases 56-r0.2 (2025) · published 30 June 2025 | Toggle abstract · src

  Stokesian Dynamics (SD) is a powerful computational framework for simulating the motion of particles in a viscous Newtonian fluid under Stokes-flow conditions. Traditional SD implementations can be computationally expensive as they rely on the inversion of large mobility matrices to determine hydrodynamic interactions. Recently, however, the simulation of thermalized systems with large numbers of particles has become feasible [Fiore and Swan, J. Fluid Mech. 878, 544 (2019)]. Their "fast Stokesian dynamics" (FSD) method leverages a saddle-point formulation to ensure overall scaling of the algorithm that is linear in the number of particles $\mathcal{O}(N)$; performance relies on dedicated graphics-processing-unit computing. Here, we present a different route toward implementing FSD, which instead leverages the Just-in-Time (JIT) compilation capabilities of Google JAX. We refer to this implementation as JFSD and perform benchmarks on it to verify that it has the right scaling and is sufficiently fast by the standards of modern computational physics. In addition, we provide a series of physical test cases that help ensure accuracy and robustness, as the code undergoes further development. Thus, JFSD is ready to facilitate the study of hydrodynamic effects in particle suspensions across the domains of soft, active, and granular matter.

  1 citation