Publications

Order by: publication date number of citations

Codebase release 0.2 for Pychastic

Radost Waszkiewicz, Maciej Bartczak, Kamil Kolasa, Maciej Lisicki

SciPost Phys. Codebases 11-r0.2 (2023) · published 3 April 2023 | · src

In the last decade, Python-powered physics simulations ecosystem has been growing steadily, allowing greater interoperability, and becoming an important tool in numerical exploration of physical phenomena, particularly in soft matter systems. Driven by the need for fast and precise numerical integration in colloidal dynamics, here we formulate the problem of Brownian Dynamics (BD) in a mathematically consistent formalism of the Itō calculus, and develop a Python package to assist numerical computations. We show that, thanks to the automatic differentiation packages, the classical truncated Taylor-Itō integrators can be implemented without the burden of computing the derivatives of the coefficient functions beforehand. Furthermore, we show how to circumvent the difficulties of BD simulations such as calculations of the divergence of the mobility tensor in the diffusion equation and discontinuous trajectories encountered when working with dynamics on $S^2$ and $SO(3)$. The resulting Python package, Pychastic, is capable of performing BD simulations including hydrodynamic interactions at speeds comparable to dedicated implementations in lower-level programming languages, but with a much simpler end-user interface.

1 citation
MatsubaraFunctions.jl: An equilibrium Green's function library in the Julia programming language

Dominik Kiese, Anxiang Ge, Nepomuk Ritz, Jan von Delft, Nils Wentzell

SciPost Phys. Codebases 24 (2024) · published 11 January 2024 | · pdf

The Matsubara Green's function formalism stands as a powerful technique for computing the thermodynamic characteristics of interacting quantum many-particle systems at finite temperatures. In this manuscript, our focus centers on introducing MatsubaraFunctions.jl, a Julia library that implements data structures for generalized n-point Green's functions on Matsubara frequency grids. The package's architecture prioritizes user-friendliness without compromising the development of efficient solvers for quantum field theories in equilibrium. Following a comprehensive introduction of the fundamental types, we delve into a thorough examination of key facets of the interface. This encompasses avenues for accessing Green's functions, techniques for extrapolation and interpolation, as well as the incorporation of symmetries and a variety of parallelization strategies. Examples of increasing complexity serve to demonstrate the practical utility of the library, supplemented by discussions on strategies for sidestepping impediments to optimal performance.

1 citation
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 | · 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.

1 citation
Codebase release 2.2 for Pyqcm

Théo N. Dionne, Alexandre Foley, Moïse Rousseau, David Sénéchal

SciPost Phys. Codebases 23-r2.2 (2023) · published 20 December 2023 | · src

Pyqcm is a Python/C++ library that implements a few quantum cluster methods with an exact diagonalization impurity solver. Quantum cluster methods are used in the study of strongly correlated electrons to provide an approximate solution to Hubbard-like models. The methods covered by this library are Cluster Perturbation Theory (CPT), the Variational Cluster Approach (VCA) and Cellular (or Cluster) Dynamical Mean Field Theory (CDMFT). The impurity solver (the technique used to compute the cluster's interacting Green function) is exact diagonalization from sparse matrices, using the Lanczos algorithm and variants thereof. The core library is written in C++ for performance, but the interface is in Python, for ease of use and inter-operability with the numerical Python ecosystem. The library is distributed under the GPL license.

1 citation
Pyqcm: An open-source Python library for quantum cluster methods

Théo N. Dionne, Alexandre Foley, Moïse Rousseau, David Sénéchal

SciPost Phys. Codebases 23 (2023) · published 20 December 2023 | · pdf

Pyqcm is a Python/C++ library that implements a few quantum cluster methods with an exact diagonalization impurity solver. Quantum cluster methods are used in the study of strongly correlated electrons to provide an approximate solution to Hubbard-like models. The methods covered by this library are Cluster Perturbation Theory (CPT), the Variational Cluster Approach (VCA) and Cellular (or Cluster) Dynamical Mean Field Theory (CDMFT). The impurity solver (the technique used to compute the cluster's interacting Green function) is exact diagonalization from sparse matrices, using the Lanczos algorithm and variants thereof. The core library is written in C++ for performance, but the interface is in Python, for ease of use and inter-operability with the numerical Python ecosystem. The library is distributed under the GPL license.

1 citation
Codebase release 1.1 for TASEPy

Luca Ciandrini, Richmond L. Crisostomo, Juraj Szavits-Nossan

SciPost Phys. Codebases 22-r1.1 (2023) · published 13 December 2023 | · src

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic lattice model for one-dimensional particle transport subject to excluded-volume interactions. Solving the inhomogeneous TASEP in which particles' hopping rate vary across the lattice is a long-standing problem. In recent years, a power series approximation (PSA) has been developed to tackle this problem, however no computer algorithm currently exists that implements this approximation. This paper addresses this issue by providing a Python-based package TASEPy that finds the steady state solution of the inhomogeneous TASEP for any set of hopping rates using the PSA truncated at a user-defined order.

1 citation
TASEPy: A Python-based package to iteratively solve the inhomogeneous exclusion process

Luca Ciandrini, Richmond L. Crisostomo, Juraj Szavits-Nossan

SciPost Phys. Codebases 22 (2023) · published 13 December 2023 | · pdf

The totally asymmetric simple exclusion process (TASEP) is a paradigmatic lattice model for one-dimensional particle transport subject to excluded-volume interactions. Solving the inhomogeneous TASEP in which particles' hopping rate vary across the lattice is a long-standing problem. In recent years, a power series approximation (PSA) has been developed to tackle this problem, however no computer algorithm currently exists that implements this approximation. This paper addresses this issue by providing a Python-based package TASEPy that finds the steady state solution of the inhomogeneous TASEP for any set of hopping rates using the PSA truncated at a user-defined order.

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 | · 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
xSPDE3: Extensible software for stochastic ordinary and partial differential equations

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

SciPost Phys. Codebases 17 (2023) · published 2 October 2023 | · pdf

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
GRIFFIN: A C++ library for electroweak radiative corrections in fermion scattering and decay processes

Lisong Chen, Ayres Freitas

SciPost Phys. Codebases 18 (2023) · published 10 October 2023 | · pdf

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

Previous 1 2 3 4 5 6 Next

Publications

1 citation

1 citation

1 citation

1 citation

1 citation

1 citation

1 citation

1 citation

1 citation

1 citation