Publications

Order by: publication date number of citations

• Robust independent validation of experiment and theory: Rivet version 4 release note

  Christian Bierlich, Andy Buckley, Jonathan Butterworth, Christian Gutschow, Leif Lonnblad, Tomasz Procter, Peter Richardson, Yoran Yeh

  SciPost Phys. Codebases 36 (2024) · published 4 November 2024 | Toggle abstract · pdf

  The Rivet toolkit is the primary mechanism for phenomenological preservation of collider-physics measurements, containing both a computational core and API for analysis implementation, and a large collection of more than a thousand preserved analyses. In this note we summarise the main changes in the new Rivet 4 major release series. These include a major generalisation and more semantically coherent model for histograms and related data objects, a thorough clean-up of inelegant and legacy observable-computation tools, and new systems for extended analysis-data, incorporation of preserved machine-learning models, and serialization for high-performance computing applications. Where these changes introduce backward-incompatible interface changes, existing analyses have been updated and indications are given on how to update new analysis routines and workflows.

  8 citations

• HEJ 2.2: W boson pairs and Higgs boson plus jet production at high energies

  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 (2023) · published 6 December 2023 | Toggle abstract · pdf

  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.

  7 citations

• Codebase release 2.0 for ALF (Algorithms for Lattice Fermions)

  F. F. Assaad, M. Bercx, F. Goth, A. Götz, J. S. Hofmann, E. Huffman, Z. Liu, F. Parisen Toldin, J. S. E. Portela, J. Schwab

  SciPost Phys. Codebases 1-r2.0 (2022) · published 22 August 2022 | Toggle abstract · src

  The Algorithms for Lattice Fermions package provides a general code for the finite-temperature and projective auxiliary-field quantum Monte Carlo algorithm. The code is engineered to be able to simulate any model that can be written in terms of sums of single-body operators, of squares of single-body operators and single-body operators coupled to a bosonic field with given dynamics. The package includes five pre-defined model classes: SU(N) Kondo, SU(N) Hubbard, SU(N) t-V and SU(N) models with long range Coulomb repulsion on honeycomb, square and N-leg lattices, as well as $Z_2$ unconstrained lattice gauge theories coupled to fermionic and $Z_2$ matter. An implementation of the stochastic Maximum Entropy method is also provided. One can download the code from our Git instance at https://git.physik.uni-wuerzburg.de/ALF/ALF/-/tree/ALF-2.0 and sign in to file issues.

  7 citations

• Codebase release 3.4 for NetKet

  Filippo Vicentini, Damian Hofmann, Attila Szabó, Dian Wu, Christopher Roth, Clemens Giuliani, Gabriel Pescia, Jannes Nys, Vladimir Vargas-Calderón, Nikita Astrakhantsev, Giuseppe Carleo

  SciPost Phys. Codebases 7-r3.4 (2022) · published 24 August 2022 | Toggle abstract · src

  We introduce version 3 of NetKet, the machine learning toolbox for many-body quantum physics. NetKet is built around neural-network quantum states and provides efficient algorithms for their evaluation and optimization. This new version is built on top of JAX, a differentiable programming and accelerated linear algebra framework for the Python programming language. The most significant new feature is the possibility to define arbitrary neural network ansätze in pure Python code using the concise notation of machine-learning frameworks, which allows for just-in-time compilation as well as the implicit generation of gradients thanks to automatic differentiation. NetKet 3 also comes with support for GPU and TPU accelerators, advanced support for discrete symmetry groups, chunking to scale up to thousands of degrees of freedom, drivers for quantum dynamics applications, and improved modularity, allowing users to use only parts of the toolbox as a foundation for their own code.

  7 citations

• LinReTraCe: The linear response transport centre

  Matthias Pickem, Emanuele Maggio, Jan M. Tomczak

  SciPost Phys. Codebases 16 (2023) · published 17 July 2023 | Toggle abstract · pdf

  We describe the "Linear Response Transport Centre" (LinReTraCe), a package for the simulation of transport properties of solids. LinReTraCe captures quantum (in)coherence effects beyond semi-classical Boltzmann techniques, while incurring similar numerical costs. The enabling algorithmic innovation is a semi-analytical evaluation of Kubo formulae for resistivities and the coefficients of Hall, Seebeck and Nernst. We detail the program's architecture, its interface and usage with electronic-structure packages such as WIEN2k, VASP, and Wannier90, as well as versatile tight-binding settings.

  7 citations

• Quantum many-body simulations with PauliStrings.jl

  Nicolas Loizeau, J. Clayton Peacock, Dries Sels

  SciPost Phys. Codebases 54 (2025) · published 24 March 2025 | Toggle abstract · pdf

  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.

  7 citations

• DFT2kp: Effective kp models from ab-initio data

  João Victor V. Cassiano, Augusto de Lelis Araújo, Paulo E. Faria Junior, Gerson J. Ferreira

  SciPost Phys. Codebases 25 (2024) · published 5 February 2024 | Toggle abstract · pdf

  The $k\cdot p$ method, combined with group theory, is an efficient approach to obtain the low energy effective Hamiltonians of crystalline materials. Although the Hamiltonian coefficients are written as matrix elements of the generalized momentum operator $\pi= p+ p_{SOC}$ (including spin-orbit coupling corrections), their numerical values must be determined from outside sources, such as experiments or ab initio methods. Here, we develop a code to explicitly calculate the Kane (linear in crystal momentum) and Luttinger (quadratic in crystal momentum) parameters of $k\cdot p$ effective Hamiltonians directly from ab initio wavefunctions provided by Quantum ESPRESSO. Additionally, the code analyzes the symmetry transformations of the wavefunctions to optimize the final Hamiltonian. This is an optional step in the code, where it numerically finds the unitary transformation $U$ that rotates the basis towards an optimal symmetry-adapted representation informed by the user. Throughout the paper, we present the methodology in detail and illustrate the capabilities of the code applying it to a selection of relevant materials. Particularly, we show a "hands-on" example of how to run the code for graphene (with and without spin-orbit coupling). The code is open source and available at https://gitlab.com/dft2kp/dft2kp.

  6 citations

• kdotpy: k·p theory on a lattice for simulating semiconductor band structures

  Wouter Beugeling, Florian Bayer, Christian Berger, Jan Böttcher, Leonid Bovkun, Christopher Fuchs, Maximilian Hofer, Saquib Shamim, Moritz Siebert, Li-Xian Wang, Ewelina M. Hankiewicz, Tobias Kießling, Hartmut Buhmann, Laurens W. Molenkamp

  SciPost Phys. Codebases 47 (2025) · published 28 January 2025 | Toggle abstract · pdf

  The software project kdotpy provides a Python application for simulating electronic band structures of semiconductor devices with k·p theory on a lattice. The application implements the widely used Kane model, capable of reliable predictions of transport and optical properties for a large variety of topological and non-topological materials with a zincblende crystal structure. The application automates the tedious steps of simulating band structures. The user inputs the relevant physical parameters on the command line, for example materials and dimensions of the device, magnetic field, and temperature. The program constructs the appropriate matrix Hamiltonian on a discretized lattice of spatial coordinates and diagonalizes it. The physical observables are extracted from the eigenvalues and eigenvectors and saved as output. The program is highly customizable with a large set of configuration options and material parameters.

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

  6 citations

• Nevanlinna.jl: A Julia implementation of Nevanlinna analytic continuation

  Kosuke Nogaki, Jiani Fei, Emanuel Gull, Hiroshi Shinaoka

  SciPost Phys. Codebases 19 (2023) · published 13 November 2023 | Toggle abstract · pdf

  We introduce a Julia implementation of the recently proposed Nevanlinna analytic continuation method. The method is based on Nevanlinna interpolants and inherently preserves the causality of a response function due to its construction. For theoretical calculations without statistical noise, this continuation method is a powerful tool to extract real-frequency information from numerical input data on the Matsubara axis. This method has been applied to first-principles calculations of correlated materials. This paper presents its efficient and full-featured open-source implementation of the method including the Hamburger moment problem and smoothing.

  6 citations