Open Access
Issue
A&A
Volume 693, January 2025
Article Number L3
Number of page(s) 7
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202452450
Published online 23 December 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Exoplanetary astronomy has had a remarkable expansion in recent years. The number of confirmed exoplanets will soon reach 5800 (NASA Exoplanet Archive, Akeson et al. 2013; accessed November 2024). Astronomers are already on track to discover potentially habitable Earth-like planets (e.g. Crossfield et al. 2015; Gillon et al. 2016; Murgas et al. 2023; Suárez Mascareño et al. 2023). The discovery of a planet in the habitable zone of Proxima, in the closest stellar system to our own (Proxima b; Anglada-Escudé et al. 2016), rocked the scientific community. It hinted that potential Earth analogues may lie in the very vicinity of the solar neighbourhood. Subsequently, considerable effort has been devoted to the analysis of these nearby systems. This led to a natural focus on the second-closest stellar system: Barnard’s star (Table 1), an M3.5V-M4V dwarf that has been thoroughly examined since its discovery (Barnard 1916).

Table 1.

Relevant parameters of Barnard’s star.

More recently, González Hernández et al. (2024), hereafter GH24, reported the detection of a sub-Earth-mass planet around Barnard’s star. This planet, designated as Barnard b, has an orbital period Pb = 3.1533 ± 0.0006 d and a minimum mass of Mbsinib = 0.37 ± 0.05M. The proximity of Barnard b to its host star (0.02 au) implies a transit probability of 3.7%, which in turn would offer an extraordinary opportunity for its atmospheric characterisation. The Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015) had obtained photometric observations of Barnard’s star between June and July 2024, and released those data as part of the Sector 80 (S80) data release in late August 2024, a mere few days after the acceptance of GH24. The present Letter serves to extend their work by exploring potential transits of Barnard b in this new and only TESS photometry.

2. Data

The TESS Science Processing Operations Center pipeline (SPOC; Jenkins et al. 2016) provides two different photometric data products: Simple Aperture Photometry (SAP) and Pre-Search Data Conditioning Simple Aperture Photometry (PDCSAP). SAP is a simple summation of flux from the optimal aperture of the target determined by SPOC. Its data may suffer from systematics of different kinds, including contamination from other sources, cosmic rays, variable incident solar radiation, foreign particles on CCD sensors, and systematics from spacecraft motion. With a similar rationale, but regarding the Kepler mission, Kinemuchi et al. (2012) argued that SAP data should be expected to carry all these systematics; analysis of subtle photometric variations on timescales larger than 1 d is likely to be affected. PDCSAP, on the other hand, aims to mitigate signal artefacts and instrumental effects in SAP data. The basic mode of operation is described in Sect. 4.2 of the same work, and has been further refined since then (see Stumpe et al. 2012; Smith et al. 2012; Stumpe et al. 2014). We discuss SAP and PDCSAP data in all the analysis that follows.

We studied the origin of flux inside the SPOC aperture of Barnard’s star through the TESS contamination tool TESS-CONT (Castro-González et al. 2024). We analysed the TESS pixel-response functions of all sources within 500 arcsec (N = 4819). We found that 99.5% of the aperture flux indeed comes from Barnard’s star. The remaining 0.5% are ascribed to fainter sources, whose small contributions add up to a negligible contamination.

3. Analysis

Two parameters define the general shape of a transit: the transit width W from first to fourth contact, and the transit depth D. Deriving some characteristic values of W and D, even if approximate, would give us clues when we look for features in photometry. Both parameters, but especially D, are sensitive to the planetary radius Rb, which in the case of Barnard b is not directly probed in GH24. Nevertheless, we can relate their measured Mbsinib to Rb for an assumed mass-radius relation and an assumed impact-parameter distribution. We used the mass-radius relation provided by SPRIGHT (Parviainen et al. 2024) to estimate Rb under the assumption sinib ≈ 1. For an impact-parameter distribution, we adopted b = 𝒰(0, 1), which would be implied by a random orbital orientation relative to the line of sight.

We then set out to estimate W and D through rapid transit curve modelling. We assumed that the limb-darkening profile of Barnard’s star follows the quadratic law

I ( ϑ ) I ( ϑ = 0 ° ) = 1 u 1 ( 1 cos ϑ ) u 2 ( 1 cos ϑ ) 2 , $$ \begin{aligned} \frac{I(\vartheta )}{I(\vartheta =0^{\circ })}=1-u_1(1-\cos \vartheta )-u_2(1-\cos \vartheta )^2, \end{aligned} $$(1)

where I is the total emergent intensity, and ϑ is the angle between the line of sight and the direction of the emergent radiation (Kopal 1950). We fed stellar-parameter measurements from Table 1 to LDTK (Parviainen & Aigrain 2015; Husser et al. 2013), which returned u1 = 0.2470 ± 0.0013 and u2 = 0.3814 ± 0.0026 for the TESS passband. We then computed 106 transit curves from the distributions of (i) R, M, Pb, Mb from GH24; (ii) Rb from SPRIGHT; (iii) u1, u2 from LDTK, assumed Gaussian; and (iv) b = 𝒰(0, 1). We calculated W and D of each light curve individually; we list their aggregate statistics together with our derived Rb in Table 2. The bulk of the distribution of Rb appeared Gaussian-like, with few outliers beyond 1 R. We expect a transit of depth 1 . 55 0.28 + 0.26 $ 1.55^{+0.26}_{-0.28}\, $ppt (parts per thousand) and of width 0 . 82 0.27 + 0.11 $ 0.82^{+0.11}_{-0.27}\, $h. GH24 provided with an additional four-planet solution that contains Barnard b and three planetary candidates. We repeated the same procedure for those solutions; we list the corresponding results in Table A.1.

Table 2.

Relevant parameters of Barnard b.

3.1. Lack of features in TESS S80 data

We took all 2 min SAP and PDCSAP measurements from the TESS S80 lightcurve, using the default cadence quality bitmask. The SAP dataset is comprised of 16 492 measurements in the range 2460479–2460507 BTJD. SAP measurements have a root mean square (rms) of 166 es−1 and a median uncertainty of 52 es−1. We show these time series in Fig. 1a, together with predicted inferior conjunctions of Barnard b from GH24 emphemerides. SAP measurements are split into several portions, and contain irregular negative trends that tend to take place at the end of those portions (e.g. near 2460485 BTJD or 2460493 BTJD). However, there is no apparent evidence of a strictly periodic decrease in flux, which would otherwise be caused by the sought planetary transit. The PDCSAP time series share the same temporal coverage, but with 4691 excluded measurements in the interval 2460485–2460492 BTJD (Fig. 1b). This exclusion affected the coverage of two predicted inferior conjunctions. PDCSAP time series contain 11 801 measurements that have a flux rms of 78.4 es−1 and a median uncertainty of 54.5 es−1.

thumbnail Fig. 1.

TESS S80 photometric data products: (a) SAP; (b) PDCSAP, with superimposed predicted times of inferior conjunction through GH24 ephemerides. The propagated inferior-conjunction uncertainty is shown as light blue bands.

We proceeded with SAP in an attempt to include data near the inferior-conjunction predictions that PDCSAP dismissed. We smoothed the SAP time series with several variants of detrending, including second-order Savitzky-Golay filters (Savitzky & Golay 1964), Tukey’s biweight time-window sliders (Mosteller & Tukey 1977), and mean time-window sliders, all available in Wōtan (Hippke et al. 2019). We globally used a window length of 0.5 d, an edge cut-off of 0.3 d, and a break tolerance of 0.2 d. For none of these forms of detrending could we find evidence of transits after generating transit least squares periodograms (TLSPs; Hippke & Heller 2019) or box least squares periodograms (BLSPs; Kovács et al. 2002) for periodicities in the range 1–15 d and transit widths in the range 0.2–2 h.

We chose to go forward with mean time-window detrending by virtue of its simplicity. Figure 2 displays this detrended SAP photometry. It contains 14 366 measurements that are characterised by a flux rms of 0.411 ppt and a median uncertainty of 0.222 ppt. The flux rms is nearly four times smaller than the expected transit depth ( 1 . 55 0.28 + 0.26 $ 1.55^{+0.26}_{-0.28} $ ppt; Fig. 2). Figure 2a compares these processed time series again with the predicted inferior conjunctions of Barnard b from GH24. Six inferior conjunctions are covered by photometry within the ephemeris precision. We could observe no significant decreases in flux at those ephemerides or anywhere else. We folded the time series at the published period and epoch of Barnard b, and binned the folded data at 5 min. Figure 2b illustrates the result of this exercise, together with the epoch uncertainty and PYTRANSIT curves from the median posteriors of Table 2 for three impact-parameter values: 0, 0.9, and 1. The measurements clearly do not follow predictions for such a liberal range of impact parameters, regardless of whether binning was involved or not. This remained the case after binning at other intervals (2–10 min). On the whole, these timeseries give an impression of notable symmetry and lack of apparent features, as if we were observing pure noise when measuring a constant-flux source. The distribution of measurements is Gaussian-like, with its right tail being asymmetrically heavier (Fig. A.1). On the other hand, transits must break the flux distribution symmetry in the opposite direction, i.e. in the left tail.

thumbnail Fig. 2.

The search of Barnard b in TESS S80 photometry. (a) Mean time-window detrended SAP photometry with superimposed predicted times of inferior conjunction through GH24 ephemerides. The median expected transit width is displayed as dark blue bands. The propagated ephemeris uncertainty is shown as light blue bands. (b) Phase-folded plot at GH24 ephemerides: before and after applying a 5 min binning (light hollow and dark filled circles, respectively). The inferior-conjunction epoch uncertainty is shown as a light blue band. We compare photometry against the expected transit curve of Barnard b for three impact-parameter values: b = 0.0 (black dotted line), b = 0.9 (black dashed line), and b = 1.0 (black solid line). (c,d) GLSP of unbinned data before and after iterative 4σ clipping. (e,f) BLSP of unbinned data before and after iterative 4σ clipping.

We looked for any out-of-transit variations through the generalised Lomb-Scargle periodogram (GLSP; Lomb 1976; Scargle 1982; Zechmeister & Kürster 2009; VanderPlas & Ivezić 2015) of the detrended SAP photometry. Figure 2c displays this exercise. The GLSP highlights two peaks near 0.32 d and 0.35 d, and assigns them a false-alarm probability (FAP; Baluev 2008) of 3% and 2%, respectively. Their periodicities do not appear to be related to the window length we used (0.5 d). Phase-folding at those periods did not reveal any apparent structures in photometry. Analysis in the spirit of Dawson & Fabrycky (2010) revealed that the common harmonics of the two peaks do not correspond to any significant peaks in the window function of the dataset, and therefore are unlikely to be aliases of one another. We found out that the iterative sigma-clipping of the SAP photometry suppressed these signals greatly (Figure 2d). Clipping within 5σ, 4σ, and 3σ resulted in GLSPs where the strongest peak had a FAP of 5%, 78%, and > 99%, respectively. Most of the rejected outliers tend to congregate around the first (2460480–2460485 BTJD) and third seasons (2460492–2460493 BTJD) in regions of higher flux dispersion. All of these observations, together with the demonstrated good behaviour of the flux distribution, are against the physical nature of the 0.32 d and 0.35 d signals. Figures 2e and 2f show the BLSPs of the detrended SAP photometry before and after sigma-clipping. In neither periodogram could we find any significant peaks that merit discussion.

3.2. Lack of plausible transit solutions

We performed transit-curve fitting through PYTRANSIT and nested sampling (Skilling 2004). For the inference of these models, we used the nested-sampling integrator provided by DYNESTY (Speagle 2020). In every inference instance of Nparam model parameters, we required 40Nparam live points and random-walk sampling. We modelled one planet on a circular orbit, using GH24 ephemerides as priors. We imposed two different priors on the planetary radius: one informed by the expected radius 𝒩(0.76, 0.04) R from Table 2, and one with the very uninformed U log ( e 5 , e + 5 ) R $ \mathcal{U}_{\mathrm{log}}(e^{-5},e^{+5})\,\mathrm{R}_{\oplus} $ (i.e. up to 148 R). We hereafter refer to these priors as the narrow and wide priors, respectively. Our impact-parameter prior was liberally set to 𝒰(0, 1.5). Quadratic limb-darkening coefficients took priors defined by the Gaussian distributions defined by LDTK, given the stellar parameters from Table 1. Their variances were scaled three-fold in the prior to ensure impartiality. All models sampled for a zero-order flux offset ΔΦ, which pinpoints the out-of-transit flux level relative to Φ = 1. Additionally, all models included a rigid jitter component of 0.25 ppt, based on the rms of the first 10 d of data.

Bayesian evidence of one-planet models reveals that the data do not support a planetary addition (ΔlnZ = −89 for narrow Rb priors, ΔlnZ = −143 for wide Rb priors). The narrow-prior model restricts impact-parameter solutions to b >  1 (3σ), which correspond to inclinations ib <  87.9° (Figure A.2). This suggests that the data are not compatible with any transiting events with the properties of Barnard b, including cases of grazing transits. Although the wide-prior model provided upper limits to potential transiting planets, we note a large degeneracy between the radius and the inclination. For ib >  88°, planetary radii up to a few tenths of the R were allowed. For less aligned transits, sensitivity decreased to the order of R, to such an extent that even 10 R planets would remain undetected for ib <  87° (Fig. A.3). We tested models with and without Gaussian processes (GPs; Rasmussen & Williams 2006), using a Matérn 5/2 kernel, as defined in S+LEAF (Delisle et al. 2020, 2022). We additionally included three more planets in the model to try to reach the four-planet solution in GH24. The inclusion of neither GPs nor additional planets had a meaningful effect.

4. Conclusion

We investigated the first and, to date, only TESS sector, S80, that contains photometry of Barnard’s star. PDCSAP and SAP photometry were shown to carry no apparent transit features within measurement precision. The flux distribution in detrended SAP photometry appears Gaussian-like. We attempted transit-curve fitting with informed and uninformed priors on the planetary radius, but this yielded no significant results. We applied a very similar methodology on detrended PDCSAP photometry of 20 s and of 2 min cadence; both proved equally unfruitful. On that account, we conclude that current TESS photometry strongly suggests that Barnard b is non-transiting. This implies a marginal constraint in orbital inclination of the planet (ib ≤ 87.9°; 3σ). We were unable to find any other periodic darkening in the photometry time series that could be ascribed to planetary transits.

Acknowledgments

AKS acknowledges the support of a fellowship from the “la Caixa” Foundation (ID 100010434). The fellowship code is LCF/BQ/DI23/11990071. AKS, JIGH, ASM, CAP, NN, and RR acknowledge financial support from the Spanish Ministry of Science and Innovation (MICINN) project PID2020-117493GB-I00 and from the Government of the Canary Islands project ProID2020010129. NN acknowledges funding from Light Bridges for the Doctoral Thesis “Habitable Earth-like planets with ESPRESSO and NIRPS”, in cooperation with the Instituto de Astrofísica de Canarias, and the use of Indefeasible Computer Rights (ICR) being commissioned at the ASTRO POC project in Tenerife, Canary Islands, Spain. The ICR-ASTRONOMY used for his research was provided by Light Bridges in cooperation with Hewlett Packard Enterprise (HPE). ACG is funded by the Spanish Ministry of Science through MCIN/AEI/10.13039/501100011033 grant PID2019-107061GB-C61. AMS acknowledges support from the Fundação para a Ciência e a Tecnologia (FCT) through the Fellowship 2020.05387.BD (DOI: 10.54499/2020.05387.BD). Funded/Co-funded by the European Union (ERC, FIERCE, 101052347). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. This work was supported by FCT – Fundação para a Ciência e a Tecnologia through national funds by these grants: UIDB/04434/2020, UIDP/04434/2020. The work of CJAPM was financed by Portuguese funds through FCT (Fundação para a Ciência e a Tecnologia) in the framework of the project 2022.04048.PTDC (Phi in the Sky, DOI 10.54499/2022.04048.PTDC). CJAPM also acknowledges FCT and POCH/FSE (EC) support through Investigador FCT Contract 2021.01214.CEECIND/CP1658/CT0001 (DOI 10.54499/2021.01214.CEECIND/CP1658/CT0001). This paper includes data collected by the TESS mission, which are publicly available from the Mikulski Archive for Space Telescopes (MAST). Funding for the TESS mission is provided by NASA’s Science Mission Directorate. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. We acknowledge the use of public TESS data from pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. This research has made use of the Exoplanet Follow-up Observation Program website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. This work made use of TESS-CONT, which builds upon TPFPLOTTER (Aller et al. 2020) and TESS_PRF (Bell & Higgins 2022). We used the following PYTHON packages for data analysis and visualisation: ASTROPY (Astropy Collaboration 2022), DYNESTY (Speagle 2020), LIGHTKURVE (Lightkurve Collaboration 2018), LDTK (Parviainen & Aigrain 2015; Husser et al. 2013), MATPLOTLIB (Hunter 2007), NIEVA (Stefanov et al., in prep.), NUMPY (Harris et al. 2020), PANDAS (McKinney 2010; The Pandas Development Team 2022), PYASTRONOMY (Czesla et al. 2019), PYTRANSIT (Parviainen 2015), SCIPY (Virtanen et al. 2020), SEABORN (Waskom 2021), S+LEAF (Delisle et al. 2020, 2022), SPRIGHT (Parviainen et al. 2024) and Wōtan (Hippke et al. 2019). This manuscript was written and compiled in OVERLEAF.

References

  1. Akeson, R. L., Chen, X., Ciardi, D., et al. 2013, PASP, 125, 989 [Google Scholar]
  2. Aller, A., Lillo-Box, J., Jones, D., Miranda, L. F., & Barceló Forteza, S. 2020, A&A, 635, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Anglada-Escudé, G., Amado, P. J., Barnes, J., et al. 2016, Nat, 536, 437 [Google Scholar]
  4. Astropy Collaboration (Price-Whelan, A. M., et al.) 2022, AJ, 935, 167 [NASA ADS] [CrossRef] [Google Scholar]
  5. Baluev, R. V. 2008, MNRAS, 385, 1279 [Google Scholar]
  6. Barnard, E. E. 1916, AJ, 29, 181 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bell, K. J., & Higgins, M. E. 2022, Astrophysics Source Code Library [record ascl:2207.008] [Google Scholar]
  8. Castro-González, A., Lillo-Box, J., Armstrong, D. J., et al. 2024, A&A, 691, A233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Crossfield, I. J. M., Petigura, E., Schlieder, J. E., et al. 2015, AJ, 804, 10 [NASA ADS] [Google Scholar]
  10. Czesla, S., Schröter, S., Schneider, C. P., et al. 2019, Astrophysics Source Code Library [record ascl:1906.010] [Google Scholar]
  11. Dawson, R. I., & Fabrycky, D. C. 2010, ApJ, 722, 937 [Google Scholar]
  12. Delisle, J. B., Hara, N., & Ségransan, D. 2020, A&A, 638, A95 [EDP Sciences] [Google Scholar]
  13. Delisle, J.-B., Unger, N., Hara, N. C., & Ségransan, D. 2022, A&A, 659, A182 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Gillon, M., Jehin, E., Lederer, S. M., et al. 2016, Nature, 533, 221 [Google Scholar]
  16. González Hernández, J. I., Suárez Mascareño, A., Silva, A. M., et al. 2024, A&A, 690, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
  18. Hippke, M., & Heller, R. 2019, A&A, 623, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Hippke, M., David, T. J., Mulders, G. D., & Heller, R. 2019, AJ, 158, 143 [Google Scholar]
  20. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  21. Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE, 9913, 99133E [NASA ADS] [Google Scholar]
  23. Kinemuchi, K., Barclay, T., Fanelli, M., et al. 2012, PASP, 124, 963 [Google Scholar]
  24. Kopal, Z. 1950, Harvard Coll. Obs. Circ., 454, 1 [Google Scholar]
  25. Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369 [Google Scholar]
  26. Lightkurve Collaboration (Cardoso, J. V. D. M., et al.) 2018, Astrophysics Source Code Library [record ascl:1812.013] [Google Scholar]
  27. Lomb, N. R. 1976, Ap&SS, 39, 447 [Google Scholar]
  28. McKinney, W. 2010, Python in Science Conference, Austin, Texas, 56 [CrossRef] [Google Scholar]
  29. Mosteller, F., & Tukey, J. W. 1977, Data Analysis and Regression: A Second Course in Statistics, Addison-Wesley Series in Behavioral Science (Reading, Mass: Addison-Wesley Pub. Co) [Google Scholar]
  30. Murgas, F., Castro-González, A., Pallé, E., et al. 2023, A&A, 677, A182 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Parviainen, H. 2015, MNRAS, 450, 3233 [Google Scholar]
  32. Parviainen, H., & Aigrain, S. 2015, MNRAS, 453, 3822 [Google Scholar]
  33. Parviainen, H., Luque, R., & Palle, E. 2024, MNRAS, 527, 5693 [Google Scholar]
  34. Rasmussen, C. E., & Williams, C. K. I. 2006, Gaussian Processes for Machine Learning, Adaptive Computation and Machine Learning (Cambridge, Mass: MIT Press) [Google Scholar]
  35. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1, 014003D [Google Scholar]
  36. Savitzky, A., & Golay, M. J. E. 1964, Anal. Chem., 36, 1627 [Google Scholar]
  37. Scargle, J. D. 1982, AJ, 263, 835 [Google Scholar]
  38. Schweitzer, A., Passegger, V. M., Cifuentes, C., et al. 2019, A&A, 625, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Skilling, J. 2004, AIP, 735, 395 [NASA ADS] [Google Scholar]
  40. Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 [Google Scholar]
  41. Speagle, J. S. 2020, MNRAS, 493, 3132 [Google Scholar]
  42. Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985 [Google Scholar]
  43. Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100 [Google Scholar]
  44. Suárez Mascareño, A., González-Álvarez, E., Zapatero Osorio, M. R., et al. 2023, A&A, 670, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. The Pandas Development Team 2022, https://doi.org/10.5281/zenodo.3509134 [Google Scholar]
  46. VanderPlas, J. T., & Ivezić, Ž. 2015, AJ, 812, 18 [NASA ADS] [Google Scholar]
  47. Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat Methods, 17, 261 [CrossRef] [PubMed] [Google Scholar]
  48. Waskom, M. 2021, JOSS, 6, 3021 [NASA ADS] [CrossRef] [Google Scholar]
  49. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]

Supplementary Material

thumbnail Fig. A.1.

Characteristics of detrended PDCSAP photometry. (a) Distribution plot, conventionally split in two into a lower distribution (blue bars) and an upper distribution (green bars). The size of these sub-distributions are given in figures of corresponding colour. A Gaussian-fuction fit to the distribution is given as a black solid line. (b) Probability plot assessing whether the full distribution (black solid line) follows the fitted Gaussian (black dashed line).

Table A.1.

Relevant parameters of Barnard b and planetary candidates from the four-planet solution in González Hernández et al. (2024).

thumbnail Fig. A.2.

rior (black, cross-hatched) and posterior (blue, dot-hatched) parameter distributions in our narrow-prior radius search of Barnard b.

thumbnail Fig. A.3.

Orbital-inclination and planetary-radius distributions in our wide-prior radius search of Barnard b. The black dots mark the location of all samples in our posterior distribution. The shaded contours highlight 20% 40%, 60%, and 80% of the enclosed probability mass.

All Tables

Table 1.

Relevant parameters of Barnard’s star.

Table 2.

Relevant parameters of Barnard b.

Table A.1.

Relevant parameters of Barnard b and planetary candidates from the four-planet solution in González Hernández et al. (2024).

All Figures

thumbnail Fig. 1.

TESS S80 photometric data products: (a) SAP; (b) PDCSAP, with superimposed predicted times of inferior conjunction through GH24 ephemerides. The propagated inferior-conjunction uncertainty is shown as light blue bands.

In the text
thumbnail Fig. 2.

The search of Barnard b in TESS S80 photometry. (a) Mean time-window detrended SAP photometry with superimposed predicted times of inferior conjunction through GH24 ephemerides. The median expected transit width is displayed as dark blue bands. The propagated ephemeris uncertainty is shown as light blue bands. (b) Phase-folded plot at GH24 ephemerides: before and after applying a 5 min binning (light hollow and dark filled circles, respectively). The inferior-conjunction epoch uncertainty is shown as a light blue band. We compare photometry against the expected transit curve of Barnard b for three impact-parameter values: b = 0.0 (black dotted line), b = 0.9 (black dashed line), and b = 1.0 (black solid line). (c,d) GLSP of unbinned data before and after iterative 4σ clipping. (e,f) BLSP of unbinned data before and after iterative 4σ clipping.

In the text
thumbnail Fig. A.1.

Characteristics of detrended PDCSAP photometry. (a) Distribution plot, conventionally split in two into a lower distribution (blue bars) and an upper distribution (green bars). The size of these sub-distributions are given in figures of corresponding colour. A Gaussian-fuction fit to the distribution is given as a black solid line. (b) Probability plot assessing whether the full distribution (black solid line) follows the fitted Gaussian (black dashed line).

In the text
thumbnail Fig. A.2.

rior (black, cross-hatched) and posterior (blue, dot-hatched) parameter distributions in our narrow-prior radius search of Barnard b.

In the text
thumbnail Fig. A.3.

Orbital-inclination and planetary-radius distributions in our wide-prior radius search of Barnard b. The black dots mark the location of all samples in our posterior distribution. The shaded contours highlight 20% 40%, 60%, and 80% of the enclosed probability mass.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.