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.

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the Schechter function, galaxies were investigated. The cross-correlation function with the surrounding galaxies was measured and compared between the AGN and blue galaxy samples. We also compared the distributions of color, absolute magnitude, and stellar mass of the galaxies around such target objects. The properties of clusters detected using surrounding galaxies selected based on a photometric redshift were examined and compared for different samples. The target AGNs were drawn from the Million Quasars (MILLIQUAS) catalog, and the blue galaxies were drawn from six redshift survey catalogs (SDSS, WiggleZ, DEEP2, VVDS, VIPERS, and PRIMUS). The galaxies used as a measure of the environment around the targets are drawn from S18a internal data released by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We found that, among the five AGN and blue galaxy samples considered, the environment of AGN-LG pairs is the most enriched with luminous galaxies. We also found an enhancement in the number of mass-selected clusters in the AGN-LG pair sample against those in the other samples. The results obtained in this study indicate that existence of multiple clusters is the major driver in the association of AGNs and LGs, rather than a single large-mass dark matter halo hosting the AGN.

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.

We investigate infrared colours and spectral energy distributions (SEDs) of 338 X-ray selected AGN from Swift-BAT 105-month survey catalogue that have AKARI detection, in order to find a new selection criteria for Compton-thick AGN. By combining data from Galaxy Evolution Explore (GALEX), Sloan Digital Sky Survey (SDSS) Data Release 14 (DR14), Two Micron All Sky Survey (2MASS), Wide-field Infrared Survey Explorer (WISE), AKARI and Herschel for the first time we perform ultraviolet (UV) to far-infrared (FIR) SED fitting 158 Swift BAT AGN by CIGALE and constrain the AGN model parameters of obscured and Compton-thick AGN. The comparison of average SEDs show while the mid-IR (MIR) SEDs are similar for the three AGN populations, optical/UV and FIR regions have differences. We measure the dust luminosity, the pure AGN luminosity and the total infrared (IR) luminosity. We examine the relationships between the measured infrared luminosities and the hard X-ray luminosity in the 14-195 keV band. We show that the average covering factor of Compton-thick AGN is higher compared to the obscured and unobscured AGN. We present a new infrared selection for Compton-thick AGN based on MIR and FIR colours ([9$\mu$m - 22$\mu$m]$ > 3.0$ and [22$\mu$m - 90$\mu$m]$ < 2.7$) from WISE and AKARI. We find two known Compton-thick AGN that are not included in the Swift-BAT sample, and conclude that MIR colours covering 9.7$\mu$m silicate absorption and MIR continuum can be a promising new tool to identify Compton-thick AGN.

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.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

We first present new Subaru narrow-band observations of the Ly$\alpha$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Ly\alpha}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$\alpha$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$\alpha$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$\alpha$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$\alpha$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~\mu$m wavelength range is successfully detected with $>5\sigma$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $\mu$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.

We report on the torus constraints of the Compton-thick AGN with double-peaked optical narrow line region (NLR) emission lines, ANEPD-CXO245, at z=0.449 in the AKARI NEP Deep Field. The unique infrared data on this field, including those from the nine-band photometry over 2-24 $\mu$m with the AKARI Infrared Camera (IRC), and the X-ray spectrum from Chandra allow us to constrain torus parameters such as the torus optical depth, X-ray absorbing column, torus angular width ($\sigma$) and viewing angle ($i$). We analyze the X-ray spectrum as well as the UV-optical-infrared spectral energy distribution (UOI-SED) with clumpy torus models in X-ray (XCLUMPY; Tanimoto et al. 2019) and infrared (CLUMPY; Nenkova et al. 2008) respectively. From our current data, the constraints on $\sigma$--$i$ from both X-rays and UOI show that the line of sight crosses the torus as expected for a type 2 AGN. We obtain a small X-ray scattering fraction (<0.1%), which suggests narrow torus openings, giving preference to the bi-polar outflow picture of the double-peaked profile. Comparing the optical depth of the torus from the UOI-SED and the absorbing column density $N_{\rm H}$ from the X-ray spectrum, we find that the gas-to-dust ratio is $\gtrsim 4$ times larger than the Galactic value.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$\sigma$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

Galaxies evolve from a blue star-forming phase into a red quiescent one by quenching their star formation activity. In high density environments, this galaxy evolution proceeds earlier and more efficiently. Therefore, local galaxy clusters are dominated by well-evolved red, elliptical galaxies. The fraction of blue galaxies in clusters monotonically declines with decreasing redshift, i.e., the Butcher-Oemler effect. In the local Universe, observed blue fractions of massive clusters are as small as $\lesssim$ 0.2. Here we report a discovery of a \lq \lq blue cluster\rq \rq, that is a local galaxy cluster with an unprecedentedly high fraction of blue star-forming galaxies yet hosted by a massive dark matter halo. The blue fraction is 0.57, which is 4.0 $\sigma$ higher than those of the other comparison clusters under the same selection and identification criteria. The velocity dispersion of the member galaxies is 510 km s$^{-1}$, which corresponds to a dark matter halo mass of 2.0$^{+1.9}_{-1.0}\times 10^{14}$ M$_{\odot}$. The blue fraction of the cluster is more than 4.7 $\sigma$ beyond the standard theoretical predictions including semi-analytic models of galaxy formation. The probability to find such a high blue fraction in an individual cluster is only 0.003\%, which challenges the current standard frameworks of the galaxy formation and evolution in the $\Lambda$CDM Universe. The spatial distribution of galaxies around the blue cluster suggests that filamentary cold gas streams can exist in massive halos even in the local Universe. However these cold streams have already disappeared in the theoretically simulated local universes.

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.

Nature of dark energy remains unknown. Especially, to constrain the time variability of the dark-energy, a new, standardisable candle that can reach more distant Universe has been awaited. Here we propose a new distance measure using fast radio bursts (FRBs), which are a new emerging population of $\sim$ ms time scale radio bursts that can reach high-$z$ in quantity. We show an empirical positive correlation between the time-integrated luminosity (L$_{\nu}$) and rest-frame intrinsic duration ($w_{\rm int,rest}$) of FRBs. The L$_{\nu}-w_{\rm int,rest}$ correlation is with a weak strength but statistically very significant, i.e., Pearson coefficient is $\sim$ 0.5 with p-value of $\sim$0.038, despite the smallness of the current sample. This correlation can be used to measure intrinsic luminosity of FRBs from the observed $w_{\rm int,rest}$. By comparing the luminosity with observed flux, we measure luminosity distances to FRBs, and thereby construct the Hubble diagram. This FRB cosmology with the L$_{\nu}-w_{\rm int,rest}$ relation has several advantages over SNe Ia, Gamma-Ray Burst (GRB), and well-known FRB dispersion measure (DM)-$z$ cosmology; (i) access to higher redshift Universe beyond the SNe Ia, (ii) high event rate that is $\sim$ 3 order of magnitude more frequent than GRBs, and (iii) it is free from the uncertainty from intergalactic electron density models, i.e., we can remove the largest uncertainty in the well-debated DM-$z$ cosmology of FRB. Our simulation suggests that the L$_{\nu}-w_{\rm int,rest}$ relation provides us with useful constraints on the time variability of the dark energy when the next generation radio telescopes start to find FRBs in quantity.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly.

The aim of this work is to create a new catalog of reliable AGN candidates selected from the AKARI NEP-Deep field. Selection of the AGN candidates was done by applying a fuzzy SVM algorithm, which allows to incorporate measurement uncertainties into the classification process. The training dataset was based on the spectroscopic data available for selected objects in the NEP-Deep and NEP-Wide fields. The generalization sample was based on the AKARI NEP-Deep field data including objects without optical counterparts and making use of the infrared information only. A high quality catalog of previously unclassified 275 AGN candidates was prepared.

A linear correlation has been proposed between the CO luminosity ($\rm{L}^{\prime}_{\rm{CO}}$) and full-width at half maximum (FWHM) for high-redshift (z > 1) submillimeter galaxies. However, the controversy concerning the $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlation seems to have been caused by the use of heterogeneous samples (e.g., different transition lines) and/or data with large measurement uncertainties. In order to avoid the uncertainty caused by using different rotational transitions, in this work we make an extensive effort to select only CO($J = 1-0$) data from the literature. We separate these wide-ranging redshift data into two samples : the low-redshift (z < 1) and high-redshift (z > 1) samples. The samples are corrected for lensing magnification factors if gravitational-lensing effects appeared in the observations. The correlation analysis shows that there exists significant $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations for both the low-redshift and high-redshift samples. A comparison of the low- and high-redshift $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations does not show strong evolution with redshift. Assuming that there is no evolution, we can use this relation to determine the model independent distances of high-redshift galaxies. We then constrain cosmological models with the calibrated high-redshift CO data and the sample of Type Ia supernovae in the Union 2.1 compilation. In the constraint for wCDM with our samples, the derived values are w_{0} = -1.02 {\pm} 0.17, {\Omega}_{m0} = 0.30{\pm}0.02, and H_{0} = 70.00 {\pm}0.60 km\,s^{-1}\,Mpc^{-1}.

Evolutionary properties of infrared (IR) luminous galaxies are important keys to understand dust-obscured star formation history and galaxy evolution. Based on the near- to mid-IR imaging with 9 continuous filters of AKARI space telescope, we present the characteristics of dusty star-forming (SF) galalxies showing polycyclic aromatic hydrocarbon (PAH) features observed by the North Ecliptic Pole (NEP) wide field survey of AKARI and Herschel. All the sample galaxies from the AKARI/NEP-Wide data are selected based both on the Herschel/SPIRE 250 {\mu}m detection and optical spectroscopic redshift data. The physical modelling of spectral energy distribution (SED) using all available data points from u to sub-mm 500 {\mu}m band, including WISE and PACS data where available, takes unique advantages of the continuous near to mid-IR coverage, reliable constraint on far-IR peak, spectroscopically determined accurate redshifts, as well as energy balance principle by MAGPHYS. This enables us to derive physically meaningful and accurate total infrared luminosity and 8 {\mu}m (or PAH) luminosity consistently. Our sample galaxies are in the redshift range z <1, and majority of them appear to be normal SF/spiral populations showing PAH features near the 8 {\mu}m. These SF galaxies showing PAHs in the mid-IR include various types from quiescent to starbursts. Some of our sample show shortage of 8 {\mu}m luminosity compared to the total IR luminosity and this PAH deficit gets severe in more luminous IR galaxies, suggesting PAH molecules in these galaxies destroyed by strong radiation field from SF region or a large amount of cold dust in ISM. The specific SFR of our sample shows mass dependent time evolution which is consistent with downsizing evolutionary pattern.

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

Revealing what fraction of galaxies harbor AGN is central in understanding black hole accretion history of the Universe. However, optical and soft X-ray surveys miss the most highly obscured AGNs. Infrared (IR), instead, is more robust against absorption. Previous IR photometric surveys, however, only had 4 or 5 filters in mid-IR. Our AKARI North Ecliptic Pole (NEP) wide field sample has 18 filters in mid-IR (9 from AKARI, 4 from WISE, and 5 from Spitzer), for the first time, allowing a sophisticated mid-IR SED fitting diagnosis for a statistical number of sources (89178 over 5.4 deg$^2$). By using a SED fitting technique, we investigate the evolution of AGN fraction as a function of redshift and IR (8-1000 $\mu$m) luminosity in an extinction-free way. We found that the AGN fraction (F$_{\rm AGN}$) shows no sign of strong redshift evolution. Instead, F$_{\rm AGN}$ increases with increasing IR luminosity in all redshifts bins ($0<z<2$).

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed.

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

We present the analysis of a new near-infrared (NIR) spectrum of a recently discovered $z=6.621$ quasar PSO J006+39 in an attempt to explore the early growth of supermassive black holes (SMBHs). This NIR (rest-frame ultraviolet, UV) spectrum shows blue continuum slope and rich metal emission lines in addition to Ly$\alpha$ line. We utilize the MgII line width and the rest frame luminosity $L_\text{3000\AA}$ to find the mass of SMBH ($M_\text{BH}$) to be $\sim 10^8 M_\odot$, making this one of the lowest mass quasars at high redshift. The power-law slope index ($\alpha_\lambda$) of the continuum emission is $-2.94\pm0.03$, significantly bluer than the slope of $\alpha_\lambda=-7/3$ predicted from standard thin disc models. We fit the spectral energy distribution (SED) using a model which can fit local SMBHs, which includes warm and hot Comptonisation powered by the accretion flow as well as an outer standard disc. The result shows that the very blue slope is probably produced by a small radial ($\sim230$ gravitational radius, $R_\text{g}$) extent of the standard accretion disc. All plausible SED models require that the source is super-Eddington ($L_\text{bol}/L_\text{Edd} \gtrsim 9$), so the apparently small disc may simply be the inner funnel of a puffed up flow, and clearly the SMBH in this quasar is in a rapid growth phase. We also utilize the rest-frame UV emission lines to probe the chemical abundance in the broad line region (BLR) of this quasar. We find that this quasar has super solar metallicity through photoionization model calculations.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9{\deg} for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1{\sigma} and with the Auger result to within 1.4{\sigma}, meaning that it is not capable to discriminate between these two hypotheses.

A detailed analysis of Herschel-PACS observations at the North Ecliptic Pole is presented. High quality maps, covering an area of 0.44 square degrees, are produced and then used to derive potential candidate source lists. A rigorous quality control pipeline has been used to create final legacy catalogues in the PACS Green 100 micron and Red 160 micron bands, containing 1384 and 630 sources respectively. These catalogues reach to more than twice the depth of the current archival Herschel/PACS Point Source Catalogue, detecting 400 and 270 more sources in the short and long wavelength bands respectively. Galaxy source counts are constructed that extend down to flux densities of 6mJy and 19mJy (50% completeness) in the Green 100 micron and Red 160 micron bands respectively. These source counts are consistent with previously published PACS number counts in other fields across the sky. The source counts are then compared with a galaxy evolution model identifying a population of luminous infrared galaxies as responsible for the bulk of the galaxy evolution over the flux range (5-100mJy) spanned by the observed counts, contributing approximate fractions of 50% and 60% to the cosmic infrared background (CIRB) at 100 microns and 160 microns respectively.

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 structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence.

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution.

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.

Star-forming galaxies display a close relation among stellar mass, metallicity and star-formation rate (or molecular-gas mass). This is known as the fundamental metallicity relation (FMR) (or molecular-gas FMR), and it has a profound implication on models of galaxy evolution. However, there still remains a significant residual scatter around the FMR. We show here that a fourth parameter, the surface density of stellar mass, reduces the dispersion around the molecular-gas FMR. In a principal component analysis of 29 physical parameters of 41,338 star-forming galaxies, the surface density of stellar mass is found to be the fourth most important parameter. The new four-dimensional fundamental relation forms a tighter hypersurface that reduces the metallicity dispersion to 50\% of that of the molecular-gas FMR. We suggest that future analyses and models of galaxy evolution should consider the FMR in a four-dimensional space that includes surface density. The dilution time scale of gas inflow and the star-formation efficiency could explain the observational dependence on surface density of stellar mass.

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshift; the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in 9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, we measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe much more precisely than previous work. In the AKARI NEP wide field, AKARI has obtained deep images in the mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage. To rectify the situation, we used the newly advent Subaru telescope's Hyper Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the AKARI NEP wide field. With this deep and wide optical data, we, for the first time, can use the entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from z=0 to z=2.2, all probed by the AKARI satellite.

Star-forming galaxies display a close relation among stellar mass, metallicity and star-formation rate (or molecular-gas mass). This is known as the fundamental metallicity relation (FMR) (or molecular-gas FMR), and it has a profound implication on models of galaxy evolution. However, there still remains a significant residual scatter around the FMR. We show here that a fourth parameter, the surface density of stellar mass, reduces the dispersion around the molecular-gas FMR. In a principal component analysis of 29 physical parameters of 41,338 star-forming galaxies, the surface density of stellar mass is found to be the fourth most important parameter. The new four-dimensional (4D) fundamental relation forms a tighter hypersurface that reduces the metallicity dispersion to 50% of that of the molecular-gas FMR. We suggest that future analyses and models of galaxy evolution should consider the FMR in a 4D space that includes surface density. The dilution time scale of gas inflow and the star-formation efficiency could explain the observational dependence on surface density of stellar mass. AKARI is expected to play an important role in shedding light on the infrared properties of the new 4D FMR.

Local infrared (IR) luminosity functions (LFs) are necessary benchmarks for high-redshift IR galaxy evolution studies. Any accurate IR LF evolution studies require accordingly accurate local IR LFs. We present infrared galaxy LFs at redshifts redshifts of $z \leq 0.3$ from AKARI space telescope, which performed an all-sky survey in six IR bands (9, 18, 65, 90, 140 and 160 micron) with 10 times better sensitivity than its precursor IRAS. Availability of 160 micron filter is critically important in accurately measuring total IR luminosity of galaxies, covering across the peak of the dust emission. By combining data from Wide-field Infrared Survey Explorer (WISE), Sloan Digital Sky Survey (SDSS) Data Release 13 (DR13), 6-degree Field Galaxy Survey (6dFGS) and the 2MASS Redshift Survey (2MRS), we created a sample of 15,638 local IR galaxies with spectroscopic redshifts, factor of 7 larger compared to previously studied AKARI -SDSS sample. After carefully correcting for volume effects in both IR and optical, the obtained IR LFs agree well with previous studies, but comes with much smaller errors. Measured local IR luminosity density is $\Omega_{IR}=$ 1.19$\pm$0.05 $\times 10^{8}$ L$_{\odot}$ Mpc$^{-3}$. The contributions from luminous infrared galaxies and ultra luminous infrared galaxies to IR are very small, 9.3 per cent and 0.9 per cent, respectively. There exists no future all sky survey in far-infrared wavelengths in the foreseeable future. The IR LFs obtained in this work will therefore remain an important benchmark for high-redshift studies for decades.

The AKARI North Ecliptic Pole (NEP) survey consists of two survey projects: NEP-Deep (0.5 sq.deg) and NEP-Wide (5.4 sq.deg), providing with tens of thousands of galaxies. A continuous filter coverage in the mid-infrared wavelengths (7, 9, 11, 15, 18 and 24 $\mu$m) is unique to diagnose the contributions from dusty star-formation activity and AGNs. Here we present current status focused on the newly obtained optical images and near-future prospects with a new X-ray telescope. Hyper Suprime-Cam on Subaru telescope is a gigantic optical camera with huge Field of View (FoV). Thanks to the wide FoV, we successfully obtained deep optical images at g, r, i, z and Y-bands covering most of the NEP-Wide field. Using the deep optical images, we identified over 5000 optical counterparts of the mid-IR sources, presumably deeply obscured galaxies in NEP-Wide field. We also investigated properties of these infrared sources with SED-fitting. eROSITA, to be launched early 2018, is a new all-sky X-ray survey telescope, and expected to conduct ultra deep 2-10 keV imaging toward NEP. We expect unprecedentedly numerous Compton-thick AGN candidates when combined with the multi-wavelength data in NEP region.

In this research, we provide a new, efficient method to select infrared (IR) active galatic nucleus (AGN). In the past, AGN selection in IR had been established by many studies using color-color diagrams. However, those methods have a problem in common that the number of bands is limited. The AKARI North Ecliptic Pole (NEP) survey was carried out by the AKARI Infrared Camera (IRC), which has 9 filters in mid-IR with a continuous wavelength coverage from 2 to 24$\mu$m$^{-1}$. Based on the intrinsic different mid-IR features of AGN and star-forming galaxies (SFGs), we performed SED fitting to separate these two populations by the best-fitting model. In the X-ray AGN sample, our method by SED fitting selects 50$\%$ AGNs, while the previous method by colour criteria recovers only 30$\%$ of them, which is a significant improvement. Furthermore, in the whole NEP deep sample, SED fitting selects two times more AGNs than the color selection. This may imply that the black hole accretion history could be more stronger than people expected before.

The recent updates of the North Ecliptic Pole deep (0.5~deg$^2$, NEP-Deep) multi-wavelength survey covering from X-ray to radio-wave is presented. The NEP-Deep provides us with several thousands of 15~$\mu$m or 18~$\mu$m selected sample of galaxies, which is the largest sample ever made at this wavelengths. A continuous filter coverage in the mid-infrared wavelength (7, 9, 11, 15, 18, and 24~$\mu$m) is unique and vital to diagnose the contributions from starbursts and AGNs in the galaxies out to $z$=2.The new goal of the project is to resolve the nature of the cosmic star formation history at the violent epoch (e.g. $z$=1--2), and to find a clue to understand its decline from $z$=1 to present universe by utilizing the unique power of the multi-wavelength survey. The progress in this context is briefly mentioned.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.

We conduct a systematic search for galaxy protoclusters at $z\sim3.8$ based on the latest internal data release (S16A) of the Hyper SuprimeCam Subaru strategic program (HSC-SSP). In the Wide layer of the HSC-SSP, we investigate the large-scale projected sky distribution of $g$-dropout galaxies over an area of $121\,\mathrm{deg^2}$, and identify 216 large-scale overdense regions ($>4\sigma$ overdensity significance) that are good protocluster candidates. Of these, 37 are located within $8\,\mathrm{arcmin}$ ($3.4\,\mathrm{physicalMpc}$) from other protocluster candidates of higher overdensity, and are expected to merge into a single massive structure by $z=0$. Therefore, we find 179 unique protocluster candidates in our survey. A cosmological simulation that includes projection effects predicts that more than 76\% of these candidates will evolve into galaxy clusters with halo masses of at least $10^{14}\,M_{\odot}$ by $z=0$. The unprecedented size of our protocluster candidate catalog allowed us to perform, for the first time, an angular clustering analysis of the systematic sample of protocluster candidates. We find a correlation length of $35.0\,h^{-1}\,\mathrm{Mpc}$. The relation between correlation length and number density of $z\sim3.8$ protocluster candidates is consistent with the prediction of the $\Lambda$CDM model, and the correlation length is similar to that of rich clusters in the local universe. This result suggests that our protocluster candidates are tracing similar spatial structures as those expected of the progenitors of rich clusters and enhances the confidence that our method to identify protoclusters at high redshifts is robust. In the coming years, our protocluster search will be extended to the entire HSC-SSP Wide sky coverage of $\sim1400\,\mathrm{deg^2}$ to probe cluster formation over a wide redshift range of $z\sim2\mathrm{-}6$.

We have developed an efficient Active Galactic Nucleus (AGN) selection method using 18-band Spectral Energy Distribution (SED) fitting in mid-infrared (mid-IR). AGNs are often obscured by gas and dust, and those obscured AGNs tend to be missed in optical, UV and soft X-ray observations. Mid-IR light can help us to recover them in an obscuration free way using their thermal emission. On the other hand, Star-Forming Galaxies (SFG) also have strong PAH emission features in mid-IR. Hence, establishing an accurate method to separate populations of AGN and SFG is important. However, in previous mid-IR surveys, only 3 or 4 filters were available, and thus the selection was limited. We combined AKARI's continuous 9 mid-IR bands with WISE and Spitzer data to create 18 mid-IR bands for AGN selection. Among 4682 galaxies in the AKARI NEP deep field, 1388 are selected to be AGN hosts, which implies an AGN fraction of 29.6$\pm$0.8$\%$ (among them 47$\%$ are Seyfert 1.8 and 2). Comparing the result from SED fitting into WISE and Spitzer colour-colour diagram reveals that Seyferts are often missed by previous studies. Our result has been tested by stacking median magnitude for each sample. Using X-ray data from Chandra, we compared the result of our SED fitting with WISE's colour box selection. We recovered more X-ray detected AGN than previous methods by 20$\%$.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

Understanding how QSO's UV radiation affects galaxy formation is vital to our understanding of reionization era. Using a custom made narrow-band filter, $NB906$, on Subaru/Suprime-Cam, we investigated the number density of Ly$\alpha$ emitters (LAE) around a QSO at z=6.4. To date, this is the highest redshift narrow-band observation, where LAEs around a luminous QSO are investigated. Due to the large field-of-view of Suprime-Cam, our survey area is $\sim$5400~cMpc$^2$, much larger than previously studies at z=5.7 ($\sim$200 cMpc$^2$). In this field, we previously found a factor of 7 overdensity of Lyman break galaxies (LBGs). Based on this, we expected to detect $\sim$100 LAEs down to $NB906$=25 ABmag. However, our 6.4 hour exposure found none. The obtained upper limit on the number density of LAEs is more than an order lower than the blank fields. Furthermore, this lower density of LAEs spans a large scale of 10 $p$Mpc across. A simple argument suggests a strong UV radiation from the QSO can suppress star-formation in halos with $M_{vir}<10^{10}M_{\odot}$ within a $p$Mpc from the QSO, but the deficit at the edge of the field (5 $p$Mpc) remains to be explained.

Mass, metallicity, and star formation rate (SFR) of a galaxy are crucial parameters in understanding galaxy formation and evolution. However, the relation among these is still a matter of debate for luminous infrared galaxies, which carry a bulk of SFR budget of the universe at $z\sim1$. We have investigated the relation among stellar mass, gas-phase oxygen abundance, and SFR of AKARI-detected mid-IR galaxies at $z\sim0.88$ in the AKARI NEP deep field. We observed about 350 AKARI sources with Subaru/FMOS NIR spectrograph, and detected secure and expected H$\alpha$ emission lines from 25 and 44 galaxies, respectively. The SFR of our sample is almost constant ($\sim 25M_{\odot}/yr$) over the stellar mass range of our sample. Compared with main-sequence (MS) galaxies at a similar redshift range, the average SFR of our detected sample is comparable for massive galaxies ($\sim10^{10.58}~M_{\odot}$), while higher by $\sim$0.6dex for less massive galaxies ($\sim 10^{10.05}~M_{\odot}$). We measure metallicities from the [NII]/H$\alpha$ emission line ratio. We find that the mass-metallicity relation of our individually measured sources agrees with that for optical-selected star-forming galaxies at $z\sim0.1$, while metallicities of stacked spectra agree with that of MS galaxies at $z\sim0.78$. Considering high SFR of individually measured sources, FMR of the IR galaxies is different from that at $z\sim0.1$. However, on the mass-metallicity plane, they are consistent with the MS galaxies, highlighting higher SFR of the IR galaxies. This suggests the evolutionary path of our IR galaxies is different from that of MS galaxies. A possible physical interpretation includes that the star-formation activities of IR galaxies at $z\sim0.88$ in our sample are enhanced by interaction and/or merger of galaxies, but the inflow of metal-poor gas is not yet induced, keeping the metallicity intact.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

Aims. We report our finding of a high excitation emission line nebula associated with an Ultra Luminous X-ray source (ULX) at $z=$ 0.027, which we found in our Chandra observation of the AKARI North Ecliptic Pole (NEP) Deep Field. Methods. We present a Chandra X-ray and Gran Telescopio Canarias (GTC) optical spectral analysis of the ULX blob. We investigate the nature of the emission line nebula by using line ratio diagnostic diagrams, and its physical properties estimated from Spectral Energy Distribution (SED) fitting. Results. The optical spectrum of this ULX blob shows emission line ratios that are located on the borderlines between star-formation and Seyfert regimes in [OIII]/H$\beta$-[OI]/H$\alpha$, [OIII]/H$\beta$-[SII]/H$\alpha$ and [OIII]/H$\beta$-[OIII]/[OII] diagnostic diagrams. These are in contrast with those of a nearby blob observed with the same slit, which occupy the HII regimes. This result suggests that the ionization of the emission line nebula associated with the ULX is significantly contributed by energy input from the accretion power of the ULX, in addition to the star formation activity in the blob, suggesting the existence of an accretion disk in the ULX emitting UV radiation, or exerting shock waves.

We investigate the star forming activity of a sample of infrared (IR)-bright dust-obscured galaxies (DOGs) that show an extreme red color in the optical and IR regime, $(i - [22])_{\rm AB} > 7.0$. Combining an IR-bright DOG sample with the flux at 22 $\mu$m $>$ 3.8 mJy discovered by Toba & Nagao (2016) with IRAS faint source catalog version 2 and AKARI far-IR (FIR) all-sky survey bright source catalog version 2, we selected 109 DOGs with FIR data. For a subsample of 7 IR-bright DOGs with spectroscopic redshift ($0.07 < z < 1.0$) that was obtained from literature, we estimated their IR luminosity, star formation rate (SFR), and stellar mass based on the spectral energy distribution fitting. We found that (i) WISE 22 $\mu$m luminosity at observed frame is a good indicator of IR luminosity for IR-bright DOGs and (ii) the contribution of active galactic nucleus (AGN) to IR luminosity increases with IR luminosity. By comparing the stellar mass and SFR relation for our DOG sample and literature, we found that most of IR-bright DOGs lie significantly above the main sequence of star-forming galaxies at similar redshift, indicating that the majority of IRAS- and/or AKARI-detected IR-bright DOGs are starburst galaxies.

The identification of high-redshift massive galaxies with old stellar populations may pose challenges to some models of galaxy formation. However, to securely classify a galaxy as quiescent, it is necessary to exclude significant ongoing star formation, something that can be challenging to achieve at high redshift. In this letter, we analyse deep ALMA/870um and SCUBA-2/450um imaging of the claimed "post-starburst" galaxy ZF-20115 at z=3.717 that exhibits a strong Balmer break and absorption lines. The far-infrared imaging reveals a luminous starburst located 0.4+/-0.1 arcsec (~3kpc in projection) from the position of the rest-frame ultra-violet/optical emission, with an obscured star-formation rate of 100 Mo/yr. This star-forming component is undetected in the rest-frame ultraviolet but contributes significantly to the lower angular resolution photometry at restframe wavelengths >3500A, significantly complicating the determination of a reliable stellar mass. Importantly, in the presence of dust obscuration, strong Balmer features are not a unique signature of a post-starburst galaxy and are indeed frequently observed in infrared-luminous galaxies. We conclude that the ZF20015 system does not pose a challenge to current models of galaxy formation and that deep sub-/millimeter observations are a prerequisite for any claims of quiescence. The multi-wavelength observations of ZF20115 unveil a complex system with an intricate and spatially-varying star-formation history. ZF20115 demonstrates that understanding high-redshift obscured starbursts will only be possible with multi-wavelength studies that include high-resolution observations, available with the JWST, at mid-infrared wavelengths.

The aim of this work is to characterize physical properties of Ultra Luminous Infrared Galaxies (ULIRGs) and Luminous Infrared Galaxies (LIRGs) detected in the far-infrared (FIR) 90um band in the AKARI Deep Field-South (ADF-S) survey. In particular, we want to estimate the AGN contribution to the [U]LIRGs' infrared emission and which types of AGNs are related to their activity. We examined 69 galaxies at z>0.05 detected at 90um by the AKARI satellite in the ADF-S, with optical counterparts and spectral coverage from the ultraviolet to the FIR. We used two independent spectral energy distribution fitting codes: one fitting the SED from FIR to FUV (CIGALE) and gray-body + power spectrum fit for the infrared part of the spectra (CMCIRSED) in order to identify a subsample of [U]LIRGs, and to estimate their properties. Based on the CIGALE SED fitting, we have found that [U]LIRGs selected at the 90um AKARI band compose ~56% of our sample (we found 17 ULIRGs and 22 LIRGs, spanning over the redshift range 0.06<z<1.23). Their physical parameters, such as stellar mass, star formation rate (SFR), and specific SFR are consistent with the ones found for other samples selected at IR wavelengths. We have detected a significant AGN contribution to the MIR luminosity for 63% of LIRGs and ULIRGs. Our LIRGs contain Type 1, Type 2, and intermediate types of AGN, whereas for ULIRGs, a majority (more than 50%) of AGN emission originates from Type 2 AGNs. The temperature--luminosity and temperature--mass relations for the dust component of ADF--S LIRGs and ULIRGs indicate that these relations are shaped by the dust mass and not by the increased dust heating. We conclude that LIRGs contain Type 1, Type 2, and intermediate types of AGNs, with an AGN contribution to the MIR emission at the median level of 13+/-3%, whereas the majority of ULIRGs contain Type 2 AGNs, with a median AGN fraction equal to 19+/-8%.

Cosmic rays in the energy range $10^{18.0}$ - $10^{18.5}$ eV are thought to have a light, probably protonic, composition. To study their origin one can search for anisotropy in their arrival directions. Extragalactic cosmic rays should be isotropic, but galactic cosmic rays of this type should be seen mostly along the galactic plane, and there should be a shortage of events coming from directions near the galactic anticenter. This is due to the fact that, under the influence of the galactic magnetic field, the transition from ballistic to diffusive behavior is well advanced, and this qualitative picture persists over the whole energy range. Guided by models of the galactic magnetic field that indicate that the enhancement along the galactic plane should have a standard deviation of about 20$^\circ$ in galactic latitude, and the deficit in the galactic anticenter direction should have a standard deviation of about 50$^\circ$ in galactic longitude, we use the data of the Telescope Array surface detector in $10^{18.0}$ to $10^{18.5}$ eV energy range to search for these effects. The data are isotropic. Neither an enhancement along the galactic plane nor a deficit in the galactic anticenter direction is found. Using these data we place an upper limit on the fraction of EeV cosmic rays of galactic origin at 1.3% at 95% confidence level.

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

We present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$\alpha$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$\alpha$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$\alpha$ emitters, which have H$\alpha$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$\alpha$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$\lambda$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$\alpha$, [NII]$\lambda \lambda$6548,6583, [SII]$\lambda \lambda$6717, 6731, and [OI]$\lambda \lambda$6300,6364 are seen.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We present mid-infrared (MIR) luminosity functions (LFs) of local star-forming (SF) galaxies in the AKARI NEP-Wide Survey field. In order to derive more accurate luminosity function, we used spectroscopic sample only. Based on the NEP-Wide point source catalogue containing a large number of infrared (IR) sources distributed over the wide (5.4 sq. deg.) field, we incorporated the spectroscopic redshift data for about 1790 selected targets obtained by optical follow-up surveys with MMT/Hectospec and WIYN/Hydra. The AKARI continuous 2 to 24 micron wavelength coverage as well as photometric data from optical u band to NIR H-band with the spectroscopic redshifts for our sample galaxies enable us to derive accurate spectral energy distributions (SEDs) in the mid-infrared. We carried out SED fit analysis and employed 1/Vmax method to derive the MIR (8, 12, and 15 micron rest-frame) luminosity functions. We fit our 8 micron LFs to the double power-law with the power index of alpha= 1.53 and beta= 2.85 at the break luminosity. We made extensive comparisons with various MIR LFs from several literatures. Our results for local galaxies from the NEP region are generally consistent with other works for different fields over wide luminosity ranges. The comparisons with the results from the NEP-Deep data as well as other LFs imply the luminosity evolution from higher redshifts towards the present epoch.

We present a method of selection of 24~$\mu$m galaxies from the AKARI North Ecliptic Pole (NEP) Deep Field down to $150 \mbox{ }\mu$Jy and measurements of their two-point correlation function. We aim to associate various 24 $\mu$m selected galaxy populations with present day galaxies and to investigate the impact of their environment on the direction of their subsequent evolution. We discuss using of Support Vector Machines (SVM) algorithm applied to infrared photometric data to perform star-galaxy separation, in which we achieve an accuracy higher than 80\%. The photometric redshift information, obtained through the CIGALE code, is used to explore the redshift dependence of the correlation function parameter ($r_{0}$) as well as the linear bias evolution. This parameter relates galaxy distribution to the one of the underlying dark matter. We connect the investigated sources to their potential local descendants through a simplified model of the clustering evolution without interactions. We observe two different populations of star-forming galaxies, at $z_{med}\sim 0.25$, $z_{med}\sim 0.9$. Measurements of total infrared luminosities ($L_{TIR}$) show that the sample at $z_{med}\sim 0.25$ is composed mostly of local star-forming galaxies, while the sample at $z_{med}\sim0.9$ is composed of luminous infrared galaxies (LIRGs) with $L_{TIR}\sim 10^{11.62}L_{\odot}$. We find that dark halo mass is not necessarily correlated with the $L_{TIR}$: for subsamples with $L_{TIR}= 10^{11.15} L_{\odot}$ at $z_{med}\sim 0.7$ we observe a higher clustering length ($r_{0}=6.21\pm0.78$ $[h^{-1} \mbox{Mpc}]$) than for a subsample with mean $L_{TIR}=10^{11.84} L_{\odot}$ at $z_{med}\sim1.1$ ($r_{0}=5.86\pm0.69$ $h^{-1} \mbox{Mpc}$). We find that galaxies at $z_{med}\sim 0.9$ can be ancestors of present day $L_{*}$ early type galaxies, which exhibit a very high $r_{0}\sim 8$~$h^{-1} \mbox{Mpc}$.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$\alpha$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$\lambda$6583/H$\alpha$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$\lambda\lambda$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$\lambda\lambda$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$\lambda$6583/H$\alpha$ vs. [S~{\sc ii}]$\lambda\lambda$6716,6731/H$\alpha$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.

The reionisation of the Universe is thought to have ended around z~6, as inferred from spectroscopy of distant bright background sources, such as quasars (QSO) and gamma-ray burst (GRB) afterglows. Furthermore, spectroscopy of a GRB afterglow provides insight in its host galaxy, which is often too dim and distant to study otherwise. We present the high S/N VLT/X-shooter spectrum of GRB130606A at z=5.913. We aim to measure the degree of ionisation of the IGM between 5.02<z<5.84 and to study the chemical abundance pattern and dust content of its host galaxy. We measured the flux decrement due to absorption at Ly$\alpha$, $\beta$ and $\gamma$ wavelength regions. The hydrogen and metal absorption lines formed in the host galaxy were fitted with Voigt profiles to obtain column densities. Our measurements of the Ly$\alpha$-forest optical depth are consistent with previous measurements of QSOs, but have a much smaller uncertainty. The analysis of the red damping wing yields a neutral fraction $x_{HI}<0.05$ (3$\sigma$). We obtain column density measurements of several elements. The ionisation corrections due to the GRB is estimated to be negligible (<0.03 dex), but larger corrections may apply due to the pre-existing radiation field (up to 0.4 dex based on sub-DLA studies). Our measurements confirm that the Universe is already predominantly ionised over the redshift range probed in this work, but was slightly more neutral at z>5.6. GRBs are useful probes of the ionisation state of the IGM in the early Universe, but because of internal scatter we need a larger statistical sample to draw robust conclusions. The high [Si/Fe] in the host can be due to dust depletion, alpha-element enhancement, or a combination of both. The very high value of [Al/Fe]=2.40+/-0.78 might connected to the stellar population history. We estimate the host metallicity to be -1.7<[M/H]<-0.9 (2%-13% of solar). (trunc.)

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$\mu$m, 12$\mu$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($\Omega_{\mathrm{IR}}$), we found that the $\Omega_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

Infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution, since their most intense stages are often obscured by dust. Japanese infrared satellite, AKARI, provided unique data sets to probe these both at low and high redshifts. The AKARI performed an all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160$\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, AKARI can measure the total infrared luminosity ($L_{TIR}$) of individual galaxies much more precisely, and thus, the total infrared luminosity density of the local Universe. In the AKARI NEP deep field, we construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2 using 4128 infrared sources. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from $z$=0 to $z$=2.2, all probed by the AKARI satellite. The next generation space infrared telescope, SPICA, will revolutionize our view of the infrared Universe with superb sensitivity of the cooled 3m space telescope. We conclude with our survey proposal and future prospects with SPICA.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

The extragalactic background suggests half the energy generated by stars reprocessed into the infrared (IR) by dust. At z$\sim$1.3, 90\% of star formation is obscured by dust. To fully understand the cosmic star formation history, it is critical to investigate infrared emission. AKARI has made deep mid-IR observation using its continuous 9-band filters in the NEP field (5.4 deg$^2$), using $\sim$10\% of the entire pointed observations available throughout its lifetime. However, there remain 11,000 AKARI's infrared sources undetected with the previous CFHT/Megacam imaging ($r\sim$25.9ABmag). Redshift and IR luminosity of these sources are unknown. These sources may contribute significantly to the cosmic star-formation rate density (CSFRD). For example, if they all lie at 1$<z<$2, the CSFRD will be twice as high at the epoch. We are carrying out deep imaging of the NEP field in 5 broad bands ($g,r,i,z,$ and $y$) using Hyper Suprime-Camera (HSC), which has 1.5 deg field of view in diameter on Subaru 8m telescope. This will provide photometric redshift information, and thereby IR luminosity for the previously-undetected 11,000 faint AKARI IR sources. Combined with AKARI's mid-IR AGN/SF diagnosis, and accurate mid-IR luminosity measurement, this will allow a complete census of cosmic star-formation/AGN accretion history obscured by dust.

Context. A possible correlation between CO luminosity (L_CO ) and its line width (FWHM) has been suggested and denied in the literature. Such claims were often based on a small, or heterogeneous sample of galaxies, and thus inconclusive. Aims. We aim to prove or dis-prove the L_CO -FWHM correlation. Methods. We compile a large sample of submm galaxies at z>2 from the literature, and investigate the L_CO-FWHM relation. Results. After carefully evaluating the selection effects and uncertainties such as inclination and magnification via gravitational lensing, we show that there exist a weak but significant correlation between L_CO and FWHM. We also discuss a feasibility to measure the cosmological distance using the correlation.

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution, since their most intense stages are often obscured by dust. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshifts. The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160$\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, AKARI can much more precisely measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe. In the AKARI NEP deep field, we construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2 using 4128 infrared sources. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from $z$=0 to $z$=2.2, all probed by the AKARI satellite.

(Abridged) We aim to study the evolution of dust attenuation in galaxies selected in the IR in the redshift range in which they are known to dominate the star formation activity in the universe. The comparison with other measurements of dust attenuation in samples selected using different criteria will give us a global picture of the attenuation at work in star-forming galaxies and its evolution with redshift. Using multiple filters of IRC instrument, we selected more than 4000 galaxies from their rest-frame emission at 8 microns, from z~0.2 to 2$. We built SEDs from the rest-frame UV to the far-IR by adding data in the optical-NIR and from GALEX and Herschel surveys. We fit SEDs with the physically-motivated code CIGALE. We test different templates for AGNs and recipes for dust attenuation and estimate stellar masses, SFRs, amount of dust attenuation, and AGN contribution to the total IR luminosity. The AGN contribution to the total IR luminosity is found to be on average approximately 10% with a slight increase with redshift. Dust attenuation in galaxies dominating the IR luminosity function is found to increase from z=0 to z=1 and to remain almost constant from z=1 to z=1.5. Conversely, when galaxies are selected at a fixed IR luminosity, their dust attenuation slightly decreases as redshift increases but with a large dispersion. The attenuation in our mid-IR selected sample is found ~ 2 mag higher than that found globally in the universe or in UV and Halpha line selections in the same redshift range. This difference is well explained by an increase of dust attenuation with the stellar mass, in global agreement with other recent studies. Starbursting galaxies do not systematically exhibit a high attenuation

We present the rest-frame optical spectral properties of 155 luminous quasars at 3.3<z<6.4 taken with the AKARI space telescope, including the first detection of H$\alpha$ emission line as far out as z~6. We extend the scaling relation between the rest-frame optical continuum and line luminosity of active galactic nuclei (AGNs) to the high luminosity, high redshift regime that has rarely been probed before. Remarkably, we find that a single log-linear relation can be applied to the 5100${\rm \AA}$ and H$\alpha$ AGN luminosities over a wide range of luminosity (10$^{42}$<$L_{5100}$<10$^{47}$ergs/s) or redshift (0<z<6), suggesting that the physical mechanism governing this relation is unchanged from z=0 to 6, over five decades in luminosity. Similar scaling relations are found between the optical and the UV continuum luminosities or line widths. Applying the scaling relations to the H$\beta$ black hole mass ($M_{\rm BH}$) estimator of local AGNs, we derive the $M_{\rm BH}$ estimators based on H$\alpha$, MgII, and CIV lines, finding that the UV-line based masses are overall consistent with the Balmer-line based, but with a large intrinsic scatter of 0.40dex for the CIV estimates. Our 43 $M_{\rm BH}$ estimates from H$\alpha$ confirm the existence of BHs as massive as ~10$^{10}M_{\odot}$ out to z~5, and provide a secure footing for previous MgII-line based studies that a rapid $M_{\rm BH}$ growth has occurred in the early universe.

We present data products from the 300 ks Chandra survey in the AKARI North Ecliptic Pole (NEP) deep field. This field has a unique set of 9-band infrared photometry covering 2-24 micron from the AKARI Infrared Camera, including mid-infrared (MIR) bands not covered by Spitzer. The survey is one of the deepest ever achieved at ~15 micron, and is by far the widest among those with similar depths in the MIR. This makes this field unique for the MIR-selection of AGN at z~1. We design a source detection procedure, which performs joint Maximum Likelihood PSF fits on all of our 15 mosaicked Chandra pointings covering an area of 0.34 square degree. The procedure has been highly optimized and tested by simulations. We provide a point source catalog with photometry and Bayesian-based 90 per cent confidence upper limits in the 0.5-7, 0.5-2, 2-7, 2-4, and 4-7 keV bands. The catalog contains 457 X-ray sources and the spurious fraction is estimated to be ~1.7 per cent. Sensitivity and 90 per cent confidence upper flux limits maps in all bands are provided as well. We search for optical MIR counterparts in the central 0.25 square degree, where deep Subaru Suprime-Cam multiband images exist. Among the 377 X-ray sources detected there, ~80 per cent have optical counterparts and ~60 per cent also have AKARI mid-IR counterparts. We cross-match our X-ray sources with MIR-selected AGN from Hanami et al. (2012). Around 30 per cent of all AGN that have MID-IR SEDs purely explainable by AGN activity are strong Compton-thick AGN candidates.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.

Previous measurements of the composition of Ultra-High Energy Cosmic Rays(UHECRs) made by the High Resolution Fly's Eye(HiRes) and Pierre Auger Observatory(PAO) are seemingly contradictory, but utilize different detection methods, as HiRes was a stereo detector and PAO is a hybrid detector. The five year Telescope Array(TA) Middle Drum hybrid composition measurement is similar in some, but not all, respects in methodology to PAO, and good agreement is evident between data and a light, largely protonic, composition when comparing the measurements to predictions obtained with the QGSJetII-03 and QGSJet-01c models. These models are also in agreement with previous HiRes stereo measurements, confirming the equivalence of the stereo and hybrid methods. The data is incompatible with a pure iron composition, for all models examined, over the available range of energies. The elongation rate and mean values of Xmax are in good agreement with Pierre Auger Observatory data. This analysis is presented using two methods: data cuts using simple geometrical variables and a new pattern recognition technique.

We present a new catalog of 118 Ultraluminous Infrared Galaxies (ULIRGs) and one Hyperluminous Infrared Galaxy (HLIRG) by crossmatching AKARI all-sky survey with the Sloan Digital Sky Survey Data Release 10 (SDSS DR10) and the Final Data Release of the Two-Degree Field Galaxy Redshift Survey (2dFGRS). 40 of the ULIRGs and one HLIRG are new identifications. We find that ULIRGs are interacting pair galaxies or ongoing/post mergers. This is consistent with the widely accepted view: ULIRGs are major mergers of disk galaxies. We confirm the previously known positive trend between the AGN fraction and IR luminosity. We show that ULIRGs have a large off-set from the 'main sequence' up to z~1; their off-set from the z~2 'main sequence' is relatively smaller. We find a consistent result with the previous studies showing that compared to local star forming SDSS galaxies of similar mass, local ULIRGs have lower oxygen abundances. We for the first time demonstrate that ULIRGs follow the fundamental metallicity relation (FMR). The scatter of ULIRGs around the FMR (0.09 dex - 0.5 dex) is comparable with the scatter of z~2-3 galaxies. Their optical colors show that ULIRGs are mostly blue galaxies and this agrees with previous findings. We provide the largest local (0.050 < z < 0.487) ULIRG catalog with stellar masses, SFRs, gas metallicities and optical colors. Our catalog provides us active galaxies analogous to high-z galaxies in the local Universe where they can be rigorously scrutinized.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We present herein galaxy number counts of the nine bands in the 2-24 micron range on the basis of the AKARI North Ecliptic Pole (NEP) surveys. The number counts are derived from NEP-deep and NEP-wide surveys, which cover areas of 0.5 and 5.8 deg2, respectively. To produce reliable number counts, the sources were extracted from recently updated images. Completeness and difference between observed and intrinsic magnitudes were corrected by Monte Carlo simulation. Stellar counts were subtracted by using the stellar fraction estimated from optical data. The resultant source counts are given down to the 80% completeness limit; 0.18, 0.16, 0.10, 0.05, 0.06, 0.10, 0.15, 0.16, and 0.44 mJy in the 2.4, 3.2, 4.1, 7, 9, 11, 15, 18 and 24 um bands, respectively. On the bright side of all bands, the count distribution is flat, consistent with the Euclidean Universe, while on the faint side, the counts deviate, suggesting that the galaxy population of the distant universe is evolving. These results are generally consistent with previous galaxy counts in similar wavebands. We also compare our counts with evolutionary models and find them in good agreements. By integrating the models down to the 80% completeness limits, we calculate that the AKARI NEP-survey revolves 20%-50% of the cosmic infrared background, depending on the wavebands.

We present the first galaxy counts at 18 microns using the Japanese AKARI satellite's survey at the North Ecliptic Pole (NEP), produced from the images from the NEP-Deep and NEP-Wide surveys covering 0.6 and 5.8 square degrees respectively. We describe a procedure using a point source filtering algorithm to remove background structure and a minimum variance method for our source extraction and photometry that delivers the optimum signal to noise for our extracted sources, confirming this by comparison with standard photometry methods. The final source counts are complete and reliable over three orders of magnitude in flux density, resulting in sensitivities (80 percent completeness) of 0.15mJy and 0.3mJy for the NEP-Deep and NEP-Wide surveys respectively, a factor of 1.3 deeper than previous catalogues constructed from this field. The differential source counts exhibit a characteristic upturn from Euclidean expectations at around a milliJansky and a corresponding evolutionary bump between 0.2-0.4 mJy consistent with previous mid-infrared surveys with ISO and Spitzer at 15 and 24 microns. We compare our results with galaxy evolution models confirming the striking divergence from the non-evolving scenario. The models and observations are in broad agreement implying that the source counts are consistent with a strongly evolving population of luminous infrared galaxies at redshifts higher than unity. Integrating our source counts down to the limit of the NEP survey at the 150 microJy level we calculate that AKARI has resolved approximately 55 percent of the 18 micron cosmic infrared background relative to the predictions of contemporary source count models.

We have searched for intermediate-scale anisotropy in the arrival directions of ultrahigh-energy cosmic rays with energies above 57~EeV in the northern sky using data collected over a 5 year period by the surface detector of the Telescope Array experiment. We report on a cluster of events that we call the hotspot, found by oversampling using 20$^\circ$-radius circles. The hotspot has a Li-Ma statistical significance of 5.1$\sigma$, and is centered at R.A.=146.7$^{\circ}$, Dec.=43.2$^{\circ}$. The position of the hotspot is about 19$^{\circ}$ off of the supergalactic plane. The probability of a cluster of events of 5.1$\sigma$ significance, appearing by chance in an isotropic cosmic-ray sky, is estimated to be 3.7$\times$10$^{-4}$ (3.4$\sigma$).

We study the behaviour of polycyclic aromatic hydrocarbon emission in galaxies at z=0.3-1.4 using 1868 samples from the revised catalogue of AKARI North Ecliptic Pole Deep survey. The continuous filter coverage at 2-24um makes it possible to measure 8um luminosity, which is dominated by polycyclic aromatic hydrocarbon emission for galaxies at up to z=2. We compare the IR8 (= LIR/L(8)) and 8um to 4.5um luminosity ratio (L(8)/L(4.5)) with the starburstiness, Rsb, defined as excess of specific star -formation rate over that of main-sequence galaxy. All AGN candidates were excluded from our sample using an SED fitting. We found L(8)/L(4.5) increases with starburstiness at log Rsb < 0.5 and stays constant at higher starburstiness. On the other hand, IR8 is constant at log Rsb < 0, while it increases with starburstiness at log Rsb > 0. This behaviour is seen in all redshift range of our study. These results indicate that starburst galaxies have deficient polycyclic aromatic hydrocarbon emission compared with main-sequence galaxies. We also find that galaxies with extremely high L(8)/L(4.5) ratio have only moderate starburstiness. These results suggest that starburst galaxies have compact star-forming regions with intense radiation, which destroys PAHs and/or have dusty HII regions resulting in a lack of ionising photons.

Aims. We present an 8-band (u*, g', r', i', z', Y, J, Ks) optical to NIR deep photometric catalog based on the observations made with MegaCam and WIRCam at CFHT, and compute photometric redshifts, zp in the North Ecliptic Pole (NEP) region. Our catalog provides us to identify the counterparts, and zp for AKARI NIR/MIR sources. Results. The estimated 4sigma detection limits within an 1" aperture radius are 26.7, 25.9, 25.1, and 24.1 mag [AB] for g', r', i', and z'-bands and 23.4, 23.0, and 22.7 mag for Y, J, and Ks-bands, respectively. There are a total of 85797 sources in the band-merged catalog. An astrometric accuracy of this catalog determined by examining coordinate offsets with regard to 2MASS is 0.013" with a root mean square offset of 0.32". We distinguish 5441 secure stars from extended sources using u*-J vs. g'-Ks colors, combined with the SExtractor stellarity index of the images. Comparing with galaxy spectroscopic redshifts, we find a photometric redshift dispersion, sigma_(dz/(1+z)), of 0.032 and catastrophic failure rate, dz/(1+z)>0.15, of 5.8% at z<1, while a dispersion of 0.117 and a catastrophic failures rate of 16.6% at z>1. We extend estimate of the zp uncertainty over the full magnitude/redshift space with a redshift probability distribution function and find that our redshift are highly accurate with z'<22 at zp<2.5 and for fainter sources with z'<24 at z<1. From the investigation of photometric properties of AKARI infrared sources (23354 sources) using the g'z'Ks diagram, <5% of AKARI sources with optical counterparts are classified as high-z (1.4<z<2.5) star-forming galaxies. Among the high-z star-forming galaxies, AKARI MIR detected sources seem to be affected by stronger dust extinction compared with sources with non-detections in the AKARI MIR bands. The full, electronic version of our catalog with zp will be available at the CDS.

The efficient selection of high-redshift emission galaxies is important for future large galaxy redshift surveys for cosmology. Here we describe the target selection methods for the FastSound project, a redshift survey for H alpha emitting galaxies at z=1.2-1.5 using Subaru/FMOS to measure the linear growth rate f\sigma 8 via Redshift Space Distortion (RSD) and constrain the theory of gravity. To select ~400 target galaxies in the 0.2 deg^2 FMOS field-of-view from photometric data of CFHTLS-Wide (u*g'r'i'z'), we test several different methods based on color-color diagrams or photometric redshift estimates from spectral energy distribution (SED) fitting. We also test the improvement in selection efficiency that can be achieved by adding near-infrared data from the UKIDSS DXS (J). The success rates of H alpha detection with FMOS averaged over two observed fields using these methods are 11.3% (color-color, optical), 13.6% (color-color, optical+NIR), 17.3% (photo-z, optical), and 15.1% (photo-z, optical+NIR). Selection from photometric redshifts tends to give a better efficiency than color-based methods, although there is no significant improvement by adding J band data within the statistical scatter. We also investigate the main limiting factors for the success rate, by using the sample of the HiZELS H alpha emitters that were selected by narrow-band imaging. Although the number density of total H alpha emitters having higher H alpha fluxes than the FMOS sensitivity is comparable with the FMOS fiber density, the limited accuracy of photometric redshift and H alpha flux estimations have comparable effects on the success rate of <~20% obtained from SED fitting.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.

We derive the stacked 1.4 GHz flux from FIRST (Faint Images of the Radio Sky at Twenty Centimeters) survey for 811 K+A galaxies selected from the SDSS DR7. For these objects we find a mean flux density of $56\pm 9$ $\mu$Jy. A similar stack of radio-quiet white dwarfs yields an upper limit of 43 $\mu$Jy at a 5$\sigma$ significance to the flux in blank regions of the sky. This implies an average star formation rate of 1.6 $\pm$ 0.3 M$_{\odot}$ year$^{-1}$ for K+A galaxies. However the majority of the signal comes from $\sim$4% of K+A fields that have aperture fluxes above the $5\sigma$ noise level of the FIRST survey. A stack of the remaining galaxies shows little residual flux consistent with an upper limit on star formation of 1.3 M$_{\odot}$ year$^{-1}$. Even for a subset of 456 `young' (spectral ages $<$ 250 Myr) K+A galaxies we find that the stacked 1.4 GHz flux is consistent with no current star formation. Our data suggest that the original starburst has been terminated in the majority of K+A galaxies, but that this may represent part of a duty cycle where a fraction of these galaxies may be active at a given moment with dusty starbursts and AGNs being present.

We present a photometric catalog of infrared (IR) sources based on the North Ecliptic PoleWide field (NEP-Wide) survey of AKARI, which is an infrared space telescope launched by Japan. The NEP-Wide survey covered 5.4 deg2 area, a nearly circular shape centered on the North Ecliptic Pole, using nine photometric filter-bands from 2 - 25 {\mu}m of the Infrared Camera (IRC). Extensive efforts were made to reduce possible false objects due to cosmic ray hits, multiplexer bleeding phenomena around bright sources, and other artifacts. The number of detected sources varied depending on the filter band: with about 109,000 sources being cataloged in the near-IR bands at 2 - 5 {\mu}m, about 20,000 sources in the shorter parts of the mid-IR bands between 7 - 11 {\mu}m, and about 16,000 sources in the longer parts of the mid-IR bands, with \sim 4,000 sources at 24 {\mu}m. The estimated 5? detection limits are approximately 21 magnitude (mag) in the 2 - 5 {\mu}m bands, 19.5 - 19 mag in the 7 - 11 {\mu}m, and 18.8 - 18.5 mag in the 15 - 24 {\mu}m bands in the AB magnitude scale. The completenesses for those bands were evaluated as a function of magnitude: the 50% completeness limits are about 19.8 mag at 3 {\mu}m, 18.6 mag at 9 {\mu}m, and 18 mag at 18 {\mu}m band, respectively. To construct a reliable source catalog, all of the detected sources were examined by matching them with those in other wavelength data, including optical and ground-based near-IR bands. The final band-merged catalog contains about 114,800 sources detected in the IRC filter bands. The properties of the sources are presented in terms of the distributions in various color-color diagrams.

The results of a deep 20 cm radio survey at 20 cm are reported of the AKARI Deep Field South (ADF-S) near the South Ecliptic Pole (SEP), using the Australia Telescope Compact Array telescope, ATCA. The survey has 1 sigma detection limits ranging from 18.7--50 microJy per beam over an area of ~1.1 sq degrees, and ~2.5 sq degrees to lower sensitivity. The observations, data reduction and source count analysis are presented, along with a description of the overall scientific objectives, and a catalogue containing 530 radio sources detected with a resolution of 6.2" x 4.9". The derived differential source counts show a pronounced excess of sources fainter than ~1 mJy, consistent with an emerging population of star forming galaxies. Cross-correlating the radio with AKARI sources and archival data we find 95 cross matches, with most galaxies having optical R-magnitudes in the range 18-24 mag, and 52 components lying within 1" of a radio position in at least one further catalogue (either IR or optical). We have reported redshifts for a sub-sample of our catalogue finding that they vary between galaxies in the local universe to those having redshifts of up to 0.825. Associating the radio sources with the Spitzer catalogue at 24 microns, we find 173 matches within one Spitzer pixel, of which a small sample of the identifications are clearly radio loud compared to the bulk of the galaxies. The radio luminosity plot and a colour-colour analysis suggest that the majority of the radio sources are in fact luminous star forming galaxies, rather than radio-loud AGN. There are additionally five cross matches between ASTE or BLAST submillimetre galaxies and radio sources from this survey, two of which are also detected at 90 microns, and 41 cross-matches with submillimetre sources detected in the Herschel HerMES survey Public Data release.