Electron dynamics surrounding the X-line in magnetopause-type asymmetric reconnection is investigated using a two-dimensional particle-in-cell simulation. We study electron properties of three characteristic regions in the vicinity of the X-line. The fluid properties, velocity distribution functions (VDFs), and orbits are studied and cross-compared. On the magnetospheric side of the X-line, the normal electric field enhances the electron meandering motion from the magnetosheath side. The motion leads to a crescent-shaped component in the electron VDF, in agreement with recent studies. On the magnetosheath side of the X-line, the magnetic field line is so stretched in the third dimension that its curvature radius is comparable with typical electron Larmor radius. The electron motion becomes nonadiabatic, and therefore the electron idealness is no longer expected to hold. Around the middle of the outflow regions, the electron nonidealness is coincident with the region of the nonadiabatic motion. Finally, we introduce a finite-time mixing fraction (FTMF) to evaluate electron mixing. The FTMF marks the magnetospheric side of the X-line, where the nonideal energy dissipation occurs.

Particle dynamics in the electron current layer in collisionless magnetic reconnection is investigated by using a particle-in-cell simulation. Electron motion and velocity distribution functions are studied by tracking self-consistent trajectories. New classes of electron orbits are discovered: figure-eight-shaped regular orbits inside the electron jet, noncrossing regular orbits on the jet flanks, noncrossing Speiser orbits, and nongyrotropic electrons in the downstream of the jet termination region. Properties of a super-Alfv\'{e}nic outflow jet are attributed to an ensemble of electrons traveling through Speiser orbits. Noncrossing orbits are mediated by the polarization electric field near the electron current layer. The noncrossing electrons are found to be non-negligible in number density. The impact of these new orbits to electron mixing, spatial distribution of energetic electrons, and observational signatures, is presented.

Kinetic aspects of the ion current layer at the center of a reconnection outflow exhaust near the X-type region are investigated by a two-dimensional particle-in-cell (PIC) simulation. The layer consists of magnetized electrons and unmagnetized ions that carry a perpendicular electric current. The ion fluid appears to be nonideal, sub-Alfvenic, and nondissipative. The ion velocity distribution functions contain multiple populations such as global Speiser ions, local Speiser ions, and trapped ions. The particle motion of the local Speiser ions in an appropriately rotated coordinate system explains the ion fluid properties very well. The trapped ions are the first demonstration of the regular orbits in the chaotic particle dynamics [Chen and Palmadesso, J. Geophys. Res., 91, 1499 (1986)] in self-consistent PIC simulations. They would be observational signatures in the ion current layer near reconnection sites.

Signatures of the dissipation region of collisionless magnetic reconnection are investigated by the Geotail spacecraft for the 15 May 2003 event. The energy dissipation in the rest frame of the electron's bulk flow is considered in an approximate form D*_e, which is validated by a particle-in-cell simulation. The dissipation measure is directly evaluated from the {plasma moments}, the electric field, and the magnetic field. Using D*_e, a compact dissipation region is successfully detected in the vicinity of the possible X-point in Geotail data. The dissipation rate is 45 pWm**{-3}. The length of the dissipation region is estimated to 1--2 local ion inertial length. The Lorentz work W, the work rate by Lorentz force to plasmas, is also introduced. It is positive over the reconnection region and it has a peak around the pileup region away from the X-point. These new measures D*_e and W provide useful information to understand the reconnection structure.

The intermediate-frequency peaked BL Lacertae (IBL) object 3C 66A is detected during 2007 - 2008 in VHE (very high energy: E > 100 GeV) gamma-rays with the VERITAS stereoscopic array of imaging atmospheric Cherenkov telescopes. An excess of 1791 events is detected, corresponding to a significance of 21.2 standard deviations (sigma), in these observations (32.8 hours live time). The observed integral flux above 200 GeV is 6% of the Crab Nebula's flux and shows evidence for variability on the time-scale of days. The measured energy spectrum is characterized by a soft power law with photon index Gamma = 4.1 +- 0.4_stat +- 0.6_sys. The radio galaxy 3C 66B is excluded as a possible source of the VHE emission.

We present results from an intensive VERITAS monitoring campaign of the high-frequency peaked BL Lac object 1ES 1218+304 in 2008/2009. Although 1ES 1218+304 was detected previously by MAGIC and VERITAS at a persistent level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite common in many nearby blazars, the case of 1ES 1218+304 (redshift z = 0.182) is particularly interesting since it belongs to a group of blazars that exhibit unusually hard VHE spectra considering their redshifts. When correcting the measured spectra for absorption by the extragalactic background light, 1ES 1218+304 and a number of other blazars are found to have differential photon indices less than 1.5. The difficulty in modeling these hard spectral energy distributions in blazar jets has led to a range of theoretical gamma-ray emission scenarios, one of which is strongly constrained by these new VERITAS observations. We consider the implications of the observed light curve of 1ES 1218+304, which shows day scale flux variations, for shock acceleration scenarios in relativistic jets, and in particular for the viability of kiloparsec-scale jet emission scenarios.

Although Galactic cosmic rays (protons and nuclei) are widely believed to be dominantly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly a century after their discovery [1]. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size, more than 50 times the diameter of similar Galactic regions, uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density [2]. The cosmic rays produced in the formation, life, and death of their massive stars are expected to eventually produce diffuse gamma-ray emission via their interactions with interstellar gas and radiation. M 82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in gamma rays [3, 4]. Here we report the detection of >700 GeV gamma rays from M 82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core of M 82, or about 500 times the average Galactic density. This result strongly supports that cosmic-ray acceleration is tied to star formation activity, and that supernovae and massive-star winds are the dominant accelerators.

We present results from a long-term monitoring campaign on the TeV binary LSI +61 303 with VERITAS at energies above 500 GeV, and in the 2-10 keV hard X-ray bands with RXTE and Swift, sampling nine 26.5 day orbital cycles between September 2006 and February 2008. The binary was observed by VERITAS to be variable, with all integrated observations resulting in a detection at the 8.8 sigma (2006/2007) and 7.3 sigma (2007/2008) significance level for emission above 500 GeV. The source was detected during active periods with flux values ranging from 5 to 20% of the Crab Nebula, varying over the course of a single orbital cycle. Additionally, the observations conducted in the 2007-2008 observing season show marginal evidence (at the 3.6 sigma significance level) for TeV emission outside of the apastron passage of the compact object around the Be star. Contemporaneous hard X-ray observations with RXTE and Swift show large variability with flux values typically varying between 0.5 and 3.0*10^-11 ergs cm^-2 s^-1 over a single orbital cycle. The contemporaneous X-ray and TeV data are examined and it is shown that the TeV sampling is not dense enough to detect a correlation between the two bands.

The VERITAS collaboration reports the detection of very-high-energy (VHE) gamma-ray emission from the high-frequency-peaked BL Lac object 1ES 1218+304 located at a redshift of z=0.182. A gamma-ray signal was detected with a statistical significance of 10.4 standard deviations (10.4 sigma) for the observations taken during the first three months of 2007, confirming the discovery of this object made by the MAGIC collaboration. The photon spectrum between ~160 GeV and ~1.8 TeV is well described by a power law with an index of Gamma = 3.08 +/- 0.34_stat +/- 0.2_sys. The integral flux is Phi(E > 200 GeV) = (12.2 +/- 2.6) X 10^-12 cm^-2 s^-1, which corresponds to ~6% of that of the Crab Nebula. The light curve does not show any evidence for VHE flux variability. Using lower limits on the density of the extragalactic background light in the near to mid-infrared we are able to limit the range of intrinsic energy spectra for 1ES 1218+304. We show that the intrinsic photon spectrum has an index that is harder than Gamma = 2.32 +/- 0.37_stat. When including constraints from the spectra of 1ES 1101-232 and 1ES 0229+200, the spectrum of 1ES 1218+304 is likely to be harder than Gamma = 1.86 +/- 0.37_stat.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

The high-frequency-peaked BL-Lacertae object \objectname{1ES 0806+524}, at redshift z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by VERITAS between November 2006 and April 2008. These data encompass the two-, and three-telescope commissioning phases, as well as observations with the full four-telescope array. \objectname{1ES 0806+524} is detected with a statistical significance of 6.3 standard deviations from 245 excess events. Little or no measurable variability on monthly time scales is found. The photon spectrum for the period November 2007 to April 2008 can be characterized by a power law with photon index $3.6 \pm 1.0_{\mathrm{stat}} \pm 0.3_{\mathrm{sys}}$ between $\sim$300 GeV and $\sim$700 GeV. The integral flux above 300 GeV is $(2.2\pm0.5_{\mathrm{stat}}\pm0.4_{\mathrm{sys}})\times10^{-12}\:\mathrm{cm}^{2}\:\mathrm{s}^{-1}$ which corresponds to 1.8% of the Crab Nebula flux. Non contemporaneous multiwavelength observations are combined with the VHE data to produce a broadband spectral energy distribution that can be reasonably described using a synchrotron-self Compton model.

VERITAS, an Imaging Atmospheric Cherenkov Telescope (IACT) system for gammma-ray astronomy in the GeV-TeV range, has recently completed its first season of observations with a full array of four telescopes. A number of astrophysical gamma-ray sources have been detected, both galactic and extragalactic, including sources previously unknown at TeV energies. We describe the status of the array and some highlight results, and assess the technical performance, sensitivity and shower reconstruction capabilities.

We report the detection of very high-energy gamma-ray emission from the intermediate-frequency-peaked BL Lacertae object W Comae (z=0.102) by VERITAS. The source was observed between January and April 2008. A strong outburst of gamma-ray emission was measured in the middle of March, lasting for only four days. The energy spectrum measured during the two highest flare nights is fit by a power-law and is found to be very steep, with a differential photon spectral index of Gamma = 3.81 +- 0.35_stat +- 0.34_syst. The integral photon flux above 200GeV during those two nights corresponds to roughly 9% of the flux from the Crab Nebula. Quasi-simultaneous Swift observations at X-ray energies were triggered by the VERITAS observations. The spectral energy distribution of the flare data can be described by synchrotron-self-Compton (SSC) or external-Compton (EC) leptonic jet models, with the latter offering a more natural set of parameters to fit the data.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy gamma-rays. The system was observed over several orbital cycles (between September 2006 and February 2007) with the VERITAS array of imaging air-Cherenkov telescopes. A signal of gamma-rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of Gamma=2.40+-0.16_stat+-0.2_sys and a flux above 300 GeV corresponding to 15-20% of the flux from the Crab Nebula.

The multiwavelength observation of the nearby radio galaxy M87 provides a unique opportunity to study in detail processes occurring in Active Galactic Nuclei from radio waves to TeV gamma-rays. Here we report the detection of gamma-ray emission above 250 GeV from M87 in spring 2007 with the VERITAS atmospheric Cherenkov telescope array and discuss its correlation with the X-ray emission. The gamma-ray emission is measured to be point-like with an intrinsic source radius less than 4.5 arcmin. The differential energy spectrum is fitted well by a power-law function: dPhi/dE=(7.4+-1.3_{stat}+-1.5_{sys})(E/TeV)^{-2.31+-0.17_{stat}+-0.2_{sys}} 10^{-9}m^{-2}s^{-1}TeV^{-1}. We show strong evidence for a year-scale correlation between the gamma-ray flux reported by TeV experiments and the X-ray emission measured by the ASM/RXTE observatory, and discuss the possible short-time-scale variability. These results imply that the gamma-ray emission from M87 is more likely associated with the core of the galaxy than with other bright X-ray features in the jet.

We present observations of the dwarf galaxies Draco and Ursa Minor, the local group galaxies M32 and M33, and the globular cluster M15 conducted with the Whipple 10m gamma-ray telescope to search for the gamma-ray signature of self-annihilating weakly interacting massive particles (WIMPs) which may constitute astrophysical dark matter (DM). We review the motivations for selecting these sources based on their unique astrophysical environments and report the results of the data analysis which produced upper limits on excess rate of gamma rays for each source. We consider models for the DM distribution in each source based on the available observational constraints and discuss possible scenarios for the enhancement of the gamma-ray luminosity. Limits on the thermally averaged product of the total self-annihilation cross section and velocity of the WIMP, <\sigma v>, are derived using conservative estimates for the magnitude of the astrophysical contribution to the gamma-ray flux. Although these limits do not constrain predictions from the currently favored theoretical models of supersymmetry (SUSY), future observations with VERITAS will probe a larger region of the WIMP parameter phase space, <\sigma v> and WIMP particle mass (m_\chi).

VERITAS is a new atmospheric Cherenkov imaging telescope array to detect very high energy gamma rays above 100 GeV. The array is located in southern Arizona, USA, at an altitude of 1268m above sea level. The array consists of four 12-m telescopes of Davies-Cotton design and structurally resembling the Whipple 10-m telescope. The four focal plane instruments are equipped with high-resolution (499 pixels) fast photo-multiplier-tube (PMT) cameras covering a 3.5 degree field of view with 0.15 degree pixel separation. Light concentrators reduce the dead-space between PMTs to 25% and shield the PMTs from ambient light. The PMTs are connected to high-speed preamplifiers allowing operation at modest anode current and giving good single photoelectron peaks in situ. Electronics in the focus box provides real-time monitoring of the anode currents for each pixel and ambient environmental conditions. A charge injection subsystem installed in the focus box allows daytime testing of the trigger and data acquisition system by injecting pulses of variable amplitude and length directly into the photomultiplier preamplifiers. A brief description of the full VERITAS focal plane instrument is given in this paper.

Gamma-ray burst (GRB) observations at very high energies (VHE, E > 100 GeV) can impose tight constraints on some GRB emission models. Many GRB afterglow models predict a VHE component similar to that seen in blazars and plerions, in which the GRB spectral energy distribution has a double-peaked shape extending into the VHE regime. VHE emission coincident with delayed X-ray flare emission has also been predicted. GRB follow-up observations have had high priority in the observing program at the Whipple 10m Gamma-ray Telescope and GRBs will continue to be high priority targets as the next generation observatory, VERITAS, comes on-line. Upper limits on the VHE emission, at late times (>~4 hours), from seven GRBs observed with the Whipple Telescope are reported here.

We report on observations of the sky region around the unidentified TeV gamma-ray source TeV J2032+4130 carried out with the Whipple Observatory 10 m atmospheric Cherenkov telescope for a total of 65.5 hrs between 2003 and 2005. The standard two-dimensional analysis developed by the Whipple collaboration for a stand-alone telescope reveals an excess in the field of view at a pre-trials significance level of 6.1 standard deviations. The measured position of this excess is alpha(2000) =20 h 32 m 27 s, delta(2000) = 41 deg 39 min 17 s. The estimated integral flux for this gamma-ray source is about 8% of the Crab-Nebula flux. The data are consistent with a point-like source. Here we present a detailed description of the standard two-dimensional analysis technique used for the analysis of data taken with the Whipple Observatory 10 m telescope and the results for the TeV J2032+4130 campaign. We include a short discussion of the physical mechanisms that may be responsible for the observed gamma-ray emission, based on possible association with known astrophysical objects, in particular Cygnus OB2.

Q^2 evolution equations are important not only for describing hadron reactions in accelerator experiments but also for investigating ultrahigh-energy cosmic rays. The standard ones are called DGLAP evolution equations, which are integrodifferential equations. There are methods for solving the Q^2 evolution equations for parton-distribution and fragmentation functions. Because the equations cannot be solved analytically, various methods have been developed for the numerical solution. We compare brute-force, Laguerre-polynomial, and Mellin-transformation methods particularly by focusing on the numerical accuracy and computational efficiency. An efficient solution could be used, for example, in the studies of a top-down scenario for the ultrahigh-energy cosmic rays.

The BL Lac object H1426+428 was recently detected as a high energy gamma-ray source by the VERITAS collaboration (Horan et al. 2002). We have reanalyzed the 2001 portion of the data used in the detection in order to examine the spectrum of H1426+428 above 250 GeV. We find that the time-averaged spectrum agrees with a power law of the shape dF/dE = 10^(-7.31 +- 0.15(stat) +- 0.16(syst)) x E^(-3.50 +- 0.35(stat) +- 0.05(syst)) m^(-2)s^(-1)TeV^(-1) The statistical evidence from our data for emission above 2.5 TeV is 2.6 sigma. With 95% c.l., the integral flux of H1426+428 above 2.5 TeV is larger than 3% of the corresponding flux from the Crab Nebula. The spectrum is consistent with the (non-contemporaneous) measurement by Aharonian et al. (2002) both in shape and in normalization. Below 800 GeV, the data clearly favours a spectrum steeper than that of any other TeV Blazar observed so far indicating a difference in the processes involved either at the source or in the intervening space.