1 of 19

A Scientific exploration only possible with Rubin-LSST: Search for hidden baryons through interstellar scintillation

With a 0.06 image/s movie of LMC/SMC at Rubin observatory

SCOC Workshop

Marc Moniez, 7 may 2025

2 of 19

Search for missing H2 turbulent galactic gas

through scintillation detection (the OSER project)

Light received by telescope varies with

- timescale ~10 min (due to the relative velocity of the gas)

- modulation of a few % (depending on distances / turbulence parameters / source extension)

3 of 19

Distance scales

4 distance scales characterize the speckle pattern

  • Diffusion radius Rdiff
    • separation (on screen) such that: σ[φ(r+Rdiff)-φ(r)] = 1 radian
    • Characterizes the strength of turbulence: strong turbulence -> small Rdiff)

  • Refraction radius Rref
  • size of the region from which most of the scattered signal, seen by an observer, originates ~ z0λ/Rdiff
  • Larger scale structures of the diffusive gas can play a role if focusing/defocusing configurations happen
  • Projected source size RS�speckle from a pointlike source�is convoluted by the source�projected profile. -> Critically impacts�the contrast of the illumination pattern

Rref

4 of 19

Illumination in Ks�by a K0V star@8kpc (mV=20.4) through a cloud@160pc (B68) with turbulence parameter Rdiff =150km

4

m = σI/I = modulation index

Simulation of a realistic speckle from a stellar source crossing a (very) turbulent gaseous structure

Simulated light curve (VT=100km/s)

tref

m

monochromatic point source

stellar source in Ks

m=0.23

m=1.08

Projection of star

Convolution within projected star disk

5 of 19

Observable 1: tref (characteristic time)

tref

m

VT = relative velocity of the cloud w/r to the line of sight

-> VT is also the apparent velocity of the illumination pattern entering the telescope

VT

telescope

is here

Rref

Distance to the cloud

Observable 2: m (modulation index)

Scintillation@λ = 1 μm

of Sun@10kpc (V~20) through a cloud@160pc

with Rdiff=1000km

Essentially depends on RS and Rref -> not on the details of the power spectrum of the fluctuations

One realization

average

6 of 19

Signature of scintillation

  • Stochastic light-curve (not random)
    • Autocorrelation (power spectrum)
    • Characteristic time tref (few minutes)
    • Modulation index m can be > a few %
      • decreases with increasing star radius
      • depends on cloud structure
  • Signatures of a propagation effect
    • Chromaticity (optical wavelengths)
      • Long time-scale variations (few min.) ~ achromatic λ−1/5
      • Short time-scale variations (sub-min.) varies with λ6/5
    • Correlation between light-curves measured by 2 telescopes decreases with their mutual distance

7 of 19

Results from IR-NTT observations toward B68

-> A star scintillating through dusty gas?

1000 x 10s exposures in Ks

Night 1

Night 2

8 of 19

Scintillation with Vera Rubin telescope�(15s exposures)

For a given value of the turbulence parameter Rdiff ->

modulation m = σI /<I> decreases with the sources’ magnitude since the apparent stellar radius increases (here MS star)

Turbulent gas structures in Galactic halo

Texp = 15s

  • Source@55Kpc (LMC) in V
  • Screen@1Kpc : halo- hidden gas

LSST precision

for Texp = 15s

Turbulence strength

Scintillation

detectable

with LSST

Vsource

Typical light-curve with m ~ 1% (simulation)

9 of 19

Search for hidden baryons�with a 0.06 image/s movie of LMC at Rubin-LSST

  • To discover scintillation effects, we need:
    • Large diameter telescope(s) (>4m) because faint stars+short exposures
    • Wide field camera (visible) because small optical depth
    • fast readout -> to get highly sampled light-curves

---> Vera Rubin telescope is the unique option

  • Extremely sensitive technique to gas clumpuscules with column density fluctuations down to 10-7 --- i.e. ~1 mm of H2 under normal (P, T) conditions --- per few 1000km
  • A few hours' filming over 2 nights (for redundancy) at the LMC in G passband will provide decisive sensitivity to the hidden gas of the galactic halo, the last unknown baryonic component�� -> Early science possible during commissioning or micro-survey

Biblio : A&A 412, 105-120 (2003); A&A 525, A108 (2011); A&A 552, A93 (2013)

Marc Moniez, SCOC workshop, 7 may 2025: moniez@ijclab.in2p3.fr

>> Other communities should be interested in this movie (transits, flares, very short timescale variabilities...)

10 of 19

Complements

11 of 19

Atmosphere, atmosphere?

  • Blurs PSF, but doesn’t affect the intensity collected by a large telescope
  • ~ 5cm size speckle due to turbulent layers at ~ 10km
  • Visible just before total solar eclipses: « shadow bands »

Polar star intensity variance as a function of the aperture (2 series of measurements)

12 of 19

Fresnel diffraction on stars is a classic observation

  • In radioastronomy:

classic technique for studying the interstellar medium

  • In optics:

diffraction during asteroid/lunar occultations, clearly distinct from atmospheric effects

13 of 19

Rdiff varies with 1/Rref

Rref

Rapid optical path variations -> small Rdiff large Rref

Slow optical path variations -> large Rdiff small Rref

14 of 19

Light-curve of an A5V-LMC star

(integral in the sliding disk)

Diffraction image of�a point-like source�through this cloud @1 kpc

Simulation of a turbulent cloud

15 of 19

Point source

Extended source (0.85 Rsun)

Monochromatic

λ=2.18

Polychromatic

(Ks passband)

m = σI/I modulation index

Illumination on Earth

- in Ks passband�- by a K0V star@1.2Kpc� (mV=16.3)

- through a cloud@160pc

- with Rdiff =100km

m=100%

m=55%

tref

m

zoom

Simulation of a realistic scintillation speckle

Simulated light curve detected by a telescope

t

m=3.3%

m=3.3%

Monochromatic

λ=2.18

Polychromatic

(Ks passband)

16 of 19

Scintillation with Vera Rubin telescope�(15s exposures)

  • Source@7Kpc (gal. plane) in V
  • Screen@160pc : B68- visible gas

For given Rdiff modulation m = σI /<I> decreases with the sources’ magnitude since the apparent stellar radius increases (here MS star)

Gas structures in Galactic plane

Gas structures in Galactic halo

LSST precision

for Texp = 15s

Rdiff = 1000 km

Rdiff = 4000 km

Rdiff = 16,000 km

m = σI /<I>

V

Texp = 15s

  • Source@55Kpc (LMC) in V
  • Screen@1Kpc : halo- hidden gas

LSST precision

for Texp = 15s

Turbulence strength

Texp = 15s

Turbulence strength

17 of 19

Maximum fraction of LMC/SMC scintillating stars : 1%

For a given detection threshold of modulation

τ(m > mthreshold) < 10-2 x f(mthreshold)

Modulation

index

fraction of gas turbulent enough to produce m > mthreshold

At maximum, 1% of the sky is covered by turbulent gas

18 of 19

Series of light curves seen by telescopes 2000 km far from each other

  • Decorrelation increases with the distance
  • But statistical behaviour constant (power spectrum)

-> Strong signature of propagation effect (not intrinsic variability) with only 2 telescopes far apart

Illumination from a 0.5xRsun star at 1Kpc

through a diffusor at 160pc with Rdiff = 1000km

19 of 19

star

Point source

Phase screen

Speckle on Earth

convolution

Impact of the size of the source