JuliaCon 2026

Testability-First Design for Few-Body Systems Physics
2026-08-12 , Room 5

Open-source software (OSS) is a practical way to ensure reproducibility and avoid reinventing the wheel in computational physics. In particular, few-body systems physics has various targets and various software requirements, so we are building an open-source ecosystem in Julia as a shared infrastructure with reusability and extensibility. In this talk, we introduce the architectural design and development roadmap derived from system-level testability and report the current development status.


Few-body systems physics is an intersection of many subfields of physics and chemistry. Its targets are various few-body systems across various scales, from hadrons and nuclei to atoms, molecules, and quantum dots. Ideally, researchers would like to calculate physical observables for arbitrary particle numbers, arbitrary interactions, and arbitrary quantum states. In practice, no single software package can satisfy all of these requirements, and researchers often develop self-tailored codebases for their specific problems.

To make this development culture scalable, we are building an open-source ecosystem as shared infrastructure that provides reusable building blocks and a flexible common interface in the Julia programming language. CI/CD with automated testing is a baseline requirement for maintainability and for accepting external contributions. Beyond its role in maintenance, testing also has a role in design, because the requirements for the system-level testability derive the architectural design.

Differential testing based on cross-method comparisons has been used in few-body calculations. For established benchmark problems, the final outputs can be validated using test oracles built from a database of results in the literature. However, for previously unexplored problems, there is no verified “correct answer”. This requirement motivates a common interface that allows multiple methods to operate on the same Hamiltonian.

The Rayleigh–Ritz method using explicitly correlated Gaussian (ECG) basis functions is widely used for solving the Schrödinger equation of quantum-mechanical few-body systems, while variational Monte Carlo (VMC) can calculate the energy for a given trial wave function. Currently, we are developing FewBodyECG.jl and MetropolisAlgorithm.jl (toward FewBodyVMC.jl) for cross-checking multiple methods. FewBodyHamiltonian.jl supports flexible construction of Hamiltonians as a common input for multiple methods. FewBodyDB.jl works as a test oracle database, and TwoBody.jl works as a prototype of FewBody.jl, the common interface of this JuliaFewBody ecosystem.

While we don't necessarily practice test-first or test-driven development, we find testability to be a compass for our overall architectural design and long-term development roadmap. We hope this ecosystem will open the gate to a "trading port" of computational methods.

I am a Ph.D. candidate at Yokohama City University and a junior research associate (JRA) at RIKEN in Japan. I launched JuliaFewBody for developing FewBody.jl, general-purpose flexible solvers for quantum mechanical few-body problems. At JuliaCon 2025, I presented several of his packages towards completing FewBody.jl.

I am the organizer of several events related to Julia and computational physics:
- Julia in Physics 2024,
- JuliaLang Japan 2025,
- Spring School on Computational Physics,
- Computational Physics Hackathon 2026 (CompPhysHack2026).

PhD student at the University of Copenhagen working on tensor networks