Search

Gamma-ray bursts are the most energetic known explosions. Despite fading rapidly, they allow to measure redshift and important properties of their host-galaxies. We report the photometric and spectroscopic study of GRB 160203A and its host-galaxy. Fine-structure absorption lines, detected in the afterglow at different epochs, allow us to investigate variability due to the strong fading background source. We obtained two optical to near-infrared spectra of the afterglow with X-shooter on ESO/VLT, 18 min and 5.7 hrs after the burst, allowing us to investigate temporal changes of fine-structure absorption lines. We measured HI column density log N(HI/cm-2)=21.75+/-0.10, and several heavy-element ions along the GRB sight-line in the host-galaxy: SiII,AlII,AlIII,CII,NiII,SiIV,CIV,ZnII,FeII, and FeII and SiII fine structure transitions from energetic levels excited by the afterglow, at a redshift z=3.518. We measured [M/H]TOT=-0.78+/-0.13 and [Zn/Fe]FIT=0.69+/-0.15, representing the total(dust-corrected) metallicity and dust depletion, respectively. We detected additional intervening systems along the line of sight at z=1.03,z=1.26,z=1.98,z=1.99,z=2.20 and z=2.83. We could not measure significant variability in the fine-structure lines throughout all the observations and determined an upper limit for the GRB distance from the absorber of d<300 pc, adopting the canonical UV pumping scenario. However, we note that the quality of our data is not sufficient to conclusively rule out collisions as an alternative mechanism. GRB 160203A belongs to a growing sample of GRBs with medium resolution spectroscopy, provided by the Swift/X-shooter legacy program, which enables detailed investigation of the interstellar medium in high-redshift GRB host-galaxies. In particular, this host galaxy shows relatively high metal enrichment and dust depletion already in place when the universe was only 1.8 Gyr old., Comment: 16 pages, 9 figures, 2 appendices, A&A accepted

Long GRB hosts at z<1 are usually low-mass, low metallicity star-forming galaxies. Here we present the until now most detailed, spatially resolved study of the host of GRB 171205A, a grand-design barred spiral galaxy at z=0.036. Our analysis includes MUSE integral field spectroscopy, complemented by high spatial resolution UV/VIS HST imaging and CO(1-0) and HI 21cm data. The GRB is located in a small star-forming region in a spiral arm of the galaxy at a deprojected distance of ~ 8 kpc from the center. The galaxy shows a smooth negative metallicity gradient and the metallicity at the GRB site is half solar, slightly below the mean metallicity at the corresponding distance from the center. Star formation in this galaxy is concentrated in a few HII regions between 5-7 kpc from the center and at the end of the bar, inwards of the GRB region, however, the HII region hosting the GRB is in the top 10% of regions with highest specific star-formation rate. The stellar population at the GRB site has a very young component (< 5 Myr) contributing a significant part of the light. Ionized and molecular gas show only minor deviations at the end of the bar. A parallel study found an asymmetric HI distribution and some additional gas near the position of the GRB, which might explain the star-forming region of the GRB site. Our study shows that long GRBs can occur in many types of star-forming galaxies, however, the actual GRB sites consistently have low metallicity, high star formation and a young population. Furthermore, gas inflow or interactions triggering the star formation producing the GRB progenitor might not be evident in ionized or even molecular gas but only in HI., Comment: 22 pages, 18 figures, Appendix with additional figures, A&A under revision

Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the $\textit{Kilonova Seekers}$ citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. $\textit{Kilonova Seekers}$ launched in July 2023 and received over 600,000 classifications from approximately 2,000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a `gold-standard' training set of 17,682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development., Comment: 20 pages, 15 figures. Accepted in MNRAS

The Epoch of Reionization (EoR) is a key period of cosmological history in which the intergalactic medium (IGM) underwent a major phase change from being neutral to almost completely ionized. Gamma-ray bursts (GRBs) are luminous and unique probes of their environments that can be used to study the timeline for the progression of the EoR. Here we present a detailed analysis of the ESO Very Large Telescope X-shooter spectrum of GRB 210905A, which resides at a redshift of z ~ 6.3. We focus on estimating the fraction of neutral hydrogen, xHI, on the line of sight to the host galaxy of GRB 210905A by fitting the shape of the Lyman-alpha damping wing of the afterglow spectrum. The X-shooter spectrum has a high signal to noise ratio, but the complex velocity structure of the host galaxy limits the precision of our conclusions. The statistically preferred model suggests a low neutral fraction with an 3-sigma upper limit of xHI < 0.15, indicating that the IGM around the GRB host galaxy is mostly ionized. We discuss complications in current analyses and potential avenues for future studies of the progression of the EoR and its evolution with redshift., Comment: 12 pages, 11 figures, submitted to Monthly Notices of the Royal Astronomical Society

Long gamma ray bursts (GRBs) are produced by the collapse of some very massive stars, which emit ultra-relativistic jets. When the jets collide with the interstellar medium they decelerate and generate the so-called afterglow emission, which has been observed to be polarised. In this work we study the polarimetric evolution of GRB 210610B afterglow, at $z = 1.1341$. This allows to evaluate the role of geometric and/or magnetic mechanisms in the GRB afterglow polarisation. We observed GRB 210610B using imaging polarimetry with CAFOS on the 2.2 m Calar Alto Telescope and FORS2 on the 4 $\times$ 8.1 m Very Large Telescope. Complementary optical spectroscopy was obtained with OSIRIS on the 10.4 m Gran Telescopio Canarias. We study the GRB light-curve from X-rays to optical bands and the Spectral Energy Distribution (SED). This allows us to strongly constrain the line-of-sight extinction. Finally, we study the GRB host galaxy using optical/NIR data to fit the SED and derive its integrated properties. GRB 210610B had a bright afterglow with a negligible line-of-sight extinction. Polarimetry was obtained at three epochs: during an early plateau phase, at the time when the light curve breaks, and after the light curve steepened. We observe an initial polarisation of $\sim 4\%$ that goes to zero at the time of the break, and then increases again to $\sim 2\%$ with a change of the position angle of $54 \pm 9$ deg. The spectrum show features with very low equivalent widths, indicating a small amount of material in the line-of-sight within the host. The lack of dust and the low amount of material on the line-of-sight to GRB 210610B allow us to study the intrinsic polarisation of the GRB optical afterglow. We find the GRB polarisation signals are consistent with ordered magnetic fields in refreshed shock or/and hydrodynamics-scale turbulent fields in the forward shock., Comment: 17 pages, 7 figures, 5 tables. Accepted for publication in Astronomy and Astrophysics 23 Feb 2024

In this paper, we present rapid follow-up observations of the short GRB 201006A, consistent with being a compact binary merger, using the LOw Frequency ARray (LOFAR). We have detected a candidate 5.6$\sigma$, short, coherent radio flash at 144 MHz at 76.6 mins post-GRB with a 3$\sigma$ duration of 38 seconds. This radio flash is 27 arcsec offset from the GRB location, which has a probability of occurring by chance of $\sim$0.05% (3.8$\sigma$) when accounting for measurement uncertainties. Despite the offset, we show that the probability of finding an unrelated transient within 40 arcsec of the GRB location is $<10^{-6}$ and conclude that this is a candidate radio counterpart to GRB 201006A. We performed image plane dedispersion and the radio flash is tentatively (2.4$\sigma$) shown to be highly dispersed, allowing a distance estimate, corresponding to a redshift of $0.58\pm0.06$. The corresponding luminosity of the event at this distance is $6.7^{+6.6}_{-4.4} \times 10^{32}$ erg s$^{-1}$ Hz$^{-1}$. If associated with GRB 201006A, this emission would indicate prolonged activity from the central engine that is consistent with being a newborn, supramassive, likely highly magnetised, millisecond spin neutron star (a magnetar)., Comment: MNRAS Accepted

Much of what is known of the chemical composition of the universe is based on emission line spectra from star forming galaxies. Emission-based inferences are, nevertheless, model-dependent and they are dominated by light from luminous star forming regions. An alternative and sensitive probe of the metallicity of galaxies is through absorption lines imprinted on the luminous afterglow spectra of long gamma ray bursts (GRBs) from neutral material within their host galaxy. We present results from a JWST/NIRSpec programme to investigate for the first time the relation between the metallicity of neutral gas probed in absorption by GRB afterglows and the metallicity of the star forming regions for the same host galaxy sample. Using an initial sample of eight GRB host galaxies at z=2.1-4.7, we find a tight relation between absorption and emission line metallicities when using the recently proposed $\hat{R}$ metallicity diagnostic (+/-0.2dex). This agreement implies a relatively chemically-homogeneous multi-phase interstellar medium, and indicates that absorption and emission line probes can be directly compared. However, the relation is less clear when using other diagnostics, such as R23 and R3. We also find possible evidence of an elevated N/O ratio in the host galaxy of GRB090323 at z=3.58, consistent with what has been seen in other $z>4$ galaxies. Ultimate confirmation of an enhanced N/O ratio and of the relation between absorption and emission line metallicities will require a more direct determination of the emission line metallicity via the detection of temperature-sensitive auroral lines in our GRB host galaxy sample., Comment: Accepted for publication in MNRAS; 24 pages, 15 figures

The composition of relativistic gamma-ray burst (GRB) jets and their emission mechanisms are still debated, and they could be matter or magnetically dominated. One way to distinguish these mechanisms arises because a Poynting flux dominated jet may produce low-frequency radio emission during the energetic prompt phase, through magnetic reconnection at the shock front. We present a search for radio emission coincident with three GRB X-ray flares with the LOw Frequency ARray (LOFAR), in a rapid response mode follow-up of long GRB 210112A (at z~2) with a 2 hour duration, where our observations began 511 seconds after the initial swift-BAT trigger. Using timesliced imaging at 120-168 MHz, we obtain upper limits at 3 sigma confidence of 42 mJy averaging over 320 second snapshot images, and 87 mJy averaging over 60 second snapshot images. LOFAR's fast response time means that all three potential radio counterparts to X-ray flares are observable after accounting for dispersion at the estimated source redshift. Furthermore, the radio pulse in the magnetic wind model was expected to be detectable at our observing frequency and flux density limits which allows us to disfavour a region of parameter space for this GRB. However, we note that stricter constraints on redshift and the fraction of energy in the magnetic field are required to further test jet characteristics across the GRB population., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

The chemical enrichment of dust and metals in the interstellar medium (ISM) of galaxies throughout cosmic time is one of the key driving processes of galaxy evolution. Here we study the evolution of the gas-phase metallicities, dust-to-gas (DTG), and dust-to-metal (DTM) ratios of 36 star-forming galaxies at $1.7 < z < 6.3$ probed by gamma-ray bursts (GRBs). We compile all GRB-selected galaxies with intermediate (R=7000) to high (R>40,000) resolution spectroscopic data for which at least one refractory (e.g. Fe) and one volatile (e.g. S or Zn) element have been detected at S/N>3. This is to ensure that accurate abundances and dust depletion patterns can be obtained. We first derive the redshift evolution of the dust-corrected, absorption-line based gas-phase metallicity [M/H]$_{\rm tot}$ in these galaxies, for which we determine a linear relation with redshift ${\rm [M/H]_{tot}}(z) = (-0.21\pm 0.04)z -(0.47\pm 0.14)$. We then examine the DTG and DTM ratios as a function of redshift and through three orders of magnitude in metallicity, quantifying the relative dust abundance both through the direct line-of-sight visual extinction $A_V$ and the derived depletion level. We use a novel method to derive the DTG and DTM mass ratios for each GRB sightline, summing up the mass of all the depleted elements in the dust-phase. We find that the DTG and DTM mass ratios are both strongly correlated with the gas-phase metallicity and show a mild evolution with redshift as well. While these results are subject to a variety of caveats related to the physical environments and the narrow pencil-beam sightlines through the ISM probed by the GRBs, they provide strong implications for studies of dust masses to infer the gas and metal content of high-redshift galaxies, and particularly demonstrate the large offset from the average Galactic value in the low-metallicity, high-redshift regime., Comment: Accepted in A&A

Luminous Fast Blue Optical Transients (LFBOTs) - the prototypical example being AT2018cow - are a rare class of events whose origins are poorly understood. They are characterised by rapid evolution, featureless blue spectra at early times, and luminous X-ray and radio emission. LFBOTs thus far have been found exclusively at small projected offsets from star-forming host galaxies. We present Hubble Space Telescope, Gemini, Chandra and Very Large Array observations of a new LFBOT, AT2023fhn. The Hubble Space Telescope data reveal a large offset (greater than 3.5 half-light radii) from the two closest galaxies, both at a redshift of 0.24. The location of AT2023fhn is in stark contrast with previous events, and demonstrates that LFBOTs can occur in a range of galactic environments., Comment: Accepted for publication in MNRASL. 7 pages, 4 figures, 2 tables

The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia, it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique insight into supernova progenitor systems. In this paper, we employ binary population synthesis to produce kinematic and photometric predictions for ejected secondary stars. We demonstrate that the unbound neutron star velocity distribution from supernovae in binaries closely traces the input kicks. Therefore, the observed distribution of neutron star velocities should be representative of their natal kicks. We evaluate the probability for any given filter, magnitude limit, minimum measurable proper motion (as a function of magnitude), temporal baseline, distance and extinction that an unbound companion can be associated with a remnant. We compare our predictions with results from previous companion searches, and demonstrate that the current sample of stars ejected by the supernova of their companion can be increased by a factor of 5-10 with Gaia data release 3. Further progress in this area is achievable by leveraging the absolute astrometric precision of Gaia, and by obtaining multiple epochs of deep, high resolution near-infrared imaging with the Hubble Space Telescope, JWST and next-generation wide-field near-infrared observatories such as Euclid or the Nancy Grace Roman Space Telescope., Comment: Accepted for publication in MNRAS. 19 pages, 17 figures

The astronomical transient AT2018cow is the closest example of the new class of luminous, fast blue optical transients (FBOTs). Liverpool Telescope RINGO3 observations of AT2018cow are reported here, which constitute the earliest polarimetric observations of an FBOT. At 5.7 days post-explosion, the optical emission of AT2018cow exhibited a chromatic polarization spike that reached ~7% at red wavelengths. This is the highest intrinsic polarization recorded for a non-relativistic explosive transient, and is observed in multiple bands and at multiple epochs over the first night of observations, before rapidly declining. The apparent wavelength dependence of the polarization may arise through depolarization or dilution of the polarized flux, due to conditions in AT~2018cow at early times. A second ``bump" in the polarization is observed at blue wavelengths at ~12 days. Such a high polarization requires an extremely aspherical geometry that is only apparent for a brief period (<1 day), such as shock breakout through an optically thick disk. For a disk-like configuration, the ratio of the thickness to radial extent must be ~10%., Comment: MNRAS Accepted, 10 pages, 8 figures

Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium., Comment: 9 pages, 4 figures, submitted to Astronomy & Astrophysics

We present JWST and Hubble Space Telescope (HST) observations of the afterglow of GRB 221009A, the brightest gamma-ray burst (GRB) ever observed. This includes the first mid-IR spectra of any GRB, obtained with JWST/NIRSPEC (0.6-5.5 micron) and MIRI (5-12 micron), 12 days after the burst. Assuming that the intrinsic spectral slope is a single power-law, with $F_{\nu} \propto \nu^{-\beta}$, we obtain $\beta \approx 0.35$, modified by substantial dust extinction with $A_V = 4.9$. This suggests extinction above the notional Galactic value, possibly due to patchy extinction within the Milky Way or dust in the GRB host galaxy. It further implies that the X-ray and optical/IR regimes are not on the same segment of the synchrotron spectrum of the afterglow. If the cooling break lies between the X-ray and optical/IR, then the temporal decay rates would only match a post jet-break model, with electron index $p<2$, and with the jet expanding into a uniform ISM medium. The shape of the JWST spectrum is near-identical in the optical/nIR to X-shooter spectroscopy obtained at 0.5 days and to later time observations with HST. The lack of spectral evolution suggests that any accompanying supernova (SN) is either substantially fainter or bluer than SN 1998bw, the proto-type GRB-SN. Our HST observations also reveal a disc-like host galaxy, viewed close to edge-on, that further complicates the isolation of any supernova component. The host galaxy appears rather typical amongst long-GRB hosts and suggests that the extreme properties of GRB 221009A are not directly tied to its galaxy-scale environment., Comment: Accepted for publication to the Astrophysical Journal Letters for the GRB 221009A Special Issue. The results of this paper are under press embargo until March 28, 18 UT. 19 pages, 8 figures, 2 tables

The study of the properties of galaxies in the first billion years after the Big Bang is one of the major topic of current astrophysics. Optical/near-infrared spectroscopy of the afterglows of long Gamma-ray bursts (GRBs) provide a powerful diagnostic tool to probe the interstellar medium (ISM) of their host galaxies and foreground absorbers, even up to the highest redshifts. We analyze the VLT/X-shooter afterglow spectrum of GRB 210905A, triggered by the Swift Neil Gehrels Observatory, and detect neutral-hydrogen, low-ionization, high-ionization, and fine-structure absorption lines from a complex system at z=6.3118, that we associate with the GRB host galaxy. We study the ISM properties of the host system, revealing the metallicity, kinematics and chemical abundance pattern. The total metallicity of the z~6.3 system is [M/H]=-1.72+/-0.13, after correcting for dust-depletion and taking into account alpha-element enhancement. In addition, we determine the overall amount of dust and dust-to-metal mass ratio (DTM) ([Zn/Fe]_fit=0.33+/-0.09, DTM=0.18+/-0.03). We find indications of nucleosynthesis due to massive stars and evidence of peculiar over-abundance of aluminium. From the analysis of fine-structure lines, we determine distances of several kpc for the low-ionization gas clouds closest to the GRB. Those farther distances are possibly due to the high number of ionizing photons. Using the HST/F140W image of the GRB field, we show the GRB host galaxy as well as multiple objects within 2" from the GRB. We discuss the galaxy structure and kinematics that could explain our observations, also taking into account a tentative detection of Lyman-alpha emission. Deep spectroscopic observations with VLT/MUSE and JWST will offer the unique possibility of combining our results with the ionized-gas properties, with the goal of better understanding how galaxies in the reionization era form and evolve., Comment: Accepted Publication (In Press on A&A) - 22 pages, 10 figures, 6 tables - Appendix: 6 figures, 3 tables

Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, the deep-embedded self-organizing map (DESOM) to this "real-bogus" classification problem. DESOM combines an autoencoder and a self-organizing map to perform clustering in order to distinguish between real and bogus detections, based on their dimensionality-reduced representations. We use 32x32 normalized detection thumbnails as the input of DESOM. We demonstrate different model training approaches, and find that our best DESOM classifier shows a missed detection rate of 6.6% with a false positive rate of 1.5%. DESOM offers a more nuanced way to fine-tune the decision boundary identifying likely real detections when used in combination with other types of classifiers, for example built on neural networks or decision trees. We also discuss other potential usages of DESOM and its limitations.

Among the large variety of astrophysical sources that we can observe, gamma-ray bursts (GRBs) are the most energetic of the whole Universe. The definition of a general picture describing the physics behind GRBs has always been a compelling task, but the results obtained so far from observations have revealed a puzzling landscape. The lack of a clear, unique paradigm calls for further observations and additional, independent techniques for this purpose. Polarimetry constitutes a very useful example as it allows us to investigate some features of the source such as the geometry of the emitting region and the magnetic field configuration. To date, only a handful of bursts detected by space telescopes have been accompanied by ground-based spectro-polarimetric follow-up, and therefore such an analysis of more GRBs is of crucial importance in order to increase the sample of bursts with multi-epoch polarisation analysis. In this work, we present the analysis of the GRB 080928 optical afterglow, with observations performed with the ESO-VLT FORS1 instrument. We find that the GRB optical afterglow was not significantly polarised on the first observing night. The polarisation degree ($P$) grew on the following night to a level of $P \sim$ 4.5%, giving evidence of polarised radiation at a 4 $\sigma$ confidence level. The GRB 080928 light curve is not fully consistent with standard afterglow models, making any comparison with polarimetric models partly inconclusive. The most conservative interpretation is that the GRB emission was characterised by a homogeneous jet and was observed at an angle of 0.6 $< \theta_{obs}/\theta_{jet} <$ 0.8. Moreover, the non-zero polarisation degree on the second night suggests the presence of a dominant locally ordered magnetic field in the emitting region., Comment: 10 pages, 7 figures. Accepted for publication in A&A

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

Large constellations of orbiting communication satellites will become an important source of noise for present and future astronomical observatories. Mitigation measures rely on high quality predictive models of the position and expected brightness of these objects. Optical linear imaging polarimetry holds promise as a quantitative tool to improve our understanding of the physics of reflection of sunlight off satellite components and through which models of expected brightness can be improved. We present the first simultaneous short-timescale linear polarimetry and optical photometry observations of a geostationary satellite, using the new MOPTOP imaging polarimeter on the 2m Liverpool Telescope. Our target, telecommunication satellite Thor-6, shows prominent short timescale glint-like features in the lightcurve, some as short as seconds. Our polarimetric observations overlap with several of these micro-glints, and have the cadence required to resolve them. We find that the polarisation lightcurve is remarkably smooth, the short time scale glints are not seen to produce strong polarimetric features in our observation. We show how short timescale polarimetry can further constrain the properties of the components responsible for these micro-glints., Comment: 6 figures, 30 pages; accepted for publication in Advances in Space Research

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a 'CSM population synthesis' with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind and ISM densities are compared with the largest sample of afterglow-derived parameters yet compiled, which we make available for the community. We find that high ISM densities around 1000/cm3 best reproduce observations. If long GRBs instead occur in typical ISM densities of approximately 1/cm3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius., Comment: 22 pages, 15 figures, 4 tables. Accepted for publication in MNRAS

It is well established that magnetars are neutron stars with extreme magnetic fields and young ages, but the evolutionary pathways to their creation are still uncertain. Since most massive stars are in binaries, if magnetars are a frequent result of core-collapse supernovae, some fraction are expected to have a bound companion at the time of observation. In this paper, we utilise literature constraints, including deep Hubble Space Telescope imaging, to search for bound stellar companions to magnetars. The magnitude and colour measurements are interpreted in the context of binary population synthesis predictions. We find two candidates for stellar companions associated with CXOU J171405.7-381031 and SGR 0755-2933, based on their J-H colours and H-band absolute magnitudes. Overall, the proportion of the Galactic magnetar population with a plausibly stellar near-infrared counterpart candidate, based on their magnitudes and colours, is between 5 and 10 per cent. This is consistent with a population synthesis prediction of 5 per cent, for the fraction of core-collapse neutron stars arising from primaries which remain bound to their companion after the supernova. These results are therefore consistent with magnetars being drawn in an unbiased way from the natal core-collapse neutron star population, but some contribution from alternative progenitor channels cannot be ruled out., Comment: 14 pages, 6 figures, 3 tables, accepted for publication in MNRAS

We report the discovery of six new magnetar counterpart candidates from deep near-infrared Hubble Space Telescope imaging. The new candidates are among a sample of nineteen magnetars for which we present HST data obtained between 2018-2020. We confirm the variability of previously established near-infrared counterparts, and newly identify candidates for PSRJ1622-4950, SwiftJ1822.3-1606, CXOUJ171405.7-381031, SwiftJ1833-0832, SwiftJ1834.9-0846 and AXJ1818.8-1559 based on their proximity to X-ray localisations. The new candidates are compared with the existing counterpart population in terms of their colours, magnitudes, and near-infrared to X-ray spectral indices. We find two candidates for AXJ1818.8-1559 which are both consistent with previously established counterparts. The other new candidates are likely to be chance alignments, or otherwise have a different origin for their near-infrared emission not previously seen in magnetar counterparts. Further observations and studies of these candidates are needed to firmly establish their nature., Comment: 11 pages, 5 figures, 3 tables, accepted for publication in MNRAS

We present a detailed follow-up of the very energetic GRB 210905A at a high redshift of z = 6.312 and its luminous X-ray and optical afterglow. We obtained a photometric and spectroscopic follow-up in the optical and near-infrared (NIR), covering both the prompt and afterglow emission from a few minutes up to 20 Ms after burst. With an isotropic gamma-ray energy release of Eiso = 1.27E54 erg, GRB 210905A lies in the top ~7% of gamma-ray bursts (GRBs) in terms of energy released. Its afterglow is among the most luminous ever observed. It starts with a shallow evolution that can be explained by energy injection, and it is followed by a steeper decay, while the spectral energy distribution is in agreement with slow cooling in a constant-density environment within the standard fireball theory. A jet break at ~ 46.2+-16.3 d (~6.3 d rest-frame) has been observed in the X-ray light curve; however, it is hidden in the H band due to the contribution from the likely host galaxy, the fourth GRB host at z > 6 known to date. We derived a half-opening angle of 8.4+-1.0 degrees, which is the highest ever measured for a z>6 burst, but within the range covered by closer events. The resulting collimation-corrected gamma-ray energy release of 1E52 erg is also among the highest ever measured. The moderately large half-opening angle argues against recent claims of an inverse dependence of the half-opening angle on the redshift. The total jet energy is likely too large to be sustained by a standard magnetar, and it suggests that the central engine of this burst was a newly formed black hole. Despite the outstanding energetics and luminosity of both GRB 210905A and its afterglow, we demonstrate that they are consistent with those of less distant bursts, indicating that the powering mechanisms and progenitors do not evolve significantly with redshift., Comment: 19 pages, 12 figures, 5 tables, accepted in Astronomy & Astrophysics

We present deep Hubble Space Telescope near-infrared (NIR) observations of the magnetar SGR 1935+2154 from June 2021, approximately 6 years after the first HST observations, a year after the discovery of fast radio burst like emission from the source, and in a period of exceptional high frequency activity. Although not directly taken during a bursting period the counterpart is a factor of ~1.5 to 2.5 brighter than seen at previous epochs with F140W(AB) = $24.65\pm0.02$ mag. We do not detect significant variations of the NIR counterpart within the course of any one orbit (i.e. on minutes-hour timescales), and contemporaneous X-ray observations show SGR 1935+2154 to be at the quiescent level. With a time baseline of 6 years from the first identification of the counter-part we place stringent limits on the proper motion of the source, with a measured proper motion of ${\mu} = 3.1\pm1.5$ mas/yr. The direction of proper motion indicates an origin of SGR 1935+2154 very close to the geometric centre of SNR G57.2+08, further strengthening their association. At an adopted distance of $6.6\pm0.7$ kpc, the corresponding tangential space velocity is ${\nu_T} = 97\pm48$ km/s (corrected for differential Galactic rotation and peculiar Solar motion), although its formal statistical determination may be compromised owing to few epochs of observation. The current velocity estimate places it at the low end of the kick distribution for pulsars, and makes it among the lowest known magnetar kicks. When collating the few-magnetar kick constraints available, we find full consistency between the magnetar kick distribution and the much larger pulsar kick sample, Comment: Accepted for publication in ApJ

The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites., Comment: 19 pages, 16 Figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Short Gamma-Ray Bursts (SGRBs) are produced by the coalescence of compact binary systems which are remnants of massive stars. GRB 160410A is classified as a short-duration GRB with extended emission and is currently the farthest SGRB with a redshift determined from an afterglow spectrum and also one of the brightest SGRBs to date. The fast reaction to the Neil Gehrels Swift Observatory alert allowed us to obtain a spectrum of the afterglow using the X-shooter spectrograph at the Very Large Telescope (VLT). The spectrum shows several absorption features at a redshift of z=1.7177, in addition, we detect two intervening systems at z=1.581 and z=1.444. The spectrum shows ly-alpha in absorption with a column density of log N(HI)=21.2+/-0.2 cm$^{-2}$ which, together with FeII, CII, SiII, AlII and OI, allow us to perform the first study of chemical abundances in a SGRB host galaxy. We determine a metallicity of [X/H]=-2.3+/-0.2 for FeII and -2.5+/-0.2 for SiII and no dust depletion. We also find no evidence for extinction in the afterglow spectral energy distribution (SED) modeling. The environment has a low degree of ionisation and the CIV and SiIV lines are completely absent. We do not detect an underlying host galaxy down to deep limits. Additionally, we compare GRB 160410A to GRB 201221D, another high-z short GRB that shows absorption lines at z=1.045 and an underlying massive host galaxy., Comment: 25 pages, 22 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The typical detection rate of $\sim1$ gamma-ray burst (GRB) per day by the \emph{Fermi} Gamma-ray Burst Monitor (GBM) provides a valuable opportunity to further our understanding of GRB physics. However, the large uncertainty of the \emph{Fermi} localization typically prevents rapid identification of multi-wavelength counterparts. We report the follow-up of 93 \emph{Fermi} GRBs with the Gravitational-wave Optical Transient Observer (GOTO) prototype on La Palma. We selected 53 events (based on favourable observing conditions) for detailed analysis, and to demonstrate our strategy of searching for optical counterparts. We apply a filtering process consisting of both automated and manual steps to 60\,085 candidates initially, rejecting all but 29, arising from 15 events. With $\approx3$ GRB afterglows expected to be detectable with GOTO from our sample, most of the candidates are unlikely to be related to the GRBs. Since we did not have multiple observations for those candidates, we cannot confidently confirm the association between the transients and the GRBs. Our results show that GOTO can effectively search for GRB optical counterparts thanks to its large field of view of $\approx40$ square degrees and its depth of $\approx20$ mag. We also detail several methods to improve our overall performance for future follow-up programs of \emph{Fermi} GRBs.

The advent of wide-field sky surveys has led to the growth of transient and variable source discoveries. The data deluge produced by these surveys has necessitated the use of machine learning (ML) and deep learning (DL) algorithms to sift through the vast incoming data stream. A problem that arises in real-world applications of learning algorithms for classification is imbalanced data, where a class of objects within the data is underrepresented, leading to a bias for over-represented classes in the ML and DL classifiers. We present a recurrent neural network (RNN) classifier that takes in photometric time-series data and additional contextual information (such as distance to nearby galaxies and on-sky position) to produce real-time classification of objects observed by the Gravitational-wave Optical Transient Observer (GOTO), and use an algorithm-level approach for handling imbalance with a focal loss function. The classifier is able to achieve an Area Under the Curve (AUC) score of 0.972 when using all available photometric observations to classify variable stars, supernovae, and active galactic nuclei. The RNN architecture allows us to classify incomplete light curves, and measure how performance improves as more observations are included. We also investigate the role that contextual information plays in producing reliable object classification., Comment: 16 pages, 12 figures, to be published in Monthly Notices of the Royal Astronomical Society

We have adapted the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines to process data from the Gravitational-Wave Optical Transient Observer (GOTO) prototype. In this paper, we describe how we used the Rubin Observatory LSST Science Pipelines to conduct forced photometry measurements on nightly GOTO data. By comparing the photometry measurements of sources taken on multiple nights, we find that the precision of our photometry is typically better than 20~mmag for sources brighter than 16 mag. We also compare our photometry measurements against colour-corrected PanSTARRS photometry, and find that the two agree to within 10~mmag (1$\sigma$) for bright (i.e., $\sim14^{\rm th}$~mag) sources to 200~mmag for faint (i.e., $\sim18^{\rm th}$~mag) sources. Additionally, we compare our results to those obtained by GOTO's own in-house pipeline, {\sc gotophoto}, and obtain similar results. Based on repeatability measurements, we measure a $5\sigma$ L-band survey depth of between 19 and 20 magnitudes, depending on observing conditions. We assess, using repeated observations of non-varying standard SDSS stars, the accuracy of our uncertainties, which we find are typically overestimated by roughly a factor of two for bright sources (i.e., $<15^{\rm th}$~mag), but slightly underestimated (by roughly a factor of 1.25) for fainter sources ($>17^{\rm th}$~mag). Finally, we present lightcurves for a selection of variable sources, and compare them to those obtained with the Zwicky Transient Factory and GAIA. Despite the Rubin Observatory LSST Science Pipelines still undergoing active development, our results show that they are already delivering robust forced photometry measurements from GOTO data., Comment: Accepted for publication in PASA

Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aware classification of transient candidates in difference imaging, and demonstrate its application to the datastream from the GOTO wide-field optical survey. Not only are candidates assigned a well-calibrated probability of being real, but also an associated confidence that can be used to prioritise human vetting efforts and inform future model optimisation via active learning. To fully realise the potential of this architecture, we present a fully-automated training set generation method which requires no human labelling, incorporating a novel data-driven augmentation method to significantly improve the recovery of faint and nuclear transient sources. We achieve competitive classification accuracy (FPR and FNR both below 1%) compared against classifiers trained with fully human-labelled datasets, whilst being significantly quicker and less labour-intensive to build. This data-driven approach is uniquely scalable to the upcoming challenges and data needs of next-generation transient surveys. We make our data generation and model training codes available to the community., Comment: 17 pages, 12 figures, resubmitted to MNRAS following reviewer comments

The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO), whose primary science objective is the optical follow-up of Gravitational Wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near real-time to fully exploit the time-domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this "off-the-shelf" pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to sub-pixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$. These values compare favourably to those obtained using GOTO's primary, in-house pipeline, GOTOPHOTO, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic "obs package" that others can build-upon should they wish to use the LSST Science Pipelines to process data from other facilities., Comment: 17 pages, 13 figures. Accepted for publication in PASA

We report on results of spectropolarimetry of the afterglow of the long gamma-ray burst GRB 191221B, obtained with SALT/RSS and VLT/FORS2, as well as photometry from two telescopes in the MASTER Global Robotic Network, at the MASTER-SAAO (South Africa) and MASTER-OAFA (Argentina) stations. Prompt optical emission was detected by MASTER-SAAO 38 s after the alert, which dimmed from a magnitude (white-light) of ~10 to 16.2 mag over a period of ~10 ks, followed by a plateau phase lasting ~10 ks and then a decline to ~18 mag after 80 ks. The light curve shows complex structure, with four or five distinct breaks in the power-law decline rate. SALT/RSS linear spectropolarimetry of the afterglow began ~2.9 h after the burst, during the early part of the plateau phase of the light curve. Absorption lines seen at ~6010 \r{A} and 5490 \r{A} are identified with the Mg II 2799 \r{A} line from the host galaxy at z=1.15 and an intervening system located at z=0.96. The mean linear polarisation measured over 3400-8000 \r{A} was ~1.5% and the mean equatorial position angle theta ~65 degrees. VLT/FORS2 spectropolarimetry was obtained ~10 h post-burst, during a period of slow decline (alpha = -0.44), and the polarisation was measured to be p = 1.2% and theta = 60 degrees. Two observations with the MeerKAT radio telescope, taken 30 and 444 days after the GRB trigger, detected radio emission from the host galaxy only. We interpret the light curve and polarisation of this long GRB in terms of a slow-cooling forward-shock., Comment: Accepted for publication in MNRAS on 19 June 2021; 14 pages, 6 figures

We predict linear polarization for a radioactively-powered kilonova following the merger of a black hole and a neutron star. Specifically, we perform 3-D Monte Carlo radiative transfer simulations for two different models, both featuring a lanthanide-rich dynamical ejecta component from numerical-relativity simulations while only one including an additional lanthanide-free disk wind component. We calculate polarization spectra for nine different orientations at 1.5, 2.5 and 3.5 d after the merger and in the $0.1-2\,\mu$m wavelength range. We find that both models are polarized at a detectable level 1.5 d after the merger while show negligible levels thereafter. The polarization spectra of the two models are significantly different. The model lacking a disk wind shows no polarization in the optical, while a signal increasing at longer wavelengths and reaching $\sim1\%-6\%$ at $2\,\mu$m depending on the orientation. The model with a disk-wind component, instead, features a characteristic "double-peak" polarization spectrum with one peak in the optical and the other in the infrared. Polarimetric observations of future events will shed light on the debated neutron richness of the disk-wind component. The detection of optical polarization would unambiguously reveal the presence of a lanthanide-free disk-wind component, while polarization increasing from zero in the optical to a peak in the infrared would suggest a lanthanide-rich composition for the whole ejecta. Future polarimetric campaigns should prioritize observations in the first $\sim48$ hours and in the $0.5-2\,\mu$m range, where polarization is strongest, but also explore shorter wavelengths/later times where no signal is expected from the kilonova and the interstellar polarization can be safely estimated., Comment: 8 pages, 4 figures; accepted for publication in MNRAS

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged], Comment: 12 pages, 9 figures. Abridged abstract. Final version published in A&A

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event., Comment: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

Context. EE Cep is one of few eclipsing binary systems with a dark, dusty disk around an invisible object similar to {\epsilon} Aur. The system is characterized by grey and asymmetric eclipses every 5.6 yr, with a significant variation in their photometric depth, ranging from ~ 0 m .5 to ~ 2 m .0. Aims. The main aim of the observational campaign of the EE Cep eclipse in 2014 was to test the model of disk precession (Galan et al. 2012). We expected that this eclipse would be one of the deepest with a depth of ~ 2 m .0. Methods. We collected multicolor observations from almost 30 instruments located in Europe and North America. This photometric data covers 243 nights during and around the eclipse. We also analyse the low- and high-resolution spectra from several instruments. Results. The eclipse was shallow with a depth of 0 m .71 in V-band. The multicolor photometry illustrates small color changes during the eclipse with a total amplitude of order ~ +0 m . 15 in B-I color index. The linear ephemeris for this system is updated by including new times of minima, measured from the three most recent eclipses at epochs E = 9, 10 and 11. New spectroscopic observations were acquired, covering orbital phases around the eclipse, which were not observed in the past and increased the data sample, filling some gaps and giving a better insight into the evolution of the H {\alpha} and NaI spectral line profiles during the primary eclipse. Conclusions. The eclipse of EE Cep in 2014 was shallower than expected 0 m .71 instead of ~ 2 m . 0. This means that our model of disk precession needs revision., Comment: 10 pages, 13 figures, 29 tables in appendix, submited to Astronomy & Astrophysics

Context: Polarimetry is a very powerful tool to uncover various properties of astronomical objects that remain otherwise hidden in standard imaging or spectroscopic observations. However, the reliable measurement of the low polarization signal from astronomical sources requires a good control of spurious instrumental polarization induced by the various components of the optical system and the detector. Aims: We perform a detailed multi-wavelength calibration study of the FORS2 instrument at the VLT operating in imaging polarimetric mode (IPOL) to characterize the spatial instrumental polarization that may affect the study of extended sources. Methods: We use imaging polarimetry of a) high signal-to-noise blank fields BVRI observations during full-moon, when the polarization is expected to be constant across the field-of-view and deviations originate from the instrument and b) a crowded star cluster in broad-band RI and narrow-band H{\alpha} filters, where individual polarization values of each star across the field can be measured. Results: We find an instrumental polarization pattern that increases radially outwards from the optical axis of the instrument reaching up to 1.4% at the edges, depending on the filter. Our results are well approximated by an elliptical paraboloid down to less than {\sim0.05\%} accuracy,and {\sim0.02\%} when using non-analytic fits. We present 2D maps to correct for this spurious instrumental polarization. We also give several tips and tricks to analyze polarimetric measurements of extended sources. Conclusions: FORS2 is a powerful instrument allowing to map the linear polarimetry of extended sources. We present and discuss a methodology to measure the polarization of such sources, and to correct for the spatial polarization induced in the optical system. This methodology could be applied to polarimetric measurements using other dual-beam polarimeters., Comment: 21 pages, 23 figures. Accepted in A&A

A small fraction of candidate tidal disruption events (TDEs) show evidence of powerful relativistic jets, which are particularly pronounced at radio wavelengths, and likely contribute non-thermal emission at a wide range of wavelengths. A non-thermal emission component can be diagnosed using linear polarimetry, even when the total received light is dominated by emission from an accretion disk or disk outflow. In this paper we present Very Large Telescope (VLT) measurements of the linear polarisation of the optical light of jetted TDE Swift J2058+0516. This is the second jetted TDE studied in this manner, after Swift J1644+57. We find evidence of non-zero optical linear polarisation, P_V ~ 8%, a level very similar to the near-infrared polarimetry of Swift J1644+57. These detections provide an independent test of the emission mechanisms of the multiwavelength emission of jetted tidal disruption events., Comment: 7 pages, 4 figures, MNRAS

We report the detections of molecular hydrogen (H$_2$), vibrationally-excited H$_2$ (H$^*_2$), and neutral atomic carbon (CI), in two new afterglow spectra of GRBs\,181020A ($z=2.938$) and 190114A ($z=3.376$), observed with X-shooter at the Very Large Telescope (VLT). Both host-galaxy absorption systems are characterized by strong damped Lyman-$\alpha$ absorbers (DLAs) and substantial amounts of molecular hydrogen with $\log N$(HI, H$_2$) = $22.20\pm 0.05,~20.40\pm 0.04$ (GRB\,181020A) and $\log N$(HI, H$_2$) = $22.15\pm 0.05,~19.44\pm 0.04$ (GRB\,190114A). The DLA metallicites, depletion levels and dust extinctions are [Zn/H] = $-1.57\pm 0.06$, [Zn/Fe] = $0.67\pm 0.03$, and $A_V = 0.27\pm 0.02$\,mag (GRB\,181020A) and [Zn/H] = $-1.23\pm 0.07$, [Zn/Fe] = $1.06\pm 0.08$, and $A_V = 0.36\pm 0.02$\,mag (GRB\,190114A). We then examine the molecular gas content of all known H$_2$-bearing GRB-DLAs and explore the physical conditions and characteristics of these systems. We confirm that H$_2$ is detected in all CI- and H$^*_2$-bearing GRB absorption systems, but that these rarer features are not necessarily detected in all GRB H$_2$ absorbers. We find that a large molecular fraction of $f_{\rm H_2} \gtrsim 10^{-3}$ is required for CI to be detected. The defining characteristic for H$^*_2$ to be present is less clear, though a large H$_2$ column density is an essential factor. We then derive the H$_2$ excitation temperatures of the molecular gas and find that they are relatively low with $T_{\rm ex} \approx 100 - 300$\,K, however, there could be evidence of warmer components populating the high-$J$ H$_2$ levels in GRBs\,181020A and 190114A. Finally, we demonstrate that the otherwise successful X-shooter GRB afterglow campaign is hampered by a significant dust bias excluding the most dust-obscured H$_2$ absorbers from identification [Abridged]., Comment: 15 pages, 15 figures + Appendix. Accepted for publication in A&A

We present X-ray and optical observations of the short duration gamma-ray burst GRB 071227 and its host at $z=0.381$, obtained using \textit{Swift}, Gemini South and the Very Large Telescope. We identify a short-lived and moderately bright optical transient, with flux significantly in excess of that expected from a simple extrapolation of the X-ray spectrum at 0.2-0.3 days after burst. We fit the SED with afterglow models allowing for high extinction and thermal emission models that approximate a kilonova to assess the excess' origins. While some kilonova contribution is plausible, it is not favoured due to the low temperature and high luminosity required, implying superluminal expansion and a large ejecta mass of $\sim 0.1$ M$_{\odot}$. We find, instead, that the transient is broadly consistent with power-law spectra with additional dust extinction of $E(B-V)\sim0.4$ mag, although a possibly thermal excess remains in the \textit{z}-band. We investigate the host, a spiral galaxy with an edge-on orientation, resolving its spectrum along its major axis to construct the galaxy rotation curve and analyse the star formation and chemical properties. The integrated host emission shows evidence for high extinction, consistent with the afterglow findings. The metallicity and extinction are consistent with previous studies of this host and indicate the galaxy is a typical, but dusty, late-type SGRB host., Comment: 16 pages, 13 figures, accepted to MNRAS

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known., Comment: Accepted for publication in MNRAS

We report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date., Comment: 17 pages, 6 figures, Version accepted by ApJ

We present a study of 21 dark gamma-ray burst (GRB) host galaxies, predominantly using X-ray afterglows obtained with the Chandra X-Ray Observatory (CXO) to precisely locate the burst in deep Hubble Space Telescope (HST) imaging of the burst region. The host galaxies are well-detected in F160W in all but one case and in F606W imaging in approx 60 per cent of cases. We measure magnitudes and perform a morphological analysis of each galaxy. The asymmetry, concentration and ellipticity of the dark burst hosts are compared against the host galaxies of optically bright GRBs. In agreement with other studies, we find that dark GRB hosts are redder and more luminous than the bulk of the GRB host population. The distribution of projected spatial offsets for dark GRBs from their host galaxy centroids is comparable to that of optically-bright bursts. The dark GRB hosts are physically larger, more massive and redder, but are morphologically similar to the hosts of bright GRBs in terms of concentration and asymmetry. Our analysis constrains the fraction of high redshift (z greater than 5) GRBs in the sample to approx 14 per cent, implying an upper limit for the whole long-GRB population of less than 4.4 per cent. If dust is the primary cause of afterglow darkening amongst dark GRBs, the measured extinction may require a clumpy dust component in order to explain the observed offset and ellipticity distributions., Comment: Accepted for publication in MNRAS

We investigate the afterglow of GRB 140713A, a gamma-ray burst (GRB) that was detected and relatively well-sampled at X-ray and radio wavelengths, but was not present at optical and near-infrared wavelengths, despite searches to deep limits. We present the emission spectrum of the likely host galaxy at $z = 0.935$ ruling out a high-redshift explanation for the absence of the optical flux detection. Modelling the GRB multi-wavelength afterglow using the radiative transfer hydrodynamics code BOXFIT provides constraints on physical parameters of the GRB jet and its environment, for instance a relatively wide jet opening angle and an electron energy distribution slope $p$ below 2. Most importantly, the model predicts an optical flux about two orders of magnitude above the observed limits. We calculated that the required host extinction to explain the observed limits in the $r$, $i$ and $z$ bands was $A^{\rm host}_{V} > 3.2$ mag, equivalent to $E(B-V)^{\rm host} > 1.0$ mag. From the X-ray absorption we derive that the GRB host extinction is $A^{\rm host}_{\rm V} = 11.6^{+7.5}_{-5.3}$ mag, equivalent to $E(B-V)^{\rm host} = 3.7^{+2.4}_{-1.7}$ mag, which is consistent with the extinction required from our BOXFIT derived fluxes. We conclude that the origin of the optical darkness is a high level of extinction in the line of sight to the GRB, most likely within the GRB host galaxy., Comment: Accepted in MNRAS. 13 pages including 8 figures and a small appendix

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes., Comment: accepted for publication in A&A, 24 pages, 10 figures, tables with the data will be available electronically

Long gamma-ray bursts mark the death of massive stars, as revealed by their association with energetic broad-lined stripped-envelope supernovae. The scarcity of nearby events and the brightness of the GRB afterglow, dominating the first days of emission, have so far prevented the study of the very early stages of the GRB-SN evolution. Here we present detailed, multi-epoch spectroscopic observations of SN 2017iuk, associated with GRB 171205A which display features at extremely high expansion velocities of $\sim$ 100,000 km s$^{-1}$ within the first day after the burst. These high-velocity components are characterized by chemical abundances different from those observed in the ejecta of SN 2017iuk at later times. Using spectral synthesis models developed for the SN 2017iuk, we explain these early features as originating not from the supernova ejecta, but from a hot cocoon generated by the energy injection of a mildly-relativistic GRB jet expanding into the medium surrounding the progenitor star. This cocoon becomes rapidly transparent and is outshone by the supernova emission which starts dominating three days after the burst. These results proves that the jet plays an important role not only in powering the GRB event but also its associated supernova., Comment: 30 pages, 11 figures, 4 tables. Original author manuscript version of a Letter published in Nature journal. Full article available at https://goo.gl/7y9ZeM

We present observations of the optical afterglow of GRB\,170817A, made by the {\it Hubble Space Telescope}, between February and August 2018, up to one year after the neutron star merger, GW170817. The afterglow shows a rapid decline beyond $170$~days, and confirms the jet origin for the observed outflow, in contrast to more slowly declining expectations for `failed-jet' scenarios. We show here that the broadband (radio, optical, X-ray) afterglow is consistent with a structured outflow where an ultra-relativistic jet, with Lorentz factor $\Gamma\gtrsim100$, forms a narrow core ($\sim5^\circ$) and is surrounded by a wider angular component that extends to $\sim15^\circ$, which is itself relativistic ($\Gamma\gtrsim5$). For a two-component model of this structure, the late-time optical decline, where $F \propto t^{-\alpha}$, is $\alpha=2.20\pm0.18$, and for a Gaussian structure the decline is $\alpha=2.45\pm0.23$. We find the Gaussian model to be consistent with both the early $\sim10$ days and late $\gtrsim290$ days data. The agreement of the optical light curve with the evolution of the broadband spectral energy distribution and its continued decline indicates that the optical flux is arising primarily from the afterglow and not any underlying host system. This provides the deepest limits on any host stellar cluster, with a luminosity $\lesssim 4000 L_\odot~(M_{\rm F606W}\gtrsim-4.3)$., Comment: 8 pages, 3 figures, accepted for publication in ApJL

The origin of the luminosity of superluminous supernovae (SLSNe) is an unresolved mystery, and a number of very different physical scenarios (including energy injection from magnetars, collision with a dense circumstellar medium and pair instability-induced explosions) have been invoked. The application of polarimetry to normal SNe has been shown to probe the three-dimensional structure of exploding stars, providing clues to the nature of the explosion mechanism. We report imaging linear polarimetry observations of the Type I SLSN 2017egm, in the galaxy NGC 3191, conducted with the Liverpool Telescope and the RINGO3 instrument. Observations were acquired at four epochs, spanning 4 - 19 days after light-curve maximum, however, polarization was not detected at a level of $>3\sigma$. At +7 and +15 days, and in the average over all epochs, we find a possible polarization signal, detected at a significance of $\approx 2\sigma$ in the "blue" channel. This signal is seen, primarily, in the Stokes $q$ parameter, with a corresponding polarization angle consistent with the orientation of the spiral arm in proximity to the position of SN 2017egm. We interpret this as indicating that any polarization, if present, originates from dust in the host galaxy rather than being intrinsic to the SN itself. Despite its apparent peculiarities, compared to other Type I SLSNe, the polarization characteristics of SN 2017egm are consistent with the previously reported low polarization of other SLSNe of this variety., Comment: 6 pages, 4 figures, MNRAS accepted

We present SPLOT, a small scale pilot survey to test the potential of snapshot (single epoch) linear imaging polarimetry as a supplementary tool to traditional transient follow-up. Transients exist in a vast volume of observational parameter space and polarimetry has the potential to highlight sources of scientific interest and add value to near real-time transient survey streams. We observed a sample of $\sim 50$ randomly selected optical transients with the EFOSC2 and SofI instruments, on the 3.6m New Technology Telescope (NTT) to test the feasibility of the survey. Our sample contained a number of interesting individual sources: a variety of supernovae, X-ray binaries, a tidal disruption event, blazar outbursts, and, by design, numerous transients of unknown nature. We discuss the results, both for the individual sources and the survey in detail. We provide an overview on the success and limitations of SPLOT and also describe a novel calibration method for removing instrumental polarisation effects from Nasymth-mounted telescopes. We find that a SPLOT-like survey would be a benefit to the large scale future transient survey streams such as LSST. The polarimetric measurements have added scientific value to a significant number of the sources and, most importantly, has shown the potential to highlight unclassified transient sources of scientific interest for further study., Comment: Accepted in MNRAS. 20 pages, 10 figures, 3 tables. Also contains a 10 page, two part appendix with unpolarised standard star measurements used for calibration (A) and individual source information and light curves (B). Corrected more typos