To study the behaviour of wave dark matter, I developed qaxions is a high-performance simulation code. It solves the Schrödinger-Poisson system of equations on a 2D or 3D lattice, capturing the quantum wave-like dynamics that make axions and similar particles distinct from heavier dark matter models. The code is written in C++ and designed to scale from a laptop to HPC clusters, with OpenMP multithreading and optional GPU acceleration via CUDA that delivers up to 40× speedups over CPU-only runs. Python bindings via pybind11 make it easy to set up simulations, process outputs, and interface with existing analysis workflows without leaving Python. qaxions is intended as a flexible and accessible entry point for researchers looking to explore ultralight dark matter dynamics on modern hardware.

Check out the online documentation and the source code

Code features

• High-performance C++ core
• OpenMP parallelism
• HDF5 and FFTW support
• GPU acceleration with CUDA and cuFFT
• Python bindings via pybind11
• Linux and macOS support