RKSTAR Survey: Mapping Stellar Radio Emission

• RKSTAR Survey is a technical framework that exploits the SKA's capabilities to systematically catalog faint stellar radio sources such as OB stars, Be stars, M-dwarf flares, and UCHII regions.
• The survey employs multi-tiered sensitivity strategies with SKA1-Mid and SKA2 configurations to achieve continuum sensitivities down to 10 nJy and high-resolution imaging for detailed population analysis.
• It integrates empirical radiative scaling relations and the Besançon Milky Way model to predict source yields and guide survey strategies for both wide-area fields and targeted deep observations.

The RKSTAR Survey is a technical framework for exploiting the sensitivity and survey capabilities of the Square Kilometer Array (SKA) to probe the radio emission from galactic stellar populations, including OB stars, Be stars, M-dwarf flares, and Ultra Compact HII (UCHII) regions. This approach utilizes high-resolution, deep imaging of both the Galactic plane and off-plane fields to systematically catalog and characterize faint stellar radio sources, using the Besançon Milky Way model for spatial distribution and tailored radiative scaling relations for flux prediction. The design emphasizes multi-tiered sensitivity strategies, spectral coverage, and temporal monitoring to enable population-scale statistical inference of radio luminosity functions, ionized wind diagnostics, and flare burst rates from diverse stellar classes (Yu et al., 2021).

1. SKA Array Configurations and Observing Setups

The RKSTAR Survey is predicated on leveraging both SKA1-Mid and SKA2 configurations for optimal coverage and sensitivity. SKA1-Mid operates between 350 MHz and 14 GHz, with survey observations concentrated at 5 GHz using a bandwidth (Δν) of 1.5 GHz (and up to ≲3.9 GHz for the full band). Synthesized beam sizes vary from 0.17″ to 23″, contingent on uv-weighting schemes, yielding continuum sensitivities of 1.4 µJy (5σ) in 1 hour (σ≈280 nJy) and σ≈60 nJy in 8 hours (5σ≈300 nJy). Natural weighting in fields with sparse sources can provide a sensitivity gain of ≲2×.

SKA2 extends these capabilities, featuring ∼2000 array elements over baselines of several thousand kilometers, and bandwidths Δν ≳ 4 GHz. The expected sensitivity for 8-hour deep integrations is ≲10 nJy (5σ, σ≈2 nJy), with synthesized beams of ≲0.1″ for high-resolution imaging. The confusion limit at 5 GHz is negligible for θ ≲ 1″, facilitating ultra-faint source detection.

Configuration	Δν (GHz)	Beam θ	σ (1 h, 5σ)	σ (8 h, 5σ)
SKA1-Mid	1.5 (3.9)	0.17″–23″	1.4 µJy	300 nJy
SKA2	≳4	≲0.1″	—	10 nJy

2. Quantitative Relations for Stellar Radio Emission

The translation of stellar parameters to expected radio flux relies on established luminosity scaling relations and radiative transfer models:

1. General flux–luminosity conversion:

$S_\nu = \dfrac{L_\nu}{4\pi d^2}$ where $S_\nu$ is the observed flux density, $L_\nu$ is the spectral luminosity, and $d$ is the source distance.

2. OB Stars (thermal free–free wind; Wright & Barlow 1975):

$$\dot{M} = \dfrac{3.01 \times 10^{-6}\,\mu}{Z\sqrt{\gamma\,g_{\rm ff}(\nu)}\,V_{\infty}} S_\nu^{3/4} d^{3/2}$$ $$g_{\rm ff}(\nu) = -1.66 + 1.27 \log \left( \dfrac{T_{\rm wind}^{3/2}}{Z\nu} \right)$$ Theoretical wind prescription (Vink et al. 2000; Vink 2011) uses a piecewise log \dot{M} formula as a function of stellar effective temperature.

3. Be Stars (Taylor et al. 1990 empirical scaling):

For SpT < B5: $$L_{14\,\text{GHz}} = 2.5\times 10^{16}\,\text{erg\,s}^{-1}\,\text{Hz}^{-1}$$; SpT ≥ B5: $$L_{14\,\text{GHz}} = 0.2\times 10^{16}\,\text{erg\,s}^{-1}\,\text{Hz}^{-1}$$ (±50% scatter);

$$S_{5\,\text{GHz}} = \left( \dfrac{5}{14} \right)^{1.3} \dfrac{L_{14}}{4\pi d^2}$$

4. Ultra-compact H II regions (Zijlstra & Pottasch 1989):

$$S_\nu / F(\mathrm{H}\beta) = 2.51\times 10^7\, T_e^{0.53}\, \nu^{-0.1}\, Y$$, where $T_e = 10^4\,\text{K}$, $Y \approx 1.1$ for He/H.

5. M-dwarf Flares (Villadsen & Hallinan 2019): The 10-min flux density PDF at 1 pc (2.8–4 GHz) is empirically sampled, flares are assigned a duration τ_flare ≈ 20 min yielding an integrated flux ≈2× the 10-min value, and fluxes are scaled to 5 GHz using the observed spectral trend.

3. Galactic Population Modeling and Simulation

Stellar population synthesis is conducted using the Besançon model, facilitating predictions for detection yields and spatial distributions:

• OB Stars: Modeled over 0–30 kpc, spanning l = 0–360°, b = ±1°, in Δl = 1°, Δb = 0.2° bins. No extinction is applied. Output yields N_OB ≃ 6.33×10^6 within the full region.
• Be Stars: Defined as a random 20% subset of B-type stars in the OB sample.
• UCHII Regions: Represented as a fraction f_UCHII ≈ 3% of OB stars, reflecting phase lifetime.
• M-dwarf Flares: Surveyed in a patch l = 0–10°, b = ±20° (~400 deg²), including SpT = M0–M9, V < 23, yielding N_M ≈ 6.82×10^7. The fraction actively flaring (f_act) is assigned per SpT from West et al. (2008): M0–M9 → [6%, 6.8%, 10%, 16.7%, 31.8%, 66.1%, 68.9%, 83.8%, 98.6%, 90.2%]. Flares are distributed per the 10-min empirical PDF and scaled for distance and duration.

4. Predicted Detection Yields and Source Statistics

Sensitivity limits dictate survey yield across the various stellar radio populations. At 5 GHz:

• F_lim = 10 nJy (5σ, 8 h, SKA2):
  • OB stars: ≃50 deg⁻²
  • Be stars: ≃1500 deg⁻²
  • M-dwarf flares: ≃4500 deg⁻²
  • UCHII regions: all detectable (≳1 µJy), ≃1000 deg⁻² toward Galactic Center (GC)
• F_lim = 100 nJy (5σ, 8 h, SKA1):
  • Results in ≃8× fewer sources per class (OB: ≃6 deg⁻², Be: ≃190 deg⁻², M flares: ≃560 deg⁻²).
• Galactic longitude dependence (|b| ≤ 1°, l = 0–10° per deg²):

Class	10 nJy	100 nJy
OB	70.6	8.8
Be	1466	185
UCHII	1031	1031
M flares	1091	136

Cumulative N(>F) curves for all classes scale as F^{-0.8}, indicating a steep faint-end slope that underscores the need for ultra-sensitive imaging (Yu et al., 2021).

5. Survey Strategy and Design Recommendations

Survey strategy is tailored for both high completeness and depth, balancing wide-area coverage against ultra-sensitive, targeted deep fields:

A. Frequency & Bandwidth:

Primary survey at 5 GHz, minimizing synchrotron confusion and maximizing sensitivity to thermal emission. Simultaneous L-band (1–2 GHz) sub-arrays are recommended for spectral characterization of M-dwarf flares.

B. Sensitivity Tiers:

• Wide survey (SKA1): 300 nJy (5σ, 8 h) across the Galactic plane (|b| ≤ 1°) to catalog UCHII regions and the brightest stars.
• Deep fields (SKA2): 10 nJy (5σ, 8 h) in selected regions (l = 0–60°) to capture faint OB/Be and flare populations.

C. Sky Coverage:

• Plane: |b| ≤ 1° for bulk UCHII and hot stars; extended to |b| ≤ 3° in l = 90–270° for warped disk structure.
• Off-plane pilot: |b| = 1°–20° (l = 0–10°) for M-dwarf flare background assessment.

D. Cadence & Imaging Depth:

Deep fields employ 8-h continuous blocks with 20-min sub-integrations to resolve transient activity. Shallow surveys can use a single 8-h exposure or split into ≤4×2 h repeats to probe variability. Imaging dynamic range must exceed 10^6 to accommodate both bright UCHII and faint stellar radio events.

E. Data Products:

• Continuum maps at σ ≲ 10 nJy (deep) / 300 nJy (wide)
• Time-resolved light curves at 20-min cadence for flare statistics.
• Source catalogs with classification tags (OB, Be, UCHII, flare).

6. Technical Justification: Key Equations and Source Yield Tables

The RKSTAR Survey is underpinned by explicit equations and tabulated yields:

Key Equations:

• (1) $S_\nu = \dfrac{L_\nu}{4π\,d^2}$
• (2) $$\dot{M} = \dfrac{3.01\times10^{-6}\,\mu}{Z\sqrt{\gamma\,g_{\rm ff}(ν)}\,V_\infty}\,S_\nu^{3/4}\,d^{3/2}$$
• (3) $$g_{\rm ff}(ν) = -1.66 + 1.27\,\log\left(\dfrac{T_{\rm wind}^{3/2}}{Z\,ν}\right)$$
• (4) $$S_{5\,\text{GHz}} = \left(\dfrac{5}{14}\right)^{1.3}\dfrac{L_{14\,\text{GHz}}}{4π\,d^2}$$
• (5) $$L_{Hβ} = L_\star\,\dfrac{h\,ν_{Hβ}}{π^4\,k\,T_\star} \dfrac{α_{Hβ}}{α_B}\,15\,G_1(T_\star)$$
• (6) $$S_\nu/F(\mathrm{H}\beta) = 2.51\times10^7\,T_e^{0.53}\,ν^{-0.1}\,Y$$

Source Yield Table:

Class	F_lim = 10 nJy	F_lim = 100 nJy
OB stars	~50	~6
Be stars	~1500	~190
M-dwarf flares	~4500	~560
UCHII regions	≳1000	≳1000

Figures (as described) include cumulative N(>F) curves (F^{-0.8} scaling) and Galactic longitude profiles (|b| ≤ 1°), delineating the spatial and flux-dependent survey returns.

A plausible implication is that SKA-enabled imaging at the specified thresholds will profoundly deepen the sample size and diversity of stellar radio sources, revealing large reservoirs of faint transients and low-luminosity wind-driven emission previously below instrumental detection limits. These capabilities are critical for population studies of stellar magnetism, mass-loss, and Galactic stellar structure (Yu et al., 2021).