Search

Long GRB hosts at z < 1 are usually low-mass, low-metallicity star-forming galaxies. Here we present the most detailed, spatially resolved study of the host of GRB 171205A so far, a grand-design barred spiral galaxy at z = 0.036. Our analysis includes MUSE integral field spectroscopy complemented with high-spatial-resolution UV/VIS HST imaging and CO(1−0) and H I 21 cm 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 H II regions between 5 and 7 kpc from the center and at the end of the bar, inwards from the GRB region; however the H II region hosting the GRB is in the top 10% of the regions with the highest specific star-formation rate. The stellar population at the GRB site has a very young component (< 5 Myr) that contributes 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 H I 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 have consistently low metallicity, high star formation rates, 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 H I. [ABSTRACT FROM AUTHOR]

Context. Gamma-ray bursts are the most energetic known explosions. Although they fade rapidly, they give us the opportunity to measure redshift and important properties of their host galaxies. We report the photometric and spectroscopic study of the Swift GRB 160203A at z = 3.518, 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. Aims. We obtained two optical to near-infrared spectra of the GRB afterglow with X-shooter on ESO/VLT, 18 minutes and 5.7 hours after the burst, allowing us to investigate temporal changes of fine-structure absorption lines. Methods. We measured H I column density log N(HI/cm–2) = 21.75 ± 0.10, and several heavy-element ions along the GRB sightline in the host galaxy, among which Si II, Al II, Al III, C II, Ni II, Si IV, C IV, Zn II and Fe II, and Fe II∗ and Si II∗ fine-structure transitions from energetic levels excited by the afterglow, at the common redshift z = 3.518. We measured [M/H]TOT = –0.78 ± 0.13 and a [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 ɀ = 1.03, ɀ = 1.26, ɀ = 1.98, ɀ = 1.99, ɀ = 2.20, and ɀ = 2.83. We could not measure significant variability in the strength of 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. Results. GRB 160203A belongs to a growing sample of GRBs with medium resolution spectroscopy, provided by the Swift/X-shooter legacy programme, which enables a 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. [ABSTRACT FROM AUTHOR]

Future detection of high-redshift gamma-ray bursts (GRBs) will be an important tool for studying the early Universe. Fast and accurate redshift estimation for detected GRBs is key for encouraging rapid follow-up observations by ground- and space-based telescopes. Low-redshift dusty interlopers pose the biggest challenge for GRB redshift estimation using broad photometric bands, as their high extinction can mimic a high-redshift GRB. To assess false alarms of high-redshift GRB photometric measurements, we simulate and fit a variety of GRBs using phozzy , a simulation code developed to estimate GRB photometric redshifts, and test the ability to distinguish between high- and low-redshift GRBs when using simultaneously observed photometric bands. We run the code with the wavelength bands and instrument parameters for the Photo-z Infrared Telescope (PIRT), an instrument designed for the Gamow mission concept. We explore various distributions of host galaxy extinction as a function of redshift, and their effect on the completeness and purity of a high-redshift GRB search with the PIRT. We find that for assumptions based on current observations, the completeness and purity range from ∼82 to 88 per cent and from ∼84 to |$\gt 99~{{\ \rm per\ cent}}$|⁠ , respectively. For the priors optimized to reduce false positives, only |$\sim 0.6~{{\ \rm per\ cent}}$| of low-redshift GRBs will be mistaken as a high-redshift one, corresponding to ∼1 false alarm per 500 detected GRBs. [ABSTRACT FROM AUTHOR]

Context. GRB 220627A is a rare burst with two distinct γ-ray emission episodes separated by almost 1000 s that triggered the Fermi Gamma-ray Burst Monitor twice. High-energy GeV emission was detected by the Fermi Large Area Telescope coincident with the first emission episode but not the second. The discovery of the optical afterglow with MeerLICHT led to MUSE observations which secured the burst redshift to z = 3.08, making this the most distant ultra-long gamma-ray burst (GRB) detected to date. Aims. The progenitors of some ultra-long GRBs have been suggested in the literature to be different to those of normal long GRBs. Our aim is to determine whether the afterglow and host properties of GRB 220627A agree with this interpretation. Methods. We performed empirical and theoretical modelling of the afterglow data within the external forward shock framework, and determined the metallicity of the GRB environment through modelling the absorption lines in the MUSE spectrum. Results. Our optical data show evidence for a jet break in the light curve at ∼1.2 days, while our theoretical modelling shows a preference for a homogeneous circumburst medium. Our forward shock parameters are typical for the wider GRB population, and we find that the environment of the burst is characterised by a sub-solar metallicity. Conclusions. Our observations and modelling of GRB 220627A do not suggest that a different progenitor compared to the progenitor of normal long GRBs is required. We find that more observations of ultra-long GRBs are needed to determine if they form a separate population with distinct prompt and afterglow features, and possibly distinct progenitors. [ABSTRACT FROM AUTHOR]

BL Lacertae objects represent a large fraction (22%) of γ-ray sources in the Third Fermi Large Area Telescope catalog (3FGL). Nearly half of the BL Lac object population remains without a redshift because of their featureless optical spectra. We aim to increase the number of BL Lac objects with a redshift measurement by using the photometric technique. For this work, we have used six Swift-UVOT filters and SDSS g′, r′, i′, and z′ optical filters mounted on the 0.65 m SARA–CTIO located in Chile and the 1.0 m SARA–ORM in the Canary Islands. A sample of 45 sources was selected from the 3FGL catalog for which photometry was performed in 10 optical and UV filters to obtain redshift measurements. We found three sources with z > 1.3, while reliable upper limits have been provided for 17 sources. The results presented here bring the total number of high-z (z > 1.3) BL Lac objects to 29. [ABSTRACT FROM AUTHOR]

We report the detections of molecular hydrogen (H2), vibrationally-excited H2 (H 2∗ 2 ∗ $ _{2}^{*} $), and neutral atomic carbon (C I), an efficient tracer of molecular gas, 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-α absorbers (DLAs) and substantial amounts of molecular hydrogen with logN(H I, H2) = 22.20 ± 0.05,  20.40 ± 0.04 (GRB 181020A) and logN(H I, H2) = 22.15 ± 0.05,  19.44 ± 0.04 (GRB 190114A). The DLA metallicites, depletion levels, and dust extinctions are within the typical regimes probed by GRBs with [Zn/H] = −1.57 ± 0.06, [Zn/Fe] = 0.67 ± 0.03, and AV = 0.27 ± 0.02 mag (GRB 181020A) and [Zn/H] = −1.23 ± 0.07, [Zn/Fe] = 1.06 ± 0.08, and AV = 0.36 ± 0.02 mag (GRB 190114A). In addition, we examine the molecular gas content of all known H2-bearing GRB-DLAs and explore the physical conditions and characteristics required to simultaneously probe C I and H 2∗ 2 ∗ $ _{2}^{*} $. We confirm that H2 is detected in all C I- and H 2∗ 2 ∗ $ _{2}^{*} $ -bearing GRB absorption systems, but that these rarer features are not necessarily detected in all GRB H2 absorbers. We find that a large molecular fraction of fH2 ≳ 10−3 is required for C I to be detected. The defining characteristic for H 2∗ 2 ∗ $ _{2}^{*} $ to be present is less clear, though a large H2 column density is an essential factor. We also find that the observed line profiles of the molecular-gas tracers are kinematically "cold", with small velocity offsets of δv <  20 km s−1 from the bulk of the neutral absorbing gas. We then derive the H2 excitation temperatures of the molecular gas and find that they are relatively low with Tex ≈ 100−300 K, however, there could be evidence of warmer components populating the high-J H2 levels in GRBs 181020A and 190114A. Finally, we demonstrate that even though the X-shooter GRB afterglow campaign has been successful in recovering several H2-bearing GRB-host absorbers, this sample is still hampered by a significant dust bias excluding the most dust-obscured H2 absorbers from identification. C I and H 2∗ 2 ∗ $ _{2}^{*} $ could open a potential route to identify molecular gas even in low-metallicity or highly dust-obscured bursts, though they are only efficient tracers for the most H2-rich GRB-host absorption systems. [ABSTRACT FROM AUTHOR]

Context. Damped Lyman-α (DLA) absorption-line systems at the redshifts of gamma-ray burst (GRB) afterglows offer a unique way to probe the physical conditions within star-forming galaxies in the early Universe. Aims. Here we built up a large sample of 22 GRBs at redshifts z > 2 observed with VLT/X-shooter in order to determine the abundances of hydrogen, metals, dust, and molecular species. This allows us to study the metallicity and dust depletion effects in the neutral interstellar medium at high redshift and to answer the question of whether (and why) there might be a lack of H2 in GRB-DLAs. Methods. We developed new methods based on the Bayesian inference package, PyMC, to FIT absorption lines and measure the column densities of different metal species as well as atomic and molecular hydrogen. The derived relative abundances are used to FIT dust depletion sequences and determine the dust-to-metals ratio and the host-galaxy intrinsic visual extinction. Additionally, we searched for the absorption signatures of vibrationally-excited H2 and carbon monoxide. Results. We find that there is no lack of H2-bearing GRB-DLAs. We detect absorption lines from molecular hydrogen in 6 out of 22 GRB afterglow spectra, with molecular fractions ranging between f ≃ 5 × 10−5 and f ≃ 0.04, and claim tentative detections in three additional cases. For the remainder of the sample, we measure, depending on S/N, spectral coverage and instrumental resolution, more or less stringent upper limits. The GRB-DLAs in our sample have on average low metallicities, [X/H]¯ ≈ −1.3 $\overline{\textrm{[X/H]}}\approx -1.3$ [ X/H ] ¯ ≈ − 1.3 , comparable to the population of extremely-strong QSO-DLAs (log N(H I) > 21.5). Furthermore, H2-bearing GRB-DLAs are found to be associated with significant dust extinction, AV > 0.1 mag, and dust-to-metals ratios DTM > 0.4 ${\mathcal{DTM}} > 0.4$ D T M > 0.4 , confirming the importance of dust grains for the production of molecules. All these systems exhibit neutral hydrogen column densities log N(H I) > 21.7. The overall fraction of H2 detections in GRB-DLAs is ≥ 27% (41% including tentative detections), which is three to four times larger than in the general QSO-DLA population. For 2 < z < 4, and considering column densities log N(H I) > 21.7, the H2 detection fraction is 60–80% in GRB-DLAs and in extremely strong QSO-DLAs. This is likely due to the fact that both GRB- and QSO-DLAs with high neutral hydrogen column densities are probed by sight-lines with small impact parameters, indicating that the absorbing gas is associated with the inner regions of the absorbing galaxy, where the gas pressure is higher and the conversion of H I to H2 takes place. In the case of GRB hosts, this diffuse molecular gas is located at distances ≳ 500 pc from the GRB and hence is unrelated to the star-forming region where the event occurred. [ABSTRACT FROM AUTHOR]

Long γ-ray bursts are associated with energetic, broad-lined, stripped-envelope supernovae1,2 and as such mark the death of massive stars. The scarcity of such events nearby and the brightness of the γ-ray burst afterglow, which dominates the emission in the first few days after the burst, have so far prevented the study of the very early evolution of supernovae associated with γ-ray bursts3. In hydrogen-stripped supernovae that are not associated with γ-ray bursts, an excess of high-velocity (roughly 30,000 kilometres per second) material has been interpreted as a signature of a choked jet, which did not emerge from the progenitor star and instead deposited all of its energy in a thermal cocoon4. Here we report multi-epoch spectroscopic observations of the supernova SN 2017iuk, which is associated with the γ-ray burst GRB 171205A. Our spectra display features at extremely high expansion velocities (around 115,000 kilometres per second) within the first day after the burst5,6. Using spectral synthesis models developed for SN 2017iuk, we show that these features are characterized by chemical abundances that differ from those observed in the ejecta of SN 2017iuk at later times. We further show that the high-velocity features originate from the mildly relativistic hot cocoon that is generated by an ultra-relativistic jet within the γ-ray burst expanding and decelerating into the medium that surrounds the progenitor star7,8. This cocoon rapidly becomes transparent9 and is outshone by the supernova emission, which starts to dominate the emission three days after the burst. Multi-epoch observations of a supernova associated with a γ-ray burst reveal spectral features at extremely high expansion velocities within the first day after the burst, indicative of a choked jet. [ABSTRACT FROM AUTHOR]

Context. Classical novae (CNe) represent the major class of supersoft X-ray sources (SSSs) in the central region of our neighbouring galaxy M 31. Aims. We performed a dedicated monitoring of the M 31 central region, which aimed to detect SSS counterparts of CNe, with XMM-Newton and Chandra between Nov. and Mar. of the years 2009/10, 2010/11, and 2011/12. Methods. We systematically searched our data for X-ray counterparts of CNe and determined their X-ray light curves and also their spectral properties in the case of XMM-Newton data. Additionally, we determined luminosity upper limits for all previously known X-ray emitting novae, which are not detected anymore, and for all CNe in our field of view with recent optical outbursts. Results. In total, we detected 24 novae in X-rays. Seven of these sources were known from previous observations, including the M 31 nova with the longest SSS phase, M31N 1996-08b, which was found to fade below our X-ray detection limit 13.8 yr after outburst. Of the new discoveries, several novae exhibit significant variability in their short-term X-ray light curves with one object showing a suspected period of about 1.3 h.We studied the SSS state of the most recent outburst of a recurrent nova, which had previously shown the shortest time ever observed between two outbursts (~5 yr). The total number of M 31 novae with X-ray counterpart was increased to 79, and we subjected this extended catalogue to detailed statistical studies. Four previously indicated correlations between optical and X-ray parameters could be confirmed and improved. Furthermore, we found indications that the multi-dimensional parameter space of nova properties might be dominated by a single physical parameter, and we provide interpretations and suggest implications. We studied various outliers from the established correlations and discuss evidence of a different X-ray behaviour of novae in the M 31 bulge and disk. Conclusions. Exploration of the multi-wavelength parameter space of optical and X-ray measurements is shown to be a powerful tool for examining properties of extragalactic nova populations. While there are hints that the different stellar populations of M 31 (bulge vs. disk) produce dissimilar nova outbursts, there is also growing evidence that the overall behaviour of an average nova might be understood in surprisingly simple terms. [ABSTRACT FROM AUTHOR]

Context. The central field of the Andromeda galaxy (M31) was monitored from 2006 to 2012 using the Chandra HRC-I detector (about 0.1-10 keV energy range) with the main aim of detecting X-rays from optical novae. We present a systematic analysis of all X-ray sources found in the 41 nova monitoring observations, along with 23 M31 central field HRC-I observations available from the Chandra data archive starting in December 1999. Aims. Based on these observations, we studied the X-ray long-term variability of the source population and especially of X-ray binaries in M31. Methods. We created a catalogue of sources detected in the 64 available observations that adds up to a total exposure time of about 1Ms. To study the variability, we developed a processing pipeline to derive long-term Chandra HRC-I light curves for each source over the 13 years of observations. We also searched for extended X-ray sources in the merged images. Results. We present a point-source catalogue containing 318 X-ray sources with detailed long-term variability information, 28 of which are published for the first time. The spatial and temporal resolution of the catalogue allows us to classify 115 X-ray binary candidates showing high X-ray variability or even outbursts, as well as 14 globular cluster X-ray binary candidates showing no significant variability. The analysis may suggest that outburst sources are less frequent in globular clusters than in the field of M31. We detected seven supernova remnants, one of which is a new candidate, and also resolved the first X-rays from a known radio supernova remnant. In addition to 33 known optical nova/X-ray source correlations, we discovered one previously unknown super-soft X-ray outburst and several new nova candidates. Conclusions. The catalogue contains a large sample of detailed long-term X-ray light curves in the M31 central field, which helps in understanding the X-ray population of our neighbouring spiral galaxy M31. [ABSTRACT FROM AUTHOR]

Context. Classical novae (CNe) represent the main class of supersoft X-ray sources (SSSs) in the central region of our neighbouring galaxy M 31. Only three confirmed novae and three SSSs have been discovered in globular clusters (GCs) of any galaxy so far, of which one nova and two SSSs (including the nova) were found in M 31 GCs. Aims. To study the SSS state of CNe we carried out a high-cadence X-ray monitoring of the M 31 central area with XMM-Newton and Chandra. This project is supplemented by regular optical monitoring programmes at various observatories. Methods. We analysed X-ray and optical monitoring data of a new transient X-ray source in the M 31 GC Bol 126, discovered serendipitously in Swift observations. Our optical data set was based on regular M 31 monitoring programmes from five different small telescopes and was reduced using a homogeneous method. Additionally, we made use of Pan-STARRS 1 data obtained during the PAndromeda survey. We extracted light curves of the source in the optical and X-rays, as well as X-ray spectra. Results. Our observations reveal that the X-ray source in Bol 126 is the third SSS in an M 31 GC and can be confirmed as the second CN in the M 31 GC system. This nova is named M31N 2010-10f. Its properties in the X-ray (high black-body temperature, short SSS phase) and optical (relatively high maximum magnitude, fast decline) regimes agree with a massive white dwarf (MWD ≳ 1.3 M⊕) in the binary system. Incorporating the data on previously found (suspected) novae in M 31 GCs we used our high-cadence X-ray monitoring observations to estimate a tentative nova rate in the M 31 GC system of 0.05 yr-1 GC-1. An optical estimate, based on the recent 10.5-year We CAPP survey, gives a lower nova rate, which is compatible with the X-ray rate on the 95% confidence level. Conclusions. Although still based on small-number statistics, there is growing evidence that the nova rate in GCs is higher than expected from primordial binary formation and under conditions as in the field. Dynamical binary formation and/or additional accretion from the intracluster medium are possible scenarios for an increased nova rate, but observational confirmation for this enhancement has been absent, so far. Regular X-ray monitoring observations of M 31 provide a promising strategy to find these novae. [ABSTRACT FROM AUTHOR]

We report on follow-up observations of 20 short-duration gamma-ray bursts (GRBs; T90 < 2 s) performed in g'r'i'z'JHKs with the Gamma-Ray Burst Optical Near-Infrared Detector (GROND) between mid-2007 and the end of 2010. This is one of the most comprehensive data sets on GRB afterglow observations of short bursts published so far. In three cases, GROND was on target within less than 10 min after the trigger, leading to the discovery of the afterglow of GRB 081226A and its faint underlying host galaxy. In addition, GROND was able to image the optical afterglow and follow the light curve evolution in five further cases: GRBs 090305, 090426, 090510, 090927, and 100117A. In all other cases, optical/near-infrared upper limits can be provided on the afterglow magnitudes. After shifting all light curves to a common redshift, we find that the optical luminosities of the six events with light curves group into two subsamples. GRBs 090426 and 090927 are situated in the regime occupied by long-duration events (collapsars), while the other four bursts occupy the parameter space typical for merger events, confirming that the short-burst population is contaminated by collapsar events. Three of the aforementioned six bursts with optical light curves show a break: GRBs 090426 and 090510 (Papers I and II) as well as GRB 090305. For GRB 090927, no break is seen in the optical/X-ray light curve until about 150 ks/600 ks after the burst. The GROND multi-color data support the view that this burst is related to a collapsar event. A decay slope of the optical afterglow of GRB 100117A could be measured. For all six GRBs a lower limit on the corresponding jet opening angle can be set. Using these data supplemented by about ten events taken from the literature, we compare the jet half-opening angles of long and short bursts. We find tentative evidence that short bursts have wider opening angles than long bursts. However, the statistics are still very poor. [ABSTRACT FROM AUTHOR]

Aims. We study the populations of massive stars in the Carina region and their energetic feedback and ejection of 26Al. Methods. We did a census of the stellar populations in young stellar clusters within a few degrees of the Carina nebula. For each star we estimated the mass, based on the spectral type and the host cluster age. We used population synthesis to calculate the energetic feedback and ejection of 26Al from the winds of the massive stars and their supernova explosions. We used 7 years of INTEGRAL observations to measure the 26Al signal from the region. Results. The INTEGRAL 26Al signal is not significant with a best-fit value of ∼1.5 ± 1.0 × 10-5 ph cm-2 s-1, approximately half of the published Compton Gamma Ray Observatory (CGRO) result, but in agreement with the latest CGRO estimates. Our analysis of the stellar populations in the young clusters leads to an expected signal of ∼half the observed value, but the results are consistent within 2σ. We find that the fraction of 26Al ejected in Wolf-Rayet winds is high, and the observed signal is unlikely to be caused by 26Al ejected in supernovae alone, indicating a strong wind ejection of 26Al. Due to the lack of prominent O stars, regions with ages ≳10 Myr are often neglected in studies of OB associations. We find that in the Carina region such clusters contribute significantly to the stellar mass and the energetics of the region. [ABSTRACT FROM AUTHOR]

ABSTRACT We report on the surprisingly high metallicity measured in two absorption systems at high redshift, detected in the Very Large Telescope spectrum of the afterglow of the gamma-ray burst (GRB) GRB 090323. The two systems, at redshift z= 3.5673 and 3.5774 (separation Δ v≈ 660 km s−1), are dominated by the neutral gas in the interstellar medium of the parent galaxies. From the singly ionized zinc and sulphur, we estimate oversolar metallicities of [Zn/H] =+0.29 ± 0.10 and [S/H] = +0.67 ± 0.34, in the blue and red absorber, respectively. These are the highest metallicities ever measured in galaxies at z > 3. We propose that the two systems trace two galaxies in the process of merging, whose star formation and metallicity are heightened by the interaction. This enhanced star formation might also have triggered the birth of the GRB progenitor. As typically seen in star-forming galaxies, the fine-structure absorption Si ii* is detected, both in z= 3.5774 ± 0.0005 and 3.5673 ± 0.0003. From the rest-frame ultraviolet emission in the GRB location, we derive a relatively high, not corrected for dust extinction, star formation rate ≈6 M⊙ yr−1. These properties suggest a possible connection between some high-redshift GRB host galaxies and high- z massive submillimetre galaxies, which are characterized by disturbed morphologies and high metallicities. Our result provides additional evidence that the dispersion in the chemical enrichment of the Universe at high redshift is substantial, with the existence of very metal-rich galaxies less than two billion years after the big bang. [ABSTRACT FROM AUTHOR]

Context. Observations of the γ-ray sky with Fermi led to significant advances towards understanding blazars, the most extreme class of active galactic nuclei. A large fraction of the population detected by Fermi is formed by BL Lacertae (BL Lac) objects, whose sample has always suffered from a severe redshift incompleteness due to the quasi-featureless optical spectra. Aims. Our goal is to provide a significant increase of the number of confirmed high-redshift BL Lac objects contained in the 2 LAC Fermi/LAT cataloge. Methods. For 103 Fermi/LAT blazars, photometric redshifts using spectral energy distribution fitting have been obtained. The photometry includes 13 broad-band filters from the far ultraviolet to the near-IR observed with Swift/UVOT and the multi-channel imager GROND at the MPG/ESO 2.2 m telescope. Data have been taken quasi-simultaneously and the remaining source-intrinsic variability has been corrected for. Results. We release the UV-to-near-IR 13-band photometry for all 103 sources and provide redshift constraints for 75 sources without previously known redshift. Out of those, eight have reliable photometric redshifts at z ≳ 1.3, while for the other 67 sources we provide upper limits. Six of the former eight are BL Lac objects, which quadruples the sample of confirmed high-redshift BL Lac. This includes three sources with redshifts higher than the previous record for BL Lac, including CRATES J0402-2615, with the best-fit solution at z ≈ 1.9. [ABSTRACT FROM AUTHOR]