The large brightness variation in the observed lightcurve of Oumuamua is probably related to its shape, i.e., to the ratio between its longest axis and its shortest axis (aspect ratio). Several approaches found the aspect ratio of Oumuamua to be unusually elongated. Moreover, the spin axis orientation has to be almost perpendicular to the observer in order to obtain such an extreme lightcurve, a configuration which is unlikely. However, interstellar Oumuamua may have different surface properties than we know in our solar system. Therefore, in this work we widen the parameter space for surface properties beyond the asteroid-like models and study its effect on the lightcurve of Oumuamua. We calculate reflection from a rotating ellipsoidal object for four models: Lambertian reflection, specular reflection, single scattering diffusive and backscatter. We then calculate the probability to obtain a lightcurve ratio larger than the observed, as a function of the object aspect ratio, assuming an isotopic spin orientation distribution. We find the elongation of Oumuamua to be less extreme for the Lambertian and specular reflection models. Consequently, the probability to observe the lightcurve ratio of Oumuamua given its unknown spin axis orientation is larger for those models. We conclude that different surface reflection properties may suggest alternatives to the extreme shape of Oumuamua, relieving the need for complicated formation scenario, extreme albedo variation, or unnatural origin. Although the models suggested here are for ideal ellipsoidal shape and ideal reflection method, the results emphasize the importance of surface properties for the derived aspect ratio.

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from $\approx150$ s to $\approx163$ days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and mm-band data reveal a distinct signature at $\approx3.9$ days, which we interpret as reverse shock (RS) emission. Our observations present the first time that a single radio-frequency spectral energy distribution can be decomposed directly into RS and forward shock (FS) components. We perform detailed modeling of the full multiwavelength data set, using Markov Chain Monte Carlo sampling to construct the joint posterior density function of the underlying physical parameters describing the RS and FS synchrotron emission. We uncover and account for all degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized ($\sigma\approx3\times10^{-3}$), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of $\Gamma_0\approx103$ for the GRB jet. Our results support the hypothesis that low-density environments are conducive to the observability of RS emission. We compare our observations to other events with strong RS detections, and find a likely observational bias selecting for longer lasting, non-relativistic reverse shocks. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multi-frequency data sets.

The Galactic Center has been under intense scrutiny in the recent years thanks to the unprecedented missions aiming at measuring the gas and star dynamics near the supermassive black hole (SMBH) and at finding gravitational wave (GW) signatures of inspiralling stellar black holes. In the crowded environment of galactic nuclei, the two-body interactions alter the distribution of stars on long timescales, making them drift in energy and angular momentum. We present a simplified analytical treatment of the scattering processes in galactic stellar nuclei, assuming all stars have the same mass. We have discussed how the interplay between two-body relaxation and gravitational wave emission modifies the slope of the inner stellar cusp within the SMBH sphere of influence, and calculated the rates of tidal disruption events (TDEs) and main-sequence extreme-mass ratio inspirals (MS-EMRIs) of stars that are tidally disrupted by the SMBH. We find that typically the ratio of the TDE and MS-EMRI rates is the square of the ratio of the tidal and Schwarzschild radii. For our Galaxy, this implies that the rate of MS-EMRIs is just about a percent of the TDE rate. We then consider the role of stars injected on highly eccentric orbits in the vicinity of the SMBH due to Hills binary disruption mechanism, and show that the MS-EMRI rate can almost approach the TDE rate if the binary fraction at the SMBH influence radius is close to unity. Finally, we discuss that physical stellar collisions affect a large area of phase space.

When the shock wave generated in a supernova explosion breaks out of the stellar envelope, the first photons, typically in the X-ray to UV range, escape to the observer. Following this breakout emission, radiation from deeper shells diffuses out of the envelope as the supernova ejecta expands. Previous studies have shown that the radiation throughout the planar phase (i.e., before the expanding envelope has doubled its radius) originates in the same mass coordinate, called the `breakout shell'. We derive a self-similar solution for the radiation inside the envelope, and show that this claim is incorrect, and that the diffusion wave propagates logarithmically into the envelope (in Lagrangian sense) rather than remaining at a fixed mass coordinate. The logarithmic correction implies that the luminosity originates in regions where the density is $\sim 10$ times higher than previously thought, where the photon production rate is increased and helps thermalization. We show that this result has significant implications on the observed temperature. In our model, the radiation emitted from blue supergiant and Wolf-Rayet explosions is still expected to be out of thermal equilibrium during the entire planar phase, but the observed temperature will decrease by two orders of magnitude, contrary to previous estimates. Considering the conditions at the end of the planar phase, we also find how the temperature and luminosity transition into the spherical phase.

The emergence of a shock from a medium with a free surface is an important process in various astrophysical phenomena. It generates the first light associated with explosions like supernovae and Gamma Ray Bursts. Most previous works considered planar or spherical geometries, where the shock front is parallel to the surface, and emerges simultaneously from all points. Here we study the hydrodynamics of an oblique planar shock breaking out from the planar surface of a uniform density ideal gas with adiabatic index $\gamma$. We obtain an analytic solution to the flow as a function of the angle between the plane of the shock and the surface $\beta$. We find steady state solutions (in a frame moving with the intersection point of the shock and the surface) up to some critical angle ($\beta_{max}=63.4$ degress for $\gamma=5/3$ and $\beta_{max}=69.3$ degrees for $\gamma=4/3$). We show how this analytic solution can be used in more complicated geometries where the shock is not planar, giving the exact profile of the outermost breakout ejecta. We apply our analytical results to a few realistic problems, such as underwater explosions, detonation under the surface of an asteroid, or off center detonations in a uniform sphere.

The light curves of type-II supernovae (SNe) are believed to be highly affected by recombination of hydrogen that takes place in their envelopes. In this work, 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 that 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 consider the diffusion of photons through the envelope while analyzing the ionization fraction and the coupling between radiation and gas, and find that the main effect of recombination is on the evolution of the observed temperature. Before recombination the temperature decreases relatively fast, while after recombination starts it significantly reduces the rate at which the observed temperature drops with time. This behaviour is the main cause for the observed flattening in the optical bands, where for a typical red supergiant explosion, the recombination wave affects the bolometric luminosity only mildly during most of the photospheric phase. Moreover, the plateau phase observed in some type-II SNe is not a generic result of recombination, and it also depends on the density structure of the progenitor. This is one possible explanation to the different light curve decay rates observed in type II (P and L) SNe.

GW170817 was the first detection of a binary neutron star merger via gravitational waves. The event was observed over a wide range of the electromagnetic spectrum, revealing a thermal kilonova dominating the optical signal during the first $\sim$15 days, and a non-thermal synchrotron emission that has continued to rise $\sim$200 days post-merger, dominating the radio and x-ray emission. At early times, when the kilonova is still dominant, the synchrotron emitting electrons can efficiently cool by up-scattering the kilonova photos through inverse-Compton. Yet, the cooling frequency is not observed up to the X-ray band. This can only be explained if the source is moving at least at a mildly relativistic velocity. We find a lower limit on the source's bulk Lorentz factor of $\Gamma > 2.1$ at $9$ days. This lower limit is model independent, and relies directly on the observed quantities, providing an additional robust evidence to the relativistic motion in this event at early times.

Tidal disruption events of stars by supermassive black holes have so far been discovered months to years after the fact. In this paper we explore the short, faint and hard burst of radiation is emitted at maximum compression, as a result of shock breakout. The detection of this burst can be used to capture tidal disruption events in real time. We verify that shock breakout from main sequence stars produces radiation in the X-ray range, but find that it is difficult to detect using all sky X-ray surveying telescopes. In the case of shock breakout from red giants, most of the radiation is emitted in the UV and visible range, which is significantly easier to detect. A similar burst of UV/optical radiation will also be emitted by stars puffed by tidal heating from a previous passage close to the central black hole. This radiation can be detected by surveys like ZTF and LSST. We calculate detection rates for different types of galactic nuclei. For the case of a very full loss cone we predict a detection rate of once per month with LSST, whereas for the case of a very empty loss cone we predict a rate of once per year with LSST. Evidence from a recent tidal disruption event, ASASSN-14li, seems to favour a very full loss cone, in which case LSST is expected to detect one such event every month.

Planets can open cavities (gaps) in the protoplanetary gaseous discs in which they are born by exerting gravitational torques. Viscosity counters these torques and limits the depletion of the gaps. We present a simple one-dimensional scheme to calculate the gas density profile inside gaps by balancing the gravitational and viscous torques. By generalizing the results of Goodman & Rafikov (2001), our scheme properly accounts for the propagation of angular momentum by density waves. This method allows us to easily study low-viscosity discs, which are challenging for full hydrodynamical simulations. We complement our numerical integration by analytical equations for the gap's steady-state depth and width as a function of the planet's to star's mass ratio $\mu$, the gas disc's aspect ratio $h$, and its Shakura & Sunyaev viscosity parameter $\alpha$. Specifically, we focus on low-mass planets ($\mu<\mu_{\rm th}\equiv h^3$) and identify a new low-viscosity regime, $\alpha<h(\mu/\mu_{\rm th})^5$, in which the classical analytical scaling relations are invalid. Equivalently, this low-viscosity regime applies to every gap that is depleted by more than a factor of $(\mu_{\rm th}/\mu)^3$ relative to the unperturbed density. We show that such gaps are significantly deeper and wider than previously thought, and consequently take a longer time to reach equilibrium.

Transit timing variations (TTVs) are a powerful tool for characterizing the properties of transiting exoplanets. However, inferring planet properties from the observed timing variations is a challenging task, which is usually addressed by extensive numerical searches. We propose a new, computationally inexpensive method for inverting TTV signals in a planetary system of two transiting planets. To the lowest order in planetary masses and eccentricities, TTVs can be expressed as a linear combination of 3 functions, which we call the \textit{TTV modes}. These functions depend only on the planets' linear ephemerides, and can be either constructed analytically, or by performing 3 orbital integrations of the three-body system. Given a TTV signal, the underlying physical parameters are found by decomposing the data as a sum of the TTV modes. We demonstrate the use of this method by inferring the mass and eccentricity of 6 \textit{Kepler} planets that were previously characterized in other studies. Finally we discuss the implications and future prospects of our new method.

We develop an analytic, steady-state model for the gas environment in quiescent galactic nuclei. We assume that the mass is constantly supplied by a spherically symmetric distribution of wind emitting stars, and that gravity is solely due to a central supermassive black hole. We show that at some finite radius, where the Keplerian velocity is comparable to the wind velocity, the bulk velocity vanishes. Matter generated below that radius will be accreted onto the black hole, while matter outside it will escape the system. Under certain conditions, the flow may become supersonic at both domains. We obtain radial profiles of the hydrodynamic variables and verify them using a time-dependent hydrodynamic simulation. We delineate the conditions under which radiative cooling can be neglected, and predict the luminosity and spectrum of the free-free X-ray emission from such a system. We discuss applications of our solution to our own Galactic Centre and other quiescent galactic nuclei.

Binary stars that are on close orbits around massive black holes (MBH) such as Sgr A* in the center of the Milky Way are liable to undergo tidal disruption and eject a hypervelocity star. We study the interaction between such a MBH and circular binaries for general binary orientations and penetration depths (i.e. binaries penetrate into the tidal radius around the BH). We show that for very deep penetrators, all binaries are disrupted when the binary rotation axis is roughly oriented toward the BH or it is in the opposite direction. The surviving chance becomes significant when the angle between the binary rotation axis and the BH direction is between \pi /4 and 3 \pi /4. The surviving chance is as high as $\sim$ 20$\%$ when the binary rotation axis is perpendicular to the BH direction. The angular dependence is opposite for very shallow penetrators where coplanar prograde orbits have the lowest surviving chance (or equivalently most vulnerable). We provide numerical fits to the disruption probability and energy gain at the the BH encounter as a function of the penetration depth. The latter can be simply rescaled in terms of binary masses, their initial separation and the binary-to-BH mass ratio to evaluate the ejection velocity of a binary members in various systems. We also investigate the disruption of coplanar, eccentric binaries by a MBH. It is shown that for highly eccentric binaries retrograde orbits have a significantly increased disruption probability and ejection velocities compared to the circular binaries.

Recent observations identify a valley in the radius distribution of small exoplanets, with planets in the range $1.5-2.0\,{\rm R}_{\oplus}$ significantly less common than somewhat smaller or larger planets. This valley may suggest a bimodal population of rocky planets that are either engulfed by massive gas envelopes that significantly enlarge their radius, or do not have detectable atmospheres at all. One explanation of such a bimodal distribution is atmospheric erosion by high-energy stellar photons. We investigate an alternative mechanism: the luminosity of the cooling rocky core, which can completely erode light envelopes while preserving heavy ones, produces a deficit of intermediate sized planets. We evolve planetary populations that are derived from observations using a simple analytical prescription, accounting self-consistently for envelope accretion, cooling and mass loss, and demonstrate that core-powered mass loss naturally reproduces the observed radius distribution, regardless of the high-energy incident flux. Observations of planets around different stellar types may distinguish between photoevaporation, which is powered by the high-energy tail of the stellar radiation, and core-powered mass loss, which depends on the bolometric flux through the planet's equilibrium temperature that sets both its cooling and mass-loss rates.

The Galactic Centre is a highly crowded stellar field and frequent unrecognized events of source confusion, which involve undetected faint stars, are expected to introduce astrometric noise on a sub-mas level. This confusion noise is the main non-instrumental effect limiting the astrometric accuracy and precision of current near-infrared imaging observations and the long-term monitoring of individual stellar orbits in the vicinity of the central supermassive black hole. We self-consistently simulate the motions of the known and the yet unidentified stars to characterize this noise component and show that a likely consequence of source confusion is a bias in estimates of the stellar orbital elements, as well as the inferred mass and distance of the black hole, in particular if stars are being observed at small projected separations from it, such as the star S2 during pericentre passage. Furthermore, we investigate modeling the effect of source confusion as an additional noise component that is time-correlated, demonstrating a need for improved noise models to obtain trustworthy estimates of the parameters of interest (and their uncertainties) in future astrometric studies.

The relaxed distribution of stars around a massive black hole is known to follow a cusp profile $\rho(r)\propto r^{-\alpha}$ with characteristic slope $\alpha=7/4$. This follows from energy conservation and a scattering rate as given by two body encounters. However, we show that injection of stars close to the black hole, i.e. a source term in the standard cusp picture, modifies this profile. In the steady-state configuration, the cusp develops a central region with typical slope $\alpha=9/4$ in which stars diffuse outward. Binary disruption by the intense tidal field of the massive black hole is among the phenomena that take place in the Galactic Center. In such disruption, one of the binary members remains bound to the black hole, thus providing a source term of stars close to the black hole. Assuming a binary fraction of $0.1$ and an orbital circularization efficiency of $0.35$, we show that this source is strong enough to modify the cusp profile within $\approx 0.07$ pc in the Galactic Center. If the binary fraction at the influence radius is of order unity and the orbits of all the captured stars are efficiently circularized, the steeper cusp extends almost as far as the radius of influence of the black hole.

The high planetary multiplicity revealed by Kepler implies that Transit Time Variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach to TTVs technique that allows expanding the TTVs catalog towards lower TTV amplitude, shorter orbital period, and shallower transit depth. In the Spectral Approach we assume that a sinusoidal TTV exists in the data and then calculate the improvement to $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 new periodic TTVs in Kepler data (an increase of ~2/3 over a previous TTV catalog); (b) Constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; (c) Identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs sample shows no deficit of short period or low amplitude transits, in contrast to previous surveys in which the detection schemes were significantly biased against such systems.

The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_\gamma$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions. We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reproduce the observed $Br_\gamma$ emission. In our model, G2 is a gas cloud that undergoes tidal disruption in a dilute ambient medium. We find that during pericenter passage, the efficient cooling of the cloud results in a vertical collapse, compressing the cloud by a factor of $\sim5000$. By properly taking into account the ionization state of the gas, we find that the cloud is UV starved and are able to reproduce the observed $Br_\gamma$ luminosity. For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the cloud fragments, due to cooling instabilities and the emitted radiation is inconsistent with observations. For lower densities, the cloud survives the pericenter passage intact and its emitted radiation matches the observed lightcurve. From the duration of $Br_\gamma$ emission which contains both redshifted and blueshifted components, we show that the cloud is not spherical but rather elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation grows as it travels towards pericenter and is consistent with observations, due to viewing angles. The simulation is also consistent with having a spherical shape at apocenter.

We present detailed multi-frequency, multi-epoch radio observations of GRB 140304A at $z=5.283$ from 1 to 86 GHz and 0.45 d to 89 d. The radio and mm data exhibit unusual multiple spectral components, which cannot be simply explained by standard forward and reverse shock scenarios. Through detailed multi-wavelength analysis spanning radio to X-rays, we constrain the forward shock parameters to $E_{\rm K, iso}\approx4.9\times10^{54}\,$erg, $A_* \approx 2.6\times10^{-2}$, $\epsilon_{\rm e}\approx2.5\times10^{-2}$, $\epsilon_{\rm B}\approx5.9\times10^{-2}$, $p\approx2.6$, and $\theta_{\rm jet}\approx1.1^{\circ}$, yielding a beaming corrected $\gamma$-ray and kinetic energy, $E_{\gamma}\approx2.3\times10^{49}\,$erg and $E_{\rm K}\approx9.5\times10^{50}\,$erg, respectively. We model the excess radio emission as due to a combination of a late-time reverse shock (RS) launched by a shell collision, which also produces a re-brightening in the X-rays at $\approx0.26\,$d, and either a standard RS or diffractive interstellar scintillation. Under the standard RS interpretation, we invoke consistency arguments between the forward and reverse shocks to derive a deceleration time, $t_{\rm dec}\approx100\,$s, the ejecta Lorentz factor, $\Gamma(t_{\rm dec})\approx300$, and a low RS magnetization, $R_{\rm B}\approx0.6$. Our observations highlight both the power of radio observations in capturing RS emission and thus constraining the properties of GRB ejecta and central engines, and the challenge presented by interstellar scintillation in conclusively identifying RS emission in GRB radio afterglows.

The orbits of many observed hot Jupiters are decaying rapidly due to tidal interaction, eventually reaching the Roche limit. We analytically study the ensuing coupled mass loss and orbital evolution during the Roche-lobe overflow and find two possible scenarios. Planets with light cores $M_c\lesssim 6M_\oplus$ (assuming a nominal tidal dissipation factor $Q\sim 10^6$ for the host star) are transformed into Neptune-mass gas planets, orbiting at a separation (relative to the stellar radius) $a/R_\star\approx 3.5$. Planets with heavier cores $M_c\gtrsim 6M_\oplus$ plunge rapidly until they are destroyed at the stellar surface. Remnant gas-Neptunes, which are stable to photo-evaporation, are absent from the observations, despite their unique transit radius ($5-10R_\oplus$). This result suggests that $M_c\gtrsim 6M_\oplus$, providing a useful constraint on the poorly-known core mass that may distinguish between different formation theories of gas giants. Alternatively, if one assumes a prior of $M_c\approx 6 M_\oplus$ from the core-accretion theory, our results suggest that $Q$ does not lie in the range $10^6\lesssim Q\lesssim 10^7$.

We consider two body relaxation in a spherical system with a loss cone. Considering two-dimensional angular momentum space, we focus on "empty loss cone" systems, where the typical scattering during a dynamical time $j_{d}$ is smaller than the size of the loss cone $j_{\rm lc}$. As a result, the occupation number within the loss cone is significantly smaller than outside. Classical diffusive treatment of this regime predict exponentially small occupation number deep in the loss cone. We revisit this classical derivation of occupancy distribution of objects in the empty loss cone regime. We emphasize the role of the rare large scatterings and show that the occupancy does not decay exponentially within the loss cone, but it is rather flat, with a typical value $\sim [(j_d/j_{\rm lc})]^2\ln^{-2}(j_{\rm lc}/j_{\min})$ compared to the occupation in circular angular momentum (where $j_{\min}$ is the smallest possible scattering). Implication are that although the loss cone for tidal break of Giants or binaries is typically empty, tidal events which occurs significantly inside the loss cone ($\beta\gtrsim 2$), are almost as common as those with $\beta\cong 1$ where $\beta$ is the ratio between the tidal radius and the periastron. The probability for event with penetration factor $>\beta$ decreases only as $\beta^{-1}$ rather than exponentially. This effect has no influence on events characterized by full loss cone, such as tidal disruption event of $\sim 1m_\odot$ main sequence star.

We consider stable mass transfer from the secondary to the primary of an extreme mass ratio binary system. We show that when the mass transfer is sufficiently fast, mass leakage occurs through the outer Lagrange point L2, in addition to the usual transfer through L1. We provide an analytical estimate for the mass leakage rate through L2 and find the conditions in which it is comparable to the mass transfer rate through L1. Focusing on a binary system of a main-sequence star and a super-massive black hole, driven by the emission of gravitational radiation, we show that it may sustain stable mass transfer, along with mass loss through L2. If such a mass-transferring system occurs at our Galactic Centre, it produces a gravitational wave signal detectable by future detectors, such as eLISA. The signal evolves according to the star's adiabatic index and cooling time. For low mass stars, the evolution is faster than the Kelvin-Helmholtz cooling rate driving the star out of the main-sequence. In some cases, the frequency and amplitude of the signal may both decrease with time, contrary to the standard chirp of a coalescing binary. Mass loss through L2, when occurs, decreases the evolution timescale of the emitted gravitational wave signal by up to a few tens of per cent. We conclude that L2 mass ejection is a crucial factor in analyzing gravitational waves signals produced by such systems.

The mass assembly history of the Milky Way can inform both theory of galaxy formation and the underlying cosmological model. Thus, observational constraints on the properties of both its baryonic and dark matter contents are sought. Here we show that hypervelocity stars (HVSs) can in principle provide such constraints. We model the observed velocity distribution of HVSs, produced by tidal break-up of stellar binaries caused by Sgr A*. Considering a Galactic Centre (GC) binary population consistent with that inferred in more observationally accessible regions, a fit to current HVS data with significance level > 5\% can only be obtained if the escape velocity from the GC to 50 kpc is $V_G < 850$ km/s, regardless of the enclosed mass distribution. When a NFW matter density profile for the dark matter halo is assumed, haloes with $V_G < 850$ km/s are in agreement with predictions in the $\Lambda$CDM model and that a subset of models around $M_{200} \sim 0.5-1.5 \times 10^{12}$ solar masses and $r_s < 35$ kpc can also reproduce Galactic circular velocity data. HVS data alone cannot currently exclude potentials with $V_G > 850$ km/s. Finally, specific constraints on the halo mass from HVS data are highly dependent on the assumed baryonic mass potentials. This first attempt to simultaneously constrain GC and dark halo properties is primarily hampered by the paucity and quality of data. It nevertheless demonstrates the potential of our method, that may be fully realised with the ESA Gaia mission.

Macronovae (kilonovae) that arise in binary neutron star mergers are powered by radioactive beta decay of hundreds of $r$-process nuclides. We derive, using Fermi's theory of beta decay, an analytic estimate of the nuclear heating rate. We show that the heating rate evolves as a power law ranging between $t^{-6/5}$ to $t^{-4/3}$. The overall magnitude of the heating rate is determined by the mean values of nuclear quantities, e.g., the nuclear matrix elements of beta decay. These values are specified by using nuclear experimental data. We discuss the role of higher order beta transitions and the robustness of the power law. The robust and simple form of the heating rate suggests that observations of the late-time bolometric light curve $\propto t^{-\frac{4}{3}}$ would be a direct evidence of a $r$-process driven macronova. Such observations could also enable us to estimate the total amount of $r$-process nuclei produced in the merger.

Using 25 years of data from uninterrupted monitoring of stellar orbits in the Galactic Center, we present an update of the main results from this unique data set: A measurement of mass of and distance to SgrA*. Our progress is not only due to the eight year increase in time base, but also due to the improved definition of the coordinate system. The star S2 continues to yield the best constraints on the mass of and distance to SgrA*; the statistical errors of 0.13 x 10^6 M_sun and 0.12 kpc have halved compared to the previous study. The S2 orbit fit is robust and does not need any prior information. Using coordinate system priors, also the star S1 yields tight constraints on mass and distance. For a combined orbit fit, we use 17 stars, which yields our current best estimates for mass and distance: M = 4.28 +/- 0.10|stat. +/. 0.21|sys. x 10^6 M_sun and R_0 = 8.32 +/- 0.07|stat. +/- 0.14|sys. kpc. These numbers are in agreement with the recent determination of R_0 from the statistical cluster parallax. The positions of the mass, of the near-infrared flares from SgrA* and of the radio source SgrA* agree to within 1mas. In total, we have determined orbits for 40 stars so far, a sample which consists of 32 stars with randomly oriented orbits and a thermal eccentricity distribution, plus eight stars for which we can explicitly show that they are members of the clockwise disk of young stars, and which have lower eccentricity orbits.

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.

Supernovae near the galactic center evolve differently from regular galactic supernovae. This is mainly due to the environment into which the supernova remnants propagate. Instead of a static, uniform density medium, SNRs near the galactic center propagate into a wind swept environment with a velocity away from the galactic center, and a graded density profile. This causes these SNRs to be non - spherical, and to evolve faster than their galactic counterparts. We develop an analytic theory for the evolution of explosions within a stellar wind, and verify it using a hydrodynamic code. We show that such explosions can evolve in one of three possible morphologies. Using these results we discuss the association between the two SNRs (SGR East and SGR A's bipolar radio/X-ray Lobes) and the two neutron stars (the cannonball and SGR J1745-2900) near the galactic center. We show that, given the morphologies of the SNR and positions of the neutron stars, the only possible association is between SGR A's bipolar radio/X-ray Lobes and SGR J1745-2900. The compact object created in the explosion of SGR East remains undetected, and the SNR of the supernova that created the cannonball has already disappeared.

We investigate the relativistic breakout of a shock wave from the surface of a star. In this process, each fluid shell is endowed with some kinetic and thermal energy by the shock, and then continues to accelerate adiabatically by converting thermal energy into kinetic energy. This problem has been previously studied for a mildly relativistic breakout, where the acceleration ends close to the surface of the star. The current work focuses on the case where the acceleration ends at distances much greater than the radius of the star. We derive an analytic description for the hydrodynamic evolution of the ejecta in this regime, and validate it using a numerical simulation. We also provide predictions for the expected light curves and spectra from such an explosion. The relevance to astrophysical explosions is discussed, and it is shown that such events require more energy than is currently believed to result from astrophysical explosions.

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.

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.

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.

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.

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.

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.

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.

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 $\tau_{\rm dep}$, creates an exterior convective region, between optical depths $L_{\rm eq}/L_{\rm dep}$ and $\tau_{\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}\tau_{\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}\tau_{\rm dep}/L_{\rm eq}\right)^\beta$, where $\beta\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\tau^\beta$. 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.

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.

The Asteroid Belt and the Kuiper Belt are relics from the formation of our solar system. Understanding the size and spin distribution of the two belts is crucial for a deeper understanding of the formation of our solar system and the dynamical process that govern it. In this paper, we investigate the effect of collisions on the evolution of the spin distribution of asteroids and KBO's. We find that the power law nature of the impactors' size distribution leads to a L\'evy distribution of the spin rates. This results in a power law tail of the spin distribution, in stark contrast to the usually quoted Maxwellian distribution. We show that for bodies larger than 10 km, collisions alone lead to spin rates peaking at 0.15-0.5 revolutions per day. Comparing that to the observed spin rates of large asteroids ($R>50$ km), we find that the spins of large asteroids, peaking at $\sim1-2$ revolutions per day, are dominated by a primordial component that reflects the formation mechanism of the asteroids. Similarly, the Kuiper Belt has undergone virtually no collisional spin evolution, assuming current day density. Collisions contribute a spin rate of $\sim0.01$ revolutions per day, thus the observed fast spin rates of KBOs are also primordial in nature.

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.

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.

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.

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.

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.

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 ($\gamma=5/3$) that the local atmospheric mass loss fraction for ground velocities $v_g < 0.25 v_{esc}$ is given by $\chi_{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} \rho_0 (\pi 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 $\rho_0$ are the atmospheric scale height, radius of the target, and its atmospheric density at the ground. 2) Smaller impactors, but above $m>4 \pi \rho_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)

We present new solutions to the strong explosion problem in a non power law density profi{}le. The unperturbed self similar solutions developed by Sedov, Taylor and Von Neumann describe strong Newtonian shocks propagating into a cold gas with a density profile falling off as $r^{-\omega}$, where $\omega\le\frac{7-\gamma}{\gamma+1}$ (filled type I solutions), and $\gamma$ is the adiabatic index of the gas. The perturbations we consider are spherically symmetric and log periodic with respect to the radius. While the unperturbed solutions are continuously self similar, the log periodicity of the density perturbations leads to a discrete self similarity of the perturbations, i.e., the solution repeats itself up to a scaling at discrete time intervals. We discuss these solutions and verify them against numerical integrations of the time dependent hydrodynamic equations. This is an extension of a previous investigation on type II solutions and helps clarifying boundary conditions for perturbations to type I self similar solutions.

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 $\rho\left(r,\theta,\phi\right)=kr^{-\omega}\left(1+\varepsilon\left(\frac{r}{r_{0}}\right)^{q}Y_{lm}\left(\theta,\phi\right)\right)$, where $\varepsilon\ll1$ and $\omega\le\frac{7-\gamma}{\gamma+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

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

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.

We present continued multi-frequency radio observations of the relativistic tidal disruption event Sw1644+57 extending to dt~600 d. The data were obtained with the JVLA and AMI Large Array. We combine these data with public Swift/XRT and Chandra X-ray observations over the same time-frame to show that the jet has undergone a dramatic transition starting at ~500 d, with a sharp decline in the X-ray flux by about a factor of 170 on a timescale of dt/t<0.2. The rapid decline rules out a forward shock origin (direct or reprocessing) for the X-ray emission at <500 d, and instead points to internal dissipation in the inner jet. On the other hand, our radio data uniquely demonstrate that the low X-ray flux measured by Chandra at ~610 d is consistent with emission from the forward shock. Furthermore, the Chandra data are inconsistent with thermal emission from the accretion disk itself since the expected temperature of 30-60 eV and inner radius of 2-10 R_s cannot accommodate the observed flux level or the detected emission at >1 keV. We associate the rapid decline with a turn off of the relativistic jet when the mass accretion rate dropped below Mdot_Edd~0.006 Msun/yr (for a 3x10^6 Msun black hole and order unity efficiency) indicating that the peak accretion rate was about 330 Mdot_Edd, and the total accreted mass by 500 d is about 0.15 Msun. From the radio data we further find significant flattening in the integrated energy of the forward shock at >250 d with E_j,iso~2x10^54 erg (E_j~10^52$ erg for a jet opening angle, theta_j=0.1) following a rise by about a factor of 15 at 30-250 d. Projecting forward, we predict that the emission in the radio and X-ray bands will evolve in tandem with similar decline rates.

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 $\gamma$-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.

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)

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.

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.

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.

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.

We present the discovery of PTF 10vgv, a Type Ic supernova detected by the Palomar Transient Factory, using the Palomar 48-inch telescope (P48). R-band observations of the PTF 10vgv field with P48 probe the supernova emission from its very early phases (about two weeks before R-band maximum), and set limits on its flux in the week prior to the discovery. Our sensitive upper limits and early detections constrain the post-shock-breakout luminosity of this event. Via comparison to numerical (analytical) models, we derive an upper-limit of R \lesssim 4.5 Rsun (R \lesssim 1 Rsun) on the radius of the progenitor star, a direct indication in favor of a compact Wolf-Rayet star. Applying a similar analysis to the historical observations of SN 1994I, yields R \lesssim 1/4 Rsun for the progenitor radius of this supernova.

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.

According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the \gamma-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.

We present continued radio observations of the tidal disruption event SwiftJ164449.3+573451 extending to \sim216 days after discovery. The data are part of a long-term program to monitor the expansion and energy scale of the relativistic outflow, and to trace the parsec-scale environment around a previously-dormant supermassive black hole (SMBH). The new observations reveal a significant change in the radio evolution starting at \sim1 month, with a brightening at all frequencies that requires an increase in the energy by about an order of magnitude, and an overall density profile around the SMBH of rho \propto r^{-3/2} (0.1-1.2 pc) with a significant flattening at r\sim0.4-0.6 pc. The increase in energy cannot be explained with continuous injection from an L \propto t^{-5/3} tail, which is observed in the X-rays. Instead, we conclude that the relativistic jet was launched with a wide range of Lorentz factors, obeying E(>Gamma) \propto Gamma^{-2.5}. The similar ratio of duration to dynamical timescale for Sw1644+57 and GRBs suggests that this result may be applicable to GRBs as well. The radial density profile may be indicative of Bondi accretion, with the inferred flattening at r\sim0.5 pc in good agreement with the Bondi radius for a \sim10^6 M_sun black hole. The density at \sim0.5 pc is about a factor of 30 times lower than inferred for the Milky Way galactic center, potentially due to a smaller number of mass-shedding massive stars. From our latest observations (\sim216 d) we find that the jet energy is E_{iso}\sim5x10^{53} erg (E_j\sim2.4x10^{51} erg for theta_j=0.1), the radius is r\sim1.2 pc, the Lorentz factor is Gamma\sim2.2, the ambient density is n\sim0.2 cm^{-3}, and the projected size is r_{proj}\sim25 microarcsec. Assuming no future changes in the observed evolution we predict that the radio emission from Sw1644+57 should be detectable with the EVLA for several decades, and will be resolvable with VLBI in a few years.

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.

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.

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.

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.

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.

Runway growth is an important stage in planet formation during which large protoplanets form, while most of the initial mass remains in small planetesimals. The amount of mass converted into large protoplanets and their resulting size distribution are not well understood. Here, we use analytic work, that we confirm by coagulation simulations, to describe runaway growth and the corresponding evolution of the velocity dispersion. We find that runaway growth proceeds as follows: Initially all the mass resides in small planetesimals, with mass surface density \sigma, and large protoplanets start to form by accreting small planetesimals. This growth continues until growth by merging large protoplanets becomes comparable to growth by planetesimal accretion. This condition sets in when \Sigma/\sigma ~\alpha^{3/4} ~ 10^{-3}, where \Sigma is the mass surface density in protoplanets in a given logarithmic mass interval and \alpha is the ratio of the size of a body to its Hill radius. From then on, protoplanetary growth and the evolution of the velocity dispersion become self-similar and \Sigma remains roughly constant, since an increase in \Sigma by accretion of small planetesimals is balanced by a decrease due to merging with large protoplanets. We show that this growth leads to a protoplanet size distribution given by N(>R) \propto R^{-3} where N(>R) is the number of objects with radii greater than R (i.e., a differential power-law index of 4). We apply our results to the Kuiper Belt, which is a relic of runaway growth. Our results successfully match the observed Kuiper belt size distribution, they illuminate the physical processes that shaped it and explain the total mass that is present in large Kuiper belt objects (KBOs) today. This work suggests that the current mass in large KBOs is primordial and that it has not been significantly depleted. Abridged

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.

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.

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.

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.

The Kuiper belt is a remnant of the primordial Solar System. Measurements of its size distribution constrain its accretion and collisional history, and the importance of material strength of Kuiper belt objects (KBOs). Small, sub-km sized, KBOs elude direct detection, but the signature of their occultations of background stars should be detectable. Observations at both optical and X-ray wavelengths claim to have detected such occultations, but their implied KBO abundances are inconsistent with each other and far exceed theoretical expectations. Here, we report an analysis of archival data that reveals an occultation by a body with a 500 m radius at a distance of 45 AU. The probability of this event to occur due to random statistical fluctuations within our data set is about 2%. Our survey yields a surface density of KBOs with radii larger than 250 m of 2.1^{+4.8}_{-1.7} x 10^7 deg^{-2}, ruling out inferred surface densities from previous claimed detections by more than 5 sigma. The fact that we detected only one event, firmly shows a deficit of sub-km sized KBOs compared to a population extrapolated from objects with r>50 km. This implies that sub-km sized KBOs are undergoing collisional erosion, just like debris disks observed around other stars.

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]

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.

We present analytic approximations to the optically thin synchrotron and synchrotron self-Compton (SSC) spectra when Klein-Nishina (KN) effects are important and pair production and external radiation fields can be neglected. This theory is useful for analytical treatment of radiation from astrophysical sources, such as gamma-ray bursts (GRBs), active galactic nuclei and pulsar wind nebula, where KN effects may be important. We consider a source with a continuous injection of relativistic electrons with a power-law energy distribution above some typical injection energy. We find that the synchrotron-SSC spectra can be described by a broken power-law, and provide analytic estimates for the break frequencies and power-law indices. In general, we show that the dependence of the KN cross-section on the energy of the upscattering electron results in a hardening of the energy distribution of fast cooling electrons and therefore in a hardening of the observed synchrotron spectrum. As a result the synchrotron spectrum of fast cooling electrons, below the typical injection energy, can be as hard as $F_\nu \propto \nu^0$, instead of the classical $\nu^{-1/2}$ when KN effects are neglected. The synchrotron energy output can be dominated by electrons with energy above the typical injection energy. We solve self-consistently for the cooling frequency and find that the transition between synchrotron and SSC cooling can result in a discontinuous variations of the cooling frequency and the synchrotron and SSC spectra. We demonstrate the application of our results to theory by applying them to prompt and afterglow emission models of GRBs.

We present a simple theoretical framework for massive galaxies at high redshift, where the main assembly and star formation occurred, and report on the first cosmological simulations that reveal clumpy disks consistent with our analysis. The evolution is governed by the interplay between smooth and clumpy cold streams, disk instability, and bulge formation. Intense, relatively smooth streams maintain an unstable dense gas-rich disk. Instability with high turbulence and giant clumps, each a few percent of the disk mass, is self-regulated by gravitational interactions within the disk. The clumps migrate into a bulge in ~10 dynamical times, or ~0.5Gyr. The cosmological streams replenish the draining disk and prolong the clumpy phase to several Gigayears in a steady state, with comparable masses in disk, bulge, and dark matter within the disk radius. The clumps form stars in dense subclumps following the overall accretion rate, ~100 Msun/yr, and each clump converts into stars in ~0.5 Gyr. While the clumps coalesce dissipatively to a compact bulge, the star-forming disk is extended because the incoming streams keep the outer disk dense and susceptible to instability and because of angular momentum transport. Passive spheroid-dominated galaxies form when the streams are more clumpy: the external clumps merge into a massive bulge and stir up disk turbulence that stabilize the disk and suppress in situ clump and star formation. We predict a bimodality in galaxy type by z~3, involving giant-clump star-forming disks and spheroid-dominated galaxies of suppressed star formation. After z~1, the disks tend to be stabilized by the dominant stellar disks and bulges. Most of the high-z massive disks are likely to end up as today's early-type galaxies.

We present new solutions to the strong explosion problem in a non-power law density profile. The unperturbed self-similar solutions discovered by Waxman & Shvarts describe strong Newtonian shocks propagating into a cold gas with a density profile falling off as $r^{-\omega}$, where $\omega>3$ (Type-II solutions). The perturbations we consider are spherically symmetric and log-periodic with respect to the radius. While the unperturbed solutions are continuously self-similar, the log-periodicity of the density perturbations leads to a discrete self-similarity of the perturbations, i.e. the solution repeats itself up to a scaling at discrete time intervals. We discuss these solutions and verify them against numerical integrations of the time dependent hydrodynamic equations. Finally we show that this method can be generalized to treat any small, spherically symmetric density perturbation by employing Fourier decomposition.

Recently, the first collisional family was discovered in the Kuiper belt. The parent body of this family, Haumea, is one of the largest objects in the Kuiper belt and is orbited by two satellites. It has been proposed that the Haumea family was created from dispersed fragments that resulted from a giant impact. This proposed origin of the Haumea family is however in conflict with the observed velocity dispersion between the family members (\sim 140 m/s) which is significantly less than the escape velocity from Haumea's surface (\sim 900 m/s). In this paper we propose a different formation scenario for Haumea's collisional family. In our scenario the family members are ejected while in orbit around Haumea. This scenario, therefore, gives naturally rise to a lower velocity dispersion among the family members than expected from direct ejection from Haumea's surface. In our scenario Haumea's giant impact forms a single moon that tidally evolves outward until it suffers a destructive collision from which the family is created. We show that this formation scenario yields a velocity dispersion of \sim 190m/s among the family members which is in good agreement with the observations. The probability for Haumea's initial giant impact in todays Kuiper belt is less than 10^{-3}. In our scenario, however, Haumea's giant impact can occur before the excitation of the Kuiper belt and the ejection of the family members afterwards. This has the advantage that one can preserve the dynamical coherence of the family and explain Haumea's original giant impact, which is several orders of magnitude more likely to have occurred in the primordial dynamically cold Kuiper belt compared to the dynamically excited Kuiper belt today. Abridged

The detection of GeV photons from gamma-ray bursts (GRBs) has important consequences for the interpretation and modelling of these most-energetic cosmological explosions. The full exploitation of the high-energy measurements relies, however, on the accurate knowledge of the distance to the events. Here we report on the discovery of the afterglow and subsequent redshift determination of GRB 080916C, the first GRB detected by the Fermi Gamma-Ray Space Telescope with high significance detection of photons at >0.1 GeV. Observations were done with 7-channel imager GROND at the 2.2m MPI/ESO telescope, the SIRIUS instrument at the Nagoya-SAAO 1.4m telescope in South Africa, and the GMOS instrument at Gemini-S. The afterglow photometric redshift of z=4.35+-0.15, based on simultaneous 7-filter observations with the Gamma-Ray Optical and Near-infrared Detector (GROND), places GRB 080916C among the top 5% most distant GRBs, and makes it the most energetic GRB known to date. The detection of GeV photons from such a distant event is rather surprising. The observed gamma-ray variability in the prompt emission together with the redshift suggests a lower limit for the Lorentz factor of the ultra-relativistic ejecta of Gamma > 1090. This value rivals any previous measurements of Gamma in GRBs and strengthens the extreme nature of GRB 080916C.

Inverse Compton (IC) scattering is one of two viable mechanisms that can produce the prompt gamma-ray emission in Gamma-Ray Bursts. IC requires low energy seed photons and a population of relativistic electrons that upscatter them. The same electrons can upscatter the gamma-ray photons to even higher energies in the TeV range. Using the current upper limits on the prompt optical emission we show here that under general conservative assumption the IC mechanism suffers from an "energy crisis". Namely, IC will over-produce a very high energy component that would carry much more energy than the observed prompt gamma-rays. Our analysis is general and it makes no assumptions on the specific mechanism that produces the relativistic electrons population.

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, $\Sigma/\sigma$. Early in the oligarchic phase, when $\Sigma/\sigma$ 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.

The extremely bright optical flash that accompanied GRB 080319B suggested, at first glance, that the prompt $\gamma$-rays in this burst were produced by Synchrotron self Compton (SSC). We analyze here the observed optical and $\gamma$ spectrum. We find that the very strong optical emission poses, due to self absorption, very strong constraints on the emission processes and put the origin of the optical emission at a very large radius, almost inconsistent with internal shock. Alternatively it requires a very large random Lorentz factor for the electrons. We find that SSC could not have produced the prompt $\gamma$-rays. We also show that the optical emission and the $\gamma$ rays could not have been produced by synchrotron emission from two populations of electron within the same emitting region. Thus we must conclude that the optical and the $\gamma$-rays were produced in different physical regions. A possible interpretation of the observations is that the $\gamma$-rays arose from internal shocks but the optical flash resulted from external shock emission. This would have been consistent with the few seconds delay observed between the optical and $\gamma$-rays signals.

We examine the evolution of an almost circular Keplerian orbit interacting with unbound perturbers. We calculate the change in eccentricity and angular momentum that results from a single encounter, assuming the timescale for the interaction is shorter than the orbital period. The orbital perturbations are incorporated into a Boltzmann equation that allows for eccentricity dissipation. We present an analytic solution to the Boltzmann equation that describes the distribution of orbital eccentricity and relative inclination as a function of time. The eccentricity and inclination of the binary do not evolve according to a normal random walk but perform a Levy flight. The slope of the mass spectrum of perturbers dictates whether close gravitational scatterings are more important than distant tidal ones. When close scatterings are important, the mass spectrum sets the slope of the eccentricity and inclination distribution functions. We use this general framework to understand the eccentricities of several Kuiper belt systems: Pluto, 2003 EL 61, and Eris. We use the model of Tholen et al (2007) to separate the non-Keplerian components of the orbits of Pluto's outer moons Nix and Hydra from the motion excited by interactions with other Kuiper belt objects. Our distribution is consistent with the observations of Nix, Hydra, and the satellites of 2003 EL 61 and Eris. We address applications of this work to objects outside of the solar system, such as extrasolar planets around their stars and millisecond pulsars.

We investigate the GeV emission from gamma-ray bursts (GRBs), using the results from the Energetic Gamma Ray Experimental Telescope (EGRET), and in view of the Gamma-ray Large Area Space Telescope (GLAST). Assuming that the conventional prompt and afterglow photons originate from synchrotron radiation, we compare an accompanying inverse-Compton component with EGRET measurements and upper limits on GeV fluence, taking Klein-Nishina feedback into account. We find that EGRET constraints are consistent with the theoretical framework of the synchrotron self-Compton model for both prompt and afterglow phases, and discuss constraints on microphysical parameters in both phases. Based on the inverse-Compton model and using EGRET results, we predict that GLAST would detect GRBs with GeV photons at a rate >~20 yr^{-1} from each of the prompt and afterglow phases. This rate applies to the high-energy tail of the prompt synchrotron emission and to the inverse-Compton component of the afterglow. Theory predicts that in a large fraction of the cases where synchrotron GeV prompt emission would be detected by GLAST, inverse-Compton photons should be detected as well at high energies >~10 GeV. Therefore GLAST will enable a more precise test of the high-energy emission mechanism. Finally, we show that the contribution of GRBs to the flux of the extragalactic gamma-ray background measured with EGRET is at least 0.01% and likely around 0.1%.

The only supernovae (SNe) to have shown early gamma-ray or X-ray emission thus far are overenergetic, broad-lined Type Ic SNe (Hypernovae - HNe). Recently, SN 2008D shows several novel features: (i) weak XRF, (ii) an early, narrow optical peak, (iii) disappearance of the broad lines typical of SNIc HNe, (iv) development of He lines as in SNeIb. Detailed analysis shows that SN 2008D was not a normal SN: its explosion energy (KE ~ 6*10^{51} erg) and ejected mass (~7 Msun) are intermediate between normal SNeIbc and HNe. We derive that SN 2008D was originally a ~30Msun star. When it collapsed a black hole formed and a weak, mildly relativistic jet was produced, which caused the XRF. SN 2008D is probably among the weakest explosions that produce relativistic jets. Inner engine activity appears to be present whenever massive stars collapse to black holes.

With the discovery of Kuiper Belt binaries that have wide separations and roughly equal masses new theories were proposed to explain their formation. Two formation scenarios were suggested by Goldreich and collaborators: In the first, dynamical friction that is generated by a sea of small bodies enables a transient binary to become bound ($L^2s$ mechanism); in the second, a transient binary gets bound by an encounter with a third body ($L^3$ mechanism). We show that these different binary formation scenarios leave their own unique signatures in the relative abundance of prograde to retrograde binary orbits. This signature is due to stable retrograde orbits that exist much further out in the Hill sphere than prograde orbits. It provides an excellent opportunity to distinguish between the different binary formation scenarios observationally. We predict that if binary formation proceeded while sub-Hill velocities prevailed, the vast majority of all comparable mass ratio binaries have retrograde orbits. This dominance of retrograde binary orbits is a result of binary formation via the $L^2s$ mechanism, or any other mechanism that dissipates energy in a smooth and gradual manner. For super-Hill velocities binary formation proceeds via the $L^3$ mechanism which produces a roughly equal number of prograde and retrograde binaries. These predictions assume that subsequent orbital evolution due to dynamical friction and dynamical stirring of the Kuiper belt did not alter the sense of the binary orbit after formation.

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.

Many small bodies in the solar system are believed to be rubble piles, a collection of smaller elements separated by voids. We propose a model for the structure of a self-gravitating rubble pile. Static friction prevents its elements from sliding relative to each other. Stresses are concentrated around points of contact between individual elements. The effective dimensionless rigidity, $\tilde\mu_{rubble}$, is related to that of a monolithic body of similar composition and size, $\tilde\mu$ by $\tilde \mu_{rubble} \sim \tilde \mu^{1/2} \epsilon_Y^{-1/2}$, where $\epsilon_Y \sim 10^{-2}$ is the yield strain. This represents a reduction in effective rigidity below the maximum radius, $R_{max}\sim [\mu\epsilon_Y/(G\rho^2)]^{1/2}\sim 10^3\km$, at which a rubble pile can exist. Densities derived for binary near-Earth asteroids imply that they are rubble piles. As a consequence, their tidal evolution proceeds $10^3$ to $10^4$ times faster than it would if they were monoliths. This accounts for both the sizes of their semimajor axes and their small orbital eccentricities. We show that our model for the rigidity of rubble piles is compatible with laboratory experiment in sand.

The discovery that a substantial fraction of Kuiper Belt objects (KBOs) exists in binaries with wide separations and roughly equal masses, has motivated a variety of new theories explaining their formation. Goldreich et al. (2002) proposed two formation scenarios: In the first, a transient binary is formed, which becomes bound with the aid of dynamical friction from the sea of small bodies (L^2s mechanism); in the second, a binary is formed by three body gravitational deflection (L^3 mechanism). Here, we accurately calculate the L^2s and L^3 formation rates for sub-Hill velocities. While the L^2s formation rate is close to previous order of magnitude estimates, the L^3 formation rate is about a factor of 4 smaller. For sub-Hill KBO velocities (v << v_H) the ratio of the L^3 to the L^2s formation rate is 0.05 (v/v_H) independent of the small bodies' velocity dispersion, their surface density or their mutual collisions. For Super-Hill velocities (v >> v_H) the L^3 mechanism dominates over the L^2s mechanism. Binary formation via the L^3 mechanism competes with binary destruction by passing bodies. Given sufficient time, a statistical equilibrium abundance of binaries forms. We show that the frequency of long-lived transient binaries drops exponentially with the system's lifetime and that such transient binaries are not important for binary formation via the L^3 mechanism, contrary to Lee et al. (2007). For the L^2s mechanism we find that the typical time, transient binaries must last, to form Kuiper Belt binaries (KBBs) for a given strength of dynamical friction, D, increases only logarithmically with D. Longevity of transient binaries only becomes important for very weak dynamical friction (i.e. D \lesssim 0.002) and is most likely not crucial for KBB formation.

On 2006 May 5, a four second duration, low energy, ~10^49 erg, Gamma-Ray Burst (GRB) was observed, spatially associated with a z=0.0894 galaxy. Here, we report the discovery of the GRB optical afterglow and observations of its environment using Gemini-south, Hubble Space Telescope (HST), Chandra, Swift and the Very Large Array. The optical afterglow of this GRB is spatially associated with a prominent star forming region in the Sc-type galaxy 2dFGRS S173Z112. Its proximity to a star forming region suggests that the progenitor delay time, from birth to explosion, is smaller than about 10 Myr. Our HST deep imaging rules out the presence of a supernova brighter than an absolute magnitude of about -11 (or -12.6 in case of ``maximal'' extinction) at about two weeks after the burst, and limits the ejected mass of radioactive Nickel 56 to be less than about 2x10^-4 solar mass (assuming no extinction). Although it was suggested that GRB 060505 may belong to a new class of long-duration GRBs with no supernova, we argue that the simplest interpretation is that the physical mechanism responsible for this burst is the same as for short-duration GRBs.

Planetesimal accretion during planet formation is usually treated as collisionless. Such accretion from a uniform and dynamically cold disk predicts protoplanets with slow retrograde rotation. However, if the building blocks of protoplanets, planetesimals, are small, of order of a meter in size, then they are likely to collide within the protoplanet's sphere of gravitational influence, creating a prograde accretion disk around the protoplanet. The accretion of such a disk results in the formation of protoplanets spinning in the prograde sense with the maximal spin rate allowed before centrifugal forces break them apart. As a result of semi-collisional accretion, the final spin of a planet after giant impacts is not completely random but is biased toward prograde rotation. The eventual accretion of the remaining planetesimals in the post giant-impact phase might again be in the semi-collisional regime and delivers a significant amount of additional prograde angular momentum to the terrestrial planets. We suggest that in our Solar System, semi-collisional accretion gave rise to the preference for prograde rotation observed in the terrestrial planets and perhaps the largest asteroids.

We examine the growth of eccentricities of a population of particles with initially circular orbits around a central massive body. Successive encounters between pairs of particles increase the eccentricities in the disk on average. As long as the epicyclic motions of the particles are small compared to the shearing motion between Keplerian orbits, there is no preferred scale for the eccentricities. The simplification due to this self-similarity allows us to find an analytic form for the distribution function; full numerical integrations of a disk with 200 planetesimals verify our analytical self-similar distribution. The shape of this non-equilibrium profile is identical to the equilibrium profile of a shear-dominated population whose mutual excitations are balanced by dynamical friction or Epstein gas drag.

Over the last decade, long-duration gamma-ray bursts (GRBs) including the subclass of X-ray flashes (XRFs) have been revealed to be a rare variety of Type Ibc supernova (SN). While all these events result from the death of massive stars, the electromagnetic luminosities of GRBs and XRFs exceed those of ordinary Type Ibc SNe by many orders of magnitude. The essential physical process that causes a dying star to produce a GRB or XRF, and not just an SN, remains the crucial open question. Here we present radio and X-ray observations of XRF 060218 (associated with SN 2006aj), the second nearest GRB identified to-date, which allow us to measure its total energy and place it in the larger context of cosmic explosions. We show that this event is 100 times less energetic but ten times more common than cosmological GRBs. Moreover, it is distinguished from ordinary Type Ibc SNe by the presence of 10^48 erg coupled to mildly-relativistic ejecta, along with a central engine (an accretion-fed, rapidly rotating compact source) which produces X-rays for weeks after the explosion. This suggests that the production of relativistic ejecta is the key physical distinction between GRBs/XRFs and ordinary SNe, while the nature of the central engine (black hole or magnetar) may distinguish typical bursts from low-luminosity, spherical events like XRF 060218.

Over the past decade our physical understanding of gamma-ray bursts (GRBs) has progressed rapidly thanks to the discovery and observation of their long-lived afterglow emission. Long-duration (T < 2 s) GRBs are associated with the explosive deaths of massive stars (``collapsars''), which produce accompanying supernovae, while the short-duration (T > 2 s) GRBs arise from a different origin, which has been argued to be the merger of two compact objects, either neutron stars or black holes. Here we present observations of GRB060614, a 100-s long burst discovered by the Swift satellite, which require the invocation of a new explosive process: either a massive ``collapsar'' that powers a GRB without any associated supernova, or a new type of engine, as long-lived as the collapsar but without any such massive stellar host. We also discuss the properties of this burst's redshift z=0.125 host galaxy, which distinguish it from other long-duration GRBs and suggest that an entirely new type of GRB progenitor may be required.

We present observations of the unusually bright and long gamma-ray burst GRB 050820A, one of the best-sampled broadband data sets in the Swift era. The gamma-ray light curve is marked by a soft precursor pulse some 200 s before the main event; the lack of any intervening emission suggests that it is due to a physical mechanism distinct from the GRB itself. The large time lag between the precursor and the main emission enabled simultaneous observations in the gamma-ray, X-ray, and optical band-passes, something only achieved for a handful of events to date. While the contemporaneous X-rays are the low-energy tail of the prompt emission, the optical does not directly track the gamma-ray flux. Instead, the early-time optical data appear mostly consistent with the forward shock synchrotron peak passing through the optical, and are therefore likely the beginning of the afterglow. On hour time scales after the burst, the X-ray and optical light curves are inconsistent with an adiabatic expansion of the shock into the surrounding region, but rather indicate that there is a period of energy injection. Observations at late times allow us to constrain the collimation angle of the relativistic outflow to theta = 6.8 - 9.3 degrees. Our estimates of both the kinetic energy of the afterglow and the prompt gamma-ray energy release make GRB 050820A one of the most energetic events for which such values could be determined.

On December 27, 2004, a giant flare from SGR 1806$-$20 was detected on earth. Its thermal spectrum and temperature suggest that the flare resulted from an energy release of about $10^{47}$ erg/sec close to the surface of a neutron star in the form of radiation and/or pairs. This plasma expanded under its own pressure producing a fireball and the observed gamma-rays escaped once the fireball became optically thin. The giant flare was followed by a bright radio afterglow, with an observable extended size, implying an energetic relativistic outflow. We revisit here the evolution of relativistic fireballs and we calculate the Lorentz factor and energy remaining in relativistic outflow once the radiation escapes. We show that pairs that arise naturally in a pure pairs-radiation fireball do not carry enough energy to account for the observed afterglow. We consider various alternatives and we show that if the relativistic outflow that causes the afterglow is related directly to the prompt flare, then the initial fireball must be loaded by baryons or Poynting flux. While we focus on parameters applicable to the giant flare and the radio afterglow of SGR 1806$-$20 the calculations presented here might be also applicable to GRBs.

The velocity dispersion, or eccentricity distribution, of protoplanets interacting with planetesimals is set by a balance between dynamical friction and viscous stirring. We calculate analytically the eccentricity distribution function of protoplanets embedded in a cold, shear-dominated planetesimal swarm. We find a distinctly non-Rayleigh distribution with a simple analytical form. The peak of the distribution lies much lower than the root-mean-squared value, indicating that while most of the bodies have similarly small eccentricities, a small subset of the population contains most of the thermal energy. We also measure the shear-dominated eccentricity distribution using numerical simulations. The numerical code treats each protoplanet explicitly and adds an additional force term to each body to represent the dynamical friction of the planetesimals. Without fitting any parameters, the eccentricity distribution of protoplanets in the N-body simulation agrees with the analytical results. This distribution function provides a useful tool for testing hybrid numerical simulations of late-stage planet formation.

We compare different examples of spherical accretion onto a gravitating mass. Limiting cases include the accretion of a collisionally dominated fluid and the accretion of collisionless particles. We derive expressions for the accretion rate and density profile for semi-collisional accretion which bridges the gap between these limiting cases. Particle crossing of the Hill sphere during the formation of the outer planets is likely to have taken place in the semi-collisional regime.

The ROTSE-IIIc telescope at the H.E.S.S. site, Namibia, obtained the earliest detection of optical emission from a Gamma-Ray Burst (GRB), beginning only 21.8 s from the onset of Swift GRB 050801. The optical lightcurve does not fade or brighten significantly over the first ~250 s, after which there is an achromatic break and the lightcurve declines in typical power-law fashion. The Swift/XRT also obtained early observations starting at 69 s after the burst onset. The X-ray lightcurve shows the same features as the optical lightcurve. These correlated variations in the early optical and X-ray emission imply a common origin in space and time. This behavior is difficult to reconcile with the standard models of early afterglow emission.

We consider a strong ultrarelativistic shock moving through a star whose envelope has a polytrope-like density profile. When the shock is close to the star's outer boundary, its behavior follows the self-similar solution given by Sari (2005) for implosions in planar geometry. Here we outline this solution and find the asymptotic solution as the shock reaches the star's edge. We then show that the motion after the shock breaks out of the star is described by a self-similar solution remarkably like the solution for the motion inside the star. In particular, the characteristic Lorentz factor, pressure, and density vary with time according to the same power laws both before and after the shock breaks out of the star. After emergence from the star, however, the self-similar solution's characteristic position corresponds to a point behind the leading edge of the flow rather than at the shock front, and the relevant range of values for the similarity variable changes. Our numerical integrations agree well with the analytic results both before and after the shock reaches the star's edge.

We derive self similar solutions for ultra-relativistic shock waves propagating into cold material of powerlaw density profile in radius rho ~ r^-k. We treat both implosions and explosions in three geometries: planar, cylindrical and spherical. For spherical explosions these are the first type solutions of Blandford and McKee for k<4 and the second type solutions found by Best and Sari for k>5-sqrt(3/4). In addition we find new, hollow (with evacuated interior), first type solutions that may be applicable for 4<k<17/4. This ``sequence'' with increasing k of first type solutions, hollow first type solutions, and then second type solutions is reminiscent of the non-relativistic sequence. However, while in the non relativistic case there is a range of k which corresponds to a ``gap'' - a range in $k$ with neither first nor second type solutions which separates the hollow first type solutions and the second type solutions, here there is an ``overlap'': a range of k for which current considerations allow for both hollow first and second type solutions. Further understanding is needed to determine which of the two solutions apply in this overlap regime. We provide similar exploration for the other geometries and for imploding configurations. Interestingly, we find a gap for imploding spherical shocks and exploding planar shocks and an overlap for imploding planar solutions. Cylindrical configurations have no hollow solutions and exhibit direct transition from first type to second type solutions, without a gap or an overlap region.

We present optical, near-IR, and radio follow up of sixteen Swift bursts, including our discovery of nine afterglows and a redshift determination for three. These observations, supplemented by data from the literature, provide an afterglow recovery rate of 60% in the optical/near-IR, much higher than in previous missions (BeppoSAX, HETE-2, INTEGRAL, and IPN). The optical/near-IR afterglows of Swift events are on average 1.7 mag fainter at t=12 hr than those of previous missions. The X-ray afterglows are similarly fainter compared to those of pre-Swift bursts. In the radio the limiting factor is the VLA threshold and the detection rate for Swift bursts is similar to that for past missions. The redshift distribution of pre-Swift bursts peaked at z~1, whereas the five Swift bursts with measured redshifts are distributed evenly between 1.3 and 3.2. From these results we conclude that (i) the pre-Swift distributions were biased in favor of bright events and low redshift events, (ii) the higher sensitivity and accurate positions of Swift result in a better representation of the true burst redshift and brightness distributions (which are higher and dimmer, respectively), and (iii) as many as 1/3 of the bursts can be optically dark, as a result of a high redshift and/or dust extinction. We remark that the apparent lack of low redshift, low luminosity Swift bursts, and the lower event rate compared to pre-launch estimates (90 vs. 150 per year), are the result of a threshold that is similar to that of BATSE. In view of these inferences, afterglow observers may find it advisable to make significant changes in follow up strategies of Swift events. [abridged]

The transits of a planet on a Keplerian orbit occur at time intervals exactly equal to the period of the orbit. If a second planet is introduced the orbit is not Keplerian and the transits are no longer exactly periodic. We compute the magnitude of these variations in the timing of the transits, dt. We investigate analytically several limiting cases: (i) interior perturbing planets with much smaller periods; (ii) exterior perturbing planets on eccentric orbits with much larger periods; (iii) both planets on circular orbits with arbitrary period ratio but not in resonance; and (iv) planets on initially circular orbits locked in resonance. Case (iv) is perhaps the most interesting case since some systems are known to be in resonances and the perturbations are the largest. As long as the perturber is more massive than the transiting planet, the timing variations would be of order of the period regardless of the perturber mass! For lighter perturbers, we show that the timing variations are smaller than the period by the perturber to transiting planet mass ratio. An earth mass planet in 2:1 resonance with a 3-day period transiting planet (e.g. HD 209458b) would cause timing variations of order 3 minutes, which would be accumulated over a year. These are easily detectable with current ground-based measurements. For the case of both planets on eccentric orbits, we compute numerically the transit timing variations for several cases of known multiplanet systems assuming they were edge-on. Transit timing measurements may be used to constrain the masses and radii of the planetary system and, when combined with radial velocity measurements, to break the degeneracy between mass and radius of the host star. (abstract truncated)

Planets form in the circumstellar disks of young stars. We review the basic physical processes by which solid bodies accrete each other and alter each others' random velocities, and we provide order-of-magnitude derivations for the rates of these processes. We discuss and exercise the two-groups approximation, a simple yet powerful technique for solving the evolution equations for protoplanet growth. We describe orderly, runaway, neutral, and oligarchic growth. We also delineate the conditions under which each occurs. We refute a popular misconception by showing that the outer planets formed quickly by accreting small bodies. Then we address the final stages of planet formation. Oligarchy ends when the surface density of the oligarchs becomes comparable to that of the small bodies. Dynamical friction is no longer able to balance viscous stirring and the oligarchs' random velocities increase. In the inner-planet system, oligarchs collide and coalesce. In the outer-planet system, some of the oligarchs are ejected. In both the inner- and outer-planet systems, this stage ends once the number of big bodies has been reduced to the point that their mutual interactions no longer produce large-scale chaos. Subsequently, dynamical friction by the residual small bodies circularizes and flattens their orbits. The final stage of planet formation involves the clean up of the residual small bodies. Clean up has been poorly explored.

We address three questions regarding solar system planets. What determined their number? Why are their orbits nearly circular and coplanar? How long did they take to form? Runaway accretion in a disk of small bodies resulted in a tiny fraction of the bodies growing much larger than all the others. This was followed by oligarchic growth during which the big bodies maintained similar masses and uniformly spaced semi-major axes. We depart from standard treatments of planet formation by assuming that as the big bodies got bigger, the small ones got smaller as the result of undergoing a collisional fragmentation cascade. It follows that oligarchy was a brief stage in solar system evolution. When the oligarchs' surface mass density matched that of the small bodies, dynamical friction was no longer able to balance viscous stirring, so their velocity dispersion increased to the extent that their orbits crossed. This marked the end of oligarchy. What happened next differed in the inner and outer parts of the planetary system. In the inner part, where the ratios of the escape velocities from the surfaces of the planets to the escape velocities from their orbits are smaller than unity, big bodies collided and coalesced after their random velocities became comparable to their escape velocities. In the outer part, where these ratios are larger than unity, the random velocities of some of the big bodies continued to rise until they were ejected. In both parts, the number density of the big bodies eventually decreased to the extent that gravitational interactions among them no longer produced large scale chaos. After that their orbital eccentricities and inclinations were damped by dynamical friction from the remaining small bodies. The last and longest stage in planet formation was the clean up of small bodies.

The observed size distribution of Kuiper belt objects (KBOs)--small icy and rocky solar system bodies orbiting beyond Neptune--is well described by a power law at large KBO sizes. However, recent work by Bernstein et al. (2003) indicates that the size spectrum breaks and becomes shallower for KBOs smaller than about 70 km in size. Here we show that we expect such a break at KBO radius ~40 km since destructive collisions are frequent for smaller KBOs. Specifically, we assume that KBOs are rubble piles with low material strength rather than solid monoliths. This gives a power-law slope q~3 where the number N(r) of KBOs larger than a size r is given by $N(r) \propto r^{1-q}$; the break location follows from this slope through a self-consistent calculation. The existence of this break, the break's location, and the power-law slope we expect below the break are consistent with the findings of Bernstein et al. (2003). The agreement with observations indicates that KBOs are effectively strengthless rubble piles.

We report 31 polarimetric observations of the afterglow of GRB 030329 with high signal-to-noise and high sampling frequency. The data imply that the afterglow magnetic field has small coherence length and is mostly random, probably generated by turbulence.

We report 31 polarimetric observations of the afterglow of GRB 030329 with high signal-to-noise and high sampling frequency. We establish the polarization light curve, detect sustained polarization at the percent level, and find significant variability of polarization degree and angle. The data imply that the afterglow magnetic field has small coherence length and is mostly random, probably generated by turbulence.

Many past studies of cosmological gamma-ray bursts (GRBs) have been limited because of the large distance to typical GRBs, resulting in faint afterglows. There has long been a recognition that a nearby GRB would shed light on the origin of these mysterious cosmic explosions, as well as the physics of their fireballs. However, GRBs nearer than z=0.2 are extremely rare, with an estimated rate of localisation of one every decade. Here, we report the discovery of bright optical afterglow emission from GRB 030329. Our prompt dissemination and the brilliance of the afterglow resulted in extensive followup (more than 65 telescopes) from radio through X-ray bands, as well as measurement of the redshift, z=0.169. The gamma-ray and afterglow properties of GRB 030329 are similar to those of cosmological GRBs (after accounting for the small distance), making this the nearest known cosmological GRB. Observations have already securely identified the progenitor as a massive star that exploded as a supernova, and we anticipate futher revelations of the GRB phenomenon from studies of this source.

Past studies suggest that long-duration gamma-ray bursts (GRBs) have a standard energy of E_gamma ~ 10^51 erg in ultra-relativistic ejecta when corrected for asymmetry ("jets"). However, recently a group of sub-energetic bursts, including the peculiar GRB 980425 associated with SN 1998bw (E_gamma ~ 10^48 erg), has been identified. Here we report radio observations of GRB 030329, the nearest burst to date, which allow us to undertake calorimetry of the explosion. Our observations require a two-component explosion: a narrow (5 degrees) ultra-relativistic component responsible for the gamma-rays and early afterglow, and a wide, mildly relativistic component responsible for the radio and optical afterglow beyond 1.5 days. While the gamma-rays are energetically minor, the total energy release, dominated by the wide component, is similar to that of other GRBs. Given the firm link of GRB 030329 with SN 2003dh our result suggests a common origin for cosmic explosions in which, for reasons not understood, the energy in the highest velocity ejecta is highly variable

Planets form in disks around young stars. Interactions with these disks cause them to migrate and thus affect their final orbital periods. We suggest that the connection between planets and disks may be deeper and involve a symbiotic evolution. By contributing to the outward transport of angular momentum, planets promote disk accretion. Here we demonstrate that planets sufficiently massive to open gaps could be the primary agents driving disk accretion. Those having masses below the gap opening threshold drift inward more rapidly than the disk material and can only play a minor role in its accretion. Eccentricity growth during gap formation may involve an even more intimate symbiosis. Given a small initial eccentricity, just a fraction of a percent, the orbital eccentricity of a massive planet may grow rapidly once a mass in excess of the planet's mass has been repelled to form a gap around the planet's orbit. Then, as the planet's radial excursions approach the gap's width, subsequent eccentricity growth slows so that the planet's orbit continues to be confined within the gap.

We employ a fireball model of the GRB explosion to constrain intrinsic and environmental parameters of four events with good broadband afterglow data; GRB970508, GRB980329, GRB980703, and GRB000926. Using standard assumptions of constant circumburst density and no evolution of the fraction of the explosion energy in the post-shock magnetic field, we investigate the uniformity of the derived explosion and shock physics parameters among the events. We find densities that range from those of the ISM to diffuse clouds, energies comparable to the total GRB gamma-ray energy, collimations from near-isotropy to 0.04 rad, substantial electron energy fractions of 10-30% with energy distribution indices of 2.1-2.9 and magnetic energy fractions from 0.2-25%. We also investigate the level to which the data constrain standard model assumptions, such as the magnetic field evolution, and the allowed density profiles of the medium. Fits generally improve slightly with an increasing magnetic energy fraction epB. Good fits can be produced with epB accumulating or decaying with the shock strength over the afterglow as epB~gamma^x; -2<=x<=+1. The data are not very sensitive to increasing density profiles, allowing good fits even with density ~r^10. Some parameter values change by up to an order of magnitude under altered assumptions; parameters of even good fits cannot be taken at face value. The data are sensitive to decreasing densities; r^-2 profiles may produce reasonable fits; steeper profiles, even r^-2.5, will not fit the data.

The two leading interpretations of achromatic breaks that are observed in the light curves of GRBs afterglow are (i) the manifestation of the edge of a jet, which has a roughly uniform energy profile and a sharp edge and (ii) a line of sight effect in jets with a variable energy profile. The first scenario requires the inner engine to produce a jet with a different opening angle each explosion, while the latter requires a standard engine. The physical structure of the jet is a crucial factor in understanding GRB progenitors, and therefore discriminating the two jet scenarios is particularly relevant. In the structured jet case, specific predictions can be made for the distribution of observed break angles $\theta_{\rm break}$, while that distribution is arbitrary in the first scenario. We derive the theoretical distribution for the structured jet model. Specifically, we predict the most common angle to be about 0.12 rad, in rough agreement with the sample. If this agreement would hold as the sample size increases, it would strengthen the case for the standard jet hypothesis. We show that a prediction of this model is that the average viewing angle is an increasing function of the survey sensitivity, and in particular that a mission like {\em Swift} will find the typical viewing angle to be about 0.3 rad. The local event rate predicted by this model is $R_{\rm GRB}(z=0)\sim 0.5 $ Gpc$^{-3}$ yr$^{-1}$.

We present the discovery of GRB 020405 made with the Inter-Planetary Network (IPN). With a duration of 60 s, the burst appears to be a typical long duration event. We observed the 75-square acrminute IPN error region with the Mount Stromlo Observatory's 50-inch robotic telescope and discovered a transient source which subsequently decayed and was also associated with a variable radio source. We identify this source as the afterglow of GRB 020405. Subsequent observations by other groups found varying polarized flux and established a redshift of 0.690 to the host galaxy. Motivated by the low redshift we triggered observations with WFPC2 on-board the Hubble Space Telescope (HST). Modeling the early ground-based data with a jet model, we find a clear red excess over the decaying optical lightcurves that is present between day 10 and day 141 (the last HST epoch). This `bump' has the spectral and temporal features expected of an underlying supernova (SN). In particular, the red color of the putative SN is similar to that of the SN associated with GRB 011121, at late time. Restricting the sample of GRBs to those with z<0.7, a total of five bursts, red bumps at late times are found in GRB 970228, GRB 011121, and GRB 020405. It is possible that the simplest idea, namely that all long duration GRBs have underlying SNe with a modest dispersion in their properties (especially peak luminosity), is sufficient to explain the non detections.

We develop a framework based on energy kicks for the evolution of high-eccentricity long-period orbits with Jacobi constant close to 3 in the restricted circular planar three-body problem where the secondary and primary masses have mass ratio mu << 1. We use this framework to explore mean-motion resonances between the test particle and the secondary mass. This approach leads to (i) a redefinition of resonance orders to reflect the importance of interactions at periapse; (ii) a pendulum-like equation describing the librations of resonance orbits; (iii) an analogy between these new fixed points and the Lagrangian points as well as between librations around the fixed points and the well known tadpole and horseshoe orbits; (iv) a condition a ~ mu^(-2/5) for the onset of chaos at large semimajor axis a; (v) the existence of a range mu < ~5 x 10^(-6) in secondary mass for which a test particle initially close to the secondary cannot escape from the system, at least in the planar problem; (vi) a simple explanation for the presence of asymmetric librations in exterior 1:N resonances.

This is a review article for the book "Compact Stellar X-Ray Sources", Editors W. Lewin and M. van der Klis, to be published by Cambridge University Press in 2003. It covers the phenomenology of cosmic gamma-ray bursts, the theory of their afterglows, and the observations and interpretation of their host galaxies.

Studies of the cosmic gamma-ray bursts (GRBs) and their host galaxies are starting to provide interesting or even unique new insights in observational cosmology. GRBs represent a new way of identifying a population of star-forming galaxies at cosmological redshifts. GRB hosts are broadly similar to the normal field galaxy populations at comparable redshifts and magnitudes, and indicate at most a mild luminosity evolution out to z ~ 1.5 - 2. GRB optical afterglows seen in absorption provide a powerful new probe of the ISM in dense, central regions of their host galaxies, complementary to the traditional studies using QSO absorbers. Some GRB hosts are heavily obscured, and provide a new way to select a population of cosmological sub-mm sources, and a novel constraint on the total obscured fraction of star formation over the history of the universe. Finally, detection of GRB afterglows at z > 6 may provide a unique way to probe the primordial star formation, massive IMF, early IGM, and chemical enrichment at the end of the cosmic reionization era.

We report our discovery and early time optical, near-infrared, and radio wavelength follow-up observations of the afterglow of the gamma-ray burst GRB 021211. Our optical observations, beginning 21 min after the burst trigger, demonstrate that the early afterglow of this burst is roughly three magnitudes fainter than the afterglow of GRB 990123 at similar epochs, and fainter than almost all known afterglows at an epoch of 1d after the GRB. Our near-infrared and optical observations indicate that this is not due to extinction. Combining our observations with data reported by other groups, we identify the signature of a reverse shock. This reverse shock is not detected to a 3-sigma limit of 110 uJy in an 8.46-GHz VLA observation at t=0.10d, implying either that the Lorentz factor of the burst gamma <~ 200, or that synchrotron self-absorption effects dominate the radio emission at this time. Our early optical observations, near the peak of the optical afterglow (forward shock), allow us to characterize the afterglow in detail. Comparing our model to flux upper limits from the VLA at later times, t >~ 1 week, we find that the late-time radio flux is suppressed by a factor of two relative to the >~ 80 uJy peak flux at optical wavelengths. This suppression is not likely to be due to synchrotron self-absorption or an early jet break, and we suggest instead that the burst may have suffered substantial radiative corrections.

We present radio observations of the afterglow of the bright gamma-ray burst GRB980703 made between one day and one year after the burst. These data are combined with published late-time radio measurements and existing optical, near-infrared (NIR) and X-ray observations to create a comprehensive broadband dataset for modeling the physical parameters of the outflow. While a wind-stratified medium cannot be ruled out statistically, it requires a high fraction of the shock energy in the electrons, and so is not favored on theoretical grounds. Instead, the data are consistent with a fireball model in which the ejecta are collimated and expanding into a constant density medium. The radio data cannot be fit with an isotropic shock but instead require a jet break at ~ 3.5 days, not seen at optical wavelengths due to the presence of a a bright host galaxy. The addition of the full radio dataset constrains the self-absorption frequency, giving an estimate of the circumburst density of n ~ 30 cm^-3, a value which differs substantially from previous estimates. This result is consistent with the growing number of GRB afterglows for which broadband modeling yields n ~= 0.1-100 cm^-3, with a typical value ~ 10 cm^-3.

Studies of the cosmic gamma-ray bursts (GRBs) and their host galaxies are now starting to provide interesting or even unique new insights in observational cosmology. Observed GRB host galaxies have a median magnitude R ~ 25 mag, and show a range of luminosities, morphologies, and star formation rates, with a median redshift z ~ 1. They represent a new way of identifying a population of star-forming galaxies at cosmological redshifts, which is mostly independent of the traditional selection methods. They seem to be broadly similar to the normal field galaxy populations at comparable redshifts and magnitudes, and indicate at most a mild luminosity evolution over the redshift range they probe. Studies of GRB optical afterglows seen in absorption provide a powerful new probe of the ISM in dense, central regions of their host galaxies, which is complementary to the traditional studies using QSO absorption line systems. Some GRB hosts are heavily obscured, and provide a new way to select a population of cosmological sub-mm sources. A census of detected optical tranistents may provide an important new way to constrain the total obscured fraction of star formation over the history of the universe. Finally, detection of GRB afterglows at high redshifts (z > 6) may provide a unique way to probe the primordial star formation, massive IMF, early IGM, and chemical enrichment at the end of the cosmic reionization era.

We report on Hubble Space Telescope WFPC2 optical and STIS near ultraviolet MAMA observations, and ground-based optical observations of GRB010222, spanning 15 hrs to 71 days. The observations are well-described by a relativistic blast-wave model with a hard electron-energy distribution, p = 1.57, and a jet transition at t_j=0.93 days. These values are slightly larger than previously found as a result of a correction for the contribution from the host galaxy to the late-time ground-based observations and the larger temporal baseline provided by the Hubble Space Telescope observations. The host galaxy is found to contain a very compact core (size <0.25 arcsec) which coincides with the position of the optical transient. The STIS near ultraviolet MAMA observations allow for an investigation of the extinction properties along the line of sight to GRB010222. We find that the far ultraviolet curvature component (c_4) is rather large. In combination with the low optical extinction A_V =0.11 mag, when compared to the Hydrogen column inferred from X-ray observations, we suggest that this is evidence for dust destruction.

The optical afterglow of gamma-ray burst 020813 was observed for 3 hours with the LRIS spectropolarimeter at the Keck-I telescope, beginning 4.7 hours after the burst was detected by HETE-2. The spectrum reveals numerous metal absorption lines that we identify with two systems at z=1.223 and z=1.255. We also detect an O II 3727 emission line at z=1.255 and we identify this galaxy as the likely host of the GRB. After a correction for Galactic interstellar polarization, the optical afterglow has a linear polarization of 1.8-2.4% during 4.7-7.9 hours after the burst. A measurement of p = 0.80% +/- 0.16% on the following night by Covino et al. demonstrates significant polarization variability over the next 14 hours. The lack of strong variability in the position angle of linear polarization indicates that the magnetic field in the jet is likely to be globally ordered rather than composed of a number of randomly oriented cells. Within the framework of afterglow models with collimated flows, the relatively low observed polarization suggests that the magnetic field components perpendicular and parallel to the shock front are only different by about 20%.

We present radio observations of the afterglow of the bright gamma-ray burst GRB 991208 at frequencies of 1.4, 4.9 and 8.5 GHz, taken between two weeks and 300 days after the burst. The well-sampled radio light curve at 8.5 GHz shows that the peak flux density peaked about 10 days after the burst and decayed thereafter as a power-law t^-1.07. This decay rate is more shallow than the optical afterglow with t^-2.2, which was measured during the first week. These late-time data are combined with extensive optical, millimeter and centimeter measurements and fitted to the standard relativistic blast wave model. In agreement with previous findings, we find that an isotropic explosion in a constant density or wind-blown medium cannot explain these broadband data without modifying the assumption of a single power-law slope for the electron energy distribution. A jet-like expansion provides a reasonable fit to the data. In this case, the flatter radio light curve compared to the optical may be due to emission from an underlying host galaxy, or due to the blastwave making a transition to non-relativistic expansion. The model that best represents the data is a free-form model in which it is assumed that the broadband emission originates from a synchrotron spectrum, while the time-evolution of the break frequencies and peak flux density are solved for explicitly. Although the decay indices for most of the synchrotron parameters are similar to the jet model, the evolution of the cooling break is unusually rapid, and therefore requires some non-standard evolution in the shock. (abridged)

It appears that at least several percent of large Kuiper belt objects are members of wide binaries. Physical collisions are too infrequent to account for their formation. Collisionless gravitational interactions are more promising. These provide two channels for binary formation. In each, the initial step is the formation of a transient binary when two large bodies penetrate each other's Hill spheres. Stabilization of a transient binary requires that it lose energy. Either dynamical friction due to small bodies or the scattering of a third large body can be responsible. Our estimates favor the former, albeit by a small margin. We predict that most objects of size comparable to those currently observed in the Kuiper belt are members of multiple systems. More specifically, we derive the probability that a large body is a member of a binary with semi-major axis of order a. The probability depends upon sigma, the total surface density, Sigma, the surface density of large bodies having radius R, and theta=10^-4, the angle subtended by the solar radius as seen from the Kuiper belt. For (sigma/Sigma)R<a< R/theta, the probability is just (Sigma/rho R)theta^-2, the optical depth of the large bodies divided by the solid angle subtended by the Sun. For R<a<r_u=(sigma/Sigma)R, it varies inversely with semimajor axis and reaches (sigma/rho R)theta^-2 at a=R. Based on current surveys of the Kuiper belt, we estimate Sigma/rho=3 10^-4cm and R=100km. We obtain sigma/rho=0.3cm by extrapolating the surface density deduced for the minimum mass solar nebula. Rough predictions are: outside of the critical separation r_u/a_odot=3'', the binary probability is 0.3%; at separations of 0.2'', comparable to current resolving capabilities, it reaches 5%, in agreement with results from the HST binary survey by Brown.

We present ground-based optical observations of GRB 020124 starting 1.6 hours after the burst, as well as subsequent Very Large Array (VLA) and Hubble Space Telescope (HST) observations. The optical afterglow of GRB 020124 is one of the faintest afterglows detected to date, and it exhibits a relatively rapid decay, $F_\nu\propto t^{-1.60\pm 0.04}$, followed by further steepening. In addition, a weak radio source was found coincident with the optical afterglow. The HST observations reveal that a positionally coincident host galaxy must be the faintest host to date, R>29.5 mag. The afterglow observations can be explained by several models requiring little or no extinction within the host galaxy, A_V~0-0.9 mag. These observations have significant implications for the statistics of the so-called dark bursts (bursts for which no optical afterglow is detected), which are usually attributed to dust extinction within the host galaxy. The faintness and relatively rapid decay of the afterglow of GRB 020124, combined with the low inferred extinction indicate that some dark bursts are intrinsically dim and not dust obscured. Thus, the diversity in the underlying properties of optical afterglows must be observationally determined before substantive inferences can be drawn from the statistics of dark bursts.

GRB 010921 was the first HETE-2 GRB to be localized via its afterglow emission. The low-redshift of the host galaxy, z=0.451, prompted us to undertake intensive multi-color observations with the Hubble Space Telescope with the goal of searching for an underlying supernova component. We do not detect any coincident supernova to a limit 1.34 mag fainter than SN 1998bw at 99.7% confidence, making this one of the most sensitive searches for an underlying SN. Analysis of the afterglow data allow us to infer that the GRB was situated behind a net extinction (Milky Way and the host galaxy) of A_V ~ 1.8 mag in the observer frame. Thus, had it not been for such heavy extinction our data would have allowed us to probe for an underlying SN with brightness approaching those of more typical Type Ib/c supernovae.

We present radio observations of the afterglow of the bright gamma-ray burst GRB980329 made between one month and several years after the burst, a re-analysis of previously published submillimeter data, and late-time optical and near-infrared (NIR) observations of the host galaxy. From the absence of a spectral break in the optical/NIR colors of the host galaxy, we exclude the earlier suggestion that GRB980329 lies at a redshift of z >~5. We combine our data with the numerous multi-wavelength observations of the early afterglow, fit a comprehensive afterglow model to the entire broadband dataset, and derive fundamental physical parameters of the blast-wave and its host environment. Models for which the ejecta expand isotropically require both a high circumburst density and extreme radiative losses from the shock. No low density model (n << 10 cm^{-3}) fits the data. A burst with a total energy of ~ 10^{51} erg, with the ejecta narrowly collimated to an opening angle of a few degrees, driven into a surrounding medium with density ~ 20 cm^{-3}, provides a satisfactory fit to the lightcurves over a range of redshifts.

Using observations from an extensive monitoring campaign with the Hubble Space Telescope we present the detection of an intermediate-time flux excess which is redder in color relative to the afterglow of GRB 011121, currently distinguished as the gamma-ray burst with the lowest known redshift. The red ``bump,'' which exhibits a spectral roll-over at ~7200 Angstrom, is well described by a redshifted Type Ic supernova that occurred approximately at the same time as the gamma-ray burst event. The inferred luminosity is about half that of the bright supernova 1998bw. These results serve as compelling evidence for a massive star origin of long-duration gamma-ray bursts. Models that posit a supernova explosion weeks to months preceding the gamma-ray burst event are excluded by these observations. Finally, we discuss the relationship between spherical core-collapse supernovae and gamma-ray bursts.

Of the cosmological gamma-ray bursts, GRB 011121 has the lowest redshift, z=0.36. More importantly, the multi-color excess in the afterglow detected in the Hubble Space Telescope (HST) light curves is compelling observational evidence for an underlying supernova. Here we present near-infrared and radio observations of the afterglow. We undertake a comprehensive modeling of these observations and those reported in the literature and find good evidence favoring a wind-fed circumburst medium. In detail, we infer the progenitor had a mass loss rate of Mdot ~ 10^-7 / v_w3 Mo/yr where v_w3 is the speed of the wind from the progenitor in units of 10^3 km/s. This mass loss rate is similar to that inferred for the progenitor of SN 1998bw which has been associated with GRB 980425. Our data, taken in conjunction with the HST results of Bloom et al. (2002), provide a consistent picture: the long duration GRB 011121 had a massive star progenitor which exploded as a supernova at about the same time as the GRB event.

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.

Observations of GRB afterglows ranging from radio to X-ray frequencies generate large data sets. Careful analysis of these broad-band data can give us insight into the nature of the GRB progenitor population by yielding such information like the total energy of the burst, the geometry of the fireball and the type of environment into which the GRB explodes. We illustrate, by example, how global, self-consistent fits are a robust approach for characterizing the afterglow emission. This approach allows a relatively simple comparison of different models and a way to determine the strengths and weaknesses of these models, since all are treated self-consistently. Here we quantify the main differences between the broad-band, self-consistent approach and the traditional approach, using GRB000301C and GRB970508 as test cases.

Of all the well localized gamma-ray bursts, GRB 000911 has the longest duration (T_90 ~ 500 s), and ranks in the top 1% of BATSE bursts for fluence. Here, we report the discovery of the afterglow of this unique burst. In order to simultaneously fit our radio and optical observations, we are required to invoke a model involving an hard electron distribution, p ~ 1.5 and a jet-break time less than 1.5 day. A spectrum of the host galaxy taken 111 days after the burst reveals a single emission line, interpreted as [OII] at a redshift z = 1.0585, and a continuum break which we interpret as the Balmer limit at this redshift. Despite the long T_90, the afterglow of GRB 000911 is not unusual in any other way when compared to the set of afterglows studied to date. We conclude that the duration of the GRB plays little part in determining the physics of the afterglow.

We present millimeter and submillimeter wavelength observations and near-infrared K-band imaging toward the bright gamma-ray burst GRB 010222. Over seven epochs the flux density of the source was constant with an average flux density 3.74+/-0.53 mJy at 350 GHz and 1.05+/-0.22 mJy at 250 GHz, giving a spectral index alpha=3.78+/-0.25. We rule out the possibility that this emission originated from the burst or its afterglow and we conclude that it is due to a dusty, high redshift starburst galaxy (SMM J14522+4301). We argue that the host galaxy of GRB 010222 is the most plausible counterpart of SMM J14522+4301, based in part on the centimeter detection of the host at the expected level. The optical/NIR properties of the host galaxy of GRB 010222 suggest that it is a blue, sub-L_*, similar to other GRB host galaxies. This contrasts with the enormous far-infrared luminosity of this galaxy based on our submillimeter detection (L_Bol~4X10^{12} L_odot). We suggest that this GRB host galaxy has a very high star formation rate, SFR 600 M_\odot/yr, most of which is unseen at optical wavelengths.

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.

We develop a model to fit the broadband afterglows of GRBs from the intrinsic parameters of the fireball's synchrotron emission, and apply it to a few well-studied events, with the goal of constraining the intrinsic variability of GRB parameters. We give an example here of fitting to the recent bright event GRB000926.

Observations of GRB host galaxies and their environments in general can provide valuable clues about the nature of progenitors. Bursts are associated with faint, <R> ~ 25 mag, galaxies at cosmological redshifts, <z> ~ 1. The host galaxies span a range of luminosities and morphologies, and appear to be broadly typical for the normal, evolving, actively star-forming galaxy populations at comparable redshifts and magnitudes, but may have somewhat elevated SFR per unit luminosity. There are also spectroscopic hints of massive star formation, from the ratios of [Ne III] and [O II] lines. The observed, unobscured star formation rates are typically a few M_sun/yr, but a considerable fraction of the total star formation in the hosts may be obscured by dust. A census of detected optical afterglows provides a powerful new handle on the obscured fraction of star formation in the universe; the current results suggest that at most a half of the massive star formation was hidden by dust.

Cosmic gamma-ray bursts are one of the great frontiers of astrophysics today. They are a playground of relativists and observers alike. They may teach us about the death of stars and the birth of black holes, the physics in extreme conditions, and help us probe star formation in the distant and obscured universe. In this review we summarise some of the remarkable progress in this field over the past few years. While the nature of the GRB progenitors is still unsettled, it now appears likely that at least some bursts originate in explosions of very massive stars, or at least occur in or near the regions of massive star formation. The physics of the burst afterglows is reasonably well understood, and has been tested and confirmed very well by the observations. Bursts are found to be beamed, but with a broad range of jet opening angles; the mean gamma-ray energies after the beaming corrections are ~ 10^51 erg. Bursts are associated with faint <R> ~ 25 mag) galaxies at cosmological redshifts, with <z> ~ 1. The host galaxies span a range of luminosities and morphologies, but appear to be broadly typical for the normal, actively star-forming galaxy populations at comparable redshifts and magnitudes. Some of the challenges for the future include: the nature of the short bursts and possibly other types of bursts and transients; use of GRBs to probe the obscured star formation in the universe, and possibly as probes of the very early universe; and their detection as sources of high-energy particles and gravitational waves.

GRB 000926 has one of the best-studied afterglows to-date, with multiple X-ray observations, as well as extensive multi-frequency optical and radio coverage. Broadband afterglow observations, spanning from X-ray to radio frequencies, provide a probe of the density structure of the circumburst medium, as well as of the ejecta energetics, geometry, and the physical parameters of the relativistic blastwave resulting from the explosion. We present an analysis of {\em Chandra X-ray Observatory} observations of this event, along with {\em Hubble Space Telescope} and radio monitoring. We combine these data with ground-based optical and IR observations and fit the synthesized afterglow lightcurve using models where collimated ejecta expand into a surrounding medium. We find that we can explain the broadband lightcurve with reasonable physical parameters only if the cooling is dominated by inverse Compton scattering. Excess X-ray emission in the broadband spectrum indicates that we are directly observing a contribution from inverse Compton scattering. It is the first time this has been observed in a GRB afterglow, and it implies that the GRB exploded in a reasonably dense (n~30 cm^{-3}) medium, consistent with a diffuse interstellar cloud environment.

As is well-known, the requirement that gamma ray bursts (GRB's) be optically thin to high energy photons yields a lower limit on the Lorentz factor (\gamma) of the expansion. In this paper, we provide a simple derivation of the lower limit on \gamma due to the annihilation of photon pairs, and correct the errors in some of the previous calculations of this limit. We also derive a second limit on \gamma due to scattering of photons by pair-created electrons and positrons. For some bursts, this limit is the more stringent. In addition, we show that a third limit on \gamma, which is obtained by considering the scattering of photons by electrons which accompany baryons, is nearly always less important than the second limit. Finally, we evaluate these limits for a number of bursts.

The fluences of gamma-ray bursts (GRBs) are measured with a variety of instruments in different detector energy ranges. A detailed comparison of the implied energy releases of the GRB sample requires, then, an accurate accounting of this diversity in fluence measurements which properly corrects for the redshifting of GRB spectra. Here, we develop a methodology to ``k-correct'' the implied prompt energy release of a GRB to a fixed co-moving bandpass. This allows us to homogenize the prompt energy release of 17 cosmological GRBs (using published redshifts, fluences, and spectra) to two common co-moving bandpasses: 20-2000 keV and 0.1 keV-10 MeV (``bolometric''). While the overall distribution of GRB energy releases does not change significantly by using a k-correction, we show that uncorrected energy estimates systematically undercounts the bolometric energy by ~5% to 600%, depending on the particular GRB. We find that the median bolometric isotropic-equivalent prompt energy release is 2.2 x 10^{53} erg with an r.m.s. scatter of 0.80 dex. The typical estimated uncertainty on a given k-corrected energy measurement is ~20%.

Gamma-ray bursts (GRBs) are the most brilliant objects in the Universe but efforts to estimate the total energy released in the explosion -- a crucial physical quantity -- have been stymied by their unknown geometry: spheres or cones. We report on a comprehensive analysis of GRB afterglows and derive their conical opening angles. We find that the gamma-ray energy release, corrected for geometry, is narrowly clustered around 5x10**50 erg. We draw three conclusions. First, the central engines of GRBs release energies that are comparable to ordinary supernovae, suggesting a connection. Second, the wide variation in fluence and luminosity of GRBs is due entirely to a distribution of opening angles. Third, only a small fraction of GRBs are visible to a given observer and the true GRB rate is at least a factor of 500 times larger than the observed rate.

We report on optical, near-infrared and centimeter radio observations of GRB000418 which allow us to follow the evolution of the afterglow from 2 to 200 days after the gamma-ray burst. In modeling these broad-band data, we find that an isotropic explosion in a constant density medium is unable to simultaneously fit both the radio and optical data. However, a jet-like outflow with an opening angle of 10-20 degress provides a good description of the data. The evidence in favor of a jet interpretation is based on the behavior of the radio light curves, since the expected jet break is masked at optical wavelengths by the light of the host galaxy. We also find evidence for extinction, presumably arising from within the host galaxy, with A(V)=0.4 mag, and host flux densities of F_R=1.1 uJy and F_K=1.7 uJy. These values supercede previous work on this burst due to the availability of a broad-band data set allowing a global fitting approach. A model in which the GRB explodes into a wind-stratified circumburst medium cannot be ruled out by these data. However, in examining a sample of other bursts (e.g. GRB990510, GRB000301C) we favor the jet interpretation for GRB000418.

We present multi-color light-curves of the optical afterglow of GRB 000926. Beginning ~1.5 days after the burst, the light-curves of this GRB steepen measurably. The existence of such achromatic breaks are usually taken to be an important observational signature that the ejecta are not expanding isotropically, but rather have a collimated jet-like geometry. If we interpret the data in this context, we derive an opening angle of 5 degrees, which reduces the energy release compared to an isotropic model by a factor of ~275, to 1.7 x 10^51 erg. To fit the data with a simple jet model requires extinction along the line of sight. The derived A_V is in the range 0.11 - 0.82 mag, depending on the adopted extinction law and whether the electrons giving rise to the optical emission are undergoing synchrotron cooling or not. Since this is in excess of the expected extinction from our Galaxy, we attribute this to the GRB host. We note that this extinction is typical of a galactic disk, and therefore the event likely took place in the disk of its host.

Gamma-ray bursts are believed to originate from internal shocks which arise in an irregular relativistic wind. The process has been thought to be inefficient, converting only a few percent of the kinetic energy into gamma-rays. We define ultra efficient internal shocks as those in which the fraction of emitted energy is larger than the fraction of energy given to the radiating electrons at each collision. We show that such a scenario is possible and even plausible. In our model, colliding shells which do not emit all their internal energy are reflected from each other and it causes subsequent collisions, allowing more energy to be emitted. As an example, we obtain about 60% overall efficiency even if the fraction of energy that goes to electrons is $\epsilon_e=0.1$ provided that the shells' Lorentz factor varies between 10 and $10^4$. The numerical temporal profile reflects well the activity of the source which ejects the shells, though numerous collisions take place in this model.

The millimeter wavelength emission from GRB 991208 is the second brightest ever detected, yielding a unique data set. We present here well-sampled spectra and light curves over more than two decades in frequency for a two-week period. This data set has allowed us for the first time to trace the evolution of the characteristic synchrotron self-absorption frequency nu_a and peak frequency nu_m, and the peak flux density F_m: we obtain nu_a \propto t^{-0.15 +- 0.12}, nu_m \propto t^{-1.7 +- 0.4}, and $_m \propto t^{-0.47 +- 0.11}. From the radio data we find that models of homogeneous or wind-generated ambient media with a spherically symmetric outflow can be ruled out. A model in which the relativistic outflow is collimated (a jet) can account for the observed evolution of the synchrotron parameters, the rapid decay at optical wavelengths, and the observed radio to optical spectral flux distributions that we present here, provided that the jet transition has not been fully completed in the first two weeks after the event. These observations provide additional evidence that rapidly decaying optical/X-ray afterglows are due to jets and that such transitions either develop very slowly or perhaps never reach the predicted asymptotic decay F(t) \propto t^{-p}.

We present broad-band radio observations of the afterglow of GRB000301C, spanning from 1.4 to 350 GHz for the period of 3 to 83 days after the burst. This radio data, in addition to measurements at the optical bands, suggest that the afterglow arises from a collimated outflow, i.e. a jet. To test this hypothesis in a self-consistent manner, we employ a global fit and find that a model of a jet, expanding into a constant density medium (ISM+jet), provides the best fit to the data. A model of the burst occurring in a wind-shaped circumburst medium (wind-only model) can be ruled out, and a wind+jet model provides a much poorer fit of the optical/IR data than the ISM+jet model. In addition, we present the first clear indication that the reported fluctuations in the optical/IR are achromatic with similar amplitudes in all bands, and possibly extend into the radio regime. Using the parameters derived from the global fit, in particular a jet break time, t_{jet}=7.5 days, we infer a jet opening angle of \theta=0.2, and consequently the estimate of the emitted energy in the GRB itself is reduced by a factor of 50 relative to the isotropic value, giving E=1.1 \times 10^{51} ergs.

We consider the effects of inverse Compton scattering of synchrotron photons from relativistic electrons in GRB afterglows. We compute the spectrum of the inverse Compton emission and find that it can dominate the total cooling rate of the afterglow for several months or even years after the initial explosion. We demonstrate that the presence of strong inverse Compton cooling can be deduced from the effect it has on the time-evolution of the cooling break in the synchrotron spectral component, and therefore on the optical and X-ray afterglow lightcurves. We then show how the physical interpretation of the observed characteristics of the synchrotron spectrum must be modified to take into consideration this extra source of cooling, and give a revised prescription for computing physical parameters characterizing the expanding shock wave from the observed quantities. We find that for a given set of observables (synchrotron break frequencies and fluxes) there is either no consistent physical interpretation or two of them. Finally we discuss the prospects of directly detecting the inverse Compton emission with Chandra. We argue that such a detection is possible for GRBs exploding in a reasonably dense (n>1 cm^-3) medium.

We discuss the spectrum arising from synchrotron emission by fast cooling (FC) electrons, when fresh electrons are continually accelerated by a strong blast wave, into a power law distribution of energies. The FC spectrum was so far described by four power law segments divided by three break frequencies $\nu_{sa} < \nu_{c} < \nu_{m}$. This is valid for a homogeneous electron distribution. However, hot electrons are located right after the shock, while most electrons are farther down stream and have cooled. This spatial distribution changes the optically thick part of the spectrum, introducing a new break frequency, $\nu_{ac} < \nu_{sa}$, and a new spectral slope, $F_{\nu }\propto \nu^{11/8}$ for $\nu_{ac} < \nu < \nu_{sa}$. The familiar $F_{\nu }\propto \nu ^{2}$ holds only for $\nu < \nu_{ac}$. This ordering of the break frequencies is relevant for typical gamma-ray burst (GRB) afterglows in an ISM environment. Other possibilities arise for internal shocks or afterglows in dense circumstellar winds. We discuss possible implications of this spectrum for GRBs and their afterglows, in the context of the internal-external shock model. Observations of $F_{\nu} \propto \nu^{11/8}$ would enable us to probe scales much smaller than the typical size of the system, and constrain the amount of turbulent mixing behind the shock.

The standard model of Gamma-Ray Bursts afterglows is based on synchrotron radiation from a blast wave produced when the relativistic ejecta encounters the surrounding medium. We reanalyze the refreshed shock scenario, in which slower material catches up with the decelerating ejecta and reenergizes it. This energization can be done either continuously or in discrete episodes. We show that such scenario has two important implications. First there is an additional component coming from the reverse shock that goes into the energizing ejecta. This persists for as long as the re-energization itself, which could extend for up to days or longer. We find that during this time the overall spectral peak is found at the characteristic frequency of the reverse shock. Second, if the injection is continuous, the dynamics will be different from that in constant energy evolution, and will cause a slower decline of the observed fluxes. A simple test of the continuously refreshed scenario is that it predicts a spectral maximum in the far IR or mm range after a few days.

We present wide-band radio observations spanning from 1.4 GHz to 350 GHz of the afterglow of GRB 991216, taken from 1 to 80 days after the burst. The optical and X-ray afterglow of this burst were fairly typical and are explained by a jet fireball. In contrast, the radio light curve is unusual in two respects: (a) the radio light curve does not show the usual rise to maximum flux on timescales of weeks and instead appears to be declining already on day 1 and (b) the power law indices show significant steepening from the radio through the X-ray bands. We show that the standard fireball model, in which the afterglow is from a forward shock, is unable to account for (b) and we conclude that the bulk of the radio emission must arise from a different source. We consider two models, neither of which can be ruled out with the existing data. In the first (conventional) model, the early radio emission is attributed to emission from the reverse shock as in the case of GRB 990123. We predict that the prompt optical emission would have been as bright (or brighter) than 8th magnitude. In the second (exotic) model, the radio emission originates from the forward shock of an isotropically energetic fireball (10^54 erg) expanding into a tenuous medium (10^-4 cm^-3). The resulting fireball would remain relativistic for months and is potentially resolvable with VLBI techniques. Finally, we note that the near-IR bump of the afterglow is similar to that seen in GRB 971214 and no fireball model can explain this bump.

Gamma-ray burst astronomy has undergone a revolution in the last three years, spurred by the discovery of fading long-wavelength counterparts. We now know that at least the long duration GRBs lie at cosmological distances with estimated electromagnetic energy release of 10**51 -- 10**53 erg, making these the brightest explosions in the Universe. In this article we review the current observational state, beginning with the statistics of X-ray, optical, and radio afterglow detections. We then discuss the insights these observations have given to the progenitor population, the energetics of the GRB events, and the physics of the afterglow emission. We focus particular attention on the evidence linking GRBs to the explosion of massive stars. Throughout, we identify remaining puzzles and uncertainties, and emphasize promising observational tools for addressing them. The imminent launch of HETE-2 and the increasingly sophisticated and coordinated ground-based and space-based observations have primed this field for fantastic growth.

The origin of GRBs have been a mystery for almost 30 years. The afterglow observed in the last few years enabled redshift determination for a handful of bursts, and the cosmological origin is now firmly established. Though the distance scale is settled, there still remains orders of magnitude uncertainty in their rate and in the total energy that is released in the explosion due to the possibility that the emission is not spherical but jet-like. Contrary to the GRB itself, the afterglow can be measured up to months and even years after the burst, and it can provide crucial information on the geometry of the ejecta. We review the theory of afterglow from jets and discuss the evidence that at least some of the bursts are not spherical. We discuss the prospects of polarization measurements, and show that this is a powerful tool in constraining the geometry of the explosion.

We investigated the spectral evolution in the 2--700 keV energy band of Gamma-Ray Bursts (GRBs) detected by the Gamma-Ray Burst Monitor (GRBM) and localized with the Wide Field Cameras (WFCs) aboard the BeppoSAX satellite before May 1998. Most of them have been followed-up with the Narrow Field Instruments aboard the same satellite. In the light of these results we discuss open issues on the GRB phenomenon. We find that the optically thin synchrotron shock model (SSM) provides an acceptable representation of most of the time-resolved GRB spectra extending down to 2 keV, except in the initial phases of several bursts and during the whole duration of the quite strong GRB970111, where a low-energy photon depletion with respect to the thin SSM spectrum is observed. We find that the X-ray afterglow starts at about 50% of the GRB duration, and that its fluence, as computed from the WFC light curve, is consistent with the decay law found from the afterglow NFI observations. We also investigate the hydrodynamical evolution of the GRB in our sample and their associated afterglow, when it was detected. We find that the photon index of the latest spectrum of the GRB prompt emission is correlated with the index of the afterglow fading law, when available, as expected on the basis of an external shock of a relativistic fireball.

The strong optical flash observed by ROTSE, as well as the radio flare associated with GRB 990123 can be attributed to the emission of the fireball ejecta, initially heated by the reverse shock. We numerically study the evolution of an adiabatic relativistic fireball interacting with an ambient uniform medium, both in the initial energy transfer stage and in its late evolution. It is shown that the Blandford-McKee solution adequately describes the evolution of the shocked shell quite early on and for as long as the fireball material has relativistic temperatures. In the case where the reverse shock is only mildly relativistic the shocked shell becomes cold almost immediately and the evolution deviates from the Blandford-McKee solution. We derive analytical expressions for the ejecta evolution in its cold regime. This solution gives a good approximation to the numerical results. We estimate the radiation from the fireball ejecta using the numerical hydrodynamic evolution in both cases: cold and hot shells. Surprisingly, we find that both evolutions give rather similar light curves, decaying approximately as $T^{-2}$ in the optical and peaking at about one day in the radio, even though the hydrodynamics is different.

We present multi-frequency radio observations from the afterglow of GRB 980519 beginning 7.2 hours after the gamma-ray burst and ending 63 days later. The fast decline in the optical and X-ray light curves for this burst has been interpreted either as afterglow emission originating from a collimated outflow -- a jet -- or the result of a blast wave propagating into a medium whose density is shaped by the wind of an evolved massive star. These two models predict divergent behavior for the radio afterglow, and therefore, radio observations are capable, in principle, of discriminating between the two. We show that a wind model describes the subsequent evolution of the radio afterglow rather well. However, we see strong modulation of the light curve, which we interpret as diffractive scintillation. These variations prevent us from decisively rejecting the jet model.

We report the discovery of a radio transient VLA 232937.2-235553, coincident with the proposed X-ray afterglow for the gamma-ray burst GRB 981226. This GRB has the highest ratio of X-ray to gamma-ray fluence of all the GRBs detected by BeppoSAX so far and yet no corresponding optical transient was detected. The radio light curve of VLA 232937.2-235553 is qualitatively similar to that of several other radio afterglows. At the sub-arcsecond position provided by the radio detection, optical imaging reveals an extended R=24.9 mag object, which we identify as the host galaxy of GRB 981226. Afterglow models which invoke a jet-like geometry for the outflow or require an ambient medium with a radial density dependence, such as that produced by a wind from a massive star, are both consistent with the radio data. Furthermore, we show that the observed properties of the radio afterglow can explain the absence of an optical transient without the need for large extinction local to the GRB.

We present multi-color optical and two-frequency radio observations of the bright SAX event, GRB 990510. The well-sampled optical decay, together with the radio observations are inconsistent with simple spherical afterglow models. The achromatic optical steepening and the decay of the radio afterglow both occuring at $t \sim 1$ day are evidence for hydrodynamical evolution of the source, and can be most easily interpreted by models where the GRB ejecta are collimated in a jet. Employing a simple jet model to interpret the observations, we derive a jet opening angle of $\theta_o = 0.08$, reducing the isotropic gamma-ray emission of $2.9 \times 10^{53}$ erg by a factor $\sim 300$. If the jet interpretation is correct, we conclude that GRB observations to-date are consistent with an energy for the central source of $E \lsim 10^{52}$ erg.

We report the discovery of a radio counterpart to GRB 990123. In contrast to previous well-studied radio afterglows which rise to peak flux on a timescale of a week and then decay over several weeks to months, the radio emission from this GRB was clearly detected one day after the burst, after which it rapidly faded away. The simplest interpretation of this ``radio flare'' is that it arises from the reverse shock. In the framework of the afterglow models discussed to date, a forward shock origin for the flare is ruled out by our data. However, at late times, some radio afterglow emission (commensurate with the observed late-time optical emission, the optical afterglow) is expected from the forward shock. The relative faintness of the observed late-time radio emission provides an independent indication for a jet-like geometry in this GRB. We use the same radio observations to constrain two key parameters of the forward shock, peak flux and peak frequency, to within a factor of two. These values are inconsistent with the notion advocated by several authors that the prompt optical emission detected by ROTSE smoothly joins the optical afterglow emission. Finally, with hindsight we now recognize another such radio flare and this suggests that one out of eight GRBs has a detectable radio flare. This abundance coupled with the reverse shock interpretation suggests that the radio flare phenomenon has the potential to shed new light into the physics of reverse shocks in GRBs.

We investigate the polarization and proper motion expected in a beamed Gamma-Ray Burst's ejecta. We find that even if the magnetic field has well defined orientation relative to the direction of motion of the shock, the polarization is not likely to exceed 20%. Taking into account the dynamics of beamed ejecta we find that the polarization rises and decays with peak around the jet break time (when the Lorentz factor of the flow is comparable to the initial opening angle of the jet). Interestingly, we find that when the offset of the observer from the center of the beam is large enough, the polarization as function of time has three peaks, and the polarization direction of the middle peak is rotated by 90 degrees relative to the two other peaks. We also show that some proper motion is expected, peaking around the jet break time. Detection of both proper motion and the direction of polarization can determine which component of the magnetic field is the dominant.

In several GRBs afterglows, rapid temporal decay is observed which is inconsistent with spherical (isotropic) blast-wave models. In particular, GRB 980519 had the most rapidly fading of the well-documented GRB afterglows, with t^{-2.05\pm 0.04} in optical as well as in X-rays. We show that such temporal decay is more consistent with the evolution of a jet after it slows down and spreads laterally, for which t^{-p} decay is expected (where p is the index of the electron energy distribution). Such a beaming model would relax the energy requirements on some of the more extreme GRBs by a factor of several hundreds. It is likely that a large fraction of the weak (or no) afterglow observations are also due to the common occurrence of beaming in GRBs, and that their jets have already transitioned to the spreading phase before the first afterglow observations were made. With this interpretation, a universal value of p~2.5 is consistent with all data.

Afterglow, or long-lived emission, has now been detected from about a dozen well-positioned gamma-ray bursts. Distance determinations made by measuring optical emission lines from the host galaxy, or absorption lines in the afterglow spectrum, place the burst sources at significant cosmological distances, with redshifts ranging from ~1--3. The energy required to produce the bright gamma-ray flashes is enormous: up to ~10^{53} erg or 10 percent of the rest mass energy of a neutron star, if the emission is isotropic. Here we present the discovery of the optical afterglow and the redshift of GRB 990123, the brightest well-localized GRB to date. With our measured redshift of >1.6, the inferred isotropic energy release exceeds the rest mass of a neutron star thereby challenging current theoretical models for the origin of GRBs. We argue that the optical and IR afterglow measurements reported here may provide the first observational evidence of beaming in a GRB, thereby reducing the required energetics to a level where stellar death models are still tenable.

We present deep images of the field of gamma-ray burst (GRB) 990123 obtained in a broad-band UV/visible bandpass with the Hubble Space Telescope, and deep near-infrared images obtained with the Keck-I 10-m telescope. Both the HST and Keck images show that the optical transient (OT) is clearly offset by 0.6 arcsec from an extended object, presumably the host galaxy. This galaxy is the most likely source of the metallic-line absorption at z = 1.6004 seen in the spectrum of the OT. With magnitudes V_{C} ~ 24.6 +/- 0.2 and K = 21.65 +/- 0.30 mag this corresponds to an L ~ 0.7 L_* galaxy, assuming that it is located at z = 1.6. The estimated unobscured star formation rate is SFR ~ 6 M_sun/yr, which is not unusually high for normal galaxies at comparable redshifts. The strength of the observed metallic absorption lines is suggestive of a relatively high metallicity of the gas, and thus of a chemically evolved system which may be associated with a massive galaxy. It is also indicative of a high column density of the gas, typical of damped Ly-alpha systems at high redshifts. We conclude that this is the host galaxy of GRB 990123. No other obvious galaxies are detected within the same projected radius from the OT. There is thus no evidence for strong gravitational lensing magnification of this burst, and some alternative explanation for its remarkable energetics may be required. The observed offset of the OT from the center of its apparent host galaxy, 5.5 +/- 0.9 proper kpc (projected) in the galaxy's rest-frame, both refutes the possibility that GRBs are related to galactic nuclear activity and supports models of GRBs which involve the death and/or merger of massive stars. Further, the HST image suggests an intimate connection of GRB 990123 and a star-forming region.

We compare the ongoing observations of the remarkable burst GRB990123, the mother of all bursts, with the predictions of the afterglow theory. We show that the observations agree with the recent prediction that a reverse shock propagating into the ejecta would produce a very strong prompt optical flash. This reverse shock has also produced the 8.46GHz radio signal, observed after one day. The forward shock, which propagates into the ISM is the origin of the classical afterglow. It has produced the prompt X-ray signal as well as the late optical and IR emission. It would most likely produce a radio emission within the next few weeks. The observations suggest that the initial Lorentz factor of the ejecta was $\sim 200$. Within factors of order unity, this crude model explains all current observations of GRB990123.

According to the internal-external shocks model for $\gamma $-ray bursts (GRBs), the GRB is produced by internal shocks within a relativistic flow while the afterglow is produced by external shocks with the ISM. We explore the early afterglow emission. For short GRBs the peak of the afterglow will be delayed, typically, by few dozens of seconds after the burst. For long GRBs the early afterglow emission will overlap the GRB signal. We calculate the expected spectrum and the light curves of the early afterglow in the optical, X-ray and $\gamma $-ray bands. These characteristics provide a way to discriminate between late internal shocks emission (part of the GRB) and the early afterglow signal. If such a delayed emission, with the characteristics of the early afterglow, will be detected it can be used both to prove the internal shock scenario as producing the GRB, as well as to measure the initial Lorentz factor of the relativistic flow. The reverse shock, at its peak, contains energy which is comparable to that of the GRB itself, but has a much lower temperature than that of the forward shock so it radiates at considerably lower frequencies. The reverse shock dominates the early optical emission, and an optical flash brighter than 15th magnitude, is expected together with the forward shock peak at x-rays or $\gamma$-rays. If this optical flash is not observed, strong limitations can be put on the baryonic contents of the relativistic shell deriving the GRBs, leading to a magnetically dominated energy density.

We calculate the expected spectrum and light curves of the early afterglow. For short GRBs the peak of the afterglow will be delayed, typically, by a few dozens of seconds after the burst. The x-ray and gamma-ray characteristics of this delayed emission provide a way to discriminate between late internal shocks emission (part of the GRB) and the early afterglow signal. Detection of this delayed emission will prove the internal shock scenario as producing the GRB, and will pinpoint the initial Lorentz factor gamma_0. In the optical band, the dominant emission arises from the reverse shock. This shock, carries an energy comparable to that of the forward shock. It radiates this energy at much lower frequencies, producing a short optical flash of 15th magnitude or brighter.

We calculate the light curve and spectra near the peak and the self absorption break, for an adiabatic blast wave described by the Blandford-McKee solution, considering the emission from the whole region behind the shock front. The expected light curve and spectra are flat near the peak. This rules out the interpretation of the sharp peak observed in the optical afterglow of GRB970508 as the expected peak of the light curve. The observed image of an afterglow is calculated for a broad range of frequencies. We show that for frequencies below the self absorption frequency the image is rather homogeneous, as opposed to the bright ring at the outer edge and dim center, which appear at higher frequencies. We fit the observed spectra of GRB970508 to the detailed theory and obtain estimates of the physical parameters of this burst.

GRB afterglow is reasonably described by synchrotron emission from relativistic blast waves at cosmological distances. We perform detailed calculations taking into account the effect of synchrotron self absorption. We consider emission from the whole region behind the shock front, and use the Blandford McKee self similar solution to describe the fluid behind the shock. We calculate the spectra and the observed image of a GRB afterglow near the self absorption frequency $\nu_a$ and derive an accurate expression for $\nu_a$. We show that the image is rather homogeneous for $\nu<\nu_a$, as opposed to the bright ring at the outer edge and dim center, which appear at higher frequencies. We compare the spectra we obtain to radio observations of GRB970508. We combine the calculations of the spectra near the self absorption frequency with other parts of the spectra and obtain revised estimates for the physical parameters of the burst: $E_{52}=0.53$, $\epsilon_e=0.57$, $\epsilon_B=0.0082$, $n_1=5.3$. These estimates are different by up to two orders of magnitude than the estimates based on an approximate spectrum.

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.

We derive a self similar solution for the propagation of an extreme relativistic (or Newtonian) radiative spherical blast wave into a surrounding cold medium. The solution is obtained under the assumption that the radiation process is fast, it takes place only in the vicinity of the shock and that it radiates away a fixed fraction of the energy generated by the shock. In the Newtonian regime these solutions generalize the Sedov-Taylor adiabatic solution and the pressure-driven fully radiative solution. In the extreme relativistic case these solutions generalize the Blandford-McKee adiabatic solution. They provide a new fully radiative extreme relativistic solution which is different from the Blandford-McKee fully radiative relativistic solution. This new solution develops a hot interior which causes it to cool faster than previous estimates. We find that the energy of the blast wave behaves as a power law of the location of the shock. The power law index depends on the fraction of the energy emitted by the shock. We obtain an analytic solution for the interior of the blast wave. These new solutions might be applicable to the study of GRB afterglow or SNRs.

We study numerically the evolution of an adiabatic relativistic fireball expanding into a cold uniform medium. We follow the stages of initial free expansion and acceleration, coasting and then deceleration and slowing down to a non-relativistic velocity. We compare the numerical results with simplified analytical estimates. We show that the relativistic self similar Blandford-McKee solution describes well the relativistic deceleration epoch. It is an excellent approximation throughout the relativistic deceleration stage, down to $\gamma \sim 5$, and a reasonable approximation even down to $\gamma \sim 2$ though the solution is rigorous only for $\gamma \gg 1$. We examine the transition into the Blandford-McKee solution, and the transition from the solution to the non-relativistic self similar Sedov-Taylor solution. These simulations demonstrate the attractive nature of the Blandford-McKee solution and its stability to radial perturbations.

The recently discovered GRB afterglow is believed to be described reasonably well by synchrotron emission from a slowing down relativistic shell that collides with an external medium. To compare theoretical models with afterglow observations we calculate here the broad band spectrum and corresponding light curve of synchrotron radiation from a power-law distribution of electrons in an expanding relativistic shock. Both the spectrum and the light curve consist of several power-law segments. The light curve is constructed under two limiting models for the hydrodynamical evolution of the shock: fully adiabatic and fully radiative. We compare the results with observations of $\gamma$-ray burst afterglows.

The detection of delayed emission in X-ray, optical and radio wave length, ``afterglow'', following a gamma-ray burst can be described by the emission of a relativistic shell decelerating upon collision with the ISM. We show that the observed radiation surface have well defined bright edges. We derive an explicit expression for the size as a function of time, and obtain the surface brightness distribution. This might be directly observed if the burst occurs at small redshift so that its radio signal can be resolved. The size and shape are relevant for detailed analysis of scintillation or microlensing. We show that the effective Lorentz factor depends on the observed frequency and it is higher for frequencies above the synchrotron typical frequency (optical and X-ray) than for low frequencies (radio). Consequently transition to non relativistic evolution, will be observed first in low frequencies and only a factor of ~2 later in the high frequencies.

The detection of delayed emission at X-ray optical and radio wavelengths, ``after-glow'', suggests that the relativistic shell which emitted the initial GRB due to internal shocks decelerates on encountering an external medium, giving rise to the after-glow. At a very early stage (few seconds), the observed bolometric luminosity increases as t^2. On longer time scales (more than about 10s), the luminosity drops as t^{-1}. If the main burst is long enough, an intermediate stage of constant luminosity will form. In this case, the after-glow overlaps the main burst, otherwise there is a time separation between the two. On the long time scale, the flow decelerate in a self similar way, reaching non relativistic velocities after about 30days. Due to the deceleration and the accumulation of ISM material, the relation between the observed time, the shock radius, and its Lorentz factor, is given by t=R/16\gamma^2c which is a factor of 8 different from the usual expression. The majority of particles are those of the original ejecta (and not the ISM) up to about 900s. These particles reach sub-relativistic velocities on a time scale of about 2hours, well before the flow becomes sub-relativistic. Therefore the ejecta particles are probably unimportant for most of the after-glow radiation. We show that even though only a small fraction of the energy is given to the electrons, most of the energy can be radiated over time. If this fraction is greater than 10% radiation losses will significantly influence the hydrodynamical evolution at early times (less than 1day).

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.

The gamma-ray burst source GRB970228 has been observed after a delay of 8--12 hours in X-rays and after one day in visible and near infrared light. This marks the first detection of emission at lower frequencies following the gamma-ray observation of a GRB and the first detection of any visible counterpart to a GRB. We consider possible delayed visible and X-ray emission mechanisms, and conclude that the intrinsic gamma-ray activity continued at a much reduced intensity for at least a day. There are hints of such continued activity in other GRB, and future observations can decide if this is true of GRB in general. The observed multi-band spectrum of GRB970228 agrees with the predictions of relativistic shock theory when the flux is integrated over a time longer than that required for a radiating electron to lose its energy.

Relativistic shock models of gamma-ray bursts may be tested by comparison of their predicted low energy asymptotic spectral indices s to observations. Synchrotron radiation theory predicts that the instantaneous spectrum has s = - 1/3 and the spectrum integrated over the radiative decay of the electrons' energies has s = 1/2, with other cases lying between these limits. We examine the spectra of 11 bursts obtained by the Large Area Detectors on BATSE. One agrees with the predicted instantaneous spectrum, as does the initial portion of a second, and three are close to the predicted integrated spectrum. All the observed asymptotic spectral slopes lie in the predicted range. This evidence for relativistic shocks is independent of detailed models of bursts and of assumptions about their distances. Radiation observed with the predicted instantaneous spectrum has a comparatively smooth time dependence, consistent with the necessarily long radiation time, while that with the predicted integrated spectrum has a spiky time dependence, consistent with the necessarily short radiation time.

The highly variable temporal structure observed in most GRBs provides us with unexpected clues. We show that variable GRBs cannot be produced by external shocks models and consequently cannot be produced by an ``explosive'' inner engine. The observed temporal structure must reflect the activity of the inner engine that must be producing unsteady and irregular ``wind'' which is converted to radiation via internal shocks.

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.

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.

A leading mechanism for producing cosmological gamma-ray bursts (GRBs) is via ultra-relativistic particles in an expanding fireball. The kinetic energy of the particles is converted into thermal energy in a forward shock and a reverse shock, when the outward flowing particles encounter the interstellar medium. The thermal energy is radiated via synchrotron emission and Comptonization. We estimate the synchrotron cooling time scale of the shocked regions for electrons of various Lorentz factors, focusing in particular on those electrons whose radiation falls within the energy detection range of the BATSE detectors. We find that in order to produce the rapid variability observed in most bursts the magnetic energy density in the shocked material must be greater than about 1\% of the thermal energy density. Additionally, the electrons must be nearly in equipartition with the protons, otherwise the radiative efficiency is low. Inverse Compton scattering can increase the cooling rate of the relevant electrons but the Comptonized emission itself is never within the BATSE range. These arguments allow us to pinpoint the conditions within the radiating regions in GRBs and to determine the important radiation processes. In addition, the model predicts that the duty cycle of intensity variations in GRB light curves should be nearly independent of burst duration, and should scale inversely as the square root of the observed photon energy. Both correlations are in agreement with observations. The model also provides a plausible explanation for the bimodal distribution of burst durations. There is no explanation, however, for the characteristic break energy in GRB spectra.

We calculate the hydrodynamic time scales for a spherical ultra-relativistic shell that is decelerated by the ISM and discuss the possible relations between these time scales and the observed temporal structure in $\gamma$-ray bursts. We suggest that the bursts' duration is related to the deceleration time, the variability is related to the ISM inhomogeneities and precursors are related to internal shocks within the shell. Good agreement can be achieved for these quantities with reasonable, not fined tuned, astrophysical parameters. The difference between Newtonian and relativistic reverse shocks may lead to the observed bimodal distribution of bursts' durations.

We calculate the hydrodynamic time scales for a spherical ultra-relativistic shell that is decelerated by the ISM and discuss the possible relations between these time scales and the observed temporal structure in $\gamma$-ray bursts. We suggest that the bursts' duration is related to the deceleration time, the variability is related to the ISM inhomogeneities and precursors are related to internal shocks within the shell. Good agreement can be achieved for these quantities with reasonable, not fined tuned, astrophysical parameters. The difference between Newtonian and relativistic reverse shocks may lead to the observed bimodal distribution of bursts' durations.