We report the detection of high-energy gamma-ray signal towards the young star-forming region, W40. Using 10-year Pass 8 data from the Fermi Large Area Telescope (Fermi-LAT), we extracted an extended gamma-ray excess region with a significance of about 18sigma. The radiation has a spectrum with a photon index of 2.49 +/- 0.01. The spatial correlation with the ionized gas content favors the hadronic origin of the gamma-ray emission. The total cosmic-ray (CR) proton energy in the gamma-ray production region is estimated to be the order of 10^47 erg. However, this could be a small fraction of the total energy released in cosmic rays (CRs) by local accelerators, presumably by massive stars, over the lifetime of the system. If so, W40, together with earlier detections of gamma-rays from Cygnus cocoon, Westerlund 1, Westerlund 2, NGC 3603, and 30 Dor C, supports the hypothesis that young star clusters are effective CR factories. The unique aspect of this result is that the gamma-ray emission is detected, for the first time, from a stellar cluster itself, rather than from the surrounding "cocoons".

Starburst galaxies are huge reservoirs of cosmic rays (CRs) and these CRs convert a significant fraction of their energy into gamma-rays by colliding with the interstellar medium (ISM). The produced GeV gamma-ray emission is temporally stable and the flux is found to correlate well with indicators of star formation rates, such as the total infrared (IR) luminosity $L_{\rm IR}$ and monochromatic radio continuum luminosity at 1.4 GHz $L_{\rm 1.4\ GHz}$, i.e., following $L_{\gamma}-L_{\rm IR}$ and $L_{\gamma}-L_{\rm 1.4\ GHz}$ relations. Recently, gamma-ray excesses are reported to be spatially coincident with two starburst galaxies NGC 3424 and UGC 11041 in the fourth \textsl{Fermi} Large Area Telescope (LAT) source catalog (4FGL). Different from other starburst galaxies detected by \textsl{Fermi}-LAT, we find that the gamma-ray emission associated with NGC 3424 and UGC 11041 show significant flux variability. With relatively weak infrared and radio emission, NGC 3424 and UGC 11041 appear as outliers of the $L_{\gamma}-L_{\rm IR}$ and $L_{\gamma}-L_{\rm 1.4\ GHz}$ relations of starburst galaxies. These results suggest that NGC 3424 and UGC 11041 may harbor obscured active galactic nuclei (AGNs) and the AGN activities provide the dominant contribution to the gamma-ray emission as compared to that provided by the starburst activities.

Neutrino annihilation process on a hyperaccreting disk is one of the leading models to explain the generation of relativistic jets of gamma-ray bursts (GRBs). However, there is still uncertainty regarding the NAE of neutron star-accretion disc (NS-disc) system because of complicated microphysics processes and effects of strong magnetic field. In this paper, we investigate the NAE by assuming that the prompt jet of GRB 070110 is driven by neutrino pair annihilation in the NS-disc system. Our calculation shows $\eta_{\nu\bar{\nu}}> 1.2 \times 10^{-3}$ under the estimated accretion rate $\dot{M}\simeq 0.04\rm M_{\odot}\cdot s^{-1}$. Independent of the detailed accretion disc models, our result shows that the magnetic field may play an important role in the neutrino annihilation process on the hyperaccreting magnetized accretion disc. Compared with the theoretical value of $\eta_{\rm \nu \bar{\nu}}$ of the nonmagnetized BH-disc system, the NAE should increase significantly in the case of the NS-disc system if the GRB is powered by magnetar-disc system.

The measurement of the speed of gravitational waves (GWs) is useful to distinguish general relativity from massive gravity. We propose a new model-independent strategy to measure the speed of GWs with the distorted pulsars. Theoretically, when the the GW frequencies from a distorted pulsar are twice the frequencies of EM pulses, they should be emitted at the same time. By measuring the arrival times of these two signals emitted at the same time, the speed of GWs can be calculated with the time difference. Specifically, when the glitches of pulsars are taken into consideration, some pulsars in our Galaxy and nearby galaxies are potential to test our new strategy at high accuracy in the foreseeable future. On the other hand, the new method is meaningful as a motivation for future design of instruments.

It has been found that there is a quasi-linear scaling relationship between the gamma-ray luminosity in GeV energies and the total infrared luminosity of star-forming galaxies, i.e. $L_{\gamma}\propto L_{\rm IR}^{\alpha}$ with $\alpha\simeq 1$. However, the origin of this linear slope is not well understood. Although extreme starburst galaxies can be regarded as calorimeters for hadronic cosmic ray interaction and thus a quasi-linear scaling may hold, it may not be the case for low star-formation-rate (SFR) galaxies, as the majority of cosmic rays in these galaxies are expected to escape. We calculate the gamma-ray production efficiency in star-forming galaxies by considering realistic galaxy properties, such as the gas density and galactic wind velocity in star-forming galaxies. We find that the slope for the relation between gamma-ray luminosity and the infrared luminosity gets steeper for low infrared luminosity galaxies, i.e. $\alpha\rightarrow 1.6$, due to increasingly lower efficiency for the production of gamma-ray emission. We further find that the measured data of the gamma-ray luminosity is compatible with such a steepening. The steepening in the slope suggests that cosmic-ray escape is very important in low-SFR galaxies.

We investigate emission signatures of binary compact star gravitational wave sources consisting of strongly magnetized neutron stars (NSs) and/or white dwarfs (WDs) in their late-time inspiral phase. Because of electromagnetic interactions between the magnetospheres of the two compact stars, a substantial amount of energy will be extracted, and the resultant power is expected to be $\sim 10^{38} - 10^{44}$ erg/s in the last few seconds before the two stars merge, when the binary system contains a NS with a surface magnetic field $10^{12}$ G. The induced electric field in the process can accelerate charged particles up to the EeV energy range. Synchrotron radiation is emitted from energetic electrons, with radiative energies reaching the GeV energy for binary NSs and the MeV energy for NS - WD or double WD binaries. In addition, a blackbody component is also presented and it peaks at several to hundreds keV for binary NSs and at several keV for NS - WD or double WD binaries. The strong angular dependence of the synchrotron radiation and the isotropic nature of the blackbody radiation lead to distinguishable modulation patterns between the two emission components. If coherent curvature radiation is presented, fast radio bursts could be produced. These components provide unique simultaneous electromagnetic signatures as precursors of gravitational wave events associated with magnetized compact star mergers and short gamma ray bursts (e.g., GRB 100717).

Star-forming regions on different scales, such as giant molecular clouds in our Galaxy and star-forming galaxies, emit GeV gamma-rays. These are thought to originate from hadronic interactions of cosmic-ray (CR) nuclei with the interstellar medium. It has recently been shown that the gamma-ray luminosity ($L_\gamma$) of star-forming galaxies is well correlated with their star formation rates (SFR). We investigated \textsl{Fermi} data of eight Galactic molecular clouds in the Gould belt and found that molecular clouds do not follow the $L_\gamma-{\rm SFR}$ correlation of star-forming galaxies. We also compared the scaling relations of gamma-ray luminosity, SFR, and the gas mass for molecular clouds and star-forming galaxies. Using a multiple-variable regression analysis, we found different dependences of gamma-ray emission on SFR or mass for molecular clouds and star-forming galaxies. This suggests that different mechanisms may govern the production of gamma-rays in these two types of sources. Specifically, the strong dependence on mass supports that gamma-ray emission of molecular clouds primarily comes from {\em \textup{passive}} interaction by diffuse Galactic CRs, whereas the strong dependence on SFR supports that gamma-ray emission of star-forming galaxies originates from CRs that are accelerated by local {\em \textup{active}} sources.

Many galaxy clusters have giant halos of non-thermal radio emission, indicating the presence of relativistic electrons in the clusters. Relativistic protons may also be accelerated by merger and/or accretion shocks in galaxy clusters. These cosmic-ray (CR) electrons and/or protons are expected to produce gamma-rays through inverse-Compton scatterings or inelastic $pp$ collisions respectively. Despite of intense efforts in searching for high-energy gamma-ray emission from galaxy clusters, conclusive evidence is still missing so far. Here we report the discovery of $\ge 200$ MeV gamma-ray emission from the Coma cluster direction with an unbinned likelihood analysis of the 9 years of {\it Fermi}-LAT Pass 8 data. The gamma-ray emission shows a spatial morphology roughly coincident with the giant radio halo, with an apparent excess at the southwest of the cluster. Using the test statistic analysis, we further find tentative evidence that the gamma-ray emission at the Coma center is spatially extended. The extended component has an integral energy flux of $\sim 2\times 10^{-12}{\rm \ erg\ cm^{-2}\ s^{-1}}$ in the energy range of 0.2 - 300 GeV and the spectrum is soft with a photon index of $\simeq-2.7$. Interpreting the gamma-ray emission as arising from CR proton interaction, we find that the volume-averaged value of the CR to thermal pressure ratio in the Coma cluster is about $\sim 2\%$. Our results show that galaxy clusters are likely a new type of GeV gamma-ray sources, and they are probably also giant reservoirs of CR protons.

We study the gravitational wave (GW) production induced by the asymmetric jets of gamma-ray bursts (GRBs). The asymmetric jets result in a recoil force acted on the central compact object, whose motion leads to emission of GW. Under reasonable assumptions and simplifications, we derive the analytic form of the produce GWs. The amplitude of emitted GWs is estimated to be relatively low, but possibility exists that they can be detected by future experiments such as the Einstein Telescope. We find the dynamical properties of the central object, which is difficult to be studied via the electromagnetic (EW) channel, can be inferred by measuring the emitted GWs. Moreover, we find the emitted GWs can be used determine whether the relativistic jets is launched by the neutrino annihilation process or the Blandford-Znajek process, which cannot be clearly distinguished by the current GRB observations. Our work manifests the importance of the GW channel in multi-messenger astronomy. The physical information encoded in the GW and EW emissions of an astrophysical object is complementary to each other; in case some physics can not be effectively investigated using the EW channel alone, including the GW channel can be very helpful.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

It has been suggested that high-energy gamma-ray emission ($>100{\rm MeV}$) of nearby star-forming galaxies may be produced predominantly by cosmic rays colliding with the interstellar medium through neutral pion decay. The pion-decay mechanism predicts a unique spectral signature in the gamma-ray spectrum, characterized by a fast rising spectrum and a spectral break below a few hundreds of MeV. We here report the evidence of a spectral break around 500 MeV in the disk emission of Large Magellanic Cloud (LMC), which is found in the analysis of the gamma-ray data extending down to 60 MeV observed by {\it Fermi}-Large Area Telescope. The break is well consistent with the pion-decay model for the gamma-ray emission, although leptonic models, such as the electron bremsstrahlung emission, cannot be ruled out completely.

The non-repeating fast radio bursts (FRBs) could arise from catastrophic stellar explosions or magnetar giant flares, so relativistic blast waves might be produced in these events. Motivated by this, we here search for GeV counterparts to all non-repeating FRBs with Fermi Large Area Telescope (LAT), including FRB 131104 that is claimed to be possibly associated with a gamma-ray transient candidate detected by Swift Burst Alert Telescope (BAT). FRB 131104 enters the field of view (FoV) of LAT ~5000s after the burst time, so we are only able to search for the GeV afterglow emission during this period, but no significant GeV emission is detected. We also perform a search for GeV emissions from other FRBs, but no significant GeV emissions are detected either. Upper limit fluences in the range of (4.7-29.2)x10^(-7) erg/cm2 are obtained, and then the upper limits of the isotropic blast wave kinetic energy of about (1-200)x10^(53) erg are inferred under certain assumptions. Although the current limits on the isotropic blast wave energy are not sufficiently stringent to rule out the connection between FRBs and GRB-like transients, future more sensitive observations with Fermi or Imaging Atmospheric Cherenkov Telescopes might be able to constrain the connection.

Dozens of high-energy neutrinos have been detected by the IceCube neutrino telescope, but no clear association with any classes of astrophysical sources has been identified so far. Recently, Kadler et al. (2016) report that a PeV cascade-like neutrino event occurs in positional and temporal coincidence with a giant gamma-ray flare of the blazar PKS B1424-418. Since IceCube track-like events have much better angular resolution, we here search for possible short-term gamma-ray flares that are associated with the IceCube track-like events with Fermi Large Area Telescope (LAT) observations. Among them, three track-like neutrino events occur within the field of view of Fermi-LAT at the time of the detection, so search for the {\em prompt} gamma-ray emission associated with neutrinos are possible. Assuming a point source origin and a single power law spectrum for the possible gamma-ray sources associated with neutrinos, a likelihood analysis of 0.2-100 GeV photons observed by Fermi-LAT on the timescales of $\sim 12$ hours and one year are performed, and for the three special neutrinos, the analysis are also performed on the timescales of thousand of seconds before and after the neutrino detection. No significant GeV excesses over the background are found and the upper limit fluxes at 95\% confidence level are obtained for different timescales. We also search for possible hard X-ray transient sources associated with the IceCube track-like neutrino events, but the search also yields null results. We discuss the implication of the non-detection of gamma-ray flares for the constraints on the neutrino source density.

Massive black holes at galaxy center may tear apart a star when the star passes occasionally within the disruption radius, which is the so-called tidal disruption event(TDE). Most TDEs radiate with thermal emission resulted from the accretion disk, but three TDEs have been detected in bright non-thermal X-ray emission, which is interpreted as arising from the relativistic jets. Search for high-energy gamma-ray emission from one relativistic TDE (Swift J164449.3+573451) with the \textsl{Fermi} Large Area Telescope (LAT) has yielded non-detection. In this paper, we report the search for high energy emission from the other two relativistic TDEs (Swift J2058.4+0516 Swift J1112.2-8238) during the flare period. No significant GeV emission is found, with an upper limit fluence in LAT energy range being less than $1\%$ of that in X-rays. Compared with gamma-ray bursts (GRBs) and blazars, these TDEs have the lowest flux ratio between GeV emission and X-ray emission. The non-detection of high-energy emission from relativistic TDEs could be due to that the high-energy emission is absorbed by soft photons in the source. Based on this hypothesis, upper limits on the bulk Lorentz factors, $\Gamma\la 30$, are then obtained for the jets in these TDEs. We also search for high-energy gamma-ray emission from the nearest TDE discovered to date, ASASSN-14li. No significant GeV emission is found and an upper limit of $L(\rm 0.1-10 GeV)\le 4.4\times 10^{42}$ erg s$^{-1}$ (at $95\%$ confidence level) is obtained for the first $10^{7}$ s after the disruption.

Cosmic rays (CRs) in starburst galaxies produce high energy gamma-rays by colliding with the dense interstellar medium (ISM). Arp 220 is the nearest ultra luminous infrared galaxy (ULIRG) that has star-formation at extreme levels, so it has long been predicted to emit high-energy gamma-rays. However, no evidence of gamma-ray emission was found despite intense efforts of search. Here we report the discovery of high-energy gamma-ray emission above 200 MeV from Arp 220 at a confidence level of $\sim 6.3 \sigma $ using 7.5 years of \textsl {Fermi} Large Area Telescope observations. The gamma-ray emission shows no significant variability over the observation period and it is consistent with the quasi-linear scaling relation between the gamma-ray luminosity and total infrared luminosity for star-forming galaxies, suggesting that these gamma-rays arise from CR interactions. As the high density medium of Arp 220 makes it an ideal CR calorimeter, the gamma-ray luminosity can be used to measure the efficiency of powering CRs by supernova (SN) remnants given a known supernova rate in Arp 220. We find that this efficiency is about $4.2\pm2.6\%$ for CRs above 1 GeV.

Gamma-ray bursts (GRBs) are powered by ultra-relativistic jets. Usually a minimum value of the Lorentz factor of the relativistic bulk motion is obtained based on the argument that the observed high energy photons ($\gg {\rm MeV}$) can escape without suffering from absorption due to pair production. The exact value, rather than a lower limit, of the Lorentz factor can be obtained if the spectral cutoff due to such absorption is detected. With the good spectral coverage of the Large Area Telescope (LAT) on \textsl{Fermi}, measurements of such cutoff become possible, and two cases (GRB~090926A and GRB~100724B) have been reported to have high-energy cutoffs or breaks. We systematically search for such high energy spectral cutoffs/breaks from the LAT and the Gamma-ray burst monitor (GBM) observations of the prompt emission of GRBs detected since August 2011. Six more GRBs are found to have cutoff-like spectral feature at energies of $\sim10-500$~MeV. Assuming that these cutoffs are caused by pair-production absorption within the source, the bulk Lorentz factors of these GRBs are obtained. We further find that the Lorentz factors are correlated with the isotropic gamma-ray luminosity of the bursts, indicating that more powerful GRB jets move faster.

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $\Gamma_{\rm ph}$ of the outflow. The derived $\Gamma_{\rm ph}$ range from $50$ to $150$ and a relation of $\Gamma_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $\Gamma_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs.

GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star.

We investigate the effect of a new triple-alpha reaction rate from Ogata et al. (2009) on helium ignition conditions on accreting neutron stars and on the properties of the subsequent type I X-ray burst. We find that the new rate leads to significantly lower ignition column density for accreting neutron stars at low accretion rates. We compare the results of our ignition models for a pure helium accretor to observations of bursts in ultra-compact X-ray binary (UCXBs), which are believed to have nearly pure helium donors. For mdot > 0.001 mdot_Edd, the new triple-alpha reaction rate from Ogata et al. (2009) predicts a maximum helium ignition column of ~ 3 x 10^9 g cm^{-2}, corresponding to a burst energy of ~ 4 x 10^{40} ergs. For mdot ~ 0.01 mdot_Edd at which intermediate long bursts occur, the predicted burst energies are at least a factor of 10 too low to explain the observed energies of such bursts in UCXBs. This finding adds to the doubts cast on the triple-alpha reaction rate of Ogata et al. (2009) by the low-mass stellar evolution results of Dotter & Paxton (2009).

Recent observations with Galaxy Evolution Explorer (GALEX) show strong unexpected UV excess in the spectrum of brightest cluster galaxies (BCGs). It is believed that the excess UV signal is produced by old and evolved core-He burning stars, and the UV flux strength could be greatly enhanced if the progenitor stars have high value of He abundance. In this work, we propose that sedimentation process can greatly enhance the He abundance in BCGs. Our model predicts that the UV flux strength is stronger in more massive, low-redshift, and dynamically relaxed BCGs. These predictions are testable with the current generation of GALEX+SDSS observations.

The uniformity of the helium-to-hydrogen abundance ratio in X-ray emitting intracluster medium (ICM) is one of the commonly adopted assumptions in X-ray analyses of galaxy clusters and cosmological constraints derived from these measurements. In this work, we investigate the effect of He sedimentation on X-ray measurements of galaxy clusters in order to assess this assumption and associated systematic uncertainties. By solving a set of flow equations for a H-He plasma, we show that the helium-to-hydrogen mass ratio is significantly enhanced in the inner regions of clusters. The effect of He sedimentation, if not accounted for, introduces systematic biases in observable properties of clusters derived using X-ray observations. We show that these biases also introduce an apparent evolution in the observed gas mass fractions of X-ray luminous, dynamically relaxed clusters and hence biases in observational constraints on the dark energy equation of state parameter, w, derived from the cluster distance-redshift relation. The Hubble parameter derived from the combination of X-ray and Sunyaev-Zel'dovich effect (SZE) measurements is affected by the He sedimentation process as well. Future measurements aiming to constrain w or H_0 to better than 10% may need to take into account the effect of He sedimentation. We propose that the evolution of gas mass fraction in the inner regions of clusters should provide unique observational diagnostics of the He sedimentation process.

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Recent observations of unstable thermonuclear burning (observed as X-ray bursts) on the surfaces of slowly accreting neutron stars ($< 0.01$ of the Eddington rate) motivate us to examine how sedimentation of CNO isotopes affects the ignition of these bursts. We further estimate the burst development using a simple one-zone model with a full reaction network. We report a region of mass accretion rates for weak H flashes. Such flashes can lead to a large reservoir of He, the unstable burning of which may explain some observed long bursts (duration $\sim 1000$ s).

We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Motivated by observations of Type I X-ray bursts from sources with extremely low persistent accretion luminosities, $L_X < 10^{36}\usp\ergspersecond (\simeq 0.01\ensuremath{L_{\mathrm{Edd}}}$), we study how sedimentation affects the distribution of isotopes and the ignition of H and He in the envelope of an accreting neutron star. For local mass accretion rates $\mdot \lesssim 10^{-2}\medd$ (for which the ignition of H is unstable), where $\medd = 8.8\times 10^{4}\nsp\gpscps$, the helium and CNO elements sediment out of the accreted fuel before reaching a temperature where H would ignite. Using one-zone calculations of the thermonuclear burning, we find a range of accretion rates for which the unstable H ignition does not trigger unstable He burning. This range depends on the emergent flux from reactions in the deep neutron star crust; for $F = 0.1\nsp\MeV(\dot{m}/\mb)$, the range is $3\times 10^{-3}\medd\lesssim\mdot\lesssim 10^{-2}\medd$. We speculate that sources accreting in this range will build up a massive He layer that later produces an energetic and long X-ray burst. At mass accretion rates lower than this range, we find that the H flash leads to a strong mixed H/He flash. Surprisingly, even at accretion rates $\mdot \gtrsim 0.1\medd$, although the H and He do not completely segregate, the H abundance at the base of the accumulated layer is still reduced. While following the evolution of the X-ray burst is beyond the scope of this introductory paper, we note that the reduced proton-to-seed ratio favors the production of \iso{12}{C}--an important ingredient for subsequent superbursts.

Recent observational and theoretical studies have raised the possibility that the collimated outflows in gamma-ray burst (GRB) sources have two distinct components: a narrow (opening half-angle $\theta_{\rm n}$), highly relativistic (initial Lorentz factor $\eta_\rmn \gtrsim 10^2$) outflow, from which the $\gamma$-ray emission originates, and a wider ($\theta_{\rm w} \lesssim 3 \theta_{\rm n}$), moderately relativistic ($\eta_{\rm w}\sim 10$) surrounding flow. Using a simple synchrotron emission model, we calculate the R-band afterglow lightcurves expected in this scenario and derive algebraic expressions for the flux ratios of the emission from the two jet components at the main transition times in the lightcurve. We apply this model to GRB sources, for explaining the structure of afterglows and source energetics, as well as to X-ray flash sources, which we interpret as GRB jets viewed at an angle $\theta_{\rm obs} > \theta_{\rm n}$. Finally, we argue that a neutron-rich hydromagnetic outflow may naturally give rise to repeated brightening episodes in the afterglow lightcurve as observed in GRB 021004 and GRB 030329.

We demonstrate that ``hot'' MHD outflows from neutron-rich black-hole debris disks can significantly alleviate the baryon-loading problem in gamma-ray burst (GRB) sources. We argue that the neutron-to-proton ratio in disk-fed outflows might be as high as ~30 and show, with the help of an exact semianalytic relativistic-MHD solution, that the neutrons can decouple at a Lorentz factor gamma_d ~ 15 even as the protons continue to accelerate to Lorentz factor ~200 and end up acquiring ~30% of the injected energy. We clarify the crucial role that the magnetic field plays in this process and prove that purely hydrodynamic outflows must have gamma_d > a few hundreds. The motion of the decoupled neutrons is not collinear with that of the decoupled protons, so, in contrast to previous suggestions based on purely hydrodynamic models, the two particle groups do not collide after decoupling. If the decoupled neutrons move at an angle > 1/gamma_d = 3.8 degr (15/gamma_d) to the line of sight to the GRB source, most of their emission after they decay into protons will remain unobservable.