Search

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: Submitted to ApJ, 21 figures, 7 tables

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 analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments., Comment: 30 pages, 17 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Comments are warmly welcomed

We present the results of a survey of CO(1-0) emission in 14 infrared luminous dusty star forming galaxies (DSFGs) at 2 < z < 4 with the NSF's Karl G. Jansky Very Large Array. All sources are detected in CO(1-0), with an ~1arcsec angular resolution. Seven sources show extended and complex structure. We measure CO luminosities of $({\mu})L'_{CO(1-0)}=0.4-2.9x10^{11}$ K km s$^{-1}$ pc$^2$, and molecular gas masses of (${\mu}$)M$_{H2}$ = 1.3 - 8.6 x 10$^{11}$ Mo, where ({\mu}) is the magnification factor. The derived molecular gas depletion times of t$_{\rm dep}$ = 40 - 460 Myr, cover the expected range of both normal star forming galaxies and starbursts. Comparing to the higher-J CO transitions previously observed for the same sources, we find CO temperature brightness ratios of r$_{32/10}$ = 0.4 - 1.4, r$_{43/10}$ = 0.4 - 1.7, and r$_{54/10}$ = 0.3 - 1.3. We find a wide range of CO spectral line energy distributions (SLEDs), in agreement with other high-z DSFGs, with the exception of three sources that are most comparable to the Cloverleaf and APM08279+5255. Based on radiative transfer modelling of the CO SLEDs we determine densities of n$_{H2}$ = 0.3 - 8.5 x 10$^3$ cm$^{-3}$ and temperatures of T$_K$ = 100 - 200 K. Lastly, four sources are detected in the continuum, three have radio emission consistent with their infrared derived star formation rates, while HerBS-70E requires an additional synchrotron radiation component from an active galactic nucleus. Overall, we find that even though the sample is similarly luminous in the infrared, by tracing the CO(1-0) emission a diversity of galaxy and excitation properties are revealed, demonstrating the importance of CO(1-0) observations in combination to higher-J transitions., Comment: 25pages , 13 figures, Accepted for publication at ApJ

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 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

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

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

Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date., Comment: ApJ, in press; more information at https://yse.ucsc.edu/

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts., Comment: 20 Pages, 8 Figures + 7 Appendix Figures + 1 animated figure. Resubmitted in order to reflect the version as printed in MNRAS

We report on INTEGRAL observations of the soft $\gamma$-ray repeater SGR 1935+2154 performed between 2020 April 28 and May 3. Several short bursts with fluence of $\sim10^{-7}-10^{-6}$ erg cm$^{-2}$ were detected by the IBIS instrument in the 20-200 keV range. The burst with the hardest spectrum, discovered and localized in real time by the INTEGRAL Burst Alert System, was spatially and temporally coincident with a short and very bright radio burst detected by the CHIME and STARE2 radio telescopes at 400-800 MHz and 1.4 GHz, respectively. Its lightcurve shows three narrow peaks separated by $\sim$29 ms time intervals, superimposed on a broad pulse lasting $\sim$0.6 s. The brightest peak had a delay of 6.5$\pm$1.0 ms with respect to the 1.4 GHz radio pulse (that coincides with the second and brightest component seen at lower frequencies). The burst spectrum, an exponentially cut-off power law with photon index $\Gamma=0.7_{-0.2}^{+0.4}$ and peak energy $E_p=65\pm5$ keV, is harder than those of the bursts usually observed from this and other magnetars. By the analysis of an expanding dust scattering ring seen in X-rays with the {\it Neil Gehrels Swift Observatory} XRT instrument, we derived a distance of 4.4$_{-1.3}^{+2.8}$ kpc for SGR 1935+2154, independent of its possible association with the supernova remnant G57.2+0.8. At this distance, the burst 20-200 keV fluence of $(6.1\pm 0.3)\times10^{-7}$ erg cm$^{-2}$ corresponds to an isotropic emitted energy of $\sim1.4\times10^{39}$ erg. This is the first burst with a radio counterpart observed from a soft $\gamma$-ray repeater and it strongly supports models based on magnetars that have been proposed for extragalactic fast radio bursts., Comment: Accepted for publication on The Astrophysical Journal Letters - revised accepted version

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

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

Thanks to its high orbit and a set of complementary detectors providing continuous coverage of the whole sky, the INTEGRAL satellite has unique capabilities for the identification and study of the electromagnetic radiation associated to gravitational waves signals and, more generally, for multi-messenger astrophysics. Here we briefly review the results obtained during the first two observing runs of the advanced LIGO/Virgo interferometers., Comment: To be published in the Proceedings of the THESEUS Workshop 2017 (http://www.isdc.unige.ch/theseus/workshop2017.html), Journal of the Italian Astronomical Society (Mem.SAIt), Editors L. Amati, E. Bozzo, M. Della Valle, D. Gotz, P. O'Brien

We report the e INTernational Gamma-ray Astrophysics Laboratory (INTEGRAL) detection of the short gamma-ray burst GRB 170817A (discovered by Fermi-GBM) with a signal-to-noise ratio of 4.6, and, for the first time, its association with the gravitational waves (GWs) from binary neutron star (BNS) merging event GW170817 detected by the LIGO and Virgo observatories. The significance of association between the gamma-ray burst observed by INTEGRAL and GW170817 is 3.2 $\sigma$, while the association between the Fermi-GBM and INTEGRAL detections is 4.2 $\sigma$. GRB 170817A was detected by the SPI-ACS instrument about 2 s after the end of the gravitational wave event. We measure a fluence of $(1.4 \pm 0.4 \pm 0.6) \times$10$^{-7}$ erg cm$^{-2})$ (75--2000 keV), where, respectively, the statistical error is given at the 1 $\sigma$ confidence level, and the systematic error corresponds to the uncertainty in the spectral model and instrument response. We also report on the pointed follow-up observations carried out by INTEGRAL, starting 19.5 h after the event, and lasting for 5.4 days. We provide a stringent upper limit on any electromagnetic signal in a very broad energy range, from 3 keV to 8 MeV, constraining the soft gamma-ray afterglow flux to $<7.1\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (80--300 keV). Exploiting the unique capabilities of INTEGRAL, we constrained the gamma-ray line emission from radioactive decays that are expected to be the principal source of the energy behind a kilonova event following a BNS coalescence. Finally, we put a stringent upper limit on any delayed bursting activity, for example from a newly formed magnetar., Comment: ApJ Letters published

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 used data from the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) to set upper-limits on the ${\gamma}$-ray and hard X-ray prompt emission associated with the gravitational wave event GW170104, discovered by the LIGO/Virgo collaboration. The unique omni-directional viewing capability of the instruments on-board INTEGRAL allowed us to examine the full 90% confidence level localization region of the LIGO trigger. Depending on the particular spectral model assumed and the specific position within this region, the upper limits inferred from the INTEGRAL observations range from F${\gamma}$=1.9x10-7 erg cm-2 to F${\gamma}$=10-6 erg cm-2 (75 keV - 2 MeV energy range). This translates into a ratio between the prompt energy released in ${\gamma}$-rays along the direction to the observer and the gravitational wave energy of E${\gamma}$/EGW <2.6x10-5 . Using the INTEGRAL results, we can not confirm the ${\gamma}$-ray proposed counterpart to GW170104 by the AGILE team with the MCAL instrument. The reported flux of the AGILE/MCAL event, E2, is not compatible with the INTEGRAL upper limits within most of the 90% LIGO localization region. There is only a relatively limited portion of the sky where the sensitivity of the INTEGRAL instruments was not optimal and the lowest allowed fluence estimated for E2 would still be compatible with the INTEGRAL results. This region was also observed independently by Fermi/GBM and AstroSAT, from which, as far as we are aware, there are no reports of any significant detection of a prompt high-energy event., Comment: accepted for publication in ApJL

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

During the first observing run of LIGO, two gravitational wave events and one lower-significance trigger (LVT151012) were reported by the LIGO/Virgo collaboration. At the time of LVT151012, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was pointing at a region of the sky coincident with the high localization probability area of the event and thus permitted us to search for its electromagnetic counterpart (both prompt and afterglow emission). The imaging instruments on-board INTEGRAL (IBIS/ISGRI, IBIS/PICsIT, SPI, and the two JEM-X modules) have been exploited to attempt the detection of any electromagnetic emission associated with LVT151012 over 3 decades in energy (from 3 keV to 8 MeV). The omni-directional instruments on-board the satellite, i.e. the SPI-ACS and IBIS monitored the entire LVT151012 localization region at energies above 75 keV. We did not find any significant transient source that was spatially and/or temporally coincident with LVT151012, obtaining tight upper limits on the associated hard X-ray and $\gamma$-ray radiation. For typical spectral models, the upper limits on the fluence of the emission from any 1 s long-lasting counterpart of LVT151012 ranges from $F_{\gamma}=$3.5$\times$10$^{-8}$ erg cm$^{-2}$ (20 - 200 keV) to $F_{\gamma}$=7.1$\times$10$^{-7}$ erg cm$^{-2}$ (75 - 2000 keV), constraining the ratio of the isotropic equivalent energy released in the electromagnetic emission to the total energy of the gravitational waves: $E_{75-2000~keV}/E_{GW}<$4.4$\times$10$^{-5}$. Finally, we provide an exhaustive summary of the capabilities of all instruments on-board INTEGRAL to hunt for $\gamma$-ray counterparts of gravitational wave events, exploiting both serendipitous and pointed follow-up observations. This will serve as a reference for all future searches., Comment: accepted for publication in A&A

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)

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands., Comment: For the main Letter, see arXiv:1602.08492

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams., Comment: For Supplement, see https://arxiv.org/abs/1604.07864

Using observations of the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL), we put upper limits on the gamma-ray and hard X-ray prompt emission associated with the gravitational wave event GW150914, discovered by the LIGO/Virgo collaboration. The omni-directional view of the INTEGRAL/SPI-ACS has allowed us to constrain the fraction of energy emitted in the hard X-ray electromagnetic component for the full high-probability sky region of LIGO trigger. Our upper limits on the hard X-ray fluence at the time of the event range from $F_{\gamma}=2 \times 10^{-8}$ erg cm$^{-2}$ to $F_{\gamma}=10^{-6}$ erg cm$^{-2}$ in the 75 keV - 2 MeV energy range for typical spectral models. Our results constrain the ratio of the energy promptly released in gamma-rays in the direction of the observer to the gravitational wave energy E$_\gamma/$E$_\mathrm{GW}<10^{-6}$. We discuss the implication of gamma-ray limits on the characteristics of the gravitational wave source, based on the available predictions for prompt electromagnetic emission., Comment: accepted for publication in ApJL

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