Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We investigated the pulsar radio luminosity ($L$), emission efficiency (ratio of radio luminosity to its spin-down power $\dot{E}$), and death line in the diagram of magnetic field (B) versus spin period (P), and found that the dependence of pulsar radio luminosity on its spin-down power ($L-\dot{E}$) is very weak, shown as $L\sim\dot{E}^{0.06}$, which deduces an equivalent inverse correlation between emission efficiency and spin-down power as $\xi\sim \dot{E}^{-0.94}$. Furthermore, we examined the distributions of radio luminosity of millisecond and normal pulsars, and found that, for the similar spin-down powers, the radio luminosity of millisecond pulsars is about one order of magnitude lower than that of the normal pulsars. The analysis of pulsar radio flux suggests that this correlations are not due to a selective effect, but are intrinsic to the pulsar radio emission physics. Their radio radiations may be dominated by the different radiation mechanisms. The cut-off phenomenon of currently observed radio pulsars in B-P diagram is usually referred as the "pulsar death line", which corresponds to $\dot{E}\approx 10^{30}$ erg/s and is obtained by the cut-off voltage of electron acceleration gap in the polar cap model of pulsar proposed by Ruderman and Sutherland. Observationally, this death line can be inferred by the actual observed pulsar flux $S\approx $1mJy and 1kpc distance, together with the maximum radio emission efficiency of 1\%. At present, the actual observed pulsar flux can reach 0.01mJy by FAST telescope, which will arise the observational limit of spin-down power of pulsar as low as $\dot{E}\approx 10^28$ erg/s. This means that the new death line is downward shifted two orders of magnitude, which might be favorably referred as the "observational limit-line", and accordingly the pulsar theoretical model for the cut-off voltage of gap should be heavily modified.

The transonic phenomenon of black hole accretion and the existence of the photon sphere are the characteristics of strong gravitational fields near a black hole horizon. In this work, we study spherical accretion flow onto a general parametrized spherically symmetric black hole spacetimes. For this purpose, we analyze the accretion process of various perfect fluids, such as the isothermal fluid of ultra-stiff, ultra-relativistic, and sub-relativistic types and polytropic fluid, respectively. The influences of extra parameters beyond the Schwarzschild black hole in the general parameterized spherically symmetric black hole on the flow behaviors of the above-mentioned test fluids are studied in detail. In addition, by studying the accretion of ideal photon gas, we further discuss the correspondence between the sonic radius of accreting photon gas and the photon sphere for the general parameterized spherically symmetric black hole. Some possible future extensions of our analysis are also discussed.

The Hubble parameter $H(z)$ is directly related to the expansion of our Universe. It can be used to study dark energy and constrain cosmology models. In this paper, we propose that $H(z)$ can be measured using fast radio bursts (FRBs) with redshift measurements. We use dispersion measures contributed by the intergalactic medium, which is related to $H(z)$, to measure Hubble parameter. We find that 500 mocked FRBs with dispersion measures and redshift information can accurately measure Hubble parameters using Monte Carlo simulation. The maximum deviation of $H(z)$ from standard $\Lambda$CDM model is about 6\% at redshift $z= 2.4$. We also test our method using Monte Carlo simulation. Kolmogorov-Smirnov (K-S) test is used to check the simulation. The $p$-value of K-S test is 0.23, which confirms internal consistency of the simulation. In future, more localizations of FRBs make it as an attractive cosmological probe.

Recently a novel four-dimensional Einstein-Gauss-Bonnet (4EGB) theory of gravity was proposed by Glavan and Lin [D. Glavan and C. Lin, Phys. Rev. Lett. 124, 081301 (2020)] which bypasses the Lovelock' s theorem with the re-scaling of the Gauss-Bonnet coupling constant $\alpha \to \alpha/(D-4)$ in the limit $D\to 4$. Here we study the physical properties of the electromagnetic radiation emitted from a thin accretion disk around the static spherically symmetric black hole in the 4EGB gravity. For this purpose, we assume the disk is in a steady-state and in hydrodynamic and thermodynamic equilibrium so that the emitted electromagnetic radiation is a black body spectrum. We study in detail the effects of the Gauss-Bonnet coupling constant $\alpha$ in 4EGB gravity on the energy flux, temperature distribution, and electromagnetic spectrum of the disk. It is shown that with the increases of the parameter $\alpha$, the energy flux, temperature distribution, and electromagnetic spectrum of the accretion disk all increase. Besides, we also show that the accretion efficiency increases as the growth of the parameter $\alpha$. Our results indicate that the thin accretion disk around the static spherically symmetric black hole in the 4EGB gravity is hotter, more luminosity, and more efficient than that around a Schwarzschild black hole with the same mass for a positive $\alpha$, while it is cooler, less luminosity, and less efficient for a negative $\alpha$.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

In this paper, we study the motion of test particles around two exact charged black-hole solutions in Einstein-{\AE}ther theory. Specifically, we first consider the quasi-periodic oscillations (QPOs) and their resonances generated by the particle moving in the Einstein-{\AE}ther black hole and then turn to study the periodic orbits of the massive particles. For QPOs, we drop the usually adopted assumptions $\nu_U=\nu_\theta$, $\nu_L=\nu_r$, and $\nu_U/\nu_L=3/2$ with $\nu_U$ ($\nu_L$) and $\nu_r$ ($\nu_\theta$) being the upper (lower) frequency of QPOs and radial (vertical) epicyclic frequency of the orbiting particles, respectively. Instead, we put-forward a new working ansatz for which the Keplerian radius is much closer to that of the innermost stable circular orbit and explore in details the effects of the {\ae}ther field on the frequencies of QPOs. We then realize good curves for the frequencies of QPOs, which fit to data of three microquasars very well by ignoring any effects of rotation and magnetic fields. For the study of the periodic orbits of massive particle moving in the two Einstein-{\AE}ther black hole solutions, we first consider the geodesic equations in the two black holes and solve the geodesic equation for one of black hole (the black hole with \ae{}ther parameter $c_{14}=0$ but $c_{123}\neq 0$) analytically. The innermost stable circular orbits are also analyzed and we find the isco radius increases with increasing $c_{13}$ for the first type black hole while decreases with increasing $c_{14}$ for the second one. We also obtain several periodic orbits and find that they share similar taxonomy schemes as the periodic equatorial orbits in the Schwarzschild/Kerr metrics. The equations for null geodesics are also briefly considered, where we study circular photon orbits and bending angles for gravitational lensing.

In this paper, we construct an effective rotating loop quantum black hole (LQBH) solution, starting from the spherical symmetric LQBH by applying the Newman-Janis algorithm modified by Azreg-A\"{i}nou's non-complexification procedure, and study the effects of loop quantum gravity (LQG) on its shadow. Given the rotating LQBH, we discuss its horizon, ergosurface, and regularity as $r \to 0$. Depending on the values of the specific angular momentum $a$ and the polymeric function $P$ arising from LQG, we find that the rotating solution we obtained can represent a regular black hole, a regular extreme black hole, or a regular spacetime without horizon (a non-black-hole solution). We also study the effects of LQG and rotation and show that, in addition to the specific angular momentum, the polymeric function also causes deformations in the size and shape of the black hole shadow. Interestingly, for a given value of $a$ and inclination angle $\theta_0$, the apparent size of the shadow monotonically decreases, and the shadow gets more distorted with increasing $P$. We also consider the effects of $P$ on the deviations from the circularity of the shadow and find that the deviation from circularity increases with increasing $P$ for fixed values of $a$ and $\theta_0$. Additionally, we explore the observational implications of $P$ in comparing with the latest Event Horizon Telescope (EHT) observation of the supermassive black hole, M87*. The connection between the shadow radius and quasinormal modes in the eikonal limit as well as the deflection of massive particles are also considered.

The properties of hydrogen at high pressure have wide implications in astrophysics and high-pressure physics. Its phase change in the liquid is variously described as a metallization, H2-dissociation, density discontinuity or plasma phase transition. It has been tacitly assumed that these phenomena coincide at a first-order liquid-liquid transition (LLT). In this work, the relevant pressure-temperature conditions are thoroughly explored with first-principles molecular dynamics. We show there is a large dependency on exchange-correlation functional and significant finite size effects. We use hysteresis in a number of measurable quantities to demonstrate a first-order transition up to a critical point, above which molecular and atomic liquids are indistinguishable. At higher temperature beyond the critical point, H2-dissociation becomes a smooth cross-over in the supercritical region that can be modelled by a pseudo-transition, where the H2-2H transformation is localized and does not cause a density discontinuity at metallization. Thermodynamic anomalies and counter-intuitive transport behavior of protons are also discovered even far beyond the critical point, making this dissociative transition highly relevant to the interior dynamics of Jovian planets. Below the critical point, simulation also reveals a dynamic H2-2H chemical equilibrium with rapid interconversion, showing that H2 and H are miscible. The predicted critical temperature lies well below the ionization temperature. Our calculations unequivocally demonstrate that there are three distinct regimes in the liquid-liquid transition of warm dense hydrogen.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

N-body simulations are essential tools in physical cosmology to understand the large-scale structure (LSS) formation of the Universe. Large-scale simulations with high resolution are important for exploring the substructure of universe and for determining fundamental physical parameters like neutrino mass. However, traditional particle-mesh (PM) based algorithms use considerable amounts of memory, which limits the scalability of simulations. Therefore, we designed a two-level PM algorithm CUBE towards optimal performance in memory consumption reduction. By using the fixed-point compression technique, CUBE reduces the memory consumption per N-body particle toward 6 bytes, an order of magnitude lower than the traditional PM-based algorithms. We scaled CUBE to 512 nodes (20,480 cores) on an Intel Cascade Lake based supercomputer with $\simeq$95\% weak-scaling efficiency. This scaling test was performed in "Cosmo-$\pi$" -- a cosmological LSS simulation using $\simeq$4.4 trillion particles, tracing the evolution of the universe over $\simeq$13.7 billion years. To our best knowledge, Cosmo-$\pi$ is the largest completed cosmological N-body simulation. We believe CUBE has a huge potential to scale on exascale supercomputers for larger simulations.

The anti- and positive correlations between X-ray photon index and Eddington-scaled X-ray luminosity were found in decay phase of X-ray binary outbursts and a sample of active galactic nuclei in former works. We further systematically investigate the evolution of X-ray spectral index, along the X-ray flux and Eddington ratio in eight outbursts of four black-hole X-ray binaries, where all selected outbursts have observational data from Rossi X-ray Timing Explorer in both rise and decay phases. In the initial rise phase, the X-ray spectral index is anti-correlated with the flux, and the X-ray spectrum quickly softens when the X-ray flux is approaching the peak value. In the decay phase, the X-ray photon index and the flux follow two different positive correlations, and they become anti-correlated again when the X-ray flux is below a critical value, where the anti-correlation part follow the same trend as that found in the initial rise phase. Compared with other X-ray binaries, GRO J1655-40 has an evident lower critical Eddington ratio for the anti- and positive transition, which suggests that its black-hole mass and distance are not well constrained or its intrinsic physic is different.

The discovery of changing-look active galactic nuclei (CL AGNs), with appearance and disappearance of broad emission lines and/or with strong variation of line-of-sight column density within a few years, challenges the AGN unification model. We explore the physical mechanisms based on the X-ray spectral evolution for a sample of 15 CL AGNs. We find that the X-ray photon index, $\Gamma$, and Eddington-scaled X-ray luminosity, $L_{\rm 2-10 keV}/L_{\rm Edd}$, follow negative and positive correlations when $L_{\rm 2-10 keV}/L_{\rm Edd}$ is lower and higher than a critical value of $\sim 10^{-3}$. This different X-ray spectral evolution is roughly consistent with the prediction of the accretion-mode transition (e.g., clumpy cold gas or cold disk to advection dominated accretion flow, or vice visa). With quasi-simultaneous X-ray and optical spectrum observations within one year, we find that the CL AGNs observed with and without broad emission lines stay in the positive and negative part of the $\Gamma-L_{\rm 2-10 keV}/L_{\rm Edd}$ correlation respectively. Our result suggest that the change of the accretion mode may be the physical reason for the CL AGNs.

The effective field theory (EFT) of inflation provides an essential picture to explore the effects of the unknown high energy physics in the single scalar field inflation models. For a generic EFT of inflation, possible high energy corrections to simple slow-roll inflation can modify both the propagating speed and dispersion relations of the cosmological scalar and tensor perturbations. With the arrival of the era of precision cosmology, it is expected that these high energy corrections become more important and have to be taken into account in the analysis with future precise observational data. In this paper we study the observational predictions of the EFT of inflation by using the third-order uniform asymptotic approximation method. We calculate explicitly the primordial power spectra, spectral indices, running of the spectral indices for both scalar and tensor perturbations, and the ratio between tensor and scalar spectra. These expressions are all written in terms of the Hubble flow parameters and the flow of four new slow-roll parameters and expanded up to the next-to-leading order in the slow-roll expansions so they represent the most accurate results obtained so far in the literature. The flow of the four new slow-roll parameters, which arise from the four new operators introduced in the action of the EFT of inflation, can affect the primordial perturbation spectra at the leading-order and the corresponding spectral indices at the next-to-leading order.

The recently published fourth Fermi Large Area Telescope source catalog (4FGL) reports 5065 gamma-ray sources in terms of direct observational gamma-ray properties. Among the sources, the largest population is the Active Galactic Nuclei (AGN), which consists of 3137 blazars, 42 radio galaxies, and 28 other AGNs. The blazar sample comprises 694 flat-spectrum radio quasars (FSRQs), 1131 BL Lac-type objects (BL Lacs), and 1312 blazar candidates of an unknown type (BCUs). The classification of blazars is difficult using optical spectroscopy given the limited knowledge with respect to their intrinsic properties, and the limited availability of astronomical observations. To overcome these challenges, machine learning algorithms are being investigated as alternative approaches. Using the 4FGL catalog, a sample of 3137 Fermi blazars with 23 parameters is systematically selected. Three established supervised machine learning algorithms (random forests (RFs), support vector machines (SVMs), artificial neural networks (ANNs)) are employed to general predictive models to classify the BCUs. We analyze the results for all of the different combinations of parameters. Interestingly, a previously reported trend the use of more parameters leading to higher accuracy is not found. Considering the least number of parameters used, combinations of eight, 12 or 10 parameters in the SVM, ANN, or RF generated models achieve the highest accuracy (Accuracy $\simeq$ 91.8\%, or $\simeq$ 92.9\%). Using the combined classification results from the optimal combinations of parameters, 724 BL Lac type candidates and 332 FSRQ type candidates are predicted; however, 256 remain without a clear prediction.

The effective field theory (EFT) of inflation provides a natural framework to study the new physical effects on primordial perturbations. Recently a healthy extension of the EFT of inflation with high-order operators has been proposed, which avoids ghosts and meanwhile leads to a nonlinear dispersion relation of the scalar perturbations. This paper is devoted to studying the effects of these high-order operators by using the uniform asymptotic approximation method. In particular, we first construct the approximate analytical solution to the mode function of the scalar perturbations. Because of the presence of the high-order operators, the perturbation modes usually experience a period of non-adiabatic evolution before they cross the Hubble radius, which could lead to the production of excited states and modifications of the primordial perturbation spectrum. However, we show that the modified power spectrum is still nearly scale-invariant and the presence of the high-order operators can only affect the overall amplitude of the spectrum. In particular, after showing explicitly the impact of these new effects on particle production rate and perturbation spectrum, we explore their origin in detail.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We construct samples of jetted active galactic nuclei (AGNs) with low-frequency radio data from the recent released TGSS AD1 catalog at 150 MHz. With these samples, we compare the properties of jet power for blazars, radio-loud narrow-line Seyfert 1 galaxies (RL-NLS1s), young radio sources (YRSs) and radio galaxies. The jet-disk connection, and the unification of jetted AGNs are explored. On the Eddington ratio --- Eddington-scaled jet power plane, jetted AGNs can be generally divided into two populations. Low power radio galaxies, low excitation FR IIs (LEG/FR IIs), and most YRSs show larger jet power than accretion power, while FSRQs, RL-NLS1s, and high excitation FR IIs (HEG/FR IIs) are on the contrary. LEG/FR IIs share similar jet power properties with HEG/FR IIs, while their accretion properties are different with the latter. These facts suggest an evolutional sequence from HEG/FR IIs, LEG/FR IIs to FR Is, where the accretion and jet activities get dimmed gradually. LEG/FR IIs are the transitional objects that accretion processes have switched off, while jets are still active. After correcting the contribution from radio core of blazars, the unification between blazars and radio galaxies is confirmed with the jet power distributions. The unification involved RL-NLS1s is more complicated than the simple scenarios of black hole growth or orientation effect. In addition, our results manifest that low synchrotron peaked BL Lacs (LBLs) contain two distinct groups on the distribution of jet power, with one group similar with FSRQs and the other similar with intermediate synchrotron peaked BL Lacs. The LBLs with higher jet powers might be the aligned counterparts of LEG/FR IIs.

We formulate and solve the problem of spherically symmetric, steady state, adiabatic accretion onto a Schwarzschild-like black hole obtained recently. We derive the general analytic expressions for the critical points, the critical velocity, the critical speed of sound, and subsequently the mass accretion rate. The case for polytropic gas is discussed in detail. We find the parameter characterizing the breaking of Lorentz symmetry will slow down the mass accretion rate, while has no effect on the gas compression and the temperature profile below the critical radius and at the event horizon.

Motivated by the next generation of gravitational wave (GW) detectors, we study the wave mechanics of a twisted light beam in the GW perturbed spacetime. We found a new gravitational dipole interaction of photons and gravitational waves. Physically, this interaction is due to coupling between the angular momentum of twisted light and the GW polarizations. We demonstrate that for the higher-order Laguerre-Gauss (LG) modes, this coupling effect makes photons undergoing dipole transitions between different orbital-angular-momentum(OAM) eigenstates, and leads to some measurable optical features in the 2-D intensity pattern. It offers an alternative way to realize precision measurements of the gravitational waves, and enables us to extract more information about the physical properties of gravitational waves than the current interferometry. With a well-designed optical setup, this dipole interaction is expected to be justified in laboratories.

The gas pixel detector (GPD) is designed and developed for high-sensitivity astronomical X-ray polarimetry, which is a new window about to open in a few years. Due to the small mass, low power, and compact geometry of the GPD, we propose a CubeSat mission Polarimeter Light (PolarLight) to demonstrate and test the technology directly in space. There is no optics but a collimator to constrain the field of view to 2.3 degrees. Filled with pure dimethyl ether (DME) at 0.8 atm and sealed by a beryllium window of 100 micron thick, with a sensitive area of about 1.4 mm by 1.4 mm, PolarLight allows us to observe the brightest X-ray sources on the sky, with a count rate of, e.g., ~0.2 counts/s from the Crab nebula. The PolarLight is 1U in size and mounted in a 6U CubeSat, which was launched into a low Earth Sun-synchronous orbit on October 29, 2018, and is currently under test. More launches with improved designs are planned in 2019. These tests will help increase the technology readiness for future missions such as the enhanced X-ray Timing and Polarimetry (eXTP), better understand the orbital background, and may help constrain the physics with observations of the brightest objects.

In the third catalog of active galactic nuclei detected by the Fermi-LAT (3LAC) Clean Sample, there are 402 blazars candidates of uncertain type (BCU). Due to the limitations of astronomical observation or intrinsic properties, it is difficult to classify blazars using optical spectroscopy. The potential classification of BCUs using machine learning algorithms is essential. Based on the 3LAC Clean Sample, we collect 1420 Fermi blazars with 8 parameters of {\gamma}-ray photon spectral index, radio flux, flux density, curve significance, the integral photon flux in 100 to 300 MeV, 0.3 to 1 GeV, 10 to 100 GeV and variability index. Here, we apply 4 different supervised machine learning (SML) algorithms (\emph{Decision trees, Random forests, support vector machines and Mclust Gaussian finite mixture models}) to evaluate the classification of BCUs based on the direct observational properties. All the 4 methods can perform exceedingly well with a more accuracy and can effective forecast the classification of Fermi BCUs. The evaluating results show the results of these methods (SML) are valid and robust, where, about 1/4 sources are FSRQs and 3/4 are BL Lacs in 400 BCUs, which are consistent with some other recent results. Although a number of factors influence the accuracy of SML, the results are stable at a fixed ratio 1:3 between FSRQs and BL Lacs, which suggests that the SML can provides an effective method to evaluate the potential classification of BCUs. Among the 4 methods, Mclust Gaussian Mixture Modelling has the highest accuracy for our training sample (4/5, seed=123).

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

While loop quantum cosmology (LQC) predicts a robust quantum bounce of the background evolution of a Friedmann-Robertson-Walker (FRW) spacetime prior to the standard slow-roll inflation, whereby the big bang singularity is resolved, there are several different quantization procedures to cosmological perturbations, for instance, {\em the deformed algebra, dressed metric, and hybrid quantizations}. This paper devotes to study the quantum bounce effects of primordial perturbations in the hybrid approach. The main discrepancy of this approach is the effective positive mass at the quantum bounce for the evolution of the background that is dominated by the kinetic energy of the inflaton field at the bounce, while this mass is always nonpositive in the dressed metric approach. It is this positivity of the effective mass that violates the adiabatic evolution of primordial perturbations at the initial moments of the quantum bounce. With the assumption that the evolution of the background is dominated by the kinetic energy of the inflaton at the bounce, we find that the effective potentials for both scalar and tensor perturbations can be well approximately described by a P\"{o}schl-Teller (PT) potential, which allows us to find analytical solutions of perturbations, and from these analytical expressions we are able to study the non-adiabatic evolution of primordial perturbations in details. In particular, we derive their quantum bounce effects and investigate their observational constraints. In addition, the impacts of quantum bounce effects on the non-Gaussinity and their implication on the explanations of observed power asymmetry in CMB have also been explored.

We construct a large sample of gamma-ray blazars with low-frequency radio data using the recent released TGSS AD1 catalog at 150 MHz. The radio/gamma connections of blazars are compared from 143 GHz to 150 MHz. The radio flux density at all radio frequencies shows strong correlation with gamma-ray flux for blazars, as well as for the two subclasses, FSRQs and BL Lacs. But the correlations get worse from high to low radio frequencies, which indicates that the low-frequency radio emission is the mixture of extended and core components for blazars. In addition, we find that the correlation between 150 MHz radio flux density and gamma-ray flux is more significant for BL Lacs than that for FSRQs. The slope for the luminosity correlation between radio and gamma-ray also get flatter than unity at 150 MHz. These results indicate that the core dominance at 150 MHz for BL Lacs is larger than that for FSRQs. We also compare the radio luminosity from direct TGSS observation and the extended radiation at 150 MHz for blazars. The results show that the ratio between core and extended component at 150 MHz is about 1:1 on average.

By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW170104, GW170814, and GW150914, while with the DM accretion, it can cover all observations well. For the higher metallicity ($Z = 0.02$), our model can not explain the mergers of high-mass binary BHs like GW170104, GW170814, and GW150914. We estimate that the merger rate of binary BHs lies between $\rm55Gpc^{-3}yr^{-1}$ to $\rm197Gpc^{-3}yr^{-1}$.

The composition and energy dissipation in jets are two of the fundamental questions of jet physics that are not fully understood. In this paper, we attempt to constrain the composition, magnetization as well as radiative efficiency for blazar with the recently released low-frequency radio catalog of the TIFR GMRT Sky Survey at 150 MHz. The jet power estimated from the low-frequency radio emission is much lower than that derived from spectral energy distribution fitting assuming one proton per electron. Assuming the jet power estimated from low-frequency radio emission is physical, the fraction of electron/positron pairs can be constrained with $n_{\rm pairs}/n_{\rm p} \sim 10$. By comparing the power carried by magnetic field and radiation with the jet power estimated from the low-frequency radio emission, we find both relatively high magnetization parameter of $\sigma \sim 0.5$ and radiative efficiency of $\eta \sim 0.4$ in the dissipation region of blazars. These results suggest that the magnetic reconnection processes may play an important role in the energy dissipation of blazars. We also explore the connection between these three parameters ($n_{\rm pairs}/n_{\rm p}$, $\sigma$, and $\eta$) and the black hole mass, disk luminosity as well as Eddington ratio. No significant correlation is found, except that $\sigma$ shows possible correlation with disk luminosity.

The production site of gamma-rays in blazar jet is an unresolved problem. We present a method to locate gamma-ray emission region in the framework of one-zone emission model. From measurements of core-shift effect, the relation between the magnetic field strengths ($B'$) in the radio cores of jet and the distances ($R$) of these radio cores from central supermassive black hole (SMBH) can be inferred. Therefore once the magnetic field strength in gamma-ray emission region ($B'_{\rm diss}$) is obtained, one can use the relation of $B'$-$R$ to derive the distance ($R_{\rm diss}$) of gamma-ray emission region from SMBH. Here we evaluate the lower limit of $B'_{\rm diss}$ by using the criteria that the optical variability timescale $t_{\rm var}$ should be longer or equal to the synchrotron radiation cooling timescale of the electrons that emit optical photons. We test the method with the observations of PSK 1510-089 and BL Lacertae, and derive $R_{\rm diss}<0.15\delta_{\rm D}^{1/3}(1+A)^{2/3}\ $pc for PSK 1510-089 with $t_{\rm var}\sim$ a few hours, and $R_{\rm diss}<0.003\delta_{\rm D}^{1/3}(1+A)^{2/3}\ $pc for BL Lacertae with $t_{\rm var}\sim$ a few minutes. Here $\delta_{\rm D} $ is the Doppler factor and $A$ is the Compton dominance (i.e., the ratio of the Compton to the synchrotron peak luminosities).

The curvature of the $\gamma$-ray spectrum in blazars may reflect the intrinsic distribution of the emitting electron distribution, which will further give some information on the possible acceleration and cooling processes in the emitting region. The $\gamma$-ray spectra of Fermi blazars are normally fitted either by a single power-law (PL) or a log-normal (call Logarithmic Parabola, LP) form. The possible reason for this differnece is not unclear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3LAC Clean sample. We find that the $\gamma$-ray flux (100 MeV-100 GeV) and variability index follow bimodal distributions for PL and LP blazars, where $\gamma$-ray flux and variability index show {a positive correlation}. However, the distributions of the $\gamma$-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of $\gamma$-ray flux for LP and PL blazars may be not intrinsic and all blazars may have an intrinsic curved $\gamma$-ray spectrum and the PL spectrum is just caused by the fitting effect due to the less photons.

We present our photometric monitoring of a flat spectrum radio quasar (FSRQ) 3C 454.3 at Yunnan observatories from 2006 to 2011. We find that the optical color of 3C 454.3 shows obvious redder-when-brighter trend, which reaches a saturation stage when the source is brighter than 15.15 mag at V band. We perform a simulation with multiple values of disk luminosity and spectral index to reproduce the magnitude-color diagram. The results show that the contamination caused by the disk radiation alone is difficult to produce the observed color variability. The variability properties during the outburst in December 2009 are also compared with $\gamma$-ray data derived from Fermi $\gamma$-ray space telescope. The flux variation of these two bands follow a linear relation with $F_{\gamma} \propto F_R^{1.14\pm0.07}$, which provides an observational evidence for external Compton process in 3C 454.3. Meanwhile, this flux correlation indicates that electron injection is the main mechanism for variability origin. We also explore the variation of the flux ratio $F_{\gamma}/F_R$ and the detailed structures in the lightcurves, and discuss some possible origins for the detailed variability behaviors.

We analytically investigate the Schwinger pair production in the de Sitter background by using {\em the uniform asymptotic approximation method}, and show that the equation of motion in general has two turning points, and the nature of these points could be single, double, real or complex, depending on the choice of the free parameters involved in the theory. Different natures of these points lead to different electric currents. In particular, when $\beta \equiv m^2/H^2-9/4$ is positive, both turning points are complex, and the electric current due to the Schwinger process is highly suppressed, where $m$ and $H$ denote, respectively, the mass of the particle and the Hubble parameter. For the turning points to be real, it is necessary to have $\beta < 0$, and the more negative of $\beta$, the easier to produce particles. In addition, when $\beta < 0$, we also study the particle production when the electric field $E$ is very weak. We find that the electric current in this case is proportional to $E^{1/2 - \sqrt{|\beta|}}$, which is strongly enhanced in the weak electric field limit when $m < \sqrt{2} H$.

NGC 4993, as a host galaxy of the electromagnetic counterpart of the first gravitational-wave detection of a binary neutron-star merger, was observed by many powerful telescopes from radio to $\gamma$-ray waveband. The weak nuclear activities of NGC 4993 suggest that it is a low-luminosity active galactic nuclear (LLAGN). We build the multi-waveband spectral energy distributions (SEDs) of NGC 4993 from literatures. We find that the radio spectrum at $\sim 100-300$ GHz is much steeper than that of low-frequency waveband (e.g., 6-100 GHz), where this break was also found in the supermassive black holes in our galaxy center (Sgr A*), and in some other nearby AGNs. The radio emission above and below this break may has different physical origins, which provide an opportunity to probe the accretion and jet properties. We model the multi-waveband SEDs of NGC 4993 with an advection-dominated accretion flow (ADAF)-jet model. We find that the high-frequency steep radio emission at millimeter waveband is consistent with the prediction of the ADAF, while the low-frequency flat radio spectrum is better fitted by the jet. Further more, the X-ray emission can be also explained by the ADAF model simultaneously. From the model fits, we estimate important parameters of the central engine (e.g., accretion rate near the horizon of the black hole and mass-loss rate in jet) for NGC 4993. This result strengthen that the millimiter, submillimer and deep X-ray observations are crucial to understand the weak or quiescent activities in supermassive black-hole systems. Further simultaneous millimeter and X-ray monitoring of this kind of LLAGNs will help us to better understand the physical origin of multiwaveband emission.

We model simultaneous or quasi-simultaneous multi-band spectral energy distributions (SEDs) for a sample of 25 blazars that have radio core-shift measurements, where a one-zone leptonic model and Markov chain Monte Carlo technique are adopted. In the SED fitting for 23 low-synchrotron-peaked (LSP) blazars, the seed photons from broad line region (BLR) and molecular torus are considered respectively in the external Compton process. We find that the SED fitting with the seed photons from the torus are better than those utilizing BLR photons, which suggests that the $\gamma$-ray emitting region may be located outside the BLR. Assuming the magnetic field strength in the $\gamma$-ray emitting region as constrained from the SED fitting follows the magnetic field distribution as derived from the radio core-shift measurements (i.e., $B(R)\simeq B_{\rm 1 pc}(R/\rm 1pc)^{-1}$, $R$ is the distance from the central engine and $B_{\rm 1 pc}$ is the magnetic field strength at 1 pc), we further calculate the location of the $\gamma$-ray emitting region, $R_{\gamma}$, for these blazars. We find that $R_{\gamma}\sim2\times10^{4}R_{\rm S}\simeq10 R_{\rm BLR}$ ($R_{\rm S}$ is the Schwarzschild radius and $R_{\rm BLR}$ is the BLR size), where $R_{\rm BLR}$ is estimated from the broad line luminosities using the empirical correlations obtained using the reverberation mapping methods.

NGC 3115 is the nearest galaxy hosting a billion solar mass black hole and is also a low-luminosity active galactic nucleus (LLAGN). X-ray observations of this LLAGN are able to spatially resolve the hot gas within the sphere of gravitational influence of the supermassive black hole. These observations make NGC 3115 an important testbed for black hole accretion theory in galactic nuclei since they constrain the outer boundary conditions of the hot accretion flow. We present a compilation of the multiwavelength spectral energy distribution (SED) of the nucleus of NGC 3115 from radio to X-rays. We report the results from modeling the observed SED with radiatively inefficient accretion flow (RIAF) models. The radio emission can be well-explained by synchrotron emission from the RIAF without the need for contribution from a relativistic jet. We obtain a tight constraint on the RIAF density profile, $\rho(r) \propto r^{-0.73_{-0.02} ^{+0.01}}$, implying that mass-loss through subrelativistic outflows from the RIAF is significant. The lower frequency radio observation requires the synchrotron emission from a nonthermal electron population in the RIAF, similarly to Sgr A*.

In this paper, we study the pre-inflationary dynamics for the power-law potential $(V(\phi) \propto \phi^n)$ with $n<2$ in the framework of loop quantum cosmology. In the case where the kinetic energy of the inflaton dominates at the initial, the evolution of the universe can always be divided into three different phases prior to preheating: {\em bouncing, transition and slow-roll inflation}. During the bouncing phase, the evolution of the expansion factor is independent not only on the initial conditions but also the inflationary potentials, and is given explicitly by an analytical solution. In contrast, for the potential energy dominated initial conditions, this universality is lost. We also obtain total number of e-folds during the slow-roll inflation, whereby physically viable models are identified. In addition, we present phase space analysis for the inflationary potentials under consideration and compare our results with the ones obtained previously for different potentials.

GRB 170817A is a weak short gamma-ray burst (GRB) accompanied by the gravitational wave (GW) event GW170817. It is believed, that an off beaming relativistic jet, produces this weak GRB. Here we use the $E_{\rm p,i}- E_{\rm iso}$ and $\Gamma - E_{\rm iso}$ relations to determine the Lorentz factor $\Gamma$ and the viewing angle from the edge of the jet $\theta'_{\rm obs}$ of GRB 170817A. Our result shows, $\Gamma = 13.4 ^{+9.8}_{-5.5}$ and $\theta'_{\rm obs} = 4.3^{+1.8}_{-1.5} \,^{\circ}$. This corresponds to the on-axis $E_{\rm p,i} = 415^{+361}_{-167}$ keV and $E_{\rm iso} = (2.4^{+1.6}_{-1.9}) \times 10^{47}$ ergs of the GRB, which is an intrinsically weak short GRB. Interestingly, the Doppler factor and the luminosity follow a universal relation from GRBs and blazars, which indicates they may share similar radiation mechanism.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

The slow-roll inflation for a single scalar field that couples to the Gauss-Bonnet (GB) term represents an important higher-order curvature correction inspired by string theory. With the arrival of the era of precision cosmology, it is expected that the high-order corrections become more and more important. In this paper we study the observational predictions of the slow-roll inflation with the GB term by using the third-order uniform asymptotic approximation method. We calculate explicitly the primordial power spectra, spectral indices, running of the spectral indices for both scalar and tensor perturbations, and the ratio between tensor and scalar spectra. These expressions are all written in terms of the Hubble and GB coupling flow parameters and expanded up to the next-to-leading order in the slow-roll expansions. The upper bounds of errors of the approximations at the third-order are $0.15\%$, so they represent the most accurate results obtained so far in the literature. We expect that the understanding of the GB corrections in the primordial spectra and their constraints by forthcoming observational data will provide clues for the UV complete theory of quantum gravity, such as the string/M-theory.

The millimeter bump, as found in high-resolution multi-waveband observations of M87, most possibly comes from the synchrotron emission of thermal electrons in advection dominated accretion flow(ADAF). It is possible to constrain the accretion rate near the horizon if both the nuclear millimeter emission and its polarization are produced by the hot plasma in the accretion flow. The jet power of M87 has been extensively explored, which is around $8_{\rm -3}^{+7}\times10^{42} {\rm erg/s}$ based on the analysis of the X-ray cavity. The black hole(BH) spin can be estimated if the jet power and the accretion rate near the horizon are known. We model the multi-wavelength spectral energy distribution (SED) of M87 with a coupled ADAF-jet model surrounding a Kerr BH, where the full set of relativistic hydrodynamical equations of the ADAF are solved. The hybrid jet formation model, as a variant of Blandford-Znajek model, is used to model the jet power. We find that the SMBH should be fast rotating with a dimensionless spin parameter $a_{*}\simeq0.98_{\rm -0.02}^{+0.012}$.

Black holes (BHs) surrounded by accretion disks are present in the Universe at different scales of masses, from microquasars up to the active galactic nuclei (AGNs). The current picture remains ad hoc due to the complexity of the magnetic field action. In addition, accretion disks at high Eddington rates can be radiation-pressure dominated and, according to some of the heating prescriptions, thermally unstable. The observational verification of their resulting variability patterns may shed light on both the role of radiation pressure and magnetic fields. We compute the structure and time evolution of an accretion disk, using the code GLADIS (which models the global accretion disk instability). We supplement this model with a modified viscosity prescription, which can to some extent describe the magnetisation of the disk. We study the results for a large grid of models, to cover the whole parameter space, and we derive conclusions separately for different scales of black hole masses, which are characteristic for various types of cosmic sources. We show the dependencies between the flare or outburst duration, its amplitude, and period, on the accretion rate and viscosity scaling. We show that if the heating rate in the accretion disk grows more rapidly with the total pressure and temperature, the instability results in longer and sharper flares. In general, we confirm that the disks around the supermassive black holes are more radiation-pressure dominated and present relatively brighter bursts. Our method can also be used as an independent tool for the black hole mass determination, which we confront now for the intermediate black hole in the source HLX-1.

The phase diagram of hydrogen is one of the most important challenges in high-pressure physics and astrophysics. Especially, the melting of dense hydrogen is complicated by dimer dissociation, metallization and nuclear quantum effect of protons, which together lead to a cold melting of dense hydrogen when above 500 GPa. Nonetheless, the variation of the melting curve at higher pressures is virtually uncharted. Here we report that using ab initio molecular dynamics and path integral simulations based on density functional theory, a new atomic phase is discovered, which gives an uplifting melting curve of dense hydrogen when beyond 2 TPa, and results in a reentrant solid-liquid transition before entering the Wigner crystalline phase of protons. The findings greatly extend the phase diagram of dense hydrogen, and put metallic hydrogen into the group of alkali metals, with its melting curve closely resembling those of lithium and sodium.

The hyperluminous X-ray source (HLX-1, the peak X-ray luminosity $\sim 10^{42}\rm erg\ s^{-1}$) near the spiral galaxy ESO 243-49 is possibly the best candidate for intermediate mass black hole (IMBH), which underwent recurrent outbursts with a period of $\sim 400$ days. The physical reason for this quasi-periodic variability is still unclear. We explore the possibility of radiation-pressure instability in accretion disk by modeling the light curve of HLX-1, and find that it can roughly reproduce the duration, period and amplitude of the recurrent outbursts HLX-1 with an IMBH of ~10^5Msun. Our result provides a possible mechanism to explain the recurrent outbursts in HLX-1. We further find a universal correlation between the outburst duration and the bolometric luminosity for the BH sources with a very broad mass range (e.g., X-ray binaries, XRBs, HLX-1 and active galactic nuclei, AGNs), which is roughly consistent with the prediction of radiation-pressure instability of the accretion disk. These results imply that "heartbeat" oscillations triggered by radiation-pressure instability may appears in different-scale BH systems.

HLX-1, currently the best intermediate-mass black hole candidate, has undergone seven violent outbursts, each with a peak X-ray luminosity of $L_{\mathrm{peak},\mathrm{X}}\sim 10^{42}\ \rm{erg\ s^{-1}}$. Interestingly, the properties of the HLX-1 outbursts evolve with time. In this work, we aim to constrain the physical parameters of the central engine of the HLX-1 outbursts in the framework of the black hole accretion. We find that the physical properties of the HLX-1 outbursts are consistent with being driven by the radiation pressure instability. This scenario can explain the evolution of the recurrent timescales of the HLX-1 outbursts as a function of the durations.

M87 is arguably the best supermassive black hole (BH) to explore the jet and/or accretion physics due to its proximity and fruitful high-resolution multi-waveband observations. We model the multi-wavelength spectral energy distribution (SED) of M87 core that observed at a scale of 0.4 arcsec ($\sim 10^5R_{\rm g}$, $R_{\rm g}$ is gravitational radius) as recently presented by Prieto et al. Similar to Sgr A*, we find that the millimeter bump as observed by Atacama Large Millimeter/submillimeter Array (ALMA) can be modeled by the synchrotron emission of the thermal electrons in advection dominated accretion flow (ADAF), while the low-frequency radio emission and X-ray emission may dominantly come from the jet. The millimeter radiation from ADAF dominantly come from the region within $10R_{\rm g}$, which is roughly consistent with the recent very long baseline interferometry observations at 230\,GHz. We further calculate the Faraday rotation measure (RM) from both ADAF and jet models, and find that the RM predicted from the ADAF is roughly consistent with the measured value while the RM predicted from the jet is much higher if jet velocity close to the BH is low or moderate (e.g., $v_{\rm jet}\lesssim0.6\,c$). With the constraints from the SED modeling and RM, we find that the accretion rate close to the BH horizon is $\sim (0.2-1)\times10^{-3}{M}_{\odot} \rm yr^{-1}\ll\dot{\it M}_{\rm B}\sim 0.2\it {M}_{\odot} \rm yr^{-1}$ ($\dot{M}_{\rm B}$ is Bondi accretion rate), where the electron density profile, $n_{\rm e}\propto r^{\sim -1}$, in the accretion flow is consistent with that determined from X-ray observation inside the Bondi radius and recent numerical simulations.

PKS 1424+240 is a distant very high energy gamma-ray BL Lac object with redshift $z=0.601$. It was found that pure synchrotron self-Compton (SSC) process normally need extreme input parameters (e.g., very low magnetic field intensity and extraordinarily large Doppler factor) to explain its multi-wavelength spectral energy distributions (SEDs). To avoid the extreme model parameters, different models have been proposed (e.g., two-zone SSC model or lepto-hadronic model). In this work, we employ the traditional one-zone leptonic model after including a weak external Compton component to re-explore the simultaneous multi-wavelength SEDs of PKS 1424+240 in both high (2009) and low (2013) states. We find that the input parameters of magnetic field and Doppler factor are roughly consistent with those of other BL Lacs if a weak external photon field from either broad line region (BLR) or the dust torus. However, the required energy density of seed photons from BLR or torus is about 3 orders of magnitude less than that constrained in luminous quasars (e.g., flat-spectrum radio quasars, FSRQs). This result suggests that the BLR/torus in BL Lacs is much weaker than that of luminous FSRQs (but not fully disappear), and the inverse-Compton of external photons from BLR/torus may still play a role even in high synchrotron peaked blazars.

Previous studies have shown that WISE-selected hyperluminous, hot dust-obscured galaxies (Hot DOGs) are powered by highly dust-obscured, possibly Compton-thick AGNs. High obscuration provides us a good chance to study the host morphology of the most luminous AGNs directly. We analyze the host morphology of 18 Hot DOGs at $z\sim3$ using Hubble Space Telescope/WFC3 imaging. We find that Hot DOGs have a high merger fraction ($62\pm 14 \%$). By fitting the surface brightness profiles, we find that the distribution of S\'ersic indices in our Hot DOG sample peaks around 2, which suggests that most of Hot DOGs have transforming morphologies. We also derive the AGN bolometric luminosity ($\sim10^{14}L_\odot$) of our Hot DOG sample by using IR SEDs decomposition. The derived merger fraction and AGN bolometric luminosity relation is well consistent with the variability-based model prediction (Hickox et al. 2014). Both the high merger fraction in IR-luminous AGN sample and relatively low merger fraction in UV/optical-selected, unobscured AGN sample can be expected in the merger-driven evolutionary model. Finally, we conclude that Hot DOGs are merger-driven and may represent a transit phase during the evolution of massive galaxies, transforming from the dusty starburst dominated phase to the unobscured QSO phase.

We first derive the primordial power spectra, spectral indices and runnings of both scalar and tensor perturbations of a flat inflationary universe to the second-order approximations of the slow-roll parameters, in the framework of loop quantum cosmology with the inverse-volume quantum corrections. This represents an extension of our previous work in which the parameter $\sigma$ was assumed to be an integer, where $\sigma$ characterizes the quantum corrections and in general can take any of values from the range $\sigma \in (0, 6]$. Restricting to the first-order approximations of the slow-roll parameters, we find corrections to the results obtained previously in the literature, and point out the causes for such errors. To our best knowledge, these represent the most accurate calculations of scalar and tensor perturbations given so far in the literature. Then, fitting the perturbations to the recently released data by Planck (2015), we obtain the most severe constraints for various values of $\sigma$. Using these constraints as our referring point, we discuss whether these quantum gravitational corrections can lead to measurable signatures in the future cosmological observations. We show that, depending on the value of $\sigma$, the scale-dependent contributions to the relativistic inflationary spectra due to the inverse-volume corrections could be well within the range of the detectability of the forthcoming generations of experiments, such as the Stage IV experiments.

Loop quantum cosmology (LQC) provides promising resolutions to the trans-Planckian issue and initial singularity arising in the inflationary models of general relativity. In general, due to different quantization approaches, LQC involves two types of quantum corrections, the holonomy and inverse-volume, to both of the cosmological background evolution and perturbations. In this paper, using {\em the third-order uniform asymptotic approximations}, we derive explicitly the observational quantities of the slow-roll inflation in the framework of LQC with these quantum corrections. We calculate the power spectra, spectral indices, and running of the spectral indices for both scalar and tensor perturbations, whereby the tensor-to-scalar ratio is obtained. We expand all the observables at the time when the inflationary mode crosses the Hubble horizon. As the upper error bounds for the uniform asymptotic approximation at the third-order are $\lesssim 0.15\%$, these results represent the most accurate results obtained so far in the literature. It is also shown that with the inverse-volume corrections, both scalar and tensor spectra exhibit a deviation from the usual shape at large scales. Then, using the Planck, BAO and SN data we obtain new constraints on quantum gravitational effects from LQC corrections, and find that such effects could be within the detection of the forthcoming experiments.

It is very controversial whether radio--X-ray correlation as defined in LH state of XRBs can extend to quiescent state (e.g., X-ray luminosity less than a critical value of $L_{\rm X,c} \sim10^{-5.5}L_{\rm Edd}$) or not. In this work, we collect a sample of XRBs and low luminosity active galactic nuclei (LLAGNs) with wide distribution of Eddington ratios to reexplore the fundamental plane between 5 GHz radio luminosity, $L_{\rm R}$, 2-10 keV X-ray luminosity, $L_{\rm X}$, and black hole (BH) mass, $M_{\rm BH}$, namely $\log L_{\rm R}=\xi_{\rm X} \log L_{\rm X}+\xi_{\rm M}\log M_{\rm BH}+\rm constant$. For the whole sample, we confirm the former fundamental plane of Merloni et al. and Falcke et al. that $\xi_{\rm X}\sim 0.6$ and $\xi_{\rm M}\sim 0.8$ even after including more quiescent BHs. The quiescent BHs follow the fundamental plane very well, and, however, FR I radio galaxies follow a steeper track comparing other BH sources. After excluding FR Is, we investigate the fundamental plane for BHs in quiescent state with $L_{\rm X}< L_{\rm X,c}$ and sub-Eddington BHs with $L_{\rm X}> L_{\rm X,c}$ respectively, and both subsamples have a similar slope, $\xi_{\rm X}\sim0.6$, which support that quiescent BHs may behave similar to those in low-hard state. We further select two subsamples of AGNs with BH mass in a narrow range (FR Is with $M_{\rm BH}=10^{8.8\pm0.4}$ and other LLAGNs with $M_{\rm BH}=10^{8.0\pm0.4}$) to simulate the behavior of a single supermassive BH evolving from sub-Eddington to quiescent state. We find that the highly sub-Eddington sources with $L_{\rm X}/L_{\rm Edd}\sim10^{-6}-10^{-9}$ still roughly stay on the extension of radio--X-ray correlation as defined by other sub-Eddington BHs. Our results are consistent with several recent observations in XRBs that the radio--X-ray correlation as defined in low-hard state can extend to highly sub-Eddington quiescent state.

Both gamma-ray bursts (GRBs) and blazars have relativistic jets pointing at a small angle from our line of sight. Several recent studies suggested that these two kinds of sources may share similar jet physics. In this work, we explore the variability properties for GRBs and blazars as a whole. We find that the correlation between minimum variability timescale (MTS) and Lorentz factor, $\Gamma$, as found only in GRBs by Sonbas et al. can be extended to blazars with a joint correlation of $\rm MTS\propto\Gamma^{-4.7\pm0.3}$. The same applies to the $\rm MTS\propto \it L_{\gamma}^{\rm -1.0\pm0.1}$ correlation as found in GRBs, which can be well extended into blazars as well. These results provide further evidence that the jets in these two kinds of sources are similar despite of the very different mass scale of their central engines. Further investigations of the physical origin of these correlations are needed, which can shed light on the nature of the jet physics.

The X-ray Timing and Polarization (XTP) is a mission concept for a future space borne X-ray observatory and is currently selected for early phase study. We present a new design of X-ray polarimeter based on the time projection gas chamber. The polarimeter, placed above the focal plane, has an additional rear window that allows hard X-rays to penetrate (a transmission of nearly 80% at 6 keV) through it and reach the detector on the focal plane. Such a design is to compensate the low detection efficiency of gas detectors, at a low cost of sensitivity, and can maximize the science return of multilayer hard X-ray telescopes without the risk of moving focal plane instruments. The sensitivity in terms of minimum detectable polarization, based on current instrument configuration, is expected to be 3% for a 1mCrab source given an observing time of 10^5 s. We present preliminary test results, including photoelectron tracks and modulation curves, using a test chamber and polarized X-ray sources in the lab.

We derive the primordial power spectra and spectral indexes of the density fluctuations and gravitational waves in the framework of loop quantum cosmology (LQC) with holonomy and inverse-volume corrections, by using the uniform asymptotic approximation method to its third-order, at which the upper error bounds are $\lesssim 0.15\%$, and accurate enough for the current and forthcoming cosmological observations. Then, using the Planck, BAO and SN data we obtain the tightest constraints on quantum gravitational effects from LQC corrections, and find that such effects could be well within the detection of the current and forthcoming cosmological observations.

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

The well-known correlation between the radio luminosity ($L_R$) and the X-ray luminosity ($L_X$) $L_R / L_X \simeq 10^{-5}$ holds for a variety of objects like active galactic nuclei, galactic black holes, solar flares and cool stars. Here we extend the relation to gamma-ray bursts (GRBs), and find the GRBs also lay on the same $L_R-L_X$ relation, with a slightly different slope as $L_R \propto L_X^{1.1}$. This relation implies the explosions in different scales may have a common underlying origin.

We fit the (quasi-)simultaneous multi-waveband spectral energy distributions (SEDs) for a sample of low-synchrotron-peaked (LSP) blazars with a one-zone leptonic model. The seed photons that predominantly come from broad line region (BLR) and infrared (IR) molecular torus are considered respectively in external Compton process. We find that the modeling with IR seed photons is systematically better than that with BLR photons based on a $\chi^2$ test, which suggest that $\gamma$-ray emitting region most possibly stay outside the BLR. The minimum electron Lorentz factor, $\gamma_{\rm min}$, is constrained from the modeling for these LSP blazars with good soft X-ray data (ranges from 5 to 160 with a median value of 55), which plays a key role in jet power estimation. Assuming one-to-one ratio of proton and electron, we find that the jet power for LSP blazars is systematically higher than that of FR II radio galaxies at given 151 MHz radio luminosity, $L_{\rm 151MHz}$, even though FR IIs are regarded as same as LSP blazars in unification scheme except the jet viewing angle. The possible reason is that there are some $e^{\pm}$ pairs in the jet of these blazars. If this is the case, we find the number density of $e^{\pm}$ pairs should be several times higher than that of $e^{-}-p$ pairs by assuming the jet power is the same for LSP blazars and FR IIs at given $L_{\rm 151MHz}$.

Recently an `outliers' track of radio-X-ray correlation was found, which is much steeper than the former universal correlation, where dual tracks were speculated to be triggered by different accretion processes. In this work, we test this issue by exploring hard X-ray spectral evolution in four black-hole X-ray binaries (XRBs) with multiple, quasi-simultaneous radio and X-ray observations. Firstly, we find that hard X-ray photon indices, $\Gamma$, are anti- and positively correlated to X-ray fluxes when the X-ray flux, $F_{\rm 3-9keV}$, is below and above a critical flux, $F_{\rm X,crit}$, which are consistent with prediction of advection dominated accretion flow (ADAF) and disk-corona model respectively. Secondly and most importantly, we find that the radio-X-ray correlations are also clearly different when the X-ray fluxes are higher and lower than the critical flux that defined by X-ray spectral evolution. The data points with $F_{\rm 3-9keV}\gtrsim F_{\rm X,crit}$ have a steeper radio-X-ray correlation ($F_{\rm X}\propto F_{\rm R}^{b}$ and $b\sim 1.1-1.4$), which roughly form the `outliers' track. However, the data points with anti-correlation of $\Gamma-F_{\rm 3-9keV}$ either stay in the universal track with $b\sim0.61$ or stay in transition track (from the universal to `outliers' tracks or vice versa). Therefore, our results support that the universal and `outliers' tracks of radio-X-ray correlations are regulated by radiatively inefficient and radiatively efficient accretion model respectively.

In recent years, it was found that there are several low/hard state of X-ray binaries (XRBs) follow an `outliers' track of radio--X-ray correlation ($L_{\rm R}\propto L_{\rm X}^{b}$ and $b\sim1.4$), which is much steeper than the former universal track with $b\sim0.6$. In this work, we compile a sample of bright radio-quiet active galactic nuclei (AGNs) and find that their hard X-ray photon indices and Eddington ratios are positively correlated, which is similar to that of `outliers' of XRBs, where both bright AGNs and `outliers' of XRBs have bolometric Eddington ratios $\gtrsim1%L_{\rm Edd}$ ($L_{\rm Edd}$ is Eddington luminosity). The Eddington-scaled radio--X-ray correlation of these AGNs is also similar to that of `outliers' of XRBs, which has a form of $L_{\rm 5 GHz}/L_{\rm Edd}\propto (L_{\rm 2-10 keV}/L_{\rm Edd})^{c}$ with $c\simeq1.59$ and 1.53 for AGNs and XRBs respectively. Both the positively correlated X-ray spectral evolution and the steeper radio--X-ray correlation can be regulated by a radiatively efficient accretion flow (e.g., disk-corona). Based on these similarities, we further present a new fundamental plane for `outliers' of XRBs and bright AGNs in black-hole (BH) mass, radio and X-ray luminosity space: $\log L_{\rm R}=1.59^{+0.28}_{-0.22} \log L_{\rm X}- 0.22^{+0.19}_{-0.20}\log M_{\rm BH}-28.97^{+0.45}_{-0.45}$ with a scatter of $\sigma_{\rm R}=0.51\rm dex$. This fundamental plane is suitable for radiatively efficient BH sources, while the former plane proposed by Merloni et al. and Falcke et al. may be most suitable for radiatively inefficient sources.

In this paper, we study static vacuum solutions of quantum gravity at a fixed Lifshitz point in (2+1) dimensions, and present all the diagonal solutions in closed forms in the infrared limit. The exact solutions represent spacetimes with very rich structures: they can represent generalized BTZ black holes, Lifshitz space-times or Lifshitz solitons, in which the spacetimes are free of any kind of space-time singularities, depending on the choices of the free parameters of the solutions. We also find several classes of exact static non-diagonal solutions, which represent similar space-time structures as those given in the diagonal case. The relevance of these solutions to the non-relativistic Lifshitz-type gauge/gravity duality is discussed.

Wu et al. (2009) suggested an empirical relation between the magnetic field strength along the line of sight (LOS) and the dispersion of Faraday rotation measure (RM) distribution in turbulent media with root-mean-square sonic Mach number $M_s \simeq 1$. In this paper, we extend the work by incorporating the Mach number dependence. Media with $\sim 0.5 < M_s < \sim 2$ are considered to cover the Mach number range of the warm ionized medium (WIM) of our Galaxy. Three-dimensional, magnetohydrodynamic isothermal turbulence simulations with solenoidal forcing are used. We suggest a new relation among the LOS magnetic field strength, the dispersion of RM distribution, and the Mach number, which approximately represents the relation for Alfv\'enic perturbations. In addition, we suggest a relation between the Mach number and the dispersion of log-normal distribution of emission measure (EM), which is basically the relation for the Mach number and the density dispersion. The relations could be used for a quick and rough estimation of the LOS magnetic field strength in the turbulent WIM.

The universal radio--X-ray correlation (F_R \propto F_X^{b}, b~0.5-0.7) has been found for a sample of black-hole X-ray binaries (BHBs) in their low/hard states, which can roughly be explained by the coupled model of jet and radiatively inefficient advection dominated accretion flow. However, more and more `outliers' were found in last few years, which evidently deviate from the universal radio-X-ray correlation and usually show a much steeper correlation with an index of ~1.4. Based on simple physical assumptions, the radiatively efficient accretion flows are speculated to exist in these `outliers'. In this work, we test this issue by modelling the `outliers' track based on the radiatively efficient disc-corona model and the hybrid jet model. We find that our model predicts a steeper radio--X-ray correlation with slopes >~1.2 for the typical viscosity parameter of \alpha~0.05-0.2. In particular, the slope is ~1.4 for the case of \alpha~0.1, which is consistent with the observational results of H1743-322 very well. Our results suggest that the `outliers' track may be regulated by the disc-corona model.

In this paper, we study the effects of parity violation on non-gaussianities of primordial gravitational waves in the framework of Ho\v{r}ava-Lifshitz theory of gravity, in which high-order spatial derivative operators, including the ones violating parity, generically appear. By calculating the three point function, we find that the leading-order contributions to the non-gaussianities come from the usual second-order derivative terms, which produce the same bispectrum as that found in general relativity. The contributions from high-order spatial n-th derivative terms are always suppressed by a factor $(H/M_*)^{n-2} \; (n \ge 3)$, where $H$ denotes the inflationary energy and $M_*$ the suppression mass scale of the high-order spatial derivative operators of the theory. Therefore, the next leading-order contributions come from the 3-dimensional gravitational Chern-Simons term. With some reasonable arguments, it is shown that this 3-dimensional operator is the only one that violates the parity and in the meantime has non-vanishing contributions to non-gaussianities.

We fit the most updated broadband spectral energy distribution from radio to X-rays for NGC 4258 with a coupled accretion-jet model that surrounding a Kerr black hole (BH), where both the jet and the warped H_2O maser disk are assumed to be triggered by a spinning BH through Blandford-Znajek mechanism and Bardeen-Petterson effect respectively. The accretion flow consists with an inner radiatively inefficient accretion flow (RIAF) and an outer truncated standard thin disk, where the transition radius R_tr~3*10^3Rg for NGC 4258 based on the width and variability of its narrow Fe K$\alpha$ line. The hybrid jet formation model, as a variant of Blandford-Znajek model, is used to model the jet power. Therefore, we can estimate the accretion rate and BH spin through the two observed quantities--X-ray emission and jet power, where the observed jet power is estimated from the low-frequency radio emission. Through this method, we find that the BH of NGC 4258 should be mildly spinning with dimensionless spin parameter a_*=0.7\pm0.2. The outer thin disk mainly radiates at near infrared waveband and the jet contributes predominantly at radio waveband. Using above estimated BH spin and the inferred accretion rate at the region of the maser disk based on the physical existence of the H_2 O maser, we find that the warp radius is ~8.6*10^4 R_g if it is driven by the Bardeen-Petterson effect, which is consistent with the observational result very well.

We study primordial gravitational waves (PGWs) in the Horava-Lifshitz (HL) theory of quantum gravity, in which high-order spatial derivative operators, including the ones violating parity, generically appear in order for the theory to be power-counting renormalizable and ultraviolet (UV) complete. Because of both parity violation and non-adiabatic evolution of the modes due to a modified dispersion relationship, a large polarization of PGWs becomes possible, and it could be well within the range of detection of the BB, TB and EB power spectra of the forthcoming cosmic microwave background (CMB) observations.

Recent observations suggest that in black hole X-ray binaries jet/outflow formation is related to the hot plasma in the vicinity of the black hole, either in the form of an advection-dominated accretion flow at low accretion rates or in a disk corona at high accretion rates. We test the viability of this scenario for supermassive black holes using two samples of active galactic nuclei distinguished by the presence (radio-strong) and absence (radio-weak) of well-collimated, relativistic jets. Each is centered on a narrow range of black hole mass but spans a very broad range of Eddington ratios, effectively simulating, in a statistical manner, the behavior of a single black hole evolving across a wide spread in accretion states. Unlike the relationship between the radio and optical luminosity, which shows an abruptly break between high- and low-luminosity sources at an Eddington ratio of ~1, the radio emission-a measure of the jet power-varies continuously with the hard X-ray (2-10 keV) luminosity, roughly as L_R \propto L_X^(0.6-0.75). This relation, which holds for both radio-weak and radio-strong active galaxies, is similar to the one seen in X-ray binaries. Jet/outflow formation appears to be closely linked to the conditions that give rise to the hot, optically thin coronal emission associated with accretion flows, both in the regime of low and high accretion rates.

In this paper, we show that the spin-0 gravitons appearing in Horava-Lifshitz gravity without the projectability condition can be eliminated by extending the gauge symmetries of the foliation-preserving diffeomorphisms to include a local U(1) symmetry. As a result, the problems of stability, ghost, strong coupling, and different speeds in the gravitational sector are automatically resolved. In addition, with the detailed balance condition softly breaking, the number of independent coupling constants can be significantly reduced (from more than 70 down to 15), while the theory is still UV complete and possesses a healthy IR limit, whereby the prediction powers of the theory are considerably improved. The strong coupling problem in the matter sector can be cured by introducing an energy scale $M_{*}$, so that $M_{*} < \Lambda_{\omega}$, where $M_{*}$ denotes the suppression energy of high order derivative terms, and $\Lambda_{\omega}$ the would-be strong coupling energy scale.

The two components of radio emission, above and below 86 GHz respectively, from the Galactic center source-Sgr A* can be naturally explained by the hybrid of thermal and nonthermal electrons in hot accretion flows (e.g., radiatively inefficient accretion flow, RIAF, e.g., Yuan et al. 2003). We further apply this model to a sample of nearby low-luminosity active galactic nuclei (LLAGNs), which are also believed to be powered by the RIAF. We selected the LLAGNs with only compact radio cores according to the high-resolution radio observations, and the sources observed with jets or jet-like features are excluded. We find that the radio emission of LLAGNs is severely underpredicted by pure RIAF model, and can be naturally explained by the RIAF model with a hybrid electron population consisting of both thermal and nonthermal particles. Our model can roughly reproduce the observed anti-correlation between the mass-corrected radio loudness and Eddington ratio for the LLAGNs in our sample. We further model the spectral energy distributions of each source in our sample, and find that roughly all sources can be well fitted if a small fraction of the steady state electron energy is ejected into the nonthermal electrons. The size of radio emission region of our model is around several thousand gravitational radii, which is also roughly consistent with the recent high-resolution VLBI observations for some nearby LLAGNs.

In this paper, we study inflation in the framework of the nonrelativistic general covariant theory of the Ho\v{r}ava-Lifshitz gravity with the projectability condition and an arbitrary coupling constant $\lambda$. We find that the Friedmann-Robterson-Walker (FRW) universe is necessarily flat in such a setup. We work out explicitly the linear perturbations of the flat FRW universe without specifying to a particular gauge, and find that the perturbations are different from those obtained in general relativity, because of the presence of the high-order spatial derivative terms. Applied the general formulas to a single scalar field, we show that in the sub-horizon regions, the metric and scalar field are tightly coupled and have the same oscillating frequencies. In the super-horizon regions, the perturbations become adiabatic, and the comoving curvature perturbation is constant. We also calculate the power spectra and indices of both the scalar and tensor perturbations, and express them explicitly in terms of the slow roll parameters and the coupling constants of the high-order spatial derivative terms. In particular, we find that the perturbations, of both scalar and tensor, are almost scale-invariant, and the spectrum indices are the same as those given in GR, but the ratio of the scalar and tensor power spectra depends on the high-order spatial derivative terms, and can be different from that of GR significantly.

The bulk Lorentz factor of the gamma-ray burst (GRB) ejecta (Gamma_0) is a key parameter to understand the GRB physics. Liang et al. have discovered a correlation between Gamma_0 and isotropic gamma-ray energy: Gamma_0\proptoE_{gamma,iso,52}^{0.25}. By including more GRBs with updated data and more methods to derive Gamma_0, we confirm this correlation and obtain Gamma_0~91E_{gamma,iso,52}^{0.29}. Evaluating the mean isotropic gamma-ray luminosities L_{gamma,iso} of the GRBs in the same sample, we discover an even tighter correlation Gamma_0~249 L_{gamma,iso,52}^{0.30}. We propose an interpretation to this later correlation. Invoking a neutrino-cooled hyperaccretion disk around a stellar mass black hole as the central engine of GRBs, we derive jet luminosity powered by neutrino annihilation and baryon loading from a neutrino-driven wind. Applying beaming correction, we finally derive Gamma_0\proptoL_{gamma,iso}^{0.22}, which is well consistent with the data. This suggests that the central engine of long GRBs is likely a stellar mass black hole surrounded by a hyper-accreting disk.

We consider an extended theory of Horava-Lifshitz gravity with the detailed balance condition softly breaking, but without the projectability condition. With the former, the number of independent coupling constants is significantly reduced. With the latter and by extending the original foliation-preserving diffeomorphism symmetry $ {{Diff}}(M, {\cal{F}})$ to include a local U(1) symmetry, the spin-0 gravitons are eliminated. Thus, all the problems related to them disappear, including the instability, strong coupling, and different speeds in the gravitational sector. When the theory couples to a scalar field, we find that the scalar field is not only stable in both the ultraviolet (UV) and infrared (IR), but also free of the strong coupling problem, because of the presence of high-order spatial derivative terms of the scalar field. Furthermore, applying the theory to cosmology, we find that due to the additional U(1) symmetry, the Friedmann-Robertson-Walker (FRW) universe is necessarily flat. We also investigate the scalar, vector, and tensor perturbations of the flat FRW universe, and derive the general linearized field equations for each kind of the perturbations.

We propose a two-component jet model consistent with the observations of several gamma ray bursts (GRBs) and active galactic nuclei (AGNs). The jet consists of inner and outer components, and they are supposed to be driven by the Blandford-Znajek (BZ) and Blandford-Payne (BP) processes, respectively. The baryons in the BP jet is accelerated centrifugally via the magnetic field anchored in the accretion disk. The BZ jet is assumed to be entrained a fraction of accreting matter leaving the inner edge of the accretion disk, and the baryons are accelerated in the conversion from electromagnetic energy to the kinetic energy. By fitting the Lorentz factors of some GRBs (GRB 030329, GRB 051221A, GRB 080413B) and AGNs (Cen A, Mkn 501 and Mkn 421) with this model, we constrain the physical parameters related to the accretion and outflow of these two kind of objects. We conclude that the spine/sheath structure of the jet from these sources can be interpreted naturally by the BZ and BP processes.

We compile 23 Gamma-ray Bursts (GRBs) and 21 blazars with estimated Doppler factors, and the Doppler factors of GRBs are estimated from their Lorentz factors by assuming their jet viewing angles \theta->0 degree. Using the conventional assumption that the prompt emission of GRBs is dominated by the synchrotron radiation, we calculate the synchrotron luminosity of GRBs from their total isotropic energy and burst duration. Intriguingly, we discover a uniform correlation between the synchrotron luminosity and Doppler factor, L_syn \propto D^3.1, for GRBs and blazars, which suggests that they may share some similar jet physics. One possible reason is that GRBs and blazars have, more or less, similar intrinsic synchrotron luminosities and both of them are strongly enhanced by the beaming effect. After Doppler and redshift-correction, we find that the intrinsic peak energy of the GRBs ranges from 0.1 to 3 keV with a typical value of 1 keV. We further correct the beaming effect for the observed luminosity of GRBs and find that there exists a positive correlation between the intrinsic synchrotron luminosity and peak energy for GRBs, which is similar to that of blazars. Our results suggest that both the intrinsic positive correlation and the beaming effect may be responsible for the observed tight correlation between the isotropic energy and the peak energy in GRBs (so called "Amati" relation).

We study the strong coupling problem in the Horava-Melby-Thompson setup of the Horava-Lifshitz gravity with an arbitrary coupling constant $\lambda$, generalized recently by da Silva, where $\lambda$ describes the deviation of the theory in the infrared from general relativity that has $\lambda_{GR} = 1$. We find that a scalar field in the Minkowski background becomes strong coupling for processes with energy higher than $\Lambda_{\omega} [\equiv (M_{pl}/c_1)^{3/2} M_{pl}|\lambda - 1|^{5/4}]$, where generically $c_1 \ll M_{pl}$. However, this problem can be cured by introducing a new energy scale $M_{*}$, so that $M_{*} < \Lambda_{\omega}$, where $M_{*}$ denotes the suppression energy of high order derivative terms of the theory.

We investigate the relation between black hole mass, M_bh, and jet power, Q_jet, for a sample of BL Lacs and radio quasars. We find that BL Lacs are separated from radio quasars by the FR I/II dividing line in M_bh-Q_jet plane, which strongly supports the unification scheme of FR I/BL Lac and FR II/radio quasar. The Eddington ratio distribution of BL Lacs and radio quasars exhibits a bimodal nature with a rough division at L_bol/L_Edd~0.01, which imply that they may have different accretion modes. We calculate the jet power extracted from advection dominated accretion flow (ADAF), and find that it require dimensionless angular momentum of black hole j~0.9-0.99 to reproduce the dividing line between FR I/II or BL Lac/radio quasar if dimensionless accretion rate mdot=0.01 is adopted, which is required by above bimodal distribution of Eddington ratios. Our results suggest that black holes in radio galaxies are rapidly spinning.

We investigate the correlation between 151 MHz radio luminosity, L_151MHz, and jet power, P_jet, for a sample of low-power radio galaxies, of which the jet power is estimated from X-ray cavities. The jet power for a sample of FR I radio galaxies is estimated with the derived empirical correlation. We find that P_jet/L_Edd is positively correlated with Lx(2-10 keV)/L_Edd for FR Is, where L_Edd is the Eddington luminosity and Lx(2-10 keV) is 2-10 keV X-ray luminosity. We calculate the jet power of a hybrid model, as a variant of Blandford-Znajek model proposed by Meier, based on the global solution of the advection-dominated accretion flow (ADAF) surrounding a Kerr black hole (BH). Our model calculations suggest that the maximal jet power is a function of mass accretion rate and the black hole spin parameter j. The hard X-ray emission is believed to be mainly from the ADAFs in FR Is, and the mass accretion rate is therefore constrained with the X-ray emission in our ADAF model calculations. We find that the dimensionless angular momentum of BH j>~0.9 is required in order to reproduce the observed relation of P_jet/L_Edd-Lx(2-10 keV)/L_Edd for FR Is. Our conclusion will be strengthened if part of the X-ray emission is contributed by the jets. Our results suggest that BHs in FR I radio galaxies are rapidly spinning, which are almost not affected by the uncertainty of the black hole mass estimates.

The stability of de Sitter spacetime in Horava-Lifshitz theory of gravity with projectability but without detailed balance condition is studied. It is found that, in contrast to the case of the Minkowski background, the spin-0 graviton now is stable for any given $\xi$, and free of ghost for $\xi \le 0$ in the infrared limit, where $\xi$ is the dynamical coupling constant.

In this paper, we consider two different issues, stability and strong coupling, raised lately in the newly-proposed Horava-Lifshitz (HL) theory of quantum gravity with projectability condition. We find that all the scalar modes are stable in the de Sitter background, due to two different kinds of effects, one from high-order derivatives of the spacetime curvature, and the other from the exponential expansion of the de Sitter space. Combining these effects properly, one can make the instability found in the Minkowski background never appear even for small-scale modes, provided that the IR limit is sufficiently closed to the relativistic fixed point. At the fixed point, all the modes become stabilized. We also show that the instability of Minkowski spacetime can be cured by introducing mass to the spin-0 graviton. The strong coupling problem is investigated following the effective field theory approach, and found that it cannot be cured by the Blas-Pujolas-Sibiryakov mechanism, initially designed for the case without projectability condition, but might be circumvented by the Vainshtein mechanism, due to the non-linear effects. In fact, we construct a class of exact solutions, and show explicitly that it reduces smoothly to the de Sitter spacetime in the relativistic limit.

We study the radion stability and radion mass in the framework of the Horava-Witten (HW) heterotic M-Theory on $S^{1}/Z_{2}$, and find that the radion is stable and its mass can be of the order of GeV. The gravity is localized on the visible brane, and the spectrum of the gravitational Kaluza-Klein (KK) modes is discrete and can have a mass gap of TeV. The corrections to the 4D Newtonian potential from the higher order gravitational KK modes are exponentially suppressed. Applying such a setup to cosmology, we find the generalized Friedmann-like equations on each of the two orbifold branes.

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

The bolometric luminosities and black hole (BH) masses are estimated by various methods for a sample of young radio galaxies with known ages. We find that the ages are positively correlated with the bolometric luminosities in these young radio galaxies, which is consistent with theoretical prediction based on radiation pressure instability of accretion disk in Czerny et al. The ages of young radio galaxies are also found to be consistent with the theoretical durations of outbursts in BH mass and accretion rate (in Eddington unit) plane, where the outbursts are assumed to be triggered by the radiation pressure instabilities. Our results provide the observational evidence for the radiation pressure instability, which causes limit-cycle behavior, as a physical mechanism that may be responsible for these short-lived young radio galaxies.

We estimate black hole masses and Eddington ratios for a sample of 81 young radio galaxies (42 CSS +39 GPS). We find that the average black hole (BH) mass of these young radio galaxies is <log Mbh>~8.3, which is less than that of radio loud QSOs and low redshift radio galaxies. The CSS/GPS sources have relatively high Eddington ratios, with an average value of <log Lbol/LEdd>=-0.75, which are similar to those of narrow line Seyfert 1 galaxies (NLS1s). This suggests that young radio galaxies may not only be in the early stages of their radio activity, but also in the early stage of their accretion activity. We find that the young radio galaxies as a class deviate systematically from M_bh-\sigma relation defined by nearby inactive galaxies, when using [O III] as a surrogate for stellar velocity dispersion, \sigma_* . We also find that the deviation of the [O III] line width is correlated with the Eddington ratio and sources with Lbol/LEdd~1 have the largest deviations, which are similar to those of radio quiet QSOs/NLS1s (radio jets is absent or weak). A similar result has been obtained for 9 linear radio Seyfert galaxies. On the basis of these results, we suggest that, in addition to the possible jet-gas interactions, accretion activities may also play an important role in shaping the kinematics of the narrow [O III] line in young radio galaxies and linear radio Seyferts.

In this paper, we study the implementation of brane worlds in type II string theory. Starting with the NS/NS sector of type II string, we first compactify the $(D+d_{+} + d_{-})$-dimensional spacetime, and reduce the corresponding action to a D-dimensional effective action, where the topologies of $M_{d_{+}}$ and $M_{d_{-}}$ are arbitrary. We further compactify one of the $(D-1)$ spatial dimensions on an $S^{1}/Z_{2}$ orbifold, and derive the gravitational and matter field equations both in the bulk and on the branes. Then, we investigate two key issues in such a setup: (i) the radion stability and radion mass; and (ii) the localization of gravity, and the corresponding Kaluza-Klein (KK) modes. We show explicitly that the radion is stable and its mass can be in the order of $GeV$. In addition, the gravity is localized on the visible brane, and its spectrum of the gravitational KK towers is discrete and can have a mass gap of $TeV$, too. The high order Yukawa corrections to the 4-dimensional Newtonian potential is exponentially suppressed, and can be negligible. Applying such a setup to cosmology, we obtain explicitly the field equations in the bulk and the generalized Friedmann equations on the branes.

The BH masses and the Eddington ratios for a sample of 65 young radio galaxies [27 GPS and 38 CSS sources] are estimated by various methods. We find that the average BH mass of these young radio galaxies is <log M_bh>=8.3, which is less than that of radio loud QSOs and low redshift radio galaxies(<log M_bh>=9.0). The CSS/GPS sources have relatively high Eddington ratios with average ratio <log L_bol/L_Edd>=-0.56, which are similar to those of narrow line Seyfert 1 galaxies (NLS1s). It suggests that the CSS/GPS sources may not only be in the early stage of radio activities, but also in the early stage of their accretion activities. We find that the young radio galaxies as a class deviate systematically from the M_bh-\sigma_* relation defined by nearby inactive galaxies, when using \sigma_[O III] as a surrogate for the stellar velocity dispersion \sigma_*. There is no significant correlation between the deviation of the [O III] emission line width, $\Delta\sigma$=\sigma_[O III]-\sigma_[pred], and the jet/accretion power, where \sigma_[pred] are calculated from the Tremaine et al. relation using the estimated BH masses. However, we find that the deviation $\Delta\sigma$ in young radio galaxies is well correlated with the Eddington ratio, and this correlation is found to be similar to that of radio quiet AGN (QSOs/NLS1s) where the radio jet is absent or weak. We suggest that the accretion activities may still play an important role in shaping the kinematics of [O III] narrow line in these young radio galaxies.

In the frame of unification schemes for radio-loud active galactic nuclei (AGNs), FR I radio galaxies are believed to be BL Lacertae (BL Lac) objects with the relativistic jet misaligned to our line of sight, and FR II radio galaxies correspond to misaligned radio quasars. The Ledlow-Owen dividing line for FR I/FR II dichotomy in the optical absolute magnitude of host galaxy-radio luminosity (M_R-L_Rad) plane can be translated to the line in the black hole mass-jet power (M_bh-Q_jet) plane by using two empirical relations: Q_jet-L_Rad and M_bh}-M_R. We use a sample of radio quasars and BL Lac objects with measured black hole masses to explore the relation of the jet power with black hole mass, in which the jet power is estimated from the extended radio emission. It is found that the BL Lac objects are clearly separated from radio quasars by the Ledlow & Owen FR I/II dividing line in the M_bh-Q_jet plane. This strongly supports the unification schemes for FR I/BL Lac object and FR II/radio quasar. We find that the Eddington ratios L_bol/L_Edd of BL Lac objects are systematically lower than those of radio quasars in the sample with a rough division at L_bol/L_Edd 0.01, and the distribution of Eddington ratios of BL Lac objects/quasars exhibits a bimodal nature, which imply that the accretion mode of BL Lac objects may be different from that of radio quasars.

We explore the relationship between the hard X-ray photon index $\Gamma$ and the Eddington ratio (\xi=L_{X}(0.5-25 keV)/L_{Edd}) in six XRBs. We find that different XRBs follow different anti-correlations between $\Gamma$ and $\xi$ when $\xi$ is less than a critical value, while they follow the same positive correlation when $\xi$ is larger than the critical value. This anti-correlation and positive correlation are also found in LLAGNs and QSOs respectively, and the anti-correlation and positive correlation of different XRBs roughly converge to the same point ($\log \xi=-2.1, \Gamma=1.5$), which may correspond to the accretion mode transition, since that the anti-correlation and positive correlation are consistent with the prediction of ADAFs and standard disk/corona system respectively. The traditional low/hard state are divided into two parts by the cross point $\log \xi\sim-2.1$, i.e., faint-hard state in the anti-correlation part and bright-hard state in the positive correlation part. The accretion process in the bright-hard state may be still the standard accretion disk as that in the high/soft state, which is consistent with that both the cold disk component and broad Fe K emission line are observed in some bright-hard state of XRBs. The ADAF is only important in the faint-hard state XRBs. Motivated by the similarities of the state transition and timing properties of the ULXs to that of XRBs, we then constrain the BH masses for seven luminous ULXs assuming that their X-ray spectral evolution is similar to that of XRBs. We find that the BH masses of these seven ULXs are around $10^{4}M_sun$, which are typical intermediate mass BHs (IMBHs). Our results are roughly consistent with the BH masses constrained from the model fitting with a multi-color disk and/or the timing properties(e.g., QPO and break frequency).

We investigate the X-ray origin in FRIs using the multi-waveband high resolution data of eight FR I sources, which have very low Eddington ratios. We fit their multi-waveband spectrum using a coupled accretion-jet model. We find that X-ray emission in the source with the highest L_X (~1.8*10^-4 L_Edd) is from the advection-dominated accretion flow (ADAF). Four sources with moderate L_X(~several*10^-6 L_Edd) are complicated. The X-ray emission of one FR I is from the jet, and the other three is from the sum of the jet and ADAF. The X-ray emission in the three least luminous sources (L_X<1.0*10^-6L_Edd) is dominated by the jet. These results roughly support the predictions of Yuan and Cui(2005) where they predict that when the X-ray luminosity of the system is below a critical value, the X-radiation will not be dominated by the emission from the ADAF any longer, but by the jet. We also find that the accretion rates in four sources must be higher than the Bondi rates, which implies that other fuel supply (e.g., stellar winds) inside the Bondi radius should be important.

We calculate the jet power of the classical Blandford-Znajek(BZ) model and hybrid model developed by Meier based on the global solutions of advection dominated accretion flows (ADAFs) surrounding Kerr black holes. We find that the jet power of the hybrid model is larger than that of the pure BZ model. The jet power will dominate over the accretion power, and the objects will enter into "jet-power-dominated advective systems", when the accretion rate is less than a critical value mdot_c=Mdot_c/Mdot_Edd, where 3*10^-4 < mdot_c < 5*10^-3 is a function of black hole spin parameter. The accretion power will be dominant when mdot<mdot_c and the objects will enter into "accretion-power-dominated advective systems." This is roughly consistent with that constrained from the low/hard-state black hole X-ray binaries (e.g., Fender et al.). We calculate the maximal jet power as a function of black hole mass with the hybrid jet formation model, and find it can roughly reproduce the dividing line of the Ledlow-Owen relation for FR I/FR II dichotomy in the jet power-black hole(BH) mass plane (Q_jet-M_BH) if the dimensionless accretion rate mdot~0.01 and BH spin parameter j~0.9-0.99 are adopted. This accretion rate mdot~0.01 is consistent with that of the critical accretion rate for the accretion mode transition of a standard disk to an ADAF constrained from the state transition of X-ray binaries. Our results imply that most FR I galaxies may be in the ADAF accretion mode similar to the low/hard-state XRBs.

Recently, in Gong {\em et al} \cite{GWW07} and Wang and Santos \cite{WS07} it was shown that the effective cosmological constant on each of the two orbifold branes can be easily lowered to its current observational value, by using the large extra dimensions in the framework of both M-Theory and string theory on $S^{1}/Z_{2}$. In this paper, we study the current acceleration of the universe, using the formulas developed in \cite{WS07}. We first construct explicitly time-dependent solution to the 10-dimensional bulk of the Neveu-Schwarz/Neveu-Schwarz sector, compactified on a 5-dimensional torus. Then, we write down the generalized Friedmann equations on each of the two dynamical branes, and fit the models to the 182 gold supernova Ia data and the BAO parameter from SDSS, using both of our MINUIT and Monte-Carlo Markov Chain (MCMC) codes. With the best fitting values of the parameters involved as initial conditions, we integrate the generalized Friedmann equations numerically and find the future evolution of the universe. We find that it depends on the choice of the radion potentials $V_{4}^{(I)} (I = 1, 2)$ of the branes. In particular, when choosing them to be the Goldberger-Wise potentials, $V_{4}^{(I)} = \lambda_{4}^{(I)} (\psi^{2} - {v_{I}}^{2})^{2}$, we find that the current acceleration of the universe driven by the effective cosmological constant is only temporary. Due to the effects of the potentials, the universe will be finally in its decelerating expansion phase again. We also study the proper distance between the two branes, and find that it remains almost constant during the whole future evolution of the universe in all the models considered.

We use the Monte-Carlo Markov Chain method to explore the dark energy property and the cosmic curvature by fitting two popular dark energy parameterizations to the observational data. The new 182 gold supernova Ia data and the ESSENCE data both give good constraint on the DE parameters and the cosmic curvature for the dark energy model $w_0+w_a z/(1+z)$. The cosmic curvature is found to be $|\Omega_k|\la 0.03$. For the dark energy model $w_0+w_a z/(1+z)^2$, the ESSENCE data gives better constraint on the cosmic curvature and we get $|\Omega_k|\leq 0.02$.

Orbifold branes are studied in the framework of the 11-dimensional Horava-Witten heterotic M-Theory. It is found that the effective cosmological constant can be easily lowered to its current observational value by the mechanism of large extra dimensions. The domination of this constant over the evolution of the universe is only temporarily. Due to the interaction of the bulk and the branes, the universe will be in its decelerating expansion phase again in the future, whereby all problems connected with a far future de Sitter universe are resolved.

We investigate the relationship between the hard X-ray photon index $\Gamma$ and the Eddington ratio ($\xi=L_{X}(0.5-25 \rm keV)/L_{Edd}$) in six X-ray binaries (XRBs) with well constrained black hole masses and distances. We find that different XRBs follow different anti-correlations between $\Gamma$ and $\xi$ when $\xi$ is less than a critical value, while $\Gamma$ and $\xi$ generally follow the same positive correlation when $\xi$ is larger than the critical value. The anti-correlation and the positive correlation may suggest that they are in different accretion modes (e.g., radiatively inefficient accretion flow (RIAF) and standard disk). We fit both correlations with the linear least-square method for individual sources, from which the crosspoint of two fitted lines is obtained. Although the anti-correlation varies from source to source, the crosspoints of all sources roughly converge to the same point with small scatter($\log \xi=-2.1\pm0.2, \Gamma=1.5\pm 0.1$), which may correspond to the transition point between RIAF and standard accretion disk. Motivated by the observational evidence for the similarity of the X-ray spectral evolution of ultraluminous X-ray sources (ULXs) to that of XRBs, we then constrain the black hole masses for seven ULXs assuming that their X-ray spectral evolution is similar to that of XRBs. We find that the BH masses of these seven luminous ULXs are around $10^{4}\msun$, which are typical intermediate-mass BHs (IMBHs). Our results are generally consistent with the BH masses constrained from the timing properties (e.g., break frequency) or the model fitting with a multi-color disk.

Although there is mounting observational evidence that the cosmic expansion is undergoing a late-time acceleration, the physical mechanism behind such a phenomenon is yet unknown. In this paper, we investigate a holographic dark energy (HDE) model with interaction between the components of the dark sector in the light of current cosmological observations. We use both the new \emph{gold} sample of 182 type Ia supernovae (SNe Ia) and the 192 SNe Ia ESSENCE data, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and the shift parameter from the three-year Wilkinson Microwave Anisotropy Probe data. In agreement with previous results, we show that these observations suggest a very weak coupling or even a noninteracting HDE. The phantom crossing behavior in the context of these scenarios is also briefly discussed.

The energy conditions give upper bounds on the luminosity distance. We apply these upper bounds to the 192 essence supernova Ia data to show that the Universe had experienced accelerated expansion. This conclusion is drawn directly from the distance modulus-reshift graph. In addition to be a very simple method, this method is also totally independent of any cosmological model. From the degeneracy of the distance modulus at low redshift, we argue that the choice of $w_0$ for probing the property of dark energy is misleading. One explicit example is used to support this argument.

We investigate the X-ray origin in FR Is using the radio, submillimetre, optical, and {\em Chandra} X-ray data of a small sample consisting of eight FR I sources. These sources are very dim, with X-ray luminosities $L_{\rm X}/L_{\rm Edd} \sim 10^{-4}-10^{-8}$ ($L_{\rm X}$ is the X-ray luminosity between 2-10 keV). We try to fit the multiwaveband spectrum using a coupled accretion-jet model. In this model, the accretion flow is described by an advection-dominated accretion flow (ADAF) while in the innermost region of ADAF a fraction of accretion flow is transferred into the vertical direction and forms a jet. We find that X-ray emission in the source with the highest $L_{\rm X}$ ($\sim 1.8 \times 10^{-4}L_{\rm Edd}$) is from the ADAF. The results for the four sources with moderate $L_{\rm X}$ ($\sim$ several $\times 10^{-6}L_{\rm Edd}$) are complicated. Two are mainly from the ADAFs, one from the jet, and the other from the sum of the jet and ADAF. The X-ray emission in the three least luminous sources ($L_{\rm X} \lesssim 1.0\times 10^{-6}L_{\rm Edd}$) is dominated by the jet although for one source it can also be interpreted by the ADAF since the quality of X-ray data is low. We conclude that these results roughly support the predictions of Yuan & Cui (2005) where they predict that when the X-ray luminosity of the system is below a critical value, the X-radiation will not be dominated by the emission from the ADAF any longer, but by the jet. We also investigate the fuel supply in these sources. We find that the accretion rates in four sources among the five in which we can have good constraints to their accretion rates must be higher than the Bondi rates. This implies that other fuel supply, such as the gas released by the stellar population inside the Bondi radius, should be important.

It is argued that there is a linear correlation between star formation rate (SFR) and accretion rate for normal bright active galactic nuclei (AGNs). However, it is still unclear whether this correlation holds for LINERs, of which the accretion rates are relatively lower than those of normal bright AGNs. The radiatively inefficient accretion flows (RIAFs) are believed to be present in these LINERs. In this work, we derive accretion rates for a sample of LINERs from their hard X-ray luminosities based on spectral calculations for RIAFs. We find that LINERs follow the same correlation between star formation rate and accretion rate defined by normal bright AGNs, when reasonable parameters are adopted for RIAFs. It means that the gases feed the black hole and star formation in these low-luminosity LINERs may follow the same way as that in normal bright AGNs, which is roughly consistent with recent numerical simulations on quasar evolution.

We use the observational data in radio, optical and X-ray wavebands, for a sample of active galactic nuclei (AGNs) with measured black hole masses, to explore the origin of radio emission from nearby low-luminosity active galactic nuclei (LLAGNs). We find that the radio luminosities are higher than the maximal luminosities expected from the ADAF model, for most sources in this sample. This implies that the radio emission is dominantly from the jets in these sources. The radio emission from a small fraction of the sources in this sample can be explained by the ADAF model. However, comparing the observed multi-band emission data with the spectra calculated for the ADAF or ADIOS cases, we find that neither ADAF nor ADIOS models can reproduce the observed multi-band emission simultaneously, with reasonable magnetic field strengths, for these radio-weak sources. A variety of other possibilities are discussed, and we suggest that the radio emission is probably dominated by jet emission even in these radio-weak LLAGNs.