We used two XMM-Newton and six Neutron Star Interior Composition Explorer (NICER) observations to investigate the fractional rms amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary 4U 1636-53. We studied, for the first time, the fractional rms amplitude of the mHz QPOs vs. energy in 4U 1636-53 down to 0.2 keV. We find that, as the energy increases from 0.2 keV to 3 keV, the rms amplitude of the mHz QPOs increases, different from the decreasing trend that has been previously observed above 3 keV. This finding has not yet been predicted by any current theoretical model, however, it provides an important observational feature to speculate whether a newly discovered mHz oscillation originates from the marginally stable nuclear burning process on the neutron star surface.

We report the discovery of year-scale X-ray variation in the nuclear region of the M87 by reanalyzing the 8 Chandra observations from 2007 to 2008. The X-ray spectra are fitted and decomposed into the disk and flaring components. This year-scale X-ray variability can be explained quite well by a simple clumpy accretion model. We conclude that the central super-massive black hole of the M87 was accreting a cloud of $\sim0.5$M$_\odot$ at that time.

Based on statistical analysis of synchrotron polarization intensity, we study the anisotropic properties of compressible magnetohydrodynamic (MHD) turbulence. The second-order normalized structure function, quadrupole ratio modulus and anisotropic coefficient are synergistically used to characterize the anisotropy of the polarization intensity. On the basis of pre-decomposition data cubes, we first explore the anisotropy of the polarization intensity in different turbulence regimes and find that the most significant anisotropy occurs in the sub-Alfv\'enic regime. Using post-decomposition data cubes in this regime, we then study the anisotropy of the polarization intensity from Alfv\'en, slow and fast modes. Statistics of polarization intensity from Alfv\'en and slow modes demonstrate the significant anisotropy while statistics of polarization intensity from fast modes show isotropic structures, which is consistent with the earlier results provided in Cho & Lazarian (2002). As a result, both quadrupole ratio modulus and anisotropic coefficient for polarization intensities can quantitatively recover the anisotropy of underlying compressible MHD turbulence. The synergistic use of the two methods helps enhance the reliability of the magnetic field measurement.

The carriers of the still (mostly) unidentified diffuse interstellar bands (DIBs) have been a long-standing mystery ever since their first discovery exactly 100 years ago. In recent years, the ubiquitous detection of a large number of DIBs in a wide range of Galactic and extragalactic environments has led to renewed interest in connecting the occurrence and properties of DIBs to the physical and chemical conditions of the interstellar clouds, with particular attention paid to whether the DIB strength is related to the shape of the interstellar extinction curve. To shed light on the nature and origin of the DIB carriers, we investigate the relation between the DIB strength and $R_V$, the total-to-selective extinction ratio, which characterizes how the extinction varies with wavelength (i.e., the shape of the extinction curve). We find that the DIB strength and $R_V$ are not related if we represent the strength of a DIB by its reddening-normalized equivalent width (EW), in contrast to the earlier finding of an anticorrelation in which the DIB strength is measured by the extinction-normalized EW. This raises a fundamental question about the appropriate normalization for the DIB EW.

Interstellar graphene could be present in the interstellar medium (ISM), resulting from the photochemical processing of polycyclic aromatic hydrocarbon (PAH) molecules and/or collisional fragmentation of graphitic particles. Indeed, by comparing the observed ultraviolet (UV) extinction and infrared (IR) emission of the diffuse ISM with that predicted for graphene, as much as $\sim$ 2% of the total interstellar carbon could have been locked up in graphene without violating the observational constraints. While the possible detection of planar C$_{24}$, a small piece of a graphene sheet, has been reported towards several Galactic and extragalactic planetary nebulae, graphene has not yet been detected in interstellar environments. In this work, we search for the characteristic IR features of C$_{24}$ at $\sim$ 6.6, 9.8, 20 $\mu m$ toward Sgr B2, a high-mass star formation region, and find a candidate target toward R.A. (J2000) = $267^{\circ}.05855$ and Decl. (J2000) = $-28^{\circ}.01479$ in Sgr B2 whose Spitzer/IRS spectra exhibit three bands peaking at $\sim$ 6.637, 9.853 and 20.050 $\mu m$ which appear to be coincident with that of C$_{24}$. Possible features of C$_{60}$ are also seen in this region. The candidate region is a warm dust environment heated by massive stars or star clusters, associated with a WISE spot (a tracer of star-formation activity), close to the HII region candidate IRAS 17450-2759, and is surrounded by seven young stellar object candidates within $\sim$ 5$^{\prime}$, suggesting that the creation and/or excitation of C$_{24}$ could be related to star formation activities.

We analysed an XMM-Newton plus a simultaneous Rossi X-ray Timing Explorer observation and a separate Suzaku observation of the neutron-star low-mass X-ray binary 4U 1728-34. We fitted the X-ray spectra with the self-consistent reflection model relxill. We found that the inclination angle of 4U 1728-34 is 49 degrees, consistent with the upper limit of 60 degrees deduced from the absence of eclipses or dips in this source. The inclination angle in the fit of the XMM-Newton/RXTE observation is larger than 85 degrees, which may be due to the possible calibration issues of the PN instrument in timing mode. We also found that the thermal emission from the accretion disc is not significant. This could be explained either by the relatively high column density of the interstellar medium along the line of sight to the source, which decreases the number of soft disc photons, or if most of the soft thermal photons from the disc are reprocessed in the corona. The ionisation parameter derived from the fits is larger than the value predicted in the framework of the standard reflection model, wherein the disc is irradiated by an X-ray source above the compact object. This inconsistency suggests that irradiation from the neutron star and its boundary layer may play an important role in the ionisation of the accretion disc, and hence in the reflection component in this source.

We investigate how to recover the spectral properties of underlying magnetohydrodynamic (MHD) turbulence using fluctuation statistics of synchrotron polarization radiation, based on the synthetic observations. Taking spatially coincident, separated, and compounded synchrotron emission and Faraday rotation regions into account, we extract the power spectrum of synchrotron polarization intensities integrated along the line of sight. Our results demonstrate that in the short wavelength range, the power spectra reflect fluctuation statistics of the perpendicular component of turbulent magnetic fields, and the spectra at long wavelengths reveal the fluctuation of the Faraday rotation density, which is a product of the parallel component of magnetic field and thermal electron density. We find that our numerical results (in the case of spatially coincident regions) are in agreement with the analytical prediction in Lazarian \& Pogosyan (2016), and this theoretical prediction is applicable to more complicated settings, i.e., the spatially separated and compounded regions. We simulate telescopic observations that incorporate the effects of telescope angular resolution and noise, and find that statistics of underlying MHD turbulence can be recovered successfully. We expect that the technique can be applied to a variety of astrophysical environments, with existing synchrotron data cubes and a large number of forthcoming data sets from such as the LOw-Frequency Array for Radio astronomy (LOFAR), the Square Kilometer Array (SKA) and the Five-hundred-meter Aperture Spherical radio Telescope (FAST).

This paper studies acceleration processes of background thermal electrons in X-ray binary jets via turbulent stochastic interactions and shock collisions. By considering turbulent magnetized jets mixed with fluctuation magnetic fields and ordered, large-scale one, and numerically solving the transport equation along the jet axis, we explore the influence of such as magnetic turbulence, electron injections, location of an acceleration region, and various cooling rates on acceleration efficiency. The results show that (1) the existence of the dominant turbulent magnetic fields in the jets is necessary to accelerate background thermal electrons to relativistic energies. (2) Acceleration rates of electrons depend on magnetohydrodynamic turbulence types, from which the turbulence type with a hard slope can accelerate electrons more effectively. (3) An effective acceleration region should be located at the distance $>10^3R_{\rm g}$ away from the central black hole ($R_{\rm g}$ being a gravitational radius). As a result of acceleration rates competing with various cooling rates, background thermal electrons obtain not only an increase in their energies but also their spectra are broadened beyond the given initial distribution to form a thermal-like distribution. (4) The acceleration mechanisms explored in this work can reasonably provide the electron maximum energy required for interpreting high-energy $\gamma$-ray observations from microquasars, but it needs to adopt some extreme parameters in order to predict a possible very high-energy $\gamma$-ray signal.

We study the properties of polarized radiation in turbulent magnetic fields from X-ray binary jets. These turbulent magnetic fields are composed of large- and small-scale configurations, which result in the polarized jitter radiation when the characteristic length of turbulence is less than the non-relativistic Larmor radius. On the contrary, the polarized synchrotron emission occurs, corresponding to a large-scale turbulent environment. We calculate the spectral energy distributions and the degree of polarization for a general microquasar. Numerical results show that turbulent magnetic field configurations can indeed provide a high degree of polarization, which does not mean that a uniform, large-scale magnetic field structure exists. The model is applied to investigate the properties of polarized radiation of black hole X-ray binary Cygnus X--1. Under the constraint of multiband observations of this source, our studies demonstrate that the model can explain the high polarization degree at MeV tail and predict the highly polarized properties at high-energy $\gamma$-ray region, and that the dominant small-scale turbulent magnetic field plays an important role for explaining the highly polarized observation at hard X-ray/soft $\gamma$-ray bands. This model can be tested by polarization observations of upcoming polarimeters at high-energy $\gamma$-ray bands.

A large number of interstellar absorption features at ~ 4000\AA\ -- 1.8 {\mu}m, known as the "diffuse interstellar bands" (DIBs), remains unidentified. Most recent works relate them to large polycyclic aromatic hydrocarbon (PAH) molecules or ultrasmall carbonaceous grains which are also thought to be responsible for the 2175 \AA\ extinction bump and/or the far ultraviolet (UV) extinction rise at $\lambda^{-1} > 5.9\ {\mu}m^{-1}$. Therefore, one might expect some relation between the UV extinction and DIBs. Such a relationship, if established, could put important constraints on the carrier of DIBs. Over the past four decades, whether DIBs are related to the shape of the UV extinction curves has been extensively investigated. However, the results are often inconsistent, partly due to the inconsistencies in characterizing the UV extinction. Here we re-examine the connection between the UV extinction curve and DIBs. We compile the extinction curves and the equivalent widths of 40 DIBs along 97 slightlines. We decompose the extinction curve into three Drude-like functions composed of the visible/near-infrared component, the 2175 \AA\ bump, and the far-UV extinction at $\lambda^{-1} > 5.9\ {\mu}m^{-1}$. We argue that the wavelength-integrated far-UV extinction derived from this decomposition technique best measures the strength of the far-UV extinction. No correlation is found between the far-UV extinction and most (~90\%) of the DIBs. We have also shown that the color excess E(1300-1700), the extinction difference at 1300 \AA\ and 1700 \AA\ often used to measure the strength of the far-UV extinction, does not correlate with DIBs. Finally, we confirm the earlier findings of no correlation between the 2175 \AA\ bump and DIBs or between the 2175 \AA\ bump and the far-UV extinction.

We refine a classic relation between the hydrogen column density (NH) and optical extinction (Av) by employing 39 Galactic Supernova Remnants (SNRs) with X-rays, optical and/or infra-red data available. We find NH = (1.69+/-0.07)*10^21 Av cm^(-2) mag^(-1) . Applying this relation to three Galactic SNRs with good historical records allows us to further constrain either their progenitor's distances or magnitudes, which is independent access to their distances.

We perform a survey of the Cr, Mn and Fe-K emission lines in young supernova remnants (SNRs) with the Japanese X-ray astronomy satellite {\sl Suzaku}. The Cr and/or Mn emission lines are detected in 3C\,397 and 0519-69.0 for the first time. We also confirm the detection of these lines in Kepler, W49B, N103B and Cas A. We derive the line parameters (i.e., the line centroid energy, flux and equivalent width [EW]) for these six sources and perform a correlation analysis for the line center energies of Cr, Mn and Fe. Also included in the correlation analysis are Tycho and G344.7-0.1 for which the Cr, Mn and Fe-K line parameters were available in the literature through {\sl Suzaku} observations. We find that the line center energies of Cr correlates very well with that of Fe and that of Mn. This confirms our previous findings that the Cr, Mn and Fe are spatially co-located, share a similar ionization state, and have a common origin in the supernova nucleo-synthesis. We find that the ratio of the EW of the Cr emission line to that of Fe ($\EWR\equiv \EWCr/\EWFe$) provides useful constraints on the SNR progenitors and on the SN explosion mechanisms: for SNRs with $\EWR > 2%$, a type Ia origin is favored (e.g., N103B, G344.7-0.1, 3C\,397 and 0519-69.0); for SNRs with $\EWR < 2%$, they could be of either core-collpase origin or carbon-deflagration Ia origin.

The diffuse interstellar bands (DIBs) are ubiquitous absorption spectral features arising from the tenuous material in the space between stars -- the interstellar medium (ISM). Since their first detection nearly nine decades ago, over 400 DIBs have been observed in the visible and near-infrared wavelength range in both the Milky Way and external galaxies, both nearby and distant. However, the identity of the species responsible for these bands remains as one of the most enigmatic mysteries in astrophysics. An equally mysterious interstellar spectral signature is the 2175 Angstrom extinction bump, the strongest absorption feature observed in the ISM. Its carrier also remains unclear since its first detection 46 years ago. Polycyclic aromatic hydrocarbon (PAH) molecules have long been proposed as a candidate for DIBs as their electronic transitions occur in the wavelength range where DIBs are often found. In recent years, the 2175 Angstrom extinction bump is also often attributed to the \pi--\pi* transition in PAHs. If PAHs are indeed responsible for both the 2175 Angstrom extinction feature and DIBs, their strengths may correlate. We perform an extensive literature search for lines of sight for which both the 2175 Angstrom extinction feature and DIBs have been measured. Unfortunately, we found no correlation between the strength of the 2175 Angstrom feature and the equivalent widths of the strongest DIBs. A possible explanation might be that DIBs are produced by small free gas-phase PAH molecules and ions, while the 2175 Angstrom bump is mainly from large PAHs or PAH clusters in condensed phase so that there is no tight correlation between DIBs and the 2175 Angstrom bump.

Using a deprojection technique, we study the X-ray properties of the galaxy cluster Abell 1835 observed with Chandra, including temperature, abundance, electron density, gas mass fraction, and total mass. A comparison with the results without deprojection shows that the properties do not change much. When we compare the results with those of XMM-Newton, the difference between the temperature profiles derived from Chandra and XMM-Newton data still exists, even if the point-spread function effect of XMM-Newton is corrected. To investigate the reasons for the difference, we used the spectra to cross-calibrate the temperatures. They show that the Chandra spectra can be fitted well with XMM-Newton temperatures. Furthermore, we derive the electron density profile from Chandra data with XMM-Newton temperatures and calculate the projected mass, which is consistent with the XMM-Newton mass and a little lower than the weak lensing mass at r_500. Thus, it seems that the temperature derived from XMM-Newton may be more reliable.

We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.

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

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

We consider the impact of thermal conduction in clusters of galaxies on the (unmagnetized) interface between a cold gaseous cloud and a hotter gas flowing over the cloud (the so-called cold front). We argue that near the stagnation point of the flow conduction creates a spatially extended layer of constant thickness $\Delta$, where $\Delta$ is of order $\sim\sqrt{kR/U}$, and $R$ is the curvature radius of the cloud, $U$ is the velocity of the flow at infinity, and $k$ is the conductivity of the gas. For typical parameters of the observed fronts, one finds $\Delta \ll R$. The formation time of such a layer is $\sim R/U$. Once the layer is formed, its thickness only slowly varies with time and the quasi-steady layer may persist for many characteristic time scales. Based on these simple arguments one can use the observed width of the cold fronts in galaxy clusters to constrain the effective thermal conductivity of the intra-cluster medium.

We present the results from a detailed deprojection analysis of Abell 1835 as observed by {\it XMM-Newton}. If we fit the spectra with an isothermal plasma model, the deprojected temperature profile is flat in the outer region around 7.6 keV and decreases to $\sim$ 5.6 keV in the center, which may be connected with the gas cooling. In the central part, a two-component thermal plasma model can fit the spectrum significantly better. Moreover, the cool component (T $\sim$ 1.8 keV) has a much lower metal abundance than the hot component (T $\sim$ 8 keV), which may be due to the longer cooling time for the cool gas with lower abundance. In addition, it was found that without a main isothermal component, the standard cooling flow model cannot fit the spectrum satisfactorily. From the isothermal model fitting results we also derived the electron density $n_e$, and fitted its radial distribution with a double-$\beta$ model. The $n_e$ profile inferred with the double-$\beta$ model and the deprojected X-ray gas temperature profile were then combined to derive the total mass and the total projected mass of the cluster. The projected mass is lower than that derived from the weak lensing method. However, assuming that the cluster extends to a larger radius $\sim15'$ as found by Clowe & Schneider (2002), the two results are consistent within the error bars. Furthermore, we calculated the projected mass within the radius of $\sim$ 153 kpc implied by the presence of a gravitational lensing arc, which is about half of the mass determined from the optical lensing.