darkmatter 0.0.2

Dark matter simulation engine — gravitational fields, particle dynamics, halo stability, and cosmological constants
Documentation
# Changelog

All notable changes to the **DarkMatter** crate will be documented in this file.

Format based on [Keep a Changelog](https://keepachangelog.com/).

---

## [0.0.1] — 2026-04-05

### Added

**Core framework** — 12 modules, 69 source files, 781+ public items

- `candidates` — WIMP annihilation and scattering, axion mass and coupling, fuzzy DM soliton profiles and BEC condensation, warm DM transfer functions, asymmetric DM, primordial black hole constraints
- `clusters` — dynamical mass estimators, intracluster medium, Sunyaev-Zel'dovich effect, cluster mass functions and scaling relations
- `constants` — physical constants (G, c, ħ, k_B, solar mass, parsec) and cosmological parameters (H₀, Ω_m, Ω_Λ, T_CMB)
- `cosmology` — Boltzmann equation solver (RK4) with freeze-out, freeze-in, coannihilation, entropy injection; CMB anisotropies; relic abundance; perturbation theory; structure formation; reionization; 21 cm signal
- `detection` — direct detection recoil rates, indirect detection fluxes (gamma, neutrino, antimatter), collider signatures (mono-X), EFT operators, dark photon mediators, nuclear response functions, gravitational wave signatures (phase transitions, binary dephasing, LISA sensitivity), annihilation/decay spectra (bb̄, τ⁺τ⁻, W⁺W⁻, μ⁺μ⁻)
- `dynamics` — Barnes-Hut octree O(N log N) gravity with leapfrog integrator, direct N-body, Chandrasekhar dynamical friction, tidal stripping, virial and Jeans stability
- `galaxy` — galaxy formation in DM halos, merger rates, rotation curve decomposition
- `halos` — NFW, Einasto, Burkert, isothermal density profiles; adiabatic contraction; density spikes around black holes; subhalo mass functions; halo rotation curves
- `large_scale` — two-point correlation functions, power spectra, cosmic web filaments, void statistics
- `physics` — gravitational potentials, lensing (convergence, shear, magnification), SIDM transfer cross-sections, Sommerfeld enhancement (Coulomb, Yukawa, resonant), phase-space distributions (Fermi-Dirac, Bose-Einstein, Tremaine-Gunn bound), thermodynamics
- `solar_system` — DM capture by solar system bodies, local density measurements, orbital perturbations
- `utils` — unit conversions, special functions (erf, Bessel, NFW), ODE solvers, root finders, statistical distributions

**Tests** — 87 integration tests across 7 test files

**Examples** — 7 runnable examples: fuzzy_dm, sommerfeld, spectra, phase_space, boltzmann, gravitational_waves, barnes_hut

### What is solid

- Full coverage of major dark matter candidate models
- Production-ready halo profiles and rotation curve fitting
- Boltzmann solver supports freeze-out, freeze-in, coannihilation, and entropy injection
- Sommerfeld enhancement covers Coulomb, Yukawa, and resonant regimes with unitarity bound
- Barnes-Hut tree gravity scales to large N
- All 87 tests pass, clippy clean with zero warnings

### Points to watch

- First release — API may evolve
- No parallelism yet in N-body routines
- Spectra limited to four annihilation channels

### Verdict

Solid foundation. Ready for research use and community feedback.