Embedded Branching Particle Process

• Embedded branching particle process is a methodological framework that integrates branching dynamics within external geometries or reduced state spaces.
• It encompasses various models—spatial trees, censored processes, measure-valued populations, and crossing trees—to capture phenomena like exclusion, selection slowdown, and density limits.
• These models reveal crossover regimes from explosive to constrained growth, offering insights applicable to evolution, epidemics, and random media analyses.

An embedded branching particle process is a branching system whose analytically relevant structure is obtained by placing branching particles in an external state space or by extracting, from a larger interacting system, a lower-dimensional process that still has an exact or asymptotically exact branching interpretation. In the literature represented here, this includes Euclidean trees with hard-core exclusion, frontier populations in branching-selection systems that evolve as censored Galton–Watson chains, branching random walks in random media represented through parabolic Anderson and spine formulas, measure-valued populations coded by Lévy excursions, and path-valued constructions whose multiscale crossing tree is itself a branching genealogy (Forgerini et al., 2011).

1. Conceptual scope and basic constructions

The term does not denote a single universally fixed model class. Rather, several closely related constructions recur.

A first class consists of genuinely spatial branching particle systems in which every particle carries a position in a metric or Euclidean space and branching depends on that position. A representative example is the discrete-time tree embedded in $\mathbb R^D$ in which a generation-$t$ node at ${\bf x}_i$ attempts to place up to two children at

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$

with ${\bf \Delta}_i$ sampled uniformly in a $D$-dimensional hypercube of side length $2$, and births suppressed if the proposed child lies within Euclidean distance $a$ of forbidden preexisting nodes (Forgerini et al., 2011).

A second class consists of embedded coordinate processes extracted from a larger interacting particle system. In the branching-selection model on $\mathbb Z$, the coordinate

$k\mapsto Y_k^N(k)$

counts particles at the rightmost possible position at time $t$0, and it evolves exactly as a censored Galton–Watson process: $t$1 Here the embedded object is not the full $t$2-particle configuration but a distinguished frontier class whose law closes autonomously under censoring at level $t$3 (Couronné et al., 2011).

A third class appears when a constrained system is represented by an auxiliary branching process with killing. For the $t$4-branching Markov process on $t$5, each birth is followed by deletion of the current leftmost particle, keeping population size fixed. In the large-$t$6 limit, this system is compared to a branching Markov process killed below a deterministic moving boundary $t$7, with survival condition

$t$8

The selected process is then approximated by the law of a single Markov particle conditioned on survival above $t$9, while the proof proceeds through the unconstrained killed branching system (Bérard et al., 2023).

A fourth class is measure-valued. In the nonlocal branching framework, individuals are indexed by a genealogical tree, but the state of the whole population is the empirical measure

${\bf x}_i$0

The embedded branching structure is the genealogical family indexed by the Ulam–Harris tree, while the observable process is the measure-valued evolution. This dual viewpoint underlies many-to-one formulas, eigenfunction tilts, and growth-fragmentation limits (Cloez, 2011).

A fifth class is genealogical rather than particle-dynamic in the usual spatial sense. In the multifractal construction based on a continuous path ${\bf x}_i$1, a level-${\bf x}_i$2 crossing decomposes into finer subcrossings, and the resulting crossing tree becomes a two-type Galton–Watson process with types corresponding to up- and down-crossings. The embedded branching process is therefore the path’s multiscale crossing genealogy rather than a cloud of independently moving particles (Decrouez et al., 2012).

These variants share a common structural idea: branching is either embedded in an external geometry or recovered as an exact or effective subsystem inside a larger stochastic evolution. This suggests that the term is best understood as a methodological category rather than a single canonical model.

2. Spatial embedding, exclusion, and geometry

The clearest directly spatial example is the Euclidean branching tree with exclusion distance ${\bf x}_i$3. The process starts from a root at the origin in ${\bf x}_i$4, time is the generation index, every node lives forever as a vertex of the tree, and only the current frontier reproduces. In the main version, conflict checks are performed against the previous generation and already created contemporaneous nodes; in a stricter variant, each attempted birth is checked against all existing nodes in the tree except the parent (Forgerini et al., 2011).

In the sparse regime, exclusion is negligible and the number of generation-${\bf x}_i$5 nodes follows the unrestricted binary law

${\bf x}_i$6

The paper identifies this as an initial “explosive” regime. At long times, the occupied region has linear size of order ${\bf x}_i$7, while the exclusion distance forces packing on scale ${\bf x}_i$8. This leads heuristically to

${\bf x}_i$9

for the previous-generation-only model. Matching the two asymptotics yields the crossover estimate

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$0

and in ${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$1 simulations give

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$2

consistent with ${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$3 (Forgerini et al., 2011).

The stricter all-existing-nodes exclusion rule changes the asymptotics because cumulative occupancy, rather than only frontier crowding, becomes the bottleneck. There the packing argument applies to total population,

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$4

and hence

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$5

In ${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$6, this implies a stationary-rate frontier,

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$7

with numerical plateau

${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$8

The contrast between the two variants is therefore structural: short-memory suppression allows continuing frontier growth, whereas all-history exclusion produces much stronger cumulative constraint (Forgerini et al., 2011).

The same paper connects these growth laws to explicit geometry. In one dimension, the previous-generation-only model exhibits an approximately symmetric step profile whose boundaries move outward with speed about ${\bf x}_{\text{new}}={\bf x}_i+{\bf \Delta}_i,$9, with bulk density about ${\bf \Delta}_i$0. In the all-existing-nodes variant, the expanding step profile applies to the distribution of all nodes, with boundary speed about ${\bf \Delta}_i$1 and bulk density about ${\bf \Delta}_i$2. This front-propagation picture explains why the packing heuristics are accurate.

When the embedding space is bounded, finite-size effects induce fluctuation-dominated regimes and extinction. On ${\bf \Delta}_i$3, the population fluctuates around ${\bf \Delta}_i$4, but extinction occurs when an entire generation produces no offspring. For ${\bf \Delta}_i$5, ${\bf \Delta}_i$6, and ${\bf \Delta}_i$7, a typical realization dies out before ${\bf \Delta}_i$8 generations, and over observation horizon ${\bf \Delta}_i$9 the numerical extinction boundary $D$0 is monotone increasing, with $D$1 for small $D$2 and divergence as $D$3 (Forgerini et al., 2011).

The network-theoretic implication is a time-driven transition from small-world to large-world structure. Because the graph is a tree built by generations, diameter $D$4 is of order age $D$5. During exponential growth,

$D$6

whereas during crowding-limited power-law growth, diameter becomes a power-law function of size. The same embedded tree therefore changes its effective network architecture over time without any external parameter change (Forgerini et al., 2011).

3. Embedded processes in selection, censoring, and pruning

A central use of embedded branching particle processes is the reduction of a complicated selected system to an exact or asymptotically exact branching subdynamics.

In the discrete-time branching-selection system on $D$7, the full state $D$8 contains $D$9 particles after each selection step. Each particle at position $2$0 is replaced independently by $2$1 children at $2$2 and $2$3 children at $2$4, after which only the $2$5 rightmost particles are kept. The asymptotic front speed is

$2$6

The key embedded process is the top-class population $2$7, which counts particles that have advanced right at every generation. Its law is exactly that of the censored Galton–Watson chain

$2$8

This identity in law is the analytical core of the model (Couronné et al., 2011).

The censored process has extinction time

$2$9

and the main asymptotic theorem states

$a$0

where $a$1 is the extinction probability of the ordinary Galton–Watson process with offspring law $a$2. The corresponding front slowdown satisfies

$a$3

Thus the large-scale speed defect of the selected frontier is controlled by the lifetime of the embedded censored process (Couronné et al., 2011).

The mechanism can be made quantitative through the last roof-visit time

$a$4

and the two-sided speed bounds

$a$5

The proof constructs comparison systems that are restarted when the embedded top-class population either vanishes or ceases to occupy the full ceiling $a$6. In this sense, the embedded process provides both exact reduction and renewal structure (Couronné et al., 2011).

A related but asymptotic embedding appears in the $a$7-branching Markov process on $a$8. There each particle moves according to a Markov process $a$9, branches at rate $\mathbb Z$0, and every branching event deletes the current leftmost particle. The hydrodynamic limit shows that the empirical c.d.f.

$\mathbb Z$1

converges uniformly to

$\mathbb Z$2

where $\mathbb Z$3 and $\mathbb Z$4 satisfies

$\mathbb Z$5

The selected system is therefore asymptotically represented by a single-particle survival-conditioned law, while the proof uses the auxiliary $\mathbb Z$6-BMP: an unconstrained branching Markov process killed at the moving boundary $\mathbb Z$7 (Bérard et al., 2023).

The minimum particle position

$\mathbb Z$8

also converges to the same deterministic boundary. Under assumptions (i)–(v), the paper proves uniform convergence of $\mathbb Z$9 to $k\mapsto Y_k^N(k)$0 on compact time intervals away from $k\mapsto Y_k^N(k)$1, with polynomial tail bounds. A plausible implication is that endogenous leftmost deletion becomes deterministic killing in the large-$k\mapsto Y_k^N(k)$2 limit (Bérard et al., 2023).

Another embedded reduction appears in semipushed fronts. In the one-dimensional dyadic branching Brownian motion with space-dependent branching rate, drift $k\mapsto Y_k^N(k)$3, and killing at $k\mapsto Y_k^N(k)$4, the decisive rare events are particles that hit a high barrier $k\mapsto Y_k^N(k)$5. The future contribution of such particles is encoded in descendant clusters sampled when they first hit $k\mapsto Y_k^N(k)$6, and the rescaled family size

$k\mapsto Y_k^N(k)$7

converges to a random variable $k\mapsto Y_k^N(k)$8 with heavy tail

$k\mapsto Y_k^N(k)$9

This breakout skeleton drives convergence of the rescaled population to an $t$00-stable continuous-state branching process on time scale $t$01 in the semipushed regime $t$02 (Tourniaire, 2021).

4. Random environments, measure-valued systems, and coded branching populations

Embedded branching particle processes also arise when a branching system is represented by analytic objects such as PDEs, expectations, or excursion-coded counting measures.

For branching random walks in random environment, the basic object is a continuous-time branching random walk on a graph $t$03, with motion given by a Markov process $t$04 and site-dependent splitting and killing rates $t$05 and $t$06. The first-moment field

$t$07

satisfies

$t$08

with Feynman–Kac representation

$t$09

Here the parabolic Anderson model is literally the first-moment equation of an underlying branching particle process. In the multitype version, expected population size is represented by a weighted Markov chain on $t$10, and higher moments are represented by finite-splitting skeletons derived through many-to-one and spine techniques (König, 2020).

The survey emphasizes three embedded viewpoints: the underlying BRWRE behind the parabolic Anderson model, the one-particle Feynman–Kac reduction, and the smaller branching skeletons used for higher moments. These are not alternative models but different embedded representations of the same branching system. They are also the basis for intermittency, localization, and optimal-type-cycle asymptotics in random environment (König, 2020).

A more explicit measure-valued embedding appears in the nonlocal branching system with empirical measure

$t$11

The first-moment operator is

$t$12

and, when $t$13 for some positive eigenfunction $t$14, the process

$t$15

is a martingale. The size-biased auxiliary process $t$16, obtained through the $t$17-transform, yields the many-to-one identity

$t$18

The long-time empirical measure is then governed by ergodic properties of $t$19, while large-population scaling produces a deterministic growth-fragmentation equation (Cloez, 2011).

In the Lévy-coded construction, a single excursion of a spectrally one-sided Lévy process already contains a branching genealogy. For the spectrally positive process with negative drift, level $t$20 is interpreted as time, and jumps that upcross $t$21 define alive particles: $t$22 Each atom records residual lifetime $t$23. As the level increases, particles drift deterministically left at unit speed and die at $t$24; new jumps create offspring at positions determined by overshoots. The resulting $t$25 is a measure-valued Borel right Markov process, and its total mass is a Crump–Mode–Jagers process (He et al., 2012).

In a rather different direction, the random-environment may itself be a branching tree. In the electrostatics model on the profinite completion $t$26 of an infinite rooted branching process, particles do not branch; instead they form a Gibbs gas on the boundary of a random genealogical space. The tree yields a random boundary measure

$t$27

and ultrametric

$t$28

This is not a branching particle system in the classical sense, but it is an interacting particle system embedded in a branching random environment (Sinclair, 2024).

5. Genealogical, pathwise, and macro-particle embeddings

Several models use embedded branching structures that are not ordinary spatial particle clouds but nevertheless preserve the essential genealogy.

In the multifractal crossing-tree construction, a continuous path $t$29 is decomposed into crossings between dyadic spatial levels. The level-$t$30 passage times are

$t$31

and the $t$32th level-$t$33 crossing is the path segment between $t$34 and $t$35. Each crossing becomes a node in a tree, its finer subcrossings become offspring, and crossing orientations $t$36 and $t$37 define a two-type Galton–Watson structure. Under Assumption 2.1,

$t$38

the canonical embedded branching process exists as a continuous process with discrete self-similarity

$t$39

Random edge weights then turn the crossing tree into a multitype branching random walk and generate a cascade time change $t$40, producing the multifractal process

$t$41

The embedded branching process is thus the path’s crossing hierarchy itself (Decrouez et al., 2012).

The paper explicitly identifies Brownian motion as a canonical embedded branching process. In that case the offspring word consists of a geometric$t$42 number of excursion pairs followed by the terminal direct crossing, so multifractal time-changed Brownian motion appears as a special case of this embedding (Decrouez et al., 2012).

A different macro-level embedding is used for branching processes in random environment with sibling dependence. At the individual level, siblings do not reproduce independently: a type-$t$43 sibling group at generation $t$44 has joint offspring vector distributed according to $t$45. Conditional independence is restored by grouping each sibling group into a “macro particle.” The associated macro process

$t$46

counts sibling groups by size, and the original population size is recovered as

$t$47

The macro process is an $t$48-type Galton–Watson process in random environment, with mean matrices related by

$t$49

The embedding therefore shifts the state space from individuals to family clusters (Vatutin et al., 2018).

A similar state-space enlargement appears in mutually catalytic branching particle systems on $t$50. The particle model itself is not reduced to an exact embedded Galton–Watson chain, but the event-time jump chain

$t$51

and especially the total particle count at jump times form an embedded discrete object used in the construction and non-explosion analysis. In finite tori, the renormalized total masses converge to the mutually catalytic diffusion

$t$52

which is a coarse-grained effective branching mechanism rather than an exact genealogical embedding (Fugenfirov et al., 2023).

These examples show that “embedded” may refer to several distinct operations: extracting a closed coordinate process, grouping dependent particles into macro-particles, reading genealogy from a path decomposition, or lifting a constrained system to a branching system in an enlarged state space.

6. Asymptotic regimes, applications, and recurrent themes

Across the models surveyed here, several recurrent asymptotic regimes define the subject.

One is the crossover from free branching to density-limited branching. In Euclidean embedded trees, this appears as the transition from

$t$53

to

$t$54

or, under all-history exclusion,

$t$55

The same model also produces extinction in bounded domains and a time-driven transition from small-world to large-world tree architecture (Forgerini et al., 2011).

A second regime is selection-induced slowdown governed by an embedded frontier population. In the censored Galton–Watson reduction of branching-selection dynamics, the slowdown satisfies

$t$56

and the extinction time of the embedded censored process is exponentially large in $t$57 after scaling by $t$58 (Couronné et al., 2011).

A third regime is hydrodynamic replacement of adaptive selection by deterministic killing. In the $t$59-branching Markov process, the large-$t$60 empirical measure converges to a survival-conditioned law above a moving boundary, and the leftmost particle converges to that boundary. This suggests that the embedded killed branching system captures the macroscopic effect of instantaneous pruning (Bérard et al., 2023).

A fourth regime is intermittency and localization in random media. In BRWRE and parabolic Anderson models, the main part of the expected mass concentrates on favorable islands, and higher moments are dominated by small embedded branching skeletons. The location and size of these islands depend strongly on the upper tail of the potential; for double-exponential tails, the survey reports bounded-size intermittent islands (König, 2020).

A fifth regime is heavy-tailed breakout-driven growth. In semipushed fronts, rare barrier-hitting particles create descendant clusters with $t$61-stable tail, and the rescaled population converges to an $t$62-stable CSBP on time scale $t$63 for $t$64. The embedded branching object there is neither the full branching Brownian motion nor a simple Galton–Watson process, but a jump skeleton of breakout events and descendant masses (Tourniaire, 2021).

Applications in the cited literature are correspondingly diverse. The Euclidean exclusion model is presented as a prototype for “tree of life” evolution and overcrowding-limited diversification (Forgerini et al., 2011). BRWRE on the hypercube is interpreted as a mutation–selection system on genotype space (König, 2020). The crossing-tree model produces multifractal time changes and includes Brownian motion as a special case (Decrouez et al., 2012). The two-type lattice BRW includes an epidemic model with infected and immunity-generated particles, for which moment asymptotics and local non-intermittency are studied (Makarova et al., 2022).

A common misconception is that an embedded branching particle process must itself be a standard branching random walk or branching Brownian motion. The surveyed literature shows a broader landscape. In some papers the branching object is exact and autonomous; in others it is auxiliary, censored, killed, weighted, or genealogically coded. In still others, such as electrostatic gases on random tree boundaries, the branching structure belongs to the environment rather than to the particle dynamics themselves (Sinclair, 2024).

The unifying principle is structural rather than taxonomic. An embedded branching particle process is a representation in which branching genealogy remains the central organizing mechanism, but the physically or analytically relevant state may be a Euclidean frontier, a censored top class, a killed auxiliary colony, a measure-valued population, a Lévy-excursion level set, a family-cluster macro-particle chain, or a multiscale crossing tree. This suggests that the concept is most useful as a bridge between branching-process theory and more complex stochastic systems, especially when exact independence is partially hidden by geometry, selection, interaction, or pathwise encoding.