Rotation Periods of Earth-approaching Main-Belt Asteroids

Study Goal

  • Determine rotation periods for 14 main-belt asteroids with Earth–MOID < 1.1AU1.1\,\text{AU} using sparse, uneven brightness data.

Core Data Sources

  • Minor Planet Center (MPC) brightness measurements (multi-observatory, multi-filter).
  • Additional Baldone Observatory data.
  • Validation sets from TESS, ATLAS, Catalina, Pan-STARRS, etc.

Essential Pre-processing

  • Select observatory/filter sets with >7070 measurements (total >10001000 for an asteroid).
  • Brightness correction for changing distances:
    m<em>ri=m</em>i5log(Δ<em>iΔ</em>0r<em>0r</em>i)m<em>{ri}=m</em>i-5\log\left(\frac{\Delta<em>i}{\Delta</em>0}\frac{r<em>0}{r</em>i}\right)
  • Light-time correction:
    JD<em>ri=JD</em>i+Δ<em>0Δ</em>icJD<em>{ri}=JD</em>i+\frac{\Delta<em>0-\Delta</em>i}{c}
  • Remove points with phase <7<7^\circ and >35>35^\circ.
  • Three-sigma clipping per night; apply linear phase correction:
    m<em>ϕ,i=m</em>ri(ϕ<em>iϕ</em>0)kjm<em>{\phi,i}=m</em>{ri}-(\phi<em>i-\phi</em>0)k_j

Lomb–Scargle Analysis

  • Frequency search range: 0.1!!50day10.1!\text{–}!50\,\text{day}^{-1} (periods 0.48!!240h0.48!\text{–}!240\,\text{h}).
  • Use Astropy L-S with 10001000 trial frequencies; refine around peaks (power w>0.2).
  • Data-set quality metric:
    R2=1N1<em>i=1N(SS</em>v)2σiR^2=\frac{1}{N-1}\sum<em>{i=1}^{N}\frac{(S-S</em>v)^2}{\sigma_i} (require R^2>0.5).

Period Selection Rules

  • Accept peaks showing:
    • Bimodal light curve (2 maxima & minima).
    • Power-spectrum peak w>0.38, Gaussian-like.
  • Combine multiple observatories via weighted average:
    P=(w<em>i/w</em>s)p<em>i+(N</em>i/N<em>s)p</em>i+(R<em>i2/R</em>s2)pi3P=\frac{\sum (w<em>i/w</em>s)p<em>i+\sum (N</em>i/N<em>s)p</em>i+\sum (R<em>i^2/R</em>s^2)p_i}{3}
  • Weighted error:
    σ<em>kv=(p</em>ipˉ)2(M1)M\sigma<em>{kv}=\sqrt{\frac{\sum(p</em>i-\bar p)^2}{(M-1)M}}

Key Results (Weighted Periods)

  • Newly determined (first time):
    • (1779) Parana 22h\approx22\,\text{h} (uncertain),
    • (1818) Brahms 5.357h5.357\,\text{h},
    • (2128) Wetherill 19.748h19.748\,\text{h},
    • (2318) Lubarsky 5.421h5.421\,\text{h},
    • (2497) Kulikovskij 77.750h77.750\,\text{h},
    • (2503) Liaoning 103.042h103.042\,\text{h},
    • (2538) Vanderlinden 53.422h53.422\,\text{h},
    • (2539) Ningxia 9.790h9.790\,\text{h},
    • (2583) Fatyanov 7.815h7.815\,\text{h}.
  • Previously known periods re-confirmed:
    • (1951) Lick 5.306h5.306\,\text{h},
    • (1963) Bezovec 18.181h18.181\,\text{h},
    • (2134) Dennispalm 4.113h4.113\,\text{h},
    • (2150) Nyctimene 6.125h6.125\,\text{h}.
  • Possible alternative shorter period for (2174) Asmodeus 4.785h4.785\,\text{h} (previous 41.05h41.05\,\text{h}).

Method Limitations

  • Ineffective for spin axes nearly perpendicular to orbital plane (low light-curve amplitude).
  • Assumes simple bimodal curves; fails on multimodal, binaries, tumblers.
  • Linear phase–magnitude assumption valid only for 7!!407^\circ!\text{–}!40^\circ; bias <0.19%0.19\% when w>0.38w>0.38.
  • Requires S/N >5 (used 0.1mag0.1\,\text{mag} accuracy ⇒ S/N 10\ge10).

Practical Checklist

  • Ensure filtered data counts (per set>7070, total>10001000).
  • Apply distance, light-time, phase corrections + 3σ3\sigma clipping.
  • Use L-S; accept only peaks meeting w>0.38 & bimodal test.
  • Cross-validate across ≥2 datasets; compute weighted PP & σkv\sigma_{kv}.

Takeaways

  • Lomb–Scargle + rigorous filtering recovers reliable periods from sparse MPC photometry.
  • Nine asteroid periods reported for the first time; four literature values confirmed; one alternate period suggested.
  • Method best for elongated, single-body asteroids with appreciable light-curve amplitude.