---
title: Generation III Bern Model Overview
url: https://www.emergentmind.com/topics/generation-3-bern-model
type: topic
---

# Generation III Bern Model Overview

Searching arXiv for papers on the Generation III Bern Model and related NGPPS studies.
The Generation III Bern model is a combined global end-to-end planetary formation and evolution model within the New Generation Planetary Population Synthesis (NGPPS) framework. Built on the core accretion paradigm, it follows planetary systems from a protoplanetary gas and planetesimal disc through concurrent multi-embryo growth, gas accretion, migration, collisions, and full N-body interactions, and then into Gigayear evolution including cooling, contraction, atmospheric escape, bloating, Roche-lobe overflow, and stellar tides. Its stated purpose is to predict not only masses and orbital distances but also radii, luminosities, magnitudes, evaporation rates, compositions, and long-term orbital evolution across a broad mass and separation range, and to do so in a form suitable for direct statistical comparison with exoplanet surveys [2007.05561]. In later NGPPS studies, the same framework is used to generate synthetic ensembles of planetary systems and to compare them with Kepler-like transit samples and the HARPS/Coralie radial-velocity survey after applying explicit observational selection effects [2105.12745, 2509.09762].

## 1. Genealogy within the Bern population-synthesis program

The Generation III Bern model is presented as the first Bern framework that unifies multi-planet formation and long-term post-formation evolution in a single coupled architecture. Earlier Bern generations treated only subsets of that problem: Generation I addressed single-planet formation up to disc dispersal; Generation Ib extended to long-term thermodynamic evolution but remained essentially single-planet; Generation II introduced multiple embryos and N-body interactions but only during the formation phase; Generation III combines these strands into one end-to-end model [2007.05561].

| Bern generation | Scope | Limitation relative to Generation III |
|---|---|---|
| Generation I | Single-planet formation up to disc dispersal | No long-term evolution; no concurrent multi-planet N-body formation |
| Generation Ib | Long-term thermodynamic evolution | Still essentially single-planet |
| Generation II | Multiple embryos and N-body interactions during formation | No full post-formation evolution |
| Generation III | Multi-planet formation plus long-term evolution | Integrates both in one framework |

This generational framing is important because the defining novelty of Generation III is not a single physical prescription, but the coupling of previously separate components. The model is explicitly described as bringing together N-body interactions among many concurrently forming embryos and full post-formation evolution including cooling, contraction, atmospheric escape, bloating, Roche-lobe overflow, and tidal decay [2007.05561]. A plausible implication is that the model’s identity lies as much in systems integration as in any individual submodule.

## 2. Physical architecture and governing components

The model is organized as a chain of coupled submodules for stellar evolution, gas-disc evolution, planetesimal-disc evolution, planetary internal structure, accretion, migration, N-body dynamics, and long-term evolution [2007.05561]. The gas disc is treated as a 1D axisymmetric viscous disc with \(\alpha\)-viscosity, internal and external photoevaporation, irradiation, and vertical thermal structure; later NGPPS studies summarize the adopted turbulence parameter as fixed to
\[
\alpha = 2\times 10^{-3},
\]
chosen to reproduce observed disc mass-accretion relations [2509.09762]. The 2020 global-model paper gives the viscous diffusion equation in the form
\[
\frac{\partial \Sigma_g}{\partial t} = \frac{1}{r}\frac{\partial}{\partial r} \left[ 3r^{1/2}\frac{\partial}{\partial r} \left(r^{1/2}\nu\Sigma_g\right) \right] - \text{photoevaporation} - \text{planet accretion},
\qquad
\nu = \alpha c_s H,
\]
and emphasizes that disc thermodynamics includes both viscous heating and stellar irradiation [2007.05561].

Solid growth proceeds through planetesimal accretion in the oligarchic regime. Planetesimals are treated as a surface-density field with dynamical state variables \(e\) and \(i\), rather than as individual bodies, and their evolution includes gas drag and stirring by embryos and other planetesimals [2007.05561]. In the nominal Generation III setup described in the original NGPPS paper, the planetesimal radius is fixed to \(300\,\mathrm{m}\), a simplification later identified as consequential when confronting RV data [2007.05561, 2509.09762]. The model tracks planetary internal structure through 1D hydrostatic equations and treats planets as having an “onion-like” structure with iron core, silicate mantle, water-ice layer, and H/He envelope; a later improvement noted in the metallicity study is that core composition is tracked self-consistently using the disc compositional model, so the iron/silicate ratio can evolve through planetesimal accretion and giant impacts [2507.09874].

Gas accretion is split into an attached, planet-limited phase and a detached, disc-limited phase. The papers consistently identify the transition to runaway gas accretion with core masses of order \(\sim 10\,M_\oplus\), after which accretion becomes limited by disc supply [2105.12745, 2509.09762]. Migration is included in both Type I and Type II regimes, with transition after the Crida et al. criterion in later NGPPS descriptions; multi-planet interactions are handled with the Mercury N-body package, allowing scattering, resonances, collisions, and dynamical instabilities [2509.09762, 2507.09874]. Long-term evolution then extends to \(10\,\mathrm{Gyr}\), including cooling and contraction, atmospheric escape, and tidal migration, which is why Generation III can be used not only for masses and periods but also for radius-distribution and radius-valley analyses [2007.05561, 2507.09874].

## 3. Population-synthesis workflow and initial conditions

The Generation III Bern model is used not merely to simulate individual systems, but to perform population synthesis by Monte Carlo sampling of initial disc and embryo properties. In the NGPPS architecture papers, the synthesis starts from a gas disc, a planetesimal disc, and multiple embryos already inserted into the disc; a representative setup places embryos of mass \(10^{-2}\,M_\oplus\) randomly in logarithmic orbital distance between the inner edge and \(40\,\mathrm{au}\), with no two embryos within \(10\) Hill radii [2007.05561]. The metallicity study states that the synthesis starts from a disk state after collapse in the Class I stage, with \(100\) embryos already inserted, each with mass \(10^{-2}\,M_\oplus\), randomly placed in logarithmic orbital distance [2507.09874]. In the “peas in a pod” study, embryo locations are described as distributed approximately uniformly in \(\log a\) between the disk inner edge and \(40\,\mathrm{au}\), with spacing by at least \(10H\) after runaway growth [2105.12745].

The disk properties varied across synthetic systems include gas disk mass, disk size, disk inner edge, disk lifetime or photoevaporation rate, and metallicity or dust-to-gas ratio [2105.12745, 2507.09874]. The metallicity prescription is given explicitly as
\[
\frac{f_{D/G}}{f_{D/G,\oplus}} = 10^{[\mathrm{Fe/H}]},
\]
with the solar dust-to-gas ratio taken as \(f_{D/G,\oplus}=0.0149\) [2507.09874]. This makes host-star metallicity a direct control variable for the initial solid reservoir, and therefore for core-growth timescales, gas-accretion onset, migration history, and dynamical excitation.

Different NGPPS studies use related but distinct synthetic ensembles. The “peas in a pod” analysis simulates \(1{,}000\) planetary systems for each of three nominal embryo counts—\(20\), \(50\), and \(100\) embryos per disk—and focuses mainly on the \(100\)-embryo population; it also removes “failed embryos” below \(0.1\,M_\oplus\) from most architectural analyses [2105.12745]. The HARPS/Coralie comparison focuses mainly on the nominal population NG76, with \(100\) embryos per disc and a \(20\,\mathrm{Myr}\) formation stage, and on NG76longshot, which extends the formation and N-body stage to \(100\,\mathrm{Myr}\) to test later dynamical instabilities [2509.09762]. Synthetic systems are then evolved individually to \(10\,\mathrm{Gyr}\), although specific observational comparisons may evaluate them at \(5\,\mathrm{Gyr}\) or use survey-specific age cuts [2509.09762, 2507.09874].

## 4. Forward modeling of observations and survey confrontation

A central feature of the Generation III Bern model in practice is that its output is not compared to observations “raw.” Instead, the synthetic populations are passed through explicit observational filters. For transit surveys, the NGPPS VI paper introduces KOBE, “Kepler Observes Bern Exoplanets,” as a forward model that takes the intrinsic Bern populations and simulates what a Kepler-like survey would detect [2105.12745]. KOBE operates in three sequential layers: `kobeshadows/` for transit geometry, `kobetransits/` for transit detectability and Kepler signal-to-noise using the FGK-star CDPP distribution, and `kobevetter/` for Robovetter completeness and reliability. A planet becomes a periodic threshold crossing event if it has at least two transits and \(S/N \ge 7.1\), following Weiss et al. in using a minimum of two transits rather than three [2105.12745]. The final output is the KOBE catalogue, from which the authors define KOBE multi-planetary systems (KMPS) for direct comparison with the California-Kepler Survey multi-planet sample.

For radial-velocity data, the HARPS/Coralie study applies an analogous biasing procedure to the synthetic populations. The observed reference is the combined HARPS + Coralie volume-limited sample of \(822\) FGK stars within \(50\,\mathrm{pc}\), containing \(169\) planets/candidates around \(105\) systems and a mean multiplicity
\[
\bar{N}_\mathrm{pl} = 1.61.
\]
For each Monte Carlo realization, \(822\) synthetic systems are drawn, a random line of sight is assigned, \(M\sin i\) is computed, and each planet receives a detection probability from the survey completeness map in \(M\sin i\) and period; planets in regions with detection probability below \(1\%\) are treated as undetected [2509.09762]. Statistical comparison then uses a two-sided Kolmogorov-Smirnov test for one-dimensional distributions, a 2D KS test for mass-period comparisons, and \(1000\) biased realizations of the synthetic population to estimate confidence intervals [2509.09762].

These forward models are not merely technical wrappers. They are structurally important because they allow three-way comparisons among the underlying Bern population, the biased synthetic detections, and real surveys. This makes it possible to distinguish intrinsic formation signatures from the distortions introduced by transit geometry, completeness, vetting, RV sensitivity, and inclination effects [2105.12745, 2509.09762].

## 5. Emergent architectures, demographic regularities, and metallicity systematics

One of the most studied outputs of the Generation III Bern model is the emergence of within-system architectural regularities. In the KOBE analysis of “peas in a pod,” adjacent planets in the same synthetic system exhibit strong size similarity and mass similarity. In the KOBE-observed multi-planet sample, the size correlation is reported as Pearson \(R \approx 0.64\), while the underlying Bern population already shows \(R \approx 0.66\); about \(65\%\) of adjacent pairs are size-ordered in KMPS, whereas only about \(48\%\) are ordered in the underlying population [2105.12745]. The model interpretation is that size similarity is largely derivative of mass similarity, because the KMPS planets display strong mass-radius monotonicity with a Spearman coefficient of about \(0.81\) [2105.12745]. The same study argues that mass similarity is primordial, arising from oligarchic growth and broadly uniform early accretion, whereas spacing and packing are not primordial but are strengthened later by orbital migration, resonance capture, mergers, ejections, and other dynamical interactions [2105.12745].

The spacing and packing results are particularly revealing about the model’s internal dynamics. For spacing similarity, measured through
\[
\frac{(P_{j+2}/P_{j+1})}{(P_{j+1}/P_j)} \approx 1,
\]
the underlying Bern population shows a stronger correlation, about \(R \approx 0.55\), than the KOBE-observed population, about \(R \approx 0.25\), while the CKS sample is reported at about \(R \approx 0.46\) [2105.12745]. Packing correlations show the same pattern: for KMPS, size-based packing gives \(R \approx 0.23\) and mass-based packing \(R \approx 0.26\), whereas the underlying population gives about \(R \approx 0.45\) and \(R \approx 0.57\), respectively [2105.12745]. The stated interpretation is that transit selection tends to dilute, rather than generate, these correlations. Systems with larger lost-planet fractions, used as a proxy for dynamical activity, show stronger spacing and packing correlations, especially in the \(50\)- and \(100\)-embryo populations [2105.12745].

The metallicity study extends this demographic program by analyzing how Generation III populations depend on \([\mathrm{Fe/H}]\). Using biased synthetic samples, it finds positive metallicity correlations for hot, warm, and cold Jupiters with \(\beta \approx 1.3\) in the occurrence-law parameterization
\[
F([\mathrm{Fe/H}],t)=C\times 10^{\beta [\mathrm{Fe/H}]},
\]
and weaker but still positive correlations for Neptune-size planets, with \(\beta \sim 0.5\) for hot Neptunes, \(\beta \sim 0.8\) for warm Neptunes, and \(\beta \sim 0.6\) for cold Neptunes [2507.09874]. For small Kepler-like planets of \(1\)–\(3.5\,R_\oplus\), the occurrence rate, system fraction, and multiplicity increase with metallicity up to \([\mathrm{Fe/H}] \sim 0.1\) and then decrease at higher metallicity, while sub-Earths exhibit a strong anti-correlation with metallicity [2507.09874].

A major success claimed for the Generation III Bern model is the reproduction of radius-valley morphology and its metallicity dependence. The metallicity paper adopts the standard solar-mass valley boundaries
\[
R^{\rm Valley}_{\rm upper} = 2.1\,R_\oplus\left(\frac{P}{10\,\mathrm{days}}\right)^{-0.09},
\qquad
R^{\rm Valley}_{\rm lower} = 1.7\,R_\oplus\left(\frac{P}{10\,\mathrm{days}}\right)^{-0.09},
\]
and reports that the radius valley becomes deeper and more prominent with increasing \([\mathrm{Fe/H}]\), the super-Earth/sub-Neptune ratio changes accordingly, and the mean radius above the valley increases with metallicity, while the lower-side mean radius remains essentially constant [2507.09874]. The quantitative fits
\[
C_{\rm valley}= 15.9^{+3.4}_{-4.6}\,[\mathrm{Fe/H}] + 7.0^{+1.6}_{-1.1},
\]
\[
A_{\rm valley}= -1.0^{+0.2}_{-0.2}\,[\mathrm{Fe/H}] + 0.12^{+0.01}_{-0.01},
\]
\[
\log_{10}\left(\frac{R^{+}_{\rm valley}}{R_\oplus}\right) = 0.06^{+0.02}_{-0.02}\,[\mathrm{Fe/H}] + 0.45^{+0.03}_{-0.02},
\]
\[
R^{-}_{\rm valley}=1.35^{+0.02}_{-0.02}R_\oplus,
\]
\[
f_{\rm NP}=0.15^{+0.04}_{-0.04}\,[\mathrm{Fe/H}] + 0.05^{+0.02}_{-0.01}
\]
are reported to be in quantitative agreement with the LAMOST-Gaia-Kepler analysis of Chen et al. (2022) [2507.09874].

## 6. Empirical performance, tensions, and current limitations

When confronted quantitatively with RV data, the Generation III Bern model reproduces several broad observed features but also shows systematic mismatches. The HARPS/Coralie comparison states that the nominal population reproduces two main groups in the mass-distance plane—close-in sub-Neptunes and distant giants—a bimodal mass function with a relative desert between Neptune and Saturn masses, mean multiplicity close to the observed value, a broadly correct metallicity trend, a broadly similar period-ratio distribution, and an eccentricity distribution that at least qualitatively resembles the survey sample [2509.09762]. The same paper, however, is unusually explicit that the nominal synthetic population has too many planets by about \(70\%\), a desert that is too deep by about \(60\%\), a relative excess of giant planets by about \(40\%\), eccentricities that are on average too low by a factor of about two, with median \(0.07\) versus \(0.15\), and a metallicity effect that is too weak; the synthetic planets are also overall too close to the star compared with the HARPS sample [2509.09762].

Several discrepancies are highlighted as structurally important. The nominal population predicts \(290^{+30}_{-28}\) detectable planets compared with \(169\) observed, but with mean multiplicity \(1.50^{+0.31}_{-0.25}\), close to the observed \(1.61\), implying that the excess lies mainly in the number of systems with detectable planets rather than in extreme over-multiplicity [2509.09762]. The giant-planet peak is too massive, near \(\sim 1000\,M_\oplus\) instead of the observed \(\sim 500\,M_\oplus\), which the paper interprets as evidence that envelope masses and detached-phase disc-limited gas accretion are overestimated [2509.09762]. Sub-Neptunes are too light and too close to the star, with observed systems typically at \(3\)–\(100\) days while synthetic ones are more often at \(1\)–\(30\) days [2509.09762]. The nominal synthetic population contains no hot Jupiters, even though the survey should contain about \(5\)–\(10\) of them given an occurrence rate of about \(0.5\)–\(1\%\) in \(1000\) stars [2509.09762].

Parameter retuning improves some observables but does not remove the underlying tension. The optimized NG192 population increases the planetesimal radius from \(300\,\mathrm{m}\) to \(2\,\mathrm{km}\), reduces both Type I and Type II migration to \(7/8\) of the nominal rate, and changes the maximum gas accretion rate from the minimum of Bondi and local reservoir rate to the minimum of Bondi rate and the radial gas flow through the disc [2509.09762]. This yields \(179\) planets in \(120\) systems with mean multiplicity \(1.49^{+0.37}_{-0.29}\), and the KS mass-function distance becomes \(1.35\), just below the threshold \(1.36\) [2509.09762]. Yet the model still fails in joint mass-period space, with \(p<0.001\) in the 2D KS test, and the planets remain too close in [2509.09762]. The paper therefore concludes that simple parameter adjustments are insufficient, and that physical processes appear to be missing.

The metallicity study reaches a closely related conclusion from a different angle. Although it reports that the Generation III Bern model successfully reproduces many observed metallicity correlations, it also states that the predicted dependences of orbital eccentricity and period on \([\mathrm{Fe/H}]\) are significantly weaker than observed [2507.09874]. For instance, the model finds a metallicity difference between planets with \(P<10\) d and \(P>10\) d of
\[
\Delta[\mathrm{Fe/H}] = 0.05^{+0.01}_{-0.01}\,\mathrm{dex},
\]
whereas the observational difference is described as around \(0.15\) dex [2507.09874]. The stated missing ingredients include the absence of a metallicity dependence of the disk inner edge, the absence of post-disc high-eccentricity migration, and the fact that only the first \(100\,\mathrm{Myr}\) of planetary interactions are modeled, so long-term secular chaos and late dynamical delivery are not captured [2507.09874]. Across both survey-comparison papers, the emerging interpretation is consistent: the Generation III Bern model captures the core physics of formation and early evolution in a physically rich and statistically useful way, but remains incomplete for processes that set the dynamical temperature of planetary systems over Gyr timescales and for channels that produce wider initial architectures, slower detached-phase gas accretion, stronger eccentricity excitation, and late-arriving close-in giants [2509.09762, 2507.09874].

Source: https://www.emergentmind.com/topics/generation-3-bern-model