The Hitomi results for the Perseus cluster have shown that accurate atomic models are essential to the success of X-ray spectroscopic missions, and just as important as knowledge on instrumental calibration and astrophysical modeling. Preparing the models requires a multifaceted approach, including theoretical calculations, laboratory measurements, and calibration using real observations. In a previous paper, we presented a calculation of the electron impact cross sections on the transitions forming the Fe-L complex. In the present work, we systematically test the calculation against cross sections of ions measured in an electron beam ion trap experiment. A two-dimensional analysis in the electron beam energies and X-ray photon energies is utilized to disentangle radiative channels following dielectronic recombination, direct electron-impact excitation, and resonant excitation processes in the experimental data. The data calibrated through laboratory measurements are further fed into global modeling of the Chandra grating spectrum of Capella. We investigate and compare the fit quality, as well as sensitivity of the derived physical parameters to the underlying atomic data and the astrophysical plasma modeling. We further list the potential areas of disagreement between the observation and the present calculations, which in turn calls for renewed efforts in theoretical calculations and targeted laboratory measurements.

We present the study on the X-ray emission for a sample of radio-detected quasars constructed from the cross-matches between SDSS, FIRST catalogs and XMM-Newton archives. A sample of radio-quiet SDSS quasars without FIRST radio detection is also assembled for comparison. We construct the optical and X-ray composite spectra normalized at rest frame $4215\,\rm \AA$ (or $2200\,\rm \AA$) for both radio-loud quasars (RLQs) and radio-quiet quasars (RQQs) at $z\le3.2$, with matched X-ray completeness of 19\%, redshift and optical luminosity. While the optical composite spectrum of RLQs is similar to that of RQQs, we find that RLQs have higher X-ray composite spectrum than RQQs, consistent with previous studies in the literature. By dividing the radio-detected quasars into radio loudness bins, we find the X-ray composite spectra are generally higher with increasing radio loudness. Moreover, a significant correlation is found between the optical-to-X-ray spectral index and radio loudness, and there is a unified multi-correlation between the radio, X-ray luminosities and radio loudness in radio-detected quasars. These results could be possibly explained with the corona-jet model, in which the corona and jet are directly related.

The origin and evolution of spectral curvature in blazar spectral energy distribution (SED) is still unclear. Since the observed SED curvature is related to an intrinsic curvature in emitting electron energy distribution (EED), we study this question by employing a log-parabolic EED with a curvature parameter and peak energy to model the quasi-simultaneous broadband SEDs of selected blazars in Fermi-LAT Bright AGN Sample (LBAS) using synchrotron and inverse Compton (IC) processes. We find that log-parabolic IC model can successfully explain the emission in all blazars in our sample. On average, FSRQs have higher magnetic field, Doppler factor, and curvature than BL Lac objects. The BL Lac objects show an anticorrelation between the curvature parameter of the EED and its peak energy, which is a signature of stochastic acceleration. FSRQs do not manifest such correlation and rather show a mild positive relationship between these parameters. This suggests that the evolution of spectral curvature in the BL Lac objects is dominated by a strong stochastic acceleration component, whereas the curvature in FSRQs evolves in a cooling dominated regime due to an additional external Compton (EC) component. The strong cooling in FSRQs not only restricts the electron peak energy but also adds extra curvature to the high energy tail of emitting EED. Since the curvature decreases from FSRQs toward high peak BL Lac objects (HBLs), opposite to peak energy, the curvature parameter can be considered a third parameter of the blazar sequence in addition to peak frequency and luminosity.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We perform a systematic search for high-redshift ($z >$ 1.5) extreme variability quasars (EVQs) using repeat spectra from the Sixteenth Data Release of Sloan Digital Sky Survey, which provides a baseline spanning up to $\sim$18 yrs in the observed frame. We compile a sample of 348 EVQs with a maximum continuum variability at rest frame 1450 Angstrom of more than 100% (i.e., $\delta$V $\equiv$ (Max$-$Min)/Mean $>$1). Among those EVQs, we discover 23 new Changing-Look Active Galactic Nuclei (CLAGNs) with CL behaviors in different UV lines, e.g., CIV, CIII], SiIV, and MgII. Importantly, spurious CLAGNs caused by SDSS problematic spectral flux calibration, e.g., fiber drop issue, are rejected. Compared with normal EVQs, we find that CLAGNs have systematically lower continuum luminosity and broader emission line widths, hence lower Eddington ratios. Statistical tests of EVQ/CLAGN properties suggest that CLAGNs are likely to be a subset of EVQs with less efficient accretion rather than a distinctive AGN population. Moreover, we find that the difference composite spectra (i.e., bright$-$faint) of EVQs and CLAGNs are both similar to those of normal quasars with a power-law slope of $\alpha_{\rm \lambda,diff}\sim-2$, further supporting that they belong to the same population with analogous intrinsic variability mechanisms, e.g., rapid accretion rate change with a possible state transition. Finally, we also confirm the anti-breathing of CIV (i.e., line width increases as luminosity increases) in EVQs, and find that in addition to $\sim$ 0.4 dex systematic uncertainty in single-epoch CIV virial black hole mass estimates, an extra scatter of $\sim$ 0.3 dex will be introduced by extreme variability.

We present the results of new radio observations carried out with the Karl G. Jansky Very Large Array C-configuration at 5.5 GHz for a sample of southern narrow-line Seyfert 1 galaxies (NLS1s). This work increases the number of known radio-detected NLS1s in the southern hemisphere, and confirms that the radio emission of NLS1s is mainly concentrated in a central region at kpc-scale and only a few sources show diffuse emission. In RQ NLS1s, the radio luminosity tends to be higher in steep-spectrum sources and be lower in flat-spectrum sources, which is opposite to RL NLS1s. This may be because the radio emission of steep NLS1s is dominated by misaligned jets, AGN-driven outflows, or star formation superposing on a compact core. Instead the radio emission of flat NLS1s may be produced by a central core which has not yet developed radio jets and outflows. We discover new NLS1s harboring kpc-scale radio jets and confirm that a powerful jet does not require a large-mass black hole to be generated. We also find sources dominated by star formation. These NLS1s could be new candidates in investigating the radio emission of different mechanisms.

In this work, we investigated the X-ray emission for a sample of young radio AGNs by combining their data from Chandra/XMM-Newton and at other wavebands. We find strong correlations between the X-ray luminosity $L_{\rm X}$ in 2$-$10 keV and the radio luminosities $L_{\rm R}$ at 5 GHz for VLBI radio-core, VLA radio-core and FIRST component, indicating that both pc- and kpc-scale radio emission strongly correlate with X-ray emission in these sources. We find approximately linear dependence of radio on X-ray luminosity in the sources with radiative efficient accretion flows (i.e., the Eddington ratio $R_{\rm edd} \gtrsim 10^{-3} $) with b $\sim$ 1 ($L_{\rm R}$ $\propto$ $L_{\rm X} ^{b}$) and $\xi_{\rm RX}$ $\sim$ 1 in fundamental plane using VLBI data, where the dependence is consistent with the re-analysed result on the previous study in \cite{2016ApJ...818..185F} at $R_{\rm edd} \gtrsim 10^{-3}$, however is significantly deviated from the theoretical prediction of accretion flow as the origin of X-ray emission. In contrast to radio-quiet quasars, there is no significant correlation between $\Gamma$ and Eddington ratio. Our results seem to indicate that the X-ray emission of high-accreting young radio AGNs may be from jet. We constructed the SEDs for 18 sources (most are in radiative efficient accretion) including 9 galaxies and 9 quasars with high-quality X-ray data, and find that the X-ray emission of most quasars is more luminous than that of normal radio-quiet quasars. This is clearly seen from the quasar composite SED, of which the X-ray emission is apparently higher than that of radio-quiet quasars, likely supporting the jet-related X-ray emission in young radio AGNs. The scenario that the X-ray emission is from self-synchrotron Compton (SSC) is discussed.

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

ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.

The reason for the difference in the composite X-ray spectrum between radio-loud quasars (RLQs) and radio-quiet quasars (RQQs) is still unclear. To study this difference, we built a new composite X-ray spectrum of RLQs by using Chandra X-ray data and Sloan Digital Sky Survey (SDSS) optical data for the sample of 3CRR quasars. We find the X-ray spectra of all 3CRR quasars except for 3C 351 have no soft X-ray excess and can be fitted with an absorbed power-law model well. Our composite X-ray spectrum is similar to that of \cite{ShangEtal2011apjs} for RLQs, showing higher hard X-ray and lower soft X-ray flux than the composite X-ray spectrum of RQQs. Most blazar-like 3CRR quasars have higher X-ray flux than the median composite X-ray spectrum, which could be related to the contribution of beamed jet emission at X-ray band. From the literature, we find that nineteen 3CRR quasars have extended X-ray emission related to radio jets, indicating inevitable contribution of jets at X-ray band. In contrast to RQQs, the X-ray photon index of 3CRR quasars does not correlate with the Eddington ratio. Our results suggest that the jet emission at X-ray band in RLQs could be related with the difference of composite X-ray spectrum between RLQs and RQQs.

We continuously monitored the blazar S5 0716+714 in the optical $g$, $r$ and $i$ bands from Nov. 10, 2017 to Jun. 06, 2019. The total number of observations is 201 nights including 26973 data points. This is a very large quasi-simultaneous multicolor sample for the blazar. The average time spans and time resolutions are 3.4 hours and 2.9 minutes per night, respectively. During the period of observations, the target source in the $r$ band brightens from $14^{\rm m}.16$ to $12^{\rm m}.29$ together with five prominent sub-flares, and then first becomes fainter to $14^{\rm m}.76$ and again brightens to $12^{\rm m}.94$ with seven prominent sub-flares. For the long-term variations, we find a strong flatter when brighter (FWB) trend at a low flux state and then a weak FWB trend at a higher flux state. A weak FWB trend at a low flux state and then a strong FWB trend at a higher flux state are also reported. Most of sub-flares show the strong FWB trends, except for two flares with a weak FWB trend. The particle acceleration and cooling mechanisms together with the superposition of different FWB-slopes from sub-flares are likely to explain the optical color behaviours. A scenario of bent jet is discussed.

The Gigahertz Peaked Spectrum (GPS) sources, Compact Steep Spectrum (CSS) radio sources, and High Frequency Peakers (HFP) radio sources are thought to be young radio AGNs, at the early stage of AGN evolution. We investigated the optical properties of the largest sample of 126 young radio AGNs based on the spectra in SDSS DR12. We find the black hole masses $M_{\rm BH}$ range from $10^{7.32}$ to $10^{9.84}$ $\rm M_{sun}$, and the Eddington ratios $R_{\rm edd}$ vary from $10^{-4.93}$ to $10^{0.37}$, suggesting that young radio AGNs have various accretion activities and not all are accreting at high accretion rate. Our young radio sources generally follow the evolutionary trend towards large-scale radio galaxies with increasing linear size and decreasing accretion rate in the radio power - linear size diagram. The radio properties of low luminosity young radio AGNs with low $R_{\rm edd}$ are discussed. The line width of [O III] $\lambda5007$ core ($\sigma_{\rm [O III]}$) is found to be a good surrogate of stellar velocity dispersion $\sigma_{\ast}$. The radio luminosity $L_{\rm 5GHz}$ correlates strongly with [O III] core luminosity $L_{\rm [O III]}$, suggesting that radio activity and accretion are closely related in young radio sources. We find one object that can be defined as narrow-line Seyfert 1 galaxies (NLS1s), representing a population of young AGNs both with young jet and early accretion activity. The optical variabilities of 15 quasars with multi-epoch spectroscopy were investigated. Our results show that the optical variability in young AGN quasars presents low variations ($\leqslant60\%$) similar to the normal radio-quiet quasars.

The spectral energy distribution of blazars around the synchrotron peak can be well described by the log-parabolic model that has three parameters: peak energy ($E_\textrm{p}$), peak luminosity ($L_\textrm{p}$) and the curvature parameter ($b$). It has been suggested that $E_\textrm{p}$ shows relations with $L_\textrm{p}$ and $b$ in several sources, which can be used to constrain the physical properties of the emitting region and/or acceleration processes of the emitting particles. We systematically study the $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations for 14 BL Lac objects using the 3-25~keV $RXTE$/PCA and 0.3-10~keV $Swift$/XRT data. Most objects (9/14) exhibit positive $E_\textrm{p}$-$L_\textrm{p}$ correlations, three sources show no correlation, and two sources display negative correlations. In addition, most targets (7/14) present no correlation between $E_\textrm{p}$ and 1$/b$, five sources pose negative correlations, and two sources demonstrate positive correlations. 1ES~1959+650 displays two different $E_\textrm{p}$-$L_\textrm{p}$ relations in 2002 and 2016. We also analyze $E_\textrm{p}$-$L_\textrm{p}$ and $E_\textrm{p}$-(1$/b$) relations during flares lasting for several days. The $E_\textrm{p}$-$L_\textrm{p}$ relation does not exhibit significant differences between flares, while the $E_\textrm{p}$-(1$/b$) relation varies from flare to flare. For the total sample, when $L_\textrm{p}$ < $\textrm{10}^\textrm{45}\ \textrm{erg}\ \textrm{s}^\textrm{-1}$, there seems to be a positive $E_\textrm{p}$-$L_\textrm{p}$ correlation. $L_\textrm{p}$ and the slope of $E_\textrm{p}$-$L_\textrm{p}$ relation present an anti-correlation, which indicates that the causes of spectral variations might be different between luminous and faint sources. $E_\textrm{p}$ shows a positive correlation with the black hole mass. We discuss the implications of these results.

Characterisation of the long-term variations of the broad line region in a luminous blazar, where Comptonization of broad line emission within a relativistic jet is the standard scenario for production of gamma-ray emission that dominates the spectral energy distribution. We analyse 10 years of optical spectroscopic data from the Steward Observatory for the blazar 3C 454.3, as well as gamma-ray data from the Fermi Large Area Telescope. The optical spectra are dominated by a highly variable non-thermal synchrotron continuum with a prominent Mg II broad emission line. The line flux is obtained by spectral decomposition including significant contribution from the Fe II pseudo-continuum. Three methods are used to characterise variations of the line flux: (1) stacking of the continuum-subtracted spectra, (2) subtracting the running mean light curves calculated for different time scales, (3) Discrete Correlation Function (DCF) for evaluating potential time delays. Despite very large variations of the gamma-ray and optical continua, the line flux changes only moderately (< 0.1 dex). The data suggest that the line flux responds to a dramatic change of the blazar activity from a very high state in 2010 to a deep low state in 2012. Two interpretations are possible: either the line flux is anti-correlated with the continuum, or increase in the line luminosity is delayed by ~600 days. If this time delay results from reverberation of poorly constrained accretion disk emission in both the broad line region (BLR) and the synchrotron emitting blazar zone within a relativistic jet, we would obtain natural estimates for the BLR radius [R_{BLR,MgII} >~ 0.28 pc] and for the supermassive black hole mass [M_SMBH ~ 8.5x10^8 M_sun]. We have not identified additional examples of short-term 'flares' of the line flux, besides the case of MJD 55512 reported by Le\'on-Tavares et al. (2013).

We present a study of a sample of 223 radio-loud quasars (up to redshift $<$0.3) in order to investigate their spectral properties. Twenty-six of these radio-loud quasars are identified as Flat Spectrum Radio Quasars (FSRQs) and fifty-four are identified as Steep Spectrum Radio Quasars (SSRQs) based on their radio spectral index. We study the [O III] line properties of these quasars to investigate the origin and properties of blue wings (shift of the profile towards lower wavelengths) and blue outliers (shift of the whole spectroscopic feature). Most of the quasars show blue wings with velocities up to 420 km $s^{-1}$. We find that around 17\% of the quasars show outliers with velocities spanning 419 to -315 km $s^{-1}$. Finally, we revisit the $\it M_{\rm BH} - \sigma$ relation of our sample using [S II]$\lambda$6716, 6731 and [O III] linewidths as surrogates for stellar velocity dispersions, $\sigma$, to investigate their location on the $\it M_{\rm BH} - \sigma$ relation for quiescent galaxies. Due to the strong blending of [S II] with $\rm H_{\alpha}$, we could estimate $\sigma_{[\rm SII]}$ of only 123 quasars. We find that the radio-loud quasars do not show a relationship between $\it M_{\rm BH}$ and $\sigma_{\rm [SII]/[OIII]}$ up to a redshift of 0.3, although they cluster around the local relation. We find an overall offset of 0.12$\pm$0.05 dex of our sample of radio-loud quasars from the $\it M_{\rm BH} - \sigma$ relation of quiescent galaxies. Quasars in our highest redshift bin (z=0.25-0.3) show a deviation of $\sim$0.33 $\pm$ 0.06 dex with respect to the local relation.

We have used an electron beam ion trap to measure electron-density-diagnostic line-intensity ratios for extreme ultraviolet lines from F XII, XIII, and XIV at wavelengths of 185-205 255-276 Angstroms. These ratios can be used as density diagnostics for astrophysical spectra and are especially relevant to solar physics. We found that density diagnostics using the Fe XIII 196.53/202.04 and the Fe XIV 264.79/274.21 and 270.52A/274.21 line ratios are reliable using the atomic data calculated with the Flexible Atomic Code. On the other hand, we found a large discrepancy between the FAC theory and experiment for the commonly used Fe XII (186.85 + 186.88)/195.12 line ratio. These FAC theory calculations give similar results to the data tabulated in CHIANTI, which are commonly used to analyze solar observations. Our results suggest that the discrepancies seen between solar coronal density measurements using the Fe XII (186.85 + 186.88)/195.12 and Fe XIII 196.54/202.04 line ratios are likely due to issues with the atomic calculations for Fe XII.

We present X-ray flux and spectral analyses of the three pointed Suzaku observations of the TeV high synchrotron peak blazar Mrk 421 taken throughout its complete operational duration. The observation taken on 5 May 2008 is, at 364.6 kiloseconds (i.e., 101.3 hours), the longest and most evenly sampled continuous observation of this source, or any blazar, in the X-ray energy 0.8 - 60 keV until now. We found large amplitude intra-day variability in all soft and hard bands in all the light curves. The discrete correction function analysis of the light curves in soft and hard bands peaks on zero lag, showing that the emission in hard and soft bands are cospatial and emitted from the same population of leptons. The hardness ratio plots imply that the source is more variable in the harder bands compared to the softer bands. The source is harder-when-brighter, following the general behavior of high synchrotron peak blazars. Power spectral densities of all three light curves are red noise dominated, with a range of power spectra slopes. If one assumes that the emission originates very close to the central super massive black hole, a crude estimate for its mass, of ~ 4 * 10^{8} M_{\odot}, can be made; but if the variability is due to perturbations arising there that are advected into the jet and are thus Doppler boosted, substantially higher masses are consistent with the quickest seen variations. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux and spectral variability.

We determined relative X-ray photon emission cross sections in Fe XVII ions that were mono-energetically excited in an electron beam ion trap. Line formation for the 3s (3s-2p) and 3d (3d-2p) transitions of interest proceeds through dielectronic recombination (DR), direct electron-impact excitation (DE), resonant excitation (RE), and radiative cascades. By reducing the electron-energy spread to a sixth of that of previous works and increasing counting statistics by three orders of magnitude, we account for hitherto unresolved contributions from DR and the little-studied RE process to the 3d transitions, and also for cascade population of the 3s line manifold through forbidden states. We found good agreement with state-of-the-art many-body perturbation theory (MBPT) and distorted-wave (DW) method for the 3s transition, while in the 3d transitions known discrepancies were confirmed. Our results show that DW calculations overestimate the 3d line emission due to DE by ~20%. Inclusion of electron-electron correlation effects through the MBPT method in the DE cross section calculations reduces this disagreement by ~11%. The remaining ~9% in 3d and ~11% in 3s/3d discrepancies are consistent with those found in previous laboratory measurements, solar, and astrophysical observations. Meanwhile, spectral models of opacity, temperature, and turbulence velocity should be adjusted to these experimental cross sections to optimize the accuracy of plasma diagnostics based on these bright soft X-ray lines of Fe XVII.

Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.

We present and analyze the variability of three Seyfert galaxies on intra-day timescales. We have analyzed in a uniform manner the 38 longest ( $>30$ ks) observations made for NGC 4051, MCG$-06-30-15$ and NGC 4151 by XMM-Newton between 2000 and 2015. The nuclei were quite active during essentially all of these observations and the overall X-ray fluxes (0.3-10 keV) varied by an order of magnitude. Most of the observations do appear to show characteristic timescales, estimated through their auto-correlation functions, ranging between $\sim 2.9$ ks and $\sim 45.3$ ks. The hard (2-10 keV) and soft (0.3-02 keV) bands are very well correlated but consideration of their hardness ratios show that the sources typically soften during flares. We also provide new estimates of the central black hole masses for these three AGN that support the hypothesis that Narrow Line Seyfert 1 galaxies have relatively small central black holes.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The stellar content of the intracluster light (ICL) provides unique insight into the hierarchical assembly process of galaxy clusters.However, the ICL is difficult to study due to its low surface brightness and large physical extent. We present optical spectra of three ICL regions in the Coma cluster, located between 100-180kpc from their nearest BCGs: NGC4889 and NGC4874. The mean surface brightness of the three ICL regions are {\mu}$_g$~25.3-26.2mag arcsec$^{-2}$. IFU spectroscopy with 13.5 hr on-source integration time were acquired as part of an ancillary program within the SDSS-IV MaNGA survey. We stacked the 127 individual fiber spectra in each IFU in order to achieve a 1{\sigma} limiting surface brightness of 27.9mag arcsec$^{-2}$, corresponding to a mean S/N in the optical of 21.6,9.6,and 11.6\AA$^{-1}$. We apply stellar population models to the stacked spectra, and measure the recession velocities, velocity dispersions ($\sigma$), stellar ages, and [Fe/H]. Our results show that the $\sigma$ of ICL regions are very high, indicating the stars are tracing the gravitational potential of the cluster, instead of any individual galaxy. The line-of-sight velocities of the three ICL regions are different from each other by ~700km/s, while the velocity of each region is similar to the closest BCG. This suggests that the ICL regions are associated with two distinct subclusters centered on NGC4889 and NGC4874.The stellar populations of these regions are old and metal poor, with ages of 7-12Gyr, and [Fe/H] of -0.8 to -0.6 dex. From the derived age and [Fe/H], the build-up of ICL in Coma is likely to be through the accretion of low mass galaxies or the tidal stripping of the outskirts of massive galaxies that have ended their star formation early on, instead of directly from major mergers of massive galaxies.

The radio-X-ray slope in the fundamental plane of radio-loud active galactic nuclei (AGNs) is found to be steeper compared with that of radio-quiet AGNs in previous works. In this work, we reinvestigate the fundamental plane in radio-loud AGNs by compiling a sample of 13 low-excitation radio galaxies (LERG) from the 3CR radio galaxies, for the reason that the accretion mode in LERG is believed to be a radiatively inefficient accretion flow. All the sources in our sample possess the data available at both the 5 GHz core radio luminosity detected by VLA/VLBI/VLBA and the core X-ray luminosity detected by Chandra/XMM-Newton. Surprisingly, we find the slope in the fundamental plane ($\log L_{\rm R}=0.52 \log L_{\rm X}+ 0.84 \log M_{\rm BH} + 10.84$) of LERG is well consistent with that reported by \citet{m2003}. However, the normalization is found to be shifted by about 0.7 dex, which can be due to the difference on magnetic field strength in different objects. A shallower slope of $L_{\rm R}-L_{\rm X}$ relation ($L_{\rm R}\sim L_{\rm X}^{0.63}$) is also given by our sample, which demonstrates that the X-ray emission in LERG may come from accretion disc instead of jets as suggested by previous works.

We found a possible ~ 1 hour quasi-periodic oscillation (QPO) in a ~ 55 ks X-ray observation of the narrow-line Seyfert 1 galaxy MCG--06--30--15 made with the XMM-Newton EPIC/pn detector in the energy range 0.3 -- 10 keV. We identify a total modulation of ~ 16% in the light curve and find a $\simeq$ 3670~s quasi-period using Lomb-Scargle periodogram (LSP) and weighted wavelet Z-transform (WWZ) techniques. Our analyses of eight light curves of MCG--06--30--15, indicated the possible presence of an oscillation during one of them. The LSP indicates a statistically significant ($\simeq$ 3$\sigma$) QPO detection. A WWZ analysis shows that the signal at this possible roughly 3670s period is present, and rather persistent, throughout the observation; however, a signal around 8735s is more persistent. We briefly discuss models that can produce X-ray QPOs with such periods in narrow line Seyfert 1 galaxies, as both other claimed QPO detections in this class of AGN had very similar periods.

A large population of ultra-diffuse galaxies (UDGs) was recently discovered in the Coma cluster. Here we present optical spectra of three such UDGs, DF7, DF44 and DF17, which have central surface brightnesses of $\mu_g \approx 24.4-25.1$ mag arcsec$^{-2}$. The spectra were acquired as part of an ancillary program within the SDSS-IV MaNGA Survey. We stacked 19 fibers in the central regions from larger integral field units (IFUs) per source. With over 13.5 hours of on-source integration we achieved a mean signal-to-noise ratio (S/N) in the optical of $9.5$\AA$^{-1}$, $7.9$\AA$^{-1}$ and $5.0$\AA$^{-1}$, respectively, for DF7, DF44 and DF17. Stellar population models applied to these spectra enable measurements of recession velocities, ages and metallicities. The recession velocities of DF7, DF44 and DF17 are $6599^{+40}_{-25}$km/s, $6402^{+41}_{-39}$km/s and $8315^{+43}_{-43}$km/s, spectroscopically confirming that all of them reside in the Coma cluster. The stellar populations of these three galaxies are old and metal-poor, with ages of $7.9^{+3.6}_{-2.5}$Gyr, $8.9^{+4.3}_{-3.3}$Gyr and $9.1^{+3.9}_{-5.5}$Gyr, and iron abundances of $\mathrm{[Fe/H]}$ $-1.0^{+0.3}_{-0.4}$, $-1.3^{+0.4}_{-0.4}$ and $-0.8^{+0.5}_{-0.5}$, respectively. Their stellar masses are $3$-$6\times10^8 M_\odot$. The UDGs in our sample are as old or older than galaxies at similar stellar mass or velocity dispersion (only DF44 has an independently measured dispersion). They all follow the well-established stellar mass$-$stellar metallicity relation, while DF44 lies below the velocity dispersion-metallicity relation. These results, combined with the fact that UDGs are unusually large for their stellar mass, suggest that stellar mass plays a more important role in setting stellar population properties for these galaxies than either size or surface brightness.

Brightest Cluster Galaxies (BCGs) in massive dark matter halos are shaped by complex merging processes. We present a detailed stellar population analysis in the central region of Abell 3827 at $z\sim0.1$, including five-nucleus galaxies involved in a BCG assembly. Based on deep spectroscopy from Multi Unit Spectroscopic Explorer (MUSE), we fit the optical spectra of 13 early-type galaxies (ETGs) in the central 70 kpc of the cluster. The stellar populations in the central $R=1$ kpc of these ETGs are old (6-10 Gyr). Their [Fe/H] increases with $\sigma_{\star}$ and stellar mass. More importantly, [$\alpha$/Fe] of galaxies close to the cluster center do not seem to depend on $\sigma_{\star}$ or stellar mass, indicating that the cluster center shapes the [$\alpha$/Fe]-$\sigma_{\star}$ and [$\alpha$/Fe]-$M_{\star}$ relations differently than other environments where [$\alpha$/Fe] is observed to increase with increasing $\sigma_{\star}$ or $M_{\star}$. Our results reveal the coordinated assembly of BCGs: their building blocks are different from the general low mass populations by their high [$\alpha$/Fe]. Massive galaxies thus grow by accreting preferentially high [$\alpha$/Fe] systems. The radial profiles also bear the imprint of the coordinated assembly. Their declining [Fe/H] and flat [$\alpha$/Fe] radial profiles confirm that the accreted systems have low metallicity and high [$\alpha$/Fe] stellar contents.

The relation between galaxies and dark matter halos reflects the combined effects of many distinct physical processes. Observations indicate that the $z=0$ stellar mass-halo mass (SMHM) relation has remarkably small scatter in stellar mass at fixed halo mass ($\lesssim$ 0.2 dex) with little dependence on halo mass. We investigate the origins of this scatter by combining N-body simulations with observational constraints on the SMHM relation. We find that at the group and cluster scale ($M_{\rm vir}>10^{14}{M_\odot}$) the scatter due purely to hierarchical assembly is $\approx0.16$ dex, which is comparable to recent direct observational estimates. At lower masses, mass buildup since $z\approx2$ is driven largely by in-situ growth. We include a model for the in-situ buildup of stellar mass and find that an intrinsic scatter in this growth channel of $0.2$ dex produces a relation between scatter and halo mass that is consistent with observations from $10^{12}{M_\odot}<M_{\rm vir}<10^{14.75}{M_\odot}$. The approximately constant scatter across a wide range of halo masses at $z=0$ thus appears to be a coincidence as it is determined largely by in-situ growth at low masses and by hierarchical assembly at high masses. These results indicate that the scatter in the SMHM relation can provide unique insight into the regularity of the galaxy formation process.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

We investigated the relationship between the accretion process and jet properties by ultilizing the VLBA and mid-infrared (MIR) data for a sample of 45 3CRR radio galaxies selected with a flux density at 178 MHz $>16.4$ Jy, 5 GHz VLA core flux density $\geq$ 7 mJy, and MIR observations. The pc-scale radio structure at 5 GHz are presented by using our VLBA observations for 21 sources in February, 2016, the analysis on the archival data for 16 objects, and directly taking the measurements for 8 radio galaxies available in literatures. The accretion mode is constrained from the Eddington ratio with a dividing value of 0.01, which is estimated from the MIR-based bolometric luminosity and the black hole masses. While most FRII radio galaxies have higher Eddington ratio than FRIs, we found that there is indeed no single correspondence between the FR morphology and accretion mode with eight FRIIs at low accretion and two FRIs at high accretion rate. There is a significant correlation between the VLBA core luminosity at 5 GHz and the Eddington ratio. Various morphologies are found in our sample, including core only, single-sided core-jet, and two-sided core-jet structures. We found that the higher accretion rate may be more likely related with the core-jet structure, thus more extended jet. These results imply that the higher accretion rates are likely able to produce more powerful jets. There is a strong correlation between the MIR luminosity at 15 $\mu$m and VLBA 5 GHz core luminosity, in favour of the tight relation between the accretion disk and jets. In our sample, the core brightness temperature ranges from $10^{9}$ to $10^{13.38}$ K with a median value of $10^{11.09}$ K indicating that systematically the beaming effect may not be significant....

We monitored BL Lacertae frequently during 2014 - 2016 when it was generally in a high state. We searched for intra-day variability for 43 nights using quasi-simultaneous measurements in the B, V, R, and I bands (totaling 143 light curves); the typical sampling interval was about eight minutes. On hour-like timescales, BL Lac exhibited significant variations during 13 nights in various optical bands. Significant spectral variations are seen during most of these nights such that the optical spectrum becomes bluer when brighter. The amplitude of variability is usually greater for longer observations but is lower when BL Lac is brighter. No evidence for periodicities or characteristic variability time-scales in the light curves was found. The color variations are mildly chromatic on long timescales.

The synchrotron peak of the X-ray bright High Energy Peaked Blazar (HBL) PKS 2155$-$304 occurs in the UV-EUV region and hence its X-ray emission (0.6--10 keV) lies mostly in the falling part of the synchrotron hump. We aim to study the X-ray emission of PKS 2155$-$304 during different intensity states in 2009$-$2014 using XMM$-$Newton satellite. We studied the spectral curvature of all of the observations to provide crucial information on the energy distribution of the non-thermal particles. Most of the observations show curvature or deviation from a single power-law and can be well modeled by a log parabola model. In some of the observations, we find spectral flattening after 6 keV. In order to find the possible origin of the X-ray excess, we built the Multi-band Spectral Energy distribution (SED). We find that the X-ray excess in PKS 2155--304 is difficult to fit in the one zone model but, could be easily reconciled in the spine/layer jet structure. The hard X-ray excess can be explained by the inverse Comptonization of the synchrotron photons (from the layer) by the spine electrons.

Blazars are classified into high, intermediate and low energy peaked sources based on the location of their synchrotron peak. This lies in infra-red/optical to ultra-violet bands for low and intermediate peaked blazars. The transition from synchrotron to inverse Compton emission falls in the X-ray bands for such sources. We present the spectral and timing analysis of 14 low and intermediate energy peaked blazars ob- served with XMMNewton spanning 31 epochs. Parametric fits to X-ray spectra helps constrain the possible location of transition from the high energy end of the syn- chrotron to the low energy end of the inverse Compton emission. In seven sources in our sample, we infer such a transition and constrain the break energy in the range 0.6 10 keV. The Lomb-Scargle periodogram is used to estimate the power spectral density (PSD) shape. It is well described by a power law in a majority of light curves, the index being flatter compared to general expectation from AGN, ranging here between 0.01 and 1.12, possibly due to short observation durations resulting in an absence of long term trends. A toy model involving synchrotron self-Compton (SSC) and exter- nal Compton (EC; disk, broad line region, torus) mechanisms are used to estimate magnetic field strength 6 0.03 - 0.88 G in sources displaying the energy break and infer a prominent EC contribution. The timescale for variability being shorter than synchrotron cooling implies steeper PSD slopes which are inferred in these sources.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

Narrow line Seyfert 1 galaxies (NLS1s) usually do not host relativistic jets and the detection of their $\gamma$-ray emissions is rarely reported. Evidences for strong relativistic beaming effect have been found in all $\gamma$-ray NLS1s to date and hence their jets are believed to be closely aligned. In this work we perform a systematic analysis of 7-year {\it Fermi} Large Area Telescope (LAT) data of 18 steep radio spectrum NLS1s (S-NLS1s) and report the {\it first} two detections of $\gamma$-ray emissions in this class, increasing the total number of known $\gamma$-ray NLS1s to 9. A strong and stable $\gamma$-ray signal (over 8$\sigma$) is produced by SDSS J144318.56+472556.7, and a significant $\gamma$-ray emission (4.5$\sigma$) of SDSS J130522.74+511640.2 has been detected in a monthly flare in November 2013. A thorough investigation of their multiwavelength properties from radio to X rays has been performed, finding no significant variability. Along with their spectral energy distributions and the radio morphology, these facts strongly suggest that these two sources harbor mis-aligned and possibly underdeveloped relativistic jets, which provides valuable targets to reveal the formation and evolution of relativistic jets under extreme circumstances.

This White Paper details the intentions and plans of the East Asian Very Long Baseline Interferometry (VLBI) community for pushing the frontiers of millimeter/submillimeter VLBI. To this end, we shall endeavor to actively promote coordinated efforts in the East Asia region. Our goal is to establish firm collaborations among the East Asia VLBI community in partnership with related institutes in North America and Europe and to expand existing global mm/submm VLBI arrays for (a) exploring the vicinity of black holes with an ultimate angular resolution down to 1~R$_{S}$ (Schwarzschild radius) and (b) investigating the dynamics of circumstellar gas in star-forming regions and late-type stars, and circumnuclear gas around active galactic nuclei (AGNs). In the first half of this White Paper, we highlight scientific accomplishments of the East Asia (EA) VLBI community. Various VLBI research results on M87, Sgr A*, blazars, narrow-line Seyfert~1 galaxies, and compact symmetric objects are described, and future visions of our VLBI science are briefly presented. Maser science of star formation, stellar evolution, and physics of accretion disks around AGNs are also discussed. A new vision for conducting multi-transition maser studies using mm/submm VLBI together with the Atacama Large Millimeter/Submillimeter Array (ALMA) is described. In the second half of this White Paper, we describe the EA community's vision for using mm/submm VLBI arrays in the framework or extended version of the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA). The accomplishment of the aforementioned goal will maximize the overall scientific outcomes of mm/submm VLBI in the world.

We studied the pc-scale core shift effect using radio light curves for three blazars, S5 0716+714, 3C 279 and BL Lacertae, which were monitored at five frequencies ($\nu$) between 4.8 GHz and 36.8 GHz using the University of Michigan Radio Astronomical Observatory (UMRAO), the Crimean Astrophysical Observatory (CrAO), and Metsahovi Radio Observatory for over 40 years. Flares were Gaussian fitted to derive time delays between observed frequencies for each flare ($\Delta t$), peak amplitude ($A$), and their half width. Using $A \propto \nu^{\alpha}$ we infer $\alpha$ in the range $-$16.67 to 2.41 and using $\Delta t \propto \nu^{1/k_r}$, we infer $k_r \sim 1$, employed in the context of equipartition between magnetic and kinetic energy density for parameter estimation. From the estimated core position offset ($\Omega_{r \nu}$) and the core radius ($r_{\rm core}$), we infer that opacity model may not be valid in all cases. The mean magnetic field strength at 1 pc ($B_1$) and at the core ($B_{\rm core}$), are in agreement with previous estimates. We apply the magnetically arrested disk model to estimate black hole spins in the range $0.15-0.9$ for these blazars, indicating that the model is consistent with expected accretion mode in such sources. The power law shaped power spectral density has slopes $-$1.3 to $-$2.3 and is interpreted in terms of multiple shocks or magnetic instabilities.

Multi-frequency Very Long Baseline Array (VLBA) polarimetry observation of the GHz-Peaked Spectrum (GPS) quasar OQ172 (J1445+0958) has been performed at 1.6, 2.2, 4.8, 8.3 and 15.3 GHz in 2005. Core-jet structures are detected in all bands with the jet strongly bent at about 3 mas from the core. The radio emission of the source is polarised at all five bands. We study the Faraday Rotation in the core and jet components at all five bands, and find good linear fits of Faraday Rotation in the core and jet components at 4.8 and 8.3 GHz. At these two frequencies, the Rotation Measure (RM) is $\sim 2000~\rm rad~m^{-2}$ in the core and $\sim 700~\rm rad~m^{-2}$ in the inner jet components and continues to decrease at the outer jet parts. We find that the depolarisation at 4.8 and 8.3 GHz might be caused by the internal medium in the source. We investigate consistency of the turnover spectra of VLBI components with the Synchrotron Self-Absorption (SSA) and Free-Free Absorption (FFA) models. Although these two models can not be easily distinguished due to the lack of low-frequency data, the physical parameters can be constrained for each model. We find that the large width of the $\rm [OIII]_{5007}$ line is likely caused by a jet interaction with a Narrow Line Region (NLR) medium. The jet bending, significant RM variations, Faraday depolarisation, spectral turnover, and broad line width of $\rm [OIII]_{5007}$ could be closely related, likely caused by the same nucleus medium, presumably NLR.

We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color index variabilities are $\bigtriangleup I=0^{\rm m}.67$, $\bigtriangleup R=0^{\rm m}.72$, $\bigtriangleup V=0^{\rm m}.68$, and $\bigtriangleup (V-R)=0^{\rm m}.25$ respectively. The largest clear IDV has an amplitude of 42% over just 5.8 minutes and the weakest detected IDV is 5.4% over 175 minutes. The BWB (bluer when brighter) chromatic trend is dominant for 3C 273 and appears at different flux levels on intraday timescales. The BWB trend exists for short-term timescales and intermediate-term timescales but different timescales have different correlations. There is no BWB trend for our whole time-series data sets. A significant anti-correlation between BWB trend and length of timescales is found. Combining with $V$-band data from previous works, we find a possible quasi-periodicity of $P=3918\pm1112$ days. The possible explanations for the observed variability, BWB chromatic trend and periodicity are discussed.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

The optical/ultraviolet (UV) variability of quasars has been discovered to be correlated with other quasar properties, such as luminosity, black hole mass and rest-frame wavelength. However,the origin of variability has been a puzzle so far. In this work, we upgrade the accretion disc model (Li & Cao 2008), which assumed the variability is caused by the change of global mass accretion rate, by constraining the disc size to match the viscous timescale of accretion disc to the variability timescale observed and by including the irradiation/X-ray reprocessing to make the emitted spectrum become steeper. We find this hybrid model can reproduce the observed bluer-when-brighter trend quite well, which is used to validate the theoretical model by several works recently. The traditional correlation between the variability amplitude and rest-frame wavelength can also be well fitted by our model. In addition, a weak positive correlation between variability amplitude and black hole mass is present, qualitatively consistent with recent observations.

We have collected near-infrared to X-ray data of 20 multi-epoch heavily reddened SDSS quasars to investigate the physical mechanism of reddening. Of these, J2317+0005 is found to be a UV cutoff quasar. Its continuum, which usually appears normal, decreases by a factor 3.5 at 3000{\AA}, compared to its more typical bright state during an interval of 23 days. During this sudden continuum cut-off, the broad emission line fluxes do not change, perhaps due to the large size of the Broad Line Region (BLR), r > 23 / (1+z) days. The UV continuum may have suffered a dramatic drop out. However, there are some difficulties with this explanation. Another possibility is that the intrinsic continuum did not change, but was temporarily blocked out, at least towards our line of sight. As indicated by X-ray observations, the continuum rapidly recovers after 42 days. A comparison of the bright state and dim states would imply an eclipse by a dusty cloud with a reddening curve having a remarkably sharp rise shortward of 3500{\AA}. Under the assumption of being eclipsed by a Keplerian dusty cloud, we characterized the cloud size with our observations, however, which is a little smaller than the 3000\AA\ continuum-emitting size inferred from accretion disk models. Therefore, we speculate this is due to a rapid outflow or inflow with a dusty cloud passing through our line-of-sight to the center.

In this paper, by cross-correlating an archive sample of 542 extragalactic radio sources with the \emph{Fermi}-LAT Third Source Catalog(3FGL), we have compiled a sample of 80 $\gamma$-ray sources and 462 non-\emph{Fermi} sources with available core dominance parameter($R_{CD}$), core and extended radio luminosity; all the parameters are directly measured or derived from available data in the literature. We found that $R_{CD}$ have significant correlations with radio core luminosities, $\gamma$-ray luminosity and $\gamma$-ray flux respectively; the \emph{Fermi} sources have on average higher $R_{CD}$ than non-\emph{Fermi} sources. These results indicate that the \emph{Fermi} sources should be more compact, and beaming effect should play a crucial role for the detection of $\gamma$-ray emission. Moreover, our results also show \emph{Fermi} sources have systematically larger radio flux than non-\emph{Fermi} sources at fixed $R_{CD}$, indicating larger intrinsic radio flux in \emph{Fermi} sources. These results show a strong connection between radio and $\gamma$-ray flux for the present sample and indicate that the non-\emph{Fermi} sources is likely due to low beaming effect, and/or the low intrinsic $\gamma$-ray flux, support a scenario in the literature: a co-spatial origin of the activity for the radio and $\gamma$-ray emission, suggesting that the origin of the seed photons for the high-energy $\gamma$-ray emission is within the jet.

We investigated the optical/ultraviolet (UV) color variations for a sample of 2169 quasars based on multi-epoch spectroscopy in the Sloan Digital Sky Survey (SDSS) data release seven (DR7) and data release nine (DR9). To correct the systematic difference between DR7 and DR9 due to the different instrumental setup, we produced a correction spectrum by using a sample of F-stars observed both in DR7 and DR9. The correction spectrum was then applied to quasars when comparing the spectra of DR7 with DR9. In each object, the color variation was explored by comparing the spectral index of the continuum power-law fit on the brightest spectrum with the faintest one, and also by the shape of their difference spectrum. In 1876 quasars with consistent color variations from two methods, we found that most sources (1755, $\sim 94\%$) show bluer-when-brighter (BWB) trend, and the redder-when-brighter (RWB) trend is only detected in 121 objects ($\sim 6\%$). The common BWB trend is supported by the bluer composite spectrum constructed from bright spectra than that from faint spectra, and also by the blue composite difference spectrum. The correction spectrum is proved to be highly reliable by comparing the composite spectrum from corrected DR9 and original DR7 spectra. Assuming that the optical/UV variability is triggered by fluctuations, RWB trend can likely be explained if the fluctuations occur firstly at outer disk region, and the inner disk region has not fully responded yet when the fluctuation being propagated inward. In contrast, the common BWB trend implies that the fluctuations are likely more often happening firstly in inner disk region.

We present the detection of compact radio structures of fourteen radio-loud narrow line Seyfert 1 (NLS1) galaxies from Very Long Baseline Array observations at 5 GHz, which were performed in 2013. While 50\% of the sources of our sample show a compact core only, the remaining 50\% exhibit a core-jet structure. The measured brightness temperatures of the cores range from $10^{8.4}$ to $10^{11.4}$ K with a median value of $10^{10.1}$ K, indicating that the radio emission is from non-thermal jets, and that, likely, most sources are not strongly beamed, then implying a lower jet speed in these radio-loud NLS1 galaxies. In combination with archival data taken at multiple frequencies, we find that seven sources show flat or even inverted radio spectra, while steep spectra are revealed in the remaining seven objects. Although all these sources are very radio-loud with $R > 100$, their jet properties are diverse, in terms of their milli-arcsecond (mas) scale (pc scale) morphology and their overall radio spectral shape. The evidence for slow jet speeds (i.e., less relativistic jets), in combination with the low kinetic/radio power, may offer an explanation for the compact VLBA radio structure in most sources. The mildly relativistic jets in these high accretion rate systems are consistent with a scenario, where jets are accelerated from the hot corona above the disk by the magnetic field and the radiation force of the accretion disk. Alternatively, a low jet bulk velocity can be explained by low spin in the Blandford-Znajek mechanism.

We observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source.

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

In a sample of 60 quasars selected from Sloan Digital Sky Survey (SDSS) with at least six-epoch spectroscopy, we investigate the variability of emission lines and continuum luminosity at various aspects. A strong anti-correlation between the variability and continuum luminosity at 2500 {\AA} is found for the sample, which is consistent with previous works. In individual sources, we find that half of the sample objects have a trend of being bluer-when-brighter trend (BWB), while the remaining half exhibit redder-when-brighter trend (RWB). Although the mechanism for RWB is unclear, the effects of host galaxy contribution due to seeing variations can not be completely ruled out. As expected from photoionization model, the positive correlations between the broad emission line and continuum luminosity are found in most individual sources, as well as for the whole sample. We confirm the Baldwin effect in most individual objects and the whole sample, while a negative Baldwin effect is also found in several quasars, which can be at least partly (if not all) due to the host galaxy contamination. We find positive correlations between the broad emission line luminosity and line width in most individual quasars, as well as the whole sample, implying a more variable line base than the line core.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

By cross-correlating an archival sample of 170 BL Lacs with 2 year \fermilat AGN sample, we have compiled a sample of 100 BL Lacs with \fermi detection (FBLs), and a sample of 70 non-\fermi BL Lacs (NFBLs). We compared various parameters of FBLs with those of NFBLs, including the redshift, the low frequency radio luminosity at 408 MHz ($L_{\rm 408MHz}$), the absolute magnitude of host galaxies ($M_{\rm host}$), the polarization fraction from NVSS survey ($P_{\rm NVSS}$), the observed arcsecond scale radio core flux at 5 GHz ($F_{\rm core}$) and jet Doppler factor; all the parameters are directly \textbf{measured} or derived from available data from literatures. We found that the Doppler factor is on average larger in FBLs than in NFBLs, and the $Fermi~ \gamma$-ray detection rate is higher in sources with higher Doppler factor. In contrast, there are no significant differences in terms of the intrinsic parameters of redshift, $ L_{\rm 408MHz}$, $ M_{\rm host}$ and $ P_{\rm NVSS}$. FBLs seem to have a higher probability of exhibiting measurable proper motion. These results strongly indicate a higher beaming effect in FBLs compared to NFBLs. The radio core flux is found to be strongly correlated with $\gamma$-ray flux, which remains after excluding the common dependence of the Doppler factor. At the fixed Doppler factor, FBLs have systematically larger radio core flux than NFBLs, implying lower $\gamma$-ray emission in NFBLs since the radio and $\gamma$-ray flux are significantly correlated. Our results indicate that the Doppler factor is an important parameter of $\gamma$-ray detection, the non-detection of $\gamma$-ray emission in NFBLs is likely due to low beaming effect, and/or low intrinsic $\gamma$-ray flux, and the gamma-rays are likely produced co-spatially with the arcsecond-scale radio core radiation and mainly through the SSC process.

In this letter, we present for the first time, the discovery of the disappearance of a narrow Mg II $\lambda\lambda2796,2803$ absorption system from the spectra of quasar SDSS J165501.31+260517.4 ($z_{\rm e}=1.8671$). This absorber is located at $z_{\rm abs} =1.7877$, and has a velocity offset of $8,423\rm ~km~s^{-1}$ with respect to the quasar. According to the velocity offset and the line variability, this narrow Mg II $\lambda\lambda2796,2803$ absorption system is likely intrinsic to the quasar. Since the corresponding UV continuum emission and the absorption lines of another narrow Mg II $\lambda\lambda2796,2803$ absorption system at $z_{\rm abs}=1.8656$ are very stable, we think that the disappearance of the absorption system is unlikely to be caused by the change in ionization of absorption gas. Instead, it likely arises from the motion of the absorption gas across the line of sight.

We have assembled a large sample of 5996 quasars at redshift 2.0=< z <= 2.4 (high-z) or 0.7=< z <= 1.1 (low-z) from SDSS data release nine and seven quasar catalogs. The spectral energy distribution (SED) of quasars were constructed by collecting WISE, UKIDSS, and GALEX photometric data in addition to SDSS, from which the IR luminosity at 1-7 \mu m and bolometric luminosity at 1100 \AA - 1 \mu m were calculated. A red tail is clearly seen in the distributions of the spectral index in 1100 \AA - 1 \mu m both for high-z and low-z sources, which is likely due to red or reddened quasars. The covering factor of dusty torus is estimated as the ratio of the IR luminosity to the bolometric luminosity. We found significant anti-correlations between the covering factor and bolometric luminosity, in both high-z and low-z quasars, however they follow different tracks. At overlapped bolometric luminosity, the covering factor of high-z quasars are systematically larger than those of low-z quasars, implying an evolution of covering factor with redshift.

Context. The steep-spectrum radio quasars (SSRQs) are powerful radio sources, with thermal emission from accretion disk and jet nonthermal emission likely both contributing in the Ultraviolet (UV)/optical luminosity, however the former may play a dominant role. While the UV/optical variability of SSRQs has been poorly studied, little is known on the mechanism of their variability. Aims. We investigate the mechanism of the UV/optical variability of SSRQs. Methods. A sample of eighteen SSRQs has been established in SDSS Stripe 82 region in our previous works, in which the flux and spectral variability have been studied. In this work, we construct the flux-flux diagram using SDSS u and i multi-epoch data for these eighteen SSRQs. The standard accretion disk model is used to fit the flux-flux variations, in order to explore the variability mechanism. Results. The model fit to flux-flux diagram are tuned with fixed black hole mass and varying accretion rate. We found that the flux-flux diagram of all our SSRQs can be qualitatively described by the standard accretion disk model with change of the accretion rate. Although non-thermal jet power-law emission can also qualitatively reproduce the variability, the reasonable accretion rates and black hole masses required to fit the flux-flux variations suggest that a disk emission with variable accretion rate is a plausible description of the data.

We compiled a sample of 60 quasars with spectroscopy on at least six epochs from the Sloan Digital Sky Survey (SDSS) to study the variabilities of the spectral shape, the continuum and the emission lines luminosity. In this paper, we present the results of SDSS J030639.57+000343.1. We found a strong anti-correlation between the continuum luminosity at $5100 \AA$ and the spectral index, implying a bluer-when-brighter trend. The luminosity of the broad $\rm H_\alpha$ line is proportion to the continuum luminosity at $5100 \AA$. Correspondingly, we did not find strong correlation between the equivalent width of broad $\rm H\alpha$ and the continuum luminosity, i.e. no baldwin effect of broad $\rm H\alpha$ in this source.

Minor mergers are thought to be important for the build-up and structural evolution of massive elliptical galaxies. In this work, we report the discovery of a system of four shell features in NGC 4889, one of the brightest members of the Coma cluster, using optical images taken with the Hubble Space Telescope and the Sloan Digital Sky Survey. The shells are well aligned with the major axis of the host and are likely to have been formed by the accretion of a small satellite galaxy. We have performed a detailed two-dimensional photometric decomposition of NGC 4889 and of the many overlapping nearby galaxies in its vicinity. This comprehensive model allows us not only to firmly detect the low-surface brightness shells, but, crucially, also to accurately measure their luminosities and colors. The shells are bluer than the underlying stars at the same radius in the main galaxy. We make use of the colors of the shells and the color-magnitude relation of the Coma cluster to infer the luminosity (or mass) of the progenitor galaxy. The shells in NGC 4889 appear to have been produced by the minor merger of a moderate-luminosity (M_I ~ -18.7 mag) disk (S0 or spiral) galaxy with a luminosity (mass) ratio of ~ 90:1 with respect to the primary galaxy. The novel methodology presented in this work can be exploited to decode the fossil record imprinted in the photometric substructure of other nearby early-type galaxies.

The spectral variability of a sample of 44 flat-spectrum radio quasars (FSRQs) and 18 steep-spectrum radio quasars (SSRQs) in the SDSS stripe 82 region is investigated. Twenty-five of 44 FSRQs show a bluer-when-brighter trend (BWB), while only one FSRQ shows a redder-when-brighter trend, which is in contrast to our previous results. Eight of 18 SSRQs display a BWB. We found an anti-correlation between the Eddington ratio and the variability amplitude in r band for SSRQs, which is similar to that in radio-quiet AGNs. This implies that the thermal emission from the accretion disk may be responsible for the variability in SSRQs. The spectral variability from the SDSS multi-epoch spectroscopy also shows BWB for several SSRQs, consistent with that from photometry.

We use a sample of radio-loud active galactic nuclei (AGNs) with measured black hole masses to explore the jet formation mechanisms in these sources. Based on the K\"{o}nigl's inhomogeneous jet model, the jet parameters, such as the bulk motion Lorentz factor, magnetic field strength, and electron density in the jet, can be estimated with the very long-baseline interferometry and X-ray data. We find a significant correlation between black hole mass and the bulk Lorentz factor of the jet components for this sample, while no significant correlation is present between the bulk Lorentz factor and the Eddington ratio. The massive black holes will be spun up through accretion, as the black holes acquire mass and angular momentum simultaneously through accretion. Recent investigation indeed suggested that most supermassive black holes in elliptical galaxies have on average higher spins than the black holes in spiral galaxies, where random, small accretion episodes (e.g., tidally disrupted stars, accretion of molecular clouds) might have played a more important role. If this is true, the correlation between black hole mass and the bulk Lorentz factor of the jet components found in this work implies that the motion velocity of the jet components is probably governed by the black hole spin. No correlation is found between the magnetic field strength at $10R_{\rm S}$ ($R_{\rm S}=2GM/c^2$ is the Schwarzschild radius) in the jets and the bulk Lorentz factor of the jet components for this sample. This is consistent with the black hole spin scenario, i.e., the faster moving jets are magnetically accelerated by the magnetic fields threading the horizon of more rapidly rotating black holes. The results imply that the Blandford-Znajek (BZ) mechanism may dominate over the Blandford-Payne (BP) mechanism for the jet acceleration at least in these radio-loud AGNs.

We present the results of EVN and MERLIN 5 GHz observations of nine ultra-high-energy synchrotron peak BL Lacs (UHBLs) selected as all BL Lacs with \textbf{log ($\nu_{\rm peak}/ \rm Hz)>20$} from Nieppola et al.. The radio structure was investigated for these sources, in combination with the available VLBA archive data. We found that the core-jet structure is detected in five sources, while four sources only have a compact core on pc scale. The core of all sources shows high brightness temperature (with mean and median values \textbf{log ($T_{\rm b} / {\rm K}) \sim11$}, which implies that the beaming effect likely present in all sources. When the multi-epoch VLBI data are available, we found no significant variations either for core or total flux density in two sources (2E 0414+0057 and EXO 0706.1+5913), and no evident proper motion in 2E 0414+0057, while the superluminal motion is likely detected in EXO 0706.1+5913. Our sources are found to be less compact than the typical HBLs in Giroletti et al, by comparing the ratio of the VLBI total flux to the core flux at arcsec scale. Combining all our results, we propose that the beaming effect might be present in the jets of UHBLs, however, it is likely weaker than that of typical HBLs. Moreover, we found that UHBLs could be less Doppler beamed versions of HBLs with similar jet power, by comparing the distribution of redshift, and radio luminosities. The results are in good consistence with the expectations from our previous work.

Context. Although a bluer-when-brighter trend is commonly observed in blazars, the opposite trend of redder-when-brighter has also been found in some blazars. Aims. We investigate the frequency of the redder-when-brighter trend in flat-spectrum radio quasars (FSRQs). Methods. We investigate the optical variability of 29 FSRQs in the SDSS Stripe 82 region using SDSS DR7 released multi-epoch data covering about nine years. We determined the spectral index by fitting a powerlaw to SDSS ugriz photometric data, and explored the relationship between the spectral index and source brightness. Results. For all FSRQs studied, we detect variations in r band flux of overall amplitude between 0.24 mag and 3.46 mag in different sources. Fourteen of 29 FSRQs display a bluer-when-brighter trend. However, only one source exhibits a redder-when-brighter trend, which implies that this behavior is rare in our FSRQ sample. In this source, the thermal emission from the accretion disk may be responsible for the redder-when-brighter trend.

While there are a lot of investigations on the optical variability of flat-spectrum radio quasars (FSRQs), not much work has been done on the optical variability of steep-spectrum radio quasars (SSRQs). We investigate the optical variability of 18 SSRQs and 15 FSRQs in the SDSS Stripe 82 region using SDSS DR7 released multi-epoch data covering about nine years. For all SSRQs studied, we detect variations in r band flux of overall amplitude between 0.22 mag and 0.92 mag in different sources. Eight of 18 SSRQs display a bluer-when-brighter (BWB) trend. In contrast, the variability amplitude of 15 FSRQs in r band ranges from 0.18 to 0.97 mag, and 11 of 15 sources show a BWB trend. We found an anti-correlation between the Eddington ratio and the variability amplitude in r band for SSRQs, which is similar to that in radio quiet AGNs. This implies that the thermal emission from accretion disk may be responsible for the variability in SSRQs, although the jet nonthermal emission cannot be excluded. In contrast, there is no correlation in FSRQs, which implies that the mechanisms of variability in FSRQs may be different from that in SSRQs. The optical continuum variability of radio loud broad absorption line quasars (BALQs) are investigated for the first time here on two sources with steep radio spectrum. Both radio loud BALQs show variations with amplitude of about 0.3 mag at r band. Their spectral variability all show a BWB trend. In a combined sample (18 SSRQs, and 44 FSRQs) of our sample with the FSRQs in Gu & Ai (2011), we found a trend of a broader line width of broad Mg II emission line with steeper radio spectral index. It implies a disc-like broad line region (BLR) geometry may present in these quasars. In these 62 quasars, we found a ~ 57% source percentage showing BWB trend in FSRQs, whereas it is ~ 44% in SSRQs.

Neutral hydrogen (HI) 21 cm absorption has been detected against more and more powerful radio jets. In this work, based on the Guppta et al. 2006a sample, we present our preliminary study of the correlations between the HI column density N(HI) and the jet power, N(HI) versus the low frequency luminosity at 408MHz, and N(HI) versus the radio luminosity at 1400MHz.

The optical variability of 29 flat spectrum radio quasars in SDSS Stripe 82 region are investigated by using DR7 released multi-epoch data. All FSRQs show variations with overall amplitude ranging from 0.24 mag to 3.46 mag in different sources. About half of FSRQs show a bluer-when-brighter trend, which is commonly observed for blazars. However, only one source shows a redder-when-brighter trend, which implies it is rare in FSRQs. In this source, the thermal emission may likely be responsible for the spectral behavior.

We examine the distribution of axis ratios of a large sample of disk galaxies hosting type 2 AGNs selected from the Sloan Digital Sky Survey and compare it with a well-defined control sample of non-active galaxies. We find them significantly different, where the type 2 AGNs show both an excess of edge-on objects and deficit of round objects. This systematical bias can not be explained by a nuclear obscurer oriented randomly with respect to the stellar disk. However, a nuclear obscurer coplanar with the stellar disk also does not fit the data very well. By assuming that the nuclear obscurer having an opening angle of ~60 degree, we find the observed axis ratio distribution can be nicely reproduced by a mean tilt angle of ~30 degree between the nuclear obscurer and the stellar disk.

We present the compact radio structure of three radio-loud narrow line Seyfert 1 galaxies from VLBA archive data at 2.3, 5 and 8.4 GHz. In RXS J16290+4007, the radio structure is mostly unresolved. The combination of compact radio structure, high brightness temperature and inverted spectrum between simultaneous 2.3 and 8.4 GHz, strongly favors jet relativistic beaming. Combining with the VLBI data at 1.6 and 8.4 GHz from literatures, we argued that RXS J16333+4718 may also harbor a relativistic jet, with resolved core-jet structure in 5 GHz. B3 1702+457 is clearly resolved with well defined jet component. The overall radio steep spectrum indicates that B3 1702+457 is likely a source optically defined as NLS1 with radio definition of compact steep spectrum sources. From these three sources, we found that radio loud NLS1s can be either intrinsically radio loud (e.g. B3 1702+457), or apparently radio loud due to jet beaming effect (e.g. RXS J16290+4007 and RXS J16333+4718).

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

We have collected all available spectra and photometric data from SDSS catalogue for bright AGNs complied from the first three months of the Fermi large area telescope all-sky survey. Based on the 106 high-confidence and 11 low-confidence association bright AGN list, the photometry data are collected from SDSS DR7 for 28 sources (12 BL Lacs and 16 FSRQs), two of which are low-confidence association bright AGNs. Among these 28 SDSS photometric sources, SDSS spectra are available for 20 sources (6 BL Lacs and 14 FSRQs). The black hole mass M_BH and the broad line region (BLR) luminosity were obtained for 14 FSRQs by measuring the line-width and strength of broad emission lines from SDSS spectra. The broad emission lines measurements of five FSRQs are presented for the first time in this work. The optical continuum emission of these 14 FSRQs is found to be likely dominated by the nonthermal jet emission through comparing the relationship between the broad Mg II line and continuum luminosity to that of radio quiet AGNs. The black hole mass of 14 FSRQs ranges from 10^8.2 to 10^9.9 solar mass, with most of sources larger than 10^9 solar mass. The Eddington ratio L_bol/L_Edd ranges from 10^-1.5 to ~ 1. This implies that the optically thin, geometrically thick accretion disk may exist in these FSRQs.

We find a significant anti-correlation between the hard X-ray photon index and the Eddington ratio L_Bol/L_Edd for a sample of Low-Ionization Nuclear Emission-line Regions (LINERs) and local Seyfert galaxies, compiled from literatures with Chandra or XMM-Newton observations. This result is in contrast with the positive correlation found in luminous active galactic nuclei (AGNs), while it is similar to that of X-ray binaries (XRBs) in low/hard state. Our result is qualitatively consistent with the spectra produced from advection dominated accretion flows (ADAFs). It implies that the X-ray emission of low-luminosity active galactic nuclei (LLAGNs) may originate from the Comptonization process in ADAF, and the accretion process in LLAGNs may be similar to that of XRBs in the low/hard state, which is different from that in luminous AGNs.

The relationship between the black hole mass and velocity dispersion indicated with [O III] line width is investigated for a sample of 87 flat-spectrum radio quasars (FSRQs) selected from SDSS DR3 quasar catalogue. We found the M_bh - sigma_[O III] relation is deviated from Tremaine et al. relation for nearby inactive galaxies, with a larger black hole mass at given velocity dispersion. There is no strong evidence of cosmology evolution in M_bh - sigma_[O III] relation up to z~0.8. A significant correlation between the [O III] luminosity and Broad Line Region (BLR) luminosity is found. When transferring the [O III] luminosity to Narrow Line Region (NLR) luminosity, the BLR luminosity is, on average, larger than NLR one by about one order of magnitude. We found a strong correlation between the synchrotron peak luminosity and NLR luminosity, which implies a tight relation between the jet physics and accretion process.

We review the accuracy of existing \ion{Fe}{17} X-ray line emission models by comparing them with an extensive analysis of \textit{Chandra} high energy transmission grating (HETG) observations of stellar coronae. We find significant discrepancies between most theoretical predictions and observations for at least some of the intensity ratios involving the six principal Fe XVII lines, 3C (15.01 {\AA}), 3D (15.26 {\AA}), 3E (15.45 {\AA}), 3F (16.78 {\AA}), 3G (17.05 {\AA}), and M2 (17.10 {\AA}). We suggest that the main problem of most previous theoretical studies to their inability to fully include electron correlation effects in the atomic structure calculations, while any deficiencies in the scattering approximation methods are of minor importance, regardless of it being close-coupling (CC) or distorted-wave (DW). An approximate method based on the many-body perturbation theory and DW approximation is proposed to include such correlation effects in the calculation of collisional excitation cross sections. The results are shown to agree with coronal observations and laboratory measurements better than most previous theories. Using the new atomic data, we then investigate the electron density sensitivity of the M2/3G intensity ratio and provide an improved density diagnostic tool for astrophysical observations.

We constructed a sample of 185 Flat Spectrum Radio Quasars (FSRQs) by cross-correlating the Shen et al.'s SDSS DR3 X-ray quasar sample with FIRST and GB6 radio catalogues. From the spectrum energy distribution (SED) constructed using multi-band (radio, UV, optical, Infrared and X-ray) data, we derived the synchrotron peak frequency and peak luminosity. The black hole mass and the broad line region (BLR) luminosity (then the bolometric luminosity) were obtained by measuring the line-width and strength of broad emission lines from SDSS spectra. We define a subsample of 118 FSRQs, of which the nonthermal jet emission is thought to be dominated over the thermal emission from accretion disk and host galaxy. For this subsample, we found 25 FSRQs having synchrotron peak frequency > 10^{15} Hz, which is higher than the typical value for FSRQs. While only a weak anti-correlation is found between the synchrotron peak frequency and peak luminosity, it becomes significant when combining with the Wu et al.'s sample of 170 BL Lac objects. At similar peak frequency, the peak luminosity of FSRQs with $\nupeak > 10^{15}$ Hz is systematically higher than that of BL Lac objects, with some FSRQs out of the range covered by BL Lac objects. Although high $\nupeak$ are found in some FSRQs, they do not reach the extreme value of BL Lacs. For the subsample of 118 FSRQs, we found significant correlations between the peak luminosity and black hole mass, the Eddington ratio, and the BLR luminosity, indicating that the jet physics may be tightly related with the accretion process.

The correlation between black hole masses and stellar velocity dispersions provides an efficient method to determine the masses of black holes in active galaxies. We obtained optical spectra of a Compact-Steep-Spectrum (CSS) galaxy 4C +29.70, using the Canada-France-Hawaii Telescope (CFHT) equipped with OSIS, in August 6, 2003. Several stellar absorption features, such as Mg I (5175\AA), Ca E band (5269\AA) and Na D (5890\AA), were detected in the spectra. The stellar velocity dispersion, $\sigma$, of the host galaxy, measured from absorption features is $\rm \approx 250 km s^{-1}$. If 4C +29.70 follows the $\rm M_{BH}-\sigma$ relation established for nearby galaxies, then its central black hole has a mass of $\rm \approx3.3\times10^{8}M_{\odot}$. In combination with the black hole masses of seven GPS galaxies in Snellen et al. (2003), we find that the average black hole mass of these eight young radio sources is smaller than that of the Bettoni et al. (2003) sample of extended radio galaxies. This may indicate that young radio sources are likely at the early evolutionary stage of radio galaxies, at which the central black holes may still undergo rapid growth. However, this needs further investigations.

Based on the K\"onigl's inhomogeneous jet model, we estimate the jet parameters, such as bulk Lorentz factor $\Gamma$, viewing angle $\theta$ and electron number density $n_{\rm e}$ from radio VLBI and X-ray data for a sample of active galactic nuclei (AGNs) assuming that the X-rays are from the jet rather than the intracluster gas. The bulk kinetic power of jets is then calculated using the derived jet parameters. We find a strong correlation between the total luminosity of broad emission lines and the bulk kinetic power of the jets. This result supports the scenario that the accretion process are tightly linked with the radio jets, though how the disk and jet are coupled is not revealed by present correlation analysis. Moreover, we find a significant correlation between the bulk kinetic power and radio extended luminosity. This implies that the emission from the radio lobes are closely related with the energy flux transported through jets from the central part of AGNs.

High resolution X-ray absorption spectroscopy is a powerful diagnostic tool for probing chemical and physical properties of the interstellar medium (ISM) at various phases. We present detections of K transition absorption lines from the low ionization ions of OI, OII, NeI, NeII, and NeIII, and the high ionization ones of OVI, OVII, OVIII, NeIX, and MgXI, as well as details of neutral absorption edges from Mg, Ne, and O in an unprecedented high quality spectrum of the low mass X-ray binary Cyg X-2. These absorption features trace the intervening interstellar medium which is indicated by the unshifted line centroids with respect to the rest frame wavelengths of the corresponding atomic transitions. We have measured the column densities of each ion. We complement these measurements with the radio HI and optical Halpha observations toward the same sight line and estimate the mean abundances of Ne, O, and Mg in the cool phase to Ne/H=0.84^{+0.13}_{-0.10}\times10^{-4}, O/H=3.83^{+0.48}_{-0.43}\times10^{-4}, and Mg/H=0.35^{+0.09}_{-0.11}\times10^{-4}, and O and Mg in the hot phase to O/H=5.81^{+1.30}_{-1.34}\times10^{-4} and Mg/H=0.33^{+0.09}_{-0.09}\times10^{-4}, respectively. These results indicate a mild depletion of oxygen into dust grains in the cool phase and little or no depletion of magnesium. We also find that absorption from highly ionized ions in the hot Galactic disk gas can account for most of the absorption observed toward the extragalactic sight lines like Mrk 421. The bulk of the observed OVI likely originates from the conductive interfaces between the cool and hot gases, from which a significant amount of NV and CIV emission is predicted.

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

We estimate the intrinsic luminosities and synchrotron peak frequency using the derived Doppler factor for a sample of 170 BL Lac objects, of which the synchrotron peak frequency are derived by fitting the SED constructed with the collected multi-band data from literatures. We find that the debeamed radio and optical core luminosities follow the same correlation found for FR I radio galaxies, which is in support of the unification of the BL Lac objects and the FR I galaxies based on orientation. For the debeamed luminosity at synchrotron peak frequency, we find a significant positive correlation between the luminosity and intrinsic synchrotron peak frequency. This implies that the more powerful sources may have the majority of jet emission at higher frequency. At synchrotron peak frequency, the intrinsic luminosity and black hole mass show strong positive correlation, while mild correlation is found in the case of jet power, indicating that the more powerful sources may have heavier black hole.

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

We present EVN observations of seven BL Lac objects selected from the RGB sample. To investigate the intrinsic radiation property of BL Lac objects, we estimated the Doppler factor with the VLA or MERLIN core and the total 408 MHz luminosity for a sample of 170 BL Lac objects. The intrinsic (comoving) synchrotron peak frequency was then calculated by using the estimated Doppler factor. Assuming a Lorentz factor of 5, the viewing angle of jets was constrained. The high-resolution VLBI images of seven sources all show a core-jet structure. We estimated the proper motions of three sources with the VLBI archive data, and find that the apparent speed increases with the distance of components to the core for all of them. In our BL Lacs sample, the Doppler factor of LBLs is systematically larger than that of IBLs and HBLs. We find a significant anti-correlation between the total 408 MHz luminosity and the intrinsic synchrotron peak frequency. However, the scatter is much larger than for the blazar sequence. Moreover, we find a significant positive correlation between the viewing angle and the intrinsic synchrotron peak frequency. The BL Lac objects show a continuous distribution on the viewing angle. While LBLs have a smaller viewing angle than that of IBLs and HBLs, IBLs are comparable to HBLs. We conclude that the intrinsic synchrotron peak frequency is not only related to the intrinsic radio power (though with a large scatter), but also to the viewing angle for the present sample.

We present a combined analysis of the X-ray emission of the Capella corona obtained with \textit{XMM-Newton} RGS, \textit{Chandra} HETGS, and LETGS. An improved atomic line database and a new differential emission measure (DEM) deconvolution method are developed for this purpose. Our new atomic database is based on the Astrophysical Plasma Emission Database and incorporates improved calculations of ionization equilibrium and line emissivities for L-shell ions of abundant elements using the Flexible Atomic Code. The new DEM deconvolution method uses a Markov Chain Monte-Carlo (MCMC) technique which differs from existing MCMC or $\chi^2$-fitting based methods. We analyze the advantages and disadvantages of each individual instrument in determining the DEM and elemental abundances. We conclude that results from either RGS or HETGS data alone are not robust enough due to their failure to constrain DEM in some temperature region or the lack of significant continuum emission in the wavelength band of the spectrometers, and that the combination of HETGS and RGS gives more stringent constraints on the DEM and abundance determinations. Using the LETGS data, we show that the recently discovered inconsistencies between the EUV and X-ray lines of Fe XVIII and XIX also exist in more highly charged iron ions, up to Fe XXIII, and that enhanced interstellar absorption due to partially ionized plasma along the Capella line of sight may explain some, but not all, of these discrepancies.

There is clear evidence that the $\gamma$-ray emission from active galactic nuclei (AGNs) is attributed to the inverse Compton scatterings in the relativistic blobs near the massive black holes. If the soft seed photons are from the regions outside the blobs, a linear relation between $(\nu F_{\nu,\gamma}/\nu F_{\nu, \rm syn} u^{*})^{1/(1+\alpha)}$ and Doppler factor $\delta$ is expected, where $\nu F_{\nu,\gamma}$ and $\nu F_{\nu, \rm syn}$ are monochromatic $\gamma$-ray and synchrotron fluxes, respectively, and $u^{*}$ is the energy density of soft seed photons \citep{D97}. We estimate the soft photon energy density in the relativistic blobs contributed by the broad line region (BLRs) in these $\gamma$-ray AGNs using their broad-line emission data. The Doppler factors $\delta$ are derived from their radio core and X-ray emission data, based on the assumption that the X-ray emission is produced through synchrotron self-Compton (SSC) scatterings. We find two nearly linear correlations: $(\nu F_{\nu,\gamma}/\nu F_{\rm opt} u^{*})^{1/(1+\alpha)} \propto \delta^{1.09}$, and $(\nu F_{\nu,\gamma}/\nu F_{\rm IR} u^{*})^{1/(1+\alpha)} \propto \delta^{1.20}$, which are roughly consistent with the linear correlation predicted by the theoretical model for external Compton scatterings. Our results imply that the soft seed photons are dominantly from the BLRs in these $\gamma$-ray AGNs.

We present the observational results of multi-colour optical monitoring of eight red blazars from 2003 September to 2004 February. The aim of our monitoring is to investigate the spectral variability as well as the flux variations at short and long time scales. The observations were carried out using the 1.0 m robotic telescope of Mt. Lemmon Optical Astronomy Observatory, in Arizona, USA, the 0.6 m telescope of Sobaeksan Optical Astronomy Observatory and the 1.8 m telescope of Bohyunsan Optical Astronomy Observatory, in the Republic of Korea. During the observations, all sources show strong flux variations with amplitudes of larger than 0.5 mag. Variations with amplitudes of over 1 mag are found in four sources. Intraday variations with amplitudes larger than 0.15 mag, and a rapid brightness increase with a rate of ~0.2 mag per day in four days, are detected in S5 0716+71. We investigate the relationship between the colour index and source brightness for each source. We find that two out of three FSRQs tend to be redder when they are brighter, and, conversely, all BL Lac objects tend to be bluer. In particular, we find a significant anti-correlation between the V-I colour index and R magnitude for 3C 454.3. This implies that the spectrum became steeper when the source was brighter, which is opposite to the common trend for blazars. In contrast, significant positive correlations are found in 3C 66A, S5 0716+71, and BL Lac. However, there are only very weak correlations for PKS 0735+17 and OJ 287. We propose that the different relative contributions of the thermal versus non-thermal radiation to the optical emission may be responsible for the different trends of the colour index with brightness in FSRQs and BL Lac objects.

We provide improved atomic calculation of wavelengths, oscillator strengths, and autoionization rates relevant to the $2\to 3$ inner-shell transitions of Fe~VI--XVI, the so-called Fe~M-shell unresolved transition array (UTA). A second order many-body perturbation theory is employed to obtain accurate transition wavelengths, which are systematically larger than previous theoretical results by 15--45~m{\AA}. For a few transitions of Fe~XVI and Fe~XV where laboratory measurements exist, our new wavelengths are accurate to within a few m{\AA}. Using these new calculations, the apparent discrepancy in the velocities between the Fe~M-shell UTA and other highly ionized absorption lines in the outflow of NGC 3783 disappears. The oscillator strengths in our new calculation agree well with the previous theoretical data, while the new autoionization rates are significantly larger, especially for lower charge states. We attribute this discrepancy to the missing autoionization channels in the previous calculation. The increased autoionization rates may slightly affect the column density analysis of the Fe~M-shell UTA for sources with high column density and very low turbulent broadening. The complete set of atomic data is provided as an electronic table.

The jet formation is thought to be closely connected with the mass of central supermassive black hole in Active Galactic Nuclei. The radio luminosity commonly used in investigating this issue is merely an indirect measure of the energy transported through the jets from the central engine, and severely Doppler boosted in core-dominated radio quasars. In this work, we investigate the relationship between the jet power and black hole mass, by estimating the jet power using extrapolated extended 151 MHz flux density from the VLA 5 GHz extended radio emission, for a sample of 146 radio loud quasars complied from literature. After removing the effect of relativistic beaming in the radio and optical emission, we find a significant intrinsic correlation between the jet power and black hole mass. It strongly implies that the jet power, so as jet formation, is closely connected with the black hole mass.To eliminate the beaming effect in the conventional radio loudness, we define a new radio loudness as the ratio of the radio extended luminosity to the optical luminosity estimated from the broad line luminosity.In a tentatively combined sample of radio quiet with our radio loud quasars, the apparent gap around the conventional radio loudness R=10 is not prominent for the new-defined radio loudness. In this combined sample, we find a significant correlation between the black hole mass and new-defined radio loudness.

We present high resolution laboratory spectra of K-shell X-ray lines from inner-shell excited and ionized ions of oxygen, obtained with a reflection grating spectrometer on the electron beam ion trap (EBIT-I) at the Lawrence Livermore National Laboratory. Only with a multi-ion model including all major atomic collisional and radiative processes, are we able to identify the observed K-shell transitions of oxygen ions from \ion{O}{3} to \ion{O}{6}. The wavelengths and associated errors for some of the strongest transitions are given, taking into account both the experimental and modeling uncertainties. The present data should be useful in identifying the absorption features present in astrophysical sources, such as active galactic nuclei and X-ray binaries. They are also useful in providing benchmarks for the testing of theoretical atomic structure calculations.

gamma Cas has long been famous for its unique hard X-ray characteristics. We report herein on a 53 ks Chandra HETGS observation of this target. An inspection of our spectrum shows that it is quite atypical for a massive star, with abnormally weak Fe XXV, XXVI lines, Ly-alpha lines of H-like species from Fe XVII, XXIII, XXIV, S XVI, Si XIV, Mg XII, Ne X, O VII, VIII, and N VII. Also, line ratios of the rif-triplet of for a few He-like ions XVII are consistent with the dominance of collisional atomic processes. Yet, the presence of Fe and Si fluorescence K features indicates that photoionization also occurs in nearby cold gas. The line profiles indicate a mean velocity at rest and a broadening of 500 km/s. A global fitting analysis of the line and continuum spectrum finds that there are 3-4 plasma emission components. The dominant hot (12 keV) component and has a Fe abundance of 0.22 solar. Some fraction of this component (10-30%) is heavily absorbed. The other 2-3 components, with temperatures 0.1, 0.4, 3 keV, are "warm," have a nearly solar composition, a lower column absorption, and are responsible for most other emission lines. The strength of the fluorescence features and the dual-column absorption model for the hot plasma component suggest the presence near the hot sites of a cold gas structure with a column density of 10^23 cm^-2. Since this value is consistent with theoretical estimates of the vertical disk column of this star, these attributes suggest that the X-rays originate near the star or disk. It is possible that the Fe anomaly in the hot component is related to the First Ionization Potential effect found in coronal structures around active cool stars. This would be yet another indication that the X-rays -rays are produced in the immediate vicinity of the Be star.

We present a detailed model of the discrete X-ray spectroscopic features expected from steady-state, low-density photoionized plasmas. We apply the Flexible Atomic Code (FAC) to calculate all of the necessary atomic data for the full range of ions relevant for the X-ray regime. These calculations have been incorporated into a simple model of a cone of ions irradiated by a point source located at its tip (now available as the XSPEC model PHOTOION). For each ionic species in the cone, photoionization is balanced by recombination and ensuing radiative cascades, and photoexcitation of resonance transitions is balanced by radiative decay. This simple model is useful for diagnosing X-ray emission mechanisms, determining photoionization/photoexcitation/recombination rates, fitting temperatures and ionic emission measures, and probing geometrical properties (covering factor/column densities/radial filling factor/velocity distributions) of absorbing/reemitting regions in photoionized plasmas. Such plasmas have already been observed in diverse astrophysical X-ray sources, including active galactic nuclei, X-ray binaries, cataclysmic variables, and stellar winds of early-type stars, and may also provide a significant contribution to the X-ray spectra of gamma-ray-burst afterglows and the intergalactic medium.

We report on the observation of OVII Doppler line broadening in a compact knot at the edge of SN 1006 detected with the Reflective Grating Spectrometer on board XMM-Newton. The observed line width of sigma = 3.4+/-0.5 eV at a line energy of 574 eV indicates an oxygen temperature of kT = 528+/-150 keV. Combined with the observed electron temperature of 1.5 keV the observed broadening is direct evidence for temperature non-equilibration in high Mach number shocks, and slow subsequent equilibration. The OVII line emission allows an accurate determination of the ionization state of the plasma, which is characterized by a relatively high forbidden line contribution, indicating log(n_e t) ~ 9.2.

High resolution spectroscopic observations of Seyfert galaxies with Chandra and XMM-Newton allows us to study the detailed ionization and thermal structures of the X-ray absorbing/emitting material in the circumnuclear environment. The vast improvement in the spectral resolving power by more than an order of magnitude enables us, for the first time, to unambiguously distinguish the dominant line emission mechanisms and to measure its dynamical properties as well. The X-ray band harbors spectral transitions from a wide range of ionization states, including valence-shell transitions in K-shell and L-shell ions from most cosmically abundant elements, as well as inner-shell transitions of iron and other mid-Z elements, which can be probed through absorption measurements. The X-ray spectrum, therefore, provides simultaneous velocity and column density constraints of highly ionized to only slightly ionized gas harbored in many of these systems. We summarize recent results that have emerged from observations of Seyfert 2 galaxies with the grating spectrometers onboard Chandra and XMM-Newton. We give particular emphasis to an empirical physical model that we have developed based on the observed spectra, and how it can be used for comparative studies with Seyfert 1 galaxies to test the AGN unification scenarios.

We present the first high-resolution, soft-X-ray spectrum of the prototypical Seyfert 2 galaxy, NGC 1068. This spectrum was obtained with the XMM-Newton Reflection Grating Spectrometer. Emission lines from H-like and He-like low-Z ions (from C to Si) and Fe-L-shell ions dominate the spectrum. Strong, narrow radiative recombination continua (RRC) for several ions are also present, implying that most of the observed soft-X-ray emission arises in low-temperature (few eV) plasma. This plasma is photoionized by the inferred nuclear continuum (obscured along our line of sight), as in the unified model of active galactic nuclei (AGN). We find excess emission (compared with pure recombination) in all resonance lines (np to 1s) up to the photoelectric edge, demonstrating the importance of photoexcitation as well. We introduce a simple model of a cone of plasma irradiated by the nuclear continuum; the line emission we observe along our line of sight perpendicular to the cone is produced through recombination/radiative cascade following photoionization and radiative decay following photoexcitation. A remarkably good fit is obtained to the H-like/He-like ionic line series, with inferred radial ionic column densities consistent with recent observations of warm absorbers in Seyfert 1 galaxies. Previous Chandra imaging revealed a large (extending out to 500 pc) ionization cone containing most of the X-ray flux, implying that the warm absorber in NGC 1068 is a large-scale outflow. To explain the ionic column densities, a broad, flat distribution in the logarithm of the ionization parameter ($\xi=L_X/n_e r^2$) is necessary, spanning $\log\xi=0$--3. This suggests either radially-stratified ionization zones or the existence of a broad density distribution (spanning a few orders of magnitude) at each radius.

We present high-resolution soft-X-ray spectra of the prototypical Seyfert 2 galaxy, NGC 1068, taken with XMM-Newton RGS and Chandra LETGS. Its rich emission-line spectrum is dominated by recombination in a warm plasma (bright, narrow radiative recombination continua provide the ``smoking gun''), which is photoionized by the inferred nuclear power-law continuum. Radiative decay following photoexcitation of resonant transitions is also significant. A self-consistent model of an irradiated cone of gas is capable of reproducing the hydrogenic/heliumlike ionic line series in detail. The radial ionic column densities we infer are consistent with absorption measurements (the warm absorber) in Seyfert 1 galaxies. This strongly suggests that the emission spectrum we observe from NGC 1068 emanates from its warm absorber. The observed extent of the ionization-cone/warm absorber in NGC 1068 of about 300 pc implies that a large fraction of the gas associated with generic warm absorbers may typically exist on the hundreds-of-parsec scale rather than much closer to the nucleus (e.g., less than a parsec). Spatially-resolved spectroscopy using the LETGS of two distinct emission regions yields two noticeably different spectra. We show that these differences are solely due to differing radial column densities. A fairly flat distribution in ionization parameter is necessary to explain the inferred radial ionic column densities of all spectra. This must primarily be due to a broad density distribution at each radius, spanning roughly 0.1-100 cm$^{-3}$. (Abridged)

The central black-hole masses of a sample of radio-loud quasars are estimated by using the data of $H_{\beta}$ line-width and the optical continuum luminosity. The vast majority of the quasars in this sample have black-hole masses larger than $10^{8} M_{\odot}$, while a few quasars may contain relatively smaller black-holes. We found a significant anti-correlation between the radio-loudness and the central black-hole mass. It might imply that the jet formation is governed by the black-hole mass.

We derive a semi-analytic solution for the structure of conduction-mediated transition layers above an X-ray illuminated accretion disk, and calculate explicitly the X-ray line radiation resulting from both resonance line scattering and radiative recombination in these layers. The vertical thermal structure of the illuminated disk is found to depend on the illuminating continuum: for a hard power law continuum, there are two stable phases connected by a single transition layer; while for a softer continuum, there may exist three stable phases connected by two separate transition layers, with an intermediate stable layer in between. We show that the structure can be written as a function of the electron scattering optical depth through these layers, which leads to unique predictions of the equivalent width of the resulting line radiation from both recombination cascades and resonance line scattering. We find that resonance line scattering plays an important role, especially for the case where there is no intermediate stable layer.

The low frequency radio luminosity is believed to be an indicator of jet power, while the optical/UV emission is probably from accretion discs in the nuclei of steep-spectrum radio quasars. We present a correlation between the ratio of radio to optical luminosities and the continuum spectral index in blue/UV bands, which might indicate that the continuum shape in blue/UV bands is related to the ratio of jet to accretion power. The results may imply that the spectra and structure of accretion discs are probably affected by the interactions between jets and discs.