We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix project with 3 different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes ($M_{sub} > 2.5 \times 10^9h^{-1}M_{\odot}$, or $V_{max} > 70 kms^{-1}$) but more abundant for less massive subhaloes in the NHD simulations. This results in different shapes in the subhalo mass/$V_{max}$ function in two different sets of simulations. At given subhalo mass, the subhaloes less massive than $10^{10} h^{-1}M_{\odot}$ have larger $V_{max}$ in the NHD than DMO simulations, while $V_{max}$ is similar for the subhaloes more massive than the mass value. This is mainly because the progenitors of present day low mass subhaloes have larger concentration parameters in the NHD than DMO simulations. The survival number fraction of the accreted low mass progenitors of the main halo at redshift 2 is about 50 percent higher in the NHD than DMO simulations.

The thermal history of cosmic gas in the Dark Ages remains largely unknown. It is important to quantify the impact of relevant physics on the IGM temperature between $z=10$ and $z \sim 30$, in order to interpret recent and oncoming observations, including results reported by EDGES. We revisit the gas heating due to structure formation shocks in this era, using a set of fixed grid cosmological hydrodynamical simulations performed by three different codes. In all our simulations, the cosmic gas is predicted to be in multiphase state since $z>30$. The gas surrounding high density peaks gradually develops a relation more sharp than $T \propto \rho^{2/3}$, approximately $T \propto \rho^{2}$, from $z=30$ to $z=11$, might due to shock heating. Meanwhile, the gas in void region tends to have a large local mach number, and their thermal state varies significantly from code to code. In the redshift range $11-20$, the mass fraction of gas shock heated above the CMB temperature in our simulations is larger than previous semi-analytical results by a factor of 2 to 8. At $z=15$, the fraction varies from $\sim 19\%$ to $52 \%$ among different codes. Between $z=11$ and $z=20$, the gas temperature $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be $\sim 10-20$ K by two codes, much higher than the adiabatic cooling model and some previous works. However, in our simulations performed by RAMSES, $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be even below the temperature required to explain result of the EDGES. Given the fact that different codes give different predictions, currently, it seems a challenge to make solid prediction on the temperature of gas at $z \sim 17$ in simulations.

In this work, we examine the possibility of using the diffuse supernova neutrino background (DSNB) to test the Chaplygin gas (CG) models of the Universe. With a typical supernova rate $R_{\mathrm{SN}}(z)$ and supernova neutrino spectrum $\mathrm{d} N(E_\nu)/\mathrm{d} E_\nu$, the DSNB flux spectrum $n(E_\nu)$ in three categories of CG models, the generalized CG (GCG), modified CG (MCG) and extended CG (ECG) models, are studied. It is found that generally the flux spectra take a form similar to a Fermi-Dirac distribution with a peak centered around 3.80-3.97 MeV. The spectrum shape and peak positions are primarily determined by $R_{\mathrm{SN}}(z)$ and $\mathrm{d} N(E_\nu)/\mathrm{d} E_\nu$ and only slightly affected by the CG models. However, the height of the spectra in each category of the CG models can vary dramatically for different models, with variances of 13.2%, 23.6% and 14.9% for GCG, MCG and ECG categories respectively. The averaged total flux in each category are also different, with the ECG model average 10.0% and 12.7% higher than that of the GCG and MCG models. These suggest that the DSNB flux spectrum height and total flux can be used to constrain the CG model parameters, and if the measured to a sub-10% accuracy, might be used to rule out some models.

Time delay in Schwarzschild spacetime for null and timelike signals with arbitrary velocity $v$ is studied. The total travel time $t_\mathrm{if}$ is evaluated both exactly and approximately in the weak field limit, with the result given as functions of signal velocity, source-lens and lens-observer distances, angular position of the source and lens mass. Two time delays, $\Delta t_v$ between signals with different velocities but coming from same side of the lens and $\Delta t_\mathrm{p}$ between signals from different sides of the lens, as well as the difference $\Delta t_{\mathrm{p}v}$ between two $\Delta t_\mathrm{p}$'s are calculated. These time delays are applied to the gravitational-lensed supernova neutrinos and gravitational waves (GW). It is shown that the $\Delta t_v$ between different mass eigenstates of supernova neutrinos can be related to the mass square difference of these eigenstates and therefore could potentially be used to discriminate neutrino mass orderings, while the difference $\Delta t_{\mathrm{p}v}$ between neutrino and optical signals can be correlated with the absolute mass of neutrinos. The formula for time delay in a general lens mass profile is derived and the result is applied to the singular isothermal sphere case. For GWs, it is found that the difference $\Delta t_{\mathrm{p}v}$ between GW and GRB can only reach $1.45\times 10^{-5}$ second for very large source distance ($2\times 10^4$ [Mpc]) and source angle (10 [as]) if $v_{GM}=(1-3\times 10^{-15})c$. This time difference is at least three order smaller than uncertainties in time measurement of the recently observed GW/GRB signals and thus calls for improvement if $\Delta t_{\mathrm{p}v}$ is to be used to further constrain the GW velocity.

We present an analysis of the metallicity and star formation activities of HII regions in the interacting system Arp 86, based on the first scientific observation of the multi-object spectroscopy on the 2.16m Telescope at Xinglong Observatory. We find that the oxygen abundance gradient in Arp 86 is flatter than that in normal disk galaxies, which confirms that gas inflows caused by tidal forces during encounters can flatten the metallicity distributions in galaxies. The companion galaxy NGC 7752 is currently experiencing a galaxy-wide starburst with higher surface density of star formation rate than the main galaxy NGC 7753, which can be explained that the companion galaxy is more susceptible to the effects of interaction than the primary. We also find that the galaxy 2MASX J23470758+2926531 has similar abundance and star formation properties to NGC 7753, and may be a part of the Arp 86 system.

Using self-consistent, physically motivated models, we investigate the X-ray obscuration in 19 Type 2 [OIII] 5007 \AA\ selected AGN, 9 of which are local Seyfert 2 galaxies and 10 of which are Type 2 quasar candidates. We derive reliable line-of-sight and global column densities for these objects, which is the first time this has been reported for an AGN sample; 4 AGN have significantly different global and line-of-sight column densities. Five sources are heavily obscured to Compton-thick. We comment on interesting sources revealed by our spectral modeling, including a candidate ``naked'' Sy2. After correcting for absorption, we find that the ratio of the rest-frame, 2-10 keV luminosity (L$_{\rm 2-10keV,in}$) to L$_{\rm [OIII]}$ is 1.54 $\pm$ 0.49 dex which is essentially identical to the mean Type 1 AGN value. The Fe K$\alpha$ luminosity is significantly correlated with L$_{\rm [OIII]}$, but with substantial scatter. Finally, we do not find a trend between L$_{\rm 2-10keV,in}$ and global or line-of-sight column density, between column density and redshift, between column density and scattering fraction or between scattering fraction and redshift.

(Abridged) We analyzed the {\it Chandra} and {\it XMM-Newton} archival observations for 72 type 2 quasars at $z<1$. These objects were was selected based on the [O III]$\lambda$5007 optical emission line which we assume to be an approximate indicator of the intrinsic AGN luminosity. We find that the means of the column density and photon index of our sample are $\log N_{\rm H}=23.0$ cm$^{-2}$ and $\Gamma=1.87$ respectively, which are consistent with results from deep X-ray surveys. The observed ratios of hard X-ray and [O III] line luminosities imply that the majority of our sample suffer significant amounts of obscuration in the hard X-ray band. A more physically realistic model which accounts for both Compton scattering and a potential partial covering of the central X-ray source was used to estimate the true absorbing column density. We find that the absorbing column density estimates based on simple power-law models significantly underestimate the actual absorption in approximately half of the sources. Eleven sources show a prominent Fe K$\alpha$ emission line, and we detect this line in the other sources through a joint fit (spectral stacking). The correlation between the Fe K$\alpha$ and [O III] fluxes and the inverse correlation of the equivalent width of Fe K$\alpha$ line with the ratio of hard X-ray and [O III] fluxes is consistent with previous results for lower luminosity Seyfert 2 galaxies. We conclude that obscuration is the cause of the weak hard X-ray emission rather than intrinsically low X-ray luminosities. We find that about half of the population of optically-selected type 2 quasars are likely to be Compton-thick. We also find no evidence that the amount of X-ray obscuration depends on the AGN luminosity (over a range of more than three orders-of-magnitude in luminosity).

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

We present a {\it Chandra} observation of IRAS 19254--7245, a nearby ULIRG also known as {\it the Superantennae}. The high spatial resolution of {\it Chandra} allows us to disentangle for the first time the diffuse starburst emission from the embedded Compton-thick AGN. The 2-10 keV spectrum of the AGN emission is fitted by a flat power-law $\Gamma=1.3$) and a He-like Fe K$\alpha$ line with EW$\sim$1.5 keV, consistent with previous observations. The Fe K$\alpha$ line profile could be resolved as a blend of a neutral 6.4 keV line and an ionized 6.7 keV (He-like) or 6.9 keV (H-like) line. Variability is detected compared with the previous {\it XMM-Newton} and {\it suzaku} observations, demonstrating the compact size of the iron line emission. We fit the spectrum of the galaxy-scale extended emission excluding the AGN and other bright point sources with a soft thermal component with kT~0.8 keV. The luminosity of the extended emission is about one order of magnitude lower than that of the AGN. The basic physical and structural properties of the extended emission are fully consistent with a galactic wind being driven by the starburst (no contribution to the feedback by the AGN is required). A candidate ultra-luminous X-ray source is detected 8\arcsec\ south of the southern nucleus. The 0.3-10 keV luminosity of this off-nuclear point source is ~$6\times 10^{40}$ erg s$^{-1}$ if the emission is isotropic and the source is associated with the Superantennae.

We have used XMM-Newton to observe six Lyman Break Analogs (LBAs): members of the rare population of local galaxies that have properties that are very similar to distant Lyman Break Galaxies. Our six targets were specifically selected because they have optical emission-line properties that are intermediate between starbursts and Type 2 (obscured) AGN. Our new X-ray data provide an important diagnostic of the presence of an AGN. We find X-ray luminosities of order 10^{42} erg/s and ratios of X-ray to far-IR luminosities that are higher than values in pure starburst galaxies by factors ranging from ~ 3 to 30. This strongly suggests the presence of an AGN in at least some of the galaxies. The ratios of the luminosities of the hard (2-10 keV) X-ray to [O III]\lambda 5007 emission-line are low by about an order-of-magnitude compared to Type 1 AGN, but are consistent with the broad range seen in Type 2 AGN. Either the AGN hard X-rays are significantly obscured or the [O III] emission is dominated by the starburst. We searched for an iron emission line at ~ 6.4 keV, which is a key feature of obscured AGN, but only detected emission at the ~ 2\sigma level. Finally, we find that the ratios of the mid-infrared (24\mu m) continuum to [O III]\lambda 5007 luminosities in these LBAs are higher than the values for Type 2 AGN by an average of 0.8 dex. Combining all these clues, we conclude that an AGN is likely to be present, but that the bolometric luminosity is produced primarily by an intense starburst. If these black holes are radiating at the Eddington limit, their masses would lie in the range of 10^5 to 10^6 M_{sun}. These objects may offer ideal local laboratories to investigate the processes by which black holes grew in the early universe.

The Copernican principle, stating that we do not occupy any special place in our universe, is usually taken for granted in modern cosmology. However recent observational data of supernova indicate that we may live in the under-dense center of our universe, which makes the Copernican principle challenged. It thus becomes urgent and important to test the Copernican principle via cosmological observations. Taking into account that unlike the cosmic photons, the cosmic neutrinos of different energies come from the different places to us along the different worldlines, we here propose cosmic neutrino background as a test of the Copernican principle. It is shown that from the theoretical perspective cosmic neutrino background can allow one to determine whether the Copernican principle is valid or not, but to implement such an observation the larger neutrino detectors are called for.

We present a geometrical methodology to interpret the periodical light curves of Soft Gamma Repeaters based on the magnetar model and the numerical arithmetic of the three-dimensional magnetosphere model for the young pulsars. The hot plasma released by the star quake is trapped in the magnetosphere and photons are emitted tangent to the local magnetic field lines. The variety of radiation morphologies in the burst tails and the persistent stages could be well explained by the trapped fireballs on different sites inside the closed field lines. Furthermore, our numerical results suggests that the pulse profile evolution of SGR 1806-20 during the 27 December 2004 giant flare is due to a lateral drift of the emitting region in the magnetosphere.

We use a modified outer gap model to study the multi-frequency phase-resolved spectra of the Crab pulsar. The emissions from both poles contribute to the light curve and the phase-resolved spectra. Using the synchrotron self-Compton mechanism and by considering the incomplete conversion of curvature photons into secondary pairs, the observed phase-averaged spectrum from 100 eV - 10 GeV can be explained very well. The predicted phase-resolved spectra can match the observed data reasonably well, too. We find that the emission from the north pole mainly contributes to Leading Wing 1. The emissions in the remaining phases are mainly dominated by the south pole. The widening of the azimuthal extension of the outer gap explains Trailing Wing 2. The complicated phase-resolved spectra for the phases between the two peaks, namely Trailing Wing 1, Bridge and Leading Wing 2, strongly suggest that there are at least two well-separated emission regions with multiple emission mechanisms, i.e. synchrotron radiation, inverse Compton scattering and curvature radiation. Our best fit results indicate that there may exist some asymmetry between the south and the north poles. Our model predictions can be examined by GLAST.

We present a modified outer gap model to study the phase-resolved spectra of the Crab pulsar. A theoretical double peak profile of the light curve containing the whole phase is shown to be consistent with the observed light curve of the Crab pulsar by shifting the inner boundary of the outer gap inwardly to $\sim 10$ stellar radii above the neutron star surface. In this model, the radial distances of the photons corresponding to different phases can be determined in the numerical calculation. Also the local electrodynamics, such as the accelerating electric field, the curvature radius of the magnetic field line and the soft photon energy, are sensitive to the radial distances to the neutron star. Using a synchrotron self-Compton mechanism, the phase-resolved spectra with the energy range from 100 eV to 3 GeV of the Crab pulsar can also be explained.

We report on the emission properties of PSR B1929+10 and its putative X-ray trail from a multi-wavelength study performed with XMM-Newton, the ESO NTT, the HST, the Effelsberg 100m Radio Telescope and the Jodrell Bank Radio Observatory. The XMM-Newton observations confirm the existence of the diffuse emission with a trail morphology lying in a direction opposite to the transverse motion of the pulsar. The trail has a length of ~15 arcmin. Its spectrum is non-thermal and produced by electron-synchrotron emission in the shock between the pulsar wind and the surrounding medium. Assuming that the electron lifetime against synchrotron cooling is comparable to the source transit time over the X-ray trail length, the magnetic field strength in the trail emitting region is inferred to be ~5 uG. Inspecting data from the Effelsberg 11cm radio continuum survey of the Galactic plane we discovered an elongated feature apparently coincident with the X-ray trail. The emission properties observed from PSR 1929+10 are found to finally challenge the picture of the emission properties of old non-recycled rotation powered pulsars. Both, the temporal and spectral X-ray emission properties of PSR 1929+10 are in excellent agreement with a non-thermal and, thus, magnetospheric radiation dominated emission scenario. A flux contribution from the thermal emission of heated polar caps of ~7% is inferred from a best fitting composite Planckian and power law spectral model. The X-ray pulse profile is found to be markedly different from the broad sinusoidal pulse profile seen in the low statistic ROSAT data. Simulations in the framework of an outer-gap emission model are able to reproduce the observed X-ray pulse profile and its phase relation relative to the radio pulse.