The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

GeV gamma ray is an important observation target of DArk Matter Particle Explorer (DAMPE) for indirect dark matter searching and high energy astrophysics. We present in this work a set of accurate instrument response functions of DAMPE (DmpIRFs) including the effective area, point-spread function and energy dispersion that are crucial for the gamma-ray data analysis based on the high statistics simulation data. A dedicated software named DmpST is developed to facilitate the scientific analyses of DAMPE gamma-ray data. Considering the limited number of photons and the angular resolution of DAMPE, the maximum likelihood method is adopted in the DmpST to better disentangle different source components. The basic mathematics and the framework regarding this software are also introduced in this paper.

The three major large-scale, diffuse $\gamma$-ray structures of the Milky Way are the Galactic disk, a bulge-like GeV excess towards the Galactic center, and the Fermi bubble. Whether such structures can also be present in other normal galaxies remains an open question. M31, as the nearest massive normal galaxy, holds promise for spatially-resolving the $\gamma$-ray emission. Based on more than 8 years of Fermi-LAT observations, we use (1) disk, (2) bulge, and (3) disk-plus-bulge templates to model the spatial distribution of the $\gamma$-ray emission from M31. Among these, the disk-plus-bulge template delivers the best-fit, in which the bulge component has a TS value 25.7 and a photon-index of $2.57\pm0.17$, providing strong evidence for a centrally-concentrated gamma-ray emission from M31, that is analogous to the Galactic center excess. The total 0.2-300 GeV gamma-ray luminosity from this bulge component is $(1.16\pm0.14)\times10^{38} erg s^{-1}$, which would require $\sim1.5\times10^5$ millisecond pulsars, if they were the dominant source. We also search for a Fermi bubble-like structure in M31 using the full dataset (pass8), but no significant evidence is found. In addition, a likelihood analysis using only photons with the most accurate reconstructed direction (i.e., PSF3-only data) reveals a 4.8 $\sigma$ point-like source located at $\sim$ 10kpc to the northwest of the M31 disk, with a luminosity of $(0.97\pm0.27)\times10^{38} erg s^{-1}$ and a photon-index of $2.31\pm0.18$. Lacking of a counterpart on the southeast side of the disk, the relation between this point-like source and a bubble-like structure remains elusive.

A tentative excess in the electron spectrum at 1.4 TeV was recently reported by the DArk Matter Particle Explorer (DAMPE). A non-astrophysical scenario in which dark matter particles annihilate or decay in a local clump has been invoked to explain the excess. If $e^\pm$ annihilation channels in the final states are mediated by left-handed leptons as a component in the $SU(2)_L$ doublet, neutrinos with similar energies should have been simultaneously produced. We demonstrate that generic dark matter models can be decisively tested by the existing IceCube data. In case of a non-detection, such models would be excluded at the $5\sigma$ level by the five-year data for a point-like source and by the ten-year data for an extended source of dark matter particles with left-handed leptons.

The tentative 1.4 TeV excess in the $e^+e^-$ spectrum measured by The DArk Matter Particle Explorer (DAMPE) motivates the possible existence of one or more local dark matter concentrated regions. In particular, Ultra-compact Micro Halos (UCMHs) seeded by large density perturbations in the early universe, allocated within ~0.3 kpc from the solar system, could provide the potential source of electrons and positrons produced from dark matter annihilation, enough to explain the DAMPE signal. Here we consider a UCMH with density profile assuming radial in-fall and explore the preferred halo parameters to explain the 1.4 TeV "DAMPE excess". We find that typical parameter space of UCMHs can easily explain the "DAMPE excess" with usual thermal-averaged annihilation cross section of WIMP. The fraction of dark matter stored in such UCMHs in the Galactic-scale halo can be reduced to as small as $O(10^{-5})$, well within the current cosmological and astrophysical constraints.

The DArk Matter Particle Explorer (DAMPE), also known as Wukong in China, launched on December 17, 2015, is a new high energy cosmic ray and {\gamma}-ray satellite-borne observatory in space. One of the main scientific goals of DAMPE is to observe GeV-TeV high energy {\gamma}-rays with accurate energy, angular, and time resolution, to indirectly search for dark matter particles and for the study of high energy astrophysics. Due to the comparatively higher fluxes of charged cosmic rays with respect to {\gamma}-rays, it is challenging to identify {\gamma}-rays with sufficiently high efficiency minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify {\gamma}-rays in DAMPE data based on Monte Carlo simulations, using the powerful electromagnetic/hadronic shower discrimination provided by the calorimeter and the veto detection of charged particles provided by the plastic scintillation detector. Monte Carlo simulations show that after this selection the number of electrons and protons that contaminate the selected {\gamma}-ray events at $\sim10$ GeV amounts to less than 1% of the selected sample. Finally, we use flight data to verify the effectiveness of the method by highlighting known {\gamma}-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

The nature of the $\gamma$-ray emission from the \emph{Fermi} bubbles is unknown. Both hadronic and leptonic models have been formulated to explain the peculiar $\gamma$-ray signal observed by the Fermi-LAT between 0.1-500~GeV. If this emission continues above $\sim$30~TeV, hadronic models of the \emph{Fermi} bubbles would provide a significant contribution to the high-energy neutrino flux detected by the IceCube observatory. Even in models where leptonic $\gamma$-rays produce the \emph{Fermi} bubbles flux at GeV energies, a hadronic component may be observable at very high energies. The combination of IceCube and HAWC measurements have the ability to distinguish these scenarios through a comparison of the neutrino and $\gamma$-ray fluxes at a similar energy scale. We examine the most recent four-year dataset produced by the IceCube collaboration and find no evidence for neutrino emission originating from the \emph{Fermi} bubbles. In particular, we find that previously suggested excesses are consistent with the diffuse astrophysical background with a p-value of 0.22 (0.05 in an extreme scenario that all the IceCube events that overlap with the bubbles come from them). Moreover, we show that existing and upcoming HAWC observations provide independent constraints on any neutrino emission from the \emph{Fermi} bubbles, due to the close correlation between the $\gamma$-ray and neutrino fluxes in hadronic interactions. The combination of these results disfavors a significant contribution from the \emph{Fermi} bubbles to the IceCube neutrino flux.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

The Ngari (Ali) prefecture of Tibet, one of the highest areas in the world, has recently emerged as a promising site for future astronomical observation. Here we use 31 years of reanalysis data from the Climate Forecast System Reanalysis (CFSR) to examine the astroclimatology of Ngari, using the recently-erected Ali Observatory at Shiquanhe (5~047~m above mean sea level) as the representative site. We find the percentage of photometric night, median atmospheric seeing and median precipitable water vapor (PWV) of the Shiquanhe site to be $57\%$, $0.8"$ and 2.5~mm, comparable some of the world's best astronomical observatories. Additional calculation supports the Shiquanhe region as one of the better sites for astronomical observations over the Tibetan Plateau. Based on the studies taken at comparable environment at Atacama, extraordinary observing condition may be possible at the few vehicle-accessible 6~000~m heights in the Shiquanhe region. Such possibility should be thoroughly investigated in future.

BICEP1 is a millimeter-wavelength telescope designed specifically to measure the inflationary B-mode polarization of the Cosmic Microwave Background (CMB) at degree angular scales. We present results from an analysis of the data acquired during three seasons of observations at the South Pole (2006 to 2008). This work extends the two-year result published in Chiang et al. (2010), with additional data from the third season and relaxed detector-selection criteria. This analysis also introduces a more comprehensive estimation of band-power window functions, improved likelihood estimation methods and a new technique for deprojecting monopole temperature-to-polarization leakage which reduces this class of systematic uncertainty to a negligible level. We present maps of temperature, E- and B-mode polarization, and their associated angular power spectra. The improvement in the map noise level and polarization spectra error bars are consistent with the 52% increase in integration time relative to Chiang et al. (2010). We confirm both self-consistency of the polarization data and consistency with the two-year results. We measure the angular power spectra at 21 <= l <= 335 and find that the EE spectrum is consistent with Lambda Cold Dark Matter (LCDM) cosmology, with the first acoustic peak of the EE spectrum now detected at 15sigma. The BB spectrum remains consistent with zero. From B-modes only, we constrain the tensor-to-scalar ratio to r = 0.03+0.27-0.23, or r < 0.70 at 95% confidence level.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

We describe the instrument concept of a high angular resolution telescope dedicated to the sub-GeV (from $\gtrsim$10 MeV to $\gtrsim$1 GeV) gamma-ray photon detection. This mission, named PANGU (PAir-productioN Gamma-ray Unit), has been suggested as a candidate for the joint small mission between the European Space Agency (ESA) and the Chinese Academy of Science (CAS). A wide range of topics of both astronomy and fundamental physics can be attacked with PANGU, covering Galactic and extragalactic cosmic-ray physics, extreme physics of a variety of extended (e.g. supernova remnants, galaxies, galaxy clusters) and compact (e.g. black holes, pulsars, gamma-ray bursts) objects, solar and terrestrial gamma-ray phenomena, and searching for dark matter decay and/or annihilation signature etc. The unprecedented point spread function can be achieved with a pair-production telescope with a large number of thin active tracking layers to precisely reconstruct the pair-produced electron and positron tracks. Scintillating fibers or thin silicon micro-strip detectors are suitable technology for such a tracker. The energy measurement is achieved by measuring the momentum of the electrons and positrons through a magnetic field. The innovated spectrometer approach provides superior photon pointing resolution, and is particular suitable in the sub-GeV range. The level of tracking precision makes it possible to measure the polarization of gamma rays, which would open up a new frontier in gamma-ray astronomy. The frequent full-sky survey at sub-GeV with PANGU's large field of view and significantly improved point spread function would provide crucial information to GeV-TeV astrophysics for current/future missions including Fermi, DAMPE, HERD, and CTA, and other multi-wavelength telescopes.

A pair of giant gamma-ray bubbles have been revealed by the {\it Fermi} LAT. In this paper we investigate their formation mechanism. Observations have indicated that the activity of the supermassive black hole located at the Galactic center, Sgr A*, was much stronger than the present time. Specifically, one possibility is that while Sgr A* was also in the hot accretion regime, the accretion rate should be $10^3-10^4$ times higher during the past $\sim 10^7$ yr. On the other hand, recent MHD numerical simulations of hot accretion flows have unambiguously shown the existence of strong winds and obtained their properties. Based on these knowledge, by performing three-dimensional hydrodynamical simulations, we show in this paper that the Fermi bubbles could be inflated by winds launched from the ``past' hot accretion flow in Sgr A*. In our model, the active phase of Sgr A* is required to last for about 10 million years and it was quenched no more than 0.2 million years ago. The Central Molecular Zone (CMZ) is included and it collimates the wind orientation towards the Galactic poles. Viscosity suppresses the Rayleigh-Taylor and Kelvin-Helmholtz instabilities and results in the smoothness of the bubble edge. The main observational features of the bubbles can be well explained. Specifically, the {\it ROSAT} X-ray features are interpreted by the shocked interstellar medium and the interaction region between winds and CMZ gas. The thermal pressure and temperature obtained in our model are in good consistency with the recent {\it Suzaku} observations.

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

Cosmic Microwave Background (CMB) polarimeters aspire to measure the faint $B$-mode signature predicted to arise from inflationary gravitational waves. They also have the potential to constrain cosmic birefringence which would produce non-zero expectation values for the CMB's $TB$ and $EB$ spectra. However, instrumental systematic effects can also cause these $TB$ and $EB$ correlations to be non-zero. In particular, an overall miscalibration of the polarization orientation of the detectors produces $TB$ and $EB$ spectra which are degenerate with isotropic cosmological birefringence, while also introducing a small but predictable bias on the $BB$ spectrum. The \bicep three-year spectra, which use our standard calibration of detector polarization angles from a dielectric sheet, are consistent with a polarization rotation of $\alpha = -2.77^\circ \pm 0.86^\circ \text{(statistical)} \pm 1.3^\circ \text{(systematic)}$. We revise the estimate of systematic error on the polarization rotation angle from the two-year analysis by comparing multiple calibration methods. We investigate the polarization rotation for the \bicep 100 GHz and 150 GHz bands separately to investigate theoretical models that produce frequency-dependent cosmic birefringence. We find no evidence in the data supporting either these models or Faraday rotation of the CMB polarization by the Milky Way galaxy's magnetic field. If we assume that there is no cosmic birefringence, we can use the $TB$ and $EB$ spectra to calibrate detector polarization orientations, thus reducing bias of the cosmological $B$-mode spectrum from leaked $E$-modes due to possible polarization orientation miscalibration. After applying this "self-calibration" process, we find that the upper limit on the tensor-to-scalar ratio decreases slightly, from $r<0.70$ to $r<0.65$ at $95\%$ confidence.

In this paper we investigate the nature and distribution of large neutral regions during the late epoch of reionization. In the "bubble model" of reionization, the mass distribution of large ionized regions ("bubbles") during the early stage of reionization is obtained by using the excursion set model, where the ionization of a region corresponds to the first up-crossing of a barrier by random trajectories. We generalize this idea, and develop a method to predict the distribution of large scale neutral regions during the late stage of reionization, taking into account the ionizing background after the percolation of HII regions. The large scale neutral regions which we call "neutral islands" are not individual galaxies or minihalos, but larger regions where fewer galaxies formed and hence ionized later, and they are identified in the excursion set model with the first down-crossings of the island barrier. Assuming that the consumption rate of ionizing background photons is proportional to the surface area of the neutral islands, we obtained the size distribution of the neutral islands. We also take the "bubbles-in-island" effect into account by considering the conditional probability of up-crossing a bubble barrier after down-crossing the island barrier. We find that this effect is very important. An additional barrier is set to avoid islands being percolated through. We find that there is a characteristic scale for the neutral islands, while the small islands are rapidly swallowed up by the ionizing background, this characteristic scale does not change much as the reionization proceeds.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.

Evidence for a spectral line in the inner Galaxy has caused a great deal of excitement over the last year, mainly because of its interpretation as a possible dark matter signal. The observation has raised important questions about statistics and suspicions about systematics, especially in photons from the Earth limb. With enough additional data, we can address these concerns. In this white paper, we summarize the current observational situation and project future sensitivities, finding that the status quo is dangerously close to leaving the issue unresolved until 2015. We advocate a change in survey strategy that more than doubles the data rate in the inner Galaxy, and is relatively non-disruptive to other survey science. This strategy will clearly separate the null hypothesis from the line signal hypothesis and provide ample limb data for systematics checks by the end of 2014. The standard survey mode may not.

Both the robust INTEGRAL 511 keV gamma-ray line and the recent tentative hint of the 135 GeV gamma-ray line from Fermi-LAT have similar signal morphologies, and may be produced from the same dark matter annihilation. Motivated by this observation, we construct a dark matter model to explain both signals and to accommodate the two required annihilation cross sections that are different by more than six orders of magnitude. In our model, to generate the low-energy positrons for INTEGRAL, dark matter particles annihilate into a complex scalar that couples to photon via a charge-radius operator. The complex scalar contains an excited state decaying into the ground state plus an off-shell photon to generate a pair of positron and electron. Two charged particles with non-degenerate masses are necessary for generating this charge-radius operator. One charged particle is predicted to be long-lived and have a mass around 3.8 TeV to explain the dark matter thermal relic abundance from its late decay. The other charged particle is predicted to have a mass below 1 TeV given the ratio of the two signal cross sections. The 14 TeV LHC will concretely test the main parameter space of this lighter charged particle.

Our recent claims of a Galactic center feature in Fermi-LAT data at approximately 130 GeV have prompted an avalanche of papers proposing explanations ranging from dark matter annihilation to exotic pulsar winds. Because of the importance of such interpretations for physics and astrophysics, a discovery will require not only additional data, but a thorough investigation of possible LAT systematics. While we do not have access to the details of each event reconstruction, we do have information about each event from the public event lists and spacecraft parameter files. These data allow us to search for suspicious trends that could indicate a spurious signal. We consider several hypotheses that might make an instrumental artifact more apparent at the Galactic center, and find them implausible. We also search for an instrumental signature in the Earth limb photons, which provide a smooth reference spectrum for null tests. We find no significant 130 GeV feature in the Earth limb sample. However, we do find a marginally significant 130 GeV feature in Earth limb photons with a limited range of detector incidence angles. This raises concerns about the 130 GeV Galactic center feature, even though we can think of no plausible model of instrumental behavior that connects the two. A modest amount of additional limb data would tell us if the limb feature is a statistical fluke. If the limb feature persists, it would raise doubts about the Pass 7 processing of E > 100 GeV events. At present we find no instrumental systematics that could plausibly explain the excess Galactic center emission at 130 GeV.

Gamma ray emission from dark matter subhalos in the Milky Way has long been sought as a sign of dark matter particle annihilation. So far, searches for gamma-ray continuum from subhalos have been unsuccessful, and line searches are difficult without prior knowledge of the line energies. Guided by recent claims of line emission at 111 GeV and 129 GeV in the Galactic center, we examine the co-added gamma-ray spectrum of unassociated point sources in the Second Fermi-LAT catalog (2FGL) using 3.9 years of LAT data. Using the SOURCE event class, we find evidence for lines at 111 GeV and 129 GeV with a local significance of 3.3 sigma based on a conservative estimate of the background at E>135 GeV. Other 2FGL sources analyzed in the same way do not show line emission at 111 GeV and 129 GeV. The line-emitting sources are mostly within 30 degrees of the Galactic plane, although this anisotropy may be a selection effect. If the double-line emission from these objects is confirmed with future data, it will provide compelling support for the hypothesis that the Galactic center line signal is indeed from dark matter annihilation.

Using 3.7 years of \Fermi-LAT data, we examine the diffuse 80-200 GeV emission in the inner Galaxy and find a resolved gamma-ray feature at $\sim 110-140$ GeV. We model the spatial distribution of this emission with a $\sim3\degree$ FWHM Gaussian, finding a best fit position $1.5\degree$ West of the Galactic Center. Even better fits are obtained for off-center Einasto and power-law profiles, which are preferred over the null (no line) hypothesis by $6.5\sigma$ ($5.0\sigma/5.4\sigma$ after trials factor correction for one/two line case) assuming an NFW density profile centered at $(\ell, b)=(-1.5\degree,0\degree)$ with a power index $\alpha=1.2$ . The energy spectrum of this structure is consistent with a single spectral line (at energy $127.0\pm 2.0$ GeV with $\chi^2=4.48$ for 4 d.o.f.). A pair of lines at $110.8\pm 4.4$ GeV and $128.8\pm 2.7$ GeV provides a marginally better fit (with $\chi^2=1.25$ for 2 d.o.f.). The total luminosity of the structure is $(3.2\pm0.6)\times 10^{35}$ erg/s, or $(1.7\pm0.4)\times 10^{36}$ photons/sec. The energies in the two-line case are compatible with a $127.3 \pm 2.7$ GeV WIMP annihilating through $\gamma \gamma$ and $\gamma Z$ (with $\chi^2=1.67$ for 3 d.o.f.). We describe a possible change to the \Fermi\ scan strategy that would accumulate S/N on spectral lines in the Galactic center 4 times as fast as the current survey strategy.

Although accretion onto supermassive black holes in other galaxies is seen to produce powerful jets in X-ray and radio, no convincing detection has ever been made of a kpc-scale jet in the Milky Way. The recently discovered pair of 10 kpc tall gamma-ray bubbles in our Galaxy may be a sign of earlier jet activity from the central black hole. In this paper, we identify a gamma-ray cocoon feature in the southern bubble, a jet-like feature along the cocoon's axis of symmetry, and another directly opposite the Galactic center in the north. Both the cocoon and jet-like feature have a hard spectrum with spectral index ~ -2 from 1 to 100 GeV, with a cocoon luminosity of (5.5 +/- 0.45) x 10^35 erg/s and luminosity of the jet-like feature of (1.8 +/- 0.35) x 10^35 erg/s at 1 to 100 GeV. If confirmed, these jets are the first resolved gamma-ray jets ever seen.

Inhomogeneous reionization gives rise to angular fluctuations in the Cosmic Microwave Background (CMB) optical depth tau(n) to the last scattering surface, correlating different spherical harmonic modes and imprinting characteristic non-Gaussianity on CMB maps. Recently the minimum variance quadratic estimator has been derived using this mode-coupling signal, and found that the optical depth fluctuations could be detected with (S/N)^2 ~ 100 in futuristic experiments like CMBPol. We first demonstrate that the non-Gaussian signal from gravitational lensing of CMB is the dominant source of contamination for reconstructing inhomogeneous reionization signals, even with 98% of its contribution removed by delensing. We then construct unbiased estimators that simultaneously reconstruct inhomogeneous reionization signals tau(n) and gravitational lensing potential phi(n). We apply our new unbiased estimators to future CMB experiment to assess the detectability of inhomogeneous reionization signals. With more physically motivated simulations of inhomogenuous reionizations that predict an order of magnitude smaller optical depth power spectrum than previous studies, we show that a CMBPol-like experiment could achieve a marginal detection of inhomogeneous reionization,(S/N)^2 ~ O(1) with this quadratic estimator to ~O(10) with the analogous maximum likelihood estimator.

Data from the Fermi-LAT reveal two large gamma-ray bubbles, extending 50 degrees above and below the Galactic center, with a width of about 40 degrees in longitude. The gamma-ray emission associated with these bubbles has a significantly harder spectrum (dN/dE ~ E^-2) than the IC emission from electrons in the Galactic disk, or the gamma-rays produced by decay of pions from proton-ISM collisions. There is no significant spatial variation in the spectrum or gamma-ray intensity within the bubbles, or between the north and south bubbles. The bubbles are spatially correlated with the hard-spectrum microwave excess known as the WMAP haze; the edges of the bubbles also line up with features in the ROSAT X-ray maps at 1.5-2 keV. We argue that these Galactic gamma-ray bubbles were most likely created by some large episode of energy injection in the Galactic center, such as past accretion events onto the central massive black hole, or a nuclear starburst in the last ~10 Myr. Dark matter annihilation/decay seems unlikely to generate all the features of the bubbles and the associated signals in WMAP and ROSAT; the bubbles must be understood in order to use measurements of the diffuse gamma-ray emission in the inner Galaxy as a probe of dark matter physics. Study of the origin and evolution of the bubbles also has the potential to improve our understanding of recent energetic events in the inner Galaxy and the high-latitude cosmic ray population.

We study the effects of instrumental systematics on the estimation of primordial non-Gaussianity using the cosmic microwave background (CMB) bispectrum from both the temperature and the polarization anisotropies. For temperature systematics we consider gain fluctuation and beam distortions. For polarization we consider effects related to known instrumental systematics: calibration, pixel rotation, differential gain, pointing, and ellipticity of the instrument beam. We consider these effects at next to leading order, which we refer to as non-linear systematic effects. We find that if the instrumental response is linearly proportional to the received CMB intensity, then only the shape of the primordial CMB bispectrum, if there is any, will be distorted. We show that the nonlinear response of the instrument can in general result in spurious non-Gaussian features on both the CMB temperature and polarization anisotropies, even if the primordial CMB is completely Gaussian. We determine the level for both the linear and non-linear systematics parameters for which they would cause no significant degradation of our ability to constrain the primordial non-Gaussianity amplitude f_{nl}. We find that the non-linear systematics are potentially bigger worry for extracting the primordial non-Gaussianity than the linear systematics. Especially because the current and near future CMB probes are optimized for CMB power-spectrum measurements which are not particularly sensitive to the non-linear instrument response. We find that if instrumental non-linearities are not controlled by dedicated calibration, the effective local non-Gaussianity can be as large as f_{nl} ~ O(10) before the corresponding non-linearities show up in the CMB dipole measurements. The higher order multipoles are even less sensitive to instrumental non-linearities.

Montage is a portable software toolkit for constructing custom, science-grade mosaics by composing multiple astronomical images. The mosaics constructed by Montage preserve the astrometry (position) and photometry (intensity) of the sources in the input images. The mosaic to be constructed is specified by the user in terms of a set of parameters, including dataset and wavelength to be used, location and size on the sky, coordinate system and projection, and spatial sampling rate. Many astronomical datasets are massive, and are stored in distributed archives that are, in most cases, remote with respect to the available computational resources. Montage can be run on both single- and multi-processor computers, including clusters and grids. Standard grid tools are used to run Montage in the case where the data or computers used to construct a mosaic are located remotely on the Internet. This paper describes the architecture, algorithms, and usage of Montage as both a software toolkit and as a grid portal. Timing results are provided to show how Montage performance scales with number of processors on a cluster computer. In addition, we compare the performance of two methods of running Montage in parallel on a grid.

Distortions in the primordial cosmic microwave background (CMB) along the line-of-sight can be modeled and described using 11 fields. These distortion fields correspond to various cosmological signals such as weak gravitational lensing of the CMB by large-scale structure, screening from patchy reionization, rotation of the plane of polarization due to magnetic fields or parity violating physics. Various instrumental systematics such as gain fluctuations, pixel rotation, differential gain, pointing, differential ellipticity are also described by the same distortion model. All these distortions produce B-mode that contaminate the primordial tensor B-modes signal. In this paper we show that apart from generating B-modes, each distortion uniquely couples different modes $(\bfl_1\ne \bfl_2)$ of the CMB anisotropies, generating $<EB>$ correlations which for the primordial CMB are zero. We describe and implement unbiased minimum variance quadratic estimators which using the off diagonal correlations in the CMB can extract the map of distortions. We perform Monte-Carlo simulations to characterize the estimators and illustrate the level of distortions that can be detected with current and future experiments. The estimators can be used to look for cosmological signals, or to check for any residual systematics in the data. As a specific example of primordial tensor B-mode diagnostics we compare the level of minimum detectable distortions using our method with maximum allowed distortion level for the B-modes detection. We show that for any experiment, the distortions will be detected at high significance using correlations before they would show up as spurious B-modes in the power spectrum.

Following Kamionkowski (2008), a quadratic estimator of the rotation of the plane of polarization of the CMB is constructed. This statistic can estimate a spatially varying rotation angle. We use this estimator to quantify the prospects of detecting such a rotation field with forthcoming experiments. For PLANCK and CMBPol we find that the estimator containing the product of the E and B components of the polarization field is the most sensitive. The variance of this EB estimator, N(L) is roughly independent of the multipole L, and is only weakly dependent on the instrumental beam. For FWHM of the beam size ~ 5'-50', and instrument noise $\Delta_p ~ 5-50 uK-arcmin, the scaling of variance N(L) can be fitted by a power law N(L)=3.3 x 10^{-7} \Delta^2_p (FWHM)^{1.3} sq-deg. For small instrumental noise \Delta_p \leq 5 uK-arcmin, the lensing B-modes become important, saturating the variance to ~10^{-6} sq-deg even for an ideal experiment. Upcoming experiments like PLANCK will be able to detect a power spectrum of the rotation angle, C^{\alpha \alpha}(L), as small as 0.01 sq-deg, while futuristic experiment like CMBPol will be able to detect rotation angle power spectrum as small as 2.5 x 10^{-5} sq-deg. We discuss the implications of such constraints, both for the various physical effects that can rotate the polarization as photons travel from the last scattering surface as well as for constraints on instrumental systematics that can also lead to a spurious rotation signal. Rotation of the CMB polarization generates B-modes which will act as contamination for the primordial B-modes detection. We discuss an application of our estimator to de-rotate the CMB to increase the sensitivity for the primordial B-modes.

In this paper, we study the effects of instrumental systematics on the reconstruction of the deflection angle power spectrum from weak lensing of Cosmic Microwave Background (CMB) temperature and polarization observations. We consider seven types of effects which are related to known instrumental systematics: calibration, rotation, pointing, spin-flip, monopole leakage, dipole leakage and quadrupole leakage. These effects can be characterized by 11 distortion fields. Each of these systematic effects can mimic the effective projected matter power spectrum and hence contaminate the lensing reconstruction. To demonstrate the effect of these instrumental systematics, we consider two types of experiments, one with a detector noise level for polarization of 9.6 uK-arcmin and FWHM of 8.0', typical of upcoming ground and balloon-based CMB experiments, and a CMBPol-like instrument with a detector noise level for polarization of 2.0 uK-arcmin and FWHM of 4.0', typical of future space-based CMB experiments. For each systematics, we consider various choices of coherence scale. Among all the 11 systematic parameters, rotation and monopole leakage place the most stringent requirements, while quadrupole leakage, pointing error, and calibration are among the least demanding. The requirements from lensing extraction are about 1-2 orders of magnitude less stringent than the requirements to measure the primordial B-modes with inflationary energy scale of 1.0*10^{16} GeV. On the other hand the requirements for lensing reconstruction are comparable or even more stringent for some systematic parameters than the requirements to detect primordial B-modes with inflationary scale E_i = 3.0*10^{16} GeV.

In this paper we analyze the constraints on the property of dark energy from cosmological observations. Together with SNe Ia Gold sample, WMAP, SDSS and 2dFGRS data, we include 69 long Gamma-Ray Bursts (GRBs) data in our study and perform global fitting using Markov Chain Monte Carlo (MCMC) technique. Dark energy perturbations are explicitly considered. We pay particular attention to the time evolution of the equation of state of dark energy parameterized as $w_{DE}=w_0+w_a(1-a)$ with $a$ the scale factor of the universe, emphasizing the complementarity of high redshift GRBs to other cosmological probes. It is found that the constraints on dark energy become stringent by taking into account high redshift GRBs, especially for $w_a$, which delineates the evolution of dark energy.

A simple speed-up cosmology model is proposed to account for the dark energy puzzle. We condense contributions from dark energy and curvature term into one effective parameter in order to reduce parameter degeneracies and to find any deviation from flat concordance $\Lambda$CDM model, by considering that the discrimination between dynamical and non-dynamical sources of cosmic acceleration as the best starting point for analyzing dark energy data sets both at present and in future. We also combine recent Type Ia Supernova (SNIa), Cosmic Microwave Background (CMB) and Baryon Oscillation (BAO) to constrain model parameter space. Degeneracies between model parameters are discussed by using both degeneracy diagram and data analysis including high redshift information from Gamma Ray Bursts (GRBs) sample. The analysis results show that our model is consistent with cosmological observations. We try to distinct the curvature effects from the specially scaling dark energy component as parameterized. We study the linear growth of large scale structure, and finally show the effective dark energy equation of state in our model and how the matter component coincidences with the dark energy numerically.

We combine recent long Gamma Ray Bursts (GRBs) sample including 52 objects out to z$=$6.3 compiled from Swift Gamma Ray Bursts by Schaefer (2006) with Type Ia Supernova (SNIa), Cosmic Microwave Background (CMB) and Baryon Oscillation (BAO) to constrain cosmologies. We study the constraints arising from GRBs alone and the complementarity of GRBs to other cosmological probes. To analyze the cosmological role of GRBs, we adopt the Hubble constant as a free parameter with a prior in the range of $H_0=72\pm8\enspace\kmsmpc$ instead of a fixed $H_0$ used in previous studies. By jointly using SNIa gold/SNLS samples and GRBs, the constraints on \Om -\Ol parameter space are dramatically improved in comparing with those from SNIa data alone. The complementarity of GRBs is mostly from the sub-sample with z$>$1.5 due to the different parameter degeneracies involved in luminosity distances at different redshifts. Including GRB data in addition to SNIa, BAO and CMB in our analysis, we find that the concordance model with \Om $\sim$0.3 and \Ol$\sim$0.7 is still well within 1$\sigma$ confidence range.