Mascon gravity + spherical harmonics
Point-mass superposition for arbitrary geometries; fully-normalized harmonic expansion as a cross-check. JAX mirror for end-to-end differentiability.
Mission 01 — In development
Mathilde recovers the interior density distribution of asteroids and comet nuclei from satellite-to-satellite tracking — the same physics that mapped the Earth's gravity field with GRACE, scaled down to a single small body and a constellation of cubesats.
The problem
We can image an asteroid's surface from a flyby. We can sometimes infer its bulk mass from a tracking arc. But the interior — whether it is a coherent monolith, a rubble pile bound by friction, or a differentiated body with a denser core — remains a guess.
That guess is everything. It determines whether a body can be safely deflected, whether a sample-return mission can anchor to it, whether prospective resources are concentrated or diffuse, and how the body formed in the first place.
Today, the only way to peek inside a small body is to send a flagship mission costing north of a billion dollars and wait a decade. That doesn't scale to a Solar System with hundreds of thousands of catalogued near-Earth and main-belt objects.
The approach
The principle is older than the satellites that carry it: a constellation in orbit is itself a gravity instrument. As the body's mass distribution bends the orbits, the inter-spacecraft ranges record the signal. The inversion problem — turning ranges back into densities — is what we do.
Three or more cubesats co-orbit the body in the rotating frame. Each pair exchanges inter-satellite range measurements continuously.
Sub-millimeter range and range-rate signals encode every inhomogeneity in the body's gravity field — exactly what GRACE and GRAIL exploited.
We discretize the body into mass concentrations. The forward model is a linear sum of point-mass potentials over the constellation's orbits.
A JAX-differentiable forward model feeds a NUTS sampler. Out comes a posterior over densities — uncertainty included, not just a point estimate.
Capabilities
Point-mass superposition for arbitrary geometries; fully-normalized harmonic expansion as a cross-check. JAX mirror for end-to-end differentiability.
RK45 integration in the body-fixed frame with Coriolis and centrifugal terms. Energy and Jacobi-integral diagnostics for trustworthy long arcs.
Torque-free Euler equations and quaternion kinematics — the body's spin state is part of the model, not an idealization.
Single-link SST and dual one-way ranging (DOWR) with imperfect clocks: bias, drift, white frequency noise.
Extended Kalman filter that recovers spacecraft trajectories and clock parameters from inter-satellite ranges alone.
Hamiltonian Monte Carlo on the JAX forward model produces a calibrated posterior over interior density — not a single fit.
Visualizations
The figures below come straight from the simulation pipeline that produces our results. No marketing renders.
Generated live in your browser using the same spherical-harmonic perturbation model the inversion uses.
A J₂+J₃ gravity field fingerprints itself in the range derivative. Real signal, simulated.
Filter accuracy as a function of constellation size. More baselines, better state recovery.
Posterior draws over interior density from NUTS. Where the data does and doesn't constrain the body.
// Live figures rendered from examples/ coming soon
Founders
Mathilde is engineered, not pitched. Every line of the inversion pipeline comes from the same hand.
Meet the foundersLogbook
Every design decision in Mathilde started as a discussion. We keep them — equations, dead ends, derivations and all — in the open.
All entriesCurrently raising
If asteroid prospecting, planetary defense, or in-orbit science is on your thesis, we'd like to share the deck and walk you through the model.