In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

We construct a theoretical framework to study Population III (Pop III) star formation in the post-reionization epoch ($z\lesssim 6$) by combining cosmological simulation data with semi-analytical models. We find that due to radiative feedback (i.e. Lyman-Werner and ionizing radiation) massive haloes ($M_{\rm halo}\gtrsim 10^{9}\ \rm M_{\odot}$) are the major ($\gtrsim 90$\%) hosts for potential Pop III star formation at $z\lesssim 6$, where dense pockets of metal-poor gas may survive to form Pop III stars, under inefficient mixing of metals released by supernovae. Metal mixing is the key process that determines not only when Pop III star formation ends, but also the total mass, $M_{\rm PopIII}$, of \textit{active} Pop III stars per host halo, which is a crucial parameter for direct detection and identification of Pop III hosts. Both aspects are still uncertain due to our limited knowledge of metal mixing during structure formation. Current predictions range from early termination at the end of reionization ($z\sim 5$) to continuous Pop III star formation extended to $z=0$ at a non-negligible rate $\sim 10^{-7}\ \rm M_{\odot}\ yr^{-1}\ Mpc^{-3}$, with $M_{\rm PopIII}\sim 10^{3}-10^{6}\ \rm M_{\odot}$. This leads to a broad range of redshift limits for direct detection of Pop III hosts, $z_{\rm PopIII}\sim 0.5-12.5$, with detection rates $\lesssim 0.1-20\ \rm arcmin^{-2}$, for current and future space telescopes (e.g. HST, WFIRST and JWST). Our model also predicts that the majority ($\gtrsim 90$\%) of the cosmic volume is occupied by metal-free gas. Measuring the volume filling fractions of this metal-free phase can constrain metal mixing parameters and Pop III star formation.

We propose a pebble-driven planet formation scenario to form giant planets with high multiplicity and large orbital distances in the early gas disk phase. We perform N-body simulations to investigate the growth and migration of low-mass protoplanets in the disk with inner viscously heated and outer stellar irradiated regions. The key feature of this model is that the giant planet cores grow rapidly by a combination of pebble accretion and planet-planet collisions. This consequently speeds up their gas accretion. Because of efficient growth, the planet transitions from rapid type I migration to slow type II migration early, reducing the inward migration substantially. Multiple giant planets can sequentially form in this way with increasing semimajor axes. Both mass growth and orbital retention are more pronounced when a large number of protoplanets are taken into account compared to the case of single planet growth. Eventually, a few numbers of giant planets form with orbital distances of a few to a few tens of AUs within $1.5{-}3$ Myr after the birth of the protoplanets. The resulting simulated planet populations could be linked to the substructures exhibited in disk observations as well as large orbital distance exoplanets observed in radial velocity and microlensing surveys.

Observationally, an anti-correlation between the X-ray photon index $\Gamma$ (obtained by fitting the X-ray spectrum between 0.5 and 10 keV with a single power law) and the X-ray luminosity $L_{\rm 0.5-10keV}$, is found in neutron star low-mass X-ray binaries (NS-LMXBs) in the range of $L_{\rm 0.5-10keV}\sim 10^{34}-10^{36}\ \rm erg\ s^{-1}$. In this paper, we explain the observed anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$ within the framework of the self-similar solution of the advection-dominated accretion flow (ADAF) around a weakly magnetized NS. The ADAF model intrinsically predicts an anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$. In the ADAF model, there is a key parameter, $f_{\rm th}$, which describes the fraction of the ADAF energy released at the surface of the NS as thermal emission to be scattered in the ADAF. We test the effect of $f_{\rm th}$ on the anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$. It is found that the value of $f_{\rm th}$ can significantly affect the anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$. Specifically, the anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$ becomes flatter with decreasing $f_{\rm th}$ as taking $f_{\rm th}=0.1, 0.03, 0.01, 0.005$, $0.003$ and $0$ respectively. By comparing with a sample of non-pulsating NS-LMXBs with well measured $\Gamma$ and $L_{\rm 0.5-10keV}$, we find that indeed only a small value of $0.003\lesssim f_{\rm th}\lesssim 0.1$ is needed to match the observed anti-correlation between $\Gamma$ and $L_{\rm 0.5-10keV}$. Finally, we argue that the small value of $f_{\rm th}\lesssim 0.1$ derived in this paper further confirms our previous conclusion that the radiative efficiency of NSs with an ADAF accretion may not be as high as $\epsilon \sim {\dot M GM\over R_{*}}/{\dot M c^2}\sim 0.2$.

We present a novel deep learning model, intelligent quasar continuum neural network (iQNet), which predicts the intrinsic continuum of any quasar in the rest-frame wavelength range 1020 Angstroms $\leq \lambda \leq$ 1600 Angstroms. We train this network using quasar spectra at low redshift ($z \sim 0.2$) from the Hubble Spectroscopic Legacy Archive, and apply it to predict quasar continua from different astronomical surveys. We introduce a standardization process to the data, reducing the absolute fractional flux error (AFFE) of the predicted continua approximately by half. We use principal component analysis and Gaussian mixture model to classify the HSLA quasar spectra into four classes and use them to synthesize mock quasar spectra create a training data set for iQNet. iQNet achieves a median AFFE of 1.31% on the training quasar spectra, approximately ten times better than traditional PCA-based prediction methods, and 4.17% on the testing quasar spectra. We apply iQNet to predict the continua of $\sim$ 3200 quasar spectra at higher redshift ($2< z \leq 5$) and measure the redshift evolution of mean transmitted flux ($< F >$) in the Ly-$\alpha$ forest region. We measure a gradual evolution of $< F >$ with redshift, which we characterize as a power-law evolution. These estimates are broadly consistent with literature, but provide a more accurate measurement as we measure the quasar continuum with minimum contamination from the Ly-$\alpha$ forest. This work proves that the iQNet model can predict the quasar continuum with high accuracy and shows the viability of such methods for quasar continuum prediction.

A long-standing question in active galactic nucleus (AGN) research is how the corona is heated up to produce X-ray radiation much stronger than that arising from the viscous heating within the corona. In this paper, we carry out detailed investigations of magnetic-reconnection heating to the corona, specifically, studying how the disc and corona are self-consistently coupled with the magnetic field, and how the emergent spectra depend on the fundamental parameters of AGN. It is shown that diverse spectral shapes and luminosities over a broad bandpass from optical to X-ray can be produced from the coupled disc and corona within a limited range of the black hole mass, accretion rate and magnetic field strength. The relative strength of X-ray emission with respect to optical/ultraviolet (UV) depends on the strength of the magnetic field in the disc, which, together with accretion rate, determines the fraction of accretion energy transported and released in the corona. This refined disc-corona model is then applied to reproduce the broad-band spectral energy distributions (SEDs) of a sample of 20 bright local AGNs observed simultaneously in X-ray and optical/UV. We find that, in general, the overall observed broad-band SEDs can be reasonably reproduced, except for rather hard X-ray spectral shapes in some objects. The radiation pressure-dominant region, as previously predicted for the standard accretion disc in AGN, disappears for strong X-ray sources, revealing that AGN accretion discs are indeed commonly stable as observed. Our study suggests the disc-corona coupling model involving magnetic fields to be a promising approach for understanding the broad-band spectra of bright AGNs.

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.

We study the dynamics of merging compact binaries near a rotating supermassive black hole (SMBH) in a hierarchical triple configuration. We include various general relativistic effects that couple the inner orbit, the outer orbit and the spin of the SMBH. During the orbital decay due to gravitational radiation, the inner binary can encounter an "apsidal precession resonance" and experience eccentricity excitation. This resonance occurs when the apsidal precession rate of the inner binary matches that of the outer binary, with the precessions driven by both Newtonian interactions and various post-Newtonian effects. The eccentricity excitation requires the outer orbit to have a finite eccentricity, and is most effective for triples with small mutual inclinations, in contrast to the well-studied Lidov-Kozai effect. The resonance and the associated eccentricity growth may occur while the binary emits gravitational waves in the low-frequency band, and may be detectable by future space-based gravitational wave detectors.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use cosmological hydrodynamic simulations to study the gravitational wave (GW) signals from high-redshift binary black holes (BBHs) formed by dynamical capture (ex-situ formation channel). We in particular focus on BHs originating from the first generation of massive, metal-poor, so-called Population~III (Pop~III) stars. An alternative (in-situ) formation pathway arises in Pop~III binary stars, whose GW signature has been intensively studied. In our optimistic model, we predict a local GW event rate density for ex-situ BBHs (formed at z>4) of ~0.04 /(yr Gpc^3). This is comparable to or even higher than the conservative predictions of the rate density for in-situ BBHs ~0.01-0.1 /(yr Gpc^3), indicating that the ex-situ formation channel may be as important as the in-situ one for producing GW events. We also evaluate the detectability of our simulated GW events for selected planned GW instruments, such as the Einstein Telescope (ET). For instance, we find the all-sky detection rate with signal-to-noise ratios above 10 to be <100 per year for the xylophone configuration of ET. However, our results are highly sensitive to the sub-grid models for BBH identification and evolution, such that the GW event efficiency (rate) is reduced by a factor of 4 (20) in the pessimistic case. The ex-situ channel of Pop III BBHs deserves further investigation with better modeling of the environments around Pop III-seeded BHs.

We demonstrate a new approach of indirectly constraining both early star and structure formation via mature galaxy clusters at cosmic noon (z~2), using the cluster XLSSC 122 as an example. With the standard Press Schechter formalism, we infer a rapid evolution of the star formation efficiency (the ratio of stellar to halo mass) from 10^-4 to 0.01 during z~20-13, based on the age distribution of stars in post-starburst galaxies of XLSSC 122, measured by HST photometry assuming no dust extinction. Here, we consider all low-mass haloes, including minihaloes, that host the first stars and galaxies (5*10^5 Msun < M_halo < 10^10 Msun). We also place new constraints on fuzzy dark matter models of m_a < 5*10^-21 eV/c^2 for the ultra-light boson mass, from the abundance of galaxies with star formation at z > 13 in XLSSC 122. Our exploratory results are consistent with existing constraints. More comprehensive results will be obtained if our approach is extended to a large sample of clusters or field post-starburst galaxies at cosmic noon, with improved modelling of halo and stellar populations.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

Context: Cygnus X-1 is a black hole X-ray binary system in which the black hole captures and accretes gas from the strong stellar wind emitted by its supergiant O9.7 companion star. The irradiation of the supergiant star essentially determines the flow properties of the stellar wind and the X-ray luminosity from the system. The results of three-dimensional hydrodynamical simulations of wind-fed X-ray binary systems reported in recent work reveal that the ionizing feedback of the X-ray irradiation leads to the existence of two stable states with either a soft or a hard spectrum. Aims: We discuss the observed radiation of Cygnus X-1 in the soft and hard state in the context of mass flow in the corona and disk, as predicted by the recent application of a condensation model. Methods: The rates of gas condensation from the corona to the disk for Cygnus X-1 are determined, and the spectra of the hard and soft radiation are computed. The theoretical results are compared with the MAXI observations of Cygnus X-1 from 2009 to 2018. In particular, we evaluate the hardness-intensity diagrams (HIDs) for its ten episodes of soft and hard states which show that Cygnus X-1 is distinct in its spectral changes as compared to those found in the HIDs of low-mass X-ray binaries. Results: The theoretically derived values of photon counts and hardness are in approximate agreement with the observed data in the HID. However, the scatter in the diagram is not reproduced. Improved agreement could result from variations in the viscosity associated with clumping in the stellar wind and corresponding changes of the magnetic fields in the disk. The observed dipping events in the hard state may also contribute to the scatter and to a harder spectrum than predicted by the model.

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$).

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

Nearly-axisymmetric gaps and rings are commonly observed in protoplanetary discs. The leading theory regarding the origin of these patterns is that they are due to dust trapping at the edges of gas gaps induced by the gravitational torques from embedded planets. If the concentration of solids at the gap edges becomes high enough, it could potentially result in planetesimal formation by the streaming instability. We test this hypothesis by performing global 1-D simulations of dust evolution and planetesimal formation in a protoplanetary disc that is perturbed by multiple planets. We explore different combinations of particle sizes, disc parameters, and planetary masses, and find that planetesimals form in all these cases. We also compare the spatial distribution of pebbles from our simulations with protoplanetary disc observations. Planets larger than one pebble isolation mass catch drifting pebbles efficiently at the edge of their gas gaps, and depending on the efficiency of planetesimal formation at the gap edges, the protoplanetary disc transforms within a few 100,000 years to either a transition disc with a large inner hole devoid of dust or to a disc with narrow bright rings. For simulations with planetary masses lower than the pebble isolation mass, the outcome is a disc with a series of weak ring patterns but no strong depletion between the rings. Lowering the pebble size artificially to 100 micrometer-sized "silt", we find that regions between planets get depleted of their pebble mass on a longer time-scale of up to 0.5 million years. These simulations also produce fewer planetesimals than in the nominal model with millimeter-sized particles and always have at least two rings of pebbles still visible after 1 Myr.

The Green Bank Telescope (GBT) Diffuse Ionized Gas Survey (GDIGS) is a fully-sampled radio recombination line (RRL) survey of the inner Galaxy at C-band (4-8 GHz). We average together ~15 Hn$\alpha$ RRLs within the receiver bandpass to improve the spectral signal-to-noise ratio. The average beam size for the RRL observations at these frequencies is ~2'. We grid these data to have spatial and velocity spacings of 30" and 0.5 km/s, respectively. Here we discuss the first RRL data from GDIGS: a six square-degree-area surrounding the Galactic HII region complex W43. We attempt to create a map devoid of emission from discrete HII regions and detect RRL emission from the diffuse ionized gas (DIG) across nearly the entire mapped area. We estimate the intensity of the DIG emission by a simple empirical model, taking only the HII region locations, angular sizes, and RRL intensities into account. The DIG emission is predominantly found at two distinct velocities: ~40 km/s and ~100 km/s. While the 100 km/s component is associated with W43 at a distance of ~6 kpc, the origin of the 40 km/s component is less clear. Since the distribution of the 40 km/s emission cannot be adequately explained by ionizing sources at the same velocity, we hypothesize that the plasma at the two velocity components is interacting, placing the 40 km/s DIG at a similar distance as the 100 km/s emission. We find a correlation between dust temperature and integrated RRL intensity, suggesting that the same radiation field that heats the dust also maintains the ionization of the DIG.

Observationally, the X-ray spectrum ($0.5-10$ keV) of low-level accreting neutron stars (NSs) can generally be well fitted by the model with two components, i.e, a thermal soft X-ray component plus a power-law component. Meanwhile, the fractional contribution of the power-law luminosity $\eta$ ($\eta\equiv L^{\rm power\ law}_{\rm 0.5-10\rm keV}/L_{\rm 0.5-10\rm keV}$) varies with the X-ray luminosity $L_{\rm 0.5-10\rm keV}$. In this paper, we systematically investigate the origin of such X-ray emission within the framework of the advection-dominated accretion flow (ADAF) around a weakly magnetized NS, in which the thermal soft X-ray component arises from the surface of the NS and the power-law component arises from the ADAF itself. We test the effects of the viscosity parameter $\alpha$ in the ADAF and thermalized parameter $f_{\rm th}$ (describing the fraction of the ADAF energy released at the surface of the NS as thermal emission) on the relation of $\eta$ versus $L_{\rm 0.5-10\rm keV}$. It is found that $\eta$ is nearly a constant ($\sim$ zero) with $L_{\rm 0.5-10\rm keV}$ for different $\alpha$ with $f_{\rm th}=1$, which is inconsistent with observations. Meanwhile, it is found that a change of $f_{\rm th}$ can significantly change the relation of $\eta$ versus $L_{\rm 0.5-10\rm keV}$. By comparing with a sample of non-pulsating NS-LMXBs probably dominated by low-level accretion onto NSs, it is found that a small value of $f_{\rm th} \lesssim 0.1$ is needed to match the observed range of $\eta \gtrsim 10\%$ in the diagram of $\eta$ versus $L_{\rm 0.5-10\rm keV}$. Finally, we argue that the small value of $f_{\rm th} \lesssim 0.1$ implies that the radiative efficiency of NSs with an ADAF accretion may not be as high as the predicted result previously of $\epsilon \sim {\dot M GM\over R_{*}}/{\dot M c^2}\sim 0.2$ despite the existence of the hard surface.

Impact craters, as "lunar fossils", are the most dominant lunar surface features and occupy most of the Moon's surface. Their formation and evolution record the history of the Solar System. Sixty years of triumphs in the lunar exploration projects accumulated a large amount of lunar data. Currently, there are 9137 existing recognized craters. However, only 1675 of them have been determined age, which is obviously not satisfactory to reveal the evolution of the Moon. Identifying craters is a challenging task due to their enormous difference in size, large variations in shape and vast presence. Furthermore, estimating the age of craters is extraordinarily difficult due to their complex and different morphologies. Here, in order to effectively identify craters and estimate their age, we convert the crater identification problem into a target detection task and crater age estimation into a taxonomy structure. From an initial small number of available craters, we progressively identify craters and estimate their age from Chang'E data by transfer learning (TL) using deep neural networks. For comprehensive identification of multi-scale craters, a two-stage craters detection approach is developed. Thus 117240 unrecognized lunar craters that range in diameter from 532 km to 1 km are identified. Then, a two-stage classification approach is developed to estimate the age of craters by simultaneously extracting their morphological features and stratigraphic information. The age of 79243 craters larger than 3 km in diameter is estimated. These identified and aged craters throughout the mid and low-latitude regions of the Moon are crucial for reconstructing the dynamic evolution process of the Solar System.

The influences and consequences of toroidal magnetic fields in shaping the visible lobes of pre planetary nebulae ("prePNe") are explored in this, the last of a series of papers of parameter studies of prePN evolution. To probe these influences we start with the steady, diverging, and field-free wind model of our previous papers and add weak to moderate toroidal fields to the winds in order to generate arrays of outcomes after 500 y, after which the structures grow almost homologously. As expected, toroidal fields in the stellar winds invariably form very thin and dense axial features whose structure is best described as a thin cold jet with an ultra-dense and neutral leading knot, or "bullet", at its tip. The speed of the leading knot depends only on the density contrast (the ratio of injected to ambient gas densities at the nozzle) and wind injection speed, but not on the field strength or opening angle. The lobes formed by the ram pressure of the winds take a variety of forms and sizes that depend primarily on the geometric structure of the injected gas and the density contrast. About 20% of the HST images of prePNe show the unique signatures of shaping by toroidal fields. Pairs of low-ionization knots seen along the major axis of fully ionized PNe, often called "FLIERs" are easily explained as the very dense, cold , and neutral remnants of magnetically formed knots.

Cold atomic hydrogen clouds are the precursors of molecular clouds. Due to self-absorption, the opacity of cold atomic hydrogen may be high, and this gas may constitute an important mass component of the interstellar medium (ISM). Atomic hydrogen gas can be cooled to temperatures much lower than found in the cold neutral medium (CNM) through collisions with molecular hydrogen. In this paper, we search for HI Narrow Self-Absorption (HINSA) features in the Large Magellanic Cloud (LMC) as an indicator of such cold HI clouds, and use the results to quantify atomic masses and atomic-to-molecular gas ratio. Our search for HINSA features was conducted towards molecular clouds in the LMC using the ATCA+Parkes HI survey and the MAGMA CO survey. HINSA features are prevalent in the surveyed sightlines. This is the first detection of HINSA in an external galaxy. The HINSA-HI/$\rm{H}_{2}$ ratio in the LMC varies from 0.5\e{-3} to 3.4\e{-3} (68\% interval), with a mean value of $(1.31 \pm 0.03)$\e{-3}, after correcting for the effect of foreground HI gas. This is similar to the Milky Way value and indicates that similar fractions of cold gas exist in the LMC and the Milky Way, despite their differing metallicities, dust content and radiation fields. The low ratio also confirms that, as with the Milky Way, the formation timescale of molecular clouds is short. The ratio shows no radial gradient, unlike the case for stellar metallicity. No correlation is found between our results and those from previous HI absorption studies of the LMC.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

IRAS F11119+3257 is a quasar-dominated Ultra-Luminous InfraRed Galaxy, with a partially obscured narrow-line seyfert 1 nucleus. In this paper, we present the NIR spectroscopy of F11119+3257, in which we find unusual Paschen emission lines, and metastable He I* $\lambda$10830 absorption associated with the previously reported atomic sodium and molecular OH mini-BAL (Broad Absorption Line) outflow. Photo-ionization diagnosis confirms previous findings that the outflows are at kilo-parsec scales. Such large-scale outflows should produce emission lines. We indeed find that high-ionization emission lines ([O III], [Ne III], and [Ne V]) are dominated by blueshifted components at similar speeds to the mini-BALs. The blueshifted components are also detected in some low-ionization emission lines, such as [O II] $\lambda$3727 and some Balmer lines (H$\alpha$, H$\beta$, and H$\gamma$), even though their cores are dominated by narrow ($FWHM_{\rm NEL} = 570\pm40$km s$^{-1}$) or broad components at the systemic redshift of $z=0.18966\pm0.00006$. The mass flow rate (230-730$~M_\odot \rm yr^{-1}$) and the kinetic luminosity ($\dot{E}_k \sim 10^{43.6-44.8} $erg s$^{-1}$) are then inferred jointly from the blueshifted emission and absorption lines. In the NIR spectrum of F11119+3257, we also find that the Paschen emission lines are unique, in which a very narrow ($FWHM=260\pm20~$km s$^{-1}$) component is shown in only Pa$\alpha$. This narrow component most probably comes from heavily obscured star formation. Based on the Pa$\alpha$ and Pa$\beta$ emissions, we obtain an extinction at the $H$ band, $A_H~>~2.1$ (or a reddenning of $E_{B-V}~>~$3.7), and a star formation rate of $SFR~>~130\rm M_\odot yr^{-1}$ that resembles the estimates inferred from the FIR emissions ($SFR_{\rm FIR} = 190\pm90$ $M_\odot$ yr$^{-1}$).

We extract the multiple-epoch Balmer-line profiles of the heavily obscured quasar 3C 445 from the spectral curves in the literature, and analyze the emission-line profiles of the H$\alpha$ and H$\beta$ lines and the profile variability in the H$\alpha$ line in the large time interval of more than three decades. The profile comparison between the H$\alpha$ and H$\beta$ lines shows that both Balmer lines share the profile with the same form, while the blue system of the H$\beta$ line is seriously weaker than that of the H$\alpha$ line. Moreover, the blue system of the H$\alpha$ line suddenly disappeared completely and then did not appear again, however the other two components did not exhibit significant variation in the velocity or the amplitude. These findings suggest that the blue system of 3C 445, as with SDSS J153636.22+044127.0 and its analogs, is probable the result of the shock-heated outflowing gases. The observation angle of almost edge-on which the previous studies suggested can easily produce the high-speed and high-temperature shock in the collision between the massive outflow and the inner surface of the dusty torus.

Quasars are high-luminosity active galactic nuclei believed to be powered by accretion of interstellar matter onto a super-massive black hole (SMBH) therein. Most of the observed energy is released in an accretion disk of inspiralling gas surrounding the SMBH. An enormous amount of fueling material is expected to be transported inwards. However, basic questions remain unanswered as to whether and how the accretion disks are supplied with external gas, since no disk-feeding inflow has hitherto been observed clearly. Here we report the discovery of highly redshifted broad absorption lines arising from neutral hydrogen and helium atoms in a small sample of quasars. Their absorption troughs show a broad range of Doppler velocities from zero extending continuously inward up to as high as $\sim 5,000$ km s$^{-1}$. We thus see through streams of cold gas moving with a radially inward velocity component that spans an immense gradient---a result of gravitational acceleration by the central SMBH. Extensive photo-ionization modeling for the archetypical object SDSS J103516.20+142200.6 indicates the inflowing gas to be dense, thick and moderately ionized, with a characteristic distance to the SMBH of $\sim 1,000$ gravitational radii, possibly overlapping or close to the outer accretion disk. Our results present the first compelling evidence for the long-sought inflow directly feeding quasars' accretion disks with external materials, likely originating from the dusty torus at a parsec scale. Our approach provides a new tool to probe the bulk of the so far elusive fueling inflows in quasars. Their studies may help address some of the fundamental questions concerning accretion physics, the onset and sustainment of quasar activity, and the SMBH growth at centers of most galaxies.

We study the general relativitic (GR) effects induced by a spinning supermassive black hole on the orbital and spin evolution of a merging black hole binary (BHB) in a hierarchical triple system. A sufficiently inclined outer orbit can excite Lidov-Kozai eccentricity oscillations in the BHB and induce its merger. These GR effects generate extra precessions on the BHB orbits and spins, significantly increasing the inclination window for mergers and producing a wide range of spin orientations when the BHB enters LIGO band. This "GR-enhanced" channel may play an important role in BHB mergers.

We develop a pebble-driven model to study the formation and evolution of planets around stars in the mass range of 0.08 and 1 solar mass. The growth and migration of a large number of individual protoplanetary embryos are simulated in a population synthesis manner. We test two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplanetary disks. Two types of disks with different turbulent viscous parameters alpha of 1e-3 and 1e-4 are investigated, to shed light on the role of outward migration of protoplanets. The forming planets are compared with the observed exoplanets in terms of masses, semimajor axes, metallicities, and water contents. We find that gas giant planets are likely to form when the characteristic disk sizes are larger, the disk accretion rates are higher, the disks are more metal-rich and/or their stellar hosts are more massive. Our model shows that 1) the characteristic mass of super-Earth is set by the pebble isolation mass. Super-Earth masses increase linearly with the mass of its stellar host, corresponding to one Earth mass around a late M-dwarf star and 20 Earth masses around a solar-mass star. 2) The low-mass planets up to 20 Earth masses can form around stars with a wide range of metallicities, while massive gas giant planets are preferred to grow around metal-rich stars. 3) Super-Earth planets that are mainly composed of silicates, with relatively low water fractions can form from protoplanetary embryos at the water ice line in weakly turbulent disks where outward migration is suppressed. However, if the embryos are formed over a wide range of radial distances, the super-Earths would end up having a distinctive, bimodal composition in water mass. Our model succeeds in quantitatively reproducing several important observed properties of exoplanets and correlations with their stellar hosts.

We present a detailed study of the gamma-ray emission from the direction of the star-forming region W30 based on a decade of the Fermi Large Area Telescope data in the 0.3-300 GeV photon energy range. The morphological and spectral analyses allow us to resolve the gamma-ray emission into two extended structures from different origins. One of them mostly overlaps with the supernova remnant (SNR) G8.7-0.1 and has a soft spectrum that resembles with the spectra of other middle-aged SNRs interacting with molecular clouds. The other shows remarkable spatial and spectral consistency with the TeV emission from HESS J1804-216, and its spectrum could be naturally explained by inverse Compton scattering of electrons like a number of TeV gamma-ray emitting pulsar wind nebulae. Thus we attribute this source to the nebula around the pulsar PSR J1803-2137.

We present the first results from the Small Magellanic Cloud portion of a new Australia Telescope Compact Array (ATCA) HI absorption survey of both of the Magellanic Clouds, comprising over 800 hours of observations. Our new HI absorption line data allow us to measure the temperature and fraction of cold neutral gas in a low metallicity environment. We observed 22 separate fields, targeting a total of 55 continuum sources against 37 of which we detected HI absorption; from this we measure a column density weighted mean average spin temperature of $<T_{s}>=150$ K. Splitting the spectra into individual absorption line features, we estimate the temperatures of different gas components and find an average cold gas temperature of $\sim{30}$ K for this sample, lower than the average of $\sim{40}$ K in the Milky Way. The HI appears to be evenly distributed throughout the SMC and we detect absorption in $67\%$ of the lines of sight in our sample, including some outside the main body of the galaxy ($N_{\text{HI}}>2\times{10^{21}}$ cm$^{-2}$). The optical depth and temperature of the cold neutral atomic gas shows no strong trend with location spatially or in velocity. Despite the low metallicity environment, we find an average cold gas fraction of $\sim{20\%}$, not dissimilar from that of the Milky Way.

We study the dynamical signatures of black hole (BH) and neutron star (NS) binary mergers via Lidov-Kozai oscillations induced by tertiary companions in hierarchical triple systems. For each type of binaries (BH-BH and BH-NS), we explore a wide range of binary/triple parameters that lead to binary mergers, and determine the distributions of eccentricity ($e_\mathrm{m}$) and spin-orbit misalignment angle ($\theta_\rm {sl}^\rm f$) when the binary enters the LIGO/VIRGO band. We use the double-averaged (over both orbits) and single-averaged (over the inner orbit) secular equations, as well as N-body integration, to evolve systems with different hierarchy levels, including the leading-order post-Newtonian effect, spin-orbit coupling and gravitational radiation. We find that for merging BH-BH binaries with comparable masses, about $7\%$ have $e_\mathrm{m}>0.1$ and $0.7\%$ have $e_{\rm m}>0.9$. The majority of the mergers have significant eccentricities in the LISA band. The BH spin evolution and $\theta_\rm {sl}^\rm f$ are correlated with the orbital evolution and $e_{\rm m}$. Mergers with $e_{\rm m} \lesssim 10^{-3}$ have a distribution of $\theta_\rm {sl}^\rm f$ that peaks around $90^\circ$ (and thus favoring a projected binary spin parameter $\chi_{\rm eff}\sim 0$), while mergers with larger $e_{\rm m}$ have a more isotropic spin-orbit misalignments. For typical BH-NS binaries, strong octuple effects lead to more mergers with non-negligible $e_\mathrm{m}$ (with $\sim 18\%$ having $e_\mathrm{m}>0.1$ and $2.5\%$ having $e_\mathrm{m}>0.9$), and the final BH spin axis tends to be randomly orientated. Measurements or constraints on eccentric mergers and $\theta_\rm {sl}^\rm f$ from LIGO/VIRGO and LISA would provide useful diagnostics on the dynamical formation of merging BH or NS binaries in triples.

Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc.Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.

Fast radio bursts (FRBs) are highly dispersed and probably extragalactic radio flashes with millisecond-duration. Recently, the Canadian Hydrogen Intensity Mapping Experiment (using the CHIME/FRB instrument) has reported detections of 13 FRBs during a pre-commissioning phase. It is more exciting that one of the 13 FRBs is a second source of repeaters which suggests that CHIME/FRB and other wide-field sensitive radio telescopes will find a substantial population of repeating FRBs. We have proposed strongly lensed repeating FRBs as a precision cosmological probe, e.g. constraining the Hubble constant and model-independently estimating the cosmic curvature. Here, we study complementary constraints on the equation of state of dark energy from strongly lensed FRBs to currently available popular probes. It is found that, in the framework of Chevalier-Polarski-Linder parametrization, adding time delay measurement of 30 strongly lensed FRB systems to cosmic microwave background radiation and type Ia supernovae can improve the dark energy figure of merit by a factor 2. In the precision cosmology era, this improvement is of great significance for studying the nature of dark energy.

We report on simulations of laboratory experiments in which magnetized supersonic flows are driven through a wire mesh. The goal of the study was to investigate the ability of such a configuration to generate supersonic, MHD turbulence. We first report on the morphological structures that develop in both magnetized and non-magnetized cases. We then analyze the flow using a variety of statistical measures, including power spectra and probability distribution functions of the density. Using these results we estimate the sonic mach number in the flows downstream of the wire mesh. We find the initially hypersonic (M=20) planar shock through the wire mesh does lead to downstream turbulent conditions. However, in both magnetized and non-magnetized cases, the resultant turbulence was marginally supersonic to transonic (M~1), and highly anisotropic in structure.

Recently, seven Earth-sized planets were discovered around the M-dwarf star TRAPPIST-1. Thanks to transit-timing variations, the masses and therefore the bulk densities of the planets have been constrained, suggesting that all TRAPPIST-1 planets are consistent with water mass fractions on the order of 10%. These water fractions, as well as the similar planet masses within the system, constitute strong constraints on the origins of the TRAPPIST-1 system. In a previous work, we outlined a pebble-driven formation scenario. In this paper we investigate this formation scenario in more detail. We used a Lagrangian smooth-particle method to model the growth and drift of pebbles and the conversion of pebbles to planetesimals through the streaming instability. We used the N-body code \texttt{MERCURY} to follow the composition of planetesimals as they grow into protoplanets by merging and accreting pebbles. This code is adapted to account for pebble accretion, type-I migration, and gas drag. In this way, we modelled the entire planet formation process (pertaining to planet masses and compositions, not dynamical configuration). We find that planetesimals form in a single, early phase of streaming instability. The initially narrow annulus of planetesimals outside the snowline quickly broadens due to scattering. Our simulation results confirm that this formation pathway indeed leads to similarly-sized planets and is highly efficient in turning pebbles into planets ($\sim$50% solids-to-planets conversion efficiency). [...] The water content of planets resulting from our simulations is on the order of 10%, and our results predict a `V-shaped' trend in the planet water fraction with orbital distance: from relatively high (innermost planets) to relatively low (intermediate planets) to relatively high (outermost planets).

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We derive new constraints on the non-gravitational baryon-dark-matter scattering (BDMS) by evaluating the mass thresholds of dark matter (DM) haloes in which primordial gas can cool efficiently to form Population III (Pop III) stars, based on the timing of the observed 21-cm absorption signal. We focus on the BDMS model with interaction cross-section $\sigma=\sigma_{1}[v/(1\ \mathrm{km\ s^{-1}})]^{-4}$, where $v$ is the relative velocity of the encounter. Our results rule out the region in parameter space with $\sigma_{1}\gtrsim 10^{-19}\ \mathrm{cm^{2}}$ and DM particle mass $m_{\chi}c^{2}\lesssim 3\times 10^{-2}\ \mathrm{GeV}$, where the cosmic number density of Pop III hosts at redshift $z\sim 20$ is at least three orders of magnitude smaller than in the standard Lambda cold DM ($\Lambda$CDM) case. In these BDMS models, the formation of Pop III stars is significantly suppressed for $z\gtrsim 20$, inconsistent with the timing of the observed global 21-cm absorption signal. For the fiducial BDMS model with $m_{\chi}c^{2}=0.3$ GeV and $\sigma_{1}=8\times 10^{-20}\ \mathrm{cm^{2}}$, capable of accommodating the measured absorption depth, the number density of Pop III hosts is reduced by a factor of $3-10$ at $z\sim 15-20$, when the 21-cm signal is imprinted, compared with the $\Lambda$CDM model. The confluence of future detailed cosmological simulations with improved 21-cm observations promises to probe the particle-physics nature of DM at the small-scale frontier of early structure formation.

The broadband spectral energy distribution (SED) of Active Galactic Nuclei (AGN) is investigated for a well-selected sample composed of $23$ Seyfert 1 galaxies observed simultaneously in the optical/UV and X-ray bands with the Neil Gehrels {\it Swift} Observatory. The optical to UV continuum spectra are modeled, for the first time, with emission from an accretion disk with a generalized radial temperature profile, in order to account for the intrinsic spectra which are found to be generally redder than the model prediction of the standard Shakura-Sunyaev disk (SSD) ($F_\nu\propto\nu^{+1/3}$). The power-law indices of the radial temperature profile ($T_{\rm eff}(R)\propto R^{-p}$, $R$ is the radius of the accretion disk) are inferred to be $p=0.5$ -- $0.75$ (a median of $0.63$), deviating from the canonical $p=0.75$ for the SSD model as widely adopted in previous studies. A marginal correlation of a flatter radial temperature profile (a smaller $p$ value) with increasing the Eddington ratio is suggested. Such a model produces generally a lower peak of accretion disk emission and thus a smaller bolometric luminosity in some of the AGN, particularly those with high Eddington ratios, than that based on the SSD model by a factor of several. The broadband SED, the bolometric correction factors and their dependence on some of the AGN parameters are re-visited. We suggest that such non-standard SSD disks may operate in AGN and are at least partly responsible for the reddened optical/UV spectra as observed. One possible explanation for these flattened temperature profiles is the mass loss process in form of disk winds/outflows.

We analyze 10 years of Fermi-LAT data towards the SS433/W50 region. With the latest source catalog and diffuse background models, the gamma-ray excess from SS433/W50 is detected with a significance of 6{\sigma} in the photon energy range of 500 MeV - 10 GeV. Our analysis indicates that an extended flat disk morphology is preferred over a point-source description, suggesting that the GeV emission region is much larger than that of the TeV emission detected by HAWC. The size of the GeV emission is instead consistent with the extent of the radio nebula W50, a supernova remnant being distorted by the jets, so we suggest that the GeV emission may originate from this supernova remnant. The spectral result of the GeV emission is also consistent with an supernova remnant origin. We also derive the GeV flux upper limits on the TeV emission region, which put moderate constrains on the leptonic models to explain the multiwavelength data.

Observationally, for neutron star low-mass X-ray binaries, so far, the correlation between the radio luminosity $L_{\rm R}$ and the X-ray luminosity $L_{\rm X}$, i.e., $L_{\rm R}\propto L_{\rm X}^{\beta}$, has been reasonably well-established only in three sources 4U 1728-34, Aql X-1 and EXO 1745-248 in their hard state. The slope $\beta$ of the radio/X-ray correlation of the three sources is different, i.e., $\beta \sim 1.4$ for 4U 1728-34, $\beta \sim 0.4$ for Aql X-1, and $\beta \sim 1.6$ for EXO 1745-248. In this paper, for the first time we explain the different radio/X-ray correlation of 4U 1728-34, Aql X-1 and EXO 1745-248 with the coupled advection-dominated accretion (ADAF)-jet model respectively. We calculate the emergent spectrum of the ADAF-jet model for $L_{\rm X}$ and $L_{\rm R}$ at different $\dot m$ ($\dot m=\dot M/\dot M_{\rm Edd}$), adjusting $\eta$ ($\eta \equiv \dot M_{\rm jet}/\dot M$, describing the fraction of the accreted matter in the ADAF transfered vertically forming the jet) to fit the observed radio/X-ray correlations. Then we derive a fitting formula of $\eta$ as a function of $\dot m$ for 4U 1728-34, Aql X-1 and EXO 1745-248 respectively. If the relation between $\eta$ and $\dot m$ can be extrapolated down to a lower value of $\dot m$, we find that in a wide range of $\dot m$, the value of $\eta$ in Aql X-1 is greater than that of in 4U 1728-34 and EXO 1745-248, implying that Aql X-1 may have a relatively stronger large-scale magnetic field, which is supported by the discovery of the coherent millisecond X-ray pulsation in Aql X-1.

Within standard $\Lambda$CDM cosmology, Population III (Pop III) star formation in minihalos of mass $M_\mathrm{halo}\gtrsim 5\times10^5$ M$_\odot$ provides the first stellar sources of Lyman$\alpha$ (Ly$\alpha$) photons. The Experiment to Detect the Global Epoch of Reionization Signature (EDGES) has measured a strong absorption signal of the redshifted 21 cm radiation from neutral hydrogen at $z\approx 17$, requiring efficient formation of massive stars before then. In this paper, we investigate whether star formation in minihalos plays a significant role in establishing the early Ly$\alpha$ background required to produce the EDGES absorption feature. We find that Pop III stars are important in providing the necessary Ly$\alpha$-flux at high redshifts, and derive a best-fitting average Pop III stellar mass of $\sim$ 750M$_\odot{}$ per minihalo, corresponding to a star formation efficiency of 0.1%. Further, it is important to include baryon-dark matter streaming velocities in the calculation, to limit the efficiency of Pop~III star formation in minihalos. Without this effect, the cosmic dawn coupling between 21 cm spin temperature and that of the gas would occur at redshifts higher than what is implied by EDGES.

It has been shown that time delay of strong gravitational lensing is not only an effective cosmological probe to constrain the Hubble constant, the matter density and the curvature of the universe but also useful for breaking the degeneracy of the dark energy equation of state and thus provide complementarity with other popular probes, such as type Ia supernovae, baryon acoustic oscillations, and cosmic microwave background radiation. Interestingly, compared to traditional strong lensing systems where quasars act as sources, strongly lensed gravitational waves (GWs) from compact binary coalescence and their electromagnetic (EM) counterparts systems have been recently proposed to be more powerful for studying cosmology since GWs and their EM counterparts are transients. Short durations of GWs and their EM counterparts are very advantageous to time delay measurement and lens profile modeling. Here, in the framework of Chevalier-Polarski-Linder (CPL) parametrization, we investigate improvement of constraining power on dark energy equation of state due to including time delay measurements of strong lensed GW systems. It is suggested that, on the basis of the third generation ground-based detector, e.g. Einstein Telescope, adding time delay of only 30 strong lensed GW systems to type Ia supernovae and cosmic microwave background radiation can improve the dark energy figure of merit by a factor 2. For the traditional standard siren method where the uncertainties of luminosity distances of GWs are $\sim10\%$, a few $\times~10^4$ events are expected to present similar constraints. In the precision cosmology era, this progress is of great significance for studying the nature of dark energy.

We use cosmological hydrodynamic zoom-in simulations to study early structure formation in two dark matter (DM) cosmologies, the standard CDM model, and a thermal warm DM (WDM) model with a particle mass of $m_{\chi}c^{2}=3\ \mathrm{keV}$. We focus on DM haloes with virial masses $M\sim 10^{10}\ M_{\odot}$. We find that the first star formation activity is delayed by $\sim 200\ \mathrm{Myr}$ in the WDM model, with similar delays for metal enrichment and the formation of the second generation of stars. However, the differences between the two models in globally-averaged properties, such as star formation rate density and mean metallicity, decrease towards lower redshifts ($z\lesssim 10$). Metal enrichment in the WDM cosmology is restricted to dense environments, while low-density gas can also be significantly enriched in the CDM case. The free-free contribution from early structure formation at redshifts $z>6$ to the cosmic radio background (CRB) is $3_{-1.5}^{+13}\%$ ($8_{-3.5}^{+33}\%$) of the total signal inferred from radio experiments such as ARCADE 2, in the WDM (CDM) model. The direct detection of the $\mathrm{H_{2}}$ emission from early structure formation ($z\gtrsim 7.2$), originating from the low-mass haloes explored here, will be challenging even with the next generation of far-infrared space telescopes, unless the signal is magnified by at least a factor of 10 via gravitational lensing or shocks. However, more massive haloes with $M\gtrsim 10^{12}\ M_{\odot}$ may be observable for $z\gtrsim 10$, even without magnification, provided that our extrapolation from the scale of our simulated haloes is valid.

We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total envelope energy is nearly constant during this time because positive energy transferred to the gas from the core particles is counterbalanced by the negative binding energy from the closer proximity of the inner layers to the plunged-in secondary. During the subsequent slow spiral-in phase, energy continues to transfer to the envelope from the red giant core and secondary core particles. We also propose that relative motion between the centre of mass of the envelope and the centre of mass of the particles could account for the offsets of planetary nebula central stars from the nebula's geometric centre.

Streaming instability is a key mechanism in planet formation, clustering pebbles into planetesimals. It is triggered at a particular disk location where the local volume density of solids exceeds that of the gas. After their formation, planetesimals can grow by feeding from other planetesimals in the birth ring as well as by accreting inwardly drifting pebbles from the outer disk. To investigate the growth of planetesimals at a single location by the streaming instability, we test the conditions under which super-Earths are able to form within the lifetime of the gaseous disk. We modify the \texttt{Mercury} N-body code to trace the growth and dynamical evolution of a swarm of planetesimals at the ice line for a solar-mass star. Three distributions of planetesimal sizes are investigated: (i) a mono-dispersed population of 400 km radius planetesimals, (ii) a poly-dispersed populations of planetesimals from 200 km up to 1000 km, (iii) a bimodal distribution with a single runaway body and a swarm of smaller, 100 km size planetesimals. The mono-disperse population of 400 km size planetesimals cannot form $\gtrsim$ Earth mass protoplanets. Their velocity dispersions are quickly excited, which suppresses both planetesimal and pebble accretion. Planets can form from the poly-dispersed and bimodal distributions. In these circumstances, the two-component nature damps the random velocity of the large embryo by small planetesimals' dynamical friction, allowing the embryo to accrete pebbles efficiently when it approaches $10^{-2}$ Earth mass. We find that super-Earth planets are preferred to form when the pebble mass flux is higher, the disk turbulence is lower, or the Stokes number of the pebbles is higher.

The Survey of Ionized Gas of the Galaxy, Made with the Arecibo telescope (SIGGMA) provides a fully-sampled view of the radio recombination line (RRL) emission from the portion of the Galactic plane visible by Arecibo. Observations use the Arecibo L-band Feed Array (ALFA), which has a FWHM beam size of 3.4 arcmin. Twelve hydrogen RRLs from H163$\alpha$ to H174$\alpha$ are located within the instantaneous bandpass from 1225 MHz to 1525 MHz. We provide here cubes of average ("stacked") RRL emission for the inner Galaxy region $32 \le \ell \le 70$ degrees, $|b|\le1.5$ degree, with an angular resolution of 6 arcmin. The stacked RRL rms at 5.1 km/s velocity resolution is $\sim0.65$ mJy beam$^{-1}$, making this the most sensitive large-scale fully-sampled RRL survey extant. We use SIGGMA data to catalogue 319 RRL detections in the direction of 244 known HII regions, and 108 new detections in the direction of 79 HII region candidates. We identify 11 Carbon RRL emission regions, all of which are spatially coincident with known HII regions. We detect RRL emission in the direction of 14 of the 32 supernova remnants (SNRs) found in the survey area. This RRL emission frequently has the same morphology as the SNRs. The RRL velocities give kinematic distances in agreement with those found in the literature, indicating that RRLs may provide an additional tool to constrain distances to SNRs. Finally, we analyze the two bright star-forming complexes: W49 and W51. We discuss the possible origins of the RRL emission in directions of SNRs W49B and W51C.

We develop a physical framework for interpreting high-resolution images of pre planetary nebule ("prePNe") with pairs of candle shaped lobes. We use hydrodynamical models to infer the historical properties of the flows injected from the nucleus that shape the lobes into standard forms. First, we find a suitable set of parameters of a fast, collimated, tapered flow that is actively reshaped by an exterior slow AGB wind and that nicely fits the basic shape, kinematics, mass, and momenta of this class of prePNe. Next we vary the most influential parameters of this "baseline" model-such as density, speed, and geometry-to see how changes in the flow parameters affect the nebular observables after 900y. Several generic conlusions emerge, such as the injected flows that create the hollow candle-shaped lobes must be light, "tapered", and injected considerably faster than the lobe expansion speed. Multi-polar and starfish prePNe probably evolve from wide angle flows in which thin-shell instabilites corrugate their leading edges. We show how the common linear relationship of Doppler shift and position along the lobe is a robust outcome the interaction of tapered diverging streamlines with the lobes' curved walls. Finally we probe how magnetic fields affect the basline model by adding a toroidal field to the injected baseline flow. Examples of prePNe and PNe that may have been magnetically shaped are listed. We conclude that the light, field-free, tapered baseline flow model is an successful and universal pardigm for unravelling the histories of lobe formation in prePNe.

We study the orbital evolution of black hole (BH) binaries in quadruple systems, where the tertiary binary excites large eccentricity in the BH binary through Lidov-Kozai (LK) oscillations, causing the binary BHs to merge via gravitational radiation. For typical BH binaries with masses $m_{1,2}\simeq 20M_\odot-30M_\odot$ and initial semimajor axis $a_0\sim100$ AU (such that the binaries have no chance of merging by themselves within $\sim10^{10}$ yrs), we show that binary-binary interactions can significantly increase the LK window for mergers (the range of companion inclinations that allows the BH binary to merge within 10~Gyrs). This increase arises from a secular resonance between the LK oscillation of the BH binary and the nodal precession of the outer (binary-binary) orbit driven by the tertiary binary. Therefore, in the presence of tertiary binary, the BH merger fraction is increased to $10-30\%$, an order of magnitude larger than the merger fraction found in similar triple systems. Our analysis (with appropriate scalings) can be easily adapted to other configurations of systems, such as relatively compact BH binaries and moderately hierarchical triples, which may generate even higher merger fractions. Since the occurrence rate of stellar quadruples in the galactic fields is not much smaller than that of stellar triples, our result suggests that dynamically induced BH mergers in quadruple systems may be an important channel of producing BH mergers observed by LIGO/VIRGO.

Using Green Bank Telescope radio recombination line (RRL) data, we analyze the role of leaking radiation from HII regions in maintaining the ionization of the interstellar medium. We observed a sample of eight Galactic HII regions of various sizes, morphologies, and luminosities. For each region the hydrogen RRL intensity decreases roughly as a power-law with distance from the center of the region. This suggests that radiation leaking from the HII region is responsible for the majority of surrounding ionized gas producing RRL emission. Our results further indicate that the hydrogen RRL intensity appears to be fundamentally related to the HII region sizes traced by their photodissociation regions, such that physically smaller HII regions show a steeper decrease in intensity with increasing distance from the region centers. As a result, giant HII regions may have a much larger effect in maintaining the ionization of the interstellar medium. For six of the eight observed HII regions we find a decrease in the 4He+/H+ abundance ratio with increasing distance, indicating that He-ionizing photons are being absorbed within the ionization front of the HII region. There is enhanced carbon RRL emission toward directions with strong continuum background, suggesting that the carbon emission is amplified by stimulated emission.

We employ the supernova fallback disk model to simulate the spin evolution of isolated young neutron stars (NSs). We consider the submergence of the NS magnetic fields during the supercritical accretion stage and its succeeding reemergence. It is shown that the evolution of the spin periods and the magnetic fields in this model is able to account for the relatively weak magnetic fields of central compact objects and the measured braking indices of young pulsars. For a range of initial parameters, evolutionary links can be established among various kinds of NS sub-populations including magnetars, central compact objects and young pulsars. Thus, the diversity of young NSs could be unified in the framework of the supernova fallback accretion model.

Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.

A thermal soft X-ray component is often detected in low-level accreting neutron stars (NSs), but is not detected in low-level accreting stellar-mass black holes (BHs). In this paper, we investigate the origin of such a thermal soft X-ray component in the framework of the self-similar solution of the advection-dominated accretion flow (ADAF) around NSs. It is assumed that a fraction, $f_{\rm th}$, of the energy transferred onto the surface of the NS is thermalized at the surface of the NS as the soft photons to be scattered in the ADAF. We self-consistently calculate the structure and the corresponding emergent spectrum of the ADAF by considering the radiative coupling between the soft photons from the surface of the NS and the ADAF itself. We show that the Compton $y$-parameter of the ADAF for NSs is systematically lower than that of BHs. Meanwhile, we find that the temperature of the thermal soft X-ray component in NSs decreases with decreasing mass accretion rate, which is qualitatively consistent with observations. We test the effect of $f_{\rm th}$ on the structure, as well as the emergent spectrum of the ADAF. It is found that a change of $f_{\rm th}$ can significantly change the temperature of the thermal soft X-ray component as well as the spectral slope in hard X-rays. Finally, it is suggested that the value of $f_{\rm th}$ can be constrained by fitting the high-quality X-ray data, such as the $\it XMM$-$\it Newton$ spectrum between 0.5-10 keV in the future work.

Black hole (BH) mergers driven by gravitational perturbations of external companions constitute an important class of formation channels for merging BH binaries detected by LIGO. We have studied the orbital and spin evolution of binary BHs in triple systems, where the tertiary companion excites large eccentricity in the inner binary through Lidov-Kozai oscillations, causing the binary to merge via gravitational radiation. Using the single-orbital-averaged and double-orbital-averaged secular dynamics equations of triples, we determine the merger window (the range of companion inclinations that allows the inner binary to merge) and the merger fraction as a function of various system parameters. For typical BH binaries, the merger fraction increases rapidly with the perturber eccentricity $e_{\rm out}$ because of the octupole perturbation, reaching $10-20\%$ at $e_{\rm out}=0.9$. We derive the analytical expressions for the merger window and merger fraction for systems with negligible octupole effect, and apply them to neutron star binary mergers in triples. We also follow the spin evolution of the BHs during the companion-induced orbital decay. Starting from aligned spin axes, a wide range of final spin-orbit misalignment angle $\theta_{\rm sl}^{\rm f}$ can be generated when the binary enters the LIGO band. For systems where the octupole effect is small (such as equal mass binaries or $e_{\rm out}\sim 0$), the distribution of $\theta_{\rm sl}^{\rm f}$ peaks around $90^\circ$. As the octuple effect increases, a more isotropic distribution of final spin axis is produced. Overall, merging BH binaries produced by Lidov-Kozai oscillations in triples exhibit a unique distribution of the effective spin parameter $\chi_{\rm eff}$; this may be used to distinguish this formation channel from others.

A new accretion picture based on a small disk surrounding a black hole is developed for the wind-fed source Cyg X-1. The hard and soft spectral states of Cyg X-1 are interpreted in terms of co-spatial two component flows for the innermost region of an accretion disk. The state transitions result from the outward expansion and inward recession of this inner disk for the hard to soft and soft to hard transition respectively. The theoretical framework for state transitions in black hole X-ray binaries with high mass companions involving a change in the inner disk size, thus, differs from systems with low mass companions involving the change in the outer disk size. This fundamental difference stems from the fact that matter captured and supplied to the black hole in wind-fed systems has low specific angular momentum and is hot essentially heated in the bow and spiral shocks, whereas it has high specific angular momentum and is cool in Roche lobe overflow systems. The existence of a weak cool disk around the ISCO region in the hard state allows for the presence of a relativistically broadened Fe K line. The small disk fed by gas condensation forms without an extensive outer disk, precluding thermal instabilities and large outbursts, resulting in the lack of large amplitude outbursts and hysteresis effects in the light curve of high mass black hole X-ray binaries. Their relatively persistent X-ray emission is attributed to their wind-fed nature.

We describe work in progress by a collaboration of astronomers and statisticians developing a suite of Bayesian data analysis tools for extrasolar planet (exoplanet) detection, planetary orbit estimation, and adaptive scheduling of observations. Our work addresses analysis of stellar reflex motion data, where a planet is detected by observing the "wobble" of its host star as it responds to the gravitational tug of the orbiting planet. Newtonian mechanics specifies an analytical model for the resulting time series, but it is strongly nonlinear, yielding complex, multimodal likelihood functions; it is even more complex when multiple planets are present. The parameter spaces range in size from few-dimensional to dozens of dimensions, depending on the number of planets in the system, and the type of motion measured (line-of-sight velocity, or position on the sky). Since orbits are periodic, Bayesian generalizations of periodogram methods facilitate the analysis. This relies on the model being linearly separable, enabling partial analytical marginalization, reducing the dimension of the parameter space. Subsequent analysis uses adaptive Markov chain Monte Carlo methods and adaptive importance sampling to perform the integrals required for both inference (planet detection and orbit measurement), and information-maximizing sequential design (for adaptive scheduling of observations). We present an overview of our current techniques and highlight directions being explored by ongoing research.

Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes.

Coagulation theory predicts that micron-sized dust grains grow into pebbles which drift inward towards the star, when they reach sizes of mm-cm. When they cross the orbit of a planet, a fraction of these drifting pebbles will be accreted. In the pebble accretion mechanism, the combined effects of the planet's gravitational attraction and gas drag greatly increase the accretion rate. We calculate the pebble accretion efficiency $\varepsilon_{2D}$ -- the probability a pebble is accreted by the planet -- in the 2D limit (pebbles reside in the midplane). In particular, we investigate the dependence of $\varepsilon_{2D}$ on the planet eccentricity and its implications for planet formation models. We conduct N-body simulations to calculate the pebble accretion efficiency in both the local frame and the global frame. With the global method we investigate the pebble accretion efficiency when the planet is on an eccentric orbit. We find that the local and the global methods generally give consistent results. The efficiency increases with the planet's eccentricity once the relative velocity between the pebble and the planet is determined by the planet's eccentric velocity. At large eccentricities, however, the relative velocity becomes too large for pebble accretion. The efficiency then drops significantly and the accretion enters the ballistic regime. We present general expressions for $\varepsilon_{2D}$. Applying the obtained formula to the formation of a secondary planet, in resonance with an already-formed giant planet, we find that the embryo grows quickly due to its larger eccentricity. The maximum $\varepsilon_{2D}$ for a planet on an eccentric orbit is several times higher than for a planet on a circular orbit, but this increase gives the planet an important headstart and boosts its following mass growth.

Turbulence plays a key role in the transport of pebble-sized particles. It also affects the ability of pebbles to be accreted by protoplanets, because it stirs pebbles out of the disk midplane. In addition, turbulence can suppress pebble accretion once the relative velocities become too large for the settling mechanism to be viable. Following Paper I, we aim to quantify these effects by calculating the pebble accretion efficiency $\varepsilon$ using three-body simulations. To model the effect of turbulence on the pebbles, we derive a stochastic equation of motion (SEOM) applicable to stratified disk configurations. In the strong coupling limit (ignoring particle inertia) the limiting form of this equation agrees with previous works. We conduct a parameter study and calculate $\varepsilon$ in 3D, varying pebble and gas (turbulence) properties and accounting for the planet inclination. We find that strong turbulence suppresses pebble accretion through turbulent diffusion, agreeing within factors of order unity with previous works. Another reduction of $\varepsilon$ occurs when the turbulent rms motions are large and the settling mechanism fails. Efficiency-wise, the outer disk regions are more affected by turbulence than the inner regions. At the location of the H$_2$O iceline, planets around low-mass stars achieve much higher efficiencies. Including the results from Paper I, we present a framework to obtain $\varepsilon$ under general circumstances.

In this paper, we explained the observed $\Gamma_{\rm 2-10keV}-L_{\rm bol}/L_{\rm Edd}$ correlation in luminous active galactic nuclei within the framework of the condensation of the corona around a supermassive black hole (Liu et al. 2015; Qiao & Liu 2017). Specifically, we systemically test the effects of black hole mass $M$, the viscosity parameter $\alpha$, and the magnetic parameter $\beta$ (with magnetic pressure $p_{\rm m}=B^2/{8\pi}=(1-\beta)p_{\rm tot}$, $p_{\rm tot}=p_{\rm gas}+p_{\rm m}$) on the structure of the accretion disc and the corona, as well as the corresponding emergent spectra. It is found that the hard X-ray photon index $\Gamma_{\rm 2-10keV}$ nearly does not change with changing black hole mass $M$, or changing magnetic parameter $\beta$. Meanwhile, it is found that the geometry of the accretion flow, i.e., the relative configuration of the disc and corona, as well as the emergent spectra can be strongly affected by changing the value of $\alpha$. By comparing with a sample composed of 29 luminous active galactic nuclei with well constrained X-ray spectra and Eddington ratios, it is found that the observed $\Gamma_{\rm 2-10keV}-L_{\rm bol}/L_{\rm Edd}$ correlation can be well matched with a relatively bigger value of $\alpha$, i.e., $\alpha \sim 1$, as previously also suggested by Narayan (1996) for luminous accreting black holes.

We complete the formulation of the standard model of first star formation by exploring the possible impact of $\mathrm{LiH}$ cooling, which has been neglected in previous simulations of non-linear collapse. Specifically, we find that at redshift $z\gtrsim 5$, the cooling by $\mathrm{LiH}$ has no effect on the thermal evolution of shocked primordial gas, and of collapsing primordial gas into minihaloes or relic HII regions, even if the primordial lithium abundance were enhanced by one order of magnitude. Adding the most important lithium species to a minimum network of primordial chemistry, we demonstrate that insufficient $\mathrm{LiH}$ is produced in all cases considered, about $[\mathrm{LiH/Li}]\sim 10^{-9}$ for $T\lesssim 100$ K. Indeed, $\mathrm{LiH}$ cooling would only be marginally significant in shocked primordial gas for the highly unlikely case that the $\mathrm{LiH}$ abundance were increased by nine orders of magnitude, implying that $all$ lithium would have to be converted into $\mathrm{LiH}$. In this study, photo-destruction processes are not considered, and the collisional disassociation rate of $\mathrm{LiH}$ is possibly underestimated, rendering our results an extreme upper limit. Therefore, the cooling by $\mathrm{LiH}$ can safely be neglected for the thermal evolution of Population~III star-forming gas.

M2-9, or the "Butterfly Nebula" is one of the most iconic outflow sources from an evolved star. In this paper we present a hydrodynamic model of M2-9 in which the nebula is formed and shaped by a steady, low-density ("light"), mildly collimated "spray" of gas injected at 200 km s^-1 that interacts with a far denser, intrinsically simple pre-existing AGB wind has slowly formed all of the complex features within M2-9's lobes (including the knot pairs N3/S3 and N4/S4 at their respective leading edges, and the radial gradient of Doppler shifts within 20" of the nucleus). We emphasize that the knot pairs are not ejected from the star but formed in situ. In addition, the observed radial speed of the knots is only indirectly related to the speed of the gas injected by the star. The model allows us to probe the early history of the wind geometry and lobe formation. We also formulate a new estimate of the nebular distance D = 1.3 kpc. The physical mechanism that accounts for the linear radial speed gradient in M2-9 applies generally to many other pre planetary nebulae whose hollow lobes exhibit similar gradients along their edges.

Proxima Centauri is known as the closest star from the Sun. Recently, radial velocity observations revealed the existence of an Earth-mass planet around it. With an orbital period of ~11 days, the surface of Proxima Centauri b is temperate and might be habitable. We took a photometric monitoring campaign to search for its transit, using the Bright Star Survey Telescope at the Zhongshan Station in Antarctica. A transit-like signal appearing on 2016 September 8th, is identified tentatively. Its midtime, $T_{C}=2,457,640.1990\pm0.0017$ HJD, is consistent with the predicted ephemeris based on RV orbit in a 1$\sigma$ confidence interval. Time-correlated noise is pronounced in the light curve of Proxima Centauri, affecting detection of transits. We develop a technique, in a Gaussian process framework, to gauge the statistical significance of potential transit detection. The tentative transit signal reported here, has a confidence level of $2.5\sigma$. Further detection of its periodic signals is necessary to confirm the planetary transit of Proxima Centauri b. We plan to monitor Proxima Centauri in next Polar night at Dome A in Antarctica, taking the advantage of continuous darkness. \citet{Kipping17} reported two tentative transit-like signals of Proxima Centauri b, observed by the Microvariability and Oscillation of Stars space Telescope in 2014 and 2015, respectively. The midtransit time of our detection is 138 minutes later than that predicted by their transit ephemeris. If all the signals are real transits, the misalignment of the epochs plausibly suggests transit timing variations of Proxima Centauri b induced by an outer planet in this system.

We report our independent GeV gamma-ray study of the young shell-type supernova remnant (SNR) Kes 73 which harbors a central magnetar, and CO-line millimeter observations toward the SNR. Using 7.6 years of Fermi-LAT observation data, we detected an extended gamma-ray source ("source A") with the centroid on the west of the SNR, with a significance of 21.6 sigma in 0.1-300 GeV and an error circle of 5.4 arcminute in angular radius. The gamma-ray spectrum cannot be reproduced by a pure leptonic emission or a pure emission from the magnetar, and thus a hadronic emission component is needed. The CO-line observations reveal a molecular cloud (MC) at V_LSR~90 km/s, which demonstrates morphological correspondence with the western boundary of the SNR brightened in multiwavelength. The 12CO (J=2-1)/12CO (J=1-0) ratio in the left (blue) wing 85-88 km/s is prominently elevated to ~1.1 along the northwestern boundary, providing kinematic evidence of the SNR-MC interaction. This SNR-MC association yields a kinematic distance 9 kpc to Kes 73. The MC is shown to be capable of accounting for the hadronic gamma-ray emission component. The gamma-ray spectrum can be interpreted with a pure hadronic emission or a magnetar+hadronic hybrid emission. In the case of pure hadronic emission, the spectral index of the protons is 2.4, very similar to that of the radio-emitting electrons, essentially consistent with the diffusive shock acceleration theory. In the case of magnetar+hadronic hybrid emission, a magnetic field decay rate >= 10^36 erg/s is needed to power the magnetar's curvature radiation.

Accretion flows toward black holes can be of a quite different nature, described as an optically thick cool gas flow in a disk for high accretion rates or as a hot coronal optically thin gas flow for low accretion rates, possibly affected by outflowing gas. The detection of broad iron emission lines in active galactic nuclei (AGN) indicates the coexistence of corona and disk. The appearance and relative strength of such flows essentially depends on their interaction. Liu et al. suggested that condensation of gas from the corona to the disk allows to understand accretion flows of comparable strength of emission. Matter inflow due to gravitational capture of gas is important for the condensation process. We discuss observational features predicted by the model. Data from simultaneous observations of AGN with {\it {Swift's}} X-ray and UV-optical telescopes are compared with the theoretical predictions. The frequent detection of broad iron K$\alpha$ emission lines and the dependence of the emitted spectra on the Eddington ratio, described by the values of the photon index $\Gamma$ and the two-point spectral index $\alpha_{\rm{ox}}$ are in approximate agreement with the predictions of the condensation model; the latter, however, with a large scatter. The model further yields a coronal emission concentrated in a narrow inner region as is also deduced from the analysis of emissivity profiles. The accretion flows in bright AGN could be described by the accretion of stellar wind or interstellar medium and its condensation into a thin disk.

The X-ray emission from bright active galactic nuclei (AGNs) is believed to originate in a hot corona lying above a cold, geometrically thin accretion disk. A highly concentrated corona located within $\sim10$ gravitational radii above the black hole is inferred from observations. Based on the accretion of interstellar medium/wind, a disk corona model has been proposed in which the corona is well coupled to the disk by radiation, thermal conduction, as well as by mass exchange \citep{Liu2015, Qiao2017}. Such a model avoids artificial energy input to the corona and has been used to interpret the spectral features observed in AGN. In this work, it is shown that the bulk emission size of the corona is very small for the extended accretion flow in our model. More than 80\% of the hard X-ray power is emitted from a small region confined within 10 Schwarzschild radii around a non-spinning black hole, which is expected to be even smaller accordingly for a spinning black hole. Here, the corona emission is more extended at higher Eddington ratios. The compactness parameter of the corona, $l={L\over R}{\sigma_{\rm T}\over m_{\rm e} c^3}$, is shown to be in the range of 1-33 for Eddington ratios of 0.02 - 0.1. Combined with the electron temperature in the corona, this indicates that electron--positron pair production is not dominant in this regime. A positive relation between the compactness parameter and photon index is also predicted. By comparing the above model predictions with observational features, we find that the model is in agreement with observations.

We study the effect of external companion on the orbital and spin evolution of merging black-hole (BH) binaries. An sufficiently close by and inclined companion can excite Lidov-Kozai (LK) eccentricity oscillations in the binary, thereby shortening its merger time. During such LK-enhanced orbital decay, the spin axis of the BH generally exhibits chaotic evolution, leading to a wide range ($0^\circ$-$180^\circ$) of final spin-orbit misalignment angle from an initially aligned configuration. For systems that do not experience eccentricity excitation, only modest ($\lesssim 20^\circ$) spin-orbit misalignment can be produced, and we derive an analytic expression for the final misalignment using the principle of adiabatic invariance. The spin-orbit misalignment directly impacts the gravitational waveform, and can be used to constrain the formation scenarios of BH binaries and dynamical influences of external companions.

Context.The Kepler mission has provided a large sample to statistically analyze the orbital properties of the super-Earth planets. We hypothesize that these planets formed early and consider the problem of matching planet formation theory to the current observations. Two scenarios, disk migration and in-situ formation, have been proposed to explain their origin. In the migration scenario planets migrate inward due to planet-disk interaction, whereas in the in-situ scenario planets assemble locally. Therefore, planets formed by migration are expected to end up in resonances, whereas those formed in-situ are expected to stay in short period ratios and in non-resonant orbits. Both predictions are at odds with observations. Aims. We investigate whether a preferred formation scenario can be identified through a comparison between the magnetospheric rebound model and the Kepler data. Methods. We conduct N-body simulations of two-planet systems during the disk dispersal phase, and make a statistical comparison between the simulations and the Kepler observations. Results. Comparing the two scenarios, we find that magnetospheric rebound tends to erase the difference in the orbital configuration that was initially presented. After disk dispersal, not all planets are in resonance in the migration scenario, whereas planets do not remain in compact configurations in the in-situ scenario. In both scenarios, the orbits of planets increase with the cavity expansion, and their period ratios have a wider distribution. Conclusions. From a statistical perspective, the magnetospheric rebound model reproduces several observed properties of Kepler planets, such as the significant number of planets are not in resonances and planet pairs can end up at large period ratios. The disparity in orbital configuration between the two formation scenarios is substantially reduced after disk dispersal.

This paper addresses maximum likelihood (ML) estimation based model fitting in the context of extrasolar planet detection. This problem is featured by the following properties: 1) the candidate models under consideration are highly nonlinear; 2) the likelihood surface has a huge number of peaks; 3) the parameter space ranges in size from a few to dozens of dimensions. These properties make the ML search a very challenging problem, as it lacks any analytical or gradient based searching solution to explore the parameter space. A population based searching method, called estimation of distribution algorithm (EDA), is adopted to explore the model parameter space starting from a batch of random locations. EDA is featured by its ability to reveal and utilize problem structures. This property is desirable for characterizing the detections. However, it is well recognized that EDAs can not scale well to large scale problems, as it consists of iterative random sampling and model fitting procedures, which results in the well-known dilemma curse of dimensionality. A novel mechanism to perform EDAs in interactive random subspaces spanned by correlated variables is proposed and the hope is to alleviate the curse of dimensionality for EDAs by performing the operations of sampling and model fitting in lower dimensional subspaces. The effectiveness of the proposed algorithm is verified via both benchmark numerical studies and real data analysis.

We study the secular dynamics of compact binaries (consisting of white dwarfs, neutron stars or black holes) with tertiary companions in hierarchical triple systems. As the inner binary (with initially negligible eccentricity) undergoes orbital decay due to gravitational radiation, its eccentricity can be excited by gravitational forcing from the tertiary. This excitation occurs when the triple system passes through an "apsidal precession resonance", when the precession rate of the inner binary, driven by the gravitational perturbation of the external companion and general relativity, matches the precession rate of the outer binary. The eccentricity excitation requires the outer companion to be on an eccentric orbit, with the mutual inclination between the inner and outer orbits less than $\sim 40^\circ$. Gravitational wave (GW) signals from the inner binary can be significantly modified as the system evolves through the apsidal precession resonance. For some system parameters (e.g., a white dwarf binary with a brown dwarf tertiary), the resonance can happen when the binary emits GWs in the $10^{-4}-10^{-1}$~Hz range (the sensitivity band of LISA).

We present a detailed hydrodynamic model that matches the present structure of the well-observed preplanetary nebula OH231.8+04.2. The purpose of the model is to present a physically justified and coherent picture of its evolutionary history from about 100 years of the start of the formation of its complex outer structures to the present. We have adopted a set of initial conditions that are heavily constrained by high-quality observations of its present structure and kinematics. The shaping of the nebula occurs while the densities of the flows are "light": less than the surrounding AGB-wind environment. The simulations show that pairs of essentially coeval clumps and sprays of the same extent and density but different outflow speeds sculpted both the pair of thin axial flow "or spine" and the bulbs. The total ejected mass and momentum in the best fit model are surprisingly large -- 3 M_sun and 2.2 x 10^41 gm cm s-1, respectively -- however, these values are reduced by up to a factor of ten in other models that fit the data almost as well. Our ultimate goal is to combine the present model results of masses, momenta, flow speeds, and flow geometries for OH231 with those of other models to be published in the future in order to find common attributes of their ejection histories.

TRAPPIST-1 is a nearby 0.08 M M-star, which was recently found to harbor a planetary system of at least seven Earth-mass planets, all within 0.1 au. The configuration confounds theorists as the planets are not easily explained by either in situ or migration models. In this Paper we present a scenario for the formation and orbital architecture of the TRAPPIST-1 system. In our model, planet formation starts at the H2O iceline, where pebble-size particles -- whose origin is the outer disk -- concentrate to trigger streaming instabilities. After their formation, planetary embryos quickly mature by pebble accretion. Planet growth stalls at Earth masses, where the planet's gravitational feedback on the disk keeps pebbles at bay. Planets are transported by Type I migration to the inner disk, where they stall at the magnetospheric cavity and end up in mean motion resonances. During disk dispersal, the cavity radius expands and the inner-most planets escape resonance. We argue that the model outlined here can also be applied to other compact systems and that the many close-in super-Earth systems are a scaled-up version of TRAPPIST-1. We also hypothesize that few close-in compact systems harbor giant planets at large distances, since they would have stopped the pebble flux from the outer disk.

We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.

In this paper we present an analysis of absorption line variability in mini-BAL quasar LBQS 1206+1052. The SDSS spectrum demonstrates that the absorption troughs can be divided into two components of blueshift velocities of $\sim$700 km s$^{-1}$ and $\sim$1400 km s$^{-1}$ relative to the quasar rest-frame. The former component shows rare Balmer absorption, which is an indicator of high density absorbing gas, thus the quasar is worth follow-up spectroscopic observations. Our follow-up optical and near-infrared spectra using MMT, YFOSC, TripleSpec and DBSP reveal that the strengths of the absorption lines vary for both of the two components, while the velocities do not change. We reproduce all of the spectral data by assuming that only the ionization state of the absorbing gas is variable and that all other physical properties are invariable. The variation of ionization is consistent with the variation of optical continuum from the V-band light-curve. Additionally, we can not interpret the data by assuming that the variability is due to a movement of the absorbing gas. Therefore, our analysis strongly indicates that the absorption line variability in LBQS 1206+1052 is photoionization-driven. As shown from photo-ionization simulations, the absorbing gas with blueshift velocity of $\sim$700 km s$^{-1}$ has a density in the range of $10^9$ to $10^{10}$ cm$^{-3}$ and a distance of $\sim$1 pc, and the gas with blueshift velocity of $\sim$1400 km s$^{-1}$ has a density of $10^3$ cm$^{-3}$ and a distance of $\sim$1 kpc.

The Kepler mission has discovered that multiple close-in super-Earth planets are common around solar-type stars, but their period ratios do not show strong pile-ups near mean motion resonances (MMRs). One scenario is that super-Earths form in a gas-rich disk, and they interact gravitationally with the surrounding gas, inducing their orbital migration. Disk migration theory predicts, however, that planets would end up at resonant orbits due to their differential migration speed. Motivated by the discrepancy between observation and theory, we seek for a mechanism that moves planets out of resonances. We examine the orbital evolution of planet pairs near the magnetospheric cavity during the gas disk dispersal phase. Our study determines the conditions under which planets can escape resonances. We perform two-planet N-body simulations, varying the planet masses, stellar magnetic field strengths, disk accretion rates and gas disk depletion timescales. As planets migrate outward with the expanding magnetospheric cavity, their dynamical configurations can be rearranged. Migration of planets is substantial (minor) in a massive (light) disk. When the outer planet is more massive than the inner planet, the period ratio of two planets increases through outward migration. On the other hand, when the inner planet is more massive, the final period ratio tends to remain similar to the initial one. Larger stellar magnetic field strengths result in planets stopping their migration at longer periods. We highlight \textit{magnetospheric rebound} as an important ingredient able to reconcile disk migration theory with observations. Even when planets are trapped into MMR during the early gas-rich stage, subsequent cavity expansion would induce substantial changes to their orbits, moving them out of resonance.

Bell's theorem states that some predictions of quantum mechanics cannot be reproduced by a local-realist theory. That conflict is expressed by Bell's inequality, which is usually derived under the assumption that there are no statistical correlations between the choices of measurement settings and anything else that can causally affect the measurement outcomes. In previous experiments, this "freedom of choice" was addressed by ensuring that selection of measurement settings via conventional "quantum random number generators" was space-like separated from the entangled particle creation. This, however, left open the possibility that an unknown cause affected both the setting choices and measurement outcomes as recently as mere microseconds before each experimental trial. Here we report on a new experimental test of Bell's inequality that, for the first time, uses distant astronomical sources as "cosmic setting generators." In our tests with polarization-entangled photons, measurement settings were chosen using real-time observations of Milky Way stars while simultaneously ensuring locality. Assuming fair sampling for all detected photons, and that each stellar photon's color was set at emission, we observe statistically significant $\gtrsim 7.31 \sigma$ and $\gtrsim 11.93 \sigma$ violations of Bell's inequality with estimated $p$-values of $ \lesssim 1.8 \times 10^{-13}$ and $\lesssim 4.0 \times 10^{-33}$, respectively, thereby pushing back by $\sim$600 years the most recent time by which any local-realist influences could have engineered the observed Bell violation.

The evolution of main sequence binaries resided in the galactic centre is influenced a lot by the central super massive black hole (SMBH). Due to this perturbation, the stars in a dense environment are likely to experience mergers or collisions through secular or non-secular interactions. In this work, we study the dynamics of the stellar binaries at galactic center, perturbed by another distant SMBH. Geometrically, such a four-body system is supposed to be decomposed into the inner triple (SMBH-star-star) and the outer triple (SMBH-stellar binary-SMBH). We survey the parameter space and determine the criteria analytically for the stellar mergers and the tidal disruption events (TDEs). For a relative distant and equal masses SMBH binary, the stars have more opportunities to merge as a result from the Lidov-Kozai(LK) oscillations in the inner triple. With a sample of tight stellar binaries, our numerical experiments reveal that a significant fraction of the binaries, ~70 per cent, experience merger eventually. Whereas the majority of the stellar TDEs are likely to occur at a close periapses to the SMBH, induced by the outer Kozai effect. The tidal disruptions are found numerically as many as ~10 per cent for a close SMBH binary that is enhanced significantly than the one without the external SMBH. These effects require the outer perturber to have an inclined orbit (>=40 degree) relatively to the inner orbital plane and may lead to a burst of the extremely astronomical events associated with the detection of the SMBH binary.

Observationally, it is found that there is a strong correlation between the hard X-ray photon index $\Gamma$ and the Compton reflection scaling factor $\Re$ in active galactic nuclei. In this paper, we propose that the $\Gamma-\Re$ correlation can be explained within the framework of the condensation of the hot corona onto the cold accretion disc around a supermassive black hole. In the model, it is presumed that, initially, a vertically extended hot gas (corona) is supplied to the central supermassive black hole by capturing the interstellar medium and stellar wind. In this scenario, when the initial mass accretion rate $\dot M/ \dot M_{\rm Edd} \gtrsim 0.01$, at a critical radius $r_{\rm d}$, part of the hot gas begins to condense onto the equatorial disc plane of the black hole, forming an inner cold accretion disc. Then the matter is accreted in the form of the disc-corona structure extending down to the innermost stable circular orbits of the black hole. The size of the inner disc is determined by the initial mass accretion rate. With the increase of the initial mass accretion rate, the size of the inner disc increases, which results in both the increase of the Compton reflection scaling factor $\Re$ and the increase of the hard X-ray photon index $\Gamma$. By comparing with a sample of Seyfert galaxies with well-fitted X-ray spectra, it is found that our model can roughly explain the observations. Finally, we discuss the possibility to apply our model to high mass X-ray binaries, which are believed to be fueled by the hot wind from the companion star.

It is believed that the hard X-ray emission in the luminous active galactic nuclei (AGNs) is from the hot corona above the cool accretion disk. However, the formation of the corona is still debated. Liu et al. investigated the spectrum of the corona heated by the reconnection of the magnetic field generated by dynamo action in the thin disk and emerging into the corona as a result of buoyancy instability. In the present paper, we improve this model to interpret the observed relation of the hard X-ray spectrum becoming softer at higher accretion rate in luminous AGNs. The magnetic field is characterized by $\beta_{\rm 0}$, i.e., the ratio of the sum of gas pressure and radiation pressure to magnetic pressure in the disk ($\beta_{\rm 0}=(P_{\rm g,d}+P_{\rm r,d})/P_{\rm B}$). Besides, both the intrinsic disk photons and reprocessed photons by the disk are included as the seed photons for inverse Compton scattering. These improvements are crucial for investigating the effect of magnetic field on the accretion disk-corona when it is not clear whether the radiation pressure or gas pressure dominates in thin disk. We change the value of $\beta_{\rm 0}$ in order to constrain the magnetic field in the accretion disk. We find that the energy fraction released in the corona ($f$) gradually increases with the decrease of $\beta_{\rm 0}$ for the same accretion rate. When $\beta_{\rm 0}$ decreases to less than 50, the structure and spectrum of the disk-corona is independent on accretion rate, which is similar to the hard spectrum found in Liu et al.(2003). Comparing with the observational results of the hard X-ray bolometric correction factor in a sample of luminous AGNs, we suggest that the value of $\beta_{\rm 0}$ is about 100-200 for $\alpha=0.3$ and the energy fraction $f$ should be larger than $30\%$ for hard X-ray emission.

Emission lines from the broad emission line region (BELR) and the narrow emission line region (NELR) of active galactic nuclei (AGNs) are extensively studied. However, between these two regions emission lines are rarely detected. We present a detailed analysis of a quasar SDSS J232444.80-094600.3 (SDSS J2324$-$0946), which is remarkable for its strong intermediate-width emission lines (IELs) with FWHM $\approx$ 1800 \kmps. The IEL component is presented in different emission lines, including the permitted lines \lya\ $\lambda$1216, \civ\ $\lambda$1549, semiforbidden line \ciii\ $\lambda$1909, and forbidden lines \oiii\ $\lambda\lambda$4959, 5007. With the aid of photo-ionization models, we found that the IELs are produced by gas with a hydrogen density of $n_{\rm H} \sim 10^{6.2}-10^{6.3}~\rm cm^{-3}$, a distance to the central ionizing source of $R \sim 35-50$ pc, a covering factor of CF $\sim$ 6\%, and a dust-to-gas ratio of $\leq 4\%$ times of SMC. We suggest that the strong IELs of this quasar are produced by nearly dust-free and intermediate-density gas located at the skin of the dusty torus. Such strong IELs, served as a useful diagnose, can provide an avenue to study the properties of gas between the BELR and the NELR.

Radial velocity and transit surveys have found that the fraction of FGKM stars with close-in super-Earth(s) ($\eta_\oplus$) is around $30 \%- 50\%$, independent of the stellar mass $M_\ast$ and metallicity $Z_\ast$. In contrast, the fraction of solar-type stars harboring one or more gas giants ($\eta_J $) with masses $M_{\rm p} > 100 \ M_\oplus $ is nearly $ 10\%-15\%$, and it appears to increase with both $M_\ast$ and $Z_\ast$. Regardless of the properties of their host stars, the total mass of some multiple super-Earth systems exceeds the core mass of Jupiter and Saturn. We suggest that both super-Earths and supercritical cores of gas giants were assembled from a population of embryos that underwent convergent type I migration from their birthplaces to a transition location between viscously heated and irradiation heated disk regions. We attribute the cause for the $\eta_\oplus$-$\eta_{\rm J}$ dichotomy to conditions required for embryos to merge and to acquire supercritical core mass ($M_c \sim 10 \ M_\oplus$) for the onset of efficient gaseous envelope accretion. We translate this condition into a critical disk accretion rate, and our analysis and simulation results show that it weakly depends on $M_\ast$ and decreases with metallicity of disk gas $Z_{\rm d}$. We find that embryos are more likely to merge into supercritical cores around relatively massive and metal-rich stars. This dependence accounts for the observed $\eta_{\rm J}$-$M_\ast$. We also consider the $Z_{\rm d}$-$Z_\ast$ dispersed relationship and reproduce the observed $\eta_J$-$Z_\ast$ correlation.

Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.

Quasi-periodic oscillation (QPO) detected in the X-ray radiation of black hole X-ray binaries (BHXBs) is thought to originate from dynamical processes in the close vicinity of the black holes (BHs), and thus carries important physical information therein. Such a feature is extremely rare in active galactic nuclei (AGNs) with supermassive BHs. Here we report on the detection of a possible X-ray QPO signal with a period of 3800\,s at a confidence level $>99.99\%$ in the narrow-line Seyfert 1 galaxy (NLS1) 1H~0707-495 in one data set in 0.2-10\,keV taken with {\it XMM-Newton}. The statistical significance is higher than that of most previously reported QPOs in AGNs. The QPO is highly coherent (quality factor $Q=\nu/\Delta\nu \geqslant 15$) with a high rms fractional variability ($\sim15\%$). A comprehensive analysis of the optical spectra of this AGN is also performed, yielding a central BH mass $5.2\times10^6\,M_{\odot}$ from the broad emission lines based on the scaling relation. The QPO follows closely the known frequency-BH mass relation, which spans from stellar-mass to supermassive BHs. The absence of the QPO in other observations of the object suggests it a transient phenomenon. We suggest that the (high-frequency) QPOs tend to occur in highly accreting BH systems, from BHXBs to supermassive BHs. Future precise estimation of the BH mass may be used to infer the BH spin from the QPO frequency.

Lidov-Kozai oscillations of planets in stellar binaries, combined with tidal dissipation, can lead to the formation of hot Jupiters (HJs) or tidal disruption of planets. Recent population synthesis studies have found that the fraction of systems resulting in HJs (F_HJ) depends strongly on the planet mass, host stellar type and tidal dissipation strength, while the total migration fraction F_mig = F_ HJ + F_dis (including both HJ formation and tidal disruption) exhibits much weaker dependence. We present an analytical method for calculating F_HJ and F_mig in the Lidov-Kozai migration scenario. The key ingredient of our method is to determine the critical initial planet-binary inclination angle that drives the planet to reach sufficiently large eccentricity for efficient tidal dissipation or disruption. This calculation includes the effects of the octupole potential and short-range forces on the planet. Our analytical method reproduces the resulting planet migration/disruption fractions from population synthesis, and can be easily implemented for various planet, stellar/companion types, and for different distributions of initial planetary semi-major axes, binary separations and eccentricities. We extend our calculations to planets in the super-Earth mass range and discuss the conditions for such planets to survive Lidov-Kozai migration and form close-in rocky planets.

The existence of intermediate-width emission line regions (IELRs) in active galactic nuclei has been discussed for over two decades. A consensus, however, is yet to be arrived at due to the lack of convincing evidence for their detection. We present a detailed analysis of the broadband spectrophotometry of the partially obscured quasar OI 287. The ultraviolet intermediate-width emission lines (IELs) are very prominent, in high contrast to the corresponding broad emission lines (BELs) which are heavily suppressed by dust reddening. Assuming that the IELR is virialized, we estimated its distance to the central black hole of $\sim 2.9$ pc, similar to the dust sublimation radius of $\sim 1.3$ pc. Photo-ionization calculations suggest that the IELR has a hydrogen density of $\sim 10^{8.8}-10^{9.4} ~ \rm cm^{-3}$, within the range of values quoted for the dusty torus near the sublimation radius. Both its inferred location and physical conditions suggest that the IELR originates from the inner surface of the dusty torus. In the spectrum of this quasar, we identified only one narrow absorption-line system associated with the dusty material. With the aid of photo-ionization model calculations, we found that the obscuring material might originate from an outer region of the dusty torus. We speculate that the dusty torus, which is exposed to the central ionizing source, may produce IELs through photo-ionization processes, while also obscure BELs as a natural "coronagraph". Such a "coronagraph" could be found in a large number of partially obscured quasars and be a useful tool to study IELRs.

Numerical simulations suggest that merging double white dwarfs (WDs) may produce a newborn neutron star surrounded by a fossil disk. We investigate the evolution of the fossil disk following the coalescence of double WDs. We demonstrate that the evolution can be mainly divided into four phases: the slim disk phase (with time $\lesssim$ 1 yr), the inner slim plus outer thin disk phase ($\sim 10-\DP{6}$ yr), the thin disk phase ($\sim \DP{2}-\DP{7}$ yr), and the inner advection-dominated accretion flow plus outer thin disk phase, given the initial disk mass $\sim 0.05-0.5\,M_{\sun}$ and the disk formation time $10^{-3}-1$ s. Considering possible wind mass loss from the disk, we present both analytic formulae and numerically calculated results for the disk evolution, which is sensitive to the condition that determines the location of the outer disk radius. The systems are shown to be very bright in X-rays in the early phase, but quickly become transient within $\lesssim$ 100 yr, with peak luminosities decreasing with time. We suggest that they might account for part of the very faint X-ray transients around the Galactic center region, which generally require a very low mass transfer rate.

We report on the identification of a new $\gamma$-ray-emitting narrow-line Seyfert 1 (NLS1) galaxy, SDSS J122222.55+041315.7, which increases the number of known objects of this remarkable but rare type of active galactic nuclei (AGN) to seven. Its optical spectrum, obtained in the Sloan Digital Sky Survey-Baryon Oscillation Spectroscopic Survey, reveals a broad H $\beta$ emission line with a width (FWHM) of 1734$\pm$104 km s$^{-1}$. This, along with strong optical Fe II multiplets [$R_{4570}=0.9$] and a weak [O III] $\lambda 5007$ emission line, makes the object a typical NLS1. On the other hand, the source exhibits a high radio brightness temperature, rapid infrared variability, and a flat X-ray spectrum extending up to $\sim$200 keV. It is associated with a luminous $\gamma$-ray source detected significantly with {\it Fermi}/LAT. Correlated variability with other wavebands has not yet been tested. The spectral energy distribution can be well modelled by a one-zone leptonic jet model. This new member is by far the most distant $\gamma$-ray-emitting NLS1, at a redshift of $z=0.966$.

We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.

Hadron emission from supernova remnant (SNR)-molecular cloud (MC) association systems has been widely regarded as a probe of the shock-accelerated cosmic-ray protons. Here, we report on the detection of a {\gamma}-ray emission source, with a significance of 24{\sigma} in 0.2-300 GeV, projected to lie to the northwest of the thermal composite SNR Kesteven 41, using 5.6 years of Fermi-Large Area Telescope (LAT) observation data. No significant long-term variability in the energy range 0.2--300 GeV is detected around this source. The 3{\sigma} error circle, 0.09 degree; in radius, covers the 1720MHz OH maser and is essentially consistent with the location of the V_{LSR} ~-50 km/s MC with which the SNR interacts. The source emission has an exponential cutoff power-law spectrum with a photon index of 1.9+/-0.1 and a cutoff energy of 4.0+/-0.9 GeV, and the corresponding 0.2-300 GeV luminosity is ~1.3*10^36 erg/s at a distance of 12 kpc. There is no radio pulsar in the 3{\sigma} circle responsible for the high {\gamma}-ray luminosity. While the inverse Compton scattering scenario would lead to a difficulty in the electron energy budget, the source emission can naturally be explained by the hadronic interaction between the relativistic protons accelerated by the shock of SNR Kesteven 41 and the adjacent northwestern MC. In this paper, we present a list of Galactic thermal composite SNRs detected at GeV {\gamma}-ray energies by Fermi-LAT is presented.

Recent studies of active galactic nuclei (AGN) found a statistical inverse linear scaling between the X-ray normalized excess variance $\sigma_{\rm rms}^2$ (variability amplitude) and the black hole mass spanning over $M_{\rm BH}=10^6- 10^9\ M_{\odot}$. Being suggested to have a small scatter, this scaling relation may provide a novel method to estimate the black hole mass of AGN. However, a question arises as to whether this relation can be extended to the low-mass regime below $\sim10^6\ M_{\odot}$. If confirmed, it would provide an efficient tool to search for AGN with low-mass black holes using X-ray variability. This paper presents a study of the X-ray excess variances for a sample of AGN with black hole masses in the range of $10^5- 10^6\ M_{\odot}$ observed with {\it XMM-Newton} and {\it ROSAT}, including data both from the archives and from newly preformed observations. It is found that the relation is no longer a simple extrapolation of the linear scaling; instead, the relation starts to flatten at $\sim10^6\ M_{\odot}$ toward lower masses. Our result is consistent with the recent finding of \citet{L15}. Such a flattening of the $M_{\rm BH}-\sigma_{\rm rms}^2$ relation is actually expected from the shape of the power spectrum density of AGN, whose break frequency is inversely scaled with the mass of black holes.

As a radio-loud narrow-line Seyfert 1 galaxy (NLS1) detected by Fermi/LAT in GeV $\gamma$-rays, 1H 0323+342 is a remarkable Active Galactic Nucleus (AGN) showing properties characteristic of both NLS1s and blazars. Here we present results of simultaneous X-ray and UV/optical monitoring observations on 1H 0323+342 taken with the UV/Optical Telescope (UVOT) and X-ray Telescope (XRT) onboard the Swift satellite over six years from 2006. Overall, the object showed statistically correlated variations in both the UV and X-ray bands on timescales of years as well as on timescales of days. A deep Suzaku observation reveals X-ray variability on timescales as short as a few tens of thousand seconds, and an X-ray spectrum typical of Seyfert galaxies. The broad-band spectral energy distribution, for which the data of UV and X-ray observations taken on 2009 July 26-27 were used, can be well modeled with a simple one-zone leptonic jet model plus accretion disk/corona emission. The latter is predominantly responsible for the UV/optical and X-ray (0.3-10 keV) emission and their observed variations. The correlated UV-X-ray variability on the timescale of days is consistent with reprocessing of the X-ray radiation by the accretion disk. The shortest timescale and large normalized excess variance of the X-ray variability detected with Suzaku suggest a relatively small black hole mass of the order of $10^7M_{\odot}$, consistent with the estimation based on the broad H$\beta$ line in the optical.

It is commonly believed that the optical/UV and X-ray emissions in luminous AGN are produced in an accretion disk and an embedded hot corona respectively. The inverse Compton scattering of disk photons by hot electrons in the corona can effectively cool the coronal gas if the mass supply is predominantly via a cool disk like flow as in BHXRBs. Thus, the application of such a model to AGNs fails to produce their observed X-ray emission. As a consequence, a fraction of disk accretion energy is usually assumed to be transferred to the corona. To avoid this assumption, we propose that gas in a vertically extended distribution is supplied to a supermassive black hole by the gravitational capture of interstellar medium or stellar wind material. In this picture, the gas partially condenses to an underlying cool disk as it flows toward the black hole, releasing accretion energy as X-ray emission and supplying mass for the disk accretion. Detailed numerical calculations reveal that the X-ray luminosity can reach a few tens of percent of the bolometric luminosity. The value of $\alpha_{\rm ox}$ varies from 0.9 to 1.2 for the mass supply rate ranging from 0.03 to 0.1 times the Eddington value. The corresponding photon index in the 2-10 keV energy band varies from 1.9 to 2.3. Such a picture provides a natural extension of the model for low luminosity AGN where condensation is absent at low mass accretion rates and no optically thick disk exists in the inner region.

We present the variability study of broad absorption lines (BALs) in a uniformly radio-selected sample of 28 BAL quasars using the archival data from the FIRST Bright Quasar Survey (FBQS) and the Sloan Digital Sky Survey (SDSS), as well as those obtained by ourselves, covering time scales $\sim 1-10$ years in the quasar's rest-frame. To our surprise, 5 quasars showing strong variations are all belong to a special subclass of overlapping iron low ionization BAL (OFeLoBAL) quasars, however, other 4 non-overlapping FeLoBALs (non-OFeLoBALs) are invariable except one case with weak optical depth change. Meanwhile, we also find 6 typical variations of high-ionization and low-ionization BALs in this BAL quasar sample. Photoionization models suggest that OFeLoBALs are formed in a relative dense ($n_e>10^6$ cm$^{-3}$) outflows at a distance from the subparsec to the dozens of parsecs from the continuum source. They differ from those of non-OFeLoBALs, which are likely produced by low-density gas, locating at a distance of hundreds to thousands parsecs. Thus, OFeLoBALs and non-OFeLoBALs, i.e., FeLoBALs with/without strong BAL variations, perhaps represent the bimodality of Fe II absorption, the former is located in the active galactic nucleus environment rather than the host galaxy. We suggest that high density and small distance are the necessary conditions that cause OFeLoBALs. As suggested in the literature, strong BAL variability is possibly due to variability of the covering factor of BAL regions caused by clouds transiting across the line of sight rather than ionization variations.

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

The observed tight radio/X-ray correlation in the low spectral state of some black hole X-ray binaries implies the strong coupling of the accretion and jet. The correlation of $L_{\rm R} \propto L_{\rm X}^{\sim 0.5-0.7}$ was well explained by the coupling of a radiatively inefficient accretion flow and a jet. Recently, however, a growing number of sources show more complicated radio/X-ray correlations, e.g., $L_{\rm R} \propto L_{\rm X}^{\sim 1.4}$ for $L_{\rm X}/L_{\rm Edd} \gtrsim 10^{-3}$, which is suggested to be explained by the coupling of a radiatively efficient accretion flow and a jet. In this work, we interpret the deviation from the initial radio/X-ray correlation for $L_{\rm X}/L_{\rm Edd} \gtrsim 10^{-3}$ with a detailed disc corona-jet model. In this model, the disc and corona are radiatively and dynamically coupled. Assuming a fraction of the matter in the accretion flow, $\eta\equiv \dot M_{\rm jet}/\dot M$, is ejected to form the jet, we can calculate the emergent spectrum of the disc corona-jet system. We calculate $L_{\rm R}$ and $L_{\rm X}$ at different $\dot M$, adjusting $\eta$ to fit the observed radio/X-ray correlation of the black hole X-ray transient H1743-322 for $L_{\rm X}/L_{\rm Edd}> 10^{-3}$. It is found that always the X-ray emission is dominated by the disc corona and the radio emission is dominated by the jet. We noted that the value of $\eta$ for the deviated radio/X-ray correlation for $L_{\rm X}/L_{\rm Edd} > 10^{-3}$, is systematically less than that of the case for $L_{\rm X}/L_{\rm Edd} < 10^{-3}$, which is consistent with the general idea that the jet is often relatively suppressed at the high luminosity phase in black hole X-ray binaries.

The Lidov-Kozai (LK) mechanism plays an important role in the secular evolution of many hierarchical triple systems. The standard LK mechanism consists of large-amplitude oscillations in eccentricity and inclination of a binary subject to the quadrupole potential from an outer perturber. Recent work has shown that when the octupole terms are included in the potential, the inner binary can reach more extreme eccentricities as well as undergo orientation flips. It is known that pericenter precessions due to short-range effects, such as General Relativity and tidal and rotational distortions, can limit the growth of eccentricity and even suppress standard (quadrupolar) LK oscillations, but their effect on the octupole-level LK mechanism has not been fully explored. In this paper, we systematically study how these short-range forces affect the extreme orbital behaviour found in octupole LK cycles. In general, the influence of the octupole potential is confined to a range of initial mutual inclinations Itot centered around 90deg (when the inner binary mass ratio is <<1), with this range expanding with increasing octupole strength. We find that, while the short-range forces do not change the width and location of this "window of influence", they impose a strict upper limit on the maximum achievable eccentricity. This limiting eccentricity can be calculated analytically, and its value holds even for strong octupole potential and for the general case of three comparable masses. Short-range forces also affect orbital flips, progressively reducing the range of Itot within which flips are possible as the intensity of these forces increases.

According to the core-accretion scenario, planets form in protostellar disks through the condensation of dust, coagulation of planetesimals, and emergence of protoplanetary embryos. At a few AU in a minimum mass nebula, embryos' growth is quenched by dynamical isolation due to the depletion of planetesimals in their feeding zone. However, embryos with masses ($M_p$) in the range of a few Earth masses ($M_\oplus$) migrate toward a transition radius between the inner viscously heated and outer irradiated regions of their natal disk. Their limiting isolation mass increases with the planetesimals surface density. When $M_p > 10 M_\oplus$, embryos efficiently accrete gas and evolve into cores of gas giants. We use numerical simulation to show that, despite streamline interference, convergent embryos essentially retain the strength of non-interacting embryos' Lindblad and corotation torque by their natal disks. In disks with modest surface density (or equivalently accretion rates), embryos capture each other in their mutual mean motion resonances and form a convoy of super Earths. In more massive disks, they could overcome these resonant barriers to undergo repeated close encounters including cohesive collisions which enable the formation of massive cores.

Nearly $15-20%$ of solar type stars contain one or more gas giant planet. According to the core-accretion scenario, the acquisition of their gaseous envelope must be preceded by the formation of super-critical cores with masses ten times or larger than that of the Earth. It is natural to link the formation probability of gas giant planets with the supply of gas and solid in their natal disks. However, a much richer population of super Earths suggests that 1) there is no shortage of planetary building-block material, 2) gas giants' growth barrier is probably associated with whether they can merge into super-critical cores, and 3) super Earths are probably failed cores which did not attain sufficient mass to initiate efficient accretion of gas before it is severely depleted. Here we construct a model based on the hypothesis that protoplanetary embryos migrated extensively before they were assembled into bona fide planets. We construct a Hermite-Embryo code based on a unified viscous-irradiation disk model and a prescription for the embryo-disk tidal interaction. This code is used to simulate 1) the convergent migration of embryos, and 2) their close encounters and coagulation. Around the progenitors of solar-type stars, the progenitor super-critical-mass cores of gas giant planets primarily form in protostellar disks with relatively high ($\gtrsim 10^{-7} M_\odot$ yr$^{-1}$) mass accretion rates whereas systems of super Earths (failed cores) are more likely to emerge out of natal disks with modest mass accretion rates, due to the mean motion resonance barrier and retention efficiency.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_\nu \propto \nu^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

Observations show that the accretion flows in low-luminosity active galactic nuclei (LLAGNs) probably have a two-component structure with an inner ADAF and an outer truncated accretion disk. As shown by Taam et al. (2012), the truncation radius as a function of mass accretion rate is strongly affected by including the magnetic field within the framework of disk evaporation model, i.e., an increase of the magnetic field results in a smaller truncation radius of the accretion disk. In this work, we calculate the emergent spectrum of an inner ADAF + an outer truncated accretion disk around a supermassive black hole based on the prediction by Taam et al. (2012). It is found that an increase of the magnetic field from $\beta=0.8$ to $\beta=0.5$ (with magnetic pressure $p_{\rm m}=B^2/{8\pi}=(1-\beta)p_{\rm tot}$, $p_{\rm tot}=p_{\rm gas}+p_{\rm m}$) results in an increase of $\sim 8.7$ times of the luminosity from the truncated accretion disk. We found that the equipartition of gas pressure to magnetic pressure, i.e., $\beta=0.5$, failed to explain the observed anti-correlation between $L_{\rm 2-10 keV}/L_{\rm Edd}$ and the bolometric correction $\kappa_{\rm 2-10 keV}$ (with $\kappa_{\rm 2-10 keV} = L_{\rm bol}/L_{\rm 2-10 keV}$). The emergent spectra for larger value $\beta=0.8$ or $\beta=0.95$ can well explain the observed $L_{\rm 2-10 keV}/L_{\rm Edd}$-$\kappa_{\rm 2-10 keV}$ correlation. We argue that in the disk evaporation model, the electrons in the corona are assumed to be heated only by a transfer of energy from the ions to electrons via Coulomb collisions, which is reasonable for the accretion with a lower mass accretion rate. Coulomb heating is the dominated heating mechanism for the electrons only if the magnetic field is strongly sub-equipartition, which is roughly consistent with observations.

A Survey of Ionized Gas in the Galaxy, made with the Arecibo telescope (SIGGMA) uses the Arecibo L-band Feed Array (ALFA) to fully sample the Galactic plane (30 < l < 75 and -2 < b < 2; 175 < l < 207 and -2 < b < 1) observable with the telescope in radio recombination lines (RRLs). Processed data sets are being produced in the form of data cubes of 2 degree (along l) x 4 degree (along b) x 151 (number of channels), archived and made public. The 151 channels cover a velocity range of 600 km/s and the velocity resolution of the survey changes from 4.2 km/s to 5.1 km/s from the lowest frequency channel to the highest frequency channel, respectively.RRL maps with 3.4 arcmin resolution and line flux density sensitivity of 0.5 mJy will enable us to identify new HII regions, measure their electron temperatures, study the physics of photodissociation regions (PDRs) with carbon RRLs, and investigate the origin of the extended low density medium (ELDM). Twelve Hn{\alpha} lines fall within the 300 MHz bandpass of ALFA; they are resampled to a common velocity resolution to improve the signal-to-noise ratio (SN) by a factor of 3 or more and preserve the line width. SIGGMA will produce the most sensitive fully sampled RRL survey to date. Here we discuss the observing and data reduction techniques in detail. A test observation toward the HII region complex S255/S257 has detected Hn{\alpha} and Cn{\alpha} lines with SN>10.

The quasar 1156+295 (4C +29.45) is one of the targets in the Urumqi monitoring program which aimed to search for evidence of annual modulation in the timescales of intra-day variable (IDV) sources. The IDV observations of 1156+295 were carried out nearly monthly from October 2007 to October 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV of total flux density in most observing sessions with variability timescales of <= 1 day at 4.8 GHz. The estimated IDV timescales seem to follow an annual cycle that can be fitted with an anisotropic interstellar scintillation (ISS) model suggesting that a significant part of the flux density variations is due to ISS. The source underwent a dramatic flare in 2008. We studied the possible consequences of the flare on the IDV of 1156+295 by comparing the changes in its IDV characteristics with the evolution of the 43 GHz Very Long Baseline Array (VLBA) core size of the source. Conclusions: The quasar 1156+295 shows evidence for an annual modulation of its IDV timescales at 4.8 GHz, the ISS-induced IDV timescales and variability strength might be affected by the overall activity state of the source core. More frequent IDV and VLBI measurements are required to confirm the relation between the IDV appearance and the core-size evolution of the source.

The magnitude of the viscosity and magnetic field parameters in hot accretion flows is investigated in low luminosity active galactic nuclei (LLAGNs). Theoretical studies show that a geometrically thin, optically thick disk is truncated at mass accretion rates less than a critical value by mass evaporated vertically from the disk to the corona, with the truncated region replaced by an advection dominated accretion flow (ADAF). The critical accretion rate for such a truncation is a function of the viscosity and magnetic field. Observations of X-ray photon indices and spectral fits of a number of LLAGNs published in the literature provide an estimate of the critical rate of mass accretion and the truncation radius respectively. By comparing the observational results with theoretical predictions, the viscosity and magnetic field parameters in the hot accretion flow region are estimated. Specifically, the mass accretion rates inferred in different sources constrain the viscosity parameter, whereas the truncation radii of the disk, as inferred from spectral fits, further constrain the magnetic field parameter. It is found that the value of the viscosity parameter in the corona/ADAF ranges from 0.17 to 0.5, with values clustered about 0.2-0.3. Magnetic pressure is required by the relatively small truncation radii for some LLAGNs and is found to be as high as its equipartition value with the gas pressure. The inferred values of the viscosity parameter are in agreement with those obtained from the observations of non-stationary accretion in stellar mass black hole X-ray transients. This consistency provides support for the paradigm that a geometrically thin disk is truncated by means of a mass evaporation process from the disk to the corona at low mass accretion rates.

The ongoing characterization of hot Jupiters has motivated a variety of circulation models of their atmospheres. Such models must be integrated starting from an assumed initial state, which is typically taken to be a wind-free, rest state. Here, we investigate the sensitivity of hot-Jupiter atmospheric circulation models to initial conditions. We consider two classes of models--shallow-water models, which have proven successful at illuminating the dynamical mechanisms at play on these planets, and full three-dimensional models similar to those being explored in the literature. Models are initialized with zonal jets, and we explore a variety of different initial jet profiles. We demonstrate that, in both classes of models, the final, equilibrated state is independent of initial condition--as long as frictional drag near the bottom of the domain and/or interaction with a specified planetary interior are included so that the atmosphere can adjust angular momentum over time relative to the interior. When such mechanisms are included, otherwise identical models initialized with vastly different initial conditions all converge to the same statistical steady state. In some cases, the models exhibit modest time variability; this variability results in random fluctuations about the statistical steady state, but we emphasize that, even in these cases, the statistical steady state itself does not depend on initial conditions. Although the outcome of hot-Jupiter circulation models depend on details of the radiative forcing and frictional drag, aspects of which remain uncertain, we conclude that the specification of initial conditions is not a source of uncertainty, at least over the parameter range explored in most current models.

The hadron-quark phase transition in the core of massive neutron stars is studied with a newly constructed two-phase model. For nuclear matter, a nonlinear Walecka type model with general nucleon-meson and meson-meson couplings, recently calibrated by Steiner, Hemper and Fischer, is taken. For quark matter, a modified Polyakov-Nambu--Jona-Lasinio (mPNJL) model, which gives consistent results with lattice QCD data, is used. Most importantly, we introduce an isoscalar-vector interaction in the description of quark matter, and we study its influence on the hadron-quark phase transition in the interior of massive neutron stars. With the constraints of neutron star observations, our calculation shows that the isoscalar-vector interaction between quarks is indispensable if massive hybrids star exist in the universe, and its strength determines the onset density of quark matter, as well as the mass-radius relations of hybrid stars. Furthermore, as a connection with heavy-ion-collision experiments we give some discussions about the strength of isoscalar-vector interaction and its effect on the signals of hadron-quark phase transition in heavy-ion collisions, in the energy range of the NICA at JINR-Dubna and FAIR at GSI-Darmstadt facilities.

Observations show that there is a positive correlation between Eddington ratio $\lambda$ and hard X-ray index $\Gamma$ for $\lambda \gtrsim 0.01$, and there is an anti-correlation between $\lambda$ and $\Gamma$ for $\lambda \lesssim 0.01$ in black hole X-ray binaries (with $\lambda=L_{\rm bol}/L_{\rm Edd}$). In this work, we theoretically investigate the correlation between $\Gamma$ and $\lambda$ within the framework of disk-corona model. We improve the model by taking into account all cooling processes including synchrotron and self-Compton radiations in the corona, Comptonization of the soft photons from the underlying accretion disk, and the Bremsstrahlung radiations. Presuming that the coronal flow above the disk can reach up to 0.1 Eddington rate at the outer region, we calculate the structure of the two-phase accretion flows and the emergent spectra for accretion rates from 0.003 to 0.1. It is found that at accretion rates larger than $\backsim 0.01$ Eddington rate a fraction of coronal gas condenses into the disk and an inner disk can be sustained by condensation. In this case, the X-ray emission is dominated by the scattering the soft photon from the underlying disk in the corona. The emission from the inner disk and corona can produce the positive correlation between $\lambda$ and $\Gamma$. While at accretion rates lower than $\backsim 0.01$ Eddington accretion rate, the inner disk vanishes completely by evaporation, the accretion is dominated by ADAF, in which the X-ray emission is produced by the Comptonization of the synchrotron and bremsstrahlung photons of ADAF itself. The emission from ADAF can produce the anti-correlation between $\lambda$ and $\Gamma$.

The truncation of an optically thick, geometrically thin accretion disk is investigated in the context of low luminosity AGN (LLAGN). We generalize the disk evaporation model used in the interpretative framework of black hole X-ray binaries by including the effect of a magnetic field in accretion disks surrounding supermassive black holes. The critical transition mass accretion rate for which the disk is truncated is found to be insensitive to magnetic effects, but its inclusion leads to a smaller truncation radius in comparison to a model without its consideration. That is, a thin viscous disk is truncated for LLAGN at an Eddington ratio less than 0.03 for a standard viscosity parameter ($\alpha = 0.3$). An increase of the viscosity parameter results in a higher critical transition mass accretion rate and a correspondingly smaller truncation distance, the latter accentuated by greater magnetic energy densities in the disk. Based on these results, the truncation radii inferred from spectral fits of LLAGN published in the literature are consistent with the disk evaporation model. The infrared emission arising from the truncated geometrically thin accretion disks may be responsible for the red bump seen in such LLAGN.

We aim to investigate the relation between the long-term flux density and the position angle (PA) evolution of inner-jet in blazars. We have carried out the elliptic Gaussian model-fit to the `core' of 50 blazars from 15 GHz VLBA data, and analyzed the variability properties of three blazars from the model-fit results. Diverse correlations between the long-term peak flux density and the PA evolution of the major axis of the `core' have been found in $\sim$ 20% of the 50 sources. Of them, three typical blazars have been analyzed, which also show quasi-periodic flux variations of a few years (T). The correlation between the peak flux density and the PA of inner-jet is positive for S5~0716+714, and negative for S4~1807+698. The two sources cannot be explained with the ballistic jet models, the non-ballistic models have been analyzed to explain the two sub-luminal blazars. A correlation between the peak flux density and the PA (with a T/4 time lag) of inner-jet is found in [HB89]~1823+568, this correlation can be explained with a ballistic precession jet model. All the explanations are based mainly on the geometric beaming effect; physical flux density variations from the jet base would be considered for more complicated situations in future, which could account for the no or less significance of the correlation between the peak flux density and the PA of inner-jet in the majority blazars of our sample.

The geometry of the accretion flow around stellar mass and supermassive black holes depends on the accretion rate. Broad iron emission lines originating from the irradiation of cool matter can indicate that there is an inner disk below a hot coronal flow.These emission lines have been detected in X-ray binaries. Observations with the Chandra X-ray Observatory, XMM Newton and Suzaku have confirmed the presence of these emission lines also in a large fraction of Seyfert-1 active galactic nuclei (AGN). We investigate the accretion flow geometry for which broad iron emission lines can arise in hard and soft spectral state. We study an ADAF-type coronal flow, where the ions are viscously heated and electrons receive their heat only by collisions from the ions and are Compton cooled by photons from an underlying cool disk. For a strong mass flow in the disk and the resulting strong Compton cooling only a very weak coronal flow is possible. This limitation allows the formation of ADAF-type coronae above weak inner disks in the hard state, but almost rules them out in the soft state. The observed hard X-ray luminosity in the soft state, of up to 10% or more of the total flux, indicates that there is a heating process that directly accelerates the electrons. This might point to the action of magnetic flares of disk magnetic fields reaching into the corona. Such flares have also been proposed by observations of the spectra of X-ray black hole binaries without a thermal cut-off around 200 keV.

We aim to search for evidence of annual modulation in the time scales of the BL Lac object S5 0716+714. The intra-day variability (IDV) observations were carried out monthly from 2005 to 2009, with the Urumqi 25m radio telescope at 4.8 GHz. The source has shown prominent IDV as well as long-term flux variations. The IDV time scale does show evidence in favor of an annual modulation, suggesting that the IDV of 0716+714 is dominated by interstellar scintillation. The source underwent a strong outburst phase between mid-2008 and mid-2009; a second intense flare was observed in late 2009, but no correlation between the total flux density and the IDV time scale is found, implying that the flaring state of the source does not have serious implications for the general characteristics of its intra-day variability. However, we find that the inner-jet position angle is changing throughout the years, which could result in an annual modulation noise in the anisotropic ISS model fit. There is also an indication that the lowest IDV amplitudes (rms flux density) correspond to the slowest time scales of IDV, which would be consistent with an ISS origin of the IDV of 0716+714.

We investigate the accretion process in high-luminosity AGNs (HLAGNs) in the scenario of the disk evaporation model. Based on this model, the thin disk can extend down to the innermost stable circular orbit (ISCO) at accretion rates higher than $0.02\dot{M}_{\rm Edd}$; while the corona is weak since part of the coronal gas is cooled by strong inverse Compton scattering of the disk photons. This implies that the corona cannot produce as strong X-ray radiation as observed in HLAGNs with large Eddington ratio. In addition to the viscous heating, other heating to the corona is necessary to interpret HLAGN. In this paper, we assume that a part of accretion energy released in the disk is transported into the corona, heating up the electrons and thereby radiated away. We for the first time, compute the corona structure with additional heating, taking fully into account the mass supply to the corona and find that the corona could indeed survive at higher accretion rates and its radiation power increases. The spectra composed of bremsstrahlung and Compton radiation are also calculated. Our calculations show that the Compton dominated spectrum becomes harder with the increase of energy fraction ($f$) liberating in the corona, and the photon index for hard X-ray($2-10 \rm keV$) is $2.2 < \Gamma < 2.7 $. We discuss possible heating mechanisms for the corona. Combining the energy fraction transported to the corona with the accretion rate by magnetic heating, we find that the hard X-ray spectrum becomes steeper at larger accretion rate and the bolometric correction factor ($L_{\rm bol}/L_{\rm 2-10keV}$) increases with increasing accretion rate for $f<8/35$, which is roughly consistent with the observational results.

Recent observations reveal that a cool disk may survive in the innermost stable circular orbit (ISCO) for some black hole X-ray binaries in the canonical low/hard state. The spectrum is characterized by a power law with a photon index $\Gamma \sim 1.5-2.1$ in the range of 2-10 keV and a weak disk component with temperature of $\sim 0.2$ keV. In this work, We revisit the formation of such a cool, optically thick, geometrically thin disk in the most inner region of black hole X-ray binaries at the low/hard state within the context of disk accretion fed by condensation of hot corona. By taking into account the cooling process associated with both Compton and conductive processes in a corona, and the irradiation of the hot corona to the disk, we calculate the structure of the corona. For viscosity parameter $\alpha=0.2$, it's found that the inner disk can exist for accretion rate ranging from $\dot M \sim 0.006-0.03 \dot M_{\rm Edd}$, over which the electron temperatures of the corona are in the range of $1-5\times 10^9\ \rm K$ producing the hard X-ray emission. We calculate the emergent spectra of the inner disk and corona for different mass accretion rates. The effect of viscosity parameter $\alpha$ and albedo $a$ ($a$ is defined as the energy ratio of reflected radiation from the surface of the thin disk to incident radiation upon it from the corona) to the emergent spectra are also presented. Our model is used to explain the recent observations of GX 339-4 and Cyg X-1, in which the thin disk may exist at ISCO region in the low/hard state at luminosity around a few percent of $L_{\rm Edd}$. It's found that the observed maximal effective temperature of the thermal component and the hard X-ray photon index $\Gamma$ can be matched well by our model.

There are mutually contradictory views in the literature of the kinematics and structure of high-ionization line (e.g. CIV) emitting regions in active galactic nuclei (AGNs). Two kinds of broad emission line region (BELR) models have been proposed, outflow and gravitationally bound BELR, which are supported respectively by blueshift of the CIV line and reverberation mapping observations. To reconcile these two apparently different models, we present a detailed comparison study between the CIV and MgII lines using a sample of AGNs selected from the Sloan Digital Sky Survey. We find that the kinematics of the CIV region is different from that of MgII, which is thought to be controlled by gravity. A strong correlation is found between the blueshift and asymmetry of the CIV profile and the Eddington ratio. This provides strong observational support for the postulation that the outflow is driven by radiation pressure. In particular, we find robust evidence that the CIV line region is largely dominated by outflow at high Eddington ratios, while it is primarily gravitationally bounded at low Eddington ratios. Our results indicate that these two emitting regions coexist in most of AGNs. The emission strength from these two gases varies smoothly with Eddington ratio in opposite ways. This explanation naturally reconciles the apparently contradictory views proposed in previous studies. Finally, candidate models are discussed which can account for both, the enhancement of outflow emission and suppression of normal BEL, in AGN with high Eddington ratios.

The phase transition of hadronic to quark matter at high baryon and isospin density is analyzed. Relativistic mean field models are used to describe hadronic matter, and the MIT bag model is adopted for quark matter. The boundaries of the mixed phase and the related critical points for symmetric and asymmetric matter are obtained. Due to the different symmetry term in the two phases, isospin effects appear to be rather significant. With increasing isospin asymmetry the binodal transition line of the (T,\rho_B) diagram is lowered to a region accessible through heavy ion collisions in the energy range of the new planned facilities, e.g. the FAIR/NICA projects. Some observable effects are suggested, in particular an "Isospin Distillation" mechanism with a more isospin asymmetric quark phase, to be seen in charged meson yield ratios, and an onset of quark number scaling of the meson/baryon elliptic flows. The presented isospin effects on the mixed phase appear to be robust with respect to even large variations of the poorly known symmetry term at high baryon density in the hadron phase. The dependence of the results on a suitable treatment of isospin contributions in effective QCD Lagrangian approaches, at the level of explicit isovector parts and/or quark condensates, is finally discussed.

Thermal conduction between a cool accretion disk and a hot inner corona can result in either evaporation of the disk or condensation of the hot corona. At low mass accretion rates, evaporation dominates and can completely remove the inner disk. At higher mass accretion rates, condensation becomes more efficient in the very inner regions, so that part of the mass accretes via a weak (initially formed) inner disk which is separated from the outer disk by a fully evaporated region at mid radii. At still higher mass accretion rates, condensation dominates everywhere, so there is a continuous cool disk extending to the innermost stable circular orbit. We extend these calculations by including the effect of irradiation by the hot corona on the disk structure. The flux which is not reflected is reprocessed in the disk, adding to the intrinsic thermal emission from gravitational energy release. This increases the seed photons for Compton cooling of the hot corona, enhancing condensation of the hot flow and re-inforcing the residual inner disk rather than evaporating it. Our calculations confirm that a residual inner disk can co-exist with a hard, coronally dominated, spectrum over a range of $0.006<\dot m<0.016$ (for $\alpha=0.2$). This provides an explanation for the weak thermal component seen recently in the low/hard state of black hole X-ray binary systems.

We report the discovery of a luminous ultra-soft X-ray excess in a radio-loud narrow-line Seyfert1 galaxy, RXJ1633+4718, from archival ROSAT observations. The thermal temperature of this emission, when fitted with a blackbody, is as low as 32.5(+8.0,-6.0)eV. This is in remarkable contrast to the canonical temperatures of ~0.1-0.2keV found hitherto for the soft X-ray excess in active galactic nuclei (AGN), and is interestingly close to the maximum temperature predicted for a postulated accretion disc in this object. If this emission is indeed blackbody in nature, the derived luminosity [3.5(+3.3,-1.5)x10^(44)ergs/s] infers a compact emitting area with a size [~5x10^(12)cm or 0.33AU in radius] that is comparable to several times the Schwarzschild radius of a black hole at the mass estimated for this AGN (3x10^6Msun). In fact, this ultra-steep X-ray emission can be well fitted as the (Compton scattered) Wien tail of the multi-temperature blackbody emission from an optically thick accretion disc, whose parameters inferred (black hole mass and accretion rate) are in good agreement with independent estimates using optical emission line spectrum. We thus consider this feature as a signature of the long-sought X-ray radiation directly from a disc around a super-massive black hole, presenting observational evidence for a black hole accretion disc in AGN. Future observations with better data quality, together with improved independent measurements of the black hole mass, may constrain the spin of the black hole.

We investigate the spectral features of accretion flows composed of an outer cool, optically thick disk and inner hot, optically thin, advection dominated accretion flows (ADAF) within the framework of disk and corona with mass evaporation (Liu et al. 2002a). In this work, both the magnetic field and Compton scattering of soft photons from the disk by electrons in the corona are included to calculate the evaporation rates at different distances. The disk is truncated at the distance where the evaporation rate equals to the accretion rate ($\dot m_{\rm evap}(r_{\rm tr})=\dot m$). For a series of accretion rates, the corresponding truncation radii are calculated out, with which we are able to calculate the emergent spectra from the inner ADAF + outer disk + corona. At very low accretion rates, the spectra are similar to that of a pure ADAF because the disk is truncated at large distances. The disk component becomes important at high accretion rates since the truncation occurs at small distances. An anti-correlation between the Eddington ratio $\xi \equiv L_{\rm 0.5-25 \,keV}/L_{\rm Edd}$ and the hard X-ray photon index $\Gamma_{\rm 3-25 \,keV}$ at low/hard states is predicted by the model. Comparing the theoretical results with observations, we find that our model can reproduce the anti-correlation between the Eddington ratio $\xi$ and the hard X-ray photon index observed for the X-ray binary XTE J1118+480.

Soft and hard spectral states of X-ray transient sources reflect two modes of accretion, accretion via a geometrically thin, optically thick disk or an advection-dominated accretion flow (ADAF). The luminosity at transition between these two states seems to vary from source to source, or even for the same source during different outbursts, as observed for GX 339-4. We investigate how the existence of an inner weak disk in the hard state affects the transition luminosity. We evaluate the structure of the corona above an outer truncated disk and the resulting disk evaporation rate for different irradiation. In some cases, recent observations of X-ray transients indicate the presence of an inner cool disk during the hard state. Such a disk can remain during quiescence after the last outburst as long as the luminosity does not drop to very low values (10^-4 to 10^-3 of the Eddington luminosity). Consequently, as part of the matter accretes via the inner disk, the hard irradiation is reduced. The hard irradiation is further reduced, occulted and partly reflected by the inner disk. This leads to a hard-soft transition at a lower luminosity if an inner disk exists below the ADAF. This seems to be supported by observations for GX 339-4.

The disk corona evaporation model extensively developed for the interpretation of observational features of black hole X-ray binaries (BHXRBs) is applied to AGNs. Since the evaporation of gas in the disk can lead to its truncation for accretion rates less than a maximal evaporation rate, the model can naturally account for the soft spectrum in high luminosity AGNs and the hard spectrum in low luminosity AGNs. The existence of two different luminosity levels describing transitions from the soft to hard state and from the hard to soft state in BHXRBs, when applied to AGNs, suggests that AGNs can be in either spectral state within a range of luminosities. For example, at a viscosity parameter, \alpha, equal to 0.3, the Eddington ratio from the hard to soft transition and from the soft to hard transition occurs at 0.027 and 0.005 respectively. When the Eddington ratio of the AGN lies below the critical value corresponding to its evolutionary state, the disk is truncated. With decreasing Eddington ratios, the inner edge of the disk increases to greater distances from the black hole with a concomitant increase in the inner radius of the broad line region, $R_{BLR}$. The absence of an optically thick inner disk at low luminosities gives rise to region in the size of borad line-luminosity plane for which the relation $R_{BLR} \propto L^{1/2}$ inferred at high luminosities is excluded. As a result, a lower limit to the accretion rate is predicted for the observability of broad emission lines, if the broad line region is associated with an optically thick accretion disk. Thus, true Seyfert 2 galaxies may exist at very low accretion rates/luminosities. The differences between BHXRBs and AGNs in the framework of the disk corona model are discussed and possible modifications to the model are briefly suggested.

We compile a blue AGN sample from SDSS and investigate the ratio of hard X-ray to bolometric luminosity in dependence on Eddington ratio and black hole mass. Our sample comprises 240 radio-quiet Seyfert 1 galaxies and QSOs. We find that the fraction of hard X-ray luminosity (log$(L_{\rm 2-10 kev}/L_{\rm bol})$) decreases with the increase of Eddington ratio. We also find that the fraction of hard X-ray luminosity is independent on the black hole mass for the radio-quiet AGNs. The relation of log$(L_{\rm 2-10 kev}/L_{\rm bol})$ decreasing with increasing Eddington ratio indicates that X-ray bolometric correction is not a constant, from a larger sample supporting the results of Vasudevan & Fabian (2007). We interpret our results by the disk corona evaporation/condensation model (Meyer et al. \cite{me200}; Liu et al. 2002a; Liu et al. 2007). In the frame of this model, the Compton cooling becomes efficient in cooling of the corona at high accretion rate (in units of Eddington rate), leading to condensation of corona gas to the disk. Consequently, the relative strength of corona to the disk becomes weaker at higher Eddington ratio. Therefore, the fraction of hard X-ray emission to disk emission and hence to the bolometric emission is smaller at higher Eddington ratio. The independence of the fraction of hard X-ray luminosity on the mass of the black hole can also be explained by the disk corona model since the corona structure and luminosity (in units of Eddington luminosity) are independent on the mass of black holes.

A wealth of Galactic accreting X-ray binaries have been observed both in low/hard state and high/soft state. The transition between these two states was often detected. Observation shows that the transition luminosity between these two states is different for different sources, ranging from 1% to 4% of the Eddington luminosity. Even for the same source the transition luminosity at different outbursts is also different. The transition can occur from 0.0069 to 0.15 Eddington luminosity. To investigate the underlying physics, we study the influence of viscosity parameter $\alpha$ on the transition luminosity on the basis of the disk-corona model for black holes. We calculate the mass evaporation rate for a wide range of viscosity parameter, $0.1\le \alpha\le 0.9$. By fitting the numerical results, we obtain fitting formulae for both the transition accretion rate and the corresponding radius as a function of $\alpha$. We find that the transition luminosity is very sensitive to the value of $\alpha$, $L/L_{\rm Edd}\propto\alpha^{2.34}$. For $0.1\le\alpha\le 0.6$, the transition luminosity varies by two orders of magnitude, from 0.001 to 0.2 Eddington luminosity. Comparing with observations we find that the transition luminosity can be fitted by adjusting the value of $\alpha$, and the model determined values of $\alpha$ are mostly in the range of observationally inferred value. Meanwhile we investigate the truncation of the disk in the low/hard state for some luminous sources. Our results are roughly in agreement with the observations.

The formation of a cool disk in the innermost regions of black hole X-ray transient systems in the low hard state is investigated. Taking into account the combined cooling associated with the Compton and conductive energy transport processes in a corona, the radial structure of a disk is described for a range of mass accretion rates. The mass flow in an optically thick inner region can be maintained by the condensation of matter from a corona with the disk temperature and luminosity varying continuously as a function of the accretion rate. Although such a disk component can be present, the contribution of the optically thick disk component to the total luminosity can be small since the mass flow due to condensation in the optically thick disk underlying the corona can be significantly less than the mass flow rate in the corona. The model is applied to the observations of the low quiescent state of the black hole source GX 339-4 at luminosities of around $0.01 L_{Edd}$ and is able to explain the temperature of the thermal component at the observed luminosities. Since conductive cooling dominates Compton cooling at low mass accretion rates, the luminosity corresponding to the critical mass accretion rate above which a weak thermal disk component can be present in the low hard state is estimated to be as low as $0.001 L_{Edd}$.

The condensation of matter from a corona to a cool, optically thick inner disk is investigated for black hole X-ray transient systems in the low hard state. A description of a simple model for the exchange of energy and mass between corona and disk originating from thermal conduction is presented, taking into account the effect of Compton cooling of the corona by photons from the underlying disk. It is found that a weak, condensation-fed inner disk can be present in the low hard state of black hole transient systems for a range of luminosities which depend on the magnitude of the viscosity parameter. For $\alpha \sim 0.1-0.4$ an inner disk can exist for luminosities in the range $\sim 0.001- 0.02$ Eddington value. The model is applied to the X-ray observations of the black hole candidate sources GX 339-4 and Swift J1753.5-0127 in their low hard state. It is found that Compton cooling is important in the condensation process, leading to the maintenance of cool inner disks in both systems. As the results of the evaporation/condensation model are independent of the black hole mass, it is suggested that such inner cool disks may contribute to the optical and ultraviolet emission of low luminosity active galactic nuclei.

The disk-corona evaporation model naturally interprets many observational phenomena in black hole X-ray binaries, such as the truncation of an accretion disk and the spectral state transitions. On the other hand, magnetic field is known to play an important role in transporting angular momentum and producing viscosity in accretion flows. In this work, we explicitly take the magnetic field in the accretion disk corona into account and numerically calculate the coronal structure on the basis of our two-temperature evaporation code. We show that the magnetic field influences the coronal structure by its contribution to the pressure, energy and radiative cooling in the corona and by decreasing the vertical heat conduction. We found that the maximal evaporation rate keeps more or less constant ($\sim 0.03$ Eddington rate) while the strength of magnetic fields changes, but that the radius corresponding to the maximal evaporation rate decreases with increasing magnetic field. This predicts that the spectral state transition always occurs at a few percent of Eddington accretion rate, while the inner edge of thin disk can be at $\sim 100 R_{\rm S} $ or even less in the hard state before the transition to the soft state. These results alleviate the problem that previous evaporation models predict too large a truncation radius, and are in better agreement with the observational results of several black hole X-ray binaries, though discrepancies remain.

Context: Accretion onto galactic and supermassive black holes occurs in different modes, which are documented in hard and soft spectral states, commonly attributed to an advection-dominated flow (ADAF) inside a truncated disk and standard disk accretion, respectively. At the times of spectral transition an intermediate state is observed, for which the accretion flow pattern is still unclear. Aims: We analyze the geometry of the accretion flow when the mass flow rate in the disk decreases (soft/hard transition) and evaporation of gas into the coronal flow leads to disk truncation. Methods: We evaluate the physics of an advection-dominated flow affected by thermal conduction to a cool accretion disk underneath. We find re-condensation of gas from the ADAF into the underlying inner disk at distances from the black hole and at rates, that depend on the properties of the hot ADAF and vary with the mass accretion rate. This sustains an inner disk for longer than a viscous decay time after the spectral transition occurred, in accordance with the spectra that indicate cool gas in the neighborhood of the accreting black hole. The model allows us to understand why Cyg X-1 does not show hysteresis in the spectral state transition luminosity that is commonly observed for X-ray transient sources. Conclusions: Our results shed new light on the complex mass flow pattern during spectral state transition.

We combine recent long Gamma Ray Bursts (GRBs) sample including 52 objects out to z$=$6.3 compiled from Swift Gamma Ray Bursts by Schaefer (2006) with Type Ia Supernova (SNIa), Cosmic Microwave Background (CMB) and Baryon Oscillation (BAO) to constrain cosmologies. We study the constraints arising from GRBs alone and the complementarity of GRBs to other cosmological probes. To analyze the cosmological role of GRBs, we adopt the Hubble constant as a free parameter with a prior in the range of $H_0=72\pm8\enspace\kmsmpc$ instead of a fixed $H_0$ used in previous studies. By jointly using SNIa gold/SNLS samples and GRBs, the constraints on \Om -\Ol parameter space are dramatically improved in comparing with those from SNIa data alone. The complementarity of GRBs is mostly from the sub-sample with z$>$1.5 due to the different parameter degeneracies involved in luminosity distances at different redshifts. Including GRB data in addition to SNIa, BAO and CMB in our analysis, we find that the concordance model with \Om $\sim$0.3 and \Ol$\sim$0.7 is still well within 1$\sigma$ confidence range.

Aims: The hard and soft spectral states of black hole accretion are understood as connected with ADAF accretion (truncated disk) and standard disk accretion, respectively. However, observations indicate the existence of cool gas in the inner region at times when the disk is already truncated outside. We try to shed light on these not yet understood intermediate states. Methods: The disk-corona model allows to understand the spectral state transitions as caused by changes of the mass flow rate in the disk and provides a picture for the accretion geometry when disk truncation starts at the time of the soft/hard transition, the formation of a gap in the disk filled by an advection-dominated flow (ADAF) at the distance where the evaporation is maximal. We study the interaction of such an ADAF with an inner thin disk below. Results: We show that, when the accretion rate is not far below the transition rate, an inner disk could exist below an ADAF, leading to an intermediate state of black hole accretion.

In neutron star and black hole X-ray binaries the transitions between the two spectral states, hard and soft, signals the change between accretion via a hot advection-dominated flow(ADAF) and disk accretion. In a few cases the hard/soft transition were observed during the rise to the nova outburst, mostly only the soft/hard transition during the luminosity decrease. Surprisingly the luminosities at the second transition is always lower by a factor of 3 to 5. A model for this hysteresis was presented in a preceding paper: It was shown that this switch in the accretion mode at different mass accretion rates and therefore different luminosities is caused by the different amount of Compton cooling or heating of the accretion disk corona as it is irradiated by hard or soft radiation from the central light source, respectively. We now give detailed results on the dependence on hardness of this radiation and on radiation efficiency. We further discuss the influence of the inclination and a possible warping of the disk on the observed hysteresis.

Many low-mass X-ray binaries show both hard and soft spectral states. For several sources the transitions between these states have been observed, mostly from the soft to the hard state during a luminosity decrease. In a few cases also the transition from the hard to the soft state was observed, coincident with an increase of the luminosity. Surprisingly this luminosity was not the same as the one during a following change back to the hard state. The values differed by a factor of about 3 to 5. We present a model for this hysteresis in the light curves of low-mass X-ray binaries (sources with neutron stars or black holes). We show that the different amount of Compton cooling or heating acting on the accretion disk corona at the time of the transition causes this switch in the accretion mode at different mass accretion rates and therefore different luminosities. The inner disk during the soft state provides a certain amount of Compton cooling which is either not present or much less if the inner region is filled with a hot advection-dominated accretion flow (ADAF) that radiates a hard spectrum.

We raise the question whether in the past a disk could have existed in our Galactic Center which has disappeared now. Our model for the interaction of a cool disk and a hot corona above (Liu et al. 2004) allows to estimate an upper limit for the mass that might have been present in a putative accretion disk after a last star forming event, but would now have evaporated by coronal action.

We study the possibility of a cool disk existing in the Galactic Center in the framework of the disk-corona evaporation/condensation model. Assuming an inactive disk, a hot corona should form above the disk since there is a continuous supply of hot gas from stellar winds of the close-by massive stars. Whether the cool disk can survive depends on the mass exchange between the disk and corona. If the disk-corona interaction is dominated by evaporation and the rate is larger than the Bondi accretion rate in the Galactic Center, the disk will be depleted within a certain time period and no persistent disk will exist. On the other hand, if the interaction results in hot gas steadily condensing into the disk, an inactive cool disk might survive. For this case we further investigate the Bremsstrahlung radiation from the hot corona and compare it with the observed X-ray luminosity. Our model shows that, for standard viscosity in the corona (alpha=0.3), the mass evaporation rate is much higher than the Bondi accretion rate and the coronal density is much larger than that inferred from Chandra observations. An inactive disk can not survive such strong evaporation. For small viscosity (alpha<0.07) we find condensation solutions. But detailed computations show that in this case there is too much X-ray radiation from the corona to be in agreement with the observations. Therefore, we conclude that there should be no thin/inactive disk presently in the Galactic Center. However, we do not exclude that the alternative non-radiative model of Nayakshin (2004) might instead be realized in nature.

We investigate a corona coupled with underlying disk through magnetic field and radiation field, and present emergent spectra calculated from Monte Carlo simulations. The spectra are determined for given black-hole mass and accretion rate. We find two types of solutions corresponding for hard spectrum and soft spectrum. In the hard-spectrum solution, the accretion energy is dominantly dissipated in the corona, supporting a strong corona above a cool disk; The hard X-ray spectral indices are around 1.1, same for different accretion rates. In the soft-spectrum solution, the accretion energy is mainly dissipated in the disk. The coronal temperature and density are quite low. Consequently, the spectra are dominated by the disk radiation peaking at UV and soft X-rays. For low-luminosity systems, there exists only the solution of hard spectra; While for high-luminosity systems, there exist both solutions of hard and soft spectra. For middle-luminosity systems, besides the hard spectra, moderately soft spectra composed of an inner soft-spectrum solution and an outer hard-spectrum solution may occur, the softness of which increases with increasing luminosity. The hard spectra are close to the observed spectra in Seyfert galaxies and radio-quiet QSOs. The composite spectra may account for the diversity of broad band spectra observed in narrow-line Seyfert 1 galaxies.

We construct a simple model for a magnetic reconnection heated corona above a thin accretion disk in AGNs and Galactic black hole candidates(GBHCs). The basic assumptions are that (1) the magnetic reconnection heat is cooled down overwhelmingly by Compton scattering in the corona and that (2) thermal conduction is dominantly cooled by evaporation of the chromospheric plasma in the disk-corona interface before Compton cooling sets in. With these two basic equations as well as equi-partition of magnetic energy with gas energy in the disk, we can consistently determine the fraction of accretion energy dissipated in the corona without free parameters, and thus determine the temperature and all other quantities in both the corona and disk for given black hole mass and accretion rate. Then, we calculate the luminosity contributed from the disk and corona and the coronal flux weighted Compton $y$ parameter. It is found that, at a low luminosity (less than 0.1 Eddington luminosity) the spectrum is hard with energy spectral index of around 1, while at a high luminosity (larger than 0.1 Eddington luminosity) the spectrum can be either soft or hard. We also find that the situation is almost same for supermassive and stellar-mass black holes. These features are consistent with observations of AGNs and GBHCs.

We extended the disk corona model (Meyer & Meyer-Hofmeister 1994; Meyer, Liu, & Meyer-Hofmeister 2000a) to the inner region of galactic nuclei by including different temperatures in ions and electrons as well as Compton cooling. We found that the mass evaporation rate and hence the fraction of accretion energy released in the corona depend strongly on the rate of incoming mass flow from outer edge of the disk, a larger rate leading to more Compton cooling, less efficient evaporation and a weaker corona. We also found a strong dependence on the viscosity, higher viscosity leading to an enhanced mass flow in the corona and therefore more evaporation of gas from the disk below. If we take accretion rates in units of the Eddington rate our results become independent on the mass of the central black hole. The model predicts weaker contributions to the hard X-rays for objects with higher accretion rate like narrow-line Seyfert 1 galaxies (NLS1s), in agreement with observations. For luminous active galactic nuclei (AGN) strong Compton cooling in the innermost corona is so efficient that a large amount of additional heating is required to maintain the corona above the thin disk.

The concept of an advection-dominated accretion flow (ADAF), with or without wind loss, was used to describe the spectra of the galactic center source Sgr A*, low-luminosity AGN and nuclei of elliptical galaxies (including M87). The spectral fits of various authors show that the transition from the geometrically thin disk to the hot flow occurs at quite different distances, apparently not uniquely related to the mass flow rate in the disk. We compare these results with the results of theoretical modeling where we determine the truncation of the thin outer disk from the efficiency of mass evaporation. The physics is the same as in the case of galactic black holes systems (Meyer et al. 2000b). For the observationally indicated mass flow rates our model predicts a truncation at 10^3 to 10^4 Schwarzschild radii. We discuss whether far inside this truncation an innermost cool thin disk could exist and affect the spectrum.

We present a model for a corona above a geometrically thin standard disk around a black hole. This corona is fed by matter from the thin disk which evaporates from the cool layers underneath. An equilibrium establishes between the cool accretion stream and the hot flow. In the inner region evaporation becomes so efficient that at low accretion rates all matter flows via the corona and proceeds towards the black hole as an advection-dominated accretion flow (ADAF). We investigate this transition for accretion disks around stellar black holes. We derive a set of four ordinary differential equations which describe the relevant processes and solve them numerically. The evaporation efficiency has a maximum value where the advection-dominated structure of the corona (at large distances from the black hole) changes to a radiation-dominated structure (at small distances). This maximum can be related to the observed spectral transitions in black hole X-ray binaries. Our solutions of coronal structure are valid for both, accretion in galactic X-ray transients and in AGN.

The theory of coronal evaporation predicts the formation of an inner hole in the cool thin accretion disk for mass accretion rates below a certain value (~ 1/50 of the Eddington mass accretion rate) and the sudden disappearance of this hole when the mass accretion rate rises above that value. The inner edge of the standard thin disk then suddenly shifts inward from about a few hundred Schwarzschild radii to the last stable orbit. This appears to quantitatively account for the observed transitions between hard and soft spectral states at critical luminosities. Due to the evaporation process the matter accreting in the geometrically thin disk changes to a hot coronal flow which proceeds towards the black hole as an advection-dominated accretion flow (ADAF; for a review see Narayan et al. 1998).

We apply the disk-corona evaporation model (Meyer & Meyer-Hofmeister) originally derived for dwarf novae to black hole systems. This model describes the transition of a thin cool outer disk to a hot coronal flow. The mass accretion rate determines the location of this transition. For a number of well studied black hole binaries we take the mass flow rates derived from a fit of the ADAF model to the observed spectra (for a review see Narayan, Mahadevan, & Quataert) and determine where the transition of accretion via a cool disk to a coronal flow/ADAF would be located for these rates. We compare with the observed location of the inner disk edge, as estimated from the maximum velocity of the $\rm H_\alpha$ emission line. We find that the transition caused by evaporation agrees with this in stellar disks. We also show that the ADAF and the ``thin outer disk + corona'' are compatible in terms of the physics in the transition region.

We model the evolution of the accretion disk of WZ Sagittae during the long quiescence. We find that the large amount of mass in the disk derived from the outburst luminosity is a severe constraint and demands values of $\alpha_c$ $\approx$ 0.001 in contradiction to some recent suggestions. We include in our computations the formation of an inner disk hole and the growth of the disk due to redistribution of angular momentum. We find a new mode of disk evolution. The disk is quasi-stationary. Only about half of the mass transfered from the companion star flows through the disk, the other half is needed to build up the steadily growing outer disk. When the 3:1 resonance radius is reached the disk growth ends. From then on all transferred matter flows inward, the surface density increases, leading to an outburst within a few years. We predict superhumps at low luminosity during this last phase. We discuss X-rays expected, the white dwarf mass and distance to WZ Sagittae.

All data available in B band for the BL Lac object MRK421 from 22 publications are used to construct a historical light curve, dating back to 1900. It is found that the light curve is very complicated and consists of a set of outbursts with very large duration. The brightness of MRK421 varies from 11.6 magnitude to more than 16 magnitude. Analyses with Jurkevich method of computing period of cyclic phenomena reveal in the light curve two kinds of behaviors. The first one is non-periodic with rapid, violent variations in intensity on time scales of hours to days. The second one is periodic with a possible period of $23.1\pm 1.1$ years. Another possible period of $15.3\pm 0.7$ years is not very significant. We have tested the robustness of the Jurkevich method. The period of about one year found in the light curves of MRK421 and of other objects is a spurious period due to the method and the observing window. We try to explain the period of $23.1 \pm1.1$ years under the thermal instability of a slim accretion disk around a massive black hole of mass of $2 *10^6 M_\odot$.