Observationally, it has been reported that the densest stellar system in the Universe does not exceed a maximum stellar surface density, $\Sigma^{\max}_{*}$ = $3\times10^5$M$_{\odot}$pc$^{-2}$, throughout a wide physical scale ranging from star cluster to galaxy. This suggests there exists a fundamental physics which regulates the star formation and stellar density. However, factors that determine this maximum limit are not clear. In this study, we show that $\Sigma^{\max}_{*}$ of galaxies is not a constant as previous work reported, but actually depends on the stellar mass. We select galaxy sample from the Sloan Digital Sky Survey Data Release 12 at $z=0.01-0.5$. In contrast to a constant maximum predicted by theoretical models, $\Sigma^{\max}_{*}$ strongly depends on stellar mass especially for less massive galaxies with $\sim10^{10}$M$_{\odot}$. We also found that a majority of high-$\Sigma_{*}$ galaxies show red colours and low star-formation rates. These galaxies probably reach the $\Sigma^{\max}_{*}$ as a consequence of the galaxy evolution from blue star forming to red quiescent by quenching star formation. One possible explanation of the stellar-mass dependency of $\Sigma^{\max}_{*}$ is a mass dependent efficiency of stellar feedback. The stellar feedback could be relatively more efficient in a shallower gravitational potential, which terminates star formation quickly before the stellar system reaches a high stellar density.

Fast radio bursts (FRBs) are mysterious radio bursts with a time scale of approximately milliseconds. Two populations of FRB, namely repeating and non-repeating FRBs, are observationally identified. However, the differences between these two and their origins are still cloaked in mystery. Here we show the time-integrated luminosity-duration ($L_{\nu}$-$w_{\rm int,rest}$) relations and luminosity functions (LFs) of repeating and non-repeating FRBs in the FRB Catalogue project. These two populations are obviously separated in the $L_{\nu}$-$w_{\rm int,rest}$ plane with distinct LFs, i.e., repeating FRBs have relatively fainter $L_{\nu}$ and longer $w_{\rm int,rest}$ with a much lower LF. In contrast with non-repeating FRBs, repeating FRBs do not show any clear correlation between $L_{\nu}$ and $w_{\rm int,rest}$. These results suggest essentially different physical origins of the two. The faint ends of the LFs of repeating and non-repeating FRBs are higher than volumetric occurrence rates of neutron-star mergers and accretion-induced collapse (AIC) of white dwarfs, and are consistent with those of soft gamma-ray repeaters (SGRs), type Ia supernovae, magnetars, and white-dwarf mergers. This indicates two possibilities: either (i) faint non-repeating FRBs originate in neutron-star mergers or AIC and are actually repeating during the lifetime of the progenitor, or (ii) faint non-repeating FRBs originate in any of SGRs, type Ia supernovae, magnetars, and white-dwarf mergers. The bright ends of LFs of repeating and non-repeating FRBs are lower than any candidates of progenitors, suggesting that bright FRBs are produced from a very small fraction of the progenitors regardless of the repetition. Otherwise, they might originate in unknown progenitors.

Investigating the Gunn-Peterson trough of high redshift quasars (QSOs) is a powerful way to reveal the cosmic reionization. As one of such attempts, we perform a series of analyses to examine the absorption lines observed with one of the highest redshift QSOs, PSO J006.1240+39.2219, which we previously discovered at z = 6.62. Using the Subaru telescope, we obtained medium-resolution spectrum with a total exposure time of 7.5 hours. We calculate the Ly$\alpha$ transmission in different redshift bins to determine the near zone radius and the optical depth at 5.6$<$z$<$6.5. We find a sudden change in the Ly$\alpha$ transmission at 5.75$<$z$<$5.86, which is consistent with the result from the literature. The near zone radius of the QSO is 5.79$\pm$0.09 $p$Mpc, within the scatter of the near zone radii of other QSOs measured in previous studies. We also analyze the dark gap distribution to probe the neutral hydrogen fractions beyond the saturation limit of the Gunn-Peterson trough. We extend the measurement of the dark gaps to 5.7$<$z$<$6.3. We find that the gap widths increase with increasing redshifts, suggesting more neutral Universe at higher redshifts. However, these measurements strongly depend on the continuum modeling. As a continuum model-free attempt, we also perform the dark-pixel counting analysis, to find the upper limit of $\langle x_{\rm H I} \rangle \sim$0.6 (0.8) at $z<$5.8 ($z>$5.8). All three analyses based on this QSO show increasingly neutral hydrogen towards higher redshifts, adding precious measurements up to z$\sim$6.5.

Faraday rotation measure at radio wavelengths is commonly used to diagnose large-scale magnetic fields. It is argued that the length-scales on which magnetic fields vary in large-scale diffuse astrophysical media can be inferred from correlations in the observed RM. RM is a variable which can be derived from the polarised radiative transfer equations in restrictive conditions. This paper assesses the usage of RMF (rotation measure fluctuation) analyses for magnetic field diagnostics in the framework of polarised radiative transfer. We use models of various magnetic field configurations and electron density distributions to show how density fluctuations could affect the correlation length of the magnetic fields inferred from the conventional RMF analyses. We caution against interpretations of RMF analyses when a characteristic density is ill-defined, e.g. in cases of log-normal distributed and fractal-like density structures. As the spatial correlations are generally not the same in the line-of-sight longitudinal direction and the sky plane direction, one also needs to clarify the context of RMF when inferring from observational data. In complex situations, a covariant polarised radiative transfer calculation is essential to capture all aspects of radiative and transport processes, which would otherwise ambiguate the interpretations of magnetism in galaxy clusters and larger-scale cosmological structures.

Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covariant. It accounts for cosmological and relativistic effects in a self-consistent manner and explicitly treats Faraday rotation, as well as Faraday conversion, emission, and absorption processes. The formulation is derived from the first principles of conservation of phase-space volume and photon number. Without loss of generality, we consider a flat Friedmann-Robertson-Walker (FRW) space-time metric and construct the corresponding polarized radiative transfer equations. We propose an all-sky CPRT calculation algorithm, based on a ray-tracing method, which allows cosmological simulation results to be incorporated and, thereby, model templates of polarization maps to be constructed. Such maps will be crucial in our interpretation of polarized data, such as those to be collected by the Square Kilometer Array (SKA). We describe several tests which are used for verifying the code and demonstrate applications in the study of the polarization signatures in different distributions of electron number density and magnetic fields. We present a pencil-beam CPRT calculation and an all-sky calculation, using a simulated galaxy cluster or a model magnetized universe obtained from GCMHD+ simulations as the respective input structures. The implications on large-scale magnetic field studies are discussed; remarks on the standard methods using rotation measure are highlighted.

The intracluster medium of galaxy clusters is permeated by {\mu}G magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic simulations. We critically assess the current Smoothed Particle Magneto-Hydrodynamics (SPMHD) schemes, especially highlighting the impact of a hyperbolic divergence cleaning scheme and artificial resistivity switch on the magnetic field evolution in cosmological simulations of the formation of a galaxy cluster using the N-body/SPMHD code gcmhd++. The impact and performance of the cleaning scheme and two different schemes for the artificial resistivity switch is demonstrated via idealized test cases and cosmological simulations. We demonstrate that the hyperbolic divergence cleaning scheme is effective at suppressing the growth of the numerical divergence error of the magnetic field and should be applied to any SPMHD simulation. Although the artificial resistivity is important in the strong field regime, it can suppress the growth of the magnetic field in the weak field regime, such as galaxy clusters. With sufficient resolution, simulations with divergence cleaning can reproduce observed magnetic fields. We conclude that the cleaning scheme alone is sufficient for galaxy cluster simulations, but our results indicate that the SPMHD scheme must be carefully chosen depending on the regime of the magnetic field.

High-energy neutrinos and photons are complementary messengers, probing violent astrophysical processes and structural evolution of the Universe. X-ray and neutrino observations jointly constrain conditions in active galactic nuclei (AGN) jets: their baryonic and leptonic contents, and particle production efficiency. Testing two standard neutrino production models for local source Cen A \citep{KT2008,BB2009}, we calculate the high-energy neutrino spectra of single AGN sources and derive the flux of high-energy neutrinos expected for the current epoch. Assuming that accretion determines both X-rays and particle creation, our parametric scaling relations predict neutrino yield in various AGN classes. We derive redshift-dependent number densities of each class, from {\it Chandra} and {\it Swift}/BAT X-ray luminosity functions \citep{SGB2008,ACS2009}. We integrate the neutrino spectrum expected from the cumulative history of AGN (correcting for cosmological and source effects, e.g. jet orientation and beaming). Both emission scenarios yield neutrino fluxes well above limits set by {\it IceCube} (by $\sim 4$--$10^6 \times$ at 1 PeV, depending on the assumed jet models for neutrino production). This implies that: (i) Cen A might not be a typical neutrino source as commonly assumed; (ii) both neutrino production models overestimate the efficiency; (iii) neutrino luminosity scales with accretion power differently among AGN classes and hence does not follow X-ray luminosity universally; (iv) some AGN are neutrino-quiet (e.g. below a power threshold for neutrino production); (v) neutrino and X-ray emission have different duty cycles (e.g. jets alternate between baryonic and leptonic flows); or (vi) some combination of the above.