We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $\tau(\lambda)$, generally follow the relation $\tau(\lambda) \propto \lambda^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $\tau(\lambda) \propto \lambda^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_\nu\propto\nu^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

We investigate whether the hard X-ray photon index (${\Gamma}$) versus accretion rate correlation for super-Eddington accreting quasars is different from that for sub-Eddington accreting quasars. We construct a sample of 113 bright quasars from the Sloan Digital Sky Survey Data Release 14 quasar catalog, including 38 quasars as the super-Eddington subsample and 75 quasars as the sub-Eddington subsample. We derive black-hole masses using a simple-epoch virial mass formula based on the ${\rm H\beta}$ lines, and we use the standard thin disk model to derive the dimensionless accretion rates ($\dot{\mathscr{M}}$) for our sample. The X-ray data for these quasars are collected from the Chandra and XMM-Newton archives. We fit the hard X-ray spectra using a single power-law model to obtain ${\Gamma}$ values. We find a statistically significant ($R_{\rm S}=0.43$, $p=7.75\times{10}^{-3}$) correlation between ${\Gamma}$ and $\dot{\mathscr{M}}$ for the super-Eddington subsample. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlation for the sub-Eddington subsample is also significant, but weaker ($R_{\rm S}=0.30$, $p=9.98\times{10}^{-3}$). Linear regression analysis shows that ${\rm \Gamma}=(0.34\pm0.11){\rm log}{\dot{\mathscr{M}}}+(1.71\pm0.17)$ and ${\rm \Gamma}=(0.09\pm0.04){\rm log}{\dot{\mathscr{M}}}+(1.93\pm0.04)$ for the super- and sub-Eddington subsamples, respectively. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlations of the two subsamples are different, suggesting different disk-corona connections in these two types of systems. We propose one qualitative explanation of the steeper ${\Gamma}$-$\dot{\mathscr{M}}$ correlation in the super-Eddington regime that involves larger seed photon fluxes received by the compact coronae from the thick disks in super-Eddington accreting quasars.

We perform a-priori validation tests of subgrid-scale (SGS) models for the turbulent transport of momentum, energy and passive scalars. To this end, we conduct two sets of high-resolution hydrodynamical simulations with a Lagrangian code: an isothermal turbulent box with rms Mach number of 0.3, 2 and 8, and the classical wind tunnel where a cold cloud traveling through a hot medium gradually dissolves due to fluid instabilities. Two SGS models are examined: the eddy diffusivity (ED) model wildly adopted in astrophysical simulations and the "gradient model" due to Clark et al. (1979). We find that both models predict the magnitude of the SGS terms equally well (correlation coefficient > 0.8). However, the gradient model provides excellent predictions on the orientation and shape of the SGS terms while the ED model predicts poorly on both, indicating that isotropic diffusion is a poor approximation of the instantaneous turbulent transport. The best-fit coefficient of the gradient model is in the range of [0.16, 0.21] for the momentum transport, and the turbulent Schmidt number and Prandtl number are both close to unity, in the range of [0.92, 1.15].

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

PSR J2021+4026 is a radio-quiet gamma-ray pulsar and the first pulsar that shows state change of the gamma-ray emission and spin-down rate. The state change of PSR J2021+4026 was first observed at 2011 October, at which the pulsar changes the state from high gamma-ray flux/low spin-down rate state to low gamma-ray flux/high spin-down rate st\ ate. In December 2014, PSR J2021+4026 recovered the state before the 2011 state change over a timescale of a few months. We report that the long term evolution of the gamma-ray flux and timing behavior suggests that PSR J2021+4026 changed the state near 2018 February 1st and entered a new low gamma-ray flux/high spin-down rate state. At the 2018 state change, the averaged flux dropped from $(1.29\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1}}$ to $(1.12\pm 0.01)\times 10^{-6} {\rm cts~cm^{-2}s^{-1\ }}$, which has the similar behavior to the case of 2011 event. The spin-down rate has increased by $\sim 3\%$ in the new state since the 2018 state change. The shapes of pulse profile and spectrum in GeV bands also changed at the 2018 event, and they are consistent with behavior at the 2011 state change. Our results probably suggest that PSR J2021+4026 is switching between different states with a timescale of several years, like some radio pulsars (e.g. PSR~B1828-11). PSR J2021+4026 will provide a unique opportunity to study the mechanism of the state switching.

We investigate the nature of the unidentified very-high-energy (VHE) gamma-ray object, HESS J1832-093, in a multi-wavelength context. Based on X-ray variability and spectral index ($\Gamma_X\sim\,1.5$), and its broad-band spectrum (which was remarkably similar to HESS J0632+057, a confirmed gamma-ray binary), HESS J1832-093 has been considered to be a strong gamma-ray binary candidate in previous works. In this work, we provide further evidence for this scenario. We obtained a spectrum of its IR counterpart using Gemini/Flamingo, finding absorption lines that are usually seen in massive stars, in particular O stars. We also obtained a rather steep ATCA spectrum ($\alpha=-1.18^{+1.04}_{-0.88}$) which prefers a gamma-ray binary over an AGN scenario. Based on spatial-spectral analysis and variability search, we found that 4FGL J1832.9-0913 is possible to be associated with SNR G22.7-0.2 rather than with HESS J1832-093 only.

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$\beta$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$\beta$, He II, He I, and Fe II lines with respect to the varying 5100\AA continuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

The long-term Swift monitoring of ESO 243-49 HLX-1 provides an opportunity to investigate the detailed timing and spectral behaviour of this hyper-luminous X-ray source. Swift has detected 7 outbursts since 2009 mid-August. Using different dynamical timing algorithms, we confirm an increasing trend for the time intervals between outbursts, which is manifest in the delays between the latest outbursts. The X-ray spectra of HLX-1 in quiescence can be described with a single power-law model while the thermal component dominates the X-ray emission during outburst. There is only marginal evidence for photon index (or spectral hardness) changes between quiescent states with about 1{\sigma} deviation. With the updated temporal and spectral features, we re-examine different scenarios to explain the origin of the quasi-periodic modulation of HLX-1. A significantly increasing trend without obvious stochastic fluctuations on the timescale of the detected quasi-period may not fully support an orbital period origin as might be due to mass transfer episodes from a donor star at periastron of an extremely eccentric orbit. The outburst profile seems to be consistent with the effect of tidal-induced-precession of an accretion disc or an oscillating wind scenario in the inner disc. Based on these models, we speculate that the true orbital period is much shorter than the detected quasi-periodicity.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

The X-ray pulsar SMC X-1 shows a superorbital modulation with an unstable cycle length in the X-ray band. We present its timing behaviors, including the spin, orbital, and superorbital modulations, beyond the end of the Rossi X-ray Timing Explorer mission. The time-frequency maps derived by the wavelet Z-transform and the Hilbert-Huang transform suggest that a new superorbital period excursion event occurred in ~MJD 57100 (2015 March). This indicates the excursion is recurrent and probably (quasi)periodic. The hardness ratio obtained with the Monitor of All-sky X-ray Image (MAXI) suggests increased absorption during the transition from the high to the low state in the superorbital cycle. Compared to the regular epochs, the superorbital profile during the excursion epochs has a shallower and narrower valley, likely caused by a flatter warp. By tracking the spin period evolution with the MAXI Gas Slit Camera in 2-20 keV, we derive an averaged spin-up rate of $\dot{\nu}=2.515(3)\times10^{-11}$ s$^{-2}$ during the period between MJD 55141 (2009 November) and 58526 (2019 February). We obtain no positive correlation between the spin frequency residual and the superorbital frequency, but a torque change accompanying the superorbital period excursion is possible. We suggest that the accretion torque on the neutron star could be changed by various mechanisms, including the change of mass accretion rate and the warp angle. We update the value of the orbital decay as $\dot{P}_{\rm{orb}}/P_{\rm{orb}}=-3.380(6)\times10^{-6}$ yr$^{-1}$. Finally, we reconfirm the detection of the superorbital modulation in the optical band and its coherence in phase with the X-ray modulation.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

We present the revised and high quality CCD time series observations in 2017 and 2018 of the W UMa-type contact binary BH~Cas. The combination of our new determinations of minimum times and those collected from the literatures reveal a long-term period increase with a rate of $\textit{dP/dt}\,=\,+3.27 \times 10^{-7} \textrm{d}\,\textrm{yr}^{-1}$, which can be accounted for by angular momentum transfer as the less-massive component loses its mass to the more massive one. One cyclical oscillation ($A'=0.00300$~days and $P'=20.09$~yrs) superimposed on the long-term period increase tendency is attributed to one additional object that orbits BH Cas. The low-resolution spectra of this source are observed at phase $\sim$0.0 and $\sim$0.5 in 2018 and 2019, and the spectral types are identified as K3V$\,\pm\,$1 and G8V$\,\pm\,$2. According to the luminosity contributions and spectra at these phases, the spectral type of the primary star is close to K3V$\,\pm\,$1. We infer a cool spot on the hotter and less-massive component to fit the asymmetry of light curves. In addition, changes in the location, temperature, and area of spots on the secondary star in different observations may indicate that the magnetic field activity of the secondary is more active than that of the primary. Such strong spot activities and optical flare are both detected at the same time from a W UMa-type star for the first time. With the \emph{XMM-Newton} observations taken in 2003 and 2012, the X-ray light curve is obtained and the X-ray spectra can be well described with a black-body or a bremsstrahlung model. The soft X-ray luminosity of BH Cas is estimated as 8.8\,--\,9.7$\,\times\,10^{30}$\,erg\,s$^{-1}$.

We present a hydro-dynamical simulation at sub-parsec and few solar mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium (ISM) and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of $\sim 900$ stellar clusters with masses above $10^{2.5}\; \rm{M_\odot}$ and a cluster mass function (CMF), which is well fitted with a power-law with a slope of $\alpha=-1.70\pm0.08$. We describe here in detail the formation of the three most massive clusters (M$_{*} \gtrsim 10^5$ M$_\odot$), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of $6$-$8$ Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of $\Sigma_\mathrm{SFR}>10\; \mathrm{M}_\odot\; \mathrm{yr}^{-1}\; \mathrm{kpc}^{-2}$ and thermal gas pressures in excess of $P_{\rm th}\sim10^7 \; \mathrm{k}_{\rm B}\;(\rm K\;\mathrm{cm}^{-3})^{-1}$. The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters. The simulation presented here provides a general model for the formation of metal-poor globular clusters in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts.

We study the destruction of interstellar dust via sputtering in supernova (SN) shocks using three-dimensional hydrodynamical simulations. With a novel numerical framework, we follow both sputtering and dust dynamics governed by direct collisions, plasma drag and betatron acceleration. Grain-grain collisions are not included and the grain-size distribution is assumed to be fixed. The amount of dust destroyed per SN is quantified for a broad range of ambient densities and fitting formulae are provided. Integrated over the grain-size distribution, nonthermal (inertial) sputtering dominates over thermal sputtering for typical ambient densities. We present the first simulations that explicitly follow dust sputtering within a turbulent multiphase interstellar medium. We find that the dust destruction timescales $\tau$ are 0.35 Gyr for silicate dust and 0.44 Gyr for carbon dust in solar neighborhood conditions. The SN environment has an important impact on $\tau$. SNe that occur in preexisting bubbles destroy less dust as the destruction is limited by the amount of dust in the shocked gas. This makes $\tau$ about 2.5 times longer than the estimate based on results from a single SN explosion. We investigate the evolution of the dust-to-gas mass ratio (DGR), and find that a spatial inhomogeneity of $\sim$ 14\% develops for scales below 10 pc. It locally correlates positively with gas density but negatively with gas temperature even in the exterior of the bubbles due to incomplete gas mixing. This leads to a $\sim$ 30\% lower DGR in the volume filling warm gas compared to that in the dense clouds.

We compiled a sample of 73 active galactic nuclei (AGNs) with reverberation mapping (RM) observations from RM campaigns including our ongoing campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs). This sample covers a large range of black hole (BH) mass $(M_{\bullet}=10^{6-9}~M_\odot)$, dimensionless accretion rates $(\dot{\mathscr{M}}=10^{-2.7}-10^{2.7})$ and 5100~\AA~luminosity $(L_{5100}=10^{42-46}~\rm erg~s^{-1})$, allowing us to systematically study the AGN variability and their relations with BH mass, accretion rates, and optical luminosity. We employed the damped random walk (DRW) model to delineate the optical variability of continuum at 5100~\AA~and obtained damped variability timescale ($\tau_{\rm d}$) and amplitude ($\sigma_{\rm d}$) using a Markov Chain Monte Carlo (MCMC) method. We also estimated the traditional variability amplitudes ($F_{\rm var}$), which provide a model-independent measure and therefore are used to test the DRW results. We found that AGN variability characteristics are generally correlated with $(M_{\bullet},\dot{\mathscr{M}},L_{5100})$. These correlations are smooth from sub-Eddington to super-Eddington accretion AGNs, probably implying that the AGN variability may be caused by the same physical mechanism.

We report the discovery of extreme X-ray variability in a type 1 quasar: SDSS J$075101.42+291419.1$. It has a black hole mass of $1.6\times 10^7~\rm M_\odot$ measured from reverberation mapping (RM), and the black hole is accreting with a super-Eddington accretion rate. Its XMM-Newton observation in 2015 May reveals a flux drop by a factor of $\sim 22$ with respect to the Swift observation in 2013 May when it showed a typical level of X-ray emission relative to its UV/optical emission. The lack of correlated UV variability results in a steep X-ray-to-optical power-law slope ($\alpha_{\rm OX}$) of -1.97 in the low X-ray flux state, corresponding to an X-ray weakness factor of 36.2 at rest-frame 2 keV relative to its UV/optical luminosity. The mild UV/optical continuum and emission-line variability also suggest that the accretion rate did not change significantly. A single power-law model modified by Galactic absorption describes well the $0.3-10$ keV spectra of the X-ray observations in general. The spectral fitting reveals steep spectral shapes with $\Gamma\approx3$. We search for active galactic nuclei (AGNs) with such extreme X-ray variability in the literature and find that most of them are narrow-line Seyfert 1 galaxies and quasars with high accretion rates. The fraction of extremely X-ray variable objects among super-Eddington accreting AGNs is estimated to be $\approx 15-24\%$. We discuss two possible scenarios, disk reflection and partial covering absorption, to explain the extreme X-ray variability of SDSS J$075101.42+291419.1$. We propose a possible origin for the partial covering absorber, which is the thick inner accretion disk and its associated outflow in AGNs with high accretion rates.

We report results of the first reverberation mapping campaign of I Zwicky 1 during $2014$-$2016$, which showed unambiguous reverberations of the broad H$\beta$ line emission to the varying optical continuum. From analysis using several methods, we obtain a reverberation lag of $\tau_{\rm H\beta}=37.2^{+4.5}_{-4.9}\,$ days. Taking a virial factor of $f_{_{\rm BLR}}=1$, we find a black hole mass of $M_{\bullet}=9.30_{-1.38}^{+1.26}\times 10^6 M_{\odot}$ from the mean spectra. The accretion rate is estimated to be $203.9_{-65.8}^{+61.0}\,L_{\rm Edd}c^{-2}$, suggesting a super-Eddington accretor, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. By decomposing {\it Hubble Space Telescope} images, we find that the stellar mass of the bulge of its host galaxy is $\log (M_{\rm bulge}/M_{\odot}) = \rm 10.92\pm 0.07$. This leads to a black hole to bulge mass ratio of $\sim 10^{-4}$, which is significantly smaller than that of classical bulges and elliptical galaxies. After subtracting the host contamination from the observed luminosity, we find that I Zw 1 follows the empirical $R_{\rm BLR}\propto L_{5100}^{1/2}$ relation.

Magnetars are strongly magnetized pulsars and they occasionally show violent radiative outbursts. They also often exhibit glitches which are sudden changes in the spin frequency. It was found that some glitches were associated with outbursts but their connection remains unclear. We present a systematic study to identify possible correlations between them. We find that the glitch size of magnetars likely shows a bimodal distribution, different from the distribution of the Vela-like recurrent glitches but consistent with the high end of that of normal pulsars. A glitch is likely a necessary condition for an outburst but not a sufficient condition because only 30\% of glitches were associated with outbursts. In the outburst cases, the glitches tend to induce larger frequency changes compared to those unassociated ones. We argue that a larger glitch is more likely to trigger the outburst mechanism, either by reconfiguration of the magnetosphere or deformation of the crust. A more frequent and deeper monitoring of magnetars is necessary for further investigation of their connection.

Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$\beta$ relative to those of the 5100 Angstrom continuum is $146.8_{-12.1}^{+8.3}$ days in the rest frame, which agrees very well with the Paschen-$\alpha$ region measured by the GRAVITY at The Very Large Telescope Interferometer. The time lag of the H$\gamma$ emission line is found to be nearly the same as for H$\beta$. The lag of the Fe II emission is $322.0_{-57.9}^{+55.5}$ days, longer by a factor of $\sim$2 than that of the Balmer lines. The velocity-resolved lag measurements of the H$\beta$ line show a complex structure which can be possibly explained by a rotation-dominated disk with some inflowing radial velocity in the H$\beta$-emitting region. Taking the virial factor of $f_{\rm BLR} = 1.3$, we derive a BH mass of $M_{\bullet} = 4.1_{-0.4}^{+0.3} \times 10^8 M_{\odot}$ and an accretion rate of $9.3\,L_{\rm Edd}\,c^{-2}$ from the H$\beta$ line. The decomposition of its $HST$ images yields a host stellar mass of $M_* = 10^{11.3 \pm 0.7} M_\odot$, and a ratio of $M_{\bullet}/M_*\approx 2.0\times 10^{-3}$ in agreement with the Magorrian relation. In the near future, it is expected to compare the geometrically-thick BLR discovered by the GRAVITY in 3C 273 with its spatially-resolved torus in order to understand the potential connection between the BLR and the torus.

Magnetars are neutron stars with extremely high surface magnetic fields. They show diverse X-ray pulse profiles in the quiescent state. We perform a systematic Fourier analysis of their soft X-ray pulse profiles. We find that most magnetars have a single-peaked profile and hence have low amplitudes of the second Fourier harmonics (A2). On the other hand, the pulsed fraction (PF) spreads over a wide range. We compared the results with theoretical profiles assuming various surface hotspot asymmetries, viewing geometries, and beaming functions. We found that a single value of the intensity ratio r between two antipodal hotspots is unable to reproduce the observed distribution of A2 and PF for all magnetars. The inferred r is probably anticorrelated with the thermal luminosity, implying that high-luminosity magnetars tend to have two symmetric hotspots. Our results are consistent with theoretical predictions, for which the existence of an evolving toroidal magnetic field breaks the symmetry of the surface temperature.

We investigate galactic winds driven by supernova (SN) explosions in an isolated dwarf galaxy using high-resolution (particle mass $m_{\rm gas} = 1{\rm M_\odot}$, number of neighbor $N_{\rm ngb} = 100$) smoothed-particle hydrodynamics simulations that include non-equilibrium cooling and chemistry, individual star formation, stellar feedback and metal enrichment. Clustered SNe lead to the formation of superbubbles which break out of the disk and vent out hot gas, launching the winds. We find much weaker winds than what cosmological simulations typically adopt at this mass scale. At the virial radius, the time-averaged loading factors of mass, momentum and energy are 3, 1 and 0.05, respectively, and the metal enrichment factor is 1.5. Winds that escape the halo consist of two populations that differ in their launching temperatures. Hot gas acquires enough kinetic energy to escape when launched while warm gas does not. However, warm gas can be further accelerated by the ram pressure of the subsequently launched hot gas and eventually escape. The strong interactions between different temperature phases highlight the caveat of extrapolating properties of warm gas to large distances based on its local conditions (e.g. the Bernoulli parameter). Our convergence study finds that wind properties converge when the cooling masses of individual SNe are resolved, which corresponds to $m_{\rm gas}=5 {\rm M_\odot}$ with an injection mass of $500 {\rm M_\odot}$. The winds weaken dramatically once the SNe become unresolved. We demonstrate that injecting the terminal momentum of SNe, a popular sub-grid model in the literature, fails to capture SN winds irrespective of the inclusion of residual thermal energy.

This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$\beta$ line measured from the mean spectrum. The virial factors for the other measures of the H$\beta$ line width are also presented.

We have started a long-term reverberation mapping project using the Wyoming Infrared Observatory 2.3 meter telescope titled "Monitoring AGNs with H\beta\ Asymmetry" (MAHA). The motivations of the project are to explore the geometry and kinematics of the gas responsible for complex H\beta\ emission-line profiles, ideally leading to an understanding of the structures and origins of the broad-line region (BLR). Furthermore, such a project provides the opportunity to search for evidence of close binary supermassive black holes. We describe MAHA and report initial results from our first campaign, from December 2016 to May 2017, highlighting velocity-resolved time lags for four AGNs with asymmetric H\beta\ lines. We find that 3C 120, Ark 120, and Mrk 6 display complex features different from the simple signatures expected for pure outflow, inflow, or a Keplerian disk. While three of the objects have been previously reverberation mapped, including velocity-resolved time lags in the cases of 3C 120 and Mrk 6, we report a time lag and corresponding black hole mass measurement for SBS 1518+593 for the first time. Furthermore, SBS 1518+593, the least asymmetric of the four, does show velocity-resolved time lags characteristic of a Keplerian disk or virialized motion more generally. Also, the velocity-resolved time lags of 3C 120 have significantly changed since previously observed, indicating an evolution of its BLR structure. Future analyses of the data for these objects and others in MAHA will explore the full diversity of H\beta\ lines and the physics of AGN BLRs.

NGC 5548 has been well spectroscopically monitored for reverberation mapping of the central kinematics by 19 campaigns. Using the maximum entropy method in this Letter, we build up a high-quality velocity-delay map of the H$\beta$ emission line in the light curves of the continuum and the line variations observed between 2015-2016. The map shows the response strength and lags of the velocity fields of the H$\beta$ emitting regions. The velocity-delay structure of the map is generally symmetric, with strong red and blue wings at time lag $\tau \leq 15$ days, a narrower velocity distribution at $\tau \geq 15$ days, and a deficit of response in the core. This is suggestive of a disk geometry of the broad-line region (BLR). The relatively weaker H$\beta$ response at the longer lags in the red side indicates anisotropic emission from the outer part of the BLR. We also recover the velocity-delay maps of NGC 5548 from the historical data of 13 years to investigate the long-term variability of its BLR. In general, the BLR of NGC 5548 was switching between the inflow and virialized phases in the past years. The resultant maps of seven years reveal inflow signatures and show decreasing lags, indicating that the changes in the BLR size are related to the infalling BLR gas. The other four maps show potential disk signatures which are similar to our map.

PSR J1119-6127 is a radio pulsar that behaved magnetar like bursts, and we performed a comprehensive investigation of this pulsar using the archival high-energy observations obtained after its outburst in 2016 July. After the 2016 outburst, specific regions on the neutron star surface were heated up to > 0.3 and > 1 keV from about 0.2 keV. A hard non-thermal spectral component with a photon index < 0.5 related to the magnetospheric emission can be resolved from the NuSTAR spectra above 10 keV. We find that the thermal emitting regions did not cool down and gradually shrank by about (20-35)% four months after the outburst. Hard X-ray pulsations were detected with NuSTAR immediately after the outburst, at 5 sigma confidence level and with a background-subtracted pulsed fraction of (40+-10)%. However, the signal became undetectable after a few days. Using Fermi data, we found that the gamma-ray emission in 0.5-300 GeV was suppressed along with the disappearance of the radio pulsations. This is likely caused by a reconfiguration of the magnetic field. We also discovered that the timing noise evolved dramatically, and the spin-down rate significantly increased after the 2016 glitch. We proposed that post-outburst temporal and spectral behaviors from radio to gamma-ray bands were caused by changes of the magnetosphere structure, pair plasma injection and the shrinking emission sites on the neutron star.

The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than the disc scale height that are resolved in modern galaxy formation simulations. We present a new statistical method to derive the evolutionary timeline of molecular clouds and star-forming regions. By quantifying the excess or deficit of the gas-to-stellar flux ratio around peaks of gas or star formation tracer emission, we directly measure the relative rarity of these peaks, which allows us to derive their lifetimes. We present a step-by-step, quantitative description of the method and demonstrate its practical application. The method's accuracy is tested in nearly 300 experiments using simulated galaxy maps, showing that it is capable of constraining the molecular cloud lifetime and feedback time-scale to $<0.1$ dex precision. Access to the evolutionary timeline provides a variety of additional physical quantities, such as the cloud-scale star formation efficiency, the feedback outflow velocity, the mass loading factor, and the feedback energy or momentum coupling efficiencies to the ambient medium. We show that the results are robust for a wide variety of gas and star formation tracers, spatial resolutions, galaxy inclinations, and galaxy sizes. Finally, we demonstrate that our method can be applied out to high redshift ($z\lesssim4$) with a feasible time investment on current large-scale observatories. This is a major shift from previous studies that constrained the physics of star formation and feedback in the immediate vicinity of the Sun.

Transient ultraluminous X-ray sources (ULXs) provide an important link bridging transient low-mass X-ray binaries and ULXs. Here we report the first discovery of both a canonical sub-Eddington outburst and an ultraluminous super-Eddington outburst from an unusual transient ULX, NGC 7793 P9 with a variability factor higher than $10^3$. Its X-ray spectrum switches between the typical high/soft state and the steep power-law state during the canonical outburst. The inner radius of the accretion disk and the disk temperature--luminosity correlation suggest that P9 harbors a stellar-mass black hole (BH). At the beginning of the ultraluminous outburst, we observe a cool outer disk with a hard Comptonized spectrum, implying a transition to the ULX regime. When the luminosity increases to $L\gtrsim3\times10^{39}$ erg/s, P9 shows a significantly curved spectrum that can be described by either a slim disk or a strongly curved Comptonized corona. The phenomenological model suggests that the hot disk observed near the peak of the ultraluminous outburst is coincidentally consistent with the extension of the thermal track. Utilizing more physical Comptonized disk models, we suggest that the corona cools down and the apparent disk-like spectrum is dominated by soft Comptonization. The significant variability above 1 keV supports this two-component scenario. The spectral evolution can also be interpreted with the supercritical accretion model. All these indicate that a canonical black hole X-ray binary can show properties of a ULX.

As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their H\beta\ and H\gamma\ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of H\beta\ and H\gamma\ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its H\beta\ and H\gamma\ lines show signatures of outflow. The H\beta\ and H\gamma\ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the H\beta\ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The H\beta\ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object.

PSR~J2021+4026 showed a sudden decrease in the gamma-ray emission at the glitch that occurred around 2011, October 16, and a relaxation of the flux to the pre-glitch state at around 2014 December. We report X-ray analysis results of the data observed by XMM-Newton on 2015 December 20 in the post-relaxation state. To examine any change in the X-ray emission, we compare the properties of the pulse profiles and spectra at the low gamma-ray flux state and at the post-relaxation state. The phase-averaged spectra for both states can be well described by a power-law component plus a blackbody component. The former is dominated by unpulsed emission and is probably originated from the pulsar wind nebula as reported by Hui et al (2015). The emission property of the blackbody component is consistent with the emission from the polar cap heated by the back-flow bombardment of the high-energy electrons or positrons that were accelerated in the magnetosphere. We found no significant change in the X-ray emission properties between two states. We suggest that the change of the X-ray luminosity is at an order of ~4%, which is difficult to measure with the current observations. We model the observed X-ray light curve with the heated polar cap emission and we speculate that the observed large pulsed fraction is owing to asymmetric magnetospheric structure.

As one of the series of papers reporting on a large reverberation mapping campaign of super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs), we present the results of 10 SEAMBHs monitored spectroscopically during 2015-2017. Six of them are observed for the first time, and have generally higher 5100 \AA\ luminosities than the SEAMBHs monitored in our campaign from 2012 to 2015; the remaining four are repeat observations to check if their previous lags change. Similar to the previous SEAMBHs, the H$\beta$ time lags of the newly observed objects are shorter than the values predicted by the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation of sub-Eddington AGNs, by factors of $\sim2-6$, depending on the accretion rate. The four previously observed objects have lags consistent with previous measurements. We provide linear regressions for the $R_{\mathrm{H\beta}}$-$L_{5100}$ relation, solely for the SEAMBH sample and for low-accretion AGNs. We find that the relative strength of Fe II and the profile of the H$\beta$ emission line can be used as proxies of accretion rate, showing that the shortening of H$\beta$ lags depends on accretion rates. The recent SDSS-RM discovery of shortened H$\beta$ lags in AGNs with low accretion rates provides compelling evidence for retrograde accretion onto the black hole. These evidences show that the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation holds only in AGNs with moderate accretion rates. At low accretion rates, it should be revised to include the effects of black hole spin, whereas the accretion rate itself becomes a key factor in the regime of high accretion rates.

We report the analysis result of UV/X-ray emission from AR~Scorpii, which is an intermediate polar (IP) composed of a magnetic white dwarf and a M-type star, with the XMM-Newton data. The X-ray/UV emission clearly shows a large variation over the orbit, and their intensity maximum (or minimum) is located at the superior conjunction (or inferior conjunction) of the M-type star orbit. The hardness ratio of the X-ray emission shows a small variation over the orbital phase, and shows no indication of the absorption by an accretion column. These properties are naturally explained by the emission from the M-type star surface rather than from the accretion column on the WD's star similar to the usual IPs. Beside, the observed X-ray emission also modulates with WD's spin with a pulse fraction of $\sim 14\%$. The peak position is aligned in the optical/UV/X-ray band. This supports the hypothesis that the electrons in AR~Scorpii are accelerated to a relativistic speed, and emit non-thermal photons via the synchrotron radiation. In the X-ray bands, the evidence of the power-law spectrum is found in the pulsed component, although the observed emission is dominated by the optically thin thermal plasma emissions with several different temperatures. It is considered that the magnetic dissipation/reconnection process on the M-type star surface heats up the plasma to a temperature of several keV, and also accelerates the electrons to the relativistic speed. The relativistic electrons are trapped in the WD's closed magnetic field lines by the magnetic mirror effect. In this model, the observed pulsed component is explained by the emissions from the first magnetic mirror point.

Type 1 active galactic nuclei display broad emission lines, regarded as arising from photoionized gas moving in the gravitational potential of a supermassive black hole. The origin of this broad-line region gas is unresolved so far, however. Another component is the dusty torus beyond the broad-line region, likely an assembly of discrete clumps that can hide the region from some viewing angles and make them observationally appear as Type 2 objects. Here we report that these clumps moving within the dust sublimation radius, like the molecular cloud G2 discovered in the Galactic center, will be tidally disrupted by the hole, resulting in some gas becoming bound at smaller radii while other gas is ejected and returns to the torus. The clumps fulfill necessary conditions to be photoionized. Specific dynamical components of tidally disrupted clumps include spiral-in gas as inflow, circularized gas, and ejecta as outflow. We calculate various profiles of emission lines from these clouds, and find they generally agree with H$\beta$ profiles of Palomar-Green quasars. We find that asymmetry, shape and shift of the profiles strongly depend on [O III], luminosity, which we interpret as a proxy of dusty torus angles. Tidally disrupted clumps from the torus may represent the source of the broad-line region gas.

We report a systematic analysis of the spin, orbital, and superorbital modulations of 4U 0114+650, a high-mass X-ray binary consists of one of the slowest spinning neutron stars. Utilizing the dynamic power spectrum, we found that the spin period varied dramatically and anti-correlated with the long-term X-ray flux variation that can be observed using RXTE ASM, Swift BAT, and MAXI. The spin-up rate over the entire dataset is consistent with previously reported values; however, the local spin-up rate is considerably higher. The corresponding local spin-up timescale is comparable to the local spin-up rate of OAO 1657-415, indicating that 4U 0114+650 could also have a transient disk. Moreover, the spin period evolution shows two ~1000-day spin-down/random-walk epochs that appeared together with depressions of the superorbital modulation amplitude. This implies that the superorbital modulation was closely related to the presence of the accretion disk, which is not favored in the spin-down/random-walk epochs because the accretion is dominated by the direct wind accretion.The orbital period is stable during the entire time span; however, the orbital profile changes with time significantly. We found that the depth of the dip near the inferior conjunction of the companion is highly variable, which disfavors the eclipsing scenario. Moreover, the dip was less obvious during the spin-down/random-walk epochs, indicating its correlation with the accretion disk. Further monitoring in both X-ray and optical bands could reveal the establishment of the accretion disk in this system.

A systematic study on the long-term periodicities of known Galactic cataclysmic variables (CVs) was conducted. Among 1580 known CVs, 344 sources were matched and extracted from the Palomar Transient Factory (PTF) data repository. The PTF light curves were combined with the Catalina Real-Time Transient Survey (CRTS) light curves and analyzed. Ten targets were found to exhibit long-term periodic variability, which is not frequently observed in the CV systems. These long-term variations are possibly caused by various mechanisms, such as the precession of the accretion disk, hierarchical triple star system, magnetic field change of the companion star, and other possible mechanisms. We discuss the possible mechanisms in this study. If the long-term period is less than several tens of days, the disk precession period scenario is favored. However, the hierarchical triple star system or the variations in magnetic field strengths are most likely the predominant mechanisms for longer periods.

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.

A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere.

We report on a timing and spectral analysis of the young, high magnetic field rotation-powered pulsar (RPP) B1509-58 using Chandra continuous-clocking mode observation. The pulsar's X-ray light curve can be fit by the two Gaussian components and the pulsed fraction shows moderate energy dependence over the Chandra band. The pulsed X-ray spectrum is well described by a power law with a photon index 1.16(4), which is harder than the values measured with RXTE PCA and NuSTAR. This result supports the log-parabolic model for the broadband X-ray spectrum. With the unprecedented angular resolution of Chandra, we clearly identified off-pulse X-ray emission from the pulsar, and its spectrum is best fit by a power law plus blackbody model. The latter component has a temperature of ~0.14 keV with a bolometric luminosity comparable to the luminosities of other young and high magnetic field RPPs, and it lies between the temperature of magnetars and typical RPPs. In addition, we found that the nonthermal X-ray emission of PSR B1509-58 is significantly softer in the off-pulse phase than in the pulsed phase, with the photon index varying between 1.0 and 1.8 and anticorrelated with the flux. This is similar to the behavior of three other young pulsars. We interpreted it as different contributions of pair-creation processes at different altitudes from the neutron star surface according to the outer-gap model.

We present high-resolution hydrodynamical simulations of isolated dwarf galaxies including self-gravity, non-equilibrium cooling and chemistry, interstellar radiation fields (IRSF) and shielding, star formation, and stellar feedback. This includes spatially and temporally varying photoelectric (PE) heating, photoionization, resolved supernova (SN) blast waves and metal enrichment. A new flexible method to sample the stellar initial mass function allows us to follow the contribution to the ISRF, the metal output and the SN delay times of individual massive stars. We find that SNe play the dominant role in regulating the global star formation rate, shaping the multi-phase interstellar medium (ISM) and driving galactic outflows. Outflow rates (with mass-loading factors of a few) and hot gas fractions of the ISM increase with the number of SNe exploding in low-density environments where radiative energy losses are low. While PE heating alone can suppress star formation slightly more (a factor of a few) than SNe alone can do, it is unable to drive outflows and reproduce the multi-phase ISM that emerges naturally when SNe are included. These results are in conflict with recent results of Forbes et al. who concluded that PE heating is the dominant process suppressing star formation in dwarfs, about an order of magnitude more efficient than SNe. Potential origins for this discrepancy are discussed. In the absence of SNe and photoionization (mechanisms to disperse dense clouds), the impact of PE heating is highly overestimated owing to the (unrealistic) proximity of dense gas to the radiation sources. This leads to a substantial boost of the infrared continuum emission from the UV-irradiated dust and a far infrared line-to-continuum ratio too low compared to observations. Though sub-dominant in regulating star formation, the ISRF controls the abundance of molecular hydrogen via photodissociation.

In this work, which is a continuation of Castell\'o-Mor et al. (2016), we present new X-ray and infrared (IR) data for a sample of active galactic nuclei (AGN) covering a wide range in Eddington ratio over a small luminosity range. In particular, we rigorously explore the dependence of the optical-to-X-ray spectral index $\alpha_{OX}$ and the IR-to-optical spectral index on the dimensionless accretion rate, $\dot{\mathcal{M}}=\dot{m}/\eta$ where $\dot{m}=L_{AGN}/L_{Edd}$ and $\eta$ is the mass-to-radiation conversion efficiency, in low and high accretion rate sources. We find that the SED of the faster accreting sources are surprisingly similar to those from the comparison sample of sources with lower accretion rate. In particular: i) the optical-to-UV AGN SED of slow and fast accreting AGN can be fitted with thin AD models. ii) The value of $\alpha_{OX}$ is very similar in slow and fast accreting systems up to a dimensionless accretion rate $\dot{\mathcal{M}}_{c}\sim$10. We only find a correlation between $\alpha_{OX}$ and $\dot{\mathcal{M}}$ for sources with $\dot{\mathcal{M}} > \dot{\mathcal{M}}_{c}$. In such cases, the faster accreting sources appear to have systematically larger $\alpha_{OX}$ values. iii) We also find that the torus in the faster accreting systems seems to be less efficient in reprocessing the primary AGN radiation having lower IR-to-optical spectral slopes. These findings, failing to recover the predicted differences between the SEDs of slim and thin ADs within the observed spectral window, suggest that additional physical processes or very special geometry act to reduce the extreme UV radiation in fast accreting AGN. This may be related to photon trapping, strong winds, and perhaps other yet unknown physical processes.

NGC 7793 P13 is an ultraluminous X-ray source harboring an accreting pulsar. We report on the detection of a ~65 d period X-ray modulation with Swift observations in this system. The modulation period found in the X-ray band is P=65.05+/-0.10 d and the profile is asymmetric with a fast rise and a slower decay. On the other hand, the u-band light curve collected by Swift UVOT confirmed an optical modulation with a period of P=64.24+/-0.13 d. We explored the phase evolution of the X-ray and optical periodicities and propose two solutions. A superorbital modulation with a period of ~2,700-4,700 d probably caused by the precession of a warped accretion disk is necessary to interpret the phase drift of the optical data. We further discuss the implication if this ~65d periodicity is caused by the superorbital modulation. Estimated from the relationship between the spin-orbital and orbital-superorbital periods of known disk-fed high-mass X-ray binaries, the orbital period of P13 is roughly estimated as 3-7 d. In this case, an unknown mechanism with a much longer time scale is needed to interpret the phase drift. Further studies on the stability of these two periodicities with a long-term monitoring could help us to probe their physical origins.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.

We report the Chandra/HRC-S and Swift/XRT observations for the 2015 outburst of the high-mass X-ray binary (HMXB) pulsar in the Small Magellanic Cloud, SMC X-2. While previous studies suggested that either an O star or a Be star in the field is the high-mass companion of SMC X-2, our Chandra/HRC-S image unambiguously confirms the O-type star as the true optical counterpart. Using the Swift/XRT observations, we extracted accurate orbital parameters of the pulsar binary through a time of arrivals (TOAs) analysis. In addition, there were two X-ray dips near the inferior conjunction, which are possibly caused by eclipses or an ionized high-density shadow wind near the companion's surface. Finally, we propose that an outflow driven by the radiation pressure from day ~10 played an important role in the X-ray/optical evolution of the outburst.

We propose a method for the flux calibration of reverberation mapping spectra based on accurate measurement of [O III] $\lambda 5007$ emission by spectral fitting. The method can achieve better accuracy than the traditional method of van Groningen & Wanders (1992), allowing reverberation mapping measurements for object with variability amplitudes as low as $\sim$ 5%. As a demonstration, we reanalyze the data of the Seyfert 1 galaxy MCG--6-30-15 taken from the 2008 campaign of the Lick AGN Monitoring Project, which previously failed to obtain a time lag for this weakly variable object owing to a relatively large flux calibration uncertainty. We detect a statistically significant rest-frame time lag of $6.38_{-2.69}^{+3.07}$ days between the H$\beta$ and $V$-band light curves. Combining this lag with FWHM(H$\beta$) = $1933\pm81$ $\rm km~s^{-1}$ and a virial coefficient of $f$ = 0.7, we derive a virial black hole mass of $3.26_{-1.40}^{+1.59}\times10^6$ $M_{\odot}$, which agrees well with previous estimates by other methods.

We report our measurements for orbital and spin parameters of X 1822-371 using its X-ray partial eclipsing profile and pulsar timing from data collected by the Rossi X-ray Timing Explorer (RXTE). Four more X-ray eclipse times obtained by the RXTE 2011 observations were combined with historical records to trace evolution of orbital period. We found that a cubic ephemeris likely better describes evolution of the X-ray eclipse times during a time span of about 34 years with a marginal second order derivative of $\ddot{P}_{orb}=(-1.05 \pm 0.59) \times 10^{-19}$ s$^{-1}$. Using the pulse arrival time delay technique, the orbital and spin parameters were obtained from RXTE observations from 1998 to 2011. The detected pulse periods show that the neutron star in X 1822-371 is continuously spun-up with a rate of $\dot{P}_{s}=(-2.6288 \pm 0.0095) \times 10^{-12}$ s s$^{-1}$. Evolution of the epoch of the mean longitude $l=\pi /2$ (i.e. $T_{\pi / 2}$) gives an orbital period derivative value consistent with that obtained from the quadratic ephemeris evaluated by the X-ray eclipse but the detected $T_{\pi / 2}$ values are significantly and systematically earlier than the corresponding expected X-ray eclipse times by $90 \pm 11$ s. This deviation is probably caused by asymmetric X-ray emissions. We also attempted to constrain the mass and radius of the neutron star using the spin period change rate and concluded that the intrinsic luminosity of X 1822-371 is likely more than $10^{38}$ ergs s$^{-1}$.

The black hole mass measurement in active galaxies is a challenge, particularly in sources where the reverberation method cannot be applied. We aim to determine the black hole mass in a very special object, RE J1034+396, one of the two AGN with QPO oscillations detected in X-rays, and a single bright AGN with optical band totally dominated by starlight. We fit the stellar content using the code starlight, and the broad band disk contribution to optical/UV/X-ray emission is modeled with optxagnf. We also determine the black hole mass using several other independent methods. Various methods give contradictory results. Most measurements of the blacc hole mass are in the range 1.e6-1.e7 Msun, and the measurements based on dynamics give higher values than measurements based on Hbeta and Mg II emission lines.

NGC 5548 is the best-observed reverberation-mapped active galactic nucleus with long-term, intensive monitoring. Here we report results from a new observational campaign between January and July, 2015. We measure the centroid time lag of the broad H$\beta$ emission line with respect to the 5100 \AA continuum and obtain $\tau_{\rm cent} = 7.20^{+1.33}_{-0.35}$ days in the rest frame. This yields a black hole mass of $M_{\bullet}=8.71^{+3.21}_{-2.61} $x$ 10^{7}M_{\odot}$ using a broad H$\beta$ line dispersion of $3124\pm302$ km s$^{-1}$ and a virial factor of $f_{_{\rm BLR}}=6.3\pm1.5$ for the broad-line region (BLR), consistent with the mass measurements from previous H$\beta$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$\beta$ line, which shows a symmetric response pattern around the line center, a plausible kinematic signature of virialized motion of the BLR. Combining all the available measurements of H$\beta$ time lags and the associated mean 5100 {\AA} luminosities over 18 campaigns between 1989 and 2015, we find that the H$\beta$ BLR size varies with the mean optical luminosity, but, interestingly, with a possible delay of $2.35_{-1.25}^{+3.47}$ yrs. This delay coincides with the typical BLR dynamical timescale of NGC 5548, indicating that the BLR undergoes dynamical changes, possibly driven by radiation pressure.

We report the optical spectroscopy of four young radio sources which are observed with the Lijiang 2.4m telescope. The Eddington ratios of these sources are similar with those of narrow-line Seyfert 1 galaxies (NLS1s). Their Fe {\sc ii} emission is strong while [O {\sc iii}] strength is weak. These results confirm the NLS1 features of young radio sources, except that the width of broad H$\beta$ of young radio sources is larger than that of NLS1s. We thus suggest that the young radio sources are the high black hole mass counterparts of steep-spectrum radio-loud NLS1s. In addition, the broad H$\beta$ component of \astrobj{4C 12.50} is the blue wing of the narrow component, but not from the broad line region.

We report a possible detection of a 55-day X-ray modulation for the ultraluminous accreting pulsar M82 X-2 from archival Chandra observations. Because M82 X-2 is known to have a 2.5-day orbital period, if the 55-day period is real, it will be the superorbital period of the system. We also investigated variabilities of other three nearby ultraluminous X-ray sources in the central region of M82 with the Chandra data and did not find any evidence of periodicities. Furthermore, we re-examined the previously reported 62-day periodicity near the central region of M82 by performing a systematic timing study with all the archival Rossi X-Ray Timing Explorer and Swift data. Using various dynamic timing analysis methods, we confirmed that the 62-day period is not stable, suggesting that it is not the orbital period of M82 X-1 in agreement with previous work.

The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft $\gamma-$ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the \emph{Fermi} Large Area Telescope (LAT), we identified an extended LAT source to the south-west of Cl* 1806-20. The centroid of the 1-50~GeV emission is consistent with that of HESS J1808-204 (until now unidentified). The LAT spectrum is best-fit by a broken power-law with the break energy $E_\mathrm{b}=297\pm15$ MeV. The index above $E_\mathrm{b}$ is $2.60\pm0.04$, and is consistent with the flux and spectral index above 100 GeV for HESS J1808-204, suggesting an association between the two sources. Meanwhile, the interacting supernova remnant SNR G9.7-0.0 is also a potential contributor to the LAT flux. A tentative flux enhancement at the MeV band during a 45-day interval (2011 Jan 21 - 2011 Mar 7) is also reported. We discuss possible origins of the extended LAT source in the context of both leptonic and hadronic scenarios.

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

Broad emission lines in active galactic nuclei (AGNs) mainly arise from gas photoionized by continuum radiation from an accretion disk around a central black hole. The shape of the broad-line profile, described by ${\cal D}_{_{\rm H\beta}}={\rm FWHM}/\sigma_{_{\rm H\beta}}$, the ratio of full width at half maximum to the dispersion of broad H$\beta$, reflects the dynamics of the broad-line region (BLR) and correlates with the dimensionless accretion rate ($\dot{\mathscr{M}}$) or Eddington ratio ($L_{\rm bol}/L_{\rm Edd}$). At the same time, $\dot{\mathscr{M}}$ and $L_{\rm bol}/L_{\rm Edd}$ correlate with ${\cal R}_{\rm Fe}$, the ratio of optical Fe II to H$\beta$ line flux emission. Assembling all AGNs with reverberation mapping measurements of broad H$\beta$, both from the literature and from new observations reported here, we find a strong bivariate correlation of the form $\log(\dot{\mathscr{M}},L_{\rm bol}/L_{\rm Edd})=\alpha+\beta{\cal D}_{_{\rm H\beta}}+\gamma{\cal R}_{\rm Fe},$ where $\alpha=(2.47,0.31)$, $\beta=-(1.59,0.82)$ and $\gamma=(1.34,0.80)$. We refer to this as the fundamental plane of the BLR. We apply the plane to a sample of $z < 0.8$ quasars to demonstrate the prevalence of super-Eddington accreting AGNs are quite common at low redshifts.

We present the time-frequency analysis results based on the Hilbert-Huang transform (HHT) for the evolution of a 4-Hz low-frequency quasi-periodic oscillation (LFQPO) around the black hole X-ray binary XTE J1550-564. The origin of LFQPOs is still debated. To understand the cause of the peak broadening, we utilized a recently developed time-frequency analysis, HHT, for tracking the evolution of the 4-Hz LFQPO from XTE J1550 564. By adaptively decomposing the ~4-Hz oscillatory component from the light curve and acquiring its instantaneous frequency, the Hilbert spectrum illustrates that the LFQPO is composed of a series of intermittent oscillations appearing occasionally between 3 Hz and 5 Hz. We further characterized this intermittency by computing the confidence limits of the instantaneous amplitudes of the intermittent oscillations, and constructed both the distributions of the QPO's high and low amplitude durations, which are the time intervals with and without significant ~4-Hz oscillations, respectively. The mean high amplitude duration is 1.45 s and 90% of the oscillation segments have lifetimes below 3.1 s. The mean low amplitude duration is 0.42 s and 90% of these segments are shorter than 0.73 s. In addition, these intermittent oscillations exhibit a correlation between the oscillation's rms amplitude and mean count rate. This correlation could be analogous to the linear rms-flux relation found in the 4-Hz LFQPO through Fourier analysis. We conclude that the LFQPO peak in the power spectrum is broadened owing to intermittent oscillations with varying frequencies, which could be explained by using the Lense-Thirring precession model.

Long-term X-ray modulations on timescales from tens to hundreds of days have been widely studied for X-ray binaries located in the Milky Way and the Magellanic Clouds. For other nearby galaxies, only the most luminous X-ray sources can be monitored with dedicated observations. We here present the first systematic study of long-term X-ray variability of four ultraluminous X-ray sources (ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, and NGC 5408 X-1) monitored with Swift. By using various dynamic techniques to analyse their light curves, we find several interesting low-frequency quasi-periodicities. Although the periodic signals may not represent any stable orbital modulations, these detections reveal that such long-term regular patterns may be related to superorbital periods and structure of the accretion discs. In particular, we show that the outburst recurrence time of ESO 243-49 HLX-1 varies over time and suggest that it may not be the orbital period. Instead, it may be due to some kinds of precession, and the true binary period is expected to be much shorter.

Here we report the results of searching millisecond pulsar (MSP) candidates from the Fermi LAT second source catalog (2FGL). Seven unassociated $\gamma-$ray sources in this catalog are identified as promising MSP candidates based on their $\gamma$-ray properties. Through the X-ray analysis, we have detected possible X-ray counterparts, localized to an arcsecond accuracy. We have systematically estimated their X-ray fluxes and compared with the corresponding $\gamma$-ray fluxes. The X-ray to $\gamma$-ray flux ratios for 2FGL J1653.6-0159 and 2FGL J1946.4-5402 are comparable with the typical value for pulsars. For 2FGL J1625.2-0020, 2FGL J1653.6-0159 and 2FGL J1946.4-5402, their candidate X-ray counterparts are bright enough for performing a detailed spectral and temporal analysis to discriminate their thermal/non thermal nature and search for the periodic signal. We have also searched for possible optical/IR counterparts at the X-ray positions. For the optical/IR source coincident with the brightest X-ray object that associated with 2FGL J1120.0-2204, its spectral energy distribution is comparable with a late-type star. Evidence for the variability has also been found by examining its optical light curve. All the aforementioned 2FGL sources resemble a pulsar in one or more aspects, which make them as the promising targets for follow-up investigations.

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$\beta$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$\beta$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm H\beta}}$) and optical luminosity at 5100\AA, $R_{_{\rm H\beta}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm H\beta}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

We have investigated the intrabinary shock emission from the redback millisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modulation in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak structure with a dip in between. The observed X-rays are non-thermal dominant which can be modeled by a power-law with a photon index of ~1.2. Intrabinary shock can be the origin of the observed X-rays. The UV light curve is resulted from the ellipsoidal modulation of the companion. Modeling the UV light curve prefers a large viewing angle. The heating effect of the UV light curve is found to be negligible which suggests the high energy radiation beam of PSR J2129-0429 does not direct toward its companion. On the other hand, no significant orbital modulation can be found in gamma-ray which suggests the majority of the gamma-rays come from the pulsar.

We have investigated the field around the radio-quiet $\gamma$-ray pulsar, PSR J2021+4026, with a ~140 ks XMM-Newton observation and a ~56 ks archival Chandra data. Through analyzing the pulsed spectrum, we show that the X-ray pulsation is purely thermal in nature which suggests the pulsation is originated from a hot polar cap with $T\sim3\times10^{6}$ K on the surface of a rotating neutron star. On the other hand, the power-law component that dominates the pulsar emission in the hard band is originated from off-pulse phases, which possibly comes from a pulsar wind nebula. In re-analyzing the Chandra data, we have confirmed the presence of bow-shock nebula which extends from the pulsar to west by ~10 arcsec. The orientation of this nebular feature suggests that the pulsar is probably moving eastward which is consistent with the speculated proper motion by extrapolating from the nominal geometrical center of the supernova remnant (SNR) G78.2+2.1 to the current pulsar position. For G78.2+2.1, our deep XMM-Newton observation also enables a study of the central region and part of the southeastern region with superior photon statistics. The column absorption derived for the SNR is comparable with that for PSR J2021+4026, which supports their association. The remnant emission in both examined regions are in an non-equilibrium ionization state. Also, the elapsed time of both regions after shock-heating is apparently shorter than the Sedov age of G78.2+2.1. This might suggest the reverse shock has reached the center not long ago. Apart from PSR J2021+4026 and G78.2+2.1, we have also serendipitously detected an X-ray flash-like event XMM J202154.7+402855 from this XMM-Newton observation.

We report multi-wavelength observations of the unidentified Fermi object 2FGL J1653.6-0159. With the help of high-resolution X-ray observation, we have identified an X-ray and optical counterpart of 2FGL J1653.6-0159. The source exhibits a periodic modulation of 75 min in optical and possibly also in X-ray. We suggest that 2FGL J1653.6-0159 is a compact binary system with an orbital period of 75 min. Combining the gamma-ray and X-ray properties, 2FGL J1653.6-0159 is potentially a black widow/redback type gamma-ray millisecond pulsar (MSP). The optical and X-ray lightcurve profiles show that the companion is mildly heated by the high-energy emission and the X-rays are from intrabinary shock. Although no radio pulsation has been detected yet, we estimated that the spin period of the MSP is ~2ms based on a theoretical model. If pulsation can be confirmed in the future, 2FGL J1653.6-0159 will become the first ultracompact rotation-powered MSP.

We report the X-ray pulsation of ~173.3 ms for the "next Geminga", PSR J1836+5925, with recent XMM-Newton investigations. The X-ray periodicity is consistent wtih the gamma-ray ephemeris at the same epoch. The X-ray folded light curve has a sinusoidal structure which is different from the double-peaked gamma-ray pulse profile. We have also analysed the X-ray phase-averaged spectra which shows the X-ray emission from PSR J1836+5925 is thermal dominant. This suggests the X-ray pulsation mainly originates from the modulated hot spot on the stellar surface.

This is the second in a series of papers reporting on a large reverberation mapping (RM) campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). The goal is to identify super-Eddington accreting massive black holes (SEAMBHs) and to use their unique properties to construct a new method for measuring cosmological distances. Based on theoretical models, the saturated bolometric luminosity of such sources is proportional to the BH mass which can be used to obtain their distance. Here we report on five new RM measurements and show that in four of the cases we can measure the BH mass and three of these sources are SEAMBHs. Together with the three sources from our earlier work, we now have six new sources of this type. We use a novel method based on a minimal radiation efficiency to identify nine additional SEAMBHs from earlier RM-based mass measurements. We use a Bayesian analysis to determine the parameters of the new distance expression, and the method uncertainties, from the observed properties of the objects in the sample. The ratio of the newly measured distances to the standard cosmological ones has a mean scatter of 0.14 dex, indicating that SEAMBHs can be use as cosmological distance probes. With their high luminosity, long period of activity and large numbers at high redshifts, SEAMBHs have a potential to extend the cosmic distance ladder beyond the range now explored by type Ia supernovae.

A tight relation between the size of the broad-line region (BLR) and optical luminosity has been established in about 50 active galactic nuclei studied through reverberation mapping of the broad Hbeta emission line. The R_blr-L relation arises from simple photoionization considerations. Using a general relativistic model of an optically thick, geometrically thin accretion disk, we show that the ionizing luminosity jointly depends on black hole mass, accretion rate, and spin. The non-monotonic relation between the ionizing and optical luminosity gives rise to a complicated relation between the BLR size and the optical luminosity. We show that the reverberation lag of Hbeta to the varying continuum depends very sensitively to black hole spin. For retrograde spins, the disk is so cold that there is a deficit of ionizing photons in the BLR, resulting in shrinkage of the hydrogen ionization front with increasing optical luminosity, and hence shortened Hbeta lags. This effect is specially striking for luminous quasars undergoing retrograde accretion, manifesting in strong deviations from the canonical R_blr-L relation. This could lead to a method to estimate black hole spins of quasars and to study their cosmic evolution. At the same time, the small scatter of the observed R_blr-L relation for the current sample of reverberation-mapped active galaxies implies that the majority of these sources have rapidly spinning black holes.

We present the smoothed-particle hydrodynamics implementation SPHGal, which combines some recently proposed improvements in GADGET. This includes a pressure-entropy formulation with a Wendland kernel, a higher order estimate of velocity gradients, a modified artificial viscosity switch with a modified strong limiter, and artificial conduction of thermal energy. With a series of idealized hydrodynamic tests we show that the pressure-entropy formulation is ideal for resolving fluid mixing at contact discontinuities but performs conspicuously worse at strong shocks due to the large entropy discontinuities. Including artificial conduction at shocks greatly improves the results. In simulations of Milky Way like disk galaxies a feedback-induced instability develops if too much artificial viscosity is introduced. Our modified artificial viscosity scheme prevents this instability and shows efficient shock capturing capability. We also investigate the star formation rate and the galactic outflow. The star formation rates vary slightly for different SPH schemes while the mass loading is sensitive to the SPH scheme and significantly reduced in our favored implementation. We compare the accretion behavior of the hot halo gas. The formation of cold blobs, an artifact of simple SPH implementations, can be eliminated efficiently with proper fluid mixing, either by conduction and/or by using a pressure-entropy formulation.

RE J1034+396, a narrow-line Seyfert-1 active galactic nucleus (AGN), is the first example of AGNs that exhibited a nearly coherent quasi-periodic oscillation (QPO) for the data collected by XMM-Newton in 2007. The spectral behaviors and timing properties of the QPO have been studied since its discovery. We present an analysis of the QPO in RE J1034+396 based on the Hilbert-Huang transform (HHT). Comparing other time-frequency analysis methods, the Hilbert spectrum reveals the variation of the QPO period in great detail. Furthermore, the empirical mode decomposition provides band-pass filtered data that can be used in the O -- C and correlation analysis. We suggest that it is better to divide the evolution of the QPO in this observation into three epochs according to their different periodicities. In addition to the periodicities, the correlations between the QPO periods and corresponding mean count rates are also different in these three epochs. Further examining the phase lags in these epochs, we found no significant phase lags between the soft and hard X-ray bands, which is also confirmed in the QPO phase-resolved spectral analysis. Finally, we discuss the indications of current models including a spotted accretion disk, diskoseismology, and oscillation of shock, according to the observed time-frequency and spectral behaviors.

Flux calibration of spectra in reverberation mapping (RM) is most often performed by assuming the flux constancy of some specified narrow emission lines, which stem from an extended region that is sometimes partially spatially resolved, in contrast to the point-like broad-line region and the central continuum source. The inhomogeneous aperture geometries used among different observation sets in a joint monitoring campaign introduce systematic deviations to the fluxes of broad lines and central continuum, and intercalibration over these data sets is required. As an improvement to the previous empirical correction performed by comparing the (nearly) contemporaneous observation points, we describe a feasible Bayesian method that obviates the need for (nearly) contemporaneous observations, naturally incorporates physical models of flux variations, and fully takes into account the measurement errors. In particular, it fits all the data sets simultaneously regardless of samplings and makes use of all of the information in the data sets. A Markov Chain Monte Carlo implementation is employed to recover the parameters and uncertainties for intercalibration. Application to the RM data sets of NGC 5548 with joint monitoring shows the high fidelity of our method.

The metallicity of active galactic nuclei (AGNs), which can be measured by emission line ratios in their broad and narrow line regions (BLRs and NLRs), provides invaluable information about the physical connection between the different components of AGNs. From the archival databases of the International Ultraviolet Explorer, the Hubble Space Telescope and the Sloan Digital Sky Survey, we have assembled the largest sample available of AGNs which have adequate spectra in both the optical and ultraviolet bands to measure the narrow line ratio [N II]/H{\alpha} and also, in the same objects, the broad-line N V/C IV ratio. These permit the measurement of the metallicities in the NLRs and BLRs in the same objects. We find that neither the BLR nor the NLR metallicity correlate with black hole masses or Eddington ratios, but there is a strong correlation between NLR and BLR metallicities. This metallicity correlation implies that outflows from BLRs carry metal-rich gas to NLRs at characteristic radial distances of ~ 1.0 kiloparsec. This chemical connection provides evidence for a kinetic feedback of the outflows to their hosts. Metals transported into the NLR enhance the cooling of the ISM in this region, leading to local star formation after the AGNs turn to narrow line LINERs. This post-AGN star formation is predicted to be observable as an excess continuum emission from the host galaxies in the near infrared and ultraviolet, which needs to be further explored.

We report first results from a large project to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). Such objects may be different from other AGNs in being powered by slim accretion disks and showing saturated accretion luminosities, but both are not yet fully understood. The results are part of a large reverberation mapping (RM) campaign using the 2.4-m Shangri-La telescope at the Yunnan Observatory in China. The goals are to investigate the gas distribution near the BH and the properties of the central accretion disks, to measure BH mass and Eddington ratios, and to test the feasibility of using such objects as a new type of cosmological candles. The paper presents results for three objects, Mrk 335, Mrk 142 and IRAS F12397+3333 with H$\beta$ time lags relative to the 5100\AA\ continuum of $10.6^{+1.7}_{-2.9}$, $6.4^{+0.8}_{-2.2}$ and $11.4^{+2.9}_{-1.9}$ days, respectively. The corresponding BH masses are $(8.3_{-3.2}^{+2.6})\times 10^6M_{\odot}$, $(3.4_{-1.2}^{+0.5})\times 10^6M_{\odot}$ and $(7.5_{-4.1}^{+4.3})\times 10^6M_{\odot}$, and the lower limits on the Eddington ratios 0.6, 2.3, and 4.6 for the minimal radiative efficiency of 0.038. Mrk 142 and IRAS F12397+333 (extinction corrected) clearly deviate from the currently known relation between H$\beta$ lag and continuum luminosity. The three Eddington ratios are beyond the values expected in thin accretion disks and two of them are the largest measured so far among objects with RM-based BH masses. We briefly discuss implications for slim disks, BH growth and cosmology.

The high-mass X-ray binary SMC X-1 is an eclipsing binary with an orbital period of 3.89 d. This system exhibits a superorbital modulation with a period varying between ~40 d and ~65 d. The instantaneous frequency and the corresponding phase of the superorbital modulation can be obtained by a recently developed time-frequency analysis technique, the Hilbert-Huang transform (HHT). We present a phase-resolved analysis of both the spectra and the orbital profiles with the superorbital phase derived from the HHT. The X-ray spectra observed by the Proportional Counter Array onboard the Rossi X-ray Timing Explorer are fitted well by a blackbody plus a Comptonized component. The plasma optical depth, which is a good indicator of the distribution of material along the line of sight, is significantly anti-correlated with the flux detected at 2.5-25 keV. However, the relationship between the plasma optical depth and the equivalent width of the iron line is not monotonic: there is no significant correlation for fluxes higher than ~35 mCrab but clear positive correlation when the intensity is lower than ~20 mCrab. This indicates that the iron line production is dominated by different regions of this binary system in different superorbital phases. To study the dependence of the orbital profile on the superorbital phase, we obtained the eclipse profiles by folding the All Sky Monitor light curve with the orbital period for different superorbital states. A dip feature, similar to the pre-eclipse dip in Her X-1, lying at orbital phase ~0.6-0.85, was discovered during the superorbital transition state. This indicates that the accretion disk has a bulge that absorbs considerable X-ray emission in the stream-disk interaction region. The dip width is anti-correlated with the flux, and this relation can be interpreted by the precessing tilted accretion disk scenario.

We report the discovery of X-ray periodicity of ~265.3 ms from a deep XMM-Newton observation of the radio-quiet gamma-ray pulsar, PSR J2021+4026, located at the edge of the supernova remnant G78.2+2.1 (gamma-Cygni).The detected frequency is consistent with the gamma-ray pulsation determined by the observation of Fermi Gamma-ray Space Telescope at the same epoch. The X-ray pulse profile resembles the modulation of hot spot on the surface of the neutron star. The phase-averaged spectral analysis also suggests that the majority of the observed X-rays have a thermal origin. This is the third member in the class of radio-quiet pulsars with the significant pulsations detected from both X-rays and gamma-ray regimes.

Super-Eddington accreting massive black holes (SEAMBHs) reach saturated luminosities above a certain accretion rate due to photon trapping and advection in slim accretion disks. We show that these SEAMBHs could provide a new tool for estimating cosmological distances if they are properly identified by hard X-ray observations, in particular by the slope of their 2-10 keV continuum. To verify this idea we obtained black hole mass estimates and X-ray data for a sample of 60 narrow line Seyfert 1 galaxies that we consider to be the most promising SEAMBH candidates. We demonstrate that the distances derived by the new method for the objects in the sample get closer to the standard luminosity distances as the hard X-ray continuum gets steeper. The results allow us to analyze the requirements for using the method in future samples of active black holes and to demonstrate that the expected uncertainty, given large enough samples, can make them into a useful, new cosmological ruler.

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 \AA. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.

For some samples, it has been shown that spectra of QSOs with low redshift are bluer during their brighter phases. For the FIRST bright QSO sample, we assemble their spectra from SDSS DR7 to investigate variability between the spectra from White et al. (2000) and from the SDSS for a long rest-frame time-lag, up to 10 years. There are 312 radio loud QSOs and 232 radio quiet QSOs in this sample, up to $z \sim 3.5$. With two-epoch variation, it is found that spectra of half of the QSOs appear redder during their brighter phases. There is no obvious difference in slope variability between sub-samples of radio quiet and radio loud QSOs. This result implies that the presence of a radio jet does not affect the slope variability on 10-year timescales. The arithmetic composite difference spectrum for variable QSOs is steep at blueward of $\sim$ 2500\AA. The variability for the region blueward of 2500 \AA\ is different to that for the region redward of 2500 \AA.

We report the results from a detailed analysis of an archival XMM-Newton observation of the X-ray source XGPS-I J183251-100106, which has been suggested as a promising magnetic cataclysmic variable candidate based on its optical properties. A single periodic signal of 1.5 hrs is detected from all EPIC cameras on board XMM-Newton. The phase-averaged X-ray spectrum can be well-modeled with a thermal bremsstrahlung of a temperature kT~50 keV. Both X-ray spectral and temporal behavior of this system suggest it as a eclipsing cataclysmic variable of AM Herculis (or polar) type.

Recently, much attention has been given to double-peaked narrow emission-line galaxies, some of which are suggested to be related with merging galaxies. We make a systematic search to build the largest sample of these sources from Data Release 7 of Sloan Digital Sky Survey (SDSS DR7). With reasonable criteria of fluxes, full-width-half-maximum of emission lines and separations of the peaks, we select 3,030 double-peaked narrow emission-lines galaxies. In light of the existence of broad Balmer lines and the locations of the two components of double-peaked narrow emission lines distinguished by the Kauffmann et al. (2003) criteria in the Baldwin-Phillips-Terlevich (BPT) diagram, we find that there are 81 type I AGN, 837 double-type II AGN (2-type II), 708 galaxies with double star forming components (2-SF), 400 with mixed star forming and type II AGN components (type II + SF) and 1,004 unknown-type objects. As a by-product, a sample of galaxies (12,582) with asymmetric or top-flat profiles of emission lines is established. After inspecting the SDSS images of the two samples visually, we find 54 galaxies with dual cores. The present samples can be used to study the dynamics of merging galaxies, the triggering mechanism of black hole activity, the hierarchical growth of galaxies and dynamics of narrow line regions driven by outflows and rotating disk.

Recent x-ray observations of Mrk 766 suggest that broad emission line region clouds cross our line of sight and produce variable x-ray absorption. Here we investigate what optical/ultraviolet spectroscopic features would be produced by such "Intervening BLR Clouds" (IBC) crossing our line of sight to the accretion disk, the source of the optical/UV continuum. Although the emission spectrum produced by intervening clouds is identical to the standard BLR model, they may produce absorption features on the optical or UV continuum. Single clouds will have little effect on the optical/UV spectrum because BLR clouds are likely to be much smaller than the accretion disk. This is unlike the X-ray case, where the radiation source is considerably smaller. However, an ensemble of intervening clouds will produce spectroscopic features in the FUV including a strong depression between the Lyman limit and Ly{\alpha}. The amount of the depression will indicate the line-of-sight covering factor of clouds, an unknown quantity that is important for the ionization of the intergalactic medium and the energy budget of AGN. Comparison with observations suggests that the SED of Mrk 766 may be affected by intervening BLR clouds and IBC may exist in most of AGNs.

(abridged) We study the consequence of star formation (SF) in an self-gravity dominated accretion disk in quasars. The warm skins of the SF disk are governed by the radiation from the inner part of the accretion disk to form Compton atmosphere (CAS). The CAS are undergoing four phases to form broad line regions. Phase I is the duration of pure accumulation supplied by the SF disk. During phase II clouds begin to form due to line cooling and sink to the SF disk. Phase III is a period of preventing clouds from sinking to the SF disk through dynamic interaction between clouds and the CAS. Finally, phase IV is an inevitable collapse of the entire CAS through line cooling. This CAS evolution drives the episodic appearance of BLRs. Geometry and dynamics of BLRs can be self-consistently derived from the thermal instability of the CAS during phases II and III by linear analysis. The metallicity gradient of SF disk gives rise to different properties of clouds from outer to inner part of BLRs. We find that clouds have column density N_H < 10^22cm^{-2} in the metal-rich regions whereas they have N_H > 10^22 cm^{-2} in the metal-poor regions. The metal-rich clouds compose the high ionization line (HIL) regions whereas the metal-poor clouds are in low ionization line (LIL) regions. Metal-rich clouds in HIL regions will be blown away by radiation pressure, forming the observed outflows. The LIL regions are episodic due to the mass cycle of clouds with the CAS in response to continuous injection by the SF disk, giving rise to different types of AGNs. Based on SDSS quasar spectra, we identify a spectral sequence in light of emission line equivalent width from Phase I to IV. A key phase in the episodic appearance of the BLRs is bright type II AGNs with no or only weak BLRs. We discuss observational implications and tests of the theoretical predictions of this model.

The high-mass X-ray binary (HMXB) SMC X-1 exhibits a superorbital modulation with a dramatically varying period ranging between ~40 d and ~60 d. This research studies the time-frequency properties of the superorbital modulation of SMC X-1 based on the observations made by the All-Sky Monitor (ASM) onboard the Rossi X-ray Timing Explorer (RXTE).We analyzed the entire ASM database collected since 1996. The Hilbert-Huang Transform (HHT), developed for non-stationary and nonlinear time series analysis, was adopted to derive the instantaneous superorbital frequency. The resultant Hilbert spectrum is consistent with the dynamic power spectrum while it shows more detailed information in both the time and frequency domains. The RXTE observations manifest that the superorbital modulation period was mostly betweenn ~50 d and ~65 d, whenas it changed to ~45 d around MJD 50,800 and MJD 54,000. Our analysis further indicates that the instantaneous frequency changed in a time scale of hundreds of days between ~MJD 51,500 and ~MJD 53,500. Based on the instantaneous phase defined by HHT, we folded the ASM light curve to derive a superorbital profile, from which an asymmetric feature and a low state with barely any X-ray emissions (lasting for ~0.3 cycles) were observed. We also calculated the correlation between the mean period and the amplitude of the superorbital modulation. The result is similar to the recently discovered relationship between the superorbital cycle length and the mean X-ray flux for Her X-1.

It has been suggested that the high metallicity generally observed in active galactic nuclei (AGNs) and quasars originates from ongoing star formation in the self-gravitating part of accretion disks around the supermassive black holes. We designate this region as the star forming (SF) disk, in which metals are produced from supernova explosions (SNexp) while at the same time inflows are driven by SNexp-excited turbulent viscosity to accrete onto the SMBHs. In this paper, an equation of metallicity governed by SNexp and radial advection is established to describe the metal distribution and evolution in the SF disk. We find that the metal abundance is enriched at different rates at different positions in the disk, and that a metallicity gradient is set up that evolves for steady-state AGNs. Metallicity as an integrated physical parameter can be used as a probe of the SF disk age during one episode of SMBH activity. In the SF disk, evaporation of molecular clouds heated by SNexp blast waves unavoidably forms hot gas. This heating is eventually balanced by the cooling of the hot gas, but we show that the hot gas will escape from the SF disk before being cooled, and diffuse into the BLRs forming with a typical rate of $\sim 1\sunmyr$. The diffusion of hot gas from a SF disk depends on ongoing star formation, leading to the metallicity gradients in BLR observed in AGNs. We discuss this and other observable consequences of this scenario.

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov pulse in radio detection originates from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulse by making far- field approximation fails when the size of elongated showers become comparable with detection distance. We investigate the Cherenkov pulse in near-field by a numerical code based on the finite- difference time-domain (FDTD) method. Our study shows that the near-field radiation exhibits very different behaviors from the far-field one and therefore can be easily recognized. For ground array neutrino detectors, the near-field radiation would provide a unique signature for ultra high energy electromagnetic showers induced by the electron neutrino charge-current interaction. This can be useful in neutrino flavor identification.

We analyze the optical spectrum of type 1 QSO SDSS J1425+3231. This ob- ject is interesting since its narrow emission lines such as [O III]{\lambda}{\lambda}4959, 5007 are double- peaked, and the line structure can be modeled well by three Gaussian components: two components for the two peaks (we refer the peaks at low/high redshift as "the blue/red component") and another one for the line wing which has the same line center as that of the blue component, but ~ 3 times broader. The separation between the blue and red components is ~ 500 km/s with blue component ~ 2 times broader than the red one. The H{\beta} emission can be separated into four components: two for the double-peaked narrow line and two for the broad line which comes from the broad line region (BLRs). The black hole mass estimated from the broad H{\beta} emission line using the typical reverberation map- ping relation is 0.85 \times 108M\odot, which is consistent with that derived from parameters of [O III]{\lambda} 5007 of the blue component. We suggest this QSO might be a dual AGN system, the broad H{\beta} emission line is mainly contributed by the primary black hole (traced by the blue component) while the broad H{\beta} component of the secondary black hole (traced by the red component) is hard to be separated out considering a resolution of ~2000 of SDSS spectra or it is totally obscured by the dusty torus.

Self-gravitating accretion disks collapse to star-forming(SF) regions extending to the inner edge of the dusty torus in active galactic nuclei (AGNs). A full set of equations including feedback of star formation is given to describe the dynamics of the regions. We explore the role of supernovae explosion (SNexp), acting to excite turbulent viscosity, in the transportation of angular momentum in the regions within 1pc scale. We find that accretion disks with typical rates in AGNs can be driven by SNexp in the regions and metals are produced spontaneously. The present model predicts a metallicity--luminosity relationship consistent with that observed in AGNs. As relics of SF regions, a ring (or belt) consisting of old stars remains for every episode of supermassive black hole activity. We suggest that multiple stellar rings with random directions interact and form a nuclear star cluster after episodes driven by star formation.

It is found that \feii emission contributes significantly to the optical and ultraviolet spectra of most active galactic nuclei. The origin of the optical/UV \feii emission is still a question open to debate. The variability of \feii would give clues to this origin. Using 7.5 yr spectroscopic monitoring data of one Palomer-Green (PG) quasi-stellar object (QSO), PG 1700+518, with strong optical \feii emission, we obtain the light curves of the continuum \lv, \feii, the broad component of \hb, and the narrow component of \hb by the spectral decomposition. Through the interpolation cross-correlation method, we calculate the time lags for light curves of \feii, the total \hb, the broad component of \hb, and the narrow component of \hb with respect to the continuum light curve. We find that the \feii time lag in PG1700+518 is $209^{+100}_{-147}$ days, and the \hb time lag cannot be determined. Assuming that \feii and \hb emission regions follow the virial relation between the time lag and the FWHM for the \hb and \feii emission lines, we can derive that the \hb time lag is $148^{+72}_{-104}$ days. The \hb time lag calculated from the empirical luminosity--size relation is 222 days, which is consistent with our measured \feii time lag. Considering the optical \feii contribution, PG 1700+518 shares the same characteristic on the spectral slope variability as other 15 PG QSOs in our previous work, i.e., harder spectrum during brighter phase.

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov signals in radio detection are originated from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulses by making Fraunhofer approximation fails when the sizes of the elongated showers become comparable with the detection distances. We present a calculation method of Cherenkov pulses based on the finite-difference time-domain (FDTD) method, and attain a satisfying effeciency via the GPU- acceleration. Our method provides a straightforward way of the near field calculation, which would be important for ultra high energy particle showers, especailly the electromagnetic showers induced by the high energy leptons produced in the neutrino charge current interactions.

We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-emitting gas. Such gas would be accelerated away from the central object if the radiation force is able to act on the entire cloud. As had previously been deduced, infall requires that the clouds have large column densities so that a substantial amount of shielded gas is present. The critical column density required for infall to occur depends on L/L_Edd, establishing a link between Eigenvector 1 and the Fe II/Hbeta ratio. We see predominantly the shielded face of the infalling clouds rather than the symmetric distribution of emitters that has been assumed. The Fe II spectrum emitted by the shielded face is in good agreement with observations thus solving several long-standing mysteries in quasar emission lines.

Double-peaked [O III]5007, profiles in active galactic nuclei (AGNs) may provide evidence for the existence of dual AGNs, but a good diagnostic for selecting them is currently lacking. Starting from $\sim$ 7000 active galaxies in SDSS DR7, we assemble a sample of 87 type 2 AGNs with double-peaked [O III]5007, profiles. The nuclear obscuration in the type 2 AGNs allows us to determine redshifts of host galaxies through stellar absorption lines. We typically find that one peak is redshifted and another is blueshifted relative to the host galaxy. We find a strong correlation between the ratios of the shifts and the double peak fluxes. The correlation can be naturally explained by the Keplerian relation predicted by models of co-rotating dual AGNs. The current sample statistically favors that most of the [O III] double-peaked sources are dual AGNs and disfavors other explanations, such as rotating disk and outflows. These dual AGNs have a separation distance at $\sim 1$ kpc scale, showing an intermediate phase of merging systems. The appearance of dual AGNs is about $\sim 10^{-2}$, impacting on the current observational deficit of binary supermassive black holes with a probability of $\sim 10^{-4}$ (Boroson & Lauer).

Motivated by Genzel et al.'s observations of high-redshift star-forming galaxies, containing clumpy and turbulent rings or disks, we build a set of equations describing the dynamical evolution of gaseous disks with inclusion of star formation and its feedback. Transport of angular momentum is due to "turbulent" viscosity induced by supernova explosions in the star formation region. Analytical solutions of the equations are found for the initial cases of a gaseous ring and the integrated form for a gaseous disk, respectively. For a ring with enough low viscosity, it evolves in a slow processes of gaseous diffusion and star formation near the initial radius. For a high viscosity, the ring rapidly diffuses in the early phase. The diffusion drives the ring into a region with a low viscosity and start the second phase undergoing pile-up of gas at a radius following the decreased viscosity torque. The third is a sharply deceasing phase because of star formation consumption of gas and efficient transportation of gas inward forming a stellar disk. We apply the model to two $z\sim 2$ galaxies BX 482 and BzK 6004, and find that they are undergoing a decline in their star formation activity.

The growth of supermassive black holes (BHs) located at the centers of their host galaxies comes mainly from accretion of gas, but how to fuel them remains an outstanding unsolved problem in quasar evolution. This issue can be elucidated by quantifying the radiative efficiency parameter ($\eta$) as a function of redshift, which also provides constraints on the average spin of the BHs and its possible evolution with time. We derive a formalism to link $\eta$ with the luminosity density, BH mass density, and duty cycle of quasars, quantities we can estimate from existing quasar and galaxy survey data. We find that $\eta$ has a strong cosmological evolution: at z~2, $\eta \approx 0.3$, and by $z\approx 0$ it has decreased by an order of magnitude, to $\eta\approx 0.03$. We interpret this trend as evolution in BH spin, and we appeal to episodic, random accretion as the mechanism for reducing the spin. The observation that the fraction of radio-loud quasars decreases with increasing redshift is inconsistent with the popular notion that BH spin is a critical factor for generating strong radio jets. In agreement with previous studies, we show that the derived history of BH accretion closely follows the cosmic history of star formation, consistent with other evidence that BHs and their host galaxies coevolve.

Tracing the star formation history in circumnuclear regions (CNRs) is a key step towards understanding the starburst-AGN connection. However, bright nuclei outshining the entire host galaxy prevent the analysis of the stellar populations of CNRs around type-I AGNs. Obscuration of the nuclei by the central torus provides an unique opportunity to study the stellar populations of AGN host galaxies. We assemble a sample of 10, 848 type-II AGNs with a redshift range of $0.03\le z\le 0.08$ from the Sloan Digital Sky Survey's Data Release 4, and measure the mean specific star formation rates (SSFRs) over the past 100Myr in the central $\sim1-2$ kpc . We find a tight correlation between the Eddington ratio ($\lambda$) of the central black hole (BH) and the mean SSFR, strongly implying that supernova explosions (SNexp) play a role in the transportation of gas to galactic centers. We outline a model for this connection by accounting for the role of SNexp in the dynamics of CNRs. In our model, the viscosity of turbulence excited by SNexp is enhanced, and thus angular momentum can be efficiently transported, driving inflows towards galactic centers. Our model explains the observed relation $\lambda \propto \rm SSFR^{1.5-2.0}$, suggesting that AGN are triggered by SNexp in CNRs.

The second moment of the H$\beta$ emission line is calculated for 329 narrow line Seyfert 1 galaxies (NLS1s) selected from the Sloan Digital Sky Survey (SDSS), which is used to calculate the central supermassive black hole (SMBHs) mass of each. We find that the second moment depends strongly on the broader component of the H$\beta$ line profile. We find that for the NLS1s requiring two Gaussians to fit the H$\beta$ line the mean value of the SMBH mass from the H$\beta$ second moment is larger by about 0.50 dex than that from the full width at half maximum (FWHM). Using the gas velocity dispersion of the core/narrow component of \oiii $\lambda$ 5007 to estimate the stellar velocity dispersion, $\sigma_{*}$, the new mass makes NLS1s fall very close to the $\mbh - \sigma_{*}$ relation for normal AGNs. By using $\sigma_{*}$ measured directly from SDSS spectra with a simple stellar population synthesis method, we find that for NLS1s with mass lower than $10^7 \msun$, they fall only marginally below the $\mbh - \sigma_{*}$ relation considering the large scatter in the mass calculation.

We report on a systematic investigation of the Hbeta and Fe II emission lines in a sample of 568 quasars within z < 0.8 selected from the Sloan Digital Sky Survey. The conventional broad Hbeta emission line can be decomposed into two components--one with intermediate velocity width and another with very broad width. The velocity shift and equivalent width of the intermediate-width component do not correlate with those of the very broad component of Hbeta, but its velocity shift and width do resemble Fe II. Moreover, the width of the very broad component is roughly 2.5 times that of the intermediate-width component. These characteristics strongly suggest the existence of an intermediate-line region, whose kinematics seem to be dominated by infall, located at the outer portion of the broad-line region.

Broad Fe II emission is a prominent feature of the optical and ultraviolet spectra of quasars. We report on a systematical investigation of optical Fe II emission in a large sample of 4037 z < 0.8 quasars selected from the Sloan Digital Sky Survey. We have developed and tested a detailed line-fitting technique, taking into account the complex continuum and narrow and broad emission-line spectrum. Our primary goal is to quantify the velocity broadening and velocity shift of the Fe II spectrum in order to constrain the location of the Fe II-emitting region and its relation to the broad-line region. We find that the majority of quasars show Fe II emission that is redshifted, typically by ~ 400 km/s but up to 2000 km/s, with respect to the systemic velocity of the narrow-line region or of the conventional broad-line region as traced by the Hbeta line. Moreover, the line width of Fe II is significantly narrower than that of the broad component of Hbeta. We show that the magnitude of the Fe II redshift correlates inversely with the Eddington ratio, and that there is a tendency for sources with redshifted Fe II emission to show red asymmetry in the Hbeta line. These characteristics strongly suggest that Fe II originates from a location different from, and most likely exterior to, the region that produces most of Hbeta. The Fe II-emitting zone traces a portion of the broad-line region of intermediate velocities whose dynamics may be dominated by infall.

We investigate the $\mbh- \sigma_*$ relation for radio-loud quasars with redshift $z<0.83$ in Data Release 3 of the Sloan Digital Sky Survey (SDSS). The sample consists of 3772 quasars with better model of H$\beta$ and \oiii lines and available radio luminosity, including 306 radio-loud quasars, 3466 radio-quiet quasars with measured radio luminosity or upper-limit of radio luminosity (181 radio-quiet quasars with measured radio luminosity). The virial supermassive black hole mass (\mbh) is calculated from the broad \hb line, the host stellar velocity dispersion ($\sigma_*$) is traced by the core \oiii gaseous velocity dispersion, and the radio luminosity and the radio loudness are derived from the FIRST catalog. Our results are follows: (1) For radio-quiet quasars, we confirm that there is no obvious deviation from the $\mbh- \sigma_*$ relation defined in inactive galaxies when \mbh uncertainties and luminosity bias are concerned. (2) We find that radio-loud quasars deviate much from the $\mbh- \sigma_*$ relation respect to that for radio-quiet quasars. This deviation is only partly due to the possible cosmology evolution of the $\mbh- \sigma_*$ relation and the luminosity bias. (3) The radio luminosity is proportional to $\mbh^{1.28^{+0.23}_{-0.16}}(\lb/\ledd)^{1.29^{+0.31}_{-0.24}}$ for radio-quiet quasars and $\mbh^{3.10^{+0.60}_{-0.70}}(\lb/\ledd)^{4.18^{+1.40}_{-1.10}}$ for radio-loud quasars. The weaker correlation of the radio luminosity dependence upon the mass and the Eddington ratio for radio-loud quasars shows that other physical effects would account for their radio luminosities, such as the black hole spin.

This study reports pulse variation analysis results for the forth discovered accretion-powered millisecond pulsar XTE J1807-294 during its 2003 outburst observed by {\it Rossi X-ray Timing Explorer}. The pulsation is significantly detected only in the first $\sim$90d out of $\sim$150d observations. The pulse phase variation is too complex to be described as an orbital motion plus a simple polynomial model. The precise orbital parameters with $P_{orb}=40.073601(8)$ min and ${\it a_x}\sin {\it i}=4.823(5)$ lt-ms were obtained after applying the trend removal to the daily observed 150s segments pulse phases folded with a constant spin frequency without Keplerian orbit included. The binary barycenter corrected pulse phases show smooth evolution and clear negative phase shifts coincident with the flares seen on the light curve and the enhancements of fractional pulse amplitude. The non-flare pulse phases for the first $\sim$60d data are well described as a fourth order polynomial implying that the neutron star was spun-up during the first $\sim$60d with a rate $\dot \nu=(1.7\pm0.3) \times 10^{-13}$ Hz/s at the beginning of the outburst. Significant soft phase lags up to $\sim$500 $\mu s$ ($\sim$10% cycle) between 2 to 20 keV were detected for the nonflare pulse phases. We conclude that the anomalous phase shifts are unlikely due to the accretion torque but could result from the ``hot spot'' moving on the surface of neutron star.

The ultra-compact Low Mass X-ray Binary (LMXB) X1916-053, composed of a neutron star and a semi-degenerated white dwarf, exhibits periodic X-ray dips with variable width and depth. We have developed new methods to parameterize the dip to systematically study its variations. This helps to further understand binary and accretion disk behaviors. The RXTE 1998 observations clearly show a 4.87d periodic variation of the dip width. This is probably due to the nodal precession of the accretion disk, although there are no significant sidebands in the spectrum from the epoch folding search. From the negative superhump model (Larwood et. al. 1996), the mass ratio can be estimated as q = 0.045. Combined with more than 24 years of historical data, we found an orbital period derivative of $\dot{P}_{orb}/P_{orb}=(1.62 \pm 0.48)\times 10^{-7} yr^{-1}$ and established a quadratic ephemeris for the X-ray dips. The period derivative seems inconsistent with the prediction of the standard model of binary orbital evolution proposed by Rappaport et. al. (1987). On the other hand, the radiation-driven model (Tavani et. al. 1991) may properly interpret the period derivative even though the large mass outflow predicted by this model has never been observed in this system. With the best ephemeris, we obtained that the standard deviation of primary dips are smaller than that of secondary dips. This means that the primary dips are more stable than the secondary dips. Thus, we conclude that the primary dips of X1916-053 occur from the bulge at the rim instead of the ring of the disk proposed by Frank et. al. (1987).

Episodic activity of quasars is driving growth of supermassive black holes (SMBHs) via accretion of baryon gas. In this Letter, we develop a simple method to analyse the duty cycle of quasars up to redshift $z\sim 6$ universe from luminosity functions (LFs). We find that the duty cycle below redshift $z\sim 2$ follows the cosmic history of star formation rate (SFR) density. Beyond $z\sim 2$, the evolutionary trends of the duty cycle are just opposite to that of the cosmic SFR density history, implying the role of feedback from black hole activity. With the duty cycle, we get the net lifetime of quasars ($z\le 5$) about $\sim 10^9$yrs. Based on the local SMBHs, the mean mass of SMBHs is obtained at any redshifts and their seeds are of $10^5\sunm$ at the reionization epoch ($z_{\rm re}$) of the universe through the conservation of the black hole number density in comoving frame. We find that primordial black holes ($\sim 10^3\sunm$) are able to grow up to the seeds via a moderate super-Eddington accretion of $\sim 30$ times of the critical rate from $z=24$ to $z_{\rm re}$. Highly super-Eddington accretion onto the primordials is not necessary.

Feedback from black hole activity is widely believed to play a key role in regulating star formation and black hole growth. A long-standing issue is the relation between the star formation and fueling the supermassive black holes in active galactic nuclei (AGNs). We compile a sample of 57 Seyfert galaxies to tackle this issue. We estimate the surface densities of gas and star formation rates in circumnuclear regions (CNRs). Comparing with the well-known Kennicutt-Schmidt (K-S) law, we find that the star formation rates in CNRs of most Seyfert galaxies are suppressed in this sample. Feedback is suggested to explain the suppressed star formation rates.

Seed black holes formed in the collapse of population III stars have been invoked to explain the presence of supermassive black holes at high redshift. It has been suggested that a seed black hole can grow up to $10^{5\sim 6}\sunm$ through highly super-Eddington accretion for a period of $\sim 10^{6\sim 7}$ yr between redshift $z=20\sim 24$. We studied the feedback of radiation pressure, Compton heating and outflow during the seed black hole growth. It is found that its surrounding medium fueled to the seed hole is greatly heated by Compton heating. For a super-critical accretion onto a $10^3\sunm$ seed hole, a Compton sphere (with a temperature $\sim 10^6$K) forms in a timescale of $1.6\times 10^3$yr so that the hole is only supplied by a rate of $10^{-3}$ Eddington limit from the Compton sphere. Beyond the Compton sphere, the kinetic feedback of the strong outflow heats the medium at large distance, this leads to a dramatical decrease of the outer Bondi accretion onto the black hole and avoid the accumulation of the matter. The highly super-critical accretion will be rapidly halted by the strong feedback. The seed black holes hardly grow up at the very early universe unless the strong feedback can be avoided.

The highly disturbed hot gas in elliptical galaxies, as revealed in many {\em Chandra} X-ray images, implies a source of energy in the galactic nucleus. In some elliptical galaxies faint X-ray ``ghost'' cavities appear without corresponding radio lobes. It has been suggested that ghost cavities are caused by short-lived activity with a timescale of $\sim 10^3-10^4$ years, but this is difficult to understand within the popular paradigm of active galactic nuclei. We suggest an episode model for ghost cavities, invoking captures of red giant stars by the black hole located at the center of the elliptical galaxies at a typical rate of $10^{-5}$yr$^{-1}$ per galaxy. The accretion of tidally disrupted red giant stars onto the black hole powers activity in a timescale of a few years. The total energy channeled into the jet/outflow during the cooling time of the hot gas is $\sim 10^{56}$ erg, which is the typical work required to form the observed cavities. In this scenario, the faint cavities are produced by the feedback following accretion of the debris of the captured red giant stars onto the black holes. We apply the present model to several elliptical galaxies and find that it can explain the formation of the ghost cavities. This model can be tested in the future by comparisons between radio and X-ray observations.