-
Galactokinetics
Authors:
Chris Hamilton,
Shaunak Modak,
Scott Tremaine
Abstract:
Galactic disks lie at the heart of many of the most pressing astrophysical puzzles. There are sophisticated kinetic theories that describe some aspects of galaxy disk dynamics, but extracting quantitative predictions from those theories has proven very difficult, meaning they have shed little light on observations/simulations of galaxies. Here, we begin to address this issue by developing a tracta…
▽ More
Galactic disks lie at the heart of many of the most pressing astrophysical puzzles. There are sophisticated kinetic theories that describe some aspects of galaxy disk dynamics, but extracting quantitative predictions from those theories has proven very difficult, meaning they have shed little light on observations/simulations of galaxies. Here, we begin to address this issue by developing a tractable theory describing fluctuations and transport in thin galactic disks. Our main conceptual advance is to split potential fluctuations into asymptotic wavelength regimes relative to orbital guiding radius and epicyclic amplitude (similar to plasma gyrokinetics), and then to treat separately the dynamics in each regime. As an illustration, we apply our results to quasilinear theory, calculating the angular-momentum transport due to a transient spiral. At each stage we verify our formulae with numerical examples. Our approach should simplify many important calculations in galactic disk dynamics.
△ Less
Submitted 6 August, 2024;
originally announced August 2024.
-
Iron Snails: non-equilibrium dynamics and spiral abundance patterns
Authors:
Neige Frankel,
David W. Hogg,
Scott Tremaine,
Adrian Price-Whelan,
Jeff Shen
Abstract:
Galaxies are not in a dynamical steady state. They continually undergo perturbations, e.g., from infalling dwarf galaxies and dark-matter substructure. After a dynamical perturbation, stars phase mix towards a new steady state; in so doing they generally form spiral structures, such as spiral density waves in galaxy disks and the Gaia Snail observed in the vertical phase-space density in the solar…
▽ More
Galaxies are not in a dynamical steady state. They continually undergo perturbations, e.g., from infalling dwarf galaxies and dark-matter substructure. After a dynamical perturbation, stars phase mix towards a new steady state; in so doing they generally form spiral structures, such as spiral density waves in galaxy disks and the Gaia Snail observed in the vertical phase-space density in the solar neighborhood. Structures in phase-space density can be hard to measure accurately, because spatially varying selection effects imprint their own patterns on the density. However, stellar labels such as metallicity, or other element abundances, or stellar masses and ages, can be measured even in the face of complex or unknown spatial selection functions. We show that if the equilibrium galaxy has phase-space gradients in these labels, any perturbation that could raise a spiral wave in the phase-space density will raise a spiral wave in the distribution of labels as well. We work out the relationship between the spiral patterns in the density and in the labels. As an example, we analyze the Gaia Snail and show that its amplitude and dynamical age as derived from elemental abundances (mainly [Mg/Fe]) follow similar patterns to those derived from the phase-space density. Our best model dates the Snail's perturbation to about 400 Myr ago although we find significant variations with angular momentum in the best-fit age. Conceptually, the ideas presented here are related to Orbital Torus Imaging, chemical tagging, and other methods that use stellar labels to trace dynamics.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Studying Binary Formation under Dynamical Friction Using Hill's Problem
Authors:
Mark Dodici,
Scott Tremaine
Abstract:
Using the equations of motion from Hill's problem, with added accelerations for different forms of dynamical friction, we provide the (to-date) broadest scale-free study of friction-driven binary formation in gaseous disks and stellar clusters. We focus mainly on binary formation between stellar-mass black holes in active galactic nuclei (AGNi), considering both gas dynamical friction from AGN dis…
▽ More
Using the equations of motion from Hill's problem, with added accelerations for different forms of dynamical friction, we provide the (to-date) broadest scale-free study of friction-driven binary formation in gaseous disks and stellar clusters. We focus mainly on binary formation between stellar-mass black holes in active galactic nuclei (AGNi), considering both gas dynamical friction from AGN disks and stellar dynamical friction from the nuclear star cluster. We first find simple, dimensionless friction coefficients that approximate the effects of standard models for gas and stellar dynamical friction. We perform extensive simulations of Hill's problem under such friction, and we present a picture of binary formation through encounters between single stars on nearby orbits, as a function of friction parameter, eccentricity, and inclination. Notably, we find that the local binary formation rate is a linear function of the friction coefficient so long as the friction is weak. Due to the dimensionless nature of our model problem, our findings are generalizable to binary formation at all scales (e.g., intermediate-mass black holes in a star cluster, planetesimals in a gaseous disk).
△ Less
Submitted 11 April, 2024;
originally announced April 2024.
-
Testing MOND on small bodies in the remote solar system
Authors:
David Vokrouhlický,
David Nesvorný,
Scott Tremaine
Abstract:
Modified Newtonian dynamics (MOND), which postulates a breakdown of Newton's laws of gravity/dynamics below some critical acceleration threshold, can explain many otherwise puzzling observational phenomena on galactic scales. MOND competes with the hypothesis of dark matter, which successfully explains the cosmic microwave background and large-scale structure. Here we provide the first solar-syste…
▽ More
Modified Newtonian dynamics (MOND), which postulates a breakdown of Newton's laws of gravity/dynamics below some critical acceleration threshold, can explain many otherwise puzzling observational phenomena on galactic scales. MOND competes with the hypothesis of dark matter, which successfully explains the cosmic microwave background and large-scale structure. Here we provide the first solar-system test of MOND that probes the sub-critical acceleration regime. Using the Bekenstein-Milgrom AQUAL formulation, we simulate the evolution of myriads of test particles (planetesimals or comets) born in the trans-Neptunian region and scattered by the giant planets over the lifetime of the Sun to heliocentric distances of $10^2$-$10^5$ au. We include the effects of the Galactic tidal field and passing stars. While Newtonian simulations reproduce the distribution of binding energies of long-period and Oort-cloud comets detectable from Earth, MOND-based simulations do not. This conclusion is robust to plausible changes in the migration history of the planets, the migration history of the Sun, the MOND transition function, effects of the Sun's birth cluster, and the fading properties of long-period comets. For the most popular version of AQUAL, characterized by a gradual transition between the Newtonian and MOND regimes, our MOND-based simulations also fail to reproduce the orbital distribution of trans-Neptunian objects in the detached disk (perihelion > 38 au). Our results do not rule out some MOND theories more elaborate than AQUAL, in which non-Newtonian effects are screened on small spatial scales, at small masses, or in external gravitational fields comparable in strength to the critical acceleration.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Scattered Disk Dynamics: The Mapping Approach
Authors:
Sam Hadden,
Scott Tremaine
Abstract:
We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map shows excellent agreement with direct numerical integrations and elucidates how the dynamics depends on perturber mass and pericenter distance. We also use the…
▽ More
We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map shows excellent agreement with direct numerical integrations and elucidates how the dynamics depends on perturber mass and pericenter distance. We also use the map to explore the onset of chaos, statistical descriptions of chaotic transport, and sticking in mean-motion resonances. We discuss implications of our mapping model for the dynamical evolution of the solar system's scattered disk and other highly eccentric trans-Neptunian objects.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
The Hamiltonian for von Zeipel-Lidov-Kozai oscillations
Authors:
Scott Tremaine
Abstract:
The Hamiltonian used in classical analyses of von Zeipel-Lidov-Kozai or ZLK oscillations in hierarchical triple systems is based on the quadrupole potential from a distant body on a fixed orbit, averaged over the orbits of both the inner and the outer bodies ("double-averaging"). This approximation can be misleading, because the corresponding Hamiltonian conserves the component of angular momentum…
▽ More
The Hamiltonian used in classical analyses of von Zeipel-Lidov-Kozai or ZLK oscillations in hierarchical triple systems is based on the quadrupole potential from a distant body on a fixed orbit, averaged over the orbits of both the inner and the outer bodies ("double-averaging"). This approximation can be misleading, because the corresponding Hamiltonian conserves the component of angular momentum of the inner binary normal to the orbit of the outer binary, thereby restricting the volume of phase space that the system can access. This defect is usually remedied by including the effects of the octopole potential, or by allowing the outer orbit to respond to variations in the inner orbit. However, in a wide variety of astrophysical systems nonlinear perturbations are comparable to or greater than these effects. The long-term effects of nonlinear perturbations are described by an additional Hamiltonian, which we call Brown's Hamiltonian. At least three different forms of Brown's Hamiltonian are found in the literature; we show that all three are related by a gauge freedom, although one is much simpler than the others. We argue that investigations of ZLK oscillations in triple systems should include Brown's Hamiltonian.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Vertical motion in the Galactic disc: unwinding the Snail
Authors:
Neige Frankel,
Jo Bovy,
Scott Tremaine,
David W. Hogg
Abstract:
The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Sn…
▽ More
The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Snail. Here, we devise and fit a simple model in which the spiral starts as a lopsided perturbation from steady state, that winds up into the present-day morphology. The winding occurs because the vertical frequency decreases with vertical action. We use data from stars in Gaia EDR3 that have measured radial velocities, pruned by simple distance and photometric selection functions. We divide the data into boxels of dynamical invariants (radial action, angular momentum); our model fits the data well in many of the boxels. The model parameters have physical interpretations: one, $A$, is a perturbation amplitude, and one, $t$, is interpretable in the simplest models as the time since the event that caused the Snail. We find trends relating the strength and age to angular momentum: (i) the amplitude $A$ is small at low angular momentum ($<1\,600\mathrm{\,kpc\,km\,s}^{-1}$ or guiding-centre radius $< 7.3\,$kpc), and over a factor of three larger, with strong variations, in the outer disc; (ii) there is no single well-defined perturbation time, with $t$ varying between 0.2 and 0.6 Gyr. Residuals between the data and the model display systematic trends, implying that the data call for more complex models.
△ Less
Submitted 22 December, 2022;
originally announced December 2022.
-
The origin and fate of the Gaia phase-space snail
Authors:
Scott Tremaine,
Neige Frankel,
Jo Bovy
Abstract:
The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar nei…
▽ More
The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar neighbourhood. In this paper we investigate the alternative hypothesis that the snail is created by many small disturbances rather than one large one, that is, by Gaussian noise in the gravitational potential. Probably most of this noise is due to substructures in the dark-matter halo. We show that this hypothesis naturally reproduces most of the properties of the snail. In particular it predicts correctly, with no free parameters, that the apparent age of the snail will be $ \sim 0.5$ Gyr. An important ingredient of this model is that any snail-like feature in the solar neighbourhood, whatever its cause, is erased by scattering from giant molecular clouds or other small-scale structure on a time-scale $\lesssim 1$ Gyr.
△ Less
Submitted 22 December, 2022;
originally announced December 2022.
-
A Criterion for the Onset of Chaos in Compact, Eccentric Multiplanet Systems
Authors:
Daniel Tamayo,
Norman Murray,
Scott Tremaine,
Joshua Winn
Abstract:
We derive a semi-analytic criterion for the presence of chaos in compact, eccentric multiplanet systems. Beyond a minimum semimajor-axis separation, below which the dynamics are chaotic at all eccentricities, we show that (i) the onset of chaos is determined by the overlap of two-body mean motion resonances (MMRs), like it is in two-planet systems; (ii) secular evolution causes the MMR widths to e…
▽ More
We derive a semi-analytic criterion for the presence of chaos in compact, eccentric multiplanet systems. Beyond a minimum semimajor-axis separation, below which the dynamics are chaotic at all eccentricities, we show that (i) the onset of chaos is determined by the overlap of two-body mean motion resonances (MMRs), like it is in two-planet systems; (ii) secular evolution causes the MMR widths to expand and contract adiabatically, so that the chaotic boundary is established where MMRs overlap at their greatest width. For closely spaced two-planet systems, a near-symmetry strongly suppresses this secular modulation, explaining why the chaotic boundaries for two-planet systems are qualitatively different from cases with more than two planets. We use these results to derive an improved angular-momentum-deficit (AMD) stability criterion, i.e., the critical system AMD below which stability should be guaranteed. This introduces an additional factor to the expression from Laskar and Petit (2017) that is exponential in the interplanetary separations, which corrects the AMD threshold toward lower eccentricities by a factor of several for tightly packed configurations. We make routines for evaluating the chaotic boundary available to the community through the open-source SPOCK package.
△ Less
Submitted 23 August, 2021; v1 submitted 28 June, 2021;
originally announced June 2021.
-
Thermal equilibrium of an ideal gas in a free-floating box
Authors:
Scott Tremaine,
Bence Kocsis,
Abraham Loeb
Abstract:
The equilibrium and fluctuations of an ideal gas in a rigid container are studied by every student of statistical mechanics. Here we study the less well-known case when the box is floating freely; in particular we determine the fluctuations of the box in velocity and position due to interactions with the gas it contains. This system is a toy model for the fluctuations in velocity and position of a…
▽ More
The equilibrium and fluctuations of an ideal gas in a rigid container are studied by every student of statistical mechanics. Here we study the less well-known case when the box is floating freely; in particular we determine the fluctuations of the box in velocity and position due to interactions with the gas it contains. This system is a toy model for the fluctuations in velocity and position of a black hole surrounded by stars at the center of a galaxy. These fluctuations may be observable in nearby galaxies.
△ Less
Submitted 24 December, 2020;
originally announced December 2020.
-
Secular Dynamics around a Supermassive Black Hole via Multipole Expansion
Authors:
Jean-Baptiste Fouvry,
Walter Dehnen,
Scott Tremaine,
Ben Bar-Or
Abstract:
In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small ti…
▽ More
In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small timesteps. An alternative approach, pioneered by Gauss, is to use the secular approximation of smearing out the $N$ stars over their Keplerian orbits, using $K$ nodes along each orbit. In this study we propose three novel improvements to this method. First, we re-formulate the discretisation of the rates of change of the variables describing the orbital states to ensure that all conservation laws are exactly satisfied. Second, we replace the pairwise sum over nodes by a multipole expansion up to order $\ell_{\mathrm{max}}$, reducing the overall computational costs from $O(N^2K^2)$ to $O(NK\ell_{\mathrm{max}}^2)$. Finally, we show that the averaged dynamical system is equivalent to $2N$ interacting unit spin vectors and provide two time integrators: a second-order symplectic scheme and a fourth-order Lie-group Runge-Kutta method, both of which are straightforward to generalize to higher order. These new simulations recover the diffusion coefficients of stellar eccentricities obtained through analytical calculations of the secular dynamics.
△ Less
Submitted 1 March, 2022; v1 submitted 3 November, 2020;
originally announced November 2020.
-
Relaxation in a Fuzzy Dark Matter Halo. II. Self-consistent kinetic equations
Authors:
Ben Bar-Or,
Jean-Baptiste Fouvry,
Scott Tremaine
Abstract:
Fuzzy dark matter (FDM) is composed of ultra-light bosons having a de Broglie wavelength that is comparable to the size of the stellar component of galaxies at typical galactic velocities. FDM behaves like cold dark matter on large scales. However, on the scale of the de Broglie wavelength, an FDM halo exhibits density fluctuations that lead to relaxation, a process similar to the two-body relaxat…
▽ More
Fuzzy dark matter (FDM) is composed of ultra-light bosons having a de Broglie wavelength that is comparable to the size of the stellar component of galaxies at typical galactic velocities. FDM behaves like cold dark matter on large scales. However, on the scale of the de Broglie wavelength, an FDM halo exhibits density fluctuations that lead to relaxation, a process similar to the two-body relaxation that occurs in classical gravitational N-body systems and is described by the Fokker-Planck equation. We derive the FDM analog of that kinetic equation, and solve it to find the evolution of the velocity distribution in a spatially homogeneous FDM halo. We also determine the dielectric function and the dispersion relation for linear waves in an FDM halo.
△ Less
Submitted 20 October, 2020;
originally announced October 2020.
-
Peas in a Pod? Radius correlations in Kepler multi-planet systems
Authors:
Lena Murchikova,
Scott Tremaine
Abstract:
We address the claim of Weiss et al. (2018) that the radii of adjacent planets in Kepler multi-planet systems are correlated. We explore two simple toy models---in the first the radii of the planets are chosen at random from a single universal distribution, and in the second we postulate several types of system with distinct radius distributions. We show that an apparent correlation between the ra…
▽ More
We address the claim of Weiss et al. (2018) that the radii of adjacent planets in Kepler multi-planet systems are correlated. We explore two simple toy models---in the first the radii of the planets are chosen at random from a single universal distribution, and in the second we postulate several types of system with distinct radius distributions. We show that an apparent correlation between the radii of adjacent planets similar to the one reported by Weiss et al. (2018) can arise in both models. In addition the second model fits the radius and signal-to-noise distribution of the observed planets. We also comment on the validity of a commonly used correction that is used to estimate intrinsic planet occurrence rates, based on weighting planets by the inverse of their detectability.
△ Less
Submitted 14 August, 2020; v1 submitted 4 March, 2020;
originally announced March 2020.
-
Resonant capture in quadruple stellar systems
Authors:
Scott Tremaine
Abstract:
Some quadruple star systems in the hierarchical 2+2 configuration exhibit orbit-orbit resonances between the two compact binaries. We show that the most important resonances occur at period ratios of 1:1, 3:2 and 2:1. We describe the conditions required for capture and show that they can be satisfied at the 3:2 and 2:1 resonances in binaries that migrate significantly in semimajor axis after circu…
▽ More
Some quadruple star systems in the hierarchical 2+2 configuration exhibit orbit-orbit resonances between the two compact binaries. We show that the most important resonances occur at period ratios of 1:1, 3:2 and 2:1. We describe the conditions required for capture and show that they can be satisfied at the 3:2 and 2:1 resonances in binaries that migrate significantly in semimajor axis after circularization, probably through magnetic braking or gravitational radiation.
△ Less
Submitted 7 April, 2020; v1 submitted 2 March, 2020;
originally announced March 2020.
-
Order-disorder phase transition in black-hole star clusters -- III. A mono-energetic cluster
Authors:
Scott Tremaine
Abstract:
Supermassive black holes at the centres of galaxies are often surrounded by dense star clusters. For a wide range of cluster properties and orbital radii the resonant relaxation times in these clusters are much shorter than the Hubble time. Since resonant relaxation conserves semimajor axes, these clusters should be in the maximum-entropy state consistent with the given semimajor axis distribution…
▽ More
Supermassive black holes at the centres of galaxies are often surrounded by dense star clusters. For a wide range of cluster properties and orbital radii the resonant relaxation times in these clusters are much shorter than the Hubble time. Since resonant relaxation conserves semimajor axes, these clusters should be in the maximum-entropy state consistent with the given semimajor axis distribution. We determine these maximum-entropy equilibria in a simplified model in which all of the stars have the same semimajor axes. We find that the cluster exhibits a phase transition from a disordered, spherical, high-temperature equilibrium to an ordered low-temperature equilibrium in which the stellar orbits have a preferred orientation or line of apsides. Here `temperature' is a measure of the non-Keplerian or self-gravitational energy of the cluster; in the spherical state, temperature is a function of the rms eccentricity of the stars. We explore a simple two-parameter model of black-hole star clusters -- the two parameters are semimajor axis and black-hole mass --- and find that clusters are susceptible to the lopsided phase transition over a range of ~100 in semimajor axis, mostly for black-hole masses less than $10^{7.5}$ solar masses.
△ Less
Submitted 12 February, 2020;
originally announced February 2020.
-
Order-disorder phase transition in black-hole star clusters. II. A scale-free cluster
Authors:
Scott Tremaine
Abstract:
The supermassive black holes found at the centers of galaxies are often surrounded by dense star clusters. The ages of these clusters are generally longer than the resonant-relaxation time and shorter than the two-body relaxation time over a wide range of radii. We explore the thermodynamic equilibria of such clusters using a simple self-similar model. We find that the cluster exhibits a phase tra…
▽ More
The supermassive black holes found at the centers of galaxies are often surrounded by dense star clusters. The ages of these clusters are generally longer than the resonant-relaxation time and shorter than the two-body relaxation time over a wide range of radii. We explore the thermodynamic equilibria of such clusters using a simple self-similar model. We find that the cluster exhibits a phase transition between a high-temperature spherical equilibrium and a low-temperature equilibrium in which the stars are on high-eccentricity orbits with nearly the same orientation. In the absence of relativistic precession, the spherical equilibrium is metastable below the critical temperature and the phase transition is first-order. When relativistic effects are important, the spherical equilibrium is linearly unstable below the critical temperature and the phase transition is continuous. A similar phase transition has recently been found in a model cluster composed of stars with a single semimajor axis. The presence of the same phenomenon in two quite different cluster models suggests that lopsided equilibria may form naturally in a wide variety of black-hole star clusters.
△ Less
Submitted 17 December, 2019; v1 submitted 11 November, 2019;
originally announced November 2019.
-
An order-disorder phase transition in black-hole star clusters
Authors:
Jihad Touma,
Scott Tremaine,
Mher Kazandjian
Abstract:
The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that,…
▽ More
The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.
△ Less
Submitted 2 July, 2019;
originally announced July 2019.
-
Relaxation in a Fuzzy Dark Matter Halo
Authors:
Ben Bar-Or,
Jean-Baptiste Fouvry,
Scott Tremaine
Abstract:
Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $λ\sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales (…
▽ More
Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $λ\sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales ($\lesssim 10\,\mathrm{kpc}$). Because of its wave nature, FDM exhibits stochastic density fluctuations on the scale of the de Broglie wavelength that never damp. The gravitational field from these fluctuations scatters stars and black holes, causing their orbits to diffuse through phase space. We show that this relaxation process can be analyzed quantitatively with the same tools used to analyze classical two-body relaxation in an $N$-body system, and can be described by treating the FDM fluctuations as quasiparticles, with effective mass $\sim 10^7 M_\odot {(1\,\mathrm{kpc}/r)}^2{(10^{-22}\,\mathrm{eV}/m_b)}^3$ in a galaxy with a constant circular speed of $200\,\mathrm{kms}$. This novel relaxation mechanism may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred pc, and can heat and expand the central regions of galaxies. These processes can be used to constrain the mass of the light bosons that might comprise FDM.
△ Less
Submitted 25 January, 2019; v1 submitted 20 September, 2018;
originally announced September 2018.
-
Measuring the mass distribution in stellar systems
Authors:
Scott Tremaine
Abstract:
One of the fundamental tasks of dynamical astronomy is to infer the distribution of mass in a stellar system from a snapshot of the positions and velocities of its stars. The usual approach to this task (e.g., Schwarzschild's method) involves fitting parametrized forms of the gravitational potential and the phase-space distribution to the data. We review the practical and conceptual difficulties w…
▽ More
One of the fundamental tasks of dynamical astronomy is to infer the distribution of mass in a stellar system from a snapshot of the positions and velocities of its stars. The usual approach to this task (e.g., Schwarzschild's method) involves fitting parametrized forms of the gravitational potential and the phase-space distribution to the data. We review the practical and conceptual difficulties with this approach and describe a novel statistical method for determining the mass distribution that does not require determining the phase-space distribution of the stars. We show that this new estimator out-performs other distribution-free estimators for the harmonic and Kepler potentials.
△ Less
Submitted 23 December, 2017;
originally announced December 2017.
-
Producing Distant Planets by Mutual Scattering of Planetary Embryos
Authors:
Kedron Silsbee,
Scott Tremaine
Abstract:
It is likely that multiple bodies with masses between those of Mars and Earth ("planetary embryos") formed in the outer planetesimal disk of the solar system. Some of these were likely scattered by the giant planets into orbits with semi-major axes of hundreds of AU. Mutual torques between these embryos may lift the perihelia of some of them beyond the orbit of Neptune, where they are no longer pe…
▽ More
It is likely that multiple bodies with masses between those of Mars and Earth ("planetary embryos") formed in the outer planetesimal disk of the solar system. Some of these were likely scattered by the giant planets into orbits with semi-major axes of hundreds of AU. Mutual torques between these embryos may lift the perihelia of some of them beyond the orbit of Neptune, where they are no longer perturbed by the giant planets so their semi-major axes are frozen in place. We conduct N-body simulations of this process, and its effect on smaller planetesimals in the region of the giant planets and the Kuiper belt. We find that (i) there is a significant possibility that one sub-Earth mass embryo, or possibly more, is still present in the outer solar system; (ii) the orbit of the surviving embryo(s) typically has perihelion of 40--70 AU, semi-major axis less than 200 AU, and inclination less than 30 degrees; (iii) it is likely that any surviving embryos could be detected by current or planned optical surveys or have a significant effect on solar-system ephemerides; (iv) whether or not an embryo has survived to the present day, their dynamical influence earlier in the history of the solar system can explain the properties of the detached disk (defined in this paper as containing objects with perihelia > 38 AU and semi-major axes between 80 and 500 AU).
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Hot Jupiters driven by high-eccentricity migration in globular clusters
Authors:
Adrian S. Hamers,
Scott Tremaine
Abstract:
Hot Jupiters (HJs) are short-period giant planets that are observed around ~ 1% of solar-type field stars. One possible formation scenario for HJs is high-eccentricity (high-e) migration, in which the planet forms at much larger radii, is excited to high eccentricity by some mechanism, and migrates to its current orbit due to tidal dissipation occurring near periapsis. We consider high-e migration…
▽ More
Hot Jupiters (HJs) are short-period giant planets that are observed around ~ 1% of solar-type field stars. One possible formation scenario for HJs is high-eccentricity (high-e) migration, in which the planet forms at much larger radii, is excited to high eccentricity by some mechanism, and migrates to its current orbit due to tidal dissipation occurring near periapsis. We consider high-e migration in dense stellar systems such as the cores of globular clusters (GCs), in which encounters with passing stars can excite planets to the high eccentricities needed to initiate migration. We study this process via Monte Carlo simulations of encounters with a star+planet system including the effects of tidal dissipation, using an efficient regularized restricted three-body code. HJs are produced in our simulations over a significant range of the stellar number density n_*. Assuming the planet is initially on a low-eccentricity orbit with semimajor axis 1 au, for n_* < 1e3 pc^{-3} the encounter rate is too low to induce orbital migration, whereas for n_* > 1e6 pc^{-3} HJ formation is suppressed because the planet is more likely ejected from its host star, tidally disrupted, or transferred to a perturbing star. The fraction of planets that are converted to HJs peaks at ~ 2% for intermediate number densities of ~ 4e4 pc^{-3}. Warm Jupiters, giant planets with periods between 10 and 100 days, are produced in our simulations with an efficiency of up to ~ 0.5%. Our results suggest that HJs can form through high-e migration induced by stellar encounters in the centers of of dense GCs, but not in their outskirts where the densities are lower.
△ Less
Submitted 29 November, 2017; v1 submitted 29 September, 2017;
originally announced October 2017.
-
Collisional Fragmentation is Not a Barrier to Close-in Planet Formation
Authors:
Joshua Wallace,
Scott Tremaine,
John Chambers
Abstract:
Collisional fragmentation is shown to not be a barrier to rocky planet formation at small distances from the host star. Simple analytic arguments demonstrate that rocky planet formation via collisions of homogeneous gravity-dominated bodies is possible down to distances of order the Roche radius ($r_\mathrm{Roche}$). Extensive N-body simulations with initial bodies ${\gtrsim}1700$ km that include…
▽ More
Collisional fragmentation is shown to not be a barrier to rocky planet formation at small distances from the host star. Simple analytic arguments demonstrate that rocky planet formation via collisions of homogeneous gravity-dominated bodies is possible down to distances of order the Roche radius ($r_\mathrm{Roche}$). Extensive N-body simulations with initial bodies ${\gtrsim}1700$ km that include plausible models for fragmentation and merging of gravity-dominated bodies confirm this conclusion and demonstrate that rocky planet formation is possible down to ${\sim}$1.1 $r_\mathrm{Roche}$. At smaller distances, tidal effects cause collisions to be too fragmenting to allow mass build-up to a final, dynamically stable planetary system. We argue that even differentiated bodies can accumulate to form planets at distances that are not much larger than $r_\mathrm{Roche}$.
△ Less
Submitted 13 September, 2017; v1 submitted 24 May, 2017;
originally announced May 2017.
-
Isotropic-Nematic Phase Transitions in Gravitational Systems
Authors:
Zacharias Roupas,
Bence Kocsis,
Scott Tremaine
Abstract:
We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital…
▽ More
We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital eccentricities or semimajor axes relax. We show that the statistical mechanics of such systems exhibit a striking resemblance to liquid crystals, with analogous ordered-nematic and disordered-isotropic phases. The ordered phase consists of bodies orbiting in a disk in both directions, with the disk thickness depending on temperature, while the disordered phase corresponds to a nearly isotropic distribution of the orbit normals. We show that below a critical value of the total angular momentum, the system undergoes a first-order phase transition between the ordered and disordered phases. At the critical point the phase transition becomes second-order while for higher angular momenta there is a smooth crossover. We also find metastable equilibria containing two identical disks with mutual inclinations between $90^{\circ}$ and $180^\circ$.
△ Less
Submitted 17 June, 2017; v1 submitted 12 January, 2017;
originally announced January 2017.
-
Galactic disc profiles and a universal angular momentum distribution from statistical physics
Authors:
Jakob Herpich,
Scott Tremaine,
Hans-Walter Rix
Abstract:
We show that the stellar surface-brightness profiles in disc galaxies---observed to be approximately exponential---can be explained if radial migration efficiently scrambles the individual stars' angular momenta while conserving the circularity of the orbits and the total mass and angular momentum. In this case the disc's distribution of specific angular momenta $j$ should be near a maximum-entrop…
▽ More
We show that the stellar surface-brightness profiles in disc galaxies---observed to be approximately exponential---can be explained if radial migration efficiently scrambles the individual stars' angular momenta while conserving the circularity of the orbits and the total mass and angular momentum. In this case the disc's distribution of specific angular momenta $j$ should be near a maximum-entropy state and therefore approximately exponential, $dN\propto\exp(-j/\langle j\rangle)dj$. This distribution translates to a surface-density profile that is generally not an exponential function of radius: $Σ(R)\propto\exp[-R/R_e(R)]/(RR_e(R))(1+d\log v_c(R)/d\log R)$, for a rotation curve $v_c(R)$ and $R_e(R)\equiv\langle j\rangle/v_c(R)$. We show that such a profile matches the observed surface-brightness profiles of disc-dominated galaxies as well as the empirical exponential profile. Disc galaxies that exhibit population gradients cannot have fully reached a maximum-entropy state but appear to be close enough that their surface-brightness profiles are well-fit by this idealized model.
△ Less
Submitted 12 June, 2017; v1 submitted 9 December, 2016;
originally announced December 2016.
-
Ultralight scalars as cosmological dark matter
Authors:
Lam Hui,
Jeremiah P. Ostriker,
Scott Tremaine,
Edward Witten
Abstract:
An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM…
▽ More
An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $λ\sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim ρλ^3$ where $ρ$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.
△ Less
Submitted 24 February, 2017; v1 submitted 26 October, 2016;
originally announced October 2016.
-
Lidov-Kozai Cycles with Gravitational Radiation: Merging Black Holes in Isolated Triple Systems
Authors:
Kedron Silsbee,
Scott Tremaine
Abstract:
We show that a black-hole binary with an external companion can undergo Lidov-Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has mass comparable to the reduced mass of the binary and the companion orbit has semi-major axis within a factor of $\sim 10$ of the binary…
▽ More
We show that a black-hole binary with an external companion can undergo Lidov-Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has mass comparable to the reduced mass of the binary and the companion orbit has semi-major axis within a factor of $\sim 10$ of the binary semi-major axis. Using a simple population-synthesis model and 3-body simulations, we estimate the rate of mergers in triple black hole systems in the field to be about six per Gpc$^3$ per year in the absence of natal kicks during black hole formation. This value is within the low end of the 90\% credible interval for the total black-hole black-hole merger rate inferred from the current LIGO results. There are many uncertainties in these calculations, the largest of which is the unknown distribution of natal kicks. Even modest natal kicks of $40\mbox{km s}^{-1}$ will reduce the merger rate by a factor of 40. A few percent of these systems will have eccentricity greater than 0.999 when they first enter the frequency band detectable by aLIGO (above 10 Hz).
△ Less
Submitted 25 January, 2017; v1 submitted 26 August, 2016;
originally announced August 2016.
-
Modeling the Nearly Isotropic Comet Population in Anticipation of LSST Observations
Authors:
Kedron Silsbee,
Scott Tremaine
Abstract:
We run simulations to determine the expected distribution of orbital elements of nearly isotropic comets (NICs) in the outer solar system, assuming that these comets originate in the Oort Cloud at thousands of AU and are perturbed into the planetary region by the Galactic tide. We show that the Large Synoptic Survey Telescope (LSST) should detect and characterize the orbits of hundreds to thousand…
▽ More
We run simulations to determine the expected distribution of orbital elements of nearly isotropic comets (NICs) in the outer solar system, assuming that these comets originate in the Oort Cloud at thousands of AU and are perturbed into the planetary region by the Galactic tide. We show that the Large Synoptic Survey Telescope (LSST) should detect and characterize the orbits of hundreds to thousands of NICs with perihelion distance outside 5 AU. Observing NICs in the outer solar system is our only way of directly detecting comets from the inner Oort Cloud, as these comets are dynamically excluded from the inner solar system by the giant planets. Thus the distribution of orbital elements constrains the spatial distribution of comets in the Oort cloud and the environment in which the solar system formed. Additionally, comet orbits can be characterized more precisely when they are seen far from the Sun as they have not been affected by non-gravitational forces.
△ Less
Submitted 26 July, 2016;
originally announced July 2016.
-
Warm Jupiters from secular planet-planet interactions
Authors:
Cristobal Petrovich,
Scott Tremaine
Abstract:
Most warm Jupiters (gas-giant planets with $0.1~{\rm AU}\lesssim a \lesssim1$ AU) have pericenter distances that are too large for significant orbital migration by tidal friction. We study the possibility that the warm Jupiters are undergoing secular eccentricity oscillations excited by an outer companion (a planet or star) in an eccentric and/or mutually inclined orbit. In this model the warm Jup…
▽ More
Most warm Jupiters (gas-giant planets with $0.1~{\rm AU}\lesssim a \lesssim1$ AU) have pericenter distances that are too large for significant orbital migration by tidal friction. We study the possibility that the warm Jupiters are undergoing secular eccentricity oscillations excited by an outer companion (a planet or star) in an eccentric and/or mutually inclined orbit. In this model the warm Jupiters migrate periodically, in the high-eccentricity phase of the oscillation when the pericenter distance is small, but are typically observed at much lower eccentricities. We show that the steady-state eccentricity distribution of the warm Jupiters migrating by this mechanism is approximately flat, which is consistent with the observed distribution if and only if we restrict the sample to warm Jupiters that have outer companions detected by radial-velocity surveys. The eccentricity distribution of warm Jupiters without companions exhibits a peak at low eccentricities ($e\lesssim 0.2$) that must be explained by a different formation mechanism. Based on a population-synthesis study we find that high-eccentricity migration excited by an outer planetary companion (i) can account for $\sim20\%$ of the warm Jupiters and most of the warm Jupiters with $e\gtrsim 0.4$, a fraction that is consistent with the observed population of warm Jupiters with outer companions; (ii) can produce most of the observed population of hot Jupiters, with a semimajor axis distribution that matches the observations, but fails to account adequately for $\sim 60\%$ of hot Jupiters with projected obliquities $\lesssim 20^\circ$. Thus $\sim 20\%$ of the warm Jupiters and $\sim 60\%$ of the hot Jupiters can be produced by high-eccentricity migration. We also provide predictions for the expected mutual inclinations and spin-orbit angles of the planetary systems with hot and warm Jupiters produced by high-eccentricity migration.
△ Less
Submitted 5 July, 2016; v1 submitted 31 March, 2016;
originally announced April 2016.
-
Adiabatic black hole growth in Sérsic models of elliptical galaxies
Authors:
Naveen Jingade,
Tarun Deep Saini,
Scott Tremaine
Abstract:
We have examined the effect of slow growth of a central black hole on spherical galaxies that obey Sérsic or $R^{1/m}$ surface-brightness profiles. During such growth the actions of each stellar orbit are conserved, which allows us to compute the final distribution function if we assume that the initial distribution function is isotropic. We find that black-hole growth leads to a central cusp or `…
▽ More
We have examined the effect of slow growth of a central black hole on spherical galaxies that obey Sérsic or $R^{1/m}$ surface-brightness profiles. During such growth the actions of each stellar orbit are conserved, which allows us to compute the final distribution function if we assume that the initial distribution function is isotropic. We find that black-hole growth leads to a central cusp or ``excess light', in which the surface brightness varies with radius as $R^{-1.3}$ (with a weak dependence on Sérsic index $m$), the line-of-sight velocity dispersion varies as $R^{-1/2}$, and the velocity anisotropy is $β\simeq -0.24$ to $-0.28$ depending on $m$. The excess stellar mass in the cusp scales approximately linearly with the black-hole mass, and is typically 0.5--0.85 times the black-hole mass. This process may strongly influence the structure of nuclear star clusters if they contain black holes.
△ Less
Submitted 14 September, 2015; v1 submitted 26 August, 2015;
originally announced August 2015.
-
Comparison of simple mass estimators for slowly rotating elliptical galaxies
Authors:
N. Lyskova,
J. Thomas,
E. Churazov,
S. Tremaine,
T. Naab
Abstract:
We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) u…
▽ More
We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $σ_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{σ_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.
△ Less
Submitted 14 April, 2015;
originally announced April 2015.
-
The statistical mechanics of planet orbits
Authors:
Scott Tremaine
Abstract:
The final "giant-impact" phase of terrestrial planet formation is believed to begin with a large number of planetary "embryos" on nearly circular, coplanar orbits. Mutual gravitational interactions gradually excite their eccentricities until their orbits cross and they collide and merge; through this process the number of surviving bodies declines until the system contains a small number of planet…
▽ More
The final "giant-impact" phase of terrestrial planet formation is believed to begin with a large number of planetary "embryos" on nearly circular, coplanar orbits. Mutual gravitational interactions gradually excite their eccentricities until their orbits cross and they collide and merge; through this process the number of surviving bodies declines until the system contains a small number of planets on well-separated, stable orbits. In this paper we explore a simple statistical model for the orbit distribution of planets formed by this process, based on the sheared-sheet approximation and the ansatz that the planets explore uniformly all of the stable region of phase space. The model provides analytic predictions for the distribution of eccentricities and semimajor axis differences, correlations between orbital elements of nearby planets, and the complete N-planet distribution function, in terms of a single parameter, the "dynamical temperature", that is determined by the planetary masses. The predicted properties are generally consistent with N-body simulations of the giant-impact phase and with the distribution of semimajor axis differences in the Kepler catalog of extrasolar planets. A similar model may apply to the orbits of giant planets if these orbits are determined mainly by dynamical evolution after the planets have formed and the gas disk has disappeared.
△ Less
Submitted 10 June, 2015; v1 submitted 5 April, 2015;
originally announced April 2015.
-
A few of Michel Henon's contributions to dynamical astronomy
Authors:
Scott Tremaine
Abstract:
This article reviews Michel Henon's contributions to a diverse set of problems in astrophysical dynamics, including violent relaxation, Saturn's rings, roundoff error in orbit integrations, and planet formation.
This article reviews Michel Henon's contributions to a diverse set of problems in astrophysical dynamics, including violent relaxation, Saturn's rings, roundoff error in orbit integrations, and planet formation.
△ Less
Submitted 19 November, 2014; v1 submitted 17 November, 2014;
originally announced November 2014.
-
Advances in exoplanet science from Kepler
Authors:
Jack J. Lissauer,
Rebekah I. Dawson,
Scott Tremaine
Abstract:
Numerous telescopes and techniques have been used to find and study extrasolar planets, but none has been more successful than NASA's Kepler Space Telescope. Kepler has discovered the majority of known exoplanets, the smallest planets to orbit normal stars, and the worlds most likely to be similar to our home planet. Most importantly, Kepler has provided our first look at typical characteristics o…
▽ More
Numerous telescopes and techniques have been used to find and study extrasolar planets, but none has been more successful than NASA's Kepler Space Telescope. Kepler has discovered the majority of known exoplanets, the smallest planets to orbit normal stars, and the worlds most likely to be similar to our home planet. Most importantly, Kepler has provided our first look at typical characteristics of planets and planetary systems for planets with sizes as small as and orbits as large as those of the Earth.
△ Less
Submitted 4 September, 2014;
originally announced September 2014.
-
The odd couple: quasars and black holes
Authors:
Scott Tremaine
Abstract:
Quasars emit more energy than any other objects in the universe, yet are not much bigger than the solar system. We are almost certain that quasars are powered by giant black holes of up to $10^{10}$ times the mass of the Sun, and that black holes of between $10^6$ and $10^{10}$ solar masses---dead quasars---are present at the centers of most galaxies. Our own galaxy contains a black hole of…
▽ More
Quasars emit more energy than any other objects in the universe, yet are not much bigger than the solar system. We are almost certain that quasars are powered by giant black holes of up to $10^{10}$ times the mass of the Sun, and that black holes of between $10^6$ and $10^{10}$ solar masses---dead quasars---are present at the centers of most galaxies. Our own galaxy contains a black hole of $4.3\times10^6$ solar masses. The mass of the central black hole appears to be closely related to other properties of its host galaxy, such as the total mass in stars, but the origin of this relation and the role that black holes play in the formation of galaxies are still mysteries.
△ Less
Submitted 23 July, 2014;
originally announced July 2014.
-
Relativistic redshifts in quasar broad lines
Authors:
Scott Tremaine,
Yue Shen,
Xin Liu,
Abraham Loeb
Abstract:
The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital S…
▽ More
The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital Sky Survey DR7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ~30-45 degrees, consistent with simple AGN unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.
△ Less
Submitted 10 June, 2014;
originally announced June 2014.
-
A numerical study of vector resonant relaxation
Authors:
Bence Kocsis,
Scott Tremaine
Abstract:
Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplecti…
▽ More
Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplectic integrator, N-ring, which follows the dynamical evolution of a cluster of N stars through vector resonant relaxation, by averaging the pairwise interactions over the orbital period and periapsis-precession timescale. We use N-ring to follow the evolution of clusters containing over 10^4 stars for tens of relaxation times. Among other results, we find that the evolution is dominated by torques among stars with radially overlapping orbits, and that resonant relaxation can be modelled as a random walk of the orbit normals on the sphere, with angular step size ranging from 0.5-1 radian. The relaxation rate in a cluster with a fixed number of stars is proportional to the RMS mass of the stars. The RMS torque generated by the cluster stars is reduced below the torque between Kepler orbits due to apsidal precession and declines weakly with the eccentricity of the perturbed orbit. However since the angular momentum of an orbit also decreases with eccentricity, the relaxation rate is approximately eccentricity-independent for e<0.7 and grows rapidly with eccentricity for e>0.8. We quantify the relaxation using the autocorrelation function of the spherical multipole moments; this decays exponentially and the e-folding time may be identified with the vector resonant relaxation timescale.
△ Less
Submitted 9 January, 2015; v1 submitted 4 June, 2014;
originally announced June 2014.
-
The statistical mechanics of self-gravitating Keplerian disks
Authors:
Jihad Touma,
Scott Tremaine
Abstract:
We describe the dynamics and thermodynamics of collisionless particle disks orbiting a massive central body, in the case where the disk mass is small compared to the central mass, the self-gravity of the disk dominates the non-Keplerian force, and the spread in semi-major axes is small. We show that with plausible approximations such disks have logarithmic two-body interactions and a compact phase…
▽ More
We describe the dynamics and thermodynamics of collisionless particle disks orbiting a massive central body, in the case where the disk mass is small compared to the central mass, the self-gravity of the disk dominates the non-Keplerian force, and the spread in semi-major axes is small. We show that with plausible approximations such disks have logarithmic two-body interactions and a compact phase space, and therefore exhibit thermodynamics that are simpler than most other gravitating systems, which require a confining box and artificial softening of the potential at small scales to be thermodynamically well-behaved. We solve for the microcanonical axisymmetric thermal equilibria and demonstrate the existence of a symmetry-breaking bifurcation into lopsided equilibria. We discuss the relation between thermal and dynamical instability in these systems and draw connections to astrophysical settings, as well as to the wider subject of the statistical mechanics of particles with logarithmic long-range interactions, such as point vortices in two-dimensional fluids.
△ Less
Submitted 11 September, 2014; v1 submitted 21 January, 2014;
originally announced January 2014.
-
Scattering outcomes of close-in planets: constraints on planet migration
Authors:
Cristobal Petrovich,
Scott Tremaine,
Roman R. Rafikov
Abstract:
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 separati…
▽ More
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.
△ Less
Submitted 21 April, 2014; v1 submitted 17 January, 2014;
originally announced January 2014.
-
Constraining Sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. II. The Population with Kinematically Offset Broad Balmer Emission Lines
Authors:
Xin Liu,
Yue Shen,
Fuyan Bian,
Abraham Loeb,
Scott Tremaine
Abstract:
A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable…
▽ More
A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]
△ Less
Submitted 13 June, 2014; v1 submitted 23 December, 2013;
originally announced December 2013.
-
The Black Hole Mass and the Stellar Ring in NGC 3706
Authors:
Kayhan Gültekin,
Karl Gebhardt,
John Kormendy,
Tod R. Lauer,
Ralf Bender,
Scott Tremaine,
Douglas O. Richstone
Abstract:
We determine the mass of the nuclear black hole ($M$) in NGC 3706, an early type galaxy with a central surface brightness minimum arising from an apparent stellar ring, which is misaligned with respect to the galaxy's major axis at larger radii. We fit new HST/STIS and archival data with axisymmetric orbit models to determine $M$, mass-to-light ratio ($Υ_V$), and dark matter halo profile. The best…
▽ More
We determine the mass of the nuclear black hole ($M$) in NGC 3706, an early type galaxy with a central surface brightness minimum arising from an apparent stellar ring, which is misaligned with respect to the galaxy's major axis at larger radii. We fit new HST/STIS and archival data with axisymmetric orbit models to determine $M$, mass-to-light ratio ($Υ_V$), and dark matter halo profile. The best-fit model parameters with 1$σ$ uncertainties are $M = (6.0^{+0.7}_{-0.9}) \times 10^8\ M_{\scriptscriptstyle \odot}$ and $Υ_V = 6.0 \pm 0.2\ M_{\scriptscriptstyle \odot}\ L_{{\scriptscriptstyle \odot},V}^{-1}$ at an assumed distance of 46 Mpc. The models are inconsistent with no black hole at a significance of $Δχ^2 = 15.4$ and require a dark matter halo to adequately fit the kinematic data, but the fits are consistent with a large range of plausible dark matter halo parameters. The ring is inconsistent with a population of co-rotating stars on circular orbits, which would produce a narrow line-of-sight velocity distribution (LOSVD). Instead, the ring's LOSVD has a small value of $|V|/σ$, the ratio of mean velocity to velocity dispersion. Based on the observed low $|V|/σ$, our orbit modeling, and a kinematic decomposition of the ring from the bulge, we conclude that the stellar ring contains stars that orbit in both directions. We consider potential origins for this unique feature, including multiple tidal disruptions of stellar clusters, a change in the gravitational potential from triaxial to axisymmetric, resonant capture and inclining of orbits by a binary black hole, and multiple mergers leading to gas being funneled to the center of the galaxy.
△ Less
Submitted 17 December, 2013;
originally announced December 2013.
-
Why do Earth satellites stay up?
Authors:
Scott Tremaine,
Tomer Yavetz
Abstract:
Satellites in low Earth orbits must accurately conserve their orbital eccentricity, since a decrease in perigee of only 5-10% would cause them to crash. However, these satellites are subject to gravitational perturbations from the Earth's multipole moments, the Moon, and the Sun that are not spherically symmetric and hence do not conserve angular momentum, especially over the tens of thousands of…
▽ More
Satellites in low Earth orbits must accurately conserve their orbital eccentricity, since a decrease in perigee of only 5-10% would cause them to crash. However, these satellites are subject to gravitational perturbations from the Earth's multipole moments, the Moon, and the Sun that are not spherically symmetric and hence do not conserve angular momentum, especially over the tens of thousands of orbits made by a typical satellite. Why then do satellites not crash? We describe a vector-based analysis of the long-term behavior of satellite orbits and apply this to several toy systems containing a single non-Keplerian perturbing potential. If only the quadrupole potential from the Earth's equatorial bulge is present, all near-circular orbits are stable. If only the octupole potential is present, all such orbits are unstable. If only the lunar or solar potential is present, all near-circular orbits with inclinations to the ecliptic exceeding 39 degrees are unstable. We describe the behavior of satellites in the simultaneous presence of all of these perturbations and show that almost all low Earth orbits are stable because of an accidental property of the dominant quadrupole potential. We also relate these results to the phenomenon of Lidov-Kozai oscillations.
△ Less
Submitted 20 September, 2013;
originally announced September 2013.
-
Dynamics of warped accretion discs
Authors:
Scott Tremaine,
Shane W. Davis
Abstract:
Accretion discs are present around both stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. A wide variety of circumstantial evidence implies that many of these discs are warped. The standard Bardeen--Petterson model attributes the shape of the warp to the competition between Lense--Thirring torque from the central black hole and viscous angular-momen…
▽ More
Accretion discs are present around both stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. A wide variety of circumstantial evidence implies that many of these discs are warped. The standard Bardeen--Petterson model attributes the shape of the warp to the competition between Lense--Thirring torque from the central black hole and viscous angular-momentum transport within the disc. We show that this description is incomplete, and that torques from the companion star (for X-ray binaries) or the self-gravity of the disc (for active galactic nuclei) can play a major role in determining the properties of the warped disc. Including these effects leads to a rich set of new phenomena. For example, (i) when a companion star is present and the warp arises from a misalignment between the companion's orbital axis and the black hole's spin axis, there is no steady-state solution of the Pringle--Ogilvie equations for a thin warped disc when the viscosity falls below a critical value; (ii) in AGN accretion discs, the warp can excite short-wavelength bending waves that propagate inward with growing amplitude until they are damped by the disc viscosity. We show that both phenomena can occur for plausible values of the black hole and disc parameters, and briefly discuss their observational implications.
△ Less
Submitted 8 August, 2013;
originally announced August 2013.
-
Co-Evolution of Galactic Nuclei and Globular Cluster Systems
Authors:
Oleg Y. Gnedin,
Jeremiah P. Ostriker,
Scott Tremaine
Abstract:
We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant ell…
▽ More
We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation.
△ Less
Submitted 13 January, 2014; v1 submitted 31 July, 2013;
originally announced August 2013.
-
Constraining sub-Parsec Binary Supermassive Black Holes in Quasars with Multi-Epoch Spectroscopy. I. The General Quasar Population
Authors:
Yue Shen,
Xin Liu,
Abraham Loeb,
Scott Tremaine
Abstract:
We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0…
▽ More
We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]
△ Less
Submitted 18 August, 2013; v1 submitted 18 June, 2013;
originally announced June 2013.
-
Planets Near Mean-Motion Resonances
Authors:
Cristobal Petrovich,
Renu Malhotra,
Scott Tremaine
Abstract:
The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipat…
▽ More
The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.
△ Less
Submitted 23 November, 2012;
originally announced November 2012.
-
Gravitational Collapse in One Dimension
Authors:
A. E. Schulz,
Walter Dehnen,
Gerard Jungman,
Scott Tremaine
Abstract:
We simulate the evolution of one-dimensional gravitating collisionless systems from non- equilibrium initial conditions, similar to the conditions that lead to the formation of dark- matter halos in three dimensions. As in the case of 3D halo formation we find that initially cold, nearly homogeneous particle distributions collapse to approach a final equilibrium state with a universal density prof…
▽ More
We simulate the evolution of one-dimensional gravitating collisionless systems from non- equilibrium initial conditions, similar to the conditions that lead to the formation of dark- matter halos in three dimensions. As in the case of 3D halo formation we find that initially cold, nearly homogeneous particle distributions collapse to approach a final equilibrium state with a universal density profile. At small radii, this attractor exhibits a power-law behavior in density, ρ(x) \propto |x|^(-γ_crit), γ_crit \simeq 0.47, slightly but significantly shallower than the value γ = 1/2 suggested previously. This state develops from the initial conditions through a process of phase mixing and violent relaxation. This process preserves the energy ranks of particles. By warming the initial conditions, we illustrate a cross-over from this power-law final state to a final state containing a homogeneous core. We further show that inhomogeneous but cold power-law initial conditions, with initial exponent γ_i > γ_crit, do not evolve toward the attractor but reach a final state that retains their original power-law behavior in the interior of the profile, indicating a bifurcation in the final state as a function of the initial exponent. Our results rely on a high-fidelity event-driven simulation technique.
△ Less
Submitted 1 June, 2012;
originally announced June 2012.
-
On the local dark matter density
Authors:
Jo Bovy,
Scott Tremaine
Abstract:
An analysis of the kinematics of 412 stars at 1-4 kpc from the Galactic mid-plane by Moni Bidin et al. (2012) has claimed to derive a local density of dark matter that is an order of magnitude below standard expectations. We show that this result is incorrect and that it arises from the assumption that the mean azimuthal velocity of the stellar tracers is independent of Galactocentric radius at al…
▽ More
An analysis of the kinematics of 412 stars at 1-4 kpc from the Galactic mid-plane by Moni Bidin et al. (2012) has claimed to derive a local density of dark matter that is an order of magnitude below standard expectations. We show that this result is incorrect and that it arises from the assumption that the mean azimuthal velocity of the stellar tracers is independent of Galactocentric radius at all heights. We substitute the assumption, supported by data, that the circular speed is independent of radius in the mid-plane. We demonstrate that the assumption of constant mean azimuthal velocity is implausible by showing that it requires the circular velocity to drop more steeply than allowed by any plausible mass model, with or without dark matter, at large heights above the mid-plane. Using the approximation that the circular velocity curve is flat in the mid-plane, we find that the data imply a local dark-matter density of 0.008 +/- 0.003 Msun/pc^3 = 0.3 +/- 0.1 GeV/cm3, fully consistent with standard estimates of this quantity. This is the most robust direct measurement of the local dark-matter density to date.
△ Less
Submitted 29 June, 2012; v1 submitted 17 May, 2012;
originally announced May 2012.
-
Density waves in debris discs and galactic nuclei
Authors:
Mir Abbas Jalali,
Scott Tremaine
Abstract:
We study the linear perturbations of collisionless near-Keplerian discs. Such systems are models for debris discs around stars and the stellar discs surrounding supermassive black holes at the centres of galaxies. Using a finite-element method, we solve the linearized collisionless Boltzmann equation and Poisson's equation for a wide range of disc masses and rms orbital eccentricities to obtain th…
▽ More
We study the linear perturbations of collisionless near-Keplerian discs. Such systems are models for debris discs around stars and the stellar discs surrounding supermassive black holes at the centres of galaxies. Using a finite-element method, we solve the linearized collisionless Boltzmann equation and Poisson's equation for a wide range of disc masses and rms orbital eccentricities to obtain the eigenfrequencies and shapes of normal modes. We find that these discs can support large-scale `slow' modes, in which the frequency is proportional to the disc mass. Slow modes are present for arbitrarily small disc mass so long as the self-gravity of the disc is the dominant source of apsidal precession. We find that slow modes are of two general types: parent modes and hybrid child modes, the latter arising from resonant interactions between parent modes and singular van Kampen modes. The most prominent slow modes have azimuthal wavenumbers $m=1$ and $m=2$. We illustrate how slow modes in debris discs are excited during a fly-by of a neighbouring star. Many of the non-axisymmetric features seen in debris discs (clumps, eccentricity, spiral waves) that are commonly attributed to planets could instead arise from slow modes; the two hypotheses can be distinguished by long-term measurements of the pattern speed of the features.
△ Less
Submitted 1 January, 2012; v1 submitted 20 October, 2011;
originally announced October 2011.
-
Super-Eccentric Migrating Jupiters
Authors:
Aristotle Socrates,
Boaz Katz,
Subo Dong,
Scott Tremaine
Abstract:
An important class of formation theories for hot Jupiters involves the excitation of extreme orbital eccentricity (e=0.99 or even larger) followed by tidal dissipation at periastron passage that eventually circularizes the planetary orbit at a period less than 10 days. In a steady state, this mechanism requires the existence of a significant population of super-eccentric (e>0.9) migrating Jupiters…
▽ More
An important class of formation theories for hot Jupiters involves the excitation of extreme orbital eccentricity (e=0.99 or even larger) followed by tidal dissipation at periastron passage that eventually circularizes the planetary orbit at a period less than 10 days. In a steady state, this mechanism requires the existence of a significant population of super-eccentric (e>0.9) migrating Jupiters with long orbital periods and periastron distances of only a few stellar radii. For these super-eccentric planets, the periastron is fixed due to conservation of orbital angular momentum and the energy dissipated per orbit is constant, implying that the rate of change in semi-major axis a is \dot a \propto a^0.5 and consequently the number distribution satisfies dN/dlog a\propto a^0.5. If this formation process produces most hot Jupiters, Kepler should detect several super-eccentric migrating progenitors of hot Jupiters, allowing for a test of high-eccentricity migration scenarios.
△ Less
Submitted 7 October, 2011;
originally announced October 2011.
-
The statistics of multi-planet systems
Authors:
Scott Tremaine,
Subo Dong
Abstract:
We describe statistical methods for measuring the exoplanet multiplicity function - the fraction of host stars containing a given number of planets - from transit and radial-velocity surveys. The analysis is based on the approximation of separability - that the distribution of planetary parameters in an n-planet system is the product of identical 1-planet distributions. We review the evidence that…
▽ More
We describe statistical methods for measuring the exoplanet multiplicity function - the fraction of host stars containing a given number of planets - from transit and radial-velocity surveys. The analysis is based on the approximation of separability - that the distribution of planetary parameters in an n-planet system is the product of identical 1-planet distributions. We review the evidence that separability is a valid approximation for exoplanets. We show how to relate the observable multiplicity function in surveys with similar host-star populations but different sensitivities. We also show how to correct for geometrical selection effects to derive the multiplicity function from transit surveys if the distribution of relative inclinations is known. Applying these tools to the Kepler transit survey and radial-velocity surveys, we find that (i) the Kepler data alone do not constrain the mean inclination of multi-planet systems; even spherical distributions are allowed by the data but only if a small fraction of host stars contain large planet populations (> 30); (ii) comparing the Kepler and radial-velocity surveys shows that the mean inclination of multi-planet systems lies in the range 0-5 degrees; (iii) the multiplicity function of the Kepler planets is not well-determined by the present data.
△ Less
Submitted 27 June, 2011;
originally announced June 2011.