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 analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

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 investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas.

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.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

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.

We present Gemini and Keck spectroscopic redshifts and velocity dispersions for twenty clusters detected via the Sunyaev-Zel'dovich (SZ) effect by the Planck space mission, with estimated masses in the range $2.3 \times 10^{14} M_{\odot} < M < 9.4 \times 10^{14} M_{\odot}$. Cluster members were selected for spectroscopic follow-up with Palomar, Gemini and Keck optical and (in some cases) infrared imaging. Seven cluster redshifts were measured for the first time with this observing campaign, including one of the most distant Planck clusters confirmed to date, at $z=0.782\pm0.010$, PSZ2 G085.95+25.23. The spectroscopic redshift catalogs of members of each confirmed cluster are included as on-line tables. We show the galaxy redshift distributions and measure the cluster velocity dispersions. The cluster velocity dispersions obtained in this paper were used in a companion paper to measure the Planck mass bias and to constrain the cluster velocity bias.

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.

Recently, the optical counterpart of a gravitational wave source GW170817 has been identified in NGC 4993 galaxy. Together with evidence from observations in electromagnetic waves, the event has been suggested as a result of a merger of two neutron stars. We analyze the multi-wavelength data to characterize the host galaxy property and its distance to examine if the properties of NGC 4993 are consistent with this picture. Our analysis shows that NGC 4993 is a bulge-dominated galaxy with reff ~ 2-3 kpc and the Sersic index of n = 3-4 for the bulge component. The spectral energy distribution from 0.15 to 24 micron indicates that this galaxy has no significant ongoing star formation, the mean stellar mass of (0.3 - 1.2) times 10^11 Msun,the mean stellar age greater than ~3 Gyr, and the metallicity of about 20% to 100% of solar abundance. Optical images reveal dust lanes and extended features that suggest a past merging activity. Overall, NGC 4993 has characteristics of normal, but slightly disturbed elliptical galaxies. Furthermore, we derive the distance to NGC 4993 with the fundamental plane relation using 17 parameter sets of 7 different filters and the central stellar velocity dispersion from literature, finding an angular diameter distance of 37.7 +- 8.7 Mpc. NGC 4993 is similar to some host galaxies of short gamma-ray bursts but much different from those of long gamma-ray bursts, supporting the picture of GW170817 as a result of a merger of two NSs.

We measure the Planck cluster mass bias using dynamical mass measurements based on velocity dispersions of a subsample of 17 Planck-detected clusters. The velocity dispersions were calculated using redshifts determined from spectra obtained at Gemini observatory with the GMOS multi-object spectrograph. We correct our estimates for effects due to finite aperture, Eddington bias and correlated scatter between velocity dispersion and the Planck mass proxy. The result for the mass bias parameter, $(1-b)$, depends on the value of the galaxy velocity bias $b_v$ adopted from simulations: $(1-b)=(0.51\pm0.09) b_v^3$. Using a velocity bias of $b_v=1.08$ from Munari et al., we obtain $(1-b)=0.64\pm 0.11$, i.e, an error of 17% on the mass bias measurement with 17 clusters. This mass bias value is consistent with most previous weak lensing determinations. It lies within $1\sigma$ of the value needed to reconcile the Planck cluster counts with the Planck primary CMB constraints. We emphasize that uncertainty in the velocity bias severely hampers precision measurements of the mass bias using velocity dispersions. On the other hand, when we fix the Planck mass bias using the constraints from Penna-Lima et al., based on weak lensing measurements, we obtain a positive velocity bias $b_v \gtrsim 0.9$ at $3\sigma$.

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.

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.

The RGB Bayer filter system consists of a mosaic of R, G, and B filters on the grid of the photo sensors which typical commercial DSLR (Digital Single Lens Reflex) cameras and CCD cameras are equipped with. Lot of unique astronomical data obtained using an RGB Bayer filter system are available, including transient objects, e.g. supernovae, variable stars, and solar system bodies. The utilization of such data in scientific research requires that reliable photometric transformation methods are available between the systems. In this work, we develop a series of equations to convert the observed magnitudes in the RGB Bayer filter system ($R_B$, $G_B$, and $B_B$) into the Johnson-Cousins BVR filter system ($B_J$, $V_J$, and $R_C$). The new transformation equations derive the calculated magnitudes in the Johnson-Cousins filters ($B_{Jcal}$, $V_{Jcal}$, and $R_{Ccal}$) as functions of RGB magnitudes and colors. The mean differences between the transformed magnitudes and original magnitudes, i.e. the residuals, are $\Delta(B_J-B_{Jcal})$ = 0.064 mag, $\Delta(V_J-V_{Jcal})$ = 0.041 mag, and $\Delta(R_C-R_{Ccal})$ = 0.039 mag. The calculated Johnson-Cousins magnitudes from the transformation equations show a good linear correlation with the observed Johnson-Cousins magnitudes.

We study the mid-infrared (MIR) properties of the galaxies in the A2199 supercluster at z = 0.03 to understand the star formation activity of galaxy groups and clusters in the supercluster environment. Using the Wide-field Infrared Survey Explorer data, we find no dependence of mass-normalized integrated SFRs of galaxy groups/clusters on their virial masses. We classify the supercluster galaxies into three classes in the MIR color-luminosity diagram: MIR blue cloud (massive, quiescent and mostly early-type), MIR star-forming sequence (mostly late-type), and MIR green valley galaxies. These MIR green valley galaxies are distinguishable from the optical green valley galaxies, in the sense that they belong to the optical red sequence. We find that the fraction of each MIR class does not depend on virial mass of each group/cluster. We compare the cumulative distributions of surface galaxy number density and cluster/group-centric distance for the three MIR classes. MIR green valley galaxies show the distribution between MIR blue cloud and MIR SF sequence galaxies. However, if we fix galaxy morphology, early- and late-type MIR green valley galaxies show different distributions. Our results suggest a possible evolutionary scenario of these galaxies: 1) Late-type MIR SF sequence galaxies -> 2) Late-type MIR green valley galaxies -> 3) Early-type MIR green valley galaxies -> 4) Early-type MIR blue cloud galaxies. In this sequence, star formation of galaxies is quenched before the galaxies enter the MIR green valley, and then morphological transformation occurs in the MIR green valley.

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.

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.

Strong Halpha Emitters (HAEs) dominate the z~4 Lyman-break galaxy population. We have identified local analogs of these HAEs using the Sloan Digital Sky Survey (SDSS). At z<0.4, only 0.04% of galaxies are classified as HAEs with Halpha equivalent widths (>500A) comparable to that of z~4 HAEs. Local HAEs have lower stellar mass and lower ultraviolet (UV) luminosity than z~4 HAEs, yet the Halpha-to-UV luminosity ratio as well as their specific star-formation rate is consistent with that of z~4 HAEs indicating that they are scaled down versions of high-z star-forming galaxies. Compared to the previously studied local analogs of z~2 Lyman break galaxies selected using rest-frame UV, local HAEs show similar UV luminosity surface density, weaker Dn(4000) breaks, lower metallicity and lower stellar mass. This supports the idea that local HAEs are less evolved galaxies than the traditional Lyman break analogs. We are not able to constrain if the star-formation history in local HAEs is powered by mergers or by cosmological cold flow accretion. However, in the stacked spectrum, local HAEs show a strong HeII4686 emission line suggesting a population of young (<10Myr), hot, massive stars similar to that seen in some Wolf-Rayet galaxies. Low [NII]/[OIII] line flux ratios imply that local HAEs are inconsistent with being systems that host bright AGN. Instead, it is highly likely that local HAEs are galaxies with an elevated ionization parameter, either due to a high electron density or large escape fraction of hydrogen ionizing photons as in the case for Wolf-Rayet galaxies.

We have studied the star-formation and AGN activity of massive galaxies in the redshift range $z=0.4-2$, which are detected in a deep survey field using the AKARI InfraRed (IR) astronomical satellite and {\em Subaru} telescope toward the North Ecliptic Pole (NEP). The AKARI/IRC Mid-InfraRed (MIR) multiband photometry is used to trace their star-forming activities with the Polycyclic-Aromatic Hydrocarbon (PAH) emissions, which is also used to distinguish star-forming populations from AGN dominated ones and to estimate the Star Formation Rate (SFR) derived from their total emitting IR (TIR) luminosities. In combination with analyses of their stellar components, we have studied the MIR SED features of star-forming and AGN-harboring galaxies.

We present the mid-infrared (MIR) properties of galaxies within a supercluster in the North Ecliptic Pole region at z?0.087 observed with the AKARI satellite. We use data from the AKARI NEP-Wide (5.4 deg2) IR survey and the CLusters of galaxies EVoLution studies (CLEVL) mission program. We show that near-IR (3 {\mu}m)-mid- IR (11 {\mu}m) color can be used as an indicator of the specific star formation rate and the presence of intermediate age stellar populations. From the MIR observations, we find that red-sequence galaxies consist not only of passively evolving red early-type galaxies, but also of 1) "weak-SFG" (disk-dominated star-forming galaxies which have star formation rates lower by \sim 4 \times than blue-cloud galaxies), and 2) "intermediate- MXG" (bulge-dominated galaxies showing stronger MIR dust emission than normal red early-type galaxies). Those two populations can be a set of transition galaxies from blue, star-forming, late-type galaxies evolving into red, quiescent, early-type ones. We find that the weak-SFG are predominant at intermediate masses (1010M\odot < M\star < 1010.5M\odot) and are typically found in local densities similar to the outskirts of galaxy clusters. As much as 40% of the supercluster member galaxies in this mass range can be classified as weak-SFGs, but their proportion decreases to < 10% at larger masses (M\star > 1010.5 M\odot) at any galaxy density. The fraction of the intermediate-MXG among red- sequence galaxies at 1010M\odot < M\star < 1011M\odot also decreases as the density and mass increase. In particular, \sim42% of the red-sequence galaxies with early-type morphologies are classified as intermediate-MXG at intermediate densities. These results suggest that the star formation activity is strongly dependent on the stellar mass, but that the morphological transformation is mainly controlled by the environment.

The WFC3 Infrared Spectroscopic Parallel Survey (WISP) uses the Hubble Space Telescope (HST) infrared grism capabilities to obtain slitless spectra of thousands of galaxies over a wide redshift range including the peak of star formation history of the Universe. We select a population of very strong emission-line galaxies with rest-frame equivalent widths higher than 200 A. A total of 176 objects are found over the redshift range 0.35 < z < 2.3 in the 180 arcmin^2 area we analyzed so far. After estimating the AGN fraction in the sample, we show that this population consists of young and low-mass starbursts with higher specific star formation rates than normal star-forming galaxies at any redshift. After spectroscopic follow-up of one of these galaxies with Keck/LRIS, we report the detection at z = 0.7 of an extremely metal-poor galaxy with 12+Log(O/H)= 7.47 +- 0.11. The nebular emission-lines can substantially affect the broadband flux density with a median brightening of 0.3 mag, with examples producing brightening of up to 1 mag. The presence of strong emission lines in low-z galaxies can mimic the color-selection criteria used in the z ~ 8 dropout surveys. In order to effectively remove low redshift interlopers, deep optical imaging is needed, at least 1 mag deeper than the bands in which the objects are detected. Finally, we empirically demonstrate that strong nebular lines can lead to an overestimation of the mass and the age of galaxies derived from fitting of their SED. Without removing emission lines, the age and the stellar mass estimates are overestimated by a factor of 2 on average and up to a factor of 10 for the high-EW galaxies. Therefore the contribution of emission lines should be systematically taken into account in SED fitting of star-forming galaxies at all redshifts.

We present evidence for unusually strong Halpha emission in galaxies with spectroscopic redshifts in the range of 3.8<z<5.0, over the Great Observatories Origins Deep Survey (GOODS) fields. Among 74 galaxies detected in the Spitzer IRAC 3.6 and 4.5um bands, more than 70% of the galaxies show clear excess at 3.6um compared to the expected flux density from stellar continuum only. We provide evidence that this 3.6um excess is due to Halpha emission redshifted into the 3.6um band, and classify these 3.6um excess galaxies to be Halpha emitter (HAE) candidates. The selection of HAE candidates using an excess in broad-band filters is sensitive to objects whose rest-frame H$\alpha$ equivalent width is larger than 350A, star formation rates of 20-500 Msun/yr. The Halpha-to-UV luminosity ratio of HAEs is on average larger than that of local starbursts. Possible reasons for such strong Halpha emission in these galaxies include different dust extinction properties, young stellar population age, continuous star formation history, low metallicity, and top-heavy stellar initial mass function. Although the correlation between UV slope beta and L(Ha)/L(UV) raises the possibility that HAEs prefer a dust extinction curve which is steeper in the UV, the most dominant factor that results in strong Halpha emission appears to be star formation history. The Halpha equivalent widths of HAEs are large despite their relatively old stellar population ages constrained by SED fitting, suggesting that at least 60% of HAEs produce stars at a constant rate. Under the assumption that the cold gas supply is sustained, HAEs are able to produce ~50% of the observed stellar mass density at z~3 that is caught in massive (M*>10^11 Msun) galaxies. This 'strong Halpha phase' of star formation plays a dominant role in early galaxy growth, being a feasible progenitors of massive red galaxies at lower redshifts.

We present the mid-infrared (MIR) observation of a nearby galaxy cluster, Abell 2255 by the AKARI space telescope. Using the AKARI's continuous wavelength coverage between 3-24 micron and the wide field of view, we investigate the properties of cluster member galaxies to see how the infall of the galaxies, the cluster substructures, and the cluster-cluster merger influence their evolution. We show that the excess of MIR (11 micron) flux is a good indicator to discriminate galaxies at different evolutionary stages, and divide galaxies into three classes accordingly : strong MIR-excess (N3-S11>0.2) galaxies that include both unobscured and obscured star-forming galaxies, weak MIR-excess (-2.0<N3-S11<-1.2) galaxies that are quiescent, old (>5 Gyr) galaxies where the MIR emission arises mainly from the circumstellar dust around AGB stars, and intermediate MIR-excess (-1.2<N3-S11<0.2) galaxies in between the two classes that are less than a few Gyrs old past the prime star formation activity. With the MIR-excess diagnostics, we investigate how local and cluster-scale environments affect the individual galaxies. We derive the total star formation rate of ~130 Msun/yr for A2255 using the strong MIR-excess galaxies, which is consistent with other clusters of galaxies at similar redshifts and with similar masses. We find no strong evidence that supports enhanced star formation neither inside the cluster nor in the substructure region. The intermediate MIR-excess galaxies, representing galaxies in transition from star-forming galaxies to quiescent galaxies, are located preferentially at the medium density region or cluster substructures. Our findings suggest that galaxies are being transformed from star-forming galaxies into red, quiescent galaxies from the infall region through near the core, which can be well-explained by the ram-pressure stripping as previous simulation suggests.

We present the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. WISP is obtaining slitless, near-infrared grism spectroscopy of ~ 90 independent, high-latitude fields by observing in the pure parallel mode with Wide Field Camera-3 on the Hubble Space Telescope for a total of ~ 250 orbits. Spectra are obtained with the G102 (lambda=0.8-1.17 microns, R ~ 210) and G141 grisms (lambda=1.11-1.67 microns, R ~ 130), together with direct imaging in the J- and H-bands (F110W and F140W, respectively). In the present paper, we present the first results from 19 WISP fields, covering approximately 63 square arc minutes. For typical exposure times (~ 6400 sec in G102 and ~ 2700 sec in G141), we reach 5-sigma detection limits for emission lines of 5 x 10^(-17) ergs s^(-1) cm^(-2) for compact objects. Typical direct imaging 5sigma-limits are 26.8 and 25.0 magnitudes (AB) in F110W and F140W, respectively. Restricting ourselves to the lines measured with highest confidence, we present a list of 328 emission lines, in 229 objects, in a redshift range 0.3 < z < 3. The single-line emitters are likely to be a mix of Halpha and [OIII]5007,4959 A, with Halpha predominating. The overall surface density of high-confidence emission-line objects in our sample is approximately 4 per arcmin^(2).These first fields show high equivalent width sources, AGN, and post starburst galaxies. The median observed star formation rate of our Halpha selected sample is 4 Msol/year. At intermediate redshifts, we detect emission lines in galaxies as faint as H_140 ~ 25, or M_R < -19, and are sensitive to star formation rates down to less than 1 Msol/year. The slitless grisms on WFC3 provide a unique opportunity to study the spectral properties of galaxies much fainter than L* at the peak of the galaxy assembly epoch.

Using an AKARI multi-wavelength mid-infrared (IR) survey, we identify luminous starburst galaxies at z> 0.5 based on the PAH luminosity, and investigate the nature of these PAH-selected starbursts. An extragalactic survey with AKARI towards the north ecliptic pole (NEP), the NEP-Deep survey, is unique in terms of a comprehensive wavelength coverage from 2 to 24um using all 9 photometric bands of the InfraRed Camera (IRC). This survey allows us to photometrically identify galaxies whose mid-IR emission is clearly dominated by PAHs. We propose a single colour selection method to identify such galaxies, using two mid-IR flux ratios at 11-to-7um and 15-to-9um (PAH-to-continuum flux ratio in the rest-frame), which are useful to identify starburst galaxies at z~0.5 and 1, respectively. We perform a fitting of the spectral energy distributions (SEDs) from optical to mid-IR wavelengths, using an evolutionary starburst model with a proper treatment of radiative transfer (SBURT), in order to investigate their nature. The SBURT model reproduces observed optical-to-mid-IR SEDs of more than a half of PAH-selected galaxies. Based on the 8um luminosity, we find ultra luminous infrared galaxies (ULIRGs) among PAH-selected galaxies. Their PAH luminosity is higher than local ULIRGs with a similar luminosity, and the PAH-to-total IR luminosity ratio is consistent with that of less luminous starburst galaxies. They are a unique galaxy population at high redshifts and we call these PAH-selected ULIRGs "PAH-luminous" galaxies. Although they are not as massive as submillimetre galaxies at z~2, they have the stellar mass of >3x10^{10} Msun and therefore moderately massive.

Dust-obscured star-formation becomes much more important with increasing intensity, and increasing redshift. We aim to reveal cosmic star-formation history obscured by dust using deep infrared observation with the AKARI. We construct restframe 8um, 12um, and total infrared (TIR) luminosity functions (LFs) at 0.15<z<2.2 using 4128 infrared sources in the AKARI NEP-Deep field. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24um) by the AKARI satellite allows us to estimate restframe 8um and 12um luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. We have found that all 8um (0.38<z<2.2), 12um (0.15<z<1.16), and TIR LFs (0.2<z<1.6), show a continuous and strong evolution toward higher redshift. In terms of cosmic infrared luminosity density (Omega_IR), which was obtained by integrating analytic fits to the LFs, we found a good agreement with previous work at z<1.2, and that the Omega_IR evolves as propto (1+z)^4.4+-1.0. When we separate contributions to Omega_IR by LIRGs and ULIRGs, we found more IR luminous sources are increasingly more important at higher redshift. We found that the ULIRG (LIRG) contribution increases by a factor of 10 (1.8) from z=0.35 to z=1.4.

We aim to reveal environmental dependence of infrared luminosity functions (IR LFs) of galaxies at z~0.8 using the AKARI satellite. We construct restframe 8um IR LFs in the cluster region RXJ1716.4+6708 at z=0.81, and compare them with a blank field using the AKARI North Ecliptic Pole deep field data at the same redshift. AKARI's wide field of view (10'x10') is suitable to investigate wide range of galaxy environments. AKARI's 15um filter is advantageous here since it directly probes restframe 8um at z~0.8, without relying on a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. We have found that cluster IR LFs at restframe 8um have a factor of 2.4 smaller L^* and a steeper faint-end slope than that of the field. Confirming this trend, we also found that faint-end slopes of the cluster LFs becomes flatter and flatter with decreasing local galaxy density. These changes in LFs cannot be explained by a simple infall of field galaxy population into a cluster. Physics that can preferentially suppress IR luminous galaxies in high density regions is required to explain the observed results.

We present the {\it AKARI} InfraRed Camera (IRC) imaging observation of early-type galaxies in A2218 at z $\simeq$ 0.175. Mid-infrared (MIR) emission from early-type galaxies traces circumstellar dust emission from AGB stars or/and residual star formation. Including the unique imaging capability at 11 and 15 $\mu$m, our {\it AKARI} data provide an effective way to investigate MIR properties of early-type galaxies in the cluster environment. Among our flux-limited sample of 22 red sequence early-type galaxies with precise dynamical and line strength measurements ($<$ 18 mag at 3 $\mu m$), we find that at least 41% have MIR-excess emission. The $N3-S11$ versus $N3$ (3 and 11 $\mu$m) color-magnitude relation shows the expected blue sequence, but the MIR-excess galaxies add a red wing to the relation especially at the fainter end. A SED analysis reveals that the dust emission from AGB stars is the most likely cause for the MIR-excess, with low level of star formation being the next possible explanation. The MIR-excess galaxies show a wide spread of $N3-S11$ colors, implying a significant spread (2--11 Gyr) in the estimated mean ages of stellar populations. We study the environmental dependence of MIR-excess early-type galaxies over an area out to a half virial radius ($\sim$1 Mpc). We find that the MIR-excess early-type galaxies are preferentially located in the outer region. From these evidences, we suggest that the fainter, MIR-excess early-type galaxies have just joined the red sequence, possibly due to the infall and subsequent morphological/spectral transformation induced by the cluster environment.

We determine the global star formation rate density at 0.7<z<1.9 using emission-line selected galaxies identified in Hubble Space Telescope Near Infrared Camera and Multi-Object Spectrograph (HST-NICMOS) grism spectroscopy observations. Observing in pure parallel mode throughout HST Cycles 12 and 13, our survey covers ~104 arcmin2 from which we select 80 galaxies with likely redshifted Ha emission lines. In several cases, a somewhat weaker [OIII] doublet emission is also detected. The Ha luminosity range of the emission-line galaxy sample is 4.4 x 10^41 < L(Ha) < 1.5 x 10^43 erg/s. In this range, the luminosity function is well described by a Schechter function with phi* = (4.24\pm3.55) x 10^-3 Mpc^-3, L* = (2.88\pm1.58) x 10^42 erg/s, and alpha = -1.39\pm0.43. We derive a volume-averaged star formation rate density of 0.138\pm0.058 Msun/yr/Mpc3 at z=1.4 without an extinction correction. Subdividing the redshift range, we find star formation rate densities of 0.088\pm0.056 Msun/yr/Mpc3 at z=1.1 and 0.265\pm0.174 Msun/yr/Mpc3 at z=1.6. The overall star formation rate density is consistent with previous studies using Ha when the same average extinction correction is applied, confirming that the cosmic peak of star formation occurs at z>1.5.

We present the survey strategy and the data characteristics of the North Ecliptic Pole (NEP) Wide Survey of AKARI. The survey was carried out for about one year starting from May 2006 with 9 passbands from 2.5 to 24 micron and the areal coverage of about 5.8 sq. degrees centered on NEP. The survey depth reaches to 21.8 AB magnitude near infrared (NIR) bands, and ~ 18.6 AB maggnitude at the mid infrared (MIR) bands such as 15 and 18 micron. The total number of sources detected in this survey is about 104,000, with more sources in NIR than in the MIR. We have cross matched infrared sources with optically identified sources in CFHT imaging survey which covered about 2 sq. degrees within NEP-Wide survey region in order to characterize the nature of infrared sources. The majority of the mid infrared sources at 15 and 18 micron band are found to be star forming disk galaxies, with smaller fraction of early type galaxies and AGNs. We found that a large fraction (60~80 %) of bright sources in 9 and 11 micron stars while stellar fraction decreases toward fainter sources. We present the histograms of the sources at mid infrared bands at 9, 11, 15 and 18 micron. The number of sources per magnitude thus varies as m^0.6 for longer wavelength sources while shorter wavelength sources show steeper variation with m, where m is the AB magnitude.

We present results from the first phase of the Seoul National University Bright Quasar Survey in Optical (SNUQSO) as well as its basic observational setup. Previous and current large-area surveys have been successful in identifying many quasars, but they could have missed bright quasars due to their survey design. In order to help complete the census of bright quasars, we have performed spectroscopic observations of new bright quasar candidates selected from various methods based on optical colors, near-infrared colors, radio, and X-ray data. In 2005/2006, we observed 55 bright quasar candidates using the Bohyunsan Optical Echelle Spectrograph (BOES) on the 1.8 m telescope at the Bohyunsan Optical Astronomy Observatory in Korea. We identify 14 quasars/Seyferts from our observation, including an optically bright quasar with i=14.98 mag at z=0.092 (SDSS J003236.59-091026.2). Non-quasar/Seyfert objects are found to be mostly stars, among which there are five M-type stars and one cataclysmic variable. Our result shows that there still exist bright quasars to be discovered. However, at the same time, we conclude that finding new bright quasars in high Galactic latitude regions is very challenging and that the existing compilation of optically bright quasars is nearly complete in the northern hemisphere.

We present the properties of 11 $\mu$m selected sources detected in the early data of the North Ecliptic Pole Deep (NEP-Deep) Survey of AKARI. The data set covers 6 wavelength bands from 2.5 to 11 $\mu$m, with the exposure time of 10 ~ 20 minutes. This field lies within the CFHT survey with four filter bands ($g', r', i',z'), enabling us to establish nearly continuous spectral energy distributions (SEDs) for wavelengths ranging from 0.4 to 11 $\mu$m. The main sample studied here consists of 71 sources whose 11 $\mu$m AB magnitudes are equal to or brighter than 18.5 (251 $\mu$Jy), which is complete to more than 90%. The 11 $\mu$m band has an advantage of sampling star forming galaxies with low to medium redshifts since the prominent PAH feature shifts into this band. As expected, we find that the majority (~68%) of 11 $\mu$m bright sources are star forming galaxies at 0.2 < z < 0.7 with $L_{IR} ~ 10^{10}$ -- $10^{12} L_{\odot}$ based on the detailed modelling of SEDs. We also find four AGNs lying at various redshifts in the main sample. In addition, we discuss a few sources which have non-typical SEDs of the main sample, including a brown dwarf candidate, a steep power-law source, flat spectrum object, and an early-type galaxy at moderate redshift.

We present deep u*-, and g-band images taken with the MegaCam on the 3.6 m Canada-France-Hawaii Telescope (CFHT) to support the extragalactic component of the Spitzer First Look Survey (hereafter, FLS). In this paper we outline the observations, present source catalogs and characterize the completeness, reliability, astrometric accuracy and number counts of this dataset. In the central 1 deg2 region of the FLS, we reach depths of g~26.5 mag, and u*~26.2 mag (AB magnitude, 5$\sigma$ detection over a 3" aperture) with ~4 hours of exposure time for each filter. For the entire FLS region (~5 deg2 coverage), we obtained u*-band images to the shallower depth of u*=25.0--25.4 mag (5$\sigma$, 3" aperture). The average seeing of the observations is 0.85" for the central field, and ~1.00" for the other fields. Astrometric calibration of the fields yields an absolute astrometric accuracy of 0.15" when matched with the SDSS point sources between 18<g<22. Source catalogs have been created using SExtractor. The catalogs are 50% complete and greater than 99.3% reliable down to g~26.5 mag and u*~26.2 mag for the central 1 deg2 field. In the shallower u*-band images, the catalogs are 50% complete and 98.2% reliable down to 24.8--25.4 mag. These images and source catalogs will serve as a useful resource for studying the galaxy evolution using the FLS data.

We investigate the properties of 1088 Lyman Break Galaxies (LBGs) at z~3 selected from a ~2.63$ deg2 sub-region of the First Look Survey field using the ground-based multi-color data and the Spitzer Space Telescope mid-infrared data at 3--8 and 24 um. With the wide area and the broad wavelength coverage, we sample a large number of ``rare'' u-band dropouts which are massive (M* > 10^11 Msun), allowing us to perform a statistical analysis of these subsets of LBGs that have not been studied in detail. Optically bright (R(AB) < 24.5 mag) LBGs detected in mid-infrared (S_{3.6um} > 6 uJy) reside at the most massive and dusty end of the LBG population, with relatively high and tight $M/L$ in rest-frame near-infrared. Most infrared-luminous LBGs (S_{24um} > 100 uJy) are dusty star-forming galaxies with star formation rates of 100--1000 Msun/yr, total infrared luminosity of > 10^12 Lsun. By constructing the UV luminosity function of massive LBGs, we estimate that the lower limit for the star formation rate density from LBGs more massive than 10^11 Msun at z~3 is > 3.3 x 10^-3 Msun/yr/Mpc^3, showing for the first time that the UV-bright population of massive galaxies alone contributes significantly to the global star formation rate density at z~3. When combined with the star formation rate densities at z < 2, our result reveals a steady increase in the contribution of massive galaxies to the global star formation from z=0 to z=3, providing strong support to the downsizing of galaxy formation.