Publications

Order by: publication date number of citations

• A Python GPU-accelerated solver for the Gross-Pitaevskii equation and applications to many-body cavity QED

  Lorenzo Fioroni, Luca Gravina, Justyna Stefaniak, Alexander Baumgärtner, Fabian Finger, Davide Dreon, Tobias Donner

  SciPost Phys. Codebases 38 (2024) · published 6 November 2024 | Toggle abstract · pdf

  TorchGPE is a general-purpose Python package developed for solving the Gross-Pitaevskii equation (GPE). This solver is designed to integrate wave functions across a spectrum of linear and non-linear potentials. A distinctive aspect of TorchGPE is its modular approach, which allows the incorporation of arbitrary self-consistent and time-dependent potentials, e.g., those relevant in many-body cavity QED models. The package employs a symmetric split-step Fourier propagation method, effective in both real and imaginary time. In our work, we demonstrate a significant improvement in computational efficiency by leveraging GPU computing capabilities. With the integration of the latter technology, TorchGPE achieves a substantial speed-up with respect to conventional CPU-based methods, greatly expanding the scope and potential of research in this field.

  1 citation

• Codebase release 3.44 for xSPDE

  Simon Kiesewetter, Ria R. Joseph, Peter D. Drummond

  SciPost Phys. Codebases 17-r3.44 (2023) · published 2 October 2023 | Toggle abstract · src

  The xSPDE toolbox treats stochastic partial and ordinary differential equations, with applications in biology, chemistry, engineering, medicine, physics and quantum technologies. It computes statistical averages, including time-step and/or sampling error estimation. xSPDE can provide higher order convergence, Fourier spectra and probability densities. The toolbox has graphical output and $\chi^{2}$ statistics, as well as weighted, projected, or forward-backward equations. It can generate input-output quantum spectra. All equations may have independent periodic, Dirichlet, and Neumann or Robin boundary conditions in any dimension, for any vector field component, and at either end of any interval.

  1 citation

• Codebase release 2.0 for sauce

  Caleb Marshall

  SciPost Phys. Codebases 37-r2.0 (2024) · published 6 November 2024 | Toggle abstract · src

  Low energy nuclear physics experiments are transitioning towards fully digital data acquisition systems. Realizing the gains in flexibility afforded by these systems relies on equally flexible data reduction techniques. In this paper, methods utilizing data frames and in-memory techniques to work with data, including data from self-triggering, digital data acquisition systems, are discussed within the context of a Python package, sauce. It is shown that data frame operations can encompass common analysis needs and allow interactive data analysis. Two event building techniques, dubbed referenced and referenceless event building, are shown to provide a means to transform raw list mode data into correlated multi-detector events. These techniques are demonstrated in the analysis of two example data sets.

  1 citation

• Data reduction for low energy nuclear physics experiments using data frames

  Caleb Marshall

  SciPost Phys. Codebases 37 (2024) · published 6 November 2024 | Toggle abstract · pdf

  Low energy nuclear physics experiments are transitioning towards fully digital data acquisition systems. Realizing the gains in flexibility afforded by these systems relies on equally flexible data reduction techniques. In this paper, methods utilizing data frames and in-memory techniques to work with data, including data from self-triggering, digital data acquisition systems, are discussed within the context of a Python package, sauce. It is shown that data frame operations can encompass common analysis needs and allow interactive data analysis. Two event building techniques, dubbed referenced and referenceless event building, are shown to provide a means to transform raw list mode data into correlated multi-detector events. These techniques are demonstrated in the analysis of two example data sets.

  1 citation

• Codebase release 0.1 for infstat

  Kyoungchul Kong, Konstantin T. Matchev, Stephen Mrenna, Prasanth Shyamsundar

  SciPost Phys. Codebases 14-r0.1 (2023) · published 7 July 2023 | Toggle abstract · src

  We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential $\varphi$. In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

  1 citation

• Codebase release r1.0 for amcheck

  Andriy Smolyanyuk, Libor Šmejkal, Igor I. Mazin

  SciPost Phys. Codebases 30-r1.0 (2024) · published 12 June 2024 | Toggle abstract · src

  Altermagnets (AM) is a recently discovered class of collinear magnets that share some properties (anomalous transport, etc) with ferromagnets, some (zero net magnetization) with antiferromagnets, while also exhibiting unique properties (spin-splitting of electronic bands and resulting spin-splitter current). Since the moment compensation in AM is driven by symmetry, it must be possible to identify them by analyzing the crystal structure directly, without computing the electronic structure. Given the significant potential of AM for spintronics, it is very useful to have a tool for such an analysis. This work presents an open-access code implementing such a direct check.

  1 citation

• Codebase release 2.2 for HEJ

  Jeppe R. Andersen, Bertrand Ducloué, Conor Elrick, Hitham Hassan, Andreas Maier, Graeme Nail, Jérémy Paltrinieri, Andreas Papaefstathiou, Jennifer M. Smillie

  SciPost Phys. Codebases 21-r2.2 (2023) · published 6 December 2023 | Toggle abstract · src

  We present version 2.2 of the High Energy Jets (HEJ) Monte Carlo event generator for hadronic scattering processes at high energies. The new version adds support for two further processes of central phenomenological interest, namely the production of a W boson pair with equal charge together with two or more jets and the production of a Higgs boson with at least one jet. Furthermore, a new prediction for charged lepton pair production with high jet multiplicities is provided in the high-energy limit. The accuracy of HEJ 2.2 can be increased further through an enhanced interface to standard predictions based on conventional perturbation theory. We describe all improvements and provide extensive usage examples. HEJ 2.2 can be obtained from https://hej.hepforge.org.

  1 citation

• Codebase release r0.1 for PanScales

  Melissa van Beekveld, Mrinal Dasgupta, Basem Kamal El-Menoufi, Silvia Ferrario Ravasio, Keith Hamilton, Jack Helliwell, Alexander Karlberg, Rok Medves, Pier Francesco Monni, Gavin P. Salam, Ludovic Scyboz, Alba Soto-Ontoso, Gregory Soyez, Rob Verheyen

  SciPost Phys. Codebases 31-r0.1 (2024) · published 25 June 2024 | Toggle abstract · src

  In this article, we document version 0.1 of the PANSCALES code for parton shower simulations. With the help of a few examples, we discuss basic usage of the code, including tests of logarithmic accuracy of parton showers. We expose some of the numerical techniques underlying the logarithmic tests and include a description of how users can implement their own showers within the framework. Some of the simpler logarithmic tests can be performed in a few minutes on a modern laptop. As an early step towards phenomenology, we also outline some aspects of a preliminary interface to PYTHIA 8.3, for access to its hard matrix elements and its hadronisation modules.

  1 citation

• Codebase release 1.5 for PauliStrings.jl

  Nicolas Loizeau, J. Clayton Peacock, Dries Sels

  SciPost Phys. Codebases 54-r1.5 (2025) · published 24 March 2025 | Toggle abstract · src

  We present the Julia package PauliStrings.jl for quantum many-body simulations, which performs fast operations on the Pauli group by encoding Pauli strings in binary. All of the Pauli string algebra is encoded into low-level logic operations on integers, and is made efficient by various truncation methods which allow for systematic extrapolation of the results. We illustrate the effectiveness of our package by (i) performing Heisenberg time evolution through direct numerical integration and (ii) by constructing a Liouvillian Krylov space. We benchmark the results against tensor network methods, and we find our package performs favorably. In addition, we show that this representation allows for easy encoding of any geometry. We present results for chaotic and integrable spin systems in 1D as well as some examples in 2D. Currently, the main limitations are the inefficiency of representing non-trivial pure states (or other low-rank operators), as well as the need to introduce dissipation to probe long-time dynamics.

  1 citation

• Codebase release 1.2 for YASTN

  Marek M. Rams, Gabriela Wójtowicz, Aritra Sinha, Juraj Hasik

  SciPost Phys. Codebases 52-r1.2 (2025) · published 26 February 2025 | Toggle abstract · src

  We present an open-source tensor network Python library for quantum many-body simulations. At its core is an Abelian-symmetric tensor, implemented as a sparse block structure managed by a logical layer on top of a dense multidimensional array backend. This serves as the basis for higher-level tensor network algorithms operating on matrix product states and projected entangled pair states. An appropriate backend, such as PyTorch, gives direct access to automatic differentiation (AD) for cost-function gradient calculations and execution on GPU and other supported accelerators. We show the library performance in simulations with infinite projected entangled-pair states, such as finding the ground states with AD and simulating thermal states of the Hubbard model via imaginary time evolution. For these challenging examples, we identify and quantify sources of the numerical advantage exploited by the symmetric-tensor implementation.

  1 citation