With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders.

We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ \mu$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models.

This paper presents a detailed description of various subsystems of CALLISTO solar radio spectrograph installed at the USO/PRL. In the front-end system, a log periodic dipole antenna (LPDA) is designed for the frequency range of 40-900 MHz. In this paper LPDA design, its modifications, and simulation results are presented. We also present some initial observations taken by CALLISTO at Udaipur.

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere.

Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and an array of plastic scintillators. We compare the reconstruction of cosmic-ray energy using radio and particle techniques and show that they are in good agreement. We also present the radiation energy of air showers measured at LOFAR. This value scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LOFAR energy scale and that of the Pierre Auger Observatory, which are shown to be in agreement. This technique demonstrates the potential to cross-calibrate the energy scales of different experiments using radiation energy.

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.

We study a complex GOES M1.1 circular ribbon flare and related pre-flare activity on 26 January 2015 [SOL26-01-2015] in solar active region NOAA 12268. This flare activity was observed by the AIA on board SDO and the RHESSI. The examination of photospheric magnetograms during the extended period, prior to the event, suggests the successive development of a so-called 'anemone' type magnetic configuration. NLFFF extrapolation reveals a fan-spine magnetic configuration with the presence of a coronal null-point. We found that the pre-flare activity in the active region starts ~15 min prior to the main flare in the form of localized bright patches at two locations. A comparison of locations and spatial structures of the pre-flare activity with magnetic configuration of the corresponding region suggests onset of magnetic reconnection at the null-point along with the low-atmosphere magnetic reconnection caused by the emergence and the cancellation of the magnetic flux. The main flare of M1.1 class is characterized by the formation of a well-developed circular ribbon along with a region of remote brightening. Remarkably, a set of relatively compact parallel ribbons formed inside the periphery of the circular ribbon which developed lateral to the brightest part of the circular ribbon. During the peak phase of the flare, a coronal jet is observed at the north-east edge of the circular ribbon which suggests interchange reconnection between large-scale field lines and low-lying closed field lines. Our investigation suggests a combination of two distinct processes in which ongoing pre-flare null-point reconnection gets further intensified as the confined eruption along with jet activity proceeded from within the circular ribbon region which results to the formation of inner parallel ribbons and corresponding post-reconnection arcade.

In this paper, we present a comprehensive study of the evolutionary phases of a major M6.6 long duration event (LDE) with special emphasize on its pre-flare phase. The event occurred in NOAA 12371 on 2015 June 22. A remarkable aspect of the event was an active pre-flare phase lasting for about an hour during which a hot EUV coronal channel was in build-up stage and displayed co-spatial hard X-ray (HXR) emission up to energies of 25 keV. As such, this is the first evidence of HXR coronal channel. The coronal magnetic field configuration based on NLFFF modeling clearly exhibited a magnetic flux rope (MFR) oriented along the polarity inversion line (PIL) and co-spatial with the coronal channel. We observed significant changes in the AR's photospheric magnetic field during an extended period of $\approx$42 hours in the form of rotation of sunspots, moving magnetic features, and flux cancellation along the PIL. Prior to the flare onset, the MFR underwent a slow rise phase ($\approx$14 km s$^{-1}$) for $\approx$12 min which we attribute to the faster build-up and activation of the MFR by tether-cutting reconnection occurring at multiple locations along the MFR itself. The sudden transition in the kinematic evolution of the MFR from the phase of slow to fast rise ($\approx$109 km s$^{-1}$ with acceleration $\approx$110 m s$^{-2}$) precisely divides the pre-flare and impulsive phase of the flare, which points toward the feedback process between the early dynamics of the eruption and the strength of the flare magnetic reconnection.

We report direct-detection constraints on light dark matter particles interacting with electrons. The results are based on a method that exploits the extremely low levels of leakage current of the DAMIC detector at SNOLAB of 2-6$\times$10$^{-22}$ A cm$^{-2}$. We evaluate the charge distribution of pixels that collect $<10~\rm{e^-}$ for contributions beyond the leakage current that may be attributed to dark matter interactions. Constraints are placed on so-far unexplored parameter space for dark matter masses between 0.6 and 100 MeV$c^{-2}$. We also present new constraints on hidden-photon dark matter with masses in the range $1.2$-$30$ eV$c^{-2}$.

In this article, we study the origin of precursor flare activity and investigate its role towards triggering the eruption of a flux rope which resulted into a dual-peak M-class flare (SOL2015-06-21T02:36) in the active region NOAA 12371. The flare evolved in two distinct phases with peak flux levels of M2.1 and M2.6 at an interval of $\approx$54 min. The active region exhibited striking moving magnetic features (MMFs) along with sunspot rotation. Non-linear force free field (NLFFF) modelling of the active region corona reveals a magnetic flux rope along the polarity inversion line in the trailing sunspot group which is observationally manifested by the co-spatial structures of an active region filament and a hot channel identified in the 304 and 94 \AA\ images, respectively, from the Atmospheric Imaging Assembly (AIA). The active region underwent a prolonged phase of flux enhancement followed by a relatively shorter period of flux cancellation prior to the onset of the flare which led to the build up and activation of the flux rope. Extreme ultra-violet (EUV) images reveal localised and structured pre-flare emission, from the region of MMFs, adjacent to the location of the main flare. Our analysis reveals strong, localised regions of photospheric currents of opposite polarities at the precursor location, thereby making the region susceptible to small-scale magnetic reconnection. Precursor reconnection activity from this location most likely induced a slipping reconnetion towards the northern leg of the hot channel which led to the destabilization of the flux rope. The application of magnetic virial theorem suggests that there was an overall growth of magnetic free energy in the active region during the prolonged pre-flare phase which decayed rapidly after the hot channel eruption and its successful transformation into a halo coronal mass ejection (CME).

This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate.

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}$.

We present a multi-wavelength analysis of the eruption of a hot coronal channel associated with an X1.0 flare (SOL2013-10-28T02:03) from the active region NOAA 11875 by combining observations from AIA/SDO, HMI/SDO, RHESSI, and HiRAS. EUV images at high coronal temperatures indicated the presence of a hot channel at the core of the active region from the early pre-flare phase evidencing the pre-existence of a quasi-stable magnetic flux rope. The hot channel underwent an activation phase after a localized and prolonged pre-flare event occurring adjacent to one of its footpoints. Subsequently, the flux rope continued to rise slowly for $\approx$16 min during which soft X-ray flux gradually built-up characterizing a distinct precursor phase. The flux rope transitioned from the state of slow rise to the eruptive motion with the onset of the impulsive phase of the X1.0 flare. The eruptive expansion of the hot channel is accompanied by a series of type III radio bursts in association with impulsive rise of strong hard X-ray non-thermal emissions that included explicit hard X-ray sources of energies up to $\approx$50 keV from the coronal loops and $\approx$100 keV from their footpoint locations. Our study contains evidence that pre-flare activity occurring within the spatial extent of a stable flux rope can destabilize it toward eruption. Moreover, sudden transition of the flux rope from the state of slow rise to fast acceleration precisely bifurcated the precursor and the impulsive phases of the flare which points toward a feedback relationship between early CME dynamics and the strength of the large-scale magnetic reconnection.

The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope.

Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR.

The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement.

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.

In this paper, we present a multi-wavelength analysis of two X-class solar eruptive flares of classes X2.2 and X9.3 that occurred in the sigmoidal active region NOAA 12673 on 2017 September 6, by combining observations of Atmospheric Imaging Assembly and Helioseismic Magnetic Imager instruments on board the Solar Dynamics Observatory. On the day of the reported activity, the photospheric structure of the active region displayed a very complex network of $\delta$-sunspots that gave rise to the formation of a coronal sigmoid observed in the hot EUV channels. Both X-class flares initiated from the core of the sigmoid sequentially within an interval of $\sim$3 hours and progressed as a single "sigmoid--to--arcade" event. Differential emission measure analysis reveals strong heating of plasma at the core of the active region right from the pre-flare phase which further intensified and spatially expanded during each event. The identification of a pre-existing magnetic null by non-force-free-field modeling of the coronal magnetic fields at the location of early flare brightenings and remote faint ribbon-like structures during the pre-flare phase, which were magnetically connected with the core region, provide support for the breakout model of solar eruption. The magnetic extrapolations also reveal flux rope structures prior to both flares which are subsequently supported by the observations of the eruption of hot EUV channels. The second X-class flare diverged from the standard flare scenario in the evolution of two sets of flare ribbons, that are spatially well separated, providing firm evidence of magnetic reconnections at two coronal heights.

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017)

In this paper, we present unique observations of a blowout coronal jet that possibly triggered a two-ribbon confined C1.2 flare in a bipolar solar active region NOAA 12615 on 2016 December 5. The jet activity initiates at chromospheric/transition-region heights with a small brightening that eventually grows in a larger volume with well developed standard morphological jet features, viz., base and spire. The spire widens up with a collimated eruption of cool and hot plasma components, observed in the 304 and 94 A channels of AIA, respectively. The speed of the plasma ejection, which forms the jet's spire, was higher for the hot component (~200 km/s) than the cooler one (~130 km/s). The NLFF model of coronal fields at pre- and post-jet phases successfully reveal opening of previously closed magnetic field lines with a rather inclined/low-lying jet structure. The peak phase of the jet emission is followed by the development of a two-ribbon flare that shows coronal loop emission in HXRs up to ~25 keV energy. The coronal magnetic fields rooted at the location of EUV flare ribbons, derived from the NLFF model, demonstrate the pre-flare phase to exhibit an "X-type" configuration while the magnetic fields at the post-flare phase are more or less parallel oriented. The comparisons of multi-wavelength measurements with the magnetic field extrapolations suggest that the jet activity likely triggered the two-ribbon flare by perturbing the field in the interior of the active region.

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.

The PICASSO dark matter search experiment operated an array of 32 superheated droplet detectors containing 3.0 kg of C$_{4}$F$_{10}$ and collected an exposure of 231.4 kgd at SNOLAB between March 2012 and January 2014. We report on the final results of this experiment which includes for the first time the complete data set and improved analysis techniques including \mbox{acoustic} localization to allow fiducialization and removal of higher activity regions within the detectors. No signal consistent with dark matter was observed. We set limits for spin-dependent interactions on protons of $\sigma_p^{SD}$~=~1.32~$\times$~10$^{-2}$~pb (90\%~C.L.) at a WIMP mass of 20 GeV/c$^{2}$. In the spin-independent sector we exclude cross sections larger than $\sigma_p^{SI}$~=~4.86~$\times$~10$^{-5 }$~pb~(90\% C.L.) in the region around 7 GeV/c$^{2}$. The pioneering efforts of the PICASSO experiment have paved the way forward for a next generation detector incorporating much of this technology and experience into larger mass bubble chambers.

We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed.

As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.

For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype.

We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. A significant increase in the number of detected UHECR is expected to be achieved by employing Earth's moon as detector, and search for short radio pulses that are emitted when a particle interacts in the lunar rock. Observation of these short pulses with current and future radio telescopes also allows to search for the even lower fluxes of neutrinos with energies above $10^{22}$ eV, that are predicted in certain Grand-Unifying-Theories (GUTs), and e.g. models for super-heavy dark matter (SHDM). In this contribution we present the initial design for such a search with the LOFAR radio telescope.

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.

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.

We consider the communication channel given by a fiber optical transmission line. We develop a method to perturbatively calculate the information capacity of a nonlinear channel, given the corresponding evolution equation. Using this technique, we compute the decrease of the channel capacity to the leading order in the perturbative parameter for fiber optics communication systems.