We investigate the evolution of supermassive binary black holes (BBHs) in galaxies with realistic property distributions and the gravitational wave (GW) radiation from the cosmic population of these BBHs. We incorporate a comprehensive treatment of the BBH dynamical interactions with their environments by including the effects of galaxy triaxial shapes and inner stellar distributions, and generate a large number of BBH evolution tracks. By combining these BBH evolution tracks, galaxy mass functions, galaxy merger rates, and supermassive black hole-host galaxy relations into our model, we obtain the statistical distributions of surviving BBHs, BBH coalescence rates, the strength of their GW radiation, and the stochastic GW background contributed by the cosmic BBH population. About ~1%-3% (or ~10%) of supermassive BHs at nearby galactic centers are expected to be binaries with mass ratio >1/ (or >1/100). The characteristic strain amplitude of the GW background at frequency 1/yr are estimated to be ~$2.0^{+1.4}_{-0.8}\times 10^{-16}$, and its upper bound of the results obtained with the different BH--host galaxy relations can be up to $5.4\times 10^{-16}$, which wait to be tested by future experiments (e.g., SKA, FAST, ngVLA). The turnover frequency of the GW background spectrum is at ~0.25nHz. The uncertainties to the above estimates and prospects for detecting individual sources are also discussed. The application of the cosmic BBH population to the LISA band provides a lower limit to the detection rate of BBHs by LISA, ~0.9/yr.

We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope.

We developed a pulsar search pipeline based on PRESTO (PulsaR Exploration and Search Toolkit). This pipeline simply runs dedispersion, FFT (Fast Fourier Transformation), and acceleration search in process-level parallel to shorten the processing time. With two parallel strategies, the pipeline can highly shorten the processing time in both the normal searches or acceleration searches. This pipeline was first tested with PMPS (Parkes Multibeam Pulsar Survery) data and discovered two new faint pulsars. Then, it was successfully used in processing the FAST (Five-hundred-meter Aperture Spherical radio Telescope) drift scan data with tens of new pulsar discoveries up to now. The pipeline is only CPU-based and can be easily and quickly deployed in computing nodes for testing purposes or data processes.

Observations and semianalytical galaxy formation and evolution models (SAMs) have suggested the existence of a stellar mass-stellar metallicity relation (MZR), which is shown to be universal for different types of galaxies over a large range of stellar masses ($M_*\sim 10^3$-$10^{11}M_\odot$) and dark matter (DM) halo masses ($M_{\rm halo}\sim 10^9$-$10^{15}h^{-1}M_\odot$). In this work, we construct a chemical evolution model to investigate the origin of the MZR, including both the effects of gas inflows and outflows in galaxies. We solve the MZR from the chemical evolution model, by assuming that the cold gas mass ($M_{\rm cold}$) and the stellar feedback efficiency ($\beta$) follow some power-law scaling relationships with $M_*$ during the growth of a galaxy, i.e., $M_{\rm cold}\propto M_*^{\alpha_{\rm gs}}$ and $\beta\propto M_*^{\alpha_{\beta{\rm s}}}$. We use the SAM to obtain these power-law scaling relations, which appear to be roughly universal over a large range of stellar masses for both satellites and central galaxies within a large range of halo masses. The range of the MZRs produced by our models is in a narrow space, which provides support to the universality of the MZRs. The formation of the MZR is a result caused jointly by that the cold gas fraction decreases with increasing $M_*$ and by that the stellar feedback efficiency decreases with increasing $M_*$ in the galaxy growth, and the exponent in the MZR is around $-\alpha_{\beta{\rm s}}$ or $1-\alpha_{\rm gs}$. The MZR represents an ''average'' evolutional track for the stellar metallicity of a galaxy. The comparison of our model with some previous models for the origin of MZRs is also discussed.

Observations reveal a universal stellar mass-stellar metallicity relation (MZR) existing in Local Group dwarfs of different types, $Z_*\propto M_*^{\alpha}$ with $\alpha=0.30\pm0.02$. In this work, we investigate the "universality" of the MZRs for both satellites and central galaxies in a large number of different host dark matter halos covering a large mass range of $10^9$-$10^{15}h^{-1}M_\odot$, by using a semianalytical galaxy formation and evolution model. We obtain the following results. (1) The exponents ($\alpha$) for the MZRs of the satellites in halos with the same mass as the Milky Way halo but different individual assembly histories are mostly $\sim$0.2-0.4, i.e., having a scatter of $\sim 0.2$; and the scatter of $\alpha$ increases with decreasing halo masses. (2) The MZR relations are changed little by the variation of halo masses and the classification between central galaxies and satellites, if many halos with the same mass are stacked together. (3) A double power law exists in the MZR relations for both central galaxies and stacked satellites, with $\alpha\sim$0.2-0.4 at $10^3M_\odot< M_*<10^{8} M_\odot$ and a relatively higher $\alpha\sim0.5$ at $10^8M_\odot<M_*<10^{11}M_\odot$. (4) The high-mass satellites ($M_*>10^8M_\odot$) existing mostly in high-mass halos can lead to an apparent increase of $\alpha$ (from $\sim0.2$ to $\sim0.4$) with increasing host halo masses shown in the single power law fitting results of stacked satellites. The universality of the MZR suggests the common physical processes in stellar formation and chemical evolution of galaxies can be unified over a large range of galaxy masses and halo masses.

The massive galaxy cluster "El Gordo" (ACT-CL J0102--4915) is a rare merging system with a high collision speed suggested by multi-wavelength observations and the theoretical modeling. Zhang et al. (2015) propose two types of mergers, a nearly head-on merger and an off-axis merger with a large impact parameter, to reproduce most of the observational features of the cluster, by using numerical simulations. The different merger configurations of the two models result in different gas motion in the simulated clusters. In this paper, we predict the kinetic Sunyaev-Zel'dovich (kSZ) effect, the relativistic correction of the thermal Sunyaev-Zel'dovich (tSZ) effect, and the X-ray spectrum of this cluster, based on the two proposed models. We find that (1) the amplitudes of the kSZ effect resulting from the two models are both on the order of $\Delta T/T\sim10^{-5}$; but their morphologies are different, which trace the different line-of-sight velocity distributions of the systems; (2) the relativistic correction of the tSZ effect around $240 {\rm\,GHz}$ can be possibly used to constrain the temperature of the hot electrons heated by the shocks; and (3) the shift between the X-ray spectral lines emitted from different regions of the cluster can be significantly different in the two models. The shift and the line broadening can be up to $\sim 25{\rm\,eV}$ and $50{\rm\,eV}$, respectively. We expect that future observations of the kSZ effect and the X-ray spectral lines (e.g., by ALMA, XARM) will provide a strong constraint on the gas motion and the merger configuration of ACT-CL J0102--4915.

X-ray variation is a ubiquitous feature of active galactic nuclei (AGNs), however, its origin is not well understood. In this paper, we show that the X-ray flux variations in some AGNs, and correspondingly the power spectral densities (PSDs) of the variations, may be interpreted as being caused by absorptions of eclipsing clouds or clumps in the broad line region (BLR) and the dusty torus. By performing Monte-Carlo simulations for a number of plausible cloud models, we systematically investigate the statistics of the X-ray variations resulting from the cloud eclipsing and the PSDs of the variations. For these models, we show that the number of eclipsing events can be significant and the absorption column densities due to those eclipsing clouds can be in the range from 10^{21} to 10^{24} cm^{-2}, leading to significant X-ray variations. We find that the PSDs obtained from the mock observations for the X-ray flux and the absorption column density resulting from these models can be described by a broken double power law, similar to those directly measured from observations of some AGNs. The shape of the PSDs depend strongly on the kinematic structures and the intrinsic properties of the clouds in AGNs. We demonstrate that the X-ray eclipsing model can naturally lead to a strong correlation between the break frequencies (and correspondingly the break timescales) of the PSDs and the masses of the massive black holes (MBHs) in the model AGNs, which can be well consistent with the one obtained from observations. Future studies of the PSDs of the AGN X-ray (and possibly also the optical-UV) flux and column density variations may provide a powerful tool to constrain the structure of the BLR and the torus and to estimate the MBH masses in AGNs.

In this paper, we investigate the quality of constraining the spin of the massive black hole (MBH) at the Galactic center (GC) by using full general relativistic simulations of the motion of a surrounding star. We obtain the dependence mapping of the spin-induced signals on any spin direction of the MBH for given example stars, which indicates the feasibility to test whether the spin direction is the same as the normal of the young stellar disk located at the GC, and, further to provide insights into the assembly history of the MBH. We demonstrate the quality of constraining the MBH spin that may be achieved, given any set of the astrometric and the redshift precisions of observational facilities. We find that in the ranges of the astrometric and the velocity precisions with 1--30$\mu$as and 0.1--10 km/s, an improvement in astrometric precision would be more effective at improving the quality of constraining the spin than an improvement in velocity precision. We obtain the parameter space of the semimajor axis and the eccentricity for the orbit of the target star that a high-precision constraint on the GC MBH spin can be obtained via the motion of the star. Our results show that the spin of the GC MBH can be constrained with 1-sigma error <~0.1 or even >~0.02 by monitoring the orbital motion of a star, if existing as expected, with semimajor axis <~300AU and eccentricity >~0.95 over a period shorter than a decade through future facilities.

Accurate estimation of the merger timescale of galaxy clusters is important to understand the cluster merger process and further the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases with increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions 0.15, compared with the timescale obtained with zero gas fractions. The baryonic effect is significant for a wide range of merger parameters and especially more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have impacts on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of "no-gas" may exist in the results obtained from the dark matter-only cosmological simulations.

The observational features of the massive galaxy cluster "El Gordo" (ACT-CL J0102-4915), such as the X-ray emission, the Sunyaev-Zel'dovich (SZ) effect, and the surface mass density distribution, indicate that they are caused by an exceptional ongoing high-speed collision of two galaxy clusters, similar to the well-known Bullet Cluster. We perform a series of hydrodynamical simulations to investigate the merging scenario and identify the initial conditions for the collision in ACT-CL J0102-4915. By surveying the parameter space of the various physical quantities that describe the two colliding clusters, including their total mass (M), mass ratio (\xi), gas fractions (f_b), initial relative velocity (V), and impact parameter (P), we find out an off-axis merger with P~800h_{70}^{-1}kpc, V~2500km/s, M~3x10^{15}Msun, and \xi=3.6 that can lead to most of the main observational features of ACT-CL J0102-4915. Those features include the morphology of the X-ray emission with a remarkable wake-like substructure trailing after the secondary cluster, the X-ray luminosity and the temperature distributions, and also the SZ temperature decrement. The initial relative velocity required for the merger is extremely high and rare compared to that inferred from currently available Lambda cold dark matter (LCDM) cosmological simulations, which raises a potential challenge to the LCDM model, in addition to the case of the Bullet Cluster.

Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor axis of ~590AU and an orbital period of ~1.2 year. The flat optical and UV continua are mainly emitted from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed sharp drop off and flux deficit at wavelength lambda ~ 4000-2500 Angstrom is due to a gap (or hole) opened by the secondary BH migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among the spectra of quasars.

The S-stars in the Galactic center (GC) are anticipated to provide unique dynamical constraint on the spin of the GC massive black hole (MBH). In this paper, we develop a fast full general relativistic method to simultaneously constrain the MBH mass, spin, and spin direction by considering both the motion of a star and the propagation of photons from the star to a distant observer. Assuming some example stars, we demonstrate that the spin-induced effects on the projected trajectory and redshift curve of a star depend on both the value and the direction of the spin. The maximum effects over a full orbit can differ by a factor upto more than one order of magnitude for cases with significantly different spin directions. Adopting the Markov Chain Monte Carlo fitting technique, we illustrate that the spin of the GC MBH is likely to be well constrained by using the motion of S0-2/S2 over a period of ~45yr if it is close to one and the astrometric and spectroscopic precisions (sigma_p,sigma_Z) can be as high as (10muas, 1km/s). In the mean time, the distance from the sun to the GC and the MBH mass can also be constrained to an unprecedented accuracy (0.01%-0.1%). If there exists a star with semimajor axis significantly smaller than that of S0-2/S2 and eccentricity larger than that of S0-2/S2, the MBH spin can be constrained with high accuracy over a period of <~10yr for (sigma_p,sigma_Z) ~ (10muas,1km/s), even if the spin is only moderately large (>~0.2).

The SkyMapper Southern Sky Survey is carrying out a search for the most metal-poor stars in the Galaxy. It identifies candidates by way of its unique filter set that allows for estimation of stellar atmospheric parameters. The set includes a narrow filter centered on the Ca II K 3933A line, enabling a robust estimate of stellar metallicity. Promising candidates are then confirmed with spectroscopy. We present the analysis of Magellan-MIKE high-resolution spectroscopy of 122 metal-poor stars found by SkyMapper in the first two years of commissioning observations. 41 stars have [Fe/H] <= -3.0. Nine have [Fe/H] <= -3.5, with three at [Fe/H] ~ -4. A 1D LTE abundance analysis of the elements Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, Ba and Eu shows these stars have [X/Fe] ratios typical of other halo stars. One star with low [X/Fe] values appears to be "Fe-enhanced," while another star has an extremely large [Sr/Ba] ratio: >2. Only one other star is known to have a comparable value. Seven stars are "CEMP-no" stars ([C/Fe] > 0.7, [Ba/Fe] < 0). 21 stars exhibit mild r-process element enhancements (0.3 <=[Eu/Fe] < 1.0), while four stars have [Eu/Fe] >= 1.0. These results demonstrate the ability to identify extremely metal-poor stars from SkyMapper photometry, pointing to increased sample sizes and a better characterization of the metal-poor tail of the halo metallicity distribution function in the future.

Observations reveal that the peaks of the X-ray map and the Sunyaev-Zel'dovich (SZ) effect map of some galaxy clusters are offset from each other. In this paper, we perform a set of hydrodynamical simulations of mergers of two galaxy clusters to investigate the spatial offset between the maxima of the X-ray and the SZ surface brightness of the merging clusters. We find that significantly large SZ-X-ray offsets (>100kpc) can be produced during the major mergers of galaxy clusters. The significantly large offsets are mainly caused by a `jump effect' occurred between the primary and secondary pericentric passages of the two merging clusters, during which the X-ray peak may jump to the densest gas region located near the center of the small cluster, but the SZ peak remains near the center of the large one. Our simulations show that merging systems with higher masses and larger initial relative velocities may result in larger offset sizes and longer offset time durations; and only nearly head-on mergers are likely to produce significantly large offsets. We further investigate the statistical distribution of the SZ-X-ray offset sizes and find that (1) the number distribution of the offset sizes is bimodal with one peak located at low offsets ~0 and the other at large offsets ~350-450kpc/h, but the objects with intermediate offsets are scarce; and (2) the probabilities of the clusters in the mass range higher than 2x10^{14}Msun/h that have offsets larger than 20, 50, 200, 300, and 500kpc/h are 34.0%, 11.1%, 8.0%, 6.5%, and 2.0% respectively at z=0.7. The probability is sensitive to the underlying pairwise velocity distribution and the merger rate of clusters. Future observations on the offsets for a large number of clusters may put strong constraints on the cosmic velocity fields on the cluster scale and the cluster merger rate. (Abridged)

In this paper, we study the chemical properties of the stars in the dwarf satellites around the MW-like host galaxies, and explore the possible effects of several baryonic processes, including supernova (SN) feedback, the reionization of the universe and H$_2$ cooling, on them and how current and future observations may put some constraints on these processes. We use a semi-analytical model to generate MW-like galaxies, for which a fiducial model can reproduce the luminosity function and the stellar metallicity--stellar mass correlation of the MW dwarfs. Using the simulated MW-like galaxies, we focus on investigating three metallicity properties of their dwarfs: the stellar metallicity--stellar mass correlation of the dwarf population, and the metal-poor and metal-rich tails of the stellar metallicity distribution in individual dwarfs. We find that (1) the slope of the stellar metallicity--stellar mass correlation is sensitive to the SN feedback strength and the reionization epoch; (2) the extension of the metal-rich tails is mainly sensitive to the SN feedback strength; (3) the extension of the metal-poor tails is mainly sensitive to the reionization epoch; (4) none of the three chemical properties are sensitive to the H$_2$ cooling process; and (5) comparison of our model results with the current observational slope of the stellar metallicity--stellar mass relation suggests that the local universe is reionized earlier than the cosmic average and local sources may have a significant contribution to the reionization in the local region, and an intermediate to strong SN feedback strength is preferred. Future observations of metal-rich and metal-poor tails of stellar metallicity distributions will put further constraints on the SN feedback and the reionization processes.

We present a detailed analysis of high-resolution, high S/N spectra for 5 Aquarius stream stars observed with the MIKE spectrograph on the Magellan Clay telescope. Our sample represents one third of the 15 known members in the stream. We find the stream is not mono-metallic: the metallicity ranges from [Fe/H] = -0.63 to -1.58. No anti-correlation in Na-O abundances is present, and we find a strong positive Mg-Al relationship, similar to that observed in the thick disk. We find no evidence that the stream is a result of a disrupted classical globular cluster, contrary to a previously published claim. High [(Na, Ni, alpha)/Fe] and low [Ba/Y] abundance ratios in the stream suggests it is not a tidal tail from a disrupted dwarf galaxy, either. The stream is chemically indistinguishable from Milky Way field stars with the exception of one candidate, C222531-145437. From its position, velocity, and detailed chemical abundances, C222531-145437 is likely a star that was tidally disrupted from omega-Centauri. We propose the Aquarius stream is Galactic in origin, and could be the result from a disk-satellite perturbation in the Milky Way thick disk on the order of a few Gyr ago: derived orbits, UVW velocities, and angular momenta of the Aquarius members offer qualitative support for our hypothesis. Assuming C222531-145437 is a tidally disrupted member of omega-Centauri, this system is the most likely disk perturber. In the absence of compelling chemical and/or dynamical evidence that the Aquarius stream is the tidal tail of a disrupted satellite, we advocate the "Aquarius group" as a more appropriate description. Like the Canis Major over-density, as well as the Hercules and Monoceros groups, the Aquarius group joins the list of kinematically-identified substructures that are not actually accreted material: they are simply part of the rich complexity of the Milky Way structure.

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.

Pulsars, if existing and detectable in the immediate vicinity of the massive black hole (MBH) in the Galactic center (GC), may be used as a superb tool to probe both the environment and the metric of the central MBH. The recent discovery of a magnetized pulsar in the GC suggests that many more pulsars should exist near the MBH. In this paper, we estimate the number and the orbital distribution of pulsars in the vicinity of the MBH in the GC by assuming that the pulsar progenitors, similar to the GC S-stars, were captured to orbits tightly bound to the MBH through the tidal breakup of stellar binaries. We use the current observations on both the GC S-stars and the hypervelocity stars to calibrate the injection rate(s) of and the dynamical model(s) for the stellar binaries. By including the relaxation processes, supernova kicks, and gravitational wave radiation in our simulations, we estimate that ~97-190 (9-14) pulsars may presently orbit the central MBH with semimajor axes <=4000AU (<=1000AU), which is compatible with the current observational constraints on the number of the GC pulsars. The semimajor axis and the pericenter distance of the pulsar closest to the central MBH are probably in the range of ~120-460AU and ~2-230AU, respectively. Future telescopes, such as the SKA, may be able to detect a significant number of pulsars with semimajor axis smaller than a few thousand AU in the GC. Long-term monitoring of these pulsars would be helpful in constraining both the environment and the metric of the central MBH. Our preferred model also results in about ten hyperfast pulsars with velocity >~1500km/s moving away from the Milky Way.

At the centre of a spherical and kinematically isotropic galaxy, the orientation of a massive binary black hole (BBH) orbit (i.e., the direction of the BBH orbital angular momentum) undergoes a random walk. If the stars in a spherical system have a non-zero total angular momentum, the BBH orbital orientation evolves towards aligning with the total stellar angular momentum direction. In this paper, we show that a triaxial galaxy has an alignment-erasing effect, that is, the alignment of the BBH orientations towards the galaxy rotation axis can be decreased significantly or erased. We also show that in a non-rotating axisymmetric galaxy, the BBH orbital orientation evolves towards the axisymmetric axis and precesses about it in a retrograde direction. Our results provide a step towards understanding the spin orientations of the final merged BH (and hence probable orientation of any jet produced) within its host galaxy, and may help to constrain the recoiling velocity of the merged BH arose from gravitational wave radiation as well.

In this paper, we investigate the link between the hypervelocity stars (HVSs) discovered in the Galactic halo and the S-stars moving in the Galactic center (GC), under the hypothesis that they are both the products of the tidal breakup of the same population of stellar binaries by the central massive black hole (MBH). By adopting several hypothetical models for binaries to be injected into the vicinity of the MBH and doing numerical simulations, we realize the tidal breakup processes of the binaries and their follow-up evolution. We find that many statistical properties of the detected HVSs and S-stars can be reproduced under some binary injecting models, and their number ratio can be reproduced if the stellar initial mass function is top-heavy (e.g., with slope ~-1.6). The total number of the captured companions is ~50 that have masses in the range ~3-7Msun and semimajor axes <~4000 AU and survive to the present within their main-sequence lifetime. The innermost one is expected to have a semimajor axis ~300-1500 AU and a pericenter distance ~10-200 AU, with a significant probability of being closer to the MBH than S2. Future detection of such a closer star would offer an important test to general relativity. The majority of the surviving ejected companions of the S-stars are expected to be located at Galactocentric distances <~20 kpc, and have heliocentric radial velocities ~-500-1500 km/s and proper motions up to ~5-20 mas/yr. Future detection of these HVSs may provide evidence for the tidal-breakup formation mechanism of the S-stars.

It is well established that stellar effective temperatures determined from photometry and spectroscopy yield systematically different results. We describe a new, simple method to correct spectroscopically derived temperatures ("excitation temperatures") of metal-poor stars based on a literature sample with -3.3<[Fe/H]<-2.5. Excitation temperatures were determined from FeI line abundances in high-resolution optical spectra in the wavelength range of ~3700 to ~7000A, although shorter wavelength ranges, up to 4750 to 6800A, can also be employed, and compared with photometric literature temperatures. Our adjustment scheme increases the temperatures up to several hundred degrees for cool red giants, while leaving the near-main-sequence stars mostly unchanged. Hence, it brings the excitation temperatures in good agreement with photometrically derived values. The modified temperature also influences other stellar parameters, as the FeI-FeII ionization balance is simultaneously used to determine the surface gravity, while also forcing no abundance trend on the absorption line strengths to obtain the microturbulent velocity. As a result of increasing the temperature, the often too low gravities and too high microturbulent velocities in red giants become higher and lower, respectively. Our adjustment scheme thus continues to build on the advantage of deriving temperatures from spectroscopy alone, independent of reddening, while at the same time producing stellar chemical abundances that are more straightforwardly comparable to studies based on photometrically derived temperatures. Hence, our method may prove beneficial for comparing different studies in the literature as well as the many high-resolution stellar spectroscopic surveys that are or will be carried out in the next few years.

We study the observational constraints on the cosmic evolution of the relationships between the massive black hole (MBH) mass (M_bh) and the stellar mass (M^*_sph; or velocity dispersion \sigma) of the host galaxy/spheroid. Assuming that the M_bh-M^*_sph (or M_bh-\sigma) relation evolves with redshift as \propto (1+z)^\Gamma, the MBH mass density can be obtained from either the observationally determined galaxy stellar mass functions or velocity dispersion distribution functions over redshift z~0-1.2 for any given \Gamma. The MBH mass density at different redshifts can also be inferred from the luminosity function of QSOs/AGNs provided known radiative efficiency \epsilon. By matching the MBH density inferred from galaxies to that obtained from QSOs/AGNs, we find that \Gamma=0.64^{+0.27}_{-0.29} for the M_bh-M^*_sph relation and \Gamma=-0.21^{+0.28}_{-0.33} for the M_bh-\sigma relation, and \epsilon=0.11^{+0.04}_{-0.03}. Our results suggest that the MBH mass growth precedes the bulge mass growth but the galaxy velocity dispersion does not increase with the mass growth of the bulge after the quench of nuclear activity, which is roughly consistent with the two-phase galaxy formation scenario proposed by Oser et al. (2012) in which a galaxy roughly double its masses after z=1 due to accretion and minor mergers while its velocity dispersion drops slightly.

In this paper we study possible signatures of binary planets or exomoons on the Rossiter-McLaughlin (R-M) effect. Our analyses show that the R-M effect for a binary planet or exomoon during its complete transit phase can be divided into two parts. The first is the conventional one similar to the R-M effect from the transit of a single planet, of which the mass and the projected area are the combinations of the binary components; and the second is caused by the orbital rotation of the binary components, which may add a sine- or linear-mode deviation to the stellar radial velocity curve. We find that the latter effect can be up to several or several ten m/s. By doing numerical simulations as well as analytical analyses, we illustrate that the distribution and dispersion of the latter effects obtained from multiple transit events can be used to constrain the dynamical configuration of the binary planet, such as, how the inner orbit of the binary planet is inclined to its orbit rotating around the central star. We find that the signatures caused by the orbital rotation of the binary components are more likely to be revealed if the two components of binary planet have different masses and mass densities, especially if the heavy one has a high mass density and the light one has a low density. Similar signature on the R-M effect may also be revealed in a hierarchical triple star system containing a dark compact binary and a tertiary star.

Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evolution according to the observationally-estimated major merger rates of galaxies and various scaling relations on the properties of galaxies and their central massive black holes. We show that our model reproduces the observed frequency and separation distribution of dual AGNs provided that significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with increasing redshift and only about 0.02%--0.06% of AGNs are dual AGNs with double-peaked narrow line features at redshifts of z 0.5-1.2. Future observations of high-redshift dual AGNs would provide a solid test for this prediction.

The proximity effect (PE) of QSOs is believed to be useful in constraining the QSO lifetime. Observations on the PE so far, however, give apparently contradictory results -- some are consistent with a long QSO lifetime (>~ a few 10^7 yr), but others appear to be only consistent with a short QSO lifetime <~ 10^6 yr. In this paper, we show that this apparent contradiction may be solved by simultaneously taking into account both the effect due to the density enhancement in the QSO near zones and that due to the obscuration of the tori associated with the QSOs, using a large number of Monte-Carlo generated synthetic Lyman alpha forest spectra. We demonstrate that the QSO properties and environment can be constrained simultaneously by the transverse PE and the line of sight PE of bright type 1 QSOs together. The current available measurements on the PEs of type 1 QSOs suggest that (1) the density is significantly enhanced in the vicinity of the QSOs; (2) the QSO lifetime is consistent with being as large as a few 10^7 yr and a substantially shorter lifetime (e.g., <~10^6 yr) is not required; and (3) the half opening angle of the tori associated with QSOs is ~60 deg, consistent with some other independent estimates. Our simulations also show that the TPE of type 2 QSOs can be significantly different from that of type 1 QSOs, which may be useful to put further constraints on the QSO properties and the QSO environment.

Hypervelocity stars (HVSs) found in the Galactic halo are probably the dynamical products of interactions between (binary) stars and the massive black hole(s) (MBH) in the Galactic center (GC). It has been shown that the detected HVSs are spatially consistent with being located on two thin disks (Lu et al.), one of which has the same orientation as the clockwise-rotating stellar disk in the GC. Here we perform a large number of three-body experiments of the interactions between the MBH and binary stars bound to it, and find that the probability of ejecting HVSs is substantially enhanced by multiple encounters between the MBH and binary stars at distances substantially larger than the initial tidal breakup radii. Assuming that the HVS progenitors are originated from the two disks, the inclination distribution of the HVSs relative to the disk planes can be reproduced by either the mechanism of tidal breakup of binary stars or the mechanism of ejecting HVSs by a hypothetical binary black hole (BBH) in the GC. However, an isotropical origination of HVS progenitors is inconsistent with the observed inclination distribution. Assuming that the HVSs were ejected out by the tidal breakup mechanism, its velocity distribution can be reproduced if their progenitors diffuse onto low angular momentum orbits slowly and most of the progenitors were broken up at relatively large distances due to multiple encounters. Assuming that the HVSs were ejected out by a BBH within the allowed parameter space in the GC, our simulations produce relatively flatter velocity spectra compared to the observed ones; however, the BBH mechanism cannot be statistically ruled out, yet. Future surveys of HVSs and better statistics of their spatial and velocity distributions should enable to distinguish the ejection mechanisms of HVSs and shed new light on the dynamical environment of the MBH.(abridged)

Hypervelocity stars (HVSs) escaping away from the Galactic halo are dynamical products of interactions of stars with the massive black hole(s) (MBH) in the Galactic Center (GC). They are mainly B-type stars with their progenitors unknown. OB stars are also populated in the GC, with many being hosted in a clockwise-rotating young stellar (CWS) disk within half a parsec from the MBH and their formation remaining puzzles. In this paper, we demonstrate that HVSs can well memorize the injecting directions of their progenitors using both analytical arguments and numerical simulations, i.e., the ejecting direction of an HVS is almost anti-parallel to the injecting direction of its progenitor. Therefore, the spatial distribution of HVSs maps the spatial distribution of the parent population of their progenitors directly. We also find that almost all the discovered HVSs are spatially consistent with being located on two thin disk planes. The orientation of one plane is consistent with that of the (inner) CWS disk, which suggests that most of the HVSs originate from the CWS disk or a previously existed disk-like stellar structure with an orientation similar to it. The rest of HVSs may be correlated with the plane of the northern arm of the mini-spiral in the GC or the plane defined by the outer warped part of the CWS disk. Our results not only support the GC origin of HVSs but also imply that the central disk (or the disk structure with a similar orientation) should persist or be frequently rejuvenated over the past 200 Myr, which adds a new challenge to the stellar disk formation and provides insights to the longstanding problem of gas fueling into massive black holes.

Growth of massive black holes (MBHs) in galactic centers comes mainly from gas accretion during their QSO/AGN phases. In this paper we apply an extended Soltan argument, connecting the local MBH mass function with the time-integral of the QSO luminosity function, to the demography of MBHs and QSOs from recent optical and X-ray surveys, and obtain robust constraints on the luminosity evolution (or mass growth history) of individual QSOs (or MBHs). We find that the luminosity evolution probably involves two phases: an initial exponentially increasing phase set by the Eddington limit and a following phase in which the luminosity declines with time as a power law (with a slope of -1.2--1.3) set by a self-similar long-term evolution of disk accretion. Neither an evolution involving only the increasing phase with a single Eddington ratio nor an exponentially declining pattern in the second phase is likely. The period of a QSO radiating at a luminosity higher than 10% of its peak value is about (2-3)x10^8 yr, during which the MBH obtains ~80% of its mass. The mass-to-energy conversion efficiency is $0.16\pm0.04 ^{+0.05}_{-0}$, with the latter error accounting for the maximum uncertainty due to Compton-thick AGNs. The expected Eddington ratios in QSOs from the constrained luminosity evolution cluster around a single value close to 0.5-1 for high-luminosity QSOs and extend to a wide range of lower values for low-luminosity ones. The Eddington ratios for high luminosity QSOs appear to conflict with those estimated from observations (~0.25) by using some virial mass estimators for MBHs in QSOs unless the estimators systematically over-estimate MBH masses by a factor of 2-4. We also infer the fraction of optically obscured QSOs ~60-80%. Further applications of the luminosity evolution of individual QSOs are also discussed.

The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center through interactions between binary stars and the central massive black hole or between single stars and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be ejected out of the Galactic center with velocities up to 10^3 km/s, while preserving their integrity, through interactions with a massive binary black hole. Binary stars are unlikely to attain such high velocities via scattering by a single massive black hole or through any other mechanisms. Based on the above theoretical prediction, we propose a search for binary systems among the hypervelocity stars. Discovery of hypervelocity binary stars, even one, is a definitive evidence of the existence of a massive binary black hole in the Galactic center.

Within a half-parsec from the Galactic center (GC), there is a population of coeval young stars which appear to reside in a coherent disk. Surrounding this dynamically-cool stellar system, there is a population of stars with a similar age and much larger eccentricities and inclinations relative to the disk. We propose a hypothesis for the origin of this dynamical dichotomy. Without specifying any specific mechanism, we consider the possibility that both stellar populations were formed within a disk some 6 Myr ago. But this orderly structure was dynamically perturbed outside-in by an intruding object with a mass ~10^4 Msun, which may be an intermediate-mass black hole (IMBH) or a dark stellar cluster hosting an IMBH. We suggest that the perturber migrated inward to ~0.15-0.3pc from the GC under the action of dynamical friction. Along the way, it captured many stars in the outer disk region into its mean-motion resonance, forced them to migrate with it, closely encountered with them, and induced the growth of their eccentricity and inclination. But stars in the inner regions of the disk retain their initial coplanar structure. We predict that some of the inclined and eccentric stars surrounding the disk may have similar Galactocentric semimajor axis. Future precision determination of their kinematic distribution of these stars will not only provide a test for this hypothesis but also evidences for the presence of an IMBH or a dark cluster at the immediate proximity of the massive black hole at the GC. (abridged)

Hypervelocity stars (HVSs) ejected by the massive black hole at the Galactic center have unique kinematic properties compared to other halo stars. Their trajectories will deviate from being exactly radial because of the asymmetry of the Milky Way potential produced by the flattened disk and the triaxial dark matter halo, causing a change of angular momentum that can be much larger than the initial small value at injection. We study the kinematics of HVSs and propose an estimator of dark halo triaxiality that is determined only by instantaneous position and velocity vectors of HVSs at large Galactocentric distances (r>~50kpc). We show that, in the case of a substantially triaxial halo, the distribution of deflection angles (the angle between the stellar position and velocity vector) for HVSs on bound orbits is spread uniformly over the range 10--180deg. Future astrometric and deep wide-field surveys should measure the positions and velocities of a significant number of HVSs, and provide useful constraints on the shape of the Galactic dark matter halo.

Observations have shown that the Eddington ratios (the ratio of the bolometric luminosity to the Eddington luminosity) in QSOs/active galactic nuclei (AGNs) cover a wide range. In this paper we connect the demography of AGNs obtained by the Sloan Digital Sky Survey with the accretion physics around massive black holes and propose that the diversity in the Eddington ratios is a natural result of the long-term evolution of accretion disks in AGNs. The observed accretion rate distribution of AGNs (with host galaxy velocity dispersion sigma~70-200 km/s) in the nearby universe (z<0.3) is consistent with the predictions of simple theoretical models in which the accretion rates evolve in a self-similar way. We also discuss the implications of the results for the issues related to self-gravitating disks, coevolution of galaxies and QSOs/AGNs, and the unification picture of AGNs.

In this paper we investigate the environment and reionization process around the highest redshift QSOs having Gunn-Peterson troughs (z>6.1). Starting with the cosmic density perturbation and structure formation theory and the fact that the highest redshift QSOs are located in rare overdense regions, we show that the halo formation, gas distribution, and star formation around QSOs are biased from those of the cosmic average. We argue that a significant fraction of hydrogen in the Str\"omgren sphere around QSOs is ionized by photons from stars and that only about several percent to at most 10%--20% of the total hydrogen is left (e.g., in minihalos, halos, or high-density subregions) to be ionized by QSO photons. The cosmic average neutral hydrogen fraction at z~6.2-6.4 should also be smaller than the upper limit of 10%--20% and may be only a few percent. We analyze the clumping property of the hydrogen ionized by QSOs and study the evolution of the Str\"omgren sphere. We find that the expected Str\"omgren radii from our models are consistent with observations if the lifetime of the highest redshift QSOs is about or longer than a few times 10^7 yr (as is the lifetime of the main population of QSOs; with comoving number density peaked at z~2-3). With such a QSO lifetime, the ages of most of the observed QSOs are long enough that the QSO photon emission is balanced by the recombination of the hydrogen ionized by QSO photons in their Str\"omgren spheres, and the expected Str\"omgren radii from the balance are independent of the detailed values of the QSO ages. (abridged)

Although the highest redshift QSOs (z>6.1) are embedded in a significantly neutral background universe (mass-averaged neutral hydrogen fraction >1%) as suggested by the Gunn-Peterson absorption troughs in their spectra, the intergalactic medium in their vicinity is highly ionized. The highly ionized region is generally idealized as spherical and called the Str\"omgren sphere. In this paper, by combining the expected evolution of the Str\"omgren sphere with the rule that the speed of light is finite, we illustrate the apparent shape of the ionization fronts around the highest redshift QSOs and its evolution, which depends on the age, luminosity evolution, and environment of the QSO (e.g., the hydrogen reionization history). The apparent shape may systematically deviate from a spherical shape, unless the QSO age is significantly long compared to the hydrogen recombination process within the ionization front and the QSO luminosity evolution is significantly slow. Effects of anisotropy of QSO emission are also discussed. The apparent shape of the "Str\"omgren sphere'' may be directly mapped by transmitted spectra of background sources behind or inside the ionized regions or by surveys of the hyperfine transition (21cm) line emission of neutral hydrogen.

Studies of massive black holes (BHs) in nearby galactic centers have revealed a tight correlation between BH mass and galactic velocity dispersion. In this paper we investigate how the BH mass versus velocity dispersion relation and the nuclear luminosity versus velocity dispersion relation in QSOs/active galactic nuclei (AGNs) are connected with the BH mass versus velocity dispersion relation in local galaxies, through the nuclear luminosity evolution of individual QSOs/AGNs and the mass growth of individual BHs. In the study we ignore the effects of BH mergers and assume that the velocity dispersion does not change significantly during and after the nuclear activity phase. Using the observed correlation in local galaxies and an assumed form of the QSO/AGN luminosity evolution and BH growth, we obtain the simulated observational appearance of the BH mass versus velocity dispersion relation in QSOs/AGNs. The simulation results illustrate how the BH accretion history (e.g., the lifetime of nuclear activity and the possibility that QSOs/AGNs accrete at a super-Eddington accretion rate at the early evolutionary stage) can be inferred from the difference between the relation in QSOs/AGNs and that in local galaxies. We also show how the difference may be weakened by the flux limit of telescopes. We expect that a large complete sample of QSOs/AGNs with accurate BH mass and velocity dispersion measurements will help to quantitatively constrain QSO/AGN luminosity evolution and BH growth models.

QSOs are believed to be powered by accretion onto massive black holes (BHs). In this paper, assuming that each central BH in nearby galaxies has experienced the QSO phase and ignoring BH mergers, we establish a relation between the QSO luminosity function (LF) and the local BH mass function (MF). The QSOLF is jointly controlled by the luminosity evolution of individual QSOs and the triggering history of the accretion onto seed BHs. By comparing the time integral of the QSOLF with that inferred from local BHs, we separate the effect of the luminosity evolution of individual QSOs from the effect of the triggering history. Assuming that the nuclear luminosity evolution includes two phases (first increasing at the Eddington luminosity with growth of BHs and then declining), we find that observations are generally consistent with the expected relation between the QSOLF and the local BHMF and obtain the constraints on QSO models and BH growth. We point out that it is hard to accurately estimate the value of the QSO lifetime from the QSOLF and/or the local BHMF, if it is longer than a certain value (e.g., four times of the Salpeter timescale in this study). We discuss the importance of accurate measurements of the intrinsic scatter in the BH mass and velocity dispersion relation of local galaxies and the scatter in the bolometric correction of QSOs. We also discuss some possible applications of the work in this paper, such as to the study of the demography of QSOs and the demography of normal galaxies at intermediate redshift (abridged).

We study three processes that eject hypervelocity (>10^3 km/s) stars from the Galactic center: (i) close encounters of two single stars; (ii) tidal breakup of binary stars by the central black hole, as originally proposed by Hills; and (iii) three-body interactions between a star and a binary black hole (BBH). Mechanism (i) expels hypervelocity stars to the solar radius at a negligible rate, ~10^{-11}/yr. Mechanism (ii) expels hypervelocity stars at a rate ~ 10^{-5}(\eta/0.1)/yr, where \eta is the fraction of stars in binaries with semimajor axis a_b<~0.3 AU. For solar-mass stars, the corresponding number of hypervelocity stars within the solar radius R_0=8 kpc is ~60(\eta/0.1)(a_b/0.1 AU)^{1/2}. For mechanism (iii), Sgr A^* is assumed to be one component of a BBH. We constrain the allowed parameter space (semimajor axis, mass ratio) of the BBH. In the allowed region (for example, semimajor axis of 0.5x10^{-3} pc and mass ratio of 0.01), the rate of ejecting hypervelocity stars can be as large as ~10^{-4}/yr and the expected number of hypervelocity stars within the solar radius can be as large as ~10^3. Hypervelocity stars may be detectable by the next generation of large-scale optical surveys.

We study the effects of stellar collisions, particularly on feeding massive black holes (BHs) and color gradients, in realistic galactic centers. We find that the mass released by stellar collisions is not sufficient to account for the present BH mass in galactic centers, especially in bright galaxies. This study, together with the study by Magorrian & Tremaine (1999) on tidal disruption of stars by massive BHs, implies that the material for BH growth (especially in galaxies brighter than ~10^9 Lsun) can only come from other sources, for example, the mass released by stellar evolution in the initial ~1 Gyr of the galaxy's lifetime, or the gas that sinks to the galactic center in a galaxy merger. We also analyze how the color of a stellar system is affected by collisions of stars. We find that collisions between main-sequence stars cannot cause observable color gradients in the visible bands at projected radius R>0.1" in M31, M32 and other nearby galactic centers. This result is consistent with the lack of an observable color gradient in M32 at R>0.1". At even smaller radii, the color differences caused by collisions between main-sequence stars are at most 0.08 mag at R=0.02". The averaged blueing due to stellar collisions in the region R<0.1" of M32 should not be larger than 0.06 mag in color index U-V and 0.02 mag in V-I. The observed blueing in the center of the galaxy M31 (in a 0.14"x0.14" box) must be caused by some mechanism other than collisions between main-sequence stars.

We study the observational constraints on the growth of massive black holes (BHs) in galactic nuclei. We use the velocity dispersions of early-type galaxies obtained by the SDSS and the relation between BH mass and velocity dispersion to estimate the local BH mass density to be 2.5x10^5 Msun/Mpc^3. We also use the QSO luminosity function from the 2dF Redshift Survey to estimate the BH mass density accreted during optically bright QSO phases. The local BH mass density is consistent with the density accreted during optically bright QSO phases if QSOs have an efficiency 0.1. By studying the continuity equation for the BH mass distribution, including the effect of BH mergers, we find relations between the local BH mass function and the QSO luminosity function. If the BH mass is assumed to be conserved during BH mergers, comparison of the predicted relations with the observations suggests that luminous QSOs (L_{bol}>10^{46} erg/s) have a high efficiency (e.g. 0.2), and the growth of high-mass BHs (>10^8 Msun) comes mainly from accretion during optically bright QSO phases, or that luminous QSOs have a super-Eddington luminosity. If luminous QSOs are not accreting with super-Eddington luminosities and the growth of low-mass BHs also occurs mainly during optically bright QSO phases, less luminous QSOs must accrete with a low efficiency <0.1; alternatively, they may accrete with high efficiency, but a significant fraction should be obscured. We estimate that the mean lifetime of luminous QSOs is (3-13)x10^7 yr, which is comparable to the Salpeter time. We also investigate the case in which total BH mass decreases during BH mergers due to gravitational radiation, and the observations again suggest that BHs in most luminous QSOs are Kerr BHs accreting with an efficiency >0.1.

Since many or most galaxies have central massive black holes (BHs), mergers of galaxies can form massive binary black holes (BBHs). In this paper, we study the evolution of massive BBHs in realistic galaxy models, using a generalization of techniques used to study tidal disruption rates around massive BHs. The evolution of BBHs depends on BH mass ratio and host galaxy type. BBHs with very low mass ratios (say, $\la$ 0.001) are hardly ever formed by mergers of galaxies because the dynamical friction timescale is too long for the smaller BH to sink into the galactic center within a Hubble time. BBHs with moderate mass ratios are most likely to form and survive in spherical or nearly spherical galaxies and in high-luminosity or high-dispersion galaxies; they are most likely to have merged in low-dispersion galaxies (line-of-sight velocity dispersion $\la$ 90 km/s) or in highly flattened or triaxial galaxies. The semimajor axes and orbital periods of surviving BBHs are generally in the range 10^{-3}-10 pc and 10-10^5 yr; and they are larger in high-dispersion galaxies than in low-dispersion galaxies, larger in nearly spherical galaxies than in highly flattened or triaxial galaxies, and larger for BBHs with equal masses than for BBHs with unequal masses. The orbital velocities of surviving BBHs are generally in the range 10^2-10^4 km/s. The methods of detecting surviving BBHs are also discussed. If no evidence of BBHs is found in AGNs, this may be either because gas plays a major role in BBH orbital decay or because nuclear activity switches on soon after a galaxy merger, and ends before the smaller BH has had time to spiral to the center of the galaxy.

We study the effects of the bulk motion of X-ray flares on the shape and equivalent width of the iron Kalpha line from an untruncated cold disk around a Kerr black hole using fully relativistic calculations. The flares are located above a cold accretion disk -- either on or off the rotation axis. For on- or off-axis flares, the upward/outward bulk motion causes a reduction of the iron Kalpha line width. To a distant observer with a low inclination angle (\theta_o \simlt 30deg.), larger upward/outward bulk velocities decrease the extension of the red wing, with little change in the location of the blue `edge'. In contrast, an observer at a large inclination angle (e.g. \theta_o=60deg.) sees both the red wing and the blue `edge' change with the bulk velocity. The equivalent width of the iron Kalpha line decreases rapidly with increasing bulk velocity of flares. However, the `narrower' line profiles observed in some objects (e.g. IC4329A and NGC4593) are difficult to produce using the out-flowing magnetic flare model with an appropriate equivalent width unless the X-ray emission is concentrated in an outer region with a radius of several tens of r_g=GM/c^2 or more. An important result is that the iron Kalpha line intensity is found to be constant even though the continuum flux varies significantly, which is true for out-flowing magnetic flares with different bulk velocities but similar intrinsic luminosities when located close to the central black hole. We find that fluctuations in the bulk velocities of out-flowing low-height flares located at the inner region (r\simlt 15r_g) can account for a constant iron Kalpha line and significant continuum variation as observered in MCG-6-30-15 and NGC5548. (Abridged)

During the recombination epoch, cosmic background photons couple not only to free electrons through Thompson scattering, but also to the neutral hydrogen through Rayleigh scattering. This latter is ~2% effect for photons near the peak of the photon energy distribution at z=800 and a ~0.2% effect at z=1100. Including Rayleigh scattering in the calculation reduces Silk damping at fixed redshift, alters the position of the surface of last scattering and alters the propagation of acoustic waves. We estimate the amplitude of these effects. For the Microwave Anisotropy Probe (MAP), Rayleigh scattering increases the anisotropy spectrum by 0.1% at the most. For the highest frequencies of the Planck Surveyor, the effects of Rayleigh scattering are much more dramatic (decreasing the anisotropy spectrum by 3% at \nu~550GHz and l~1000). The relative difference between the spectra of low and high frequencies is imposed by an oscillation with a function of multipole l and the oscillation amplitude is up to 0.5% between 100 and 550GHz. Rayleigh scattering also slows the decoupling between radiation and matter, but the effect is undetectably small.

We assembled a sample of Seyfert 1 galaxies, QSOs and Low-Luminosity Active Galactic Nuclei (LLAGNs) observed by ASCA, whose central black hole masses have been measured. We found that the X-ray variability (which is quantified by the ``excess variance'' sigma^2_{rms}) is significantly anti-correlated with the central black hole mass, and there likely exists a linear relationship sigma^2_{rms}\propto M^{-1}_{bh}. This can be interpreted that the short time-scale X-ray variability is caused by some global coherent variations in the X-ray emission region which is scaled by the size of the central black hole. Hence, the central black hole mass is the driving parameter of the previously established relation between X-ray variability and luminosity. This findings favor the hypothesis that the Narrow Line Seyfert 1 galaxies and QSOs harbor smaller black holes than the broad line objects, and can also easily explain the observational fact that high redshift QSOs have greater variability than local AGNs at a given luminosity. Further investigations are needed to confirm our findings, and a large sample X-ray variability investigation can give constraints on the physical mechanism and evolution of AGNs.

Hierarchical mergers of galaxies can form binary black holes (BBHs) since many or most galaxies have central massive black holes (BHs). It is possible that some BBHs exist in active galactic nuclei (AGNs). We argue that each BH may be surrounded by an accretion disc with a different inclination angle to the line of sight (due to different BH spin directions and the Bardeen-Petterson effect). The observed Fe K_alpha line profile from a BBH system is a combination of the lines from the inner regions of the two discs, which is significantly affected by the inclination angles of the two discs. The Fe K_alpha line profile associated with BBHs may have an unusual shape with double or more peaks as well as short-term variability, which can be distinguished from the Fe K_alpha line properties of some other possible models. We suggest that with the improvement of resolution in X-ray astronomy, Fe K_alpha line profiles be a potential tool to probe the existence of massive BBHs in AGNs. The Fe K_alpha line profile associated with BBHs may also provide a tool to investigate the dynamics in strong gravitation field (e.g. providing evidence of the Bardeen-Petterson effect).

We investigate resonant capture of small bodies by planets that migrate inwards, using analytic arguments and three-body integrations. If the orbits of the planet and the small body are initially circular and coplanar, the small body is captured when it crosses the 2:1 resonance with the planet. As the orbit shrinks it becomes more eccentric, until by the time its semimajor axis has shrunk by a factor of four, its eccentricity reaches nearly unity (1-e<<10^{-4}). In typical planetary systems, bodies in this high-eccentricity phase are likely to be consumed by the central star. If they can avoid this fate, as migration continues the inclination flips from 0 to i=180 degrees; thereafter the eccentricity declines until the semimajor axis is a factor of nine smaller than at capture, at which point the small body is released from the 2:1 resonance on a nearly circular retrograde orbit. Small bodies captured into resonance from initially inclined or eccentric orbits can also be ejected from the system, or released from the resonance on highly eccentric polar orbits (i\simeq 90 degrees) that are stabilized by a secular resonance. We conclude that migration can drive much of the inner planetesimal disk into the star, and that post-migration multi-planet systems may not be coplanar.

We use large-scale cosmological simulations to estimate the mass-to-light ratio of galaxy systems as a function of scale, and compare the results with observations of galaxies, groups, clusters, and superclusters of galaxies. We find remarkably good agreement between observations and simulations. Specifically, we find that the simulated mass-to-light ratio increases with scale on small scales and flattens to a constant value on large scales, as suggested by observations. We find that while mass typically follows light on large scales, high overdensity regions --- such as rich clusters and superclusters of galaxies --- exhibit higher M/L_B values than average, while low density regions exhibit lower M/L_B values; high density regions are thus antibiased in M/L_B, with mass more strongly concentrated than blue light. This is true despite the fact that the galaxy mass density is unbiased or positively biased relative to the total mass density in these regions. The M/L_B antibias is likely due to the relatively old age of the high density regions, where light has declined significantly since their early formation time, especially in the blue band which traces recent star formation. Comparing the simulated results with observations, we place a powerful constraint on the mass density of the universe; using, for the first time, the entire observed mass-to-light function, from galaxies to superclusters, we find Omega =0.16+/-0.05.

We suggest that polar rings and/or counter-rotating disks in flattened galaxies can be formed from stars captured at the Binney resonance, where the rate of precession of the angular momentum vector of a disk star equals the pattern speed of a triaxial halo. If the halo pattern speed is initially retrograde and slowly decays to zero, stars can be trapped as the Binney resonance sweeps past them, and levitated into polar orbits. If the halo pattern speed is initially retrograde and slowly changes to prograde, trapped stars can evolve from prograde to retrograde disk orbits. The stellar components of polar rings formed by this process should consist of two equal, counter-rotating star streams.

The mass of the central black hole in Seyfert galaxies and QSOs can be determined from the broad emission lines and the reverberation method. Using the measured black hole mass and the bolometric or ionizing luminosity, the accretion rate can be estimated. Compiling a sample of Seyfert 1 galaxies and QSOs with reliable central masses, estimated ionizing luminosities and X-ray spectral slopes, we find that the X-ray spectral slope strongly correlates with the accretion rate. The objects in the sample are found to be distributed in two distinct classes in the spectral index versus Log(L_ion/L_Edd) plane. We argue that these two classes may correspond to ADAF and thin disk accretion. The observations of a ``two-state'' Seyfert 1 galaxy, 1H0419-577, confirm our results. Detailed fitting of the spectra of individual Seyfert 1 galaxies and QSOs using ADAF and/or thin disk models should further clarify the two-class classification.

Previous calculations of Fe K-alpha line profiles are based on axisymmetric emissivity laws. In this paper, we show line profiles driven by non-axial symmetric illumination which results from an off-axis X-ray point source. We find that source location and motion have significant effects on the red wing and blue horn of the line profiles. The disk region under the source will receive more flux, which is the most important factor to affect the line profiles. We suggest that at least part of the variation in Fe K-alpha line profiles is caused by the motion of X-ray sources. Future observations of Fe K-alpha line profiles will provide more information about the distribution and motion of the X-ray sources around black holes, and hence the underlying physics.

Plutinos are Kuiper-belt objects that share the 3:2 Neptune resonance with Pluto. The long-term stability of Plutino orbits depends on their eccentricity. Plutinos with eccentricities close to Pluto (fractional eccentricity difference |e-e_p|/e_p<=0.1) can be stable because the longitude difference librates, in a manner similar to the tadpole and horseshoe libration in coorbital satellites. Plutinos with |e-e_p|/e_p>=0.3 can also be stable; the longitude difference circulates and close encounters are possible, but the effects of Pluto are weak because the encounter velocity is high. Orbits with intermediate eccentricity differences are likely to be unstable over the age of the solar system, in the sense that encounters with Pluto drive them out of the 3:2 Neptune resonance and thus into close encounters with Neptune. This mechanism may be a source of Jupiter-family comets.