Source author record

Boaz Katz

Boaz Katz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2021arXiv

Polarization signatures of the head-on collision model for Type Ia supernovae: How much asymmetry is too much?

In a previous paper, we showed that the asymmetric ejecta produced by (zero impact parameter) head-on collisions of carbon-oxygen white dwarfs allow these progenitor models for Type Ia supernovae (SNe Ia) to cover the observed two-dimensional (2D) distribution of Si II line depths (Branch plot). In this paper, we study the polarization signature associated with the 2D asymmetric ejecta of the collision model and a double-detonation model using similar TARDIS radiative transfer simulations along different lines of sight with a spherical photosphere, combined with a new 3D Monte Carlo polarization code. We show that the polarization $Q$ can be parametrized as a product $Q=Q_{\max}Q_{\rm{x}}$ of a radial structure component $Q_{\max}$ which is insensitive to the model specifics and is shown to be universally around $Q_{\max}\sim 5\%$, and a cancellation component $Q_{\rm{x}}$ which depends on the asymmetry details. The continuum polarization is found to be low for both the collision and double-detonation models with $Q\sim 0.5\%$. However, the irregular Si distribution in the 2D head-on collision model results in Si II line polarization reaching $Q\sim 3\%$ ($Q_{\rm{x}} \lesssim 50\%$) in tension with observations (mostly $\lesssim 1.2\%$). In contrast, we show that the double-detonation model also covers the Branch plot, and yet results in low line polarization $Q\lesssim 0.7\%$ ($Q_{\rm{x}} \sim 10\%$) consistent with previous results and most SNe Ia. These results strengthen the case for asymmetric explosions as progenitors of SNe Ia, emphasizing an additional requirement for large polarization cancellations to account for the low observed line polarizations.

preprint2020arXiv

A Simple Random-Walk Model Explains the Disruption Process of Hierarchical, Eccentric 3-Body Systems

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.

preprint2020arXiv

An accurate and efficient numerical calculation of detonation waves in multidimensional supernova simulations using a burning limiter and adaptive quasi-statistical equilibrium

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.

preprint2020arXiv

An asymmetric explosion mechanism may explain the diversity of Si II line widths in Type Ia supernovae

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.

preprint2019arXiv

Comments on "Numerical Stability of Detonations in White Dwarf Simulations"

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.

preprint2016arXiv

Double-Averaging Can Fail to Characterize the Long-Term Evolution of Lidov-Kozai Cycles & Derivation of an Analytical Correction

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.

preprint2016arXiv

Rapidly Rising Transients in the Supernova - Superluminous Supernova Gap

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.

preprint2016arXiv

The importance of 56Ni in shaping the light curves of type II supernovae

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.

preprint2016arXiv

Type II supernovae progenitor and ejecta properties from the total emitted light, ET

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.

preprint2015arXiv

Failure of a neutrino-driven explosion after core-collapse may lead to a thermonuclear supernova

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.

preprint2015arXiv

Type Ia Supernovae with Bi-Modal Explosions Are Common -- Possible Smoking Gun for Direct Collisions of White-Dwarfs

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 (Δ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.

preprint2014arXiv

Directly Imaging Tidally Powered Migrating Jupiters

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.

preprint2014arXiv

Low X-ray emission challenges supernovae remnants as the source of cosmic-ray electrons

The X-ray synchrotron emission of each of the young supernova-remnants (SNRs) SN1006, Kepler, Tycho, RCW86 and Cas A, is roughly given by $νL_ν\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νL_ν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.

preprint2014arXiv

Luminosity function suggests up to 100 white dwarfs within 20 pc may be hiding in multiple systems

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.

preprint2014arXiv

OGLE-LMC-ECL-11893: The discovery of a long-period eclipsing binary with a circumstellar disk

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 Δ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.

preprint2013arXiv

AMS02 results support the secondary origin of cosmic ray positrons

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.

preprint2013arXiv

An exact integral relation between the Ni56 mass and the bolometric light curve of a type Ia supernova

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.

preprint2013arXiv

Early Hydrodynamic Evolution of a Stellar Collision

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).

preprint2013arXiv

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

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.

preprint2013arXiv

Non-relativistic radiation mediated shock breakouts: III. Spectral properties of SN shock breakout

The spectrum of radiation emitted following shock breakout from a star's surface with a power-law density profile $ρ\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 $ρ_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 $ρ_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.

preprint2013arXiv

The energy production rate density of cosmic rays in the local universe is $\sim10^{44-45}\rm erg~Mpc^{-3}~yr^{-1}$ at all particle energies

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.

preprint2013arXiv

Warm Jupiters Need Close "Friends" for High-Eccentricity Migration -- A Stringent Upper Limit on the Perturber's Separation

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.

preprint2012arXiv

Exploring a Stream of Highly-Eccentric Binaries with Kepler

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 (Δ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.

preprint2012arXiv

Physical Basis for a Constant Lag Time

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.

preprint2012arXiv

Q in Other Solar Systems

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 τ, 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.

preprint2012arXiv

The rate of WD-WD head-on collisions may be as high as the SNe Ia rate

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.

preprint2011arXiv

Exponential growth of eccentricity in secular theory

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.

preprint2011arXiv

Long-Term Cycling of Kozai-Lidov Cycles: Extreme Eccentricities and Inclinations Excited by a Distant Eccentric Perturber

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.

preprint2011arXiv

Non-relativistic radiation mediated shock breakouts: I. Exact bolometric planar breakout solutions

The problem of a non-steady planar radiation mediated shock (RMS) breaking out from a surface with a power-law density profile, ρ\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, τ>> 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.

preprint2011arXiv

Non-relativistic radiation mediated shock breakouts: II. Bolometric properties of SN shock breakout

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, κ, 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 ρ_0 respectively, and κ. The total breakout energy, E_BO, and the maximal ejecta velocity, v_max, are shown to be E_BO=8.0πR^2κ^-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/κρ_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, Δ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(Δ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}.

preprint2011arXiv

Radius and magnetic field from Synchrotron-self-absorbed radio and Inverse Compton X-ray observations of Supernovae

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.

preprint2011arXiv

Super-Eccentric Migrating Jupiters

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.

preprint2011arXiv

X-rays, gamma-rays and neutrinos from collisoinless shocks in supernova wind breakouts

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.

preprint2010arXiv

Fast radiation mediated shocks and supernova shock breakouts

We present a simple analytic model for the structure of non-relativistic and relativistic radiation mediated shocks. At shock velocities β_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 β_s\approx 0.2. At higher velocities, β_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 β_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).

preprint2010arXiv

Long wavelength unstable modes in the far upstream of relativistic collisionless shocks

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 Γ_{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.

preprint2010arXiv

Magnetic fields and cosmic rays in clusters of galaxies

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π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.

preprint2010arXiv

Relativistic Radiation Mediated Shocks

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 $Γ_u$ in the range $6\leΓ_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 $ Γ$ drops from $ Γ_u $ to $ \sim 1 $, is characterized by high plasma temperatures $ T\sim Γ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 Γ^2 m_ec^2$, respectively.Photon scattering is dominated by e$^\pm$ pairs, with pair to proton density ratio reaching $\approx10^2Γ_u$. The width of the deceleration region, in terms of Thomson optical depths for upstream going photons, is large, $Δτ\simΓ_u^2$ ($Δτ\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, $νI_ν\proptoν^0$, and an energy cutoff at $ \sim Γ_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.

preprint2010arXiv

The spectrum of Cosmic Rays escaping from relativistic shocks

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 ~ Γ^2 m_pc^2 where Γ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 (ε_B/0.1)(n/1 cm^-3)^(1/6)(E_tot/10^51 erg)^(1/3) eV where ε_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.

preprint2010arXiv

What can we really learn from positron flux 'anomalies'?

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.

preprint2009arXiv

The energy production rate & the generation spectrum of UHECRs

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^-α(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)(α-1) 10^44 erg Mpc^-3 yr^-1 at E<10^19.5 eV with αroughly confined to 2\lesseqα<2.7. The uncertainty is dominated by the systematic and statistic errors in the experimental determination of individual CR event energy, (ΔE/E)_{sys} (ΔE/E)_{stat} ~20%. At lower energy, d\dot{n}/dE is uncertain due to the unknown Galactic contribution. Simple models in which α\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, α\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.

preprint2006arXiv

Self-Similar Collisionless Shocks

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 μof momentum p in diffusive shock acceleration models must satisfy D_{μμ}(p,D) = D^{-1}D'_{μμ}(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).