---
title: 'WASP-156 b: Aligned Exo-Neptune'
url: https://www.emergentmind.com/topics/wasp-156-b
type: topic
---

# WASP-156 b: Aligned Exo-Neptune

WASP-156 b is a short-period exo-Neptune whose orbital geometry has been reassessed through high-resolution transit spectroscopy and joint photometric–radial-velocity analysis. The principal result is that the planet does not follow a polar orbit; instead, it is consistent with an aligned configuration, with a projected obliquity measured from ESPRESSO data of $λ=-8\pm16^\circ$ [2605.27291]. In the same study, the system parameters were updated, the stellar projected rotation was found to be substantially lower than previously reported, and long-term radial-velocity monitoring was used to exclude Jupiter-mass companions within 5 au. Collectively, the measured aligned and circular orbit, together with the absence of nearby massive companions, are consistent with in situ formation or early disc-driven migration [2605.27291].

## 1. System classification and physical parameters

WASP-156 b is described as a warm “ridge” Neptune with $P\approx3.84$ d, $R\approx6.15\,R_\oplus$, $M\approx42\,M_\oplus$, and $\rho\approx1$ g/cm$^3$ [2605.27291]. In the terminology used for the short-period Neptune population, it resides between the hot-Neptune “desert” ($P\lesssim3$ d) and the longer-period “savanna” [2605.27291]. This placement is relevant because the population of short-period exo-Neptunes is thought to be shaped by an interplay between orbital migration, tidal effects, and photoevaporation.

The stellar parameters derived from ESPRESSO spectra with the PAWS spectral-synthesis pipeline are $T_{\rm eff}=5036\pm104$ K, $\log g=4.46\pm0.08$, and ${\rm [Fe/H]}=+0.32\pm0.15$ [2605.27291]. The same analysis gives a projected stellar rotation of $v\sin i_\star\leq2$ km/s as an upper limit, superseding the literature SOPHIE-based value $v\sin i_\star\simeq3.8\pm0.9$ km/s [2605.27291].

From a joint fit to all TESS and NGTS photometry together with archival out-of-transit radial velocities from CORALIE, SOPHIE, HARPS, and HIRES, the orbital period was updated to $P=3.8361629\pm0.0000006$ d and the transit epoch to $T_0=2459058.6113\pm0.0001$ BJD [2605.27291]. The same fit yielded $R_{\rm p}/R_\star=0.0682\pm0.0007$, impact parameter $b=0.40\pm0.05$, stellar density $\rho_\star=2.25\pm0.14$ g cm$^{-3}$, and radial-velocity semi-amplitude $K=18.6\pm0.8$ m/s [2605.27291]. Eccentricity was fixed to $e=0$, with $e<0.16$ at $3σ$ [2605.27291].

The derived planetary properties are $M_{\rm p}=42\pm2\,M_\oplus$, $R_{\rm p}=6.15\pm0.12\,R_\oplus$, mean density $\rho_{\rm p}=0.99\pm0.07$ g cm$^{-3}$, and semi-major axis $a=0.046\pm0.001$ au [2605.27291].

## 2. Observational basis

The reassessment of the system relied on a multi-instrument campaign centered on transit spectroscopy. Two transits were observed with ESPRESSO on ESO’s VLT in high-resolution mode at $R\approx138\,000$, on 2022-09-02 and 2022-09-29 [2605.27291]. Each night produced 21 spectra, with signal-to-noise ratios of approximately 85 per pixel at 703 nm for the first night and approximately 70 for the second; fibre A tracked the star and fibre B monitored sky, with exposure times of 600–700 s [2605.27291].

A third transit was observed with MAROON-X on Gemini North on 2021-11-04 at $R\approx85\,000$, spanning 5000–9200 Å [2605.27291]. Thirty-eight spectra were obtained, 17 of them in transit, using 400 s exposures; the signal-to-noise ratio increased from approximately 60 to 100 toward egress [2605.27291].

Simultaneous ground-based photometry was obtained with NGTS during the two ESPRESSO transits [2605.27291]. NGTS used 4–5 20 cm telescopes, a custom 520–890 nm filter, and 10 s exposures, achieving approximately 1 ppt precision after airmass detrending [2605.27291].

The ESPRESSO and MAROON-X spectra were reduced to 2D blaze-corrected, sky-subtracted spectra, and disc-integrated radial velocities were extracted both via cross-correlation functions and via template matching with SERVAL to maximize precision [2605.27291]. These high-cadence, high-S/N in-transit radial-velocity time series contain the Rossiter–McLaughlin anomaly, produced because the transiting planet blocks light from different line-of-sight velocity regions across the rotating stellar disk [2605.27291].

## 3. Rossiter–McLaughlin measurement and projected obliquity

Two independent approaches were used to measure the sky-projected obliquity $λ$ from the Rossiter–McLaughlin signal [2605.27291]. The first was a classical RM analysis based on disc-integrated radial velocities, implemented with ironman/rmfit and using the ESPRESSO SERVAL radial velocities only, because they were more precise than the CCF radial velocities [2605.27291].

In that framework, the in-transit velocity anomaly was modeled as

$$
\Delta V_{RM}(t)\;=\;-\,\frac{\iint_{\rm occulted} v_{\rm rot}(x,y)\,I(\mu)\,dA}{\iint_{\,\rm unocculted} I(\mu)\,dA},
$$

with

$$
v_{\rm rot}(x,y)= (v\sin i_\star)\frac{-x\sin \lambda + y\cos \lambda}{R_\star}.
$$

The fit adopted Gaussian priors on orbital parameters from Table 1, uninformative priors on $λ$ and $v\sin i_\star$, specifically $U[-180^\circ,180^\circ]$ and $U[0,2]$ km/s, and a prior on the total intrinsic broadening $\beta=3.6\pm1.9$ km/s [2605.27291]. The joint fit to both ESPRESSO nights yielded $λ=-15^{\circ\,+17^\circ}_{\ \ \ -16^\circ}$ and $v\sin i_\star=0.44\pm0.08$ km/s [2605.27291].

The second approach was the reloaded RM technique, using spatially resolved local CCFs following Cegla et al. 2016 [2605.27291]. For each exposure, the out-of-transit master-out CCF was subtracted from the in-transit CCF, after normalization by the instantaneous photometric flux, to isolate the local CCF from the blocked stellar patch [2605.27291]. Gaussian fits to these local CCFs provided local radial velocities $v_{\rm local}(\phi)$, which were then modeled under a rigid-body rotation law,

$$
v_{\rm model}(\phi) = (v\sin i_\star)\,\frac{x(\phi)\cos \lambda + y(\phi)\sin \lambda}{R_\star}.
$$

An MCMC fit with emcee, using uniform priors $λ\in[-180^\circ,180^\circ]$ and $v\sin i_\star\in[0,100]$ km/s, gave $λ=-8^\circ\pm16^\circ$ and $v\sin i_\star=0.40\pm0.11$ km/s from the ESPRESSO data taken over both nights [2605.27291].

These two analyses are mutually consistent and both indicate an aligned orbit rather than a polar one. The lower measured $v\sin i_\star$ is central to that inference, because the study explicitly identifies a reduction of the $λ$–$v\sin i_\star$ degeneracy when the stellar projected rotation is constrained near 0.4 km/s rather than near 4 km/s [2605.27291].

## 4. Orbital architecture and companion constraints

The updated global fit supports a circular short-period orbit, with eccentricity fixed to zero and a $3σ$ upper bound of $e<0.16$ [2605.27291]. Long-term radial-velocity residuals from the juliet analysis, spanning approximately 7.8 yr and combining CORALIE, SOPHIE, HARPS, and HIRES, show no additional trends or periodicities [2605.27291].

A Monte-Carlo injection–recovery analysis rules out, at 5σ confidence, any companion $\gtrsim1\,M_{\rm Jup}$ interior to approximately 5 au [2605.27291]. Independent high-resolution imaging constraints are also reported: Gemini/Zorro contrast curves at 562/832 nm exclude stellar companions with $\Delta{\rm mag}\lesssim5$–7 from 0.2″–1.2″, corresponding to approximately 25–150 au, while Shane adaptive optics excludes companions of approximately 3–7 mag from approximately 60–500 au [2605.27291]. Gaia DR3 RUWE of approximately 1.1 and the El-Badry et al. wide-binary catalog similarly show no stellar-mass companion out to 1 pc [2605.27291].

These results define the current observational picture of the system’s architecture: a close-in Neptune on an aligned, nearly circular orbit, without evidence for massive nearby perturbers on either planetary or stellar scales. A plausible implication is that mechanisms requiring a present external driver are not supported by the available data.

## 5. Reassessment of the earlier polar-orbit claim

A previous study by Bourrier et al. (2023), based on CARMENES data, reported $λ=105.7^\circ\pm14.4^\circ$ and therefore described the system as having a “polar orbit” [2605.27291]. That analysis adopted $v\sin i_\star=3.8\pm0.9$ km/s from SOPHIE line-profile calibrations [2605.27291].

The later reassessment identifies several reasons for the discrepancy. The CARMENES RM fit is described as suffering from low signal-to-noise, since the observations were obtained with the 3.5 m Calar Alto telescope rather than 8 m class telescopes; from strong telluric contamination and airmass-correlated systematics; and from a restrictive but inaccurate prior on $v\sin i_\star$ based on an empirical calibration outside its valid $T_{\rm eff}$ range [2605.27291]. By contrast, the ESPRESSO and MAROON-X spectra, combined with PAWS spectral modelling, indicate $v\sin i_\star\lesssim2$ km/s, while both the classical and reloaded RM analyses yield values near 0.4 km/s [2605.27291].

The same paper also notes an independent classical RM analysis by Jiang et al. (2026), using the first ESPRESSO night and CCF radial velocities with old system parameters, which found $λ=-5^{\circ\,+33^\circ}_{\ \ \ -35^\circ}$ [2605.27291]. That result is also consistent with alignment, although less precise.

The central methodological point is therefore not merely that a different best-fit obliquity was obtained, but that the inferred obliquity changed once the stellar projected rotation was re-estimated using higher-quality spectroscopy and once the RM signal was analyzed with updated system parameters. This suggests that the earlier polar interpretation was sensitive to prior assumptions and data quality limitations.

## 6. Dynamical interpretation and evolutionary implications

The aligned projected obliquity, the nearly circular orbit, and the lack of nearby massive companions are interpreted as being consistent with either in situ formation or early disc-driven migration within a primordial protoplanetary disc aligned with the stellar spin [2605.27291]. In contrast, high-eccentricity migration scenarios such as planet–planet scattering or Kozai–Lidov forcing would require an external companion capable of exciting $e\gg0$ and would generically produce misaligned orbits [2605.27291]. The observational constraints reported for the system do not support such a configuration.

The study further evaluates the true three-dimensional obliquity using the Fabrycky & Winn (2009) relation,

$$
\cos \psi\;=\;\sin i_\star\,\cos \lambda\,\sin i_p \;+\;\cos i_\star\,\cos i_p.
$$

When the distribution is simulated from $λ$ and random $\cos i_\star$, the resulting median is $\psi\sim35^\circ$, which makes a true polar configuration “extremely unlikely” [2605.27291]. The paper also reports that the stellar tidal-realignment timescale is $\tau_{\rm align}\gtrsim10^4$ Gyr, much longer than the system age of approximately 5–6 Gyr, so tides cannot erase a primordial misalignment [2605.27291]. This strengthens the interpretation that the currently observed alignment is likely primordial rather than the result of later tidal damping.

Photoevaporative evolution is also discussed. The study states that any low-density Neptune migrating into the desert or ridge would likely be eroded, and that WASP-156 b, with $\rho\approx1$ g/cm$^3$, lies near the “density brink,” suggesting that it survived moderate atmospheric loss [2605.27291]. He I 1083 nm and Ly$\alpha$ observations show no detectable escape to date [2605.27291]. This suggests a scenario in which the planet’s present-day density reflects some atmospheric sculpting, but not the level of mass loss required to remove the envelope entirely.

Within the broader short-period Neptune population, the principal significance of the revised result is taxonomic as well as dynamical. WASP-156 b is moved from a tentative cluster of close-in Neptunes in polar orbits to the group of aligned Neptunes [2605.27291]. That reclassification matters because stellar obliquity is used as a diagnostic of formation and migration history; for WASP-156 b, the currently available evidence is most naturally read as favoring formation in, and migration through, a well-aligned disc, with moderate photoevaporation shaping the current bulk density [2605.27291].

Source: https://www.emergentmind.com/topics/wasp-156-b