Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$\sigma$. The ten years of data post-trial significance is 4.1$\sigma$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented.

Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger.

Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.

We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function.

The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper.

The Wide-Field InfraRed Survey Telescope (WFIRST) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point-spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to the cosmological landscape after the initial years of observations from other concurrent Stage IV dark energy experiments. This potential requires equally careful planning and requirements validation as the mission prepares to enter its construction phase. We present a suite of image simulations based on GalSim that are used to construct a complex, synthetic WFIRST weak lensing survey that incorporates realistic input galaxies and stars, relevant detector non-idealities, and the current reference five-year WFIRST survey strategy. We present a first study to empirically validate the existing WFIRST weak lensing requirements flowdown using a suite of 12 matched image simulations, each representing a different perturbation to the wavefront or image motion model. These are chosen to induce a range of potential static and low- and high-frequency time-dependent PSF model errors. We analyze the measured shapes of galaxies from each of these simulations and compare them to a reference, fiducial simulation to infer the response of the shape measurement to each of these modes in the wavefront model. We then compare this to existing analytic flowdown requirements, and find general agreement between the empirically derived response and that predicted by the analytic model.

We present the large-scale structure over more than 50 comoving Mpc scale at z $\sim$ 0.9 where the CL1604 supercluster, which is one of the largest structures ever known at high redshifts, is embedded. The wide-field deep imaging survey by the Subaru Strategic Program with Hyper Suprime-Cam reveals that the already-known CL1604 supercluster is a mere part of larger-scale structure extending to both the north and the south. We confirm that there are galaxy clusters at three slightly different redshifts in the northern and southern sides of the supercluster by determining the redshifts of 55 red-sequence galaxies and 82 star-forming galaxies in total by the follow-up spectroscopy with Subaru/FOCAS and Gemini-N/GMOS. This suggests that the structure ever known as the CL1604 supercluster is the tip of the iceberg. We investigate stellar population of the red-sequence galaxies using 4000 A break and Balmer H$\delta$ absorption line. Almost all of the red-sequence galaxies brighter than 21.5 mag in $z$-band show an old stellar population with $\gtrsim2$ Gyr. The comparison of composite spectra of the red-sequence galaxies in the individual clusters show that the galaxies at a similar redshift have similar stellar population age, even if they are located $\sim$50 Mpc apart from each other. However, there could be a large variation in the star formation history. Therefore, it is likely that galaxies associated with the large-scale structure at 50 Mpc scale formed at almost the same time, have assembled into the denser regions, and then have evolved with different star formation history along the hierarchical growth of the cosmic web.

We investigate the locations of satellite galaxies in the z = 0 redshift slice of the hydrodynamical Illustris-1 simulation. As expected from previous work, the satellites are distributed anisotropically in the plane of the sky, with a preference for being located near the major axes of their hosts. Due to misalignment of mass and light within the hosts, the degree of anisotropy is considerably less when satellite locations are measured with respect to the hosts' stellar surface mass density than when they are measured with respect to the hosts' dark matter surface mass density. When measured with respect to the hosts' dark matter surface mass density, the mean satellite location depends strongly on host stellar mass and luminosity, with the satellites of the faintest, least massive hosts showing the greatest anisotropy. When measured with respect to the hosts' stellar surface mass density, the mean satellite location is essentially independent of host stellar mass and luminosity. In addition, the satellite locations are largely insensitive to the amount of stellar mass used to define the hosts' stellar surface mass density, as long as at least 50% to 70% of the hosts' total stellar mass is used. The satellite locations are dependent upon the stellar masses of the satellites, with the most massive satellites having the most anisotropic distributions.

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.

Improved statistics and mass-composition-sensitive observation are required to clarify the origin of ultra-high energy cosmic rays (UHECRs). Inevitably in the future, the UHECR observatories will have to be expanded due to the small flux; however, the upgrade will be expensive with the detectors currently in use. Hence, we are developing a new fluorescence detector for UHECR observation. The proposed fluorescence detector, called cosmic ray air fluorescence Fresnel-lens telescope (CRAFFT), has an extremely simple structure and can observe the longitudinal development of an air shower. Furthermore, CRAFFT has the potential to significantly reduce costs for the realization of a huge observatory for UHECR research. We deployed four CRAFFT detectors at the Telescope Array site and conducted the test observation. We have successfully observed ten air-shower events using CRAFFT. Thus, CRAFFT can be a solution to realize the next generation of UHECR observatories.

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.

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.

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.

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}$).

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.

We present new results from near-infrared spectroscopy with Keck/MOSFIRE of [OIII]-selected galaxies at $z\sim3.2$. With our $H$ and $K$-band spectra, we investigate the interstellar medium (ISM) conditions, such as ionization states and gas metallicities. [OIII] emitters at $z\sim3.2$ show a typical gas metallicity of $\mathrm{12+log(O/H) = 8.07\pm0.07}$ at $\mathrm{log(M_*/M_\odot) \sim 9.0-9.2}$ and $\mathrm{12+log(O/H) = 8.31\pm0.04}$ at $\mathrm{log(M_*/M_\odot) \sim 9.7-10.2}$ when using the empirical calibration method. We compare the [OIII] emitters at $z\sim3.2$ with UV-selected galaxies and Ly$\alpha$ emitters at the same epoch and find that the [OIII]-based selection does not appear to show any systematic bias in the selection of star-forming galaxies. Moreover, comparing with star-forming galaxies at $z\sim2$ from literature, our samples show similar ionization parameters and gas metallicities as those obtained by the previous studies using the same calibration method. We find no strong redshift evolution in the ISM conditions between $z\sim3.2$ and $z\sim2$. Considering that the star formation rates at a fixed stellar mass also do not significantly change between the two epochs, our results support the idea that the stellar mass is the primary quantity to describe the evolutionary stages of individual galaxies at $z>2$.

We present the molecular gas mass fraction ($f_\mathrm{H_2}$) and star-formation efficiency (SFE) of local galaxies on the basis of our new CO($J=1-0$) observations with the Nobeyama 45m radio telescope, combined with the COLDGASS galaxy catalog, as a function of galaxy environment defined as the local number density of galaxies measured with SDSS DR7 spectroscopic data. Our sample covers a wide range in the stellar mass and SFR, and covers wide environmental range over two orders of magnitude. This allows us to conduct the first, systematic study of environmental dependence of molecular gas properties in galaxies from the lowest- to the highest-density environments in the local universe. We confirm that both $f_\mathrm{H_2}$ and SFE have strong positive correlations with the SFR offset from the star-forming main sequence ($\Delta$MS), and most importantly, we find that these correlations are universal across all environments. Our result demonstrates that star-formation activity within individual galaxies is primarily controlled by their molecular gas content, regardless of their global environment. Therefore, we claim that one always needs to be careful about the $\Delta$MS distribution of the sample when investigating the environmental effects on the H$_2$ gas content in galaxies.

We present the environmental dependence of colour, stellar mass, and star formation (SF) activity in H-alpha-selected galaxies along the large-scale structure at z=0.4 hosting twin clusters in DEEP2-3 field, discovered by Subaru Strategic Programme of Hyper Suprime-Cam (HSC SSP). By combining photo-z selected galaxies and H-alpha emitters selected with broad-band and narrow-band (NB) data of the recent data release of HSC SSP (DR1), we confirm that galaxies in higher-density environments or galaxies in the cluster central regions show redder colours. We find that there still remains a possible colour-density and colour-radius correlation even if we restrict the sample to H-alpha-selected galaxies, likely due to the presence of massive H-alpha emitters in denser regions. We also find a hint of increased star formation rates (SFR) amongst H-alpha emitters towards the highest-density environment, again primarily driven by the excess of red/massive H-alpha emitters in high-density environment, while their specific SFR does not significantly change with environment. This work demonstrates the power of the HSC SSP NB data to study SF galaxies across environment in the distant universe.

We have carried out a panoramic Ly$\alpha$ narrowband imaging with Suprime-Cam on Subaru towards the known protocluster USS1558--003 at $z=2.53$. Our previous narrowband imaging at near-infrared has identified multiple dense groups of H$\alpha$ emitters (HAEs) within the protocluster. We have now identified the large-scale structures across a $\sim$50 comoving Mpc scale traced by Ly$\alpha$ emitters (LAEs) in which the protocluster traced by the HAEs is embedded. On a smaller scale, however, there are remarkably few LAEs in the regions of HAE overdensities. Moreover, the stacking analyses of the images show that HAEs in higher-density regions show systematically lower escape fractions of Ly$\alpha$ photons than those of HAEs in lower-density regions. These phenomena may be driven by the extra depletion of Ly$\alpha$ emission lines along our line of sight by more intervening cold circumgalactic/intergalactic medium and/or dust existing in the dense core. We also caution that all the past high-$z$ protocluster surveys using LAEs as the tracers would have largely missed galaxies in the very dense cores of the protoclusters where we would expect to see any early environmental effects.

We conducted a deep narrow-band imaging survey with the Subaru Prime Focus Camera on the Subaru Telescope and constructed a sample of Ly$\alpha$ emitters (LAEs) at z=2.53 in the UDS-CANDELS field where a sample of H$\alpha$ emitters (HAEs) at the same redshift is already obtained from our previous narrow-band observation at NIR. The deep narrow-band and multi broadband data allow us to find LAEs of stellar masses and star-formation rates (SFRs) down to $\gtrsim$$10^8$ M$_\odot$ and $\gtrsim$0.2 M$_\odot$/yr, respectively. We show that the LAEs are located along the same mass-SFR sequence traced by normal star-forming galaxies such as HAEs, but towards a significantly lower mass regime. Likewise, LAEs seem to share the same mass--size relation with typical star-forming galaxies, except for the massive LAEs, which tend to show significantly compact sizes. We identify a vigorous mass growth in the central part of LAEs: the stellar mass density in the central region of LAEs increases as their total galaxy mass grows. On the other hand, we see no Ly$\alpha$ line in emission for most of the HAEs. Rather, we find that the Ly$\alpha$ feature is either absent or in absorption (Ly$\alpha$ absorbers; LAAs), and its absorption strength may increase with reddening of the UV continuum slope. We demonstrate that a deep Ly$\alpha$ narrow-band imaging like this study is able to search for not only LAEs but also LAAs in a certain redshift slice. This work suggests that LAEs trace normal star-forming galaxies in the low-mass regime, while they remain as a unique population because the majority of HAEs are not LAEs.

We investigate a correlation between star-formation rate (SFR) and stellar mass for Halpha emission line galaxies (HAEs) in one of the richest proto-clusters ever known at z~2.5, USS 1558-003 proto-cluster. This study is based on a 9.7-hour narrow-band imaging data with MOIRCS on the Subaru telescope. We are able to construct a sample, in combination with additional H-band data taken with WFC3 on Hubble Space Telescope (HST), of 100 HAEs reaching the dust-corrected SFRs down to 3 Msun/yr and the stellar masses down to $10^{8.0}$ Msun. We find that while the star-forming galaxies with >$10^{9.3}$ Msun are located on the universal SFR-mass main sequence irrespective of the environment, less massive star-forming galaxies with <$10^{9.3}$ Msun show a significant upward scatter from the main sequence in this proto-cluster. This suggests that some less massive galaxies are in a starburst phase, although we do not know yet if this is due to environmental effects.

We investigate the properties of $z$=2.23 H$\alpha$ and [OIII]$\lambda$5007 emitters using the narrow-band-selected samples obtained from the High-$z$ Emission Line Survey (HiZELS: Sobral et al. 2013). We construct two samples of the H$\alpha$ and [OIII] emitters and compare their integrated physical properties. We find that the distribution of stellar masses, dust extinction, star formation rates (SFRs), and specific SFRs, is not statistically different between the two samples. When we separate the full galaxy sample into three subsamples according to the detections of the H$\alpha$ and/or [OIII] emission lines, most of the sources detected with both H$\alpha$ and [OIII] show ${\rm log(sSFR_{UV})}$$\gtrsim$-9.5. The comparison of the three subsamples suggests that sources with strong [OIII] line emission tend to have the highest star-forming activity out all galaxies that we study. We argue that the [OIII] emission line can be used as a tracer of star-forming galaxies at high redshift, and that it is especially useful to investigate star-forming galaxies at $z$$>$3, for which H$\alpha$ emission is no longer observable from the ground.

We present the physical properties of [OIII] emission line galaxies at z>3 as the tracers of active galaxies at 1Gyr before the peak epoch at z~2. We have performed deep narrow-band imaging surveys in the Subaru/XMM-Newton Deep Survey Field with MOIRCS on the Subaru Telescope and have constructed coherent samples of 34 [OIII] emitters at z=3.2 and 3.6, as well as 107 H$\alpha$ emitters at z=2.2 and 2.5. We investigate their basic physical quantities, such as stellar masses, star formation rates (SFRs), and sizes using the publicly available multi-wavelength data and high resolution images by the Hubble Space Telescope. The stellar masses and SFRs show a clear correlation known as the "main sequence" of star-forming galaxies. It is found that the location of the main sequence of the [OIII] emitters at z=3.2 and 3.6 is almost identical to that of the H$\alpha$ emitters at z=2.2 and 2.5. Also, we investigate their mass-size relation and find that the relation does not change between the two epochs. When we assume that the star-forming galaxies at z=3.2 grow simply along the same main sequence down to z=2.2, galaxies with $M_* = 10^{9}$-$10^{11} M_{\odot}$ increase their stellar masses significantly by a factor of 10-2. They climb up the main sequence, and their star formation rates also increase a lot as their stellar masses grow. This indicates that star formation activities of galaxies are accelerated from z>3 towards the peak epoch of galaxy formation at z~2.

In the redshift interval of $2<z<3$, the physical conditions of the inter-stellar medium (ISM) in star-forming galaxies are likely to be different from those in the local Universe because of lower gaseous metallicities, higher gas fractions, and higher star formation activities. In fact, observations suggest that higher electron densities, higher ionization parameters, and harder UV radiation fields are common. In this paper, based on the spectra of H$\alpha$-selected star-forming galaxies at $z=2.5$ taken with Multi-Object Spectrometer for InfraRed Exploration (MOSFIRE) on Keck-1 telescope, we measure electron densities ($n_e$) using the oxygen line ratio ( [OII]$\lambda\lambda$3726,3729), and investigate the relationships between the electron density of ionized gas and other physical properties. As a result, we find that the specific star formation rate (sSFR) and the surface density of SFR ($\Sigma_\mathrm{SFR}$) are correlated with the electron density at $z=2.5$ for the first time. The $\Sigma_\mathrm{SFR}-n_e$ relation is likely to be linked to the star formation law in HII regions (where star formation activity is regulated by interstellar pressure). Moreover, we discuss the mode of star formation in those galaxies. The correlation between sSFR and $\Sigma_\mathrm{SFR}$ suggests that highly star-forming galaxies (with high sSFR) tend to be characterized by higher surface densities of star formation ($\Sigma_\mathrm{SFR}$) and thus higher $n_e$ values as well.

We present a first application of the subhalo abundance matching (SHAM) method to describe the redshift-space clustering of galaxies including the non-linear redshift-space distortion, i.e., the Fingers-of-God. We find that the standard SHAM connecting the luminosity of galaxies to the maximum circular velocity of subhalos well reproduces the luminosity dependence of redshift-space clustering of galaxies in the Sloan Digital Sky Survey in a wide range of scales from 0.3 to 40 Mpc/h. The result indicates that the SHAM approach is very promising for establishing a theoretical model of redshift-space galaxy clustering without additional parameters. We also test color abundance matching using two different proxies for colors: subhalo age and local dark matter density following the method by Masaki et al. (2013b). Observed clustering of red galaxies exhibits much stronger Fingers-of-God effect than blue galaxies. We find that the subhalo age model describes the observed color-dependent redshift-space clustering much better than the local dark matter density model. The result infers that the age of subhalos is a key ingredient to determine the color of galaxies.