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 statistically examine the gamma-ray burst (GRB) photon indices obtained by the Fermi-GBM and Fermi-LAT observations and compare the LAT GRB photon indices to the GBM GRB photon indices. We apply the jitter radiation to explain the GRB spectral diversities in the high-energy bands. In our model, the jitter radiative spectral index is determined by the spectral index of the turbulence. We classify GRBs into three classes depending on the shape of the GRB high-energy spectrum when we compare the GBM and LAT detections: the GRB spectrum is concave (GRBs turn out to be softer and are labeled as S-GRBs), the GRB spectrum is convex (GRBs turn out to be harder and are labeled as H-GRBs), and the GRBs have no strong spectral changes (labeled as N-GRBs). A universal Kolmogorov index 7/3 in the turbulent cascade is consistent with the photon index of the N-GRBs. The S-GRB spectra can be explained by the turbulent cascade due to the kinetic magnetic reconnection with the spectral index range of the turbulence from 8/3 to 3.0. The H-GRB spectra originate from the inverse turbulent cascade with the spectral index range of the turbulence from 2.0 to 3.5 that occurred during the large lengthscale magnetic reconnection. Thus, the GRB radiative spectra are diversified because the turbulent cascade modifies the turbulent energy spectrum. More observational samples are expected in the future to further identify our suggestions.

We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

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.

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.

We report observations of a transient source \fermi\ from radio to \grs. \fermi\ was discovered by the {\it Fermi-LAT} in May 2017. Follow-up {\it Swift-XRT} observations revealed three flaring episodes through March 2018, and the peak X-ray flux is about $10^3$ higher than the {\it ROSAT all-sky survey (RASS)} flux upper limit. Optical spectral measurements taken by the {\it Magellan 6.5-m telescope} and the {\it Lick-Shane telescope} both show a largely featureless spectrum, strengthening the BL Lac interpretation first proposed by \citet{Bruni18}. The optical and mid-infrared (MIR) emission goes to a higher state in 2018, when the flux in high energies goes down to a lower level. Our {\it RATAN-600m} measurements at 4.8~GHz and 8.2~GHz do not indicate any significant radio flux variation over the monitoring seasons in 2017 and 2018, nor deviate from the archival {\it NVSS} flux level. During GeV flaring times, the spectrum is very hard ($\Gamma_\gamma\sim$1.7) in the GeV band and at times also very hard (($\Gamma_{\rm X}\lesssim2$) in the X-rays, similar to a high-synchrotron-peak (or even an extreme) BL Lac object, making \fermi\ a good target for ground-based {\it Cherenkov telescopes}.

Context. Obscuration events caused by outflowing clumps or streams of high column density, low ionisation gas, heavily absorbing the X-ray continuum, have been witnessed in a number of Seyfert galaxies. Aims. We report on the X-ray spectral-timing analysis of the December 2016 obscuration event in NGC 3783, aimed at probing variability of the X-ray obscurer on the shortest possible timescales. The main goals of this study are to obtain independent constraints on the density, and ultimately on the distance of the obscuring gas, as well as to characterise the impact of variable obscuration on the observed X-ray spectral-timing characteristics of Seyfert galaxies. Methods. We carried out a comparative analysis of NGC 3783 during unobscured (using archival 2000-2001 XMM-Newton data) and obscured states (using XMM-Newton and NuSTAR data from the 2016 observational campaign). The timescales analysed range between ten hours and about one hour. This study was then generalized to discuss the signatures of variable obscuration in the X-ray spectral-timing characteristics of Seyfert galaxies as a function of the physical properties of the obscuring gas. Results. The X-ray obscurer in NGC 3783 is found to vary on timescales between about one hour to ten hours. This variability is incoherent with the variations of the X-ray continuum. A fast response (on timescales shorter than about 1.5 ks) of the ionisation state of the obscuring gas to the short timescale variability of the primary X-ray continuum provides a satisfactory interpretation of all the observed X-ray spectral-timing properties. This study enabled us to put independent constraints on the density and location of the obscuring gas. We found the gas to have a density of $n_{e}> 7.1 \times 10^7 \rm{cm^{-3}}$, consistent with being part of the broad line region.

Emission and absorption features from C-like ions serve as temperature and density diagnostics of astrophysical plasmas. $R$-matrix electron-impact excitation data sets for C-like ions in the literature merely cover a few ions, and often only for the ground configuration. Our goal is to obtain level-resolved effective collision strength over a wide temperature range for C-like ions from \ion{N}{II} to \ion{Kr}{XXXI} (i.e., N$^{+}$ to Kr$^{30+}$) with a systematic set of $R$-matrix calculations. We also aim to assess their accuracy. For each ion, we included a total of 590 fine-structure levels in both the configuration interaction target and close-coupling collision expansion. These levels arise from 24 configurations $2l^3 nl^{\prime}$ with $n=2-4$, $l=0-1$, and $l^{\prime}=0-3$ plus the three configurations $2s^22p5l$ with $l=0-2$. The AUTOSTRUCTURE code was used to calculate the target structure. Additionally, the $R$-matrix intermediate coupling frame transformation method was used to calculate the collision strengths. We compare the present results of selected ions with archival databases and results in the literature. The comparison covers energy levels, transition rates, and effective collision strengths. We illustrate the impact of using the present results on an \ion{Ar}{xiii} density diagnostic for the solar corona. The electron-impact excitation data is archived according to the Atomic Data and Analysis Structure (ADAS) data class adf04 and will be available in OPEN-ADAS. The data will be incorporated into spectral codes, such as CHIANTI and SPEX, for plasma diagnostics.

We present the SPEX software package for modeling and fitting X-ray spectra. Our group has developed spectral models, atomic data and code for X-ray applications since the 1970's. Since the 1990's these are further developed in the public SPEX package. In the last decades, X-ray spectroscopy has been revolutionized by the high-resolution grating spectrometers aboard XMM-Newton and Chandra. Currently, new high-resolution detectors aboard the Hitomi mission, and future missions, like XRISM and Athena, will provide another major step forward in spectral resolution. This poses challenges for us to increase the atomic database substantially, while keeping model calculation times short. In this paper, we summarize our efforts to prepare the SPEX package for the next generation of X-ray observatories.

Accurate chemical abundance measurements of X-ray emitting atmospheres pervading massive galaxies, galaxy groups, and clusters provide essential information on the star formation and chemical enrichment histories of these large scale structures. Although the collisionally ionised nature of the intracluster medium (ICM) makes these abundance measurements relatively easy to derive, underlying spectral models can rely on different atomic codes, which brings additional uncertainties on the inferred abundances. Here, we provide a simple, yet comprehensive comparison between the codes SPEXACT v3.0.5 (cie model) and AtomDB v3.0.9 (vapec model) in the case of moderate, CCD-like resolution spectroscopy. We show that, in cool plasmas ($kT \lesssim 2$ keV), systematic differences up to $\sim$20% for the Fe abundance and $\sim$45% for the O/Fe, Mg/Fe, Si/Fe, and S/Fe ratios may still occur. Importantly, these discrepancies are also found to be instrument-dependent, at least for the absolute Fe abundance. Future improvements in these two codes will be necessary to better address questions on the ICM enrichment.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

Gamma-ray bursts (GRBs) can be divided into three subclasses: X-ray flash (XRF), X-ray rich (XRR), and classical GRB (C-GRB). An X-ray flare is the rebrightening emission shown in the early X-ray afterglow of some GRBs. In this paper, we comprehensively examine the X-ray flare properties among XRF, XRR, and C-GRB subclasses. We utilize the XRF, XRR, and C-GRB subclass samples obtained from the Swift-BAT3 catalog, and the X-ray flare observational properties are collected from Falcone et al., Chincarini et al., and Yi et al. We find that XRFs and XRRs have more bright X-ray flares than C-GRBs. The ratio of the X-ray flare fluence to the prompt emission fluence has different distributions between XRF and C-GRB subclasses. The linear correlation between the duration and the peak time of the X-ray flares is also different between XRF and C-GRB subclasses. We are inclined to identify the GRBs with the bright X-ray flares as XRFs or XRRs. We discuss some issues that are related to the XRF/XRR/C-GRB classification. We also caution the selection effects and the instrument bias in our investigation. Large samples are required in the future to further confirm our results.

Aims. We aim to investigate and characterise the photoionised X-ray emission line regions within NGC 7469. Methods. We apply the photoionisation model, PION, within the spectral fitting code SPEX to analyse the 640 ks RGS spectrum of NGC 7469 gathered during an XMM-Newton observing campaign in 2015. Results. We find the emission line region in NGC 7469 to be multiphased, consisting of two narrow components with ionisation parameters of $\log \xi = 0.4$ and 1.6. A third, broad emission component, with a broadening velocity of $v_b \sim 1400$ km s\textsuperscript{-1} and an outflow velocity of $v_{out} \sim -4500$ km s\textsuperscript{-1}, is required to fit the residuals in the O VII triplet, at around 22 \AA. Assuming a volume filling factor of 0.1, the lower distance limits of the narrow emission line region components are estimated for the first time at 2.6 and 2.5 pc from the central black hole, whereas the broad component has an estimated lower bound distance between 0.004 to 0.03 pc, depending on the assumed plasma parameters. The collisionally ionised plasma from the star burst region in NGC 7469 has a plasma temperature of 0.32 keV and outflow velocity of $-280$ km s\textsuperscript{-1}, consistent with previous results in this campaign. In addition, we model the photoionised plasma of the warm absorber (WA) in NGC 7469, and find that it consists of three photoionised phases, with different values of $\xi$, $N_H$ and $v_{out}$. The upper bound distances of these WA components are 1.9, 0.3 and 0.6 pc, respectively, consistent with archival results. Conclusions. The environment of NGC 7469 is a complex mix of plasma winds absorbing and emitting X-rays. We find the picture painted by our results can be attributed to line emitting plasma located at distances ranging from near the black hole to the torus and beyond the ionised outflows.

Event rate of long Gamma-Ray Bursts (GRBs) is expected to be an useful tracer of the cosmic star-formation history. For this purpose, it is necessary to understand what kind of star formations/galaxies are traced by GRBs. Here we report rest-frame far-infrared (FIR) continuum detections of GRB 070521 and 080207 host galaxies at $z\sim2$ with ALMA band 8 and 9. The FIR photometries provide the reliable star-formation rates (SFRs), because FIR emission is free from dust extinction and possible radio contamination from long-lived afterglows of GRBs. The spectral energy distribution fittings indicate 49.85$^{+72.33}_{-2.86}$ and 123.4$^{+25.19}_{-21.78}$ $M_{\odot}$ yr$^{-1}$ for the 070521 and 080207 hosts, respectively. The derived SFRs place them on the \lq \lq main sequence\rq \rq\ of normal star-forming galaxies at $z\sim2$. The derived SFRs are significantly lower than that of radio observations. It is inferred that the observed radio fluxes in a previous study are contaminated by the afterglows. ALMA marginally detected [C~{\sc ii}]\,158\,$\mu$m emission line from the GRB 080207 host galaxy with S/N $\sim$ 4. This is the first detection of [C~{\sc ii}]\,158\,$\mu$m of a GRB host at $z>2$, and the second detection among known GRBs. The luminosity ratio of [C~{\sc~ii}]\,158$\mu$m to FIR is 7.5$\times 10^{-4}$, which is one of the smallest values among galaxies at $z\sim 1-2$ with the same FIR luminosity. The \lq \lq [C~{\sc ii}] deficit\rq \rq\ could be a new physical property to characterise GRB hosts at $z\sim1-2$. Possible parameters controlling the deficit include the metallicity, initial mass function, and gas density.

The high-resolution X-ray spectrum of the Perseus galaxy cluster observed with the Hitomi satellite challenges astrophysical collisional ionized plasma models that are widely used in the community. Although Hitomi spun out of control, several Hitomi-level missions have been proposed and some funded. The spectrometers aboard these future missions have a broader energy range and/or a higher spectral resolution to achieve different scientific goals. Accurate plasma models and atomic data are crucial for plasma diagnostics of high-quality spectra. Here, we present a few cases where improvement of plasma diagnostics will be decisive to better understand celestial bodies and their physical processes at play. We focus on collisional ionized and photoionized astrophysical plasmas in the context of developments of plasma models, as well as current and future generations of spectrometers.

Clusters of galaxies are the largest known gravitationally-bound structures in the Universe. When clusters collide, they create merger shocks on cosmological scales, which transform most of the kinetic energy carried by the cluster gaseous halos into heat. Observations of merger shocks provide key information of the merger dynamics, and enable insights into the formation and thermal history of the large-scale structures. Nearly all of the merger shocks are found in systems where the clusters have already collided, knowledge of shocks in the pre-merger phase is a crucial missing ingredient. Here we report on the discovery of a unique shock in a cluster pair 1E 2216 and 1E 2215. The two clusters are observed at an early phase of major merger. Contrary to all the known merger shocks observed ubiquitously on merger axes, the new shock propagates outward along the equatorial plane of the merger. This discovery uncovers an important epoch in the formation of massive clusters, when the rapid approach of the cluster pair leads to strong compression of gas along the merger axis. Current theoretical models predict that the bulk of the shock energy might be dissipated outside the clusters, and eventually turn into heat of the pristine gas in the circum-cluster space.

It has been reported that some X-ray spectra of gamma-ray burst (GRB) afterglows cannot be fitted by a simple power law. A blackbody component is added to precisely fit the thermal feature in these spectra. Alternatively, we propose that bremsstrahlung radiation can also be one possible mechanism to explain the thermal component of the GRB X-ray afterglow. In particular, we examine the X-ray afterglow of the ultra-long GRB 130925A in this paper. By our calculation, we find that the X-ray thermal component observed by both Swift-XRT and NuSTAR can be well explained by the bremsstrahlung radiation. Our results indicate that the GRBs with the bremsstrahlung emission in the X-ray afterglow could be born in a metal-rich and dusty environment.

The Hitomi results on the Perseus cluster lead to improvements in our knowledge of atomic physics which are crucial for the precise diagnostic of hot astrophysical plasma observed with high-resolution X-ray spectrometers. However, modeling uncertainties remain, both within but especially beyond Hitomi's spectral window. A major challenge in spectral modeling is the Fe-L spectrum, which is basically a complex assembly of n>2 to n=2 transitions of Fe ions in different ionization states, affected by a range of atomic processes such as collisional excitation, resonant excitation, radiative recombination, dielectronic recombination, and innershell ionization. In this paper we perform a large-scale theoretical calculation on each of the processes with the flexible atomic code (FAC), focusing on ions of Fe XVII to Fe XXIV that form the main body of the Fe-L complex. The new data are found to be consistent within 20% with the recent individual R-matrix calculations for the main Fe-L lines. By further testing the new FAC calculations with the high-quality RGS data from 15 elliptical galaxies and galaxy clusters, we note that the new model gives systematically better fits than the current SPEX v3.04 code, and the mean Fe abundance decreases by 12%, while the O/Fe ratio increases by 16% compared with the results from the current code. Comparing the FAC fit results to those with the R-matrix calculations, we find a temperature-dependent discrepancy of up to ~10% on the Fe abundance between the two theoretical models. Further dedicated tests with both observed spectra and targeted laboratory measurements are needed to resolve the discrepancies, and ultimately, to get the atomic data ready for the next high-resolution X-ray spectroscopy mission.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

High-resolution X-ray spectroscopy has advanced our understanding of the hot Universe by revealing physical properties like kinematics, temperature, and abundances of the astrophysical plasmas. Despite the technical and scientific achievements, the lack of scientific products at a level higher than count spectra is hampering full scientific exploitation of high-quality data. This paper introduces the Catalogue of Ionized Emission Lines Observed by the Reflection Grating Spectrometer (CIELO-RGS) onboard the XMM-Newton space observatory. The CIELO-RGS catalogue aims to facilitate the exploitation of emission features in the public RGS spectra archive, in particular, to perform the correlation between X-ray spectral diagnostics parameters with measurements at other wavelengths. This paper focuses on the methodology of catalogue generation, describing the automated line detection algorithm. A moderate sample (~2400 observations) of high-quality RGS spectra available at XMM-Newton Science Archive is used as our starting point. A list of potential emission lines is selected based on a multi-scale peak detection algorithm in a uniform and automated way without prior assumption on the underlying astrophysical model. The candidate line list is validated via spectral fitting with simple continuum and line profile models. We also compare the catalogue content with published literature results on a small number of exemplary sources. We generate a catalogue of emission lines ~12000 detected in ~1600 observations toward stars, X-ray binaries, supernovae remnants, active galactic nuclei, and groups and clusters of galaxies. For each line, we report the observed wavelength, broadening, energy and photon flux, and equivalent width, etc.

The stability of the flow under the magnetic force is one of the classical problems in fluid mechanics. In this paper, the flow in a rectangular duct with different Hartmann (Ha) number is simulated. The finite volume method and the SIMPLE algorithm are used to solve a system of equations and the energy gradient theory is then used to study the (associated) stability of magnetohydrodynamics (MHD). The flow stability of MHD flow for different Hartmann (Ha) number, from Ha=1 to 40, at the fixed Reynolds number, Re=190 are investigated. The simulation is validated firstly against the simulation in literature. The results show that, with the increasing Ha number, the centerline velocity of the rectangular duct with MHD flow decreases and the absolute value of the gradient of total mechanical energy along the streamwise direction increases. The maximum of K appears near the wall in both coordinate axis of the duct. According to the energy gradient theory, this position of the maximum of K would initiate flow instability (if any) than the other positions. The higher the Hartmann number is, the smaller the K value becomes, which means that the fluid becomes more stable in the presence of higher magnetic force. As the Hartmann number increases, the K value in the parallel layer decreases more significantly than in the Hartmann layer. The most dangerous position of instability tends to migrate towards wall of the duct as the Hartmann number increases. Thus, with the energy gradient theory, the stability or instability in the rectangular duct can be controlled by modulating the magnetic force.

Our Swift monitoring program triggered two joint XMM-Newton, NuSTAR and HST observations on 11 and 21 December 2016 targeting NGC 3783, as its soft X-ray continuum was heavily obscured. Consequently, emission features, including the O VII radiative recombination continuum, stand out above the diminished continuum. We focus on the photoionized emission features in the December 2016 RGS spectra and compare them to the time-averaged RGS spectrum obtained in 2000--2001 when the continuum was unobscured. A two-phase photoionized plasma is required to account for the narrow emission features. These narrow emission features are weakly varying between 2000--2001 and December 2016. We also find a statistically significant broad emission component in the time-averaged RGS spectrum in 2000--2001. This broad emission component is significantly weaker in December 2016, suggesting that the obscurer is farther away than the X-ray broad-line region. In addition, by analyzing the archival high-resolution X-ray spectra, we find that nine photoionized absorption components with different ionization parameters and kinematics are required for the warm absorber in X-rays.

The radio luminosity functions (RLFs) of active galactic nuclei (AGNs) are traditionally measured based on total emission, which doesn't reflect the current activity of the central black hole. The increasing interest in compact radio cores of AGNs requires determination of the RLF based on core emission (i.e., core RLF). In this work we have established a large sample (totaling 1207) of radio-loud AGNs, mainly consisting of radio galaxies (RGs) and steep-spectrum radio quasars (SSRQs). Based on the sample, we explore the relationship between core luminosity ($L_c$) and total luminosity ($L_t$) via a powerful statistical tool called "Copula". The conditional probability distribution $p(\log L_{c} \mid \log L_{t})$ is obtained. We derive the core RLF as a convolution of $p(\log L_{c} \mid \log L_{t})$ with the total RLF which was determined by previous work. We relate the separate RG and SSRQ core RLFs via a relativistic beaming model and find that SSRQs have an average Lorentz factor of $\gamma=9.84_{-2.50}^{+3.61}$, and that most are seen within $8^{\circ} \lesssim \theta \lesssim 45^{\circ}$ of the jet axis. Compared with the total RLF which is mainly contributed by extended emission, the core RLF shows a very weak luminosity-dependent evolution, with the number density peaking around $z\thicksim 0.8$ for all luminosities. Differences between core and total RLFs can be explained in a framework involving a combination of density and luminosity evolutions where the cores have significantly weaker luminosity evolution than the extended emission.

Following previous work, we further confirm that the cosmic evolution of steep-spectrum radio-loud AGNs (active galactic nuclei) can be reproduced by a simple combination of density evolution (DE) and luminosity evolution (LE). This mixture evolution scenario can naturally explain the luminosity-dependent evolution of radio-loud AGNs. Our models successfully fitted a large amount of data on radio luminosity functions (RLFs) of steep-spectrum sources and multi-frequency source counts. The modeling indicates that the DE slowly increase as $(1+z)^{0.3 \thicksim 1.3}$ out to $z \thicksim 0.8$, and then rapidly decreases as $(1+z)^{-6.8 \thicksim -5.7}$, while the LE rapidly increase as $(1+z)^{4.8}$ out to a higher redshift (at least $z>3.5$). We find a high-redshift decline (i.e. redshift cutoff) in the number density of steep-spectrum radio sources, but we cannot conclude whether such decline is sharp or shallow. We believe that whether a redshift cutoff occurs or not depends mainly on DE, while its steepness is decided by LE, which, however, cannot be well constrained due to the lack of high-redshift samples. Most intriguingly, according to our mixture evolution scenario, there appears to be no need for different evolution for the low- and high-power radio-loud AGNs. Both types of sources experience the same combined evolution of DE and LE.

To understand the nature of transient obscuring outflows in active galactic nuclei, we observed the Seyfert 1 galaxy NGC 3783 on two occasions in December 2016 triggered by Swift monitoring indicating strong soft X-ray absorption in November. We obtained ultraviolet spectra using COS on HST and optical spectra using FEROS on the MPG/ESO 2.2-m telescope that were simultaneous with X-ray spectra from XMM-Newton and NuSTAR. We find new components of broad, blue-shifted absorption associated with Ly$\alpha$, \ion{N}{v}, \ion{Si}{iv}, and \ion{C}{iv} in our COS spectra. The absorption extends from velocities near zero in the rest-frame of the host galaxy to $-6200$ $\rm km~s^{-1}$. These features appear for the first time in NGC 3783 at the same time as the heavy soft X-ray absorption seen in the XMM-Newton X-ray spectra. The X-ray absorption has a column density of $\sim 10^{23}~\rm cm^{-2}$, and it partially covers the X-ray continuum source. The X-ray absorption becomes more transparent in the second observation, as does the UV absorption. Combining the X-ray column densities with the UV spectral observations yields an ionization parameter for the obscuring gas of log $\xi =1.84^{+0.4}_{-0.2}$ $\rm erg~cm~s^{-1}$. Despite the high intensity of the UV continuum in NGC 3783, F(1470 \AA)=$8 \times 10^{-14}~\rm erg~cm^{-2}~s^{-1}~\AA^{-1}$, the well known narrow UV absorption lines are deeper than in earlier observations in unobscured states, and low ionization species such as \ion{C}{iii} appear, indicating that the narrow-line gas is more distant from the nucleus and is being shadowed by the gas producing the obscuration. Despite the high continuum flux levels in our observations of NGC 3783, moderate velocities in the UV broad line profiles have substantially diminished. We suggest that a collapse of the broad line region has led to the outburst and triggered the obscuring event.

We build a comprehensive sample to statistically describe the properties of X-ray flashes (XRFs) and X-ray riches (XRRs) from the latest third Swift Burst Alert Telescope (BAT3) catalog of Gamma-ray bursts (GRBs). We obtain 81 XRFs, 540 XRRs, and 394 classical GRBs (C-GRBs). We statistically explore the different properties of the $\gamma$-ray prompt emission, the X-ray emission, the X-ray lightcurve type, the association with supernovae (SNe), and the host galaxy properties for these sources. We confirm that most XRFs/XRRs are long GRBs with low values of peak energy $E^{obs}_{peak}$ and they are low-luminosity GRBs. XRFs, XRRs, and C-GRBs follow the same $E_{X,iso}$-$E_{\gamma,iso}$-$E_{peak,z}$ correlations. Compared to the classical GRBs, XRFs are favorable to have the association with SN explosions. We do not find any significant differences of redshift distribution and host galaxy properties among XRFs, XRRs, and C-GRBs. We also discuss some observational biases and selection effects that may affect on our statistical results. The GRB detectors with wide energy range and low energy threshold are expected for the XRF/XRR research in the future.

Chemical abundances in the X-ray halos (also known as the intracluster medium, ICM) of clusters and groups of galaxies can be measured via prominent emission line features in their X-ray spectra. Elemental abundances are footprints of time-integrated yields of various stellar populations that have left their specific abundance patterns prior to and during the cluster and group evolution. We aim to constrain nitrogen abundances in the CHEmical Evolution RGS Sample (CHEERS), which contains 44 nearby groups and clusters of galaxies, in order to have a better understanding of their chemical enrichment. We examine the high-resolution spectra of the CHEERS sample carefully and take into account various systematic effects in the spectral modelling. We compare the observed abundance ratios with those in the Galactic stellar populations, as well as predictions from stellar yields (low- and intermediate-mass stars, massive stars and degenerate stars). The nitrogen abundance can only be well constrained ($\gtrsim3\sigma$) in one cluster of galaxies and seven groups of galaxies. The [O/Fe] -- [Fe/H] relation of the ICM is comparable to that for the Galaxy, while both [N/Fe] and [N/O] ratios of the ICM are higher than in the Galaxy. Future studies on nitrogen radial distributions are required to tell whether the obtained higher [N/Fe] and [N/O] ratios are biased due to the small extraction region ($r/r_{500}\lesssim0.05$) that we adopt here. Since abundances of odd-$Z$ elements are more sensitive to the initial metallicity of stellar populations, accurate abundance measurements of N, Na and Al are required to better constrain the chemical enrichment in the X-ray halos of clusters and groups of galaxies.

Theoretically long gamma-ray bursts (GRBs) are expected to happen in low-metallicity environments, because in a single massive star scenario, low iron abundance prevents loss of angular momentum through stellar wind, resulting in ultra-relativistic jets and the burst. In this sense, not just a simple metallicity measurement but also low iron abundance ([Fe/H]<-1.0) is essentially important. Observationally, however, oxygen abundance has been measured more often due to stronger emission. In terms of oxygen abundance, some GRBs have been reported to be hosted by high-metallicity star-forming galaxies, in tension with theoretical predictions. Here we compare iron and oxygen abundances for the first time for GRB host galaxies (GRB 980425 and 080517) based on the emission-line diagnostics. The estimated total iron abundances, including iron in both gas and dust, are well below the solar value. The total iron abundances can be explained by the typical value of theoretical predictions ([Fe/H]<-1.0), despite high oxygen abundance in one of them. According to our iron abundance measurements, the single massive star scenario still survives even if the oxygen abundance of the host is very high, such as the solar value. Relying only on oxygen abundance could mislead us on the origin of the GRBs. The measured oxygen-to-iron ratios, [O/Fe], can be comparable to the highest values among the iron-measured galaxies in the Sloan Digital Sky Survey. Possible theoretical explanations of such high [O/Fe] include the young age of the hosts, top-heavy initial mass function, and fallback mechanism of the iron element in supernova explosions.

We present optical spectroscopic and photometric observations of the nearby type Ic supernova (SN Ic) SN 2014L. This SN was discovered by the Tsinghua-NAOC Transient Survey (TNTS) in the nearby type-Sc spiral galaxy M99 (NGC 4254). Fitting to the early-time light curve indicates that SN 2014L was detected at only a few hours after the shock breakout, and it reached a peak brightness of $M_{\rm V} = -17.73 \pm 0.28$ mag ($L= [2.06\pm 0.50] \times 10^{42}$ erg s$^{-1}$) approximately 13 days later. SN 2014L shows a close resemblance to SN 2007gr in the photometric evolution, while it shows stronger absorption features of intermediate-mass elements (especially \CaII) in the early-time spectra. Based on simple modeling of the observed light curves, we derived the mass of synthesized $^{56}$Ni as $M_{\rm Ni} = 0.075 \pm 0.025 M_{\odot}$, and the mass and total energy of the ejecta as $M_{\rm ej} = 1.00 \pm 0.20 M_\odot$ and $E_{\rm ej}=1.45\pm0.25$ foe, respectively. Given these typical explosion parameters, the early detection, and the extensive observations, we suggest that SN 2014L could be a template sample for the investigation of SNe Ic.

The hot intracluster medium (ICM) pervading galaxy clusters and groups is rich in metals, which were synthesised by billions of supernovae and have accumulated in cluster gravitational wells for several Gyrs. Since the products of both Type Ia and core-collapse supernovae - expected to explode over different time scales - are found in the ICM, constraining accurately the chemical composition these hot atmospheres can provide invaluable information on the history of the enrichment of large-scale structures. Recently, Hitomi observations reported solar abundance ratios in the core of the Perseus cluster, in tension with previous XMM-Newton measurements obtained for 44 cool-core clusters, groups, and massive ellipticals (the CHEERS sample). In this work, we revisit the CHEERS results by using an updated version of the spectral code used to fit the data (SPEXACT v3), the same as was used to obtain the Hitomi measurements. Despite limitations in the spectral resolution, the average Cr/Fe and Ni/Fe ratios are now found to be remarkably consistent with unity and in excellent agreement with the Hitomi results. Our updated measurements suggest that the solar composition of the ICM of Perseus is a very common feature in nearby cool-core systems.

X-ray measurements find systematically lower Fe abundances in the X-ray emitting haloes pervading groups ($kT\lesssim1.7$ keV) than in clusters of galaxies. These results have been difficult to reconcile with theoretical predictions. However, models using incomplete atomic data or the assumption of isothermal plasmas may have biased the best fit Fe abundance in groups and giant elliptical galaxies low. In this work, we take advantage of a major update of the atomic code in the spectral fitting package SPEX to re-evaluate the Fe abundance in 43 clusters, groups, and elliptical galaxies (the CHEERS sample) in a self-consistent analysis and within a common radius of 0.1$r_{500}$. For the first time, we report a remarkably similar average Fe enrichment in all these systems. Unlike previous results, this strongly suggests that metals are synthesised and transported in these haloes with the same average efficiency across two orders of magnitude in total mass. We show that the previous metallicity measurements in low temperature systems were biased low due to incomplete atomic data in the spectral fitting codes. The reasons for such a code-related Fe bias, also implying previously unconsidered biases in the emission measure and temperature structure, are discussed.

Relativistic thermal electrons moving in a large-scale magnetic field can produce synchrotron radiation. Linear synchrotron polarization can also be produced by the relativistic thermal electrons. In this paper, we utilize a hybrid thermal-nonthermal electron energy distribution to calculate circular synchrotron polarization. We further compute the radiative transfer of the synchrotron polarization in the optical and radio bands when we consider the contribution of the thermal electrons. We attempt to apply the polarization results to some astrophysical objects, such as kilonova like AT 2017gfo/GW170817, the fast radio burst (FRB), the Gamma-ray burst (GRB) afterglow, and the supernova remnant (SNR). The large optical depth of radiative transfer effects the small polarization degrees of these populations when the media surrounding the synchrotron sources take heavy absorption to the polarized photons. We need a strong magnetic field in our model to reproduce the linear and circular polarization properties that were observed in FRB 140514. This indicates that FRBs have a neutron star origin.

The occurrence rates of bright X-ray flares in z<1 gamma-ray bursts (GRBs) with or without observed supernovae (SNe) association were compared. Our Sample I: the z<1 long GRBs (LGRBs) with SNe association (SN-GRBs) and with early Swift/X-Ray Telescope (XRT) observations, consists of 18 GRBs, among which only two GRBs have bright X-ray flares. Our Sample II: for comparison, all the z<1 LGRBs without observed SNe association and with early Swift/XRT observations, consists of 45 GRBs, among which 16 GRBs present bright X-ray flares. Thus, the study indicates a lower occurrence rate of bright X-ray flares in Sample I (11.1%) than in Sample II (35.6%). In addition, if dim X-ray fluctuations are included as flares, then 16.7% of Sample I and 55.6% of Sample II are found to have flares, again showing the discrepancy between these two samples. We examined the physical origin of these bright X-ray flares and found that most of them are probably related to the central engine reactivity. To understand the discrepancy, we propose that such a lower occurrence rate of flares in the SN-GRB sample may hint at an energy partition among the GRB, SNe, and X-ray flares under a saturated energy budget of massive star explosion.

Obscuration of the continuum emission from active galactic nuclei by streams of gas with relatively high velocity (> 1000 km/s) and column density (>3E25 per m2) has been seen in a few Seyfert galaxies. This obscuration has a transient nature. In December 2016 we have witnessed such an event in NGC 3783. The frequency and duration of these obscuration events is poorly known. Here we study archival data of NGC 3783 in order to constrain this duty cycle. We use archival Chandra/NuSTAR spectra taken in August 2016. We also study the hardness ratio of all Swift XRT spectra taken between 2008-2017. In August 2016, NGC 3783 also showed evidence for obscuration. While the column density of the obscuring material is ten times lower than in December 2016, the opacity is still sufficient to block a significant fraction of the ionising X-ray and EUV photons. From the Swift hardness ratio behaviour we find several other epochs with obscuration. Obscuration with columns >1E26 per m2 may take place in about half of the time. Also in archival X-ray data taken by ASCA in 1993 and 1996 we find evidence for obscuration. Obscuration of the ionising photons in NGC 3783 occurs more frequently than previously thought. This may not always have been recognised due to low spectral resolution observations, too limited spectral bandwidth or confusion with underlying continuum variations.

Short gamma-ray bursts may originate from the merger of double neutron stars (NS) or that of a black hole (BH) and an NS. We propose that the bright X-ray flare related to the central engine reactivity may hint a BH-NS merger, since such a merger can provide more fall-back materials and therefore a more massive accretion disk than the NS-NS merger. Based on the observed 49 short bursts with Swift/X-ray Telescope follow-up observations, we find that three bursts have bright X-ray flares, among which three flares from two bursts are probably related to the central engine reactivity. We argue that these two bursts may originate from the BH-NS merger rather than the NS-NS merger. Our suggested link between the central engine-powered bright X-ray flare and the BH-NS merger event can be checked by the future gravitational wave detections from advanced LIGO and Virgo.

We investigate the physical structure of the AGN wind in the Seyfert-1 galaxy NGC 7469 through high-resolution X-ray spectroscopy with Chandra HETGS and photoionisation modelling. Contemporaneous data from Chandra, HST, and Swift are used to model the optical-UV-X-ray continuum and determine the spectral energy distribution (SED) at two epochs, 13 years apart. For our investigation we use new observations taken in December 2015-January 2016, and historical ones taken in December 2002. We study the impact of a change in the SED shape, seen between the two epochs, on the photoionisation of the wind. The HETGS spectroscopy shows that the AGN wind in NGC 7469 consists of four ionisation components, with their outflow velocities ranging from -400 to -1800 km/s. From our modelling we find that the change in the ionising continuum shape between the two epochs results in some variation in the ionisation state of the wind components. However, for the main ions detected in X-rays, the sum of their column densities over all four components, remains in practice unchanged. For two of the four components, which are found to be thermally unstable in both epochs, we obtain 2 < r < 31 pc and 12 < r < 29 pc using the cooling and recombination timescales. For the other two thermally stable components, we obtain r < 31 pc and r < 80 pc from the recombination timescale. The results of our photoionisation modelling and thermal stability analysis suggest that the absorber components in NGC 7469 are consistent with being a thermally-driven wind from the AGN torus. Finally, from analysis of the zeroth-order ACIS/HETG data, we discover that the X-ray emission between 0.2-1 keV is spatially extended over 1.5-12". This diffuse soft X-ray emission is explained by coronal emission from the nuclear starburst ring in NGC 7469.

We conducted a multi-wavelength six-month campaign to observe the Seyfert galaxy NGC~7469, using the space-based observatories \textit{HST}, \textit{Swift}, \textit{XMM-Newton} and \textit{NuSTAR}. Here we report the results of the spectral analysis of the 7 simultaneous \textit{XMM-Newton} and \textit{NuSTAR} observations. The sources shows significant flux variability within each observation, but the average flux is less variable among the different pointings of our campaign. Our spectral analysis reveals a prominent narrow neutral \ion{Fe} K$\alpha$ emission line in all the spectra, with weaker contributions from Fe K$\beta$, neutral Ni K$\alpha$ and ionised iron. We find no evidence for variability or relativistic effects acting on the emission lines, which indicates that they originate from distant material. Analysing jointly \textit{XMM-Newton} and \textit{NuSTAR} data a constant photon index is found ($\Gamma$=$1.78\pm0.02$), together with a high energy cut-off $E_{\rm{cut}}=170^{+60}_{-40}$ keV. Adopting a self-consistent Comptonization model, these values correspond to an average coronal electron temperature of kT=$45^{+15}_{-12}$ keV and, assuming a spherical geometry, an optical depth $\tau=2.6\pm0.9$. The reflection component is consistent with being constant, with a reflection fraction in the range $R=0.3-0.6$. A prominent soft excess dominates the spectra below 4 keV. This is best fit with a second Comptonization component, arising from a \virg{warm corona} with an average $kT=0.67\pm0.03$ keV and a corresponding optical depth $\tau=9.2\pm0.2$.

Relativistic electrons accelerated by both the first-order and the second-order Fermi accelerations in some synchrotron sources have a hybrid shape of thermal and nonthermal energy distribution. This particle acceleration result is supported by some recent numerical simulations. We calculate the synchrotron polarization by applying this electron energy distribution. The polarization degrees in the cases of active galactic nucleus (AGN) jets and gamma-ray bursts (GRBs) are given as examples. The possible application for the polarization study of Sgr A* is also mentioned. We finally suggest high-energy polarization measurements for these synchrotron sources to test our results.

The X-ray narrow emission line region (NELR) of the archetypal Seyfert 1 galaxy NGC\,5548 has been interpreted as a single-phase photoionized plasma that is absorbed by some of the warm absorber components. This scenario requires those overlaying warm absorber components to have larger distance (to the central engine) than the X-ray NELR, which is not fully consistent with the distance estimates found in the literature. Therefore, we reanalyze the high-resolution spectra obtained in 2013--2014 with the Reflection Grating Spectrometer (RGS) aboard \textit{XMM}-Newton to provide an alternative interpretation of the X-ray narrow emission features. We find that the X-ray narrow emission features in NGC\,5548 can be described by a two-phase photoionized plasma with different ionization parameters ($\log \xi=1.3$ and $0.1$) and kinematics ($v_{\rm out}=-50$ and $-400~{\rm km~s^{-1}}$), and no further absorption by the warm absorber components. The X-ray and optical NELR might be the same multi-phase photoionized plasma. Both X-ray and optical NELR have comparable distances, asymmetric line profiles, and the underlying photoionized plasma is turbulent and compact in size. The X-ray NELR is not the counterpart of the UV/X-ray absorber outside the line of sight because their distances and kinematics are not consistent. In addition, X-ray broad emission features that we find in the spectrum can be accounted for by a third photoionized emission component. The RGS spectrum obtained in 2016 is analyzed as well, where the luminosity of most prominent emission lines (the \ion{O}{vii} forbidden line and \ion{O}{viii} Ly$\alpha$ line) are the same (at a 1 $\sigma$ confidence level) as in 2013--2014.

We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "kilonova" or "macronova". We report the detection of linear optical polarization P = (0.50 +/- 0.07)% at 1.46 days after detection of the GWs from GW170817, a double neutron star merger associated with an optical macronova counterpart and a short GRB. The optical emission from a macronova is expected to be characterized by a blue, rapidly decaying, component and a red, more slowly evolving, component due to material rich of heavy elements, the lanthanides. The polarization measurement was made when the macronova was still in its blue phase, during which there is an important contribution from a lanthanide-free outflow. The low degree of polarization is consistent with intrinsically unpolarized emission scattered by Galactic dust, suggesting a symmetric geometry of the emitting region and low inclination of the merger system. Stringent upper limits to the polarization degree from 2.45 - 9.48 days post-burst are consistent with the lanthanides-rich macronova interpretation.

We present the spectroscopic evolution of AT 2017gfo, the optical counterpart of the first binary neutron star (BNS) merger detected by LIGO and Virgo, GW170817. While models have long predicted that a BNS merger could produce a kilonova (KN), we have not been able to definitively test these models until now. From one day to four days after the merger, we took five spectra of AT 2017gfo before it faded away, which was possible because it was at a distance of only 39.5 Mpc in the galaxy NGC 4993. The spectra evolve from blue ($\sim6400$K) to red ($\sim3500$K) over the three days we observed. The spectra are relatively featureless --- some weak features exist in our latest spectrum, but they are likely due to the host galaxy. However, a simple blackbody is not sufficient to explain our data: another source of luminosity or opacity is necessary. Predictions from simulations of KNe qualitatively match the observed spectroscopic evolution after two days past the merger, but underpredict the blue flux in our earliest spectrum. From our best-fit models, we infer that AT 2017gfo had an ejecta mass of $0.03M_\odot$, high ejecta velocities of $0.3c$, and a low mass fraction $\sim10^{-4}$ of high-opacity lanthanides and actinides. One possible explanation for the early excess of blue flux is that the outer ejecta is lanthanide-poor, while the inner ejecta has a higher abundance of high-opacity material. With the discovery and follow-up of this unique transient, combining gravitational-wave and electromagnetic astronomy, we have arrived in the multi-messenger era.

We report on SALT low resolution optical spectroscopy and optical/IR photometry undertaken with other SAAO telescopes (MASTER-SAAO and IRSF) of the kilonova AT 2017gfo (aka SSS17a) in the galaxy NGC4993 during the first 10 days of discovery. This event has been identified as the first ever electromagnetic counterpart of a gravitational wave event, namely GW170817, which was detected by the LIGO and Virgo gravitational wave observatories. The event is likely due to a merger of two neutron stars, resulting in a kilonova explosion. SALT was the third telescope to obtain spectroscopy of AT 2017gfo and the first spectrum, 1.2 d after the merger, is quite blue and shows some broad features, but no identifiable spectral lines and becomes redder over time. We compare the spectral and photometric evolution with recent kilonova simulations and conclude that they are in qualitative agreement for post-merger wind models with proton: nucleon ratios of $Y_e$ = 0.25$-$0.30. The blue colour of the first spectrum is consistent with the lower opacity of the Lathanide-free r-process elements in the ejecta. Differences between the models and observations are likely due to the choice of system parameters combined with the absence of atomic data for more elements in the ejecta models.

Though theoretically expected, the charge exchange emission from galaxy clusters has not yet been confidently detected. Accumulating hints were reported recently, including a rather marginal detection with the Hitomi data of the Perseus cluster. As suggested in Gu et al. (2015), a detection of charge exchange line emission from galaxy clusters would not only impact the interpretation of the newly-discovered 3.5 keV line, but also open up a new research topic on the interaction between hot and cold matter in clusters. We aim to perform the most systematic search for the O VIII charge exchange line in cluster spectra using the RGS on board XMM. We introduce a sample of 21 clusters observed with the RGS. The dominating thermal plasma emission is modeled and subtracted with a two-temperature CIE component, and the residuals are stacked for the line search. The systematic uncertainties in the fits are quantified by refitting the spectra with a varying continuum and line broadening. By the residual stacking, we do find a hint of a line-like feature at 14.82 A, the characteristic wavelength expected for oxygen charge exchange. This feature has a marginal significance of 2.8 sigma, and the average equivalent width is 2.5E-4 keV. We further demonstrate that the putative feature can be hardly affected by the systematic errors from continuum modelling and instrumental effects, or the atomic uncertainties of the neighbouring thermal lines. Assuming a realistic temperature and abundance pattern, the physical model implied by the possible oxygen line agrees well with the theoretical model proposed previously to explain the reported 3.5 keV line. If the charge exchange source indeed exists, we would expect that the oxygen abundance is potentially overestimated by 8-22% in previous X-ray measurements which assumed pure thermal lines.

Ionized outflows in Active Galactic Nuclei (AGN) are thought to influence their nuclear and local galactic environment. However, the distance of the outflows with respect to the central engine is poorly constrained, which limits our understanding of their kinetic power as a cosmic feedback channel. Therefore, the impact of AGN outflows on their host galaxies is uncertain. However, when the density of the outflows is known, their distance can be immediately obtained from their modelled ionization parameter. With the new self-consistent PhotoIONization (PION) model in the SPEX code, we are able to calculate detailed level populations, including the ground and metastable levels. This enables us to determine under what physical conditions the metastable levels are significantly populated. We then identify characteristic lines from these metastable levels in the 1 -- 2000 {\AA} wavelength range. In the large density range of $n_H \in (10^6, 10^{20} m^{-3}$, the metastable levels 2s 2p $(^3P_{0-2})$ in Be-like ions can be significantly populated. For B-like ions, merely the first excited level 2s$^2$ 2p $(^2P_{3/2})$ can be used as a density probe. For C-like ions, the first two excited levels 2s$^2$ 2p$^2$ ($^3P_1$ and $^3P_2$) are better density probes than the next two excited levels 2s$^2$ 2p$^2$ ($^1S_0$ and $^1D_2$). Different ions in the same isoelectronic sequence cover not only a wide range of ionization parameter, but also a wide range of density. On the other hand, within the same isonuclear sequence, less ionized ions probe lower density and smaller ionization parameter. Finally, we re-analyzed the high-resolution grating spectra of NGC 5548 observed with Chandra in January 2002, using a set of PION components to account for the ionized outflow. We derive lower (or upper) limits of plasma density in five out of six PION components, based on the presence (or absence) of the metastable absorption lines.

We present multi-wavelength follow-up campaigns by the AstroSat-CZTI and GROWTH collaborations to search for an electromagnetic counterpart to the gravitational wave event GW170104. At the time of the GW170104 trigger, the AstroSat CZTI field-of-view covered 50.3\% of the sky localization. We do not detect any hard X-ray (>100 keV) signal at this time, and place an upper limit of $\approx 4.5 \times 10^{-7}~{\rm erg~cm}^{-2}{\rm~s}^{-1}$ for a 1\,s timescale. Separately, the ATLAS survey reported a rapidly fading optical source dubbed ATLAS17aeu in the error circle of GW170104. Our panchromatic investigation of ATLAS17aeu shows that it is the afterglow of an unrelated long, soft GRB~170105A, with only a fortuitous spatial coincidence with GW170104. We then discuss the properties of this transient in the context of standard long GRB afterglow models.

In 2016 we carried out a Swift monitoring program to track the X-ray hardness variability of eight type-I AGN over a year. The purpose of this monitoring was to find intense obscuration events in AGN, and thereby study them by triggering joint XMM-Newton, NuSTAR, and HST observations. We successfully accomplished this for NGC 3783 in December 2016. We found heavy X-ray absorption produced by an obscuring outflow in this AGN. As a result of this obscuration, interesting absorption features appear in the UV and X-ray spectra, which are not present in the previous epochs. Namely, the obscuration produces broad and blue-shifted UV absorption lines of Ly$\alpha$, C IV, and N V, together with a new high-ionisation component producing Fe XXV and Fe XXVI absorption lines. In soft X-rays, only narrow emission lines stand out above the diminished continuum as they are not absorbed by the obscurer. Our analysis shows that the obscurer partially covers the central source with a column density of few $10^{23}$ cm$^{-2}$, outflowing with a velocity of few thousand km s$^{-1}$. The obscuration in NGC 3783 is variable and lasts for about a month. Unlike the commonly-seen warm-absorber winds at pc-scale distances from the black hole, the eclipsing wind in NGC 3783 is located at about 10 light days. Our results suggest the obscuration is produced by an inhomogeneous and clumpy medium, consistent with clouds in the base of a radiatively-driven disk wind at the outer broad-line region of the AGN.

We present some recent highlights and prospects for the study of hot astrophysical plasmas. Hot plasmas can be studied primarily through their X-ray emission and absorption. Most astrophysical objects, from solar system objects to the largest scale structures of the Universe, contain hot gas. In general we can distinguish collisionally ionised gas and photoionised gas. We introduce several examples of both classes and show where the frontiers of this research in astrophysics can be found. We put this also in the context of the current and future generation of X-ray spectroscopy satellites. The data coming from these missions challenge the models that we have for the calculation of the X-ray spectra.

We carefully examine the depolarization feature of blazars in the optical and near-infrared bands using the sample of Mead et al. Magnetohydrodynamics (MHD) turbulence could be one possible reason for the depolarization of optical/infrared blazars when we apply the theoretical analysis of Lazarian and Pogosyan. We further identify in the sample that the depolarization results shown in most blazars roughly obey the form of the three-dimensional anisotropic Kolmogorov scaling. The effective Faraday rotation window length scale is not small enough to resolve the polarization correlation length scale in the blazar sample. The depolarization and the related turbulent features show diversities in different blazar sources. We suggest more simultaneous observations in both the optical/infrared and the high-energy bands for the study of the blazar polarization.

The hot intra-cluster medium (ICM) permeating galaxy clusters and groups is not pristine, as it is continuously enriched by metals synthesised in Type Ia (SNIa) and core-collapse (SNcc) supernovae since the major epoch of star formation (z ~ 2-3). The cluster/group enrichment history and the mechanisms responsible for releasing and mixing the metals can be probed via the radial distribution of SNIa and SNcc products within the ICM. In this paper, we use deep XMM-Newton/EPIC observations from a sample of 44 nearby cool-core galaxy clusters, groups, and ellipticals (CHEERS) to constrain the average radial O, Mg, Si, S, Ar, Ca, Fe, and Ni abundance profiles. The radial distributions of all these elements, averaged over a large sample for the first time, represent the best constrained profiles available currently. We find an overall decrease of the Fe abundance with radius out to ~$0.9 r_{500}$ and ~$0.6 r_{500}$ for clusters and groups, respectively, in good agreement with predictions from the most recent hydrodynamical simulations. The average radial profiles of all the other elements (X) are also centrally peaked and, when rescaled to their average central X/Fe ratios, follow well the Fe profile out to at least ~0.5$r_{500}$. Using two sets of SNIa and SNcc yield models reproducing well the X/Fe abundance pattern in the core, we find that, as predicted by recent simulations, the relative contribution of SNIa (SNcc) to the total ICM enrichment is consistent with being uniform at all radii, both for clusters and groups. In addition to implying that the central metal peak is balanced between SNIa and SNcc, our results suggest that the enriching SNIa and SNcc products must share the same origin, and that the delay between the bulk of the SNIa and SNcc explosions must be shorter than the timescale necessary to diffuse out the metals.

Linear and circular polarizations of gamma-ray bursts (GRBs) have been detected during recent years. We adopt a simplified model to investigate GRB polarization characteristics in this paper. A compressed two-dimensional turbulent slab containing stochastic magnetic fields is considered, and jitter radiation can produce the linear polarization under this special magnetic field topology. Turbulent Faraday rotation measure (RM) of this slab makes strong wavelength-dependent depolarization. The jitter photons can also scatter with those magnetic clumps inside the turbulent slab, and a nonzero variance of the Stokes parameter $V$ can be generated. Furthermore, the linearly and circularly polarized photons in the optical and radio bands may suffer heavy absorptions from the slab. Thus we consider the polarized jitter radiation transfer processes. Finally, we compare our model results with the optical detections of GRB 091018, GRB 121024A, and GRB 131030A. We suggest simultaneous observations of GRB multi-wavelength polarization in the future.

Charge exchange emission is known to provide a key diagnostic to the interface between hot and cold matter in many astrophysical environments. Most of the recent charge exchange studies focus on its emission in the X-ray band, but few on the UV part, although the latter can also provide a powerful probe of the charge exchange process. An atomic calculation, as well as an application to observed data, are presented to explore and describe the potential use of the UV data for the study of cosmic charge exchange. Using the newest charge exchange model in the SPEX code v3.03, we re-analyze an archival Hubble STIS data of the central region of NGC 1275. The NGC 1275 spectrum shows hints for three possible weak lines at about 1223.6~{\AA}, 1242.4~{\AA}, and 1244.0~{\AA}, each with a significance of about $2-3\sigma$. The putative features are best explained by charge exchange between highly ionized hydrogen, neon, and sulfur with neutral matter. The wavelengths of the charge exchange lines are found robustly with uncertainties $\leq 0.3$~{\AA}. The possible charge exchange emission shows a line-of-sight velocity offset of about $-3400$ km s$^{-1}$ with respect to the NGC 1275 nucleus, which resembles one of the Ly$\alpha$ absorbers reported in Baum et al. (2005). This indicates that the charge exchange lines might be emitted as the same position of the absorber, which could be ascribed to outflowing gas from the nucleus.

For photoionized plasmas, electron energy-loss rates due to radiative recombination (RR) are required for thermal equilib- rium calculations, which assume a local balance between the energy gain and loss. While many calculations of total and/or partial RR rates are available from literature, specific calculations of associated RR electron energy-loss rates are lacking. Here we focus on electron energy-loss rates due to radiative recombination of H-like to Ne-like ions for all the elements up to and including zinc (Z = 30), over a wide temperature range. We use the AUTOSTRUCTURE code to calculate the level-resolved photoionization cross section and modify the ADASRR code so that we can simultaneously obtain level-resolved RR rate coefficients and associated RR electron energy-loss rate coefficients. The total RR rates and electron energy-loss rates of H i and He i are compared with those found in literature. Furthermore, we utilize and parameterize the weighted electron energy-loss factors (dimensionless) to characterize total electron energy-loss rates due to RR. The RR electron energy-loss data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf48. The RR electron energy-loss data are also incorporated into the SPEX code for detailed modelling of photoionized plamsas.

We present a multi-wavelength study of the unidentified Fermi object, 3FGL J0212.1+5320. Within the 95% error ellipse, Chandra detects a bright X-ray source (i.e., F(0.5-7keV) = 1.4e-12 erg/cm^2/s), which has a low-mass optical counterpart (M < 0.4 Msun and T ~ 6000 K). A clear ellipsoidal modulation is shown in optical/infrared at 20.87 hours. The gamma-ray properties of 3FGL J0212.1+5320 are all consistent with that of a millisecond pulsar, suggesting that it is a gamma-ray redback millisecond pulsar binary with a low-mass companion filling > 64% of the Roche-lobe. If confirmed, it will be a redback binary with one of the longest orbital periods known. Spectroscopic data taken in 2015 from the Lijiang observatory show no evidence of strong emission lines, revealing that the accretion is currently inactive (the rotation-powered pulsar state). This is consistent with the low X-ray luminosities (Lx ~ 10^32 erg/s) and the possible X-ray modulation seen by Chandra and Swift. Considering that the X-ray luminosity and the high X-ray-to-gamma-ray flux ratio (8%) are both comparable to that of the two known gamma-ray transitional millisecond pulsars, we suspect that 3FGL J0212.1+5320 could be a potential target to search for future transition to the accretion active state.

We build a simple physical model to study the high-redshift active galactic Nucleus (AGN) evolution within the co-evolution framework of central black holes (BHs) and their host galaxies. The correlation between the circular velocity of a dark halo $V_c$ and the velocity dispersion of a galaxy $\sigma$ is used to link the dark matter halo mass and BH mass. The dark matter halo mass function is converted to the BH mass function for any given redshift. The high-redshift optical AGN luminosity functions (LFs) are constructed. At $z\sim 4$, the flattening feature is not shown at the faint end of the optical AGN LF. This is consistent with observational results. If the optical AGN LF at $z\sim 6$ can be reproduced in the case in which central BHs have the Eddington-limited accretion, it is possible for the AGN lifetime to have a small value of $2\times 10^5$ yrs. The X-ray AGN LFs and X-ray AGN number counts are also calculated at $2.0<z<5.0$ and $z>3$, respectively, using the same parameters adopted in the calculation for the optical AGN LF at $z\sim 4$. It is estimated that about 30 AGNs per $\rm{deg}^2$ at $z>6$ can be detected with a flux limit of $3\times 10^{-17}~\rm{erg~cm^{-2}~s^{-1}}$ in the $0.5-2$ keV band. Additionally, the cosmic reionization is also investigated. The ultraviolet photons emitted from the high-redshift AGNs mainly contribute to the cosmic reionization, and the central BHs of the high-redshift AGNs have a mass range of $10^6-10^8M_\odot$. We also discuss some uncertainties in both the AGN LFs and AGN number counts originating from the $M_{\rm{BH}}-\sigma$ relation, Eddington ratio, AGN lifetime, and X-ray attenuation in our model.

The hot intra-cluster medium (ICM) is rich in metals, which are synthesized by supernovae (SNe) explosions and accumulate over time into the deep gravitational potential well of clusters of galaxies. Since most of the elements visible in X-rays are formed by type Ia (SNIa) and/or core-collapse (SNcc) supernovae, measuring their abundances gives us direct information on the nucleosynthesis products of billions of SNe since the epoch of the star formation peak (z ~ 2-3). In this study, we use the EPIC and RGS instruments onboard XMM-Newton to measure the abundances of 9 elements (O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni) from a sample of 44 nearby cool-core galaxy clusters, groups, and elliptical galaxies. We find that the Fe abundance shows a large scatter (~20-40%) over the sample, within 0.2$r_{500}$ and, especially, 0.05$r_{500}$. Unlike the absolute Fe abundance, the abundance ratios (X/Fe) are quite uniform over the considered temperature range (~0.6-8 keV), and with a limited scatter. In addition to a unprecedented treatment of systematic uncertainties, we provide the most accurate abundance ratios measured so far in the ICM, including Cr/Fe and Mn/Fe that we firmly detect (>4{\sigma} with MOS and pn independently). We find that Cr/Fe, Mn/Fe and Ni/Fe, differ significantly from the proto-solar values. However, the large uncertainties in the proto-solar abundances prevent us from making a robust comparison between the local and the intra-cluster chemical enrichments. We also note that, interestingly, and despite the large net exposure time (~4.5 Ms) of our dataset, no line emission feature is seen around ~3.5 keV.

The hot intra-cluster medium (ICM) is rich in metals, which are synthesised by supernovae (SNe) and accumulate over time into the deep gravitational potential well of clusters of galaxies. Since most of the elements visible in X-rays are formed by type Ia (SNIa) and/or core-collapse (SNcc) supernovae, measuring their abundances gives us direct information on the nucleosynthesis products of billions of SNe since the epoch of the star formation peak (z~2-3). In this study, we compare the most accurate average X/Fe abundance ratios (compiled in a previous work from XMM-Newton EPIC and RGS observations of 44 galaxy clusters, groups, and ellipticals), representative of the chemical enrichment in the nearby ICM, to various SNIa and SNcc nucleosynthesis models found in the literature. The use of a SNcc model combined to any favoured standard SNIa model (deflagration or delayed-detonation) fails to reproduce our abundance pattern. In particular, the Ca/Fe and Ni/Fe ratios are significantly underestimated by the models. We show that the Ca/Fe ratio can be reproduced better, either by taking a SNIa delayed-detonation model that matches the observations of the Tycho supernova remnant, or by adding a contribution from the Ca-rich gap transient SNe, whose material should easily mix into the hot ICM. On the other hand, the Ni/Fe ratio can be reproduced better by assuming that both deflagration and delayed-detonation SNIa contribute in similar proportions to the ICM enrichment. In either case, the fraction of SNIa over the total number of SNe (SNIa+SNcc) contributing to the ICM enrichment ranges within 29-45%. This fraction is found to be systematically higher than the corresponding SNIa/SNe fraction contributing to the enrichment of the proto-solar environnement (15-25%). We also discuss and quantify two useful constraints on both SNIa and SNcc that can be inferred from the ICM abundance ratios.

We present the results of $Suzaku$ and $XMM-Newton$ X-ray observations of the cluster pair 1E2216.0-0401 and 1E2215.7-0404. We discover an X-ray bridge between the clusters. $Suzaku$ and $XMM-Newton$ observations revealed that each cluster hosts gas with moderate temperature of $kT_{1E2216.0-0401}=$4.8$\pm$0.1 keV and $kT_{1E2215.7-0404}=$5.8$\pm$0.2 keV, respectively. On the other hand, the bridge region shows a remarkably high temperature ({\it kT}=6.6$\pm$0.5 keV). Furthermore, at the position of the bridge, we detected an enhancement in the wavelet-decomposed soft-band (0.5-4.0 keV) $XMM-Newton$ image at 3 sigma significance, this is most likely due to a compression of the intracluster medium (ICM) as a consequence of the merging activity. This X-ray intensity and temperature enhancement are not consistent with those expected from a late phase, but are in agreement with the predictions by numerical simulations of an early phase merger. From the temperature jump at the location of the bridge, the Mach number is estimated to be ${\cal M}=1.4\pm0.1$, which corresponds to a shock propagation velocity of about 1570 km/s. From the shock properties, we estimate that core-passage will occur in 0.3-0.6 Gyr and that the age of the shock structure is 50--100 Myr. Based on the measured properties of the ICM at the bridge and estimation of timescales, we find indications for non-equilibrium ionization. We also discover possible diffuse radio emission located between the merging clusters. Combining the radio, X-ray, and optical image data, we speculate that the detected radio sources are most likely related to the merger event. Thus, 1E2216.0-0401 and 1E2215.7-0404 is a new example of an early phase cluster merger with remarkable characteristics.

By revisiting the Suzaku and XMM-Newton data of the North Polar Spur, we discovered that the spectra are inconsistent with the traditional model consisting of pure thermal emission and neutral absorption. The most prominent discrepancies are the enhanced O VII and Ne IX forbidden-to-resonance ratios, and a high O VIII Ly$\beta$ line relative to other Lyman series. A collisionally ionized absorption model can naturally explain both features, while a charge exchange component can only account for the former. By including the additional ionized absorption, the plasma in the North Polar Spur can be described by a single-phase CIE component with temperature of 0.25 keV, and nitrogen, oxygen, neon, magnesium, and iron abundances of $0.4-0.8$ solar. The abundance pattern of the North Polar Spur is well in line with those of the Galactic halo stars. The high nitrogen-to-oxygen ratio reported in previous studies can be migrated to the large transmission of the O VIII Ly$\alpha$ line. The ionized absorber is characterized by a balance temperature of $0.17-0.20$ keV and a column density of $3-5 \times 10^{19}$ cm$^{-2}$. Based on the derived abundances and absorption, we speculate that the North Polar Spur is a structure in the Galactic halo, so that the emission is mostly absorbed by Galactic ISM in the line of sight.

In this paper, we scrutinize the effect of spectral index distribution on estimating the AGN (active galactic nucleus) radio luminosity function (RLF) by a Monte Carlo method. We find that the traditional bivariate RLF estimators can cause bias in varying degree. The bias is especially pronounced for the flat-spectrum radio sources whose spectral index distribution is more scattered. We believe that the bias is caused because the $K$-corrections complicate the truncation boundary on the $L-z$ plane of the sample, but the traditional bivariate RLF estimators have difficulty in dealing with this boundary condition properly. We suggest that the spectral index distribution should be incorporated into the RLF analysis process to obtain a robust estimation. This drives the need for a trivariate function of the form $\Phi(\alpha,z,L)$ which we show provides an accurate basis for measuring the RLF.

We examine the scenario that the Doppler factor determines the observational differences of blazars in this paper. Significantly negative correlations are found between the observational synchrotron peak frequency and the Doppler factor. After correcting the Doppler boosting, the intrinsic peak frequency further has a tightly linear relation with the Doppler factor. It is more interesting that this relation is consistent with the scenario that the black hole mass governs both the bulk Lorentz factor and the synchrotron peak frequency. In addition, the distinction of the kinetic jet powers between BL Lacs and FSRQs disappears after the boosting factor $\delta^2$ is considered. The negative correlation between the peak frequency and the observational isotropic luminosity, known as the blazar sequence, also disappears after the Doppler boosting is corrected. We also find that the correlation between the Compton dominance and the Doppler factor exists for all types of blazars. Therefore, this correlation is unsuitable to examine the external Compton emission dominance.

We propose a mixture evolution scenario to model the evolution of the radio luminosity function (RLF) of steep spectrum AGNs (active galactic nuclei) based on a Bayesian method. In this scenario, the shape of the RLF is determined by both the density and luminosity evolution. Our models indicate that the density evolution is positive until a redshift of $\thicksim 0.9$, at which point it becomes negative, while the luminosity evolution is positive to a higher redshift ($z \thicksim 5$ for model B and $z \thicksim 3.5$ for model C), where it becomes negative. Our mixture evolution model works well, and the modeled RLFs are in good agreement with previous determinations. The mixture evolution scenario can naturally explain the luminosity-dependent evolution of the RLFs.

The level-resolved radiative recombination (RR) rate coefficients for H-like to Na-like ions from H (Z = 1) up to and including Zn (Z = 30) are studied here. For H-like ions, the quantum-mechanical exact photoionization cross sections for nonrelativistic hydrogenic systems are used to calculate the RR rate coefficients under the principle of detailed balance, while for He-like to Na-like ions, the archival data on ADAS are adopted. Parameterizations are made for the direct capture rates in a wide temperature range. The fitting accuracies are better than 5% for about 99% of the 30, 000 levels considered here. The 1% exceptions include levels from low- charged many-electron ions, and/or high-shell (n>=4) levels are less important in terms of interpreting X-ray emitting astrophysical plasmas. The RR data will be incorporated into the high-resolution spectral analysis package SPEX.

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($\theta$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

X-ray photons coming from an X-ray point source not only arrive at the detector directly, but also can be strongly forward-scattered by the interstellar dust along the line of sight (LOS), leading to a detectable diffuse halo around the X-ray point source. The geometry of small angle X-ray scattering is straightforward, namely, the scattered photons travel longer paths and thus arrive later than the unscattered ones; thus the delay time of X-ray scattered halo photons can reveal information of the distances of the interstellar dust and the point source. Here we present a study of the X-ray scattered around IGR J17544-2619, which is one of the so-called supergiant fast X-ray transients. IGR J17544-2619 underwent a striking outburst when observed with Chandra on 2004 July 3, providing a near delta-function lightcurve. We find that the X-ray scattered halo around IGR J17544-2619 is produced by two interstellar dust clouds along the LOS. The one which is closer to the observer gives the X-ray scattered at larger observational angles; whereas the farther one, which is in the vicinity of the point source, explains the halo with a smaller angular size. By comparing the observational angle of the scattered halo photons with that predicted by different dust grain models, we are able to determine the normalized dust distance. With the delay times of the scattered halo photons, we can determine the point source distance, given a dust grain model. Alternatively we can discriminate between the dust grain models, given the point source distance.

We present multiwavelength observations of the afterglow of GRB 130427A, the brightest (in total fluence) gamma-ray burst of the past 29 years. Optical spectroscopy from Gemini-North reveals the redshift of the GRB to be z=0.340, indicating that its unprecedented brightness is primarily the result of its relatively close proximity to Earth; the intrinsic luminosities of both the GRB and its afterglow are not extreme in comparison to other bright GRBs. We present a large suite of multiwavelength observations spanning from 300 s to 130 d after the burst and demonstrate that the afterglow shows relatively simple, smooth evolution at all frequencies with no significant late-time flaring or rebrightening activity. The entire dataset from 1 GHz to 10 GeV can be modeled as synchrotron emission from a combination of reverse and forward shocks in good agreement with the standard afterglow model, providing strong support to the applicability of the underlying theory and clarifying the nature of the GeV emission observed to last for minutes to hours following other very bright GRBs. A tenuous, wind-stratified circumburst density profile is required by the observations, suggesting a massive-star progenitor with a low mass-loss rate, perhaps due to low metallicity. GRBs similar in nature to GRB 130427A, inhabiting low-density media and exhibiting strong reverse shocks, are probably not uncommon but may have been difficult to recognize in the past due to their relatively faint late-time radio emission; more such events should be found in abundance by the new generation of sensitive radio and millimeter instruments.

A high-degree of polarization of gamma-ray burst (GRB) prompt emission has been confirmed in recent years. In this paper, we apply jitter radiation to study the polarization feature of GRB prompt emission. In our framework, relativistic electrons are accelerated by turbulent acceleration. Random and small-scale magnetic fields are generated by turbulence. We further determine that the polarization property of GRB prompt emission is governed by the configuration of the random and small-scale magnetic fields. A two-dimensional compressed slab, which contains stochastic magnetic fields, is applied in our model. If the jitter condition is satisfied, the electron deflection angle in the magnetic field is very small and the electron trajectory can be treated as a straight line. A high-degree of polarization can be achieved when the angle between the line of sight and the slab plane is small. Moreover, micro-emitters with mini-jet structure are considered to be within a bulk GRB jet. The jet "off-axis" effect is intensely sensitive to the observed polarization degree. We discuss the depolarization effect on GRB prompt emission and afterglow. We also speculate that the rapid variability of GRB prompt polarization may be correlated with the stochastic variability of the turbulent dynamo or the magnetic reconnection of plasmas.

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

We use an airborne pulsed integrated path differential absorption lidar to make spectroscopic measurements of the pressure-induced line broadening and line center shift of atmospheric CO2 at the 1572.335 nm absorption line. We measure the absorption lineshape in the vertical column between the aircraft and ground. A comparison of our measured absorption lineshape to calculations based on HITRAN shows excellent agreement with the peak optical depth accurate to within 0.3%. Additionally, we measure changes in the line center position to within 5.2 MHz of calculations, and the absorption linewidth to within 0.6% of calculations.

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.

We present our successful Chandra program designed to identify, with sub-arcsecond accuracy, the X-ray afterglow of the short GRB 111117A, which was discovered by Swift and Fermi. Thanks to our rapid target of opportunity request, Chandra clearly detected the X-ray afterglow, though no optical afterglow was found in deep optical observations. The host galaxy was clearly detected in the optical and near-infrared band, with the best photometric redshift of z=1.31_{-0.23}^{+0.46} (90% confidence), making it one of the highest known short GRB redshifts. Furthermore, we see an offset of 1.0 +- 0.2 arcseconds, which corresponds to 8.4 +- 1.7 kpc, between the host and the afterglow position. We discuss the importance of using Chandra for obtaining sub-arcsecond X-ray localizations of short GRB afterglows to study GRB environments.

Jitter radiation, the emission of relativistic electrons in a random and small-scale magnetic field, has been applied to explain the gamma-ray burst (GRB) prompt emission. The seed photons produced from jitter radiation can be scattered by thermal/nonthermal electrons to the high-energy bands. This mechanism is called jitter self-Compton (JSC) radiation. GRB 100728A, which was simultaneously observed by the Swift and Fermi, is a great example to constrain the physical processes of jitter and JSC. In our work, we utilize jitter/JSC radiation to reproduce the multiwavelength spectrum of GRB 100728A. In particular, due to JSC radiation, the powerful emission above the GeV band is the result of those jitter photons in X-ray band scattered by the relativistic electrons with a mixed thermal-nonthermal energy distribution. We also combine the geometric effect of microemitters to the radiation mechanism, such that the "jet-in-jet" scenario is considered. The observed GRB duration is the result of summing up all of the contributions from those microemitters in the bulk jet.

We present the analysis of a large sample of gamma-ray burst (GRB) X-ray light curves in the rest frame to characterise their intrinsic properties in the context of different theoretical scenarios. We determine the morphology, time scales, and energetics of 64 long GRBs observed by \emph{Swift}/XRT \emph{without} flaring activity. We furthermore provide a one-to-one comparison to the properties of GRBs \emph{with} X-ray flares. We find that the steep decay morphology and its connection with X-ray flares favour a scenario in which a central engine origin. We show that this scenario can also account for the shallow decay phase, provided that the GRB progenitor star has a self-similar structure with a constant envelope-to-core mass ratio $\sim 0.02-0.03$. However, difficulties arise for very long duration ($t_p\gtrsim10^4$ s) shallow phases. Alternatively, a spinning-down magnetar whose emitted power refreshes the forward shock can quantitatively account for the shallow decay properties. In particular we demonstrate that this model can account for the plateau luminosity vs. end time anticorrelation.

Context. Multi-wavelength observations of gamma-ray burst (GRB) afterglows provide important information about the activity of their central engines and their environments. In particular, the short timescale variability, such as bumps and/or rebrightening features visible in the multi-wavelength light curves, is still poorly understood. Aims. We analyze the multi-wavelength observations of the GRB100219A afterglow at redshift 4.7. In particular, we attempt to identify the physical origin of the late achromatic flares/bumps detected in the X-ray and optical bands. Methods. We present ground-based optical photometric data and Swift X-ray observations on GRB100219A. We analyzed the temporal behavior of the X-ray and optical light curves, as well as the X-ray spectra. Results. The early flares in the X-ray and optical light curves peak simultaneously at about 1000 s after the burst trigger, while late achromatic bumps in the X-ray and optical bands appear at about 20000 s after the burst trigger. These are uncommon features in the afterglow phenomenology. Considering the temporal and spectral properties, we argue that both optical and X-ray emissions come from the same mechanism. The late flares/bumps may be produced by late internal shocks from long-lasting activity of the central engine. An off-axis origin for a structured jet model is also discussed to interpret the bump shapes. The early optical bump can be interpreted as the afterglow onset, while the early X-ray flare could be caused by the internal activity. GRB 100219A exploded in a dense environment as revealed by the strong attenuation of X-ray emission and the optical-to-X-ray spectral energy distribution.

We present an analysis of the unusual optical light curve of the gamma-ray burst GRB 081029, a long-soft burst with a redshift of z = 3.8479. We combine X-ray and optical observations from the Swift X-Ray Telescope and the Swift UltraViolet/Optical Telescope with ground-based optical and infrared data obtained using the REM, ROTSE, and CTIO 1.3-m telescopes to construct a detailed data set extending from 86 s to approximately 100,000 s after the BAT trigger. Our data cover a wide energy range, from 10 keV to 0.77 eV (1.24 to 16,000 Angstrom). The X-ray afterglow shows a shallow initial decay followed by a rapid decay starting at about 18,000 s. The optical and infrared afterglow, however, shows an uncharacteristic rise at about 3000 s that does not correspond to any feature in the X-ray light curve. Our data are not consistent with synchrotron radiation from a jet interacting with an external medium, a two-component jet, or continuous energy injection from the central engine. We find that the optical light curves can be broadly explained by a collision between two ejecta shells within a two-component jet. A growing number of gamma-ray burst afterglows are consistent with complex jets, which suggests that some (or all) gamma-ray burst jets are complex and will require detailed modelling to fully understand them.

We revisit the radiation mechanism of relativistic electrons in the stochastic magnetic field and apply it to the high-energy emissions of gamma-ray bursts (GRBs). We confirm that jitter radiation is a possible explanation for GRB prompt emission in the condition of a large electron deflection angle. In the turbulent scenario, the radiative spectral property of GRB prompt emission is decided by the kinetic energy spectrum of turbulence. The intensity of the random and small-scale magnetic field is determined by the viscous scale of the turbulent eddy. The microphysical parameters $\epsilon_e$ and $\epsilon_B$ can be obtained. The acceleration and cooling timescales are estimated as well. Due to particle acceleration in magnetized filamentary turbulence, the maximum energy released from the relativistic electrons can reach a value of about $10^{14}$ eV. The GeV GRBs are possible sources of high-energy cosmic-ray.

We present an analysis of the unusual optical light curve of the gamma-ray burst GRB~081029, which occurred at a redshift of z = 3.8479$. We combine X-ray and optical observations from the Swift X-Ray Telescope and the Swift UltraViolet/Optical Telescope with optical and infrared data obtained using the REM and ROTSE telescopes to construct a detailed data set extending from 86 s to approximately 100,000 s after the BAT trigger. Our data also cover a wide energy range, from 10 keV to 0.77 eV (1.24 Angstrom to 16,000 Angstrom). The X-ray afterglow shows a shallow initial decay followed by a rapid decay starting at about 18,000s. The optical and infrared afterglow, however, shows an uncharacteristic rise at about 5000 s that does not correspond to any feature in the X-ray light curve. Our data are not consistent with synchrotron radiation from a single-component jet interacting with an external medium. We do, however, find that the observed light curve can be explained using multi-component model for the jet.

In this paper we shed light on late time (i.e. with peak time t_{pk} \gtrsim 1000 s) flaring activity. We address the morphology and energetic of flares in the window \sim 10^3-10^6 s to put constraints on the temporal evolution of the flare properties and to identify possible differences in the mechanism producing the early and late time flaring emission, if any. This requires the complete understanding of the observational biases affecting the detection of X-ray flares superimposed on a fading continuum at t > 1000 s. We consider all the Swift GRBs that exhibit late time flares. Our sample consists of 36 flares, 14 with redshift measurements. We inherit the strategy of data analysis from Chincarini et al. (2010) in order to make a direct comparison with the early time flare properties. The morphology of the flare light curve is the same for both early time and late time flares, while they differ energetically. The width of late time flares increases with time similarly to the early time flares. Simulations confirmed that the increase of the width with time is not due to the decaying statistics, at least up to 10^4 s. The energy output of late time flares is one order of magnitude lower than the early time flare one, being \sim 1% E_{prompt}. The evolution of the peak luminosity as well as the distribution of the peak flux-to-continuum ratio for late time flares indicate that the flaring emission is decoupled from the underlying continuum, differently from early time flares/steep decay. A sizable fraction of late time flares are compatible with afterglow variability. The internal shock origin seems the most promising explanation for flares. However, some differences that emerge between late and early time flares suggest that there could be no unique explanation about the nature of late time flares.

Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.

Long-duration gamma-ray bursts(LGRBs) are believed to be linked with the star formation. We adopt a galactic evolution model, in which the star formation process inside the virialized dark halo at given redshift can be achieved. In this paper, the gamma-ray burst(GRB) host galaxies are assumed to be the star-forming galaxies within the small dark halos. The star formation rates(SFRs) in the host galaxies of LGRBs at different redshifts have been derived from our model with the galactic evolutionary time about a few times of $10^7$ yr and the dark halo mass of about $5\times 10^{11}M_\odot$. The related stellar masses, luminosities and metallicities of these hosts are estimated as well. We further calculate the X-ray and optical absorption of GRB afterglow emission. From our model calculation, at higher redshift, the SFR of host galaxy is larger, the absorption in X-ray band and optical band of GRB afterglow is stronger, in the condition that the dust and metal components are released locally, surrounding the GRB environment. These model predictions are compared with the {\it Swift} and other observational data. At lower redshift $z<1$, as the merger and interaction of some host galaxies are involved, one monolithic physical process is not sufficient to fully explain all kinds of observed phenomena.

The temporal and spectral analysis of 9 bright X-ray flares out of a sample of 113 flares observed by Swift reveals that the flare phenomenology is strictly analogous to the prompt gamma-ray emission: high energy flare profiles rise faster, decay faster and peak before the low energy emission. However, flares and prompt pulses differ in one crucial aspect: flares evolve with time. As time proceeds flares become wider, with larger peak lag, lower luminosities and softer emission. The flare spectral peak energy E_{p,i} evolves to lower values following an exponential decay which tracks the decay of the flare flux. The two flares with best statistics show higher than expected isotropic energy E_{iso} and peak luminosity L_{p,iso} when compared to the E_{p,i}-E_{iso} and E_{p,i}-L_{iso} prompt correlations. E_{p,i} is found to correlate with L_{iso} within single flares, giving rise to a time resolved E_{p,i}(t)-L_{iso}(t). Like prompt pulses, flares define a lag-luminosity relation: L_{p,iso}^{0.3-10 keV} t_{lag}^{-0.95+/-0.23}. The lag-luminosity is proven to be a fundamental law extending 5 decades in time and 5 in energy. Moreover, this is direct evidence that GRB X-ray flares and prompt gamma-ray pulses are produced by the same mechanism. Finally we establish a flare-afterglow morphology connection: flares are preferentially detected superimposed to one-break or canonical X-ray afterglows.

We present an updated catalog of 113 X-ray flares detected by Swift in the ~33% of the X-ray afterglows of Gamma-Ray Bursts (GRB). 43 flares have a measured redshift. For the first time the analysis is performed in 4 different X-ray energy bands, allowing us to constrain the evolution of the flare temporal properties with energy. We find that flares are narrower at higher energies: their width follows a power-law relation w~E^{-0.5} reminiscent of the prompt emission. Flares are asymmetric structures, with a decay time which is twice the rise time on average. Both time scales linearly evolve with time, giving rise to a constant rise-to-decay ratio: this implies that both time scales are stretched by the same factor. As a consequence, the flare width linearly evolves with time to larger values: this is a key point that clearly distinguishes the flare from the GRB prompt emission. The flare 0.3-10 keV peak luminosity decreases with time, following a power-law behaviour with large scatter: L_{pk}~ t_{pk}^{-2.7}. When multiple flares are present, a global softening trend is established: each flare is on average softer than the previous one. The 0.3-10 keV isotropic energy distribution is a log-normal peaked at 10^{51} erg, with a possible excess at low energies. The flare average spectral energy distribution (SED) is found to be a power-law with spectral energy index beta~1.1. These results confirmed that the flares are tightly linked to the prompt emission. However, after considering various models we conclude that no model is currently able to account for the entire set of observations.

Gamma-ray bursts (GRBSs) are produced by rare types of massive stellar explosions. Their rapidly fading afterglows are often bright enough at optical wavelengths, that they are detectable up to cosmological distances. Hirtheto, the highest known redshift for a GRB was z=6.7, for GRB 080913, and for a galaxy was z=6.96. Here we report observations of GRB 090423 and the near-infrared spectroscopic measurement of its redshift z=8.1^{+0.1}_{-0.3}. This burst happened when the Universe was only ~4% of its current age. Its properties are similar to those of GRBs observed at low/intermediate redshifts, suggesting that the mechanisms and progenitors that gave rise to this burst about 600 million years after the Big Bang are not markedly different from those producing GRBs ~10 billion years later.

We constrain the physical composition of the outflows of GRBs 080916C and 090510 with the prompt emission data and find that the former is likely magnetic while the latter may be baryonic. The X-ray and optical afterglow emission of both GRBs can be reasonably fitted using the standard external shock model but the density profiles of the circum-burst medium are different. We also propose a simple method to estimate the number of the seed photons suppose the GeV afterglow photons are due to the inverse Compton radiation of external forward shock electrons. The seed photons needed in the modeling are too many to be realistic for both events. The synchrotron radiation of the forward shock seems able to account for the GeV afterglow data.

Swift captured for the first time a smoothly rising X-ray re-brightening of clear non-flaring origin after the steep decay in a long gamma-ray burst (GRB): GRB081028. A rising phase is likely present in all GRBs but is usually hidden by the prompt tail emission and constitutes the first manifestation of what is later to give rise to the shallow decay phase. Contemporaneous optical observations reveal a rapid evolution of the injection frequency of a fast cooling synchrotron spectrum through the optical band, which disfavours the afterglow onset (start of the forward shock emission along our line of sight when the outflow is decelerated) as the origin of the observed re-brightening. We investigate alternative scenarios and find that the observations are consistent with the predictions for a narrow jet viewed off-axis. The high on-axis energy budget implied by this interpretation suggests different physical origins of the prompt and (late) afterglow emission. Strong spectral softening takes place from the prompt to the steep decay phase: we track the evolution of the spectral peak energy from the gamma-rays to the X-rays and highlight the problems of the high latitude and adiabatic cooling interpretations. Notably, a softening of both the high and low spectral slopes with time is also observed. We discuss the low on-axis radiative efficiency of GRB081028 comparing its properties against a sample of Swift long GRBs with secure E_gamma,iso measurements.

We present a detailed study of GRB090111, focusing on its extra soft power-law photon index Gamma>5 at the very steep decay phase emission (power-law index alpha=5.1, steeper than 96% of GRBs detected by Swift) and the following peculiar X-ray re-brightening. Our spectral analysis supports the hypothesis of a comoving Band spectrum with the the peak of the nu*F_nu spectrum evolving with time to lower values: a period of higher temporal variability in the 1-2 keV light-curve ends when the E_peak evolves outside the energy band. The X-ray re-brightening shows extreme temporal properties when compared to a homogeneous sample of 82 early flares detected by Swift. While an internal origin cannot be excluded, we show these properties to be consistent with the energy injection in refreshed shocks produced by slow shells colliding with the fastest ones from behind, well after the internal shocks that are believed to give rise to the prompt emission have ceased.

X-ray flashes (XRFs) are a class of gamma-ray bursts (GRBs) with the peak energy of the time-integrated spectrum, Ep, below 30 keV, whereas classical GRBs have Ep of a few hundreds keV. Apart from Ep and the lower luminosity, the properties of XRFs are typical of the classical GRBs. Yet, the nature of XRFs and the differences from that of GRBs are not understood. In addition, there is no consensus on the interpretation of the shallow decay phase observed in most X-ray afterglows of both XRFs and GRBs. We examine in detail the case of XRF 080330 discovered by Swift at the redshift of 1.51. This burst is representative of the XRF class and exhibits an X-ray shallow decay. The rich and broadband (from NIR to UV) photometric data set we collected across this phase makes it an ideal candidate to test the off-axis jet interpretation proposed to explain both the softness of XRFs and the shallow decay phase. We present prompt gamma-ray, early and late IR/visible/UV and X-ray observations of the XRF 080330. We derive a SED from NIR to X-ray bands across the plateau phase with a power-law index of 0.79 +- 0.01 and negligible rest-frame dust extinction. The multi-wavelength evolution of the afterglow is achromatic from ~10^2 s out to ~8x10^4 s. We describe the temporal evolution of the multi-wavelength afterglow within the context of the standard afterglow model and show that a single-component jet viewed off-axis explains the observations (abriged).

Short GRBs are commonly thought to originate from the merging of double compact object binaries but direct evidence for this scenario is still missing. Optical observations of short GRBs allow us to measure redshifts, firmly identify host galaxies, characterize their properties, and accurately localize GRBs within them. Multiwavelength observations of GRB afterglows provide useful information on the emission mechanisms at work. These are all key issues that allow one to discriminate among different models of these elusive events. We carried out photometric observations of the short/hard GRB 051227, GRB 061006, and GRB 071227 with the ESO-VLT starting from several hours after the explosion down to the host galaxy level several days later. For GRB 061006 and GRB 071227 we also obtained spectroscopic observations of the host galaxy. We compared the results obtained from our optical observations with the available X-ray data of these bursts. For all the three above bursts, we discovered optical afterglows and firmly identified their host galaxies. About half a day after the burst, the optical afterglows of GRB 051227 and GRB 061006 present a decay significatly steeper than in the X-rays. In the case of GRB 051227, the optical decay is so steep that it likely indicates different emission mechanisms in the two wavelengths ranges. The three hosts are blue, star forming galaxies at moderate redshifts and with metallicities comparable to the Solar one. The projected offsets of the optical afterglows from their host galaxies centers span a wide range, but all afterglows lie within the light of their hosts and present evidence for local absorption in their X-ray spectra. We discuss our findings in light of the current models of short GRB progenitors.

After considering the effects of negative feedback on the process of star formation, we explore the relationship between star formation process and the associated feedback, by investigating how the mechanical feedback from supernovae(SNe) and radiative feedback from luminous objects regulate the star formation rate and therefore affect the cosmic reionization.Based on our present knowledge of the negative feedback theory and some numerical simulations, we construct an analytic model in the framework of the Lambda cold dark matter model. In certain parameter regions, our model can explain some observational results properly. In large halos(T_vir>10000 K), both mechanical and radiative feedback have a similar behavior: the relative strength of negative feedback reduces as the redshift decreases. In contrast, in small halos (T_vir<10000 K$) that are thought to breed the first stars at early time, the radiative feedback gets stronger when the redshift decreases. And the star formation rate in these small halos depends very weakly on the star-formation efficiency. Our results show that the radiative feedback is important for the early generation stars. It can suppress the star formation rate considerably. But the mechanical feedback from the SNe explosions is not able to affect the early star formation significantly. The early star formation in small-halo objects is likely to be self-regulated. The radiative and mechanical feedback dominates the star formation rate of the PopII/I stars all along. The feedback from first generation stars is very strong and should not be neglected. However, their effects on the cosmic reionization are not significant, which results in a small contribution to the optical depth of Thomson scattering.

We illustrate some of the preliminary results obtained with a new sample of flares and a new analysis. In these proceedings we deal mainly with the analysis related to the flare energy and describe the work in progress to measure the average flare luminosity curve. We discuss in brief GRB050904 and GRB050724 for matters relevant to this work. In particular we measure the contribution given to the flares by GRB050904 and give a new interpretation for the decaying early XRT light curve of GRB050724. We briefly illustrate the first evidence that the early decay is given by the subsequent emission of events with Width/TPeak < 1 and the total energy of these events is larger than the energy emitted during the prompt emission spike showing, indeed, that not only the central engine may still be active after hundreds of seconds of the first spike but that this may still be part of the prompt emission.

The flare activity that is observed in GRBs soon after the prompt emission with the XRT (0.3-10 KeV) instrument on board of the Swift satellite is leading to important clues in relation to the physical characteristics of the mechanism generating the emission of energy in Gamma Ray Bursts. We will briefly refer to the results obtained with the recent analysis and and discuss the preliminary results we obtained with a new larger sample of GRBs [limited to early flares] based on fitting of the flares using the Norris 2005 profile. We find, in agreement with previous results, that XRT flares follow the main characteristics observed in Norris 2005 for the prompt emission spikes. The estimate of the flare energy for the subsample with redshift is rather robust and an attempt is made, using the redshisft sample, to estimate how the energy emitted in flares depends on time. We used a $H_0=70 km/s/Mpc$, $\Omega_\Lambda=0.7$, $\Omega_m=0.3$ cosmology.

We performed the first systematic search for the minimum variability time scale between 0.3 and 10 keV studying the 28 brightest early (<3000 s) afterglows detected by Swift-XRT up to March 2008. We adopt the power spectrum analysis in the time domain: unlike the Fourier spectrum, this is suitable to study the rms variations at different time-scales. We find that early XRT afterglows show variability in excess of the Poissonian noise level on time-scales as short as about 1 s (rest frame value), with the shortest t_{min} associated with the highest energy band. The gamma-ray prompt emission of GRB080319B shows a characteristic average variability time-scale t_{var} of about 1s; this parameter undergoes a remarkable evolution during the prompt emission (BAT observation).

In this Letter, we propose that a microphysical process takes a vital role in the shocked region in which the prompt emission of GRBs is emitted. The turbulent energy is included in the internal energy transferred by the kinetic energy of the shock. It dissipates through stochastic acceleration for the electrons to supply the early X-ray emission in the phase of shallow decay. We put the constraints on the time evolution of microphysical parameters. The early X-ray fluxes can be obtained by this scenario and these results are consistent with the Swift observation

The emission of relativistic electrons moving in the random and small-scale magnetic field is presented by diffusive synchrotron radiation (DSR). In this Letter, we revisit the perturbative treatment of DSR. We propose that random and small-scale magnetic field might be generated by the turbulence. As an example, multi-band radiation of the knot in Cen A comes from the electrons with energy $\gamma_e\sim 10^3-10^4$ in the magnetic field of $10^{-3}G$. The multi-band spectrum of DSR is well determined by the feature of stochastic magnetic field. These results put strong constraint to the models of particle acceleration.

We explore the onset of star formation in the early Universe, exploiting the observations of high-redshift Lyman-break galaxies (LBGs) and Lyman alpha emitters (LAEs), in the framework of the galaxy formation scenario elaborated by Granato et al. (2004) already successfully tested against the wealth of data on later evolutionary stages. Complementing the model with a simple, physically plausible, recipe for the evolution of dust attenuation in metal poor galaxies we reproduce the luminosity functions (LFs) of LBGs and of LAEs at different redshifts. This recipe yields a much faster increase with galactic age of attenuation in more massive galaxies, endowed with higher star formation rates. These objects have therefore shorter lifetimes in the LAE and LBG phases, and are more easily detected in the dusty submillimeter bright phase. The short UV bright lifetimes of massive objects strongly mitigate the effect of the fast increase of the massive halo density with decreasing redshift, thus accounting for the weaker evolution of the LBG LF, compared to that of the halo mass function, and the even weaker evolution between z~6 and z~3 of the LAE LF. LAEs are on the average expected to be younger, with lower stellar masses, and associated to less massive halos than LBGs. Finally, we show that the intergalactic medium can be completely reionized at redshift z~6-7 by massive stars shining in protogalactic spheroids with halo masses from a few 10^10 to a few 10^11 M_sun, showing up as faint LBGs with magnitude in the range -17<M_1350<-20, without resorting to any special stellar initial mass function.

We measure the luminosity function of QSOs in the redshift range 3.5 < z < 5.2 for the absolute magnitude interval -21 < M_{145} < -28. Suitable criteria are defined to select faint QSOs in the GOODS fields, checking their effectiveness and completeness in detail. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. Using a Monte-Carlo technique we estimate the properties of the luminosity function. Our results show that models based on pure density evolution show better agreement with observation than models based on pure luminosity evolution, even if a different break magnitude with respect to z ~ 2.1 is required at 3.5 < z < 5.2. According to our modeling a faint-end slope steeper than low-redshift observations is required to reproduce the data, moreover models with a steep bright-end slope score a higher probability than models with a bright-end flattening. Determining the faint-end of the luminosity function at these redshifts provides important constraints on models of the joint evolution of galaxies and AGNs.

Aims: In this work the luminosity function of QSOs is measured in the redshift range 3.5<z<5.2 for the absolute magnitude interval -21<M_{145}<-28. The determination of the faint end of the luminosity function at these redshifts provides important constraints on models of joint evolution of galaxies and AGNs. Methods: We have defined suitable criteria to select faint QSOs in the GOODS fields, checking in detail their effectiveness and completeness. Spectroscopic follow-up of the resulting QSO candidates has been carried out. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. We have used a Monte-Carlo technique to estimate the properties of the luminosity function, checking various parameterizations for its shape and evolution. Results: Models based on Pure Density Evolution show better agreement with observation than models based on Pure Luminosity Evolution. However a different break magnitude with respect to z~2.1 is required at 3.5<z<5.2. Models with a steeper faint end score a higher probability. We do not find any evidence for a flattening of the bright end at redshift z>3.5. Conclusions: The estimated space density evolution of QSOs indicates a suppression of the formation and/or feeding of Supermassive Black Holes at these redshifts. The QSO contribution to the UV background is insufficient to ionize the IGM at 3.5<z<5.2.

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

The Great Observatories Origins Deep Survey provides significant constraints on the space density of less luminous QSOs at high redshift, which is particularly important to understand the interplay between the formation of galaxies and super-massive black holes and to measure the QSO contribution to the UV ionizing background. We present the results of a search for high-z QSOs, identified in the two GOODS fields on the basis of deep imaging in the optical (with HST) and X-ray (Chandra), and discuss the allowed space density of QSOs in the early universe.

The new results of the evolution of the synchrotron peak for Mrk 421 are mostly likely related to the particle acceleration process. In order to account for the above results, we present a model of blazar variability during the flare in which the emission comes from accelerating electrons. A diffusion advection equation of the electron energy distribution is derived to calculate the spectrum and light curve of synchrotron radiation. We present that the observed shifts of the synchrotron peak moving to higher energies during the flare are caused by shock acceleration. The observed relation between changes in the fluxes at specified frequency ranges and shifts of the peak position is fitted to constrain the physical parameters of the dissipation region.