Search

Fast X-ray transients (FXTs) are extragalactic bursts of X-rays first identified in archival X-ray data, and now routinely discovered by the Einstein Probe in real time, which is continuously surveying the night sky in the soft ($0.5 - 4$ keV) X-ray regime. In this Letter, we report the discovery of the second optical counterpart (AT2024gsa) to an FXT (EP240414a). EP240414a is located at a projected radial separation of 27 kpc from its likely host galaxy at $z = 0.4018 \pm 0.0010$. The optical light curve of AT2024gsa displays three distinct components. The initial decay from our first observation is followed by a re-brightening episode, displaying a rapid rise in luminosity to an absolute magnitude of $M_r \sim -21$ after two rest-frame days. While the early optical luminosity and decline rate is similar to luminous fast blue optical transients, the colour temperature of AT2024gsa is distinctly red and we show that the peak flux is inconsistent with a thermal origin. The third component peaks at $M_i \sim -19$ at $\gtrsim 16$ rest-frame days post-FXT, and is compatible with an emerging supernova. We fit the $riz$-band data with a series of power laws and find that the decaying components are in agreement with gamma-ray burst afterglow models, and that the re-brightening may originate from refreshed shocks. By considering EP240414a in context with all previously reported known-redshift FXT events, we propose that Einstein Probe FXT discoveries may all result from high-redshift gamma-ray bursts, and thus are distinct from the previously discovered lower redshift, lower luminosity population of FXTs., Comment: To be submitted to ApJL

Quasi-periodic Eruptions (QPEs) are luminous bursts of soft X-rays from the nuclei of galaxies, repeating on timescales of hours to weeks. The mechanism behind these rare systems is uncertain, but most theories involve accretion disks around supermassive black holes (SMBHs), undergoing instabilities or interacting with a stellar object in a close orbit. It has been suggested that this disk could be created when the SMBH disrupts a passing star, implying that many QPEs should be preceded by observable tidal disruption events (TDEs). Two known QPE sources show long-term decays in quiescent luminosity consistent with TDEs, and two observed TDEs have exhibited X-ray flares consistent with individual eruptions. TDEs and QPEs also occur preferentially in similar galaxies. However, no confirmed repeating QPEs have been associated with a spectroscopically confirmed TDE or an optical TDE observed at peak brightness. Here we report the detection of nine X-ray QPEs with a mean recurrence time of approximately 48 hours from AT2019qiz, a nearby and extensively studied optically-selected TDE. We detect and model the X-ray, ultraviolet and optical emission from the accretion disk, and show that an orbiting body colliding with this disk provides a plausible explanation for the QPEs.

We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction with hydrogen-free circumstellar material. The second peak is fainter (M$_{r}$ = -17.29 $\pm$ 0.03 mag), and spectroscopically similar to an evolved stripped-envelope SNe, with strong blended forbidden [Ca II] and [O II] features. No other known double-peak SN exhibits a light curve similar to that of SN 2020acct. We find the likelihood of two individual SNe occurring in the same star-forming region within that time to be highly improbable, while an implausibly fine-tuned configuration would be required to produce two SNe from a single binary system. We find that the peculiar properties of SN2020acct match models of pulsational pair instability (PPI), in which the initial peak is produced by collisions of shells of ejected material, shortly followed by a terminal explosion. Pulsations from a star with a 72 M$_{\odot}$ helium core provide an excellent match to the double-peaked light curve. The local galactic environment has a metallicity of 0.4 Z$_{\odot}$, a level where massive single stars are not expected retain enough mass to encounter the PPI. However, late binary mergers or a low-metallicity pocket may allow the required core mass. We measure the rate of SN 2020acct-like events to be $<3.3\times10^{-8}$ Mpc$^{-3}$ yr$^{-1}$ at z = 0.07, or <0.1% of the total core-collapse SN rate., Comment: 31 pages, 14 figures. Submitted to ApJL, comments welcome

We present a comprehensive analysis of the photometric and spectroscopic evolution of SN~2021foa, unique among the class of transitional supernovae for repeatedly changing its spectroscopic appearance from hydrogen-to-helium-to-hydrogen-dominated (IIn-to-Ibn-to-IIn) within 50 days past peak brightness. The spectra exhibit multiple narrow ($\approx$ 300--600~km~s$^{-1}$) absorption lines of hydrogen, helium, calcium and iron together with broad helium emission lines with a full-width-at-half-maximum (FWHM) of $\sim 6000$~km~s$^{-1}$. For a steady, wind-mass loss regime, light curve modeling results in an ejecta mass of $\sim 8$ M$_{\odot}$ and CSM mass below 1 M$_{\odot}$, and an ejecta velocity consistent with the FWHM of the broad helium lines. We obtain a mass-loss rate of $\approx 2$ M$_{\odot} {\rm yr}^{-1}$. This mass-loss rate is three orders of magnitude larger than derived for normal Type II SNe. We estimate that the bulk of the CSM of SN~2021foa must have been expelled within half a year, about 15 years ago. Our analysis suggests that SN~2021foa had a helium rich ejecta which swept up a dense shell of hydrogen rich CSM shortly after explosion. At about 60 days past peak brightness, the photosphere recedes through the dense ejecta-CSM region, occulting much of the red-shifted emission of the hydrogen and helium lines, which results in observed blue-shift ($\sim -3000$~km~s$^{-1}$). Strong mass loss activity prior to explosion, such as those seen in SN~2009ip-like objects and SN~2021foa as precursor emission, are the likely origin of a complex, multiple-shell CSM close to the progenitor star., Comment: Revised. Accepted in ApJ

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

We present LAISS (Lightcurve Anomaly Identification and Similarity Search), an automated pipeline to detect anomalous astrophysical transients in real-time data streams. We deploy our anomaly detection model on the nightly ZTF Alert Stream via the ANTARES broker, identifying a manageable $\sim$1-5 candidates per night for expert vetting and coordinating follow-up observations. Our method leverages statistical light-curve and contextual host-galaxy features within a random forest classifier, tagging transients of rare classes (spectroscopic anomalies), of uncommon host-galaxy environments (contextual anomalies), and of peculiar or interaction-powered phenomena (behavioral anomalies). Moreover, we demonstrate the power of a low-latency ($\sim$ms) approximate similarity search method to find transient analogs with similar light-curve evolution and host-galaxy environments. We use analogs for data-driven discovery, characterization, (re-)classification, and imputation in retrospective and real-time searches. To date we have identified $\sim$50 previously known and previously missed rare transients from real-time and retrospective searches, including but not limited to: SLSNe, TDEs, SNe IIn, SNe IIb, SNe Ia-CSM, SNe Ia-91bg-like, SNe Ib, SNe Ic, SNe Ic-BL, and M31 novae. Lastly, we report the discovery of 325 total transients, all observed between 2018-2021 and absent from public catalogs ($\sim$1% of all ZTF Astronomical Transient reports to the Transient Name Server through 2021). These methods enable a systematic approach to finding the "needle in the haystack" in large-volume data streams. Because of its integration with the ANTARES broker, LAISS is built to detect exciting transients in Rubin data., Comment: 44 pages (68 pages with Appendix), 15 figures, accepted to ApJ

We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($\delta t < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist on a characteristic timescale ($t_{\rm IIn}$) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early-time IIn-like features in addition to 35 "comparison" SNe with no evidence of early-time IIn-like features, all with ultraviolet observations. The total sample consists of 50 unpublished objects with 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both $t_{\rm IIn}$ and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through matching of peak multiband absolute magnitudes, rise times, $t_{\rm IIn}$ and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: $t_{\rm IIn} \approx 3.8[\dot{M}/(0.01 \textrm{M}_{\odot} \textrm{yr}^{-1})]$ days., Comment: 58 pages, 24 figures, submitted to ApJ. Supplementary figures available on Github (https://github.com/wynnjacobson-galan/Flash_Spectra_Sample). Data release following publication

We present a pan-chromatic study of AT2017bcc, a nuclear transient that was discovered in 2017 within the skymap of a reported burst-like gravitational wave candidate, G274296. It was initially classified as a superluminous supernova, and then reclassified as a candidate tidal disruption event. Its optical light curve has since shown ongoing variability with a structure function consistent with that of an active galactic nucleus, however earlier data shows no variability for at least 10 years prior to the outburst in 2017. The spectrum shows complex profiles in the broad Balmer lines: a central component with a broad blue wing, and a boxy component with time-variable blue and red shoulders. The H$\alpha$ emission profile is well modelled using a circular accretion disc component, and a blue-shifted double Gaussian which may indicate a partially obscured outflow. Weak narrow lines, together with the previously flat light curve, suggest that this object represents a dormant galactic nucleus which has recently been re-activated. Our time-series modelling of the Balmer lines suggests that this is connected to a disturbance in the disc morphology, and we speculate this could involve a sudden violent event such as a tidal disruption event involving the central supermassive black hole, though this cannot be confirmed, and given an estimated black hole mass of $\gtrsim10^7-10^8$ M$_\odot$ instabilities in an existing disc may be more likely. Although we find that the redshifts of AT2017bcc ($z=0.13$) and G274296 ($z>0.42$) are inconsistent, this event adds to the growing diversity of both nuclear transients and multi-messenger contaminants., Comment: Submitted to MNRAS

GW190425 is the second of only two binary neutron star (BNS) merger events to be significantly detected by the LIGO-Virgo- Kagra gravitational wave detectors. With a detection only in LIGO Livingston, the skymap containing the source was large and no plausible electromagnetic counterpart was found in real time searching in 2019. Here we summarise our ATLAS and Pan-STARRS wide-field optical coverage of the skymap beginning within 1 hour and 3 hours respectively of the GW190425 merger time. More recently, a potential coincidence between GW190425 and a fast radio burst FRB 190425 has been suggested, given their spatial and temporal coincidence. The smaller sky localisation area of FRB 190425 and its dispersion measure have led to the identification of a likely host galaxy, UGC 10667 at a distance of 141 +/- 10 Mpc. Our optical imaging covered the galaxy 6.0 hrs after GW190425 was detected and 3.5 hrs after the FRB 190425. No optical emission was detected and further imaging at +1.2 and +13.2 days also revealed no emission. If the FRB 190425 and GW190425 association were real, we highlight our limits on kilonova emission from a BNS merger in UGC 10667. The model for producing FRB 190425 from a BNS merger involves a supramassive magnetised neutron star spinning down by dipole emission on the timescale of hours. We show that magnetar enhanced kilonova emission is ruled out by optical upper limits. The lack of detected optical emission from a kilonova in UGC 10667 disfavours, but does not disprove, the FRB-GW link for this source., Comment: Submitted to MNRAS, 20th Sept 2023, 9 pages

We systematically study and classify scotogenic models with a local U(1) gauge symmetry. These models give rise to radiative neutrino masses and a stable dark matter candidate, but avoid the theoretical problems of global and discrete symmetries. We restrict the dark sector particle content to up to four scalar or fermionic SU(2) singlets, doublets or triplets and use theoretical arguments based on anomaly freedom, Lorentz and gauge symmetry to find all possible charge assignments of these particles. The U(1) symmetry can be broken by a new Higgs boson to a residual discrete symmetry, that still stabilizes the dark matter candidate. We list the particle content and charge assignments of all non-equivalent models. Specific examples in our class of models that have been studied previously in the literature are the U(1)$_D$ scotogenic and singlet-triplet scalar models breaking to $Z_2$. We also briefly discuss the new phenomenological aspects of our model arising from the presence of a new massless dark photon or massive $Z'$ boson as well as the additional Higgs boson., Comment: Version accepted for publication in JHEP

We present the discovery and extensive follow-up of a remarkable fast-evolving optical transient, AT2022aedm, detected by the Asteroid Terrestrial impact Last Alert Survey (ATLAS). AT2022aedm exhibited a rise time of $9\pm1$ days in the ATLAS $o$-band, reaching a luminous peak with $M_g\approx-22$ mag. It faded by 2 magnitudes in $g$-band during the next 15 days. These timescales are consistent with other rapidly evolving transients, though the luminosity is extreme. Most surprisingly, the host galaxy is a massive elliptical with negligible current star formation. X-ray and radio observations rule out a relativistic AT2018cow-like explosion. A spectrum in the first few days after explosion showed short-lived He II emission resembling young core-collapse supernovae, but obvious broad supernova features never developed; later spectra showed only a fast-cooling continuum and narrow, blue-shifted absorption lines, possibly arising in a wind with $v\approx2700$ km s$^{-1}$. We identify two further transients in the literature (Dougie in particular, as well as AT2020bot) that share similarities in their luminosities, timescales, colour evolution and largely featureless spectra, and propose that these may constitute a new class of transients: luminous fast-coolers (LFCs). All three events occurred in passive galaxies at offsets of $\sim4-10$ kpc from the nucleus, posing a challenge for progenitor models involving massive stars or massive black holes. The light curves and spectra appear to be consistent with shock breakout emission, though usually this mechanism is associated with core-collapse supernovae. The encounter of a star with a stellar mass black hole may provide a promising alternative explanation., Comment: Accepted in ApJL

We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for $\sim$8 days with respect to first light, at which time Doppler broadened features from the fastest SN ejecta form, suggesting a reduction in CSM density at $r \gtrsim 10^{15}$ cm. The early-time light curve of SN2023ixf shows peak absolute magnitudes (e.g., $M_{u} = -18.6$ mag, $M_{g} = -18.4$ mag) that are $\gtrsim 2$ mag brighter than typical type II supernovae, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light curve and multi-epoch spectral models from the non-LTE radiative transfer code CMFGEN and the radiation-hydrodynamics code HERACLES suggests dense, solar-metallicity, CSM confined to $r = (0.5-1) \times 10^{15}$ cm and a progenitor mass-loss rate of $\dot{M} = 10^{-2}$ M$_{\odot}$yr$^{-1}$. For the assumed progenitor wind velocity of $v_w = 50$ km s$^{-1}$, this corresponds to enhanced mass-loss (i.e., ``super-wind'' phase) during the last $\sim$3-6 years before explosion., Comment: 19 pages, 8 figures. Accepted for publication in ApJL

The formation process of fossil groups (FGs) is still under debate, and large samples of such objects are still missing. The aim of this paper is to increase the sample of known FGs, and to analyse the properties of their brightest group galaxies (BGG) and compare them with a control sample of non-FG BGGs. Based on the Tinker spectroscopic catalogue of haloes and galaxies, we extract 87 FG and 100 non-FG candidates. For all the objects with data available in UNIONS in the u and r bands, and/or in an extra r-band processed to preserve all low surface brightness features (rLSB), we made a 2D photometric fit of the BGG with GALFIT with one or two Sersic components and analysed how the subtraction of intracluster light contribution modifies the BGG properties. From the SDSS spectra available for the BGGs of 65 FGs and 82 non-FGs, we extracted the properties of their stellar populations with Firefly. We also investigated the origin of the emission lines in a nearby FG, NGC 4104, that has an AGN. A single Sersic profile can fit most objects in the u band, while two Sersics are needed in the r and rLSB bands, both for FGs and non-FGs. Non-FG BGGs cover a larger range of Sersic index. FG BGGs follow the Kormendy relation derived for almost one thousand brightest cluster galaxies (BCGs) by Chu et al. (2022) while non-FGs BGGs are mostly located below this relation, suggesting that FG BGGs have evolved similarly to BCGs, while non-FG BGGs have evolved differently. The above properties can be strongly modified by the subtraction of intracluster light contribution. The stellar populations of FG and non-FG BGGs do not differ significantly. Our results suggest FG and non-FG BGGs have had different formation histories, but it is not possible to trace differences in their stellar populations or large scale distributions., Comment: Accepted for publication in Astronomy & Astrophysics

We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-law $f(t) \propto t^{-1.556\pm0.002}$ best describes the light curve. In addition to the high foreground Milky Way dust extinction along this line of sight, the data favour additional extinction to consistently model the optical to X-ray flux with optically thin synchrotron emission. We fit the X-ray-derived power-law to the optical light curve and find good agreement with the measured data up to $5-6$\,days. Thereafter we find a flux excess in the $riy$ bands which peaks in the observer frame at $\sim20$\,days. This excess shares similar light curve profiles to the type Ic broad-lined supernovae SN~2016jca and SN~2017iuk once corrected for the GRB redshift of $z=0.151$ and arbitrarily scaled. This may be representative of a supernova emerging from the declining afterglow. We measure rest-frame absolute peak AB magnitudes of $M_g=-19.8\pm0.6$ and $M_r=-19.4\pm0.3$ and $M_z=-20.1\pm0.3$. If this is an SN component, then Bayesian modelling of the excess flux would imply explosion parameters of $M_{\rm ej}=7.1^{+2.4}_{-1.7}$ M$_{\odot}$, $M_{\rm Ni}=1.0^{+0.6}_{-0.4}$ M$_{\odot}$, and $v_{\rm ej}=33,900^{+5,900}_{-5,700} kms^{-1}$, for the ejecta mass, nickel mass and ejecta velocity respectively, inferring an explosion energy of $E_{\rm kin}\simeq 2.6-9.0\times10^{52}$ ergs., Comment: Submitted to ApJL on 23rd January 2023, for consideration for publication in the special issue on GRB 221009A. Accepted on 2nd March 2023. The results of this paper are under press embargo until 28th March 2023. 15 pages, 7 figures, 1 table

Thanks to the advent of large-scale optical surveys, a diverse set of flares from the nuclear regions of galaxies has recently been discovered. These include the disruption of stars by supermassive black holes at the centers of galaxies - nuclear transients known as tidal disruption events (TDEs). Active galactic nuclei (AGN) can show extreme changes in the brightness and emission line intensities, often referred to as changing-look AGN (CLAGN). Given the physical and observational similarities, the interpretation and distinction of nuclear transients as CLAGN or TDEs remains difficult. One of the obstacles of making progress in the field is the lack of well-sampled data of long-lived nuclear outbursts in AGN. Here, we study PS16dtm, a nuclear transient in a Narrow Line Seyfert 1 (NLSy1) galaxy, which has been proposed to be a TDE candidate. Our aim is to study the spectroscopic and photometric properties of PS16dtm, in order to better understand the outbursts originating in NLSy1 galaxies. Our extensive multiwavelength follow-up that spans around 2000 days includes photometry and spectroscopy in the UV/optical, as well as mid-infrared (MIR) and X-ray observations. Furthermore, we improved an existing semiempirical model in order to reproduce the spectra and study the evolution of the spectral lines. The UV/optical light curve shows a double peak at $\sim50$ and $\sim100$ days after the first detection, and it declines and flattens afterward, reaching preoutburst levels after 2000 days of monitoring. The MIR light curve rises almost simultaneously with the optical, but unlike the UV/optical which is approaching the preoutburst levels in the last epochs of our observations, the MIR emission is still rising at the time of writing. The optical spectra show broad Balmer features and the strongest broad Fe II emission ever detected in a nuclear transient. [abridged], Comment: Accepted for publication in A&A, 19 pages and 18 figures

We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time., Comment: Accepted to ApJS; 64 pages; 35 figures; 10 tables

We present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of 800 km/s; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outflowing slower than a typical Wolf-Rayet wind velocity of $>$1000 km/s. We identify helium in NIR spectra obtained two weeks after maximum and in optical spectra at three weeks, demonstrating that the CSM is not fully devoid of helium. We never detect broad spectral features from SN ejecta, including in spectra extending to the nebular phase, a unique characteristic among SNe~Icn. Compared to other SNe Icn, SN 2022ann has a low luminosity, with a peak o-band absolute magnitude of -17.7, and evolves slowly. We model the bolometric light curve and find it is well-described by 1.7 M_Sun of SN ejecta interacting with 0.2 M_sun of CSM. We place an upper limit of 0.04 M_Sun of Ni56 synthesized in the explosion. The host galaxy is a dwarf galaxy with a stellar mass of 10^7.34 M_Sun (implied metallicity of log(Z/Z_Sun) $\approx$ 0.10) and integrated star-formation rate of log(SFR) = -2.20 M_sun/yr; both lower than 97\% of the galaxies observed to produce core-collapse supernovae, although consistent with star-forming galaxies on the galaxy Main Sequence. The low CSM velocity, nickel and ejecta masses, and likely low-metallicity environment disfavour a single Wolf-Rayet progenitor star. Instead, a binary companion star is likely required to adequately strip the progenitor before explosion and produce a low-velocity outflow. The low CSM velocity may be indicative of the outer Lagrangian points in the stellar binary progenitor, rather than from the escape velocity of a single Wolf-Rayet-like massive star., Comment: 21 pages, 13 figures, to be submitted to MNRAS

Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emitted from accreting gaseous discs close to the BHs. Without this light, they are difficult to detect. Tidal disruption events (TDEs), the luminous flares produced when a star strays close to a BH and is shredded, are a direct way to probe massive BHs. The rise times of these flares theoretically correlate with the BH mass. Here we present AT2020neh, a fast rising TDE candidate, hosted by a dwarf galaxy. AT2020neh can be described by the tidal disruption of a main sequence star by a 10$^{4.7} - 10^{5.9} M_{\odot}$ BH. We find the observable rate of fast rising nuclear transients like AT2020neh to be rare, at $\lesssim 2 \times 10^{-8}$ events Mpc$^{-3}$ yr$^{-1}$. Finding non-accreting BHs in dwarf galaxies is important to determine how prevalent BHs are within these galaxies, and constrain models of BH formation. AT2020neh-like events may provide a galaxy-independent method of measuring IMBH masses., Comment: Accepted for publication in Nature Astronomy

We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns., Comment: By the ENGRAVE collaboration (engrave-eso.org). 35 pages, 20 figures, final version accepted by A&A

We present photometric and spectroscopic data on three extragalactic luminous red novae (LRNe): AT2018bwo, AT2021afy, and AT2021blu. AT2018bwo was discovered in NGC45 (at 6.8 Mpc) a few weeks after the outburst onset. During the monitoring period, the transient reached a peak luminosity of 10^40 erg/s. AT2021afy, hosted by UGC10043 (49.2 Mpc), showed a double-peaked light curve, with the two peaks reaching a similar luminosity of 2.1(+-0.6)x10^41 erg/s. For AT2021blu in UGC5829, (8.6 Mpc), the pre-outburst phase was well-monitored by several photometric surveys, and the object showed a slow luminosity rise before the outburst. The light curve of AT2021blu was sampled with an unprecedented cadence until the object disappeared behind the Sun, and it was then recovered at late phases. The light curve of AT2021blu shows a double peak, with a prominent early maximum reaching a luminosity of 6.5x10^40 erg/s, which is half of that of AT2021afy. The spectra of AT2021afy and AT2021blu display the expected evolution for LRNe: a blue continuum dominated by prominent Balmer lines in emission during the first peak, and a redder continuum consistent with that of a K-type star with narrow absorption metal lines during the second, broad maximum. The spectra of AT2018bwo are markedly different, with a very red continuum dominated by broad molecular features in absorption. As these spectra closely resemble those of LRNe after the second peak, AT2018bwo was probably discovered at the very late evolutionary stages. This would explain its fast evolution and the spectral properties compatible with that of an M-type star. From the analysis of deep frames of the LRN sites years before the outburst, and considerations of the light curves, the quiescent progenitor systems of the three LRNe were likely massive, with primaries ranging from 13Mo for AT2018bwo, to 13-18Mo for AT2021blu, and over 40Mo for AT2021afy., Comment: 33 pages, 19 figures, 7 tables (plus 3 available at the CDS). Accepted for publication in A&A

We present an observational study of the luminous red nova (LRN) AT\,2021biy in the nearby galaxy NGC\,4631. The field of the object was routinely imaged during the pre-eruptive stage by synoptic surveys, but the transient was detected only at a few epochs from $\sim 231$\,days before maximum brightness. The LRN outburst was monitored with unprecedented cadence both photometrically and spectroscopically. AT\,2021biy shows a short-duration blue peak, with a bolometric luminosity of $\sim 1.6 \times 10^{41}$\,erg\,s$^{-1}$, followed by the longest plateau among LRNe to date, with a duration of 210\,days. A late-time hump in the light curve was also observed, possibly produced by a shell-shell collision. AT\,2021biy exhibits the typical spectral evolution of LRNe. Early-time spectra are characterised by a blue continuum and prominent H emission lines. Then, the continuum becomes redder, resembling that of a K-type star with a forest of metal absorption lines during the plateau phase. Finally, late-time spectra show a very red continuum ($T_{\mathrm{BB}} \approx 2050$ K) with molecular features (e.g., TiO) resembling those of M-type stars. Spectropolarimetric analysis indicates that AT\,2021biy has local dust properties similar to those of V838\,Mon in the Milky Way Galaxy. Inspection of archival {\it Hubble Space Telescope} data taken on 2003 August 3 reveals a $\sim 20$\,\msun\ progenitor candidate with log\,$(L/{\rm L}_{\odot}) = 5.0$\,dex and $T_{\rm{eff}} = 5900$\,K at solar metallicity. The above luminosity and colour match those of a luminous yellow supergiant. Most likely, this source is a close binary, with a 17--24\,\msun\ primary component., Comment: 21 pages, 14 figures. Accepted by Astronomy and Astrophysics

We present the results from a dense multi-wavelength (optical/UV, near-infrared (IR), and X-ray) follow-up campaign of the nuclear transient AT2017gge, covering a total of 1698 days from the transient's discovery. The bolometric lightcurve, the black body temperature and radius, the broad H and He I $\lambda$5876 emission lines and their evolution with time, are all consistent with a tidal disruption event (TDE) nature. A soft X-ray flare is detected with a delay of $\sim$200 days with respect to the optical/UV peak and it is rapidly followed by the emergence of a broad He II $\lambda$4686 and by a number of long-lasting high ionization coronal emission lines. This indicate a clear connection between a TDE flare and the appearance of extreme coronal line emission (ECLEs). An IR echo, resulting from dust re-radiation of the optical/UV TDE light is observed after the X-ray flare and the associated near-IR spectra show a transient broad feature in correspondence of the He I $\lambda$10830 and, for the first time in a TDE, a transient high-ionization coronal NIR line (the [Fe XIII] $\lambda$10798) is also detected. The data are well explained by a scenario in which a TDE occurs in a gas and dust rich environment and its optical/UV, soft X-ray, and IR emission have different origins and locations. The optical emission may be produced by stellar debris stream collisions prior to the accretion disk formation, which is instead responsible for the soft X-ray flare, emitted after the end of the circularization process., Comment: Accepted for publication in MNRAS

We use the SMASH survey to obtain unprecedented deep photometry reaching down to the oldest main sequence turn-offs in the colour-magnitude diagrams (CMDs) of the Small Magellanic Cloud (SMC) and quantitatively derive its star formation history (SFH) using CMD fitting techniques. We identify five distinctive peaks of star formation in the last 3.5 Gyr, at $\sim $3, $\sim$2, $\sim$1.1, $\sim $0.45 Gyr ago, and one presently. We compare these to the SFH of the Large Magellanic Cloud (LMC) finding unequivocal synchronicity, with both galaxies displaying similar periods of enhanced star formation over the past $\sim$3.5 Gyr. The parallelism between their SFHs indicates that tidal interactions between the MCs have recurrently played an important role in their evolution for at least the last $\sim$3.5 Gyr, tidally truncating the SMC and shaping the LMC's spiral arm. We show, for the first time, an SMC-LMC correlated SFH at recent times in which enhancements of star formation are localised in the northern spiral arm of the LMC, and globally across the SMC. These novel findings should be used to constrain not only the orbital history of the MCs but also how star formation should be treated in simulations., Comment: 6 pages, 4 figures. Accepted for publication to MNRAS Letters

Using data from a year-long dedicated campaign to observe bright stars, we study the crosstalk channels present in the GPC1 camera. By analyzing these data, we construct a dataset that checks source stars on almost every CCD of every chip within the camera against all possible crosstalk destinations. We use a clustering algorithm to find potential crosstalk occurrences, and then also check all possible combinations (driven by the hardware layout) by eye. This results in a total of 640 rules, with a flux attenuation factor ranging from 2.5x10$^{2}$ for the bright end to 2.5$\times$10$^{4}$ at the faint end. The average value of m$_{cross}$-m$_{src}\approx$-10.25 corresponds to an attenuating factor of 1.25x10$^{4}$, which produces crosstalk ghosts with an average signal-to-noise ratio of 0.64$\pm$0.1 on the bright images. We find no evidence of crosstalk signals between CCDs not connected in the hardware setup. The distribution of attenuation factors is also found to be dependent on crosstalk movement. A clear dependence on cell column offsets is found, consistent with the idea that the source star charge is progressively attenuated during the traversal of cell readout lines. While we can see the trends, the uncertainties on aperture magnitude measurements are large at this stage., Comment: Accepted by PASP

We present multi-wavelength follow-up observations of the ATLAS discovered faint Iax supernova SN 2020kyg that peaked at an absolute magnitude of $M_g \approx -14.9 \pm 0.2$, making it another member of the faint Iax supernova population. The bolometric light curve requires only $\approx 7 \times 10^{-3}$ M$_{\odot}$ of radioactive $^{56}$Ni, with an ejected mass of $M_{\rm ej} \sim 0.4$ M$_{\odot}$ and a low kinetic energy of $E \approx 0.05 \pm 0.02 \times 10^{51}$ erg. We construct a homogeneous volume-limited sample of 902 transients observed by ATLAS within 100 Mpc during a 3.5 year span. Using this sample, we constrain the rates of faint Iax ($M_r \gtrsim -16$) events within 60 Mpc at $12^{+14}_{-8}\%$ of the SN Ia rate. The overall Iax rate, at $15^{+17}_{-9}\%$ of the Ia rate, is dominated by the low-luminosity events, with luminous SNe Iax ($M_r \lesssim -17.5$) like 2002cx and 2005hk accounting for only $0.9^{+1.1}_{-0.5}\%$ of the Ia rate (a 2$\sigma$ upper limit of approximately 3\%). We favour the hybrid CONe WD + He star progenitor channel involving a failed deflagration of a near Chandrasekhar mass white dwarf, expected to leave a bound remnant and a surviving secondary companion, as a candidate explanation for faint Iax explosions. This scenario requires short delay times, consistent with the observed environments of SNe Iax. Furthermore, binary population synthesis calculations have suggested rates of $1-18\%$ of the SN Ia rate for this channel, consistent with our rate estimates., Comment: Accepted to MNRAS after minor revision

Radiative seesaw models have the attractive property of providing dark matter candidates in addition to the generation of neutrino masses. Here we present a study of neutrino signals from the annihilation of dark matter particles that have been gravitationally captured in the Sun in the framework of the scotogenic model. We compute expected event rates in the IceCube detector in its 86-string configuration. As fermionic dark matter does not scatter off nucleons due to its singlet nature and therefore does not accumulate in the Sun, we study the case of scalar dark matter with a scan over the parameter space. Due to a naturally small mass splitting between the two neutral scalar components, inelastic scattering processes with nucleons can occur. We find that for most of the parameter space, i.e. for mass splittings below 500 keV, inelastic scattering in the Sun yields IceCube event rates above 10 events per year, whereas direct detection on Earth is sensitive only to 250 keV. Consequently, a detailed analysis with IceCube could lead to a lower limit on the scalar coupling $\lambda_5\gtrsim1.6\cdot10^{-5}\cdot m_{DM}$/TeV. For larger mass splittings, only elastic scattering occurs in the Sun. In this case, XENON1T limits only allow for models with expected event rates of up to O(0.1) per year. Some of these models, in particular those with large DM mass and fermion coannihilation, could also be tested with a dedicated IceCube analysis of DM annihilation in the Galactic Center., Comment: 24 pages, 9 figures, 2 tables

The DECam Local Volume Exploration survey (DELVE) is a 126-night survey program on the 4-m Blanco Telescope at the Cerro Tololo Inter-American Observatory in Chile. DELVE seeks to understand the characteristics of faint satellite galaxies and other resolved stellar substructures over a range of environments in the Local Volume. DELVE will combine new DECam observations with archival DECam data to cover ~15000 deg$^2$ of high-Galactic-latitude (|b| > 10 deg) southern sky to a 5$\sigma$ depth of g,r,i,z ~ 23.5 mag. In addition, DELVE will cover a region of ~2200 deg$^2$ around the Magellanic Clouds to a depth of g,r,i ~ 24.5 mag and an area of ~135 deg$^2$ around four Magellanic analogs to a depth of g,i ~ 25.5 mag. Here, we present an overview of the DELVE program and progress to date. We also summarize the first DELVE public data release (DELVE DR1), which provides point-source and automatic aperture photometry for ~520 million astronomical sources covering ~5000 deg$^2$ of the southern sky to a 5$\sigma$ point-source depth of g=24.3, r=23.9, i=23.3, and z=22.8 mag. DELVE DR1 is publicly available via the NOIRLab Astro Data Lab science platform., Comment: 30 pages, 14 figures, 6 tables; updated to match published version; public data release at https://datalab.noirlab.edu/delve

Dark matter (DM) scattering and its subsequent capture in the Sun can boost the local relic density, leading to an enhanced neutrino flux from DM annihilations that is in principle detectable at neutrino telescopes. We calculate the event rates expected for a radiative seesaw model containing both scalar triplet and singlet-doublet fermion DM candidates. In the case of scalar DM, the absence of a spin dependent scattering on nuclei results in a low capture rate in the Sun, which is reflected in an event rate of less than one per year in the current IceCube configuration with 86 strings. For singlet-doublet fermion DM, there is a spin dependent scattering process next to the spin independent one, which significantly boosts the event rate and thus makes indirect detection competitive with respect to the direct detection limits imposed by PICO-60. Due to a correlation between both scattering processes, the limits on the spin independent cross section set by XENON1T exclude also parts of the parameter space that can be probed at IceCube. Previously obtained limits by ANTARES, IceCube and Super-Kamiokande from the Sun and the Galactic Center are shown to be much weaker., Comment: 23 pages, 10 figures, 2 tables

We present the results of our monitoring campaigns of the luminous red novae (LRNe) AT 2020hat in NGC 5068 and AT 2020kog in NGC 6106. The two objects were imaged (and detected) before their discovery by routine survey operations. They show a general trend of slow luminosity rise, lasting at least a few months. The subsequent major LRN outbursts were extensively followed in photometry and spectroscopy. The light curves present an initial short-duration peak, followed by a redder plateau phase. AT 2020kog is a moderately luminous event peaking at ~7 x 10^40 erg/s, while AT 2020hat is almost one order of magnitude fainter than AT 2020kog, although it is still more luminous than V838 Mon. In analogy with other LRNe, the spectra of AT 2020kog change significantly with time. They resemble those of type IIn supernovae at early phases, then they become similar to those of K-type stars during the plateau, and to M-type stars at very late phases. In contrast, AT 2020hat already shows a redder continuum at early epochs, and its spectrum shows the late appearance of molecular bands. A moderate-resolution spectrum of AT 2020hat taken at +37 d after maximum shows a forest of narrow P Cygni lines of metals with velocities of 180 km/s, along with an Halpha emission with a full-width at half-maximum velocity of 250 km/s. For AT 2020hat, a robust constraint on its quiescent progenitor is provided by archival images of the Hubble Space Telescope. The progenitor is clearly detected as a mid-K type star, with an absolute magnitude of MF606W = -3.33+-0.09 mag and a colour of F606W-F814W = 1.14+-0.05 mag, which are inconsistent with the expectations from a massive star that could later produce a core-collapse supernova. Although quite peculiar, the two objects nicely match the progenitor versus light curve absolute magnitude correlations discussed in the literature., Comment: 14 pages, 9 figures, 2 tables. Accepted for publication in A&A

With the KATRIN experiment, the determination of the absolute neutrino mass scale down to cosmologically favored values has come into reach. We show that this measurement provides the missing link between the Standard Model and the dark sector in scotogenic models, where the suppression of the neutrino masses is economically explained by their only indirect coupling to the Higgs field. We determine the linear relation between the electron neutrino mass and the scalar coupling $\lambda_5$ associated with the dark neutral scalar mass splitting to be $\lambda_5=3.1\times10^{-9}\ m_{\nu_e}/$eV. This relation then induces correlations among the DM and new scalar masses and their Yukawa couplings. Together, KATRIN and future lepton flavor violation experiments can then probe the fermion DM parameter space, irrespective of the neutrino mass hierarchy and CP phase., Comment: 5 pages, 5 figures, published version

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample., Comment: 28 pages, 23 figures. MNRAS, in press

The GD-1 stream is one of the longest and coldest stellar streams discovered to date, and one of the best objects for constraining the dark matter properties of the Milky Way. Using data from {\it Gaia} DR2 we study the proper motions, distance, morphology and density of the stream to uncover small scale perturbations. The proper motion cleaned data shows a clear distance gradient across the stream, ranging from 7 to 12 kpc. However, unlike earlier studies that found a continuous gradient, we uncover a distance minimum at $\varphi_{1}\approx$-50 deg, after which the distance increases again. We can reliably trace the stream between -85$<\varphi_{1}<$15 deg, showing an even further extent to GD-1 beyond the earlier extension of \citet{Price-Whelan18a}. We constrain the stream track and density using a Boolean matched filter approach and find three large under densities and find significant residuals in the stream track lining up with these gaps. In particular, a gap is visible at $\varphi_{1}$=-3 deg, surrounded by a clear sinusoidal wiggle. We argue that this wiggle is due to a perturbation since it has the wrong orientation to come from a progenitor. We compute a total initial stellar mass of the stream segment of 1.58$\pm$0.07$\times$10$^{4}$ M$_{\odot}$. With the extended view of the spur in this work, we argue that the spur may be unrelated to the adjacent gap in the stream. Finally, we show that an interaction with the Sagittarius dwarf can create features similar to the spur., Comment: Submitted to MNRAS. comments are very welcome

Tidal streams are highly sensitive to perturbations from passing dark matter (DM) subhalos and thus provide a means of measuring their abundance. In a recent paper, we analyzed the distribution of stars along the GD-1 stream with a combination of data from the Gaia satellite and the Pan-STARRS survey, and we demonstrated that the population of DM subhalos predicted by the cold dark matter (CDM) paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. In this paper, we use the measurements of the subhalo mass function (SHMF) from the GD-1 data combined with a similar analysis of the Pal 5 stream to provide novel constraints on alternative DM scenarios that predict a suppression of the SHMF on scales smaller than the mass of dwarf galaxies, marginalizing over uncertainties in the slope and normalization of the unsuppressed SHMF and the susceptibility of DM subhalos in the inner Milky Way to tidal disruption. In particular, we derive a 95% lower limit on the mass of warm dark matter (WDM) thermal relics $m_{\rm WDM}>3.6\,\mathrm{keV}$ from streams alone that strengthens to $m_{\rm WDM}>6.2\,\mathrm{keV}$ when adding dwarf satellite counts. Similarly, we constrain the axion mass in ultra-light ("fuzzy") dark matter (FDM) models to be $m_{\rm FDM}>1.4\times10^{-21}\,\mathrm{eV}$ from streams alone or $m_{\rm FDM}>2.2\times10^{-21}\,\mathrm{eV}$ when adding dwarf satellite counts. Because we make use of simple approximate forms of the streams' SHMF measurement, our analysis is easy to replicate with other alternative DM models that lead to a suppression of the SHMF., Comment: JCAP, in press. Code available to fully replicate the analysis at https://github.com/jobovy/dark-matter-constraints-from-stellar-streams

New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses $\approx 10^{7}$ to $10^9 \ M_{\rm{\odot}}$. We infer a total abundance of dark subhalos normalised to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ ($68 \%$), which corresponds to a mass fraction contained in the subhalos $f_{\rm{sub}} = 0.14 ^{+0.11}_{-0.07} \%$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons., Comment: 20 pages, 14 figures. Also see companion paper : "Novel constraints on the particle nature of dark matter from stellar streams". Accepted in MNRAS. Modified main text slightly. Updated Appendix B and added a plot. Main results and conclusions remain unchanged

We report a new star formation history for the Tucana dwarf spheroidal galaxy, obtained from a new look at a deep HST/ACS colour-magnitude diagram. We combined information from the main sequence turn-off and the horizontal branch to resolve the ancient star formation rates on a finer temporal scale than previously possible. We show that Tucana experienced three major phases of star formation, two very close together at ancient times and the last one ending between 6 and 8 Gyr ago. We show that the three discrete clumps of stars on the horizontal branch are linked to the distinct episodes of star formation in Tucana. The spatial distribution of the clumps reveals that each generation of stars presents a higher concentration than the previous one. The simultaneous modelling of the horizontal branch and the main sequence turn-off also allows us to measure the amount of mass lost by red giant branch stars in Tucana with unprecedented precision, confirming dwarf spheroidals to be excellent laboratories to study the advanced evolution of low-mass stars., Comment: Accepted for publication in A&A. 10 pages, 11 figures

Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data is supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters with a complete density profile covering the full tidal radius (and beyond) for each cluster. We model the density profiles using a set of single-mass models ranging from King and Wilson models to generalised lowered isothermal limepy models and the recently introduced spes models, which allow for the inclusion of potential escapers. We find that both King and Wilson models are too simple to fully reproduce the density profiles, with King (Wilson) models on average underestimating(overestimating) the radial extent of the clusters. The truncation radii derived from the limepy models are similar to estimates for the Jacobi radii based on the cluster masses and their orbits. We show clear correlations between structural and environmental parameters, as a function of Galactocentric radius and integrated luminosity. Notably, the recovered fraction of potential escapers correlates with cluster pericentre radius, luminosity and cluster concentration. The ratio of half mass over Jacobi radius also correlates with both truncation parameter and PE fraction, showing the effect of Roche lobe filling., Comment: 23 pages, 23 figures. Accepted in MNRAS

Using deep wide-field photometry three-year data (Y3) from the Dark Energy Survey (DES), we present a panoramic study of the Fornax dwarf spheroidal galaxy. The data presented here -- a small subset of the full survey -- uniformly covers a region of 25 square degrees centered on the galaxy to a depth of g ~ 23.5. We use this data to study the structural properties of Fornax, overall stellar population, and its member stars in different evolutionary phases. We also search for possible signs of tidal disturbance. Fornax is found to be significantly more spatially extended than what early studies suggested. No statistically significant distortions or signs of tidal disturbances were found down to a surface brightness limit of ~32.1 mag/arcsec^2. However, there are hints of shell-like features located ~30' - 40' from the center of Fornax that may be stellar debris from past merger events. We also find that intermediate age and young main-sequence populations show different orientation at the galaxy center and have many substructures. The deep DES Y3 data allows us to characterize the age of those substructures with great accuracy, both those previously known and those newly discovered in this work, on the basis of their color-magnitude diagram morphology. We find that the youngest overdensities are all found on the Eastern side of Fornax, where the Fornax field population itself is slightly younger than in the West. The high quality DES Y3 data reveals that Fornax has many rich structures, and provides insights into its complex formation history., Comment: 15 pages, 10 figures, 3 tables, submitted to ApJ

We present a new method that incorporates the horizontal branch morphology into synthetic colour-magnitude diagram based star formation history determinations. This method, we call MORGOTH, self-consistently takes into account all the stellar evolution phases up to the early asymptothic giant branch, flexibly modelling red giant branch mass loss. We test MORGOTH on a range of synthetic populations, and find that the inclusion of the horizontal branch significantly increases the precision of the resulting star formation histories. When the main sequence turn-off is detected, MORGOTH can fit the star formation history and the red giant branch mass loss at the same time, efficiently breaking this degeneracy. As part of testing MORGOTH, we also model the observed colour-magnitude diagram of the well studied Sculptor dwarf spheroidal galaxy. We recover a new more detailed star formation history for this galaxy. Both the new star formation history and the red giant branch mass loss we determined for Sculptor with MORGOTH are in good agreement with previous analyses, thus demonstrating the power of this new approach., Comment: 12 pages, 11 figures, accepted for publication in MNRAS

Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos., Comment: submitted to MNRAS, comments welcome

We report the discovery of a thin stellar stream - which we name the Jet stream - crossing the constellations of Hydra and Pyxis. The discovery was made in data from the SLAMS survey, which comprises deep $g$ and $r$ imaging for a $650$ square degree region above the Galactic disc performed by the CTIO Blanco + DECam. SLAMS photometric catalogues will be made publicly available. The stream is approximately 0.18 degrees wide and 10 degrees long, though it is truncated by the survey footprint. Its colour-magnitude diagram is consistent with an old, metal-poor stellar population at a heliocentric distance of approximately 29 kpc. We corroborate this measurement by identifying a spatially coincident overdensity of likely blue horizontal branch stars at the same distance. There is no obvious candidate for a surviving stream progenitor., Comment: MNRAS accepted version

We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $\sigma$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born out of a more turbulent ISM at high redshift and retain that larger velocity dispersion till present day - consistent with recent IFU studies. To ascertain the dominant mechanism and better understand the impact of instabilities and feedback, we develop models based on observed SFHs of these Local Group galaxies in order to create an evolutionary formalism which describes the ISM velocity dispersion due to a galaxy's evolving gas fraction. These empirical models relax the common assumption that the stars are born from gas which has constant velocity dispersion at all redshifts. Using only the observed SFHs as input, the ISM velocity dispersion and a mid-plane scattering model fits the observed AVRs of low mass galaxies without fine tuning. Higher mass galaxies above $M_{vir} > 10^{11}$ need a larger contribution from latent dynamical heating processes (for example minor mergers), in excess of the ISM model. Using the SFHs we also find that supernovae feedback does not appear to be a dominant driver of the gas velocity dispersion compared to gravitational instabilities - at least for dispersions $\sigma \gtrsim 25$ km/s. Together our results point to stars being born with a velocity dispersion close to that of the gas at the time of their formation, with latent dynamical heating operating with a galaxy mass-dependent efficiency. These semi-empirical relations may help constrain the efficiency of feedback and its impact on the physics of disk settling in galaxy formation simulations., Comment: 21 pages, 10 figures. accepted for publication in MNRAS

We present results from deep and very spatially extended CTIO/DECam $g$ and $r$ photometry (reaching out to $\sim$ 2 mag below the oldest MSTO and covering $\sim$ 20 deg$^2$) around the Sextans dSph. We use this data-set to study the structural properties of Sextans overall stellar population and its different stellar evolutionary phases, as well as to search for signs of tidal disturbance from the MW, which would indicate departure from dynamical equilibrium. We perform the most accurate structural analysis to-date of Sextans' stellar components by applying Bayesian MCMC methods to the individual stars' positions. Surface density maps are built by decontaminating the sample through a matched filter analysis of the CMD, and then analysed for departures from axisymmetry. Sextans is found to be considerably less spatially extended than early studies suggested. No significant distortions or tidal disturbances are found down to a surface brightness of $\sim$ 31.8 mag/arcsec$^{-2}$ in V-band. We identify an overdensity in the central regions that may correspond to previously reported kinematic substructure(s). In agreement with previous findings, old $\&$ metal-poor stars such as BHB stars cover a much larger area than stars in other evolutionary phases, and bright BSs are less spatially extended than faint ones. However, the different spatial distribution of bright and faint BSs appears consistent with the general age/metallicity gradients found in Sextans' stellar component. This is compatible with BSs having formed by evolution of binaries and not necessarily due to the presence of a central disrupted globular cluster, as suggested in the literature. We provide structural parameters for the various populations analyzed and make publicly available the photometric catalogue of point-sources as well as a catalogue of literature spectroscopic measurements with updated membership probabilities., Comment: 21 pages, 14 figures, 9 tables, accepted for publication in Astronomy & Astrophysics (A&A), the associated photometric and spectroscopic catalogues will be available at CDS, abstract abridged for arXiv

We present the first detailed study of the behaviour of the stellar proper motion across the entire Galactic Anti-centre area visible in the Sloan Digital Sky Survey data. We use recalibrated SDSS astrometry in combination with positions from {\it Gaia} DR1 to provide tangential motion measurements with a systematic uncertainty $<$5 kms$^{-1}$ for the Main Sequence stars at the distance of the Monoceros Ring. We demonstrate that Monoceros members rotate around the Galaxy with azimuthal speeds of ~230 kms$^{-1}$, only slightly lower than that of the Sun. Additionally, both vertical and azimuthal components of their motion are shown to vary considerably but gradually as a function of Galactic longitude and latitude. The stellar over-density in the Anti-centre region can be split into two components, the narrow, stream-like ACS and the smooth Ring. According to our analysis, these two structures show very similar but clearly distinct kinematic trends, which can be summarised as follows: the amplitude of the velocity variation in $v_{\phi}$ and $v_z$ in the ACS is higher compared to the Ring, whose velocity gradients appear to be flatter. Currently, no model available can explain the entirety of the data in this area of the sky. However, the new accurate kinematic map introduced here should provide strong constraints on the genesis of the Monoceros Ring and the associated sub-structure., Comment: Submitted to MNRAS. comments welcome

We present the spatially resolved star formation history (SFH) of the Carina dwarf spheroidal galaxy, obtained from deep, wide-field g,r imaging and a metallicity distribution from the literature. Our photometry covers $\sim2$ deg$^2$, reaching up to $\sim10$ times the half-light radius of Carina with a completeness higher than $50\%$ at $g\sim24.5$, more than one magnitude fainter than the oldest turnoff. This is the first time a combination of depth and coverage of this quality has been used to derive the SFH of Carina, enabling us to trace its different populations with unprecedented accuracy. We find that Carina's SFH consists of two episodes well separated by a star formation temporal gap. These episodes occurred at old ($>10$ Gyr) and intermediate ($2$-$8$ Gyr) ages. Our measurements show that the old episode comprises the majority of the population, accounting for $54\pm5\%$ of the stellar mass within $1.3$ times the King tidal radius, while the total stellar mass derived for Carina is $1.60\pm0.09\times 10^{6} M_{\rm{\odot}}$, and the stellar mass-to-light ratio $1.8\pm0.2$. The SFH derived is consistent with no recent star formation which hints that the observed blue plume is due to blue stragglers. We conclude that the SFH of Carina evolved independently of the tidal field of the Milky Way, since the frequency and duration of its star formation events do not correlate with its orbital parameters. This result is supported by the age/metallicity relation observed in Carina, and the gradients calculated indicating that outer regions are older and more metal poor., Comment: Accepted in ApJ (22 pages, 13 figures)

We present [C/Fe] and [N/Fe] abundance ratios and CH({\lambda}4300) and S({\lambda}3883) index measurements for 94 red giant branch (RGB) stars in the Sculptor dwarf spheroidal galaxy from VLT/VIMOS MOS observations at a resolving power R= 1150 at 4020 {\AA}. This is the first time that [N/Fe] abundances are derived for a large number of stars in a dwarf spheroidal. We found a trend for the [C/Fe] abundance to decrease with increasing luminosity on the RGB across the whole metallicity range, a phenomenon observed in both field and globular cluster giants, which can be interpreted in the framework of evolutionary mixing of partially processed CNO material. Both our measurements of [C/Fe] and [N/Fe] are in good agreement with the theoretical predictions for stars at similar luminosity and metallicity. We detected a dispersion in the carbon abundance at a given [Fe/H], which cannot be ascribed to measurement uncertainties alone. We interpret this observational evidence as the result of the contribution of different nucleosynthesis sources over time to a not well-mixed interstellar medium. We report the discovery of two new carbon-enhanced, metal-poor stars. These are likely the result of pollution from material enriched by asymptotic giant branch stars, as indicated by our estimates of [Ba/Fe]> +1. We also attempted a search for dissolved globular clusters in the field of the galaxy by looking for the distinctive C-N pattern of second population globular clusters stars in a previously detected, very metal-poor, chemodynamical substructure. We do not detect chemical anomalies among this group of stars. However, small number statistics and limited spatial coverage do not allow us to exclude the hypotheses that this substructure forms part of a tidally shredded globular cluster., Comment: 18 pages, 14 figures, 3 tables. Accepted to A&A

We have determined the detailed star formation history and total mass of the globular clusters in the Fornax dwarf spheroidal using archival HST WFPC2 data. Colour magnitude diagrams are constructed in the F555W and F814W bands and corrected for the effect of Fornax field star contamination, after which we use the routine Talos to derive the quantitative star formation history as a function of age and metallicity. The star formation history of the Fornax globular clusters shows that Fornax 1, 2, 3 and 5 are all dominated by ancient (>10 Gyr) populations. Cluster Fornax 1,2 and 3 display metallicities as low as [Fe/H]=-2.5 while Fornax 5 is slightly more metal-rich at [Fe/H]=-1.8, consistent with resolved and unresolved metallicity tracers. Conversely, Fornax 4 is dominated by a more metal-rich~([Fe/H]=-1.2) and younger population at 10 Gyr, inconsistent with the other clusters. A lack of stellar populations overlapping with the main body of Fornax argues against the nucleus cluster scenario for Fornax 4. The combined stellar mass in globular clusters as derived from the SFH is (9.57$\pm$0.93)$\times$10$^{5}$ M$_{\odot}$ which corresponds to 2.5$\pm$0.2 percent of the total stellar mass in Fornax. The mass of the four most metal-poor clusters can be further compared to the metal-poor Fornax field to yield a mass fraction of 19.6$\pm$3.1 percent. Therefore, the SFH results provide separate supporting evidence for the unusually high mass fraction of the GCs compared to the Fornax field population., Comment: 11 pages, 11 figures. Accepted in A&A