In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), where the inner boundary is modeled by removing mass within a finite zone, the inner density profile typically becomes anomalously depleted. Here we show how the boundary condition should be applied in such codes, via a simple modification of the usual approach: when one removes mass, one must speed up the remaining material so that the disk's angular momentum is unchanged. We show with both 1D and 2D moving-mesh (AREPO) simulations that this scheme works as desired in viscously evolving disks. It produces no spurious density depletions and is independent of the mass removal rate, provided that the disk is adequately resolved and that the mass removal rate is not so extreme as to trigger instabilities. This "torque-free" mass removal technique permits the use of quasi-Lagrangian codes to simulate viscously evolving disks, while including a variety of additional effects. As an example, we apply our scheme to a 2D simulation of an accretion disk perturbed by a very massive planet, in which the disk is evolved to viscous steady state.

We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.

We investigate how planets interact with viscous accretion disks, in the limit that the disk is sufficiently low mass that the planet migrates more slowly than the disk material. In that case, the disk's surface density profile is determined by the disk being in viscous steady state (VSS), while overflowing the planet's orbit. We compute the VSS profiles with 2D hydrodynamical simulations, and show that disk material piles up behind the planet, with the planet effectively acting as a leaky dam. Previous 2D hydrodynamical simulations missed the pileup effect because of incorrect boundary conditions, while previous 1D models greatly overpredicted the pileup due to the neglect of non-local deposition. Our simulations quantify the magnitude of the pileup for a variety of planet masses and disk viscosities. We also calculate theoretically the magnitude of the pileup for moderately deep gaps, showing good agreement with simulations. For very deep gaps, current theory is inadequate, and we show why and what must be understood better. The pileup is important for two reasons. First, it is observable in directly imaged protoplanetary disks, and hence can be used to diagnose the mass of a planet that causes it or the viscosity within the disk. And second, it determines the planet's migration rate. Our simulations determine a new Type-II migration rate (valid for low mass disks), and show how it connects continuously with the well-verified Type-I rate.

Accretion disks can be eccentric: they support $m=1$ modes that are global and slowly precessing. But whether the modes remain trapped in the disk---and hence are long-lived---depends on conditions at the outer edge of the disk. Here we show that in disks with realistic boundaries, in which the surface density drops rapidly beyond a given radius, eccentric modes are trapped and hence long-lived. We focus on pressure-only disks around a central mass, and show how this result can be understood with the help of a simple second-order WKB theory. We show that the longest lived mode is the zero-node mode in which all of the disk's elliptical streamlines are aligned, and that this mode decays coherently on the viscous timescale of the disk. Hence such a mode, once excited, will live for the lifetime of the disk. It may be responsible for asymmetries seen in recent images of protoplanetary disks.

Saturn is ringing weakly. Exquisite data from the Cassini mission reveal the presence of f-mode oscillations as they excite density waves in Saturn's rings. These oscillations have displacement amplitudes of order a metre on Saturn's surface. We propose that they result from large impacts in the past. Experiencing little dissipation inside Saturn on account of its weak luminosity, f-modes may live virtually forever; but the very ring waves that reveal their existence also remove energy from them, in 10^4 to 10^7 yrs for the observed f-modes (spherical degree 2-10). We find that the largest impacts that arrive during these times excite the modes to their current levels, with the exception of the few lowest degree modes. To explain the latter, either a fortuitously large impact in the recent past, or a new source of stochastic excitation, is needed. We extend this scenario to Jupiter which has no substantial rings. With an exceedingly long memory of past bombardments, Jovian f-modes and p-modes can acquire much higher amplitudes, possibly explaining past reports of radial-velocity detections, and detectable by the Juno spacecraft.

Accretion disks around stars, or other central massive bodies, can support long-lived, slowly precessing $m=1$ disturbances in which the fluid motion is nearly Keplerian with non-zero eccentricity. We study such `slow modes' in disks that are subject to both pressure and self-gravity forces. We derive a second-order WKB dispersion relation that describes the dynamics quite accurately, and show that the apparently complicated nature of the various modes can be understood in a simple way with the help of a graphical method. We also solve the linearized fluid equations numerically, and show that the results agree with the theory. We find that when self-gravity is weak ($Q\gtrsim 1/h$, where $Q$ is Toomre's parameter, and $h$ is the disk aspect ratio) the modes are pressure dominated. But when self-gravity is strong ($1<Q\lesssim 1/h$), two kinds of gravity-dominated modes appear: one is an aligned elliptical pattern and the other is a one-armed spiral. In the context of protoplanetary disks, we suggest that if the radial eccentricity profile can be measured, it could be used to determine the total disk mass.

We derive a criterion for the onset of chaos in systems consisting of two massive, eccentric, coplanar planets. Given the planets' masses and separation, the criterion predicts the critical eccentricity above which chaos is triggered. Chaos occurs where mean motion resonances overlap, as in Wisdom (1980)'s pioneering work. But whereas Wisdom considered only nearly circular planets, and hence examined only first order resonances, we extend his results to arbitrarily eccentric planets (up to crossing orbits) by examining resonances of all orders. We thereby arrive at a simple expression for the critical eccentricity. We do this first for a test particle in the presence of a planet, and then generalize to the case of two massive planets, based on a new approximation to the Hamiltonian (Hadden, in prep). We then confirm our results with detailed numerical simulations. Finally, we explore the extent to which chaotic two-planet systems eventually result in planetary collisions.

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

We conduct a uniform analysis of the transit timing variations (TTVs) of 145 planets from 55 Kepler multiplanet systems to infer planet masses and eccentricities. Eighty of these planets do not have previously reported mass and eccentricity measurements. We employ two complementary methods to fit TTVs: Markov chain Monte Carlo simulations based on N-body integration and an analytic fitting approach. Mass measurements of 49 planets, including 12 without previously reported masses, meet our criterion for classification as robust. Using mass and radius measurements, we infer the masses of planets' gaseous envelopes for both our TTV sample as well as transiting planets with radial velocity observations. Insight from analytic TTV formulae allows us to partially circumvent degeneracies inherent to inferring eccentricities from TTV observations. We find that planet eccentricities are generally small, typically a few percent, but in many instances are non-zero.

Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was placed at >~1AU distance from the stellar host. This finding is at odds with formation theories that invoke the migration of such planets from semi-major axes >~1AU due to secular dynamical processes (e.g., secular chaos, Lidov-Kozai oscillations) coupled with tidal dissipation. In fact, the results of N-body integrations show that the evolution of dynamically unstable systems does not lead to tidal migration but rather to planet ejections and collisions with the host star. This and other arguments lead us to suggest that most of the observed planets with a companion could not have been transported from further out through secular migration processes. More generally, by using a combination of numerical and analytic techniques we show that the high-e Lidov-Kozai migration scenario can only account for less than 10% of all gas giants observed between 0.1-1 AU. Simulations of multi-planet systems support this result. Our study indicates that rather than starting on highly eccentric orbits with orbital periods above one year, these "warm" Jupiters are more likely to have reached the region where they are observed today without having experienced significant tidal dissipation.

Large Kuiper Belt Objects are conventionally thought to have formed out of a massive planetesimal belt that is a few thousand times its current mass. Such a picture, however, is incompatible with multiple lines of evidence. Here, we present a new model for the conglomeration of Cold Classical Kuiper belt objects, out of a solid belt only a few times its current mass, or a few percent of the solid density in a Minimum Mass Solar Nebula. This is made possible by depositing most of the primordial mass in grains of size centimetre or smaller. These grains collide frequently and maintain a dynamically cold belt out of which large bodies grow efficiently: an order-unity fraction of the solid mass can be converted into large bodies, in contrast to the ~0.1% efficiency in conventional models. Such a light belt may represent the true outer edge of the Solar system, and it may have effectively halted the outward migration of Neptune. In addition to the high efficiency, our model can also produce a mass spectrum that peaks at an intermediate size, similar to the observed Cold Classicals, if one includes the effect of cratering collisions. In particular, the observed power-law break observed at ~30 km for Cold Classicals, one that has been interpreted as a result of collisional erosion, may be primordial in origin.

We develop and apply methods to extract planet masses and eccentricities from observed transit time variations (TTVs). First, we derive simple analytic expressions for the TTV that include the effects of both first- and second-order resonances. Second, we use N-body Markov chain Monte Carlo (MCMC) simulations, as well as the analytic formulae, to measure the masses and eccentricities of ten planets discovered by Kepler that have not previously been analyzed. Most of the ten planets have low densities. Using the analytic expressions to partially circumvent degeneracies, we measure small eccentricities of a few percent or less.

We study the efficiency of forming large bodies, starting from a sea of equal-sized planetesimals. This is likely one of the earlier steps of planet formation and relevant for the formation of the asteroid belt, the Kuiper belt and extra-solar debris disks. Here we consider the case that the seed planetesimals do not collide frequently enough for dynamical collisional to be important (the collisionless limit), using a newly constructed conglomeration code, and by carefully comparing numerical results with analytical scalings. In the absence of collisional cooling, as large bodies grow by accreting small bodies, the velocity dispersion of the small bodies ($u$) is increasingly excited. Growth passes from the well-known run-away stage (when $u$ is higher than the big bodies' hill velocity) to the newly discovered trans-hill stage (when $u$ and big bodies both grow, but $u$ remains at the big bodies' hill velocity). We find, concurring with the analytical understandings developed in Lithwick (2014), as well as previous numerical studies, that a size spectrum $dn/dR \propto R^{-4}$ results, and that the formation efficiency, defined as mass fraction in bodies much larger than the initial size, is $\sim {\mathrm{ a \,\, few}}\times R_\odot/a$, or $\sim 10^{-3}$ at the distance of the Kuiper belt. We argue that this extreme inefficiency invalidates the conventional conglomeration model for the formation of both our Kuiper belt and extra-solar debris disks. New theories, possibly involving direct gravitational collapse, or strong collisional cooling of small planetesimals, are required.

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield $0.75\pm 0.06$. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.

A transiting planet exhibits sinusoidal transit-time-variations (TTVs) if perturbed by a companion near a mean-motion-resonance (MMR). We search for sinusoidal TTVs in more than 2600 Kepler candidates, using the publicly available Kepler light-curves (Q0-Q12). We find that the TTV fractions rise strikingly with the transit multiplicity. Systems where four or more planets transit enjoy roughly five times higher TTV fraction than those where a single planet transits, and about twice higher than those for doubles and triples. In contrast, models in which all transiting planets arise from similar dynamical configurations predict comparable TTV fractions among these different systems. One simple explanation for our results is that there are at least two different classes of Kepler systems, one closely packed and one more sparsely populated.

We extract densities and eccentricities of 139 sub-Jovian planets by analyzing transit time variations (TTVs) obtained by the Kepler mission through Quarter 12. We partially circumvent the degeneracies that plague TTV inversion with the help of an analytical formula for the TTV. From the observed TTV phases, we find that most of these planets have eccentricities of order a few percent. More precisely, the r.m.s. eccentricity is 0.018^{+0.005}_{-0.004}, and planets smaller than 2.5R_\earth are around twice as eccentric as those bigger than 2.5R_\earth. We also find a best-fit density-radius relationship \rho ~3 g/cm^3 \times (R/3R_\earth)^{-2.3} for the 56 planets that likely have small eccentricity and hence small statistical correction to their masses. Many planets larger than 2.5R_\earth are less dense than water, implying that their radii are largely set by a massive hydrogen atmosphere.

In the inner solar system, the planets' orbits evolve chaotically, driven primarily by secular chaos. Mercury has a particularly chaotic orbit, and is in danger of being lost within a few billion years. Just as secular chaos is reorganizing the solar system today, so it has likely helped organize it in the past. We suggest that extrasolar planetary systems are also organized to a large extent by secular chaos. A hot Jupiter could be the end state of a secularly chaotic planetary system reminiscent of the solar system. But in the case of the hot Jupiter, the innermost planet was Jupiter- (rather than Mercury-) sized, and its chaotic evolution was terminated when it was tidally captured by its star. In this contribution, we review our recent work elucidating the physics of secular chaos and applying it to Mercury and to hot Jupiters. We also present new results comparing the inclinations of hot Jupiters thus produced with observations.

Tidally distorted rotating stars and gaseous planets are subject to a well-known linear fluid instability -- the elliptical instability. It has been proposed that this instability might drive enough energy dissipation to solve the long-standing problem of the origin of tidal dissipation in stars and planets. But the nonlinear outcome of the elliptical instability has yet to be investigated in the parameter regime of interest, and the resulting turbulent energy dissipation has not yet been quantified. We do so by performing three dimensional hydrodynamical simulations of a small patch of a tidally deformed fluid planet or star subject to the elliptical instability. We show that when the tidal deformation is weak, the nonlinear outcome of the instability leads to the formation of long-lived columnar vortices aligned with the axis of rotation. These vortices shut off the elliptical instability, and the net result is insufficient energy dissipation to account for tidal dissipation. However, further work is required to account for effects neglected here, including magnetic fields, turbulent convection, and realistic boundary conditions.

We investigate whether the elliptical instability is important for tidal dissipation in gaseous planets and stars. In a companion paper, we found that the conventional elliptical instability results in insufficient dissipation because it produces long-lived vortices that then quench further instability. Here, we study whether the addition of a magnetic field prevents those vortices from forming, and hence leads to enhanced dissipation. We present results from magnetohydrodynamical simulations that evolve the elliptical instability in a local patch of a rotating planet or star, in the presence of a weak magnetic field. We find that magnetic fields do indeed prevent vortices from forming, and hence greatly enhance the steady state dissipation rate. In addition, the resulting turbulence acts as a small-scale dynamo, amplifying the initially weak field. The inferred tidal dissipation is potentially important at short orbital periods. For example, it can circularise hot Jupiters with orbital periods shorter than 2.5 days, and synchronise their spins with their orbits out to 6 days. However, it appears unable to account for the hot Jupiters that appear to have been circularised out to six to ten day orbital periods. It also cannot account for the inferred circularisation of many close binary stars.

We analyze the transit timing variations obtained by the Kepler mission for 22 sub-jovian planet pairs (17 published, 5 new) that lie close to mean motion resonances. We find that the TTV phases for most of these pairs lie close to zero, consistent with an eccentricity distribution that has a very low RMS value of e ~ 0.01; but about a quarter of the pairs possess much higher eccentricities, up to 0.1 - 0.4. For the low-eccentricity pairs, we are able to statistically remove the effect of eccentricity to obtain planet masses from TTV data. These masses, together with those measured by radial velocity, yield a best fit mass-radius relation M~3 M_E (R/R_E). This corresponds to a constant surface escape velocity of 20km/s. We separate the planets into two distinct groups, "mid-sized" (those greater than 3 R_E), and "compact" (those smaller). All mid-sized planets are found to be less dense than water and therefore contain extensive H/He envelopes, likely comparable in mass to that of their cores. We argue that these planets have been significantly sculpted by photoevaporation. Surprisingly, mid-sized planets are discovered exclusively around stars more massive than 0.8 M_sun. The compact planets, on the other hand, are often denser than water. Combining our density measurements with those from radial velocity studies, we find that hotter compact planets tend to be denser, with the hottest ones reaching rock density. Moreover, hotter planets tend to be smaller in sizes. These results can be explained if the compact planets are made of rocky cores overlaid with a small amount of hydrogen, < 1% in mass, with water contributing little to their masses or sizes. Photoevaporation has exposed bare rocky cores in cases of the hottest planets. Our conclusion that these planets are likely not water-worlds contrasts with some previous studies.

When planetesimals begin to grow by coagulation, they enter an epoch of runaway, during which the biggest bodies grow faster than all the others. The questions of how runaway ends and what comes next have not been answered satisfactorily. Here we show that runaway is followed by a `trans-hill stage' that commences once the bodies become trans-hill, i.e. once the Hill velocity of the bodies that dominate viscous stirring matches the random speed of the small bodies they accrete. Subsequently, the small bodies' speed grows in lockstep with the big bodies' sizes, such that the bodies remain in the trans-hill state. Trans-hill growth is crucial for determining the efficiency of growing big bodies, as well as their growth timescale and size spectrum. We work out the properties of trans-hill growth analytically and confirm these numerically. Trans-hill growth has two sub-stages. In the earlier one, collisional cooling is irrelevant, in which case the efficiency of forming big bodies remains very low (0.1% in the Kuiper belt) and their mass spectrum is flat. This explains results from previous coagulation simulations for both the Kuiper belt and the asteroid belt. The second sub-stage commences when small bodies begin to collide with one another. Collisional cooling then controls the evolution, in which case the efficiency of forming big bodies rises, and their size spectrum becomes more top-heavy. Trans-hill growth terminates in one of two ways, depending on parameters. First, mutual accretion of big bodies can become significant and conglomeration proceeds until half the total mass is converted into big bodies. This mode of growth may explain the size distributions of minor bodies in the Solar System, and is explored in forthcoming work. Second, if big bodies become separated by their Hill radius, oligarchy commences. This mode likely precedes the formation of fully-fledged planets.

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

We study the radius evolution of close-in extra-solar jupiters under Ohmic heating, a mechanism that was recently proposed to explain the large observed sizes of many of these planets. Planets are born with high entropy and they subsequently cool and contract. We focus on two cases: first, that ohmic heating commences when the planet is hot (high entropy); and second, that it commences after the planet has cooled. In the former case, we use analytical scalings and numerical experiments to confirm that Ohmic heating is capable of suspending the cooling as long as a few percent of the stellar irradiation is converted into Ohmic heating, and the planet has a surface wind that extends to pressures of ~10 bar or deeper. For these parameters, the radii at which cooling is stalled are consistent with (or larger than) the observed radii of most planets. The only two exceptions are WASP-17b and HAT-P-32b. In contrast to the high entropy case, we show that Ohmic heating cannot significantly re-inflate planets after they have already cooled. This leads us to suggest that the diversity of radii observed in hot jupiters may be partially explained by the different epochs at which they are migrated to their current locations.

Most planet pairs in the Kepler data that have measured transit time variations (TTV) are near first-order mean-motion resonances. We derive analytical formulae for their TTV signals. We separate planet eccentricity into free and forced parts, where the forced part is purely due to the planets' proximity to resonance. This separation yields simple analytical formulae. The phase of the TTV depends sensitively on the presence of free eccentricity: if the free eccentricity vanishes, the TTV will be in phase with the longitude of conjunctions. This effect is easily detectable in current TTV data. The amplitude of the TTV depends on planet mass and free eccentricity, and it determines planet mass uniquely only when the free eccentricity is sufficiently small. We proceed to analyze the TTV signals of six short period Kepler pairs. We find that three of these pairs (Kepler-18,24,25) have TTV phase consistent with zero. The other three (Kepler-23,28,32) have small TTV phases, but ones that are distinctly non-zero. We deduce that the free eccentricities of the planets are small, < 0.01, but not always vanishing. Furthermore, as a consequence of this, we deduce that the true masses of the planets are fairly accurately determined by the TTV amplitudes, within a factor < 2. The smallness of the free eccentricities suggests that the planets have experienced substantial dissipation. This is consistent with the hypothesis that the observed pile-up of Kepler pairs near mean-motion resonances is caused by resonant repulsion. But the fact that some of the planets have non-vanishing free eccentricity suggests that after resonant repulsion occurred there was a subsequent phase in the planets' evolution when their eccentricities were modestly excited, perhaps by interplanetary interactions.

Planetary systems discovered by the Kepler space telescope exhibit an intriguing feature. While the period ratios of adjacent low-mass planets appear largely random, there is a significant excess of pairs that lie just wide of resonances and a deficit on the near side. We demonstrate that this feature naturally arises when two near-resonant planets interact in the presence of weak dissipation that damps eccentricities. The two planets repel each other as orbital energy is lost to heat. This moves near-resonant pairs just beyond resonance, by a distance that reflects the integrated dissipation they experienced over their lifetimes. We find that the observed distances may be explained by tides if tidal dissipation is unexpectedly efficient (tidal quality factor ~10). Once the effect of resonant repulsion is accounted for, the initial orbits of these low mass planets show little preference for resonances. This could constrain their origin.

Modelling the formation of super-km-sized planetesimals by gravitational collapse of regions overdense in small particles requires numerical algorithms capable of handling simultaneously hydrodynamics, particle dynamics and particle collisions. While the initial phases of radial contraction are dictated by drag forces and gravity, particle collisions become gradually more significant as filaments contract beyond Roche density. Here we present a new numerical algorithm for treating momentum and energy exchange in collisions between numerical superparticles representing a high number of physical particles. We adopt a Monte Carlo approach where superparticle pairs in a grid cell collide statistically on the physical collision time-scale. Collisions occur by enlarging particles until they touch and solving for the collision outcome, accounting for energy dissipation in inelastic collisions. We demonstrate that superparticle collisions can be consistently implemented at a modest computational cost. In protoplanetary disc turbulence driven by the streaming instability, we argue that the relative Keplerian shear velocity should be subtracted during the collision calculation. If it is not subtracted, density inhomogeneities are too rapidly diffused away, as bloated particles exaggerate collision speeds. Local particle densities reach several thousand times the mid-plane gas density. We find efficient formation of gravitationally bound clumps, with a range of masses corresponding to contracted radii from 100 to 400 km when applied to the asteroid belt and 150 to 730 km when applied to the Kuiper belt, extrapolated using a constant self-gravity parameter. The smaller planetesimals are not observed at low resolution, but the masses of the largest planetesimals are relatively independent of resolution and treatment of collisions.

We investigate the collapse and internal structure of dark matter halos. We consider halo formation from initially scale-free perturbations, for which gravitational collapse is self-similar. Fillmore and Goldreich (1984) and Bertschinger (1985) solved the one dimensional (i.e. spherically symmetric) case. We generalize their results by formulating the three dimensional self-similar equations. We solve the equations numerically and analyze the similarity solutions in detail, focusing on the internal density profiles of the collapsed halos. By decomposing the total density into subprofiles of particles that collapse coevally, we identify two effects as the main determinants of the internal density structure of halos: adiabatic contraction and the shape of a subprofile shortly after collapse; the latter largely reflects the triaxiality of the subprofile. We develop a simple model that describes the results of our 3D simulations. In a companion paper, we apply this model to more realistic cosmological fluctuations, and thereby explain the origin of the nearly universal (NFW-like) density profiles found in N-body simulations.

We study the dynamical evolution of a test particle that orbits a star in the presence of an exterior massive planet, considering octupole-order secular interactions. In the standard Kozai mechanism (SKM), the planet's orbit is circular, and so the particle conserves vertical angular momentum. As a result, the particle's orbit oscillates periodically, exchanging eccentricity for inclination. However, when the planet's orbit is eccentric, the particle's vertical angular momentum varies and its Kozai oscillations are modulated on longer timescales---we call this the eccentric Kozai mechanism (EKM). The EKM can lead to behavior that is dramatically different from the SKM. In particular, the particle's orbit can flip from prograde to retrograde and back again, and it can reach arbitrarily high eccentricities given enough time. We map out the conditions under which this dramatic behavior (flipping and extreme eccentricities) occurs, and show that when the planet's eccentricity is sufficiently high, it occurs quite generically. For example, when the planet's eccentricity exceeds a few percent of the ratio of semi-major axes (outer to inner), around half of randomly oriented test particle orbits will flip and reach extreme eccentricities. The SKM has often been invoked for bringing pairs of astronomical bodies (star-star, planet-star, compact-object pairs) close together. Including the effect of the EKM will enhance the rate at which such matchmaking occurs.

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

We study the chaotic orbital evolution of planetary systems, focusing on secular (i.e., orbit-averaged) interactions, because these often dominate on long timescales. We first focus on the evolution of a test particle that is forced by multiple massive planets. To linear order in eccentricity and inclination, its orbit precesses with constant frequencies. But nonlinearities can shift the frequencies into and out of secular resonance with the planets' eigenfrequencies, or with linear combinations of those frequencies. The overlap of these nonlinear secular resonances drive secular chaos in planetary systems. We quantify the resulting dynamics for the first time by calculating the locations and widths of nonlinear secular resonances. When results from both analytical calculations and numerical integrations are displayed together in a newly developed "map of the mean momenta" (MMM), the agreement is excellent. This map is particularly revealing for non-coplanar planetary systems and demonstrates graphically that chaos emerges from overlapping secular resonances. We then apply this newfound understanding to Mercury. Previous numerical simulations have established that Mercury's orbit is chaotic, and that Mercury might even collide with Venus or the Sun. We show that Mercury's chaos is primarily caused by the overlap between resonances that are combinations of four modes, the Jupiter-dominated eccentricity mode, the Venus-dominated inclination mode and Mercury's free eccentricity and inclination. Numerical integration of the Solar system confirms that a slew of these resonant angles alternately librate and circulate. We are able to calculate the threshold for Mercury to become chaotic: Jupiter and Venus must have eccentricity and inclination of a few percent. Mercury appears to be perched on the threshold for chaos.

In a planetary system with two or more well-spaced, eccentric, inclined planets, secular interactions may lead to chaos. The innermost planet may gradually become very eccentric and/or inclined, as a result of the secular degrees of freedom drifting towards equipartition of angular momentum deficit. Secular chaos is known to be responsible for the eventual destabilization of Mercury in our own Solar System. Here we focus on systems with three giant planets. We characterize the secular chaos and demonstrate the criterion for it to occur, but leave a detailed understanding of secular chaos to a companion paper (Lithwick & Wu, 2010). After an extended period of eccentricity diffusion, the inner planet's pericentre can approach the star to within a few stellar radii. Strong tidal interactions and ensuing tidal dissipation extracts orbital energy from the planet and pulls it inward, creating a hot Jupiter. In contrast to other proposed channels for the production of hot Jupiters, such a scenario (which we term "secular migration") explains a range of observations: the pile-up of hot Jupiters at 3-day orbital periods, the fact that hot Jupiters are in general less massive than other RV planets, that they may have misaligned inclinations with respect to stellar spin, and that they have few easily detectable companions (but may have giant companions in distant orbits). Secular migration can also explain close-in planets as low in mass as Neptune; and an aborted secular migration can explain the "warm Jupiters" at intermediate distances. In addition, the frequency of hot Jupiters formed via secular migration increases with stellar age. We further suggest that secular chaos may be responsible for the observed eccentricities of giant planets at larger distances, and that these planets could exhibit significant spin-orbit misalignment.

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics.

Two dimensional hydrodynamical disks are nonlinearly unstable to the formation of vortices. Once formed, these vortices essentially survive forever. What happens in three dimensions? We show with pseudospectral simulations that in 3D a vortex in a short box forms and survives just as in 2D. But a vortex in a tall box is unstable and is destroyed. In our simulation, the unstable vortex decays into a transient turbulent-like state that transports angular momentum outward at a nearly constant rate for hundreds of orbital times. The 3D instability that destroys vortices is a generalization of the 2D instability that forms them. We derive the conditions for these nonlinear instabilities to act by calculating the coupling between linear modes, and thereby derive the criterion for a vortex to survive in 3D as it does in 2D: the azimuthal extent of the vortex must be larger than the scale height of the accretion disk. When this criterion is violated, the vortex is unstable and decays. Because vortices are longer in azimuthal than in radial extent by a factor that is inversely proportional to their excess vorticity, a vortex with given radial extent will only survive in a 3D disk if it is sufficiently weak. This counterintuitive result explains why previous 3D simulations always yielded decaying vortices: their vortices were too strong. We conclude that in protoplanetary disks weak vortices can trap dust and serve as the nurseries of planet formation. Decaying strong vortices might be responsible for the outwards transport of angular momentum that is required to make accretion disks accrete.

Pluto's recently discovered minor moons, Nix and Hydra, have almost circular orbits, and are nearly coplanar with Charon, Pluto's major moon. This is surprising because tidal interactions with Pluto are too weak to damp their eccentricities. We consider an alternative possibility: that Nix and Hydra circularize their orbits by exciting Charon's eccentricity via secular interactions, and Charon in turn damps its own eccentricity by tidal interaction with Pluto. The timescale for this process can be less than the age of the Solar System, for plausible tidal parameters and moon masses. However, as we show numerically and analytically, the effects of the 2:1 and 3:1 resonant forcing terms between Nix and Charon complicate this picture. In the presence of Charon's tidal damping, the 2:1 term forces Nix to migrate outward and the 3:1 term changes the eccentricity damping rate, sometimes leading to eccentricity growth. We conclude that this mechanism probably does not explain Nix and Hydra's current orbits. Instead, we suggest that they were formed in-situ with low eccentricities. We also show that an upper limit on Nix's migration speed sets a lower limit on Pluto-Charon's tidal circularization timescale of >10^5 yrs. Moreover, Hydra's observed proper eccentricity may be explained by the 3:2 forcing by Nix.

How did Pluto's recently discovered minor moons form? Ward and Canup propose an elegant solution in which Nix and Hydra formed in the collision that produced Charon, then were caught into corotation resonances with Charon, and finally were transported to their current location as Charon migrated outwards. We show with numerical integrations that, if Charon's eccentricity is judiciously chosen, this scenario works beautifully for either Nix or Hydra. However, it cannot work for both Nix and Hydra simultaneously. To transport Nix, Charon's eccentricity must satisfy e_C< 0.024; otherwise, the second order Lindblad resonance at 4:1 overlaps with the corotation resonance, leading to chaos. To transport Hydra, e_C > 0.7 R_p/a_C > 0.04; otherwise migration would be faster than libration, and Hydra would slip out of resonance. These two restrictions conflict. Having ruled out this scenario, we suggest an alternative: that many small bodies were captured from the nebular disk, and they were responsible for forming, migrating and damping Nix and Hydra. If this is true, small moons could be common around large Kuiper belt objects.

We examine how perturbed shear flows evolve in two-dimensional, incompressible, inviscid hydrodynamical fluids, with the ultimate goal of understanding the dynamics of accretion disks. To linear order, vorticity waves are swung around by the background shear, and their velocities are amplified transiently before decaying. It has been speculated that sufficiently amplified modes might couple nonlinearly, leading to turbulence. Here we show how nonlinear coupling occurs in two dimensions. This coupling is remarkably simple because it only lasts for a short time interval, when one of the coupled modes is in mid-swing. We focus on the interaction between a swinging and an axisymmetric mode. There is instability provided that k_{y,swing}/k_{x,axi} < omega/q, i.e., that the ratio of wavenumbers is less than the ratio of the axisymmetric mode's vorticity to the background vorticity. If this is the case, then when the swinging mode is in mid-swing it couples with the axisymmetric mode to produce a new leading swinging mode that has larger vorticity than itself; this new mode in turn produces an even larger leading mode, etc. Therefore all axisymmetric modes, regardless of how small in amplitude, are unstable to perturbations with sufficiently large azimuthal wavelength. We show that this shear instability occurs whenever the momentum transported by a perturbation has the sign required for it to diminish the background shear; only when this occurs can energy be extracted from the mean flow and hence added to the perturbation. For an accretion disk, this means that the instability transports angular momentum outwards while it operates.

We describe the diffusion and random velocities of solid particles due to stochastic forcing by turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas to particle diffusivity) is Sc = 1 + (Omega t_stop)^2, in terms of the particle stopping time, t_stop, and the orbital frequency, Omega. The standard result, Sc = 1 + t_stop/t_eddy, in terms of the eddy turnover time, t_eddy, is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with stopping time. Consequently, particles larger than ~ 10 cm in protoplanetary disks will suffer less radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients will help interpret numerical simulations. We confirm previous results for the vertical stirring of particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to orbital times. The particle layer becomes thinner for t_eddy > 1/Omega, with the strength of turbulent diffusion held fixed. We use two analytic techniques -- the Hinze-Tchen formalism and the Fokker-Planck equation with velocity diffusion -- with identical results when the regimes of validity overlap. We include simple physical arguments for the scaling of our results.

The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects (KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations from the solar neighborhood. We critically review the present-day belt's observed properties and the theories designed to explain them. Theories are organized according to a possible time-line of events. In chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by Neptune-mass planets (``oligarchs''), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) collisional comminution of the smallest KBOs. Recent work underscores how small, collisional, primordial planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass planets whose orbits cross into the Kuiper belt for up to ~20 Myr help generate the high-perihelion members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made of how little primordial mass might reside in bodies having sizes of order 100 km. These and other frontier issues in trans-Neptunian space are discussed quantitatively.

We present new high-resolution observations of Sagittarius A* at wavelengths of 17.4 to 23.8 cm with the Very Large Array in A configuration with the Pie Town Very Long Baseline Array antenna. We use the measured sizes to calibrate the interstellar scattering law and find that the major axis size of the scattering law is smaller by ~6% than previous estimates. Using the new scattering law, we are able to determine the intrinsic size of Sgr A* at wavelengths from 0.35 cm to 6 cm using existing results from the VLBA. The new law increases the intrinsic size at 0.7 cm by ~20% and <5% at 0.35 cm. The intrinsic size is 13^{+7}_{-3} Schwarzschild radii at 0.35 cm and is proportional to lambda^gamma, where gamma is in the range 1.3 to 1.7.

We present a phenomenological model of imbalanced MHD turbulence in an incompressible magnetofluid. The steady-state cascades, of waves traveling in opposite directions along the mean magnetic field, carry unequal energy fluxes to small length scales, where they decay due to viscous and resistive dissipation. The inertial-range scalings are well-understood when both cascades are weak. We study the case when both cascades are, in a sense, strong. The inertial-range of this imbalanced cascade has the following properties: (i) the ratio of the r.m.s. Elsasser amplitudes is independent of scale, and is equal to the ratio of the corresponding energy fluxes; (ii) in common with the balanced strong cascade, the energy spectra of both Elsasser waves are of the anisotropic Kolmogorov form, with their parallel correlation lengths equal to each other on all scales, and proportional to the two-thirds power of the transverse correlation length; (iii) the equality of cascade time and waveperiod (critical balance) that characterizes the strong balanced cascade does not apply to the Elsasser field with the larger amplitude. Instead, the more general criterion that always applies to both Elsasser fields is that the cascade time is equal to the correlation time of the straining imposed by oppositely-directed waves. Our results are particularly relevant for turbulence in the solar wind. Spacecraft measurements have established that, in the inertial range of solar wind turbulence, waves travelling away from the sun have higher amplitudes than those travelling towards it. Result (i) allows us to infer the turbulent flux ratios from the amplitude ratios, thus providing insight into the origin of the turbulence.

We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and can be used to simulate, for the first time, the interaction of planets with disks over many viscous times. Though the algorithm is implemented in two dimensions-i.e., the motions of bodies need only be tracked in a plane-it captures the behavior of fully three-dimensional disks in which collisions maintain inclinations that are comparable to random eccentricities. We subject the algorithm to a battery of tests for the case of an isolated, narrow, circular ring. Numerical simulations agree with analytic theory with regards to how particles' random velocities equilibrate; how the ring viscously spreads; and how energy dissipation, angular momentum transport, and material transport are connected. We derive and measure the critical value of the coefficient of restitution above which viscous stirring dominates inelastic damping and the particles' velocity dispersion runs away.

We report the discovery of a transient radio source 2.7 arcsec (0.1 pc projected distance) South of the Galactic Center massive black hole, Sagittarius A*. The source flared with a peak of at least 80 mJy in March 2004. The source was resolved by the Very Large Array into two components with a separation of ~0.7 arcsec and characteristic sizes of ~0.2 arcsec. The two components of the source faded with a power-law index of 1.1 +/- 0.1. We detect an upper limit to the proper motion of the Eastern component of ~3 x 10^3 km s^-1 relative to Sgr A*. We detect a proper motion of ~10^4 km s^-1 for the Western component relative to Sgr A*. The transient was also detected at X-ray wavelengths with the Chandra X-ray Observatory and the XMM-Newton telescope and given the designation CXOGC J174540.0-290031. The X-ray source falls in between the two radio components. The maximum luminosity of the X-ray source is ~10^36 erg s^-1, significantly sub-Eddington. The radio jet flux density predicted by the X-ray/radio correlation for X-ray binaries is orders of magnitude less than the measured flux density. We conclude that the radio transient is the result of a bipolar jet originating in a single impulsive event from the X-ray source and interacting with the dense interstellar medium of the Galactic Center.

The problem of accretion in the Trojan 1:1 resonance is akin to the standard problem of planet formation, transplanted from a star-centered disk to a disk centered on the Lagrange point. The newly discovered class of Neptune Trojans promises to test theories of planet formation by coagulation. Neptune Trojans resembling the prototype 2001 QR322 (``QR'')--whose radius of ~100 km is comparable to that of the largest Jupiter Trojan--may outnumber their Jovian counterparts by a factor of ~10. We discover that seeding the 1:1 resonance with debris from planetesimal collisions and having the seed particles accrete in situ naturally reproduces the inferred number of QR-sized Trojans. We analyze accretion in the Trojan sub-disk by applying the two-groups method, accounting for kinematics specific to the resonance. We find that a Trojan sub-disk comprising decimeter-sized seed particles and having a surface density 1e-3 that of the local minimum-mass disk produces ~10 QR-sized objects in ~1 Gyr, in accord with observation. Further growth is halted by collisional diffusion of seed particles out of resonance. In our picture, the number and sizes of the largest Neptune Trojans represent the unadulterated outcome of dispersion-dominated, oligarchic accretion. Large Neptune Trojans, perhaps the most newly accreted objects in our Solar System, may today have a dispersion in orbital inclination of less than ~10 degrees, despite the existence of niches of stability at higher inclinations. Such a vertically thin disk, born of a dynamically cold environment necessary for accretion, and raised in minimal contact with external perturbations, contrasts with the thick disks of other minor body belts.

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.

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.

MHD turbulence consists of waves that propagate along magnetic fieldlines, in both directions. When two oppositely directed waves collide, they distort each other, without changing their respective energies. In weak MHD turbulence, a given wave suffers many collisions before cascading. "Imbalance" means that more energy is going in one direction than the other. In general, MHD turbulence is imbalanced. A number of complications arise for the imbalanced cascade that are unimportant for the balanced one. We solve weak MHD turbulence that is imbalanced. Of crucial importance is that the energies going in both directions are forced to equalize at the dissipation scale. We call this the "pinning" of the energy spectra. It affects the entire inertial range. Weak MHD turbulence is particularly interesting because perturbation theory is applicable. Hence it can be described with a simple kinetic equation. Galtier et al. (2000) derived this kinetic equation. We present a simpler, more physical derivation, based on the picture of colliding wavepackets. In the process, we clarify the role of the zero-frequency mode. We also explain why Goldreich & Sridhar claimed that perturbation theory is inapplicable, and why this claim is wrong. (Our "weak" is equivalent to Goldreich & Sridhar's "intermediate.") We perform numerical simulations of the kinetic equation to verify our claims. We construct simplified model equations that illustrate the main effects. Finally, we show that a large magnetic Prandtl number does not have a significant effect, and that hyperviscosity leads to a pronounced bottleneck effect.

Radio-wave scintillation observations reveal a nearly Kolmogorov spectrum of density fluctuations in the ionized interstellar medium. Although this density spectrum is suggestive of turbulence, no theory relevant to its interpretation exists. We calculate the density spectrum in turbulent magnetized plasmas by extending the theory of incompressible MHD turbulence given by Goldreich & Sridhar to include the effects of compressibility and particle transport. Our most important results are as follows. (1) Density fluctuations are due to the slow mode and the entropy mode. Both modes are passively mixed by the cascade of shear Alfven waves. Since the shear Alfven waves have a Kolmogorov spectrum, so do the density fluctuations. (2) Observed density fluctuation amplitudes imply either that the magnetic and gas pressures are comparable, or that the outer scale of the turbulence is very small. (3) A high degree of ionization is required for the cascade to survive damping by neutrals and thereby to extend to small lengthscales. Regions that are insufficiently ionized produce density fluctuations only on lengthscales larger than the neutral damping scale. These regions may account for the excess of power that is found on large scales. (4) Both the entropy mode and the slow mode are damped on lengthscales below that at which protons can diffuse across an eddy during the eddy's turnover time. Consequently, eddies whose extents along the magnetic field are smaller than the proton collisional mean free path do not contribute to the density spectrum. However, in MHD turbulence eddies are highly elongated along the magnetic field. From an observational perspective, the relevant lengthscale is that transverse to the magnetic field. Thus the cut-off lengthscale for density fluctuations is significantly smaller than the proton mean free path.

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.

Extremely strong magnetic fields change the vacuum index of refraction. This induces a lensing effect that is not unlike the lensing phenomenon in strong gravitational fields. The main difference between the two is the polarization dependency of the magnetic lensing, a behaviour that induces a handful of interesting effects. The main prediction is that the thermal emission of neutron stars with extremely strong magnetic fields is polarized - up to a few percent for the largest fields known. This potentially allows a direct method for measuring their magnetic fields.