Cellular Automation of Galactic Habitable Zone

We present a preliminary results of our Galactic Habitable Zone (GHZ) 2D probabilistic cellular automata models. The relevant time-scales (emergence of life, it’s diversification and evolution influenced with the global risk function) are modeled as the probability matrix elements and are chosen in accordance with the Copernican principle to be well-represented by the data inferred from the Earth’s fossil record. With Fermi’s paradox as a main boundary condition the resulting histories of astrobiological landscape are discussed.

💡 Research Summary

The paper introduces a two‑dimensional probabilistic cellular automaton (CA) framework to model the temporal and spatial evolution of the Galactic Habitable Zone (GHZ). The authors discretize the Milky Way’s disk into a 500 × 500 grid, each cell representing a small region of the galaxy. Every cell can occupy one of four biological states: no life, unicellular life, multicellular life, and intelligent (technological) life. Transitions between these states are governed by a probability matrix P(i → j) whose elements are derived from three main ingredients: (1) the characteristic timescales for each evolutionary step, (2) the local metallicity (Z) that determines whether the chemistry needed for life is feasible, and (3) a “global risk function” R(t, x, y) that captures the suppressive effect of astrophysical hazards such as supernovae, gamma‑ray bursts, and heightened activity near the galactic centre.

The characteristic timescales are calibrated using Earth’s fossil record under the Copernican principle, i.e., the Earth is assumed to be a typical example. The authors adopt approximate mean durations of 1 Gyr for the emergence of unicellular organisms, 2 Gyr for the transition to multicellular complexity, and 3 Gyr for the appearance of intelligence. These values are not simply inverted to obtain probabilities; instead, they are combined with the risk function to produce a time‑dependent transition rate that varies across the galactic plane. Metallicities are taken from observed radial gradients, and only cells with Z ≥ 0.1 Z⊙ are allowed to host life. The risk function is modeled as an exponential attenuation of transition probabilities, with higher hazard levels producing stronger suppression.

Initial conditions are set by seeding a very low density of primitive life (≈10⁻⁶ of the total stellar mass) across the disk, weighted by the present‑day star formation rate (SFR) profile. The simulation proceeds in 10 Myr steps for a total of 10 Gyr, updating each cell’s state according to the locally computed transition probabilities.

Key results emerge from the ensemble of runs. First, the “sweet spot” for life – where metallicity is sufficient and astrophysical hazards are modest – lies in an annulus roughly 8–12 kpc from the galactic centre. Within this region, unicellular and multicellular life proliferate relatively robustly. Intelligent life, however, is far more constrained: it tends to appear only in a narrower band (≈9–11 kpc) that avoids both the high‑risk inner galaxy and the low‑metallicity outer disk. Second, when the authors impose the Fermi paradox as a hard boundary condition – namely, that the number of observable technological civilizations in the Milky Way must be ≤ 1 – they must dramatically lower the probability of the final transition (multicellular → intelligent). In practice, the model forces the intelligent‑life transition probability down to the order of 10⁻⁸ per cell per Gyr, rendering intelligent civilizations extremely rare on galactic scales.

The discussion acknowledges several limitations. The two‑dimensional representation neglects vertical structure and stellar orbital migration, which could redistribute habitable zones over time. The probability matrix is calibrated solely on Earth data, so alternative evolutionary pathways that might be common elsewhere are not captured. The risk function includes only known astrophysical threats and ignores potential unknown or stochastic events. Despite these simplifications, the CA approach provides a computationally inexpensive way to explore a high‑dimensional parameter space and to visualise how different assumptions shape the astrobiological landscape.

Future work outlined by the authors includes extending the model to three dimensions, incorporating up‑to‑date metallicity and SFR maps from surveys such as Gaia and APOGEE, testing alternative life‑emergence pathways (e.g., panspermia, radiation‑driven chemistry), and refining the Fermi‑paradox constraint by modelling signal propagation, detection biases, and observational selection effects.

In conclusion, the study demonstrates that a probabilistic cellular automaton can serve as a powerful tool for quantifying the distribution and timing of life and intelligence within the GHZ. By explicitly embedding the Copernican principle and the Fermi paradox into the model, the authors provide a coherent, quantitative framework that supports the view that intelligent, technologically capable civilizations are likely to be exceedingly sparse across the Milky Way.