Halting Migration: Numerical Calculations of Corotation Torques in the Weakly Nonlinear Regime

Planets in their formative years can migrate due to the influence of gravitational torques in the protoplanetary disk they inhabit. For low-mass planets in an isothermal disk, it is known that there is a strong negative torque on the planet due to its linear perturbation to the disk, causing fast inward migration. The current investigation demonstrates that in these same isothermal disks, for intermediate-mass planets, there is a strong positive nonlinear corotation torque due to the effects of gas being pulled through a gap on horseshoe orbits. For intermediate-mass planets, this positive torque can partially or completely cancel the linear (Type I) torque, leading to slower or outward migration, even in an isothermal disk. The effect is most significant for Super-Earth and Sub-Jovian planets, during the transition from a low-mass linear perturber to a non-linear gap-opening planet, when the planet has opened a so-called 'partial gap'. In this study, numerical calculations of planet-disk interactions calculate these torques explicitly, and scalings are empirically constructed for migration rates in this weakly nonlinear regime. These results find outward migration is possible for planets with masses in the 20 - 100 Earth Mass range.Comment: ApJ Accepte

A One-Dimensional Model for Rayleigh-Taylor Instability in Supernova Remnants

This study presents a method for approximating the multidimensional effects of Rayleigh-Taylor instability as a modification of the one-dimensional hydro equations. This modification is similar to the Shakura-Sunyaev {\alpha} prescription for modeling the coarse-grained effects of turbulence in astrophysical disks. The model introduces several dimensionless tunable parameters that are calibrated by comparing with high-resolution two-dimensional axisymmetric numerical calculations of Rayleigh-Taylor unstable flows. A complete description of the model is presented, along with a handful of test problems that demonstrate the extent to which the one-dimensional model is able to reproduce multidimensional effects.Comment: ApJ Accepte

A Simple Analytical Model for Gaps in Protoplanetary Disks

An analytical model is presented for calculating the surface density as a function of radius $\Sigma(r)$ in protoplanetary disks in which a planet has opened a gap. This model is also applicable to circumbinary disks with extreme binary mass ratios. The gap profile can be solved for algebraically, without performing any numerical integrals. In contrast with previous one-dimensional gap models, this model correctly predicts that low-mass (sub-Jupiter) planets can open gaps in sufficiently low-viscosity disks, and it correctly recovers the power-law dependence of gap depth on planet-to-star mass ratio $q$, disk aspect ratio $h/r$, and dimensionless viscosity $\alpha$ found in previous numerical studies. Analytical gap profiles are compared with numerical calculations over a range of parameter space in $q$, $h/r$, and $\alpha$, demonstrating accurate reproduction of the "partial gap" regime, and general agreement over a wide range of parameter space

Shock Corrugation by Rayleigh-Taylor Instability in GRB Afterglow Jets

Afterglow jets are Rayleigh-Taylor unstable and therefore turbulent during the early part of their deceleration. There are also several processes which actively cool the jet. In this letter, we demonstrate that if cooling significantly increases the compressibility of the flow, the turbulence collides with the forward shock, destabilizing and corrugating it. In this case, the forward shock is turbulent enough to produce the magnetic fields responsible for synchrotron emission via small scale turbulent dynamo. We calculate light curves assuming the magnetic field is in energy equipartition with the turbulent kinetic energy and discover that dynamic magnetic fields are well-approximated by a constant magnetic-to- thermal energy ratio of 1%, though there is a sizeable delay in the time of peak flux as the magnetic field turns on only after the turbulence has activated. The reverse shock is found to be significantly more magnetized than the forward shock, with a magnetic-to-thermal energy ratio of order 10%. This work motivates future Rayleigh-Taylor calculations using more physical cooling models.Comment: ApJ Accepte

Synchrotron Magnetic Fields from Rayleigh-Taylor Instability in Supernovae

Synchrotron emission from a supernova necessitates a magnetic field, but it is unknown how strong the relevant magnetic fields are, and what mechanism generates them. In this study, we perform high-resolution numerical gas dynamics calculations to determine the growth of turbulence due to Rayleigh-Taylor instability, and the resulting kinetic energy in turbulent fluctuations, to infer the strength of magnetic fields amplified by this turbulence. We find that Rayleigh-Taylor instability can produce turbulent fluctuations strong enough to amplify magnetic fields to a few percent of equipartition with the thermal energy. This turbulence stays concentrated near the reverse shock, but averaging this magnetic energy throughout the shocked region (weighting by emissivity) sets the magnetic fields at a minimum of 0.3 percent of equipartition. This suggests a minimum effective magnetic field strength ($\epsilon_B > 0.003$) which should be present in all interacting supernovae

On the Deceleration and Spreading of Relativistic Jets I: Jet Dynamics

Jet breaks in gamma ray burst (GRB) afterglows provide a direct probe of their collimation angle. Modeling a jet break requires an understanding of the "jet spreading" process, whereby the jet transitions from a collimated outflow into the spherical Sedov-Taylor solution at late times. Currently, direct numerical calculations are the most accurate way to capture the deceleration and spreading process, as analytical models have previously given inaccurate descriptions of the dynamics. Here (in paper I) we present a new, semi-analytical model built empirically by performing relativistic numerical jet calculations and measuring the relationship between Lorentz factor and opening angle. We then calculate the Lorentz factor and jet opening angle as a function of shock radius and compare to the numerical solutions. Our analytic model provides an efficient means of computing synthetic GRB afterglow light curves and spectra, which is the focus of paper II.Comment: ApJ Submitte

Eccentric Jupiters via Disk-Planet Interactions

Numerical hydrodynamics calculations are performed to determine conditions under which giant planet eccentricities can be excited by parent gas disks. Unlike in other studies, Jupiter-mass planets are found to have their eccentricities amplified --- provided their orbits start eccentric. We disentangle the web of co-rotation, co-orbital, and external resonances to show that this finite-amplitude instability is consistent with that predicted analytically. Ellipticities can grow until they reach of order the disk's aspect ratio, beyond which the external Lindblad resonances that excite eccentricity are weakened by the planet's increasingly supersonic epicyclic motion. Forcing the planet to still larger eccentricities causes catastrophic eccentricity damping as the planet collides into gap walls. For standard parameters, the range of eccentricities for instability is modest; the threshold eccentricity for growth ($\sim$$0.04$) is not much smaller than the final eccentricity to which orbits grow ($\sim$$0.07$). If this threshold eccentricity can be lowered (perhaps by non-barotropic effects), and if the eccentricity driving documented here survives in 3D, it may robustly explain the low-to-moderate eccentricities $\lesssim 0.1$ exhibited by many giant planets (including Jupiter and Saturn), especially those without planetary or stellar companions.Comment: Accepted to ApJ with added references and minor revision

A "Boosted Fireball" Model for Structured Relativistic Jets

We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor \gamma_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum frame the fireball expands isotropically, converting its internal energy into radially expanding flow with asymptotic Lorentz factor \eta_0 ~ E/M. In the lab frame the flow is beamed, expanding with Lorentz factor \Gamma = 2 \eta_0 \gamma_B in the direction of its initial bulk motion and with characteristic opening angle \theta_0 ~ 1/\gamma_B. The flow is jet-like with \Gamma \theta_0 ~ 2 \eta_0 such that jets with \Gamma > 1/\theta_0 are naturally produced. The choice \eta_0 ~ \gamma_B ~ 10 yields a jet with \Gamma ~ 200 on-axis and angular structure characterized by opening angle \theta_0 ~ 0.1 of relevance for cosmological GRBs, while \gamma_B >~ 1 may be relevant for low-luminosity GRBs. The model produces a family of outflows, of relevance for different relativistic phenomena with structures completely determined by \eta_0 and \gamma_B. We calculate the energy per unit solid angle for the model and use it to compute light curves for comparison with the widely used top-hat model. The jet break in the boosted fireball light curve is greatly subdued when compared to the top-hat model because the edge of the jet is smoother than for a top-hat. This may explain missing jet breaks in afterglow light curves.Comment: ApJ Accepte

From Engine to Afterglow: Collapsars Naturally Produce Top-Heavy Jets and Early-Time Plateaus in Gamma Ray Burst Afterglows

We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through the late afterglow phase. The resultant outflow includes a jet core which is highly relativistic after breaking out of the star, but becomes baryon-loaded after colliding with a massive outer shell, corresponding to mass from the stellar atmosphere of the progenitor star which became trapped in front of the jet core at breakout. The prompt emission produced before or during this collision would then have the signature of a high Lorentz factor jet, but the afterglow is produced by the amalgamated post-collision ejecta which has more inertia than the original highly relativistic jet core and thus has a delayed deceleration. This naturally explains the early light curve behavior discovered by Swift, including a steep decay and a long plateau, without invoking late-time energy injection from the central engine. The numerical simulation is performed continuously from engine to afterglow, covering a dynamic range of over ten orders of magnitude in radius. Light curves calculated from the numerical output demonstrate that this mechanism reproduces basic features seen in early afterglow data. Initial steep decays are produced by internal shocks, and the plateau corresponds to the coasting phase of the outflow.Comment: ApJ Accepte

Gap Opening by Extremely Low Mass Planets in a Viscous Disk

By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ~10^{-4}) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the low-viscosity limit, we find no minimum mass necessary to open a gap for a planet held on a fixed orbit. In particular, a Neptune-mass planet will open a gap in a minimum mass solar nebula with suitably low viscosity (\alpha <10^{-4}). We find that the mass threshold scales as the square root of viscosity in the low mass regime. This is because the gap width for critical planet masses in this regime is a fixed multiple of the scale height, not of the Hill radius of the planet.Comment: ApJ accepte