As my loyal reader knows, our radial-velocity companion-search code called The Joker is a brute-force least-square fitting (plus some Bayes) of an elliptical orbit at every possible period, orientation, and phase. In astrophysics, the Lomb–Scargle periodogram is a workhorse tool that, under the hood, is a brute-force least-square fitting of a sinusoid at every possible period and phase. So these two ideas are fundamentally incredibly similar. And indeed, today Winston Harris (MTSU) demonstrated quantitatively that these are the same, and will become identical as eccentricities go to zero. That's interesting, because if we are doing companion search and we don't mind making the approximation that eccentricities are zero, the fitting of sinusoids (even at arbitrary phase, because: trig identities) is way, way faster than the fitting of Kepler functions.
I think the "Bayesian Kepler periodogram"that Tom Loredo and Phil Gregory have developed are relevant and perhaps better tools:
ReplyDeleteA Bayesian Kepler periodogram detects a second planet in HD 208487
P. C. Gregory
Monthly Notices of the Royal Astronomical Society, Volume 374, Issue 4, February 2007, Pages 1321–1333, https://doi.org/10.1111/j.1365-2966.2006.11240.x
Interested in the complex Fourier transform, not just the periodogram? I have code ... not quite brute force, but straightforward ... (jeffscargle@gmail.com).
ReplyDelete