Aims. Using HARPS spectroscopic data obtained by the RedDots campaign, as well as archival data from HARPS and CARMENES, supplemented with ASH2 and T90 photometry, we aim to search for additional planets around the three M dwarfs GJ 832, GJ 674, and Ross 128. We also aim to determine limits on possible undetected, habitable planets. We investigate (i) the reliability of the recovered orbital eccentricities and (ii) the reliability of Bayesian evidence as a diagnostic for selecting the best model. Methods. We employed Markov-chain Monte Carlo, nested sampling, and Gaussian process (GP) analyses to fit a total of 20 different models. We used the residuals to create grids for injection-recovery simulations to obtain detection limits on potentially undiscovered planets. Results. Our refined orbital elements for GJ 832 b, GJ 674 b, and Ross 128 b confirm (GJ 832, GJ 674) or increase (Ross 128) prior eccentricity determinations. No additional planets were found in any of the systems. The detection limits obtained for all three systems are between 30 and 50 cm/s for orbital periods in the range of 1 to 10 000 days. Using N-body simulations, we find that undiscovered secondary planets are unlikely (Ross 128) or incapable (GJ 674) of having caused the observed eccentricities of the known planets. We find that the eccentricity of GJ 832 b is not significantly different from zero. Conclusions. GJ 832 b, GJ 674 b, and Ross 128 b retain their status as hosting lonely and (for the latter two) eccentric planets. Finally, our results show that Bayesian evidence, when used in conjunction with GP, is not a robust diagnostic for selecting the best model in cases of low-activity stars. In such cases, we advise an inspection of the shapes of the posterior distributions and to ensure that relevant simulations are performed to assess the validity of the perceived best model.

The radial velocity (RV) method of exoplanet detection requires mitigation of nuisance signals arising from stellar activity. Using analytic cool and facular spot models, we explore the use of central line moments (CLMs) for recovering and monitoring rotation induced RV variability. Different spot distribution patterns, photosphere-spot contrast ratios, and the presence or absence of the convective blueshift lead to differences in CLM signals between M and G dwarfs. Harmonics of the rotation period are often recovered with the highest power in standard periodogram analyses. By contrast, we show the true stellar rotation may be more reliably recovered with string length minimization. For solar minimum activity levels, recovery of the stellar rotation signal from CLMs is found to require unfeasibly high signal-to-noise observations. The stellar rotation period can be recovered at solar maximum activity levels from CLMs for reasonable cross-correlation function (CCF) signal-to-noise ratios >1000-5000. The CLMs can be used to recover and monitor stellar activity through their mutual correlations and correlations with RV and bisector inverse span. The skewness of a CCF, a measure of asymmetry, is described by the third CLM, M-3. Our noise-free simulations indicate the linear RV versus M-3 correlation is up to 10 per cent higher than the RV versus bisector inverse span correlation. We find a corresponding similar to 5 per cent increase in linear correlation for CARMENES observations of the M star, AU Mic. We also assess the effectiveness of the time derivative of the second CLM, M-2, for monitoring stellar activity.

We present a series of nine Doppler images of the magnetically active K component of the RS CVn-type binary SZ Piscium (SZ Psc), based on high-resolution spectroscopic data collected from 2014 to 2018. We apply least-squares deconvolution (LSD) to all spectra to extract the average profiles with high signal-to-noise ratios for Doppler imaging. The surface maps of the K subgiant show starspots widely distributed along its latitude and longitude. A prominent, nonaxisymmetric polar spot around phase 0 is revealed by all images with sufficient phase coverage, which may be a stable feature on the K component. The starspots evolve over a timescale of 1 month. We have determined the surface shear rate of the K component from starspot maps reconstructed 10 days apart in 2017 November-December through the cross-correlation method. The surface differential rotation parameters are Omega eq = 1.591 +/- 0.002 rad day-1 and Delta Omega = 0.035 +/- 0.003 rad day-1. The absorption lines contributed from the tertiary component of SZ Psc are detected in all LSD profiles, and we measure the radial velocity of the binary system and the tertiary component to derive an elliptical orbit with a period of 1530 +/- 3 days and a mass of 0.75 +/- 0.06 M circle dot for the tertiary component.

Context. The nearby (d < 5 pc) M dwarfs GJ 832, GJ 674, and Ross 128 each host a single exoplanet, with Ross 128 b located within the optimistic habitable zone. Due to their low mass and close proximity, these three systems are prime candidates for further characterization studies. Aims. Using HARPS spectroscopic data obtained by the RedDots campaign, as well as archival data from HARPS and CARMENES, supplemented with ASH2 and T90 photometry, we aim to search for additional planets in the three systems. We also aim to determine limits on possible undetected, habitable planets. We investigate (i) the reliability of the recovered orbital eccentricities and (ii) the reliability of Bayesian evidence as a diagnostic for selecting the best model. Methods. We employed Markov-chain Monte Carlo, nested sampling, and Gaussian process (GP) analyses to fit a total of 20 different models comprising 0-2 Keplerian signals and three different GP kernels for stellar activity. We used the residuals to create grids for injection-recovery simulations to obtain detection limits on potentially undiscovered planets. Results. Our refined orbital elements for GJ 832 b, GJ 674 b, and Ross 128 b confirm (GJ 832, GJ 674) or increase (Ross 128) prior eccentricity determinations. No additional planets were found in any of the systems. The detection limits obtained for all three systems are between 30 and 50 cm s(-1) for orbital periods in the range of 1-10 000 days. This corresponds to habitable planet masses of <1.5M(circle plus)for GJ 832 and < 1M(circle plus) for GJ 674 and Ross 128. Using N-body simulations, we find that undiscovered secondary planets are unlikely (Ross 128) or incapable (GJ 674) of having caused the observed eccentricities of the known planets. We find that the eccentricity of GJ 832b is not significantly different from zero. Conclusions. GJ 832 b, GJ 674 b, and Ross 128 b retain their status as hosting lonely and (for the latter two) eccentric planets (e = 0.04, 0.24, 0.21; respectively). This is unexpected in classical planet formation scenarios, which favor circular orbits and multiplanet configurations, demonstrating that planet formation in these cases is more complicated than traditionally thought. Additionally, the eccentricity of Ross 128 indicates that it spends some of its orbit outside of the optimistic habitable zone. Finally, our results show that Bayesian evidence, when used in conjunction with GP, is not a robust diagnostic for selecting the best model in cases of low- activity stars. In such cases, we advise an inspection of the shapes of the posterior distributions and to ensure that relevant simulations are performed to assess the validity of the perceived best model.

Discovering exoplanets orbiting young Suns can provide insight into the formation and early evolution of our own solar system, but the extreme magnetic activity of young stars obfuscates exoplanet detection. Here we monitor the long-term magnetic field and chromospheric activity variability of the young solar analogue V889 Her, model the activity-induced radial velocity variations and evaluate the impacts of extreme magnetism on exoplanet detection thresholds. We map the magnetic field and surface brightness for 14 epochs between 2004 and 2019. Our results show potential 3-4 yr variations of the magnetic field which evolves from weak and simple during chromospheric activity minima to strong and complex during activity maxima but without any polarity reversals. A persistent, temporally-varying polar spot coexists with weaker, short-lived lower-latitude spots. Due to their different decay time-scales, significant differential rotation and the limited temporal coverage of our legacy data, we were unable to reliably model the activity-induced radial velocity using Gaussian Process regression. Doppler Imaging can be a useful method for modelling the magnetic activity jitter of extremely active stars using data with large phase gaps. Given our data and using Doppler Imaging to filter activity jitter, we estimate that we could detect Jupiter-mass planets with orbital periods of ∼3 d. A longer baseline of continuous observations is the best observing strategy for the detection of exoplanets orbiting highly active stars.

ABSTRACT We have created an up-to-date catalogue of 214 brown dwarfs (BDs) in binaries with P < 104 d. This allows us to examine the population statistics of the BD desert. We searched Gaia third data release (DR3) non-single star (NSS) results for orbital inclinations of BD candidates, deriving 12 new masses. Three remain as desert BDs whereas nine candidates are found to be low-mass stars. We improved the radial velocity (RV) solutions for three previously studied BD candidates. A further 19 BD masses with periods less than ∼1200 d were identified in the DR3 binary_masses data base. We confirm a valley in the mass distribution with a minimum around 30–35 Mjup, and find that periods <100 d are still underpopulated in comparison with longer periods. The updated mass and eccentricity distribution of BDs still shows a marginally statistically significant split into high- and low-mass BD populations. This hints at two different parent distributions, and two potential origins – either akin to planetary formation, or stellar. There are no low-eccentricity BDs at periods around 100 d. The mass–metallicity distribution of BDs indicates that core accretion is not the dominant formation mechanism for BDs as they do not follow the same trends that giant exoplanets do with metallicity. We identify a diagonal envelope bounding the Gaia BDs in the mass–period plane due to the detection thresholds of the currently available NSS solutions from 34 months of data.

We present radial velocity measurements of the very bright ($V\sim5.7$) nearby F star, DMPP-4 (HD 184960). The anomalously low Ca II H&K emission suggests mass loss from planets orbiting a low activity host star. Periodic radial velocity variability with $\sim 10$ ms$^{-1}$ amplitude is found to persist over a $>4$ year timescale. Although the non-simultaneous photometric variability in four TESS sectors supports the view of an inactive star, we identify periodic photometric signals and also find spectroscopic evidence for stellar activity. We used a posterior sampling algorithm that includes the number of Keplerian signals, $N_\textrm{p}$, as a free parameter to test and compare (1) purely Keplerian models (2) a Keplerian model with linear activity correlation and (3) Keplerian models with Gaussian processes. A preferred model, with one Keplerian and quasi-periodic Gaussian process indicates a planet with a period of $P_\textrm{b} = 3.4982^{+0.0015}_{-0.0027}$ d and corresponding minimum mass of $m_\textrm{b}\,\textrm{sin}\,i = 12.2^{+1.8}_{-1.9}$ M$_\oplus$. Without further high time resolution observations over a longer timescale, we cannot definitively rule out the purely Keplerian model with 2 candidates planets with $P_\textrm{b} = 2.4570^{+0.0026}_{-0.0462}$ d, minimum mass $m_\textrm{b}\,\textrm{sin}\,i = 8.0^{+1.1}_{-1.5}$ M$_\oplus$ and $P_\textrm{c} = 5.4196^{+0.6766}_{-0.0030}$ d and corresponding minimum mass of $m_\textrm{b}\,\textrm{sin}\,i = 12.2^{+1.4}_{-1.6}$ M$_\oplus$. The candidate planets lie in the region below the lower-envelope of the Neptune Desert. Continued mass loss may originate from the highly irradiated planets or from an as yet undetected body in the system.

The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active research areas from stellar astrophysics and exoplanet studies to Galactic and Local Group archaeology. HRMOS fills a gap in capabilities amongst the landscape of future instrumentation planned for the next decade. The key characteristics of HRMOS will be high spectral resolution (R = 60000 - 80000) combined with multi-object (20-100) capabilities and long term stability that will provide excellent radial velocity precision and accuracy (10m/s). Initial designs predict that a SNR 100 will be achievable in about one hour for a star with mag(AB) = 15, while with the same exposure time a SNR 30 will be reached for a star with mag(AB) = 17. The combination of high resolution and multiplexing with wavelength coverage extending to relatively blue wavelengths (down to 380 nm), makes HRMOS a spectrograph that will push the boundaries of our knowledge and that is envisioned as a workhorse instrument in the future. The science cases presented in this White Paper include topics and ideas developed by the Core Science Team with the contributions from the astronomical community, also through the wide participation in the first HRMOS Workshop (https://indico.ict.inaf.it/event/1547/) that took place in Firenze (Italy) in October 2021.

Precision radial velocity (RV) measurements continue to be a key tool to detect and characterise extrasolar planets. While instrumental precision keeps improving, stellar activity remains a barrier to obtain reliable measurements below 1-2 m/s accuracy. Using simulations and real data, we investigate the capabilities of a Deep Neural Network approach to produce activity free Doppler measurements of stars. As case studies we use observations of two known stars (Eps Eridani and AUMicroscopii), both with clear signals of activity induced RV variability. Synthetic data using the starsim code are generated for the observables (inputs) and the resulting RV signal (labels), and used to train a Deep Neural Network algorithm. We identify an architecture consisting of convolutional and fully connected layers that is adequate to the task. The indices investigated are mean line-profile parameters (width, bisector, contrast) and multi-band photometry. We demonstrate that the RV-independent approach can drastically reduce spurious Doppler variability from known physical effects such as spots, rotation and convective blueshift. We identify the combinations of activity indices with most predictive power. When applied to real observations, we observe a good match of the correction with the observed variability, but we also find that the noise reduction is not as good as in the simulations, probably due to the lack of detail in the simulated physics. We demonstrate that a model-driven machine learning approach is sufficient to clean Doppler signals from activity induced variability for well known physical effects. There are dozens of known activity related observables whose inversion power remains unexplored indicating that the use of additional indicators, more complete models, and more observations with optimised sampling strategies can lead to significant improvements in our detrending capabilities.

ABSTRACT We present additional HARPS radial velocity observations of the highly eccentric (e ∼ 0.6) binary system DMPP-3AB, which comprises a K0V primary and a low-mass companion at the hydrogen burning limit. The binary has a 507 d orbital period and a 1.2 au semimajor axis. The primary component harbours a known 2.2 M⊕ planet, DMPP-3A b, with a 6.67-d orbit. New HARPS measurements constrain periastron passage for the binary orbit and add further integrity to previously derived solutions for both companion and planet orbits. Gaia astrometry independently confirms the binary orbit and establishes the inclination of the binary is 63.89 ± 0.78°. We performed dynamical simulations that establish that the previously identified ∼800 d RV signal cannot be attributed to an orbiting body. The additional observations, a deviation from strict periodicity, and our new analyses of activity indicators suggest the ∼800 d signal is caused by stellar activity. We conclude that there may be long-period planet ‘detections’ in other systems, which are similar misinterpreted stellar activity artefacts. Without the unusual eccentric binary companion to the planet-hosting star, we could have accepted the ∼800 d signal as a probable planet. Further monitoring of DMPP-3 will reveal which signatures can be used to most efficiently identify these imposters. We also report a threshold detection (0.2 per cent FAP) of a ∼2.26 d periodicity in the RVs, potentially attributed to an Earth-mass S-type planet interior to DMPP-3A b.

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution.

Context. Gliese-832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 ± 9.3 days) is close to the orbital period of putative planet c (35.68 ± 0.03 days). Aims. We aim to confirm or dismiss the planetary nature of the RV signature attributed to GJ 832c, by adding 119 new RV data points, new photometric data, and an analysis of the spectroscopic stellar activity indicators. Additionally, we update the orbital parameters of the planetary system and search for additional signals. Methods. We performed a frequency content analysis of the RVs to search for periodic and stable signals. Radial velocity time series were modelled with Keplerians and Gaussian process (GP) regressions alongside activity indicators to subsequently compare them within a Bayesian framework. Results. We updated the stellar rotational period of GJ 832 from activity indicators, obtaining 37.5+1.4-1.5 days, improving the precision by a factor of 6. The new photometric data are in agreement with this value. We detected an RV signal near 18 days (FAP < 4.6%), which is half of the stellar rotation period. Two Keplerians alone fail at modelling GJ 832b and a second planet with a 35-day orbital period. Moreover, the Bayesian evidence from the GP analysis of the RV data with simultaneous activity indices prefers a model without a second Keplerian, therefore negating the existence of planet c.

Current exoplanet surveys using the radial velocity (RV) technique are targeting M dwarfs because any habitable zone terrestrial-mass planets will induce a high RV and orbit on shorter periods than for more massive stars. One of the main caveats is that M dwarfs show a wide range of activity levels from inactive to very active, which can induce an asymmetry in the line profiles and, consequently, a spurious RV measurement. We aim to benchmark the impact of stellar activity on high-precision RV measurements using regular-cadence CARMENES visible and near-infrared observations of the active M3.5 dwarf EV Lac. We used the newly developed technique of low-resolution Doppler imaging to determine the centre-of-light, or spot-induced RV component, for eight observational epochs. We confirm a statistically significant and strong correlation between the independently measured centre-of-light and the chromatic index, which is a measure of the amplitude variation with wavelength of the RVs. We also find circular closed-loop relations of several activity indices with RV for a subset of data that covers only several rotation periods. We also investigate the implications of large phase gaps in the periodograms of activity indicators. Finally, by removing the spot-induced RV component we improve the planet-mass sensitivity by a factor of at least three. We conclude that for active M stars, a regular-cadence observing strategy is the most efficient way to identify and eliminate sources of correlated noise.

Context. Gliese 832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 $\pm$ 9.3 days) is close to the orbital period of putative planet c (35.68 $\pm$ 0.03 days). Aims. We aim to confirm or dismiss the planetary nature of the RV signature attributed to GJ 832c, by adding 119 new RV data points, new photometric data, and an analysis of the spectroscopic stellar activity indicators. Additionally, we update the orbital parameters of the planetary system and search for additional signals. Methods. We performed a frequency content analysis of the RVs to search for periodic and stable signals. Radial velocity time series were modelled with Keplerians and Gaussian process (GP) regressions alongside activity indicators to subsequently compare them within a Bayesian framework. Results. We updated the stellar rotational period of GJ 832 from activity indicators, obtaining $37.5^{+1.4}_{-1.5}$ days, improving the precision by a factor of 6. The new photometric data are in agreement with this value. We detected an RV signal near 18 days (FAP < 4.6%), which is half of the stellar rotation period. Two Keplerians alone fail at modelling GJ 832b and a second planet with a 35-day orbital period. Moreover, the Bayesian evidence from the GP analysis of the RV data with simultaneous activity indices prefers a model without a second Keplerian, therefore negating the existence of planet c.

We present a readout scheme for CMOS image sensors that can be used to achieve arbitrarily high dynamic range (HDR) in principle. The linear full well capacity (LFWC) in high signal regions was extended 50 times from 20 to 984 ke $^{-}$ via an interlaced row-wise readout order, while the noise floor remained unchanged in low signal regions, resulting in a 34-dB increase in DR. The peak signal-to-noise ratio (PSNR) is increased in a continuous fashion from 43 to 60 dB. This was achieved by summing user-selected rows that were read out multiple times. Centroiding uncertainties were lowered when template-fitting a projected pattern, compared to the standard readout scheme. Example applications are aimed at scientific imaging due to the linearity and PSNR increase.

Ariel ’s ambitious goal to survey a quarter of known exoplanets will transform our knowledge of planetary atmospheres. Masses measured directly with the radial velocity technique are essential for well determined planetary bulk properties. Radial velocity masses will provide important checks of masses derived from atmospheric fits or alternatively can be treated as a fixed input parameter to reduce possible degeneracies in atmospheric retrievals. We quantify the impact of stellar activity on planet mass recovery for the Ariel mission sample using Sun-like spot models scaled for active stars combined with other noise sources. Planets with necessarily well-determined ephemerides will be selected for characterisation with Ariel . With this prior requirement , we simulate the derived planet mass precision as a function of the number of observations for a prospective sample of Ariel targets. We find that quadrature sampling can significantly reduce the time commitment required for follow-up RVs, and is most effective when the planetary RV signature is larger than the RV noise. For a typical radial velocity instrument operating on a 4 m class telescope and achieving 1 m s −1 precision, between ~17% and ~ 37% of the time commitment is spent on the 7% of planets with mass M p < 10 M ⊕ . In many low activity cases, the time required is limited by asteroseismic and photon noise. For low mass or faint systems, we can recover masses with the same precision up to ~3 times more quickly with an instrumental precision of ~10 cm s −1 .

Small low-mass stars are favourable targets for the detection of rocky habitable planets. In particular, planetary systems in the solar neighbourhood are interesting and suitable for precise characterisation. The Red Dots campaigns seek to discover rocky planets orbiting nearby low-mass stars. The 2018 campaign targeted GJ 1061, which is the 20$^{\rm th}$ nearest star to the Sun. For three consecutive months we obtained nightly, high-precision radial velocity measurements with the HARPS spectrograph. We analysed these data together with archival HARPS data. We report the detection of three planet candidates with periods of $3.204\pm 0.001$, $6.689\pm 0.005$ and $13.03\pm 0.03$ days, which is close to 1:2:4 period commensurability. After several considerations related to the properties of the noise and sampling, we conclude that a 4$^{\rm th}$ signal is most likely explained by stellar rotation, although it may be due to a planet. The proposed three-planet system (and the potential four-planet solution) is long-term dynamically stable. Planet-planet gravitational interactions are below our current detection threshold. The minimum masses of the three planets range from $1.4\pm 0.2$ to $1.8\pm 0.3$ M$_\oplus$. Planet d, with $m \sin i = 1.68\pm0.25$ M$_\oplus$, receives a similar amount of energy as Earth receives from the Sun. Consequently it lies within the liquid-water habitable zone of the star and has a similar equilibrium temperature to Earth. GJ 1061 has very similar properties to Proxima Centauri but activity indices point to lower levels of stellar activity.

We analyze TESS photometry of DMPP-1 (HD 38677; TIC 66560666), a nearby F8V star hosting hot super-Earth planets and a warm Neptune. Using the Transit Least Squares algorithm and other methods we find a transit signal at days with depth ppm and false alarm probability 1.6%. This is shallower than hitherto published TESS discoveries. The 3.285 day signal is recovered for several, but not all, methods for detrending stellar astrophysical variability. Further observations are needed to improve the significance of the detection. If this transit were due to an Earth-like rocky planet it would have been detected in the RV data, but it is not. The TESS data cover seven individual transits, one of which is consistent with zero depth. The insolation of the putative planet is 990 S⊕, typical of fluxes experienced by the three known catastrophically disintegrating exoplanets (CDEs). The transits can be self-consistently attributed to a CDE with a mass below the RV detection threshold. We searched for transits of the known RV planets, finding null results and detection thresholds of <100 ppm, which we quantify for each. The DMPP-1 planetary system was discovered as a consequence of circumstellar gas attributed to ablation of hot planets. The RV planets may have been ablated to near-pure iron cores. We place limits on the orbital inclinations of the RV planets where the expected transit depth exceeds the detection threshold. If the 3.2854 day transit detection is confirmed, e.g., with CHEOPS photometry, DMPP-1 would be a first-rate target for James Webb Space Telescope spectroscopy.

The closet exoplanets to the Sun provide opportunities for detailed characterization of planets outside the Solar System. We report the discovery, using radial velocity measurements, of a compact multiplanet system of super-Earth exoplanets orbiting the nearby red dwarf star GJ 887. The two planets have orbital periods of 9.3 and 21.8 days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8-day planet is ~350 kelvin. The planets are interior to, but close to the inner edge of, the liquid-water habitable zone. We also detect an unconfirmed signal with a period of ~50 days, which could correspond to a third super-Earth in a more temperate orbit. Our observations show that GJ 887 has photometric variability below 500 parts per million, which is unusually quiet for a red dwarf.

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.