Search

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data., Comment: Accepted for publication in Physical Review D

We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final kinetic energy and nickel-56 yields conform better to observational values than is the case for the "classical" GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel), the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that are necessary for the Zel'dovich gradient mechanism to produce a detonation. We demonstrate feasibility of this scenario with three 3-dimensional (3D), full star simulations of this model using the FLASH code. We characterized the simulations by the energy released during the deflagration phase, which ranged from 38% to 78% of the white dwarf's binding energy. We show that the necessary conditions for detonation are achieved in all three of the models., Comment: 22 pages, 8 figures; ApJ

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background., Comment: 28 pages, 20 figures; accepted for publication in Phys. Rev. D; Contact authors: A. Cuoco, E. Komatsu, T. Linden, M. N. Mazziotta, J. Siegal-Gaskins, V. Vitale

GRB 091127 is a bright gamma-ray burst (GRB) detected by Swift at a redshift z=0.49 and associated with SN 2009nz. We present the broadband analysis of the GRB prompt and afterglow emission and study its high-energy properties in the context of the GRB/SN association. While the high luminosity of the prompt emission and standard afterglow behavior are typical of cosmological long GRBs, its low energy release, soft spectrum and unusual spectral lag connect this GRB to the class of sub-energetic bursts. We discuss the suppression of high-energy emission in this burst, and investigate whether this behavior could be connected with the sub-energetic nature of the explosion., Comment: 12 pages, 7 figures, submitted to ApJ

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel., Comment: 29 pages, 5 figures, 3 tables; Accepted by ApJ; Contact authors: Alex Drlica-Wagner, Ping Wang, Louie Strigari, and Elliott Bloom

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009., Comment: published in A&A, 9pages, 4figures. Due to some style difficulties using the standard style file with the long author list, the abstract has been included manually. A version with the correct style can be found in the website http://www.lsw.uni-heidelberg.de/users/skaufman/publications/aa16170-10published.pdf

Spectral lag, which is defined as the difference in time of arrival of high and low energy photons, is a common feature in Gamma-ray Bursts (GRBs). Previous investigations have shown a correlation between this lag and the isotropic peak luminosity for long duration bursts. However, most of the previous investigations used lags extracted in the observer-frame only. In this work (based on a sample of 43 Swift long GRBs with known redshifts), we present an analysis of the lag-luminosity relation in the GRB source-frame. Our analysis indicates a higher degree of correlation -0.82 +/- 0.05 (chance probability of ~ 5.5 x 10^-5) between the spectral lag and the isotropic peak luminosity, Liso, with a best-fit power-law index of -1.2 +/- 0.2, such that Liso proportional to lag^-1.2. In addition, there is an anti-correlation between the source-frame spectral lag and the source-frame peak energy of the burst spectrum, E_pk(1+z)., Comment: 11 pages, 6 figures, 6 table; Accepted for publication in MNRAS

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy., Comment: 5 figures, 1 table, revtex 4.1, updated to match PRL published version

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors., Comment: 6 pages, 2 figures; Contact authors: Johann Cohen-Tanugi, Jan Conrad, and Maja Llena Garde

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance., Comment: Accepted for publications by ApJ. Contact authors: M. Giroletti, V. Pavlidou, A. Reimer

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample., Comment: 22 pages, 11 figures, 4 tables; accepted for publication in the Astrophysical Journal

We develop, implement and characterise an enhanced data reduction approach which delivers precise, accurate, radial velocities from moderate resolution spectroscopy with the fibre-fed VLT/FLAMES+GIRAFFE facility. This facility, with appropriate care, delivers radial velocities adequate to resolve the intrinsic velocity dispersions of the very faint dSph dwarf galaxies. Importantly, repeated measurements let us reliably calibrate our individual velocity errors ($0.2 \leq \delta_V\leq 5$ km s$^{-1}$) and directly detect stars with variable radial velocities. We show, by application to the Bootes-1 dwarf spheroidal, that the intrinsic velocity dispersion of this system is significantly below 6.5\,km/s reported by previous studies. Our data favor a two-population model of Bootes-1, consisting of a majority `cold' stellar component, with velocity dispersion $2.4^{+0.9}_{-0.5}$\,km/s, and a minority `hot' stellar component, with velocity dispersion $\sim 9$\,km/s, although we can not completely rule out a single component distribution with velocity dispersion $4.6^{0.8}_{-0.6}$\,km/s. We speculate this complex velocity distribution actually reflects the distribution of velocity anisotropy in Bootes-1, which is a measure of its formation processes., Comment: ApJ in press. Anticipated publication July 1, 2011

We present a detailed chemical composition analysis of 35 red giant stars in the globular cluster M22. High resolution spectra for this study were obtained at five observatories, and analyzed in a uniform manner. We have determined abundances of representative light proton-capture, alpha, Fe-peak and neutron-capture element groups. Our aim is to better understand the peculiar chemical enrichment history of this cluster, in which two stellar groups are characterized by a different content in iron, neutron capture elements Y, Zr and Ba, and alpha element Ca. The principal results of this study are: (i) substantial star-to-star metallicity scatter (-2.0<[Fe/H]<-1.6); (ii) enhancement of s-process/r-process neutron-capture abundance ratios in a fraction of giants, positively correlated with metallicity; (iii) sharp separation between the s-process rich and s-process poor groups by [La/Eu] ratio; (iv) possible increase of [Cu/Fe] ratios with increasing [Fe/H], suggesting that this element also has a significant s-process component; and (v) presence of Na-O and C-N anticorrelations in both the stellar groups., Comment: 26 pages, 19 figures. Accepted for publication in A&A

(Abridged) Carbon-enhanced metal-poor (CEMP) stars in the halo components of the Milky Way are explored, based on accurate determinations of the carbon-to-iron ([C/Fe]) abundance ratios and kinematic quantities for over 30000 calibration stars from the Sloan Digital Sky Survey (SDSS). Using our present criterion that low-metallicity stars exhibiting [C/Fe] ratios ("carbonicity") in excess of [C/Fe]$ = +0.7$ are considered CEMP stars, the global frequency of CEMP stars in the halo system for \feh\ $< -1.5$ is 8%; for \feh\ $< -2.0$ it is 12%; for \feh\ $<-2.5$ it is 20%. We also confirm a significant increase in the level of carbon enrichment with declining metallicity, growing from $<$[C/Fe]$>$ $\sim +1.0$ at \feh\ $= -1.5$ to $<$[C/Fe]$>$ $\sim +1.7$ at \feh\ $= -2.7$. The nature of the carbonicity distribution function (CarDF) changes dramatically with increasing distance above the Galactic plane, $|$Z$|$. For $|$Z$|$ $< 5$ kpc, relatively few CEMP stars are identified. For distances $|$Z$|$ $> 5$ kpc, the CarDF exhibits a strong tail towards high values, up to [C/Fe] $>$ +3.0. We also find a clear increase in the CEMP frequency with $|$Z$|$. For stars with $-2.0 <$ [Fe/H] $< -$1.5, the frequency grows from 5% at $|$Z$|$ $\sim 2$ kpc to 10% at $|$Z$|$ $\sim 10$ kpc. For stars with [Fe/H] $< -$2.0, the frequency grows from 8% at $|$Z$|$ $\sim 2$ kpc to 25% at $|$Z$|$ $\sim 10$ kpc. For stars with $-2.0 <$ [Fe/H] $< -$1.5, the mean carbonicity is $<$[C/Fe]$>$ $\sim +1.0$ for 0 kpc $<$ $|$Z$|$ $<$ 10 kpc, with little dependence on $|$Z$|$; for [Fe/H] $< -$2.0, $<$[C/Fe]$>$ $\sim +1.5$, again roughly independent of $|$Z$|$., Comment: Accepted for publication in the Astrophysical Journal, 32 pages, 15 figures

We have studied power density spectra (PDS) of 206 long Gamma-Ray Bursts (GRBs). We fitted the PDS with a simple power-law and extracted the exponent of the power-law (alpha) and the noise-crossing threshold frequency (f_th). We find that the distribution of the extracted alpha peaks around -1.4 and that of f_th around 1 Hz. In addition, based on a sub-set of 58 bursts with known redshifts, we show that the redshift-corrected threshold frequency is positively correlated with the isotropic peak luminosity. The correlation coefficient is 0.57 +/- 0.03., Comment: 9 pages, 17 figures, 1 table; Accepted for publication in MNRAS

We determine the metallicity distribution function (MDF) of the Galactic halo based on metal-poor main-sequence turnoff-stars (MSTO) which were selected from the Hamburg/ESO objective-prism survey (HES) database. Corresponding follow-up moderateresolution observations (R ~ 2000) of some 682 stars (among which 617 were accepted program stars) were carried out with the 2.3m telescope at the Siding Spring Observatory (SSO). Corrections for the survey volume covered by the sample stars were quantitatively estimated and applied to the observed MDF. The corrections are quite small, when compared with those for a previously studied sample of metal-poor giants. The corrected observational MDF of the turnoff sample was then compared with that of the giants, as well as with a number of theoretical predictions of Galactic chemical evolution, including the mass-loss modified Simple Model. Although the survey-volume corrected MDFs of the metal-poor turnoff and the halo giants notably differ in the region of [Fe/H] > -2.0, below [Fe/H] ~ -2.0, (the region we scientifically focus on most) both MDFs show a sharp drop at [Fe/H] ~ -3.6 and present rather similar distributions in the low-metallicity tail. Theoretical models can fit some parts of the observed MDF, but none is found to simultaneously reproduce the peak as well as the features in the metal-poor region with [Fe/H] between -2.0 to -3.6. Among the tested models only the GAMETE model, when normalized to the tail of the observed MDF below [Fe/H] ~ -3.0, and with Z_{cr} = 10^{-3.4}Z_{\odot}, is able to predict the sharp drop at [Fe/H] ~ -3.6., Comment: 10 pages, 11 figures, accepted for publication in A&A

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

One of the most important factors in determining the stellar lithium abundance is the effective temperature. In a previous study by the authors, new effective temperatures Teff for sixteen metal-poor halo dwarfs were derived using a local thermodynamic equilibrium (LTE) description of the formation of Fe lines. This new Teff scale reinforced the discrepancy. For six of the stars from our previous study we calculate revised temperatures using a non-local thermodynamic equilibrium (NLTE) approach. These are then used to derive a new mean primordial lithium abundance in an attempt to solve the lithium discrepancy. Using the code MULTI we calculate NLTE corrections to the LTE abundances for the Fe I lines measured in the six stars, and determine new Teff's. We keep other physical parameters, i.e. log g, [Fe/H] and xi, constant at the values calculated in Paper I. With the revised Teff scale we derive new Li abundances. We compare the NLTE values of Teff with the photometric temperatures of Ryan et al. (1999, ApJ, 523, 654), the infrared flux method (IRFM) temperatures of Melendez & Ramirez (2004, ApJ, 615, 33), and the Balmer line wing temperatures of Asplund et al. (2006, ApJ, 644, 229). We find that our temperatures are hotter than both the Ryan et al. and Asplund et al. temperatures by typically ~ 110 K - 160 K, but are still cooler than the temperatures of Melendez & Ramirez by typically ~ 190 K. The temperatures imply a primordial Li abundance of 2.19 dex or 2.21 dex, depending on the magnitude of collisions with hydrogen in the calculations, still well below the value of 2.72 dex inferred from WMAP + BBN. We discuss the effects of collisions on trends of 7Li abundances with [Fe/H] and Teff, as well as the NLTE effects on the determination of log g through ionization equilibrium, which imply a collisional scaling factor SH > 1 for collisions between Fe and H atoms., Comment: 12 pages, 10 Figures

The volume of a Wiener sausage constructed from a diffusion process with periodic, mean-zero, divergence-free velocity field, in dimension 3 or more, is shown to have a non-random and positive asymptotic rate of growth. This is used to establish the existence of a homogenized limit for such a diffusion when subject to Dirichlet conditions on the boundaries of a sparse and independent array of obstacles. There is a constant effective long-time loss rate at the obstacles. The dependence of this rate on the form and intensity of the obstacles and on the velocity field is investigated. A Monte Carlo algorithm for the computation of the volume growth rate of the sausage is introduced and some numerical results are presented for the Taylor--Green velocity field., Comment: 19 pages, 7 figures, To appear in Bingham, N. H., and Goldie, C. M. (eds), Probability and Mathematical Genetics: Papers in Honour of Sir John Kingman. London Math. Soc. Lecture Note Series. Cambridge: Cambridge Univ. Press

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals., Comment: 29 pages, 7 figures, to be published in JCAP. Corresponding authors: Jan Conrad, Michael Gustafsson, Alexander Sellerholm, Gabrijela Zaharijas

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed., Comment: 6 pages, 1 figure. Accepted for publication by The Physical Review Letters

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged), Comment: 25 pages, 4 figures, accepted to ApJ, Corresponding authors: J. Cohen-Tanugi, C. Farnier, T.E. Jeltema, E. Nuss, and S. Profumo

We set out a general procedure which allows the approximation of certain Markov chains by the solutions of differential equations. The chains considered have some components which oscillate rapidly and randomly, while others are close to deterministic. The limiting dynamics are obtained by averaging the drift of the latter with respect to a local equilibrium distribution of the former. Some general estimates are proved under a uniform mixing condition on the fast variable which give explicit error probabilities for the fluid approximation. Mitzenmacher, Prabhakar and Shah [In Proc. 43rd Ann. Symp. Found. Comp. Sci. (2002) 799-808, IEEE] introduced a variant with memory of the "join the shortest queue" or "supermarket" model, and obtained a limit picture for the case of a stable system in which the number of queues and the total arrival rate are large. In this limit, the empirical distribution of queue sizes satisfies a differential equation, while the memory of the system oscillates rapidly and randomly. We illustrate our general fluid limit estimate by giving a proof of this limit picture., Comment: Published in at http://dx.doi.org/10.1214/12-AAP861 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Spectral lag, the time difference between the arrival of high-energy and low-energy photons, is a common feature in Gamma-ray Bursts (GRBs). Norris et al. 2000 reported a correlation between the spectral lag and the isotropic peak luminosity of GRBs based on a limited sample. More recently, a number of authors have provided further support for this correlation using arbitrary energy bands of various instruments. In this paper we report on a systematic extraction of spectral lags based on the largest Swift sample to date of 31 GRBs with measured redshifts. We extracted the spectral lags for all combinations of the standard Swift hard x-ray energy bands: 15-25 keV, 25-50 keV, 50-100 keV and 100-200 keV and plotted the time dilation corrected lag as a function of isotropic peak luminosity. The mean value of the correlation coefficient for various channel combinations is -0.68 with a chance probability of ~ 0.7 x 10^{-3}. In addition, the mean value of the power-law index is 1.4 +/- 0.3. Hence, our study lends support for the existence of a lag-luminosity correlation, albeit with large scatter., Comment: 19 Pages, 11 Figures and 5 Tables; Accepted to The Astrophysical Journal

One of the most prominent, yet controversial associations derived from the ensemble of prompt-phase observations of gamma-ray bursts (GRBs) is the apparent correlation in the source frame between the peak energy Epeak) of the nu-F(nu) spectrum and the isotropic radiated energy, Eiso. Since most gamma-ray bursts (GRBs) have Epeak above the energy range (15-150 keV) of the Burst Alert Telescope (BAT) on Swift, determining accurate Epeak values for large numbers of Swift bursts has been difficult. However, by combining data from Swift/BAT and the Suzaku Wide-band All-Sky Monitor (WAM), which covers the energy range from 50-5000 keV, for bursts which are simultaneously detected, one can accurately fit Epeak and Eiso and test the relationship between them for the Swift sample. Between the launch of Suzaku in July 2005 and the end of April 2009, there were 48 gamma-ray bursts (GRBs) which triggered both Swift/BAT and WAM and an additional 48 bursts which triggered Swift and were detected by WAM, but did not trigger. A BAT-WAM team has cross-calibrated the two instruments using GRBs, and we are now able to perform joint fits on these bursts to determine their spectral parameters. For those bursts with spectroscopic redshifts, we can also calculate the isotropic energy. Here we present the results of joint Swift/BAT-Suzaku/WAM spectral fits for 91 of the bursts detected by the two instruments. We show that the distribution of spectral fit parameters is consistent with distributions from earlier missions and confirm that Swift bursts are consistent with earlier reported relationships between Epeak and isotropic energy. We show through time-resolved spectroscopy that individual burst pulses are also consistent with this relationship., Comment: Accepted for publication in the Astrophysical Journal

We apply the optimal filter technique to Sloan Digital Sky Survey photometry around Segue 1 and find that the outer parts of the cluster are distorted. There is strong evidence for ~ 1 degree elongations of extra-tidal stars, extending both eastwards and southwestwards of the cluster. The extensions have similar differential Hess diagrams to Segue 1. A Kolmogorov-Smirnov test suggests a high probability that both come from the same parent distribution. The location of Segue 1 is close to crossings of the tidal wraps of the Sagittarius stream. By extracting blue horizontal branch stars from Sloan's spectral database, two kinematic features are isolated and identified with different wraps of the Sagittarius stream. We show that Segue 1 is moving with a velocity that is close to one of the wraps. At this location, we estimate that there are enough Sagittarius stars, indistinguishable from Segue 1 stars, to inflate the velocity dispersion and hence the mass-to-light ratio. All the available evidence is consistent with the interpretation that Segue 1 is a star cluster, originally from the Sagittarius galaxy, and now dissolving in the Milky Way., Comment: MNRAS, in press, 13 figures, 2 tables

The LAT instrument on the Fermi mission will reveal the rich spectral and temporal gamma-ray burst phenomena in the > 100 MeV band. The synergy with Fermi's GBM detectors will link these observations to those in the well explored 10-1000 keV range; the addition of the > 100 MeV band observations will resolve theoretical uncertainties about burst emission in both the prompt and afterglow phases. Trigger algorithms will be applied to the LAT data both onboard the spacecraft and on the ground. The sensitivity of these triggers will differ because of the available computing resources onboard and on the ground. Here we present the LAT's burst detection methodologies and the instrument's GRB capabilities., Comment: Accepted by ApJ

The integration-by-parts formula discovered by Malliavin for the Ito map on Wiener space is proved using the two-parameter stochastic calculus. It is also shown that the solution of a one-parameter stochastic differential equation driven by a two-parameter semimartingale is itself a two-parameter semimartingale., Comment: 19 pages

Gamma-ray bursts (GRBs) have been an enigma since their discovery forty years ago. However, considerable progress unraveling their mysteries has been made in recent years. Developments in observations, theory, and instrumentation have prepared the way so that the next decade can be the one in which we finally answer the question, "What are gamma-ray bursts?" This question encompasses not only what the progenitors are that produce the GRBs, but also how the enormous luminosity of the GRBs, concentrated in gamma rays, is achieved. Observations across the electromagnetic spectrum, from both the ground and space, will be required to fully tackle this important question. This white paper, mostly distilled from a recent study commissioned by the Division of Astrophysics of the American Physical Society, focuses on what very high energy (~100 GeV and above) gamma-ray observations can contribute. Very high energy gamma rays probe the most extreme high energy particle populations in the burst environment, testing models of lepton and proton acceleration in GRBs and constraining the bulk Lorentz factor and opacity of the outflow. Sensitivity improvements of more than an order of magnitude in the very high energy gamma-ray band can be achieved early in the next decade, in order to contribute to this science., Comment: White paper contributed to the Stars and Stellar Evolution Science Frontier Panel of the Astro2010 Decadal Survey

The discovery of the Spite plateau in the abundances of 7Li for metal-poor stars led to the determination of an observationally deduced primordial lithium abundance. However, after the success of the Wilkinson Microwave Anisotropy Probe (WMAP) in determining the baryon density, OmegaBh^2, there was a discrepancy between observationally determined and theoretically determined abundances in the case of 7Li. One of the most important uncertain factors in the calculation of the stellar 7Li abundance is the effective temperature, Teff. We use sixteen metal-poor halo dwarfs to calculate new Teff values using the excitation energy method. With this temperature scale we then calculate new Li abundances for this group of stars in an attempt to resolve the 7Li discrepancy. Using high signal-to-noise (S/N ~ 100) spectra of 16 metal-poor halo dwarfs, obtained with the UCLES spectrograph on the AAT, measurements of equivalent widths from a set of unblended FeI lines are made. These equivalent widths are then used to calculate new Teff values with the use of the single line radiative transfer program WIDTH6, where we have constrained the gravity using either theoretical isochrones or the Hipparcos parallax, rather than the ionization balance. The lithium abundances of the stars are calculated with these temperatures. The physical parameters are derived for the 16 programme stars, and two standards. These include Teff, log g, [Fe/H], microturbulence and 7Li abundances. A comparison between the temperature scale of this work and those adopted by others has been undertaken. We find good consistency with the temperatures derived from the Halpha line by Asplund et al. (2006), but not with the hotter scale of Melendez & Ramirez (2004)., Comment: 13 pages, 9 figures

We report on observations of GRB 080503, a short gamma-ray burst with very bright extended emission (about 30 times the gamma-ray fluence of the initial spike) in conjunction with a thorough comparison to other short Swift events. In spite of the prompt-emission brightness, however, the optical counterpart is extraordinarily faint, never exceeding 25 mag in deep observations starting at ~1 hr after the BAT trigger. The optical brightness peaks at ~1 day and then falls sharply in a manner similar to the predictions of Li & Paczynski (1998) for supernova-like emission following compact-binary mergers. However, a shallow spectral index and similar evolution in X-rays inferred from Chandra observations are more consistent with an afterglow interpretation. The extreme faintness of this probable afterglow relative to the bright gamma-ray emission argues for a very low-density medium surrounding the burst (a "naked" GRB), consistent with the lack of a coincident host galaxy down to 28.5 mag in deep HST imaging. Our observations reinforce the notion that short GRBs generally occur outside regions of active star formation, but demonstrate that in some cases the luminosity of the extended prompt emission can greatly exceed that of the short spike, which may constrain theoretical interpretation of this class of events. Because most previous BAT short bursts without observed extended emission are too faint for this signature to have been detectable even if it were present at typical level, conclusions based solely on the observed presence or absence of extended emission in the existing Swift sample are premature. (abridged), Comment: 16 pages, 11 figures; accepted for publication in ApJ

This is a report on the findings of the gamma ray burst working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe gamma ray bursts at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that major advances are possible and that the detection of very high energy emission would have strong implications for GRB models, as well as cosmic ray origin., Comment: report from the Gamma Ray Burst Working group of the APS commissioned White paper on ground-based TeV gamma ray astronomy

We determine the metallicity distribution function (MDF) of the Galactic halo by means of a sample of 1638 metal-poor stars selected from the Hamburg/ESO objective-prism survey (HES). The sample was corrected for minor biases introduced by the strategy for spectroscopic follow-up observations of the metal-poor candidates, namely "best and brightest stars first". [...] We determined the selection function of the HES, which must be taken into account for a proper comparison between the HES MDF with MDFs of other stellar populations or those predicted by models of Galactic chemical evolution. The latter show a reasonable agreement with the overall shape of the HES MDF for [Fe/H] > -3.6, but only a model of Salvadori et al. (2007) with a critical metallicity for low-mass star formation of Z_cr = 10^{-3.4} * Z_Sun reproduces the sharp drop at [Fe/H] ~-3.6 present in the HES MDF. [...] A comparison of the MDF of Galactic globular clusters and of dSph satellites to the Galaxy shows qualitative agreement with the halo MDF, derived from the HES, once the selection function of the latter is included. However, statistical tests show that the differences between these are still highly significant. [ABSTRACT ABRIDGED], Comment: Accepted for publication in A&A

Using a sample of 248 metal-poor stars (RR Lyraes, red giants and RHB stars) which is remarkable for the accuracy of its 6-D kinematical data, we find a new component for the local halo which has an axial ratio c/a ~ 0.2, a similar flattening to the thick disk. It has a small prograde rotation but is supported by velocity anisotropy, and contains more intermediate-metallicity stars (with -1.5 < [Fe/H] < -1.0) than the rest of our sample. We suggest that this component was formed quite late, during or after the formation of the disk. It formed either from the gas that was accreted by the last major mergers experienced by the Galaxy, or by dynamical friction of massive infalling satellite(s) with the halo and possibly the stellar disk or thick disk. The remainder of the stars in our sample exhibit a clumpy distribution in energy and angular momentum, suggesting that the early, chaotic conditions under which the inner halo formed were not violent enough to erase the record of their origins. The clumpy structure suggests that a relatively small number of progenitors were responsible for building up the inner halo, in line with theoretical expectations. We find a difference in mean binding energy between the RR Lyrae variables and the red giants in our sample, suggesting that more of the RR Lyraes in the sample belong to the outer halo, and that the outer halo may be somewhat younger, as first suggested by Searle and Zinn (1978). We also find that the RR Lyrae mean rotation is more negative than the red giants, which is consistent with the recent result of Carollo et al.(2007) that the outer halo has a retrograde rotation and with the difference in kinematics seen between RR Lyraes and BHB stars by Kinman et al.(2007)., Comment: 16 pages, 10 figures, this version accepted by ApJ

This is a short report on the preliminary findings of the gamma ray burst (GRB) working group for the white paper on the status and future of very high energy (VHE; >50 GeV) gamma-ray astronomy. The white paper discusses the status of past and current attempts to observe GRBs at GeV-TeV energies, including a handful of low-significance, possible detections. The white paper concentrates on the potential of future ground-based gamma-ray experiments to observe the highest energy emission ever recorded for GRBs, particularly for those that are nearby and have high Lorentz factors in the GRB jet. It is clear that the detection of VHE emission would have strong implications for GRB models, as well as cosmic ray origin. In particular, the extended emission phase (including both afterglow emission and possible flaring) of nearby long GRBs could provide the best possibility for detection. The difficult-to-obtain observations during the prompt phase of nearby long GRBs and short GRBs could also provide particularly strong constraints on the opacity and bulk Lorentz factors surrounding the acceleration site. The synergy with upcoming and existing observatories will, of course, be critical for both identification of GRBs and for multiwavelength/multimessenger studies., Comment: to appear in Proceedings of Gamma Ray Bursts 2007 Meeting, Santa Fe, New Mexico; 5 pages, 4 figures

We have gathered optical photometry data from the literature on a large sample of Swift-era gamma-ray burst (GRB) afterglows including GRBs up to September 2009, for a total of 76 GRBs, and present an additional three pre-Swift GRBs not included in an earlier sample. Furthermore, we publish 840 additional new photometry data points on a total of 42 GRB afterglows, including large data sets for GRBs 050319, 050408, 050802, 050820A, 050922C, 060418, 080413A and 080810. We analyzed the light curves of all GRBs in the sample and derived spectral energy distributions for the sample with the best data quality, allowing us to estimate the host galaxy extinction. We transformed the afterglow light curves into an extinction-corrected z=1 system and compared their luminosities with a sample of pre-Swift afterglows. The results of a former study, which showed that GRB afterglows clustered and exhibited a bimodal distribution in luminosity space, is weakened by the larger sample. We found that the luminosity distribution of the two afterglow samples (Swift-era and pre-Swift) are very similar, and that a subsample for which we were not able to estimate the extinction, which is fainter than the main sample, can be explained by assuming a moderate amount of line-of-sight host extinction. We derived bolometric isotropic energies for all GRBs in our sample, and found only a tentative correlation between the prompt energy release and the optical afterglow luminosity at one day after the GRB in the z=1 system. A comparative study of the optical luminosities of GRB afterglows with echelle spectra (which show a high number of foreground absorbing systems) and those without reveals no indication that the former are statistically significantly more luminous. (abridged), Comment: ApJ, in press; 65 pages in journal format; 20 pages main text, 18 pages Appendix, 5 pages references, 6 tables (21 pages), 9 figures, 840 original data points; v4: Updated references and acknowledgements, corrected mistake in table 1

The origins of Gamma-ray Burst prompt emission are currently not well understood and in this context long, well-observed events are particularly important to study. We present the case of GRB 070616, analysing the exceptionally long-duration multipeaked prompt emission, and later afterglow, captured by all the instruments on-board Swift and by Suzaku WAM. The high energy light curve remained generally flat for several hundred seconds before going into a steep decline. Spectral evolution from hard to soft is clearly taking place throughout the prompt emission, beginning at 285 s after the trigger and extending to 1200 s. We track the movement of the spectral peak energy, whilst observing a softening of the low energy spectral slope. The steep decline in flux may be caused by a combination of this strong spectral evolution and the curvature effect. We investigate origins for the spectral evolution, ruling out a superposition of two power laws and considering instead an additional component dominant during the late prompt emission. We also discuss origins for the early optical emission and the physics of the afterglow. The case of GRB 070616 clearly demonstrates that both broadband coverage and good time resolution are crucial to pin down the origins of the complex prompt emission in GRBs., Comment: 13 pages, 11 figures (2 in colour), MNRAS accepted

We validate the performance and accuracy of the current SEGUE (Sloan Extension for Galactic Understanding and Exploration) Stellar Parameter Pipeline (SSPP), which determines stellar atmospheric parameters (effective temperature, surface gravity, and metallicity) by comparing derived overall metallicities and radial velocities from selected likely members of three globular clusters (M 13, M 15, and M 2) and two open clusters (NGC 2420 and M 67) to the literature values. Spectroscopic and photometric data obtained during the course of the original Sloan Digital Sky Survey (SDSS-I) and its first extension (SDSS-II/SEGUE) are used to determine stellar radial velocities and atmospheric parameter estimates for stars in these clusters. Based on the scatter in the metallicities derived for the members of each cluster, we quantify the typical uncertainty of the SSPP values, sigma([Fe/H]) = 0.13 dex for stars in the range of 4500 K < Teff < 7500 K and 2.0 < log g < 5.0, at least over the metallicity interval spanned by the clusters studied (-2.3 < [Fe/H] < 0). The surface gravities and effective temperatures derived by the SSPP are also compared with those estimated from the comparison of the color-magnitude diagrams with stellar evolution models; we find satisfactory agreement. At present, the SSPP underestimates [Fe/H] for near-solar-metallicity stars, represented by members of M 67 in this study, by about 0.3 dex., Comment: 56 pages, 8 Tables, 15 figures, submitted to the Astronomical Journal

We describe the development and implementation of the SEGUE (Sloan Extension for Galactic Exploration and Understanding) Stellar Parameter Pipeline (SSPP). The SSPP derives, using multiple techniques, radial velocities and the fundamental stellar atmospheric parameters (effective temperature, surface gravity, and metallicity) for AFGK-type stars, based on medium-resolution spectroscopy and $ugriz$ photometry obtained during the course of the original Sloan Digital Sky Survey (SDSS-I) and its Galactic extension (SDSS-II/SEGUE). The SSPP also provides spectral classification for a much wider range of stars, including stars with temperatures outside of the window where atmospheric parameters can be estimated with the current approaches. This is Paper I in a series of papers on the SSPP; it provides an overview of the SSPP, and initial tests of its performance using multiple data sets. Random and systematic errors are critically examined for the current version of the SSPP, which has been used for the sixth public data release of the SDSS (DR-6)., Comment: 64 pages, 8 tables, 12 figures, submitted to the Astronomical Journal

We formulate some simple conditions under which a Markov chain may be approximated by the solution to a differential equation, with quantifiable error probabilities. The role of a choice of coordinate functions for the Markov chain is emphasised. The general theory is illustrated in three examples: the classical stochastic epidemic, a population process model with fast and slow variables, and core-finding algorithms for large random hypergraphs., Comment: Published in at http://dx.doi.org/10.1214/07-PS121 the Probability Surveys (http://www.i-journals.org/ps/) by the Institute of Mathematical Statistics (http://www.imstat.org)

The halo of the Milky Way provides unique elemental abundance and kinematic information on the first objects to form in the Universe, which can be used to tightly constrain models of galaxy formation and evolution. Although the halo was once considered a single component, evidence for its dichotomy has slowly emerged in recent years from inspection of small samples of halo objects. Here we show that the halo is indeed clearly divisible into two broadly overlapping structural components -- an inner and an outer halo -- that exhibit different spatial density profiles, stellar orbits and stellar metallicities (abundances of elements heavier than helium). The inner halo has a modest net prograde rotation, whereas the outer halo exhibits a net retrograde rotation and a peak metallicity one-third that of the inner halo. These properties indicate that the individual halo components probably formed in fundamentally different ways, through successive dissipational (inner) and dissipationless (outer) mergers and tidal disruption of proto-Galactic clumps., Comment: Two stand-alone files in manuscript, concatenated together. The first is for the main paper, the second for supplementary information. The version is consistent with the version published in Nature

Swift triggered on a precursor to the main burst of GRB 061121 (z=1.314), allowing observations to be made from the optical to gamma-ray bands. Many other telescopes, including Konus-Wind, XMM-Newton, ROTSE and the Faulkes Telescope North, also observed the burst. The gamma-ray, X-ray and UV/optical emission all showed a peak ~75s after the trigger, although the optical and X-ray afterglow components also appear early on - before, or during, the main peak. Spectral evolution was seen throughout the burst, with the prompt emission showing a clear positive correlation between brightness and hardness. The Spectral Energy Distribution (SED) of the prompt emission, stretching from 1eV up to 1MeV, is very flat, with a peak in the flux density at ~1keV. The optical-to-X-ray spectra at this time are better fitted by a broken, rather than single, power-law, similar to previous results for X-ray flares. The SED shows spectral hardening as the afterglow evolves with time. This behaviour might be a symptom of self-Comptonisation, although circumstellar densities similar to those found in the cores of molecular clouds would be required. The afterglow also decays too slowly to be accounted for by the standard models. Although the precursor and main emission show different spectral lags, both are consistent with the lag-luminosity correlation for long bursts. GRB 061121 is the instantaneously brightest long burst yet detected by Swift. Using a combination of Swift and Konus-Wind data, we estimate an isotropic energy of 2.8x10^53 erg over 1keV - 10MeV in the GRB rest frame. A probable jet break is detected at ~2x10^5s, leading to an estimate of ~10^51 erg for the beaming-corrected gamma-ray energy., Comment: 17 pages, 11 figures (1 in colour), accepted for publication in ApJ

We present techniques for the estimation of stellar atmospheric parameters (Teff,logg,[Fe/H]) for stars from the SDSS/SEGUE survey. The atmospheric parameters are derived from the observed medium-resolution (R=2000) stellar spectra using non-linear regression models trained either on (1) pre-classified observed data or (2) synthetic stellar spectra. In the first case we use our models to automate and generalize parametrization produced by a preliminary version of the SDSS/SEGUE Spectroscopic Parameter Pipeline (SSPP). In the second case we directly model the mapping between synthetic spectra (derived from Kurucz model atmospheres) and the atmospheric parameters, independently of any intermediate estimates. After training, we apply our models to various samples of SDSS spectra to derive atmospheric parameters, and compare our results with those obtained previously by the SSPP for the same samples. We obtain consistency between the two approaches, with RMS deviations of 150K in Teff, 0.35dex in logg, and 0.22dex in [Fe/H]. The models are applied to pre-processed spectra, either via Principal Components Analysis or a Wavelength Range Selection method, which employs a subset of the full 3850-9000A spectral range. This is both for computational reasons, and because it delivers higher accuracy. From an analysis of cluster candidates with available SDSS spectroscopy (M15, M13, M2, and NGC2420), we find evidence for small systematic offsets in Teff and/or logg for the atmospheric parameter estimates from the model trained on real data with the SSPP. Thus, this model turns out to derive more precise, but less accurate, atmospheric parameters than the model trained on synthetic data., Comment: 17 pages, 13 figures, accepted for publication in A&A

We numerically investigate abundance properties of the Galactic globular clusters (GCs) by adopting a new ``external pollution'' scenario. In this framework, GCs are assumed to originate in forming low-mass dwarfs embedded in dark matter subhalos at very high redshifts (z) and thus be chemically influenced by field AGB stars of the dwarfs during early GC formation processes. In this external pollution scenario, the ratio of the total mass of infalling gas to that of AGB ejecta during GC formation in a dwarf (s) and the time scale of gas infall (sigma_I) are the most important key parameters that can determine abundance properties of GCs. We mainly investigate the abundance inhomogeneity among light elements (e.g., C, N, O, Na, and Al) of stars in GCs by using the latest stellar yield models of metal-poor AGB stars with and without third dredge-up. Our principal results for the models with no third dredge-up, which are more consistent with observations, are as follows. Both [N/Fe] and [C/Fe] can be diverse among stars within a GC owing to chemical pollution from field AGB stars. [N/Fe] distributions in some GCs can clearly show bimodality whereas [C/Fe] is monomodal in most models. [N/Fe] distributions depend on s such that models with smaller $s$ (i.e., larger mass fraction of AGB ejecta used for GC formation) show the [N/Fe] bimodality more clearly. N-rich, C-poor stars in GCs also have higher He abundances owing to pollution from massive AGB stars with He-rich ejecta. The number fraction of He-rich stars (Y >0.30) is higher for the models with smaller s and shorter sigma_I for 3 < s <24 and 10^5 < sigma_I < 10^7 yr.Our models are in strong disagreement with the observed O-Na anticorrelation., Comment: 18 pages, 17 figures, accepted by MNRAS

GRB061007 is the brightest gamma-ray burst (GRB) to be detected by Swift and is accompanied by an exceptionally luminous afterglow that had a V-band magnitude <11.1 at 80s after the prompt emission. From the start of the Swift observations the afterglow decayed as a power law with a slope of \alpha_X=1.66+/-0.01 in the X-ray and \alpha_{opt}=1.64+/-0.01 in the UV/optical, up to the point that it was no longer detected above background in the optical or X-ray bands. The brightness of this GRB and the similarity in the decay rate of the X-ray, optical and gamma-ray emission from 100s after the trigger distinguish this burst from others and present a challenge to the fireball model. The lack of a cooling or jet break in the afterglow up to \~10^5s constrains any model that can produce the large luminosity observed in GRB061007, which we found to require either an excessively large kinetic energy or highly collimated outflow. Analysis of the multi-wavelength spectral and high-resolution temporal data taken with Swift suggest an early time jet-break to be a more plausible scenario. This must have occurred within 80s of the prompt emission, which places an upper limit on the jet opening angle of \theta_j=0.8deg. Such a highly collimated outflow resolves the energy budget problem presented in a spherical emission model, reducing the isotropic equivalent energy of this burst to E_{\gamma}^{corr}=10^{50} ergs; consistent with other GRBs., Comment: Accepted by MNRAS, 11 pages, 3 figures

Gamma ray bursts (GRBs) are known to come in two duration classes, separated at ~2 s. Long bursts originate from star forming regions in galaxies, have accompanying supernovae (SNe) when near enough to observe and are likely caused by massive-star collapsars. Recent observations show that short bursts originate in regions within their host galaxies with lower star formation rates consistent with binary neutron star (NS) or NS - black hole (BH) mergers. Moreover, although their hosts are predominantly nearby galaxies, no SNe have been so far associated with short GRBs. We report here on the bright, nearby GRB 060614 that does not fit in either class. Its ~102 s duration groups it with long GRBs, while its temporal lag and peak luminosity fall entirely within the short GRB subclass. Moreover, very deep optical observations exclude an accompanying supernova, similar to short GRBs. This combination of a long duration event without accompanying SN poses a challenge to both a collapsar and merging NS interpretation and opens the door on a new GRB classification scheme that straddles both long and short bursts., Comment: 13 pages, 2 figures, accepted in Nature