Open Access
Issue
A&A
Volume 684, April 2024
Article Number A83
Number of page(s) 24
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202348813
Published online 05 April 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

One of the most uncommon types of exoplanets discovered to date are ultra-short-period (USP) planets. USP planets are defined arbitrarily as planets with orbital periods shorter than one day (Sahu et al. 2006; Sanchis-Ojeda et al. 2013). Owing to their short planet-to-star distance, they experience levels of irradiation on the order of hundreds or thousands of times stronger than Earth’s, and are thus subject to extreme temperatures. These planets typically have small sizes, with planetary radii up to Rp < 2 R (Winn et al. 2018), and appear to have mostly Earthlike compositions (Dai et al. 2019), some of them even exhibiting iron-enhanced densities (e.g., Rappaport et al. 2013; Santerne et al. 2018; Price & Rogers 2020). USP planets have certain characteristics that put them in a distinct category of exoplanets. For example, unlike hot Jupiters, the presence of USP planets does not strongly correlate with the metallicity of their parent star (Winn et al. 2017). Another characteristic that separates USP planets from hot Jupiters is that the former are often found in multi-planetary systems (Sanchis-Ojeda et al. 2014; Winn et al. 2018). Furthermore, the abundance of USP planets seems to depend on the spectral type of their host star; Sanchis-Ojeda et al. (2014) computed the USP planet occurrence rates for different stellar spectral types and found 0.15 ± 0.05% for F-type stars, 0.51 ± 0.07% for G-type stars, and 1.1 ± 0.4% for M-type stars. Although the measured occurrence of USP planets is higher for M dwarfs than for F dwarfs, that difference is of modest statistical significance (about 2σ). More significant is the higher occurrence around G dwarfs compared to F dwarfs (about 4σ), making it seem reasonable that the occurrence decreases as a function of stellar mass, in general.

The formation pathway of USP planets is still unknown. Because of the high irradiation levels they experience, USP planets typically lie within the region where dust is sublimated, which means that they probably had to form farther away from their current orbit. Proposed formation mechanisms include photo-evaporation of sub-Neptune planets (e.g., Lundkvist et al. 2016; Lee & Chiang 2017) and migration of rocky planets (e.g., Petrovich et al. 2019).

Highly irradiated rocky USP planets may offer a unique opportunity to study the surface and internal composition of exoplanets with current instrumentation. Planets with surface temperatures above ~1100 K will probably have large molten regions since this temperature range marks the melting point of rocks for Earth-like compositions (Katz et al. 2003). This surface temperature range can be achieved by USP planets owing to their short orbital periods and because they are most probably tidally locked to their star. Although USP rocky planets probably lose their primordial atmospheres in the early stages of their formation owing to the effects of the stellar activity and winds of their host star, some volatile elements associated with the magma oceans located in the planet’s dayside may be detected in emission (Ito et al. 2015; Nguyen et al. 2020, 2022).

Currently, no confirmed detection of a USP rocky planet atmosphere has been reported; however, the canonical USP rocky exoplanet 55 Cnc e (P = 0.73 days, Rp = 1.87 R, and Mp = 7.99 M; McArthur et al. 2004; Dawson & Fabrycky 2010; Winn et al. 2011; Demory et al. 2011) is a promising candidate for the detection of a secondary atmosphere with an origin associated with its magma ocean. The depths of the secondary eclipses of 55 Cnc e have shown evidence of variability at both optical (Meier Valdés et al. 2022, Meier Valdés et al. 2023; Demory et al. 2023) and infrared wavelengths (Demory et al. 2016; Tamburo et al. 2018), while the primary transit depths do not exhibit any significant variability (Tamburo et al. 2018; Meier Valdés et al. 2023). Heng (2023) proposed an explanation for 55 Cnc e observations in which a secondary atmosphere is produced by geochemical outgassing events on the dayside of the planet; this gas envelope is transient owing to atmospheric escape processes. In this scenario an eclipse with a depth consistent with zero is produced when there is no secondary atmosphere; when the outgassed atmosphere accumulates in the dayside, this then has the effect of increasing the eclipse depth at optical wavelengths owing to a combination of Rayleigh scattering and thermal emission; and finally fluctuations in temperature are responsible for the variable infrared eclipse depth. Future JWST observations may detect evidence for this proposed scenario for 55 Cnc e; however, more observations are needed to explore whether other rocky USP planets also have variable infrared secondary eclipse depths, such observations perhaps opening the door to the study of secondary atmospheres in exoplanets.

To further improve our understanding of the formation and evolution of USP planets, it is important to increase the sample of this type of object with accurately measured radii and masses. In this paper we report the detection of a USP planet around the star Wolf 327. This star is an M2.5-type star located at d = 28.5 pc. The transit events were discovered by the Transiting Exoplanet Survey Satellite (TESS) as part of its search for transiting planets. This paper is organized as follows: in Sects. 2 and 3 we describe TESS and the ground-based data of this target; Sect. 4 presents the stellar properties of the host star; Sect. 5 describes our data analysis methods and results; and Sect. 6 discusses our results and puts Wolf 327b in the context of known USP planets; and finally in Sect. 7 we present our conclusions.

2 TESS observations

The transiting planet candidate was discovered by TESS (Ricker et al. 2014). TESS takes photometric data continuously for ~27 days and sends data to Earth every ~13.7 days. Wolf 327 (TIC 4918918, TOI-5747) was observed at 2-minute cadence in sector 21 from 2020 January 21 to 2020 February 18 and in sector 48 from 2022 January 28 to 2022 February 26. Once the data are received, the TESS Science Processing Operations Center (SPOC; Jenkins et al. 2016) at the NASA Ames Research Center generates simple aperture-photometry (SAP; Morris et al. 2020) light curves that are then processed to remove systematic effects using the Presearch Data Conditioning (PDC) pipeline module (Smith et al. 2012; Stumpe et al. 2012,Stumpe et al. 2014). The light curves are then scanned for transit-like signals with a wavelet-based adaptive noise-compensating matched filter (Jenkins 2002; Jenkins et al. 2010, Jenkins et al. 2020). A limb-darkened transiting planet model is fitted (Li et al. 2019), and a suite of false-positive detection tests are applied to rule out certain non-planetary scenarios (Twicken et al. 2018). Once this transit search and model-fitting process is completed, the candidate reports are vetted by the TESS Science Office (TSO) at Massachusetts Institute of Technology (MIT), and the planet candidates are announced to the community and assigned a TESS object of interest (TOI) number. On 2022 September 1, TSO announced the detection of a transiting planet candidate around Wolf 327. A transit signature with an orbital period of 0.573 day (~13.7 h) was detected in the SPOC search of the combined light curve from the two sectors on 2022 May 26. The transit signature passed all the diagnostic tests presented in the SPOC Data Validation reports. The SPOC Data Validation process performs difference imaging (out-of-transit baseline image minus in-transit image) and centroid offset analysis as described by Twicken et al. (2018) to establish the location of the transit source with respect to the location of the target star. A mean difference image is produced in each sector (with transits) for a given threshold crossing event and then a centroid based on the pixel response function is computed to determine the location of the transit source. Afterwards a centroid offset is determined by subtracting the location of the target star, corrected for proper motion. Wolf 327 was observed in two sectors (sectors 21 and 48); for targets observed in multiple sectors SPOC computes a mean centroid offset by robustly averaging over the centroid offsets associated with the individual sectors. The mean offset located the source of the transit at 5.13 ± 4.20″ from the coordinates of Wolf 327. There are no other TESS Input Catalog (TIC) objects within 3σ of the mean centroid offset, hence all TIC objects are excluded as potential sources of the transit signal. After passing all the diagnostics tests performed by SPOC, the candidate was assigned TOI number TOI-5747.01 by the TESS Science Office.

In this study, we have analyzed TESS PDCSAP photometry of Wolf 327 from sectors 21 and 48, both taken with a 2-min cadence. Figure 1 shows the TESS target pixel file (TPF) produced with tpfplotter1, the figure shows the field around the target, highlighting in red the apertures used to compute the light curves. TESS apertures were computed by SPOC’s compute optimal aperture, which selects the photometric aperture pixels for each target and sector. The code also calculates a “crowding metric” for each photometric aperture; this is the fraction of light in the photometric aperture due to the target star after background subtraction. In sector 21 the crowding metric was 0.9259, therefore 7.41% of the median flux level was subtracted on each cadence to account for light from other nearby stars in the PDC-SAP light curve. In sector 48 the crowding metric was 0.9060, hence 9.40% of the median flux level was subtracted to account for dilution. Otherwise, the transit depths would be artificially small. This dilution correction is based on a scalar average of each sector and is not dynamic.

Figure 2 shows the TESS light curves for sectors 21 and 48; both time series clearly show the eclipses of the nearby star TIC 4918919 (V = 12.1 mag, J = 10.7 mag). TIC 4918919 is located at a projected separation of r = 55″ from Wolf 327. This nearby star is a known eclipsing binary classified by Gaia (Gaia DR3 796185373590623360), and TESS observes an eclipse every ~3.89 days. This eclipsing binary is a background star (Gaia DR3 parallax π = 1.53 ± 0.02 mas) and is not bound to Wolf 327 (Gaia DR3 parallax π = 35.00 ± 0.03 mas). Although this star is positioned outside TESS’s photometric aperture for both sectors (see Fig. 1) owing to its brightness, part of its point spread function (PSF) is included inside the aperture used to integrate the flux of Wolf 327. As previously mentioned, the constant flux contamination level of TIC 4918919 is taken into account by the PDCSAP pipeline; however, this dilution correction is not dynamic and the eclipses are still seen in the final TESS light curve. In Sect. 5.2, we describe the steps taken to account for the contributions of the transiting candidate and eclipsing binary present in both TESS time series.

thumbnail Fig. 1

Tpfplotter (Aller et al. 2020) target pixel file (TPF) images of Wolf 327 for sectors 21 and 48. The red squares correspond to the TESS aperture used to compute the photometry, the size of the red circles represent the Gaia DR3 (Gaia Collaboration 2023) magnitudes of the stars, and the gray arrows shows the Gaia DR3 proper motion directions of each star inside the field of view. Wolf 327b is positioned at the center (star 1), star 7 (TIC 4918919) is a known eclipsing binary.

3 Ground-based observations

3.1 Seeing-limited photometry

The TESS pixel scale is ~21″ pixel−1 and photometric apertures typically extend out to roughly 1′, generally causing multiple stars to blend in the TESS photometric aperture. To determine the true source of the TESS detection, we acquired ground-based time-series follow-up photometry of the field around Wolf 327 as part of the validation process of the candidate.

3.1.1 Long-term photometric monitoring

We searched for long-term photometric data of Wolf 327 in archival databases. We found a time series taken by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014). ASAS-SN is a survey dedicated to finding bright supernovae and uses an array of 24 robotic telescopes located in several countries (USA, Chile, South Africa, China). The project provides light curves for stars with magnitudes between 9 and 18 in g and V filters through their Sky Patrol portal2. We searched the ASAS-SN portal for Wolf 327 photometry and found ~200 epochs taken with the V filter, spanning a time baseline of ~2060 days. The median values of the photometric measurements are V¯=12.98$\bar V = 12.98$ mag with a median photometric uncertainty of σV = 0.02 mag; the standard deviation of the time series is σdev = 0.016 mag. During the period of time in which the star was visible, the cadence of the observations was roughly one visit every 3 days.

A long term photometric campaign of Wolf 327 was carried out from April to December 2023 with the 0.8 m Joan Oró telescope (TJO; Colomé et al. 2010) at the Montsec Observatory in Lleida, Spain. We obtained a total of 353 images on 61 different nights with an exposure time of 60 s each using the Johnson R filter of the LAIA imager, a 4k×4k with a field of view of 30′ and a scale of 0.″4 pixel−1. The images were calibrated with darks, bias, and flat fields with the ICAT pipeline (Colome & Ribas 2006) of the TJO. The differential photometry was extracted with AstroImageJ (Collins et al. 2017) using the aperture size that minimized the root-mean-square (rms) of the resulting relative fluxes, and a selection of the ten brightest comparison stars in the field that did not show variability. Then, we used our own pipelines to remove outliers and measurements affected by poor observing conditions or with a low signal-to-noise ratio (S/N). The resulting rms of the differential photometry from the TJO in the R filter is ~11 ppt.

Owing to the relatively large photometric uncertainties and the cadence of the observations, it is difficult to detect the transit events in the ASAS-SN and TJO data; however, the long-term monitoring of Wolf 327 provides insights into the level of stellar activity of the star and helps to establish the rotation period of Wolf 327 (see Sect. 4.2).

3.1.2 MuSCAT2 photometry

Wolf 327 was observed on the night of 2022 November 6 with the 1.5 m Telescopio Carlos Sánchez (TCS) at Teide Observatory, Spain. The data were acquired by the TCS MuSCAT2 imager (Narita et al. 2019). This instrument is equipped with four cameras and is capable of obtaining simultaneous images in g′, r′, i′, and z-short with little overhead time.

During the observations the exposure times were set to 10 s for the four bands. Standard data reduction and calibrations were performed by the MuSCAT2 pipeline (Parviainen et al. 2019); in addition to basic reduction, this pipeline also optimizes the photometric aperture while fitting a transit model (including systematic effects) to the time series.

Since the data covered an egress of the transit and not a full event, they were not used in the final fit. However, the MuS-CAT2 observations were useful to detect hints of an egress on Wolf 327, confirming the transit event on the target star and ruling out the presence of hitherto unrecognized eclipsing binaries located inside TESS’s photometric aperture that could mimic a signal similar to that caused by a planetary transit.

thumbnail Fig. 2

TESS PDCSAP photometry for sectors 21 and 48. The deep transit events occurring with a period of ~3.89 days are from the known eclipsing binary TIC 4918919. The solid black line shows the best-fitting model from the joint data fit described in Sect. 5.2. The vertical lines represent the central time of the transit of the USP Wolf 327b (orange) and the eclipsing times of the EB TIC 4918919 (blue).

3.1.3 LCOGT photometry

Wolf 327 was included as part of the TESS Follow-up Observing Program (TFOP; Collins 2019)3. The on-target follow-up light curves are also used to place constraints on the transit depth and the TESS ephemeris. To schedule our transit observations we used the TESS Transit Finder, which is a customized version of the Tapir software package (Jensen 2013).

We observed three full transit windows of Wolf 327b on 2022 November 19, 2022 December 2 and 2022 December 12 in the Pan-STARRS z-short band and two in the Sloan i′ band, respectively, using the Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013) 1.0 m network nodes at Teide Observatory on the island of Tenerife and at McDonald Observatory near Fort Davis, Texas, United States. The images were calibrated with the standard LCOGT BANZAI pipeline (McCully et al. 2018) and differential photometric data were extracted using AstroImageJ. We used circular photometric apertures with a radius of 4.″7 for the z-short band observations and 12.″4 and 5.″8 for the two i′ band observations, respectively, to extract the differential photometry. The smaller target star apertures excluded all of the flux from the nearest known neighbor in the Gaia DR3 catalog (Gaia DR3 796185407950360320), which is ~9″ southwest of Wolf 327. The observational mode on the night of 2022 December 2 used greater defocusing and thus resulted in a larger photometric aperture that includes flux from the nearest known neighbor. However the neighbor is more than a hundred times fainter than Wolf 327 and insignificantly affects the transit depth. The three light curves are included in the global modeling described in Sect. 5.2.

Two full transit windows and two ingress windows of Wolf 327b were observed simultaneously in Sloan g′, r′, i′, and Pan-STARRS z-short bands on 2023 January 19, 2023 January 23, 2023 March 25, and 2023 April 25 using the LCOGT 2m Faulkes Telescope North at Haleakala Observatory on Maui, Hawai’i. The telescope is equipped with the MuSCAT3 multi-band imager (Narita et al. 2020). The images were calibrated using the standard LCOGT BANZAI pipeline (McCully et al. 2018), and photometric data were extracted using AstroImageJ (Collins et al. 2017). We used circular apertures with a radius of ~5.″0 to extract the differential photometry. The two ingress window observations on 2023 January 23 and 2025 March 25 were not included in the global modeling.

3.2 High resolution imaging

As part of our standard process for validating transiting exoplanets to assess the possible contamination by bound or unbound companions on the derived planetary radii (Ciardi et al. 2015), we observed Wolf 327 with high-resolution near-infrared adaptive optics (AO) imaging at Keck Observatory and with optical speckle observations at Sternberg Astronomical Institute.

thumbnail Fig. 3

Adaptive optics and speckle imaging of Wolf 327. Left panel: Keck NIR AO imaging and sensitivity curve for Wolf 327 taken with the Kcont filter. The black points represent the 5σ limits and are separated in steps of 1 FWHM; the purple area represents the 1σ azimuthal dispersion of the contrast determinations (see text). The inset image of the primary target shows no additional close-in companions. Right panel: Ic band speckle sensitivity curve of Wolf 327 taken with the SPeckle Polarimeter (SPP) at the Caucasian Observatory. No stellar companions are detected with magnitudes brighter than ΔIc = 3.9 and 6.2 at ρ = 0.″25 and 1.″0 respectively. The inset shows the speckle autocorrelation function.

3.2.1 Near-infrared AO at Keck

Observations of Wolf 327 were taken with the NIRC2 instrument on Keck-II behind the natural guide star AO system (Wizinowich et al. 2000) on 2023 April 26 in the standard 3-point dither pattern that is used with NIRC2 to avoid the left lower quadrant of the detector, which is typically noisier than the other three quadrants. The dither pattern step size was 3″ and was repeated twice, with each dither offset from the previous dither by 0.″5. NIRC2 was used in the narrow-angle mode with a full field of view of ~10″ and a pixel scale of approximately 0.″0099442 per pixel. The Keck observations were made in the Kcont filter (λo = 2.2706; Δλ = 0.0296 μm) with an integration time in each filter of 20 s for a total of 180 s. Flat fields were generated from a median average of dark-subtracted dome flats. Sky frames were generated from the median average of the nine dithered science frames; each science image was then sky-subtracted and flat-fielded. The reduced science frames were combined into a single combined image using an intra-pixel interpolation that conserves flux and shifts the individual dithered frames by the appropriate fractional pixels; the final resolution of the combined dithers was determined from the full-width half-maximum (FWHM) of the point spread function (0.″100). To within the limits of the AO observations, no stellar companions were detected. The final 5σ limit at each separation was determined from the average of all of the determined limits at that separation and the uncertainty on the limit was set by the root-mean-square dispersion of the azimuthal slices at a given radial distance (Fig. 3, left panel).

3.2.2 Optical speckle at SAI

Wolf 327 was observed on 2023 January 17 with the speckle Polarimeter on the 2.5 m telescope at the Caucasian Observatory of Sternberg Astronomical Institute (SAI) of Lomonosov Moscow State University. The speckle Polarimeter uses a high-speed low-noise CMOS detector Hamamatsu ORCA–quest (Strakhov et al. 2023). The atmospheric dispersion compensator was active, which allowed use of the Ic band. The respective angular resolution is 0.″083, while the long-exposure atmospheric seeing was 0.″8. We did not detect any stellar companions brighter than ΔIc = 3.9 and 6.2 at ρ = 0.″25 and 1.″0, respectively, where ρ is the separation between the source and the potential companion. The sensitivity curve obtained is shown in Fig. (right panel).

3.2.3 Gaia assessment

In addition to employing high-resolution imaging, we have utilized Gaia DR3 to identify possible wide companions bound to Wolf 327. Typically, these stars are already cataloged in the TESS Input Catalog, and their impact on the transit events has been factored into the analysis of the transits and the associated derived parameters. From a comparison of parallaxes and proper motions (e.g., Mugrauer & Michel 2020, 2021; Mugrauer et al. 2022, 2023), our investigation reveals no additional widely separated companions identified by Gaia.

Furthermore, the astrometry from Gaia offers supplementary insights into the potential existence of close companions that might have evaded detection by both Gaia and high-resolution imaging. The Gaia Renormalized Unit Weight Error (RUWE) serves as a metric, akin to a reduced chi-square, with values approximately ≲ 1.4 indicating that the astrometric solution is compatible with a solitary star. Conversely, RUWE values ≳ 1.4 suggest an excess of astrometric noise, potentially attributed to the presence of an unseen companion (e.g., Ziegler et al. 2020). Wolf 327, with a RUWE of 1.3, aligns with this single-star model, suggesting no unseen companions causing excess astrometric noise.

3.3 CARMENES radial velocity measurements

Wolf 327 was observed with the CARMENES spectrograph as part of an RV follow-up program of TESS candidates (22B-3.5-006; PI: E. Pallé). CARMENES is installed at the 3.5m telescope at the Calar Alto Observatory in Almeria, Spain. The instrument has two channels: the visible one (VIS, spectral range 0.52–0.96 μm) and the near-infrared one (NIR, spectral range 0.96–1.71 μm); with an average spectral resolution for both channels of R = 94 600 and R = 80 400, respectively (Quirrenbach et al. 2014, 2018).

We collected a total of 22 CARMENES spectra from 2023 January 6 until 2023 February 6, covering a time baseline of 31 days. The exposure time used to acquire the spectra was 1800 s. The data reduction was done with the CARACAL pipeline (Caballero et al. 2016b), which performs basic data reduction (bias, flat, and cosmic ray corrections) and extracts the spectra using the FOX optimal extraction algorithm (Zechmeister et al. 2014). The wavelength calibration is done following the algorithms described in Bauer et al. (2015). The radial velocity measurements were computed with SERVAL4 (Zechmeister et al. 2018). SERVAL produces a template spectrum by co-adding and shifting the observed spectra and computes the RV shift relative to this template using a χ2 minimization with the RV shift as a free parameter. The RV measurements were corrected for other effects such as barycentric motion, secular acceleration, instrumental drifts, and nightly zero-points (e.g., Trifonov et al. 2018; Tal-Or et al. 2018).

SERVAL also computes several activity indices, such as spectral line indices (the Na I doublet λλ589.0 nm, 589.6nm; Hα λ656.2 nm, and the Ca II infrared triplet λλλ849.8 nm, 854.2 nm, 866.2 nm), and other indicators such as the RV chromatic index (CRX) and the differential line width (dLW). These indicators are useful for monitoring the stellar activity of the star and its effect on the RV measurements (cf. Zechmeister et al. 2018).

The uncertainties of the RV measurements derived by SERVAL have a median value of 1.88 m s−1 and a standard deviation of σdev = 1.43 m s−1 for the VIS channel. In the NIR channel, the median RV uncertainty value is larger at 9.69 m s−1, with a standard deviation of σdev = 4.95 m s−1. The median overall signal-to-noise ratio for the spectra taken with the VIS channel was 329, with a standard deviation of 99 (minimum S/N = 88, maximum S/N = 469). For the NIR spectra, the median S/N was 446, with a standard deviation of 138 (min. S/N = 111, max. S/N = 602).

The SERVAL RV measurements for both channels are shown in Fig. 4. In this study, we used the RV data extracted using the VIS channel, as this dataset presented lower scatter than the NIR measurements (see Bauer et al. 2020 for a comparison between CARMENES VIS and NIR RVs). The CARMENES VIS channel velocity measurements and activity indicators used in this study can be found in Tables A.1 and A.2, respectively.

4 Stellar properties

4.1 Stellar parameters

The coordinates, photometric magnitudes, astrometric values, and stellar parameters of Wolf 327 are listed in Table 1. The stellar parameters were computed using a combined CARMENES stellar spectrum. The effective temperature (Teff), surface gravity (log g), and iron abundance ([Fe/H]) of the star were derived using SteParSyn5 code (Tabernero et al. 2022). SteParSyn makes use of the line list (for the visible and near-infrared wavelength ranges available in CARMENES data) and model grid described by Marfïl et al. (2021). The error in Teff is the quadratic sum of the measurement error given by SteParSyn (8 K) and the weighted average of the spreads of the differences to literature values (69 K) given in Marfïl et al. (2021). The target’s spectral type was determined with ±0.5 dex accuracy from the color-, absolute magnitude-, and luminosity-spectral type relations of Cifuentes et al. (2020). The stellar mass (M*) and radius (R*) were determined following Schweitzer et al. (2019), and the stellar luminosity (L) was computed following Cifuentes et al. (2020). We determine that Wolf 327 is an M2.5 V dwarf with an effective temperature of Teff = 3542 ± 70 K, and a stellar mass and radius of M* = 0.405 + 0.019 M and R* = 0.406 ± 0.015 R, respectively.

We searched for common proper-motion companions to Wolf 327 in the Gaia DR3 catalog up to a radius of 1°. We imposed a restriction on the parallax with a range of 35±5 mas, bracketing the parallax of Wolf 327 (see Table 1). The query returned only Wolf 327, which thus does not seem to share motion and distance with any other star within one degree down to the depth of Gaia. At the distance of Wolf 327 (d = 28.57 ± 1.02 pc), Gaia is sensitive to all stellar companions and even the highest-mass brown dwarfs (Reylé et al. 2021; Lodieu et al. 2019).

The galactocentric velocities (U, V, W) of the star presented in Table 1 were computed following Cortés Contreras (2017), using Gaia DR3 astrometry and adopting a right-handed system with no solar motion correction. The derived Galactic velocities indicate that Wolf 327 belongs to the thin disk (Montes et al. 2001), suggesting an age > 1 Gyr. According to BANYAN Σ6 (Gagné et al. 2018), the U, V, W values suggest that Wolf 327 does not belong to any young moving group, a finding supporting a relatively old age for this star.

thumbnail Fig. 4

Wolf 327 CARMENES radial velocity measurements obtained using SERVAL for the visible (orange circles) and near-infrared (blue squares) channels.

4.2 Stellar rotation from seeing-limited photometry

We used the long-term photometric V-band data from ASAS-SN and the R-band data from TJO to search for periodic flux variations attributable to the rotation period of the star. We attempted to fit the available datasets jointly using a unified modeling approach. However, when adopting a common modeling procedure, the final model for the TJO photometry exhibited clear signs of overfitting. Consequently, we decided to rely solely on the ASAS-SN dataset to establish the rotation period of the star.

In this section, we describe our adopted results for the ASAS-SN dataset, while our results for TJO are presented in Appendix B.

To model the flux variations, we used a combination of a linear function in addition to a periodic term using Gaussian processes (GPs; e.g., Rasmussen & Williams 2010; Gibson et al. 2012) of the form

L(t)=zpt+αt+kij QP,$L(t) = {z_{{\rm{pt}}}} + \alpha t + {k_{ij\,{\rm{QP}}}},$(1)

where zpt and α are the zero point and slope of the linear function, respectively, and kij QP is a GP periodic kernel. For the periodic term, we used the implementation provided by the GP package George (Ambikasaran et al. 2015) and chose the quasi-periodic kernel described by

kij QP=Aexp[ (| titj |)22l2Γ2sin2(π| titj |Prot) ],${k_{ij\,{\rm{QP}}}} = A\,\exp \,\left[ {{{ - {{\left( {\left| {{t_i} - {t_j}} \right|} \right)}^2}} \over {2{l^2}}} - {{\rm{\Gamma }}^2}{{\sin }^2}\,\left( {{{\pi \left| {{t_i} - {t_j}} \right|} \over {{P_{{\rm{rot}}}}}}} \right)} \right]{\rm{,}}$(2)

where |titj| is the difference between two epochs or observations, A is the variance, l is the evolution timescale, Γ is the scale factor and Prot is the period of the sinusoidal variation, that is, the rotational period of star (for a discussion of the physical interpretation of these terms see Nicholson & Aigrain 2022). We also added a jitter term (σphot jitter) as a free parameter to fit the white noise component present in the time series.

The fitting procedure was performed in two steps: first we optimized a posterior probability function using PyDE7 and then we used the optimal set of parameters to start exploring the parameter space using emcee (Foreman-Mackey et al. 2013). We let emcee run for 3000 iterations as a burn-in period and then ran the main Markov chain Monte Carlo (MCMC) procedure for 10 000 iterations. The final values of the fitted parameters were obtained by computing the median and 1σ percentiles of the posterior distributions. Table 2 presents the values and 1σ uncertainties of our fit.

Figure 5 shows the photometric ASAS-SN V-band time series and its periodogram. The left panel shows the median model (black line) computed using the posterior distributions of the fitted parameters; the blue-shaded area shows 1σ uncertainty of the model. The right panel shows the Generalized Lomb-Scargle (GLS) periodogram computed using the python implementation by Zechmeister & Kürster (2009).8 The root-mean-square of the residuals of the fit is rms = 0.013 mag.

The periodogram for the ASAS-SN observations presents a significant peak around 44 days, which we attribute to the stellar rotation of the star. With our fitting procedure we find a stellar rotational period of Prot = 44.4 ± 0.4 days, which will be our final adopted value. This value is in agreement with the results found for the TJO data set (Prot=42.15.8+3.5${P_{{\rm{rot}}}} = 42.1_{ - 5.8}^{ + 3.5}$days, see Appendix B).

Finally, we used the age–rotation period relation described in Angus et al. (2015) to get an estimate of the age of the star. Using our rotation period value (including uncertainties) and Eq. (1) of Angus et al. (2015), we find an age of 4.12.4+3.2$4.1_{ - 2.4}^{ + 3.2}$ Gyr. Considering the uncertainties, this value is in agreement with the estimate obtained using the age-rotation relation for M dwarfs derived by Engle & Guinan (2018); with this relation we obtain an age of 2.7 ± 1.1 Gyr.

Table 1

Stellar parameters of Wolf 327.

thumbnail Fig. 5

Long-term photometric archival data of Wolf 327 from ASAS-SN (V band). Left panel: photometric time-series and best-fitted model computed using GP adopting a quasi-periodic kernel (black line). The 1σ uncertainty regions of the fit are shown in light blue. The residuals of the fit are shown in the bottom panel. Right panel: GLS periodogram (Zechmeister & Kürster 2009) for ASAS-SN V-band photometry. The vertical green line shows the fitted rotational period of the star. The horizontal lines represent the false alarm probability (FAP) levels of 10% (red dotted line) and 1% (blue dashed line).

Table 2

Prior functions and fitted parameters values for the ASAS-SN V-band photometry of Wolf 327.

4.3 Stellar activity indicators

We investigated the spectral line diagnostics of stellar activity across both the VIS and NIR channels of CARMENES. In the CARMENES wavelength range we computed a total of sixteen activity indices following the methods described in Schöfer et al. (2019) and Jeffers et al. (2018). Because the spectra of M dwarfs do not exhibit extensive continuum regions, the spectral line activity indices were computed using the pseudo-equivalent width (pEW) of the lines by integrating the fluxes in regions centered at the position of the line and in two nearby regions (see Table 1 of Schöfer et al. 2019). When dealing with the absorption bands of TiO and VO present in the spectrum, the pEW proves ineffective as a measure for quantifying strength. This inadequacy arises due to the potential extension of these bands across multiple spectral orders, where they can become blended with other spectral features. Instead of using pEWs, the TiO and VO activity levels were computed by measuring the mean fluxes in two regions and taking the ratio between these two quantities (see Table 3 of Schöfer et al. 2019). The same approach was employed for the iron hydride (FeH) Wing-Ford band (Wing & Ford 1969; Schiavon et al. 1997; Schöfer 2021). The regions where the mean fluxes were computed are defined by 989.8 – 990.7 nm and 988.7 – 989.6 nm (see Table 2.3 of Schöfer 2021). The photospheric and chromospheric line activity indices studied here were: log(L/Lbol), pEW(Hα), He D3, Na D, Ca IRT-a,-b,-c, He 1 10830, Paβ, CaH2, CaH3, TiO 7050, TiO 8430, TiO 8860, VO 7436, VO 7942, and FeH Wing-Ford (see Table 3).

In addition to the absorption line indicators, we also investigated the RV chromatic (CRX) and the differential line width (dLW) indices. These indicators serve the purpose of monitoring the stellar activity of the star and assessing its impact on the RV measurements. The CRX index is the slope of the radial velocity (computed in each echelle spectral order) as a function of the logarithm of the wavelength (Zechmeister et al. 2018). The dLW index is used as a proxy to track differential changes in the absorption line widths of the spectrum. This index is computed by comparing the second derivative of the high signal-to-noise spectrum template to the residuals from the fit for the best RVs. This comparison is performed in each echelle spectral order and then averaged taking into account the uncertainties of the spectra (Zechmeister et al. 2018). As is the case for the RVs, we use only the CRX and dLW values derived using the CARMENES VIS channel. Considering all the stellar activity indices, we studied a total of twenty activity diagnostic parameters.

Our results show that Wolf 327 is a relatively stable and Hα inactive star. Stars with positive pEW(Hα) values have Hα lines in absorption in contrast to Hα-active stars with Hα in emission and negative pEW(Hα) values. The low activity level is consistent with the ~44 day rotation period of Wolf 327.

We also investigated whether there was a statistical correlation between the twenty activity diagnostics and the RV values using Pearson’s r, as previously described by Jeffers et al. (2020). None of the 20 diagnostics had a statistically significant correlation with the RV which is defined as Pearson’s r having a value r > |0.7|. Given the low number of RV measurements, any further statistical analysis would be misleading.

5 Analysis and results

5.1 Radial velocity and activity indices periodogram analysis

We computed the GLS periodogram of the CARMENES radialvelocity and activity-indicator measurements for Wolf 327. Figure 6 shows the periodograms for the CARMENES RV measurements (VIS channel), the RVs after subtracting a linear trend and the planetary signal (see Sect. 5.2), and the activity indices provided by SERVAL, the vertical black line marks the orbital period of the transiting candidate around Wolf 327.

The RV periodogram shows a broad peak consistent with the period of the transiting candidate detected by TESS. The peak is slightly above the 1% false alarm probability (FAP) level, although this peak is not the one with the highest power in the periodogram. The RV signature of the rotation period of the star is not present in any of the periodograms, probably because the RV measurements were taken with a time baseline shorter than the rotation period of the star (31 days versus Prot = 44.4 ± 0.4 days from photometric measurements) and the observation strategy, as the majority of the RV data acquisition was made during two full nights of continuous observations.

Table 3

Wolf 327 activity indicators computed from a coadded CARMENES VIS and NIR channel spectrum.

5.2 TESS, ground-based light curves, and CARMENES data

To derive the planetary orbital parameters we decided to use TESS, ground-based light curves, and CARMENES RV measurements to perform a joint fit of the data. We use the same method as described in previous papers (e.g., Murgas et al. 2022, 2023). To model the transits, we use PyTransit9 (Parviainen 2015) and we adopted a quadratic limb darkening (LD) law. The fitted LD coefficients were compared through a likelihood function to the predicted values computed by LDTK10 (Parviainen & Aigrain 2015) using the stellar parameters presented in Table 1. During the fitting process we used the LD parameterization of Kipping (2013). The CARMENES radial velocity time series was modeled using RadVel11 (Fulton et al. 2018).

For the analysis, we selected only ground-based observations that covered a full transit of the planet candidate and presented a detection of the transit event (for details see Appendix C). In the end we used the observations taken on 2022 November 19, 2022 December 2, and 2022 December 12 by LCOGT 1 m telescopes; and the 2023 January 19 and 2023 April 25 time series taken with LCOGT 2 m telescope using the MuSCAT3 instrument. We used ground-based light curves that were corrected by observational and instrumental effects (i.e., detrended light curves) provided by the ExoFOP team.

To model the systematic effects of unknown origin present in TESS and ground-based photometry we used GPs. For all the light curves analyzed here we adopted the commonly used Matérn 3/2 kernel:

kij=ck2(1+3| titj |τk)exp(3| titj |τk),${k_{ij}} = c_k^2\,\left( {1 + {{\sqrt 3 \,\left| {{t_i} - {t_j}} \right|} \over {{\tau _k}}}} \right)\,\exp \left( { - {{\sqrt 3 \left| {{t_i} - {t_j}} \right|} \over {{\tau _k}}}} \right),$(3)

where |titj| is the time between epochs in the series, and the hyperparameters, ck and τk, were allowed to be free, with k indicating the different time series (i.e. individual TESS sectors, ground-based observations of transit events).

We performed an initial test, fitting TESS data and CARMENES measurements simultaneously using RV models with and without GPs, and a zero-point offset for each observing night. We found that all the methodologies tested deliver similar results in the determination of the radial velocity semi-amplitude (see Fig. D.1). For the case in which we included GPs in the RVs to model the systematic noise present in the measurements, we adopted a Gaussian kernel of the form

kij RV=crv2exp((titj)2τrv2),${k_{ij\,{\rm{RV}}}} = c_{rv}^2\,\exp \,\left( { - {{{{\left( {{t_i} - {t_j}} \right)}^2}} \over {\tau _{rv}^2}}} \right),$(4)

where titj is the time between epochs in the series, and the hyperparameters, crv and τrv, were set free. We decided to use this kernel since it has few free parameters and assumes that points closer in time will be more correlated than observations separated by a long time baseline. With this approach we expect to model the effect on the RVs of transient stellar activity events on the star. The use of a more complex kernel, such as the popular quasi-periodic kernel, is not warranted since the RV periodogram did not show any evidence of a strong correlation with the rotation period of the star.

To model the planetary transits and RV variations we set as free parameters the planet-to-star radius ratio (Rp/R*), the LD coefficients (q1, q2 following Kipping 2013), the central time of the transit (Tc), the planetary orbital period (P), the stellar density (ρ*), the transit impact parameter (b), the RV semi-amplitude (KRV), the host star systemic velocity (γ0), the long-term linear trend of the radial velocity variation (γ˙${\dot \gamma }$) and the RV jitter (σRV jitter). Owing to the short orbital period of the candidate, we decided to assume a circular orbit and fix the eccentricity to 0. As shown in Fig. 2, TESS light curves show deep eclipses caused by the nearby star TIC 4918919 (see Sect. 2); we model the contribution of this EB with its own light curve model by setting as free parameters the transit depth, period, central time of transit, stellar density, and impact parameter. The LD coefficients for TIC 4918919 were fixed to the values computed with LDTK using the TESS Input Catalog (Stassun et al. 2018, 2019) stellar parameter values for this star.

The fitting procedure started with the optimization of a posterior probability function with PyDE. After the optimization finished, we used the optimal parameter vector as a start of an MCMC procedure using emcee to explore the parameter space.

We ran the MCMC using a number of chains that were at least 5.5 times the number of free parameters for 35 000 iterations as a burn-in stage to ensure convergence. We then ran the main MCMC for another 15 000 iterations. We computed the final values of the fitted parameters using the median and 1σ limits of the posterior distributions.

We compared the results of different modeling approaches for the joint fit of all the datasets available (i.e. TESS, ground-based photometry; CARMENES RVs) through the Bayesian information criterion metric (BIC, Schwarz 1978). Because of the photometric variability present in TESS and ground-based light curves, we decided to keep the use of GPs for these datasets. In the case of the RV measurements we tested models without and with GPs. We computed the BIC metric for three circular orbit models using different modeling approaches for the RVs: Model 1) RVs without GPs, Model 2) RVs including GPs, and Model 3) RVs using the zero-point offset method (no GPs). The computed BIC differences between the models were ΔBICM1M2=10.4${\rm{\Delta BI}}{{\rm{C}}_{{{\rm{M}}_1} - {{\rm{M}}_2}}} = - 10.4$ and ΔBICM1M3=74.4${\rm{\Delta BI}}{{\rm{C}}_{{{\rm{M}}_1} - {{\rm{M}}_3}}} = - 74.4$, meaning that Model 1 presented the lowest BIC value. Therefore, according to the BIC metric, the best model was the simplest one, i.e. transit analysis including GPs for TESS and ground-based observations and no GPs for the RV measurements. Hence, we adopted the parameter values of Model 1 as our final results. Using this approach we had a total of 50 free parameters with which to model the planet candidate (transits and RVs), the eclipsing binary TIC 4918919 in TESS light curves, and systematic effects.

Table 4 presents the median and 1σ uncertainty range for the fitted and derived parameters for Wolf 327b. In Fig. D.2 shows the posterior distributions of the transit and RV measurements fitted parameters (without the GP hyperparameters for easy viewing). Figures 79 present the best-fitting model for TESS and ground-based light curves, and RV measurements, respectively. Figure 10 shows the phase-folded TESS data and best-fitting model for the eclipsing binary TIC 4918919. The radial velocity residuals (see Fig. 9 middle panel) shows a linear trend with time. This trend may be caused by another massive object in the system; however, the relatively short baseline of the RV observations (31 days) does not allow us to put limits on the mass of the object. The potential presence of other planets in Wolf 327 system is in agreement with what is known of the USP population (Sanchis-Ojeda et al. 2014; Winn et al. 2018); however, more RV observations are needed to confirm the presence of more objects in the system and rule out certain stellar activity induced signal that may be the cause of this linear trend.

Using our fitted values and the stellar parameters presented in Table 1, we find that Wolf 327b has a radius of Rp = 1.24 ± 0.06 R and a mass of Mp = 2.53 ± 0.46 M. This planet orbits its host star with a period of Ρ = 0.57347 days (13.76 h) and has an equilibrium temperature of Teq = 996 ± 22 K.

thumbnail Fig. 6

GLS periodograms of the CARMENES RV measurements taken with the optical channel, RV residuals (i.e., linear trend and planet signal removed), spectral activity indices, and window function. The horizontal lines represent the FAP levels of 10% (purple dotted line) and 1% (blue dash line). The vertical black line marks the orbital period of Wolf 327b (P = 0.5734 day).

thumbnail Fig. 7

TESS phase-folded light curve of Wolf 327b after subtracting the photometric variations from the time series (top panel). The best-fit transit model is shown in black, the circles are TESS binned observations. Bottom panel: residuals of the fit.

Table 4

Prior functions, fitted, and derived parameters of Wolf 327b.

thumbnail Fig. 8

Ground-based individual transit observations of Wolf 327b used in the joint fit, systematic effects were removed from the light curves. The points show the individual observations, the circles represent binned data points, and the best-fit model is shown by the black line.

thumbnail Fig. 9

Wolf 327b radial velocity measurements taken with CARMENES. (a) RV time series and best-fitting model (blue line). The best-fitting model was computed using the median values of the posterior distribution of the fitted parameters, (b) Residuals of the fit after subtracting the single-planet Keplerian, the velocity linear trend seen here was included in the RV modeling, (c) Phase-folded RV measurements.

6 Discussion

6.1 Comparison with other exoplanets

Once we confirmed the planetary nature of Wolf 327b, we studied how it compares to the known population of USP planets. Figure 11 shows a period versus radius diagram of known exoplanets with short orbital periods (P < 1 day), similar to the figures presented in the studies of Adams et al. (2021) and Morello et al. (2023). The planets with radius and mass measurements are represented with circles, the size and color of the symbols being proportional to the mass of the planet; planets with only radius measurements are marked by the green squares. The figure highlights the range of sizes presented by USP planets: from Jupiter-sized planets, to a less populated sub-Jovian region, and with a more populated region inhabited by small (Rp < 3 R) planets. The radius of Wolf 327b (Rp = 1.24 + 0.06 R) places it well below the sub-Jovian desert and in a region where the typical mass range goes up to a few tens of Earth masses, in agreement with our measured planetary mass.

Figure 12 shows the mass and radius values of known transiting exoplanets (data taken from TEPcat12 catalog; Southworth 2011). In the figure we also show the composition models of Zeng et al. (2016, 2019). As expected for small USP planets, Wolf 327b’s position in this diagram indicates that this planet probably has a rocky composition, with little to no atmosphere, and with a higher density than that corresponding to an Earth-like composition for the planet’s range of mass and radius.

In terms of mass and radius, the USP planet with the most resemblance to Wolf 327b is K2-229b. This planet was discovered by Santerne et al. (2018), who measured an orbital of period of Ρ = 0.58 day (13.9 h) and estimated a planetary radius of Rp = 1.17 + 0.07 R, a mass of Mp = 2.59 + 0.43 Μ and a bulk density of ρp = 8.9 + 2.1 g cm−3. K2-229b is the inner planet of a planetary system, with a second planet in a 8.32 days orbit (K2-229C, Rp = 2.12 + 0.11 R, Mp < 21.3 Μ) and a possible third transiting candidate in a 31 days orbit (K2-229d, Rp = 2.65 + 0.24 R, Mp < 25.1 Μ). Incidentally, given the trend observed in the CARMENES RV data, Wolf 327b might also be part of a multiplanet system. Although K2-229b is similar in terms of size and mass, this planet has an equilibrium temperature of Teq = 1960 + 40 Κ because it is orbiting a star hotter than Wolf 327, a K0 dwarf (Teff = 5185 + 32 K). Santerne et al. (2018) modeled the planetary interior composition of K2-229b with only a metallic core and a silicate mantle, and found a core-mass fraction of 68%, making K2-229b a planet with an enlarged core in comparison with its mantle, similar to Mercury. In the following section we explore the composition of Wolf 327b using theoretical planetary interior models.

thumbnail Fig. 10

Top panel: TESS phase-folded light curve of the eclipsing binary TIC 4918919. The photometric variations modeled by the GPs were removed from the time series. The best-fit transit model is shown in black, the circles are TESS binned observations. Bottom panel: residuals of the fit.

6.2 Planet interior modeling

We utilized ExoMDN13 (Baumeister & Tosi 2023) to model the interior of Wolf 327b, using as input our derived values for the planetary mass, radius, and equilibrium temperature while taking into account their uncertainties. ExoMDN is an exo-planet interior inference model that adopts a machine-learning approach. The program uses a set of mixture density networks (MDN) trained on 5.6 million synthetic planet models computed with the ΤΑΤ00ΙΝΕ code (Baumeister et al. 2020; MacKenzie et al. 2023). The planets from the training set have masses below 25 Μ, equilibrium temperatures in the range 100–1000 K, and have compositions consisting of an iron core, a silicate mantle, a water and high-pressure ice layer, and an H/He atmosphere. Since our derived equilibrium temperature for Wolf 327b (Teq = 996 ± 23 K) places the planet at the edge of temperature range of ExoMDN’s training dataset, we recommend that the results obtained by this model should be considered just as an approximation.

Figures 13 and 14 present the thickness and mass fraction of the interior layers of Wolf 327b computed with ExoMDN. The interior model predicts that the planet is dominated in terms of layer size and mass fraction by a large planetary iron core; Wolf 327b’s core could take up as much as 78% of the planet’s size and 93% of its mass. For the other layers the model predicts a relatively small mantle, and a negligible fraction of water and gas layers. The predicted lack of an atmosphere is in agreement with what is expected for highly irradiated planets of relatively small sizes. Interestingly, ExoMDN’s prediction of a relatively large planetary core is in agreement with the simulations made for K2-229b in Santerne et al. (2018), thus confirming the similarities between these two planets.

thumbnail Fig. 11

Orbital period versus radius diagram for ultra-short-period planets. Confirmed planets with only radius measurements are shown with green squares, and planets with mass and radius values are shown with circles. The color and size of the circles represent the mass of the planet. The position of Wolf 327b is marked by the red star. Data taken from NASA Exoplanet Archive (https://exoplanetarchive.ipac.caltech.edu/).

thumbnail Fig. 12

Mass-radius diagram for Wolf 327b (red star) and known transiting planets with mass determinations with a precision better than 30% (parameters taken from the TEPcat database; Southworth 2011). Planets orbiting M-type stars (Teff ≤ 4000 K) are marked with circles, orange indicate planets with periods P > 1 day, and red circles indicate USP planets (P ≤ 1 day) around M dwarfs. The lines in the mass–-radius diagram represent the composition models of Zeng et al. (2016, 2019) for planets with pure iron cores (100% Fe, solid brown line), Earth-like rocky compositions (32.5% Fe plus 67.5% MgSiO3, dashed green line) and pure rock (100% MgSiO3, dash-dotted purple line).

6.3 Prospects for atmospheric characterization

Studying the prerequisites for habitability and the orbital limits of the habitable zone requires the characterization of both temperate and hot planets (e.g., Turbet et al. 2019; Schlecker et al. 2023). With Sp = 233.9 ± 18.3 S, Wolf 327b is well within the latter regime (Kopparapu et al. 2013, 2014). If the planet accreted significant amounts of water or other volatiles during its formation, it probably went through a runaway greenhouse phase with an inflated atmosphere (Ingersoll 1969; Kasting 1988; Turbet et al. 2020; Mousis et al. 2020), in which water is removed through photolysis and subsequent hydrodynamic escape of the hydrogen (Wordsworth & Pierrehumbert 2013, 2014). The expected past evolution of the planet’s instellation (Baraffe et al. 2015) makes it likely that the planet never left this phase before being completely desiccated (Luger & Barnes 2015). Our density estimate of Wolf 327b (see Fig. 12) is consistent with a post-runaway greenhouse state without an extended atmosphere (Turbet et al. 2020), which may have been further depleted through intense irradiation at the orbital distance of the planet (Owen & Wu 2013).

To evaluate the possibility of further detailed studies of this planet we computed the transmission spectroscopy metric (TSM) and emission spectroscopy metric (ESM) proposed by Kempton et al. (2018). Both metrics are used to evaluate the feasibility of detecting atmospheric signals during a primary transit (TSM) or detecting the eclipse (ESM) with JWST. For Wolf 327b we find a TSM of 13, Kempton et al. (2018) recommends a threshold of TSM > 10 for planets with radii of Rp < 1.5 R which marks this planet as a potential candidate for future transmission spec-troscopy studies. However, owing to the high levels of irradiation suffered by Wolf 327b, it is possible that this planet has lost its original gas envelope. As for the possibility of eclipse detection, we find an ESM of 12.5; the recommended ESM threshold for rocky planets is of 7.5 which indicates that this planet is also a good candidate for eclipse spectroscopy with JWST. Both emission and transmission JWST observations would be a good way to test the hypothesis that the planet has little or no atmosphere.

At the moment of submission, no JWST results on USP rocky planets had been reported. Early JWST transmission spec-troscopy studies have focused on small planets with orbital periods longer than 1 day and equilibrium temperatures below 1000 K. These planets are: TRAPPIST-1b (P = 1.51 days, Rp = 1.11 R, and Teq ≈ 400 K; Greene et al. 2023; Ih et al. 2023), TRAPPIST-1c (P = 2.42 days, Rp = 1.09 R, and Teq ≈ 340 K; Zieba et al. 2023; Lincowski et al. 2023), LHS 475b (P = 2.03 days, Rp = 0.99 R, and Teq ≈ 590 K; Lustig-Yaeger et al. 2023), and GJ 486b (P = 1.46 days, Rp = 1.30 R, and Teq ≈ 700 K; Moran et al. 2023). These planets are not considered to be USP planets and were expected to have some gaseous envelope. The observations for TRAPPIST-1b and 1c are consistent with a thin atmosphere or a bare rock surface, for LHS 475b JWST found a featureless transmission spectrum, and for GJ 486b the transmission spectrum can be explained by either the presence of water or by the effects of unocculted spots in the transmission spectrum. As discussed in Sect. 1, the USP rocky planet 55 Cnc e is a prime candidate for JWST observations and has two approved Cycle 1 programs. Brandeker et al. (2021) will study whether the planet is in a 3:2 spin-orbit resonance to explain the secondary eclipse depth variability, and Hu et al. (2021) will try to detect any volatiles in the atmosphere of 5 Cnc e using the spectral features of H2O, CO, CO2, and SiO.

During the reviewing process of this article, Zhang et al. (2024) announced the first JWST results on a USP planet. The target of the study is GJ 367b (Lam et al. 2021), a USP sub-Earth with a radius and mass of Rp = 0.699 ± 0.024 R and Mp = 0.633 ± 0.050 R (Goffo et al. 2023), respectively, and an orbital period of 0.322 day (7.7 h). Zhang et al. (2024) measured the phase curve of GJ 367b in the mid-infrared wavelength range (5–12 μm) using the Mid-Infrared Instrument (MIRI; Kendrew et al. 2015) on board of JWST. The team measured an eclipse depth of 79 ± 4 ppm, with a lack of any detectable atmosphere, and featuring a dark surface (i.e., low albedo) within the MIRI bandpass. The measured emission spectrum is consistent with a blackbody spectrum. The measured low albedo value for GJ 367b is consistent with what is expected for lava worlds with molten surfaces (Essack et al. 2020).

thumbnail Fig. 13

Predicted thickness of the interior layers of Wolf 327b computed with ExoMDN.

thumbnail Fig. 14

Predicted mass fraction of the interior layers of Wolf 327b computed with ExoMDN.

6.4 Planet-star interaction

In this section, we study the prospects for detecting radio emission from magnetic star-planet interaction between Wolf 327 and its newly discovered planet. The phenomenon behind this emission is the electron cyclotron maser (ECM) instability (Melrose & Dulk 1982), which can generate auroral radio emission in both the planet and its host star. The characteristic frequency of this emission is given by the electron gyrofrequency, vG = 2.8 B (in MHz), with B being the local magnetic field (in G). The resulting nonthermal emission has a coherent nature, with a high degree of circular polarization (close to 100% in some instances), and can have a large bandwidth (Δν ~ vG/2).

We note that any potential radio emission coming directly from the planet is most probably impossible to detect. Since the magnetic field of a planet is merely a few gauss, any resulting emission would fall well below the ionosphere cut-off. However, if the planet is in the sub-Alfvénic regime, where the speed of the plasma wind relative to the planet is smaller than the Alfvén speed, then energy and momentum can be transferred back to the star via Alfvén wings, and coherent electron-cyclotron emission can originate close to the star. Since the global magnetic fields of M-dwarf stars are typically larger than 100 G, sometimes reaching several kG, the electron gyrofrequency would be of a few hundred MHz, or even a few GHz, making the resulting auroral emission potentially detectable with current radio telescopes.

We estimated the stellar magnetic field of Wolf 327 using the relations in Reiners et al. (2022), which depend on the value of the Rossby number, Ro. From the two different Ro estimates in Wright et al. (2018), we get a magnetic field of the star, B* ≈ 292 G. To estimate the radio emission arising from star–planet interaction, we used the models in Pérez-Torres et al. (2021), where we consider two different magnetic field geometries: a closed dipole and an open Parker spiral. We assumed that the stellar wind was isothermal, with T = 3 × 106 K, and that the solid angle covered by the ECM emission was 1.6 steradians. We considered two different emission models, Saur/Turnpenney s and Zarka/Lanza s, which are respectively shown as orange and blue in the plots. For each model we consider efficiency values ranging from 0.2 to 2%. Since the planet is very close to the star, the interaction does indeed take place in the sub-Alfvénic regime (see top panels of Fig. 15).

We show in Figs. 15a,b the radio flux density from sub-Alfvénic star-planet interaction as a function of the mass-loss rate of the star and find that the emission is strong and very similar in both a close dipolar and an open Parker spiral geometry of the magnetic field. Assuming a conservative threshold value of 100 μJy for the rms obtained in a given observation, we should be able to detect radio emission from star–planet interaction, regardless of the geometry, efficiency, model, or magnetic field of the planet.

thumbnail Fig. 15

Flux density arising from starplanet interaction as a function of the stellar mass-loss rate (panels a and b) and the magnetic field of the planet (panels c and d). (a) Expected flux density for a closed dipolar geometry. (b) Expected flux density for an open Parker spiral geometry. (c) Expected flux density for a closed dipolar geometry. (d) Expected flux density for an open Parker spiral geometry. For (c) and (d) we assume = 0.1 . Note that the resulting radio emission is unaffected by the magnetic field of the planet, unless the latter is very large (~4 G for closed dipole or even > 10 G for a Parker spiral). The emission expected from Saur/Turnpenney’s model is shown in orange, and the emission expected from Zarka/Lanza’s model is shown in blue. The overlap between both models is shown in brown. The dashed black line represents the assumed detection threshold of 100 μJy.

7 Conclusions

We report the discovery of Wolf 327b, an ultra-short-period super-Earth orbiting around an M2.5 V dwarf with a period of 0.573 day. The detection of the transit events was led by NASA’s TESS mission, and ground-based observations confirmed that the transit events occur on the star Wolf 327. Follow-up RV measurements taken with the CARMENES spectrograph allowed us to establish the mass of the candidate and confirm its planetary nature.

From a coadded CARMENES mean stellar spectrum we determined stellar parameters for the host star. We estimate that Wolf 327 has an effective temperature of Teff = 3542 ± 70 K, an iron abundance of [Fe/H] = −0.17 ± 0.04 dex, a derived stellar mass of M* = 0.405 ± 0.019 M, and a derived stellar radius of R* = 0.406 ± 0.015 R. Using long-term ground-based photometric observations, we estimate a stellar rotation period of Prot = 44.4 ± 0.4 days. Gaia DR3 measurements indicate the host star does not present any proper-motion companions, based on the star’s Galactic velocity components it is likely that it does not belong to any young moving group association.

We analyzed photometric data taken by TESS (sectors 21 and 48), five ground-based individual full-transit observations, and 22 RV measurements taken with the CARMENES spectrograph. All the datasets were fitted simultaneously using an MCMC procedure that included the contribution of red noise sources for each time series modeled by Gaussian processes. We note that after subtracting the signal of the USP planet from the RV measurements, the residuals present a slope that might indicate the presence of an additional companion in the system. From our joint data fit we find that Wolf 327b is a super-Earth with a radius of Rp = 1.24 ± 0.06 R and mass of Mp = 2.53 ± 0.46 M. These measurements imply that this planet has an iron-like bulk density of ρp = 7.24 ± 1.66 g cm−3. Owing to its short orbital period (P = 0.57347 day), the planet has an equilibrium temperature of Teq = 996 ± 22 K. Given its density and equilibrium temperature, it is likely that Wolf 327b possesses a rocky composition with little to no atmosphere.

In terms of its mass and radius Wolf 327b, is very similar to the known transiting planet K2-229b. Previous studies of K2-229b indicate that this planet probably possesses a large iron core, with an internal structure similar to that of Mercury. Planet interior modeling applied to Wolf 327b arrives at a similar result: the internal structure of this planet is probably characterized by a substantial iron core, which dominates both in terms of layer size and mass fraction. Surrounding this core is a comparatively small mantle layer, accompanied by a negligible atmosphere (if any).

Overall, Wolf 327b is an interesting addition to the known USP planets around M dwarfs. This planet presents an interesting internal structure and favorable metrics for secondary transit detection with the JWST. Furthermore, since USP planets are usually found in multi-planetary systems, Wolf 327 also presents a good opportunity for additional RV follow-up to search for planets in wider orbits (as the data so far seem to indicate).

Acknowledgements

CARMENES is an instrument at the Centro Astronómico Hispano en Andalucía (CAHA) at Calar Alto (Almería, Spain), operated jointly by the Junta de Andalucía and the Instituto de Astrofísica de Andalucía (CSIC). CARMENES was funded by the Max-Planck-Gesellschaft (MPG), the Consejo Superior de Investigaciones Científicas (CSIC), the Ministerio de Economía y Competitividad (MINECO) and the European Regional Development Fund (ERDF) through projects FICTS-2011-02, ICTS-2017-07-CAHA-4, and CAHA16-CE-3978, and the members of the CARMENES Consortium (Max-Planck-Institut für Astronomie, Instituto de Astrofísica de Andalucía, Landessternwarte Königstuhl, Institut de Ciències de l’Espai, Institut für Astrophysik Göttingen, Universidad Complutense de Madrid, Thüringer Landessternwarte Tautenburg, Instituto de Astrofísica de Canarias, Hamburger Sternwarte, Cen-tro de Astrobiología and Centro Astronómico Hispano-Alemán), with additional contributions by the MINECO, the Deutsche Forschungsgemeinschaft (DFG) through the Major Research Instrumentation Programme and Research Unit FOR2544 “Blue Planets around Red Stars”, the Klaus Tschira Stiftung, the states of Baden-Württemberg and Niedersachsen, and by the Junta de Andalucía. 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 makes use of observations from the LCOGT network. Part of the LCOGT telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF. The Joan Oró Telescope (TJO) of the Montsec Observatory (OdM) is owned by the Catalan Government and operated by the Institute for Space Studies of Catalonia (IEEC). This work is partly supported by JSPS KAKENHI Grant Number JP18H05439, JST CREST Grant Number JPMJCR1761. This article is based on observations made with the MuSCAT2 instrument, developed by ABC, at Telescopio Carlos Sánchez operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide. This paper is based on observations made with the MuSCAT3 instrument, developed by the Astrobiology Center and under financial supports by JSPS KAKENHI (JP18H05439) and JST PRESTO (JPMJPR1775), at Faulkes Telescope North on Maui, HI, operated by the Las Cumbres Observatory. Some of the data presented herein were obtained at Keck Observatory, which is a private 501(c)3 nonprofit organization operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the Native Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. This research has made use of the Exoplanet Follow-up Observation Program (ExoFOP; DOI: 10.26134/Exo-FOP5) website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. We acknowledge financial support from the Agencia Estatal de Investigatión (AEI/10.13039/501100011033) of the Ministerio de Ciencia e Innovatión and the ERDF “A way of making Europe” through projects PID2022-137241NB-C4[1:4], PID2021-125627OB-C3[1:2], PID2019-109522GB-C5[1:4], PID2019-107061GB-C61, and RYC2021-031640-I, and the Centre of Excellence “Severo Ochoa” and “María de Maeztu” awards to the Insti-tuto de Astrofísica de Canarias (CEX2019-000920-S), Instituto de Astrofísica de Andalucía (CEX2021-001131-S) and Institut de Ciències de l’Espai (CEX2020-001058-M). This work was also funded by the Generalitat de Catalunya/CERCA programme; the M.V. Lomonosov Moscow State University Program of Development; the Massachusetts Institute of Technology through the TESS mission via subaward s3449; the University of La Laguna under the EU Next Generation funds through the Margarita Salas Fellowship from the Ministerio de Universidades ref. UNI/551/2021-May 26; the DFG through the priority program SPP 1992 “Exploring the Diversity of Extrasolar Planets” (JE 701/5-1); the DFG through the Major Research Instrumentation Programme and Research Unit FOR2544 “Blue Planets around Red Stars”; the Jet Propulsion Laboratory, California Institute of Technology, under contract 80NM0018D0004 with the National Aeronautics and Space Administration; the Francqui Foundation and the F.R.S.-FNRS through grant T.0109.20; and IPAC through a Visiting Research Fellowship. The results reported herein benefited from collaborations and/or information exchange within NASA’s Nexus for Exoplanet System Science (NExSS) research coordination network sponsored by NASA’s Science Mission Directorate under Agreement No. 80NSSC21K0593 for the program “Alien Earths”.

Appendix A Radial velocity data

Table A.1 present the RV measurements of Wolf 327 taken with the visible channel of CARMENES. Table presents A.2 the line activity indicators computed with SERVAL.

Table A.1

CARMENES VIS channel corrected radial velocities, instrumental drift, and nightly zero-point values computed by SERVAL. The RVs presented here were corrected by instrumental drifts and nightly zero-point offsets.

Table A.2

CARMENES spectroscopic line activity indicators for Wolf 327 computed by SERVAL.

Appendix B Stellar rotation from TJO photometry

thumbnail Fig. B.1

Long-term photometric observations of Wolf 327 taken with the Joan Oró telescope (TJO). Left panel: normalized R-band flux (orange points) versus time. The best-fitted model is shown by the black line, the 1σ uncertainty regions of the fit are shown in light blue. Right panel: GLS periodogram (Zechmeister & Kürster 2009) for TJO R-band photometry. The vertical green line shows the best-fitted rotation period (Prot=42.15.8+3.5 d${P_{{\rm{rot}}}} = 42.1_{ - 5.8}^{ + 3.5}\,{\rm{d}}$). The horizontal lines represent the false alarm probability (FAP) levels of 10% (red dotted line) and 1% (blue dashed line).

In Fig. B.1 we show the ground-based photometric observations of Wolf 327 taken with the 0.8 m Joan Oró telescope (TJO). On the left panel we present the normalized R-band flux versus time, in the right side panel the we show the GLS periodogram of the time-series. The periodogram presents a maximum peak at Pmax = 40.6 d. We attempted to fit the light curve using the quasi-periodic kernel described in #2 Sec. 4.2, however we found that by using this procedure the kernel was overfitting our data. We decided to model the light curve using a combination of a linear and a sinusoidal function:

L(t)=zpt+αt+Asin(2πtP+ϕ),$L(t) = {z_{{\rm{pt}}}} + \alpha t + A\sin \left( {{{2\pi t} \over P} + \phi } \right),$(B.1)

where t is the epoch of the measurements, the free parameters zpt and α are the zero point and slope of the linear function, respectively. For the sinusoidal function the free parameters are the amplitude (A), the rotation period (Prot), and a phase angle (φ). To model the red noise present in the data, we used the GP implementation by Celerite (Foreman-Mackey et al. 2017) and chose a simple exponential kernel of the form

kijexp=a2expc(| titj |),${k_{ij\,\exp }} = {a^2}{\exp ^{ - c\left( {\left| {{t_i} - {t_j}} \right|} \right)}},$(B.2)

where |titj| is the time between epochs in the series, and the hyperparameters, a and c, were allowed to be free.

The fitting procedure was the same as the one described in Sec. 4.2. In Table B.1 we present the median and 1σ uncertainties of the fitted parameters. With this procedure we find a rotation period of Prot=42.15.8+3.5d${P_{{\rm{rot}}}} = 42.1_{ - 5.8}^{ + 3.5}\,d$, consistent with the rotation period found by fitting the ASAS-SN V-band data.

Table B.1

Prior functions and fitted parameters values for the TJO R-band photometry of Wolf 327.

Appendix C Ground-based transit light curves

In this section we describe our initial fitting process and selection for the ground-based light curves considered for the joint fit of the data. Figure C.1 shows all the ground-based transit observations available in ExoFOP and a time series taken with MuSCAT2, along with a transit model for each individual light curve (black line). The transit model was computed using normal priors in the central time of the transit and orbital period of the candidate, the values for these parameters were taken from the results of a test joint fit of the TESS and CARMENES datasets. To model some of the systematics present in the data we use a transit baseline taking into account the position of the star, the FWHM, and airmass during the observations:

=f0+cXXp+cYYp+cfFw+caAir,${\cal B} = {f_0} + {c_X}{X_p} + {c_Y}{Y_p} + {c_f}{F_w} + {c_a}{A_{ir}},$(C.1)

where the free parameters are f0 as a flux offset, the coefficients cX and cY model the flux variations caused by changes in the star’s position on the detector (Xp, Yp), cf is the coefficient associated with the changes of the PSF’s FWHM (Fw) during the observations, and ca is the free parameter used to model flux changes depending on the airmass (Air). For the individual transit fits we did not use GPs to model the red noise component present in each of the time series.

thumbnail Fig. C.1

Ground-based transit light curves of Wolf 327b. The best fitted model (black line) includes systematic effects but no red noise component (i.e., without GPs).

From Figure C.1 we can see that the observations that cover a full transit event are from 2022 November 19, 2022 December 2 and 12, 2023 January 19, and 2023 May 25. These full transit observations were reanalyzed with AstrolmageJ using some transit parameter values as priors taken from a TESS and CARMENES joint fit. Then systematic instrumental effects were removed from the data, which were used for the final joint fit.

Appendix D Joint fit results

During our data fitting test phase we compared the results of a joint fit performed on TESS photometric data and CARMENES RV measurements but using different approaches to deal with the RVs. In one test we used GPs to model the red noise of both datasets (TESS and CARMENES), for a second model we used GPs for TESS and no GPs for CARMENES, and for the third model we used GPs only for TESS data and no GPs for CARMENES but assuming different zero-point offsets for each observed night. The zero-point offset method assumes that each night has a different radial velocity zero-point offset, meaning that we added five free γ parameters corresponding to each night that the target was observed. Figure D.1 shows the posterior distribution of the radial velocity semi-amplitude K obtained from all the models we tested, from the distributions we can see that either using or not using GPs, or the zero-point offset method, we get the same posterior distribution for K; meaning the that determination of the mass of the planet is unaffected by the use of GPs to model the systematics present in the RV measurements.

thumbnail Fig. D.1

Posterior distribution of the radial velocity semi-amplitude (K) for a joint fit using only TESS and CARMENES RV observations. We show the results for a fit using GPs to model the red noise of both datasets (light yellow), GPs for TESS data and no GPs for CARMENES measurements (light blue), and a fit of TESS photometry including GPs and RVs without GPs but using a different zero-point (ZP) offset for each night (light purple).

thumbnail Fig. D.2

Correlation plot for the fitted transit and orbital parameter for Wolf 327b. The fitted orbital parameters of the eclipsing binary TIC 4918919 and the free parameters used to model the systematic effects were intentionally left out for easy viewing. The blue lines mark the median values of the distribution. For plotting purposes the distributions for the central time of the transit and orbital period were offset by T1 = 2459252 d.

Figure D.2 presents the posterior distributions for the fitted orbital parameters of Wolf 327b using all the data available (TESS, ground-based transit observations, and RV measurements from CARMENES). The fitted parameters for the eclipsing binary TIC 4918919 and Gaussian processes components were left out intentionally from the plot for easy viewing.

References

  1. Adams, E. R., Jackson, B., Johnson, S., et al. 2021, Planet. Sci. J., 2, 152 [NASA ADS] [CrossRef] [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. Ambikasaran, S., Foreman-Mackey, D., Greengard, L., Hogg, D. W., & O’Neil, M. 2015, IEEE Trans. Pattern Anal. Mach. Intell., 38, 252 [Google Scholar]
  4. Angus, R., Aigrain, S., Foreman-Mackey, D., & McQuillan, A. 2015, MNRAS, 450, 1787 [Google Scholar]
  5. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147 [Google Scholar]
  6. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bauer, F. F., Zechmeister, M., & Reiners, A. 2015, A&A, 581, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bauer, F. F., Zechmeister, M., Kaminski, A., et al. 2020, A&A, 640, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Baumeister, P., & Tosi, N. 2023, A&A, 676, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Baumeister, P., Padovan, S., Tosi, N., et al. 2020, ApJ, 889, 42 [Google Scholar]
  11. Brandeker, A., Alibert, Y., Bourrier, V., et al. 2021, Is it raining lava in the evening on 55 Cancri e?, JWST Proposal. Cycle 1, #2084 [Google Scholar]
  12. Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031 [Google Scholar]
  13. Caballero, J. A., Cortés-Contreras, M., Alonso-Floriano, F. J., et al. 2016a, in 19th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (CS19), 148 [Google Scholar]
  14. Caballero, J. A., Guàrdia, J., López del Fresno, M., et al. 2016b, SPIE Conf. Ser., 9910, 99100E [Google Scholar]
  15. Ciardi, D. R., Beichman, C. A., Horch, E. P., & Howell, S. B. 2015, ApJ, 805, 16 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cifuentes, C., Caballero, J. A., Cortés-Contreras, M., et al. 2020, A&A, 642, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Collins, K. 2019, in Am. Astron. Soc. Meeting Abstracts #233, 140.05 [Google Scholar]
  18. Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77 [Google Scholar]
  19. Colome, J., & Ribas, I. 2006, IAU Special Session, 6, 11 [Google Scholar]
  20. Colomé, J., Casteels, K., Ribas, I., & Francisco, X. 2010, SPIE Conf. Ser., 7740, 77403K [Google Scholar]
  21. Cortés Contreras, M. 2017, PhD Thesis, Complutense University of Madrid, Spain [Google Scholar]
  22. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog: II/328 [Google Scholar]
  23. Dai, F., Masuda, K., Winn, J. N., & Zeng, L. 2019, ApJ, 883, 79 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dawson, R. I., & Fabrycky, D. C. 2010, ApJ, 722, 937 [Google Scholar]
  25. Demory, B. O., Gillon, M., Deming, D., et al. 2011, A&A, 533, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Demory, B.-O., Gillon, M., Madhusudhan, N., & Queloz, D. 2016, MNRAS, 455, 2018 [NASA ADS] [CrossRef] [Google Scholar]
  27. Demory, B. O., Sulis, S., Meier Valdés, E., et al. 2023, A&A, 669, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Engle, S. G., & Guinan, E. F. 2018, RNAAS, 2, 34 [NASA ADS] [Google Scholar]
  29. Essack, Z., Seager, S., & Pajusalu, M. 2020, ApJ, 898, 160 [CrossRef] [Google Scholar]
  30. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  31. Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, AJ, 154, 220 [Google Scholar]
  32. Fulton, B. J., Petigura, E. A., Blunt, S., & Sinukoff, E. 2018, PASP, 130, 044504 [Google Scholar]
  33. Gagné, J., Mamajek, E. E., Malo, L., et al. 2018, ApJ, 856, 23 [Google Scholar]
  34. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gibson, N. P., Aigrain, S., Roberts, S., et al. 2012, MNRAS, 419, 2683 [NASA ADS] [CrossRef] [Google Scholar]
  36. Goffo, E., Gandolfi, D., Egger, J. A., et al. 2023, ApJ, 955, L3 [NASA ADS] [CrossRef] [Google Scholar]
  37. Greene, T. P., Bell, T. J., Ducrot, E., et al. 2023, Nature, 618, 39 [NASA ADS] [CrossRef] [Google Scholar]
  38. Heng, K. 2023, ApJ, 956, L20 [NASA ADS] [CrossRef] [Google Scholar]
  39. Hu, R., Brandeker, A., Damiano, M., et al. 2021, Determining the Atmospheric Composition of the Super-Earth 55 Cancri e, JWST Proposal. Cycle 1, #1952 [Google Scholar]
  40. Ih, J., Kempton, E. M. R., Whittaker, E. A., & Lessard, M. 2023, ApJ, 952, L4 [CrossRef] [Google Scholar]
  41. Ingersoll, A. P. 1969, J. Atmos. Sci., 26, 1191 [NASA ADS] [CrossRef] [Google Scholar]
  42. Ito, Y., Ikoma, M., Kawahara, H., et al. 2015, ApJ, 801, 144 [NASA ADS] [CrossRef] [Google Scholar]
  43. Jeffers, S. V., Schöfer, P., Lamert, A., et al. 2018, A&A, 614, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Jeffers, S. V., Dreizler, S., Barnes, J. R., et al. 2020, Science, 368, 1477 [Google Scholar]
  45. Jenkins, J. M. 2002, ApJ, 575, 493 [Google Scholar]
  46. Jenkins, J. M., Chandrasekaran, H., McCauliff, S. D., et al. 2010, SPIE Conf. Ser., 7740, 77400D [Google Scholar]
  47. Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E [Google Scholar]
  48. Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search, Kepler Science Document KSCI-19081-003 [Google Scholar]
  49. Jensen, E. 2013, Astrophysics Source Code Library [record ascl:1306.007] [Google Scholar]
  50. Kasting, J. F. 1988, Icarus, 74, 472 [NASA ADS] [CrossRef] [Google Scholar]
  51. Katz, R. F., Spiegelman, M., & Langmuir, C. H. 2003, Geochem. Geophys. Geosyst., 4, 1073 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kempton, E. M. R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401 [Google Scholar]
  53. Kendrew, S., Scheithauer, S., Bouchet, P., et al. 2015, PASP, 127, 623 [NASA ADS] [CrossRef] [Google Scholar]
  54. Kipping, D. M. 2013, MNRAS, 435, 2152 [Google Scholar]
  55. Kopparapu, R. K., Ramirez, R., Kasting, J. F., et al. 2013, ApJ, 765, 131 [NASA ADS] [CrossRef] [Google Scholar]
  56. Kopparapu, R. K., Ramirez, R. M., SchottelKotte, J., et al. 2014, ApJ, 787, L29 [Google Scholar]
  57. Lam, K. W. F., Csizmadia, S., Astudillo-Defru, N., et al. 2021, Science, 374, 1271 [NASA ADS] [CrossRef] [Google Scholar]
  58. Lee, E. J., & Chiang, E. 2017, ApJ, 842, 40 [Google Scholar]
  59. Lépine, S., & Shara, M. M. 2005, AJ, 129, 1483 [Google Scholar]
  60. Li, J., Tenenbaum, P., Twicken, J. D., et al. 2019, PASP, 131, 024506 [NASA ADS] [CrossRef] [Google Scholar]
  61. Lincowski, A. P., Meadows, V. S., Zieba, S., et al. 2023, ApJ, 955, L7 [NASA ADS] [CrossRef] [Google Scholar]
  62. Lodieu, N., Smart, R. L., Pérez-Garrido, A., & Silvotti, R. 2019, A&A, 623, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Luger, R., & Barnes, R. 2015, Astrobiology, 15, 119 [Google Scholar]
  64. Lundkvist, M. S., Kjeldsen, H., Albrecht, S., et al. 2016, Nat. Commun., 7, 11201 [Google Scholar]
  65. Lustig-Yaeger, J., Fu, G., May, E. M., et al. 2023, Nat. Astron., 7, 1317 [NASA ADS] [CrossRef] [Google Scholar]
  66. MacKenzie, J., Grenfell, J. L., Baumeister, P., et al. 2023, A&A, 671, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Marfil, E., Tabernero, H. M., Montes, D., et al. 2021, A&A, 656, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. McArthur, B. E., Endl, M., Cochran, W. D., et al. 2004, ApJ, 614, L81 [Google Scholar]
  69. McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, SPIE Conf. Ser., 10707, 107070K [Google Scholar]
  70. Meier Valdés, E. A., Morris, B. M., Wells, R. D., Schanche, N., & Demory, B. O. 2022, A&A, 663, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Meier Valdés, E. A., Morris, B. M., Demory, B. O., et al. 2023, A&A, 677, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Melrose, D. B., & Dulk, G. A. 1982, ApJ, 259, 844 [Google Scholar]
  73. Montes, D., López-Santiago, J., Gálvez, M. C., et al. 2001, MNRAS, 328, 45 [NASA ADS] [CrossRef] [Google Scholar]
  74. Moran, S. E., Stevenson, K. B., Sing, D. K., et al. 2023, ApJ, 948, L11 [NASA ADS] [CrossRef] [Google Scholar]
  75. Morello, G., Parviainen, H., Murgas, F., et al. 2023, A&A, 673, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Morris, R. L., Twicken, J. D., Smith, J. C., et al. 2020, Kepler Data Processing Handbook: Photometric Analysis, Kepler Science Document KSCI-19081-003 [Google Scholar]
  77. Mousis, O., Deleuil, M., Aguichine, A., et al. 2020, ApJ, 896, L22 [NASA ADS] [CrossRef] [Google Scholar]
  78. Mugrauer, M., & Michel, K.-U. 2020, Astron. Nachr., 341, 996 [NASA ADS] [CrossRef] [Google Scholar]
  79. Mugrauer, M., & Michel, K.-U. 2021, Astron. Nachr., 342, 840 [NASA ADS] [CrossRef] [Google Scholar]
  80. Mugrauer, M., Zander, J., & Michel, K.-U. 2022, Astron. Nachr., 343, e24017 [NASA ADS] [Google Scholar]
  81. Mugrauer, M., Rück, J., & Michel, K. U. 2023, Astron. Nachr., 344, e20230055 [Google Scholar]
  82. Murgas, F., Nowak, G., Masseron, T., et al. 2022, A&A, 668, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Murgas, F., Castro-González, A., Pallé, E., et al. 2023, A&A, 677, A182 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Narita, N., Fukui, A., Kusakabe, N., et al. 2019, J. Astron. Telesc. Instrum. Syst., 5, 015001 [NASA ADS] [Google Scholar]
  85. Narita, N., Fukui, A., Yamamuro, T., et al. 2020, SPIE Conf. Ser., 11447, 114475K [NASA ADS] [Google Scholar]
  86. Nguyen, T. G., Cowan, N. B., Banerjee, A., & Moores, J. E. 2020, MNRAS, 499, 4605 [NASA ADS] [CrossRef] [Google Scholar]
  87. Nguyen, T. G., Cowan, N. B., Pierrehumbert, R. T., Lupu, R. E., & Moores, J. E. 2022, MNRAS, 513, 6125 [NASA ADS] [CrossRef] [Google Scholar]
  88. Nicholson, B. A., & Aigrain, S. 2022, MNRAS, 515, 5251 [NASA ADS] [CrossRef] [Google Scholar]
  89. Owen, J. E., & Wu, Y. 2013, ApJ, 775, 1 [CrossRef] [Google Scholar]
  90. Parviainen, H. 2015, MNRAS, 450, 3233 [Google Scholar]
  91. Parviainen, H., & Aigrain, S. 2015, MNRAS, 453, 3821 [Google Scholar]
  92. Parviainen, H., Tingley, B., Deeg, H. J., et al. 2019, A&A, 630, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Pérez-Torres, M., Gómez, J. F., Ortiz, J. L., et al. 2021, A&A, 645, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Petrovich, C., Deibert, E., & Wu, Y. 2019, AJ, 157, 180 [NASA ADS] [CrossRef] [Google Scholar]
  95. Price, E. M., & Rogers, L. A. 2020, ApJ, 894, 8 [NASA ADS] [CrossRef] [Google Scholar]
  96. Quirrenbach, A., Amado, P. J., Caballero, J. A., et al. 2014, Proc. SPIE, 9147, 91471F [Google Scholar]
  97. Quirrenbach, A., Amado, P. J., Ribas, I., et al. 2018, SPIE Conf. Ser., 10702, 107020W [Google Scholar]
  98. Rappaport, S., Sanchis-Ojeda, R., Rogers, L. A., Levine, A., & Winn, J. N. 2013, ApJ, 773, L15 [NASA ADS] [CrossRef] [Google Scholar]
  99. Rasmussen, C., & Williams, C. 2010, The MIT Press, 122, 935 [Google Scholar]
  100. Reiners, A., Shulyak, D., Käpylä, P. J., et al. 2022, A&A, 662, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Reylé, C., Jardine, K., Fouqué, P., et al. 2021, A&A, 650, A201 [Google Scholar]
  102. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, SPIE Conf. Ser., 9143, 914320 [Google Scholar]
  103. Sahu, K. C., Casertano, S., Bond, H. E., et al. 2006, Nature, 443, 534 [Google Scholar]
  104. Sanchis-Ojeda, R., Rappaport, S., Winn, J. N., et al. 2013, ApJ, 774, 54 [Google Scholar]
  105. Sanchis-Ojeda, R., Rappaport, S., Winn, J. N., et al. 2014, ApJ, 787, 47 [Google Scholar]
  106. Santerne, A., Brugger, B., Armstrong, D. J., et al. 2018, Nat. Astron., 2, 393 [Google Scholar]
  107. Schiavon, R. P., Barbuy, B., & Singh, P. D. 1997, ApJ, 484, 499 [NASA ADS] [CrossRef] [Google Scholar]
  108. Schlecker, M., Apai, D., Lichtenberg, T., et al. 2023, https://doi.org/10.5281/zenodo.8251078 [Google Scholar]
  109. Schöfer, P. 2021, PhD Thesis, Georg August University of Gottingen, Germany [Google Scholar]
  110. Schöfer, P., Jeffers, S. V., Reiners, A., et al. 2019, A&A, 623, A44 [Google Scholar]
  111. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  112. Schweitzer, A., Passegger, V. M., Cifuentes, C., et al. 2019, A&A, 625, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  113. Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48 [Google Scholar]
  114. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  115. Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 [Google Scholar]
  116. Southworth, J. 2011, MNRAS, 417, 2166 [Google Scholar]
  117. Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 102 [Google Scholar]
  118. Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138 [Google Scholar]
  119. Strakhov, I. A., Safonov, B. S., & Cheryasov, D. V. 2023, Astrophys. Bull., 78, 234 [NASA ADS] [CrossRef] [Google Scholar]
  120. Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985 [Google Scholar]
  121. Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100 [Google Scholar]
  122. Tabernero, H. M., Marfil, E., Montes, D., & González Hernández, J. I. 2022, A&A, 657, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Tal-Or, L., Zechmeister, M., Reiners, A., et al. 2018, A&A, 614, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Tamburo, P., Mandell, A., Deming, D., & Garhart, E. 2018, AJ, 155, 221 [NASA ADS] [CrossRef] [Google Scholar]
  125. Trifonov, T., Kürster, M., Zechmeister, M., et al. 2018, A&A, 609, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  126. Turbet, M., Ehrenreich, D., Lovis, C., Bolmont, E., & Fauchez, T. 2019, A&A, 628, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Turbet, M., Bolmont, E., Ehrenreich, D., et al. 2020, A&A, 638, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  128. Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502 [Google Scholar]
  129. Wing, R. F., & Ford, W. Kent, J. 1969, PASP, 81, 527 [CrossRef] [Google Scholar]
  130. Winn, J. N., Matthews, J. M., Dawson, R. I., et al. 2011, ApJ, 737, L18 [Google Scholar]
  131. Winn, J. N., Sanchis-Ojeda, R., Rogers, L., et al. 2017, AJ, 154, 60 [NASA ADS] [CrossRef] [Google Scholar]
  132. Winn, J. N., Sanchis-Ojeda, R., & Rappaport, S. 2018, New Astron. Rev., 83, 37 [CrossRef] [Google Scholar]
  133. Wizinowich, P., Acton, D. S., Shelton, C., et al. 2000, PASP, 112, 315 [NASA ADS] [CrossRef] [Google Scholar]
  134. Wolf, M. 1919, Astron. Nachr., 209, 11 [NASA ADS] [CrossRef] [Google Scholar]
  135. Wordsworth, R. D., & Pierrehumbert, R. T. 2013, ApJ, 778, 154 [Google Scholar]
  136. Wordsworth, R., & Pierrehumbert, R. 2014, ApJ, 785, L20 [CrossRef] [Google Scholar]
  137. Wright, N. J., Newton, E. R., Williams, P. K. G., Drake, J. J., & Yadav, R. K. 2018, MNRAS, 479, 2351 [Google Scholar]
  138. Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44 [Google Scholar]
  139. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
  140. Zechmeister, M., Anglada-Escudé, G., & Reiners, A. 2014, A&A, 561, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  141. Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Zeng, L., Sasselov, D. D., & Jacobsen, S. B. 2016, ApJ, 819, 127 [Google Scholar]
  143. Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, PNAS, 116, 9723 [Google Scholar]
  144. Zhang, M., Hu, R., Inglis, J., et al. 2024, ApJ, 961, L44 [NASA ADS] [CrossRef] [Google Scholar]
  145. Zieba, S., Kreidberg, L., Ducrot, E., et al. 2023, Nature, 620, 746 [NASA ADS] [CrossRef] [Google Scholar]
  146. Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19 [Google Scholar]

All Tables

Table 1

Stellar parameters of Wolf 327.

Table 2

Prior functions and fitted parameters values for the ASAS-SN V-band photometry of Wolf 327.

Table 3

Wolf 327 activity indicators computed from a coadded CARMENES VIS and NIR channel spectrum.

Table 4

Prior functions, fitted, and derived parameters of Wolf 327b.

Table A.1

CARMENES VIS channel corrected radial velocities, instrumental drift, and nightly zero-point values computed by SERVAL. The RVs presented here were corrected by instrumental drifts and nightly zero-point offsets.

Table A.2

CARMENES spectroscopic line activity indicators for Wolf 327 computed by SERVAL.

Table B.1

Prior functions and fitted parameters values for the TJO R-band photometry of Wolf 327.

All Figures

thumbnail Fig. 1

Tpfplotter (Aller et al. 2020) target pixel file (TPF) images of Wolf 327 for sectors 21 and 48. The red squares correspond to the TESS aperture used to compute the photometry, the size of the red circles represent the Gaia DR3 (Gaia Collaboration 2023) magnitudes of the stars, and the gray arrows shows the Gaia DR3 proper motion directions of each star inside the field of view. Wolf 327b is positioned at the center (star 1), star 7 (TIC 4918919) is a known eclipsing binary.

In the text
thumbnail Fig. 2

TESS PDCSAP photometry for sectors 21 and 48. The deep transit events occurring with a period of ~3.89 days are from the known eclipsing binary TIC 4918919. The solid black line shows the best-fitting model from the joint data fit described in Sect. 5.2. The vertical lines represent the central time of the transit of the USP Wolf 327b (orange) and the eclipsing times of the EB TIC 4918919 (blue).

In the text
thumbnail Fig. 3

Adaptive optics and speckle imaging of Wolf 327. Left panel: Keck NIR AO imaging and sensitivity curve for Wolf 327 taken with the Kcont filter. The black points represent the 5σ limits and are separated in steps of 1 FWHM; the purple area represents the 1σ azimuthal dispersion of the contrast determinations (see text). The inset image of the primary target shows no additional close-in companions. Right panel: Ic band speckle sensitivity curve of Wolf 327 taken with the SPeckle Polarimeter (SPP) at the Caucasian Observatory. No stellar companions are detected with magnitudes brighter than ΔIc = 3.9 and 6.2 at ρ = 0.″25 and 1.″0 respectively. The inset shows the speckle autocorrelation function.

In the text
thumbnail Fig. 4

Wolf 327 CARMENES radial velocity measurements obtained using SERVAL for the visible (orange circles) and near-infrared (blue squares) channels.

In the text
thumbnail Fig. 5

Long-term photometric archival data of Wolf 327 from ASAS-SN (V band). Left panel: photometric time-series and best-fitted model computed using GP adopting a quasi-periodic kernel (black line). The 1σ uncertainty regions of the fit are shown in light blue. The residuals of the fit are shown in the bottom panel. Right panel: GLS periodogram (Zechmeister & Kürster 2009) for ASAS-SN V-band photometry. The vertical green line shows the fitted rotational period of the star. The horizontal lines represent the false alarm probability (FAP) levels of 10% (red dotted line) and 1% (blue dashed line).

In the text
thumbnail Fig. 6

GLS periodograms of the CARMENES RV measurements taken with the optical channel, RV residuals (i.e., linear trend and planet signal removed), spectral activity indices, and window function. The horizontal lines represent the FAP levels of 10% (purple dotted line) and 1% (blue dash line). The vertical black line marks the orbital period of Wolf 327b (P = 0.5734 day).

In the text
thumbnail Fig. 7

TESS phase-folded light curve of Wolf 327b after subtracting the photometric variations from the time series (top panel). The best-fit transit model is shown in black, the circles are TESS binned observations. Bottom panel: residuals of the fit.

In the text
thumbnail Fig. 8

Ground-based individual transit observations of Wolf 327b used in the joint fit, systematic effects were removed from the light curves. The points show the individual observations, the circles represent binned data points, and the best-fit model is shown by the black line.

In the text
thumbnail Fig. 9

Wolf 327b radial velocity measurements taken with CARMENES. (a) RV time series and best-fitting model (blue line). The best-fitting model was computed using the median values of the posterior distribution of the fitted parameters, (b) Residuals of the fit after subtracting the single-planet Keplerian, the velocity linear trend seen here was included in the RV modeling, (c) Phase-folded RV measurements.

In the text
thumbnail Fig. 10

Top panel: TESS phase-folded light curve of the eclipsing binary TIC 4918919. The photometric variations modeled by the GPs were removed from the time series. The best-fit transit model is shown in black, the circles are TESS binned observations. Bottom panel: residuals of the fit.

In the text
thumbnail Fig. 11

Orbital period versus radius diagram for ultra-short-period planets. Confirmed planets with only radius measurements are shown with green squares, and planets with mass and radius values are shown with circles. The color and size of the circles represent the mass of the planet. The position of Wolf 327b is marked by the red star. Data taken from NASA Exoplanet Archive (https://exoplanetarchive.ipac.caltech.edu/).

In the text
thumbnail Fig. 12

Mass-radius diagram for Wolf 327b (red star) and known transiting planets with mass determinations with a precision better than 30% (parameters taken from the TEPcat database; Southworth 2011). Planets orbiting M-type stars (Teff ≤ 4000 K) are marked with circles, orange indicate planets with periods P > 1 day, and red circles indicate USP planets (P ≤ 1 day) around M dwarfs. The lines in the mass–-radius diagram represent the composition models of Zeng et al. (2016, 2019) for planets with pure iron cores (100% Fe, solid brown line), Earth-like rocky compositions (32.5% Fe plus 67.5% MgSiO3, dashed green line) and pure rock (100% MgSiO3, dash-dotted purple line).

In the text
thumbnail Fig. 13

Predicted thickness of the interior layers of Wolf 327b computed with ExoMDN.

In the text
thumbnail Fig. 14

Predicted mass fraction of the interior layers of Wolf 327b computed with ExoMDN.

In the text
thumbnail Fig. 15

Flux density arising from starplanet interaction as a function of the stellar mass-loss rate (panels a and b) and the magnetic field of the planet (panels c and d). (a) Expected flux density for a closed dipolar geometry. (b) Expected flux density for an open Parker spiral geometry. (c) Expected flux density for a closed dipolar geometry. (d) Expected flux density for an open Parker spiral geometry. For (c) and (d) we assume = 0.1 . Note that the resulting radio emission is unaffected by the magnetic field of the planet, unless the latter is very large (~4 G for closed dipole or even > 10 G for a Parker spiral). The emission expected from Saur/Turnpenney’s model is shown in orange, and the emission expected from Zarka/Lanza’s model is shown in blue. The overlap between both models is shown in brown. The dashed black line represents the assumed detection threshold of 100 μJy.

In the text
thumbnail Fig. B.1

Long-term photometric observations of Wolf 327 taken with the Joan Oró telescope (TJO). Left panel: normalized R-band flux (orange points) versus time. The best-fitted model is shown by the black line, the 1σ uncertainty regions of the fit are shown in light blue. Right panel: GLS periodogram (Zechmeister & Kürster 2009) for TJO R-band photometry. The vertical green line shows the best-fitted rotation period (Prot=42.15.8+3.5 d${P_{{\rm{rot}}}} = 42.1_{ - 5.8}^{ + 3.5}\,{\rm{d}}$). The horizontal lines represent the false alarm probability (FAP) levels of 10% (red dotted line) and 1% (blue dashed line).

In the text
thumbnail Fig. C.1

Ground-based transit light curves of Wolf 327b. The best fitted model (black line) includes systematic effects but no red noise component (i.e., without GPs).

In the text
thumbnail Fig. D.1

Posterior distribution of the radial velocity semi-amplitude (K) for a joint fit using only TESS and CARMENES RV observations. We show the results for a fit using GPs to model the red noise of both datasets (light yellow), GPs for TESS data and no GPs for CARMENES measurements (light blue), and a fit of TESS photometry including GPs and RVs without GPs but using a different zero-point (ZP) offset for each night (light purple).

In the text
thumbnail Fig. D.2

Correlation plot for the fitted transit and orbital parameter for Wolf 327b. The fitted orbital parameters of the eclipsing binary TIC 4918919 and the free parameters used to model the systematic effects were intentionally left out for easy viewing. The blue lines mark the median values of the distribution. For plotting purposes the distributions for the central time of the transit and orbital period were offset by T1 = 2459252 d.

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.