Fixed correlated HARPS noise moves 60-75-day GJ 581 forecast scores toward zero in 54 later velocities
We test whether the disputed GJ 581 velocity signal continues in a later instrument. Four circular terms are trained on all 250 released HARPS rows. All 54 later CARMENES rows are scored without fitting their signal parameters. Orthonormal contrasts remove the unknown instrument offset. The forecast includes uncertainty in coefficients, frequencies, and added HARPS noise. Both hypotheses use the same HARPS model. They differ only in whether one trained term continues. The primary white-noise model gives a log predictive density change of -6.25 nats for the 66.8-day term. Fixed 133-day and 10-day correlated HARPS noise move this score to -1.00 and -2.44 nats with the same test noise. Four training models and six fixed test covariances give 24 comparisons for this term. Most scores are adverse. One weak positive result gives +0.675 nat in a 32000-draw check. A secondary period mode near 71 days has 13 percent of the training posterior but supplies 95 percent of its full forecast density. The three established planet signals improve all 72 control comparisons. This shows that the adverse forecast for the disputed term depends on the noise model. The package contains exact public data, source hashes, code, posterior draws, all scores, integration checks, simulations, figures, and reproduction steps. This is a retrospective conditional forecast test. It does not determine whether GJ 581 d exists or prove an activity origin. Epoch and instrument change together. The fixed noise grid does not cover all stellar noise models.
GJ 581 · radial velocity · predictive density · signal continuation · stellar activity · HARPS · CARMENES · public astronomical data
References
9- usesdoi:10.1051/0004-6361/202449375
This paper releases the HARPS and CARMENES velocity tables in CDS catalog J/A+A/688/A112. It supplies the data and scientific context for the forecast test.
- usesdoi:10.26093/cds/vizier
The CDS service hosts the exact public source tables. This is the service DOI, not a version-specific dataset DOI.
- usesA period average reduces the GJ 887 candidate forecast gain under fixed noise assumptions sci:sha256:f44c4414fade1d53edad73d53323c3456a9934e8d657df95e67d1e2ddd7b7c43
This GJ 887 study informs coefficient integration and period uncertainty in a velocity forecast. We study a different star and test signal continuation. No code is copied.
- usesTwo 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
This GJ 887 study informs joint instrument forecasts and offset removal with orthonormal contrasts. It also informs planet controls and conditional simulations.
- doi:10.1126/science.1253253
This paper gives an activity interpretation of the disputed 66.7-day GJ 581 signal.
- doi:10.1126/science.1260796
This comment questions residual analysis in the earlier activity interpretation. The present study fits all training velocity terms at the same time.
- doi:10.1093/mnras/sts476
This paper shows that correlated noise affects weak GJ 581 signals. It motivates the separate 10-day exponential noise check.
- doi:10.3847/2515-5172/ad1de4
This paper reports a photometric rotation period near 149 days. It differs from the 133-day period used in one check. We do not resolve that difference.
- usesdoi:10.1086/670067
The emcee ensemble sampler is used for the training frequency and noise posterior.
Cited by
2 publicationsPublic catalog entries on Sciscape.
- A 41.78-day term raises the total prediction score for omitted years in 204 GJ 3998 velocities under fixed noise
Informs the use of joint prediction scores, coefficient uncertainty, and stated fixed-noise checks. No code or data are reused.
- 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
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.