Papers for Friday, Feb 05 2021

Alongside future observations with the new European Extremely Large Telescope (ELT), optimised instruments on the 8-10m generation of telescopes will still be competitive at 'ground UV' wavelengths (3000-4000 A). The near UV provides a wealth of unique information on the nucleosynthesis of iron-peak elements, molecules, and neutron-capture elements. In the context of development of the near-UV CUBES spectrograph for ESO's Very Large Telescope (VLT), we are investigating the impact of spectral resolution on the ability to estimate chemical abundances for beryllium and more than 30 iron-peak and heavy elements. From work ahead of the Phase A conceptual design of CUBES, here we present a comparison of the elements observable at the notional resolving power of CUBES (R~20,000) to those with VLT-UVES (R~40,000). For most of the considered lines signal-to-noise is a more critical factor than resolution. We summarise the elements accessible with CUBES, several of which (e.g. Be, Ge, Hf) are now the focus of quantitative simulations as part of the ongoing Phase A study.

Without additional heating, radiative cooling of gas in the halos of massive galaxies (Milky Way and above) produces cold gas or stars in excess of that observed. Previous work suggested that AGN jets are likely required, but the form of jet energy required to quench remains unclear. This is particularly challenging for galaxy simulations, in which the resolution is orders of magnitude coarser than necessary to form and evolve the jet. On such scales, the uncertain parameters include: jet energy form (kinetic, thermal, and cosmic ray (CR) energy), energy, momentum, and mass flux, magnetic field strength and geometry, jet precession angle and period, opening-angle, and duty cycle. We investigate all of these parameters in a $10^{14}\,{\rm M}_{\odot}$ halo using high-resolution non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios match observational constraints and show that CR-dominated jets can most efficiently quench the central galaxy through a combination of CR pressure support and a modification of the thermal instability. Jets with most energy in mildly relativistic ($\sim$ MeV or $\sim10^{10}$ K) thermal plasma work, but require a factor $\sim 10$ larger energy input. For a fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo size, kinetic jets are less efficient in quenching unless they have wide opening or precession angles. Magnetic fields play a minor role except when the magnetic flux reaches $\gtrsim 10^{44}$ erg s$^{-1}$ in a kinetic jet model, which causes the jet cocoon to significantly widen, and the quenching to become explosive. We conclude that the criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with long enough cooling time at the cooling radius.

Sterile neutrinos with masses in the keV range are well-motivated extensions to the Standard Model that could explain the observed neutrino masses while also making up the dark matter (DM) of the Universe. If sterile neutrinos are DM then they may slowly decay into active neutrinos and photons, giving rise to the possibility of their detection through narrow spectral features in astrophysical X-ray data sets. In this work, we perform the most sensitive search to date for this and other decaying DM scenarios across the mass range from 5 to 16 keV using archival XMM-Newton data. We reduce 547 Ms of data from both the MOS and PN instruments using observations taken across the full sky and then use this data to search for evidence of DM decay in the ambient halo of the Milky Way. We determine the instrumental and astrophysical baselines with data taken far away from the Galactic Center, and use Gaussian Process modeling to capture additional continuum background contributions. No evidence is found for unassociated X-ray lines, leading us to produce the strongest constraints to date on decaying DM in this mass range.

To determine the importance of merging galaxies to galaxy evolution, it is necessary to design classification tools that can identify different types and stages of merging galaxies. Previously, using GADGET-3/SUNRISE simulations of merging galaxies and linear discriminant analysis (LDA), we created an accurate merging galaxy classifier from imaging predictors. Here, we develop a complementary tool based on stellar kinematic predictors derived from the same simulation suite. We design mock stellar velocity and velocity dispersion maps to mimic the specifications of the Mapping Nearby Galaxies at Apache Point (MaNGA) integral field spectroscopy (IFS) survey and utilize an LDA to create a classification based on a linear combination of 11 kinematic predictors. The classification varies significantly with mass ratio; the major (minor) merger classifications have a mean statistical accuracy of 80% (70%), a precision of 90% (85%), and a recall of 75% (60%). The major mergers are best identified by predictors that trace global kinematic features, while the minor mergers rely on local features that trace a secondary stellar component. While the kinematic classification is less accurate than the imaging classification, the kinematic predictors are better at identifying post-coalescence mergers. A combined imaging + kinematic classification has the potential to reveal more complete merger samples from imaging and IFS surveys like MaNGA. We note that since the suite of simulations used to train the classifier covers a limited range of galaxy properties (i.e., the galaxies are intermediate mass and disk-dominated), the results may not be applicable to all MaNGA galaxies.

Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib $\&$ Lin (2021), we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the "tail" of the velocity distribution, and the consistency of the data with different underlying models. In this paper, we apply this method to the second data release of Gaia and find that a model with at least two components is preferred. Based on a fit with three bound components to account for the disk, relaxed halo, and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be $v_{\rm{esc}}= 484.6^{+17.8}_{-7.4}$ km/s. Assuming a Navarro-Frenck-White dark matter profile, and taken in conjunction with a recent measurement of the circular velocity at the solar position of $v_c = 230 \pm 10$ km/s, we find a Milky Way concentration of $c_{200} = 13.8^{+6.0}_{-4.3}$ and a mass of $M_{200} = 7.0^{+1.9}_{-1.2} \times 10^{11} M_{\odot}$, which is considerably lighter than previous measurements.

We describe a new algorithm to solve the time dependent, frequency integrated radiation transport (RT) equation implicitly, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using {\sf Athena++}. The radiation filed is represented by specific intensities along discrete rays, which are evolved using a conservative finite volume approach for both cartesian and curvilinear coordinate systems. All the terms for spatial transport of photons and interactions between gas and radiation are calculated implicitly together. An efficient Jacobi-like iteration scheme is used to solve the implicit equations. This removes any time step constrain due to the speed of light in RT. We evolve the specific intensities in the lab frame to simplify the transport step. The lab-frame specific intensities are transformed to the co-moving frame via Lorentz transformation when the source term is calculated. Therefore, the scheme does not need any expansion in terms of $v/c$. The radiation energy and momentum source terms for the gas are calculated via direct quadrature in the angular space. The time step for the whole scheme is determined by the normal Courant -- Friedrichs -- Lewy condition in the MHD module. We provide a variety of test problems for this algorithm including both optically thick and thin regimes, and for both gas and radiation pressure dominated flows to demonstrate its accuracy and efficiency.

We find that, under certain conditions, protoplanetary disks may spontaneously generate multiple, concentric gas rings without an embedded planet through an eccentric cooling instability. Using both linear theory and non-linear hydrodynamics simulations, we show that a variety of background states may trap a slowly precessing, one-armed spiral mode that becomes unstable when a gravitationally-stable disk rapidly cools. The angular momentum required to excite this spiral comes at the expense of non-uniform mass transport that generically results in multiple rings. For example, one long-term hydrodynamics simulation exhibits four long-lived, axisymmetric gas rings. We verify the instability evolution and ring formation mechanism from first principles with our linear theory, which shows remarkable agreement with the simulation results. Dust trapped in these rings may produce observable features consistent with observed disks. Additionally, direct detection of the eccentric gas motions may be possible when the instability saturates, and any residual eccentricity leftover in the rings at later times may also provide direct observational evidence of this mechanism.

Context. The paths followed by the known extreme trans-Neptunian objects (ETNOs) effectively avoid direct gravitational perturbations from the four giant planets, yet their orbital eccentricities are in the range between 0.69-0.97. Solar system dynamics studies show that such high values of the eccentricity can be produced via close encounters or secular perturbations. In both cases, the presence of yet-to-be-discovered trans-Plutonian planets is required. Aims. If the high eccentricities of the known ETNOs are the result of relatively recent close encounters with putative planets, the mutual nodal distances of sizeable groups of ETNOs with their assumed perturber may still be small enough to be identifiable geometrically. In order to confirm or reject this possibility, we used Monte Carlo random search techniques. Methods. Two arbitrary orbits may lead to close encounters when their mutual nodal distance is sufficiently small. We generated billions of random planetary orbits with parameters within the relevant ranges and computed the mutual nodal distances with a set of randomly generated orbits with parameters consistent with those of the known ETNOs and their uncertainties. We monitored which planetary orbits had the maximum number of potential close encounters with synthetic ETNOs and we studied the resulting distributions. Results. We provide narrow ranges for the orbital parameters of putative planets that may have experienced orbit-changing encounters with known ETNOs. Conclusions. Our calculations suggest that more than one perturber is required if scattering is the main source of orbital modification for the known ETNOs. Perturbers might not be located farther than 600 AU and they have to follow moderately eccentric and inclined orbits to be capable of experiencing close encounters with multiple known ETNOs.

We report the discovery of two short-period massive giant planets from NASA's Transiting Exoplanet Survey Satellite (TESS). Both systems, TOI-558 (TIC 207110080) and TOI-559 (TIC 209459275), were identified from the 30-minute cadence Full Frame Images and confirmed using ground-based photometric and spectroscopic follow-up observations from TESS's Follow-up Observing Program Working Group. We find that TOI-558 b, which transits an F-dwarf ($M_{\star}=1.349^{+0.064}_{-0.065}\ M_{\odot}$, $R_{\star}=1.496^{+0.042}_{-0.040}\ R_{\odot}$, $T_{eff}=6466^{+95}_{-93}$ K, age $1.79^{+0.91}_{-0.73}$ Gyr) with an orbital period of 14.574 days, has a mass of $3.61\pm0.15\ M_J$, a radius of $1.086^{+0.041}_{-0.038}\ R_J$, and an eccentric (e=$0.300^{+0.022}_{-0.020}$) orbit. TOI-559 b transits a G-dwarf ($M_{\star}=1.026\pm0.057\ M_{\odot}$, $R_{\star}=1.233^{+0.028}_{-0.026}\ R_{\odot}$, $T_{eff}=5925^{+85}_{-76}$ K, age $1.79^{+0.91}_{-0.73}$ Gyr) in an eccentric (e=$0.151\pm0.011$) 6.984-day orbit with a mass of $6.01^{+0.24}_{-0.23}\ M_J$ and a radius of $1.091^{+0.028}_{-0.025}\ R_J$. Our spectroscopic follow-up also reveals a long-term radial velocity trend for TOI-559, indicating a long-period companion. The statistically significant orbital eccentricity measured for each system suggests that these planets migrated to their current location through dynamical interactions. Interestingly, both planets are also massive ($>3\ M_J$), adding to the population of massive hot Jupiters identified by TESS. Prompted by these new detections of high-mass planets, we analyzed the known mass distribution of hot Jupiters but find no significant evidence for multiple populations. TESS should provide a near magnitude-limited sample of transiting hot Jupiters, allowing for future detailed population studies.

Long duration gamma-ray bursts (GRBs) are among the least understood astrophysical transients powering the high-energy universe. To date, various mechanisms have been proposed to explain the observed electromagnetic GRB emission. In this work, we show that, although different jet models may be equally successful in fitting the observed electromagnetic spectral energy distributions, the neutrino production strongly depends on the adopted emission and dissipation model. To this purpose, we compute the neutrino production for a benchmark high-luminosity GRB in the internal shock model, including a dissipative photosphere as well as three emission components, in the jet model invoking internal-collision-induced magnetic reconnection and turbulence (ICMART), in the case of a magnetic jet with gradual dissipation, and in a jet with dominant proton synchrotron radiation. We find that the expected neutrino fluence can vary up to 1-1.5 orders of magnitude in amplitude and peak at energies ranging from $10^4$ to $10^8$ GeV. For our benchmark input parameters, none of the explored GRB models is excluded by the targeted searches carried out by the IceCube and ANTARES Collaborations. However, our work highlights the potential of high-energy neutrinos of pinpointing the underlying GRB emission mechanism and the importance of relying on different jet models for unbiased stacking searches.

(abridged) Data and results from the WISH key program are summarized, designed to provide a legacy data set to address its physics and chemistry. WISH targeted ~80 sources along the two axes of luminosity and evolutionary stage: from low- to high-mass protostars and from pre-stellar cores to protoplanetary disks. Lines of H2O, HDO, OH, CO and [O I] were observed with the HIFI and PACS instruments, complemented by molecules that probe UV, X-ray or grain chemistry. Most of the far-infrared water emission from protostars is found to be compact, originating from warm outflowing and shocked gas at high density and temperature in at least two physical components. This gas is not probed by low-J CO lines, only by J>14. Water is a significant, but not dominant, coolant. Its abundance is universally low, of order H2O/H2=2E-6, pointing to shock and outflow cavity models that include UV radiation at 100-1000 times the ISRF. In cold quiescent pre-stellar cores and envelopes, the water abundance structure is accurately probed through velocity-resolved line profiles, confirming basic chemistry networks. The gaseous HDO/H2O ratio of 0.025, much higher than that of bulk ice, is representative of the outer photodesorbed ice layers and cold chemistry. Water abundances in the inner hot cores are high, but with variations from 5E-6 to 2E-4. Combined analyses of water gas and ice show that up to 50% of the oxygen budget may be missing, with possible explanations discussed. Water vapor emission from disks is weak, indicating that water ice is locked up in larger pebbles early on and that these pebbles have settled and drifted inward by the Class II stage. Quantitatively, many oceans of water ice are available. Extragalactic low-J H2O emission is mostly compact and collisionally excited. Prospects for future mid- to far-infrared missions are given.

We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public MOSFiT package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly ($B>10^{10}$) favours the inclusion of an additional luminosity source during the first day, well fit by a shock-heated cocoon, alongside lanthanide-poor dynamical ejecta. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is $q=0.92\pm0.07$ (90% credible interval). We place tight constraints on the maximum stable NS mass, $M_{\rm TOV}=2.17^{+0.08}_{-0.11}$ M$_\odot$. For a uniform prior in tidal deformability, the radius of a 1.4 M$_\odot$ NS is $R_{1.4}\sim 10.7$ km. Re-weighting with a prior based on equations of state that support our credible range in $M_{\rm TOV}$, we derive a final measurement $R_{1.4}=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally-detected neutron star merger, GW190425, we estimate that an associated kilonova would have been fainter (by $\sim0.7$ mag at one day post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.

The coherent low-frequency radio emission detected by LOFAR from Gliese 1151, a quiescent M4.5 dwarf star, has radio emission properties consistent with theoretical expectations of star-planet interactions for an Earth-sized planet on a 1-5 day orbit. New near-infrared radial velocities from the Habitable-zone Planet Finder (HPF) spectrometer on the 10m Hobby-Eberly Telescope at McDonald Observatory, combined with previous velocities from HARPS-N, reveal a periodic Doppler signature consistent with an $m\sin i = 2.5 \pm 0.5 M_\oplus$ exoplanet on a 2.02-day orbit. Precise photometry from the Transiting Exoplanet Survey Satellite (TESS) shows no flares or activity signature, consistent with a quiescent M dwarf. While no planetary transit is detected in the TESS data, a weak photometric modulation is detectable in the photometry at a $\sim2$ day period. This independent detection of a candidate planet signal with the Doppler radial-velocity technique adds further weight to the claim of the first detection of star-exoplanet interactions at radio wavelengths, and helps validate this emerging technique for the detection of exoplanets.

We explore the stellar mass density and colour profiles of 118 low redshift, massive, central galaxies, selected to have assembled 90 percent of their stellar mass 6 Gyr ago, finding evidence of the minor merger activity expected to be the driver behind the size growth of quiescent galaxies. We use imaging data in the $g, r, i, z, y$ bands from the Subaru Hyper Suprime-Cam survey and perform SED fitting to construct spatially well-resolved radial profiles in colour and stellar mass surface density. Our visual morphological classification reveals that $\sim 42$ percent of our sample displays tidal features, similar to previous studies, $\sim 43$ percent of the remaining sample display a diffuse stellar halo and only $\sim 14$ percent display no features, down to a limiting $\mu_{r\mathrm{-band}}$ $\sim$ 28 mag arcsec$^{-2}$. We find good agreement between the stacked colour profiles of our sample to those derived from previous studies and an expected smooth, declining stellar mass surface density profile in the central regions (< 3 R$_{\mathrm{e}}$). However, we also see a flattening of the profile ($\Sigma_* \sim 10^{7.5}$ M$_\odot$ kpc$^{-2}$) in the outskirts (up to 10 R$_{\mathrm{e}}$), which is revealed by our method of specifically targeting tidal/accretion features. We find similar levels of tidal features and behaviour in the stellar mass surface density profiles in a younger comparison sample, however a lack of diffuse haloes. We also apply stacking techniques, similar to those in previous studies, finding such procedures wash out tidal features and thereby produces smooth declining profiles. The stellar material in the outskirts contributes on average $\sim 10^{10}$ M$_\odot$ or a few percent of the total stellar mass and has similar colours to SDSS satellites of similar stellar mass.

The presence of relativistic electrons within the diffuse gas phase of galaxy clusters is now well established, but their detailed origin remains unclear. Cosmic ray protons are also expected to accumulate during the formation of clusters and would lead to gamma-ray emission through hadronic interactions within the thermal gas. Recently, the detection of gamma-ray emission has been reported toward the Coma cluster with Fermi-LAT. Assuming that this gamma-ray emission arises from hadronic interactions in the ICM, we aim at exploring the implication of this signal on the cosmic ray content of the Coma cluster. We use the MINOT software to build a physical model of the cluster and apply it to the Fermi-LAT data. We also consider contamination from compact sources and the impact of various systematic effects. We confirm that a significant gamma-ray signal is observed within the characteristic radius $\theta_{500}$ of the Coma cluster, with a test statistic TS~27 for our baseline model. The presence of a possible point source may account for most of the observed signal. However, this source could also correspond to the peak of the diffuse emission of the cluster itself and extended models match the data better. We constrain the cosmic ray to thermal energy ratio within $R_{500}$ to $X_{\rm CRp}=1.79^{+1.11}_{-0.30}$\% and the slope of the energy spectrum of cosmic rays to $\alpha=2.80^{+0.67}_{-0.13}$. Finally, we compute the synchrotron emission associated with the secondary electrons produced in hadronic interactions assuming steady state. This emission is about four times lower than the overall observed radio signal, so that primary cosmic ray electrons or reacceleration of secondary electrons is necessary to explain the total emission. Assuming an hadronic origin of the signal, our results provide the first quantitative measurement of the cosmic ray proton content in a cluster.[Abridged]

Ground-based Cosmic Microwave Background (CMB) experiments (except those located at the South Pole) take advantage of sky rotation to achieve multiple crossing angles (i.e. scans across each pixel in multiple directions). This crossing-angle coverage is important for controlling the magnitude of some systematics, and for the removal of 1/f noise during the map-making process. Typically, these experiments scan the sky using constant elevation scans, which puts strong constraints on the crossing-angle coverage that is achievable, and therefore on the advantage gained from sky rotation. In this work we elucidate the relationship between scanning elevation and crossing angle, and derive some general "rules of thumb" for scheduling constant elevation scans to maximise the crossing-angle coverage. We derive some bounds on the quality of crossing-angle coverage that can be achieved independently of detailed scheduling choices, which are an important consideration for map-making algorithms if upcoming CMB surveys are to achieve their aims, and we discuss how these results might relate to field selection. These bounds can be used to forecast the effect of scan-coupled systematics for upcoming surveys. We also show quantitatively how some simple choices of boresight rotation improve the possible systematics mitigation from purely sky rotation. Our results are relevant for other surveys that perform constant elevation scans and may have scan-coupled systematics, such as intensity mapping surveys.

We present a recent ALMA observation of the CO(1-0) line emission in the central galaxy of the Zw 3146 galaxy cluster ($z=0.2906$). We also present updated X-ray cavity measurements from archival Chandra observations. The $5\times 10^{10}\,M_{\odot}$ supply of molecular gas, which is confined to the central 4 kpc, is marginally resolved into three extensions that are reminiscent of the filaments observed in similar systems. No velocity structure that would be indicative of ordered motion is observed. The three molecular extensions all trail X-ray cavities, and are potentially formed from the condensation of intracluster gas lifted in the wakes of the rising bubbles. Many cycles of feedback would be require to account for the entire molecular gas reservoir. The molecular gas and continuum source are mutually offset by 2.6 kpc, with no detected line emission coincident with the continuum source. It is the molecular gas, not the continuum source, that lies at the gravitational center of the brightest cluster galaxy. As the brightest cluster galaxy contains possible tidal features, the displaced continuum source may correspond to the nucleus of a merging galaxy. We also discuss the possibility that a gravitational wave recoil following a black hole merger may account for the displacement.

During the early stages of an impact a small amount material may be jetted and ejected at speeds exceeding the impact velocity. Jetting is an important process for producing melt during relatively low velocity impacts. How impact angle affects the jetting process has yet to be fully understood. Here, we simulate jetting during oblique impacts using the iSALE shock physics code. Assuming both the target and impactor have the same composition (dunite), we examine the jetted material which exceeds the impact velocity. Our results show that oblique impacts always produce more jetted ejecta than vertical impacts, except for grazing impacts with impact angles $< 15^{\circ}$. A 45$^{\circ}$ impact with an impact velocity of 3 km/s produces jetted material equal to $\sim$ 7 \% of the impactor mass. This is 6 times the jetted mass produced by a vertical impact with similar impact conditions. We also find that the origin of jetted ejecta depends on impact angle; for impact angles less than 45$^{\circ}$, most of the jet is composed of impactor material, while at higher impact angles the jet is dominated by target material. Our findings are consistent with previous experimental work. In all cases, jetted materials are preferentially distributed downrange of the impactor.

Waiting time distributions of solar flares and {\sl coronal mass ejections (CMEs)} exhibit power law-like distribution functions with slopes in the range of $\alpha_{\tau} \approx 1.4-3.2$, as observed in annual data sets during 4 solar cycles (1974-2012). We find a close correlation between the waiting time power law slope $\alpha_\tau$ and the {\sl sunspot number (SN)}, i.e., $\alpha_\tau$ = 1.38 + 0.01 $\times$ SN. The waiting time distribution can be fitted with a Pareto-type function of the form $N(\tau) = N_0$ $(\tau_0 + \tau)^{-\alpha_{\tau}}$, where the offset $\tau_0$ depends on the instrumental sensitivity, the detection threshold of events, and pulse pile-up effects. The time-dependent power law slope $\alpha_{\tau}(t)$ of waiting time distributions depends only on the global solar magnetic flux (quantified by the sunspot number) or flaring rate, independent of other physical parameters of {\sl self-organized criticality (SOC)} or {\sl magneto-hydrodynamic (MHD)} turbulence models. Power law slopes of $\alpha_{\tau}\approx 1.2-1.6$ were also found in solar wind switchback events, as observed with the {\sl Parker Solar Probe (PSP)}. We conclude that the annual variability of switchback events in the heliospheric solar wind is modulated by flare and CME rates originating in the photosphere and lower corona.

We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB, we consider two different reionization histories for our mock observations: fiducial and late-reionization. For both models we find that bias on the inferred parameter means and 1-$\sigma$ credible intervals can be substantially reduced by using the isosceles bispectrum (calculated for a wide range of scales and triangle shapes) together with the power spectrum (as opposed to just using one of the statistics). We find that making the simplifying assumption of a Gaussian likelihood with a diagonal covariance matrix does not notably bias parameter constraints for the three-parameter reionization model and basic instrumental effects considered here. This is true even if we use extreme (unlikely) initial conditions which would be expected to amplify biases. We also find that using the cosmic variance error calculated with Monte-Carlo simulations using the fiducial model parameters whilst assuming the late-reionization model for the simulated data also does not strongly bias the inference. This implies we may be able to sparsely sample and interpolate the cosmic variance error over the parameter space, substantially reducing computational costs. All codes used in this work are publicly-available.

Although quasi-Keplerian discs are among the most common astrophysical structures, computation of secular angular momentum transport within them routinely presents a considerable practical challenge. In this work, we investigate the secular small-inclination dynamics of a razor-thin particle disc as the continuum limit of a discrete Lagrange-Laplace secular perturbative theory and explore the analogy between the ensuing secular evolution -- including non-local couplings of self-gravitating discs -- and quantum mechanics. We find the 'quantum' Hamiltonian that describes the time evolution of the system and demonstrate the existence of a conserved inner product. The lowest-frequency normal modes are numerically approximated by performing a Wick rotation on the equations of motion. These modes are used to quantify the accuracy of a much simpler local-coupling model, revealing that it predicts the shape of the normal modes to a high degree of accuracy, especially in narrow annuli, even though it fails to predict their eigenfrequencies.

The origin of the inner dust cavities observed in transition discs remains unknown. The segregation of dust and size of the cavity is expected to vary depending on which clearing mechanism dominates grain evolution. We present the results from the Discs Down Under program, an 8.8 mm continuum Australia Telescope Compact Array (ATCA) survey targeting 15 transition discs with large (> 20 au) cavities, and compare the resulting dust emission to Atacama Large millimetre/sub-millimetre Array (ALMA) observations. Our ATCA observations resolve the inner cavity for 8 of the 14 detected discs. We fit the visibilities and reconstruct 1D radial brightness models for 10 sources with a S/N > 5sigma. We find that, for sources with a resolved cavity in both wavebands, the 8.8 mm and sub-mm brightness distributions peak at the same radius from the star. We suggest that a similar cavity size for 8.8 mm and sub-mm dust grains is due to a dust trap induced by the presence of a companion.

We present infrared photometry of all 36 potential JWST calibrators for which there is archival Spitzer IRAC data. This photometry can then be used to inform stellar models necessary to provide absolute calibration for all JWST instruments. We describe in detail the steps necessary to measure IRAC photometry from archive retrieval to photometric corrections. To validate our photometry we examine the distribution of uncertainties from all detections in all four IRAC channels as well as compare the photometry and its uncertainties to those from models, ALLWISE, and the literature. 75% of our detections have standard deviations per star of all observations within each channel of less than three percent. The median standard deviations are 1.2, 1.3, 1.1, and 1.9% in [3.6] - [8.0] respectively. We find less than 8% standard deviations in differences of our photometry with ALLWISE, and excellent agreement with literature values (less than 3% difference) lending credence to our measured fluxes. JWST is poised to do ground-breaking science, and accurate calibration and cross-calibration with other missions will be part of the underpinnings of that science.

SPIRou is a near-infrared (nIR) spectropolarimeter at the CFHT, covering the YJHK nIR spectral bands ($980-2350\,\mathrm{nm}$). We describe the development and current status of the SPIRou wavelength calibration in order to obtain precise radial velocities (RVs) in the nIR. We make use of a UNe hollow-cathode lamp and a Fabry-P\'erot \'etalon to calibrate the pixel-wavelength correspondence for SPIRou. Different methods are developed for identifying the hollow-cathode lines, for calibrating the wavelength dependence of the Fabry-P\'erot cavity width, and for combining the two calibrators. The hollow-cathode spectra alone do not provide a sufficiently accurate wavelength solution to meet the design requirements of an internal error of $\mathrm{<0.45\,m\,s^{-1}}$, for an overall RV precision of $\mathrm{1\,m\,s^{-1}}$. However, the combination with the Fabry-P\'erot spectra allows for significant improvements, leading to an internal error of $\mathrm{\sim 0.15\,m\,s^{-1}}$. We examine the inter-night stability, intra-night stability, and impact on the stellar RVs of the wavelength solution.

A new era of exploration of the low radio frequency Universe from the Moon will soon be underway with landed payload missions facilitated by NASA's Commercial Lunar Payload Services (CLPS) program. CLPS landers are scheduled to deliver two radio science experiments, ROLSES to the nearside and LuSEE to the farside, beginning in 2021. These instruments would be pathfinders for a 10-km diameter interferometric array, FARSIDE, composed of 128 pairs of dipole antennas proposed to be delivered to the lunar surface later in the decade. ROLSES and LuSEE, operating at frequencies from 100 kHz to a few tens of MHz, will investigate the plasma environment above the lunar surface and measure the fidelity of radio spectra on the surface. Both use electrically-short, spiral-tube deployable antennas and radio spectrometers based upon previous flight models. ROLSES will measure the photoelectron sheath density to better understand the charging of the lunar surface via photoionization and impacts from the solar wind, charged dust, and current anthropogenic radio frequency interference. LuSEE will measure the local magnetic field and exo-ionospheric density, interplanetary radio bursts, Jovian and terrestrial natural radio emission, and the galactic synchrotron spectrum. FARSIDE, and its precursor risk-reduction six antenna-node array PRIME, would be the first radio interferometers on the Moon. FARSIDE would break new ground by imaging radio emission from Coronal Mass Ejections (CME) beyond 2 solar radii, monitor auroral radiation from the B-fields of Uranus and Neptune (not observed since Voyager), and detect radio emission from stellar CMEs and the magnetic fields of nearby potentially habitable exoplanets.

The impact of stellar rotation on the morphology of star cluster colour-magnitude diagrams is widely acknowledged. However, the physics driving the distribution of the equatorial rotation velocities of main-sequence turn-off (MSTO) stars is as yet poorly understood. Using Gaia Data Release 2 photometry and new Southern African Large Telescope medium-resolution spectroscopy, we analyse the intermediate-age ($\sim1\,$Gyr-old) Galactic open clusters NGC 3960, NGC 6134 and IC 4756 and develop a novel method to derive their stellar rotation distributions based on SYCLIST stellar rotation models. Combined with literature data for the open clusters NGC 5822 and NGC 2818, we find a tight correlation between the number ratio of slow rotators and the clusters' binary fractions. The blue-main-sequence stars in at least two of our clusters are more centrally concentrated than their red-main-sequence counterparts. The origin of the equatorial stellar rotation distribution and its evolution remains as yet unidentified. However, the observed correlation in our open cluster sample suggests a binary-driven formation mechanism.

We investigate the possibility of having massive compact stars, $(M>2M_{\odot})$, and at the same time fulfilling the theoretical bounds on the speed of sound $c_s$ in dense matter suggesting that the conformal limit of $c_s^2=1/3$ is approached from below as the density increases. This is possible if two families of stars exist: hadronic stars and quark stars, the latter being entirely composed by quark matter. By using astrophysical data on electromagnetic and gravitational waves signals of a few sources interpreted as quark stars, we show, within a Bayesian analysis framework, that the posterior distribution of $c_s^2$ is peaked around 0.3, and the maximum mass of the most probable equation of state is $\sim 2.13 M_{\odot}$. We finally discuss also the possibility that the maximum mass is larger than $2.6M_{\odot}$ as it could be the case if the secondary component of GW190814 is a compact star and not a black hole.

An origin of Earth life on Mars would resolve significant inconsistencies between the inferred history of life and Earth's geologic history. Life as we know it utilizes amino acids, nucleic acids, and lipids for the metabolic, informational, and compartment-forming subsystems of a cell. Such building blocks may have formed simultaneously from cyanosulfidic chemical precursors in a planetary surface scenario involving ultraviolet light, wet-dry cycling, and volcanism. However, early Earth was a water world, and the timing of the rise of oxygen on Earth is inconsistent with final fixation of the genetic code in response to oxidative stress. A cyanosulfidic origin of life could have taken place on Mars via photoredox chemistry, facilitated by orders of magnitude more sub-aerial crust than early Earth, and an earlier transition to oxidative conditions. Meteoritic bombardment may have generated transient habitable environments and ejected and transferred life to Earth. The Mars 2020 Perseverance Rover offers an unprecedented opportunity to confirm or refute evidence consistent with a cyanosulfidic origin of life on Mars, search for evidence of ancient life, and constrain the evolution of Mars' oxidation state over time. We should seek to prove or refute a Martian origin for life on Earth alongside other possibilities.

Recently, Meneghetti et al. reported an excess of small-scale gravitational lenses in galaxy clusters, compared to simulations of standard cold dark matter (CDM). We propose a self-interacting dark matter (SIDM) scenario, where a population of subhalos in the clusters experiences gravothermal collapse. Using controlled N-body simulations, we show the presence of early-type galaxies in substructures accelerates gravothermal evolution and a collapsed SIDM subhalo has a steeper density profile than its CDM counterpart, leading to a larger radial galaxy-galaxy strong lensing cross section and more lens images, in better agreement with the observations. Our results indicate that strong gravitational lensing can provide a promising test of the self-interacting nature of dark matter.

The BICEP/Keck Array experiment is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background polarization from the South Pole in search of a primordial $B$-mode signature. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We use high-fidelity, in-situ measurements of the beam response to estimate the temperature-to-polarization (T $\rightarrow$ P) leakage in our latest data including observations from 2016 through 2018. This includes three years of BICEP3 observing at 95 GHz, and multifrequency data from Keck Array. Here we present band-averaged far-field beam maps, differential beam mismatch, and residual beam power (after filtering out the leading difference modes via deprojection) for these receivers. We show preliminary results of "beam map simulations," which use these beam maps to observe a simulated temperature (no $Q/U$) sky to estimate T $\rightarrow$ P leakage in our real data.

In the image data collected by astronomical surveys, stars and galaxies often overlap. Deblending is the task of distinguishing and characterizing individual light sources from survey images. We propose StarNet, a fully Bayesian method to deblend sources in astronomical images of crowded star fields. StarNet leverages recent advances in variational inference, including amortized variational distributions and the wake-sleep algorithm. Wake-sleep, which minimizes forward KL divergence, has significant benefits compared to traditional variational inference, which minimizes a reverse KL divergence. In our experiments with SDSS images of the M2 globular cluster, StarNet is substantially more accurate than two competing methods: Probablistic Cataloging (PCAT), a method that uses MCMC for inference, and a software pipeline employed by SDSS for deblending (DAOPHOT). In addition, StarNet is as much as $100,000$ times faster than PCAT, exhibiting the scaling characteristics necessary to perform fully Bayesian inference on modern astronomical surveys.

In $2016-17$, the Galactic transient black hole candidate GRS 1716-249 exhibited an outburst event after a long period of quiescence of almost 23 years. The source remained in the outbursting phase for $\sim 9$ months. We study the spectral and temporal properties of the source during this outburst using archival data from four astronomy satellites, namely MAXI, Swift, NuSTAR and AstroSat. Initial spectral analysis is done using combined disk black body and power-law models. For a better understanding of the accretion flow properties, we studied spectra with the physical two component advective flow (TCAF) model. Accretion flow parameters are extracted directly from the spectral fits with the TCAF model. Low frequency quasi-periodic oscillations are also observed in the Swift/XRT and AstroSat/LAXPC data. From the nature of the variation of the spectral and temporal properties, we find the source remains in hard state during the entire outburst. It never had a transition to other states which makes this event a `failed' outburst. An absence of the softer spectral states is consistent with the class of short orbital period objects, where the source belongs to. From the spectral fit, we also estimate the probable mass of GRS~1716-249 to be in the range of $4.50-5.93 M_\odot$ or $5.02^{+0.91}_{-0.52} M_\odot$.

Spectroscopic observations at extreme and far ultraviolet wavelengths have revealed systematic upflows in the solar transition region and corona. These upflows are best seen in the network structures of the quiet Sun and coronal holes, boundaries of active regions, and dimming regions associated with coronal mass ejections. They have been intensively studied in the past two decades because they are highly likely to be closely related to the formation of the solar wind and heating of the upper solar atmosphere. We present an overview of the characteristics of these upflows, introduce their possible formation mechanisms, and discuss their potential roles in the mass and energy transport in the solar atmosphere. Though past investigations have greatly improved our understanding of these upflows, they have left us with several outstanding questions and unresolved issues that should be addressed in the future. New observations from the Solar Orbiter mission, the Daniel K. Inouye Solar Telescope and the Parker Solar Probe will likely provide critical information to advance our understanding of the generation, propagation and energization of these upflows.

Abridged for arXiv: In this work, we apply a powerful new technique in order to observationally derive accurate assembly histories through a self-consistent combined stellar dynamical and population galaxy model. We present this approach for three edge-on lenticular galaxies from the Fornax3D project -- FCC 153, FCC 170, and FCC 177 -- in order to infer their mass assembly histories individually and in the context of the Fornax cluster. The method was tested on mock data from simulations to quantify its reliability. We find that the galaxies studied here have all been able to form dynamically-cold (intrinsic vertical velocity dispersion $\sigma_z \lesssim 50\ {\rm km}\ {\rm s}^{-1}$) stellar disks after cluster infall. Moreover, the pre-existing (old) high angular momentum components have retained their angular momentum (orbital circularity $\lambda_z > 0.8$) through to the present day. Comparing the derived assembly histories with a comparable galaxy in a low-density environment -- NGC 3115 -- we find evidence for cluster-driven suppression of stellar accretion and merging. We measured the intrinsic stellar age--velocity-dispersion relation and find that the shape of the relation is consistent with galaxies in the literature across redshift. There is tentative evidence for enhancement in the luminosity-weighted intrinsic vertical velocity dispersion due to the cluster environment. But importantly, there is an indication that metallicity may be a key driver of this relation. We finally speculate that the cluster environment is responsible for the S0 morphology of these galaxies via the gradual external perturbations, or `harassment', generated within the cluster.

The energy range between about 100 keV and 1 GeV is of interest for a vast class of astrophysical topics. In particular, (1) it is the missing ingredient for understanding extreme processes in the multi-messenger era; (2) it allows localizing cosmic-ray interactions with background material and radiation in the Universe, and spotting the reprocessing of these particles; (3) last but not least, gamma-ray emission lines trace the formation of elements in the Galaxy and beyond. In addition, studying the still largely unexplored MeV domain of astronomy would provide for a rich observatory science, including the study of compact objects, solar- and Earth-science, as well as fundamental physics. The technological development of silicon microstrip detectors makes it possible now to detect MeV photons in space with high efficiency and low background. During the last decade, a concept of detector ("ASTROGAM") has been proposed to fulfil these goals, based on a silicon hodoscope, a 3D position-sensitive calorimeter, and an anticoincidence detector. In this paper we stress the importance of a medium size (M-class) space mission, dubbed "ASTROMEV", to fulfil these objectives.

The large enhancement of the primordial power spectrum of the curvature perturbation can seed the formation of primordial black hole, that can play as a dark matter component in the Universe. In multi-filed inflation models, the curved trajectory of the scalar fields in the field space can generate a peak in the power spectrum on small scales due to the existence of the isocurvature perturbation. Here we show that a potential can be reconstructed from a given power spectrum, which is made of a scale-invariant one on large scales and the other function with a peak on small scales. In multi-field inflation models the reconstructed potential may not be unique and we can find different potentials from a given power spectrum.

We investigated the radio spectra of two magnetars, PSR J1622$-$4950 and 1E 1547.0$-$5408, using observations from the Australia Telescope Compact Array and the Atacama Large Millimeter/submillimeter Array taken in 2017. Our observations of PSR J1622$-$4950 show a steep spectrum with a spectral index of $-$1.3 $\pm$ 0.2 in the range of 5.5-45 GHz during its re-activating X-ray outburst in 2017. By comparing the data taken at different epochs, we found significant enhancement in the radio flux density. The spectrum of 1E 1547.0$-$5408 was inverted in the range of 43-95 GHz, suggesting a spectral peak at a few hundred gigahertz. Moreover, we obtained the X-ray and radio data of radio magnetars, PSR J1622$-$4950 and SGR J1745$-$2900, from literature and found two interesting properties. First, radio emission is known to be associated with X-ray outburst but has different evolution. We further found that the rising time of the radio emission is much longer than that of the X-ray during the outburst. Second, the radio magnetars may have double peak spectra at a few GHz and a few hundred GHz. This could indicate that the emission mechanism is different in the cm and the sub-mm bands. These two phenomenons could provide a hint to understand the origin of radio emission and its connection with the X-ray properties.

AstroSat/UVIT carries two gratings in the FUV channel and a single grating in the NUV channel. These gratings are useful for low resolution, slitless spectroscopy in the far and near UV bands of a variety of cosmic sources such as hot stars, interacting binaries, active galactic nuclei, etc. We present the calibration of these gratings using observations of UV standards NGC40 and HZ4. We perform wavelength and flux calibration and derive effective areas for different grating orders. We find peak effective areas of 18.7cm^2 at 2325 Angstrom for the -1 order of NUV-Grating, 4.5cm^2 at 1390 Angstrom for the -2 order of FUV-Grating1, and 4.3cm^2 at 1500 Angstrom for the -2 order of FUV-Grating2. The FWHM spectral resolution of the FUV gratings is 14.6 Angstrom in the -2 order. The -1 order of NUV grating has an FWHM resolution of 33 Angstrom. We find excellent agreement in flux measurements between the FUV/NUV gratings and all broadband filters. We have generated spectral responses of the UVIT gratings and broadband filters that can directly be used in the spectral fitting packages such as XSPEC, Sherpa, and ISIS, thus allowing spectral analysis of UVIT data either separately or jointly with X-ray data from AstroSat or other missions.

A Rb deficiency by a factor two with respect to the Sun has been found in M dwarfs of solar metallicity. This deficiency is difficult to understand from both the observational and nucleosynthesis point of views. To test the reliability of this Rb deficiency, we study the Rb and Zr abundances in a sample of KM-type giant stars in a similar metallicity range extracted from the AMBRE Project. We derive Rb and Zr abundances in 54 giant stars with metallicity close to solar by spectral synthesis in LTE and NLTE. The impact of the Zeeman broadening in the RbI line is also studied. The LTE analysis results in a Rb deficiency in giant stars smaller than that obtained in M dwarfs, but the NLTE [Rb/Fe] ratios are very close to solar in the full metallicity range. This contrasts with the figure found in M dwarfs. We investigate the effect of gravitational settling and magnetic activity as possible causes of the Rb deficiency found in M dwarfs. While, the former phenomenon has a negligible impact on the surface Rb abundance, the existence of an average magnetic field with intensity typical of that observed in M dwarfs may result in systematic Rb abundance underestimations if the Zeeman broadening is not considered in the spectral synthesis. The new [Rb,Zr/Fe] vs. [Fe/H] relationships can be explained when the Rb production by rotating massive stars and low-and-intermediate mass stars are considered, without the need of any deviation from the standard s-process nucleosynthesis in AGB stars as previously suggested.

A small fraction of GRBs with available data down to soft X-rays ($\sim0.5$ keV) have been shown to feature a spectral break in the low-energy part ($\sim1-10$ keV) of their prompt emission spectrum. The overall spectral shape is consistent with optically thin synchrotron emission from a population of particles that have cooled on a timescale comparable to the dynamic time to energies that are still much higher than their rest mass energy (marginally fast cooling regime). We consider a hadronic scenario and investigate if the prompt emission of these GRBs can originate from relativistic protons that radiate synchrotron in the marginally fast cooling regime. Using semi-analytical methods, we derive the source parameters, such as magnetic field strength and proton luminosity, and calculate the high-energy neutrino emission expected in this scenario. We also investigate how the emission of secondary pairs produced by photopion interactions and $\gamma\gamma$ pair production affect the broadband photon spectrum. We support our findings with detailed numerical calculations. Strong modification of the photon spectrum below the break energy due to the synchrotron emission of secondary pairs is found, unless the bulk Lorentz factor is very large ($\Gamma \gtrsim 10^3$). Moreover, this scenario predicts unreasonably high Poynting luminosities because of the strong magnetic fields ($10^6-10^7$ G) that are necessary for the incomplete proton cooling. Our results strongly disfavor marginally fast cooling protons as an explanation of the low-energy spectral break in the prompt GRB spectra.

We use Gamma-ray burst (GRB) spectra total continuum absorption to estimate the key intergalactic medium (IGM) properties of hydrogen column density ($\mathit{N}\textsc{hxigm}$), metallicity, temperature and ionisation parameter over a redshift range of $1.6 \leq z \leq 6.3$, using photo-ionisation (PIE) and collisional ionisation equilibrium (CIE) models for the ionised plasma. We use more realistic host metallicity, dust corrected where available, in generating the host absorption model, assuming that the host intrinsic hydrogen column density is equal to the measured ionisation corrected intrinsic neutral column from UV spectra ($\textit{N}\textsc{h}\/\ \textsc{i,ic}$). We find that the IGM property results are similar, regardless of whether the model assumes all PIE or CIE. The $\mathit{N}\textsc{hxigm}$ scales as $(1 + z)^{1.0\/\ -\/\ 1.9}$, with equivalent hydrogen mean density at $z = 0$ of $n_0 = 1.8^{+1.5}_{-1.2} \times 10^{-7}$ cm$^{-3}$. The metallicity ranges from $\sim0.1Z\sun$ at $z \sim 2$ to $\sim0.001Z\sun$ at redshift $z > 4$. The PIE model implies a less rapid decline in average metallicity with redshift compared to CIE. Under CIE, the temperature ranges between $5.0 <$ log$(T/$K$)<\/\ 7.1$. For PIE the ionisation parameter ranges between $0.1 <$ log$(\xi) < 2.9$. Using our model, we conclude that the IGM contributes substantially to the total absorption seen in GRB spectra and that this contribution rises with redshift, explaining why the hydrogen column density inferred from X-rays is substantially in excess of the intrinsic host contribution measured in UV.

We report about observations of the solar U+III bursts on 5 June of 2020 by means of a new active antenna designed to receive radiation in 4-70 MHz. This instrument can serve as a prototype of the ultra-long-wavelength radiotelescope for observations on the farside of the Moon. Our analysis of experimental data is based on simultaneous records obtained with the antenna arrays GURT and NDA in high frequency and time resolution, e-Callisto network as well as by using the space-based observatories STEREO and WIND. The results from this observational study confirm the model of Reid and Kontar (2017).

So far, only transient Gravitational waves (GWs) produced by catastrophic events of extra-galactic origin have been detected. However, it is generally believed that there should be also continuous sources of GWs within our galaxy, such as accreting neutron stars (NSs). In fact, in accreting NSs, centrifugal forces can be so strong to break the neutron star crust (causing a starquake), thus producing a quadrupole moment responsible for the continuous emission of GWs. At equilibrium, the angular momentum gained by accretion and lost via GWs emission should balance each other, stopping the stellar spin-up. We hereinafter investigate the above physical picture within the framework of a Newtonian model describing compressible, non-magnetized, and self-gravitating NSs. In particular, we calculate the rotational frequency need to break the stellar crust of an accreting pulsar and we estimate the upper limit for the ellipticity due to this event. Depending on the equation of state (EoS) and on the mass of the star, we calculated that the starquake-induced ellipticity ranges from $10^{-9}$ to $10^{-5}$. The corresponding equilibrium frequency that we find is in good agreement with observations and, for all the scenarios, it is below the observational limit frequency of $716.36$ Hz. Finally, we also discuss possible observational constraints on the ellipticity upper limit of accreting pulsars.

The oscillation frequencies observed in Sun-like stars are susceptible to being shifted by magnetic activity effects. The measured shifts depend on a complex relationship involving the mode type, the field strength and spatial distribution of activity, as well as the inclination angle of the star. Evidence of these shifts is also present in frequency separation ratios which are often used when inferring global properties of stars in order to avoid surface effects. However, one assumption when using frequency ratios for this purpose is that there are no near-surface perturbations that are non-spherically symmetric. In this work, we studied the impact on inferred stellar properties when using frequency ratios that are influenced by non-homogeneous activity distributions. We generate several sets of artificial oscillation frequencies with various amounts of shift and determine stellar properties using two separate pipelines. We find that for asteroseismic observations of Sun-like targets we can expect magnetic activity to affect mode frequencies which will bias the results from stellar modelling analysis. Although for most stellar properties this offset should be small, typically less than 0.5% in mass, estimates of age and central hydrogen content can have an error of up to 5% and 3% respectively. We expect a larger frequency shift and therefore larger bias for more active stars. We also warn that for stars with very high or low inclination angles, the response of modes to activity is more easily observable in the separation ratios and hence will incur a larger bias.

Stars in the mass range from 8 to 10 solar masses are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes, the observables of these SNe are, therefore, expected to be similar. In this study we explore the differences in these two types of SNe. Specifically, we investigate three different progenitor models: a solar-metallicity ECSN progenitor with an initial mass of 8.8 solar masses, a zero-metallicity progenitor with 9.6 solar masses, and a solar-metallicity progenitor with 9 solar masses, carrying out radiative transfer simulations for these progenitors. We present the resulting light curves for these models. The models exhibit very low photospheric velocity variations of about 2000 km/s, therefore, this may serve as a convenient indicator of low-mass SNe. The ECSN has very unique light curves in broad bands, especially the U band, and does not resemble any currently observed SN. This ECSN progenitor being part of a binary will lose its envelope for which reason the light curve becomes short and undetectable. The SN from the 9.6 solar masses progenitor exhibits also quite an unusual light curve, explained by the absence of metals in the initial composition. The artificially iron polluted 9.6 solar masses model demonstrates light curves closer to normal SNe IIP. The SN from the 9 solar masses progenitor remains the best candidate for so-called low-luminosity SNe IIP like SN 1999br and SN 2005cs.

Galaxies in clusters undergo several phenomena such as ram pressure stripping and tidal interactions, that can trigger or quench their star formation and, in some cases, lead to galaxies acquiring unusual shapes and long tails. We searched for jellyfish galaxy candidates in a sample of 40 clusters from the DAFT/FADA and CLASH surveys covering the redshift range 0.2<z<0.9. In MACS J0717.5+3745 (MACS0717), our large spatial coverage and abundant sampling of spectroscopic redshifts allowed us to pursue a detailed analysis of jellyfish galaxy candidates in this cluster and its extended filament. We looked at the Hubble Space Telescope images of all the cluster galaxies (based on redshifts), and classified them as a function of their likeliness to be jellyfish galaxies, and give catalogues of jellyfish candidates with positions, redshifts, magnitudes, and projected distance to the respective cluster centre. We found 81 jellyfish candidates in the extended region around MACS0717, and 97 in 22 other clusters. Jellyfish galaxy candidates in MACS0717 tend to avoid the densest regions of the cluster, while this does not appear to be the case in the other clusters. For 79 galaxies in MACS0717 and 31 in other clusters, we computed the best stellar population fits with LePhare through the GAZPAR interface. We find that jellyfish candidates tend to be star forming objects, with blue colours, young ages, high star formation rates and specific star formation rates. In a SFR versus stellar mass diagram, jellyfish galaxy candidates appear to have somewhat larger SFRs than non-jellyfish star forming galaxies Based on several arguments, the jellyfish candidates identified in MACS0717 seem to have fallen rather recently into the cluster. A very rough estimate of the proportions of jellyfish galaxies in the studied clusters is about 10%.

We study the decay phase of solar flares in several spectral bands using a method based on that successfully applied to white light flares observed on an M4 dwarf. We selected and processed 102 events detected in the Sun-as-a-star flux obtained with SDO/AIA images in the 1600~{\AA} and 304~{\AA} channels and 54 events detected in the 1700~{\AA} channel. The main criterion for the selection of time profiles was a slow, continuous flux decay without significant new bursts. The obtained averaged time profiles were fitted with analytical templates, using different time intervals, that consisted of a combination of two independent exponents or a broken power law. The average flare profile observed in the 1700~{\AA} channel decayed more slowly than the average flare profile observed on the M4 dwarf. As the 1700~{\AA} emission is associated with a similar temperature to that usually ascribed to M dwarf flares, this implies that the M dwarf flare emission comes from a more dense layer than solar flare emission in the 1700~{\AA} band. The cooling processes in solar flares were best described by the two exponents model, fitted over the intervals t1=[0, 0.5]$t_{1/2}$ and t2=[3, 10]$t_{1/2}$ where $t_{1/2}$ is time taken for the profile to decay to half the maximum value. The broken power law model provided a good fit to the first decay phase, as it was able to account for the impact of chromospheric plasma evaporation, but it did not successfully fit the second decay phase.

Precise and accurate measurements of distances to Galactic X-ray binaries (XRBs) reduce uncertainties in the determination of XRB physical parameters. We have cross-matched the XRB catalogues of Liu et al. (2006, 2007) to the results of Gaia Data Release 2. We identify 86 X-ray binaries with a Gaia candidate counterpart, of which 32 are low-mass X-ray binaries (LMXBs) and 54 are high-mass X-ray binaries (HMXBs). Distances to Gaia candidate counterparts are, on average, consistent with those measured by Hipparcos and radio parallaxes. When compared to distances measured by Gaia candidate counterparts, distances measured using Type I X-ray bursts are systematically larger, suggesting that these bursts reach only 50% of the Eddington limit. However, these results are strongly dependent on the prior assumptions used for estimating distance from the Gaia parallax measurements. Comparing positions of Gaia candidate counterparts for XRBs in our sample to positions of spiral arms in the Milky Way, we find that HMXBs exhibit mild preference for being closer to spiral arms; LMXBs exhibit mild preference for being closer to inter-arm regions. LMXBs do not exhibit any preference for leading or trailing their closest spiral arm. HMXBs exhibit a mild preference for trailing their closest spiral arm. The lack of a strong correlation between HMXBs and spiral arms may be explained by star formation occurring closer to the midpoint of the arms, or a time delay between star formation and HMXB formation manifesting as a spatial separation between HMXBs and the spiral arm where they formed.

We analyse from an observational perspective the formation history and kinematics of a Milky Way-like galaxy from a high-resolution zoom-in cosmological simulation that we compare to those of our Galaxy as seen by Gaia DR2 to better understand the origin and evolution of the Galactic thin and thick discs. The cosmological simulation was carried out with the GADGET-3 TreePM+SPH code using the MUlti Phase Particle Integrator (MUPPI) model. We disentangle the complex overlapping of stellar generations that rises from the top-down and inside-out formation of the galactic disc. We investigate cosmological signatures in the phase-space of mono-age populations and highlight features stemming from past and recent dynamical perturbations. In the simulation, we identify a satellite with a stellar mass of $1.2 \times 10^9$ M$_\odot$, i.e. stellar mass ratio $\Delta \sim 5.5$ per cent at the time, accreted at $z \sim 1.6$, which resembles the major merger Gaia-Sausage-Enceladus that produced the Galactic thick disc, i.e. $\Delta \sim 6$ per cent. We found at $z \sim 0.5-0.4$ two merging satellites with a stellar mass of $8.8 \times 10^8$ M$_\odot$ and $5.1 \times 10^8$ M$_\odot$ that are associated to a strong starburst in the Star Formation History, which appears fairly similar to that recently found in the Solar Neighbourhood. Our findings highlight that detailed studies of coeval stellar populations kinematics, which are made available by current and future Gaia data releases and in synergy with simulations, are fundamental to unravel the formation and evolution of the Milky Way discs.

Some observations of warm carbon chain chemistry (WCCC) cores indicate that they are often located near the edges of molecular clouds. This finding may suggest that WCCC is promoted in star-forming cores exposed to radiation from the interstellar medium. We aim to investigate the chemistry of carbon chains in such a core. A chemical simulation of a gas parcel in a low-mass star-forming core with a full level of irradiation by interstellar photons and cosmic rays was compared to a simulation of a core receiving only one-tenth of such irradiation. In the full irradiation model, the abundances of carbon chains were found to be higher by a factor of few to few hundred, compared to the model with low irradiation. Higher carbon-chain abundances in the prestellar stage and, presumably, in the extended circumstellar envelope, arise because of irradiation of gas and dust by interstellar photons and cosmic rays. A full standard rate of cosmic-ray induced ionization is essential for a high carbon-chain abundance peak to occur in the circumstellar envelope, which is heated by the protostar (the "true" WCCC phenomenon). The full irradiation model has lower abundances of complex organic molecules than the low-irradiation model. We conclude that WCCC can be caused by exposure of a star-forming core to interstellar radiation, or even just to cosmic rays. The Appendix describes an updated accurate approach for calculating the rate of cosmic-ray induced desorption.

The star formation activity of the host galaxies of active galactic nuclei (AGNs) provides valuable insights into the complex interconnections between black hole growth and galaxy evolution. A major obstacle arises from the difficulty of estimating accurate star formation rates in the presence of a strong AGN. Analyzing the $1-500\, \mu m$ spectral energy distributions and high-resolution mid-infrared spectra of low-redshift ($z < 0.5$) Palomar-Green quasars with bolometric luminosity $\sim 10^{44.5}-10^{47.5}\rm\,erg\,s^{-1}$, we find, from comparison with an independent star formation rate indicator based on [Ne II] 12.81$\, \mu m$ and [Ne III] 15.56$\, \mu m$, that the torus-subtracted, total infrared ($8-1000\, \mu m$) emission yields robust star formation rates in the range $\sim 1-250\,M_\odot\,{\rm yr^{-1}}$. Combined with available stellar mass estimates, the vast majority ($\sim 75\%-90\%$) of the quasars lie on or above the main sequence of local star-forming galaxies, including a significant fraction ($\sim 50\%-70\%$) that would qualify as starburst systems. This is further supported by the high star formation efficiencies derived from the gas content inferred from the dust masses. Inspection of high-resolution Hubble Space Telescope images reveals a wide diversity of morphological types, including a number of starbursting hosts that have not experienced significant recent dynamical perturbations. The origin of the high star formation efficiency is unknown.

We study preheating in the Palatini formalism with a quadratic inflaton potential and an added $\alpha R^2$ term. In such models, the oscillating inflaton field repeatedly returns to the plateau of the Einstein frame potential, on which the tachyonic instability fragments the inflaton condensate within less than an e-fold. We find that tachyonic preheating takes place when $\alpha \gtrsim 10^{13}$ and that the energy density of the fragmented field grows with the rate $\Gamma/H \approx 0.011 \times \alpha^{0.31}$. The model extends the family of plateau models with similar preheating behaviour. Although it contains non-canonical quartic kinetic terms in the Einstein frame, we show that, in the first approximation, these can be neglected during both preheating and inflation.

We calculate the finite-temperature r-mode spectrum of a superfluid neutron star accounting for both muons in the core and the entrainment between neutrons and protons. We show that the standard perturbation scheme, considering the rotation rate as an expansion parameter, breaks down in this case. We develop an original perturbation scheme which circumvents this problem by treating both the perturbations due to rotation and (weak) entrainment simultaneously. Applying this scheme, we propose a simple method for calculating the superfluid r-mode eigenfrequency in the limit of vanishing rotation rate. We also calculate the r-mode spectrum at finite rotation rate for realistic microphysics input (adopting, however, the Newtonian framework and Cowling approximation when considering perturbed oscillation equations) and show that the normal r-mode exhibits resonances with superfluid r-modes at certain values of temperatures and rotation frequencies in the parameter range relevant to neutron stars in low-mass X-ray binaries (LMXBs). This turns the recently suggested phenomenological model of resonance r-mode stabilization into a quantitative theory, capable of explaining observations. A strong dependence of resonance rotation rates and temperatures on the neutron superfluidity model allows us to constrain the latter by confronting our calculations with the observations of neutron stars in LMXBs.

We have demonstrated the advantage of combining multi-wavelength observations, from the ultraviolet (UV) to near-infrared, to study Kron 3, a massive star cluster in the Small Magellanic Cloud. We have estimated the radius of the cluster Kron 3 to be 2.'0 and for the first time, we report the identification of NUV-bright red clump (RC) stars and the extension of the RCin colour and magnitude in the NUV vs (NUV-optical) colour-magnitude diagram (CMD). We found that extension of the RC is an intrinsic property of the cluster and it is not due to contamination of field stars or differential reddening across the field. We studied the spectral energy distribution of the RC stars and estimated a small range in temperature ~5000 - 5500K, luminosity ~60 - 90 Land radius ~8.0 - 11.0 Supporting their RC nature. The range of UV magnitudes amongst the RC stars (~23.3 to 24.8 mag) is likely caused by the combined effect of variable mass loss, variation in initial helium abundance (Y_ini=0.23 to 0.28), and a small variation in age (6.5-7.5 Gyr) and metallicity ([Fe/H]=-1.5 to -1.3). Spectroscopic follow-up observations of RC stars in Kron 3 are necessary to confirm the cause of the extended RC.

The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object.

Stellar evolution models require an initial isotopic abundance set as input, but these abundances have not been thoroughly established outside of our sun. Nucleosynthesis studies require reliable isotopic abundances which are challenging to infer from elemental observations and are galaxy specific. Despite the challenges, accurate GCE models for dSphs can provide significant insight on the galactic hierarchy. We present an isotopic history model for the Sculptor dSph galaxy based on astrophysical processes, which is a complementary approach to GCE models. We estimate the isotopic composition of Sculptor's late stage evolution using OMEGA, and we use BBN as the other boundary condition. Each astrophysical process was assigned a parametric function with free parameters according to the underlying physics that dictate their average chemical evolution. The isotopic yields were summed into elemental yields and fit to observational Sculptor abundance data to fix the free parameters. This procedure gives an average isotopic history of Sculptor for massive star, Type 1a SNe, main s-process peak, and r-process contributions, which can be compared with the chemical history of the MW. We find that Type 1a SNe contribute approximately 86 per cent to the late stage evolution Fe abundance, which is greater than typical MW solar values of approximately 70 per cent, and in agreement with other dSph chemical evolution studies. The model also finds that NSMs only contribute approximately 30 per cent to the late stage evolution Eu abundance, further suggesting that CCSNe are the more dominant r-process progenitor in dSphs.

W33A is a well-known example of a high-mass young stellar object showing evidence of a circumstellar disc. We revisited the $K$-band NIFS/Gemini North observations of the W33A protostar using principal components analysis tomography and additional post-processing routines. Our results indicate the presence of a compact rotating disc based on the kinematics of the CO absorption features. The position-velocity diagram shows that the disc exhibits a rotation curve with velocities that rapidly decrease for radii larger than 0\farcs1 ($\sim$250 AU) from the central source, suggesting a structure about four times more compact than previously reported. We derived a dynamical mass of 10.0$^{+4.1}_{-2.2}$ M$_\odot$ for the "disc+protostar" system, about $\sim$33% smaller than previously reported, but still compatible with high-mass protostar status. A relatively compact H$_2$ wind was identified at the base of the large-scale outflow of W33A, with a mean visual extinction of $\sim$63 mag. By taking advantage of supplementary near-infrared maps, we identified at least two other point-like objects driving extended structures in the vicinity of W33A, suggesting that multiple active protostars are located within the cloud. The closest object (Source B) was also identified in the NIFS field of view as a faint point-like object at a projected distance of $\sim$7,000 AU from W33A, powering extended $K$-band continuum emission detected in the same field. Another source (Source C) is driving a bipolar H$_2$ jet aligned perpendicular to the rotation axis of W33A.

Context: The evolution of solar active regions is still not fully understood. The growth and decay of active regions have mostly been studied in case-by-case studies. Aims: Instead of studying the evolution of active regions case by case, we performed a large-scale statistical study to find indications for the statistically most frequent scenario. Methods: We studied a large sample of active regions recorded by the Helioseismic and Magnetic Imager instrument. The sample was split into two groups: forming (367 members) and decaying (679 members) active regions. We tracked individual dark features (i.e. those that are assumed to be intensity counterparts of magnetised fragments from small objects to proper sunspots) and followed their evolution. We investigated the statistically most often locations of fragment merging and splitting as well as their properties. Results: Our results confirm that statistically, sunspots form by merging events of smaller fragments. The coalescence process is driven by turbulent diffusion in a process similar to random-walk, where supergranular flows seem to play an important role. The number of appearing fragments does not seem to significantly correlate with the number of sunspots formed. The formation seems to be consistent with the magnetic field accumulation. Statistically, the merging occurs most often between a large and a much smaller object. The decay of the active region seems to take place preferably by a process similar to the erosion.

Modeling and forecasting the solar wind-driven global magnetic field perturbations is an open challenge. Current approaches depend on simulations of computationally demanding models like the Magnetohydrodynamics (MHD) model or sampling spatially and temporally through sparse ground-based stations (SuperMAG). In this paper, we develop a Deep Learning model that forecasts in Spherical Harmonics space 2, replacing reliance on MHD models and providing global coverage at one minute cadence, improving over the current state-of-the-art which relies on feature engineering. We evaluate the performance in SuperMAG dataset (improved by 14.53%) and MHD simulations (improved by 24.35%). Additionally, we evaluate the extrapolation performance of the spherical harmonics reconstruction based on sparse ground-based stations (SuperMAG), showing that spherical harmonics can reliably reconstruct the global magnetic field as evaluated on MHD simulation.

It has recently been proved that, in the presence of vortex flows, the fluctuation dynamics of a rotating photon-fluid model is governed by the Klein-Gordon equation of an effective massive scalar field in a $(2+1)$-dimensional acoustic black-hole spacetime. Interestingly, it has been demonstrated numerically that the rotating acoustic black hole, like the familiar Kerr black-hole spacetime, may support spatially regular stationary density fluctuations (linearized acoustic scalar `clouds') in its exterior regions. In particular, it has been shown that the composed rotating-acoustic-black-hole-stationary-scalar-field configurations of the photon-fluid model exist in the narrow dimensionless regime $\alpha\equiv\Omega_0/m\Omega_{\text{H}}\in(1,\alpha_{\text{max}})$ with $\alpha_{\text{max}}\simeq 1.08$ [here $\Omega_{\text{H}}$ is the angular velocity of the black-hole horizon and $\{\Omega_0,m\}$ are respectively the effective proper mass and the azimuthal harmonic index of the acoustic scalar field]. In the present paper we use analytical techniques in order to explore the physical and mathematical properties of the acoustic scalar clouds of the photon-fluid model in the regime $\Omega_{\text{H}}r_{\text{H}}\gg1$ of rapidly-spinning central supporting acoustic black holes. In particular, we derive a remarkably compact analytical formula for the discrete resonance spectrum $\{\Omega_0(\Omega_{\text{H}},m;n)\}$ which characterizes the stationary bound-state acoustic scalar clouds of the photon-fluid model. Interestingly, it is proved that the critical (maximal) mass parameter $\alpha_{\text{max}}$, which determines the regime of existence of the composed acoustic-black-hole-stationary-bound-state-massive-scalar-field configurations, is given by the exact dimensionless relation $\alpha_{\text{max}}=\sqrt{{{32}\over{27}}}$.

We perform an observational test of no-hair theorem using quasi-periodic oscillations within the relativistic precession model. Two well motivated metrics we apply are Kerr-Q and Hartle-Thorne metrics in which the quadrupole is the parameter that possibly encodes deviations from the Kerr black hole. The expressions for the quasi-periodic frequencies are derived before comparing the models with the observation. We encounter a degeneracy in constraining spin and quadrupole parameters that makes it difficult to measure their values. In particular, we here propose a novel test of no-hair theorem by adapting the Hartle-Thorne metric. It turns out that a Kerr black hole is a good description of the central object in GRO J1655$-$40 given the present observational precisions.

In this work by using a numerical analysis, we investigate in a quantitative way the late-time dynamics of scalar coupled $f(R,\mathcal{G})$ gravity. Particularly, we consider a Gauss-Bonnet term coupled to the scalar field coupling function $\xi(\phi)$, and we study three types of models, one with $f(R)$ terms that are known to provide a viable late-time phenomenology, and two Einstein-Gauss-Bonnet types of models. Our aim is to write the Friedmann equation in terms of appropriate statefinder quantities frequently used in the literature, and we numerically solve it by using physically motivated initial conditions. In the case that $f(R)$ gravity terms are present, the contribution of the Gauss-Bonnet related terms is minor, as we actually expected. This result is robust against changes in the initial conditions of the scalar field, and the reason is the dominating parts of the $f(R)$ gravity sector at late times. In the Einstein-Gauss-Bonnet type of models, we examine two distinct scenarios, firstly by choosing freely the scalar potential and the scalar Gauss-Bonnet coupling $\xi(\phi)$, in which case the resulting phenomenology is compatible with the latest Planck data and mimics the $\Lambda$-Cold-Dark-Matter model. In the second case, since there is no fundamental particle physics reason for the graviton to change its mass, we assume that primordially the tensor perturbations propagate with the speed equal to that of light's, and thus this constraint restricts the functional form of the scalar coupling function $\xi(\phi)$, which must satisfy the differential equation $\ddot{\xi}=H\dot{\xi}$.

Mutual entrainment effects in cold neutron-proton mixtures are studied in the framework of the self-consistent nuclear energy-density functional theory. Exact expressions for the mass currents, valid for both homogeneous and inhomogeneous systems, are directly derived from the time-dependent Hartree-Fock equations with no further approximation. The equivalence with the Fermi-liquid expression is also demonstrated. Focusing on neutron-star cores, a convenient and simple analytical formulation of the entrainment matrix in terms of the isovector effective mass is found, thus allowing to relate entrainment phenomena in neutron stars to isovector giant dipole resonances in finite nuclei. Results obtained with different functionals are presented. These include the Brussels-Montreal functionals, for which unified equations of state of neutron stars have been recently calculated.

This paper analyzes the experiment presented in 2019 by the Event Horizon Telescope (EHT) Collaboration that revealed the first image of the supermassive black hole at the center of galaxy M87. The very first question asked by the EHT Collaboration is: What is the compact object at the center of galaxy M87? Does it have a horizon? Is it a Kerr black hole? In order to answer these questions, the EHT Collaboration first endorses the working hypothesis that the central object is a black hole described by the Kerr metric, i.e. a spinning Kerr black hole as predicted by classical general relativity. They choose this hypothesis based on previous research and observations of the galaxy M87. After having adopted the Kerr black hole hypothesis, the EHT Collaboration proceeds to test it. They confront this hypothesis with the data collected in the 2017 EHT experiment. They then compare the Kerr rotating black hole hypothesis with alternative explanations and finally find that their hypothesis is consistent with the data. In this paper I describe the complex methods used to test the spinning Kerr black hole hypothesis. I conclude this paper with a discussion of the implications of the findings presented here with respect to Hawking radiation.

A new potential energy surface (PES) and dynamical study are presented of the reactive process between H2CO + OH towards the formation of HCO + H2O and HCOOH + H. In this work a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, what prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many body term has been used fitted with a novel artificial neural networks scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing a good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated to excitations of the normal modes of the ring polymer.

We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, which can be exploited to efficiently produce waveform models that faithfully encode the properties of EMRI systems. We look at examples for a range of different orbital geometries: spherical orbits, equatorial eccentric orbits, and one example of generic (inclined and eccentric) orbits. To our knowledge, this is the first calculation of a generic EMRI waveform that uses strong-field radiation reaction. We examine waveforms in both the time and frequency domains. Although EMRIs evolve slowly enough that the stationary phase approximation (SPA) to the Fourier transform is valid, the SPA calculation must be done to higher order for some voices, since their instantaneous frequency can change from chirping forward ($\dot f > 0$) to chirping backward ($\dot f < 0$). The approach we develop can eventually be extended to more complete EMRI waveform models, for example to include effects neglected by the adiabatic approximation such as the conservative self force and spin-curvature coupling.