Gravitational coupling between young planets and their parent disks is often explored using numerical simulations, which typically treat the disk thermodynamics in a highly simplified manner. In particular, many studies adopt the locally isothermal approximation, in which the disk temperature is a fixed function of the stellocentric distance. We explore the dynamics of planet-driven density waves in disks with more general thermodynamics, in which the temperature is relaxed towards an equilibrium profile on a finite cooling timescale $t_{\rm c}$. We use both linear perturbation theory and direct numerical simulations to examine the global structure of density waves launched by planets in such disks. A key diagnostic used in this study is the behavior of the wave angular momentum flux (AMF), which directly determines the evolution of the underlying disk. The AMF of free waves is constant for slowly cooling (adiabatic) disks, but scales with the disk temperature for rapidly cooling (and locally isothermal) disks. However, cooling must be extremely fast, with $\beta = \Omega t_{\rm c} \lesssim 10^{-3}$ for the locally isothermal approximation to provide a good description of density wave dynamics in the linear regime (relaxing to $\beta \lesssim 10^{-2}$ when nonlinear effects are important). For intermediate cooling timescales, density waves are subject to a strong linear damping. This modifies the appearance of planet-driven spiral arms and the characteristics of axisymmetric structures produced by massive planets: in disks with $\beta \approx 0.1$ -- $1$, a near-thermal mass planet opens only a single wide gap around its orbit, in contrast to the several narrow gaps produced when cooling is either faster or slower.

Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced in collisions into the corresponding mass bins results in O(N^3) calculation, which leads to significant computational overhead when N is large. Here we formulate a more efficient explicit O(N^2) fragmentation algorithm, which works when the size spectrum of fragments produced in an individual collision has a self-similar shape with only a single characteristic mass scale (which can have arbitrary dependence on the energy and masses of colliding objects). Fragment size spectra used in existing fragmentation codes typically possess this property. We also show that our O(N^2) approach can be easily extended to work with non-self-similar fragment size distributions, for which we provide a worked example. This algorithm offers a substantial speedup of fragmentation calculations for large N> 100, even over the implicit methods, making it an attractive tool for studying collisionally evolving systems.

In many astrophysical problems involving discs (gaseous or particulate) orbiting a dominant central mass, gravitational potential of the disc plays an important dynamical role. Its impact on the motion of external objects, as well as on the dynamics of the disc itself, can usually be studied using secular approximation. This is often done using softened gravity to avoid singularities arising in calculation of the orbit-averaged potential --- disturbing function --- of a razor-thin disc using classical Laplace-Lagrange theory. We explore the performance of several softening formalisms proposed in the literature in reproducing the correct eccentricity dynamics in the disc potential. We identify softening models that, in the limit of zero softening, give results converging to the expected behavior exactly, approximately or not converging at all. We also develop a general framework for computing secular disturbing function given an arbitrary softening prescription for a rather general form of the interaction potential. Our results demonstrate that numerical treatments of the secular disc dynamics, representing the disc as a collection of $N$ gravitationally interacting annuli, are rather demanding: for a given value of the (dimensionless) softening parameter, $\varsigma\ll 1$, accurate representation of eccentricity dynamics requires $N \sim C\varsigma^{-\chi}\gg 1$, with $C\sim O(10)$, $1.5\lesssim \chi\lesssim 2$. In discs with sharp edges a very small value of the softening parameter $\varsigma$ ($\lesssim 10^{-3}$) is required to correctly reproduce eccentricity dynamics near the disc boundaries; this finding is relevant for modelling planetary rings.

The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive black hole) until gravitational wave emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH-BH, NS-BH and NS-NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several pc. For NS-NS binaries we find merger rates in the range $0.01-0.07$ Gpc$^{-3}$ yr$^{-1}$ from globular clusters and $0.1-0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters. For NS-BH and BH-BH binaries we find small merger rates from globular clusters, but a rate of $0.1 - 0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several per cent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future gravitational wave detectors.

Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $\Gamma$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $\Gamma$ is varied. We find that for $\Gamma > 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $\Gamma = 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $\gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters.

Orbital evolution of binary systems in dense stellar clusters is important in a variety of contexts: origin of blue stragglers, progenitors of compact object mergers, millisecond pulsars, and so on. Here we consider the general problem of secular evolution of the orbital elements of a binary system driven by the smooth tidal field of an axisymmetric stellar cluster (globular, nuclear, etc.) in which the binary orbits. We derive a secular Hamiltonian (averaged over both the inner Keplerian orbit of the binary and its outer orbit within the cluster) valid to quadrupole order for an arbitrary cluster potential and explore its characteristics. This doubly-averaged 'tidal' Hamiltonian depends on just two parameters, which fully absorb the information about the background cluster potential and the binary's orbit within it: a dimensional parameter $A$ setting the secular timescale, and a dimensionless parameter $\Gamma$ which determines the phase portrait of the binary's inner orbital evolution. We examine the dependence of $A$ and $\Gamma$ on cluster potential (both spherical and axisymmetric) and on the binary orbit within the cluster. Our theory reproduces known secular results - such as Lidov-Kozai evolution and the effect of the Galactic tide on Oort Cloud comets - in appropriate limits, but is more general. It provides a universal framework for understanding dynamical evolution of various types of binaries driven by the smooth tidal field of any axisymmetric potential. In a companion paper (Hamilton & Rafikov 2019b) we provide a detailed exploration of the resulting orbital dynamics.

Massive planetary cores embedded in protoplanetary discs are believed to accrete extended atmospheres, providing a pathway to forming gas giants and gas-rich super-Earths. The properties of these atmospheres strongly depend on the nature of the coupling between the atmosphere and the surrounding disc. We examine the formation of gaseous envelopes around massive planetary cores via three-dimensional inviscid and isothermal hydrodynamic simulations. We focus the changes in the envelope properties as the core mass varies from low (sub-thermal) to high (super-thermal) values, a regime relevant to close-in super-Earths. We show that global envelope properties such as the amount of rotational support or turbulent mixing are mostly sensitive to the ratio of the Bondi radius of the core to its physical size. High-mass cores are fed by supersonic inflows arriving along the polar axis and shocking on the densest parts of the envelope, driving turbulence and mass accretion. Gas flows out of the core's Hill sphere in the equatorial plane, describing a global mass circulation through the envelope. The shell of shocked gas atop the core surface delimits regions of slow (inside) and fast (outside) material recycling by gas from the surrounding disc. While recycling hinders the runaway growth towards gas giants, the inner regions of protoplanetary atmospheres, more immune to mixing, may remain bound to the planet.

Measurements of exoplanetary masses and radii have revealed a population of massive super-Earths --- planets sufficiently large that, according to one dimensional models, they should have turned into gas giants. To better understand the origin of these objects, we carry out hydrodynamical simulations of planetary cores embedded in a nascent protoplanetary disk. In this first paper of a series, to gain intuition as well as to develop useful diagnostics, we focus on two-dimensional simulations of the flow around protoplanetary cores. We use the pluto code to study isothermal and adiabatic envelopes around cores of sub- to super-thermal masses, fully resolving the envelope properties down to the core surface. Owing to the conservation of vortensity, envelopes acquire a substantial degree of rotational support when the core mass increases beyond the thermal mass, suggesting a limited applicability of one-dimensional models for describing the envelope structure. The finite size of the core (relatively large for super-Earths) also controls the amount of rotational support in the entire envelope. Steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores, triggering mass accretion and turbulent mixing in their interiors. We also examine the influence of the gas self-gravity on the envelope structure. Although it only weakly alters the properties of the envelopes, the gas gravity has significant effect on the properties of the density waves triggered by the core in the protoplanetary disk.

It has been recently suggested that the multiple concentric rings and gaps discovered by ALMA in many protoplanetary disks may be produced by a single planet, as a result of the complex propagation and dissipation of the multiple spiral density waves it excites in the disk. Numerical efforts to verify this idea have largely utilized the so-called locally isothermal approximation with a prescribed disk temperature profile. However, in protoplanetary disks this approximation does not provide an accurate description of the density wave dynamics on scales of tens of au. Moreover, we show that locally isothermal simulations tend to overestimate the contrast of ring and gap features, as well as misrepresent their positions, when compared to simulations in which the energy equation is evolved explicitly. This outcome is caused by the non-conservation of the angular momentum flux of linear perturbations in locally isothermal disks. We demonstrate this effect using simulations of locally isothermal and adiabatic disks (with essentially identical temperature profiles) and show how the dust distributions, probed by mm wavelength observations, differ between the two cases. Locally isothermal simulations may thus underestimate the masses of planets responsible for the formation of multiple gaps and rings on scales of tens of au observed by ALMA. We suggest that caution should be exercised in using the locally isothermal simulations to explore planet-disk interaction, as well as in other studies of wave-like phenomena in astrophysical disks.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Recent observations of protoplanetary disks, as well as simulations of planet-disk interaction, have suggested that a single planet may excite multiple spiral arms in the disk, in contrast to the previous expectations based on linear theory (predicting a one-armed density wave). We re-assess the origin of multiple arms in the framework of linear theory, by solving for the global two-dimensional response of a non-barotropic disk to an orbiting planet. We show that the formation of a secondary arm in the inner disk, at about half of the orbital radius of the planet, is a robust prediction of linear theory. This arm becomes stronger than the primary spiral at several tenths of the orbital radius of the planet. Several additional, weaker spiral arms may also form in the inner disk. On the contrary, a secondary spiral arm is unlikely to form in the outer disk. Our linear calculations, fully accounting for the global behavior of both the phases and amplitudes of perturbations, generally support the recently proposed WKB phase argument for the secondary arm origin (as caused by the intricacy of constructive interference of azimuthal harmonics of the perturbation at different radii). We provide analytical arguments showing that the process of a single spiral wake splitting up into multiple arms is a generic linear outcome of wave propagation in differentially rotating disks. It is not unique to planet-driven waves and occurs also in linear calculations of spiral wakes freely propagating with no external torques. These results are relevant for understanding formation of multiple rings and gaps in protoplanetary disks.

Observations of the first interstellar minor object 1I/2017 'Oumuamua did not reveal direct signs of outgassing that would have been natural if it had volatile-rich composition. However, a recent measurement by Micheli et al (2018) of a substantial non-gravitational acceleration affecting the orbit of this object has been interpreted as resulting from its cometary activity, which must be rather vigorous. Here we critically re-assess this interpretation by exploring the implications of measured non-gravitational acceleration for the 'Oumuamua's rotational state. We show that outgassing torques should drive rapid evolution of 'Oumuamua's spin (on a timescale of a few days), assuming torque asymmetry typical for the Solar System comets. However, given the highly elongated shape of the object, its torque asymmetry is likely higher, implying even faster evolution. This would have resulted in rapid rotational fission of 'Oumuamua during its journey through the Solar System and is clearly incompatible with the relative stability of its rotational state inferred from photometric variability. Based on these arguments, as well as the lack of direct signs of outgassing, we conclude that the classification of 'Oumuamua as a comet (invoked to explain its claimed anomalous acceleration) is questionable.

Motion of many comets is affected by non-gravitational forces caused by outgassing from their surfaces. Outgassing also produces reactive torques resulting in cometary spin evolution. We propose that the two processes are correlated and show that the change of cometary spin rate over its heliocentric orbit scales linearly with the amplitude of its non-gravitational acceleration. The proportionality constant depends on the comet size and orbital elements (semi-major axis and eccentricity) and on the (dimensionless) lever arm parameter $\zeta$ that relates the outgassing-induced torque and acceleration. We determine $\zeta$ for 7 comets for which both non-gravitational acceleration and change of spin period $\Delta P$ were measured and verify this relation. This sample spanning almost 4 decades in $\Delta P$ yields $\log\zeta=-2.21\pm 0.54$, surprisingly small value and spread. We then apply our framework to 209 comets with measured non-gravitational accelerations and determine the objects most likely to exhibit large spin period changes, $\Delta P\gtrsim 20$ min per orbit assuming rotation period of 10 hr and $\zeta$ comparable to our control sample. These objects should be primary targets for future studies of cometary spin variability, further constraining distribution of $\zeta$. Lack of comets with very high expected spin rate changes (which is not equivalent to having the highest non-gravitational acceleration) suggests that (1) cometary fission due to outgassing-driven spin-up must be an important process and (2) the distribution of $\zeta$ has a lower limit $\sim 10^{-3}$.

Massive eccentric disks (gaseous or particulate) orbiting a dominant central mass appear in many astrophysical systems, including planetary rings, protoplanetary and accretion disks in binaries, and nuclear stellar disks around supermassive black holes in galactic centers. We present an analytical framework for treating the nearly Keplerian secular dynamics of test particles driven by the gravity of an eccentric, apsidally aligned, zero-thickness disk with arbitrary surface density and eccentricity profiles. We derive a disturbing function describing the secular evolution of coplanar objects, which is explicitly related (via one-dimensional, convergent integrals) to the disk surface density and eccentricity profiles without using any ad hoc softening of the potential. Our analytical framework is verified via direct orbit integrations, which show it to be accurate in the low-eccentricity limit for a variety of disk models (for disk eccentricity < 0.1-0.2). We find that free precession in the potential of a disk with a smooth surface density distribution can naturally change from prograde to retrograde within the disk. Sharp disk features - edges and gaps - are the locations where this tendency is naturally enhanced, while the precession becomes very fast. Radii where free precession changes sign are the locations where substantial (formally singular) growth of the forced eccentricity of the orbiting objects occurs. Based on our results, we formulate a self-consistent analytical framework for computing an eccentricity profile for an aligned, eccentric disk (with a prescribed surface density profile) capable of precessing as a solid body under its own self-gravity.

We conducted high-contrast polarimetry observations of T Tau in the H-band, using the HiCIAO instrument mounted on the Subaru Telescope, revealing structures as near as 0.$\arcsec$1 from the stars T Tau N and T Tau S. The whole T Tau system is found to be surrounded by nebula-like envelopes, and several outflow-related structures are detected in these envelopes. We analyzed the detailed polarization patterns of the circumstellar structures near each component of this triple young star system and determined constraints on the circumstellar disks and outflow structures. We suggest that the nearly face-on circumstellar disk of T Tau N is no larger than 0.$\arcsec$8, or 117 AU, in the northwest, based on the existence of a hole in this direction, and no larger than 0.$\arcsec$27, or 40 AU, in the south. A new structure "N5" extends to about 0.$\arcsec$42, or 59 AU, on the southwest of the star, believed to be part of the disk. We suggest that T Tau S is surrounded by a highly inclined circumbinary disk with a radius of about 0.$\arcsec$3, or 44 AU, with a position angle of about 30$^\circ$, that is misaligned with the orbit of the T Tau S binary. After analyzing the positions and polarization vector patterns of the outflow-related structures, we suggest that T Tau S should trigger the well-known E-W outflow, and is also likely to be responsible for a southwest precessing outflow "coil" and a possible south outflow.

Discovery of the first interstellar asteroid (ISA) - 1I/2017 'Oumuamua - raised a number of questions regarding its origin. Many of them relate to its lack of cometary activity, suggesting refractory composition of 'Oumuamua. Here we explore the possibility that 'Oumuamua-like ISAs are produced in tidal disruption events (TDEs) of refractory planetoids (asteroids, dwarf planets, etc.) by the white dwarfs (WDs). This idea is supported by existing spectroscopic observations of metal-polluted WDs, hinting at predominantly volatile-poor composition of accreted material. We show that such TDEs sourced by realistic planetary systems (including a population of >1000 km planetoids and massive perturbers - Neptune-to-Saturn mass planets) can eject to interstellar space up to 30% of planetary mass involved in them. Collisional fragmentation, caused by convergent vertical motion of the disrupted planetoid's debris inside the Roche sphere of the WD, channels most of the original mass into 0.1-1 km fragments, similar to 'Oumuamua. Such size spectrum of ISAs (very different from the top-heavy distributions expected in other scenarios) implies that planetary TDEs can account for a significant fraction (up to ~30% under optimistic assumptions) of the ISAs. This figure is based on existing observations of WD metal pollution and accounts for observational biases by using realistic models of circum-WD planetary systems. ISAs should exhibit kinematic characteristics similar to old, dynamically hot Galactic populations; we interpret 'Oumuamua's slow Galactic motion as a statistical fluctuation. ISA ejection in individual planetary TDEs is highly anisotropic, resulting in large fluctuations of their space density. We also show that other ISA production mechanisms involving stellar remnants - direct ejection by massive planets around WDs and SN explosions - have difficulty explaining 'Oumuamua-like ISAs.

Spectroscopic observations of some metal-rich white dwarfs (WDs), believed to be polluted by planetary material, reveal the presence of compact gaseous metallic disks orbiting them. The observed variability of asymmetric, double-peaked emission line profiles in about half of such systems could be interpreted as the signature of precession of an eccentric gaseous debris disk. The variability timescales --- from decades down to $1.4$ yr (recently inferred for the debris disk around HE 1349--2305) --- are in rough agreement with the rate of general relativistic (GR) precession in the test particle limit. However, it has not been demonstrated that this mechanism can drive such a fast, coherent precession of a radially extended (out to $1 R_\odot$) gaseous disk mediated by internal stresses (pressure). Here we use the linear theory of eccentricity evolution in hydrodynamic disks to determine several key properties of eccentric modes in gaseous debris disks around WDs. We find a critical dependence of both the precession period and radial eccentricity distribution of the modes on the inner disk radius, $r_\mathrm{in}$. For small inner radii, $r_\mathrm{in} \lesssim (0.2 - 0.4) R_\odot$, the modes are GR-driven, with periods of $\approx 1 - 10$ yr. For $r_\mathrm{in} \gtrsim (0.2 - 0.4) R_\odot$, the modes are pressure-dominated, with periods of $\approx 3 - 20$ yr. Correspondence between the variability periods and inferred inner radii of the observed disks is in general agreement with this trend. In particular, the short period of HE 1349--2305 is consistent with its small $r_\mathrm{in}$. Circum-WD debris disks may thus serve as natural laboratories for studying the evolution of eccentric gaseous disks.

Spiral density waves are known to exist in many astrophysical disks, potentially affecting disk structure and evolution. We conduct a numerical study of the effects produced by a density wave, evolving into a shock, on the characteristics of the underlying disk. We measure the deposition of angular momentum in the disk by spiral shocks of different strength and verify the analytical prediction of Rafikov (2016) for the behavior of this quantity, using shock amplitude (which is potentially observable) as the input variable. Good agreement between the theory and numerics is found as we vary shock amplitude (including highly nonlinear shocks), disk aspect ratio, equation of state, radial profiles of the background density and temperature, and pattern speed of the wave. We show that high numerical resolution is required to properly capture shock-driven transport, especially at low wave amplitudes. We also demonstrate that relating local mass accretion rate to shock dissipation in rapidly evolving disks requires accounting for the time-dependent contribution to the angular momentum budget, caused by the time dependence of the radial pressure support. We provide a simple analytical prescription for the behavior of this contribution and demonstrate its excellent agreement with the simulation results. Using these findings we formulate a theoretical framework for studying one-dimensional (in radius) evolution of the shock-mediated accretion disks, which can be applied to a variety of astrophysical systems.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

Radial transport of particles, elements and fluid driven by internal stresses in three-dimensional (3D) astrophysical accretion disks is an important phenomenon, potentially relevant for the outward dust transport in protoplanetary disks, origin of the refractory particles in comets, isotopic equilibration in the Earth-Moon system, etc. To gain better insight into these processes, we explore the dependence of meridional circulation in 3D disks with shear viscosity on their thermal stratification, and demonstrate strong effect of the latter on the radial flow. Previous locally isothermal studies have normally found a pattern of the radial outflow near the midplane, switching to inflow higher up. Here we show, both analytically and numerically, that a flow, which is inward at all altitudes, is possible in disks with entropy and temperature steeply increasing with height. Such thermodynamic conditions may be typical in the optically thin, viscously heated accretion disks. Disks in which these conditions do not hold should feature radial outflow near the midplane, as long as their internal stress is provided by the shear viscosity. Our results can also be used for designing hydrodynamical disk simulations with a prescribed pattern of the meridional circulation.

Protoplanetary disks are believed to evolve on Myr timescales in a diffusive (viscous) manner as a result of angular momentum transport driven by internal stresses. Here we use a sample of 26 protoplanetary disks resolved by ALMA with measured (dust-based) masses and stellar accretion rates to derive the dimensionless $\alpha$-viscosity values for individual objects, with the goal of constraining the angular momentum transport mechanism. We find that the inferred values of $\alpha$ do not cluster around a single value, but instead have a broad distribution extending from $10^{-4}$ to $0.04$. Moreover, they correlate with neither the global disk parameters (mass, size, surface density) nor the stellar characteristics (mass, luminosity, radius). However, we do find a strong linear correlation between $\alpha$ and the central mass accretion rate $\dot M$. This correlation is unlikely to result from the direct physical effect of $\dot M$ on disk viscosity on global scales. Instead, we suggest that it is caused by the decoupling of stellar $\dot M$ from the global disk characteristics in one of the following ways. (1) The behavior (and range) of $\alpha$ is controlled by a yet unidentified parameter (e.g. ionization fraction, magnetic field strength, or geometry), ultimately driving the variation of $\dot M$. (2) The central $\dot M$ is decoupled from the global viscous mass accretion rate as a result of an instability or mass accumulation (or loss) in the inner disk. (3) Perhaps the most intriguing possibility is that angular momentum in protoplanetary disks is transported non-viscously, e.g. via magnetohydrodynamic winds or spiral density waves.

By performing non-masked polarization imaging with Subaru/HiCIAO, polarized scattered light from the inner region of the disk around the GG Tau A system was successfully detected in the $H$ band with a spatial resolution of approximately 0.07$\arcsec$, revealing the complicated inner disk structures around this young binary. This paper reports the observation of an arc-like structure to the north of GG Tau Ab and part of a circumstellar structure that is noticeable around GG Tau Aa extending to a distance of approximately 28 AU from the primary star. The speckle noise around GG Tau Ab constrains its disk radius to <13 AU. Based on the size of the circumbinary ring and the circumstellar disk around GG Tau Aa, the semi-major axis of the binary's orbit is likely to be 62 AU. A comparison of the present observations with previous ALMA and near-infrared (NIR) H$_2$ emission observations suggests that the north arc could be part of a large streamer flowing from the circumbinary ring to sustain the circumstellar disks. According to the previous studies, the circumstellar disk around GG Tau Aa has enough mass and can sustain itself for a duration sufficient for planet formation; thus, our study indicates that planets can form within close (separation $\lesssim$ 100 AU) young binary systems.

Supermassive black hole binaries (SMBHs) are expected to result from galaxy mergers, and thus are natural byproducts (and probes) of hierarchical structure formation in the Universe. They are also the primary expected source of low-frequency gravitational wave emission. We search for binary BHs using time-variable velocity shifts in broad Mg II emission lines of quasars with multi-epoch observations. First, we inspect velocity shifts of the binary SMBH candidates identified in Ju et al. (2013), using SDSS spectra with an additional epoch of data that lengthens the typical baseline to ~10 yr. We find variations in the line-of-sight velocity shifts over 10 years that are comparable to the shifts observed over 1-2 years, ruling out the binary model for the bulk of our candidates. We then analyze 1438 objects with 8 yr median time baselines, from which we would expect to see velocity shifts >1000 km/s from sub-pc binaries. We find only one object with an outlying velocity of 448 km/s, indicating, based on our modeling, that ~< 1 per cent (the value varies with different assumptions) of SMBHs that are active as quasars reside in binaries with ~0.1 pc separations. Binaries either sweep through these small separations rapidly or stall at larger radii.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

Several classes of stellar binaries with post-main sequence (post-MS) components - millisecond pulsars with the white dwarf companions (MSP+WD) and periods P_b~30 d, binaries hosting post-asymptotic giant branch (post-AGB) stars or barium stars with P_b~several yr - feature high eccentricities (up to 0.4) despite the expectation of their efficient tidal circularization during the AGB phase. It was suggested that the eccentricities of these binaries can be naturally excited by their tidal coupling to the circumbinary disk, formed by the material ejected from the binary. Here we critically reassess this idea using simple arguments rooted in the global angular momentum conservation of the disk+binary system. Compared to previous studies we (1) fully account for the viscous spreading of the circumbinary disk, (2) consider the possibility of reaccretion from the disk onto the binary (in agreement with simulations and empirical evidence), and (3) allow for the reduced viscosity after the disk expands, cools, and forms dust. These ingredients conspire to significantly lower the efficiency of eccentricity excitation by the disk tides. We find that explaining eccentricities of the post-AGB binaries is difficult and requires massive ($> 10^{-2}M_\odot$), long-lived ($> 10^5$ yr) circumbinary disks, that do not reaccrete. Reaccretion is detrimental for the eccentricity growth also in the MSP+WD systems. Reduced efficiency of the disk-driven excitation motivates study of alternative mechanisms for producing the peculiar eccentricities of the post-MS binaries.

Decretion (or external) disks are gas disks freely expanding to large radii due to their internal stresses. They are expected to naturally arise in tidal disruption events, around Be stars, in mass-losing post main sequence binaries, as a result of supernova fallback, etc. Their evolution is theoretically understood in two regimes: when the central object does not exert torque on the disk (a standard assumption for conventional accretion disks) or when no mass inflow (or outflow) occurs at the disk center. However, many astrophysical objects - circumbinary disks, Be stars, neutron stars accreting in a propeller regime, etc. - feature non-zero torque simultaneously with the non-zero accretion (or ejection of mass) at the disk center. We provide a general description for the evolution of such disks (both linear and non-linear) in the self-similar regime, to which the disk should asymptotically converge with time. We identify a similarity parameter $\lambda$, which is uniquely related to the degree, to which the central mass accretion is suppressed by the non-zero central torque. The known decretion disk solutions correspond to the two discrete values of $\lambda$, while our new solutions cover a continuum of its physically allowed values, corresponding to either accretion or mass ejection by the central object. A direct relationship between $\lambda$ and central $\dot M$ and torque is also established. We describe the time evolution of the various disk characteristics for different $\lambda$, and show that the observable properties (spectrum and luminosity evolution) of the decretion disks are in general different from the standard accretion disks with no central torque.

Many galaxies are expected to harbor binary supermassive black holes (SMBHs) in their centers. Their interaction with the surrounding gas results in accretion and exchange of angular momentum via tidal torques, facilitating binary inspiral. Here we explore the non-trivial coupling between these two processes and analyze how the global properties of externally supplied circumbinary disks depend on the binary accretion rate. By formulating our results in terms of the angular momentum flux driven by internal stresses, we come up with a very simple classification of the possible global disk structures, which differ from the standard constant $\dot M$ accretion disk solution. Suppression of accretion by the binary tides, leading to a significant mass accumulation in the inner disk, accelerates binary inspiral. We show that once the disk region strongly perturbed by the viscously transmitted tidal torque exceeds the binary semi-major axis, the binary can merge in less than its mass-doubling time due to accretion. Thus, unlike the inspirals driven by stellar scattering, the gas-assisted merger can occur even if the binary is embedded in a relatively low mass disk (lower than its own mass). This is important for resolving the "last parsec" problem for SMBH binaries and understanding powerful gravitational wave sources in the Universe. We argue that the enhancement of accretion by the binary found in some recent simulations cannot persist for a long time and should not affect the long-term orbital inspiral. We also review the existing simulations of the SMBH binary-disk coupling and propose a numerical setup, which is particularly well suited for verifying our theoretical predictions.

High-resolution imaging of some protoplanetary disks in scattered light reveals presence of the global spiral arms of significant amplitude, likely excited by massive planets or stellar companions. Assuming that these arms are density waves, evolving into spiral shocks, we assess their effect on the thermodynamics, accretion, and global evolution of the disk. We derive analytical expressions for the direct (irreversible) heating, angular momentum transport, and mass accretion rate induced by the disk shocks of arbitrary strength. We find these processes to be very sensitive to the shock amplitude. Focusing on the waves of moderate strength (density jump at the shock $\Delta\Sigma/\Sigma\sim 1$) we show the associated disk heating to be negligible (contributing at $\sim 1\%$ level to the energy budget) in passive, irradiated protoplanetary disks on $\sim 100$ AU scales, but becoming important within several AU from the star. At the same time, shock heating can be a significant (or even dominant) energy source in disks of cataclysmic variables, stellar X-ray binaries, and supermassive black hole binaries, heated mainly by viscous dissipation. Mass accretion induced by the global spiral shocks is comparable to (or exceeds) the mass inflow due to viscous stresses. Protoplanetary disks featuring prominent global spirals must be evolving rapidly, in $\lesssim 0.5$ Myr at $\sim 100$ AU. A direct upper limit on the disk evolution timescale can be established via the measurement of the gravitational torque due to the spiral arms from the imaging data. Our findings suggest that, regardless of their origin, global spiral waves must be important agents of the protoplanetary disk evolution. They may serve as an effective mechanism of disk dispersal and could be related to the transitional disk phenomenon.

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

Disk accretion at high rate onto a white dwarf or a neutron star has been suggested to result in the formation of a spreading layer (SL) - a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g. $\alpha$-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective re-arrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin-Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal non-uniformity of mixing in accreting white dwarfs may be linked to the observed bipolar morphology of classical novae ejecta.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

Recent discoveries of circumbinary planets by Kepler mission provide motivation for understanding their birthplaces - protoplanetary disks around stellar binaries with separations <1 AU. We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscous evolution of the disk properties, which are verified and calibrated using 1D numerical calculations with realistic inputs. Injection of angular momentum by the central binary suppresses mass accretion onto the binary and causes radial distribution of the viscous angular momentum flux F_J to be different from that in a standard accretion disk around a single star with no torque at the center. Disks with no mass accretion at the center develop F_J profile which is flat in radius. Radial profiles of temperature and surface density are also quite different from those in disks around single stars. Damping of the density waves driven by the binary and viscous dissipation dominate heating of the inner disk (within 1-2 AU), pushing the iceline beyond 3-5 AU, depending on disk mass and age. Irradiation by the binary governs disk thermodynamics beyond ~10 AU. However, self-shadowing by the hot inner disk may render central illumination irrelevant out to ~20 AU. Spectral energy distribution of a circumbinary disk exhibits a distinctive bump around 10 micron, which may facilitate identification of such disks around unresolved binaries. Efficient tidal coupling to the disk drives orbital inspiral of the binary and may cause low-mass and compact binaries to merge into a single star within the disk lifetime. We generally find that circumbinary disks present favorable sites for planet formation (despite wider zone of volatile depletion), in agreement with the statistics of Kepler circumbinary planets.

Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.

The Kepler mission has discovered about a dozen circumbinary planetary systems, all containing planets on ~ 1 AU orbits. We place bounds on the locations in the circumbinary protoplanetary disk, where these planets could have formed through collisional agglomeration starting from small (km-sized or less) planetesimals. We first present a model of secular planetesimal dynamics that accounts for the (1) perturbation due to the eccentric precessing binary, as well as the (2) gravity and (3) gas drag from a precessing eccentric disk. Their simultaneous action leads to rich dynamics, with (multiple) secular resonances emerging in the disk. We derive analytic results for size-dependent planetesimal eccentricity, and demonstrate the key role of the disk gravity for circumbinary dynamics. We then combine these results with a simple model for collisional outcomes and find that in systems like Kepler 16, planetesimal growth starting with 10-100 m planetesimals is possible outside a few AU. The exact location exterior to which this happens is sensitive to disk eccentricity, density and precession rate, as well as to the size of the first generation of planetesimals. Strong perturbations from the binary in the inner part of the disk, combined with a secular resonance at a few AU inhibit the growth of km-sized planetesimals within 2 - 4 AU of the binary. In situ planetesimal growth in the Kepler circumbinary systems is possible only starting from large (few-km-sized) bodies in a low-mass disk with surface density less than 500 g/cm^2 at 1 AU.

Gravitationally unstable accretion disks emerge in a variety of astrophysical contexts - giant planet formation, FU Orioni outbursts, feeding of AGNs, and the origin of Pop III stars. When a gravitationally unstable disk is unable to cool rapidly it settles into a quasi-stationary, fluctuating gravitoturbulent state, in which its Toomre Q remains close to a constant value Q_0~1. Here we develop an analytical formalism describing the evolution of such a disk, which is based on the assumptions of Q=Q_0 and local thermal equilibrium. Our approach works in the presence of additional sources of angular momentum transport (e.g. MRI), as well as external irradiation. Thermal balance dictates a unique value of the gravitoturbulent stress \alpha_{gt} driving disk evolution, which is a function of the local surface density and angular frequency. We compare this approach with other commonly used gravitoturbulent viscosity prescriptions, which specify the explicit dependence of stress \alpha_{gt} on Toomre Q in an ad hoc fashion, and identify the ones that provide consistent results. We nevertheless argue that our Q=Q_0 approach is more flexible, robust, and straightforward, and should be given preference in applications. We illustrate this with a couple of analytical calculations - locations of the snow line and of the outer edge of the dead zone in a gravitoturbulent protoplanetary disk - which clearly show the simplicity and versatility of the Q=Q_0 approach.

Detections of planets in eccentric, close (separations of ~20 AU) binary systems such as \alpha Cen or \gamma Cep provide an important test of planet formation theories. Gravitational perturbations from the companion are expected to excite high planetesimal eccentricities resulting in destruction, rather than growth, of objects with sizes of up to several hundred km in collisions of similar-size bodies. It was recently suggested that gravity of a massive axisymmetric gaseous disk in which planetesimals are embedded drives rapid precession of their orbits, suppressing eccentricity excitation. However, disks in binaries are themselves expected to be eccentric, leading to additional planetesimal excitation. Here we develop secular theory of eccentricity evolution for planetesimals perturbed by the gravity of an elliptical protoplanetary disk (neglecting gas drag) and the companion. For the first time we derive an expression for the disturbing function due to an eccentric disk, which can be used for a variety of other astrophysical problems. We obtain explicit analytical solutions for planetesimal eccentricity evolution and delineate four different regimes of dynamical excitation. We show that in systems with massive (>10^{-2}M_Sun) disks, planetesimal eccentricity is usually determined by the gravity of the eccentric disk alone, and is comparable to the disk eccentricity. As a result, the latter imposes a lower limit on collisional velocities of solids, making their growth problematic. This fragmentation barrier can be removed if the gaseous disk rapidly precesses or if its own self-gravity is efficient at lowering disk eccentricity.

About $20\%$ of exoplanets discovered by radial velocity surveys reside in stellar binaries. To clarify their origin one has to understand the dynamics of planetesimals in protoplanetary disks within binaries. The standard description, accounting for only gas drag and gravity of the companion star has been challenged recently, as the gravity of the protoplanetary disk was shown to play a crucial role in planetesimal dynamics. An added complication is the tendency of protoplanetary disks in binaries to become eccentric, giving rise to additional excitation of planetesimal eccentricity. Here, for the first time, we analytically explore secular dynamics of planetesimals in binaries such as $\alpha$ Cen and $\gamma$ Cep under the combined action of (1) gravity of the eccentric protoplanetary disk, (2) perturbations due to the (coplanar) eccentric companion, and (3) gas drag. We derive explicit solutions for the behavior of planetesimal eccentricity ${\bf e}_p$ in non-precessing disks (and in precessing disks in certain limits). We obtain the analytical form of the distribution of relative velocities of planetesimals, which is a key input for understanding their collisional evolution. Disk gravity strongly influences relative velocities and tends to push sizes of planetesimals colliding with comparable objects at the highest speed to small values, $\sim 1$ km. We also find that planetesimals in eccentric protoplanetary disks apsidally aligned with the binary orbit collide at lower relative velocities than in mis-aligned disks. Our results highlight a decisive role that disk gravity plays in planetesimal dynamics in binaries.

Planet formation in small-separation (~20 AU) eccentric binaries such as gamma Cephei or alpha Centauri is believed to be adversely affected by the presence of the stellar companion. Strong dynamical excitation of planetesimals by the eccentric companion can result in collisional destruction (rather than growth) of 1-100 km objects, giving rise to the "fragmentation barrier" for planet formation. We revise this issue using a novel description of secular dynamics of planetesimals in binaries, which accounts for the gravity of the eccentric, coplanar protoplanetary disk, as well as gas drag. By studying planetesimal collision outcomes we show, in contrast to many previous studies, that planetesimal growth and subsequent formation of planets (including gas giants) in AU-scale orbits within ~20 AU separation binaries may be possible, provided that the protoplanetary disks are massive (>10^{-2}M_\odot) and only weakly eccentric (disk eccentricity <0.01). These requirements are compatible with both the existence of massive (several M_J) planets in gamma Cep-like systems and the results of recent simulations of gaseous disks in eccentric binaries. Terrestrial and Neptune-like planets can also form in lower-mass disks at small (sub-AU) radii. We find that fragmentation barrier is less of a problem in eccentric disks which are apsidally aligned with the binary orbit. Alignment gives rise to special locations, where (1) relative planetesimal velocities are low and (2) the timescale of their drag-induced radial drift is long. This causes planetesimal pileup at such locations in the disk and promotes their growth.

Many exoplanets in close-in orbits are observed to have relatively high eccentricities and large stellar obliquities. We explore the possibility that these result from planet-planet scattering by studying the dynamical outcomes from a large number of orbit integrations in systems with two and three gas-giant planets in close-in orbits (0.05 AU < a < 0.15 AU). We find that at these orbital separations, unstable systems starting with low eccentricities and mutual inclinations ($e\lesssim0.1$, $i\lesssim0.1$) generally lead to planet-planet collisions in which the collision product is a planet on a low-eccentricity, low-inclination orbit. This result is inconsistent with the observations. We conclude that eccentricity and inclination excitation from planet-planet scattering must precede migration of planets into short-period orbits. This result constrains theories of planet migration: the semi-major axis must shrink by 1-2 orders of magnitude without damping the eccentricity and inclination.

We explore secular dynamics of a recently discovered hierarchical triple system consisting of the radio pulsar PSR J0337+1715 and two white dwarfs (WDs). We show that three body interactions endow the inner binary with a large forced eccentricity and suppress its apsidal precession, to about 24% of the rate due to the general relativity. However, precession rate is still quite sensitive to the non-Newtonian effects and may be used to constrain gravity theories if measured accurately. Small value of the free eccentricity of the inner binary $e_{i}^{free}\approx 2.6\times 10^{-5}$ and vanishing forced eccentricity of the outer, relatively eccentric binary naturally result in their apsidal near-alignment. In addition, this triple system provides a unique opportunity to explore excitation of both eccentricity and inclination in neutron star-WD binaries, e.g. due to random torques caused by convective eddies in the WD progenitor. We show this process to be highly anisotropic and more effective at driving eccentricity rather than inclination. The outer binary eccentricity as well as $e_{i}^{free}$ exceed by more than an order of magnitude the predictions of the eccentricity-period relation of Phinney (1992), which is not uncommon. We also argue that the non-zero mutual inclination of the two binaries emerges at the end of the Roche lobe overflow of the outer (rather than the inner) binary.

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{\oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<\Omega t_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $\Omega t_{s}\sim 1$. The dust concentration is highly elongated in the $\phi$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($\Omega t_{s}\gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.

Understanding the dominant brown dwarf and giant planet formation processes, and finding out whether these processes rely on completely different mechanisms or share common channels represents one of the major challenges of astronomy and remains the subject of heated debates. It is the aim of this review to summarize the latest developments in this field and to address the issue of origin by confronting different brown dwarf and giant planet formation scenarios to presently available observational constraints. As examined in the review, if objects are classified as "Brown Dwarfs" or "Giant Planets" on the basis of their formation mechanism, it has now become clear that their mass domains overlap and that there is no mass limit between these two distinct populations. Furthermore, while there is increasing observational evidence for the existence of non-deuterium burning brown dwarfs, some giant planets, characterized by a significantly metal enriched composition, might be massive enough to ignite deuterium burning in their core. Deuterium burning (or lack of) thus plays no role in either brown dwarf or giant planet formation. Consequently, we argue that the IAU definition to distinguish these two populations has no physical justification and brings scientific confusion. In contrast, brown dwarfs and giant planets might bear some imprints of their formation mechanism, notably in their mean density and in the physical properties of their atmosphere. Future direct imaging surveys will undoubtedly provide crucial information and perhaps provide some clear observational diagnostics to unambiguously distinguish these different astrophysical objects.

The Oort Cloud remains one of the most poorly explored regions of the Solar System. We propose that its properties can be constrained by studying a population of dust grains produced in collisions of comets in the outer Solar System. We explore the dynamics of micron-size grains outside the heliosphere (beyond ~250 AU), which are affected predominantly by the magnetic field of the interstellar medium (ISM) flow past the Sun. We derive analytic models for the production and motion of small particles as a function of their birth location in the Cloud and calculate the particle flux and velocity distribution in the inner Solar System. These models are verified by direct numerical simulations. We show that grains originating in the Oort Cloud have a unique distribution of arrival directions, which should easily distinguish them from both interplanetary and interstellar dust populations. We also demonstrate that the distribution of particle arrival velocities is uniquely determined the mass distribution and dust production rate in the Cloud. Cometary collisions within the Cloud produce a flux of micron-size grains in the inner Solar System of up to several per square meter per year. The next-generation dust detectors may be sensitive enough to detect and constrain this dust population, which will illuminate us about the Oort Cloud's properties. We also show that the recently-detected mysterious population of large (micron-size) unbound particles, which seems to arrive with the ISM flow is unlikely to be generated by collisions of comets in the Oort Cloud.

Supermassive black hole (SMBH) binaries are expected in a Lambda CDM cosmology given that most (if not all) massive galaxies contain a massive black hole at their center. So far, however, direct evidence for such binaries has been elusive. We use cross-correlation to search for temporal velocity shifts in the MgII broad emission lines of 0.36 < z < 2 quasars with multiple observations in the Sloan Digital Sky Survey. For ~ 10^9 Msun BHs in SMBH binaries, we are sensitive to velocity drifts for binary separations of ~ 0.1 pc with orbital periods of ~100 years. We find seven candidate sub-pc--scale binaries with velocity shifts > 3.4 sigma ~ 280 km/s, where sigma is our systematic error. Comparing the detectability of SMBH binaries with the number of candidates (N < 7), we can rule out that most 10^9 Msun BHs exist in ~ 0.03-0.2 pc scale binaries, in a scenario where binaries stall at sub-pc scales for a Hubble time. We further constrain that < one-third of quasars host SMBH binaries after considering gas-assisted sub-pc evolution of SMBH binaries, although this result is very sensitive to the assumed size of the broad line region. We estimate the detectability of SMBH binaries with ongoing or next-generation surveys (e.g., BOSS, Subaru Prime Focus Spectrograph), taking into account the evolution of the sub-parsec binary in circumbinary gas disks. These future observations will provide longer time baselines for searches similar to ours and may in turn constrain the evolutionary scenarios of SMBH binaries.

Eclipsing binary DI Herculis (DI Her) is known to exhibit anomalously slow apsidal precession, below the rate predicted by the general relativity. Recent measurements of the Rossiter-McLauglin effect indicate that stellar spins in DI Her are almost orthogonal to the orbital angular momentum, which explains the anomalous precession in agreement with the earlier theoretical suggestion by Shakura. However, these measurements yield only the projections of the spin-orbit angles onto the sky plane, leaving the spin projection onto our line of sight unconstrained. Here we describe a method of determining the full three-dimensional spin orientation of the binary components relying on the use of the gravity darkening effect, which is significant for the rapidly rotating stars in DI Her. Gravity darkening gives rise to nonuniform brightness distribution over the stellar surface, the pattern of which depends on the stellar spin orientation. Using archival photometric data obtained during multiple eclipses spread over several decades we are able to constrain the unknown spin angles in DI Her with this method, finding that spin axes of both stars lie close to the plane of the sky. Our procedure fully accounts for the precession of stellar spins over the long time span of observations.

We study wakes and gap opening by low mass planets in gaseous protoplanetary disks threaded by net vertical magnetic fields which drive magnetohydrodynamical (MHD) turbulence through the magnetorotational instabilty (MRI), using three dimensional simulations in the unstratified local shearing box approximation. The wakes, which are excited by the planets, are damped by shocks similar to the wake damping in inviscid hydrodynamic (HD) disks. Angular momentum deposition by shock damping opens gaps in both MHD turbulent disks and inviscid HD disks even for low mass planets, in contradiction to the "thermal criterion" for gap opening. To test the "viscous criterion", we compared gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective \alpha within the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. We also confirmed the large excess torque close to the planet in MHD disks, and found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. The comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant \alpha-viscosity to model gaps in protoplanetary disks.

We perform global unstratified 3D magnetohydrodynamic simulations of an astrophysical boundary layer (BL) -- an interface region between an accretion disk and a weakly magnetized accreting object such as a white dwarf -- with the goal of understanding the effects of magnetic field on the BL. We use cylindrical coordinates with an isothermal equation of state and investigate a number of initial field geometries including toroidal, vertical, and vertical with zero net flux. Our initial setup consists of a Keplerian disk attached to a non-rotating star. In a previous work, we found that in hydrodynamical simulations, sound waves excited by shear in the BL were able to efficiently transport angular momentum and drive mass accretion onto the star. Here we confirm that in MHD simulations, waves serve as an efficient means of angular momentum transport in the vicinity of the BL, despite the magnetorotational instability (MRI) operating in the disk. In particular, the angular momentum current due to waves is at times larger than the angular momentum current due to MRI. Our results suggest that angular momentum transport in the BL and its vicinity is a global phenomenon occurring through dissipation of waves and shocks. This point of view is quite different from the standard picture of transport by a local anomalous turbulent viscosity. In addition to angular momentum transport, we also study magnetic field amplification within the BL. We find that the field is indeed amplified in the BL, but only by a factor of a few and remains subthermal.

Circumbinary planetary systems recently discovered by Kepler represent an important testbed for planet formation theories. Planetesimal growth in disks around binaries has been expected to be inhibited interior to ~10 AU by secular excitation of high relative velocities between planetesimals, leading to their collisional destruction (rather than agglomeration). Here we show that gravity of the gaseous circumbinary disk in which planets form drives fast precession of both the planetesimal and binary orbits, resulting in strong suppression of planetesimal eccentricities beyond 2-3 AU and making possible growth of 1-100 km objects in this region. The precise location of the boundary of accretion-friendly region depends on the size of the inner disk cavity cleared by the binary torques and on the disk mass (even 0.01 M_Sun disk strongly suppresses planetesimal excitation), among other things. Precession of the orbit of the central binary, enhanced by the mass concentration naturally present at the inner edge of a circumbinary disk, plays key role in this suppression, which is a feature specific to the circumbinary planet formation.

Existence of planets is binaries with relatively small separations (around 20 AU), such as \alpha Centauri or \gamma Cephei poses severe challenges to standard planet formation theories. The problem lies in the vigorous secular excitation of planetesimal eccentricities at separations of several AU, where some of the planets are found, by the massive, eccentric stellar companions. High relative velocities of planetesimals preclude their growth in mutual collisions for a wide range of sizes, from below 1 km up to several hundred km, resulting in fragmentation barrier to planet formation. Here we show that rapid apsidal precession of planetesimal orbits, caused by the gravity of the circumstellar protoplanetary disk, acts to strongly reduce eccentricity excitation, lowering planetesimal velocities by an order of magnitude or even more at 1 AU. By examining the details of planetesimal dynamics we demonstrate that this effect eliminates fragmentation barrier for in-situ growth of planetesimals as small as < 10 km even at separations as wide as 2.6 AU (semi-major axis of the giant planet in HD 196885), provided that the circumstellar protoplanetary disk is relatively massive, ~0.1 M_Sun.

The nature of angular momentum transport in the boundary layers of accretion disks has been one of the central and long-standing issues of accretion disk theory. In this work we demonstrate that acoustic waves excited by supersonic shear in the boundary layer serve as an efficient mechanism of mass, momentum and energy transport at the interface between the disk and the accreting object. We develop the theory of angular momentum transport by acoustic modes in the boundary layer, and support our findings with 3D hydrodynamical simulations, using an isothermal equation of state. Our first major result is the identification of three types of global modes in the boundary layer. We derive dispersion relations for each of these modes that accurately capture the pattern speeds observed in simulations to within a few percent. Second, we show that angular momentum transport in the boundary layer is intrinsically nonlocal, and is driven by radiation of angular momentum away from the boundary layer into both the star and the disk. The picture of angular momentum transport in the boundary layer by waves that can travel large distances before dissipating and redistributing angular momentum and energy to the disk and star is incompatible with the conventional notion of local transport by turbulent stresses. Our results have important implications for semianalytical models that describe the spectral emission from boundary layers.

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

Disk accretion onto a weakly magnetized central object, e.g. a star, is inevitably accompanied by the formation of a boundary layer near the surface, in which matter slows down from the highly supersonic orbital velocity of the disk to the rotational velocity of the star. We perform high resolution 2D hydrodynamical simulations in the equatorial plane of an astrophysical boundary layer with the goal of exploring the dynamics of non-axisymmetric structures that form there. We generically find that the supersonic shear in the boundary layer excites non-axisymmetric quasi-stationary acoustic modes that are trapped between the surface of the star and a Lindblad resonance in the disk. These modes rotate in a prograde fashion, are stable for hundreds of orbital periods, and have a pattern speed that is less than and of order the rotational velocity at the inner edge of the disk. The origin of these intrinsically global modes is intimately related to the operation of a corotation amplifier in the system. Dissipation of acoustic modes in weak shocks provides a universal mechanism for angular momentum and mass transport even in purely hydrodynamic (i.e. non-magnetized) boundary layers. We discuss the possible implications of these trapped modes for explaining the variability seen in accreting compact objects.

We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.

A number of metal-rich white dwarfs (WDs) are known to host compact, dense particle disks, which are thought to be responsible for metal pollution of these stars. In many such systems the inner radii of disks inferred from their spectra are so close to the WD that particles directly exposed to starlight must be heated above 1500 K and are expected to be unstable against sublimation. To reconcile this expectation with observations we explore particle sublimation in H-poor debris disks around WDs. We show that because of the high metal vapor pressure the characteristic sublimation temperature in these disks is 300-400 K higher than in their protoplanetary analogues, allowing particles to survive at higher temperatures. We then look at the structure of the inner edges of debris disks and show that they should generically feature superheated inner rims directly exposed to starlight with temperatures reaching 2500-3500 K. Particles migrating through the rim towards the WD (and rapidly sublimating) shield the disk behind them from strong stellar heating, making the survival of solids possible close to the WD. Our model agrees well with observations of WD+disk systems provided that disk particles are composed of Si-rich material such as olivine, and have sizes in the range ~(0.03-30) cm.

It is generally believed that gaseous disks around supermassive black hole (SMBH) binaries in centers of galaxies can facilitate binary merger and give rise to observational signatures both in electromagnetic and gravitational wave domains. We explore general properties of circumbinary disks by reformulating standard equations for the viscous disk evolution in terms of the viscous angular momentum flux F_J. In steady state F_J is a linear function of the specific angular momentum, which is a generalization of (but is not equivalent to) the standard constant \dot M disk solution. If the torque produced by the central binary is effective at stopping gas inflow and opening a gap (or cavity) in the disk, then the inner part of the circumbinary disk can be approximated as a constant F_J disk. We compute properties of such disks in different physical regimes relevant for SMBH binaries and use these results to understand the gas-assisted evolution of SMBH pairs starting at separations 10^{-4}-10^{-2} pc. We find the following. (1) Pile-up of matter at the inner edge of the disk leads to continuous growth of the torque acting on the binary and can considerably accelerate its orbital evolution compared to the gravitational wave-driven decay. (2) Torque on the binary is determined non-locally and does not in general reflect the disk properties in the vicinity of the binary. (3) Binary evolution depends on the past history of the disk evolution. (4) Eddington limit can be important in circumbinary disks even if they accrete at sub-Eddington rates at late stages of binary evolution. (5) Circumbinary disk self-consistently evolved under the action of the binary torque emits more power and has spectrum different from the spectrum of constant \dot M disk - it is steeper (\nu F_\nu\propto \nu^{12/7}) and extends to shorter wavelength, facilitating its detection.

Gravitational coupling between a gaseous disk and an orbiting perturber leads to angular momentum exchange between them which can result in gap opening by planets in protoplanetary disks and clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks. Understanding the co-evolution of the disk and the orbit of the perturber in these circumstances requires knowledge of the spatial distribution of the torque exerted by the latter on a highly nonuniform disk. Here we explore disk-satellite interaction in disks with gaps in linear approximation both in Fourier and in physical space, explicitly incorporating the disk non-uniformity in the fluid equations. Density gradients strongly displace the positions of Lindblad resonances in the disk (which often occur at multiple locations), and the waveforms of modes excited close to the gap edge get modified compared to the uniform disk case. The spatial distribution of the excitation torque density is found to be quite different from the existing prescriptions: most of the torque is exerted in a rather narrow region near the gap edge where Lindblad resonances accumulate, followed by an exponential fall-off with the distance from the perturber. Despite these differences, for a given gap profile the full integrated torque exerted on the disk agrees with the conventional uniform disk theory prediction at the level of ~10%. The nonlinearity of the density wave excited by the perturber is shown to decrease as the wave travels out of the gap, slowing down its nonlinear evolution and damping. Our results suggest that gap opening in protoplanetary disks and gas clearing around SMBH binaries can be more efficient than the existing theories predict. They pave the way for self-consistent calculations of the gap structure and the orbital evolution of the perturber using accurate prescription for the torque density behavior.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

A growing sample of white dwarfs (WDs) with metal-enriched atmospheres are accompanied by excess infrared emission, indicating that they are encircled by a compact dusty disk of solid debris. Such `WD debris disks' are thought to originate from the tidal disruption of asteroids or other minor bodies, but the precise mechanism(s) responsible for transporting matter to the WD surface remains unclear, especially in those systems with the highest inferred metal accretion rates dM_Z/dt ~ 1e8-1e10 g/s. Here we present global time-dependent calculations of the coupled evolution of the gaseous and solid components of WD debris disks. Solids transported inwards (initially due to PR drag) sublimate at tens of WD radii, producing a source of gas that accretes onto the WD surface and viscously spreads outwards in radius, where it overlaps with the solid disk. If the aerodynamic coupling between the solids and gaseous disks is sufficiently strong (and/or the gas viscosity sufficiently weak), then gas builds up near the sublimation radius faster than it can viscously spread away. Since the rate of drag-induced solid accretion increases with gas density, this results in a runaway accretion process, during which the WD accretion rate reaches values orders of magnitude higher than can be achieved by PR drag alone. We explore the evolution of WD debris disks across a wide range of physical conditions and calculate the predicted distribution of observed accretion rates dM_Z/dt, finding reasonable agreement with the current sample. Although the conditions necessary for runaway accretion are at best marginally satisfied given the minimal level of aerodynamic drag between circular gaseous and solid disks, the presence of other stronger forms of solid-gas coupling---such as would result if the gaseous disk is only mildly eccentric---substantially increase the likelihood of runaway accretion.

Disk accretion onto weakly magnetized astrophysical objects often proceeds via a boundary layer that forms near the object's surface, in which the rotation speed of the accreted gas changes rapidly. Here we study the initial stages of formation for such a boundary layer around a white dwarf or a young star by examining the hydrodynamical shear instabilities that may initiate mixing and momentum transport between the two fluids of different densities moving supersonically with respect to each other. We find that an initially laminar boundary layer is unstable to two different kinds of instabilities. One is an instability of a supersonic vortex sheet (implying a discontinuous initial profile of the angular speed of the gas) in the presence of gravity, which we find to have a growth rate of order (but less than) the orbital frequency. The other is a sonic instability of a finite width, supersonic shear layer, which is similar to the Papaloizou-Pringle instability. It has a growth rate proportional to the shear inside the transition layer, which is of order the orbital frequency times the ratio of stellar radius to the boundary layer thickness. For a boundary layer that is thin compared to the radius of the star, the shear rate is much larger than the orbital frequency. Thus, we conclude that sonic instabilities play a dominant role in the initial stages of nonmagnetic boundary layer formation and give rise to very fast mixing between disk gas and stellar fluid in the supersonic regime.

Tidal interaction between a gaseous disk and a massive orbiting perturber is known to result in angular momentum exchange between them. Understanding astrophysical manifestations of this coupling such as gap opening by planets in protoplanetary disks or clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks requires knowledge of the spatial distribution of the torque exerted on the disk by a perturber. Recent hydrodynamical simulations by Dong et al (2011) have shown evidence for the tidal torque density produced in a uniform disk to change sign at the radial separation of $\approx 3.2$ scale heights from the perturber's orbit, in clear conflict with the previous studies. To clarify this issue we carry out a linear calculation of the disk-satellite interaction putting special emphasis on understanding the behavior of the perturbed fluid variables in physical space. Using analytical as well as numerical methods we confirm the reality of the negative torque density phenomenon and trace its origin to the overlap of Lindblad resonances in the vicinity of the perturber's orbit - an effect not accounted for in previous studies. These results suggest that calculations of the gap and cavity opening in disks by planets and binary SMBHs should rely on more realistic torque density prescriptions than the ones used at present.

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

Recent infrared observations have revealed presence of compact (radii < R_Sun) debris disks around more than a dozen of metal-rich white dwarfs (WD), likely produced by tidal disruption of asteroids. Accretion of high-Z material from these disks may account for the metal contamination of these WDs. It was previously shown using local calculations that the Poynting-Robertson (PR) drag acting on the dense, optically thick disk naturally drives metal accretion onto the WD at the typical rate \dot M_{PR} \approx 10^8 g/s. Here we extend this local analysis by exploring global evolution of the debris disk under the action of the PR drag for a variety of assumptions about the disk properties. We find that massive disks (mass > 10^{20} g), which are optically thick to incident stellar radiation inevitably give rise to metal accretion at rates \dot M > 0.2\dot M_{PR}. The magnitude of \dot M and its time evolution are determined predominantly by the initial pattern of the radial distribution of the debris (i.e. ring-like vs. disk-like) but not by the total mass of the disk. The latter determines only the disk lifetime, which can be several Myr or longer. Evolution of an optically thick disk generically results in the development of a sharp outer edge of the disk. We also find that the low mass (< 10^{20} g), optically thin disks exhibit \dot M << \dot M_{PR} and evolve on characteristic timescale \sim 10^5-10^6 yr, independent of their total mass.

Collisions resulting in fragmentation are important in shaping the mass spectrum of minor bodies in the asteroid belt, the Kuiper belt, and debris disks. Models of fragmentation cascades typically find that in steady-state, the solution for the particle mass distribution is a power law in the mass. However, previous studies have typically assumed that the mass of the largest fragment produced in a collision with just enough energy to shatter the target and disperse half its mass to infinity is directly proportional to the target mass. We show that if this assumption is not satisfied, then the power law solution for the steady-state particle mass distribution is modified by a multiplicative factor, which is a slowly varying function of the mass. We derive analytic solutions for this correction factor and confirm our results numerically. We find that this correction factor proves important when extrapolating over many orders of magnitude in mass, such as when inferring the number of large objects in a system based on infrared observations. In the course of our work, we have also discovered an unrelated type of non-power law behavior: waves can persist in the mass distribution of objects even in the absence of upper or lower cutoffs to the mass distribution or breaks in the strength law.

It was recently proposed that metal-rich white dwarfs (WDs) accrete their metals from compact debris disks found to exist around more than a dozen of them. At the same time, elemental abundances measured in atmospheres of some WDs imply vigorous metal accretion at rates up to $10^{11}$ g/s, far in excess of what can be supplied solely by Poynting-Robertson drag acting on such debris disks. To explain this observation we propose a model, in which rapid transport of metals from the disk onto the WD naturally results from interaction between this particulate disk and spatially coexisting disk of metallic gas. The latter is fed by evaporation of debris particles at the sublimation radius located at several tens of WD radii. Because of pressure support gaseous disk orbits WD slower than particulate disk. Resultant azimuthal drift between them at speed ~1 m/s causes aerodynamic drag on the disk of solids and drives inward migration of its constituent particles. Upon reaching the sublimation radius particles evaporate, enhancing the density of metallic gaseous disk and leading to positive feedback. Under favorable circumstances (low viscosity in the disk of metallic gas and efficient aerodynamic coupling between the disks) system evolves in a runaway fashion, destroying debris disk on time scale of $\sim 10^5$ yr, and giving rise to high metal accretion rates up to $10^{10}-10^{11}$ g/s, in agreement with observations.

Recent discoveries of compact (sizes $<R_\odot$) debris disks around more than a dozen of metal-rich white dwarfs (WDs) suggest that pollution of these stars with metals may be caused by accretion of high-Z material from the disk. But the mechanism responsible for efficient transfer of mass from a particulate disk to the WD atmosphere has not yet been identified. Here we demonstrate that radiation of the WD can effectively drive accretion of matter through the disk towards the sublimation radius (located at several tens of WD radii), where particles evaporate, feeding a disk of metal gas accreting onto the WD. We show that, contrary to some previous claims, Poynting-Robertson (PR) drag on the debris disk is effective at providing metal accretion rate $\dot M_{PR}\sim 10^8$ g/s and higher, scaling quadratically with WD effective temperature. We compare our results with observations and show that, as expected, no WD hosting a particulate debris disk shows evidence of metal accretion rate below that produced by the PR drag. Existence of WDs accreting metals at rates significantly higher than $\dot M_{PR}$ suggests that another mechanism in addition to PR drag drives accretion of high-Z elements in these systems.

We study the dynamics of large dust grains >1 micron with orbits outside of the heliosphere (beyond 250 AU). Motion of the Solar System through the interstellar medium (ISM) at a velocity of 26 km/s subjects these particles to gas and Coulomb drag (grains are expected to be photoelectrically charged) as well as the Lorentz force and the electric force caused by the induction electric field. We show that to zeroth order the combined effect of these forces can be well described in the framework of the classical Stark problem: particle motion in a Keplerian potential subject to an additional constant force. Based on this analogy, we elucidate the circumstances in which the motion becomes unbound, and show that under local ISM conditions dust grains smaller than ~100 microns originating in the Oort Cloud (e.g. in collisions of comets) beyond 10000 AU are ejected from the Solar System under the action of the electric force. Orbital motion of larger, bound grains is described analytically using the orbit-averaged Hamiltonian approach and consists of orbital plane precession at a fixed semi-major axis, accompanied by the periodic variations of the inclination and eccentricity (the latter may approach unity in some cases). A more detailed analysis of the combined effect of gas and Coulomb drag shows it is possible to reduce particle semi-major axes, but that the degree of orbital decay is limited (a factor of several at best) by passages through atomic and molecular clouds, which easily eject small particles.

The issue of giant planet formation by core instability (CI) far from the central star is rather controversial because the growth of massive solid core necessary for triggering the CI can take longer than the lifetime of the protoplanetary disk. In this work we assess the range of separations at which the CI may operate by (1) allowing for arbitrary (physically meaningful) rate of planetesimal accretion by the core and (2) properly taking into account the dependence of the critical mass for the CI on the planetesimal accretion luminosity. This self-consistent approach distinguishes our work from similar studies in which only a specific planetesimal accretion regime was explored and/or the critical core mass was fixed at some arbitrary level. We demonstrate that the largest separation at which the CI can occur within 3 Myr corresponds to the surface density of solids in the disk higher than 0.1 g cm^{-2} and is 40-50 AU in the minimum mass Solar nebula. This limiting separation is achieved when the planetesimal accretion proceeds at the fastest possible rate, even though the high associated accretion luminosity increases the critical core mass delaying the onset of the CI. Our constraints are independent of the mass of the central star and vary only weakly with the core density and its atmospheric opacity. We also discuss various factors which can strengthen or weaken our limits on the operation of the CI.

We study the dynamics of a vertically thin, dispersion-dominated disk of planetesimals with eccentricities $e$ and inclinations $i$ (normalized in Hill units) satisfying $e >> 1$, $i << e^{-2} << 1$. This situation may be typical for e.g. a population of protoplanetary cores in the end of the oligarchic phase of planet formation. In this regime of orbital parameters planetesimal scattering has an anisotropic character and strongly differs from scattering in thick ($i ~ e$) disks. We derive analytical expressions for the planetesimal scattering coefficients and compare them with numerical calculations. We find significant discrepancies in the inclination scattering coefficients obtained by the two approaches and ascribe this difference to the effects not accounted for in the analytical calculation: multiple scattering events (temporary captures, which may be relevant for the production of distant planetary satellites outside the Hill sphere) and distant interaction of planetesimals prior to their close encounter. Our calculations show that the inclination of a thin, dispersion-dominated planetesimal disk grows exponentially on a very short time scale implying that (1) such disks must be very short-lived and (2) planetesimal accretion in this dynamical phase is insignificant. Our results are also applicable to the dynamics of shear-dominated disks switching to the dispersion-dominated regime.

We explore the properties of cold gravitoturbulent accretion disks - non-fragmenting disks hovering on the verge of gravitational instability - using a realistic prescription for the effective viscosity caused by gravitational torques. This prescription is based on a direct relationship between the angular momentum transport in a thin accretion disk and the disk cooling in a steady state. Assuming that opacity is dominated by dust we are able to self-consistently derive disk properties for a given $\dot M$ assuming marginal gravitational stability. We also allow external irradiation of the disk and account for a non-zero background viscosity which can be due to the MRI. Spatial transitions between different co-existing disk states (e.g. between irradiated and self-luminous or between gravitoturbulent and viscous) are described and the location of the boundary at which disk must fragment is determined in a variety of situations. We demonstrate in particular that at low enough $\dot M$ external irradiation stabilizes gravitoturbulent disk against fragmentation all the way to infinity thus providing means of steady mass transport to the central object. Implications of our results for the possibility of planet formation by gravitational instability in protoplanetary disks and star formation in the Galactic Center and for the problem of feeding supermassive black holes in galactic nuclei are discussed.

Duration and period of transits in extrasolar planetary systems can exhibit long-term variations for a variety of reasons. Here we investigate how systemic proper motion, which steadily re-orients planetary orbit with respect to our line of sight, affects the timing of transits. We find that in a typical system with a period of several days proper motion at the level of 100 mas/yr makes transit duration vary at a rate ~10-100 ms/yr. In some isolated systems this variation is at the measurable level (can be as high as 0.6 s/yr for GJ436) and may exceed all other transit timing contributions (due to the general relativity, stellar quadrupole, etc.). In addition, proper motion causes evolution of the observed orbital period via the Shklovskii effect at a rate $\gtrsim 10$ $\mu$s/yr for the nearby transiting systems (0.26 ms/yr in GJ436), which in some cases exceeds all other contributions to $\dot P$. Earth's motion around the Sun gives rise to additional periodic timing signal (even for systems with zero intrinsic proper motion) allowing a full determination of the spatial orientation of the planetary orbit. Unlike most other timing effects the proper motion signatures persist even in systems with zero eccentricity and get stronger as the planetary period increases. They should be the dominant cause of transit timing variations in isolated wide separation (periods of months) systems that will be sought by Kepler.

Young stars are expected to gain most of their mass by accretion from a disk that forms around them as a result of angular momentum conservation in the collapsing protostellar cloud. Accretion initially proceeds at high rates of 10^{-6}-10^{-5} M_Sun/yr resulting in strong irradiation of the stellar surface by the hot inner portion of the disk and leading to the suppression of the intrinsic stellar luminosity. Here we investigate how this luminosity suppression affects evolution of the protostellar properties. Using simple model based on the energy balance of accreting star we demonstrate that disk irradiation causes only a slight increase of the protostellar radius, at the level of several per cent. Such a weak effect is explained by a minor role played by the intrinsic stellar luminosity (at the time when it is significantly altered by irradiation) in the protostellar energy budget compared to the stellar deuterium burning luminosity and the inflow of the gravitational potential energy brought in by the freshly accreted material. Our results justify the neglect of irradiation effects in previous studies of the protostellar growth via disk accretion. Evolution of some other actively accreting objects such as young brown dwarfs and planets should also be only weakly sensitive to the effects of disk irradiation.

In the initial formation stages young stars must acquire a significant fraction of their mass by accretion from a circumstellar disk that forms in the center of a collapsing protostellar cloud. Throughout this period mass accretion rates through the disk can reach 10^{-6}-10^{-5} M_Sun/yr leading to substantial energy release in the vicinity of stellar surface. We study the impact of irradiation of the stellar surface produced by the hot inner disk on properties of accreting fully convective low-mass stars, and also look at objects such as young brown dwarfs and giant planets. At high accretion rates irradiation raises the surface temperature of the equatorial region above the photospheric temperature T_0 that a star would have in the absence of accretion. The high-latitude (polar) parts of the stellar surface, where disk irradiation is weak, preserve their temperature at the level of T_0. In strongly irradiated regions an almost isothermal outer radiative zone forms on top of the fully convective interior, leading to the suppression of the local internal cooling flux derived from stellar contraction (similar suppression occurs in irradiated ``hot Jupiters''). Properties of this radiative zone likely determine the amount of thermal energy that gets advected into the convective interior of the star. Total intrinsic luminosity integrated over the whole stellar surface is reduced compared to the non-accreting case, by up to a factor of several in some systems (young brown dwarfs, stars in quasar disks, forming giants planets), potentially leading to the retardation of stellar contraction. Stars and brown dwarfs irradiated by their disks tend to lose energy predominantly through their cool polar regions while young giant planets accreting through the disk cool through their whole surface.

Precision timing is the key ingredient of ongoing pulsar-based gravitational wave searches and tests of general relativity using binary pulsars. The conventional approach to timing explicitly assumes that the radio emitting region is located at the center of the pulsar, while polarimetric observations suggest that radio emission is in fact produced at altitudes ranging from tens to thousands of kilometers above the neutron star surface. Here we present a calculation of the effects of finite emission height on the timing of binary pulsars using a simple model for the emitting region geometry. Finite height of emission changes the propagation path of radio photons through the binary and gives rise to a large spin velocity of the emission region co-rotating with the neutron star. Under favorable conditions these two effects introduce corrections to the conventional time delays at the microsecond level (for a millisecond pulsar in a double neutron star binary with a period of several hours and assuming the emission height of 100 km). Exploiting the dependence of the emission height on frequency (radius-to-frequency mapping) and using multi-frequency observations one should be able to detect these timing corrections even though they are formally degenerate with conventional time delays. Although even in the most accurately timed systems the magnitude of the finite emission height effects is currently somewhat below timing precision, longer-term observations and future facilities like SKA will make measurement of these effects possible, providing an independent check of existing emission height estimates.

Gravitationally unstable disks can fragment and form bound objects provided that their cooling time is short. In protoplanetary disks radiative cooling is likely to be too slow to permit formation of planets by fragmentation within several tens of AU from the star. Recently, convection has been suggested as a faster means of heat loss from the disk but here we demonstrate that it is only marginally more efficient than the radiative cooling. The crucial factor is the rate at which energy can be radiated from the disk photosphere, which is robustly limited from above in the convective case by the adiabatic temperature gradient (given a certain midplane temperature). Thus, although vigorous convection is definitely possible in disks, the inefficiency of radiative loss from the photosphere may create a bottleneck limiting the ability of the disk to form self-gravitating objects. Based on this argument we derive a set of analytical constraints which diagnose the susceptibility of an unstable disk to fragmentation and show that the formation of giant planets by fragmentation of protoplanetary disks is unlikely to occur at distances of tens of AU. At the same time these constraints do not preclude the possibility of fragmentation and star formation in accretion disks around supermassive black holes.

The measurement of the Shapiro time delay in binary pulsar systems with highly inclined orbit can be affected both by the motion of the pulsar's companion because of the finite time it takes a photon to cross the binary, and by the gravitational light bending if the orbit is sufficiently edge-on relative to the line of sight. Here we calculate the effect of retardation due to the companion's motion on various time delays in pulsar binaries, including the Shaipro delay, the geometric lensing delay, and the lens-induced delays associated with the pulsar rotation. Our results can be applied to systems so highly inclined that near conjunction gravitational lensing of the pulsar radiation by the companion becomes important (the recently discovered double pulsar system J0737-3039 may exemplify such a system). To the leading order, the effect of retardation is to shift all the delay curves backward in time around the orbit conjunction, without affecting the shape and amplitude of the curves. The time shift is of order the photon orbit crossing time, and ranges from a second to a few minutes for the observed binary pulsar systems. In the double pulsar system J0737-3039, the motion of the companion may also affect the interpretation of the recent correlated interstellar scintillation measurements. Finally, we show that lensing sets an upper limit on the magnitude of the frame-dragging time delay caused by the companion's spin, and makes this delay unobservable in stellar-mass binary pulsar systems.

We study the effect of pulsar rotation on timing of binary pulsars, with particular emphasis on the double pulsar system J0737-3039. Special relativistic aberration due to the orbital motion of pulsar changes both the longitude and colatitude of the emission direction with respect to the pulsar spin axis. The former gives rise to a shift of the arrival time of the pulse centroid (this is the conventional ``longitudinal'' aberration delay), the latter results in a distortion (contraction or dilation) of the pulse profile on the orbital timescale. The amplitude of pulse distortion depends inversely on the variation of polarization position angle across the pulse. For small angle between the pulsar magnetic and spin axes, as inferred for PSR J0737-3039A from polarimetric observations, the pulse distortion is significant (\sim 1%) and the associated ``latitudinal'' aberration delay is much larger than the longitudinal one. By monitoring the arrival time of separate pulse components as a function of pulsar orbital phase, the latitudinal aberration delay may be easily measured with the current timing precision, thus constraining the spin geometry of the system. The latitudinal delay can also be detected by monitoring system's orbital parameters on the geodetic precession timescale. Because of the near edge-on orbital orientation of the PSR J0737-3039 system, general relativistic bending of pulsar A's radio beam near its superior conjunction also introduces spin-dependent time delays of similar order of magnitude as the aberration delays. In addition, light bending splits the pulse profile into two variable components, corresponding to two gravitationally lensed images of the source. Detection of lensing effects is challenging, but may be possible with existing technology.

In the high density environments of circumstellar disks dust grains are expected to grow to large sizes by coagulation. Somewhat unexpectedly, recent near-IR observations of PAH features from disks around Herbig Ae/Be stars demonstrate that substantial amount of dust mass in these disks (up to several tens of per cent of the total carbon content) can be locked up in particles with sizes ranging from several to tens of nanometers. We investigate the possibility of detecting the electric dipole emission produced by these nanoparticles as they spin at thermal rates (tens of GHz) in cold gas. We show that such emission peaks in the microwave range and dominates over the thermal disk emission at \nu < 50 GHz typically by a factor of several if > 5 % of the total carbon abundance is locked up in nanoparticles. We test the sensitivity of this prediction to various stellar and disk parameters and show that if the potential contamination of the spinning dust component by the free-free and/or synchrotron emission can be removed, then the best chances of detecting this emission would be in disks with small opacity, having SEDs with steep sub-mm slopes (which minimizes thermal disk emission at GHz frequencies). Detection of the spinning dust emission would provide important evidence for the existence, properties, and origin of the population of small dust particles in protoplanetary disks, with possible ramifications for planet formation.

Structure and spectral energy distributions (SEDs) of externally irradiated circumstellar disks are often computed on the basis of the two-temperature model of Chiang & Goldreich. We refine these calculations by using a more realistic temperature profile which is continuous at all optical depths and thus goes beyond the two-temperature model. It is based on the approximate solution of the radiation transfer in the disk obtained from the frequency-integrated moment equations in the Eddington approximation. We come up with a simple procedure (``constant g_z approximation'') for treating the vertical structure of the disk in regions where its optical depth to stellar radiation is high. This allows us to obtain expressions for the vertical profiles of density and pressure at every point in the disk and to determine the shape of its surface. Armed with these analytical results we calculate the full radial structure of the disk and demonstrate that it favorably agrees with the results of direct numerical calculations. We also describe a simple and efficient way of the SED calculation based on our adopted temperature profile. Resulting spectra provide very good match (especially at short wavelengths) to the results of more detailed (but also more time-consuming) SED calculations solving the full frequency- and angle-dependent radiation transfer within the disk.

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

We investigate the effects of gravitational lensing in the binary pulsar system J0737-3039. Current measurement of the orbital inclination allows the millisecond pulsar (A) to pass very close (at R_{min}=4000 km) in projection to the companion pulsar (B), with R_{min} comparable to the Einstein radius (2600 km). For this separation at the conjunction, lensing causes small (about 10%) magnification of the pulsar A signal on a timescale of several seconds, and displaces the pulsar image on the sky plane by about 1200 km. More importantly, lensing introduces a correction (of several microsec) to the conventional Shapiro delay formula used in pulsar timing analysis, and gives rise to a geometric time delay together with the delays associated with the pulsar spin period. These lensing effects can influence the determination of the system parameters by both timing and scintillation studies. Given the current uncertainty in the orbital inclination, more extreme manifestations of lensing (e.g. magnification by a factor of several) are possible. We compare our predictions with the existing observations and discuss the possibility of detecting gravitational lensing signatures in the system. The anomalously high point in A's lightcurve close to superior conjunction might be caused by gravitational lensing.

Gravitational instability has been invoked as a possible mechanism of giant planet formation in protoplanetary disks. Here we critically revise its viability by noting that for the direct production of giant planets it is not enough for protoplanetary disks to be gravitationally unstable. They must also be able to cool efficiently (on a timescale comparable to the local disk orbital period) to allow the formation of bound clumps by fragmentation. Combination of dynamical and thermal constraints puts very stringent lower limits on the surface density and temperature of disks capable of fragmenting into self-gravitating objects: for the gravitational instability to form giant planets at 10 AU gas temperature at this location must exceed 10^3 K for a minimum disk mass of 0.7 M_Sun and minimum disk luminosity of 40 L_Sun. Although these requirements relax in more distant parts of the disk, masses of bound objects formed as a result of instability even at 100 AU are too large (\sim 10 M_Jupiter) to explain characteristics of known extrasolar giant planets. Such protoplanetary disks (and planets formed in them) have very unusual observational properties and this severely constrains the possibility of giant planet formation by direct gravitational instability.

We explore the dynamics of small planetesimals coexisting with massive protoplanetary cores in a gaseous nebula. Gas drag strongly affects the motion of small bodies leading to the decay of their eccentricities and inclinations, which are excited by the gravity of protoplanetary cores. Drag acting on larger ($\gtrsim 1$ km), high velocity planetesimals causes a mere reduction of their average random velocity. By contrast, drag qualitatively changes the dynamics of smaller ($\lesssim 0.1-1$ km), low velocity objects: (1) small planetesimals sediment towards the midplane of the nebula forming vertically thin subdisk; (2) their random velocities rapidly decay between successive passages of the cores and, as a result, encounters with cores typically occur at the minimum relative velocity allowed by the shear in the disk. This leads to a drastic increase in the accretion rate of small planetesimals by the protoplanetary cores, allowing cores to grow faster than expected in the simple oligarchic picture, provided that the population of small planetesimals contains more than roughly 1% of the solid mass in the nebula. Fragmentation of larger planetesimals ($\gtrsim 1$ km) in energetic collisions triggered by the gravitational scattering by cores can easily channel this amount of material into small bodies on reasonable timescales ($< 1$ Myr in the outer Solar System), providing a means for the rapid growth (within several Myr at 30 AU) of rather massive protoplanetary cores. Effects of inelastic collisions between planetesimals and presence of multiple protoplanetary cores are discussed.

We study quasi-static atmospheres of accreting protoplanetary cores for different opacity behaviors and realistic planetesimal accretion rates in various parts of protoplanetary nebula. Atmospheres segregate into those having outer convective zone which smoothly merges with the nebular gas, and those having almost isothermal outer radiative region decoupling atmospheric interior from the nebula. Specific type of atmosphere depends only on the relations between the Bondi radius of the core, photon mean free path in the nebular gas, and the luminosity radius (roughly the size of the sphere which can radiate luminosity of the core at effective temperature equal to the nebular temperature). Cores in the inner parts of protoplanetary disk (within roughly 0.3 AU from the Sun) have large luminosity radii resulting in the atmospheres of the first type, while cores in the giant planet region (beyond several AU) have small luminosity radii and always accumulate massive atmospheres of the second type. Critical core mass for nucleated instability is found to vary as a function of distance from the Sun. It is 5-20 M_Earth at 0.1-1 AU which is too large to permit the formation of ``hot Jupiters'' by nucleated instability near the cores that have grown in situ. In the region of giant planets critical mass is 20-60 M_Earth (for opacity 0.1 cm^2/g) if planetesimal accretion was fast enough for protoplanetary cores to form prior to the nebular gas dissipation. This might indicate that giant planets in the Solar System have gained their atmospheres by nucleated instability only after their cores have accumulated most of the mass in solids during the epoch of oligarchic growth, subsequent to which planetesimal accretion slowed down and cores became supercritical.

Formation of terrestrial planets by agglomeration of planetesimals in protoplanetary disks sensitively depends on the velocity evolution of planetesimals. We describe a novel semi-analytical approach to the treatment of planetesimal dynamics incorporating the gravitational scattering by massive protoplanetary bodies. Using this method we confirm that planets grow very slowly in the outer Solar System if gravitational scattering is the only process determining planetesimal velocities, making it hard for giant planets to acquire their massive gaseous envelopes within less than 10 Myr. We put forward several possibilities for alleviating this problem.

The current picture of terrestrial planet formation relies heavily on our understanding of the dynamical evolution of planetesimals -- asteroid-like bodies thought to be planetary building blocks. In this study we investigate the growth of eccentricities and inclinations of planetesimals in spatially homogeneous protoplanetary disks using methods of kinetic theory. We explore disks with a realistic mass spectrum of planetesimals evolving in time, similar to that obtained in self-consistent simulations of planetesimal coagulation. We calculate the behavior of planetesimal random velocities as a function of the planetesimal mass spectrum both analytically and numerically; results obtained by the two approaches agree quite well. Scaling of random velocity with mass can always be represented as a combination of power laws corresponding to different velocity regimes (shear- or dispersion-dominated) of planetesimal gravitational interactions. For different mass spectra we calculate analytically the exponents and time dependent normalizations of these power laws, as well as the positions of the transition regions between different regimes. It is shown that random energy equipartition between different planetesimals can only be achieved in disks with very steep mass distributions (differential surface number density of planetesimals falling off steeper than m^{-4}), or in the runaway tails. In systems with shallow mass spectra (shallower than m^{-3}) random velocities of small planetesimals turn out to be independent of their masses. We also discuss the damping effects of inelastic collisions between planetesimals and of gas drag, and their importance in modifying planetesimal random velocities.

We consider the growth of a protoplanetary embryo embedded in a planetesimal disk. We take into account the dynamical evolution of the disk caused by (1) planetesimal-planetesimal interactions, which increase random motions and smooth gradients in the disk, and (2) gravitational scattering of planetesimals by the embryo, which tends to heat up the disk locally and repels planetesimals away. The embryo's growth is self-consistently coupled to the planetesimal disk dynamics. We demonstrate that details of the evolution depend on only two dimensionless parameters incorporating all the physical characteristics of the problem: the ratio of the physical radius to the Hill radius of any solid body in the disk and the number of planetesimals inside the annulus of the disk with width equal to the planetesimal Hill radius. The results of exploration in the framework of our model of several situations typical for protosolar nebula can be summarized as follows: initially, the planetesimal disk dynamics is not affected by the presence of the embryo and the growth of the embryo's mass proceeds very rapidly in the runaway regime. Later on, when the embryo starts being dynamically important, its accretion slows down similar to the ``oligarchic'' growth picture. The scenario of orderly growth suggested by Safronov (1972) is never realized in our calculations; scenario of runaway growth suggested by Wetherill & Stewart (1989) is only realized for a limited range in mass. Slow character of the planetesimal accretion on the oligarchic stage of the embryo's accumulation leads to a considerable increase of the protoplanetary formation timescale compared to that following from a simple runaway accretion picture valid in the homogeneous planetesimal disks.

We investigate the process of an inhomogeneous planetesimal disk evolution caused by the planetesimal-planetesimal gravitational scattering. We develop a rather general approach based on the kinetic theory which self-consistently describes the evolution in time and space of both the disk surface density and its kinematic properties -- dispersions of eccentricity and inclination. The gravitational scattering of planetesimals is assumed to be in the dispersion-dominated regime which considerably simplifies analytical treatment. The resultant equations are of advection-diffusion type. Distance dependent scattering coefficients entering these equations are calculated analytically under the assumption of two-body scattering in the leading order in Coulomb logarithm. They are essentially nonlocal in nature. Our approach allows one to explore the dynamics of nonuniform planetesimal disks with arbitrary mass and random velocity distributions. It can also naturally include other physical mechanisms which are important for the evolution of such disks -- gas drag, migration, and so on.

The process of gravitational scattering of planetesimals by a massive protoplanetary embryo is explored theoretically. We propose a method to describe the evolution of the disk surface density, eccentricity, and inclination caused by the embryo-planetesimal interaction. It relies on the analytical treatment of the scattering in two extreme regimes of the planetesimal epicyclic velocities: shear-dominated (dynamically ``cold'') and dispersion-dominated (dynamically ``hot''). In the former, planetesimal scattering can be treated as a deterministic process. In the latter, scattering is mostly weak because of the large relative velocities of interacting bodies. This allows one to use the Fokker-Planck approximation and the two-body approximation to explore the disk evolution. We compare the results obtained by this method with the outcomes of the direct numerical integrations of planetesimal orbits and they agree quite well. In the intermediate velocity regime an approximate treatment of the disk evolution is proposed based on interpolation between the two extreme regimes. We also calculate the rate of embryo's mass growth in an inhomogeneous planetesimal disk and demonstrate that it is in agreement with both the simulations and earlier calculations. Finally we discuss the question of the direction of the embryo-planetesimal interaction in the dispersion-dominated regime and demonstrate that it is repulsive. This means that the embryo always forms a gap in the disk around it, which is in contrast with the results of other authors. The machinery developed here will be applied to realistic protoplanetary systems in future papers.

Gap formation in a gas disk triggered by disk-planet tidal interaction is considered. Density waves launched by the planet are assumed to be damped as a result of their nonlinear evolution leading to shock formation and its subsequent dissipation. As a consequence wave angular momentum is transferred to the disk,leading to evolution of its surface density. Planetary migration is an important ingredient of the theory; effects of the planet-induced surface density perturbations on the migration speed are considered. A gap is assumed to form when a stationary solution for the surface density profile is no longer possible in the frame of reference migrating with the planet. An analytical limit on the planetary mass necessary to open a gap in an inviscid disk is derived. The critical mass turns out to be smaller than mass M_1 for which planetary Hill's radius equals disk scaleheight by a factor of at least Q^{5/7} (Q is the Toomre stability parameter) depending on the strength of the migration feedback. In viscous disks the critical planetary mass could vary from about 0.1M_1 to M_1, depending on the disk viscosity. This implies that a gap could be formed by a planet with mass 1-10 times bigger than the Earth mass depending on the disk aspect ratio, viscosity, and planet's location in the nebula.

The propagation and evolution of planet-generated density waves in protoplanetary disks is considered. The evolution of waves, leading to the shock formation and wake dissipation, is followed in the weakly nonlinear regime. The local approach of Goodman & Rafikov (2001) is extended to include the effects of surface density and temperature variations in the disk as well as the disk cylindrical geometry and nonuniform shear. Wave damping due to shocks is demonstrated to be a nonlocal process spanning a significant fraction of the disk. Torques induced by the planet could be significant drivers of disk evolution on timescales of the order 1-10 Myr even in the absence of strong background viscosity. A global prescription for angular momentum deposition is developed which could be incorporated into the study of gap formation in a gaseous disk around the planet.

We discuss measurements of the properties of about 10,000 asteroids detected in 500 deg2 of sky in the Sloan Digital Sky Survey (SDSS) commissioning data. The moving objects are detected in the magnitude range 14 < r < 21.5, with a baseline of 5 minutes. Extensive tests show that the sample is at least 98% complete, with the contamination rate of less than 3%. We find that the size distribution of asteroids resembles a broken power-law, independent of the heliocentric distance: D^{-2.3} for 0.4 km < D < 5 km, and D^{-4} for 5 km < D < 40 km. As a consequence of this break, the number of asteroids with r < 21.5 is ten times smaller than predicted by extrapolating the power-law relation observed for brighter asteroids (r < 18). The observed counts imply that there are about 530,000 objects with D>1 km in the asteroid belt, or about four times less than previous estimates. The distribution of main belt asteroids in the 4-dimensional SDSS color space is bimodal, and the two groups can be associated with S (rocky) and C (carbonaceous) asteroids. A strong bimodality is also seen in the heliocentric distribution of asteroids and suggests the existence of two distinct belts: the inner rocky belt, about 1 AU wide (FWHM) and centered at R~2.8 AU, and the outer carbonaceous belt, about 0.5 AU wide and centered at R~3.2 AU. The colors of Hungarias, Mars crossers, and near-Earth objects are more similar to the C-type than to S-type asteroids, suggesting that they originate in the outer belt. (abridged).

The process of gap formation by a growing planetary embryo embedded in a planetesimal disk is considered. It is shown that there exists a single parameter characterizing this process, which represents the competition between the gravitational influence of the embryo and planetesimal-planetesimal scattering. For realistic assumptions about the properties of the planetesimal disk and the planetary embryo, a gap is opened long before the embryo can accrete all the bodies within its region of influence. The implication of this result is that the embryo stops growing and, thus, large bodies formed during the coagulation stage should be less massive than is usually assumed. For conditions expected at 1 AU in the solar protoplanetary disk, gap formation is expected to occur around bodies of mass < 10^24 g. The effect of protoplanetary radial migration is also discussed.

It has been proposed that the Galaxy might contain a population of cold clouds in numbers sufficient to account for a substantial fraction of the total mass of the Galaxy. These clouds would have masses of the order of 10^{-3} Solar mass and sizes of the order of 10 AU. We consider here the lensing effects of such clouds on the light from background stars. A semianalytical formalism for calculation of the magnification event rate produced by such gaseous lensing is developed, taking into account the spatial distribution of the dark matter in the Galaxy, the velocity distribution of the lensing clouds and source stars, and motion of the observer. Event rates are calculated for the case of gaseous lensing of stars in the Large Magellanic Cloud and results are directly compared with the results of the search for gravitational microlensing events undertaken by the MACHO collaboration. The MACHO experiment strongly constrains the properties of the proposed molecular clouds, but does not completely rule them out. Future monitoring programs will either detect or more strongly constrain this proposed population.

We revisit the idea that density-wave wakes of planets drive accretion in protostellar disks. The effects of many small planets can be represented as a viscosity if the wakes damp locally, but the viscosity is proportional to the damping length. Damping occurs mainly by shocks even for earth-mass planets. The excitation of the wake follows from standard linear theory including the torque cutoff. We use this as input to an approximate but quantitative nonlinear theory based on Burger's equation for the subsequent propagation and shock. Shock damping is indeed local but weakly so. If all metals in a minimum-mass solar nebula are invested in planets of a few earth masses each, dimensionless viscosities [alpha] of order dex(-4) to dex(-3) result. We compare this with observational constraints. Such small planets would have escaped detection in radial-velocity surveys and could be ubiquitous. If so, then the similarity of the observed lifetime of T Tauri disks to the theoretical timescale for assembling a rocky planet may be fate rather than coincidence.

Purely gravitational perturbations are considered in a thin rotating disk composed of several gas and stellar components. The dispersion relation for the axisymmetric density waves propagating through the disk is found and the criterion for the local axisymmetric stability of the whole system is formulated. In the appropriate limit of two-component gas we confirm the findings of Jog & Solomon (1984) and extend consideration to the case when one component is collisionless. Gravitational stability of the Galactic disk in the Solar neighborhood based on the multicomponent instability condition is explored using recent measurements of the stellar composition and kinematics in the local Galactic disk obtained by Hipparcos satellite.

The possibility of the magnetic flux expulsion from the Galaxy in the superbubble (SB) explosions, important for the Alpha-Omega dynamo, is considered. Special emphasis is put on the investigation of the downsliding of the matter from the top of the shell formed by the SB explosion which is able to influence the kinematics of the shell. It is shown that either Galactic gravity or the development of the Rayleigh-Taylor instabilities in the shell limit the SB expansion, thus, making impossible magnetic flux expulsion. The effect of the cosmic rays in the shell on the sliding is considered and it is shown that it is negligible compared to Galactic gravity. Thus, the question of possible mechanism of flux expulsion in the Alpha-Omega dynamo remains open.

The two-fluid effects on the radial outflow of relativistic electron-positron plasma are considered. It is shown that for large enough Michel magnetization parameter (1969) and multiplication parameter, the one-fluid MHD approximation remains correct in the whole region where |E| < |B|. In the case in which the longitudinal electric current is smaller than the Goldreich-Julian one, the acceleration of particles is determined near the light surface where |E| = |B|. It is shown that, as in the previously considered cylindrical geometry (Beskin, Gurevich & Istomin 1983), almost all electromagnetic energy is transformed into the energy of particles in the narrow boundary layer. The relative thickness of this layer is shown to be of the order of the inverse of multiplication parameter.