We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

Gamma-ray bursts (GRBs) are spectacularly energetic events, with the potential to inform on the early universe and its evolution, once their redshifts are known. Unfortunately, determining redshifts is a painstaking procedure requiring detailed follow-up multi-wavelength observations often involving various astronomical facilities, which have to be rapidly pointed at these serendipitous events. Here we use Machine Learning algorithms to infer redshifts from a collection of observed temporal and spectral features of GRBs. We obtained a very high correlation coefficient ($0.96$) between the inferred and the observed redshifts, and a small dispersion (with a mean square error of $0.003$) in the test set. The addition of plateau afterglow parameters improves the predictions by $61.4\%$ compared to previous results. The GRB luminosity function and cumulative density rate evolutions, obtained from predicted and observed redshift are in excellent agreement indicating that GRBs are effective distance indicators and a reliable step for the cosmic distance ladder.

We explore the general question of correlations among different waveband luminosities in a flux-limited multiband observational data set. Such correlations, often observed for astronomical sources, may either be intrinsic or induced by the redshift evolution of the luminosities and the data truncation due to the flux limits. We first address this question analytically. We then use simulated flux-limited data with three different known intrinsic luminosity correlations, and prescribed luminosity functions and evolution similar to the ones expected for quasars. We explore how the intrinsic nature of luminosity correlations can be deduced, including exploring the efficacy of partial correlation analysis with redshift binning in determining whether luminosity correlations are intrinsic and finding the form of the intrinsic correlation. By applying methods that we have developed in recent works, we show that we can recover the true cosmological evolution of the luminosity functions and the intrinsic correlations between the luminosities. Finally, we demonstrate the methods for determining intrinsic luminosity correlations on actual observed samples of quasars with mid-infrared, radio, and optical fluxes and redshifts, finding that the luminosity-luminosity correlation is significantly stronger between mid-infrared and optical than that between radio and optical luminosities, supporting the canonical jet-launching and heating model of active galaxies.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

A knowledge of the particle escape time from the acceleration regions of many space and astrophysical sources is of critical importance in the analysis of emission signatures produced by these particles and in the determination of the acceleration and transport mechanisms at work. This paper addresses this general problem, in particular in solar flares, where in addition to scattering by turbulence, the magnetic field convergence from the acceleration region towards its boundaries also influences the particle escape. We test an (approximate) analytic relation between escape and scattering times, and the field convergence rate, based on the work of Malyshkin and Kulsrud (2001), valid for both strong and weak diffusion limits and isotropic pitch angle distribution of the injected particles, with a numerical model of particle transport. To this end, a kinetic Fokker-Planck transport model of particles is solved with a stochastic differential equation scheme assuming different initial pitch angle distributions. This approach enables further insights into the phase-space dynamics of the transport process, which would otherwise not be accessible. We find that in general the numerical results agree well with the analytic equation for the isotropic case, however, there are significant differences weak diffusion regime for non-isotopic cases, especially for distributions beamed along the magnetic field lines. The results are important in the interpretation of observations of energetic particles in solar flares, and other similar space and astrophysical acceleration sites and for the determination of acceleration-transport coefficients, commonly used in Fokker-Planck type kinetic equations.

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

The Fermi gamma-ray Space Telescope (Fermi) has detected hard X-ray (HXR) and gamma-ray photons from three flares, which according to \stereo occurred in active regions behind the limb of the Sun as delineated by near Earth instruments. For two of these flares \r has provided HXR images with sources located just above the limb, presumably from the loop top (LT) region of a relatively large loop. Fermi-Gamma-ray Burst Monitor has detected HXRs and gamma-rays, and RSTN has detected microwaves emissions with similar light curves. This paper presents a quantitative analysis of these multi-wavelength observations assuming that HXRs and microwaves are produced by electrons accelerated at the LT source, with emphasize on the importance of the proper treatment of escape of the particles from the acceleration-source region and the trans-relativistic nature of the analysis. The observed spectra are used to determine the magnetic field and relativistic electron spectra. It is found that a simple power-law in momentum (with cut off above a few 100 MeV) agrees with all observations, but in energy space a broken power law spectrum (steepening at rest mass energy) may be required. It is also shown that the production of the $>100$ MeV photons detected by The Fermi-Large Area Telescope at the LT source would require more energy compared to photospheric emission. These energies are smaller than that required for electrons, so that the possibility that all the emissions originate in the LT cannot be ruled out on energetic grounds. However, the differences in the light curves and emission centroids of HXRs and $>100$ MeV gamma-rays favour a different source for the latter.

The discovery of the Fermi bubbles---a huge bilobular structure seen in GeV gamma-rays above and below the Galactic center---implies the presence of a large reservoir of high energy particles at $\sim 10 \, \text{kpc}$ from the disk. The absence of evidence for a strong shock coinciding with the edge of the bubbles, and constraints from multi-wavelength observations point towards stochastic acceleration by turbulence as a likely mechanism of acceleration. We have investigated the time-dependent acceleration of electrons in a large-scale outflow from the Galactic centre. For the first time, we present a detailed numerical solution of the particle kinetic equation that includes the acceleration, transport and relevant energy loss processes. We also take into account the addition of shock acceleration of electrons at the bubble's blast wave. Fitting to the observed spectrum and surface brightness distribution of the bubbles allows determining the transport coefficients, thereby shedding light on the origin of the Fermi bubbles.

The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6\sigma$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2.

We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (CMEs) and strong solar energetic particle events. We present updated localizations of the >100 MeV photon emission, hard X-ray (HXR)and EUV images, and broadband spectra from 10 keV to 10 GeV, as well as microwave spectra. We also provide a comparison of the BTL flares detected by Fermi-LAT with three on-disk flares and present a study of some of the significant quantities of these flares as an attempt to better understand the acceleration mechanisms at work during these occulted flares. We interpret the HXR emission to be due to electron bremsstrahlung from a coronal thin-target loop top with the accelerated electron spectra steepening at semirelativistic energies. The >100 MeV gamma-rays are best described by a pion-decay model resulting from the interaction of protons (and other ions) in a thick-target photospheric source. The protons are believed to have been accelerated (to energies >10 GeV) in the CME environment and precipitate down to the photosphere from the downstream side of the CME shock and landed on the front side of the Sun, away from the original flare site and the HXR emission.

Flares close to the solar limb, where the footpoints are occulted, can reveal the spectrum and structure of the coronal loop-top source in X-rays. We aim at studying the properties of the corresponding energetic electrons near their acceleration site, without footpoint contamination. To this end, a statistical study of partially occulted flares observed with RHESSI is presented here, covering a large part of solar cycles 23 and 24. We perform a detailed spectra, imaging and light curve analysis for 116 flares and include contextual observations from SDO and STEREO when available, providing further insights into flare emission that was previously not accessible. We find that most spectra are fitted well with a thermal component plus a broken power-law, non-thermal component. A thin-target kappa distribution model gives satisfactory fits after the addition of a thermal component. X-rays imaging reveals small spatial separation between the thermal and non-thermal components, except for a few flares with a richer coronal source structure. A comprehensive light curve analysis shows a very good correlation between the derivative of the soft X-ray flux (from GOES) and the hard X-rays for a substantial number of flares, indicative of the Neupert effect. The results confirm that non-thermal particles are accelerated in the corona and estimated timescales support the validity of a thin-target scenario with similar magnitudes of thermal and non-thermal energy fluxes.

At hard X-ray energies, the bright footpoint emission from solar flare loops often prevents a detailed analysis of the weaker loop-top source morphology due to the limited dynamic range available for X-ray imaging. Here, we study the X1.3 April 25, 2014 flare with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). This partially occulted limb flare allows the analysis of the loop-top emission in isolation. We present results on the flare light curve at different energies, the source morphology from X-ray imaging and a detailed spectral analysis of the different source components by imaging spectroscopy. The loop-top source, a likely site of particle acceleration, shows a clear composition of different emission components. The results indicate the opportunities that detailed imaging of hard X-rays can provide to learn about particle acceleration, transport and heating processes in solar flares.

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in MgII h&k, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DST/IBIS) in H$\alpha$ 6563 \AA\ and CaII 8542 \AA, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multi-threaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated MgII line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of $\approx$800 km below the transition region can turn the simulated MgII line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Observations relating the characteristics of electrons seen near Earth (SEPs) and those producing flare radiation show that in certain (prompt) events the origin of both population appears to be the flare site, which show strong correlation between the number and spectral index of SEP and hard X-ray radiating electrons, but in others(delayed), which are associated with fast CMEs, this relation is complex and SEPs tend to be harder. Prompt event spectral relation disagrees with that expected in thick or thin target models. We show that using a a more accurate treatment of the transport of the accelerated electrons to the footpoints and to the Earth can account for this discrepancy. Our results are consistent with those found by Chen and Petrosian (2013) for two flares using non-parametric inversion methods, according to which we have weak diffusion conditions, and trapping mediated by magnetic field convergence. The weaker correlations and harder spectra of delayed events can come about by re-acceleration of electrons in the CME shock environment. We describe under what conditions such a hardening can be achieved. Using this (acceleration at the flare and re-acceleration in the CME) scenario we show that we can describe the similar dichotomy that exists between the so called impulsive, highly enriched ($^3$He and heavy ions) and softer SEP events, and stronger more gradual SEP events with near normal ionic abundances and harder spectra. These methods can be used to distinguish the acceleration mechanisms and to constrain their characteristics.

We investigate heating and acceleration of protons from a thermal gas with a generic diffusion and acceleration model, and subject to Coulomb scattering and energy loss, as was carried out in Petrosian & East (2008) for electrons. As protons gain energy their loss to electrons becomes important. Thus, we need to solve the coupled proton-electron kinetic equation. We numerically solve the coupled Fokker-Plank equations and computes the time evolution of the spectra of both particles. We show that this can lead to a quasi-thermal component plus a high energy nonthermal tail. We determine the evolution of nonthermal tail and the quasi-thermal component. The results may be used to explore the possibility of inverse bremsstrahlung radiation as a source of hard X-ray emissions from hot sources such as solar flares, accretion disk coronas and the intracluster medium of galaxy clusters. We find that emergence of nonthermal protons is accompanied by excessive heating of the entire plasma, unless the turbulence needed for scattering and acceleration is steeper than Kolmogorov and the acceleration parameters, the duration of the acceleration, and/or the initial distributions are significantly fine-tuned. These results severely constraint the feasibility of nonthermal inverse bremsstrahlung process producing hard X-ray emissions. However the nonthermal tail may be the seed particles for further re-acceleration to relativistic energies, say by a shock. In the Appendix we present some tests of the integrity of the algorithm used and present a new formula for the energy loss rate due to inelastic proton-proton interactions.

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $\gamma$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3\deg that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $\gamma$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle \sigma v \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $\gamma$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

Solar flares involve complex processes that are coupled and span a wide range of temporal, spatial, and energy scales. Modeling such processes self-consistently has been a challenge in the past. Here we present results from simulations that couple particle kinetics with hydrodynamics of the atmospheric plasma. We combine the Stanford unified Fokker-Planck code that models particle acceleration and transport with the RADYN hydrodynamic code that models the atmospheric response to collisional heating by accelerated electrons through detailed radiative transfer calculations. We perform simulations using two different electron spectra, one an {\it ad hoc} power law and the other predicted by the model of stochastic acceleration by turbulence or plasma waves. Surprisingly, the later model, even with energy flux $\ll 10^{10}$ erg s$^{-1}$ cm$^{-2}$, can cause "explosive" chromospheric evaporation and drive stronger up- and downflows (and hydrodynamic shocks). This is partly because our acceleration model, like many others, produces a spectrum consisting of a quasi-thermal component plus a power-law tail. We synthesize emission line profiles covering different heights in the lower atmosphere, including H$\alpha$ 6563 \AA, HeII 304 \AA, CaII K 3934 \AA\ and SiIV 1393 \AA. One interesting result is the unusual high temperature (up to a few $10^5$ K) of the formation site of HeII 304 \AA, which is expected due to photoionization-recombination under flare conditions, compared to those in the quiet Sun dominated by collisional excitation. When compared with observations, our results can constrain the properties of non-thermal electrons and thus the poorly understood particle acceleration mechanism.

Fermi-LAT >30 MeV observations have increased the number of detected solar flares by almost a factor of 10 with respect to previous space observations. These sample both the impulsive and long duration phases of GOES M and X class flares. Of particular interest is the recent detections of three solar flares whose position behind the limb was confirmed by the STEREO-B spacecraft. While gamma-ray emission up to tens of MeV resulting from proton interactions has been detected before from occulted solar flares, the significance of these particular events lies in the fact that these are the first detections of >100 MeV gamma-ray emission from footpoint-occulted flares. We will present the Fermi-LAT, RHESSI and STEREO observations of these flares and discuss the various emission scenarios for these sources and implications for the particle acceleration mechanisms.

We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $\sigma$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $\rho=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $\rho=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists.

We report the first detection of >100 MeV gamma rays associated with a behind-the-limb solar flare, which presents a unique opportunity to probe the underlying physics of high-energy flare emission and particle acceleration. On 2013 October 11 a GOES M1.5 class solar flare occurred ~ 9.9 degrees behind the solar limb as observed by STEREO-B. RHESSI observed hard X-ray emission above the limb, most likely from the flare loop-top, as the footpoints were occulted. Surprisingly, the Fermi Large Area Telescope (LAT) detected >100 MeV gamma-rays for ~30 minutes with energies up to GeV. The LAT emission centroid is consistent with the RHESSI hard X-ray source, but its uncertainty does not constrain the source to be located there. The gamma-ray spectra can be adequately described by bremsstrahlung radiation from relativistic electrons having a relatively hard power-law spectrum with a high-energy exponential cutoff, or by the decay of pions produced by accelerated protons and ions with an isotropic pitch-angle distribution and a power-law spectrum with a number index of ~3.8. We show that high optical depths rule out the gamma rays originating from the flare site and a high-corona trap model requires very unusual conditions, so a scenario in which some of the particles accelerated by the CME shock travel to the visible side of the Sun to produce the observed gamma rays may be at work.

Solar flares involve impulsive energy release, which results in enhanced radiation over a broad spectral range and a wide range of heights. In particular, line emission from the chromosphere can provide critical diagnostics of plasma heating processes. Thus, a direct comparison between high-resolution spectroscopic observations and advanced numerical modeling results could be extremely valuable, but has not yet been attempted. In this paper, we present in this paper such a self-consistent investigation of an M3.0 flare observed by the Dunn Solar Telescope's Interferometric Bi-dimensional Spectrometer (IBIS) on 2011 September 24 which we have modeled with the radiative hydrodynamic code RADYN. We obtained images and spectra of the flaring region with IBIS in H$\alpha$ 6563 \AA\ and Ca II 8542 \AA, and with the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in X-rays. The latter observations were used to infer the non-thermal electron population, which was passed to RADYN to simulate the atmospheric response to electron collisional heating. We then synthesized spectral lines and compared their shapes and intensities to those observed by IBIS and found a general agreement. In particular, the synthetic Ca II 8542 \AA\ profile fits well to the observed profile, while the synthetic H$\alpha$ profile is fainter in the core than for the observation. This indicates that H$\alpha$ emission is more responsive to the non-thermal electron flux than the Ca II 8542 \AA\ emission. We suggest that it is necessary to refine the energy input and other processes is necessary to resolve this discrepancy.

Gamma-ray bursts (GRBs) by virtue of their high luminosities can be detected up to very high redshifts and therefore can be excellent probes of the early universe. This task is hampered by the fact that most of their characteristics have a broad range so that we first need to obtain an accurate description of the distribution of these characteristics, and specially, their cosmological evolution. We use a sample of about 200 \swift long GRBs with known redshift to determine the luminosity and formation rate evolutions and the general shape of the luminosity function. In contrast to most other forward fitting methods of treating this problem we use the Efron Petrosian methods which allow a non-parametric determination of above quantities. We find a relatively strong luminosity evolution, a luminosity function that can be fitted to a broken power law, and an unusually high rate of formation rate at low redshifts, a rate more than one order of magnitude higher than the star formation rate (SFR). On the other hand, our results seem to agree with the almost constant SFR in redshifts 1 to 3 and the decline above this redshift.

During its first six years of operation, the Fermi Large Area Telescope (LAT) has detected >30 MeV gamma-ray emission from more than 40 solar flares, nearly a factor of 10 more than those detected by EGRET. These include detections of impulsive and sustained emissions, extending up to 20 hours in the case of the 2012 March 7 X-class flares. We will present an overview of solar flare detections with LAT, highlighting recent results and surprising features, including the detection of >100 MeV emission associated with flares located behind the limb. Such flares may shed new light on the relationship between the sites of particle acceleration and gamma-ray emission.

Active galactic nuclei jets are detected via their radio and/or gamma-ray emissions while the accretion disks are detected by their optical and UV radiation. Observations of the radio and optical luminosities show a strong correlation between the two luminosities. However, part of this correlation is due to the redshift or distances of the sources that enter in calculating the luminosities from the observed fluxes and part of it could be due to the differences in the cosmological evolution of luminosities. Thus, the determination of the intrinsic correlations between the luminosities is not straightforward. It is affected by the observational selection effects and other factors that truncate the data, sometimes in a complex manner (e.g. Antonucci (2011) and Pavildou et al. (2010)). In this paper we describe methods that allow us to determine the evolution of the radio and optical luminosities, and determine the true intrinsic correlation between the two luminosities. We find a much weaker correlation than observed and sub-linear relations between the luminosities. This has a significant implication for jet and accretion disk physics.

We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $\kappa$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $\kappa$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $\kappa$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $\kappa$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $\kappa$-distribution with $10 < \kappa_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $\kappa$-distribution with $\kappa_e < 2$ at quasi-perpendicular shocks. For such $\kappa$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays.

We have developed an inversion method for determination of the characteristics of the acceleration mechanism directly and non-parametrically from observations, in contrast to the usual forward fitting of parametric model variables to observations. In two recent papers (Petrosian & Chen 2010, Chen & Petrosian 2013), we demonstrate the efficacy of this inversion method by its application to acceleration of electrons in solar flares based on stochastic acceleration by turbulence. Here we explore its application for determining the characteristics of shock acceleration in supernova remnants (SNRs) based on the electron spectra deduced from the observed nonthermal radiation from SNRs and the spectrum of the cosmic ray electrons observed near the Earth. These spectra are related by the process of escape of the electrons from SNRs and energy loss during their transport in the galaxy. Thus, these observations allow us to determine spectral characteristics of the momentum and pitch angle diffusion coefficients, which play crucial roles in both direct acceleration by turbulence and in high Mach number shocks. Assuming that the average electron spectrum deduced from a few well known SNRs is representative of those in the solar neighborhood we find interesting discrepancies between our deduced forms for these coefficients and those expected from well known wave-particle interactions. This may indicate that the standard assumptions made in treatment of shock acceleration need revision. In particular, the escape of particles from SNRs may be more complex than generally assumed.

We present the redshift evolutions and distributions of the gamma-ray luminosity and photon spectral index of flat spectrum radio quasar (FSRQ) type blazars, using non-parametric methods to obtain the evolutions and distributions directly from the data. The sample we use for analysis consists of almost all FSRQs observed with a greater than approximately 7 sigma detection threshold in the first year catalog of the Fermi Gamma-ray Space Telescope's Large Area Telescope, with redshfits as determined from optical spectroscopy by Shaw et al. We find that FSQRs undergo rapid gamma-ray luminosity evolution, but negligible photon index evolution, with redshift. With these evolutions accounted for we determine the density evolution and luminosity function of FSRQs, and calculate their total contribution to the extragalactic gamma-ray background radiation, resolved and unresolved, which is found to be 16(+10/-4)%, in agreement with previous studies.

Following our recent paper (Petrosian & Chen 2010), we have developed an inversion method to determine the basic characteristics of the particle acceleration mechanism directly and non-parametrically from observations under the leaky box framework. In the above paper, we demonstrated this method for obtaining the energy dependence of the escape time. Here, by converting the Fokker-Planck equation to its integral form, we derive the energy dependences of the energy diffusion coefficient and direct acceleration rate for stochastic acceleration in terms of the accelerated and escaping particle spectra. Combining the regularized inversion method of Piana et al. 2007 and our procedure, we relate the acceleration characteristics in solar flares directly to the count visibility data from RHESSI. We determine the timescales for electron escape, pitch angle scattering, energy diffusion, and direct acceleration at the loop top acceleration region for two intense solar flares based on the regularized electron flux spectral images. The X3.9 class event shows dramatically different energy dependences for the acceleration and scattering timescales, while the M2.1 class event shows a milder difference. The M2.1 class event could be consistent with the stochastic acceleration model with a very steep turbulence spectrum. A likely explanation of the X3.9 class event could be that the escape of electrons from the acceleration region is not governed by a random walk process, but instead is affected by magnetic mirroring, in which the scattering time is proportional to the escape time and has an energy dependence similar to the energy diffusion time.

Gamma-ray bursts (GRBs), which have been observed up to redshifts z approx 9.5 can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the rest frame time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic and describe how we recover it and how this can be used to constrain physical models of the plateau emission, whose origin is still unknown. The present result could help clarifing the debated issue about the nature of the plateau emission.

Gamma-ray bursts (GRBs), which have been observed up to redshifts $z \approx 9.5$ can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the \underline{rest frame} time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to estimate the luminosity and time evolution and to correct for this effects to determine the intrinsic nature of this correlation.

The Fermi Large Area Telescope (LAT) is the most sensitive instrument ever deployed in space for observing gamma-ray emission >100 MeV. This sensitivity has enabled the LAT to detect gamma-ray emission from the Sun during quiescent periods from pions produced by cosmic-ray protons interacting in the solar atmosphere and from cosmic-ray electrons interacting with solar optical photons. The LAT has detected high-energy gamma-ray emission associated with GOES M-class and X-class X-ray flares accompanied by coronal mass ejections and solar energetic particle events. In a number of cases, LAT has detected gamma rays with energies up to several hundreds of MeV during the impulsive phase and gamma rays up to GeV energies sustained for several hours after the impulsive phase. This presentation focuses on observations in the impulsive emission phase in solar flares, including the modest GOES M2-class flare at SOL2010-06-12T0057 and more recent detections, such as the bright X-class flares of March 2012.

Where particle acceleration and plasma heating take place in relation to magnetic reconnection is a fundamental question for solar flares. We report analysis of an M7.7 flare on 2012 July 19 observed by SDO/AIA and RHESSI. Bi-directional outflows in forms of plasmoid ejections and contracting cusp-shaped loops originate between an erupting flux rope and underlying flare loops at speeds of typically 200-300 km/s up to 1050 km/s. These outflows are associated with spatially separated double coronal X-ray sources with centroid separation decreasing with energy. The highest temperature is located near the nonthermal X-ray loop-top source well below the original heights of contracting cusps near the inferred reconnection site. These observations suggest that the primary loci of particle acceleration and plasma heating are in the reconnection outflow regions, rather than the reconnection site itself. In addition, there is an initial ascent of the X-ray and EUV loop-top source prior to its recently recognized descent, which we ascribe to the interplay among multiple processes including the upward development of reconnection and the downward contractions of reconnected loops. The impulsive phase onset is delayed by 10 minutes from the start of the descent, but coincides with the rapid speed increases of the upward plasmoids, the individual loop shrinkages, and the overall loop-top descent, suggestive of an intimate relation of the energy release rate and reconnection outflow speed.

We determine the radio and optical luminosity evolutions and the true distribution of the radio loudness parameter R, defined as the ratio of the radio to optical luminosity, for a set of more than 5000 quasars combining SDSS optical and FIRST radio data. We apply the method of Efron and Petrosian to access the intrinsic distribution parameters, taking into account the truncations and correlations inherent in the data. We find that the population exhibits strong positive evolution with redshift in both wavebands, with somewhat greater radio evolution than optical. With the luminosity evolutions accounted for, we determine the density evolutions and local radio and optical luminosity functions. The intrinsic distribution of the radio loudness parameter R is found to be quite different than the observed one, and is smooth with no evidence of a bi-modality in radio loudness. The results we find are in general agreement with the previous analysis of Singal et al. 2011 which used POSS-I optical and FIRST radio data.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

We present a determination of the distributions of photon spectral index and gamma-ray flux - the so called LogN-LogS relation - for the 352 blazars detected with a greater than approximately seven sigma detection threshold and located above +/- 20 degrees Galactic latitude by the Large Area Telescope of the Fermi Gamma-ray Space Telescope in its first year catalog. Because the flux detection threshold depends on the photon index, the observed raw distributions do not provide the true LogN-LogS counts or the true distribution of the photon index. We use the non-parametric methods developed by Efron and Petrosian to reconstruct the intrinsic distributions from the observed ones which account for the data truncations introduced by observational bias and includes the effects of the possible correlation between the two variables. We demonstrate the robustness of our procedures using a simulated data set of blazars and then apply these to the real data and find that for the population as a whole the intrinsic flux distribution can be represented by a broken power law with high and low indexes of -2.37 +/- 0.13 and -1.70 +/- 0.26, respectively, and the intrinsic photon index distribution can be represented by a Gaussian with mean of 2.41 +/- 0.13 and width of 0.25 +/- 0.03. We also find the intrinsic distributions for the sub-populations of BL Lac and FSRQs type blazars separately. We then calculate the contribution of Fermi blazars to the diffuse extragalactic gamma-ray background radiation. Under the assumption that the flux distribution of blazars continues to arbitrarily low fluxes, we calculate the best fit contribution of all blazars to the total extragalactic gamma-ray output to be 60%, with a large uncertainty.

The subject of this paper is stochastic acceleration by plasma turbulence, a process akin to the original model proposed by Fermi. We review the relative merits of different acceleration models, in particular the so called first order Fermi acceleration by shocks and second order Fermi by stochastic processes, and point out that plasma waves or turbulence play an important role in all mechanisms of acceleration. Thus, stochastic acceleration by turbulence is active in most situations. We also show that it is the most efficient mechanism of acceleration of relatively cool non relativistic thermal background plasma particles. In addition, it can preferentially accelerate electrons relative to protons as is needed in many astrophysical radiating sources, where usually there are no indications of presence of shocks. We also point out that a hybrid acceleration mechanism consisting of initial acceleration by turbulence of background particles followed by a second stage acceleration by a shock has many attractive features. It is demonstrated that the above scenarios can account for many signatures of the accelerated electrons, protons and other ions, in particular $^3$He and $^4$He, seen directly as Solar Energetic Particles and through the radiation they produce in solar flares.

[Abridged]We present analysis of a pair of unusually energetic coronal hard X-ray (HXR) sources detected by RHESSI during the impulsive phase of an X3.9 class solar flare on 2003 November 3, which simultaneously shows two intense footpoint (FP) sources. A distinct loop top (LT) coronal source is detected up to ~150 keV and a second (upper) coronal source up to ~80 keV. These photon energies are much higher than commonly observed in coronal sources and pose grave modeling challenges. The LT source in general appears higher in altitude with increasing energy and exhibits a more limited motion compared to the expansion of the thermal loop. The high energy LT source shows an impulsive time profile and its nonthermal power law spectrum exhibits soft-hard-soft evolution during the impulsive phase, similar to the FP sources. The upper coronal source exhibits an opposite spatial gradient and a similar spectral slope compared to the LT source. These properties are consistent with the model of stochastic acceleration of electrons by plasma waves or turbulence. However, the LT and FP spectral index difference (varying from ~0-1) is much smaller than commonly measured and than that expected from a simple stochastic acceleration model. Additional confinement or trapping mechanisms of high energy electrons in the corona are required. Comprehensive modeling including both kinetic effects and the macroscopic flare structure may shed light on this behavior. These results highlight the importance of imaging spectroscopic observations of the LT and FP sources up to high energies in understanding electron acceleration in solar flares. Finally, we show that the electrons producing the upper coronal HXR source may very likely be responsible for the type III radio bursts at the decimetric/metric wavelength observed during the impulsive phase of this flare.

We calculate simultaneously the radio and optical luminosity evolutions of quasars, and the distribution in radio loudness R defined as the ratio of radio and optical luminosities, using a flux limited data set containing 636 quasars with radio and optical fluxes from White et al. We first note that when dealing with multivariate data it is imperative to first determine the true correlations among the variables, not those introduced by the observational selection effects, before obtaining the individual distributions of the variables. We use the methods developed by Efron and Petrosian which are designed to obtain unbiased correlations, distributions, and evolution with redshift from a data set truncated due to observational biases. It is found that the population of quasars exhibits strong positive correlation between the radio and optical luminosities. With this correlation, whether intrinsic or observationally induced accounted for, we find that there is a strong luminosity evolution with redshift in both wavebands, with significantly higher radio than optical evolution. We also construct the local radio and optical luminosity functions and the density evolution. Finally, we consider the distribution of the radio loudness parameter R obtained from careful treatment of the selection effects and luminosity evolutions with that obtained from the raw data without such considerations. We find a significant difference between the two distributions and no clear sign of bi-modality in the true distribution for the range of R values considered. Our results indicate therefore, somewhat surprisingly, that there is no critical switch in the efficiency of the production of disk outflows/jets between very radio quiet and very radio loud quasars, but rather a smooth transition. Also, this efficiency seems higher for the high-redshift and more luminous sources in the considered sample.

We examine the effects of time dilation on the temporal profiles of gamma-ray burst (GRB) pulses. By using prescriptions for the shape and evolution of prompt gamma-ray spectra, we can generate a simulated population of single pulsed GRBs at a variety of redshifts and observe how their light curves would appear to a gamma-ray detector here on Earth. We find that the observer frame duration of individual pulses does not increase as a function of redshift as one would expect from the cosmological expansion of a Friedman-Lemaitre-Robertson-Walker Universe. In fact, the duration of individual pulses is seen to decrease as their signal-to-noise decreases with increasing redshift, as only the brightest portion of a high redshift GRB's light curve is accessible to the detector. The results of our simulation are consistent with the fact that a systematic broadening of GRB durations as a function of redshift has not materialized in either the Swift or Fermi detected GRBs with known redshift. We show that this fundamental duration bias implies that the measured durations and associated Eiso estimates for GRBs detected near an instrument's detection threshold should be considered lower limits to their true values. We conclude by predicting that the average peak-to-peak time for a large number of multi-pulsed GRBs as a function of redshift may eventually provide the evidence for time dilation that has so far eluded detection.

We review basic theoretical concepts in particle acceleration, with particular emphasis on processes likely to occur in regions of magnetic reconnection. Several new developments are discussed, including detailed studies of reconnection in three-dimensional magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell approaches are used and results compared. While these studies show considerable promise in accounting for the various observational manifestations of solar flares, they are limited by a number of factors, mostly relating to available computational power. Not the least of these issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated particles themselves on the environment in which they are accelerated. A brief prognosis for future advancement is offered.

We investigate different scenarios for the origin of the extragalactic radio background. The surface brightness of the background, as reported by the ARCADE 2 collaboration, is several times higher than that which would result from currently observed radio sources. We consider contributions to the background from diffuse synchrotron emission from clusters and the intergalactic medium, previously unrecognized flux from low surface brightness regions of radio sources, and faint point sources below the flux limit of existing surveys. By examining radio source counts available in the literature, we conclude that most of the radio background is produced by radio point sources that dominate at sub microJy fluxes. We show that a truly diffuse background produced by electrons far from galaxes is ruled out because such energetic electrons would overproduce the obserevd X-ray/gamma-ray background through inverse Compton scattering of the other photon fields. Unrecognized flux from low surface brightness regions of extended radio sources, or moderate flux sources missed entirely by radio source count surveys, cannot explain the bulk of the observed background, but may contribute as much as 10 per cent. We consider both radio supernovae and radio quiet quasars as candidate sources for the background, and show that both fail to produce it at the observed level because of insufficient number of objects and total flux, although radio quiet quasars contribute at the level of at least a few percent. We conclude that if the radio background is at the level reported, a majority of the total surface brightness would have to be produced by ordinary starforming galaxies above redshift 1 characterized by an evolving radio far-infrared correlation, which changes toward the radio loud with redshift.

The model of stochastic acceleration of particles by turbulence has been successful in explaining many observed features of solar flares. Here we demonstrate a new method to obtain the accelerated electron spectrum and important acceleration model parameters from the high resolution hard X-ray observations provided by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). In our model, electrons accelerated at or very near the loop top produce thin target bremsstrahlung emission there and then escape downward producing thick target emission at the loop footpoints. Based on the electron flux spectral images obtained by the regularized inversion of the RHESSI count visibilities, we derive several important parameters for the acceleration model. We apply this procedure to the 2003 November 03 solar flare, which shows a loop top source up to 100--150 keV in hard X-ray with a relatively flat spectrum in addition to two footpoint sources. The results imply presence of strong scattering and a high density of turbulence energy with a steep spectrum in the acceleration region.

The Fermi Gamma-ray Space Telescope observed the bright and long GRB090902B, lying at a redshift of z = 1.822. Together the Large Area Telescope (LAT) and the Gamma-ray Burst Monitor (GBM) cover the spectral range from 8 keV to >300 GeV is covered. Here we show that the prompt burst spectrum is consistent with emission from the jet photosphere combined with non-thermal emission described by a single power-law with photon index -1.9. The photosphere gives rise to a strong quasi-blackbody spectrum which is somewhat broader than a single Planck function and has a characteristic temperature of ~290keV. We derive the photospheric radius Rph = (1.1 \pm 0.3) x 10^12 Y^{1/4} cm and the bulk Lorentz factor of the flow, which is found to vary by a factor of two and has a maximal value of Gamma = 750 Y^{1/4}. Here Y is the ratio between the total fireball energy and the energy emitted in the gamma-rays. We find that during the first quarter of the prompt phase the photospheric emission dominates, which explains the delayed onset of the observed flux in the LAT compared to the GBM. We model the photospheric emission with a multi-color blackbody and its shape indicates that the photospheric radius increases at higher latitudes. We interpret the broad band emission as synchrotron emission at R ~ 4x10^15 cm. Our analysis emphasize the importance of having high temporal resolution when performing spectral analysis on GRBs, since there is strong spectral evolution.

Recent observations of cosmic ray electrons from several instruments have revealed various degrees of deviation in the measured electron energy distribution from a simple power-law, in a form of an excess around TeV energies. An even more prominent deviation has been observed in the fraction of cosmic ray positrons around 100 GeV energies. In this paper we show that the observed excesses in the electron spectrum may be easily re-produced without invoking any unusual sources other than the general diffuse Galactic components of cosmic rays. The primary physical effect involved is the Klein-Nishina suppression of the electron cooling rate around TeV energies. With a very reasonable choice of the model parameters characterizing the local interstellar medium, we can reproduce the most recent observations by Fermi and HESS experiments. We also find that high positron fraction increasing with energy, as claimed by the PAMELA experiment, cannot be explained in our model with the conservative set of the model parameters. We are able, however, to reproduce the PAMELA results assuming high values of the starlight and interstellar gas densities, which would be more appropriate for vicinities of supernova remnants. A possible solution to this problem may be that cosmic rays undergo most of their interactions near their sources due to the efficient trapping in the far upstream of supernova shocks by self-generated, cosmic ray-driven turbulence.

In recent years there has been considerable activity in using gamma-ray bursts as cosmological probes for determining global cosmological parameters complementing results from type Ia supernovae and other methods. This requires a characteristics of the source to be a standard candle. We show that contrary to earlier indications the accumulated data speak against this possibility. Another method would be to use correlation between a distance dependent and a distance independent variable to measure distance and determine cosmological parameters as is done using Cepheid variables and to some extent Type Ia supernovae. Many papers have dealt with the use of so called Amati relation, first predicted by Lloyd, Petrosian and Mallozzi, or the Ghirlanda relation for this purpose. We have argued that these procedure involve many unjustified assumptions which if not true could invalidate the results. In particular, we point out that many evolutionary effects can affect the final outcome. In particular, we demonstrate that the existing data from Swift and other earlier satellites show that the gamma-ray burst may have undergone luminosity evolution. Similar evolution may be present for other variables such as the peak photon energy of the total radiated energy. Another out come of our analysis is determination of the luminosity function and the comoving rate evolution of gamma-ray bursts which does not seem to agree with the cosmic star formation rate. We caution however, that the above result are preliminary and includes primarily the effect of detection threshold. Other selection effects, perhaps less important than this, are also known to be present and must be accounted for. We intend to address these issues in future publications.

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares. We present here self-consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the NRL flux tube code. We calculated the collisional heating rate from the particle transport code, which is more accurate than those based on approximate analytical solutions. We used a realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a nonthermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (vs. chromosphere), a larger downward conductive flux, and thus a stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced X-ray emission.

It is well-recognized that the presence of magnetic fields will lead to anisotropic energy cascade and dissipation of astrophysical turbulence. With the diffusion approximation and linear dissipation rates, we study the cascade and damping of Alfv\'en-cyclotron fluctuations in solar plasmas numerically. For an isotropic case the steady-state turbulence spectra are nearly isotropic in the inertial range and can be fitted by a single power-law function with a spectral index of -3/2, similar to the Iroshnikov-Kraichnan phenomenology. Beyond the MHD regime the kinetic effects make the spectrum softer at higher wavenumbers. In the dissipation range the turbulence spectrum cuts off at the wavenumber, where the damping rate becomes comparable to the cascade rate, and the cutoff wavenumber changes with the wave propagation direction. The angle averaged turbulence spectrum of the isotropic model resembles a broken power-law. Taking into account the Doppler effects, the model naturally reproduces the broken power-law turbulence spectra observed in the solar wind and predicts that a higher break frequency always comes along with a softer dissipation range spectrum that may be caused by the increase of the turbulence intensity, the reciprocal of the plasma \beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by detailed comparisons with more accurate observations.

We present a detailed imaging and spectroscopic study of the conjugate hard X-ray (HXR) footpoints (FPs) observed with RHESSI in the 2003 October 29 X10 flare. The double FPs first move toward and then away from each other, mainly parallel and perpendicular to the magnetic neutral line, respectively. The transition of these two phases of FP unshearing motions coincides with the direction reversal of the motion of the loop-top (LT) source, and with the minima of the estimated loop length and LT height. The FPs show temporal correlations between HXR flux, spectral index, and magnetic field strength. The HXR flux exponentially correlates with the magnetic field strength, which also anti-correlates with the spectral index before the second HXR peak's maximum, suggesting that particle acceleration sensitively depends on the magnetic field strength and/or reconnection rate. Asymmetries are observed between the FPs: on average, the eastern FP is 2.2 times brighter in HXR flux and 1.8 times weaker in magnetic field strength, and moves 2.8 times faster away from the neutral line than the western FP; the estimated coronal column density to the eastern FP from the LT source is 1.7 times smaller. The two FPs have marginally different spectral indexes. The eastern-to-western FP HXR flux ratio and magnetic field strength ratio are anti-correlated only before the second HXR peak's maximum. Neither magnetic mirroring nor column density alone can explain the totality of these observations, but their combination, together with other transport effects, might provide a full explanation. We have also developed novel techniques to remove particle contamination from HXR counts and to estimate effects of pulse pileup in imaging spectroscopy, which can be applied to other RHESSI flares in similar circumstances.

Solar Energetic Particles (SEPs) show a rich variety of spectra and relative abundances of many ionic species and their isotopes. A long standing puzzle has been the extreme enrichments of 3He ions. The most extreme enrichments are observed in low fluence, the so-called impulsive, events which are believed to be produced at the flare site in the solar corona with little scattering and acceleration during transport to the Earth. In two earlier papers (Liu et al. 2004 and 2006) we showed how such extreme enrichments can result in the model developed by Petrosian and Liu (2004), where ions are accelerated stochastically by plasma waves or turbulence. In this paper we address the relative distributions of the fluences of 3He and 4He ions presented by Ho et al. (2005) which show that while the distribution of 4He fluence like many other extensive characteristics of solar flare, is fairly broad, the 3He fluence is limited to a narrow range. Moreover, the ratio of the fluences shows a strong correlation with the 4He fluence. One of the predictions of our model was presence of steep variation of the fluence ratio with the level of turbulence or the rate of acceleration. We show here that this feature of the model can reproduce the observed distribution of the fluences with very few free parameters. The primary reason for the success of the model in both fronts is because fully ionized 3He ion, with its unique charge to mass ratio, can resonantly interact with more plasma modes and accelerate more readily than 4He. Essentially in most flares, all background 3He ions are accelerated to few MeV/nucleon range, while this happens for 4He ions only in very strong events. A much smaller fraction of 4He ions reach such energies in weaker events.

We discuss possible mechanisms of acceleration of particles in solar flares and show that turbulence plays an important role in all the mechanism. It is also argued that stochastic particle acceleration by turbulent plasma waves is the most likely mechanism for production of the high energy electrons and ions responsible for observed radiative signatures of solar flares and for solar energetic particle or SEPs, and that the predictions of this model agrees well with many past and recent high spectral and temporal observations of solar flares. It is shown that, in addition, the model explains many features of SEPs that accompany flares. In particular we show that it can successfully explain the observed extreme enhancement, relative to photospheric values, of $^3$He ions and the relative spectra of $^3$He and $^4$He. It has also the potential of explaining the relative abundances of most ions including the increasing enhancements of heavy ions with ion mass or mass-to-charge ratio.

We address the problem of particle acceleration in solar flares by fast modes which may be excited during the reconnection and undergo cascade and are subjected to damping. We extend the calculations beyond quasilinear approximation and compare the acceleration and scattering by transit time damping and gyroresonance interactions. We find that the acceleration is dominated by the so called transit time damping mechanism. We estimate the total energy transferred into particles, and show that our approach provides sufficiently accurate results We compare this rate with energy loss rate. Scattering by fast modes appears to be sufficient to prevent the protons from escaping the system during the acceleration. Confinement of electrons, on the other hand, requires the existence of plasma waves. Electrons can be accelerated to GeV energies through the process described here for solar flare conditions.

Here we investigate some aspects of stochastic acceleration of ultrarelativistic electrons by magnetic turbulence. In particular, we discuss the steady-state energy spectra of particles undergoing momentum diffusion due to resonant interactions with turbulent MHD modes, taking rigorously into account direct energy losses connected with different radiative cooling processes. For the magnetic turbulence we assume a given power spectrum of the type $W(k) \propto k^{-q}$. In contrast to the previous approaches, however, we assume a finite range of turbulent wavevectors $k$, consider a variety of turbulence spectral indexes $1 =< q =< 2$, and concentrate on the case of a very inefficient particle escape from the acceleration site. We find that for different cooling and injection conditions, stochastic acceleration processes tend to establish a modified ultrarelativistic Maxwellian distribution of radiating particles, with the high-energy exponential cut-off shaped by the interplay between cooling and acceleration rates. For example, if the timescale for the dominant radiative process scales with the electron momentum as $\propto p^r$, the resulting electron energy distribution is of the form $n_e(p) \propto p^2 exp[ - (1 / a) (p / p_eq)^a]$, where $a = 2-q-r$, and $p_eq$ is the equilibrium momentum defined by the balance between stochastic acceleration and energy losses timescales. We also discuss in more detail the synchrotron and inverse-Compton emission spectra produced by such an electron energy distribution, taking into account Klein-Nishina effects. We point out that the curvature of the high frequency segments of these spectra, even though being produced by the same population of electrons, may be substantially different between the synchrotron and inverse-Compton components.

In this paper we investigate the feasibility of bremsstrahlung radiation from `nonthermal' electrons as a source of hard X-rays from the intracluster medium of clusters of galaxies. With an exact treatment of the Coulomb collisions in a Fokker-Planck analysis of the electron distribution we find that the severe difficulties with lifetimes of `nonthermal' particles found earlier by Petrosian (2001) using a cold target model remain problematic. We then address possible acceleration of background electrons into a nonthermal tail. We assume a simplified but generic acceleration rate and determine the expected evolution of an initially Maxwellian distribution of electrons. We find that strong nonthermal components arise only for rapid rate of acceleration which also heats up the entire plasma. These results confirm the conclusion that if the observed `nonthermal' excesses are due to some process accelerating the background thermal electrons this process must be short lived.

With the diffusion approximation, we study the cascade and damping of Alfv\'{e}n-cyclotron fluctuations in solar plasmas numerically. Motivated by wave-wave couplings and nonlinear effects, we test several forms of the diffusion tensor. For a general locally anisotropic and inhomogeneous diffusion tensor in the wave vector space, the turbulence spectrum in the inertial range can be fitted with power-laws with the power-law index varying with the wave propagation direction. For several locally isotropic but inhomogeneous diffusion coefficients, the steady-state turbulence spectra are nearly isotropic in the absence of damping and can be fitted by a single power-law function. However, the energy flux is strongly polarized due to the inhomogeneity that leads to an anisotropic cascade. Including the anisotropic thermal damping, the turbulence spectrum cuts off at the wave numbers, where the damping rates become comparable to the cascade rates. The combined anisotropic effects of cascade and damping make this cutoff wave number dependent on the wave propagation direction, and the propagation direction integrated turbulence spectrum resembles a broken power-law, which cuts off at the maximum of the cutoff wave numbers or the $^4$He cyclotron frequency. Taking into account the Doppler effects, the model can naturally reproduce the broken power-law wave spectra observed in the solar wind and predicts that a higher break frequency is aways accompanied with a greater spectral index change that may be caused by the increase of the Alfv\'{e}n Mach number, the reciprocal of the plasma beta, and/or the angle between the solar wind velocity and the mean magnetic field. These predictions can be tested by future observations.

The Warm-Hot Intergalactic Medium (WHIM) is thought to contribute about 40-50 % to the baryonic budget at the present evolution stage of the universe. The observed large scale structure is likely to be due to gravitational growth of density fluctuations in the post-inflation era. The evolving cosmic web is governed by non-linear gravitational growth of the initially weak density fluctuations in the dark energy dominated cosmology. Non-linear structure formation, accretion and merging processes, star forming and AGN activity produce gas shocks in the WHIM. Shock waves are converting a fraction of the gravitation power to thermal and non-thermal emission of baryonic/leptonic matter. They provide the most likely way to power the luminous matter in the WHIM. The plasma shocks in the WHIM are expected to be collisionless. Collisionless shocks produce a highly non-equilibrium state with anisotropic temperatures and a large differences in ion and electron temperatures. We discuss the ion and electron heating by the collisionless shocks and then review the plasma processes responsible for the Coulomb equilibration and collisional ionisation equilibrium of oxygen ions in the WHIM. MHD-turbulence produced by the strong collisionless shocks could provide a sizeable non-thermal contribution to the observed Doppler parameter of the UV line spectra of the WHIM.

In this paper we review the possible radiation mechanisms for the observed non-thermal emission in clusters of galaxies, with a primary focus on the radio and hard X-ray emission. We show that the difficulty with the non-thermal, non-relativistic Bremsstrahlung model for the hard X-ray emission, first pointed out by Petrosian (2001) using a cold target approximation, is somewhat alleviated when one treats the problem more exactly by including the fact that the background plasma particle energies are on average a factor of 10 below the energy of the non-thermal particles. This increases the lifetime of the non-thermal particles, and as a result decreases the extreme energy requirement, but at most by a factor of three. We then review the synchrotron and so-called inverse Compton emission by relativistic electrons, which when compared with observations can constrain the value of the magnetic field and energy of relativistic electrons. This model requires a low value of the magnetic field which is far from the equipartition value. We briefly review the possibilities of gamma-ray emission and prospects for GLAST observations. We also present a toy model of the non-thermal electron spectra that are produced by the acceleration mechanisms discussed in an accompanying paper.

We present the work of an international team at the International Space Science Institute (ISSI) in Bern that worked together to review the current observational and theoretical status of the non-virialised X-ray emission components in clusters of galaxies. The subject is important for the study of large-scale hierarchical structure formation and to shed light on the "missing baryon" problem. The topics of the team work include thermal emission and absorption from the warm-hot intergalactic medium, non-thermal X-ray emission in clusters of galaxies, physical processes and chemical enrichment of this medium and clusters of galaxies, and the relationship between all these processes. One of the main goals of the team is to write and discuss a series of review papers on this subject. These reviews are intended as introductory text and reference for scientists wishing to work actively in this field. The team consists of sixteen experts in observations, theory and numerical simulations.

We review the possible mechanisms for production of non-thermal electrons which are responsible for non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We briefly review acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We consider two scenarios for production of non-thermal radiation. The first is hard X-ray emission due to non-thermal Bremsstrahlung by nonrelativistic particles. Non-thermal tails are produced by accelerating electrons from the background plasma with an initial Maxwellian distribution. However, these tails are accompanied by significant heating and they are present for a short time of <10^6 yr, which is also the time that the tail will be thermalised. Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.

We present data analysis and interpretation of an M1.4-class flare observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) on April 30, 2002. This event, with its footpoints occulted by the solar limb, exhibits a rarely observed, but theoretically expected, double-source structure in the corona. The two coronal sources, observed over the 6-30 keV range, appear at different altitudes and show energy-dependent structures with the higher-energy emission being closer together. Spectral analysis implies that the emission at higher energies in the inner region between the two sources is mainly nonthermal, while the emission at lower energies in the outer region is primarily thermal. The two sources are both visible for about 12 minutes and have similar light curves and power-law spectra above about 20 keV. These observations suggest that the magnetic reconnection site lies between the two sources. Bi-directional outflows of the released energy in the form of turbulence and/or particles from the reconnection site can be the source of the observed radiation. The spatially resolved thermal emission below about 15 keV, on the other hand, indicates that the lower source has a larger emission measure but a lower temperature than the upper source. This is likely the result of the differences in the magnetic field and plasma density of the two sources.

Stochastic acceleration of electrons interacting resonantly with a turbulent magnetic field in a small accretion torus appears to be the likely mechanism responsible for much of Sagittarius A*'s millimeter and shorter wavelength spectrum. The longer wavelength radiation is produced at larger radii by electrons either diffusing from smaller scales or accelerated in situ. An important prediction of this model is the ejection of a significant flux of relativistic protons from a magnetic-field-dominated acceleration site into the wind-shocked medium surrounding the black hole. Recently, several air Cerenkov telescopes, notably HESS, have detected TeV emission from the Galactic center, with characteristics hinting at a p-p-induced pion decay process for the \gamma-ray emission. Given (1) the size of this acceleration region measured in the radio band and (2) the wind-injected ISM mapped with Chandra using the diffuse X-rays, it is feasible to test the idea that protons accelerated within \~20 Schwarzschild radii of the black hole produce the TeV emission farther out. We show a fraction of TeV protons scattering about once within ~3 pc of Sagittarius A* and the proton power (~10^37 erg s^-1) produced in concert with the 7 mm radio emission matches the TeV luminosity well. This model explains why the TeV source does not vary on a timescale of a year or less. The particle cascade generated by the p-p scatterings also produces bremsstrahlung, inverse Compton, and synchrotron emission at longer wavelengths from secondary particles. We compare these with current measurements and demonstrate that GLAST will detect this source during its one-year all-sky survey.

The near-IR and X-ray flares in Sagittarius A* are believed to be produced by relativistic electrons via synchrotron and synchrotron self-Comptonization, respectively. These electrons are likely energized by turbulent plasma waves through second-order Fermi acceleration that, in combination with the radiative cooling processes, produces a relativistic Maxwellian distribution in the steady state. This model has four principal parameters, namely the magnetic field B, the electron density n and temperature \gamma_c mec2, and the size of the flare region R. In the context of stochastic acceleration, the quantities Rn^1/2B and \gamma_cRn should remain nearly constant in time. Therefore, simultaneous spectroscopic observations in the NIR and X-ray bands can readily test the model, which, if proven to be valid, may be used to determine the evolution of the plasma properties during an eruptive event with spectroscopic observations in either band or simultaneous flux density measurements in both bands. The formulae can be applied to other isolated or confined systems, where electrons are accelerated to relativistic energies by plasma wave turbulence and produce most of the emission via synchrotron processes.

Evidence for non-thermal activity in clusters of galaxies is well established from radio observations of synchrotron emission by relativistic electrons. New windows in the Extreme Ultraviolet and Hard X-ray ranges have provided for more powerful tools for the investigation of this phenomenon. Detection of hard X-rays in the 20 to 100 keV range have been reported from several clusters of galaxies, notably from Coma and others. Based on these earlier observations we identified the relatively high redshift cluster 1E0657-56 (also known as RX J0658-5557) as a good candidate for hard X-ray observations. This cluster, also known as the bullet cluster, has many other interesting and unusual features, most notably that it is undergoing a merger, clearly visible in the X-ray images. Here we present results from a successful RXTE observations of this cluster. We summarize past observations and their theoretical interpretation which guided us in the selection process. We describe the new observations and present the constraints we can set on the flux and spectrum of the hard X-rays. Finally we discuss the constraints one can set on the characteristics of accelerated electrons which produce the hard X-rays and the radio radiation.

Flares in Sagittarius A* are produced by hot plasmas within a few Schwarzschild radii of the supermassive black hole at the Galactic center. The recent detection of a correlation between the spectral index and flux during a near infrared (NIR) flare provides a means to conduct detailed investigations of the plasma heating and radiation processes. We study the evolution of the electron distribution function under the influence of a turbulent magnetic field in a hot collisionless plasma. The magnetic field, presumably generated through instabilities in the accretion flow, can both heat the plasma via resonant wave-particle coupling and cool the electrons via radiation. The electron distribution can generally be approximated as relativistic Maxwellian. To account for the observed correlation, we find that the magnetic field needs to be anti-correlated with the electron ''temperature''. NIR and X-ray light curves are produced for a cooling and a heating phase. The model predicts simultaneous flare activity in the NIR and X-ray bands, which can be compared with observations. These results can be applied to MHD simulations to study the radiative characteristics of collisionless plasmas, especially accretion flows in low-luminosity AGNs.

We present analyses of the spatial and spectral evolution of hard X-ray emission observed by {\it RHESSI} during the impulsive phase of an M1.7 flare on 2003 November 13. In general, as expected, the loop top (LT) source dominates at low energies while the footpoint (FP) sources dominate the high energy emission. At intermediate energies, both the LT and FPs may be seen, but during certain intervals emission from the legs of the loop dominates, in contrast to the commonly observed LT and FP emission. The hard X-ray emission tends to rise above the FPs and eventually merge into a single LT source. This evolution starts first at low energies and proceeds to higher energies. The spectrum of the resultant LT source becomes more and more dominated by a thermal component with an increasing emission measure as the flare proceeds. The soft and hard X-rays show a Neupert-type behavior. With a non-thermal bremsstrahlung model the brightness profile along the loop is used to determine the density profile and its evolution, which reveals a gradual increase of the gas density in the loop. These results are evidence for chromospheric evaporation and are consistent with the qualitative features of hydrodynamic simulations of this phenomenon. However, some observed source morphology and its evolution cannot be accounted for by previous simulations. Therefore simulations with more realistic physical conditions are required to explain the results and the particle acceleration and plasma heating processes.

(Abridged) We describe the cascade of plasma waves or turbulence injected, presumably by reconnection, at scales comparable to the size of a solar flare loop to scales comparable to particle gyroradii, and evaluate their damping by various mechanisms. We show that the classical viscous damping is unimportant for magnetically dominated or low beta plasmas and the primary damping mechanism is the collisionless damping by the background particles. We show that the damping rate is proportional to the total random momentum density of the particles. For solar flare conditions this means that in most flares, except the very large ones, the damping is dominated by thermal background electrons. For large flares one requires acceleration of essentially all background electrons into a nonthermal distribution so that the accelerated electrons can be important in the damping of the waves. In general, damping by thermal or nonthermal protons is negligible compared to that of electrons except for quasi-perpendicular propagating waves or for rare proton dominated flares with strong nuclear gamma-ray line emission. Using the rate for damping we determine the critical scale below which the damping becomes important and the spectrum of the turbulence steepens. This critical scale, however, has strong dependence on the angle of propagation with respect to the magnetic field direction. The waves can cascade down to very small scales, such as the gyroradii of the particles at small angles (quasi-parallel propagation) and possibly near 90 degree (quasi-perpendicular propagation) giving rise to a highly anisotropic spectral distribution.

We present a study of the spatial and spectral evolution of the loop-top (LT) sources in a sample of 6 flares near the solar limb observed by {\it RHESSI}. A distinct coronal source, which we identify as the LT source, was seen in each of these flares from the early ``pre-heating'' phase through the late decay phase. Spectral analyses reveal an evident steep power-law component in the pre-heating and impulsive phases, suggesting that the particle acceleration starts upon the onset of the flares. In the late decay phase the LT source has a thermal spectrum and appears to be confined within a small region near the top of the flare loop, and does not spread throughout the loop, as is observed at lower energies. The total energy of this source decreases usually faster than expected from the radiative cooling but much slower than that due to the classical Spitzer conductive cooling along the flare loop. These results indicate the presence of a distinct LT region, where the thermal conductivity is suppressed significantly and/or there is a continuous energy input. We suggest that plasma wave turbulence could play important roles in both heating the plasma and suppressing the conduction during the decay phase of solar flares. With a simple quasi-steady loop model we show that the energy input in the gradual phase can be comparable to that in the impulsive phase and demonstrate how the observed cooling and confinement of the LT source can be used to constrain the wave-particle interaction.

Petrosian magnitudes were designed to help with the difficult task of determining a galaxy's total light. Although these magnitudes (taken here as the flux within 2R_P, with the inverted Petrosian index 1/eta(R_P)=0.2) can represent most of an object's flux, they do of course miss the light outside of the Petrosian aperture (2R_P). The size of this flux deficit varies monotonically with the shape of a galaxy's light-profile, i.e., its concentration. In the case of a de Vaucouleurs R^{1/4} profile, the deficit is 0.20 mag; for an R^{1/8} profile this figure rises to 0.50 mag. Here we provide a simple method for recovering total (Sersic) magnitudes from Petrosian magnitudes using only the galaxy concentration (R_90/R_50 or R_80/R_20) within the Petrosian aperture. The corrections hold to the extent that Sersic's model provides a good description of a galaxy's luminosity profile. We show how the concentration can also be used to convert Petrosian radii into effective half-light radii, enabling a robust measure of the mean effective surface brightness. Our technique is applied to the SDSS DR2 Petrosian parameters, yielding good agreement with the total magnitudes, effective radii, and mean effective surface brightnesses obtained from the NYU--VAGC Sersic R^{1/n} fits by Blanton et al. (2005). Although the corrective procedure described here is specifically applicable to the SDSS DR2 and DR3, it is generally applicable to all imaging data where any Petrosian index and concentration can be constructed.

Recent near-IR (NIR) and X-ray observations of Sagittarius A*'s spectrum have yielded several strong constraints on the transient energization mechanism, justifying a re-examination of the stochastic acceleration model proposed previously for these events. We here demonstrate that the new results are fully consistent with the acceleration of electrons via the transit-time damping process. But more importantly, these new NIR and X-ray flares now can constrain the source size, the gas density, the magnetic field, and the wave energy density in the turbulent plasma. Future simultaneous multi-wavelength observations with good spectral information will, in addition, allow us to study their temporal evolution, which will eventually lead to an accurate determination of the behavior of the plasma just minutes prior to its absorption by the black hole.

We study the acceleration in solar flares of $^3$He and $^4$He from a thermal background by parallel propagating plasma waves with a general broken power-law spectrum that takes into account the turbulence generation processes at large scales and the thermal damping effects at small scales. The exact dispersion relation for a cold plasma is used to describe the relevant wave modes. Because low-energy $\alpha$-particles only interact with small scale waves in the $^4$He-cyclotron branch, where the wave frequencies are below the $\alpha$-particle gyro-frequency, their pitch angle averaged acceleration time is at least one order of magnitude longer than that of $^3$He ions, which mostly resonate with relatively higher frequency waves in the proton-cyclotron (PC) branch. The $\alpha$-particle acceleration rate starts to approach that of $^3$He beyond a few tens of keV nucleon$^{-1}$, where $\alpha$-particles can also interact with long wavelength waves in the PC branch. However, the $^4$He acceleration rate is always smaller than that of $^3$He. Consequently, the acceleration of $^4$He is suppressed significantly at low energies, and the spectrum of the accelerated $\alpha$-particles is always softer than that of $^3$He. The model gives reasonable account of the observed low-energy $^3$He and $^4$He fluxes and spectra in the impulsive solar energetic particle events observed with the {\it Advanced Composition Explorer}. We explore the model parameter space to show how observations may be used to constrain the model.

Evidence for nonthermal activity in clusters of galaxies is well established from radio observations of synchrotron emission by relativistic electrons, and new windows (in EUV and Hard X-ray ranges) have provided more powerful tools for its investigation. The hard X-ray observations, notably from Coma, are summarized and results of a new RXTE observations of a high redshift cluster are presented. It is shown that the most likely emission mechanisms for these radiations is the inverse Compton scattering of the cosmic microwave background photons by the same electrons responsible for the radio radiation. Various scenarios for acceleration of the electrons are considered and it is shown that the most likely model is episodic acceleration by shocks or turbulence, presumably induced by merger activity, of high energy electrons injected into the intercluster medium by galaxies or active galactic nuclei.

We present data analysis and interpretation of a simple X-class flare observed with RHESSI on November 3, 2003. In contrast to other X-class flares observed previously, this flare shows a very simple morphology with well defined looptop (LT) and footpoint (FP) sources. The almost monotonic upward motion of the LT source and increase in separation of the two FP sources are consistent with magnetic reconnection models proposed for solar flares. In addition, we find that the source motions are relatively slower during the more active phases of hard X-ray emission; the emission centroid of the LT source shifts toward higher altitudes with the increase of energy; the separation between the LT emission centroids at two different photon energies is anti-correlated with the FP flux. Non-uniformity of the reconnecting magnetic fields could be a possible explanation of these features.

Stochastic acceleration of electrons and protons by waves propagating parallel to the large scale magnetic fields of magnetized plasmas is studied with emphasis on the feasibility of accelerating particles from a thermal background to relativistic energies and with the aim of determining the relative acceleration of the two species in one source. In general, the stochastic acceleration by these waves results in two distinct components in the particle distributions, a quasi-thermal and a hard nonthermal, with the nonthermal one being more prominent in hotter plasmas and/or with higher level turbulence. This can explain many of the observed features of solar flares. Regarding the proton to electron ratio, we include the effects of $^4$He in the dispersion relation and the damping of the turbulent waves by the thermal background plasma. The relative acceleration of protons and electrons is very sensitive to the ratio of the electron plasma frequency to its gyro-frequency. Protons are preferentially accelerated in weakly magnetized plasmas. The formalism developed here is applicable to the acceleration of other ion species and to other astrophysical systems.

The recent detection of variable infrared emission from Sagittarius A*, combined with its previously observed flare activity in X-rays, provides compelling evidence that at least a portion of this object's emission is produced by nonthermal electrons. We show here that acceleration of electrons by plasma wave turbulence in hot gases near the black hole's event horizon can account both for Sagittarius A*'s mm and shorter wavelengths emission in the quiescent state, and for the infrared and X-ray flares, induced either via an enhancement of the mass accretion rate onto the black hole or by a reorganization of the magnetic field coupled to the accretion gas. The acceleration model proposed here produces distinct flare spectra that may be compared with future coordinated multi-wavelength observations. We further suggest that the diffusion of high energy electrons away from the acceleration site toward larger radii might be able to account for the observed characteristics of Sagittarius A*'s emission at cm and longer wavelengths.

Stochastic acceleration of $^3$He and $^4$He from a thermal background by parallel propagating turbulent plasma waves with a single power-law spectrum of the wavenumber is studied. In the model, both ions interact with several resonant modes. When one of these modes dominates, the acceleration rate is reduced considerably. At low energies, this happens for $^4$He, but not for $^3$He where contributions from the two stronger modes are comparable so that acceleration of $^3$He is very efficient. As a result, the acceleration of $^4$He is suppressed by a barrier below $\sim 100$ keV nucleon$^{-1}$ and there is a prominent quasi-thermal component in the $^4$He spectra, while almost all the injected $^3$He ions are accelerated to high energies. This accounts for the large enrichment of $^3$He at high energies observed in impulsive solar energetic particle events. With reasonable plasma parameters this also provides a good fit to the spectra of both ions. Beyond $\sim 1$ MeV nucleon$^{-1}$, the spectrum of $^3$He is softer than that of $^4$He, which is consistent with the observed decrease of the $^3$He to $^4$He ratio with energy. This study also indicates that the acceleration, Coulomb losses and diffusive escape of the particles from the acceleration site {\it all} play important roles in shaping the ion spectra. This can explain the varied spectral shapes observed recently by the {\it Advanced Composition Explorer}.

We present spectra of the optical transient of GRB021004 obtained with the Hobby-Eberly telescope starting 15.48, 20.31 hours, and 4.84 days after the burst and a spectrum obtained with the H. J. Smith 2.7 m Telescope starting 14.31 hours after the burst. GRB021004 is the first afterglow whose spectrum is dominated by absorption lines from high ionization species with multiple velocity components separated by up to 3000 km/s. We argue that these lines are likely to come from shells around a massive progenitor star. The high velocities and high ionizations arise from a combination of acceleration and flash-ionization by the burst photons and the wind velocity and steady ionization by the progenitor. We also analyze the broad-band spectrum and the light curve. We distinguish six components along the line of sight: (1) The z~2.293 absorption lines arise from the wind of a massive star. For a mass loss rate of ~6 x 10^{-5} solar masses per year, this component also provides the external medium to create the afterglow light. (2) A second shell produces absorption lines with a relative velocity of 560 km/s, and this is associated with the shell created by the fast massive star wind blowing a bubble in the preceding slow wind at a radial distance of order 10 pc. (3) More distant clouds within the host galaxy lie between 30-2500 pc, where they have been ionized by the burst. (4-6) The massive star wind has clumps with radii and over-densities of 0.022, 0.063, and 0.12 parsecs and 50%, 10%, and 10% respectively. The immediate progenitor of the burster could either be a WC-type Wolf-Rayet star or a highly evolved star whose original mass was just too small for it to become a WN-type Wolf-Rayet star.

Understanding of the nature of gamma-ray bursts (GRBs) is one of the challenging problem facing the astrophysics community. The field of gamma ray bursts is a rapidly developing one and we expect that new missions, like SWIFT and GLAST will bring the field to a new quantitative level. The detail of the explosion, the formation of the fireball, its propagation, the generation of the shocks, the source of the magnetic fields and soft photons, the particle acceleration process, and the details of the radiation process are all still outstanding questions. The goal of this review is to attract the attention of the community to a number of physical processes (MHD turbulence and particle acceleration, magnetic reconnection, etc.) that seem to be very relevant for these sources.

Recent observations of excess radiation at extreme ultraviolet and hard X-ray energies straddling the well known thermal soft X-ray emission have provided new tools and puzzles for investigation of the acceleration of nonthermal particles in the intercluster medium of clusters of galaxies. It is shown that these radiations can be produced by the inverse Compton upscattering of the cosmic microwave background photons by the same population of relativistic electrons that produce the well known diffuse radio radiation via the synchrotron mechanism. It is shown that the commonly discussed discrepancy between the value of the magnetic field required for the production of these radiation with that obtained from Faraday rotation measures could be resolved by more realistic models and by considerations of observational selection effects. In a brief discussion of the acceleration process it is argued that the most likely scenario is reacceleration of injected relativistic electrons involving shocks and turbulence. The seed electrons cannot be the background thermal electrons because of energetic considerations, and a steady state situation may not agree with the details of the observed spectra. Episodic generation of shocks and turbulence or episodic injection of relativistic electrons is a more likely scenario for acceleration.

The discovery of hard X-rays from tops of flaring loops by the HXT of YOHKOH represents a significant progress in the understanding of solar flares. This report describes the properties of 20 limb flares observed by YOHKOH from October 1991 to August 1998, 15 of which show detectable impulsive looptop emission. Considering the finite dynamic range (about a decade) of the detection it can be concluded that looptop emission is a common feature of all flares. The light curves and images of a representative flare are presented and the statistical properties of the footpoint and looptop fluxes and spectral indexes are summarized. The importance of these observations, and those expected from RHESSI with its superior angular, spectral and temporal resolution, in constraining the acceleration models and parameters is discussed. briefly.

Several patterns have been discovered in how some spectral characteristics change during the decaying phase of long (> few s.) GRB pulses. We compare these observed signatures with those expected from a relativistically expanding shell. Within the internal shock model and assuming a short cooling time, we show that the angular dependence in arrival time can explain the general characteristics of long GRB pulses. This includes the pulse shape, with a fast rise and a slower decay, ~ (1+t/tau)^2, and the spectral evolution, which can be described by the hardness-intensity correlation (HIC), with the intensity being proportional to the square of the hardness. A variation of the relevant time scales involved (the angular spreading and the dynamic) can explain the broad, observed dispersion of the HIC index. Reasonable estimates of physical parameters lead to situations where the HIC relation deviates from a pure power law; features that are indeed present in the observations. Depending on the relative values of the rise and decay times of the intrinsic light curve, the spectral/temporal behavior, as seen by an observer, will produce hard-to-soft or tracking pulses. The observed spectrum is a superposition of many intrinsic spectra arriving from different parts of the shell with varying spectral shifts. Therefore, it will be broader than the emitted spectrum and its spectral parameters could have complex relations with the intrinsic ones. Furthermore, we show that the softening of the low-energy power-law index, that has been observed in some pulses, can be explained by geometric effects and does not need to be an intrinsic behavior.

The aim of this paper is to describe new statistical methods for determination of the correlations among and distributions of physical parameters from a multivariate data with general and arbitrary truncations and selection biases. These methods, developed in collaboration with B. Efron of Department of Statistics at Stanford, can be used for analysis of combined data from many surveys with different and varied observational selection criteria. For clarity we will use the luminosity function of AGNs and its evolution to demonstrate the methods. We will first describe the general features of data truncation and present a brief review of past methods of analysis. Then we will describe the new methods and results from simulations testing their accuracy. Finally we will present the results from application of the methods to a sample of quasars.

The discovery of hard X-ray sources near the top of a flaring loop by the HXT instrument on board the YOHKOH satellite represents a significant progress towards the understanding of the basic processes driving solar flares. In this paper we extend the previous study of limb flares by Masuda (1994) by including all YOHKOH observations up through August 1998. We report that from October 1991 to August 1998, YOHKOH observed 20 X-ray bright limb flares (where we use the same selection criteria as Masuda), of which we have sufficient data to analyze 18 events, including 8 previously unanalyzed flares. Of these 18 events, 15 show detectable impulsive looptop emission. Considering that the finite dynamic range (about a decade) of the detection introduces a strong bias against observing comparatively weak looptop sources, we conclude that looptop emission is a common feature of all flares. We summarize the observations of the footpoint to looptop flux ratio and the spectral indices. We present light curves and images of all the important newly analyzed limb flares. Whenever possible we present results for individual pulses in multipeak flares and for different loops for multiloop flares. We then discuss the statistics of the fluxes and spectral indices of the looptop and footpoint sources taking into account observational selection biases. The importance of these observations (and those expected from the scheduled HESSI satellite with its superior angular spectral and temporal resolution) in constraining acceleration models and parameters is discussed briefly.

In this paper, we explore time resolved Gamma-Ray Burst (GRB) spectra in the context of the synchrotron emission model presented in Lloyd and Petrosian (2000; LP00). First, we show that our model - which involves three distinct emission regimes - can provide excellent fits to the time resolved spectra of GRBs, and we present these results for a few bursts. We then describe how the phenomenological Band spectrum (Band et al., 1993) can be interpreted in the context of our models based on the value of the low energy photon index $\alpha$. We discuss the types of correlations one would expect to observe among the Band parameters if these models are correct. We then compare these predictions to the existing data, combining a sample of 2,026 time resolved spectra (from approximately 80 bursts). We show that the correlations found in the data are consistent with the models, and discuss the constraints they place on the emission physics. In particular, we find a 4-sigma negative correlation between the peak of the nu-Fnu spectrum, Ep, and the low energy photon index $\alpha$ for bursts with $-2/3 < \alpha < 0$, in contrast to what is predicted by the instrumental effect discussed in LP00. We suggest that this correlation is simply due to the mechanism responsible for producing $\alpha$'s above the value of -2/3 - namely, a decreasing mean pitch angle of the electrons. We also show that Ep is correlated with the photon flux, and interpret this as a result of changing magnetic field or characteristic electron energy between emission episodes. Finally, we discuss the implications our results have on particle acceleration in GRBs, and prospects for further testing these models with the anticipated data from HETE-2, Swift and GLAST.

Some clusters of galaxies in addition to thermal bremsstrahlung (TB), emit diffuse radiation from the intercluster medium (ICM) at radio, EUV and hard x-ray (HXR) ranges. The radio radiation is due to synchrotron by relativistic electrons, and the inverse Compton (IC) scattering by the cosmic microwave background radiation of the same electrons is the most natural source for the HXR and perhaps the EUV emissions. However, simple estimates give a weaker magnetic field than that suggested by Faraday rotation measurements. Consequently, non-thermal bremsstrahlung (NTB) and TB have also been suggested as sources of these emissions. We show that NTB cannot be the source of the HXRs and that the difficulty with the low magnetic field in the IC model is alleviated if we take into account the effects of observational bias, nonisotropic pitch angle distribution and spectral breaks. We derive a spectrum for the radiating electrons and discuss acceleration scenarios. We show that continuous and in situ acceleration in the ICM of the background thermal electrons requires unreasonably high energy input and acceleration of injected relativistic electrons gives rise to a much flatter spectrum than desired, unless a large fraction of electrons escape the ICM, in which case one obtains EUV and HXR emissions extending well beyond the boundaries of the cluster. A continuous emission by a cooling spectrum resulting from interaction with ICM of electrons accelerated elsewhere also suffers from similar shortcomings. The most likely scenario appears to be an episodic injection-acceleration model, whereby one obtains a time dependent spectrum that for certain phases of its evolution satisfies all the requirements.

We investigate synchrotron emission models as the source of gamma-ray burst spectra. We show that including the possibility for synchrotron self-absorption, a ``smooth cutoff'' to the electron energy distribution, and an anisotropic distribution for the electron pitch angles produces a whole range of low energy spectral behavior. In addition, we show that the procedure of spectral fitting to GRB data over a finite bandwidth can introduce a spurious correlation between spectral parameters - in particular, the value of the peak of the nu F_nu spectrum, E_p, and the low energy photon spectral index alpha (the lower E_p is, the lower (softer) the fitted value of alpha will be). From this correlation and knowledge of the E_p distribution, we show how to derive the expected distribution of alpha. We show that optically thin synchrotron models with an isotropic electron pitch angle distribution can explain the distribution of alpha below alpha=-2/3. This agreement is achieved if we relax the unrealistic assumption of the presence of a sharp low energy cutoff in the spectrum of accelerated electrons, and allow for a more gradual break. We show that this low energy portion of the electron spectrum can be at most flat. We also show that optically thin synchrotron models with an anisotropic electron pitch angle distribution can explain all bursts with -2/3 < alpha <= 0$. The very few bursts with low energy spectral indices that fall above alpha=0 may be due the presence of a the synchrotron self-absorption frequency entering the lower end of the BATSE window. Our results also predict a particular relationship between alpha and E_p during the temporal evolution of a GRB. We give examples of spectral evolution in GRBs and discuss how the behavior are consistent with the above models.

The time profiles of many gamma-ray bursts consist of distinct pulses, which offers the possibility of characterizing the temporal structure of these bursts using a relatively small set of pulse shape parameters. We have used a pulse decomposition procedure to analyze the Time-to-Spill (TTS) data for all bursts observed by BATSE up through trigger number 2000, in all energy channels for which TTS data is available. We obtain amplitude, rise and decay timescales, a pulse shape parameter, and the fluences of individual pulses in all of the bursts. We investigate the correlations between brightness measures (amplitude and fluence) and timescale measures (pulse width and separation) which may result from cosmological time dilation of bursts, or from intrinsic properties of burst sources or from selection effects. The effects of selection biases are evaluated through simulations. The correlations between these parameters among pulses within individual bursts give a measure of the intrinsic effects while the correlations among bursts could result both from intrinsic and cosmological effects. We find that timescales tend to be shorter in bursts with higher peak fluxes, as expected from cosmological time dilation effects, but also find that there are non-cosmological effects contributing to this inverse correlation. We find that timescales tend to be longer in bursts with higher total fluences, contrary to what is expected from cosmological effects. We also find that peak fluxes and total fluences of bursts are uncorrelated, indicating that they cannot both be good distance indicators for bursts.

The time profiles of many gamma-ray bursts consist of distinct pulses, which offers the possibility of characterizing the temporal structure of these bursts using a relatively small set of pulse shape parameters. This pulse decomposition analysis has previously been performed on a small sample of bright long bursts using binned data from BATSE, which comes in several data types, and on a sample of short bursts using the BATSE Time-Tagged Event (TTE) data type. We have developed an interactive pulse-fitting program using the phenomenological pulse model of Norris, et al. and a maximum-likelihood fitting routine. We have used this program to analyze the Time-to-Spill (TTS) data for all bursts observed by BATSE up through trigger number 2000, in all energy channels for which TTS data is available. We present statistical information on the attributes of pulses comprising these bursts, including relations between pulse characteristics in different energy channels and the evolution of pulse characteristics through the course of a burst. We carry out simulations to determine the biases that our procedures may introduce. We find that pulses tend to have shorter rise times than decay times, and tend to be narrower and peak earlier at higher energies. We also find that pulse brightness, pulse width, and pulse hardness ratios do not evolve monotonically within bursts, but that the ratios of pulse rise times to decay times tend to decrease with time within bursts.

We investigate synchrotron emission models as the source of gamma ray burst spectra. We show that allowing for synchrotron self absorption and a ``smooth cutoff'' to the electron energy distribution produces a wide range of low energy spectral behavior. We show that there exists a correlation between the value of the peak of the $\nu F_{\nu}$ spectrum, $E_{p}$, and the low energy spectral index $\alpha$ as determined by spectral fits over a finite bandwidth. Finally, we discuss the implications of synchrotron emission from internal shocks for GRB spectral evolution.

We test the models of synchrotron emission presented in Part I of this series (Lloyd & Petrosian, these proceedings) against the distributions and evolution of GRB spectral parameters (particularly the low energy index, $\alpha$). With knowledge of the $E_{p}$ distribution and the correlation between $\alpha$ and $E_{p}$ presented in Part I, we show how to derive the expected distribution of $\alpha$ from fits to optically thin synchrotron spectra, and compare this with the observed distribution. We show that there is no difficulty explaining bursts below the ``line of death'', $\alpha < -2/3$, and that these bursts indicate that the spectrum of accelerated electrons must flatten or decline at low energies. Bursts with low energy spectral indices that fall above this limit are explained by the synchrotron self-absorption frequency entering the lower end of the BATSE window. Finally, we discuss a variety of spectral evolution behavior seen in GRBs and explain this behavior in the context of synchrotron emission from internal shocks.

Observational aspects of solar flares relevant to the acceleration process of electrons and protons are reviewed and it is shown that most of these observations can be explained by the interaction with flare plasma of a power law energy distribution of electrons (and protons) that are injected at the top of a flaring loop, in the so-called thick target model. Some new observations that do not agree with this model are described and it is shown that these can be explained most naturally if most of the energy released by the reconnection process goes first into the generation of plasma turbulence, which accelerates, scatters and traps the ambient electrons near the top of the loop stochastically. The resultant bremsstrahlung photon spectral and spatial distributions agree with the new observations. This model is also justified by some theoretical arguments. Results from numerical evaluation of the spectra of the accelerated electrons and their bremsstrahlung emission are compared with observations and shown how one can constrain the model parameters describing the flare plasma and the spectrum and the energy density of the turbulence.

We present results of correlation studies, examining the association between the peak of the nu F_nu spectrum of gamma ray bursts, E_p, with the burst's energy fluence and photon peak flux. We discuss methods to account for data truncation in E_p and fluence or flux when performing the correlation analyses. However, because bursts near the detector threshold are not usually able to provide reliable spectral parameters, we focus on results for the brightest bursts in which we can better understand the selection effects relevant to E_p and burst strength. We find that there is a strong correlation between total fluence and E_p. We discuss these results in terms of both cosmological and intrinsic effects. In particular, we show that for realistic distributions of the burst parameters, cosmological expansion alone cannot account for the correlation between E_p and total fluence; the observed correlation is likely a result of an intrinsic relation between the burst rest-frame peak energy and the total radiated energy. We investigate this latter scenario in the context of synchrotron radiation from external and internal shock models of GRBs. We find that the internal shock model is consistent with our interpretation of the correlation, while the external shock model cannot easily explain this intrinsic relation between peak energy and burst radiated energy.

The aim of this paper is to investigate the spatial structure of the impulsive phase hard X-ray emission from solar flares. This work is motivated by the YOHKOH and the forthcoming HESSI observations. Summarizing past results, it is shown that the transport effects can account for the observations by inhomogeneous loops where there is a strong field convergence and/or density enhancement at the top of the flaring loop. Scattering by plasma turbulence at the acceleration site or pancake type pitch angle distribution of the accelerated electrons can also give rise to enhanced emission at the loop tops. These could be a natural consequence of acceleration by plasma waves. This paper considers a general case of stochastic scattering and acceleration that leads to an isotropic pitch angle distribution and an enhanced emission from the loop tops or the acceleration site. Following the formalism developed in earlier papers the strength and the spectrum of the radiation expected from the acceleration site and the foot points are evaluated and their dependence on the parameters describing the acceleration process and the flare plasma are determined. The theoretical ratio of these two intensities and relative values of their spectral indices are compared with the YOHKOH observations, demonstrating that the above mentioned parameters can be constrained with such observations. It is shown that future high spatial and spectral resolution observations, for example those expected from HESSI, can begin to distinguish between different models and constrain their parameters.

There have been several suggestions of the existence of cosmological redshift signatures in the temporal and spectral characteristics of gamma-ray bursts. However, recent discoveries of afterglows and redshift measurements indicate the presence of broad ``luminosity functions'' which may overwhelm such weaker cosmological signatures. The primary goal of this paper is to determine if the intrinsic and cosmological dispersions can be separated. We have expanded the search for cosmological signatures to several other temporal and spectral features, have determined correlations which could arise from the cosmological redshift of the sources and have carried out tests to determine if the observed correlations can be due to cosmology. We find that the intrinsic dispersions are the dominant factors and that detection of cosmological signatures must await accumulation of a much larger number of identification with galaxies and measurements of redshifts.

Truncated data plays an important role in the statistical analysis of astronomical observations as well as in survival analysis. The motivating example for this paper concerns a set of measurements on quasars in which there is double truncation. That is, the quasars are only observed if their luminosity occurs within a certain finite interval, bounded at both ends, with the interval varying for different observations. Nonparametric methods for the testing and estimation of doubly truncated data are developed. These methods extend some known techniques for data that is only truncated on one side, in particular Lynden-Bell's estimator and the truncated version of Kendall's tau statistic. However the kind of hazard function arguments that underlie the one-sided methods fail for two-sided truncation. Bootstrap and Markov Chain Monte Carlo techniques are used here in their place. Finally, we apply these techniques to the quasar data, answering a question about their long-term luminosity evolution.

We investigate the cosmological evolution of GRBs, using the total gamma ray fluence as a measure of the burst strength. This involves an understanding of the distributions of the spectral parameters of GRBs as well as the total fluence distribution - both of which are subject to detector selection effects. We present new non-parametric statistical techniques to account for these effects, and use these methods to estimate the true distribution of the peak of the nu F_nu spectrum, E_p, from the raw distribution. The distributions are obtained from four channel data and therefore are rough estimates. Here, we emphasize the methods and present qualitative results. Given its spectral parameters, we then calculate the total fluence for each burst, and compute its cumulative and differential distributions. We use these distributions to estimate the cosmological rate evolution of GRBs, for three cosmological models. Our two main conclusions are the following: 1) Given our estimates of the spectral parameters, we find that there may exist a significant population of high E_p bursts that are not detected by BATSE, 2) We find a GRB co-moving rate density quite different from that of other extragalactic objects; in particular, it is different from the recently determined star formation rate.

We describe the classical cosmological tests, such as the Log$N$-Log$S$, redshift-magnitude and angular diameter tests, and propose some new tests of the evolution of galaxies and the universe. Most analyses of these tests treat the problem in terms of a luminosity function and its evolution which can lead to incorrect conclusions when dealing with high redshift sources. We develop a proper treatment in three parts. In the first part we describe these tests based on the isophotal values of the quantities such as flux, size or surface brightness. We show the shortcomings of the simple point source approximation based solely on the luminosity function and consideration of the flux limit. We emphasize the multivariate nature of the problem and quantify the effects of other selection biases due to the surface brightness and angular size limitations. In these considerations the surface brightness profile plays a critical role. In the second part we show that considerable simplification over the complicated isophotal scheme is achieved if these test are carried out in some sort of metric scheme, for example that suggested by Petrosian (1976). This scheme, however, is limited to well resolved sources. Finally, we describe the new tests, which use the data to a fuller extent than the isophotal or metric based tests, and amount to simply counting the pixels or adding their intensities as a function of the pixel surface brightness, instead of dealing with surface brightness, sizes and fluxes of individual galaxies. We show that the data analysis and its comparison with the theoretical models of the distributions and evolution of galaxies has the simplicity of the metric test and utilizes the data more fully than the isophotal test.

We use several quasar samples to determine the density and luminosity evolution of quasars. Combining different samples and accounting for varying selection criteria require tests of correlation and the determination of distributions for multiply truncated data. We describe new non-parametric techniques for accomplishing these tasks, which have been developed recently by Efron and Petrosian (1998). We use matter dominated cosmological models with either zero cosmological constant or zero spatial curvature. Of the two most common models of luminosity evolution, L = exp(k t(z)) and L = (1+z)^k', we find the second model to be a better description of the data at all luminosities; we find k' = 2.58 ([2.14,2.91] one sigma region) for the Einstein - de Sitter model. Using this form of luminosity evolution we determine a global luminosity function and the evolution of the co-moving density for the two types of cosmological models. For the Einstein - de Sitter cosmological model we find a relatively strong increase in co-moving density up to a redshift of about 2, at which point the density peaks and begins to decrease rapidly. This is in agreement with results from high redshift surveys. We find some co-moving density evolution for all cosmological models, with the rate of evolution lower for models with lower matter density. We find that the local cumulative luminosity function exhibits the usual double power law behavior. The luminosity density (i.e. the total rate of energy output of quasars) is found to increase rapidly at low redshift and to peak at around z = 2. Our results for the luminosity density are compared to results from high redshift surveys and to the variation of the star formation rate with redshift.

The distributions of spectral characteristics and their correlations with fluence, peak flux, or duration, are essential in understanding the nature of GRBs. However, the selection effects involved in detecting GRBs can distort these distributions. Here, we discuss how to deal with selection effects involving the peak energy E_p of the GRB nu F_nu spectrum, which suffers from both an upper and lower threshold. We describe a new method to account for this double-sided truncation, and show that the true distribution of E_p is significantly different from the observed distribution.

We argue that recent observations of afterglows and some theoretical considerations indicate that the gamma-ray fluence may be a better measure of strength of a burst than its peak flux. We discuss how the distribution of the fluence, or any physical quantity related to the peak counts (or flux), can be obtained with proper correction for threshhold effects. Using this method, we compare the cumulative and differential distributions of flux and fluence. We find some differences between these distributions, but more remarkable are the similarities between the distributions. The other striking feature is how different these distributions are from those of other extragalactic objects. Using the fluence distributions, we derive the expected comoving density rate evolution for different assumed total luminosity; we compare this with the GRB rate expected from the recently determined star formation rate.

This paper extends our earlier work on the acceleration of low-energy electrons by plasma turbulence to include the effects of finite temperature of the plasma. We consider the resonant interaction of thermal electrons with the whole transverse branch of plasma waves propagating along the magnetic field. We show that our earlier published results for acceleration of low-energy electrons can be applied to the case of finite temperature if a sufficient level of turbulence is present. From comparison of the acceleration rate of the thermal particles with the decay rate of the waves with which they interact, we determine the required energy density of the waves as a fraction of the magnetic energy density, so that a substantial fraction of the background plasma electrons can be accelerated. The dependence of this value on the plasma parameter alpha = omega_pe / Omega_e (the ratio of electron plasma frequency to electron gyrofrequency), plasma temperature, and turbulence spectral parameters is quantified. We show that the result is most sensitive to the plasma parameter alpha. We estimate the required level of total turbulence by calculating the level of turbulence required for the initial acceleration of thermal electrons and that required for further acceleration to higher energies.

We investigate the possibility of stochastic acceleration of background low-energy electrons by turbulent plasma waves. We consider the resonant interaction of the charged particles with all branches of the transverse plasma waves propagating parallel to a uniform magnetic field. Numerical results and asymptotic analytic solutions valid at non-relativistic and ultra-relativistic energies are obtained for the acceleration and scattering times of electrons. These times have a strong dependence on plasma parameter alpha = Omega_pe / Omega_e (the ratio of electron plasma frequency to electron gyrofrequency) and on the spectral index of plasma turbulence. It is shown that particles with energies above certain critical value may interact with higher frequency electromagnetic plasma waves and this interaction is allowed only in plasmas with alpha < 1. We show that for non-relativistic and semi-relativistic electrons in low-alpha plasmas the ratio of the acceleration time to the scattering time can be less than unity for a wide range of energies. From this we conclude that the transport equation derived for cosmic rays which requires this ratio to be much larger than one is not applicable at these energies. An approximate "critical" value of particle energy above which the dynamics of charged particles may be described by this transport equation is determined as a function of plasma parameters. We propose new transport equation for the opposite limit (energies less than this critical value) when the acceleration rate is much faster than the pitch angle scattering rate. This equation is needed to describe the electron dynamics in plasmas with alpha <= 0.1.

We perform several nonparametric correlation tests on the BATSE 3B data to search for evidence of cosmological time dilation. These tests account for the effects of data truncation due to threshold effects in both limiting brightness and limiting duration, enabling us to utilize a larger number of bursts than in previous analyses. We find little significant evidence for correlation between various measures of peak intensity and duration, but the tests cannot conclusively rule out time dilation factors of 2 or less without more data. There is stronger evidence for a positive correlation between fluence and duration, which if confirmed would rule out simple no-evolution cosmological models unless there is a strong intrinsic correlation between the total radiant energy and the duration of bursts.

In the use and interpretation of $\log{N}$--$\log{S}$ distributions for gamma-ray bursts, burst peak flux has typically been used for $S$. We consider here the use of the fluence as a measure of $S$, which may be a more appropriate quantity than the peak flux in such highly variable sources. We demonstrate how using the BATSE trigger data we can determine the selection effects on fluence. Then using techniques developed elsewhere to account for the important threshold effects and correlations. Applying the appropriate corrections to the distributions, we obtain a fluence distribution which shows a somewhat sharper break than the peak flux distribution, implying a possibly narrower fluence luminosity distribution. If bursts are at cosmological distances, these observations together indicate that evolution of the luminosity function is required.

Many of the important conclusions about Gamma-Ray Bursts follow from the distributions of various quantities such as peak flux or duration. We show that for astrophysical transients such as bursts, multiple selection thresholds can lead to various forms of data truncation, which can strongly affect the distributions obtained from the data if not accounted for properly. Thus the data should be considered to form a multivariate distribution. We also caution that if the variables forming the multivariate distribution are not statistically independent of each other, further biases can result. A general method is described to properly account for these effects, and as a specific example we extract the distributions of flux and duration from the BATSE 3B Gamma-Ray Burst data. It is shown that properly accounting for the aforementioned biases tends to increase the slope of the $\log{N}$-$\log{S}$ relation at low values of $S$, and dramatically increases the number of short duration bursts.