Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter search to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The "standard" EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition to it, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) range. The first is due to bremsstrahlung of electrons scattered on neutral atoms (so-called neutral bremsstrahlung, NBrS). The second, responsible for electron avalanche scintillations in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a SiPM matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillations effects. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter search. It should be emphasized that this is the first practical application of the NBrS effect in detection science.

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $\beta$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $\gamma$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector.

We report measurements of the charged daughter fraction of $^{218}$Po as a result of the $^{222}$Rn alpha decay, and the mobility of $^{218}$Po$^+$ ions, using radon-polonium coincidences from the $^{238}$U chain identified in 532 live-days of DarkSide-50 WIMP-search data. The fraction of $^{218}$Po that is charged is found to be 0.37$\pm$0.03 and the mobility of $^{218}$Po$^+$ is (8.6$\pm$0.1)$\times$10$^{-4}$$\frac{\text{cm}^2}{\text{Vs}}$.

We have experimentally determined the production rate of $^{39}$Ar and $^{37}$Ar from cosmic ray neutron interactions in argon at sea level. Understanding these production rates is important for argon-based dark matter experiments that plan to utilize argon extracted from deep underground because it is imperative to know what the ingrowth of $^{39}$Ar will be during the production, transport, and storage of the underground argon. These measurements also allow for the prediction of $^{39}$Ar and $^{37}$Ar concentrations in the atmosphere which can be used to determine the presence of other sources of these isotopes. Through controlled irradiation with a neutron beam that mimics the cosmic ray neutron spectrum, followed by direct counting of $^{39}$Ar and $^{37}$Ar decays with sensitive ultra-low background proportional counters, we determined that the production rate from cosmic ray neutrons at sea-level is expected to be $(759 \pm 128)$ atoms/kg$_\text{Ar}$/day for $^{39}$Ar, and $(51.0 \pm 7.4)$ atoms/kg$_\text{Ar}$/day for $^{37}$Ar. We also performed a survey of the alternate production mechanisms based on the state-of-knowledge of the associated cross-sections to obtain a total sea-level cosmic ray production rate of $(1048 \pm 133)$ atoms/kg$_\text{Ar}$/day for $^{39}$Ar, $(56.7 \pm 7.5)$ atoms/kg$_\text{Ar}$/day for $^{37}$Ar in underground argon, and $(92 \pm 13)$ atoms/kg$_\text{Ar}$/day for $^{37}$Ar in atmospheric argon.

The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

Intermediate-mass black holes (IMBHs) of mass $M_{\bullet} \approx 10^{2} - 10^{5}$ solar masses, $M_{\odot}$, are the long-sought missing link between stellar black holes, born of supernovae, and massive black holes, tied to galaxy evolution by the empirical $M_{\bullet}/\sigma_{\star}$ correlation. We show that low-mass black hole seeds that accrete stars from locally dense environments in galaxies following a universal $M_{\bullet}/\sigma_{\star}$ relation grow over the age of the Universe to be above ${\mathcal{M}}_{0}\approx3\times10^{5}M_{\odot}$ ($5\%$ lower limit), independent of the unknown seed masses and formation processes. The mass ${\mathcal{M}}_{0}$ depends weakly on the uncertain formation redshift, and sets a universal minimal mass scale for present-day black holes. This can explain why no IMBHs have yet been found, and it implies that present-day galaxies with ${\sigma_{\star}<{\mathcal{S}}_{0}\approx40\,\mathrm{km\,s}^{-1}}$ lack a central black hole, or formed it only recently. A dearth of IMBHs at low redshifts has observable implications for tidal disruptions and gravitational wave mergers.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

Extreme mass ratio inspiral (EMRI) events are vulnerable to perturbations by the stellar background, which can abort them prematurely by deflecting EMRI orbits to plunging ones that fall directly into the massive black hole (MBH), or to less eccentric ones that no longer interact strongly with the MBH. A coincidental hierarchy between the collective resonant Newtonian torques due to the stellar background, and the relative magnitudes of the leading-order post-Newtonian precessional and radiative terms of the general relativistic 2-body problem, allows EMRIs to decouple from the background and produce semi-periodic gravitational wave signals. I review the recent theoretical developments that confirm this conjectured fortunate coincidence, and briefly discuss the implications for EMRI rates, and show how these dynamical effects can be probed locally by stars near the Galactic MBH.

Most galactic nuclei harbor a massive black hole (MBH), whose birth and evolution are closely linked to those of its host galaxy. The unique conditions near the MBH: high velocity and density in the steep potential of a massive singular relativistic object, lead to unusual modes of stellar birth, evolution, dynamics and death. A complex network of dynamical mechanisms, operating on multiple timescales, deflect stars to orbits that intercept the MBH. Such close encounters lead to energetic interactions with observable signatures and consequences for the evolution of the MBH and its stellar environment. Galactic nuclei are astrophysical laboratories that test and challenge our understanding of MBH formation, strong gravity, stellar dynamics, and stellar physics. I review from a theoretical perspective the wide range of stellar phenomena that occur near MBHs, focusing on the role of stellar dynamics near an isolated MBH in a relaxed stellar cusp.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

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.

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.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

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 $\nu_{\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 $\nu_{GR}/j^{2}$, is exponentially suppressed for $j<j_{b}$, where $\nu_{GR}/j_{b}^{2}\sim\nu_{\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.

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.

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.

As part of the DarkSide program of direct dark matter searches using liquid argon TPCs, a prototype detector with an active volume containing 10 kg of liquid argon, DarkSide-10, was built and operated underground in the Gran Sasso National Laboratory in Italy. A critically important parameter for such devices is the scintillation light yield, as photon statistics limits the rejection of electron-recoil backgrounds by pulse shape discrimination. We have measured the light yield of DarkSide-10 using the readily-identifiable full-absorption peaks from gamma ray sources combined with single-photoelectron calibrations using low-occupancy laser pulses. For gamma lines of energies in the range 122-1275 keV, we get consistent light yields averaging 8.887+-0.003(stat)+-0.444(sys) p.e./keVee. With additional purification, the light yield measured at 511 keV increased to 9.142+-0.006(stat) p.e./keVee.

A long-term spectroscopic and photometric survey of the most luminous and massive stars in the vicinity of the super-massive black hole Sgr A* revealed two new binaries; a long-period Ofpe/WN9 binary, GCIRS 16NE, with a modest eccentricity of 0.3 and a period of 224 days and an eclipsing Wolf-Rayet binary with a period of 2.3 days. Together with the already identified binary GCIRS 16SW, there are now three confirmed OB/WR binaries in the inner 0.2\,pc of the Galactic Center. Using radial velocity change upper limits, we were able to constrain the spectroscopic binary fraction in the Galactic Center to $F_{\rm SB}=0.27^{+0.29}_{-0.19}$ at a confidence level of 95%, a massive binary fraction similar to that observed in dense clusters. The fraction of eclipsing binaries with photometric amplitudes $\Delta m>0.4$ is $F^{\rm GC}_{\rm EB}=3\pm2%$, which is consistent with local OB star clusters ($F_{\rm EB}=1%$). Overall the Galactic Center binary fraction seems to be close to the binary fraction in comparable young clusters.

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.

We report on the performance and commissioning of a cryogenic distillation column for low radioactivity underground argon at Fermi National Accelerator Laboratory. The distillation column is designed to accept a mixture of argon, helium, and nitrogen and return pure argon with a nitrogen contamination less than 10 ppm. In the first commissioning, we were able to run the distillation column in a continuous mode and produce argon that is 99.9% pure. After running in a batch mode, the argon purity was increased to 99.95%, with 500 ppm of nitrogen remaining. The efficiency of collecting the argon from the gas mixture was between 70% and 81%, at an argon production rate of 0.84-0.98 kg/day.

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.

The ultra-luminous intermediate-mass black-hole system HLX-1 in the ESO 243-49 galaxy exhibits variability with a possible recurrence time of a few hundred days. Finding the origin of this variability would constrain the still largely unknown properties of this extraordinary object. Since it exhibits an intensity-hardness behavior characteristic of black-hole X-ray transients, we have analyzed the variability of HLX-1 in the framework of the disk instability model that explains outbursts of such systems. We find that the long-term variability of HLX-1 is unlikely to be explained by a model in which outbursts are triggered by thermal-viscous instabilities in an accretion disc. Possible alternatives include the instability in a radiation-pressure dominated disk but we argue that a more likely explanation is a modulated mass-transfer due to tidal stripping of a star on an eccentric orbit around the intermediate-mass black hole. We consider an evolutionary scenario leading to the creation of such a system and estimate the probability of its observation. We conclude, using a simplified dynamical model of the post-collapse cluster, that no more than 1/100 to 1/10 of Mbh < 10^4 Msun IMBHs - formed by run-away stellar mergers in the dense collapsed cores of young clusters - could have a few times 1 Msun Main-Sequence star evolve to an AGB on an orbit eccentric enough for mass transfer at periapse, while avoiding collisional destruction or being scattered into the IMBH by 2-body encounters. The finite but low probability of this configuration is consistent with the uniqueness of HLX-1. We note, however, that the actual response of a standard accretion disk to bursts of mass transfer may be too slow to explain the observations unless the orbit is close to parabolic (and hence even rarer) and/or additional heating, presumably linked to the highly time-dependent gravitational potential, are invoked.

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.

We present new observations of the nuclear star cluster in the central parsec of the Galaxy with the adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our work allows the spectroscopic detection of early and late type stars to m_K >= 16, more than 2 magnitudes deeper than our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We find that most of these Wolf Rayet (WR), O- and B- stars reside in two strongly warped disks between 0.8" and 12" from SgrA*, as well as a central compact concentration (the S-star cluster) centered on SgrA*. The later type B stars (m_K>15) in the radial interval between 0.8" and 12" seem to be in a more isotropic distribution outside the disks. The observed dearth of late type stars in the central few arcseconds is puzzling, even when allowing for stellar collisions. The stellar mass function of the disk stars is extremely top heavy with a best fit power law of dN/dm ~ m^(-0.45+/-0.3). Since at least the WR/O-stars were formed in situ in a single star formation event ~6 Myrs ago, this mass function probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and of the early-type stars at distances beyond 12" are compatible with a standard Salpeter/Kroupa IMF (best fit power law of dN/dm ~ m^(-2.15+/-0.3)).

The maser disk around the massive black hole (MBH) in active galaxy NGC 4258 exhibits an O(10 deg) warp on the O(0.1 pc) scale. The physics driving the warp are still debated. Suggested mechanisms include torquing by relativistic frame dragging or by radiation pressure. We propose here a new warping mechanism: resonant torquing of the disk by stars in the dense cusp around the MBH. We show that resonant torquing can induce such a warp over a wide range of observed and deduced physical parameters of the maser disk.

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 use N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (<0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration; and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The N-body simulations include a dense cluster of 10 M_sun stellar black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent at the ~95 % level with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20 % - 30 % of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with stellar black holes.

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

Electromagnetic observations over the last 15 years have yielded a growing appreciation for the importance of supermassive black holes (SMBH) to the evolution of galaxies, and for the intricacies of dynamical interactions in our own Galactic center. Here we show that future low-frequency gravitational wave observations, alone or in combination with electromagnetic data, will open up unique windows to these processes. In particular, gravitational wave detections in the 10^{-5}-10^{-1} Hz range will yield SMBH masses and spins to unprecedented precision and will provide clues to the properties of the otherwise undetectable stellar remnants expected to populate the centers of galaxies. Such observations are therefore keys to understanding the interplay between SMBHs and their environments.

We present the results of 16 years of monitoring stellar orbits around the massive black hole in center of the Milky Way using high resolution NIR techniques. This work refines our previous analysis mainly by greatly improving the definition of the coordinate system, which reaches a long-term astrometric accuracy of 300 microarcsecond, and by investigating in detail the individual systematic error contributions. The combination of a long time baseline and the excellent astrometric accuracy of adaptive optics data allow us to determine orbits of 28 stars, including the star S2, which has completed a full revolution since our monitoring began. Our main results are: all stellar orbits are fit extremely well by a single point mass potential to within the astrometric uncertainties, which are now 6 times better than in previous studies. The central object mass is (4.31 +- 0.06|stat +- 0.36|R0) * 10^6 M_sun where the fractional statistical error of 1.5 percent is nearly independent from R0 and the main uncertainty is due to the uncertainty in R0. Our current best estimate for the distance to the Galactic Center is R0 = 8.33 +- 0.35 kpc. The dominant errors in this value is systematic. The mass scales with distance as (3.95 +- 0.06) * 10^6 M_sun * (R0/8kpc)^2.19. The orientations of orbital angular momenta for stars in the central arcsecond are random. We identify six of the stars with orbital solutions as late type stars, and six early-type stars as members of the clockwise rotating disk system, as was previously proposed. We constrain the extended dark mass enclosed between the pericenter and apocenter of S2 at less than 0.066, at the 99% confidence level, of the mass of Sgr A*. This is two orders of magnitudes larger than what one would expect from other theoretical and observational estimates.

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 \Delta=4 N_H M_H^2 /[N_L M_L^2(3+M_H/M_L)]. When the heavy stars are relatively common (\Delta>>1), they scatter frequently on each other. This efficient self-coupling leads to weak mass segregation, where the stars form n \propto r^{-\alpha_M} mass-dependent cusps near the MBH, with indices \alpha_H=7/4 for the heavy stars and 3/2<\alpha_L<7/4 for the light stars (i.e. \max(\alpha_H-\alpha_L)~=1/4). However, when the heavy stars are relatively rare (\Delta<<1), they scatter mostly on light stars, sink to the center by dynamical friction and settle into a much steeper cusp with 2~<\alpha_H<11/4, while the light stars form a 3/2<\alpha_L<7/4 cusp, resulting in strong segregation (i.e. \max(\alpha_H-\alpha_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 \Delta<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.

The massive Galactic black hole and the stars around it are a unique laboratory for studying how relaxation processes lead to close interactions of stars and compact remnants with the central massive black hole, in particular those leading to the emission of gravitational waves. I review new results on the processes of strong mass segregation and loss-cone refilling by massive perturbers and resonant relaxation; describe observational evidence that these processes play a role in the Galactic Center and can be studied there; and discuss some of the implications for Extreme Mass Ration Inspiral event rates and their properties.

Recent observations of the Galactic center revealed a nuclear disk of young OB stars near the massive black hole (MBH), in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). In addition, observations show the existence of young B stars (the 'S-cluster') in an isotropic distribution in the close vicinity of the MBH ($<0.04$ pc). We use extended N-body simulations to probe the dynamical evolution of these two populations. We show that the stellar disk could have evolved to its currently observed state from a thin disk of stars formed in a gaseous disk, and that the dominant component in its evolution is the interaction with stars in the cusp around the MBH. We also show that the currently observed distribution of the S-stars could be consistent with a capture origin through 3-body binary-MBH interactions. In this scenario the stars are captured at highly eccentric orbits, but scattering by stellar black holes could change their eccentricity distribution to be consistent with current observations.

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.

We show that dynamical relaxation in the aftermath of a galactic merger and the ensuing formation and decay of a binary massive black hole (MBH), are dominated by massive perturbers (MPs) such as giant molecular clouds or clusters. MPs accelerate relaxation by orders of magnitude relative to 2-body stellar relaxation alone, and efficiently scatter stars into the binary MBH's orbit. The 3-body star-binary MBH interactions shrink the binary MBH to the point where energy losses from the emission of gravitational waves (GW) lead to rapid coalescence. We model this process based on observed and simulated MP distributions and take into account the decreased efficiency of the star-binary MBH interaction due to acceleration in the galactic potential. We show that mergers of gas-rich galactic nuclei lead to binary MBH coalescence well within the Hubble time. Moreover, lower-mass binary MBHs (<10^8 Msun) require only a few percent of the typical gas mass in a post-merger nucleus to coalesce in a Hubble time. The fate of a binary MBH in a gas poor galactic merger is less certain, although massive stellar structures (e.g. clusters, stellar rings) could likewise lead to efficient coalescence. These coalescence events are observable by their strong GW emission. MPs thus increase the cosmic rate of such GW events, lead to a higher mass deficit in the merged galactic core and suppress the formation of triple MBH systems and the resulting ejection of MBHs into intergalactic space.

Recent observations of the Galactic center revealed a nuclear disk of young OB stars, in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). Binaries in such nuclear disks, if they exist in non-negligible fractions, could have a major role in the evolution of the disks through binary heating of this stellar system. We suggest that interactions with/in binaries may explain some (or all) of the observed outlying young stars in the Galactic center. Such stars could have been formed in a disk, and later on kicked out from it through binary related interactions, similar to ejection of high velocity runaway OB stars in young clusters throughout the galaxy.

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.

The massive black hole (MBH) in the Galactic Center and the stars around it form a unique stellar dynamics laboratory for studying how relaxation processes affect the distribution of stars and compact remnants and lead to close interactions between them and the MBH. Recent theoretical studies suggest that processes beyond "minimal" two-body relaxation may operate and even dominate relaxation and its consequences in the Galactic Center. I describe loss-cone refilling by massive perturbers, strong mass segregation and resonant relaxation; review observational evidence that these processes play a role in the Galactic Center; and discuss some cosmic implications for the rates of gravitational wave emission events from compact remnants inspiraling into MBHs, and the coalescence timescales of binary MBHs.

We classify 329 late-type giants within 1 parsec of Sgr A*, using the adaptive optics integral field spectrometer SINFONI on the VLT. These observations represent the deepest spectroscopic data set so far obtained for the Galactic Center, reaching a 50% completeness threshold at the approximate magnitude of the helium-burning red clump (Ks ~ 15.5 mag.). Combining our spectroscopic results with NaCo H and Ks photometry, we construct an observed Hertzsprung-Russell diagram, which we quantitatively compare to theoretical distributions of various star formation histories of the inner Galaxy, using a chi-squared analysis. Our best-fit model corresponds to continuous star formation over the last 12 Gyr with a top-heavy initial mass function (IMF). The similarity of this IMF to the IMF observed for the most recent epoch of star formation is intriguing and perhaps suggests a connection between recent star formation and the stars formed throughout the history of the Galactic Center.

We present high-resolution seeing limited and AO NIR imaging observations of the stellar cluster within about one parsec of Sgr A*, the massive black hole at the centre of the Milky Way. Stellar number counts and the diffuse background light density were extracted from these observations in order to examine the structure of the nuclear stellar cluster.Our findings are as follows: (a) A broken-power law provides an excellent fit to the overall structure of the GC nuclear cluster. The power-law slope of the cusp is $\Gamma=0.19\pm0.05$, the break radius is $R_{\rm break} = 6.0'' \pm 1.0''$ or $0.22\pm0.04$ pc, and the cluster density decreases with a power-law index of $\Gamma=0.75\pm0.1$ outside of $R_{\rm break}$. (b) Using the best velocity dispersion measurements from the literature, we derive higher mass estimates for the central parsec than assumed until now. The inferred density of the cluster at the break radius is $2.8\pm1.3\times 10^{6} {\rm M_{\odot} pc^{-3}}$. This high density agrees well with the small extent and flat slope of the cusp. Possibly, the mass of the stars makes up only about 50% of the total cluster mass. (c) Possible indications of mass segregation in the cusp are found (d) The cluster appears not entirely homogeneous. Several density clumps are detected that are concentrated at projected distances of $R=3''$ and $R=7''$ from Sgr A*.(e) There appears to exist an under-density of horizontal branch/red clump stars near $R=5''$, or an over-density of stars of similar brightness at $R=3''$ and $R=7''$. (f) The extinction map in combination with cometary-like features in an L'-band image may provide support for the assumption of an outflow from Sgr A*.

We investigate the role of massive perturbers, such as giant molecular clouds or stellar clusters, in supplying low-angular momentum stars that pass very close to the central massive black hole (MBH) or fall into it. We show that massive pe rturbers can play an important role in supplying both binaries and single stars to the vicinity of the MBH. We discuss possible implications for the ejection of high velocity stars; for the capture of stars on tight orbits around the MBH; for the emission of gravitational waves from low-eccentricity inspiraling stars; and for the origin of the young main sequence B stars observed very near the Gal actic MBH. Massive perturbers may also enhance the the growth rate of MBHs, and may accelerate the dynamical orbital decay of coalescing binary MBHs.

High-precision radial-velocity techniques, which enabled the detection of extrasolar planets are now sensitive to relativistic effects in the data of spectroscopic binary stars (SBs). We show how these effects can be used to derive the absolute masses of the components of eclipsing single-lined SBs and double-lined SBs from Doppler measurements alone. High-precision stellar spectroscopy can thus substantially increase the number of measured stellar masses, thereby improving the mass-radius and mass-luminosity calibrations.

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.

The coherent torques between stars on orbits near massive black holes (MBHs) lead to resonant angular momentum relaxation. Due to the fact that orbits are Keplerian to good approximation, the torques efficiently change the magnitude of the angular momenta and rotate the orbital inclinations. As a result the stars are rapidly randomized. The galactic MBH is a good system for the observational study of resonant relaxation. The age of the young B-stars at a distance of $\sim0.01 pc$ from the MBH is comparable to the resonant relaxation time, implying that resonant relaxation may have played an important role in their dynamical structure. In contrast, the O-stars in the stellar disks at $\sim0.1 pc$ are younger than the resonant relaxation time, as required by their dynamical coherence. Resonant relaxation dynamics dominates the event rate of gravitational wave (GW) emission from inspiraling stars into MBHs of masses comparable to the Galactic MBH. Resonant relaxation leads to rates $\lesssim 10$ times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting these events by future GW detectors, such as LISA.

Resonant relaxation (RR) of orbital angular momenta occurs near massive black holes (MBHs) where the stellar orbits are nearly Keplerian and so do not precess significantly. The resulting coherent torques efficiently change the magnitude of the angular momenta and rotate the orbital inclination in all directions. As a result, many of the tightly bound stars very near the MBH are rapidly destroyed by falling into the MBH on low-angular momentum orbits, while the orbits of the remaining stars are efficiently randomized. We solve numerically the Fokker-Planck equation in energy for the steady state distribution of a single mass population with a RR sink term. We find that the steady state current of stars, which sustains the accelerated drainage close to the MBH, can be up to ~10 times larger than that due to non-coherent 2-body relaxation alone. RR mostly affects tightly bound stars, and so it increases only moderately the total tidal disruption rate, which is dominated by stars originating from less bound orbits farther away. We show that the event rate of gravitational wave (GW) emission from inspiraling stars, originating much closer to the MBH, is dominated by RR dynamics. The GW event rate depends on the uncertain efficiency of RR. The efficiency indicated by the few available simulations implies rates ~10 times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting such events by future GW detectors, such as LISA. However, a higher, but still plausible RR efficiency can lead to the drainage of all tightly bound stars and strong suppression of GW events from inspiraling stars. We apply our results to the Galactic MBH, and show that the observed dynamical properties of stars there are consistent with RR.

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.

We report the definite spectroscopic identification of 41 OB supergiants, giants and main sequence stars in the central parsec of the Galaxy. Detection of their absorption lines have become possible with the high spatial and spectral resolution and sensitivity of the adaptive optics integral field spectrometer SPIFFI/SINFONI on the ESO VLT. Several of these OB stars appear to be helium and nitrogen rich. Almost all of the ~80 massive stars now known in the central parsec (central arcsecond excluded) reside in one of two somewhat thick (<|h|/R>~0.14) rotating disks. These stellar disks have fairly sharp inner edges (R~1") and surface density profiles that scale as R^{-2}. We do not detect any OB stars outside the central 0.5 pc. The majority of the stars in the clockwise system appear to be on almost circular orbits, whereas most of those in the `counter-clockwise' disk appear to be on eccentric orbits. Based on its stellar surface density distribution and dynamics we propose that IRS 13E is an extremely dense cluster (core density > 3x10^8 sunmass/pc^3), which has formed in the counter-clockwise disk. The stellar contents of both systems are remarkably similar, indicating a common age of ~6+/-2 Myr. The K-band luminosity function of the massive stars suggests a top-heavy mass function and limits the total stellar mass contained in both disks to ~1.5x10^4 sunmass. Our data strongly favor in situ star formation from dense gas accretion disks for the two stellar disks. This conclusion is very clear for the clockwise disk and highly plausible for the counter-clockwise system.

Stars closely approaching the massive black hole in the center of the Galaxy provide a unique opportunity to probe post-Newtonian physics in a yet unexplored regime of celestial mechanics. Recent advances in infrared stellar spectroscopy allow the precise measurement of stellar Doppler shift curves and thereby the detection of beta-squared post-Newtonian effects (gravitational redshift in the black hole's potential and the transverse Doppler shift). We formulate a detection procedure in terms of a simplified post-Newtonian parametrization. We then use simulations to show that these effects can be decisively detected with existing instruments after about a decade of observations. We find that neglecting these effects can lead to statistically significant systematic errors in the derived black hole mass and distance.

We study the orbital parameters distribution of stars that are scattered into nearly radial orbits and then spiral into a massive black hole (MBH) due to dissipation, in particular by emission of gravitational waves (GW). This is important for GW detection, e.g. by the Laser Interferometer Space Antenna (LISA). Signal identification requires knowledge of the waveforms, which depend on the orbital parameters. We use analytical and Monte Carlo methods to analyze the interplay between GW dissipation and scattering in the presence of a mass sink during the transition from the initial scattering-dominated phase to the final dissipation-dominated phase of the inspiral. Our main results are (1) Stars typically enter the GW-emitting phase with high eccentricities. (2) The GW event rate per galaxy is a few per Gyr for typical central stellar cusps, almost independently of the relaxation time or the MBH mass. (3) For intermediate mass black holes (IBHs) of ~a thousand solar masses such as may exist in dense stellar clusters, the orbits are very eccentric and the inspiral is rapid, so the sources are very short-lived.

A massive black hole resides in the center of most, perhaps all galaxies. The one in the center of our home galaxy, the Milky Way, provides a uniquely accessible laboratory for studying in detail the connections and interactions between a massive black hole and the stellar system in which it grows; for investigating the effects of extreme density, velocity and tidal fields on stars; and for using stars to probe the central dark mass and probe post-Newtonian gravity in the weak- and strong-field limits. Recent results, open questions and future prospects are reviewed in the wider context of the theoretical framework and physical processes that underlie them. Contents: [1] Introduction (1.1) Astrophysical context (1.2) Science questions (1.3) Scope and connections to related topics [2] Observational overview: Stars in the Galactic center (2.1) The central 100 parsecs (2.2) The central parsec [3] Stellar dynamics at extreme densities (3.1) Physical processes and scales (3.2) The stellar cusp in the Galactic center (3.3) Mass segregation (3.4) Stellar Collisions [4] Probing the dark mass with stellar dynamics (4.1) Weighing and pinpointing the dark mass (4.2) Constraints on non-BH dark mass alternatives (4.3) Limits on MBH binarity (4.4) High-velocity runaway stars [5] Probing post-Newtonian gravity near the MBH (5.1) Relativistic orbital effects (5.2) Gravitational lensing [6] Strong star-MBH interactions (6.1) Tidal disruption (6.2) Dissipative interactions with the MBH [7] The riddle of the young stars (7.1) The difficulties of forming or importing stars near a MBH (7.2) Proposed solutions (7.3) Feeding the MBH with stellar winds [8] Outlook (8.1) Progress report (8.2) Future directions

We report 75 milli-arcsec resolution, near-IR imaging spectroscopy within the central 30 light days of the Galactic Center [...]. To a limiting magnitude of K~16, 9 of 10 stars in the central 0.4 arcsec, and 13 of 17 stars out to 0.7 arcsec from the central black hole have spectral properties of B0-B9, main sequence stars. [...] all brighter early type stars have normal rotation velocities, similar to solar neighborhood stars. We [...] derive improved 3d stellar orbits for six of these S-stars in the central 0.5 arcsec. Their orientations in space appear random. Their orbital planes are not co-aligned with those of the two disks of massive young stars 1-10 arcsec from SgrA*. We can thus exclude [...] that the S-stars as a group inhabit the inner regions of these disks. They also cannot have been located/formed in these disks [...]. [...] we conclude that the S-stars were most likely brought into the central light month by strong individual scattering events. The updated estimate of distance to the Galactic center from the S2 orbit fit is Ro = 7.62 +/- 0.32 kpc, resulting in a central mass value of 3.61 +/- 0.32 x 10^6 Msun. We happened to catch two smaller flaring events from SgrA* [...]. The 1.7-2.45 mum spectral energy distributions of these flares are fit by a featureless, red power law [...]. The observed spectral slope is in good agreement with synchrotron models in which the infrared emission comes from [...] radiative inefficient accretion flow in the central R~10 Rs region.

We propose a dynamical mechanism for capturing stars around a massive black hole (MBH), which is based on the accumulation there of a very dense cluster of compact stellar remnants. This study is motivated by the presence of ~10 young massive stars (M*~3-15 Mo, ~B9V--O8V) less than 0.04 pc from the MBH in the Galactic center (GC). Their existence in the extreme environment so close to a MBH is a challenge for theories of star formation and stellar dynamics. We show that young stars, which formed far from the MBH and were then scattered into eccentric orbits, repeatedly cross a cluster of stellar black holes (SBHs), where they may undergo rare direct 3-body exchanges with a MBH-SBH "binary". The interaction between two objects of similar mass ejects the SBH and captures the star on a tight orbit around the MBH. Such captures can naturally explain trends observed in the orbits of the young stars. We derive the capture cross-section, validate it by Monte-Carlo simulations, and calculate the number of captured stars in the GC using the currently uncertain estimates of the numbers of SBHs inside 0.04 pc and of young stars in the inner few pc of the GC. We find that under favorable conditions 3-body exchange can account for \~25% of the observed stars, mostly at the fainter end of the observed range. We discuss additional effects that possibly increase the capture efficiency. Future observations will establish whether there are enough SBHs and young stars for exchange captures to singly account for the central young stars. We estimate that there are also ~35 lower mass stars (M*~1-3 Mo, ~G2V--A0V) inside 0.04 pc similarly captured by exchanges with neutron stars. Ongoing replacement of compact remnants by main sequence stars may regulate the accumulation of compact remnants near the MBH.

The nature of ultraluminous X-ray sources (ULXs) is presently unknown. A possible explanation is that they are accreting intermediate mass black holes (IBHs) that are fed by Roche lobe overflow from a tidally captured stellar companion. We show that a star can circularize around an IBH without being destroyed by tidal heating (in contrast to the case of M_bh> 10^6 M_sun massive black holes in galactic centers, where survival is unlikely). We find that the capture and circularization rate is of the order of 5 \times 10^-8 yr^-1, almost independently of the cluster's relaxation time. We follow the luminosity evolution of the binary system during the main sequence Roche lobe overflow phase and show it can maintain ULX-like luminosities for >10 Myr. In particular, we show that the ULX in the young cluster MGG-11 in star-burst galaxy M82, which possibly harbors an IBH, is well explained by this mechanism, and we predict that \gtrsim 10% of similar clusters with IBHs have a tidally captured circularized star. The cluster can evaporate on a time-scale shorter than the lifetime of the binary. This raises the possibility of a ULX that outlives its host cluster, or even lights up only after the cluster has evaporated, in agreement with observations of host-less ULXs.

We present a study of the intrinsic X-ray and far-ultraviolet absorption in the Seyfert 1.5 galaxy Markarian 279 using simultaneous observations from the Chandra X-ray Observatory, the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope, and the Far Ultraviolet Spectroscopic Explorer (FUSE). We also present FUSE observations made at three additional epochs. We detect the Fe K-alpha emission line in the Chandra spectrum, and its flux is consistent with the low X-ray continuum flux level of Mrk 279 at the time of the observation. Due to low signal-to-noise ratios in the Chandra spectrum, no O VII or O VIII absorption features are observable in the Chandra data, but the UV spectra reveal strong and complex absorption from HI and high-ionization species such as O VI, N V, and C IV, as well as from low-ionization species such as C III, N III, C II, and N II in some velocity components. The far-UV spectral coverage of the FUSE data provides information on high-order Lyman series absorption, which we use to calculate the optical depths and line and continuum covering fractions in the intrinsic HI absorbing gas in a self-consistent fashion. The UV continuum flux of Mrk 279 decreases by a factor of ~7.5 over the time spanning these observations and we discuss the implications of the response of the absorption features to this change. From arguments based on the velocities, profile shapes, covering fractions and variability of the UV absorption, we conclude that some of the absorption components, particularly those showing prominent low-ionization lines, are likely associated with the host galaxy of Mrk 279, and possibly with its interaction with a close companion galaxy, while the remainder arises in a nuclear outflow.

Recent measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a 3.6-million-solar-mass black hole. Sgr A* is remarkably faint in all wavebands other than the radio region, however, which challenges current theories of matter accretion and radiation surrounding black holes. The black hole's rotation rate is not known, and therefore neither is the structure of space-time around it.Here we report high-resolution infrared observations of Sgr A* that reveal 'quiescent' emission and several flares. The infrared emission originates from within a few milliarcseconds of the black hole, and traces very energetic electrons or moderately hot gas within the innermost accretion region. Two flares exhibit a 17-minute quasi-periodic variability. If the periodicity arises from relativistic modulation of orbiting gas, the emission must come from just outside the event horizon, and the black hole must be rotating at about half of the maximum possible rate.

We report new astrometric and spectroscopic observations of the star S2 orbiting the massive black hole in the Galactic Center, which were taken at the ESO VLT with the adaptive optics assisted, near-IR camera NAOS/CONICA and the near-IR integral field spectrometer SPIFFI. We use these data to determine all orbital parameters of the star with high precision, including the Sun-Galactic Center distance, which is a key parameter for calibrating stellar standard candles and an important rung in the extragalactic distance ladder. Our deduced value of Ro = 8.0 +/- 0.4 kpc is the most accurate primary distance measurement to the center of the Milky Way and has minimal systematic uncertainties of astrophysical origin. It is in excellent agreement with other recent determinations of Ro.

We present proper motions for $>$40 stars at projected distances $\leq1.2''$ from Sagittarius A* (Sgr A*). We find evidence on a $\geq2\sigma$ level for radial anisotropy of the cluster of stars within $1''$ of Sgr A*. We find no evidence for a stationary source or variable source at the position of Sgr A*. We confirm/find accelerated motion for 6 stars, with 4 stars having passed the pericenter of their orbits during the observed time span. We calculated/constrained the orbital parameters of these stars. All orbits have moderate to high eccentricities. We discuss the possible bias in detecting preferentially orbits with high eccentricities. We find that the center of acceleration for all the orbits coincides with the radio position of Sgr A*. From the orbit of the star S2, the currently most tightly constrained one, we determine the mass of Sgr A* to be $3.3\pm0.7\times10^{6}$M$_{\odot}$ and its position to $2.0\pm2.4$ mas East and $2.7\pm4.5$ mas South of the nominal radio position. The mass estimate for the central dark mass from the orbit of S2 is fully consistent with the mass estimate of $3.4\pm0.5\times10^{6}$M$_{\odot}$ obtained from stellar proper motions within $1.2''$ of Sgr A* using a Leonard-Merritt mass estimator. We find that radio astronomical observations of the proper motion of Sgr A* in combination with its intrinsic source size place at the moment the tightest constraints on the mass density of Sgr A*, which must exceed $\rho_{\mathrm{Sgr A*}}>3\times10^{19}\mathrm{M}_{\odot}\mathrm{pc}^{-3}$.

We analyze deep near-IR adaptive optics imaging as well as new proper motion data of the nuclear star cluster of the Milky Way. The surface density distribution of faint stars peaks within 0.2" of the black hole candidate SgrA*. The radial density distribution of this stellar 'cusp' follows a power law of exponent 1.3-1.4. The K-band luminosity function of the overall nuclear stellar cluster (within 9" of SgrA*) resembles that of the large scale, Galactic bulge, but shows an excess of stars at K<14. We find that most of the massive early type stars at distances 1-10" from SgrA* are located in two rotating and geometrically thin disks. These disks are inclined at large angles and counter-rotate with respect to each other. Their stellar content is essentially the same, indicating that they formed at the same time. The star closest to SgrA* in 2002, S2, exhibits a 3.8 micron excess. We propose that the mid-IR emission either comes from the accretion flow around the black hole itself, or from dust in the accretion flow that is heated by the ultra-violet emission of S2.

We propose that there exists a class of transient sources, "squeezars", which are stars caught in highly eccentric orbits around a massive (m<10^8 Mo) black hole (MBH), whose atypically high luminosity (up to a significant fraction of their Eddington luminosity) is powered by tidal interactions with the MBH. Their existence follows from the presence of a mass sink, the MBH, in the galactic center, which drives a flow of stars into nearly radial orbits to replace those it has destroyed. We consider two limits for the stellar response to tidal heating: surface heating with radiative cooling ("hot squeezars") and bulk heating with adiabatic expansion ("cold squeezars"), and calculate the evolution of the squeezar orbit, size, luminosity and effective temperature. The squeezar formation rate is only ~0.05 that of tidal disruption flares, but squeezar lifetimes are many orders of magnitude longer, and so future observations of squeezars in nearby galaxies can probe the tidal process that feeds MBHs and the effects of extreme tides on stars. The mean number of squeezars orbiting the Galactic MBH is estimated at 0.1-1.

A massive black hole (MBH) in a galactic center drives a flow of stars into nearly radial orbits to replace those it destroyed. Stars whose orbits cross the event horizon r_s or the tidal disruption radius r_t are promptly destroyed in an orbital period P. Stars with orbital periapse r_p slightly larger than the sink radius q=max(r_s,r_t) may slowly spiral in due to dissipative interactions with the MBH, e.g. gravitational wave emission, tidal heating or accretion disk drag, with observable consequences and implications for the MBH growth rate. Unlike prompt destruction, the inspiral time is typically >>P. This time is limited by the same scattering process that initially deflected the star into its eccentric orbit, since it can deflect it again to a wider orbit where dissipation is inefficient. The ratio between slow and prompt event rates is therefore much smaller than that implied by the ratio of cross-sections, ~r_p/q, and so only prompt disruption contributes significantly to the mass of the MBH. Conversely, most stars that scatter off the MBH survive the extreme tidal interaction ("tidal scattering"). We derive general expressions for the inspiral event rate and the mean number of inspiraling stars, and show that the survival probability of tidally scattered stars is ~1, and that the number of tidally heated stars ("squeezars") and gravity wave emitting stars in the Galactic Center is ~0.1-1.

Close tidal interactions of stars with a central massive black hole (MBH) or with other stars in the high density cusp around it can affect a significant fraction of the stellar population within the MBH radius of influence. We consider three strong processes that have the potential of modifying stellar structure and evolution there. (1) Tidal spin-up by hyperbolic star-star encounters. (2) Tidal scattering of stars on the MBH. (3) Tidal heating of inspiraling stars - "squeezars" - that were tidally captured by the MBH. We discuss the implications for stellar populations near MBHs and for the growth of MBHs by tidal disruption of stars, and the possible observational signature of such processes near the MBH in the Galactic Center. We compare the event rates of prompt tidal encounters (tidal disruption and tidal scattering) and slow inspiral events (squeezars / tidal capture), and find that tidal capture is at least an order of magnitude less efficient than prompt disruption. This means that past studies, which assigned similar weights to prompt disruption and tidal capture, over-estimated the contribution of tidal disruption to the growth of the MBH by at least a factor of two.

We report on the first results obtained using adaptive optics measurements of the Galactic Centre done with NAOS/CONICA.

Measurements of stellar velocities and variable X-ray emission near the centre of the Milky Way have provided the strongest evidence so far that the dark mass concentrations seen in many galactic nuclei are likely supermassive black holes, but have not yet excluded several alternative configurations. Here we report ten years of high resolution astrometric imaging that allow us to trace two thirds of the orbit of the star currently closest to the compact radio source and massive black hole candidate SgrA*. In particular, we have observed both peri- and apocentre passages. Our observations show that the star is on a bound, highly elliptical Keplerian orbit around SgrA*, with an orbital period of 15.2 years and a peri-centre distance of only 17 light hours. The orbital elements require an enclosed point mass of 3.7+-1.5x10^6 solar masses. The data exclude with high confidence that the central dark mass consists of a cluster of astrophysical objects or massive, degenerate fermions, and strongly constrain the central density structure.

This is a pedagogical review of recent results on the interactions of central massive black holes with stars very near them, focused on the black hole in the center of the Milky Way. Table of contents: [1] Introduction [2] Stellar dynamics near a black hole [2.1] Physical scales [2.2] A relaxed stellar system around a MBH [3] The stellar collider in the Galactic Center [3.1] The case for a dense stellar cusp in the Galactic Center [3.2] Tidal spin-up [3.3] Tidal scattering [4] The gravitational telescope in the Galactic Center [4.1] Gravitational lensing by a point mass [4.2] Pinpointing the MBH with lensed images [4.3] The detection of gravitational lensing [4.4] Magnification bias [4.5] Beyond the point mass lens approximation [5] Summary

Some of the mass that feeds the growth of a massive black hole (BH) in a galactic center is supplied by tidal disruption of stars that approach it on unbound, low angular momentum orbits. For each star that is disrupted, others narrowly escape after being subjected to extreme tidal distortion, spin-up, mixing and mass-loss, which may affect their evolution and appearance. We show that it is likely that a significant fraction of the stars around massive BHs in galactic centers have undergone such extreme tidal interactions and survived subsequent total disruption, either by being deflected off their orbit or by missing the BH due to its Brownian motion. We discuss possible long-term observable consequences of this process, which may be relevant for understanding the nature of stars in galactic centers, and may provide a signature of the existence of massive BHs there.

We reassess constraints on the cosmological baryon density from observations of the mean decrement and power spectrum of the Lyman-alpha forest, taking into account uncertainties in all free parameters in the simplest gravitational instability model. The uncertainty is dominated by that of the photoionizing background, but incomplete knowledge of the thermal state of the intergalactic medium also contributes significantly to the error-budget. While current estimates of the baryon fraction from the forest do prefer values that are somewhat higher than the big bang nucleosynthesis value of \Omega_b h^2 = 0.02 \pm 0.001, the discrepancy is at best about 3 \sigma. For instance, assuming the highest estimate of the ionizing background, as indicated by recent measurements of a large escape fraction from Lyman-break galaxies by Steidel, Pettini & Adelberger, we find \Omega_b h^2 = 0.045 \pm 0.008. A recent measurement of the ionizing background from the proximity effect by Scott et al., on the other hand, implies \Omega_b h^2 = 0.03 \pm 0.01. We provide an expression from which future likelihoods for \Omega_b h^2 can be derived as measurements of the ionizing background improve -- consistency among constraints from the forest, nucleosynthesis and the microwave background will provide a powerful test of the gravitational instability model. We also develop a formalism which treats lower bounds on the baryon density in a statistical manner, which is appropriate if only a lower bound on the ionizing background is known. Finally, we discuss the implications of the escape fraction measurement for the age, structure and stellar content of Lyman-break galaxies.

The effect of stars on the lensing properties of the supermassive black hole in the Galactic Center is similar to the effect of planets on microlensing by a star. We show that the dense stellar cluster around SgrA* increases by factors of a few the probability of high-magnification lensing events of a distant background source by the black hole. Conversely, the gravitational shear of the black hole changes and enhances the microlensing properties of the individual stars. The effect is largest when the source image lies near the Einstein radius of the black hole 1.75"+/-0.20" for a source at infinity). We estimate that the probability of observing at least one distant background star which is magnified by a factor >5 in any infrared snapshot of the inner ~ 2" of the Galactic Center is 1% with a K-band detection threshold of 20 mag. The largest source of uncertainty in this estimate is the luminosity function of the background stars. The gravitational shear of the black hole lengthens the duration of high-magnification events near the Einstein radius up to a few months, and introduces a large variety of lightcurve shapes that are different from those of isolated microlenses. Identification of such events by image subtraction can be used to probe the mass function, density and velocity distributions of faint stars near the black hole, which are not detectable otherwise.

A new statistical method for pinpointing the massive black hole (BH) in the Galactic Center on the IR grid is presented and applied to astrometric IR observations of stars close to the BH. This is of interest for measuring the IR emission from the BH, in order to constrain accretion models; for solving the orbits of stars near the BH, in order to measure the BH mass and to search for general relativistic effects; and for detecting the fluctuations of the BH away from the dynamical center of the stellar cluster, in order to study the stellar potential. The BH lies on the line connecting the two images of any background source it gravitationally lenses, and so the intersection of these lines fixes its position. A combined search for a lensing signal and for the BH shows that the most likely point of intersection coincides with the center of acceleration of stars orbiting the BH. This statistical detection of lensing by the BH has a random probability of ~0.01. It can be verified by deep IR stellar spectroscopy, which will determine whether the most likely lensed image pair candidates (listed here) have identical spectra.

We show that main-sequence stars in dense stellar cusps around massive black holes are likely to rotate at a significant fraction of the centrifugal breakup velocity due to spin-up by hyperbolic tidal encounters. We use realistic stellar structure models to calculate analytically the tidal spin-up in soft encounters, and extrapolate these results to close and penetrating collisions using smoothed particle hydrodynamics simulations. We find that the spin-up falls off only slowly with distance from the black hole because the increased tidal coupling in slower collisions at larger distances compensates for the decrease in the stellar density. We apply our results to the stars near the massive black hole in the Galactic Center. Over their lifetime, ~1 Msol main sequence stars in the inner 0.3 pc of the Galactic Center are spun-up on average to ~10%--30% of the centrifugal breakup limit. Such rotation is ~20--60 times higher than is usual for such stars and may affect their subsequent evolution and their observed properties.

The ISO-SWS 2.5-45 um infrared spectroscopic observations of the nucleus of the Seyfert 2 galaxy NGC 1068 (see companion paper) are combined with a compilation of UV to IR narrow emission line data to determine the spectral energy distribution (SED) of the obscured extreme-UV continuum that photoionizes the narrow line emitting gas in the active galactic nucleus. We search a large grid of gas cloud models and SEDs for the combination that best reproduces the observed line fluxes and NLR geometry. Our best fit model reproduces the observed line fluxes to better than a factor of 2 on average and is in general agreement with the observed NLR geometry. It has two gas components that are consistent with a clumpy distribution of dense outflowing gas in the center and a more extended distribution of less dense and more clumpy gas farther out that has no net outflow. The best fit SED has a deep trough at ~4 Ryd, which is consistent with an intrinsic Big Blue Bump that is partially absorbed by ~6x10^19 cm^-2 of neutral hydrogen interior to the NLR.

We present ISO-SWS spectroscopy of NGC 1068 for the wavelength range 2.4 to 45um, detecting a total of 36 emission lines. Most of the observed transitions are fine structure and recombination lines originating in the narrow line region. We compare the line profiles of optical lines and reddening-insensitive infrared lines to constrain the dynamical structure and extinction properties of the NLR. The considerable differences found are most likely explained by two effects. (1) The spatial structure of the NLR is a combination of a highly ionized outflow cone and lower excitation extended emission. (2) Parts of the NLR, mainly in the receding part at velocities above systemic, are subject to extinction that is significantly suppressing optical emission. Line asymmetries and net blueshifts remain, however, even for infrared fine structure lines suffering very little obscuration. This may be either due to an intrinsic asymmetry of the NLR, or due to a very high column density obscuring component which is hiding part of the NLR even from infrared view. Mid-infrared emission of molecular hydrogen in NGC 1068 arises in a dense molecular medium at temperatures of a few hundred Kelvin that is most likely closely related to the warm and dense components seen in the near-infrared H2 transitions, and in millimeter wave tracers of molecular gas. Any emission of the putative pc-scale molecular torus is likely overwhelmed by this larger scale emission.

We report infrared spectroscopy of the prototypical Seyfert 2 galaxy NGC 1068, aiming at detection of broad components of hydrogen recombination lines that originate in the obscured broad-line region. Using the Short Wavelength Spectrometer on board the Infrared Space Observatory, we have observed for the first time the regions of Brackett beta 2.626um and Pfund alpha 7.460um, and present improved data for Brackett alpha 4.052um. No significant broad components are detected, implying an equivalent visual extinction to the broad-line region of at least 50 magnitudes and an obscuring column density of at least 10^23 cm^-2. While consistent with a highly obscured broad-line region, as required by the classical unified scenario, these limits are not yet significant enough to discriminate strongly between different torus models or to constrain properties of the gas causing the very large X-ray obscuration. We discuss the systematic limitations of infrared broad-line region searches and suggest that Brackett alpha may often be the most favorable transition for future searches.

We analyze three sets of infrared star counts in the inner ~0.5 pc of the Galactic Center. We perform statistical tests on the star counts and model in detail the extinction field and the effects of dwarf-giant collisions on the luminosity function. We find that both the star counts and the depletion of the brightest stars in the inner ~0.05 pc can be explained by a r^-\alpha stellar cusp with \alpha in the range 3/2 to 7/4, in which the envelopes of the brightest giants are destroyed by stellar collisions. Such a cusp is consistent with the Bahcall-Wolf solution for the distribution of stars that have undergone two-body relaxation around a black hole. We show that systematic uncertainties due to variable extinction and unrelaxed stars are probably small, but deeper star counts are required to confirm these results. We estimate that the tidal disruption rate of cusp stars by the black hole is few x 10^-5 yr^-1.

We investigate microlensing amplification of faint stars in the dense stellar cluster in the Galactic Center (GC) by the super-massive black hole (BH). Such events would appear very close to the position of the radio source SgrA*, which is thought to coincide with the BH, and could be observed during the monitoring of stellar motions in the GC. We use the observed K-band (2.2 um) luminosity function (KLF) in the GC and in Baade's Window, as well as stellar population synthesis computations, to construct KLF models for the inner 300 pc of the Galaxy. These, and the observed dynamical properties of this region, are used to compute the rates of microlensing events, which amplify stars above specified detection thresholds. We present computations of the lensing rates and amplifications as functions of the event durations (weeks to years), for a range of detection thresholds. We find that short events dominate the total rate and that long events tend to have large amplifications. For the current detection limit of K=17 mag, the total microlensing rate is 0.003 1/yr, and the rate of events with durations >1 yr is 0.001 1/yr. Recent GC proper motion studies have revealed the possible presence of one or two variable K-band sources very close to SgrA* (Genzel et al 97; Ghez et al 98). These sources may have attained peak brightnesses of K~15 mag, about 1.5-2 mag above the observational detection limits, and appear to have varied on a timescale of ~1 yr. This behavior is consistent with long-duration microlensing of faint stars by the BH. However, we estimate that the probability that such an event could have been detected during the course of the recent proper motion studies is \~0.5%. A ten-fold improvement in the detection limit and 10 yr of monthly monitoring would increase the total detection probability to ~20%. (Abridged)

We present ISO-SWS and ISOPHOT-S spectroscopy of the Seyfert galaxy NGC 4151. We detect a total of 17 fine structure emission lines emitted by a wide range of low- and high-excitation ions, two rotational lines of molecular hydrogen, and the Br beta HI line. We find that the mid-IR fine structure line profiles display blue asymmetries which are very similar to those observed in the optical lines produced in the narrow line region. Because the mid-infrared lines are much less sensitive to extinction than are the optical lines this similarity places strong constraints on scenarios which have been invoked to explain the optical line asymmetries. For example, we are able to rule out the simplest radial-motion-plus-dust scenarios for the production of the line asymmetries. Our preferred model is that of a central, geometrically thin but optically thick, obscuring screen of sub-arcsecond extension, enclosing a total hydrogen gas mass of >~5E6 Msun. This mass may be molecular. The weakness of `PAH' emission features in the low resolution spectrum is evidence that star formation plays a minor role in the circumnuclear region of NGC 4151. In a companion paper (Alexander et al., accepted by ApJ), we use the rich set of mid-infrared lines to determine the obscured photoionizing continuum produced by the active galactic nucleus.

The ISO-SWS infrared spectroscopic observations of the nucleus of Seyfert galaxy NGC4151, which are described in a companion paper, are used together with a compilation of UV to IR narrow emission line data to determine the spectral shape of the obscured extreme-UV continuum that photoionizes the narrow line emitting gas in the active galactic nucleus. We present a new method to determine the best fitting photoionizing continuum and emission line cloud model from a heterogeneous set of emission line data. For NGC4151, we find a best-fit model that reproduces the observed line fluxes to within a factor of 2 on average, and which is consistent with the observed geometry of the optical Narrow Line Region (NLR). Our model consists of a clumpy, optically thick (ionization bounded) gas distribution, with a hydrogen gas density of ~1000 cm^-3 and a volume filling factor of 6.5*10^-4. Our best fitting spectral energy distribution (SED) falls sharply beyond the Lyman limit and then rises sharply again towards 100 eV. In particular, it does not display a `Big Blue Bump' signature of a hot accretion disk. We find that this SED, which best reproduces the NLR line emission on the 100--500 pc scale, does not have enough UV photons to produce the observed BLR recombination emission from the inner 1 pc. This suggests that the BLR is photoionized by the intrinsic continuum source, which does have a strong UV component (perhaps a Big Blue Bump), but that this UV component is absorbed by material located between the NLR and BLR. Our analysis suggests that the absorber consists of \~5*10^19 cm^-2 of neutral hydrogen. Such an absorber was independently detected by UV absorption lines (Kriss et al. 1992, 1995).

We describe a model-independent method of assessing the uncertainties in cross-correlation lags determined from AGN light curves, and use this method to investigate the reality of lags between UV and optical continuum variations in well-studied AGNs. Our results confirm the existence of such lags in NGC 7469. We find that the continuum variations at 1825 A, 4845 A, and 6962 A follow those at 1315A by 0.22^{+0.12}_{-0.13} days, 1.25^{+0.48}_{-0.35} days, and 1.84^{+0.93}_{-0.94} days, respectively, based on the centroids of the cross-correlation functions; the error intervals quoted correspond to 68% confidence levels, and each of these lags is greater than zero at no less than 97% confidence. We do not find statistically significant interband continuum lags in NGC 5548, NGC 3783, or Fairall 9. Wavelength-dependent continuum lags may be marginally detected in the case of NGC 4151. However, on the basis of theoretical considerations, wavelength-dependent continuum lags in sources other than NGC 7469 are not expected to have been detectable in previous experiments. We also confirm the existence of a statistically significant lag between X-ray and UV continuum variations in the blazar PKS 2155-304.

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.

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.

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 $\Lambda = 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.