GASTON-GP: Source catalogue and millimetre variability of massive protostellar objects

This study analyzes four years of millimetre variability data from the GASTON-GP survey of massive protostellar objects and reports no robust detections of variable protostars, a null result attributed to observational limitations and the absence of extreme luminosity bursts during the monitoring period.

The Gist

Imagine the universe as a giant, cosmic construction site. For decades, astronomers have been trying to figure out how the "buildings" of this site—stars—are constructed. Specifically, they want to know: Do stars grow up slowly and steadily, like a child eating a steady diet, or do they have sudden, massive growth spurts, like a teenager hitting puberty all at once?

This paper, titled "GASTON-GP: Source catalogue and millimetre variability of massive protostellar objects," is a report from a team of astronomers who went looking for those growth spurts in the most massive baby stars (protostars) in our galaxy.

Here is the story of their search, explained simply:

1. The Big Question: The "Luminosity Problem"

Astronomers have long noticed a puzzle. If you look at baby stars, they seem much dimmer than they "should" be if they were growing at a constant, steady rate. It's like seeing a teenager who eats a whole pizza every day but somehow stays the same size.

The leading theory to solve this is Episodic Accretion. The idea is that these stars spend most of their time sleeping (growing slowly), but occasionally, they have a massive "feast" where they gulp down huge amounts of gas and dust in a short burst. This feast makes them shine incredibly bright for a while before they go back to sleep.

We've seen these feasts in small, low-mass stars. But what about the giants? The massive stars that will become the most powerful beacons in the universe? Do they have even bigger, more dramatic feasts?

2. The Detective Work: The GASTON-GP Survey

To find out, the team used the IRAM 30-meter telescope in Spain, equipped with a super-sensitive camera called NIKA2. Think of this camera as a high-speed, multi-spectral eye that can see in "millimeter waves" (a type of light that passes through the thick dust clouds where stars are born).

They pointed this camera at a specific patch of the Milky Way (the Galactic Plane) for four years. They didn't just take one picture; they took 11 snapshots over that time, like a time-lapse video of a construction site.

• The Target: They looked at a region about 2.4 square degrees (roughly the size of 10 full moons) centered on a specific point in the sky.
• The Goal: To see if any of the thousands of baby stars in that patch suddenly got brighter or dimmer between the snapshots.

3. The Challenge: Finding a Needle in a Haystack

The team faced two main problems:

1. The Noise: The telescope data was a bit "noisy" (like a radio with static). The conditions changed slightly every time they took a picture. To fix this, they had to invent a special mathematical "tuning" method to make sure that if a star looked brighter, it was actually brighter, and not just because the telescope was having a good day.
2. The Crowd: The region they looked at was incredibly crowded. Imagine trying to watch a single person in a packed stadium from a mile away. If that person stands up (has a burst), can you see them? Or are they hidden by the thousands of people around them?

4. The Results: The Great Silence

After analyzing 2,925 compact sources (potential baby stars) at one wavelength and 1,713 at another, the team looked for the "growth spurts."

• The Finding: They found zero massive baby stars that showed a clear, dramatic burst of light during those four years.
• The One Oddball: They did find one object that changed brightness significantly. However, after investigating, they realized it wasn't a baby star at all. It was likely a weird, isolated object (maybe a dying star or a planetary nebula) that had nothing to do with star formation.

In short: The massive baby stars in this region were remarkably calm. They didn't have the massive feasts the astronomers were hoping to catch.

5. Why Didn't They See Anything?

The authors explain this silence with two main reasons:

• The "Crowded Room" Effect: The telescope's view is a bit blurry (like looking through a foggy window). When they look at a baby star, they are actually looking at a whole cluster of them mixed together. If one star has a burst, its light gets diluted by the light of its neighbors, making the burst too faint to notice. It's like trying to hear a single person whisper in a roaring crowd.
• The Rarity of the Feasts: Computer simulations suggest that for massive stars, these huge bursts (where the star gets 100 times brighter) are actually very rare. They might only happen once every few thousand years. Since the team only watched for four years, they might have just missed the party.

6. The Conclusion: Keep Watching

The paper concludes that while they didn't find the "smoking gun" of a massive star burst, they didn't disprove the theory either. They just didn't have the right tools to see it clearly yet.

The Takeaway:
To catch these massive stars in the act of feasting, we need:

1. Sharper eyes: Telescopes with higher resolution (to separate the stars in the crowd).
2. Faster cameras: Taking pictures more often to catch the short bursts.

For now, the massive stars in this part of the galaxy seem to be growing quietly, keeping their secrets hidden in the dust. But the hunt continues!

Technical Summary

Here is a detailed technical summary of the paper "GASTON-GP: Source catalogue and millimetre variability of massive protostellar objects" by Zhou et al. (2023).

1. Problem Statement

The mechanisms governing protostellar mass growth, particularly the role of episodic accretion, remain a subject of intense debate. While "luminosity bursts" (accretion events increasing luminosity by factors of 10–100) have been observed in low-mass protostars and a handful of massive protostars (e.g., S255IR NIR3, NGC6334-MM1), the overall statistics for massive objects are limited.

• The Challenge: Massive protostars are rare, embedded in crowded environments, and their variability is difficult to characterize at (sub-)millimetre wavelengths where the flux response to luminosity changes is diluted compared to infrared wavelengths.
• The Gap: Previous surveys (e.g., SCUBA-2 Transient Survey) focused largely on low-mass regions. There was a lack of large-scale, multi-epoch millimetre data for massive protostellar sources in the Galactic plane to constrain the frequency and amplitude of accretion bursts.

2. Methodology

The authors utilized data from the GASTON-GP (Galactic Star Formation with NIKA2 - Galactic Plane) programme, a large programme on the IRAM 30m telescope.

• Observations:

  • Instrument: NIKA2 (New IRAM KID Array 2), a dual-band camera observing at 1.15 mm (260 GHz) and 2.00 mm (150 GHz) simultaneously.
  • Field: A $\sim$2.4 deg$^2$ region of the Galactic plane centered at $l=24^\circ$.
  • Cadence: 11 observing runs (epochs) conducted between October 2017 and November 2021 (4-year baseline).
  • Resolution: 11.6$^{\prime\prime}$ at 1.15 mm and 18$^{\prime\prime}$ at 2.00 mm.

• Data Processing & Source Extraction:

  • Reduction: Used the piic software to convert time-series data into calibrated continuum images, removing sky noise and instrumental instabilities.
  • Extended Emission Removal: Applied the Nebuliser tool with a 60-arcsec spatial filter to isolate compact sources (clumps) from large-scale diffuse emission.
  • Extraction: Used the Astrodendro algorithm on Signal-to-Noise Ratio (SNR) maps to identify hierarchical structures. They defined "leaves" (smallest structures) as compact sources.
  • Catalogues: Generated catalogues of 2,925 sources at 1.15 mm and 1,713 sources at 2.00 mm.

• Physical Property Derivation:

  • Velocities: Derived using NH$_3$ (RAMPS survey) and CO isotopologues (FUGIN survey). A decision tree assigned flags (0–7) based on tracer reliability (NH$_3$ being most reliable).
  • Distances: Calculated using the Bayesian Distance Calculator (BDC) with kinematic distances, resolved using maser parallax associations where available.
  • Mass & Size: Derived from integrated flux, dust temperature (from Herschel colour maps), and distance, assuming optically thin dust emission.

• Variability Analysis:

  • Calibration: Developed a dedicated relative calibration scheme to correct for run-to-run flux variations. This involved selecting 36–37 bright, stable sources as calibrators and calculating a "flattening parameter" for each run to normalize the data.
  • Target Sample: Focused on $\sim$200 high-SNR sources (SNR > 5 in at least 8 runs) to construct light curves.
  • Criteria for Variability: A source was considered a candidate if:
    1. Variability significance $s_{var} \ge 3$ (deviation from mean relative to uncertainty).
    2. Correlated variability at both 1.15 mm and 2.00 mm (Pearson correlation $p < 0.05$).

3. Key Contributions

1. First Large-Scale Millimetre Catalogue: Provided the first deep, multi-epoch catalogues of massive protostellar clumps in the Galactic plane, containing nearly 3,000 sources with derived physical properties (mass, size, distance, temperature).
2. Advanced Calibration Technique: Demonstrated a robust relative calibration method that reduced flux uncertainty from $\sim$14% (absolute) to 3.4% at 1.15 mm and 2.6% at 2.00 mm, enabling the detection of subtle variability.
3. Variability Constraints: Conducted the first systematic search for millimetre variability in a large sample of massive protostars, providing strong constraints on the frequency of accretion bursts.

4. Results

• Source Properties:

  • Median distance: 5.2 kpc.
  • Median mass: 20 M$_\odot$ (ranging from 0.014 to $8.7 \times 10^3$M$_\odot$).
  • Median radius: 0.31 pc.
  • The sample is dominated by low-mass clumps, reflecting the high sensitivity of the survey.

• Variability Findings:

  • No Robust Protostellar Variables: Despite analyzing $\sim$200 high-SNR sources, no robust variable protostellar sources were detected.
  • False Positives: Three sources initially appeared variable ($s_{var} \ge 3$), but their light curves correlated with known poor-quality observing runs, indicating they were calibration residuals rather than intrinsic variability.
  • Non-Protostellar Variable: One source (G024.485+0.614) showed a strong variability ($s_{var} = 6.38$) at 2.00 mm. However, spectral energy distribution (SED) analysis (negative spectral index, lack of surrounding dust) confirmed it is not a protostar (likely a planetary nebula or radio source).

• Comparison with Known Variables:

  • The authors cross-matched known variable sources (e.g., G24.32+0.14, which showed a mid-IR burst) with GASTON-GP data. No corresponding millimetre burst was detected, highlighting the difficulty of detecting variability at these wavelengths.

5. Significance and Discussion

• Absence of Bursts: The lack of detected variability is attributed to two main factors:

  1. Observational Limitations: The single-dish beam size ($\sim$11.6$^{\prime\prime}$) dilutes the flux variation. If a massive clump fragments into many cores (as suggested by ALMA), a burst in one core is diluted by the emission of non-variable neighbors.
  2. Burst Amplitude: The data suggests that bursts with luminosity increases $< 100$-fold are undetectable with current sensitivity and resolution. Simulations suggest that while frequent, smaller bursts are common, the rare, massive bursts ($>100\times$) required for detection in this dataset are statistically unlikely to occur within the 4-year window for the specific sample size.

• Implications for Theory:

  • The results are consistent with isolated core collapse simulations (e.g., Meyer et al. 2021) which predict that only a small fraction of the stellar mass is accreted during short, high-amplitude bursts.
  • The study highlights that millimetre variability is significantly harder to detect than infrared variability due to the Rayleigh-Jeans tail of the SED, where flux changes are smaller for a given luminosity change.

• Future Outlook:

  • Detecting accretion bursts in massive protostars requires high-resolution (to resolve individual cores within clumps) and high-cadence surveys.
  • Future facilities like AtLAST (single-dish) and interferometers like NOEMA and ALMA are necessary to overcome beam dilution and sensitivity limits to fully constrain the accretion histories of massive stars.

In conclusion, the GASTON-GP survey provides a comprehensive static view of massive star-forming clumps but demonstrates that current single-dish millimetre surveys are insufficient to capture the episodic accretion events predicted by theory, necessitating a shift toward higher-resolution interferometric monitoring.