Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$.

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

Motivated by the current search for exomoons, this paper considers the stability of tidal equilibrium for hierarchical three-body systems containing a star, a planet, and a moon. In this treatment, the energy and angular momentum budgets include contributions from the planetary orbit, lunar orbit, stellar spin, planetary spin, and lunar spin. The goal is to determine the optimized energy state of the system subject to the constraint of constant angular momentum. Due to the lack of a closed form solution for the full three-body problem, however, we must use use an approximate description of the orbits. We first consider the Keplerian limit and find that the critical energy states are saddle points, rather than minima, so that these hierarchical systems have no stable tidal equilibrium states. We then generalize the calculation so that the lunar orbit is described by a time-averaged version of the circular restricted three-body problem. In this latter case, the critical energy state is a shallow minimum, so that a tidal equilibrium state exists. In both cases, however, the lunar orbit for the critical point lies outside the boundary (roughly half the Hill radius) where (previous) numerical simulations indicate dynamical instability. These results suggest that star-planet-moon systems have no viable long-term stable states analogous to those found for two-body systems.

Motivated by the possibility that different versions of the laws of physics could be realized within other universes, this paper delineates the galactic parameters that allow for habitable planets and revisits constraints on the amplitude $Q$ of the primordial density fluctuations. Previous work indicates that large values of $Q$ lead to galaxies so dense that planetary orbits cannot survive long enough for life to develop. Small values of $Q$ lead to delayed star formation, loosely bound galaxies, and compromised heavy element retention. This work generalizes previous treatments: [A] We consider models for the internal structure of galaxies and find the fraction of galactic real estate that allows stable, long-lived planetary orbits. [B] We perform a large ensemble of numerical simulations to estimate cross sections for the disruption of planetary orbits due to interactions with passing stars. [C] We consider disruption due to the background radiation fields produced by the galaxies. [D] One consequence of intense galactic background radiation fields is that some portion of the galaxy, denoted as the Galactic Habitable Zone, will provide the right flux levels to support habitable planets for essentially any planetary orbit. As $Q$ increases, the fraction of stars in a galaxy that allow for habitable planets decreases due to both orbital disruption and the intense background radiation. However, the outer parts of the galaxy always allow for habitable planets, so that the value of $Q$ does not have a well-defined upper limit. Moreover, some Galactic Habitable Zones are large enough to support more potentially habitable planets than the galaxies found in our universe. These results suggest that the possibilities for habitability in other universes are somewhat more favorable and far more diverse than previously imagined.

This paper considers the stability of tidal equilibria for planetary systems in which stellar rotation provides a significant contribution to the angular momentum budget. We begin by applying classic stability considerations for two bodies to planetary systems --- where one mass is much smaller than the other. The application of these stability criteria to a subset of the Kepler sample indicates that the majority of the systems are not in a stable equilibrium state. Motivated by this finding, we generalize the stability calculation to include the quadrupole moment for the host star. In general, a stable equilibrium requires that the total system angular momentum exceeds a minimum value (denoted here as $L_X$) and that the orbital angular momentum of the planet exceeds a minimum fraction of the total. Most, but not all, of the observed planetary systems in the sample have enough total angular momentum to allow an equilibrium state. Even with the generalizations of this paper, however, most systems have too little orbital angular momentum (relative to the total) and are not in an equilibrium configuration. Finally, we consider the time evolution of these planetary systems; the results constrain the tidal quality factor of the stars and suggest that $10^6\le{Q_\ast}\le10^7$.

Intermediate mass stars and stellar remnants often host planets, and these dynamical systems evolve because of mass loss and tides. This paper considers the combined action of stellar mass loss and tidal dissipation on planetary orbits in order to determine the conditions required for planetary survival. Stellar mass loss is included using a so-called Jeans model, described by a dimensionless mass loss rate \gamma and an index \beta. We use an analogous prescription to model tidal effects, described here by a dimensionless dissipation rate \Gamma and two indices (q,p). The initial conditions are determined by the starting value of angular momentum parameter \eta (equivalently, the initial eccentricity) and the phase \theta of the orbit. Within the context of this model, we derive an analytic formula for the critical dissipation rate \Gamma, which marks the boundary between orbits that spiral outward due to stellar mass loss and those that spiral inward due to tidal dissipation. This analytic result \Gamma=\Gamma(\gamma,\beta,q,p,\eta,\theta) is essentially exact for initially circular orbits and holds to within an accuracy of 50% over the entire multi-dimensional parameter space, where the individual parameters vary by several orders of magnitude. For stars that experience mass loss, the stellar radius often displays quasi-periodic variations, which produce corresponding variations in tidal forcing; we generalize the calculation to include such pulsations using a semi-analytic treatment that holds to the same accuracy as the non-pulsating case. These results can be used in many applications, e.g., to predict/constrain properties of planetary systems orbiting white dwarfs.

Observations indicate that intermediate mass stars, binary stars, and stellar remnants often host planets; a complete explanation of these systems requires an understanding of how planetary orbits evolve as their central stars lose mass. Motivated by these dynamical systems, this paper generalizes in two directions previous studies of orbital evolution in planetary systems with stellar mass loss: [1] Many previous treatments focus on constant mass loss rates and much of this work is carried out numerically. Here we study a class of single planet systems where the stellar mass loss rate is time dependent. The mass loss rate can be increasing or decreasing, but the stellar mass always decreases monotonically. For this class of models, we develop analytic approximations to specify the final orbital elements for planets that remain bound after the epoch of mass loss, and find the conditions required for the planets to become unbound. We also show that for some mass loss functions, planets become unbound only in the asymptotic limit where the stellar mass vanishes. [2] We consider the chaotic evolution for two planet systems with stellar mass loss. Here we focus on a model consisting of analogs of Jupiter, Saturn, and the Sun. By monitoring the divergence of initially similar trajectories through time, we calculate the Lyapunov exponents of the system. This analog solar system is chaotic in the absence of mass loss with Lyapunov time in the range 5 - 10 Myr; we find that the Lyapunov time decreases with increasing stellar mass loss rate, with a nearly linear relationship between the two time scales. Taken together, the results of this paper help provide an explanation for a wide range of dynamical evolution that occurs in solar systems with stellar mass loss.

Hill's equations arise in a wide variety of physical problems, and are specified by a natural frequency, a periodic forcing function, and a forcing strength parameter. This classic problem is generalized here in two ways: [A] to Random Hill's equations which allow the forcing strength q_k, the oscillation frequency \lambda_k, and the period \tau_k of the forcing function to vary from cycle to cycle, and [B] to Stochastic Hill's equations which contain (at least) one additional term that is a stochastic process \xi. This paper considers both random and stochastic Hill's equations with small parameter variations, so that p_k=q_k-<q_k>, \ell_k=\lambda_k-<\lambda_k>, and \xi are all O(\epsilon), where \epsilon<<1. We show that random Hill's equations and stochastic Hill's equations have the same growth rates when the parameter variations p_k and \ell_k obey certain constraints given in terms of the moments of \xi. For random Hill's equations, the growth rates for the solutions are given by the growth rates of a matrix transformation, under matrix multiplication, where the matrix elements vary from cycle to cycle. Unlike classic Hill's equations where the parameter space (the \lambda-q plane) displays bands of stable solutions interlaced with bands of unstable solutions, random Hill's equations are generically unstable. We find analytic approximations for the growth rates of the instability; for the regime where Hill's equation is classically stable, and the parameter variations are small, the growth rate \gamma = O(\epsilon^2). Using the relationship between the (\ell_k,p_k) and the \xi, this result for \gamma can be used to find growth rates for stochastic Hill's equations.

Motivated by the large number of extrasolar planetary systems that are near mean motion resonances, this paper explores a related type of dynamical behavior known as "nodding". Here, the resonance angle of a planetary system executes libration (oscillatory motion) for several cycles, circulates for one or more cycles, and then enters once again into libration. This type of complicated dynamics can affect our interpretation of observed planetary systems that are in or near mean motion resonance. This work shows that planetary systems in (near) mean motion resonance can exhibit nodding behavior, and outlines the portion of parameter space where it occurs. This problem is addressed using both full numerical integrations of the planetary systems and via model equations obtained through expansions of the disturbing function. In the latter approach, we identify the relevant terms that allow for nodding. The two approaches are in agreement, and show that nodding often occurs when a small body is in an external mean motion resonance with a larger planet. As a result, the nodding phenomenon can be important for interpreting observations of transit timing variations, where the existence of smaller bodies is inferred through their effects on larger, observed transiting planets. For example, in actively nodding planetary systems, both the amplitude and frequency of the transit timing variations depend on the observational time window.

The observed population of Hot Jupiters displays a stunning variety of physical properties, including a wide range of densities and core sizes for a given planetary mass. Motivated by the observational sample, this paper studies the accretion of rocky planets by Hot Jupiters, after the Jovian planets have finished their principal migration epoch and become parked in $\sim4$-day orbits. In this scenario, rocky planets form later and then migrate inward due to torques from the remaining circumstellar disk, which also damps the orbital eccentricity. This mechanism thus represents one possible channel for increasing the core masses and metallicities of Hot Jupiters. This paper determines probabilities for the possible end states for the rocky planet: collisions with the Jovian planets, accretion onto the star, ejection from the system, and long-term survival of both planets. These probabilities depend on the mass of the Jovian planet and its starting orbital eccentricity, as well as the eccentricity damping rate for the rocky planet. Since these systems are highly chaotic, a large ensemble ($N\sim10^3$) of simulations with effectively equivalent starting conditions is required. Planetary collisions are common when the eccentricity damping rate is sufficiently low, but are rare otherwise. For systems that experience planetary collisions, this work determines the distributions of impact velocities -- both speeds and impact parameters -- for the collisions. These velocity distributions help determine the consequences of the impacts, e.g., where energy and heavy elements are deposited within the giant planets.

Pairs of migrating extrasolar planets often lock into mean motion resonance as they drift inward. This paper studies the convergent migration of giant planets (driven by a circumstellar disk) and determines the probability that they are captured into mean motion resonance. The probability that such planets enter resonance depends on the type of resonance, the migration rate, the eccentricity damping rate, and the amplitude of the turbulent fluctuations. This problem is studied both through direct integrations of the full 3-body problem, and via semi-analytic model equations. In general, the probability of resonance decreases with increasing migration rate, and with increasing levels of turbulence, but increases with eccentricity damping. Previous work has shown that the distributions of orbital elements (eccentricity and semimajor axis) for observed extrasolar planets can be reproduced by migration models with multiple planets. However, these results depend on resonance locking, and this study shows that entry into -- and maintenance of -- mean motion resonance depends sensitively on migration rate, eccentricity damping, and turbulence.

This paper derives expressions for the growth rates for the random 2 x 2 matrices that result from solutions to the random Hill's equation. The parameters that appear in Hill's equation include the forcing strength and oscillation frequency. The development of the solutions to this periodic differential equation can be described by a discrete map, where the matrix elements are given by the principal solutions for each cycle. Variations in the forcing strength and oscillation frequency lead to matrix elements that vary from cycle to cycle. This paper presents an analysis of the growth rates including cases where all of the cycles are highly unstable, where some cycles are near the stability border, and where the map would be stable in the absence of fluctuations. For all of these regimes, we provide expressions for the growth rates of the matrices that describe the solutions.

This paper presents a generalized treatment of Type I planetary migration in the presence of stochastic perturbations. In many planet-forming disks, the Type I migration mechanism, driven by asymmetric torques, acts on a short time scale and compromises planet formation. If the disk also supports MHD instabilities, however, the corresponding turbulent fluctuations produce additional stochastic torques that modify the steady inward migration scenario. This work studies the migration of planetary cores in the presence of stochastic fluctuations using complementary methods, including a Fokker-Planck approach and iterative maps. Stochastic torques have two main effects: [1] Through outward diffusion, a small fraction of the planetary cores can survive in the face of Type I inward migration. [2] For a given starting condition, the result of any particular realization of migration is uncertain, so that results must be described in terms of the distributions of outcomes. In addition to exploring different regimes of parameter space, this paper considers the effects of the outer disk boundary condition, varying initial conditions, and time-dependence of the torque parameters. For disks with finite radii, the fraction of surviving planets decreases exponentially with time. We find the survival fractions and decay rates for a range of disk models, and find the expected distribution of locations for surviving planets. For expected disk properties, the survival fraction lies in the range $0.01 < p_S < 0.1$.

This paper considers random Hill's equations in the limit where the periodic forcing function becomes a Dirac delta function. For this class of equations, the forcing strength $q_k$, the oscillation frequency $\af_k$, and the period are allowed to vary from cycle to cycle. Such equations arise in astrophysical orbital problems in extended mass distributions, in the reheating problem for inflationary cosmologies, and in periodic Schr{\"o}dinger equations. The growth rates for solutions to the periodic differential equation can be described by a matrix transformation, where the matrix elements vary from cycle to cycle. Working in the delta function limit, this paper addresses several coupled issues: We find the growth rates for the $2 \times 2$ matrices that describe the solutions. This analysis is carried out in the limiting regimes of both large $q_k \gg 1$ and small $q_k \ll 1$ forcing strength parameters. For the latter case, we present an alternate treatment of the dynamics in terms of a Fokker-Planck equation, which allows for a comparison of the two approaches. Finally, we elucidate the relationship between the fundamental parameters $(\af_k,q_k)$ appearing in the stochastic differential equation and the matrix elements that specify the corresponding discrete map. This work provides analytic -- and accurate -- expressions for the growth rates of these stochastic differential equations in both the $q_k \gg1 $ and the $q_k \ll 1$ limits.

This paper continues previous work on the effects of turbulence on mean motion resonances in extrasolar planetary systems. Turbulence is expected to arise in the disks that form planets, and these fluctuations act to compromise resonant configurations. This paper extends previous work by considering how interactions between the planets and possible damping effects imposed by the disk affect the outcomes. These physical processes are studied using three approaches: numerical integrations of the 3-body problem with additional forcing due to turbulence, model equations that reduce the problem to stochastically driven oscillators, and Fokker-Planck equations that describe the time evolution of an ensemble of systems. With this combined approach, we elucidate the physics of how turbulence can remove extrasolar planetary systems from mean motion resonance. As expected, systems with sufficiently large damping (dissipation) can maintain resonance, in spite of turbulent forcing. In the absence of strong damping, ensembles of these systems exhibit two regimes of behavior, where the fraction of the bound states decreases as a power-law or as an exponential. Both types of behavior can be understood through the model developed herein. For systems with weak interactions between planets, the model reduces to a stochastic pendulum, and the fraction of bound states decreases as a power-law. For highly interactive systems, the dynamics are more complicated and the fraction of bound states decreases exponentially. We show how planetary interactions lead to drift terms in the Fokker-Planck equation and account for this exponential behavior. In addition to clarifying the physical processes involved, this paper strengthens the finding that turbulence implies that mean motions resonances should be rare.

This paper considers the effects of turbulence on mean motion resonances in extrasolar planetary systems and predicts that systems rarely survive in a resonant configuration. A growing number of systems are reported to be in resonance, which is thought to arise from the planet migration process. If planets are brought together and moved inward through torques produced by circumstellar disks, then disk turbulence can act to prevent planets from staying in a resonant configuration. This paper studies this process through numerical simulations and via analytic model equations, where both approaches include stochastic forcing terms due to turbulence. We explore how the amplitude and forcing time intervals of the turbulence affect the maintenance of mean motion resonances. If turbulence is common in circumstellar disks during the epoch of planet migration, with the amplitudes indicated by current MHD simulations, then planetary systems that remain deep in mean motion resonance are predicted to be rare. More specifically, the fraction of resonant systems that survive over a typical disk lifetime of 1 Myr is of order 0.01. If mean motion resonances are found to be common, their existence would place tight constraints on the amplitude and duty cycle of turbulent fluctuations in circumstellar disks. These results can be combined by expressing the expected fraction of surviving resonant systems in the approximate form P_b = C / N_{orb}^{1/2}, where the dimensionless parameter C = 10 - 50 and where N_{orb} is the number of orbits for which turbulence is active.

The rosette-shaped motion of a particle in a central force field is known to be classically solvable by quadratures. We present a new approach of describing and characterizing such motion based on the eccentricity vector of the two body problem. In general, this vector is not an integral of motion. However, the orbital motion, when viewed from the nonuniformly rotating frame defined by the orientation of the eccentricity vector, can be solved analytically and will either be a closed periodic circulation or libration. The motion with respect to inertial space is then given by integrating the argument of periapsis with respect to time. Finally we will apply the decomposition to a modern central potential, the spherical Hernquist-Newton potential, which models dark matter halos of galaxies with central black holes.

Motivated by a class of orbit problems in astrophysics, this paper considers solutions to Hill's equation with forcing strength parameters that vary from cycle to cycle. The results are generalized to include period variations from cycle to cycle. The development of the solutions to the differential equation is governed by a discrete map. For the general case of Hill's equation in the unstable limit, we consider separately the case of purely positive matrix elements and those with mixed signs; we then find exact expressions, bounds, and working estimates for the growth rates. We also find exact expressions, estimates, and bounds for the infinite products of several $2 \times 2$ matrices with random variables in the matrix elements. In the limit of sharply spiked forcing terms (the delta function limit), we find analytic solutions for each cycle and for the discrete map that matches solutions from cycle to cycle; for this case we find the growth rates and the condition for instability in the limit of large forcing strength, as well as the widths of the stable/unstable zones.

This paper constructs an analytic form for a triaxial potential that describes the dynamics of a wide variety of astrophysical systems, including the inner portions of dark matter halos, the central regions of galactic bulges, and young embedded star clusters. Specifically, this potential results from a density profile of the form $\rho (m) \propto m^{-1}$, where the radial coordinate is generalized to triaxial form so that $m^2 = x^2/a^2 + y^2/b^2 + z^2/c^2 $. Using the resulting analytic form of the potential, and the corresponding force laws, we construct orbit solutions and show that a robust orbit instability exists in these systems. For orbits initially confined to any of the three principal planes, the motion in the perpendicular direction can be unstable. We discuss the range of parameter space for which these orbits are unstable, find the growth rates and saturation levels of the instability, and develop a set of analytic model equations that elucidate the essential physics of the instability mechanism. This orbit instability has a large number of astrophysical implications and applications, including understanding the formation of dark matter halos, the structure of galactic bulges, the survival of tidal streams, and the early evolution of embedded star clusters.

This paper constructs an analytic framework for calculating the assembly of galactic disks from the collapse of gas within dark matter halos, with the goal of determining the surface density profiles. Gas parcels (baryons) fall through the potentials of dark matter halos on nearly ballistic, zero energy orbits and collect in a rotating disk. The dark matter halos have a nearly universal form, as determined previously through numerical simulations. The calculation is first carried out for a variety of pre-collapse mass distributions and rotation profiles, including polytropic spheres in hydrostatic equilibrium with the halo potential. The resulting disk surface density profiles have nearly power-law forms, with well-defined edges. This idealized scenario is generalized to include non-spherical starting states and multiple accretion events (due to gas being added to the halo via merger events). This latter complication is explored in detail and considers a log-normal distribution for the angular momenta of the pre-collapse states. If this distribution is wide, then the composite surface density approaches a universal power-law form, independent of the shape of the constituent profiles. When the angular momentum distribution has an intermediate width, the composite surface density attains a nearly exponential form, roughly consistent with profiles of observed galaxies.

This paper explores orbits in extended mass distributions and develops an analytic approximation scheme based on epicycloids (spirograph patterns). We focus on the Hernquist potential which provides a good model for many astrophysical systems, including elliptical galaxies, dark matter halos, and young embedded star clusters. For a given potential, one can readily calculate orbital solutions as a function of energy and angular momentum using numerical methods. In contrast, this paper presents a number of analytic results for the Hernquist potential and proves a series of general constraints showing that orbits have similar properties for any extended mass distribution (including, e.g., the NFW profile). We discuss circular orbits, radial orbits, zero energy orbits, different definitions of eccentricity, analogs of Kepler's law, the definition of orbital elements, and the relation of these orbits to spirograph patterns (epicycloids). Over much of parameter space the orbits can be adequately described (with accuracy better than 10%) using the parametric equations of epicycloids, thereby providing an analytic description of the orbits. As an application of this formal development, we find a solution for the orbit of the Large Magellanic Cloud in the potential of our Galaxy.