BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Re'em Sari

Re'em Sari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.EPastro-phastro-ph.GAastro-ph.IMphysics.comp-phastro-ph.COastro-ph.SRphysics.flu-dyngr-qcphysics.geo-ph

Catalog footprint

What is connected

49works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

49 published item(s)

preprint2022arXiv

A new perspective on interiors of ice-rich planets: Ice-rock mixture instead of ice on top of rock

Ice-rich planets are formed exterior to the water ice-line and thus are expected to contain a substantial amount of ices. The high ice content leads to unique conditions in the interior, under which the structure of a planet is affected by ice interaction with other metals. We apply experimental data of ice-rock interaction at high pressure, and calculate detailed thermal evolution for possible interior configurations of ice-rich planets, in the mass range of super-Earth to Neptunes (5-15 Earth masses). We model the effect of migration inward on the ice-rich interior by including the influences of stellar flux and envelope mass loss. We find that ice and rock are expected to remain mixed, due to miscibility at high pressure, in substantial parts of the planetary interior for billions of years. We also find that the deep interior of planetary twins that have migrated to different distances from the star are usually similar, if no mass loss occurs. Significant mass loss results in separation of the water from the rock on the surface and emergence of a volatile atmosphere of less than 1 percent of the planet's mass. The mass of the atmosphere of water/steam is limited by the ice-rock interaction. We conclude that when ice is abundant in planetary interiors the planet structure may differ significantly from the standard layered structure of a water shell on top of a rocky core. Similar structure is expected in both close-in and further-out planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Formation of Intermediate Mass Black Holes in Galactic Nuclei

Most stellar evolution models predict that black holes (BHs) should not exist above approximately $50-70$ M$_\odot$, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate mass black holes (IMBHs), can form in galactic nuclei through collisions between stellar-mass black holes and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as $10^4$ M$_\odot$. This upper limit assumes that (1) the BHs accrete a substantial fraction of the stellar mass captured during each collision and (2) that the rate at which new stars are introduced into the region near the SMBH is high enough to offset depletion by stellar disruptions and star-star collisions. We discuss deviations from these key assumptions in the text. Our results suggest that BHs in the pair-instability mass gap and IMBHs may be ubiquitous in galactic centers. This formation channel has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from x-ray binaries and tidal disruption events. Additionally, formed through this channel, both black holes in the mass gap and IMBHs can merge with the supermassive black hole at the center of a galactic nucleus through gravitational waves. These gravitational wave events are extreme and intermediate mass ratio inspirals (EMRIs and IMRIs, respectively).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Partial Stellar Explosions -- Ejected Mass and Minimal Energy

Many massive stars appear to undergo enhanced mass loss during late stages of their evolution. In some cases, the ejected mass likely originates from non-terminal explosive outbursts, rather than continuous winds. Here we study the dependence of the ejecta mass, $m_{\rm ej}$, on the energy budget $E$ of an explosion deep within the star, using both analytical arguments and numerical hydrodynamics simulations. Focusing on polytropic stellar models, we find that for explosion energies smaller than the stellar binding energy, the ejected mass scales as $m_{\rm ej} \propto E^{\varepsilon_{m}}$, where $\varepsilon_m = 2.4-3.0$ depending on the polytropic index. The loss of energy due to shock breakout emission near the stellar edge leads to the existence of a minimal mass-shedding explosion energy, corresponding to a minimal ejecta mass. For a wide range of progenitors, from Wolf-Rayet stars to red supergiants, we find a similar limiting energy of $E_{\rm min} \approx 10^{46}-10^{47} \rm \, erg$, almost independent of the stellar radius. The corresponding minimal ejecta mass varies considerably across different progenitors, ranging from $\sim \! 10^{-8} \, \rm M_\odot$ in compact stars, up to $\sim \! 10^{-2} \, \rm M_\odot$ in red supergiants. We discuss implications of our results for pre-supernova outbursts driven by wave heating, and complications caused by the non-constant opacity and adiabatic index of realistic stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Blackbody Radiation from Isolated Neptunes

Recent analyses of the orbits of some Kuiper Belt objects hypothesize the presence of an undiscovered Neptune-size planet at a very large separation from the Sun. The energy budget of Neptunes on such distant orbits is dominated by the internal heat released by their cooling rather than solar irradiation (making them effectively "isolated"). The blackbody radiation that these planets emit as they cool may provide the means for their detection. Here we use an analytical toy model to study the cooling and radiation of isolated Neptunes. This model can translate a detection (or a null detection) to a constraint on the size and composition of the hypothesized "Planet Nine". Specifically, the thick gas atmosphere of Neptune-like planets serves as an insulating blanket which slows down their cooling. Therefore, a measurement of the blackbody temperature, $T_{\rm eff}\sim 50\textrm{K}$, at which a Neptune emits can be used to estimate the mass of its atmosphere, $M_{\rm atm}$. Explicitly, we find the relation $T_{\rm eff}\propto M_{\rm atm}^{1/12}$. Despite this weak relation, a measurement of the flux at the Wien tail can constrain the atmospheric mass, at least to within a factor of a few, and provide useful limits to possible formation scenarios of these planets. Finally, we constrain the size and composition of Planet Nine by combining our model with the null results of recent all-sky surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extended Heat Deposition in Hot Jupiters: Application to Ohmic Heating

Many giant exoplanets in close orbits have observed radii which exceed theoretical predictions. One suggested explanation for this discrepancy is heat deposited deep inside the atmospheres of these "hot Jupiters". Here, we study extended power sources which distribute heat from the photosphere to the deep interior of the planet. Our analytical treatment is a generalization of a previous analysis of localized "point sources". We model the deposition profile as a power law in the optical depth and find that planetary cooling and contraction halt when the internal luminosity (i.e. cooling rate) of the planet drops below the heat deposited in the planet's convective region. A slowdown in the evolutionary cooling prior to equilibrium is possible only for sources which do not extend to the planet's center. We estimate the Ohmic dissipation resulting from the interaction between the atmospheric winds and the planet's magnetic field, and apply our analytical model to Ohmically heated planets. Our model can account for the observed radii of most inflated planets which have equilibrium temperatures $\approx 1500\textrm{K}-2500\textrm{K}$, and are inflated to a radius $\approx 1.6 R_J$. However, some extremely inflated planets remain unexplained by our model. We also argue that Ohmically inflated planets have already reached their equilibrium phase, and no longer contract. Following Wu & Lithwick who argued that Ohmic heating could only suspend and not reverse contraction, we calculate the time it takes Ohmic heating to re-inflate a cold planet to its equilibrium configuration. We find that while it is possible to re-inflate a cold planet, the re-inflation timescales are longer by a factor of $\approx 30$ than the cooling time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grid Noise in Moving Mesh Codes: Fixing The Volume Inconsistency Problem

Current Voronoi based moving mesh hydro codes suffer from "grid noise". We identify the cause of this noise as the volume inconsistency error, where the volume that is transferred between cells is inconsistent with the hydrodynamical calculations. As a result, the codes do not achieve second order convergence. In this paper we describe how a simple fix allows Voronoi based moving mesh codes to attain second order convergence. The fix is based on the understanding that the volume exchanged between cells should be consistent with the hydrodynamical calculations. We benchmark our fix with three test problems and show that it can significantly improve the computational accuracy. We also examine the effect of initial mesh initialization and present an improved model for the Green-Gauss based gradient estimator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Super-Earth Atmospheres: Self-Consistent Gas Accretion and Retention

Some recently discovered short-period Earth to Neptune sized exoplanets (super Earths) have low observed mean densities which can only be explained by voluminous gaseous atmospheres. Here, we study the conditions allowing the accretion and retention of such atmospheres. We self-consistently couple the nebular gas accretion onto rocky cores and the subsequent evolution of gas envelopes following the dispersal of the protoplanetary disk. Specifically, we address mass-loss due to both photo-evaporation and cooling of the planet. We find that planets shed their outer layers (dozens of percents in mass) following the disk's dispersal (even without photo-evaporation), and their atmospheres shrink in a few Myr to a thickness comparable to the radius of the underlying rocky core. At this stage, atmospheres containing less particles than the core (equivalently, lighter than a few % of the planet's mass) can be blown away by heat coming from the cooling core, while heavier atmospheres cool and contract on a timescale of Gyr at most. By relating the mass-loss timescale to the accretion time, we analytically identify a Goldilocks region in the mass-temperature plane in which low-density super Earths can be found: planets have to be massive and cold enough to accrete and retain their atmospheres, while not too massive or cold, such that they do not enter runaway accretion and become gas giants (Jupiters). We compare our results to the observed super-Earth population and find that low-density planets are indeed concentrated in the theoretically allowed region. Our analytical and intuitive model can be used to investigate possible super-Earth formation scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Spectral Sharpness Angle of Gamma-ray Bursts

We explain the results of Yu et al. (2015b) of the novel sharpness angle measurement to a large number of spectra obtained from the Fermi gamma-ray burst monitor. The sharpness angle is compared to the values obtained from various representative emission models: blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. It is found that more than 91% of the high temporally and spectrally resolved spectra are inconsistent with any kind of optically thin synchrotron emission model alone. It is also found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up to 58+23 -18% of the peak flux. These results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain a large fraction of the observed spectra. Since any combination of physically possible synchrotron spectra added together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tidal Heating of Young Super-Earth Atmospheres

Short-period Earth to Neptune size exoplanets (super-Earths) with voluminous gas envelopes seem to be very common. These gas atmospheres are thought to have originated from the protoplanetary disk in which the planets were embedded during their first few Myr. The accretion rate of gas from the surrounding nebula is determined by the ability of the gas to cool and radiate away its gravitational energy. Here we demonstrate that heat from the tidal interaction between the star and the young (and therefore inflated) planet can inhibit the gas cooling and accretion. Quantitatively, we find that the growth of super-Earth atmospheres halts for planets with periods of about 10 days, provided that their initial eccentricities are of the order of 0.2. Thus, tidal heating provides a robust and simple mechanism that can simultaneously explain why these planets did not become gas giants and account for the deficit of low-density planets closer to the star, where the tides are even stronger. We suggest that tidal heating may be as important as other factors (such as the nebula's lifetime and atmosphere evaporation) in shaping the observed super-Earth population.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotic Steady State Solution to a Bow Shock with an Infinite Mach Number

The problem of a cold gas flowing past a stationary object is considered. It is shown that at large distances from the obstacle the shock front forms a parabolic solid of revolution. The interior of the shock front is obtained by solution of the hydrodynamic equations in parabolic coordinates. The results are verified with a hydrodynamic simulation. The drag force and expected spectra are calculated for such shock, both in case of an optically thin and thick media. Finally, relations to astrophysical bow shocks and other analytic works on oblique shocks are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy Injection in Gamma-ray Burst Afterglows

We present multi-wavelength observations and modeling of Gamma-ray Bursts (GRBs) that exhibit a simultaneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming corrected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both z~1 and z>6, suggesting that the jet launching mechanism and environment of these events are similar to that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simultaneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of energy injection, forming the most extreme events in their class. Future X-ray observations of GRB afterglows with Swift and its successors will thus likely discover several more such events, while radio follow-up and multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hot-Jupiter Inflation due to Deep Energy Deposition

Some extrasolar giant planets in close orbits---"hot Jupiters"---exhibit larger radii than that of a passively cooling planet. The extreme irradiation $L_{\rm eq}$ these hot Jupiters receive from their close in stars creates a thick isothermal layer in their envelopes, which slows down their convective cooling, allowing them to retain their inflated size for longer. This is yet insufficient to explain the observed sizes of the most inflated planets. Some models invoke an additional power source, deposited deep in the planet's envelope. Here we present an analytical model for the cooling of such irradiated, and internally heated gas giants. We show that a power source $L_{\rm dep}$, deposited at an optical depth $τ_{\rm dep}$, creates an exterior convective region, between optical depths $L_{\rm eq}/L_{\rm dep}$ and $τ_{\rm dep}$, beyond which a thicker isothermal layer exists, which in extreme cases may extend to the center of the planet. This convective layer, which occurs only for $L_{\rm dep}τ_{\rm dep}>L_{\rm eq}$, further delays the cooling of the planet. Such a planet is equivalent to a planet irradiated with $L_{\rm eq}\left(1+L_{\rm dep}τ_{\rm dep}/L_{\rm eq}\right)^β$, where $β\approx 0.35$ is an effective power-law index describing the radiative energy density as function of the optical depth for a convective planet $U\proptoτ^β$. Our simple analytical model reproduces the main trends found in previous numerical works, and provides an intuitive understanding. We derive scaling laws for the cooling rate of the planet, its central temperature, and radius. These scaling laws can be used to estimate the effects of tidal or Ohmic dissipation, wind shocks, or any other mechanism involving energy deposition, on sizes of hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The afterglow of a relativistic shock breakout and low-luminosity GRBs

The prompt emission of low-luminosity gamma-ray bursts (llGRBs) indicates that these events originate from a relativistic shock breakout. In this case we can estimate, based on the properties of the prompt emission, the energy distribution of the ejecta. We develop a general formalism to estimate the afterglow produced by synchrotron emission from the forward shock resulting from the interaction of this ejecta with the circum-burst matter. We assess whether this emission can produce the observed radio and X-ray afterglows of the available sample of 4 llGRBs. All 4 radio afterglows can be explained within this model, providing further support for shock breakouts being the origin of llGRBs. We find that in one of the llGRBs (GRB 031203) the predicted X-ray emission, using the same parameters that fit the radio, can explain the observed one. In another one (GRB 980425) the observed X-rays can be explained if we allow for a slight modification of the simplest model. For the last two cases (GRBs 060218 and 100316D), we find that, as is the case for previous attempts to model these afterglows, the simplest model that fits the radio emission underpredicts the observed X-ray afterglows. Using general arguments, we show that the most natural location of the X-ray source is, like the radio source, within the ejecta-external medium interaction layer but that emission is due to a different population of electrons or to a different emission process.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The sharpness of gamma-ray burst prompt emission spectra

We aim to obtain a measure of the curvature of time-resolved spectra that can be compared directly to theory. This tests the ability of models such as synchrotron emission to explain the peaks or breaks of GBM prompt emission spectra. We take the burst sample from the official Fermi GBM GRB time-resolved spectral catalog. We re-fit all spectra with a measured peak or break energy in the catalog best-fit models in various energy ranges, which cover the curvature around the spectral peak or break, resulting in a total of 1,113 spectra being analysed. We compute the sharpness angles under the peak or break of the triangle constructed under the model fit curves and compare to the values obtained from various representative emission models: blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. We find that 35% of the time-resolved spectra are inconsistent with the single-electron synchrotron function, and 91% are inconsistent with the Maxwellian synchrotron function. The single temperature, single emission time and location blackbody function is found to be sharper than all the spectra. No general evolutionary trend of the sharpness angle is observed, neither per burst nor for the whole population. It is found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up to $58^{+23}_{-18}$% of the peak flux. Our results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain a large fraction of the observed GRB prompt spectra. Because of the fact that any combination of physically possible synchrotron spectra added together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular Density Perturbations to Filled Type I Strong Explosions

In this paper we extend the Sedov - Taylor - Von Neumann model for a strong explosion to account for small angular and radial variations in the density. We assume that the density profile is given by $ρ\left(r,θ,ϕ\right)=kr^{-ω}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(θ,ϕ\right)\right)$, where $\varepsilon\ll1$ and $ω\le\frac{7-γ}{γ+1}$. In order to verify our results we compare them to analytical approximations and full hydrodynamic simulations. We demonstrate how this method can be used to describe arbitrary (not just self similar) angular perturbations. This work complements our previous analysis on radial, spherically symmetric perturbations, and allows one to calculate the response of an explosion to arbitrary perturbations in the upstream density. Together, they settle an age old controversy about the inner boundary conditions

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atmospheric Mass Loss During Planet Formation: The Importance of Planetesimal Impacts

We quantify the atmospheric mass loss during planet formation by examining the contributions to atmospheric loss from both giant impacts and planetesimal accretion. Giant impacts cause global motion of the ground. Using analytic self-similar solutions and full numerical integrations we find (for isothermal atmospheres with adiabatic index ($γ=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $χ_{loss}=(1.71 v_g/v_{esc})^{4.9}$, where $v_{esc}$ is the escape velocity from the target. Yet, the global atmospheric mass loss is a weaker function of the impactor velocity $v_{Imp}$ and mass $m_{Imp}$ and given by $X_{loss} ~ 0.4x+1.4x^2-0.8x^3$ (isothermal atmosphere) and $X_{loss} ~ 0.4x+1.8x^2-1.2x^3$ (adiabatic atmosphere), where $x=(v_{Imp}m/v_{esc}M)$. Atmospheric mass loss due to planetesimal impacts proceeds in two different regimes: 1) Large enough impactors $m > \sqrt{2} ρ_0 (πh R)^{3/2}$ (25~km for the current Earth), are able to eject all the atmosphere above the tangent plane of the impact site, which is $h/2R$ of the whole atmosphere, where $h$, $R$ and $ρ_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 πρ_0 h^3$ (1~km for the current Earth) are only able to eject a fraction of the atmospheric mass above the tangent plane. We find that the most efficient impactors (per unit impactor mass) for atmospheric loss are planetesimals just above that lower limit and that the current atmosphere of the Earth could have resulted from an equilibrium between atmospheric erosion and volatile delivery to the atmosphere from planetesimals. We conclude that planetesimal impacts are likely to have played a major role in atmospheric mass loss over the formation history of the terrestrial planets. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Balancing the load: A Voronoi based scheme for parallel computations

The use of numerical simulations in science is ever increasing and with it the computational size. In many cases single processors are no longer adequate and simulations are run on multiple core machines or supercomputers. One of the key issues when running a simulation on multiple CPUs is maintaining a proper load balance throughout the run and minimizing communications between CPUs. We propose a novel method of utilizing a Voronoi diagram to achieve a nearly perfect load balance without the need of any global redistributions of data. As a show case, we implement our method in RICH, a 2D moving mesh hydrodynamical code, but it can be extended trivially to other codes in 2D or 3D. Our tests show that this method is indeed efficient and can be used in a large variety of existing hydrodynamical codes as well as other applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recombination Effects on Supernovae Light-Curves

Supernovae of type IIP are marked by the long plateau seen in their optical light curves. The plateau is believed to be the result of a recombination wave that propagates through the outflowing massive hydrogen envelope. Here, we analytically investigate the transition from a fully ionized envelope to a partially recombined one and its effects on the SN light curve. The motivation is to establish the underlying processes which dominate the evolution at late times when recombination takes place in the envelope, yet early enough so that $^{56}$Ni decay is a negligible source of energy. We assume a simple, yet adequate, hydrodynamic profile of the envelope and study the mechanisms which dominate the energy emission and the observed temperature. We consider the diffusion of photons through the envelope while analyzing the ionization fraction and the coupling between radiation and gas. We find that once recombination starts, the observed temperature decreases slowly in time. However, in a typical red supergiant (RSG) explosion, the recombination wave does not affect the bolometric luminosity immediately. Only at later times, the cooling wave may reach layers that are deep enough to affect the luminosity. We find that the plateau is not a generic result of a recombination process in expanding gas. Instead it depends on the density profile of the parts of the envelope which undergo recombination. Our results are useful to investigate the light curves of RSG explosions. We show the resulting light curves of two examples of RSG explosions according to our model and discuss their compatibility with observations. In addition, we improve the analytical relations between the plateau luminosity and plateau duration to the properties of the pre-explosion progenitor (Arnett 1980; Popov 1993).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh

We present here RICH, a state of the art 2D hydrodynamic code based on Godunov's method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time advancement scheme as well as a novel parallelization scheme based on Voronoi tessellation. Using our code we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time advancement scheme is more accurate and robust than AREPO's, when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Llyod iterations) and it effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way, and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. Moreover, we show that Voronoi based moving mesh schemes suffer from an error, that is resolution independent, due to inconsistencies between the flux calculation and change in the area of a cell. Our code is publicly available as open source and designed in an object oriented, user friendly way that facilitates incorporation of new algorithms and physical processes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Galactic Center cloud G2 and its gas streamer

We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Velocity Distribution of Hypervelocity Stars

We consider the process of stellar binaries tidally disrupted by a supermassive black hole. For highly eccentric orbits, as one star is ejected from the three-body system, the companion remains bound to the black hole. Hypervelocity stars (HVSs) observed in the Galactic halo and S-stars observed orbiting the central black hole may originate from such mechanism. In this paper, we predict the velocity distribution of the ejected stars of a given mass, after they have travelled out of the Galactic potential. We use both analytical methods and Monte Carlo simulations. We find that each part of the velocity distribution encodes different information. At low velocities < 800 km/s, the Galactic Potential shapes universally the observed distribution, which rises towards a peak, related to the Galactic escape velocity. Beyond the peak, the velocity distribution depends on binary mass and separation distributions. Finally, the finite star life introduces a break related to their mass. A qualitative comparison of our models with current observations shows the great potential of HVSs to constrain bulge and Galactic properties. Standard choices for parameter distributions predict velocities below and above ~800 km/s with equal probability, while none are observed beyond ~700 km/s and the current detections are more clustered at low velocities 300-400 km/s. These features may indicate that the separation distribution of binaries that reach the tidal sphere is not flat in logarithmic space, as observed in more local massive binaries, but has more power towards larger separations, enhancing smaller velocities. In addition, the binary formation/evolution process or the injection mechanism might also induce a cut-off at ~ 10 solar radii in the separation distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new rare type of supernovae: hypervelocity stellar collisions at galactic centers

When a binary star system is tidally disrupted by a supermassive black hole at a galactic nucleus, one star is ejected at a high speed while the other remains in a tightly bound orbit around the black hole. The cluster of tightly bound stars builds over time, eventually creating a steady state in which the rate of collisions between these stars is similar to the rate of capturing new stars. A large fraction of the collisions occur near the periapsis of the orbits around the black hole, where the kinetic energies are sufficient to generate an explosive disruption of the two stars involved. The typical flare brightens for several days, with a peak luminosity that is comparable to the lower-luminosity end of known supernovae. The explosion lightcurve is followed by a longer flare due to accretion of ejected matter onto the black hole. Dedicated searches in the near universe could observe several such "collisional-supernovae" per year.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Consequences of Strong Compression in Tidal Disruption Events

The tidal disruption of a star by a supermassive black hole (SMBH) is a highly energetic event with consequences dependent on the degree to which the star plunges inside the SMBH's tidal sphere. We introduce a new analytic model for tidal disruption events (TDEs) to analyze the dependence of these events on beta, the ratio of the tidal radius to the orbital pericenter. We find, contrary to most previous work, that the spread in debris energy for a TDE is largely constant for all beta. This result has important consequences for optical transient searches targeting TDEs, which we discuss. We quantify leading-order general relativistic corrections to this spread in energy and find that they are small. We also examine the role of stellar spin, and find that a combination of spin-orbit misalignment, rapid rotation, and high beta may increase the spread in debris energy. Finally, we quantify for the first time the gravitational wave emission due to the strong compression of a star in a high-beta TDE. Although this signal is unlikely to be detectable for disruptions of main sequence stars, the tidal disruption of a white dwarf by an intermediate mass black hole can produce a strong signal visible to Advanced LIGO at tens of megaparsecs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ejection and Capture Dynamics in Restricted Three-Body Encounters

We study the tidal disruption of binaries by a massive point mass (e.g. the black hole at the Galactic center), and we discuss how the ejection and capture preference between unequal-mass binary members depends on which orbit they approach the massive object. We show that the restricted three-body approximation provides a simple and clear description of the dynamics. The orbit of a binary with mass m around a massive object M should be almost parabolic with an eccentricity |1-e| < (m/M)^{1/3} << 1 for a member to be captured, while the other is ejected. Indeed, the energy change of the members obtained for a parabolic orbit can be used to describe non-parabolic cases. If a binary has an encounter velocity much larger than (M/m)^{1/3} times the binary rotation velocity, it would be abruptly disrupted, and the energy change at the encounter can be evaluated in a simple disruption model. We evaluate the probability distributions for the ejection and capture of circular binary members and for the final energies. In principle, for any hyperbolic (elliptic) orbit, the heavier member has more chance to be ejected (captured), because it carries a larger fraction of the orbital energy. However, if the orbital energy is close to zero, the difference between the two members becomes small, and there is practically no ejection and capture preference. The preference becomes significant when the orbital energy is comparable to the typical energy change at the encounter. We discuss its implications to hypervelocity stars and irregular satellites around giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Measuring the Abundance of sub-kilometer sized Kuiper Belt Objects using Stellar Occultations

We present here the analysis of about 19,500 new star hours of low ecliptic latitude observations (|b| < 20 deg) obtained by the Hubble Space Telescope's FGS over a time span of more than nine years; which is an addition to the 12,000 star hours previously analyzed by Schlichting et al. (2009). Our search for stellar occultations by small Kuiper belt objects (KBOs) yielded one new candidate event corresponding to a body with a 530 +/-70m radius at a distance of about 40AU. Using bootstrap simulations, we estimate a probability of approx 5%, that this event is due to random statistical fluctuations within the new data set. Combining this new event with the single KBO occultation reported by Schlichting et al. (2009) we arrive at the following results: 1) The ecliptic latitudes of 6.6 deg and 14.4 deg of the two events are consistent with the observed inclination distribution of larger, 100km-sized KBOs. 2) Assuming that small, sub-km sized KBOs have the same ecliptic latitude distribution as their larger counterparts, we find an ecliptic surface density of KBOs with radii larger than 250m of N(r>250m) = 1.1^{+1.5}_{-0.7} x 10^7 deg^{-2}; if sub-km sized KBOs have instead a uniform ecliptic latitude distribution for -20 deg < b< 20 deg then N(r>250m) = 4.4^{+5.8}_{-2.8} x 10^6 deg^{-2}. This is the best measurement of the surface density of sub-km sized KBOs to date. 3) Assuming the KBO size distribution can be well described by a single power law given by N(>r) \propto r^{1-q}, where N(>r) is the number of KBOs with radii greater than r, and q is the power law index, we find q=3.8+/-0.2 for a KBO ecliptic latitude distribution that follows the observed distribution for larger, 100-km sized KBOs. 4) Regardless of the exact power law, our results suggest that small KBOs are numerous enough to satisfy the required supply rate for the Jupiter family comets. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optical to X-rays supernovae light curves following shock breakout through a thick wind

Recent supernovae (SNe) observations have motivated renewed interest in SN shock breakouts from stars surrounded by thick winds. In such events the interaction with the wind powers the observed luminosity, and predictions include observable hard X-rays. Wind breakouts on timescales of a day or longer are currently the most probable for detection. Here we study the signal that follows such events. We start from the breakout of the radiation mediated shock, finding that the breakout temperature can vary significantly from one event to another ($10^4 - 5 \times 10^6$ K) due to possible deviation from thermal equilibrium. In general, events with longer breakout pulse duration, $t_{bo}$, are softer. We follow the observed radiation through the evolution of the collisionless shock which forms after the breakout of the radiation mediated shock. We restrict the study of the collisionless shock evolution to cases where the breakout itself is in thermal equilibrium, peaking in optical/UV. In these cases the post-breakout emission contains two spectral components - soft (optical/UV) and hard (X-rays and possibly soft $γ$-rays). Right after the breakout pulse X-rays are strongly suppressed, and they carry only a small fraction of the total luminosity. The hard component becomes harder and its luminosity rises quickly afterwards, gaining dominance at $\sim 10-50 \,t_{bo}$. The ratio of the peak optical/UV to the peak X-ray luminosity depends mostly on the breakout time. In terms of prospects for X-ray and soft gamma-ray detections, it is best to observe 100-500 days after explosions with breakout timescales between a week and a month.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Short vs Long and Collapsars vs. non-Collapsar: a quantitative classification of GRBs

Gamma-Ray Bursts (GRBs) are traditionally divided to long and short according to their durations (>/< 2 sec). It was generally believed that this reflects a different physical origin: Collapsars (long) and non-Collapsars (short). We have recently shown that the duration distribution of Collapsars is flat, namely independent of the duration, at short durations. Using this model for the distribution of Collapsars we determine the duration distribution of non-Collapsars and estimate the probability that a burst with a given duration (and hardness) is a Collapsar or not. We find that this probability depends strongly on the spectral window of the observing detector. While the commonly used limit of 2 sec is conservative and suitable for BATSE bursts, 40% of swift's bursts shorter than 2 sec are Collapsars and division >/<0.8 sec is more suitable for swift. We find that the duration overlap of the two populations is very large. On the one hand there is a non-negligible fraction of non-Collapsars longer than 10 sec, while on the other hand even bursts shorter than 0.5 sec in the swift sample have a non-negligible probability to be Collapsars. Our results enable the construction of non-Collapsar samples while controlling the Collapsar contamination. They also highlight that no firm conclusions can be drawn based on a single burst and they have numerous implications concerning previous studies of non-Collapsar properties that were based on the current significantly contaminated swift samples of localized short GRBs. Specifically: (i) all known short bursts with z>1 are most likely Collapsars, (ii) the only short burst with a clear jet break is most likely a Collapsar, indicating our lack of knowledge concerning non-Collapsar beaming (iii) the existence of non-Collapsars with durations up to 10 sec impose new challenges to non-Collapsar models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sterile and Fertile Planetary Systems - Statistical Analysis of Multi-Planet Systems in Kepler's data

The Kepler mission has discovered a large number of planetary systems. We analyze the implications of the discovered single/multi-exoplanet systems from Kepler's data. As done in previous works, we test a simple model in which the intrinsic occurrence of planet is an independent process, and with equal probability around all planet producing stars. This leads to a Poisson distribution for the intrinsic number of planets around each host. However, the possibility of zero/low mutual inclination is taken into account, creating a correlation between detecting different planets in a given stellar system, leading to a non Poisson distribution for the number of transiting planets per system. Comparing the model's predictions with the observations made by Kepler, we find that the correlation produced by planarity is insufficient and a higher correlation is needed; either the formation of one planet in the system enhances the likelihood of other planets to form, and/or that some stars are considerably more fertile than others. Kepler's data presents evidences that both correlations might play a part, in particular a significant dependency in the radial distribution of planets in multi-planet systems is shown. Followup observations on Kepler planet's hosts can help pinpoint the physical nature of this correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Long, the Short and the Weak - the origin of GRBs

The origin of Gamma-Ray Bursts is one of the most interesting puzzles in recent astronomy. During the last decade a consensus formed that long GRBs (LGRBs) arise from the collapse of massive stars and that short GRBs (SGRBs) have a different origin, most likely neutron star mergers. A key ingredient of the Collapsar model that explains how the collapse of massive stars produces a GRB is the emergence of a relativistic jet that penetrates the stellar envelope. The condition that the emerging jet penetrates the envelope poses strong constraints on the system. Using these constraints we show that: (i) Low luminosity GRBs (llGRBs), a sub population of GRBs with a very low luminosities (and other peculiar properties: single peaked, smooth and soft) cannot be formed by Collapsars. llGRBs must have a different origin (most likely a shock breakout). (ii) On the other hand regular LGRBs must be formed by Collapsars. (iii) While for BATSE the dividing duration between Collapsars and non-Collapsar is indeed at $\sim 2$ sec, the dividing duration is different for other GRBs detectors. In particular most Swift bursts longer than 0.8 sec are of a Collapsar origin. This last results requires a revision of many conclusions concerning the origin of Swift SGRBs which were based on the commonly used 2 sec limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An observational imprint of the Collapsar model of long Gamma Ray Bursts

The Collapsar model provides a theoretical framework for the well known association between long gamma-ray bursts (GRBs) and collapsing massive stars. A bipolar relativistic jet, launched at the core of a collapsing star, drills its way through the stellar envelope and breaks out of the surface before producing the observed gamma-rays. While a wealth of observations associate GRBs with the death of massive stars, as yet there is no direct evidence for the Collapsar model itself. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all three major satellites. These findings provide an evidence that directly supports the Collapsar model. Additionally, it suggests the existence of a large population of choked (failed) GRBs and that the 2 s duration commonly used to separate Collapsars and non-Collapasars is inconsistent with the duration distributions of {\it Swift} and Fermi GRBs and only holds for BATSE GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

New insights on the duration distribution of long GRBs from Collapsars

According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the γ-rays after they break out. This model provides a theoretical framework for the well known association between LGRBs and massive stars. However although this association is supported by a wealth of observations, to this date there is no direct observational evidence for the emergence of the jet from the star. In other words there is no direct evidence for the Collapsar model. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all major GRB satellites, and provides a direct evidence supporting the Collapsar model. It also enables us to place limits on the sizes and masses of LGRB progenitors; suggests the existence of a large population of choked (failed) GRBs; and indicates that the 2 s duration commonly used to separate Collapsars and non-Collapsars holds for BATSE and possibly Fermi GBM GRBs, but it is inconsistent with the duration distributions of Swift GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Relativistic shock breakouts - a variety of gamma-ray flares: from low luminosity gamma-ray bursts to type Ia supernovae

A large variety of explosions result in mildly or ultra relativistic shock breakouts. Here we calculate the luminosity and spectrum that these breakouts produce. In order to do so we improve an analytic description of relativistic radiation mediated shocks and follow the system from the breakout itself, through the planar phase and into the spherical phase. We limit our calculation to cases where the post breakout acceleration of the gas ends during the planar phase. (i.e., the final gas Lorentz factor <~30). We find that spherical relativistic breakouts produce a flash of gamma-rays with energy, temperature, and duration that always satisfy a specific relation between them as well as providing the breakout radius and Lorentz factor. The breakout flare is typically followed, on longer time scales, by X-rays that carry a comparable energy. We apply our model to a variety of explosions, including Ia and .Ia supernovae (SNe), AIC, energetic SNe and gamma-ray bursts (GRBs). We find that all these events produce detectable gamma-ray signals, some of which may have already been seen. Some particular examples are: (i) Relativistic shock breakouts provide a natural explanation to the energy, temperature and time scales, as well as many other observed features, of all low luminosity GRBs. (ii) Nearby broad-line Ib/c (like SN 2002ap) may produce a detectable gamma-ray signal. (iii) Galactic Ia SNe produce detectable gamma-ray flares, if their progenitors are single degenerate. We conclude that relativistic shock breakout is a generic process for the production of gamma-ray flares, which opens a new window for the study and detection of a variety of stellar explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Slightly Two or Three Dimensional Self-Similar Solutions

Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The propagation of relativistic jets in external media

Relativistic jets are ubiquitous in astrophysical systems that contain compact objects. They transport large amounts of energy to large distances from the source, and their interaction with the ambient medium has a crucial effect on the evolution of the system. The propagation of the jet is characterized by the formation of a shocked "head" at the front of the jet which dissipates the jet's energy and a cocoon that surrounds the jet and potentially collimates it. We present here a self consistent, analytic model that follows the evolution of the jet and its cocoon, and describes their interaction. We show that the critical parameter that determines the properties of the jet-cocoon system is the dimensionless ratio between the jet's energy density and the rest-mass energy density of the ambient medium. This parameter, together with the jet's injection angle, also determines whether the jet is collimated by the cocoon or not. The model is applicable to relativistic, unmagnetized, jets on all scales and may be used to determine the conditions in AGNs jets as well as in GRBs or microquasars. It shows that AGN and microquasar jets are hydrodynamically collimated due to the interaction with the ambient medium, while GRB jets can be collimated only inside a star and become uncollimated once they breakout.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Unified Theory for the Effects of Stellar Perturbations and Galactic Tides on Oort Cloud Comets

We examine the effects of passing field stars on the angular momentum of a nearly radial orbit of an Oort cloud comet bound to the Sun. We derive the probability density function (PDF) of the change in angular momentum from one stellar encounter, assuming a uniform and isotropic field of perturbers. We show that the total angular momentum follows a Levy flight, and determine its distribution function. If there is an asymmetry in the directional distribution of perturber velocities, the marginal probability distribution of each component of the angular momentum vector can be different. The constant torque attributed to Galactic tides arises from a non-cancellation of perturbations with an impact parameter of order the semimajor axis of the comet. When the close encounters are rare, the angular momentum is best modeled by the stochastic growth of stellar encounters. If trajectories passing between the comet and sun occur frequently, the angular momentum exhibits the coherent growth attributed to the Galactic tides.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Binary YORP and Evolution of Binary Asteroids

The rotation states of kilometer sized near earth asteroids are known to be affected by the YORP effect. In a related effect, Binary YORP (BYORP) the orbital properties of a binary asteroid evolves under a radiation effect mostly acting on a tidally locked secondary. The BYORP effect can alter the orbital elements in $\sim 10^{4-5}$ years for a $D_{p}=2\;km$ primary with a $D_{s}=0.4\; km$ secondary at $1\; AU$. It can either separate the binary components or cause them to collide. In this paper we devise a simple approach to calculate the YORP effect on asteroids and BYORP effect on binaries including $J_2$ effects due to primary oblateness and the sun. We apply this to asteroids with known shapes as well as a set of randomly generated bodies with various degrees of smoothness. We find a strong correlation between the strengths of an asteroid's YORP and BYORP effects. Therefore, a statistical knowledge on one, could be used to estimate the effect of the other. We show that the action of BYORP preferentially shrinks rather than expands the binary orbit and that YORP preferentially slows down asteroids. This conclusion holds for the two extremes of thermal conductivities studied in this work and assuming the asteroid reaches a stable point, but may break down for moderate thermal conductivity. The YORP and BYORP effects are shown to be smaller than what could be naively expected due to near cancellation of the effects on small scales. Taking this near cancellation into account, a simple order of magnitude estimate of the YORP and BYORP effects as function of the sizes and smoothness of the bodies is calculated. Finally, we provide a simple proof showing that there is no secular effect due to absorption of radiation in BYORP.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Early supernovae light-curves following the shock-breakout

The first light from a supernova (SN) emerges once the SN shock breaks out of the stellar surface. The first light, typically a UV or X-ray flash, is followed by a broken power-law decay of the luminosity generated by radiation that leaks out of the expanding gas sphere. Motivated by recent detection of emission from very early stages of several SNe, we revisit the theory of shock breakout and the following emission. We derive analytic light curves, paying special attention to the photon-gas coupling and deviations from thermal equilibrium. We then consider the breakout from several SNe progenitors. We find that for more compact progenitors, white dwarfs, Wolf-Rayet stars (WRs) and possibly more energetic blue-supergiant explosions, the observed radiation is out of thermal equilibrium at the breakout, during the planar phase (i.e., before the expanding gas doubles its radius), and during the early spherical phase. Therefore, during these phases we predict significantly higher temperatures than previous analysis that assumed equilibrium. When thermal equilibrium prevails, we find the location of the thermalization depth and its temporal evolution. Our results are useful for interpretation of early SN light curves. Some examples are: (i) Red supergiant SNe have an early bright peak in optical and UV flux, less than an hour after breakout. It is followed by a minimum at the end of the planar phase (about 10 hr), before it peaks again once the temperature drops to the observed frequency range. In contrast WRs show only the latter peak in optical and UV. (ii) Bright X-ray flares are expected from all core-collapse SNe types. (iii) The light curve and spectrum of the initial breakout pulse holds information on the explosion geometry and progenitor wind opacity. Its spectrum in compact progenitors shows a (non-thermal) power-law.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Equilibrium Configurations of Synchronous Binaries: Numerical Solutions and Application to Kuiper-Belt Binary 2001 QG298

We present numerical computations of the equilibrium configurations of tidally-locked homogeneous binaries, rotating in circular orbits. Unlike the classical Roche approximations, we self-consistently account for the tidal and rotational deformations of both components, and relax the assumptions of ellipsoidal configurations and Keplerian rotation. We find numerical solutions for mass ratios q between 1e-3 and 1, starting at a small angular velocity for which tidal and rotational deformations are small, and following a sequence of increasing angular velocities. Each series terminates at an appropriate ``Roche limit'', above which no equilibrium solution can be found. Even though the Roche limit is crossed before the ``Roche lobe'' is filled, any further increase in the angular velocity will result in mass-loss. For close, comparable-mass binaries, we find that local deviations from ellipsoidal forms may be as large as 10-20%, and departures from Keplerian rotation are significant. We compute the light curves that arise from our equilibrium configurations, assuming their distance is >>1 AU (e.g. in the Kuiper Belt). We consider both backscatter (proportional to the projected area) and diffuse (Lambert) reflections. Backscatter reflection always yields two minima of equal depths. Diffuse reflection, which is sensitive to the surface curvature, generally gives rise to unequal minima. We find detectable intensity differences of up to 10% between our light curves and those arising from the Roche approximations. Finally, we apply our models to Kuiper Belt binary 2001 QG298, and find a nearly edge-on binary with a mass ratio q = 0.93 ^{+0.07}_{-0.03}, angular velocity Omega^2/G rho = 0.333+/-0.001 (statistical errors only), and pure diffuse reflection. For the observed period of 2001 QG298, these parameters imply a bulk density, rho = 0.72 +/- 0.04 g cm^-3.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Transits of Transparent Planets - Atmospheric Lensing Effects

Light refracted by the planet's atmosphere is usually ignored in analysis of planetary transits. Here we show that refraction can add shoulders to the transit light curve, i.e., an increase in the observed flux, mostly just before and after transit. During transit, light may be refracted away from the observer. Therefore, even completely transparent planets will display a very similar signal to that of a standard transit, i.e., that of an opaque planet. We provide analytical expression for the amount of additional light deflected towards the observer before the transit, and show that the effect may be as large as $10^{-4}$ of the stellar light and therefore measurable by current instruments. By observing this effect we can directly measure the scale height of the planet's atmosphere. We also consider the attenuation of starlight in the planetary atmosphere due to Rayleigh scattering and discuss the conditions under which the atmospheric lensing effect is most prominent. We show that, for planets on orbital periods larger than about 70 days, the size of the transit is determined by refraction effects, and not by absorption within the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Co-orbital Oligarchy

We present a systematic examination of the changes in semi-major axis caused by the mutual interactions of a group of massive bodies orbiting a central star in the presence of eccentricity dissipation. For parameters relevant to the oligarchic stage of planet formation, dynamical friction keeps the typical eccentricities small and prevents orbit crossing. Interactions at impact parameters greater than several Hill radii cause the protoplanets to repel each other; if the impact parameter is instead much less than the Hill radius, the protoplanets shift slightly in semi-major axis but remain otherwise unperturbed. If the orbits of two or more protoplanets are separated by less than a Hill radius, they are each pushed towards an equilibrium spacing between their neighbors and can exist as a stable co-orbital system. In the shear-dominated oligarchic phase of planet formation we show that the feeding zones contain several oligarchs instead of only one. Growth of the protoplanets in the oligarchic phase drives the disk to an equilibrium configuration that depends on the mass ratio of protoplanets to planetesimals, $Σ/σ$. Early in the oligarchic phase, when $Σ/σ$ is low, the spacing between rows of co-orbital oligarchs are about 5 Hill radii wide, rather than the 10 Hill radii cited in the literature. It is likely that at the end of oligarchy the average number of co-orbital oligarchs is greater than unity. In the outer solar system this raises the disk mass required to form the ice giants. In the inner solar system this lowers the mass of the final oligarchs and requires more giant impacts than previously estimated. This result provides additional evidence that Mars is not an untouched leftover from the oligarchic phase, but must be composed of several oligarchs assembled through giant impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Discrete Self-Similarity in Ultra-Relativistic Type-II Strong Explosions

A solution to the ultra-relativistic strong explosion problem with a non-power law density gradient is delineated. We consider a blast wave expanding into a density profile falling off as a steep radial power-law with small, spherically symmetric, and log-periodic density perturbations. We find discretely self-similar solutions to the perturbation equations and compare them to numerical simulations. These results are then generalized to encompass small spherically symmetric perturbations with arbitrary profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hyper Velocity Stars and the Restricted Parabolic 3-body Problem

Motivated by detections of hypervelocity stars that may originate from the Galactic Center, we revist the problem of a binary disruption by a passage near a much more massive point mass. The six order of magnitude mass ratio between the Galactic Center black hole and the binary stars allows us to formulate the problem in the restricted parabolic three-body approximation. In this framework, results can be simply rescaled in terms of binary masses, its initial separation and binary-to-black hole mass ratio. Consequently, an advantage over the full three-body calculation is that a much smaller set of simulations is needed to explore the relevant parameter space. Contrary to previous claims, we show that, upon binary disruption, the lighter star does not remain preferentially bound to the black hole. In fact, it is ejected exactly in 50% of the cases. Nonetheless, lighter objects have higher ejection velocities, since the energy distribution is independent of mass. Focusing on the planar case, we provide the probability distributions for disruption of circular binaries and for the ejection energy. We show that even binaries that penetrate deeply into the tidal sphere of the black hole are not doomed to disruption, but survive in 20% of the cases. Nor do these deep encounters produce the highest ejection energies, which are instead obtained for binaries arriving to 0.1-0.5 of the tidal radius in a prograde orbit. Interestingly, such deep-reaching binaries separate widely after penetrating the tidal radius, but always approach each other again on their way out from the black hole.[shortened]

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Composite self-similar solutions for relativistic shocks: the transition to cold fluid temperatures

The flow resulting from a strong ultrarelativistic shock moving through a stellar envelope with a polytrope-like density profile has been studied analytically and numerically at early times while the fluid temperature is relativistic--that is, just before and just after the shock breaks out of the star. Such a flow should expand and accelerate as its internal energy is converted to bulk kinetic energy; at late enough times, the assumption of relativistic temperatures becomes invalid. Here we present a new self-similar solution for the post-breakout flow when the accelerating fluid has bulk kinetic Lorentz factors much larger than unity but is cooling through $p/n$ of order unity to subrelativistic temperatures. This solution gives a relation between a fluid element's terminal Lorentz factor and that element's Lorentz factor just after it is shocked. Our numerical integrations agree well with the solution. While our solution assumes a planar flow, we show that corrections due to spherical geometry are important only for extremely fast ejecta originating in a region very close to the stellar surface. This region grows if the shock becomes relativistic deeper in the star.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Eccentricity Evolution for Planets in Gaseous Disks

We investigate the hypothesis that interactions between a giant planet and the disk from which it forms promote eccentricity growth. These interactions are concentrated at discrete Lindblad and corotation resonances. Interactions at principal Lindblad resonances cause the planet's orbit to migrate and open a gap in the disk if the planet is sufficiently massive. Those at first order Lindblad and corotation resonances change the planet's orbital eccentricity. Eccentricity is excited by interactions at external Lindblad resonances which are located on the opposite side of corotation from the planet, and damped by co-orbital Lindblad resonances which overlap the planet's orbit. If the planet clears a gap in the disk, the rate of eccentricity damping by co-orbital Lindblad resonances is reduced. Density gradients associated with the gap activate eccentricity damping by corotation resonances at a rate which initially marginally exceeds that of eccentricity excitation by external Lindblad resonances. But the corotation torque drives a mass flux which reduces the density gradient near the resonance. Sufficient partial saturation of corotation resonances can tip the balance in favor of eccentricity excitation. A minimal initial eccentricity of a few percent is required to overcome viscous diffusion which acts to unsaturate corotation resonances by reestablishing the large scale density gradient. Thus eccentricity growth is a finite amplitude instability. Formally, interactions at the apsidal resonance, which is a special kind of co-orbital Lindblad resonance, appears to damp eccentricity faster than external Lindblad resonances can excite it. However, apsidal waves have such long wavelengths that they do not propagate in protoplanetary disks. This reduces eccentricity damping by the apsidal resonance to a modest level.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Shape of Spectral Breaks in GRB Afterglows

Gamma-Ray Burst (GRB) afterglows are well described by synchrotron emission from relativistic blast waves expanding into an external medium. The blast wave is believed to amplify the magnetic field and accelerate the electrons into a power law distribution of energies promptly behind the shock. These electrons then cool both adiabatically and by emitting synchrotron and inverse Compton radiation. The resulting spectra is known to consist several power law segments, which smoothly join at certain break frequencies. Here, we give a complete description of all possible spectra under those assumptions, and find that there are 5 possible regimes, depending on the ordering of the break frequencies. The flux density is calculated by integrating over the contributions from all the shocked region, using the Blandford McKee solution. This allows us to calculate more accurate expressions for the value of these break frequencies, and describe the shape of the spectral breaks around them. This also provides the shape of breaks in the light curves caused by the passage of a break frequency through the observed band. These new, more exact, estimates are different from more simple calculations by up to a factor of about 70, and describe some new regimes which where previously ignored.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Images and Spectra From the Interior of a Relativistic Fireball

The detection of an afterglow, following a gamma-ray burst (GRB), can be described reasonably well by synchrotron emission from a relativistic spherical expanding blast wave, driven by an expanding fireball. We perform detailed calculations considering the emission from the whole region behind the shock front. We use the Blandford McKee self similar adiabatic solution to describe the fluid behind the shock. Using this detailed model, we derive expressions for the peak flux, and the peak frequency at a given observed time. These expressions provide important numerical corrections to previous, more simplified models. We calculate the observed light curve and spectra for several magnetic field models. We show that both the light curve and the spectra are flat near the peak. This rules out the interpretation of the optical peak of GRB970508 as the peak of the light curve, predicted by the existing fireball models. We calculate the observed image of a GRB afterglow. The observed image is bright near the edge and dimmer at the center, thus creating a ring. The contrast between the edge and the center is larger at high frequencies and the width of the ring is smaller.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Can internal shocks produce the variability in GRBs?

We discuss the possibility that gamma-ray bursts result from internal shocks in an ultra-relativistic matter. Using a simple model we calculate the temporal structure and we estimate the efficiency of this process. In this model the ultra-relativistic matter flow is represented by a succession of shells with random values of the Lorentz factor. We calculate the shocks that take place between those shells and we estimate the resulting emission. Internal shocks can produce the highly variable temporal structure observed in most of the bursts provided that the source emitting the relativistic flow is highly variable. The observed peaks are in almost one to one correlation to the activity of the emitting source. A large fraction of the kinetic energy is converted to radiation. The most efficient case is when an inner engine produces shells with comparable energy but very different Lorentz factors. It also gives the most preferable temporal structure.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Variability in GRBs - A Clue

We show that external shocks cannot produce a variable GRB, unless they are produced by an extremely narrow jets (angular opening of < ~10^{-4}) or if only a small fraction of the shell emits the radiation and the process is very inefficient. Internal shocks can produce the observed complex temporal structure provided that the source itself is variable. In this case, the observed temporal structure reflects the activity of the ``inner engine'' that drives the bursts. This sets direct constraints on it.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Cosmological GRBs: Internal vs. External Shocks

Internal (IS) or external (ES) shocks in relativistic expanding shells are currently the best known mechanism for producing GRBs. We calculate the hydrodynamic conditions and the cooling processes in internal shocks, that occur when one layer overtakes another. We compare these conditions with the conditions in external shocks that occur when the shell encounters the ISM. We find that the complicated temporal structure observed in GRBs could be easily explained in IS scenario. Unlike ES, the observed temporal structure simply reflects the temporal behavior of the internal engine that produces the shell. We also find that opacity and in particular the opacity for pair production poses strong constraints on the parameter space and consequently bursts with very narrow peaks are expected to be optically thick and would not contain GeV photons. Finally we find that as in ES synchrotron is the most likely radiation process.

0 citations0 reviewsNo code linkedOpen access