Publications

Order by: publication date number of citations

• Codebase release 1.1 for QDsim

  Valentina Gualtieri, Charles Renshaw-Whitman, Vinicius Hernandes, Eliska Greplova

  SciPost Phys. Codebases 46-r1.1 (2025) · published 23 January 2025 | Toggle abstract · src

  We introduce QDsim, a Python package tailored for the rapid generation of charge stability diagrams in large-scale quantum dot devices, extending beyond traditional double or triple dots. QDsim is founded on the constant interaction model from which we rephrase the task of finding the lowest energy charge configuration as a convex optimization problem. Therefore, we can leverage the existing package CVXPY, in combination with an appropriate powerful solver, for the convex optimization which streamlines the creation of stability diagrams and polytopes. Through multiple examples, we demonstrate how QDsim enables the generation of large-scale dataset that can serve a basis for the training of machine-learning models for automated tuning algorithms. While the package currently does not support quantum effects beyond the constant interaction model, QDsim is a tool that directly addresses the critical need for cost-effective and expeditious data acquisition for better tuning algorithms in order to accelerate the development of semiconductor quantum devices.

  1 citation

• Codebase release 0.2 for Carlo.jl

  Lukas Weber

  SciPost Phys. Codebases 49-r0.2 (2025) · published 7 February 2025 | Toggle abstract · src

  Carlo.jl is a Monte Carlo simulation framework written in Julia. It provides MPI-parallel scheduling (including parallel tempering), organized storage of input, checkpoint, and output files, as well as statistical postprocessing. With a minimalist design, it aims to aid the development of high-quality Monte Carlo codes, especially for demanding applications in condensed matter and statistical physics. This hands-on user guide shows how to implement a simple code with Carlo.jl, how to run an existing stochastic series expansion code based on Carlo.jl, and provides benchmarks to show its efficacy.

  1 citation

• Kolya: An open-source package for inclusive semileptonic B decays

  Matteo Fael, Ilija S. Milutin, K. Keri Vos

  SciPost Phys. Codebases 55 (2025) · published 25 March 2025 | Toggle abstract · pdf

  We introduce the code kolya, an open-source tool for phenomenological analyses of inclusive semileptonic $B\to X_c l\bar \nu_l$ decays. It contains a library to compute predictions for the total rate and various kinematic moments within the framework of the heavy quark expansion, utilizing the so-called kinetic scheme. The library currently includes power corrections up to $1/m_b^5$. All available QCD perturbative corrections are implemented via interpolation grids for fast numerical evaluation. We also include effects from new physics parameterized as Wilson coefficients of dimension-six operators in the weak effective theory below the electroweak scale. The library is interfaced to CRunDec for easy evaluation of the quark masses and strong coupling constant at different renormalization scales. The library is developed in Python and does not require compilation. The code is flexible and can be extended to other kinds of inclusive decays such as $B$ meson lifetimes, $B \to X_u l\bar \nu_l$ and $D$ meson decays.

  1 citation

• QDsim: A user-friendly toolbox for simulating large-scale quantum dot devices

  Valentina Gualtieri, Charles Renshaw-Whitman, Vinicius Hernandes, Eliska Greplova

  SciPost Phys. Codebases 46 (2025) · published 23 January 2025 | Toggle abstract · pdf

  We introduce QDsim, a Python package tailored for the rapid generation of charge stability diagrams in large-scale quantum dot devices, extending beyond traditional double or triple dots. QDsim is founded on the constant interaction model from which we rephrase the task of finding the lowest energy charge configuration as a convex optimization problem. Therefore, we can leverage the existing package CVXPY, in combination with an appropriate powerful solver, for the convex optimization which streamlines the creation of stability diagrams and polytopes. Through multiple examples, we demonstrate how QDsim enables the generation of large-scale dataset that can serve a basis for the training of machine-learning models for automated tuning algorithms. While the package currently does not support quantum effects beyond the constant interaction model, QDsim is a tool that directly addresses the critical need for cost-effective and expeditious data acquisition for better tuning algorithms in order to accelerate the development of semiconductor quantum devices.

  1 citation

• Codebase release 0.5 for HarmonicBalance.jl

  Jan Košata, Javier del Pino, Toni L. Heugel, Oded Zilberberg

  SciPost Phys. Codebases 6-r0.5 (2022) · published 24 August 2022 | Toggle abstract · src

  HarmonicBalance.jl is a publicly available Julia package designed to simplify and solve systems of periodic time-dependent nonlinear ordinary differential equations. Time dependence of the system parameters is treated with the harmonic balance method, which approximates the system's behaviour as a set of harmonic terms with slowly-varying amplitudes. Under this approximation, the set of all possible steady-state responses follows from the solution of a polynomial system. In HarmonicBalance.jl, we combine harmonic balance with contemporary implementations of symbolic algebra and the homotopy continuation method to numerically determine all steady-state solutions and their associated fluctuation dynamics. For the exploration of involved steady-state topologies, we provide a simple graphical user interface, allowing for arbitrary solution observables and phase diagrams. HarmonicBalance.jl is a free software available at https://github.com/NonlinearOscillations/HarmonicBalance.jl.

  1 citation

• Codebase release 5.3 for EDIpack

  Lorenzo Crippa, Igor Krivenko, Samuele Giuli, Gabriele Bellomia, Alexander Kowalski, Francesco Petocchi, Alberto Scazzola, Markus Wallerberger, Giacomo Mazza, Luca de Medici, Giorgio Sangiovanni, Massimo Capone, Adriano Amaricci

  SciPost Phys. Codebases 58-r5.3 (2025) · published 24 September 2025 | Toggle abstract · src

  We present a next-generation version of EDIpack, a flexible, high-performance numerical library using Lanczos-based exact diagonalization to solve generic quantum impurity problems, such as those introduced in Dynamical Mean-Field Theory to describe extended strongly correlated materials. This new release efficiently solves impurity problems allowing for different broken-symmetry solutions, including superconductivity, featuring local spin-orbit coupling and/or electron-phonon coupling. It provides quick access to dynamical correlation functions on the entire complex frequency plane at zero and low-temperatures. The modular architecture of the software not only provides Fortran APIs but also includes bindings to C/C++, interfaces with Python and Julia or with TRIQS and w2dynamics research platforms, thus ensuring unprecedented level of inter-operability. The outlook includes further extensions to study quantum materials and cold atoms quantum simulators, as well as quantum information applications.

  1 citation

• The ITransverse.jl library for transverse tensor network contractions

  Stefano Carignano

  SciPost Phys. Codebases 63 (2026) · published 9 January 2026 | Toggle abstract · pdf

  Transverse contraction methods are extremely promising tools for the efficient contraction of tensor networks associated with the time evolution of quantum many-body systems, allowing in some cases to circumvent the entanglement barrier that would normally prevent the study of quantum dynamics with classical resources. We present here the ITransverse.jl package, written in Julia and based on ITensors.jl, containing several of these high-level algorithms, including novel prescriptions for efficient truncations of temporal matrix product states.

  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

• Codebase release 1.0 for GRIFFIN

  Lisong Chen, Ayres Freitas

  SciPost Phys. Codebases 18-r1.0 (2023) · published 10 October 2023 | Toggle abstract · src

  This paper describes a modular framework for the description of electroweak scattering and decay processes, including but not limited to Z-resonance physics. The framework consistently combines a complex-pole expansion near an s-channel resonance with a regular fixed-order perturbative description away from the resonance in a manifestly gauge-invariant scheme. Leading vertex correction contributions are encapsulated in form factors that can be predicted or treated as numerical fit parameters. This framework has been implemented in the publicly available object-oriented C++ library GRIFFIN. Version 1.0 of this library provides Standard Model predictions for the IR-subtracted matrix elements for the process $f\bar{f} \to f'\bar{f}'$ with full NNLO and leading higher-order contributions on the Z-resonance, and with NLO corrections off-resonance. The library can straightforwardly be extended to include higher-order corrections, should they become available, or predictions for new physics models. It can be interfaced with Monte-Carlo programs to account for QED and QCD initial-state and final-state radiation.

  1 citation

• Codebase release 1.0 for NumDynInflation

  Siddharth S. Bhatt, Swagat S. Mishra, Soumen Basak, Surya N. Sahoo

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

  Numerical simulations of the inflationary dynamics are presented here for a single canonical scalar field minimally coupled to gravity. We spell out the basic equations governing the inflationary dynamics in terms of cosmic time $t$ and define a set of dimensionless variables convenient for numerical analysis. We then provide a link to our simple numerical Python code on GitHub that can be used to simulate the background dynamics as well as the evolution of linear perturbations during inflation. The code computes both scalar and tensor power spectra for a given inflaton potential $V(\phi)$. We discuss a concrete algorithm to use the code for various purposes, especially for computing the enhanced scalar power spectrum in the context of Primordial Black Holes and scalar-induced Gravitational Waves. We also compare the efficiency of different variables used in the literature to compute the scalar fluctuations. We intend to extend the framework to simulate the dynamics of a number of different quantities, including the computation of scalar-induced second-order tensor power spectrum in the near future.

  1 citation