Rimu.jl: Random integrators for many-body quantum systems
arXiv:2601.19505
Abstract
Rimujl is a Julia package for solving many-body quantum problems. The core of the package is a matrix-free implementation of Hamiltonians and other operators and compact representation of Fock states, which together allow for efficient methods suitable for high-performance computing. Rimujl includes a Julia implementation of the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm which is a type of projector QMC algorithm for stochastically solving the time-independent Schrödinger equation. It also includes many well-known model Hamiltonians, and an interface for exact diagonalisation based on external eigenvalue solvers. Both the stochastic and exact diagonalisation methods are accessed with a CommonSolvejl interface. We describe the FCIQMC algorithm and how to obtain estimators of observables as well as the key features of the implementation.
39 pages, 10 figures