QPOs are seen as peak features in the X-ray power spectral density of stellar mass black holes and neutron stars, of which frequencies can be measured with high precision. These QPO frequencies are believed to be related to fundamental frequencies of test particles, which are mainly determined by the background metric. We consider the metric introduced in Ghasemi-Nodehi and Bambi (Eur Phys J C 76:290, 2016). The fundamental frequencies in this metric are sensitive to some of the metric parameters but insensitive to other parameters, which means, the differences in fundamental frequencies in this metric and those in the Kerr ones can be significant for small changes of some but not all parameters around the Kerr value. By comparing with the QPO observations of GRO J1655-40, we find that only one parameter of the Ghasemi-Nodehi-Bambi metric can be strongly constrained, but other parameters cannot. We also use nested algorithm to investigate whether better constraints on the Ghasemi-Nodehi-Bambi metric parameters can be obtained from QPO observations of multiple objects by simulations. We find that four parameters can be strongly constrained while other parameters cannot . Our results suggest that QPOs may be important tools for testing the Kerr metric.

The spin distribution of massive black holes (MBHs) contains rich information on their assembly history. However, only limited information can be extracted from currently available spin measurements of MBHs owing to the small sample size and large measurement uncertainties. Upcoming X-ray telescopes with improved spectral resolution and larger effective area are expected to provide new insights into the growth history of MBHs. Here we investigate, at a proof of concept level, how stringent constraints can be placed on the accretion history of MBHs by the spin measurements from future X-ray missions. We assume a toy model consisting of a two-phase accretion history composed of an initial coherent phase with a constant disk orientation, followed by a chaotic phase with random disk orientations in each accretion episode. By utilizing mock spin data generated from such models and performing Bayesian Markov Chain Monte Carlo simulations, we find that most accretion models of MBHs can be reconstructed provided that $\gtrsim100$ MBH spins are measured with an accuracy of $\lesssim0.1$. We also quantify the precision of the reconstructed parameters by adopting various combinations of sample sizes and spin accuracies, and find that the sample size is more crucial to model reconstruction once the spin accuracy reaches $\sim 0.1$. To some extent, a better spin accuracy will compensate for a small sample size and vice versa. Future X-ray missions such as the Advanced Telescope for High Energy Astrophysics and the enhanced X-ray Timing and Polarimetry mission, may provide spin measurements of $\gtrsim100$ MBHs with an uncertainty of $\sim0.04-0.1$ and will thus put strong constraints on the MBH growth history.

We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the \startrack\ binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for \gaia\ and \lamost\ surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$--$340$ nBHBs to be observed by \gaia, although the numbers may drop to $\sim10$--$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For \lamost\ we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

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.

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of \textit{Kepler} data, we recover stellar properties, including surface gravity for red giant stars (with an uncertainty of $\lesssim$ 0.06 dex), and rotation period for main sequence stars (with an uncertainty of $\lesssim$ 5.2 days, and unbiased from $\approx$5 to 40 days). Shortening the \textit{Kepler} data to a 27-day TESS-like baseline, we recover stellar properties with a small decrease in precision, $\sim$0.07 dex for log $g$ and $\sim$5.5 days for $P_{\rm rot}$, unbiased from $\approx$5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal feature engineering, and can be generalized to other surveys and datasets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys.

We explore the prospects of detecting of Galactic double white dwarf (DWD) binaries with the space-based gravitational wave (GW) observatory TianQin. In this work we analyse both a sample of currently known and a realistic synthetic population of DWDs to assess the number of guaranteed detections and the full capacity of the mission. We find that TianQin can detect 12 out of $\sim100$ known DWDs; GW signals of these binaries can be modelled in detail ahead of the mission launch, and therefore they can be used as verification sources. Besides we estimate that TianQin has potential to detect as many as $10^4$ DWDs in the Milky Way. TianQin is expected to measure their orbital periods and amplitudes with accuracy of $\sim10^{-7}$ and $\sim0.2$ respectively, and to localize on the sky a large fraction (39%) of the detected population to better than 1 deg$^2$. We conclude that TianQin has the potential to significantly advance our knowledge on Galactic DWDs by increasing the sample up to 2 orders of magnitude, and will allow their multi-messenger studies in combination with electromagnetic telescopes. We also test the possibilities of different configuration of TianQin: 1) the same mission with a different orientation, 2) two perpendicular constellations combined into a network, and 3) the combination of the network with the ESA-lead Laser Interferometer Space Antenna. We find that the network of detectors boosts the accuracy on the measurement of source parameters by $1-2$ orders of magnitude, with the improvement on sky localization being the most significant.

The exact relationship between the long gamma-ray burst (LGRB) rate and the cosmic star formation rate (CSFR) is essential for using LGRBs as cosmological probes. In this work, we collect a large sample composed of 371 Swift LGRBs with known redshifts and prompt emission properties. We first compare the rest-frame prompt properties of these bursts in different redshift bins, finding negligible redshift evolution of the luminosity of LGRBs with $L_{\mathrm{iso}}\gtrsim10^{51}\,\mathrm{erg\, s^{-1}}$ between $z\sim1$ and $z\sim4$. Then, by utilizing the CSFR obtained from the large-scale cosmological hydrodynamical simulation, the Illustris simulation, we calculate the cumulative redshift distribution of LGRBs under different metallicity thresholds. After comparing with our sample, we find that the predictions with a moderate threshold between $0.3\,Z_{\odot}\leqslant Z_{\mathrm{th}}\leqslant1.0\,Z_{\odot}$ are consistent with the sample between redshift $0<z<3$, while at higher redshifts between $3<z<5$, all metallicity thresholds fit the data well. When changing to an empirical model based on observations, the predictions show similar results as well. After comparing with the metallicity distribution of the observed LGRB host galaxies between $0<z<1$, we confirm that the production of LGRBs in galaxies with super-solar metallicity is suppressed. Nevertheless, considering that a significant fraction of stars are born in sub-solar metallicity environments at $z\gtrsim3$, we suggest that, as a first approximation, LGRBs can be used as direct tracers of the CSFR in this redshift range.

The possibility that in the mass range around $10^{-12}\ M_\odot$ most of dark matter constitutes of primordial black holes (PBHs) is a very interesting topic. To produce PBHs with this mass, the primordial scalar power spectrum needs to be enhanced to the order of 0.01 at the scale $k\sim 10^{12}\ \text{Mpc}^{-1}$. The enhanced power spectrum also produces large secondary gravitational waves at the mHz band. A phenomenological delta function power spectrum is usually used to discuss the production of PBHs and secondary gravitational waves. Based on G and k inflations, we propose a new mechanism to enhance the power spectrum at small scales by introducing a non-canonical kinetic term $[1-2G(\phi)]X$ with the function $G(\phi)$ having a peak. Away from the peak, $G(\phi)$ is negligible and we recover the usual slow-roll inflation which is constrained by the cosmic microwave background anisotrpy observations. Around the peak, the slow-roll inflation transiently turns to ultra slow-roll inflation. The enhancement of the power spectrum can be obtained with generic potentials, and there is no need to fine tune the parameters in $G(\phi)$. The energy spectrum $\Omega_{GW}(f)$ of secondary gravitational waves have the characteristic power law behaviour $\Omega_{GW}(f)\sim f^{n}$ and is testable by pulsar timing array and space based gravitational wave detectors.

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

We report the discovery of a Tidal Disruption Event (TDE) candidate occurring in the Active Galactic Nucleus SDSS J022700.77-042020.6. A sudden increase in flux of J0227-0420 during the second half of 2009 is clearly shown in the long-term optical, UV, and NIR light curves. A plateau phase, following an initial decline, is seen in the NUV and optical u, g, r, i light curves. Moreover, we find possible evidence that the plateau phase in the NUV band may lag behind the optical ones by approximately 70-80 days with also a much shorter duration, i.e. $\sim$7-15 days against $\sim$40-50 days. The long-term NUV/optical (after the plateau phase), NIR and MIR light curves can be well characterized with a form of $L(t)\propto t^{-\beta}$, consistent with the expectation of a TDE. The plateaus can be explained if the stellar streams collide with the pre-existing AGN disk at different radii. Though the overall fallback rate decreases, the material in the outer disk gradually drifts inward and increases the local accretion rate at the inner region, producing the optical and UV plateaus. The possible lag between the optical and NUV plateaus can then be attributed to viscosity delay. The index $\beta$ of the NIR $J, H, K_s$ bands ($\sim1.4-3.3$) is steeper than that of the UV/optical ($\sim0.7-1.3$) and MIR bands ($\sim0.9-1.8$), which may suggest that a certain fraction of the dust in the inner region of the dusty torus may be sublimated during the TDE phase. Our results indicate that, due to collisions between stellar debris and pre-existing disk, the light curves of TDEs occurring in AGN may show distinctive features, which may shed new light on the accretion process.

The Moon's optical characteristics in visible and long-wavelength infrared (LWIR) have long been observed with our eyes or with instruments. What the mid-infrared (MIR) Moon looks like is still a mystery. For the first time we present detailed appearance of the MIR Moon observed by a high-resolution geostationary satellite and reveal the essence behind its appearance. The appearance of the MIR Moon is opposite to its normal visible appearance. In addition the MIR Moon shows limb darkening. Both the absolute and the relative brightness distribution of the MIR lunar disk changes with the solar incidence angle. The signatures of the MIR Moon are controlled by both the reflection and emission of the lunar surface. We also show first-ever brightness temperature maps of the lunar disk without needing a mosaic, which better show the temperature variation across the lunar disk. They reveal that the relationship between brightness temperature and solar incidence angle i is cos1/bi, and the power parameter is smaller than the Lambertian temperature model of cos1/4i observed for lunar orbit-based measurements. The slower decrease of the brightness temperature when moving away from the sub-solar point than the Lambertian model is due to topographic effects. The brightness temperature is dominated by albedo and the solar incidence angle and influenced by the topography. Our results indicate that the Moon in the MIR exhibits many interesting phenomena which were previously unknown, and contains abundant information about lunar reflection and thermal emission for future study.

The discovery of gravitational-wave (GW) signals, produced by the coalescence of stellar-mass binary black holes (SBBHs), opens a new window to study the astrophysical origins and dynamical evolutions of compact binaries. In addition, these GW events can be treated as the standard sirens to constrain various cosmological parameters. Both issues require the host identification for these GW events, with help of the spatial resolution of GW detector networks. In this paper, we investigate the capabilities of various detector networks for identifying the SBBHs' host galaxy groups, rather than their host galaxies, which can overcome the influence of galaxies' proper motions in dark matter halos for measuring the cosmological parameters. In our analysis, the group catalog of SDSS DR7 with redshift $z\in(0.01,0.1)$ is considered as an example of the application. We find that for the second-generation (2G) detector network, the host galaxy groups of around $17\%$ $30-30\ M_{\odot}$ SBBHs can be identified. For the 3G detector network, this fraction becomes $\sim 99\%$. We also investigate the potential constraint on the Hubble constant $H_0$ by these GW events, if their redshift information is extracted from the candidates of host galaxy groups. We find that, by five-year observations, 2G detector network is expected to give a constraint of $\Delta H_0/H_0\sim (1\%,4\%)$, which can be more than two order smaller if considering the 3G detector network.

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

Primordial black holes and secondary gravitational waves can be used to probe the small scale physics at very early time. For secondary gravitational waves produced after the horizon reentry, we derive an analytical formula for the time integral of the source and analytical behavior of the time dependence of the energy density of induced gravitational waves is obtained. By proposing a piecewise power law parametrization for the power spectrum of primordial curvature perturbations, we use the observational constraints on primordial black hole dark matter to obtain an upper bound on the power spectrum, and discuss the test of the model with future space based gravitational wave antenna.

Supermassive binary black holes (SMBBHs) in galactic centers may radiate gravitational wave (GW) in the nano-Hertz frequency band, which are expected to be detected by pulsar timing arrays (PTAs) in the near future. GW signals from individual SMBBHs at cosmic distances, if detected by PTAs, are potentially powerful standard sirens that can be used to independently measure distances and thus put constraints on cosmological parameters. In this paper, we investigate the constraint that may be obtained on the equation of state ($w$) of dark energy by using those SMBBHs, expected to be detected by the PTAs in the Square Kilometre Array (SKA) era. By considering both the currently available SMBBH candidates and mock SMBBHs in the universe resulting from a simple galaxy major merger model, we find that $\sim 200$ to $3000$ SMBBHs with chirp mass $>10^9M_\odot$ are expected to be detected with signal-to-noise ratio $>10$ by SKA-PTA with conservative and optimistic settings and they can be used to put a constraint on $w$ to an uncertainty of $\Delta w\sim 0.02-0.1$. If further information on the mass and mass ratio of those SMBBHs can be provided by electromagnetic observations (e.g., chirp mass uncertainty $\lesssim 50\%$), the constraint may be further improved to $\lesssim 0.01$ level, as many more SMBBHs will be detected by SKA-PTA with relatively better distance measurements and can be used as the standard sirens.

Optical periodicity QSOs found by transient surveys are suggested to be sub-parsec supermassive binary black holes (BBHs). An intriguing interpretation for the periodicity of some of those QSOs is that the continuum is radiated from the accretion disk associated with the BBH secondary component and modulated by the periodical rotation of the secondary via Doppler-boost effect. Close to edge-on orbital orientation can lead to more significant Doppler-boost effect and thus are preferred for these systems, which is distinct from those normal type-1 QSOs with more or less face-on orientations. Therefore, the profiles of broad lines emitted from these Doppler-modulated systems may be significantly different from other systems that are not Doppler-modulated. We investigate the properties of the broad emission lines of optical-periodicity QSOs, including both a sample of QSOs that can be interpreted by the Doppler-modulated effects and a sample that cannot. We find that there is no obvious difference in the profiles and other properties of various (stacked) broad emission lines of these two samples, though a simple broad line region model would suggest significant differences. Our finding raises a challenge to the Doppler boost hypothesis for some of those BBHs candidates with optical periodicity.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

We present the observational results from a detailed timing analysis of the black hole candidate Swift J1658.2-4242 during its 2018 outburst with the observations of Hard X-ray Modulation Telescope (Insight-HXMT), Neutron Star Interior Composition Explorer (NICER) and AstroSat in 0.1-250keV. The evolution of intensity, hardness and integrated fractional root mean square (rms) observed by Insight-HXMT and NICER are presented in this paper. Type-C quasi-periodic oscillations (QPOs) observed by NICER (0.8-3.5Hz) and Insight-HXMT (1-1.6Hz) are also reported in this work. The features of the QPOs are analysed with an energy range of 0.5-50keV. The relations between QPO frequency and other characteristics such as intensity, hardness and QPO rms are carefully studied. The timing and spectral properties indicate that Swift J1658.2-4242 is a black hole binary system. Besides, the rms spectra of the source calculated from the simultaneous observation of Insight-HXMT, NICER and AstroSat support the Lense-Thirring origin of the QPOs. The relation between QPO phase lag and the centroid frequency of Swift J1658.2-4242 reveals a near zero constant when < 4Hz and a soft phase lag at 6.68 Hz. This independence follows the same trend as the high inclination galactic black hole binaries such as MAXI J1659-152.

Using the archive data from the Rossi X-ray Timing Explorer ({\sl RXTE}), we have studied the evolution of the X-ray profile of the Crab pulsar in a time span of 11 years. The X-ray profiles, as characterized by a few parameters, changed slightly but significantly in these years: the separation of the two peaks increased with a rate $0.88\pm0.20\,\textordmasculine$\,per century, the flux ratio of the second pulse to the first pulse decreased with $(3.64\pm0.86)\times10^{-2}$\,per century, and the pulse widths of the two pulses descended with $1.44\pm0.15\,\textordmasculine$\, and $1.09\pm0.73\,\textordmasculine$\,per century, respectively. The evolutionary trends of the above parameters are similar to the radio results, but the values are different. We briefly discussed the constraints of these X-ray properties on the geometry of the emission region of this pulsar.

With large-scale galaxy surveys, we can observe hundreds of thousands of galaxies or more, up to billions with upcoming experiments such as WFIRST, Euclid and LSST. While such surveys cannot obtain spectra for all observed galaxies, we have access to the galaxy magnitudes in color filters. This data set behaves like a high-dimensional nonlinear surface, making it an excellent target for machine learning methods. In this work, we use a lightcone of semianalytic galaxies tuned to match CANDELS observations from Lu et al. (2014) to train a set of neural networks on a set of galaxy properties (stellar mass, metallicity, and average star formation rate) using the truth values from the semianalytic catalogs. We also demonstrate the effect of adding simulated observational noise to the simulated magnitudes, and then use neural networks trained on the noisy data to predict stellar masses, metallicities, and average star formation rates on real CANDELS galaxies, comparing our results to the physical parameters obtained from SED fitting. On semianalytic galaxies alone, we are nearly competitive with template-fitting methods. For the observed CANDELS data, our results are not as accurate, with indications that this inaccuracy is due to a combination of different assumptions in template-fitting and differences between the semianalytic models and the observed galaxies, particularly in the noise properties. Our results show that stellar mass, metallicity, and star formation rate can in principle be measured with neural networks at a competitive degree of accuracy and precision relative to physically motivated template-fitting methods if an appropriate training set can be obtained.

We studied the long-term evolution of the spectral-temporal correlated properties of the black hole candidate Swift J1753.5-0127 from the onset of its outburst until 2011 with the Rossi X-ray Timing Explorer (RXTE). The source stayed most of its lifetime during hard state, with occasionally transitioned to the hard intermediate state. Similar to typical black hole transients, Swift J1753.5-0127 traces a clear hard line in absolute rms - intensity diagram during the low hard state, with expected highest absolute rms, while shows a clear turn during the hard intermediate state, accompanied by lower absolute rms. Different from Cyg X-1, we found that frequency-dependent time lag increased significantly in the 0.02 -- 3.2 Hz band during state transition in this source. The X-ray time lags in 0.02 -- 3.2 Hz can therefore be used as indicators of state transition in this source. Type-C QPO frequency is positively related with its fractional rms and X-ray photon index, suggesting a moving inwards disk/corona scenario. We discussed the physical interpretation of our results in our paper.

Accretion of tidally disrupted stars (TDSs) is expected to contribute significantly to the growth of massive black holes (MBHs) with mass $\sim 10^6 M_{\odot}$ in galactic centers. In this paper, we quantitatively investigate the effect of the TDS accretion on the spin evolution of these relatively small MBHs, by also considering the accretion of gas-clouds with (many) chaotic episodes. We find that the accretion of TDSs can play an important role or even a dominant role in shaping the spin distribution of $\sim 10^6 M_{\odot}$ MBHs, depending on the contribution fraction ($f_{\rm TDE}$) of the TDS accretion to the MBH growth. If $f_{\rm TDE}$ is as large as $\gtrsim 0.9$, most $\sim 10^6M_{\odot}$ MBHs have low spins ($|a| \lesssim 0.3$); if $f_{\rm TDE}$ is as small as $\lesssim 0.1$, most $\sim 10^6 M_{\odot}$ MBHs have high spins ($|a| \gtrsim 0.7$). We also find that (1) the fraction of highly spinning $\sim 10^6 M_{\odot}$ MBHs in the TDS accretion states is smaller than that in the gas-cloud (AGN) accretion states, which is a consequence of more rapid spin decrease during the period of consecutive TDS accretion than the spin increase during the AGN periods when the spin is large; (2) the fraction of retrograde spin accretion in the TDS accretion states is almost the same as that of prograde spin accretion, while it is negligible in the gas-cloud (AGN) accretion states. Current scarce sample of AGNs ($\sim 10^6 M_{\odot}$) with spin measurements hints an insignificant contribution from TDS accretion to MBH growth. Future measurements on spins of $\sim 10^6 M_{\odot}$ MBHs may provide stronger constraints on the importance of both AGN and TDS accretion states in their growth history.

We study the shape of the gas-phase mass-metallicity relation (MZR) of a combined sample of present-day dwarf and high-mass star-forming galaxies using IZI, a Bayesian formalism for measuring chemical abundances presented in Blanc et al. 2015. We observe a characteristic stellar mass scale at $M_* \simeq 10^{9.5}$M$_{\odot}$, above which the ISM undergoes a sharp increase in its level of chemical enrichment. In the $10^{6}-10^{9.5}$M$_{\odot}$ range the MZR follows a shallow power-law ($Z\propto M^{\alpha}_*$) with slope $\alpha=0.14\pm0.08$. At approaching $M_* \simeq 10^{9.5}$M$_{\odot}$ the MZR steepens significantly, showing a slope of $\alpha=0.37\pm0.08$ in the $10^{9.5}-10^{10.5}$M$_{\odot}$ range, and a flattening towards a constant metallicity at higher stellar masses. This behavior is qualitatively different from results in the literature that show a single power-law MZR towards the low mass end. We thoroughly explore systematic uncertainties in our measurement, and show that the shape of the MZR is not induced by sample selection, aperture effects, a changing N/O abundance, the adopted methodology used to construct the MZR, secondary dependencies on star formation activity, nor diffuse ionized gas (DIG) contamination, but rather on differences in the method used to measure abundances. High resolution hydrodynamical simulations can qualitatively reproduce our result, and suggest a transition in the ability of galaxies to retain their metals for stellar masses above this threshold. The MZR characteristic mass scale also coincides with a transition in the scale height and clumpiness of cold gas disks, and a typical gas fraction below which the efficiency of star formation feedback for driving outflows is expected to decrease sharply.

Radiative efficiency is an important physical parameter that describes the fraction of accretion material converted to radiative energy for accretion onto massive black holes (MBHs). With the simplest So{\l}tan argument, the radiative efficiency of MBHs can be estimated by matching the mass density of MBHs in the local universe to the accreted mass density by MBHs during AGN/QSO phases. In this paper, we estimate the local MBH mass density through a combination of various determinations of the correlations between the masses of MBHs and the properties of MBH host galaxies, with the distribution functions of those galaxy properties. We also estimate the total energy density radiated by AGNs and QSOs by using various AGN/QSO X-ray luminosity functions in the literature. We then obtain several hundred estimates of the mean radiative efficiency of AGNs/QSOs. Under the assumption that those estimates are independent of each other and free of systematic effects, we apply the median statistics as described by Gott et al.\cite{got01} and find the mean radiative efficiency of AGNs/QSOs is $\epsilon=0.105^{+0.006}_{-0.008}$, which is consistent with the canonical value $\sim0.1$. Considering that about $20\%$ Compton-thick objects may be missed from current available X-ray surveys, the true mean radiative efficiency may be actually $\sim 0.12$.

The spin distribution of massive black holes (MBHs) contains rich information on the MBH growth history. In this paper, we investigate the spin evolution of MBHs by assuming that each MBH experiences two-phase accretion, with an initial phase of coherent-accretion via either the standard thin disc or super-Eddington disc, followed by a chaotic-accretion phase composed of many episodes with different disc orientations. If the chaotic-phase is significant to the growth of an MBH, the MBH spin quickly reaches the maximum value because of the initial coherent-accretion, then changes to a quasi-equilibrium state and fluctuates around a value mainly determined by the mean ratio of the disc to the MBH mass ($M_{\bullet}$) in the chaotic-accretion episodes, and further declines due to late chaotic-accretion if $M_\bullet \gtrsim (1-3) \times 10^8 M_\odot$. The turning point to this decline is determined by the equality of the disc warp radius and disc size. By matching the currently available spin measurements with mock samples generated from the two-phase model(s) on the spin-mass plane, we find that MBHs must experience significant chaotic-accretion phase with many episodes and the mass accreted in each episode is roughly 1-2 percent of M_bh or less. MBHs with $M_{\bullet}\gtrsim 10^8 M_{\odot}$ appear to have intermediate-to-high spins ($\sim 0.5-1$), while lighter MBHs have higher spins ($\gtrsim 0.8$). The best matches also infer that (1) the radiative efficiencies ($\eta$) of those active MBHs appear to slightly decrease with $M_{\bullet}$; however, the correlation between $\eta$ and $M_{\bullet}$, if any, is weak; (2) the mean radiative efficiency of active MBHs is $<\eta> \sim 0.09-0.15$, consistent with the global constraints.

Three basic ingredients are required to generate a simple stellar population (SSP) library, i.e., an initial mass function (IMF), a stellar evolution model/isochrones, and an empirical/theoretical stellar spectral library. However, there are still some uncertainties to the determination and understanding of these ingredients. We perform the spectral fitting to test the relative parameter offsets between these uncertainties using two different stellar population models, two different empirical stellar libraries, two different isochrones, and the Salpeter and Chabrier IMFs. Based on these setups, we select five SSP libraries generated with the Galaxev/STELIB and Vazdekis/MILES models, and apply them to the pPXF full-spectrum fitting of both MaNGA and mock spectra. We find that: 1) Compared to the Galaxev/STELIB model, spectral fitting qualities with the Vazdekis/MILES model have significant improvements for those metal-rich (especially over-solar) spectra, which cause better reduced $\chi^2$ distributions and more precisely fitted absorption lines. This might due to the lack of metal rich stars in the empirical STELIB library, or code improvement of the Vazdekis model. 2) When applying the Vazdekis/MILES model for spectral fitting, the IMF variation will lead to not only a systematic offset in $M_*/L_r$, but also offsets in age and metallicity, and these offsets increase with increasing stellar population ages. However, the IMF-variation caused metallicity offsets disappear in the case of Galaxev/STELIB based libraries. 3) The Padova2000 model provides a better match to the MaNGA galaxy spectra at [M/H]$_L<-1.0$, while the BaSTI model match the local galaxy spectra better at [M/H]$_L>-1.0$. Current tests suggest that spectral fitting with the Vazdekis/MILES+BaSTI combination would be a better choice for local galaxies.

We investigate the chemical evolution of a collapsing core that starts from a hydrostatic core and finally form a low-mass protostar. New multiphase gas-grain models that include bulk diffusion and photon penetration are simulated by the macroscopic Monte Carlo method in order to derive the chemical evolution. There are two types of species in the ice bulk in the new multiphase models. Interstitial species can diffuse and sublime at their own sublimation temperatures while normal species are locked in the ice bulk. Photodissociation rates of icy species are reduced by the exponential decay of UV flux within the ice mantle. Two-phase models and basic multiphase models without bulk diffusion and photon penetration are also simulated for comparison. Our physical model for the collapsing core is base on a one-dimensional radiation hydrodynamics model. Abundant icy radicals are produced at around 10 K in the new multiphase models. Interstitial radicals can diffuse inside ice mantles to form complex organic molecules (COMS) upon warming-up. Thus, COMs produced by radical recombination at higher temperatures in the new multiphase models are more than one order of magnitude higher than those in the two-phase and basic multiphase models. Moreover, COMs produced at around 10 K in the new multiphase models are about one order of magnitude higher than those in the two-phase model. Our model shows a reasonable agreement with observations toward low-mass protostars. Moreover, molecular oxygen abundances predicted by our new multiphase models agree reasonably well with that found in cometary materials.

Space weathering is an important surface process occurring on the Moon and other airless bodies, especially those that have no magnetic field. The optical effects of the Moon's space weathering have been largely investigated in the laboratory for lunar samples and lunar analogues. However, duplication of the pristine regolith on Earth is not possible. Here we report the space weathering from the unique perspective of the Chang'E-3's (CE-3) "Yutu" rover, building on our previous work (Wang et al. 2017; Wu and Hapke 2018). Measurement of the visually undisturbed uppermost regolith as well as locations that have been affected by rocket exhaust from the spacecraft by the Visible-Near Infrared Spectrometer (VNIS) revealed that the returned samples bring a biased information about the pristine lunar regolith. The uppermost surficial regolith is much more weathered than the regolith immediately below, and the finest fraction is rich in space weathered products. These materials are very dark and attenuated throughout the visible and near-infrared (VNIR) wavelengths, hence reduce the reflectance and mask the absorption features. The effects on the spectral slope caused by space weathering are wavelength-dependent: the visible and near-infrared continuum slope (VNCS) increases while the visible slope (VS) decreases. In the visible wavelength, the optical effects of space weathering and TiO2 are identical: both reduce albedo and blue the spectra. This suggests that developing new TiO2 abundance algorithm is needed. Optical maturity indices are composition related and hence only locally meaningful. Since optical remote sensing can only sense the uppermost few microns of regolith and since this surface tends to be very weathered, the interpretation of surface composition using optical remote sensing data needs to be carefully evaluated. Sampling the uppermost surface is suggested.

The mineralogy of mare basalts reflects the chemical composition of the magma source, as well as the physical and chemical environment of the rock's formation. It is significant for understanding the thermal evolution of the Moon. In this study, the spatial distribution of the mineralogy of the lunar northern nearside basalts was mapped using the Moon Mineralogy Mapper (M3) data. The study area, which is an elongated mare, Mare Frigoris and northern Mare Imbrium, was mapped and characterized into 27 units based on multi-source data, including spectrum, terrain, and element abundance. We extracted 177 M3 spectra from fresh craters. Spectral parameters such as absorption center and band area ratio (BAR) were obtained through data processing. The variation of mafic mineral of this region was acquired by analyzing these parameters. The basaltic units in eastern Mare Frigoris, which are older, have been found to be dominated by clinopyroxene with lower CaO compared to the returned lunar samples; this is similar to older basaltic units in Mare Imbrium. The basaltic units of western Mare Frigoris and Sinus Roris which are younger have been found to be rich in olivine. The late-stage basalts in Oceanus Procellarum and Mare Imbrium show the same feature. These wide-spread olivine-rich basalts suggest the uniqueness in the evolution of the Moon. Geographically speaking, Mare Frigoris is an individual mare, but the magma source region have connections with surrounding maria in consideration of mineral differences between western and eastern Frigoris, as well as mineral similarities with maria at the same location.

Gravitational waves emitted from compact binary coalescence can be subject to wave diffraction if they are gravitationally lensed by an intervening mass clump whose Schwarzschild timescale matches the wave period. Waves in the ground-based frequency band $f\sim 10$--$10^3\,$Hz are sensitive to clumps with masses $M_E \sim 10^2$--$10^3\,M_\odot$ enclosed within the impact parameter. These can be the central parts of low mass $M_L \sim 10^3$--$10^6\,M_\odot$ dark matter halos, which are predicted in Cold Dark Matter scenarios but are challenging to observe. Neglecting finely-tuned impact parameters, we focus on lenses aligned generally on the Einstein scale for which multiple lensed images may not form in the case of an extended lens. In this case, diffraction induces amplitude and phase modulations whose sizes $\sim 10\%$--$20\%$ are small enough so that standard matched filtering with unlensed waveforms do not degrade, but are still detectable for events with high signal-to-noise ratio. We develop and test an agnostic detection method based on dynamic programming, which does not require a detailed model of the lensed waveforms. For pseudo-Jaffe lenses aligned up to the Einstein radius, we demonstrate that a pair of fully upgraded aLIGO/Virgo detectors can extract diffraction imprints from binary black hole mergers out to $z_s \sim 0.2$--$0.3$. The prospect will improve dramatically for a third-generation detector for which binary black hole mergers out to $z_s \sim 2$--$4$ will all become valuable sources.

We analyze 168 Swift monitoring observations of the nearest broad absorption line quasar Mrk 231 in the UV and X-ray bands, where we detect significant variability in the UV ($\sim$2246\AA) light curve with a null probability of $4.3\times10^{-10}$ for a constant model. Separately, from an archival sample of Swift observed active galactic nuclei (AGN), we measure the relation between UV excess variance and luminosity, finding that the normalized UV excess variance decreases with luminosity. Comparing to this mean relation, the normalized UV excess variance of Mrk 231 is smaller, however within the scatter characterising the full population. The upper limit of the X-ray excess variance is consistent with other AGN. The power spectrum density of the UV light curve can be well fit by a power law model with a slope of $1.82\pm0.14$ between $10^{-7.5}$ and $10^{-6}$ Hz, consistent with those for typical AGN, with no obvious quasi-periodical oscillation peaks. The UV variability and its power spectrum suggest that a significant amount of the UV emission of Mrk 231 is from the accretion disk. The consistencies in the normalized UV variability and the shape of the power spectrum density between Mrk 231 and other normal AGN suggest that the origin of UV variability of broad absorption line quasars is similar to other AGN, and dust scattering at large scales such as the torus is not a dominating process for the UV emission of Mrk 231. Significant scattering, if present, is constrained to smaller than $\sim$10 light days. We perform lagged correlation analysis between the UV and X-ray light curves and find the correlation insignificant within the present data.

In this paper, we report the discovery of a new late-B type unbound hyper-runaway star (LAMOST-HVS4) from the LAMOST spectroscopic surveys. According to its atmospheric parameters, it is either a B-type main sequence (MS) star or a blue horizontal branch (BHB) star. Its Galactocentric distance and velocity are 30.3 +/- 1.6 kpc and 586 +/- 7 km/s if it is an MS star, and they are 13.2 +/- 3.7 kpc and 590 +/- 7 km/s if a BHB star. We track its trajectories back, and find that the trajectories intersect with the Galactic disk and the Galactic center lies outside of the intersection region at the 3 sigma confidence level. We investigate a number of mechanisms that could be responsible for the ejection of the star, and find that it is probably ejected from the Galactic disk by supernova explosion or multiple-body interactions in dense young stellar clusters.

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

Using mock spectra based on Vazdekis/MILES library fitted within the wavelength region 3600-7350\AA, we analyze the bias and scatter on the resulting physical parameters induced by the choice of fitting algorithms and observational uncertainties, but avoid effects of those model uncertainties. We consider two full-spectrum fitting codes: pPXF and STARLIGHT, in fitting for stellar population age, metallicity, mass-to-light ratio, and dust extinction. With pPXF we find that both the bias in the population parameters and the scatter in the recovered logarithmic values follows the expected trend. The bias increases for younger ages and systematically makes recovered ages older, $M_*/L_r$ larger and metallicities lower than the true values. For reference, at S/N=30, and for the worst case ($t=10^8$yr), the bias is 0.06 dex in $M_*/L_r$, 0.03 dex in both age and [M/H]. There is no significant dependence on either E(B-V) or the shape of the error spectrum. Moreover, the results are consistent for both our 1-SSP and 2-SSP tests. With the STARLIGHT algorithm, we find trends similar to pPXF, when the input E(B-V)<0.2 mag. However, with larger input E(B-V), the biases of the output parameter do not converge to zero even at the highest S/N and are strongly affected by the shape of the error spectra. This effect is particularly dramatic for youngest age, for which all population parameters can be strongly different from the input values, with significantly underestimated dust extinction and [M/H], and larger ages and $M_*/L_r$. Results degrade when moving from our 1-SSP to the 2-SSP tests. The STARLIGHT convergence to the true values can be improved by increasing Markov Chains and annealing loops to the "slow mode". For the same input spectrum, pPXF is about two order of magnitudes faster than STARLIGHT's "default mode" and about three order of magnitude faster than STARLIGHT's "slow mode".

We propose a novel method to constrain the missing fraction of galaxies using galaxy clustering measurements in the galaxy conditional stellar mass function (CSMF) framework, which is applicable to surveys that suffer significantly from sample selection effects. The clustering measurements, which are not sensitive to the random sampling (missing fraction) of galaxies, are widely used to constrain the stellar-halo mass relation (SHMR). By incorporating a missing fraction (incompleteness) component into the CSMF model (ICSMF), we use the incomplete stellar mass function and galaxy clustering to simultaneously constrain the missing fractions and the SHMRs. Tests based on mock galaxy catalogs with a few typical missing fraction models show that this method can accurately recover the missing fraction and the galaxy SHMR, and hence provides us reliable measurements of the galaxy stellar mass functions. We then apply it to the Baryon Oscillation Spectroscopic Survey (BOSS) over the redshift range of 0.1<z<0.8 for galaxies of $M_*>10^{11}M_\odot$. We find the sample completeness for BOSS is over 80% at z<0.6, but decreases at higher redshifts to about 30%. After taking these completeness factors into account, we provide accurate measurements of the stellar mass functions for galaxies with $10^{11}M_\odot<M_*<10^{12}M_\odot$, as well as the SHMRs, over the redshift range 0.1<z<0.8 in this largest galaxy redshift survey.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

A stationary and spherically symmetric black hole (For example, Reissner-Nordstrom black hole or Kerr-Newman black hole) has at most one singularity and two horizons. One horizon is the outer event horizon and the other is the inner Cauchy horizon. Can we construct static and spherically symmetric black hole solutions with N horizons and M singularities? De Sitter cosmos has only one apparent horizon. Can we construct cosmos solutions with N horizons? In this article, we present the static and spherically symmetric black hole and cosmos solutions with N horizons and M singularities in the vector-tensor theories. Following these motivations, we also construct the black hole solutions with a firewall. The deviation of these black hole solutions from the usual ones can be potentially tested by future measurements of gravitational waves.

In this paper, we study the strong gravitational lensing of gravitational waves (GWs) from a statistical perspective, with particular focus on the high frequency GWs from stellar binary black hole coalescences. These are most promising targets for ground-based detectors such as Advanced Laser Interferometer Gravitational Wave Observatory (aLIGO) and the proposed Einstein Telescope (ET) and can be safely treated under the geometrical optics limit for GW propagation. We perform a thorough calculation of the lensing rate, by taking account of effects caused by the ellipticity of lensing galaxies, lens environments, and magnification bias. We find that in certain GW source rate scenarios, we should be able to observe strongly lensed GW events once per year ($\sim1~\text{yr}^{-1}$) in the aLIGO survey at its design sensitivity; for the proposed ET survey, the rate could be as high as $\sim80~\text{yr}^{-1}$. These results depend on the estimate of GW source abundance, and hence can be correspondingly modified with an improvement in our understanding of the merger rate of stellar binary black holes. We also compute the fraction of four-image lens systems in each survey, predicting it to be $\sim30$ per cent for the aLIGO survey and $\sim6$ per cent for the ET survey. Finally, we evaluate the possibility of missing some images due to the finite survey duration, by presenting the probability distribution of lensing time delays. We predict that this selection bias will be insignificant in future GW surveys, as most of the lens systems ($\sim90$ per cent) will have time delays less than $\sim1$ month, which will be far shorter than survey durations.

The redshifted 21~cm signal from neutral hydrogen (HI) is potentially a very powerful probe for cosmology, but a difficulty in its observation is that it is much weaker than foreground radiation from the Milky Way as well as extragalactic radio sources. The foreground radiation at different frequencies are however coherent along one line of sight, and various methods of foreground subtraction based on this property have been proposed. In this paper, we present a new method based on the Robust Principal Component Analysis (RPCA) to subtract foreground and extract 21~cm signal, which explicitly uses both the low-rank property of the frequency covariance matrix (i.e. frequency coherence) of the foreground and the sparsity of the frequency covariance matrix of the 21~cm signal. The low-rank property of the foregrounds frequency covariance has been exploited in many previous works on foreground subtraction, but to our knowledge the sparsity of the frequency covariance of the 21~cm signal is first explored here. By exploiting both properties in the RPCA method, in principle, the foreground and signal may be separated without the signal loss problem. Our method is applicable to both small patch of sky with the flat-sky approximation, and to large area of sky where the sphericity has to be considered. It is also easy to be extended to deal with more complex conditions such as sky map with defects.

We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $\Sigma_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $\Sigma_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $\Sigma_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $\Sigma_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve.

The Hawking-Penrose singularity theorem states that a singularity forms inside a black hole in general relativity. To remove this singularity one must resort to a more fundamental theory. Using a corrected dynamical equation arising in loop quantum cosmology and braneworld models, we study the gravitational collapse of a perfect fluid sphere with a rather general equation of state. In the frame of an observer comoving with this fluid, the sphere pulsates between a maximum and a minimum size, avoiding the singularity. The exterior geometry is also constructed. There are usually an outer and an inner apparent horizon, resembling the Reissner-Nordstr\"om situation. For a distant observer the {horizon} crossing occurs in an infinite time and the pulsations of the black hole quantum "beating heart" are completely unobservable. However, it may be observable if the black hole is not spherical symmetric and radiates gravitational wave due to the quadrupole moment, if any.

Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged]

Understanding the host galaxy properties of stellar binary black hole (SBBH) mergers is important for revealing the origin of the SBBH gravitational-wave sources detected by advanced LIGO and helpful for identifying their electromagnetic counterparts. Here we present a comprehensive analysis of the host galaxy properties of SBBHs by implementing semi-analytical recipes for SBBH formation and merger into cosmological galaxy formation model. If the time delay between SBBH formation and merger ranges from $\la$\,Gyr to the Hubble time, SBBH mergers at redshift $z\la0.3$ occur preferentially in big galaxies with stellar mass $M_*\ga2\times10^{10}\msun$ and metallicities $Z$ peaking at $\sim0.6Z_\odot$. However, the host galaxy stellar mass distribution of heavy SBBH mergers ($M_{\bullet\bullet}\ga50\msun$) is bimodal with one peak at $\sim10^9\msun$ and the other peak at $\sim2\times10^{10}\msun$. The contribution fraction from host galaxies with $Z\la0.2Z_\odot$ to heavy mergers is much larger than that to less heavy mergers. If SBBHs were formed in the early universe (e.g., $z>6$), their mergers detected at $z\la0.3$ occur preferentially in even more massive galaxies with $M_*>3\times10^{10}\msun$ and in galaxies with metallicities mostly $\ga0.2Z_\odot$ and peaking at $Z\sim0.6Z_\odot$, due to later cosmic assembly and enrichment of their host galaxies. SBBH mergers at $z\la0.3$ mainly occur in spiral galaxies, but the fraction of SBBH mergers occur in elliptical galaxies can be significant if those SBBHs were formed in the early universe; and about two thirds of those mergers occur in the central galaxies of dark matter halos. We also present results on the host galaxy properties of SBBH mergers at higher redshift.

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.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

We apply a halo-based group finder to four large redshift surveys, the 2MRS, 6dFGS, SDSS and 2dFGRS, to construct group catalogs in the low-redshift Universe. The group finder is based on that of Yang et al. but with an improved halo mass assignment so that it can be applied uniformly to various redshift surveys of galaxies. Halo masses are assigned to groups according to proxies based on the stellar mass/luminosity of member galaxies. The performances of the group finder in grouping galaxies according to common halos and in halo mass assignments are tested using realistic mock samples constructed from hydrodynamical simulations and empirical models of galaxy occupation in dark matter halos. Our group finder finds $\sim 94\%$ of the correct true member galaxies for $90-95\%$ of the groups in the mock samples; the halo masses assigned by the group finder are un-biased with respect to the true halo masses, and have a typical uncertainty of $\sim0.2\,{\rm dex}$. The properties of group catalogs constructed from the observational samples are described and compared with other similar catalogs in the literature.

One of the major unsolved questions on the understanding of the AGN population is the origin of the dichotomy between radio-quiet and radio-loud quasars. The most promising explanation is provided by the spin paradigm, which suggests radio-loud quasars have higher black hole spin. However, the measurement of black hole spin remains extremely challenging. We here aim at comparing the mean radiative efficiencies of carefully matched samples of radio-loud and radio-quiet SDSS quasars at 0.3<z<0.8. We use the [OIII] luminosity as an indirect average tracer of the ionizing continuum in the extreme-UV regime where differences in the SED due to black hole spin are most pronounced. We find that the radio-loud sample shows an enhancement in [OIII] line strength by a factor of at least 1.5 compared to a radio-quiet sample matched in redshift, black hole mass and optical continuum luminosity or accretion rate. We argue that this enhancement is caused by differences in the SED, suggesting higher average bolometric luminosities at fixed accretion rate in the radio-loud population. This suggests that the radio-loud quasar population has on average systematically larger radiative efficiencies and therefore higher black hole spin than the radio-quiet population, providing observational support for the black hole spin paradigm.

We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ($ M_r < -12.3 $). We define a Milky Way analog based on $K$-band luminosity and local environment. Here, we present satellite luminosity functions for 8 Milky Way analog galaxies between 20 to 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to $ r_o < 20.75 $ using low redshift $gri$ color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky Way analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are five Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

Matching the number counts and redshift distribution of submillimeter galaxies (SMGs) without invoking modifications to the initial mass function (IMF) has proved challenging for semi-analytic models (SAMs) of galaxy formation. We adopt a previously developed SAM that is constrained to match the z=0 galaxy stellar mass function and makes various predictions that agree well with observational constraints; we do not recalibrate the SAM for this work. We implement three prescriptions to predict the submillimeter flux densities of the model galaxies; two depend solely on star formation rate, whereas the other also depends on the dust mass. By comparing the predictions of the models, we find that taking into account the dust mass, which affects the dust temperature and thus influences the far-infrared spectral energy distribution, is crucial for matching the number counts and redshift distribution of SMGs. Moreover, despite using a standard IMF, our model can match the observed SMG number counts and redshift distribution reasonably well, which contradicts the conclusions of some previous studies that a top-heavy IMF, in addition to taking into account the effect of dust mass, is needed to match these observations. Although we have not identified the key ingredient that is responsible for our model matching the observed SMG number counts and redshift distribution without IMF variation- which is challenging given the different prescriptions for physical processes employed in the SAMs of interest-our results demonstrate that in SAMs, IMF variation is degenerate with other physical processes, such as stellar feedback.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

We study the X-ray phases of the Crab pulsar utilizing the 11-year observations from the Rossi X-ray Timing Explorer, 6-year radio observations from the Nanshan Telescope, and the ephemeris from Jodrell Bank Observatory. It is found that the X-ray phases in different energy bands and the radio phases from Nanshan Telescope show similar behaviors, including long-time evolution and short-time variations. Such strong correlations between the X-ray and radio phases imply that the radio and X-ray timing noises are both generated from the pulsar spin that cannot be well described by the the monthly ephemeris from the Jodrell Bank observatory. When using the Nanshan phases as references to study the X-ray timing noise, it has a significantly smaller variation amplitude and shows no long-time evolution, with a change rate of $(-1.1\pm1.1)\times10^{-7}$ periods per day. These results show that the distance of the X-ray and radio emission regions on the Crab pulsar has no detectable secular change, and it is unlikely that timing-noises resulted from any unique physical processes in the radio or X-ray emitting regions. The similar behaviors of the X-ray and radio timing noises also imply that the variation of the interstellar medium is not the origin of the Crab pulsar's timing noises, which is consistent with the results obtained from the multi-frequency radio observations of PSR B1540-06.

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.

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy connection exists beyond the local universe. To this purpose, we only need a statistically representative, rather than a complete, sample of low-mass galaxies, which enables our study to $z\gtrsim1.5$. For each low-mass galaxy, we measure the projected distance ($d_{proj}$) to its nearest massive neighbor (M*$>10^{10.5} M_\odot$) within a redshift range. At a given redshift and M*, the environmental quenching effect is considered to be observed if the $d_{proj}$ distribution of QGs ($d_{proj}^Q$) is significantly skewed toward lower values than that of star-forming galaxies ($d_{proj}^{SF}$). For galaxies with $10^{8} M_\odot < M* < 10^{10} M_\odot$, such a difference between $d_{proj}^Q$ and $d_{proj}^{SF}$ is detected up to $z\sim1$. Also, about 10\% of the quenched galaxies in our sample are located between two and four virial radii ($R_{Vir}$) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, $d_{proj}^Q / R_{Vir}$, decreases with satellite M* at $M* \lesssim 10^{9.5} M_\odot$, but increases with satellite M* at $M* \gtrsim 10^{9.5} M_\odot$. This trend suggests a smooth, if any, transition of the quenching timescale around $M* \sim 10^{9.5} M_\odot$ at $0.5<z<1.0$.

Radio-loud high-redshift quasars (HRQs), although only a few of them are known to date, are crucial for the studies of the growth of supermassive black holes (SMBHs) and the evolution of active galactic nuclei (AGN) at early cosmological epochs. Radio jets offer direct evidence of SMBHs, and their radio structures can be studied with the highest angular resolution using Very Long Baseline Interferometry (VLBI). Here we report on the observations of three HRQs (J0131-0321, J0906+6930, J1026+2542) at z>5 using the Korean VLBI Network and VLBI Exploration of Radio Astrometry Arrays (together known as KaVA) with the purpose of studying their pc-scale jet properties. The observations were carried out at 22 and 43 GHz in 2016 January among the first-batch open-use experiments of KaVA. The quasar J0906+6930 was detected at 22 GHz but not at 43 GHz. The other two sources were not detected and upper limits to their compact radio emission are given. Archival VLBI imaging data and single-dish 15-GHz monitoring light curve of J0906+6930 were also acquired as complementary information. J0906+6930 shows a moderate-level variability at 15 GHz. The radio image is characterized by a core-jet structure with a total detectable size of ~5 pc in projection. The brightness temperature, 1.9x10^{11} K, indicates relativistic beaming of the jet. The radio properties of J0906+6930 are consistent with a blazar. Follow-up VLBI observations will be helpful for determining its structural variation.

We perform the first spatially-resolved stellar population study of galaxies in the early universe (z = 3.5 - 6.5), utilizing the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) imaging dataset over the GOODS-S field. We select a sample of 418 bright and extended galaxies at z = 3.5 - 6.5 from a parent sample of ~ 8000 photometric-redshift selected galaxies from Finkelstein et al. (2015). We first examine galaxies at 3.5< z < 4.0 using additional deep K-band survey data from the HAWK-I UDS and GOODS Survey (HUGS) which covers the 4000A break at these redshifts. We measure the stellar mass, star formation rate, and dust extinction for galaxy inner and outer regions via spatially-resolved spectral energy distribution fitting based on a Markov Chain Monte Carlo algorithm. By comparing specific star formation rates (sSFRs) between inner and outer parts of the galaxies we find that the majority of galaxies with the high central mass densities show evidence for a preferentially lower sSFR in their centers than in their outer regions, indicative of reduced sSFRs in their central regions. We also study galaxies at z ~ 5 and 6 (here limited to high spatial resolution in the rest-frame ultraviolet only), finding that they show sSFRs which are generally independent of radial distance from the center of the galaxies. This indicates that stars are formed uniformly at all radii in massive galaxies at z ~ 5 - 6, contrary to massive galaxies at z < 4.

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

In this paper, the dependence of the broad Fe K$\alpha$ line on the physical parameters of AGN, such as the black hole mass $M_{BH}$, accretion rate (equivalently represented by Eddington ratio $\lambda_{Edd}$), and optical classification, is investigated by applying the X-ray spectra stacking method to a large sample of AGN which have well measured optical parameters. A broad line feature is detected ($>3\sigma$) in the stacked spectra of the high $\lambda_{Edd}$ sub-sample ($\log\lambda_{Edd}>-0.9$). The profile of the broad line can be well fitted with relativistic broad line model, with the line energy consistent with highly ionized Fe K$\alpha$ line (i.e. Fe xxvi). A model consisting of multiple narrow lines cannot be ruled out, however. We found hints that the Fe K line becomes broader as the $\lambda_{Edd}$ increases. No broad line feature is shown in the sub-sample of broad-line Seyfert 1 (BLS1) galaxies and in the full sample, while a broad line might be present, though at low significance, in the sub-sample of narrow-line Seyfert 1 (NLS1) galaxies. We find no strong dependence of the broad line on black hole masses. Our results indicate that the detection/properties of the broad Fe K$\alpha$ line may strongly depend on $\lambda_{Edd}$, which can be explained if the ionization state and/or truncation radius of the accretion disc changes with $\lambda_{Edd}$. The non-detection of the broad line in the BLS1 sub-sample can be explained if the the average EW of the relativistic Fe K$\alpha$ line is weak or/and the fraction of sources with relativistic Fe K$\alpha$ line is small in BLS1 galaxies.

Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

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.

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 \AA\ ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $\alpha=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $\alpha=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

It is widely believed that the large discrepancy between the observed number of satellite galaxies and the predicted number of dark subhalos can be resolved via a variety of baryonic effects which suppress star formation in low mass halos.Supporting this hypothesis, numerous high resolution simulations with star formation, and associated feedback have been shown to reproduce the satellite luminosity function around Milky Way-mass simulated galaxies at redshift zero. However, a more stringent test of these models is their ability to simultaneously match the satellite luminosity functions of a range of host halo masses and redshifts. In this work we measure the luminosity function of faint (sub-Small Magellanic Cloud luminosity) satellites around hosts with stellar masses 10.5$<\log_{10}$M$_*$/M$_\odot<11.5$ to an unprecedented redshift of 1.5. This new measurement of the satellite luminosity function provides powerful new constraining power; we compare these results with predictions from four different simulations and show that although the models perform similarly over-all, no one model reproduces the satellite luminosity function reliably at all redshifts and host stellar masses. This result highlights the continued need for improvement in understanding the fundamental physics that governs satellite galaxy evolution.

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.

We propose a geometrically thick, super-Eddington accretion disk model, where an optically thick wind is not necessary, to understand ultraluminous supersoft sources (ULSs). For high mass accretion rates $\dot M \ga 30\dot M_{\rm Edd}$ and not small inclination angles $\theta \ga 25^{\circ}$, where $\dot M_{\rm Edd}$ is the Eddington accretion rate, the hard photons from the hot inner region may be shaded by the geometrically thick inner disk, and therefore only the soft photons from the outer thin disk and the outer photosphere of the thick disk can reach the observer. Our model can naturally explain the approximate relation between the typical thermal radius and the thermal temperature, $R_{\rm bb} \propto T_{\rm bb}^{-2}$. Moreover, the thick disk model can unify ULSs and normal ultraluminous X-ray sources, where the different observational characteristics are probably related to the inclination angle and the mass accretion rate. By comparing our model with the optically thick outflow model, we find that less mass accretion rate is required in our model.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

The formation of relativistic jets by an accreting compact object is one of the fundamental mysteries of astrophysics. While the theory is poorly understood, observations of relativistic jets from systems known as microquasars have led to a well-established phenomenology. Relativistic jets are not expected from sources with soft or supersoft X-ray spectra, although two such systems are known to produce relatively low-velocity bipolar outflows. Here we report optical spectra of an ultraluminous supersoft X-ray source (ULS) in the nearby galaxy M81 (M81 ULS-1) showing blueshifted broad H\alpha\ emission lines, characteristic of baryonic jets with relativistic speeds. The time variable jets have projected velocities ~17 per cent of the speed of light, and seem similar to those in the prototype microquasar SS 433. Such relativistic jets are not expected to be launched from white dwarfs, but an origin from a black hole or neutron star in M81 ULS-1 is hard to reconcile with its constant soft X-rays. The completely unexpected presence of relativistic jets in a ULS challenges the canonical theories for jet formation, but may possibly be explained by a long speculated super-critically accreting black hole with optically thick outflows

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.

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3\sigma$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

Massive luminous red galaxies (LRGs) are believed to be evolving passively and can be used as cosmic chronometers to estimate the Hubble constant. However, different LRGs may be located in different environments. The environmental effects may limit the use of the LRGs as cosmic chronometers. We aim to investigate the environmental and mass dependence of the formation of 'quiescent' LRGs selected from the Sloan Digital Sky Survey Date Release 8 and to pave the way for using LRGs as cosmic chronometers. Using the population synthesis software STARLIGHT, we derive the stellar populations in each LRG through the full spectrum fitting and obtain the mean age distribution and the mean star formation history (SFH) of those LRGs. We find that there is no apparent dependence of the mean age and the SFH of quiescent LRGs on their environment, while the ages of those quiescent LRGs depend weakly on their mass. We compare the SFHs of the SDSS LRGs with those obtained from a semi-analytical galaxy formation model, and find that they are roughly consistent with each other if we consider the errors in the STARLIGHT-derived ages. We find that a small fraction of later star formation in LRGs leads to a systematical overestimation (~28 %) of the Hubble constant by the differential age method, and the systematical errors in the STARLIGHT-derived ages may lead to an underestimation (~ 16 %) of the Hubble constant. However, these errors can be corrected by a detailed study of the mean SFH of those LRGs and by calibrating the STARLIGHT-derived ages to those obtained independently by other methods. The environmental effects do not play significant role in the age estimates of quiescent LRGs, and the quiescent LRGs as a population can be used securely as cosmic chronometers.

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

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).

A scheme for estimating atmospheric parameters T$_{eff}$, log$~g$, and [Fe/H] is proposed on the basis of Least Absolute Shrinkage and Selection Operator (LASSO) algorithm and Haar wavelet. The proposed scheme consists of three processes. A spectrum is decomposed using the Haar wavelet transform and low-frequency components at the fourth level are considered as candidate features. Then, spectral features from the candidate features are detected using the LASSO algorithm to estimate the atmospheric parameters. Finally, atmospheric parameters are estimated from the extracted spectral features using the support-vector regression (SVR) method. The proposed scheme was evaluated using three sets of stellar spectra respectively from Sloan Digital Sky Survey (SDSS), Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST), and Kurucz's model, respectively. The mean absolute errors are as follows: for 40~000 SDSS spectra, 0.0062 dex for log~T$_{eff}$ (85.83 K for T$_{eff}$), 0.2035 dex for log$~g$ and 0.1512 dex for [Fe/H]; for 23963 LAMOST spectra, 0.0074 dex for log~T$_{eff}$ (95.37 K for T$_{eff}$), 0.1528 dex for log~$g$, and 0.1146 dex for [Fe/H]; and for 10469 synthetic spectra, 0.0010 dex for log T$_{eff}$(14.42K for T$_{eff}$), 0.0123 dex for log~$g$, and 0.0125 dex for [Fe/H].

Hyper-velocity stars are believed to be ejected out from the Galactic center through dynamical interactions between (binary) stars and the central massive black hole(s). In this paper, we report 19 low mass F/G/K type hyper-velocity star candidates from over one mil- lion stars of the first data release of the LAMOST general survey. We determine the unbound probability for each candidate using a Monte-Carlo simulation by assuming a non-Gaussian proper-motion error distribution, Gaussian heliocentric distance and radial velocity error dis- tributions. The simulation results show that all the candidates have unbound possibilities over 50% as expected, and one of them may even exceed escape velocity with over 90% probabili- ty. In addition, we compare the metallicities of our candidates with the metallicity distribution functions of the Galactic bulge, disk, halo and globular cluster, and conclude that the Galactic bulge or disk is likely the birth place for our candidates.

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

Large-scale and deep sky survey missions are rapidly collecting a large amount of stellar spectra, which necessitate the estimation of atmospheric parameters directly from spectra and makes it feasible to statistically investigate latent principles in a large dataset. We present a technique for estimating parameters $T_{eff}$, log$~g$ and [Fe/H] from stellar spectra. With this technique, we first extract features from stellar spectra using the LASSO algorithm; then, the parameters are estimated from the extracted features using the SVR. On a subsample of 20~000 stellar spectra from SDSS with reference parameters provided by SDSS/SEGUE Pipeline SSPP, estimation consistency are 0.007458 dex for log$~T_{eff}$ (101.609921 K for $T_{eff}$), 0.189557 dex for log$~g$ and 0.182060 for [Fe/H], where the consistency is evaluated by mean absolute error. Prominent characteristics of the proposed scheme are sparseness, locality, and physical interpretability. In this work, every spectrum consists of 3821 fluxes, and 10, 19, and 14 typical wavelength positions are detected respectively for estimating $T_{eff}$, log$~g$ and [Fe/H]. It is shown that the positions are related to typical lines of stellar spectra. This characteristic is important in investigating physical indications from analysis results. Then, stellar spectra can be described by the individual fluxes on the detected positions (PD) or local integration of fluxes near them (LI). The abovementioned consistency is the result based on features described by LI. If features are described by PD, consistency are 0.009092 dex for log$~T_{eff}$ (124.545075 K for $T_{eff}$), 0.198928 dex for log$~g$, and 0.206814 dex for [Fe/H].

We describe a scheme to extract linearly supporting (LSU) features from stellar spectra to automatically estimate the atmospheric parameters $T_{eff}$, log$~g$, and [Fe/H]. "Linearly supporting" means that the atmospheric parameters can be accurately estimated from the extracted features through a linear model. The successive steps of the process are as follow: first, decompose the spectrum using a wavelet packet (WP) and represent it by the derived decomposition coefficients; second, detect representative spectral features from the decomposition coefficients using the proposed method Least Absolute Shrinkage and Selection Operator (LARS)$_{bs}$; third, estimate the atmospheric parameters $T_{eff}$, log$~g$, and [Fe/H] from the detected features using a linear regression method. One prominent characteristic of this scheme is its ability to evaluate quantitatively the contribution of each detected feature to the atmospheric parameter estimate and also to trace back the physical significance of that feature. This work also shows that the usefulness of a component depends on both wavelength and frequency. The proposed scheme has been evaluated on both real spectra from the Sloan Digital Sky Survey (SDSS)/SEGUE and synthetic spectra calculated from Kurucz's NEWODF models. On real spectra, we extracted 23 features to estimate $T_{eff}$, 62 features for log$~g$, and 68 features for [Fe/H]. Test consistencies between our estimates and those provided by the Spectroscopic Sarameter Pipeline of SDSS show that the mean absolute errors (MAEs) are 0.0062 dex for log$~T_{eff}$ (83 K for $T_{eff}$), 0.2345 dex for log$~g$, and 0.1564 dex for [Fe/H]. For the synthetic spectra, the MAE test accuracies are 0.0022 dex for log$~T_{eff}$ (32 K for $T_{eff}$), 0.0337 dex for log$~g$, and 0.0268 dex for [Fe/H].

Using a semi-analytical approach we investigate the characteristics of predictions for the masses and metallicities of the baryonic matter in and around galaxies made by three galaxy formation models. These models represent three different feedback scenarios: one model with purely ejective feedback, one model with ejective feedback with reincorporation of ejected gas, and one preventative model. We find that, when the model parameters are adjusted to predict the correct stellar masses for a range of halo masses between 10^{10} to 10^{12}Msun, these three scenarios have very different predictions for the masses and metallicities of the interstellar and circum-galactic media. Compared with current observational data, the model implementing preventative feedback has a large freedom to match a broad range of observational data, while the ejective models have difficulties to match a number of observational constraints simultaneously, independent of how the ejection and reincorporation are implemented. Our results suggest that the feedback process which regulates the amounts of stars and cold gas in low-mass galaxies is preventative in nature.

We use a set of observational data for galaxy cold gas mass fraction and gas phase metallicity to constrain the content, inflow and outflow of gas in central galaxies hosted by halos with masses between $10^{11} M_{\odot}$ to $10^{12} M_{\odot}$. The gas contents in high redshift galaxies are obtained by combining the empirical star formation histories of Lu et al. (2014) and star formation models that relate star formation rate with the cold gas mass in galaxies. We find that the total baryon mass in low-mass galaxies is always much less than the universal baryon mass fraction since $z = 2$, regardless of star formation model adopted. The data for the evolution of the gas phase metallicity require net metal outflow at $z\lesssim 2$, and the metal loading factor is constrained to be about $0.01$, or about $60\%$ of the metal yield. Based on the assumption that galactic outflow is more enriched in metal than both the interstellar medium and the material ejected at earlier epochs, we are able to put stringent constraints on the upper limits for both the net accretion rate and the net mass outflow rate. The upper limits strongly suggest that the evolution of the gas phase metallicity and gas mass fraction for low-mass galaxies at $z < 2$ is not compatible with strong outflow. We speculate that the low star formation efficiency of low-mass galaxies is owing to some preventative processes that prevent gas from accreting into galaxies in the first place.

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

Reliable halo mass estimation for a given galaxy system plays an important role both in cosmology and galaxy formation studies. Here we set out to find the way that can improve the halo mass estimation for those galaxy systems with limited brightest member galaxies been observed. Using four mock galaxy samples constructed from semi-analytical formation models, the subhalo abundance matching method and the conditional luminosity functions, respectively, we find that the luminosity gap between the brightest and the subsequent brightest member galaxies in a halo (group) can be used to significantly reduce the scatter in the halo mass estimation based on the luminosity of the brightest galaxy alone. Tests show that these corrections can significantly reduce the scatter in the halo mass estimations by $\sim 50\%$ to $\sim 70\%$ in massive halos depending on which member galaxies are considered. Comparing to the traditional ranking method, we find that this method works better for groups with less than five members, or in observations with very bright magnitude cut.

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

While a broad profile of the Fe K$\alpha$ emission line is frequently found in the X-ray spectra of typical Seyfert galaxies, the situation is unclear in the case of Narrow Line Seyfert 1 galaxies (NLS1s)---an extreme subset which are generally thought to harbor less massive black holes with higher accretion rates. In this paper, the ensemble property of the Fe K$\alpha$ line in NLS1s is investigated by stacking the X-ray spectra of a large sample of 51 NLS1s observed with {\it XMM-Newton}. The composite X-ray spectrum reveals a prominent, broad emission feature over 4--7 keV, characteristic of the broad Fe K$\alpha$ line. In addition, there is an indication for a possible superimposing narrow (unresolved) line, either emission or absorption, corresponding to Fe XXVI or Fe XXV, respectively. The profile of the broad emission feature can well be fitted with relativistic broad-line models, with the line energy consistent either with 6.4 keV (i.e., neutral Fe) or with 6.67 keV (i.e., highly ionized Fe), in the case of the narrow line being emission and absorption, respectively. Interestingly, there are tentative indications for low or intermediate values of the average spins of the black holes ($a<0.84$), as inferred from the profile of the composite broad line. If the observed feature is indeed a broad line rather than resulting from partial covering absorption, our results suggest that a relativistic Fe line may in fact be common in NLS1s; and there are tentative indications that black holes in NLS1s may not spin very fast in general.

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

A large inflationary tensor-to-scalar ratio $r_\mathrm{0.002} = 0.20^{+0.07}_{-0.05}$ is reported by the BICEP2 team based on their B-mode polarization detection, which is outside of the $95\%$ confidence level of the Planck best fit model. We explore several possible ways to reduce the tension between the two by considering a model in which $\alpha_\mathrm{s}$, $n_\mathrm{t}$, $n_\mathrm{s}$ and the neutrino parameters $N_\mathrm{eff}$ and $\Sigma m_\mathrm{\nu}$ are set as free parameters. Using the Markov Chain Monte Carlo (MCMC) technique to survey the complete parameter space with and without the BICEP2 data, we find that the resulting constraints on $r_\mathrm{0.002}$ are consistent with each other and the apparent tension seems to be relaxed. Further detailed investigations on those fittings suggest that $N_\mathrm{eff}$ probably plays the most important role in reducing the tension. We also find that the results obtained from fitting without adopting the consistency relation do not deviate much from the consistency relation. With available Planck, WMAP, BICEP2 and BAO datasets all together, we obtain $r_{0.002} = 0.14_{-0.11}^{+0.05}$, $n_\mathrm{t} = 0.35_{-0.47}^{+0.28}$, $n_\mathrm{s}=0.98_{-0.02}^{+0.02}$, and $\alpha_\mathrm{s}=-0.0086_{-0.0189}^{+0.0148}$; if the consistency relation is adopted, we get $r_{0.002} = 0.22_{-0.06}^{+0.05}$.

We explore a cyclic universe due to phantom and quintessence fields. We find that, in every cycle of the evolution of the universe, the phantom dominates the cosmic early history and quintessence dominates the cosmic far future. In this model of universe, there are infinite cycles of expansion and contraction. Different from the inflationary universe, the corresponding cosmic space-time is geodesically complete and quantum stable. But similar to the Cyclic Model, the flatness problem, the horizon problem and the large scale structure of the universe can be explained in this cyclic universe.

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

We present a new technique for overcoming confusion noise in deep far-infrared \Herschel space telescope images making use of prior information from shorter $\lambda<2$\micron wavelengths. For the deepest images obtained by \Herschels, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated \Herschel PACS-160\micron images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict FIR flux of sources based on ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of $3\sigma_{c} $=2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For \Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$\sigma$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

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.

Combining the latest Planck, Wilkinson Microwave Anisotropy Probe (WMAP), and baryon acoustic oscillation (BAO) data, we exploit the recent cosmic microwave background (CMB) B-mode power spectra data released by the BICEP2 collaboration to constrain the cosmological parameters of the $\Lambda$CDM model, esp. the primordial power spectra parameters of the scalar and the tensor modes, $n_s,\alpha_s, r, n_t$. We obtained constraints on the parameters for a lensed $\Lambda$CDM model using the Markov Chain Monte Carlo (MCMC) technique, the marginalized $68\%$ bounds are $r = 0.1043\pm_{0.0914}^{0.0307}$, $n_s = 0.9617\pm_{0.0061}^{ 0.0061}$, $\alpha_s = -0.0175\pm_{0.0097}^{0.0105}$, $n_t = 0.5198\pm_{0.4579}^{0.4515}$. We found that a blue tilt for $n_t$ is favored slightly, but it is still well consistent with flat or even red tilt. Our $r$ value is slightly smaller than the one obtained by the BICEP team, as we allow $n_t$ as a free parameter without imposing the single-field slow roll inflation consistency relation. If we impose this relation, $r=0.2130\pm_{0.0609}^{0.0446}$. For most other parameters, the best fit values and measurement errors are not altered much by the introduction of the BICEP2 data.

We fit the spectral energy distributions (SEDs) of a GeV-TeV FSRQ sample with the leptonic model. Their gamma_min of the relativistic electron distributions, which significantly affect the estimates of the jet properties, are constrained, with a typical value of 48. Their jet power, magnetized parameter, radiation efficiency, and jet production/radiation rates per central black hole (BH) mass are derived and compared to that of BL Lacs. We show that he FSRQ jets may be dominated by the Poynting flux and have a high radiation efficiency, whereas the BL Lac jets are likely dominated by particles and have a lower radiation efficiency than FSRQs. Being different from BL Lacs, the jet powers of FSRQs are proportional to their central BH masses. The jet production and radiation rates of the FSRQs distribute in narrow ranges and are correlated with each other, whereas no similar feature is found for the BL Lacs. We also show that the jet power is correlated with the cavity kinetic power, and the magnetic field energy in the jets may provide the cavity kinetic energy of FSRQs and the kinetic energy of cold protons in the jets may be crucial for cavity kinetic energy of BL Lacs. We suggest that the dominating formation mechanism of FSRQ jets may be the BZ process, but BL Lac jets may be produced via the BP and/or BZ processes, depending on the structures and accretion rates of accretion disks.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

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.

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.

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

We estimate the radiative efficiency (epsilon) of individual type 1 SDSS QSOs by using their bolometric luminosities (L_{bol}) and accretion rates (\dot{M}_{\bullet,acc}), which may be related to the assembly histories and spins of the central massive black holes (MBHs). We estimate L_{bol} by using the empirical spectral energy distributions of QSOs and \dot{M}_{\bullet,acc} by fitting the observed optical luminosity(/-ies) with the thin accretion disk model, assuming the MBH masses given by the virial mass estimator(s) (M_{\bullet,vir}). We find an apparent correlation between epsilon and M_{\bullet,vir}, which is strong at redshift z \la 1.8, weak at z \ga 2, and consistent with that found by Davis & Laor (2011) for 80 PG QSOs at z \leq 0.5. To investigate whether this correlation is intrinsic or not, we construct a mock sample of QSOs according to the true MBH mass and Eddington ratio distributions given in Kelly & Shen (2013). By comparing the results obtained from the mock sample with that from the SDSS sample, we demonstrate that the apparent epsilon-M_{\bullet,vir} correlation can be produced by and mainly due to the selection effects of the SDSS sample and the bias induced by the usage of M_{\bullet,vir} as the true MBH mass. The mean values of epsilon of those SDSS QSOs are consistent with being a constant \simeq 0.11-0.16 over the redshift range of 0.3 \la z \la 4. We conclude that the current SDSS QSO data is consistent with no strong intrinsic correlation between radiative efficiency and true MBH mass and no significant redshift evolution of radiative efficiencies.

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.

We estimate the radiative efficiency (epsilon) of individual type 1 SDSS QSOs by using their bolometric luminosities (L_{bol}) and accretion rates (\dot{M}_{\bullet,acc}), which may be related to the assembly histories and spins of the central massive black holes (MBHs). We estimate L_{bol} by using the empirical spectral energy distributions of QSOs and \dot{M}_{\bullet,acc} by fitting the observed optical luminosity(/-ies) with the thin accretion disk model, assuming the MBH masses given by the virial mass estimator(s) (M_{\bullet,vir}). We find an apparent correlation between epsilon and M_{\bullet,vir}, which is strong at redshift z \la 1.8, weak at z \ga 2, and consistent with that found by Davis & Laor (2011) for 80 PG QSOs at z \leq 0.5. To investigate whether this correlation is intrinsic or not, we construct a mock sample of QSOs according to the true MBH mass and Eddington ratio distributions given in Kelly & Shen (2013). By comparing the results obtained from the mock sample with that from the SDSS sample, we demonstrate that the apparent epsilon-M_{\bullet,vir} correlation can be produced by and mainly due to the selection effects of the SDSS sample and the bias induced by the usage of M_{\bullet,vir} as the true MBH mass. The mean values of epsilon of those SDSS QSOs are consistent with being a constant \simeq 0.11-0.16 over the redshift range of 0.3 \la z \la 4. We conclude that the current SDSS QSO data is consistent with no strong intrinsic correlation between radiative efficiency and true MBH mass and no significant redshift evolution of radiative efficiencies.

Many schemes have been proposed to perform a model-independent constraint on cosmological dynamics, such as nonparametric dark energy equation of state (EoS) \omega(z) or the deceleration parameter q(z). These methods usually contain derivative processes with respect to observational data with noise. However, it still remains remarkably uncertain when one estimates the numerical differentiation, especially the corresponding truncation errors. In this work, we introduce a global numerical differentiation method, first formulated by Reinsch(1967), which is smoothed by cubic spline functions. The optimal solution is obtained by minimizing the functional \Phi(f). To investigate the potential of the algorithm further, we apply it to the estimation of the transition redshift z_{t} with simulated expansion rate E(z) based on observational Hubble parameter data(OHD). An effective method to determine the free parameter S appearing in Reinsch Splines is provided.

Gamma-ray bursts (GRBs) and GeV-TeV selected radio loud Active Galactic Nuclei (AGNs) are compared based on our systematic modeling of the observed spectral energy distributions of a sample of AGNs with a single-zone leptonic model. We show that the correlation between the jet power (P_{jet}) and the prompt gamma-ray luminosity (L_{jet}) of GRBs is consistent, within the uncertainties, with the correlation between jet power and the synchrotron peak luminosity (L_{s, jet}) of flat spectrum radio quasars (FSRQs). Their radiation efficiencies (\varepsilon) are also comparable (>10% for most sources), which increase with the bolometric jet luminosity (L_{bol,jet}) for FSRQs and with the L_{jet} for GRBs with similar power-law indices. BL Lacs do not follow the P_{jet}-L_{s, jet} relation of FSRQs. They have lower \varepsilon and L_{bol, jet} values than FSRQs, and a tentative L_{bol, jet}-\varepsilon relation is also found, with a power-law index being different from that of the FSRQs. The magnetization parameters (\sigma) of FSRQs are averagely larger than that of BL Lacs. They are anti-correlated with $\varepsilon$ for the FSRQs, but positive correlated with \varepsilon for the BL Lacs. GeV Narrow-line Seyfert 1 galaxies potentially share similar properties with FSRQs. Based on the analogy between GRBs and FSRQs, we suggest that the prompt gamma-ray emission of GRBs is likely produced by synchrotron process in a magnetized jet with high radiation efficiency, similar to FSRQs. The jets of BL Lacs, on the other hand, are less efficient and are likely more matter dominated.

Many approaches to obtaining cosmological constraints rely on the connection between galaxies and dark matter. However, the distribution of galaxies is dependent on their formation and evolution as well as the cosmological model, and galaxy formation is still not a well-constrained process. Thus, methods that probe cosmology using galaxies as a tracer for dark matter must be able to accurately estimate the cosmological parameters without knowing the details of galaxy formation a priori. We apply this reasoning to the method of obtaining $\Omega_m$ and $\sigma_8$ from galaxy clustering combined with the mass-to-number ratio of galaxy clusters. To test the sensitivity of this method to variations due to galaxy formation, we consider several different models applied to the same cosmological dark matter simulation. The cosmological parameters are then estimated using the observables in each model, marginalizing over the parameters of the Halo Occupation Distribution (HOD). We find that for models where the galaxies can be well represented by a parameterized HOD, this method can successfully extract the desired cosmological parameters for a wide range of galaxy formation prescriptions.

The standard cosmological model based on cold dark matter (CDM) predicts a large number of subhalos for each galaxy-size halo. It is well known that matching the subhalos to the observed properties of luminous satellites of galaxies in the local universe poses a significant challenge to our understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and host mass dependence of the luminosity function of satellites provides a powerful new diagnostic to disentangle astrophysical effects from variations in the underlying dark matter mass function. We illustrate this by comparing the results of recent observations of satellites out to $z=0.8$ based on Hubble Space Telescope images with the predictions of three different sets of state-of-the art semi-analytic models with underlying CDM power spectra and one semi-analytic model with an underlying Warm Dark Matter (WDM) power spectrum. We find that even though CDM models provide a reasonable fit to the local luminosity function of satellites around galaxies comparable or slightly larger than the Milky Way, they do not reproduce the data as well for different redshift and host galaxy stellar mass. This tension indicates that further improvements are likely to be needed in the description of star formation if the models are to be reconciled with the data. The WDM model matches the observed mass dependence and redshift evolution of satellite galaxies more closely than any of the CDM models, indicating that a modification of the underlying power spectrum may offer an alternative solution to this tension. We conclude by presenting predictions for the color magnitude relation of satellite galaxies to demonstrate how future observations will be able to further distinguish between these models and help constrain baryonic and non-baryonic physics.

In this paper, we select a sample of quiescent luminous red galaxies (LRGs) from Sloan Digital Sky Survey Data Release 7 (SDSS DR7) with signal-to-noise ratio (S/N) substantially high to study the consistency of the full spectrum fitting method by using different packages, mainly, ULySS and STARLIGHT. The spectrum of each galaxy in the sample is fitted by the full spectrum fitting packages ULySS and STARLIGHT, respectively. We find: (1) for spectra with higher S/Ns, the ages of stellar populations obtained from ULySS are slightly older than that from STARLIGHT, and metallicities derived from ULySS are slightly richer than that from STARLIGHT. In general, both packages can give roughly consistent fitting results. (2) for low S/N spectra, it is possible that the fitting by ULySS can be trapped some local minimum parametric regions during execution and thus may give unreliable results, while STARLIGHT can still give reliable results. Based on the fitting results of LRGs, we further analyse their star formation history (SFH) and the relation between their age and velocity dispersion, and find that those highly agrees with conclusions from others' previous work.

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

We study the evolution of the number density, as a function of the size, of passive early-type galaxies with a wide range of stellar masses 10^10<M*/Msun<10^11.5) from z~3 to z~1, exploiting the unique dataset available in the GOODS-South field, including the recently obtained WFC3 images as a part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). In particular, we select a sample of 107 massive (M*>10^10 M_sun), passive (SSFR<10^-2 Gyr^-1) and morphologically spheroidal galaxies at 1.2<z<3, taking advantage of the panchromatic dataset available for GOODS, including VLT, CFHT, Spitzer, Chandra and HST ACS+WFC3 data. We find that at 1<z<3 the passively evolving early-type galaxies are the reddest and most massive objects in the Universe, and we prove that a correlation between mass, morphology, color and star-formation activity is already in place at that epoch. We measure a significant evolution in the mass-size relation of passive early-type galaxies (ETGs) from z~3 to z~1, with galaxies growing on average by a factor of 2 in size in a 3 Gyr timescale only. We witness also an increase in the number density of passive ETGs of 50 times over the same time interval. We find that the first ETGs to form at z>2 are all compact or ultra-compact, while normal sized ETGs (meaning ETGs with sizes comparable to those of local counterparts of the same mass) are the most common ETGs only at z<1. The increase of the average size of ETGs at 0<z<1 is primarily driven by the appearance of new large ETGs rather than by the size increase of individual galaxies.

Four well-sampled spectral energy distributions (SEDs) of PMN J0948+0022 are fitted with the syn+SSC+EC model to derive the physical parameters of its jets and to investigate the spectral variations of its SEDs. A tentative correlation between the peak luminosity L_c and peak frequency nu_c of its inverse Compton (IC) bump is found in both the observer and co-moving frames, indicating that the variations of luminosity are accompanied with the spectral shift. A correlation between L_c and delta is found, and thus the magnification of the external photon field by the bulk motion of the radiation regions is an essential reason for the spectral variation since the IC bump of PMN J0948+0022 is dominated by the EC process.

We present a new submm/mm galaxy counterpart identification technique which builds on the use of Spitzer IRAC colors as discriminators between likely counterparts and the general IRAC galaxy population. Using 102 radio- and SMA-confirmed counterparts to AzTEC sources across three fields (GOODS-N, GOODS-S, and COSMOS), we develop a non-parametric IRAC color-color characteristic density distribution (CDD), which, when combined with positional uncertainty information via likelihood ratios, allows us to rank all potential IRAC counterparts around SMGs and calculate the significance of each ranking via the reliability factor. We report all robust and tentative radio counterparts to SMGs, the first such list available for AzTEC/COSMOS, as well as the highest ranked IRAC counterparts for all AzTEC SMGs in these fields as determined by our technique. We demonstrate that the technique is free of radio bias and thus applicable regardless of radio detections. For observations made with a moderate beamsize (~18"), this technique identifies ~85 per cent of SMG counterparts. For much larger beamsizes (>30"), we report identification rates of 33-49 per cent. Using simulations, we demonstrate that this technique is an improvement over using positional information alone for observations with facilities such as AzTEC on the LMT and SCUBA-2 on JCMT.

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.

The relative age of galaxies at different redshifts can be used to infer the Hubble parameter and put constraints on cosmological models. We select 23,883 quiescent luminous red galaxies (LRGs) from the SDSS DR7 and divide them into four sub-samples according to their velocity dispersions and each sub-sample is further divided into 12 redshift bins. The spectra of the LRGs in each redshift and velocity bin are co-added in order to obtain a combined spectrum with relatively high $S/N$. Adopting the GalexEV/SteLib model, we estimate the mean ages of the LRGs from these combined spectra by the full-spectrum fitting method. We check the reliability of the estimated age by using Monte-Carlo simulations and find that the estimates are robust and reliable. Assuming that the LRGs in each sub-sample and each redshift bin were on average formed at the same time, the Hubble parameter at the present time $H_0$ is estimated from the age--redshift relation obtained for each sub-sample, which is compatible with the $H_0$ value measured by other methods. We demonstrate that a systematic bias (up to $\sim 20%$) may be introduced to the $H_0$ estimation because of recent star formation in the LRGs due to the later major mergers at $z\la 0.4$, but this bias may be negligible for those sub-samples with large velocity dispersions. Using the age--redshift relations obtained from the sub-sample with the largest velocity dispersion or the two sub-samples with high velocity dispersions, we find $H_0= 65^{+7}_{-3}\kmsmpc$ or $H_0= 74^{+5}_{-4}\kmsmpc$ by assuming a spatially flat $\Lambda$CDM cosmology. With upcoming surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), even larger samples of quiescent massive LRGs may be obtained, and thus the Hubble parameter can be measured with high accuracy through the age--redshift relation.

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

Hypervelocity stars are believed to be ejected out from the Galactic center through dynamical interactions of (binary) stars with the central massive black hole(s). In this letter, we report 13 metal-poor F-type hypervelocity star candidates selected from 370,000 stars of the data release 7 of the Sloan Digital Sky Survey. With a detailed analysis of the kinematics of these stars, we find that seven of them were likely ejected from the Galactic center (GC) or the Galactic disk, four neither originated from the GC nor the Galactic disk, and the other two were possibly ejected from either the Galactic disk or other regions. Those candidates which unlikely originated from the GC or the Galactic disk, may be explained by other mechanisms, like the tidal disruption of the Milky Way's dwarf galaxies in the Galactic potential, or the gravitational interactions with a massive black hole at the center of M31 or M32.

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $\Lambda$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

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.

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)

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.

The RXTE and Swift observations on the 2009 outburst from a newly discovered transient and black-hole candidate XTE J1652$-$453 are analyzed. The source was observed by RXTE to behave a sequence of spectral states that are typical to the black hole XRBs. During the first 7 observations, the source is diagnosed to stay in a high/soft state, from the spectrum dominated by soft thermal component which contributes an average of $\sim$ 85% to the X-ray flux at 2-20keV, and from the hardness $\sim$ 0.1 showing up in the hardness-intensity diagram (HID). For the last 20 observations the spectral state is classified as {\sl low/hard} according to an average hardness of $\sim$ 0.8 and a balance between the thermal and the non-thermal components: a power-law component takes ${\geq}80%$ of total 2-20keV flux. Located in between is an {\sl intermediate} state that the source might have experienced. The usual relationship between rms and hardness presents in HRD as well. Throughout the outburst no QPOs are found in XTE J1652$-$453.

We investigate the non-Gaussian features of the IGM at redshift $z\sim 5 - 6$ using Ly$\alpha$ transmitted flux of quasar absorption spectra and cosmological hydrodynamic simulation of the concordance $\Lambda$CDM universe. We show that the neutral hydrogen mass density field and Ly$\alpha$ transmitted flux fluctuations possess all the non-Gaussian features predicted by the log-Poisson hierarchy, which depends only on two dimensionless parameters $\beta$ and $\gamma$, describing, respectively, the intermittence and singularity of the random fields. We find that the non-Gaussianity of the Ly$\alpha$ transmitted flux of quasars from $z=4.9$ to $z=6.3$ can be well reconstructed by the hydrodynamical simulation samples. Although the Gunn-Peterson optical depth and its variance underwent a significant evolution in the redshift range of $5 - 6$, the intermittency measured by $\beta$ is almost redshift-independent in this range. More interesting, the intermittency of quasar's absorption spectra on physical scales $0.1-1$ h$^{-1}$Mpc in redshift $5 - 6$ are found to be about the same as that on physical scales $1-10$ h$^{-1}$Mpc at redshifts $2 - 4$. Considering the Jeans length is less than 0.1 h$^{-1}$Mpc at $z\sim 5$, and $1$ h$^{-1}$Mpc at $z\sim 2$, these results imply that the nonlinear evolution in high and low redshifts will lead the cosmic baryon fluid to a state similar to fully developed turbulence. The log-Poisson high order behavior of current high redshift data of quasar's spectrum can be explained by uniform UV background in the redshift range considered. We also studied the log-Poisson non-Gaussianity by considering inhomogeneous background. With several simplified models of inhomogeneous background, we found the effect of the inhomogeneous background on the log-Poisson non-Gaussianity is not larger than 1-sigma.

It has been argued that the long gamma-ray burst (GRB) of GRB 060614 without an associated supernova (SN) has challenged the current classification and fuel model for long GRBs, and thus a tidal disruption model has been proposed to account for such an event. Since it is difficult to detect SNe for long GRBs at high redshift, the absence of an SN association cannot be regarded as the solid criterion for a new classification of long GRBs similar to GRB 060614, called GRB 060614-type bursts. Fortunately, we now know that there is an obvious periodic substructure observed in the prompt light curve of GRB 060614. We thus use such periodic substructure as a potential criterion to categorize some long GRBs into a new class of bursts, which might have been fueled by an intermediate-mass black hole (IMBH) gulping a star, rather than a massive star collapsing to form a black hole. Therefore, the second criterion to recognize for this new class of bursts is whether they fit the tidal disruption model. From a total of 328 Swift GRBs with accurately measured durations and without SN association, we find 25 GRBs satisfying the criteria for GRB 060614-type bursts: seven of them are with known redshifts and 18 with unknown redshifts. These new bursts are ~6% of the total Swift GRBs, which are clustered into two subclasses: Type I and Type II with considerably different viscous parameters of accretion disks formed by tidally disrupting their different progenitor stars. We suggest that the two different kinds of progenitors are solar-type stars and white dwarfs: the progenitors for four Type I bursts with viscous parameter of around 0.1 are solar-type stars, and the progenitors for 21 Type II bursts with viscous parameter of around 0.3 are white dwarfs. The potential applications of this new class of GRBs as cosmic standard candles are discussed briefly

We identified a large sample of radio quasars, including those with complex radio morphology, from the Sloan Digital Sky Survey (SDSS) and the Faint Images of Radio Sky at Twenty-cm (FIRST). Using this sample, we inspect previous radio quasar samples for selection effects resulting from complex radio morphologies and adopting positional coincidence between radio and optical sources alone. We find that 13.0% and 8.1% radio quasars do not show a radio core within 1.2 and 2 arcsecs of their optical position, and thus are missed in such samples. Radio flux is under-estimated by a factor of more than 2 for an additional 8.7% radio quasars. These missing radio extended quasars are more radio loud with a typical radio-to-optical flux ratio namely radio loudness RL >100, and radio power P >10^{25} W/Hz. They account for more than one third of all quasars with RL>100. The color of radio extended quasars tends to be bluer than the radio compact quasars. This suggests that radio extended quasars are more radio powerful sources, e.g., Fanaroff-Riley type 2 (FR-II) sources, rather than the compact ones viewed at larger inclination angles. By comparison with the radio data from the NRAO VLA Sky Survey (NVSS), we find that for sources with total radio flux less than 3 mJy, low surface brightness components tend to be underestimated by FIRST, indicating that lobes in these faint radio sources are still missed.

In the unification scheme, narrow-lined (type 2) active galactic nuclei (AGN) are intrinsically similar to broad-lined (type 1) AGN with the exception that the line of sight to the broad emission line region and accretion disk is blocked by a dusty torus. The fraction of type 1 AGN measures the average covering factor of the torus. In this paper, we explore the dependence of this fraction on nuclear properties for a sample of low redshift (z <0.35) radio strong (P_{1.4GHz} >10^{23}W/Hz) AGN selected by matching the spectroscopic catalog of Sloan Digital Sky Survey and the radio source catalog of Faint Image of Radio Sky at Twenty cm. After correcting for several selection effects, we find that : (1) type 1 fraction $f_1$ keeps at a constant of ~20 per cent in the [O III] 5007 luminosity range of 40.7< log(L_{[O III]}/ erg/s) <43.5 . This result is significantly different from previous studies, and the difference can be explained by extinction correction and different treatment of selection effects. (2) $f_1$ rises with black hole mass from ~20 per cent (M_bh below 10^8 Msun) to ~30 per cent (M_bh above that). This coincides with the decrease of the fraction of highly-inclined disk galaxies with black hole mass, implying a population of Seyfert galaxies seen as type 2 due to galaxy-scale obscuration in disk when the host galaxy type transfer from bulge-dominant to disk-dominant. (3) $f_1$ is independent of the Eddington ratio for its value between 0.01 and 1; (4) $f_1$ ascends from 15 per cent to 30 per cent in the radio power range of 23< log(P_{1.4GHz}/ W/Hz) <24, then remain a constant at ~30 per cent up to 10^{26} W/Hz.

We present a study of the centroid frequencies and phase lags of the quasi-periodic oscillations (QPOs) as functions of photon energy for GRS 1915+105. It is found that the centroid frequencies of the 0.5-10 Hz QPOs and their phase lags are both energy dependent, and there exists an anti-correlation between the QPO frequency and phase lag. These new results challenge the popular QPO models, because none of them can fully explain the observed properties. We suggest that the observed QPO phase lags are partially due to the variation of the QPO frequency with energy, especially for those with frequency higher than 3.5 Hz.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

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.

This paper aims to study the configuration of two components caused by rotational and tidal distortions in the model of a binary system. The potentials of the two distorted components can be approximated to 2nd-degree harmonics. Furthermore, both the accretion luminosity ($\sigma_{i}$) and the irradiative luminosity are included in stellar structure equations. The equilibrium structure of rotationally and tidally distorted star is exactly a triaxial ellipsoids. A formula describing the isobars is presented, and the rotational velocity and the gravitational acceleration at the primary surface simulated. The results show the distortion at the outer layers of the primary increases with temporal variation and system evolution. Besides, it was observed that the luminosity accretion is unstable, and the curve of the energy-generation rate fluctuates after the main sequence in rotation sequences. The luminosity in rotation sequences is slightly weaker than that in non-rotation sequences. As a result, the volume expands slowly. Polar ejection is intensified by the tidal effect. The ejection of an equatorial ring may be favoured by both the opacity effect and the $g_{e}(\theta,\varphi)$-effect in the binary system.

We demonstrate that the feedback from stellar bulges can play an essential role in shaping the halo gas of galaxies with substantial bulge components by conducting 1-D hydrodynamical simulations. The feedback model we consider consists of two distinct phases: 1) an early starburst during the bulge formation and 2) a subsequent long-lasting mass and energy injection from stellar winds of low-mass stars and Type Ia SNe. An energetic outward blastwave is initiated by the starburst and is maintained and enhanced by the long-lasting stellar feedback. For a MW-like galactic bulge, this blastwave sweeps up the halo gas in the proto-galaxy and heats up the surrounding medium to a scale much beyond the virial radius of the halo, thus the accretion of the halo hot gas can be completely stopped. In addition, the long-lasting feedback in the later phase powers a galactic bulge wind that is reverse-shocked at a large radius in the presence of surrounding intergalactic medium and hence maintains a hot gaseous halo. As the mass and energy injection decreases with time, the feedback evolves to a subsonic and quasi-stable outflow, which is enough to prevent halo gas from cooling. The two phases of the feedback thus re-enforce each-other's impact on the gas dynamics. The simulation results demonstrate that the stellar bulge feedback may provide a plausible solution to the long-standing problems in understanding the MW type galaxies, such as the "missing stellar feedback" problem and the "over-cooling" problem. The simulations also show that the properties of the hot gas in the subsonic outflow state depend sensitively on the environment and the formation history of the bulge. This dependence and variance may explain the large dispersion in the X-ray to B-band luminosity ratio of the low $L_X/L_B$ Es.

we study, in this paper, the non-Gaussian features of the mass density field of neutral hydrogen fluid and the Ly-alpha transmitted flux of QSO absorption spectrum from the point-of-view of self-similar log-Poisson hierarchy. It has been shown recently that, in the scale range from the onset of nonlinear evolution to dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the She-Leveque's scaling formula, which is due to the log-Poisson hierarchical cascade. Since the mass density ratio between ionized hydrogen to total hydrogen is not uniform in space, the mass density field of neutral hydrogen component is not given by a similar mapping of total baryon fluid. Nevertheless, we show, with hydrodynamic simulation samples of the concordance $\Lambda$CDM universe, that the mass density field of neutral hydrogen, is also well described by the log-Poisson hierarchy. We then investigate the field of Ly$\alpha$ transmitted flux of QSO absorption spectrum. Due to redshift distortion, Ly$\alpha$ transmitted flux fluctuations are no longer to show all features of the log-Poisson hierarchy. However, some non-Gaussian features predicted by the log-Poisson hierarchy are not affected by the redshift distortion. We test these predictions with the high resolution and high S/N data of quasars Ly$\alpha$ absorption spectra. All results given by real data, including $\beta$-hierarchy, high order moments and scale-scale correlation, are found to be well consistent with the log-Poisson hierarchy. We compare the log-Poisson hierarchy with the popular log-normal model of the Ly$\alpha$ transmitted flux. The later is found to yield too strong non-Gaussianity at high orders, while the log-Poisson hierarchy is in agreement with observed data.

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 combination of a long duration and the absence of any accompanying supernova clearly shows that GRB 060614 can not be grouped into the two conventional classes of gamma-ray bursts, i.e. the long/soft bursts deemed to be collapsars and the short/hard bursts deemed to be merging binary compact stars. A new progenitor model is required for this anomalous gamma-ray burst. We propose that GRB 060614 might be produced through the tidal disruption of a star by an intermediate mass black hole. In this scenario, the long duration and the lack of any associated supernova are naturally expected. The theoretical energy output is also consistent with observations. The observed 9-s periodicity in the $\gamma$-ray light curve of GRB 060614 can also be satisfactorily explained.

About 10TeV gamma-ray emission within 10 pc region from the Galactic Center had been reported by 4 independent groups. Considering that this TeV gamma-ray emission is produced via a hadronic model, and the relativistic protons came from the tidal disruption of stars by massive black holes, we investigate the spectral nature of the injected relativistic protons required by the hadronic model. The calculation was carried on the tidal disruption of the different types of stars and the different propagation mechanisms of protons in the interstellar medium. Compared with the observation data from HESS, we find for the best fitting that the power-law index of the spectrum of the injected protons is about -1.9, when a red giant star is tidally disrupted, and the effective confinement of protons diffusion mechanism is adopted.

During a core-collapse supernova, absorption of anti-nu_e emitted from the proto-neutron star by protons in the hydrogen envelope produces neutrons and positrons. Neutron capture on protons and positron annihilation then produce gamma rays of 2.22 and 0.511 MeV, respectively. We calculate the fluxes of these gamma rays expected from a supernova with an 11 M_sun progenitor. The flux from neutron capture on protons exponentially decays on a timescale of 564 s, which is determined by neutron decay and capture on protons and 3He nuclei. The peak flux is 2.38x10^{-7}/cm^2/s for a supernova at a distance of 1 kpc. In contrast, the gamma-ray flux from positron annihilation follows the time evolution of the anti-nu_e luminosity and lasts for ~10 s. The peak flux in this case is 6.8x10^{-5}/cm^2/s for a supernova at a distance of 1 kpc. Detection of the above gamma-ray fluxes is beyond the capability of current instruments, and perhaps even those planned for the near future. However, if such fluxes can be detected, they not only constitute a new kind of signals that occur during the gap of several hours between the neutrino signals and the optical display of a supernova, but may also provide a useful probe of the conditions in the surface layers of the supernova progenitor.

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)

Due to the relativistic motion of gamma-ray burst remnant and its deceleration in the circumburst medium, the equal arrival time surfaces at any moment are not spherical, but should be distorted ellipsoids. This will leave some imprints in the afterglows. In this article, we study the effect of equal arrival time surfaces numerically under various conditions, i.e., for isotropic fireballs, collimated jets, density jump conditions, and energy injection events. For each condition, direct comparison between the two instances when the effect is and is not included, is presented. For isotropic fireballs and jets viewed on axis, the effect slightly hardens the spectra and postpones the peak time of afterglows, but does not change the shapes of the spectra and light curves significantly. In the cases when a density jump or an energy injection is involved, the effect smears the variability of the afterglows markedly.

The X-ray flares of NGC 5905, RX J1242.6-1119A, and RX J1624.9+7554 observed by Chandra in 2001 and 2002 have been suggested as the candidate tidal disruption events. The distinct features observed from these events may be used to determine the type of a star tidally disrupted by a massive black hole. We investigate these three events, focusing on the differences for the tidal disruption of a giant star and a main sequence, resulted from their different relation between the mass and the radius. We argue that their X-ray flare properties could be modeled by the partial stripping of the outer layers of a solar type star. The tidal disruption of a giant star is excluded completely. This result may be useful for understanding the growth of a supermassive black hole by capturing stars, versus the growth mode through continuous mass accretion.

This proceedings contribution summarizes first results obtained from a systematic study of 2011 NLS1 galaxies, as presented in Zhou et al. (2006, ApJS 166, 128). The sample was compiled by examining the spectral parameters of galaxies and QSOs derived from SDSS DR3 data. We discuss some preliminary results on the statistic properties of the sample, such as the fraction of NLS1, the properties of broad and narrow emission lines, and emission in other wavebands. The black hole mass--velocity dispersion relation for NLS1s was re-examined using the velocity dispersion values estimated from the stellar absorption spectra of the host galaxies. Preliminary result from an X-ray study for a small subset using data obtained by XMM-Newton is briefly discussed.

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter halos in the $\Lambda$CDM paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though halos and disks may assemble simultaneously. We further demonstrate the validity of the simple approximation for $\Lambda$CDM halos with isotropic velocity distributions using three-dimensional N-body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low-angular momentum orbits in the vicinity of the disk and the growth time scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modeling the response of dark matter halos to the growth of a disk.

(abridged) We use a one-dimensional hydrodynamical code to investigate the effects of preheating on gas accretion and cooling in cold dark matter halos. In the absence of radiative cooling, preheating reduces the amount of gas that can be accreted into a halo, and the accreted gas fraction is determined by the ratio of the initial specific entropy of the gas to the virial entropy of the halo. In the presence of radiative cooling, preheating affects the gas fraction that can cool in two different ways. For small halos with masses <10^12Msun, preheating suppresses gas accretion, but most of the accreted gas can cool. For more massive halos, preheating not only reduces the amount of accreted gas, but also reduces the cooling efficiency. For both small and massive halos, gas cooling is delayed by preheating and in an inside-out fashion if the halo gas is assumed to be a single-phase medium. However, cooling can occur over a wider range of redshifts and radii, if we assume a multi-phase medium. As examples, two specific preheating cases are investigated. In the first case, the preheating entropy is assumed to be proportional to the virial entropy of the halo, as expected from AGN feedback. Such preheating effectively suppresses radiative cooling in halos with M>10^13Msun. We suggest that this may be the reason why the stellar mass function of galaxies breaks sharply at the massive end. Such preheating also helps create the hot diffused halos within which the "radio mode" feedback of AGNs can act effectively. In the second case, we assume the intergalactic medium is warm. Here the total amount of gas that can cool in a halo scales with halo mass as ~M^2, as would be required to match the observed stellar- and HI-mass functions in the current CDM model at the small mass end.

Prompt optical emission from the $\gamma$-ray burst of GRB 041219A has been reported by Vestrand et al. There was a fast rise of optical emission simultaneous with the dominant $\gamma$-ray pulse, and a tight correlation with the prompt $\gamma$-ray emission has been displayed. These indicate that the prompt optical emission and $\gamma$-ray emission would naturally have a common origin. We propose that this optical component can be modeled by considering the Comptonization of $\gamma$-ray photons by an electron cloud. As a result of this mechanism, the arrival time of the optical photons is delayed compared with that of the $\gamma$-rays. We restrict that the lag time to be shorter than 10 s, within which the prompt optical emission is considered to vary simultaneously with the prompt $\gamma$-ray emission. Taking the observations of GRB 041219A into account, we derive the number density of the surrounding electron cloud required for Comptonization. The redshift of GRB 041219A is predicted to be $z\lesssim 0.073$ as well.

This is the first paper in a series dedicated to the study of the emission-line and continuum properties of narrow line Seyfert 1 galaxies (NLS1s). We carried out a systematic search for NLS1s from objects assigned as ``QSOs'' or ``galaxies'' in the spectroscopic sample of the Sloan Digital Sky Survey Data Release 3 (SDSS DR3) by a careful modeling of their emission lines and continua. The result is a uniform sample comprising $\sim 2~000$ NLS1. This sample dramatically increases the number of known NLS1 by a factor of $\sim 10$ over previous compilations. This paper presents the parameters of the prominent emission lines and continua, which were measured accurately with typical uncertainties $<10%$. Taking advantage of such an unprecedented large and uniform sample with accurately measured spectral parameters, we carried out various statistical analysis, some of which were only possible for the first time (abbreviated).

Whether radio intermediate quasars possess relativistic jets as radio-loud quasars is an important issue in the understanding of the origin of radio emission in quasars. In this letter, using the two-epoch radio data obtained during Faint Image of Radio Sky at Twenty centimeter sky (FIRST) and NOAO VLA Sky Survey (NVSS), we identify 89 radio variable sources in the Sloan Digital Sky Survey. Among them, more than half are radio intermediate quasars (RL=f2cm/f2500A<250). For all objects with available multiple band radio observations, the radio spectra are either flat or inverted. The brightness temperature inferred from the variability is larger than the synchrotron-self Compton limit for a stationary source in 87 objects, indicating of relativistic beaming. Considering the sample selection and viewing angle effect, we conclude that relativistic jets probably exist in a substantianl fraction of radio intermediate quasars.

Whether radio intermediate quasars possess relativistic jets as radio-loud quasars is an important issue in the understanding of the origin of radio emission in quasars. In this letter, using the two-epoch radio data obtained during Faint Image of Radio Sky at Twenty centimeter sky (FIRST) and NOAO VLA Sky Survey (NVSS), we identify 89 radio variable sources in the Sloan Digital Sky Survey. Among them, more than half are radio intermediate quasars (RL=f2cm/f2500A<250). For all objects with available multiple band radio observations, the radio spectra are either flat or inverted. The brightness temperature inferred from the variability is larger than the synchrotron-self Compton limit for a stationary source in 87 objects, indicating of relativistic beaming. Considering the sample selection and viewing angle effect, we conclude that relativistic jets probably exist in a substantianl fraction of radio intermediate quasars.

N-body simulations predict that CDM halo-assembly occurs in two phases: 1) a fast accretion phase with a rapidly deepening potential well; and 2) a slow accretion phase characterised by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM halos of the NFW form and provides physical insight into the properties of the mass accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropised by fluctuations in the gravitational potential during the fast accretion phase, we show that gravitational collapse in this phase leads to an inner profile rho(r) ~ r^{-1}. Slow accretion onto an established potential well leads to an outer profile with rho(r) ~ r^{-3}. The concentration of a halo is determined by the fraction of mass that is accreted during the fast accretion phase. Using an ensemble of realistic mass accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time, and hence the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N-body simulations. Using a simple analytic model that captures much of the important physics we show that the inner r^{-1} profile of CDM halos is a natural result of hierarchical mass assembly with a initial phase of rapid accretion.

Non-thermal TeV $\gamma$-ray emission within a multiparsec has been observed from the center region of our Galaxy. We argue that these $\gamma$-rays are the result of transient activity of the massive black hole Sgr A$^*$ that resides at the Galactic Center. Several thousand years ago, the black hole may have experienced an active phase by capturing a red giant star and forming an accretion disk, temporarily behaving like an active galactic nucleus. A powerful jet, which contains plenty of high speed protons, was launched during the process. These runaway protons interact with the dense ambient medium, producing TeV $\gamma$-ray emission through the $\pi^0$-decay process. We show that the total energy deposited in this way is large enough to account for observations. The diffusion length of protons is also consistent with the observed size of the TeV source.

Using the radio observations by FIRST and NVSS, we build a sample of 151 radio variable quasars selected from the Sloan Digital Sky Survey Data Release 3 (SDSS DR3). Six (probably another two) among them are classified as broad absorption line (BAL) quasars, with radio flux variations of a few 10 percent within 1.5-5 years. Such large amplitudes of the variations imply brightness temperatures much higher than the inverse Compton limits (10$^{12}$ K) in all the BAL quasars, suggesting the presence of relativistic jets beaming toward the observer. The angle between the outflow and the jet is constrained to be less than $\sim 20^{\circ}$. Such BAL quasars with polar outflows are beyond the simple unification models of BAL quasars and non-BAL quasars, which hypothesize that BAL quasars are normal quasars seen nearly edge-on.

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.

Double-peaked Balmer lines have been observed in about 150 AGNs and were interpreted preferably as emission from relativistic accretion disks. In this paper, we report the discovery of extreme double-peaked lines in SDSS J0942+0900 and SDSS J1417+6141. The FWHM of the Halpha line ~40,600 km/s in the first object is almost twice as large as the broadest one previously known. By comparing the line profile with accretion disk models, we find that most of the line flux is emitted from a narrow annulus around 100Rg in SDSS J0942+0900, and from a disk of radii between 100 and 400Rg in SDSS J1417+6141. This is the first time that an accretion disk at radii below 100Rg can be directly probed through optical emission lines. A certain asymmetry in the disk is required for both objects. Another much weaker broad Halpha component (EW~20A, and FWHM 4000km/s) is also detected in both objects. Both objects show weak radio emission and strong high-ionization narrow lines.

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)

The observed optical luminosity in 5100 angstorm and black hole mass correlation is used to probe the magnetic field of black holes harbored in active galactic nuclei(AGNs). The model is based on the assumption that the disk is heated by energy injection due to the magnetic coupling(MC) process and the gravitational dissipation due to accretion. The MC process can transfer energy and angular momentum from a rotating Kerr black hole to its surrounding disk. The relation of optical luminosity in 5100 angstorm and black hole mass as functions of the spin and magnetic field of the black hole is modelled. The model predicts that optical luminosity in 5100 angstorm emitted from the disk is sensitive to the strength of the poloidal component of the magnetic field on the BH horizon. Based on the observations of optical luminosity in 5100 angstorm for 143 AGN sources, we obtain the correlation between mass and magnetic field of black hole. And we compared out result with the approximate result between mass and magnetic field of black hole derived from the condition in the standard accretion disc theory.

We show that changes in the orientation of the inner accretion disk of Cygnus X-1 affect the shape of the broad Fe K$\alpha$ emission line emitted from this object, in such a way that eV-level spectral resolution observations (such as those that will be carried out by the {\it ASTRO-E2} satellite) can be used to analyze the dynamics of the disk. We here present a potential diagnostic tool, supported by numerical simulations, by which a few observations of Cygnus X-1, separated in time, can determine whether its accretion disk actually precesses, and if so, determine its period and precession angle. This approach could also be used for similar studies in other microquasar systems.

A maser curvature emission mechanism in the presence of curvature drift is used to explain the optical pulsations from anomalous X-ray pulsars. For the source of AXP0142+61,the optical pulsation occurs at the radial distance $R(\nu_M)\sim 4.75\times 10^9$ cm to the neutron star. The corresponding curvature maser frequency is about $\nu_M\approx1.39\times 10^{14}$ Hz. The result is consistent with the observation of the optical pulsations from the anomalous X-ray pulsar 4U0142+61.

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 propose the study of long-term variations in the Fe Kalpha line profile as a discriminator for binary black hole (BH) systems. The existence of a binary BH in the center of an active galaxy will produce a particular signature in the evolution of the line profile, as a result of disk precession. This signature is a periodic variation of the position of the blue edge of the profile, accompanied by periodic change of its intensity. We show that detection of the former is clearly within the observational capabilities of planned X-ray missions. Detecting a periodic variation of line profiles would provide the first direct evidence for precessing discs in active galactic nuclei, as opposed to the existing evidence supporting only the precession of jets. We apply these ideas to 3C 273.

In this paper, the spectra of the persistent emission from anomalous X-ray pulsars (AXPs) and their variation with spin-down rate $\dot{\Omega}$ is considered. Firstly, based on an accretion-powered model, the influences of both magnetic field and mass accretion rate on the spectra properties of AXPs are addressed. Subsequently, the relation between the spectral property of AXPs and mass accretion rate $\dot{M}$ is investigated. The result shows that there exists a linear correlation between the photon index and mass accretion rate, and the spectral hardness increases with increasing $\dot{M}$. A possible emission mechanism for the explanation of spectral properties of AXPs is also discussed.

A possible accretion model associated with the ionization instability of quasar disks is proposed to address the growth of the central black hole harbored in the host galaxy.The mass ratio between black hole and its host galactic bulge is a nature consequence of our model.

This paper continues the study of the properties of an accretion disc rotating around a non-baryonic (assumed super-massive) compact object. This kind of objects, generically known as boson stars, were earlier proposed as a possible alternative scenario to the existence of super-masive black holes in the center of every galaxy. A dilute boson star has also been proposed as a large part of the non-baryonic dark matter, flattening galactic rotational velocities curves. In this contribution, we compute the profile of the emission lines of Iron; its shape has been for long known as a useful diagnosis of the space-time geometry. We compare with the case of a Schwarzschild black hole, concluding that the differences are observationally distinguishable.

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)

Controversy regarding whether Narrow Line Seyfert 1 galaxies (NLS1s) follow the same black hole --- galactic bulge relation as normal galaxies has been raised by Mathur et al. (2001) and Ferrarese et al. (2001). In this paper, the correlation between the black hole mass and the velocity dispersion of the narrow line region, indicated by the [OIII] line width for a heterogeneous sample of NLS1s, is examined. We show that the non-virial component subtracted [OIII] width, which may well represent the stellar velocity dispersion ($\sigma$), correlates with the estimated black hole mass $M_bh$, made using the empirical law for broad line region sizefrom the reverberation mapping measurements. Considering measurement errors, we find that the relation between $M_bh$ and the [OIII] width in active galactic nuclei (both NLS1s in this paper and normal Seyferts in Nelson 2000) is consistent with that defined in nearby hot galaxies (ellipticals and spiral bulge) but with more scatter. It seems that there is no clear difference in the relation between $M_bh$ and $\sigma$ (represented by the [OIII] width for AGN) in both NLS1s and normal broad line active galactic nuclei from the same relation defined by nearby hot galaxies.

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.

A possible relation between the high dispersion in metallicity of metal-poor halo stars and the minor merger processes in the history of the Galaxy is presented. The foreign populations of stars in the satellites through minor merger processes, make considerable contributions to the Galactic halo, therefore alter the intrinsic distribution of metallicities of the halo stars. Our model reproduces the observed trends in the metallicity distribution of the extremely metal-poor halo field stars in the Galaxy. Taken all the parameters for our model, a density of $1.9\times 10^{-3}kpc^{-3}$ of extremely metal-poor halo stars in $-5<[Fe/H]<-4$ is predicted, which explains why no such stars have been observed so far.

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 compile the multi-band spectral energy distribution (SED) of Seyfert 1 galaxy NGC4593, and investigate its optical-to-ultraviolet (OUV) continuum and iron K\alpha line profile. No strong ``big blue bump'' is found in OUV band. The iron K\alpha line profile is not as broad as expected from a disk which extends down to the last stable orbit. Both of them can be modeled by a truncated thin disk, and support the existence of a truncated thin disk in the system. It is of much interest that NGC4593 may be a typical object with an accretion rate \dot{m}=\dot{M}/M_{Edd} \sim 0.01-0.1 and harbors a two component accretion flow: an outer thin disk and an inner hot flow--possibly ADAF, analogous to the low and intermediate state of low mass X-ray binaries (LMXB).

Project GRAND is an extensive air shower array utilizing position sensitive detectors of proportional wire chambers. The 64 detectors deployed in a field 100 m x 100 m are located at 86.2 deg W and 41.7 deg N, at 220 m above sea level. The project was completed about two years ago and has been taking data, simultaneously on two triggers: 1) multiple-hit coincidence triggers which collect data on extensive air showers, and 2) 200 Hz triggers which collect the single tracks stored in each station during the last 5 msec.

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.

We suggest that an extreme Kerr black hole with a mass $\sim 10^6M_\odot$, a dimensionless angular momentum $A\sim 1$ and a marginal stable orbital radius $r_{ms}\sim 3r_s\sim 10^{12}M_6 cm$ located in a normal galaxy, may produced a Gamma-ray Burst by capturing and disrupting a star. During this period, a transient accretion disk is formed and a strong transient magnetic field $\sim 2.4\times 10^9M_6^{-1/2}$ Gauss, lasting for $r_{ms}/c\sim 30 M_6 s$, may be produced in the inner boundary of the accretion disk. A large amount of rotational energy of the black hole is extracted and released in the ultra relativistic jet with a bulk Lorentz factor $\Gamma$ larger than $10^3$ via Blandford-Znajek process. The relativistic jet energy can be converted into $\gamma$-ray radiation via internal shock mechanism. The gamma-ray burst (GRB) duration should be the same as that of the life time of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be $r_{ms}/h\sim (10 - 10^2)$ because the disk material is likely broken into pieces with the size about the thickness of the disk $h$ at the cusp ($2r_s\le r \le 3r_s$). The shortest rising time of the burst estimated from this model is $\sim h/\Gamma c\sim 3\times 10^{-4}\Gamma^{-1}_3(h/r)_{-2}M_6$ s. The model gamma-ray burst density rate is also estimated.

Anisotropic emission of the ionizing continuum is a general prediction of the accretion disk models. In this paper, we present the results of correlation analysis of the UV emission line and UV to X-ray continuum properties for a large sample of broad emission line AGNs observed with ROSAT, IUE and HST. We find strong correlations between the CIV/$Ly\alpha$ ratio, the equivalent width of CIV, and the UV to soft X-ray spectral slope. The results are in good agreement with the photoionization calculation, suggesting that the overall ionizing continuum can well match the observed UV to soft X-ray spectrum. These results are consistent with the assumption of isotropic ionizing continuum shape. Our analysis suggests a small range for the ``big blue bump'' cutoff energy for the objects in this sample, consistent with the similar results of Laor et al. (1997) and Walter & Fink (1993) based on the continuum properties . The mean UV-to-X-ray spectral slope is similar to the soft X-ray spectral slope. This similarity also holds for radio-loud and radio-quiet objects separately. This suggests that the two might be drawn from the same distribution. The two spectral slopes are only weakly correlated. The UV to X-ray spectral index is correlated with absolute optical magnitude. This result confirms the earlier suggestion that the ionizing continua are softer for higher luminosity objects.