We study the disruption process of hierarchical 3-body systems with bodies of comparable mass. Such systems have long survival times that vary by orders of magnitude depending on the initial conditions. By comparing with 3-body numerical integrations, we show that the evolution and disruption of such systems can be statistically described as a simple random-walk process in the outer-orbit's energy, where the energy-exchange per pericenter passage (step-size) is calculated from the initial conditions. In our derivation of the step-size, we use previous analytic results for parabolic encounters, and average over the (Kozai-Lidov) oscillations in orbital parameters, which are faster then the energy diffusion timescale. While similar random-walk models were studied before, this work differs in two manners: (a) this is the first time that the Kozai-Lidov averaged step-size is derived from first principles and demonstrated to reproduce the statistical evolution of numerical ensembles without fitting parameters, and (b) it provides a characteristic life-time, instead of answering the binary question (stable/unstable), set by case-specific criteria.

Resolving the small length-scale of thermonuclear detonation waves (TNDWs) in supernovae is currently not possible in multidimensional full-star simulations. Additionally, multidimensional simulations usually use small, oversimplistic reaction networks and adopt an ad hoc transition criterion to nuclear statistical equilibrium (NSE). The errors due to the applied approximations are not well understood. We present here a new accurate and efficient numerical scheme that accelerates the calculations by orders of magnitudes and allows the structure of TNDWs to be resolved. The numerical scheme has two important ingredients: (1) a burning limiter that broadens the width of the TNDW while accurately preserving its internal structure, and (2) an adaptive separation of isotopes into groups that are in nuclear statistical quasi-equilibrium, which resolves the time-consuming burning calculation of reactions that are nearly balanced out. Burning is calculated in situ employing the required large networks without the use of post-processing or pre-describing the conditions behind the TNDW. In particular, the approach to and deviation from NSE are calculated self-consistently. The scheme can be easily implemented in multidimensional codes. We test our scheme against accurate solutions of the structure of TNDWs and against homogeneous expansion from NSE. We show that with resolutions that are typical for multidimensional full-star simulations, we reproduce the accurate thermodynamic trajectory (density, temperature, etc.) to an accuracy that is better than a percent for the resolved scales (where the burning limiter is not applied), while keeping the error for unresolved scales (broadened by the burning limiter) within a few percent.

Katz & Zingale (2019, KZ19) recently studied a one-dimensional test problem, intended to mimic the process of detonation ignition in head-on collisions of two carbon--oxygen (CO) white dwarfs. They do not obtain ignition of a detonation in pure CO compositions unless the temperature is artificially increased or 5% He is included. In both of these cases they obtain converged ignition only for spatial resolutions better than 0.1 km, which are beyond the capability of multidimensional simulations. This is in a contradiction with the claims of Kushnir et al. (2013, K13), that a convergence to $\sim10\%$ is achieved for a resolution of a few km. Using Eulerian and Lagrangian codes we show that a converged and resolved ignition is obtained for pure CO in this test problem without the need for He or increasing the temperature. The two codes agree to within 1% and convergence is obtained at resolutions of several km. We calculate the case that includes He and obtain a similar slow convergence, but find that it is due to a boundary numerical artifact that can (and should) be avoided. Correcting the boundary conditions allows convergence with resolution of $\sim10\,\textrm{km}$ in an agreement with the claims of K13. It is likely that the slow convergence obtained by KZ19 in this case is because of a similar boundary numerical artifact, but we are unable to verify this. KZ19 further recommended to avoid the use of the burning limiter introduced by K13. We show that their recommendation is not justified.

Near maximum brightness, the spectra of type Ia supernovae (SNe Ia) present typical absorption features of Silicon II observed at roughly 6100A and 5750A. The 2-D distribution of the pseudo-equivalent widths (pEWs) of these features is a useful tool for classifying SNe Ia spectra (Branch plot). Comparing the observed distribution of SNe on the Branch plot to results of simulated explosion models, we find that 1-D models fail to cover most of the distribution. In contrast, we find that TARDIS radiative transfer simulations of the WD head-on collision models along different lines of sight almost fully cover the distribution. We use several simplified approaches to explain this result. We perform order-of-magnitude analysis and model the opacity of the Si lines using LTE and NLTE approximations. Introducing a simple toy model of spectral feature formation, we show that the pEW is a good tracer for the extent of the absorption region in the ejecta. Using radiative transfer simulations of synthetic SNe ejecta, we reproduce the observed Branch plot distribution by varying the luminosity of the SN and the Si density profile of the ejecta. We deduce that the success of the collision model in covering the Branch plot is a result of its asymmetry, which allows for a significant range of Si density profiles along different viewing angles, uncorrelated with a range of $^{56}$Ni yields that cover the observed range of SNe Ia luminosity. We use our results to explain the shape and boundaries of the Branch plot distribution.

The diffusion approximation is often used to study supernovae light-curves around peak light, where it is applicable. By analytic arguments and numerical studies of toy models, Khatami & Kasen (2019) recently argued for a new approximate relation between peak bolometric Luminosity, $L_p$, and the time of peak since explosion, $t_p$, for transients involving homologous expansion: $L_p=2/(\beta t_p)^2\int_0 ^{\beta t_{p}} t'Q(t')dt'$, where $Q(t)$ is the heating rate of the ejecta, and $\beta$ is an order unity parameter that is calibrated from numerical calculations. Khatami & Kasen (2019) demonstrated its validity using Monte-Carlo radiation transfer simulations of ejecta with homogenous density and (for most cases considered) constant opacity. Interestingly, constant values of $\beta$ accurately reproduce the numerical calculations for different heating distributions and over a wide range of energy release times. Here we show that the diffusion and the adiabatic loss of energy in homologous expansion is equivalent to a static diffusion equation and provide an analytic solution for the case of uniform density and opacity (extending the results of Pinto & Eastman 2000). Our accurate analytical solutions reproduce and extend the results of Khatami & Kasen (2019) for this case, allowing clarification for the universality of their peak time-luminosity relation as well as new limitations to its use.

We report comprehensive multi-wavelength observations of a peculiar Type Ia-like supernova ("SN Ia-pec") ASASSN-15pz. ASASSN-15pz is a spectroscopic "twin" of SN 2009dc, a so-called "Super-Chandrasekhar-mass" SN, throughout its evolution, but it has a peak luminosity M_B,peak = -19.69 +/- 0.12 mag that is \approx 0.6 mag dimmer and comparable to the SN 1991T sub-class of SNe Ia at the luminous end of the normal width-luminosity relation. The synthesized Ni56 mass of M_Ni56 = 1.13 +/- 0.14 M_sun is also substantially less than that found for several 2009dc-like SNe. Previous well-studied 2009dc-like SNe have generally suffered from large and uncertain amounts of host-galaxy extinction, which is negligible for ASASSN-15pz. Based on the color of ASASSN-15pz, we estimate a host extinction for SN 2009dc of E(B-V)_host=0.12 mag and confirm its high luminosity (M_B, peak[2009dc] \approx -20.3 mag). The 2009dc-like SN population, which represents ~1% of SNe Ia, exhibits a range of peak luminosities, and do not fit onto the tight width-luminosity relation. Their optical light curves also show significant diversity of late-time (>~ 50 days) decline rates. The nebular-phase spectra provide powerful diagnostics to identify the 2009dc-like events as a distinct class of SNe Ia. We suggest referring to these sources using the phenomenology-based "2009dc-like SN Ia-pec" instead of "Super-Chandrasekhar SN Ia," which is based on an uncertain theoretical interpretation.

While the width-luminosity relation (WLR) among type Ia supernovae (slower is brighter) has been extensively studied, its physical basis has not been convincingly identified. In particular, the 'width' has not been quantitatively linked yet to a physical time scale. We demonstrate that there are two robust fundamental time scales that 1. can be calculated based on integral quantities of the ejecta, with little dependence on radiation transfer modeling and 2. can be inferred from observations. The first is the gamma-ray escape time $t_0$, which determines the long-term evolution of the bolometric light curve and is studied in this Paper I. The second is the recombination time of $^{56}$Fe and $^{56}$Co, which sets the long-term color evolution of the emitted light and is studied in Paper II. Here we show that the gamma-ray escape time $t_0$ can be derived with $\sim 15\%$ accuracy from bolometric observations based on first principles. When applied to a sample of supernovae, the observed values of $t_0$ span a narrow range of $30-45$ days for the wide range of observed $^{56}\rm Ni$ masses $0.1M_{\odot}\lesssim M_{^{56}\rm Ni}\lesssim 1 M_{\odot}$. %This (trivial) bolometric WLR This narrow range of the gamma-ray escape time across the range of luminosities (a trivial WLR) is consistent with central detonations and direct collisions of sub-Chandrasekhar mass white dwarfs (WDs) but not with delayed detonation models for explosions of Chandrasekhar mass WDs, which are therefore disfavored as the primary channel for the population of type Ia supernovae. Computer codes for extracting $t_0$ from observations and models and for calculating gamma-ray transfer in 1D-3D are provided.

While the width-luminosity relation (WLR) among type Ia supernovae (slower is brighter) is one of the best studied properties of this type of events, its physical basis has not been identified convincingly. The 'luminosity' is known to be related to a clear physical quantity -- the amount of $^{56}$Ni synthesized, but the 'width' has not been quantitatively linked yet to a physical time scale. We show that the recombination time of $^{56}$Fe and $^{56}$Co from doubly to singly ionized states causes the typical observed break in the color curve B-V due to a cliff in the mean opacities, and is a robust width measure of the light curve, which is insensitive to radiation transfer uncertainties. A simple photospheric model is shown to predict the recombination time to an accuracy of $\sim5$ days, allowing a quantitative understanding of the color WLR. Two physical times scales of the width luminosity relation are shown to be set by two column densities -- the total column density which sets the gamma-ray escape time $t_0$ (previous Paper I) and the $^{56}$Ni column density which sets the recombination time (this Paper II). Central detonations of sub-$\rm M_{ch}$ WDs and direct WD collision models have gamma-ray escape times and recombination times which are consistent with observations across the luminosity range of type Ia's. Delayed detonation Chandrasekhar mass models have recombination times that are broadly consistent with observations, with tension at the bright end of the luminosity range and inconsistent gamma-ray escape times at the faint end.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

We report a luminous Type II supernova, ASASSN-15nx, with a peak luminosity of M_V=-20 mag, that is between typical core-collapse supernovae and super-luminous supernovae. The post-peak optical light curves show a long, linear decline with a steep slope of 2.5 mag/100 days (i.e., an exponential decline in flux), through the end of observations at phase ~260 days. In contrast, the light curves of hydrogen rich supernovae (SNe II-P/L) always show breaks in their light curves at phase ~100 days, before settling onto Co56 radioactive decay tails with a decline rate of about 1 mag/100 days. The spectra of ASASSN-15nx do not exhibit the narrow emission-line features characteristic of Type IIn SNe, which can have a wide variety of light-curve shapes usually attributed to strong interactions with a dense circumstellar medium (CSM). ASASSN-15nx has a number of spectroscopic peculiarities, including a relatively weak and triangularly-shaped H-alpha emission profile with no absorption component. The physical origin of these peculiarities is unclear, but the long and linear post-peak light curve without a break suggests a single dominant powering mechanism. Decay of a large amount of Ni56 (M_Ni56 = 1.6 +/- 0.2 M_sun) can power the light curve of ASASSN-15nx, and the steep light-curve slope requires substantial gamma-ray escape from the ejecta, which is possible given a low-mass hydrogen envelope for the progenitor. Another possibility is strong CSM interactions powering the light curve, but the CSM needs to be sculpted to produce the unique light-curve shape and to avoid producing SN IIn-like narrow emission lines.

We present nebular-phase spectra of the Type Ia supernova (SN Ia) 2016brx, a member of the 1991bg-like subclass that lies at the faint end of the SN Ia luminosity function. Nebular spectra are available for only three other 1991bg-like SNe, and their Co line centers are all within <~ 500 km/s of each other. In contrast, the nebular Co line center of SN 2016brx is blue-shifted by >1500 km/s compared to them and by ~1200 km/s compared to the rest frame. This is a significant shift relative to the narrow nebular line velocity dispersion of <~ 2000 km/s of these SNe. The large range of nebular line shifts implies that the Ni56 in the ejecta of SN 1991bg-like events is off-center by ~1000 km/s rather than universally centrally confined as previously suggested. With the addition of SN 2016brx, the Co nebular line shapes of 1991bg-like objects appear to connect with the brighter SNe Ia that show double-peak profiles, hinting at a continuous distribution of line profiles among SNe Ia. One class of models to produce both off-center and bi-modal Ni56 distributions is collisions of white dwarfs with unequal and equal masses.

We study close approaches in hierarchical triple systems with comparable masses using full N-body simulations, motivated by a recent model for type Ia supernovae involving direct collisions of white dwarfs (WDs). For stable hierarchical systems where the inner binary components have equal masses, we show that the ability of the inner binary to achieve very close approaches, where the separation between the components of the inner binary reaches values which are orders of magnitude smaller than the semi-major axis, can be analytically predicted from initial conditions. The rate of close approaches is found to be roughly linear with the mass of the tertiary. The rate increases in systems with unequal inner binaries by a marginal factor of $\lesssim 2$ for mass ratios ${0.5<m_1/m_2<1}$ relevant for the inner white-dwarf binaries. For an average tertiary mass of $\sim 0.3 M_{\odot}$ which is representative of typical M-dwarfs, the chance for clean collisions is $\sim 1$% setting challenging constraints on the collisional model for type Ia's.

The earliest supernova (SN) emission is produced when the optical depth of the plasma lying ahead of the shock, which ejects the envelope, drops below c/v, where v is the shock velocity. This "breakout" may occur when the shock reaches the edge of the star, producing a bright X-ray/UV flash on time scales of seconds to a fraction of an hour, followed by UV/optical "cooling" emission from the expanding cooling envelope on a day time-scale. If the optical depth of circumstellar material (CSM) ejected from the progenitor star prior to the explosion is larger than c/v, the breakout will take place at larger radii, within the CSM, extending its duration to days time scale. The properties of the early, breakout and cooling, emission carry unique signatures of the structure of the progenitor star (e.g. its radius and surface composition) and of its mass-loss history. The recent progress of wide-field transient surveys enable SN detections on a day time scale, and are being used to set unique constraints on the progenitors of SNe of all types. This chapter includes a pedagogical description of SN breakout theory, and a concise overview of what we have learned from observations so far, and of advances in observational capabilities that are required in order to make further significant progress.

We find that $\sim 15-20$ per cent of A-type stars or red giants are bound with a massive companion ($M_{\rm secondary} > 1M_{\odot}$) in an intermediate wide orbit ($0.5<P<5000\mbox{ yr}$). These massive binaries are expected to form wide-orbit, double-degenerate systems (WODDs) within $\lesssim10\mbox{ Gyr}$ implying that $\sim10$ per cent of white dwarfs (WDs) are expected to be part of a WODD with a lighter WD companion. These findings are based on an analysis of previous adaptive optics observations of A-type stars and radial velocity measurements of red giants and shed light on the claimed discrepancy between the seemingly high multiplicity function of stars and the rather low number of detected double degenerates. We expect that GAIA will find $\sim 10$ new WODDs within $20\mbox{ pc}$ from the sun. These results put a stringent constraint on the collision model of type Ia supernovae in which triple stellar systems that include a WODD as the inner binary are required to be abundant.

What intrinsic properties shape the light curves of Type II supernovae (SNe)? To address this question we derive observational measures that are robust (i.e., insensitive to detailed radiative transfer) and constrain the contribution from $^{56}$Ni, as well as a combination of the envelope mass, progenitor radius, and explosion energy. By applying our methods to a sample of type II SNe from the literature we find that $^{56}$Ni contribution is often significant. In our sample its contribution to the time weighted integrated luminosity during the photospheric phase ranges between 8% and 72% with a typical value of 30%. We find that the $^{56}$Ni relative contribution is anti-correlated with the luminosity decline rate. When added to other clues, this in turn suggests that the flat plateaus often observed in type II SNe are not a generic feature of the cooling envelope emission, and that without $^{56}$Ni many of the SNe that are classified as II-P would have shown a decline rate that is steeper by up to 1 mag/100 d. Nevertheless, we find that the cooling envelope emission, and not $^{56}$Ni contribution, is the main driver behind the observed range of decline rates. Furthermore, contrary to previous suggestions, our findings indicate that fast decline rates are not driven by lower envelope masses. We therefore suggest that the difference in observed decline rates is mainly a result of different density profiles of the progenitors.

It was recently shown that the bolometric light curves of type II supernovae (SNe) allow an accurate and robust measurement of the product of the radiation energy in the ejecta, $E_r$, and the time since the explosion, $t$, at early phases ($t\lesssim 10d$) of the homologous expansion. This observable, denoted here $ET \equiv E_rt$ is constant during that time and depends only on the progenitor structure and explosion energy. We use a 1D hydrodynamic code to find $ET$ of simulated explosions of 145 red supergiant progenitors obtained using the stellar evolution code MESA, and relate this observable to the properties of the progenitor and the explosion energy. We show that $ET$ probes only the properties of the envelope (velocity, mass and initial structure), similarly to other observables that rely on the photospheric phase emission. Nevertheless, for explosions where the envelope dominates the ejected mass, $M_{env}/M_{ej} \gtrsim 0.6$, $ET$ is directly related to the explosion energy $E_{exp}$ and ejected mass $M_{ej}$ through the relation $ET \approx 0.15 E_{exp}^{1/2} R_* M_{ej}^{1/2}$, where $R_*$ is the progenitor radius, to an accuracy better than $30\%$. We also provide relations between $ET$ and the envelope properties that are accurate (to within 20\%) for all the progenitors in our sample, including those that lost most of their envelope. We show that when the envelope velocity can be reasonably measured by line shifts in observed spectra, the envelope is directly constrained from the bolometric light curve (independent of $E_{exp}$). We use that to compare observations of 11 SNe with measured $ET$ and envelope velocity to our sample of numerical progenitors. This comparison suggests that many SNe progenitors have radii that are $\lesssim 500~R_\odot$. In the framework of our simulations this indicates, most likely, a rather high value of the mixing length parameter.

The double-averaging (DA) approximation is widely employed as the standard technique in studying the secular evolution of the hierarchical three-body system. We show that effects stemmed from the short-timescale oscillations ignored by DA can accumulate over long timescales and lead to significant errors in the long-term evolution of the Lidov-Kozai cycles. In particular, the conditions for having an orbital flip, where the inner orbit switches between prograde and retrograde with respect to the outer orbit and the associated extremely high eccentricities during the switch, can be modified significantly. The failure of DA can arise for a relatively strong perturber where the mass of the tertiary is considerable compared to the total mass of the inner binary. This issue can be relevant for astrophysical systems such as stellar triples, planets in stellar binaries, stellar-mass binaries orbiting massive black holes and moons of the planets perturbed by the Sun. We derive analytical equations for the short-term oscillations of the inner orbit to the leading order for all inclinations, eccentricities and mass ratios. Under the test particle approximation, we derive the "corrected double-averaging" (CDA) equations by incorporating the effects of short-term oscillations into the DA. By comparing to N-body integrations, we show that the CDA equations successfully correct most of the errors of the long-term evolution under the DA approximation for a large range of initial conditions. We provide an implementation of CDA that can be directly added to codes employing DA equations.

We present observations of four rapidly rising (t_{rise}~10d) transients with peak luminosities between those of supernovae (SNe) and superluminous SNe (M_{peak}~-20) - one discovered and followed by the Palomar Transient Factory (PTF) and three by the Supernova Legacy Survey (SNLS). The light curves resemble those of SN 2011kl, recently shown to be associated with an ultra-long-duration gamma ray burst (GRB), though no GRB was seen to accompany our SNe. The rapid rise to a luminous peak places these events in a unique part of SN phase space, challenging standard SN emission mechanisms. Spectra of the PTF event formally classify it as a Type II SN due to broad Halpha emission, but an unusual absorption feature, which can be interpreted as either high velocity Halpha (though deeper than in previously known cases) or Si II (as seen in Type Ia SNe), is also observed. We find that existing models of white dwarf detonations, CSM interaction, shock breakout in a wind (or steeper CSM) and magnetar spindown can not readily explain the observations. We consider the possibility that a "Type 1.5 SN" scenario could be the origin of our events. More detailed models for these kinds of transients and more constraining observations of future such events should help better determine their nature.

We demonstrate that $\sim10\,\textrm{s}$ after the core-collapse of a massive star, a thermonuclear explosion of the outer shells is possible for some (tuned) initial density and composition profiles, assuming that the neutrinos failed to explode the star. The explosion may lead to a successful supernova, as first suggested by Burbidge et al. We perform a series of one-dimensional (1D) calculations of collapsing massive stars with simplified initial density profiles (similar to the results of stellar evolution calculations) and various compositions (not similar to 1D stellar evolution calculations). We assume that the neutrinos escaped with a negligible effect on the outer layers, which inevitably collapse. As the shells collapse, they compress and heat up adiabatically, enhancing the rate of thermonuclear burning. In some cases, where significant shells of mixed helium and oxygen are present with pre-collapsed burning times of $\lesssim100\,\textrm{s}$ ($\approx10$ times the free-fall time), a thermonuclear detonation wave is ignited, which unbinds the outer layers of the star, leading to a supernova. The energy released is small, $\lesssim10^{50}\,\textrm{erg}$, and negligible amounts of synthesized material (including $^{56}$Ni) are ejected, implying that these 1D simulations are unlikely to represent typical core-collapse supernovae. However, they do serve as a proof of concept that the core-collapse-induced thermonuclear explosions are possible, and more realistic two-dimensional and three-dimensional simulations are within current computational capabilities.

We discover clear doubly-peaked line profiles in 3 out of ~20 type Ia supernovae (SNe Ia) with high-quality nebular-phase spectra. The profiles are consistently present in three well-separated Co/Fe emission features. The two peaks are respectively blue-shifted and red-shifted relative to the host galaxies and are separated by ~5000 km/s. The doubly-peaked profiles directly reflect a bi-modal velocity distribution of the radioactive Ni56 in the ejecta that powers the emission of these SNe. Due to their random orientations, only a fraction of SNe with intrinsically bi-modal velocity distributions will appear as doubly-peaked spectra. Therefore SNe with intrinsic bi-modality are likely common, especially among the SNe in the low-luminosity part on the Philips relation (\Delta m15(B) >~ 1.3; ~40% of all SNe Ia). Such bi-modality is naturally expected from direct collisions of white dwarfs (WDs) due to the detonation of both WDs and is demonstrated in a 3D 0.64 M_Sun-0.64 M_Sun WD collision simulation. In the future, with a large sample of nebular spectra and a comprehensive set of numerical simulations, the collision model can be unambiguously tested as the primary channel for type Ia SNe, and the distribution of nebular line profiles will either be a smoking gun or rule it out.

We report the serendipitous discovery of a disk-eclipse system OGLE-LMC-ECL-11893. The eclipse occurs with a period of 468 days, a duration of about 15 days and a deep (up to \Delta I ~1.5), peculiar and asymmetric profile. A possible origin of such an eclipse profile involves a circumstellar disk. The presence of the disk is confirmed by the H-alpha line profile from the follow-up spectroscopic observations, and the star is identified as Be/Ae type. Unlike the previously known disk-eclipse candidates (Epsilon Aurigae, EE Cephei, OGLE-LMC-ECL-17782, KH 15D), the eclipses of OGLE-LMC-ECL-11893 retain the same shape throughout the span of ~17 years (13 orbital periods), indicating no measurable orbital precession of the disk.

The X-ray synchrotron emission of each of the young supernova-remnants (SNRs) SN1006, Kepler, Tycho, RCW86 and Cas A, is roughly given by $\nu L_{\nu}\sim 10^{45}\rm{erg}/t$, where $t$ is the remnant's age. The electrons emitting the X-ray emission cool fast, implying that the X-ray emission is calorimetric and equal to half of the cosmic ray (CR) electron acceleration efficiency (per logarithmic interval of particle energies, at multi TeV energies). Assuming Sedov-Taylor expansion, the resulting CR electron yield per SNR is estimated to be $E^2dN_e/dE\approx 6\nu L_{\nu}t \sim 10^{46}\rm erg$. This is about two orders of magnitudes below the required amount for explaining the observed electron CRs at $E\sim 10\rm GeV$. Possible resolutions are 1. a soft acceleration spectrum allowing much more energy at $E\sim 10\rm GeV$ compared to $E\sim 10\rm TeV$, 2. an increased acceleration efficiency at later phases of the SNR evolution (unlikely), or 3. SNRs are not the source of CR electrons.

We examine the luminosity function of white dwarfs (WDs) in the local ``complete'' WD sample ($d<20$ pc) of Holberg et. al. 2008. We find that the fraction of bright and young WDs is anomalously high among the WDs detected in multiple systems with main sequence (MS) companions compared to that of the single WDs and theoretical expectations. This indicates a significant observation bias against finding relatively faint WDs in multiple systems. At the bright end ($M_V<11.5$), the amount of WDs with MS companions is approximately equal to that of single white dwarfs, indicating that $\gtrsim 50\%$ of WDs have MS companions, consistent with the high multiplicity fraction of early-type MS stars. If true, a significant fraction of WDs in multiple systems within 20 pc may have not been detected yet, and the number density of WDs in the solar neighborhood and elsewhere may be up to twice as much as presently believed.

Upcoming direct-imaging experiments may detect a new class of long-period, highly luminous, tidally powered extrasolar gas giants. Even though they are hosted by ~ Gyr-"old" main-sequence stars, they can be as "hot" as young Jupiters at ~100 Myr, the prime targets of direct-imaging surveys. They are on years-long orbits and presently migrating to "feed" the "hot Jupiters." They are expected from "high-e" migration mechanisms, in which Jupiters are excited to highly eccentric orbits and then shrink semi-major axis by a factor of ~10-100 due to tidal dissipation at close periastron passages. The dissipated orbital energy is converted to heat, and if it is deposited deep enough into the atmosphere, the planet likely radiates steadily at luminosity L ~ 100-1000 L_Jup(2 x 10-7-2 x 10-6 L_Sun) during a typical ~ Gyr migration timescale. Their large orbital separations and expected high planet-to-star flux ratios in IR make them potentially accessible to high-contrast imaging instruments on 10 m class telescopes. ~10 such planets are expected to exist around FGK dwarfs within ~50 pc. Long-period radial velocity planets are viable candidates, and the highly eccentric planet HD 20782b at maximum angular separation ~0.''08 is a promising candidate. Directly imaging these tidally powered Jupiters would enable a direct test of high-e migration mechanisms. Once detected, the luminosity would provide a direct measurement of the migration rate, and together with mass (and possibly radius) estimate, they would serve as a laboratory to study planetary spectral formation and tidal physics.

We propose a stringent observational test on the formation of warm Jupiters (gas-giant planets with 10 d <~ P <~ 100 d) by high-eccentricity (high-e) migration mechanisms. Unlike hot Jupiters, the majority of observed warm Jupiters have pericenter distances too large to allow efficient tidal dissipation to induce migration. To access the close pericenter required for migration during a Kozai-Lidov cycle, they must be accompanied by a strong enough perturber to overcome the precession caused by General Relativity (GR), placing a strong upper limit on the perturber's separation. For a warm Jupiter at a ~ 0.2 AU, a Jupiter-mass (solar-mass) perturber is required to be <~ 3 AU (<~ 30 AU) and can be identified observationally. Among warm Jupiters detected by Radial Velocities (RV), >~ 50% (5 out of 9) with large eccentricities (e >~ 0.4) have known Jovian companions satisfying this necessary condition for high-e migration. In contrast, <~ 20 % (3 out of 17) of the low-e (e <~ 0.2) warm Jupiters have detected additional Jovian companions, suggesting that high-e migration with planetary perturbers may not be the dominant formation channel. Complete, long-term RV follow-ups of the warm-Jupiter population will allow a firm upper limit to be put on the fraction of these planets formed by high-e migration. Transiting warm Jupiters showing spin-orbit misalignments will be interesting to apply our test. If the misalignments are solely due to high-e migration as commonly suggested, we expect that the majority of warm Jupiters with low-e (e <~0.2) are not misaligned, in contrast with low-e hot Jupiters.

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.

We show that the recent AMS02 positron fraction measurement is consistent with a secondary origin for positrons, and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production (Katz et al 2009), obtained by neglecting radiative losses. This coincidence, which will be further tested by upcoming AMS02 data at higher energy, is a compelling indication for a secondary source. Within the secondary model the AMS02 data imply a cosmic ray propagation time in the Galaxy of < Myr and an average traversed interstellar matter density of order 1/cc, comparable to the density of the Milky Way gaseous disk, at a rigidity of 300 GV.

The early phase of the hydrodynamic evolution following collision of two stars is analyzed. Two strong shocks propagate at a constant velocity (which is a small fraction of the velocity of the approaching stars) from the contact surface toward the center of each star. The shocked region near the contact surface has a planar symmetry and a uniform pressure. The density vanishes at the (Lagrangian) surface of contact and the speed of sound diverges there. The temperature, however, reaches a finite value, since as the density vanishes, the finite pressure is radiation dominated. For Carbon-Oxygen white dwarfs collisions this temperature is too low for any appreciable nuclear burning at early times. The divergence of the speed of sound limits numerical studies of stellar collisions, as it makes convergence tests exceedingly expensive unless dedicated schemes are used. We provide a new one-dimensional Lagrangian numerical scheme to achieve this. Self-similar planar solutions are derived for zero-impact parameter collisions between two identical stars, under some simplifying assumptions. These solutions provide rough approximations that capture the main features of the flow and allow a general study as well as a detailed numerical verification test problem. The self-similar solution in the upstream frame is the planar version of previous piston problems that were studied in cylindrical and spherical symmetries. We found it timely to present a global picture of self similar piston problems. In particular, we derive new results regarding the non trivial transition to accelerating shocks at sufficiently declining densities (not relevant for collisions).

The energy output (per logarithmic interval of particle energies) of Cosmic Rays (CRs) with energies $10{\rm GeV}\lesssim\varepsilon_p\lesssim100{\rm GeV}$ is $\sim 10^{47}\rm erg$ per solar mass of star$-$formation, based on the CR production rate in the Milky Way and in starburst galaxies, implying a generation rate of $\varepsilon_p^2Q\sim 10^{45}\rm erg~Mpc^{-3}~yr^{-1}$ in the local universe. It is only $\sim 10$ times larger than the output, $\varepsilon_p^2 Q=0.5\pm0.2\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$, of Ultra High Energy CRs (UHECRs) at energies $10^{10.5}{\rm GeV}<\varepsilon_p<10^{12}\rm GeV$ (obtained assuming they are mostly protons), which in turn is comparable to the lower limit of $\varepsilon_p^2 Q\ge 0.5\times 10^{44}\rm erg~Mpc^{-3}~yr^{-1}$ of high energy CRs with $10^6{\rm GeV}\lesssim\varepsilon_p\lesssim 10^{8}\rm GeV$ implied by the saturation of the Waxman-Bahcall bound by the neutrino excess recently discovered by IceCube. These similarities are consistent with a flat production spectrum, $\varepsilon_p^2 Q\sim \text{const}$ for CRs at all observed energies. If a flat production spectrum is generated by our galaxy, the observed CR flux in the range $10^{6.5}-10^{9.5}\rm GeV$, above the "knee", is suppressed compared to lower energies due to propagation effects rather than acceleration upper limits. As suggested by Parizot and Aublin, the most exciting possibility is that cosmic rays at all energies are emitted from a single type of (unknown) sources, which can not be supernova remnants.

The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $\rho \propto x^n$ is investigated. Assuming planar geometry, local Compton equilibrium and bremsstrahlung emission as the dominant photon production mechanism, numerical solutions are obtained for the photon number density and temperature profiles as a function of time, for hydrogen-helium envelopes. The temperature solutions are determined by the breakout shock velocity $v_0$ and the pre-shock breakout density $\rho_0$, and depend weakly on the value of n. Fitting formulas for the peak surface temperature at breakout as a function of $v_0$ and $\rho_0$ are provided, with $T_{peak}\approx 9.44\exp{[12.63(v_0/c)^{1/2}]}$ eV, and the time dependence of the surface temperature is tabulated. The time integrated emitted spectrum is a robust prediction of the model, determined by $\mathcal{T}_{\rm peak}$ and $v_0$ alone and insensitive to details of light travel time or slight deviations from spherical symmetry. Adopting commonly assumed progenitor parameters, breakout luminosities of ~10^45 erg/s and ~10^44 erg/s in the 0.3-10 keV band are expected for BSG and RSG/He-WR progenitors respectively (T_{peak} is well below the band for RSGs, unless their radius is ~10^13 cm). >30 detections of SN1987A-like (BSG) breakouts are expected over the lifetime of ROSAT and XMM-Newton. An absence of such detections would imply that either the typical parameters assumed for BSG progenitors are grossly incorrect or that their envelopes are not hydrostatic. The observed spectrum and duration of XRF 080109/SN2008D are in tension with a non-relativistic breakout from a stellar surface interpretation.

An exact relation between the Ni56 mass and the bolometric light curve of a type Ia supernova can be derived as follows, using the following excellent approximations: 1. the emission is powered solely by Ni56-> Co56 ->Fe56; 2. each mass element propagates at a non-relativistic velocity which is constant in time (free coasting); and 3. the internal energy is dominated by radiation. Under these approximations, the energy E(t) carried by radiation in the ejecta satisfies: dE/dt=-E(t)/t-L(t)+Q(t), where Q(t) is the deposition of energy by the decay which is precisely known and L(t) is the bolometric luminosity. By multiplying this equation by time and integrating over time we find: E(t)*t=\int_0^t Q(t')t'dt' -\int_0^t L(t')t'dt'. At late time, t>> t_peak, the energy inside the ejecta decreases rapidly due to its escape, and thus we have \int_0^t Q(t')t'dt'=\int_0^t L(t')t'dt'. This relation is correct regardless of the opacities, density distribution or Ni56 deposition distribution in the ejecta and is very different from "Arnett's rule", L_peak ~ Q(t_peak). By comparing \int_0^t Q(t')t'dt' with \int_0^t L(t')t'dt' at t~40 day after the explosion, the mass of Ni56 can be found directly from UV, optical and infrared observations with modest corrections due to the unobserved gamma-rays and due to the small residual energy in the ejecta, E(t)*t>0.

We show that a White Dwarf-White Dwarf (WD-WD) binary with semi-major axis a=1-300 AU, which is orbited by a stellar mass outer perturber with a moderate pericenter r_{p, out} \sim 3-10 x a, has a few percent chance of experiencing a head-on collision within ~5 Gyr. Such a perturber is sufficiently distant to allow the triple system to remain intact for millions of orbits while efficiently exchanging angular momentum with the WD-WD binary. In ~ 5% of the initial orientations, the inner orbit efficiently scans the (equal energy) phase space in the region of zero angular momentum. In these systems, the binary experiences increasingly closer, stochastic, pericenter approaches r_p ~ a/2N with the increasing number (N) of orbits elapsed. Within N~10^5(a/30AU) orbits, a collision is likely to occur. This is shown by performing \simten thousand 3-body integrations and is explained by simple analytic arguments. The collisions are conservatively restricted to "clean" collisions in which all passages prior to the collision are greater than 4R_WD=4x10^9cm. In particular, within the last single orbit, the pericenter changes from r_p>4R_WD to a collision value of r_p<2R_WD. The effects of tidal deformations and General Relativistic (GR) corrections are negligible in these scenarios. The WDs approach each other with a high velocity >3000 km/s and the collision is likely to detonate the WDs leading to a type Ia SNe. If a significant fraction of WDs reside in such triples, the rate of such collisions is as high as the SNe Ia rate, and it is possible that some or all type Ia SNe occur in this way. Such SNe have a unique gravitational wave signature, which will allow a decisive identification in the future.

We show that the constant time lag prescription for tidal dissipation follows directly from the equations of motion of a tidally-forced viscous fluid body, given some basic assumptions. They are (i) dissipation results from a viscous force that is proportional to the velocity of the tidal flow (ii) tidal forcing and dissipation are weak and non-resonant (iii) the equilibrium structure of the forced body is spherically-symmetric. The lag time is an intrinsic property of the tidally-forced body and is independent of the orbital configuration.

A significant fraction of the hot Jupiters with final circularized orbital periods of less than 5 days are thought to form through the channel of high-eccentricity migration. Tidal dissipation at successive periastron passages removes orbital energy of the planet, which has the potential for changes in semi-major axis of a factor of ten to a thousand. In the equilibrium tide approximation we show that, in order for high-eccentricity migration to take place, the relative level of tidal dissipation in Jupiter analogues must be at least 10 times higher than the upper-limit attributed to the Jupiter-Io interaction. While this is not a severe problem for high-e migration, it contradicts the results of several previous calculations. We show that these calculations of high-e migration inadvertently over-estimated the strength of tidal dissipation by three to four orders of magnitude. These discrepancies were obscured by the use of various parameters, such as lag time \tau, tidal quality factor Q and viscous time t_V. We provide the values of these parameters required for the Jupiter-Io interaction, tidal circularization and high-e migration. Implications for tidal theory as well as models of the inflated radii of hot Jupiters are discussed. Though the tidal Q is not, in general, well-defined, we derive a formula for it during high-eccentricity migration where Q is approximately constant throughout evolution.

We report the discovery and detailed monitoring of X-ray emission associated with the Type IIb SN 2011dh using data from the Swift and Chandra satellites, placing it among the best studied X-ray supernovae to date. We further present millimeter and radio data obtained with the SMA, CARMA, and EVLA during the first three weeks after explosion. Combining these observations with early optical photometry, we show that the panchromatic dataset is well-described by non-thermal synchrotron emission (radio/mm) with inverse Compton scattering (X-ray) of a thermal population of optical photons. In this scenario, the shock partition fractions deviate from equipartition by a factor, (e_e/e_B) ~ 30. We derive the properties of the shockwave and the circumstellar environment and find a shock velocity, v~0.1c, and a progenitor mass loss rate of ~6e-5 M_sun/yr. These properties are consistent with the sub-class of Type IIb SNe characterized by compact progenitors (Type cIIb) and dissimilar from those with extended progenitors (Type eIIb). Furthermore, we consider the early optical emission in the context of a cooling envelope model to estimate a progenitor radius of ~1e+11 cm, in line with the expectations for a Type cIIb SN. Together, these diagnostics are difficult to reconcile with the extended radius of the putative yellow supergiant progenitor star identified in archival HST observations, unless the stellar density profile is unusual. Finally, we searched for the high energy shock breakout pulse using X-ray and gamma-ray observations obtained during the purported explosion date range. Based on the compact radius of the progenitor, we estimate that the breakout pulse was detectable with current instruments but likely missed due to their limited temporal/spatial coverage. [Abridged]

With 16-month Kepler data, 14 long-period (40 d - 265 d) eclipsing binaries on highly eccentric orbits (minimum e between 0.5 and 0.85) are recognized from their closely separated primary and secondary eclipses (\Delta t_I,II = 3 d - 10 d). These systems confirm the existence of a previously hinted binary population situated near a constant angular momentum track at P(1-e^2)^(3/2) ~ 15 d, close to the tidal circularization period P_circ. They may be presently migrating due to tidal dissipation and form a steady-state stream (~1% of stars) feeding the close-binary population (few percent of stars). If so, future Kepler data releases will reveal a growing number (dozens) of systems at longer periods, following dN/dlgP \propto P^(1/3) with increasing eccentricities reaching e -> 0.98 for P -> 1000d. Radial-velocity follow up of long-period eclipsing binaries with no secondary eclipses could offer a significantly larger sample. Orders of magnitude more (hundreds) may reveal their presence from periodic "eccentricity pulses", such as tidal ellipsoidal variations, near pericenter passages. Several new few-day-long eccentricity-pulse candidates with long period (P = 25 d - 80 d) are reported.

An important class of formation theories for hot Jupiters involves the excitation of extreme orbital eccentricity (e=0.99 or even larger) followed by tidal dissipation at periastron passage that eventually circularizes the planetary orbit at a period less than 10 days. In a steady state, this mechanism requires the existence of a significant population of super-eccentric (e>0.9) migrating Jupiters with long orbital periods and periastron distances of only a few stellar radii. For these super-eccentric planets, the periastron is fixed due to conservation of orbital angular momentum and the energy dissipated per orbit is constant, implying that the rate of change in semi-major axis a is \dot a \propto a^0.5 and consequently the number distribution satisfies dN/dlog a\propto a^0.5. If this formation process produces most hot Jupiters, Kepler should detect several super-eccentric migrating progenitors of hot Jupiters, allowing for a test of high-eccentricity migration scenarios.

Simple expressions for the radius and magnetic field of a system emitting Synchrotron-self-absorbed radio and Inverse Compton X-rays are derived from first principles which involve observable quantities only. These expressions are useful for analyzing observations of Supernova blastwaves interacting with dense circumstellar material at early times.

Kozai-Lidov oscillations of Jupiter-mass planets, excited by comparable planetary or brown dwarf mass perturbers were recently shown in numerical experiments to be slowly modulated and to exhibit striking features, including extremely high eccentricities and the generation of retrograde orbits with respect to the perturber. Here we solve this problem analytically for the case of a test particle orbiting a host star and perturbed by a distant companion whose orbit is eccentric and highly inclined. We give analytic expressions for the conditions that produce retrograde orbits and high eccentricities. This mechanism likely operates in various systems thought to involve Kozai-Lidov oscillations such as tight binaries, mergers of compact objects, irregular moons of planets and many others. In particular, it could be responsible for exciting eccentricities and inclinations of exo-planetary orbits and be important for understanding the spin-orbit (mis)alignment of hot Jupiters.

We show that a collisionless shock necessarily forms during the shock breakout of a supernova (SN) surrounded by an optically thick wind. An intense non-thermal flash of <~ MeV gamma rays, hard X-rays and multi-TeV neutrinos is produced simultaneously with and following the soft X-ray breakout emission, carrying similar or larger energy than the soft emission. The non-thermal flash is detectable by current X-ray telescopes and may be detectable out to 10's of Mpc by km-scale neutrino telescopes.

The Kozai mechanism for exponentially exciting eccentricity of a Keplerian orbit by a distant perturber is extended to a general perturbing potential. In particular, the case of an axisymmetric potential is solved analytically. The analysis is applied to orbits around an oblate central object with a distant perturber. If the equatorial plane of the central object is aligned with the orbit of the distant perturber (axisymmetric potential), a single instability zone, in which eccentricity grows exponentially, is found between two critical inclinations; if misaligned (non-axisymmetric potential), a rich set of critical inclinations separating stable and unstable zones is obtained (Vashkoviak 1974). The analysis is also applied to a general quadratic potential. Similarly, for non-axisymmetric cases, multiple stability and instability zones are obtained. Here eccentricity can reach very high values in the instability zones even when the potential's deviation from axisymmetry is small.

Exact bolometric light curves of supernova shock breakouts are derived based on the universal, non relativistic, planar breakout solutions (Sapir et al. 2011), assuming spherical symmetry, constant Thomson scattering opacity, \kappa, and angular intensity corresponding to the steady state planar limit. These approximations are accurate for progenitors with a scale height much smaller than the radius. The light curves are insensitive to the density profile and are determined by the progenitor radius R, and the breakout velocity and density, v_0 and \rho_0 respectively, and \kappa. The total breakout energy, E_BO, and the maximal ejecta velocity, v_max, are shown to be E_BO=8.0\pi R^2\kappa^-1cv_0 and v_max=2.0v_0 respectively, to an accuracy of about 10%. The calculated light curves are valid up to the time of transition to spherical expansion, t_sph\approx R/4v_0. Approximate analytic expressions for the light curves are provided for breakouts in which the shock crossing time at breakout, t_0=c/\kappa\rho_0v_0^2, is << R/c (valid for R<10^14 cm). Modifications of the flux angular intensity distribution and differences in shock arrival times to the surface, \Delta t_asym, due to moderately asymmetric explosions, affect the early light curve but do not affect v_max and E_BO. For 4v_0<<c, valid for large (RSG) progenitors, L\propto t^{-4/3} at max(\Delta t_asym,R/c)< t<t_sph and R may be accurately estimated from R\approx 2*10^13 (L/10^43 erg s^-1)^{2/5}(t/1 hr)^{8/15}.

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a power-law density profile, \rho\propto x^n, is numerically solved in the approximation of diffusion with constant opacity. For an appropriate choice of time, length and energy scales, determined by the breakout opacity, velocity and density, the solution is universal, i.e. depends only on the density power law index n. The resulting luminosity depends weakly on the value of n. An approximate analytic solution, based on the self-similar hydrodynamic solutions and on the steady RMS solutions, is constructed and shown to agree with the numerical solutions as long as the shock is far from the surface, \tau>> c/v_{sh}. Approximate analytic expressions, calibrated based on the exact solutions, are provided, that describe the escaping luminosity as a function of time. These results can be used to calculate the bolometric properties of the bursts of radiation produced during supernova (SN) shock breakouts. For completeness, we also use the exact breakout solutions to provide an analytic approximation for the maximum surface temperature for fast (v_{sh}>~0.1) non-thermal breakouts, and show that it is few times smaller than inferred based on steady-state RMS solutions.

We derive expressions for the time integrated spectrum of Cosmic Rays (CRs) that are accelerated in a decelerating relativistic shock wave and escape ahead of the shock. It is assumed that at any given time the CRs have a power law form, carry a constant fraction of the energy E_tot of the shocked plasma, and escape continuously at the maximal energy attainable. The spectrum of escaping particles is highly sensitive to the instantaneous spectral index due to the fact that the minimal energy, E_min ~ \Gamma^2 m_pc^2 where \Gamma is the shock Lorentz factor, changes with time. In particular, the escaping spectrum may be considerably harder than the canonical N(E)\propto E^-2 spectrum. For a shock expanding into a plasma of density n, a spectral break is expected at the maximal energy attainable at the transition to non relativistic velocities, E ~ 10^19 (\epsilon_B/0.1)(n/1 cm^-3)^(1/6)(E_tot/10^51 erg)^(1/3) eV where \epsilon_B is the fraction of the energy flux carried by the magnetic field. If ultra-high energy CRs are generated in decelerating relativistic blast waves arising from the explosion of stellar mass objects, their generation spectrum may therefore be different than the canonical N(E)\propto E^-2.

We present a simple analytic model for the structure of non-relativistic and relativistic radiation mediated shocks. At shock velocities \beta_s\equiv v_s/c\gtrsim 0.1, the shock transition region is far from thermal equilibrium, since the transition crossing time is too short for the production of a black-body photon density (by Bremsstrahlung emission). In this region, electrons and photons (and positrons) are in Compton (pair) equilibrium at temperatures T_s significantly exceeding the far downstream temperature, T_s\gg T_d\approx 2(\varepsilon n_u \hbar^3c^3)^{1/4}. T_s\gtrsim 10 keV is reached at shock velocities \beta_s\approx 0.2. At higher velocities, \beta_s\gtrsim0.6, the plasma is dominated in the transition region by e^\pm pairs and 60 keV\lesssim T_s \lesssim 200 keV. We argue that the spectrum emitted during the breaking out of supernova shocks from the stellar envelopes (or the surrounding winds) of Blue Super Giants and Wolf-Rayet stars, which reach \beta_s>0.1 for reasonable stellar parameters, may include a hard component with photon energies reaching tens or even hundreds of keV. This may account for the X-ray outburst associated with SN2008D, and possibly for other SN-associated outbursts with spectra not extending beyond few 100 keV (e.g. XRF060218/SN2006aj).

The growth rate of long wavelength kinetic instabilities arising due to the interaction of a collimated beam of relativistic particles and a cold unmagnetized plasma are calculated in the ultra relativistic limit. For sufficiently culminated beams, all long wave-length modes are shown to be Weibel-unstable, and a simple analytic expression for their growth rate is derived. For large transverse velocity spreads, these modes become stable. An analytic condition for stability is given. These analytic results, which generalize earlier ones given in the literature, are shown to be in agreement with numerical solutions of the dispersion equation and with the results of novel PIC simulations in which the electro-magnetic fields are restricted to a given k-mode. The results may describe the interaction of energetic cosmic rays, propagating into the far upstream of a relativistic collisionless shock, with a cold unmagnetized upstream. The long wavelength modes considered may be efficient in deflecting particles and could be important for diffusive shock acceleration. It is shown that while these modes grow in relativistic shocks propagating into electron-positron pair plasmas, they are damped in relativistic shocks propagating into electron-proton plasmas with moderate Lorenz factors \Gamma_{sh}\lesssim 100. If these modes dominate the deflection of energetic cosmic rays in electron-positron shocks, it is argued that particle acceleration is suppressed at shock frame energies that are larger than the downstream thermal energy by a factor greater than the shock Lorentz factor.

The structure of relativistic radiation mediated shocks (RRMS) propagating into a cold electron-proton plasma is calculated and analyzed. A qualitative discussion of the physics of relativistic and non relativistic shocks, including order of magnitude estimates for the relevant temperature and length scales, is presented. Detailed numerical solutions are derived for shock Lorentz factors $\Gamma_u$ in the range $6\le\Gamma_u\le30$, using a novel iteration technique solving the hydrodynamics and radiation transport equations (the protons, electrons and positrons are argued to be coupled by collective plasma processes and are treated as a fluid). The shock transition (deceleration) region, where the Lorentz factor $ \Gamma $ drops from $ \Gamma_u $ to $ \sim 1 $, is characterized by high plasma temperatures $ T\sim \Gamma m_ec^2 $ and highly anisotropic radiation, with characteristic shock-frame energy of upstream and downstream going photons of a few~$\times\, m_ec^2$ and $\sim \Gamma^2 m_ec^2$, respectively.Photon scattering is dominated by e$^\pm$ pairs, with pair to proton density ratio reaching $\approx10^2\Gamma_u$. The width of the deceleration region, in terms of Thomson optical depths for upstream going photons, is large, $\Delta\tau\sim\Gamma_u^2$ ($\Delta\tau\sim1$ neglecting the contribution of pairs) due to Klein Nishina suppression of the scattering cross section. A high energy photon component, narrowly beamed in the downstream direction, with a nearly flat power-law like spectrum, $\nu I_\nu\propto\nu^0$, and an energy cutoff at $ \sim \Gamma_u^2 m_ec^2 $ carries a fair fraction of the energy flux at the end of the deceleration region. An approximate analytic model of RRMS, reproducing the main features of the numerical results, is provided.

We present a critical analysis of the observational constraints on, and of the theoretical modeling of, aspects of cosmic ray (CR) generation and propagation in the Galaxy, which are relevant for the interpretation of recent positron and anti-proton measurements. We give simple, analytic, model independent expressions for the secondary pbar flux, and an upper limit for the secondary e+ flux, obtained by neglecting e+ radiative losses, e+/(e+ + e-)<0.2\pm0.1 up to ~300 GeV. These expressions are completely determined by the rigidity dependent grammage, which is measured from stable CR secondaries up to ~150 GeV/nuc, and by nuclear cross sections measured in the laboratory. pbar and e+ measurements, available up to ~100 GeV, are consistent with these estimates, implying that there is no need for new, non-secondary, pbar or e+ sources. The radiative loss suppression factor f_{s,e+} of the e+ flux depends on the e+ propagation in the Galaxy, which is not understood theoretically. A rough, model independent estimate of f_{s,e+} 1/3 can be obtained at a single energy, E\sim20 GeV, from unstable secondary decay and is found to be consistent with e+ measurements, including the positron fraction measured by PAMELA. We show that specific detailed models, that agree with compositional CR data, agree with our simple expressions for the e+ and pbar flux, and that the claims that the positron fraction measured by PAMELA requires new primary e+ sources are based on assumptions, that are not supported by observations. If PAMELA results are correct, they suggest that f_{s,e+} is slightly increasing with energy, which provides an interesting constraint on CR propagation models. We argue that measurements of the e+ to pbar ratio are more useful for challenging secondary production models than the positron fraction.

We argue that the observed correlation between the radio luminosity and the X-ray luminosity in radio emitting galaxy clusters implies that the radio emission is due to secondary electrons that are produced by p-p interactions and lose their energy by emitting synchrotron radiation in a strong magnetic field, B>(8\pi a T_{CMB}^4)^{1/2}\simeq 3\muG. We construct a simple model that naturally explains the correlation, and show that the observations provide stringent constraints on cluster magnetic fields and cosmic rays (CRs): Within the cores of clusters, the ratio beta_{core} between the CR energy (per logarithmic particle energy interval) and the thermal energy is beta_{core}\sim 2*10^{-4}; The source of these CRs is most likely the cluster accretion shock, which is inferred to deposit in CRs ~ 0.1 of the thermal energy it generates; The diffusion time of 100 GeV CRs over scales \gtrsim100 kpc is not short compared to the Hubble time; Cluster magnetic fields are enhanced by mergers to \gtrsim 1 % of equipartition, and decay (to <1 muG) on 1 Gyr time scales. The inferred value of beta_{core} implies that high energy gamma-ray emission from secondaries at cluster cores will be difficult to detect with existing and planned instruments.

We derive simple analytic expressions for the flux and spectrum of ultra-high energy cosmic-rays (UHECRs) predicted in models where the CRs are protons produced by extra-Galactic sources. For a power-law scaling of the CR production rate with redshift and energy, d\dot{n} /dE\propto E^-\alpha (1+z)^m, our results are accurate at high energy, E>10^18.7 eV, to better than 15%, providing a simple and straightforward method for inferring d\dot{n}/dE from the observed flux at E. We show that current measurements of the UHECR spectrum, including the latest Auger data, imply E^2d\dot{n}/dE(z=0)=(0.45\pm0.15)(\alpha-1) 10^44 erg Mpc^-3 yr^-1 at E<10^19.5 eV with \alpha roughly confined to 2\lesseq\alpha<2.7. The uncertainty is dominated by the systematic and statistic errors in the experimental determination of individual CR event energy, (\Delta E/E)_{sys} (\Delta E/E)_{stat} ~20%. At lower energy, d\dot{n}/dE is uncertain due to the unknown Galactic contribution. Simple models in which \alpha\simeq 2 and the transition from Galactic to extra-Galactic sources takes place at the "ankle", E ~10^19 eV, are consistent with the data. Models in which the transition occurs at lower energies require a high degree of fine tuning and a steep spectrum, \alpha\simeq 2.7, which is disfavored by the data. We point out that in the absence of accurate composition measurements, the (all particle) energy spectrum alone cannot be used to infer the detailed spectral shapes of the Galactic and extra-Galactic contributions.

We study relativistic unmagnetized collisionless shocks using unprecedentedly large particle-in-cell simulations of two-dimensional pair plasma. High energy particles accelerated by the shock are found to drive magnetic field evolution on a timescale >10^4 plasma times. Progressively stronger magnetic fields are generated on larger scales in a growing region around the shock. Shock-generated magnetic fields and accelerated particles carry >1% and >10% of the downstream energy flux, respectively. Our results suggest limits on the magnetization of relativistic astrophysical flows.

For more than fifty years, it has been believed that cosmic ray (CR) nuclei are accelerated to high energies in the rapidly expanding shockwaves created by powerful supernova explosions. Yet observational proof of this conjecture is still lacking. Recently, Uchiyama and collaborators reported the detection of small-scale X-ray flares in one such supernova remnant, dubbed 'RX J1713-3946' (a.k.a. G347.3-0.5), which also emits very energetic, TeV (10^12 eV) range, gamma-rays. They contend that the variability of these X-ray 'hotspots' implies that the magnetic field in the remnant is about a hundred times larger than normally assumed; and this, they say, means that the detected TeV range photons were produced in energetic nuclear interactions, providing 'a strong argument for acceleration of protons and nuclei to energies of 1 PeV (10^15 eV) and beyond in young supernova remnants.' We point out here that the existing multiwavelength data on this object certainly do not support such conclusions. Though intriguing, the small-scale X-ray flares are not the long sought-after 'smoking gun' of nucleonic CR acceleration in SNRs.

We present a simple analytic model for the various contributions to the non-thermal emission from shell type SNRs, and show that this model's results reproduce well the results of previous detailed calculations. We show that the \geq 1 TeV gamma ray emission from the shell type SNRs RX J1713.7-3946 and RX J0852.0-4622 is dominated by inverse-Compton scattering of CMB photons (and possibly infra-red ambient photons) by accelerated electrons. Pion decay (due to proton-proton collisions) is shown to account for only a small fraction, \lesssim10^-2, of the observed flux, as assuming a larger fractional contribution would imply nonthermal radio and X-ray synchrotron emission and thermal X-ray Bremsstrahlung emission that far exceed the observed radio and X-ray fluxes. Models where pion decay dominates the \geq 1 TeV flux avoid the implied excessive synchrotron emission (but not the implied excessive thermal X-ray Bremsstrahlung emission) by assuming an extremely low efficiency of electron acceleration, K_ep \lesssim 10^-4 (K_ep is the ratio of the number of accelerated electrons and the number of accelerated protons at a given energy). We argue that observations of SNRs in nearby galaxies imply a lower limit of K_ep \gtrsim 10^-3, and thus rule out K_ep values \lesssim 10^-4 (assuming that SNRs share a common typical value of K_ep). It is suggested that SNRs with strong thermal X-ray emission, rather than strong non-thermal X-ray emission, are more suitable candidates for searches of gamma rays and neutrinos resulting from proton-proton collisions. In particular, it is shown that the neutrino flux from the SNRs above is probably too low to be detected by current and planned neutrino observatories (Abridged).

We derive constraints that must be satisfied by the sources of ~10^{15} to ~10^{18} eV cosmic rays, under the assumption that the sources are Galactic. We show that while these constraints are not satisfied by ordinary supernovae (SNe), which are believed to be the sources of <10^{15} eV cosmic rays, they may be satisfied by the recently discovered class of trans-relativistic supernovae (TRSNe), which were observed in association with gamma-ray bursts. We define TRSNe as SNe that deposit a large fraction, f_R>10^{-2}, of their kinetic energy in mildly relativistic, \gamma\beta>1, ejecta. The high velocity ejecta enable particle acceleration to ~10^{18} eV, and the large value of f_R (compared to f_R~10^{-7} for ordinary SNe) ensures that if TRSNe produce the observed ~10^{18} eV cosmic ray flux, they do not overproduce the flux at lower energies. This, combined with the estimated rate and energy production of TRSNe, imply that Galactic TRSNe may be the sources of cosmic rays with energies up to ~10^{18}eV .

Observations of gamma-ray burst afterglows suggest that the correlation length of magnetic field fluctuations downstream of relativistic non-magnetized collisionless shocks grows with distance from the shock to scales much larger than the plasma skin depth. We argue that this indicates that the plasma properties are described by a self-similar solution, and derive constraints on the scaling properties of the solution. For example, we find that the scaling of the characteristic magnetic field amplitude with distance from the shock is B \propto D^{s_B} with -1<s_B<=0, that the spectrum of accelerated particles is dn/dE \propto E^{-2/(s_B+1)}, and that the scaling of the magnetic correlation function is <B_i(x)B_j(x+dx)> \propto x^{2s_B} (for x>>D). We show that the plasma may be approximated as a combination of two self-similar components: a kinetic component of energetic particles and an MHD-like component representing "thermal" particles. We argue that the latter may be considered as infinitely conducting, in which case s_B=0 and the scalings are completely determined (e.g. dn/dE \propto E^{-2} and B \propto D^0). Similar claims apply to non- relativistic shocks such as in supernova remnants, if the upstream magnetic field can be neglected. Self-similarity has important implications for any model of particle acceleration and/or field generation. For example, we show that the diffusion function in the angle \mu of momentum p in diffusive shock acceleration models must satisfy D_{\mu\mu}(p,D) = D^{-1}D'_{\mu\mu}(p/D), and that a previously suggested model for the generation of large scale magnetic fields through a hierarchical merger of current-filaments should be generalized. A numerical experiment testing our analysis is outlined (Abridged).