Effective Satellite Counts: Galaxies & Constellations

• Effective Number of Satellites is a measure that corrects raw satellite counts by accounting for biases, incompleteness, contamination, and kinematic criteria.
• Empirical scaling relations and methods like background correction, radial profile integration, and luminosity function integration are used to quantify satellite populations.
• Applications span galaxy formation studies in ΛCDM cosmology to optimizing satellite constellation designs for continuous coverage and reliable connectivity.

The effective number of satellites quantifies the statistically corrected count of physically associated (e.g., luminous or gaseous) satellites for a given host galaxy or system, taking into account selection effects, incompleteness, contamination, and, when relevant, kinematic association. This concept underpins quantitative studies of galaxy formation, low-mass halo occupation, and the hierarchical assembly of structure in both observational and simulation-based contexts. In satellite constellation engineering, the effective number prescribes the absolute minimum number of satellites needed to meet predefined coverage and connectivity requirements subject to the explicit constraints of geometry and hardware.

1. Empirical Scaling Relations for Luminous Satellite Counts

State-of-the-art ΛCDM cosmological simulations and deep observations reveal a power-law relationship between the effective number of luminous satellites and the mass of the host galaxy. Specifically, the median cumulative number of satellites above a fixed stellar mass threshold ($M_* > 10^5\,M_\odot$) scales as

$$N_{\rm sat}(M_{\rm host}, M_* > 10^5\,M_\odot) = N_0 \left( \frac{M_{\rm host}}{M_{\rm ref}} \right)^\alpha$$

with $M_{\rm ref} = 10^{11}\,M_\odot$, $N_0 \simeq 3.0$, and $\alpha \simeq 1.0 \pm 0.2$ as inferred from the Auriga simulations (Pardy et al., 2019). Equivalently, in terms of host stellar mass ($M_*^{\rm pri}$):

$$N_{\rm sat}(M_*^{\rm pri}, M_* > 10^5\,M_\odot) \simeq 3.0 \bigg( \frac{M_*^{\rm pri}}{2 \times 10^9\,M_\odot} \bigg)^{1.0}$$

Scatter about the median is approximately a factor of two. The observed satellite counts in the Magellanic system (SMC, Fornax, Carina; $M_*>10^5\,M_\odot$) match these predictions, validating the statistical procedure for effective satellite enumeration under $\Lambda$CDM (Pardy et al., 2019).

2. Methodologies for Deriving the Effective Number

Galaxy satellite studies employ a suite of methodologies to obtain the effective number, depending on sample selection (photometric, spectroscopic, or multi-wavelength):

• Background-corrected direct counts: For a set of isolated hosts, satellites are counted within specified radial and luminosity/mass windows, subtracting projected background counts using control annuli (Lares et al., 2010).
• Radial projection profile integration: Measured or fitted surface-density profiles $\Sigma(r_p) = A r_p^{-\alpha}$ are integrated over the defined annulus to yield

$$N_{\rm eff} = 2\pi A \frac{R_{\rm max}^{2-\alpha} - R_{\rm min}^{2-\alpha}}{2-\alpha}$$

(Lares et al., 2010).

• Luminosity function integration: The differential luminosity function (Schechter or power law) is integrated down to the detection threshold:

$$N_{\rm eff} = \int_{L_{\rm min}}^{\infty}\Phi(L)dL = \Phi^*\Gamma(1+\alpha_L,\, L_{\rm min}/L^*)$$

where $\alpha_L$ is the faint-end slope (typically $-1.3\pm0.2$ for satellites), and $\Gamma$ is the incomplete gamma function (Lares et al., 2010, Nierenberg et al., 2012).

• Membership and completeness correction: For HI or multi-phase satellites, the effective number is given by

$$N_{\rm eff} = \sum_{i\in \text{candidates}} \frac{p_i}{f_{\rm det}(M_i)}$$

where $p_i$ is the membership probability and $f_{\rm det}(M_i)$ is detection completeness as a function of mass and linewidth (Zhu et al., 2023). For example, with raw observed counts $N_{\rm obs}\simeq2$ and completeness $f_{\rm det}\simeq 0.5$, $N_{\rm eff}\simeq 4$.

In all cases, uncertainties are propagated from Poisson errors and sample variance.

3. Kinematic and Phase-Space Selection Criteria

For robust assignment of satellites to a host, especially in the context of the Milky Way–LMC system, effective number calculations incorporate dynamical constraints:

1. Angular (sky) proximity: Satellites are required to be within a certain angular separation from the host, e.g., $\leq 30^\circ$ (Pardy et al., 2019).
2. Phase-space consistency: The sign and magnitude of line-of-sight velocity relative to the host must be compatible with the expected motion (e.g., for satellites on first approach).
3. Orbital angular momentum alignment: The angle $\alpha$ between satellite and host orbital angular momentum vectors is constrained, typically $\alpha \leq 30^\circ$:

$$\cos\alpha = \frac{\mathbf{L}_{\rm host}\cdot\mathbf{L}_{\rm sat}}{|\mathbf{L}_{\rm host}||\mathbf{L}_{\rm sat}|}$$

Satellites passing all three criteria constitute the kinematically defined effective sample (Pardy et al., 2019).

This step is essential in disentangling physically associated satellites from foreground/background interlopers, particularly in cosmological surveys where projection effects are significant.

4. Application to LEO Satellite Constellations

In satellite constellation engineering, the "effective number" denotes the minimum number of satellites, $T_{\rm min}$, required to guarantee continuous global or regional coverage and maintain uninterrupted inter-satellite connectivity:

• Coverage constraint: For a given arc coverage $\lambda^*$ (angular radius, set by payload boresight or target elevation), every point in the area of interest must be within $\lambda^*$ of some satellite at all times.
  • Bounded Voronoi tessellation and APC (Access-Profile-Constellation) decomposition are analytic/global geometric guarantees for coverage (Jeon et al., 2024).
• Connectivity constraint: Ensures angular separations $[\theta_{\rm min}, \theta_{\rm max}]$ between adjacent-plane satellites fall within ISL hardware constraints, governed by closed-form analytic expressions for the Walker–Δ and related constellations (Jeon et al., 2024).
• Optimization: Binary Integer Linear Programming (BILP) is used to minimize $T$ subject to coverage (via circulant seed-access matrices) and ISL constraints.

For example, for single-fold coverage and ISL continuity, BILP-optimized common ground-track constellations can achieve $T_{\rm min}=31$ for the same coverage where Walker–Δ requires $40$, but with greater ISL variation (Jeon et al., 2024).

5. Observational Trends and Physical Interpretation

The effective number varies systematically with host characteristics and redshift:

• Mass dependence: Early-type galaxies with $\log_{10}M_*>11$ exhibit higher values of $\Sigma_0$ and $N_{\rm eff}$ compared to late-types of the same mass; for example, $\Sigma_0\simeq1.7$ (early) versus $\simeq0.8$ (late) for $0.1Nierenberg et al., 2012).
• Luminosity dependence: For primaries with $M_r<-21.5$, the mean number of satellites down to $M_r=-14.5$ is $\sim6$, while for $-21.5<M_r<-20.5$, $\langle N_{\rm sat}\rangle<1$ (Lares et al., 2010).
• Morphology and alignment: Satellites of early-type hosts display flattened, aligned angular distributions; late-type hosts show isotropic angular patterns (Nierenberg et al., 2012).
• Redshift evolution: No statistically significant redshift dependence detected over $0.1Nierenberg et al., 2012).
• Membership probability and detection completeness: HI-selected satellites require explicit correction for detection completeness, often a factor $\sim2$ at $M_{\rm HI}\sim7\times10^6\,M_\odot$ at 10 Mpc (Zhu et al., 2023).

6. Practical Calculation Workflow

A typical end-to-end workflow for obtaining the effective number of satellites (luminous or gaseous):

1. Estimate host mass ($M_{200}$ or $M_*^{\rm pri}$) via abundance matching or kinematic indicators.
2. Select candidate satellites meeting radial, luminosity/mass, and phase-space/kinematic criteria.
3. Apply corrections for background/foreground contamination, detection completeness, and (when appropriate) membership probability.
4. Integrate spatial and luminosity functions over targeted domains, normalizing as required.
5. Compare effective counts against theoretical, simulation, or alternative survey predictions to evaluate hierarchical galaxy formation or satellite quenching scenarios.

7. Significance in Broader Contexts

The effective number of satellites is central to testing $\Lambda$CDM predictions at dwarf-galaxy scales, constraining baryonic processes (e.g., quenching, feedback), and guiding the design of space-based distributed systems. Discrepancies between predicted and observed effective numbers—after all selection effects are accounted for—indicate possible gaps in the theoretical understanding of halo occupation, environmental processing, or observational methodology. Quantitative alignment between simulation prediction and observed effective numbers, as documented for the LMC and its satellites, provides compelling evidence for the hierarchical assembly paradigm (Pardy et al., 2019). In satellite constellation engineering, optimization of the effective number directly informs mission cost, latency, and operability (Jeon et al., 2024).