BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Tal Alexander appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
I briefly review advances in the understanding and modeling of relativistic stellar dynamics around massive black holes (MBHs) in galactic nuclei, following the inclusion of coherent relaxation and of secular processes in a new formal analytic description of the dynamics.
A massive black hole (MBH) consumes stars whose orbits evolve into the small phase-space volume of unstable orbits, the "loss-cone", which take them directly into the MBH, or close enough to interact strongly with it. The resulting phenomena: tidal heating and tidal disruption, binary capture and hyper-velocity star ejection, gravitational wave (GW) emission by inspiraling compact remnants, or hydrodynamical interactions with an accretion disk, are of interest as they can produce observable signatures and thereby reveal the existence of the MBH, affect its mass and spin evolution, probe strong gravity, and provide information on stars and gas near the MBH. The continuous loss of stars and the processes that resupply them shape the central stellar distribution. We investigate relativistic stellar dynamics near the loss-cone of a non-spinning MBH in steady-state analytically and by Monte Carlo simulations of the diffusion of the orbital parameters. These take into account Newtonian mass precession due to enclosed stellar mass, in-plane precession due to general relativity, dissipation by GW, uncorrelated two-body relaxation, correlated resonant relaxation (RR) and adiabatic invariance due to secular precession, using a rigorously derived description of correlated post-Newtonian dynamics in the diffusion limit. We argue that general maximal entropy considerations strongly constrain orbital diffusion in steady-state, irrespective of the relaxation mechanism. We identify the exact phase-space separatrix between plunges and inspirals, predict their steady-state rates, and verify they are robust under a wide range of assumptions. We derive the dependence of the rates on the mass of the MBH, show that the contribution of RR is small, and discuss special cases where unquenched RR in restricted volumes of phase-space may affect the steady-state substantially.
Mass accretion by black holes (BHs) is typically capped at the Eddington rate, when radiation's push balances gravity's pull. However, even exponential growth at the Eddington-limited e-folding time t_E ~ few x 0.01 billion years, is too slow to grow stellar-mass BH seeds into the supermassive luminous quasars that are observed when the universe is 1 billion years old. We propose a dynamical mechanism that can trigger supra-exponential accretion in the early universe, when a BH seed is trapped in a star cluster fed by the ubiquitous dense cold gas flows. The high gas opacity traps the accretion radiation, while the low-mass BH's random motions suppress the formation of a slowly-draining accretion disk. Supra-exponential growth can thus explain the puzzling emergence of supermassive BHs that power luminous quasars so soon after the Big Bang.
Stars around a massive black hole (MBH) move on nearly fixed Keplerian orbits, in a centrally-dominated potential. The random fluctuations of the discrete stellar background cause small potential perturbations, which accelerate the evolution of orbital angular momentum by resonant relaxation. This drives many phenomena near MBHs, such as extreme mass-ratio gravitational wave inspirals, the warping of accretion disks, and the formation of exotic stellar populations. We present here a formal statistical mechanics framework to analyze such systems, where the background potential is described as a correlated Gaussian noise. We derive the leading order, phase-averaged 3D stochastic Hamiltonian equations of motion, for evolving the orbital elements of a test star, and obtain the effective Fokker-Planck equation for a general correlated Gaussian noise, for evolving the stellar distribution function. We show that the evolution of angular momentum depends critically on the temporal smoothness of the background potential fluctuations. Smooth noise has a maximal variability frequency $ν_{\max}$. We show that in the presence of such noise, the evolution of the normalized angular momentum $j=\sqrt{1-e^{2}}$ of a relativistic test star, undergoing Schwarzschild (in-plane) General Relativistic precession with frequency $ν_{GR}/j^{2}$, is exponentially suppressed for $j<j_{b}$, where $ν_{GR}/j_{b}^{2}\simν_{\max}$, due to the adiabatic invariance of the precession against the slowly varying random background torques. This results in an effective Schwarzschild precession-induced barrier in angular momentum. When $j_{b}$ is large enough, this barrier can have significant dynamical implications for processes near the MBH.
Massive black holes (MBHs) in galactic nuclei are believed to be surrounded by a high density stellar cluster, whose mass is mostly in hard-to-detect faint stars and compact remnants. Such dark cusps dominate the dynamics near the MBH: a dark cusp in the Galactic center (GC) of the Milky Way would strongly affect orbital tests of General Relativity there; on cosmic scales, dark cusps set the rates of gravitational wave emission events from compact remnants that spiral into MBHs, and they modify the rates of tidal disruption events, to list only some implications. A recently discovered long-period massive young binary (P_12 <~ 1 yr, M_12 ~ O(100 M_sun), T_12 ~ 6x10^6 yr), only ~0.1 pc from the Galactic MBH (Pfuhl et al 2013), sets a lower bound on the 2-body relaxation timescale there, min t_rlx ~ (P_12/M_12)^(2/3)T_12 ~ 10^7 yr, and correspondingly, an upper bound on the stellar number density, max n ~ few x 10^8/<M_star^2> 1/pc^3, based on the binary's survival against evaporation by the dark cusp. However, a conservative dynamical estimate, the drain limit, implies t_rlx > O(10^8) yr. Such massive binaries are thus too short-lived and tightly bound to constrain a dense relaxed dark cusp. We explore here in detail the use of longer-period, less massive and longer-lived binaries (P_12 ~ few yr, M_12 ~ 2-4 M_sun, T_12 ~ 10^8-10^10 yr), presently just below the detection threshold, for probing the dark cusp, and develop the framework for translating their future detections among the giants in the GC into dynamical constraints.
The cross-correlation function (CCF) is commonly employed in the study of AGN, where it is used to probe the structure of the broad line region by line reverberation, to study the continuum emission mechanism by correlating multi-waveband light curves and to seek correlations between the variability and other AGN properties. The z -transformed discrete correlation function (ZDCF) is a new method for estimating the CCF of sparse, unevenly sampled light curves. Unlike the commonly used interpolation method, it does not assume that the light curves are smooth and it does provide errors on its estimates. The ZDCF corrects several biases of the discrete correlation function method of Edelson & Krolik (1988) by using equal population binning and Fisher's z -transform. These lead to a more robust and powerful method of estimating the CCF of sparse light curves of as few as 12 points. Two examples of light curve analysis with the ZDCF are presented. 1) The ZDCF estimate of the auto-correlation function is used to uncover a correlation between AGN magnitude and variability time scale in a small simulated sample of very sparse and irregularly sampled light curves. 2) A maximum likelihood function for the ZDCF peak location is used to estimate the time-lag between two light curves. Fortran 77 and 95 code implementations of the ZDCF and the maximum likelihood peak location (PLIKE) algorithms are freely available (see http://www.weizmann.ac.il/weizsites/tal/research/software/).
The disruption of a star by the tidal field of a massive black hole is the final outcome of a chain of complex dynamical processes in the host galaxy. I introduce the "loss cone problem", and describe the many theoretical and numerical challenges on the path of solving it. I review various dynamical channels by which stars can be supplied to a massive black hole, and the relevant dynamical relaxation / randomization mechanisms. I briefly mention some "exotic" tidal disruption scenarios, and conclude by discussing some new dynamical results that are changing our understanding of dynamics near a massive black hole, and may well be relevant for tidal disruption dynamics.
[abridged] Energy relaxation around a massive black hole (MBH) is key to establishing the dynamical state of galactic nuclei, and the nature of close stellar interactions with the MBH. The standard description of relaxation as diffusion provides a perturbative 2nd-order solution in the weak two-body interaction limit. We run N-body simulations and find that this solution fails to describe the non-Gaussian relaxation on short timescale, which is strongly influenced by extreme events even in the weak limit, and is thus difficult to characterize and measure. We derive a non-perturbative solution for relaxation as an anomalous diffusion process, and develop a robust estimation technique to measure it in simulations. These enable us to analyze and model our numerical results, and validate in detail, for the first time, this model of energy relaxation around an MBH on all timescales. We derive the relation between the energy diffusion time, t_E, and the time for a small perturbation to return to steady state, t_r, in a relaxed, single mass cusp around a MBH. We constrain the contribution of strong encounters, measure that of the weakest encounters, determine the value of the Coulomb logarithm, and provide a robust analytical estimate for t_E in a finite nuclear stellar cusp. We find that t_r ~ 10t_E ~(5/32)Q^2P_h/N_h log Q, where Q=M_bh/M_* is the MBH to star mass ratio, the orbital period P_h and number of stars N_h are evaluated at the energy scale corresponding to the MBH's sphere of influence, E_h=sigma_inf^2, where sigma_inf is the velocity dispersion far from the MBH. We conclude, using the observed cosmic M_bh/sigma correlation, that cusps around lower-mass MBHs (M_bh<10^7 Mo), which evolved passively over a Hubble time, should be relaxed. We consider the effects of anomalous energy diffusion on orbital perturbations of stars observed near the Galactic MBH.
An accreting massive black hole (MBH) in a galactic nucleus is surrounded by a dense stellar cluster. We analyze and simulate numerically the evolution of a thin accretion disk due to its internal viscous torques, due to the frame-dragging torques of a spinning MBH (the Bardeen-Petterson effect) and due to the orbit-averaged gravitational torques by the stars (Resonant Relaxation). We show that the evolution of the MBH mass accretion rate, the MBH spin growth rate, and the covering fraction of the disk relative to the central ionizing continuum source, are all strongly coupled to the stochastic fluctuations of the stellar potential via the warps that the stellar torques excite in the disk. These lead to fluctuations by factors of up to a few in these quantities over a wide range of timescales, with most of the power on timescales >~(M_bh/M_d)P(R_d), where M_bh and M_d are the masses of the MBH and disk, and P is the orbital period at the disk's mass-weighted mean radius R_d. The response of the disk is stronger the lighter it is and the more centrally concentrated the stellar cusp. As proof of concept, we simulate the evolution of the low-mass maser disk in NGC4258, and show that its observed O(10 deg) warp can be driven by the stellar torques. We also show that the frame-dragging of a massive AGN disk couples the stochastic stellar torques to the MBH spin and can excite a jitter of a few degrees in its direction relative to that of the disk's outer regions.
Inspiral of compact stellar remnants into massive black holes (MBHs) is accompanied by the emission of gravitational waves at frequencies that are potentially detectable by space-based interferometers. Event rates computed from statistical (Fokker-Planck, Monte-Carlo) approaches span a wide range due to uncertaintities about the rate coefficients. Here we present results from direct integration of the post-Newtonian N-body equations of motion descrbing dense clusters of compact stars around Schwarzschild MBHs. These simulations embody an essentially exact (at the post-Newtonian level) treatment of the interplay between stellar dynamical relaxation, relativistic precession, and gravitational-wave energy loss. The rate of capture of stars by the MBH is found to be greatly reduced by relativistic precession, which limits the ability of torques from the stellar potential to change orbital angular momenta. Penetration of this "Schwarzschild barrier" does occasionally occur, resulting in capture of stars onto orbits that gradually inspiral due to gravitational wave emission; we discuss two mechanisms for barrier penetration and find evidence for both in the simulations. We derive an approximate formula for the capture rate, which predicts that captures would be strongly disfavored from orbits with semi-major axes below a certain value; this prediction, as well as the predicted rate, are verified in the N-body integrations. We discuss the implications of our results for the detection of extreme-mass-ratio inspirals from galactic nuclei with a range of physical properties.
I discuss four key questions about Galactic Center dynamics, their implications for understanding both the environment of the Galactic MBH and galactic nuclei in general, and the progress made in addressing them. The questions are (1) Is the stellar system around the MBH relaxed? (2) Is there a "dark cusp" around the MBH? (3) What is the origin of the stellar disk(s)?, and (4) What is the origin of the S-stars?
We analyze the distribution of stars of arbitrary mass function xi(m) around a massive black hole (MBH). Unless xi is strongly dominated by light stars, the steady-state distribution function approaches a power-law in specific energy x=-E/(m*sigma^2)<x_max with index p=m/4M_0, where E is the energy, sigma is the typical velocity dispersion of unbound stars, and M_0 is the mass averaged over m*xi*x_{max}^p. For light-dominated xi, p can grow as large as 3/2 - much steeper than previously thought. A simple prescription for the stellar density profile around MBHs is provided. We illustrate our results by applying them to stars around the MBH in the Milky Way.
The spin and quadrupole moment of the supermassive black hole at the Galactic center can in principle be measured via astrometric monitoring of stars orbiting at milliparsec (mpc) distances, allowing tests of general relativistic "no-hair" theorems (Will 2008). One complicating factor is the presence of perturbations from other stars, which may induce orbital precession of the same order of magnitude as that due to general relativistic effects. The expected number of stars in this region is small enough that full N-body simulations can be carried out. We present the results of a comprehensive set of such simulations, which include a post-Newtonian treatment of spin-orbit effects. A number of possible models for the distribution of stars and stellar remnants are considered. We find that stellar perturbations are likely to obscure the signal due to frame-dragging for stars beyond ~0.5 mpc from the black hole, while measurement of the quadrupole moment is likely to require observation of stars inside ~0.2 mpc. A high fraction of stellar remnants, e.g. 10-Solar-mass black holes, in this region would make tests of GR problematic at all radii. We discuss the possibility of separating the effects of stellar perturbations from those due to GR.
In recent years several hypervelocity stars (HVSs) have been observed in the halo of our Galaxy. Such HVSs have possibly been ejected from the Galactic center and then propagated in the Galactic potential up to their current position. The recent survey for candidate HVSs show an asymmetry in the kinematics of candidate HVSs (position and velocity vectors), where more outgoing stars than ingoing stars (i.e. positive Galactocentric velocities vs. negative ones) are observed. We show that such kinematic asymmetry, which is likely due to the finite lifetime of the stars and Galactic potential structure, could be used in a novel method to probe and constrain the Galactic potential, identify the stellar type of the stars in the survey and estimate the number of HVSs. Kinematics-independent identification of the stellar types of the stars in such surveys (e.g. spectroscopic identification) could further improve these results. We find that the observed asymmetry between ingoing and outgoing stars favors specific Galactic potential models. It also implies a lower limit of ~54+-8 main sequence HVSs in the survey sample (>=648+-96 in the Galaxy), assuming that all of the main sequence stars in the survey originate from the Galactic center. The other stars in the survey are likely to be hot blue horizontal branch stars born in the halo rather than stars ejected from the Galactic center.
We show that the mass-segregation solution for the steady state distribution of stars around a massive black hole (MBH) has two branches: the known weak segregation solution (Bahcall & Wolf 1977), and a newly discovered strong segregation solution, presented here. The nature of the solution depends on the heavy-to-light stellar mass ratio M_H/M_L and on the unbound population number ratio N_H/N_L, through the relaxational coupling parameter Δ=4 N_H M_H^2 /[N_L M_L^2(3+M_H/M_L)]. When the heavy stars are relatively common (Δ>>1), they scatter frequently on each other. This efficient self-coupling leads to weak mass segregation, where the stars form n \propto r^{-α_M} mass-dependent cusps near the MBH, with indices α_H=7/4 for the heavy stars and 3/2<α_L<7/4 for the light stars (i.e. \max(α_H-α_L)~=1/4). However, when the heavy stars are relatively rare (Δ<<1), they scatter mostly on light stars, sink to the center by dynamical friction and settle into a much steeper cusp with 2~<α_H<11/4, while the light stars form a 3/2<α_L<7/4 cusp, resulting in strong segregation (i.e. \max(α_H-α_L)~=1). We show that the present-day mass function of evolved stellar populations (coeval or continuously star forming) with a universal initial mass function, separate into two distinct mass scales, ~1 Mo of main sequence and compact dwarfs, and ~10 Mo of stellar black holes (SBHs), and have Δ<0.1. We conclude that it is likely that many relaxed galactic nuclei are strongly segregated. We review indications of strong segregation in observations of the Galactic Center and in results of numeric simulations, and briefly list some possible implications of a very high central concentration of SBHs around a MBH.
Resonant relaxation (RR) is a rapid relaxation process that operates in the nearly-Keplerian potential near a massive black hole (MBH). RR dominates the dynamics of compact remnants that inspiral into a MBH and emit gravitational waves (extreme mass ratio inspiral events, EMRIs). RR can either increase the EMRI rate, or strongly suppress it, depending on its still poorly-determined efficiency. We use small-scale Newtonian N-body simulations to measure the RR efficiency and to explore its possible dependence on the stellar number density profile around the MBH, and the mass-ratio between the MBH and a star (a single-mass stellar population is assumed). We develop an efficient and robust procedure for detecting and measuring RR in N-body simulations. We present a suite of simulations with a range of stellar density profiles and mass-ratios, and measure the mean RR efficiency in the near-Keplerian limit. We do not find a statistically significant dependence on the density profile or the mass-ratio. Our numerical determination of the RR efficiency in the Newtonian, single-mass population approximations, suggests that RR will likely enhance the EMRI rate by a factor of a few over the rates predicted assuming only slow stochastic two-body relaxation.
The new photometric space-borne survey missions CoRoT and Kepler will be able to detect minute flux variations in binary stars due to relativistic beaming caused by the line-of-sight motion of their components. In all but very short period binaries (P>10d), these variations will dominate over the ellipsoidal and reflection periodic variability. Thus, CoRoT and Kepler will discover a new observational class: photometric beaming binary stars. We examine this new category and the information that the photometric variations can provide. The variations that result from the observatory heliocentric velocity can be used to extract some spectral information even for single stars.
We study the role of massive perturbers (MPs) in deflecting stars and binaries to almost radial ("loss-cone") orbits, where they pass near the central massive black hole (MBH), interact with it at periapse, and are ultimately destroyed. MPs dominate dynamical relaxation when the ratio of the 2nd moments of the MP and star mass distributions, mu_2=(N_p*<M_p^2>)/(N_s*<M_s^2>), satisfies mu_2>>1. We compile the MP mass function from published observations, and show that MPs in the nucleus of the Galaxy (mainly giant molecular clouds), and plausibly in late type galaxies generally, have 10^2<mu_2<10^8. MPs thus shorten the relaxation timescale by 10^1-10^7 relative to 2-body relaxation by stars alone. We show this increases by 10-1000 the rate of \emph{large}-periapse interactions with the MBH, where loss-cone refilling by stellar 2-body relaxation is inefficient. We extend the Fokker-Planck loss-cone formalism to approximately account for relaxation by rare encounters with MPs. We show that binary stars--MBH exchanges driven by MPs can explain the origin of the young main sequence B stars that are observed very near the Galactic MBH, and increase by orders of magnitude the ejection rate of hyper-velocity stars. In contrast, the rate of \emph{small}-periapse interactions of single stars with the MBH, such as tidal disruption, is only increased by a factor of a few. We suggest that MP-driven relaxation plays an important role in the 3-body exchange capture of single stars on very tight orbits around the MBH. These captured stars may later be disrupted by the MBH via tidal orbital decay or direct scattering into the loss cone; captured compact objects may inspiral into the MBH by the emission of gravitational waves from zero-eccentricity orbits.
Gravitational waves (GWs) from the inspiral of compact remnants (CRs) into massive black holes (MBHs) will be observable to cosmological distances. While a CR spirals in, 2-body scattering by field stars may cause it to fall into the MBH before reaching a short period orbit that would give an observable signal. As a result, only CRs very near (~0.01 pc) the MBH can spiral in successfully. In a multi-mass stellar population, the heaviest objects sink to the center, where they are more likely to slowly spiral into the MBH without being swallowed prematurely. We study how mass-segregation modifies the stellar distribution and the rate of GW events. We find that the inspiral rate per galaxy for white dwarfs is 30 per Gyr, for neutron stars 6 per Gyr, and for stellar black holes (SBHs) 250 per Gyr. The high rate for SBHs is due to their extremely steep density profile, n_{BH}(r)\propto r^{-2}. The GW detection rate will be dominated by SBHs.
The Bloated Stars Scenario proposes that AGN broad line emission originates in the winds or envelopes of bloated stars (BS). Alexander and Netzer (1994) established that ~ 5e4 BSs with dense, decelerating winds can reproduce the observed emission line spectrum and avoid rapid collisional destruction. Here, we use the observed properties of AGN line profiles to further constrain the model parameters. In the BS model, the origin of the broad profiles is the stellar velocity field in the vicinity of the central black hole. We use a detailed photoionization code and a model of the stellar distribution function to calculate the BS emission line profiles and compare them to a large sample of AGN CIV, CIII] and MgII profiles. We find that the BSs can reproduce the general shape and width of typical AGN profiles as well as the line ratios if (i) The ionizing luminosity to black hole mass ratio is low enough. (ii) The broad line region size is limited by some cutoff mechanism. (iii) The fraction of the BSs in the stellar population falls off roughly as r^-2. (iv) The wind density and/or velocity are correlated with the black hole mass and ionizing luminosity. Under these conditions the strong tidal forces near the black hole play an important role in determining the line emission properties of the BSs. Some discrepancies remain: the calculated BS profiles tend to have weaker wings than the observed ones, and the differences between the profiles of different lines are somewhat smaller than those observed.
The `Bloated Stars Scenario' proposes that AGN broad line emission originates in the winds or envelopes of bloated stars (BS). Alexander and Netzer (1994, 1996) established that ~ 5e4 BSs with dense, decelerating winds can reproduce the observed emission line spectrum and line profiles while avoiding rapid collisional destruction. Here, I investigate a third prediction of the model, related to the size of the broad line region, by deriving the emission line response to variations in the ionizing continuum (`line reverberation') and comparing it to observations. The expected time lags, as well as the order of response of the various lines, strongly depend on the typical variability time scale of the ionizing continuum. The BS model studied here corresponds to a bright Syfert 1 galaxy (Sy1) or a low luminosity QSO. I find that the BS model is consistent with the observed correlation between the Balmer lines time lags and the AGN luminosity, which at present is the only line reverberation information available for this luminosity class. The model displays also some of the anti-correlation between the time lags of the metal emission lines and their degree of ionization that has been observed in a few low-luminosity Sy1s. However, the bright Sy1 model results differs from the low-luminosity Sy1 data in that the Hb time lag is relatively shorter and the CIV time lag relatively longer.
A recently discovered inverse correlation between QSO redshift and long-term continuum variability timescales was suggested to be the signature of microlensing on cosmological scales (Hawkins 1993). A general theoretical method for calculating such correlations is presented and applied to various lensing scenarios in the framework of $Λ= 0$ Friedmann cosmologies. It is shown that the observed timescales can be strongly influenced by the observational limitations: the finite duration of the monitoring campaign and the finite photometric sensitivity. In most scenarios the timescales increase with source redshift, $z_s$, although slower than the $1+z_s$ time dilation expected of intrinsic variability. A decrease can be obtained for an extended source observed with moderate sensitivity. In this case, only lenses no further away than several hundreds Mpc participate in the lensing. The resulting optical depth is too small to explain the common long-term QSO variability unless an extremely high local lens density is assumed. These results do not support the idea that the reported inverse correlation can be attributed to microlensing of a uniform QSO sample by a uniform distribution of lenses. The possibility of using observations at various wavelengths and QSO samples at various positions to identify microlensing in QSO variability is also discussed.
The Bloated Stars Scenario proposes that AGN broad line emission originates in the winds or envelopes of bloated stars (BS). Its main advantage over BLR cloud models is the gravitational confinement of the gas and its major difficulty the large estimated number of BSs and resulting high mass loss rate. We calculate the emission line spectrum by a detailed numerical photoionization code for a wide range of wind structures and a detailed QSO nucleus model with L(ion)=7E45 erg/s, M(bh)=8E7 Mo. The size and boundary density of the wind are determined by various processes: Comptonization by the central ionizing source (calculated self consistently), tidal disruption by the black hole and the limit set by the wind's finite mass. We find that the emission spectrum is mainly determined by the conditions at the boundary of the line emitting fraction of the wind rather than by its internal structure. Comptonization results in very high ionization parameters at the boundary which produces an excess of unobserved broad high excitation forbidden lines. The finite mass constraint limits the wind size, increases the boundary density and thus improves the results. Slow decelerating mass-constrained flows with high gas densities (1E8-1E12 cm^-3) are as success- ful as clouds in reproducing the observed lines. Only 5E4 dense wind BSs are required in the inner 1/3 pc. This small fraction approaches that of super- giants (SG) in the solar neighborhood and the mass loss is consistent with the observational constraints. The BS number is a very sensitive function of the wind structure. SG-like BSs are ruled out by their high mass loss. BSs with dense winds can reproduce the BLR line spectrum and be supported by the stellar population without excessive mass loss and collisional destruction rates.