Search

CHEOPS is a space telescope specifically designed to monitor transiting exoplanets orbiting bright stars. In September 2023, CHEOPS completed its nominal mission and remains in excellent operational conditions. The mission has been extended until the end of 2026. Scientific and instrumental data have been collected throughout in-orbit commissioning and nominal operations, enabling a comprehensive analysis of the mission's performance. In this article, we present the results of this analysis with a twofold goal. First, we aim to inform the scientific community about the present status of the mission and what can be expected as the instrument ages. Secondly, we intend for this publication to serve as a legacy document for future missions, providing insights and lessons learned from the successful operation of CHEOPS. To evaluate the instrument performance in flight, we developed a comprehensive monitoring and characterisation programme. It consists of dedicated observations that allow us to characterise the instrument's response. In addition to the standard collection of nominal science and housekeeping data, these observations provide input for detecting, modelling, and correcting instrument systematics, discovering and addressing anomalies, and comparing the instrument's actual performance with expectations. The precision of the CHEOPS measurements has enabled the mission objectives to be met and exceeded. Careful modelling of the instrumental systematics allows the data quality to be significantly improved during the light curve analysis phase, resulting in more precise scientific measurements. CHEOPS is compliant with the driving scientific requirements of the mission. Although visible, the ageing of the instrument has not affected the mission's performance., Comment: Accepted for publication in Astronomy and Astrophysics

Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and UV emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars' habitable zone. Aims. We used the high-cadence, high-photometric capabilities of the CHEOPS and TESS space telescopes to study the detailed morphology of white-light flares occurring in a sample of 130 late-K and M stars, and compared our findings with results obtained at a lower cadence. We developed dedicated software for this purpose. Results. Multi-peak flares represent a significant percentage ($\gtrsim 30$\%) of the detected outburst events. Our findings suggest that high-impulse flares are more frequent than suspected from lower-cadence data, so that the most impactful flux levels that hit close-in exoplanets might be more time-limited than expected. We found significant differences in the duration distributions of single-peak and complex flare components, but not in their peak luminosity. A statistical analysis of the flare parameter distributions provides marginal support for their description with a log-normal instead of a power-law function, leaving the door open to several flare formation scenarios. We tentatively confirmed previous results about quasi-periodic pulsations in high-cadence photometry, report the possible detection of a pre-flare dip, and did not find hints of photometric variability due to an undetected flare background. Conclusions. The high-cadence study of stellar hosts might be crucial to evaluate the impact of their flares on close-in exoplanets, as their impulsive phase emission might otherwise be incorrectly estimated. Future telescopes such as PLATO and Ariel will help in this respect., Comment: 28 pages, 25 figures, 4 tables, to be published in Astronomy & Astrophysics

We aim to constrain the internal structure and composition of HD 15337 b and c, two short-period planets situated on opposite sides of the radius valley, using new transit photometry and radial velocity data. We acquire 6 new transit visits with the CHaracterising ExOPlanet Satellite (CHEOPS) and 32 new radial velocity measurements from the High Accuracy Radial Velocity Planet Searcher (HARPS) to improve the accuracy of the mass and radius estimates for both planets. We reanalyse light curves from TESS sectors 3 and 4 and analyse new data from sector 30, correcting for long-term stellar activity. Subsequently, we perform a joint fit of the TESS and CHEOPS light curves, and all available RV data from HARPS and the Planet Finder Spectrograph (PFS). Our model fits the planetary signals, the stellar activity signal and the instrumental decorrelation model for the CHEOPS data simultaneously. The stellar activity was modelled using a Gaussian-process regression on both the RV and activity indicators. We finally employ a Bayesian retrieval code to determine the internal composition and structure of the planets. We derive updated and highly precise parameters for the HD 15337 system. Our improved precision on the planetary parameters makes HD 15337 b one of the most precisely characterised rocky exoplanets, with radius and mass measurements achieving a precision better than 2\% and 7\%, respectively. We are able to improve the precision of the radius measurement of HD 15337 c to 3\%. Our results imply that the composition of HD 15337 b is predominantly rocky, while HD 15337 c exhibits a gas envelope with a mass of at least $0.01\ M_\oplus$.Our results lay the groundwork for future studies, which can further unravel the atmospheric evolution of these exoplanets and give new insights into their composition and formation history and the causes behind the radius gap., Comment: 17 pages, including appendix

Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$\sigma$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared., Comment: accepted for publication in A&A

Context - Around the year 2000, Triton's south pole experienced an extreme summer solstice that occurs every about 650 years, when the subsolar latitude reached about 50{\deg}. Bracketing this epoch, a few occultations probed Triton's atmosphere in 1989, 1995, 1997, 2008 and 2017. A recent ground-based stellar occultation observed on 6 October 2022 provides a new measurement of Triton's atmospheric pressure which is presented here. Aims- The goal is to constrain the Volatile Transport Models (VTMs) of Triton's atmosphere that is basically in vapor pressure equilibrium with the nitrogen ice at its surface. Methods - Fits to the occultation light curves yield Triton's atmospheric pressure at the reference radius 1400 km, from which the surface pressure is induced. Results - The fits provide a pressure p_1400= 1.211 +/- 0.039 microbar at radius 1400 km (47 km altitude), from which a surface pressure of p_surf= 14.54 +/- 0.47 microbar is induced (1-sigma error bars). To within error bars, this is identical to the pressure derived from the previous occultation of 5 October 2017, p_1400 = 1.18 +/- 0.03 microbar and p_surf= 14.1 +/- 0.4 microbar, respectively. Based on recent models of Triton's volatile cycles, the overall evolution over the last 30 years of the surface pressure is consistent with N2 condensation taking place in the northern hemisphere. However, models typically predict a steady decrease in surface pressure for the period 2005-2060, which is not confirmed by this observation. Complex surface-atmosphere interactions, such as ice albedo runaway and formation of local N2 frosts in the equatorial regions of Triton could explain the relatively constant pressure between 2017 and 2022., Comment: 8 pages, 4 figures, accepted for publication in Astronomy and Astrophysics

HD 139139 (a.k.a. 'The Random Transiter') is a star that exhibited enigmatic transit-like features with no apparent periodicity in K2 data. The shallow depth of the events ($\sim$200 ppm -- equivalent to transiting objects with radii of $\sim$1.5 R$_\oplus$ in front of a Sun-like star), and their non-periodicity, constitutes a challenge for the photometric follow-up of this star. The goal of this study is to confirm with independent measurements the presence of shallow, non-periodic transit-like features on this object. We performed observations with CHEOPS, for a total accumulated time of 12.75 d, distributed in visits of roughly 20 h in two observing campaigns in years 2021 and 2022. The precision of the data is sufficient to detect 150 ppm features with durations longer than 1.5 h. We use the duration and times of the events seen in the K2 curve to estimate how many should have been detected in our campaigns, under the assumption that their behaviour during the CHEOPS observations would be the same as in the K2 data of 2017. We do not detect events with depths larger than 150 ppm in our data set. If the frequency, depth, and duration of the events were the same as in the K2 campaign, we estimate the probability of having missed all events due to our limited observing window would be 4.8 %. We suggest three different scenarios to explain our results: 1) Our observing window was not long enough, and the events were missed with the estimated 4.8 % probability. 2) The events recorded in the K2 observations were time critical, and the mechanism producing them was either not active in the 2021 and 2022 campaigns or created shallower events under our detectability level. 3) The enigmatic events in the K2 data are the result of an unidentified and infrequent instrumental noise in the original data set or its data treatment., Comment: 7 pages, 3 figures, accepted for publication in A&A. Language-corrected version

55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-curve modulation in 55 Cnc e. To this end, we used the CHaracterising ExOPlanet Satellite (CHEOPS), whose exquisite photometric precision provides an opportunity to characterise minute changes in the phase curve from one orbit to the next. CHEOPS observed 29 individual visits of 55 Cnc e between March 2020 and February 2022. Based on these observations, we investigated the different processes that could be at the origin of the observed modulation. In particular, we built a toy model to assess whether a circumstellar torus of dust driven by radiation pressure and gravity might match the observed flux variability timescale. We find that the phase-curve amplitude and peak offset of 55 Cnc e do vary between visits. The sublimation timescales of selected dust species reveal that silicates expected in an Earth-like mantle would not survive long enough to explain the observed phase-curve modulation. We find that silicon carbide, quartz, and graphite are plausible candidates for the circumstellar torus composition because their sublimation timescales are long. The extensive CHEOPS observations confirm that the phase-curve amplitude and offset vary in time.We find that dust could provide the grey opacity source required to match the observations. However, the data at hand do not provide evidence that circumstellar material with a variable grain mass per unit area causes the observed variability. Future observations with the James Webb Space Telescope promise exciting insights into this iconic super-Earth., Comment: 27 pages, 22 figures. Accepted for publication on A&A

TOI-1055 is a Sun-like star known to host a transiting Neptune-sized planet on a 17.5-day orbit (TOI-1055 b). Radial velocity (RV) analyses carried out by two independent groups using nearly the same set of HARPS spectra have provided measurements of planetary masses that differ by $\sim$ 2$\sigma$. Our aim in this work is to solve the inconsistency in the published planetary masses by significantly extending the set of HARPS RV measurements and employing a new analysis tool that is able to account and correct for stellar activity. Our further aim was to improve the precision on measurements of the planetary radius by observing two transits of the planet with the CHEOPS space telescope. We fit a skew normal (SN) function to each cross correlation function extracted from the HARPS spectra to obtain RV measurements and hyperparameters to be used for the detrending. We evaluated the correlation changes of the hyperparameters along the RV time series using the breakpoint technique. We performed a joint photometric and RV analysis using a Markov chain Monte Carlo (MCMC) scheme to simultaneously detrend the light curves and the RV time series. We firmly detected the Keplerian signal of TOI-1055 b, deriving a planetary mass of $M_b=20.4_{-2.5}^{+2.6} M_{\oplus}$ ($\sim$12%). This value is in agreement with one of the two estimates in the literature, but it is significantly more precise. Thanks to the TESS transit light curves combined with exquisite CHEOPS photometry, we also derived a planetary radius of $R_b=3.490_{-0.064}^{+0.070} R_{\oplus}$ ($\sim$1.9%). Our mass and radius measurements imply a mean density of $\rho_b=2.65_{-0.35}^{+0.37}$ g cm$^{-3}$ ($\sim$14%). We further inferred the planetary structure and found that TOI-1055 b is very likely to host a substantial gas envelope with a mass of $0.41^{+0.34}_{-0.20}$ M$_\oplus$ and a thickness of $1.05^{+0.30}_{-0.29}$ R$_\oplus$., Comment: 13 pages, 6 figures, 5 tables. Accepted for publication in A&A

The planetary system around the naked-eye star $\nu^2$ Lupi (HD 136352; TOI-2011) is composed of three exoplanets with masses of 4.7, 11.2, and 8.6 Earth masses. The TESS and CHEOPS missions revealed that all three planets are transiting and have radii straddling the radius gap separating volatile-rich and volatile-poor super-earths. Only a partial transit of planet d had been covered so we re-observed an inferior conjunction of the long-period 8.6 Earth-mass exoplanet $\nu^2$ Lup d with the CHEOPS space telescope. We confirmed its transiting nature by covering its whole 9.1 h transit for the first time. We refined the planet transit ephemeris to P = 107.1361 (+0.0019/-0.0022) days and Tc = 2,459,009.7759 (+0.0101/-0.0096) BJD_TDB, improving by ~40 times on the previously reported transit timing uncertainty. This refined ephemeris will enable further follow-up of this outstanding long-period transiting planet to search for atmospheric signatures or explore the planet's Hill sphere in search for an exomoon. In fact, the CHEOPS observations also cover the transit of a large fraction of the planet's Hill sphere, which is as large as the Earth's, opening the tantalising possibility of catching transiting exomoons. We conducted a search for exomoon signals in this single-epoch light curve but found no conclusive photometric signature of additional transiting bodies larger than Mars. Yet, only a sustained follow-up of $\nu^2$ Lup d transits will warrant a comprehensive search for a moon around this outstanding exoplanet., Comment: 16 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

We report the discovery of a third planet transiting the star TOI-1260, previously known to host two transiting sub-Neptune planets with orbital periods of 3.127 and 7.493 days, respectively. The nature of the third transiting planet with a 16.6-day orbit is supported by ground-based follow-up observations, including time-series photometry, high-angular resolution images, spectroscopy, and archival imagery. Precise photometric monitoring with CHEOPS allows to improve the constraints on the parameters of the system, improving our knowledge on their composition. The improved radii of TOI-1260b, TOI-1260c are $2.36 \pm 0.06 \rm R_{\oplus}$, $2.82 \pm 0.08 \rm R_{\oplus}$, respectively while the newly discovered third planet has a radius of $3.09 \pm 0.09 \rm R_{\oplus}$. The radius uncertainties are in the range of 3\%, allowing a precise interpretation of the interior structure of the three planets. Our planet interior composition model suggests that all three planets in the TOI-1260 system contains some fraction of gas. The innermost planet TOI-1260b has most likely lost all of its primordial hydrogen-dominated envelope. Planets c and d were also likely to have experienced significant loss of atmospheric through escape, but to a lesser extent compared to planet b., Comment: 18 pages, 10 figures, accepted for publication in MNRAS

Tidal orbital decay is suspected to occur especially for hot Jupiters, with the only observationally confirmed case of this being WASP-12 b. By examining this effect, information on the properties of the host star can be obtained using the so-called stellar modified tidal quality factor $Q_*'$, which describes the efficiency with which kinetic energy of the planet is dissipated within the star. This can help to get information about the interior of the star. In this study, we aim to improve constraints on the tidal decay of the KELT-9, KELT-16 and WASP-4 systems, to find evidence for or against the presence of this particular effect. With this, we want to constrain each star's respective $Q_*'$ value. In addition to that, we also aim to test the existence of the transit timing variations (TTVs) in the HD 97658 system, which previously favoured a quadratic trend with increasing orbital period. Making use of newly acquired photometric observations from CHEOPS and TESS, combined with archival transit and occultation data, we use Markov chain Monte Carlo (MCMC) algorithms to fit three models, a constant period model, an orbital decay model, and an apsidal precession model, to the data. We find that the KELT-9 system is best described by an apsidal precession model for now, with an orbital decay trend at over 2 $\sigma$ being a possible solution as well. A Keplerian orbit model with a constant orbital period fits the transit timings of KELT-16 b the best due to the scatter and scale of their error bars. The WASP-4 system is represented the best by an orbital decay model at a 5 $\sigma$ significance, although apsidal precession cannot be ruled out with the present data. For HD 97658 b, using recently acquired transit observations, we find no conclusive evidence for a previously suspected strong quadratic trend in the data., Comment: 19 pages, 11 figures, 9 tables, accepted for publication in A&A

The HD108236 system was first announced with the detection of four small planets based on TESS data. Shortly after, the transit of an additional planet with a period of 29.54d was serendipitously detected by CHEOPS. In this way, HD108236 (V=9.2) became one of the brightest stars known to host five small transiting planets (R$_p$<3R$_{\oplus}$). We characterize the planetary system by using all the data available from CHEOPS and TESS space missions. We use the flexible pointing capabilities of CHEOPS to follow up the transits of all the planets in the system, including the fifth transiting body. After updating the host star parameters by using the results from Gaia eDR3, we analyzed 16 and 43 transits observed by CHEOPS and TESS, respectively, to derive the planets physical and orbital parameters. We carried out a timing analysis of the transits of each of the planets of HD108236 to search for the presence of transit timing variations. We derived improved values for the radius and mass of the host star (R$_{\star}$=0.876$\pm$0.007 R$_{\odot}$ and M$_{\star}$=0.867$_{-0.046}^{+0.047}$ M$_{\odot}$). We confirm the presence of the fifth transiting planet f in a 29.54d orbit. Thus, the system consists of five planets of R$_b$=1.587$\pm$0.028, R$_c$=2.122$\pm$0.025, R$_d$=2.629$\pm$0.031, R$_e$=3.008$\pm$0.032, and R$_f$=1.89$\pm$0.04 [R$_{\oplus}$]. We refine the transit ephemeris for each planet and find no significant transit timing variations for planets c, d, and e. For planets b and f, instead, we measure significant deviations on their transit times (up to 22 and 28 min, respectively) with a non-negligible dispersion of 9.6 and 12.6 min in their time residuals. We confirm the presence of planet f and find no significant evidence for a potential transiting planet in a 10.9d orbital period, as previously suggested. Full abstract in the PDF file., Comment: 18 Figures and 25 pages. Accepted for publication in A&A

Much remains to be understood about the nature of exoplanets smaller than Neptune, most of which have been discovered in compact multi-planet systems. With its inner ultra-short period planet b aligned with the star and two larger outer planets d-c on polar orbits, the multi-planet system HD 3167 features a peculiar architecture and offers the possibility to investigate both dynamical and atmospheric evolution processes. To this purpose we combined multiple datasets of transit photometry and radial velocimetry (RV) to revise the properties of the system and inform models of its planets. This effort was spearheaded by CHEOPS observations of HD 3167b, which appear inconsistent with a purely rocky composition despite its extreme irradiation. Overall the precision on the planetary orbital periods are improved by an order of magnitude, and the uncertainties on the densities of the transiting planets b and c are decreased by a factor of 3. Internal structure and atmospheric simulations draw a contrasting picture between HD 3167d, likely a rocky super-Earth that lost its atmosphere through photo-evaporation, and HD 3167c, a mini-Neptune that kept a substantial primordial gaseous envelope. We detect a fourth, more massive planet on a larger orbit, likely coplanar with HD 3167d-c. Dynamical simulations indeed show that the outer planetary system d-c-e was tilted, as a whole, early in the system history, when HD 3167b was still dominated by the star influence and maintained its aligned orbit. RV data and direct imaging rule out that the companion that could be responsible for the present-day architecture is still bound to the HD\,3167 system. Similar global studies of multi-planet systems will tell how many share the peculiar properties of the HD3167 system, which remains a target of choice for follow-up observations and simulations., Comment: 22 pages, 23 pages, accepted for publication in A&A (18 August 2022). Updated author list in new version

Stellar occultation is a powerful technique that allows the determination of some physical parameters of the occulting object. The result depends on the photometric accuracy, the temporal resolution, and the number of chords obtained. Space telescopes can achieve high photometric accuracy as they are not affected by atmospheric scintillation. Using ESA's CHEOPS space telescope, we observed a stellar occultation by the Transneptunian object (50000) Quaoar. We compare the obtained chord with previous occultations by this object and determine its astrometry with sub-milliarcsecond precision. Also, we determine upper limits to the presence of a global methane atmosphere on the occulting body. We predicted and observed a stellar occultation by Quaoar using the CHEOPS space telescope. We measured the occultation light curve from this data-set and determined the dis- and re-appearance of the star behind the occulting body. Furthermore, a ground-based telescope in Australia was used to constrain Quaoar's limb. Combined with results from previous works, these measurements allow us to obtain a precise position of Quaoar at the occultation time. We present results obtained from the first stellar occultation by a Transneptunian object (TNO) using space telescope orbiting Earth. It was the occultation by Quaoar observed on 2020 June 11. We used the CHEOPS light curve to obtain a surface pressure upper limit of 85 nbar for the detection of a global methane atmosphere. Also, combining this observation with a ground-based observation we fit Quaoar's limb to determine its astrometric position with an uncertainty below 1.0 mas. This observation is a first of its kind, and it shall be considered as a proof of concept of stellar occultation observations of Transneptunian objects with space telescopes orbiting Earth. Moreover, it shows significant prospects for the James Webb Space Telescope., Comment: 9 pages, 10 figures

We present the discovery of two small planets transiting HD 93963A (TOI-1797), a G0\,V star (M$_*$=1.109\,$\pm$\,0.043\,M$_\odot$, R$_*$=1.043\,$\pm$\,0.009\,R$_\odot$) in a visual binary system. We combined TESS and CHEOPS space-borne photometry with data from MuSCAT 2, `Alopeke, PHARO, TRES, FIES, and SOPHIE. We validated and spectroscopically confirmed the outer transiting planet HD 93963 Ac, a sub-Neptune with an orbital period of P$_c \approx$ 3.65 d, reported as a TESS object of interest (TOI) shortly after the release of Sector 22 data. HD 93963 Ac has a mass of M$_c = 19.2 \pm 4.1$ M$_{\oplus}$ and a radius of R$_c = 3.228 \pm 0.059$ R$_{\oplus}$, implying a mean density of $\rho_c=3.1\pm0.7$ gcm$^{-3}$. The inner object, HD 93963 Ab, is a validated 1.04 d ultra-short period (USP) transiting super-Earth that we discovered in the TESS light curve and that was not listed as a TOI, owing to the low significance of its signal (TESS signal-to-noise ratio $\approx$ 6.7, TESS $+$ CHEOPS combined transit depth D$_b=141.5 \pm 8.5$ ppm). We intensively monitored the star with CHEOPS by performing nine transit observations to confirm the presence of the inner planet and validate the system. HD 93963 Ab is the first small (R$_b = 1.35 \pm 0.042$ R$_{\oplus}$) USP planet discovered and validated by TESS and CHEOPS. Unlike planet c, HD 93963 Ab is not significantly detected in our radial velocities (M$_b = 7.8 \pm 3.2$ M$_{\oplus}$). We also discovered a linear trend in our Doppler measurements, suggesting the possible presence of a long-period outer planet. With a V-band magnitude of 9.2, HD 93963 A is among the brightest stars known to host a USP planet, making it one of the most favourable targets for precise mass measurement via Doppler spectroscopy and an important laboratory to test formation, evolution, and migration models of planetary systems hosting ultra-short period planets., Comment: Accepted in Astronomy & Astrophysics

We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($\rho_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies., Comment: 22 pages, 19 figures. Accepted for publication in MNRAS

CHEOPS(CHaracterising ExOPlanet Satellite) is an ESA S-class mission that observes bright stars at high cadence from low-Earth orbit. The main aim of the mission is to characterize exoplanets that transit nearby stars using ultrahigh precision photometry. Here we report the analysis of transits observed by CHEOPS during its Early Science observing programme for four well-known exoplanets: GJ436b, HD106315b, HD97658b and GJ1132b. The analysis is done using pycheops, an open-source software package we have developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques that are fully described herein. We show that the precision of the transit parameters measured using CHEOPS is comparable to that from larger space telescopes such as Spitzer Space Telescope and Kepler. We use the updated planet parameters from our analysis to derive new constraints on the internal structure of these four exoplanets., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review

55 Cnc e is a transiting super-Earth (radius $1.88\rm\,R_\oplus$ and mass $8\rm\, M_\oplus$) orbiting a G8V host star on a 17-hour orbit. Spitzer observations of the planet's phase curve at 4.5 $\mu$m revealed a time-varying occultation depth, and MOST optical observations are consistent with a time-varying phase curve amplitude and phase offset of maximum light. Both broadband and high-resolution spectroscopic analyses are consistent with either a high mean molecular weight atmosphere or no atmosphere for planet e. A long term photometric monitoring campaign on an independent optical telescope is needed to probe the variability in this system. We seek to measure the phase variations of 55 Cnc e with a broadband optical filter with the 30 cm effective aperture space telescope CHEOPS and explore how the precision photometry narrows down the range of possible scenarios. We observed 55 Cnc for 1.6 orbital phases in March of 2020. We designed a phase curve detrending toolkit for CHEOPS photometry which allows us to study the underlying flux variations of the 55 Cnc system. We detected a phase variation with a full-amplitude of $72 \pm 7$ ppm but do not detect a significant secondary eclipse of the planet. The shape of the phase variation resembles that of a piecewise-Lambertian, however the non-detection of the planetary secondary eclipse, and the large amplitude of the variations exclude reflection from the planetary surface as a possible origin of the observed phase variations. They are also likely incompatible with magnetospheric interactions between the star and planet but may imply that circumplanetary or circumstellar material modulate the flux of the system. Further precision photometry of 55 Cnc from CHEOPS will measure variations in the phase curve amplitude and shape over time this year., Comment: 15 pages, accepted by A&A

The accuracy of theoretical mass, radius and effective temperature values for M-dwarf stars is an active topic of debate. Differences between observed and theoretical values have raised the possibility that current theoretical stellar structure and evolution models are inaccurate towards the low-mass end of the main sequence. To explore this issue we use the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions. We use these light curves combined with the spectroscopic orbit for the solar-type companion to measure the mass, radius and effective temperature of the M-dwarf star. Here we present the analysis of three eclipsing binaries. We use the pycheops data analysis software to fit the observed transit and eclipse events of each system. Two of our systems were also observed by the TESS satellite -- we similarly analyse these light curves for comparison. We find consistent results between CHEOPS and TESS, presenting three stellar radii and two stellar effective temperature values of low-mass stellar objects. These initial results from our on-going observing programme with CHEOPS show that we can expect to have ~24 new mass, radius and effective temperature measurements for very low mass stars within the next few years., Comment: 12 pages, 8 figures, accepted for publication in MNRAS

Black-hole binary (BHB) systems comprise a stellar-mass black hole and a closely orbiting companion star. Matter is transferred from the companion to the black hole, forming an accretion disk, corona and jet structures. The resulting release of gravitational energy leads to emission of X-rays. The radiation is affected by special/general relativistic effects, and can serve as a probe of the properties of the black hole and surrounding environment, if the accretion geometry is properly identified. Two competing models describe the disk-corona geometry for the hard spectral state of BHBs, based on spectral and timing measurements. Measuring the polarization of hard X-rays reflected from the disk allows the geometry to be determined. The extent of the corona differs between the two models, affecting the strength of relativistic effects (e.g., enhancement of polarization fraction and rotation of polarization angle). Here, we report observational results on linear polarization of hard X-ray (19-181 keV) emission from a BHB, Cygnus X-1, in the hard state. The low polarization fraction, <8.6% (upper limit at 90% confidence level), and the alignment of the polarization angle with the jet axis show that the dominant emission is not influenced by strong gravity. When considered together with existing spectral and timing data, our result reveals that the accretion corona is either an extended structure, or is located far from the black hole in the hard state of Cygnus X-1., Comment: 44 pages, 15 figures, 1 table

The linear polarisation fraction and angle of the hard X-ray emission from the Crab provide unique insight into high energy radiation mechanisms, complementing the usual imaging, timing and spectroscopic approaches. Results have recently been presented by two missions operating in partially overlapping energy bands, PoGO+ (18-160 keV) and AstroSat CZTI (100-380 keV). We previously reported PoGO+ results on the polarisation parameters integrated across the light-curve and for the entire nebula-dominated off-pulse region. We now introduce finer phase binning, in light of the AstroSat CZTI claim that the polarisation fraction varies across the off-pulse region. Since both missions are operating in a regime where errors on the reconstructed polarisation parameters are non-Gaussian, we adopt a Bayesian approach to compare results from each mission. We find no statistically significant variation in off-pulse polarisation parameters, neither when considering the mission data separately nor when they are combined. This supports expectations from standard high-energy emission models., Comment: 6 pages, 3 figures, 5 tables. Accepted for publication in MNRAS Letters

Strong magnetic fields, synchrotron emission, and Compton scattering are omnipresent in compact celestial X-ray sources. Emissions in the X-ray energy band are consequently expected to be linearly polarized. X-ray polarimetry provides a unique diagnostic to study the location and fundamental mechanisms behind emission processes. The polarization of emissions from a bright celestial X-ray source, the Crab, is reported here for the first time in the hard X-ray band (~20-160 keV). The Crab is a complex system consisting of a central pulsar, a diffuse pulsar wind nebula, as well as structures in the inner nebula including a jet and torus. Measurements are made by a purpose-built and calibrated polarimeter, PoGO+. The polarization vector is found to be aligned with the spin axis of the pulsar for a polarization fraction, PF = (20.9 $\pm$ 5.0)%. This is higher than that of the optical diffuse nebula, implying a more compact emission site, though not as compact as, e.g., the synchrotron knot. Contrary to measurements at higher energies, no significant temporal evolution of phase-integrated polarisation parameters is observed. The polarization parameters for the pulsar itself are measured for the first time in the X-ray energy band and are consistent with observations at optical wavelengths., Comment: Accepted for publication in Scientific Reports, www.nature.com/srep

We have measured the linear polarisation of hard X-ray emission from the Crab in a previously unexplored energy interval, 20-120 keV. The introduction of two new observational parameters, the polarisation fraction and angle stands to disentangle geometrical and physical effects, thereby providing information on the pulsar wind geometry and magnetic field environment. Measurements are conducted using the PoGOLite Pathfinder - a balloon-borne polarimeter. Polarisation is determined by measuring the azimuthal Compton scattering angle of incident X-rays in an array of plastic scintillators housed in an anticoincidence well. The polarimetric response has been characterised prior to flight using both polarised and unpolarised calibration sources. We address possible systematic effects through observations of a background field. The measured polarisation fraction for the integrated Crab light-curve is ($18.4^{+9.8}_{-10.6}$)%, corresponding to an upper limit (99% credibility) of 42.4%, for a polarisation angle of ($149.2\pm16.0)^\circ$., Comment: 6 pages, 5 figures. Accepted for publication in MNRAS Letters (In this revision the polarisation angle is correctly stated as 149.2$\pm$16.0 degrees and not 142.2$\pm$16.0 degrees (typo corrected in abstract and on page 4)

In the 50 years since the advent of X-ray astronomy there have been many scientific advances due to the development of new experimental techniques for detecting and characterising X-rays. Observations of X-ray polarisation have, however, not undergone a similar development. This is a shortcoming since a plethora of open questions related to the nature of X-ray sources could be resolved through measurements of the linear polarisation of emitted X-rays. The PoGOLite Pathfinder is a balloon-borne hard X-ray polarimeter operating in the 25 - 240 keV energy band from a stabilised observation platform. Polarisation is determined using coincident energy deposits in a segmented array of plastic scintillators surrounded by a BGO anticoincidence system and a polyethylene neutron shield. The PoGOLite Pathfinder was launched from the SSC Esrange Space Centre in July 2013. A near-circumpolar flight was achieved with a duration of approximately two weeks. The flight performance of the Pathfinder design is discussed for the three Crab observations conducted. The signal-to-background ratio for the observations is shown to be 0.25$\pm$0.03 and the Minimum Detectable Polarisation (99% C.L.) is (28.4$\pm$2.2)%. A strategy for the continuation of the PoGOLite programme is outlined based on experience gained during the 2013 maiden flight., Comment: Accepted for publication in Experimental Astronomy. 26 pages, 20 figures

Context. The binary B5V star BD+31?643 exhibits a disk-like structure detected at optical wavelengths. Even though the feature is well centered on the star, it has been argued, based on Spitzer observations, that the feature is a filament not directly associated to the binary star. Aims. The purpose of the present paper is to investigate whether polarization imaging may provide evidence either for or against the disk hypothesis. In addition, we aim at clarifying whether there might be any additional close companion to the binary star. Methods. We used the coronagraph PolCor in its polarization mode in combination with an EMCCD camera allowing short unit exposure times. As a result of shift-and-add and frame selection, the spatial resolution is improved compared to traditional CCD imaging. In order to possibly reveal an additional stellar companion, we used high resolution spectroscopy in the optical and high spatial resolution imaging in the near-IR. Results. The disk/filament is much better seen in polarization; it is narrow and a line drawn along the ridge passes within a second of arc from the star. The degree of polarization is high (?50% after correction for the extended component of the reflection nebula) which means that the disk/filament must be approximately at the same distance as the star. Although we confirm that the feature is much brighter south-east than north-west of the star, the evidence that the feature is physically connected to the star is strengthened and suggests that we are witnessing the destruction process of an accretion disk. Our spectroscopy shows that at least one of the stars is a spectroscopic binary. We were, however, not able to spatially resolve any stellar component in addition to the two well separated stars., Comment: 9 pages, 19 figures

The Odin satellite, launched in Feb. 2001, is equipped with a 1.1-m submillimetre telescope. Odin was used to observe the 557 GHz line of water with high spectral resolution in 12 comets between 2001 and 2005. Line shapes and spatial mapping provide information on the anisotropy of the outgassing and constraints on water excitation, enabling accurate measurements of the water production rate. Five comets were regularly observed over periods of more than one month to monitor the variation of their water outgassing rate with heliocentric distance. Observing campaigns have been generally coordinated with ground-based observations of molecular lines at Nan\c{c}ay, CSO or IRAM 30-m telescopes to obtain molecular abundances relative to water. Thanks to Odin's frequency coverage, it was also possible to detect the H_2^18O 548GHz line, first in comet 153P/Ikeya-Zhang in April 2002 (Lecacheux et al., 2003) and then in comets C/2002 T7 (LINEAR), C/2001 Q4 (NEAT) and C/2004 Q2 (Machholz). The ^16O/^18O isotopic ratio (\approx450) is consistent with the terrestrial value. Ammonia has been searched for in three comets through its J_K = 1_0-0_0 line at 572 GHz and was tentatively detected in C/2001 Q4 and C/2002 T7. The derived abundances of NH_3 relative to water are 0.5 % and 0.3 %, respectively, similar to values obtained in other comets with different techniques., Comment: accepted for publication in Planetary and Space Science

Comet 9P/Tempel 1 was the target of a multi-wavelength worldwide investigation in 2005. The NASA Deep Impact mission reached the comet on 4.24 July 2005, delivering a 370 kg impactor which hit the comet at 10.3 km/s. Following this impact, a cloud of gas and dust was excavated from the comet nucleus. The comet was observed in 2005 prior to and after the impact, at 18-cm wavelength with the Nan\c{c}ay radio telescope, in the millimetre range with the IRAM and CSO radio telescopes, and at 557 GHz with the Odin satellite. Post-impact observations at IRAM and CSO did not reveal a significant change of the outgassing rates and relative abundances, with the exception of CH_3OH which may have been more abundant by up to one order of magnitude in the ejecta. Most other variations are linked to the intrinsic variability of the comet. The Odin satellite monitored nearly continuously the H_2O line at 557 GHz during the 38 hours following the impact on the 4th of July, in addition to weekly monitoring. Once the periodic variations related to the nucleus rotation are removed, a small increase of outgassing related to the impact is present, which corresponds to the release of \approx 5000+/-2000 tons of water. Two other bursts of activity, also observed at other wavelengths, were seen on 23 June and 7 July; they correspond to even larger releases of gas., Comment: accepted for publication in Icarus

Aims. The aim of this paper is to characterise the star formation activity in the poorly studied embedded cluster Serpens/G3-G6, located ~ 45' (3 pc) to the south of the Serpens Cloud Core, and to determine the luminosity and mass functions of its population of Young Stellar Objects (YSOs). Methods. Multi-wavelength broadband photometry was obtained to sample the near and mid-IR spectral energy distributions to separate YSOs from field stars and classify the YSO evolutionary stage. ISOCAM mapping in the two filters LW2 (5-8.5 um) and LW3 (12-18 um) of a 19' x 16' field was combined with JHKs data from 2MASS, Ks data from Arnica/NOT, and L' data from SIRCA/NOT. Continuum emission at 1.3 mm (IRAM) and 3.6 cm (VLA) was mapped to study the cloud structure and the coldest/youngest sources. Deep narrow band imaging at the 2.12 um S(1) line of H2 from NOTCam/NOT was obtained to search for signs of bipolar outflows. Results. We have strong evidence for a stellar population of 31 Class II sources, 5 flat-spectrum sources, 5 Class I sources, and two Class 0 sources. Our method does not sample the Class III sources. The cloud is composed of two main dense clumps aligned along a ridge over ~ 0.5 pc plus a starless core coinciding with absorption features seen in the ISOCAM maps. We find two S-shaped bipolar collimated flows embedded in the NE clump, and propose the two driving sources to be a Class 0 candidate (MMS3) and a double Class I (MMS2). For the Class II population we find a best age of ~ 2 Myr and compatibility with recent Initial Mass Functions (IMFs) by comparing the observed Class II luminosity function (LF), which is complete to 0.08 L_sun, to various model LFs with different star formation scenarios and input IMFs., Comment: 18 pages, 16 figures, 3 online tables, accepted by A&A