Search

A dedicated gamma-ray burst (GRB) afterglow observing program was performed between 2007 and 2016 with GROND, a seven-channel optical and near-infrared imager at the 2.2m telescope of the Max-Planck Society at ESO/La Silla. In this first of a series of papers, we describe the GRB observing plan, providing first readings of all so far unpublished GRB afterglow measurements and some observing statistics. In total, we observed 514 GRBs with GROND, including 434 Swift-detected GRBs, representing 81\% of the observable Swift sample. For GROND-observations within 30 min of the GRB trigger, the optical/NIR afterglow detection rate is 81\% for long- and 57\% for short-duration GRBs. We report the discovery of ten new GRB afterglows plus one candidate, along with redshift estimates (partly improved) for four GRBs and new host detections for seven GRBs. We identify the (already known) afterglow of GRB 140209A as the sixth GRB exhibiting a 2175 Angstroem dust feature. As a side result, we identified two blazars, with one at a redshift of z=3.8 (in the GRB 131209A field)., Comment: 20 pages, 31 figures, accepted for publication in Astronomy & Astrophysics

PICO bubble chambers have exceptional sensitivity to inelastic dark matter-nucleus interactions due to a combination of their extended nuclear recoil energy detection window from a few keV to $O$(100 keV) or more and the use of iodine as a heavy target. Inelastic dark matter-nucleus scattering is interesting for studying the properties of dark matter, where many theoretical scenarios have been developed. This study reports the results of a search for dark matter inelastic scattering with the PICO-60 bubble chambers. The analysis reported here comprises physics runs from PICO-60 bubble chambers using CF$_{3}$I and C$_{3}$F$_{8}$. The CF$_{3}$I run consisted of 36.8 kg of CF$_{3}$I reaching an exposure of 3415 kg-day operating at thermodynamic thresholds between 7 and 20 keV. The C$_{3}$F$_{8}$ runs consisted of 52 kg of C$_{3}$F$_{8}$ reaching exposures of 1404 kg-day and 1167 kg-day running at thermodynamic thresholds of 2.45 keV and 3.29 keV, respectively. The analysis disfavors various scenarios, in a wide region of parameter space, that provide a feasible explanation of the signal observed by DAMA, assuming an inelastic interaction, considering that the PICO CF$_{3}$I bubble chamber used iodine as the target material., Comment: 7 pages, 3 figures

The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. The PICO Collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C$_3$F$_8$ from calibration data taken with 5 distinct neutron spectra at various thermodynamic thresholds ranging from 2.1 keV to 3.9 keV. Instead of assuming any particular functional forms for the nuclear recoil efficiency, a generalized piecewise linear model is proposed with systematic errors included as nuisance parameters to minimize model-introduced uncertainties. A Markov-Chain Monte-Carlo (MCMC) routine is applied to sample the nuclear recoil efficiency for fluorine and carbon at 2.45 keV and 3.29 keV thermodynamic thresholds simultaneously. The nucleation efficiency for fluorine was found to be $\geq 50\, \%$ for nuclear recoils of 3.3 keV (3.7 keV) at a thermodynamic Seitz threshold of 2.45 keV (3.29 keV), and for carbon the efficiency was found to be $\geq 50\, \%$ for recoils of 10.6 keV (11.1 keV) at a threshold of 2.45 keV (3.29 keV). Simulated data sets are used to calculate a p-value for the fit, confirming that the model used is compatible with the data. The fit paradigm is also assessed for potential systematic biases, which although small, are corrected for. Additional steps are performed to calculate the expected interaction rates of WIMPs in the PICO-60 detector, a requirement for calculating WIMP exclusion limits., Comment: 17 pages, 22 figures, 5 tables

Many compelling models predict dark matter coupling to the electromagnetic current through higher multipole interactions, while remaining electrically neutral. Different multipole couplings have been studied, among them anapole moment, electric and magnetic dipole moments, and millicharge. This study sets limits on the couplings for these photon-mediated interactions using non-relativistic contact operators in an effective field theory framework. Using data from the PICO-60 bubble chamber leading limits for dark matter masses between 2.7 GeV/c$^2$ and 24 GeV/c$^2$ are reported for the coupling of these photon-mediated dark matter-nucleus interactions. The detector was filled with 52 kg of C$_3$F$_8$ operating at thermodynamic thresholds of 2.45 keV and 3.29 keV, reaching exposures of 1404 kg-day and 1167 kg-day, respectively.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution., Comment: 9 pages, 3 figures, 1 table. To appear in Eur.Phys.J. C

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$., Comment: 18 pages, 12 figures

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude., Comment: 9 pages, 8 figures

The identification of non-signal events is a major hurdle to overcome for bubble chamber dark matter experiments such as PICO-60. The current practice of manually developing a discriminator function to eliminate background events is difficult when available calibration data is frequently impure and present only in small quantities. In this study, several different discriminator input/preprocessing formats and neural network architectures are applied to the task. First, they are optimized in a supervised learning context. Next, two novel semi-supervised learning algorithms are trained, and found to replicate the Acoustic Parameter (AP) discriminator previously used in PICO-60 with a mean of 97% accuracy., Comment: 27 pages, 10 figures

Twenty years ago, GRB 980425/SN 1998bw revealed that long Gamma-Ray Bursts (GRBs) are physically associated with broad-lined type Ic supernovae. Since then more than 1000 long GRBs have been localized to high angular precision, but only in about 50 cases the underlying supernova (SN) component was identified. Using the multi-channel imager GROND (Gamma-Ray Burst Optical Near-Infrared Detector) at ESO/La Silla, during the last ten years we have devoted a substantial amount of observing time to reveal and to study SN components in long-GRB afterglows. Here we report on four more GRB-SNe (associated with GRBs 071112C, 111228A, 120714B, and 130831A) which were discovered and/or followed-up with GROND and whose redshifts lie between z=0.4 and 0.8. We study their afterglow light curves, follow the associated SN bumps over several weeks, and characterize their host galaxies. Using SN 1998bw as a template, the derived SN explosion parameters are fully consistent with the corresponding properties of the so-far known GRB-SN ensemble, with no evidence for an evolution of their properties as a function of redshift. In two cases (GRB 120714B/SN 2012eb at z=0.398 and GRB 130831A/SN 2013fu at z=0.479) additional Very Large Telescope (VLT) spectroscopy of the associated SNe revealed a photospheric expansion velocity at maximum light of about 40 000 and 20 000 km/s, respectively. For GRB 120714B, which was an intermediate-luminosity burst, we find additional evidence for a blackbody component in the light of the optical transient at early times, similar to what has been detected in some GRB-SNe at lower redshifts., Comment: submitted to Astron. Astroph., revised version

Long gamma-ray bursts give us the chance to study both their extreme physics and the star-forming galaxies in which they form. GRB 100418A, at a z = 0.6239, had a bright optical and radio afterglow, and a luminous star-forming host galaxy. This allowed us to study the radiation of the explosion as well as the interstellar medium of the host both in absorption and emission. We collected photometric data from radio to X-ray wavelengths to study the evolution of the afterglow and the contribution of a possible supernova and three X-shooter spectra obtained during the first 60 hr. The light curve shows a very fast optical rebrightening, with an amplitude of 3 magnitudes, starting 2.4 hr after the GRB onset. This cannot be explained by a standard external shock model and requires other contributions, such as late central-engine activity. Two weeks after the burst we detect an excess in the light curve consistent with a SN with peak absolute magnitude M_V = -18.5 mag, among the faintest GRB-SNe detected to date. The host galaxy shows two components in emission, with velocities differing by 130 km s^-1, but otherwise having similar properties. While some absorption and emission components coincide, the absorbing gas spans much higher velocities, indicating the presence of gas beyond the star-forming regions. The host has a star-formation rate of 12.2 M_sol yr^-1, a metallicity of 12 + log(O/H) = 8.55 and a mass of 1.6x10^9 M_sol. GRB 100418A is a member of a class of afterglow light curves which show a steep rebrightening in the optical during the first day, which cannot be explained by traditional models. Its very faint associated SN shows that GRB-SNe can have a larger dispersion in luminosities than previously seen. Furthermore, we have obtained a complete view of the host of GRB 100418A owing to its spectrum, which contains a remarkable number of both emission and absorption lines., Comment: 23 pages, 14 figures, accepted for publication in A&A

Afterglows of Gamma-Ray Bursts (GRBs) are simple in the most basic model, but can show many complex features. The ultra-long duration GRB 111209A, one of the longest GRBs ever detected, also has the best-monitored afterglow in this rare class of GRBs. We want to address the question whether GRB 111209A was a special event beyond its extreme duration alone, and whether it is a classical GRB or another kind of high-energy transient. The afterglow may yield significant clues. We present afterglow photometry obtained in seven bands with the GROND imager as well as in further seven bands with the UVOT telescope on-board the Neil Gehrels Swift Observatory. The light curve is analysed by multi-band modelling and joint fitting with power-laws and broken power-laws, and we use the contemporaneous GROND data to study the evolution of the spectral energy distribution. We compare the optical afterglow to a large ensemble we have analysed in earlier works, and especially to that of another ultra-long event, GRB 130925A. We furthermore undertake a photometric study of the host galaxy. We find a strong, chromatic rebrightening event at approx 0.8 days after the GRB, during which the spectral slope becomes redder. After this, the light curve decays achromatically, with evidence for a break at about 9 days after the trigger. The afterglow luminosity is found to not be exceptional. We find that a double-jet model is able to explain the chromatic rebrightening. The afterglow features have been detected in other events and are not unique. The duration aside, the GRB prompt emission and afterglow parameters of GRB 111209A are in agreement with the known distributions for these parameters. While the central engine of this event may differ from that of classical GRBs, there are multiple lines of evidence pointing to GRB 111209A resulting from the core-collapse of a massive star with a stripped envelope., Comment: A&A, in press, expanded discussion on energetics, and updated otherwise. 32 pages, 17 pages main paper, 2 pages Appendix, 11 pages data tables. arXiv admin note: substantial text overlap with arXiv:1606.06791

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results., Comment: 6 pages, 6 figures

GRB 111209A, one of the longest Gamma-Ray Bursts (GRBs) ever observed, is linked to SN 2011kl, the most luminous GRB-Supernova (SN) detected so far, which shows evidence for being powered by a magnetar central engine. We place SN 2011kl into the context of large samples of SNe, addressing in more detail the question of whether it could be radioactively powered, and whether it represents an extreme version of a GRB-SN or an underluminous Superluminous SN (SLSN). We model SN 2011kl using SN 1998bw as a template and derive a bolometric light curve including near-infrared data. We compare the properties of SN 2011kl to literature results on stripped-envelope and superluminous supernovae. Comparison in the k,s context, i.e., comparing it to SN 1998bw templates in terms of luminosity and light-curve stretch, clearly shows SN 2011kl is the most luminous GRB-SN to date, and it is spectrally very dissimilar to other events, being significantly bluer/hotter. Although SN 2011kl does not reach the classical luminosity threshold of SLSNe and evolves faster than any of them, it resembles SLSNe more than the classical GRB-associated broad-lined Type Ic SNe in several aspects. GRB 111209A was a very energetic event, both at early (prompt emission) and at very late (SN) times. We have shown in a further publication that with the exception of the extreme duration, the GRB and afterglow parameters are in agreement with the known distributions for these parameters. SN 2011kl, on the other hand, is exceptional both in luminosity and spectral characteristics, indicating that GRB 111209A was likely not powered by a standard-model collapsar central engine, further supporting our earlier conclusions. Instead, it reveals the possibility of a direct link between GRBs and SLSNe., Comment: A&A, in press; v4: updated, text compressed, references removed upon request of the referee; 21 pages: 17 pages article, 2 pages tables

New data are reported from a second run of the 2-liter PICO-2L C$_3$F$_8$ bubble chamber with a total exposure of 129$\,$kg-days at a thermodynamic threshold energy of 3.3$\,$keV. These data show that measures taken to control particulate contamination in the superheated fluid resulted in the absence of the anomalous background events observed in the first run of this bubble chamber. One single nuclear-recoil event was observed in the data, consistent both with the predicted background rate from neutrons and with the observed rate of unambiguous multiple-bubble neutron scattering events. The chamber exhibits the same excellent electron-recoil and alpha decay rejection as was previously reported. These data provide the most stringent direct detection constraints on weakly interacting massive particle (WIMP)-proton spin-dependent scattering to date for WIMP masses $<$ 50$\,$GeV/c$^2$.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out., Comment: v3 to reflect published version

[Abridged] We present data and initial results from VLT/X-Shooter emission-line spectroscopy of 96 GRB-selected galaxies at 0.1

Aims. The aim of the study is to constrain the physics of gamma-ray bursts (GRBs) by analysing the multi-wavelength afterglow data set of GRB 121024A that covers the full range from radio to X-rays. Methods. Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measured the three characteristic break frequencies of the synchrotron spectrum. We used six epochs of combined XRT and GROND data to constrain the temporal slopes, the dust extinction, the X-ray absorption, and the spectral slope with high accuracy. Two more epochs of combined data from XRT, GROND, APEX, CARMA, and EVLA were used to set constraints on the break frequencies and therefore on the micro-physical and dynamical parameters. Results. The XRT and GROND light curves show a simultaneous and achromatic break at around 49 ks. As a result, the crossing of the synchrotron cooling break is no suitable explanation for the break in the light curve. The multi wavelength data allow us to test two plausible scenarios explaining the break: a jet break, and the end of energy injection. The jet-break scenario requires a hard electron spectrum, a very low cooling break frequency, and a non-spreading jet. The energy injection avoids these problems, but requires $\epsilon_e > 1 (k = 2)$, spherical outflow, and $\epsilon_B < 10^{-9}$. Conclusions. In light of the extreme microphysical parameters required by the energy-injection model, we favour a jet-break scenario where $\nu_m < \nu_{sa}$ to explain the observations. This scenario gives physically meaningful microphysical parameters, and it also naturally explains the reported detection of linear and circular polarisation., Comment: Final accepted version

New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering., Comment: 6 pages, 6 figures, v2 to match published version

After the discovery of the first connection between GRBs and SNe almost two decades ago, tens of SN-like rebrightenings have been discovered and about seven solid associations have been spectroscopically confirmed to date. Using GROND optical/NIR data and Swift X-ray/UV data, we estimate the intrinsic extinction, luminosity, and evolution of three SN rebrightenings in GRB afterglow light curves at z~0.5. The SNe 2008hw, 2009nz, and 2010ma exhibit 0.80, 1.15, and 1.78 times the optical (r band) luminosity of SN 1998bw, respectively. While SN 2009nz evolves similarly to SN 1998bw, SNe 2008hw and 2010ma show earlier peak times. The quasi-bolometric light curves were corrected for the contribution of the NIR bands using data available in the literature and blackbody fits. The large luminosity of SN 2010ma (1.4x10^43 erg/s) is confirmed, while SNe 2008hw and 2009nz reached a peak luminosity closer to SN 1998bw. Physical parameters of the SN explosions, such as synthesised nickel mass, ejecta mass, and kinetic energy, are estimated using Arnett's analytic approach, which resulted in nickel masses of around 0.4-0.5 Msun. By means of the a very comprehensive data set, we found that the luminosity and the nickel mass of SNe 2008hw, 2009nz, and 2010ma resembles those of other known GRB-associated SNe. This findings strengthens previous claims of GRB-SNe being brighter than type-Ic SNe unaccompanied by GRBs., Comment: 11 pages, 9 figures, accepted for publication in Astronomy & Astrophysics, abstract abridged