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.

We compare two observations of gamma-rays before, during, and after lightning flashes initiated by upward leaders from a tower during low-altitude winter thunderstorms on the western coast of Honshu, Japan. While the two leaders appear similar, one produced a terrestrial gamma-ray flash (TGF) so bright that it paralyzed the gamma-ray detectors while it was occurring, and could be observed only via the weaker flux of neutrons created in its wake, while the other produced no detectable TGF gamma-rays at all. The ratio between the indirectly derived gamma-ray fluence for the TGF and the 95% confidence gamma-ray upper limit for the gamma-ray quiet flash is a factor of $1\times10^7$. With the only two observations of this type providing such dramatically different results -- a TGF probably as bright as those seen from space and a powerful upper limit -- we recognize that weak, sub-luminous TGFs in this situation are probably not common, and we quantify this conclusion. While the gamma-ray quiet flash appeared to have a faster leader and more powerful initial continuous current pulse than the flash that produced a TGF, the TGF-producing flash occurred during a weak gamma-ray "glow", while the gamma-ray quiet flash did not, implying a higher electric field aloft when the TGF was produced. We suggest that the field in the high-field region approached by a leader may be more important for whether a TGF is produced than the characteristics of the leader itself.

Production of runaway electron avalanches and gamma rays originating inside Martian dust storms are studied using Monte Carlo simulations. In the absence of in situ measurements, we use theoretical predictions of electric fields inside dust storms. Electrons are produced through the relativistic runaway electron avalanches process, and energetic photons are results of the bremsstrahlung scattering of the electrons with the air. Characteristic lengths of the runaway electron avalanche for different electric fields and the energy spectrum of electrons are derived and compared to their terrestrial counterparts. It is found that it is possible for Martian dust storms to develop energetic electron avalanches and produce large fluxes of gamma ray photons similar to terrestrial gamma ray flashes from Earth's thunderstorms. The phenomenon could be called Martian gamma ray flash, and due to the very thin atmosphere on Mars, it can be observed by both ground-based instruments or satellites orbiting the planet.

Terrestrial gamma-ray flashes (TGFs) are short intense flashes of gamma rays associated with lightning activity in thunderstorms. Using Monte Carlo simulations of the relativistic runaway electron avalanche (RREA) process, theoretical predictions for the temporal and spectral evolution of TGFs are compared to observations made with the Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope. Assuming a single source altitude of 15 km, a comparison of simulations to data is performed for a range of empirically chosen source electron variation time scales. The data exhibit a clear softening with increased source distance, in qualitative agreement with theoretical predictions. The simulated spectra follow this trend in the data, but tend to underestimate the observed hardness. Such a discrepancy may imply that the basic RREA model is not sufficient. Alternatively, a TGF beam that is tilted with respect to the zenith could produce an evolution with source distance that is compatible with the data. Based on these results, we propose that the source electron distributions of TGFs observed by GBM vary on time scales of at least tens of microseconds, with an upper limit of approx. 100 microseconds.

The Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope has triggered on over 300 terrestrial gamma-ray flashes (TGFs) since its launch in June 2008. With 14 detectors, GBM collects on average ~100 counts per triggered TGF, enabling unprecedented studies of the time profiles of TGFs. Here we present the first rigorous analysis of the temporal properties of a large sample of TGFs (278), including the distributions of the rise and fall times of the individual pulses and their durations. A variety of time profiles are observed with 19 of TGFs having multiple pulses separated in time and 31 clear cases of partially overlapping pulses. The effect of instrumental dead time and pulse pileup on the temporal properties are also presented. As the observed gamma ray pulse structure is representative of the electron flux at the source, TGF pulse parameters are critical to distinguish between relativistic feedback discharge and lightning leader models. We show that at least 67% of TGFs at satellite altitudes are significantly asymmetric. For the asymmetric pulses, the rise times are almost always shorter than the fall times. Those which are not are consistent with statistical fluctuations. The median rise time for asymmetric pulses is ~3 times shorter than for symmetric pulses while their fall times are comparable. The asymmetric shapes observed are consistent with the relativistic feedback discharge model when Compton scattering of photons between the source and Fermi is included, and instrumental effects are taken into account.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of high-energy radiation from thunderstorms and lightning. For a summary, we refer to the paper.

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GWs) and high-energy neutrinos (HENs). A network of GW detectors such as LIGO and Virgo can determine the direction/time of GW bursts while the IceCube and ANTARES neutrino telescopes can also provide accurate directional information for HEN events. Requiring the consistency between both, totally independent, detection channels shall enable new searches for cosmic events arriving from potential common sources, of which many extra-galactic objects.

We present direct upper limits on gravitational wave emission from the Crab pulsar using data from the first nine months of the fifth science run of the Laser Interferometer Gravitational-wave Observatory (LIGO). These limits are based on two searches. In the first we assume that the gravitational wave emission follows the observed radio timing, giving an upper limit on gravitational wave emission that beats indirect limits inferred from the spin-down and braking index of the pulsar and the energetics of the nebula. In the second we allow for a small mismatch between the gravitational and radio signal frequencies and interpret our results in the context of two possible gravitational wave emission mechanisms.

We present a coincidence search method for astronomical events using gravitational wave detectors in conjunction with other astronomical observations. We illustrate our method for the specific case of the LIGO gravitational wave detector and the IceCube neutrino detector. LIGO trigger-events and IceCube events which occur within a given time window are selected as time-coincident events. Then the spatial overlap of the reconstructed event directions is evaluated using an unbinned maximum likelihood method. Our method was tested with Monte Carlo simulations based on realistic LIGO and IceCube event distributions. We estimated a typical false alarm rate for the analysis to be 1 event per 435 years. This is significantly smaller than the false alarm rates of the ndividual detectors.

We present spectra of the inverse-Compton and extragalactic components of the high-energy gamma-radiation based on an analysis of the emission at high galactic latitudes ($\vert$b$\vert$ $\geq$ $30\deg$). We correlate the gamma-ray intensity with a model consisting of an isotropic component, a component proportional to the 408 MHz synchrotron radiation, and a component proportional to atomic hydrogen (H I) column density with different emissivities in eight galactic octants. The spectrum of gamma-radiation that is correlated with the H I column density indicates that this component originates in cosmic-ray/matter interactions. The cosmic-ray electrons which produce the 408 MHz radio continuum emission also produce gamma-radiation through inverse-Compton interactions with interstellar photons. By correlating the gamma-radiation with the 408 MHz continuum, we measure the spectrum and absolute intensity of the IC emission. The isotropic component gives us the spectrum and intensity of the extragalactic gamma-radiation. We discuss interpretation of the extragalactic spectrum.