Planetary Transits Research Articles

The GAPS programme at TNG

Context. Atmospheric characterisation plays a key role in the study of exoplanetary systems, giving hints about the current and past conditions of the planets. The information retrieved from the analysis of pivotal lines such as the Hα and He I triplet allow us to constrain the evolutionary path of the planets due to atmospheric photo-evaporation. After focussing for many years on ultra-hot Jupiters, atmospheric characterisation is slowly moving towards smaller and colder planets, which are harder to study due to the difficulties in extracting the planetary signal and which require more precise analysis. Aims. We aim to characterise the atmosphere of TOI-5398 b (P ~ 10.59 days), the outer warm Saturn orbiting a young (~650 Myr) G-type star that also hosts the small inner planet TOI-5398 c (P ~ 4.77 days). Both planets are suitable for atmospheric probing due to the closeness to their host star, which results in strong photo-evaporation processes, especially the larger outer one with an estimated transmission spectroscopy metric of 288 (higher than those of several well-known hot Jupiters). Methods. We investigated the atmosphere of planet b, analysing the data collected during a transit with HARPS-N and GIANO-B high-resolution spectrographs, employing both cross-correlation and single-line analysis to study the presence of atomic species. Incidentally, we recorded the simultaneous transit of planet c, and hence we also focussed on discerning the origin of the signal. We expect planet b to be the cause of the detected signal, since, according to existing evaporation models, it is currently expected to lose more mass than planet c. Results. We detected the presence of Hα and He I triplets, two markers of the photo-evaporation processes predicted for the system, retrieving a height in the atmosphere of 2.33 Rp and 1.65 Rp, respectively. We confirmed these predictions by employing the models computed with the ATES software, which predict a He I absorption arising from planet b comparable with the observed one. Moreover, the ATES models suggested an He/H ratio of 1/99 to match our observations. The investigation of atomic species led to the detection of an Na I doublet via single-line analysis, while the cross-correlation did not return a detection for any of the atomic species investigated.

Detecting and sizing the Earth with PLATO: A feasibility study based on solar data

The PLAnetary Transits and Oscillations of stars (PLATO) mission will observe the same area of the sky continuously for at least two years in an effort to detect transit signals of an Earth-like planet orbiting a solar-like star. We aim to study how short-term solar-like variability caused by oscillations and granulation would affect PLATO's ability to detect and size Earth if PLATO were to observe the Solar System itself. We also compare different approaches to mitigate noise caused by short-term solar-like variability and perform realistic transit fitting of transit signals in PLATO-like light curves. We injected Earth-like transit signals onto real solar data taken by the Helioseismic and Magnetic Imager (HMI) instrument on board the Solar Dynamics Observatory (SDO). We isolated short-term stellar variability in the HMI observations by removing any variability with characteristic timescales longer than five hours using a smooth Savitzky-Golay filter. We then added a noise model for a variety of different stellar magnitudes computed by assuming an observation by all 24 normal cameras. We first compared four different commonly used treatments of correlated noise in the time domain by employing them in a transit fitting scheme. We then tried to recover pairs of transit signals using an algorithm similar to the transit least squares algorithm. Finally, we performed transit fits using realistic priors on planetary and stellar parameters and assessed how accurately the pair of two injected transits was recovered. We find that short-term solar-like variability affects the correct retrieval of Earth-like transit signals in PLATO data. Variability models accounting for variations with typical timescales at the order of one hour are sufficient to mitigate these effects. We find that when the limb-darkening coefficients of the host star are properly constrained, the impact parameter does not negatively affect the detectability of a transit signal or the retrieved transit parameters, except for high values ($b > 0.8$). For bright targets (8.5 - 10.5 mag), the transit signal of an Earth analogue can reliably be detected in PLATO data. For faint targets a detection is still likely, though the results of transit search algorithms have to be verified by transit-fitting algorithms to avoid false positive detections being flagged. For bright targets (V-mag leq 9.5), the radius of an Earth-like planet orbiting a solar-like star can be correctly determined at a precision of 3<!PCT!> or less, assuming that at least two transit events are observed and the characteristics of the host star are well understood.

Probing the magnetic fields of starspots with transit mapping

Starspots, regions of strong magnetic fields, serve as indicators of stellar activity and the dynamo mechanism at play in the interior of stars. The magnetic fields of main-sequence stars play a crucial role in driving stellar activity. An effective approach to better understanding stellar magnetic fields and activity lies in the detailed characterisation of starspot properties. We propose a new method for estimating the magnetic fields of starspots that employs modelling techniques of planetary transit mapping, which provides estimates of the size, intensity, and location of spots on the stellar photosphere. A starspot's maximum magnetic field was calculated using the linear relationship with the spot flux deficit, $ spot $ (the spot's brightness times its area) and the well-characterised relation for sunspots determined in this work, $B_ spot (G). Applying this relationship to previously mapped spots on the photospheres of 14 FGK and M stars yields spot maximum magnetic fields ranging from 2700 G to 4600 G, with an overall average of $3900 400$ G. We looked for correlations between starspot magnetic fields and stellar properties.\ We did not find any correlation between a spot's mean extreme magnetic field and effective temperature, nor the differential shear. However, a weak anti-correlation is seen between the spots' magnetic field and stellar age as well as between the magnetic field and the rotation period. When compared with previous results of small-scale magnetic field measurements, the B values obtained here are basically constant and near the saturation limit found for rapid rotators. This implies that it is not the intensity of the magnetic field of starspots that decreases with age but rather the filling factor. This result offers a unique window into the magnetic dynamo of stars.

TESS Investigation—Demographics of Young Exoplanets (TI-DYE). II. A Second Giant Planet in the 17 Myr System HIP 67522

The youngest (<50 Myr) planets are vital to understand planet formation and early evolution. The 17 Myr system HIP 67522 is already known to host a giant (≃10R ⊕) planet on a tight orbit. In their discovery paper, Rizzuto et al. reported a tentative single-transit detection of an additional planet in the system using TESS. Here, we report the discovery of HIP 67522c, a 7.9 R ⊕ planet that matches with that single-transit event. We confirm the signal with ground-based multiwavelength photometry from Sinistro and MuSCAT4. At a period of 14.33 days, planet c is close to a 2:1 mean-motion resonance with b (6.96 days or 2.06:1). The light curve shows distortions during many of the transits, which are consistent with spot-crossing events and/or flares. Fewer stellar activity events are seen in the transits of planet b, suggesting that planet c is crossing a more active latitude. Such distortions, combined with systematics in the TESS light-curve extraction, likely explain why planet c was previously missed.

Measuring stellar surface rotation and activity with the PLATO mission

The Planetary Transits and Oscillations of stars mission (PLATO) will allow us to measure surface rotation and monitor photometric activity of tens of thousands of main sequence solar-type and subgiant stars. This paper is the first of a series dedicated to the preparation of the analysis of stellar surface rotation and photospheric activity with the near-future PLATO data. We describe in this work the strategy that will be implemented in the PLATO pipeline to measure stellar surface rotation, photometric activity, and long-term modulations. The algorithms are applied on both noise-free and noisy simulations of solar-type stars, which include activity cycles, latitudinal differential rotation, and spot evolution. PLATO simulated systematics are included in the noisy light curves. We show that surface rotation periods can be recovered with confidence for most of the stars with only six months of observations and that the recovery rate of the analysis significantly improves as additional observations are collected. This means that the first PLATO data release will already provide a substantial set of measurements for this quantity, with a significant refinement on their quality as the instrument obtains longer light curves. Measuring the Schwabe-like magnetic activity cycle during the mission will require that the same field be observed over a significant timescale (more than four years). Nevertheless, PLATO will provide a vast and robust sample of solar-type stars with constraints on the activity-cycle length. Such a sample is lacking from previous missions dedicated to space photometry.

Differential Rotation of CoRoT Stars and a Kepler Binary Star from Starspot Transit Mapping

Abstract While direct magnetic field measurements are rare for slowly rotating stars, spots observed during planetary transits provide a potential indicator of magnetic activity on stellar surfaces. Moreover, the rotation of the stellar surface can be probed by monitoring the spots’ position with time in subsequent transits. This study investigates the dynamic interplay of rotational shear and stellar rotation rate in six stars of spectral types F to M, all hosting exoplanets observed by the CoRoT space mission, except for one binary star from Kepler. The analysis, facilitated by the ECLIPSE code, unveils the physical properties of stellar spots, including radius, intensity, temperature, and position. The five CoRoT stars exhibit spot characteristics consistent with those observed in solar type stars. The determination of a spot longitude during different transits allows for the inference of the star’s differential rotation profile, revealing a decreasing trend of rotational shear with the mean stellar rotation period, given by Δ Ω ∝ P ¯ − 0.9 . This implies that, at least for slow rotators (mean rotation period &gt;5 days), as stars age, their differential rotation decreases. Additionally, the six stars analyzed here seem to fall into two categories, according to their spot flux deficit: those with ΔF spot &gt; 0.005 and those with ΔF spot &lt; 0.002.

Occultation records in the Royal Frankish Annals for A.D. 807: Knowledge transfer from Arabia to Frankia?

The Royal Frankish Annals (RFA) report under A.D. 807 three lunar and one solar eclipse (complete and dated correctly), a lunar occultation of Jupiter, and “a small spot” on the Sun interpreted as Mercury transit. This transmission is embedded between texts on the landing of diplomatic delegates from Baghdad and Jerusalem in Venice and on their meeting with Charlemagne. Technology transfer is explicitly mentioned (a mechanical water clock). The interpretation of a spot on the Sun as planetary transit is attested in Arabic reports since A.D. 840. The Jupiter occultation by the Moon is dated “Luna 17,” correctly counted in Muslim calendar practice since the first crescent detection. Eclipse predictions and records by astronomers from the cAbbasid caliphate using Ptolemy’s Almagest are extant from the early ninth century. There are time-spans of 6 months each between the three reported lunar eclipses (806 Sep 1/2, 807 Feb 25/26, Aug 21/22), a period of possible repetition given in the Almagest, but not known to the Carolingians. The 6-month rule for lunar eclipses might be preserved in the Carolingian Seven Book Computus. We present multiple evidence that visitors from Baghdad brought new insight and fostered the Carolingian progress in astronomy.

Assessing the residual biases in high-resolution transit absorption spectra correction

Context. In recent years it has become common practice to divide observed transit absorption spectra by synthetic absorption spectra computed for the transit of an atmosphere-less planet. This action supposedly corrects the observed absorption spectrum, leaving the sole atmospheric absorption signature free from the biases induced by stellar rotation and centre-to-limb variations. Aims. We aim to show that while this practice is beneficial, it does not completely correct the absorption spectrum from the stellar distortions and that some residual biases remain, leaving a possibly altered atmospheric signature. Methods. By reducing the problem to its most basic form, we show that dividing the observed absorption spectrum by a synthetic absorption spectrum of the planet does not isolate the pure atmospheric absorption signature. We also used simulated synthetic transit observations to assess the magnitude of these residual biases in typical transit observations. Results. We show that dividing the observed absorption spectrum by the planetary absorption spectrum results in an atmospheric signature modulated by the ratio of the flux behind the atmosphere and the flux behind the planet. Depending on the non-homogeneity of the stellar spectrum, this leads to distorted atmospheric signatures. Eventually, directly analysing these biased signatures will lead to wrong estimates of planetary atmosphere properties.

NotPlaNET: Removing False Positives from Planet Hunters TESS with Machine Learning

Differentiating between real transit events and false-positive signals in photometric time-series data is a bottleneck in the identification of transiting exoplanets, particularly long-period planets. This differentiation typically requires visual inspection of a large number of transit-like signals to rule out instrumental and astrophysical false positives that mimic planetary transit signals. We build a one-dimensional convolutional neural network (CNN) to separate eclipsing binaries and other false positives from potential planet candidates, reducing the number of light curves that require human vetting. Our CNN is trained using the TESS light curves that were identified by Planet Hunters citizen scientists as likely containing a transit. We also include the background flux and centroid information. The light curves are visually inspected and labeled by project scientists and are minimally preprocessed, with only normalization and data augmentation taking place before training. The median percentage of contaminants flagged across the test sectors is 18% with a maximum of 37% and a minimum of 10%. Our model keeps 100% of the planets for 16 of the 18 test sectors, while incorrectly flagging one planet candidate (0.3%) for one sector and two (0.6%) for the remaining sector. Our method shows potential to reduce the number of light curves requiring manual vetting by up to a third with minimal misclassification of planet candidates.

Study of stellar properties under TESS satellite observations

This study employs data from the TESS satellite to investigate the properties of stars. TESS, the Transiting Exoplanet Survey Satellite, offers insights into star activity, periodic variations, and planetary transits by monitoring changes in stellar brightness. Our analysis of TESS data explores critical star attributes, including mass, luminosity, activity levels, and dynamic changes. By examining the relationship between these stellar properties and the presence of planets, we aim to uncover the intricate interplay between celestial bodies and unveil the universe's mysteries. This research encompasses data collection methods, processing techniques, and data analysis approaches, addressing star classifications, parameters, and periodic stellar activity. We also delve into flare energy, which reveals the energy released during stellar flares, and photoperiodic analysis, a tool for understanding periodic light changes in celestial bodies. The study concludes by highlighting the significance of newly discovered planets and the diversity of stars observed by TESS. It underscores the importance of comprehending the connections between stars and planets, shaping our understanding of celestial habitability. In summary, this study enriches our knowledge of star and planet properties, contributing to a broader comprehension of the universe. Through our analysis, we seek to illuminate the complex relationships between celestial entities and unveil the universe's enigmas.

Interpretation of the Transit Light Curve in the Presence of Principal Main Minimum with Allowance for the Eccentricity of the Transit (Planet) Orbit

Interpretation of the Transit Light Curve in the Presence of Principal Main Minimum with Allowance for the Eccentricity of the Transit (Planet) Orbit

Latitudinal Asymmetry in the Dayside Atmosphere of WASP-43b

We present two-dimensional near-infrared temperature maps of the canonical hot Jupiter WASP-43b using a phase-curve observation with JWST NIRSpec/G395H. From the white-light planetary transit, we improve constraints on the planet’s orbital parameters and measure a planet-to-star radius ratio of 0.15883−0.00053+0.00056 . Using the white-light phase curve, we measure a longitude of maximum brightness of 6.9−0.°5+0.°5 east of the substellar point and a phase-curve offset of 10.0−0.°8+0.°8 . We also find a ≈4σ detection of a latitudinal hotspot offset of −13.4−1.°7+3.°2 , the first significant detection of a nonequatorial hotspot in an exoplanet atmosphere. We show that this detection is robust to variations within planetary parameter uncertainties, but only if the transit is used to improve constraints, showing the importance of transit observations to eclipse mapping. Maps retrieved from the NRS1 and NRS2 detectors are similar, with hotspot locations consistent between the two detectors at the 1σ level. Our JWST data show brighter (hotter) nightsides and a dimmer (colder) dayside at the shorter wavelengths relative to fits to Spitzer 3.6 and 4.5 μm phase curves. Through comparison between our phase curves and a set of general circulation models, we find evidence for clouds on the nightside and atmospheric drag or high metallicity reducing the eastward hotspot offset.

New ephemerides and detection of transit-timing variations in the K2-138 system using high-precision CHEOPS photometry

Multi-planet systems are a perfect laboratory for constraining planetary formation models. A few of these systems present planets that come very close to mean motion resonance, potentially leading to significant transit-timing variations (TTVs) due to their gravitational interactions. Of these systems represents a excellent laboratory for studying the dynamics of its six small planets (with radii ranging between $ as the five innermost planets are in a near 3:2 resonant chain. In this work, we aim to constrain the orbital properties of the six planets in the system by monitoring their transits with CHaracterising ExOPlanets Satellite ( We also seek to use this new data to lead a TTV study on this system. We obtained twelve light curves of the system with transits of planets $d$, $e$, $f,$ and $g$. With these data, we were able to update the ephemerides of the transits for these planets and search for timing transit variations. With our measurements, we reduced the uncertainties in the orbital periods of the studied planets, typically by an order of magnitude. This allowed us to correct for large deviations, on the order of hours, in the transit times predicted by previous studies. This is key to enabling future reliable observations of the planetary transits in the system. We also highlight the presence of potential TTVs ranging from 10 minutes to as many as 60 minutes for planet $d$.

Beyond Circular Eclipsers (BeyonCE) light curve modelling

Context. Time series photometry offers astronomers the tools to study time-dependent astrophysical phenomena, from stellar activity to fast radio bursts and exoplanet transits. Transit events, in particular, are focussed primarily on planetary transits and eclipsing binaries with eclipse geometries that can be parameterised with a few variables. However, more complex light curves caused by the substructure within the transiting object would require a customised analysis code. Aims. We present Beyond Circular Eclipsers (BeyonCE), which reduces the parameter space encompassed by the transit of circumsec-ondary disc (CSD) systems with azimuthally symmetric, non-uniform optical-depth profiles. By rejecting disc geometries that are not able to reproduce the measured gradients within their light curves, we can constrain the size and orientation of discs with a complex sub-structure. Methods. We mapped out all the possible geometries of a disc and calculated the gradients for rings crossing the star. We then rejected those configurations where the measured gradient of the light curve is greater than the theoretical gradient from the given disc orientation. Results. We present the fitting code BeyonCE and demonstrate its effectiveness in considerably reducing the parameter space of discs that contain an azimuthally symmetric structure. We used the code to analyse the light curves seen towards J1407 and PDS 110, attributed to CSD transits.

HD 110067 c has an aligned orbit

Planetary systems in mean motion resonances hold a special place among the planetary population. They allow us to study planet formation in great detail as dissipative processes are thought to have played an important role in their existence. Additionally, planetary masses in bright resonant systems can be independently measured via both radial velocities and transit timing variations. In principle, they also allow us to quickly determine the inclination of all planets in the system since, for the system to be stable, they are likely all in coplanar orbits. To describe the full dynamical state of the system, we also need the stellar obliquity, which provides the orbital alignment of a planet with respect to the spin of its host star and can be measured thanks to the Rossiter–McLaughlin effect. It was recently discovered that HD 110067 harbors a system of six sub-Neptunes in resonant chain orbits. We here analyze an ESPRESSO high-resolution spectroscopic time series of HD 110067 during the transit of planet c. We find the orbit of HD 110067 c to be well aligned, with a sky-projected obliquity of λ =6+24-26 deg. This result indicates that the current architecture of the system was reached through convergent migration without any major disruptive events. Finally, we report transit-timing variation in this system as we find a significant offset of 19 ± 4 min in the center of the transit compared to the published ephemeris.

The AstraLux-TESS high spatial resolution imaging survey

Context. Chance-aligned sources or blended companions can cause false positives in planetary transit detections or simply bias the determination of the candidate properties. In the era of high-precision space-based photometers, the need for high spatial resolution images has been demonstrated to be critical for validating and confirming transit signals. This already applied to theKeplermission, is now applicable to the TESS survey, and will be critical for the PLATO mission.Aims. In this paper we present the results of the AstraLux-TESS survey, a catalog of high spatial resolution images obtained with the AstraLux instrument at the Calar Alto observatory (Almería, Spain) in the context of the TESS Follow-up Observing Program.Methods. We used the lucky imaging technique to obtain high spatial resolution images from planet candidate hosts included mostly in two relevant regimes: exoplanet candidates belonging to the level one requirement of the TESS mission (planets with radiiR &lt;4R⊕) and TESS planet candidates around intermediate-mass main-sequence stars.Results. Among the 185 planet host candidate stars observed, we found 13 (7%) to be accompanied by additional sources within a separation of 2.2 arcsec. Among them, six are not associated with sources in theGaiaDR3 catalog, thus contaminating the TESS light curve. Even if no contaminants have been detected, we can provide upper limits and probabilities to the possible existence of field contaminants through the sensitivity limits of our images. Among the isolated hosts, we can discard hazardous companions (bright enough to mimic a planetary transit signals) with an accuracy below 1% for all their planets.Conclusions. The results from this catalog are key to the statistical validation of small planets (prime targets of the TESS mission) and planets around intermediate-mass stars in the main sequence. These two populations of planets are difficult to confirm with the radial velocity technique because of the shallow amplitude of small planets and the high rotational velocities and low number of available spectral lines in the intermediate stellar mass regime. Our results also demonstrate the importance of this type of follow-up observation for future transit missions such as PLATO, even in theGaiaera.

Nuance: Efficient Detection of Planets Transiting Active Stars

The detection of planetary transits in the light curves of active stars, featuring correlated noise in the form of stellar variability, remains a challenge. Depending on the noise characteristics, we show that the traditional technique that consists of detrending a light curve before searching for transits alters their signal-to-noise ratio and hinders our capability to discover exoplanets transiting rapidly rotating active stars. We present nuance, an algorithm to search for transits in light curves while simultaneously accounting for the presence of correlated noise, such as stellar variability and instrumental signals. We assess the performance of nuance on simulated light curves as well as on the Transiting Exoplanet Survey Satellite light curves of 438 rapidly rotating M dwarfs. For each data set, we compare our method to five commonly used detrending techniques followed by a search with the Box-Least-Squares algorithm. Overall, we demonstrate that nuance is the most performant method in 93% of cases, leading to both the highest number of true positives and the lowest number of false-positive detections. Although simultaneously searching for transits while modeling correlated noise is expected to be computationally expensive, we make our algorithm tractable and available as the JAX-powered Python package nuance, allowing its use on distributed environments and GPU devices. Finally, we explore the prospects offered by the nuance formalism and its use to advance our knowledge of planetary systems around active stars, both using space-based surveys and sparse ground-based observations.

High-resolution Spectroscopic Reconnaissance of a Temperate Sub-Neptune

The study of temperate sub-Neptunes is the new frontier in exoplanetary science. A major development in the past year has been the first detection of carbon-bearing molecules in the atmosphere of a temperate sub-Neptune, K2-18 b, a possible Hycean world, with the James Webb Space Telescope (JWST). The JWST is poised to characterize the atmospheres of several other such planets, with important implications for planetary processes in the temperate regime. Meanwhile, ground-based high-resolution spectroscopy has been highly successful in detecting chemical signatures of giant exoplanets, though low-mass planets have remained elusive. In the present work, we report the atmospheric reconnaissance of a temperate sub-Neptune, TOI-732 c, using ground-based high-resolution transmission spectroscopy. The long orbital period and the low systemic velocity result in a low planetary radial velocity during transit, making this system a valuable test bed for high-resolution spectroscopy of temperate sub-Neptunes. We observe high-resolution time-series spectroscopy in the H and K bands during the planetary transit with the IGRINS instrument (R ∼ 45,000) on Gemini-South. Using observations from a single transit, we find marginal evidence (2.2σ) for the presence of methane (CH4) in the atmosphere and no evidence for ammonia (NH3) despite its strong detectability for a cloud-free H2-rich atmosphere. We assess our findings using injection tests with different atmospheric scenarios and find them to be consistent with a high CH4/NH3 ratio and/or the presence of high-altitude clouds. Our results demonstrate the capability of Gemini-S/IGRINS for atmospheric characterization of temperate sub-Neptunes and the complementarity between space- and ground-based facilities in this planetary regime.

Doppler tomography as a tool for characterising exoplanet atmospheres II: an analysis of HD 179949 b

Abstract High-resolution Doppler spectroscopy provides an avenue to study the atmosphere of both transiting and non-transiting planets. This powerful method has also yielded some of the most robust atmospheric detections to date. Currently, high-resolution Doppler spectroscopy detects atmospheric signals by cross-correlating observed data with a model atmospheric spectrum. This technique has been successful in detecting various molecules such as H2O, CO, HCN and TiO, as well as several atomic species. Here we present an alternative method of performing high-resolution Doppler spectroscopy, using a technique known as Doppler tomography. We present an analysis of HD 179949 b using Doppler tomography and provide Doppler tomograms confirming previous detections of CO at 2.3 μm and H2O at both 2.3 μm and 3.5 μmwithin the atmosphere of HD 179949 b, showing significantly lower background noise levels when compared to cross-correlation methods applied to the same data. We also present a novel detection of H2O at 2.1 μm as well as a tentative detection of CO on the night side of the planet at 2.3 μm. This represents the first observational evidence for molecular absorption in the night-side emission spectrum of an exoplanet using Doppler spectroscopy.

Transit spectroscopy of K2-33b with subaru/IRD: Spin-Orbit alignment and tentative atmospheric helium

ABSTRACT Exoplanets in their infancy are ideal targets to probe the formation and evolution history of planetary systems, including the planet migration and atmospheric evolution and dissipation. In this paper, we present spectroscopic observations and analyses of two planetary transits of K2-33b, which is known to be one of the youngest transiting planets (age ≈ 8–11 Myr) around a pre-main-sequence M-type star. Analysing K2-33’s near-infrared spectra obtained by the IRD instrument on Subaru, we investigate the spin-orbit angle and transit-induced excess absorption for K2-33b. We attempt both classical modelling of the Rossiter–McLaughlin (RM) effect and Doppler-shadow analyses for the measurements of the projected stellar obliquity, finding a low angle of $\lambda =-6_{-58}^{+61}$ deg (for RM analysis) and $\lambda =-10_{-24}^{+22}$ deg (for Doppler-shadow analysis). In the modelling of the RM effect, we allow the planet-to-star radius ratio to float freely to take into account the possible smaller radius in the near infrared, but the constraint we obtain ($R_p/R_s=0.037_{-0.017}^{+0.013}$) is inconclusive due to the low radial-velocity precision. Comparison spectra of K2-33 of the 1083 nm triplet of metastable ortho-He I obtained in and out of the 2021 transit reveal excess absorption that could be due to an escaping He-rich atmosphere. Under certain conditions on planet mass and stellar XUV emission, the implied escape rate is sufficient to remove an Earth-mass H/He in ∼1 Gyr, transforming this object from a Neptune to a super-Earth.