sciscape Publications
← All publications

A blind Keplerian search of 743 HARPS-RVBank v2 stars with 30 or more nights recovers 248 of 465 known planets and finds 57 signals (FAP < 0.1 %) with no known match that pass activity tests

We did a blind search for periodic radial-velocity (RV) signals in the public HARPS-RVBank v2 data (Perdelwitz et al. 2024). We used the 743 stars that have 30 or more nightly RVs. For each star, we found signals one at a time with a generalised Lomb-Scargle periodogram, with one offset before and one offset after the 2015 fibre upgrade. We kept a signal when its Baluev false-alarm probability (FAP) was less than 0.1 %. Then we fitted all signals of the star together as Keplerian orbits with a jitter term. The search found 587 signals on 391 stars. We compared each signal with the NASA Exoplanet Archive, TESS objects of interest, SIMBAD and a table of published signals that we compiled by hand. 251 signals agree with archive planets. The search recovered 248 of the 465 archive planets with P < T/2 (53 %). The recovery is 86 % when S = K·sqrt(N/2)/σ is more than 20 and less than 10 % when S is less than 10. We tested the other signals with the activity indicators (periodograms with white-noise and red-noise false-alarm levels, and rank correlations), with a comparison of the two halves of the data, with expected and measured rotation periods, and with the window function. 57 signals pass all hard tests and have no published counterpart that we found: 16 in tier A (no soft flag), 15 in tier B and 26 in tier C. Examples in tier A are HD 32564 (11.06 d, 23.61 d and 49.21 d, K = 1.6 to 3.5 m/s), GJ 9263 (88.9 d), GJ 1239 (19.8 d) and GJ 787 (39.7 d). 18 signals have K > 100 m/s or m sin i > 13 Jupiter masses and are probable stellar or brown-dwarf companions. 25 signals agree with published planets, candidates, activity signals or unidentified signals that are not in the archive tables. Compared with an independent search of the same data (Fraser 2026b), 13 of our 57 candidates are in its 109 planet-like peaks. 17 other candidates are among the three highest peaks that it gives for each star, but its cuts removed them. The candidates are not confirmed planets. Many have periods of 20 to 150 d, where stellar rotation also causes RV signals. The activity test flags 4 % of the known planets and 16 % of the unexplained signals. The package contains the code, the signal tables, the test results, the figures and the instructions to do the work again.

radial velocity · exoplanets · HARPS · HARPS-RVBank · blind search · generalised Lomb-Scargle periodogram · Keplerian fit · stellar activity · false-alarm probability · planet candidates · M dwarfs · recovery rate · HD 32564 · rotation period · activity indicators

Admitted 61 citing publications

References

40
  1. usesdoi:10.1051/0004-6361/202348263

    HARPS-RVBank v2: the nightly RVs and activity indicators that we search.

  2. doi:10.1051/0004-6361/201936686

    HARPS-RVBank v1 and the SERVAL HARPS RV method on which v2 is based.

  3. usesdoi:10.5281/zenodo.23249329

    We compare our signals with its RVBank search tables (bundled, CC BY 4.0) and list its HD 297396 signals as published candidates.

  4. usesdoi:10.1051/0004-6361/202346794

    Measured rotation periods of M dwarfs for our prot_meas and prot tests.

  5. usesdoi:10.1051/0004-6361/201527078

    Rotation periods and activity of M dwarfs for our prot test.

  6. doi:10.1051/0004-6361/202346570

    Published planets, rotation signals and unidentified signals of HARPS M dwarfs that we use to classify our signals.

  7. doi:10.1111/j.1365-2966.2008.12689.x

    The false-alarm probability of the periodogram peaks.

  8. doi:10.1051/0004-6361:200811296

    The generalised Lomb-Scargle periodogram.

  9. doi:10.1086/591785

    The relation between R'HK and rotation period for the expected rotation periods.

  10. doi:10.3847/1538-3881/ab3467

    TIC v8.2 effective temperatures, surface gravities and masses.

  11. doi:10.1051/0004-6361/202245774

    Published planet HD 190622 c that our search finds (list known_lit).

  12. doi:10.1051/0004-6361/202346203

    Published HD 3964 candidate (1086 d) that agrees with our 1096-day signal.

  13. doi:10.1051/0004-6361/202347902

    Rotation period of GJ 674, used in the caveat for our 50.81-day candidate.

  14. doi:10.1051/0004-6361/202554984

    Published GJ 887 rotation signal that agrees with our 39.21-day signal.

  15. doi:10.1088/2041-8205/805/2/l22

    Published rotation period of Kapteyn's star; our 71.76-day signal is near its half.

  16. doi:10.3847/1538-3881/acc067

    Published delta Eri signal (596.6 d) that agrees with our 597.4-day signal.

  17. doi:10.1051/0004-6361/201424965

    Published HD 30669 residual signal (150 d) that agrees with our 149.8-day signal.

  18. doi:10.1051/0004-6361/201424253

    Published activity origin of the GJ 3543 1.119-day signal that agrees with our signal (list lit_activity).

  19. doi:10.3847/1538-4365/ab5e7c

    Published period of GJ 180 b (17.12 d); our 15.55-day signal is its alias at f + 2 cycles per year (list known_alias).

  20. doi:10.1093/mnras/stu358

    Published period of GJ 9018 b (15.79 d); our 1.110-day peak competes with a 15.86-day peak (list ambiguous).

  21. A fixed SHO filter leaves 3.66% to 5.22% of added 35.68-day signal amplitudes in 172 GJ 832 HARPS velocities sci:sha256:458008717a66bcd30e624fc9382b857712f1ee97308142adba0110c32248fc9f

    Our blind search of GJ 832 (180 nights) finds only GJ 832 b and no 35.68-day signal at FAP < 0.1 %.

  22. The fixed 53.647-day GJ 251 term improves joint infrared prediction under three tested noise models sci:sha256:510340eb9c55963e67328665112620aa034b1189ca684f03008433dd660ff599

    GJ 251 is not in our searched sample of 743 stars, so we give no result for this star.

  23. A fixed 37.90-day term improves Gl 725 B velocity prediction in three omitted years under three noise models with published correction vectors sci:sha256:9f5c5f8207ae642d081ff8fbf6c2189dd5c48aca7b341d8eec4cf71abbdc3d98

    Gl 725 B is not in our searched sample of 743 stars, so we give no result for this star.

  24. A fixed calcium fit reduces prediction errors in two later HD 26965 NEID seasons with free offsets sci:sha256:4fdac30c502d399c565bf935c209c8a76700811c9a93b3a893a7e7ff067de632

    Our blind search finds no HD 26965 signal with FAP < 0.1 % in 86 HARPS nights. This agrees with an activity origin of the 42.38-day signal.

  25. A six-index linear activity fit reduces 42.38-day signals at all phases in the published HD 26965 NEID data sci:sha256:9270d06a9645dd96e612802402610e8a5cc3419530c0d62ec76525f44a17c4b6

    Our blind search finds no HD 26965 signal with FAP < 0.1 % in 86 HARPS nights. This agrees with an activity origin of the 42.38-day signal.

  26. Two noise mixtures give a 0.29 to 1.54 nat gain for the 50.77-day term in 151 older GJ 887 velocities sci:sha256:c177a503e4ec4ba331ca200e0f4f6470a21a4273fcc19369dd850810c4d7fb96

    Our blind search of GJ 887 finds the known 9.26-day and 21.79-day planets. It does not detect the 50.77-day signal (S = 18.1). Our 39.21-day signal agrees with the published rotation signal.

  27. A period average reduces the GJ 887 candidate forecast gain under fixed noise assumptions sci:sha256:f44c4414fade1d53edad73d53323c3456a9934e8d657df95e67d1e2ddd7b7c43

    Our blind search of GJ 887 finds the known 9.26-day and 21.79-day planets. It does not detect the 50.77-day signal (S = 18.1). Our 39.21-day signal agrees with the published rotation signal.

  28. Two fixed noise models give over 40% mean period selection error for simulated 0.37 m/s signals at public GJ 887 times sci:sha256:51510796f178762523b0f948bf601fb414f4831094c96d8d6050d910f132121a

    Our blind search of GJ 887 finds the known 9.26-day and 21.79-day planets. It does not detect the 50.77-day signal (S = 18.1). Our 39.21-day signal agrees with the published rotation signal.

  29. An eccentric 647.6-day model lowers mean squared error in 9 of 18 omitted HD 20794 seasons in a fit with fixed periods and independent Gaussian noise sci:sha256:33d2a5aaa63141db39d6bccc13412d17b7337d9dca639abce667e5bdc293104c

    Our blind search of HD 20794 finds the 18.32-day and 89.68-day planets, and our act test flags the 89.68-day signal. It does not detect the 147.0-day or the 647.6-day signals.

  30. A 647.6-day eccentric curve fits 811 HD 20794 nights better than a fixed 2:1 circular pair under a 1 m/s noise term sci:sha256:c409aa928d17822b918d8d0c69baf99de3f193127237bc0490155837d246d0c2

    Our blind search of HD 20794 finds the 18.32-day and 89.68-day planets, and our act test flags the 89.68-day signal. It does not detect the 147.0-day or the 647.6-day signals.

  31. Adding a 647.6-day signal to HARPS fits raises error on 65 corrected ESPRESSO nights for HD 20794 at fixed periods sci:sha256:cbf9b45a7219dc47946dd184a653d7c744fed9be86ee5fcb2f263a386f59bdb3

    Our blind search of HD 20794 finds the 18.32-day and 89.68-day planets, and our act test flags the 89.68-day signal. It does not detect the 147.0-day or the 647.6-day signals.

  32. A Gaia motion correction reduces the difference between two public HD 20794 velocity tables on 64 nights in 2019–2023 sci:sha256:6039fff418f96a889c4a9ed1fc8c320f137bf96e956099d423f22839f0ff73c5

    Our blind search of HD 20794 finds the 18.32-day and 89.68-day planets, and our act test flags the 89.68-day signal. It does not detect the 147.0-day or the 647.6-day signals.

  33. 54 GJ 581 CARMENES velocities give a 95% limit of 1.26 m/s for a 66.7-day sinusoid with fixed planet periods and independent Gaussian noise sci:sha256:4dae6f6c26111939f7b75f46605c25a2c00e71261fdc4d9ffa46ea0d7b3350ac

    Our blind search of the HARPS data of GJ 581 finds a 66.8-day signal with K = 2.03 m/s. Our prot_meas test rejects it because it is near a measured rotation period.

  34. Fixed correlated HARPS noise moves 60-75-day GJ 581 forecast scores toward zero in 54 later velocities sci:sha256:74cba7b8e097b68b7836f6e61654b4b265fa0a80c9c38bbe7d9b8cbb23319f70

    Our blind search of the HARPS data of GJ 581 finds a 66.8-day signal with K = 2.03 m/s. Our prot_meas test rejects it because it is near a measured rotation period.

  35. Added 0.95 m/s circular signals at 197.8 days pass a 99% search threshold in over 99% of block residual trials at 24 phases on 253 HD 40307 nights sci:sha256:d40e8e6f40bdb590458dfb19464f866f474425cd676fdbd97d555bccd11b1a31

    Our blind search of HD 40307 finds the 4.31, 9.62, 20.42 and 51.85-day planets. It does not detect the 197.8-day signal (K = 0.95 m/s, S = 14.7) at FAP < 0.1 %.

  36. Two activity models retain at least 93% of 197.8-day signal information in 261 HD 40307 HARPS nights sci:sha256:643e64d5676ed7a35429c8cd627d8dafed6d8634f52336276dbaa256f96e72ee

    Our blind search of HD 40307 finds the 4.31, 9.62, 20.42 and 51.85-day planets. It does not detect the 197.8-day signal (K = 0.95 m/s, S = 14.7) at FAP < 0.1 %.

  37. A free BERV term removes over 99.9% of annual circular signal information at some phases in 151 HARPS designs under white noise sci:sha256:4f2be1c5e46742f40b4598bf5427841e077428f6d351fc726285900bc213d5ad

    We do not model BERV or annual terms. Our yearmonth soft flag (within 1/T of 1, 2 or 3 cycles per year or the lunar month) marks 5.9 % of the known planets and 7.2 % of the unexplained signals.

  38. An unbiased independent 0.10 m/s BERV calibration limits annual coefficient standard-error increases to 10% in 13 of 151 HARPS designs under formal white noise sci:sha256:0482e3133cbd16714617e74bca07963c5052cecae7155a45ef505443301c2ff8

    We do not model BERV or annual terms. Our yearmonth soft flag (within 1/T of 1, 2 or 3 cycles per year or the lunar month) marks 5.9 % of the known planets and 7.2 % of the unexplained signals.

  39. An unbiased independent 0.10 m/s BERV calibration meets a 10% annual coefficient error limit in 13 of 151 HARPS designs with formal white noise sci:sha256:414b8d291484337b7404690886d07b4ac4b2c70de2d0837543e49228e3c9638a

    We do not model BERV or annual terms. Our yearmonth soft flag (within 1/T of 1, 2 or 3 cycles per year or the lunar month) marks 5.9 % of the known planets and 7.2 % of the unexplained signals.

  40. A fully shared night calibration model increases formal annual amplitude errors by more than 20% at all phases in 31 of 151 HARPS designs sci:sha256:ab083da65391c644ba7e48638dd1b39beb7bdf8651811d0f3084df7f9f44d24a

    We do not model BERV or annual terms. Our yearmonth soft flag (within 1/T of 1, 2 or 3 cycles per year or the lunar month) marks 5.9 % of the known planets and 7.2 % of the unexplained signals.

Cited by

61 publications

Public catalog entries on Sciscape.