Harvesting primordial black holes from stochastic trees with FOREST

We introduce a novel framework to implement stochastic inflation on stochastic trees, modelling the inflationary expansion as a branching process. Combined with the δ N formalism, this allows us to generate real-space maps of the curvature perturbation that fully capture quantum diffusion and its non-perturbative backreaction during inflation. Unlike lattice methods, trees do not proceed on a fixed background since new spacetime units emerge dynamically as trees unfold, naturally incorporating metric fluctuations. The recursive structure of stochastic trees also offers remarkable numerical efficiency, and we develop the FOrtran Recursive Exploration of Stochastic Trees (FOREST) tool and demonstrate its performance. We show how primordial black holes blossom at unbalanced nodes of the trees, and how their mass distribution can be obtained while automatically accounting for the “cloud-in-cloud” effect. In the “quantum-well” toy model, we find broad mass distributions, with mild power laws terminated by exponential tails. We finally compare our results with existing approximations in the literature and discuss several prospects.