Search

Modeling based on differentiable programming holds great promise for astronomy, as it can employ techniques such as Hamiltonian Monte Carlo, gradient-based optimization, and other machine learning techniques. This new programming paradigm has motivated us to develop the first auto-differentiable spectrum model of exoplanets and brown dwarfs, ExoJAX (Kawahara et al. 2022). ExoJAX is designed to directly calculate cross-sections as functions of temperature and pressure, rather than interpolating tabulated data, to minimize errors in high-dispersion spectra modeling. However, its application was primarily proof-of-concept and limited to narrowband high-dispersion emission spectroscopy. In this paper, we have enhanced the differentiable opacity calculation using a new fast and memory-efficient algorithm, and have developed differentiable radiative transfer schemes, including emission, transmission, and reflection spectroscopy. These enhancements significantly expand the range of applications, as demonstrated through actual atmospheric retrievals: high-dispersion emission spectra of the brown dwarf GL229 B, medium-dispersion transmission spectra of the hot Saturn WASP-39 b from JWST, and high-dispersion reflection spectra of Jupiter. We obtained a C/O ratio for GL229 B consistent with its host star, constrained WASP-39 b's radial velocity from molecular fine structures at original resolution ($R \sim 3,000$), and estimated Jupiter's metallicity consistent with previous studies., Comment: 32 pages, 12 figures, submitted. ExoJAX is available at https://github.com/HajimeKawahara/exojax. ExoJAX2 is scheduled to be released after the review process

Optical phase curves of hot Jupiters can reveal global scattering properties. We implement a Bayesian inference framework for optical phase curves with flux contributions from: reflected light from a potentially inhomogeneous atmosphere, thermal emission, ellipsoidal variations, Doppler beaming, and stellar rotation via a Gaussian process in the time domain. We probe for atmospheric homogeneity and time-variability using the reflected light inferences for highly precise Kepler light curves of five hot Jupiters. We also investigate the scattering properties which constrain the most likely condensates in the inhomogeneous atmospheres. Cross validation prefers inhomogeneous albedo distributions for Kepler-7 b and Kepler-41 b, and a weak preference for inhomogeneity for KOI-13 b. None of the five planets exhibit significant variations in geometric albedo on one-year timescales, in agreement with theoretical expectations. We show that analytic reflected light phase curves with isotropic multiple scattering are in excellent agreement with full Rayleigh multiple scattering calculations, allowing for accelerated and analytic inference. In a case study of Kepler-41 b, we identify perovskite, forsterite, and enstatite as possible scattering species consistent with the reflected light phase curves, with condensate particle radii in the range 0.01-0.1 micron., Comment: Accepted in A&A

Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.

Ultra-hot Jupiters are tidally locked with their host stars dividing their atmospheres into a hot dayside and a colder nightside. As the planet moves through transit, different regions of the atmosphere rotate into view revealing different chemical regimes. High-resolution spectrographs can observe asymmetries and velocity shifts, and offer the possibility for time-resolved spectroscopy. In this study, we search for other atoms and molecules in the planet`s transmission spectrum and investigate asymmetric signals. We analyse and combine eight transits of the ultra-hot Jupiter WASP-189 b taken with the HARPS, HARPS-N, ESPRESSO and MAROON-X high-resolution spectrographs. Using the cross-correlation technique, we search for neutral and ionised atoms, and oxides and compare the obtained signals to model predictions. We report significant detections for H, Na, Mg, Ca, Ca+, Ti, Ti+, TiO, V, Cr, Mn, Fe, Fe+, Ni, Sr, Sr+, and Ba+. Of these, Sr, Sr+, and Ba+ are detected for the first time in the transmission spectrum of WASP-189 b. In addition, we robustly confirm the detection of titanium oxide based on observations with HARPS and HARPS-N using the follow-up observations performed with MAROON-X and ESPRESSO. By fitting the orbital traces of the detected species by means of time-resolved spectroscopy using a Bayesian framework, we infer posterior distributions for orbital parameters as well as lineshapes. Our results indicate that different species must originate from different regions of the atmosphere to be able to explain the observed time dependence of the signals. Throughout the course of the transit, most signal strengths are expected to increase due to the larger atmospheric scale height at the hotter trailing terminator. For some species, however, the signals are instead observed to weaken due to ionisation for atoms and their ions, or the dissociation of molecules on the dayside., Comment: 38 pages, 34 figures, published in A&A on October 24, 2023

We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $\rho=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$., Comment: Published in ApJL

The distribution of small-scale magnetic fields in stellar photospheres is an important ingredient in our understanding of the magnetism of low mass stars. Their spatial distribution connects the field generated in the stellar interior with the outer corona and the large scale field, and thereby affects the space weather of planets. Unfortunately, we lack techniques that can locate them on most low-mass stars. One strategy is to localize field concentrations using the flares that occur in their vicinity. We explore a new method that adapts the spot simulation software fleck to study the modulation of flaring times as a function of active latitude. We use empirical relations to construct flare light curves similar to those available from Kepler and the Transiting Exoplanet Survey Satellite (TESS), search them for flares, and use the waiting times between flares to determine the location of active latitudes. We find that the mean and standard deviation of the waiting time distribution provide a unique diagnostic of flaring latitudes as a function of the number of active regions. Latitudes are best recovered when stars have three or less active regions that flare repeatedly, and active latitude widths below 20 deg; when either increases, the information about the active latitude location is gradually lost. We demonstrate our technique on a sample of flaring G dwarfs observed with the Kepler satellite, and furthermore suggest that combining ensemble methods for spots and flares could overcome the limitations of each individual technique for the localization of surface magnetic fields., Comment: 15 pages, 10 figures. Accepted to MNRAS. Repository with source code: https://github.com/ekaterinailin/flare-locations-ensembles-science Repository with data: https://zenodo.org/record/7996929

Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the magnetic flux carried by flare-associated CMEs results in planetary interior heating by Ohmic dissipation and leads to a variety of interior--exterior interactions. We construct a physical model to study this effect and apply it to the TRAPPIST-1 star whose flaring activity has been constrained by Kepler observations. Our model is posed in a stochastic manner to account for uncertainty and variability in input parameters. Particularly for the innermost planets, our results suggest that the heat dissipated in the silicate mantle is both of sufficient magnitude and longevity to drive geological processes and hence facilitate volcanism and outgassing of the TRAPPIST-1 planets. Furthermore, our model predicts that Joule heating can further be enhanced for planets with an intrinsic magnetic field compared to those without. The associated volcanism and outgassing may continuously replenish the atmosphere and thereby mitigate the erosion of the atmosphere caused by the direct impact of flares and CMEs. To maintain consistency of atmospheric and geophysical models, the impact of stellar flares and CMEs on atmospheres of close-in exoplanetary systems needs to be studied in conjunction with the effect on planetary interiors.

The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project., Comment: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

Fast yellow pulsating supergiants (FYPS) are a recently-discovered class of evolved massive pulsator. As candidate post-red supergiant objects, and one of the few classes of pulsating evolved massive stars, these objects have incredible potential to change our understanding of the structure and evolution of massive stars. Here we examine the lightcurves of a sample of 126 cool supergiants in the Magellanic Clouds observed by the Transiting Exoplanet Survey Satellite (\tess~) in order to identify pulsating stars. After making quality cuts and filtering out contaminant objects, we examine the distribution of pulsating stars in the Hertzprung-Russel (HR) diagram, and find that FYPS occupy a region above $\log L/L_\odot \gtrsim 5.0$. This luminosity boundary corresponds to stars with initial masses of $\sim$18-20 $M_\odot$, consistent with the most massive red supergiant progenitors of supernovae (SNe) II-P, as well as the observed properties of SNe IIb progenitors. This threshold is in agreement with the picture that FYPS are post-RSG stars. Finally, we characterize the behavior of FYPS pulsations as a function of their location in the HR diagram. We find low frequency pulsations at higher effective temperatures, higher frequency pulsations at lower temperatures, with a transition between the two behaviors at intermediate temperatures. The observed properties of FYPS make them fascinating objects for future theoretical study., Comment: Consistent with published version which contains significantly improved detection and rejection of contaminant objects. Comments welcome

Understanding magnetic activity on the surface of stars other than the Sun is important for exoplanet analyses to properly characterize an exoplanet's atmosphere and to further characterize stellar activity on a wide range of stars. Modeling stellar surface features of a variety of spectral types and rotation rates are key to understanding of the magnetic activity of these stars. Using data from Kepler, we use the starspot modeling program STarSPot (STSP) to measure the position and size of spots for KOI-340 which is an eclipsing binary consisting of a subgiant star ($\rm T_{\rm eff} = 5593 \pm 27 K; R_{*} = 1.98 \pm 0.05 R_{\odot}$) with an M-dwarf companion ($M_{*} = 0.214 \pm 0.006 M_{\odot}$). STSP uses a novel technique to measure the spot positions and radii by using the transiting secondary to study and model individual active regions on the stellar surface using high-precision photometry. We find the average size of spot features on KOI-340's primary is $\sim$10% the radius of the star, i.e. two times larger than the mean size of Solar-maximum sunspots. The spots on KOI-340 are present at every longitude and show possible signs of differential rotation. The minimum fractional spotted area of KOI-340's primary is $2\substack{+12\\ -2} \%$ while the spotted area of the Sun is at most 0.2%. One transit of KOI-340 shows a signal in the transit consistent with a plage; this plage occurs right before a dark spot indicating the plage and spot might be co-located on the surface of the star., Comment: 16 pages, 17 figures; Accepted to AJ

General circulation models (GCMs) provide context for interpreting multi-wavelength, multi-phase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K) and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR+HELIOS WASP-43b GCM (horizontal resolution of about 4 degrees on the sphere and with 40 radial points) with multi-wavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864,000 time steps) takes about 19-26 days to complete depending on the type of GPU., Comment: 31 pages, 24 figures, accepted for publication at MNRAS

We explore the transit timing variations (TTVs) of the young (22 Myr) nearby AU Mic planetary system. For AU Mic b, we introduce three Spitzer (4.5 $\mu$m) transits, five TESS transits, 11 LCO transits, one PEST transit, one Brierfield transit, and two transit timing measurements from Rossiter-McLaughlin observations; for AU Mic c, we introduce three TESS transits. We present two independent TTV analyses. First, we use EXOFASTv2 to jointly model the Spitzer and ground-based transits and to obtain the midpoint transit times. We then construct an O--C diagram and model the TTVs with Exo-Striker. Second, we reproduce our results with an independent photodynamical analysis. We recover a TTV mass for AU Mic c of 10.8$^{+2.3}_{-2.2}$ M$_{E}$. We compare the TTV-derived constraints to a recent radial-velocity (RV) mass determination. We also observe excess TTVs that do not appear to be consistent with the dynamical interactions of b and c alone, and do not appear to be due to spots or flares. Thus, we present a hypothetical non-transiting "middle-d" candidate exoplanet that is consistent with the observed TTVs, the candidate RV signal, and would establish the AU Mic system as a compact resonant multi-planet chain in a 4:6:9 period commensurability. These results demonstrate that the AU Mic planetary system is dynamically interacting producing detectable TTVs, and the implied orbital dynamics may inform the formation mechanisms for this young system. We recommend future RV and TTV observations of AU Mic b and c to further constrain the masses and to confirm the existence of possible additional planet(s)., Comment: Corrected typos; revised Section 3, 4, and 5 to reflect reanalysis, results unchanged. Submitted to AAS Journals Nov 11th, 2020; favorable referee report received Jan 3rd; final draft accepted for publication in the AJ Apr 19th

Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature., Comment: Invited review in press at RASTI. Based on the ExoPAG SAG21 report (arXiv:2201.09905v1) and refined via feedback from three reviewers. 75 pages, 30 figures, 5 tables