• We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.
• We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
• A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($\delta$, $B$ and $R$), in addition to the parameters related to the particle population.
• We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.
• We present the results of 71.6 hours of observations of the Geminga pulsar (PSR J0633+1746) with the VERITAS very-high-energy gamma-ray telescope array. Data taken with VERITAS between November 2007 and February 2013 were phase-folded using a Geminga pulsar timing solution derived from data recorded by the XMM-\emph{Newton} and \emph{Fermi}-LAT space telescopes. No significant pulsed emission above 100 GeV is observed, and we report upper limits at the 95% confidence level on the integral flux above 135 GeV (spectral analysis threshold) of 4.0$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ and 1.7$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ for the two principal peaks in the emission profile. These upper limits, placed in context with phase-resolved spectral energy distributions determined from five years of data from the \emph{Fermi}-LAT, constrain possible hardening of the Geminga pulsar emission spectra above $\sim$50 GeV.
• Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
• ### Super-Sample CMB Lensing(1401.7992)

Jan. 30, 2014 astro-ph.CO
Lensing of the CMB by modes that are larger than the size of the survey dilates intrinsic scales in the temperature and polarization fields and coherently shifts their observed power spectra with respect to the ensemble or all-sky mean. The effect can be simply encapsulated as a contribution to the power spectrum covariance matrix in accordance with the lensing trispectrum or as an additional parameter, the mean convergence in the field, for parameter estimation. It should be included for upcoming surveys that precisely measure acoustic polarization features deep into the damping tail at multipoles of $\ell \gtrsim 1500$ with less than $10\%$ of sky. Its omission may lead to seemingly conflicting values for the angular scale of the sound horizon which may then provide erroneous cosmological parameters when compared to baryon acoustic oscillation measurements.
• ### Physical properties of simulated galaxy populations at z=2 - I. Effect of metal-line cooling and feedback from star formation and AGN(1211.1021)

Aug. 7, 2013 astro-ph.CO
We use hydrodynamical simulations from the OWLS project to investigate the dependence of the physical properties of galaxy populations at redshift 2 on metal-line cooling and feedback from star formation and active galactic nuclei (AGN). We find that if the sub-grid feedback from star formation is implemented kinetically, the feedback is only efficient if the initial wind velocity exceeds a critical value. This critical velocity increases with galaxy mass and also if metal-line cooling is included. This suggests that radiative losses quench the winds if their initial velocity is too low. If the feedback is efficient, then the star formation rate is inversely proportional to the amount of energy injected per unit stellar mass formed (which is proportional to the initial mass loading for a fixed wind velocity). This can be understood if the star formation is self-regulating, i.e. if the star formation rate (and thus the gas fraction) increase until the outflow rate balances the inflow rate. Feedback from AGN is efficient at high masses, while increasing the initial wind velocity with gas pressure or halo mass allows one to generate galaxy-wide outflows at all masses. Matching the observed galaxy mass function requires efficient feedback. In particular, the predicted faint-end slope is too steep unless we resort to highly mass loaded winds for low-mass objects. Such efficient feedback from low-mass galaxies (M_* << 10^10 Msun) also reduces the discrepancy with the observed specific star formation rates, which are higher than predicted unless the feedback transitions from highly efficient to inefficient just below the observed stellar mass range.
• ### Parity violation constraints using 2006-2007 QUaD CMB polarization spectra(0811.0618)

April 2, 2009 astro-ph, hep-ph
We constrain parity-violating interactions to the surface of last scattering using spectra from the QUaD experiment's second and third seasons of observations by searching for a possible systematic rotation of the polarization directions of CMB photons. We measure the rotation angle due to such a possible "cosmological birefringence" to be 0.55 deg. +/- 0.82 deg. (random) +/- 0.5 deg. (systematic) using QUaD's 100 and 150 GHz TB and EB spectra over the multipole range 200 < l < 2000, consistent with null, and constrain Lorentz violating interactions to < 2^-43 GeV (68% confidence limit). This is the best constraint to date on electrodynamic parity violation on cosmological scales.
• In this paper we present a parameter estimation analysis of the polarization and temperature power spectra from the second and third season of observations with the QUaD experiment. QUaD has for the first time detected multiple acoustic peaks in the E-mode polarization spectrum with high significance. Although QUaD-only parameter constraints are not competitive with previous results for the standard 6-parameter LCDM cosmology, they do allow meaningful polarization-only parameter analyses for the first time. In a standard 6-parameter LCDM analysis we find the QUaD TT power spectrum to be in good agreement with previous results. However, the QUaD polarization data shows some tension with LCDM. The origin of this 1 to 2 sigma tension remains unclear, and may point to new physics, residual systematics or simple random chance. We also combine QUaD with the five-year WMAP data set and the SDSS Luminous Red Galaxies 4th data release power spectrum, and extend our analysis to constrain individual isocurvature mode fractions, constraining cold dark matter density, alpha(cdmi)<0.11 (95 % CL), neutrino density, alpha(ndi)<0.26 (95 % CL), and neutrino velocity, alpha(nvi)<0.23 (95 % CL), modes. Our analysis sets a benchmark for future polarization experiments.
• ### Measurements of cosmic-ray secondary nuclei at high energies with the first flight of the CREAM balloon-borne experiment(0808.1718)

Aug. 12, 2008 astro-ph
We present new measurements of heavy cosmic-ray nuclei at high energies per- formed during the first flight of the balloon-borne cosmic-ray experiment CREAM (Cosmic-Ray Energetics And Mass). This instrument uses multiple charge detectors and a transition radiation detector to provide the first high accuracy measurements of the relative abundances of elements from boron to oxygen up to energies around 1 TeV/n. The data agree with previous measurements at lower energies and show a relatively steep decline (~E$^-0.6$ to E$^-0.5$) at high energies. They further show the source abundance of nitrogen relative to oxygen is ~10% in the TeV/n region.
• ### On relaxation processes in collisionless mergers(astro-ph/0609612)

Nov. 30, 2006 astro-ph
We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.
• ### Type Ia Supernova Explosion: Gravitationally Confined Detonation(astro-ph/0405163)

May 8, 2004 astro-ph
We present a new mechanism for Type Ia supernova explosions in massive white dwarfs. The proposed scenario follows from relaxing the assumption of symmetry in the model and involves a detonation created in an unconfined environment. The explosion begins with an essentially central ignition of stellar material initiating a deflagration. This deflagration results in the formation of a buoyantly-driven bubble of hot material that reaches the stellar surface at supersonic speeds. The bubble breakout forms a strong pressure wave that laterally accelerates fuel-rich outer stellar layers. This material, confined by gravity to the white dwarf, races along the stellar surface and is focused at the location opposite to the point of the bubble breakout. These streams of nuclear fuel carry enough mass and energy to trigger a detonation just above the stellar surface. The flow conditions at that moment support a detonation that will incinerate the white dwarf and result in an energetic explosion. The stellar expansion following the deflagration redistributes stellar mass in a way that ensures production of intermediate mass and iron group elements consistent with observations. The ejecta will have a strongly layered structure with a mild amount of asymmetry following from the early deflagration phase. This asymmetry, combined with the amount of stellar expansion determined by details of the evolution (principally the energetics of deflagration, timing of detonation, and structure of the progenitor), can be expected to create a family of mildly diverse Type Ia supernova explosions.
• ### Lorentz invariance violation in top-down scenarios of ultrahigh energy cosmic ray creation(hep-ph/0306288)

Jan. 28, 2004 astro-ph, hep-ph
The violation of Lorentz invariance (LI) has been invoked in a number of ways to explain issues dealing with ultrahigh energy cosmic ray (UHECR) production and propagation. These treatments, however, have mostly been limited to examples in the proton-neutron system and photon-electron system. In this paper we show how a broader violation of Lorentz invariance would allow for a series of previously forbidden decays to occur, and how that could lead to UHECR primaries being heavy baryonic states or Higgs bosons.
• ### The Quiescent X-Ray Spectrum of the Neutron Star in Cen X-4 Observed with Chandra/ACIS-S(astro-ph/0012400)

Dec. 19, 2000 astro-ph
We report on spectral and intensity variability analysis from a Chandra/ACIS-S observation of the transient, type-I X-ray bursting low-mass X-ray binary Cen X-4. The quiescent X-ray spectrum during this observation is statistically identical to one observed previously with Beppo/SAX, and close, but not identical, to one observed previously with ASCA. The X-ray spectrum is best described as a pure Hydrogen atmosphere thermal spectrum plus a power-law component that dominates the spectrum above 2 keV. The best-fit radius of the neutron star is r=12.9+/-2.6 (d/1.2 kpc) km if the interstellar absorption is fixed at the value implied by the optical reddening. Allowing the interstellar absorption to be a free parameter yields r=19+45-10 (d/1.2 kpc) km (90% confidence). The thermal spectrum from the neutron star surface is inconsistent with a solar metallicity. We find a 3sigma upper-limit of root-mean-square variability <18% (0.2-2.0 keV; 0.0001-1 Hz) during the observation. On the other hand, the 0.5-10.0 keV luminosity decreased by 40+/-8% in the 4.9 years between the Asca and Chandra observations. This variability can be attributed to the power-law component. Moreover, we limit the variation in thermal temperature to <10% over these 4.9 years. The stability of the thermal temperature and emission area radius supports the interpretation that the quiescent thermal emission is due to the hot neutron star core.
• ### Anisotropy in the Microwave Sky at 90 GHz: Results From Python III(astro-ph/9606175)

June 28, 1996 astro-ph
The third year of observations with the Python microwave background experiment densely sample a $5.5^o\times 22^o$ region of sky that includes the fields measured during the first two years of observations with this instrument. The sky is sampled in two multipole bands centered at $l \approx 92$ and $l \approx 177$. These two data sets are analyzed to place limits on fluctuations in the microwave sky at 90 GHz. Interpreting the observed fluctuations as anisotropy in the cosmic microwave background we find flat band power estimates of $\delta T_l\equiv \sqrt{l(l+1)C_l/2\pi }% =54_{-12}^{+14}\mu K$ at $l=92_{-24}^{+37}$ and $\delta T_l=58_{-13}^{+15}\mu K$ at $l=177_{-58}^{+66}$. Combining the entire three year set of Python observations, we find that the angular power spectrum of fluctuations has a spectral index $m=.02_{-.16}^{+.18}$ and an amplitude $% \delta T_{l_e}=53_{-11}^{+13}\mu K$ at $l_e=145_{-74}^{+59}$ for the functional form $\delta T_l=\delta T_{l_e}(l/l_e)^m$. The stated uncertainties in the amplitudes and spectral index represent 1$\sigma$ confidence intervals in the likelihood added in quadrature with a 20\% calibration uncertainty and an estimate of the effects introduced due to imperfect overlap of the beams on the sky. The limits of $l$ are determined from the full width at half maximum of the window functions.
• ### Constraints on the Strength of Primordial Magnetic Fields from Big Bang Nucleosynthesis Revisited(astro-ph/9606163)

June 26, 1996 astro-ph, hep-ph, gr-qc, nucl-th
In this paper, we revisit in detail the effects of primordial magnetic fields on big bang nucleosynthesis (BBN) including a discussion of the magnetic field geometry and the anomalous magnetic moment. The presence of magnetic fields affects BBN by (1) increasing the weak reaction rates; (2) increasing the electron density due to changes to the electron phase space; and (3) by increasing the expansion rate of the universe, due both to the magnetic field energy density and to the modified electron energy density. Of the effects considered, the increase in the expansion rate due to the magnetic field energy is the most significant for the interests of BBN. The allowed magnetic field intensity at the end of nucleosynthesis (0.01 MeV) is about $2 \times 10^{9}$G and corresponds to an upper limit on the magnetic field energy density of about 28\% of the neutrino energy density ($\rho_B \le 0.28 \rho_\nu$).