As established by ground based surveys, white dwarfs with hydrogen atmospheres pulsate as they cool across the temperature range, $12500\,\mathrm{K} \gtrsim T_{\mathrm{eff}} \gtrsim 10800\,\mathrm{K}$. Known as DAVs or ZZ Ceti stars, their oscillations are attributed to overstable g-modes excited by convective driving. The effective temperature at the blue edge of the instability strip is slightly lower than that at which a surface convection zone appears. The temperature at the red edge is a two-decade old puzzle. Recently, {\it Kepler} discovered a number of cool DAVs which pulsate at higher frequencies and with much smaller photometric amplitudes than expected based on trends extrapolated from DAVs found by ground based observations. Remarkably, some of them exhibit sporadic outbursts separated by days, each lasting several hours, and releasing $\sim 10^{33}-10^{34}\,\mathrm{erg}$. We provide quantitative explanations for both the red edge and the outbursts. The minimal frequency for overstable modes rises abruptly near the red edge. Although high frequency overstable modes exist below the red edge, their photometric amplitudes are generally too small to be detected by ground based observations. Nevertheless, these overstable parent modes can manifest themselves through nonlinear mode couplings to damped daughter modes which generate limit cycles giving rise to photometric outbursts.

The properties of supersonic isothermal turbulence influence a variety of astrophysical phenomena, including the structure and evolution of star forming clouds. This work presents a simple model for the structure of dense regions in turbulence in which the density distribution behind isothermal shocks originates from rough hydrostatic balance between the pressure gradient behind the shock and its deceleration from ram pressure applied by the background fluid. Using simulations of supersonic isothermal turbulence and idealized waves moving through a background medium, we show that the structural properties of dense, shocked regions broadly agree with our analytical model. Our work provides a new conceptual picture for describing the dense regions, which complements theoretical efforts to understand the bulk statistical properties of turbulence and attempts to model the more complex features of star forming clouds like magnetic fields, self-gravity, or radiative properties.

Rubble piles are a common feature of solar system bodies. They are composed of monolithic elements of ice or rock bound by gravity. Voids occupy a significant fraction of the volume of a rubble pile. They can exist up to pressure $P\approx \epsy\mu$, where $\epsy$ is the monolithic material's yield strain and $\mu$ its rigidity. At low $P$, contacts between neighboring elements are confined to a small fraction of their surface areas. As a result, the effective thermal conductivity of a rubble pile, $\kcon\approx k(P/(\epsy\mu))^{1/2}$, can be orders of magnitude smaller than, $k$, the thermal conductivity of its monolithic elements. In a fluid-free environment, only radiation can transfer energy across voids. It contributes an additional component, $\krad=16\ell\sigma T^3/3$, to the total effective conductivity, $\keff=\kcon +\krad$. Here $\ell$, the inverse of the opacity per unit volume, is of order the size of the elements and voids. An important distinction between $\kcon$ and $\krad$ is that the former is independent of the size of the elements whereas the latter is proportional to it. Our expression for $\keff$ provides a good fit to the depth dependence of thermal conductivity in the top $140\,\mathrm{cm}$ of the lunar regolith. It also offers a good starting point for detailed modeling of thermal inertias for asteroids and satellites. Measurement of the response of surface temperature to variable insolation is a valuable diagnostic of a regolith. There is an opportunity for careful experiments under controlled laboratory conditions to test models of thermal conductivity such as the one we outline.

The pulsar triple system, J0337+1715, is remarkably regular and highly hierarchical. Secular dynamics controls its long term evolution with orbital commensurability having negligible effect. Secular interactions transfer angular momentum between inner and outer orbits unless their apsidal lines are parallel or anti-parallel. These choices correspond, respectively, to orthogonal eigenmodes p and a. Each is characterized by an eccentricity ratio set by the masses and semi-major axes, i.e., $e_{p, 1}/e_{p, 2}\sim a_1/a_2$ while $e_{a, 1}/e_{a, 2}\sim (a_1/a_2)^{-3/2}(m_2/m_1)$. Mode p dominates the system's current state so $e_1/e_2$ always remains close to $e_{p, 1}/e_{p, 2}$. A small contribution by Mode a causes $e_1$ and $e_2$ to oscillate with period $\sim 10^3\,\yr$. Orbital changes should be apparent in a few years. These will reveal the forcing of the apsidal precession of the inner orbit by general relativity (GR), and possibly also the smaller contribution due to the tidal and rotational distortion of the inner white dwarf (WD). Phinney (1992) proposes that the epicyclic energy of a WD-pulsar binary reaches equipartition with the kinetic energy of a single convective eddy when the WD's progenitor fills its Roche lobe. We extend Phinney's theory to apply to modes rather than individual orbits. Thus we predict that Mode p and Mode a achieved equipartition with eddies in the giant envelopes of the progenitors of the outer and inner WD, respectively. The most effective eddies are those with lifetimes closest to the orbit period. These were more energetic in the progenitor of the outer WD. This explains why Mode p overwhelms Mode a, and also why the inner binary's orbit is far more eccentric than other WD-pulsar binaries with similar orbit periods. Mode a's small but finite amplitude places a lower bound of $Q\sim 10^6$ on the tidal quality parameter of the inner WD.

Fast Radio Bursts (FRBs) are isolated, $\ms$ radio pulses with dispersion measure (DM) of order $10^3\DMunit$. Galactic candidates for the DM of high latitude bursts detected at $\GHz$ frequencies are easily dismissed. DM from bursts emitted in stellar coronas are limited by free-free absorption and those from HII regions are bounded by the nondetection of associated free-free emission at radio wavelengths. Thus, if astronomical, FRBs are probably extra-galactic. FRB 110220 has a scattering tail of $\sim 5.6\pm 0.1 \ms$. If the electron density fluctuations arise from a turbulent cascade, the scattering is unlikely to be due to propagation through the diffuse intergalactic plasma. A more plausible explanation is that this burst sits in the central region of its host galaxy. Pulse durations of order $\ms$ constrain the sizes of FRB sources implying high brightness temperatures that indicates coherent emission. Electric fields near FRBs at cosmological distances would be so strong that they could accelerate free electrons from rest to relativistic energies in a single wave period.

We assess the multi-planet systems discovered by the Kepler satellite in terms of current ideas about orbital migration and eccentricity damping due to planet-disk interactions. Our primary focus is on mean motion resonances. Only a few percent of planet pairs are in close proximity to a resonance. However, predicted migration rates (parameterized by $\tau_n=n/{|\dot n|}$) imply that during convergent migration most planets would have been captured into first order resonances. Eccentricity damping (parameterized by $\tau_e=e/{|\dot e|}$) offers a plausible resolution. Estimates suggest $\tau_e/\tau_n\sim (h/a)^2\sim 10^{-2}$, where $h/a$ is the ratio of disk thickness to radius. Together, eccentricity damping and orbital migration give rise to an equilibrium eccentricity, $e_{eq}\sim(\tau_e/\tau_n)^{1/2}$. Capture is permanent provided $e_{eq}\lesssim \mu^{1/3}$, where $\mu$ denotes the planet to star mass ratio. But for $e_{eq}\gtrsim \mu^{1/3}$, capture is only temporary because librations around equilibrium are overstable and lead to passage through resonance on timescale $\tau_e$. Most Kepler planet pairs have $e_{eq}>\mu^{1/3}$. Since $\tau_n>> \tau_e$ is the timescale for migration between neighboring resonances, only a modest percentage of pairs end up trapped in resonances after the disk disappears. Planet pairs close to a mean motion resonance typically exhibit period ratios 1-2% larger than those for exact resonance. The direction of this shift undoubtedly reflects the same asymmetry that requires convergent migration for resonance capture. Permanent resonance capture at these separations from exact resonance would demand $\mu (\tau_n/\tau_e)^{1/2}\gtrsim 0.01$, a value that estimates of $\mu$ from transit data and $(\tau_e/\tau_n)^{1/2}$ from theory are insufficient to match. Plausible alternatives involve eccentricity damping during or after disk dispersal. (Abridged)

Turbulence influences the behavior of many astrophysical systems, frequently by providing non-thermal pressure support through random bulk motions. Although turbulence is commonly studied in systems with constant volume and mean density, turbulent astrophysical gases often expand or contract under the influence of pressure or gravity. Here, we examine the behavior of turbulence in contracting volumes using idealized models of compressed gases. Employing numerical simulations and an analytical model, we identify a simple mechanism by which the turbulent motions of contracting gases "adiabatically heat", experiencing an increase in their random bulk velocities until the largest eddies in the gas circulate over a "Hubble" time of the contraction. Adiabatic heating provides a mechanism for sustaining turbulence in gases where no large-scale driving exists. We describe this mechanism in detail and discuss some potential applications to turbulence in astrophysical settings.

A number of synchronous moons are thought to harbor water oceans beneath their outer ice shells. A subsurface ocean frictionally decouples the shell from the interior. This has led to proposals that a weak tidal or atmospheric torque might cause the shell to rotate differentially with respect to the synchronously rotating interior. As a result of centrifugal and tidal forces, the ocean would assume an ellipsoidal shape with its long axis aligned toward the parent planet. Any displacement of the shell away from its equilibrium position would induce strains thereby increasing its elastic energy and giving rise to an elastic restoring torque. We compare the elastic torque with the tidal torque acting on Europa and the atmospheric torque acting on Titan. For Europa, the tidal torque is far too weak to produce stresses that could fracture the ice shell, thus refuting a widely advocated idea. Instead, we suggest that cracks arise from time-dependent stresses due to non-hydrostatic gravity anomalies from tidally driven, episodic convection in the interior. Two years of Cassini RADAR observations of Titan's surface are interpreted as implying an angular displacement of ~0.24 degrees relative to synchroneity. Compatibility of the amplitude and phase of the observed non-synchronous rotation with estimates of the atmospheric torque requires that Titan's shell be decoupled from its interior. We find that the elastic torque balances the atmospheric torque at an angular displacement <0.05 degrees, thus coupling the shell to the interior. Moreover, if Titan's surface were spinning faster than synchronous, the tidal torque tending to restore synchronous rotation would certainly be larger than the atmospheric torque. There must either be a problem with the interpretation of the radar observations, or with our understanding of Titan's atmosphere and/or interior.

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 atmospheres of Jupiter and Saturn exhibit strong and stable zonal winds. How deep the winds penetrate unabated into each planet is unknown. Our investigation favors shallow winds. It consists of two parts. The first part makes use of an Ohmic constraint; Ohmic dissipation associated with the planet's magnetic field cannot exceed the planet's net luminosity. Application to Jupiter (J) and Saturn (S) shows that the observed zonal winds cannot penetrate below a depth at which the electrical conductivity is about six orders of magnitude smaller than its value at the molecular-metallic transition. Measured values of the electrical conductivity of molecular hydrogen yield radii of maximum penetration of 0.96R_J and 0.86R_S, with uncertainties of a few percent of R. At these radii, the magnetic Reynolds number based on the zonal wind velocity and the scale height of the magnetic diffusivity is of order unity. These limits are insensitive to difficulties in modeling turbulent convection. They permit complete penetration along cylinders of the equatorial jets observed in the atmospheres of Jupiter and Saturn. The second part investigates how deep the observed zonal winds actually do penetrate. Truncation of the winds in the planet's convective envelope would involve breaking the Taylor-Proudman constraint on cylindrical flow. This would require a suitable nonpotential acceleration which none of the obvious candidates appears able to provide. Accelerations arising from entropy gradients, magnetic stresses, and Reynolds stresses appear to be much too weak. These considerations suggest that strong zonal winds are confined to shallow, stably stratified layers, with equatorial jets being the possible exception.

Saturn's magnetic field is remarkably axisymmetric. Its dipole axis is inclined by less than 0.2 deg with respect to its rotation axis. Rotationally driven convection of magnetospheric plasma breaks the axisymmetry of its external magnetic field. Field aligned currents transfer angular momentum from the planet to a tongue of outflowing plasma. This transfer slows the rate of rotation of the ionosphere relative to that of the underlying atmosphere. The currents are the source for the non-axisymmetric components of the field. The common rotation rates of these components and Saturn's kilometric radio (SKR) bursts is that of the plasma near the orbit of Enceladus, and by extension the rotation rate in the ionosphere to which this plasma is coupled. That rate tells us nothing about the rotation rate of Saturn's deep interior. Of that we remain ignorant. Magnetic perturbations with magnitudes similar to those observed by Cassini are produced for Mdot ~ 10^4 g/s, a value similar to estimates for the rate of production of plasma from Saturn's E-ring. Enhancement of the SKR occurs in a narrow range of longitudes where the tip of the outgoing plasma stream connects to the auroral ionosphere via field lines that are bowed outwards by currents that supply the plasma's centripetal acceleration. (abridged)

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

A strong case has been made that radio waves from sources within about half a degree of the Galactic Center undergo extreme diffractive scattering. However, problems arise when standard (``Kolmogorov'') models of electron density fluctuations are employed to interpret the observations of scattering in conjunction with those of free-free radio emission. Specifically, the outer scale of a Kolmogorov spectrum of electron density fluctuations is constrained to be so small that it is difficult to identify an appropriate astronomical setting. Moreover, an unacceptably high turbulent heating rate results if the outer scale of the velocity field coincides with that of the density fluctuations. We propose an alternative model based on folded magnetic field structures that have been reported in numerical simulations of small-scale dynamos. Nearly isothermal density variations across thin current sheets suffice to account for the scattering. There is no problem of excess turbulent heating because the outer scale for the velocity fluctuations is much larger than the widths of the current sheets. We speculate that interstellar magnetic fields could possess geometries that reflect their origins: fields maintained by the galactic dynamo could have large correlation lengths, whereas those stirred by local energetic events might exhibit folded structures.

We revisit the dynamics of Prometheus and Pandora, two small moons flanking Saturn's F ring. Departures of their orbits from freely precessing ellipses result from mutual interactions via their 121:118 mean motion resonance. Motions are chaotic because the resonance is split into four overlapping components. Orbital longitudes were observed to drift away from Voyager predictions, and a sudden jump in mean motions took place close to the time at which the orbits' apses were antialigned in 2000. Numerical integrations reproduce both the longitude drifts and the jumps. The latter have been attributed to the greater strength of interactions near apse antialignment (every 6.2 years), and it has been assumed that this drift-jump behavior will continue indefinitely. We re-examine the dynamics by analogy with that of a nearly adiabatic, parametric pendulum. In terms of this analogy, the current value of the action of the satellite system is close to its maximum in the chaotic zone. Consequently, at present, the two separatrix crossings per precessional cycle occur close to apse antialignment. In this state libration only occurs when the potential's amplitude is nearly maximal, and the 'jumps' in mean motion arise during the short intervals of libration that separate long stretches of circulation. Because chaotic systems explore the entire region of phase space available to them, we expect that at other times the system would be found in states of medium or low action. In a low action state it would spend most of the time in libration, and separatrix crossings would occur near apse alignment. We predict that transitions between these different states can happen in as little as a decade. Therefore, it is incorrect to assume that sudden changes in the orbits only happen near apse antialignment.

Goertz and Morfill (1983) propose that spokes on Saturn's rings form under radially moving plasma clouds produced by meteoroid impacts. We demonstrate that the speed at which a plasma cloud can move relative to the ring material is bounded from above by the difference between the Keplerian and corotation velocities. The radial orientation of new spokes requires radial speeds that are at least an order of magnitude larger than this upper limit, thus the model advanced by Goertz and Morfill fails to make radial spokes.

We argue that eclipses of radio emission from the millisecond pulsar A in the double pulsar system J0737-3039 are due to synchrotron absorption by plasma in the closed field line region of the magnetosphere of its normal pulsar companion B. A's radio beam only illuminates B's magnetosphere for about 10 minutes surrounding the time of eclipse. During this time it heats particles at r\gtrsim 10^9 cm to relativistic energies and enables extra plasma to be trapped by magnetic mirroring. An enhancement of the plasma density by a factor \sim 10^2 is required to match the duration and optical depth of the observed eclipses. The extra plasma might be supplied by a source near B through B\gamma pair creation by energetic photons produced in B's outer gap. Excitation of pairs' gyrational motions by cyclotron absorption of A's radio beam can result in their becoming trapped between conjugate mirror points in B's magnetosphere. Because the trapping efficiency decreases with increasing optical depth, the plasma density enhancement saturates even under steady state illumination. The result is an eclipse with finite, frequency dependent, optical depth. After illumination by A's radio beam ceases, the trapped particles cool and are lost. The entire cycle repeats every orbital period. We speculate that the asymmetries between eclipse ingress and egress result in part from the magnetosphere's evolution toward a steady state when illuminated by A's radio beam. We predict that A's linear polarization will vary with both eclipse phase and B's rotational phase.

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.

In most books the Delaunay and Lagrange equations for the orbital elements are derived by the Hamilton-Jacobi method: one begins with the 2-body Hamilton equations, performs a canonical transformation to the orbital elements, and obtains the Delaunay system. A standard trick is then used to generalise the approach to the N-body case. We re-examine this step and demonstrate that it contains an implicit condition which restricts the dynamics to a 9(N-1)-dimensional submanifold of the 12(N-1)-dimensional space spanned by the elements and their time derivatives. The tacit condition is equivalent to the so-called Lagrange constraint. It is the condition of the orbital elements being osculating, i.e., of the instantaneous ellipse (or hyperbola) being always tangential to the physical velocity. Imposure of any condition different from the Lagrange constraint (but compatible with the equations of motion) is legitimate and will not alter the physical trajectory or velocity (though will alter the mathematical form of the planetary equations). This freedom of nomination of the supplementary constraint reveals a gauge-type internal symmetry of the celestial-mechanics equations and has consequences for the stability of numerical integrators. Another important aspect of this freedom is that any gauge different from that of Lagrange makes the Delaunay system non-canonical. In a more general setting, when the disturbance depends not only upon positions but also upon velocities, there is a "generalised Lagrange gauge" wherein the Delaunay system is symplectic. This gauge renders orbital elements that are osculating in the phase space. It coincides with the regular Lagrange gauge when the perturbation is velocity-independent.

Cosmic rays scatter off magnetic irregularities (Alfven waves) with which they are resonant, that is waves of wavelength comparable to their gyroradii. These waves may be generated either by the cosmic rays themselves, if they stream faster than the Alfven speed, or by sources of MHD turbulence. Waves excited by streaming cosmic rays are ideally shaped for scattering, whereas the scattering efficiency of MHD turbulence is severely diminished by its anisotropy. We show that MHD turbulence has an indirect effect on cosmic ray propagation by acting as a damping mechanism for cosmic ray generated waves. The hot (``coronal'') phase of the interstellar medium is the best candidate location for cosmic ray confinement by scattering from self-generated waves. We relate the streaming velocity of cosmic rays to the rate of turbulent dissipation in this medium, for the case in which turbulent damping is the dominant damping mechanism. We conclude that cosmic rays with up to 10^2 GeV could not stream much faster than the Alfven speed, but that 10^6 GeV cosmic rays would stream unimpeded by self-generated waves unless the coronal gas were remarkably turbulence-free.

We demonstrate that the chaotic orbits of Prometheus and Pandora are due to interactions associated with the 121:118 mean motion resonance. Differential precession splits this resonance into a quartet of components equally spaced in frequency. Libration widths of the individual components exceed the splitting resulting in resonance overlap which causes the chaos. A single degree of freedom model captures the essential features of the chaotic dynamics. Mean motions of Prometheus and Pandora wander chaotically in zones of width 1.8 deg/yr and 3.1 deg/yr, respectively.

We summarise research reported in (Efroimsky 2002, 2003; Efroimsky and Goldreich 2003a,b) and develop its application to planetary equations in non-inertial frames. We provide a practical example illustrating how the gauge formalism considerably simplifies the calculation of satellite motion about an oblate precessing planet.

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.

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.

Recent HST images of the Saturnian satellites Prometheus and Pandora show that their longitudes deviate from predictions of ephemerides based on Voyager images. Currently Prometheus is lagging and Pandora leading these predictions by somewhat more than 20 degrees. We show that these discrepancies are fully accounted for by gravitational interactions between the two satellites. These peak every 24.8 days at conjunctions and excite chaotic perturbations. The Lyapunov exponent for the Prometheus-Pandora system is of order 0.35 inverse years for satellite masses based on a nominal density of 1.3 gm/cm^3. Interactions are strongest when the orbits come closest together. This happens at intervals of 6.2 years when their apses are anti-aligned. In this context we note the sudden changes of opposite signs in the mean motions of Prometheus and Pandora at the end of 2000 occured shortly after their apsidal lines were anti-aligned.

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.

Two planets with orbital period ratio approximately 10:1 have been discovered around the star HD 83443. The inner and more massive planet, HD 83443b, has the smallest semi-major axis among all currently known exoplanets. Unlike other short period exoplanets, it maintains a substantial orbital eccentricity, e_1=0.079+/- 0.008, in spite of efficient tidal damping. This is a consequence of its secular interactions with HD 83443c whose orbital eccentricity e_2=0.42+/-0.06. Dissipation, associated with tides the star raises in the inner planet, removes energy but not angular momentum from its orbit, while secular interactions transfer angular momentum but not energy from the inner to the outer planet's orbit. The outward transfer of angular momentum decreases the tidal decay rate of the inner planet's orbital eccentricity while increasing that of the outer planet. The alignment of the apsides of the planets' orbits is another consequence of tidal and secular interactions. In this state the ratio of their orbital eccentricities, e_1/e_2, depends upon the secular perturbations the planets exert on each other and on additional perturbations that enhance the inner planet's precession rate. Tidal and rotational distortions of the inner planet along with general relativity provide the most important of these extra precessional perturbations, each of which acts to reduce e_1/e_2. Provided the planets' orbits are coplanar, the observed eccentricity ratio uniquely relates sin i and C=(k_2/k_{2J})(R_1/R_J)^5, where the tidal Love number, k_2, and radius, R_1, of the inner planet are scaled by their Jovian equivalents.

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.

We simulate incompressible MHD turbulence in the presence of a strong background magnetic field. Our major conclusions are: 1) MHD turbulence is most conveniently described in terms of counter propagating shear Alfven and slow waves. Shear Alfven waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfven waves, as does a passive scalar. 2) MHD turbulence is anisotropic with energy cascading more rapidly along k_perp than along k_parallel, where k_perp and k_parallel refer to wavevector components perpendicular and parallel to the local magnetic field. Anisotropy increases with increasing k_perp. 3) MHD turbulence is generically strong in the sense that the waves which comprise it suffer order unity distortions on timescales comparable to their periods. Nevertheless, turbulent fluctuations are small deep inside the inertial range compared to the background field. 4) Decaying MHD turbulence is unstable to an increase of the imbalance between the flux of waves propagating in opposite directions along the magnetic field. 5) Items 1-4 lend support to the model of strong MHD turbulence by Goldreich & Sridhar (GS). Results from our simulations are also consistent with the GS prediction gamma=2/3. The sole notable discrepancy is that 1D power law spectra, E(k_perp) ~ k_perp^{-alpha}, determined from our simulations exhibit alpha ~ 3/2, whereas the GS model predicts alpha = 5/3.

The boundaries of the Uranian epsilon, alpha, and beta rings can be fitted by Keplerian ellipses. The pair of ellipses that outline a given ring share a common line of apsides. Apse alignment is surprising because the quadrupole moment of Uranus induces differential precession. We propose that rigid precession is maintained by a balance of forces due to ring self-gravity, planetary oblateness, and interparticle collisions. Collisional impulses play an especially dramatic role near ring edges. Pressure-induced accelerations are maximal near edges because there (1) velocity dispersions are enhanced by resonant satellite perturbations, and (2) the surface density declines steeply. Remarkably, collisional forces felt by material in the last 100 m of a 10 km wide ring can increase equilibrium masses up to a factor of 100. New ring surface densities are derived which accord with Voyager radio measurements. In contrast to previous models, collisionally modified self-gravity appears to allow for both negative and positive eccentricity gradients; why all narrow planetary rings exhibit positive eccentricity gradients remains an open question.

ZZ Ceti stars exhibit small amplitude photometric pulsations in multiple gravity-modes. We demonstrate that parametric instability, a form of resonant 3-mode coupling, limits overstable modes to amplitudes similar to those observed. In particular, it reproduces the observed trend that longer period modes have larger amplitudes. Parametric instability involves the destabilization of a pair of stable daughter modes by an overstable parent mode. The 3-modes must satisfy exact angular selection rules and approximate frequency resonance. The lowest instability threshold for each parent mode is provided by the daughter pair that minimizes $(\delta\omega^2+\gamma_d^2)/\kappa^2$, where $\kappa$ is the nonlinear coupling constant, $\delta\omega$ is the frequency mismatch, and $\gamma_d$ is the energy damping rate of the daughter modes. The overstable mode's amplitude is maintained at close to the instability threshold value. Although parametric instability defines an upper envelope for the amplitudes of overstable modes in ZZ Ceti stars, other nonlinear mechanisms are required to account for the irregular distribution of amplitudes of similar modes and the non-detection of modes with periods longer than $1,200\s$. Resonant 3-mode interactions involving more than one excited mode may account for the former. Our leading candidate for the latter is Kelvin-Helmholtz instability of the mode-driven shear layer below the convection zone.

When a self-gravitating spherical gas cloud collapses or accretes onto a central mass, the inner region of the cloud develops a density profile $\rho\propto r^{-3/2}$ and the velocity approaches free-fall. We show that in this region, nonspherical perturbations grow with decreasing radius. In the linear regime, the tangential velocity perturbation increases as $r^{-1}$, while the Lagrangian density perturbation, $\Delta\rho/\rho$, grows as $r^{-1/2}$. Faster growth occurs if the central collapsed object maintains a finite multiple moment, in which case $\Delta\rho/\rho$ increases as $r^{-l}$, where $l$ specifies the angular degree of the perturbation. These scaling relations are different from those obtained for the collapse of a homogeneous cloud. Our numerical calculations indicate that nonspherical perturbations are damped in the subsonic region, and that they grow and approach the asymptotic scalings in the supersonic region. The implications of our results to asymmetric supernova collapse and to black hole accretion are briefly discussed.

We compute spectral energy distributions (SEDs) for passive T Tauri disks viewed at arbitrary inclinations. Semi-analytic models of disks in radiative and hydrostatic equilibrium are employed. Over viewing angles for which the flared disk does not occult the central star, the SED varies negligibly with inclination. For such aspects, the SED shortward of ~80 microns is particularly insensitive to orientation, since short wavelength disk emission is dominated by superheated surface layers which are optically thin. The SED of a nearly edge-on disk is that of a class I source. The outer disk occults inner disk regions, and emission shortward of ~30 microns is dramatically extinguished. Spectral features from dust grains may appear in absorption. However, millimeter wavelength fluxes decrease by at most a factor of 2 from face-on to edge-on orientations. We present illustrative applications of our SED models. The class I source 04108+2803B is considered a T Tauri star hidden from view by an inclined circumstellar disk. Fits to its observed SED yield model-dependent values for the disk mass of ~0.015 solar masses and a disk inclination of ~65 degrees relative to face-on. The class II source GM Aur represents a T Tauri star unobscured by its circumstellar disk. Fitted parameters include a disk mass of \~0.050 solar masses and an inclination of ~60 degrees.

We report on numerical calculations of nonadiabatic eigenvalues and eigenfunctions for g-modes in ZZ Ceti variables. The spectrum of overstable $l=1$ modes delineates the instability strip. Its blue edge occurs where $\omega \tau_c \approx 1$ for the $n=1$ mode. Here $\omega$ is radian frequency and $\tau_c$ is about four times the thermal timescale at the bottom of the surface convection zone. As a ZZ Ceti cools, its convection zone deepens, longer period modes become overstable, but the critical value of $\omega\tau_c$ separating overstable and damped modes rises. The latter is a consequence of enhanced radiative damping for modes which propagate immediately below the convection zone. The critical value of $\omega\tau_c$ is of observational significance because modes with the smallest value of $\omega\tau_c$ are most observable photometrically. Maximum periods for overstable modes predicted for our cooler model envelopes are about a factor two longer than the observational upper limit of $1,200\s$. We assess a number of plausible resolutions for this discrepancy among which convective overshoot and nonlinear saturation look promising. The nonadiabatic eigenfunctions enable us to predict relative amplitudes and phases of photospheric variations of flux and velocity, quantities made accessible by recent observations. We also present asymptotic formula for damping rates of high order modes, a result of consequence for future investigations of nonlinear saturation of the amplidues of overstable modes.

We investigate dynamical interactions between turbulent convection and g-mode pulsations in ZZ Ceti variables (DAVs). Since our understanding of turbulence is rudimentary, we are compelled to settle for order of magnitude results. A key feature of these interactions is that convective response times are much shorter than pulsation periods. Thus the dynamical interactions enforce near uniform horizontal velocity inside the convection zone. They also give rise to a narrow shear layer in the region of convective overshoot at the top of the radiative interior. Turbulent damping inside the convection zone is negligible for all modes, but that in the region of convective overshoot may be significant for a few long period modes near the red edge of the instability strip. These conclusions are in accord with those reached earlier by Brickhill. Our major new result concerns nonlinear damping arising from the Kelvin-Helmholtz instability of the aforementioned shear layer. Amplitudes of overstable modes saturate where dissipation due to this instability balances excitation by convective driving. This mechanism of amplitude saturation is most effective for long period modes, and it may play an important role in defining the red edge of the instability strip.

We analyze the stability of g-modes in variable white dwarfs with hydrogen envelopes. In these stars, the radiative layer contributes to mode damping because its opacity decreases upon compression and the amplitude of the Lagrangian pressure perturbation increases outward. The overlying convective envelope is the seat of mode excitation because it acts as an insulating blanket with respect to the perturbed flux that enters it from below. A crucial point is that the convective motions respond to the instantaneous pulsational state. Driving exceeds damping by as much as a factor of two provided $\omega\tau_c\geq 1$, where $\omega$ is the radian frequency of the mode and $\tau_c\approx 4\tau_{th}$ with $\tau_{th}$ being the thermal time constant evaluated at the base of the convective envelope. As a white dwarf cools, its convection zone deepens, and modes of lower frequency become overstable. However, the deeper convection zone impedes the passage of flux perturbations from the base of the convection zone to the photosphere. Thus the photometric variation of a mode with constant velocity amplitude decreases. These factors account for the observed trend that longer period modes are found in cooler DAVs. The linear growth time, ranging from hours for the longest period observed modes ($P\approx 20$ minutes) to thousands of years for those of shortest period ($P\approx 2 $ minutes), probably sets the time-scale for variations of mode amplitude and phase. This is consistent with observations showing that longer period modes are more variable than shorter period ones. Our investigation confirms many results obtained by Brickhill in his pioneering studies of ZZ Cetis.

We derive hydrostatic, radiative equilibrium models for passive disks surrounding T Tauri stars. Each disk is encased by an optically thin layer of superheated dust grains. This layer re-emits directly to space about half the stellar energy it absorbs. The other half is emitted inward and regulates the interior temperature of the disk. The heated disk flares. As a consequence, it absorbs more stellar radiation, especially at large radii, than a flat disk would. The portion of the spectral energy distribution contributed by the disk is fairly flat throughout the thermal infrared. At fixed frequency, the contribution from the surface layer exceeds that from the interior by about a factor 3 and is emitted at more than an order of magnitude greater radius. Spectral features from dust grains in the superheated layer appear in emission if the disk is viewed nearly face-on.

Kraichnan (1965) proposed that MHD turbulence occurs as a result of collisions between oppositely directed Alfv\'en wave packets. Recent work has generated some controversy over the nature of non linear couplings between colliding Alfv\'en waves. We find that the resolution to much of the confusion lies in the existence of a new type of turbulence, intermediate turbulence, in which the cascade of energy in the inertial range exhibits properties intermediate between those of weak and strong turbulent cascades. Some properties of intermediate MHD turbulence are: (i) in common with weak turbulent cascades, wave packets belonging to the inertial range are long lived; (ii) however, components of the strain tensor are so large that, similar to the situation in strong turbulence, perturbation theory is not applicable; (iii) the breakdown of perturbation theory results from the divergence of neighboring field lines due to wave packets whose perturbations in velocity and magnetic fields are localized, but whose perturbations in displacement are not; (iv) 3--wave interactions dominate individual collisions between wave packets, but interactions of all orders $n\geq 3$ make comparable contributions to the intermediate turbulent energy cascade; (v) successive collisions are correlated since wave packets are distorted as they follow diverging field lines; (vi) in common with the weak MHD cascade, there is no parallel cascade of energy, and the cascade to small perpendicular scales strengthens as it reaches higher wave numbers; (vii) For an appropriate weak excitation, there is a natural progression from a weak, through an intermediate, to a strong cascade.