We study the orbital stability of a non-zero mass, close-in circular orbit planet around an eccentric orbit binary for various initial values of the binary eccentricity, binary mass fraction, planet mass, planet semi--major axis, and planet inclination by means of numerical simulations that cover $5 \times 10^4$ binary orbits. For small binary eccentricity, the stable orbits that extend closest to the binary (most stable orbits) are nearly retrograde and circulating. For high binary eccentricity, the most stable orbits are highly inclined and librate near the so-called generalised polar orbit which is a stationary orbit that is fixed in the frame of the binary orbit. For more extreme mass ratio binaries, there is a greater variation in the size of the stability region (defined by initial orbital radius and inclination) with planet mass and initial inclination, especially for low binary eccentricity. For low binary eccentricity, inclined planet orbits may be unstable even at large orbital radii (separation $> 5 \,a_{\rm b}$). The escape time for an unstable planet is generally shorter around an equal mass binary compared with an unequal mass binary. Our results have implications for circumbinary planet formation and evolution and will be helpful for understanding future circumbinary planet observations.

We investigate the formation mechanism for the observed nearly polar aligned (perpendicular to the binary orbital plane) debris ring around the eccentric orbit binary 99 Herculis. An initially inclined nonpolar debris ring or disc will not remain flat and will not evolve to a polar configuration, due to the effects of differential nodal precession that alter its flat structure. However, a gas disc with embedded well coupled solids around the eccentric binary may evolve to a polar configuration as a result of pressure forces that maintain the disc flatness and as a result of viscous dissipation that allows the disc to increase its tilt. Once the gas disc disperses, the debris disc is in a polar aligned state in which there is little precession. We use three-dimensional hydrodynamical simulations, linear theory, and particle dynamics to study the evolution of a misaligned circumbinary gas disc and explore the effects of the initial disc tilt, mass, and size. We find that for a wide range of parameter space, the polar alignment timescale is shorter than the lifetime of the gas disc. Using the observed level of alignment of 3 deg. from polar, we place an upper limit on the mass of the gas disc of about 0.014 M_sun at the time of gas dispersal. We conclude that the polar debris disc around 99 Her can be explained as the result of an initially moderately inclined gas disc with embedded solids. Such a disc may provide an environment for the formation of polar planets.

We investigate the evolution of a multi--planet--disc system orbiting one component of a binary star system. The planet--disc system is initially coplanar but misaligned to the binary orbital plane. The planets are assumed to be giants that open gaps in the disc. We first study the role of the disc in shaping the mutual evolution of the two planets using a secular model for low initial tilt. In general we find that the planets and the disc do not remain coplanar, in agreement with Lubow & Martin (2016) for the single planet case. Instead, the planets and the disc undergo tilt oscillations. A high mass disc between the two planets causes the planets and the disc to nodally precess at the same average rate but they are generally misaligned. The amplitude of the tilt oscillations between the planets is larger while the disc is present. We then consider higher initial tilts using hydrodynamical simulations and explore the possibility of the formation of eccentric Kozai-Lidov (KL) planets. We find that the inner planet's orbit undergoes eccentricity growth for a large range of disc masses and initial misalignments. For a low disc mass and large initial misalignment both planets and the disc can undergo KL oscillations. Furthermore, we find that sufficiently massive discs can cause the inner planet to increase its inclination beyond $90^{\circ}$ and therefore to orbit the binary in a retrograde fashion. The results have important implications for the explanation of very eccentric planets and retrograde planets observed in multi-planet systems.

The Hubble Space Telescope (HST) has obtained multi-epoch observations providing the opportunity for a comprehensive variability search aiming to uncover new variables. We have therefore undertaken the task of creating a catalog of variable sources based on version 3 of the Hubble Source Catalog (HSC), which relies on publicly available images obtained with the WFPC2, ACS, and WFC3 instruments onboard the HST. We adopted magnitude-dependent thresholding in median absolute deviation (a robust measure of light curve scatter) combined with sophisticated preprocessing techniques and visual quality control to identify and validate variable sources observed by Hubble with the same instrument and filter combination five or more times. The Hubble Catalog of Variables (HCV) includes 84,428 candidate variable sources (out of 3.7 million HSC sources that were searched for variability) with $V \leq 27$ mag; for 11,115 of them the variability is detected in more than one filter. The data points in the light curves of the variables in the HCV catalog range from five to 120 points (typically having less than ten points); the time baseline ranges from under a day to over 15 years; while $\sim$8% of all variables have amplitudes in excess of 1 mag. Visual inspection performed on a subset of the candidate variables suggests that at least 80% of the candidate variables that passed our automated quality control are true variable sources rather than spurious detections resulting from blending, residual cosmic rays, and calibration errors. The HCV is the first, homogeneous catalog of variable sources created from archival HST data and currently is the deepest catalog of variables available. The catalog includes variable stars in our Galaxy and nearby galaxies, as well as transients and variable active galactic nuclei. (abbreviated)

Martin & Lubow (2017) found that an initially sufficiently misaligned low mass protoplanetary disc around an eccentric binary undergoes damped nodal oscillations of tilt angle and longitude of ascending node. Dissipation causes evolution towards a stationary state of polar alignment in which the disc lies perpendicular to the binary orbital plane with angular momentum aligned to the eccentricity vector of the binary. We use hydrodynamic simulations and analytic methods to investigate how the mass of the disc affects this process. The simulations suggest that a disc with nonzero mass settles into a stationary state in the frame of the binary, the generalised polar state, at somewhat lower levels of misalignment with respect to the binary orbital plane, in agreement with the analytic model. Provided that discs settle into this generalised polar state, the observational determination of the misalignment angle and binary properties can be used to determine the mass of a circumbinary disc. We apply this constraint to the circumbinary disc in HD 98800. We obtain analytic criteria for polar alignment of a circumbinary ring with mass that approximately agree with the simulation results. Very broad misaligned discs undergo breaking, but the inner regions at least may still evolve to a polar state. The long term evolution of the disc depends on the evolution of the binary eccentricity that we find tends to decrease. Although the range of parameters required for polar alignment decreases somewhat with increasing disc mass, such alignment appears possible for a broad set of initial conditions expected in protostellar circumbinary discs.

Astrometric calibration of images with a small field of view is often inferior to the internal accuracy of the source detections due to the small number of accessible guide stars. One important experiment with such challenges is the Hubble Space Telescope (HST). A possible solution is to cross-calibrate overlapping fields instead of just relying on standard stars. Following the approach of \citet{2012ApJ...761..188B}, we use infinitesimal 3D rotations for fine-tuning the calibration but devise a better objective that is robust to a large number of false candidates in the initial set of associations. Using Bayesian statistics, we accommodate bad data by explicitly modeling the quality, which yields a formalism essentially identical to an $M$-estimation in robust statistics. Our results on simulated and real catalogs show great potentials for improving the HST calibration, and those with similar challenges.

We investigate the dynamics of a nonzero mass, circular orbit planet around an eccentric orbit binary for various values of the binary eccentricity, binary mass fraction, planet mass, and planet semi--major axis by means of numerical simulations. Previous studies investigated the secular dynamics mainly by approximate analytic methods. In the stationary inclination state, the planet and binary precess together with no change in relative tilt. For both prograde and retrograde planetary orbits, we explore the conditions for planetary orbital libration versus circulation and the conditions for stationary inclination. As was predicted by analytic models, for sufficiently high initial inclination, a prograde planet's orbit librates about the stationary tilted state. For a fixed binary eccentricity, the stationary angle is a monotonically decreasing function of the ratio of the planet--to--binary angular momentum $j$. The larger $j$, the stronger the evolutionary changes in the binary eccentricity and inclination. We also calculate the critical tilt angle that separates the circulating from the librating orbits for both prograde and retrograde planet orbits. The properties of the librating orbits and stationary angles are quite different for prograde versus retrograde orbits. The results of the numerical simulations are in very good quantitative agreement with the analytic models. Our results have implications for circumbinary planet formation and evolution.

We investigate the misalignment of the circumbinary disk around the binary HD 98800 BaBb with eccentricity $e\simeq 0.8$. \cite{Kennedy2019} observed the disk to be either at an inclination of $48^{\circ}$ or polar aligned to the binary orbital plane. Their simulations showed that alignment from $48^\circ$ to a polar configuration can take place on a shorter timescale than the age of this system. We perform hydrodynamical numerical simulations in order to estimate the cavity size carved by the eccentric binary for different disk inclinations as an independent check of polar alignment. Resonance theory suggests that torques on the inner parts of a polar disk around such a highly eccentric binary are much weaker than in the coplanar case, indicating a significantly smaller central cavity than in the coplanar case. We show that the inferred inner radius (from carbon monoxide measurements) of the accretion disk around BaBb can exclude the possibility of it being mildly inclined with respect to the binary orbital plane and therefore confirm the polar configuration. This study constitutes an important diagnostic for misaligned circumbinary disks, since it potentially allows us to infer the disk inclination from observed gas disk inner radii.

We analyse the evolution of a mildly inclined circumbinary disc that orbits an eccentric orbit binary by means of smoother particle hydrodynamic (SPH) simulations and linear theory. We show that the alignment process of an initially misaligned circumbinary disc around an eccentric orbit binary is significantly different than around a circular orbit binary and involves tilt oscillations. The more eccentric the binary, the larger the tilt oscillations and the longer it takes to damp these oscillations. A circumbinary disc that is only mildly inclined may increase its inclination by a factor of a few before it moves towards alignment. The results of the SPH simulations agree well with those of linear theory. We investigate the properties of the circumbinary disc/ring around KH 15D. We determine disc properties based on the observational constraints imposed by the changing binary brightness. We find that the inclination is currently at a local minimum and will increase substantially before setting to coplanarity. In addition, the nodal precession is currently near its most rapid rate. The recent observations that show a reappearance of Star B impose constraints on the thickness of the layer of obscuring material. Our results suggest that disc solids have undergone substantial inward drift and settling towards to disc midplane. For disc masses $\sim 0.001 M_\odot$, our model indicates that the level of disc turbulence is low $\alpha \ll 0.001$. Another possibility is that the disc/ring contains little gas.

The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(\alpha < H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R<\alpha \ll 1$ it evolves diffusively. Here $\alpha$ is the viscosity parameter and $H/R$ the disk aspect ratio. Currently there is no simple set of equations which describes the evolution in both regimes. In this paper we describe a possible solution to this problem and introduce a set of one-dimensional equations that describe the evolution of a warped disk that are applicable in both high and low viscosity regimes for arbitrary tilts, but small warps.

We investigate the formation and evolution of misaligned accretion discs around the secondary component of a binary through mass transfer driven by Kozai-Lidov oscillations of the circumprimary disc's eccentricity and inclination. We perform SPH simulations to study the amount of mass transferred to the secondary star as a function of both the disc and binary parameters. For the range of parameters we explore, we find that increasing the disc aspect ratio, viscosity parameter and initial inclination as well as decreasing the binary mass ratio leads to larger amount of mass transfer, up to a maximum of about ten per cent of the initial mass of the primary disc. The circumsecondary disc forms with a high eccentricity and a high inclination and is also able to undergo KL oscillations. The circumsecondary disc oscillations have a shorter period than those in the disc around the primary. We find that some of the material that escapes the Roche-lobe of the two components forms a narrow misaligned circumbinary accretion disc. This study has implications for disc evolution in young binary star systems.

In the past few decades, observations have revealed signatures of metals polluting the atmospheres of white dwarfs. The diffusion timescale for metals to sink from the atmosphere of a white dwarf is of the order of days for a hydrogen-dominated atmosphere. Thus, there must be a continuous supply of metal-rich material accreting onto these white dwarfs. We investigate the role of secular resonances that excite the eccentricity of asteroids allowing them to reach star-grazing orbits leading them to tidal disruption and the formation of a debris disc. Changes in the planetary system during the evolution of the star lead to a change in the location of secular resonances. In our Solar System, the engulfment of the Earth will cause the $\nu_6$ resonance to shift outwards which will force previously stable asteroids to undergo secular resonant perturbations. With analytic models and $N$--body simulations we show that secular resonances driven by two outer companions can provide a source of continuous pollution. Secular resonances are a viable mechanism for the pollution of white dwarfs in a variety of exoplanetary system architectures.

In a recent paper Martin & Lubow showed that a circumbinary disc around an eccentric binary can undergo damped nodal oscillations that lead to the polar (perpendicular) alignment of the disc relative to the binary orbit. The disc angular momentum vector aligns to the eccentricity vector of the binary. We explore the robustness of this mechanism for a low mass disc (0.001 of the binary mass) and its dependence on system parameters by means of hydrodynamic disc simulations. We describe how the evolution depends upon the disc viscosity, temperature, size, binary mass ratio, orbital eccentricity and inclination. We compare results with predictions of linear theory. We show that polar alignment of a low mass disc may occur over a wide range of binary-disc parameters. We discuss the application of our results to the formation of planetary systems around eccentric binary stars.

The Hubble Source Catalog (HSC) combines lists of sources detected on images obtained with the WFPC2, ACS and WFC3 instruments aboard the Hubble Space Telescope (HST) available in the Hubble Legacy Archive. The catalog contains time-domain information with about two million of its sources detected with the same instrument and filter in at least five HST visits. The Hubble Catalog of Variables (HCV) project aims to identify HSC sources showing significant brightness variations. A magnitude-dependent threshold in the median absolute deviation of photometric measurements (an outlier-resistant measure of lightcurve scatter) is adopted as the variability-detection statistic. It is supplemented with a cut in $\chi_{\rm red}^2$ that removes sources with large photometric errors. A pre-processing procedure involving bad image identification, outlier rejection and computation of local magnitude zero-point corrections is applied to HSC lightcurves before computing the variability detection statistic. About 52000 HSC sources are identified as candidate variables, among which 7800 show variability in more than one filter. Visual inspection suggests that $\sim 70\%$ of the candidates detected in multiple filters are true variables while the remaining $\sim 30\%$ are sources with aperture photometry corrupted by blending, imaging artifacts or image processing anomalies. The candidate variables have AB magnitudes in the range 15-27$^{m}$ with the median 22$^{m}$. Among them are the stars in our own and nearby galaxies as well as active galactic nuclei.

The Hubble Catalog of Variables (HCV) project aims to identify the variable sources in the Hubble Source Catalog (HSC), which includes about 92 million objects with over 300 million measurements detected by the WFPC2, ACS and WFC3 cameras on board of the Hubble Space Telescope (HST), by using an automated pipeline containing a set of detection and validation algorithms. All the HSC sources with more than a predefined number of measurements in a single filter/instrument combination are pre-processed to correct systematic effect and to remove the bad measurements. The corrected data are used to compute a number of variability indexes to determine the variability status of each source. The final variable source catalog will contain variables stars, active galactic nuclei (AGNs), supernovae (SNs) or even new types of variables, reaching an unprecedented depth (V$\leq$27 mag). At the end of the project, the first release of the HCV will be available at the Mikulski Archive for Space Telescopes (MAST) and the ESA Hubble Science Archives. The HCV pipeline will be deployed at the Space Telescope Science Institute (STScI) so that an updated HCV may be generated following future releases of HSC.

Martin & Lubow (2017) showed through simulations that an initially tilted disc around an eccentric binary can evolve to polar alignment in which the disc lies perpendicular to the binary orbital plane. We apply linear theory to show both analytically and numerically that a nearly polar aligned low mass circumbinary disc evolves to polar alignment and determine the alignment timescale. Significant disc evolution towards the polar state around moderately eccentric binaries can occur for typical protostellar disc parameters in less than a typical disc lifetime for binaries with orbital periods of order 100 years or less. Resonant torques are much less effective at truncating the inner parts of circumbinary polar discs than the inner parts of coplanar discs. For polar discs, they vanish for a binary eccentricity of unity. The results agree with the simulations in showing that discs can evolve to a polar state. Circumbinary planets may then form in such discs and reside on polar orbits.

We aim to construct an exceptionally deep (V ~< 27) catalog of variable objects in selected Galactic and extragalactic fields visited multiple times by the Hubble Space Telescope (HST). While HST observations of some of these fields were searched for specific types of variables before (most notably, the extragalactic Cepheids), we attempt a systematic study of the population of variable objects of all types at the magnitude range not easily accessible with ground-based telescopes. The variability timescales that can be probed range from hours to years depending on how often a particular field has been visited. For source extraction and cross-matching of sources between visits we rely on the Hubble Source Catalog which includes 10^7 objects detected with WFPC2, ACS, and WFC3 HST instruments. The lightcurves extracted from the HSC are corrected for systematic effects by applying local zero-point corrections and are screened for bad measurements. For each lightcurve we compute variability indices sensitive to a broad range of variability types. The indices characterize the overall lightcurve scatter and smoothness. Candidate variables are selected as having variability index values significantly higher than expected for objects of similar brightness in the given set of observations. The Hubble Catalog of Variables will be released in 2018.

Recent results by Martin et al. (2014) showed in 3D SPH simulations that tilted discs in binary systems can be unstable to the development of global, damped Kozai-Lidov (KL) oscillations in which the discs exchange tilt for eccentricity. We investigate the linear stability of KL modes for tilted inviscid discs under the approximations that the disc eccentricity is small and the disc remains flat. By using 1D equations, we are able to probe regimes of large ratios of outer to inner disc edge radii that are realistic for binary systems of hundreds of AU separations and are not easily probed by multi-dimensional simulations. For order unity binary mass ratios, KL instability is possible for a window of disc aspect ratios H/r in the outer parts of a disc that roughly scale as (n_b/n)^2 < H/r < n_b/n, for binary orbital frequency n_b and orbital frequency n at the disc outer edge. We present a framework for understanding the zones of instability based on the determination of branches of marginally unstable modes. In general, multiple growing eccentric KL modes can be present in a disc. Coplanar apsidal-nodal precession resonances delineate instability branches. We determine the range of tilt angles for unstable modes as a function of disc aspect ratio. Unlike the KL instability for free particles that involves a critical (minimum) tilt angle, disc instability is possible for any nonzero tilt angle depending on the disc aspect ratio.

We use three-dimensional hydrodynamical simulations to show that an initially mildly misaligned circumbinary accretion disk around an eccentric binary can evolve to an orientation that is perpendicular to the orbital plane of the binary (polar alignment). As the disk evolves to the perpendicular state, it undergoes nodal libration oscillations of the tilt angle and the longitude of the ascending node. Dissipation within the disk causes the oscillations to damp. The process operates above a critical initial misalignment angle that depends upon the eccentricity of the binary and the mass of the disk. For binary eccentricity of 0.5, the process operates typically for disk masses smaller than a few percent of the binary mass and initial tilt angle of more than 40 degrees. This evolution has important implications for planet formation around eccentric binary star systems.

We analyze the gravitational instability (GI) of a locally isothermal inclined disk around one component of a binary system. Such a disk can undergo global Kozai-Lidov (KL) cycles if the initial disk tilt is above the critical KL angle (of about 40 degrees). During these cycles, an initially circular disk exchanges its inclination for eccentricity, and vice versa. Self-gravity may suppress the cycles under some circumstances. However, with hydrodynamic simulations including self-gravity we show that for a sufficiently high initial disk tilts and for certain disk masses, disks can undergo KL oscillations and fragment due to GI, even when the Toomre Q value for an equivalent undisturbed disk is well within the stable regime (Q > 2). We suggest that KL triggered disk fragmentation provides a mechanism for the efficient formation of giant planets in binary systems and may enhance fragmentation of disks in massive black hole binaries.

We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow (2011). We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low (2009), the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.

With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central star, observed for some exoplanets. The planet is assumed to be massive enough to open a gap in the disc. The planet's tilt relative to the binary orbital plane is subject to two types of oscillations. The first type, present at even small inclination angles relative to the binary orbital plane, is due to the interaction of the planet with the disc and binary companion and is amplified by a secular resonance. The second type of oscillation is the KL oscillation that operates on both the planet and disc at larger binary inclination angles. We find that for a sufficiently massive disc, even a relatively low inclination planet-disc system can force a planet to an inclination above the critical KL angle, as a consequence of the first type of tilt oscillation, allowing it to undergo the second type of oscillation. We conclude that the hydrodynamical evolution of a sufficiently massive and inclined disc in a binary system broadens the range of systems that form eccentric and misaligned giant planets to include a wide range of initial misalignment angles (20 to 160 degrees).

The Hubble Source Catalog is designed to help optimize science from the Hubble Space Telescope by combining the tens of thousands of visit-based source lists in the Hubble Legacy Archive into a single master catalog. Version 1 of the Hubble Source Catalog includes WFPC2, ACS/WFC, WFC3/UVIS, and WFC3/IR photometric data generated using SExtractor software to produce the individual source lists. The catalog includes roughly 80 million detections of 30 million objects involving 112 different detector/filter combinations, and about 160 thousand HST exposures. Source lists from Data Release 8 of the Hubble Legacy Archive are matched using an algorithm developed by Budavari & Lubow (2012). The mean photometric accuracy for the catalog as a whole is better than 0.10 mag, with relative accuracy as good as 0.02 mag in certain circumstances (e.g., bright isolated stars). The relative astrometric residuals are typically within 10 mas, with a value for the mode (i.e., most common value) of 2.3 mas. The absolute astrometric accuracy is better than $\sim$0.1 arcsec for most sources, but can be much larger for a fraction of fields that could not be matched to the PanSTARRS, SDSS, or 2MASS reference systems. In this paper we describe the database design with emphasis on those aspects that enable the users to fully exploit the catalog while avoiding common misunderstandings and potential pitfalls. We provide usage examples to illustrate some of the science capabilities and data quality characteristics, and briefly discuss plans for future improvements to the Hubble Source Catalog.

We determine the evolution of a giant planet-disk system that orbits a member of a binary star system and is mildly inclined with respect to the binary orbital plane. The planet orbit and disk are initially mutually coplanar. We analyze the evolution of the planet and the disk by analytic means and hydrodynamic simulations. We generally find that the planet and the disk do not remain coplanar unless the disk mass is very large or the gap that separates the planet from the disk is very small. The relative planet-disk tilt undergoes secular oscillations whose initial amplitudes are typically of order the initial disk tilt relative to the binary orbital plane for disk masses ~1% of the binary mass or less. The effects of a secular resonance and the disk tilt decay enhance the planet-disk misalignment. The secular resonance plays an important role for disk masses greater than the planet mass. At later times, the accretion of disk gas by the planet causes its orbit to evolve towards alignment, if the disk mass is sufficiently large. The results have several implications for the evolution of massive planets in binary systems.

Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai-Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the effects of self-gravity, which are expected to suppress the KL oscillations for sufficiently massive disks, were ignored. Here, we analyze the effects of disk self-gravity by means of hydrodynamic simulations and compare the results with the expectations of analytic theory. The disk mass required for suppression in the simulations is a few percent of the mass of the central star and this roughly agrees with an analytical estimate. The conditions for suppression of the KL oscillations in the simulations are close to requiring that the disk be gravitationally unstable. We discuss some implications of our results for the dynamics of protoplanetary disks and the related planet formation.

Martin et al. (2014b) showed that a substantially misaligned accretion disk around one component of a binary system can undergo global damped Kozai-Lidov oscillations. During these oscillations, the inclination and eccentricity of the disk are periodically exchanged. However, the robustness of this mechanism and its dependence on the system parameters were unexplored. In this paper, we use three-dimensional hydrodynamical simulations to analyze how various binary and disk parameters affect the Kozai-Lidov mechanism in hydrodynamical disks. The simulations include the effect of gas pressure and viscosity, but ignore the effects of disk self-gravity. We describe results for different numerical resolutions, binary mass ratios and orbital eccentricities, initial disk sizes, initial disk surface density profiles, disk sound speeds, and disk viscosities. We show that the Kozai-Lidov mechanism can operate for a wide range of binary-disk parameters. We discuss the applications of our results to astrophysical disks in various accreting systems.

We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case.

We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.

We extend previous studies of the tidal truncation of coplanar disks in binary systems to the more general case of noncoplanar disks. As in the prograde coplanar case, Lindblad resonances play a key role in tidal truncation. We analyze the tidal torque acting on a misaligned nearly circular disk in a circular orbit binary system. We concentrate on the 2:1 inner Lindblad resonance associated with the m=2 tidal forcing (for azimuthal wavenumber m) that plays a major role in the usual coplanar case. We determine the inclination dependence of this torque, which is approximately cos^8(i/2) for misalignment angle i. Compared to the prograde coplanar case (i=0), this torque decreases by a factor of about 2 for i = pi/6 and by a factor of about 20 for i=pi/2. The Lindblad torque decreases to zero for a tilt angle of pi (counter-rotation), consistent with previous investigations. The effects of higher order resonances associated with m>2 tidal forcing may contribute somewhat, but are much more limited than in the i=0 case. These results suggest that misaligned disks in binary systems can be significantly extended compared to their coplanar counterparts. In cases where a disk is sufficiently inclined and viscous, it can overrun all Lindblad resonances and overflow the Roche lobe of the disk central object.

We carried out two-dimensional high-resolution simulations to study the effect of dust feedback on the evolution of vortices induced by massive planets in protoplanetary disks. Various initial dust to gas disk surface density ratios ($0.001$ -- $0.01$) and dust particle sizes (Stokes number $4\times10^{-4}$ -- $0.16$) are considered. We found that while dust particles migrate inwards, vortices are very effective in collecting them. When dust density becomes comparable to gas density within the vortex, a dynamical instability is excited and it alters the coherent vorticity pattern and destroys the vortex. This dust feedback effect is stronger with higher initial dust/gas density ratio and larger dust grain. Consequently, we found that the disk vortex lifetime can be reduced up to a factor of 10. We discuss the implications of our findings on the survivability of vortices in protoplanetary disks and planet formation.

We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged periodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.

Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.

Recent observations of large-scale asymmetric features in protoplanetary disks suggest that large-scale vortices exist in such disks. Massive planets are known to be able to produce deep gaps in protoplanetary disks. The gap edges could become hydrodynamically unstable to the Rossby wave/vortex instability and form large-scale vortices. In this study we examine the long term evolution of these vortices by carrying out high-resolution two dimensional hydrodynamic simulations that last more than $10^4$ orbits (measured at the planet's orbit). We find that the disk viscosity has a strong influence on both the emergence and lifetime of vortices. In the outer disk region where asymmetric features are observed, our simulation results suggest that the disk viscous $\alpha$ needs to be low $\sim 10^{-5 }$ - $10^{-4}$ to sustain vortices to thousands and up to $10^{4}$ orbits in certain cases. The chance of finding a vortex feature in a disk then decreases with smaller planet orbital radius. For $\alpha \sim 10^{-3}$ or larger, even planets with masses of 5 Jupiter-masses will have difficulty either producing or sustaining vortices. We have also studied the effects of different disk temperatures and planet masses. We discuss the implications of our findings on current and future protoplanetary disk observations.

Zhu, Stone, and Rafikov (2012) found in 3D shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with results of Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k_y > 1/h (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.

We consider the evolution of accretion discs that contain some turbulence within a disc dead zone, a region about the disc midplane of a disc that is not sufficiently ionised for the magneto-rotational instability (MRI) to drive turbulence. In particular, we determine whether additional sources of turbulence within a dead zone are capable of suppressing gravo-magneto (GM) disc outbursts that arise from a rapid transition from gravitationally produced to MRI produced turbulence. With viscous $\alpha$ disc models we consider two mechanisms that may drive turbulence within the dead zone. First, we examine a constant $\alpha$ parameter within the dead zone. This may be applicable to hydrodynamical instability, such as baroclinic instability, where the turbulence level is independent of the MRI active surface layer properties. In this case, we find that the disc will not become stable to the outbursts unless the dead zone turbulent viscosity is comparable to that in the MRI active surface layers. Under such conditions, the accretion rate through the dead zone must be larger than that through the MRI active layers. In a second model, we assume that the accretion flow though the dead zone is a constant fraction (less than unity) of that through the active layers. This may be applicable to turbulence driven by hydrodynamic waves that are excited by the MRI active surface layers. We find that the instability is hardly affected by the viscous dead zone. In both cases however, we find that the triggering of the MRI during the outburst may be due to the heating from the viscosity in the dead zone, rather than self-gravity. Thus, neither mechanism for generating turbulence within the dead zone can significantly stabilise a disc or the resulting outburst behaviour.

Discs that contain dead zones are subject to the Gravo-Magneto (GM) instability that arises when the turbulence shifts from gravitational to magnetic. We have previously described this instability through a local analysis at some radius in the disc in terms of a limit cycle. A disc may be locally unstable over a radial interval. In this paper, we consider how the local instability model can describe global disc outbursts. The outburst is triggered near the middle of the range of locally unstable radii. The sudden increase in turbulence within high surface density material causes a snow plough of density that propagates both inwards and outwards. All radii inside of the trigger radius become unstable, as well as locally unstable radii outside the trigger radius. In addition, a locally stable region outside of the trigger radius may also become unstable as the gravitational instability is enhanced by the snow plough. For the circumstellar disc model we consider, we find that a quarter of the disc mass is accreted on to the central object during the outburst. The radius out to which the disc is globally unstable is twice that for which it is locally unstable.

Standard, planar accretion discs operate through a dissipative mechanism, usually thought to be turbulent, and often modelled as a viscosity. This acts to take energy from the radial shear, enabling the flow of mass and angular momentum in the radial direction. In a previous paper we discussed observational evidence for the magnitude of this viscosity, and pointed out discrepancies between these values and those obtained in numerical simulations. In this paper we discuss the observational evidence for the magnitude of the dissipative effects which act in non--planar discs, both to transfer and to eliminate the non--planarity. Estimates based on the model by Ogilvie (1999), which assumes a small--scale, isotropic viscosity, give alignment timescales for fully ionized discs which are apparently too short by a factor of a few compared with observations, although we emphasize that more detailed computations as well as tighter observational constraints are required to verify this conclusion. For discs with low temperature and conductivity, we find that the timescales for disc alignment based on isotropic viscosity are too short by around two orders of magnitude. This large discrepancy suggests that our understanding of viscosity in quiescent discs is currently inadequate.

Motivated by recent observations of short-timescale variations in the infrared emission of circumstellar disks, we propose that coronal mass ejections can remove dust grains on timescales as short as a few days. Continuous monitoring of stellar activity, coupled to infrared observations, can place meaningful constraints on the proposed mechanism.

Object cross-identification in multiple observations is often complicated by the uncertainties in their astrometric calibration. Due to the lack of standard reference objects, an image with a small field of view can have significantly larger errors in its absolute positioning than the relative precision of the detected sources within. We present a new general solution for the relative astrometry that quickly refines the World Coordinate System of overlapping fields. The efficiency is obtained through the use of infinitesimal 3-D rotations on the celestial sphere, which do not involve trigonometric functions. They also enable an analytic solution to an important step in making the astrometric corrections. In cases with many overlapping images, the correct identification of detections that match together across different images is difficult to determine. We describe a new greedy Bayesian approach for selecting the best object matches across a large number of overlapping images. The methods are developed and demonstrated on the Hubble Legacy Archive, one of the most challenging data sets today. We describe a novel catalog compiled from many Hubble Space Telescope observations, where the detections are combined into a searchable collection of matches that link the individual detections. The matches provide descriptions of astronomical objects involving multiple wavelengths and epochs. High relative positional accuracy of objects is achieved across the Hubble images, often sub-pixel precision in the order of just a few milli-arcseconds. The result is a reliable set of high-quality associations that are publicly available online.

Regular satellites in the solar system are thought to form within circumplanetary discs. We consider a model of a layered circumplanetary disc that consists of a nonturbulent midplane layer and and strongly turbulent disc surface layers. The dead zone provides a favorable site for satellite formation. It is a quiescent environment that permits the growth of solid bodies. Viscous torques within the disc cause it to expand to a substantial fraction of its Hill radius (~0.4 R_H) where tidal torques from the central star remove its angular momentum. For certain parameters, the dead zone develops into a high density substructure well inside the Hill sphere. The radial extent of the dead zone may explain the compactness of the regular satellites orbits for Jupiter and Saturn. The disc temperatures can be low enough to be consistent with the high ice fractions of Ganymede and Callisto.

We present high spatial resolution maps of ro-vibrational molecular hydrogen emission from the environment of the GG Tau A binary component in the GG Tau quadruple system. The H2 v= 1-0 S(1) emission is spatially resolved and encompasses the inner binary, with emission detected at locations that should be dynamically cleared on several hundred-year timescales. Extensions of H2 gas emission are seen to ~100 AU distances from the central stars. The v = 2-1 S(1) emission at 2.24 microns is also detected at ~30 AU from the central stars, with a line ratio of 0.05 +/- 0.01 with respect to the v = 1-0 S(1) emission. Assuming gas in LTE, this ratio corresponds to an emission environment at ~1700 K. We estimate that this temperature is too high for quiescent gas heated by X-ray or UV emission from the central stars. Surprisingly, we find that the brightest region of H2 emission arises from a spatial location that is exactly coincident with a recently revealed dust "streamer" which seems to be transferring material from the outer circumbinary ring around GG Tau A into the inner region. As a result, we identify a new excitation mechanism for ro-vibrational H2 stimulation in the environment of young stars. The H2 in the GG Tau A system appears to be stimulated by mass accretion infall as material in the circumbinary ring accretes onto the system to replenish the inner circumstellar disks. We postulate that H2 stimulated by accretion infall could be present in other systems, particularly binaries and "transition disk" systems which have dust cleared gaps in their circumstellar environments.

We perform global time-dependent simulations of an accretion disc around a young stellar object with a dead zone (a region where the magneto-rotational instability cannot drive turbulence because the material is not sufficiently ionised). For infall accretion rates on to the disc of around 10^-7 Msun/yr, dead zones occur if the critical magnetic Reynolds number is larger than about 10^4. We model the collapse of a molecular gas cloud. At early times when the infall accretion rate is high, the disc is thermally ionised and fully turbulent. However, as the infall accretion rate drops, a dead zone may form if the critical magnetic Reynolds number is sufficiently large, otherwise the disc remains fully turbulent. With a dead zone the disc can become unstable to the gravo-magneto instability. The mass of the star grows in large accretion outbursts that may explain FU Orionis events. At late times there is not sufficient mass in the disc for outbursts to occur but the dead zone becomes even more prominent as the disc cools. Large inner dead zones in the later stages of disc evolution may help to explain observations of transition discs with an inner hole.

We describe a model for the long term evolution of a circumplanetary disk that is fed mass from a circumstellar disk and contains regions of low turbulence (dead zones). We show that such disks can be subject to accretion driven outbursts, analogous to outbursts previously modeled in the context of circumstellar disks to explain FU Ori phenomena. Circumplanetary disks around a proto-Jupiter can undergo outbursts for infall accretion rates onto the disks in the range ~10^{-9} to 10^{-7} M_sun/yr, typical of accretion rates in the T Tauri phase. During outbursts, the accretion rate and disk luminosity increases by several orders of magnitude. Most of the planet mass growth during planetary gas accretion may occur via disk outbursts involving gas that is considerably hotter than predicted by steady state models. For low infall accretion rates less than ~10^{-10} M_sun/yr that occur in late stages of disk accretion, disk outbursts are unlikely to occur, even if dead zones are present. Such conditions are favorable for the formation of icy satellites.

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations, strong torques by the binary can maintain a gap of radius $\simeq 2a$. Streams curve inward from $r \simeq 2a$ toward the binary; some of their mass passes through the inner boundary, while the remainder swings back out to the disk. However, we also find that near its inner edge the disk develops both a strong $m=1$ asymmetry and growing orbital eccentricity. Because the MHD stresses introduce more matter into the gap, the total torque per unit disk mass is $\simeq 14$ times larger than found previously. The inner boundary accretion rate per unit disk mass is $\simeq 40$ times greater than found from previous hydrodynamical calculations. The implied binary shrinkage rate is determined by a balance between the rate at which the binary gains angular momentum by accretion and loses it by gravitational torque. The large accretion rate brings these two rates nearly into balance, but in net, we find that $\dot a/a < 0$, and its magnitude is about 2.7 times larger than predicted by the earlier hydrodynamic simulations. If the binary comprises two massive black holes, the accretion rate may be great enough for one or both to be AGN, with consequences for the physical state of the gas both in the disk body and in its inner gap.

Angular momentum is transported outwards through an accretion disc by magnetohydrodynamical (MHD) turbulence thus allowing material to accrete on to the central object. The magneto-rotational instability (MRI) requires a minimum ionisation fraction to drive turbulence in a disc. The inner parts of the disc around a young stellar object are sufficiently hot to be thermally ionised. Further out, cosmic rays ionise the surface layers and a dead zone forms at the mid-plane where the disc is too cool for the MRI to operate. The surface density in the turbulent active layer is often assumed to be constant with radius because the cosmic rays penetrate a constant layer. However, if a critical magnetic Reynolds number, Re_{M,crit}, is used to determine the extent of the dead zone, the surface density in the layer generally increases with radius. For small critical magnetic Reynolds number of order 1, the constant layer approximation may be a reasonable fit. However, MHD simulations suggest the critical magnetic Reynolds number may be much larger, of order 10^4. Analytical fits for the surface density in the magnetic active layer show that \Sigma_m \propto Re_{M,crit}^{-2} R^{9/2} T^{\,2}$, at temperature T and radius R, are a good fit for higher critical magnetic Reynolds number. For the metallicity variation between our galaxy, the LMC and the SMC, there should be no significant difference in the extent of the dead zone. Observations suggest an increase in the lifetime of the disc with decreasing metallicity that cannot be explained by the dead zone structure (ignoring possible differences in dust abundances).

Previous theoretical studies have found that repeating outbursts can occur in certain regions of an accretion disk, due to sudden transitions in time from gravitationally produced turbulence to magnetically produced turbulence. We analyze the disk evolution in a state diagram that plots the mass accretion rate versus disk surface density. We determine steady state accretion branches that involve gravitational and magnetic sources of turbulence. Using time-dependent numerical disk simulations, we show that cases having outbursts track along a nonsteady 'dead zone' branch and some steady state accretion branches. The outburst is the result of a rapid inter-branch transition. The gravo-magneto outbursts are then explained on this diagram as a limit cycle that is analogous to the well-known S-curve that has been applied to dwarf nova outbursts. The diagram and limit cycle provide a conceptual framework for understanding the nature of the outbursts that may occur in accretion disks of all scales, from circumplanetary to protoplanetary to AGN accretion disks.

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed edge on. Transitions between the two implies that the discs may be warped and precessing. Type II X-ray outbursts are thought to occur when the warped disc interacts with the fast stellar wind. We suggest that tides from a misaligned companion neutron star can cause the observed effects. We make numerical models of a Be star decretion disc in which the spin of the Be star is misaligned with the orbital axis of a neutron star companion. Tidal torques from the neutron star truncate the disc at a radius small enough that the neutron star orbit does not intersect the disc unless the eccentricity or misalignment is very large. A magnetic torque from the Be star that is largest at the equator, where the rotation is fastest, is approximated by an inner boundary condition. There are large oscillations in the mass and inclination of the disc as it moves towards a steady state. These large variations may explain the observed changes from Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find the tidal timescale on which the disc warps, precesses and reaches a steady state to be around a year up to a few hundred years. If present, the oscillations in mass and disc inclination occur on a fraction of this timescale depending on the orbital parameters of the binary. The timescales associated with the tidal torque for observed Be star binaries suggest that these effects are important in all but the longest period binaries.

We analyse some properties of circumplanetary discs. Flow through such discs may provide most of the mass to gas giant planets, and such discs are likely sites for the formation of regular satellites. We model these discs as accretion discs subject to the tidal forces of the central star. The tidal torques from the star remove the disc angular momentum near the disc outer edge and permit the accreting disc gas to lose angular momentum at the rate appropriate for steady accretion. Circumplanetary discs are truncated near the radius where periodic ballistic orbits cross, where tidal forces on the disc are strong. This radius occurs at approximately 0.4 r_H for the planet Hill radius r_H. During the T Tauri stage of disc accretion, the disc is fairly thick with aspect ratio H/r > 0.2 and the disc edge tapering occurs over a radial scale ~ H ~ 0.1 r_H. For a circular or slightly eccentric orbit planet, no significant resonances lie within the main body of the disc. Tidally driven waves involving resonances nonetheless play an important role in truncating the disc, especially when it is fairly thick. We model the disc structure using one dimensional time-dependent and steady-state models and also two dimensional SPH simulations. The circumplanetary disc structure depends on the variation of the disc turbulent viscosity with radius and is insensitive to the angular distribution of the accreting gas. Dead zones may occur within the circumplanetary disc and result in density structures. If the disc is turbulent throughout, the predicted disc structure near the location of the regular Jovian and Saturnian satellites is smooth with no obvious feature that would favor formation at their current locations.

We determine an expression for the Type I planet migration torque involving a locally isothermal disk, with moderate turbulent viscosity (~0.0005 < alpha < ~0.05), based on three-dimensional nonlinear hydrodynamical simulations. The radial gradients (in a dimensionless logarithmic form) of density and temperature are assumed to be constant near the planet. We find that the torque is roughly equally sensitive to the surface density and temperature radial gradients. Both gradients contribute to inward migration when they are negative. Our results indicate that two-dimensional calculations with a smoothed planet potential, used to account for the effects of the third dimension, do not accurately determine the effects of density and temperature gradients on the three-dimensional torque. The results suggest that substantially slowing or stopping planet migration by means of changes in disk opacity or shadowing is difficult and appears unlikely for a disk that is locally isothermal. The scalings of the torque and torque density with planet mass and gas sound speed follow the expectations of linear theory. We also determine an improved formula for the torque density distribution that can be used in one-dimensional long-term evolution studies of planets embedded in locally isothermal disks. This formula can be also applied in the presence of mildly varying radial gradients and of planets that open gaps. We illustrate its use in the case of migrating super-Earths and determine some conditions sufficient for survival.

Planet migration is the process by which a planet's orbital radius changes in time. The main agent for causing gas giant planet migration is the gravitational interaction of the young planet with the gaseous disk from which it forms. We describe the migration rates resulting from these interactions based on a simple model for disk properties. These migration rates are higher than is reasonable for planet survival. We discuss some proposed models for which the migration rates are lower. There are major uncertainties in migration rates due to a lack of knowledge about the detailed physical properties of disks. We also describe some additional forms of migration.

We consider discs that orbit a central object and are tidally perturbed by a circular orbit companion. Such discs are sometimes subject to an eccentric instability due to the effects of certain resonances. Eccentric instabilities may be present in planetary rings perturbed by satellites, protostellar discs perturbed by planets, and discs in binary star systems. Although the basic mechanism for eccentric instability is well understood, the detailed response of a gaseous disc to such an instability is not understood. We apply a linear eccentricity evolution equation developed by Goodchild and Ogilvie. We explore how the eccentricity is distributed in such a disc and how the distribution in turn affects the instability growth rate for a range of disc properties. We identify a disc mode, termed the superhump mode, that is likely at work in the superhump binary star case. The mode results from the excitation of the fundamental free precession mode. We determine an analytic expression for the fundamental free mode precession rate that is applicable to a sufficiently cool disc. Depending on the disc sound speed and disc edge location, other eccentric modes can grow faster than the superhump mode and dominate.

We carry out 2-D high resolution numerical simulations of type I planet migration with different disk viscosities. We find that the planet migration is strongly dependent on disk viscosities. Two kinds of density wave damping mechanisms are discussed. Accordingly, the angular momentum transport can be either viscosity dominated or shock dominated, depending on the disk viscosities. The long term migration behavior is different as well. Influences of the Rossby vortex instability on planet migration are also discussed. In addition, we investigate very weak shock generation in inviscid disks by small mass planets and compare the results with prior analytic results.

We describe 2D hydrodynamic simulations of the migration of low-mass planets ($\leq 30 M_{\oplus}$) in nearly laminar disks (viscosity parameter $\alpha < 10^{-3}$) over timescales of several thousand orbit periods. We consider disk masses of 1, 2, and 5 times the minimum mass solar nebula, disk thickness parameters of $H/r = 0.035$ and 0.05, and a variety of $\alpha$ values and planet masses. Disk self-gravity is fully included. Previous analytic work has suggested that Type I planet migration can be halted in disks of sufficiently low turbulent viscosity, for $\alpha \sim 10^{-4}$. The halting is due to a feedback effect of breaking density waves that results in a slight mass redistribution and consequently an increased outward torque contribution. The simulations confirm the existence of a critical mass ($M_{cr} \sim 10 M_{\oplus}$) beyond which migration halts in nearly laminar disks. For $\alpha \ga 10^{-3}$, density feedback effects are washed out and Type I migration persists. The critical masses are in good agreement with the analytic model of Rafikov (2002). In addition, for $\alpha \la 10^{-4}$ steep density gradients produce a vortex instability, resulting in a small time-varying eccentricity in the planet's orbit and a slight outward migration. Migration in nearly laminar disks may be sufficiently slow to reconcile the timescales of migration theory with those of giant planet formation in the core accretion model.

We analyze the orbital and mass evolution of planets that undergo run-away gas accretion by means of 2D and 3D hydrodynamic simulations. The disk torque distribution per unit disk mass as a function of radius provides an important diagnostic for the nature of the disk-planet interactions. We first consider torque distributions for nonmigrating planets of fixed mass and show that there is general agreement with the expectations of resonance theory. We then present results of simulations for mass-gaining, migrating planets. For planets with an initial mass of 5 Earth masses, which are embedded in disks with standard parameters and which undergo run-away gas accretion to one Jupiter mass (Mjup), the torque distributions per unit disk mass are largely unaffected by migration and accretion for a given planet mass. The migration rates for these planets are in agreement with the predictions of the standard theory for planet migration (Type I and Type II migration). The planet mass growth occurs through gas capture within the planet's Bondi radius at lower planet masses, the Hill radius at intermediate planet masses, and through reduced accretion at higher planet masses due to gap formation. During run-away mass growth, a planet migrates inwards by only about 20% in radius before achieving a mass of ~1 Mjup. For the above models, we find no evidence of fast migration driven by coorbital torques, known as Type III migration. We do find evidence of Type III migration for a fixed mass planet of Saturn's mass that is immersed in a cold and massive disk. In this case the planet migration is assumed to begin before gap formation completes. The migration is understood through a model in which the torque is due to an asymmetry in density between trapped gas on the leading side of the planet and ambient gas on the trailing side of the planet.

The exoplanets discovered so far have been mostly around relatively nearby and bright stars. As a result, the host stars are mostly (i) in the Galactic disk, (ii) relatively massive, and (iii) relatively metal rich. The aim of the SWEEPS project is to extend our knowledge to stars which (i) are in a different part of the Galaxy, (ii) have lower masses, and (iii) have a large range of metallicities. To achieve this goal, we used the Hubble Space Telescope to search for transiting planets around F, G, K, and M dwarfs in the Galactic bulge. We photometrically monitored 180,000 stars in a dense bulge field continuously for 7 days. We discovered 16 candidate transiting extrasolar planets with periods of 0.6 to 4.2 days, including a new class of ultra-short period planets (USPPs) with P < 1.2 days. Radial-velocity observations of the two brightest candidates support their planetary nature. These results suggest that planets are as abundant in the Galactic bulge as they are in the solar neighborhood, and they are equally abundant around low-mass stars (within a factor 2). The planet frequency increases with metallicity even for the stars in the Galactic bulge. All the USPP hosts are low-mass stars, suggesting either that close-in planets around higher-mass stars are irradiatively evaporated, or that the planets can migrate to close-in orbits only around such old and low-mass stars.

There is evidence for the existence of massive planets at orbital radii of several hundred AU from their parent stars where the timescale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.

We present the results of a coronagraphic scattered-light imaging survey of six young disk candidate stars using the Hubble Space Telescope Advanced Camera for Surveys. The observations made use of the 1.8" occulting spot through the F606W (broad V) filter. Circumstellar material was imaged around HD 97048, a Herbig Ae/Be star located in the Chamaeleon I dark cloud at a distance of 180 pc. The material is seen between ~2" (360 AU) and ~4" (720 AU) from the star in all directions. A V-band azimuthally-averaged radial surface brightness profile peaks at r = 2" with a value of 19.6 +/- 0.2 mag arcsec^-2 and smoothly decreases with projected distance from the star as r^(-3.3 +/- 0.5). An integrated flux of 16.8 +/- 0.1 mag is measured between 2" and 4", corresponding to a scattered-light fractional luminosity lower limit of 8.4 x 10^-4. Filamentary structure resembling spiral arms similar to that seen in Herbig Ae/Be disks is observed. Such structure has been attributed to the influence of orbiting planets or stellar encounters. Average surface brightness upper limits are determined for the five non-detections: HD 34282, HD 139450, HD 158643, HD 159492, and HD 195627. Possible reasons for the non-detections are disks that are too faint or disks hidden by the occulter.

More than 200 extrasolar planets have been discovered around relatively nearby stars, primarily through the Doppler line shifts owing to the reflex motions of their host stars, and more recently through transits of some planets across the face of the host stars. The detection of planets with the shortest known periods, 1.2 to 2.5 days, has mainly resulted from transit surveys which have generally targeted stars more massive than 0.75 M_sun. Here we report the results from a planetary transit search performed in a rich stellar field towards the Galactic bulge. We discovered 16 candidates with orbital periods between 0.4 and 4.2 days, five of which orbit stars of 0.44 to 0.75 M_sun. In two cases, radial-velocity measurements support the planetary nature of the companions. Five candidates have orbital periods below 1.0 day, constituting a new class of ultra-short-period planets (USPPs), which occur only around stars of less than 0.88 M_sun. This indicates that those orbiting very close to more luminous stars might be evaporatively destroyed, or that jovian planets around lower-mass stars might migrate to smaller radii.

We investigate the interaction between a giant planet and a viscous circumstellar disk by means of high-resolution, two-dimensional hydrodynamical simulations. We consider planet masses that range from 1 to 3 Jupiter masses (Mjup) and initial orbital eccentricities that range from 0 to 0.4. We find that a planet can cause eccentricity growth in a disk region adjacent to the planet's orbit, even if the planet's orbit is circular. Disk-planet interactions lead to growth in a planet's orbital eccentricity. The orbital eccentricities of a 2 Mjup and a 3 Mjup planet increase from 0 to 0.11 within about 3000 orbits. Over a similar time period, the orbital eccentricity of a 1 Mjup planet grows from 0 to 0.02. For a case of a 1 Mjup planet with an initial eccentricity of 0.01, the orbital eccentricity grows to 0.09 over 4000 orbits. Radial migration is directed inwards, but slows considerably as a planet's orbit becomes eccentric. If a planet's orbital eccentricity becomes sufficiently large, e > ~0.2, migration can reverse and so be directed outwards. The accretion rate towards a planet depends on both the disk and the planet orbital eccentricity and is pulsed over the orbital period. Planet mass growth rates increase with planet orbital eccentricity. For e~0.2 the mass growth rate of a planet increases by approximately 30% above the value for e=0. For e > ~0.1, most of the accretion within the planet's Roche lobe occurs when the planet is near the apocenter. Similar accretion modulation occurs for flow at the inner disk boundary which represents accretion toward the star.

The migration of a planet through a gaseous disc causes the locations of their resonant interactions to drift and can alter the torques exerted between the planet and the disc. We analyse the time-dependent dynamics of a non-coorbital corotation resonance under these circumstances. The ratio of the resonant torque in a steady state to the value given by Goldreich & Tremaine (1979) depends essentially on two dimensionless quantities: a dimensionless turbulent diffusion time-scale and a dimensionless radial drift speed. When the drift speed is comparable to the libration speed and the viscosity is small, the torque can become much larger than the unsaturated value in the absence of migration, but is still proportional to the large-scale vortensity gradient in the disc. Fluid that is trapped in the resonance and drifts with it acquires a vortensity anomaly relative to its surroundings. If the anomaly is limited by viscous diffusion in a steady state, the resulting torque is inversely proportional to the viscosity, although a long time may be required to achieve this state. A further, viscosity-independent, contribution to the torque comes from fluid that streams through the resonant region. In other cases, torque oscillations occur before the steady value is achieved. We discuss the significance of these results for the evolution of eccentricity in protoplanetary systems. We also describe the possible application of these findings to the coorbital region and the concept of runaway (or type III) migration. [Abridged]

We analyze the gas accretion flow through a planet-produced gap in a protoplanetary disk. We adopt the alpha disk model and ignore effects of planetary migration. We develop a semi-analytic, one-dimensional model that accounts for the effects of the planet as a mass sink and also carry out two-dimensional hydrodynamical simulations of a planet embedded in a disk. The predictions of the mass flow rate through the gap based on the semi-analytic model generally agree with the hydrodynamical simulations at the 25% level. Through these models, we are able to explore steady state disk structures and over large spatial ranges. The presence of an accreting Jupiter-mass planet significantly lowers the density of the disk within a region of several times the planet's orbital radius. The mass flow rate across the gap (and onto the central star) is typically 10% to 25% of the mass accretion rate outside the orbit of the planet, for planet-to-star mass ratios that range from 5e-5 to 1e-3.

We consider the alignment torque between a spinning black hole and an accretion disc whose angular momenta are misaligned. This situation must hold initially in almost all gas accretion events on to supermassive black holes, and may occur in binaries where the black hole receives a natal supernova kick. We show that the torque always acts to align the hole's spin with the total angular momentum without changing its magnitude. The torque acts dissipatively on the disc, reducing its angular momentum, and aligning it with the hole if and only if the angle theta between the angular momenta J_d of the disc and J_h of the hole satisfies the inequality cos theta > -J_d / 2 J_h. If this condition fails, which requires both theta > pi/2 and J_d < 2 J_h, the disc counteraligns.

We study the dynamical origin of the structures observed in the scattered-light images of the resolved debris disk around HD 141569A. We explore the roles of radiation pressure from the central star, gas drag from the gas disk, and the tidal forces from two nearby stars in creating and maintaining these structures. We use a simple one-dimensional axisymmetric model to show that the presence of the gas helps confine the dust and that a broad ring of dust is produced if a central hole exists in the disk. This model also suggests that the disk is in a transient, excited dynamical state, as the observed dust creation rate applied over the age of the star is inconsistent with submillimeter mass measurements. We model in two dimensions the effects of a fly-by encounter between the disk and a binary star in a prograde, parabolic, coplanar orbit. We track the spatial distribution of the disk's gas, planetesimals, and dust. We conclude that the surface density distribution reflects the planetesimal distribution for a wide range of parameters. Our most viable model features a disk of initial radius 400 AU, a gas mass of 50 M_earth, and beta = 4 and suggests that the system is being observed within 4000 yr of the fly-by periastron. The model reproduces some features of HD 141569A's disk, such as a broad single ring and large spiral arms, but it does not reproduce the observed multiple spiral rings or disk asymmetries nor the observed clearing in the inner disk. For the latter, we consider the effect of a 5 M_Jup planet in an eccentric orbit on the planetesimal distribution of HD 141569A.

We investigate the migration rates of high-mass protoplanets embedded in accretion discs via two and three-dimensional hydrodynamical simulations. The simulations follow the planet's radial motion and employ a nested-grid code that allows for high resolution close to the planet. We concentrate on the possible role of the coorbital torques in affecting migration rates. We analyse two cases: (a) a Jupiter-mass planet in a low-mass disc and (b) a Saturn-mass planet in a high-mass disc. Planet migration in case (b) is much more susceptible to coorbital torques than in case (a). We find that the coorbital torques in both cases do not depend sensitively on whether the planet is allowed to migrate through the disc or is held on a fixed orbit. We also examine the dependence of the planet's migration rate on the numerical resolution near the planet. For case (a), numerical convergence is relatively easy to obtain, even when including torques arising from deep within the planet's Hill sphere. The migration rate in this case is numerically on order of the Type II migration rate and much smaller than the Type I rate, if the disc has 0.01 solar-masses inside 26 AU. Torques from within the Hill sphere provide a substantial opposing contribution to the migration rate. In case (b), the gas mass within the Hill sphere is larger than the planet's mass and convergence is more difficult to obtain. If the torques within the Hill sphere are ignored, convergence is more easily achieved but the migration rate is artificially large. At our highest resolution, the migration rate for case (b) is much less than the Type I rate, but somewhat larger than the Type II rate.

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from one Earth mass (1 M$_e$) to one Jupiter mass (1 M$_J$) by using the ZEUS hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses $M \gsim 0.1$ M$_J$ produce significant perturbations in the disc's surface density. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc midplane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the planet's orbit. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M$_J$ and is highly efficient, occurring at the local viscous rate. The migration timescales for planets of mass less than 0.1 M$_J$, based on torques from disc material outside the planets' Roche lobes, are in excellent agreement with the linear theory of Type I (non-gap) migration for three-dimensional discs. The transition from Type I to Type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M$_J$ with little migration. Planets with final masses of order 10 M$_J$ would undergo large migration, which makes formation and survival difficult.

We determine the torque exerted in a steady state by an external potential on a three-dimensional gaseous disk at a non-coorbital corotation resonance. Our model accounts for the feedback of the torque on the surface density and vorticity in the corotation region, and assumes that the disk has a barotropic equation of state and a nonzero effective viscosity. The ratio of the torque to the value given by the formula of Goldreich & Tremaine depends essentially on a single dimensionless parameter, which quantifies the extent to which the resonance is saturated. We discuss the implications for the eccentricity evolution of young planets.

We present three-dimensional numerical simulations of the interaction of a circular-orbit planet with a protoplanetary disk. We calculate the flow pattern, the accretion rate, and torques on the planet. We consider planet masses ranging from 1 Earth mass to 1 Jupiter mass.

We consider the evolution of a warped disc around a Kerr black hole, under conditions such that the warp propagates in a wavelike manner. This occurs when the dimensionless effective viscosity, alpha, that damps the warp is less than the characteristic angular semi-thickness, H/R, of the disc. We adopt linearized equations that are valid for warps of sufficiently small amplitude in a Newtonian disc, but also account for the apsidal and nodal precession that occur in the Kerr metric. Through analytical and time-dependent studies, we confirm the results of Demianski & Ivanov, and of Ivanov & Illarionov, that such a disc takes on a characteristic warped shape. The inner part of the disc is not necessarily aligned with the equator of the hole, even in the presence of dissipation. We draw attention to the fact that this might have important implications for the directionality of jets emanating from discs around rotating black holes.

We investigate the migration of massive extrasolar planets due to gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 AU at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log (a), where a is the orbital radius. Between 1 AU and 3 AU, the expected number of planets per logarithmic interval in orbital radius roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, about 15 percent of stars targeted by existing radial velocity searches are predicted to possess planets with masses 0.3 M_Jupiter < M_p sin (i) < 10 M_Jupiter, and radii 0.1 AU < a < 5 AU. A third of these planets (around 5 percent of the target stars) lie at the radii most amenable to detection via microlensing. A further 5-10 percent of stars could have planets at radii of 5 AU < a < 8 AU that have migrated outwards. We discuss the probability of forming a system (akin to the Solar System) in which significant radial migration of the most massive planet does not occur. About 10-15 percent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low mass planets than high mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic ($\gamma = 5/3$). Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance ${\cal M}_{\rm r} \ga 0.01$ undergo shocks within a distance of order the resonance radius.

Doppler spectroscopy of rho1 Cnc has detected evidence of a companion with an orbital period of 14.65 days and a minimum mass of 0.88 Jupiter masses. Astrometric observations performed with the Hubble Space Telescope Fine Guidance Sensor 1r using a novel new observing technique have placed an upper limit on the astrometric reflex motion of rho1 Cnc in a time period of only one month. These observations detected no reflex motion induced by the 14.65 day period radial velocity companion, allowing us to place a 3-sigma upper limit of \~0.3 milliarcseconds on the semi-major axis of this motion, ruling out the preliminary Hipparcos value of 1.15 milliarcseconds. The corresponding upper limit on the true mass of the companion is ~30 Mj, confirming that it is a sub-stellar object.

A planet of low mass orbiting in a two-dimensional gaseous disc generates a one-armed spiral wake. We explain this phenomenon as the result of constructive interference between wave modes in the disc, somewhat similar to the Kelvin wedge produced in the wake of a ship. The same feature is not expected in a three-dimensional disc with thermal stratification.

Withdrawn in current form.

In a planetary system, a secular particle resonance occurs at a location where the precession rate of a test particle (e.g. an asteroid) matches the frequency of one of the precessional modes of the planetary system. We investigate the secular interactions of a system of mutually inclined planets with a gaseous protostellar disk that may contain a secular nodal particle resonance. We determine the normal modes of some mutually inclined planet-disk systems. The planets and disk interact gravitationally, and the disk is internally subject to the effects of gas pressure, self-gravity, and turbulent viscosity. The behavior of the disk at a secular resonance is radically different from that of a particle, owing mainly to the effects of gas pressure. The resonance is typically broadened by gas pressure to the extent that global effects, including large-scale warps, dominate. The standard resonant torque formula is invalid in this regime. Secular interactions cause a decay of the inclination at a rate that depends on the disk properties, including its mass, turbulent viscosity, and sound speed. For a Jupiter-mass planet embedded within a minimum-mass solar nebula having typical parameters, dissipation within the disk is sufficient to stabilize the system against tilt growth caused by mean-motion resonances.

In a recent paper, Elmegreen (2000) has made a cogent case, from an observational point of view, that the lifetimes of molecular clouds are comparable to their dynamical timescales. If so, this has important implications for the mechanisms by which molecular clouds form. In particular we consider the hypothesis that molecular clouds may form not by {\it in situ} cooling of atomic gas, but rather by the agglomeration of the dense phase of the interstellar medium (ISM), much, if not most, of which is already in molecular form.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

We consider the dynamics of a protostellar disk surrounding a star in a circular-orbit binary system. Our aim is to determine whether, if the disk is initially tilted with respect to the plane of the binary orbit, the inclination of the system will increase or decrease with time. The problem is formulated in the binary frame in which the tidal potential of the companion star is static. We consider a steady, flat disk that is aligned with the binary plane and investigate its linear stability with respect to tilting or warping perturbations. The dynamics is controlled by the competing effects of the m=0 and m=2 azimuthal Fourier components of the tidal potential. In the presence of dissipation, the m=0 component causes alignment of the system, while the m=2 component has the opposite tendency. We find that disks that are sufficiently large, in particular those that extend to their tidal truncation radii, are generally stable and will therefore tend to alignment with the binary plane on a time-scale comparable to that found in previous studies. However, the effect of the m=2 component is enhanced in the vicinity of resonances where the outer radius of the disk is such that the natural frequency of a global bending mode of the disk is equal to twice the binary orbital frequency. Under such circumstances, the disk can be unstable to tilting and acquire a warped shape, even in the absence of dissipation. The outer radius corresponding to the primary resonance is always smaller than the tidal truncation radius. For disks smaller than the primary resonance, the m=2 component may be able to cause a very slow growth of inclination through the effect of a near resonance that occurs close to the disk center. We discuss these results in the light of recent observations of protostellar disks in binary systems.

We analyze the nonlinear, two-dimensional response of a gaseous, viscous protoplanetary disk to the presence of a planet of one Jupiter mass (1 M_J) and greater that orbits a 1 solar mass star by using the ZEUS hydrodynamics code with high resolution near the planet's Roche lobe. The planet is assumed to be in a circular orbit about the central star and is not allowed to migrate. A gap is formed about the orbit of the planet, but there is a nonaxisymmetric flow through the gap and onto the planet. The gap partitions the disk into an inner (outer) disk that extends inside (outside) the planet's orbit. For a 1 M_J planet and typical disk parameters, the accretion through the gap onto the planet is highly efficient. For typical disk parameters, the mass doubling time scale is less than 10^5 years, considerably shorter than the disk lifetime. Following shocks near the L1 and L2 Lagrange points, disk material enters the Roche lobe in the form of two gas streams. Shocks occur within the Roche lobe as the gas streams collide, and shocks lead to rapid inflow towards the planet within much of planet's Roche lobe. Shocks also propagate in the inner and outer disks that orbit the star. For higher mass planets (of order 6 M_J), the flow rate onto the planet is considerably reduced, which suggests an upper mass limit to planets in the range of 10 M_J. This rate reduction is related to the fact that the gap width increases relative to the Roche (Hill sphere) radius with increasing planetary mass. The flow in the gap affects planetary migration. For the 1 M_J planet case, mass can penetrate from the outer disk to the inner disk, so that the inner disk is not depleted. The results suggest that most of the mass in gas giant planets is acquired by flows through gaps.

We extend our analysis of the three-dimensional response of a vertically polytropic disk to tidal forcing at Lindblad resonances by including the effects of a disk atmosphere. The atmosphere is modeled as an isothermal layer that joins smoothly on to an underlying polytropic layer. The launched wave progressively enters the atmosphere as it propagates away from the resonance. The wave never propagates vertically, however, and the wave energy rises to a (finite) characteristic height in the atmosphere. The increase of wave amplitude associated with this process of wave channeling is reduced by the effect of the atmosphere. For waves of large azimuthal mode number m generated by giant planets embedded in a disk, the increase in wave amplitude is still substantial enough to be likely to dissipate the wave energy by shocks for even modest optical depths (tau greater than about 10) over a radial distance of a few times the disk thickness. For low-m waves generated in circumstellar disks in binary stars, the effects of wave channeling are less important and the level of wave nonlinearity increases by less than a factor of 10 in going from the disk edge to the disk center. For circumbinary disks, the effects of wave channeling remain important, even for modest values of optical depth.

We analyze the linear, 3D response to tidal forcing of a disk that is thin and thermally stratified in the direction normal to the disk plane. We model the vertical disk structure locally as a polytrope which represents a disk of high optical depth. We solve the 3D gas-dynamic equations semi-analytically in the neighborhood of a Lindblad resonance. These solutions match asymptotically on to those valid away from resonances and provide solutions valid at all radii. We obtain the following results. 1) A variety of waves are launched at resonance. However, the f mode carries more than 95% of the torque exerted at the resonance. 2) These 3D waves collectively transport exactly the amount of angular momentum predicted by the 2D torque formula. 3) Near resonance, the f mode occupies the full vertical extent of the disk. Away from resonance, the f mode becomes confined near the surface of the disk, and, in the absence of other dissipation mechanisms, damps via shocks. The radial length scale for this process is roughly r_L/m (for resonant radius r_L and azimuthal wavenumber m), independent of the disk thickness H. This wave channeling process is due to the variations of physical quantities in r and is not due to wave refraction. 4) However, the inwardly propagating f mode launched from an m=2 inner Lindblad resonance experiences relatively minor channeling. We conclude that for binary stars, tidally generated waves in highly optically thick circumbinary disks are subject to strong nonlinear damping by the channeling mechanism, while those in circumstellar accretion disks are subject to weaker nonlinear effects. We also apply our results to waves excited by young planets for which m is approximately r/H and conclude that the waves are damped on the scale of a few H.

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiter's mass around a solar-like star could have dynamically survived in an orbit with a period as short as $\sim10\,$hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun.