-
High-resolution atmospheric retrievals of WASP-76b transmission spectroscopy with ESPRESSO: Monitoring limb asymmetries across multiple transits
Authors:
Cathal Maguire,
Neale P. Gibson,
Stevanus K. Nugroho,
Mark Fortune,
Swaetha Ramkumar,
Siddharth Gandhi,
Ernst de Mooij
Abstract:
Direct atmospheric retrievals of exoplanets at high-resolution have recently allowed for a more detailed characterisation of their chemistry and dynamics from the ground. By monitoring the longitudinal distribution of species, as well as the varying vertical temperature structure and dynamics between the limbs of WASP-76b, across multiple transits, we aim to enhance our understanding of the 3D nat…
▽ More
Direct atmospheric retrievals of exoplanets at high-resolution have recently allowed for a more detailed characterisation of their chemistry and dynamics from the ground. By monitoring the longitudinal distribution of species, as well as the varying vertical temperature structure and dynamics between the limbs of WASP-76b, across multiple transits, we aim to enhance our understanding of the 3D nature and chemical/dynamical evolution of such objects over the timescales of months/years. We present retrievals of three VLT/ESPRESSO observations of the ultra-hot Jupiter WASP-76b, including one not yet reported in the literature, from which we constrain the atmospheric abundances, vertical temperature structure, and atmospheric dynamics for both the leading and trailing limbs of the atmosphere separately, via novel rotational broadening kernels. We confirm the presence of VO recently reported in the atmosphere of WASP-76b. We find a uniform longitudinal distribution of Fe and Mg across the limbs of the atmosphere, for each of our transits, consistent with previous works as well as with stellar values. We constrain substellar Na/Fe and Cr/Fe ratios across each of our transits, consistent with previous studies of WASP-76b. Where constrained, V/Fe and VO/Fe ratios were also found to be broadly consistent between the limbs of the atmosphere for each of the transits, as well as with previous studies. However, for two of the transits, both V and VO were unconstrained in the leading limb, suggesting possible depletion due to recombination and condensation. The consistency of our constraints across multiple high-resolution observations, as well as with previous studies using varying modelling/retrieval frameworks and/or instruments, affirms the efficacy of high-resolution ground-based retrievals of exoplanetary atmospheres.
△ Less
Submitted 17 April, 2024; v1 submitted 16 April, 2024;
originally announced April 2024.
-
How do wavelength correlations affect transmission spectra? Application of a new fast and flexible 2D Gaussian process framework to transiting exoplanet spectroscopy
Authors:
Mark Fortune,
Neale P. Gibson,
Daniel Foreman-Mackey,
Thomas M. Evans-Soma,
Cathal Maguire,
Swaetha Ramkumar
Abstract:
The use of Gaussian processes (GPs) is a common approach to account for correlated noise in exoplanet time series, particularly for transmission and emission spectroscopy. This analysis has typically been performed for each wavelength channel separately, with the retrieved uncertainties in the transmission spectrum assumed to be independent. However, the presence of noise correlated in wavelength…
▽ More
The use of Gaussian processes (GPs) is a common approach to account for correlated noise in exoplanet time series, particularly for transmission and emission spectroscopy. This analysis has typically been performed for each wavelength channel separately, with the retrieved uncertainties in the transmission spectrum assumed to be independent. However, the presence of noise correlated in wavelength could cause these uncertainties to be correlated, which could significantly affect the results of atmospheric retrievals. We present a method that uses a GP to model noise correlated in both wavelength and time simultaneously for the full spectroscopic dataset while avoiding the use of a 'common-mode' correction. To make this analysis computationally tractable, we introduce a new fast and flexible GP method that can analyse 2D datasets when the input points lie on a (potentially non-uniform) 2D grid - in our case a time by wavelength grid - and the kernel function has a Kronecker product structure. This simultaneously fits all light curves and enables the retrieval of the covariance matrix of the transmission spectrum. By testing on synthetic datasets, we demonstrate that our new approach can reliably recover atmospheric features contaminated by noise correlated in time and wavelength. In contrast, fitting each spectroscopic light curve separately performed poorly when wavelength-correlated noise was present. It frequently underestimated the uncertainty of the scattering slope and overestimated the uncertainty in the strength of sharp absorption peaks in transmission spectra. Two archival VLT/FORS2 transit observations of WASP-31b were re-analysed, with our method strongly constraining the presence of wavelength-correlated noise in both datasets and recovering significantly different constraints on atmospheric features such as the scattering slope and strength of sodium and potassium features.
△ Less
Submitted 4 March, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
High-resolution emission spectroscopy retrievals of MASCARA-1b with CRIRES+: Strong detections of CO, H$_2$O and Fe emission lines and a C$/$O consistent with solar
Authors:
Swaetha Ramkumar,
Neale P. Gibson,
Stevanus K. Nugroho,
Cathal Maguire,
Mark Fortune
Abstract:
The characterization of exoplanet atmospheres has proven to be successful using high-resolution spectroscopy. Phase curve observations of hot/ultra-hot Jupiters can reveal their compositions and thermal structures, thereby allowing the detection of molecules and atoms in the planetary atmosphere using the cross-correlation technique. We present pre-eclipse observations of the ultra-hot Jupiter, MA…
▽ More
The characterization of exoplanet atmospheres has proven to be successful using high-resolution spectroscopy. Phase curve observations of hot/ultra-hot Jupiters can reveal their compositions and thermal structures, thereby allowing the detection of molecules and atoms in the planetary atmosphere using the cross-correlation technique. We present pre-eclipse observations of the ultra-hot Jupiter, MASCARA-1b, observed with the recently upgraded CRIRES+ high-resolution infrared spectrograph at the VLT. We report a detection of $\rm Fe$ ($\approx$8.3$σ$) in the K-band and confirm previous detections of $\rm CO$ (>15$σ$) and $\rm H_2O$ (>10$σ$) in the day-side atmosphere of MASCARA-1b. Using a Bayesian inference framework, we retrieve the abundances of the detected species and constrain planetary orbital velocities, $T$-$P$ profiles, and the carbon-to-oxygen ratio ($\rm C/O$). A free retrieval results in an elevated $\rm CO$ abundance ($\log_{10}$($χ_{\rm{{}^{12}CO}}$) = $-2.85^{+0.57}_{-0.69}$), leading to a super-solar $\rm C/O$ ratio. More realistically, allowing for vertically-varying chemistry in the atmosphere by incorporating a chemical-equilibrium model results in a $\rm C/O$ of $0.68^{+0.12}_{-0.22}$ and a metallicity of $[\rm M/H] = 0.62^{+0.28}_{-0.55}$, both consistent with solar values. Finally, we also report a slight offset of the $\rm Fe$ feature in both K$_{\rm p}$ and v$_{\rm sys}$ that could be a signature of atmospheric dynamics. Due to the 3D structure of exoplanet atmospheres and the exclusion of time/phase dependence in our 1D forward models, further follow-up observations and analysis are required to confirm or refute this result.
△ Less
Submitted 14 August, 2023;
originally announced August 2023.
-
Retrieval survey of metals in six ultra-hot Jupiters: Trends in chemistry, rain-out, ionisation and atmospheric dynamics
Authors:
Siddharth Gandhi,
Aurora Kesseli,
Yapeng Zhang,
Amy Louca,
Ignas Snellen,
Matteo Brogi,
Yamila Miguel,
Núria Casasayas-Barris,
Stefan Pelletier,
Rico Landman,
Cathal Maguire,
Neale P. Gibson
Abstract:
Ground-based high-resolution spectroscopy (HRS) has detected numerous chemical species and atmospheric dynamics in exoplanets, most notably ultra-hot Jupiters (UHJs). However, quantitative estimates on abundances have been challenging but are essential for accurate comparative characterisation and to determine formation scenarios. In this work we retrieve the atmospheres of six UHJs (WASP-76~b, MA…
▽ More
Ground-based high-resolution spectroscopy (HRS) has detected numerous chemical species and atmospheric dynamics in exoplanets, most notably ultra-hot Jupiters (UHJs). However, quantitative estimates on abundances have been challenging but are essential for accurate comparative characterisation and to determine formation scenarios. In this work we retrieve the atmospheres of six UHJs (WASP-76~b, MASCARA-4~b, MASCARA-2~b, WASP-121~b, HAT-P-70~b and WASP-189~b) with ESPRESSO and HARPS-N/HARPS observations, exploring trends in eleven neutral species and dynamics. While Fe abundances agree well with stellar values, Mg, Ni, Cr, Mn and V show more variation, highlighting the difficulty in using a single species as a proxy for metallicity. We find that Ca, Na, Ti and TiO are under-abundant, potentially due to ionisation and/or night-side rain-out. Our retrievals also show that relative abundances between species are more robust, consistent with previous works. We perform spatially- and phase-resolved retrievals for WASP-76~b and WASP-121~b given their high signal-to-noise observations, and find the chemical abundances in each of the terminator regions are broadly consistent. We additionally constrain dynamics for our sample through Doppler shifts and broadening of the planetary signals during the primary eclipse, with median blue shifts between $\sim$0.9-9.0~km/s due to day-night winds. Furthermore, we constrain spectroscopic masses for MASCARA-2~b and HAT-P-70~b consistent with their known upper limits, but we note that these may be biased due to degeneracies. This work highlights the importance of future HRS studies to further probe differences and trends between exoplanets.
△ Less
Submitted 26 May, 2023;
originally announced May 2023.
-
High-resolution atmospheric retrievals of WASP-121b transmission spectroscopy with ESPRESSO: Consistent relative abundance constraints across multiple epochs and instruments
Authors:
Cathal Maguire,
Neale P. Gibson,
Stevanus K. Nugroho,
Swaetha Ramkumar,
Mark Fortune,
Stephanie R. Merritt,
Ernst de Mooij
Abstract:
Recent progress in high-resolution transmission spectroscopy has offered new avenues in which to characterise the atmospheres of transiting exoplanets. High-resolution cross-correlation spectroscopy allows for the unambiguous detection of molecules/atoms. It has also been used to map both atmospheric dynamics and longitudinal variations in the abundance of species across the morning and evening li…
▽ More
Recent progress in high-resolution transmission spectroscopy has offered new avenues in which to characterise the atmospheres of transiting exoplanets. High-resolution cross-correlation spectroscopy allows for the unambiguous detection of molecules/atoms. It has also been used to map both atmospheric dynamics and longitudinal variations in the abundance of species across the morning and evening limbs. We present multiple VLT/ESPRESSO observations of the ultra-hot Jupiter WASP-121b, from which we constrain relative abundances of various neutral metals consistently across all observations, whilst accounting for the distortion of the exoplanet's signal caused by traditional data processing techniques. We also constrain planetary orbital velocities and $T$-$P$ profiles. We compare our abundance constraints with previous constraints using VLT/UVES transmission spectroscopy of WASP-121b, and find our results to be consistent between observations, and also in agreement with stellar values for species previously detected in the atmosphere of WASP-121b. Our retrieval framework can also be used to identify potential exospheric species, resulting in extended absorption features beyond the transit equivalent Roche limit of WASP-121b ($R_{\rm eqRL}$ $\sim$ 1.3 $R_{\rm p}$). H$α$, Fe II, and Ca II were found to extend to high altitudes ($1.54\pm0.04$ $R_{\rm p}$, $1.17\pm0.01$ $R_{\rm p}$, and $2.52\pm0.34$ $R_{\rm p}$, respectively), which are broadly consistent with literature values. The consistency of our constraints across multiple high-resolution observations is a strong validation of our model filtering and retrieval framework, as well as the stability of the atmosphere over the timescales of months/years, and could allow for planet formation processes to be inferred from future ground-based observations of exoplanetary atmospheres.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
Relative abundance constraints from high-resolution optical transmission spectroscopy of WASP-121b, and a fast model-filtering technique for accelerating retrievals
Authors:
Neale P. Gibson,
Stevanus K. Nugroho,
Joshua Lothringer,
Cathal Maguire,
David K. Sing
Abstract:
High-resolution Doppler-resolved spectroscopy has presented new opportunities for studying the atmospheres of exoplanets. While the 'classical' cross-correlation approach has proven to be efficient at finding atmospheric species, it is unable to perform direct atmospheric retrievals. Recent work has shown that retrievals are possible using a direct likelihood evaluation or likelihood 'mappings'. T…
▽ More
High-resolution Doppler-resolved spectroscopy has presented new opportunities for studying the atmospheres of exoplanets. While the 'classical' cross-correlation approach has proven to be efficient at finding atmospheric species, it is unable to perform direct atmospheric retrievals. Recent work has shown that retrievals are possible using a direct likelihood evaluation or likelihood 'mappings'. The unique aspect of high-resolution retrievals is that the data-processing methods required to remove the stellar and telluric lines also distort the underlying exoplanet's signal and therefore the forward model must be pre-processed to match this filtering. This was the key remaining limitation in our previously published framework. This paper directly addresses this by introducing a simple and fast model-filtering technique that can replicate the processing performed by algorithms such as SysRem and PCA. This enables retrievals to be performed without having to perform expensive injection and pre-processing steps for every model. We show that we can reliably constrain quantitative measures of the atmosphere from transmission spectra including the temperature-pressure profile, relative abundances, planetary velocities and rotational broadening parameters. Finally, we demonstrate our framework using UVES transmission spectroscopy of WASP-121b. We constrain the temperature-pressure profile and relative abundances of Fe, Cr, and V to be $\log_{10}(χ_{\rm Fe}/χ_{\rm Cr})$=1.66$\pm$0.28, $\log_{10}(χ_{\rm Fe}/χ_{\rm V})$=3.78$\pm$0.29 and $\log_{10}(χ_{\rm Fe}/χ_{\rm Mg})$=-1.26$\pm$0.60. The relative abundances are consistent with solar values, with the exception of Fe/Mg, where the large Mg abundance is probably explained by the escaping atmosphere of WASP-121b that is not accounted for in our atmospheric model.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
First Results from the REAL-time Transient Acquisition backend (REALTA) at the Irish LOFAR station
Authors:
P. C. Murphy,
P. Callanan,
J. McCauley,
D. J. McKenna,
D. Ó Fionnagáin,
C. K. Louis,
M. P. Redman,
L. A. Cañizares,
E. P. Carley,
S. A. Maloney,
B. Coghlan,
M. Daly,
J. Scully,
J. Dooley,
V. Gajjar,
C. Giese,
A. Brennan,
E. F. Keane,
C. A. Maguire,
J. Quinn,
S. Mooney,
A. M. Ryan,
J. Walsh,
C. M. Jackman,
A. Golden
, et al. (5 additional authors not shown)
Abstract:
Modern radio interferometers such as the LOw Frequency ARray (LOFAR) are capable of producing data at hundreds of gigabits to terabits per second. This high data rate makes the analysis of radio data cumbersome and computationally expensive. While high performance computing facilities exist for large national and international facilities, that may not be the case for instruments operated by a sing…
▽ More
Modern radio interferometers such as the LOw Frequency ARray (LOFAR) are capable of producing data at hundreds of gigabits to terabits per second. This high data rate makes the analysis of radio data cumbersome and computationally expensive. While high performance computing facilities exist for large national and international facilities, that may not be the case for instruments operated by a single institution or a small consortium. Data rates for next generation radio telescopes are set to eclipse those currently in operation, hence local processing of data will become all the more important. Here, we introduce the REAL-time Transient Acquisition backend (REALTA), a computing backend at the Irish LOFAR station (I-LOFAR) which facilitates the recording of data in near real-time and post-processing. We also present first searches and scientific results of a number of radio phenomena observed by I-LOFAR and REALTA, including pulsars, fast radio bursts (FRBs), rotating radio transients (RRATs), the search for extraterrestrial intelligence (SETI), Jupiter, and the Sun.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
LOFAR observations of a jet-driven piston shock in the low solar corona
Authors:
Ciara A. Maguire,
Eoin P. Carley,
Pietro Zucca,
Nicole Vilmer,
Peter T. Gallagher
Abstract:
The Sun produces highly dynamic and eruptive events that can drive shocks through the corona. These shocks can accelerate electrons, which result in plasma emission in the form of a type II radio burst. Despite the large number of type II radio bursts observations, the precise origin of coronal shocks is still subject to investigation. Here we present a well observed solar eruptive event that occu…
▽ More
The Sun produces highly dynamic and eruptive events that can drive shocks through the corona. These shocks can accelerate electrons, which result in plasma emission in the form of a type II radio burst. Despite the large number of type II radio bursts observations, the precise origin of coronal shocks is still subject to investigation. Here we present a well observed solar eruptive event that occurred on 16 October 2015, focusing on a jet observed in the extreme ultraviolet (EUV) by the Atmospheric Imaging Assembly (SDO/AIA), a streamer observed in white-light by the Large Angle and Spectrometric Coronagraph (SOHO/LASCO), and a metric type II radio burst observed by the LOw Frequency Array (LOFAR). LOFAR interferometrically imaged the fundamental and harmonic sources of a type II radio burst and revealed that the sources did not appear to be co-spatial, as would be expected from the plasma emission mechanism. We correct for the separation between the fundamental and harmonic using a model which accounts for scattering of radio waves by electron density fluctuations in a turbulent plasma. This allows us to show the type II radio sources were located $\sim$0.5 R$_\odot$ above the jet and propagated at a speed of $\sim$1000 kms$^{-1}$, which was significantly faster than the jet speed of $\sim$200 kms$^{-1}$. This suggests that the type II burst was generated by a piston shock driven by the jet in the low corona.
△ Less
Submitted 14 January, 2021;
originally announced January 2021.
-
Evolution of the Alfvén Mach number associated with a coronal mass ejection shock
Authors:
Ciara A. Maguire,
Eoin P. Carley,
Joseph McCauley,
Peter T. Gallagher
Abstract:
The Sun regularly produces large-scale eruptive events, such as coronal mass ejections (CMEs) that can drive shock waves through the solar corona. Such shocks can result in electron acceleration and subsequent radio emission in the form of a type II radio burst. However, the early-phase evolution of shock properties and its relationship to type II burst evolution is still subject to investigation.…
▽ More
The Sun regularly produces large-scale eruptive events, such as coronal mass ejections (CMEs) that can drive shock waves through the solar corona. Such shocks can result in electron acceleration and subsequent radio emission in the form of a type II radio burst. However, the early-phase evolution of shock properties and its relationship to type II burst evolution is still subject to investigation. Here we study the evolution of a CME-driven shock by comparing three commonly used methods of calculating the Alfvén Mach number ($M_A$), namely: shock geometry, a comparison of CME speed to a model of the coronal Alfvén speed, and the type II band-splitting method. We applied the three methods to the 2017 September 2 event, focusing on the shock wave observed in extreme ultraviolet (EUV) by the Solar Ultraviolet Imager (SUVI) on board GOES-16, in white-light by the Large Angle and Spectrometric Coronagraph (LASCO) on board SOHO, and the type II radio burst observed by the Irish Low Frequency Array (I-LOFAR). We show that the three different methods of estimating shock $M_A$ yield consistent results and provide a means of relating shock property evolution to the type II emission duration. The type II radio emission emerged from near the nose of the CME when $M_A$ was in the range 1.4-2.4 at a heliocentric distance of $\sim$1.6 $R_\odot$. The emission ceased when the CME nose reached $\sim$2.4 $R_\odot$, despite an increasing Alfvén Mach number (up to 4). We suggest the radio emission cessation is due to the lack of quasi-perpendicular geometry at this altitude, which inhibits efficient electron acceleration and subsequent radio emission.
△ Less
Submitted 6 December, 2019; v1 submitted 4 December, 2019;
originally announced December 2019.