Massive planet migration: Theoretical predictions and comparison with observations

We quantify the utility of large radial velocity surveys for constraining theoretical models of Type II migration and protoplanetary disk physics. We describe a theoretical model for the expected radial distribution of extrasolar planets that combines an analytic description of migration with an empirically calibrated disk model. The disk model includes viscous evolution and mass loss via photoevaporation. Comparing the predicted distribution to a uniformly selected subsample of planets from the Lick / Keck / AAT planet search programs, we find that a simple model in which planets form in the outer disk at a uniform rate, migrate inward according to a standard Type II prescription, and become stranded when the gas disk is dispersed, is consistent with the radial distribution of planets for orbital radii 0.1 AU < a < 2.5 AU and planet masses greater than 1.65 Jupiter masses. Some variant models are disfavored by existing data, but the significance is limited (~95%) due to the small sample of planets suitable for statistical analysis. We show that the favored model predicts that the planetary mass function should be almost independent of orbital radius at distances where migration dominates the massive planet population. We also study how the radial distribution of planets depends upon the adopted disk model. We find that the distribution can constrain not only changes in the power-law index of the disk viscosity, but also sharp jumps in the efficiency of angular momentum transport that might occur at small radii.Comment: ApJ, in press. References updated to match published versio

Scale-free equilibria of self-gravitating gaseous disks with flat rotation curves

We introduce exact analytical solutions of the steady-state hydrodynamic equations of scale-free, self-gravitating gaseous disks with flat rotation curves. We express the velocity field in terms of a stream function and obtain a third-order ordinary differential equation (ODE) for the angular part of the stream function. We present the closed-form solutions of the obtained ODE and construct hydrodynamical counterparts of the power-law and elliptic disks, for which self-consistent stellar dynamical models are known. We show that the kinematics of the Large Magellanic Cloud can well be explained by our findings for scale-free elliptic disks.Comment: AAS preprint format, 21 pages, 8 figures, accepted for publication in The Astrophysical Journa

Evolution of Giant Planets in Eccentric Disks

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

Constraining the Mass Profiles of Stellar Systems: Schwarzschild Modeling of Discrete Velocity Datasets

(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations that can only address continuous data in the form of (the moments of) velocity distributions, thus avoiding potentially important losses of information due to data binning. Furthermore, the code can handle any combination of available velocity components, i.e., only line-of-sight velocities, only proper motions, or a combination of both. It can also handle a combination of discrete and continuous data. The code finds the distribution function (DF, a function of the three integrals of motion E, Lz, and I3) that best reproduces the available kinematic and photometric observations in a given axisymmetric gravitational potential. The fully numerical approach ensures considerable freedom on the form of the DF f(E,Lz,I3). This allows a very general modeling of the orbital structure, thus avoiding restrictive assumptions about the degree of (an)isotropy of the orbits. We describe the implementation of the discrete code and present a series of tests of its performance based on the modeling of simulated datasets generated from a known DF. We find that the discrete Schwarzschild code recovers the original orbital structure, M/L ratios, and inclination of the input datasets to satisfactory accuracy, as quantified by various statistics. The code will be valuable, e.g., for modeling stellar motions in Galactic globular clusters, and those of individual stars, planetary nebulae, or globular clusters in nearby galaxies. This can shed new light on the total mass distributions of these systems, with central black holes and dark matter halos being of particular interest.Comment: ApJ, in press; 51 pages, 11 figures; manuscript revised following comments by refere

Measurement of an AGN Central Mass on Centiparsec Scales: Results of Long-Term Optical Monitoring of Arp 102B

The optical spectrum of the broad-line radio galaxy Arp 102B has been monitored for more than thirteen years to investigate the nature of the source of its broad, double-peaked hydrogen Balmer emission lines. The shape of the lines varied subtly; there was an interval during which the variation in the ratio of the fluxes of the two peaks appeared to be sinusoidal, with a period of 2.16 years and an amplitude of about 16% of the average value. The variable part of the broad H-alpha line is well fit by a model in which a region of excess emission (a quiescent ``hot spot'') within an accretion disk (fitted to the non-varying portion of the double-peaked line) completes at least two circular orbits and eventually fades. Fits to spectra from epochs when the hot spot is not present allow determination of the disk inclination, while fits for epochs when it is present provide a measurement of the radius of the hot spot's orbit. From these data and the period of variation, we find that the mass within the hot spot's orbit is 2.2 +0.2/-0.7 times 10^8 solar masses, within the range of previous estimates of masses of active galactic nuclei. Because this mass is determined at a relatively small distance (~1000 AU) from the central body, it is extremely difficult to explain without assuming that a supermassive black hole lies within Arp 102B. The lack of any systematic change in the velocity of the blue peak over time yields a lower limit on the combined mass of the two bodies in a binary black hole model like that of Gaskell (1983) of 10^10 solar masses.Comment: 29 pages, including 6 figures; to appear in the Astrophysical Journal 199

A Photometric Method for Quantifying Asymmetries in Disk Galaxies

A photometric method for quantifying deviations from axisymmetry in optical images of disk galaxies is applied to a sample of 32 face-on and nearly face-on spirals. The method involves comparing the relative fluxes contained within trapezoidal sectors arranged symmetrically about the galaxy center of light, excluding the bulge and/or barred regions. Such a method has several advantages over others, especially when quantifying asymmetry in flocculent galaxies. Specifically, the averaging of large regions improves the signal-to-noise in the measurements; the method is not strongly affected by the presence of spiral arms; and it identifies the kinds of asymmetry that are likely to be dynamically important. Application of this "method of sectors" to R-band images of 32 disk galaxies indicates that about 30% of spirals show deviations from axisymmetry at the 5-sigma level.Comment: 17 pages, 2 tables and 6 figures, uses psfig and AAS LaTex; to appear in A

Neutralino dark matter vs galaxy formation

Neutralino dark matter may be incompatible with current cold dark matter models with cuspy dark halos, because excessive synchrotron radiation may originate from neutralino annihilations close to the black hole at the galactic center.Comment: 6 pages, 3 figures, talk given at "Sources and detection of dark matter in the Universe", Marina del Rey, CA, February 23-25, 200

Numerical simulation of small perturbation on an accretion disk due to the collision of a star with the disk near the black hole

In this paper, perturbations of an accretion disk by a star orbiting around a black hole are studied. We report on a numerical experiment, which has been carried out by using a parallel-machine code originally developed by D\"{o}nmez (2004). An initially steady state accretion disk near a non-rotating (Schwarzschild) black hole interacts with a "star", modeled as an initially circular region of increased density. Part of the disk is affected by the interaction. In some cases, a gap develops and shock wave propagates through the disk. We follow the evolution for order of one dynamical period and we show how the non-axisymetric density perturbation further evolves and moves downwards where the material of the disk and the star become eventually accreted onto the central body. When the star perturbs the steady state accretion disk, the disk around the black hole is destroyed by the effect of perturbation. The perturbed accretion disk creates a shock wave during the evolution and it loses angular momentum when the gas hits on the shock waves. Colliding gas with the shock wave is the one of the basic mechanism of emitting the $X-$rays in the accretion disk. The series of supernovae occurring in the inner disk could entirely destroy the disk in that region which leaves a more massive black hole behind, at the center of galaxies.Comment: 20pages, 8 figures, accepted for publication in Astrophysics and Space Scienc

A Magnetic Alpha-Omega Dynamo in Active Galactic Nuclei Disks: I. The Hydrodynamics of Star-Disk Collisions and Keplerian Flow

A magnetic field dynamo in the inner regions of the accretion disk surrounding the supermassive black holes in AGNs may be the mechanism for the generation of magnetic fields in galaxies and in extragalactic space. We argue that the two coherent motions produced by 1) the Keplerian motion and 2) star-disk collisions, numerous in the inner region of AGN accretion disks, are both basic to the formation of a robust, coherent dynamo and consequently the generation of large scale magnetic fields. They are frequent enough to account for an integrated dynamo gain, e^{10^{9}} at 100 gravitational radii of a central black hole, many orders of magnitude greater than required to amplify any seed field no matter how small. The existence of extra-galactic, coherent, large scale magnetic fields whose energies greatly exceed all but massive black hole energies is recognized. In paper II (Pariev, Colgate, and Finn 2006) we argue that in order to produce a dynamo that can access the free energy of black hole formation and produce all the magnetic flux in a coherent fashion the existence of these two coherent motions in a conducting fluid is required. The differential winding of Keplerian motion is obvious, but the disk structure depends upon the model of "alpha", the transport coefficient of angular momentum chosen. The counter rotation of driven plumes in a rotating frame is less well known, but fortunately the magnetic effect is independent of the disk model. Both motions are discussed in this paper, paper I. The description of the two motions are preliminary to two theoretical derivations and one numerical simulation of the alpha-omega dynamo in paper II. (Abridged)Comment: 34 pages, 1 figure, accepted by Ap