# OS-0001: The star that glowed before it dimmed

> In 2018 this Sun-like star started glowing in the infrared. Three years later it began to fade. Something hot came first.

- **Star:** ASASSN-21qj (2MASS J08152329-3859234, Gaia DR3 5539970601632026752, ASASSN -21qj)
- **Position (J2000, deg):** RA 123.84708180613582, Dec -38.98980679315472
- **Status on the discovery ladder:** reproduced (observation → anomaly → reproduced → interpreted → validated; exits: artifact, closed)
- **Web page:** https://oddsky.ai/case/OS-0001
- **All data (JSON, ~0.1-0.5 MB):** https://oddsky.ai/data/cases/OS-0001.json
- **Case directory in git:** https://github.com/oddsky-org/oddsky/tree/main/cases/OS-0001

## Gaia DR3

| quantity | value |
|---|---|
| Gaia G mag | 13.371 |
| BP−RP | 0.815 |
| Teff (K) | 5761.088 |
| distance (pc) | 567.172 |
| absolute G | 4.602 |
| RUWE | 1.104 |
| non_single_star | 0 |

## Where the argument stands

A Sun-like star whose infrared suddenly brightened in 2017-18 with ~1,000 K dust, followed about 3 years later by ragged optical dips of up to 2 mag. The dips are more strongly coloured than interstellar dust, so small grains are blocking the light. The infrared excess faded as the optical recovered, and by 2025-26 the star is back to its pre-event brightness and flat to 0.005 mag; the best explanation is a one-off cloud of debris, plausibly from a planetary collision.

## Light curves

Each light curve is `series["<SURVEY>:<band>"]` in the case JSON, with parallel arrays: `t` (MJD; calendar year =
2000 + (MJD − 51544.5)/365.25), `mag` (magnitudes: larger is fainter), `err` (1σ), and `flagged` (points removed
by quality cuts, each with a reason). Dense series are binned; `binned_days` says how. `context.gaia` holds
Gaia DR3 parameters, `context.neighbors` the Gaia sources nearby, `skeptic.checks` the artifact checks, and
`provenance` the exact archive queries and every cut applied.

| series | points | span | flagged |
|---|---|---|---|
| `AAVSO:B` | 2076 | 2021-08-28 to 2026-10-09 | 4 |
| `AAVSO:I` | 2108 | 2021-08-28 to 2026-10-25 | 4 |
| `AAVSO:R` | 815 | 2024-08-18 to 2026-10-25 | 4 |
| `AAVSO:V` | 2116 | 2021-08-28 to 2026-10-25 | 4 |
| `NEOWISE:W1` | 21 | 2014-05-11 to 2024-04-20 | 4 |
| `NEOWISE:W2` | 21 | 2014-05-11 to 2024-04-20 | 4 |
| `TESS:TESS` | 167 | 2019-02-03 to 2025-03-12 | 4 |

## The Skeptic (deterministic artifact checks)

0 fail, 5 warn, 9 pass, 0 skipped

- **WARN** Neighbours inside the NEOWISE beam (6.1"): 1 Gaia source(s) within 6.1"; brightest contributes 3.8% of the target's G flux at 3.8"; 1 is much redder than the target and may matter more in the infrared
- **PASS** Neighbours inside the TESS beam (42.0"): 5 Gaia source(s) within 42.0"; brightest contributes 3.8% of the target's G flux at 3.8"
- **PASS** Neighbours inside the AAVSO beam (5.0"): 1 Gaia source(s) within 5.0"; brightest contributes 3.8% of the target's G flux at 3.8"
- **PASS** WISE saturation: all visits fainter than the W1/W2 saturation limits
- **PASS** NEOWISE single-exposure quality cuts: 417 exposures; cuts removed up to 84 (20%) per band; 0 visits kept with fewer than 3 exposures
- **WARN** Independent optical confirmation: optical data from 1 independent source(s): AAVSO observer HMB (7115). One instrument or observer is a single point of failure; seek a second. TESS is not counted: it is normalised per sector, so it confirms an event only if a sector overlaps it, and never its absolute depth.
- **WARN** Does AAVSO start at the star's normal brightness?: first 60 days of V: median 14.67; Gaia DR3 (2014-17 mean) implies V = 13.52; difference +1.15 mag. The series begins mid-event: its own median is not the star's normal state.
- **WARN** Colour of the variability (all simultaneous AAVSO data): dB/dV = 1.55 (ISM dust 1.32, r = 0.98), dR/dV = 0.66 (ISM dust 0.81, r = 0.50), dI/dV = 0.46 (ISM dust 0.59, r = 0.94). At least one band does not track V: suspect calibration, a blend, or a band-specific artifact.
- **PASS** Was the start seen? (AAVSO B): the 3 deepest dip(s) ( 0.1 mag vs a 600-day baseline) all have an observed ingress
- **PASS** Was the start seen? (AAVSO I): the 3 deepest dip(s) ( 0.1 mag vs a 600-day baseline) all have an observed ingress
- **PASS** Was the start seen? (AAVSO V): the 3 deepest dip(s) ( 0.1 mag vs a 600-day baseline) all have an observed ingress
- **WARN** Periodicity vs survey rhythms (AAVSO V): strongest period 1.00 d coincides with: matches spectral-window peak at 1.00 d; sidereal day; solar day; best period above 2 d: 889.7 d
- **PASS** Periodicity vs survey rhythms (TESS TESS): strongest period 4.13 d is not a known survey/Earth rhythm (not yet evidence it is real)
- **PASS** Gaia astrometric wobble (RUWE): RUWE = 1.10 (single-star-like astrometry)

## Claims so far

- **C1** `observed` Before late 2017 the infrared was boringly photospheric: 8 NEOWISE visits from 2014-05 to 2017-11 have W1 = 11.654 (median, scatter ~0.01 mag) and W1-W2 between -0.05 and -0.09, i.e. no detectable warm dust. (receipts: K3, K4, K14)
- **C2** `computed` The infrared brightened abruptly between consecutive NEOWISE visits on 2017-11-15 (MJD 58072) and 2018-04-19 (MJD 58227): W1 by 0.34 mag and W2 by 0.78 mag, and W1-W2 jumped from -0.09 to +0.35. Over 2018-04 to 2020-04 the median excess vs. 2014-17 was -0.41 mag in W1 (20 sigma) and -0.80 mag in W2 (39 sigma). (receipts: K3, K4, K5, K6, K14, K15)
- **C3** `inferred` The excess alone (total minus the 2014-17 photosphere) is +45% of the stellar flux in W1 and +110% in W2, giving an excess W1-W2 of about 0.9 mag (Vega); a single blackbody with that 3.4/4.6 micron ratio has T ~ 1000 K (range ~850-1050 K across visits). Hot dust, far hotter than a cold debris belt. (Hand arithmetic on tool outputs; no tool computes excess temperature.) (receipts: K14, K15, K4)
- **C4** `computed` The infrared excess faded over 2020-2023: W1-W2 = +0.32 (2020-04), +0.22 (2020-11), +0.17 (2021-11), +0.16 (2022-04), -0.01 (2022-11), -0.09 (2023-11); W1 tracks the colour (r = -0.98). The last three visits (2023-04 to 2024-04) match the 2014-17 level to within 0.01 mag in W1 and W2. NEOWISE ends in 2024, so there is no infrared after 2024-04. (receipts: K4, K16, K17)
- **C5** `computed` When AAVSO coverage begins (2021-08-26) the star is already deep in an event: the first 60 days have median V = 14.67 vs. V = 13.52 implied by the Gaia DR3 2014-17 mean (+1.15 mag). Within the event, dips reach up to 1.8 mag below a 300-day running baseline on timescales of days to weeks (faintest V 15.05, faintest B 16.51), with 18 excursions, all between 2021-08 and 2022-06. (receipts: K60, K18, K8, K59)
- **C6** `inferred` TESS sectors 8 (2019-02) and 34-35 (2021-01 to 2021-03) contain no excursion deeper than 0.05 mag and only ~1-2% day-to-day wiggles, so the deep, fast dips were not happening then, even though the infrared excess had been present since 2018-04. The infrared therefore led the dipping by at least ~2.9 years. Caveat: TESS is normalised per sector, so a smooth constant dimming during a sector would be invisible. (receipts: K29, K31, K32, K3)
- **C7** `computed` The dimming is strongly chromatic and steeper than interstellar reddening: from 2000+ near-simultaneous pairs, dB/dV = 1.55 (ISM 1.32), dI/dV = 0.45 (ISM 0.59), dB/dI = 2.86 (ISM 2.24), with r = 0.90-0.98. (receipts: K36, K37, K38)
- **C8** `inferred` A colour slope steeper than ISM dust means the occulter is not grey: it is dominated by small (sub-micron) grains, not large grains or an opaque solid body. Caveat: these are a single observer's magnitudes and the crossband slopes could absorb uncorrected colour terms at the ~10% level. (receipts: K36, K37, K38)
- **C9** `computed` Recovery: the V seasonal median went from 13.587 (2022-09 to 2023-07) to 13.506 (2024-06 to 2025-07) to 13.489 (2025-08 to 2026-06), reaching ~13.47-13.48 by mid-2023. The optical recovery (2022-23) coincides with the disappearance of the infrared excess (2022-11 to 2023-11). (receipts: K7, K44, K19, K20, K4)
- **C10** `computed` Since early 2025 the star has been flat: same-season medians Feb-Jun 2025 vs Feb-Jun 2026 differ by -0.0045 mag in V (1.5 sigma), 0.000 in B and -0.0035 in I. No excursion beyond the 4 sigma (0.14 mag) threshold after 2022-06. The current season (Oct 2026, V = 13.489, n = 16) matches the first half of 2026 (13.486). V is now equal to the Gaia-implied 13.52 within the ~0.03-0.05 mag uncertainty of a G-to-V transformation. (receipts: K55, K56, K57, K58, K8)
- **C11** `computed` A step occurs at MJD ~60721.7 (2025-02-16), across a 10-day gap: I brightens by 0.120 mag (15.7 sigma), R by 0.047, B by 0.021 and V by 0.013. The brightening grows toward the red, the opposite of dust clearing. (receipts: K46, K47, K48, K49, K50)
- **C12** `inferred` The 2025-02 step is almost certainly a calibration change by the single observer (comparison stars, aperture or filter) rather than the star. It also makes the apparent 2024/25 to 2025/26 brightening in I (-0.08) and R (-0.05) unreliable; only B and V, which barely moved, should be used to judge late recovery. (receipts: K46, K47, K48, K49, K50, K26, K27)
- **C13** `observed` Apparent 2026-06 brightening is not astrophysical: on MJD 61196.98, V is 0.37 mag and I 0.45 mag bright while B and R, taken within ~5 minutes, are normal. On 61197.98, V, B and R are 0.3-0.5 mag bright and I is only 0.12 mag bright. Real brightness changes would show up in all four filters together. These are end-of-season, low-altitude frames. (receipts: K10, K11, K12, K13)
- **C14** `inferred` The 3.85 arcsec neighbour (G = 16.91, BP-RP = 1.84, 3.8% of the target G flux) cannot produce the optical dips (losing it entirely would dim the blend by 0.04 mag). It sits inside the WISE beam: assuming G-W1 ~ 2.6 for its colour, it is ~9% of the W1 flux, so to fake the IR excess it would have to brighten ~5x in W1 and ~12x in W2 for four years, then fade exactly as the target recovers. Possible in principle (e.g. an eruptive red star), but contrived. (receipts: K60, K40)
- **C15** `observed` All optical evidence for the event in this case comes from one AAVSO observer (HMB). No TESS sector overlaps the deep dips (sectors 34-35 end 2021-03, sector 61 starts 2023-01), and ASAS-SN is unavailable. So the Skeptic's 'independent confirmation: pass' refers to the dataset, not to the event itself. (receipts: K59, K60, K29)
- **C16** `computed` TESS shows a 4.13 d periodicity (LS power 0.17, side peaks at 4.03 and 4.22 d that look like sector-gap aliases) not matching any spacecraft or Earth rhythm. AAVSO V shows no periodicity in 2-400 d beyond the trend (top peak 323 d, power 0.08, close to the yearly observing window). (receipts: K35, K61)
- **C17** `inferred` Timing tension for any single-cloud picture: blackbody grains at 1000 K around a roughly solar-luminosity star (Teff 5760 K, abs G 4.6) sit at r ~ (278/1000)^2 ~ 0.08 au, where the orbital period is ~8 days. Dust that close should cross the line of sight within weeks of appearing, not ~3 years later. So either the emitters are superheated small grains much farther out, or the hot emitter and the later occulter are not the same material (e.g. a hot body or cloud on an eccentric orbit, or dust spreading along an orbit). (receipts: K14, K15, K59)
- **C18** `computed` The new irexcess tool reproduces C3, which was hand arithmetic: over the 2014-17 photosphere, the 2018-20 excess is +44% in W1 and +108% in W2, with a W1/W2 colour temperature of 998 K. (receipts: K96)

## Hypotheses

- **H1** (collision_debris, favored): A violent dust-production event (e.g. a collision between planetary bodies) around 2018 made a hot (~1000 K) dust cloud, seen first as a sudden infrared excess. Small-grain debris from it later (2021-22) crossed our line of sight, producing ragged, strongly chromatic dips. The material then dispersed or was removed, so the infrared and optical both returned to normal by 2023-24 and have stayed there.
  - *Would falsify:* (a) Resolved infrared photometry showing the 2018-23 excess belonged to the 3.85 arcsec neighbour; (b) a repeat of the same dip pattern on a fixed period (a permanent orbiting structure, not a one-off cloud); (c) a renewed optical dimming without any new infrared excess, of similar depth and colour; (d) spectra showing the star is a young accretor with gas-rich inner disk.
- **H2** (dust_occultation, live): Generic circumstellar dust (a clump, a warp or a pre-existing structure, not necessarily collision-made) passed in front of the star in 2021-22. The infrared event is related only loosely or by coincidence.
  - *Would falsify:* The tight coupling seen here, an abrupt infrared turn-on from a photospheric baseline followed by optical dips while the IR excess decays and both recover together, already argues against a pre-existing structure. A second, independent infrared excess event with no dips (or dips with no IR) would separate them; recurrence on a period would turn this into a companion/disk model.
- **H3** (eclipse_companion, disfavored): A long-period companion carrying a dusty disk or ring system eclipsed the star (epsilon Aur / J1407-like), producing the 2021-22 dimming.
  - *Would falsify:* Already strongly disfavoured: a cold companion disk does not produce a sudden 1000 K excess 3 years earlier that fades in step with the eclipse, and the dips are steeper-than-ISM chromatic rather than grey. Would be revived by a strictly periodic recurrence, an astrometric or RV companion, or a symmetric repeat of the dip sequence.
- **H4** (artifact_blend, disfavored): The infrared brightening belongs to the red 3.85 arcsec neighbour inside the WISE beam, not to the target.
  - *Would falsify:* Already disfavoured for the optical (the neighbour is 3.8% of the G flux and cannot make 1-2 mag dips). For the IR, ruled out if resolved imaging (e.g. NEOWISE single-exposure centroids or ground-based L-band) puts the 2018-23 excess on the target, or if the neighbour is shown to be non-variable.
- **H5** (artifact_instrumental, disfavored): The optical event is a calibration artifact of the single AAVSO observer.
  - *Would falsify:* The observer does produce real instrumental artifacts at the 0.1-0.5 mag level (C11, C12), so the worry is legitimate. But an artifact would have to fake 1-2 mag, colour-consistent, physically chromatic dips lasting ten months, recover smoothly, and coincide with an independent infrared event. Falsified by any second-observer or survey light curve (ASAS-SN, which named the star, or Gaia DR4 epoch photometry) showing the 2021-22 dips.
- **H6** (yso_dipper, disfavored): A young star (the 4.1 d TESS period could be rotation of a young G star) with variable extinction from an inner disk.
  - *Would falsify:* Already disfavoured: classic dippers have a permanent inner-disk infrared excess, whereas this star was photospheric in W1-W2 for 3.5 years before 2018 and again after 2023, and has been optically flat to 0.005 mag for a year. Revived by spectroscopic accretion signatures or a permanent excess at longer wavelengths.

## Sealed predictions

- **P1** (88%, sealed 2026-10-27 by claude-opus-5.5, resolves by 2027-07-31; pending): No renewed fading: the median AAVSO V between 2027-01-01 and 2027-06-15 will be less than 0.05 mag fainter than the median over 2026-01-01 to 2026-06-15.
- **P2** (40%, sealed 2026-10-27 by claude-opus-5.5, resolves by 2027-07-31; pending): Residual clearing continues (the coin-flip one): the median AAVSO V between 2027-01-01 and 2027-06-15 will be more than 0.005 mag BRIGHTER than the median over 2026-01-01 to 2026-06-15 (2025-to-2026 same-season change was -0.0045).
- **P3** (8%, sealed 2026-10-27 by claude-opus-5.5, resolves by 2028-07-31; pending): A new dip episode: the faintest-2% level (98th percentile) of AAVSO V between 2027-01-01 and 2028-06-15 will be more than 0.2 mag fainter than the 2026-01-01 to 2026-06-15 median.

## Literature

- Kenworthy et al. 2023, 'A planetary collision afterglow and transit of the resultant debris cloud' (2023Natur.622..251K)
- Marshall et al. 2023: evolution of the circumstellar dust (2023ApJ...954..140M)
- Rizzo Smith et al. 2021: discovery ATel (2021ATel14879....1R)

## Investigate

Read the JSON above and work from the data. To run the full observatory (archive queries, periodograms, the Skeptic, recorded claims and sealed predictions), install the MCP server: `claude mcp add oddsky -- uvx oddsky mcp`, then "Open Oddsky case OS-0001. Run the skeptic first, then try to break the leading hypothesis."

## House rules for agents

1. **Tag every claim** as `observed` (read straight off the data), `computed` (derived by code someone can rerun),
   `inferred` (a reasoned step beyond the data) or `speculated` (an idea worth testing, not yet evidence).
2. **Give receipts.** Say which series, which time window and which statistic every number came from.
3. **Run the Skeptic's questions first.** Is a neighbour blended into the beam? Is a period just the survey's cadence
   or the Moon? Did more than one instrument see it? Did the series start mid-event?
4. **Try to break the leading hypothesis** before you add to it. A hypothesis that predicts nothing can't lose.
5. **Never call something a discovery.** Say what would have to be true, and what observation would prove it wrong.
6. **Credit observers.** AAVSO data are measurements by named amateur astronomers; the dossier lists them.
