Showing papers in "arXiv: High Energy Astrophysical Phenomena in 2013"

The Man Behind the Curtain: X-rays Drive the UV through NIR Variability in the 2013 AGN Outburst in NGC 2617

Abstract: After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ~70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such "changing look Active Galactic Nuclei (AGN)" are rare and provide us with important insights about AGN physics. Based on the Hbeta line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 +/- 1) x 10^7 M_sun. When we cross-correlate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2-3 days) to the NIR (6-9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

Abstract: This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Radiative Transfer Simulations for Neutron Star Merger Ejecta

Abstract: The merger of binary neutron stars (NSs) is among the most promising gravitational wave (GW) sources. Next-generation GW detectors are expected to detect signals from the NS merger within 200 Mpc. Detection of electromagnetic wave (EM) counterpart is crucial to understand the nature of GW sources. Among possible EM emission from the NS merger, emission powered by radioactive r-process nuclei is one of the best targets for follow-up observations. However, prediction so far does not take into account detailed r-process element abundances in the ejecta. We perform radiative transfer simulations for the NS merger ejecta including all the r-process elements from Ga to U for the first time. We show that the opacity in the NS merger ejecta is about kappa = 10 cm^2 g^{-1}, which is higher than that of Fe-rich Type Ia supernova ejecta by a factor of ~ 100. As a result, the emission is fainter and longer than previously expected. The spectra are almost featureless due to the high expansion velocity and bound-bound transitions of many different r-process elements. We demonstrate that the emission is brighter for a higher mass ratio of two NSs and a softer equation of states adopted in the merger simulations. Because of the red color of the emission, follow-up observations in red optical and near-infrared (NIR) wavelengths will be the most efficient. At 200 Mpc, expected brightness of the emission is i = 22 - 25 AB mag, z = 21 - 23 AB mag, and 21 - 24 AB mag in NIR JHK bands. Thus, observations with wide-field 4m- and 8m-class optical telescopes and wide-field NIR space telescopes are necessary. We also argue that the emission powered by radioactive energy can be detected in the afterglow of nearby short gamma-ray bursts.

Leptonic and Hadronic Modeling of Fermi-Detected Blazars

Abstract: We describe new implementations of leptonic and hadronic models for the broadband emission from relativistic jets in AGN in a temporary steady state. For the leptonic model, a temporary equilibrium between particle injection/acceleration, radiative cooling, and escape from a spherical emission region is evaluated, and the self-consistent radiative output is calculated. For the hadronic model, a temporary equilibrium between particle injection/acceleration, radiative and adiabatic cooling, and escape is evaluated for both primary electrons and protons. A new, semi-analytical method to evaluate the radiative output from cascades initiated by internal gamma-gamma pair production is presented. We use our codes to fit snap-shot spectral energy distributions of a representative set of Fermi-LAT detected blazars. We find that the leptonic model provides acceptable fits to the SEDs of almost all blazars with parameters close to equipartition between the magnetic field and the relativistic electron population. However, the hard gamma-ray spectrum of AO 0235+164, in contrast to the very steep IR-optical-UV continuum, poses a severe problem for the leptonic model. If charge neutrality in leptonic models is provided by cold protons, the kinetic energy carried by the jet should be dominated by protons. We find satisfactory representations of the snapshot SEDs of most blazars in our sample with the hadronic model presented here. However, in the case of two quasars the characteristic break at a few GeV energies can not be well modelled. All of our hadronic model fits require powers in relativistic protons in the range L_p ~ 1e47 - 1e49 erg/s.

The McGill Magnetar Catalog

Abstract: We present a catalog of the 26 currently known magnetars and magnetar candidates. We tabulate astrometric and timing data for all catalog sources, as well as their observed radiative properties, particularly the spectral parameters of the quiescent X-ray emission. We show histograms of the spatial and timing properties of the magnetars, comparing them with the known pulsar population, and we investigate and plot possible correlations between their timing, X-ray, and multiwavelength properties. We find the scale height of magnetars to be in the range 20-31 pc, assuming they are exponentially distributed. This range is smaller than that measured for OB stars, providing evidence that magnetars are born from the most massive O stars. From the same fits, we find that the Sun lies ~13-22 pc above the Galactic plane, consistent with previous measurements. We confirm previously identified correlations between quiescent X-ray luminosity L_X and magnetic field B, as well as X-ray spectral power-law index Gamma and B, and show evidence for an excluded region in a plot of L_X vs. Gamma. We also present an updated kT versus characteristic age plot, showing magnetars and high-B radio pulsars are hotter than lower-B neutron stars of similar age. Finally, we observe a striking difference between magnetars detected in the the hard X-ray and radio bands; there is a clear correlation between the hard and soft X-ray flux, whereas the radio-detected magnetars all have low soft X-ray flux suggesting, if anything, that the two bands are anti-correlated. An online version of the catalog is located at this http URL

Effect of a High Opacity on the Light Curves of Radioactively Powered Transients from Compact Object Mergers

Abstract: The coalescence of compact objects are a promising astrophysical sources of gravitational wave (GW) signals. The ejection of r-process material from such mergers may lead to a radioactively-powered electromagnetic counterpart which, if discovered, would enhance the science return of a GW detection. As very little is known about the optical properties of heavy r-process elements, previous light curve models have adopted opacities similar to those of iron group elements. Here we report that the presence of heavier elements, particularly the lanthanides, increase the ejecta opacity by several orders of magnitude. We include these higher opacities in time dependent, multi-wavelength radiative transport calculations to predict the broadband light curves of one-dimensional models over a range of parameters (ejecta masses from 0.001 to 0.1 solar masses and velocities from 0.1 to 0.3c). We find that the higher opacities lead to much longer duration light curves which can last a week or more. The emission is shifted toward the infrared bands due to strong optical line blanketing, and the colors at later times are representative of a blackbody near the recombination temperature of the lanthanides (T ~ 2500K). We further consider the case in which a second mass outflow, composed of 56Ni, is ejected from a disk wind, and show that the net result is a distinctive two component spectral energy distribution, with a bright optical peak due to 56Ni and an infrared peak due to r-process ejecta. We briefly consider the prospects for detection and identification of these transients.

Super Luminous Ic Supernovae: catching a magnetar by the tail

Abstract: We report extensive observational data for five of the lowest redshift Super-Luminous Type Ic Supernovae (SL-SNe Ic) discovered to date, namely PTF10hgi, SN2011ke, PTF11rks, SN2011kf and SN2012il. Photometric imaging of the transients at +50 to +230 days after peak combined with host galaxy subtraction reveals a luminous tail phase for four of these SL-SNe. A high resolution, optical and near infrared spectrum from xshooter provides detection of a broad He I $\lambda$10830 emission line in the spectrum (+50d) of SN2012il, revealing that at least some SL-SNe Ic are not completely helium free. At first sight, the tail luminosity decline rates that we measure are consistent with the radioactive decay of \co, and would require 1-4M of i to produce the luminosity. These i masses cannot be made consistent with the short diffusion times at peak, and indeed are insufficient to power the peak luminosity. We instead favour energy deposition by newborn magnetars as the power source for these objects. A semi-analytical diffusion model with energy input from the spin-down of a magnetar reproduces the extensive lightcurve data well. The model predictions of ejecta velocities and temperatures which are required are in reasonable agreement with those determined from our observations. We derive magnetar energies of $0.4\lesssim E$($10^{51}$erg) $\lesssim6.9$ and ejecta masses of $2.3\lesssim M_{ej}$(\M) $\lesssim 8.6$. The sample of five SL-SNe Ic presented here, combined with SN 2010gx - the best sampled SL-SNe Ic so far - point toward an explosion driven by a magnetar as a viable explanation for all SL-SNe Ic.

Simulations of Ion Acceleration at Non-relativistic Shocks. I. Acceleration Efficiency

Abstract: We use 2D and 3D hybrid (kinetic ions - fluid electrons) simulations to investigate particle acceleration and magnetic field amplification at non-relativistic astrophysical shocks. We show that diffusive shock acceleration operates for quasi-parallel configurations (i.e., when the background magnetic field is almost aligned with the shock normal) and, for large sonic and Alfv\'enic Mach numbers, produces universal power-law spectra proportional to p^(-4), where p is the particle momentum. The maximum energy of accelerated ions increases with time, and it is only limited by finite box size and run time. Acceleration is mainly efficient for parallel and quasi-parallel strong shocks, where 10-20% of the bulk kinetic energy can be converted to energetic particles, and becomes ineffective for quasi-perpendicular shocks. Also, the generation of magnetic turbulence correlates with efficient ion acceleration, and vanishes for quasi-perpendicular configurations. At very oblique shocks, ions can be accelerated via shock drift acceleration, but they only gain a factor of a few in momentum, and their maximum energy does not increase with time. These findings are consistent with the degree of polarization and the morphology of the radio and X-ray synchrotron emission observed, for instance, in the remnant of SN 1006. We also discuss the transition from thermal to non-thermal particles in the ion spectrum (supra-thermal region), and we identify two dynamical signatures peculiar of efficient particle acceleration, namely the formation of an upstream precursor and the alteration of standard shock jump conditions.

Turbulent, Extreme Multi-Zone Model for Simulating Flux and Polarization Variability in Blazars

Abstract: The author presents a model for variability of the flux and polarization of blazars in which turbulent plasma flowing at a relativistic speed down a jet crosses a standing conical shock. The shock compresses the plasma and accelerates electrons to energies up to gamma(max) > 1E4 times their rest-mass energy, with the value of gamma(max) determined by the direction of the magnetic field relative to the shock front. The turbulence is approximated in a computer code as many cells, each with a uniform magnetic field whose direction is selected randomly. The density of high-energy electrons in the plasma changes randomly with time in a manner consistent with the power spectral density of flux variations derived from observations of blazars. The variations in flux and polarization are therefore caused by continuous noise processes rather than by singular events such as explosive injection of energy at the base of the jet. Sample simulations illustrate the behavior of flux and linear polarization versus time that such a model produces. The variations in gamma-ray flux generated by the code are often, but not always, correlated with those at lower frequencies, and many of the flares are sharply peaked. The mean degree of polarization of synchrotron radiation is higher and its time-scale of variability shorter toward higher frequencies, while the polarization electric vector sometimes randomly executes apparent rotations. The slope of the spectral energy distribution exhibits sharper breaks than can arise solely from energy losses. All of these results correspond to properties observed in blazars.

Measurement of the Radius of Neutron Stars with High S/N Quiescent Low-mass X-ray Binaries in Globular Clusters

Abstract: This paper presents the measurement of the neutron star (NS) radius using the thermal spectra from quiescent low-mass X-ray binaries (qLMXBs) inside globular clusters (GCs). Recent observations of NSs have presented evidence that cold ultra dense matter -- present in the core of NSs -- is best described by "normal matter" equations of state (EoSs). Such EoSs predict that the radii of NSs, Rns, are quasi-constant (within measurement errors, of ~10%) for astrophysically relevant masses (Mns > 0.5 Msun). The present work adopts this theoretical prediction as an assumption, and uses it to constrain a single Rns value from five qLMXB targets with available high signal-to-noise X-ray spectroscopic data. Employing a Markov-Chain Monte-Carlo approach, we produce the marginalized posterior distribution for Rns, constrained to be the same value for all five NSs in the sample. An effort was made to include all quantifiable sources of uncertainty into the uncertainty of the quoted radius measurement. These include the uncertainties in the distances to the GCs, the uncertainties due to the Galactic absorption in the direction of the GCs, and the possibility of a hard power-law spectral component for count excesses at high photon energy, which are observed in some qLMXBs in the Galactic plane. Using conservative assumptions,we found that the radius, common to the five qLMXBs and constant for a wide range of masses, lies in the low range of possible NS radii, Rns=9.1(+1.3)(-1.5) km (90%-confidence). Such a value is consistent with low-res equations of state. We compare this result with previous radius measurements of NSs from various analyses of different types of systems. In addition, we compare the spectral analyses of individual qLMXBs to previous works.

X-ray reflected spectra from accretion disk models. III. A complete grid of ionized reflection calculations

Abstract: We present a new and complete library of synthetic spectra for modeling the component of emission that is reflected from an illuminated accretion disk The spectra were computed using an updated version of our code XILLVER that incorporates new routines and a richer atomic data base We offer in the form of a table model an extensive grid of reflection models that cover a wide range of parameters Each individual model is characterized by the photon index \Gamma of the illuminating radiation, the ionization parameter \xi at the surface of the disk (ie, the ratio of the X-ray flux to the gas density), and the iron abundance A_{Fe} relative to the solar value The ranges of the parameters covered are: 12 \leq \Gamma \leq 34, 1 \leq \xi \leq 10^4, and 05 \leq A_{Fe} \leq 10 These ranges capture the physical conditions typically inferred from observations of active galactic nuclei, and also stellar-mass black holes in the hard state This library is intended for use when the thermal disk flux is faint compared to the incident power-law flux The models are expected to provide an accurate description of the Fe K emission line, which is the crucial spectral feature used to measure black hole spin A total of 720 reflection spectra are provided in a single FITS file{\url{this http URL}} suitable for the analysis of X-ray observations via the atable model in XSPEC Detailed comparisons with previous reflection models illustrate the improvements incorporated in this version of XILLVER

Measuring Black Hole Spin using X-ray Reflection Spectroscopy

Abstract: I review the current status of X-ray reflection (a.k.a. broad iron line) based black hole spin measurements. This is a powerful technique that allows us to measure robust black hole spins across the mass range, from the stellar-mass black holes in X-ray binaries to the supermassive black holes in active galactic nuclei. After describing the basic assumptions of this approach, I lay out the detailed methodology focusing on "best practices" that have been found necessary to obtain robust results. Reflecting my own biases, this review is slanted towards a discussion of supermassive black hole (SMBH) spin in active galactic nuclei (AGN). Pulling together all of the available XMM-Newton and Suzaku results from the literature that satisfy objective quality control criteria, it is clear that a large fraction of SMBHs are rapidly-spinning, although there are tentative hints of a more slowly spinning population at high (M>5*10^7Msun) and low (M<2*10^6Msun) mass. I also engage in a brief review of the spins of stellar-mass black holes in X-ray binaries. In general, reflection-based and continuum-fitting based spin measures are in agreement, although there remain two objects (GROJ1655-40 and 4U1543-475) for which that is not true. I end this review by discussing the exciting frontier of relativistic reverberation, particularly the discovery of broad iron line reverberation in XMM-Newton data for the Seyfert galaxies NGC4151, NGC7314 and MCG-5-23-16. As well as confirming the basic paradigm of relativistic disk reflection, this detection of reverberation demonstrates that future large-area X-ray observatories such as LOFT will make tremendous progress in studies of strong gravity using relativistic reverberation in AGN.

Spectroscopy of The Largest Ever Gamma-ray Selected BL Lac Sample

Abstract: We report on spectroscopic observations covering most of the 475 BL Lacs in the 2nd Fermi LAT catalog of AGN. Including archival measurements (correcting several erroneous literature values) we now have spectroscopic redshifts for 44% of the BL Lacs. We establish firm lower redshift limits via intervening absorption systems and statistical lower limits via searches for host galaxies for an additional 51% of the sample leaving only 5% of the BL Lacs unconstrained. The new redshifts raise the median spectroscopic z from 0.23 to 0.33 and include redshifts as large as z=2.471. Spectroscopic redshift minima from intervening absorbers have ~ z= 0.70, showing a substantial fraction at large z and arguing against strong negative evolution. We find that detected BL Lac hosts are bright ellipticals with black hole masses M_\bullet ~ 10^{8.5-9}, substantially larger than the mean of optical AGN and LAT Flat Spectrum Radio Quasar samples. A slow increase in M_\bullet with z may be due to selection bias. We find that the power-law dominance of the optical spectrum extends to extreme values, but this does not strongly correlate with the gamma-ray properties, suggesting that strong beaming is the primary cause of the range in continuum dominance.

A new population of ultra-long duration gamma-ray bursts

Abstract: We present comprehensive multiwavelength observations of three gamma-ray bursts (GRBs) with durations of several thousand seconds. We demonstrate that these events are extragalactic transients; in particular we resolve the long-standing conundrum of the distance of GRB 101225A (the "Christmas-day burst"), finding it to have a redshift z=0.847, and showing that two apparently similar events (GRB 111209A and GRB 121027A) lie at z=0.677 and z=1.773 respectively. The systems show extremely unusual X-ray and optical lightcurves, very different from classical GRBs, with long lasting highly variable X-ray emission and optical light curves that exhibit little correlation with the behaviour seen in the X-ray. Their host galaxies are faint, compact, and highly star forming dwarf galaxies, typical of "blue compact galaxies". We propose that these bursts are the prototypes of a hitherto largely unrecognized population of ultra-long GRBs, that while observationally difficult to detect may be astrophysically relatively common. The long durations may naturally be explained by the engine driven explosions of stars of much larger radii than normally considered for GRB progenitors which are thought to have compact Wolf-Rayet progenitor stars. However, we cannot unambiguously identify supernova signatures within their light curves or spectra. We also consider the alternative possibility that they arise from the tidal disruption of stars by supermassive black holes.

Analytical Light Curve Models of Super-Luminous Supernovae: chi^2-Minimizations of Parameter Fits

Abstract: We present fits of generalized semi-analytic supernova (SN) light curve (LC) models for a variety of power inputs including Ni-56 and Co-56 radioactive decay, magnetar spin-down, and forward and reverse shock heating due to supernova ejecta-circumstellar matter (CSM) interaction. We apply our models to the observed LCs of the H-rich Super Luminous Supernovae (SLSN-II) SN 2006gy, SN 2006tf, SN 2008am, SN 2008es, CSS100217, the H-poor SLSN-I SN 2005ap, SCP06F6, SN 2007bi, SN 2010gx and SN 2010kd as well as to the interacting SN 2008iy and PTF09uj. Our goal is to determine the dominant mechanism that powers the LCs of these extraordinary events and the physical conditions involved in each case. We also present a comparison of our semi-analytical results with recent results from numerical radiation hydrodynamics calculations in the particular case of SN 2006gy in order to explore the strengths and weaknesses of our models. We find that CS shock heating produced by ejecta-CSM interaction provides a better fit to the LCs of most of the events we examine. We discuss the possibility that collision of supernova ejecta with hydrogen-deficient CSM accounts for some of the hydrogen-deficient SLSNe (SLSN-I) and may be a plausible explanation for the explosion mechanism of SN 2007bi, the pair-instability supernova (PISN) candidate. We characterize and discuss issues of parameter degeneracy.

NuSTAR discovery of a 3.76-second transient magnetar near Sagittarius A*

Abstract: We report the discovery of 3.76-s pulsations from a new burst source near Sgr A* observed by the NuSTAR Observatory. The strong signal from SGR J1745-29 presents a complex pulse profile modulated with pulsed fraction 27+/-3 % in the 3-10 keV band. Two observations spaced 9 days apart yield a spin-down rate of Pdot = (6.5+/-1.4)x10^{-12}. This implies a magnetic field B = 1.6x10^14 G, spin-down power Edot = 5x10^33 erg/s, and characteristic age P/2Pdot = 9x10^3 yr, for the rotating dipole model. However, the current Pdot may be erratic, especially during outburst. The flux and modulation remained steady during the observations and the 3-79 keV spectrum is well fitted by a combined blackbody plus power-law model with temperature kT_BB = 0.96+/-0.02 keV and photon index = 1.5+/-0.4, respectively. The neutral hydrogen column density (nH ~ 1.4x10^23 cm^{-2}) measured by NuSTAR and Swift suggests that SGR J1745-29 is located at or near the Galactic Center. The lack of an X-ray counterpart in the published Chandra survey catalog sets a quiescent 2-8 keV luminosity limit of Lx ~< 10^32 erg/s. The bursting, timing, and spectral properties indicate a transient magnetar undergoing an outburst with 2-79 keV luminosity up to 3.5x10^35 erg/s for a distance of 8 kpc. SGR J1745-29 joins a growing subclass of transient magnetars, indicating that many magnetars in quiescence remain undetected in the X-ray band or have been detected as high-B radio pulsars. The peculiar location of SGR J1745-29 has important implications for the formation and dynamics of neutron stars in the Galactic Center region.

Head-on collisions of white dwarfs in triple systems could explain type Ia supernova

Abstract: Type Ia supernovae (SNe Ia), thermonuclear explosions of carbon-oxygen white dwarfs (CO-WDs), are currently the best cosmological "standard candles", but the triggering mechanism of the explosion is unknown. It was recently shown that the rate of head-on collisions of typical field CO-WDs in triple systems may be comparable to the SNe Ia rate. Here we provide evidence supporting a scenario in which the majority of SNe Ia are the result of such head-on collisions of CO-WDs. In this case, the nuclear detonation is due to a well understood shock ignition, devoid of commonly introduced free parameters such as the deflagration velocity or transition to detonation criteria. By using two-dimensional hydrodynamical simulations with a fully resolved ignition process, we show that zero-impact-parameter collisions of typical CO-WDs with masses $0.5-1\,M_{\odot}$ result in explosions that synthesize $^{56}$Ni masses in the range of $\sim0.1-1\,M_{\odot}$, spanning the wide distribution of yields observed for the majority of SNe Ia. All collision models yield the same late-time ($>60$ days since explosion) bolometric light curve when normalized by $^{56}$Ni masses (to better than $30\%$), in agreement with observations. The calculated widths of the $^{56}$Ni-mass-weighted-line-of-sight velocity distributions are correlated with the calculated $^{56}$Ni yield, agreeing with the observed correlation. The strong correlation, shown here for the first time, between $^{56}$Ni yield and total mass of the colliding CO-WDs (insensitive to their mass ratio), is suggestive as the source for the continuous distribution of observed SN Ia features, possibly including the Philips relation.

A panchromatic view of the restless SN2009ip reveals the explosive ejection of a massive star envelope

Abstract: The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

Black hole triple dynamics: breakdown of the orbit average approximation and implications for gravitational wave detections

Abstract: Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.

Demographics of the Galaxies Hosting Short-duration Gamma-Ray Bursts

Abstract: We present observations of the afterglows and host galaxies of three short-duration gamma-ray bursts (GRBs): 100625A, 101219A and 110112A. We find that GRB 100625A occurred in a z=0.452 early-type galaxy with a stellar mass of 4.6e9 M_Sun and a stellar population age of 0.7 Gyr, and GRB 101219A originated in a star-forming galaxy at z=0.718 with a stellar mass of 1.4e9 M_Sun, a star formation rate of 16 M_Sun yr^-1, and a stellar population age of 50 Myr. We also report the discovery of the optical afterglow of GRB 110112A, which lacks a coincident host galaxy to i>26 mag and we cannot conclusively identify any field galaxy as a possible host. The bursts have inferred circumburst densities of ~1e-4-1 cm^-3, and isotropic-equivalent gamma-ray and kinetic energies of 1e50-1e51 erg. These events highlight the diversity of galaxies that host short GRBs. To quantify this diversity, we use the sample of 36 Swift short GRBs with robust associations to an environment (~1/2 of 68 short bursts detected by Swift to May 2012) and classify them as follows: late-type (50%), early-type (15%), inconclusive (20%), and host-less (lacking a coincident host galaxy to limits of >26 mag; 15%). To find likely ranges for the true late- and early-type fractions, we assign each of the host-less bursts to the late- or early-type category using probabilistic arguments, and consider the scenario that all hosts in the inconclusive category are early-type galaxies to set an upper bound on the early-type fraction. The most likely ranges for the late- and early-type fractions are ~60-80% and ~20-40%, respectively. We find no clear trend between gamma-ray duration and host type, and no change to the fractions when excluding events recently claimed as possible contaminants from the long GRB/collapsar population. Our reported demographics are consistent with a short GRB rate driven by both stellar mass and star formation.

The Origin of Galactic Cosmic Rays

Abstract: One century ago Viktor Hess carried out several balloon flights that led him to conclude that the penetrating radiation responsible for the discharge of electroscopes was of extraterrestrial origin. One century from the discovery of this phenomenon seems to be a good time to stop and think about what we have understood about Cosmic Rays. The aim of this review is to illustrate the ideas that have been and are being explored in order to account for the observable quantities related to cosmic rays and to summarize the numerous new pieces of observation that are becoming available. In fact, despite the possible impression that development in this field is somewhat slow, the rate of new discoveries in the last decade or so has been impressive, and mainly driven by beautiful pieces of observation. At the same time scientists in this field have been able to propose new, fascinating ways to investigate particle acceleration inside the sources, making use of multifrequency observations that range from the radio, to the optical, to X-rays and gamma rays. These ideas can now be confronted with data. I will mostly focus on supernova remnants as the most plausible sources of Galactic cosmic rays, and I will review the main aspects of the modern theory of diffusive particle acceleration at supernova remnant shocks, with special attention for the dynamical reaction of accelerated particles on the shock and the phenomenon of magnetic field amplification at the shock. Cosmic ray escape from the sources is discussed as a necessary step to determine the spectrum of cosmic rays at the Earth. In the end of this review I will also discuss the phenomenon of cosmic ray acceleration at shocks propagating in partially ionized media and the implications of this phenomenon in terms of width of the Balmer line emission.

Binary Black Hole Accretion From a Circumbinary Disk: Gas Dynamics Inside the Central Cavity

Abstract: We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

The Locations of Short Gamma-ray Bursts as Evidence for Compact Object Binary Progenitors

Abstract: We present a detailed investigation of Hubble Space Telescope (HST) rest-frame UV/optical observations of 22 short gamma-ray burst (GRB) host galaxies and sub-galactic environments. Utilizing the high angular resolution and depth of HST, we characterize the host galaxy morphologies, measure precise projected physical and host-normalized offsets between the bursts and host centers, and calculate the locations of the bursts with respect to their host light distributions (rest-frame UV and optical). We calculate a median short GRB projected physical offset of 4.5 kpc, about 3.5 times larger than that for long GRBs, and find that ~25% of short GRBs have offsets of >10 kpc. When compared to their host sizes, the median offset is 1.5 half-light radii (r_e), about 1.5 times larger than the values for long GRBs, core-collapse supernovae, and Type Ia supernovae. In addition, ~20% of short GRBs having offsets of >5 r_e, and only ~25% are located within 1 r_e. We further find that short GRBs severely under-represent their hosts' rest-frame optical and UV light, with ~30-45% of the bursts located in regions of their host galaxies that have no detectable stellar light, and ~55% in the regions with no UV light. Therefore, short GRBs do not occur in regions of star formation or even stellar mass. This demonstrates that the progenitor systems of short GRBs must migrate from their birth sites to their eventual explosion sites, a signature of kicks in compact object binary systems. Utilizing the full sample of offsets, we estimate natal kick velocities of ~20-140 km s^-1. These independent lines of evidence provide the strongest support to date that short GRBs result from the merger of compact object binaries (NS-NS/NS-BH).

Short GRB 130603B: Discovery of a jet break in the optical and radio afterglows, and a mysterious late-time X-ray excess

Abstract: We present radio, optical/NIR, and X-ray observations of the afterglow of the short-duration 130603B, and uncover a break in the radio and optical bands at 0.5 d after the burst, best explained as a jet break with an inferred jet opening angle of 4-8 deg. GRB 130603B is only the third short GRB with a radio afterglow detection to date, and the first time that a jet break is evident in the radio band. We model the temporal evolution of the spectral energy distribution to determine the burst explosion properties and find an isotropic-equivalent kinetic energy of (0.6-1.7) x 10^51 erg and a circumburst density of 5 x 10^-3-30 cm^-3. From the inferred opening angle of GRB 130603B, we calculate beaming-corrected energies of Egamma (0.5-2) x 10^49 erg and EK (0.1-1.6) x 10^49 erg. Along with previous measurements and lower limits we find a median short GRB opening angle of 10 deg. Using the all-sky observed rate of 10 Gpc^-3 yr^-1, this implies a true short GRB rate of 20 yr^-1 within 200 Mpc, the Advanced LIGO/VIRGO sensitivity range for neutron star binary mergers. Finally, we uncover evidence for significant excess emission in the X-ray afterglow of GRB 130603B at >1 d and conclude that the additional energy component could be due to fall-back accretion or spin-down energy from a magnetar formed following the merger.

The First Fermi-LAT Catalog of Sources Above 10 GeV

Abstract: We present a catalog of gamma-ray sources at energies above 10 GeV based on data from the Large Area Telescope (LAT) accumulated during the first 3 years of the Fermi Gamma-ray Space Telescope mission. This catalog complements the Second Fermi-LAT Catalog, which was based on 2 years of data extending down to 100 MeV and so included many sources with softer spectra below 10 GeV. The First Fermi-LAT Catalog of >10 GeV sources (1FHL) has 514 sources, and includes their locations, spectra, a measure of their variability, and associations with cataloged sources at other wavelengths. We found that 449 (87%) could be associated with known sources, of which 393 (76% of the 1FHL sources) are active galactic nuclei. We also highlight the subset of the 1FHL sources that are the best candidates for detection at energies above 50 GeV with ground-based gamma-ray observatories.

Neutron Star Masses and Radii from Quiescent Low-Mass X-ray Binaries

Abstract: We perform a systematic analysis of neutron star radius constraints from five quiescent low-mass X-ray binaries and examine how they depend on measurements of their distances and amounts of intervening absorbing material, as well as their assumed atmospheric compositions. We construct and calibrate to published results a semi-analytic model of the neutron star atmosphere which approximates these effects for the predicted masses and radii. Starting from mass and radius probability distributions established from hydrogen-atmosphere spectral fits of quiescent sources, we apply this model to compute alternate sets of probability distributions. We perform Bayesian analyses to estimate neutron star mass-radius curves and equation of state (EOS) parameters that best-fit each set of distributions, assuming the existence of a known low-density neutron star crustal EOS, a simple model for the high-density EOS, causality, and the observation that the neutron star maximum mass exceeds $2~M_\odot$. We compute the posterior probabilities for each set of distance measurements and assumptions about absorption and composition. We find that, within the context of our assumptions and our parameterized EOS models, some absorption models are disfavored. We find that neutron stars composed of hadrons are favored relative to those with exotic matter with strong phase transitions. In addition, models in which all five stars have hydrogen atmospheres are found to be weakly disfavored. Our most likely models predict neutron star radii that are consistent with current experimental results concerning the nature of the nucleon-nucleon interaction near the nuclear saturation density.

Helium-ignited violent mergers as a unified model for normal and rapidly declining Type Ia Supernovae

Abstract: The progenitors of Type Ia Supernovae (SNe Ia) are still unknown, despite significant progress during the last years in theory and observations. Violent mergers of two carbon--oxygen (CO) white dwarfs (WDs) are one candidate suggested to be responsible for at least a significant fraction of normal SNe Ia. Here, we simulate the merger of two CO WDs using a moving-mesh code that allows for the inclusion of thin helium (He) shells (0.01\,\msun) on top of the WDs, at an unprecedented numerical resolution. The accretion of He onto the primary WD leads to the formation of a detonation in its He shell. This detonation propagates around the CO WD and sends a converging shock wave into its core, known to robustly trigger a second detonation, as in the well-known double-detonation scenario for He-accreting CO WDs. However, in contrast to that scenario where a massive He shell is required to form a detonation through thermal instability, here the He detonation is ignited dynamically. Accordingly the required He-shell mass is significantly smaller, and hence its burning products are unlikely to affect the optical display of the explosion. We show that this scenario, which works for CO primary WDs with CO- as well as He-WD companions, has the potential to explain the different brightness distributions, delay times and relative rates of normal and fast declining SNe Ia. Finally, we discuss extensions to our unified merger model needed to obtain a comprehensive picture of the full observed diversity of SNe Ia.

Simulating radio emission from air showers with CoREAS

Abstract: CoREAS is a Monte Carlo code for the simulation of radio emission from extensive air showers. It implements the endpoint formalism for the calculation of electromagnetic radiation directly in CORSIKA. As such, it is parameter-free, makes no assumptions on the emission mechanism for the radio signals, and takes into account the complete complexity of the electron and positron distributions as simulated by CORSIKA. In this article, we illustrate the capabilities of CoREAS with simulations carried out in different frequency ranges from tens of MHz up to GHz frequencies, and describe in particular the emission characteristics at high frequencies due to Cherenkov effects arising from the varying refractive index of the atmosphere.

Towards the event horizon -- the supermassive black hole in the Galactic Center

Abstract: The center of our Galaxy hosts the best constrained supermassive black hole in the universe, Sagittarius A* (Sgr A*). Its mass and distance have been accurately determined from stellar orbits and proper motion studies, respectively, and its high-frequency radio, and highly variable near-infrared and X-ray emission originate from within a few Schwarzschild radii of the event horizon. The theory of general relativity (GR) predicts the appearance of a black hole shadow, which is a lensed image of the event horizon. This shadow can be resolved by very long baseline radio interferometry and test basic predictions of GR and alternatives thereof. In this paper we review our current understanding of the physical properties of Sgr A*, with a particular emphasis on the radio properties, the black hole shadow, and models for the emission and appearance of the source. We argue that the Galactic Center holds enormous potential for experimental tests of black hole accretion and theories of gravitation in their strong limits.

Progenitor models of the electromagnetic transient associated with the short GRB 130603B

Abstract: An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.